Genome-Scale Characterization of Toxicity-Induced Metabolic Alterations in Primary Hepatocytes

Context-specific GEnome-scale metabolic Network REconstructions (GENREs) provide a means to understand cellular metabolism at a deeper level of physiological detail. Here, we use transcriptomics data from chemically-exposed rat hepatocytes to constrain a GENRE of rat hepatocyte metabolism and predict biomarkers of liver toxicity using the Transcriptionally Inferred Metabolic Biomarker Response algorithm. We profiled alterations in cellular hepatocyte metabolism following in vitro exposure to four toxicants (acetaminophen, carbon tetrachloride, 2,3,7,8-tetrachlorodibenzodioxin, and trichloroethylene) for six hour. TIMBR predictions were compared with paired fresh and spent media metabolomics data from the same exposure conditions. Agreement between computational model predictions and experimental data led to the identification of specific metabolites and thus metabolic pathways associated with toxicant exposure. Here, we identified changes in the TCA metabolites citrate and alpha-ketoglutarate along with changes in carbohydrate metabolism and interruptions in ATP production and the TCA Cycle. Where predictions and experimental data disagreed, we identified testable hypotheses to reconcile differences between the model predictions and experimental data. The presented pipeline for using paired transcriptomics and metabolomics data provides a framework for interrogating multiple omics datasets to generate mechanistic insight of metabolic changes associated with toxicological responses.

Influence of lifestyle and genetic variants in the aldo-keto reductase 1C3 rs12529 polymorphism in high-risk prostate cancer detection variability assessed between US and New Zealand cohorts

Introduction

The prostate-specific antigen (PSA) based prostate cancer (PC) screening is currently being debated. The current assessment is to understand the variability of detecting high-risk PC in a NZ cohort in comparison to a US cohort with better PSA screening facilities. Aldo-keto reductase 1C3 (AKR1C3) is known for multiple functions with a potential to regulate subsequent PSA levels. Therefore, we wish to understand the influence of tobacco smoking and the AKR1C3 rs12529 gene polymorphism in this variability.

Method

NZ cohort (n = 376) consisted of 94% Caucasians while the US cohort consisted of African Americans (AA), n = 202, and European Americans (EA), n = 232. PSA level, PC grade and stage at diagnosis were collected from hospital databases for assigning high-risk PC status. Tobacco smoking status and the AKR1C3 rs12529 SNP genotype were considered as confounding variables. Variation of the cumulative % high-risk PC (outcome variable) with increasing PSA intervals (exposure factor) was compared between the cohorts using the Kolmogorov-Smirnov test. Comparisons were carried out with and without stratifications made using confounding variables.

Results

NZ cohort has been diagnosed at a significantly higher mean age (66.67± (8.08) y) compared to both AA (62.65±8.17y) and EA (64.83+8.56y); median PSA (NZ 8.90ng/ml compared to AA 6.86ng/ml and EA 5.80ng/ml); and Gleason sum (NZ (7) compared EA (6)) (p<0.05). The cumulative % high-risk PC detection shows NZ cohort with a significantly lower diagnosis rates at PSA levels between >6 - <10ng/ml compared to both US groups (p<0.05). These were further compounded significantly by smoking status and genetics.

Conclusions

High-risk PCs recorded at higher PSA levels in NZ could be due to factors including lower levels of PSA screening and subsequent specialist referrals for biopsies. These consequences could be pronounced among NZ ever smokers carrying the AKR1C3 rs12529 G alleles making them a group that requires increased PSA screening attention.

CES1 genetic variation affects the activation of angiotensin-converting enzyme inhibitors

The aim of the study was to determine the effect of carboxylesterase 1 (CES1) genetic variation on the activation of angiotensin-converting enzyme inhibitor (ACEI) prodrugs. In vitro incubation study of human liver, intestine and kidney s9 fractions demonstrated that the ACEI prodrugs enalapril, ramipril, perindopril, moexipril and fosinopril are selectively activated by CES1 in the liver. The impact of CES1/CES1VAR and CES1P1/CES1P1VAR genotypes and diplotypes on CES1 expression and activity on enalapril activation was investigated in 102 normal human liver samples. Neither the genotypes nor the diplotypes affected hepatic CES1 expression and activity. Moreover, among several CES1 nonsynonymous variants studied in transfected cell lines, the G143E (rs71647871) was a loss-of-function variant for the activation of all ACEIs tested. The CES1 activity on enalapril activation in human livers with the 143G/E genotype was approximately one-third of that carrying the 143G/G. Thus, some functional CES1 genetic variants (for example, G143E) may impair ACEI activation, and consequently affect therapeutic outcomes of ACEI prodrugs.

Pilot study of CYP2B6 genetic variation to explore the contribution of nitrosamine activation to lung carcinogenesis

We explored the contribution of nitrosamine metabolism to lung cancer in a pilot investigation of genetic variation in CYP2B6, a high-affinity enzymatic activator of tobacco-specific nitrosamines with a negligible role in nicotine metabolism. Previously we found that variation in CYP2A6 and CHRNA5-CHRNA3-CHRNB4 combined to increase lung cancer risk in a case-control study in European American ever-smokers (n = 860). However, these genes are involved in the pharmacology of both nicotine, through which they alter smoking behaviours, and carcinogenic nitrosamines. Herein, we separated participants by CYP2B6 genotype into a high- vs. low-risk group (*1/*1 + *1/*6 vs. *6/*6). Odds ratios estimated through logistic regression modeling were 1.25 (95% CI 0.68-2.30), 1.27 (95% CI 0.89-1.79) and 1.56 (95% CI 1.04-2.31) for CYP2B6, CYP2A6 and CHRNA5-CHRNA3-CHRNB4, respectively, with negligible differences when all genes were evaluated concurrently. Modeling the combined impact of high-risk genotypes yielded odds ratios that rose from 2.05 (95% CI 0.39-10.9) to 2.43 (95% CI 0.47-12.7) to 3.94 (95% CI 0.72-21.5) for those with 1, 2 and 3 vs. 0 high-risk genotypes, respectively. Findings from this pilot point to genetic variation in CYP2B6 as a lung cancer risk factor supporting a role for nitrosamine metabolic activation in the molecular mechanism of lung carcinogenesis.