Lessons on gene regulation learnt from the Drosophila melanogaster bithorax complex.

Homeotic or Hox genes determine the anterior-posterior body axis in all bilaterians. As expected, Hox genes are highly conserved across bilaterians. Interestingly, however, the peculiar organization of Hox genes in the form of clusters where the order of occurrence of genes in the genome corresponds to the order in which they regulate segmental identity of anterior-posterior body axis is also conserved. The relation between collinearity of arrangement of genes on the chromosomes and spatial function along the body axis has attracted attention to exploring its relevance in the precise regulation of Hox genes. Conservation of genes and their arrangement suggests a linkage between co-regulation and the higher order chromatin organization of the entire complex. To this end, we and others have used Drosophila as the model system to understand the cis-and trans-regulatory components of Hox genes. A number of chromatin-level regulatory elements, like chromatin domain boundaries, and Polycomb Response Elements (PREs) have been discovered in this process. Interestingly, much of what has emerged from the study of homeotic genes, the cis-elements and protein factors, have relevance across the genome in a large number of regulatory events beyond the Hox genes. Here, we review our findings and discuss their genome-wide implications in complex regulatory processes.

Long-Read Genome Sequencing and Assembly of Leptopilina boulardi: A Specialist Drosophila Parasitoid.

Leptopilinaboulardi (Hymenoptera: Figitidae) is a specialist parasitoid of Drosophila The Drosophila-Leptopilina system has emerged as a suitable model for understanding several aspects of host-parasitoid biology. However, a good quality genome of the wasp counterpart was lacking. Here, we report a whole-genome assembly of L. boulardi to bring it in the scope of the applied and fundamental research on Drosophila parasitoids with access to epigenomics and genome editing tools. The 375Mb draft genome has an N50 of 275Kb with 6315 scaffolds >500bp and encompasses >95% complete BUSCOs. Using a combination of ab-initio and RNA-Seq based methods, 25259 protein-coding genes were predicted and 90% (22729) of them could be annotated with at least one function. We demonstrate the quality of the assembled genome by recapitulating the phylogenetic relationship of L. boulardi with other Hymenopterans. The key developmental regulators like Hox genes and sex determination genes are well conserved in L. boulardi, and so is the basic toolkit for epigenetic regulation. The search for epigenetic regulators has also revealed that L. boulardi genome possesses DNMT1 (maintenance DNA methyltransferase), DNMT2 (tRNA methyltransferase) but lacks the de novo DNA methyltransferase (DNMT3). Also, the heterochromatin protein 1 family appears to have expanded as compared to other hymenopterans. The draft genome of L. boulardi (Lb17) will expedite the research on Drosophila parasitoids. This genome resource and early indication of epigenetic aspects in its specialization make it an interesting system to address a variety of questions on host-parasitoid biology.

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.

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.

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.

Chromatin domain boundary element search tool for Drosophila.

Chromatin domain boundary elements prevent inappropriate interaction between distant or closely spaced regulatory elements and restrict enhancers and silencers to correct target promoters. In spite of having such a general role and expected frequent occurrence genome wide, there is no DNA sequence analysis based tool to identify boundary elements. Here, we report chromatin domain Boundary Element Search Tool (cdBEST), to identify boundary elements. cdBEST uses known recognition sequences of boundary interacting proteins and looks for 'motif clusters'. Using cdBEST, we identified boundary sequences across 12 Drosophila species. Of the 4576 boundary sequences identified in Drosophila melanogaster genome, >170 sequences are repetitive in nature and have sequence homology to transposable elements. Analysis of such sequences across 12 Drosophila genomes showed that the occurrence of repetitive sequences in the context of boundaries is a common feature of drosophilids. We use a variety of genome organization criteria and also experimental test on a subset of the cdBEST boundaries in an enhancer-blocking assay and show that 80% of them indeed function as boundaries in vivo. These observations highlight the role of cdBEST in better understanding of chromatin domain boundaries in Drosophila and setting the stage for comparative analysis of boundaries across closely related species.

Transcriptional activation by GAGA factor is through its direct interaction with dmTAF3.

The GAGA factor (GAF), encoded by the Trithorax like gene (Trl) is a multifunctional protein involved in gene activation, Polycomb-dependent repression, chromatin remodeling and is a component of chromatin domain boundaries. Although first isolated as transcriptional activator of the Drosophila homeotic gene Ultrabithorax (Ubx), the molecular basis of this GAF activity is unknown. Here we show that dmTAF3 (also known as BIP2 and dTAF(II)155), a component of TFIID, interacts directly with GAF. We generated mutations in dmTAF3 and show that, in Trl mutant background, they affect transcription of Ubx leading to enhancement of Ubx phenotype. These results reveal that the gene activation pathway involving GAF is through its direct interaction with dmTAF3.

Trl-GAGA directly interacts with lola like and both are part of the repressive complex of Polycomb group of genes.

Epigenetic inheritance to maintain the expression state of the genome is essential during development. In Drosophila, the cis regulatory elements, called the Polycomb Response Elements (PREs) function to mark the epigenetic cellular memory of the corresponding genomic region with the help of PcG and trxG proteins. While the PcG genes code for the repressor proteins, the trxG genes encode activator proteins. The observations that some proteins may function both as PcG and trxG member and that both these group of proteins act upon common cis elements indicate at least a partial functional overlap among these proteins. Trl-GAGA was initially identified as a trxG member but later was shown to be essential for PcG function on several PREs. In order to understand how Trl-GAGA functions in PcG context, we have looked for the interactors of this protein. We identified lola like, aka batman, as a strong interactor of GAGA factor in a yeast two-hybrid screen. lolal also interacts with polyhomeotic and, like Trl, both lolal and ph are needed for iab-7PRE mediated pairing dependent silencing of mini-white transgene. These observations suggest a possible mechanism of how Trl-GAGA plays a role in maintaining the repressed state of target genes involving lolal, which may function as a mediator to recruit PcG complexes.