Michelle R. Campbell

Identification of novel NRF2-regulated genes by ChIP-Seq: influence on retinoid X receptor alpha

Cellular oxidative and electrophilic stress triggers a protective response in mammals regulated by NRF2 (nuclear factor (erythroid-derived) 2-like; NFE2L2) binding to deoxyribonucleic acid-regulatory sequences near stress-responsive genes. Studies using Nrf2-deficient mice suggest that hundreds of genes may be regulated by NRF2. To identify human NRF2-regulated genes, we conducted chromatin immunoprecipitation (ChIP)-sequencing experiments in lymphoid cells treated with the dietary isothiocyanate, sulforaphane (SFN) and carried out follow-up biological experiments on candidates. We found 242 high confidence, NRF2-bound genomic regions and 96% of these regions contained NRF2-regulatory sequence motifs. The majority of binding sites were near potential novel members of the NRF2 pathway. Validation of selected candidate genes using parallel ChIP techniques and in NRF2-silenced cell lines indicated that the expression of about two-thirds of the candidates are likely to be directly NRF2-dependent including retinoid X receptor alpha (RXRA). NRF2 regulation of RXRA has implications for response to retinoid treatments and adipogenesis. In mouse, 3T3-L1 cells’ SFN treatment affected Rxra expression early in adipogenesis, and knockdown of Nrf2-delayed Rxra expression, both leading to impaired adipogenesis.

Discovery and verification of functional single nucleotide polymorphisms in regulatory genomic regions: Current and developing technologies

The most common form of genetic variation, single nucleotide polymorphisms or SNPs, can affect the way an individual responds to the environment and modify disease risk. Although most of the millions of SNPs have little or no effect on gene regulation and protein activity, there are many circumstances where base changes can have deleterious effects. Non-synonymous SNPs that result in amino acid changes in proteins have been studied because of their obvious impact on protein activity. It is well known that SNPs within regulatory regions of the genome can result in disregulation of gene transcription. However, the impact of SNPs located in putative regulatory regions, or rSNPs, is harder to predict for two primary reasons. First, the mechanistic roles of non-coding genomic sequence remain poorly defined. Second, experimental validation of the functional consequences of rSNPs is often slow and laborious. In this review, we summarize traditional and novel methodologies for candidate rSNPs selection, in particular in silico techniques that aid in candidate rSNP selection. Additionally we will discuss molecular biological techniques that assess the impact of rSNPs on binding of regulatory machinery, as well as functional consequences on transcription. Standard techniques such as EMSA and luciferase reporter constructs are still widely used to assess effects of rSNPs on binding and gene transcription; however, these protocols are often bottlenecks in the discovery process. Therefore, we highlight novel and developing high-throughput protocols that promise to aid in shortening the process of rSNP validation. Given the large amount of genomic information generated from a multitude of re-sequencing and genome-wide SNP array efforts, future focus should be to develop validation techniques that will allow greater understanding of the impact these polymorphisms have on human health and disease.

A polymorphic p53 response element in KIT ligand influences cancer risk and has undergone natural selection.

The ability of p53 to regulate transcription is crucial for tumor suppression and implies that inherited polymorphisms in functional p53-binding sites could influence cancer. Here, we identify a polymorphic p53 responsive element and demonstrate its influence on cancer risk using genome-wide data sets of cancer susceptibility loci, genetic variation, p53 occupancy, and p53-binding sites. We uncover a single-nucleotide polymorphism (SNP) in a functional p53-binding site and establish its influence on the ability of p53 to bind to and regulate transcription of the KITLG gene. The SNP resides in KITLG and associates with one of the largest risks identified among cancer genome-wide association studies. We establish that the SNP has undergone positive selection throughout evolution, signifying a selective benefit, but go on to show that similar SNPs are rare in the genome due to negative selection, indicating that polymorphisms in p53-binding sites are primarily detrimental to humans.

Probing the functional impact of sequence variation on p53-DNA interactions using a novel microsphere assay for protein-DNA binding with human cell extracts.

The p53 tumor suppressor regulates its target genes through sequence-specific binding to DNA response elements (REs). Although numerous p53 REs are established, the thousands more identified by bioinformatics are not easily subjected to comparative functional evaluation. To examine the relationship between RE sequence variation—including polymorphisms—and p53 binding, we have developed a multiplex format microsphere assay of protein-DNA binding (MAPD) for p53 in nuclear extracts. Using MAPD we measured sequence-specific p53 binding of doxorubicin-activated or transiently expressed p53 to REs from established p53 target genes and p53 consensus REs. To assess the sensitivity and scalability of the assay, we tested 16 variants of the p21 target sequence and a 62-multiplex set of single nucleotide (nt) variants of the p53 consensus sequence and found many changes in p53 binding that are not captured by current computational binding models. A group of eight single nucleotide polymorphisms (SNPs) was examined and binding profiles closely matched transactivation capability tested in luciferase constructs. The in vitro binding characteristics of p53 in nuclear extracts recapitulated the cellular in vivo transactivation capabilities for eight well-established human REs measured by luciferase assay. Using a set of 26 bona fide REs, we observed distinct binding patterns characteristic of transiently expressed wild type and mutant p53s. This microsphere assay system utilizes biologically meaningful cell extracts in a multiplexed, quantitative, in vitro format that provides a powerful experimental tool for elucidating the functional impact of sequence polymorphism and protein variation on protein/DNA binding in transcriptional networks.

An African-specific polymorphism in the TP53 gene impairs p53 tumor suppressor function in a mouse model

A nonsynonymous single-nucleotide polymorphism at codon 47 in TP53 exists in African-descent populations (P47S, rs1800371; referred to here as S47). Here we report that, in human cell lines and a mouse model, the S47 variant exhibits a modest decrease in apoptosis in response to most genotoxic stresses compared with wild-type p53 but exhibits a significant defect in cell death induced by cisplatin. We show that, compared with wild-type p53, S47 has nearly indistinguishable transcriptional function but shows impaired ability to transactivate a subset of p53 target genes, including two involved in metabolism:Gls2(glutaminase 2) and Sco2 We also show that human and mouse cells expressing the S47 variant are markedly resistant to cell death by agents that induce ferroptosis (iron-mediated nonapoptotic cell death). We show that mice expressing S47 in homozygous or heterozygous form are susceptible to spontaneous cancers of diverse histological types. Our data suggest that the S47 variant may contribute to increased cancer risk in individuals of African descent, and our findings highlight the need to assess the contribution of this variant to cancer risk in these populations. These data also confirm the potential relevance of metabolism and ferroptosis to tumor suppression by p53.

Novel Hematopoietic Target Genes in the NRF2-Mediated Transcriptional Pathway

Nuclear factor- (erythroid-derived 2) like 2 (NFE2L2, NRF2) is a key transcriptional activator of the antioxidant response pathway and is closely related to erythroid transcription factor NFE2. Under oxidative stress, NRF2 heterodimerizes with small Maf proteins and binds cis-acting enhancer sequences found near oxidative stress response genes. Using the dietary isothiocyanate sulforaphane (SFN) to activate NRF2, chromatin immunoprecipitation sequencing (ChIP-seq) identified several hundred novel NRF2-mediated targets beyond its role in oxidative stress. Activated NRF2 bound the antioxidant response element (ARE) in promoters of several known and novel target genes involved in iron homeostasis and heme metabolism, including known targets FTL and FTH1, as well as novel binding in the globin locus control region. Five novel NRF2 target genes were chosen for followup: AMBP, ABCB6, FECH, HRG-1 (SLC48A1), and TBXAS1. SFN-induced gene expression in erythroid K562 and lymphoid cells were compared for each target gene. NRF2 silencing showed reduced expression in lymphoid, lung, and hepatic cells. Furthermore, stable knockdown of NRF2 negative regulator KEAP1 in K562 cells resulted in increased NQO1, AMBP, and TBXAS1 expression. NFE2 binding sites in K562 cells revealed similar binding profiles as lymphoid NRF2 sites in all potential NRF2 candidates supporting a role for NRF2 in heme metabolism and erythropoiesis.

Beyond antioxidant genes in the ancient Nrf2 regulatory network

Nrf2, a basic leucine zipper transcription factor encoded by the gene NFE2L2, is a master regulator of the transcriptional response to oxidative stress. Nrf2 is structurally and functionally conserved from insects to humans, and it heterodimerizes with the small MAF transcription factors to bind a consensus DNA sequence (the antioxidant response element, or ARE) and regulate gene expression. We have used genome-wide chromatin immunoprecipitation and gene expression data to identify direct Nrf2 target genes in Drosophila and humans. These data have allowed us to construct the deeply conserved ancient Nrf2 regulatory network-target genes that are conserved from Drosophila to human. The ancient network consists of canonical antioxidant genes, as well as genes related to proteasomal pathways and metabolism and a number of less expected genes. We have also used enhancer reporter assays and electrophoretic mobility-shift assays to confirm Nrf2-mediated regulation of ARE activity at a number of these novel target genes. Interestingly, the ancient network also highlights a prominent negative feedback loop; this, combined with the finding that Nrf2-mediated regulatory output is tightly linked to the quality of the ARE it is targeting, suggests that precise regulation of nuclear Nrf2 concentration is necessary to achieve proper quantitative regulation of distinct gene sets. Together, these findings highlight the importance of balance in the Nrf2-ARE pathway and indicate that Nrf2-mediated regulation of xenobiotic metabolism, glucose metabolism, and proteostasis has been central to this pathway since its inception.

Human single-nucleotide polymorphisms alter p53 sequence-specific binding at gene regulatory elements

p53 coordinates the expression of an intricate network of genes in response to stress signals. Sequence-specific DNA binding is essential for p53-mediated tumor suppression. We evaluated the impact of single-nucleotide polymorphisms (SNPs) in p53 response elements (p53RE) on DNA binding and gene expression in response to DNA damage. Using a bioinformatics approach based on incorporating p53 binding strength into a position weight matrix, we selected 32 SNPs in putative and validated p53REs. The microsphere assay for protein–DNA binding (MAPD) and allele-specific expression analysis was employed to assess the impact of SNPs on p53-DNA binding and gene expression, respectively. Comparing activated p53 binding in nuclear extracts from doxorubicin- or ionizing radiation (IR)-treated human cells, we observed little difference in binding profiles. Significant p53 binding was observed for most polymorphic REs and several displayed binding comparable to the p21 RE. SNP alleles predicted to lower p53 binding indeed reduced binding in 25 of the 32 sequences. Chromatin immunoprecipitation-sequencing in lymphoblastoid cells confirmed p53 binding to seven polymorphic p53 REs in response to doxorubicin. In addition, five polymorphisms were associated with altered gene expression following doxorubicin treatment. Our findings demonstrate an effective strategy to identify and evaluate SNPs that may alter p53-mediated stress responses.

Interactions of Chromatin Context, Binding Site Sequence Content, and Sequence Evolution in Stress-Induced p53 Occupancy and Transactivation

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10−7) and correlated with stronger p53RE sequences (p<10−110) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

Distinct epigenetic effects of tobacco smoking in whole blood and among leukocyte subtypes

Tobacco smoke exposure dramatically alters DNA methylation in blood cells and may mediate smoking-associated complex diseases through effects on immune cell function. However, knowledge of smoking effects in specific leukocyte subtypes is limited. To better characterize smoking–associated methylation changes in whole blood and leukocyte subtypes, we used Illumina 450K arrays and Reduced Representation Bisulfite Sequencing (RRBS) to assess genome-wide DNA methylation. Differential methylation analysis in whole blood DNA from 172 smokers and 81 nonsmokers revealed 738 CpGs, including 616 previously unreported CpGs, genome-wide significantly associated with current smoking (p <1.2x10-7, Bonferroni correction). Several CpGs (MTSS1, NKX6-2, BTG2) were associated with smoking duration among heavy smokers (>22 cigarettes/day, n = 86) which might relate to long-term heavy-smoking pathology. In purified leukocyte subtypes from an independent group of 20 smokers and 14 nonsmokers we further examined methylation and gene expression for selected genes among CD14+ monocytes, CD15+ granulocytes, CD19+ B cells, and CD2+ T cells. In 10 smokers and 10 nonsmokers we used RRBS to fine map differential methylation in CD4+ T cells, CD8+ T cells, CD14+, CD15+, CD19+, and CD56+ natural killer cells. Distinct cell-type differences in smoking-associated methylation and gene expression were identified. AHRR (cg05575921), ALPPL2 (cg21566642), GFI1 (cg09935388), IER3 (cg06126421) and F2RL3 (cg03636183) showed a distinct pattern of significant smoking-associated methylation differences across cell types: granulocytes> monocytes>> B cells. In contrast GPR15 (cg19859270) was highly significant in T and B cells and ITGAL (cg09099830) significant only in T cells. Numerous other CpGs displayed distinctive cell-type responses to tobacco smoke exposure that were not apparent in whole blood DNA. Assessing the overlap between these CpG sites and differential methylated regions (DMRs) with RRBS in 6 cell types, we confirmed cell-type specificity in the context of DMRs. We identified new CpGs associated with current smoking, pack-years, duration, and revealed unique profiles of smoking-associated DNA methylation and gene expression among immune cell types, providing potential clues to hematopoietic lineage-specific effects in disease etiology.