A MYC inhibitor selectively alters the MYC and MAX cistromes and modulates the epigenomic landscape to regulate target gene expression.

MYC regulates multiple gene programs, raising questions about the potential selectivity and downstream transcriptional consequences of MYC inhibitors as cancer therapeutics. Here, we examined the effect of a small-molecule MYC inhibitor, MYCi975, on the MYC/MAX cistromes, epigenome, transcriptome, and tumorigenesis. Integrating these data revealed three major classes of MYCi975-modulated gene targets: type 1 (down-regulated), type 2 (up-regulated), and type 3 (unaltered). While cell cycle and signal transduction pathways were heavily targeted by MYCi, RNA biogenesis and core transcriptional pathway genes were spared. MYCi975 altered chromatin binding of MYC and the MYC network family proteins, and chromatin accessibility and H3K27 acetylation alterations revealed MYCi975 suppression of MYC-regulated lineage factors AR/ARv7, FOXA1, and FOXM1. Consequently, MYCi975 synergistically sensitized resistant prostate cancer cells to enzalutamide and estrogen receptor-positive breast cancer cells to 4-hydroxytamoxifen. Our results demonstrate that MYCi975 selectively inhibits MYC target gene expression and provide a mechanistic rationale for potential combination therapies.

Epigenomic and enhancer dysregulation in uterine leiomyomas.

BACKGROUND: Uterine leiomyomas, also known as uterine fibroids or myomas, are the most common benign gynecological tumors and are found in women of reproductive and postmenopausal age. There is an exceptionally high prevalence of this tumor in women by the age of 50 years. Black women are particularly affected, with an increased incidence, earlier age of onset, larger and faster growing fibroids and greater severity of symptoms as compared to White women. Although advances in identifying genetic and environmental factors to delineate these fibroids have already been made, only recently has the role of epigenomics in the pathogenesis of this disease been considered. OBJECTIVE AND RATIONALE: Over recent years, studies have identified multiple epigenomic aberrations that may contribute to leiomyoma development and growth. This review will focus on the most recent discoveries in three categories of epigenomic changes found in uterine fibroids, namely aberrant DNA methylation, histone tail modifications and histone variant exchange, and their translation into altered target gene architecture and transcriptional outcome. The findings demonstrating how the altered 3D shape of the enhancer can regulate gene expression from millions of base pairs away will be discussed. Additionally, translational implications of these discoveries and potential roadblocks in leiomyoma treatment will be addressed. SEARCH METHODS: A comprehensive PubMed search was performed to identify published articles containing keywords relevant to the focus of the review, such as: uterine leiomyoma, uterine fibroids, epigenetic alterations, epigenomics, stem cells, chromatin modifications, extracellular matrix [ECM] organization, DNA methylation, enhancer, histone post-translational modifications and dysregulated gene expression. Articles until September 2021 were explored and evaluated to identify relevant updates in the field. Most of the articles focused on in the discussion were published between 2015 and 2021, although some key discoveries made before 2015 were included for background information and foundational purposes. We apologize to the authors whose work was not included because of space restrictions or inadvertent omission. OUTCOMES: Chemical alterations to the DNA structure and of nucleosomal histones, without changing the underlying DNA sequence, have now been implicated in the phenotypic manifestation of uterine leiomyomas. Genome-wide DNA methylation analysis has revealed subsets of either suppressed or overexpressed genes accompanied by aberrant promoter methylation. Furthermore, differential promoter access resulting from altered 3D chromatin structure and histone modifications plays a role in regulating transcription of key genes thought to be involved in leiomyoma etiology. The dysregulated genes function in tumor suppression, apoptosis, angiogenesis, ECM formation, a variety of cancer-related signaling pathways and stem cell differentiation. Aberrant DNA methylation or histone modification is also observed in altering enhancer architecture, which leads to changes in enhancer-promoter contact strength, producing novel explanations for the overexpression of high mobility group AT-hook 2 and gene dysregulation found in mediator complex subunit 12 mutant fibroids. While many molecular mechanisms and epigenomic features have been investigated, the basis for the racial disparity observed among those in the Black population remains unclear. WIDER IMPLICATIONS: A comprehensive understanding of the exact pathogenesis of uterine leiomyoma is lacking and requires attention as it can provide clues for prevention and viable non-surgical treatment. These findings will widen our knowledge of the role epigenomics plays in the mechanisms related to uterine leiomyoma development and highlight novel approaches for the prevention and identification of epigenome targets for long-term non-invasive treatment options of this significantly common disease.

The long noncoding RNA H19 regulates tumor plasticity in neuroendocrine prostate cancer.

Neuroendocrine (NE) prostate cancer (NEPC) is a lethal subtype of castration-resistant prostate cancer (PCa) arising either de novo or from transdifferentiated prostate adenocarcinoma following androgen deprivation therapy (ADT). Extensive computational analysis has identified a high degree of association between the long noncoding RNA (lncRNA) H19 and NEPC, with the longest isoform highly expressed in NEPC. H19 regulates PCa lineage plasticity by driving a bidirectional cell identity of NE phenotype (H19 overexpression) or luminal phenotype (H19 knockdown). It contributes to treatment resistance, with the knockdown of H19 re-sensitizing PCa to ADT. It is also essential for the proliferation and invasion of NEPC. H19 levels are negatively regulated by androgen signaling via androgen receptor (AR). When androgen is absent SOX2 levels increase, driving H19 transcription and facilitating transdifferentiation. H19 facilitates the PRC2 complex in regulating methylation changes at H3K27me3/H3K4me3 histone sites of AR-driven and NEPC-related genes. Additionally, this lncRNA induces alterations in genome-wide DNA methylation on CpG sites, further regulating genes associated with the NEPC phenotype. Our clinical data identify H19 as a candidate diagnostic marker and predictive marker of NEPC with elevated H19 levels associated with an increased probability of biochemical recurrence and metastatic disease in patients receiving ADT. Here we report H19 as an early upstream regulator of cell fate, plasticity, and treatment resistance in NEPC that can reverse/transform cells to a treatable form of PCa once therapeutically deactivated.

Epigenomic tensor predicts disease subtypes and reveals constrained tumor evolution.

Understanding the epigenomic evolution and specificity of disease subtypes from complex patient data remains a major biomedical problem. We here present DeCET (decomposition and classification of epigenomic tensors), an integrative computational approach for simultaneously analyzing hierarchical heterogeneous data, to identify robust epigenomic differences among tissue types, differentiation states, and disease subtypes. Applying DeCET to our own data from 21 uterine benign tumor (leiomyoma) patients identifies distinct epigenomic features discriminating normal myometrium and leiomyoma subtypes. Leiomyomas possess preponderant alterations in distal enhancers and long-range histone modifications confined to chromatin contact domains that constrain the evolution of pathological epigenomes. Moreover, we demonstrate the power and advantage of DeCET on multiple publicly available epigenomic datasets representing different cancers and cellular states. Epigenomic features extracted by DeCET can thus help improve our understanding of disease states, cellular development, and differentiation, thereby facilitating future therapeutic, diagnostic, and prognostic strategies.

Altered chromatin landscape and enhancer engagement underlie transcriptional dysregulation in MED12 mutant uterine leiomyomas.

Uterine leiomyomas (fibroids) are a major source of gynecologic morbidity in reproductive age women and are characterized by the excessive deposition of a disorganized extracellular matrix, resulting in rigid benign tumors. Although down regulation of the transcription factor AP-1 is highly prevalent in leiomyomas, the functional consequence of AP-1 loss on gene transcription in uterine fibroids remains poorly understood. Using high-resolution ChIP-sequencing, promoter capture Hi-C, and RNA-sequencing of matched normal and leiomyoma tissues, here we show that modified enhancer architecture is a major driver of transcriptional dysregulation in MED12 mutant uterine leiomyomas. Furthermore, modifications in enhancer architecture are driven by the depletion of AP-1 occupancy on chromatin. Silencing of AP-1 subunits in primary myometrium cells leads to transcriptional dysregulation of extracellular matrix associated genes and partly recapitulates transcriptional and epigenetic changes observed in leiomyomas. These findings establish AP-1 driven aberrant enhancer regulation as an important mechanism of leiomyoma disease pathogenesis.

Interferon-γ signaling is associated with BRCA1 loss-of-function mutations in high grade serous ovarian cancer.

Loss-of-function mutations of the breast cancer type 1 susceptibility protein (BRCA1) are associated with breast (BC) and ovarian cancer (OC). To identify gene signatures regulated by epigenetic mechanisms in OC cells carrying BRCA1 mutations, we assessed cellular responses to epigenome modifiers and performed genome-wide RNA- and chromatin immunoprecipitation-sequencing in isogenic OC cells UWB1.289 (carrying a BRCA1 mutation, BRCA1-null) and UWB1.289 transduced with wild-type BRCA1 (BRCA1+). Increased sensitivity to histone deacetylase inhibitors (HDACi) was observed in BRCA1-null vs. BRCA1+ cells. Gene expression profiles of BRCA1-null vs. BRCA1+ cells and treated with HDACi were integrated with chromatin mapping of histone H3 lysine 9 or 27 acetylation. Gene networks activated in BRCA1-null vs. BRCA1 + OC cells related to cellular movement, cellular development, cellular growth and proliferation, and activated upstream regulators included TGFß1, TNF, and IFN-γ. The IFN-γ pathway was altered by HDACi in BRCA1+ vs. BRCA1-null cells, and in BRCA1-mutated/or low vs. BRCA1-normal OC tumors profiled in the TCGA. Key IFN-γ-induced genes upregulated at baseline in BRCA1-null vs. BRCA1+OC and BC cells included CXCL10, CXCL11, and IFI16. Increased localization of STAT1 in the promoters of these genes occurred in BRCA1-null OC cells, resulting in diminished responses to IFN-γ or to STAT1 knockdown. The IFN-γ signature was associated with improved survival among OC patients profiled in the TCGA. In all, our results support that changes affecting IFN-γ responses are associated with inactivating BRCA1 mutations in OC. This signature may contribute to altered responses to anti-tumor immunity in BRCA1-mutated cells or tumors.

Small-Molecule MYC Inhibitors Suppress Tumor Growth and Enhance Immunotherapy.

Small molecules that directly target MYC and are also well tolerated in vivo will provide invaluable chemical probes and potential anti-cancer therapeutic agents. We developed a series of small-molecule MYC inhibitors that engage MYC inside cells, disrupt MYC/MAX dimers, and impair MYC-driven gene expression. The compounds enhance MYC phosphorylation on threonine-58, consequently increasing proteasome-mediated MYC degradation. The initial lead, MYC inhibitor 361 (MYCi361), suppressed in vivo tumor growth in mice, increased tumor immune cell infiltration, upregulated PD-L1 on tumors, and sensitized tumors to anti-PD1 immunotherapy. However, 361 demonstrated a narrow therapeutic index. An improved analog, MYCi975 showed better tolerability. These findings suggest the potential of small-molecule MYC inhibitors as chemical probes and possible anti-cancer therapeutic agents.

Decreased expression of microRNA-29 family in leiomyoma contributes to increased major fibrillar collagen production.

OBJECTIVE: To determine the expression and function of the microRNA-29 family (miRNA-29a, miRNA-29b, miRNA-29c) in human leiomyoma and myometrium. DESIGN: Basic science experimental design. SETTING: Academic medical center. PATIENT(S): Women undergoing surgery for symptomatic uterine fibroids. INTERVENTION(S): Overexpression and knockdown of miRNA-29a, miRNA-29b, and miRNA-29c in primary leiomyoma and myometrial cells. MAIN OUTCOME MEASURE(S): [1] Expression of the miRNA-29 family members in vivo in leiomyoma versus myometrium; [2] Major fibrillar collagen (I, II, III) expression in leiomyoma and myometrial cells with manipulation of miRNA-29 species. RESULT(S): Members of the miRNA-29 family (29a, 29b, 29c) are all down-regulated in leiomyoma versus myometrium in vivo. The expression of the miRNA-29 family can be successfully modulated in primary leiomyoma and myometrial cells. Overexpression of the miRNA-29 family in leiomyoma cells results in down-regulation of the major fibrillar collagens. Down-regulation of the miRNA-29 species in myometrium results in an increase in collagen type III deposition. CONCLUSION(S): The miRNA-29 family is consistently down-regulated in leiomyoma compared to matched myometrial tissue. This down-regulation contributes to the increased collagen seen in leiomyomas versus myometrium. When miRNA-29 members are overexpressed in leiomyoma cells, protein levels of all of the major fibrillar collagens decrease. The miRNA-29 members are potential therapeutic targets in this highly prevalent condition.

Genomic Determinants of THAP11/ZNF143/HCFC1 Complex Recruitment to Chromatin.

The THAP11 and ZNF143 transcription factors recognize overlapping DNA sequences and are reported to exhibit signs of both competitive and cooperative binding. HCFC1 serves as a scaffold protein, bridging interactions between transcription factors, including THAP11 and ZNF143, and transcriptional coregulators. The exact mechanism of how DNA sequences guide the recruitment of the THAP11/ZNF143/HCFC1 complex to chromatin is still controversial. In this study, we use chromosomally integrated synthetic constructs and clustered regularly interspaced short palindromic repeat (CRISPR)-Cas9-mediated approaches in intact cells to elucidate the role of the DNA sequence in the recruitment of this complex and to establish its biological relevance. We show that the ACTACA submotif, shared by both THAP11 and ZNF143, directs the recruitment of THAP11 and HCFC1 to ZNF143-occupied loci. Importantly, its position, spacing, and orientation relative to the ZNF143 core motif are critical for this action. CRISPR-Cas9-mediated alterations of the ACTACA submotif at endogenous promoters recapitulated results obtained with synthetic constructs and resulted in altered gene transcription and histone modifications at targeted promoters. Our in vivo approaches provide strong evidence for the molecular role of the ACTACA submotif in THAP11, ZNF143, and HCFC1 cooperative recruitment to chromatin and its biological role in target gene expression.

Host cell factor-1 recruitment to E2F-bound and cell-cycle-control genes is mediated by THAP11 and ZNF143.

Host cell factor-1 (HCF-1) is a metazoan transcriptional coregulator essential for cell-cycle progression and cell proliferation. Current models suggest a mechanism whereby HCF-1 functions as a direct coregulator of E2F proteins, facilitating the expression of genes necessary for cell proliferation. In this report, we show that HCF-1 recruitment to numerous E2F-bound promoters is mediated by the concerted action of zinc finger transcription factors THAP11 and ZNF143, rather than E2F proteins directly. THAP11, ZNF143, and HCF-1 form a mutually dependent complex on chromatin, which is independent of E2F occupancy. Disruption of the THAP11/ZNF143/HCF-1 complex results in altered expression of cell-cycle control genes and leads to reduced cell proliferation, cell-cycle progression, and cell viability. These data establish a model in which a THAP11/ZNF143/HCF-1 complex is a critical component of the transcriptional regulatory network governing cell proliferation.

A role for WDR5 in integrating threonine 11 phosphorylation to lysine 4 methylation on histone H3 during androgen signaling and in prostate cancer.

Upon androgen stimulation, PKN1-mediated histone H3 threonine 11 phosphorylation (H3T11P) promotes AR target gene activation. However, the underlying mechanism is not completely understood. Here, we show that WDR5, a subunit of the SET1/MLL complex, interacts with H3T11P, and this interaction facilitates the recruitment of the MLL1 complex and subsequent H3K4 tri-methylation (H3K4me3). Using ChIP-seq, we find that androgen stimulation results in a 6-fold increase in the number of H3T11P-marked regions and induces WDR5 colocalization to one third of H3T11P-enriched promoters, thus establishing a genome-wide relationship between H3T11P and recruitment of WDR5. Accordingly, PKN1 knockdown or chemical inhibition severely blocks WDR5 chromatin association and H3K4me3 on AR target genes. Finally, WDR5 is critical in prostate cancer cell proliferation and is hyperexpressed in human prostate cancers. Together, these results identify WDR5 as a critical epigenomic integrator of histone phosphorylation and methylation and as a major driver of androgen-dependent prostate cancer cell proliferation.

A transcriptional regulatory role of the THAP11-HCF-1 complex in colon cancer cell function.

The recently identified Thanatos-associated protein (THAP) domain is an atypical zinc finger motif with sequence-specific DNA-binding activity. Emerging data suggest that THAP proteins may function in chromatin-dependent processes, including transcriptional regulation, but the roles of most THAP proteins in normal and aberrant cellular processes remain largely unknown. In this work, we identify THAP11 as a transcriptional regulator differentially expressed in human colon cancer. Immunohistochemical analysis of human colon cancers revealed increased THAP11 expression in both primary tumors and metastases. Knockdown of THAP11 in SW620 colon cancer cells resulted in a significant decrease in cell proliferation, and profiling of gene expression in these cells identified a novel gene set composed of 80 differentially expressed genes, 70% of which were derepressed by THAP11 knockdown. THAP11 was found to associate physically with the transcriptional coregulator HCF-1 (host cell factor 1) and recruit HCF-1 to target promoters. Importantly, THAP11-mediated gene regulation and its chromatin association require HCF-1, while HCF-1 recruitment at these genes requires THAP11. Collectively, these data provide the first characterization of THAP11-dependent gene expression in human colon cancer cells and suggest that the THAP11-HCF-1 complex may be an important transcriptional and cell growth regulator in human colon cancer.

Chromatin binding of SRp20 and ASF/SF2 and dissociation from mitotic chromosomes is modulated by histone H3 serine 10 phosphorylation.

Histone H3 serine 10 phosphorylation is a hallmark of mitotic chromosomes, but its full function remains to be elucidated. We report here that two SR protein splicing factors, SRp20 and ASF/SF2, associate with interphase chromatin, are released from hyperphosphorylated mitotic chromosomes, but reassociate with chromatin late in M-phase. Inhibition of Aurora B kinase diminished histone H3 serine 10 phosphorylation and increased SRp20 and ASF/SF2 retention on mitotic chromosomes. Unexpectedly, we also found that HP1 proteins interact with ASF/SF2 in mitotic cells. Strikingly, siRNA-mediated knockdown of ASF/SF2 caused retention of HP1 proteins on mitotic chromatin. Finally, ASF/SF2-depleted cells released from a mitotic block displayed delayed G0/G1 entry, suggesting a functional consequence of these interactions. These findings underscore the evolving role of histone H3 phosphorylation and demonstrate a direct, functional, and histone-modification-regulated association of SRp20 and ASF/SF2 with chromatin.

Hepatic expression of heme oxygenase-1 and antioxidant response element-mediated genes following administration of ethinyl estradiol to rats.

Heme oxygenase-1 (HO-1) is one of several enzymes induced by hepatotoxicants, and is thought to have an important protective role against cellular stress during liver inflammation and injury. The objective of the present study was to evaluate the role of HO-1 in estradiol-induced liver injury. A single dose of ethinyl estradiol (500 mg/kg, po) resulted in mild liver injury. Repeated administration of ethinyl estradiol (500 mg/kg/day for 4 days, po) resulted in no detectable liver injury or dysfunction. Using RT-PCR analysis, we demonstrate that HO-1 gene expression in whole liver tissue is elevated (>20-fold) after the single dose of ethinyl estradiol. The number and intensity of HO-1 immunoreactive macrophages were increased after the single dose of ethinyl estradiol. HO-1 expression was undetectable in hepatic parenchymal cells from rats receiving Methocel control or a single dose of ethinyl estradiol, however cytosolic HO-1 immunoreactivity in these cells after repeated dosing of ethinyl estradiol was pronounced. The increases in HO-1 mRNA and HO-1 immunoreactivity following administration of a single dose of ethinyl estradiol suggested that this enzyme might be responsible for the observed protection of the liver during repeated dosing. To investigate the effect of HO-1 expression on ethinyl estradiol-induced hepatotoxicity, rats were pretreated with hemin (50 micromol/kg, ip, a substrate and inducer of HO-1), with tin protoporphyrin IX (60 micromol/kg, ip, an HO-1 inhibitor), or with gadolinium chloride (10 mg/kg, iv, an inhibitor/toxin of Kupffer cells) 24 h before ethinyl estradiol treatment. Pretreatment with modulators of HO-1 expression and activity had generally minimal effects on ethinyl estradiol-induced liver injury. These data suggest that HO-1 plays a limited role in antioxidant defense against ethinyl estradiol-induced oxidative stress and hepatotoxicity, and suggests that other coordinately induced enzymes are responsible for protection observed with repeated administration of high doses of this compound.

Predictive toxicogenomics approaches reveal underlying molecular mechanisms of nongenotoxic carcinogenicity.

Toxicogenomics technology defines toxicity gene expression signatures for early predictions and hypotheses generation for mechanistic studies, which are important approaches for evaluating toxicity of drug candidate compounds. A large gene expression database built using cDNA microarrays and liver samples treated with over one hundred paradigm compounds was mined to determine gene expression signatures for nongenotoxic carcinogens (NGTCs). Data were obtained from male rats treated for 24 h. Training/testing sets of 24 NGTCs and 28 noncarcinogens were used to select genes. A semiexhaustive, nonredundant gene selection algorithm yielded six genes (nuclear transport factor 2, NUTF2; progesterone receptor membrane component 1, Pgrmc1; liver uridine diphosphate glucuronyltransferase, phenobarbital-inducible form, UDPGTr2; metallothionein 1A, MT1A; suppressor of lin-12 homolog, Sel1h; and methionine adenosyltransferase 1, alpha, Mat1a), which identified NGTCs with 88.5% prediction accuracy estimated by cross-validation. This six genes signature set also predicted NGTCs with 84% accuracy when samples were hybridized to commercially available CodeLink oligo-based microarrays. To unveil molecular mechanisms of nongenotoxic carcinogenesis, 125 differentially expressed genes (P<0.01) were selected by Student's t-test. These genes appear biologically relevant, of 71 well-annotated genes from these 125 genes, 62 were overrepresented in five biochemical pathway networks (most linked to cancer), and all of these networks were linked by one gene, c-myc. Gene expression profiling at early time points accurately predicts NGTC potential of compounds, and the same data can be mined effectively for other toxicity signatures. Predictive genes confirm prior work and suggest pathways critical for early stages of carcinogenesis.

Thanatos-associated protein 7 associates with template activating factor-Ibeta and inhibits histone acetylation to repress transcription.

The posttranslational modifications of histones on chromatin or a lack thereof is critical in transcriptional regulation. Emerging studies indicate a role for histone-binding proteins in transcriptional activation and repression. We have previously identified template-activating factor-Ibeta (TAF-Ibeta, also called PHAPII, SET, and I(2)(pp2A)) as a component of a cellular complex called inhibitor of acetyltransferases (INHAT) that masks histone acetylation in vitro and blocks histone acetyltransferase (HAT)-dependent transcription in living cells. TAF-Ibeta has also been shown to associate with transcription factors, including nuclear receptors, to regulate their activities. To identify novel interactors of TAF-Ibeta, we employed a yeast two-hybrid screen and identified a previously uncharacterized human protein called thanatos-associated protein-7 (THAP7), a member of a large family of THAP domain-containing putative DNA-binding proteins. In this study we demonstrate that THAP7 associates with TAF-Ibeta in vitro and map their association domains to a C-terminal predicted coiled-coil motif on THAP7 and the central region of TAF-Ibeta. Similarly, stably transfected THAP7 associates with endogenous TAF-Ibeta in intact cells. Like TAF-Ibeta, THAP7 associates with histone H3 and histone H4 and inhibits histone acetylation. The histone-interacting domain of THAP7 is sufficient for this activity in vitro. Promoter-targeted THAP7 can also recruit TAF-Ibeta and silencing mediator of retinoid and thyroid receptors/nuclear hormone receptor corepressor (NCoR) proteins to promoters, and knockdown of TAF-Ibeta by small interfering RNA relieves THAP7-mediated repression, indicating that, like nuclear hormone receptors, THAP7 may represent a novel class of transcription factor that uses TAF-Ibeta as a corepressor to maintain histones in a hypoacetylated, repressed state.

Glucocorticoid receptor-dependent gene regulatory networks.

While the molecular mechanisms of glucocorticoid regulation of transcription have been studied in detail, the global networks regulated by the glucocorticoid receptor (GR) remain unknown. To address this question, we performed an orthogonal analysis to identify direct targets of the GR. First, we analyzed the expression profile of mouse livers in the presence or absence of exogenous glucocorticoid, resulting in over 1,300 differentially expressed genes. We then executed genome-wide location analysis on chromatin from the same livers, identifying more than 300 promoters that are bound by the GR. Intersecting the two lists yielded 53 genes whose expression is functionally dependent upon the ligand-bound GR. Further network and sequence analysis of the functional targets enabled us to suggest interactions between the GR and other transcription factors at specific target genes. Together, our results further our understanding of the GR and its targets, and provide the basis for more targeted glucocorticoid therapies.

Drug-induced oxidative stress in rat liver from a toxicogenomics perspective.

Macrophage activators (MA), peroxisome proliferators (PP), and oxidative stressors/reactive metabolites (OS/RM) all produce oxidative stress and hepatotoxicity in rats. However, these three classes of hepatotoxicants give three distinct gene transcriptional profiles on cDNA microarrays, an indication that rat hepatocytes respond/adapt quite differently to these three classes of oxidative stressors. The differential gene responses largely reflect differential activation of transcription factors: MA activate Stat-3 and NFkB, PP activate PPARa, and OS/RM activate Nrf2. We have used gene signature profiles for each of these three classes of hepatotoxicants to categorize over 100 paradigm (and 50+ in-house proprietary) compounds as to their oxidative stress potential in rat liver. In addition to a role for microarrays in predictive toxicology, analyses of small subsets of these signature profiles, genes within a specific pathway, or even single genes often provide important insights into possible mechanisms involved in the toxicities of these compounds.

A gene expression signature for oxidant stress/reactive metabolites in rat liver.

Formation of free radicals and other reactive molecules is responsible for the adverse effects produced by a number of hepatotoxic compounds. cDNA microarray technology was used to compare transcriptional profiles elicited by training and testing sets of 15 oxidant stressors/reactive metabolite treatments to those produced by approximately 85 other paradigm compounds (mostly hepatotoxicants) to determine a shared signature profile for oxidant stress-associated hepatotoxicity. Initially, 100 genes were chosen that responded significantly different to oxidant stressors/reactive metabolites (OS/RM) compared to other samples in the database, then a 25-gene subset was selected by multivariate analysis. Many of the selected genes (e.g., aflatoxin aldehyde reductase, diaphorase, epoxide hydrolase, heme oxgenase and several glutathione transferases) are well-characterized oxidant stress/Nrf-2-responsive genes. Less than 10 other compounds co-cluster with our training and testing set compounds and these are known to generate OS/RMs as part of their mechanisms of toxicity. Using OS/RM signature gene sets, compounds previously associated with macrophage activation formed a distinct cluster separate from OS/RM and other compounds. A 69-gene set was chosen to maximally separate compounds in control, macrophage activator, peroxisome proliferator and OS/RM classes. The ease with which these 'oxidative stressor' classes can be separated indicates a role for microarray technology in early prediction and classification of hepatotoxicants. The ability to rapidly screen the oxidant stress potential of compounds may aid in avoidance of some idiosyncratic drug reactions as well as overtly toxic compounds.

Inverse gene expression patterns for macrophage activating hepatotoxicants and peroxisome proliferators in rat liver.

Macrophage activation contributes to adverse effects produced by a number of hepatotoxic compounds. Transcriptional profiles elicited by two macrophage activators, LPS and zymosan A, were compared to those produced by 100 paradigm compounds (mostly hepatotoxicants) using cDNA microarrays. Several hepatotoxicants previously reported to activate liver macrophages produced transcriptional responses similar to LPS and zymosan, and these were used to construct a gene signature profile for macrophage activators in the liver. Measurement of cytokine mRNAs in the same liver samples by RT-PCR independently confirmed that these compounds are associated with macrophage activation. In addition to expected effects on acute phase proteins and metabolic pathways that are regulated by LPS and inflammation, a strong induction was observed for many endoplasmic reticulum-associated stress/chaperone proteins. Additionally, many genes in our macrophage activator signature profile were well-characterized PPARalpha-induced genes which were repressed by macrophage activators. A shared gene signature profile for peroxisome proliferators was determined using a training set of clofibrate, WY 14643, diethylhexylphthalate, diisononylphthalate, perfluorodecanoic acid, perfluoroheptanoic acid, and perfluorooctanoic acid. The signature profile included macrophage activator-induced genes that were repressed by peroxisome proliferators. NSAIDs comprised an interesting pharmacological class in that some compounds, notably diflunisal, co-clustered with peroxisome proliferators whereas several others co-clustered with macrophage activators, possibly due to endotoxin exposure secondary to their adverse effects on the gastrointestinal system. While much of these data confirmed findings from the literature, the transcriptional patterns detected using this toxicogenomics approach showed relationships between genes and biological pathways requiring complex analysis to be discerned.