Hierarchical mechanisms for direct reprogramming of fibroblasts to neurons.

Direct lineage reprogramming is a promising approach for human disease modeling and regenerative medicine, with poorly understood mechanisms. Here, we reveal a hierarchical mechanism in the direct conversion of fibroblasts into induced neuronal (iN) cells mediated by the transcription factors Ascl1, Brn2, and Myt1l. Ascl1 acts as an "on-target" pioneer factor by immediately occupying most cognate genomic sites in fibroblasts. In contrast, Brn2 and Myt1l do not access fibroblast chromatin productively on their own; instead, Ascl1 recruits Brn2 to Ascl1 sites genome wide. A unique trivalent chromatin signature in the host cells predicts the permissiveness for Ascl1 pioneering activity among different cell types. Finally, we identified Zfp238 as a key Ascl1 target gene that can partially substitute for Ascl1 during iN cell reprogramming. Thus, a precise match between pioneer factors and the chromatin context at key target genes is determinative for transdifferentiation to neurons and likely other cell types.

Highly transcribed RNA polymerase II genes are impediments to replication fork progression in Saccharomyces cerevisiae.

Replication forks face multiple obstacles that slow their progression. By two-dimensional gel analysis, yeast forks pause at stable DNA protein complexes, and this pausing is greatly increased in the absence of the Rrm3 helicase. We used a genome-wide approach to identify 96 sites of very high DNA polymerase binding in wild-type cells. Most of these binding sites were not previously identified pause sites. Rather, the most highly represented genomic category among high DNA polymerase binding sites was the open reading frames (ORFs) of highly transcribed RNA polymerase II genes. Twice as many pause sites were identified in rrm3 compared with wild-type cells, as pausing in this strain occurred at both highly transcribed RNA polymerase II genes and the previously identified protein DNA complexes. ORFs of highly transcribed RNA polymerase II genes are a class of natural pause sites that are not exacerbated in rrm3 cells.

Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position.

We describe an assay for transposase-accessible chromatin using sequencing (ATAC-seq), based on direct in vitro transposition of sequencing adaptors into native chromatin, as a rapid and sensitive method for integrative epigenomic analysis. ATAC-seq captures open chromatin sites using a simple two-step protocol with 500-50,000 cells and reveals the interplay between genomic locations of open chromatin, DNA-binding proteins, individual nucleosomes and chromatin compaction at nucleotide resolution. We discovered classes of DNA-binding factors that strictly avoided, could tolerate or tended to overlap with nucleosomes. Using ATAC-seq maps of human CD4(+) T cells from a proband obtained on consecutive days, we demonstrated the feasibility of analyzing an individual's epigenome on a timescale compatible with clinical decision-making.

X chromosome repression by localization of the C. elegans dosage compensation machinery to sites of transcription initiation.

Among organisms with chromosome-based mechanisms of sex determination, failure to equalize expression of X-linked genes between the sexes is typically lethal. In C. elegans, XX hermaphrodites halve transcription from each X chromosome to match the output of XO males. Here, we mapped the binding location of the condensin homolog DPY-27 and the zinc finger protein SDC-3, two components of the C. elegans dosage compensation complex (DCC). We observed strong foci of DCC binding on X, surrounded by broader regions of localization. Binding foci, but not adjacent regions of localization, were distinguished by clusters of a 10-bp DNA motif, suggesting a recruitment-and-spreading mechanism for X recognition. The DCC was preferentially bound upstream of genes, suggesting modulation of transcriptional initiation and polymerase-coupled spreading. Stronger DCC binding upstream of genes with high transcriptional activity indicated a mechanism for tuning DCC activity at specific loci. These data aid in understanding how proteins involved in higher-order chromosome dynamics can regulate transcription at individual loci.

Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity.

The human body contains thousands of unique cell types, each with specialized functions. Cell identity is governed in large part by gene transcription programs, which are determined by regulatory elements encoded in DNA. To identify regulatory elements active in seven cell lines representative of diverse human cell types, we used DNase-seq and FAIRE-seq (Formaldehyde Assisted Isolation of Regulatory Elements) to map "open chromatin." Over 870,000 DNaseI or FAIRE sites, which correspond tightly to nucleosome-depleted regions, were identified across the seven cell lines, covering nearly 9% of the genome. The combination of DNaseI and FAIRE is more effective than either assay alone in identifying likely regulatory elements, as judged by coincidence with transcription factor binding locations determined in the same cells. Open chromatin common to all seven cell types tended to be at or near transcription start sites and to be coincident with CTCF binding sites, while open chromatin sites found in only one cell type were typically located away from transcription start sites and contained DNA motifs recognized by regulators of cell-type identity. We show that open chromatin regions bound by CTCF are potent insulators. We identified clusters of open regulatory elements (COREs) that were physically near each other and whose appearance was coordinated among one or more cell types. Gene expression and RNA Pol II binding data support the hypothesis that COREs control gene activity required for the maintenance of cell-type identity. This publicly available atlas of regulatory elements may prove valuable in identifying noncoding DNA sequence variants that are causally linked to human disease.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Regulation of nucleosome stability as a mediator of chromatin function.

Interactions among regulatory proteins, enzymes and DNA are required to use the information encoded in genomes. In eukaryotes, the location and timing of interactions between these proteins and DNA is determined in large part by whether DNA at a given locus is packaged into a nucleosome. Given the central role of nucleosome formation in regulating genome function, many mechanisms have evolved to control nucleosome stability at loci across the genome. New conclusions about the role of the DNA sequence itself in the regulation of nucleosome stability have come from two new types of experiment: high-resolution mapping of in vivo nucleosome occupancy on a genomic scale; and in vivo versus in vitro comparisons of nucleosome stability on natural DNA templates. These new types of data raise intriguing questions about the evolutionary constraints on DNA sequence with regard to nucleosome formation, and at long last might enable the derivation of DNA sequence-based rules that produce reliable predictions of nucleosome behavior.

Isolation of active regulatory elements from eukaryotic chromatin using FAIRE (Formaldehyde Assisted Isolation of Regulatory Elements).

The binding of sequence-specific regulatory factors and the recruitment of chromatin remodeling activities cause nucleosomes to be evicted from chromatin in eukaryotic cells. Traditionally, these active sites have been identified experimentally through their sensitivity to nucleases. Here we describe the details of a simple procedure for the genome-wide isolation of nucleosome-depleted DNA from human chromatin, termed FAIRE (Formaldehyde Assisted Isolation of Regulatory Elements). We also provide protocols for different methods of detecting FAIRE-enriched DNA, including use of PCR, DNA microarrays, and next-generation sequencing. FAIRE works on all eukaryotic chromatin tested to date. To perform FAIRE, chromatin is crosslinked with formaldehyde, sheared by sonication, and phenol-chloroform extracted. Most genomic DNA is crosslinked to nucleosomes and is sequestered to the interphase, whereas DNA recovered in the aqueous phase corresponds to nucleosome-depleted regions of the genome. The isolated regions are largely coincident with the location of DNaseI hypersensitive sites, transcriptional start sites, enhancers, insulators, and active promoters. Given its speed and simplicity, FAIRE has utility in establishing chromatin profiles of diverse cell types in health and disease, isolating DNA regulatory elements en masse for further characterization, and as a screening assay for the effects of small molecules on chromatin organization.

Identification of human CD4+ T cell populations with distinct antitumor activity.

How naturally arising human CD4+ T helper subsets affect cancer immunotherapy is unclear. We reported that human CD4+CD26high T cells elicit potent immunity against solid tumors. As CD26high T cells are often categorized as TH17 cells for their IL-17 production and high CD26 expression, we posited these populations would have similar molecular properties. Here, we reveal that CD26high T cells are epigenetically and transcriptionally distinct from TH17 cells. Of clinical importance, CD26high and TH17 cells engineered with a chimeric antigen receptor (CAR) regressed large human tumors to a greater extent than enriched TH1 or TH2 cells. Only human CD26high T cells mediated curative responses, even when redirected with a suboptimal CAR and without aid by CD8+ CAR T cells. CD26high T cells cosecreted effector cytokines, produced cytotoxic molecules, and persisted long term. Collectively, our work underscores the promise of CD4+ T cell populations to improve durability of solid tumor therapies.

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8+ T cell exhaustion and CD4+ T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Chromatin Accessibility Landscape of Cutaneous T Cell Lymphoma and Dynamic Response to HDAC Inhibitors.

Here, we define the landscape and dynamics of active regulatory DNA in cutaneous T cell lymphoma (CTCL) by ATAC-seq. Analysis of 111 human CTCL and control samples revealed extensive chromatin signatures that distinguished leukemic, host, and normal CD4+ T cells. We identify three dominant patterns of transcription factor (TF) activation that drive leukemia regulomes, as well as TF deactivations that alter host T cells in CTCL patients. Clinical response to histone deacetylase inhibitors (HDACi) is strongly associated with a concurrent gain in chromatin accessibility. HDACi causes distinct chromatin responses in leukemic and host CD4+ T cells, reprogramming host T cells toward normalcy. These results provide a foundational framework to study personal regulomes in human cancer and epigenetic therapy.

Identification of a molecular signature of sarcopenia.

Investigating the molecular mechanisms underlying sarcopenia in humans with the use of microarrays has been complicated by low sample size and the variability inherent in human gene expression profiles. We have conducted a study using Affymetrix GeneChips to identify a molecular signature of aged skeletal muscle. The molecular signature was defined as the set of expressed genes that best distinguished the vastus lateralis muscle of young (n = 10) and older (n = 12) male subjects, when a k-nearest neighbor supervised classification method was used in conjunction with a signal-to-noise ratio gene selection method and a holdout cross-validation procedure. The age-specific expression signature was comprised of 45 genes; 27 were upregulated and 18 were downregulated. This signature also correctly classified 75% of the muscle samples from young and older subjects published by an independent laboratory, based on their expression profiles. The signature revealed increased expression of several genes involved in mediating cellular responses to inflammation and apoptosis, including complement component C1QA, Galectin-1, C/EBP-beta, and FOXO3A, among others. The increased expressions of genes that regulate pre-mRNA splicing, localization, and modification of RNA comprise markers of the aging signature. Downregulated genes in the signature were the glutamine transporter SLC38A1, a TRAF-6 inhibitory zinc finger protein, and membrane-bound transcription factor protease S2P, among others. The sarcopenia signature developed here will be useful as a molecular model to judge the effectiveness of exercise and other therapeutic treatments aimed at ameliorating the effects of muscle loss associated with aging.

Transcriptional profile of a myotube starvation model of atrophy.

Skeletal muscle wasting is a pervasive phenomenon that can result from a wide range of pathological conditions as well as from habitual muscular inactivity. The present work describes a cell-culture condition that induces significant atrophy in skeletal muscle C2C12 myotubes. The failure to replenish differentiation media in mature myotubes leads to rapid atrophy (53% in diameter), which is referred to here as starvation. Affymetrix microarrays were used to develop a transcriptional profile of control (fed) vs. atrophied (nonfed) myotubes. Myotube starvation was characterized by an upregulation of genes involved in translational inhibition, amino acid biosynthesis and transport, and cell cycle arrest/apoptosis, among others. Downregulated genes included several structural and regulatory elements of the extracellular matrix as well as several elements of Wnt/frizzled and TGF-beta signaling pathways. Interestingly, the characteristic transcriptional upregulation of the ubiquitin-proteasome system, calpains, and cathepsins known to occur in multiple in vivo models of atrophy were not seen during myotube starvation. With the exception of the downregulation of extracellular matrix genes, serine protease inhibitor genes, and the upregulation of the translation initiation factor PHAS-I, this model of atrophy in cell culture has a transcriptional profile quite distinct from any study published to date with atrophy in whole muscle. These data show that, although the gross morphology of atrophied muscle fibers may be similar in whole muscle vs. myotube culture, the processes by which this phenotype is achieved differ markedly.

A novel ATAC-seq approach reveals lineage-specific reinforcement of the open chromatin landscape via cooperation between BAF and p63.

BACKGROUND: Open chromatin regions are correlated with active regulatory elements in development and are dysregulated in diseases. The BAF (SWI/SNF) complex is essential for development, and has been demonstrated to remodel reconstituted chromatin in vitro and to control the accessibility of a few individual regions in vivo. However, it remains unclear where and how BAF controls the open chromatin landscape to regulate developmental processes, such as human epidermal differentiation. RESULTS: Using a novel "on-plate" ATAC-sequencing approach for profiling open chromatin landscapes with a low number of adherent cells, we demonstrate that the BAF complex is essential for maintaining 11.6 % of open chromatin regions in epidermal differentiation. These BAF-dependent open chromatin regions are highly cell-type-specific and are strongly enriched for binding sites for p63, a master epidermal transcription factor. The DNA sequences of p63 binding sites intrinsically favor nucleosome formation and are inaccessible in other cell types without p63 to prevent ectopic activation. In epidermal cells, BAF and p63 mutually recruit each other to maintain 14,853 open chromatin regions. We further demonstrate that BAF and p63 cooperatively position nucleosomes away from p63 binding sites and recruit transcriptional machinery to control tissue differentiation. CONCLUSIONS: BAF displays high specificity in controlling the open chromatin landscape during epidermal differentiation by cooperating with the master transcription factor p63 to maintain lineage-specific open chromatin regions.

Individuality and variation of personal regulomes in primary human T cells.

Here we survey variation and dynamics of active regulatory elements genome-wide using longitudinal samples from human individuals. We applied Assay of Transposase Accessible Chromatin with sequencing (ATAC-seq) to map chromatin accessibility in primary CD4+ T cells isolated from standard blood draws of 12 healthy volunteers over time, from cancer patients, and during T cell activation. Over 4,000 predicted regulatory elements (7.2%) showed reproducible variation in accessibility between individuals. Gender was the most significant attributable source of variation. ATAC-seq revealed previously undescribed elements that escape X chromosome inactivation and predicted gender-specific gene regulatory networks across autosomes, which coordinately affect genes with immune function. Noisy regulatory elements with personal variation in accessibility are significantly enriched for autoimmune disease loci. Over one third of regulome variation lacked genetic variation in cis, suggesting contributions from environmental or epigenetic factors. These results refine concepts of human individuality and provide a foundational reference for comparing disease-associated regulomes.

STAT3 acts through pre-existing nucleosome-depleted regions bound by FOS during an epigenetic switch linking inflammation to cancer.

BACKGROUND: Transient induction of the Src oncoprotein in a non-transformed breast cell line can initiate an epigenetic switch to a cancer cell via a positive feedback loop that involves activation of the signal transducer and activator of transcription 3 protein (STAT3) and NF-κB transcription factors. RESULTS: We show that during the transformation process, nucleosome-depleted regions (defined by formaldehyde-assisted isolation of regulatory elements (FAIRE)) are largely unchanged and that both before and during transformation, STAT3 binds almost exclusively to previously open chromatin regions. Roughly, a third of the transformation-inducible genes require STAT3 for the induction. STAT3 and NF-κB appear to drive the regulation of different gene sets during the transformation process. Interestingly, STAT3 directly regulates the expression of NFKB1, which encodes a subunit of NF-κB, and IL6, a cytokine that stimulates STAT3 activity. Lastly, many STAT3 binding sites are also bound by FOS and the expression of several AP-1 factors is altered during transformation in a STAT3-dependent manner, suggesting that STAT3 may cooperate with AP-1 proteins. CONCLUSIONS: These observations uncover additional complexities to the inflammatory feedback loop that are likely to contribute to the epigenetic switch. In addition, gene expression changes during transformation, whether driven by pre-existing or induced transcription factors, occur largely through pre-existing nucleosome-depleted regions.

A detailed protocol for formaldehyde-assisted isolation of regulatory elements (FAIRE).

Nucleosome displacement is a key event in the regulation of gene expression in the eukaryotic genome. This unit details an approach called Formaldehyde-Assisted Isolation of Regulatory Elements (FAIRE) for isolating nucleosome-depleted regions. FAIRE does not rely on the use of antibodies or enzymes, and has proven successful in most eukaryotic cells and tissues. The set of regulatory elements enriched by FAIRE is similar to those identified through DNase hypersensitivity. The enriched fragments can be detected by quantitative PCR, tiling DNA microarrays, or next-generation sequencing. Although the signal-to-noise ratio is typically lower than that observed for DNase assays, FAIRE has high sample-to-sample reproducibility, requires very low amounts of input material, is inexpensive, is amenable to high-throughput adaptations, and is a relatively simple procedure with a high rate of success, even for those without extensive experience in molecular biology protocols.

Using formaldehyde-assisted isolation of regulatory elements (FAIRE) to isolate active regulatory DNA.

Eviction or destabilization of nucleosomes from chromatin is a hallmark of functional regulatory elements in eukaryotic genomes. Historically identified by nuclease hypersensitivity, these regulatory elements are typically bound by transcription factors or other regulatory proteins. FAIRE (formaldehyde-assisted isolation of regulatory elements) is an alternative approach to identify these genomic regions and has proven successful in a multitude of eukaryotic cell and tissue types. Cells or dissociated tissues are cross-linked briefly with formaldehyde, lysed and sonicated. Sheared chromatin is subjected to phenol/chloroform extraction and the isolated DNA, typically encompassing 1-3% of the human genome, is purified. We provide guidelines for quantitative analysis by PCR, microarrays or next-generation sequencing. Regulatory elements enriched by FAIRE have high concordance with those identified by nuclease hypersensitivity or chromatin immunoprecipitation (ChIP), and the entire procedure can be completed in 3 d. FAIRE has low technical variability, which allows its usage in large-scale studies of chromatin from normal or diseased tissues.

ZINBA integrates local covariates with DNA-seq data to identify broad and narrow regions of enrichment, even within amplified genomic regions.

ZINBA (Zero-Inflated Negative Binomial Algorithm) identifies genomic regions enriched in a variety of ChIP-seq and related next-generation sequencing experiments (DNA-seq), calling both broad and narrow modes of enrichment across a range of signal-to-noise ratios. ZINBA models and accounts for factors that co-vary with background or experimental signal, such as G/C content, and identifies enrichment in genomes with complex local copy number variations. ZINBA provides a single unified framework for analyzing DNA-seq experiments in challenging genomic contexts.