Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Investigation of genomic and transcriptomic risk factors in clopidogrel response in African Americans

Guang Yang, Cristina Alarcon, Catherine Chanfreau, Norman H. Lee, Paula Friedman, View ORCID ProfileEdith Nutescu, View ORCID ProfileMatthew Tuck, Travis O’Brien, Li Gong, View ORCID ProfileTeri E. Klein, View ORCID ProfileKyong-Mi Chang, View ORCID ProfilePhilip S. Tsao, David O. Meltzer, Million Veteran Program, View ORCID ProfileSony Tuteja, View ORCID ProfileMinoli A Perera
doi: https://doi.org/10.1101/2023.12.05.23299140
Guang Yang
1Department of Pharmacology, Center for Pharmacogenomics, Fienberg School of Medicine, Northwestern University, Chicago IL
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cristina Alarcon
1Department of Pharmacology, Center for Pharmacogenomics, Fienberg School of Medicine, Northwestern University, Chicago IL
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Catherine Chanfreau
2VA Salt Lake City Heath Care System, Salt Lake City, UT
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Norman H. Lee
3Department of Pharmacology and Physiology, George Washington University, 2300 I Street NW, Washington, DC, 20037, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paula Friedman
1Department of Pharmacology, Center for Pharmacogenomics, Fienberg School of Medicine, Northwestern University, Chicago IL
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Edith Nutescu
4Department of Pharmacy Practice and Center for Pharmacoepidemiology and Pharmacoeconomic Research, University of Illinois Chicago, College of Pharmacy, Chicago, IL
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Edith Nutescu
Matthew Tuck
5Washington DC VA Medical Center, Washington, DC and The George Washington University, Washington, DC
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Matthew Tuck
Travis O’Brien
3Department of Pharmacology and Physiology, George Washington University, 2300 I Street NW, Washington, DC, 20037, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Li Gong
6Department of Biomedical Data Science, Stanford University, Stanford, CA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Teri E. Klein
7Department of Biomedical Data Science and Department of Medicine, Stanford University, Stanford, CA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Teri E. Klein
Kyong-Mi Chang
8Corporal Michael J. Crescenz VA Medical Center, Philadelphia, PA
9University of Pennsylvania Perelman School of Medicine, Philadelphia, PA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Kyong-Mi Chang
Philip S. Tsao
10VA Palo Alto Healthcare System and Stanford University, Palo Alto, CA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Philip S. Tsao
David O. Meltzer
11Section of Hospital Medicine, Department of Medicine, University of Chicago, Chicago, IL
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sony Tuteja
8Corporal Michael J. Crescenz VA Medical Center, Philadelphia, PA
9University of Pennsylvania Perelman School of Medicine, Philadelphia, PA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sony Tuteja
Minoli A Perera
1Department of Pharmacology, Center for Pharmacogenomics, Fienberg School of Medicine, Northwestern University, Chicago IL
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Minoli A Perera
  • For correspondence: minoli.perera{at}northwestern.edu
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

Abstract

Clopidogrel, an anti-platelet drug, used to prevent thrombosis after percutaneous coronary intervention. Clopidogrel resistance results in recurring ischemic episodes, with African Americans suffering disproportionately. The aim of this study was to identify biomarkers of clopidogrel resistance in African American patients.

We conducted a genome-wide association study, including local ancestry adjustment, in 141 African Americans on clopidogrel to identify associations with high on-treatment platelet reactivity (HTPR). We validated genome-wide and suggestive hits in an independent cohort of African American clopidogrel patients (N = 823) from the Million Veteran’s Program (MVP) along with in vitro functional follow up. We performed differential gene expression (DGE) analysis in whole blood with functional follow-up in MEG-01 cells.

We identified rs7807369, within thrombospondin 7A (THSD7A), as significantly associated with increasing risk of HTPR (p = 4.56 x 10−9). Higher THSD7A expression was associated with HTPR in an independent gene expression cohort of clopidogrel treated patients (p = 0.004) and supported by increased gene expression on THSD7A in primary human endothelial cells carrying the risk haplotype. Two SNPs (rs1149515 and rs191786) were validated in the MVP cohort. DGE analysis identified an association with decreased LAIR1 expression to HTPR. LAIR1 knockdown in a MEG-01 cells resulted in increased expression of SYK and AKT1, suggesting an inhibitory role of LAIR1 in the Glycoprotein VI pathway. Notably, the CYP2C19 variants showed no association with clopidogrel response in the discovery or MVP cohorts.

In summary, these finding suggest that other variants outside of CYP2C19 star alleles play an important role in clopidogrel response in African Americans.

Key Points SNPs within THSD7A were associated to clopidogrel resistance in African Americans.

THSD7A gene expression is higher in patients with clopidogrel resistance and in endothelial cells carrying the risk variants.

No association to high on therapy platelet reactivity or major adverse cardiac events was seen with CYP2C19 star alleles in either the African American cohort investigated.

Decreased LAIR1 gene expression in blood was associated to risk of HTPR and showed downstream effect on SKY and AKT.

Introduction

Clopidogrel is an anti-platelet agent commonly prescribed to prevent ischemic events in patients with coronary artery disease (CAD), acute coronary syndrome (ACS), and in patients undergoing percutaneous coronary intervention (PCI)1. However, the response to clopidogrel can vary widely between individuals, with both pharmacodynamic response (as measured by P2Y12 reaction units [PRU]) or clinical response showing significant interindividual differences2. High on-treatment platelet reactivity (HTPR),3 has been linked to a greater incidence of major adverse cardiovascular events (MACE) in clinical trials4.

This variability in response to clopidogrel is, at least in part, heritable, with genetic factors accounting for up to 70% of the observed variability5,6. Inadequate clopidogrel response (also known as clopidogrel resistance) is associated with an increased risk of recurrence of MACE, such as myocardial infarction, stroke, and even death7–9.

The first pharmacogenomics study of clopidogrel response published by Mega. et al. showed that carrier of CYP2C19 loss of function (LOF) variants were associated with significantly decreased clopidogrel response by reducing active metabolite of clopidogrel10. The first genome-wide association study (GWAS) of clopidogrel response in an European-ancestry cohort was published by Shuldiner et al. and similarly found that LOF alleles in CYP2C19 (CYP2C19*2) significantly decreased clopidogrel activity, which resulted in increased risk of post-PCI thrombotic events6. This finding has been validated by other studies11–13.

It is important to note that most pharmacogenomic studies of clopidogrel have largely been performed in populations of European ancestry, and it remains unclear whether the findings are generalizable to other populations. Additionally, while CYP2C19 has been shown to explain a significant portion (∼12%) of the variability in clopidogrel response, heritability studies suggest that up to ∼70% of the variance in clopidogrel response is genetic in nature6, indicating that other genetic variants likely play a role. Therefore, there is still a need to identify additional genetic variants that contribute to clopidogrel response to improve the efficacy and safety of clopidogrel treatment.

The African American (AA) population is an admixed population which contains more than one ancestry genome. A study done in 24 USA hospitals showed that African Americans (AAs) have a much higher one-year mortality rate (∼7.2%) compared with European-ancestry (∼3.6%) who were treated with clopidogrel14 due to the poor antiplatelet response15,16. Current GWAS studies often exclude individuals of AA ancestry,17 which can lead to a lack of knowledge regarding population-specific genetic variants that may be associated with clopidogrel response in this population. This knowledge gap may limit our ability to fully understand the genetic factors that influence drug response and to develop personalized treatment strategies for AA patients.

The aim of this study is to identify novel biomarkers of clopidogrel response in AAs through GWAS and transcriptomics analysis followed by in vitro functional analysis of implicated variants and genes.

Materials and Methods

Study populations and sample collection

Participants were recruited through AA Cardiovascular Pharmacogenomics Consortium (ACCOuNT) as previously described18. The study protocol was approved by the University of Chicago, University of Illinois, Chicago, Northwestern University, George Washington University, and the DC Veterans Administration Hospital Institutional Review Boards. All participants provided written consent prior to the start of study procedures. In this study, 154 participants were administered a daily dose of 75 mg of clopidogrel for a minimum of 15 days. Blood samples were collected into sodium citrate-vacutainers and analyzed within 4 hours of collection at Northwestern Memorial Hospital diagnostic laboratory, the Washington, DC Veterans Affairs Medical Center or at George Washington University. Samples were tested using the VerifyNow platform, which provided a P2Y12 reaction units (PRU) value for each sample. Participants with a PRU value of 230 or higher were classified as having high on-treatment platelet reactivity (HTPR). This well-established definition of HTPR has been used in clinical practice and research to identify individuals who may be at increased risk of MACE19. Other blood samples were collected into EDTA-vacutainers for DNA isolation and into PAXgene-vacutainers for RNA isolation.

Genotyping and imputation

DNA was isolated from blood sample using the Gentra Puregene Bood kit from Qiagen (Maryland, USA). DNA samples were genotyped using the Illumina Infinium Multi-Ethnic Global Kit (MEGA) at the University of Chicago Functional Genomics Core using standard protocols. Genotyping outputs were created by Genome Studio using 0.15 GenCall score cutoff.

Standard quality control (QC) procedures were implemented and included: removing the SNPs with an excessive missingness (> 5%) and individuals with missing genotypes (> 5%) or high heterozygosity rates (three standard deviations (SD) away from the mean rate). Hardy-Weinberg equilibrium (HWE) was calculated for all genotyped SNPs. Any SNPs that deviate from HWE were flagged but not removed. Gender misspecification was checked using X chromosome heterozygosity. Any individuals who did not match known sample data were excluded. Identity-by-descent (IBD) checks were performed to identify sample duplicates, contaminated samples, and cryptic relationships. For each pair of samples with estimated IBD coefficients greater than 0.185, only the sample with the highest call rate were retained. Principal component analysis (PCA) was performed to facilitate verification of ancestry and correction for population substructure including HapMap samples from three populations (CEPH (Caucasians), Asians and Yorubas (Africans)) as references.

After the quality Control measures, genotyping data was imputed using the Trans-Omics of Precision Medicine Program (TOPMed) imputation server (https://imputation.biodatacatalyst.nhlbi.nih.gov/#!), using TOPMed imputation panel20. SNPs with minor allele frequencies (MAF) below 0.05 and imputation scores under 0.8 were removed. This resulted in 8,599,814 SNPs used in the GWAS analysis.

Genome-wide association studies (GWAS) analyses

Multivariable logistic regression GWAS was performed to estimate the odds ratios (ORs) and 95% confidence intervals (CI) of each variant for the risk of HTPR (case/control analysis) adjusting for age, gender, principal component 1 (PC1), and PC2 using SNPTEST v2 software. Additionally, multivariable linear regression GWAS was performed to identify the association of SNPs and the change of PRU value as a continuous phenotype.

We also conducted local ancestry inference (LAI) GWAS using Tractor. Briefly, local ancestry at an individual level were estimated used RFMix v221 with a window size of 0.2cM, number of generations equal to 8 and the number of trees to generate per Random Forest set to 100. 1000 Genome African ancestry (YRI: Yoruba in Ibadan, Nigeria), and European ancestry (CEU: Northern Europeans from Utah) populations were used as reference panel. The RFMix v2 ancestry calls were converted to Tractor format, which includes genotype dosage and haplotype count information for African, and European ancestry at each SNP position. The Tractor local ancestry GWAS was then performed using PLINK. This analysis allowed us to use the LAI at each SNPs position as a SNP-based covariate in the GWAS for this admixed cohort. All clinical covariates that showed association to PRU or HTPR were included as covariates in the analysis. The deconvoluted model within TRACTOR was used in this analysis22. Meta-analyses on the deconvoluted AFR, and EUR summary statistics was conducted using METAL using a random-effect model.

SNPs with p-value ≤ 5×10−8 were considered as genome-wide significance, whereas SNPs with p-value ≤ 5×10−6 were considered as suggestive level of significance.

Validation in Million Veterans Program Clopidogrel cohort

To validate our genome-wide suggestive hit, we used data collected from the Million Veteran’s Program (MVP)23 for a cohort of African American patients treated with clopidogrel after a PCI performed at VA catheterization facilities (n=823). Compared to the ACCOuNT cohort, the only clinically relevant phenotype to clopidogrel response available was the Major Adverse Cardiovascular Events (MACE) composite phenotype defined as all-cause death, recurrent MI, and stroke following the index PCI, which was collected from electronic medical records data as previously described.24 We extracted genotypes from the MVP Release 4 imputed genomic data; quality of African ancestry imputation in this version of the MVP dataset was improved using the African Genomic Resources (AGR) reference panel.25 The frequencies of the alleles and the imputation quality for the genomic hits tested are provided in Table Supplement 2. We tested the association of the carrier status of candidate hits (coded as an indicator for carrying two copies of the effect allele versus 1 or 0 copy as the reference) using a Cox proportional hazard model adjusted for age at procedure, type of stent received (bare metal versus drug eluting), presence of chronic kidney disease (yes/no), prior history of PCI (yes/no) and timing of PCI relative to enrollment in the MVP study (before versus after). As a sensitivity analysis we tested the inclusion of an indicator for PCI indication (acute coronary syndrome, or ACS, versus stable ischemic heart disease) in the model. The effect of ACS was found not significant in the model, and we did not retain the variable in the final model. Statistical analyses were performed in R using the packages ‘survival’ and ‘survminer’24.

HUVEC cell culture, sample collection, and analysis of SNP association to THSD7A gene

HUVEC cell lines from one white donor and seven African American (AA) donors were purchased from several vendors: one white and 5 AA from Lonza (Walkersville, MD, USA); 1 AA from ACGT (Manassas, VA, USA); and 1 AA from PromoCell (Heidelberg, Germany).

All HUVEC cell lines were cultured in EGM-Plus medium with supplements (Lonza Walkersville) at 37° C in a 5% CO2/95% air cell culture incubator.

DNA was isolated from these 8 cell lines and sequenced to identify their rs7792678 and rs7811849 genotypes. DNA isolation was done using the PureLink™ Genomic DNA Mini Kit (ThermoFisher). Then, DNA was amplified using the Platinum™ SuperFi II Green PCR Master Mix (ThermoFisher) and the PCR product sequenced by ACGT Inc. (Wheeling, IL). The primers used for both PCR amplification and sequencing of THSD7A were: tgtgtcaaaacccaaaaggc (forward) and tcccaccatcaactttgggta (reverse).

Total protein and total RNA were isolated from cells cultured in 10-cm plates. The cells from half of one plate were used for protein extraction, and the cells from the other half of the same plate were used for total RNA extraction.

Differential gene expression analysis of THSD7A in GEO data

To investigate the association of expression of THSD7A gene and risk of HTPR, the GSE3222626 gene expression dataset was downloaded from NCBI GEO (www.ncbi.nlm.nih.gov/geo).27 The GSE32226 dataset included exon array data in peripheral blood cells from 26 CAD patients treated with clopidogrel (75 mg/day). All samples were normalized by robust multiarray average (RMA) procedure. For this study, we selected PRU value ≥ 230 as clopidogrel non-response group, while those PRU value ≤ 230 as clopidogrel response group (as defined in our study). The GEO2R interactive web tool was used for difference gene expression analysis. The raw CEL data was analysis by limma28 (Linear Models for Microarray Analysis). Only THSD7A gene expression data was extracted.

RNA-sequencing and Quality control analysis (QC)

Blood total RNA was isolated using the PAXgene Blood RNA kit from PreAnalytiX (Qiagen). RNA samples were then depleted of globin mRNA with the GLOBINclear kit (Invitrogen, Thermo Fisher) and RNA quality measured by Bioanalyzer or TapeStation (Agilent, Santa Clara, CA).

Samples with RNA integrity number (RIN) score > 7 were selected for sequencing. TruSeq Stranded mRNA Library Prep Kits, and TruSeq RNA Indexes plate for 96 samples (Illumina, San Diego, CA) were used to prepare libraries for sequencing. The libraries were sequenced using the NovaSeq (Illumina) platform by Novogene (Sacramento, CA) following the manufacturer’s instructions, generating paired-end 50 bp reads with approximately 50 million reads per sample.

Sequencing data quality was assessed using FastQC. The sequencing reads were trimmed to remove adaptor-related sequences and reads with ≥ 3 consecutive bases of aggregate quality score < 20 using Trimmomatic (version 0.39)29. The high-quality reads were then quantified at transcript level against the human reference transcriptome (Homo_sapiens.GRCh38.cdna.all.fa) using Salmon30.

Differential gene expression analysis

The differential gene expression (DGE) analysis was performed using R package “DESeq2”31. Gene expression values were filtered based on expression thresholds of at least 10 reads in at least 20% of samples. Only protein coding genes and lncRNAs were included in the analysis. The expression values for each gene were normalized using regularized log transformation. Surrogate variable analysis (SVA) was used to account for unmeasured covariates, and all relevant surrogate variables were added to DESeq2 and generalized linear regressions were fit for each gene32. Genes with a fold change greater than 2.5 and false discovery rate (FDR) less than 0.05 between HTPR and non-HTPR groups were considered differentially expressed.

shRNA knockdown of LAIR1 in MEG-01 cells

The MEG-01 cells were cultured in RPMI 1640 with 10% fetal bovine serum (FBS) and 2 mM l-glutamine plus 100 U/mL penicillin and 100 μg/mL streptomycin in a 37°C humidified atmosphere containing 5% CO2/95% air. LAIR1 shRNA lentivirus clones were purchased from the RNAi Consortium (Sigma-Aldrich). MEG-01 cells were infected by LAIR1 shRNA lentivirus clones. The pLKO.1-puro non-Mammalian control transduction particle (SHC002V) was used as negative control. Puromycin (2ug/ml) was used to select the infected cells. After selection with puromycin (2 ug/ml) for 4 days, real-time quantitative polymerase chain reaction (RT-qPCR) was performed to test the knockdown efficiency.

RT-qPCR analysis

Total RNA was isolated from MEG-01 cells and HUVECs using the miRNeasy Mini kit (Qiagen, Valencia, CA). Then, RT-qPCR analysis was performed using the Power SYBR® Green RNA-to-CTTM 1-Step Kit (Thermo Fisher, Waltham, MA, USA), in an Applied Biosystems QuantStudio 7 system (Thermo Fisher). PCR primers for human genes LAIR1 (leukocyte associated immunoglobulin like receptor 1), SRC (SRC proto-oncogene, non-receptor tyrosine kinase), SYK (spleen associated tyrosine kinase), and AKT1 (AKT serine/threonine kinase 1) were purchased from Sigma-Aldrich (Saint Louis, MO, USA). The primers sequence for each gene is shown in Table S1. PCR primers for the human THSD7A gene were predesigned (Bio Rad, qHsaCID0007366).

Relative gene expression was quantified using the comparative Ct method (2−ΔΔCT), with glyceraldehyde 3-phosphate de-hydrogenase (GAPDH) used as an internal control.

Western blotting and ELISA analyses

Total protein was extracted using a routine procedure. Briefly, total protein from non-LAIR1 knockdown and LAIR1 knockdown MEG-01 cells was extracted using RIPA buffer containing protease and phosphatase inhibitors. Total protein from HUVECs was extracted using RIPA buffer with protease inhibitors. Protein concentration was determined using a Bradford protein assay kit (ThermoFisher Scientific). Lysates were loaded into NuPAGE™ 4–12%, 1.5 mm, 10 well Mini Protein Bis-Tris gels (ThermoFisher Scientific; NP0335BOX), and run at 100 V for 90 min. Gels were then wet transferred to Immun-Blot® PVDF Membrane (Bio-Rad, Hercules, CA) using NuPAGE™ transfer buffer (ThermoFisher Scientific; NP00061) with 10% methanol for 1 hour and blocked for 1 hour at room temperature in Blotting-Grade Blocker Buffer (Bio-Rad, Hercules, CA). Primary antibodies were diluted with their respective blocking buffers and incubated overnight at 4 °C. Washes were performed with TBS 0.1% Tween-20 (TBST) before the addition of secondary antibody for one hour at room temperature. Washes were performed with 1× TBST before imaging on Odyssey Fc. Protein detection was performed using Image Studio Ver 5.2. Most primary antibodies were purchased from Santa Cruz Biotechnology (Dallas, TX), including LAIR-1, Src, Phospho-Src (Tyr419), ERK ½, Phospho-ERK ½, and Cell Signaling (Danvers, MA), including phospho-Akt (Ser473), Akt, phospho-Syk (Tyr525/526), Syk, and GAPDH. THSD7A protein expression in HUVECs was also confirmed by ELISA (MyBioSource, San Diego, CA).

All values were represented as mean ± standard deviation (SD) from at least three independent experiments. Statistical differences between groups were measured by using Student’s T-test. P-values in all experiments were considered significant at less than or equal to 0.05.

Results

The overview of the study can be seen in Figure 1. This study included 154 subjects from the ACCOuNT cohort with genotype data, of which 99 had corresponding RNA-seq data. After QC for genotyping and RNA-seq, 141 samples passed QC for genotype data (case 38, Controls 103 – Table 2). Among 99 RNA-seq data, 93 included PRU value (22 cases and 71 controls). In addition, we validated our finding in a cohort from 823 African American MVP participants on clopidogrel after a PCI. The characteristics of both patient cohorts are summarized in the Table 1.

Figure 1:
  • Download figure
  • Open in new tab
Figure 1: Overview of Study

Infographic showing the multiple stages of this study including discovery cohort analysis, functional in vitro follow up and validation cohorts for both the genomic (blue) and transcriptomic (green) findings.

View this table:
  • View inline
  • View popup
Table 1. Demographic table
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2. Top SNPs of GWAS analysis of HTPR

GWAS identified association at the THSD7A gene locus

Our logistic model GWAS in AA patients treated with clopidogrel identified one locus that attained genome-wide significance level (p<5* 10−8). The top SNP rs7807369 on chromosome 7 was associated with an increased risk of HTPR (OR: 9.11, 95% confidence interval (CI): 3.66 – 22.67, p = 5.57 * 10−9) (Figure 2, Table 2). Our GWAS association to PRU (as a continuous phenotype) did not identify any genome-wide significant signals. One signal located on chromosome 18 reached the suggestive significance level (p<1×10−6). The top SNP rs2226859 was associated with decreasing the PRU value (β = −0.61, 95%CI: −0.72 - −0.49, P = 8.41*10−7). This SNP is in an intergenic region and is not eQTL to any gene in the GTEx database. Our LAI GWAS of HTPR did not find any genome-wide significant hits. However, 10 SNPs (in 5 loci) reached nominal significance. (Table 3)

Figure 2.
  • Download figure
  • Open in new tab
Figure 2. Miami and QQ plots from GWAS

Miami plot representing the association between single nucleotide polymorphisms (SNPs) and High on Treatment Platelet Reactivity (HTPR) (Top) and the association between SNPs and the change of PRU value (Bottom). The x-axis shows the position of each SNP on the chromosome, while the y-axis displays the association (-log10 P-value) between the SNP and the trait of interest. The red line represents the threshold for genome-wide significance (p = 5*10−8), and the blue line represents the suggestive level of significance (p = 1*10−6). The QQ plot displays the observed p-values of the genome-wide association study (GWAS) summary data on the y-axis, compared to the expected p-values under the null hypothesis of no association on the x-axis. Deviations from the diagonal line indicate departures from the null hypothesis, where points above the line suggest more significant associations than expected, while points below the line suggest fewer significant associations than expected. (A) The GWAS identified one locus (THSD7A on chromosome 7) to be associated with higher risk of HTPR at genome-wide significance level (p < 5* 10−8). (B) QQ plot for GWAS. The genomic control value for the GWAS was 1.06 for logistic regression model GWAS (Left) and 1.00 for linear regression model GWAS.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 3: Top SNPs of LAI GWAS analysis of HTPR

The regional plot shows that the top SNP and variants in high linkage disequilibrium (LD) from the logistic regression analysis are located at the intronic region of THSD7A gene. Based on the output of UCSC Genome browser, two SNPs are in high LD (r2 ≥ 0.8) with rs7807369 and are located within a potential enhancer/promoter region (Figure 3). Hi-C interactions analysis from Vanno Portal corroborated these findings with SNPs rs7792678 and rs7811849 affect the THSD7A expression in Human umbilical vein endothelial (HUVEC) cell line by regulating the activation of the promoter (Supplement Figure 1). To investigate the association of THSD7A gene expression to the risk of HTPR, we compared the difference in THSD7A gene expression using public gene expression data (GSE32226)26 comparing HTPR and non-HTPR groups (as defined in our GWAS cohort). From this analysis, higher expression of THSD7A was associated with higher risk of HTPR (p = 0.004) (Figure 4A).

Figure 3.
  • Download figure
  • Open in new tab
Figure 3. Regional plot of top SNP

The colors of the circles refer to linkage disequilibrium (LD) (r2) between top SNP rs7807369 (purple diamond) and nearby SNPs (based on pairwise r2 values from the 1000 Genomes Project reference panel). The blue line and right y-axis show the estimated recombination rate. The x-axis represents the genomic position in chromosome 7 and the left y-axis represents the -log10 p of association with HTPR in discovery cohort. Region plot showed THSD7A variant rs7792678 is high LD with top SNP rs7807369. Both are located at the intronic region of THSD7A gene. SNPs rs7792678 and rs7811849 located at the potential enhancer/promoter region based on ENCODE annotation displayed in the UCSC Genome Browser dataset.

Figure 4.
  • Download figure
  • Open in new tab
Figure 4.
  • Download figure
  • Open in new tab
Figure 4. THSD7A expression associated with higher risk of HTPR

(A) The violin box plot of the expression of THSD7A gene in HTPR (PRU >230) and non-HTPR groups from the GEO dataset GSE32226. The THSD7A gene expression is significant higher in HTPR group compared with non-HTPR group (p = 0.004). (B) Fold-change in expression of THSD7A gene in carriers of the effect alleles compared with carriers wild-type alleles in HUVEC cell cultures. Bar heights indicate the mean expression of THSD7A gene in groups of white wild-type, AA wild-type, and AA effect alleles carriers. Error bars indicate 95% confidence interval estimates of the mean expressions. (*) = p<0.05, (**) = p<0.01, (***) = p<0.001, (****) = p<0.0001, ns = non-significant. (C) Quantified difference in protein expression from the THSD7A ELISA show in the bar chart for the same groups.

To investigate the potential function of rs7792678 and rs7811849 on THSD7A gene expression, we sequenced seven HUVEC (primary endothelial) cell cultures from AA donors and one from a European Ancestry donor. The European and 2 of the AA HUVEC cultures were wild type for the associated haplotype, while 5 AA HUVEC cultures were heterozygous for the associated haplotype. RT-qPCR found that the cultures that carried the effect haplotype showed increased expression of THSD7A (2.2-fold compared with the wild type, p = 0.0006) (Figure 4B). Using ELISA, we showed that THSD7A protein levels significantly increased by a 50% (p = 0.004) in the cultures carrying the effect haplotype as compared with the protein level in the cells with the wild-type haplotype (Figure 4C).

External Validation of top association in the MVP cohort

We examined the association of the top hits with the occurrence of MACE in a sample of AA MVP participants receiving clopidogrel after a PCI. We restricted this analysis to 823 AA participants under the age of 65 at time of the PCI due to loss of outcomes data in older individuals who are more likely to receive care outside of the VA healthcare system. Distribution of genotypes and MACE events are presented in Table S3. We performed a time-to-event analysis for the occurrence of MACE using a Cox proportional hazard model adjusted for age at procedure, type of stent used, existence of chronic kidney disease, prior history of PCI and timing of PCI and enrollment in the MVP study. We were able validate rs191786 (found in the LAI GWAS of HTPR, Table 3, OR = 1.85, p = 6.21×10−6), with a significantly elevated risk of MACE for the homozygous carriers of the C allele in the MVP AA cohort (OR = 1.77, p = 0.049; Table 4). No significant association was found for rs11983224. Homozygous carriers of the C allele for rs1149515 were found at significantly lower risk of MACE (OR=0.54, p=0.0179). Of note the association of CYP2C19 loss-of-function with MACE (p=0.738) was not replicated in this sample of African American patients in contrast to many of the European studies to date (Table 4)

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 4: Validation of Top SNPs in the MVP cohort using an adjusted Cox model for the risk of Major Adverse Cardiac Event (MACE) among African American patients receiving clopidogrel after a percutaneous coronary intervention (PCI).

Low expression of LAIR1 was associated with an increased risk of HTPR

The differential gene expression analysis using RNA-seq data from blood samples identified LAIR1 as associated with HTPR. Lower expression of LAIR1 was associated with an increased risk of HTPR (Figure 5). To assess the role of LAIR1 in clopidogrel response, shRNA knockdown of LAIR1 by lentivirus shRNA clones was performed in MEG−01 cell (megakaryocyte cell line). The expression of LAIR1 was significantly reduced with LAIR1 knockdown. The knockdown of LAIR1 significantly increased the expression of SYK, AKT1 but not SRC. (Figure 6).

Figure 5.
  • Download figure
  • Open in new tab
Figure 5. Volcano plot of the differential gene expression analysis

The volcano plot illustrates the statistical significance and fold change for each gene in the RNA-seq dataset from whole blood. Each point represents a gene, with the x-axis representing the log2 fold change and the y-axis representing the negative logarithm of the p-value. The red color point represents the gene with significant p-value (< 4*10−6) and log2 fold change.

Figure 6.
  • Download figure
  • Open in new tab
Figure 6. RT-qPCR

Fold-change expression of target genes LAIR1 (A), SRC (B), SYK (C), and AKT (D), between MEG−01 cells after LAIR1 knockdown and wild-type cells. Bar heights indicate the mean expression of relative genes. Error bars indicate 95% confidence interval estimates of the mean expressions. (*) = p<0.05, (**) = p<0.01, (***) = p<0.001, (****) = p<0.0001, ns = non-significant. Lack of LAIR1 expression increased the expression of SYK and AKT significantly with no change in SRC gene expression.

The effect of LAIR1 knockdown in protein expression was investigated by Western blot analysis. The results showed that the LAIR1 knockdown cells exhibited a significant decrease in LAIR1 protein compared with the control group. Moreover, the knockdown of LAIR1 led to an increase in the protein expression levels of Syk and Akt. There is no significant difference for protein level of SRC between LAIR1 knockdown and control groups. Taken together, these data demonstrate that LAIR1 inhibits Syk expression and decreases the downstream AKT pathway (Figure 7).

Figure 7.
  • Download figure
  • Open in new tab
Figure 7.
  • Download figure
  • Open in new tab
Figure 7. LAIR1 knockdown cells increases AKT pathway.

(A) Whole-cell lysates of MEG−01 from wild-type and LAIR1 knockdown groups were Western blotted with anti-LAIR1, anti-SRC, anti-p-SRC, anti-Syk, anti-p-Syk, anti-Akt, anti-p-Akt, and GAPDH antibodies. Representative images are shown. (B) Quantification of western blotting images. Bar heights indicate the mean value of relative proteins. Error bars indicate 95% confidence interval estimates of the mean values. (*) = p<0.05, (**) = p<0.01, (***) = p<0.001, (****) = p<0.0001, ns = non-significant. The LAIR1 knockdown MEG−01 cells show significant lower LAIR1 protein level compared with wild-type MEG−01 cells. The LAIR1 knockdown cells increased the phosphorylation levels of Syk and Akt proteins.

Discussion

This is the first GWAS study to investigate biomarkers of clopidogrel response in AAs. We identified a potential functional variant, rs7792678, on the potential enhancer/promoter region that is associated with a higher risk of HTPR. According to the gnomAD browser, the frequency of the effect allele for rs7792678 is 0.414 among AA populations, compared to 0.135 among European populations. This may explain the finding of that AAs have a higher clopidogrel resistance compared to Europeans33,34. Within the ICPC European clopidogrel study this SNP was found at a frequency of 0.105 and was not associated with PRU.35 This study did not investigate genetic associations to HTPR. While we were not able to validate the association of rs7792678 to MACE in the AA MVP cohort, we did find an association to a second SNP (rs1149515) at a nearby locus in Chromosome 7 that did show association to MACE in the AA MVP cohort (Table 5). rs1149515 is in moderate LD with an eQTL to the TSPAN13, however no known function has been found for this variant.

Using HEUVEC cells derived from AAs, we were able to validate the effect of the C allele of rs7792678 with higher expression of thrombospondin 7A (THSD7A). Additionally, the gene expression data from within a multiethnic cohort of clopidogrel patients showed higher expression of THSD7A in patient with HTPR as compared with the non-HTPR group. These validation efforts point to a potential role for THSD7A in clopidogrel response. THSD7A has been implicated in several disease process, including venous thromboembolism36, and CAD37, both of which are linked to platelet function. THSD7A encodes a protein that play a crucial role in endothelial cell migration and angiogenesis38. Studies have shown that THSD7A can inhibit primary endothelial cell migration and tube formation39, which would result in the dysfunction of the endothelial cells. Previous studies have shown endothelial dysfunction was significantly associated with HTPR in patients who were treated with clopidogrel40–42.

Notable, we did not detect an association of clopidogrel response with the CYP2C19*2 allele. CYP2C9*2 was present in our cohort (MAF = 0.177). However, it showed no association to either HTPR or PRU (P = 0.60 and 0.62). CYP2C19*2 carriers, a LOF variant of CYP2C19, were reported to have a decreasing response to clopidogrel in European and Asian populations6,43. While it is clear that the CYP2C19*2 variant has functional consequences on CYP2C19 enzyme function44, our lack of association suggests that additional factors (including additional SNPs) may modify this effect in African ancestry populations. More in-depth sequencing is needed to better understand the role of this clinically actionable SNP on clopidogrel response in African ancestry populations, who carry genetic variation not found in other populations. Moreover, the Cox regression analysis also showed no association of CYP2C19 LOF alleles to MACE in this independent African American cohort. More in-depth sequencing is needed to better understand the role of this clinically actionable SNPs on clopidogrel response in African ancestry populations.

This study identified LAIR1 expression is associated with a decreased risk of HTPR. LAIR1 is a transmembrane glycoprotein that is expressed on hematopoietic and megakaryocytes cells (the precursor cell to platelets)45,46. GPVI is a platelet specific collagen receptor that plays a crucial role in platelet activation and aggregation47. Both LAIR1 and the gene encoding GPVI (glycoprotein VI) are located on chromosome 19 – making up the leukocyte receptor complex (LRC)48. The collagen-binding site on LAIR1 is similar to that of GPVI, a well-known platelet receptor that initiates platelet activation and aggregation45. Unlike GPVI, LAIR1 is an inhibitory receptor that binds collagen49. A recent study by Smith et al.50 demonstrated that mice with LAIR1 deficiency showed a 25% increase in platelet counts, a longer platelet half-life, and increased platelet formation. Our study provides evidence that knockdown of LAIR1 in human megakaryoblast cells results in increased activation of the GPVI pathway. Specifically, we found that SYK and AKT gene and protein expression levels were significantly increased with LAIR1-knockdown. SYK and AKT are two important factors in the platelet aggregation pathway.51 These findings suggest that LAIR1 and GPVI may co-regulate platelet aggregation in human cells. Based on the abovementioned evidence, LAIR1 is a promising target gene that regulates platelet aggregation.

Limitations

One of the major limitations of this study is the small sample size used for analysis. The limited sample size may have reduced the statistical power to detect small to medium effect sizes and rare variants, potentially limiting the generalizability and precision of the findings. We have tried to remedy this limitation through functional validation of our finding in in vitro models as well as in independent data sets. The second limitation is the lack of replication study to validate our GWAS result. Replication is unfortunately not possible as no other independent AA clopidogrel dataset exists. The ICPC cohort is a European cohort which showed that most of our significant SNPs were below 5% allele frequency in their cohort and hence not included in the analysis. Of the SNPs that were also tested in this cohort, none were associated to PRU. They may be due to differences in LD structure between Europeans and African populations as well as differences in phenotyping (no HTPR associations were tested in the ICPC cohort). Future studies with larger AA and multi-ancestry cohorts are warranted to validate and expand upon the findings reported here.

In conclusion, we have performed the first GWAS of clopidogrel response in an African American cohort. These data will be useful as we move toward in the implementation of genotypes-guided clopidogrel therapy in diverse clinical populations. Our work shows that the effect of known variants differs by population and potentially by ancestry.

Data Availability

Genomic and de-identified clinical data for the ACCOuNT participants will be shared on the GEN3 platform at https://gen3.datacommons.io/discovery/ACCOuNT_Clopidogrel_Arm through a control access application. All relevant genetic data for the MVP cohort are available in the main paper and supplemental data. Individual data cannot be shared publicly according to the Data Access Policy of the Million Veteran Program in the VA Office of R&D in Veterans Health Administration.

Sources of Funding

Funding/support: Funding for ACCOuNT was provided by 1U54MD010723−01 from NIHMD (G.Y, C.A, N.H.L, E.N, M. T, T.O, L.G, T. K, D.M, M. A. P). This research is based on data from the Million Veteran Program, Office of Research and Development, Veterans Health Administration, and was supported by award #I01-BX003362 for MVP−003/028 to K.M.C., P.S.T. This work was supported using resources and facilities of the Department of Veterans Affairs (VA) Informatics and Computing Infrastructure (VINCI), VA HSR RES 13-457.

Disclaimer

This publication does not represent the views of the Department of Veterans Affairs or the United States Government.

Data Sharing

Genomic and de-identified clinical data for the ACCOuNT participants will be shared on the GEN3 platform at https://gen3.datacommons.io/discovery/ACCOuNT_Clopidogrel_Arm through a control access application. All relevant genetic data for the MVP cohort are available in the main paper and supplemental data. Individual data cannot be shared publicly according to the Data Access Policy of the Million Veteran Program in the VA Office of R&D in Veterans Health Administration.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S1. Primers for RT-qPCR
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S2: Frequencies of alleles and quality of imputation of the Top SNPs in the MVP dataset (Release 4)
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S3: Distribution of genotypes and MACE rate for the Top SNPs in the MVP sample of African American patients receiving clopidogrel after a PCI

Acknowledgements

Footnotes

  • ↵* A list of members and their affiliations appears in the Supplementary Information.

References

  1. 1.↵
    Nguyen TA, Diodati JG, Pharand C. Resistance to clopidogrel: a review of the evidence. J Am Coll Cardiol. Apr 19 2005;45(8):1157–64. doi:10.1016/j.jacc.2005.01.034
    OpenUrlFREE Full Text
  2. 2.↵
    Nagareddy P, Smyth SS. Inflammation and thrombosis in cardiovascular disease. Curr Opin Hematol. Sep 2013;20(5):457–63. doi:10.1097/MOH.0b013e328364219d
    OpenUrlCrossRefPubMed
  3. 3.↵
    Spiliopoulos S, Pastromas G. Current status of high on-treatment platelet reactivity in patients with coronary or peripheral arterial disease: Mechanisms, evaluation and clinical implications. World J Cardiol. Dec 26 2015;7(12):912–21. doi:10.4330/wjc.v7.i12.912
    OpenUrlCrossRef
  4. 4.↵
    Sharifi H, Habibi V, Emami Zeydi A. Platelet reactivity unit (PRU) in patients undergoing elective PCI: Rethinking the optimal cut point. Anatol J Cardiol. Aug 2017;18(2):164. doi:10.14744/AnatolJCardiol.2017.7939
    OpenUrlCrossRef
  5. 5.↵
    Lewis JP, Shuldiner AR. Clopidogrel pharmacogenetics: Beyond candidate genes and genome-wide association studies. Clin Pharmacol Ther. Mar 2017;101(3):323–325. doi:10.1002/cpt.519
    OpenUrlCrossRef
  6. 6.↵
    Shuldiner AR, O’Connell JR, Bliden KP, et al. Association of cytochrome P450 2C19 genotype with the antiplatelet effect and clinical efficacy of clopidogrel therapy. JAMA. Aug 26 2009;302(8):849–57. doi:10.1001/jama.2009.1232
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    Geisler T, Langer H, Wydymus M, et al. Low response to clopidogrel is associated with cardiovascular outcome after coronary stent implantation. Eur Heart J. Oct 2006;27(20):2420–5. doi:10.1093/eurheartj/ehl275
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.
    Zhang YJ, Li DJ, Li ZY, et al. Influence of genetic polymorphisms in P2Y12 receptor signaling pathway on antiplatelet response to clopidogrel in coronary heart disease. BMC Cardiovasc Disord. Dec 30 2022;22(1):575. doi:10.1186/s12872-022-02988-w
    OpenUrlCrossRef
  9. 9.↵
    Wisman PP, Roest M, Asselbergs FW, et al. Platelet-reactivity tests identify patients at risk of secondary cardiovascular events: a systematic review and meta-analysis. J Thromb Haemost. May 2014;12(5):736–47. doi:10.1111/jth.12538
    OpenUrlCrossRefPubMed
  10. 10.↵
    Mega JL, Close SL, Wiviott SD, et al. Cytochrome p-450 polymorphisms and response to clopidogrel. N Engl J Med. Jan 22 2009;360(4):354–62. doi:10.1056/NEJMoa0809171
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    Pereira NL, Rihal CS, So DYF, et al. Clopidogrel Pharmacogenetics. Circ Cardiovasc Interv. Apr 2019;12(4):e007811. doi:10.1161/CIRCINTERVENTIONS.119.007811
    OpenUrlCrossRefPubMed
  12. 12.
    Stone GW, Witzenbichler B, Weisz G, et al. Platelet reactivity and clinical outcomes after coronary artery implantation of drug-eluting stents (ADAPT-DES): a prospective multicentre registry study. Lancet. Aug 17 2013;382(9892):614-23. doi:10.1016/S0140-6736(13)61170-8
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    Lee CR, Luzum JA, Sangkuhl K, et al. Clinical Pharmacogenetics Implementation Consortium Guideline for CYP2C19 Genotype and Clopidogrel Therapy: 2022 Update. Clin Pharmacol Ther. Nov 2022;112(5):959–967. doi:10.1002/cpt.2526
    OpenUrlCrossRef
  14. 14.↵
    Cresci S, Depta JP, Lenzini PA, et al. Cytochrome p450 gene variants, race, and mortality among clopidogrel-treated patients after acute myocardial infarction. Circ Cardiovasc Genet. Jun 2014;7(3):277–86. doi:10.1161/CIRCGENETICS.113.000303
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    Pendyala LK, Torguson R, Loh JP, et al. Racial disparity with on-treatment platelet reactivity in patients undergoing percutaneous coronary intervention. Am Heart J. Aug 2013;166(2):266–72. doi:10.1016/j.ahj.2013.04.008
    OpenUrlCrossRefPubMed
  16. 16.↵
    Long Y, Gracely EJ, Newschaffer CJ, Liu L. Analysis of the prevalence of cardiovascular disease and associated risk factors for European-American and African-American populations in the state of Pennsylvania 2005-2009. Am J Cardiol. Jan 01 2013;111(1):68–72. doi:10.1016/j.amjcard.2012.08.045
    OpenUrlCrossRefPubMed
  17. 17.↵
    Coram MA, Fang H, Candille SI, Assimes TL, Tang H. Leveraging Multi-ethnic Evidence for Risk Assessment of Quantitative Traits in Minority Populations. Am J Hum Genet. 10 05 2017;101(4):638. doi:10.1016/j.ajhg.2017.09.005
    OpenUrlCrossRef
  18. 18.↵
    Friedman PN, Shaazuddin M, Gong L, et al. The ACCOuNT Consortium: A Model for the Discovery, Translation, and Implementation of Precision Medicine in African Americans. Clin Transl Sci. May 2019;12(3):209–217. doi:10.1111/cts.12608
    OpenUrlCrossRef
  19. 19.↵
    Collet JP. Current Concepts in the Clinical Utility of Platelet Reactivity Testing. Interv Cardiol. Aug 2013;8(2):100–106. doi:10.15420/icr.2013.8.2.100
    OpenUrlCrossRef
  20. 20.↵
    Kowalski MH, Qian H, Hou Z, et al. Use of >100,000 NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium whole genome sequences improves imputation quality and detection of rare variant associations in admixed African and Hispanic/Latino populations. PLoS Genet. 12 2019;15(12):e1008500. doi:10.1371/journal.pgen.1008500
    OpenUrlCrossRefPubMed
  21. 21.↵
    Maples BK, Gravel S, Kenny EE, Bustamante CD. RFMix: a discriminative modeling approach for rapid and robust local-ancestry inference. Am J Hum Genet. Aug 08 2013;93(2):278–88. doi:10.1016/j.ajhg.2013.06.020
    OpenUrlCrossRefPubMed
  22. 22.↵
    Atkinson EG, Maihofer AX, Kanai M, et al. Tractor uses local ancestry to enable the inclusion of admixed individuals in GWAS and to boost power. Nat Genet. 02 2021;53(2):195–204. doi:10.1038/s41588-020-00766-y
    OpenUrlCrossRefPubMed
  23. 23.↵
    Gaziano JM, Concato J, Brophy M, et al. Million Veteran Program: A mega-biobank to study genetic influences on health and disease. J Clin Epidemiol. Feb 2016;70:214–23. doi:10.1016/j.jclinepi.2015.09.016
    OpenUrlCrossRefPubMed
  24. 24.↵
    Chanfreau-Coffinier C, Anglin-Foote T, Lee KM, et al. Abstract 16710: CYP2C19 Polymorphisms and Clinical Outcomes Following Percutaneous Coronary Intervention (PCI) in the Million Veterans Program (MVP). Circulation. 2020;142(Suppl_3):A16710–A16710. doi:doi:10.1161/circ.142.suppl_3.16710
    OpenUrlCrossRef
  25. 25.↵
    O’Connell J, Yun T, Moreno M, et al. A population-specific reference panel for improved genotype imputation in African Americans. Communications Biology. 2021/11/05 2021;4(1):1269. doi:10.1038/s42003-021-02777-9
    OpenUrlCrossRef
  26. 26.↵
    Luchessi AD, Silbiger VN, Hirata RD, et al. Pharmacogenomics of anti-platelet therapy focused on peripheral blood cells of coronary arterial disease patients. Clin Chim Acta. Oct 21 2013;425:9–17. doi:10.1016/j.cca.2013.06.021
    OpenUrlCrossRef
  27. 27.↵
    Barrett T, Wilhite SE, Ledoux P, et al. NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res. Jan 2013;41(Database issue):D991-5. doi:10.1093/nar/gks1193
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    Ritchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. Apr 20 2015;43(7):e47. doi:10.1093/nar/gkv007
    OpenUrlCrossRefPubMed
  29. 29.↵
    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. Aug 01 2014;30(15):2114–20. doi:10.1093/bioinformatics/btu170
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. Apr 2017;14(4):417–419. doi:10.1038/nmeth.4197
    OpenUrlCrossRefPubMed
  31. 31.↵
    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. doi:10.1186/s13059-014-0550-8
    OpenUrlCrossRefPubMed
  32. 32.↵
    Leek JT, Storey JD. Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet. Sep 2007;3(9):1724–35. doi:10.1371/journal.pgen.0030161
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    Nguyen AB, Cavallari LH, Rossi JS, Stouffer GA, Lee CR. Evaluation of race and ethnicity disparities in outcome studies of. Front Cardiovasc Med. 2022;9:991646. doi:10.3389/fcvm.2022.991646
    OpenUrlCrossRef
  34. 34.↵
    Yang G, Alarcon C, Friedman P, et al. The Role of Global and Local Ancestry on Clopidogrel Response in African Americans. Pac Symp Biocomput. 2023;28:221–232.
    OpenUrl
  35. 35.↵
    Verma SS, Bergmeijer TO, Gong L, et al. Genomewide Association Study of Platelet Reactivity and Cardiovascular Response in Patients Treated With Clopidogrel: A Study by the International Clopidogrel Pharmacogenomics Consortium. Clin Pharmacol Ther. 11 2020;108(5):1067–1077. doi:10.1002/cpt.1911
    OpenUrlCrossRef
  36. 36.↵
    de Visser MC, van Minkelen R, van Marion V, et al. Genome-wide linkage scan in affected sibling pairs identifies novel susceptibility region for venous thromboembolism: Genetics In Familial Thrombosis study. J Thromb Haemost. Aug 2013;11(8):1474–84. doi:10.1111/jth.12313
    OpenUrlCrossRef
  37. 37.↵
    Li Y, Wang DW, Chen Y, et al. Genome-Wide Association and Functional Studies Identify. Arterioscler Thromb Vasc Biol. 04 2018;38(4):964–975. doi:10.1161/ATVBAHA.117.310594
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    Kuo MW, Wang CH, Wu HC, Chang SJ, Chuang YJ. Soluble THSD7A is an N-glycoprotein that promotes endothelial cell migration and tube formation in angiogenesis. PLoS One. 2011;6(12):e29000. doi:10.1371/journal.pone.0029000
    OpenUrlCrossRefPubMed
  39. 39.↵
    Wang CH, Su PT, Du XY, et al. Thrombospondin type I domain containing 7A (THSD7A) mediates endothelial cell migration and tube formation. J Cell Physiol. Mar 2010;222(3):685–94. doi:10.1002/jcp.21990
    OpenUrlCrossRefPubMed
  40. 40.↵
    Boehlen F, Clemetson KJ. Platelet chemokines and their receptors: what is their relevance to platelet storage and transfusion practice? Transfus Med. Dec 2001;11(6):403–17. doi:10.1046/j.1365-3148.2001.00340.x
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.
    Schober A, Manka D, von Hundelshausen P, et al. Deposition of platelet RANTES triggering monocyte recruitment requires P-selectin and is involved in neointima formation after arterial injury. Circulation. Sep 17 2002;106(12):1523–9. doi:10.1161/01.cir.0000028590.02477.6f
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    Reiner MF, Akhmedov A, Stivala S, et al. Ticagrelor, but not clopidogrel, reduces arterial thrombosis via endothelial tissue factor suppression. Cardiovasc Res. Jan 2017;113(1):61–69. doi:10.1093/cvr/cvw233
    OpenUrlCrossRefPubMed
  43. 43.↵
    Sun Y, Lu Q, Tao X, Cheng B, Yang G. *2 Polymorphism Related to Clopidogrel Resistance in Patients With Coronary Heart Disease, Especially in the Asian Population: A Systematic Review and Meta-Analysis. Front Genet. 2020;11:576046. doi:10.3389/fgene.2020.576046
    OpenUrlCrossRef
  44. 44.↵
    Jeong YH, Tantry US, Kim IS, et al. Effect of CYP2C19*2 and *3 loss-of-function alleles on platelet reactivity and adverse clinical events in East Asian acute myocardial infarction survivors treated with clopidogrel and aspirin. Circ Cardiovasc Interv. Dec 01 2011;4(6):585–94. doi:10.1161/CIRCINTERVENTIONS.111.962555
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    Steevels TA, Westerlaken GH, Tijssen MR, et al. Co-expression of the collagen receptors leukocyte-associated immunoglobulin-like receptor-1 and glycoprotein VI on a subset of megakaryoblasts. Haematologica. Dec 2010;95(12):2005–12. doi:10.3324/haematol.2010.026120
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    Verbrugge A, de Ruiter T, Geest C, Coffer PJ, Meyaard L. Differential expression of leukocyte-associated Ig-like receptor-1 during neutrophil differentiation and activation. J Leukoc Biol. Apr 2006;79(4):828–36. doi:10.1189/jlb.0705370
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    Perrella G, Nagy M, Watson SP, Heemskerk JWM. Platelet GPVI (Glycoprotein VI) and Thrombotic Complications in the Venous System. Arterioscler Thromb Vasc Biol. Nov 2021;41(11):2681–2692. doi:10.1161/ATVBAHA.121.316108
    OpenUrlCrossRef
  48. 48.↵
    Schwartz JC, Sanderson ND, Bickhart DM, Smith TPL, Hammond JA. The Structure, Evolution, and Gene Expression Within the Caprine Leukocyte Receptor Complex. Front Immunol. 2019;10:2302. doi:10.3389/fimmu.2019.02302
    OpenUrlCrossRef
  49. 49.↵
    Brondijk TH, de Ruiter T, Ballering J, et al. Crystal structure and collagen-binding site of immune inhibitory receptor LAIR-1: unexpected implications for collagen binding by platelet receptor GPVI. Blood. Feb 18 2010;115(7):1364–73. doi:10.1182/blood-2009-10-246322
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    Smith CW, Thomas SG, Raslan Z, et al. Mice Lacking the Inhibitory Collagen Receptor LAIR-1 Exhibit a Mild Thrombocytosis and Hyperactive Platelets. Arterioscler Thromb Vasc Biol. May 2017;37(5):823–835. doi:10.1161/ATVBAHA.117.309253
    OpenUrlAbstract/FREE Full Text
  51. 51.↵
    Rayes J, Watson SP, Nieswandt B. Functional significance of the platelet immune receptors GPVI and CLEC-2. J Clin Invest. Jan 02 2019;129(1):12–23. doi:10.1172/JCI122955
    OpenUrlCrossRefPubMed
Back to top
PreviousNext
Posted December 07, 2023.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Investigation of genomic and transcriptomic risk factors in clopidogrel response in African Americans
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Investigation of genomic and transcriptomic risk factors in clopidogrel response in African Americans
Guang Yang, Cristina Alarcon, Catherine Chanfreau, Norman H. Lee, Paula Friedman, Edith Nutescu, Matthew Tuck, Travis O’Brien, Li Gong, Teri E. Klein, Kyong-Mi Chang, Philip S. Tsao, David O. Meltzer, Million Veteran Program, Sony Tuteja, Minoli A Perera
medRxiv 2023.12.05.23299140; doi: https://doi.org/10.1101/2023.12.05.23299140
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Investigation of genomic and transcriptomic risk factors in clopidogrel response in African Americans
Guang Yang, Cristina Alarcon, Catherine Chanfreau, Norman H. Lee, Paula Friedman, Edith Nutescu, Matthew Tuck, Travis O’Brien, Li Gong, Teri E. Klein, Kyong-Mi Chang, Philip S. Tsao, David O. Meltzer, Million Veteran Program, Sony Tuteja, Minoli A Perera
medRxiv 2023.12.05.23299140; doi: https://doi.org/10.1101/2023.12.05.23299140

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Pharmacology and Therapeutics
Subject Areas
All Articles
  • Addiction Medicine (349)
  • Allergy and Immunology (668)
  • Allergy and Immunology (668)
  • Anesthesia (181)
  • Cardiovascular Medicine (2648)
  • Dentistry and Oral Medicine (316)
  • Dermatology (223)
  • Emergency Medicine (399)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (942)
  • Epidemiology (12228)
  • Forensic Medicine (10)
  • Gastroenterology (759)
  • Genetic and Genomic Medicine (4103)
  • Geriatric Medicine (387)
  • Health Economics (680)
  • Health Informatics (2657)
  • Health Policy (1005)
  • Health Systems and Quality Improvement (985)
  • Hematology (363)
  • HIV/AIDS (851)
  • Infectious Diseases (except HIV/AIDS) (13695)
  • Intensive Care and Critical Care Medicine (797)
  • Medical Education (399)
  • Medical Ethics (109)
  • Nephrology (436)
  • Neurology (3882)
  • Nursing (209)
  • Nutrition (577)
  • Obstetrics and Gynecology (739)
  • Occupational and Environmental Health (695)
  • Oncology (2030)
  • Ophthalmology (585)
  • Orthopedics (240)
  • Otolaryngology (306)
  • Pain Medicine (250)
  • Palliative Medicine (75)
  • Pathology (473)
  • Pediatrics (1115)
  • Pharmacology and Therapeutics (466)
  • Primary Care Research (452)
  • Psychiatry and Clinical Psychology (3432)
  • Public and Global Health (6527)
  • Radiology and Imaging (1403)
  • Rehabilitation Medicine and Physical Therapy (814)
  • Respiratory Medicine (871)
  • Rheumatology (409)
  • Sexual and Reproductive Health (410)
  • Sports Medicine (342)
  • Surgery (448)
  • Toxicology (53)
  • Transplantation (185)
  • Urology (165)