A European Spectrum of Pharmacogenomic Biomarkers: Implications for Clinical Pharmacogenomics

Cost-effectiveness analysis of CYP3A5 genotype-guided tacrolimus dosing in solid organ transplantation using real-world data

The objective of this study was to estimate the cost-effectiveness of CYP3A5 genotype-guided tacrolimus dosing in kidney, liver, heart, and lung transplant recipients relative to standard of care (SOC) tacrolimus dosing, from a US healthcare payer perspective. We developed decision-tree models to compare economic and clinical outcomes between CYP3A5 genotype-guided and SOC tacrolimus therapy in the first six months post-transplant. We derived inputs for CYP3A5 phenotype frequencies and physician use of genotype test results to inform clinical care from literature; tacrolimus exposure [high vs low tacrolimus time in therapeutic range using the Rosendaal algorithm (TAC TTR-Rosendaal)] and outcomes (incidences of acute tacrolimus nephrotoxicity, acute cellular rejection, and death) from real-world data; and costs from the Medicare Fee Schedule and literature. We calculated cost per avoided event and performed sensitivity analyses to evaluate the robustness of the results to changes in inputs. Incremental costs per avoided event for CYP3A5 genotype-guided vs SOC tacrolimus dosing were $176,667 for kidney recipients, $364,000 for liver recipients, $12,982 for heart recipients, and $93,333 for lung recipients. The likelihood of CYP3A5 genotype-guided tacrolimus dosing leading to cost-savings was 19.8% in kidney, 32.3% in liver, 51.8% in heart, and 54.1% in lung transplant recipients. Physician use of genotype results to guide clinical care and the proportion of patients with a high TAC TTR-Rosendaal were key parameters driving the cost-effectiveness of CYP3A5 genotype-guided tacrolimus therapy. Relative to SOC, CYP3A5 genotype-guided tacrolimus dosing resulted in a slightly greater benefit at a higher cost. Further economic evaluations examining intermediary outcomes (e.g., dose modifications) are needed, particularly in populations with higher frequencies of CYP3A5 expressers.

Initial Low-Dose Hydroxyurea and Anagrelide Combination in Essential Thrombocythemia: Comparable Response with Lower Toxicity

Background and Objectives: Essential thrombocythemia (ET) is a myeloproliferative neoplasm that overproduces platelets and is associated with life-threatening thrombosis. Medical cytoreduction either with hydroxyurea (HU) or anagrelide (AG) is widely used, but drug intolerance or resistance are major concerns. Low-dose combination of HU and AG as an alternative strategy has been explored in various studies. It showed comparable response with acceptable toxicity in second-line settings for patients who experienced side effects from prior monotherapy. In this study, we evaluated the efficacy and safety of upfront combination for ET patients. Materials and Methods: From January 2018 to June 2022, a total of 241 ET patients with intermediate to high risk were enrolled. We identified 21 patients with initial drug combinations and compared treatment outcomes and adverse events (AEs) between combination and monotherapy groups. Results: The median age was 62 years old (range, 26-87) and median platelet count was 912 × 109/L (range, 520-1720). Overall treatment response did not exhibit significant differences between the groups, although there was a trend towards a lower response rate in patients treated with AG alone at 3 months post-treatment (AG + HU, 85.7% vs. AG alone, 75.4%, p = 0.068). AEs of any grade occurred in 52.3% of the combination group, 44.3% of the HU monotherapy group, and 43.4% of the AG single group, respectively. Of note was that the HU plus AG combination group suffered a lower incidence of grade 3-4 AEs compared to the other two groups, with statistical significance (p = 0.008 for HU monotherapy vs. combination therapy and p < 0.01 for AG monotherapy vs. combination therapy). Conclusions: Our findings demonstrated that the upfront low-dose combination approach showed feasible clinical outcomes with significantly lower severe AEs compared to conventional monotherapy. These results may offer valuable insights to clinicians for future prospective investigations.

Pharmacogenomics-assisted schizophrenia management: A hybrid type 2 effectiveness-implementation study protocol to compare the clinical utility, cost-effectiveness, and barriers

OBJECTIVES

The response to antipsychotic therapy is highly variable. Pharmacogenomic (PGx) factors play a major role in deciding the effectiveness and safety of antipsychotic drugs. A hybrid type 2 effectiveness-implementation research will be conducted to evaluate the clinical utility (safety and efficacy), cost-effectiveness, and facilitators and barriers in implementing PGx-assisted management compared to standard of care in patients with schizophrenia attending a tertiary care hospital in eastern India.

METHODS

In part 1, a randomized controlled trial will be conducted. Adult patients with schizophrenia will be randomized (2: 1) to receive PGx-assisted treatment (drug and regimen selection depending on the results of single-nucleotide polymorphisms in genes DRD2, HTR1A, HTR2C, ABCB1, CYP2D6, CYP3A5, and CYP1A2) or the standard of care. Serum drug levels will be measured. The patients will be followed up for 12 weeks. The primary endpoint is the difference in the Udvalg for Kliniske Undersøgelser Side-Effect Rating Scale score between the two arms. In part 2, the cost-effectiveness of PGx-assisted treatment will be evaluated. In part 3, the facilitators and barriers to implementing PGx-assisted treatment for schizophrenia will be explored using a qualitative design.

EXPECTED OUTCOME

The study findings will help in understanding whether PGx-assisted management has a clinical utility, whether it is cost-effective, and what are the facilitators and barriers to implementing it in the management of schizophrenia.

TRIAL REGISTRATION

The study has been registered with the Clinical Trials Registry-India (CTRI/2023/08/056210).

The Frequency of CYP2D6 and CYP3A4/5 Genotypes and The Impact of Their Allele Translation and Phenoconversion-Predicted Enzyme Activity on Risperidone Pharmacokinetics in Saudi Children with Autism

Data on the role of CYP2D6 and CYP3A4/5 polymorphisms in relation to risperidone (RIS) pharmacokinetics (PK) in children are relatively limited and inconsistent. This is partially attributable to the limited coverage of CYP2D6 and CYP3A4/5 metabolizer phenotypes, particularly those of poor and ultrarapid metabolizers (PMs and UMs), which has led to calls for studies of populations with a non-European background that may carry variants that are less frequent in Europeans. Children ≤ 18 years old with at least 8 weeks of a RIS-based regimen were recruited from three autism centers in Riyadh, Saudi Arabia. The primary outcomes measured were plasma concentrations of RIS and 9-hydroxyrisperidone (9-OH-RIS) and their dose-adjusted (C/D) ratios as a function of phenotypes and activity score (AS). For accurate DNA genotyping, targeted pharmacogenomic testing with the Axiom PharmacoFocus Array was performed via examination of a broad collection of probesets targeting CYP2D6 and CYP3A4/5 variants. The frequency of genotypes/phenotypes and the impact of their allele translation and phenoconversion-predicted enzyme activity were examined. The final cohort included 83 individuals. The most common CYP2D6 phenotype in our population was normal metabolizers (NMs, 66.3%). Inconsistent with some previous studies, the three phenotypes of intermediate metabolizers (IMs), NMs, and UMs were significantly different in terms of RIS concentration, the RIS/9-OH-RIS ratio, the RIS C/D ratio and the 9-OH-RIS C/D ratio. According to AS analyses, there were statistically significant differences in the RIS concentration (P = 0.013), RIS/9-OH-RIS ratio (P < 0.001) and RIS C/D ratio (P = 0.030) when patients were categorized into AS ≤ 1 vs. AS > 1. None of the CYP3A4/5 star allele translated phenotypes revealed a significant influence on any of the RIS PK parameters. Notably, neither CYP2D6 nor CYP3A4/5 phenotyping demonstrated a significant impact on the total active moiety, suggesting that other gene variants could modulate RIS PK. The study confirmed the previously reported partial impact of the CYP2D6 gene on RIS PK. However, future studies using contemporary genotyping techniques targeting a wide range of variants in other candidate genes must be conducted to further examine their interactive effects on RIS PK and the clinical response.

Knowledge and attitudes of medical and pharmacy students about pharmacogenomics: a systematic review and meta-analysis

Pharmacogenomics (PGx) is rapidly growing branch of molecular genetics with high potentials to influence therapeutics. This review evaluates knowledge and attitudes of medical and pharmacy students about PGx. A literature search was conducted in electronic databases and studies were selected by following precise eligibility criteria. After quality assessment, studies were reviewed systematically, and meta-analyses of proportions were performed to estimate response rates of students. Fifteen studies (5509 students; 69% [95% confidence interval (CI): 60%, 77%] females) were included. Among students, 28% [95%CI: 12, 46] had adequate PGx knowledge; 65% [95%CI: 55, 75] were willing to have PGx test for their own risk assessment; 78% [95%CI: 71, 84] had intention to incorporate PGx in future practice; and 32% [95%CI: 21, 43] were satisfied with current PGx component of curriculum. Age, advanced year of educational program, and more time spent in PGx education were positively associated with PGx knowledge and positive attitudes.

Can Pharmacogenetic Variants in TPMT, MTHFR and SLCO1B1 Genes Be Used as Potential Markers of Outcome Prediction in Systemic Sclerosis Patients?

Systemic sclerosis (SSc) is a rare connective tissue disorder with highest morbidity and mortality among rheumatologic diseases. Disease progression is highly heterogeneous between patients, implying a strong need for individualization of therapy. Four pharmacogenetic variants, namely TPMT rs1800460, TPMT rs1142345, MTHFR rs1801133 and SLCO1B1 rs4149056 were tested for association with severe disease outcomes in 102 patients with SSc from Serbia treated either with immunosuppressants azathioprine (AZA) and methotrexate (MTX) or with other types of medications. Genotyping was performed using PCR-RFLP and direct Sanger sequencing. R software was used for statistical analysis and development of polygenic risk score (PRS) model. Association was found between MTHFR rs1801133 and higher risk for elevated systolic pressure in all patients except those prescribed with MTX, and higher risk for kidney insufficiency in patients prescribed with other types of drugs. In patients treated with MTX, variant SLCO1B1 rs4149056 was protective against kidney insufficiency. For patients receiving MTX a trend was shown for having a higher PRS rank and elevated systolic pressure. Our results open a door wide for more extensive research on pharmacogenomics markers in patients with SSc. Altogether, pharmacogenomics markers could predict the outcome of patients with SSc and help in prevention of adverse drug reactions.

Integrative toxicogenomics: Advancing precision medicine and toxicology through artificial intelligence and OMICs technology

More information about a person's genetic makeup, drug response, multi-omics response, and genomic response is now available leading to a gradual shift towards personalized treatment. Additionally, the promotion of non-animal testing has fueled the computational toxicogenomics as a pivotal part of the next-gen risk assessment paradigm. Artificial Intelligence (AI) has the potential to provid new ways analyzing the patient data and making predictions about treatment outcomes or toxicity. As personalized medicine and toxicogenomics involve huge data processing, AI can expedite this process by providing powerful data processing, analysis, and interpretation algorithms. AI can process and integrate a multitude of data including genome data, patient records, clinical data and identify patterns to derive predictive models anticipating clinical outcomes and assessing the risk of any personalized medicine approaches. In this article, we have studied the current trends and future perspectives in personalized medicine & toxicology, the role of toxicogenomics in connecting the two fields, and the impact of AI on personalized medicine & toxicology. In this work, we also study the key challenges and limitations in personalized medicine, toxicogenomics, and AI in order to fully realize their potential.

Global distribution of functionally important CYP2C9 alleles and their inferred metabolic consequences

BACKGROUND

Genetic variability in the cytochrome P450 CYP2C9 constitutes an important predictor for efficacy and safety of various commonly prescribed drugs, including coumarin anticoagulants, phenytoin and multiple non-steroidal anti-inflammatory drugs (NSAIDs). A global map of CYP2C9 variability and its inferred functional consequences has been lacking.

RESULTS

Frequencies of eight functionally relevant CYP2C9 alleles (*2, *3, *5, *6, *8, *11, *13 and *14) were analyzed. In total, 108 original articles were identified that included genotype data from a total of 81,662 unrelated individuals across 70 countries and 40 unique ethnic groups. The results revealed that CYP2C9*2 was most abundant in Europe and the Middle East, whereas CYP2C9*3 was the main reason for reduced CYP2C9 activity across South Asia. Our data show extensive variation within superpopulations with up to tenfold differences between geographically adjacent populations in Malaysia, Thailand and Vietnam. Translation of genetic CYP2C9 variability into functional consequences indicates that up to 40% of patients in Southern Europe and the Middle East might benefit from warfarin and phenytoin dose reductions, while 3% of patients in Southern Europe and Israel are recommended to reduce starting doses of NSAIDs.

CONCLUSIONS

This study provides a comprehensive map of the genetic and functional variability of CYP2C9 with high ethnogeographic resolution. The presented data can serve as a useful resource for CYP2C9 allele and phenotype frequencies and might guide the optimization of genotyping strategies, particularly for indigenous and founder populations with distinct genetic profiles.

The genetic landscape of major drug metabolizing cytochrome P450 genes-an updated analysis of population-scale sequencing data

Genes encoding cytochrome P450 enzymes (CYPs) are extremely polymorphic and multiple CYP variants constitute clinically relevant biomarkers for the guidance of drug selection and dosing. We previously reported the distribution of the most relevant CYP alleles using population-scale sequencing data. Here, we update these findings by making use of the increasing wealth of data, incorporating whole exome and whole genome sequencing data from 141,614 unrelated individuals across 12 human populations. We furthermore extend our previous studies by systematically considering also uncharacterized rare alleles and reveal that they contribute between 1.5% and 17.5% to the overall genetically encoded functional variability. By using established guidelines, we aggregate and translate the available sequencing data into population-specific patterns of metabolizer phenotypes. Combined, the presented data refine the worldwide landscape of ethnogeographic variability in CYP genes and aspire to provide a relevant resource for the optimization of population-specific genotyping strategies and precision public health.

Considerations into pharmacogenomics of COVID-19 pharmacotherapy: Hope, hype and reality

COVID-19 medicines, such as molnupiravir are beginning to emerge for public health and clinical practice. On the other hand, drugs display marked variability in their efficacy and safety. Hence, COVID-19 medicines, as with all drugs, will be subject to the age-old maxim “one size prescription does not fit all”. In this context, pharmacogenomics is the study of genome-by-drug interactions and offers insights on mechanisms of patient-to-patient and between-population variations in drug efficacy and safety. Pharmacogenomics information is crucial to tailoring the patients' prescriptions to achieve COVID-19 preventive and therapeutic interventions that take into account the host biology, patients’ genome, and variable environmental exposures that collectively influence drug efficacy and safety. This expert review critically evaluates and summarizes the pharmacogenomics and personalized medicine aspects of the emerging COVID-19 drugs, and other selected drug interventions deployed to date. Here, we aim to sort out the hope, hype, and reality and suggest that there are veritable prospects to advance COVID-19 medicines for public health benefits, provided that pharmacogenomics is considered and implemented adequately. Pharmacogenomics is an integral part of rational and evidence-based medical practice. Scientists, health care professionals, pharmacists, pharmacovigilance practitioners, and importantly, patients stand to benefit by expanding the current pandemic response toolbox by the science of pharmacogenomics, and its applications in COVID-19 medicines and clinical trials.

Pharmacogenes that demonstrate high association evidence according to CPIC, DPWG, and PharmGKB

BACKGROUND

Different levels of evidence related to the variable responses of individuals to drug treatment have been reported in various pharmacogenomic (PGx) databases. Identification of gene-drug pairs with strong association evidence can be helpful in prioritizing the implementation of PGx guidelines and focusing on a gene panel. This study aimed to determine the pharmacogenes with the highest evidence-based association and to indicate their involvement in drug-gene interactions.

METHODOLOGY

The publicly available datasets CPIC, DPWG, and PharmGKB were selected to determine the pharmacogenes with the highest drug outcome associations. The upper two levels of evidence rated by the three scoring methods were specified (levels A-B in CPIC, 3-4 in DPWG, or 1-2 levels in PharmGKB). The identified pharmacogenes were further ranked in this study based on the number of medications they interacted with.

RESULTS

Fifty pharmacogenes, with high to moderately high evidence of associations with drug response alterations, with potential influence on the therapeutic and/or toxicity outcomes of 152 drugs were identified. CYP2D6, CYP2C9, CYP2C19, G6PD, HLA-B, SLCO1B1, CACNA1S, RYR1, MT-RNR1, and IFNL4 are the top 10 pharmacogenes, where each is predicted to impact patients' responses to ≥5 drugs.

CONCLUSION

This study identified the most important pharmacogenes based on the highest-ranked association evidence and their frequency of involvement in affecting multiple drugs. The obtained data is useful for customizing a gene panel for PGx testing. Identifying the strength of scientific evidence supporting drug-gene interactions aids drug prescribers in making the best clinical decision.

PharmaKoVariome database for supporting genetic testing

Pharmacogenomics (PGx) provides information about routine precision medicine, based on the patient's genotype. However, many of the available information about human allele frequencies, and about clinical drug-gene interactions, is based on American and European populations. PharmaKoVariome database was constructed to support genetic testing for safe prescription and drug development. It consolidated and stored 2507 diseases, 11 459 drugs and 61 627 drug-target or druggable genes from public databases. PharmaKoVariome precomputed ethnic-specific abundant variants for approximately 120 M single-nucleotide variants of drug-target or druggable genes. A user can search by gene symbol, drug name, disease and reference SNP ID number (rsID) to statistically analyse the frequency of ethnical variations, such as odds ratio and P-values for related genes. In an example study, we observed five Korean-enriched variants in the CYP2B6 and CYP2D6 genes, one of which (rs1065852) is known to be incapable of metabolizing drug. It is also shown that 4-6% of North and East Asians have risk factors for drugs metabolized by the CYP2D6 gene. Therefore, PharmaKoVariome is a useful database for pharmaceutical or diagnostic companies for developing diagnostic technologies that can be applied in the Asian PGx industry. Database URL: http://www.pharmakovariome.com/.

Drug genetic associations with COVID-19 manifestations: a data mining and network biology approach

Available drugs have been used as an urgent attempt through clinical trials to minimize severe cases of hospitalizations with Coronavirus disease (COVID-19), however, there are limited data on common pharmacogenomics affecting concomitant medications response in patients with comorbidities. To identify the genomic determinants that influence COVID-19 susceptibility, we use a computational, statistical, and network biology approach to analyze relationships of ineffective concomitant medication with an adverse effect on patients. We statistically construct a pharmacogenetic/biomarker network with significant drug-gene interactions originating from gene-disease associations. Investigation of the predicted pharmacogenes encompassing the gene-disease-gene pharmacogenomics (PGx) network suggests that these genes could play a significant role in COVID-19 clinical manifestation due to their association with autoimmune, metabolic, neurological, cardiovascular, and degenerative disorders, some of which have been reported to be crucial comorbidities in a COVID-19 patient.

CYP2C19 and CYP2D6 Genotypes and Metabolizer Status Distribution in a Bulgarian Psychiatric Cohort

CYP2D6 and CYP2C19 are enzymes of essential significance for the pharmacokinetics of a multitude of commonly used antidepressants, antipsychotics, antiemetics, β-blockers, opioids, antiestrogen, antacids, etc. Polymorphisms in the respective genes are well established as resulting in functional differences, which in turn can impact safety and efficacy. Importantly, the prevalence of genetic CYP2D6 and CYP2C19 variability differs drastically between populations. Drawing on the limited information concerning genotype frequencies in Bulgaria, we here analyzed 742 Bulgarian psychiatric patients predominantly diagnosed with depression and/or anxiety. Specifically, we analyzed frequencies of CYPC19*2, *4 and *17, as well as of CYP2D6*2, *3, *4, *5, *6, *10 and *41. In total, 571 out of 742 patients (77%) carried at least one variant which impacts metabolizer status. Overall, 48.6% of the studied individuals were classified as non-normal metabolizers of CYP2D6 with most exhibiting reduced function (38.2% intermediate metabolizers and 6.6% poor metabolizers). In contrast, for CYP2C19, the majority of non-normal metabolizers showed increased functionality (28.9% rapid and 5.5% ultrarapid metabolizers), while reduced activity metabolizer status accounted for 25.6% (23.8% intermediate and 1.8% poor metabolizers). These results provide an important resource to assess the genetically encoded functional variability of CYP2D6 and CYP2C19 which may have significant implications for precision medicine in Bulgarian psychiatry practice.

Impact of SLCO1B1*5 on Flucloxacillin and Co-Amoxiclav-Related Liver Injury

Background: Idiosyncratic drug-induced liver injury (DILI) is a serious uncommon disease that may develop as a result of the intake of certain drugs such as the antimicrobials flucloxacillin and co-amoxiclav. The reported cases showed significant associations between DILI and various human leukocyte (HLA) markers. The solute carrier organic anion transporter 1B1 (SLCO1B1), a non-HLA candidate gene, was previously reported as a risk factor for liver injury induced by rifampin and methimazole. This study presumed that SLCO1B1 may play a general role in the DILI susceptibility and therefore investigated the association of rs4149056 (SLCO1B1*5, T521C) polymorphism with flucloxacillin- and co-amoxiclav-induced liver injury. Methodology: We recruited 155 and 165 DILI cases of white ancestral origin from various European countries but mainly from the United Kingdom owing to flucloxacillin and co-amoxiclav, respectively. Only adult patients (≥18 years) who were diagnosed with liver injury and who showed i) clinical jaundice or bilirubin >2x the upper limit of normal (ULN), ii) alanine aminotransferase (ALT) >5x ULN or iii) alkaline phosphatase (ALP) >2x ULN and bilirubin > ULN were selected. The population reference sample (POPRES), a European control group (n = 282), was used in comparison with the investigated cases. TaqMan SNP genotyping custom assay designed by Applied Biosystems was used to genotype both DILI cohorts for SLCO1B1 polymorphism (rs4149056). Allelic discrimination analysis was performed using a step one real-time PCR machine. Genotype differences between cases and controls were examined using Fisher's exact test. GraphPad Prism version 5.0 was used to determine the p-value, odds ratio, and 95% confidence interval. Compliance of the control group with Hardy-Weinberg equilibrium was proven using a web-based calculator available at https://wpcalc.com/en/equilibrium-hardy-weinberg/. Results: A small number of cases failed genotyping in each cohort. Thus, only 149 flucloxacillin and 162 co-amoxiclav DILI cases were analyzed. Genotyping of both DILI cohorts did not show evidence of association with the variant rs4149056 (T521C) (OR = 0.71, 95% CI = 0.46-1.12; p = 0.17 for flucloxacillin cases and OR = 0.87, 95% CI = 0.56-1.33; p = 0.58 for co-amoxiclav), although slightly lower frequency (22.8%) of positive flucloxacillin cases was noticed than that of POPRES controls (29.4%). Conclusion: Carriage of the examined allele SLCO1B1*5 is not considered a risk factor for flucloxacillin DILI or co-amoxiclav DILI as presumed. Testing a different allele (SLCO1B1*1B) and another family member gene (SLCO1B3) may still be needed to provide a clearer role of SLCO1B drug transporters in DILI development-related to the chosen antimicrobials.

Prevalence of exposure to pharmacogenetic drugs by the Saudis treated at the health care centers of the Ministry of National Guard

Background

The drugs impacted by genetic variants are known as pharmacogenetic (PGx) drugs. Patients’ responses to these drugs may vary according to the variability in patients’ genetic makeup. Hence, exploring the pharmacogenes that affect drug treatment is vital to ensure optimal therapy and patients’ safety. This study aimed to describe the usage rate of PGx drugs and the frequency of relevant variants in the Saudi population.

Methodology

Prescription patterns over seven years (2015–2021) for Saudi patients on PGx drugs treated at the Ministry of National Guard-Health Affairs (MNG-HA) were investigated. Only registered drugs in the MNG-HA formulary (n = 78) were included. The patients were subgrouped into four age groups: ≤24, 25–44, 45–64, and ≥65 years. Further subgrouping was made according to gender and drugs’ therapeutic categories following anatomical therapeutic chemical (ATC) classification.

Furthermore, an online searching was carried out to identify the pharmacogenes reported in the literature among healthy Saudis. The search included 45 genes that may affect drug outcomes based on evidence rated by either CPIC (A-B levels) or PharmGKB (1–2 levels).

Results

The screened patients were 1,483,905. Patients on PGx drugs accounted for 46.7% (n = 693,077 patients). The analgesic group was the most prescribed drug category (47%), which included ibuprofen (20.5%), celecoxib (6.3%), tramadol (5.8%), and others. Cardiovascular agents were the second-most utilized drug class (24.4%). Omeprazole was the second most commonly used medication (11.1%) but ranked third as a class (gastroenterology). Females used PGx drugs more frequently than males (53.5% versus 46.5%) and a higher usage rate by patients aged 45–64 years (31.3%) was noted. The cytochrome P450 genes (CYP2C9, CYP2C19, and CYP2D6) were estimated to impact responses of 54.3% (n = 1,156,113) of the used drugs (27.2% are possibly affected by CYP2C9, 12.8% by CYP2C19, and 14.3% by CYP2D6). Thirty-five pharmacogenes that characterize Saudi population and their variants’ allele frequencies were identified from previous reports. This study presents the largest reported number of genes that may affect drug therapies among Saudis.

Conclusion

This study confirmed that a high percentage of Saudi patients use PGx drugs and various genotypes of certain pharmacogenes are inherited by the Saudi population.

Population pharmacogenomics: an update on ethnogeographic differences and opportunities for precision public health

Both safety and efficacy of medical treatment can vary depending on the ethnogeographic background of the patient. One of the reasons underlying this variability is differences in pharmacogenetic polymorphisms in genes involved in drug disposition, as well as in drug targets. Knowledge and appreciation of these differences is thus essential to optimize population-stratified care. Here, we provide an extensive updated analysis of population pharmacogenomics in ten pharmacokinetic genes (CYP2D6, CYP2C19, DPYD, TPMT, NUDT15 and SLC22A1), drug targets (CFTR) and genes involved in drug hypersensitivity (HLA-A, HLA-B) or drug-induced acute hemolytic anemia (G6PD). Combined, polymorphisms in the analyzed genes affect the pharmacology, efficacy or safety of 141 different drugs and therapeutic regimens. The data reveal pronounced differences in the genetic landscape, complexity and variant frequencies between ethnogeographic groups. Reduced function alleles of CYP2D6, SLC22A1 and CFTR were most prevalent in individuals of European descent, whereas DPYD and TPMT deficiencies were most common in Sub-Saharan Africa. Oceanian populations showed the highest frequencies of CYP2C19 loss-of-function alleles while their inferred CYP2D6 activity was among the highest worldwide. Frequencies of HLA-B*15:02 and HLA-B*58:01 were highest across Asia, which has important implications for the risk of severe cutaneous adverse reactions upon treatment with carbamazepine and allopurinol. G6PD deficiencies were most frequent in Africa, the Middle East and Southeast Asia with pronounced differences in variant composition. These variability data provide an important resource to inform cost-effectiveness modeling and guide population-specific genotyping strategies with the goal of optimizing the implementation of precision public health.

Pharmacogenetic distinction of the Croatian population from the European average

AIM

To compare the Croatian and European population in terms of allele frequencies of clinically relevant polymorphisms in drug absorption, distribution, metabolism, and excretion (ADME) genes.

METHODS

In 429 Croatian participants, we genotyped 27 loci in 20 ADME genes. The obtained frequencies were merged with the published frequencies for the Croatian population by sample size weighting. The study sample obtained in this way was compared with the average data for the European population from the gnomAD database.

RESULTS

Variant allele frequencies in the Croatian population were higher in three and lower in two polymorphisms (Benjamini-Hochberg-corrected P values: 0.0027 for CYP2B6*4 rs2279343, CYP2C9*2 rs1799853, and VKORC1 rs9923231; 0.0297 for GSTP1 rs1695; 0.0455 for CYP2A6 rs1801272) compared with the European population. The most marked difference was observed for CYP2B6*4 (9.3% in Europe vs 24.3% in Croatia). The most clinically relevant findings were higher variant allele frequencies in two polymorphisms related to lower warfarin requirements: VKORC1*2 (34.9% in Europe vs 40.1% in Croatia) and CYP2C9*2 (12.3% in Europe vs 14.7% in Croatia). This indicates that three-quarters of Croatian people have at least one variant allele at these loci. Variants in genes GSTP1 and CYP2A6 were significantly less frequently observed in Croatia.

CONCLUSIONS

Croatian population has a higher bleeding and over-anticoagulation risk, which is why we recommend the prescription of lower doses of anticoagulation drugs such as warfarin and acenocoumarol. Lower phenytoin, and higher bupropion and efavirenz doses are also recommended in the Croatian population.

Ethnic Diversity and Warfarin Pharmacogenomics

Warfarin has remained the most commonly prescribed vitamin K oral anticoagulant worldwide since its approval in 1954. Dosing challenges including having a narrow therapeutic window and a wide interpatient variability in dosing requirements have contributed to making it the most studied drug in terms of genotype-phenotype relationships. However, most of these studies have been conducted in Whites or Asians which means the current pharmacogenomics evidence-base does not reflect ethnic diversity. Due to differences in minor allele frequencies of key genetic variants, studies conducted in Whites/Asians may not be applicable to underrepresented populations such as Blacks, Hispanics/Latinos, American Indians/Alaska Natives and Native Hawaiians/other Pacific Islanders. This may exacerbate health inequalities when Whites/Asians have better anticoagulation profiles due to the existence of validated pharmacogenomic dosing algorithms which fail to perform similarly in the underrepresented populations. To examine the extent to which individual races/ethnicities are represented in the existing body of pharmacogenomic evidence, we review evidence pertaining to published pharmacogenomic dosing algorithms, including clinical utility studies, cost-effectiveness studies and clinical implementation guidelines that have been published in the warfarin field.

Pharmacogenetics of Drugs Used in the Treatment of Cancers

Pharmacogenomics is based on the understanding of the individual differences in drug use, the response to drug therapy (efficacy and toxicity), and the mechanisms underlying variable drug responses. The identification of DNA variants which markedly contribute to inter-individual variations in drug responses would improve the efficacy of treatments and decrease the rate of the adverse side effects of drugs. This review focuses only on the impact of polymorphisms within drug-metabolizing enzymes on drug responses. Anticancer drugs usually have a very narrow therapeutic index; therefore, it is very important to use appropriate doses in order to achieve the maximum benefits without putting the patient at risk of life-threatening toxicities. However, the adjustment of the appropriate dose is not so easy, due to the inheritance of specific polymorphisms in the genes encoding the target proteins and drug-metabolizing enzymes. This review presents just a few examples of such polymorphisms and their impact on the response to therapy.

Genophenotypic Factors and Pharmacogenomics in Adverse Drug Reactions

Adverse drug reactions (ADRs) rank as one of the top 10 leading causes of death and illness in developed countries. ADRs show differential features depending upon genotype, age, sex, race, pathology, drug category, route of administration, and drug–drug interactions. Pharmacogenomics (PGx) provides the physician effective clues for optimizing drug efficacy and safety in major problems of health such as cardiovascular disease and associated disorders, cancer and brain disorders. Important aspects to be considered are also the impact of immunopharmacogenomics in cutaneous ADRs as well as the influence of genomic factors associated with COVID-19 and vaccination strategies. Major limitations for the routine use of PGx procedures for ADRs prevention are the lack of education and training in physicians and pharmacists, poor characterization of drug-related PGx, unspecific biomarkers of drug efficacy and toxicity, cost-effectiveness, administrative problems in health organizations, and insufficient regulation for the generalized use of PGx in the clinical setting. The implementation of PGx requires: (i) education of physicians and all other parties involved in the use and benefits of PGx; (ii) prospective studies to demonstrate the benefits of PGx genotyping; (iii) standardization of PGx procedures and development of clinical guidelines; (iv) NGS and microarrays to cover genes with high PGx potential; and (v) new regulations for PGx-related drug development and PGx drug labelling.

Protocol for a pharmacogenetic study of antidepressants: characterization of drug-metabolizing profiles of cytochromes CYP2D6 and CYP2C19 in a Sardinian population of patients with major depressive disorder

The effectiveness of antidepressants shows high interindividual variability ranging from full symptomatologic remission to treatment-resistant depression. Many factors can determine the variation in the clinical response, but a fundamental role is played by genetic variation within the genes encoding for the enzymes most involved in the metabolism of antidepressant drugs: the CYP2D6 and CYP2C19 isoforms of the cytochrome P450 system. This study is poised to clarify whether the different metabolizing phenotypes related to CYP2D6 and CYP2C19 could have an impact on the clinical efficacy of antidepressants and whether the frequency of these phenotypes of metabolization shows differences in the population of Sardinian patients compared to other Caucasian populations. The sample is being recruited from patients followed-up and treated at the Psychiatric Unit of the Department of Medical Science and Public Health, University of Cagliari and the University Hospital Agency of Cagliari (Italy). The study design includes three approaches: (1) a pharmacogenetic analysis of 80 patients diagnosed with MDD resistant to antidepressant treatment compared to 80 clinically responsive or remitted patients; (2) a prospective arm (N = 30) of the study where we will test the impact of genetic variation within the CYP2D6 and CYP2C19 genes on clinical response to antidepressants and on their serum levels and (3) the assessment of the socio-economic impact of antidepressant therapies, and estimation of the cost-effectiveness of the pharmacogenetic test based on CYP genes.

Pharmacogene Sequencing of a Gabonese Population with Severe Plasmodium falciparum Malaria Reveals Multiple Novel Variants with Putative Relevance for Antimalarial Treatment

Malaria remains one of the deadliest diseases in Africa, particularly for children. While successful in reducing morbidity and mortality, antimalarial treatments are also a major cause of adverse drug reactions (ADRs). Host genetic variation in genes involved in drug disposition or toxicity constitutes an important determinant of ADR risk and can prime for parasite drug resistance. Importantly, however, the genetic diversity in Africa is substantial, and thus, genetic profiles in one population cannot be reliably extrapolated to other ethnogeographic groups. Gabon is considered a high-transmission country, with more than 460,000 malaria cases per year. Yet the pharmacogenetic landscape of the Gabonese population or its neighboring countries has not been analyzed. Using targeted sequencing, here, we profiled 21 pharmacogenes with importance for antimalarial treatment in 48 Gabonese pediatric patients with severe Plasmodium falciparum malaria. Overall, we identified 347 genetic variants, of which 18 were novel, and each individual was found to carry 87.3 ± 9.2 (standard deviation [SD]) variants across all analyzed genes. Importantly, 16.7% of these variants were population specific, highlighting the need for high-resolution pharmacogenomic profiling. Between one in three and one in six individuals harbored reduced-activity alleles of CYP2A6, CYP2B6, CYP2D6, and CYP2C8 with important implications for artemisinin, chloroquine, and amodiaquine therapy. Furthermore, one in three patients harbored at least one G6PD-deficient allele, suggesting a considerably increased risk of hemolytic anemia upon exposure to aminoquinolines. Combined, our results reveal the unique genetic landscape of the Gabonese population and pinpoint the genetic basis for interindividual differences in antimalarial drug responses and toxicity.

Association of Genetic Polymorphisms in Oxidative Stress and Inflammation Pathways with Glaucoma Risk and Phenotype

Oxidative stress and neuroinflammation are involved in the pathogenesis and progression of glaucoma. Our aim was to evaluate the impact of selected single-nucleotide polymorphisms in inflammation and oxidative stress genes on the risk of glaucoma, the patients’ clinical characteristics and the glaucoma phenotype. In total, 307 patients with primary open-angle glaucoma or ocular hypertension were enrolled. The control group included 339 healthy Slovenian blood donors. DNA was isolated from peripheral blood. Genotyping was performed for SOD2 rs4880, CAT rs1001179, GPX1 rs1050450, GSTP1 rs1695, GSTM1 gene deletion, GSTT1 gene deletion, IL1B rs1143623, IL1B rs16944, IL6 rs1800795 and TNF rs1800629. We found a nominally significant association of GSTM1 gene deletion with decreased risk of ocular hypertension and a protective role of IL1B rs16944 and IL6 rs1800629 in the risk of glaucoma. The CT and TT genotypes of GPX1 rs1050450 were significantly associated with advanced disease, lower intraocular pressure and a larger vertical cup–disc ratio. In conclusion, genetic variability in IL1B and IL6 may be associated with glaucoma risk, while GPX and TNF may be associated with the glaucoma phenotype. In the future, improved knowledge of these pathways has the potential for new strategies and personalised treatment of glaucoma.

Pharmacogenetics Approach for the Improvement of COVID-19 Treatment

The treatment of coronavirus disease 2019 (COVID-19) has been a challenge. The efficacy of several drugs has been evaluated and variability in drug response has been observed. Pharmacogenetics could explain this variation and improve patients’ outcomes with this complex disease; nevertheless, several disease-related issues must be carefully reviewed in the pharmacogenetic study of COVID-19 treatment. We aimed to describe the pharmacogenetic variants reported for drugs used for COVID-19 treatment (remdesivir, oseltamivir, lopinavir, ritonavir, azithromycin, chloroquine, hydroxychloroquine, ivermectin, and dexamethasone). In addition, other factors relevant to the design of pharmacogenetic studies were mentioned. Variants in CYP3A4, CYP3A5, CYP2C8, CY2D6, ABCB1, ABCC2, and SLCO1B1, among other variants, could be included in pharmacogenetic studies of COVID-19 treatment. Besides, nongenetic factors such as drug–drug interactions and inflammation should be considered in the search for personalized therapy of COVID-19.

Prevalence of pharmacogenomic variants in 100 pharmacogenes among Southeast Asian populations under the collaboration of the Southeast Asian Pharmacogenomics Research Network (SEAPharm)

Pharmacogenomics can enhance the outcome of treatment by adopting pharmacogenomic testing to maximize drug efficacy and lower the risk of serious adverse events. Next-generation sequencing (NGS) is a cost-effective technology for genotyping several pharmacogenomic loci at once, thereby increasing publicly available data. A panel of 100 pharmacogenes among Southeast Asian (SEA) populations was resequenced using the NGS platform under the collaboration of the Southeast Asian Pharmacogenomics Research Network (SEAPharm). Here, we present the frequencies of pharmacogenomic variants and the comparison of these pharmacogenomic variants among different SEA populations and other populations used as controls. We investigated the different types of pharmacogenomic variants, especially those that may have a functional impact. Our results provide substantial genetic variations at 100 pharmacogenomic loci among SEA populations that may contribute to interpopulation variability in drug response phenotypes. Correspondingly, this study provides basic information for further pharmacogenomic investigations in SEA populations.

Genetic Diversity of Drug-Related Genes in Native Americans of the Brazilian Amazon

Introduction

The genetic admixture of the Brazilian population has considerable relevance to the implementation of the principles of pharmacogenomics (PGx), as it may compromise the extrapolation of data obtained in more homogeneous world populations.

Purpose

This study aims to investigate a panel of 117 polymorphisms in 35 pharmacogenes, which contains label recommendations or clinical evidence by international drug regulatory agencies, in Amazonian Native American populations, and compare the results obtained with continental population data from the 1000 Genomes Project Consortium.

Patients and Methods

The study population is composed of 109 Native American individuals from three Brazilian Amazon groups. The genotyping of the PGx polymorphisms was performed by allelic discrimination using TaqMan^® OpenArray Genotyping with a panel of 120 customized assays on the QuantStudio™ 12K Flex Real-Time PCR System.

Results

Statistical differences within the Native American populations were observed regarding both genotypes and phenotypes of some genes of the CYP family. The discriminant analysis of principal components (DAPCs) between the NAM group and the continental populations of the 1000 Genomes Project resulted in the clustering of the three Native American populations. Additionally, in general, the NAM group was determined to be closely situated between East Asia, America, and South Asia groups, which enabled us to infer a genetic similarity between these populations. The DAPC analysis further demonstrated that eight polymorphisms and six polymorphisms were more relevant in differentiating the NAM from the continental populations and the NAM populations among themselves, respectively.

Conclusion

Some investigated polymorphisms show differences among world populations, particularly with populations of European origin, for whom precision medicine protocols are primarily designed. The accumulated knowledge regarding these variations may assist in the design of specific protocols for Native American populations and populations admixed with them.

Cartographic atlas of frequency variation for 45 pharmacogenetic markers in populations of Russia and its neighbor states

The lack of information about the frequency of pharmacogenetic markers in Russia impedes the adoption of personalized treatment algorithms originally developed for West European populations. The aim of this paper was to study the distribution of some clinically significant pharmacogenetic markers across Russia. A total of 45 pharmacogenetic markers were selected from a few population genetic datasets, including ADME, drug target and hemostasis-controlling genes. The total number of donors genotyped for these markers was 2,197. The frequencies of these markers were determined for 50 different populations, comprised of 137 ethnic and subethnic groups. A comprehensive pharmacogenetic atlas was created, i.e. a systematic collection of gene geographic maps of frequency variation for 45 pharmacogenetic DNA markers in Russia and its neighbor states. The maps revealed 3 patterns of geographic variation. Clinal variation (a gradient change in frequency along the East-West axis) is observed in the pharmacogenetic markers that follow the main pattern of variation for North Eurasia (13% of the maps). Uniform distribution singles out a group of markers that occur at average frequency in most Russian regions (27% of the maps). Focal variation is observed in the markers that are specific to a certain group of populations and are absent in other regions (60% of the maps). The atlas reveals that the average frequency of the marker and its frequency in individual populations do not indicate the type of its distribution in Russia: a gene geographic map is needed to uncover the pattern of its variation.

Variation in 100 relevant pharmacogenes among emiratis with insights from understudied populations

Genetic variations have an established impact on the pharmacological response. Investigating this variation resulted in a compilation of variants in "pharmacogenes". The emergence of next-generation sequencing facilitated large-scale pharmacogenomic studies and exhibited the extensive variability of pharmacogenes. Some rare and population-specific variants proved to be actionable, suggesting the significance of population pharmacogenomic research. A profound gap exists in the knowledge of pharmacogenomic variants enriched in some populations, including the United Arab Emirates (UAE). The current study aims to explore the landscape of variations in relevant pharmacogenes among healthy Emiratis. Through the resequencing of 100 pharmacogenes for 100 healthy Emiratis, we identified 1243 variants, of which 63% are rare (minor allele frequency ≤ 0.01), and 30% were unique. Filtering the variants according to Pharmacogenomics Knowledge Base (PharmGKB) annotations identified 27 diplotypes and 26 variants with an evident clinical relevance. Comparison with global data illustrated a significant deviation of allele frequencies in the UAE population. Understudied populations display a distinct allelic architecture and various rare and unique variants. We underscored pharmacogenes with the highest variation frequencies and provided investigators with a list of candidate genes for future studies. Population pharmacogenomic studies are imperative during the pursuit of global pharmacogenomics implementation.

Population structure and pharmacogenomic risk stratification in the United States

BACKGROUND

Pharmacogenomic (PGx) variants mediate how individuals respond to medication, and response differences among racial/ethnic groups have been attributed to patterns of PGx diversity. We hypothesized that genetic ancestry (GA) would provide higher resolution for stratifying PGx risk, since it serves as a more reliable surrogate for genetic diversity than self-identified race/ethnicity (SIRE), which includes a substantial social component. We analyzed a cohort of 8628 individuals from the United States (US), for whom we had both SIRE information and whole genome genotypes, with a focus on the three largest SIRE groups in the US: White, Black (African-American), and Hispanic (Latino). Our approach to the question of PGx risk stratification entailed the integration of two distinct methodologies: population genetics and evidence-based medicine. This integrated approach allowed us to consider the clinical implications for the observed patterns of PGx variation found within and between population groups.

RESULTS

Whole genome genotypes were used to characterize individuals' continental ancestry fractions-European, African, and Native American-and individuals were grouped according to their GA profiles. SIRE and GA groups were found to be highly concordant. Continental ancestry predicts individuals' SIRE with > 96% accuracy, and accordingly, GA provides only a marginal increase in resolution for PGx risk stratification. In light of the concordance between SIRE and GA, taken together with the fact that information on SIRE is readily available to clinicians, we evaluated PGx variation between SIRE groups to explore the potential clinical utility of race and ethnicity. PGx variants are highly diverged compared to the genomic background; 82 variants show significant frequency differences among SIRE groups, and genome-wide patterns of PGx variation are almost entirely concordant with SIRE. The vast majority of PGx variation is found within rather than between groups, a well-established fact for almost all genetic variants, which is often taken to argue against the clinical utility of population stratification. Nevertheless, analysis of highly differentiated PGx variants illustrates how SIRE partitions PGx variation based on groups' characteristic ancestry patterns. These cases underscore the extent to which SIRE carries clinically valuable information for stratifying PGx risk among populations, albeit with less utility for predicting individual-level PGx alleles (genotypes), supporting the concept of population pharmacogenomics.

CONCLUSIONS

Perhaps most interestingly, we show that individuals who identify as Black or Hispanic stand to gain far more from the consideration of race/ethnicity in treatment decisions than individuals from the majority White population.

Costing Methods as a Means to Measure the Costs of Pharmacogenomics Testing

Clinical implementation of pharmacogenomics and personalized medicine interventions relies on addressing important financial aspects of the delivery of genetic testing to the patients, be it from public or private providers. Details on how to determine the cost items of the genetic testing are often limited. The goal of this study is to present a costing methodology in order to estimate and measure the costs as far as the technical process of pharmacogenomics testing is concerned. Moreover, an overall cost mindset strategy based on the selective genotyping workflow to guide specialized laboratories of interest effectively is provided. We particularly accounted for the resources consumed within the laboratory premises such as cost of reagents for DNA isolation, cost of consumables, cost of personnel, while costs associated with patient recruitment, blood sample collection and maintenance, administration costs in the hospital, and costs of blood sample shipment were not taken into consideration. Our article presents the first-time detailed information on a costing framework for pharmacogenomic testing that could be employed to laboratories involved in routine clinical implementation of pharmacogenomics.

Pharmacogenetics in Psychiatry: An Update on Clinical Usability

Using pharmacogenetics in guiding drug therapy experiences a steady increase in uptake, although still leads to discussions as to its clinical use. Psychiatry constitutes a field where pharmacogenomic testing might help in guiding drug therapy. To address current challenges, this minireview provides an update regarding genotyping (SNP analysis/arrays/NGS), structural variant detection (star-alleles/CNVs/hybrid alleles), genotype-to-phenotype translations, cost-effectiveness, and actionability of results (FDA/CPIC/PharmGKB) regarding clinical importance of pre-emptive pharmacogenomic testing for prescription of antidepressants and antipsychotics.

Clinical pharmacology applications in clinical drug development and clinical care: A focus on Saudi Arabia

Drug development, from preclinical to clinical studies, is a lengthy and complex process. There is an increased interest in the Kingdom of Saudi Arabia (KSA) to promote innovation, research and local content including clinical trials (Phase I-IV). Currently, there are over 650 registered clinical trials in Saudi Arabia, and this number is expected to increase. An important part of drug development and clinical trials is to assure the safe and effective use of drugs. Clinical pharmacology plays a vital role in informed decision making during the drug development stage as it focuses on the effects of drugs in humans. Disciplines such as pharmacokinetics, pharmacodynamics and pharmacogenomics are components of clinical pharmacology. It is a growing discipline with a range of applications in all phases of drug development, including selecting optimal doses for Phase I, II and III studies, evaluating bioequivalence and biosimilar studies and designing clinical studies. Incorporating clinical pharmacology in research as well as in the requirements of regulatory agencies will improve the drug development process and accelerate the pipeline. Clinical pharmacology is also applied in direct patient care with the goal of personalizing treatment. Tools such as therapeutic drug monitoring, pharmacogenomics and model informed precision dosing are used to optimize dosing for patients at an individual level. In KSA, the science of clinical pharmacology is underutilized and we believe it is important to raise awareness and educate the scientific community and healthcare professionals in terms of its applications and potential. In this review paper, we provide an overview on the use and applications of clinical pharmacology in both drug development and clinical care.

The effect of the VKORC1 promoter variant on warfarin responsiveness in the Saudi WArfarin Pharmacogenetic (SWAP) cohort

Warfarin is a frequently prescribed oral anticoagulant with a narrow therapeutic index, requiring careful dosing and monitoring. However, patients respond with significant inter-individual variability in terms of the dose and responsiveness of warfarin, attributed to genetic polymorphisms within the genes responsible for the pharmacokinetics and pharmacodynamics of warfarin. Extensive warfarin pharmacogenetic studies have been conducted, including studies resulting in genotype-guided dosing guidelines, but few large scale studies have been conducted with the Saudi population. In this study, we report the study design and baseline characteristics of the Saudi WArfarin Pharmacogenomics (SWAP) cohort, as well as the association of the VKORC1 promoter variants with the warfarin dose and the time to a stable INR. In the 936 Saudi patients recruited in the SWAP study, the minor allele C of rs9923231 was significantly associated with a 8.45 mg higher weekly warfarin dose (p value = 4.0 × 10^-46), as well as with a significant delay in achieving a stable INR level. The addition of the rs9923231 status to the model, containing all the significant clinical variables, doubled the warfarin dose explained variance to 31%. The SWAP cohort represents a valuable resource for future research with the objective of identifying rare and prevalent genetic variants, which can be incorporated in personalized anticoagulation therapy for the Saudi population.

Exome-Wide Analysis of the DiscovEHR Cohort Reveals Novel Candidate Pharmacogenomic Variants for Clinical Pharmacogenomics

Recent advances in next-generation sequencing technology have led to the production of an unprecedented volume of genomic data, thus further advancing our understanding of the role of genetic variation in clinical pharmacogenomics. In the present study, we used whole exome sequencing data from 50,726 participants, as derived from the DiscovEHR cohort, to identify pharmacogenomic variants of potential clinical relevance, according to their occurrence within the PharmGKB database. We further assessed the distribution of the identified rare and common pharmacogenomics variants amongst different GnomAD subpopulations. Overall, our findings show that the use of publicly available sequence data, such as the DiscovEHR dataset and GnomAD, provides an opportunity for a deeper understanding of genetic variation in pharmacogenes with direct implications in clinical pharmacogenomics.

Pharmacogenomic Markers of Methotrexate Response in the Consolidation Phase of Pediatric Acute Lymphoblastic Leukemia Treatment

Methotrexate (MTX) is one of the staples of pediatric acute lymphoblastic leukemia (ALL) treatment. MTX targets the folate metabolic pathway (FMP). Abnormal function of the enzymes in FMP, due to genetic aberrations, leads to adverse drug reactions. The aim of this study was to investigate variants in pharmacogenes involved in FMP and their association with MTX pharmacokinetics (MTX elimination profile) and toxicity in the consolidation therapy phase of pediatric ALL patients. Eleven variants in the thymidylate synthetase (TYMS), methylenetetrahydrofolate reductase (MTHFR), dihydrofolate reductase (DHFR), SLC19A1 and SLCO1B genes were analyzed in 148 patients, using PCR- and sequencing-based methodology. For the Serbian and European control groups, data on allele frequency distribution were extracted from in-house and public databases. Our results show that the A allele of SLC19A1 c.80 variant contributes to slow MTX elimination. Additionally, the AA genotype of the same variant is a predictor of MTX-related hepatotoxicity. Patients homozygous for TYMS 6bp deletion were more likely to experience gastrointestinal toxicity. No allele frequency dissimilarity was found for the analyzed variants between Serbian and European populations. Statistical modelling did not show a joint effect of analyzed variants. Our results indicate that SLC19A1 c.80 variant and TYMS 6bp deletion are the most promising pharmacogenomic markers of MTX response in pediatric ALL patients.

The Identification of Admixture Patterns Could Refine Pharmacogenetic Counseling: Analysis of a Population-Based Sample in Mexico

Pharmacogenetic analysis has generated translational data that could be applied to guide treatments according to individual genetic variations. However, pharmacogenetic counseling in some mestizo (admixed) populations may require tailoring to different patterns of admixture. The identification and clustering of individuals with related admixture patterns in such populations could help to refine the practice of pharmacogenetic counseling. This study identifies related groups in a highly admixed population-based sample from Mexico, and analyzes the differential distribution of actionable pharmacogenetic variants. A subsample of 1728 individuals from the Mexican Genomic Database for Addiction Research (MxGDAR/Encodat) was analyzed. Genotyping was performed with the commercial PsychArray BeadChip, genome-wide ancestry was estimated using EIGENSOFT, and model-based clustering was applied to defined admixture groups. Actionable pharmacogenetic variants were identified with a query to the Pharmacogenomics Knowledge Base (PharmGKB) database, and functional prediction using the Variant Effect Predictor (VEP). Allele frequencies were compared with chi-square tests and differentiation was estimated by F_ST. Seven admixture groups were identified in Mexico. Some, like Group 1, Group 4, and Group 5, were found exclusively in certain geographic areas. More than 90% of the individuals, in some groups (Group 1, Group 4 and Group 5) were found in the Central-East and Southeast region of the country. MTRR p.I49M, ABCG2 p.Q141K, CHRNA5 p.D398N, SLCO2B1 rs2851069 show a low degree of differentiation between admixture groups. ANKK1 p.G318R and p.H90R, had the lowest allele frequency of Group 1. The reduction in these alleles reduces the risk of toxicity from anticancer and antihypercholesterolemic drugs. Our analysis identified different admixture patterns and described how they could be used to refine the practice of pharmacogenetic counseling for this admixed population.

Documentation of clinically relevant genomic biomarker allele frequencies in the next-generation FINDbase worldwide database

FINDbase (http://www.findbase.org) is a comprehensive data resource recording the prevalence of clinically relevant genomic variants in various populations worldwide, such as pathogenic variants underlying genetic disorders as well as pharmacogenomic biomarkers that can guide drug treatment. Here, we report significant new developments and technological advancements in the database architecture, leading to a completely revamped database structure, querying interface, accompanied with substantial extensions of data content and curation. In particular, the FINDbase upgrade further improves the user experience by introducing responsive features that support a wide variety of mobile and stationary devices, while enhancing computational runtime due to the use of a modern Javascript framework such as ReactJS. Data collection is significantly enriched, with the data records being divided in a Public and Private version, the latter being accessed on the basis of data contribution, according to the microattribution approach, while the front end was redesigned to support the new functionalities and querying tools. The abovementioned updates further enhance the impact of FINDbase, improve the overall user experience, facilitate further data sharing by microattribution, and strengthen the role of FINDbase as a key resource for personalized medicine applications and personalized public health.

Pharmacogenomic genotypes define genetic ancestry in patients and enable population-specific genomic implementation

The importance of genetic ancestry characterization is increasing in genomic implementation efforts, and clinical pharmacogenomic guidelines are being published that include population-specific recommendations. Our aim was to test the ability of focused clinical pharmacogenomic SNP panels to estimate individual genetic ancestry (IGA) and implement population-specific pharmacogenomic clinical decision-support (CDS) tools. Principle components and STRUCTURE were utilized to assess differences in genetic composition and estimate IGA among 1572 individuals from 1000 Genomes, two independent cohorts of Caucasians and African Americans (AAs), plus a real-world validation population of patients undergoing pharmacogenomic genotyping. We found that clinical pharmacogenomic SNP panels accurately estimate IGA compared to genome-wide genotyping and identify AAs with ≥70 African ancestry (sensitivity >82%, specificity >80%, PPV >95%, NPV >47%). We also validated a new AA-specific warfarin dosing algorithm for patients with ≥70% African ancestry and implemented it at our institution as a novel CDS tool. Consideration of IGA to develop an institutional CDS tool was accomplished to enable population-specific pharmacogenomic guidance at the point-of-care. These capabilities were immediately applied for guidance of warfarin dosing in AAs versus Caucasians, but also provide a real-world model that can be extended to other populations and drugs as actionable genomic evidence accumulates.

Pharmacogenomics landscape of COVID-19 therapy response in Serbian population and comparison with worldwide populations

Background

Since there are no certified therapeutics to treat COVID-19 patients, drug repurposing became important. With lack of time to test individual pharmacogenomics markers, population pharmacogenomics could be helpful in predicting a higher risk of developing adverse reactions and treatment failure in COVID-19 patients. Aim of our study was to identify pharmacogenes and pharmacogenomics markers associated with drugs recommended for COVID-19 treatment, chloroquine/hydroxychloroquine, azithromycin, lopinavir and ritonavir, in population of Serbia and other world populations.

Methods

Genotype information of 143 individuals of Serbian origin was extracted from database previously obtained using TruSight One Gene Panel (Illumina). Genotype data of individuals from different world populations were extracted from the 1000 Genome Project. Fisher's exact test was used for comparison of allele frequencies.

Results

We have identified 11 potential pharmacogenomics markers in 7 pharmacogenes relevant for COVID-19 treatment. Based on high alternative allele frequencies in population and the functional effect of the variants, ABCB1 rs1045642 and rs2032582 could be relevant for reduced clearance of azithromycin, lopinavir and ritonavir drugs and UGT1A7 rs17868323 for hyperbilirubinemia in ritonavir treated COVID-19 patients in Serbian population. SLCO1B1 rs4149056 is a potential marker of lopinavir response, especially in Italian population. Our results confirmed that pharmacogenomics profile of African population is different from the rest of the world.

Conclusions

Considering population specific pharmacogenomics landscape, preemptive testing for pharmacogenes relevant for drugs used in COVID-19 treatment could contribute to better understanding of the inconsistency in therapy response and could be applied to improve the outcome of the COVID-19 patients.

Defining screening panel of functional variants of CYP1A1, CYP2C9, CYP2C19, CYP2D6, and CYP3A4 genes in Serbian population

Plethora of drugs and toxic substances is metabolized by cytochrome P450 enzymes (CYP450). These enzymes are coded by highly variable genes abundant with single nucleotide variants (SNVs) and small insertions/deletions (indels) that affect the functionality of the enzymes, increasing or decreasing their activity. CYP genes genotyping, followed by haplotype inference, provides substrate specific metabolic phenotype prediction. This is crucial in pharmacogenetics and applicable in molecular autopsy. However, high number of alleles in CYP450 superfamily and interethnic variability in frequency distribution require precise gene panel customization. To estimate informativeness of SNVs and alleles in CYP gene families 1, 2, and 3, associated with metabolic alterations, 500 unrelated individuals from 5 regions of Serbia were genotyped using TaqMan assays to determine frequencies of CYP2C9 *2 and *3, CYP2C19 *2 and *17 alleles, four variants in CYP2D6 (rs3892097, rs1065852, rs28371725, rs28371706) gene, and CYP3A4*1B allele. In addition, CYP1A1 rs4646903 and rs1048943 (m1 and m2) variants were genotyped by RFLP. Our results showed that frequencies of tested variants in Serbian population corresponded to general European population and somewhat differed from neighboring populations. SNV rs1065852, the main contributor to non-functional CYP2D6 *4, significantly departed from Hardy-Weinberg equilibrium. With the exception of rs28371706 in CYP2D6 and rs2740574 in CYP3A4, which were very rare in our sample, all other tested variants in CYP2 family are informative and appropriate for pharmacogenetic testing, molecular autopsy, and medico-legal genetic analyses.

Unveiling the guidance heterogeneity for genome-informed drug treatment interventions among regulatory bodies and research consortia

Pharmacogenomics and personalized medicine interventions hold promise to optimize drug treatment modalities and hence, improve the quality of life of the patients by minimizing the occurrence of adverse drug reactions and/or maximizing drug treatment efficacy. To this end, proper guidance for accurately prescribing the correct drug at the right dose is empowered by major regulatory bodies, namely the U.S. Food and Drug Administration (FDA) and the European Medicine Agency (EMA), and well-recognized research consortia, like the Clinical Pharmacogenetics Implementation Consortium (CPIC), that propose therapeutic recommendations after the thorough evaluation of the existing scientific evidence base. In this context, the consistency of these recommendations is crucial for smoothly integrating pharmacogenomics into the clinic. Here, we collected all of the important and clinically actionable pharmacogenomics information provided by the aforementioned renowned sources and documented it in order to assess potential similarities and, most importantly, differences. Our data show that the level of concordance regarding the guidance provided for the same drug-gene association pairs varies significantly, despite the fact that it all derives from a single evidence base. In particular, apart from the expected similarities in a number of association pairs, especially the ones related to cancer genomics, there are still major discrepancies that create confusion as to which guidance should be followed in order to properly inform drug prescribing. This regulatory deficiency calls for the fruitful engagement of the regulatory agencies involved with the contribution of other experts engaged in the field of pharmacogenomics in an effort to harmonize the existing arsenal of guidance for genome-informed drug prescription. The achievement of harmonization would in turn expedite bringing personalized medicine closer to clinical fruition.

Southeast Asian Pharmacogenomics Research Network (SEAPharm): Current Status and Perspectives

Pharmacogenomics (PGx) is increasingly being recognized as a potential tool for improving the efficacy and safety of drug therapy. Therefore, several efforts have been undertaken globally to facilitate the implementation process of PGx into routine clinical practice. Part of these efforts include the formation of PGx working groups working on PGx research, synthesis, and dissemination of PGx data and creation of PGx implementation strategies. In Asia, the Southeast Asian Pharmacogenomics Research Network (SEAPharm) is established to enable and strengthen PGx research among the various PGx communities within but not limited to countries in SEA; with the ultimate goal to support PGx implementation in the region. From the perspective of SEAPharm member countries, there are several key elements essential for PGx implementation at the national level. They include pharmacovigilance database, PGx research, health economics research, dedicated laboratory to support PGx testing for both research and clinical use, structured PGx education, and supportive national health policy. The status of these essential elements is presented here to provide a broad picture of the readiness for PGx implementation among the SEAPharm member countries, and to strengthen the PGx research network and practice in this region.

Frequencies of clinically important CYP2C19 and CYP2D6 alleles are graded across Europe

CYP2C19 and CYP2D6 are important drug-metabolizing enzymes that are involved in the metabolism of around 30% of all medications. Importantly, the corresponding genes are highly polymorphic and these genetic differences contribute to interindividual and interethnic differences in drug pharmacokinetics, response, and toxicity. In this study we systematically analyzed the frequency distribution of clinically relevant CYP2C19 and CYP2D6 alleles across Europe based on data from 82,791 healthy individuals extracted from 79 original publications and, for the first time, provide allele confidence intervals for the general population. We found that frequencies of CYP2D6 gene duplications showed a clear South-East to North-West gradient ranging from <1% in Sweden and Denmark to 6% in Greece and Turkey. In contrast, an inverse distribution was observed for the loss-of-function alleles CYP2D6*4 and CYP2D6*5. Similarly, frequencies of the inactive CYP2C19*2 allele were graded from North-West to South-East Europe. In important contrast to previous work we found that the increased activity allele CYP2C19*17 was most prevalent in Central Europe (25–33%) with lower prevalence in Mediterranean-South Europeans (11–24%). In summary, we provide a detailed European map of common CYP2C19 and CYP2D6 variants and find that frequencies of the most clinically relevant alleles are geographically graded reflective of Europe’s migratory history. These findings emphasize the importance of generating pharmacogenomic data sets with high spatial resolution to improve precision public health across Europe.

Use of Pharmacogenetic Drugs by the Dutch Population

Introduction

The Dutch Pharmacogenetics Working Group (DPWG) indicated a list of actionable genotypes that affect patients’ response to more 50 drugs; these drugs which show variable effects based on patients’ genetic traits were named as pharmacogenetics (PGX) drugs. Preemptive genetic testing before using these drugs may protect certain patients from serious adverse reactions and could help in avoiding treatment failures. The objectives of this study include identifying the rate of PGX drug usage among Dutch population, estimating the level of users who carry the actionable genotypes and determining the main genes involved in drug’s effect variability.

Methods

Usage of PGX drugs over 2011–2017 by the insured population (an average of 11.4 million) in outpatient clinics in Netherlands was obtained from the publically available GIP databank. The data of 45 drugs were analyzed and their interactions with selected pharmacogenes were estimated. Frequency of actionable genotypes of 249 Dutch parents was obtained from the public database: Genome of Netherlands (GoNL), to identify the pattern of genetic characteristics of Dutch population.

Results

Over a 7 year period, 51.3 million exposures of patients to PGX drugs were reported with an average of 5.3 exposures per each drug user. One quarterof the exposures (12.4 million) are predicted to be experienced by individuals with actionable genotypes (risky exposures). Up to 60% of the risky exposures (around 7.5 million) were related to drugs metabolized by CYP2D6. SLCO1B1, and CYP2C19 were also identified among the top genes affecting response of drugs users (involved in about 22 and 12.4% of the risky exposures, respectively). Cardiovascular medications were the top prescribed PGX drug class (43%), followed by gastroenterology (29%) and psychiatry/neurology medications (15%). Women use more PGX drugs than men (55.8 vs. 44.2%, respectively) with the majority (84%) of users in both sexes are above 45 years.

Conclusion

PGX drugs are commonly used in Netherlands. Preemptive panel testing for CYP2D6, SLCO1B1, and CYP2C19 only could be useful to predict 95% of vulnerable patients’ exposures to PGX drugs. Future studies to assess the economic impact of preemptive panel testing on patients of older age are suggested.

Inconsistency in race and ethnic classification in pharmacogenetics studies and its potential clinical implications

Introduction

Racial and ethnic categories are frequently used in pharmacogenetics literature to stratify patients; however, these categories can be inconsistent across different studies. To address the ongoing debate on the applicability of traditional concepts of race and ethnicity in the context of precision medicine, we aimed to review the application of current racial and ethnic categories in pharmacogenetics and its potential impact on clinical care.

Methods

One hundred and three total pharmacogenetics papers involving the CYP2C9, CYP2C19, and CYP2D6 genes were analyzed for their country of origin, racial, and ethnic categories used, and allele frequency data. Correspondence between the major continental racial categories promulgated by National Institutes of Health (NIH) and those reported by the pharmacogenetics papers was evaluated.

Results

The racial and ethnic categories used in the papers we analyzed were highly heterogeneous. In total, we found 66 different racial and ethnic categories used which fall under the NIH race category “White”, 47 different racial and ethnic categories for “Asian”, and 62 different categories for “Black”. The number of categories used varied widely based on country of origin: Japan used the highest number of different categories for “White” with 17, Malaysia used the highest number for “Asian” with 24, and the US used the highest number for “Black” with 28. Significant variation in allele frequency between different ethnic subgroups was identified within 3 major continental racial categories.

Conclusion

Our analysis showed that racial and ethnic classification is highly inconsistent across different papers as well as between different countries. Evidence-based consensus is necessary for optimal use of self-identified race as well as geographical ancestry in pharmacogenetics. Common taxonomy of geographical ancestry which reflects specifics of particular countries and is accepted by the entire scientific community can facilitate reproducible pharmacogenetic research and clinical implementation of its results.