Baozhu Wei
Yang Liu3†
Zhi Luo- 1Department of Cardiology, Zhongnan Hospital of Wuhan University, Wuhan University, Wuhan, China
- 2Institute of Myocardial Injury and Repair, Wuhan University, Wuhan, China
- 3Department of Endocrinology, China Resources and WISCO General Hospital, Wuhan, China
- 4Department of Gerontology, Zhongnan Hospital of Wuhan University, Wuhan University, Wuhan, China
Background: Several 9p21.3 variants, such as rs1333049, rs4977574, rs10757274, rs10757278, and rs10811661, identified from recent genome-wide association studies (GWASs) are reported to be associated with coronary artery disease (CAD) susceptibility but independent of dyslipidemia. This study investigated whether these 9p21.3 variants influenced lipid profiles.
Methods and results: By searching the PubMed and Cochrane databases, 101,099 individuals were included in the analysis. The consistent finding for the rs1333049 C allele on lipid profiles increased the triglyceride (TG) levels. Moreover, the rs4977574 G allele and the rs10757274 G allele, respectively, increased low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C) levels. However, the rs10811661 C allele largely reduced LDL-C levels. Subgroup analyses indicated that the effects of the rs1333049 C allele, rs4977574 G allele, and rs10757274 G allele on lipid profiles were stronger in Whites compared with Asians. In contrast, the effect of the rs10811661 C allele on lipid profiles was stronger in Asians compared with Whites.
Conclusion: The rs1333049 C allele, rs4977574 G allele, and rs10757274 G allele of lncRNA, and the rs10811661 G allele of CDKN2A/2B had a significant influence on lipid levels, which may help the understanding of the underlying mechanisms between 9p21.3 variants and CAD.
Introduction
The variant (mutation) of the 9p21.3 allele is prevalent in the general population at a rate of 50% (1–4). One copy of the 9p21.3 allele increases the risk of coronary artery disease (CAD) by 25%, whereas two copies augment the risk by 50% (1–4).
The 9p21.3 allele contains a long non-coding RNA (also known as lncRNA, ANRIL, or CDKN2BAS) and two protein-encoding genes (CDKN2A and CDKN2B). rs1333049, a genetic variant in lncRNA, is formed by a nucleotide substitution from guanine (G) to cytosine (C). In addition, rs4977574, rs10757274, and rs10757278, genetic variants in lncRNA, are formed by a nucleotide substitution from adenine (A) to guanine (G). Moreover, rs10811661, a genetic variant in CDKN2A/2B, is formed by a nucleotide substitution from thymine (T) to cytosine (C).
A series of basic studies tested the effects of gene expression levels of lncRNA and CDKN2A/2B on lipid profiles, but the results were inconsistent. For instance, the knockdown of lncRNA inhibits lipid uptake and accumulation and promotes macrophage reverse cholesterol transport (mRCT) (5), while the overexpression of lncRNA significantly increases intracellular lipid accumulation in macrophage-derived foam cells (5). In contrast, the knockdown of CDKN2B induced C/EBPα expression and lipid accumulation in 3T3-L1 cells (6). However, knockout of the 9p21 risk locus (chr4Δ70kb/Δ70kb, ApoE–/– on 129 background) (7) or CDKN2A (8) did not affect lipid levels in deficient mice. Similarly, a large body of clinical trials have detected the effects of these five variants on lipid profiles, but the results are controversial. For instance, the C allele of rs1333049 increased triglyceride (TG) (9), total cholesterol (TC) (10), and low-density lipoprotein cholesterol (LDL-C) (9, 10) and lowered high-density lipoprotein cholesterol (HDL-C) levels (9). In contrast, the G allele of rs4977574 increased TC (11) and lowered the HDL-C levels (11), the G allele of rs10757274 increased the HDL-C levels (12), the G allele of rs10757278 increased the TC (13) and LDL-C levels (13), and the C allele of rs10811661 decreased the TG (14–16), TC (14, 15), and LDL-C (15) levels and increased the HDL-C levels (15, 16). Notably, an observational study conducted by Ahmed et al. (17) indicated that the C allele of rs1333049 remodeled lipid metabolism, therefore influencing myocardial infarction risk. Moreover, a set of elegant lipidomics studies (NPHSII, Northwick Park Heart Study II) conducted by Meckelmann et al. (18) indicated that the G allele of rs10757274 remodeled CAD risk by modulating lipid metabolism. However, the results obtained from other clinical trials did not favor these findings (19–23). Hence, a meta-analysis is needed to clarify whether these variants affected lipid profiles to resolve these discrepancies.
It has been well documented that rs1333049 C allele (1), rs4977574 G allele (24), rs10757274 G allele (2), rs10757278 G allele (3), and rs10811661 T allele (25) increased the risk of CAD by 30–40% in Whites (1–3, 25) and East Asians (24, 25). Since dyslipidemia is one of the most important risk factors for CAD and accounts for at least 50% of the population-attributable risk (26), it is tempting to speculate that the increased CAD risk caused by these variants may stem from a remodeled lipid profile. Surprisingly, the increase in CAD risk attributed to these risk alleles was reported to be independent of dyslipidemia in four GWAS studies (1–4). Since dyslipidemia is closely related to the pathogenesis of CAD (26), these five variants indeed have the potential to affect lipid levels (9–18). It is reasonable to speculate that these five variants may affect lipid profiles after enhancing statistical power. Therefore, we conducted this meta-analysis in a large sample size (101,099 individuals) to identify whether these five variants remodeled lipid metabolism and to increase our understanding of the underlying mechanisms between 9p21.3 variants and CAD.
Due to diverse constraints, most independent or single clinical trials are imperfect or flawed. For instance, the sample size is too small, the gender is not balanced, the coverage of the age group is too narrow, or it is impossible to obtain richer and more accurate experimental results due to the limitations of the examining technology or experimental condition (27–29). Since these constraints are inevitable and usually result in severe deviations or heterogeneity, a common truth may be masked by these independent clinical trials (27–29). Therefore, we attempted to utilize the light of evidence-based medicine by a meta-analysis to examine the differences and to identify the sources of heterogeneity across these independent studies to reveal a common truth (whether variants of rs1333049, rs4977574, rs10757274, rs10757278, and rs10811661 in 9p21.3 are statistically impacted or with a trend to influence lipid levels). Moreover, the specific reason to select these five variants rather than other variants loci in 9p21.3, such as LINC-PINT (30), LINC00599 (30), or rs1537373 (31), is due to only these five variants with sufficient data to execute meta-analysis.
In the past few decades, circulating TG levels have been widely reported to be associated with the occurrence, progress, and prognosis of CAD. For instance, a meta-analysis of the Asia-Pacific region indicated that TG was an independent predictor of CAD (32). This is consistent with the findings of a Chinese multi-provincial cohort study (33). Moreover, another meta-analysis of 29 prospective studies (34) indicated that TG was moderately or highly associated with CAD risk. Intriguingly, it was verified by the Copenhagen City Heart Study (35), whereby increased TG levels were associated with an increased risk of myocardial infarction, ischemic heart disease, and death. Notably, the Bezafibrate Infarction Prevention trial (36) further revealed that increased TG levels were independently associated with increased mortality in patients with CAD, indicating that circulating TG was a critical risk factor for CAD and hypertriglyceridemia should not be ignored in CAD intervention.
Here, we systematically analyzed the effects of the rs1333049 C allele, rs4977574 G allele, rs10757274 G allele, rs10757278 G allele, and rs10811661 C allele in 9p21.3 on lipid profiles in 101,099 individuals by a meta-analysis.
Materials and methods
The present meta-analysis follows the Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) (37).
Literature search
A comprehensive search of the literature was executed from January 5, 2021 to February 15, 2022, by using the PubMed and Cochrane databases. The following keywords were used in the search: (“long non-coding RNA,” “lncRNA,” “ANRIL,” “CDKN2BAS,” “CDKN2A/2B,” “9p21.3,” “rs1333049,” “rs4977574,” “rs10757274,” “rs10757278,” or “rs10811661”), (“mutation,” “variation,” “mutant,” “variant,” or “polymorphism”), and (“lipids,” “lipid metabolism,” “lipoprotein,” “cholesterol,” “circulating lipids,” “blood lipids,” “plasma lipids,” “serum lipids,” or “lipid profile”).
Inclusion and exclusion criteria
The inclusion criteria were as follows: (1) Articles that detected the effect of the rs1333049 C allele, rs4977574 G allele, rs10757274 G allele, rs10757278 G allele, or rs10811661 C allele on lipid profiles; (2) articles that provided at least one of the four parameters in lipid profiles [TG, total cholesterol (TC), LDL-C, and HDL-C]; (3) articles that provided genotype frequencies of variants of rs1333049, rs4977574, rs10757274, rs10757278, and rs1081166; (4) articles that offered mean lipid levels with standard deviation (SD) or standard errors (SE) by the genotypes; (5) the interventional studies that provided pre-intervention data; and (6) the language of eligible studies restricted to English and Chinese. The exclusion criteria were as follows: (1) articles not related to rs1333049, rs4977574, rs10757274, rs10757278, and rs10811661; (2) articles in which human subjects used lipid-lowering drugs; (3) articles that did not present genotype counts; (4) studies that provided invalid data; (5) pedigree articles; (6) overlapping articles; and (7) abstract, review, case report, meta-analysis, and animal articles.
Subgroup analysis
Subgroup analysis was conducted on ethnicity/race and disease status. Ethnicity/race includes White, Asian, and other ethnicities. Disease status includes CAD, type 2 diabetes mellitus (T2DM), and healthy subjects. In some studies, the subjects were divided into more than one subpopulation (e.g., the subjects originated from a different gender or a different race). Each subpopulation was regarded as an independent comparison in this study.
Other items
Since data extraction and analysis, heterogeneity processing, sensitivity analysis, the risk bias test, and the publication bias test were adopted from previous methods, to avoid redundant descriptions, previous publication by Liu et al. (38) provides more details. Moreover, a study conducted by Phani et al. (23) only offered raw lipid data by the genotypes of rs10811661 in their Supplementary material. Therefore, we downloaded and analyzed those raw data by performing a one-way ANOVA (the Kruskal–Wallis test) in SPSS software (version 23.0, Inc., Chicago, IL, United States).
Results
Study selection
By searching the PubMed and Cochrane databases, 4,307 articles were identified. After the screening, 4,131 articles were excluded by their title and abstract. Next, 55 articles were further estimated by their contents, of which 5 articles provided lipid data by the genotypes of rs1333049 (39, 40), rs4977574 (41), rs10757274 (42), and rs10811661 (43) but expressed as a median and interquartile range (IQR), 3 articles provided lipid data by the genotypes of rs1333049 but human subjects used lipid-lowering drugs (44–46), 2 articles provided the percentage change of lipid data by the genotypes of rs1333049 (47, 48), 1 article did not present genotype counts of rs1333049 (17), and 1 article (49) provided lipid levels by the genotypes of rs1333049 but in an aberrant genetic model [(CG + GG) vs. CC]. Therefore, 12 articles were further excluded. Finally, 43 articles involving a total of 101,099 individuals were included in the present study (Figure 1). Of the 43 articles, 15 articles (22,513 individuals), 9 articles (31,762 individuals), 7 articles (36,636 individuals), 6 articles (5,156 individuals), and 8 articles (6,028 individuals) were identified for the effects of rs1333049 C allele, rs4977574 G allele, rs10757274 G allele, rs10757278 G allele, and rs10811661 C allele on lipid profiles, respectively.
Figure 1. A flow diagram of the study’s selection process.
Characteristics of the included studies are presented in Supplementary Table 1. Circulating lipid levels by the genotypes rs1333049, rs4977574, rs10757274, rs10757278, and rs10811661 are presented in Supplementary Tables 2–6, respectively.
Effect of rs1333049 C allele on lipid profile
All the results stated below were statistically analyzed from the included studies that eliminated heterogeneity (please see “recalculated results that eliminated heterogeneity” in Supplementary Tables 7–11 for more details). The consistent finding for the effects of the rs1333049 C allele on lipid metabolism (Supplementary Table 7 and Supplementary Figure 1) was a slight increase in TG levels (Figure 2). The subgroup analysis indicated that the effect of the rs1333049 C allele on TG levels was observed in Whites (Supplementary Table 7), indicating that Whites with the rs1333049 C allele had an increased risk of CAD. Intriguingly, this speculation was supported by the present analysis results, whereby the C allele of rs1333049 significantly increased the TG levels in patients with CAD (Supplementary Table 7). Meanwhile, the C allele of rs1333049 showed a statistical influence on the TC levels in patients with CAD (Supplementary Table 7).
Figure 2. A forest plot of lncRNA rs1333049 variant with circulating TG levels.
Effect of rs4977574 G allele on lipid profile
The G allele of rs4977574 slightly increased the LDL-C (Figure 3) levels and lowered the TG levels (Supplementary Figure 2). Subgroup analyses indicated that the effects of the rs4977574 G allele on the LDL-C, HDL-C, and TG levels were noted in Whites and healthy subjects (Supplementary Table 8). This indicates that Caucasians and healthy subjects with the G allele of rs4977574 were at a high risk of CAD due to the increased LDL-C levels.
Figure 3. A forest plot of lncRNA rs4977574 variant with circulating LDL-C levels.
Effect of rs10757274 G allele on lipid profile
The G allele of rs10757274 slightly increased HDL-C (Figure 4) and lowered TG levels (Supplementary Figure 3). Subgroup analyses indicated that the effects of the rs10757274 G allele on the HDL-C and TG levels were noted in Whites and healthy subjects (Supplementary Table 9). This indicates that Whites and healthy subjects with the G allele of rs10757274 may have a reduced CAD risk.
Figure 4. A forest plot of lncRNA rs10757274 variant with circulating HDL-C levels.
Effect of rs10757278 G allele on lipid profile
The present study indicated that the G allele of rs10757278 did not statistically influence lipid levels (Supplementary Figure 4 and Supplementary Table 10).
Effect of rs10811661 C allele on lipid profile
The C allele of rs10811661 significantly reduced the LDL-C (Figure 5) and TC levels (Supplementary Figure 5). Subgroup analyses indicated that the significant effects of the rs10811661 C allele on the LDL-C and TC levels were observed in Asians and in patients with T2DM (Supplementary Table 11). This indicates that Asians and patients with T2DM with the C allele of rs10811661 were at low risk of CAD. Moreover, the C allele of rs10811661 significantly reduced the TG levels in Whites (Supplementary Table 11), indicating that Whites with the rs10811661 C allele may have a reduced risk of CAD.
Figure 5. A forest plot of lncRNA rs10811661 variant with circulating LDL-C levels.
Evaluation of heterogeneity
Significant heterogeneity was detected in the analysis of the effects of the rs1333049 C allele, the rs4977574 G allele, and the rs10811661 G allele on lipid profiles (Supplementary Tables 7–10). The SMD and 95% CIs of TC and HDL-C modulated by the rs1333049 C allele did not change substantially after eliminating heterogeneity (see Supplementary Table 7 for more details). However, the standardized mean difference (SMD) and 95% CI of TG and LDL-C affected by the rs4977574 G allele and the rs10811661 C allele changed significantly after eliminating heterogeneity (as shown in Supplementary Tables 8, 11 for more details).
Sensitivity analysis
Sensitivity analysis showed that one comparison [Shakhtshneider et al. (19)] may affect the effect of the rs1333049 C allele on the LDL-C levels (Supplementary Figure 6), one comparison [Hindy et al. (50)] may affect the effects of the rs4977574 G allele on the LDL-C and HDL-C levels (Supplementary Figure 7), and one comparison [Mehramiz et al. (16)] may affect the effect of the rs10811661 C allele on the LDL-C levels (Supplementary Figure 8). However, the effects of the rs1333049 C allele, rs4977574 G allele, and rs10811661 C allele on lipid profiles did not change substantially after omitting these comparisons. This indicates that the synthetic results were robust.
Risk bias test
While analyzing the effects of the rs1333049 C allele, rs4977574 G allele, rs10757274 G allele, and rs10811661 C allele on lipid profiles, some concerns were observed in the randomization process (10–25%). However, the overall results showed a low risk of bias (75–90%) among the included studies (Figure 6). Consequently, the studies included in the meta-analysis were of relatively high quality.
Figure 6. The risk bias plot of lncRNA variant with circulating lipid levels [(A) rs1333049 with TG; (B) rs4977574 with LDL-C; (C) rs10757274 with HDL-C; and (D) rs10811661 with LDL-C].
Publication bias test
In the present study, Begg’s test did not find any publication bias. This was confirmed by Egger’s regression test (as seen in Supplementary Figures 9–13 for more details).
Discussion
In sharp contrast to previous GWAS studies (1–4), our study indicated that the rs1333049 C allele, rs4977574 G allele, rs10757274 G allele, and rs10811661 C allele were robustly associated with lipid profiles in 101,099 individuals. This may help the understanding of the underlying mechanisms between these 9p21.3 variants and CAD.
The precise mechanisms underlying the effect of these five 9p21.3 variants on lipid profiles have not been clarified yet. However, several putative regulatory pathways could be proposed to interpret its mechanisms. (1) By modulating the expression of CDKN2A/2B. CDKN2A/2B plays a crucial role in regulating reverse cholesterol transport (RCT) (5) and in maintaining lipid metabolism homeostasis (6). However, the expression of CDKN2A/2B was largely determined by the C allele of rs1333049 (51) and rs10811661 (52) and the G allele of rs4977574 (53), rs10757274 (54) and rs10757278 (55, 56). Therefore, these five variants may indirectly influence the lipid levels by impacting CDKN2A/2B expression. (2) By modulating the expression of lipid metabolism-related genes. lncRNA is a known regulator of multiple genes (57), which may remodel lipid metabolism by regulating lipid metabolism-related genes (58). Therefore, some genetic variants in 9p21.3, such as rs1333049 C allele (59, 60), rs4977574 G allele (60), rs10757274 G allele (60), rs10757278 G allele (60, 61), and rs10811661 C allele (52), may affect the lipid levels by influencing ANRIL expression. Intriguingly, this speculation was supported by a lipidomics study (18), in which the G allele of rs10757274 remodeled the lipid metabolism by altering the expression of lysophospholipids (LysoPLs), lysophosphatidic acid (LysoPA), and autotaxin (ATX). Notably, other lipid metabolism-related genes, such as, low density lipoprotein receptor (LDLR) (62), very low density lipoprotein receptor (VLDLR) (63), ATP-binding cassette transporter A1 (ABCA1) (64), apolipoprotein C2 (APOC2) (65), apolipoprotein A-I (APOA1) (66), and HMG-CoA reductase (HMGCR) (67) may also be regulated by these variants.
In the present study, the C allele of rs1333049 showed a statistical influence on TG and TC levels in patients with CAD (Supplementary Table 7), indicating the significant association between the rs1333049 C allele and increased CAD risk (1) was mediated, at least partly, by increased TG and TC levels. However, the G allele of rs4977574 simultaneously increased the LDL-C (harmful) levels and lowered the TG (beneficial) levels (Supplementary Table 8), indicating that the rs4977574 G allele had an ambiguous influence on lipid profiles. When combined with a previous GWAS study (24), whereby the G allele of rs4977574 largely increased the risk of CAD, the increased levels of LDL-C at least partly mediated the correlation between the rs4977574 G allele and increased CAD risk.
The G allele of rs10757274 was identified as a risk allele for CAD in the McPherson et al. (2) study. Since the rs10757274 G allele may interact with cigarette smoking, alcohol consumption, the presence of hypertension, the presence of diabetes, and a family history of CAD to modulate CAD risk (68, 69), and the G allele of rs10757274 slightly increased the HDL-C levels and reduced the TG levels in the present study (Supplementary Table 9). This indicates that the correlation between the rs10757274 G allele and increased CAD risk (2) was more likely mediated by other cardiovascular risk factors, such as cigarette smoking, alcohol consumption, the presence of hypertension, the presence of diabetes, and a family history of CAD (68, 69), but not ameliorated lipid metabolism (Supplementary Table 9). Intriguingly, this speculation was favored by a large-sale GWAS study conducted by Angelakopoulou et al. (70), in which the G allele of rs10757274 significantly increased the risk of CAD but independent of dyslipidemia. More large-scale clinical trials are needed to verify this speculation.
However, we did not observe the effects of the rs10757278 G allele on lipid profiles (Supplementary Table 10). Since the G allele of rs10757278 may regulate the expression of CDKN2A/2B (55, 56) and/or ANRIL (60, 61) and was closely related to the onset of myocardial infarction (3), it was reasonable to speculate that the G allele of rs10757278 very likely impacted lipid metabolism. One plausible explanation that can be proposed to interpret this phenomenon is that the sample size included for rs10757278 lipid association analysis was relatively small (as seen in Supplementary Table 10 for more details), which largely reduced the statistical power; therefore, future large-scale clinical trials are needed to verify or correct our findings.
Strikingly, the C allele of rs10811661 substantially reduced the LDL-C and TC levels in Asians (Supplementary Table 11), indicating that the correlation between the rs10811661 C allele and reduced CAD risk in Asians (25) was mediated, at least partly, by decreased LDL-C and TC levels. According to American College of Cardiology/American Heart Association (ACC/AHA) (71), European Society of Cardiology/European Atherosclerosis Society (ESC/EAS) (72), and the adult treatment panel III (ATP III) cholesterol guidelines (73), LDL-C was considered the major cause of CAD and was treated as the primary target for therapy, while other lipids were used as the secondary or supplementary therapeutic targets. Since rs10811661 C allele largely reduced the LDL-C levels, the C allele of rs10811661 may be a potential marker for dyslipidemia and/or CAD.
According to the JDSC study (74), TG was considered a risk factor for CAD comparable with LDL-C (a 1-mmol/L increase in the baseline TG and LDL-C levels were associated with 63 and 64% higher risk of CAD, respectively). In the present study, the C allele of rs1333049, the G allele of rs4977574, and the G allele of rs10757274 consistently affected TG levels (Supplementary Tables 7–9). This indicates the effects of the rs1333049 C allele, rs4977574 G allele, and rs10757274 G allele on lipid profiles predominantly in the TG levels. Since TG levels were closely related to the pathogenesis of CAD (30–34, 71), these three alleles of the three variants may be a potential marker for CAD.
Strengths and limitations
To the best of our knowledge, this is the first reliable evidence that demonstrates that the rs1333049 C allele, rs4977574 G allele, rs10757274 G allele, and rs10811661 C allele in 9p21.3 had a statistical influence on lipid profiles. Several strengths of the present study should be noted. For instance, the clinical lipid data of 101,099 individuals were included in the analysis, which increased the reliability of synthetic results due to high statistical power. Moreover, data analyses were performed after eliminating the studies with heterogeneity, which further advanced the preciseness of conclusions drawn in our study. Most importantly, our findings may help the understanding of the underlying mechanisms between variants of rs1333049, rs4977574, rs10757274, and rs10811661 in 9p21.3 and CAD. However, a large number of genes as well as some environmental factors are involved in dyslipidemia. Our study has not investigated the interaction of the 9p21.3 variant with other variant locus or environmental factors on lipid profiles due to the lack of original data from the included studies.
Conclusion
The rs1333049 C allele, rs4977574 G allele, rs10757274 G allele of lncRNA and the rs10811661 G allele of CDKN2A/2B had a significant influence on lipid levels, which may help the understanding of the underlying mechanisms between 9p21.3 variants and CAD.
Data availability statement
The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
ZL, BW, and YL conceived and designed this study as well as drafted the manuscript. HL and YP carried out the searches and collected the data. ZL performed the statistical analyses. All authors reviewed and approved the final manuscript.
Conflict of interest
The 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/fcvm.2022.946289/full#supplementary-material
References
1. Samani NJ, Erdmann J, Hall AS, Hengstenberg C, Mangino M, Mayer B, et al. Genomewide association analysis of coronary artery disease. N Engl J Med. (2007) 357:443–53. doi: 10.1056/NEJMoa072366
2. McPherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, Cox DR, et al. A common allele on chromosome 9 associated with coronary heart disease. Science. (2007) 316:1488–91. doi: 10.1126/science.1142447
3. Helgadottir A, Thorleifsson G, Manolescu A, Gretarsdottir S, Blondal T, Jonasdottir A, et al. A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science. (2007) 316:1491–3. doi: 10.1126/science.1142842
4. The Wellcome Trust Case Control Consortium. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. (2007) 447:661–78. doi: 10.1038/nature05911
5. Ou M, Li X, Zhao S, Cui S, Tu J. Long non-coding RNA CDKN2B-AS1 contributes to atherosclerotic plaque formation by forming RNA-DNA triplex in the CDKN2B promoter. EbioMedicine. (2020) 5:102694. doi: 10.1016/j.ebiom.2020.102694
6. Horswell SD, Fryer LG, Hutchison CE, Zindrou D, Speedy HE, Town MM, et al. CDKN2B expression in adipose tissue of familial combined hyperlipidemia patients. J Lipid Res. (2013) 54:3491–505. doi: 10.1194/jlr.M041814
7. Kojima Y, Ye J, Nanda V, Wang Y, Flores AM, Jarr KU, et al. Knockout of the murine ortholog to the human 9p21 coronary artery disease locus leads to smooth muscle cell proliferation, vascular calcification, and advanced atherosclerosis. Circulation. (2020) 141:1274–6. doi: 10.1161/CIRCULATIONAHA.119.043413
8. Kim JB, Deluna A, Mungrue IN, Vu C, Pouldar D, Civelek M, et al. Effect of 9p21.3 coronary artery disease locus neighboring genes on atherosclerosis in mice. Circulation. (2012) 126:1896–906. doi: 10.1161/CIRCULATIONAHA.111.064881
9. Abid K, Mili D, Kenani A. Polymorphism on chromosome 9p21.3 is associated with severity and early-onset CAD in type 2 diabetic Tunisian population. Dis Markers. (2015) 2015:792679. doi: 10.1155/2015/792679
10. Dandona S, Stewart AF, Chen L, Williams K, So D, O’Brien E, et al. Gene dosage of the common variant 9p21 predicts severity of coronary artery disease. J Am Coll Cardiol. (2010) 56:479–86. doi: 10.1016/j.jacc.2009.10.092
11. Temel şG, Ergören MÇ. The association between the chromosome 9p21 CDKN2B-AS1 gene variants and the lipid metabolism: a pre-diagnostic biomarker for coronary artery disease. Anatol J Cardiol. (2019) 21:31–8. doi: 10.14744/AnatolJCardiol.2018.90907
12. Zhao Q, Liao S, Wei H, Liu D, Li J, Zhang X, et al. CDKN2BAS polymorphisms are associated with coronary heart disease risk a Han Chinese population. Oncotarget. (2016) 7:82046–54. doi: 10.18632/oncotarget.12575
13. Bi J, Yang L, Liu D, Wu J, Tong X, Cen S, et al. Sequence variants on chromosome 9p21 are associated with ischemic stroke and the lipids level in Chinese Han population. J Stroke Cerebrovasc Dis. (2015) 24:894–900.
14. Chen Y, Zhao Y, Liu GF. Associations of polymorphisms of SLC30A8, CDKN2A/2B, HHEX and TCF7L2 gene with type 2 diabetes mellitus in Chinese Han population in northeast China. Chin J Gerontol. (2015) 35:4201–4.
15. Liu J, Wang SZ, Wang QL, Du JG, Wang BB. The correlation between blood lipid, blood glucose levels of gestational diabetes mellitus and CDKN2A/2B gene polymorphism. Labeled Immunoassays Clin Med. (2018) 25:1360–3.
16. Mehramiz M, Ghasemi F, Esmaily H, Tayefi M, Hassanian SM, Sadeghzade M, et al. Interaction between a variant of CDKN2A/B-gene with lifestyle factors in determining dyslipidemia and estimated cardiovascular risk: a step toward personalized nutrition. Clin Nutr. (2018) 37:254–61. doi: 10.1016/j.clnu.2016.12.018
17. Ahmed W, Ali IS, Riaz M, Younas A, Sadeque A, Niazi AK, et al. Association of ANRIL polymorphism (rs1333049:C>G) with myocardial infarction and its pharmacogenomic role in hypercholesterolemia. Gene. (2013) 515:416–20. doi: 10.1016/j.gene.2012.12.044
18. Meckelmann SW, Hawksworth JI, White D, Andrews R, Rodrigues P, O’Connor A, et al. Metabolic dysregulation of the lysophospholipid/autotaxin axis in the chromosome 9p21 gene SNP rs10757274. Circ Genom Precis Med. (2020) 13:e002806. doi: 10.1161/CIRCGEN.119.002806
19. Shakhtshneider E, Orlov P, Semaev S, Ivanoshchuk D, Malyutina S, Gafarov V, et al. Analysis of polymorphism rs1333049 (located at 9P21.3) in the white population of western Siberia and associations with clinical and biochemical markers. Biomolecules. (2019) 9:290. doi: 10.3390/biom9070290
20. Matsuoka R, Abe S, Tokoro F, Arai M, Noda T, Watanabe S, et al. Association of six genetic variants with myocardial infarction. Int J Mol Med. (2015) 35:1451–9. doi: 10.3892/ijmm.2015.2115
21. Nambi V, Boerwinkle E, Lawson K, Brautbar A, Chambless L, Franeschini N, et al. The 9p21 genetic variant is additive to carotid intima media thickness and plaque in improving coronary heart disease risk prediction in white participants of the atherosclerosis risk in communities (ARIC) study. Atherosclerosis. (2012) 222:135–7. doi: 10.1016/j.atherosclerosis.2012.01.028
22. Patel RS, Su S, Neeland IJ, Ahuja A, Veledar E, Zhao J, et al. The chromosome 9p21 risk locus is associated with angiographic severity and progression of coronary artery disease. Eur Heart J. (2010) 31:3017–23. doi: 10.1093/eurheartj/ehq272
23. Phani NM, Adhikari P, Nagri SK, D’Souza SC, Satyamoorthy K, Rai PS. Replication and relevance of multiple susceptibility loci discovered from genome wide association studies for type 2 diabetes in an Indian population. PLoS One. (2016) 11:e0157364. doi: 10.1371/journal.pone.0157364
24. Lee JY, Lee BS, Shin DJ, Woo Park K, Shin YA, Joong Kim K, et al. A genome-wide association study of a coronary artery disease risk variant. J Hum Genet. (2013) 58:120–6. doi: 10.1038/jhg.2012.124
25. Cheng X, Shi L, Nie S, Wang F, Li X, Xu C, et al. The same chromosome 9p21.3 locus is associated with type 2 diabetes and coronary artery disease in a Chinese Han population. Diabetes. (2011) 60:680–4. doi: 10.2337/db10-0185
26. Yusuf S, Hawken S, Ounpuu S, Dans T, Avezum A, Lanas F, et al. Effect of potentially modifiable risk factors associated with myocardial infarction in 52 countries (the INTERHEART study): case-control study. Lancet. (2004) 364:937–52. doi: 10.1016/S0140-6736(04)17018-9
27. Lipsey MW, Wilson DB. Practical Meta-Analysis. Thousand Oaks, CA: SAGE publications, Inc (2001).
28. Card NA. Applied Meta-Analysis for Social Science Research. New York, NY: Guilford Publications (2015).
29. Borenstein M, Hedges LV, Higgins JP, Rothstein HR. Introduction to Meta-Analysis. Hoboken, NJ: John Wiley & Sons (2011).
30. Sun M, Cao Y, Wang T, Liu T, An F, Wu H, et al. Association between LINC-PINT and LINC00599 gene polymorphism and the risk of steroid-induced osteonecrosis of the femoral head in the population of northern China. Steroids. (2021) 173:108886. doi: 10.1016/j.steroids.2021.108886
31. Li Y, Zhang D, Zhang Y, Xu X, Bi L, Zhang M, et al. Association of lncRNA polymorphisms with triglyceride and total cholesterol levels among myocardial infarction patients in Chinese population. Gene. (2020) 724:143684. doi: 10.1016/j.gene.2019.02.085
32. Patel A, Barzi F, Jamrozik K, Lam TH, Ueshima H, Whitlock G, et al. Serum triglycerides as a risk factor for cardiovascular diseases in the Asia-Pacific region. Circulation. (2004) 110:2678–86. doi: 10.1161/01.CIR.0000145615.33955.83
33. Liu J, Wang W, Wang M, Sun J, Liu J, Li Y, et al. Impact of diabetes, high triglycerides and low HDL cholesterol on risk for ischemic cardiovascular disease varies by LDL cholesterol level: a 15-year follow-up of the Chinese multi-provincial cohort study. Diabetes Res Clin Pract. (2012) 96:217–24. doi: 10.1016/j.diabres.2011.12.018
34. Sarwar N, Danesh J, Eiriksdottir G, Sigurdsson G, Wareham N, Bingham S, et al. Triglycerides and the risk of coronary heart disease: 10,158 incident cases among 262,525 participants in 29 Western prospective studies. Circulation. (2007) 115:450–8. doi: 10.1161/CIRCULATIONAHA.106.637793
35. Freiberg JJ, Tybjaerg-Hansen A, Jensen JS, Nordestgaard BG. Nonfasting triglycerides and risk of ischemic stroke in the general population. JAMA. (2008) 300:2142–52. doi: 10.1001/jama.2008.621
36. Klempfner R, Erez A, Sagit BZ, Goldenberg I, Fisman E, Kopel E, et al. Elevated triglyceride level is independently associated with increased all-cause mortality in patients with established coronary heart disease: twenty-two-year follow-up of the bezafibrate infarction prevention study and registry. Circ Cardiovasc Qual Outcomes. (2016) 9:100–8. doi: 10.1161/CIRCOUTCOMES.115.002104
37. Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gøtzsche PC, Ioannidis JP, et al. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate health care interventions: explanation and elaboration. J Clin Epidemiol. (2009) 62:1–34. doi: 10.1016/j.jclinepi.2009.06.006
38. Liu F, Wang S, Luo Z. Associations of the miRNA-146a rs2910164 and the miRNA-499a rs3746444 polymorphisms with plasma lipid levels: a meta-analysis. Front Genet. (2021) 12:746686. doi: 10.3389/fgene.2021.746686
39. Reichert S, Seitter L, Schaller HG, Schlitt A, Schulz S. ANRIL polymorphisms (rs1333049 and rs3217992) in relation to plasma CRP levels among in-patients with CHD. Cytokine. (2020) 127:154932. doi: 10.1016/j.cyto.2019.154932
40. Aschauer S, Mittermayer F, Wagner CC, Schmidt WM, Brunner M, Haslacher H, et al. Forearm vasodilator reactivity in homozygous carriers of the 9p21.3 rs1333049 G>C polymorphism. Eur J Clin Invest. (2010) 40:700–5. doi: 10.1111/j.1365-2362.2010.02321.x
41. Hua L, Yuan JX, He S, Zhao CH, Jia QW, Zhang J, et al. Analysis on the polymorphisms of site RS4977574, and RS1333045 in region 9p21 and the susceptibility of coronary heart disease in Chinese population. BMC Med Genet. (2020) 21:36. doi: 10.1186/s12881-020-0965-x
42. Paynter NP, Chasman DI, Buring JE, Shiffman D, Cook NR, Ridker PM. Cardiovascular disease risk prediction with and without knowledge of genetic variation at chromosome 9p21.3. Ann Intern Med. (2009) 150:65–72. doi: 10.7326/0003-4819-150-2-200901200-00003
43. Totomoch-Serra A, Muñoz ML, Burgueño J, Revilla-Monsalve MC, Diaz-Badillo A. Association of common polymorphisms in the VEGFA and SIRT1 genes with type 2 diabetes-related traits in Mexicans. Arch Med Sci. (2018) 14:1361–73. doi: 10.5114/aoms.2018.74757
44. Farzaneh-Far R, Na B, Schiller NB, Whooley MA. Lack of association of chromosome 9p21.3 genotype with cardiovascular structure and function in persons with stable coronary artery disease: the heart and soul study. Atherosclerosis. (2009) 205:492–6. doi: 10.1016/j.atherosclerosis.2008.12.026
45. Muendlein A, Saely CH, Rhomberg S, Sonderegger G, Loacker S, Rein P, et al. Evaluation of the association of genetic variants on the chromosomal loci 9p21.3, 6q25.1, and 2q36.3 with angiographically characterized coronary artery disease. Atherosclerosis. (2009) 205:174–80. doi: 10.1016/j.atherosclerosis.2008.10.035
46. Vakalis K, Bechlioulis A, Naka KK, Chatzikyriakidou A, Gartzonika K, Vezyraki P, et al. Role of 9p21 and 2q36 variants and arterial stiffness in the prediction of coronary artery disease. Eur J Clin Invest. (2014) 44:784–94. doi: 10.1111/eci.12295
47. Wang W, Peng W, Zhang X, Lu L, Zhang R, Zhang Q, et al. Chromosome 9p21.3 polymorphism in a Chinese Han population is associated with angiographic coronary plaque progression in non-diabetic but not in type 2 diabetic patients. Cardiovasc Diabetol. (2010) 9:33. doi: 10.1186/1475-2840-9-33
48. Dutta A, Henley W, Lang IA, Murray A, Guralnik J, Wallace RB, et al. The coronary artery disease-associated 9p21 variant and later life 20-year survival to cohort extinction. Circ Cardiovasc Genet. (2011) 4:542–8. doi: 10.1161/CIRCGENETICS.111.960146
49. Ellis KL, Pilbrow AP, Frampton CM, Doughty RN, Whalley GA, Ellis CJ, et al. A common variant at chromosome 9P21.3 is associated with age of onset of coronary disease but not subsequent mortality. Circ Cardiovasc Genet. (2010) 3:286–93. doi: 10.1161/CIRCGENETICS.109.917443
50. Hindy G, Ericson U, Hamrefors V, Drake I, Wirfält E, Melander O, et al. The chromosome 9p21 variant interacts with vegetable and wine intake to influence the risk of cardiovascular disease: a population based cohort study. BMC Med Genet. (2014) 15:1220. doi: 10.1186/s12881-014-0138-x
51. Congrains A, Kamide K, Oguro R, Yasuda O, Miyata K, Yamamoto E, et al. Genetic variants at the 9p21 locus contribute to atherosclerosis through modulation of ANRIL and CDKN2A/B. Atherosclerosis. (2012) 220:449–55. doi: 10.1016/j.atherosclerosis.2011.11.017
52. Kong Y, Sharma RB, Ly S, Stamateris RE, Jesdale WM, Alonso LC. CDKN2A/B T2D genome-wide association study risk SNPs impact locus gene expression and proliferation in human islets. Diabetes. (2018) 67:872–84. doi: 10.2337/db17-1055
53. Taheri M, Pouresmaeili F, Omrani MD, Habibi M, Sarrafzadeh S, Noroozi R, et al. Association of ANRIL gene polymorphisms with prostate cancer and benign prostatic hyperplasia in an Iranian population. Biomark Med. (2017) 11:413–22. doi: 10.2217/bmm-2016-0378
54. Zhuang J, Peng W, Li H, Wang W, Wei Y, Li W, et al. Methylation of p15INK4b and expression of ANRIL on chromosome 9p21 are associated with coronary artery disease. PLoS One. (2012) 7:e47193. doi: 10.1371/journal.pone.0047193
55. Harismendy O, Notani D, Song X, Rahim NG, Tanasa B, Heintzman N, et al. 9p21 DNA variants associated with coronary artery disease impair interferon-γ signalling response. Nature. (2011) 470:264–8. doi: 10.1038/nature09753
56. Svensson PA, Wahlstrand B, Olsson M, Froguel P, Falchi M, Bergman RN, et al. CDKN2B expression and subcutaneous adipose tissue expandability: possible influence of the 9p21 atherosclerosis locus. Biochem Biophys Res Commun. (2014) 446:1126–31. doi: 10.1016/j.bbrc.2014.03.075
57. Milligan MJ, Lipovich L. Pseudogene-derived lncRNAs: emerging regulators of gene expression. Front Genet. (2015) 5:476. doi: 10.3389/fgene.2014.00476
58. Zeng Y, Ren K, Zhu X, Zheng Z, Yi G. Long noncoding RNAs: advances in lipid metabolism. Adv Clin Chem. (2018) 87:1–36.
59. Yang J, Gu L, Guo X, Huang J, Chen Z, Huang G, et al. LncRNA ANRIL expression and ANRIL gene polymorphisms contribute to the risk of ischemic stroke in the Chinese Han population. Cell Mol Neurobiol. (2018) 38:1253–69. doi: 10.1007/s10571-018-0593-6
60. Shanker J, Arvind P, Jambunathan S, Nair J, Kakkar V. Genetic analysis of the 9p21.3 CAD risk locus in Asian Indians. Thromb Haemost. (2014) 111:960–9. doi: 10.1160/TH13-08-0706
61. Zhang W, Chen Y, Liu P, Chen J, Song L, Tang Y, et al. Variants on chromosome 9p21.3 correlated with ANRIL expression contribute to stroke risk and recurrence in a large prospective stroke population. Stroke. (2012) 43:14–21. doi: 10.1161/STROKEAHA.111.625442
62. Mitchel K, Theusch E, Cubitt C, Dosé AC, Stevens K, Naidoo D, et al. RP1-13D10.2 Is a novel modulator of statin-induced changes in cholesterol. Circ Cardiovasc Genet. (2016) 9:223–30. doi: 10.1161/CIRCGENETICS.115.001274
63. Liu H, Zhou T, Wang B, Li L, Ye D, Yu S. Identification and functional analysis of a potential key lncRNA involved in fat loss of cancer cachexia. J Cell Biochem. (2018) 119:1679–88. doi: 10.1002/jcb.26328
64. Hu YW, Yang JY, Ma X, Chen ZP, Hu YR, Zhao JY, et al. A lincRNA-DYNLRB2-2/GPR119/GLP-1R/ABCA1-dependent signal transduction pathway is essential for the regulation of cholesterol homeostasis. J Lipid Res. (2014) 55:681–97. doi: 10.1194/jlr.M044669
65. Li P, Ruan X, Yang L, Kiesewetter K, Zhao Y, Luo H, et al. A liver-enriched long non-coding RNA, lncLSTR, regulates systemic lipid metabolism in mice. Cell Metab. (2015) 21:455–67. doi: 10.1016/j.cmet.2015.02.004
66. Halley P, Kadakkuzha BM, Faghihi MA, Magistri M, Zeier Z, Khorkova O, et al. Regulation of the apolipoprotein gene cluster by a long noncoding RNA. Cell Rep. (2014) 6:222–30. doi: 10.1016/j.celrep.2013.12.015
67. Sallam T, Jones MC, Gilliland T, Zhang L, Wu X, Eskin A, et al. Feedback modulation of cholesterol metabolism by the lipid-responsive non-coding RNA LeXis. Nature. (2016) 534:124–8. doi: 10.1038/nature17674
68. Nawaz SK, Noreen A, Rani A, Yousaf M, Arshad M. Association of the rs10757274 SNP with coronary artery disease in a small group of a Pakistani population. Anatol J Cardiol. (2015) 15:709–15. doi: 10.5152/akd.2014.5470
69. Hu WL, Li SJ, Liu DT, Wang Y, Niu SQ, Yang XC, et al. Genetic variants on chromosome 9p21 and ischemic stroke in Chinese. Brain Res Bull. (2009) 79:431–5. doi: 10.1016/j.brainresbull.2009.04.001
70. Angelakopoulou A, Shah T, Sofat R, Shah S, Berry DJ, Cooper J, et al. Comparative analysis of genome-wide association studies signals for lipids, diabetes, and coronary heart disease: cardiovascular biomarker genetics collaboration. Eur Heart J. (2012) 33:393–407. doi: 10.1093/eurheartj/ehr225
71. Grundy SM, Stone NJ, Bailey AL, Beam C, Birtcher KK, Blumenthal RS, et al. 2018 AHA/ACC/AACVPR/AAPA/ABC/ACPM/ADA/AGS/APhA/ASPC/NLA/PCNA guideline on the management of blood cholesterol: a report of the American college of cardiology/American heart association task force on clinical practice guidelines. Circulation. (2019) 139:1082–143. doi: 10.1161/CIR.0000000000000624
72. Mach F, Baigent C, Catapano AL, Koskinas KC, Casula M, Badimon L, et al. 2019 ESC/EAS guidelines for the management of dyslipidaemias: lipid modification to reduce cardiovascular risk. Eur Heart J. (2020) 41:111–88. doi: 10.15829/1560-4071-2020-3826
73. National Cholesterol Education Program [NCEP] expert panel on detection, evaluation, and treatment of high blood cholesterol in adults (adult treatment panel Iii). Third report of the national cholesterol education program (NCEP) expert panel on detection, evaluation, and treatment of high blood cholesterol in adults (adult treatment panel III) final report. Circulation. (2002) 106:3143–421. doi: 10.1161/circ.106.25.3143
74. Sone H, Tanaka S, Tanaka S, Iimuro S, Oida K, Yamasaki Y, et al. Serum level of triglycerides is a potent risk factor comparable to LDL cholesterol for coronary heart disease in Japanese patients with type 2 diabetes: subanalysis of the Japan Diabetes Complications Study (JDCS). J Clin Endocrinol Metab. (2011) 96:3448–56. doi: 10.1210/jc.2011-0622
Keywords: lncRNA, CDKN2A/2B, variant, dyslipidemia, coronary artery disease
Citation: Wei B, Liu Y, Li H, Peng Y and Luo Z (2022) Effect of 9p21.3 (lncRNA and CDKN2A/2B) variant on lipid profile. Front. Cardiovasc. Med. 9:946289. doi: 10.3389/fcvm.2022.946289
Received: 17 May 2022; Accepted: 16 August 2022;
Published: 07 September 2022.
Edited by:
Prashanth N. Suravajhala, Amrita Vishwa Vidyapeetham University, IndiaReviewed by:
Aldo Grefhorst, Amsterdam University Medical Center, NetherlandsRobert N. Helsley, University of Kentucky, United States
Corina Rosales, Houston Methodist Research Institute, United States
Copyright © 2022 Wei, Liu, Li, Peng and Luo. 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.
*Correspondence: Baozhu Wei, wbz9810@163.com; Zhi Luo, 2020103030002@whu.edu.cn
†These authors have contributed equally to this work and share first authorship