Visual Overview
Abstract
Unraveling the molecular mechanisms by which genetic variants of cytochrome P450 2A6 lead to different metabolic phenotypes remains a long-standing but important challenge. CYP2A6 is an enzyme involved in the metabolism of several clinical drugs as well as the metabolic activation of carcinogenic nitrosamines. Herein, CYP2A6 genotypes and phenotypes, as indicated by protein content [by liquid chromatography–mass spectrometry (MS)/MS] and metabolic activities [Vmax, clearance (CL)], were determined for 90 human liver samples. We determined the median, range, and interindividual and intraindividual variation of CYP2A6 content and activity at the microsomal, liver tissue, and whole liver level and predicted hepatic in vivo clearance by in vitro–in vivo extrapolation based on CYP2A6-mediated coumarin metabolism by each CYP2A6 genotype. These results reveal how different CYP2A6 genotypes yield different phenotypic traits in protein content and enzyme activity. For relative Vmax, CL, and protein content, the intraindividual percentage coefficients of variation (ICVs) were 41.0% (18.8%–125.1%), 28.5% (2.39%–133.5%), and 27.8% (2.68%–88.0%), respectively. The high ICVs implied large intraindividual variation at different levels, sometimes in a genotype-dependent manner. Intergenotype analysis revealed that the CYP2A6*4 allele demonstrated the most obvious effect on phenotypic outcomes, both in protein content and in metabolic activity. Indeed, decreased CYP2A6 protein content with the CYP2A6*4 genotype might explain the decreased metabolic activity from the molecular to the organismal level. These findings may allow useful predictions for CYP2A6-mediated drug metabolism on an individual patient basis in accord with the goal of achieving personalized medicine.
SIGNIFICANCE STATEMENT We provide the median, range, and interindividual and intraindividual variation in CYP2A6 content at the microsomal, liver tissue, and whole liver level by liquid chromatography–mass spectrometry (MS)/MS as well as activities at the protein, microsomal, liver tissue, and whole liver level both in vitro and at the organismal level based on CYP2A6-mediated coumarin metabolism with each CYP2A6 genotype, thereby allowing us to elucidate how different CYP2A6 genotypes yield differing phenotypic traits (protein content and enzyme activity), facilitating the development of personalized medicine.
Introduction
A recent study (Tennessen et al., 2012) yielding the sequence of 15,585 human protein-coding genes showed that 2.3% of the 13,595 single nucleotide variants involving ∼313 genes were predicted to impact protein function. Cytochrome P450 enzymes are characterized by abundant gene polymorphisms and great interindividual variations in drug-metabolizing activity (Zanger and Schwab, 2013). Our previous studies carried out with human liver samples derived from both individuals who were normal or had hepatocellular carcinoma indicated that various genetic variants of CYP450 (Gao et al., 2016b; Fang et al., 2018a,b) combined with other factors, including CYP450 reductase (Tian et al., 2016; Fang et al., 2018b) and CYPb5 (Zhang et al., 2015a), generated significant variation in metabolic phenotypes (Zhang et al., 2016; Gao et al., 2017c).
More importantly, differences in CYP activity (Gao et al., 2016a; Zhou et al., 2016) are associated with increased risk for diseases like hepatofibrosis (Gao et al., 2017b) and hepatocarcinogenesis (Gao et al., 2018a,b). Thus, studies on CYP genotype-phenotype associations may have important disease implications. However, with regard to metabolic phenotypes for CYPs, research in this area has focused primarily on the effect of CYP genetic variants on enzyme-catalyzed drug metabolism in vitro (Palma et al., 2010; Mürdter et al., 2012; Dai et al., 2014; Werk and Cascorbi, 2014) or on drug clearance in vivo (Mega et al., 2009; Zordoky and El-Kadi, 2010; Zhu et al., 2011; Shuker et al., 2016). It has long been known that genetic effects on physiologic phenotypes often depend on the cellular and environmental context, such as culture media, cell type, or tissue. For this reason, a more comprehensive process for CYP genotype-phenotype association is needed.
In the past several years, CYP2A6 has gained importance in numerous studies on the more than 45 known allelic variants (http://www.imm.ki.se/CYPalleles/cyp2a6.htm) and their roles in the metabolism of coumarin and a number of pharmaceuticals. Along with these pharmaceuticals, toxins and procarcinogens, such as aflatoxin B1 and tobacco-specific nitrosamines, may exhibit differences in blood concentration and duration of effect as a result of CYP2A6 polymorphisms (Sadeque et al., 1997; White et al., 2001; Xu et al., 2002) that are metabolically activated by CYP2A6. Development of a comprehensive genotype-phenotype model for CYP2A6 has the potential to improve our understanding of phenotypic variation and enhance the prediction of clinical outcomes for drugs metabolized by this cytochrome P450.
Clinically, studies on CYP2A6 genotype-phenotype associations may have important clinical implications because of CYP2A6’s major role in the metabolism of nicotine. However, studies on CYP2A6 genotypes determined in vitro (Fukami et al., 2005; Wang et al., 2006) (e.g., by site-directed mutagenesis and liver microsomes) and/or in vivo (e.g., by phenotyping using substrates) are not always consistent. Take CYP2A6*1B as example: It was characterized with either elevated catalytic activity (Yoshida et al., 2003; Wang et al., 2006) or normal activity (Ujjin et al., 2002; Peamkrasatam et al., 2006; Wang et al., 2011), making it difficult to determine the genotype-phenotype association and predict phenotypic outcomes.
One recent report focused on the CYP2A6 genotype-phenotype correlation and showed large interindividual variation in phenotypes (Tanner et al., 2017). Unfortunately, insufficient background information for more than half of the donors made accurate and rigorous evaluation of the genotype-to-phenotype association difficult. Moreover, phenotypes including protein content and activity measured in that study were only determined at one level in vitro. The extrapolation from phenotype in vitro to phenotype in vivo is complex and may involve multiple intermediates. Simple genotype-to-phenotype prediction may not be sufficient to accurately determine the phenotype. As a result, efforts to identify the influential intermediate steps in the process from genotype to ultimate phenotype for CYP2A6 are still an urgent need.
With regard to the characterization of CYP2A6 phenotypes, as mentioned in a previous study (Gao et al., 2017a), they can be divided into content phenotypes and metabolic phenotypes, as evidenced by CYP2A6 protein expression and activity from the molecular level to the cellular level and then to organismal level.
We generated a genotype-phenotype model for CYP2A6 based on three gene polymorphisms of CYP2A6 (CYP2A6*1B, *4, and *9) with frequencies greater than 1% in Chinese subjects. We determined protein expression and metabolic activity from the molecular to the cellular (microsomal, liver tissue, and liver) and organismal level and presented the intraindividual and interindividual variation in metabolic activity for CYP2A6 with different genotypes. Consideration of differences in enzyme activity at the molecular level and the effect on clearance at the organismal level may provide insight into how CYP2A6 molecular genetic variation influences phenotypic variation.
Material and Methods
Human Liver Samples
Normal human liver tissues were obtained from 90 Chinese patients undergoing liver surgical resection at Henan Provincial People’s Hospital or Henan Provincial Tumor Hospital during 2012 and 2014. This research was performed in accordance with the Declaration of Helsinki. All experiments were approved by the Medical Ethics Committee of Zhengzhou University, and written informed consent was obtained from all patients. All liver samples were obtained from patients who were undergoing hepatic surgery and free of infectious disease, with medical diagnoses mostly of cavernous hemangioma of liver (61 cases); a few had metastatic cancers (seven cases), cholelithiasis (eight cases), gallbladder cancer (four cases), hepatic cholangiocarcinoma (six cases), or hepatocellular carcinoma (four cases). Detailed demographic information for the liver samples are presented in Supplemental Table 1. For patients who were tumor-bearing, normal liver samples were obtained 2 cm distant from the tumor tissue. Liver function tests, histopathological analysis, and imaging by ultrasonography or computed tomography were used to confirm that only liver tissue with normal liver function and normal histologic appearance was collected. Liver samples were frozen immediately after removal and stored in liquid nitrogen until use. There was no history of exposure to known CYP-inducing or CYP-inhibiting agents, and only routine anesthetics were used with the patients. Detailed donor characteristics of the human liver samples and human liver microsomes (HLMs) and the method of preparation by differential centrifugation are described in a previous study (Zhang et al., 2015b). Total HLM protein concentration was measured by the Bradford method (Bradford, 1976).
Genotyping of CYP2A6
Genomic DNA was isolated from human liver tissues using a genomic DNA purification kit (Beijing ComWin Biotech Co., Ltd., China). A two-step polymerase chain reaction method (Oscarson et al., 1999) was used to detect CYP2A6*1B and CYP2A6*4. Another single nucleotide polymorphism CYP2A6*9 (rs28399433) was genotyped by polymerase chain reaction sequencing. Genotyping errors were detected by regenotyping with a subsample, and reproducibility was routinely greater than 99%.
Quantification of CYP2A6 Protein
Protein quantification of CYP2A6 was performed by nano–liquid chromatography (LC)–mass spectrometry (MS)/MS as previously reported (Zhang et al., 2016). A QconCAT protein consisting of 57 stable isotope-labeled peptides from 21 drug-metabolizing enzymes (including CYP2A6) in which two or three peptides were selected for each targeted protein was employed to quantify these proteins in HLMs (Li et al., 2015). The concentration of CYP2A6 protein was determined by nano-LC–multiple reaction monitoring MS using an easy nano-LC (Thermo Fisher Scientific Inc., Waltham, MA) coupled to a TSQVantage Triple Quadrupole mass spectrometer (Thermo Fisher Scientific Inc.).
Determination of Phenotype at Different Levels
Molecular phenotype is defined as protein content (Content-P) and activity [Vmax-P, clearance (CL)-P] at the protein level. Cellular phenotype can be divided into microsomal, tissue, and organ levels; accordingly, protein expression at the cellular level in the liver is defined as Content-M, Content-LT, and Content-L, and activity is defined as Vmax-M, Vmax-LT, and Vmax-L and CL-M, CL-LT, and CL-L, respectively. For organismal phenotype, this refers to drug clearance from the body at the whole organism level (in vivo).
Microsomal Phenotype.
The subcellular phenotype, termed the microsomal phenotype, is defined as protein content (CYP-M), Vmax (Vmax-M), and clearance (CL-M). Vmax-M and CL-M were determined by measuring the rate of coumarin 7-hydroxylation by high-performance liquid chromatography with eight coumarin concentrations (0.156–20 µM) according to a previously described method (Zhang et al., 2016). Incubation mixtures contained HLMs, 100 mM phosphate buffer (pH 7.4) with 1 mM NADPH, and substrate. Optimal protein concentration was 0.3 mg/ml protein, and incubation time was 30 minutes. Perchloric acid (10 µl) was added to terminate the reaction. The metabolite was identified by high-performance LC–UV. All experiments included two replicates. The Michaelis-Menten constant (Km) and maximum reaction rate (Vmax-M) of CYP2A6 were determined by nonlinear regression analysis using GraphPad Prism 6.0, and CL-M was calculated based on the ratio of Vmax-M to Km.(1)
Molecular Phenotype.
The metabolic phenotype at the molecular level is defined as Vmax (Vmax-P) and clearance (CL-P) of drug per picomole of CYP2A6. The Vmax-P and CL-P for each individual were obtained by dividing each individual microsomal phenotype parameter by the corresponding individual CYP2A6 protein content (Content-M).(2)
Tissue Phenotype.
The tissue phenotype is defined as protein content (Content-LT), Vmax (Vmax-LT), and clearance (CL-LT). The Content-LT, Vmax-LT, and CL-LT for each individual were obtained by multiplying each individual microsomal phenotype parameter (Content-M, Vmax-M, and CL-M) by the corresponding individual microsomal protein per gram of liver (MPPGL). MPPGL is the microsomal protein per gram of liver and determined as previously reported (Zhang et al., 2015b). The mean MPPGL value was 42.0 (9.9–127.9) mg/g in this cohort of 90 liver samples.(3)
Organ Phenotype.
The organ phenotype is defined as content (Content-L), Vmax (Vmax-L), and clearance (CL-L) in liver. Content-L, Vmax-L, and CL-L were obtained by multiplying each individual liver weight (LW)/body weight (BW) by the individual tissue phenotype parameter, also known as Content-LT, Vmax-LT, and CL-LT, respectively. Detailed values of LW and BW were previously reported (Gao et al., 2017a). BW is the actual body weight for each individual, and LW is the liver weight, which is calculated by multiplying the liver volume (LV) by the liver density, wherein LV (ml) = 12.5 × BW (kg) + 536.4 (Wang et al., 2008), and the liver density is 1.001 g/ml (Yuan et al., 2008). The mean LW/BW value was 21.2 (18.4–30.4) g/kg in this cohort of 90 liver samples.
Organismal Phenotype.
The organismal phenotype is defined as clearance (CL-H) in vivo. According to our previous report, the Bias-Corrected Conventional in vitro–in vivo extrapolation method (Gao et al., 2016a, 2017a) was adopted for extrapolating the CL-H for different genotypes on the basis of QH, fu,p, and RB. CO, the cardiac output (CO), determined based on age and gender, was 5.45 (4.92–6.65) l/min. QH, the hepatic blood flow, determined as 24.5% of the CO (Barter et al., 2013), was 1334.9 (1205.4–1629.3) ml/min. The plasma unbound fraction, fu,p, determined as 0.055 for coumarin and RB, ratio of the drug concentration in blood to plasma, and as 1 for coumarin, was obtained from our previous study (Gao et al., 2016a). Detailed data and the computing method for QH, fu,p, and RB were reported in our previous studies (Zhang et al., 2015b; Gao et al., 2016a). The clearance for CYP2A6 in vivo (CL-H) was calculated using the well-stirred model:(4)
The average fold error (AFE) was introduced to assess the overall accuracy of prediction, which was calculated as follows:(5)
The individual fold error (IFE) was introduced to evaluate the individual accuracy, which was calculated as:(6)
N refers to the number of separate reports in the literature concerning coumarin intravenous clearance.
The definitions and corresponding units for parameters at each level were listed in Table 1.
Definitions of phenotypes for CYP2A6 at different levels
Molecular level, or phenotype at the protein level, Vmax-P (pmol/min per picomole), CL-P (μl/min per picomole). Subcellular level, or phenotype at microsomal level, Vmax-M (pmol/min per milligram), CL-M (μl/min per milligram), Content-M (pmol/mg). Tissue level, or phenotype at liver tissue level, Vmax-LT (pmol/min per gram), CL-LT (μl/min per gram), Content-LT (pmol/g). Organ level, or phenotype at liver level, Vmax-L (nmol/min per kilogram), CL-L (μl/min per kilogram), Content-L (nmol/kg). Organism level, or phenotype (clearance) at in vivo level, CL-H (ml/min).
Intraindividual Variation of CYP2A6 at Different Levels
To make CYP2A6 phenotypes at different levels comparable, relative Vmax (R-Vmax), relative CL (R-CL), and relative Content (R-Content) were computed. For one specific level, R-Vmax, R-CL, and R-Content were calculated as the individual Vmax, CL, and Content divided by the median for that level. Variation in R-Vmax, R-CL, and R-Content among different levels for each person was expressed as the intraindividual percentage coefficient of variation (ICV).(7)
(8)
For one specific relative phenotype (R-Vmax, R-CL or R-Content), the ICV was calculated as the S.D. divided by the corresponding average value of the relative phenotype among different levels.
To compare the variation trend of Vmax and CL from the molecular to the organismal level, the percentage rank change of individual CYP2A6 phenotypes was developed. Ninety samples were ranked in ascending order according to the value of each metabolic phenotype parameter: Vmax-P, Vmax-M, Vmax-LT, Vmax-L, CL-P, CL-M, CL-LT, CL-L, and CL-H. The rank change for each individual was determined by comparison with phenotype at the molecular level (Vmax-P or CL-P), calculation of the absolute difference value between the rank of Vmax or CL at different levels (Vmax-M, Vmax-LT, Vmax-L, or CL-M, CL-LT, CL-L, and CL-H), and the rank of corresponding Vmax-P or CL-P. The percentage of rank change was computed for the sum of the 90 samples divided by the rank change absolute value of each individual. For CYP2A6 content, content phenotype at the molecular level could not be obtained, and so analysis of its rank change could not be determined.
Intergenotype Variation of CYP2A6 at Different Levels
Vmax-P, Vmax-M, Vmax-LT, Vmax-L, Content-M, Content-LT, Content-L, CL-P, CL-M, CL-LT, CL-L, and CL-H were sorted by CYP2A6 genotype. As compared with wild type, fold changes of genotype were calculated by dividing the absolute value of the median for the phenotype by the absolute value of median for the corresponding wild type. The fold change for all phenotypes with wild type was defined as one.
Statistical Analysis
Statistical analysis was performed using SPSS 22.0 software (SPSS Inc., Chicago, IL). Kolmogorov-Smirnov and Shapiro-Wilk tests were introduced to assess the normality of the data distribution. One-way ANOVA and Bonferroni tests were used for the comparison of relative phenotype potency. Since most data sets of phenotypes, including kinetic parameters and protein content, were not normally distributed, nonparametric methods were generally used. Mann-Whitney U and Kruskal-Wallis H tests were used for pairwise comparison and multiple pairwise comparisons of phenotypes with different genetic polymorphisms. Nonparametric Spearman rank correlation analysis was conducted to calculate the correlation coefficient (r) among different phenotypes. P < 0.05 (two-tailed) was considered statistically significant. Graphs were generated using Graph Pad Prism software version 6.02. Both the genotyping and phenotyping tests were derived from three separate and independent replicates to obtain an accurate value. Experimental results were averaged to obtain a single value for that experimental series.
Results
Genotype of CYP2A6
The three most common genetic variants with an allele frequency of more than 1% in the Chinese population were analyzed for a total of 90 samples. The allele CYP2A6*1B has a gene conversion located in the 3′–untranslated region, with the frequency of the wild-type, heterozygous, and homozygous genotypes being 12.22%, 56.67%, and 17.78%, respectively. CYP2A6*4 is a CYP2A6 gene-deletion variant, with the frequency of the wild type and heterozygotes being 86.67% and 13.33%, respectively, and with no homozygous genotype detected. The CYP2A6*9 (−48T>G) is located in the TATA box, with the frequency of the wild-type, heterozygous, and homozygous genotypes being 51.11%, 32.22%, and 3.33%, respectively.
Phenotype of CYP2A6
Phenotypes at Different Levels.
Substantial interindividual variations in the metabolic phenotypes (Vmax, CL) and protein phenotypes (Content) for CYP2A6 from the molecular to cellular level and then to organismal levels are depicted in Table 2, with the greatest difference demonstrated in CL (454-fold) at microsomal level, 144.8 (1.20∼544.7) μl/min per milligram. The CL phenotype exhibited the most remarkable variations, followed by Vmax, and CYP2A6 content was relatively smaller at most levels. Variations of Vmax showed an increased trend from molecular (63.8-fold) to the tissue level (247-fold) and then to the organ level (134-fold).
CYP2A6 phenotypes at different levels (n = 90)
Relative phenotypes, R-Vmax, R-CL, R-Content, referring to relative Vmax, CL, and Content, respectively, were computed to make CYP2A6 phenotypes at different levels comparable. Each level’s relative Vmax, CL, and Content (R-Vmax, R-CL, and R-Content) were calculated as the individual Vmax, CL, and Content divided by the median of corresponding parameter for that level. Molecular level, or phenotype at the protein level, Vmax-P (pmol/min per picomole), CL-P (μl/min per picomole). Subcellular level, or phenotype at microsomal level, Vmax-M (pmol/min per milligram), CL-M (μl/min per milligram), Content-M (pmol/mg). Tissue level, or phenotype at liver tissue level, Vmax-LT (pmol/min per gram), CL-LT (μl/min per gram), Content-LT (pmol/g). Organ level, or phenotype at liver level, Vmax-L(nmol/min per kilogram), CL-L (μl/min per kilogram), Content-L (nmol/kg). Organism level, or phenotype (clearance) at in vivo level, CL-H (ml/min).
To make CYP2A6 phenotypes at different levels comparable, relative Vmax (R-Vmax), CL (R-CL), and Content (R-Content) were computed. Each level’s R-Vmax, R-CL, and R-Content were calculated as the individual Vmax, CL, and Content divided by the median for that level. To sort each level’s relative phenotypes in a strong-to-weak sequence, we evaluated the relative potencies at each level (R-Vmax, R-CL, and R-Content) in 90 samples. It is noteworthy that the R-Vmax potencies at the organ (liver) level, compared with other levels, were always higher. Analysis for the values of R-Vmax of four different levels revealed that the R-Vmax for the organ level was significantly higher than that of the other three levels (Table 2). For values of R-CL and R-Content, statistically significant differences were not demonstrated among different levels.
Intraindividual Variation among Different Levels.
The relative Vmax (R-Vmax), CL (R-CL), and Content (R-Content) values for different levels are depicted in Fig. 1. The ICVs of R-Vmax, R-CL, and R-Content within different levels were 41.0% (18.8%–125.1%), 28.5% (2.39%–133.5%), and 27.8% (2.68%–88.0%) for 90 cases (Fig. 1A). ICVs with the values >100%, 50%–100%, and <50% are indicative of high, moderate, and low variation, respectively. In these 90 cases, the number (percentage) of high, moderate, and low ICVs for the values of R-Vmax within four levels was 4 (4.44%), 10 (11.1%), and 76 (84.4%); for R-CL, they were 4 (4.44%), 24 (26.7%), and 62 (68.9%); and for R-Content, they were 15 (16.7%) and 75 (83.3%) for moderate and low ICV, respectively, with no high ICV detected. The results indicated surprisingly high intraindividual variation for CYP2A6 phenotypes among different levels. When relative phenotypes were classified by CYP2A6 genotypes, the ICVs of R-Vmax, R-CL, and R-Content within different levels varied (Fig. 1, B–D). One notable example was for the CYP2A6*4 polymorphism, with R-Vmax being 39.1% (18.8%–100.8%) and 50.0% (25.9%–125.0%) and R-CL being 26.7% (2.39%–102.9%) and 41.6% (11.1%–133.5%) for *1/*1 and *1/*4 genotypes, respectively. For R-Content, the genotype-dependent effect was less obvious.
ICVs of relative Vmax, CL, and Content of CYP2A6 at different levels (n = 90). ICVs of relative Vmax, CL, and Content of CYP2A6 for total samples (A) and ICVs of relative Vmax (B), CL (C), and Content (D) sorted by different CYP2A6 genotypes. ICVs were calculated as the S.D. of relative Vmax, CL, and Content divided by the corresponding average value for that case. Each level’s relative Vmax, CL, and Content were calculated as the individual Vmax, CL, and Content divided by the median for that level. Horizontal lines represent medians with the interquartile range. *1A/*1A, n = 11; *1A/*1B, n = 51; *1B/*1B, n = 16. *1/*1, n = 78; *1*4, n = 12.*1/*1, n = 46; *1/*9, n = 29; *9/*9, n = 3.
To obtain more detailed information on how the alteration of CYP2A6 phenotypes among different levels originated, values for MPPGL, LW/BW, and QH (hepatic blood flow) were determined for individuals with CYP2A6 genetic differences (Fig. 2). Metabolic phenotypes at the liver tissue level (Vmax-LT, CL-LT) were obtained by multiplying each individual MPPGL by the corresponding individual Vmax or CL at the microsomal level. LW/BW and QH were introduced to determine phenotypes at the organ and organismal levels, respectively. However, no significant effect of CYP2A6 genotype was found.
Effect of CYP2A6 genotypes on MPPGL (A), LW/BW (B), and QH (C). QH, hepatic blood flow. Horizontal lines represent medians with the interquartile range. *1A/*1A, n = 11; *1A/*1B, n = 51; *1B/*1B, n = 16.*1/*1, n = 78; *1*4, n = 12. *1/*1, n = 46; *1/*9, n = 29; *9/*9, n = 3.
Genotype-Phenotype Association
No significant effect of demographic factors, including gender, age, smoking status, drinking habit, and medical diagnoses (Supplemental Table 1), on phenotypes was detected (Supplemental Figs. 1–5).
CYP2A6*1B.
No significant differences were found for Vmax and CL at different levels among genetically different individuals with respect to CYP2A6*1B polymorphisms, nor was a significant effect on the Km value detected. As for the content of CYP2A6 at different levels, the effect varied, with a minor impact on Content-LT and Content-L, as demonstrated by significantly higher Content-LT and Content-L for individuals carrying CYP2A6*1B/*1B homozygote as compared with *1A/*1B heterozygotes; however, there was no significant difference as compared with wild type. It is worth mentioning that strongly significant correlations were charted between Km and Vmax at different levels among most of the genetic variants and were dependent upon CYP2A6 variant. For instance, for Vmax at microsomal level (Vmax-M), extremely high correlation was detected with Km for all CYP2A6*1B-related variants (Fig. 3); however, in terms of Vmax at the isoform level (Vmax-P), no significant correlation was obtained with Km for CYP2A6*1A/*1A wild type, and moderate correlation was obtained for *1A/*1B and *1B/*1B subjects. The mechanism underlying this observation remained to be determined.
Genotype-phenotype associations varied by CYP2A6*1B polymorphisms. A1∼A4 (Vmax-P, Vmax-M, Vmax-LT, Vmax-L), B1∼B4 (Km, Content-M, Content-LT, Content-L), C1∼C4 (CL-P, CL-M, CL-LT, CL-L) were divided by *1A/*1A, *1A/*1B, and *1B/*1B genotypes (*1A/*1A, n = 11; *1A/*1B, n = 51; *1B/*1B, n = 16). D1∼D4, E1∼E4, G1∼G4, correlation between Km with Vmax-P, Vmax-M, Vmax-LT, Vmax-L in *1A/*1A, *1A/*1B and *1B/*1B genotypes, respectively. Vmax-P (pmol/min per picomole), CL-P (μl/min per picomole), the molecular phenotype or phenotype at the protein level. Vmax-M (pmol/min per milligrams), CL-M (μl/min per milligram), Content-M (pmol/mg), the subcellular phenotype or phenotype at microsomal level. Vmax-LT (nmol/min per gram), CL-LT (ml/min per gram), Content-LT (pmol/g), the tissue phenotype or phenotype at liver tissue level. Vmax-L (μmol/min per kilogram), CL-L (ml/min per kilogram), Content-L (pmol/kg), the organ phenotype or phenotype at liver level. Km (μM), the Michaelis-Menten constant.
CYP2A6*4.
CYP2A6*4 was significantly associated with decreased Vmax, and the effect on phenotype at the cellular level was more markedly significant than that at the molecular level. The Vmax-P values of CYP2A6*1/*4 individuals (15.24 pmol/min per picomole) were dramatically lower than *1/*1 (26.45 pmol/min per picomole), and the Vmax-M, Vmax-LT, and Vmax-L of CYP2A6*1/*4 individuals were significantly less than a third of *1/*1 genotypes (Fig. 4). Subjects carrying CYP2A6*1/*4 genotypes showed lower content-M (10.62 pmol/mg) and content-LT (370 pmol/g) when compared with content-M (18.92 pmol/mg) and content-LT (629 pmol/g) for subjects with *1/*1. Observably lower values of content-M and content-LT for *1/*4 individuals than those in *1/*1 individuals might explain significantly greater lowered Vmax-M and Vmax-LT in *1/*4 individuals than *1/*1 subjects. Consistent with this, we found that the Vmax-L values of *1/*4 individuals (110.2 nmol/min per kilogram) were significantly lower than *1/*1 (366.2 nmol/min per kilogram) individuals.
Genotype-phenotype associations varied by CYP2A6*4 polymorphisms. A1∼A4 (Vmax-P, Vmax-M, Vmax-LT, Vmax-L), B1∼B4 (Km, Content-M, Content-LT, Content-L), C1∼C4 (CL-P, CL-M, CL-LT, CL-L) were divided by *1/*1 and *1/*4 genotypes (*1/*1, n = 78; *1*4, n = 12). D1∼D4, E1∼E4, correlation between Km with Vmax-P, Vmax-M, Vmax-LT, Vmax-L in *1/*1 and *1/*4 genotypes, respectively. Vmax-P (pmol/min per picomole), CL-P (μl/min per picomole), the molecular phenotype or phenotype at the protein level. Vmax-M (pmol/min per milligram), CL-M (μl/min per milligram), Content-M (pmol/mg), the subcellular phenotype or phenotype at microsomal level. Vmax-LT (nmol/min per gram), CL-LT (ml/min per gram), Content-LT (pmol/g), the tissue phenotype or phenotype at liver tissue level. Vmax-L (μmol/min per kilogram), CL-L (ml/min per kilogram), Content-L (pmol/kg), the organ phenotype or phenotype at liver level. Km (μM), the Michaelis-Menten constant.
It should be stated that individuals with the CYP2A6*1/*4 genotype showed lower Km than those with the *1/*1 genotype. No statistical difference was found for CL-P, which was obtained by dividing each individual Vmax-P by the individual Km, between the*1/*1 (10.76 µl/min per picomole) and *1/*4 (10.83 µl/min per picomole) genotypes, and this might be due to the simultaneously lowered Vmax-P and Km in individuals and the effect on Vmax-P partially counteracting the negative effect on Km, as evidenced by significant correlation between Km and Vmax-P (r = 0.457, P = 2.570 × 10−5) for *1/*1 individuals; however, no such correlation was seen for *1/*4 individuals. The CL-M (75.05 µl/min per milligram), CL-LT (3.70 ml/min per gram), and CL-L (76.97 ml/min per kilogram) of CYP2A6 for subjects with *1/*4 were dramatically lower than those for *1/*1 individuals (150.11 µl/min per milligram, 6.12 ml/min per gram, and 130.12 ml/min per kilogram, respectively). This could be due to the simultaneously lowered Vmax-P, Vmax-M, Vmax-LT, Vmax-L, and Km in individuals, as evidenced by a significant correlation between Km and Vmax at different levels, especially for the extremely strong correlation between Km and Vmax-M (r = 0.908, P = 2.092 × 10−30), Km and Vmax-LT (r = 0.751, P = 2.285 × 10−15), and Km and Vmax-L (r = 0.751, P = 1.024 × 10−13) for *1/*1 individuals; however, no significant correlation was seen for *1/*4 individuals.
CYP2A6*9.
Consistent with CYP2A6*1B, the effect of CYP2A6*9 on phenotypes was not as significant as with CYP2A6*4 (Fig. 5). Only subjects with CYP2A6*1/*1 wild type and *1/*9 heterozygotes were taken into analysis, without considering *9/*9 homozygotes since only three individuals carried the *9/*9 genotype. No significant difference was found in Vmax-P for *1/*1 carriers as compared with *1/*9 carriers, and there was a dramatically lower Km for individuals with *1/*9 (2.07 µM) than for*1/*1 (3.03 µM) genotypes. The CL-P values, acquired by dividing individual Vmax-P by each Km as described above, of *1/*9 (12.72 µl/min per picomole) individuals seemed to be higher than*1/*1 (8.71 µl/min per picomole) subjects; however, this difference was not statistically significant. For both Vmax-M and Vmax-LT, *1/*9 was associated with decreased activity as compared with wild type, which might be related to seemingly lower content-M and content-LT for *1/*9 individuals in spite of lack of statistical difference. Moreover, the CL-M for 2A6*1/*9 was significantly lower than that of *1/*1 individuals, which might be due to simultaneously lowered Vmax-M and Content-M and a strongly significant association between Km and Vmax-M. For Vmax-L, Content-LT, Content-L, CL-LT, and CL-L, no significant alteration was related with *9 polymorphisms, even with a lower tendency of Vmax-L, Content-LT, and Content-L for *1/*9 individuals. We note that Km was strongly positively associated with Vmax within different levels for both CYP2A6*1/*1 and CYP2A6*1/*9 genotypes except for Vmax-P in CYP2A6*9 wild-type carriers.
Genotype-phenotype associations varied by CYP2A6*9 polymorphisms. A1∼A4 (Vmax-P, Vmax-M, Vmax-LT, Vmax-L), B1∼B4 (Km,Content-M, Content-LT, Content-L), C1∼C4 (CL-P, CL-M, CL-LT, and CL-L) were divided by *1/*1, *1/*9, and *9/*9 genotypes (*1/*1, n = 46; *1/*9, n = 29; *9/*9, n = 3). D1∼D4, E1∼E4, correlation between Km with Vmax-P, Vmax-M, Vmax-LT, Vmax-L in *1/*1 and *1/*9 genotypes, respectively. Vmax-P (pmol/min per picomole), CL-P (μl/min per picomole), the molecular phenotype or phenotype at the protein level. Vmax-M (pmol/min per milligram), CL-M (μl/min per milligram), Content-M (pmol/mg), the subcellular phenotype or phenotype at microsomal level. Vmax-LT (nmol/min per gram), CL-LT (ml/min per gram), Content-LT (pmol/g), the tissue phenotype or phenotype at the liver tissue level. Vmax-L (μmol/min per kilogram), CL-L (ml/min per kilogram), Content-L (pmol/kg), the organ phenotype or phenotype at liver level. Km (μM), the Michaelis-Menten constant.
Intergenotype Variation
Intergenotype individual variations of CYP2A6 phenotypes at different levels were charted, and significant differences in variability among genetically different individuals are presented as seen in Fig. 6. Generally speaking, the CYP2A6*4 polymorphism led to the most remarkable alteration, with decreased Vmax, Content, and CL at the majority of levels, whereas CYP2A6*1B and *9 led to less obvious alteration either with increased or decreased phenotypes at some levels. Compared with wild type, various phenotypes including Vmax-P, Vmax-M, Vmax-LT, Vmax-L, Content-M, Content-LT, CL-M, CL-LT, CL-L, and CL-H of *1/*4 showed an approximate decrease of 51%, 85%, 80%, 88%, 51%, 51%, 71%, 59%, 62%, and 57%, respectively, whereas the Content-L and CL-P were basically unchanged. For CYP2A6*1B, Content-LT, Content-L, CL-LT, CL-L, and CL-H of mutant homozygotes were significantly higher than that of heterozygotes, whereas heterozygotes and mutant homozygotes were basically unchanged when compared with wild type. Most notably, CYP2A6*1B led to distinctly increased Content-L with heterozygotes and mutant homozygotes showing 2.80-fold and 6.87-fold increases as compared with wild type, respectively. However, CL-L did not confirm this trend of change. With regard to CYP2A6*9, lowered Vmax-P, Vmax-M, Vmax-LT, Content-L, and CL-L of mutant heterozygotes were detected when compared with wild type, whereas Content-M, CL-P, and CL-M showed a reverse trend.
The diversity in CYP2A6 phenotypes at different levels. A1∼A3 (fold changes of Vmax-P, Vmax-M, Vmax-LT, Vmax-L), A4 (fold changes of Km), B1∼B3 (fold changes of Content-M, Content-LT, Content-L), C1∼C3 (fold changes of CL-P, CL-M, CL-LT, CL-L, CL-H)were sorted by different CYP2A6 genotypes. As compared with wild type, fold changes of other genotypes were calculated by dividing the absolute value of the median for the phenotype by the absolute value of median for the corresponding wild type. The fold change for all phenotypes with wild type were defined as one. CYP2A6*1B (*1A/*1A, n = 11; *1A/*1B, n = 51; *1B/*1B, n = 16); CYP2A6*4 (*1/*1, n = 78; *1*4, n = 12); CYP2A6*9 (*1/*1, n = 46; *1/*9, n = 29; *9/*9, n = 3). *P < 0.05, **P < 0.01, ***P < 0.001 vs. wild type; #P < 0.05, ##P < 0.01 vs. mutant heterozygotes.
Moreover, with regard to Km, in line with other phenotypes, the most remarkable alteration was with *1/*4 genotypes, which showed a 49% reduction, whereas *9 polymorphisms led to distinctly lower Km, with heterozygotes and mutant homozygotes showing 32% and 47% decreases as compared with wild type, respectively. For CYP2A6*1B, Km of *1B carriers was basically unchanged for heterozygotes and homozygotes.
Discussion
In this present study of 90 human liver samples, the systematic delineation of a CYP2A6 genotype-phenotype model including protein content and metabolic activity of CYP2A6 at various levels such as molecular, cellular (microsomal, tissue, organ), and organismal levels was generated. We determined the median and range for the protein content, Vmax, and CL as well as relative phenotypes (R-Content, R-Vmax, and R-CL) at each level with every CYP2A6 genotype to elucidate how different CYP2A6 genotypes yield phenotypic traits (protein content and enzyme activity) at different levels. The high relative ICVs of phenotypes indicate large intraindividual variation within different levels, sometimes in a genotype-dependent manner. Intergenotype analysis indicated that CYP2A6*4 had the most obvious effect on phenotypes, and protein content might be related to the CYP2A6*4 allele. The correlation of phenotypes from the molecular level to the cellular level and then to the organismal level demonstrates how genotypes may affect phenotypes and could be considered for better understanding of phenotypic variation and assessment of phenotypic outcome with specific genotypes.
We did find that the alteration of phenotypes from the molecular level to the cellular level and then to organismal level varied greatly, which agrees with our previous studies (Zhang et al., 2015b; Gao et al., 2017a). However, our previous study on the CYP2D6 genotype-phenotype association (Gao et al., 2017a) indicated the effect of the CYP2D6 genotype on the metabolic phenotype at each level but did not consider the alteration of phenotypes among different levels. In this study, the ICVs of R-Vmax, R-CL, and R-Content within different levels indicated noticeable intraindividual variation and varied with the CYP2A6 genotype. Take the CYP2A6*4 allele, for example; the ICVs of R-CL for CYP2A6*1/*4, 41.6% (11.1%–133.5%), was nearly 1.56-fold of that for the *1/*1 genotype, 26.7% (2.39%–102.9%), whereas for R-Vmax, the difference was 1.28-fold. The reason for these differences might be the superimposed effect of genotype on phenotype since the effect on the molecular level might influence the cellular level and then the organismal level. In line with this speculation, the effect on CL might be the combined effect on Vmax and Km since CL was obtained by dividing Vmax by Km. However, the precise mechanism underlying this needs further exploration.
The most impactful genotype-on-phenotype effect was with CYP2A6*4, which is a gene-deletion genotype, in contrast to the *1B genotype (located in 3′ untranslated region) and the *9 mutation (located in the “TATA box”), being the least impactful phenotype as measured herein. Results are in agreement with a series of previous studies (Ujjin et al., 2002; Yoshida et al., 2003; Wang et al., 2006, 2011; Tanner et al., 2017), though those studies were conducted either in vitro by site-directed mutagenesis or in vivo by phenotyping using substrates. The location and type of mutation is clearly important when evaluating the phenotypic impact of different genotypes.
It is noteworthy that protein content might play an important role for extrapolation from the genotype of CYP2A6 to the different levels of phenotypes, especially for CYP2A6*4. CYP2A6*4, a coding-deletion genotype, was associated with significant decreased protein and acitvity in our study and also with the same effect in previous reports (Peamkrasatam et al., 2006; Wang et al., 2011). However, one such trend was detected in one study demonstrating protein content correlated well with activity (Tanner et al., 2017). The reason for the inconsistencies might be that liver samples conducted in that study were from two liver banks with different ethnicities, and only one CYP2A6*1/*4 sample was identified, and so the effect of the CYP2A6*4 genotype on phenotype could not be well evaluated. Our study focused on known but less common polymorphisms that were present at higher levels in the Chinese population. As depicted in Fig. 4, the markedly lower Content-M (10.62 pmol/mg) and Content-LT (3.70 ml/min per gram) in *1/*4 genotypes might help explain the much greater decrease in Vmax-M, Vmax-LT, and Vmax-L compared with Vmax-P correlated with CYP2A6*4. This was consistent with the view that protein content and activity might play a role in determining metabolic phenotypes and susceptibility to disease (Tian et al., 2016; Gao et al., 2018a,b,c).
Most notably, one interesting finding was the highly significant positive correlation between Km and Vmax at each level in the*1/*1 wild type. Given that, we speculate that the significantly lowered CL-M, CL-LT, and CL-L in *4 carriers might be due to the more prominent decreased Vmax when compared with Km at each corresponding level, whereas for the phenotype at the molecular level, the almost uniform extent of change for Vmax-P and Km might be the reason for the insignificant alteration of CL-P related to the *4 variation. However, the explicit mechanism underlying this remains unknown and needs further investigation.
The interindividual variation of CYP2A6 phenotypes in relation to different genotypes at various levels may play a role in interethnic or interindividual variability for CYP2A6 phenotypes and the risk for CYP2A6-related disease. CYP2A6*4, in particular, correlated with lower 2A6 phenotypes at a number of levels. For example, the median values of Vmax-M, Content-M, and CL-M were decreased to 15%, 49%, and 29%, respectively, in *4 carriers when compared with *1 wild-type carriers. This likely contributes to the known slower cotinine metabolism in Chinese-Americans than Caucasians (Benowitz et al., 2002) since the frequency of CYP2A6*4 in Chinese was 6.6%∼15.1%, much higher than the Caucasian level of 0.5%∼4.9% (Xu et al., 2002).
The interethnic as well as interindividual variation in metabolic phenotypes might have several important clinical implications; for instance, CYP2A6-related metabolic activity (e.g., nicotine metabolism), behaviors (Malaiyandi et al., 2006; Gao et al., 2016a) (i.e., smoking), and disease (Fujieda et al., 2004) (i.e., lung cancer). Genetic studies on CYP2A6 in Japanese indicated a decreased rate of lung cancer in inactive CYP2A6 variants carriers in comparison with wild type (Miyamoto et al., 1999). In addition to the cigarette smoking–related lung cancer, CYP2A6 is also involved in the metabolism of many other procarcinogens. However, a clearer understanding of the role for genotypes of CYP2A6 in relation to disease is needed before conclusions can be reached.
Since liver samples of perfectly healthy persons with comprehensive background information are difficult to obtain, the liver samples used here were from patients undergoing hepatic surgery. However, the donors were with clear medical diagnoses mostly of cavernous hemangioma of liver and mitigated by techniques and confirmed by analysis (Supplemental Fig. 5) to minimize impact of tumors and other demographic factors (Supplemental Figs. 1–4). Moreover, the total of 90 samples studied here with comprehensive background information is a large number of human samplings and is certainly adequate. Even so, high variability existed along with relatively low numbers per grouping and might still affect the ability to differentiate smaller differences. Many more samples are needed to rule out idiosyncratic variation as much as possible. Since it was not easy to simultaneously obtain in vivo and in vitro results for the same cohort of samples, no actual in vivo measurements were made in our study. We used in vitro measurements at different levels to predict hepatic clearance by in vitro–in vivo extrapolation to the in vivo level. After correction, our in vivo extrapolation by calculation yielded high accuracy, as indicated by the predicted CL of CYP2A6 matching 82.8% of the cases within 2-fold error when compared with actual observed CLin vivo determined in other studies (Ritschel et al., 1977)(Gao et al., 2016a). Even so, actual in vivo measurements are additional efforts needed to draw more clear and convincing conclusions.
In conclusion, we have developed a comprehensive CYP2A6 genotype-to-phenotype model. In this model, a role for the CYP2A6 genotypes in influencing phenotypic variation at various levels was evaluated, although clearly it may not be sufficient to explain variation in all phenotypic outcomes. It does provide the basal data for CYP2A6 content and metabolic activities at different levels with each CYP2A6 genotype as well as fundamental hierarchical relationships between CYP2A6 genotypes and phenotypes. We identified a key role for decreased protein content for the CYP2A6*4 genotype that is related to decreased metabolic phenotypes. These findings might help open new approaches for considering how genetic variation is translated to the phenotype of individuals, which might be a further step on the way to making useful predictions of outcomes for a particular individual as a part of personalized medicine and disease risk.
Acknowledgments
The authors would like to thank professor Song-He Shi and Xue-Zhong Shi for their valuable suggestions during the statistical analysis of data for this manuscript.
Authorship Contributions
Participated in research design: Qiao.
Conducted experiments: Fang, Wang.
Contributed new reagents or analytic tools: Guo, Zhang, Wen.
Performed data analysis: Fang, Wang.
Wrote or contributed to the writing of the manuscript: Fang, Xing, Gao, Qiao.
Footnotes
- Received October 9, 2019.
- Accepted December 23, 2019.
This work was supported by the National Natural Science Foundation of China [Grants 81473279, 81673507, 8187 2931, and U1804169].
↵
This article has supplemental material available at jpet.aspetjournals.org.
Abbreviations
- BW
- body weight
- CL
- clearance
- HLM
- human liver microsome
- ICV
- intraindividual percentage coefficient of variation
- LC
- liquid chromatography
- LW
- liver weight
- MPPGL
- microsomal protein per gram of liver
- MS
- mass spectrometry
- Copyright © 2020 by The American Society for Pharmacology and Experimental Therapeutics