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

Decoding transcriptional enhancers: Evolving from annotation to functional interpretation

Deciphering the intricate molecular processes that orchestrate the spatial and temporal regulation of genes has become an increasingly major focus of biological research. The differential expression of genes by diverse cell types with a common genome is a hallmark of complex cellular functions, as well as the basis for multicellular life. Importantly, a more coherent understanding of gene regulation is critical for defining developmental processes, evolutionary principles and disease etiologies. Here we present our current understanding of gene regulation by focusing on the role of enhancer elements in these complex processes. Although functional genomic methods have provided considerable advances to our understanding of gene regulation, these assays, which are usually performed on a genome-wide scale, typically provide correlative observations that lack functional interpretation. Recent innovations in genome editing technologies have placed gene regulatory studies at an exciting crossroads, as systematic, functional evaluation of enhancers and other transcriptional regulatory elements can now be performed in a coordinated, high-throughput manner across the entire genome. This review provides insights on transcriptional enhancer function, their role in development and disease, and catalogues experimental tools commonly used to study these elements. Additionally, we discuss the crucial role of novel techniques in deciphering the complex gene regulatory landscape and how these studies will shape future research.

Simultaneous Pathway Activity Inference and Gene Expression Analysis Using RNA Sequencing

Reporter gene assays are a venerable tool for studying signaling pathways, but they lack the throughput and complexity necessary to contribute to a systems-level understanding of endogenous signaling networks. We present a parallel reporter assay, transcription factor activity sequencing (TF-seq), built on synthetic DNA enhancer elements, which enables parallel measurements in primary cells of the transcriptome and transcription factor activity from more than 40 signaling pathways. Using TF-seq in Myd88(-/-) macrophages, we captured dynamic pathway activity changes underpinning the global transcriptional changes of the innate immune response. We also applied TF-seq to investigate small molecule mechanisms of action and find a role for NF-κB activation and coordination of the STAT1 response in the macrophage reaction to the anti-inflammatory natural product halofuginone. Simultaneous TF-seq and global gene expression profiling represent an integrative approach for gaining mechanistic insight into pathway activity and transcriptional changes that result from genetic and small molecule perturbations.

The Dynamic Regulatory Genome of Capsaspora and the Origin of Animal Multicellularity

The unicellular ancestor of animals had a complex repertoire of genes linked to multicellular processes. This suggests that changes in the regulatory genome, rather than in gene innovation, were key to the origin of animals. Here, we carry out multiple functional genomic assays in Capsaspora owczarzaki, the unicellular relative of animals with the largest known gene repertoire for transcriptional regulation. We show that changing chromatin states, differential lincRNA expression, and dynamic cis-regulatory sites are associated with life cycle transitions in Capsaspora. Moreover, we demonstrate conservation of animal developmental transcription-factor networks and extensive network interconnection in this premetazoan organism. In contrast, however, Capsaspora lacks animal promoter types, and its regulatory sites are small, proximal, and lack signatures of animal enhancers. Overall, our results indicate that the emergence of animal multicellularity was linked to a major shift in genome cis-regulatory complexity, most notably the appearance of distal enhancer regulation.

Developmental Acquisition of Regulomes Underlies Innate Lymphoid Cell Functionality

Innate lymphoid cells (ILCs) play key roles in host defense, barrier integrity, and homeostasis and mirror adaptive CD4(+) T helper (Th) cell subtypes in both usage of effector molecules and transcription factors. To better understand the relationship between ILC subsets and their Th cell counterparts, we measured genome-wide chromatin accessibility. We find that chromatin in proximity to effector genes is selectively accessible in ILCs prior to high-level transcription upon activation. Accessibility of these regions is acquired in a stepwise manner during development and changes little after in vitro or in vivo activation. Conversely, dramatic chromatin remodeling occurs in naive CD4(+) T cells during Th cell differentiation using a type-2-infection model. This alteration results in a substantial convergence of Th2 cells toward ILC2 regulomes. Our data indicate extensive sharing of regulatory circuitry across the innate and adaptive compartments of the immune system, in spite of their divergent developing pathways.

Medial HOXA genes demarcate haematopoietic stem cell fate during human development

Pluripotent stem cells (PSC) may provide a potential source of haematopoietic stem/progenitor cells (HSPCs) for transplantation; however, unknown molecular barriers prevent the self-renewal of PSC-HSPCs. Using two-step differentiation, human embryonic stem cells (hESCs) differentiated in vitro into multipotent haematopoietic cells that had CD34⁺CD38^−/loCD90⁺CD45⁺GPI-80⁺ foetal liver (FL) HSC immunophenotype, but displayed poor expansion potential and engraftment ability. Transcriptome analysis of immunophenotypic hESC-HSPCs revealed that, despite their molecular resemblance to FL-HSPCs, medial HOXA genes remained suppressed. Knockdown of HOXA7 disrupted FL-HSPC function and caused transcriptome dysregulation that resembled hESC-derived progenitors. Overexpression of medial HOXA genes prolonged FL-HSPC maintenance but was insufficient to confer self-renewal to hESC-HSPCs. Stimulation of retinoic acid signalling during endothelial-to-haematopoietic transition induced the HOXA cluster and other HSC/definitive haemogenic endothelium genes, and prolonged HSPC maintenance in culture. Thus, retinoic acid signalling-induced medial HOXA gene expression marks the establishment of the definitive HSC fate and controls HSC identity and function.

Age-Dependent Pancreatic Gene Regulation Reveals Mechanisms Governing Human β Cell Function

Intensive efforts are focused on identifying regulators of human pancreatic islet cell growth and maturation to accelerate development of therapies for diabetes. After birth, islet cell growth and function are dynamically regulated; however, establishing these age-dependent changes in humans has been challenging. Here, we describe a multimodal strategy for isolating pancreatic endocrine and exocrine cells from children and adults to identify age-dependent gene expression and chromatin changes on a genomic scale. These profiles revealed distinct proliferative and functional states of islet α cells or β cells and histone modifications underlying age-dependent gene expression changes. Expression of SIX2 and SIX3, transcription factors without prior known functions in the pancreas and linked to fasting hyperglycemia risk, increased with age specifically in human islet β cells. SIX2 and SIX3 were sufficient to enhance insulin content or secretion in immature β cells. Our work provides a unique resource to study human-specific regulators of islet cell maturation and function.

MNase titration reveals differences between nucleosome occupancy and chromatin accessibility

Chromatin accessibility plays a fundamental role in gene regulation. Nucleosome placement, usually measured by quantifying protection of DNA from enzymatic digestion, can regulate accessibility. We introduce a metric that uses micrococcal nuclease (MNase) digestion in a novel manner to measure chromatin accessibility by combining information from several digests of increasing depths. This metric, MACC (MNase accessibility), quantifies the inherent heterogeneity of nucleosome accessibility in which some nucleosomes are seen preferentially at high MNase and some at low MNase. MACC interrogates each genomic locus, measuring both nucleosome location and accessibility in the same assay. MACC can be performed either with or without a histone immunoprecipitation step, and thereby compares histone and non-histone protection. We find that changes in accessibility at enhancers, promoters and other regulatory regions do not correlate with changes in nucleosome occupancy. Moreover, high nucleosome occupancy does not necessarily preclude high accessibility, which reveals novel principles of chromatin regulation.

Cell freezing protocol suitable for ATAC-Seq on motor neurons derived from human induced pluripotent stem cells

In recent years, the assay for transposase-accessible chromatin using sequencing (ATAC-Seq) has become a fundamental tool of epigenomic research. However, it is difficult to perform this technique on frozen samples because freezing cells before extracting nuclei can impair nuclear integrity and alter chromatin structure, especially in fragile cells such as neurons. Our aim was to develop a protocol for freezing neuronal cells that is compatible with ATAC-Seq; we focused on a disease-relevant cell type, namely motor neurons differentiated from induced pluripotent stem cells (iMNs) from a patient affected by spinal muscular atrophy. We found that while flash-frozen iMNs are not suitable for ATAC-Seq, the assay is successful with slow-cooled cryopreserved cells. Using this method, we were able to isolate high quality, intact nuclei, and we verified that epigenetic results from fresh and cryopreserved iMNs quantitatively agree.

Distinct Gene Regulatory Pathways for Human Innate versus Adaptive Lymphoid Cells

Innate lymphoid cells (ILCs) serve as sentinels in mucosal tissues, sensing release of soluble inflammatory mediators, rapidly communicating danger via cytokine secretion, and functioning as guardians of tissue homeostasis. Although ILCs have been extensively studied in model organisms, little is known about these "first responders" in humans, especially their lineage and functional kinships to cytokine-secreting T helper (Th) cell counterparts. Here, we report gene regulatory circuitries for four human ILC-Th counterparts derived from mucosal environments, revealing that each ILC subset diverges as a distinct lineage from Th and circulating natural killer cells but shares circuitry devoted to functional polarization with their Th counterparts. Super-enhancers demarcate cohorts of cell-identity genes in each lineage, uncovering new modes of regulation for signature cytokines, new molecules that likely impart important functions to ILCs, and potential mechanisms for autoimmune disease SNP associations within ILC-Th subsets.

Identification of plant transcription factor target sequences

Regulation of gene expression depends on specific cis-regulatory sequences located in the gene promoter regions. These DNA sequences are recognized by transcription factors (TFs) in a sequence-specific manner, and their identification could help to elucidate the regulatory networks that underlie plant physiological responses to developmental programs or to environmental adaptation. Here we review recent advances in high throughput methodologies for the identification of plant TF binding sites. Several approaches offer a map of the TF binding locations in vivo and of the dynamics of the gene regulatory networks. As an alternative, high throughput in vitro methods provide comprehensive determination of the DNA sequences recognized by TFs. These advances are helping to decipher the regulatory lexicon and to elucidate transcriptional network hierarchies in plants in response to internal or external cues. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

OVOL2 Maintains the Transcriptional Program of Human Corneal Epithelium by Suppressing Epithelial-to-Mesenchymal Transition

In development, embryonic ectoderm differentiates into neuroectoderm and surface ectoderm using poorly understood mechanisms. Here, we show that the transcription factor OVOL2 maintains the transcriptional program of human corneal epithelium cells (CECs), a derivative of the surface ectoderm, and that OVOL2 may regulate the differential transcriptional programs of the two lineages. A functional screen identified OVOL2 as a repressor of mesenchymal genes to maintain CECs. Transduction of OVOL2 with several other transcription factors induced the transcriptional program of CECs in fibroblasts. Moreover, neuroectoderm derivatives were found to express mesenchymal genes, and OVOL2 alone could induce the transcriptional program of CECs in neural progenitors by repressing these genes while activating epithelial genes. Our data suggest that the difference between the transcriptional programs of some neuroectoderm- and surface ectoderm-derivative cells may be regulated in part by a reciprocally repressive mechanism between epithelial and mesenchymal genes, as seen in epithelial-to-mesenchymal transition.

Network-based approaches for analysis of complex biological systems

Cells function and respond to changes in their environment by the coordinated activity of their molecular components, including mRNAs, proteins and metabolites. At the heart of proper cellular function are molecular networks connecting these components to process extra-cellular environmental signals and drive dynamic, context-specific cellular responses. Network-based computational approaches aim to systematically integrate measurements from high-throughput experiments to gain a global understanding of cellular function under changing environmental conditions. We provide an overview of recent methodological developments toward solving two major computational problems within this field in the past two years (2013-2015): network reconstruction and network-based interpretation. Looking forward, we envision development of methods that can predict phenotypes with high accuracy as well as provide biologically plausible mechanistic hypotheses.

Evolution of Gene Regulation in Humans

As a species, we possess unique biological features that distinguish us from other primates. Here, we review recent efforts to identify changes in gene regulation that drove the evolution of novel human phenotypes. We discuss genotype-directed comparisons of human and nonhuman primate genomes to identify human-specific genetic changes that may encode new regulatory functions. We also review phenotype-directed approaches, which use comparisons of gene expression or regulatory function in homologous human and nonhuman primate cells and tissues to identify changes in expression levels or regulatory activity that may be due to genetic changes in humans. Together, these studies are beginning to reveal the landscape of regulatory innovation in human evolution and point to specific regulatory changes for further study. Finally, we highlight two novel strategies to model human-specific regulatory functions in vivo: primate induced pluripotent stem cells and the generation of humanized mice by genome editing.

Impact of the X Chromosome and sex on regulatory variation

The X Chromosome, with its unique mode of inheritance, contributes to differences between the sexes at a molecular level, including sex-specific gene expression and sex-specific impact of genetic variation. Improving our understanding of these differences offers to elucidate the molecular mechanisms underlying sex-specific traits and diseases. However, to date, most studies have either ignored the X Chromosome or had insufficient power to test for the sex-specific impact of genetic variation. By analyzing whole blood transcriptomes of 922 individuals, we have conducted the first large-scale, genome-wide analysis of the impact of both sex and genetic variation on patterns of gene expression, including comparison between the X Chromosome and autosomes. We identified a depletion of expression quantitative trait loci (eQTL) on the X Chromosome, especially among genes under high selective constraint. In contrast, we discovered an enrichment of sex-specific regulatory variants on the X Chromosome. To resolve the molecular mechanisms underlying such effects, we generated chromatin accessibility data through ATAC-sequencing to connect sex-specific chromatin accessibility to sex-specific patterns of expression and regulatory variation. As sex-specific regulatory variants discovered in our study can inform sex differences in heritable disease prevalence, we integrated our data with genome-wide association study data for multiple immune traits identifying several traits with significant sex biases in genetic susceptibilities. Together, our study provides genome-wide insight into how genetic variation, the X Chromosome, and sex shape human gene regulation and disease.

Modeling the causal regulatory network by integrating chromatin accessibility and transcriptome data

Cell packs a lot of genetic and regulatory information through a structure known as chromatin, i.e. DNA is wrapped around histone proteins and is tightly packed in a remarkable way. To express a gene in a specific coding region, the chromatin would open up and DNA loop may be formed by interacting enhancers and promoters. Furthermore, the mediator and cohesion complexes, sequence-specific transcription factors, and RNA polymerase II are recruited and work together to elaborately regulate the expression level. It is in pressing need to understand how the information, about when, where, and to what degree genes should be expressed, is embedded into chromatin structure and gene regulatory elements. Thanks to large consortia such as Encyclopedia of DNA Elements (ENCODE) and Roadmap Epigenomic projects, extensive data on chromatin accessibility and transcript abundance are available across many tissues and cell types. This rich data offer an exciting opportunity to model the causal regulatory relationship. Here, we will review the current experimental approaches, foundational data, computational problems, interpretive frameworks, and integrative models that will enable the accurate interpretation of regulatory landscape. Particularly, we will discuss the efforts to organize, analyze, model, and integrate the DNA accessibility data, transcriptional data, and functional genomic regions together. We believe that these efforts will eventually help us understand the information flow within the cell and will influence research directions across many fields.

Single-cell epigenomics: powerful new methods for understanding gene regulation and cell identity

Emerging single-cell epigenomic methods are being developed with the exciting potential to transform our knowledge of gene regulation. Here we review available techniques and future possibilities, arguing that the full potential of single-cell epigenetic studies will be realized through parallel profiling of genomic, transcriptional, and epigenetic information.

Cell type-specific epigenome profiling using affinity-purified nuclei

The development of a multicellular organism from a single zygote depends on the differentiation of progenitor cells to specialized cell types. The differentiation of these cell types is associated with changes in gene expression and the underlying chromatin landscape. To understand how these processes are regulated, it is desirable to understand how the chromatin features that constitute the epigenome differ between cell types at any given time during development. INTACT, a method for the cell type-specific purification of nuclei that can be used for the isolation of both RNA and chromatin, has emerged as a powerful tool to simultaneously study gene expression and chromatin profiles specifically in cell types of interest. In this review, we focus on the application of INTACT to different model organisms and discuss its potential for profiling cell types in their developmental context.

Strategies for targeting primate neural circuits with viral vectors

Understanding how the brain works requires understanding how different types of neurons contribute to circuit function and organism behavior. Progress on this front has been accelerated by optogenetics and chemogenetics, which provide an unprecedented level of control over distinct neuronal types in small animals. In primates, however, targeting specific types of neurons with these tools remains challenging. In this review, we discuss existing and emerging strategies for directing genetic manipulations to targeted neurons in the adult primate central nervous system. We review the literature on viral vectors for gene delivery to neurons, focusing on adeno-associated viral vectors and lentiviral vectors, their tropism for different cell types, and prospects for new variants with improved efficacy and selectivity. We discuss two projection targeting approaches for probing neural circuits: anterograde projection targeting and retrograde transport of viral vectors. We conclude with an analysis of cell type-specific promoters and other nucleotide sequences that can be used in viral vectors to target neuronal types at the transcriptional level.

The mechanisms shaping the single-cell transcriptional landscape

Recent technological and computational advances in understanding the transcriptional and chromatin features of single cells have begun answering longstanding questions in the extent and impact of biological heterogeneity. Here, we outline the intrinsic and extrinsic mechanisms that underlie the transcriptional and functional diversity within superficially homogeneous populations, and we discuss how fascinating new studies have afforded novel insight into each mechanism. The studies are chosen in part to include initial reports of novel functional genomics tools where the eventual applications will clearly have profound impact on our understanding the dynamics of cell-to-cell transcriptional variation-from individual cells to whole organisms.

Foxc1 reinforces quiescence in self-renewing hair follicle stem cells

Stem cell quiescence preserves the cell reservoir by minimizing cell division over extended periods of time. Self-renewal of quiescent stem cells (SCs) requires the reentry into the cell cycle. In this study, we show that murine hair follicle SCs induce the Foxc1 transcription factor when activated. Deleting Foxc1 in activated, but not quiescent, SCs causes failure of the cells to reestablish quiescence and allows premature activation. Deleting Foxc1 in the SC niche of gene-targeted mice leads to loss of the old hair without impairing quiescence. In self-renewing SCs, Foxc1 activates Nfatc1 and bone morphogenetic protein (BMP) signaling, two key mechanisms that govern quiescence. These findings reveal a dynamic, cell-intrinsic mechanism used by hair follicle SCs to reinforce quiescence upon self-renewal and suggest a unique ability of SCs to maintain cell identity.

Principles of metadata organization at the ENCODE data coordination center

The Encyclopedia of DNA Elements (ENCODE) Data Coordinating Center (DCC) is responsible for organizing, describing and providing access to the diverse data generated by the ENCODE project. The description of these data, known as metadata, includes the biological sample used as input, the protocols and assays performed on these samples, the data files generated from the results and the computational methods used to analyze the data. Here, we outline the principles and philosophy used to define the ENCODE metadata in order to create a metadata standard that can be applied to diverse assays and multiple genomic projects. In addition, we present how the data are validated and used by the ENCODE DCC in creating the ENCODE Portal (https://www.encodeproject.org/).

Database URL:www.encodeproject.org

Assay for transposase-accessible chromatin and circularized chromosome conformation capture, two methods to explore the regulatory landscapes of genes in zebrafish

Accurate transcriptional control of genes is fundamental for the correct functioning of organs and developmental processes. This control depends on the interplay between the promoter of genes and other noncoding sequences, whose interaction is mediated by 3D chromatin arrangements. Thus, the detailed description of transcriptional regulatory landscapes is essential to understand the mechanisms of transcriptional regulation. However, to achieve that, two important challenges have to be faced: (1) the identification of the noncoding sequences that contribute to gene transcription and (2) the association of these sequences to the respective genes they control. In this chapter, we describe two protocols that allow overcoming these important challenges: the assay for transposase-accessible chromatin using sequencing (ATAC-seq) and circularized chromosome conformation capture (4C-seq). ATAC-seq is a very efficient technique that, using a very low number of cells as starting material, allows the identification of active chromatin regions genome wide, whereas 4C-seq detects the subset of sequences that interact specifically with the promoter of a given gene. When combined, both techniques provide a comprehensive snapshot of the regulatory landscapes of developmental genes. The protocols we present here have been optimized for teleost fish samples, zebrafish and medaka, allowing the in-depth study of transcriptional regulation in these two emerging animal models. Given the amenability and easy genetic manipulation of these two experimental systems, we anticipate that they will be important in revealing general principles of the vertebrate regulatory genome.

C/EBPα creates elite cells for iPSC reprogramming by upregulating Klf4 and increasing the levels of Lsd1 and Brd4

Reprogramming somatic cells into induced pluripotent stem cells (iPSCs) is typically inefficient and has been explained by elite-cell and stochastic models. We recently reported that B cells exposed to a pulse of C/EBPα (Bα' cells) behave as elite cells, in that they can be rapidly and efficiently reprogrammed into iPSCs by the Yamanaka factors OSKM. Here we show that C/EBPα post-transcriptionally increases the abundance of several hundred proteins, including Lsd1, Hdac1, Brd4, Med1 and Cdk9, components of chromatin-modifying complexes present at super-enhancers. Lsd1 was found to be required for B cell gene silencing and Brd4 for the activation of the pluripotency program. C/EBPα also promotes chromatin accessibility in pluripotent cells and upregulates Klf4 by binding to two haematopoietic enhancers. Bα' cells share many properties with granulocyte/macrophage progenitors, naturally occurring elite cells that are obligate targets for leukaemic transformation, whose formation strictly requires C/EBPα.

Understanding Celiac Disease by Genomics

Celiac disease (CeD) is a complex immune-mediated disease. Genetic studies have implicated 43 predisposing loci that collectively explain some 50% of the genetic variance in CeD. More than ∼90% of CeD-associated single nucleotide polymorphisms (SNPs) localize to the non-coding genome, which we need to better understand to translate genetic knowledge into clinical practice. New genomic technologies and resources are permitting a systematic analysis of the functional elements in the non-coding part of the genome. Here we explain how investigating the regulatory and epigenomic landscape will help to pinpoint the cell types involved in CeD, and the driver genes and gene regulatory networks that are affected by CeD-associated SNPs.

Epigenomic landscapes of retinal rods and cones

Rod and cone photoreceptors are highly similar in many respects but they have important functional and molecular differences. Here, we investigate genome-wide patterns of DNA methylation and chromatin accessibility in mouse rods and cones and correlate differences in these features with gene expression, histone marks, transcription factor binding, and DNA sequence motifs. Loss of NR2E3 in rods shifts their epigenomes to a more cone-like state. The data further reveal wide differences in DNA methylation between retinal photoreceptors and brain neurons. Surprisingly, we also find a substantial fraction of DNA hypo-methylated regions in adult rods that are not in active chromatin. Many of these regions exhibit hallmarks of regulatory regions that were active earlier in neuronal development, suggesting that these regions could remain undermethylated due to the highly compact chromatin in mature rods. This work defines the epigenomic landscapes of rods and cones, revealing features relevant to photoreceptor development and function.

DOI:http://dx.doi.org/10.7554/eLife.11613.001

Vision in humans is made possible by a light-sensing sheet of cells at the back of the eye called the retina. The surface of the retina is populated by specialized sensory cells, known as rods and cones. The rod cells detect very dim light, while the cones are less sensitive to light but are used to detect color. Together, the rods and cones gather the information needed to create a picture that is then transmitted to the brain.

Rods and cones have been studied for decades, and genetic analyses have revealed the patterns of gene expression that lead a cell to develop into either a rod or a cone. Researchers have also identified several key regulatory genes that control these patterns, but less is known about the role of other factors that control the expression of genes.

Chemical modifications to DNA or modifications to the proteins associated with DNA – which are collectively called epigenetic modifications – can either promote or inhibit the activation of nearby genes. Now, Mo et al. have shown that rods and cones from mice have very different patterns of epigenetic modifications. The experiments also revealed that many sections of DNA that are marked to promote gene activation contain known rod-specific or cone-specific genes; and that rod cells need a known regulatory gene to develop their specific pattern of epigenetic modifications. Finally, Mo et al. showed that epigenetic regulation differed between brain cells and rods and cones.

These insights into epigenetic regulation of rod and cone genes may help explain why some people with eye diseases caused by the same genetic mutation may develop symptoms at different ages or lose vision at different rates. The new information about gene regulation may also help scientists to reprogram stem cells to become healthy rods or cones that could be transplanted into people with eye disease to restore their vision.

DOI:http://dx.doi.org/10.7554/eLife.11613.002

Smash and DASH with Cas9

DeRisi and colleagues present a creative application for Cas9 in vitro, using it to deplete unwanted sequence from DNA libraries. It seems plausible that the in vitro use of CRISPR/Cas9 has unrealized potential to revolutionize the practice of molecular biology well beyond genome editing.

Tunicates: exploring the sea shores and roaming the open ocean. A tribute to Thomas Huxley

This review is a tribute to the remarkable contributions of Thomas Huxley to the biology of tunicates, the likely sister group of vertebrates. In 1851, the great biologist and philosopher published two landmark papers on pelagic tunicates in the Philosophical Transactions of the Royal Society. They were dedicated to the description of the adult anatomy and life cycle of thaliaceans and appendicularians, the pelagic relatives of ascidians. In the first part of this review, we discuss the novel anatomical observations and evolutionary hypotheses made by Huxley, which would have a lasting influence on tunicate biology. We also briefly comment on the more philosophical reflections of Huxley on individuality. In the second part, we stress the originality and relevance of past and future studies of tunicates in the resolution of major biological issues. In particular, we focus on the complex relationship between genotype and phenotype and the phenomenon of developmental system drift. We propose that more than 150 years after Huxley's papers, tunicate embryos are still worth studying in their own right, independently of their evolutionary proximity to vertebrates, as they provide original and crucial insights into the process of animal evolution. Tunicates are still at the forefront of biological research.

4C-ker: A Method to Reproducibly Identify Genome-Wide Interactions Captured by 4C-Seq Experiments

4C-Seq has proven to be a powerful technique to identify genome-wide interactions with a single locus of interest (or "bait") that can be important for gene regulation. However, analysis of 4C-Seq data is complicated by the many biases inherent to the technique. An important consideration when dealing with 4C-Seq data is the differences in resolution of signal across the genome that result from differences in 3D distance separation from the bait. This leads to the highest signal in the region immediately surrounding the bait and increasingly lower signals in far-cis and trans. Another important aspect of 4C-Seq experiments is the resolution, which is greatly influenced by the choice of restriction enzyme and the frequency at which it can cut the genome. Thus, it is important that a 4C-Seq analysis method is flexible enough to analyze data generated using different enzymes and to identify interactions across the entire genome. Current methods for 4C-Seq analysis only identify interactions in regions near the bait or in regions located in far-cis and trans, but no method comprehensively analyzes 4C signals of different length scales. In addition, some methods also fail in experiments where chromatin fragments are generated using frequent cutter restriction enzymes. Here, we describe 4C-ker, a Hidden-Markov Model based pipeline that identifies regions throughout the genome that interact with the 4C bait locus. In addition, we incorporate methods for the identification of differential interactions in multiple 4C-seq datasets collected from different genotypes or experimental conditions. Adaptive window sizes are used to correct for differences in signal coverage in near-bait regions, far-cis and trans chromosomes. Using several datasets, we demonstrate that 4C-ker outperforms all existing 4C-Seq pipelines in its ability to reproducibly identify interaction domains at all genomic ranges with different resolution enzymes.

Multifunctional role of the transcription factor Blimp-1 in coordinating plasma cell differentiation

The transcription factor Blimp-1 is necessary for the generation of plasma cells. Here we studied its functions in plasmablast differentiation by identifying regulated Blimp-1 target genes. Blimp-1 promoted the migration and adhesion of plasmablasts. It directly repressed genes encoding several transcription factors and Aicda (which encodes the cytidine deaminase AID) and thus silenced B cell-specific gene expression, antigen presentation and class-switch recombination in plasmablasts. It directly activated genes, which led to increased expression of the plasma cell regulator IRF4 and proteins involved in immunoglobulin secretion. Blimp-1 induced the transcription of immunoglobulin genes by controlling the 3' enhancers of the loci encoding the immunoglobulin heavy chain (Igh) and κ-light chain (Igk) and, furthermore, regulated the post-transcriptional expression switch from the membrane-bound form of the immunoglobulin heavy chain to its secreted form by activating Ell2 (which encodes the transcription-elongation factor ELL2). Notably, Blimp-1 recruited chromatin-remodeling and histone-modifying complexes to regulate its target genes. Hence, many essential functions of plasma cells are under the control of Blimp-1.

Active DNA demethylation at enhancers during the vertebrate phylotypic period

The vertebrate body plan and organs are shaped during a conserved embryonic phase called the phylotypic stage. However, the mechanisms that guide the epigenome through this transition and their evolutionary conservation remain elusive. Here we report widespread DNA demethylation of enhancers during the phylotypic period in zebrafish, Xenopus tropicalis and mouse. These enhancers are linked to developmental genes that display coordinated transcriptional and epigenomic changes in the diverse vertebrates during embryogenesis. Binding of Tet proteins to (hydroxy)methylated DNA and enrichment of 5-hydroxymethylcytosine in these regions implicated active DNA demethylation in this process. Furthermore, loss of function of Tet1, Tet2 and Tet3 in zebrafish reduced chromatin accessibility and increased methylation levels specifically at these enhancers, indicative of DNA methylation being an upstream regulator of phylotypic enhancer function. Overall, our study highlights a regulatory module associated with the most conserved phase of vertebrate embryogenesis and suggests an ancient developmental role for Tet dioxygenases.

GATA3-dependent cellular reprogramming requires activation-domain dependent recruitment of a chromatin remodeler

Background

Transcription factor-dependent cellular reprogramming is integral to normal development and is central to production of induced pluripotent stem cells. This process typically requires pioneer transcription factors (TFs) to induce de novo formation of enhancers at previously closed chromatin. Mechanistic information on this process is currently sparse.

Results

Here we explore the mechanistic basis by which GATA3 functions as a pioneer TF in a cellular reprogramming event relevant to breast cancer, the mesenchymal to epithelial transition (MET). In some instances, GATA3 binds previously inaccessible chromatin, characterized by stable, positioned nucleosomes where it induces nucleosome eviction, alters local histone modifications, and remodels local chromatin architecture. At other loci, GATA3 binding induces nucleosome sliding without concomitant generation of accessible chromatin. Deletion of the transactivation domain retains the chromatin binding ability of GATA3 but cripples chromatin reprogramming ability, resulting in failure to induce MET.

Conclusions

These data provide mechanistic insights into GATA3-mediated chromatin reprogramming during MET, and suggest unexpected complexity to TF pioneering. Successful reprogramming requires stable binding to a nucleosomal site; activation domain-dependent recruitment of co-factors including BRG1, the ATPase subunit of the SWI/SNF chromatin remodeling complex; and appropriate genomic context. The resulting model provides a new conceptual framework for de novo enhancer establishment by a pioneer TF.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-0897-0) contains supplementary material, which is available to authorized users.

Ensembl regulation resources

New experimental techniques in epigenomics allow researchers to assay a diversity of highly dynamic features such as histone marks, DNA modifications or chromatin structure. The study of their fluctuations should provide insights into gene expression regulation, cell differentiation and disease. The Ensembl project collects and maintains the Ensembl regulation data resources on epigenetic marks, transcription factor binding and DNA methylation for human and mouse, as well as microarray probe mappings and annotations for a variety of chordate genomes. From this data, we produce a functional annotation of the regulatory elements along the human and mouse genomes with plans to expand to other species as data becomes available. Starting from well-studied cell lines, we will progressively expand our library of measurements to a greater variety of samples. Ensembl’s regulation resources provide a central and easy-to-query repository for reference epigenomes. As with all Ensembl data, it is freely available at http://www.ensembl.org, from the Perl and REST APIs and from the public Ensembl MySQL database server at ensembldb.ensembl.org.

Database URL: http://www.ensembl.org

Transcription factors LRF and BCL11A independently repress expression of fetal hemoglobin

Genes encoding human β-type globin undergo a developmental switch from embryonic to fetal to adult-type expression. Mutations in the adult form cause inherited hemoglobinopathies or globin disorders, including sickle cell disease and thalassemia. Some experimental results have suggested that these diseases could be treated by induction of fetal-type hemoglobin (HbF). However, the mechanisms that repress HbF in adults remain unclear. We found that the LRF/ZBTB7A transcription factor occupies fetal γ-globin genes and maintains the nucleosome density necessary for γ-globin gene silencing in adults, and that LRF confers its repressive activity through a NuRD repressor complex independent of the fetal globin repressor BCL11A. Our study may provide additional opportunities for therapeutic targeting in the treatment of hemoglobinopathies.

7SK-BAF axis controls pervasive transcription at enhancers

RNA functions at enhancers remain mysterious. Here we show that the 7SK small nuclear RNA (snRNA) inhibits enhancer transcription by modulating nucleosome position. 7SK occupies enhancers and super enhancers genome-wide in mouse and human cells, and 7SK is required to limit eRNA initiation and synthesis in a manner distinct from promoter pausing. Clustered elements at super enhancers uniquely require 7SK to prevent convergent transcription and DNA damage signaling. 7SK physically interacts with the BAF chromatin remodeling complex, recruit BAF to enhancers, and inhibits enhancer transcription by modulating chromatin structure. In turn, 7SK occupancy at enhancers coincides with Brd4 and is exquisitely sensitive to the bromodomain inhibitor JQ1. Thus, 7SK employs distinct mechanisms to counteract diverse consequences of pervasive transcription that distinguish super enhancers, enhancers, and promoters.

CD44+ Cells in Head and Neck Squamous Cell Carcinoma Suppress T-Cell-Mediated Immunity by Selective Constitutive and Inducible Expression of PD-L1

PURPOSE

Human tumors consist of heterogeneous populations of cells with distinct marker expression and functional properties. In squamous cell carcinoma of the head and neck (SCCHN), CD44 is a well-characterized marker of a resilient subpopulation of cells associated with increased tumorigenesis, radioresistance, and chemoresistance. Evidence indicates that these cells have an immunosuppressive phenotype; however, mechanisms have been elusive.

EXPERIMENTAL DESIGN

Using primary human SCCHN tumor samples and patient-derived xenografts, we examined the phenotypes of subsets of tumor cells and investigated mechanisms regulating their immunogenicity.

RESULTS

CD44(+) cells in primary human SCCHN were found to have an epithelial-to-mesenchymal (EMT) phenotype and were less immunogenic than CD44(-) cells when cultured with autologous CD8(+) tumor-infiltrating T cells. Selective expression of the programmed death-ligand 1 (PD-L1) was observed on CD44(+) cells compared with CD44(-) cells and was associated with constitutive phosphorylation of STAT3 on CD44(+) cells. Importantly, inhibition of STAT3 decreased expression of PD-L1 on CD44(+) cells. IFNγ treatment preferentially induced even further PD-L1 expression on CD44(+) cells and was associated with enhanced IFNγ receptor expression and phosphorylation of STAT1. Finally, the decreased immunogenicity of CD44(+) cells was partially reversed by antibody blockade of the programmed death 1 (PD-1) receptor, indicating that the differences in PD-L1 expression between CD44(+) and CD44(-) cells are biologically and clinically relevant.

CONCLUSIONS

Our findings provide a mechanism by which long-lived CD44(+) tumor-initiating cells can selectively evade host immune responses and provide rationale for targeting the PD-1 pathway in the adjuvant therapy setting of SCCHN. Clin Cancer Res; 22(14); 3571-81. ©2016 AACR.

Characterization of chromatin accessibility with a transposome hypersensitive sites sequencing (THS-seq) assay

Chromatin accessibility captures in vivo protein-chromosome binding status, and is considered an informative proxy for protein-DNA interactions. DNase I and Tn5 transposase assays require thousands to millions of fresh cells for comprehensive chromatin mapping. Applying Tn5 tagmentation to hundreds of cells results in sparse chromatin maps. We present a transposome hypersensitive sites sequencing assay for highly sensitive characterization of chromatin accessibility. Linear amplification of accessible DNA ends with in vitro transcription, coupled with an engineered Tn5 super-mutant, demonstrates improved sensitivity on limited input materials, and accessibility of small regions near distal enhancers, compared with ATAC-seq.

Which Genetics Variants in DNase-Seq Footprints Are More Likely to Alter Binding?

Large experimental efforts are characterizing the regulatory genome, yet we are still missing a systematic definition of functional and silent genetic variants in non-coding regions. Here, we integrated DNaseI footprinting data with sequence-based transcription factor (TF) motif models to predict the impact of a genetic variant on TF binding across 153 tissues and 1,372 TF motifs. Each annotation we derived is specific for a cell-type condition or assay and is locally motif-driven. We found 5.8 million genetic variants in footprints, 66% of which are predicted by our model to affect TF binding. Comprehensive examination using allele-specific hypersensitivity (ASH) reveals that only the latter group consistently shows evidence for ASH (3,217 SNPs at 20% FDR), suggesting that most (97%) genetic variants in footprinted regulatory regions are indeed silent. Combining this information with GWAS data reveals that our annotation helps in computationally fine-mapping 86 SNPs in GWAS hit regions with at least a 2-fold increase in the posterior odds of picking the causal SNP. The rich meta information provided by the tissue-specificity and the identity of the putative TF binding site being affected also helps in identifying the underlying mechanism supporting the association. As an example, the enrichment for LDL level-associated SNPs is 9.1-fold higher among SNPs predicted to affect HNF4 binding sites than in a background model already including tissue-specific annotation.

A single three-dimensional chromatin compartment in amphioxus indicates a stepwise evolution of vertebrate Hox bimodal regulation

The HoxA and HoxD gene clusters of jawed vertebrates are organized into bipartite three-dimensional chromatin structures that separate long-range regulatory inputs coming from the anterior and posterior Hox-neighboring regions. This architecture is instrumental in allowing vertebrate Hox genes to pattern disparate parts of the body, including limbs. Almost nothing is known about how these three-dimensional topologies originated. Here we perform extensive 4C-seq profiling of the Hox cluster in embryos of amphioxus, an invertebrate chordate. We find that, in contrast to the architecture in vertebrates, the amphioxus Hox cluster is organized into a single chromatin interaction domain that includes long-range contacts mostly from the anterior side, bringing distant cis-regulatory elements into contact with Hox genes. We infer that the vertebrate Hox bipartite regulatory system is an evolutionary novelty generated by combining ancient long-range regulatory contacts from DNA in the anterior Hox neighborhood with new regulatory inputs from the posterior side.

A zebrafish melanoma model reveals emergence of neural crest identity during melanoma initiation

The "cancerized field" concept posits that cancer-prone cells in a given tissue share an oncogenic mutation, but only discreet clones within the field initiate tumors. Most benign nevi carry oncogenic BRAF(V600E) mutations but rarely become melanoma. The zebrafish crestin gene is expressed embryonically in neural crest progenitors (NCPs) and specifically reexpressed in melanoma. Live imaging of transgenic zebrafish crestin reporters shows that within a cancerized field (BRAF(V600E)-mutant; p53-deficient), a single melanocyte reactivates the NCP state, revealing a fate change at melanoma initiation in this model. NCP transcription factors, including sox10, regulate crestin expression. Forced sox10 overexpression in melanocytes accelerated melanoma formation, which is consistent with activation of NCP genes and super-enhancers leading to melanoma. Our work highlights NCP state reemergence as a key event in melanoma initiation.

The Histone Variant MacroH2A1.2 Is Necessary for the Activation of Muscle Enhancers and Recruitment of the Transcription Factor Pbx1

Histone variants complement and integrate histone post-translational modifications in regulating transcription. The histone variant macroH2A1 (mH2A1) is almost three times the size of its canonical H2A counterpart, due to the presence of an ∼25 kDa evolutionarily conserved non-histone macro domain. Strikingly, mH2A1 can mediate both gene repression and activation. However, the molecular determinants conferring these alternative functions remain elusive. Here, we report that mH2A1.2 is required for the activation of the myogenic gene regulatory network and muscle cell differentiation. H3K27 acetylation at prospective enhancers is exquisitely sensitive to mH2A1.2, indicating a role of mH2A1.2 in imparting enhancer activation. Both H3K27 acetylation and recruitment of the transcription factor Pbx1 at prospective enhancers are regulated by mH2A1.2. Overall, our findings indicate a role of mH2A1.2 in marking regulatory regions for activation.

Genome-wide nucleosome specificity and function of chromatin remodellers in ES cells

ATP-dependent chromatin remodellers allow access to DNA for transcription factors and the general transcription machinery, but whether mammalian chromatin remodellers target specific nucleosomes to regulate transcription is unclear. Here we present genome-wide remodeller-nucleosome interaction profiles for the chromatin remodellers Chd1, Chd2, Chd4, Chd6, Chd8, Chd9, Brg1 and Ep400 in mouse embryonic stem (ES) cells. These remodellers bind one or both full nucleosomes that flank micrococcal nuclease (MNase)-defined nucleosome-free promoter regions (NFRs), where they separate divergent transcription. Surprisingly, large CpG-rich NFRs that extend downstream of annotated transcriptional start sites are nevertheless bound by non-nucleosomal or subnucleosomal histone variants (H3.3 and H2A.Z) and marked by H3K4me3 and H3K27ac modifications. RNA polymerase II therefore navigates hundreds of base pairs of altered chromatin in the sense direction before encountering an MNase-resistant nucleosome at the 3' end of the NFR. Transcriptome analysis after remodeller depletion reveals reciprocal mechanisms of transcriptional regulation by remodellers. Whereas at active genes individual remodellers have either positive or negative roles via altering nucleosome stability, at polycomb-enriched bivalent genes the same remodellers act in an opposite manner. These findings indicate that remodellers target specific nucleosomes at the edge of NFRs, where they regulate ES cell transcriptional programs.

Epigenetic States of nephron progenitors and epithelial differentiation

In mammals, formation of new nephrons ends perinatally due to consumption of mesenchymal progenitor cells. Premature depletion of progenitors due to prematurity or postnatal loss of nephrons due to injury causes chronic kidney disease and hypertension. Intensive efforts are currently invested in designing regenerative strategies to form new nephron progenitors from pluripotent cells, which upon further differentiation provide a potential source of new nephrons. To know if reprogramed renal cells can maintain their identity and fate requires knowledge of the epigenetic states of native nephron progenitors and their progeny. In this article, we summarize current knowledge and gaps in the epigenomic landscape of the developing kidney. We now know that Pax2/PTIP/H3K4 methyltransferase activity provides the initial epigenetic specification signal to the metanephric mesenchyme. During nephrogenesis, the cap mesenchyme housing nephron progenitors is enriched in bivalent chromatin marks; as tubulogenesis proceeds, the tubular epithelium acquires H3K79me2. The latter mark is uniquely induced during epithelial differentiation. Analysis of histone landscapes in clonal metanephric mesenchyme cell lines and in Wilms tumor and normal fetal kidney has revealed that promoters of poised nephrogenesis genes carry bivalent histone signatures in progenitors. Differentiation or stimulation of Wnt signaling promotes resolution of bivalency; this does not occur in Wilms tumor cells consistent with their developmental arrest. The use of small cell number ChIP-Seq should facilitate the characterization of the chromatin landscape of the metanephric mesenchyme and various nephron compartments during nephrogenesis. Only then we will know if stem and somatic cell reprogramming into kidney progenitors recapitulates normal development.

High-density genotyping of immune-related loci identifies new SLE risk variants in individuals with Asian ancestry

Systemic lupus erythematosus (SLE) has a strong but incompletely understood genetic architecture. We conducted an association study with replication in 4,478 SLE cases and 12,656 controls from six East Asian cohorts to identify new SLE susceptibility loci and better localize known loci. We identified ten new loci and confirmed 20 known loci with genome-wide significance. Among the new loci, the most significant locus was GTF2IRD1-GTF2I at 7q11.23 (rs73366469, Pmeta = 3.75 × 10(-117), odds ratio (OR) = 2.38), followed by DEF6, IL12B, TCF7, TERT, CD226, PCNXL3, RASGRP1, SYNGR1 and SIGLEC6. We identified the most likely functional variants at each locus by analyzing epigenetic marks and gene expression data. Ten candidate variants are known to alter gene expression in cis or in trans. Enrichment analysis highlights the importance of these loci in B cell and T cell biology. The new loci, together with previously known loci, increase the explained heritability of SLE to 24%. The new loci share functional and ontological characteristics with previously reported loci and are possible drug targets for SLE therapeutics.

Genomic Views of Transcriptional Enhancers: Essential Determinants of Cellular Identity and Activity-Dependent Responses in the CNS

Sprinkled throughout the genome are a million regulatory sequences called transcriptional enhancers that activate gene promoters in the right cells, at the right time. Enhancers endow the brain with its incredible diversity of cell types and also translate neural activity into gene induction. Thanks to rapid advances in genomic technologies, it is now possible to identify thousands of enhancers rapidly, test their transcriptional function en masse, and address their neurobiological functions via genome editing. Enhancers also promise to be a great technological opportunity for neuroscience, offering the potential for cell-type-specific genetic labeling and manipulation without the need for transgenesis. The objective of this review and the accompanying 2015 SfN mini-symposium is to highlight the use of new and emerging genomic technologies to probe enhancer function in the nervous system.

SIGNIFICANCE STATEMENT

Transcriptional enhancers turn on genes in the right cells, at the right time. Enhancers are also the genomic sequences that encode the incredible diversity of cell types in the brain and enable the brain to turn genes on in response to new experiences. New technology enables enhancers to be found and manipulated. The study of enhancers promises to inform our understanding of brain development and function. The application of enhancer technology holds promise in accelerating basic neuroscience research and enabling gene therapies to be targeted to specific cell types in the brain.

Cis-regulatory programs in the development and evolution of vertebrate paired appendages

Differential gene expression is the core of development, mediating the genetic changes necessary for determining cell identity. The regulation of gene activity by cis-acting elements (e.g., enhancers) is a crucial mechanism for determining differential gene activity by precise control of gene expression in embryonic space and time. Modifications to regulatory regions can have profound impacts on phenotype, and therefore developmental and evolutionary biologists have increasingly focused on elucidating the transcriptional control of genes that build and pattern body plans. Here, we trace the evolutionary history of transcriptional control of three loci key to vertebrate appendage development (Fgf8, Shh, and HoxD/A). Within and across these regulatory modules, we find both complex and flexible regulation in contrast with more fixed enhancers that appear unchanged over vast timescales of vertebrate evolution. The transcriptional control of vertebrate appendage development was likely already incredibly complex in the common ancestor of fish, implying that subtle changes to regulatory networks were more likely responsible for alterations in phenotype rather than the de novo addition of whole regulatory domains. Finally, we discuss the dangers of relying on inter-species transgenesis when testing enhancer function, and call for more controlled regulatory swap experiments when inferring the evolutionary history of enhancer elements.