Genetic polymorphism analysis of cytochrome P4502E1 (CYP2E1) in Chinese Han populations from four different geographic areas of Mainland China

CYP2E1 is one of a superfamily of enzymes that play a central role in activating and detoxifying many xenobiotics and endogenous compounds thought to be involved in the development of several human diseases. Among other factors, individual susceptibility to developing these pathologies relies on genetic polymorphisms, which are related to ethnic differences, since the frequency of mutant genotypes varies in different populations. The aim of this study was to investigate the genetic basis of CYP2E1 polymorphisms in the populations of four different geographical locations of China. Twenty-two different CYP2E1 polymorphisms, including six novel variants in promoter regions and a novel nonsense mutation, were identified. The frequencies of some polymorphisms and genotypes demonstrated significant differences among the four populations. Linkage disequilibrium analysis and tag SNP selection were performed. Haplotypes were analyzed within the selected tag SNPs. Tag SNP selection and haplotype distributions showed differences across the four populations.     

Genetic Variation in TLR Genes in Ugandan and South African Populations and Comparison with HapMap Data

Genetic epidemiological studies of complex diseases often rely on data from the International HapMap Consortium for identification of single nucleotide polymorphisms (SNPs), particularly those that tag haplotypes. However, little is known about the relevance of the African populations used to collect HapMap data for study populations conducted elsewhere in Africa. Toll-like receptor (TLR) genes play a key role in susceptibility to various infectious diseases, including tuberculosis. We conducted full-exon sequencing in samples obtained from Uganda (n = 48) and South Africa (n = 48), in four genes in the TLR pathway: TLR2, TLR4, TLR6, and TIRAP. We identified one novel TIRAP SNP (with minor allele frequency [MAF] 3.2%) and a novel TLR6 SNP (MAF 8%) in the Ugandan population, and a TLR6 SNP that is unique to the South African population (MAF 14%). These SNPs were also not present in the 1000 Genomes data. Genotype and haplotype frequencies and linkage disequilibrium patterns in Uganda and South Africa were similar to African populations in the HapMap datasets. Multidimensional scaling analysis of polymorphisms in all four genes suggested broad overlap of all of the examined African populations. Based on these data, we propose that there is enough similarity among African populations represented in the HapMap database to justify initial SNP selection for genetic epidemiological studies in Uganda and South Africa. We also discovered three novel polymorphisms that appear to be population-specific and would only be detected by sequencing efforts.     

Novel Single Nucleotide Polymorphisms in Interleukin 6 Affect Tacrolimus Metabolism in Liver Transplant Patients

Background

Tacrolimus is the first-line immunosuppressant after organ transplantation. It is mainly metabolized by cytochrome P450, family 3, subfamily A (CYP3A) enzymes, but there are large individual differences in metabolism. Interleukin 6 (IL6) has been shown to cause a pan-suppression of mRNA levels of ten major CYP enzymes in human hepatocyte cultures. IL6 has been shown to provide hepatoprotection in various models of liver injury. Rs1800796 is a locus in the IL6 gene promoter region which regulates cytokine production. We speculated that IL6 rs1800796 polymorphisms may lead to individual differences in tacrolimus metabolism by affecting CYP3A enzymes levels and liver function after liver transplantation.

Methodology/Principal Findings

Ninety-six liver transplant patients receiving tacrolimus were enrolled in the study. Two single nucleotide polymorphisms (SNP), CYP3A5 rs776746 and IL6 rs1800796, were genotyped in both donors and recipients. The effects of SNPs on tacrolimus concentration/dose (C/D ratio) at four weeks after transplantation were studied, as well as the effects of donor IL6 rs1800796 polymorphisms on liver function. Both donor and recipient CYP3A5 rs776746 allele A showed association with lower C/D ratios, while donor IL6 rs1800796 allele G showed an association with higher C/D ratios. Donor CYP3A5 rs776746 allele A, IL6 rs1800796 allele C, and recipient CYP3A5 rs776746 allele A were associated with fast tacrolimus metabolism. With increasing numbers of these alleles, patients were found to have increasingly lower tacrolimus C/D ratios at time points after transplantation. Donor IL6 rs1800796 allele G carriers showed an association with higher glutamic-pyruvic transaminase (GPT) levels.

Conclusions

Combined analysis of donor CYP3A5 rs776746, IL6 rs1800796, and recipient CYP3A5 rs776746 polymorphisms may distinguish tacrolimus metabolism better than CYP3A5 rs776746 alone. IL6 may lead to individual differences in tacrolimus metabolism mainly by affecting liver function.     

Haplotype structure,adaptive history and associations with exploratory behaviour of the DRD4 gene region in four great tit (Parus major) populations

The assessment of genetic architecture and selection history in genes for behavioural traits is fundamental to our understanding of how these traits evolve. The dopamine receptor D4 (DRD4) gene is a prime candidate for explaining genetic variation in novelty seeking behaviour, a commonly assayed personality trait in animals. Previously, we showed that a single nucleotide polymorphism in exon 3 of this gene is associated with exploratory behaviour in at least one of four Western European great tit (Parus major) populations. These heterogeneous association results were explained by potential variable linkage disequilibrium (LD) patterns between this marker and the causal variant or by other genetic or environmental differences among the populations. Different adaptive histories are further hypothesized to have contributed to these population differences. Here, we genotyped 98 polymorphisms of the complete DRD4 gene including the flanking regions for 595 individuals of the four populations. We show that the LD structure, specifically around the original exon 3 SNP is conserved across the four populations and does not explain the heterogeneous association results. Study‐wide significant associations with exploratory behaviour were detected in more than one haplotype block around exon 2, 3 and 4 in two of the four tested populations with different allele effect models. This indicates genetic heterogeneity in the association between multiple DRD4 polymorphisms and exploratory behaviour across populations. The association signals were in or close to regions with signatures of positive selection. We therefore hypothesize that variation in exploratory and other dopamine‐related behaviour evolves locally by occasional adaptive shifts in the frequency of underlying genetic variants.     

Polymorphisms and phenotypic analysis of cytochrome P450 2D6 in the Tibetan population

The cytochrome P450 2D6 enzyme (CYP2D6) metabolizes about 25% of prescribed drugs in the endoplasmic reticulum, and genetic polymorphisms in CYP2D6 can greatly affect its activity and lead to differences among individuals in drug efficacy and adverse drug reactions. To investigate genetic polymorphisms in CYP2D6 among Tibetan Chinese, we directly sequenced the whole gene in 96 unrelated, healthy Tibetans from The Tibet Autonomous Region of China and screened for genetic variants in the promoter, exons, introns, and 3′UTR. We detected fifty-one genetic polymorphisms in CYP2D6, and 16 of them are novel. The allele frequencies of CYP2D6*1, *2, *5, *10, *41, and *49 were 0.25, 0.43, 0.02, 0.29, 0.02, and 0.01, respectively. The frequency of CYP2D6*10, a putative poor-metabolizer allele, was lower in our sample population compared with that in the Han Chinese population (p < 0.001). In addition, haplotype analysis allowed 15 CYP2D6 haplotypes to be classified into three groups. In conclusion, our results provide basic information about CPY2D6 alleles in Tibetans and suggest that the enzymatic activities of CYP2D6 may differ among the diverse ethnic populations of China. Our results provide a basis for safer drug administration and better therapeutic treatment among Tibetans.     

Ecologically and evolutionarily important SNPs identified in natural populations

Evolution by natural selection acts on natural populations amidst migration, gene-by-environmental interactions, constraints, and tradeoffs, which affect the rate and frequency of adaptive change. We asked how many and how rapidly loci change in populations subject to severe, recent environmental changes. To address these questions, we used genomic approaches to identify randomly selected single nucleotide polymorphisms (SNPs) with evolutionarily significant patterns in three natural populations of Fundulus heteroclitus that inhabit and have adapted to highly polluted Superfund sites. Three statistical tests identified 1.4-2.5% of SNPs that were significantly different from the neutral model in each polluted population. These nonneutral patterns in populations adapted to highly polluted environments suggest that these loci or closely linked loci are evolving by natural selection. One SNP identified in all polluted populations using all tests is in the gene for the xenobiotic metabolizing enzyme, cytochrome P4501A (CYP1A), which has been identified previously as being refractory to induction in the three highly polluted populations. Extrapolating across the genome, these data suggest that rapid evolutionary change in natural populations can involve hundreds of loci, a few of which will be shared in independent events.     

A Systematically Combined Genotype and Functional Combination Analysis of CYP2E1, CYP2D6, CYP2C9, CYP2C19 in Different Geographic Areas of Mainland China – A Basis for Personalized Therapy

The cytochrome P450 is the major enzyme involved in drug metabolism. Single CYP genotypes and metabolic phenotypes have been widely studied, but no combination analysis has been conducted in the context of specific populations and geographical areas. This study is the first to systematically analyze the combined genotypes and functional combinations of 400 samples of major CYP genes—CYP2E1, CYP2D6, CYP2C9, and CYP2C19 in four geographical areas of mainland China. 167 different genotype combinations were identified, of which 25 had a greater than 1% frequency in the Chinese Han population. In addition, phenotypes of the four genes for each sample were in line with the predictions of previous studies of the four geographical areas. On the basis of the genotype classification, we were able to produce a systemic functional combinations analysis for the population. 25 of the combinations detected had at least two non-wild phenotypes and four showed a frequency above 1%. A bioinformatics analysis of the relationship between particular drugs and multi-genes was conducted. This is the first systematic study to analyze genotype combinations and functional combinations across whole Chinese population and could make a significant contribution in the field of personalized medicine and therapy.     

TAGster: efficient selection of LD tag SNPs in single or multiple populations

Genetic association studies increasingly rely on the use of linkage disequilibrium (LD) tag SNPs to reduce genotyping costs. We developed a software package TAGster to select, evaluate and visualize LD tag SNPs both for single and multiple populations. We implement several strategies to improve the efficiency of current LD tag SNP selection algorithms: (1) we modify the tag SNP selection procedure of Carlson et al. to improve selection efficiency and further generalize it to multiple populations. (2) We propose a redundant SNP elimination step to speed up the exhaustive tag SNP search algorithm proposed by Qin et al. (3) We present an additional multiple population tag SNP selection algorithm based on the framework of Howie et al., but using our modified exhaustive search procedure. We evaluate these methods using resequenced candidate gene data from the Environmental Genome Project and show improvements in both computational and tagging efficiency. AVAILABILITY: The software Package TAGster is freely available at http://www.niehs.nih.gov/research/resources/software/tagster/     

Somatic cell mapping, polymorphism, and linkage analysis of bovine prolactin-related proteins and placental lactogen.

The bovine prolactin gene family includes novel members expressed in the fetal placenta that are distinct from placental lactogen. In this study, we investigated the genetic organization of four members of this gene family (PRP1, PRP3, PRP6, and PRP10) as well as placental lactogen (PL). Using a bovine-rodent hybrid somatic cell panel, all five genes were assigned to bovine chromosome 23, which contains prolactin and the major histocompatibility group (BOLA). Restriction fragment length polymorphisms were detected by all probes in breeding populations with the restriction enzyme MspI, whereas no polymorphisms were detected with BamHI. EcoRI, HindIII, TaqI, and PstI produced polymorphic fragments with some but not all of the probes tested. A PRP10 polymorphism, which is apparently the result of a insertion/deletion event, detected polymorphism frequency differences between Bos indicus and Bos taurus. No recombinational events were observed with these probes and prolactin using linkage analysis involving 91 American Holsteins. The bovine prolactin gene family was incorporated into a linkage group containing CYP21. Our studies demonstrate that members of the bovine prolactin gene family have a close physical association with each other, and all members demonstrate genetic variability in the breeding population.     

Reduced Levels of DNA Polymorphism and Fixed between-Population Differences in the Centromeric Region of Drosophila Ananassae

We have estimated DNA sequence variation within and between two populations of Drosophila ananassae, using six-cutter restriction site variation at vermilion (v) and furrowed (fw). These two gene regions are located close to the centromere on the left and right X chromosome arms, respectively. In the fw region, no DNA polymorphism was detected within each population. In the v region, average heterozygosity per nucleotide was very low in both populations (pi = 0.0005 in the Burma population, and 0.0009 in the India population). These estimates are significantly lower than those from loci in more distal gene regions. The distribution of DNA polymorphisms between both populations was also striking. At fw, three fixed differences between the Burma and India populations were detected (two restriction site differences and one insertion/deletion of approximately 2 kb). At v, each DNA polymorphism in high frequency in the total sample was nearly fixed in one or the other population, although none of them reached complete fixation. The observed pattern of reduced variation within populations and fixed differences between populations appears to correlate with recombination rate. We conclude that recent hitchhiking associated with directional selection is the best explanation for this pattern. The data indicate that different selective sweeps have occurred in the two populations. The possible role of genetic hitchhiking in rapid population differentiation in gene regions of restricted recombination is discussed.     

Frequencies of four genetic polymorphisms in the CYP1A2 gene in Turkish population

Cytochrome P450 (CYP) 1A2 gene is involved in the metabolic activation of several carcinogens and altered metabolization of some clinically used drugs. We aimed to investigate the distributions of genetic polymorphisms -3860 (G/A)(CYP1A2*1C) and -2467 (T/del)(CYP1A2*1D) in the 5'-flanking region and -739 (T/G)(CYP1A2*1E) and -163(C/A)(CYP1A2*1F) in the first intron of the CYP1A2 gene in 110 unrelated healthy Turkish volunteers by PCR-RFLP technique. The frequencies of each polymorphism in Turkish population were found as 0.04, 0.92, 0.01, 0.27 for CYP1A2*1C, CYP1A2*1D, CYP1A2*1E, CYP1A2*1F, respectively. Compared with other populations, CYP1A2*1D has been found to be significantly increased in Turkish population. On the other hand, in general, the frequencies of the other polymorphisms were concordant with those in the Egyptian and Caucasian populations, and were different from those in the Japanese, Chinese and Ethiopian populations. Our results suggest that due to increased frequency of CYP1A2*1D in Turkish population, functional significance of CYP1A2*1D should be evaluated. It might be screened to determine the relationship between CYP1A2*1D and CYP1A2 related drug metabolisms in associated groups.     

CYP2A6 polymorphism reveals differences in Japan and the existence of a specific variant in Ovambo and Turk populations

CYP2A6 is a polymorphic enzyme, and CYP2A6 genotype has been shown to be associated with smoking habits and lung cancer. We investigated CYP2A6 polymorphism in Japanese from four different geographic areas of Japan and in the Ovambo and Turk populations. Using two polymerase chain reaction restriction fragment length polymorphisms (PCR-RFLPs), we identified the functionally important variants of CYP2A6: *1A, *1B, *1F, *1G, *4A, and *4D. In the Japanese population the highest frequencies of the CYP2A6*1A allele were observed in subjects from the Fukuoka (Kyushu Island) and Ehime (Shikoku Island) prefectures, whereas subjects in Shimane and Tottori (both located on the Japan Sea side of Honshu Island) showed the highest frequencies of the CYP2A6*1B allele. In the Tottori and Shimane groups no subject was homozygous for the CYP2A6*4A allele, a whole gene deletion type that is prevalent among Asians. In the Ovambo and Turk populations the CYP2A6*1A allele was predominant. Furthermore, two alleles undetected in the Japanese were observed in these latter two ethnic groups: CYP2A6*1G was found solely in the Ovambos, and CYP2A6*1F was found solely in the Turks. The present study is the first to show interprefecture differences in CYP2A6 polymorphism in Japanese who live in relatively close but distinct geographic areas; this is also the first study to evaluate CYP2A6 variations among these Japanese and the Ovambo and Turk populations. The distribution results of these alleles could help to define the true significance of CYP2A6 polymorphism as a genetic susceptibility marker in worldwide populations.     

Frequency of single nucleotide polymorphisms of CYP2D6 in the Czech population

CYP2D6 is a member of cytochrome P450 enzymes that metabolise over 25% of commonly used drugs. Genetic polymorphisms can cause insufficient drug efficacy at usually administered doses or can be the cause of adverse drug reaction. CYP2D6 genotyping can be used to predict CYP2D6 phenotype and thereby explain some abnormalities in drug response and thus optimize pharmacotherapy. The aim of this study was to investigate the frequency of functionally important variant alleles of the CYP2D6 gene throughout the Czech population to predict the prevalence of ultra-rapid and poor metabolizer phenotypes. The DNA of 223 unrelated, healthy volunteers was analysed to detect the presence of CYP2D6*6, *5, *4, *3 and gene duplication. The variant allele frequencies in our population were 0.22%, 3.14%, 22.87%, 1.12% and 3.14% for CYP2D6*6, CYP2D6*5, CYP2D6*4, CYP2D6*3 and CYP2D6*MxN, respectively. Fifteen subjects carried two variant alleles leading to predicted poor type of metabolism, 84 subjects were heterozygous extensive metabolizers (het-EM). The full-text contains detailed comparison with European white populations. The distribution of variant alleles complies with the Hardy-Weinberg equilibrium. The frequencies of functional variant alleles of CYP2D6 in Czech population are in concordance with other Caucasian populations.     

Data-mining methods as useful tools for predicting individual drug response: application to CYP2D6 data

OBJECTIVES: Selecting a maximally informative subset of polymorphisms to predict a clinical outcome, such as drug response, requires appropriate search methods due to the increased dimensionality associated with looking at multiple genotypes. In this study, we investigated the ability of several pattern recognition methods to identify the most informative markers in the CYP2D6 gene for the prediction of CYP2D6 metabolizer status. METHODS: Four data-mining tools were explored: decision trees, random forests, artificial neural networks, and the multifactor dimensionality reduction (MDR) method. Marker selection was performed separately in eight population samples of different ethnic origin to evaluate to what extent the most informative markers differ across ethnic groups. RESULTS: Our results show that the number of polymorphisms required to predict CYP2D6 metabolic phenotype with a high accuracy can be dramatically reduced owing to the strong haplotype block structure observed at CYP2D6. MDR and neural networks provided nearly identical results and performed the best. CONCLUSION: Data-mining methods, such as MDR and neural networks, appear as promising tools to improve the efficiency of genotyping tests in pharmacogenetics with the ultimate goal of pre-screening patients for individual therapy selection with minimum genotyping effort.     

DNA haplotype-dependent differences in the amino acid sequence of debrisoquine 4-hydroxylase (CYP2D6): evidence for two major allozymes in extensive metabolisers

The molecular basis for DNA haplotype-dependent debrisoquine 4-hydroxylase (CYP2D6) expression was explored by sequencing all of the nine exons of the CYP2D6 gene. Two distinct exon sequence frameworks of the CYP2D6 gene were found, each associated with specific BamHI-defined DNA haplotypes of the CYP2D cluster. They corresponded to Arg²⁹⁶/Cys²⁹⁶ and Ser⁴⁸⁶/Thr⁴⁸⁶ amino acid polymorphisms in the CYP2D6 enzyme, and occurred in almost equal frequency among the Caucasians examined. These two major allozymes with amino acid differences in the presumed substrate recognition region and in the vicinity of the heme binding site could be the source of the observed DNA haplotype-dependent variation in phenotypic expression.     

Development of Genotyping Methods for Single Nucleotide Polymorphism in the Human Pancreatic Ribonuclease Gene (RNASE1) and Their Application to Population Studies

Two potential single nucleotide polymorphisms [SNPs; rs1804215 (G979T) and rs11545379 (G1169T)] have been identified in the human pancreatic ribonuclease, RNase 1, gene (RNASE1) that could give rise to an amino acid substitution in the protein, but relevant population data are not available. We have developed genotyping methods for each SNP using the mismatched PCR-restriction fragment length polymorphism technique. These methods are advantageous in comparison with other SNP genotyping methods because they are technically simpler and do not require specialized instruments. We applied these genotyping methods to examine the genotype distribution of each SNP in four populations, including Japanese populations living in two prefectures, an Ovambo population, and a Turkish population. In all the populations studied, however, only a single genotype for each SNP was found. Therefore, irrespective of differences in ethnic groups, RNASE1 might show markedly low heterogeneity in its genetic structure with regard to these SNPs.     

Comparative study of the linkage disequilibrium of an ENCODE region, chromosome 7p15, in Korean, Japanese, and Han Chinese samples

The extent and pattern of linkage disequilibrium (LD) in the human genome provide important information for disease gene mapping. Previous studies have shown that LDs vary depending on chromosomal regions and populations. As the Asian samples of the International HapMap Project consisted of Japanese and Chinese populations, it was of interest whether we could use the HapMap data as a reference to carry out association studies of common complex diseases in a closely related population, such as Koreans. We have compared the LD and recombination patterns defined by single-nucleotide polymorphisms (SNPs) in ENCODE region ENm010, chromosome 7p15.2, in Korean, Japanese, and Chinese samples and further tested the robustness of tagSNPs among the Asian samples. We genotyped 792 SNPs in 500 kb (chromosome 7: 26699793-27199792, NCBI build 34) from 90 unrelated Koreans by fluorescence polarization detection and compared the data with Asian data from the HapMap project. Despite some differences in the position of high LD region boundaries, the overall patterns of LD were remarkably similar across the three samples, reflecting strong genetic affinities among them. Furthermore, the haplotype tag SNP transferability across the three samples was greater than 90%. Our results support the initial suggestion that the populations genotyped in the HapMap project might serve as reference populations for the selection of tagSNPs in association studies.     

Genome-wide selection of tag SNPs using multiple-marker correlation

MOTIVATIONS: The tag SNP approach is a valuable tool in whole genome association studies, and a variety of algorithms have been proposed to identify the optimal tag SNP set. Currently, most tag SNP selection is based on two-marker (pairwise) linkage disequilibrium (LD). Recent literature has shown that multiple-marker LD also contains useful information that can further increase the genetic coverage of the tag SNP set. Thus, tag SNP selection methods that incorporate multiple-marker LD are expected to have advantages in terms of genetic coverage and statistical power. RESULTS: We propose a novel algorithm to select tag SNPs in an iterative procedure. In each iteration loop, the SNP that captures the most neighboring SNPs (through pair-wise and multiple-marker LD) is selected as a tag SNP. We optimize the algorithm and computer program to make our approach feasible on today's typical workstations. Benchmarked using HapMap release 21, our algorithm outperforms standard pair-wise LD approach in several aspects. (i) It improves genetic coverage (e.g. by 7.2% for 200 K tag SNPs in HapMap CEU) compared to its conventional pair-wise counterpart, when conditioning on a fixed tag SNP number. (ii) It saves genotyping costs substantially when conditioning on fixed genetic coverage (e.g. 34.1% saving in HapMap CEU at 90% coverage). (iii) Tag SNPs identified using multiple-marker LD have good portability across closely related ethnic groups and (iv) show higher statistical power in association tests than those selected using conventional methods. AVAILABILITY: A computer software suite, multiTag, has been developed based on this novel algorithm. The program is freely available by written request to the author at ke_hao@merck.com     

Systematic functional characterization of cytochrome P450 2E1 promoter variants in the Chinese Han population

CYP2E1 promoter polymorphisms can lead to significant interindividual differences in expression of CYP2E1. Using a database of CYP2E1 gene polymorphisms established in 2010, our study aimed to functionally characterize the single nucleotide polymorphisms (SNPs) of the promoter region and corresponding haplotypes in the Chinese Han population. Six novel SNPs and seven haplotypes with a frequency equal to or greater than 0.01 were constructed on a luciferase reporter system on the basis of site-directed mutagenesis. Dual luciferase reporter systems were used to analyze regulatory activity. The constructs including single novel SNP mutations exhibited insignificant change in luciferase activity, whereas, the activity produced by Haplo1(GTTGCTATAT), Haplo2 (CTTGCTATAT) and Haplo7 (GAGCTCACAT), containing a -333T>A polymorphism was significantly greater than for the wild type in Hep G2 cells (p<0.05), being 1.5-, 2.0- and 1.4- times greater respectively. These findings suggest the possibility of significant clinical prediction of adverse drug reaction and the facilitation of personalized medicine.