key: cord-0798354-i0mpum4e authors: He, Chunjuan; Peng, Linna; Xing, Shishi; Li, Dandan; Wang, Li; Jin, Tianbo title: Population Genetic Difference of Pharmacogenomic VIP Variants in the Tibetan Population date: 2021-08-16 journal: Pharmgenomics Pers Med DOI: 10.2147/pgpm.s316711 sha: 79a6719fd4213ba040673c8a13bdd6f37b5c6de8 doc_id: 798354 cord_uid: i0mpum4e BACKGROUND: Genetic variation influences drug reaction or adverse prognosis. The purpose of this research was to genotype very important pharmacogenetic (VIP) variants in the Tibetan population. METHODS AND MATERIALS: Blood samples from 200 Tibetans were randomly collected and 59 VIP variants were genotyped, and then compared our data to 26 other populations in the 1000 project to further analyze and identify significant difference. RESULTS: The results showed that on comparing with most of the 26 populations from the 1000 project, rs4291 (ACE), rs1051296 (SLC19A1) and rs1065852 (CYP2D6) significantly differed in the Tibetan population. Furthermore, three significant loci were related to drug response. In addition, the allele frequency of Tibetans least differed from that of East Asian populations, and most differed from that of Americans. CONCLUSION: Three significant loci of variation ACE rs4291, SLC19A1 rs1051296 and CYP2D6 rs1065852 were associated with drug response. This result will contribute to improving the information of the Tibetan in the pharmacogenomics database, and providing a theoretical basis for clinical individualised drug use in Tibetans. A major challenge in current drug clinical practice, drug development and drug regulation is the huge difference in drug treatment 1 which is reasoned to the individual differences for the same drug. 2 Single nucleotide polymorphism (SNP) is an important determinant of individual differences in drug therapy. Pharmacogenomics mainly involves knowledge of pharmacokinetic, pharmacodynamic and focuses on the inheritance of individual variation in the course of drug treatment, clarifies the influence of genomic variation on drug distribution and function, and provides guidance for individualised precise medical treatment 3,4 so on. The Pharmacogenomics Knowledge Base (PharmGKB: http://www.pharmgkb.org) is a useful resource that aims to explain the gene-drug-disease relationship. In recent years, pharmacogenomics research has largely focused on genetic variations that may be related to drug response or metabolism. 5 These variations-very important genetic (VIP) variants are mainly concentrated in the pharmacogenomics database. They assessed the relationship between the VIP variants and specific drugs, and provide individuals with appropriate drugs at proper doses. 6 There are many differences in genetic heterogeneity and genetic polymorphism among different ethnic groups. 7 An important research hotspot in pharmacogenomics is to compare the efficacy of drugs among different races. Tibetan is one of the largest of the 56 ethnic Chinese nation, with a long history. Most Tibetans reside in the Tibet Autonomous Region on the Qinghai-Tibet Plateau, and some people also live in Qinghai, Gansu, Sichuan, Yunnan and other regions of China. The Tibetans have gradually formed unique ethnic custom, lifestyles and eating habits based on the climate and altitude of their residing areas. A study of the genomic variants in the Tibetan population showed that most of the Tibetan gene pool may be differ from the gene pool of the Han people. 8 This study aimed to assess the genotype frequencies of 59 VIP variants between Tibetans and the other 26 ethnicities from the 1000 genome project. The results will contribute to improving the information of the Tibetan pharmacogenomics database, and providing a theoretical basis for drug management in Tibetans. 200 unrelated Tibetans were recruited from the Tibet Autonomous Region. According to the paternal lineage of at least past three generations, all candidates were confirmed to be Tibetan. Based on the results of physical examination, people with chronic diseases were excluded. The protocol of this study was known to all participants and they have signed informed consent. According to the ethics committee of Xizang Minzu University (2019-12), all participants provided blood samples. The procedures met the agency's ethical standards and the 1964 Helsinki Declaration and its subsequent amendments or similar ethical standards. The 59 VIP variants were selected from the PharmGKB database (http://www.pharmgkb.org). The genomic DNA of participants was extracted according to the extraction method of the blood DNA extraction kit (GoldMag Ltd. Xi'an, China), and then its DNA purity and concentration were inspected. Multiplexed SNP MassEXTEND assays were designed by the Agena MassARRAY Assay Design 3.0 software, and the Agena MassARRAY RS1000 was used to genotype the 59 VIP variants. 9 Data arrangement and analysis was performed by Agena Typer 4.0 software. 10, 11 To determine whether the genotype frequencies of variants were in Hardy-Weinberg equilibrium, we used Microsoft Excel and SPSS 22.0 (SPSS, Chicago, IL) for statistical analysis. The genotype frequencies of 59 variants in the Tibetan population were separately compared with those of the other super-populations downloaded from the 1000 genomes project (http://www.internationalgenome.org/), including ①African: the African Caribbeans in Barbados the Utah residents with Northern and Western European ancestry (CEU); the Finnish in Finland (FIN); the British in England and Scotland (GBR); the Iberian populations in Spain (IBS); the Toscani in Italy (TSI); ⑤South Asian: the Bengali in Bangladesh (BEB); the Gujarati Indian in Houston, Texas (GIH); the Indian Telugu in the UK (ITU); the Punjabi in Lahore, Pakistan (PJL) and the Sri Lankan Tamil in the UK (STU). 12 Also, Bonferroni's multiple adjustments were used for significance assessment, and all p values were double-sided when p < 0.05 and p < 0.05/(59 × 27), the p value was considered to be statistically significant. 13 In the Before adjustments, ACE rs4291, SLC19A1 rs1051296 and CYP2D6 rs1065852 significantly differed in the Tibetan population when compared with the other 26 nationalities. The SULT1A1 rs750155, SLC19A1 rs1051298 and SLC19A1 rs1131596 in the Tibetan population were obviously different from those of the other nationalities. Bonferroni adjustments have been made (p < 0.05/(59 × 27)). The results showed that ACE rs4291 and CYP2D6 rs1065852 were significantly different in the Tibetan population when compared to the other 26 nationalities. Except for CHS, SLC19A1 rs1051296 significantly differed in Tibetans when compared with the other 25 ethnic groups. Moreover, SULT1A1 rs750155 and SLC19A1 rs1131596 in the Tibetan population differed from those of the other 22 nationalities. In addition, on the basis of Pharmgkb database (https:// www.pharmgkb.org/), Table 2 showed the drug-related information of significant difference. The rs4291 carriers with different genotypes have different responses to drugs (sertraline, captopril, aspirin and amlodipinechlorthalidonelisinoprilj). Rs1065852 was found to be related to alphahydroxymetoprolol, citalopram escitalopram and iloperidone. As for rs1051298, compared to the allele A individuals, individuals with allele G were associated with increased progression free survival after bevacizumab and pemetrexed treatment. Compared with allele C, rs750155 T allele was not associated with ABT-751 pharmacokinetic parameters in cancer patients receiving ABT-751 treatment. And when comparing to the allele A, rs1131596 G allele was not correlated with the response to methotrexate in children with precursor cell lymphoblastic leukaemia lymphoma and patients with rheumatoid arthritis. The correlation between rs1051296 and clinical drugs has not been reported to date. The MAF distribution map of SNPs ACE (rs4291 T allele), SLC19A1 (rs1051296 A allele), CYP2D6 (rs1065852 A allele), SULT1A1 (rs750155 T allele), SLC19A1 (rs1051298 G allele) and SLC19A1 (rs1131596 A allele) were compared with significant differences between the Tibetans and the other super-populations. As shown in Supplementary Figure 1 , the allele frequency of the Tibetans least differed from the East Asian populations, and most differed from Americans. Among them, the frequencies of rs4291-T and rs1065852-A were the highest in the East Asian population, whereas the frequency of rs1051296-A was the highest in the African population. The frequency of rs1051298-G, rs750155-T and rs1131596-A were the highest in the American population. In recent years, many studies have focused on the efficacy comparison of drug reactions among different races, which lays a foundation for clinical individualised drug use. This study identified 59 VIP loci in Tibetans and compared them with the other 26 different populations. Three significant loci of variation were inferred from the genotyping results, ACE rs4291, SLC19A1 rs1051296 and CYP2D6 rs1065852, which were associated with drug response. Angiotensin-converting enzyme (ACE) encodes an enzyme involved in blood pressure regulation and electrolyte balance. It catalyses the conversion of angiotensin I into the physiologically active peptide angiotensin II. Angiotensin II promotes hypertension through vasoconstriction and salt and water retention, and is helpful for heart remodelling, inflammation, thrombosis and plaque rupture. 14 At present, ACE inhibitors are the first-line treatment for hypertension, which can favourably affect the remodelling of patients with myocardial infarction and heart failure, and reduce the incidence rate and mortality. 15 Therefore, ACE-encoded enzyme is the drug target. 16, 17 Polymorphic loci have been reported to play a role. 18, 19 Martínez-Rodríguez et al 20 reported that rs4291 was associated with an increased risk of hypertension after adjusting for age, gender, BMI, triglyceride, drinking and smoking and with the elevation of ACE enzyme levels. In addition, the function of ACE polymorphisms in Alzheimer's Disease cannot be ignored. The study by de Oliveira et al 21 showed that ACE inhibitors, not related to blood pressure, can delay cognitive decline in ACE haplotype carriers with rs1800764-T and rs4291-A. In 2016, they also observed that each A allele of rs4291 resulted in an increase of 3.074 mg/dL urea and 0.044 mg/dl creatinine per year. The use of ACE inhibitors had a protective effect on the change of creatinine, and had no effect on the change of blood pressure. 22 The above research showed that rs4291 carriers with different genotypes had different drug sensitivity. The CYP2D6 gene encoded Cytochrome P450 Family 2 Subfamily D Member 6 that catalyses drug metabolism and the synthesis of cholesterol, steroids and other lipids. One report has shown that the gene is highly polymorphic in the population. Lee et al 23 26 In our study, rs1065852-A frequency was the highest in the South East Asian populations. Therefore, to evaluate the effect of CYP2D6 (rs1065852) on drug metabolism, we should not ignore the ethnic factors, especially in the Asian populations. The membrane protein Folate carrier protein 1 is involved in the regulation of intracellular folate concentration, and is encoded by SLC19A1 (Solute Carrier Family 19 Member 1) gene. 27 The mutant allele of SLC19A1 −43T>C was reported to be associated with low folate levels. 28 At present, the treatment response to pemetrexed has been proven to be individual specific. Zhang et al investigated the genetic characteristics of pemetrexed response in 203 Han patients with advanced non-small cell lung cancer (NSCLC). 29 The participants who received pemetrexed alone, the SNP rs1051298 of SLC19A1 gene increased the risk of all adverse reactions in different cycles of the treatment. Therefore, rs1051298 may be a marker associated with adverse reactions and efficacy of pemetrexed related therapy in Chinese Han patients. In addition, some researchers have found that rs1051296 was also related to drug reactions. Wang et al explored the influence of the miRNA binding site polymorphism (rs1051296) in SLC19A1 on the serum methotrexate (MTX) concentration in children with acute lymphoblastic leukaemia (ALL) after receiving MTX treatment. 30 In comparison with the GT and TT carriers, the MTX concentration in GG carriers was higher than the treatment threshold. Compared with GT and TT carriers, the delayed elimination of MTX was observed in the GG carriers. Rs1051296 G>T correlated with the MTX plasma concentration. In the research, the frequency of rs1051296 in Tibetans differed from that in the super-populations, leading to different genotypes of carriers such as MTX, for some drug efficacy. The results will improve the pharmacogenomic information of the Tibetan population, and deepen the study on the differences in some genetic polymorphisms between the Tibetan population and the other 26 populations in the world. Some limitations of this study are: a small sample size and lack of experiment verification. Further studies are required on the drug efficacy of drug-related gene polymorphism loci in the Tibetan patients, especially for significant loci. In conclusion, ACE rs4291, SLC19A1 rs1051296 and CYP2D6 rs1065852 were significant in Tibetans compared to the other 26 nationalities. This study supplemented the knowledge of Tibetan pharmacology, and also showed the relationship between the SNPs with significant differences and drugs more perfectly, thus making the clinical medication safer and personalized for the Tibetan population. The datasets used or analysed during the current study are available from the corresponding author upon reasonable request. This study was performed in accordance with the World Medical Association Declaration of Helsinki and was approved by the Ethics Committee of the Affiliated Hospital of Xizang Minzu University. Written informed consent was obtained from all of the subjects before participation. The authors have declared that they agreed to publish. 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. This journal is indexed on the American Chemical Society's Chemical Abstracts Service (CAS). The manuscript management system is completely online and includes a very quick and fair peer-review system, which is all easy to use. Visit http://www.dovepress.com/testimonials.php to read real quotes from published authors. Pharmacogenetics, pharmacogenomics, and individualized medicine Genetic polymorphisms of pharmacogenomic VIP variants in the Uygur population from northwestern China Pharmacogenomics in the clinic Genetic polymorphisms study of pharmacogenomic VIP variants in Han ethnic of China's Shaanxi province Pharmacogenomics-drug disposition, drug targets, and side effects Genetic polymorphisms of pharmacogenomic VIP variants in the Dai population from Yunnan province Identifying novel mutations of NKX2-5 congenital heart disease patients of Chinese minority groups Sequencing of 50 human exomes reveals adaptation to high altitude SNP genotyping using the Sequenom MassARRAY iPLEX platform Genetic polymorphisms of pharmacogenomic VIP variants in the lhoba population of southwest China Polymorphisms and phenotypic analysis of cytochrome P450 3A4 in the Uygur population in northwest China The International HapMap Project Composite variables: when and how The ACE2 gene: its potential as a functional candidate for cardiovascular disease Progression of cardiovascular damage: the role of renin-angiotensin system blockade Angiotensin-converting enzyme 2 and angiotensin 1-7: novel therapeutic targets Focus on receptors for coronaviruses with special reference to angiotensin-Converting enzyme 2 as a potential drug target -A perspective Association of angiotensin-converting enzyme 2 gene polymorphism and enzymatic activity with essential hypertension in different gender: a case-control study Hypertension and hypertensive left ventricular hypertrophy are associated with ACE2 genetic polymorphism Single nucleotide polymorphisms of the angiotensin-converting enzyme (ACE) gene are associated with essential hypertension and increased ACE enzyme levels in Mexican individuals Pharmacogenetics of angiotensin-converting enzyme inhibitors in patients with Alzheimer's disease dementia Pharmacogenetic effects of angiotensin-converting enzyme inhibitors over age-related urea and creatinine variations in patients with dementia due to Alzheimer disease Association of polymorphisms of cytochrome P450 2D6 with blood hydroxychloroquine levels in patients with systemic lupus erythematosus Influence of genetic variants of CYP2D6, CYP2C9, CYP2C19 and CYP3A4 on antiepileptic drug metabolism in pediatric patients with refractory epilepsy Five genetic polymorphisms of cytochrome P450 enzymes in the Czech non-Roma and Czech Roma population samples. Drug Metab Pers Ther Population genetic difference of pharmacogenomic VIP gene variants in the Lisu population from Yunnan Province SLC19A1 pharmacogenomics summary Distinct association of SLC19A1 polymorphism −43T>C with red cell folate levels and of MTHFR polymorphism 677C>T with plasma folate levels Discovery of novel biomarkers of therapeutic responses in han chinese pemetrexed-based treated advanced NSCLC patients Effects of a microRNA binding site polymorphism in SLC19A1 on methotrexate concentrations in Chinese children with acute lymphoblastic leukemia We are grateful to the participants for providing blood samples and extend our thanks to the professionals responsible for blood collection. We also thank the reviewers and editors for their patience. The authors report no conflicts of interest in this work.