Next Article in Journal
Ezetimibe Induces Paraptosis through Niemann–Pick C1-like 1 Inhibition of Mammalian-Target-of-Rapamycin Signaling in Hepatocellular Carcinoma Cells
Next Article in Special Issue
Autism Spectrum Disorder with Epilepsy: A Research Protocol for a Clinical and Genetic Study
Previous Article in Journal
DNAJB11 Mutation in ADPKD Patients: Clinical Characteristics in a Monocentric Cohort
Previous Article in Special Issue
Navigating the Complex Landscape of Fibrodysplasia Ossificans Progressiva: From Current Paradigms to Therapeutic Frontiers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Variants at the 9p21.3 Locus Are Associated with Risk for Non-Compressible Artery Disease: Results from the ARTPER Study

by
Marc Via
1,2,*,
Guillem Pera
3,
Rosa Forés
3,
Anna Costa-Garrido
3,
Antonio Heras
3,4,
José Miguel Baena-Díez
5,6,
Edurne Pedrosa
7,
Inmaculada C. Clemente
1,2,
Noemí Lamonja-Vicente
3,8,
Maria Mataró
2,8,
Pere Torán-Montserrat
3 and
M. Teresa Alzamora
3,4
1
Brainlab-Grup de Recerca en Neurociència Cognitiva, Departament de Psicologia Clínica i Psicobiologia, Institut de Neurociències, Universitat de Barcelona, 08035 Barcelona, Spain
2
Institut de Recerca Sant Joan de Déu, 08950 Esplugues de Llobregat, Spain
3
Unitat de Suport a la Recerca Metropolitana Nord, Fundació Institut Universitari per a la recerca a l’Atenció Primària de Salut Jordi Gol i Gurina (IDIAPJGol), 08303 Mataró, Spain
4
Centre d’Atenció Primària Riu Nord-Riu Sud Santa Coloma de Gramenet, Direcció d’Atenció Primària Barcelonés Nord i Maresme, Institut Català de la Salut, 08921 Santa Coloma de Gramenet, Spain
5
Centre d’Atenció Primària la Marina, Direcció d’Atenció Primària Barcelona Ciutat, Institut Català de la Salut, 08038 Barcelona, Spain
6
Institut Universitari d’Investigació en Atenció Primària Jordi Gol (IDIAP Jordi Gol), 08007 Barcelona, Spain
7
IGTP-HUGTP Biobank, Germans Trias i Pujol Research Institute (IGTP), 08916 Badalona, Spain
8
Departament de Psicologia Clínica i Psicobiologia, Institut de Neurociències, Universitat de Barcelona, 08035 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Submission received: 14 November 2023 / Revised: 10 December 2023 / Accepted: 15 December 2023 / Published: 19 December 2023
(This article belongs to the Special Issue Feature Papers in Human Genomics and Genetic Diseases 2023)

Abstract

:
Peripheral artery disease (PAD) and non-compressible artery disease (NCAD) constitute predictors of subclinical atherosclerosis easily assessed through the ankle brachial index (ABI). Although both diseases show substantial genetic influences, few genetic association studies have focused on the ABI and PAD, and none have focused on NCAD. To overcome these limitations, we assessed the role of several candidate genes on the ABI, both in its continuous distribution and in the clinical manifestations associated to its extreme values: PAD and NCAD. We examined 13 candidate genomic regions in 1606 participants from the ARTPER study, a prospective population-based cohort, with the ABI assessed through ultrasonography. Association analyses were conducted independently for individuals with PAD (ABI < 0.9) or with NCAD (ABI > 1.4) vs. healthy participants. After including potential covariates and correction for multiple testing, minor alleles in the genetic markers rs10757278 and rs1333049, both in the 9p21.3 region, were significantly associated with a decreased risk of NCAD. Associations with the ABI showed limited support to these results. No significant associations were detected for PAD. The locus 9p21.3 constitutes the first genetic locus associated with NCAD, an assessment of subclinical atherosclerosis feasible for implementation in primary healthcare settings that has been systematically neglected from genetic studies.

1. Introduction

Atherosclerotic cardiovascular diseases (CVDs) continue to be the main cause of death in developed countries, although mortality associated to them has decreased in most regions of the world in the last decades [1]. In the United States, they are estimated to represent an economic burden of USD 11.3 trillion until 2050 [2]. Among other CVDs, peripheral artery disease (PAD) is a common manifestation of atherosclerosis characterized by a narrowing of peripheral arteries that increases cardiovascular and cerebrovascular morbid-mortality incidence [3]. Around the world, 202 million people have PAD. The prevalence of PAD [assessed as an ankle brachial index (ABI) < 0.9] ranges from 3.7 to 7% in the general population of countries with a low cardiovascular risk [4,5] with an overall increase during the last decades. Worldwide, the prevalence of this disease has increased more than 17% between 2010 and 2015 [6]. Conversely, a high ABI (ABI ≥ 1.4) result indicates non-compressible artery disease (NCAD), suggesting the presence of arterial calcification in the lower limbs [7], and is also associated with an increased risk of morbid-mortality [8,9,10]. The prevalence of NCAD is 6.2% in the same population to which the present study refers [11] and is twice as frequent in males as in females. In this sense, the ABI constitutes a predictor of subclinical atherosclerosis feasible to be systematically assessed in primary healthcare settings.
The main factors associated to an increased risk of PAD are environmental factors also associated to atherosclerosis, especially smoking, age, sex, type 2 diabetes, dyslipidemia, hypertension, or obesity, among others. However, several studies indicate a remarkable genetic component in the etiology of PAD. A family history of PAD and its clinical symptoms have been identified as independent risk factors for PAD [12,13], indicating a role for genetic or shared environmental factors. Moreover, twin studies have estimated that 58% and 48% of the observed variability in PAD and in ABI values, respectively, could be attributed to additive genetic effects [14,15]. Estimates of genetic influences for arterial calcification have yielded similar results, with heritability values ranging between 42 and 51% depending on the study and the type of calcification analyzed, although only abdominal aortic and coronary atherosclerotic calcification (CAC) have been examined (reviewed in [16]).
Despite this important genetic component, the specific genetic variants involved in pathological ABI values remain elusive. In contrast to thousands of association studies for other atherosclerotic pathologies, such as coronary artery disease or stroke, few reports on the ABI and PAD, and none on NCAD, have been published. For years, studies on the genetic bases of PAD used a candidate gene approach and focused on pathways relevant in atherosclerosis, such as lipid metabolism, endothelial function, or inflammation, among others. However, a meta-analysis from the candidate gene association resource (CARe) consortium for PAD and continuous ABI did not identify any robust association after examining 11 different studies [17]. With the advent of genotyping array technologies, several genome-wide association studies (GWASs) were conducted without a clear consensus on their findings. A meta-analysis of 21 GWASs published in 2012 only identified one genome-wide significant locus (9p21.3) with the ABI and none with categorical PAD [18]. There has not been much advance in the genetics of PAD until two recent efforts, which have shed some light on the genetic determinants of PAD by combining data from hundreds of thousands of individuals. Klarin and colleagues identified 19 loci associated with PAD in the Million Veteran Program (MVP) and the UK Biobank datasets [19], and a genome-wide meta-analysis just confirmed five of these associations in an independent sample of 449,548 individuals (12,086 PAD cases) [20]. One of these studies also included ABI data in the analyses, and locus 9p21.3 was the only genome-wide significant region associated with continuous ABI [19].
No previous studies have examined genetic determinants of NCAD assessed through the ABI. Moreover, genetic studies on the ABI have actively excluded individuals with values higher than 1.4 from their analyses (e.g., [17,18,19]). As for arterial calcification, some studies have focused on coronary artery calcification (CAC), another marker of subclinical atherosclerosis. In this sense, the most comprehensive meta-analysis of genome-wide association studies for CAC identified 11 independent risk loci as influencing CAC [21]. In addition, different association studies have screened the implication of specific candidate genes on arterial calcification revealing promising markers in specific pathways such as endothelial function, cell adhesion, and extracellular matrix remodeling (as a review, see [16]).
In the present study, we aimed at clarifying the relationship between variants in several candidate genes and differences in the ABI, both in its continuous distribution and in the clinical manifestations associated to its extreme values: PAD and NCAD. To do so, we explored fifteen genetic markers in candidate regions in the well-characterized subjects of the population-based Peripheral Arterial Disease (PERART/ARTPER) study. These candidate regions include regions previously associated to the ABI and genes associated to other atherosclerotic pathologies but not specifically tested for PAD or NCAD.

2. Materials and Methods

2.1. Subjects

A detailed description of the ARTPER study has been published previously [4,22]. Briefly, the ARTPER study is an ongoing prospective, multicenter, population-based cohort of 3786 subjects older than 49 years with the aim of determining the prevalence of peripheral artery disease and related cardiovascular risk factors in the general population. Individuals were recruited between 2006 and 2008, with periodic phone and in situ follow-up interviews. Later, between 2011–2012, 1810 participants provided a blood sample to be used in genetic studies.
The study was approved by the local ethical review board (IDIAP Jordi Gol Foundation of Investigation in Primary Care (P11/61)), and all participants provided written informed consent. Moreover, the recommendations of the World Medical Association Declaration of Helsinki were followed.

2.2. Data Collection

The determination of the ankle brachial index (ABI) of the participants was carried under standardized conditions by two trained healthcare professionals using a standardized Doppler Ultrasonic device (Mini-Dopplex D 900-P, Huntleigh Healthcare, 8 MHz). A detailed description of the procedure has been previously published [4]. A threshold of an ABI < 0.9 was used to define PAD since it presents a high sensitivity and specificity for the detection of stenosis ≥50% in peripheral arteries [23,24]. Conversely, an ABI ≥ 1.4 was used to define the presence of calcification (NCAD) in the peripheral arteries [9].
A comprehensive questionnaire collected demographic, clinical, and lifestyle variables, including smoking status, physical activity, previous cardiovascular pathology (self-reported and in the clinical history), anthropometric measurements (height, weight, and waist circumference), hypertension, diabetes mellitus, and blood test results (total cholesterol, HDL-cholesterol, triglycerides, and glycemia), among others. The demographic and clinical characteristics of the final sample are described in Table 1.

2.3. Genetic Analyses

Genomic DNA was isolated from the whole peripheral blood by automated DNA purification based into M-PVA Magnetic Bead Technology using the Chemagic Magnetic Separation Module I instrument (Revvity, Inc., Waltham, MA, USA) at the IGTP-HUGTP Biobank (Spanish National Registry of Biobanks B.0000643) integrated in the Platform ISCIII Biobanks and Biomodels (PT20/00050). We determined DNA quality by the absorbance method. The 260/230 nm ratio was 1.7–2.1, and the 260/280 nm ratio was higher, 1.8, demonstrating a good DNA accuracy. Samples were processed following standard operating procedures with the appropriate approval of the Ethical and Scientific Committees.
Fifteen single-nucleotide polymorphisms (SNPs) in 13 different chromosomal regions were genotyped using KASPar assays, an allele-specific PCR technology, by an external genotyping core facility (Progenika Biopharma S.A., Derio, Spain). SNPs were selected based on two criteria: location in a gene or chromosomal region previously associated with PAD or with other related cardiovascular phenotypes and a minor allele frequency (MAF) higher than 10% in the Iberian population from the 1000 Genomes Project (except for the polymorphisms in the APOE gene) to maximize statistical power to detect genetic effects [25]. These included markers in gene regions involved in endothelial function, inflammation and immunity, or lipid metabolism, among other pathways. Details on the selected SNPs (identifier, chromosomal location, MAF, etc.) are described in Table 2.

2.4. Statistical Analyses

Quality control (QC) analyses were performed on the dataset. Forty-one negative controls and nine positive controls (i.e., sample duplicates) were included. Missingness by individual and by marker and deviations from Hardy–Weinberg equilibrium expectations were assessed to check data quality. After QC, the potential presence of population stratification effects between cases and controls was assessed through the calculation of the inflation factor (λ).
Association analyses were conducted by means of logistic regression models for binary traits (i.e., PAD and NCAD) and by linear regression models for continuous variables (i.e., ABI). All analyses compared, in different models, individuals with PAD or with NCAD vs. healthy participants. The analyses estimated effect sizes (odds ratio (OR) or β (β) for logistic and linear regressions, respectively) with their 95% confidence intervals under an additive genetic model. The analyses were conducted without covariates (unadjusted) and including all potential covariates identified in previous studies as influencing the ABI and PAD: sex, age, cardiovascular pathology, hypertension, diabetes, dyslipidemia, smoking, physical activity, central obesity, and abdominal obesity. Significant results were further evaluated under different models of genetic effect (i.e., dominant and recessive) in an attempt to clarify biological effects. To correct the rate of type I error, a false discovery rate (FDR) controlling procedure and a permutation test (1000 permutations) were implemented on significant results. The linkage disequilibrium was quantified between the markers in the 9p21.3 chromosomal region. All QC and association analyses were conducted using PLINK v.1.9 (https://www.cog-genomics.org/plink/1.9/ (accessed on 7 October 2019) [26]), specific for genetic association analyses.
Finally, we used the software QUANTO v.1.2.4 to estimate the minimum size of genetic effect that our cohort has statistical power to detect [27] (https://preventivemedicine.usc.edu/download-quanto/ (accessed on 15 September 2019)). We run post hoc power calculations based on our sample size and the minor allele frequency of the included SNPs and assuming additive effects, a statistical power of 80%, and a two-sided significance threshold of 0.05. Power calculations for the ABI as a quantitative outcome used a mean value of 1.039 and an SD of 0.130.

3. Results

3.1. Sample

From the original 1810 subject dataset, 58 samples were removed for low quality of genetic results, and 4 samples presented low call rate and were also removed. All genetic markers fit Hardy–Weinberg equilibrium expectations (p > 0.1), and all 50 controls (both positive and negative) were successfully identified. After quality control filters, genetic data for 1748 individuals were included, and the genotyping call rate was >99%.
Subsequently, 119 subjects were removed because of discrepancies in data. Twenty-three subjects declared non-European ancestry and were not included to avoid potential population stratification effects in the study. Thus, the final sample included 1606 participants: 112 subjects with PAD (ABI < 0.9), 127 with arterial calcification (ABI > 1.4), and 1367 healthy controls (0.9 < ABI < 1.4) (Table 1).
We did not detect signs of population stratification in our analyses, as estimated by inflation factors (λ) lower than 1.1 (λ of 1.00 and 1.08 for the adjusted and unadjusted analyses, respectively). Power calculations estimated that our cohort was powered to identify alleles than conferred increases in the risk of developing PAD or NCAD (GRR) higher than 1.40–1.43 and increases in ABI higher than 0.010–0.014 units, depending on the minor allele frequency of the analyzed SNPs.

3.2. Association Analyses

We did not find any significant association between the analyzed genetic variants and PAD (Table 2). Moreover, association analyses on the PAD severity, assessed as mild, moderate, or severe PAD (ABI values of 0.7–0.9, 0.4–0.7, and <0.4, respectively [13]), did not find significant associations. In contrast, polymorphisms in the 9p21.3 chromosomal region were significantly associated with the presence of calcified arteries in the lower limbs (Table 3). Minor alleles in SNPs rs10757278 and rs1333049 were significantly associated with a decreased risk of NCAD both in the unadjusted regression models and in the models adjusted for all covariates (p < 0.01 in all cases). Given the differences across groups for some environmental risk factors (Table 1), we further explored if the exclusion of some covariates, specifically cardiovascular pathology and hypertension, from the adjusted models modified the observations, but the association results remained very stable across adjusted, unadjusted, and different partially adjusted models.
In the post hoc analyses, we analyzed different genetic models to clarify the biological mechanism of the observed associations. Thus, we detected a higher degree of association between variants in 9p21.3 and NCAD under a recessive genetic model (results shown in Table 4), but there were not associations under a dominant model (p > 0.10).
Individuals carrying two minor alleles in any of the SNPs explored in locus 9p21.3 had less than half the risk of presenting calcified arteries either under an unadjusted model (OR = 0.45; 95%CI = 0.26/0.77; p = 0.004) or after adjusting for all potential confounding variables (OR = 0.47; 95%CI = 0.27/0.81; p = 0.007). The unadjusted models resisted both a correction for multiple comparisons using an FDR method (p = 0.029) and a permutation test (p = 0.036–0.038), but the adjusted models only resisted the FDR control (p = 0.049) and not the permutation test (p = 0.085).
Association analyses of ABI as a continuous variable showed limited support to the results found in the analyses of binary traits. Under an additive model, no associations were identified between SNPs in the candidate genes and ABI neither when we analyzed individuals with PAD and controls nor among individuals with NCAD and controls (Supplementary Tables S1 and S2). We further conducted an association test between ABI and SNPs in the 9p21.3 region under a recessive genetic model since those were the analyses that showed the highest degree of association with arterial calcification (Table 4). Only SNP rs10757278 was associated with ABI at the standard nominal threshold (p = 0.040 in an unadjusted model and p = 0.053 in an adjusted model), with carriers of two copies of the minor allele (A) showing a lower ABI (β = −0.012; 95%CI: −0.024/−0.001). However, none of the comparisons held a correction for multiple comparisons or a permutation test (p > 0.3 in all cases).

4. Discussion

In this study, we analyzed fifteen polymorphic markers in 13 different genomic regions to identify genetic factors modulating the ABI, a well-known prognostic marker of CVD, and the risk to develop PAD or NCAD. We found that genetic variants at the 9p21.3 chromosomal region are significantly associated with an increased risk of NCAD in what constitutes, to our knowledge, the first genetic locus associated with arterial calcification assessed through ABI.
There is a broad consensus that ABI values higher than 1.4 are good predictors of subclinical atherosclerosis and cardiovascular risk assessment [7]. Other assessments of CV risk, such as CAC, may provide slightly better discrimination of risk, but its routine assessment is not feasible in primary healthcare settings since it is quantified through computed tomography [28,29]. Quick and non-invasive estimators of subclinical atherosclerosis, such as ABI, constitute practical alternatives to be implemented in primary care centers. However, clinical investigations on this complex phenotype have been limited compared with other CVD markers. Moreover, research on its genetic determinants has focused only on the lower tail of the distribution of ABI values, clinically associated to PAD, while systematically ignoring high values, associated to NCAD.
In our study, we conducted separate analyses for both tails of the distribution of ABI values. Thus, we found a remarkable protective effect of the minor allele at both SNPs in the 9p21.3 region (rs10757278-A and rs1333049-G), which are in tight linkage disequilibrium. Homozygotes for the minor allele at any of these two polymorphisms, which account for almost a quarter of our population, showed a risk of developing NCAD that was 40–110% lower than carriers of major alleles (OR = 0.69 and 0.46 for the minor allele under additive and recessive models, respectively).
Chromosome 9p21.3 harbors no known protein-coding genes, but it expresses a long noncoding RNA called ANRIL (antisense noncoding RNA in the INK4 locus) that regulates genes in the atherogenic network [30]. Among other processes, ANRIL participates in the process of phenotypic switching of vascular smooth muscle cells to an adverse proliferative phenotype via microRNAs and the RAS/RAF/ERK signaling pathway (for a recent review, check [31]). The involvement of genetic variants at the 9p21.3 region in the pathophysiology of CVD was first identified by three independent GWASs in 2007 and has since been implicated in other CVD complications such as stroke or aneurysms. In this sense, our results extend to NCAD previous evidence on the implication of this chromosomal region in other endophenotypes of subclinical atherosclerosis, such as the ABI [18,19] and CAC [21]. This novel association, however, cannot be discussed with other reports on the genetics of NCAD, since the present work constitutes the first study published on the topic. So far studies have mostly focused on PAD as a binary outcome, and the few reports that have analyzed genetic determinants of the ABI as a continuous trait have actively excluded individuals with ABI > 1.4 from the analyses to remove heterogeneity from the analyses (e.g., [17,18,19]).
In contrast to coronary heart disease (CHD), very few genetic variants that influence susceptibility to PAD or ABI had been discovered until two recent genome-wide meta-analyses slightly broadened the list of associated loci [19,20]. Summarizing those and other previous works, the only locus that has been consistently associated with both PAD and ABI (limited to values lower than 1.4) is 9p21.3. In our study, we did not find any association between the markers in 9p21.3 or in any of the other 12 genomic regions analyzed and the risk to develop PAD or variation in ABI values. At most, a marginally significant association was observed between an SNP in the 9p21.3 region and ABI values among individuals with NCAD and controls but only under a recessive genetic model, the one with the strongest association with NCAD. Some of the difficulties identifying genetic factors associated to PAD and ABI might be attributed to the lower number of studies invested in the peripheral system and limited sample sizes, until some recent efforts. However, some authors have postulated that it might be also partly due to a greater degree of clinical and genetic heterogeneity in PAD [32]. In this sense, subdivision of precise phenotypes has been proposed as an approach to reduce heterogeneity and identify genetic markers that confer susceptibility to subclinical atherosclerosis [16]. The ABI can be considered an especially sensitive assessment of subclinical atherosclerosis since values at both extremes of its distribution are associated to different clinical manifestations: PAD and NCAD. This characteristic also implies that it needs special heed from an analytical point of view compared with other clinical measurements. Thus, a practical (but limited) approach has led researchers to remove high values from the analyses and focus only on PAD. Here we showed that separate analyses at both extremes of ABI values provide additional information on two different clinical entities related to subclinical atherosclerosis.
Other limitations of our study include the lack of a replication cohort to confirm our findings. Moreover, our analyses focused only on a restricted number of candidate genes instead of a genome-wide design, and thus, our conclusions were limited to a hypothesis-driven approach. Finally, our sample size had limited statistical power to detect small size effects, which tend to be the most common in genetic association studies of complex phenotypes, such as the ABI. Research in the field is starting to overcome this limitation by analyzing data from big cohorts, such as the UK Biobank or the MVP, but has also to be extended into generating joint consortia efforts, especially if multi-omics approaches are going to be implemented. In this sense, promising approaches to understand the pathophysiology of arterial calcification should consider the inclusion of system biology approaches combining at least epigenetic and genetic analyses. In the case of 9p21.3, for example, this locus exerts its functions through the expression of a non-coding RNA, which in turn is regulated by DNA methylation at the promoter region and by existing alleles [31]. Moreover, increasing evidence offers support for a relevant role of specific microRNAs in the pathophysiology, but also as diagnostic tools, of arterial calcification, PAD, and other atherosclerotic complications [33,34].

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes15010002/s1: Table S1: Linear regression results on the ankle brachial index (ABI) in PAD and control individuals; Table S2: Linear regression results on the ankle brachial index (ABI) in NCAD and controls.

Author Contributions

Conceptualization, M.V., M.M., P.T.-M. and M.T.A.; formal analysis, M.V., G.P. and N.L.-V.; investigation, I.C.C. and E.P.; resources, E.P.; data curation, R.F., A.H. and J.M.B.-D.; writing—original draft preparation, M.V. and G.P.; writing—review and editing, M.V., G.P. and M.T.A.; visualization, A.C.-G.; supervision, M.V. and M.T.A.; funding acquisition, M.T.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Carlos III Institute of Health, Ministry of Economy and Competitiveness (Spain), co-funded with the European Union ERDF funds, grant numbers PI11/00765 and PT20/00050. This work was partially supported by ICREA under the ICREA Academia program to M.M.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of the IDIAP Jordi Gol Foundation of Investigation in Primary Care (P11/61).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to consent limitations.

Acknowledgments

The study group would like to thank all the members of the ARTPER cohort for their participation in the study. We also would like to thank the editor and three anonymous reviewers for their contributions to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Moran, A.E.; Forouzanfar, M.H.; Roth, G.A.; Mensah, G.A.; Ezzati, M.; Murray, C.J.L.; Naghavi, M. Temporal Trends in Ischemic Heart Disease Mortality in 21 World Regions, 1980 to 2010. Circulation 2014, 129, 1483–1492. [Google Scholar] [CrossRef] [PubMed]
  2. Chen, S.; Kuhn, M.; Prettner, K.; Bloom, D.E. The Macroeconomic Burden of Noncommunicable Diseases in the United States: Estimates and Projections. PLoS ONE 2018, 13, e0206702. [Google Scholar] [CrossRef] [PubMed]
  3. Shu, J.; Santulli, G. Update on peripheral artery disease: Epidemiology and evidence-based facts. Atherosclerosis 2018, 275, 379–381. [Google Scholar] [CrossRef] [PubMed]
  4. Alzamora, M.T.; Forés, R.; Baena-Díez, J.M.; Pera, G.; Toran, P.; Sorribes, M.; Vicheto, M.; Reina, M.D.; Sancho, A.; Albaladejo, C.; et al. The Peripheral Arterial Disease Study (PERART/ARTPER): Prevalence and Risk Factors in the General Population. BMC Public Health 2010, 10, 38. [Google Scholar] [CrossRef]
  5. Forés, R.; Alzamora, M.T.; Boixadera-Planas, E.; Vázquez, A.; Pera, G.; Torán, P. Evolución de la prevalencia de arteriopatía periférica en la práctica clínica: Estudio descriptivo poblacional con bases de datos reales (SIDIAP-CMBD) [Evolution of the prevalence of peripheral artery disease in clinical practice: A descriptive population study with real databases (SIDIAP-CMBD)]. Aten. Primaria 2022, 54, 102437. [Google Scholar] [CrossRef]
  6. Song, P.; Rudan, D.; Zhu, Y.; Fowkes, F.J.I.; Rahimi, K.; Fowkes, F.G.R.; Rudan, I. Global, regional, and national prevalence and risk factors for peripheral artery disease in 2015: An updated systematic review and analysis. Lancet Glob. Health 2019, 7, e1020–e1030. [Google Scholar] [CrossRef]
  7. Aboyans, V.; Criqui, M.H.; Abraham, P.; Allison, M.A.; Creager, M.A.; Diehm, C.; Fowkes, F.G.R.; Hiatt, W.R.; Jönsson, B.; Lacroix, P.; et al. Measurement and Interpretation of the Ankle-Brachial Index: A Scientific Statement from the American Heart Association. Circulation 2012, 126, 2890–2909. [Google Scholar] [CrossRef]
  8. Nishimura, H.; Miura, T.; Minamisawa, M.; Ueki, Y.; Abe, N.; Hashizume, N.; Mochidome, T.; Harada, M.; Shimizu, K.; Shoin, W.; et al. Clinical Characteristics and Outcomes of Patients with High Ankle-Brachial Index from the IMPACT-ABI Study. PLoS ONE 2016, 11, e0167150. [Google Scholar] [CrossRef]
  9. Velescu, A.; Clara, A.; Martí, R.; Ramos, R.; Perez-Fernandez, S.; Marcos, L.; Grau, M.; Degano, I.R.; Marrugat, J.; Elosua, R. Abnormally High Ankle–Brachial Index Is Associated with All-Cause and Cardiovascular Mortality: The REGICOR Study. Eur. J. Vasc. Endovasc. Surg. 2017, 54, 370–377. [Google Scholar] [CrossRef]
  10. Gu, X.; Man, C.; Zhang, H.; Fan, Y. High ankle-brachial index and risk of cardiovascular or all-cause mortality: A meta-analysis. Atherosclerosis 2019, 282, 29–36. [Google Scholar] [CrossRef]
  11. Alzamora, M.T.; Forés, R.; Torán, P.; Pera, G.; Baena-Díez, J.M.; López, B.; Sierra, M.V.; Fabregat, A. Prevalence of arterial calcification and related risk factors. The multicenter population-based ARTPER study. Gac. Sanit. 2012, 26, 74–77. [Google Scholar] [CrossRef] [PubMed]
  12. Khaleghi, M.; Isseh, I.N.; Bailey, K.R.; Kullo, I.J. Family history as a risk factor for peripheral arterial disease. Am. J. Cardiol. 2014, 114, 928–932. [Google Scholar] [CrossRef] [PubMed]
  13. Firnhaber, J.M.; Powell, C.S. Lower Extremity Peripheral Artery Disease: Diagnosis and Treatment. Am. Fam. Physician 2019, 99, 362–369. [Google Scholar] [PubMed]
  14. Carmelli, D.; Fabsitz, R.R.; Swan, G.E.; Reed, T.; Miller, B.; Wolf, P.A. Contribution of Genetic and Environmental Influences to Ankle-Brachial Blood Pressure Index in the NHLBI Twin Study. National Heart, Lung, and Blood Institute. Am. J. Epidemiol. 2000, 151, 452–458. [Google Scholar] [CrossRef] [PubMed]
  15. Wahlgren, C.M.; Magnusson, P.K.E. Genetic Influences on Peripheral Arterial Disease in a Twin Population. Arterioscler. Thromb. Vasc. Biol. 2011, 31, 678–682. [Google Scholar] [CrossRef] [PubMed]
  16. Rutsch, F.; Nitschke, Y.; Terkeltaub, R. Genetics in Arterial Calcification: Pieces of a Puzzle and Cogs in a Wheel. Circ. Res. 2011, 109, 578–592. [Google Scholar] [CrossRef] [PubMed]
  17. Wassel, C.L.; Lamina, C.; Nambi, V.; Coassin, S.; Mukamal, K.J.; Ganesh, S.K.; Jacobs, D.R.; Franceschini, N.; Papanicolaou, G.J.; Gibson, Q.; et al. Genetic Determinants of the Ankle-Brachial Index: A Meta-Analysis of a Cardiovascular Candidate Gene 50K SNP Panel in the Candidate Gene Association Resource (CARe) Consortium. Atherosclerosis 2012, 222, 138–147. [Google Scholar] [CrossRef]
  18. Murabito, J.M.; White, C.C.; Kavousi, M.; Sun, Y.V.; Feitosa, M.F.; Nambi, V.; Lamina, C.; Schillert, A.; Coassin, S.; Bis, J.C.; et al. Association between Chromosome 9p21 Variants and the Ankle-Brachial Index Identified by a Meta-Analysis of 21 Genome-Wide Association Studies. Circ. Cardiovasc. Genet. 2012, 5, 100–112. [Google Scholar] [CrossRef]
  19. Klarin, D.; Lynch, J.; Aragam, K.; Chaffin, M.; Assimes, T.L.; Huang, J.; Lee, K.M.; Shao, Q.; Huffman, J.E.; Natarajan, P.; et al. Genome-Wide Association Study of Peripheral Artery Disease in the Million Veteran Program. Nat. Med. 2019, 25, 1274–1279. [Google Scholar] [CrossRef]
  20. Van Zuydam, N.R.; Stiby, A.; Abdalla, M.; Austin, E.; Dahlström, E.H.; McLachlan, S.; Vlachopoulou, E.; Ahlqvist, E.; Di Liao, C.; Sandholm, N.; et al. Genome-Wide Association Study of Peripheral Artery Disease. Circ. Genom. Precis. Med. 2021, 14, e002862. [Google Scholar] [CrossRef]
  21. Kavousi, M.; Bos, M.M.; Barnes, H.J.; Lino Cardenas, C.L.; Wong, D.; Lu, H.; Hodonsky, C.J.; Landsmeer, L.P.L.; Turner, A.W.; Kho, M.; et al. Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification. Nat. Genet. 2023, 55, 1651–1664. [Google Scholar] [CrossRef]
  22. Alzamora, M.T.; Baena-Díez, J.M.; Sorribes, M.; Forés, R.; Toran, P.; Vicheto, M.; Pera, G.; Reina, M.D.; Albaladejo, C.; Llussà, J.; et al. Peripheral Arterial Disease Study (PERART): Prevalence and Predictive Values of Asymptomatic Peripheral Arterial Occlusive Disease Related to Cardiovascular Morbidity and Mortality. BMC Public Health 2007, 7, 348. [Google Scholar] [CrossRef] [PubMed]
  23. Guo, X.; Li, J.; Pang, W.; Zhao, M.; Luo, Y.; Sun, Y.; Hu, D. Sensitivity and Specificity of Ankle-Brachial Index for Detecting Angiographic Stenosis of Peripheral Arteries. Circ. J. Off. J. Jpn. Circ. Soc. 2008, 72, 605–610. [Google Scholar] [CrossRef] [PubMed]
  24. Yao, S.T. Pressure Measurement in the Extremity. In Vascular Diagnosis; Bernstein, E.F., Ed.; Mosby: St. Louis, MI, USA, 1993; pp. 169–175. ISBN 0-8016-6557-4. [Google Scholar]
  25. 1000 Genomes Project Consortium; Auton, A.; Brooks, L.D.; Durbin, R.M.; Garrison, E.P.; Kang, H.M.; Korbel, J.O.; Marchini, J.L.; McCarthy, S.; McVean, G.A.; et al. A Global Reference for Human Genetic Variation. Nature 2015, 526, 68–74. [Google Scholar] [CrossRef] [PubMed]
  26. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-Generation PLINK: Rising to the Challenge of Larger and Richer Datasets. GigaScience 2015, 4, 7. [Google Scholar] [CrossRef] [PubMed]
  27. Gauderman, W.J. Sample Size Requirements for Matched Case-Control Studies of Gene-Environment Interaction. Stat. Med. 2002, 21, 35–50. [Google Scholar] [CrossRef]
  28. Geisel, M.H.; Bauer, M.; Hennig, F.; Hoffmann, B.; Lehmann, N.; Möhlenkamp, S.; Kröger, K.; Kara, K.; Müller, T.; Moebus, S.; et al. Comparison of Coronary Artery Calcification, Carotid Intima-Media Thickness and Ankle-Brachial Index for Predicting 10-Year Incident Cardiovascular Events in the General Population. Eur. Heart J. 2017, 38, 1815–1822. [Google Scholar] [CrossRef]
  29. Mahabadi, A.A.; Möhlenkamp, S.; Lehmann, N.; Kälsch, H.; Dykun, I.; Pundt, N.; Moebus, S.; Jöckel, K.-H.; Erbel, R. Heinz Nixdorf Recall Study Investigators CAC Score Improves Coronary and CV Risk Assessment Above Statin Indication by ESC and AHA/ACC Primary Prevention Guidelines. JACC Cardiovasc. Imaging 2017, 10, 143–153. [Google Scholar] [CrossRef]
  30. Hannou, S.A.; Wouters, K.; Paumelle, R.; Staels, B. Functional Genomics of the CDKN2A/B Locus in Cardiovascular and Metabolic Disease: What Have We Learned from GWASs? Trends Endocrinol. Metab. TEM 2015, 26, 176–184. [Google Scholar] [CrossRef]
  31. Razeghian-Jahromi, I.; Karimi Akhormeh, A.; Zibaeenezhad, M.J. The Role of ANRIL in Atherosclerosis. Dis. Markers 2022, 2022, 8859677. [Google Scholar] [CrossRef]
  32. Leeper, N.J.; Kullo, I.J.; Cooke, J.P. Genetics of Peripheral Artery Disease. Circulation 2012, 125, 3220–3228. [Google Scholar] [CrossRef] [PubMed]
  33. Nanoudis, S.; Pikilidou, M.; Yavropoulou, M.; Zebekakis, P. The Role of MicroRNAs in Arterial Stiffness and Arterial Calcification. Front. Genet. 2017, 8, 209. [Google Scholar] [CrossRef] [PubMed]
  34. Vogiatzi, G.; Oikonomou, E.; Deftereos, S.; Siasos, G.; Tousoulis, D. Peripheral Artery Disease: A Micro-RNA-Related Condition? Curr. Opin. Pharmacol. 2018, 39, 105–112. [Google Scholar] [CrossRef] [PubMed]
Table 1. Demographic and clinical characteristics of participants. Data are presented as mean (SD) for continuous variables and as absolute (%) frequencies for binary variables.
Table 1. Demographic and clinical characteristics of participants. Data are presented as mean (SD) for continuous variables and as absolute (%) frequencies for binary variables.
Controls (0.9 < ABI < 1.4) PAD
(ABI < 0.9)
NCAD
(ABI ≥ 1.4)
N1367112127
Age 67.9 (7.6)72.5 (8.3) ***70.2 (8.0) **
Male sex 567 (41%)66 (59%) ***86 (68%) ***
Cardiovascular pathology 1146 (11%)40 (36%) ***17 (13%)
Hypertension 756 (55%)86 (77%) ***72 (57%)
Diabetes276 (20%)42 (38%) ***36 (28%) *
Dyslipidemia878 (64%)84 (75%) *81 (64%)
BMI 29.3 (4.6)29.8 (4.87)30.38 (4.7) *
Abdominal obesity 2850 (63%)71 (64%)81 (64%)
Ever smoker591 (43%)73 (65%) ***66 (52%)
Physical activity 36982 (4828)5402 (3869) ***7278 (5605)
1 Cardiovascular pathology defined as AIM, angina, stroke, or TIA. 2 Abdominal obesity defined as waist circumference ≥102/88 cm (male/female). 3 Physical activity defined as METS expenditure in 14 days performing leisure activity (walking, sport, shopping, cleaning, or gardening). *: p < 0.05; **: p < 0.01; ***: p < 0.001, in the comparisons with controls using Student’s t and chi-square tests for continuous and binary variables, respectively.
Table 2. Logistic regression results on PAD. Adjusted models included as covariates: age, sex, cardiovascular pathology, hypertension, diabetes, dyslipidemia, smoking, physical activity, central obesity, and abdominal obesity. A1/A2: minor/major allele; MAF: minor allele frequency; L95/U95: lower/upper bounds of the 95% confidence interval. Analyses were conducted under an additive model.
Table 2. Logistic regression results on PAD. Adjusted models included as covariates: age, sex, cardiovascular pathology, hypertension, diabetes, dyslipidemia, smoking, physical activity, central obesity, and abdominal obesity. A1/A2: minor/major allele; MAF: minor allele frequency; L95/U95: lower/upper bounds of the 95% confidence interval. Analyses were conducted under an additive model.
Unadjusted ModelAdjusted Model
Gene/LocusChrSNPPositionA1A2MAFnORL95U95pnORL95U95p
TNFRSF91rs24530217929506AG0.3614730.970.721.290.80913881.110.811.510.507
MTHFR1rs180113311796321TC0.4114760.970.731.280.81913900.810.601.110.190
OSBPL103rs190234131754078AG0.4214740.970.731.280.80813880.980.731.330.915
VEGF-A6rs69994743768652AC0.4714770.810.621.070.13413910.820.611.100.176
IL67rs180079522727026CG0.3514721.180.881.560.26613871.170.861.590.317
9p21.39rs1075727822124478AG0.4814510.850.641.110.22813660.910.681.220.516
9p21.39rs133304922125504GC0.4814730.820.631.080.15513870.860.641.160.330
BDNF11rs626527658369TC0.2214781.010.731.410.95113921.040.731.480.845
MMP111rs1799750102799765GT0.5014701.030.791.360.82213841.070.801.450.640
MMP311rs35068180102845216CA0.4714670.810.621.070.14413820.790.581.070.126
CCL217rs102461134252769GA0.2314711.150.841.580.38113861.100.781.550.573
CCL517rs210753835880776TC0.1714760.940.651.370.75913900.940.631.400.772
ICAM-119rs549810285007GA0.5014740.790.601.040.09313880.780.571.060.107
APOE19rs42935844908684CT0.1014650.790.481.300.35713800.730.431.250.252
APOE19rs741244908822TC0.0514730.770.391.510.44213870.810.391.690.568
Table 3. Logistic regression results on NCAD. Adjusted models included as covariates: age, sex, cardiovascular pathology, hypertension, diabetes, dyslipidemia, smoking, physical activity, central obesity, and abdominal obesity. A1/A2: minor/major allele; L95/U95: lower/upper bounds of the 95% confidence interval. Analyses were conducted under an additive model.
Table 3. Logistic regression results on NCAD. Adjusted models included as covariates: age, sex, cardiovascular pathology, hypertension, diabetes, dyslipidemia, smoking, physical activity, central obesity, and abdominal obesity. A1/A2: minor/major allele; L95/U95: lower/upper bounds of the 95% confidence interval. Analyses were conducted under an additive model.
Unadjusted ModelAdjusted Model
Gene/LocusSNPA1A2nORL95U95pnORL95U95p
TNFRSF9rs2453021AG14861.020.771.330.91214671.060.801.410.700
MTHFRrs1801133TC14900.930.711.220.60214690.950.721.240.690
OSBPL10rs1902341AG14881.040.801.360.76114670.990.751.300.934
VEGF-Ars699947AC14920.920.721.190.53614710.950.731.240.725
IL6rs1800795CG14860.910.681.200.48314650.920.691.230.560
9p21.3rs10757278AG14670.700.540.910.00714460.690.530.910.008
9p21.3rs1333049GC14870.690.530.900.00614660.690.530.910.009
BDNFrs6265TC14931.240.921.670.15714721.200.881.630.261
MMP1rs1799750G-14851.270.981.640.07414641.280.981.680.068
MMP3rs35068180-A14811.020.791.330.87014610.980.741.280.861
CCL2rs1024611GA14850.810.591.120.20814640.830.591.160.269
CCL5rs2107538TC14920.940.661.320.70514710.940.651.340.713
ICAM-1rs5498GA14890.850.651.100.22314680.840.641.100.201
APOErs429358CT14810.650.391.070.09114600.660.401.110.116
APOErs7412TC14880.830.441.540.54414670.840.441.620.598
Table 4. Post hoc analyses under a recessive genetic model of associations between 9p21.3 markers and NCAD and ABI between NCAD and controls. Adjusted models included as covariates: age, sex, cardiovascular pathology, hypertension, diabetes, dyslipidemia, smoking, physical activity, central obesity, and abdominal obesity. A1/A2: minor/major allele; MAF: minor allele frequency; L95/U95: lower/upper bounds of the 95% confidence interval.
Table 4. Post hoc analyses under a recessive genetic model of associations between 9p21.3 markers and NCAD and ABI between NCAD and controls. Adjusted models included as covariates: age, sex, cardiovascular pathology, hypertension, diabetes, dyslipidemia, smoking, physical activity, central obesity, and abdominal obesity. A1/A2: minor/major allele; MAF: minor allele frequency; L95/U95: lower/upper bounds of the 95% confidence interval.
NCAD Unadjusted ModelAdjusted Model
Gene/LocusSNPA1A2MAFnORL95U95pnORL95U95p
9p21.3rs10757278AG0.4814670.450.260.770.00414460.470.270.810.007
9p21.3rs1333049GC0.4814870.450.260.780.00414660.470.270.810.007
ABI Unadjusted ModelAdjusted Model
Gene/LocusSNPA1A2MAFnβL95U95pnβL95U95p
9p21.3rs10757278AG0.481467−0.012−0.024−0.0010.0401446−0.012−0.0230.0000.053
9p21.3rs1333049GC0.481487−0.009−0.0200.0020.1101466−0.008−0.0200.0030.149
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Via, M.; Pera, G.; Forés, R.; Costa-Garrido, A.; Heras, A.; Baena-Díez, J.M.; Pedrosa, E.; Clemente, I.C.; Lamonja-Vicente, N.; Mataró, M.; et al. Genetic Variants at the 9p21.3 Locus Are Associated with Risk for Non-Compressible Artery Disease: Results from the ARTPER Study. Genes 2024, 15, 2. https://doi.org/10.3390/genes15010002

AMA Style

Via M, Pera G, Forés R, Costa-Garrido A, Heras A, Baena-Díez JM, Pedrosa E, Clemente IC, Lamonja-Vicente N, Mataró M, et al. Genetic Variants at the 9p21.3 Locus Are Associated with Risk for Non-Compressible Artery Disease: Results from the ARTPER Study. Genes. 2024; 15(1):2. https://doi.org/10.3390/genes15010002

Chicago/Turabian Style

Via, Marc, Guillem Pera, Rosa Forés, Anna Costa-Garrido, Antonio Heras, José Miguel Baena-Díez, Edurne Pedrosa, Inmaculada C. Clemente, Noemí Lamonja-Vicente, Maria Mataró, and et al. 2024. "Genetic Variants at the 9p21.3 Locus Are Associated with Risk for Non-Compressible Artery Disease: Results from the ARTPER Study" Genes 15, no. 1: 2. https://doi.org/10.3390/genes15010002

APA Style

Via, M., Pera, G., Forés, R., Costa-Garrido, A., Heras, A., Baena-Díez, J. M., Pedrosa, E., Clemente, I. C., Lamonja-Vicente, N., Mataró, M., Torán-Montserrat, P., & Alzamora, M. T. (2024). Genetic Variants at the 9p21.3 Locus Are Associated with Risk for Non-Compressible Artery Disease: Results from the ARTPER Study. Genes, 15(1), 2. https://doi.org/10.3390/genes15010002

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop