Genome-Wide Association Study of Genetic Variants Associated with Lower Extremity Amputation Risk in Peripheral Artery Disease

Rajashekar Korutla; Tanisha Garg; Michael P Wilczek; Elsie G Ross; Saeed Amal

doi:10.20944/preprints202602.0764.v1

Results

ClinVar Analysis Results

The targeted GWAS using the ClinVar dataset did not identify any variants reaching genome-wide significance (p < 5×10⁻⁸). However, three variants showed suggestive associations (p < 1×10⁻⁵) as shown in the Table 1.

The genomic inflation factor (λ) was 1.046.

The Manhattan plot in Figure 1 shows peaks of suggestive significance, particularly on chromosomes 6 and 22, while the QQ plot in Figure 2 showed some deviation from the expected distribution at the tail end, suggesting potential associations.

ACAF Analysis Results

The comprehensive GWAS using the ACAF dataset tested 7,784,837 variants in 804 samples. No variants reached genome-wide significance. 35 variants showed suggestive associations (p < 1×10⁻⁵). The genomic inflation factor (λ) was 1.017. The suggestive variants can be seen in Table 2.

Figure 3. Manhattan Plot of ACAF analysis.

Figure 4. QQ Plot of ACAF analysis.

Functional and Regulatory Annotation of Suggestive Variants

Variant Characterization and Genomic Context

Comprehensive functional annotation of the 38 suggestive variants revealed a characteristic GWAS architecture dominated by non-coding variation. Variant Effect Predictor (VEP) analysis demonstrated that 94.7% of variants were non-coding, with intergenic variants comprising the largest category at 44.7% (17/38), followed by intronic variants at 34.2% (13/38), upstream gene variants at 7.9% (3/38), 3' UTR variants at 5.3% (2/38), and one each of 5' UTR variant, synonymous variant, and splice region variant.

For intergenic variants, nearest gene analysis revealed that all variants were located within 500kb of at least one gene, with 18 variants (47.4%) residing within gene boundaries (0 bp distance) and 20 variants (52.6%) located near genes at distances ranging from 3,424 to 51,750 bp. Notably, nine variants clustered near the pseudogene RN7SL776P on chromosome 4 (distances 26-44kb), and five variants were associated with ENSG00000235965 on chromosome 21.

Population Genetics and Evolutionary Context

Population frequency analysis via gnomAD v4 revealed remarkably high allele frequencies across all variants. All 38 variants (100%) were present in gnomAD, with 37/38 (97.4%) being common variants (MAF ≥ 5%). The mean allele frequency was 30.96%, with a median of 29.60%. According to ACMG guidelines, these 37 variants meet the BA1 criterion for stand-alone benign classification based on allele frequency alone.

The sole exception was rs62420977 (chr6:129165442G>A), with an allele frequency of 2.15% in exomes and 1.88% in genomes. Population stratification analysis revealed striking enrichment in the Amish population (13.71%), with particularly high frequency in Amish females (17.23%) compared to males (9.95%). This variant showed highest frequencies in Finnish (2.37%) and Non-Finnish European (2.95%) populations, while being nearly absent in East Asian (0.02%) and having low frequency in African (0.47%) populations.

Pathogenicity Prediction

CADD v1.7 analysis uniformly indicated benign consequences. All 38 variants had PHRED-scaled scores below 10, with a mean of 2.06 ± 2.20 and maximum of 8.66. The highest-scoring variant (chr4:167367142T>C, PHRED 8.66) remained well below the threshold of 20 typically used to identify potentially pathogenic variants. These uniformly low scores as demonstrated in Figure 5 and Figure 6, suggest the variants are unlikely to cause disease through protein-damaging mechanisms.

Splicing Analysis

The splice region variant rs267607236 in MLC1 was evaluated using SpliceAI. Despite its VEP annotation as a splice_region_variant, SpliceAI predicted no significant splicing impact, with all delta scores below 0.1 (maximum Δ = 0.06), well below the 0.2 threshold for predicted splicing effects.

Expression Quantitative Trait Loci Analysis

GTEx v10 analysis was limited by database coverage, with only 7/38 variants (18.4%) present in the database. Among these, three variants (7.9%) demonstrated significant quantitative trait loci effects:

rs6010165 (MLC1): Significant eQTL in whole blood (p = 2.9×10⁻¹³, NES = -0.18), affecting expression of MLC1, MOV10L1, and PANX2. Also showed splicing QTL activity.
rs62420977 (LAMA2): Splicing QTL in thyroid tissue (p = 1.8×10⁻⁵, NES = -0.56).
rs2460694 (B3GAT2): eQTL in tibial nerve (p = 3.5×10⁻⁴, NES = 0.13).

Four additional variants (rs62332335, rs12108386, rs28806580, rs9653673) were present in GTEx but showed no significant eQTL associations. Notably, no variants showed eQTLs in arterial tissues (aorta, coronary, tibial arteries).

Regulatory Element Analysis

VEP identified only one variant, rs58661519, as overlapping a regulatory element (enhancer ENSR10_93T3ZW). However, ENCODE SCREEN analysis found no candidate cis-regulatory elements (cCREs) at any of the three positions queried (chr17:21269385, chr22:50076844, chr11:20368196) across 1,518 biosamples examined.

RegulomeDB analysis was successful for only 15/38 variants (39.5%), with two achieving score 1f (rs1045164134 and rs6010165), indicating strong regulatory evidence including eQTL data. However, the limited coverage suggests many variants occur in poorly characterized regulatory regions.

Chromatin Accessibility Analysis

ENCODE chromatin accessibility data was available for all variants but was limited to five non-vascular cell types: HCT116 cells, activated and resting T-cells, A673 cells, and LPS-treated macrophages. No PAD-relevant vascular tissue data was available, preventing tissue-specific interpretation.

Comparison of Results Across Analyses

Unified Genetic Architecture

The integration of multiple analytical approaches revealed a remarkably consistent genetic architecture underlying PAD amputation risk. Conditional analysis across seven genomic regions containing multiple suggestive variants demonstrated uniform absence of secondary signals. These regions encompassed substantial genomic territory, with variant counts ranging from 682 in the chromosome 17 region to 1,056 in the chromosome 11 region. After conditioning on the lead variant in each region, all association signals were effectively eliminated, with conditional p-values showing dramatic increases typically of several orders of magnitude compared to the original associations.

This pattern held true even for the chromosome 4 region at 167.3-167.4Mb, which contained nine of our suggestive variants representing the highest density in our dataset. The complete elimination of secondary signals across all regions indicates that each associated locus contains a single causal variant, with all other associations merely reflecting linkage disequilibrium. This relatively simple genetic architecture contrasts with many complex traits that show extensive allelic heterogeneity and suggests that identifying causal variants through fine-mapping may be more tractable than initially anticipated.

Critical Impact of Study Design

The alternative control analysis provided compelling validation of our extreme phenotype design, revealing how profoundly study design choices influence genetic discovery. When we compared our matched design utilizing 399 cases and 405 controls against an analysis using all available PAD patients without amputation as controls, comprising 399 cases and 7,139 controls, the results were striking. Despite the 18-fold increase in control sample size, which traditional power calculations would suggest should enhance discovery, we observed a paradoxical and dramatic loss of statistical power. The systematic loss of signal when using all controls versus our extreme phenotype design is visualized in Figure 7, demonstrating how the traditional approach dilutes amputation-specific genetic effects.

Of the 38 variants that achieved suggestive significance in our extreme phenotype design, 36 variants lost this significance threshold when tested with all controls, representing a 94.7% failure rate. The mean p-values increased by 2.17 orders of magnitude, approximately 148-fold, while mean odds ratios attenuated from 1.42 to 1.16, representing a 14.6% reduction in apparent effect size. Several variants showed particularly dramatic deterioration in statistical evidence. The variant rs59299975 on chromosome 11 saw its p-value deteriorate from 2.00×10⁻⁶ to 0.028, a staggering 14,000-fold increase that would lead to its complete dismissal in a conventional analysis. Similarly, rs57745152 shifted from 7.00×10⁻⁶ to 0.026, and rs1834436109 moved from 6.87×10⁻⁶ to 0.012.

Intriguingly, only one variant, rs60516505 on chromosome 19, showed modest improvement with the all-controls design, with its p-value changing from 2.54×10⁻⁶ to 6.44×10⁻⁷. This unique behavior suggests this variant may capture genetic risk for general PAD severity rather than the specific progression to amputation, highlighting how different study designs can capture distinct aspects of disease biology.

Robustness to Demographic Confounding

The stability of our genetic findings to demographic adjustment provided important validation of their biological relevance. Our covariate adjustment analysis successfully tested all 38 candidate variants in 859 individuals with complete phenotype and genotype data, comprising 424 amputation cases and 435 controls. When comparing models incorporating only ancestry principal components against models that additionally included age at enrollment, the results demonstrated remarkable consistency.

The p-values from the baseline and age-adjusted models showed near-perfect correlation with r = 0.999, indicating that age adjustment had minimal impact on the genetic associations. All 38 variants that showed suggestive association in the baseline model maintained this level of significance after age adjustment, with no variant losing its suggestive status. The mean -log₁₀(p-value) across all variants increased slightly from 4.857 in the baseline model to 4.887 in the age-adjusted model, suggesting that accounting for age, if anything, slightly strengthened rather than weakened the genetic associations. The minimal impact of age adjustment on genetic associations is shown in Figure 8, with points clustering tightly along the diagonal.

Individual variants showed minimal change, with the strongest association at rs1256321231 in the HTATIP2 gene improving from p = 1.04×10⁻⁶ to p = 8.89×10⁻⁷ after age adjustment. Among the top associations, the maximum change in p-value was less than one order of magnitude, with most variants showing changes of less than 10% in their -log₁₀(p-value).

While we attempted to additionally adjust for sex, these models failed to converge due to numerical instability, likely reflecting the strong sex imbalance with approximately 62% male representation among amputation cases. However, given the extreme stability observed with age adjustment and the fact that ancestry principal components partially capture population-level sex differences, this limitation is unlikely to substantially impact our conclusions.

Functional Annotation Integration

The comparison of results across multiple functional annotation platforms revealed both the current limitations of functional genomics resources for vascular diseases and the convergent evidence supporting a small subset of our variants. Database coverage showed marked disparities that highlight the fragmented nature of current functional annotation resources. While ENCODE provided chromatin accessibility data for all 38 variants, this apparent complete coverage was misleading as the data consistently came from the same five non-vascular cell types across all queries, including colorectal cancer cells, T-cells, Ewing sarcoma cells, and macrophages. This uniformity prevented any tissue-specific interpretation relevant to peripheral artery disease.

The Genome Aggregation Database provided complete coverage and revealed that 37 of 38 variants were common in human populations with allele frequencies exceeding 5%, definitively excluding rare variant models of disease. GTEx offered more limited coverage at only 18.4%, successfully querying just 7 of 38 variants, but yielded the most functionally relevant findings with three confirmed expression quantitative trait loci. RegulomeDB showed intermediate coverage at 39.5%, annotating 15 variants with two achieving the highest regulatory evidence scores. Most strikingly, ENCODE SCREEN failed to identify any candidate cis-regulatory elements for the specific positions we queried, despite examining 1,518 biosamples.

Despite this limited overlap across platforms, certain variants emerged with multiple lines of supporting evidence. The variant rs6010165 demonstrated convergent support as a blood eQTL with highly significant effects on gene expression, achieved a RegulomeDB score of 1f indicating comprehensive regulatory evidence, and represented a synonymous coding change that might still affect regulation. The variant rs62420977 stood out as both a splicing QTL affecting LAMA2 and the sole low-frequency variant in our set with a minor allele frequency of 2.15%, suggesting potential functional importance. The variant rs58661519 was notable for its location within an enhancer element and its intronic position within ADAM12, a metalloproteinase gene relevant to vascular remodeling.

The complete absence of arterial eQTLs across aorta, coronary, and tibial artery tissues in GTEx represents a critical gap, as these tissues would provide the most direct functional evidence for variants affecting peripheral artery disease. Similarly, the lack of vascular-specific chromatin accessibility data from ENCODE means we may be missing crucial tissue-specific regulatory mechanisms that only manifest in disease-relevant cell types.

Regional Architecture Patterns

The regional association analyses corroborated our conditional analysis findings while revealing interesting patterns of variant clustering. Analysis of the top six loci demonstrated that five regions showed isolated association peaks with minimal flanking signals, consistent with single causal variants driving the associations. The lead variants in these regions rose clearly above background noise with no accompanying signals approaching significance, supporting a simple genetic architecture. Regional association plots for the six most suggestive loci is shown in Figure 9. Each panel shows one locus with ±500kb (or ±1Mb for chr4:167Mb region) flanking the lead SNP (purple diamond). Colors indicate approximate linkage disequilibrium based on physical distance from lead variant. Horizontal lines mark genome-wide (red) and suggestive (orange) significance thresholds.

The LINC02477 locus on chromosome 4 represented the sole exception to this pattern, displaying two variants that reached suggestive significance within approximately 1.6 kilobases of each other. These variants, rs62332335 at position 160,553,934 and rs377194319 at position 160,552,306, initially suggested possible allelic heterogeneity. Similarly, the LINC02752 locus contained two variants (rs59299975 and rs59583539) separated by only 10 base pairs with identical association statistics, likely representing the same signal captured twice or variants in perfect linkage disequilibrium. Our conditional analysis definitively demonstrated that both loci represent single association signals rather than independent effects.

The concentration of variants in certain genomic regions revealed intriguing patterns that warrant further investigation. Nine variants clustered near the pseudogene RN7SL776P on chromosome 4, with distances ranging from 26 to 44 kilobases, raising questions about whether this represents a regulatory hotspot affecting nearby genes or potentially reflects technical artifacts from repetitive sequences. Similarly, five variants associated with ENSG00000235965 on chromosome 21 suggest another potential regulatory domain. The identification of variants affecting long non-coding RNAs, specifically LINC02477 and LINC02752, aligns with growing recognition of lncRNA roles in vascular biology and disease, though their specific functions remain largely uncharacterized.

Discussion

This study represents the first genome-wide investigation specifically targeting genetic determinants of amputation risk in PAD patients. Our approach, employing both targeted GWAS using Clinvar analysis and comprehensive GWAS using ACAF methodologies, has identified several suggestive genetic loci that may contribute to amputation risk in this population.

Interpretation of Key Findings

This first genome-wide investigation of genetic determinants specific to amputation risk in PAD patients identified 38 suggestive loci, though none achieved genome-wide significance. The absence of genome-wide significant findings (p < 5×10⁻⁸) with our limited sample size likely reflects limited statistical power rather than absence of genetic effects, given the observed effect sizes ranging from OR 0.48 to 3.14.

The predominance of non-coding variants at 94.7% aligns with the established GWAS paradigm that regulatory variation, rather than protein-coding changes, drives most complex trait associations. Among our 38 variants, only five showed potential for direct functional impact: one splice region variant (rs267607236 in MLC1), 1 synonymous variant (rs6010165), two 3' UTR variants (rs1565356411 in FAR1 and rs57745152 in NATD1), and one 5' UTR variant (rs1834436109 in TUBB8). However, computational predictions suggested limited functional consequences, with the MLC1 splice variant showing SpliceAI scores below significance thresholds and all variants having CADD scores below 10.

The functional evidence was strongest for three variants confirmed as expression quantitative trait loci. The rs6010165 variant in MLC1 demonstrated a highly significant blood eQTL (p = 2.9×10⁻¹³) affecting multiple genes including MLC1, MOV10L1, and PANX2, highlighting potential systemic inflammatory mechanisms in amputation risk. The rs62420977 variant in LAMA2, encoding a critical vascular basement membrane component, showed splicing QTL activity that could affect protein isoforms important for vascular integrity. Additionally, rs2460694 in B3GAT2 showed a weak eQTL signal in GTEx Nerve - Tibial tissue (p = 3.5×10⁻⁴, NES = 0.13).

Several genes implicated by our analysis have plausible biological roles in PAD progression. ADAM12 encodes a metalloproteinase involved in extracellular matrix remodeling, CUBN is involved in vitamin B12 absorption potentially affecting homocysteine levels, and FAR1 catalyzes fatty alcohol synthesis relevant to lipid metabolism. The identification of protective variants, such as rs58440816 near POGZ with OR=0.56, suggests some genetic factors may confer resilience against progression to amputation.

The remarkably high population frequencies, with 37 of 38 variants showing MAF above 5% and a mean frequency of 31%, definitively exclude Mendelian inheritance patterns. This frequency distribution, combined with uniformly low CADD scores, supports a polygenic architecture where multiple common variants with individually small effects collectively influence amputation risk.

Methodological Considerations

Our extreme phenotype design comparing amputation cases to matched PAD controls proved crucial for discovery. The alternative analysis demonstrated that using all 7,139 available PAD patients as controls, despite providing 18-fold more samples, would have missed 94.7% of our associations. This paradox occurs because including mild disease cases dilutes genetic signals specific to severe outcomes. Mean p-values deteriorated by 2.17 orders of magnitude in the conventional analysis, with some variants showing 14,000-fold increases in p-values.

The uniform absence of secondary signals in conditional analysis across all seven tested regions simplifies interpretation and future fine-mapping efforts. Each region appears to harbor a single causal variant, with other associations reflecting linkage disequilibrium. This contrasts with many complex traits showing multiple independent signals per locus.

The remarkable stability of associations with age adjustment, showing correlation of 0.999 between models, indicates these genetic factors operate independently of chronological age. This is particularly noteworthy given age is a major clinical risk factor for amputation, suggesting these variants influence biological pathways rather than simply correlating with disease duration.

Several technical constraints affected our analyses. Models adjusting for sex failed to converge due to the imbalanced sex distribution. The cross-sectional design prevented analysis of genetic effects on time to amputation. Most critically, functional annotation was severely limited by the scarcity of vascular-specific data in public databases, with no arterial eQTLs found in GTEx despite the vascular nature of PAD, and ENCODE chromatin accessibility data lacking vascular cell types for our queried regions.

Biological and Clinical Implications

Our findings support a model where amputation risk results from modest perturbations across multiple biological systems rather than severe disruption of single pathways. The two strongest eQTLs suggest mechanisms involving systemic inflammation (MLC1 blood eQTL, p=2.9×10⁻¹³) and vascular structural integrity (LAMA2 splicing QTL, p=1.8×10⁻⁵), while a nominal association with B3GAT2 in nerve tissue (p=3.5×10⁻⁴) suggests possible peripheral nerve involvement. However, the limited functional evidence for most variants indicates that current genomic databases inadequately capture vascular-specific regulatory mechanisms.

The high population frequencies of associated variants suggest they are well-tolerated under normal conditions but may lower the threshold for irreversible ischemic damage when combined with environmental stressors and traditional risk factors. The mix of protective and risk variants, with odds ratios ranging from 0.48 to 3.14, indicates complex balancing of genetic factors influencing progression to amputation.

For clinical translation, these findings suggest polygenic risk scores rather than single variant testing will be needed for risk stratification. However, the current variants explain limited phenotypic variance and would require validation before clinical implementation. The population-specific enrichment of rs62420977 in Amish populations (13.7% vs 2.15% in gnomAD exomes) highlights the need for ancestry-aware approaches in precision medicine.

Comparison with Previous Literature

As the first genome-wide association study specifically examining genetic determinants of amputation risk in peripheral artery disease patients, our findings lack direct comparisons in the existing literature. Previous genetic studies of PAD have focused primarily on disease susceptibility rather than progression to severe outcomes such as limb loss.

The most relevant prior work includes the PAD GWAS by Klarin et al. (2019) in the Million Veteran Program, which identified 19 loci associated with PAD susceptibility in over 31,000 cases. However, their phenotype definition included all PAD patients without stratification by amputation status. Similarly, Matsukura et al. (2015) conducted a GWAS in a Japanese population identifying susceptibility loci for PAD, but did not examine amputation as an outcome. Our approach of comparing amputation cases to matched PAD controls without amputation represents a methodological advance in studying genetic factors specific to disease progression rather than initiation.

The absence of overlap between our suggestive loci and previously reported PAD susceptibility variants supports the hypothesis that genetic architecture underlying disease complications differs from that influencing disease susceptibility. This aligns with observations in other complex diseases where progression genetics show distinct patterns from susceptibility genetics.

Our findings of predominantly common variants with modest effect sizes (OR 0.48-3.14) parallel the genetic architecture observed in other cardiovascular traits. The high population frequencies we observed (mean MAF 31%) are consistent with the polygenic model of complex disease, similar to findings in coronary artery disease GWAS where most associated variants are common with small individual effects.

The lack of genome-wide significant findings with our sample size mirrors early GWAS efforts in other vascular diseases. Initial PAD GWAS studies similarly required meta-analyses of multiple cohorts to achieve genome-wide significance, suggesting our sample of 804 individuals represents an important first step requiring expansion through collaborative efforts.

Strengths and Limitations

This study's primary strength lies in its novel phenotype definition specifically targeting amputation risk rather than PAD susceptibility. The extreme phenotype design maximized statistical power from a moderate sample size by comparing phenotypic extremes. Comprehensive multi-platform functional annotation, though revealing limited functional evidence, provided a thorough characterization attempt. The All of Us cohort offers greater ancestral diversity than traditional European-focused studies. Multiple sensitivity analyses, including conditional analysis and alternative control comparisons, confirmed the robustness of our approach.

However, significant limitations must be acknowledged. The sample size of 804 individuals limited power to achieve genome-wide significance or detect variants with smaller effects. No independent replication cohort was used to validate findings. Functional annotation was severely constrained by the absence of vascular-specific data in public databases, with no arterial eQTLs or relevant chromatin accessibility data. We could not adjust for important clinical variables including diabetes, smoking, or medications. The cross-sectional design prevented temporal analysis of progression. Additionally, database coverage was incomplete, with only 18.4% of variants in GTEx and 39.5% in RegulomeDB.

Future Directions

These findings establish a foundation for larger collaborative efforts needed to achieve genome-wide significance and enable robust replication. Functional validation of the two eQTL variants (rs6010165 and rs62420977) in disease-relevant cell types - including vascular cells for LAMA2 and blood/immune cells for MLC1 - represents an immediate priority Development of vascular-specific functional genomics resources, including arterial eQTL datasets and endothelial cell chromatin accessibility maps, is critical for understanding these associations. Longitudinal studies tracking progression from PAD to amputation would clarify temporal relationships between genetic factors and clinical outcomes. Integration with clinical risk factors in predictive models could assess the added value of genetic information for patient stratification.

Conclusions

This comprehensive genome-wide association study represents the first systematic investigation of genetic factors specifically influencing amputation risk in peripheral artery disease patients. Through careful extreme phenotype design, we identified 38 suggestive genetic associations that implicate diverse biological pathways in disease progression. While these variants did not achieve genome-wide significance and show limited functional annotation with current resources, the study design successfully isolated genetic factors specific to amputation risk rather than general PAD susceptibility. The simple genetic architecture with single variants per locus and the age-independent nature of associations provide optimism for future fine-mapping efforts. These results establish that progression to amputation involves genetic factors distinct from disease susceptibility, supporting precision medicine approaches that consider disease stage. Larger collaborative studies, functional validation, and development of vascular-specific genomic resources represent critical next steps in translating these findings to improved clinical outcomes for PAD patients at risk for limb loss.

Funding

None.

Data Availability Statement

The linked genotype and phenotype data analyzed in this study are from the National Institutes of Health's All of Us Research Program and are available through its controlled-access Researcher Workbench. Due to the terms of the Data Use Agreement designed to protect participant privacy and data security, the raw and individual-level data cannot be downloaded or deposited into external public repositories. The full summary-level statistics from our genome-wide association study are available upon reasonable request to the corresponding author. All qualified researchers can apply for access to the same raw dataset by registering and completing the required ethics training via the All of Us Research Hub at https://www.researchallofus.org/.

Acknowledgments

We gratefully acknowledge All of Us participants for their contributions, without whom this research would not have been possible. We also thank the National Institute of Health’s All of Us Research Program for making available the participant data [and/or samples and/or cohort] examined in this study.

Conflicts of Interest

No Conflicts of Interest.

References

Conte, M. S.; Bradbury, A. W.; Kolh, P.; White, J. V.; Dick, F.; Fitridge, R.; Mills, J. L.; Ricco, J.; Suresh, K. R.; Murad, M. H.; Aboyans, V.; Aksoy, M.; Alexandrescu, V.; Armstrong, D.; Azuma, N.; Belch, J.; Bergoeing, M.; Bjorck, M.; Chakfé, N.; Wang, S. Global Vascular Guidelines on the Management of Chronic Limb-Threatening Ischemia. European Journal of Vascular and Endovascular Surgery 2019, 58(1), S1–S109.e33. [Google Scholar] [CrossRef]
Firth, D. Bias reduction of maximum likelihood estimates. Biometrika 1993, 80(1), 27–38. [Google Scholar] [CrossRef]
Fowkes, F. G. R.; Rudan, D.; Rudan, I.; Aboyans, V.; Denenberg, J. O.; McDermott, M. M.; Norman, P. E.; Sampson, U. K.; Williams, L. J.; Mensah, G. A.; Criqui, M. H. Comparison of global estimates of prevalence and risk factors for peripheral artery disease in 2000 and 2010: a systematic review and analysis. The Lancet 2013, 382(9901), 1329–1340. [Google Scholar] [CrossRef]
Qvist, P.; Christensen, J. H.; Vardya, I.; Rajkumar, A. P.; Mørk, A.; Paternoster, V.; Füchtbauer, E.; Pallesen, J.; Fryland, T.; Dyrvig, M.; Hauberg, M. E.; Lundsberg, B.; Fejgin, K.; Nyegaard, M.; Jensen, K.; Nyengaard, J. R.; Mors, O.; Didriksen, M.; Børglum, A. D. The Schizophrenia-Associated BRD1 gene regulates behavior, neurotransmission, and expression of schizophrenia risk enriched gene sets in mice. Biological Psychiatry 2016, 82(1), 62–76. [Google Scholar] [CrossRef]
Goodney, P. P.; Tarulli, M.; Faerber, A. E.; Schanzer, A.; Zwolak, R. M. Fifteen-Year trends in lower limb amputation, revascularization, and preventive measures among Medicare patients. JAMA Surgery 2015, 150(1), 84. [Google Scholar] [CrossRef]
Klarin, D.; Lynch, J.; Aragam, K.; Chaffin, M.; Assimes, T. L.; Huang, J.; Lee, K. M.; Shao, Q.; Huffman, J. E.; Natarajan, P.; Arya, S.; Small, A.; Sun, Y. V.; Vujkovic, M.; Freiberg, M. S.; Wang, L.; Chen, J.; Saleheen, D.; Lee, J. S.; Damrauer, S. M. Genome-wide association study of peripheral artery disease in the Million Veteran Program. Nature Medicine 2019, 25(8), 1274–1279. [Google Scholar] [CrossRef] [PubMed]
Ma, C.; Blackwell, T.; Boehnke, M.; Scott, L. J. Recommended joint and Meta-Analysis strategies for Case-Control Association testing of Single Low-Count variants. Genetic Epidemiology 2013, 37(6), 539–550. [Google Scholar] [CrossRef]
Mackness, M.; Mackness, B. Human paraoxonase-1 (PON1): Gene structure and expression, promiscuous activities and multiple physiological roles. Gene 2015, 567(1), 12–21. [Google Scholar] [CrossRef]
Malik, R.; Chauhan, G.; Traylor, M.; Sargurupremraj, M.; Okada, Y.; Mishra, A.; Rutten-Jacobs, L.; Giese, A.; Van Der Laan, S. W.; Gretarsdottir, S.; Anderson, C. D.; Chong, M.; Adams, H. H. H.; Ago, T.; Almgren, P.; Amouyel, P.; Ay, H.; Bartz, T. M.; Benavente, O. R.; Dichgans, M. Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. Nature Genetics 2018, 50(4), 524–537. [Google Scholar] [CrossRef] [PubMed]
Matsukura, M.; Ozaki, K.; Takahashi, A.; Onouchi, Y.; Morizono, T.; Komai, H.; Shigematsu, H.; Kudo, T.; Inoue, Y.; Kimura, H.; Hosaka, A.; Shigematsu, K.; Miyata, T.; Watanabe, T.; Tsunoda, T.; Kubo, M.; Tanaka, T. Genome-Wide Association study of peripheral arterial disease in a Japanese population. PLoS ONE 2015, 10(10), e0139262. [Google Scholar] [CrossRef] [PubMed]
Mustapha, J. A.; Katzen, B. T.; Neville, R. F.; Lookstein, R. A.; Zeller, T.; Miller, L. E.; Jaff, M. R. Determinants of Long-Term Outcomes and Costs in the Management of Critical Limb Ischemia: A Population-Based Cohort Study. Journal of the American Heart Association 2018, 7(16). [Google Scholar] [CrossRef]
A comprehensive 1000 Genomes–based genome-wide association meta-analysis of coronary artery disease. (2015). Nature Genetics 47(10), 1121–1130. [CrossRef]
Norgren, L.; Hiatt, W.; Dormandy, J.; Nehler; Harris, K.; Fowkes, F. Inter-Society Consensus for the Management of Peripheral Arterial Disease (TASC II). Journal of Vascular Surgery 2007, 45(1), S5–S67. [Google Scholar] [CrossRef]
Jones, S.; Swaminathan, A.; Patel, M.; Vemulapalli, S. Lower extremity amputation in peripheral artery disease: improving patient outcomes. Vascular Health and Risk Management 2014, 417. [Google Scholar] [CrossRef]
Tam, V.; Patel, N.; Turcotte, M.; Bossé, Y.; Paré, G.; Meyre, D. Benefits and limitations of genome-wide association studies. Nature Reviews Genetics 2019, 20(8), 467–484. [Google Scholar] [CrossRef]
Ma, Q.; Jacobi, P. M.; Emmer, B. T.; Kretz, C. A.; Ozel, A. B.; McGee, B.; Kimchi-Sarfaty, C.; Ginsburg, D.; Li, J. Z.; Desch, K. C. Genetic variants in ADAMTS13 as well as smoking are major determinants of plasma ADAMTS13 levels. Blood Advances 2017, 1(15), 1037–1046. [Google Scholar] [CrossRef] [PubMed]
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.; Yanek, L. R.; Van Der Harst, P.; Ferguson, J. F.; Crawford, D. C.; Waite, L. L.; Allison, M. A.; Criqui, M. H.; McDermott, M. M.; Mehra, R.; Murabito, J. M. 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(1), 138–147. [Google Scholar] [CrossRef]
Yang, J.; Weedon, M. N.; Purcell, S.; Lettre, G.; Estrada, K.; Willer, C. J.; Smith, A. V.; Ingelsson, E.; O’Connell, J. R.; Mangino, M.; Mägi, R.; Madden, P. A.; Heath, A. C.; Nyholt, D. R.; Martin, N. G.; Montgomery, G. W.; Frayling, T. M.; Hirschhorn, J. N.; McCarthy, M. I.; Visscher, P. M. Genomic inflation factors under polygenic inheritance. European Journal of Human Genetics 2011, 19(7), 807–812. [Google Scholar] [CrossRef]
Yao, Y. Laminin: loss-of-function studies. Cellular and Molecular Life Sciences 2016, 74(6), 1095–1115. [Google Scholar] [CrossRef] [PubMed]
Fox, C. S.; Golden, S. H.; Anderson, C.; Bray, G. A.; Burke, L. E.; De Boer, I. H.; Deedwania, P.; Eckel, R. H.; Ershow, A. G.; Fradkin, J.; Inzucchi, S. E.; Kosiborod, M.; Nelson, R. G.; Patel, M. J.; Pignone, M.; Quinn, L.; Schauer, P. R.; Selvin, E.; Vafiadis, D. K. Update on prevention of cardiovascular disease in adults with Type 2 diabetes mellitus in light of recent evidence: A scientific statement from the American Heart Association and the American Diabetes Association. Diabetes Care 2015, 38(9), 1777–1803. [Google Scholar] [CrossRef] [PubMed]
Gutierrez, J. A.; Mulder, H.; Jones, W. S.; Rockhold, F. W.; Baumgartner, I.; Berger, J. S.; Blomster, J. I.; Fowkes, F. G. R.; Held, P.; Katona, B. G.; Mahaffey, K. W.; Norgren, L.; Hiatt, W. R.; Patel, M. R. Polyvascular disease and risk of major adverse cardiovascular events in peripheral artery disease. JAMA Network Open 2018, 1(7), e185239. [Google Scholar] [CrossRef]
McDermott, M. M.; Kibbe, M.; Guralnik, J. M.; Pearce, W. H.; Tian, L.; Liao, Y.; Zhao, L.; Criqui, M. H. Comparative effectiveness study of self-directed walking exercise, lower extremity revascularization, and functional decline in peripheral artery disease. Journal of Vascular Surgery 2013, 57(4), 990–996.e1. [Google Scholar] [CrossRef]
Ricco, J.; Bartelink, M.; Björck, M.; Brodmann, M.; Cohnert, T.; Czerny, M.; De Carlo, M.; Debus, S.; Espinola-Klein, C.; Kahan, T.; Kownator, S.; Mazzolai, L.; Naylor, A. R.; Röther, J.; Sprynger, M.; Tendera, M.; Tepe, G.; Venermo, M.; Vlachopoulos, C.; Obiekezie, A. 2017 ESC Guidelines on the Diagnosis and Treatment of Peripheral Arterial Diseases, in collaboration with the European Society for Vascular Surgery (ESVS). European Heart Journal 2017, 39(9), 763–816. [Google Scholar]
Gerhard-Herman, M. D.; Gornik, H. L.; Barrett, C.; Barshes, N. R.; Corriere, M. A.; Drachman, D. E.; Fleisher, L. A.; Fowkes, F. G. R.; Hamburg, N. M.; Kinlay, S.; Lookstein, R.; Misra, S.; Mureebe, L.; Olin, J. W.; Patel, R. A.; Regensteiner, J. G.; Schanzer, A.; Shishehbor, M. H.; Stewart, K. J.; Walsh, M. E. 2016 AHA/ACC Guideline on the Management of Patients with Lower Extremity Peripheral Artery Disease: Executive summary: A report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines. Circulation 2016, 135(12). [Google Scholar] [CrossRef] [PubMed]
Kullo, I. J.; Leeper, N. J. The genetic basis of peripheral arterial disease: current knowledge, challenges, and future directions. Circulation research 2015, 116(9), 1551–1560. [Google Scholar] [CrossRef]
Leeper, N. J.; Kullo, I. J.; Cooke, J. P. Genetics of peripheral artery Disease. Circulation 2012, 125(25), 3220–3228. [Google Scholar] [CrossRef] [PubMed]
Morley, R. L.; Sharma, A.; Horsch, A. D.; Hinchliffe, R. J. Peripheral artery disease. BMJ 2018, j5842. [Google Scholar] [CrossRef] [PubMed]
Korutla, R.; Tedder, D.; Brogan, K.; Milosevic, M.; Wilczek, M. P.; Shehadeh, N.; Shara, N.; Ross, E. G.; Amal, S. Analysis and Visualization of Confounders and Treatment Pathways Leading to Amputation and Non-Amputation in Peripheral Artery Disease Patients Using Sankey Diagrams: Enhancing Explainability. Biomedicines 2025, 13(2), 258. [Google Scholar] [CrossRef]
Barraza-Flores, P.; Bates, C. R.; Oliveira-Santos, A.; Burkin, D. J. Laminin and integrin in LAMA2-Related Congenital Muscular Dystrophy: From Disease to Therapeutics. Frontiers in Molecular Neuroscience 2020, 13. [Google Scholar] [CrossRef]
Brignone, M. S.; Lanciotti, A.; Camerini, S.; De Nuccio, C.; Petrucci, T. C.; Visentin, S.; Ambrosini, E. MLC1 protein: a likely link between leukodystrophies and brain channelopathies. Frontiers in Cellular Neuroscience 2015, 09. [Google Scholar] [CrossRef]
Zhang, M.; Bai, J.; Yuan, H.; Duan, X.; Yu, L.; Li, Y.; Li, K.; Rile, S.; Wang, X.; Wang, H.; Liu, P.; Yan, J.; Wang, C. BRD1 deficiency affects SREBF1-related lipid metabolism through regulating H3K9ac/H3K9me3 transition to inhibit HCC progression. Cell Death and Disease 2025, 16(1). [Google Scholar] [CrossRef] [PubMed]

Figure 1. Manhattan Plot on Clinvar Analysis.

Figure 2. QQ Plot on Clinvar Analysis.

Figure 5. Ranked CADD PHRED scores for 38 PAD-associated variants showing all scores below pathogenicity thresholds (PHRED=15 orange line, PHRED=20 red line), with maximum score of 8.66.

Figure 6. Distribution of CADD pathogenicity scores. Left: Histogram demonstrating all variants cluster below clinical significance thresholds. Right: Box plot showing median PHRED score of 1.23 (IQR: 0.50-3.12) with all variants in the benign range.

Figure 7. Alternative control analysis comparing extreme phenotype design versus all PAD controls. Left: Distribution of p-value differences (Δ-log₁₀(p)) showing 37/38 variants with improved significance using extreme phenotype design (positive values), with mean improvement of 2.17 orders of magnitude. Right: Odds ratio comparison demonstrating systematic attenuation of effect sizes when using all controls, with red points indicating variants that lost significance and the single blue point (rs60516505) showing improved association with all controls.

Figure 8. impact of age adjustment on genetic associations.

Figure 9. Regional association plots for the six most suggestive loci.

Table 1. 3 Suggestive SNPs from Clinvar analysis.

Chr:Position	rsID	Ref>Alt	Beta	P-value	OR	Gene/Nearest Gene	Consequence	Distance (bp)	MAF (%)	eQTL	CADD PHRED
chr6:129165442	rs62420977	G>A	0.88	1.00×10⁻⁶	2.41	LAMA2	intron_variant	0	2.15	Yes	0.4
chr22:50076844	rs6010165	G>A	0.484	2.00×10⁻⁶	1.62	MLC1	synonymous_variant	0	11.4	Yes	0.54
chr22:50076841	rs267607236	T>C	0.484	2.00×10⁻⁶	1.62	MLC1	splice_region_variant	0	10.8	No	0.02

Table 2. 35 Suggestive SNPs from ACAF analysis.

Chr:Position	rsID	Ref>Alt	Beta	P-value	OR	Gene/Nearest Gene	Consequence	Distance (bp)	MAF (%)	eQTL (presence)	CADD PHRED
chr4:160553934	rs62332335	G>C	-0.729	1.07×10⁻⁶	0.48	LINC02477	intron_variant	0	29.8	No	7.09
chr10:126157589	rs58661519	C>T	0.756	1.13×10⁻⁶	2.13	ADAM12	intron_variant	0	15.5	No	5.3
chr11:13732125	rs1565356411	TAA>T	1.128	1.58×10⁻⁶	3.09	FAR1	3_prime_UTR_variant	0	6.4	No	0.27
chr1:151481028	rs58440816	G>C	-0.588	1.80×10⁻⁶	0.56	POGZ	intergenic_variant	21,534	46.2	No	1.3
chr11:11089988	rs59583539	T>C	0.557	2.00×10⁻⁶	1.75	LINC02752	intron_variant	0	12.4	No	7.29
chr11:11089978	rs59299975	T>C	0.557	2.00×10⁻⁶	1.75	LINC02752	intron_variant	0	12.4	No	0.02
chr10:16831081	rs57775857	C>T	0.612	2.06×10⁻⁶	1.84	CUBN	intron_variant	0	33.9	No	2
chr4:167358505	rs10390166	A>T	0.517	2.17×10⁻⁶	1.68	RN7SL776P	intergenic_variant	41,990	30.6	No	4.06
chr4:167355663	rs58299204	G>A	0.515	2.23×10⁻⁶	1.67	RN7SL776P	intergenic_variant	44,832	29	No	1.03
chr19:7375984	rs60516505	A>G	0.519	2.54×10⁻⁶	1.68	ARHGEF18-AS1, ARHGEF18	intron_variant	0	27.4	No	0.08
chr4:160552306	rs377194319	T>TCTAA	-0.691	2.69×10⁻⁶	0.5	LINC02477	intron_variant	0	29.8	No	1.74
chr4:167358166	rs12108386	C>A	0.51	2.82×10⁻⁶	1.67	RN7SL776P	intergenic_variant	42,329	32.5	No	0.68
chr4:167372657	rs17601234	T>A	0.5	3.02×10⁻⁶	1.65	RN7SL776P	intergenic_variant	27,838	29	No	1.11
chr21:44596298	rs17004490	T>C	1.144	3.15×10⁻⁶	3.14	TSPEAR, KRTAP10-7, KRTAP10-6	intron_variant	0	14.9	No	4.69
chr4:167373408	rs111484533	A>ACAGT	0.494	3.43×10⁻⁶	1.64	RN7SL776P	intergenic_variant	27,087	31.1	No	2.88
chr21:19075049	rs143785792	G>A	-0.587	3.45×10⁻⁶	0.56	ENSG00000235965	intergenic_variant	27,365	37.2	No	3.79
chr4:167362902	rs60353699	G>T	0.504	3.73×10⁻⁶	1.66	RN7SL776P	intergenic_variant	37,593	29	No	1.62
chr4:167373717	rs59233454	G>A	0.491	4.13×10⁻⁶	1.63	RN7SL776P	intergenic_variant	26,778	29.1	No	1.77
chr5:133657373	rs59870923	T>C	0.516	4.40×10⁻⁶	1.68	FSTL4	intergenic_variant	44,832	51.8	No	3.21
chr21:19051267	rs1985175965	G>T	-0.534	5.00×10⁻⁶	0.59	ENSG00000235965	upstream_gene_variant	3,583	33.8	No	0.32
chr4:167374290	rs58934032	C>T	0.488	5.09×10⁻⁶	1.63	RN7SL776P	intergenic_variant	26,205	29	No	1.77
chr4:167367142	rs59244366	T>C	0.492	5.52×10⁻⁶	1.64	RN7SL776P	intergenic_variant	33,353	29	No	8.66
chr5:133661151	rs1182945568	T>C	0.515	5.56×10⁻⁶	1.67	FSTL4	intergenic_variant	48,610	52.8	No	0.75
chr17:21269327	rs2508018174	G>A	-0.504	5.76×10⁻⁶	0.6	ENSG00000289453	intergenic_variant	6,570	26.6	No	0.46
chr7:73141070	rs1801996601	C>CCACA	0.959	5.92×10⁻⁶	2.61	SPDYE10	intron_variant	0	15.9	No	0.12
chr17:21269385	rs1045164134	A>AT	-0.502	6.01×10⁻⁶	0.61	ENSG00000289453	intergenic_variant	6,512	24.4	No	1.17
chr6:70891807	rs2460694	T>C	-0.551	6.35×10⁻⁶	0.58	B3GAT2	intron_variant	0	30.1	Yes	0.18
chr10:49442	rs1834436109	AGCTCAGGTGTCCTT>A	-1.096	6.87×10⁻⁶	0.33	TUBB8	5_prime_UTR_variant	0	13.2	No	3.22
chr17:21240761	rs57745152	G>A	0.469	7.00×10⁻⁶	1.6	NATD1, TMEM11-DT	3_prime_UTR_variant	0	10.7	No	0.66
chr21:19056133	rs28806580	A>G	-0.509	7.05×10⁻⁶	0.6	ENSG00000235965	intergenic_variant	8,449	26.4	No	0.78
chr21:19051108	rs9653673	G>A	-0.515	7.26×10⁻⁶	0.6	ENSG00000235965	upstream_gene_variant	3,424	26.3	No	4.62
chr5:133664291	rs56555387	G>A	0.508	7.28×10⁻⁶	1.66	FSTL4	intergenic_variant	51,750	53.5	No	1.95
chr11:20368197	rs1256321231	T>A	0.505	8.68×10⁻⁶	1.66	HTATIP2	intron_variant	0	12.8	No	0.49
chr21:19051807	rs1568913903	T>TTA	-0.503	9.12×10⁻⁶	0.6	ENSG00000235965	upstream_gene_variant	4,123	26.4	No	1.28
chr7:37714274	rs1215147035	T>C	-0.513	9.65×10⁻⁶	0.6	GPR141	intron_variant	0	29.9	No	0.9

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.

Genome-Wide Association Study of Genetic Variants Associated with Lower Extremity Amputation Risk in Peripheral Artery Disease

Abstract

Keywords:

Subject:

Introduction

Methods

Study Population

Patient Cohorts

Genotyping and Quality Control

Statistical Analysis

Bioinformatic Analysis

Variant Annotation and Functional Prediction

Expression and Regulatory Analysis

Statistical Refinement Analyses

Regional Visualization