Decoding a signature-based model of transcription cofactor recruitment dictated by cardinal cis-regulatory elements in proximal promoter regions

Q-rich activation domains: flexible 'rulers' for transcription start site selection?

Recent findings broadened the function of RNA polymerase II (Pol II) proximal promoter motifs from quantitative regulators of transcription to important determinants of transcription start site (TSS) position. These motifs are recognized by transcription factors (TFs) that we propose to term 'ruler' TFs (rTFs), such as NRF1, NF-Y, YY1, ZNF143, BANP, and members of the SP, ETS, and CRE families, sharing as a common feature a glutamine-rich (Q-rich) effector domain also enriched in valine, isoleucine, and threonine (QVIT-rich). We propose that rTFs guide TSS location by constraining the position of the pre-initiation complex (PIC) during its promoter recognition phase through a specialized, and still enigmatic, class of activation domains.

ZNF143 binds DNA and stimulates transcription initiation to activate and repress direct target genes

Transcription factors bind to sequence motifs and act as activators or repressors. Transcription factors interface with a constellation of accessory cofactors to regulate distinct mechanistic steps to regulate transcription. We rapidly degraded the essential and pervasively expressed transcription factor ZNF143 to determine its function in the transcription cycle. ZNF143 facilitates RNA polymerase initiation and activates gene expression. ZNF143 binds the promoter of nearly all its activated target genes. ZNF143 also binds near the site of genic transcription initiation to directly repress a subset of genes. Although ZNF143 stimulates initiation at ZNF143-repressed genes (i.e. those that increase transcription upon ZNF143 depletion), the molecular context of binding leads to cis repression. ZNF143 competes with other more efficient activators for promoter access, physically occludes transcription initiation sites and promoter-proximal sequence elements, and acts as a molecular roadblock to RNA polymerases during early elongation. The term context specific is often invoked to describe transcription factors that have both activation and repression functions. We define the context and molecular mechanisms of ZNF143-mediated cis activation and repression.

ZNF143 binds DNA and stimulates transcripstion initiation to activate and repress direct target genes

Transcription factors bind to sequence motifs and act as activators or repressors. Transcription factors interface with a constellation of accessory cofactors to regulate distinct mechanistic steps to regulate transcription. We rapidly degraded the essential and ubiquitously expressed transcription factor ZNF143 to determine its function in the transcription cycle. ZNF143 facilitates RNA Polymerase initiation and activates gene expression. ZNF143 binds the promoter of nearly all its activated target genes. ZNF143 also binds near the site of genic transcription initiation to directly repress a subset of genes. Although ZNF143 stimulates initiation at ZNF143-repressed genes (i.e. those that increase expression upon ZNF143 depletion), the molecular context of binding leads to cis repression. ZNF143 competes with other more efficient activators for promoter access, physically occludes transcription initiation sites and promoter-proximal sequence elements, and acts as a molecular roadblock to RNA Polymerases during early elongation. The term context specific is often invoked to describe transcription factors that have both activation and repression functions. We define the context and molecular mechanisms of ZNF143-mediated cis activation and repression.

Predicting the impact of sequence motifs on gene regulation using single-cell data

The binding of transcription factors at proximal promoters and distal enhancers is central to gene regulation. Identifying regulatory motifs and quantifying their impact on expression remains challenging. Using a convolutional neural network trained on single-cell data, we infer putative regulatory motifs and cell type-specific importance. Our model, scover, explains 29% of the variance in gene expression in multiple mouse tissues. Applying scover to distal enhancers identified using scATAC-seq from the developing human brain, we identify cell type-specific motif activities in distal enhancers. Scover can identify regulatory motifs and their importance from single-cell data where all parameters and outputs are easily interpretable.

Targeting the LSD1/KDM1 Family of Lysine Demethylases in Cancer and Other Human Diseases

Lysine-specific demethylase 1 (LSD1) was the first histone demethylase discovered and the founding member of the flavin-dependent lysine demethylase family (KDM1). The human KDM1 family includes KDM1A and KDM1B, which primarily catalyze demethylation of histone H3K4me1/2. The KDM1 family is involved in epigenetic gene regulation and plays important roles in various biological and disease pathogenesis processes, including cell differentiation, embryonic development, hormone signaling, and carcinogenesis. Malfunction of many epigenetic regulators results in complex human diseases, including cancers. Regulators such as KDM1 have become potential therapeutic targets because of the reversibility of epigenetic control of genome function. Indeed, several classes of KDM1-selective small molecule inhibitors have been developed, some of which are currently in clinical trials to treat various cancers. In this chapter, we review the discovery, biochemical, and molecular mechanisms, atomic structure, genetics, biology, and pathology of the KDM1 family of lysine demethylases. Focusing on cancer, we also provide a comprehensive summary of recently developed KDM1 inhibitors and related preclinical and clinical studies to provide a better understanding of the mechanisms of action and applications of these KDM1-specific inhibitors in therapeutic treatment.

Differential dependencies of human RNA polymerase II promoters on TBP, TAF1, TFIIB and XPB

The effects of rapid acute depletion of components of RNA polymerase II (Pol II) general transcription factors (GTFs) that are thought to be critical for formation of preinitiation complexes (PICs) and initiation in vitro were quantified in HAP1 cells using precision nuclear run-on sequencing (PRO-Seq). The average dependencies for each factor across >70 000 promoters varied widely even though levels of depletions were similar. Some of the effects could be attributed to the presence or absence of core promoter elements such as the upstream TBP-specificity motif or downstream G-rich sequences, but some dependencies anti-correlated with such sequences. While depletion of TBP had a large effect on most Pol III promoters only a small fraction of Pol II promoters were similarly affected. TFIIB depletion had the largest general effect on Pol II and also correlated with apparent termination defects downstream of genes. Our results demonstrate that promoter activity is combinatorially influenced by recruitment of TFIID and sequence-specific transcription factors. They also suggest that interaction of the preinitiation complex (PIC) with nucleosomes can affect activity and that recruitment of TFIID containing TBP only plays a positive role at a subset of promoters.

N-Acetylcysteine Reduces miR-146a and NF-κB p65 Inflammatory Signaling Following Cadmium Hepatotoxicity in Rats

We performed a thorough screening and analysis of the impact of cadmium chloride (CdCl2) and N-acetylcysteine (NAC) on the miR146a/NF-κB p65 inflammatory pathway and mitochondrial biogenesis dysfunction in male albino rats. A total of 24 male albino rats were divided into three groups: a control group, a CdCl2-treated group (3 mg/kg, orally), and a CdCl2 + NAC-treated group (200 mg/kg of NAC, 1 h after CdCl2 treatment), for 60 consecutive days. Real-time quantitative PCR was used to analyze the expression of miR146a, Irak1, Traf6, Nrf1, Nfe2l2, Pparg, Prkaa, Stat3, Tfam, Tnfa, and Il1b, whereas tumor necrosis factor-α, interleukin-1β, and cyclooxygenase-2 protein levels were assessed using ELISA, and NF-κB p65 was detected using western blotting. A significant restoration of homeostatic inflammatory processes as well as mitochondrial biogenesis was observed after NAC and CdCl2 treatment. Decreased miR146a and NF-κB p65 were also found after treatment with NAC and CdCl2 compared with CdCl2 treatment alone. Collectively, our findings demonstrate that CdCl2 caused mtDNA release because of Tfam loss, leading to NF-κB p65 activation. Co-treatment with NAC could alleviate Cd-induced genotoxicity in liver tissue. We concluded that adding NAC to CdCl2 resulted in a decreased signaling of the NF-κB p65 signaling pathway.

Chemical Inhibition of Apurinic-Apyrimidinic Endonuclease 1 Redox and DNA Repair Functions Affects the Inflammatory Response via Different but Overlapping Mechanisms

The presence of oxidized DNA lesions, such as 7,8-dihydro-8-oxoguanine (8-oxoG) and apurinic/apyrimidinic sites (AP sites), has been described as epigenetic signals that are involved in gene expression control. In mammals, Apurinic-apyrimidinic endonuclease 1/Redox factor-1 (APE1/Ref-1) is the main AP endonuclease of the base excision repair (BER) pathway and is involved in active demethylation processes. In addition, APE1/Ref-1, through its redox function, regulates several transcriptional factors. However, the transcriptional control targets of each APE1 function are not completely known. In this study, a transcriptomic approach was used to investigate the effects of chemical inhibition of APE1/Ref-1 redox or DNA repair functions by E3330 or methoxyamine (MX) in an inflammatory cellular model. Under lipopolysaccharide (LPS) stimulation, both E3330 and MX reduced the expression of some cytokines and chemokines. Interestingly, E3330 treatment reduced cell viability after 48 h of the treatment. Genes related to inflammatory response and mitochondrial processes were downregulated in both treatments. In the E3330 treatment, RNA processing and ribosome biogenesis genes were downregulated, while they were upregulated in the MX treatment. Furthermore, in the E3330 treatment, the cellular stress response was the main upregulated process, while the cellular macromolecule metabolic process was observed in MX-upregulated genes. Nuclear respiratory factor 1 (NRF1) was predicted to be a master regulator of the downregulated genes in both treatments, while the ETS transcription factor ELK1 (ELK1) was predicted to be a master regulator only for E3330 treatment. Decreased expression of ELK1 and its target genes and a reduced 28S/18S ratio were observed, suggesting impaired rRNA processing. In addition, both redox and repair functions can affect the expression of NRF1 and GABPA target genes. The master regulators predicted for upregulated genes were YY1 and FLI1 for the E3330 and MX treatments, respectively. In summary, the chemical inhibition of APE1/Ref-1 affects gene expression regulated mainly by transcriptional factors of the ETS family, showing partial overlap of APE1 redox and DNA repair functions, suggesting that these activities are not entirely independent. This work provides a new perspective on the interaction between APE1 redox and DNA repair activity in inflammatory response modulation and transcription.

Analysis of Clonal Composition in Human iPSC and ESC and Derived 2D and 3D Differentiated Cultures

Human induced pluripotent and embryonic stem cell cultures (hiPSC/hESC) are phenotypically heterogeneous and prone to clonal deviations during subculturing and differentiation. Clonal deviations often emerge unnoticed, but they can change the biology of the cell culture with a negative impact on experimental reproducibility. Here, we describe a computational workflow to profile the bulk clonal composition in a hiPSC/hESC culture that can also be used to infer clonal deviations. This workflow processes data obtained with two versions of the same method. The two versions-epigenetic and transcriptomic-rely on a mechanism of stochastic H3K4me3 deposition during hiPSC/hESC derivation. This mechanism generates a signature of ten or more H3K4me3-enriched clustered protocadherin (PCDH) promoters distinct in every single cell. The aggregate of single-cell signatures provides an identificatory feature in every hiPSC/hESC line. This feature is stably transmitted to the cell progeny of the culture even after differentiation unless there is a clonal deviation event that changes the internal balance of single-cell signatures. H3K4me3 signatures can be profiled by chromatin immunoprecipitation and next-generation sequencing (ChIP-seq). Alternatively, an equivalent PCDH-expression version can be profiled by RNA-seq in PCDH-expressing hiPSC/hESC-derived cells (such as neurons, astrocytes, and cardiomyocytes; and, in long-term cultures, such as cerebral organoids). Notably, our workflow can also distinguish genetically identical hiPSC/hESC lines derived from the same patient or generated in the same editing process. Together, we propose a method to improve data sharing and reproducibility in the hiPSC and hESC fields.

INTACT vs. FANS for Cell-Type-Specific Nuclei Sorting: A Comprehensive Qualitative and Quantitative Comparison

Increasing numbers of studies seek to characterize the different cellular sub-populations present in mammalian tissues. The techniques “Isolation of Nuclei Tagged in Specific Cell Types” (INTACT) or “Fluorescence-Activated Nuclei Sorting” (FANS) are frequently used for isolating nuclei of specific cellular subtypes. These nuclei are then used for molecular characterization of the cellular sub-populations. Despite the increasing popularity of both techniques, little is known about their isolation efficiency, advantages, and disadvantages or downstream molecular effects. In our study, we compared the physical and molecular attributes of sfGFP+ nuclei isolated by the two methods—INTACT and FANS—from the neocortices of Arc-CreERT2 × CAG-Sun1/sfGFP animals. We identified differences in efficiency of sfGFP+ nuclei isolation, nuclear size as well as transcriptional (RNA-seq) and chromatin accessibility (ATAC-seq) states. Therefore, our study presents a comprehensive comparison between the two widely used nuclei sorting techniques, identifying the advantages and disadvantages for both INTACT and FANS. Our conclusions are summarized in a table to guide researchers in selecting the most suitable methodology for their individual experimental design.

Super-Enhancer Redistribution as a Mechanism of Broad Gene Dysregulation in Repeatedly Drug-Treated Cancer Cells

Cisplatin is an antineoplastic drug administered at suboptimal and intermittent doses to avoid life-threatening effects. Although this regimen shortly improves symptoms in the short term, it also leads to more malignant disease in the long term. We describe a multilayered analysis ranging from chromatin to translation-integrating chromatin immunoprecipitation sequencing (ChIP-seq), global run-on sequencing (GRO-seq), RNA sequencing (RNA-seq), and ribosome profiling-to understand how cisplatin confers (pre)malignant features by using a well-established ovarian cancer model of cisplatin exposure. This approach allows us to segregate the human transcriptome into gene modules representing distinct regulatory principles and to characterize that the most cisplatin-disrupted modules are associated with underlying events of super-enhancer plasticity. These events arise when cancer cells initiate without ultimately ending the program of drug-stimulated death. Using a PageRank-based algorithm, we predict super-enhancer regulator ISL1 as a driver of this plasticity and validate this prediction by using CRISPR/dCas9-KRAB inhibition (CRISPRi) and CRISPR/dCas9-VP64 activation (CRISPRa) tools. Together, we propose that cisplatin reprograms cancer cells when inducing them to undergo near-to-death experiences.

Analysis of Brugada syndrome loci reveals that fine-mapping clustered GWAS hits enhances the annotation of disease-relevant variants

Genome-wide association studies (GWASs) are instrumental in identifying loci harboring common single-nucleotide variants (SNVs) that affect human traits and diseases. GWAS hits emerge in clusters, but the focus is often on the most significant hit in each trait- or disease-associated locus. The remaining hits represent SNVs in linkage disequilibrium (LD) and are considered redundant and thus frequently marginally reported or exploited. Here, we interrogate the value of integrating the full set of GWAS hits in a locus repeatedly associated with cardiac conduction traits and arrhythmia, SCN5A-SCN10A. Our analysis reveals 5 common 7-SNV haplotypes (Hap1-5) with 2 combinations associated with life-threatening arrhythmia-Brugada syndrome (the risk Hap1/1 and protective Hap2/3 genotypes). Hap1 and Hap2 share 3 SNVs; thus, this analysis suggests that assuming redundancy among clustered GWAS hits can lead to confounding disease-risk associations and supports the need to deconstruct GWAS data in the context of haplotype composition.

ESR1 ChIP-Seq Identifies Distinct Ligand-Free ESR1 Genomic Binding Sites in Human Hepatocytes and Liver Tissue

The estrogen receptor alpha (ESR1) is an important gene transcriptional regulator, known to mediate the effects of estrogen. Canonically, ESR1 is activated by its ligand estrogen. However, the role of unliganded ESR1 in transcriptional regulation has been gaining attention. We have recently shown that ligand-free ESR1 is a key regulator of several cytochrome P450 (CYP) genes in the liver, however ligand-free ESR1 has not been characterized genome-wide in the human liver. To address this, ESR1 ChIP-Seq was conducted in human liver samples and in hepatocytes with or without 17beta-estradiol (E2) treatment. We identified both ligand-dependent and ligand-independent binding sites throughout the genome. These two ESR1 binding categories showed different genomic localization, pathway enrichment, and cofactor colocalization, indicating different ESR1 regulatory function depending on ligand availability. By analyzing existing ESR1 data from additional human cell lines, we uncovered a potential ligand-independent ESR1 activity, namely its co-enrichment with the zinc finger protein 143 (ZNF143). Furthermore, we identified ESR1 binding sites near many gene loci related to drug therapy, including the CYPs. Overall, this study shows distinct ligand-free and ligand-bound ESR1 chromatin binding profiles in the liver and suggests the potential broad influence of ESR1 in drug metabolism and drug therapy.

Deep neural networks identify sequence context features predictive of transcription factor binding

Transcription factors (TFs) bind DNA by recognizing specific sequence motifs, typically of length 6-12bp. A motif can occur many thousands of times in the human genome, but only a subset of those sites are actually bound. Here we present a machine learning framework leveraging existing convolutional neural network architectures and model interpretation techniques to identify and interpret sequence context features most important for predicting whether a particular motif instance will be bound. We apply our framework to predict binding at motifs for 38 TFs in a lymphoblastoid cell line, score the importance of context sequences at base-pair resolution, and characterize context features most predictive of binding. We find that the choice of training data heavily influences classification accuracy and the relative importance of features such as open chromatin. Overall, our framework enables novel insights into features predictive of TF binding and is likely to inform future deep learning applications to interpret non-coding genetic variants.

Altered Enhancer and Promoter Usage Leads to Differential Gene Expression in the Normal and Failed Human Heart

BACKGROUND

The failing heart is characterized by changes in gene expression. However, the regulatory regions of the genome that drive these gene expression changes have not been well defined in human hearts.

METHODS

To define genome-wide enhancer and promoter use in heart failure, cap analysis of gene expression sequencing was applied to 3 healthy and 4 failed human hearts to identify promoter and enhancer regions used in left ventricles. Healthy hearts were derived from donors unused for transplantation and failed hearts were obtained as discarded tissue after transplantation.

RESULTS

Cap analysis of gene expression sequencing identified a combined potential for ≈23 000 promoters and ≈5000 enhancers active in human left ventricles. Of these, 17 000 promoters and 1800 enhancers had additional support for their regulatory function. Comparing promoter usage between healthy and failed hearts highlighted promoter shifts which altered aminoterminal protein sequences. Enhancer usage between healthy and failed hearts identified a majority of differentially used heart failure enhancers were intronic and primarily localized within the first intron, revealing this position as a common feature associated with tissue-specific gene expression changes in the heart.

CONCLUSIONS

This data set defines the dynamic genomic regulatory landscape underlying heart failure and serves as an important resource for understanding genetic contributions to cardiac dysfunction. Additionally, regulatory changes contributing to heart failure are attractive therapeutic targets for controlling ventricular remodeling and clinical progression.

Mitochondrial Dysfunction in Kidney Disease and Uremic Sarcopenia

Recently, there has been an increased focus on the influences of mitochondrial dysfunction on various pathologies. Mitochondria are major intracellular organelles with a variety of critical roles, such as adenosine triphosphate production, metabolic modulation, generation of reactive oxygen species, maintenance of intracellular calcium homeostasis, and the regulation of apoptosis. Moreover, mitochondria are attracting attention as a therapeutic target in several diseases. Additionally, a lot of existing agents have been found to have pharmacological effects on mitochondria. This review provides an overview of the mitochondrial change in the kidney and skeletal muscle, which is often complicated with sarcopenia and chronic kidney disease (CKD). Furthermore, the pharmacological effects of therapeutics for CKD on mitochondria are explored.

Chromatin establishes an immature version of neuronal protocadherin selection during the naive-to-primed conversion of pluripotent stem cells

In the mammalian genome, the clustered protocadherin (cPcdh) locus is a paradigm of stochastic gene expression with the potential to generate a unique cPcdh combination in every neuron. Here, we report a chromatin-based mechanism emerging during the transition from the naive to the primed states of cell pluripotency that reduces by orders of magnitude the combinatorial potential in the human cPcdh locus. This mechanism selectively increases the frequency of stochastic selection of a small subset of cPcdh genes after neuronal differentiation in monolayers, months-old organoids, and engrafted cells in the rat spinal cord. Signs of these frequent selections can be observed in the brain throughout fetal development and disappear after birth, unless there is a condition of delayed maturation such as Down Syndrome. We therefore propose that a pattern of limited cPcdh diversity is maintained while human neurons still retain fetal-like levels of maturation.

Short and long-term cultures of human stem cell-derived neurons reveal that a pattern of restricted selection of clustered protocadherin isoforms, pre-established in pluripotent cells, distinguishes immature from mature neurons.

LSD1 suppresses invasion, migration and metastasis of luminal breast cancer cells via activation of GATA3 and repression of TRIM37 expression

LSD1 (KDM1A) is a histone demethylase that plays both oncogenic and tumor suppressor roles in breast cancer. However, the exact contexts under which it plays these opposite functions remain largely elusive. By characterizing its role in luminal breast epithelial cells, here we show that inhibition of LSD1 by both genetic and pharmacological approaches increases their invasion and migration, whereas its inhibition by genetic approach, but not by pharmacological approach, impairs their proliferation/survival. Induced loss of LSD1 in luminal cells in a mouse model of luminal breast cancer, MMTV-PyMT, leads to a profound increase in lung metastasis. Mechanistically, LSD1 interacts with GATA3, a key luminal-specific transcription factor (TF), and their common target genes are highly related to breast cancer. LSD1 positively regulates GATA3 expression. It also represses expression of TRIM37, a breast epithelial oncogene encoding a histone H2A ubiquitin ligase, and ELF5, a key TF gene for luminal progenitors and alveolar luminal cells. LSD1-loss also leads to reduced expression of several cell-cell adhesion genes (e.g., CDH1, VCL, CTNNA1), possibly via TRIM37-upregulation and subsequently TRIM37-mediated repression. Collectively, our data suggest LSD1 largely plays a tumor suppressor role in luminal breast cancer and the oncogenic program associated with LSD1-inhibition may be suppressed via TRIM37-inhibition.

Coordinated Transcriptional Regulation of Cytochrome P450 3As by Nuclear Transcription Factor Y and Specificity Protein 1

The cytochrome P450 3A subfamily plays vital roles in the metabolism of endogenous chemicals and xenobiotics. Understanding the basal expression of CYP3A in humans and pigs is crucial for drug evaluation. In this study, we demonstrated that the basal transcriptional regulation of CYP3A genes in hepatocytes is evolutionarily conserved between humans and pigs. The basal expression of CYP3A genes is transactivated by two cis-acting elements, the CCAAT and GC boxes, located a constant distance apart in the proximal promoter region of six CYP3A genes. Mutation analysis of these two cis-acting elements suggested that they play important roles in mediating basal expression, but to different extents because of the nucleotide variations in the elements. Two transcription factors, nuclear transcription factor Y (NF-Y) and specificity protein 1 (Sp1), directly bind to these cis-acting elements in CYP3A proximal promoters in HepG2 cells and porcine hepatocytes. Furthermore, changing the distance between the NF-Y and Sp1 binding sites resulted in decreases in the promoter activity of CYP3A genes. Conclusively, our results show that human and porcine CYP3A genes are regulated by NF-Y and Sp1 in a coordinated manner, and that the distance between these two cis-acting elements is crucial for constitutive CYP3A expression.

Bioinformatic and functional analysis of promoter region of human SLC25A13 gene

The human SLC25A13 gene encodes the liver type aspartate/glutamate carrier isoform 2 (AGC2, commonly named as citrin), which plays a key role in the main NADH-shuttle of human hepatocyte. Biallelic SLC25A13 mutations result in Citrin deficiency (CD). In order to identify the important regulatory region of SLC25A13 gene and elucidate the way how potential promoter mutations affect the citrin expression, we performed promoter deletion analysis and established the reporter constructs of luciferase gene-carrying SLC25A13 promoter containing several mutations located in putative transcription factor-binding sites. The luciferase activities of all promoter constructs were measured using a Dual-Luciferase Reporter Assay System. Bioinformatic analysis showed that the promoter of SLC25A13 gene lacks TATA box and obviously typical initiator element, but contains a CCAAT box and two GC box. Promoter deletion analysis confirmed the region from -221 to -1 upstream ATG was essential for SLC25A13 to maintain the promoter activity. We utilized dual-luciferase reporter system as function analytical model to tentatively assess the effect of artificially constructed promoter mutations on citrin expression, and our analysis revealed that mutated putative CCAAT box and GC box could significantly affect the citrin expression. Our study confirmed the important SLC25A13 promoter regions that influenced citrin expression in HL7702 cells, and constructed a function analytical model. This work may be useful to further identify the pathogenic mutations leading to CD in the promoter region.

Identification of factors associated with duplicate rate in ChIP-seq data

Chromatin immunoprecipitation and sequencing (ChIP-seq) has been widely used to map DNA-binding proteins, histone proteins and their modifications. ChIP-seq data contains redundant reads termed duplicates, referring to those mapping to the same genomic location and strand. There are two main sources of duplicates: polymerase chain reaction (PCR) duplicates and natural duplicates. Unlike natural duplicates that represent true signals from sequencing of independent DNA templates, PCR duplicates are artifacts originating from sequencing of identical copies amplified from the same DNA template. In analysis, duplicates are removed from peak calling and signal quantification. Nevertheless, a significant portion of the duplicates is believed to represent true signals. Obviously, removing all duplicates will underestimate the signal level in peaks and impact the identification of signal changes across samples. Therefore, an in-depth evaluation of the impact from duplicate removal is needed. Using eight public ChIP-seq datasets from three narrow-peak and two broad-peak marks, we tried to understand the distribution of duplicates in the genome, the extent by which duplicate removal impacts peak calling and signal estimation, and the factors associated with duplicate level in peaks. The three PCR-free histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq data had about 40% duplicates and 97% of them were within peaks. For the other datasets generated with PCR amplification of ChIP DNA, as expected, the narrow-peak marks have a much higher proportion of duplicates than the broad-peak marks. We found that duplicates are enriched in peaks and largely represent true signals, more conspicuous in those with high confidence. Furthermore, duplicate level in peaks is strongly correlated with the target enrichment level estimated using nonredundant reads, which provides the basis to properly allocate duplicates between noise and signal. Our analysis supports the feasibility of retaining the portion of signal duplicates into downstream analysis, thus alleviating the limitation of complete deduplication.

Nuclear Respiratory Factor 1 Acting as an Oncoprotein Drives Estrogen-Induced Breast Carcinogenesis

We have previously shown nuclear respiratory factor 1 (NRF1)-mediated transcriptional programming of mitobiogenesis contributes to estrogen-induced breast cancer through modulating cell cycle progression. In this study, we report a new role of NRF1 that goes beyond that of programming mitobiogenesis. Specifically, we report a novel oncogenic function of NRF1 supporting its causative role in breast cancer development and progression. The gain of NRF1 and/or treatment with 17β-estradiol (E2) produced heterogeneous breast cancer stem cell (BCSC)-like subsets composed of more than 10 distinct cell sub-populations. Flow sorting combined with confocal imaging of markers for pluripotency, epithelial mesenchymal transition (EMT), and BCSCs phenotypically confirmed that the BCSC-like subset arise from cell re-programming. Thus, we determined the molecular actions of NRF1 on its target gene CXCR4 because of its known role in the acquisition of the BCSC-like subset through EMT. CXCR4 was activated by NRF1 in a redox-dependent manner during malignant transformation. An NRF1-induced BCSC-like subset was able to form xenograft tumors in vivo, while inhibiting transcription of CXCR4 prevented xenograft tumor growth. Consistent with our observation of NRF1-driven breast tumorigenesis in the experimental model, higher protein levels of NRF1 were also found in human breast cancer tissue specimens. This highly novel role of NRF1 in the stochastic acquisition of BCSC-like subsets and their progression to a malignant phenotype may open an entirely new research direction targeting NRF1 signaling in invasive breast cancer. Our discovery of targeting transcriptional activation of CXCR4 to inhibit NRF1-induced oncogenic transformation provides a mechanistic explanation for estrogen-dependent breast carcinogenesis and opens new avenues in strategic therapeutics to fight breast cancer.

LSD1-ERRα complex requires NRF1 to positively regulate transcription and cell invasion

Lysine-specific demethylase 1 (LSD1) exerts dual effects on histone H3, promoting transcriptional repression via Lys4 (H3K4) demethylation or transcriptional activation through Lys9 (H3K9) demethylation. These activities are often exerted at transcriptional start sites (TSSs) and depend on the type of enhancer-bound transcription factor (TFs) with which LSD1 interacts. In particular, the Estrogen-Receptor Related α (ERRα) TF interacts with LSD1 and switches its activities toward H3K9 demethylation, resulting in transcriptional activation of a set of common target genes. However, how are the LSD1-TF and, in particular LSD1-ERRα, complexes determined to act at TSSs is not understood. Here we show that promoter-bound nuclear respiratory factor 1 (NRF1), but not ERRα, is essential to LSD1 recruitment at the TSSs of positive LSD1-ERRα targets. In contrast to ERRα, NRF1 does not impact on the nature of LSD1 enzymatic activity. We propose a three factor model, in which the LSD1 histone modifier requires a TSS tethering factor (NRF1) as well as an activity inducer (ERRα) to transcriptionally activate common targets. The relevance of this common network is illustrated by functional data, showing that all three factors are required for cell invasion in an MMP1 (Matrix MetalloProtease 1)-dependent manner, the expression of which is regulated by NRF1/LSD1/ERRα-mediated H3K9me2 demethylation.

Increased expression of mitochondrial transcription factor A and nuclear respiratory factor-1 predicts a poor clinical outcome of breast cancer

Nuclear respiratory factor-1 (Nrf1) and mitochondrial transcription factor A (TFAM) are involved in the regulation of a variety of mitochondrial functional genes, which are associated with decreased sensitivity of tumor cells to chemotherapy. However, the expression status of Nrf1 and TFAM, as well as their clinical significance in breast cancer, is unknown. In the present study, tumor tissues and corresponding adjacent normal tissues were collected from 336 patients with breast cancer, and Nrf1 and TFAM expression was analyzed by immunohistochemistry using a tissue microarray. Expression of Nrf1 and TFAM was significantly increased in breast cancer tissue compared with adjacent normal tissues. In addition, patients positive for Nrf1 or TFAM had a poorer clinical prognosis than patients who were negative, and those positive for Nrf1 and TFAM had the shortest survival time. These results suggest that Nrf1 and TFAM are potential biomarkers for the determination of individualized therapy and the prognosis of breast cancer, and molecular targeting of Nrf1 and TFAM is a promising strategy for the sensitization of breast cancer cells to chemotherapeutics.

HTLV-1 bZIP factor suppresses TDP1 expression through inhibition of NRF-1 in adult T-cell leukemia

Adult T-cell leukemia (ATL) is an aggressive T-cell malignancy caused by human T-cell leukemia virus type 1 (HTLV-1). We recently reported that abacavir, an anti-HIV-1 drug, potently and selectively kills ATL cells. This effect was attributed to the reduced expression of tyrosyl-DNA-phosphodiesterase 1 (TDP1), a DNA repair enzyme, in ATL cells. However, the molecular mechanism underlying the downregulation of TDP1 in ATL cells remains elusive. Here we identified the core promoter of the TDP1 gene, which contains a conserved nuclear respiratory factor 1 (NRF-1) binding site. Overexpression of NRF-1 increased TDP1-promoter activity, whereas the introduction of dominant-negative NRF-1 repressed such activity. Overexpression of NRF-1 also upregulated endogenous TDP-1 expression, while introduction of shNRF-1 suppressed TDP1 in Jurkat T cells, making them susceptible to abacavir. These results indicate that NRF-1 is a positive transcriptional regulator of TDP1-gene expression. Importantly, we revealed that HTLV-1 bZIP factor (HBZ) protein which is expressed in all ATL cases physically interacts with NRF-1 and inhibits the DNA-binding ability of NRF-1. Taken together, HBZ suppresses TDP1 expression by inhibiting NRF-1 function in ATL cells. The HBZ/NRF-1/TDP1 axis provides new therapeutic targets against ATL and might explain genomic instability leading to the pathogenesis of ATL.

Krüppel-like factors compete for promoters and enhancers to fine-tune transcription

Krüppel-like factors (KLFs) are a family of 17 transcription factors characterized by a conserved DNA-binding domain of three zinc fingers and a variable N-terminal domain responsible for recruiting cofactors. KLFs have diverse functions in stem cell biology, embryo patterning, and tissue homoeostasis. KLF1 and related family members function as transcriptional activators via recruitment of co-activators such as EP300, whereas KLF3 and related members act as transcriptional repressors via recruitment of C-terminal Binding Proteins. KLF1 directly activates the Klf3 gene via an erythroid-specific promoter. Herein, we show KLF1 and KLF3 bind common as well as unique sites within the erythroid cell genome by ChIP-seq. We show KLF3 can displace KLF1 from key erythroid gene promoters and enhancers in vivo. Using 4sU RNA labelling and RNA-seq, we show this competition results in reciprocal transcriptional outputs for >50 important genes. Furthermore, Klf3-/- mice displayed exaggerated recovery from anemic stress and persistent cell cycling consistent with a role for KLF3 in dampening KLF1-driven proliferation. We suggest this study provides a paradigm for how KLFs work in incoherent feed-forward loops or networks to fine-tune transcription and thereby control diverse biological processes such as cell proliferation.

ERRα induces H3K9 demethylation by LSD1 to promote cell invasion

Lysine Specific Demethylase 1 (LSD1) removes mono- and dimethyl groups from lysine 4 of histone H3 (H3K4) or H3K9, resulting in repressive or activating (respectively) transcriptional histone marks. The mechanisms that control the balance between these two antagonist activities are not understood. We here show that LSD1 and the orphan nuclear receptor estrogen-related receptor α (ERRα) display commonly activated genes. Transcriptional activation by LSD1 and ERRα involves H3K9 demethylation at the transcriptional start site (TSS). Strikingly, ERRα is sufficient to induce LSD1 to demethylate H3K9 in vitro. The relevance of this mechanism is highlighted by functional data. LSD1 and ERRα coregulate several target genes involved in cell migration, including the MMP1 matrix metallo-protease, also activated through H3K9 demethylation at the TSS. Depletion of LSD1 or ERRα reduces the cellular capacity to invade the extracellular matrix, a phenomenon that is rescued by MMP1 reexpression. Altogether our results identify a regulatory network involving a direct switch in the biochemical activities of a histone demethylase, leading to increased cell invasion.

A key role of nuclear factor Y in the refeeding response of fatty acid synthase in adipocytes

Fatty acid synthase (Fasn) is a key component of energy metabolism that is dynamically induced by food intake. Although extensive studies have revealed a number of transcription factors involved in the fasting/refeeding transition of Fasn expression in hepatocytes, much less evidence is available for adipocytes. Using the in vivo Ad-luc analytical system, we identified the inverted CCAAT element (ICE) around -100 nucleotides in the Fasn promoter as a critical cis-element for the refeeding response in adipocytes. Electrophoretic mobility shift assays and chromatin immunoprecipitation show that nuclear factor Y (NF-Y) binds to ICE specifically in refeeding states. Notably, the NF-Y binding to ICE is differently regulated between adipocytes and hepatocytes. These findings provide insights into the specific mechanisms controlling energy metabolism in adipocytes.

A high definition look at the NF-Y regulome reveals genome-wide associations with selected transcription factors

NF-Y is a trimeric transcription factor (TF), binding the CCAAT box element, for which several results suggest a pioneering role in activation of transcription. In this work, we integrated 380 ENCODE ChIP-Seq experiments for 154 TFs and cofactors with sequence analysis, protein–protein interactions and RNA profiling data, in order to identify genome-wide regulatory modules resulting from the co-association of NF-Y with other TFs. We identified three main degrees of co-association with NF-Y for sequence-specific TFs. In the most relevant one, we found TFs having a significant overlap with NF-Y in their DNA binding loci, some with a precise spacing of binding sites with respect to the CCAAT box, others (FOS, Sp1/2, RFX5, IRF3, PBX3) mostly lacking their canonical binding site and bound to arrays of well spaced CCAAT boxes. As expected, NF-Y binding also correlates with RNA Pol II General TFs and with subunits of complexes involved in the control of H3K4 methylations. Co-association patterns are confirmed by protein–protein interactions, and correspond to specific functional categorizations and expression level changes of target genes following NF-Y inactivation. These data define genome-wide rules for the organization of NF-Y-centered regulatory modules, supporting a model of distinct categorization and synergy with well defined sets of TFs.

LSD1-Mediated Demethylation of H3K4me2 Is Required for the Transition from Late Progenitor to Differentiated Mouse Rod Photoreceptor

Epigenetic modifiers can work in concert with transcription factors to control the transition of cells from proliferating progenitors into quiescent terminally differentiated cells. This transition involves changes in histone methylation and one of the key regulators of this is the H3K4me2/1 histone demethylase LSD1. Here, we show that the highest expression of LSD1 occurs in postmitotic retinal cells during the peak period of rod photoreceptor differentiation. Pharmacological inhibition of LSD1 in retinal explants cultured from PN1 to PN8 had three major effects. It prevented the normal decrease in expression of genes associated with progenitor function, it blocked rod photoreceptor development, and it increased expression of genes associated with other retinal cell types. The maintained expression of progenitor genes was associated with a maintained level of H3K4me2 over the gene and its promoter. Among the genes whose expression was maintained was Hes1, a repressor known to block rod photoreceptor development. The inhibition of rod photoreceptor gene expression occurred in spite of the normal expression of transcription factors CRX and NRL, and the normal accumulation of H3K4me2 marks over the promoter and gene body. We suggest that LSD1 acts in concert with a series of nuclear receptors to modify chromatin structure and repress progenitor genes as well as to inhibit ectopic patterns of gene expression in the differentiating postmitotic retinal cells.

Genome-wide identification of hypoxia-inducible factor-1 and -2 binding sites in hypoxic human macrophages alternatively activated by IL-10

Macrophages (MΦ) often accumulate in hypoxic areas, where they significantly influence disease progression. Anti-inflammatory cytokines, such as IL-10, generate alternatively activated macrophages that support tumor growth. To understand how alternative activation affects the transcriptional profile of hypoxic macrophages, we globally mapped binding sites of hypoxia-inducible factor (HIF)-1α and HIF-2α in primary human monocyte-derived macrophages prestimulated with IL-10. 713 HIF-1 and 795 HIF-2 binding sites were identified under hypoxia. Pretreatment with IL-10 altered the binding pattern, with 120 new HIF-1 and 188 new HIF-2 binding sites emerging. HIF-1 binding was most prominent in promoters, while HIF-2 binding was more abundant in enhancer regions. Comparison of ChIP-seq data obtained in other cells revealed a highly cell type specific binding of HIF. In MΦ HIF binding occurred preferentially in already active enhancers or promoters. To assess the roles of HIF on gene expression, primary human macrophages were treated with siRNA against HIF-1α or HIF-2α, followed by genome-wide gene expression analysis. Comparing mRNA expression to the HIF binding profile revealed a significant enrichment of hypoxia-inducible genes previously identified by ChIP-seq. Analysis of gene expression under hypoxia alone and hypoxia/IL-10 showed the enhanced induction of a set of genes including PLOD2 and SLC2A3, while another group including KDM3A and ADM remained unaffected or was reduced by IL-10. Taken together IL-10 influences the DNA binding pattern of HIF and the level of gene induction.

LSD1 promotes oxidative metabolism of white adipose tissue

Exposure to environmental cues such as cold or nutritional imbalance requires white adipose tissue (WAT) to adapt its metabolism to ensure survival. Metabolic plasticity is prominently exemplified by the enhancement of mitochondrial biogenesis in WAT in response to cold exposure or β3-adrenergic stimulation. Here we show that these stimuli increase the levels of lysine-specific demethylase 1 (LSD1) in WAT of mice and that elevated LSD1 levels induce mitochondrial activity. Genome-wide binding and transcriptome analyses demonstrate that LSD1 directly stimulates the expression of genes involved in oxidative phosphorylation (OXPHOS) in cooperation with nuclear respiratory factor 1 (Nrf1). In transgenic (Tg) mice, increased levels of LSD1 promote in a cell-autonomous manner the formation of islets of metabolically active brown-like adipocytes in WAT. Notably, Tg mice show limited weight gain when fed a high-fat diet. Taken together, our data establish LSD1 as a key regulator of OXPHOS and metabolic adaptation in WAT.