High FOXA1 levels induce ER transcriptional reprogramming, a pro-metastatic secretome, and metastasis in endocrine-resistant breast cancer

Aberrant activation of the forkhead protein FOXA1 is observed in advanced hormone-related cancers. However, the key mediators of high FOXA1 signaling remain elusive. We demonstrate that ectopic high FOXA1 (H-FOXA1) expression promotes estrogen receptor-positive (ER+) breast cancer (BC) metastasis in a xenograft mouse model. Mechanistically, H-FOXA1 reprograms ER-chromatin binding to elicit a core gene signature (CGS) enriched in ER+ endocrine-resistant (EndoR) cells. We identify Secretome14, a CGS subset encoding ER-dependent cancer secretory proteins, as a strong predictor for poor outcomes of ER+ BC. It is elevated in ER+ metastases vs. primary tumors, irrespective of ESR1 mutations. Genomic ER binding near Secretome14 genes is also increased in mutant ER-expressing or mitogen-treated ER+ BC cells and in ER+ metastatic vs. primary tumors, suggesting a convergent pathway including high growth factor receptor signaling in activating pro-metastatic secretome genes. Our findings uncover H-FOXA1-induced ER reprogramming that drives EndoR and metastasis partly via an H-FOXA1/ER-dependent secretome.

HNF4A modulates glucocorticoid action in the liver

The glucocorticoid receptor (GR) is a nuclear receptor critical to the regulation of energy metabolism and inflammation. The actions of GR are dependent on cell type and context. Here, we demonstrate the role of liver lineage-determining factor hepatocyte nuclear factor 4A (HNF4A) in defining liver specificity of GR action. In mouse liver, the HNF4A motif lies adjacent to the glucocorticoid response element (GRE) at GR binding sites within regions of open chromatin. In the absence of HNF4A, the liver GR cistrome is remodeled, with loss and gain of GR recruitment evident. Loss of chromatin accessibility at HNF4A-marked sites associates with loss of GR binding at weak GRE motifs. GR binding and chromatin accessibility are gained at sites characterized by strong GRE motifs, which show GR recruitment in non-liver tissues. The functional importance of these HNF4A-regulated GR sites is indicated by an altered transcriptional response to glucocorticoid treatment in the Hnf4a-null liver.

ERα and ERβ Homodimers in the Same Cellular Context Regulate Distinct Transcriptomes and Functions

The two estrogen receptors ERα and ERβ are nuclear receptors that bind estrogen (E2) and function as ligand-inducible transcription factors. They are homologues and can form dimers with each other and bind to the same estrogen-response element motifs in the DNA. ERα drives breast cancer growth whereas ERβ has been reported to be anti-proliferative. However, they are rarely expressed in the same cells, and it is not fully investigated to which extent their functions are different because of inherent differences or because of different cellular context. To dissect their similarities and differences, we here generated a novel estrogen-dependent cell model where ERα homodimers can be directly compared to ERβ homodimers within the identical cellular context. By using CRISPR-cas9 to delete ERα in breast cancer MCF7 cells with Tet-Off-inducible ERβ expression, we generated MCF7 cells that express ERβ but not ERα. MCF7 (ERβ only) cells exhibited regulation of estrogen-responsive targets in a ligand-dependent manner. We demonstrated that either ER was required for MCF7 proliferation, but while E2 increased proliferation via ERα, it reduced proliferation through a G2/M arrest via ERβ. The two ERs also impacted migration differently. In absence of ligand, ERβ increased migration, but upon E2 treatment, ERβ reduced migration. E2 via ERα, on the other hand, had no significant impact on migration. RNA sequencing revealed that E2 regulated a transcriptome of around 800 genes via each receptor, but over half were specific for either ERα or ERβ (417 and 503 genes, respectively). Functional gene ontology enrichment analysis reinforced that E2 regulated cell proliferation in opposite directions depending on the ER, and that ERβ specifically impacted extracellular matrix organization. We corroborated that ERβ bound to cis-regulatory chromatin of its unique proposed migration-related direct targets ANXA9 and TFAP2C. In conclusion, we demonstrate that within the same cellular context, the two ERs regulate cell proliferation in the opposite manner, impact migration differently, and each receptor also regulates a distinct set of target genes in response to E2. The developed cell model provides a novel and valuable resource to further complement the mechanistic understanding of the two different ER isoforms.

Epigenetic and transcriptional analysis reveals a core transcriptional program conserved in clonal prostate cancer metastases

The epigenomic regulation of transcriptional programs in metastatic prostate cancer is poorly understood. We studied the epigenomic landscape of prostate cancer drivers using transcriptional profiling and ChIP‐seq in four clonal metastatic tumors derived from a single prostate cancer patient. Our epigenomic analyses focused on androgen receptor (AR), which is a key oncogenic driver in prostate cancer, the AR pioneer factor FOXA1, chromatin insulator CCCTC‐Binding Factor, as well as for modified histones H3K27ac and H3K27me3. The vast majority of AR binding sites were shared among healthy prostate, primary prostate cancer, and metastatic tumor samples, signifying core AR‐driven transcriptional regulation within the prostate cell lineage. Genes associated with core AR‐binding events were significantly enriched for essential genes in prostate cancer cell proliferation. Remarkably, the metastasis‐specific active AR binding sites showed no differential transcriptional output, indicating a robust transcriptional program across metastatic samples. Combined, our data reveal a core transcriptional program in clonal metastatic prostate cancer, despite epigenomic differences in the AR cistrome.

A temporal hierarchy underpins the transcription factor-DNA interactome of the maize UPR

Adverse environmental conditions reduce crop productivity and often increase the load of unfolded or misfolded proteins in the endoplasmic reticulum (ER). This potentially lethal condition, known as ER stress, is buffered by the unfolded protein response (UPR), a set of signaling pathways designed to either recover ER functionality or ignite programmed cell death. Despite the biological significance of the UPR to the life of the organism, the regulatory transcriptional landscape underpinning ER stress management is largely unmapped, especially in crops. To fill this significant knowledge gap, we performed a large-scale systems-level analysis of the protein-DNA interaction (PDI) network in maize (Zea mays). Using 23 promoter fragments of six UPR marker genes in a high-throughput enhanced yeast one-hybrid assay, we identified a highly interconnected network of 262 transcription factors (TFs) associated with significant biological traits and 831 PDIs underlying the UPR. We established a temporal hierarchy of TF binding to gene promoters within the same family as well as across different families of TFs. Cistrome analysis revealed the dynamic activities of a variety of cis-regulatory elements (CREs) in ER stress-responsive gene promoters. By integrating the cistrome results into a TF network analysis, we mapped a subnetwork of TFs associated with a CRE that may contribute to UPR management. Finally, we validated the role of a predicted network hub gene using the Arabidopsis system. The PDIs, TF networks, and CREs identified in our work are foundational resources for understanding transcription-regulatory mechanisms in the stress responses and crop improvement.

Darolutamide antagonizes androgen signaling by blocking enhancer and super-enhancer activation

Prostate cancer (PCa) is one of the most frequent tumor types in the male Western population. Early-stage PCa and late-stage PCa are dependent on androgen signaling, and inhibitors of the androgen receptor (AR) axis represent the standard therapy. Here, we studied in detail the global impact of darolutamide, a newly approved AR antagonist, on the transcriptome and AR-bound cistrome in two PCa cell models. Darolutamide strongly depleted the AR from gene regulatory regions and abolished AR-driven transcriptional signaling. Enhancer activation was blocked at the chromatin level as evaluated by H3K27 acetylation (H3K27ac), H3K4 monomethylation (H3K4me1), and FOXA1, MED1, and BRD4 binding. We identified genomic regions with high affinities for the AR in androgen-stimulated, but also in androgen-depleted conditions. A similar AR affinity pattern was observed in healthy and PCa tissue samples. High FOXA1, BRD4, H3K27ac, and H3K4me1 levels were found to mark regions showing AR binding in the hormone-depleted setting. Conversely, low FOXA1, BRD4, and H3K27ac levels were observed at regulatory sites that responded strongly to androgen stimulation, and AR interactions at these sites were blocked by darolutamide. Beside marked loss of AR occupancy, FOXA1 recruitment to chromatin was also clearly reduced after darolutamide treatment. We furthermore identified numerous androgen-regulated super-enhancers (SEs) that were associated with hallmark androgen and cell proliferation-associated gene sets. Importantly, these SEs are also active in PCa tissues and sensitive to darolutamide treatment in our models. Our findings demonstrate that darolutamide is a potent AR antagonist blocking genome-wide AR enhancer and SE activation, and downstream transcription. We also show the existence of a dynamic AR cistrome that depends on the androgen levels and on high AR affinity regions present in PCa cell lines and also in tissue samples.

Human Endometrial Transcriptome and Progesterone Receptor Cistrome Reveal Important Pathways and Epithelial Regulators

CONTEXT

Poor uterine receptivity is one major factor leading to pregnancy loss and infertility. Understanding the molecular events governing successful implantation is hence critical in combating infertility.

OBJECTIVE

To define Progesterone Receptor (PGR)-regulated molecular mechanisms and epithelial roles in receptivity.

DESIGN

RNA-sequencing and PGR-ChIP-seq were conducted in parallel to identify PGR-regulated pathways during the Window of implantation (WOI) in endometrium of fertile women.

SETTING

Endometrial biopsies from the proliferative and mid-secretory phases were analyzed.

PATIENTS OR OTHER PARTICIPANTS

Participants were fertile, reproductive aged (18-37 years) women with normal cycle length, and without any history of dysmenorrhea, infertility, or irregular cycles. In total, 42 endometrial biopsies obtained from 42 women were analyzed in this study.

INTERVENTIONS

There were no interventions during this study.

MAIN OUTCOME MEASURES

Here we measured the alterations in gene expression and PGR occupancy in the genome during the WOI, based on the hypothesis that PGR binds uterine chromatin cycle dependently to regulate genes involved in uterine cell differentiation and function.

RESULTS

653 genes were identified with regulated PGR binding and differential expression during the WOI. These were involved in regulating inflammatory response, xenobiotic metabolism, epithelial mesenchymal transition, cell death, interleukin/Signal Transducer And Activator Of Transcription (STAT) signaling, estrogen response, and Mammalian target of rapamycin complex 1 (MTORC1) response. Transcriptome of the epithelium identified 3052 differentially expressed genes, of which 658 were uniquely regulated. Transcription factors Interferon Regulatory Factor 8 (IRF8) and Myocyte Enhancer Factor 2C (MEF2C) were found to be regulated in the epithelium during the WOI at the protein level, suggesting potentially important functions that are previously unrecognized.

CONCLUSION

PGR binds the genomic regions of genes regulating critical processes in uterine receptivity and function.

Modulation of androgen receptor DNA binding activity through direct interaction with the ETS transcription factor ERG

Progress in studying the androgen receptor (AR), the primary drug target in prostate cancer, has been hampered by challenges in expressing and purifying active multidomain AR for use in cell-free biochemical reconstitution assays. Here we successfully express full-length and truncated AR variants and demonstrate that the oncogenic ETS protein ERG, responsible for half of all prostate cancers, enhances the ability of AR to bind DNA through direct interaction with AR. In addition to providing a biochemical system to evaluate AR activity on different DNA templates, our findings provide insight into why ERG-positive prostate cancers have an expanded AR cistrome.

The androgen receptor (AR) is a type I nuclear hormone receptor and the primary drug target in prostate cancer due to its role as a lineage survival factor in prostate luminal epithelium. In prostate cancer, the AR cistrome is reprogrammed relative to normal prostate epithelium and particularly in cancers driven by oncogenic ETS fusion genes. The molecular basis for this change has remained elusive. Using purified proteins, we report a minimal cell-free system that demonstrates interdomain cooperativity between the ligand (LBD) and DNA binding domains (DBD) of AR, and its autoinhibition by the N terminus of AR. Furthermore, we identify ERG as a cofactor that activates AR’s ability to bind DNA in both high and lower affinity contexts through direct interaction within a newly identified AR-interacting motif (AIM) in the ETS domain, independent of ERG’s own DNA binding ability. Finally, we present evidence that this interaction is conserved among ETS factors whose expression is altered in prostate cancer. Our work highlights, at a biochemical level, how tumor-initiating ETS translocations result in reprogramming of the AR cistrome.

Deciphering essential cistromes using genome-wide CRISPR screens

Systematically dissecting the function of a large set of cis-regulatory elements or transcription factor binding sites (cistromes) has been technically challenging. Using genome-wide CRISPR screens, we profiled over 10,000 FOXA1 and CTCF binding sites for their roles in regulating the fitness of breast and prostate cancer cells, and accordingly developed a model to predict essentiality for cis-elements. These efforts not only reveal how the key transcription factors and their cistromes regulate cell essentiality in hormone-dependent cancers but also highlight an efficient approach to investigate the functions of noncoding regions of the genome.

Although millions of transcription factor binding sites, or cistromes, have been identified across the human genome, defining which of these sites is functional in a given condition remains challenging. Using CRISPR/Cas9 knockout screens and gene essentiality or fitness as the readout, we systematically investigated the essentiality of over 10,000 FOXA1 and CTCF binding sites in breast and prostate cancer cells. We found that essential FOXA1 binding sites act as enhancers to orchestrate the expression of nearby essential genes through the binding of lineage-specific transcription factors. In contrast, CRISPR screens of the CTCF cistrome revealed 2 classes of essential binding sites. The first class of essential CTCF binding sites act like FOXA1 sites as enhancers to regulate the expression of nearby essential genes, while a second class of essential CTCF binding sites was identified at topologically associated domain (TAD) boundaries and display distinct characteristics. Using regression methods trained on our screening data and public epigenetic profiles, we developed a model to predict essential cis-elements with high accuracy. The model for FOXA1 essentiality correctly predicts noncoding variants associated with cancer risk and progression. Taken together, CRISPR screens of cis-regulatory elements can define the essential cistrome of a given factor and can inform the development of predictive models of cistrome function.

The Lineage Determining Factor GRHL2 Collaborates with FOXA1 to Establish a Targetable Pathway in Endocrine Therapy-Resistant Breast Cancer

Notwithstanding the positive clinical impact of endocrine therapies in estrogen receptor-alpha (ERα)-positive breast cancer, de novo and acquired resistance limits the therapeutic lifespan of existing drugs. Taking the position that resistance is nearly inevitable, we undertook a study to identify and exploit targetable vulnerabilities that were manifest in endocrine therapy-resistant disease. Using cellular and mouse models of endocrine therapy-sensitive and endocrine therapy-resistant breast cancer, together with contemporary discovery platforms, we identified a targetable pathway that is composed of the transcription factors FOXA1 and GRHL2, a coregulated target gene, the membrane receptor LYPD3, and the LYPD3 ligand, AGR2. Inhibition of the activity of this pathway using blocking antibodies directed against LYPD3 or AGR2 inhibits the growth of endocrine therapy-resistant tumors in mice, providing the rationale for near-term clinical development of humanized antibodies directed against these proteins.

Factors that influence the androgen receptor cistrome in benign and malignant prostate cells

The androgen receptor (AR) plays key roles in the development of prostate tissue and the development and progression of prostate cancer (PC). AR guides cytodifferentiation and homeostasis in benign luminal epithelial cells; however, in PC, AR instead drives the uncontrolled proliferation of these cells. This ‘AR malignancy shift’ (AMS) is a central event in tumorigenesis. Using a ChIP‐seq approach in primary human tissues, cell lines, and mouse models, we demonstrate that the AMS occurs in every sample analyzed, suggesting that it is necessary for PC development. Using molecular and genetic techniques, we demonstrate that forkhead box (FOX)A1, HOXB13, GATA2, and c‐JUN are involved in the regulation of the AMS. AR‐binding sites (ARBS) are enriched for FOX, HOX, and GATA motifs in PC cells but not for c‐JUN motifs in benign cells. We show that the SPOP mutation commonly found in localized PCs can cause the AMS but is not transformative on its own and must be coupled to another mutation to transform cells. We show that the AMS occurs in mouse models of PC as well and that chronic low T, which is associated with increased PC risk and aggressiveness in humans, also causes the AMS in mice. We have discovered a previously unrecognized, fundamental tenet of PC, one which explains how and why AR signaling is different in cancer and benign cells. Our work has the potential to be used to stratify patients with localized PC for specific treatments. Furthermore, our work suggests that the AMS is a novel target for the treatment and/or prevention of PC.

Androgen receptor enhancer usage and the chromatin regulatory landscape in human prostate cancers

The androgen receptor (AR) is commonly known as a key transcription factor in prostate cancer development, progression and therapy resistance. Genome-wide chromatin association studies revealed that transcriptional regulation by AR mainly depends on binding to distal regulatory enhancer elements that control gene expression through chromatin looping to gene promoters. Changes in the chromatin epigenetic landscape and DNA sequence can locally alter AR-DNA-binding capacity and consequently impact transcriptional output and disease outcome. The vast majority of reports describing AR chromatin interactions have been limited to cell lines, identifying numerous other factors and interacting transcription factors that impact AR chromatin interactions. Do these factors also impact AR cistromics - the genome-wide chromatin-binding landscape of AR - in vivo? Recent technological advances now enable researchers to identify AR chromatin-binding sites and their target genes in human specimens. In this review, we provide an overview of the different factors that influence AR chromatin binding in prostate cancer specimens, which is complemented with knowledge from cell line studies. Finally, we discuss novel perspectives on studying AR cistromics in clinical samples.

Gene networks orchestrated by MeGI: a single-factor mechanism underlying sex determination in persimmon

Separating male and female sex organs is one of the main strategies used to maintain genetic diversity within a species. However, the genetic determinants and their regulatory mechanisms have been identified in only a few species. In dioecious persimmons, the homeodomain transcription factor, MeGI, which is the target of a Y chromosome-encoded small-RNA, OGI, can determine floral sexuality. The basic features of this system are conserved in the monoecious hexaploid Oriental persimmon, in which an additional epigenetic regulation of MeGI determines floral sexuality. The downstream regulatory pathways of MeGI remain uncharacterized. In this study, we examined transcriptomic data for male and female flowers from monoecious persimmon cultivars to unveil the gene networks orchestrated by MeGI. A network visualization and cistrome assessment suggested that class-1 KNOTTED-like homeobox (KNOX)/ovate family protein (OFP)/growth regulating factors (GRFs) and short vegetative phase (SVP) genes mediate the differences in gynoecium and androecium development between male and female flowers, respectively. The expression of these genes is directly controlled by MeGI. The gene networks also suggested that some cytokinin, auxin, and gibberellin signaling genes function cooperatively in the KNOX/OFP/GRF pathway during gynoecium differentiation. Meanwhile, SVP may repress PI expression in developing androecia. Overall, our results suggest that MeGI evolved the ability to promote gynoecium development and suppress androecium development by regulating KNOX/OFP/GRF and SVP expression levels, respectively. These insights may help to clarify the molecular mechanism underlying the production of unisexual flowers, while also elucidating the physiological background enabling a single-factor system to establish dioecy in plants.

Src promotes castration-recurrent prostate cancer through androgen receptor-dependent canonical and non-canonical transcriptional signatures

Progression of prostate cancer (PC) to castration-recurrent growth (CRPC) remains dependent on sustained expression and transcriptional activity of the androgen receptor (AR). A major mechanism contributing to CRPC progression is through the direct phosphorylation and activation of AR by Src-family (SFK) and ACK1 tyrosine kinases. However, the AR-dependent transcriptional networks activated by Src during CRPC progression have not been elucidated. Here, we show that activated Src (Src527F) induces androgen-independent growth in human LNCaP cells, concomitant with its ability to induce proliferation/survival genes normally induced by dihydrotestosterone (DHT) in androgen-dependent LNCaP and VCaP cells. Src induces additional gene signatures unique to CRPC cell lines, LNCaP-C4-2 and CWR22Rv1, and to CRPC LuCaP35.1 xenografts. By comparing the Src-induced AR-cistrome and/or transcriptome in LNCaP to those in CRPC and LuCaP35.1 tumors, we identified an 11-gene Src-regulated CRPC signature consisting of AR-dependent, AR binding site (ARBS)-associated genes whose expression is altered by DHT in LNCaP[Src527F] but not in LNCaP cells. The differential expression of a subset (DPP4, BCAT1, CNTNAP4, CDH3) correlates with earlier PC metastasis onset and poorer survival, with the expression of BCAT1 required for Src-induced androgen-independent proliferation. Lastly, Src enhances AR binding to non-canonical ARBS enriched for FOXO1, TOP2B and ZNF217 binding motifs; cooperative AR/TOP2B binding to a non-canonical ARBS was both Src- and DHT-sensitive and correlated with increased levels of Src-induced phosphotyrosyl-TOP2B. These data suggest that CRPC progression is facilitated via Src-induced sensitization of AR to intracrine androgen levels, resulting in the engagement of canonical and non-canonical ARBS-dependent gene signatures.

Allele-Specific Chromatin Recruitment and Therapeutic Vulnerabilities of ESR1 Activating Mutations

Estrogen receptor α (ER) ligand-binding domain (LBD) mutations are found in a substantial number of endocrine treatment-resistant metastatic ER-positive (ER⁺) breast cancers. We investigated the chromatin recruitment, transcriptional network, and genetic vulnerabilities in breast cancer models harboring the clinically relevant ER mutations. These mutants exhibit both ligand-independent functions that mimic estradiol-bound wild-type ER as well as allele-specific neomorphic properties that promote a pro-metastatic phenotype. Analysis of the genome-wide ER binding sites identified mutant ER unique recruitment mediating the allele-specific transcriptional program. Genetic screens identified genes that are essential for the ligand-independent growth driven by the mutants. These studies provide insights into the mechanism of endocrine therapy resistance engendered by ER mutations and potential therapeutic targets.

Embryonic transcription factor SOX9 drives breast cancer endocrine resistance

The estrogen receptor (ER) drives the growth of most luminal breast cancers and is the primary target of endocrine therapy. Although ER blockade with drugs such as tamoxifen is very effective, a major clinical limitation is the development of endocrine resistance especially in the setting of metastatic disease. Preclinical and clinical observations suggest that even following the development of endocrine resistance, ER signaling continues to exert a pivotal role in tumor progression in the majority of cases. Through the analysis of the ER cistrome in tamoxifen-resistant breast cancer cells, we have uncovered a role for an RUNX2-ER complex that stimulates the transcription of a set of genes, including most notably the stem cell factor SOX9, that promote proliferation and a metastatic phenotype. We show that up-regulation of SOX9 is sufficient to cause relative endocrine resistance. The gain of SOX9 as an ER-regulated gene associated with tamoxifen resistance was validated in a unique set of clinical samples supporting the need for the development of improved ER antagonists.

The ubiquitous transcription factor CTCF promotes lineage-specific epigenomic remodeling and establishment of transcriptional networks driving cell differentiation

Cell differentiation relies on tissue-specific transcription factors (TFs) that cooperate to establish unique transcriptomes and phenotypes. However, the role of ubiquitous TFs in these processes remains poorly defined. Recently, we have shown that the CCCTC-binding factor (CTCF) is required for adipocyte differentiation through epigenomic remodelling of adipose tissue-specific enhancers and transcriptional activation of Peroxisome proliferator-activated receptor gamma (PPARG), the main driver of the adipogenic program (PPARG), and its target genes. Here, we discuss how these findings, together with the recent literature, illuminate a functional role for ubiquitous TFs in lineage-determining transcriptional networks.

Biological functional annotation of retinoic acid alpha and beta in mouse liver based on genome-wide binding

Retinoic acid (RA) has diverse biological effects. The liver stores vitamin A, generates RA, and expresses receptors for RA. The current study examines the hepatic binding profile of two RA receptor isoforms, RARA (RARα) and RARB (RARβ), in response to RA treatment in mouse livers. Our data uncovered 35,521, and 14,968 genomic bindings for RARA and RARB, respectively. Each expressed unique and common bindings, implying their redundant and specific roles. RARB has higher RA responsiveness than RARB. RA treatment generated 18,821 novel RARB bindings but only 14,798 of RARA bindings, compared with the control group. RAR frequently bound the consensus hormone response element [HRE; (A/G)G(G/T)TCA], which often contained the motifs assigned to SP1, GABPA, and FOXA2, suggesting potential interactions between those transcriptional factors. Functional annotation coupled with principle component analysis revealed that the function of RAR target genes were motif dependent. Taken together, the cistrome of RARA and RARB revealed their extensive biological roles in the mouse liver. RAR target genes are enriched in various biological processes. The hepatic RAR genome-wide binding data can help us understand the global molecular mechanisms underlying RAR and RA-mediated gene and pathway regulation.

Using MACS to identify peaks from ChIP-Seq data

Model-based Analysis of ChIP-Seq (MACS) is a command-line tool designed by X. Shirley Liu and colleagues to analyze data generated by ChIP-Seq experiments in eukaryotes, especially mammals. MACS can be used to identify transcription factor binding sites and histone modification-enriched regions if the ChIP-Seq data, with or without control samples, are given. This unit describes two basic protocols that provide detailed information on how to use MACS to identify either the binding sites of a transcription factor or the enriched regions of a histone modification with broad peaks. Furthermore, the basic ideas for the MACS algorithm and its appropriate usage are discussed.