Liver eQTL meta-analysis illuminates potential molecular mechanisms of cardiometabolic traits.

Understanding the molecular mechanisms of complex traits is essential for developing targeted interventions. We analyzed liver expression quantitative-trait locus (eQTL) meta-analysis data on 1,183 participants to identify conditionally distinct signals. We found 9,013 eQTL signals for 6,564 genes; 23% of eGenes had two signals, and 6% had three or more signals. We then integrated the eQTL results with data from 29 cardiometabolic genome-wide association study (GWAS) traits and identified 1,582 GWAS-eQTL colocalizations for 747 eGenes. Non-primary eQTL signals accounted for 17% of all colocalizations. Isolating signals by conditional analysis prior to coloc resulted in 37% more colocalizations than using marginal eQTL and GWAS data, highlighting the importance of signal isolation. Isolating signals also led to stronger evidence of colocalization: among 343 eQTL-GWAS signal pairs in multi-signal regions, analyses that isolated the signals of interest resulted in higher posterior probability of colocalization for 41% of tests. Leveraging allelic heterogeneity, we predicted causal effects of gene expression on liver traits for four genes. To predict functional variants and regulatory elements, we colocalized eQTL with liver chromatin accessibility QTL (caQTL) and found 391 colocalizations, including 73 with non-primary eQTL signals and 60 eQTL signals that colocalized with both a caQTL and a GWAS signal. Finally, we used publicly available massively parallel reporter assays in HepG2 to highlight 14 eQTL signals that include at least one expression-modulating variant. This multi-faceted approach to unraveling the genetic underpinnings of liver-related traits could lead to therapeutic development.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies.

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Genetic effects on liver chromatin accessibility identify disease regulatory variants.

Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.

Identification of Reduced ERAP2 Expression and a Novel HLA Allele as Components of a Risk Score for Susceptibility to Liver Injury Due to Amoxicillin-Clavulanate.

BACKGROUND & AIMS: Drug-induced liver injury (DILI) due to amoxicillin-clavulanate (AC) has been associated with HLA-A∗02:01, HLA-DRB1∗15:01, and rs2476601, a missense variant in PTPN22. The aim of this study was to identify novel risk factors for AC-DILI and to construct a genetic risk score (GRS). METHODS: Transcriptome-wide association study and genome-wide association study analyses were performed on 444 AC-DILI cases and 10,397 population-based controls of European descent. Associations were confirmed in a validation cohort (n = 133 cases and 17,836 population-based controls). Discovery and validation AC-DILI cases were also compared with 1358 and 403 non-AC-DILI cases. RESULTS: Transcriptome-wide association study revealed a significant association of AC-DILI risk with reduced liver expression of ERAP2 (P = 3.7 × 10-7), coding for an aminopeptidase involved in antigen presentation. The lead eQTL single nucleotide polymorphism, rs1363907 (G), was associated with AC-DILI risk in the discovery (odds ratio [OR], 1.68; 95% CI, 1.23-1.66; P = 1.7 × 10-7) and validation cohorts (OR, 1.2; 95% CI, 1.04-2.05; P = .03), following a recessive model. We also identified HLA-B∗15:18 as a novel AC-DILI risk factor in both discovery (OR, 4.19; 95% CI, 2.09-8.36; P = 4.9 × 10-5) and validation (OR, 7.78; 95% CI, 2.75-21.99; P = .0001) cohorts. GRS, incorporating rs1363907, rs2476601, HLA-B∗15:18, HLA-A∗02:01, and HLA-DRB1∗15:01, was highly predictive of AC-DILI risk when cases were analyzed against both general population and non-AC-DILI control cohorts. GRS was the most significant predictor in a regression model containing known AC-DILI clinical risk characteristics and significantly improved the predictive model. CONCLUSIONS: We identified novel associations of AC-DILI risk with ERAP2 low expression and with HLA-B∗15:18. GRS based on the 5 risk variants may assist AC-DILI causality assessment and risk management.

Sequencing of genes of drug response in tumor DNA and implications for precision medicine in cancer patients.

Tumor DNA sequencing is becoming standard-of-care for patient treatment decisions. We evaluated genotype concordance between tumor DNA and genomic DNA from blood and catalogued functional effects of somatic mutations in 21 drug response genes in 752 solid tumor patients. Using a threshold of 10% difference between tumor and blood DNA variant allele fraction (VAF), concordance for heterogenous genotype calls was 78% and increased to 97.5% using a 30% VAF threshold. Somatic mutations were observed in all 21 drug response genes, and 44% of patients had at least one somatic mutation in these genes. In tumor DNA, eight patients had a frameshift mutation in CYP2C8, which metabolizes taxanes. Overall, somatic copy number losses were more frequent than gains, including for CYP2C19 and CYP2D6 which had the most frequent copy number losses. However, copy number gains in TPMT were more than four times as common as losses. Seven % of patients had copy number gains in ABCB1, a multidrug resistance transporter of anti-cancer agents. These results demonstrate tumor-only DNA sequencing might not be reliable to call germline genotypes of drug response variants.

Genome-wide association studies of survival in 1520 cancer patients treated with bevacizumab-containing regimens.

Germline variants might predict cancer progression. Bevacizumab improves overall survival (OS) in patients with advanced cancers. No biomarkers are available to identify patients that benefit from bevacizumab. A meta-analysis of genome-wide association studies (GWAS) was conducted in 1,520 patients from Phase III trials (CALGB 80303, 40503, 80405 and ICON7), where bevacizumab was randomized to treatment without bevacizumab. We aimed to identify genes and single nucleotide polymorphisms (SNPs) associated with survival independently of bevacizumab treatment or through interaction with bevacizumab. A cause-specific Cox model was used to test the SNP-OS association in both arms combined (prognostic), and the effect of SNPs-bevacizumab interaction on OS (predictive) in each study. The SNP effects across studies were combined using inverse variance. Findings were tested for replication in advanced colorectal and ovarian cancer patients from The Cancer Genome Atlas (TGCA). In the GWAS meta-analysis, patients with rs680949 in PRUNE2 experienced shorter OS compared to patients without it (P = 1.02 × 10-7 , hazard ratio [HR] = 1.57, 95% confidence interval [CI] 1.33-1.86), as well as in TCGA (P = .0219, HR = 1.58, 95% CI 1.07-2.35). In the GWAS meta-analysis, patients with rs16852804 in BARD1 experienced shorter OS compared to patients without it (P = 1.40 × 10-5 , HR = 1.51, 95% CI 1.25-1.82) as well as in TCGA (P = 1.39 × 10-4 , HR = 3.09, 95% CI 1.73-5.51). Patients with rs3795897 in AGAP1 experienced shorter OS in the bevacizumab arm compared to the nonbevacizumab arm (P = 1.43 × 10-5 ). The largest GWAS meta-analysis of bevacizumab treated patients identified PRUNE2 and BARD1 (tumor suppressor genes) as prognostic genes of colorectal and ovarian cancer, respectively, and AGAP1 as a potentially predictive gene that interacts with bevacizumab with respect to patient survival.

Plasma levels of angiopoietin-2, VEGF-A, and VCAM-1 as markers of bevacizumab-induced hypertension: CALGB 80303 and 90401 (Alliance).

Hypertension is a common toxicity induced by bevacizumab and other antiangiogenic drugs. There are no biomarkers to predict the risk of bevacizumab-induced hypertension. This study aimed to identify plasma proteins related to the function of the vasculature to predict the risk of severe bevacizumab-induced hypertension. Using pretreated plasma samples from 398 bevacizumab-treated patients in two clinical trials (CALGB 80303 and 90401), the levels of 17 proteins were measured via ELISA. The association between proteins and grade 3 bevacizumab-induced hypertension was performed by calculating the odds ratio (OR) from logistic regression adjusting for age, sex, and clinical trial. Using the optimal cut-point of each protein, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) for hypertension were estimated. Five proteins showed no difference in levels between clinical trials and were used for analyses. Lower levels of angiopoietin-2 (p = 0.0013, OR 3.41, 95% CI 1.67-7.55), VEGF-A (p = 0.0008, OR 4.25, 95% CI 1.93-10.72), and VCAM-1 (p = 0.0067, OR 2.68, 95% CI 1.34-5.63) were associated with an increased risk of grade 3 hypertension. The multivariable model suggests independent effects of angiopoietin-2 (p = 0.0111, OR 2.71, 95% CI 1.29-6.10), VEGF-A (p = 0.0051, OR 3.66, 95% CI 1.54-9.73), and VCAM-1 (p = 0.0308, OR 2.27, 95% CI 1.10-4.92). The presence of low levels of 2-3 proteins had an OR of 10.06 (95% CI 3.92-34.18, p = 1.80 × 10-5) for the risk of hypertension, with sensitivity of 89.7%, specificity of 53.5%, PPV of 17.3%, and NPV of 97.9%. This is the first study providing evidence of plasma proteins with potential value to predict patients at risk of developing bevacizumab-induced hypertension.Clinical trial registration: ClinicalTrials.gov Identifier: NCT00088894 (CALGB 80303); and NCT00110214 (CALGB 90401).

Integration of DNA sequencing with population pharmacokinetics to improve the prediction of irinotecan exposure in cancer patients.

BACKGROUND: Irinotecan (CPT-11) is an anticancer agent widely used to treat adult solid tumours. Large interindividual variability in the clearance of irinotecan and SN-38, its active and toxic metabolite, results in highly unpredictable toxicity. METHODS: In 217 cancer patients treated with intravenous irinotecan single agent or in combination, germline DNA was used to interrogate the variation in 84 genes by next-generation sequencing. A stepwise analytical framework including a population pharmacokinetic model with SNP- and gene-based testing was used to identify demographic/clinical/genetic factors that influence the clearance of irinotecan and SN-38. RESULTS: Irinotecan clearance was influenced by rs4149057 in SLCO1B1, body surface area, and co-administration of 5-fluorouracil/leucovorin/bevacizumab. SN-38 clearance was influenced by rs887829 in UGT1A1, pre-treatment total bilirubin, and EGFR rare variant burden. Within each UGT1A1 genotype group, elevated pre-treatment total bilirubin and/or presence of at least one rare variant in EGFR resulted in significantly lower SN-38 clearance. The model reduced the interindividual variability in irinotecan clearance from 38 to 34% and SN-38 clearance from 49 to 32%. CONCLUSIONS: This new model significantly reduced the interindividual variability in the clearance of irinotecan and SN-38. New genetic factors of variability in clearance have been identified.

Bevacizumab-induced hypertension and proteinuria: a genome-wide study of more than 1000 patients.

BACKGROUND: Hypertension and proteinuria are common bevacizumab-induced toxicities. No validated biomarkers are available for identifying patients at risk of these toxicities. METHODS: A genome-wide association study (GWAS) meta-analysis was performed in 1039 bevacizumab-treated patients of European ancestry in four clinical trials (CALGB 40502, 40503, 80303, 90401). Grade ≥2 hypertension and proteinuria were recorded (CTCAE v.3.0). Single-nucleotide polymorphism (SNP)-toxicity associations were determined using a cause-specific Cox model adjusting for age and sex. RESULTS: The most significant SNP associated with hypertension with concordant effect in three out of the four studies (p-value <0.05 for each study) was rs6770663 (A > G) in KCNAB1, with the G allele increasing the risk of hypertension (p-value = 4.16 × 10-6). The effect of the G allele was replicated in ECOG-ACRIN E5103 in 582 patients (p-value = 0.005). The meta-analysis of all five studies for rs6770663 led to p-value = 7.73 × 10-8, close to genome-wide significance. The most significant SNP associated with proteinuria was rs339947 (C > A, between DNAH5 and TRIO), with the A allele increasing the risk of proteinuria (p-value = 1.58 × 10-7). CONCLUSIONS: The results from the largest study of bevacizumab toxicity provide new markers of drug safety for further evaluations. SNP in KCNAB1 validated in an independent dataset provides evidence toward its clinical applicability to predict bevacizumab-induced hypertension. ClinicalTrials.gov Identifier: NCT00785291 (CALGB 40502); NCT00601900 (CALGB 40503); NCT00088894 (CALGB 80303) and NCT00110214 (CALGB 90401).

PIK3R5 genetic predictors of hypertension induced by VEGF-pathway inhibitors.

No biomarkers are available to predict patients at risk of developing hypertension induced by VEGF-pathway inhibitors. This study aimed to identify predictive biomarkers of hypertension induced by these drugs using a discovery-replication approach. The discovery set included 140 sorafenib-treated patients (TARGET study) genotyped for 973 SNPs in 56 genes. The most statistically significant SNPs associated with grade ≥2 hypertension were tested for association with grade ≥2 hypertension in the replication set of a GWAS of 1039 bevacizumab-treated patients from four clinical trials (CALGB/Alliance). In the discovery set, rs444904 (G > A) in PIK3R5 was associated with an increased risk of sorafenib-induced hypertension (p = 0.006, OR = 3.88 95% CI 1.54-9.81). In the replication set, rs427554 (G > A) in PIK3R5 (in complete linkage disequilibrium with rs444904) was associated with an increased risk of bevacizumab-induced hypertension (p = 0.008, OR = 1.39, 95% CI 1.09-1.78). This study identified a predictive marker of drug-induced hypertension that should be evaluated for other VEGF-pathway inhibitors.ClinicalTrials.gov Identifier:NCT00073307 (TARGET).

KDR genetic predictor of toxicities induced by sorafenib and regorafenib.

No biomarkers are available to predict toxicities induced by VEGFR TKIs. This study aimed to identify markers of toxicities induced by these drugs using a discovery-validation approach. The discovery set included 140 sorafenib-treated cancer patients (TARGET study) genotyped for SNPs in 56 genes. The most significant SNPs associated with grade ≥2 hypertension, diarrhea, dermatologic toxicities, and composite toxicity (any one of the toxicities) were tested for association with grade ≥2 toxicity in a validation set of 201 sorafenib-treated patients (Alliance/CALGB 80802). The validated SNP was tested for association with grade ≥2 toxicity in 107 (LCCC 1029) and 82 (Italian cohort) regorafenib-treated patients. SNP-toxicity associations were evaluated using logistic regression, and a meta-analysis between the studies was performed by inverse variance. Variant rs4864950 in KDR increased the risk of grade ≥2 composite toxicity in TARGET, Alliance/CALGB 80802, and the Italian cohort (meta-analysis p = 6.79 × 10-4, OR = 2.01, 95% CI 1.34-3.01). We identified a predictor of toxicities induced by VEGFR TKIs. CLINICALTRIALS.GOV IDENTIFIER: NCT00073307 (TARGET), NCT01015833 (Alliance/CALGB 80802), and NCT01298570 (LCCC 1029).

An initial genetic analysis of gemcitabine-induced high-grade neutropenia in pancreatic cancer patients in CALGB 80303 (Alliance).

OBJECTIVES: One of the standard of care regimens for advanced pancreatic cancer is gemcitabine-based chemotherapy. The efficacy of gemcitabine is limited by dose-limiting hematologic toxicities especially neutropenia. Uncovering the variability of these toxicities attributed to germline DNA variation is of great importance. PATIENTS AND METHODS: CALGB 80303 was a randomized study in advanced pancreatic cancer patients treated with gemcitabine with or without bevacizumab. The study protocol included genotyping of genes of gemcitabine disposition (CDA, DCTD, SLC29A1, SLC28A1, and SLC29A2), as well as a genome-wide analysis. The clinical phenotype was time to early high-grade neutropenia event accounting for progression or death or other treatment-terminating adverse events as competing for informative events. The inference was carried out on the basis of the association between genotype and cause-specific hazard of a neutropenic event. RESULTS: The primary analyses were carried out on the basis of 294 genetically estimated European pancreatic cancer patients. For CDA rs2072671 (A>C), AC and CC patients had a lower risk of neutropenia than AA patients (P=0.01, hazard ratio: 0.61, 95% confidence interval: 0.41-0.89). For SLC28A1 rs3825876 (G>A), AA patients have a higher risk of neutropenia than GA and GG patients (P=0.02, hazard ratio: 1.51, 95% confidence interval: 1.06-2.16). CDA rs2072671 was associated with increased mRNA expression in whole blood in three studies (P=2.7e-14, 6.61e-62, and 9.70e-65). In the genome-wide analysis, variants in TGFB2 were among the top hits (lowest P=1.62e-06) but had no effect in luciferase assays. CONCLUSION: This is the first genetic analysis of gemcitabine-induced neutropenia using a competing risk model in a prospective randomized clinical study has proposed a potentially novel mechanism of the protective effect of the CDA rs2072671 variant. Further confirmation is needed.

Genomic Characterization of Metformin Hepatic Response.

Metformin is used as a first-line therapy for type 2 diabetes (T2D) and prescribed for numerous other diseases. However, its mechanism of action in the liver has yet to be characterized in a systematic manner. To comprehensively identify genes and regulatory elements associated with metformin treatment, we carried out RNA-seq and ChIP-seq (H3K27ac, H3K27me3) on primary human hepatocytes from the same donor treated with vehicle control, metformin or metformin and compound C, an AMP-activated protein kinase (AMPK) inhibitor (allowing to identify AMPK-independent pathways). We identified thousands of metformin responsive AMPK-dependent and AMPK-independent differentially expressed genes and regulatory elements. We functionally validated several elements for metformin-induced promoter and enhancer activity. These include an enhancer in an ataxia telangiectasia mutated (ATM) intron that has SNPs in linkage disequilibrium with a metformin treatment response GWAS lead SNP (rs11212617) that showed increased enhancer activity for the associated haplotype. Expression quantitative trait locus (eQTL) liver analysis and CRISPR activation suggest that this enhancer could be regulating ATM, which has a known role in AMPK activation, and potentially also EXPH5 and DDX10, its neighboring genes. Using ChIP-seq and siRNA knockdown, we further show that activating transcription factor 3 (ATF3), our top metformin upregulated AMPK-dependent gene, could have an important role in gluconeogenesis repression. Our findings provide a genome-wide representation of metformin hepatic response, highlight important sequences that could be associated with interindividual variability in glycemic response to metformin and identify novel T2D treatment candidates.

Discovery and functional assessment of gene variants in the vascular endothelial growth factor pathway.

Angiogenesis is a host-mediated mechanism in disease pathophysiology. The vascular endothelial growth factor (VEGF) pathway is a major determinant of angiogenesis, and a comprehensive annotation of the functional variation in this pathway is essential to understand the genetic basis of angiogenesis-related diseases. We assessed the allelic heterogeneity of gene expression, population specificity of cis expression quantitative trait loci (eQTLs), and eQTL function in luciferase assays in CEU and Yoruba people of Ibadan, Nigeria (YRI) HapMap lymphoblastoid cell lines in 23 resequenced genes. Among 356 cis-eQTLs, 155 and 174 were unique to CEU and YRI, respectively, and 27 were shared between CEU and YRI. Two cis-eQTLs provided mechanistic evidence for two genome-wide association study findings. Five eQTLs were tested for function in luciferase assays and the effect of two KRAS variants was concordant with the eQTL effect. Two eQTLs found in each of PRKCE, PIK3C2A, and MAP2K6 could predict 44%, 37%, and 45% of the variance in gene expression, respectively. This is the first analysis focusing on the pattern of functional genetic variation of the VEGF pathway genes in CEU and YRI populations and providing mechanistic evidence for genetic association studies of diseases for which angiogenesis plays a pathophysiologic role.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD).

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD).

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Metabolism and disposition of [14C]n-butyl-p-hydroxybenzoate in male and female Harlan Sprague Dawley rats following oral administration and dermal application.

n-Butyl-p-hydroxybenzoate (n-butylparaben, BPB) is an antioxidant used in foods, pharmaceuticals and cosmetics. This study investigated the disposition of ring-labelled [(14)C]BPB in Harlan Sprague Dawley rats, and in rat and human hepatocytes. BPB was rapidly cleared in hepatocytes from rat (t(1/2) = 3-4 min) and human (t(1/2) = 20-30 min). The major metabolites detected in rat hepatocytes were hydroxybenzoic acid and in human hepatocytes were hydroxybenzoic acid and hydroxyhippuric acid. [(14)C]BPB was administered to male rats orally at 10, 100 or 1000 mg/kg, intravenously at 10 mg/kg and dermally at 10 and 100 mg/kg; female rats were administered oral doses at 10 mg/kg. Oral doses of BPB were well-absorbed (>83%) and eliminated chiefly in urine (83-84%); ≤ 1% of the radioactivity remained in tissues at 24 h or 72 h after dosing. About 4% and 8%, respectively, of 100 mg/kg dermal doses were absorbed in 24 h and 72 h, and about 50% of a 10 mg/kg dose was absorbed in 72 h. Metabolites detected in urine included those previously reported, BPB-glucuronide, BPB-sulfate, hydroxybenzoic acid and hydroxyhippuric acid, but also novel metabolites arising from ring hydroxylation followed by glucuronidation and sulfation.

Metabolism and disposition of 2-methoxy-4-nitroaniline in male and female Harlan Sprague Dawley rats and B6C3F1/N mice.

The disposition of 2-Methoxy-4-nitroaniline (MNA) was investigated in male and female Harlan Sprague Dawley rats and B6C3F(1)/N mice following oral, intravenous, and dermal exposure to [(14)C]MNA at 2, 15, or 150 mg/kg. Clearance of MNA was investigated in male and female rat, mouse, and human hepatocytes. MNA was cleared slowly in hepatocytes from rat (t(1/2) = 152-424 min) and human (t(1/2) = 118-403 min) but faster in mouse (t(1/2)= 70-106 min). MNA was well-absorbed in rats and mice following oral administration and eliminated chiefly in urine (rats, 75-79%; mice, 55-68%) 72 h post dosing. Less than 1% of the radioactivity remained in tissues at 72 h. MNA was poorly absorbed following dermal application in rats (5.5%) and mice (10%) over 24 h. The major pathway of metabolism of MNA was via hydroxylation of the phenyl ring to form 6-hydroxy MNA; major metabolites detected were sulfate and glucuronide conjugates of 6-hydroxy MNA. Following oral administration, the percent of total radioactivity bound in tissues bound was highest in liver (43%) and red blood cells (30%), whereas the radioactivity bound to DNA was highest in cecum (160 pmol/mg DNA).