Discovering enhancers directly by activity

DNA sequences that regulate the timing, tissue-specificity, and level of gene expression are critical determinants of normal organismal development and differentiation1. Changes in their function likely underly biological processes ranging from disease susceptibility2-4 to speciation5. However, our ignorance of a grammar by which to identify and decipher such regulatory sequences has impeded discovery and manipulation of these important genetic elements. Two papers in this issue of Nature Methods report novel methods for discovering regulatory sequences by directly monitoring their activity during the initial screening assay: functional identification of regulatory elements within active chromatin (FIREWACh6) and site-specific integration fluoresence activated cell sorting followed by sequencing (SIF-seq7). These methods open new avenues for discovery of regulatory sequences. The critical roles of regulatory sequences fostered decades of research into their structures and mechanisms of action. Most regulatory regions are modular, comprised of multiple binding sites for transcription factors (TFs). The TF binding site motifs direct binding by the TFs, but such short (frequently 6-8 bp) sequences do not provide sufficient discriminatory information to explain specific TF binding genome-wide. Regulatory regions control genes on the same chromosome (in cis), but they can be quite far away and need not target the nearest gene. The limited information in motifs and large distances over which enhancers act have impeded the formulation of a regulatory “grammar”, i.e. rules for interpreting regulatory information in DNA sequences of complex organisms. Thus discovery of such cis-regulatory modules (CRMs) based solely on DNA sequence information is daunting8. Phylogenetic conservation or epigenetic features such as chromatin accessibility, histone modifications, and binding by TFs and co-activators can be used to improve CRM prediction. However, many of the DNA segments marked by these features are not independently active in enhancer assays8. These apparent “false positives” can be identified by experimental tests, and several groups have introduced massively parallel reporter assays (MPRAs) for high throughput testinge.g.,9. The new methods reported in this issue take a different approach to identifying CRMs, using their ability to boost gene expression as part of the assay to identify the candidates. Thus they are not biased toward evolutionarily constrained regions or to DNA sequences in chromatin with known activating epigenetic marks. FIREWACh builds upon the fact that almost all CRMs are found in accessible chromatin. Digestion of nuclei with frequently cutting restriction endonucleases releases at least a portion of the accessible chromatin, thereby enriching substantially for regulatory sequences (Figure). Fragments of the released DNA are then cloned into a lentiviral vector upstream of a green fluorescent protein (GFP) reporter gene, and a library of viruses carrying different DNA segments is transduced into an appropriate cell line6. In SIF-seq, DNA fragments are cloned upstream of a yellow fluorescent protein (YFP) gene and inserted into a vector for homologous recombination into the Hprt locus of mouse embryonic stem (ES) cells7. In both methods, cells carrying an active enhancer upstream of the fluorescent reporter gene are isolated by fluorescence activated cell sorting (FACS). The positive cells from each technique contain a single integrant carrying a candidate enhancer. Candidate enhancers can be located by sequencing the integrated DNA from the pool of positive cells and mapping the reads to the genome or target locus. The SIF-seq approach was effective not only in ES cells, but it also was used to discover enhancers active in cardiomyocytes or neural progenitor cells after in vitro differentiation of the ES cells. Figure Two methods for identifying enhancers directly by their activity. FIREWACh starts with DNA fragments cleaved from accessible chromatin, whereas SIF-seq begins with DNA segments from a locus containing a gene of interest. In both methods, the isolated ... Given that the candidate enhancers were discovered by an increased expression of a reporter gene, one expects these new methods to have a very high success rate in identifying active enhancers. This expectation was met by both approaches. Subsequent, independent enhancer assays validated the function of candidate enhancers in 78% of the tested FIREWACh positives and all of the tested SIF-seq positives. This is substantially higher than the results reported when using a MPRA approach10 for enhancer discovery based on histone modifications and motif instances (25% to 41%) or the roughly 50% positive rate of predicted enhancers in moderate throughput assays8. Importantly, several DNA segments associated with epigenetic features indicative of enhancement (such as binding by EP300 or acetylation of histone H3K27) that were inactive in SIF-seq were confirmed to be inactive in an independent assay. While these inactive regions could reflect “opportunistic” binding by TFs and recruitment of chromatin modifiers that does not impact gene regulation, they could also be DNA segments that cooperate with other CRMs in gene regulation but are not independently active. The new methods do have limitations, e.g. they were not designed to be comprehensive. SIF-seq was developed to interrogate in detail regulatory regions around specific loci, using as the input genomic clones in bacterial artificial chromosomes. FIREWACh was not targeted to specific loci, but coverage of all accessible chromatin would require lentiviral libraries larger than is practical. Each method was successful in achieving the goals for which it was designed. In contrast, comprehensive prediction of CRMs still relies on genome-wide maps of epigenetic features associated with regulation, but those candidate CRMs require functional assays. Perhaps future developments will reveal ways to use these activity-based assays in series with the epigenetic maps to accomplish more comprehensive coverage while maintaining high specificity.