HLA-B*46:01:01:01 and HLA-DRB1*09:01:02:01 are associated with anti-rHuEPO-induced pure red cell aplasia

Treatment of anemia in patients with chronic kidney disease (CKD) with recombinant human erythropoietin (rHuEPO) can be disrupted by a severe complication, anti-rHuEPO-induced pure red cell aplasia (PRCA). Specific HLA genotypes may have played a role in the high incidence of PRCA in Thai patients (1.7/1,000 patient years vs. 0.03/10,000 patient years in Caucasians). We conducted a case-control study in 157 CKD patients with anti-rHuEPO-induced PRCA and 56 controls. The HLA typing was determined by sequencing using a highly accurate multiplex single-molecule, real-time, long-read sequencing platform. Four analytical models were deployed: Model 1 (additive: accounts for the number of alleles), Model 2 (dominant: accounts for only the presence or absence of alleles), Model 3 (adjusted additive with rHuEPO types) and Model 4 (adjusted dominant with rHuEPO types). HLA-B*46:01:01:01 and DRB1*09:01:02:01 were found to be independent risk markers for anti-rHuEPO-induced PRCA in all models [OR (95%CI), p-values for B*46:01:01:01: 4.58 (1.55-13.51), 0.006; 4.63 (1.56-13.75), 0.006; 5.72 (1.67-19.67), 0.006; and 5.81 (1.68-20.09), 0.005; for DRB1*09:01:02:01: 3.99 (1.28-12.49), 0.017, 4.50 (1.32-15.40), 0.016, 3.42 (1.09-10.74), 0.035, and 3.75 (1.08-13.07), 0.038, in Models 1-4, respectively. HLA-B*46:01:01:01 and DRB1*09:01:02:01 are susceptible alleles for anti-rHuEPO-induced PRCA. These findings support the role of HLA genotyping in helping to monitor patients receiving rHuEPO treatment.

Effects of polymorphisms in the MTHFR gene on 5-FU hematological toxicity and efficacy in Thai colorectal cancer patients

Background

The two common methylenetetrahydrofolate reductase (MTHFR) polymorphisms 677G>A and 1298A>C may have been affecting 5-FU toxicity in cancer patients for decades. Drug efficacy has also been shown by previous studies to be affected. In this study, we investigated the effects of these polymorphisms on 5-FU hematological toxicity and treatment efficacy, to provide enhanced pharmacological treatment for cancer patients.

Methods

This is a retrospective study involving 52 Thai colorectal cancer patients who were treated with 5-FU based therapy, using TaqMAN real-time PCR to genotype the MTHFR polymorphisms (677G>A and 1298A>C). The toxicity and response rate were assessed using standardized measures.

Results

Neutropenia was significantly more likely to be experienced (P=0.049, OR=7.286, 95% CI=0.697-76.181) by patients with the MTHFR 677G>A polymorphism, in the same way as leukopenia (P =0.036, OR=3.333, 95%CI=2.183-5.090) and thrombocytopenia (P<0.001, OR=3.917, 95%CI=2.404-6.382). The MTHFR 1298A>C polymorphism had no statistical association with hematological toxicity in 5-FU treatment. The response rate to 5-FU was not significantly affected by these two polymorphisms.

Conclusion

The MTHFR polymorphism 677G>A is a significant risk factor for developing leukopenia, neutropenia and thrombocytopenia as toxic effects of 5-FU therapy in cancer patients. Therefore, patients receiving 5-FU-based therapy should be aware of their polymorphisms as one risk factor for experiencing severe toxicity.

Associations between UGT1A1 and SLCO1B1 polymorphisms and susceptibility to neonatal hyperbilirubinemia in Thai population

Hyperbilirubinemia is the main mechanism that causes neonatal jaundice, and genetics is one of the risk factors of hyperbilirubinemia. Therefore, this study aims to explore the correlation between two genes, UGT1A1 and SLCO1B1, and hyperbilirubinemia in Thai neonates. One hundred thirty seven neonates were recruited from Division of Clinical Chemistry, Ramathibodi Hospital. UGT1A1*28 and *6 were determined by pyrosequencing whereas, SLCO1B1 388A > G and 521 T > C genetic variants were determined by TaqMan® real-time polymerase chain reaction. Neonates carrying with homozygous (AA) and heterozygous (GA) variants in UGT1A1*6 were significantly related to hyperbilirubinemia development compared with wild type (GG; P < 0.001). To the combined of UGT1A1, total bilirubin levels in homozygous variant were higher significantly than heterozygous variant and wild type (P = 0.002, P = 0.003, respectively). Moreover, SLCO1B1 combination was significant differences between the hyperbilirubinemia and the control group (P = 0.041). SLCO1B1 521 T > C variant provide protection for Thai neonatal hyperbilirubinemia (P = 0.041). There are no significant differences in UGT1A1*28 and SLCO1B1 388A > G for the different severity of hyperbilirubinemia. The combined UGT1A1*28 and *6 polymorphism is a strong risk factor for the development of severe hyperbilirubinemia in Thai neonates. Therefore, we suggest neonates with this gene should be closely observed to avoid higher severities of bilirubin.

Genotypic and phenotypic landscapes of 51 pharmacogenes derived from whole-genome sequencing in a Thai population

Differences in drug responses in individuals are partly due to genetic variations in pharmacogenes, which differ among populations. Here, genome sequencing of 171 unrelated Thai individuals from all regions of Thailand was used to call star alleles of 51 pharmacogenes by Stargazer, determine allele and genotype frequencies, predict phenotype and compare high-impact variant frequencies between Thai and other populations. Three control genes, EGFR, VDR, and RYR1, were used, giving consistent results. Every individual had at least three genes with variant or altered phenotype. Forty of the 51 pharmacogenes had at least one individual with variant or altered phenotype. Moreover, thirteen genes had at least 25% of individuals with variant or altered phenotype including SLCO1B3 (97.08%), CYP3A5 (88.3%), CYP2C19 (60.82%), CYP2A6 (60.2%), SULT1A1 (56.14%), G6PD (54.39%), CYP4B1 (50.00%), CYP2D6 (48.65%), CYP2F1 (46.41%), NAT2 (40.35%), SLCO2B1 (28.95%), UGT1A1 (28.07%), and SLCO1B1 (26.79%). Allele frequencies of high impact variants from our samples were most similar to East Asian. Remarkably, we identified twenty predicted high impact variants which have not previously been reported. Our results provide information that contributes to the implementation of pharmacogenetic testing in Thailand and other Southeast Asian countries, bringing a step closer to personalized medicine.

Influence of SULT1A1*2 Polymorphism on Plasma Efavirenz Concentration in Thai HIV-1 Patients

Purpose

Plasma efavirenz (EFV) concentrations within therapeutic levels are essential to successfully treat patients suffering from human immunodeficiency virus (HIV) type 1. In addition to the drug-metabolizing enzyme CYP2B6, other phase II drug-metabolizing enzymes and transporters may have an important role in the pharmacokinetics of EFV. Thus, the influence of phase II drug-metabolizing enzymes and drug transporters on plasma EFV levels was investigated in Thai HIV patients receiving EFV.

Patients and Methods

Genotyping was performed by TaqMan^® real-time PCR in 149 HIV-infected Thai adults, and plasma efavirenz concentration was measured by a validated high-performance liquid chromatography in 12 hours after dosing steady-state plasma samples at week 12 and 24.

Results

Patients with three or more copies of SULT1A1 had significantly lower median plasma EFV concentrations than those carrying two copies at week 12 (p=0.046) and SULT1A1*2 (c.638G>A) carriers had significantly lower median plasma EFV concentrations compared to those not carrying the variant at week 24 (p=0.048). However, no significant association was found after adjusting for CYP2B6 genotype.

Conclusion

Genetic variation in a combination of SULT1A1*2 and SULT1A1 copy number may contribute to variability in EFV metabolism and thereby may impact drug response. The influence of a combination between the SULT1A1 and CYP2B6 genotype on EFV pharmacokinetics should be further investigated in a larger study population.

Phenotype prediction and characterization of 25 pharmacogenes in Thais from whole genome sequencing for clinical implementation

Publicly available pharmacogenomics (PGx) databases enable translation of genotype data into clinically actionable information. As variation within pharmacogenes is population-specific, this study investigated the spectrum of 25 clinically relevant pharmacogenes in the Thai population (n = 291) from whole genome sequencing. The bioinformatics tool Stargazer was used for phenotype prediction, through assignment of alleles and detection of structural variation. Known and unreported potentially deleterious PGx variants were identified. Over 25% of Thais carried a high-risk diplotype in CYP3A5, CYP2C19, CYP2D6, NAT2, SLCO1B1, and UGT1A1. CYP2D6 structural variants accounted for 83.8% of all high-risk diplotypes. Of 39 known PGx variants identified, six variants associated with adverse drug reactions were common. Allele frequencies of CYP3A5*3 (rs776746), CYP2B6*6 (rs2279343), and NAT2 (rs1041983) were significantly higher in Thais than East-Asian and global populations. 121 unreported variants had potential to exert clinical impact, majority were rare and population-specific, with 60.3% of variants absent from gnomAD database. This study demonstrates the population-specific variation in clinically relevant pharmacogenes, the importance of CYP2D6 structural variation detection in the Thai population, and potential of unreported variants in explaining drug response. These findings are essential in development of dosing guidelines, PGx testing, clinical trials, and drugs.

CYP2D6 Predicts Plasma Donepezil Concentrations in a Cohort of Thai Patients with Mild to Moderate Dementia

Purpose

Donepezil, a drug frequently used to treat dementia, is mainly metabolized by cytochrome P450 2D6 (CYP2D6). This study investigated the relationships between CYP2D6 genotype and activity scores as well as predicted phenotype of plasma donepezil concentrations in 86 Thai dementia participants.

Materials and Methods

CYP2D6 was genotyped using bead-chip technology (Luminex xTAG^® v.3). Steady-state trough plasma donepezil concentrations were measured using high-performance liquid chromatography.

Results

Sixteen genotypes were found but the most frequent genotypes detected among our participants were CYP2D6*10/*10 (27.9%) and *1/*10 (26.7%). One-third of the participants had an activity score of 1.25 which predicted that they were normal metabolizers. The overall median (interquartile range) of plasma donepezil concentration was 51.20 (32.59–87.24) ng/mL. Normal metabolizers (NMs) had lower plasma donepezil concentrations compared to intermediate metabolizers (IMs) (41.15 (28.44–67.65) ng/mL vs 61.95 (35.25–97.00) ng/mL). Multivariate analysis showed that CYP2D6 activity score (r² = 0.50) and the predicted phenotype (independent of dose) could predict the plasma donepezil concentration (r² = 0.49).

Conclusion

Plasma donepezil concentration in NMs was lower compared to IMs. Additional studies with larger sample size and use of next-generation sequencing as well as its outcomes are warranted to confirm the benefit of using pharmacogenetic-guided treatment for donepezil.

SLCO1B1 and ABCG2 Gene Polymorphisms in a Thai Population

Introduction

Genetic polymorphisms of drug transporters influence drug transporter activity and alter pharmacokinetic profiles of the drugs. Organic anion transporting polypeptide 1B1 (OATP1B1) and breast cancer resistance protein (BCRP) are important transporters encoded by solute carrier organic anion transporter family member 1B1 (SLCO1B1) gene and ATP-binding cassette subfamily G member 2 (ABCG2) gene, respectively. Polymorphisms in these genes are associated with increased plasma statins concentrations, statin-induced myopathy and poor response to allopurinol treatment.

Purpose

We explored allele and genotype frequencies of SLCO1B1 and ABCG2 genes including their predicted phenotypes in 53 Thai participants. Of these, 17 had chronic kidney disease and were on statins.

Materials and Methods

Genotyping analysis for SLCO1B1 c.521T>C (rs4149056), c.388A>G (rs2306283), g.-11187G>A (rs4149015), and ABCG2 c.421C>A (rs2231142) was done by using TaqMan^® Real time PCR. All were tested for Hardy–Weinberg Equilibrium.

Results

Most of the participants (80%) had normal function haplotypes SLCO1B1 (*1A and *1B) while decreased (*5, *15, and *17) and unknown (*21) function haplotypes were less observed. Four phenotypes of SLCO1B1 were observed: 69.81% had normal function (*1A/*1A,*1A/*1B, and *1B/*1B), 13.21% had intermediate function (*1A/*17, *1B/*15 and *1B/*17), 9.43% had indeterminate function (*1A/*21 and *1B/*21) and 7.55% had low function (*5/*15, *15/*15, and *15/*17). ABCG2 c.421A allele frequency was 25%. The frequency of ABCG2 c.421CA and AA phenotypes were 37.7% and 5.7%, respectively. The allele and genotype frequencies observed are consistent with reports in Asians. However, there were differences in major allele distributions between Asians and Caucasians for SLCO1B1 c.388A>G; SLCO1B1 c.388G were highly found in Asians, but c.388A were more observed in Caucasians.

Conclusion

This study showed that in the Thai population, there were 4 SNPs of SLCO1B1 and ABCG2 genes. This finding may be clinically applied to minimize inter-individual variability of drugs such as statins and allopurinol. Further study with a larger sample size is needed to assess the drug profiles and responses to treatment.

Pharmacogene Variation in Thai Plasmodium vivax Relapse Patients Treated with a Combination of Primaquine and Chloroquine

Purpose

Pharmacogenes have an influence on biotransformation pathway and clinical outcome of primaquine and chloroquine which are often prescribed to treat Plasmodium vivax infection. Genetic variation may impact enzyme activity and/or transporter function and thereby contribute to relapse. The aim of the study was to assess allele, genotype frequencies and the association between pharmacogenes variation and primaquine response in Thai patients infected with Plasmodium vivax.

Patients and Methods

Fifty-one patients were genotyped for 74 variants in 18 genes by Sequenom MassARRAY^® and Taqman^® SNP Real-Time PCR.

Results

SNP frequencies were not significantly different between relapse (n=4) and non-relapse (n=47) patients. However, the CYP2C19 c.681G>A, the frequency of the A-allele that defines the non-functional CYP2C19*2 haplotype was significantly higher compared to the G-allele (OR=5.14, p=0.021). Patients heterozygous for ABCG2 c.421C>A had a higher odds ratio (OR=8.75, p=0.071) and the frequency of the G-allele of UGT2B7 c.372G>A was higher compared to the A-allele (OR=3.75, p=0.081). CYP2C19, ABCG2 and UGT2B7 emerged as potential high priority genes.

Conclusion

Decreased activity of CYP2C19, ABCG2 and UGT2B7 in combination with CYP2D6 intermediate or poor metabolizer status may expose patients to a higher risk of Plasmodium vivax relapse. Further investigations are warranted to substantiate these findings.

CYP2D6 genotype analysis of a Thai population: platform comparison

The highly polymorphic CYP2D6 gene locus leads to a wide range of enzyme activity. Since there are limited data for Thai, the major aim was to investigate CYP2D6 genetic variation in a large Thai population (n = 920). CYP2D6 genotyping was performed using four different platforms. Genotype call rates of the Luminex xTAG® and AmpliChip CYP450 test were 96.5% and 87.4%, respectively. Based on Luminex xTAG® data, the most common alleles and genotypes were *1 0 (49.6%), *1 (24.6%), *2 (10.8%), *5 (6.7%), *41 (6.5%) and *1/*10 (23.9%), *10/*10 (21.5%), *2/*10 (9.4%), *5/*10 (6.9%), *10/*41 (5.7%), respectively. Challenges and limitations of the platforms evaluated are discussed. These data add to our knowledge regarding interethnic variability in CYP2D6 activity and contribute to improving drug therapy in the Thai.

Pharmacogenomics and Efficacy of Risperidone Long-Term Treatment in Thai Autistic Children and Adolescents

The purpose of this study was to evaluate the association of pharmacogenomic factors and clinical outcome in autistic children and adolescents who were treated with risperidone for long periods. Eighty-two autistic subjects diagnosed with DSM-IV and who were treated with risperidone for more than 1 year were recruited. Pharmacogenomics and clinical outcome (CGI-I, aggressive, overactivity and repetitive score) were evaluated. Almost all patients showed stable symptoms on aggressive behaviour (89.02%), overactivity (71.95%), repetitive (70.89%) behaviour and all clinical symptoms (81.71%). Only 4.48% of patients showed minimally worse CGI-I score. Patients in the non-stable symptom group had DRD2 Taq1A non-wild-type (TT and CT) frequencies higher than the clinically stable group (p = 0.04), whereas other gene polymorphisms showed no significant association. Haplotype ACCTCAT (rs6311, rs1045642, rs1128503, rs1800497, rs4436578, rs1799978, rs6280) showed a significant association with non-stable clinical outcome (χ²  = 6.642, p = 0.010). Risperidone levels showed no association with any clinical outcome. On the other hand, risperidone dose, 9-OH risperidone levels and prolactin levels were significantly higher in the non-stable compared to the stable symptom group (p = 0.013, p = 0.044, p = 0.030). Increased appetite was the most common adverse drug reaction and associated with higher body-weight, whereas it was not significantly associated with genetic variations and non-genetic information. In conclusion, risperidone showed efficacy to control autism, especially aggressive symptoms in long-term treatment. However, Taq1A T - carrier of dopamine 2 receptor gene - is associated with non-stable response in risperidone-treated patients. This study supports pharmacogenomics testing for personalized therapy with risperidone in autistic children and adolescents.

Pharmacogenomics and Efficacy of Risperidone Long-term Treatment in Thai Autistic Children and Adolescence

Risperidone, atypical antipsychotics, was approved for irritability in autistic disorder. However, some patients had minimal improvement or no response to this treatment. The purpose of this study was to evaluate the association of pharmacogenomics factors and clinical outcomes in autistic children and adolescence who treated with risperidone for long periods. Sixty-seven autistic subjects diagnosed with DSM-IV criteria and treated with risperidone more than 1 year were evaluated clinical symptom by CGI, aggressive, over activity, and repetitive score. Polymorphisms of ABCB1, CYP2D6, DRD2, DRD3, and HTR2A were analyzed. Almost patients showed stable symptom on aggressive (91.04%), over activity (73.13%), repetitive (68.25%) behavior, and all clinical symptoms (82.09%). Only 4.48% of patients showed minimally worse on CGI-I score. Patients in non-stable of all symptom group had DRD2 Taq1A non-wildtype (TT and CT) frequencies higher than clinical stable group (P = 0.046), whereas other genes polymorphism showed no significant association. Interestingly, there was no patient with HTR2A-1438G > A wildtype in all non-stable symptoms. However, there was no significant association due to small sample sizes. Drug levels (RIS, 9OH-RIS, and active moiety) did not show the association with any clinical outcome. Increased appetite was the common ADRs, which associated with high body weight, whereas there was not significantly associated with genetic variations and non-genetic information. In conclusion, risperidone showed efficacy to control autism, especially aggressive symptom in long-term treatment. However, dopamine 2 gene variation affect to non-stable in risperidone treated patients. This study supports pharmacogenomics testing for personalized therapeutics of risperidone in autistic disorder.

Disclosure of interest

The authors have not supplied their declaration of competing interest.

Pharmacogenomic Study Reveals New Variants of Drug Metabolizing Enzyme and Transporter Genes Associated with Steady-State Plasma Concentrations of Risperidone and 9-Hydroxyrisperidone in Thai Autism Spectrum Disorder Patients

The present study sought to investigate the genetic variants in drug metabolizing enzyme and transporter (DMET) genes associated with steady-state plasma concentrations of risperidone among Thai autism spectrum disorder (ASD) patients. ASD patients taking risperidone for at least 1 month were enrolled for this pharmacogenomic study. Genotyping profile was obtained using Affymetrix DMET Plus array interrogating 1931 variants in 231 genes. Steady-state plasma risperidone and 9-hydroxyrisperidone were measured using liquid chromatography/tandem mass spectrometry assay. The final analysis included 483 markers for 167 genes. Six variants, ABCB11 (c.3084A > G, c.^∗420A > G, c.^∗368G > A, and c.^∗236G > A) and ADH7 (c.690G > A and c.-5360G > A), were found to be associated with plasma concentrations of risperidone. 9-Hydroxyrisperidone and the total active-moiety levels were associated with six gene variants, SCLO1B1 (c.-11187G > A and c.521T > C), SLCO1B3 (c.334G > T, c.699A > G, and c.1557G > A), and SLC7A5 c.^∗438C > G. Polymorphisms in UGT2B4 c.^∗448A > G and CYP2D6 (c.1661G > C, c.4180G > C, and c.-2178G > A) showed considerable but not significant associations with metabolic ratio. This pharmacogenomic study identifies new genetic variants of DMET genes in monitoring risperidone therapy.

CYP2D6 polymorphisms and their influence on risperidone treatment

Cytochrome P450 enzyme especially CYP2D6 plays a major role in biotransformation. The interindividual variations of treatment response and toxicity are influenced by the polymorphisms of this enzyme. This review emphasizes the effect of CYP2D6 polymorphisms in risperidone treatment in terms of basic knowledge, pharmacogenetics, effectiveness, adverse events, and clinical practice. Although the previous studies showed different results, the effective responses in risperidone treatment depend on the CYP2D6 polymorphisms. Several studies suggested that CYP2D6 polymorphisms were associated with plasma concentration of risperidone, 9-hydroxyrisperidone, and active moiety but did not impact on clinical outcomes. In addition, CYP2D6 poor metabolizer showed more serious adverse events such as weight gain and prolactin than other predicted phenotype groups. The knowledge of pharmacogenomics of CYP2D6 in risperidone treatment is increasing, and it can be used for the development of personalized medication in term of genetic-based dose recommendation. Moreover, the effects of many factors in risperidone treatment are still being investigated. Both the CYP2D6 genotyping and therapeutic drug monitoring are the important steps to complement the genetic-based risperidone treatment.