Characterization of the extended Myb-like DNA-binding domain of trithorax group protein Zeste

Comparative interactome analysis of the PRE DNA-binding factors: purification of the Combgap-, Zeste-, Psq-, and Adf1-associated proteins

The Polycomb group (PcG) and Trithorax group (TrxG) proteins are key epigenetic regulators controlling the silenced and active states of genes in multicellular organisms, respectively. In Drosophila, PcG/TrxG proteins are recruited to the chromatin via binding to specific DNA sequences termed polycomb response elements (PREs). While precise mechanisms of the PcG/TrxG protein recruitment remain unknown, the important role is suggested to belong to sequence-specific DNA-binding factors. At the same time, it was demonstrated that the PRE DNA-binding proteins are not exclusively localized to PREs but can bind other DNA regulatory elements, including enhancers, promoters, and boundaries. To gain an insight into the PRE DNA-binding protein regulatory network, here, using ChIP-seq and immuno-affinity purification coupled to the high-throughput mass spectrometry, we searched for differences in abundance of the Combgap, Zeste, Psq, and Adf1 PRE DNA-binding proteins. While there were no conspicuous differences in co-localization of these proteins with other functional transcription factors, we show that Combgap and Zeste are more tightly associated with the Polycomb repressive complex 1 (PRC1), while Psq interacts strongly with the TrxG proteins, including the BAP SWI/SNF complex. The Adf1 interactome contained Mediator subunits as the top interactors. In addition, Combgap efficiently interacted with AGO2, NELF, and TFIID. Combgap, Psq, and Adf1 have architectural proteins in their networks. We further investigated the existence of direct interactions between different PRE DNA-binding proteins and demonstrated that Combgap-Adf1, Psq-Dsp1, and Pho-Spps can interact in the yeast two-hybrid assay. Overall, our data suggest that Combgap, Psq, Zeste, and Adf1 are associated with the protein complexes implicated in different regulatory activities and indicate their potential multifunctional role in the regulation of transcription.

Comprehensive analysis of 1R- and 2R-MYBs reveals novel genic and protein features, complex organisation, selective expansion and insights into evolutionary tendencies

Myeloblastosis (MYB) family, the largest plant transcription factor family, has been subcategorised based on the number and type of repeats in the MYB domain. In spite of several reports, evolution of MYB genes and repeats remains enigmatic. Brassicaceae members are endowed with complex genomes, including dysploidy because of its unique history with multiple rounds of polyploidisation, genomic fractionations and rearrangements. The present study is an attempt to gain insights into the complexities of MYB family diversity, understand impacts of genome evolution on gene families and develop an evolutionary framework to understand the origin of various subcategories of MYB gene family. We identified and analysed 1129 MYBs that included 1R-, 2R-, 3R- and atypical-MYBs across sixteen species representing protists, fungi, animals and plants and exclude MYB identified from Brassicaceae except Arabidopsis thaliana; in addition, a total of 1137 2R-MYB genes from six Brassicaceae species were also analysed. Comparative analysis revealed predominance of 1R-MYBs in protists, fungi, animals and lower plants. Phylogenetic reconstruction and analysis of selection pressure suggested ancestral nature of R1-type repeat containing 1R-MYBs that might have undergone intragenic duplication to form multi-repeat MYBs. Distinct differences in gene structure between 1R-MYB and 2R-MYBs were observed regarding intron number, the ratio of gene length to coding DNA sequence (CDS) length and the length of exons encoding the MYB domain. Conserved as well as novel and lineage-specific intron phases were identified. Analyses of physicochemical properties revealed drastic differences indicating functional diversification in MYBs. Phylogenetic reconstruction of 1R- and 2R-MYB genes revealed a shared structure-function relationship in clades which was supported when transcriptome data was analysed in silico. Comparative genomics to study distribution pattern and mapping of 2R-MYBs revealed congruency and greater degree of synteny and collinearity among closely related species. Micro-synteny analysis of genomic segments revealed high conservation of genes that are immediately flanking the surrounding tandemly organised 2R-MYBs along with instances of local duplication, reorganisations and genome fractionation. In summary, polyploidy, dysploidy, reshuffling and genome fractionation were found to cause loss or gain of 2R-MYB genes. The findings need to be supported with functional validation to understand gene structure-function relationship along the evolutionary lineage and adaptive strategies based on comparative functional genomics in plants.

Structure of Zeste-DNA Complex Reveals a New Modality of DNA Recognition by Homeodomain-Like Proteins

Drosophila Zeste is a DNA binding protein important for chromatin-targeted regulation of gene expression. It is best studied in the context of transvection-a mechanism of interallelic gene regulation involving paired chromosomes-and repression of the expression of white by Zeste mutants. Both of these functions depend on the DNA binding and self-association properties of Zeste, but the underlying structural basis remains unknown. Here we report the crystal structure of the DNA binding domain of Zeste in complex with a 19-bp DNA duplex containing the consensus recognition sequence motif. The structure reveals a helix-turn-helix Myb/homeodomain-like fold with the Zeste-specific insertion sequence forming a short helix and a long loop. Direct base contacts by the major groove binding helix principally account for the sequence-specific recognition, and backbone contacts via the Zeste-specific insertion are mainly responsible for the length requirement and the orientation of DNA. Our structural and biochemical characterizations of the DNA binding property of Zeste uncover an altered DNA binding modality of homeodomain-like proteins, and the structural information should facilitate the unraveling of the intricate mechanism of Zeste in regulation of gene expression.

High quality reference genome of drumstick tree (Moringa oleifera Lam.), a potential perennial crop

The drumstick tree (Moringa oleifera Lam.) is a perennial crop that has gained popularity in certain developing countries for its high-nutrition content and adaptability to arid and semi-arid environments. Here we report a high-quality draft genome sequence of M. oleifera. This assembly represents 91.78% of the estimated genome size and contains 19,465 protein-coding genes. Comparative genomic analysis between M. oleifera and related woody plant genomes helps clarify the general evolution of this species, while the identification of several species-specific gene families and positively selected genes in M. oleifera may help identify genes related to M. oleifera's high protein content, fast-growth, heat and stress tolerance. This reference genome greatly extends the basic research on M. oleifera, and may further promote applying genomics to enhanced breeding and improvement of M. oleifera.

Optimized position weight matrices in prediction of novel putative binding sites for transcription factors in the Drosophila melanogaster genome

Position weight matrices (PWMs) have become a tool of choice for the identification of transcription factor binding sites in DNA sequences. DNA-binding proteins often show degeneracy in their binding requirement and thus the overall binding specificity of many proteins is unknown and remains an active area of research. Although existing PWMs are more reliable predictors than consensus string matching, they generally result in a high number of false positive hits. Our previous study introduced a promising approach to PWM refinement in which known motifs are used to computationally mine putative binding sites directly from aligned promoter regions using composition of similar sites. In the present study, we extended this technique originally tested on single examples of transcription factors (TFs) and showed its capability to optimize PWM performance to predict new binding sites in the fruit fly genome. We propose refined PWMs in mono- and dinucleotide versions similarly computed for a large variety of transcription factors of Drosophila melanogaster. Along with the addition of many auxiliary sites the optimization includes variation of the PWM motif length, the binding sites location on the promoters and the PWM score threshold. To assess the predictive performance of the refined PWMs we compared them to conventional TRANSFAC and JASPAR sources. The results have been verified using performed tests and literature review. Overall, the refined PWMs containing putative sites derived from real promoter content processed using optimized parameters had better general accuracy than conventional PWMs.

Architectural and functional diversity of polycomb group response elements in Drosophila

Polycomb group response elements (PREs) play an essential role in gene regulation by the Polycomb group (PcG) repressor proteins in Drosophila. PREs are required for the recruitment and maintenance of repression by the PcG proteins. PREs are made up of binding sites for multiple DNA-binding proteins, but it is still unclear what combination(s) of binding sites is required for PRE activity. Here we compare the binding sites and activities of two closely linked yet separable PREs of the Drosophila engrailed (en) gene, PRE1 and PRE2. Both PRE1 and PRE2 contain binding sites for multiple PRE-DNA-binding proteins, but the number, arrangement, and spacing of the sites differs between the two PREs. These differences have functional consequences. Both PRE1 and PRE2 mediate pairing-sensitive silencing of mini-white, a functional assay for PcG repression; however, PRE1 requires two binding sites for Pleiohomeotic (Pho), whereas PRE2 requires only one Pho-binding site for this activity. Furthermore, for full pairing-sensitive silencing activity, PRE1 requires an AT-rich region not found in PRE2. These two PREs behave differently in a PRE embryonic and larval reporter construct inserted at an identical location in the genome. Our data illustrate the diversity of architecture and function of PREs.

Molecular genetic analysis of Suppressor 2 of zeste identifies key functional domains

The Su(z)2 complex contains Posterior sex combs (Psc) and Suppressor 2 of zeste [Su(z)2], two paralogous genes that likely arose by gene duplication. Psc encodes a Polycomb group protein that functions as a central component of the PRC1 complex, which maintains transcriptional repression of a wide array of genes. Although much is known about Psc, very little is known about Su(z)2, the analysis of which has been hampered by a dearth of alleles. We have generated new alleles of Su(z)2 and analyzed them at the genetic and molecular levels. Some of these alleles display negative complementation in that they cause lethality when heterozygous with the gain-of-function Su(z)2(1) allele but are hemizygous and, in some cases, homozygous viable. Interestingly, alleles of this class identify protein domains within Su(z)2 that are highly conserved in Psc and the mammalian Bmi-1 and Mel-18 proteins. We also find several domains of intrinsic disorder in the C-terminal regions of both Psc and Su(z)2 and suggest that these domains may contribute to the essential functions of both proteins.

In vivo identification of human small ubiquitin-like modifier polymerization sites by high accuracy mass spectrometry and an in vitro to in vivo strategy

The length and precise linkage of polyubiquitin chains is important for their biological activity. Although other ubiquitin-like proteins have the potential to form polymeric chains their identification in vivo is challenging and their functional role is unclear. Vertebrates express three small ubiquitin-like modifiers, SUMO-1, SUMO-2, and SUMO-3. Mature SUMO-2 and SUMO-3 are nearly identical and contain an internal consensus site for sumoylation that is missing in SUMO-1. Combining state-of-the-art mass spectrometry with an "in vitro to in vivo" strategy for post-translational modifications, we provide direct evidence that SUMO-1, SUMO-2, and SUMO-3 form mixed chains in cells via the internal consensus sites for sumoylation in SUMO-2 and SUMO-3. In vitro, the chain length of SUMO polymers could be influenced by changing the relative amounts of SUMO-1 and SUMO-2. The developed methodology is generic and can be adapted for the identification of other sumoylation sites in complex samples.

An Sp1/KLF binding site is important for the activity of a Polycomb group response element from the Drosophila engrailed gene

Polycomb-group response elements (PREs) are DNA elements through which the Polycomb-group (PcG) of transcriptional repressors act. Many of the PcG proteins are associated with two protein complexes that repress gene expression by modifying chromatin. Both of these protein complexes specifically associate with PREs in vivo, however, it is not known how they are recruited or held at the PRE. PREs are complex elements, made up of binding sites for many proteins. Our laboratory has been working to define all the sequences and DNA binding proteins required for the activity of a 181 bp PRE from the Drosophila engrailed gene. Here we show that one of the sites necessary for PRE activity, Site 2, can be bound by members of the Sp1/KLF family of zinc finger proteins. There are 10 Sp1/KLF family members in Drosophila, and nine of them bind to Site 2. We derive a consensus binding site for the Sp1/KLF Drosophila family members and show that this consensus sequence is present in most of the molecularly characterized PREs. These data suggest that one or more Sp1/KLF family members play a role in PRE function in Drosophila.

Profilin regulates the activity of p42POP, a novel Myb-related transcription factor

Profilins, regulators of cytoplasmic actin dynamics, also bind to several nuclear proteins but the significance of these interactions is mostly unclear. Here, we describe a novel Myb-related transcription factor, p42POP, as a new ligand for profilin and show that profilin regulates its activity. p42POP comprises a unique combination of domains and is widely expressed in mouse tissues. In contrast to many other Myb proteins, it contains only one functional tryptophan-cluster motif. This is followed by an acidic domain, a leucine zipper that mediates dimerization and functional nuclear import and export signals that can direct p42POP to either the nuclear or the cytoplasmic compartment. Binding to profilins is mediated by a proline-rich cluster. p42POP-profilin complexes can be precipitated from cell lysates. In transfected cells displaying p42POP in the nucleus, nuclear profilin is markedly increased. When p42POP is anchored at mitochondrial membranes, profilin is targeted to this location. Hence, in a cellular environment, p42POP and profilin are found in the same protein complex. In luciferase assays, p42POP acts as repressor and this activity is substantially reduced by profilins, indicating that profilin can regulate p42POP activity and is therefore involved in gene regulation.

Harbinger transposons and an ancient HARBI1 gene derived from a transposase

In this study we report main properties of Harbinger DNA transposons identified in protists, plants, insects, worms, and vertebrates. This is the first superfamily of eukaryotic DNA transposons where all autonomous transposons, even those that are hosted by species from different kingdoms, encode two proteins: a superfamily-specific transposase and a DNA-binding protein characterized by the presence of the conserved SANT/myb/trihelix motif. The last motif is known to be important for the DNA binding by different transcription regulators. Therefore, we suggest that this protein is necessary for coordinated expression of the Harbinger transposase. Although mammalian genomes are free of recognizable remnants of Harbingers, we identified a widely expressed HARBI1 gene encoding a 350-aa protein entirely derived from a Harbinger transposase some 450-500 million years ago. The HARBI1 proteins are conserved in humans, rats, mice, cows, pigs, chickens, frogs, and various bony fish, as well as other extremely important proteins, including RAG1 and RAG2. Conserved motifs detected in the Harbinger transposases are also well preserved in the HARBI1 proteins. Therefore, the HARBI1 proteins are expected to be nucleases important for functioning of bony vertebrates. We also found that the protein most similar to HARBI1 is encoded by an autonomous Harbinger 3_DR transposon that was transpositionally active in the zebrafish genome a few million years ago. Nonautonomous transposons derived from Harbinger3_DR are characterized by a striking preference for a 17-bp target site never seen previously in any other DNA transposon. Based on this observation, we suggest that the hypothetical HARBI1 nucleases are also characterized by a strong DNA-target specificity.

GAGA facilitates binding of Pleiohomeotic to a chromatinized Polycomb response element

Polycomb response elements (PREs) are chromosomal elements, typically comprising thousands of base pairs of poorly defined sequences that confer the maintenance of gene expression patterns by Polycomb group (PcG) repressors and trithorax group (trxG) activators. Genetic studies have indicated a synergistic requirement for the trxG protein GAGA and the PcG protein Pleiohomeotic (PHO) in silencing at several PREs. However, the molecular basis of this cooperation remains unknown. Here, using DNaseI footprinting analysis, we provide a high-resolution map of sites for the sequence- specific DNA-binding PcG protein PHO, trxG proteins GAGA and Zeste and the gap protein Hunchback (HB) on the 1.6 kb Ultrabithorax (Ubx) PRE. Although these binding elements are present throughout the PRE, they display clear patterns of clustering, suggestive of functional collaboration at the level of PRE binding. We found that while GAGA could efficiently bind to a chromatinized PRE, PHO alone was incapable of binding to chromatin. However, PHO binding to chromatin, but not naked DNA, was strongly facilitated by GAGA, indicating interdependence between GAGA and PHO already at the level of PRE binding. These results provide a biochemical explanation for the in vivo cooperation between GAGA and PHO and suggest that PRE function involves the integrated activities of genetically antagonistic trxG and PcG proteins.