Human loci involved in drug biotransformation: worldwide genetic variation, population structure, and pharmacogenetic implications

Joint Analysis of Phenotypic and Genomic Diversity Sheds Light on the Evolution of Xenobiotic Metabolism in Humans

Variation in genes involved in the absorption, distribution, metabolism, and excretion of drugs (ADME) can influence individual response to a therapeutic treatment. The study of ADME genetic diversity in human populations has led to evolutionary hypotheses of adaptation to distinct chemical environments. Population differentiation in measured drug metabolism phenotypes is, however, scarcely documented, often indirectly estimated via genotype-predicted phenotypes. We administered seven probe compounds devised to target six cytochrome P450 enzymes and the P-glycoprotein (P-gp) activity to assess phenotypic variation in four populations along a latitudinal transect spanning over Africa, the Middle East, and Europe (349 healthy Ethiopian, Omani, Greek, and Czech volunteers). We demonstrate significant population differentiation for all phenotypes except the one measuring CYP2D6 activity. Genome-wide association studies (GWAS) evidenced that the variability of phenotypes measuring CYP2B6, CYP2C9, CYP2C19, and CYP2D6 activity was associated with genetic variants linked to the corresponding encoding genes, and additional genes for the latter three. Instead, GWAS did not indicate any association between genetic diversity and the phenotypes measuring CYP1A2, CYP3A4, and P-gp activity. Genome scans of selection highlighted multiple candidate regions, a few of which included ADME genes, but none overlapped with the GWAS candidates. Our results suggest that different mechanisms have been shaping the evolution of these phenotypes, including phenotypic plasticity, and possibly some form of balancing selection. We discuss how these contrasting results highlight the diverse evolutionary trajectories of ADME genes and proteins, consistent with the wide spectrum of both endogenous and exogenous molecules that are their substrates.

Relevance of CYP2D6 Gene Variants in Population Genetic Differentiation

A significant portion of the variability in complex features, such as drug response, is likely caused by human genetic diversity. One of the highly polymorphic pharmacogenes is CYP2D6, encoding an enzyme involved in the metabolism of about 25% of commonly prescribed drugs. In a directed search of the 1000 Genomes Phase III variation data, 86 single nucleotide polymorphisms (SNPs) in the CYP2D6 gene were extracted from the genotypes of 2504 individuals from 26 populations, and then used to reconstruct haplotypes. Analyses were performed using Haploview, Phase, and Arlequin softwares. Haplotype and nucleotide diversity were high in all populations, but highest in populations of African ancestry. Pairwise F_ST showed significant results for eleven SNPs, six of which were characteristic of African populations, while four SNPs were most common in East Asian populations. A principal component analysis of CYP2D6 haplotypes showed that African populations form one cluster, Asian populations form another cluster with East and South Asian populations separated, while European populations form the third cluster. Linkage disequilibrium showed that all African populations have three or more haplotype blocks within the CYP2D6 gene, while other world populations have one, except for Chinese Dai and Punjabi in Pakistan populations, which have two.

Genetic Polymorphism of Drug Metabolic Gene CYPs, VKORC1, NAT2, DPYD and CHST3 of Five Ethnic Minorities in Heilongjiang Province, Northeast China

Introduction

Genetic variability in genes encoding drug-metabolizing enzymes may contribute to the heterogeneity of drug responses in different populations. Extensive research in pharmacogenomics in major populations around the world provides us with a great deal of information about drug-related genetic polymorphisms.

Objective

The purpose of this study was to detect the genetic variation of drug-metabolism-related genes in the five ethnic minorities Daur, Hezhen, Ewenki, Mongolian and Manchu in China, and to analyze the distribution differences among ethnic groups.

Methods

We genotyped 32 SNPs of drug metabolism genes in 882 healthy Chinese volunteers from five ethnic groups. The genotype frequency and allele frequency of the five ethnic groups were calculated, and the different variants among the five ethnic groups were compared by chi-square test. Genetic parameters were analyzed using Popgene software. The genetic structure of five ethnic minorities was analyzed by principal component analysis, and compared with 26 populations.

Results

We found that SNPs of genes related to drug metabolism existed diversity in different populations. Among them, rs8192766 and rs9419082 in CYP2E1 showed statistical differences between Daur and Manchu, and NAT2 rs1801280 showed statistical differences between Hezhen and Mongolian. In addition, the five populations we studied had the smallest differences with EAS populations. There was haplotype diversity in CHST3, VKORC1, CYP1A2 and CYP2E1 genes in the five ethnic minorities, and these haplotype polymorphisms were related to the use of corresponding drug doses. Cluster analysis shows that the five ethnic minorities in Heilongjiang Province are clustered together with the EAS populations.

Conclusion

These results suggest that understanding the diversity of drug-related genetic markers is critical for individualized drug gene therapy programs in ethnic minorities in China as well as in populations highly mixed with these ethnic groups.

The influence of evolutionary history on human health and disease

Nearly all genetic variants that influence disease risk have human-specific origins; however, the systems they influence have ancient roots that often trace back to evolutionary events long before the origin of humans. Here, we review how advances in our understanding of the genetic architectures of diseases, recent human evolution and deep evolutionary history can help explain how and why humans in modern environments become ill. Human populations exhibit differences in the prevalence of many common and rare genetic diseases. These differences are largely the result of the diverse environmental, cultural, demographic and genetic histories of modern human populations. Synthesizing our growing knowledge of evolutionary history with genetic medicine, while accounting for environmental and social factors, will help to achieve the promise of personalized genomics and realize the potential hidden in an individual's DNA sequence to guide clinical decisions. In short, precision medicine is fundamentally evolutionary medicine, and integration of evolutionary perspectives into the clinic will support the realization of its full potential.

Genetic ancestry plays a central role in population pharmacogenomics

Recent studies have pointed out the essential role of genetic ancestry in population pharmacogenetics. In this study, we analyzed the whole-genome sequencing data from The 1000 Genomes Project (Phase 3) and the pharmacogenetic information from Drug Bank, PharmGKB, PharmaADME, and Biotransformation. Here we show that ancestry-informative markers are enriched in pharmacogenetic loci, suggesting that trans-ancestry differentiation must be carefully considered in population pharmacogenetics studies. Ancestry-informative pharmacogenetic loci are located in both protein-coding and non-protein-coding regions, illustrating that a whole-genome analysis is necessary for an unbiased examination over pharmacogenetic loci. Finally, those ancestry-informative pharmacogenetic loci that target multiple drugs are often a functional variant, which reflects their importance in biological functions and pathways. In summary, we develop an efficient algorithm for an ultrahigh-dimensional principal component analysis. We create genetic catalogs of ancestry-informative markers and genes. We explore pharmacogenetic patterns and establish a high-accuracy prediction panel of genetic ancestry. Moreover, we construct a genetic ancestry pharmacogenomic database Genetic Ancestry PhD (http://hcyang.stat.sinica.edu.tw/databases/genetic_ancestry_phd/).

Hsin-Chou Yang et al. examine population structure in several genomic databases and identify that pharmacogenetic loci are enriched for markers of genetic ancestry. Their results suggest that genetic ancestry must be carefully considered in population pharmacogenetics studies.

Brazilian cohort and genes encoding for drug-metabolizing enzymes and drug transporters

Background & aim: Genetic variability in drug absorption, distribution, metabolism and excretion (ADME) genes contributes to the high heterogeneity of drug responses. The present study investigated polymorphisms of ADME genes frequencies and compared the findings with populations from other continents, available in the 1000 Genome Project (1 KGP) and the Exome Aggregation Consortium (ExAC) databases. Methodology & results: We conducted a study of 100 patients in Brazil and a total of 2003 SNPs were evaluated by targeted next-generation sequencing in 148 genes, including Phase I enzymes (n = 50), Phase II enzymes (n = 38) and drug transporters (n = 60). Overall, the distribution of minor allele frequency (MAF) suggests that the distribution of 2003 SNPs is similar between Brazilian cohort, 1 KGP and ExAC; however, we found moderate SNP allele-frequency divergence between Brazilian cohort and both 1000 KGP and ExAC. These differences were observed in several relevant genes including CYP3A4, NAT2 and SLCO1B1. Conclusion: We concluded that the Brazilian population needs clinical assessment of drug treatment based on individual genotype rather than ethnicity.

Beyond drugs: the evolution of genes involved in human response to medications

The genetic variation of our species reflects human demographic history and adaptation to diverse local environments. Part of this genetic variation affects individual responses to exogenous substances, such as food, pollutants and drugs, and plays an important role in drug efficacy and safety. This review provides a synthesis of the evolution of loci implicated in human pharmacological response and metabolism, interpreted within the theoretical framework of population genetics and molecular evolution. In particular, I review and discuss key evolutionary aspects of different pharmacogenes in humans and other species, such as the relationship between the type of substrates and rate of evolution; the selective pressure exerted by landscape variables or dietary habits; expected and observed patterns of rare genetic variation. Finally, I discuss how this knowledge can be translated directly or after the implementation of specific studies, into practical guidelines.

Characterization of ADME genes variation in Roma and 20 populations worldwide

The products of the polymorphic ADME genes are involved in Absorption, Distribution, Metabolism, and Excretion of drugs. The pharmacogenetic data have been studied extensively due to their clinical importance in the appropriate drug prescription, but such data from the isolated populations are rather scarce. We analyzed the distribution of 95 polymorphisms in 31 core ADME genes in 20 populations worldwide and in newly genotyped samples from the Roma (Gypsy) population living in Croatia. Global distribution of ADME core gene loci differentiated three major clusters; (1) African, (2) East Asian, and (3) joint European, South Asian and South American cluster. The SLCO1B3 (rs4149117) and CYP3A4 (rs2242480) genes differentiated at the highest level the African group of populations, while NAT2 gene loci (rs1208, rs1801280, and rs1799929) and VKORC1 (rs9923231) differentiated East Asian populations. The VKORC1 rs9923231 was among the investigated loci the one with the largest global minor allele frequency (MAF) range; its MAF ranged from 0.027 in Nigeria to 0.924 in Han Chinese. The distribution of the investigated gene loci positions Roma population within the joined European and South Asian clusters, suggesting that their ADME gene pool is a combination of ancestral (Indian) and more recent (European) surrounding, as it was already implied by other genetic markers. However, when compared to the populations worldwide, the Croatian Roma have extreme MAF values in 10 out of the 95 investigated ADME core gene loci. Among loci which have extraordinary MAFs in Roma population two have strong proof of clinical importance: rs1799853 (CYP2C9) for warfarin dosage, and rs12248560 (CYP2C19) for clopidogrel dosage, efficacy and toxicity. This finding confirms the importance of taking the Roma as well as the other isolated populations`genetic profiles into account in pharmaco-therapeutic practice.

Interleukin-4 and STAT6 promoter polymorphisms but not interleukin-10 or 13 are essential for schistosomiasis and associated disease burden among Nigerian children

Schistosomiasis is endemic in many parts of rural Africa, with previous reports showing interleukin-13 polymorphisms as drivers of infectivity and disease severity in West Africa while IL-13/IL-4 polymorphisms contributes to patterns of reinfection in East Africa. We have shown that there is a genetic delineation in susceptibility to and severity of infectious diseases in Africa, in addition to sub-continental differences in disease pattern. Therefore, which immunoregulatory biomarkers are essential in driving S. haematobium infection or regulate disease burden among Nigerian school children? One hundred and thirty one age and sex-matched schistosome-infected children and 275 uninfected controls, of same ethnicity, recruited from southwestern Nigeria, were screened for variability of cytokine genes, IL-10 (rs1800872), IL-13 (rs7719175), IL-4 (rs2243250) and STAT6 (rs3024974), utilizing a polymerase chain reaction-restriction fragment length polymorphism assay. We found no difference in genotypic or allelic frequencies of IL-10 and IL-13 promoter polymorphisms alone or in association with disease. Contrariwise, we report significant differences in the frequencies of IL-4 and STAT6 variants between groups. For IL-4, the rs2243250 T/T variant was significantly different for genotypes (71.6% versus 51.2%; p < .0004) and alleles (82.6% versus 71.1%; p < .001) between disease and control groups respectively. For STAT6 (rs3024974), the frequencies of genotypes C/C and C/T are 75.4% and 24.6%, both showing an association with disease; none of the infected subjects had the T/T variant. Despite minor differences in disease covariates, we found no association between IL-4 and STAT6 variants with age, gender or anemia. However, mean egg count (indicative of disease burden), was regulated based on IL-4 variants, with highest burden in infected subjects with rs2243250 T/T variant (mean egg count: 207.5 eggs/10 ml of urine) versus rs2243250 C/T heterozygotes (mean egg count: 84.3 eggs/10 ml of urine) versus rs2243250 C/C (mean egg count: 127.9 eggs/10 ml of urine). Comparing rs2243250 C/T versus rs2243250 T/T (p < .008) or rs2243250 C/C + C/T versus rs2243250 T/T (p < .016) reveals an association with disease burden. We conclude that the IL-4 promoter gene is a susceptibility factor for schistosomiasis, and essential to regulate disease burden, with worse disease among carriers of the rs2243250 T/T variant. The absence of the STAT6, rs3024974T/T variant among infected subjects reveal the necessity of the STAT6 promoter gene in driving susceptibility to schistosomiasis in Nigeria.

Interethnic DNA methylation difference and its implications in pharmacoepigenetics

AIM

This is the first systematic study to examine the population differentiation effect of DNA methylation on the treatment response and drug absorption, distribution, metabolism and excretion in multiple tissue types and cancer types.

MATERIALS & METHODS

We analyzed the whole methylome and transcriptome data of primary tumor tissues of four cancer types (breast, colon, head & neck and uterine corpus) and lymphoblastoid cell lines for African and European ancestry populations.

RESULTS

Ethnicity-associated CpG sites exhibited similar methylation patterns in the two studied populations, but the patterns differed between tumor tissues and lymphoblastoid cell lines. Ethnicity-associated CpG sites may have triggered gene expression, influenced drug absorption, distribution, metabolism and excretion, and showed tumor-specific patterns of methylation and gene regulation.

CONCLUSION

Ethnicity should be carefully accounted for in future pharmacoepigenetics research.

A European Spectrum of Pharmacogenomic Biomarkers: Implications for Clinical Pharmacogenomics

Pharmacogenomics aims to correlate inter-individual differences of drug efficacy and/or toxicity with the underlying genetic composition, particularly in genes encoding for protein factors and enzymes involved in drug metabolism and transport. In several European populations, particularly in countries with lower income, information related to the prevalence of pharmacogenomic biomarkers is incomplete or lacking. Here, we have implemented the microattribution approach to assess the pharmacogenomic biomarkers allelic spectrum in 18 European populations, mostly from developing European countries, by analyzing 1,931 pharmacogenomics biomarkers in 231 genes. Our data show significant inter-population pharmacogenomic biomarker allele frequency differences, particularly in 7 clinically actionable pharmacogenomic biomarkers in 7 European populations, affecting drug efficacy and/or toxicity of 51 medication treatment modalities. These data also reflect on the differences observed in the prevalence of high-risk genotypes in these populations, as far as common markers in the CYP2C9, CYP2C19, CYP3A5, VKORC1, SLCO1B1 and TPMT pharmacogenes are concerned. Also, our data demonstrate notable differences in predicted genotype-based warfarin dosing among these populations. Our findings can be exploited not only to develop guidelines for medical prioritization, but most importantly to facilitate integration of pharmacogenomics and to support pre-emptive pharmacogenomic testing. This may subsequently contribute towards significant cost-savings in the overall healthcare expenditure in the participating countries, where pharmacogenomics implementation proves to be cost-effective.

Genetic diversity of disease-associated loci in Turkish population

Many consortia and international projects have investigated the human genetic variation of a large number of ethno-geographic groups. However, populations with peculiar genetic features, such as the Turkish population, are still absent in publically available datasets. To explore the genetic predisposition to health-related traits of the Turkish population, we analyzed 34 genes associated with different health-related traits (for example, lipid metabolism, cardio-vascular diseases, hormone metabolism, cellular detoxification, aging and energy metabolism). We observed relevant differences between the Turkish population and populations with non-European ancestries (that is, Africa and East Asia) in some of the investigated genes (that is, AGT, APOE, CYP1B1, GNB3, IL10, IL6, LIPC and PON1). As most complex traits are highly polygenic, we developed polygenic scores associated with different health-related traits to explore the genetic diversity of the Turkish population with respect to other human groups. This approach showed significant differences between the Turkish population and populations with non-European ancestries, as well as between Turkish and Northern European individuals. This last finding is in agreement with the genetic structure of European and Middle East populations, and may also agree with epidemiological evidences about the health disparities of Turkish communities in Northern European countries.

Genetic architectures of ADME genes in five Eurasian admixed populations and implications for drug safety and efficacy

BACKGROUND

Drug absorption, distribution, metabolism and excretion (ADME) contribute to the high heterogeneity of drug responses in humans. However, the same standard for drug dosage has been applied to all populations in China although genetic differences in ADME genes are expected to exist in different ethnic groups. In particular, the ethnic minorities in northwestern China with substantial ancestry contribution from Western Eurasian people might violate such a single unified standard.

METHODS

In this study, we used Affymetrix SNP Array 6.0 to investigate the genetic diversity of 282 ADME genes in five northwestern Chinese minority populations, namely, Tajik, Uyghur, Kazakh, Kirgiz and Hui, and attempted to identify the highly differential SNPs and haplotypes and further explore their clinical implications.

RESULTS

We found that genetic diversity of many ADME genes in the five minority groups was substantially different from those in the Han Chinese population. For instance, we identified 10 functional SNPs with substantial allele frequency differences, 14 functional SNPs with highly different heterozygous states and eight genes with significant haplotype differences between these admixed minority populations and the Han Chinese population. We further confirmed that these differences mainly resulted from the European gene flow, that is, this gene flow increased the genetic diversity in the admixed populations.

CONCLUSIONS

These results suggest that the ADME genes vary substantially among different Chinese ethnic groups. We suggest it could cause potential clinical risk if the same dosage of substances (eg, antitumour drugs) is used without considering population stratification.

Applying genome-wide gene-based expression quantitative trait locus mapping to study population ancestry and pharmacogenetics

BACKGROUND

Gene-based analysis has become popular in genomic research because of its appealing biological and statistical properties compared with those of a single-locus analysis. However, only a few, if any, studies have discussed a mapping of expression quantitative trait loci (eQTL) in a gene-based framework. Neither study has discussed ancestry-informative eQTL nor investigated their roles in pharmacogenetics by integrating single nucleotide polymorphism (SNP)-based eQTL (s-eQTL) and gene-based eQTL (g-eQTL).

RESULTS

In this g-eQTL mapping study, the transcript expression levels of genes (transcript-level genes; T-genes) were correlated with the SNPs of genes (sequence-level genes; S-genes) by using a method of gene-based partial least squares (PLS). Ancestry-informative transcripts were identified using a rank-score-based multivariate association test, and ancestry-informative eQTL were identified using Fisher's exact test. Furthermore, key ancestry-predictive eQTL were selected in a flexible discriminant analysis. We analyzed SNPs and gene expression of 210 independent people of African-, Asian- and European-descent. We identified numerous cis- and trans-acting g-eQTL and s-eQTL for each population by using PLS. We observed ancestry information enriched in eQTL. Furthermore, we identified 2 ancestry-informative eQTL associated with adverse drug reactions and/or drug response. Rs1045642, located on MDR1, is an ancestry-informative eQTL (P = 2.13E-13, using Fisher's exact test) associated with adverse drug reactions to amitriptyline and nortriptyline and drug responses to morphine. Rs20455, located in KIF6, is an ancestry-informative eQTL (P = 2.76E-23, using Fisher's exact test) associated with the response to statin drugs (e.g., pravastatin and atorvastatin). The ancestry-informative eQTL of drug biotransformation genes were also observed; cross-population cis-acting expression regulators included SPG7, TAP2, SLC7A7, and CYP4F2. Finally, we also identified key ancestry-predictive eQTL and established classification models with promising training and testing accuracies in separating samples from close populations.

CONCLUSIONS

In summary, we developed a gene-based PLS procedure and a SAS macro for identifying g-eQTL and s-eQTL. We established data archives of eQTL for global populations. The program and data archives are accessible at http://www.stat.sinica.edu.tw/hsinchou/genetics/eQTL/HapMapII.htm. Finally, the results from our investigations regarding the interrelationship between eQTL, ancestry information, and pharmacodynamics provide rich resources for future eQTL studies and practical applications in population genetics and medical genetics.