Pharmacogenomic implications of the differential distribution of CYP3A5 metabolic phenotypes among Latin American populations.

This study shows that the distribution of CYP3A5 alleles (*1, *3, *6 and *7) and genotype-predicted CYP3A5 phenotypes vary significantly across Latin American cohorts (Brazilians and the One Thousand Genomes Admixed American superpopulation), as well as among subcohorts comprising individuals with the highest proportions of Native, European or sub-Saharan African ancestry. Differences in biogeographical ancestry across the study groups are the likely explanation for these results. The differential distribution of CYP3A5 phenotypes has major pharmacogenomic implications, affecting the proportion of individuals carrying high risk CYP3A5 phenotypes for the immunosuppressant tacrolimus and the number of patients that would need to be genotyped to prevent acute rejection in kidney transplant recipients under tacrolimus treatment.

NUDT15 and TPMT polymorphisms in three distinct native populations of the Brazilian Amazon.

This is the first report of the distribution of TPMT and NUDT15 single nucleotide polymorphisms and metabolic phenotypes associated with cytotoxicity of thiopurine drugs, in indigenous groups of Brazilian Amazon: Munduruku, Paiter-Suruí and Yanomami. The minor allele frequency (MAF) of NUDT15 rs116855232 did not differ significantly across the groups; TPMT rs1800462 was absent, while rs1800460 and rs1142345 were in strong linkage disequilibrium, and 10- and 30-fold more common in Paiter-Suruí. Indeed, the MAFs in Paiter-Surui (0.193 and 0.188) are the largest report globally. The distribution of combined NUDT15/TPMT metabolic phenotypes differed significantly (p < 0.0001) and largely (Cramér´s V = 0.37) across cohorts. This has important pharmacogenetic implications: the Clinical Pharmacogenetics Implementation Consortium recommendations to reduce or consider reduction of thiopurine dose applies to 4.4% Yanomami, 5.6% Munduruku, versus 41% Paiter-Suruí. The proportion of Paiter-Suruí at risk of thiopurine intolerance is 3- to 4-fold higher than any other population worldwide.

Intestinal P-gp activity is reduced in postmenopausal women under breast cancer therapy.

Intestinal P-glycoprotein (P-gp) activity plays a crucial role in modulating the oral bioavailability of its substrates. Fexofenadine has commonly been used as a P-gp probe, although it is important to note the involvement of other drug transporters like, OATP1B1, OATP1B3, and OATP2B1. In vitro studies demonstrated an upregulation of P-gp protein in response to exposure to pregnancy-related hormones. The objective of this study was to investigate how intestinal P-gp activity is impacted by menopausal status. This study sampled fexofenadine plasma concentrations over 0-12 h after probe drug administration from two groups of patients with breast cancer: premenopausal (n = 20) and postmenopausal (n = 20). Fexofenadine plasma concentrations were quantified using liquid-chromatography tandem mass spectrometry. Area under the plasma concentration-time curve from zero to infinity (AUCinf ) was calculated through limited sampling strategies equation. Multiple linear regression was applied on AUCinf , maximum plasma concentration (Cmax ), and time to Cmax . Postmenopausal patients showed a significant increase in Cmax (geometric mean and 95% confidence interval [CI] 143.54, 110.95-176.13 vs. 223.54 ng/mL, 161.02-286.06 and in AUCinf 685.55, 534.98-878.50 vs. 933.54 ng·h/mL 735.45-1184.99) compared to premenopausal patients. The carriers of the ABCB1 3435 allele T displayed higher Cmax values of 166.59 (95% CI: 129.44-214.39) compared to the wild type at 147.47 ng/mL (95% CI: 111.91-194.34, p = 0.02). In postmenopausal individuals, the decrease in P-gp activity of ~40% may lead to an increased plasma exposure of orally administered P-gp substrates.

Pharmacogenomic implications of the differential distribution of CYP2C9 metabolic phenotypes among Latin American populations.

The CYP2C9 gene encodes the major drug metabolism enzyme CYP2C9. This gene is highly polymorphic, and no-function (CYP2C9*3) plus decreased function (CYP2C9*2, *5, *8 and *11) star alleles (haplotypes) are commonly used to predict CYP2C9 metabolic phenotypes. This study explores the pharmacogenomic implications of the differential distribution of genotype-predicted CYP2C9 phenotypes across Latin American populations. Data from 1,404 individuals from the South American countries Brazil, Colombia and Peru, from Puerto Rico in the Caribbean and from persons with Mexican ancestry living in North America were analysed. The results showed that the distribution of CYP2C9 alleles and diplotypes, and diplotype-predicted CYP2C9 phenotypes vary significantly across the distinct country cohorts, as well as among self-identified White, Brown and Black Brazilians. Differences in average proportions of biogeographical ancestry across the study groups, especially Native American and African ancestry, are the likely explanation for these results. The differential distribution of genotype-predicted CYP2C9 phenotypes has potentially clinically-relevant pharmacogenomic implications, through its influence on the proportion of individuals at high risk for adverse response to medications that are CYP2C9 substrates, the proportion on individuals with CPIC therapeutic recommendations for dosing and choice of nonsteroidal antinflammatory drugs (NSAIDs) and the number of individuals that need to be genotyped in order to prevent adverse effects of NSAIDs. Collectively, these findings are likely to impact the perceived benefits, cost-effectiveness and clinical adoption of pharmacogenomic screening for drugs that are predominantly metabolized by CYP2C9.

Characterization of pharmacogenomic variants in a Brazilian admixed cohort of elderly individuals based on whole-genome sequencing data.

Introduction: Research in the field of pharmacogenomics (PGx) aims to identify genetic variants that modulate response to drugs, through alterations in their pharmacokinetics (PK) or pharmacodynamics (PD). The distribution of PGx variants differs considerably among populations, and whole-genome sequencing (WGS) plays a major role as a comprehensive approach to detect both common and rare variants. This study evaluated the frequency of PGx markers in the context of the Brazilian population, using data from a population-based admixed cohort from Sao Paulo, Brazil, which includes variants from WGS of 1,171 unrelated, elderly individuals. Methods: The Stargazer tool was used to call star alleles and structural variants (SVs) from 38 pharmacogenes. Clinically relevant variants were investigated, and the predicted drug response phenotype was analyzed in combination with the medication record to assess individuals potentially at high-risk of gene-drug interaction. Results: In total, 352 unique star alleles or haplotypes were observed, of which 255 and 199 had a frequency < 0.05 and < 0.01, respectively. For star alleles with frequency > 5% (n = 97), decreased, loss-of-function and unknown function accounted for 13.4%, 8.2% and 27.8% of alleles or haplotypes, respectively. Structural variants (SVs) were identified in 35 genes for at least one individual, and occurred with frequencies >5% for CYP2D6, CYP2A6, GSTM1, and UGT2B17. Overall 98.0% of the individuals carried at least one high risk genotype-predicted phenotype in pharmacogenes with PharmGKB level of evidence 1A for drug interaction. The Electronic Health Record (EHR) Priority Result Notation and the cohort medication registry were combined to assess high-risk gene-drug interactions. In general, 42.0% of the cohort used at least one PharmGKB evidence level 1A drug, and 18.9% of individuals who used PharmGKB evidence level 1A drugs had a genotype-predicted phenotype of high-risk gene-drug interaction. Conclusion: This study described the applicability of next-generation sequencing (NGS) techniques for translating PGx variants into clinically relevant phenotypes on a large scale in the Brazilian population and explores the feasibility of systematic adoption of PGx testing in Brazil.

Distribution of a novel CYP2C haplotype in Native American populations.

The CYP2C19 gene, located in the CYP2C cluster, encodes the major drug metabolism enzyme CYP2C19. This gene is highly polymorphic and no-function (CYP2C19*2 and CYP2C19*3), reduced function (CYP2C19*9) and increased function (CYP2C19*17) star alleles (haplotypes) are commonly used to predict CYP2C19 metabolic phenotypes. CYP2C19*17 and the genotype-predicted rapid (RM) and ultrarapid (UM) CYP2C19 metabolic phenotypes are absent or rare in several Native American populations. However, discordance between genotype-predicted and pharmacokinetically determined CYP2C19 phenotypes in Native American cohorts have been reported. Recently, a haplotype defined by rs2860840T and rs11188059G alleles in the CYP2C cluster has been shown to encode increased rate of metabolism of the CYP2C19 substrate escitalopram, to a similar extent as CYP2C19*17. We investigated the distribution of the CYP2C:TG haplotype and explored its potential impact on CYP2C19 metabolic activity in Native American populations. The study cohorts included individuals from the One Thousand Genomes Project AMR superpopulation (1 KG_AMR), the Human Genome Diversity Project (HGDP), and from indigenous populations living in Brazil (Kaingang and Guarani). The frequency range of the CYP2C:TG haplotype in the study cohorts, 0.469 to 0.598, is considerably higher than in all 1 KG superpopulations (range: 0.014-to 0.340). We suggest that the high frequency of the CYP2C:TG haplotype might contribute to the reported discordance between CYP2C19-predicted and pharmacokinetically verified CYP2C19 metabolic phenotypes in Native American cohorts. However, functional studies involving genotypic correlations with pharmacokinetic parameters are warranted to ascertain the importance of the CYP2C:TG haplotype.

Assessment of small in-frame indels and C-terminal nonsense variants of BRCA1 using a validated functional assay.

BRCA1 (Breast Cancer 1, early onset) is linked to breast and ovarian cancer predisposition. Still, the risks conferred by a significant portion of BRCA1 variants identified in the population remains unknown. Most of these variants of uncertain significance are missense alterations. However, the functional implications of small in-frame deletions and/or insertions (indels) are also difficult to predict. Our group has previously evaluated the functional impact of 347 missense variants using an extensively validated transcriptional activity assay. Here we show a systematic assessment of 30 naturally occurring in-frame indels located at the C-terminal region of BRCA1. We identified positions sensitive and tolerant to alterations, expanding the knowledge of structural determinants of BRCA1 function. We further designed and assessed the impact of four single codon deletions in the tBRCT linker region and six nonsense variants at the C-terminus end of BRCA1. Amino acid substitutions, deletions or insertions in the disordered region do not significantly impact activity and are not likely to constitute pathogenic alleles. On the other hand, a sizeable fraction of in-frame indels at the BRCT domain significantly impact function. We then use a Bayesian integrative statistical model to derive the probability of pathogenicity for each variant. Our data highlights the importance of assessing the impact of small in-frame indels in BRCA1 to improve risk assessment and clinical decisions for carriers.

Functional Restoration of BRCA1 Nonsense Mutations by Aminoglycoside-Induced Readthrough.

BRCA1 is a major tumor suppressor that functions in the accurate repair of DNA double-strand breaks via homologous recombination (HR). Nonsense mutations in BRCA1 lead to inactive truncated protein products and are associated with high risk of breast and ovarian cancer. These mutations generate premature termination codons (PTCs). Different studies have shown that aminoglycosides can induce PTC suppression by promoting stop codon readthrough and restoring full-length (FL) protein expression. The use of these compounds has been studied in clinical trials for genetic diseases such as cystic fibrosis and Duchenne muscular dystrophy, with encouraging results. Here we show proof-of-concept data demonstrating that the aminoglycoside G418 can induce BRCA1 PTC readthrough and restore FL protein synthesis and function. We first demonstrate that G418 treatment restores BRCA1 FL protein synthesis in HCC1395, a human breast tumor cell line carrying the R1751X mutation. HCC1395 cells treated with G418 also recover HR DNA repair and restore cell cycle checkpoint activation. A set of naturally occurring BRCA1 nonsense variants encoding different PTCs was evaluated in a GFP C-terminal BRCA1 construct model and BRCA1 PTC readthrough levels vary depending on the stop codon context. Because PTC readthrough could generate FL protein carrying pathogenic missense mutations, variants representing the most probable acquired amino acid substitutions in consequence of readthrough were functionally assessed by a validated transcription activation assay. Overall, this is the first study that evaluates the readthrough of PTC variants with clinical relevance in the breast and ovarian cancer-predisposing gene BRCA1.

Single Nucleotide Variants as Proxies for HLA-A*31:01 in Native American Populations.

Carbamazepine triggers dermatologic hypersensitivity reactions, associated with specific human leukocyte antigens (HLAs), especially HLA-B*15:02 and HLA-A*31:01. Previous efforts to identify single nucleotide variants (SNVs) with high predictive value as HLA proxies, revealed that rs1061235 and rs17179220 fulfill these requirements for HLA-A*31:01 in some but not all populations. Herein we explored the predictive performance of rs1061235 and rs17179220 as HLA-A*31:01 tags in populations of Native American ancestry, which are largely underrepresented in pharmacogenomic studies. The study cohorts comprised the overall Admixed American superpopulation of the 1000 Genomes Project (1 KG_AMR), a subcohort of individuals with predominant Native American ancestry (1 KG_NAT), the Native American population of the Human Genome Diversity Project (HGDP), plus Kaingang (KRC) and Guarani (GRC and GKW) adults from indigenous reservation areas in Brazil. The diversity of cohorts is reflected in the range of frequencies of HLA-A*31:01 (0.02-0.65), rs1061235 (0.03-0.13) and rs17179220 (0.12-0.66), as well as in the predictive performance of these SNVs as HLA-A*31:01 proxies. NPV (negative predictive value), the metric of primary interest for pharmacogenetic-informed carbamazepine prescription was maximal (NPV = 1.0) for both SNVs in 1 KG_AMR and 1 KG_NAT, for rs17179220, but not rs1061235 (NPV = 0.91) in HGDP, and for rs17179220 in GWK, but not GRC (NPV = 0.88) or KRC (NPV = 0.80). Collectively, the data support the notion that rs1061235 and rs17179220 are not optimal proxies for HLA-A*31:01 across populations of Native American ancestry.

Pharmacogenetic predictors of variability in efavirenz pharmacokinetics in an admixed Brazilian HIV cohort.

AIMS: To investigate the influence of pharmacogenetic polymorphisms on efavirenz (EFV) exposure and metabolism in HIV-infected Brazilians under treatment with EFV-containing antiretroviral (ART) regimens. METHODS: HIV-positive adults (n = 82) on stable ART regimens containing 600 mg EFV once daily for at least 6 months were recruited at 2 university hospitals. Blood samples collected at mid-dose interval were used to quantify the plasma concentrations of EFV (denoted [EFV]), its major metabolite 8-OH-EFV ([8-OH-EFV]) and [8-OH-EFV]/[EFV] metabolic ratio, and to genotype single nucleotide polymorphisms in CYP2B6 (rs3745274, c.516G > T; rs28399499, c.983 T > C) and ABCB1 (rs3842, c.4036G > A). CYP2B6 metabolic phenotypes were inferred from the CYP2B6 diplotypes. Linear regression modelling was applied to identify sociodemographic, clinical and pharmacogenetic predictors of [EFV] and [8-OH-EFV]/[EFV] metabolic ratio. RESULTS: Wide (50-fold) interindividual variation in [EFV], [8-OH-EFV] and [8-OH-EFV]/[EFV] was observed; 69.5% of participants had [EFV] within the nominal therapeutic range (1000-4000 ng/mL), while 19.5 and 11.0% had [EFV] below and above this range, respectively. Multiple regression modelling retained only CYP2B6 metabolic phenotypes or the combined rs3745274 and rs28399499 genotypes, as significant predictors of [EFV] and [8-OH-EFV]/[EFV]. CONCLUSION: EFV exposure and disposition varied widely among HIV-infected Brazilians under stable treatment with EFV-containing ART regimens. About 1/10 of the participants had [EFV] exceeding nominal supratherapeutic concentration (4000 ng/mL), but reported tolerance to the ARV regimens, while 1/5 of participants had nominal subtherapeutic [EFV] (<1000 ng/mL) but adequate virological response. Genotype for the 2 CYP2B6 single nucleotide polymorphisms studied explained 48% of variation in [EFV] and 35% of variation in [8-OH-EFV]/[EFV].

Pharmacogenetic differentiation across Latin America.

Aim: To explore the pharmacogenetic differentiation across Latin American populations, using the fixation index statistics (FST). Materials & methods: FST analyses were applied to 1519 pharmacogenetic markers in the 1000 Genomes admixed American superpopulation (1KG_AMR) and an admixed Brazilian sample. Results: Allele-specific FST values for the overall cohort point to little overall pharmacogenetic differentiation (average FST = 0.017); however, moderate differentiation (FST = 0.05-0.15) was observed for 83 markers, while large differentiation (FST = 0.15-0.25) was restricted to three markers. Pairwise FST analysis identified three markers with very large differentiation (FST >0.25). Conclusion: The present study verifies and extends previous reports of little overall pharmacogenetic divergence across Latin America, although a number of markers display substantial differentiation.

Impact of CYP2D6 Genetic Variation on Radical Cure of Plasmodium vivax Malaria.

Plasmodium vivax (P. vivax) is the most widespread human malaria parasite, with 2.5 billion people at risk of infection worldwide. P. vivax forms liver hypnozoites, which trigger further symptomatic episodes (relapses) weeks or months after the initial episode. Radical cure of vivax malaria requires hypnozoitocide therapy to prevent relapses. The two US Food and Drug Administration (FDA)-approved hypnozoiticides for human use, primaquine, and tafenoquine, are pro-drugs, that require in vivo conversion into metabolites with redox activity. This mini-review focuses on the association between CYP2D6-mediated hydroxylation and hypnozoitocide efficacy of primaquine and tafenoquine. Studies in murine models show that the antimalarial activity of primaquine and tafenoquine is abolished by CYP2D knock-out and partially restored by knock-in of humanized CYP2D6. Human studies explored the impact of CYP2D6 genetic variation and genotype-inferred CYP2D6 phenotype on anti-relapse efficacy. Most, but not all, studies with primaquine report higher rates of relapse in patients with decreased CYP2D6 activity (activity scores (AS) ≤ 1) compared to normal activity (AS ≥ 1.5). Potential factors for discordance among studies include risk of reinfection in endemic areas, adherence to primaquine-treatment, assignment of CYP2D6 phenotypes based on CYP2D6 polymorphism and choice of AS values for dichotomizing the study cohorts. Tafenoquine anti-relapse efficacy did not differ between patients with AS < 1 vs. AS ≥ 1.5 in 2 studies. Absence/small number of poor CYP2D6 metabolizers in AS ≤ 1 groups, combined with lesser dependence of tafenoquine on CYP2D6-mediated conversion into active redox metabolites may account for this result. Additional tafenoquine studies with larger representation of poor CYP2D6 metabolizers are warranted.

Identification of ancestry proportions in admixed groups across the Americas using clinical pharmacogenomic SNP panels.

We evaluated the performance of three PGx panels to estimate biogeographical ancestry: the DMET panel, and the VIP and Preemptive PGx panels described in the literature. Our analysis indicate that the three panels capture quite well the individual variation in admixture proportions observed in recently admixed populations throughout the Americas, with the Preemptive PGx and DMET panels performing better than the VIP panel. We show that these panels provide reliable information about biogeographic ancestry and can be used to guide the implementation of PGx clinical decision-support (CDS) tools. We also report that using these panels it is possible to control for the effects of population stratification in association studies in recently admixed populations, as exemplified with a warfarin dosing GWA study in a sample from Brazil.