Geographic variation of mutagenic exposures in kidney cancer genomes

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

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.

Identification of adverse drug reactions that may be related to pharmacogenetics in a public hospital in the South of Brazil

BACKGROUND

Adverse drug reactions (ADRs) are of great concern in clinical practice. Pharmacogenetics can identify individuals and groups at increased risk of developing ADRs, enabling treatment adjustments to improve outcomes. The study aimed to determine the prevalence of ADRs related to drugs with pharmacogenetic evidence level 1A in a public hospital in Southern Brazil.

RESEARCH DESIGN AND METHODS

ADR information was collected from the pharmaceutical registries from 2017 to 2019. Drugs that have pharmacogenetic evidence level 1A were selected. Public genomic databases were used to estimate the genotypes/phenotypes frequency.

RESULTS

During the period, 585 ADRs were spontaneously notified. Most were moderate (76.3%), whereas severe reactions accounted for 33.8%. Additionally, 109 ADRs caused by 41 drugs presented pharmacogenetic evidence level 1A, representing 18.6% of all notified reactions. Depending on the drug-gene pair, up to 35% of individuals from Southern Brazil could be at risk of developing ADRs.

CONCLUSIONS

Relevant amount of ADRs were related to drugs with pharmacogenetic recommendations on drug labels and/or guidelines. Genetic information could guide and improve clinical outcomes, decreasing ADR incidence and reducing treatment costs.

Long-read HiFi sequencing of NUDT15: Phased full-gene haplotyping and pharmacogenomic allele discovery

To determine the phase of NUDT15 sequence variants for more comprehensive star (*) allele diplotyping, we developed a novel long-read single-molecule real-time HiFi amplicon sequencing method. A 10.5 kb NUDT15 amplicon assay was validated using reference material positive controls and additional samples for specimen type and blinded accuracy assessment. Triplicate NUDT15 HiFi sequencing of two reference material samples had nonreference genotype concordances of >99.9%, indicating that the assay is robust. Notably, short-read genome sequencing of a subset of samples was unable to determine the phase of star (*) allele-defining NUDT15 variants, resulting in ambiguous diplotype results. In contrast, long-read HiFi sequencing phased all variants across the NUDT15 amplicons, including a *2/*9 diplotype that previously was characterized as *1/*2 in the 1000 Genomes Project v3 data set. Assay throughput was also tested using 8.5 kb amplicons from 100 Ashkenazi Jewish individuals, which identified a novel NUDT15 *1 suballele (c.-121G>A) and a rare likely deleterious coding variant (p.Pro129Arg). Both novel alleles were Sanger confirmed and assigned as *1.007 and *20, respectively, by the PharmVar Consortium. Taken together, NUDT15 HiFi amplicon sequencing is an innovative method for phased full-gene characterization and novel allele discovery, which could improve NUDT15 pharmacogenomic testing and subsequent phenotype prediction.

Sales of "COVID kit" drugs and adverse drug reactions reported by the Brazilian Health Regulatory Agency

Off-label use of azithromycin, hydroxychloroquine, and ivermectin (the "COVID kit") has been suggested for COVID-19 treatment in Brazil without clinical or scientific evidence of efficacy. These drugs have known adverse drug reactions (ADR). This study aimed to analyze if the sales of drugs in the "COVID kit" are correlated to the reported number of ADR after the COVID-19 pandemic began. Data was obtained from the Brazilian Health Regulatory Agency (Anvisa) website on reported sales and ADRs for azithromycin, hydroxychloroquine, and ivermectin for all Brazilian states. The period from March 2019 to February 2020 (before the pandemic) was compared to that from March 2020 to February 2021 (during the pandemic). Trend adjustment was performed for time series data and cross-correlation analysis to investigate correlation between sales and ADR within the same month (lag 0) and in the following months (lag 1 and lag 2). Spearman's correlation coefficient was used to assess the magnitude of the correlations. After the pandemic onset, sales of all investigated drugs increased significantly (69.75% for azithromycin, 10,856,481.39% for hydroxychloroquine, and 12,291,129.32% for ivermectin). ADR levels of all medications but azithromycin were zero before the pandemic, but increased after its onset. Cross-correlation analysis was significant in lag 1 for all drugs nationwide. Spearman's correlation was moderate for azithromycin and hydroxychloroquine but absent for ivermectin. Data must be interpreted cautiously since no active search for ADR was performed. Our results show that the increased and indiscriminate use of "COVID kit" during the pandemic correlates to an increased occurrence of ADRs.

Development and Analytical Validation of a 29 Gene Clinical Pharmacogenetic Genotyping Panel: Multi-Ethnic Allele and Copy Number Variant Detection

To develop a novel pharmacogenetic genotyping panel, a multidisciplinary team evaluated available evidence and selected 29 genes implicated in interindividual drug response variability, including 130 sequence variants and additional copy number variants (CNVs). Of the 29 genes, 11 had guidelines published by the Clinical Pharmacogenetics Implementation Consortium. Targeted genotyping and CNV interrogation were accomplished by multiplex single‐base extension using the MassARRAY platform (Agena Biosciences) and multiplex ligation‐dependent probe amplification (MRC Holland), respectively. Analytical validation of the panel was accomplished by a strategic combination of > 500 independent tests performed on 170 unique reference material DNA samples, which included sequence variant and CNV accuracy, reproducibility, and specimen (blood, saliva, and buccal swab) controls. Among the accuracy controls were 32 samples from the 1000 Genomes Project that were selected based on their enrichment of sequence variants included in the pharmacogenetic panel (VarCover.org). Coupled with publicly available samples from the Genetic Testing Reference Materials Coordination Program (GeT‐RM), accuracy validation material was available for the majority (77%) of interrogated sequence variants (100% with average allele frequencies > 0.1%), as well as additional structural alleles with unique copy number signatures (e.g., CYP2D6*5, *13, *36, *68; CYP2B6*29; and CYP2C19*36). Accuracy and reproducibility for both genotyping and copy number were > 99.9%, indicating that the optimized panel platforms were precise and robust. Importantly, multi‐ethnic allele frequencies of the interrogated variants indicate that the vast majority of the general population carries at least one of these clinically relevant pharmacogenetic variants, supporting the implementation of this panel for pharmacogenetic research and/or clinical implementation programs.

PharmVar GeneFocus: CYP2C19

The Pharmacogene Variation Consortium (PharmVar) catalogues star (*) allele nomenclature for the polymorphic human CYP2C19 gene. CYP2C19 genetic variation impacts the metabolism of many drugs and has been associated with both efficacy and safety issues for several commonly prescribed medications. This GeneFocus provides a comprehensive overview and summary of CYP2C19 and describes how haplotype information catalogued by PharmVar is utilized by the Pharmacogenomics Knowledgebase and the Clinical Pharmacogenetics Implementation Consortium (CPIC).

Structural variation at the CYP2C locus: Characterization of deletion and duplication alleles

The human CYP2C locus harbors the polymorphic CYP2C18, CYP2C19, CYP2C9, and CYP2C8 genes, and of these, CYP2C19 and CYP2C9 are directly involved in the metabolism of ~15% of all medications. All variant CYP2C19 and CYP2C9 star (*) allele haplotypes currently cataloged by the Pharmacogene Variation (PharmVar) Consortium are defined by sequence variants. To determine if structural variation also occurs at the CYP2C locus, the 10q23.33 region was interrogated across deidentified clinical chromosomal microarray (CMA) data from 20,642 patients tested at two academic medical centers. Fourteen copy number variants that affected the coding region of CYP2C genes were detected in the clinical CMA cohorts, which ranged in size from 39.2 to 1,043.3 kb. Selected deletions and duplications were confirmed by MLPA or ddPCR. Analysis of the clinical CMA and an additional 78,839 cases from the Database of Genomic Variants (DGV) and ClinGen (total n = 99,481) indicated that the carrier frequency of a CYP2C structural variant is ~1 in 1,000, with ~1 in 2,000 being a CYP2C19 full gene or partial-gene deletion carrier, designated by PharmVar as CYP2C19*36 and *37, respectively. Although these structural variants are rare in the general population, their detection will likely improve metabolizer phenotype prediction when interrogated for research and/or clinical testing.

Identification of environmental and genetic factors that influence warfarin time in therapeutic range

Warfarin is an oral anticoagulant prescribed to prevent and treat thromboembolic disorders. It has a narrow therapeutic window and must have its effect controlled. Prothrombin test, expressed in INR value, is used for dose management. Time in therapeutic range (TTR) is an important outcome of quality control of anticoagulation therapy and is influenced by several factors. The aim of this study was to identify genetic, demographic, and clinical factors that can potentially influence TTR. In total,422 patients using warfarin were investigated. Glibenclamide co-medication and presence of CYP2C9*2 and/or *3 alleles were associated with higher TTR, while amiodarone, acetaminophen and verapamil co-medication were associated with lower TTR. Our data suggest that TTR is influenced by co-medication and genetic factors. Thus, individuals in use of glibenclamide may need a more careful monitoring and genetic testing (CYP2C9*2 and/or *3 alleles) may improve the anticoagulation management. In addition, in order to reach and maintain the INR in the target for a longer period, it is better to discuss dose adjustment in office instead of by telephone assessment. Other studies are needed to confirm these results and to find more variables that could contribute to this important parameter.

Integrated CYP2D6 interrogation for multiethnic copy number and tandem allele detection

AIM

To comprehensively interrogate CYP2D6 by integrating genotyping, copy number analysis and novel strategies to identify CYP2D6*36 and characterize CYP2D6 duplications.

METHODS

Genotyping of 16 CYP2D6 alleles, multiplex ligation-dependent probe amplification (MLPA) and CYP2D6*36 and duplication allele-specific genotyping were performed on 427 African-American, Asian, Caucasian, Hispanic, and Ashkenazi Jewish individuals.

RESULTS

A novel PCR strategy determined that almost half of all CYP2D6*10 (100C>T) alleles are actually *36 (isolated or in tandem with *10) and all identified duplication alleles were characterized. Integrated results from all testing platforms enabled the refinement of genotype frequencies across all studied populations.

CONCLUSION

The polymorphic CYP2D6 gene requires comprehensive interrogation to characterize allelic variation across ethnicities, which was enabled in this study by integrating multiplexed genotyping, MLPA copy number analysis, novel PCR strategies and duplication allele-specific genotyping.

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.

Multi-ethnic SULT1A1 copy number profiling with multiplex ligation-dependent probe amplification

AIM

To develop a SULT1A1 multiplex ligation-dependent probe amplification assay and to investigate multi-ethnic copy number variant frequencies.

METHODS

A novel multiplex ligation-dependent probe amplification assay was developed and tested on 472 African-American, Asian, Caucasian, Hispanic and Ashkenazi Jewish individuals.

RESULTS

The frequencies of atypical total copy number (i.e., greater or less than two) were 38.7% for Hispanics, 38.9% for Ashkenazi Jewish, 43.2% for Caucasians, 53.6% for Asians and 64.1% for African-Americans. Heterozygous SULT1A1 deletion carriers (slow sulfators) were most common among Caucasians (8.4%), whereas African-Americans had the highest frequencies of three or more copies (rapid sulfators; 60.9%).

CONCLUSION

Different ethnic and racial populations have varying degrees of SULT1A1-mediated sulfation activity, which warrants further research and that may have utility for drug response prediction among SULT1A1-metabolized medications.

Institutional profile: translational pharmacogenomics at the Icahn School of Medicine at Mount Sinai

For almost 50 years, the Icahn School of Medicine at Mount Sinai has continually invested in genetics and genomics, facilitating a healthy ecosystem that provides widespread support for the ongoing programs in translational pharmacogenomics. These programs can be broadly cataloged into discovery, education, clinical implementation and testing, which are collaboratively accomplished by multiple departments, institutes, laboratories, companies and colleagues. Focus areas have included drug response association studies and allele discovery, multiethnic pharmacogenomics, personalized genotyping and survey-based education programs, pre-emptive clinical testing implementation and novel assay development. This overview summarizes the current state of translational pharmacogenomics at Mount Sinai, including a future outlook on the forthcoming expansions in overall support, research and clinical programs, genomic technology infrastructure and the participating faculty.

Sequencing the CYP2D6 gene: from variant allele discovery to clinical pharmacogenetic testing

CYP2D6 is one of the most studied enzymes in the field of pharmacogenetics. The CYP2D6 gene is highly polymorphic with over 100 catalogued star (*) alleles, and clinical CYP2D6 testing is increasingly accessible and supported by practice guidelines. However, the degree of variation at the CYP2D6 locus and homology with its pseudogenes make interrogating CYP2D6 by short-read sequencing challenging. Moreover, accurate prediction of CYP2D6 metabolizer status necessitates analysis of duplicated alleles when an increased copy number is detected. These challenges have recently been overcome by long-read CYP2D6 sequencing; however, such platforms are not widely available. This review highlights the genomic complexities of CYP2D6, current sequencing methods and the evolution of CYP2D6 from allele discovery to clinical pharmacogenetic testing.

Genome-wide association study of warfarin maintenance dose in a Brazilian sample

AIM

Extreme discordant phenotype and genome-wide association (GWA) approaches were combined to explore the role of genetic variants on warfarin dose requirement in Brazilians.

METHODS

Patients receiving low (≤ 20 mg/week; n = 180) or high stable warfarin doses (≥ 42.5 mg/week; n = 187) were genotyped with Affymetrix Axiom(®) Biobank arrays. Imputation was carried out using data from the combined 1000 Genomes project.

RESULTS

Genome-wide signals (p ≤ 5 × 10(-8)) were identified in the well-known VKORC1 (lead SNP, rs749671; OR: 20.4; p = 1.08 × 10(-33)) and CYP2C9 (lead SNP, rs9332238, OR: 6.8 and p = 4.4 × 10(-13)) regions. The rs9332238 polymorphism is in virtually perfect LD with CYP2C9*2 (rs1799853) and CYP2C9*3 (rs1057910). No other genome-wide significant regions were identified in the study.

CONCLUSION

We confirmed the important role of VKORC1 and CYP2C9 polymorphisms in warfarin dose. Original submitted 14 January 2015; Revision submitted 26 May 2015.

Influence of VKORC1 gene polymorphisms on the effect of oral vitamin K supplementation in over-anticoagulated patients

Significant inter-individual variability on the effect of vitamin K to reverse overanticoagulation has been identified. Genetic polymorphisms of the vitamin K epoxide reductase complex subunit 1 (VKORC1) gene might explain in part this variability. The objective of this study was to evaluate the influence of VKORC1 -1639G>A and 3730G>A polymorphisms on the effect of oral vitamin K supplementation in overanticoagulated patients. We performed an interventional trial of oral vitamin K supplementation in over-anticoagulated outpatients (international normalized ratio [INR] ≥ 4). Subjects received vitamin K (2.5-5.0 mg) according to baseline INR and were genotyped by real time polymerase chain reaction (PCR). INR values were determined at 3, 6, 24 and 72 h after supplementation. We evaluated 33 outpatients, 61 % were males, with a mean age of 62 ± 12 years old. There was a significant decrease in INR values over time for both polymorphisms after oral vitamin K. At 3 h after supplementation, patients carrying the G allele for the -1639G>A polymorphism had a greater decrease in INR values compared to AA patients (p < 0.05 for difference among groups; p < 0.001 for time variation; p = 0.001 for time × group interaction), with differences of -1.01 for GG versus AA (p = 0.003) and -0.84 for GA versus AA (p = 0.024). Mean INR value at 24 h was 1.9 ± 0.6 and at 72 h was 2.1 ± 0.7, with no differences among genotypes. No significant interaction was identified between the 3730G>A polymorphism and vitamin K supplementation. Our study indicated that the VKORC1 -1639G>A polymorphism plays a role in the response to acute vitamin K supplementation in over-anticoagulated patients, with faster decrease of INR value in patients carrying the G allele.

Pharmacogenetics of coumarin anticoagulants in Brazilians

INTRODUCTION

Coumarin vitamin K antagonists are the mainstay of anticoagulant therapy in atrial fibrillation, prosthetic heart valves and thromboembolic conditions. Concerns with these drugs include large inter-individual variability in dose requirements, narrow therapeutic index and need to monitor prothrombin time repeatedly.

AREAS COVERED

Pharmacogenetic studies and dosing algorithms for warfarin and phenprocoumon.

EXPERT OPINION

Gene candidate studies in Brazilian patients verified consistently the association of warfarin and pheprocoumon stable dose requirements with CYP2C9 and VKORC1 polymorphisms, and minor or no influence of other pharmacogenes (e.g., CYP4F2 and F7). The predictive power of warfarin and phenprocoumon dosing algorithms developed for Brazilians compares favorably with those reported for other populations. A warfarin dosing algorithm derived for an admixed cohort performed equally well in self-reported White and Black patients, in marked contrast with the considerably poorer performance of other warfarin algorithms in patients of African descent compared to those of European ancestry. This discrepancy is ascribed to the extensive European/African admixture among Brazilians. Prospective studies of clinical utility of coumarin dosing algorithms, in the context of the Brazilian Public Health System, would represent an important counterpart to recently published trials in European and North American cohorts with predominant or exclusive European ancestry.

Pharmacogenomics of warfarin in populations of African descent

Warfarin is the most commonly prescribed oral anticoagulant worldwide despite its narrow therapeutic index and the notorious inter- and intra-individual variability in dose required for the target clinical effect. Pharmacogenetic polymorphisms are major determinants of warfarin pharmacokinetic and dynamics and included in several warfarin dosing algorithms. This review focuses on warfarin pharmacogenomics in sub-Saharan peoples, African Americans and admixed Brazilians. These 'Black' populations differ in several aspects, notably their extent of recent admixture with Europeans, a factor which impacts on the frequency distribution of pharmacogenomic polymorphisms relevant to warfarin dose requirement for the target clinical effect. Whereas a small number of polymorphisms in VKORC1 (3673G > A, rs9923231), CYP2C9 (alleles *2 and *3, rs1799853 and rs1057910, respectively) and arguably CYP4F2 (rs2108622), may capture most of the pharmacogenomic influence on warfarin dose variance in White populations, additional polymorphisms in these, and in other, genes (e.g. CALU rs339097) increase the predictive power of pharmacogenetic warfarin dosing algorithms in the Black populations examined. A personalized strategy for initiation of warfarin therapy, allowing for improved safety and cost-effectiveness for populations of African descent must take into account their pharmacogenomic diversity, as well as socio-economical, cultural and medical factors. Accounting for this heterogeneity in algorithms that are 'friendly' enough to be adopted by warfarin prescribers worldwide requires gathering information from trials at different population levels, but demands also a critical appraisal of racial/ethnic labels that are commonly used in the clinical pharmacology literature but do not accurately reflect genetic ancestry and population diversity.

Relationship between polymorphisms in thrombophilic genes and preeclampsia in a Brazilian population

The purpose of this study was to evaluate the thrombophilic genes in pregnant women with and without preeclampsia independently or in combination. In a prospective case-control study, we investigated four polymorphisms in thrombophilic genes in 75 women with mild or severe preeclampsia and 145 women with normal pregnancy. The genotype frequencies were assessed and the odds ratio (OR) calculated. When we analyzed the polymorphisms independently and the development of preeclampsia, no association was observed [methylenetetrahydrofolate reductase (MTHFR) 677TT genotype, OR 2.07, 95% confidence interval (CI) 0.99-4.30; prothrombin mutation (F II) (GA or AA genotypes) OR 8.11, 95% CI 0.89-73.92; factor V Leiden (FV Leiden) OR 3.94, 95% CI 0.35-44.23; plasminogen activator inhibitor (PAI-1) 4G/4G genotype, OR 1.63, 95% CI 0.87-3.05] not even with severe preeclampsia subgroup analysis. However, when we investigated a possible interaction among these polymorphisms on the development of the preeclampsia, the OR for having one risk genotype, one or two genotype risk factors and two genotype risk factors compared to those without genotype risk factors were 1.97 (95% CI 1.08-3.59), 2.21 (95% CI 1.25-3.92) and 4.27 (95% CI 1.3-13.9), respectively. In conclusion, in the population analyzed, the presence of the genotype risk factors alone does not seem to be associated with the development of preeclampsia even in the severe presentation form. However, an interaction among the MTHFR, F II, FV and PAI-1 gene polymorphisms on the development of the preeclampsia was indicated.