Genome-scale functional characterization of Drosophila developmental enhancers in vivo

Ever-Changing Landscapes: Transcriptional Enhancers in Development and Evolution

A class of cis-regulatory elements, called enhancers, play a central role in orchestrating spatiotemporally precise gene-expression programs during development. Consequently, divergence in enhancer sequence and activity is thought to be an important mediator of inter- and intra-species phenotypic variation. Here, we give an overview of emerging principles of enhancer function, current models of enhancer architecture, genomic substrates from which enhancers emerge during evolution, and the influence of three-dimensional genome organization on long-range gene regulation. We discuss intricate relationships between distinct elements within complex regulatory landscapes and consider their potential impact on specificity and robustness of transcriptional regulation.

Regulatory Enhancer-Core-Promoter Communication via Transcription Factors and Cofactors

Gene expression is regulated by genomic enhancers that recruit transcription factors and cofactors to activate transcription from target core promoters. Over the past years, thousands of enhancers and core promoters in animal genomes have been annotated, and we have learned much about the domain structure in which regulatory genomes are organized in animals. Enhancer-core-promoter targeting occurs at several levels, including regulatory domains, DNA accessibility, and sequence-encoded core-promoter specificities that are likely mediated by different regulatory proteins. We review here current knowledge about enhancer-core-promoter targeting, regulatory communication between enhancers and core promoters, and the protein factors involved. We conclude with an outlook on open questions that we find particularly interesting and that will likely lead to additional insights in the upcoming years.

Plant Enhancers: A Call for Discovery

Higher eukaryotes typically contain many different cell types, displaying different cellular functions that are influenced by biotic and abiotic cues. The different functions are characterized by specific gene expression patterns mediated by regulatory sequences such as transcriptional enhancers. Recent genome-wide approaches have identified thousands of enhancers in animals, reviving interest in enhancers in gene regulation. Although the regulatory roles of plant enhancers are as crucial as those in animals, genome-wide approaches have only very recently been applied to plants. Here we review characteristics of enhancers at the DNA and chromatin level in plants and other species, their similarities and differences, and techniques widely used for genome-wide discovery of enhancers in animal systems that can be implemented in plants.

Genome-wide identification of Drosophila dorso-ventral enhancers by differential histone acetylation analysis

BACKGROUND

Drosophila dorso-ventral (DV) patterning is one of the best-understood regulatory networks to date, and illustrates the fundamental role of enhancers in controlling patterning, cell fate specification, and morphogenesis during development. Histone acetylation such as H3K27ac is an excellent marker for active enhancers, but it is challenging to obtain precise locations for enhancers as the highest levels of this modification flank the enhancer regions. How to best identify tissue-specific enhancers in a developmental system de novo with a minimal set of data is still unclear.

RESULTS

Using DV patterning as a test system, we develop a simple and effective method to identify tissue-specific enhancers de novo. We sample a broad set of candidate enhancer regions using data on CREB-binding protein co-factor binding or ATAC-seq chromatin accessibility, and then identify those regions with significant differences in histone acetylation between tissues. This method identifies hundreds of novel DV enhancers and outperforms ChIP-seq data of relevant transcription factors when benchmarked with mRNA expression data and transgenic reporter assays. These DV enhancers allow the de novo discovery of the relevant transcription factor motifs involved in DV patterning and contain additional motifs that are evolutionarily conserved and for which the corresponding transcription factors are expressed in a DV-biased fashion. Finally, we identify novel target genes of the regulatory network, implicating morphogenesis genes as early targets of DV patterning.

CONCLUSIONS

Taken together, our approach has expanded our knowledge of the DV patterning network even further and is a general method to identify enhancers in any developmental system, including mammalian development.

Genome-Wide Ultrabithorax Binding Analysis Reveals Highly Targeted Genomic Loci at Developmental Regulators and a Potential Connection to Polycomb-Mediated Regulation

Hox homeodomain transcription factors are key regulators of animal development. They specify the identity of segments along the anterior-posterior body axis in metazoans by controlling the expression of diverse downstream targets, including transcription factors and signaling pathway components. The Drosophila melanogaster Hox factor Ultrabithorax (Ubx) directs the development of thoracic and abdominal segments and appendages, and loss of Ubx function can lead for example to the transformation of third thoracic segment appendages (e.g. halters) into second thoracic segment appendages (e.g. wings), resulting in a characteristic four-wing phenotype. Here we present a Drosophila melanogaster strain with a V5-epitope tagged Ubx allele, which we employed to obtain a high quality genome-wide map of Ubx binding sites using ChIP-seq. We confirm the sensitivity of the V5 ChIP-seq by recovering 7/8 of well-studied Ubx-dependent cis-regulatory regions. Moreover, we show that Ubx binding is predictive of enhancer activity as suggested by comparison with a genome-scale resource of in vivo tested enhancer candidates. We observed densely clustered Ubx binding sites at 12 extended genomic loci that included ANTP-C, BX-C, Polycomb complex genes, and other regulators and the clustered binding sites were frequently active enhancers. Furthermore, Ubx binding was detected at known Polycomb response elements (PREs) and was associated with significant enrichments of Pc and Pho ChIP signals in contrast to binding sites of other developmental TFs. Together, our results show that Ubx targets developmental regulators via strongly clustered binding sites and allow us to hypothesize that regulation by Ubx might involve Polycomb group proteins to maintain specific regulatory states in cooperative or mutually exclusive fashion, an attractive model that combines two groups of proteins with prominent gene regulatory roles during animal development.

Chromatin remodeling effects on enhancer activity

During organism development, a diversity of cell types emerges with disparate, yet stable profiles of gene expression with distinctive cellular functions. In addition to gene promoters, the genome contains enhancer regulatory sequences, which are implicated in cellular specialization by facilitating cell-type and tissue-specific gene expression. Enhancers are DNA binding elements characterized by highly sophisticated and various mechanisms of action allowing for the specific interaction of general and tissue-specific transcription factors (TFs). However, eukaryotic organisms package their genetic material into chromatin, generating a physical barrier for TFs to interact with their cognate sequences. The ability of TFs to bind DNA regulatory elements is also modulated by changes in the chromatin structure, including histone modifications, histone variants, ATP-dependent chromatin remodeling, and the methylation status of DNA. Furthermore, it has recently been revealed that enhancer sequences are also transcribed into a set of enhancer RNAs with regulatory potential. These interdependent processes act in the context of a complex network of chromatin interactions, which together contributes to a renewed vision of how gene activation is coordinated in a cell-type-dependent manner. In this review, we describe the interplay between genetic and epigenetic aspects associated with enhancers and discuss their possible roles on enhancer function.

Automatic Segmentation of Drosophila Neural Compartments Using GAL4 Expression Data Reveals Novel Visual Pathways

Identifying distinct anatomical structures within the brain and developing genetic tools to target them are fundamental steps for understanding brain function. We hypothesize that enhancer expression patterns can be used to automatically identify functional units such as neuropils and fiber tracts. We used two recent, genome-scale Drosophila GAL4 libraries and associated confocal image datasets to segment large brain regions into smaller subvolumes. Our results (available at https://strawlab.org/braincode) support this hypothesis because regions with well-known anatomy, namely the antennal lobes and central complex, were automatically segmented into familiar compartments. The basis for the structural assignment is clustering of voxels based on patterns of enhancer expression. These initial clusters are agglomerated to make hierarchical predictions of structure. We applied the algorithm to central brain regions receiving input from the optic lobes. Based on the automated segmentation and manual validation, we can identify and provide promising driver lines for 11 previously identified and 14 novel types of visual projection neurons and their associated optic glomeruli. The same strategy can be used in other brain regions and likely other species, including vertebrates.

•

Genome-scale enhancer expression patterns can be used to predict brain structure

•

Automated clustering of images finds known structures such as olfactory glomeruli

•

Results identify GAL4 lines with strong expression in the predicted structures

•

We validate novel predictions to reveal previously undescribed optic glomeruli

Regulation of the BMP Signaling-Responsive Transcriptional Network in the Drosophila Embryo

The BMP signaling pathway has a conserved role in dorsal-ventral axis patterning during embryonic development. In Drosophila, graded BMP signaling is transduced by the Mad transcription factor and opposed by the Brinker repressor. In this study, using the Drosophila embryo as a model, we combine RNA-seq with Mad and Brinker ChIP-seq to decipher the BMP-responsive transcriptional network underpinning differentiation of the dorsal ectoderm during dorsal-ventral axis patterning. We identify multiple new BMP target genes, including positive and negative regulators of EGF signaling. Manipulation of EGF signaling levels by loss- and gain-of-function studies reveals that EGF signaling negatively regulates embryonic BMP-responsive transcription. Therefore, the BMP gene network has a self-regulating property in that it establishes a balance between its activity and that of the antagonistic EGF signaling pathway to facilitate correct patterning. In terms of BMP-dependent transcription, we identify key roles for the Zelda and Zerknüllt transcription factors in establishing the resulting expression domain, and find widespread binding of insulator proteins to the Mad and Brinker-bound genomic regions. Analysis of embryos lacking the BEAF-32 insulator protein shows reduced transcription of a peak BMP target gene and a reduction in the number of amnioserosa cells, the fate specified by peak BMP signaling. We incorporate our findings into a model for Mad-dependent activation, and discuss its relevance to BMP signal interpretation in vertebrates.

Embryogenesis involves the patterning of many different cell fates by a limited number of types of signals. One way that these signals promote a particular cell fate is through the induction of a complex, yet highly reproducible, gene expression programme that instructs changes in the cell. For example, there is a conserved role for BMP signals in specifying cell fates during dorsal-ventral axis patterning. Here, we have used genomics approaches to identify the gene expression programme implemented in response to BMP signaling during axis patterning in the Drosophila embryo. Part of the gene network downstream of BMP signaling includes members of the EGF signaling pathway, with our data highlighting reciprocal interactions between these two pathways. We have also determined genome-wide binding of BMP-responsive transcription factors to gain new insights into how the BMP gene network is activated. Our data reveal roles for specific transcription factors and insulator binding proteins, with the latter traditionally associated with the separation of transcriptional domains. Overall, our data will provide a platform for exploiting the tractability of the Drosophila embryo to determine which features of the network are critical drivers of BMP-induced cell fate changes during embryogenesis.

Gene regulation at a distance: From remote enhancers to 3D regulatory ensembles

Large-scale identification of elements associated with gene expression revealed that many of them are located extremely far from gene transcriptional start sites. We review here the growing evidence that show that distal cis-acting elements provide key instructions to genes, as genetic variation affecting them is growingly identified as an importance source of phenotypic diversity and disease. We discuss the different mechanisms that allow these elements to exert their regulatory functions, in a robust and specific manner, despite the large genomic distances separating them from their target genes. We particularly focus on the role of the structural organization of the genome in guiding such regulatory interactions.

Chromatin Heterogeneity and Distribution of Regulatory Elements in the Late-Replicating Intercalary Heterochromatin Domains of Drosophila melanogaster Chromosomes

Late-replicating domains (intercalary heterochromatin) in the Drosophila genome display a number of features suggesting their organization is quite unique. Typically, they are quite large and encompass clusters of functionally unrelated tissue-specific genes. They correspond to the topologically associating domains and conserved microsynteny blocks. Our study aims at exploring further details of molecular organization of intercalary heterochromatin and has uncovered surprising heterogeneity of chromatin composition in these regions. Using the 4HMM model developed in our group earlier, intercalary heterochromatin regions were found to host chromatin fragments with a particular epigenetic profile. Aquamarine chromatin fragments (spanning 0.67% of late-replicating regions) are characterized as a class of sequences that appear heterogeneous in terms of their decompactization. These fragments are enriched with enhancer sequences and binding sites for insulator proteins. They likely mark the chromatin state that is related to the binding of cis-regulatory proteins. Malachite chromatin fragments (11% of late-replicating regions) appear to function as universal transitional regions between two contrasting chromatin states. Namely, they invariably delimit intercalary heterochromatin regions from the adjacent active chromatin of interbands. Malachite fragments also flank aquamarine fragments embedded in the repressed chromatin of late-replicating regions. Significant enrichment of insulator proteins CP190, SU(HW), and MOD2.2 was observed in malachite chromatin. Neither aquamarine nor malachite chromatin types appear to correlate with the positions of highly conserved non-coding elements (HCNE) that are typically replete in intercalary heterochromatin. Malachite chromatin found on the flanks of intercalary heterochromatin regions tends to replicate earlier than the malachite chromatin embedded in intercalary heterochromatin. In other words, there exists a gradient of replication progressing from the flanks of intercalary heterochromatin regions center-wise. The peculiar organization and features of replication in large late-replicating regions are discussed as possible factors shaping the evolutionary stability of intercalary heterochromatin.

Enhancer Control of Transcriptional Bursting

Transcription is episodic, consisting of a series of discontinuous bursts. Using live-imaging methods and quantitative analysis, we examine transcriptional bursting in living Drosophila embryos. Different developmental enhancers positioned downstream of synthetic reporter genes produce transcriptional bursts with similar amplitudes and duration but generate very different bursting frequencies, with strong enhancers producing more bursts than weak enhancers. Insertion of an insulator reduces the number of bursts and the corresponding level of gene expression, suggesting that enhancer regulation of bursting frequency is a key parameter of gene control in development. We also show that linked reporter genes exhibit coordinated bursting profiles when regulated by a shared enhancer, challenging conventional models of enhancer-promoter looping.

High-throughput cis-regulatory element discovery in the vector mosquito Aedes aegypti

Background

Despite substantial progress in mosquito genomic and genetic research, few cis-regulatory elements (CREs), DNA sequences that control gene expression, have been identified in mosquitoes or other non-model insects. Formaldehyde-assisted isolation of regulatory elements paired with DNA sequencing, FAIRE-seq, is emerging as a powerful new high-throughput tool for global CRE discovery. FAIRE results in the preferential recovery of open chromatin DNA fragments that are not bound by nucleosomes, an evolutionarily conserved indicator of regulatory activity, which are then sequenced. Despite the power of the approach, FAIRE-seq has not yet been applied to the study of non-model insects. In this investigation, we utilized FAIRE-seq to profile open chromatin and identify likely regulatory elements throughout the genome of the human disease vector mosquito Aedes aegypti. We then assessed genetic variation in the regulatory elements of dengue virus susceptible (Moyo-S) and refractory (Moyo-R) mosquito strains.

Results

Analysis of sequence data obtained through next generation sequencing of FAIRE DNA isolated from A. aegypti embryos revealed >121,000 FAIRE peaks (FPs), many of which clustered in the 1 kb 5’ upstream flanking regions of genes known to be expressed at this stage. As expected, known transcription factor consensus binding sites were enriched in the FPs, and of these FoxA1, Hunchback, Gfi, Klf4, MYB/ph3 and Sox9 are most predominant. All of the elements tested in vivo were confirmed to drive gene expression in transgenic Drosophila reporter assays. Of the >13,000 single nucleotide polymorphisms (SNPs) recently identified in dengue virus-susceptible and refractory mosquito strains, 3365 were found to map to FPs.

Conclusion

FAIRE-seq analysis of open chromatin in A. aegypti permitted genome-wide discovery of CREs. The results of this investigation indicate that FAIRE-seq is a powerful tool for identification of regulatory DNA in the genomes of non-model organisms, including human disease vector mosquitoes.

Electronic supplementary material

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

An ancient yet flexible cis-regulatory architecture allows localized Hedgehog tuning by patched/Ptch1

The Hedgehog signaling pathway is part of the ancient developmental-evolutionary animal toolkit. Frequently co-opted to pattern new structures, the pathway is conserved among eumetazoans yet flexible and pleiotropic in its effects. The Hedgehog receptor, Patched, is transcriptionally activated by Hedgehog, providing essential negative feedback in all tissues. Our locus-wide dissections of the cis-regulatory landscapes of fly patched and mouse Ptch1 reveal abundant, diverse enhancers with stage- and tissue-specific expression patterns. The seemingly simple, constitutive Hedgehog response of patched/Ptch1 is driven by a complex regulatory architecture, with batteries of context-specific enhancers engaged in promoter-specific interactions to tune signaling individually in each tissue, without disturbing patterning elsewhere. This structure—one of the oldest cis-regulatory features discovered in animal genomes—explains how patched/Ptch1 can drive dramatic adaptations in animal morphology while maintaining its essential core function. It may also suggest a general model for the evolutionary flexibility of conserved regulators and pathways.

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

Ig Superfamily Ligand and Receptor Pairs Expressed in Synaptic Partners in Drosophila

Information processing relies on precise patterns of synapses between neurons. The cellular recognition mechanisms regulating this specificity are poorly understood. In the medulla of the Drosophila visual system, different neurons form synaptic connections in different layers. Here, we sought to identify candidate cell recognition molecules underlying this specificity. Using RNA sequencing (RNA-seq), we show that neurons with different synaptic specificities express unique combinations of mRNAs encoding hundreds of cell surface and secreted proteins. Using RNA-seq and protein tagging, we demonstrate that 21 paralogs of the Dpr family, a subclass of immunoglobulin (Ig)-domain containing proteins, are expressed in unique combinations in homologous neurons with different layer-specific synaptic connections. Dpr interacting proteins (DIPs), comprising nine paralogs of another subclass of Ig-containing proteins, are expressed in a complementary layer-specific fashion in a subset of synaptic partners. We propose that pairs of Dpr/DIP paralogs contribute to layer-specific patterns of synaptic connectivity.

Stability-driven nonnegative matrix factorization to interpret spatial gene expression and build local gene networks

Spatial gene expression patterns enable the detection of local covariability and are extremely useful for identifying local gene interactions during normal development. The abundance of spatial expression data in recent years has led to the modeling and analysis of regulatory networks. The inherent complexity of such data makes it a challenge to extract biological information. We developed staNMF, a method that combines a scalable implementation of nonnegative matrix factorization (NMF) with a new stability-driven model selection criterion. When applied to a set ofDrosophilaearly embryonic spatial gene expression images, one of the largest datasets of its kind, staNMF identified 21 principal patterns (PP). Providing a compact yet biologically interpretable representation ofDrosophilaexpression patterns, PP are comparable to a fate map generated experimentally by laser ablation and show exceptional promise as a data-driven alternative to manual annotations. Our analysis mapped genes to cell-fate programs and assigned putative biological roles to uncharacterized genes. Finally, we used the PP to generate local transcription factor regulatory networks. Spatially local correlation networks were constructed for six PP that span along the embryonic anterior-posterior axis. Using a two-tail 5% cutoff on correlation, we reproduced 10 of the 11 links in the well-studied gap gene network. The performance of PP with theDrosophiladata suggests that staNMF provides informative decompositions and constitutes a useful computational lens through which to extract biological insight from complex and often noisy gene expression data.

Enhancer-promoter interactions are encoded by complex genomic signatures on looping chromatin

Discriminating the gene target of a distal regulatory element from other nearby transcribed genes is a challenging problem with the potential to illuminate the causal underpinnings of complex diseases. We present TargetFinder, a computational method that reconstructs regulatory landscapes from diverse features along the genome. The resulting models accurately predict individual enhancer-promoter interactions across multiple cell lines with a false discovery rate up to 15 times smaller than that obtained using the closest gene. By evaluating the genomic features driving this accuracy, we uncover interactions between structural proteins, transcription factors, epigenetic modifications, and transcription that together distinguish interacting from non-interacting enhancer-promoter pairs. Most of this signature is not proximal to the enhancers and promoters but instead decorates the looping DNA. We conclude that complex but consistent combinations of marks on the one-dimensional genome encode the three-dimensional structure of fine-scale regulatory interactions.

EMERGE: a flexible modelling framework to predict genomic regulatory elements from genomic signatures

Regulatory DNA elements, short genomic segments that regulate gene expression, have been implicated in developmental disorders and human disease. Despite this clinical urgency, only a small fraction of the regulatory DNA repertoire has been confirmed through reporter gene assays. The overall success rate of functional validation of candidate regulatory elements is low. Moreover, the number and diversity of datasets from which putative regulatory elements can be identified is large and rapidly increasing. We generated a flexible and user-friendly tool to integrate the information from different types of genomic datasets, e.g. ATAC-seq, ChIP-seq, conservation, aiming to increase the ease and success rate of functional prediction. To this end, we developed the EMERGE program that merges all datasets that the user considers informative and uses a logistic regression framework, based on validated functional elements, to set optimal weights to these datasets. ROC curve analysis shows that a combination of datasets leads to improved prediction of tissue-specific enhancers in human, mouse and Drosophila genomes. Functional assays based on this prediction can be expected to have substantially higher success rates. The resulting integrated signal for prediction of functional elements can be plotted in a build-in genome browser or exported for further analysis.

Lineage mapping identifies molecular and architectural similarities between the larval and adult Drosophila central nervous system

Neurogenesis in Drosophila occurs in two phases, embryonic and post-embryonic, in which the same set of neuroblasts give rise to the distinct larval and adult nervous systems, respectively. Here, we identified the embryonic neuroblast origin of the adult neuronal lineages in the ventral nervous system via lineage-specific GAL4 lines and molecular markers. Our lineage mapping revealed that neurons born late in the embryonic phase show axonal morphology and transcription factor profiles that are similar to the neurons born post-embryonically from the same neuroblast. Moreover, we identified three thorax-specific neuroblasts not previously characterized and show that HOX genes confine them to the thoracic segments. Two of these, NB2-3 and NB3-4, generate leg motor neurons. The other neuroblast is novel and appears to have arisen recently during insect evolution. Our findings provide a comprehensive view of neurogenesis and show how proliferation of individual neuroblasts is dictated by temporal and spatial cues.

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

Fruit flies undergo a process called metamorphosis in which they change from a maggot or larva into an adult fly. These two life stages look and behave differently and appear to have strikingly different nervous systems. The relationship between the two nervous systems has been most extensively studied in the ventral nerve cord (which is the equivalent to the spinal cord in humans). Although the ventral nerve cords of a larva and an adult fly look quite different, they are generated by the same set of stem cells known as neuroblasts. This is made possible because the neuroblasts proliferate in two separate phases: the first phase occurs in the embryo to generate the neurons of the larval nervous system, and the second phase occurs in the larva to generate neurons for the adult’s nervous system.

Now, Lacin and Truman have paired each of the neurons in the adult fruit fly’s nerve cord with their corresponding neurons in the nerve cords of fruit fly larvae. This involved identifying the original neuroblasts that gave rise to each of the neurons in both larval and adult fruit flies. The results suggest that most neurons that arise from a given neuroblast produce a similar set of molecules and extend similar nerve fibers, even though they work in two different nervous systems. Since neuroblasts in non-metamorphosing insects proliferate continuously, these findings also suggest that, when metamorphosis evolved, a pause was introduced to create the two separate phases of proliferation without a big effect on the types of neurons generated.

Lacin and Truman then went on to discover three neuroblasts that appear to be unique to the middle (or thoracic) segments of a fruit fly. The experiments reveal that the presence of these neuroblasts depended on specific genes that control the development of animal body plans. Two of these neuroblasts generate the so-called motor neurons that control the movement of a fly’s legs. Flies only have legs on their thoracic segments, so this indicates that the development of new neurons is coordinated with the development of the body plan at the stem cell level. The third neuroblast generates neurons that connect with the leg motor neurons, and Lacin and Truman propose that this neuroblast arose from a copy of a neighboring stem cell. The resulting extra neurons may have enabled finer control over the leg movements required for activities such as walking and grooming.

Following on from this work, it is now possible to investigate how molecular events that occur from the embryonic to the adult stages of a fruit fly’s life control the formation and function of its nervous system.

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

Streamlined scanning for enhancer elements in Drosophila melanogaster

Enhancer elements in most eukaryotic organisms are often positioned at a great distance away from the transcription start site of the gene they regulate. Complex three-dimensional chromatin organization and insulators usually guide and limit the range of an enhancer's regulatory activity to a specific genetic locus. Rigorous testing of an entire genomic locus is often required in order to uncover the complete set of cis-regulatory modules (CRMs) regulating a gene, especially those with complex and dynamic expression patterns. Here we report a fast and efficient method for enhancer element identification by scanning large genomic regions using transgenic reporter genes.

A genome-wide resource for the analysis of protein localisation in Drosophila

The Drosophila genome contains >13000 protein-coding genes, the majority of which remain poorly investigated. Important reasons include the lack of antibodies or reporter constructs to visualise these proteins. Here, we present a genome-wide fosmid library of 10000 GFP-tagged clones, comprising tagged genes and most of their regulatory information. For 880 tagged proteins, we created transgenic lines, and for a total of 207 lines, we assessed protein expression and localisation in ovaries, embryos, pupae or adults by stainings and live imaging approaches. Importantly, we visualised many proteins at endogenous expression levels and found a large fraction of them localising to subcellular compartments. By applying genetic complementation tests, we estimate that about two-thirds of the tagged proteins are functional. Moreover, these tagged proteins enable interaction proteomics from developing pupae and adult flies. Taken together, this resource will boost systematic analysis of protein expression and localisation in various cellular and developmental contexts.

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

The fruit fly Drosophila melanogaster is a popular model organism in biological research. Studies using Drosophila have led to important insights into human biology, because related proteins often fulfil similar roles in flies and humans. Thus, studying the role of a protein in Drosophila can teach us about what it might do in a human.

To fulfil their biological roles, proteins often occupy particular locations inside cells, such as the cell’s nucleus or surface membrane. Many proteins are also only found in specific types of cell, such as neurons or muscle cells. A protein’s location thus provides clues about what it does, however cells contain many thousands of proteins and identifying the location of each one is a herculean task.

Sarov et al. took on this challenge and developed a new resource to study the localisation of all Drosophila proteins during this animal’s development. First, genetic engineering was used to tag thousands of Drosophila proteins with a green fluorescent protein, so that they could be tracked under a microscope.

Sarov et al. tagged about 10000 Drosophila proteins in bacteria, and then introduced almost 900 of them into flies to create genetically modified flies. Each fly line contains an extra copy of the tagged gene that codes for one tagged protein. About two-thirds of these tagged proteins appeared to work normally after they were introduced into flies. Sarov et al. then looked at over 200 of these fly lines in more detail and observed that many of the proteins were found in particular cell types and localized to specific parts of the cells. Video imaging of the tagged proteins in living fruit fly embryos and pupae revealed the proteins’ movements, while other techniques showed which proteins bind to the tagged proteins, and may therefore work together in protein complexes.

This resource is openly available to the community, and so researchers can use it to study their favourite protein and gain new insights into how proteins work and are regulated during Drosophila development. Following on from this work, the next challenge will be to create more flies carrying tagged proteins, and to swap the green fluorescent tag with other experimentally useful tags.

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

Genome-wide profiling of nucleosome sensitivity and chromatin accessibility in Drosophila melanogaster

Nucleosomal DNA is thought to be generally inaccessible to DNA-binding factors, such as micrococcal nuclease (MNase). Here, we digest Drosophila chromatin with high and low concentrations of MNase to reveal two distinct nucleosome types: MNase-sensitive and MNase-resistant. MNase-resistant nucleosomes assemble on sequences depleted of A/T and enriched in G/C-containing dinucleotides, whereas MNase-sensitive nucleosomes form on A/T-rich sequences found at transcription start and termination sites, enhancers and DNase I hypersensitive sites. Estimates of nucleosome formation energies indicate that MNase-sensitive nucleosomes tend to be less stable than MNase-resistant ones. Strikingly, a decrease in cell growth temperature of about 10°C makes MNase-sensitive nucleosomes less accessible, suggesting that observed variations in MNase sensitivity are related to either thermal fluctuations of chromatin fibers or the activity of enzymatic machinery. In the vicinity of active genes and DNase I hypersensitive sites nucleosomes are organized into periodic arrays, likely due to 'phasing' off potential barriers formed by DNA-bound factors or by nucleosomes anchored to their positions through external interactions. The latter idea is substantiated by our biophysical model of nucleosome positioning and energetics, which predicts that nucleosomes immediately downstream of transcription start sites are anchored and recapitulates nucleosome phasing at active genes significantly better than sequence-dependent models.

Single-base pair differences in a shared motif determine differential Rhodopsin expression

The final identity and functional properties of a neuron are specified by terminal differentiation genes, which are controlled by specific motifs in compact regulatory regions. To determine how these sequences integrate inputs from transcription factors that specify cell types, we compared the regulatory mechanism of Drosophila Rhodopsin genes that are expressed in subsets of photoreceptors to that of phototransduction genes that are expressed broadly, in all photoreceptors. Both sets of genes share an 11-base pair (bp) activator motif. Broadly expressed genes contain a palindromic version that mediates expression in all photoreceptors. In contrast, each Rhodopsin exhibits characteristic single-bp substitutions that break the symmetry of the palindrome and generate activator or repressor motifs critical for restricting expression to photoreceptor subsets. Sensory neuron subtypes can therefore evolve through single-bp changes in short regulatory motifs, allowing the discrimination of a wide spectrum of stimuli.

High throughput preparation of fly genomic DNA in 96-well format using a paint-shaker

Sample homogenization is an essential step for genomic DNA extraction, with multiple downstream applications in Molecular Biology. Genotyping hundreds or thousands of samples requires an automation of this homogenization step, and high throughput homogenizer equipment currently costs 7000 euros or more. We present an apparatus for homogenization of individual Drosophila adult flies in 96-well micro-titer dishes, which was built from a small portable paint-shaker (F5 portable paint-shaker, Ushake). Single flies are disrupted in each well that contains extraction buffer and a 4-mm metal ball. Our apparatus can hold up to five 96-well micro-titer plates. Construction of the homogenizer apparatus takes about 3–4 days, and all equipment can be obtained from a home improvement store. The total material cost is approximately 700 euros including the paint-shaker. We tested the performance of our apparatus using the ZR-96 Quick-gDNA™ kit (Zymo Research) homogenization buffer and achieved nearly complete tissue homogenization after 15 minutes of shaking. PCR tests did not detect any cross contamination between samples of neighboring wells. We obtained on average 138 ng of genomic DNA per fly, and DNA quality was adequate for standard PCR applications. In principle, our tissue homogenizer can be used for isolation of DNA suitable for library production and high throughput genotyping by Multiplexed Shotgun Genotyping (MSG), as well as RNA isolation from single flies. The sample adapter can also hold and shake other items, such as centrifuge tubes (15–50 mL) or small bottles.

Using FlyBase, a Database of Drosophila Genes and Genomes

For nearly 25 years, FlyBase (flybase.org) has provided a freely available online database of biological information about Drosophila species, focusing on the model organism D. melanogaster. The need for a centralized, integrated view of Drosophila research has never been greater as advances in genomic, proteomic, and high-throughput technologies add to the quantity and diversity of available data and resources.FlyBase has taken several approaches to respond to these changes in the research landscape. Novel report pages have been generated for new reagent types and physical interaction data; Drosophila models of human disease are now represented and showcased in dedicated Human Disease Model Reports; other integrated reports have been established that bring together related genes, datasets, or reagents; Gene Reports have been revised to improve access to new data types and to highlight functional data; links to external sites have been organized and expanded; and new tools have been developed to display and interrogate all these data, including improved batch processing and bulk file availability. In addition, several new community initiatives have served to enhance interactions between researchers and FlyBase, resulting in direct user contributions and improved feedback.This chapter provides an overview of the data content, organization, and available tools within FlyBase, focusing on recent improvements. We hope it serves as a guide for our diverse user base, enabling efficient and effective exploration of the database and thereby accelerating research discoveries.

A Guide to Genome-Wide In Vivo RNAi Applications in Drosophila

RNAi technologies enable the testing of gene function in a cell-type- and stage-specific manner in Drosophila. The development of genome-wide RNAi libraries has allowed expansion of this approach to the genome scale and supports identification of most genes required for a given process in a cell type of choice. However, a large-scale RNAi approach also harbors many potential pitfalls that can complicate interpretation of the results. Here, we summarize published screens and provide a guide on how to optimally plan and perform a large-scale, in vivo RNAi screen. We highlight the importance of assay design and give suggestions on how to optimize the assay conditions by testing positive and negative control genes. These genes are used to estimate false-negative and false-positive rates of the screen data. We discuss the planning and logistics of a large-scale screen in detail and suggest bioinformatics platforms to identify and select gene groups of interest for secondary assays. Finally, we review various options to confirm RNAi knock-down specificity and thus identify high confidence genes for more detailed case-by-case studies in the future.

Shadow Enhancers Are Pervasive Features of Developmental Regulatory Networks

Embryogenesis is remarkably robust to segregating mutations and environmental variation; under a range of conditions, embryos of a given species develop into stereotypically patterned organisms. Such robustness is thought to be conferred, in part, through elements within regulatory networks that perform similar, redundant tasks. Redundant enhancers (or "shadow" enhancers), for example, can confer precision and robustness to gene expression, at least at individual, well-studied loci. However, the extent to which enhancer redundancy exists and can thereby have a major impact on developmental robustness remains unknown. Here, we systematically assessed this, identifying over 1,000 predicted shadow enhancers during Drosophila mesoderm development. The activity of 23 elements, associated with five genes, was examined in transgenic embryos, while natural structural variation among individuals was used to assess their ability to buffer against genetic variation. Our results reveal three clear properties of enhancer redundancy within developmental systems. First, it is much more pervasive than previously anticipated, with 64% of loci examined having shadow enhancers. Their spatial redundancy is often partial in nature, while the non-overlapping function may explain why these enhancers are maintained within a population. Second, over 70% of loci do not follow the simple situation of having only two shadow enhancers-often there are three (rols), four (CadN and ade5), or five (Traf1), at least one of which can be deleted with no obvious phenotypic effects. Third, although shadow enhancers can buffer variation, patterns of segregating variation suggest that they play a more complex role in development than generally considered.

Morphogen rules: design principles of gradient-mediated embryo patterning

The Drosophila blastoderm and the vertebrate neural tube are archetypal examples of morphogen-patterned tissues that create precise spatial patterns of different cell types. In both tissues, pattern formation is dependent on molecular gradients that emanate from opposite poles. Despite distinct evolutionary origins and differences in time scales, cell biology and molecular players, both tissues exhibit striking similarities in the regulatory systems that establish gene expression patterns that foreshadow the arrangement of cell types. First, signaling gradients establish initial conditions that polarize the tissue, but there is no strict correspondence between specific morphogen thresholds and boundary positions. Second, gradients initiate transcriptional networks that integrate broadly distributed activators and localized repressors to generate patterns of gene expression. Third, the correct positioning of boundaries depends on the temporal and spatial dynamics of the transcriptional networks. These similarities reveal design principles that are likely to be broadly applicable to morphogen-patterned tissues.

CBP binding outside of promoters and enhancers in Drosophila melanogaster

BACKGROUND

CREB-binding protein (CBP, also known as nejire) is a transcriptional co-activator that is conserved in metazoans. CBP plays an important role in embryonic development and cell differentiation and mutations in CBP can lead to various diseases in humans. In addition, CBP and the related p300 protein have successfully been used to predict enhancers in both humans and flies when they occur with monomethylation of histone H3 on lysine 4 (H3K4me1).

RESULTS

Here, we compare CBP chromatin immunoprecipitation sequencing data from Drosophila S2 cells with modENCODE data and show that CBP is bound at genomic sites with a wide range of functions. As expected, we find that CBP is bound at active promoters and enhancers. In addition, we find that the strongest CBP sites in the genome are found at Polycomb response elements embedded in histone H3 lysine 27 trimethylated (H3K27me3) chromatin, where they correlate with binding of the Pho repressive complex. Interestingly, we find that CBP also binds to most insulators in the genome. At a subset of these, CBP may regulate insulating activity, measured as the ability to prevent repressive H3K27 methylation from spreading into adjacent chromatin.

CONCLUSIONS

We conclude that CBP could be involved in a much wider range of functions than has previously been appreciated, including Polycomb repression and insulator activity. In addition, we discuss the possibility that a common role for CBP at all functional elements may be to regulate interactions between distant chromosomal regions and speculate that CBP is controlling higher order chromatin organization.

Dissecting neural pathways for forgetting in Drosophila olfactory aversive memory

Recent studies have identified molecular pathways driving forgetting and supported the notion that forgetting is a biologically active process. The circuit mechanisms of forgetting, however, remain largely unknown. Here we report two sets of Drosophila neurons that account for the rapid forgetting of early olfactory aversive memory. We show that inactivating these neurons inhibits memory decay without altering learning, whereas activating them promotes forgetting. These neurons, including a cluster of dopaminergic neurons (PAM-β'1) and a pair of glutamatergic neurons (MBON-γ4>γ1γ2), terminate in distinct subdomains in the mushroom body and represent parallel neural pathways for regulating forgetting. Interestingly, although activity of these neurons is required for memory decay over time, they are not required for acute forgetting during reversal learning. Our results thus not only establish the presence of multiple neural pathways for forgetting in Drosophila but also suggest the existence of diverse circuit mechanisms of forgetting in different contexts.

Gene regulation during development in the light of topologically associating domains

During embryonic development, complex transcriptional programs govern the precision of gene expression. Many key developmental genes are regulated via cis-regulatory elements that are located far away in the linear genome. How sequences located hundreds of kilobases away from a promoter can influence its activity has been the subject of numerous speculations, which all underline the importance of the 3D-organization of the genome. The recent advent of chromosome conformation capture techniques has put into focus the subdivision of the genome into topologically associating domains (TADs). TADs may influence regulatory activities on multiple levels. The relative invariance of TAD limits across cell types suggests that they may form fixed structural domains that could facilitate and/or confine long-range regulatory interactions. However, most recent studies suggest that interactions within TADs are more variable and dynamic than initially described. Hence, different models are emerging regarding how TADs shape the complex 3D conformations, and thereafter influence the networks of cis-interactions that govern gene expression during development. For further resources related to this article, please visit the WIREs website.

Demystifying the secret mission of enhancers: linking distal regulatory elements to target genes

Enhancers are short regulatory sequences bound by sequence-specific transcription factors and play a major role in the spatiotemporal specificity of gene expression patterns in development and disease. While it is now possible to identify enhancer regions genomewide in both cultured cells and primary tissues using epigenomic approaches, it has been more challenging to develop methods to understand the function of individual enhancers because enhancers are located far from the gene(s) that they regulate. However, it is essential to identify target genes of enhancers not only so that we can understand the role of enhancers in disease but also because this information will assist in the development of future therapeutic options. After reviewing models of enhancer function, we discuss recent methods for identifying target genes of enhancers. First, we describe chromatin structure-based approaches for directly mapping interactions between enhancers and promoters. Second, we describe the use of correlation-based approaches to link enhancer state with the activity of nearby promoters and/or gene expression. Third, we describe how to test the function of specific enhancers experimentally by perturbing enhancer–target relationships using high-throughput reporter assays and genome editing. Finally, we conclude by discussing as yet unanswered questions concerning how enhancers function, how target genes can be identified, and how to distinguish direct from indirect changes in gene expression mediated by individual enhancers.

cis-regulatory analysis of the Drosophila pdm locus reveals a diversity of neural enhancers

Background

One of the major challenges in developmental biology is to understand the regulatory events that generate neuronal diversity. During Drosophila embryonic neural lineage development, cellular temporal identity is established in part by a transcription factor (TF) regulatory network that mediates a cascade of cellular identity decisions. Two of the regulators essential to this network are the POU-domain TFs Nubbin and Pdm-2, encoded by adjacent genes collectively known as pdm. The focus of this study is the discovery and characterization of cis-regulatory DNA that governs their expression.

Results

Phylogenetic footprinting analysis of a 125 kb genomic region that spans the pdm locus identified 116 conserved sequence clusters. To determine which of these regions function as cis-regulatory enhancers that regulate the dynamics of pdm gene expression, we tested each for in vivo enhancer activity during embryonic development and postembryonic neurogenesis. Our screen revealed 77 unique enhancers positioned throughout the noncoding region of the pdm locus. Many of these activated neural-specific gene expression during different developmental stages and many drove expression in overlapping patterns. Sequence comparisons of functionally related enhancers that activate overlapping expression patterns revealed that they share conserved elements that can be predictive of enhancer behavior. To facilitate data accessibility, the results of our analysis are catalogued in cisPatterns, an online database of the structure and function of these and other Drosophila enhancers.

Conclusions

These studies reveal a diversity of modular enhancers that most likely regulate pdm gene expression during embryonic and adult development, highlighting a high level of temporal and spatial expression specificity. In addition, we discovered clusters of functionally related enhancers throughout the pdm locus. A subset of these enhancers share conserved elements including sequences that correspond to known TF DNA binding sites. Although comparative analysis of the nubbin and pdm-2 encoding sequences indicate that these two genes most likely arose from a duplication event, we found only partial evidence of sequence duplication between their enhancers, suggesting that after the putative duplication their cis-regulatory DNA diverged at a higher rate than their coding sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1897-2) contains supplementary material, which is available to authorized users.

Connecting Neural Codes with Behavior in the Auditory System of Drosophila

Brains are optimized for processing ethologically relevant sensory signals. However, few studies have characterized the neural coding mechanisms that underlie the transformation from natural sensory information to behavior. Here, we focus on acoustic communication in Drosophila melanogaster and use computational modeling to link natural courtship song, neuronal codes, and female behavioral responses to song. We show that melanogaster females are sensitive to long timescale song structure (on the order of tens of seconds). From intracellular recordings, we generate models that recapitulate neural responses to acoustic stimuli. We link these neural codes with female behavior by generating model neural responses to natural courtship song. Using a simple decoder, we predict female behavioral responses to the same song stimuli with high accuracy. Our modeling approach reveals how long timescale song features are represented by the Drosophila brain and how neural representations can be decoded to generate behavioral selectivity for acoustic communication signals.

TRAP-rc, Translating Ribosome Affinity Purification from Rare Cell Populations of Drosophila Embryos

Measuring levels of mRNAs in the process of translation in individual cells provides information on the proteins involved in cellular functions at a given point in time. The protocol dubbed Translating Ribosome Affinity Purification (TRAP) is able to capture this mRNA translation process in a cell-type-specific manner. Based on the affinity purification of polysomes carrying a tagged ribosomal subunit, TRAP can be applied to translatome analyses in individual cells, making it possible to compare cell types during the course of developmental processes or to track disease development progress and the impact of potential therapies at molecular level. Here we report an optimized version of the TRAP protocol, called TRAP-rc (rare cells), dedicated to identifying engaged-in-translation RNAs from rare cell populations. TRAP-rc was validated using the Gal4/UAS targeting system in a restricted population of muscle cells in Drosophila embryos. This novel protocol allows the recovery of cell-type-specific RNA in sufficient quantities for global gene expression analytics such as microarrays or RNA-seq. The robustness of the protocol and the large collections of Gal4 drivers make TRAP-rc a highly versatile approach with potential applications in cell-specific genome-wide studies.

Zelda overcomes the high intrinsic nucleosome barrier at enhancers during Drosophila zygotic genome activation

The Drosophila genome activator Vielfaltig (Vfl), also known as Zelda (Zld), is thought to prime enhancers for activation by patterning transcription factors (TFs). Such priming is accompanied by increased chromatin accessibility, but the mechanisms by which this occurs are poorly understood. Here, we analyze the effect of Zld on genome-wide nucleosome occupancy and binding of the patterning TF Dorsal (Dl). Our results show that early enhancers are characterized by an intrinsically high nucleosome barrier. Zld tackles this nucleosome barrier through local depletion of nucleosomes with the effect being dependent on the number and position of Zld motifs. Without Zld, Dl binding decreases at enhancers and redistributes to open regions devoid of enhancer activity. We propose that Zld primes enhancers by lowering the high nucleosome barrier just enough to assist TFs in accessing their binding motifs and promoting spatially controlled enhancer activation if the right patterning TFs are present. We envision that genome activators in general will utilize this mechanism to activate the zygotic genome in a robust and precise manner.

Genome-wide errant targeting by Hairy

Metazoan transcriptional repressors regulate chromatin through diverse histone modifications. Contributions of individual factors to the chromatin landscape in development is difficult to establish, as global surveys reflect multiple changes in regulators. Therefore, we studied the conserved Hairy/Enhancer of Split family repressor Hairy, analyzing histone marks and gene expression in Drosophila embryos. This long-range repressor mediates histone acetylation and methylation in large blocks, with highly context-specific effects on target genes. Most strikingly, Hairy exhibits biochemical activity on many loci that are uncoupled to changes in gene expression. Rather than representing inert binding sites, as suggested for many eukaryotic factors, many regions are targeted errantly by Hairy to modify the chromatin landscape. Our findings emphasize that identification of active cis-regulatory elements must extend beyond the survey of prototypical chromatin marks. We speculate that this errant activity may provide a path for creation of new regulatory elements, facilitating the evolution of novel transcriptional circuits.

JAK/STAT and Hox Dynamic Interactions in an Organogenetic Gene Cascade

Organogenesis is controlled by gene networks activated by upstream selector genes. During development the gene network is activated stepwise, with a sequential deployment of successive transcription factors and signalling molecules that modify the interaction of the elements of the network as the organ forms. Very little is known about the steps leading from the early specification of the cells that form the organ primordium to the moment when a robust gene network is in place. Here we study in detail how a Hox protein induces during early embryogenesis a simple organogenetic cascade that matures into a complex gene network through the activation of feedback and feed forward interaction loops. To address how the network organization changes during development and how the target genes integrate the genetic information it provides, we analyze in Drosophila the induction of posterior spiracle organogenesis by the Hox gene Abdominal-B (Abd-B). Initially, Abd-B activates in the spiracle primordium a cascade of transcription factors and signalling molecules including the JAK/STAT signalling pathway. We find that at later stages STAT activity feeds back directly into Abd-B, initiating the transformation of the Hox cascade into a gene-network. Focusing on crumbs, a spiracle downstream target gene of Abd-B, we analyze how a modular cis regulatory element integrates the dynamic network information set by Abd-B and the JAK/STAT signalling pathway during development. We describe how a Hox induced genetic cascade transforms into a robust gene network during organogenesis due to the repeated interaction of Abd-B and one of its targets, the JAK/STAT signalling cascade. Our results show that in this network STAT functions not just as a direct transcription factor, but also acts as a "counter-repressor", uncovering a novel mode for STAT directed transcriptional regulation.

Organogenesis is controlled by gene networks activated by upstream selector genes. To address how the network organization changes during development and how the target genes integrate the genetic information it provides, we analyze in Drosophila the induction of posterior spiracle organogenesis by the Hox gene Abdominal-B (Abd-B). Initially, Abd-B activates in the spiracle primordium a cascade of transcription factors and signalling molecules including the JAK/STAT pathway. We find that at later stages STAT activity feeds back into Abd-B, initiating the transformation of the Hox cascade into a gene-network. Focusing on a spiracle downstream target gene of Abd-B, we analyze how its cis regulatory elements integrate the dynamic network information set by Abd-B and the JAK/STAT signalling pathway during development. Our results also show that the well known transcription factor STAT can control gene expression as a “counter-repressor”, uncovering an alternative novel mode for STAT directed transcriptional regulation.

Genome-Wide Mapping of Collier In Vivo Binding Sites Highlights Its Hierarchical Position in Different Transcription Regulatory Networks

Collier, the single Drosophila COE (Collier/EBF/Olf-1) transcription factor, is required in several developmental processes, including head patterning and specification of muscle and neuron identity during embryogenesis. To identify direct Collier (Col) targets in different cell types, we used ChIP-seq to map Col binding sites throughout the genome, at mid-embryogenesis. In vivo Col binding peaks were associated to 415 potential direct target genes. Gene Ontology analysis revealed a strong enrichment in proteins with DNA binding and/or transcription-regulatory properties. Characterization of a selection of candidates, using transgenic CRM-reporter assays, identified direct Col targets in dorso-lateral somatic muscles and specific neuron types in the central nervous system. These data brought new evidence that Col direct control of the expression of the transcription regulators apterous and eyes-absent (eya) is critical to specifying neuronal identities. They also showed that cross-regulation between col and eya in muscle progenitor cells is required for specification of muscle identity, revealing a new parallel between the myogenic regulatory networks operating in Drosophila and vertebrates. Col regulation of eya, both in specific muscle and neuronal lineages, may illustrate one mechanism behind the evolutionary diversification of Col biological roles.

Getting Down to Specifics: Profiling Gene Expression and Protein-DNA Interactions in a Cell Type-Specific Manner

The majority of multicellular organisms are comprised of an extraordinary range of cell types, with different properties and gene expression profiles. Understanding what makes each cell type unique and how their individual characteristics are attributed are key questions for both developmental and neurobiologists alike. The brain is an excellent example of the cellular diversity expressed in the majority of eukaryotes. The mouse brain comprises of approximately 75 million neurons varying in morphology, electrophysiology, and preferences for synaptic partners. A powerful process in beginning to pick apart the mechanisms that specify individual characteristics of the cell, as well as their fate, is to profile gene expression patterns, chromatin states, and transcriptional networks in a cell type-specific manner, i.e., only profiling the cells of interest in a particular tissue. Depending on the organism, the questions being investigated, and the material available, certain cell type-specific profiling methods are more suitable than others. This chapter reviews the approaches presently available for selecting and isolating specific cell types and evaluates their key features.

Recent advances in functional assays of transcriptional enhancers

In this special edition of Genomics, we present reviews of the current state of the field in identifying and functionally understanding transcriptional enhancers in cells and developing tissues. Typically several enhancers coordinate the expression of an individual target gene, each controlling that gene's expression in specific cell types at specific times. Until recently, identifying each gene's enhancers had been challenging because enhancers do not occupy prescribed locations relative to their target genes. Recently there have been powerful advances in DNA sequencing and other technologies that make it possible to identify the majority of enhancers in virtually any cell type of interest. The reviews in this edition of Genomics highlight some of these new and powerful approaches.

STARR-seq - principles and applications

Differential gene expression is the basis for cell type diversity in multicellular organisms and the driving force of development and differentiation. It is achieved by cell type-specific transcriptional enhancers, which are genomic DNA sequences that activate the transcription of their target genes. Their identification and characterization is fundamental to our understanding of gene regulation. Features that are associated with enhancer activity, such as regulatory factor binding or histone modifications can predict the location of enhancers. Nonetheless, enhancer activity can only be assessed by transcriptional reporter assays. Over the past years massively parallel reporter assays have been developed for large scale testing of enhancers. In this review we focus on the principles and applications of STARR-seq, a functional assay that quantifies enhancer strengths in complex candidate libraries and thus allows activity-based enhancer identification in entire genomes. We explain how STARR-seq works, discuss current uses and give an outlook to future applications.

Using transgenic reporter assays to functionally characterize enhancers in animals

Enhancers or cis-regulatory modules play an instructive role in regulating gene expression during animal development and in response to the environment. Despite their importance, we only have an incomplete map of enhancers in the genome and our understanding of the mechanisms governing their function is still limited. Recent advances in genomics provided powerful tools to generate genome-wide maps of potential enhancers. However, most of these methods are based on indirect measures of enhancer activity and have to be followed by functional testing. Animal transgenesis has been a valuable method to functionally test and characterize enhancers in vivo. In this review I discuss how different transgenic strategies are utilized to characterize enhancers in model organisms focusing on studies in Drosophila and mouse. I will further discuss recent large-scale transgenic efforts to systematically identify and catalog enhancers as well as highlight the challenges and future directions in the field.

Optimized tools for multicolor stochastic labeling reveal diverse stereotyped cell arrangements in the fly visual system

We describe the development and application of methods for high-throughput neuroanatomy in Drosophila using light microscopy. These tools enable efficient multicolor stochastic labeling of neurons at both low and high densities. Expression of multiple membrane-targeted and distinct epitope-tagged proteins is controlled both by a transcriptional driver and by stochastic, recombinase-mediated excision of transcription-terminating cassettes. This MultiColor FlpOut (MCFO) approach can be used to reveal cell shapes and relative cell positions and to track the progeny of precursor cells through development. Using two different recombinases, the number of cells labeled and the number of color combinations observed in those cells can be controlled separately. We demonstrate the utility of MCFO in a detailed study of diversity and variability of Distal medulla (Dm) neurons, multicolumnar local interneurons in the adult visual system. Similar to many brain regions, the medulla has a repetitive columnar structure that supports parallel information processing together with orthogonal layers of cell processes that enable communication between columns. We find that, within a medulla layer, processes of the cells of a given Dm neuron type form distinct patterns that reflect both the morphology of individual cells and the relative positions of their arbors. These stereotyped cell arrangements differ between cell types and can even differ for the processes of the same cell type in different medulla layers. This unexpected diversity of coverage patterns provides multiple independent ways of integrating visual information across the retinotopic columns and implies the existence of multiple developmental mechanisms that generate these distinct patterns.

Identifying transcriptional cis-regulatory modules in animal genomes

Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.