Functional anatomy of distant-acting mammalian enhancers

Construction of an atlas of transcription factor binding during mouse development identifies popular regulatory regions

Gene regulators physically associate with the genome, in a combinatorial fashion, to drive tissue-specific gene expression. Uncovering the genome-wide activity of all gene regulators across tissues is therefore needed to understand gene regulation during development. Here, we take a first step towards this goal. Using CUT&RUN, we systematically mapped genome-wide binding profiles of key transcription factors and co-factors that mediate ontogenetically relevant signaling pathways in select mouse tissues at two developmental stages. Computation of the datasets unveiled tissue- and time-specific activity for each gene regulator. We identified 'popular' regulatory regions that are bound by a multitude of regulators, which tend to be more evolutionarily conserved. Consistently, they lie near the transcription start site of genes for which dysregulation results in early embryonic lethality. Moreover, the human homologs of these regions are similarly bound by many gene regulators and are highly conserved, indicating a retained relevance for human development. This work constitutes a decisive step towards understanding how the genome is simultaneously read and used by gene regulators in a holistic fashion to drive embryonic development.

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

An association between copy number variation of enhancer involved in craniofacial development and biogeographic ancestry

Human facial morphology is a combination of many complex traits and is determined by a large number of genes and enhancers. Here, we report a Copy Number Variation (CNV) study of enhancer hs1431 in populations of Central European and South Siberian ancestry. Central European samples included 97 Poles, while South Siberian samples included 78 Buryats and 27 Tuvinians. CNVs were detected by real-time PCR, using ViiA™ 7 Real-Time PCR System (Applied Biosystems). We revealed significant differences in CNV of hs1431 enhancer between Polish and Buryat population (p=0.0378), but not between Central European and South Siberian population (p=0.1225). Our results suggest that an increase in copy number variation of hs1431 enhancer is associated with biogeographic ancestry. However, this result needs extending and replicating in larger cohorts. This is the first study revealing the presence of copy number variation of enhancer hs1431 in humans.

Regulation of the Late Onset alzheimer's Disease Associated HLA-DQA1/DRB1 Expression

(Genome-wide Association Studies) GWAS have identified ∼42 late-onset Alzheimer's disease (LOAD)-associated loci, each of which contains multiple single nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) and most of these SNPs are in the non-coding region of human genome. However, how these SNPs regulate risk gene expression remains unknown. In this work, by using a set of novel techniques, we identified 6 functional SNPs (fSNPs) rs9271198, rs9271200, rs9281945, rs9271243, and rs9271247 on the LOAD-associated HLA-DRB1/DQA1 locus and 42 proteins specifically binding to five of these 6 fSNPs. As a proof of evidence, we verified the allele-specific binding of GATA2 and GATA3, ELAVL1 and HNRNPA0, ILF2 and ILF3, NFIB and NFIC, as well as CUX1 to these five fSNPs, respectively. Moreover, we demonstrate that all these nine proteins regulate the expression of both HLA-DQA1 and HLA-DRB1 in human microglial cells. The contribution of HLA class II to the susceptibility of LOAD is discussed.

Post-GWAS functional studies reveal an RA-associated CD40-induced NF-kB signal transduction and transcriptional regulation network targeted by class II HDAC inhibitors

Currently, it remains difficult to identify which single nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWAS) are functional and how various functional SNPs (fSNPs) interact and contribute to disease susceptibility. GWAS have identified a CD40 locus that is associated with rheumatoid arthritis (RA). We previously used two techniques developed in our laboratory, single nucleotide polymorphism-next-generation sequencing (SNP-seq) and flanking restriction enhanced DNA pulldown-mass spectrometry (FREP-MS), to determine that the RA risk gene RBPJ regulates CD40 expression via a fSNP at the RA-associated CD40 locus. In the present work, by applying the same approach, we report the identification of six proteins that regulate RBPJ expression via binding to two fSNPs on the RA-associated RBPJ locus. Using these findings, together with the published data, we constructed an RA-associated signal transduction and transcriptional regulation network (STTRN) that functionally connects multiple RA-associated risk genes via transcriptional regulation networks (TRNs) linked by CD40-induced nuclear factor kappa B (NF-kB) signaling. Remarkably, this STTRN provides insight into the potential mechanism of action for the histone deacetylase inhibitor givinostat, an approved therapy for systemic juvenile idiopathic arthritis. Thus, the generation of disease-associated STTRNs based on post-GWAS functional studies is demonstrated as a novel and effective approach to apply GWAS for mechanistic studies and target identification.

Ultraconserved enhancer function does not require perfect sequence conservation

Ultraconserved enhancer sequences show perfect conservation between human and rodent genomes, suggesting that their functions are highly sensitive to mutation. However, current models of enhancer function do not sufficiently explain this extreme evolutionary constraint. We subjected 23 ultraconserved enhancers to different levels of mutagenesis, collectively introducing 1,547 mutations, and examined their activities in transgenic mouse reporter assays. Overall, we find that the regulatory properties of ultraconserved enhancers are robust to mutation. Upon mutagenesis, nearly all (19/23, 83%) still functioned as enhancers at one developmental stage, as did most of those tested again later in development (5/9, 56%). Replacement of endogenous enhancers with mutated alleles in mice corroborated results of transgenic assays, including the functional resilience of ultraconserved enhancers to mutation. Our findings show that the currently known activities of ultraconserved enhancers do not necessarily require the perfect conservation observed in evolution and suggest that additional regulatory or other functions contribute to their sequence constraint.

Giant lungfish genome elucidates the conquest of land by vertebrates

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

A 6.5-kb intergenic structural variation enhances P450-mediated resistance to pyrethroids in malaria vectors lowering bed net efficacy

Elucidating the complex evolutionary armory that mosquitoes deploy against insecticides is crucial to maintain the effectiveness of insecticide-based interventions. Here, we deciphered the role of a 6.5-kb structural variation (SV) in driving cytochrome P450-mediated pyrethroid resistance in the malaria vector, Anopheles funestus. Whole-genome pooled sequencing detected an intergenic 6.5-kb SV between duplicated CYP6P9a/b P450s in pyrethroid-resistant mosquitoes through a translocation event. Promoter analysis revealed a 17.5-fold higher activity (p < .0001) for the SV- carrying fragment than the SV- free one. Quantitative real-time PCR expression profiling of CYP6P9a/b for each SV genotype supported its role as an enhancer because SV+/SV+ homozygote mosquitoes had a significantly greater expression for both genes than heterozygotes SV+/SV- (1.7- to 2-fold) and homozygotes SV-/SV- (4-to 5-fold). Designing a PCR assay revealed a strong association between this SV and pyrethroid resistance (SV+/SV+ vs. SV-/SV-; odds ratio [OR] = 2,079.4, p < .001). The 6.5-kb SV is present at high frequency in southern Africa (80%-100%) but absent in East/Central/West Africa. Experimental hut trials revealed that homozygote SV mosquitoes had a significantly greater chance to survive exposure to pyrethroid-treated nets (OR 27.7; p < .0001) and to blood feed than susceptible mosquitoes. Furthermore, mosquitoes homozygote-resistant at the three loci (SV+/CYP6P9a_R/CYP6P9b_R) exhibited a higher resistance level, leading to a far superior ability to survive exposure to nets than those homozygotes susceptible at the three loci, revealing a strong additive effect. This study highlights the important role of structural variations in the development of insecticide resistance in malaria vectors and their detrimental impact on the effectiveness of pyrethroid-based nets.

An osmolality/salinity-responsive enhancer 1 (OSRE1) in intron 1 promotes salinity induction of tilapia glutamine synthetase

Euryhaline tilapia (Oreochromis mossambicus) are fish that tolerate a wide salinity range from fresh water to > 3× seawater. Even though the physiological effector mechanisms of osmoregulation that maintain plasma homeostasis in fresh water and seawater fish are well known, the corresponding molecular mechanisms that control switching between hyper- (fresh water) and hypo-osmoregulation (seawater) remain mostly elusive. In this study we show that hyperosmotic induction of glutamine synthetase represents a prominent part of this switch. Proteomics analysis of the O. mossambicus OmB cell line revealed that glutamine synthetase is transcriptionally regulated by hyperosmolality. Therefore, the 5' regulatory sequence of O. mossambicus glutamine synthetase was investigated. Using an enhancer trapping assay, we discovered a novel osmosensitive mechanism by which intron 1 positively mediates glutamine synthetase transcription. Intron 1 includes a single, functional copy of an osmoresponsive element, osmolality/salinity-responsive enhancer 1 (OSRE1). Unlike for conventional enhancers, the hyperosmotic induction of glutamine synthetase by intron 1 is position dependent. But irrespective of intron 1 position, OSRE1 deletion from intron 1 abolishes hyperosmotic enhancer activity. These findings indicate that proper intron 1 positioning and the presence of an OSRE1 in intron 1 are required for precise enhancement of hyperosmotic glutamine synthetase expression.

A sequential methodology for the rapid identification and characterization of breast cancer-associated functional SNPs

GWAS cannot identify functional SNPs (fSNP) from disease-associated SNPs in linkage disequilibrium (LD). Here, we report developing three sequential methodologies including Reel-seq (Regulatory element-sequencing) to identify fSNPs in a high-throughput fashion, SDCP-MS (SNP-specific DNA competition pulldown-mass spectrometry) to identify fSNP-bound proteins and AIDP-Wb (allele-imbalanced DNA pulldown-Western blot) to detect allele-specific protein:fSNP binding. We first apply Reel-seq to screen a library containing 4316 breast cancer-associated SNPs and identify 521 candidate fSNPs. As proof of principle, we verify candidate fSNPs on three well-characterized loci: FGFR2, MAP3K1 and BABAM1. Next, using SDCP-MS and AIDP-Wb, we rapidly identify multiple regulatory factors that specifically bind in an allele-imbalanced manner to the fSNPs on the FGFR2 locus. We finally demonstrate that the factors identified by SDCP-MS can regulate risk gene expression. These data suggest that the sequential application of Reel-seq, SDCP-MS, and AIDP-Wb can greatly help to translate large sets of GWAS data into biologically relevant information.

Interrogating the Grainyhead-like 2 (Grhl2) genomic locus identifies an enhancer element that regulates palatogenesis in mouse

The highly-conserved Grainyhead-like (Grhl) transcription factors are critical regulators of embryogenesis that regulate cellular survival, proliferation, migration and epithelial integrity, especially during the formation of the craniofacial skeleton. Family member Grhl2 is expressed throughout epithelial tissues during development, and loss of Grhl2 function leads to significant defects in neurulation, abdominal wall closure, formation of the face and fusion of the maxilla/palate. Whereas numerous downstream target genes of Grhl2 have been identified, very little is known about how this crucial developmental transcription factor itself is regulated. Here, using in silico and in utero expression analyses and functional deletion in mice, we have identified a novel 2.4 ​kb enhancer element (mm1286) that drives reporter gene expression in a pattern that strongly recapitulates endogenous Grhl2 in the craniofacial primordia, modulates Grhl2 expression in these tissues, and augments Grhl2-mediated closure of the secondary palate. Deletion of this genomic element, in the context of inactivation of one allele of Grhl2 (through generation of double heterozygous Grhl2^+/-;mm1286^+/- mice), results in a significant predisposition to palatal clefting at birth. Moreover, we found that a highly conserved 325 bp region of mm1286 is both necessary and sufficient for mediating the craniofacial-specific enhancer activity of this region, and that an extremely well-conserved 12-bp sequence within this element (CTGTCAAACAGGT) substantially determines full enhancer function. Together, these data provide valuable new insights into the upstream genomic regulatory landscape responsible for transcriptional control of Grhl2 during palatal closure.

100 Years of evolving gene-disease complexities and scientific debutants

It's been over 100 years since the word `gene' is around and progressively evolving in several scientific directions. Time-to-time technological advancements have heavily revolutionized the field of genomics, especially when it's about, e.g. triple code development, gene number proposition, genetic mapping, data banks, gene-disease maps, catalogs of human genes and genetic disorders, CRISPR/Cas9, big data and next generation sequencing, etc. In this manuscript, we present the progress of genomics from pea plant genetics to the human genome project and highlight the molecular, technical and computational developments. Studying genome and epigenome led to the fundamentals of development and progression of human diseases, which includes chromosomal, monogenic, multifactorial and mitochondrial diseases. World Health Organization has classified, standardized and maintained all human diseases, when many academic and commercial online systems are sharing information about genes and linking to associated diseases. To efficiently fathom the wealth of this biological data, there is a crucial need to generate appropriate gene annotation repositories and resources. Our focus has been how many gene-disease databases are available worldwide and which sources are authentic, timely updated and recommended for research and clinical purposes. In this manuscript, we have discussed and compared 43 such databases and bioinformatics applications, which enable users to connect, explore and, if possible, download gene-disease data.

TCF/LEF dependent and independent transcriptional regulation of Wnt/β-catenin target genes

During canonical Wnt signalling, the activity of nuclear β-catenin is largely mediated by the TCF/LEF family of transcription factors. To challenge this view, we used the CRISPR/Cas9 genome editing approach to generate HEK 293T cell clones lacking all four TCF/LEF genes. By performing unbiased whole transcriptome sequencing analysis, we found that a subset of β-catenin transcriptional targets did not require TCF/LEF factors for their regulation. Consistent with this finding, we observed in a genome-wide analysis that β-catenin occupied specific genomic regions in the absence of TCF/LEF Finally, we revealed the existence of a transcriptional activity of β-catenin that specifically appears when TCF/LEF factors are absent, and refer to this as β-catenin-GHOST response. Collectively, this study uncovers a previously neglected modus operandi of β-catenin that bypasses the TCF/LEF transcription factors.

Spatial Principles of Chromatin Architecture Associated With Organ-Specific Gene Regulation

Packaging of the genome in the nucleus is a non-random process that is thought to directly contribute to cell type-specific transcriptomes, although this hypothesis remains untested. Epigenome architecture, as assayed by chromatin conformation capture techniques, such as Hi-C, has recently been described in the mammalian cardiac myocyte and found to be remodeled in the setting of heart failure. In the present study, we sought to determine whether the structural features of the epigenome are conserved between different cell types by investigating Hi-C and RNA-seq data from heart and liver. Investigation of genes with enriched expression in heart or liver revealed nuanced interaction paradigms between organs: first, the log₂ ratios of heart:liver (or liver:heart) intrachromosomal interactions are higher in organ-specific gene sets (p = 0.009), suggesting that organ-specific genes have specialized chromatin structural features. Despite similar number of total interactions between cell types, intrachromosomal interaction profiles in heart but not liver demonstrate that genes forming promoter-to-transcription-end-site loops in the cardiac nucleus tend to be involved in cardiac-related pathways. The same analysis revealed an analogous organ-specific interaction profile for liver-specific loop genes. Investigation of A/B compartmentalization (marker of chromatin accessibility) revealed that in the heart, 66.7% of cardiac-specific genes are in compartment A, while 66.1% of liver-specific genes are found in compartment B, suggesting that there exists a cardiac chromatin topology that allows for expression of cardiac genes. Analyses of interchromosomal interactions revealed a relationship between interchromosomal interaction count and organ-specific gene localization (p = 2.2 × 10^-16) and that, for both organs, regions of active or inactive chromatin tend to segregate in 3D space (i.e., active with active, inactive with inactive). 3D models of topologically associating domains (TADs) suggest that TADs tend to interact with regions of similar compartmentalization across chromosomes, revealing trans structural interactions contributing to genomic compartmentalization at distinct structural scales. These models reveal discordant nuclear compaction strategies, with heart packaging compartment A genes preferentially toward the center of the nucleus and liver exhibiting preferential arrangement toward the periphery. Taken together, our data suggest that intra- and interchromosomal chromatin architecture plays a role in orchestrating tissue-specific gene expression.

Chromosomal Conformations and Epigenomic Regulation in Schizophrenia

Chromosomal conformations, including promoter-enhancer loops, provide a critical regulatory layer for the transcriptional machinery. Therefore, schizophrenia, a common psychiatric disorder associated with broad changes in neuronal gene expression in prefrontal cortex and other brain regions implicated in psychosis, could be associated with alterations in higher-order chromatin. Here, we review early studies on spatial genome organization in the schizophrenia postmortem brain and discuss how integrative approaches using cell culture and animal model systems could gain deeper insight into the potential roles of higher-order chromatin for the neurobiology of and novel treatment avenues for common psychiatric disease.

Human evolution: the non-coding revolution

What made us human? Gene expression changes clearly played a significant part in human evolution, but pinpointing the causal regulatory mutations is hard. Comparative genomics enabled the identification of human accelerated regions (HARs) and other human-specific genome sequences. The major challenge in the past decade has been to link diverged sequences to uniquely human biology. This review discusses approaches to this problem, progress made at the molecular level, and prospects for moving towards genetic causes for uniquely human biology.

Quantifying deleterious effects of regulatory variants

The majority of genome-wide association study (GWAS) risk variants reside in non-coding DNA sequences. Understanding how these sequence modifications lead to transcriptional alterations and cell-to-cell variability can help unraveling genotype-phenotype relationships. Here, we describe a computational method, dubbed CAPE, which calculates the likelihood of a genetic variant deactivating enhancers by disrupting the binding of transcription factors (TFs) in a given cellular context. CAPE learns sequence signatures associated with putative enhancers originating from large-scale sequencing experiments (such as ChIP-seq or DNase-seq) and models the change in enhancer signature upon a single nucleotide substitution. CAPE accurately identifies causative cis-regulatory variation including expression quantitative trait loci (eQTLs) and DNase I sensitivity quantitative trait loci (dsQTLs) in a tissue-specific manner with precision superior to several currently available methods. The presented method can be trained on any tissue-specific dataset of enhancers and known functional variants and applied to prioritize disease-associated variants in the corresponding tissue.

Spatial genome organization and cognition

Nonrandom chromosomal conformations, including promoter-enhancer loopings that bypass kilobases or megabases of linear genome, provide a crucial layer of transcriptional regulation and move vast amounts of non-coding sequence into the physical proximity of genes that are important for neurodevelopment, cognition and behaviour. Activity-regulated changes in the neuronal '3D genome' could govern transcriptional mechanisms associated with learning and plasticity, and loop-bound intergenic and intronic non-coding sequences have been implicated in psychiatric and adult-onset neurodegenerative disease. Recent studies have begun to clarify the roles of spatial genome organization in normal and abnormal cognition.

Quantitative perturbation-based analysis of gene expression predicts enhancer activity in early Drosophila embryo

Enhancers constitute one of the major components of regulatory machinery of metazoans. Although several genome-wide studies have focused on finding and locating enhancers in the genomes, the fundamental principles governing their internal architecture and cis-regulatory grammar remain elusive. Here, we describe an extensive, quantitative perturbation analysis targeting the dorsal-ventral patterning gene regulatory network (GRN) controlled by Drosophila NF-κB homolog Dorsal. To understand transcription factor interactions on enhancers, we employed an ensemble of mathematical models, testing effects of cooperativity, repression, and factor potency. Models trained on the dataset correctly predict activity of evolutionarily divergent regulatory regions, providing insights into spatial relationships between repressor and activator binding sites. Importantly, the collective predictions of sets of models were effective at novel enhancer identification and characterization. Our study demonstrates how experimental dataset and modeling can be effectively combined to provide quantitative insights into cis-regulatory information on a genome-wide scale.

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

DNA contains regions known as genes, which may be “transcribed” to produce the RNA molecules that act as templates for building proteins and regulate cell activity. Proteins called transcription factors can bind to specific sequences of DNA to influence whether nearby genes are transcribed. For example, so-called enhancer regions of DNA contain several binding sites for transcription factors, and this binding activates gene transcription. Little is known about how the transcription factor binding sites are organized in enhancer regions, which makes it difficult to use DNA sequence information alone to predict the regulation of genes.

A transcription factor called Dorsal controls the activity of a network of genes that plays a crucial role in the development of fruit fly embryos. Dorsal binds to the enhancer region of a gene called rhomboid, which has been well studied and is known to be a fairly typical example of an enhancer region.

To understand the regulatory information encoded in the DNA sequences of enhancers, Sayal, Dresch et al. have now used a technique called perturbation analysis to investigate the interactions that are likely to occur between Dorsal and other transcription factors as they bind to the rhomboid enhancer. This technique involves systematically mutating the enhancer to remove different combinations of transcription factor binding sites and quantitatively investigating the effect this has on gene activity. A large set of mathematical models were then trained using this data and shown to correctly predict the activity of a range of other gene regulatory regions. The collective predictions of the models identified new enhancer regions and revealed details about how different types of transcription factor binding sites are arranged within enhancers.

As we enter an era where the DNA sequences of entire human populations are increasingly accessible, we would like to know the functional significance of changes in gene regulatory regions. Sayal, Dresch et al. show that the regulatory properties of specific control proteins are accessible by employing quantitative experiments and mathematical models. Similar studies will be required to learn how mutations found across the genome may alter gene expression, leading to better diagnosis and treatment of disease.

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

A Mammalian enhancer trap resource for discovering and manipulating neuronal cell types

There is a continuing need for driver strains to enable cell-type-specific manipulation in the nervous system. Each cell type expresses a unique set of genes, and recapitulating expression of marker genes by BAC transgenesis or knock-in has generated useful transgenic mouse lines. However, since genes are often expressed in many cell types, many of these lines have relatively broad expression patterns. We report an alternative transgenic approach capturing distal enhancers for more focused expression. We identified an enhancer trap probe often producing restricted reporter expression and developed efficient enhancer trap screening with the PiggyBac transposon. We established more than 200 lines and found many lines that label small subsets of neurons in brain substructures, including known and novel cell types. Images and other information about each line are available online (enhancertrap.bio.brandeis.edu).

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

Scientists can track and even alter the activity of different kinds of neurons, as well as the connections between neurons, by manipulating their genes. However, most genes are active in many different kinds of cells in many different places in the brain, making it difficult to track or target only a particular neuron or brain area.

Enhancers are sections of DNA that can regulate the activity of nearby genes so that they are only active in very specific cell types, and an “enhancer trap” is a genetic approach that essentially hijacks enhancers to express artificial genes in those same cell types. The technique relies on inserting a genetic marker, which can be easily tracked, into random locations in the genome. If this marker then interacts with an enhancer, it is activated and the effect of the enhancer on gene expression can be assessed.

This method has been used in fruit flies and fish to identify enhancers that specifically restrict gene expression to a small subset of cells. Now, Shima et al. show that enhancer traps can be used successfully in mammals too. The experiments produced over 200 different strains of mice, many with the fluorescent marker only in specific brain areas or in specific kinds of brain cells. Some of the types of brain cells uncovered by these experiments are new, and the labelling of specific brain cells and brain areas in different strains makes these mice a useful resource for future work. Furthermore, it will be relatively straightforward to produce many more strains of these mice, because it would simply involve crossbreeding mice. It is likely that some of these to-be-discovered strains will be useful tools for research as well.

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

Deciphering the importance of the palindromic architecture of the immunoglobulin heavy-chain 3' regulatory region

The IgH 3' regulatory region (3'RR) controls class switch recombination (CSR) and somatic hypermutation (SHM) in B cells. The mouse 3'RR contains four enhancer elements with hs1,2 flanked by inverted repeated sequences and the centre of a 25-kb palindrome bounded by two hs3 enhancer inverted copies (hs3a and hs3b). hs4 lies downstream of the palindrome. In mammals, evolution maintained this unique palindromic arrangement, suggesting that it is functionally significant. Here we report that deconstructing the palindromic IgH 3'RR strongly affects its function even when enhancers are preserved. CSR and IgH transcription appear to be poorly dependent on the 3'RR architecture and it is more or less preserved, provided 3'RR enhancers are present. By contrast, a ‘palindromic effect' significantly lowers V_H germline transcription, AID recruitment and SHM. In conclusion, this work indicates that the IgH 3'RR does not simply pile up enhancer units but also optimally exposes them into a functional architecture of crucial importance.

The IgH 3' regulatory region contains an evolutionarily conserved palindromic sequence flanking important enhancer elements. Here the authors show that the palindrome is required for generating antibody diversity.

Genes expressed in mouse cortical progenitors are enriched in Pax, Lhx, and Sox transcription factor putative binding sites

Considerable progress has been made in the understanding of molecular and cellular mechanisms controlling the development of the mammalian cortex. The proliferative and neurogenic properties of cortical progenitors located in the ventricular germinal zone start being understood. Little is known however on the cis-regulatory control that finely tunes gene expression in these progenitors. Here, we undertook an in silico-based approach to address this question, followed by some functional validation. Using the Eurexpress database, we established a list of 30 genes specifically expressed in the cortical germinal zone, we selected mouse/human conserved non-coding elements (CNEs) around these genes and we performed motif-enrichment search in these CNEs. We found an over-representation of motifs corresponding to binding sites for Pax, Sox, and Lhx transcription factors, often found as pairs and located within 100bp windows. A small subset of CNEs (n=7) was tested for enhancer activity, by ex-vivo and in utero electroporation assays. Two showed strong enhancer activity in the germinal zone progenitors. Mutagenesis experiments on a selected CNE showed the functional importance of the Pax, Sox, and Lhx TFBS for conferring enhancer activity to the CNE. Overall, from a cis-regulatory viewpoint, our data suggest an input from Pax, Sox and Lhx transcription factors to orchestrate corticogenesis. These results are discussed with regards to the known functional roles of Pax6, Sox2 and Lhx2 in cortical development.

Human Enhancers Are Fragile and Prone to Deactivating Mutations

To explore the underlying mechanisms whereby noncoding variants affect transcriptional regulation, we identified nucleotides capable of disrupting binding of transcription factors and deactivating enhancers if mutated (dubbed candidate killer mutations or KMs) in HepG2 enhancers. On average, approximately 11% of enhancer positions are prone to KMs. A comparable number of enhancer positions are capable of creating de novo binding sites via a single-nucleotide mutation (dubbed candidate restoration mutations or RSs). Both KM and RS positions are evolutionarily conserved and tend to form clusters within an enhancer. We observed that KMs have the most deleterious effect on enhancer activity. In contrast, RSs have a smaller effect in increasing enhancer activity. Additionally, the KMs are strongly associated with liver-related Genome Wide Association Study traits compared with other HepG2 enhancer regions. By applying our framework to lymphoblastoid cell lines, we found that KMs underlie differential binding of transcription factors and differential local chromatin accessibility. The gene expression quantitative trait loci associated with the tissue-specific genes are strongly enriched in KM positions. In summary, we conclude that the KMs have the greatest impact on the level of gene expression and are likely to be the causal variants of tissue-specific gene expression and disease predisposition.

Conserved higher-order chromatin regulates NMDA receptor gene expression and cognition

Three-dimensional chromosomal conformations regulate transcription by moving enhancers and regulatory elements into spatial proximity with target genes. Here we describe activity-regulated long-range loopings bypassing up to 0.5 Mb of linear genome to modulate NMDA glutamate receptor GRIN2B expression in human and mouse prefrontal cortex. Distal intronic and 3' intergenic loop formations competed with repressor elements to access promoter-proximal sequences, and facilitated expression via a "cargo" of AP-1 and NRF-1 transcription factors and TALE-based transcriptional activators. Neuronal deletion or overexpression of Kmt2a/Mll1 H3K4- and Kmt1e/Setdb1 H3K9-methyltransferase was associated with higher-order chromatin changes at distal regulatory Grin2b sequences and impairments in working memory. Genetic polymorphisms and isogenic deletions of loop-bound sequences conferred liability for cognitive performance and decreased GRIN2B expression. Dynamic regulation of chromosomal conformations emerges as a novel layer for transcriptional mechanisms impacting neuronal signaling and cognition.

In pursuit of design principles of regulatory sequences

Instructions for when, where and to what level each gene should be expressed are encoded within regulatory sequences. The importance of motifs recognized by DNA-binding regulators has long been known, but their extensive characterization afforded by recent technologies only partly accounts for how regulatory instructions are encoded in the genome. Here, we review recent advances in our understanding of regulatory sequences that influence transcription and go beyond the description of motifs. We discuss how understanding different aspects of the sequence-encoded regulation can help to unravel the genotype-phenotype relationship, which would lead to a more accurate and mechanistic interpretation of personal genome sequences.

Evolution of transcription factor binding in metazoans - mechanisms and functional implications

Differences in transcription factor binding can contribute to organismal evolution by altering downstream gene expression programmes. Genome-wide studies in Drosophila melanogaster and mammals have revealed common quantitative and combinatorial properties of in vivo DNA binding, as well as marked differences in the rate and mechanisms of evolution of transcription factor binding in metazoans. Here, we review the recently discovered rapid 're-wiring' of in vivo transcription factor binding between related metazoan species and summarize general principles underlying the observed patterns of evolution. We then consider what might explain the differences in genome evolution between metazoan phyla and outline the conceptual and technological challenges facing this research field.

Rapid and pervasive changes in genome-wide enhancer usage during mammalian development

Enhancers are distal regulatory elements that can activate tissue-specific gene expression and are abundant throughout mammalian genomes. Although substantial progress has been made toward genome-wide annotation of mammalian enhancers, their temporal activity patterns and global contributions in the context of developmental in vivo processes remain poorly explored. Here we used epigenomic profiling for H3K27ac, a mark of active enhancers, coupled to transgenic mouse assays to examine the genome-wide utilization of enhancers in three different mouse tissues across seven developmental stages. The majority of the ∼90,000 enhancers identified exhibited tightly temporally restricted predicted activity windows and were associated with stage-specific biological functions and regulatory pathways in individual tissues. Comparative genomic analysis revealed that evolutionary conservation of enhancers decreases following midgestation across all tissues examined. The dynamic enhancer activities uncovered in this study illuminate rapid and pervasive temporal in vivo changes in enhancer usage that underlie processes central to development and disease.