Chromatin domain boundary element search tool for Drosophila

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Genomic organization of the autonomous regulatory domain of eyeless locus in Drosophila melanogaster

In Drosophila, expression of eyeless (ey) gene is restricted to the developing eyes and central nervous system. However, the flanking genes, myoglianin (myo), and bent (bt) have different temporal and spatial expression patterns as compared to the ey. How distinct regulation of ey is maintained is mostly unknown. Earlier, we have identified a boundary element intervening myo and ey genes (ME boundary) that prevents the crosstalk between the cis-regulatory elements of myo and ey genes. In the present study, we further searched for the cis-elements that define the domain of ey and maintain its expression pattern. We identify another boundary element between ey and bt, the EB boundary. The EB boundary separates the regulatory landscapes of ey and bt genes. The two boundaries, ME and EB, show a long-range interaction as well as interact with the nuclear architecture. This suggests functional autonomy of the ey locus and its insulation from differentially regulated flanking regions. We also identify a new Polycomb Response Element, the ey-PRE, within the ey domain. The expression state of the ey gene, once established during early development is likely to be maintained with the help of ey-PRE. Our study proposes a general regulatory mechanism by which a gene can be maintained in a functionally independent chromatin domain in gene-rich euchromatin.

Insulators in Plants: Progress and Open Questions

The genomes of higher eukaryotes are partitioned into topologically associated domains or TADs, and insulators (also known as boundary elements) are the key elements responsible for their formation and maintenance. Insulators were first identified and extensively studied in Drosophila as well as mammalian genomes, and have also been described in yeast and plants. In addition, many insulator proteins are known in Drosophila, and some have been investigated in mammals. However, much less is known about this important class of non-coding DNA elements in plant genomes. In this review, we take a detailed look at known plant insulators across different species and provide an overview of potential determinants of plant insulator functions, including cis-elements and boundary proteins. We also discuss methods previously used in attempts to identify plant insulators, provide a perspective on their importance for research and biotechnology, and discuss areas of potential future research.

Homeotic Genes: Clustering, Modularity, and Diversity

Hox genes code for transcription factors and are evolutionarily conserved. They regulate a plethora of downstream targets to define the anterior-posterior (AP) body axis of a developing bilaterian embryo. Early work suggested a possible role of clustering and ordering of Hox to regulate their expression in a spatially restricted manner along the AP axis. However, the recent availability of many genome assemblies for different organisms uncovered several examples that defy this constraint. With recent advancements in genomics, the current review discusses the arrangement of Hox in various organisms. Further, we revisit their discovery and regulation in Drosophila melanogaster. We also review their regulation in different arthropods and vertebrates, with a significant focus on Hox expression in the crustacean Parahyale hawaiensis. It is noteworthy that subtle changes in the levels of Hox gene expression can contribute to the development of novel features in an organism. We, therefore, delve into the distinct regulation of these genes during primary axis formation, segment identity, and extra-embryonic roles such as in the formation of hair follicles or misregulation leading to cancer. Toward the end of each section, we emphasize the possibilities of several experiments involving various organisms, owing to the advancements in the field of genomics and CRISPR-based genome engineering. Overall, we present a holistic view of the functioning of Hox in the animal world.

MOCCA: a flexible suite for modelling DNA sequence motif occurrence combinatorics

BACKGROUND

Cis-regulatory elements (CREs) are DNA sequence segments that regulate gene expression. Among CREs are promoters, enhancers, Boundary Elements (BEs) and Polycomb Response Elements (PREs), all of which are enriched in specific sequence motifs that form particular occurrence landscapes. We have recently introduced a hierarchical machine learning approach (SVM-MOCCA) in which Support Vector Machines (SVMs) are applied on the level of individual motif occurrences, modelling local sequence composition, and then combined for the prediction of whole regulatory elements. We used SVM-MOCCA to predict PREs in Drosophila and found that it was superior to other methods. However, we did not publish a polished implementation of SVM-MOCCA, which can be useful for other researchers, and we only tested SVM-MOCCA with IUPAC motifs and PREs.

RESULTS

We here present an expanded suite for modelling CRE sequences in terms of motif occurrence combinatorics-Motif Occurrence Combinatorics Classification Algorithms (MOCCA). MOCCA contains efficient implementations of several modelling methods, including SVM-MOCCA, and a new method, RF-MOCCA, a Random Forest-derivative of SVM-MOCCA. We used SVM-MOCCA and RF-MOCCA to model Drosophila PREs and BEs in cross-validation experiments, making this the first study to model PREs with Random Forests and the first study that applies the hierarchical MOCCA approach to the prediction of BEs. Both models significantly improve generalization to PREs and boundary elements beyond that of previous methods-including 4-spectrum and motif occurrence frequency Support Vector Machines and Random Forests-, with RF-MOCCA yielding the best results.

CONCLUSION

MOCCA is a flexible and powerful suite of tools for the motif-based modelling of CRE sequences in terms of motif composition. MOCCA can be applied to any new CRE modelling problems where motifs have been identified. MOCCA supports IUPAC and Position Weight Matrix (PWM) motifs. For ease of use, MOCCA implements generation of negative training data, and additionally a mode that requires only that the user specifies positives, motifs and a genome. MOCCA is licensed under the MIT license and is available on Github at https://github.com/bjornbredesen/MOCCA .

Interplay of pericentromeric genome organization and chromatin landscape regulates the expression of Drosophila melanogaster heterochromatic genes

Background

Transcription of genes residing within constitutive heterochromatin is paradoxical to the tenets of epigenetic code. The regulatory mechanisms of Drosophila melanogaster heterochromatic gene transcription remain largely unknown. Emerging evidence suggests that genome organization and transcriptional regulation are inter-linked. However, the pericentromeric genome organization is relatively less studied. Therefore, we sought to characterize the pericentromeric genome organization and understand how this organization along with the pericentromeric factors influences heterochromatic gene expression.

Results

Here, we characterized the pericentromeric genome organization in Drosophila melanogaster using 5C sequencing. Heterochromatic topologically associating domains (Het TADs) correlate with distinct epigenomic domains of active and repressed heterochromatic genes at the pericentromeres. These genes are known to depend on the heterochromatic landscape for their expression. However, HP1a or Su(var)3-9 RNAi has minimal effects on heterochromatic gene expression, despite causing significant changes in the global Het TAD organization. Probing further into this observation, we report the role of two other chromatin proteins enriched at the pericentromeres-dMES-4 and dADD1 in regulating the expression of a subset of heterochromatic genes.

Conclusions

Distinct pericentromeric genome organization and chromatin landscapes maintained by the interplay of heterochromatic factors (HP1a, H3K9me3, dMES-4 and dADD1) are sufficient to support heterochromatic gene expression despite the loss of global Het TAD structure. These findings open new avenues for future investigations into the mechanisms of heterochromatic gene expression.

Genomic organization of Polycomb Response Elements and its functional implication in Drosophila and other insects

The epigenetic memory is an essential aspect of multicellular organisms to maintain several cell types and their gene expression pattern. This complex process uses a number of protein factors and specific DNA elements within the developmental cues to achieve this. The protein factors involved in the process are the Polycomb group (PcG) members, and, accordingly, the DNA sequences that interact with these proteins are called Polycomb Response Elements (PREs). Since the PcG proteins are highly conserved among higher eukaryotes, including insects, and function at thousands of sites in the genomes, it is expected that PREs mayalso be present across the genome.However, the studies on PREs in insect species, other thanDrosophila, is currently lacking.We took a bioinformatics approach to develop an inclusive PRE prediction tool, 'PRE Mapper', to address this need. By applying this tool on the Drosophila melanogaster genome, we predicted greater than 20,000 PREs.When comparedwith the available PRE prediction methods, this tool shows far better performance by correctly identifying the in vivo binding sites of PcG proteins, identified by genome-scale ChIP experiments. Further analysis of the predicted PREs shows their cohabitation with chromatin domain boundary elements at several places in the Drosophila genome, possibly defining a composite epigenetic module.We analysed 10 insect genomes in this context and find several conserved features in PREs across the insect species with some variations in their occurrence frequency. These analyses leading to the identification of PREin insect genomes contribute to our understanding of epigenetic mechanisms in these organisms.

O-GlcNAcylation of boundary element associated factor (BEAF 32) in Drosophila melanogaster correlates with active histone marks at the promoters of its target genes

Boundary Element-Associated Factor 32 (BEAF 32) is a sequence specific DNA binding protein involved in functioning of chromatin domain boundaries in Drosophila. Several studies also show it to be involved in transcriptional regulation of a large number of genes, many of which are annotated to have cell cycle, development and differentiation related function. Since post-translational modifications (PTMs) of proteins add to their functional capacity, we investigated the PTMs on BEAF 32. The protein is known to be phosphorylated and O-GlcNAcylated. We mapped O-GlcNAc site at T91 of BEAF 32 and showed that it is linked to the deposition of active histone (H3K4me3) marks at transcription start site (TSS) of associated genes. Its role as a boundary associated factor, however, does not depend on this modification. Our study shows that by virtue of O-GlcNAcylation, BEAF 32 is linked to epigenetic mechanisms that activate a subset of associated genes.

The commonly used eye-specific sev-GAL4 and GMR-GAL4 drivers in Drosophila melanogaster are expressed in tissues other than eyes also

The binary GAL4-UAS system of conditional gene expression is widely used by Drosophila geneticists to target expression of the desired transgene in tissue of interest. In many studies, a preferred target tissue is the Drosophila eye, for which the sev-GAL4 and GMR-GAL4 drivers are most widely used since they are believed to be expressed exclusively in the developing eye cells. However, several reports have noted lethality following expression of certain transgenes under these GAL4 drivers notwithstanding the fact that eye is not essential for survival of the fly. Therefore, to explore the possibility that these drivers may also be active in tissues other than eye, we examined the expression of UAS-GFP reporter driven by the sev-GAL4 or GMR-GAL4 drivers. We found that both these drivers are indeed expressed in additional tissues, including a common set of specific neuronal cells in larval and pupal ventral and cerebral ganglia. Neither sev nor glass gene has so far been reported to be expressed in these neuronal cells. Expression pattern of sev-GAL4 driver parallels that of the endogenous Sevenless protein. In addition to cells in which sev-GAL4 is expressed, the GMR-GAL4 is expressed in several other larval cell types also. Further, two different GMR-GAL4 lines also show some specific differences in their expression domains outside the eye discs. These findings emphasize the need for a careful confirmation of the expression domains of a GAL4 driver being used in a given study, rather than relying only on the empirically claimed expression domains.

Genome-Wide Mapping Targets of the Metazoan Chromatin Remodeling Factor NURF Reveals Nucleosome Remodeling at Enhancers, Core Promoters and Gene Insulators

NURF is a conserved higher eukaryotic ISWI-containing chromatin remodeling complex that catalyzes ATP-dependent nucleosome sliding. By sliding nucleosomes, NURF is able to alter chromatin dynamics to control transcription and genome organization. Previous biochemical and genetic analysis of the specificity-subunit of Drosophila NURF (Nurf301/Enhancer of Bithorax (E(bx)) has defined NURF as a critical regulator of homeotic, heat-shock and steroid-responsive gene transcription. It has been speculated that NURF controls pathway specific transcription by co-operating with sequence-specific transcription factors to remodel chromatin at dedicated enhancers. However, conclusive in vivo demonstration of this is lacking and precise regulatory elements targeted by NURF are poorly defined. To address this, we have generated a comprehensive map of in vivo NURF activity, using MNase-sequencing to determine at base pair resolution NURF target nucleosomes, and ChIP-sequencing to define sites of NURF recruitment. Our data show that, besides anticipated roles at enhancers, NURF interacts physically and functionally with the TRF2/DREF basal transcription factor to organize nucleosomes downstream of active promoters. Moreover, we detect NURF remodeling and recruitment at distal insulator sites, where NURF functionally interacts with and co-localizes with DREF and insulator proteins including CP190 to establish nucleosome-depleted domains. This insulator function of NURF is most apparent at subclasses of insulators that mark the boundaries of chromatin domains, where multiple insulator proteins co-associate. By visualizing the complete repertoire of in vivo NURF chromatin targets, our data provide new insights into how chromatin remodeling can control genome organization and regulatory interactions.

Supervised learning method for predicting chromatin boundary associated insulator elements

In eukaryotic cells, the DNA material is densely packed inside the nucleus in the form of a DNA-protein complex structure called chromatin. Since the actual conformation of the chromatin fiber defines the possible regulatory interactions between genes and their regulatory elements, it is very important to understand the mechanisms governing folding of chromatin. In this paper, we show that supervised methods for predicting chromatin boundary elements are much more effective than the currently popular unsupervised methods. Using boundary locations from published Hi-C experiments and modEncode tracks as features, we can tell the insulator elements from randomly selected background sequences with great accuracy. In addition to accurate predictions of the training boundary elements, our classifiers make new predictions. Many of them correspond to the locations of known insulator elements. The key features used for predicting boundary elements do not depend on the prediction method. Because of its miniscule size, chromatin state cannot be measured directly, we need to rely on indirect measurements, such as ChIP-Seq and fill in the gaps with computational models. Our results show that currently, at least in the model organisms, where we have many measurements including ChIP-Seq and Hi-C, we can make accurate predictions of insulator positions.

Intrinsic Toxicity of Unchecked Heterochromatin Spread Is Suppressed by Redundant Chromatin Boundary Functions in Schizosacchromyces pombe

Effective boundary mechanisms halt the spread of repressive histone methylation. In the fission yeast Schizosacchromyces pombe, two factors/elements required for boundary function have been described, the jmjC protein Epe1 and binding sites for the RNA polymerase III transcription factor TFIIIC. Perplexingly, individual mutation of Epe1 or TFIIIC sites produces only mild boundary defects, and no other boundary factors have been identified. To approach these issues, we developed a synthetic reporter gene tool that uses a tethered Clr4 histone H3K9 methyltransferase and monitors the ability of a DNA element to block heterochromatin spread. The inverted repeat (IR) that flanks the mat2/3 silent mating-type cassette region demonstrates strong boundary activity compared to sequences that flank pericentromeric heterochromatic repeats. Rather than acting in the same inhibitory pathway, Epe1 and TFIIIC sites mediate boundary function of the IR via the two parallel and largely redundant pathways. We also use the system to demonstrate that HP1/Swi6 promotes boundary activity in addition to promoting silencing and acts in the same pathway as Epe1. Inhibition of heterochromatin spread at the endogenous IR element also requires either Epe1 or TFIIIC sites. Strikingly, mutation of both mechanisms results in growth inhibition that is associated with the spread of heterochromatin over many kilobases to the nearest essential gene and the near-complete silencing of several intervening euchromatic genes. The growth defect is reversed by deletion of clr4+, indicating that the redundant boundary mechanisms protect cells from intrinsic toxicity caused by the spread of heterochromatin.

Genome-wide mapping of matrix attachment regions in Drosophila melanogaster

Background

Eukaryotic genome acquires functionality upon proper packaging within the nucleus. This process is facilitated by the structural framework of Nuclear Matrix, a nucleo-proteinaceous meshwork. Matrix Attachment Regions (MARs) in the genome serve as anchoring sites to this framework.

Results

Here we report direct sequencing of the MAR preparation from Drosophila melanogaster embryos and identify >7350 MARs. This amounts to ~2.5% of the fly genome and often coincide with AT rich non-coding regions. We find significant association of MARs with the origins of replication, transcription start sites, paused RNA Polymerase II sites and exons, but not introns, of highly expressed genes. We also identified sequence motifs and repeats that constitute MARs.

Conclusion

Our data reveal the contact points of genome to the nuclear architecture and provide a link between nuclear functions and genomic packaging.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1022) contains supplementary material, which is available to authorized users.

Mining the cis-regulatory elements of Hox clusters

Hox clusters have served as a favorite system to study the role of cis-regulatory elements at multiple layers of gene regulation. Organization and regulation of Hox genes show remarkable conservation and determine the anterior-posterior body axis across the bilaterians. Identification of a variety of regulatory regions within the complex and around it, embedded primarily in the noncoding part of the corresponding genomic region that can spread 100-150 kb, is a challenging problem. Multiple experimental and computational tools need to be employed to investigate functional features of such elements. Here we discuss parallel approaches to mine the most plausible regulatory information from the noncoding sequences of Hox clusters, among diverse species.

Transposable element 'roo' attaches to nuclear matrix of the Drosophila melanogaster

The genome of eukaryotes is organized into structural units of chromatin loops. This higher order organization is supported by a nuclear skeleton called the nuclear matrix. The genomic DNA associated with the nuclear matrix is called the matrix associated region (MAR). Only a few genome-wide screens have been attempted, although many studies have characterized locusspecific MAR DNA sequences. In this study, a MAR DNA library was prepared from the Drosophila melanogaster Meigen (Diptera: Drosophilidae) genome. One of the sequences identified as a MAR was from a long terminal repeat region of 'roo' retrotransposon (roo MAR). Sequence analysis of roo MAR showed its distribution across the D. melanogaster genome. roo MAR also showed high sequence similarity with a previously identified MAR in Drosophila, namely the 'gypsy' retrotransposon. Analysis of the genes flanking roo MAR insertions in the Drosophila genome showed that genes were co-ordinately expressed. The results from the present study in D. melanogaster suggest this sequence plays an important role in genome organization and function. The findings point to an evolutionary role of retrotransposons in shaping the genomic architecture of eukaryotes.

Conserved boundary elements from the Hox complex of mosquito, Anopheles gambiae

The conservation of hox genes as well as their genomic organization across the phyla suggests that this system of anterior–posterior axis formation arose early during evolution and has come under strong selection pressure. Studies in the split Hox cluster of Drosophila have shown that proper expression of hox genes is dependent on chromatin domain boundaries that prevent inappropriate interactions among different types of cis-regulatory elements. To investigate whether boundary function and their role in regulation of hox genes is conserved in insects with intact Hox clusters, we used an algorithm to locate potential boundary elements in the Hox complex of mosquito, Anopheles gambiae. Several potential boundary elements were identified that could be tested for their functional conservation. Comparative analysis revealed that like Drosophila, the bithorax region in A. gambiae contains an extensive array of boundaries and enhancers organized into domains. We analysed a subset of candidate boundary elements and show that they function as enhancer blockers in Drosophila. The functional conservation of boundary elements from mosquito in fly suggests that regulation of hox genes involving chromatin domain boundaries is an evolutionary conserved mechanism and points to an important role of such elements in key developmentally regulated loci.