The solution structure of a specific GAGA factor-DNA complex reveals a modular binding mode

C2H2 proteins: Evolutionary aspects of domain architecture and diversification

The largest group of transcription factors in higher eukaryotes are C2H2 proteins, which contain C2H2-type zinc finger domains that specifically bind to DNA. Few well-studied C2H2 proteins, however, demonstrate their key role in the control of gene expression and chromosome architecture. Here we review the features of the domain architecture of C2H2 proteins and the likely origin of C2H2 zinc fingers. A comprehensive investigation of proteomes for the presence of proteins with multiple clustered C2H2 domains has revealed a key difference between groups of organisms. Unlike plants, transcription factors in metazoans contain clusters of C2H2 domains typically separated by a linker with the TGEKP consensus sequence. The average size of C2H2 clusters varies substantially, even between genomes of higher metazoans, and with a tendency to increase in combination with SCAN, and especially KRAB domains, reflecting the increasing complexity of gene regulatory networks.

GAGA zinc finger transcription factor searches chromatin by 1D-3D facilitated diffusion

To elucidate how eukaryotic sequence-specific transcription factors (TFs) search for gene targets on chromatin, we used multi-color smFRET and single-particle imaging to track the diffusion of purified GAGA-Associated Factor (GAF) on DNA and nucleosomes. Monomeric GAF DNA-binding domain (DBD) bearing one zinc finger finds its cognate site by 1D or 3D diffusion on bare DNA and rapidly slides back-and-forth between naturally clustered motifs for seconds before escape. Multimeric, full-length GAF also finds clustered motifs on DNA by 1D-3D diffusion, but remains locked on target for longer periods. Nucleosome architecture effectively blocks GAF-DBD 1D-sliding into the histone core but favors retention of GAF-DBD when targeting solvent-exposed sites by 3D-diffusion. Despite the occlusive power of nucleosomes, 1D-3D facilitated diffusion enables GAF to effectively search for clustered cognate motifs in chromatin, providing a mechanism for navigation to nucleosome and nucleosome-free sites by a member of the largest TF family.

Molecular Evolution of RAMOSA1 (RA1) in Land Plants

RAMOSA1 (RA1) is a Cys2-His2-type (C2H2) zinc finger transcription factor that controls plant meristem fate and identity and has played an important role in maize domestication. Despite its importance, the origin of RA1 is unknown, and the evolution in plants is only partially understood. In this paper, we present a well-resolved phylogeny based on 73 amino acid sequences from 48 embryophyte species. The recovered tree topology indicates that, during grass evolution, RA1 arose from two consecutive SUPERMAN duplications, resulting in three distinct grass sequence lineages: RA1-like A, RA1-like B, and RA1; however, most of these copies have unknown functions. Our findings indicate that RA1 and RA1-like play roles in the nucleus despite lacking a traditional nuclear localization signal. Here, we report that copies diversified their coding region and, with it, their protein structure, suggesting different patterns of DNA binding and protein-protein interaction. In addition, each of the retained copies diversified regulatory elements along their promoter regions, indicating differences in their upstream regulation. Taken together, the evidence indicates that the RA1 and RA1-like gene families in grasses underwent subfunctionalization and neofunctionalization enabled by gene duplication.

NMR structure verifies the eponymous zinc finger domain of transcription factor ZNF750

ZNF750 is a nuclear transcription factor that activates skin differentiation and has tumor suppressor roles in several cancers. Unusually, ZNF750 has only a single zinc-finger (ZNF) domain, Z*, with an amino acid sequence that differs markedly from the CCHH family consensus. Because of its sequence differences Z* is classified as degenerate, presumed to have lost the ability to bind the zinc ion required for folding. AlphaFold predicts an irregular structure for Z* with low confidence. Low confidence predictions are often inferred to be intrinsically disordered regions of proteins, which would be the case if Z* did not bind Zn2+. We use NMR and CD spectroscopy to show that a 25-51 segment of ZNF750 corresponding to the Z* domain folds into a well-defined antiparallel ββα tertiary structure with a pM dissociation constant for Zn2+ and a thermal stability >80 °C. Of three alternative Zn2+ ligand sets, Z* uses a CCHC rather than the expected CCHH ligating motif. The switch in the last ligand maintains the folding topology and hydrophobic core of the classical ZNF motif. CCHC ZNFs are typically associated with protein-protein interactions, raising the possibility that ZNF750 interacts with DNA through other proteins rather than directly. The structure of Z* provides context for understanding the function of the domain and its cancer-associated mutations. We expect other ZNFs currently classified as degenerate could be CCHC-type structures like Z*.

WITHDRAWN: NMR structure verifies the eponymous zinc finger domain of transcription factor ZNF750

This article was initially published in the Journal of Structural Biology, instead of the Journal of Structural Biology: X, due to a publisher error. We regret the inconvenience. The link to the article published in Journal of Structural Biology: X is presented below: https://www.sciencedirect.com/science/article/pii/S2590152423000090. The full Elsevier Policy on Article Withdrawal can be found at https://www.elsevier.com/about/policies/article-withdrawal.

Kinetic principles underlying pioneer function of GAGA transcription factor in live cells

How pioneer factors interface with chromatin to promote accessibility for transcription control is poorly understood in vivo. Here, we directly visualize chromatin association by the prototypical GAGA pioneer factor (GAF) in live Drosophila hemocytes. Single-particle tracking reveals that most GAF is chromatin bound, with a stable-binding fraction showing nucleosome-like confinement residing on chromatin for more than 2 min, far longer than the dynamic range of most transcription factors. These kinetic properties require the full complement of GAF's DNA-binding, multimerization and intrinsically disordered domains, and are autonomous from recruited chromatin remodelers NURF and PBAP, whose activities primarily benefit GAF's neighbors such as Heat Shock Factor. Evaluation of GAF kinetics together with its endogenous abundance indicates that, despite on-off dynamics, GAF constitutively and fully occupies major chromatin targets, thereby providing a temporal mechanism that sustains open chromatin for transcriptional responses to homeostatic, environmental and developmental signals.

Zinc finger structure determination by NMR: Why zinc fingers can be a handful

Zinc fingers can be loosely defined as protein domains containing one or more tetrahedrally-co-ordinated zinc ions whose role is to stabilise the structure rather than to be involved in enzymatic chemistry; such zinc ions are often referred to as "structural zincs". Although structural zincs can occur in proteins of any size, they assume particular significance for very small protein domains, where they are often essential for maintaining a folded state. Such small structures, that sometimes have only marginal stability, can present particular difficulties in terms of sample preparation, handling and structure determination, and early on they gained a reputation for being resistant to crystallisation. As a result, NMR has played a more prominent role in structural studies of zinc finger proteins than it has for many other types of proteins. This review will present an overview of the particular issues that arise for structure determination of zinc fingers by NMR, and ways in which these may be addressed.

[Lack of GAGA protein in Trl mutants causes massive cell death in Drosophila spermatogenesis and oogenesis]

Белок дрозофилы GAGA (GAF) является фактором эпигенетической регуляции транскрипции большой группы генов с широким разнообразием клеточных функций. GAF кодируется геном Trithorax-like (Trl), который экспрессируется в различных органах и тканях на всех стадиях онтогенеза дрозофилы. Мутации этого гена вызывают множественные нарушения развития. В предыдущих работах мы показали, что этот белок необходим для развития половой системы как самцов, так и самок дрозофилы. Снижение экспрессии гена Trl приводило ко множественным нарушениям спермато- и оогенеза. Одно из значительных нарушений было связано с массовой деградацией и потерей клеток зародышевого пути, что позволило предположить, что этот белок вовлечен в регуляцию клеточной гибели. В представленной работе мы провели более детальное цитологическое исследование, чтобы определить, какой тип гибели клеток зародышевого пути характерен для Trl-мутантов, и происходят ли нарушения или изменения этого процесса по сравнению с нормой. Полученные результаты показали, что недостаток белка GAF вызывает массовую гибель клеток зародышевого пути как у самок, так и самцов дрозофилы, но проявляется эта гибель в зависимости от пола по-разному. У самок, мутантных по гену Trl, фенотипически этот процесс не отличается от нормы и в гибнущих яйцевых камерах выявлены признаки апоптоза и аутофагии клеток зародышевого пути. У самцов, мутантных по гену Trl, в отличие от самок, не обнаружены признаки апоптоза. У самцов мутации Trl индуцируют массовую гибель клеток через аутофагию, что не характерно для сперматогенеза дрозофилы и не описано ранее ни в норме, ни у мутаций по другим генам. Таким образом, недостаток GAF у мутантов Trl приводит к усилению апоптотической и аутофагической гибели клеток зародышевого пути. Эктопическая клеточная гибель и атрофия зародышевой линии, вероятно, связаны с нарушением экспрессии генов-мишеней GAGAфактора, среди которых есть гены, регулирующие как апоптоз, так и аутофагию.

Liquid condensation of reprogramming factor KLF4 with DNA provides a mechanism for chromatin organization

Expression of a few master transcription factors can reprogram the epigenetic landscape and three-dimensional chromatin topology of differentiated cells and achieve pluripotency. During reprogramming, thousands of long-range chromatin contacts are altered, and changes in promoter association with enhancers dramatically influence transcription. Molecular participants at these sites have been identified, but how this re-organization might be orchestrated is not known. Biomolecular condensation is implicated in subcellular organization, including the recruitment of RNA polymerase in transcriptional activation. Here, we show that reprogramming factor KLF4 undergoes biomolecular condensation even in the absence of its intrinsically disordered region. Liquid–liquid condensation of the isolated KLF4 DNA binding domain with a DNA fragment from the NANOG proximal promoter is enhanced by CpG methylation of a KLF4 cognate binding site. We propose KLF4-mediated condensation as one mechanism for selectively organizing and re-organizing the genome based on the local sequence and epigenetic state.

KLF4, OCT4, SOX2 and MYC cooperate to reorganize chromatin during somatic cell reprogramming. Here the authors show that KLF4 forms a liquid-like biomolecular condensate that recruits OCT4 and SOX2, and that condensation of the isolated KLF4 DNA binding domain with DNA is enhanced by CpG methylation

Cell-free genomics reveal intrinsic, cooperative and competitive determinants of chromatin interactions

Metazoan transcription factors distinguish their response elements from a large excess of similar sequences. We explored underlying principles of DNA shape read-out and factor cooperativity in chromatin using a unique experimental system. We reconstituted chromatin on Drosophila genomes in extracts of preblastoderm embryos, mimicking the naïve state of the zygotic genome prior to developmental transcription activation. We then compared the intrinsic binding specificities of three recombinant transcription factors, alone and in combination, with GA-rich recognition sequences genome-wide. For MSL2, all functional elements reside on the X chromosome, allowing to distinguish physiological elements from non-functional 'decoy' sites. The physiological binding profile of MSL2 is approximated through interaction with other factors: cooperativity with CLAMP and competition with GAF, which sculpts the profile by occluding non-functional sites. An extended DNA shape signature is differentially read out in chromatin. Our results reveal novel aspects of target selection in a complex chromatin environment.

Kinetoplastid kinetochore proteins KKT2 and KKT3 have unique centromere localization domains

The kinetochore is the macromolecular protein complex that assembles onto centromeric DNA and binds spindle microtubules. Evolutionarily divergent kinetoplastids have an unconventional set of kinetochore proteins. It remains unknown how kinetochores assemble at centromeres in these organisms. Here, we characterize KKT2 and KKT3 in the kinetoplastid parasite Trypanosoma brucei. In addition to the N-terminal kinase domain and C-terminal divergent polo boxes, these proteins have a central domain of unknown function. We show that KKT2 and KKT3 are important for the localization of several kinetochore proteins and that their central domains are sufficient for centromere localization. Crystal structures of the KKT2 central domain from two divergent kinetoplastids reveal a unique zinc-binding domain (termed the CL domain for centromere localization), which promotes its kinetochore localization in T. brucei. Mutations in the equivalent domain in KKT3 abolish its kinetochore localization and function. Our work shows that the unique central domains play a critical role in mediating the centromere localization of KKT2 and KKT3.

GAGA factor: a multifunctional pioneering chromatin protein

The Drosophila GAGA factor (GAF) is a multifunctional protein implicated in nucleosome organization and remodeling, activation and repression of gene expression, long distance enhancer-promoter communication, higher order chromosome structure, and mitosis. This broad range of activities poses questions about how a single protein can perform so many seemingly different and unrelated functions. Current studies argue that GAF acts as a "pioneer" factor, generating nucleosome-free regions of chromatin for different classes of regulatory elements. The removal of nucleosomes from regulatory elements in turn enables other factors to bind to these elements and carry out their specialized functions. Consistent with this view, GAF associates with a collection of chromatin remodelers and also interacts with proteins implicated in different regulatory functions. In this review, we summarize the known activities of GAF and the functions of its protein partners.

Accelerated evolution of a minimal 63-amino acid dual transcription factor

Transcription factors control gene expression in all life. This raises the question of what is the smallest protein that can support such activity. In nature, Cro from bacteriophage λ is one of the smallest known repressors (66 amino acids), and activators are typically much larger (e.g., λ cI, 237 amino acids). Previous efforts to engineer a minimal activator from λ Cro resulted in no activity in vivo in cells. In this study, we show that directed evolution results in a new Cro activator-repressor that functions as efficiently as λ cI in vivo. To achieve this, we develop phagemid-assisted continuous evolution (PACEmid). We find that a peptide as small as 63 amino acids functions efficiently as an activator and/or repressor. To our knowledge, this is the smallest protein activator that enables polymerase recruitment, highlighting the capacity of transcription factors to evolve from very short peptide sequences.

Basal hsp70 expression levels do not explain adaptive variation of the warm- and cold-climate O3 + 4 + 7 and OST gene arrangements of Drosophila subobscura

Background

Drosophila subobscura exhibits a rich inversion polymorphism, with some adaptive inversions showing repeatable spatiotemporal patterns in frequencies related to temperature. Previous studies reported increased basal HSP70 protein levels in homokaryotypic strains for a warm-climate arrangement compared to a cold-climate one. These findings do not match the similar hsp70 genomic organization between arrangements, where gene expression levels are expected to be similar. In order to test this hypothesis and understand the molecular basis for hsp70 expression, we compared basal hsp70 mRNA levels in males and females, and analysed the 5′ and 3′ regulatory regions of hsp70 genes in warm- and cold-climate isochromosomal O_{3 + 4 + 7} and O_ST lines of D. subobscura.

Results

We observed comparable mRNA levels between the two arrangements and a sex-biased hsp70 gene expression. The number of heat-shock elements (HSEs) and GAGA sites on the promoters were identical amongst the O_ST and O_{3 + 4 + 7} lines analysed. This is also true for 3′ AU-rich elements where most A and B copies of hsp70 have, respectively, two and one element in both arrangements. Beyond the regulatory elements, the only notable difference between both arrangements is the presence in 3′ UTR of a 14 bp additional fragment after the stop codon in the hsp70A copy in five O_{3 + 4 + 7} lines, which was not found in any of the six O_ST lines.

Conclusions

The equivalent hsp70 mRNA amounts in O_ST and O_{3 + 4 + 7} arrangements provide the first evidence of a parallelism between gene expression and genetic organization in D. subobscura lines having these arrangements. This is reinforced by the lack of important differential features in the number and structure of regulatory elements between both arrangements, despite the genetic differentiation observed when the complete 5′ and 3′ regulatory regions were considered. Therefore, the basal levels of hsp70 mRNA cannot account, in principle, for the adaptive variation of the two arrangements studied. Consequently, further studies are necessary to understand the intricate molecular mechanisms of hsp70 gene regulation in D. subobscura.

CoCUN, a Novel Ubiquitin Binding Domain Identified in N4BP1

Ubiquitin binding domains (UBDs) are modular elements that bind non-covalently to ubiquitin and act as downstream effectors and amplifiers of the ubiquitination signal. With few exceptions, UBDs recognize the hydrophobic path centered on Ile44, including residues Leu8, Ile44, His68, and Val70. A variety of different orientations, which can be attributed to specific contacts between each UBD and surface residues surrounding the hydrophobic patch, specify how each class of UBD specifically contacts ubiquitin. Here, we describe the structural model of a novel ubiquitin-binding domain that we identified in NEDD4 binding protein 1 (N4BP1). By performing protein sequence analysis, mutagenesis, and nuclear magnetic resonance (NMR) spectroscopy of the ¹⁵N isotopically labeled protein, we demonstrate that a Phe-Pro motif in N4BP1 recognizes the canonical hydrophobic patch of ubiquitin. This recognition mode resembles the molecular mechanism evolved in the coupling of ubiquitin conjugation to endoplasmic-reticulum (ER) degradation (CUE) domain family, where an invariant proline, usually following a phenylalanine, is required for ubiquitin binding. Interestingly, this novel UBD, which is not evolutionary related to CUE domains, shares a 40% identity and 47% similarity with cullin binding domain associating with NEDD8 (CUBAN), a protein module that also recognizes the ubiquitin-like NEDD8. Based on these features, we dubbed the region spanning the C-terminal 50 residues of N4BP1 the CoCUN domain, for Cousin of CUBAN. By performing circular dichroism and ¹⁵N NMR chemical shift perturbation of N4BP1 in complex with ubiquitin, we demonstrate that the CoCUN domain lacks the NEDD8 binding properties observed in CUBAN. We also show that, in addition to mediating the interaction with ubiquitin and ubiquitinated substrates, both CUBAN and CoCUN are poly-ubiquitinated in cells. The structural and the functional characterization of this novel UBD can contribute to a deeper understanding of the molecular mechanisms governing N4BP1 function, providing at the same time a valuable tool for clarifying how the discrimination between ubiquitin and the highly related NEDD8 is achieved.

Specific interaction of zinc finger protein Com with RNA and the crystal structure of a self-complementary RNA duplex recognized by Com

The bacteriophage Mu Com is a small zinc finger protein that binds to its cognate mom mRNA and activates its translation. The Mom protein, in turn, elicits a chemical modification (momification) of the bacteriophage genome, rendering the DNA resistant to cleavage by bacterial restriction endonucleases, and thereby protecting it from defense mechanisms of the host. We examined the basis of specificity in Com-RNA interactions by in vitro selection and probing of RNA structure. We demonstrated that Com recognizes a sequence motif within a hairpin-loop structure of its target RNA. Our data support the model of Com interaction with mom mRNA, in which Com binds to the short hairpin structure proximal to the so-called translation inhibition structure. We also observed that Com binds its target motif weakly if it is within an RNA duplex. These results suggest that the RNA structure, in addition to its sequence, is crucial for Com to recognize its target and that RNA conformational changes may constitute another level of Mom regulation. We determined a crystal structure of a Com binding site variant designed to form an RNA duplex preferentially. Our crystal model forms a 19-mer self-complementary double helix composed of the canonical and non-canonical base pairs. The helical parameters of crystalized RNA indicate why Com may bind it more weakly than a monomeric hairpin form.

C2H2-Type Zinc Finger Proteins: Evolutionarily Old and New Partners of the Nuclear Hormone Receptors

Nuclear hormone receptors (NRs) are evolutionarily conserved ligand-dependent transcription factors. They are essential for human life, mediating the actions of lipophilic molecules, such as steroid hormones and metabolites of fatty acid, cholesterol, and external toxic compounds. The C2H2-type zinc finger proteins (ZNFs) form the largest family of the transcription factors in humans and are characterized by multiple, tandemly arranged zinc fingers. Many of the C2H2-type ZNFs are conserved throughout evolution, suggesting their involvement in preserved biological activities, such as general transcriptional regulation and development/differentiation of organs/tissues observed in the early embryonic phase. However, some C2H2-type ZNFs, such as those with the Krüppel-associated box (KRAB) domain, appeared relatively late in evolution and have significantly increased family members in mammals including humans, possibly modulating their complicated transcriptional network and/or supporting the morphological development/functions specific to them. Such evolutional characteristics of the C2H2-type ZNFs indicate that these molecules influence the NR functions conserved through evolution, whereas some also adjust them to meet with specific needs of higher organisms. We review the interaction between NRs and C2H2-type ZNFs by focusing on some of the latter molecules.

DNA-binding miniproteins based on zinc fingers. Assessment of the interaction using nanopores

We report a synthetic miniprotein that combines zinc finger modules of the transcription factor GAGA with the AT-hook peptide. This designed chimera binds to extended DNA sites with high affinity and selectivity, as shown by nanopore force spectroscopy.

Obtaining artificial proteins that mimic the DNA binding properties of natural transcription factors could open new ways of manipulating gene expression at will. In this context it is particularly interesting to develop simple synthetic systems. Inspired by the modularity of natural transcription factors, we have designed synthetic miniproteins that combine the zinc finger module of the transcription factor GAGA and AT-hook peptide domains. These constructs are capable of binding to composite DNA sequences of up to 14 base pairs with high affinity and good selectivity. In particular, we have synthesized three different chimeras and characterized their DNA binding properties by electrophoresis and fluorescence anisotropy. We have also used, for the first time in the study of peptide-based DNA binders, nanopore force spectroscopy to obtain further data on the DNA interaction.

Differential Occupancy of Two GA-Binding Proteins Promotes Targeting of the Drosophila Dosage Compensation Complex to the Male X Chromosome

Little is known about how variation in sequence composition alters transcription factor occupancy to precisely recruit large transcription complexes. A key model for understanding how transcription complexes are targeted is the Drosophila dosage compensation system in which the male-specific lethal (MSL) transcription complex specifically identifies and regulates the male X chromosome. The chromatin-linked adaptor for MSL proteins (CLAMP) zinc-finger protein targets MSL to the X chromosome but also binds to GA-rich sequence elements throughout the genome. Furthermore, the GAGA-associated factor (GAF) transcription factor also recognizes GA-rich sequences but does not associate with the MSL complex. Here, we demonstrate that MSL complex recruitment sites are optimal CLAMP targets. Specificity for CLAMP binding versus GAF binding is driven by variability in sequence composition within similar GA-rich motifs. Therefore, variation within seemingly similar cis elements drives the context-specific targeting of a large transcription complex.

Vertebrate GAF/ThPOK: emerging functions in chromatin architecture and transcriptional regulation

GAGA factor of Drosophila melanogaster (DmGAF) is a multifaceted transcription factor with diverse roles in chromatin regulation. Recently, ThPOK/c-Krox was identified as its vertebrate homologue (vGAF), which has a basic domain structure similar to DmGAF and is decorated with a number of post-translationally modified residues. In vertebrate genomes, vGAF associates with purine-rich GAGA sequences and performs diverse chromatin-mediated functions, viz., gene activation, repression and enhancer blocking. Expansion of regulatory chromatin proteins with the acquisition of PTMs appears to be the general trend that facilitated the evolution of complexity in vertebrates. Here, we compare the structural and functional features of vGAF with those of DmGAF and also assess the possible functional redundancy among paralogues of vGAF. We also discuss the underlying mechanisms which aid in the diverse and context-dependent functions of this protein.

The conserved basic residues and the charged amino acid residues at the α-helix of the zinc finger motif regulate the nuclear transport activity of triple C2H2 zinc finger proteins

Zinc finger (ZF) motifs on proteins are frequently recognized as a structure for DNA binding. Accumulated reports indicate that ZF motifs contain nuclear localization signal (NLS) to facilitate the transport of ZF proteins into nucleus. We investigated the critical factors that facilitate the nuclear transport of triple C2H2 ZF proteins. Three conserved basic residues (hot spots) were identified among the ZF sequences of triple C2H2 ZF proteins that reportedly have NLS function. Additional basic residues can be found on the α-helix of the ZFs. Using the ZF domain (ZFD) of Egr-1 as a template, various mutants were constructed and expressed in cells. The nuclear transport activity of various mutants was estimated by analyzing the proportion of protein localized in the nucleus. Mutation at any hot spot of the Egr-1 ZFs reduced the nuclear transport activity. Changes of the basic residues at the α-helical region of the second ZF (ZF2) of the Egr-1 ZFD abolished the NLS activity. However, this activity can be restored by substituting the acidic residues at the homologous positions of ZF1 or ZF3 with basic residues. The restored activity dropped again when the hot spots at ZF1 or the basic residues in the α-helix of ZF3 were mutated. The variations in nuclear transport activity are linked directly to the binding activity of the ZF proteins with importins. This study was extended to other triple C2H2 ZF proteins. SP1 and KLF families, similar to Egr-1, have charged amino acid residues at the second (α2) and the third (α3) positions of the α-helix. Replacing the amino acids at α2 and α3 with acidic residues reduced the NLS activity of the SP1 and KLF6 ZFD. The reduced activity can be restored by substituting the α3 with histidine at any SP1 and KLF6 ZFD. The results show again the interchangeable role of ZFs and charge residues in the α-helix in regulating the NLS activity of triple C2H2 ZF proteins.

The Drosophila CLAMP protein associates with diverse proteins on chromatin

Gaining new insights into gene regulation involves an in-depth understanding of protein-protein interactions on chromatin. A powerful model for studying mechanisms of gene regulation is dosage compensation, a process that targets the X-chromosome to equalize gene expression between XY males and XX females. We previously identified a zinc finger protein in Drosophila melanogaster that plays a sex-specific role in targeting the Male-specific lethal (MSL) dosage compensation complex to the male X-chromosome, called the Chromatin-Linked Adapter for MSL Proteins (CLAMP). More recently, we established that CLAMP has non-sex-specific roles as an essential protein that regulates chromatin accessibility at promoters genome-wide. To identify associations between CLAMP and other factors in both male and female cells, we used two complementary mass spectrometry approaches. We demonstrate that CLAMP associates with the transcriptional regulator complex Negative Elongation Factor (NELF) in both sexes and determine that CLAMP reduces NELF recruitment to several target genes. In sum, we have identified many new CLAMP-associated factors and provide a resource for further study of this little understood essential protein.

Enhanced chromatin accessibility of the dosage compensated Drosophila male X-chromosome requires the CLAMP zinc finger protein

The essential process of dosage compensation is required to equalize gene expression of X-chromosome genes between males (XY) and females (XX). In Drosophila, the conserved Male-specific lethal (MSL) histone acetyltransferase complex mediates dosage compensation by increasing transcript levels from genes on the single male X-chromosome approximately two-fold. Consistent with its increased levels of transcription, the male X-chromosome has enhanced chromatin accessibility, distinguishing it from the autosomes. Here, we demonstrate that the non-sex-specific CLAMP (Chromatin-linked adaptor for MSL proteins) zinc finger protein that recognizes GA-rich sequences genome-wide promotes the specialized chromatin environment on the male X-chromosome and can act over long genomic distances (~14 kb). Although MSL complex is required for increasing transcript levels of X-linked genes, it is not required for enhancing global male X-chromosome chromatin accessibility, and instead works cooperatively with CLAMP to facilitate an accessible chromatin configuration at its sites of highest occupancy. Furthermore, CLAMP regulates chromatin structure at strong MSL complex binding sites through promoting recruitment of the Nucleosome Remodeling Factor (NURF) complex. In contrast to the X-chromosome, CLAMP regulates chromatin and gene expression on autosomes through a distinct mechanism that does not involve NURF recruitment. Overall, our results support a model where synergy between a non-sex-specific transcription factor (CLAMP) and a sex-specific cofactor (MSL) creates a specialized chromatin domain on the male X-chromosome.

Transcriptomic data from panarthropods shed new light on the evolution of insulator binding proteins in insects : Insect insulator proteins

Background

Body plan development in multi-cellular organisms is largely determined by homeotic genes. Expression of homeotic genes, in turn, is partially regulated by insulator binding proteins (IBPs). While only a few enhancer blocking IBPs have been identified in vertebrates, the common fruit fly Drosophila melanogaster harbors at least twelve different enhancer blocking IBPs. We screened recently compiled insect transcriptomes from the 1KITE project and genomic and transcriptomic data from public databases, aiming to trace the origin of IBPs in insects and other arthropods.

Results

Our study shows that the last common ancestor of insects (Hexapoda) already possessed a substantial number of IBPs. Specifically, of the known twelve insect IBPs, at least three (i.e., CP190, Su(Hw), and CTCF) already existed prior to the evolution of insects. Furthermore we found GAF orthologs in early branching insect orders, including Zygentoma (silverfish and firebrats) and Diplura (two-pronged bristletails). Mod(mdg4) is most likely a derived feature of Neoptera, while Pita is likely an evolutionary novelty of holometabolous insects. Zw5 appears to be restricted to schizophoran flies, whereas BEAF-32, ZIPIC and the Elba complex, are probably unique to the genus Drosophila. Selection models indicate that insect IBPs evolved under neutral or purifying selection.

Conclusions

Our results suggest that a substantial number of IBPs either pre-date the evolution of insects or evolved early during insect evolution. This suggests an evolutionary history of insulator binding proteins in insects different to that previously thought. Moreover, our study demonstrates the versatility of the 1KITE transcriptomic data for comparative analyses in insects and other arthropods.

Electronic supplementary material

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

Transcription factors GAF and HSF act at distinct regulatory steps to modulate stress-induced gene activation

The coordinated regulation of gene expression at the transcriptional level is fundamental to development and homeostasis. Inducible systems are invaluable when studying transcription because the regulatory process can be triggered instantaneously, allowing the tracking of ordered mechanistic events. Here, we use precision run-on sequencing (PRO-seq) to examine the genome-wide heat shock (HS) response in Drosophila and the function of two key transcription factors on the immediate transcription activation or repression of all genes regulated by HS. We identify the primary HS response genes and the rate-limiting steps in the transcription cycle that GAGA-associated factor (GAF) and HS factor (HSF) regulate. We demonstrate that GAF acts upstream of promoter-proximally paused RNA polymerase II (Pol II) formation (likely at the step of chromatin opening) and that GAF-facilitated Pol II pausing is critical for HS activation. In contrast, HSF is dispensable for establishing or maintaining Pol II pausing but is critical for the release of paused Pol II into the gene body at a subset of highly activated genes. Additionally, HSF has no detectable role in the rapid HS repression of thousands of genes.

A designed DNA binding motif that recognizes extended sites and spans two adjacent major grooves

We report the rational design of a DNA-binding peptide construct composed of the DNA-contacting regions of two transcription factors (GCN4 and GAGA) linked through an AT-hook DNA anchor. The resulting chimera, which represents a new, non-natural DNA binding motif, binds with high affinity and selectivity to a long composite sequence of 13 base pairs (TCAT-AATT-GAGAG).

Solution structure of the 13th C2H2 Zinc Finger of Miz-1

Miz-1 is a BTB/POZ transcription factor that contains 13C2H2 Zinc Finger domains (ZF). Miz-1 transactivates and represses the transcription of a myriad of genes involved in many aspects of the biology of the cell. The detailed molecular interactions through which Miz-1 controls transcription, including its specific DNA binding via its ZF domains, remain to be understood and documented. In our effort to shed light into the structural biology of Miz-1, we have undertaken the determination of the structure of all its ZF and the characterization of their interactions with cognate DNA. The structure of ZF 1 to 10 have already been solved and characterized. Here, we present the structure of the synthetic Miz-1 ZF13 determined by 2D (1)H-(1)H NMR.

The prokaryotic zinc-finger: structure, function and comparison with the eukaryotic counterpart

Classical zinc finger (ZF) domains were thought to be confined to the eukaryotic kingdom until the transcriptional regulator Ros protein was identified in Agrobacterium tumefaciens. The Ros Cys2 His2 ZF binds DNA in a peculiar mode and folds in a domain significantly larger than its eukaryotic counterpart consisting of 58 amino acids (the 9-66 region) arranged in a βββαα topology, and stabilized by a conserved, extensive, 15-residue hydrophobic core. The prokaryotic ZF domain, then, shows some intriguing new features that make it interestingly different from its eukaryotic counterpart. This review will focus on the prokaryotic ZFs, summarizing and discussing differences and analogies with the eukaryotic domains and providing important insights into their structure/function relationships.

The AT-Hook motif as a versatile minor groove anchor for promoting DNA binding of transcription factor fragments

We report the development of chimeric DNA binding peptides comprising a DNA binding fragment of natural transcription factors (the basic region of a bZIP protein or a monomeric zinc finger module) and an AT-Hook peptide motif. The resulting peptide conjugates display high DNA affinity and excellent sequence selectivity. Furthermore, the AT-Hook motif also favors the cell internalization of the conjugates.

Spatial proximity of homologous alleles and long noncoding RNAs regulate a switch in allelic gene expression

Physiological processes rely on the regulation of total mRNA levels in a cell. In diploid organisms, the transcriptional activation of one or both alleles of a gene may involve trans-allelic interactions that provide a tight spatial and temporal level of gene expression regulation. The mechanisms underlying such interactions still remain poorly understood. Here, we demonstrate that lipopolysaccharide stimulation of murine macrophages rapidly resulted in the actin-mediated and transient homologous spatial proximity of Tnfα alleles, which was necessary for the mono- to biallelic switch in gene expression. We identified two new complementary long noncoding RNAs transcribed from the TNFα locus and showed that their knockdown had opposite effects in Tnfα spatial proximity and allelic expression. Moreover, the observed spatial proximity of Tnfα alleles depended on pyruvate kinase muscle isoform 2 (PKM2) and T-helper-inducing POZ-Krüppel-like factor (ThPOK). This study suggests a role for lncRNAs in the regulation of somatic homologous spatial proximity and allelic expression control necessary for fine-tuning mammalian immune responses.

GAGA factor maintains nucleosome-free regions and has a role in RNA polymerase II recruitment to promoters

Previous studies have shown that GAGA Factor (GAF) is enriched on promoters with paused RNA Polymerase II (Pol II), but its genome-wide function and mechanism of action remain largely uncharacterized. We assayed the levels of transcriptionally-engaged polymerase using global run-on sequencing (GRO-seq) in control and GAF-RNAi Drosophila S2 cells and found promoter-proximal polymerase was significantly reduced on a large subset of paused promoters where GAF occupancy was reduced by knock down. These promoters show a dramatic increase in nucleosome occupancy upon GAF depletion. These results, in conjunction with previous studies showing that GAF directly interacts with nucleosome remodelers, strongly support a model where GAF directs nucleosome displacement at the promoter and thereby allows the entry Pol II to the promoter and pause sites. This action of GAF on nucleosomes is at least partially independent of paused Pol II because intergenic GAF binding sites with little or no Pol II also show GAF-dependent nucleosome displacement. In addition, the insulator factor BEAF, the BEAF-interacting protein Chriz, and the transcription factor M1BP are strikingly enriched on those GAF-associated genes where pausing is unaffected by knock down, suggesting insulators or the alternative promoter-associated factor M1BP protect a subset of GAF-bound paused genes from GAF knock-down effects. Thus, GAF binding at promoters can lead to the local displacement of nucleosomes, but this activity can be restricted or compensated for when insulator protein or M1BP complexes also reside at GAF bound promoters.

Transcriptional regulation is critical for proper gene expression in response to environmental changes and developmental programs. Eukaryotes have evolved multiple mechanisms by which transcription factors regulate transcription. One mechanism is the reorganization of chromatin to allow Pol II recruitment. Another is the release of promoter-proximal paused Pol II, where Pol II transcription that is halted 20–60 bases downstream of the transcription start site (TSS) is allowed to enter into productive elongation through the gene body. The Drosophila transcription factor GAF binds to genes that undergo pausing and interacts with nucleosome remodelers and the pausing factor NELF. Thus, GAF can regulate multiple points necessary for transcription, but its mechanistic role is not fully understood genome-wide. We depleted GAF from cells and examined the genome-wide changes in Pol II and nucleosome distributions across genes. We found that GAF depletion reduces polymerase density at genes where GAF binds just upstream of the TSS, and results in nucleosomes moving into the promoter region. Our results show that GAF is important for maintaining the promoter accessibility, allowing Pol II to be recruited to promoters and enter the pause sites downstream of the TSS. Thus, GAF is critical for providing the chromatin environment necessary for the proper control of gene expression.

Structure of the cyclic peptide [W8S]contryphan Vn: effect of the tryptophan/serine substitution on trans-cis proline isomerization

The structural characterization of [W8S]contryphan Vn, an analogue of Contryphan Vn with tryptophan 8 substituted with a serine residue (W8S), was performed by NMR spectroscopy, molecular dynamics simulations and fluorescence spectroscopy. Contryphan Vn, a bioactive cyclic peptide from the venom of the cone snail Conus ventricosus, contains an S-S bridge between two cysteines and a D-tryptophan. Like other Contryphans, [W8S]contryphan Vn has proline 7 isomerized trans, while the proline 4 has nearly equivalent populations of cis and trans configurations. The thermodynamic and kinetic parameters of the trans-cis isomerization of proline 4 were measured. The isomers of [W8S]contryphan Vn with proline 4 in cis and trans show structural differences. The absence of the salt bridge between the same Asp2 and Lys6, present in Contryphan Vn, may be attributed to the lack of the hydrophobic side chain of Trp8 where it likely protects the electrostatic interactions. These results may contribute to identifying, in these cyclic peptides, the structural determinants of the mechanism of proline trans-cis isomerization, this being also an important step in protein folding.

Towards understanding the molecular recognition process in prokaryotic zinc-finger domain

Eukaryotic Cys2His2 zinc finger domain is one of the most common and important structural motifs involved in protein-DNA interaction. The recognition motif is characterized by the tetrahedral coordination of a zinc ion by conserved cysteine and histidine residues. We have characterized the prokaryotic Cys2His2 zinc finger motif, included in the DNA binding region (Ros87) of Ros protein from Agrobacterium tumefaciens, demonstrating that, although possessing a similar zinc coordination sphere, this domain presents significant differences from its eukaryotic counterpart. Furthermore, basic residues flanking the zinc binding region on either side have been demonstrated, by Electrophoretic Mobility Shift Assay (EMSA) experiments, to be essential for Ros DNA binding. In spite of this wealth of knowledge, the structural details of the mechanism through which the prokaryotic zinc fingers recognize their target genes are still unclear. Here, to gain insights into the molecular DNA recognition process of prokaryotic zinc finger domains we applied a strategy in which we performed molecular docking studies using a combination of Nuclear Magnetic Resonance (NMR) and Molecular Dynamics (MD) simulations data. The results demonstrate that the MD ensemble provides a reasonable picture of Ros87 backbone dynamics in solution. The Ros87-DNA model indicates that the interaction involves the first two residue of the first α-helix, and several residues located in the basic regions flanking the zinc finger domain. Interestingly, the prokaryotic zinc finger domain, mainly with the C-terminal tail that is wrapped around the DNA, binds a more extended recognition site than the eukaryotic counterpart. Our analysis demonstrates that the introduction of the protein flexibility in docking studies can improve, in terms of accuracy, the quality of the obtained models and could be particularly useful for protein showing high conformational heterogeneity as well as for computational drug design applications.

Diversification of transcription factor paralogs via noncanonical modularity in C2H2 zinc finger DNA binding

A major challenge in obtaining a full molecular description of evolutionary adaptation is to characterize how transcription factor (TF) DNA-binding specificity can change. To identify mechanisms of TF diversification, we performed detailed comparisons of yeast C2H2 ZF proteins with identical canonical recognition residues that are expected to bind the same DNA sequences. Unexpectedly, we found that ZF proteins can adapt to recognize new binding sites in a modular fashion whereby binding to common core sites remains unaffected. We identified two distinct mechanisms, conserved across multiple Ascomycota species, by which this molecular adaptation occurred. Our results suggest a route for TF evolution that alleviates negative pleiotropic effects by modularly gaining new binding sites. These findings expand our current understanding of ZF DNA binding and provide evidence for paralogous ZFs utilizing alternate modes of DNA binding to recognize unique sets of noncanonical binding sites.

Regulatory modules controlling maize inflorescence architecture

Genetic control of branching is a primary determinant of yield, regulating seed number and harvesting ability, yet little is known about the molecular networks that shape grain-bearing inflorescences of cereal crops. Here, we used the maize (Zea mays) inflorescence to investigate gene networks that modulate determinacy, specifically the decision to allow branch growth. We characterized developmental transitions by associating spatiotemporal expression profiles with morphological changes resulting from genetic perturbations that disrupt steps in a pathway controlling branching. Developmental dynamics of genes targeted in vivo by the transcription factor RAMOSA1, a key regulator of determinacy, revealed potential mechanisms for repressing branches in distinct stem cell populations, including interactions with KNOTTED1, a master regulator of stem cell maintenance. Our results uncover discrete developmental modules that function in determining grass-specific morphology and provide a basis for targeted crop improvement and translation to other cereal crops with comparable inflorescence architectures.

Control of transcriptional elongation

Elongation is becoming increasingly recognized as a critical step in eukaryotic transcriptional regulation. Although traditional genetic and biochemical studies have identified major players of transcriptional elongation, our understanding of the importance and roles of these factors is evolving rapidly through the recent advances in genome-wide and single-molecule technologies. Here, we focus on how elongation can modulate the transcriptional outcome through the rate-liming step of RNA polymerase II (Pol II) pausing near promoters and how the participating factors were identified. Among the factors we describe are the pausing factors--NELF (negative elongation factor) and DSIF (DRB sensitivity-inducing factor)--and P-TEFb (positive elongation factor b), which is the key player in pause release. We also describe the high-resolution view of Pol II pausing and propose nonexclusive models for how pausing is achieved. We then discuss Pol II elongation through the bodies of genes and the roles of FACT and SPT6, factors that allow Pol II to move through nucleosomes.

GAGA factor repression of transcription is a rare event but the negative regulation of Trl is conserved in Drosophila species

GAGA is a highly conserved Drosophila transcription factor encoded by the Trithorax-like (Trl) gene. While GAGA usually activates transcription, it represses its own promoter. Here we show that GAGA-mediated repression of Trl is conserved between two distant Drosophila species. A detailed promoter study showed that GAGA repressive activity can't be attributed to any discrete element in the Trl promoter. Genome-wide analysis of the transcriptome in S2 cells indicated that repression of Trl is very likely unique, being GAGA factor a transactivator for all the other promoters. Taken together, our results suggest a new mechanism to explain GAGA-mediated repression that involves a dose-dependent change in the architecture of the Trl promoter.

Sequence-selective DNA recognition with peptide-bisbenzamidine conjugates

Transcription factors (TFs) are specialized proteins that play a key role in the regulation of genetic expression. Their mechanism of action involves the interaction with specific DNA sequences, which usually takes place through specialized domains of the protein. However, achieving an efficient binding usually requires the presence of the full protein. This is the case for bZIP and zinc finger TF families, which cannot interact with their target sites when the DNA binding fragments are presented as isolated monomers. Herein it is demonstrated that the DNA binding of these monomeric peptides can be restored when conjugated to aza-bisbenzamidines, which are readily accessible molecules that interact with A/T-rich sites by insertion into their minor groove. Importantly, the fluorogenic properties of the aza-benzamidine unit provide details of the DNA interaction that are eluded in electrophoresis mobility shift assays (EMSA). The hybrids based on the GCN4 bZIP protein preferentially bind to composite sequences containing tandem bisbenzamidine-GCN4 binding sites (TCAT⋅AAATT). Fluorescence reverse titrations show an interesting multiphasic profile consistent with the formation of competitive nonspecific complexes at low DNA/peptide ratios. On the other hand, the conjugate with the DNA binding domain of the zinc finger protein GAGA binds with high affinity (KD≈12 nM) and specificity to a composite AATTT⋅GAGA sequence containing both the bisbenzamidine and the TF consensus binding sites.

New insights into DNA recognition by zinc fingers revealed by structural analysis of the oncoprotein ZNF217

Classical zinc fingers (ZFs) are one of the most abundant and best characterized DNA-binding domains. Typically, tandem arrays of three or more ZFs bind DNA target sequences with high affinity and specificity, and the mode of DNA recognition is sufficiently well understood that tailor-made ZF-based DNA-binding proteins can be engineered. We have shown previously that a two-zinc finger unit found in the transcriptional coregulator ZNF217 recognizes DNA but with an affinity and specificity that is lower than other ZF arrays. To investigate the basis for these differences, we determined the structure of a ZNF217-DNA complex. We show that although the overall position of the ZFs on the DNA closely resembles that observed for other ZFs, the side-chain interaction pattern differs substantially from the canonical model. The structure also reveals the presence of two methyl-π interactions, each featuring a tyrosine contacting a thymine methyl group. To our knowledge, interactions of this type have not previously been described in classical ZF-DNA complexes. Finally, we investigated the sequence specificity of this two-ZF unit and discuss how ZNF217 might discriminate its target DNA sites in the cell.

Solution structure of an archaeal DNA binding protein with an eukaryotic zinc finger fold

While the basal transcription machinery in archaea is eukaryal-like, transcription factors in archaea and their viruses are usually related to bacterial transcription factors. Nevertheless, some of these organisms show predicted classical zinc fingers motifs of the C2H2 type, which are almost exclusively found in proteins of eukaryotes and most often associated with transcription regulators. In this work, we focused on the protein AFV1p06 from the hyperthermophilic archaeal virus AFV1. The sequence of the protein consists of the classical eukaryotic C2H2 motif with the fourth histidine coordinating zinc missing, as well as of N- and C-terminal extensions. We showed that the protein AFV1p06 binds zinc and solved its solution structure by NMR. AFV1p06 displays a zinc finger fold with a novel structure extension and disordered N- and C-termini. Structure calculations show that a glutamic acid residue that coordinates zinc replaces the fourth histidine of the C2H2 motif. Electromobility gel shift assays indicate that the protein binds to DNA with different affinities depending on the DNA sequence. AFV1p06 is the first experimentally characterised archaeal zinc finger protein with a DNA binding activity. The AFV1p06 protein family has homologues in diverse viruses of hyperthermophilic archaea. A phylogenetic analysis points out a common origin of archaeal and eukaryotic C2H2 zinc fingers.

Microcephaly gene links trithorax and REST/NRSF to control neural stem cell proliferation and differentiation

Microcephaly is a neurodevelopmental disorder causing significantly reduced cerebral cortex size. Many known microcephaly gene products localize to centrosomes, regulating cell fate and proliferation. Here, we identify and characterize a nuclear zinc finger protein, ZNF335/NIF-1, as a causative gene for severe microcephaly, small somatic size, and neonatal death. Znf335 null mice are embryonically lethal, and conditional knockout leads to severely reduced cortical size. RNA-interference and postmortem human studies show that ZNF335 is essential for neural progenitor self-renewal, neurogenesis, and neuronal differentiation. ZNF335 is a component of a vertebrate-specific, trithorax H3K4-methylation complex, directly regulating REST/NRSF, a master regulator of neural gene expression and cell fate, as well as other essential neural-specific genes. Our results reveal ZNF335 as an essential link between H3K4 complexes and REST/NRSF and provide the first direct genetic evidence that this pathway regulates human neurogenesis and neuronal differentiation.

Conservation of the three-dimensional structure in non-homologous or unrelated proteins

In this review, we examine examples of conservation of protein structural motifs in unrelated or non-homologous proteins. For this, we have selected three DNA-binding motifs: the histone fold, the helix-turn-helix motif, and the zinc finger, as well as the globin-like fold. We show that indeed similar structures exist in unrelated proteins, strengthening the concept that three-dimensional conservation might be more important than the primary amino acid sequence.

Exploring translocation of proteins on DNA by NMR

While an extensive body of knowledge has accumulated on the structures of transcription factors, DNA and their complexes from both NMR and crystallography, much less is known at a molecular level regarding the mechanisms whereby transcription factors locate their specific DNA target site within an overwhelming sea of non-specific DNA sites. Indirect kinetic data suggested that three processes are involved in the search procedure: jumping by dissociation of the protein from the DNA followed by re-association at another site, direct transfer from one DNA molecule or segment to another, and one-dimensional sliding. In this brief perspective I summarize recent NMR developments from our laboratory that have permitted direct characterization of the species and molecular mechanisms involved in the target search process, including the detection of highly transient sparsely-populated states. The main tool in these studies involves the application of paramagnetic relaxation enhancement, supplemented by z-exchange spectroscopy, lineshape analysis and residual dipolar couplings. These studies led to the first direct demonstration of rotation-coupled sliding of a protein along the DNA and the direct transfer of a protein from one DNA molecule to another without dissociating into free solution.

Analysis and recognition of the GAGA transcription factor binding sites in Drosophila genes

The transcription factor GAGA, encoded by the gene Trl, controls expression of many Drosophila melanogaster genes. We have compiled the presently largest sample (120 sites) of published nucleotide sequences with experimentally confirmed binding to GAGA protein. Analysis of the sample has demonstrated that despite an apparent structural diversity of the GAGA sites, they fall into four distinct groups, namely, (1) the sites containing two GAG trinucleotides with no more than one nucleotide substitution in each and separated by spacers with a length of 1 or 3 nucleotides (GAGnGAG and GAGnnnGAG); (2) the sites containing a single GAGAG motif; (3) (GA)(3-9) microsatellite repeats; and (4) the sites corresponding to three and more direct repeats of GAG trinucleotide homolog and its inverted repeats separated by spacers of various lengths. Using the software package SITECON, the methods were elaborated for recognizing the sites of GAGnGAG (method 1) and GAGnnnGAG (method 2) types in DNA sequences. Experimental verification confirmed the ability to interact with the GAGA factor for 72% of the sites predicted using method 1 and 94.5% of the sites predicted by method 2. Application of the experimentally verified methods to analyzing the localization of potential GAGA binding sites in the target genes of this transcription factor has demonstrated that the 5'-untranslated regions (5'UTRs) and first introns are enriched for these sites (two-threefold relative to the average occurrence frequency in the D. melanogaster genome) as compared with a moderate enrichment (not exceeding 1.5-fold) of promoter regions (-4000/+200 bp or -1000/+100 bp).