key: cord-0024754-62h156lw authors: Wang, Yuliang; Peng, Linna; Lu, Hongyan; Zhang, Zhanhao; Xing, Shishi; Li, Dandan; He, Chunjuan; Jin, Tianbo; Wang, Li title: Genetic Polymorphisms of Very Important Pharmacogene Variants in the Blang Population from Yunnan Province in China date: 2021-12-17 journal: Pharmgenomics Pers Med DOI: 10.2147/pgpm.s327313 sha: 848d9dd5b118215372e02f96c27a6185848b7729 doc_id: 24754 cord_uid: 62h156lw BACKGROUND: We aimed to enrich the pharmacogenomic information of a Blang population (BP) from Yunnan Province in China. METHODS: We genotyped 55 very important pharmacogene (VIP) variants from the PharmGKB database and compared their genotype distribution (GD) in a BP with that of 26 populations by the χ(2) test. The minor allele frequency (MAF) distribution of seven significantly different single-nucleotide polymorphisms (SNPs) was conducted to compare the difference between the BP and 26 other populations. RESULTS: Compared with the GD of 55 loci in the BP, among 26 studied populations, GWD, YRI, GIH, ESN, MSL, TSI, PJL, ACB, FIN and IBS were the top-10 populations, which showed a significantly different GD >35 loci. CHB, JPT, CDX, CHS, and KHV populations had a significantly different GD <20 loci. A GD difference of 27–34 loci was found between the BP and 11 populations (LWK, CEU, ITU, STU, PUR, CLM, GBR, ASW, BEB, MXL and PEL). The GD of five loci (rs750155 (SULT1A1), rs4291 (ACE), rs1051298 (SLC19A1), rs1131596 (SLC19A1) and rs1051296 (SLC19A1)) were the most significantly different in the BP as compared with that of the other 26 populations. The genotype frequency of rs1800764 (ACE) and rs1065852 (CYP2D6) was different in all populations except for PEL and LWK, respectively. MAFs of rs1065852 (CYP2D6) and rs750155 (SULT1A1) showed the largest fluctuation between the BP and SAS, EUR, AFR and AMR populations. CONCLUSION: Our data can provide theoretical guidance for safe and efficacious personalized drug use in the Blang population. The use of drugs should be different among diverse ethnic groups because of differences in ethnicity, age, sex, environmental factors and genetic factors. If these differences are ignored, then drug sensitivity, metabolic rate, and adverse reactions are affected, which influences the curative effect of drugs and aggravates the illness of patients. Genetic factors can explain up to 20-95% of the variability in drug response. 1 Variations in genes can affect the pharmacokinetics/pharmacodynamics of drugs, as well as their absorption and metabolism. "Pharmacogenes" are genes that decide the fate of drug pharmacology in a biological system. In general, pharmacogenes correspond to specific gene "superfamilies". Among numerous gene superfamilies, the cytochrome P450 superfamily is the most widely researched in pharmacogenomics studies. It has been reported that polymorphisms of cytochrome P450 account for the most frequent variations in the phase-I metabolism of drugs. 2 Variations in most gene superfamilies can affect the metabolism of drugs and disease risk. The single-nucleotide polymorphism (SNP) is the most common variation of very important pharmacogenes (VIPs). Usually, SNPs are employed to analyze the pharmacogenomic information in different populations. 3, 4 "Pharmacogenomics" is an emerging approach to "precision medicine". Pharmacogenomics plays a major part in precision medicine by "tailoring" the selection and dosing to the patient's genetic features. 5 The Pharmacogenomics Knowledge Base (PharmGKB; www.pharmgkb.org/) is one of the most commonly used databases on primary pharmacogenomics. PharmGKB contains information on gene-variant annotations, drug-centered pathways, VIPs and diverse diseases. PharmGKB aims to share genotype, phenotype, or other data on genetic variations among researchers. 6 It has been demonstrated that pharmacogenomic analysis of a specific population can aid the efficacious, accurate use of drugs in a population. 7, 8 For example, Bader et al found that variants of the vitamin K epoxide reductase complex gene (VKORC) and cytochrome P450 family 2 subfamily C member 9 gene (CYP2C9), which encode enzymes for warfarin metabolism, were the strongest predictors of variability in the warfarin dose among different populations in Middle East and North Africa. 7 In addition, Kim et al demonstrated that adverse drug reactions could be avoided if preemptive genotyping was employed in a South Korean population. 8 The US Food and Drug Administration (www.fda.gov/) have recognized >250 biomarkers with known pharmacogenomic value, and provided recommendations for therapeutic management. 9 Recently, pharmacogenomics information on increasing numbers of ethnic minorities in China has been explored. For example, Liu et al found that, compared with 11 populations in a dataset from the International HapMap Project (www.genome.gov/), differences in expression between the rs2070676 of the cytochrome P450 family 2 subfamily E member 1 gene (CYP2E1) and rs1065852 of cytochrome P450 family 2 subfamily D member 6 gene (CYP2D6) in people of Zhuang nationality were the greatest according to genotyping of samples of 105 people of Zhuang nationality. 3 Besides, He et al concluded that expression of rs4291 of the angiotensin I-converting enzyme gene (ACE), rs1051296 of the solute carrier family 19 member 1 gene (SLC19A1) and rs1065852 of CYP2D6 differed significantly in a Tibetan population compared with that of 26 other populations after genotyping of 200 samples from a Tibetan population. They also found that the allele frequency in this Tibetan population differed least from that of an East Asian population, and differed most from that of a North American population. 4 China has 56 ethnic groups. The Blang ethnic group is found in Yunnan Province in China. According to the Sixth National Census in 2010, the total number of people of Blang ethnicity was 119,639. Among them, >30,000 people live in Mount Blang, Xiding, Bada, Daluo, Mengman, Menggang and other towns in Menghai County in Xishuangbanna Dai Autonomous Prefecture. 10 People of Blang ethnicity live in mountainous areas with a mild climate and abundant rainfall, which is very conducive to plant growth. The area in which Blang populations live is one of the main raw material-producing areas of "Pu'er tea" and "Mengku tea". Even though genetic studies on Blang populations have been conducted, [10] [11] [12] pharmacogenomics information of the Blang population is lacking. Cheng et al explored the pharmacogenomics information of a Blang population. 13 Here, we shed light on the pharmacogenomic information of a Blang population by genotyping 55 different loci of 27 VIPs using 200 samples from Yunnan Province. These samples are different from those investigated by Cheng and collaborators. We also compared the distribution of genotype frequency and minor allele frequency (MAF) differences (55 loci of 27 VIPs) with a Blang population and 26 other populations. The genetic variations of the 15 gene superfamilies involved in the present study were related mainly to changes in drug metabolism and disease risk. 2, [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] We wished to enrich the pharmacogenomics information of a Blang population and provide a theoretical foundation for promoting the development of personalized precise medication for Blang populations in the future. The study protocol was approved by the Clinical Research Ethics Committee of Xizang Minzu University (Xianyang, China). Written informed consent was obtained from each study participant before a blood sample was given. Participants Two-hundred randomly selected healthy, unrelated individuals of Blang ethnicity from Yunnan Province were recruited. Whole-blood samples were collected according to the study protocol. Candidate participants were healthy individuals and had exclusive Blang ancestry for ≥3 previous generations. People suffering from cancer, infectious diseases, drug/alcohol addiction, severe dysfunction of the heart, liver, or kidney or immune disorders were excluded, as were women who were pregnant or lactating. Thus, the recruited individuals were representative of a Blang population. PharmGKB was used for selection of genetic variants from published polymorphisms associated with VIP variants. Assays for the loci of 55 genetic variants in 27 VIPs were designed. Loci that could not be designed for an assay were excluded. We extracted the genomic DNA from the peripheral blood of participants using the GoldMag ® -Mini Whole Blood Genomic DNA Purification Kit (GoldMag. Xi'an, China) according to manufacturer protocols. The DNA concentration was measured using the NanoDrop 2000C spectrophotometer (Thermo Scientific, Waltham, MA, USA). MassARRAY Assay Design 3.0 (Sequenom, San Diego, CA, USA) was employed to design multiplexed SNP MassEXTEND assays. 28 SNP genotyping was done using MassARRAY RS1000 (Sequenom) according to manufacturer protocols. Sequenom Typer 4.0 was employed to manage and analyze the data on SNP genotyping. 29 The basic information on the selected 55 loci related to 27 VIPs of the Blang population are listed in Table 1 . The polymerase chain reaction (PCR) primers designed for the selected SNPs are shown in Supplemental Table 1 . The basic information comprised the gene name, SNP ID, positions, functional consequence, genotype frequencies and MAF in the Blang population. All samples from the Blang population were genotyped with respect to these variants. PharmGKB was also used for the clinical and variant annotations for seven significantly different SNPs in the Blang population compared with 26 other populations. The genotype data of individuals from 26 populations was obtained from the International HapMap Project An exact test was used to test the frequency validity of each VIP variant by assessing the departure from the Hardy-Weinberg equilibrium. The comparison of genotype frequencies between the Blang population and 26 other populations was conducted using the χ 2 test. SPSS 17.0 (Armonk, NY, USA) and Excel™ (Microsoft, Redmond, WA, USA) were used to analyze the distribution of genotypes and MAFs. The Bonferroni correction was applied to p < 0.05 (two-sided). The VIPs corresponding to 55 loci could be classified into 15 gene superfamilies (Table 1) : cytochrome P450 superfamily; dihydropyrimidine dehydrogenase; prostaglandinendoperoxide synthase; calcium voltage-gated channel; ryanodine receptor; alcohol dehydrogenase; potassium voltage-gated ion channel; N-acetyltransferase; angiotensin I-converting enzyme; potassium inwardly rectifying channel; G-protein coupled receptor family; solute carrier organic anion transporter family; nuclear receptor family; sulfotransferase family; solute carrier family. The sequence function of these 55 loci was classified mainly Among 55 loci, after Bonferroni correction between the Blang population and 26 other populations, the distribution of genotype frequencies was significantly different in five loci: rs750155 of sulfotransferase family 1A member gene (SULT1A1), rs4291 of ACE, rs1051298, rs1131596 and rs1051296 of SLC19A1. Besides, the genotype distribution of rs1800764 (ACE) and rs1065852 (CYP2D6) was different in all populations except for PEL and LWK, respectively. Conversely, the genotype distribution of rs1801028 of the dopamine receptor D2 gene (DRD2) was significantly different only in the GIH population compared with that in the Blang population. In addition to the eight loci mentioned above, the genotype distribution of the remaining loci in the Blang population also showed a significant difference compared with that in the other 26 populations, but to different degrees. The MAF distribution of seven significantly different SNPs is shown in Table 3 and Figure 1 . The MAFs of rs1065852 (CYP2D6) and rs750155 (SULT1A1) showed the greatest similarities among SAS, EUR, AFR and AMR populations, but also showed the largest fluctuation between the Blang population and SAS, EUR, AFR and AMR populations. The MAFs of rs1800764 (ACE) and rs1131596 (SLC19A1) among the seven subpopulations of AFR showed distinct differences when compared with those of the Blang population. However, the MAFs of rs4291 (ACE), rs1051298 (SLC19A1), and rs1051296 (SLC19A1) showed relatively less fluctuation between the Blang population and the other 26 populations. Besides, the MAFs of rs1800764 (ACE) and rs750155 (SULT1A1) in the Blang population were close to those of the PEL population, even though most of other populations showed distinct differences on it. To better observe the phenotypes of these seven significantly different SNPs in the Blang population, their clinical and variant annotations were retrieved from PharmGKB (Supplemental Table 3 and Supplemental Table 4 , respectively). We genotyped 55 VIP variants from PharmGKB and compared the genotype distribution and MAF of variants in a Blang population with those of 26 other populations. Among 55 loci, the genotype distribution of five SNPs (rs750155 (SULT1A1), rs4291 (ACE), rs1051298 (SLC19A1), rs1051296 (SLC19A1) and rs1131596 (SLC19A1)) was significantly different in the Blang population compared with that in the other 26 populations. Two SNPs (rs1800764 (ACE) and rs1065852 (CYP2D6)) showed a significantly different genotype distribution in the Blang population compared with that in the other 25 populations but, compared with PEL and LWK populations, respectively, a significant difference was not observed. In addition, the MAFs of rs1065852 (CYP2D6) and rs750155 (SULT1A1) showed the greatest fluctuation between the Blang population and SAS, EUR, AFR and AMR populations. SULT1A1, encoded by SULT1A1, is an isoform of sulfotransferases. The latter are phase-II detoxification enzymes and have a crucial role in the metabolism of several xenobiotics and endogenous compounds (eg, tamoxifen). 30, 31 High polymorphism of SULT1A1 has been reported among Caucasian, Chinese, African-American and Korean populations. 32, 33 Moyer et al reported that the genetic variation in SULT1A1, including rs750155, which is located in the promoter region (the short arm of chromosome 16) of SULT1A1, could explain (at least in part) the interindividual variability in the onset of menopause and symptoms before initiation of hormone therapy, and may represent a step towards individualizing decisions for hormone therapy. 34 Besides, Innocenti et al demonstrated that allele T of rs750155 is not associated with the pharmacokinetic parameters of ABT-751 (novel anticancer agent) in people with neoplasms as compared with allele C. 35 In our study, the genotype frequency distribution of rs750155 (SULT1A1) in the Blang population was significantly different from that of the other 26 populations. Also, the MAF distribution of rs750155 (SULT1A1) showed the greatest difference between the Blang population and SAS, EUR, AFR and AMR populations. Besides, the allele T frequency of rs750155 was far higher than that of allele C [T (76.7%) vs C (23.3%)], which indicated that the T allele of rs750155 in members of the Blang population with neoplasms could metabolize ABT-751 more readily. ACE, encoded by ACE, is an enzyme that can affect the renin-angiotensin system and regulation of blood pressure. 36, 37 ACE inhibitors are first-line treatment for hypertension. They can favorably affect the vascular remodeling of patients with myocardial infarction and heart failure, and reduce its risk and mortality. 38 The functional SNPs rs1800764 (ACE) and rs4291 (ACE) are located in the promoter region (chromosome 17) of ACE. 39 Linkage disequilibrium has been identified between these two SNPs in ACE in multiple populations. 40, 41 These two SNPs possess the same pharmacokinetic characteristics and are associated with the risk of breast cancer, endstage renal disease and Alzheimer's disease. [42] [43] [44] The SNPs rs1800764 (ACE) and rs4291 (ACE) show different drug responses in different populations. [45] [46] [47] In the present study, the genotype frequency distribution of SNPs rs1800764 (ACE) and rs4291 (ACE) in the Blang population was different from that of the other populations studied, even though rs1800764 (ACE) was not significantly different in the Blang population compared with that in the PEL population. Besides, the MAF of rs1800764 (ACE) in the AFR population showed a distinct difference compared with that in the Blang population. However, rs4291 (ACE) showed relatively less fluctuation of MAF between the Blang population and the other 26 populations. Although the association between SNPs rs1800764 (ACE) and rs4291 (ACE) and the risk of breast cancer, end-stage renal disease and Alzheimer's disease have not been elucidated in the Blang population, our pharmacogenomics study of the SNPs rs1800764 (ACE) and rs4291 (ACE) in the Blang population is important for disease prevention and safe use of drugs. Reduced folate carrier protein 1 (RFC1), encoded by SLC19A1, is a high-capacity, bidirectional transporter of 5-methyl-tetrahydrofolate and thiamine monophosphate. RFC1 is involved in the uptake, homeostasis, folate deficiency as well as the transportation and sensitivity of antifolate chemotherapeutic agents, such as methotrexate. [48] [49] [50] The SNPs rs1051298 and rs1051296 are intron variants, and rs1131596 is the missense variant of SLC19A1. Scholars have postulated genotype (AA + AG) of rs1051298 to be associated with reduced overall survival upon treatment with pemetrexed in people with non-small-cell lung cancer or mesothelioma compared with that with genotype GG. 51 In addition, allele G of rs1051298 has been reported to be associated with longer progression-free survival after treatment with bevacizumab and pemetrexed in patients with lung neoplasms compared with that with allele A of rs1051298. 52 Besides, the SNP rs1051296 is associated with higher plasma concentrations of methotrexate in pediatric patients with acute lymphoblastic leukemia. 53 Evidence suggests that rs1131596 variants have a positive effect on methotrexate toxicity. 54 Research has shown that the SNP rs1131596-G is not associated with alteration of the concentration or side-effects of methotrexate treatment compared with that of the SNP rs1131596-A in Chinese children with precursor cell lymphoblastic leukemia/lymphoma and people with rheumatoid arthritis. 55 In our study, the genotype distribution of rs1051298, rs1051296, and rs1131596 in the Blang population was significantly different from that of the other 26 populations. MAF analyses showed that rs1051298 (SLC19A1), and rs1051296 (SLC19A1) showed relatively less fluctuation between the Blang population and the other 26 populations, even though the MAF of rs1131596 (SLC19A1) in the AFR population showed a distinct difference when compared with that of the Blang population. These observations suggested that pharmacogenomic research of variants of rs1051298, rs1051296 and rs1131596 may help to provide guidance for individualized drug use for the Blang population. CYP2D6, encoded by CYP2D6, is an enzyme of the cytochrome P450 superfamily. It is involved in the metabolism of ≤25% of drugs in common use in the clinic. 56 Debrisoquine and sparteine are CYP2D6 variation-related drugs. 57 The genetic variation of CYP2D6 has been reported to be closely related to the metabolism of antipsychotic, antiarrhythmic and antiepileptic drugs. [58] [59] [60] The SNP rs1065852 is an intron variant of CYP2D6. It is related to alteration of the encoded amino acids of CYP2D6 protein, reduction of CYP2D6 activity and to have a "poor metabolizer" phenotype. 61 In addition, the genotype GG of rs1065852 (CYP2D6) is a factor of increased corrected QT (QTc) interval after treatment with iloperidone in people suffering from schizophrenia. 61 The distribution of rs1065852 (CYP2D6) has been shown to be significantly different in a Zhuang population as compared with that in 11 other ethnic groups by Liu et al. 3 In the present study, the genotype distribution of rs1065852 (CYP2D6) was different in the Blang population when compared with that in all other ethnic groups except for the LWK population, and the MAF distribution showed the largest fluctuation between the Blang population and SAS, EUR, AFR and AMR populations. Hence, the different corrected QTc interval may occur in schizophrenia patients of Blang ethnicity upon treatment with iloperidone. All the above evidence indicated the non-negligible roles of CYP2D6 (rs1065852) in effective drug usage and normal drug metabolism in Blang individuals. We provided information on the genetic polymorphisms of VIP variants in the Blang population from Yunnan Province. Nevertheless, the sample size was small: a much larger sample size is needed to verify our results. The genotype distribution of five SNPs (rs750155 (SULT1A1), rs4291 (ACE), rs1051298 (SLC19A1), rs1051296 (SLC19A1) and rs1131596 (SLC19A1)) was significantly different in the Blang population compared with that in the other 26 populations tested. Two SNPs (rs1800764 (ACE) and rs1065852 (CYP2D6)) showed a significantly different genotype distribution in the Blang population as compared with all other populations tested except for PEL and LWK populations, respectively. The MAF of rs1065852 (CYP2D6) and rs750155 (SULT1A1) showed the largest fluctuation between the Blang population and SAS, EUR, AFR and AMR populations. Our data can provide theoretical guidance for safe and efficacious personalized drug use in the Blang population. VIPs, very important pharmacogenes Blang population GD, genotype distribution Western European ancestry; FIN Toscani in Italy; ACB, African Caribbeans in Barbados Esan in Nigeria; GWD, Gambian in Western Divisions MTX, methotrexate; SULT1A1, sulfotransferase family 1A member 1 ACE, angiotensin I-converting enzyme SLC19A1,solute carrier family 19 Member 1; CYP2D6, cytochrome P450 family 2 subfamily D member 6 CYP2C9, cytochrome P450 family 2 subfamily C member 9 DRD2, dopamine receptor D2; RFC1, reduced folate carrier protein 1 PCR, polymerase chain reaction SNP, single-nucleotide polymorphism; PharmGKB, Pharmacogenomics Knowledge Base Hypothesis: comparisons of interand intra-individual variations can substitute for twin studies in drug research Polymorphism of human cytochrome P450 enzymes and its clinical impact Genetic variation of pharmacogenomic VIP variants in Zhuang nationality of southern China Population genetic difference of pharmacogenomic VIP variants in the Tibetan population Pharmacogenomics and personalized medicine PharmGKB: a worldwide resource for pharmacogenomic information The impact of genetic and non-genetic factors on warfarin dose prediction in MENA region: a systematic review Role of preemptive genotyping in preventing serious adverse drug events in south Korean patients. Drug Safe Pharmacogenomics biomarker discovery and validation for translation in clinical practice HLA-DPB1 polymorphism in Blang and Puyi ethnic groups of southwest China inferred from sequence-based typing The frequency distribution of flavin-containing monooxygenase 3 mutant alleles in 28 populations from Yunnan. Zhonghua Yi Xue Yi Chuan Xue Za Zhi Allele frequencies of 20 autosomal STR loci for 207 unrelated individuals of the Blang people in China Genetic analysis of pharmacogenomic VIP variants in the Blang population from Yunnan Province of China A candidate gene study of capecitabine-related toxicity in colorectal cancer identifies new toxicity variants at DPYD and a putative role for ENOSF1 rather than TYMS Functional characterization of cyclooxygenase-2 polymorphisms Correlations of calcium voltage-gated channel subunit alpha1 A (CACNA1A) gene polymorphisms with benign paroxysmal positional vertigo Clinical and genetic findings in a large cohort of patients with ryanodine receptor 1 gene-associated myopathies Genetic variation in alcohol dehydrogenase is associated with neurocognition in men with HIV and history of alcohol use disorder: preliminary findings Genetic variants of potassium voltage-gated channel genes (KCNQ1, KCNH2, and KCNE1) affected the risk of atrial fibrillation in elderly patients N-acetyltransferase 2 polymorphism is associated with bladder cancer risk: an updated meta-analysis based on 54 case-control studies Human G-protein-coupled inwardly rectifying potassium channel (GIRK1) gene (KCNJ3): localization to chromosome 2 and identification of a simple tandem repeat polymorphism Mutations in the pH-sensing G-protein-coupled receptor GPR68 cause amelogenesis imperfecta Influence of uridine diphosphate glucuronosyltransferase family 1 member A1 and solute carrier organic anion transporter family 1 member B1 polymorphisms and efavirenz on bilirubin disposition in healthy volunteers Nuclear receptor variants in liver disease Sulfotransferase 1A1 Arg(213)His polymorphism and prostate cancer risk Associations between polymorphisms in the solute carrier family 6 member 3 and the myelin basic protein gene and posttraumatic stress disorder Angiotensin I converting enzyme gene polymorphism and coronary heart disease SNP genotyping using the Sequenom massARRAY iPLEX platform High-throughput oncogene mutation profiling in human cancer Pharmacogenetics of human cytosolic sulfotransferases Pharmacogenetics of SULT1A1 Sulfation pharmacogenetics: SULT1A1 and SULT1A2 allele frequencies in Caucasian, Chinese and African-American subjects Genetic polymorphisms and linkage disequilibrium of sulfotransferase SULT1A1 and SULT1A2 in a Korean population: comparison of other ethnic groups Influence of SULT1A1 genetic variation on age at menopause, estrogen levels, and response to hormone therapy in recently postmenopausal white women Preclinical discovery of candidate genes to guide pharmacogenetics during Phase I development: the example of the novel anticancer agent ABT-751 Brainpenetrating angiotensin-converting enzyme inhibitors and cognitive change in patients with dementia due to Alzheimer's disease ACE insertion/deletion (I/D) polymorphism and diabetic nephropathy Progression of cardiovascular damage: the role of renin-angiotensin system blockade Pharmacogenetics of angiotensin-converting enzyme inhibitors in patients with Alzheimer's disease dementia Haplotype structure of five SNPs within the ACE gene in the Tunisian population The linkage disequilibrium pattern of the angiotensin converting enzyme gene in Arabic and Asian population groups Synergistic association of six well-characterized polymorphisms in three genes of the renin-angiotensin system with breast cancer among Han Chinese women Gene-gene interactions in renin-angiotensin-aldosterone system contributes to end-stage renal disease susceptibility in a Han Chinese population Angiotensin-converting enzyme levels and activity in Alzheimer's disease: differences in brain and CSF ACE and association with ACE1 genotypes Association of angiotensin I-converting enzyme gene polymorphisms with aspirin intolerance in asthmatics Identification of specific angiotensin-converting enzyme variants and haplotypes that confer risk and protection against type 2 diabetic nephropathy Pharmacogenetic effects of angiotensin-converting enzyme inhibitors over age-related urea and creatinine variations in patients with dementia due to Alzheimer disease. Colombia medica A novel cDNA restores reduced folate carrier activity and methotrexate sensitivity to transport deficient cells Carrier-and receptor-mediated transport of folate antagonists targeting folate-dependent enzymes: correlates of molecular-structure and biological activity The reduced folate carrier (RFC) is cytotoxic to cells under conditions of severe folate deprivation. RFC as a double edged sword in folate homeostasis Pharmacogenetics of pemetrexed combination therapy in lung cancer: pathway analysis reveals novel toxicity associations Phase II trial of pemetrexed plus bevacizumab for second-line therapy of patients with advanced non-small-cell lung cancer: NCCTG and SWOG study N0426 MiR-595 suppresses the cellular uptake and cytotoxic effects of methotrexate by targeting SLC19A1 in CEM/C1 cells Transcription regulatory polymorphism −43T>C in the 5ʹ-flanking region of SLC19A1 gene could affect rheumatoid arthritis patient response to methotrexate therapy Polymorphisms in methotrexate transporters and their relationship to plasma methotrexate levels, toxicity of high-dose methotrexate, and outcome of pediatric acute lymphoblastic leukemia Influence of cytochrome P450 polymorphisms on drug therapies: pharmacogenetic, pharmacoepigenetic and clinical aspects The genetic control of sparteine and debrisoquine metabolism in man with new methods of analysing bimodal distributions Effect of routine cytochrome P450 2D6 and 2C19 genotyping on antipsychotic drug persistence in patients with schizophrenia: a randomized clinical trial Psychiatric aspects of chloroquine and hydroxychloroquine treatment in the wake of Coronavirus disease-2019: psychopharmacological interactions and neuropsychiatric sequelae Influence of genetic variants of CYP2D6, CYP2C9, CYP2C19 and CYP3A4 on antiepileptic drug metabolism in pediatric patients with refractory epilepsy A thorough QTc study of 3 doses of iloperidone including metabolic inhibition via CYP2D6 and/or CYP3A4 and a comparison to quetiapine and ziprasidone Pharmacogenomics and Personalized Medicine is an international, peer-reviewed, open access journal characterizing the influence of genotype on pharmacology leading to the development of personalized treatment programs and individualized drug selection for improved safety, efficacy and sustainability We express our thanks to all study participants. We also thank the clinicians and hospital staff who worked on sample/data collection in this study. The study was supported by the Talent The authors declare that they have no competing interests. All relevant data are available within the manuscript. Scholars interested in other information from this study should contact the corresponding author. All experiments were conducted in accordance with the Declaration of Helsinki 1964 and its later amendments. Each participant provided written informed consent before study commencement. The study protocol was approved (2019-12) by the Ethics Committee of Xizang Minzu University. This work was performed in collaboration between all authors. YLW and LNP carried out the draft and improvement of the manuscript. HYL, ZHZ and SSX designed the tables and figures. DDL and CJH performed the SNP genotyping analysis. TBJ and LW conceived of the study, worked on associated data collection and statistical analysis, participated in the coordination and funded of the study. All authors contributed to data analysis, drafting or revising the article, have agreed on the journal to which the article will be submitted, gave final approval of the version to be published, and agree to be accountable for all aspects of the work. YLW and LNP contributed equally to this article. Yuliang Wang and Linna Peng are co-first authors.