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.

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.

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 two 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 (SNPs) in CYP2B6 (rs3745274, c.516G>T; rs28399499, c.983T>C) and ABCB1 (rs3842, c.4036G>A). CYP2B6 metabolic phenotypes were inferred from the CYP2B6 diplotypes. Linear regression modeling 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 modeling 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 one-tenth of the participants had [EFV] exceeding nominal supratherapeutic concentration (4000 ng/mL), but reported tolerance to the ARV regimens, while one fifth of participants had nominal subtherapeutic [EFV] (<1000 ng/ml) but adequate virological response. Genotype for the two CYP2B6 SNPs studied explained 48% of variation in [EFV] and 35% of variation in [8-OH-EFV]/[EFV].

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 differentiation across Latin America

Aim: To explore the pharmacogenetic differentiation across Latin American populations, using the fixation index statistics (F_ST). Materials & methods: F_ST analyses were applied to 1519 pharmacogenetic markers in the 1000 Genomes Admixed American superpopulation (1KG_AMR) and an admixed Brazilian sample. Results: Allele-specific F_ST values for the overall cohort point to little overall pharmacogenetic differentiation (average F_ST = 0.017); however, moderate differentiation (F_ST = 0.05-0.15) was observed for 83 markers, while large differentiation (F_ST = 0.15-0.25) was restricted to three markers. Pairwise F_ST analysis identified three markers with very large differentiation (F_ST >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.

Impact of Plasmodium vivax malaria and antimalarial treatment on cytochrome P450 activity in Brazilian patients

AIMS

To investigate the impact of Plasmodium vivax malaria and chloroquine-primaquine chemotherapy on CYP2D6 and CYP2C19 activity in patients from the Brazilian Amazon.

METHODS

Adult patients (n = 30) were given subtherapeutic doses of CYP2D6 and CYP2C19 phenotypic probes metoprolol (10 mg) and omeprazole (2 mg) in three different stages of vivax malaria illness: acute disease (study phase 1), post chemotherapy (phase 2) and convalescence (stage 3). Plasma concentrations of probes and CYP-hydroxylated metabolites (α-OH metoprolol and 5-OH omeprazole) were measured using LC/MS/MS. Two pharmacokinetic metrics were used to estimate CYP activity: (a) ratio of plasma concentrations of probe/metabolite at 240 minutes after administration of the probes and (b) ratio of areas under the time-concentration curves for probe/metabolite (AUC_0-12h ). For statistical analysis, the pharmacokinetic metrics were normalized to the respective values in phase 3. Taqman assays were used for CYP2D6 and CYP2C19 genotyping. Cytokines levels were measured using cytometric bead array.

RESULTS

Both pharmacokinetic metrics for metoprolol and omeprazole, and plasma concentrations of cytokines IL-6, IL-8 and IL-10 varied significantly across the three study phases (ANOVA P < 0.0001). Post hoc tests showed greater metoprolol:α-OH metoprolol ratios in phases 1 and 2 compared to phase 3, larger omeprazole:5-OH omeprazole ratios in phase 1 than in phases 2 and 3, and higher circulating IL-6, IL-8 and IL-10 in phase 1 than in phases 2 and 3.

CONCLUSION

P. vivax malaria and treatment altered CYP2D6 and CYP2C19 metabolic phenotypes. CYP2C19 inhibition is attributed to a higher level of circulating proinflammatory cytokines, while suppression of CYP2D6 is ascribed mainly to chloroquine exposure.

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.

Implementation of a pharmacogenomic program in a Brazilian public institution

This narrative review describes implementation, current status and perspectives of a pharmacogenomic (PGx) program at the Brazilian National Cancer Institute (INCA), targeting the cancer chemotherapeutic drugs – fluoropyrimidines, irinotecan and thiopurines. This initiative, designed as a research project, was supported by a grant from the Brazilian Ministry of Health. A dedicated task force developed standard operational procedures from recruitment of patients to creating PGx reports with dosing recommendations, which were successfully applied to test 100 gastrointestinal cancer INCA outpatients and 162 acute lymphoblastic leukemia pediatric patients from INCA and seven other hospitals. The program has been subsequently expanded to include gastrointestinal cancer patients from three additional cancer treatment centers. We anticipate implementation of routine pre-emptive PGx testing at INCA but acknowledge challenges associated with this transition, such as continuous financing support, availability of trained personnel, adoption of the PGx-informed prescription by the clinical staff and, ultimately, evidence of cost–effectiveness.

4th ESPT Conference: pharmacogenomics and personalized medicine - research progress and clinical implementation

The Fourth European Society of Pharmacogenomics and Personalized Therapy biennial conference was organized in collaboration with the Italian Society of Personalized Medicine (SIMeP) and was held at Benedictine Monastery of San Nicolò l'Arena in Catania, Sicily (Italy) on 4-7 October 2017. The congress addressed the research progress and clinical implementation in pharmacogenomics and personalized medicine. The Fourth European Society of Pharmacogenomics and Personalized Therapy congress brought together leading international scientists and healthcare professionals actively working in the fields of pharmacogenomics and personalized therapy. Altogether, 25 speakers in 15 session comprehensively covered broad spectrum of pharmacogenetics and pharmacogenomics research, clinical applications in different clinical disciplines attended by 270 delegates.

A background subtraction approach for determination of endogenous cortisol and 6β-hydroxycortisol in urine by UPLC-MS/MS with application in a within-day variability study in HIV-infected pregnant women

Different methods have been used for CYP3A phenotyping, such as probe drugs or the urinary index 6β-hydroxycortisol/cortisol ratio (6β-OHF:C). This work describes a simple and affordable method for the simultaneous determination of the endogenous compounds cortisol and 6β-hydroxycortisol in urine using a background subtraction approach. The method was applied to investigate the CYP3A activity in HIV-infected pregnant women (n = 9) in the third trimester and postpartum periods. Also, the within-day variability in the 6β-OHF:C index was also evaluated. The sample preparation consists of a pre-cleanup with acetonitrile followed by liquid-liquid extraction with ethyl acetate. The analytes were resolved by employing an Acquity UPLC®BEH C18 column with a mobile phase that consisted of a mixture of acetonitrile containing 0.1% formic acid and 0.1% formic acid in gradient mode. The method presented linearities of 1-1.000 ng/mL and 2-1.000 ng/mL for C and 6β-OHF, respectively, and presented acceptable precision and accuracy. Qualitative and quantitative matrix effects tests were also performed. A high 6β-OHF:C within-day variability was observed in both phases. In the third trimester period, the 6β-OHF:C ranged from 2.57 to 51.69, with a mean ± standard deviation (SD) of 15.12 ± 5.41 (n = 9). Similar values were obtained in the postpartum period, with 6β-OHF:C ranging from 3.48 to 44.54 with a mean ± SD of 14.37 ± 5.73 (n = 7). Even though the 6β-OHF:C is a non-invasive index for CYP3A phenotyping, its use is susceptible to high within-day variability.

Limited Sampling Modeling for Estimation of Phenotypic Metrics for CYP Enzymes and the ABCB1 Transporter Using a Cocktail Approach

Plasma concentration data points (n = 2,640) from 16 healthy adults were used to develop and validate limited sampling strategies (LSS) for estimation of phenotypic metrics for CYP enzymes and the ABCB1 transporter, using a cocktail of subtherapeutic doses of the selective probes caffeine (CYP1A2), metoprolol (CYP2D6), midazolam (CYP3A), losartan (CYP2C9), omeprazole (CYP2C19), and fexofenadine (ABCB1). All-subsets linear regression modelling was applied to estimate the AUC_0–12h for caffeine, fexofenadine, and midazolam, and the AUC_0–12h ratio of metoprolol: α-OH metoprolol and omeprazole:5-OH omeprazole. LSS-derived metrics were compared with the parameters’ ‘best estimates’ obtained by non-compartmental analysis using all plasma concentration data points. The correlation coefficient (R²) was used to identify the LSS equations that provided the best fit for n timed plasma samples, and the jack-knife statistics was used as an additional validation procedure for the LSS models. Single time-point LSS models provided R² values greater than 0.95 (R² > 0.95) for the AUC_0–12h ratio of metoprolol:α-OH metoprolol and omeprazole:5-OH omeprazole, whereas 2 time-point models were required for R² > 0.95 for the AUC_0–12h of caffeine, fexofenadine, and midazolam. Increasing the number of sampling points to three led to minor increases in R² and/or the bias or prediction of the estimates. In conclusion, the LSS models provided accurate prediction of phenotypic indices for CYP1A2, CYP2C19, CYP2D6, CYP3A, and ABCB1, when using subtherapeutic doses of selective probes for these enzymes and transporter.

Hormonal status affects plasma exposure of tamoxifen and its main metabolites in tamoxifen-treated breast cancer patients

Background

Tamoxifen is considered a prodrug of its active metabolite endoxifen, which is dependent on the CYP2D6 and CYP3A enzymes. Tamoxifen pharmacokinetic variability influences endoxifen exposure and, consequently, its clinical outcome. This study investigated the impact of hormonal status on the pharmacokinetics of tamoxifen and its metabolites in TAM-treated breast cancer patients.

Methods

TAM-treated breast cancer patients (n = 40) previously believed to have CYP3A activity within the normal range based on oral midazolam and phenotyped as CYP2D6 normal metabolizers using oral metoprolol were divided into two groups according to premenopausal (n = 20; aged 35–50 years) or postmenopausal (n = 20; aged 60–79 years) status. All patients were treated with 20 mg/day tamoxifen for at least three months. Serial plasma samples were collected within the 24 h dose interval for analysis of unchanged tamoxifen, endoxifen, 4-hydroxytamoxifen and N-desmethyltamoxifen quantified by LC-MS/MS. CYP activities were assessed using midazolam apparent clearance (CYP3A) and the metoprolol/alfa-hydroxymetoprolol plasma metabolic ratio (CYP2D6). CYP3A4, CYP3A5 and CYP2D6 SNPs and copy number variation were investigated using TaqMan assays.

Results

Postmenopausal status increased steady-state plasma concentrations (Css) of tamoxifen (116.95 vs 201.23 ng/mL), endoxifen (8.01 vs 18.87 ng/mL), N-desmethyltamoxifen (485.16 vs 843.88 ng/mL) and 4-hydroxytamoxifen (2.67 vs 4.11 ng/mL). The final regression models included hormonal status as the only predictor for Css of tamoxifen [β-coef ± SE, p-value (75.03 ± 17.71, p = 0.0001)] and 4-hydroxytamoxifen (1.7822 ± 0.4385, p = 0.0002), while endoxifen Css included hormonal status (8.578 ± 3.402, p = 0.02) and race (11.945 ± 2.836, p = 0.007). For N-desmethyltamoxifen Css, the final model was correlated with hormonal status (286.259 ± 76.766, p = 0.0007) and weight (− 8.585 ± 3.060, p = 0.008).

Conclusion

The premenopausal status was associated with decreased endoxifen plasma concentrations by 135% compared to postmenopausal status. Thus, the endoxifen plasma concentrations should be monitored mainly in the premenopausal period to maintain plasma levels above the efficacy threshold value.

Trial registration

RBR-7tqc7k.

Abstract 1740: Characterization of CDK9 role in the DNA damage response

Introduction: The DNA damage response (DDR) is a well-coordinated pathway capable of sensing and repairing different types of DNA damage. Double-strand break (DSB) lesions activate the ATM-CHK2 axis, which is responsible for the cell cycle arrest, prompting cells for DNA repair. DSBs are primarily repaired through two distinct pathways: homology-directed recombination (HR) and non-homologous end-joining (NHEJ). BRCA1 and 53BP1 play an important role in DDR by orchestrating the decision between HR and NHEJ, but the precise mechanisms regarding both pathways are not entirely understood. Previously, our group identified the putative interaction of the cyclin-dependent kinase 9 (CDK9) with BRCA1 and BARD1 (BRCA1-associated RING domain 1). CDK9 is a component of the positive transcription elongation complex and has been implicated in genome integrity maintenance associated with the replication stress response. More recently, we characterized CDK9 as a DDR player by modulating BRCA1 recruitment to DNA damaged sites and consequently the HR pathway. We also showed that cells lacking CDK9 failure to form RAD51, but not 53BP1, ionizing radiation-induced foci (IRIF) and are characterized by an increased sensitivity to DSBs

Methods and Results: Here we show that the CDK9-BRCA1 interaction occurs through BRCA1 N-terminal region in a BRCA1-BARD1 heterodimer formation-dependent manner. We also identified a CDK9 mutant (E369A) that abrogates the CDK9-BRCA1 interaction (through BRCA1 C-terminal region). Moreover, we characterized CHK2 as a novel CDK9-interaction partner, through ectopic GST-pulldown and constitutive co-immunoprecipitation assays. Further, cells lacking CDK9 are deficient in promoting G2/M arrest after IR-induced DNA damage. We also identified new CDK9 interactions in a DDR-related context using a TAP-MS approach.

Conclusion: Collectively, our data suggest a role for CDK9 in the DDR through DNA damage repair and cell cycle checkpoint regulation

Citation Format: Thales C. Nepomuceno, Guilherme Suarez-Kurtz, Alvaro N. Monteiro, Marcelo A. Carvalho. Characterization of CDK9 role in the DNA damage response [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1740.

Effect of CYP4F2, VKORC1, and CYP2C9 in Influencing Coumarin Dose: A Single-Patient Data Meta-Analysis in More Than 15,000 Individuals

The cytochrome P450 (CYP)4F2 gene is known to influence mean coumarin dose. The aim of the present study was to undertake a meta-analysis at the individual patients level to capture the possible effect of ethnicity, gene-gene interaction, or other drugs on the association and to verify if inclusion of CYP4F2*3 variant into dosing algorithms improves the prediction of mean coumarin dose. We asked the authors of our previous meta-analysis (30 articles) and of 38 new articles retrieved by a systematic review to send us individual patients' data. The final collection consists of 15,754 patients split into a derivation and validation cohort. The CYP4F2*3 polymorphism was consistently associated with an increase in mean coumarin dose (+9% (95% confidence interval (CI) 7-10%), with a higher effect in women, in patients taking acenocoumarol, and in white patients. The inclusion of the CYP4F2*3 in dosing algorithms slightly improved the prediction of stable coumarin dose. New pharmacogenetic equations potentially useful for clinical practice were derived.

Impact of amino acid substitutions at secondary structures in the BRCT domains of the tumor suppressor BRCA1: Implications for clinical annotation

Genetic testing for BRCA1, a DNA repair protein, can identify carriers of pathogenic variants associated with a substantially increased risk for breast and ovarian cancers. However, an association with increased risk is unclear for a large fraction of BRCA1 variants present in the human population. Most of these variants of uncertain clinical significance lead to amino acid changes in the BRCA1 protein. Functional assays are valuable tools to assess the potential pathogenicity of these variants. Here, we systematically probed the effects of substitutions in the C terminus of BRCA1: the N- and C-terminal borders of its tandem BRCT domain, the BRCT-[N-C] linker region, and the α1 and α'1 helices in BRCT-[N] and -[C]. Using a validated transcriptional assay based on a fusion of the GAL4 DNA-binding domain to the BRCA1 C terminus (amino acids 1396-1863), we assessed the functional impact of 99 missense variants of BRCA1. We include the data obtained for these 99 missense variants in a joint analysis to generate the likelihood of pathogenicity for 347 missense variants in BRCA1 using VarCall, a Bayesian integrative statistical model. The results from this analysis increase our understanding of BRCA1 regions less tolerant to changes, identify functional borders of structural domains, and predict the likelihood of pathogenicity for 98% of all BRCA1 missense variants in this region recorded in the population. This knowledge will be critical for improving risk assessment and clinical treatment of carriers of BRCA1 variants.

Pharmacogenomics research and clinical implementation in Brazil

We searched PubMed entries and the Lattes database of Brazilian Pharmacogenetics Network investigators, for pharmacogenetic/genomic (PGx) studies in the Brazilian population, focusing on the drugs and genes included in the Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines. Warfarin was the most extensively studied drug in a PGx context: a genomewide association study targeting warfarin stable dose identified significant signals in VKORC1 and CYP2C9, several PGx dosing algorithms were developed based on these and other genes, and the implications of population admixture on extrapolation of dosing recommendations in the CPIC guidelines were examined. A study in renal transplanted patients disclosed association of CYP3A5*6 and CYP3A5*7 with tacrolimus dosing, which led to addition of these variants to CYP3A5*3 in the CPIC tacrolimus guideline. Studies verified predisposition of HIV-positive carriers of UGT1A1*28 to severe atazanavir-induced hyperbilirubinaemia, intolerance to 5-fluorouracyl in gastrointestinal cancer patients with deleterious DPYD variants, failure of HCV-infected carriers of IFNL3 rs12979860 to obtain a sustained viral response to PEG-IFN-α, and hypersensitivity reactions to abacavir in HIV-positive carriers of HLA-B*57:01. No prospective analyses of drug therapy outcomes or cost-effectiveness assessments of PGx-guided therapy were found. In conclusion, the limited adoption of PGx-informed drug prescription in Brazil reflects combination of recognized barriers to PGx implementation worldwide plus factors specific to the Brazilian population. The latter include rarity/absence of genetic variants on which international PGx guidelines are based (eg HLA-B*15.02 for phenytoin and carbamazepine) and the caveat of extrapolating to the admixed Brazilian population, guidelines based on categorical variables, such as continental ancestry (eg warfarin guidelines), "race" or ethnicity.

Conference report: pharmacogenomics in special populations at WCP2018

The 18th World Congress of Basic and Clinical Pharmacology (WCP2018), coordinated by IUPHAR and hosted by the Japanese Pharmacological Society and the Japanese Society of Clinical Pharmacology and Therapeutics, was held in July 2018 at the Kyoto International Conference Center, in Kyoto, Japan. Having as its main theme 'Pharmacology for the Future: Science, Drug Development and Therapeutics', WCP2018 was attended by over 4500 delegates, representing 78 countries. The present report is an overview of a symposium at WCP2018, entitled Pharmacogenomics in Special Populations, organized by IUPHAR´s Pharmacogenetics/Genomics (PGx) section. The PGx section congregates distinguished scientists from different continents, covering expertise from basic research, to clinical implementation and ethical aspects of PGx, and one of its major activities is the coordination of symposia and workshops to foster exchange of PGx knowledge (https://iuphar.org/sections-subcoms/pharmacogenetics-genomics/). The symposium attracted a large audience to listen to presentations covering various areas of research and clinical adoption of PGx in Oceania, Africa, Latin America and Asia.

Clinical Pharmacogenetics Implementation Consortium Guideline for Thiopurine Dosing Based on TPMT and NUDT15 Genotypes: 2018 Update

Thiopurine methyltransferase (TPMT) activity exhibits a monogenic codominant inheritance and catabolizes thiopurines. TPMT variant alleles are associated with low enzyme activity and pronounced pharmacologic effects of thiopurines. Loss-of-function alleles in the NUDT15 gene are common in Asians and Hispanics and reduce the degradation of active thiopurine nucleotide metabolites, also predisposing to myelosuppression. We provide recommendations for adjusting starting doses of azathioprine, mercaptopurine, and thioguanine based on TPMT and NUDT15 genotypes (updates on www.cpicpgx.org).

Pharmacogenetic testing in oncology: a Brazilian perspective

Pharmacogenetics, a major component of individualized or precision medicine, relies on human genetic diversity. The remarkable developments in sequencing technologies have revealed that the number of genetic variants modulating drug action is much higher than previously thought and that a true personalized prediction of drug response requires attention to rare mutations (minor allele frequency, MAF<1%) in addition to polymorphisms (MAF>1%) in pharmacogenes. This has major implications for the conceptual development and clinical implementation of pharmacogenetics. Drugs used in cancer treatment have been major targets of pharmacogenetics studies, encompassing both germline polymorphisms and somatic variants in the tumor genome. The present overview, however, has a narrower scope and is focused on germline cancer pharmacogenetics, more specifically, on drug/gene pairs for which pharmacogenetics-informed prescription guidelines have been published by the Clinical Pharmacogenetics Implementation Consortium and/or the Dutch Pharmacogenetic Working Group, namely, thiopurines/TPMT, fluoropyrimidines/UGT1A1, irinotecan/UGT1A1 and tamoxifen/CYP2D6. I begin by reviewing the general principles of pharmacogenetics-informed prescription, pharmacogenetics testing and the perceived barriers to the adoption of routine pharmacogenetics testing in clinical practice. Then, I highlight aspects of the pharmacogenetics testing of the selected drug-gene pairs and finally present pharmacogenetics data from Brazilian studies pertinent to these drug-gene pairs. I conclude with the notion that pharmacogenetics testing has the potential to greatly benefit patients by enabling precision medicine applied to drug therapy, ensuring better efficacy and reducing the risk of adverse effects.

Population Diversity in Pharmacogenetics: A Latin American Perspective

Pharmacogenetics/pharmacogenomics (PGx) relies on human genetic diversity. In this review we initially examine the PGx implications of human demographic history and genetic diversity, and highlight results from recent studies on the worldwide distribution of common and rare variants in pharmacogenes. The abundance of rare variants implies that a substantial effort will be required to identify their putative functional effects and to develop reliable algorithms for PGx-guided prescription. Furthermore, variants in all pharmacogenes relevant to a drug treatment must be considered. This implies a shift of the current paradigm of PGx-informed prescription based on genotyping a few common variants in selected genes toward comprehensive sequencing approaches. The following sections deal with the impact of population admixture on PGx diversity focusing on Latin America, where a kaleidoscopic combination of individual proportions of Native American, European, and sub-Saharan African ancestries prevails. We illustrate this diversity by contrasting Brazil and Mexico, the two most populous countries in Latin America, and show that population average admixture proportions are not predictive of the corresponding proportions at the individual level. As a consequence of admixture, the genetic differentiation of common pharmacogenetic variants in Latin Americans is much attenuated in comparison to their most relevant ancestral populations. Finally, we review data for tacrolimus and warfarin to illustrate the opportunities and challenges presented by Latin American populations for PGx studies and clinical implementation.

Influence of pharmacogenetic polymorphisms and demographic variables on metformin pharmacokinetics in an admixed Brazilian cohort

AIMS

To identify pharmacogenetic and demographic variables that influence the systemic exposure to metformin in an admixed Brazilian cohort.

METHODS

The extreme discordant phenotype was used to select 106 data sets from nine metformin bioequivalence trials, comprising 256 healthy adults. Eleven single-nucleotide polymorphisms in SLC22A1, SLC22A2, SLC47A1 SLC47A2 and in transcription factor SP1 were genotyped and a validated panel of ancestry informative markers was used to estimate the individual proportions of biogeographical ancestry. Two-step (univariate followed by multivariate) regression modelling was developed to identify covariates associated with systemic exposure to metformin, accessed by the area under the plasma concentration-time curve, between 0 and 48 h (AUC_0-48h ), after single oral doses of metformin (500 or 1000 mg).

RESULTS

The individual proportions of African, Amerindian and European ancestry varied widely, as anticipated from the structure of the Brazilian population The dose-adjusted, log-transformed AUC_0-48h 's (ng h ml^-1  mg^-1 ) differed largely in the two groups at the opposite ends of the distribution histogram, namely 0.82, 0.79-0.85 and 1.08, 1.06-1.11 (mean, 95% confidence interval; P = 6.10^-26 , t test). Multivariate modelling revealed that metformin AUC_0-48h increased with age, food and carriage of rs12208357 in SLC22A1 but was inversely associated with body surface area and individual proportions of African ancestry.

CONCLUSIONS

A pharmacogenetic marker in OCT1 (SLC22A1 rs12208357), combined with demographic covariates (age, body surface area and individual proportion of African ancestry) and a food effect explained 29.7% of the variability in metformin AUC_0-48h .

CYP2D6 activity and the risk of recurrence of Plasmodium vivax malaria in the Brazilian Amazon: a prospective cohort study

Background

CYP2D6 pathway mediates the activation of primaquine into active metabolite(s) in hepatocytes. CYP2D6 is highly polymorphic, encoding CYP2D6 isoforms with normal, reduced, null or increased activity. It is hypothesized that Plasmodium vivax malaria patients with defective CYP2D6 function would be at increased risk for primaquine failure to prevent recurrence. The aim of this study was to investigate the association of CYP2D6 polymorphisms and inferred CYP2D6 phenotypes with malaria recurrence in patients from the Western Brazilian Amazon, following chloroquine/primaquine combined therapy.

Methods

The prospective cohort consisted of P. vivax malaria patients who were followed for 6 months after completion of the chloroquine/primaquine therapy. Recurrence was defined as one or more malaria episodes, 28–180 days after the initial episode. Genotyping for nine CYP2D6 SNPs and copy number variation was performed using TaqMan assays in a Fast 7500 Real-Time System. CYP2D6 star alleles (haplotypes), diplotypes and CYP2D6 phenotypes were inferred, and the activity score system was used to define the functionality of the CYP2D6 diplotypes. CYP2D6 activity scores (AS) were dichotomized at ≤ 1 (gPM, gIM and gNM-S phenotypes) and ≥ 1.5 (gNM-F and gUM phenotypes).

Results

Genotyping was successfully performed in 190 patients (44 with recurrence and 146 without recurrences). Recurrence incidence was higher in individuals presenting reduced activity CYP2D6 phenotypes (adjusted relative risk = 1.89, 95% CI 1.01–3.70; p = 0.049). Attributable risk and population attributable fraction were 11.5 and 9.9%, respectively. The time elapsed from the first P. vivax malaria episode until the recurrence did not differ between patients with AS of ≤ 1 versus ≥ 1.5 (p = 0.917).

Conclusions

The results suggest that CYP2D6 polymorphisms are associated with increased risk of recurrence of vivax malaria, following chloroquine–primaquine combined therapy. This association is interpreted as the result of reduced conversion of primaquine into its active metabolites in patients with reduced CYP2D6 enzymatic activity.

Abstract A04: Structural and functional caracterization of BARD1/CDK13 interaction

DNA damage response (DDR) is crucial for cell survivor through genomic integrity maintenance. Dysfunctions in proteins implicated in DDR prompt cell to diseases, such as cancer. BRCA1 is a central protein involved in DDR and, once damaged, plays a pivotal role in cancer development. BRCA1 structure comprises two major domains, a RING finger domain and two tandem BRCT domains (tBRCT). This structure is shared by BRCA1 major interaction partner: BARD1. The tBRCT domain is commonly found in DDR-associated proteins. To study interactions in DDR, our group charted the tBRCT interactome for 7 different proteins, using three different strategies: tandem affinity purification followed by mass spectrometry (TAP-MS), yeast two hybrid (Y2H) screening, and literature curation. CDK13 was identified as a putative BARD1 tBRCT interactor by both TAP-MS and Y2H strategies. CDK13 is described as involved in several processes, such as chromatin remodeling, transcription regulation, and splicing coordination; however, its biological role remains unclear. To confirm BARD1/CDK13 interaction, whole cell extracts obtained from MCF-7 cells were used in coimmunoprecipitation assays using anti-CDK13 and anti-BARD1. Constructs encoding the HA-tagged regions 1-706aa (N-terminal), 706-982aa (kinase domain), and 1006-1452aa (C-terminal) of CDK13 were ectopically coexpressed with Flag-tagged full-length BARD1 or GFP-tagged BARD1 tBRCT in HEK293Tcells; whole cellular extracts were used in coimmunoprecipitation assays. BARD1 and CDK13 interaction was confirmed both in ectopic and constitutive expression approaches. CDK13 silencing was performed using lentiviral constructs encoding five different short hairpin RNAs (shRNA); MCF-7 cells were transduced and submitted to puromycin selection; extracts were used to evaluate CDK13 silencing by immunoblotting. Our data suggest that BARD1 tBRCT and the CDK13 N-terminal region, as well as the kinase domain, are critical for the formation of BARD1/CDK13 complex. MCF-7 cells lacking CDK13 expression will be used to investigate the kinase role during DDR.

Citation Format: Vanessa C. Fernandes, Thales C. Nepomuceno, Renato S. Carvalho, Guilherme Suarez-Kurtz, Alvaro N. Monteiro, Marcelo A. Carvalho. Structural and functional caracterization of BARD1/CDK13 interaction [abstract]. In: Proceedings of the AACR International Conference held in cooperation with the Latin American Cooperative Oncology Group (LACOG) on Translational Cancer Medicine; May 4-6, 2017; São Paulo, Brazil. Philadelphia (PA): AACR; Clin Cancer Res 2018;24(1_Suppl):Abstract nr A04.

The Role of PALB2 in the DNA Damage Response and Cancer Predisposition

The deoxyribonucleic acid (DNA) damage response (DDR) is a major feature in the maintenance of genome integrity and in the suppression of tumorigenesis. PALB2 (Partner and Localizer of Breast Cancer 2 (BRCA2)) plays an important role in maintaining genome integrity through its role in the Fanconi anemia (FA) and homologous recombination (HR) DNA repair pathways. Since its identification as a BRCA2 interacting partner, PALB2 has emerged as a pivotal tumor suppressor protein associated to hereditary cancer susceptibility to breast and pancreatic cancers. In this review, we discuss how other DDR proteins (such as the kinases Ataxia Telangiectasia Mutated (ATM) and ATM- and Rad3-Related (ATR), mediators BRCA1 (Breast Cancer 1)/BRCA2 and effectors RAD51/DNA Polymerase η (Polη) interact with PALB2 to orchestrate DNA repair. We also examine the involvement of PALB2 mutations in the predisposition to cancer and the role of PALB2 in stimulating error-free DNA repair through the FA/HR pathway.

DNA repair genes PAXIP1 and TP53BP1 expression is associated with breast cancer prognosis

Despite remarkable advances in diagnosis, prognosis and treatment, advanced or recurrent breast tumors have limited therapeutic approaches. Many treatment strategies try to explore the limitations of DNA damage response (DDR) in tumor cells to selectively eliminate them. BRCT (BRCA1 C-terminal) domains are present in a superfamily of proteins involved in cell cycle checkpoints and the DDR. Tandem BRCT domains (tBRCT) represent a distinct class of these domains. We investigated the expression profile of 7 tBRCT genes (BARD1, BRCA1, LIG4, ECT2, MDC1, PAXIP1/PTIP and TP53BP1) in breast cancer specimens and observed a high correlation between PAXIP1 and TP53BP1 gene expression in tumor samples. Tumors with worse prognosis (tumor grade 3 and triple negative) showed reduced expression of tBRCT genes, notably, PAXIP1 and TP53BP1. Survival analyses data indicated that tumor status of both genes may impact prognosis. PAXIP1 and 53BP1 protein levels followed gene expression results, i.e., are intrinsically correlated, and also reduced in more advanced tumors. Evaluation of both genes in triple negative breast tumor samples which were characterized for their BRCA1 status showed that PAXIP1 is overexpressed in BRCA1 mutant tumors. Taken together our findings indicate that PAXIP1 status correlates with breast cancer staging, in a manner similar to what has been characterized for TP53BP1.

BRCA1 recruitment to damaged DNA sites is dependent on CDK9

Double strand break lesions, the most toxic type of DNA damage, are repaired primarily through 2 distinct pathways: homology-directed recombination (HR) and non-homologous end-joining (NHEJ). BRCA1 and 53BP1, 2 proteins containing the BRCT modular domain, play an important role in DNA damage response (DDR) by orchestrating the decision between HR and NHEJ, but the precise mechanisms regarding both pathways are not entirely understood. Previously, our group identified a putative interaction between BRCA1 and BARD1 (BRCA1-associated RING domain 1) and the cyclin-dependent kinase (CDK9). CDK9 is a component of the positive transcription elongation complex and has been implicated in genome integrity maintenance associated with the replication stress response. Here we show that CDK9 interacts with endogenous BRCA1 and BARD1 mediated by their RING finger and BRCT domains, and describe CDK9 ionizing radiation-induced foci (IRIF) formation and its co-localization with BRCA1 in DNA damage sites. Cells lacking CDK9 are characterized by an altered γ-H2AX foci dynamics after DNA damage, a reduced efficiency in HR but not in NHEJ repair, failure to form BRCA1 and RAD51 IRIF and increased sensitivity to genotoxic agents. These data indicate that CDK9 is a player in the DDR and is consistent with its participation in HR pathway by modulating BRCA1 response.

Limited sampling strategy for determining metformin area under the plasma concentration-time curve

AIM

The aim was to develop and validate limited sampling strategy (LSS) models to predict the area under the plasma concentration-time curve (AUC) for metformin.

METHODS

Metformin plasma concentrations (n = 627) at 0-24 h after a single 500 mg dose were used for LSS development, based on all subsets linear regression analysis. The LSS-derived AUC(0,24 h) was compared with the parameter 'best estimate' obtained by non-compartmental analysis using all plasma concentration data points. Correlation between the LSS-derived and the best estimated AUC(0,24 h) (r(2) ), bias and precision of the LSS estimates were quantified. The LSS models were validated in independent cohorts.

RESULTS

A two-point (3 h and 10 h) regression equation with no intercept estimated accurately the individual AUC(0,24 h) in the development cohort: r(2)  = 0.927, bias (mean, 95% CI) -0.5, -2.7-1.8% and precision 6.3, 4.9-7.7%. The accuracy of the two point LSS model was verified in study cohorts of individuals receiving single 500 or 1000 mg (r(2)  = -0.933-0.934) or seven 1000 mg daily doses (r(2)  = 0.918), as well as using data from 16 published studies covering a wide range of metformin doses, demographics, clinical and experimental conditions (r(2)  = 0.976). The LSS model reproduced previously reported results for effects of polymorphisms in OCT2 and MATE1 genes on AUC(0,24 h) and renal clearance of metformin.

CONCLUSIONS

The two point LSS algorithm may be used to assess the systemic exposure to metformin under diverse conditions, with reduced costs of sampling and analysis, and saving time for both subjects and investigators.

A phase I study of mTOR inhibitor everolimus in association with cisplatin and radiotherapy for the treatment of locally advanced cervix cancer: PHOENIX I

BACKGROUND

Cervix cancer (CC) represents the fourth most common cancer in women. Treatment involving cisplatin and radiotherapy has been the standard for locally advanced disease. Everolimus inhibits the aberrant activity of mTOR that is part of carcinogenesis in CC. Further everolimus inactivates the HPV E7 oncoprotein and inhibits its proliferation. Preclinical models have suggested that everolimus sensitizes tumoral cells and vasculature to cisplatin and radiotherapy.

METHODS

In a 3 + 3 design, the trial aimed to treat three dose levels of at least three patients with daily doses of everolimus (2.5, 5 and 10 mg/day), cisplatin and radiotherapy delivered in a 9-week interval in CC patients, stage IIB, IIIA or IIIB. Patients received everolimus from day -7 up to the last day of brachytherapy. Primary objective was to evaluate safety, toxicity and the maximum-tolerated dose (MTD) of everolimus in association with cisplatin and radiotherapy. Pharmacokinetic (PK) parameters and response rates were analyzed as secondary objectives.

RESULTS

Thirteen patients were enrolled, 6 at 2.5 mg, 3 at 5 mg and 4 at 10 mg. Four patients did not complete the planned schedule, 1 at 2.5 mg presented grade 4 acute renal failure interpreted as dose-limiting toxicity (DLT) and 3 at 10 mg: 1 with disease progression, and 2 with DLTs-1 grade 3 rash and 1 grade 4 neutropenia. PK results were characterized by dose-dependent increases in AUC and C max.

CONCLUSIONS

The MTD of everolimus in combination with cisplatin and radiotherapy has been defined as 5 mg/day. The data regarding safety and response rates support further studies.

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.

Impact of visceral leishmaniasis and curative chemotherapy on cytochrome P450 activity in Brazilian patients

AIMS

The aim of the present study was to investigate the impact of human visceral leishmaniasis (VL) and curative chemotherapy on the activity of cytochrome P450 (CYP) 3A, CYP2C9 and CYP2C19 in patients from an endemic region in Brazil.

METHODS

Adult patients with parasitologically confirmed VL were given a CYP phenotyping cocktail, comprising midazolam, omeprazole and losartan, immediately before (Study phase 1), 2-3 days (phase 2) and 3-6 months (phase 3) after curative VL chemotherapy. CYP activity was assessed by the apparent clearance of midazolam (CYP3A), omeprazole/5-hydroxyomeprazol ratio in plasma (CYP2C19) and losartan/E3174 ratio in urine (CYP2C9).

RESULTS

Mean values (95% confidence interval) in phases 1, 2 and 3 were, respectively: log apparent midazolam clearance, 1.21 (1.10-1.31), 1.45 (1.32-1.57) and 1.35 (1.26-1.44) ml min(-1)  kg(-1) ; omeprazole/5-hydroxyomeprazole ratio, 0.78 (0.61-0.94), 0.45 (0.27-0.63) and 0.37 (0.20-0.55); losartan/E3174 ratio, 0.66 (0.39-0.92), 0.35 (0.20-0.50) and 0.35 (0.16-0.53). Analysis of variance revealed significant differences in CYP3A (P = 0.018) and CYP2C19 (P = 0.008), but not CYP2C9 (P = 0.11) phenotypic activity, across the three study phases.

CONCLUSION

The phenotypic activities of CYP3A4 and CYP2C19 were significantly reduced during acute VL compared with post-chemotherapy. We propose that increased plasma concentrations of proinflammatory cytokines during active disease account for the suppression of CYP activity. The failure to detect significant changes in CYP2C9 activity in the overall cohort may reflect differential effects of the inflammatory process on the expression of CYP isoforms, although the possibility of insufficient statistical power cannot be dismissed.

Effect of UGT1A1, UGT1A3, DIO1 and DIO2 polymorphisms on L-thyroxine doses required for TSH suppression in patients with differentiated thyroid cancer

AIM

To evaluate the impact of genetic polymorphisms in uridine 5'-glucuronosylytansferases UGT1A1 and UGT1A3 and iodothyronine-deiodinases types 1 and 2 on levothyroxine (T4 ; 3,5,3',5'-triiodo-L-thyronine) dose requirement for suppression of thyrotropin (TSH) secretion in patients with differentiated thyroid cancer (DTC).

METHODS

Patients (n = 268) submitted to total thyroidectomy and ablation by (131) I, under T4 therapy for at least 6 months were recruited in three public institutions in Brazil. Multivariate regression modelling was applied to assess the association of T4 dosing with polymorphisms in UGT1A1 (rs8175347), UGT1A3 (rs3806596 and rs1983023), DIO1 (rs11206244 and rs2235544) and DIO2 (rs225014 and rs12885300), demographic and clinical variables.

RESULTS

A regression model including UGT1A haplotypes, age, gender, body weight and serum TSH concentration accounted for 39% of the inter-individual variation in the T4 dosage. The association of T4 dose with UGT1A haplotype is attributed to reduced UGT1A1 expression and T4 glucuronidation in liver of carriers of low expression UGT1A1 rs8175347 alleles. The DIO1 and DIO2 genotypes had no influence of T4 dosage.

CONCLUSION

UGT1A haplotypes associate with T4 dosage in DTC patients, but the effect accounts for only 2% of the total variability and recommendation of pre-emptive UGT1A genotyping is not warranted.

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.

Distribution of CYP2D6 alleles and phenotypes in the Brazilian population

The CYP2D6 enzyme is one of the most important members of the cytochrome P450 superfamily. This enzyme metabolizes approximately 25% of currently prescribed medications. The CYP2D6 gene presents a high allele heterogeneity that determines great inter-individual variation. The aim of this study was to evaluate the variability of CYP2D6 alleles, genotypes and predicted phenotypes in Brazilians. Eleven single nucleotide polymorphisms and CYP2D6 duplications/multiplications were genotyped by TaqMan assays in 1020 individuals from North, Northeast, South, and Southeast Brazil. Eighteen CYP2D6 alleles were identified in the Brazilian population. The CYP2D6*1 and CYP2D6*2 alleles were the most frequent and widely distributed in different geographical regions of Brazil. The highest number of CYPD6 alleles observed was six and the frequency of individuals with more than two copies ranged from 6.3% (in Southern Brazil) to 10.2% (Northern Brazil). The analysis of molecular variance showed that CYP2D6 is homogeneously distributed across different Brazilian regions and most of the differences can be attributed to inter-individual differences. The most frequent predicted metabolic status was EM (83.5%). Overall 2.5% and 3.7% of Brazilians were PMs and UMs respectively. Genomic ancestry proportions differ only in the prevalence of intermediate metabolizers. The IM predicted phenotype is associated with a higher proportion of African ancestry and a lower proportion of European ancestry in Brazilians. PM and UM classes did not vary among regions and/or ancestry proportions therefore unique CYP2D6 testing guidelines for Brazilians are possible and could potentially avoid ineffective or adverse events outcomes due to drug prescriptions.

Probing structure-function relationships in missense variants in the carboxy-terminal region of BRCA1

Germline inactivating variants in BRCA1 lead to a significantly increased risk of breast and ovarian cancers in carriers. While the functional effect of many variants can be inferred from the DNA sequence, determining the effect of missense variants present a significant challenge. A series of biochemical and cell biological assays have been successfully used to explore the impact of these variants on the function of BRCA1, which contribute to assessing their likelihood of pathogenicity. It has been determined that variants that co-localize with structural or functional motifs are more likely to disrupt the stability and function of BRCA1. Here we assess the functional impact of 37 variants chosen to probe the functional impact of variants in phosphorylation sites and in the BRCT domains. In addition, we perform a meta-analysis of 170 unique variants tested by the transcription activation assays in the carboxy-terminal domain of BRCA1 using a recently developed computation model to provide assessment for functional impact and their likelihood of pathogenicity.

Characterization of LGALS3 (galectin-3) as a player in DNA damage response

DNA damage repair (DDR) is an orchestrated process encompassing the injury detection to its complete resolution. DNA double-strand break lesions are repaired mainly by two distinct mechanisms: the error-free homologous recombination (HR) and the error-prone non-homologous end-joining. Galectin-3 (GAL3) is the unique member of the chimeric galectins subfamily and is reported to be involved in several cancer development and progression related events. Recently our group described a putative protein interaction between GAL3 and BARD1, the main partner of breast and ovarian cancer susceptibility gene product BRCA1, both involved in HR pathway. In this report we characterized GAL3/BARD1 protein interaction and evaluated the role of GAL3 in DDR pathways using GAL3 silenced human cells exposed to different DNA damage agents. In the absence of GAL3 we observed a delayed DDR response activation, as well as a decrease in the G 2/M cell cycle checkpoint arrest associated with HR pathway. Moreover, using a TAP-MS approach we also determined the protein interaction network of GAL3.