Molecular analysis of the CYP2F1 gene: identification of a frequent non-functional allelic variant

Complementarity of two proteomic data analysis tools in the identification of drug-metabolising enzymes and transporters in human liver

Several software packages are available for the analysis of proteomic LC-MS/MS data, including commercial (e.g. Mascot/Progenesis LC-MS) and open access software (e.g. MaxQuant). In this study, Progenesis and MaxQuant were used to analyse the same data set from human liver microsomes (n = 23). Comparison focussed on the total number of peptides and proteins identified by the two packages. For the peptides exclusively identified by each software package, distribution of peptide length, hydrophobicity, molecular weight, isoelectric point and score were compared. Using standard cut-off peptide scores, we found an average of only 65% overlap in detected peptides, with surprisingly little consistency in the characteristics of peptides exclusively detected by each package. Generally, MaxQuant detected more peptides than Progenesis, and the additional peptides were longer and had relatively lower scores. Progenesis-specific peptides tended to be more hydrophilic and basic relative to peptides detected only by MaxQuant. At the protein level, we focussed on drug-metabolising enzymes (DMEs) and transporters, by comparing the number of unique peptides detected by the two packages for these specific proteins of interest, and their abundance. The abundance of DMEs and SLC transporters showed good correlation between the two software tools, but ABC showed less consistency. In conclusion, in order to maximise the use of MS datasets, we recommend processing with more than one software package. Together, Progenesis and MaxQuant provided excellent coverage, with a core of common peptides identified in a very robust way.

Genotypic and phenotypic landscapes of 51 pharmacogenes derived from whole-genome sequencing in a Thai population

Differences in drug responses in individuals are partly due to genetic variations in pharmacogenes, which differ among populations. Here, genome sequencing of 171 unrelated Thai individuals from all regions of Thailand was used to call star alleles of 51 pharmacogenes by Stargazer, determine allele and genotype frequencies, predict phenotype and compare high-impact variant frequencies between Thai and other populations. Three control genes, EGFR, VDR, and RYR1, were used, giving consistent results. Every individual had at least three genes with variant or altered phenotype. Forty of the 51 pharmacogenes had at least one individual with variant or altered phenotype. Moreover, thirteen genes had at least 25% of individuals with variant or altered phenotype including SLCO1B3 (97.08%), CYP3A5 (88.3%), CYP2C19 (60.82%), CYP2A6 (60.2%), SULT1A1 (56.14%), G6PD (54.39%), CYP4B1 (50.00%), CYP2D6 (48.65%), CYP2F1 (46.41%), NAT2 (40.35%), SLCO2B1 (28.95%), UGT1A1 (28.07%), and SLCO1B1 (26.79%). Allele frequencies of high impact variants from our samples were most similar to East Asian. Remarkably, we identified twenty predicted high impact variants which have not previously been reported. Our results provide information that contributes to the implementation of pharmacogenetic testing in Thailand and other Southeast Asian countries, bringing a step closer to personalized medicine.

Considering smoking status, coexpression network analysis of non-small cell lung cancer at different cancer stages, exhibits important genes and pathways

Non-small cell lung cancer (NSCLC) is the most common subtype of lung cancer among smokers, nonsmokers, women, and young individuals. Tobacco smoking and different stages of the NSCLC have important roles in cancer evolution and require different treatments. Existence of poorly effective therapeutic options for the NSCLC brings special attention to targeted therapies by considering genetic alterations. In this study, we used RNA-Seq data to compare expression levels of RefSeq genes and to find some genes with similar expression levels. We utilized the "Weighted Gene Co-expression Network Analysis" method for three different datasets to create coexpressed genetic modules having relations with the smoking status and different stages of the NSCLC. Our results indicate seven important genetic modules having important associations with the smoking status and cancer stages. Based on investigated genetic modules and their biological explanation, we then identified 13 newly candidate genes and 7 novel transcription factors in association with the NSCLC, the smoking status, and cancer stages. We then examined those results using other datasets and explained our results biologically to illustrate some important genes in relation with the smoking status and metastatic stage of the NSCLC that can bring some crucial information about cancer evolution. Our genetic findings also can be used as some therapeutic targets for different clinical conditions of the NSCLC.

Xenobiotic metabolism in differentiated human bronchial epithelial cells

Differentiated human bronchial epithelial cells in air liquid interface cultures (ALI-PBEC) represent a promising alternative for inhalation studies with rodents as these 3D airway epithelial tissue cultures recapitulate the human airway in multiple aspects, including morphology, cell type composition, gene expression and xenobiotic metabolism. We performed a detailed longitudinal gene expression analysis during the differentiation of submerged primary human bronchial epithelial cells into ALI-PBEC to assess the reproducibility and inter-individual variability of changes in transcriptional activity during this process. We generated ALI-PBEC cultures from four donors and focussed our analysis on the expression levels of 362 genes involved in biotransformation, which are of primary importance for toxicological studies. Expression of various of these genes (e.g., GSTA1, ADH1C, ALDH1A1, CYP2B6, CYP2F1, CYP4B1, CYP4X1 and CYP4Z1) was elevated following the mucociliary differentiation of airway epithelial cells into a pseudo-stratified epithelial layer. Although a substantial number of genes were differentially expressed between donors, the differences in fold changes were generally small. Metabolic activity measurements applying a variety of different cytochrome p450 substrates indicated that epithelial cultures at the early stages of differentiation are incapable of biotransformation. In contrast, mature ALI-PBEC cultures were proficient in the metabolic conversion of a variety of substrates albeit with considerable variation between donors. In summary, our data indicate a distinct increase in biotransformation capacity during differentiation of PBECs at the air-liquid interface and that the generation of biotransformation competent ALI-PBEC cultures is a reproducible process with little variability between cultures derived from four different donors.

Hypothesis-based weight-of-evidence evaluation and risk assessment for naphthalene carcinogenesis

Inhalation of naphthalene causes olfactory epithelial nasal tumors in rats (but not in mice) and benign lung adenomas in mice (but not in rats). The limited available human data have not identified an association between naphthalene exposure and increased respiratory cancer risk. Assessing naphthalene's carcinogenicity in humans, therefore, depends entirely on experimental evidence from rodents. We evaluated the respiratory carcinogenicity of naphthalene in rodents, and its potential relevance to humans, using our Hypothesis-Based Weight-of-Evidence (HBWoE) approach. We systematically and comparatively reviewed data relevant to key elements in the hypothesized modes of action (MoA) to determine which is best supported by the available data, allowing all of the data from each realm of investigation to inform interpretation of one another. Our analysis supports a mechanism that involves initial metabolism of naphthalene to the epoxide, followed by GSH depletion, cytotoxicity, chronic inflammation, regenerative hyperplasia, and tumor formation, with possible weak genotoxicity from downstream metabolites occurring only at high cytotoxic doses, strongly supporting a non-mutagenic threshold MoA in the rat nose. We also conducted a dose-response analysis, based on the likely MoA, which suggests that the rat nasal MoA is not relevant in human respiratory tissues at typical environmental exposures. Our analysis illustrates how a thorough WoE evaluation can be used to support a MoA, even when a mechanism of action cannot be fully elucidated. A non-mutagenic threshold MoA for naphthalene-induced rat nasal tumors should be considered as a basis to determine human relevance and to guide regulatory and risk-management decisions.

Metabolic polymorphisms and biomarkers of effect in the biomonitoring of occupational exposure to low-levels of benzene: state of the art

Current levels of occupational exposure to benzene, a genotoxic human carcinogen, in Western countries are reduced by two-three orders of magnitude (from ppm to ppb) as compared to the past. However, as benzene toxicity is strongly dependent on biotransformation and recent evidence underlines a higher efficiency of bio-activation pathways at lower levels of exposure, toxic effects at low doses could be higher than expected, particularly in susceptible individuals. Currently, biological monitoring can allow accurate exposure assessment, relying on sensitive and specific enough biomarkers of internal dose. The availability of similarly reliable biomarkers of early effect or susceptibility could greatly improve the risk assessment process to such an extent that risk could even be assessed at the individual level. As to susceptibility biomarkers, functional genetic polymorphisms of relevant biotransformation enzymes may modulate the risk of adverse effects (NQO1) and the levels of biomarkers of internal dose, in particular S-phenylmercapturic acid (GSTM1, GSTT1, GSTA1). Among biomarkers of early effect, genotoxicity indicators, although sensitive in some cases, are too aspecific for routine use in occupational health surveillance programmes. Currently only the periodical blood cell count seems suitable enough to be applied in the longitudinal monitoring of effects from benzene exposure. Novel biomarkers of early effect are expected from higher collaboration among toxicologists and clinicians, also using advanced "omics" techniques.

Polymorphisms of cytochrome p450 genes in three ethnic groups from Russia

OBJECTIVE

To determine the prevalence of the most common allelic variants of CYP1A1, CYP1A2, CYP1B1, CYP2C9, CYP2E1, CYP2F1, CYP2J2 and CYP2S1 in a representative sample of the three ethnic groups (Russians, Tatars and Bashkirs) from Republic of Bashkortostan (Russia), and compare the results with existing data published for other populations.

MATERIAL AND METHODS

CYPs genotypes were determined in 742 DNA samples of healthy unrelated individuals representative of three ethnic groups. The CYPs gene polymorphisms were examined using the PCR-RLFP method.

RESULTS

Analysis of the CYP1A1 (rs1048943, rs4646903), CYP1A2 (rs762551), CYP2E1 (rs2031920) allele, genotype and haplotype frequencies revealed significant differences among healthy residents of the Republic of Bashkortostan of different ethnicities. Distribution of allele and genotype frequencies of CYP1A2 (rs35694136), CYP1B1 (rs1056836), CYP2C9 (rs1799853, rs1057910), CYP2F1 (rs11399890), CYP2J2 (rs890293), CYP2S1 (rs34971233, rs338583) genes were similar in Russians, Tatars, and Bashkirs. Analysis of the CYPs genes allele frequency distribution patterns among the ethnic groups from the Republic of Bashkortostan in comparison with the different populations worldwide was conducted.

CONCLUSION

The peculiarities of the allele frequency distribution of CYPs genes in the ethnic groups of the Republic of Bashkortostan should be taken into consideration in association and pharmacogenetic studies. The results of the present investigation will be of great help in elucidating the genetic background of drug response, susceptibility to cancer and complex diseases, as well as in determining the toxic potentials of environmental pollutants in our region.

Predicting the effects of frameshifting indels

Each human has approximately 50 to 280 frameshifting indels, yet their implications are unknown. We created SIFT Indel, a prediction method for frameshifting indels that has 84% accuracy. The percentage of human frameshifting indels predicted to be gene-damaging is negatively correlated with allele frequency. We also show that although the first frameshifting indel in a gene causes loss of function, there is a tendency for the second frameshifting indel to compensate and restore protein function. SIFT Indel is available at http://sift-dna.org/www/SIFT_indels2.html

Xenobiotic metabolism and disposition in human lung: transcript profiling in non-tumoral and tumoral tissues

The lung is directly exposed to a wide variety of inhaled toxicants and carcinogens. In order to improve our knowledge of the cellular processing of these compounds in the respiratory tract, we investigated the mRNA expression level of 380 genes encoding xenobiotic-metabolizing enzymes (XME), transporters, nuclear receptors and transcription factors, in pulmonary parenchyma (PP), bronchial mucosa (BM) and tumoral lung tissues from 12 patients with non-small cell lung cancer (NSCLC). Using a high throughput quantitative real-time RT-PCR method, we found that ADH1B, CYP4B1, CES1 and GSTP1 are the major XME genes expressed both in BM and PP. Our results also documented the predominant role played by the xenosensor AhR in human lung. The gene expression profiles were different for BM and PP, with a tendency toward increased mRNA levels of phase I and phase II XME genes in BM, suggesting major differences in the initial stages of xenobiotic metabolism. Some of the significantly overexpressed genes in BM (i.e. CYP2F1, CYP2A13, CYP2W1, NQO1…) encode proteins involved in the bioactivation of procarcinogens, pointing out distinct susceptibility to xenobiotics and their toxic effects between these two tissue types. Additionally, interindividual differences in transcript levels observed for some genes may be of genetic origin and may contribute to the variability in response to environmental exposure and, consequently, in the risk of developing lung diseases. A global decrease in gene expression was observed in tumoral specimens. Some of the proteins are involved in the metabolism or transport of anti-cancer drugs and their influence in the response of tumors to chemotherapy should be considered. In conclusion, the present study provides an overview of the cellular response to toxicants and drugs in healthy and cancerous human lung tissues, and thus improves our understanding of the mechanisms of chemical carcinogenesis as well as cellular resistance to chemotherapy.

Genetic polymorphism of CYP2U1, a cytochrome P450 involved in fatty acids hydroxylation

The human cytochrome P450 2U1 (CYP2U1) has been described as a novel extrahepatic P450. CYP2U1 is a highly conserved gene mainly expressed in brain and thymus, but also at lower levels in kidney, lung or heart. This selective tissue distribution suggests important endogenous functions, in particular in the conversion of arachidonic acid into two bioactive compounds, the 19- and 20-HETE. To investigate the extent of CYP2U1 genetic polymorphism in 70 French individuals, a screening for sequence variations in the 5'-flanking and protein encoding regions was performed using PCR-SSCP and sequencing strategies. Four polymorphisms were identified and correspond to -204C>A and -241T>C in the 5'-flanking region, -37G>A in the 5'-untranslated region, and IVS2-17T>C in the intron 2. The most frequent mutations, -241T>C (59.7%) and IVS2-17T>C (66.0%), did not seem to alter CYP2U1 lung expression. These results suggest that CYP2U1 exhibits few genetic variations and support a probable role in endogenous processes.

The pneumotoxin 3-methylindole is a substrate and a mechanism-based inactivator of CYP2A13, a human cytochrome P450 enzyme preferentially expressed in the respiratory tract

3-Methylindole (3MI), a respiratory tract toxicant, can be metabolized by a number of cytochromes P450 (P450), primarily through either dehydrogenation or epoxidation of the indole. In the present study, we assessed the bioactivation of 3MI by recombinant CYP2A13, a human P450 predominantly expressed in the respiratory tract. Four metabolites were detected, and the two principal ones were identified as indole-3-carbinol (I-3-C) and 3-methyloxindole (MOI). Bioactivation of 3MI by CYP2A13 was verified by the observation of three glutathione (GSH) adducts designated as GS-A1 (glutathione adduct 1), GS-A2 (glutathione adduct 2), and GS-A3 (glutathione adduct 3) in a NADPH- and GSH-fortified reaction system. GS-A1 and GS-A2 gave the same molecular ion at m/z 437, an increase of 305 Da over 3MI. Their structures are assigned to be 3-glutathionyl-S-methylindole and 3-methyl-2-glutathionyl-S-indole, respectively, on the basis of the mass fragmentation data obtained by high-resolution mass spectrometry. Kinetic parameters were determined for the formation of I-3-C (V(max) = 1.5 nmol/min/nmol of P450; K(m) = 14 muM), MOI (V(max) = 1.9 nmol/min/nmol of P450; K(m) = 15 muM) and 3-glutathionyl-S-methylindole (V(max) = 0.7 nmol/min/nmol of P450; K(m) = 13 muM). The structure of GS-A3, a minor adduct with a protonated molecular ion at m/z 453, is proposed to be 3-glutathionyl-S-3-methyloxindole. We also discovered that 3MI is a mechanism-based inactivator of CYP2A13, given that it produced a time-, cofactor-, and 3MI concentration-dependent loss of activity toward 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone, with a relatively low K(I) value of approximately 10 muM and a k(inact) of 0.046 min(-1). Thus, CYP2A13 metabolizes 3MI through multiple bioactivation pathways, and the process can lead to a suicide inactivation of CYP2A13.

CYP2F1 genetic polymorphism: identification of interethnic variations

Since human cytochrome P450 2F1 (CYP2F1) is predominantly expressed in lung tissue and is involved in the metabolism of various pneumotoxicants with potential carcinogenic effects, variations in the nucleotidic sequence of its gene may contribute to interindividual and interethnic differences in the susceptibility to lung tumorigenesis. The aim of the current study was to compare the frequency of a previously reported frameshift mutation, namely c.14_15insC, responsible for the synthesis of a severely truncated protein, between several populations of different ethnic origins. The frequencies of this polymorphism were 26.1, 51.6, 42.7 and 22.9% in French, Gabonese, Senegalese, and Tunisian population samples, respectively, thereby representing a substantial inter ethnic variation in the CYP2F1 gene. These findings provide data for further studies that investigate the potential association of CYP2F1 haplotypes with an incidence of lung cancer genesis in respect of ethnicity.