Dioxin-like compound exposures and DNA methylation in the Anniston Community Health Survey Phase II.

The Anniston Community Health Survey (ACHS-I) was initially conducted from 2005 to 2007 to assess polychlorinated biphenyl (PCB) exposures in Anniston, Alabama residents. In 2014, a follow-up study (ACHS-II) was conducted to measure the same PCBs as in ACHS-I and additional compounds e.g., polychlorinated dibenzo-p-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs), and dioxin-like non-ortho (cPCBs) substituted PCBs. In this epigenome-wide association study (EWAS), we examined the associations between PCDD, PCDF, and PCB exposures and DNA methylation. Whole blood DNA methylation was measured using Illumina EPIC arrays (n=292). We modeled lipid-adjusted toxic equivalencies (TEQs) for: ΣDioxins (sum of 28 PCDDs, PCDFs, cPCBs, and mPCBs), PCDDs, PCDFs, cPCBs, and mPCBs using robust multivariable linear regression adjusting for age, race, sex, smoking, bisulfite conversion batch, and estimated percentages of six blood cell types. Among all exposures we identified 10 genome-wide (Bonferroni p≤6.74E-08) and 116 FDR (p≤5.00E-02) significant associations representing 10 and 113 unique CpGs, respectively. Of the 10 genome-wide associations, seven (70%) occurred in the PCDDs and four (40%) of these associations had an absolute differential methylation ≥1.00%, based on the methylation difference between the highest and lowest exposure quartiles. Most of the associations (six, 60%) represented hypomethylation changes. Of the 10 unique CpGs, eight (80%) were in genes shown to be associated with dioxins and/or PCBs based on data from the 2019 Comparative Toxicogenomics Database. In this study, we have identified a set of CpGs in blood DNA that may be particularly susceptible to dioxin, furan, and dioxin-like PCB exposures.

Polychlorinated biphenyl exposure and DNA methylation in the Anniston Community Health Survey.

Anniston, Alabama was home to a major polychlorinated biphenyl (PCB) production facility from 1929 until 1971. The Anniston Community Health Survey I and II (ACHS-I 2005-2007, ACHS-II 2013-2014) were conducted to explore the effects of PCB exposures. In this report we examined associations between PCB exposure and DNA methylation in whole blood using EPIC arrays (ACHS-I, n = 518; ACHS-II, n = 299). For both cohorts, 35 PCBs were measured in serum. We modelled methylation versus PCB wet-weight concentrations for: the sum of 35 PCBs, mono-ortho substituted PCBs, di-ortho substituted PCBs, tri/tetra-ortho substituted PCBs, oestrogenic PCBs, and antiestrogenic PCBs. Using robust multivariable linear regression, we adjusted for age, race, sex, smoking, total lipids, and six blood cell-type percentages. We carried out a two-stage analysis; discovery in ACHS-I followed by replication in ACHS-II. In ACHS-I, we identified 28 associations (17 unique CpGs) at p ≤ 6.70E-08 and 369 associations (286 unique CpGs) at FDR p ≤ 5.00E-02. A large proportion of the genes have been observed to interact with PCBs or dioxins in model studies. Among the 28 genome-wide significant CpG/PCB associations, 14 displayed replicated directional effects in ACHS-II; however, only one in ACHS-II was statistically significant at p ≤ 1.70E-04. While we identified many novel CpGs significantly associated with PCB exposures in ACHS-I, the differential methylation was modest and the effect was attenuated seven years later in ACHS-II, suggesting a lack of persistence of the associations between PCB exposures and altered DNA methylation in blood cells.

Mitigation of the effect of variability in digital PCR assays through use of duplexed reference assays for normalization.

Digital PCR has been promoted as a technique for obtaining absolute measures of the amount of nucleic acid target sequence in a sample, but still lacks standardization in data reporting. The initial method of representing data as copies per microliter produced inconsistent results and made inter-assay comparisons difficult. Normalizing copies to amount of nucleic acid gives more uniform results, but factors influencing the effective concentration of nucleic acid in the final digital PCR assay must be considered. Using droplet digital PCR and previously validated reference genes duplexed with target genes, a method of normalization was developed to estimate the amount of input nucleic acid in individual assays, subsequently reporting the number of copies of target gene relative to this amount. Correcting for the actual amount of amplifiable nucleic acid present demonstrated a higher correlation between various dilutions of sample mRNA and allowed more accurate comparisons of digital PCR results.

Gene expression changes in immune response pathways following oral administration of tetrabromobisphenol A (TBBPA) in female Wistar Han rats.

Tetrabromobisphenol A (TBBPA) is a brominated flame retardant used globally at high volumes, primarily in the epoxy resin of circuit boards. It has been detected in the environment and in humans. The National Toxicology Program found that chronic oral TBBPA treatment of 250mg/kg and higher caused an increased incidence of uterine lesions in female Wistar Han rats. The present laboratory has previously reported changes in gene expression associated with estrogen homeostasis in liver and uterine tissue of adult female Wistar Han rats after five days of gavage with 250mg/kg of TBBPA. Microarray analysis of tissue from these same TBBPA-treated rats was performed to detect additional pathways perturbed by TBBPA. Microarray analysis of uterine tissue detected downregulation of genes in pathways of the immune response following TBBPA treatment. These results, along with validation of associated gene expression changes using droplet digital PCR, are reported here. Our findings suggest mechanisms that may be related to estrogen-mediated immunosuppression.

Disruption of estrogen homeostasis as a mechanism for uterine toxicity in Wistar Han rats treated with tetrabromobisphenol A.

Chronic oral treatment of tetrabromobisphenol A (TBBPA) to female Wistar Han rats resulted in increased incidence of cell proliferation at 250mg/kg and tumor formation in the uterus at higher doses. The present study was designed to test the hypothesis that disruption of estrogen homeostasis was a major mode-of-action for the observed effects. Biological changes were assessed in serum, liver, and the proximal (nearest the cervix) and distal (nearest the ovaries) sections of the uterine horn of Wistar Han rats 24h following administration of the last of five daily oral doses of 250mg/kg. Expression of genes associated with receptors, biosynthesis, and metabolism of estrogen was altered in the liver and uterus. TBBPA treatment also resulted in changes in expression of genes associated with cell division and growth. Changes were also observed in the concentration of thyroxine in serum and in expression of genes in the liver and uterus associated with thyroid hormone receptors. Differential expression of some genes was tissue-dependent or specific to tissue location in the uterus. The biological responses observed in the present study support the hypothesis that perturbation of estrogen homeostasis is a major mode-of-action for TBBPA-mediated cell proliferation and tumorigenesis previously observed in the uterus of TBBPA-treated Wistar Han rats.

RBCK1, an E3 ubiquitin ligase, interacts with and ubiquinates the human pregnane X receptor.

The pregnane X receptor (PXR, NR1I2) plays a pivotal role in the disposition and detoxification of numerous foreign and endogenous chemicals by increasing transcription of numerous target genes, including phase I and II drug-metabolizing enzymes and transporters. In the present study, yeast two-hybrid screening identified an E3 ubiquitin ligase, RBCK1 (Ring-B-box-coiled-coil protein interacting with protein kinase C-1), as a human pregnane X receptor (hPXR)-interacting protein. Coimmunoprecipitation studies confirmed the interaction between RBCK1 and hPXR when both were ectopically expressed in AD-293 cells. Domain mapping studies showed that the interaction between RBCK1 and hPXR involves all RBCK1 domains. We further demonstrate that RBCK1 ubiquitinates hPXR, and this may target hPXR for degradation by the ubiquitin-proteasome pathway. Simultaneous ectopic overexpression of RBCK1 and PXR decreased PXR levels in AD-293 cells, and this decrease was inhibited by the proteasomal inhibitor MG-132 (carbobenzoxy-Leu-Leu-leucinal). Furthermore, overexpression of RBCK1 decreased endogenous levels of PXR in HepG2 cells. Of importance, ectopic overexpression and silencing of endogenous RBCK1 in primary human hepatocytes resulted in a decrease and increase, respectively, in endogenous PXR protein levels and in the induction of PXR target genes by rifampicin. These results suggest that RBCK1 is important for the ubiquitination of PXR and may play a role in its proteasomal degradation.

Human CYP2C8 is post-transcriptionally regulated by microRNAs 103 and 107 in human liver.

The CYP2C genes are extensively regulated at the transcriptional stage. The present study shows for the first time that CYP2Cs are also regulated post-transcriptionally by microRNAs (miRNAs). By using online search engines, we found potential miRNA response elements (MREs) in the 3'-untranslated region (3'-UTR) of the CYP2C mRNAs. Among these were a MRE for the miRNAs miR-103 and miR-107 in the 3'-UTR of human CYP2C8. CYP2C8 protein levels (measured through immunoblot analyses) did not correlate with CYP2C8 mRNA levels (measured through quantitative polymerase chain reaction analyses) in human liver samples. The translation efficiency (protein/mRNA ratio) for CYP2C8 was inversely correlated with the expression of miR-103 and miR-107. When three copies of the putative MRE from CYP2C8 were inserted downstream from a luciferase expression reporter, transfection with precursors for miR-103 or miR-107 decreased luciferase activity in primary hepatocytes, whereas transfection with antisense oligonucleotides (AsOs) for miR-103/miR-107 increased luciferase activity. As expected, there was no effect of the precursors or AsOs when three copies of the putative MRE were inserted in the reverse orientation. When precursors for miR-103/miR-107 were transfected into primary human hepatocytes, CYP2C8 protein levels were decreased, whereas AsOs increased CYP2C8 protein levels. Neither precursors nor AsOs affected CYP2C8 mRNA levels, which indicated that the effect was post-transcriptional. Putative MRE motifs were also found in the 3'-UTRs of CYP2C9 and CYP2C19, which suggested that the same miRNAs could regulate translation of other members of the CYP2C family, although to a lesser degree than CYP2C8. These results clearly show that CYP2Cs are regulated post-transcriptionally by miR-103 and miR-107.

Identification of human CYP2C8 as a retinoid-related orphan nuclear receptor target gene.

Retinoid-related orphan nuclear receptors (RORs) alpha and gamma (NR1F1, -3) are highly expressed in liver, adipose tissue, thymus, and brain and are involved in many physiological processes, such as circadian rhythm and immune function. Enzymes in the cytochrome P450 2C subfamily metabolize many clinically important drugs and endogenous compounds, such as the anticancer drug paclitaxel and arachidonic acid, and are highly expressed in liver. Here, we present the first evidence that RORs regulate the transcription of human CYP2C8. Overexpression of RORalpha and RORgamma in HepG2 cells significantly enhanced the activity of the CYP2C8 promoter but not that of the CYP2C9 or CYP2C19 promoters. Computer analyses, promoter deletion studies, gel shift assays, and mutational analysis identified an essential ROR-responsive element at -2045 base pairs in the CYP2C8 promoter that mediates ROR transactivation. Adenoviral overexpression of RORalpha and -gamma significantly induced endogenous CYP2C8 transcripts in both HepG2 cells and human primary hepatocytes. Knockdown of endogenous RORalpha and -gamma expression in HepG2 cells by RNA interference decreased the expression of endogenous CYP2C8 mRNA by approximately 50%. These data indicate that RORs transcriptionally up-regulate CYP2C8 in human liver and, therefore, may be important modulators of the metabolism of drugs and physiologically active endogenous compounds by this enzyme in liver and possibly extrahepatic tissues where RORs are expressed.

The discovery of new coding alleles of human CYP26A1 that are potentially defective in the metabolism of all-trans retinoic acid and their assessment in a recombinant cDNA expression system.

OBJECTIVES: Retinoic acid (RA) is a critical regulator of gene expression during embryonic development and in the maintenance of adult epithelial tissues. This study was undertaken to identify genetic polymorphisms of CYP26A1 which might affect these processes. We sequenced CYP26A1 in racially diverse individuals and assessed the metabolism of retinoic acid by newly identified coding alleles of CYP26A1 in a recombinant system. METHODS: CYP26A1 was sequenced in 24 Caucasians, 24 African-Americans, 24 Asians, and 20 individuals of unknown racial origin. cDNA constructs for wild-type and coding alleles of CYP26A1 were constructed in a pcDNA3.1 expression vector and expressed in Cos-1 cells. A FLAG tag at the C-terminal end of the cDNA was used to quantitate the recombinant CYP26A1 proteins. RESULTS: A total of 13 single nucleotide polymorphisms (SNPs) were identified in CYP26A1. Three SNPs produced coding changes: R173S, F186L, and C358R. These alleles were termed as CYP26A1*2, CYP26A1*3, and CYP26A1*4, respectively, by the Human Cytochrome P450 (CYP) Allele Nomenclature Committee at http://www.cypalleles.ki.se/. Wild type CYP26A1 protein metabolized all-trans-retinoic acid (at-RA) to 4-oxo-RA, 4-OH-RA as well as water-soluble metabolites. CYP26A1.3 (F186L) and CYP26A1.4 (C358R) allelic proteins exhibited significantly lower metabolism (40-80%) of at-RA than wild-type CYP26A1.1 protein. CONCLUSION: This is the first study to identify coding alleles of CYP26A1. Two coding alleles, CYP26A1*3 and CYP26A1*4, are predicted to be defective based on the metabolism of at-RA by the recombinant proteins. These studies suggest the need for future clinical studies of polymorphisms of CYP26A1 in embryonic development.

Detection of human CYP2C8, CYP2C9, and CYP2J2 in cardiovascular tissues.

The cytochrome P450 (P450) enzymes CYP2C8, CYP2C9, and CYP2J2 metabolize arachidonic acid to epoxyeicosatrienoic acids, which are known to be vital in regulation of vascular tone and cardiovascular homeostasis. Because there is limited information regarding the relative expression of these P450 enzymes in cardiovascular tissues, this study examined the expression of CYP2C8, CYP2C9, and CYP2J2 mRNA and protein in human heart, aorta, and coronary artery samples by real-time polymerase chain reaction, immunoblotting, and immunohistochemistry. CYP2J2 and CYP2C9 mRNA levels were highly variable in human hearts, whereas CYP2C8 mRNA was present in lower abundance. CYP2J2 mRNA was approximately 10(3) times higher than CYP2C9 or CYP2C8 in human heart. However, CYP2C9 mRNA was more abundant than CYP2J2 or CYP2C8 in one ischemic heart. In human aorta, mean CYP2C9 mRNA levels were approximately 50 times higher than that of CYP2J2 and 5-fold higher than that of CYP2C8. In human coronary artery, mean values for CYP2C9 mRNA were approximately 2-fold higher than that of CYP2J2 mRNA and 6-fold higher than that of CYP2C8 mRNA. Immunoblotting results show relatively high levels of CYP2J2 and CYP2C8 protein in human hearts, which was confirmed by immunohistochemistry. CYP2C9 protein was also detected at high levels in one ischemic heart by immunoblotting. CYP2C9 was present at higher levels than CYPJ2 in aorta and coronary artery, whereas CYP2C8 protein was below the limits of detection. The expression of CYP2J2 and CYP2C8 in human heart, and CYPC9 and CYP2J2 in aorta and coronary artery is consistent with a physiological role for these enzymes in these tissues.

Functional characterization of novel allelic variants of CYP2C9 recently discovered in southeast Asians.

CYP2C9 was recently resequenced in 150 Asian subjects from Singapore. Several new coding variants were reported, and these variants are now named CYP2C9*14 (R125H), CYP2C9*15 (S162X), CYP2C9*16 (T299A), CYP2C9*17 (P382S), CYP2C9*18 (D397A), and CYP2C9*19 (Q454H). The CYP2C9*18 variant also contained an I359L change previously associated with the CYP2C9*3 allele. In this study, we assessed the functional consequences of the new coding changes. cDNAs containing each of the new coding changes were constructed by site-directed mutagenesis and expressed in a bacterial cDNA expression system, the allelic proteins were partially purified, and their ability to hydroxylate a prototype CYP2C9 substrate was assayed. Expression of cDNAs in Escherichia coli containing either the D397A change or the S162X (premature stop codon) could not be detected either spectrally or at the apoprotein level. CYP2C9.14 and CYP2C9.16 exhibited 80 to 90% lower catalytic activity toward tolbutamide at two substrate concentrations compared with wild-type CYP2C9.1. Kinetic analysis confirmed that CYP2C9.14 and CYP2C9.16 have a higher Km and a >90% lower intrinsic clearance of tolbutamide compared with wild-type CYP2C9.1. Both CYP2C9.17 and CYP2C9.19 proteins exhibited modest 30 to 40% decreases in catalytic activity toward tolbutamide. Thus, CYP2C9*15 and CYP2C9*18 may represent null alleles, whereas CYP2C9*14 and CYP2C9*16 allelic variants produce proteins that are clearly catalytically defective in vitro, indicating the existence of new defective putative alleles of CYP2C9 in Asians.

Recombinant CYP3A4*17 is defective in metabolizing the hypertensive drug nifedipine, and the CYP3A4*17 allele may occur on the same chromosome as CYP3A5*3, representing a new putative defective CYP3A haplotype.

Genetic variation in CYP3A activity may influence the rate of the metabolism and elimination of CYP3A substrates in humans. We previously reported four new CYP3A4 coding variants in three different racial groups. In the present study, we examined metabolism of nifedipine by the recombinant forms of these allelic variants. Metabolism of nifedipine by the L293P (CYP3A4*18), M445T (CYP3A4*3), and P467S (CYP3A4*19) allelic variants was not significantly different from wild-type CYP3A4*1. However, F189S (CYP3A4*17) exhibited a >99% decrease in both V(max) and CL(max) of nifedipine compared with CYP3A4*1. Of 72 racially diverse individuals, CYP3A4*17 was identified in 1 of 24 Caucasian samples [1:5 Eastern European (Adygei ethnic group)]. Genotyping of an extended set of 276 genomic DNAs of Caucasians (100 from the Coriell Repository and an additional 176 from the United States) for CYP3A4*17 detected no additional individuals containing the CYP3A4*17 allele. However, additional genotyping of four more Adygei samples available from Coriell detected an additional individual carrying the CYP3A4*17 allele. New specific polymerase chain reaction-restriction fragment length polymorphism genotyping procedures were developed for the major splice variant of CYP3A5 (CYP3A5*3) and CYP3A4*17. Genotyping revealed that the two individuals carrying CYP3A4*17 were either homozygous or heterozygous for the more frequent CYP3A5*3 allele, suggesting that the two alleles may exist on the same chromosome as a new putative CYP3A poor metabolizer haplotype. We predict that individuals who are homozygous for defective alleles of both of these genes would metabolize CYP3A substrates poorly. The new genetic tests will be useful in future clinical studies to investigate genotype/phenotype associations.

Discovery of new potentially defective alleles of human CYP2C9.

CYP2C9 is a clinically important enzyme, responsible for the metabolism of numerous clinically important therapeutic drugs. In the present study, we discovered 38 single nucleotide polymorphisms in CYP2C9 by resequencing of genomic DNA from 92 individuals from three different racial groups. Haplotype analysis predicted that there are at least 21 alleles of CYP2C9 in this group of individuals. Six new alleles were identified that contained coding changes: L19I (CYP2C9*7), R150H (CYP2C9*8), H251R (CYP2C9*9), E272G (CYP2C9*10), R335W(CYP2C9*11) and P489S (CYP2C9*12). When expressed in a bacterial cDNA expression system, several alleles exhibited altered catalytic activity. CYP2C9*11 appeared to be a putative poor metabolizer allele, exhibiting a three-fold increase in the Km and more than a two-fold decrease in the intrinsic clearance for tolbutamide. Examination of the crystal structure of human CYP2C9 reveals that R335 is located in the turn between the J and J' helices and forms a hydrogen-bonding ion pair with D341 from the J' helix. Abolishing this interaction in CYP2C9*11 individuals could destabilize the secondary structure and alter the substrate affinity. This new putative poor metabolizer (PM) allele was found in Africans. A second potentially PM allele CYP2C9*12 found in a racially unidentified sample also exhibited a modest decrease in the Vmax and the intrinsic clearance for tolbutamide in a recombinant system. Further clinical studies are needed to determine the effect of these new polymorphisms on the metabolism of CYP2C9 substrates.

CYP2C44, a new murine CYP2C that metabolizes arachidonic acid to unique stereospecific products.

The human CYP2Cs have been studied extensively with respect to the metabolism of clinically important drugs and endogenous chemicals such as arachidonic acid (AA). Five members of the mouse CYP2C family have previously been described that metabolize arachidonic acid into regio- and stereospecific epoxyeicosatrienoic acids (EETs) and hydroxyeicosatetraenoic acids, which have many important physiological roles. Herein, we describe the cloning and characterization of a new mouse cytochrome P450 (P450), CYP2C44, which has the lowest homology with other known mouse CYP2Cs. Western blotting and real-time polymerase chain reaction detected CYP2C44 mRNA and protein in liver >> kidney > adrenals. Kidney contained approximately 10% of the CYP2C44 mRNA content of liver. CYP2C44 metabolized AA to unique stereospecific products, 11R,12S-EET and 8R, 9S-EET, which are similar to those produced by rat CYP2C23. CY2C23 is highly expressed in rat kidney and has been suggested to be important in producing compensatory renal artery vasodilation in response to salt-loading in this species. Immunohistochemistry showed the presence of CYP2C44 in hepatocytes, biliary cells of the liver, and the proximal tubules of the kidney. Unlike mouse CYP2C29, CYP2C38, and CYP2C39, CYP2C44 did not metabolize the common CYP2C substrate tolbutamide. CYP2C44 was not induced by phenobarbital or pregnenolone-16alpha-carbonitrile, two prototypical inducers of hepatic P450s. The presence of CYP2C44 in mouse liver, kidney, and adrenals and the unique stereospecificity of its arachidonic acid metabolites are consistent with the possibility that it may have unique physiological roles within these tissues, such as modulation of electrolyte transport or vascular tone.

Identification and functional characterization of new potentially defective alleles of human CYP2C19.

CYP2C19 is a clinically important enzyme responsible for the metabolism of a number of therapeutic drugs, such as mephenytoin, omeprazole, diazepam, proguanil, propranolol and certain antidepressants. Genetic polymorphisms in this enzyme result in poor metabolizers of these drugs. There are racial differences in the incidence of the poor metabolizer trait, which represents 13-23% of Asians but only 3-5% of Caucasians. In this study, single nucleotide polymorphisms (SNPs) in CYP2C19 were identified by direct sequencing of genomic DNA from 92 individuals from three different racial groups of varied ethnic background, including Caucasians, Asians and blacks. Several new alleles were identified containing the coding changes Arg114 His (CYP2C19*9), Pro227 Leu (CYP2C19*10), Arg150 His (CYP2C19*11), stop491 Cys (CYP2C19*12), Arg410 Cys (CYP2C19*13), Leu17 Pro (CYP2C19*14) and Ile19 Leu (CYP2C19*15). When expressed in a bacterial cDNA expression system, CYP2C19*9 exhibited a modest decrease in the V(max) for 4'-hydroxylation of -mephenytoin, and no alteration in its affinity for reductase. CYP2C19*10 exhibited a dramatically higher K(m) and lower V(max) for mephenytoin. CYP2C19*12was unstable and expressed poorly in a bacterial cDNA expression system. Clinical studies will be required to confirm whether this allele is defective in vivo. CYP2C19*9, CYP2C19*10 and CYP2C19*12 all occurred in African-Americans, or individuals of African descent, and represent new potentially defective alleles of CYP2C19 which are predicted to alter risk of these populations to clinically important drugs.