Implementing a Pharmacogenomic-driven Algorithm to Guide Antiplatelet Therapy among Caribbean Hispanics: A non-randomized prospective cohort study

Background

After percutaneous coronary intervention (PCI), clopidogrel resistant patients are at an increased risk of major adverse cardiovascular and cerebrovascular events (MACCEs). We aimed to assess whether genotype-guided selection of oral antiplatelet drugs using a clinical decision support (CDS) algorithm reduces the occurrence of these ischemic events and improves outcomes among Caribbean Hispanic patients from Puerto Rico, who are underrepresented in clinical pharmacogenomic (PGx)-guided implementation studies.

Methods

Individual platelet function testing (PRU) measures, CYP2C19*2 and PON1 rs662 genotypes, clinical and demographic data from 8 medical facilities were included. Patients were separated into standard of care (SoC) and genotype-guided groups (150 each). Risk scores were calculated based on a previously developed CDS risk prediction algorithm designed to make actionable treatment recommendations for each patient. Alternative therapy with ticagrelor was recommended for patients with a high risk score ≥2. Statistical associations between patient time free of MACCEs and predictor variables (i.e., treatment groups, risk scores) were tested in this population using Kaplan-Meier survival analyses and Cox proportional-hazards regression models.

Results

Median age of participants is 67 years; BMI: 27.8; 48% women; 14% smokers; 59% with type-2 diabetes mellitus (T2DM). Among patients with high-risk scores who were free from MACCE events 6 months after coronary stenting, genotype-driven guidance of antiplatelet therapy showed superiority over SoC in terms of reducing the incidence rate of atherothrombotic events.

Conclusions

The clinical utility of our PGx-driven CDS algorithm to reduce the incidence rate of MACCEs among post-PCI Caribbean Hispanic patients on clopidogrel was externally demonstrated.

Clinical Trial Registration Unique Identifier

NCT03419325.

Discovery of Ancestry-specific Variants Associated with Clopidogrel Response among Caribbean Hispanics

Background

High on-treatment platelet reactivity (HTPR) with clopidogrel is predictive of ischemic events in adults with coronary artery disease. Despite strong data suggesting HTPR varies with ethnicity, including clinical and genetic variables, no genome-wide association study (GWAS) of clopidogrel response has been performed among Caribbean Hispanics. This study aimed to identify genetic predictors of HTPR in a cohort of Caribbean Hispanic cardiovascular patients from Puerto Rico.

Methods

Local Ancestry inference (LAI) and traditional GWASs were performed on a cohort of 511 clopidogrel-treated patients, stratified based on their P2Y12 reaction units (PRU) into responders and non-responders (HTPR).

Results

The LAI GWAS identified variants within the CYP2C19 region associated with HTPR, predominantly driven by individuals of European ancestry and absent in those with native ancestry. Incorporating local ancestry adjustment notably enhanced our ability to detect associations. While no loci reached traditional GWAS significance, three variants showed suggestive significance at chromosomes 3, 14 and 22 (OSBPL10 rs1376606, DERL3 rs5030613, and RGS6 rs9323567). In addition, a variant in the UNC5C gene on chromosome 4 was associated with an increased risk of HTPR. These findings were not identified in other cohorts, highlighting the unique genetic landscape of Caribbean Hispanics.

Conclusion

This is the first GWAS of clopidogrel response in Hispanics, confirming the relevance of the CYP2C19 cluster, particularly among those with European ancestry, and also identifying novel markers in a diverse patient population. Further studies are warranted to replicate our findings in other diverse cohorts and meta-analyses.

Pharmacogenomic polygenic risk score for clopidogrel responsiveness among Caribbean Hispanics: A candidate gene approach

This multicenter clinical study was aimed at conducting a targeted pharmacogenomic association analysis of residual on‐clopidogrel platelet reactivity in 474 Caribbean Hispanic patients. Platelet reactivity was measured using the VerifyNow P2Y12 assay and clopidogrel resistance was defined as P2Y12 reaction units (PRUs) greater than or equal to 208. Genotyping was performed using the whole‐genome Infinium MEGA BeadChip array. An ancestry‐adjusted, weighted polygenic risk score (wPGxRS) was developed to account for the effect of multiple variants on PRU and compared between clopidogrel responders and nonresponders. The mean PRU across the study cohort was 173.8 ± 68.5 and 33.5% of patients were defined as clopidogrel resistant. Multivariate linear regression showed that 19% of PRU variability was attributed to nine independent predictors, with CYP2C19*2 (rs4244285) accounting for ~ 7% of observed PRU variation (p < 0.001). PON1 rs662, ABCB1/MDR1 rs2032582, PEAR1 rs12041331 carrier status, and the interaction between African ancestry and rs12041331 carriers also predicted PRU among the participants (p ≤ 0.05). A clear gene‐dose effect was detected between PRU and CYP2C19*2 genotype, consistent with previous studies in European patient populations, as well as rs12777823. Importantly, a significant positive correlation was detected between our novel wPGxRS (4 variants) and PRU among the Hispanic patient population (r_p = 0.35, p < 0.001). Moreover, the wPGxRS discriminated between nonresponders and responders (p = 0.003), indicating that this multigene‐based score is a useful predictor of clopidogrel resistance among Caribbean Hispanics. Taken together, these results help close the gap of knowledge on clopidogrel pharmacogenomics and its potential clinical implementation in this under‐represented population.

Genetic Factors Influencing Warfarin Dose in Black-African Patients: A Systematic Review and Meta-Analysis

Warfarin is the most commonly used oral anticoagulant in sub-Saharan Africa. Dosing is challenging due to a narrow therapeutic index and high interindividual variability in dose requirements. To evaluate the genetic factors affecting warfarin dosing in black-Africans, we performed a meta-analysis of 48 studies (2,336 patients). Significant predictors for CYP2C9 and stable dose included rs1799853 (CYP2C9*2), rs1057910 (CYP2C9*3), rs28371686 (CYP2C9*5), rs9332131 (CYP2C9*6), and rs28371685 (CYP2C9*11) reducing dose by 6.8, 12.5, 13.4, 8.1, and 5.3 mg/week, respectively. VKORC1 variants rs9923231 (-1639G>A), rs9934438 (1173C>T), rs2359612 (2255C>T), rs8050894 (1542G>C), and rs2884737 (497T>G) decreased dose by 18.1, 21.6, 17.3, 11.7, and 19.6 mg/week, respectively, whereas rs7294 (3730G>A) increased dose by 6.9 mg/week. Finally, rs12777823 (CYP2C gene cluster) was associated with a dose reduction of 12.7 mg/week. Few studies were conducted in Africa, and patient numbers were small, highlighting the need for further work in black-Africans to evaluate genetic factors determining warfarin response.

Pharmacogenetic association study on clopidogrel response in Puerto Rican Hispanics with cardiovascular disease: a novel characterization of a Caribbean population

Introduction

High on-treatment platelet reactivity (HTPR) to clopidogrel imparts an increased risk for ischemic events in adults with coronary artery disease. Platelet reactivity varies with ethnicity and is influenced by both clinical and genetic variables; however, no clopidogrel pharmacogenetic studies with Puerto Rican patients have been reported. Therefore, we sought to identify clinical and genetic determinants of on-treatment platelet reactivity in a cohort of Puerto Rican patients with cardiovascular disease.

Methods

We performed a retrospective study of 111 patients on 75 mg/day maintenance dose of clopidogrel. Patients were allocated into 2 groups: Group I, without HTPR; and Group II, with HTPR. Platelet function was measured ex vivo using the VerifyNow® P2Y12 assay and HTPR was defined as P2Y12 reaction units (PRU) ≥230. Genotyping testing was performed using Taqman^® Genotyping Assays.

Results

The mean PRU across the cohort was 203±61 PRU (range 8–324), and 42 (38%) patients had HTPR. Multiple logistic regression showed that 27% of the total variation in PRU was explained by a history of diabetes mellitus, hematocrit, CYP2C19*2, and PON1 p.Q192R. Body mass index (odds ratio [OR]=1.15; 95% CI: 1.03–1.27), diabetes mellitus (OR=3.46; 95% CI: 1.05–11.43), hematocrit (OR=0.75; 95% CI: 0.65–0.87), and CYP2C19*2 (OR=4.44; 95% CI: 1.21–16.20) were the only independent predictors of HTPR.

Conclusion

Moreover, we propose a predictive model to determine PRU values as measured by VerifyNow P2Y12 assay for the Puerto Rican Hispanic population. This model has the potential to identify Hispanic patients at higher risk for adverse events on clopidogrel.

Clinical Relevant Polymorphisms Affecting Clopidogrel Pharmacokinetics and Pharmacodynamics: Insights from the Puerto Rico Newborn Screening Program

Background: Variations in several clopidogrel-pharmacogenes have been linked to clopidogrel response variability and clinical outcomes. We aimed to determine the frequency distribution of major polymorphisms on CYP2C19, PON1, ABCB1 and P2RY12 pharmacogenes in Puerto Ricans. Methods: This was a cross-sectional, population-based study of 200 unrelated "Guthrie" cards specimens from newborns registered in the Puerto Rican newborn screening program (PRNSP) between 2004 and 2014. Taqman® SNP assay techniques were used for genotyping. Results: Minor allele frequencies (MAF) were 46% for PON1 (rs662), 41% for ABCB1 (rs1045642), 14% for CYP2C19*17, 13% for CYP2C19*2, 12% for P2RY12-H2 and 0.3% for CYP2C19*4. No carriers of the CYP2C19*3 variants were detected. All alleles and genotype proportions were found to be in Hardy⁻Weinberg equilibrium (HWE). Overall, there were no significant differences between MAFs of these variants in Puerto Ricans and the general population (n = 453) of the 1000 Genome project, except when comparisons to each individual parental group were performed (i.e., Africans, Europeans and East-Asians; p < 0.05). As expected, the prevalence of these markers in Puerto Ricans most resembled those in the 181 subjects from reference populations of the Americas. Conclusions: These prevalence data provide a necessary groundwork for future clinical studies of clopidogrel pharmacogenetics in Caribbean Hispanics.

Role of treatment-modifying MTHFR677C>T and 1298A>C polymorphisms in metformin-treated Puerto Rican patients with type-2 diabetes mellitus and peripheral neuropathy

BACKGROUND

The study was conducted to investigate potential association between MTHFR genotypes and diabetic peripheral neuropathy (DPN) in Puerto Ricans with type-2 diabetes mellitus (T2DM) treated with metformin. The prevalence of major MTHFR polymorphisms in this cohort was also ascertained.

METHODS

DNAs from 89 metformin-treated patients with T2DM and DPN were genotyped using the PCR-based RFLP assay for MTHFR677C>T and 1298A>C polymorphisms. Frequency distributions of these variants in the study cohort were compared to those reported for three reference populations (HapMap project) and controls (400 newborn specimens). Chi-square (or Fischer's exact) tests and odds ratios (OR) were used to assess association with DPN susceptibility risk (patients vs. controls) and biochemical markers (wild types vs. carriers).

RESULTS

Sixty-seven percent (67%) of participants carry at least one of these MTHFR polymorphisms. No deviations from Hardy-Weinberg equilibrium were detected. The genotype and allele frequencies showed statistically significant differences between participants and controls (p<0.0001 and p=0.03, respectively). Results suggest that 1298A>C but not 677C>T is associated with DPN susceptibility in this cohort (p=0.018). Different patterns of allelic dissimilarities are observed when comparing our cohort vs. the three parental ancestries. After sorting individuals by their carrier status, no significant associations were observed between these genetic variants (independently or combined) and any of the biochemical markers (HbA1c, folate, vitamin B12, homocysteine).

CONCLUSIONS

Prevalence of major MTHFR variants in Puerto Rican patients with T2DM is first time ever reported. The study provides further evidence on the use of this genetic marker as an independent risk factor for DPN.

Pharmacogenetics of healthy volunteers in Puerto Rico

Puerto Ricans are a unique Hispanic population with European, Native American (Taino), and higher West African ancestral contributions than other non-Caribbean Hispanics. In admixed populations, such as Puerto Ricans, genetic variants can be found at different frequencies when compared to parental populations and uniquely combined and distributed. Therefore, in this review, we aimed to collect data from studies conducted in healthy Puerto Ricans and to report the frequencies of genetic polymorphisms with major relevance in drug response. Filtering for healthy volunteers or individuals, we performed a search of pharmacogenetic studies in academic literature databases without limiting the period of the results. The search was limited to Puerto Ricans living in the island, excluding those studies performed in mainland (United States). We found that the genetic markers impacting pharmacological therapy in the areas of cardiovascular, oncology, and neurology are the most frequently investigated. Coincidently, the top causes of mortality in the island are cardiovascular diseases, cancer, diabetes, Alzheimer's disease, and stroke. In addition, polymorphisms in genes that encode for members of the CYP450 family (CYP2C9, CYP2C19, and CYP2D6) are also available due to their relevance in the metabolism of drugs. The complex genetic background of Puerto Ricans is responsible for the divergence in the reported allele frequencies when compared to parental populations (Africans, East Asians, and Europeans). The importance of reporting the findings of pharmacogenetic studies conducted in Puerto Ricans is to identify genetic variants with potential utility among this genetically complex population and eventually move forward the adoption of personalized medicine in the island.

A Novel Admixture-Based Pharmacogenetic Approach to Refine Warfarin Dosing in Caribbean Hispanics

Aim

This study is aimed at developing a novel admixture-adjusted pharmacogenomic approach to individually refine warfarin dosing in Caribbean Hispanic patients.

Patients & Methods

A multiple linear regression analysis of effective warfarin doses versus relevant genotypes, admixture, clinical and demographic factors was performed in 255 patients and further validated externally in another cohort of 55 individuals.

Results

The admixture-adjusted, genotype-guided warfarin dosing refinement algorithm developed in Caribbean Hispanics showed better predictability (R² = 0.70, MAE = 0.72mg/day) than a clinical algorithm that excluded genotypes and admixture (R² = 0.60, MAE = 0.99mg/day), and outperformed two prior pharmacogenetic algorithms in predicting effective dose in this population. For patients at the highest risk of adverse events, 45.5% of the dose predictions using the developed pharmacogenetic model resulted in ideal dose as compared with only 29% when using the clinical non-genetic algorithm (p<0.001). The admixture-driven pharmacogenetic algorithm predicted 58% of warfarin dose variance when externally validated in 55 individuals from an independent validation cohort (MAE = 0.89 mg/day, 24% mean bias).

Conclusions

Results supported our rationale to incorporate individual’s genotypes and unique admixture metrics into pharmacogenetic refinement models in order to increase predictability when expanding them to admixed populations like Caribbean Hispanics.

Trial Registration

ClinicalTrials.gov NCT01318057

Frequencies of Functional Polymorphisms in Three Pharmacokinetic Genes of Clinical Interest within the Admixed Puerto Rican Population

OBJECTIVE

This cross-sectional study was aimed at determining the allele frequencies for the CYP2C19*2, CYP2C19*3, CYP2D6*10 and PON1 (rs662) polymorphisms in the Puerto Rican population. The CYP2C19, CYP2D6 and PON1 genes are known to be associated with functional changes in drug metabolism and activation. Individuals carrying the aforementioned polymorphisms are at a higher risk of suffering from drug-induced adverse events and/ or unresponsiveness from a variety of drugs that includes antidepressants, atypical antipsychotics and antiplatelet compounds. Information on the frequency of these polymorphisms is more commonly found on homogeneous populations, but is scarce in highly heterogeneous populations like Hispanics, as in the case of Puerto Ricans.

METHOD

Genotyping was carried out in 100 genomic DNA samples from dried blood spots supplied by the Puerto Rican Newborn Screening program using Taqman® Genotyping Assays.

RESULTS

The Minor Allele Frequencies (MAF) obtained were 9% for CYP2C19*2 and CYP2D6*10, 50% for PON1 (rs662), while the CYP2C19*3 variant was not detected in our study. Furthermore, Hardy Weinberg equilibrium analysis was assessed as well as a comparison between Puerto Rico and other reference populations using a Z-test for proportions.

CONCLUSION

The observed allele and genotype frequencies on these relevant pharmacogenes in Puerto Ricans were more closely related to those early reported in two other reference populations of Americans (Mexicans and Colombians).

Genetic testing for oculocutaneous albinism type 1 and 2 and Hermansky-Pudlak syndrome type 1 and 3 mutations in Puerto Rico

Hermansky-Pudlak syndrome (HPS) (MIM #203300) is a heterogeneous group of autosomal recessive disorders characterized by oculocutaneous albinism (OCA), bleeding tendency, and lysosomal dysfunction. HPS is very common in Puerto Rico (PR), particularly in the northwest part of the island, with a frequency of approximately 1:1,800. Two HPS genes and mutations have been identified in PR, a 16-base pair (bp) duplication in HPS1 and a 3,904-bp deletion in HPS3. In Puerto Ricans with more typical OCA, the most common mutation of the tyrosinase (TYR) (human tyrosinase (OCA1) gene) gene was G47D. We describe screening 229 Puerto Rican OCA patients for these mutations, and for mutations in the OCA2 gene. We found the HPS1 mutation in 42.8% of cases, the HPS3 deletion in 17%, the TYR G47D mutation in 3.0%, and a 2.4-kb deletion of the OCA2 gene in 1.3%. Among Puerto Rican newborns, the frequency of the HPS1 mutation is highest in northwest PR (1:21; 4.8%) and lower in central PR (1:64; 1.6%). The HPS3 gene deletion is most frequent in central PR (1:32; 3.1%). Our findings provide insights into the genetics of albinism and HPS in PR, and provide the basis for genetic screening for these disorders in this minority population.

The cDNA-derived amino acid sequence of hemoglobin II from Lucina pectinata

Hemoglobin II from the clam Lucina pectinata is an oxygen-reactive protein with a unique structural organization in the heme pocket involving residues Gln65 (E7), Tyr30 (B10), Phe44 (CD1), and Phe69 (E11). We employed the reverse transcriptase-polymerase chain reaction (RT-PCR) and methods to synthesize various cDNA(HbII). An initial 300-bp cDNA clone was amplified from total RNA by RT-PCR using degenerate oligonucleotides. Gene-specific primers derived from the HbII-partial cDNA sequence were used to obtain the 5' and 3' ends of the cDNA by RACE. The length of the HbII cDNA, estimated from overlapping clones, was approximately 2114 bases. Northern blot analysis revealed that the mRNA size of HbII agrees with the estimated size using cDNA data. The coding region of the full-length HbII cDNA codes for 151 amino acids. The calculated molecular weight of HbII, including the heme group and acetylated N-terminal residue, is 17,654.07 Da.

Characterization of hemoglobin Hotel Dieu in a Puerto Rican adolescent

PURPOSE

Hemoglobin Hotel Dieu (HbHD) is a high-oxygen affinity variant of HbA never before reported in a Hispanic patient. This Hb variant was first reported in 1981 by Blouquit et al. in a white person with erythrocytosis with a substitution in the beta 99 aspartic acid residue by glycine.

METHODS

A 13-year-old Puerto Rican boy had pain in his chest, headaches, easy fatigability, and high Hb (as high as 19.1 g/dl). Protein analysis was performed by cellulose acetate, citrate agar, and isoelectric focusing electrophoresis and high-pressure liquid chromatography (HPLC), polymerase chain reaction (PCR) amplification, and DNA sequencing of the second exon of the beta gene in samples obtained from the mother, father, and the patient, and DNA fingerprinting to determine paternity.

RESULTS

The variant found in the patient migrated on cellulose acetate electrophoresis to a cathodic position relative to HbF, and a band cathodal to HbA and close to HbF on isoelectric focusing electrophoresis. The patient showed an abnormal well-resolved peak on HPLC with a retention time slightly shorter than that for HbS. DNA analysis by direct sequencing of the PCR product demonstrated heterozygosity for codon 99 (GAT-->GGT) in the patient but not in either parent. DNA fingerprinting by multiplex PCR amplification of three simple tandem repeat loci showed that the patient shared alleles in all three loci with both parents, ruling out nonpaternity.

CONCLUSIONS

The protein and DNA analysis indicate that the erythrocytosis is caused by the presence of HbHD in this Hispanic adolescent.

Longitudinal studies on maternal HIV-1 variants by biological phenotyping, sequence analysis and viral load

In this study, the HIV-1 variant viruses from ten pregnant women and their infants were isolated and characterized longitudinally in order to determine the role that viral envelope (gp120-V3 loop) gene variation and viral tropism play in vertical transmission. Biological phenotyping of each HIV variant was accomplished by growth in MT-2, and macrophages from healthy and non-HIV-infected donors. Genetic characterization of the variants was accomplished by DNA sequence analysis. All the women enrolled in this study received ZDV therapy. Virus was cultured from eight out of ten env V3-PCR positive mothers. HIV-1 isolates were all non-syncitium inducing variants. None of the mothers were found to transmit HIV, as determined by DNA PCR and quantitative co-cultures on their infants which were seronegative for HIV-1 through one year after birth. Viral cultures from infant blood samples were negative and infants were all healthy. However, nested env V3-PCR detected proviral DNA in five out of ten infants. In contrast, conventional gag-PCR was negative in the same five infants. Sequences of the five maternal-infant pairs were different, suggesting unique infant HIV-1 variants. The three highest maternal viral load values corresponded to infants that were env V3-PCR positive. These results suggest that HIV-1 particles are transmitted from ZDV-treated mothers to infants. Infant follow up is recommended to determine if HIV-1 has been inhibited by the immune system of the infants.