The zinc finger-associated domain of the Drosophila transcription factor grauzone is a novel zinc-coordinating protein-protein interaction module

Escalation of genome defense capacity enables control of an expanding meiotic driver

From RNA interference to chromatin silencing, diverse genome defense pathways silence selfish genetic elements to safeguard genome integrity. Despite their diversity, different defense pathways share a modular organization, where numerous specificity factors identify diverse targets and common effectors silence them. In the PIWI-interacting RNA (piRNA) pathway, target RNAs are first identified by complementary base pairing with piRNAs and then silenced by PIWI-clade nucleases. Such a binary architecture allows the defense systems to be readily adaptable, where new targets can be captured via innovation of specificity factors. Thus, our current understanding of genome defense against lineage-specific selfish genes has been largely limited to specificity factor innovations, while it remains poorly understood whether other types of innovations are required. Here, we describe a new type of innovation, which escalates the genome defense capacity to control a recently expanded selfish gene in Drosophila melanogaster. Through a targeted RNAi screen for repressors of Stellate-a recently evolved meiotic driver-we identified a defense factor, Trailblazer. Trailblazer is a transcription factor that promotes the expression of two PIWI-clade nucleases, Aub and AGO3, to match Stellate in abundance. Recent innovation in the DNA-binding domain of Trailblazer enabled it to elevate Aub and AGO3 expression, thereby escalating the silencing capacity of piRNA pathway to tame expanded Stellate and safeguard fertility. As copy-number expansion is a recurrent feature of diverse selfish genes across the tree of life, we envision that augmenting the defense capacity to quantitatively match selfish genes is a repeatedly employed defense strategy in evolution.

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.

A C2H2 zinc finger transcription factor of the whitefly Bemisia tabaci interacts with the capsid proteins of begomoviruses and inhibits virus retention

Begomoviruses are a group of ssDNA viruses exclusively transmitted by the whitefly Bemisia tabaci and constrain vegetable production in the old and new worlds. Although multiple molecular determinants governing the transmission of begomoviruses by whiteflies have been unravelled, factors critical for transmission majorly remain unknown. In this study, a whitefly C2H2 zinc finger (ZF) protein, 100% identical to the vascular endothelial ZF-like gene (vezf) protein was confirmed to interact with the CP of both old- and new-world begomoviruses. This was achieved by a yeast two-hybrid (Y2H) system screening of a whitefly cDNA library using capsid protein (CP) of TYLCV as a bait. In silico annotation of vezf protein revealed that it contains a N-terminal ZF-associated domain (ZAD) alongside multiple C2H2 ZF domains on the C-terminal end. ZAD-ZF proteins form the most abundant class of transcription factors within insects. Herein, we validated the interaction of vezf with four diverse begomoviruses and its functional role in begomovirus transmission. Silencing of the vezf gene of B. tabaci led to increased retention of three diverse begomoviruses tested. Vezf is the first insect transcription factor identified to interact with plant viruses and can be crucial to understand the possible mechanisms by which plant viruses modulate transcription of their insect vectors during transmission.

AoSAP8-P encoding A20 and/or AN1 type zinc finger protein in asparagus officinalis L. Improving stress tolerance in transgenic Nicotiana sylvestris

The full length CDS of an A20 and AN1 type zinc finger gene (named AoSAP8-P), located nearby the male specific Y chromosome (MSY) region of Asparagus officinalis (garden asparagus) was amplified by RT-PCR from purple passion. This gene, predicted as the stress associated protein (SAPs) gene families, encodes 172 amino acids with multiple cis elements including light, stress response box, MYB and ERF binding sites on its promoter. To analyze its function, the gene expression of different organs in different asparagus gender were analyzed and the overexpressed transgenic Nicotiana sylvestris lines were generated. The results showed the gene was highly expressed in both flower and root of male garden asparagus; the germination rate of seeds of the T2 transgenic lines (T2-5-4 and T2-7-1) under the stress conditions of 125 mM NaCl and 150 mM mannitol were significantly higher than the wild type (WT) respectively. When the potted T2-5-4, T2-7-1 lines and WT were subjected to drought stress for 24 days and the leaf discs immerged into 20 % PEG6000 and 300 mM NaCl solution for 48 h respectively, the T2-5-4 and T2-7-1 with AoSAP8-P expression showed stronger drought, salt and osmotic stress tolerance. When compared, the effects of AoSAP8-P overexpression on productive development showed that the flowering time of transgenic lines, were ∼ 9 day earlier with larger but fewer pollens than its WT counterparts. However, there were no significant differences in anthers, stigmas and pollen viability between the transgenic lines and WT. Our results suggested that, the AoSAP8-P gene plays a role in improving the stress resistance and shortening seeds generation time for perianal survival during the growth and development of garden asparagus.

Protein-protein docking and molecular dynamics studies of sericin and cocoonase of silkworm: an insight for cocoon softening

Cocoonase is known to digest the sericin protein that encapsulates the silkworm cocoon's fibroin protein. Silk fibroin and sericin are two types of proteins that make up silk, and accounts for around 20-30% of the overall cocoon weight. The aim of the study was to see the protein-protein interaction (PPI) and molecular dynamic study of sericin, cocoonase and protein-protein docked complex of silkworm by computational approaches. Here motif analysis, phylogenetic analysis, principal component analysis, root-mean-square deviation (RMSD), root mean square fluctuation, radius of gyration, structural and functional study of cocoonase and sericin as well as molecular docking study were carried out. The 33 amino acid residues of cocoonase shows interaction with 38 aa residues of sericin involving 4 disulphide bonds, 22 hydrogen bonds and 319 non-bonded contacts. The confirmational stability and flexibility of both the proteins as well as protein-protein complex were achieved at 70 ns of MD simulation study. RMSD-based data indicated that cocoonase is more stable than sericin and complex, and complex has a greater fluctuation with more compact (higher Rg) value than cocoonase and sericin, inferring higher conformational stability and flexibility of protein-protein complex than cocoonase and sericin. This study provides a new dimension for PPI study by computational approaches.Communicated by Ramaswamy H. Sarma.

Twenty-seven ZAD-ZNF genes of Drosophila melanogaster are orthologous to the embryo polarity determining mosquito gene cucoid

The C2H2 zinc finger gene cucoid establishes anterior-posterior (AP) polarity in the early embryo of culicine mosquitoes. This gene is unrelated to genes that establish embryo polarity in other fly species (Diptera), such as the homeobox gene bicoid, which serves this function in the traditional model organism Drosophila melanogaster. The cucoid gene is a conserved single copy gene across lower dipterans but nothing is known about its function in other species, and its evolution in higher dipterans, including Drosophila, is unresolved. We found that cucoid is a member of the ZAD-containing C2H2 zinc finger (ZAD-ZNF) gene family and is orthologous to 27 of the 91 members of this family in D. melanogaster, including M1BP, ranshi, ouib, nom, zaf1, odj, Nnk, trem, Zif, and eighteen uncharacterized genes. Available knowledge of the functions of cucoid orthologs in Drosophila melanogaster suggest that the progenitor of this lineage specific expansion may have played a role in regulating chromatin. We also describe many aspects of the gene duplication history of cucoid in the brachyceran lineage of D. melanogaster, thereby providing a framework for predicting potential redundancies among these genes in D. melanogaster.

Resolving the zinc binding capacity of honey bee vitellogenin and locating its putative binding sites

The protein vitellogenin (Vg) plays a central role in lipid transportation in most egg-laying animals. High Vg levels correlate with stress resistance and lifespan potential in honey bees (Apis mellifera). Vg is the primary circulating zinc-carrying protein in honey bees. Zinc is an essential metal ion in numerous biological processes, including the function and structure of many proteins. Measurements of Zn²⁺ suggest a variable number of ions per Vg molecule in different animal species, but the molecular implications of zinc-binding by this protein are not well-understood. We used inductively coupled plasma mass spectrometry to determine that, on average, each honey bee Vg molecule binds 3 Zn²⁺ -ions. Our full-length protein structure and sequence analysis revealed seven potential zinc-binding sites. These are located in the β-barrel and α-helical subdomains of the N-terminal domain, the lipid binding site, and the cysteine-rich C-terminal region of unknown function. Interestingly, two potential zinc-binding sites in the β-barrel can support a proposed role for this structure in DNA-binding. Overall, our findings suggest that honey bee Vg bind zinc at several functional regions, indicating that Zn²⁺ -ions are important for many of the activities of this protein. In addition to being potentially relevant for other egg-laying species, these insights provide a platform for studies of metal ions in bee health, which is of global interest due to recent declines in pollinator numbers.

Safeguarding Drosophila female germ cell identity depends on an H3K9me3 mini domain guided by a ZAD zinc finger protein

H3K9me3-based gene silencing is a conserved strategy for securing cell fate, but the mechanisms controlling lineage-specific installation of this epigenetic mark remain unclear. In Drosophila, H3K9 methylation plays an essential role in securing female germ cell fate by silencing lineage inappropriate phf7 transcription. Thus, phf7 regulation in the female germline provides a powerful system to dissect the molecular mechanism underlying H3K9me3 deposition onto protein coding genes. Here we used genetic studies to identify the essential cis-regulatory elements, finding that the sequences required for H3K9me3 deposition are conserved across Drosophila species. Transposable elements are also silenced by an H3K9me3-mediated mechanism. But our finding that phf7 regulation does not require the dedicated piRNA pathway components, piwi, aub, rhino, panx, and nxf2, indicates that the mechanisms of H3K9me3 recruitment are distinct. Lastly, we discovered that an uncharacterized member of the zinc finger associated domain (ZAD) containing C2H2 zinc finger protein family, IDENTITY CRISIS (IDC; CG4936), is necessary for H3K9me3 deposition onto phf7. Loss of idc in germ cells interferes with phf7 transcriptional regulation and H3K9me3 deposition, resulting in ectopic PHF7 protein expression. IDC's role is likely to be direct, as it localizes to a conserved domain within the phf7 gene. Collectively, our findings support a model in which IDC guides sequence-specific establishment of an H3K9me3 mini domain, thereby preventing accidental female-to-male programming.

The Drosophila ZAD zinc finger protein Kipferl guides Rhino to piRNA clusters

RNA interference systems depend on the synthesis of small RNA precursors whose sequences define the target spectrum of these silencing pathways. The Drosophila Heterochromatin Protein 1 (HP1) variant Rhino permits transcription of PIWI-interacting RNA (piRNA) precursors within transposon-rich heterochromatic loci in germline cells. Current models propose that Rhino's specific chromatin occupancy at piRNA source loci is determined by histone marks and maternally inherited piRNAs, but also imply the existence of other, undiscovered specificity cues. Here, we identify a member of the diverse family of zinc finger associated domain (ZAD)-C2H2 zinc finger proteins, Kipferl, as critical Rhino cofactor in ovaries. By binding to guanosine-rich DNA motifs and interacting with the Rhino chromodomain, Kipferl recruits Rhino to specific loci and stabilizes it on chromatin. In kipferl mutant flies, Rhino is lost from most of its target chromatin loci and instead accumulates on pericentromeric Satellite arrays, resulting in decreased levels of transposon targeting piRNAs and impaired fertility. Our findings reveal that DNA sequence, in addition to the H3K9me3 mark, determines the identity of piRNA source loci and provide insight into how Rhino might be caught in the crossfire of genetic conflicts.

Transcriptional Regulators of Ecdysteroid Biosynthetic Enzymes and Their Roles in Insect Development

Steroid hormones are responsible for coordinating many aspects of biological processes in most multicellular organisms, including insects. Ecdysteroid, the principal insect steroid hormone, is biosynthesized from dietary cholesterol or plant sterols. In the last 20 years, a number of ecdysteroidogenic enzymes, including Noppera-bo, Neverland, Shroud, Spook/Spookier, Cyp6t3, Phantom, Disembodied, Shadow, and Shade, have been identified and characterized in molecular genetic studies using the fruit fly Drosophila melanogaster. These enzymes are encoded by genes collectively called the Halloween genes. The transcriptional regulatory network, governed by multiple regulators of transcription, chromatin remodeling, and endoreplication, has been shown to be essential for the spatiotemporal expression control of Halloween genes in D. melanogaster. In this review, we summarize the latest information on transcriptional regulators that are crucial for controlling the expression of ecdysteroid biosynthetic enzymes and their roles in insect development.

An RNA-interference screen in Drosophila to identify ZAD-containing C2H2 zinc finger genes that function in female germ cells

The zinc finger-associated domain (ZAD) is present in over 90 C2H2 zinc finger (ZNF) proteins. Despite their abundance, only a few ZAD-ZNF genes have been characterized to date. Here, we systematically analyze the function of 68 ZAD-ZNF genes in Drosophila female germ cells by performing an in vivo RNA-interference screen. We identified eight ZAD-ZNF genes required for oogenesis, and based on further characterization of the knockdown phenotypes, we uncovered defects broadly consistent with functions in germ cell specification and/or survival, early differentiation, and egg chamber maturation. These results provide a candidate pool for future studies aimed at functionalization of this large but poorly characterized gene family.

Structural basis of diversity and homodimerization specificity of zinc-finger-associated domains in Drosophila

In arthropods, zinc finger-associated domains (ZADs) are found at the N-termini of many DNA-binding proteins with tandem arrays of Cys2-His2 zinc fingers (ZAD-C2H2 proteins). ZAD-C2H2 proteins undergo fast evolutionary lineage-specific expansion and functional diversification. Here, we show that all ZADs from Drosophila melanogaster form homodimers, but only certain ZADs with high homology can also heterodimerize. CG2712, for example, is unable to heterodimerize with its paralog, the previously characterized insulator protein Zw5, with which it shares 46% homology. We obtained a crystal structure of CG2712 protein's ZAD domain that, in spite of a low sequence homology, has similar spatial organization with the only known ZAD structure (from Grauzone protein). Steric clashes prevented the formation of heterodimers between Grauzone and CG2712 ZADs. Using detailed structural analysis, site-directed mutagenesis, and molecular dynamics simulations, we demonstrated that rapid evolutionary acquisition of interaction specificity was mediated by the more energy-favorable formation of homodimers in comparison to heterodimers, and that this specificity was achieved by multiple amino acid substitutions resulting in the formation or breaking of stabilizing interactions. We speculate that specific homodimerization of ZAD-C2H2 proteins is important for their architectural role in genome organization.

CTCF As an Example of DNA-Binding Transcription Factors Containing Clusters of C2H2-Type Zinc Fingers

In mammals, most of the boundaries of topologically associating domains and all well-studied insulators are rich in binding sites for the CTCF protein. According to existing experimental data, CTCF is a key factor in the organization of the architecture of mammalian chromosomes. A characteristic feature of the CTCF is that the central part of the protein contains a cluster consisting of eleven domains of C2H2-type zinc fingers, five of which specifically bind to a long DNA sequence conserved in most animals. The class of transcription factors that carry a cluster of C2H2-type zinc fingers consisting of five or more domains (C2H2 proteins) is widely represented in all groups of animals. The functions of most C2H2 proteins still remain unknown. This review presents data on the structure and possible functions of these proteins, using the example of the vertebrate CTCF protein and several well- characterized C2H2 proteins in Drosophila and mammals.

Innovation of heterochromatin functions drives rapid evolution of essential ZAD-ZNF genes in Drosophila

Contrary to dogma, evolutionarily young and dynamic genes can encode essential functions. We find that evolutionarily dynamic ZAD-ZNF genes, which encode the most abundant class of insect transcription factors, are more likely to encode essential functions in Drosophila melanogaster than ancient, conserved ZAD-ZNF genes. We focus on the Nicknack ZAD-ZNF gene, which is evolutionarily young, poorly retained in Drosophila species, and evolves under strong positive selection. Yet we find that it is necessary for larval development in D. melanogaster. We show that Nicknack encodes a heterochromatin-localizing protein like its paralog Oddjob, also an evolutionarily dynamic yet essential ZAD-ZNF gene. We find that the divergent D. simulans Nicknack protein can still localize to D. melanogaster heterochromatin and rescue viability of female but not male Nicknack-null D. melanogaster. Our findings suggest that innovation for rapidly changing heterochromatin functions might generally explain the essentiality of many evolutionarily dynamic ZAD-ZNF genes in insects.

Identification and effect of Zf-AD-containing C2H2 zinc finger genes on BmNPV replication in the silkworm (Bombyx mori)

Zf-AD-containing C2H2 zinc -finger genes (ZAD) are uniquely present and have lineage-specific expansion in arthropods. Arthropods are also the hosts of Baculoviruses. We studied the possible relationship between the lineage-specific expansion of ZAD genes and arthropod-Baculovirus co-evolution. We used the silkworm (Bombyx mori) as a model. We identified 73 ZAD genes (BmZAD) in the silkworm. Sequence-based similarity analysis showed that nine clusters involving 28 BmZADs may have undergone species-specific expansion in the silkworm. Expression pattern analysis showed that the BmZADs were divided into five groups. Group I comprised 10 genes with high expression in multiple tissues, suggesting that BmZADs may play roles in the development of various tissues. We identified six BmZADs that could be induced by the Nucleopolyhedrovirus (BmNPV). Among them, BmZAD69 expression is capable of responding to BmNPV infection, and the ZAD domain is indispensable for the function of BmZAD69 in BmNPV replication. We also detected a 3 bp deletion at 1.7 kb upstream of BmZAD69, which may make it more sensitive to BmNPV infection, and thus elevate the BmNPV resistance in Qiufeng_N, a strain with strong virus resistance. These data suggest that BmZADs may be involved in BmNPV infection and that ZAD genes may play a role in arthropod-Baculovirus co-evolution.

Small Drosophila zinc finger C2H2 protein with an N-terminal zinc finger-associated domain demonstrates the architecture functions

Recently, the concept has arisen that a special class of architectural proteins exists, which are responsible not only for global chromosome architecture but also for the local regulation of enhancer-promoter interactions. Here, we describe a new architectural protein, with a total size of only 375 aa, which contains an N-terminal zinc finger-associated domain (ZAD) and a cluster of five zinc finger C2H2 domains at the C-terminus. This new protein, named ZAD and Architectural Function 1 protein (ZAF1 protein), is weakly and ubiquitously expressed, with the highest expression levels observed in oocytes and embryos. The cluster of C2H2 domains recognizes a specific 15-bp consensus site, located predominantly in promoters, near transcription start sites. The expression of ZAF1 by a tissue-specific promoter led to the complete blocking of the eye enhancer when clusters of ZAF1 binding sites flanked the eye enhancer in transgenic lines, suggesting that the loop formed by the ZAF1 protein leads to insulation. The ZAF1 protein also supported long-range interactions between the yeast GAL4 activator and the white promoter in transgenic Drosophila lines. A mutant protein lacking the ZAD failed to block the eye enhancer or to support distance interactions in transgenic lines. Taken together, these results suggest that ZAF1 is a minimal architectural protein that can be used to create a convenient model for studying the mechanisms of distance interactions.

Opbp is a new architectural/insulator protein required for ribosomal gene expression

A special class of poorly characterized architectural proteins is required for chromatin topology and enhancer–promoter interactions. Here, we identify Opbp as a new Drosophila architectural protein, interacting with CP190 both in vivo and in vitro. Opbp binds to a very restrictive set of genomic regions, through a rare sequence specific motif. These sites are co-bound by CP190 in vivo, and generally located at bidirectional promoters of ribosomal protein genes. We show that Opbp is essential for viability, and loss of opbp function, or destruction of its motif, leads to reduced ribosomal protein gene expression, indicating a functional role in promoter activation. As characteristic of architectural/insulator proteins, the Opbp motif is sufficient for distance-dependent reporter gene activation and enhancer-blocking activity, suggesting an Opbp-mediated enhancer–promoter interaction. Rather than having a constitutive role, Opbp represents a new type of architectural protein with a very restricted, yet essential, function in regulation of housekeeping gene expression.

Cooperative Control of Ecdysone Biosynthesis in Drosophila by Transcription Factors Séance, Ouija Board, and Molting Defective

Ecdysteroids are steroid hormones that control many aspects of development and physiology. During larval development, ecdysone is synthesized in an endocrine organ called the prothoracic gland through a series of ecdysteroidogenic enzymes encoded by the Halloween genes. The expression of the Halloween genes is highly restricted and dynamic, indicating that their spatiotemporal regulation is mediated by their tight transcriptional control. In this study, we report that three zinc finger-associated domain (ZAD)-C₂H₂ zinc finger transcription factors—Séance (Séan), Ouija board (Ouib), and Molting defective (Mld)—cooperatively control ecdysone biosynthesis in the fruit fly Drosophila melanogaster. Séan and Ouib act in cooperation with Mld to positively regulate the transcription of neverland and spookier, respectively, two Halloween genes. Remarkably, loss-of-function mutations in séan, ouib, or mld can be rescued by the expression of neverland, spookier, or both, respectively. These results suggest that the three transcription factors have distinct roles in coordinating the expression of just two genes in Drosophila. Given that neverland and spookier are located in constitutive heterochromatin, Séan, Ouib, and Mld represent the first example of a transcription factor subset that regulates genes located in constitutive heterochromatin.

piragua encodes a zinc finger protein required for development in Drosophila

We isolated and characterized embryonic lethal mutations in piragua (prg). The prg locus encodes a protein with an amino terminus Zinc Finger-Associated-Domain (ZAD) and nine C₂H₂ zinc fingers (ZF). prg mRNA and protein expression during embryogenesis is dynamic with widespread maternal contribution, and subsequent expression in epithelial precursors. About a quarter of prg mutant embryos do not develop cuticle, and from those that do a small fraction have cuticular defects. Roughly half of prg mutants die during embryogenesis. prg mutants have an extended phenocritical period encompassing embryogenesis and first instar larval stage, since the other half of prg mutants die as first or second instar larvae. During dorsal closure, time-lapse high-resolution imaging shows defects arising out of sluggishness in closure, resolving at times in failures of closure. prg is expressed in imaginal discs, and is required for imaginal development. prg was identified in imaginal tissue in a cell super competition screen, together with other genes, like flower. We find that flower mutations are also embryonic lethal with a similar phenocritical period and strong embryonic mutant phenotypes (head involution defects, primarily). The two loci interact genetically in the embryo, as they increase embryonic mortality to close to 90% with the same embryonic phenotypes (dorsal closure and head involution defects, plus lack of cuticle). Mutant prg clones generated in developing dorsal thorax and eye imaginal tissue have strong developmental defects (lack of bristles and ommatidial malformations). prg is required in several developmental morphogenetic processes.

Architectural proteins Pita, Zw5,and ZIPIC contain homodimerization domain and support specific long-range interactions in Drosophila

According to recent models, as yet poorly studied architectural proteins appear to be required for local regulation of enhancer-promoter interactions, as well as for global chromosome organization. Transcription factors ZIPIC, Pita and Zw5 belong to the class of chromatin insulator proteins and preferentially bind to promoters near the TSS and extensively colocalize with cohesin and condensin complexes. ZIPIC, Pita and Zw5 are structurally similar in containing the N-terminal zinc finger-associated domain (ZAD) and different numbers of C2H2-type zinc fingers at the C-terminus. Here we have shown that the ZAD domains of ZIPIC, Pita and Zw5 form homodimers. In Drosophila transgenic lines, these proteins are able to support long-distance interaction between GAL4 activator and the reporter gene promoter. However, no functional interaction between binding sites for different proteins has been revealed, suggesting that such interactions are highly specific. ZIPIC facilitates long-distance stimulation of the reporter gene by GAL4 activator in yeast model system. Many of the genomic binding sites of ZIPIC, Pita and Zw5 are located at the boundaries of topologically associated domains (TADs). Thus, ZAD-containing zinc-finger proteins can be attributed to the class of architectural proteins.

Transcription factors, chromatin proteins and the diversification of Hemiptera

Availability of complete genomes provides a means to explore the evolution of enormous developmental, morphological, and behavioral diversity among insects. Hemipterans in particular show great diversity of both morphology and life history within a single order. To better understand the role of transcription regulators in the diversification of hemipterans, using sequence profile searches and hidden Markov models we computationally analyzed transcription factors (TFs) and chromatin proteins (CPs) in the recently available Rhodnius prolixus genome along with 13 other insect and 4 non-insect arthropod genomes. We generated a comprehensive collection of TFs and CPs across arthropods including 303 distinct types of domains in TFs and 139 in CPs. This, along with the availability of two hemipteran genomes, R. prolixus and Acyrthosiphon pisum, helped us identify possible determinants for their dramatic morphological and behavioral divergence. We identified five domain families (i.e. Pipsqueak, SAZ/MADF, THAP, FLYWCH and BED finger) as having undergone differential patterns of lineage-specific expansion in hemipterans or within hemipterans relative to other insects. These expansions appear to be at least in part driven by transposons, with the DNA-binding domains of transposases having provided the raw material for emergence of new TFs. Our analysis suggests that while R. prolixus probably retains a state closer to the ancestral hemipteran, A. pisum represents a highly derived state, with the emergence of asexual reproduction potentially favoring genome duplication and transposon expansion. Both hemipterans are predicted to possess active DNA methylation systems. However, in the course of their divergence, aphids seem to have expanded the ancestral hemipteran DNA methylation along with a distinctive linkage to the histone methylation system, as suggested by expansion of SET domain methylases, including those fused to methylated CpG recognition domains. Thus, differential use of DNA methylation and histone methylation might have played a role in emergence of polyphenism and cyclic parthenogenesis from the ancestral hemipteran.

The Drosophila Zinc Finger Transcription Factor Ouija Board Controls Ecdysteroid Biosynthesis through Specific Regulation of spookier

Steroid hormones are crucial for many biological events in multicellular organisms. In insects, the principal steroid hormones are ecdysteroids, which play essential roles in regulating molting and metamorphosis. During larval and pupal development, ecdysteroids are synthesized in the prothoracic gland (PG) from dietary cholesterol via a series of hydroxylation and oxidation steps. The expression of all but one of the known ecdysteroid biosynthetic enzymes is restricted to the PG, but the transcriptional regulatory networks responsible for generating such exquisite tissue-specific regulation is only beginning to be elucidated. Here, we report identification and characterization of the C₂H₂-type zinc finger transcription factor Ouija board (Ouib) necessary for ecdysteroid production in the PG in the fruit fly Drosophila melanogaster. Expression of ouib is predominantly limited to the PG, and genetic null mutants of ouib result in larval developmental arrest that can be rescued by administrating an active ecdysteroid. Interestingly, ouib mutant animals exhibit a strong reduction in the expression of one ecdysteroid biosynthetic enzyme, spookier. Using a cell culture-based luciferase reporter assay, Ouib protein stimulates transcription of spok by binding to a specific ~15 bp response element in the spok PG enhancer element. Most remarkable, the developmental arrest phenotype of ouib mutants is rescued by over-expression of a functionally-equivalent paralog of spookier. These observations imply that the main biological function of Ouib is to specifically regulate spookier transcription during Drosophila development.

Steroid hormones are crucial for development and reproduction in multicellular organisms. The spatially-restricted expression of almost all steroid biosynthesis genes is key to the specialization of steroid producing cells. In the last decade, insects have become the focus for research on the biosynthesis of the principal steroid hormones, ecdysteroids. However, the transcriptional regulatory mechanisms controlling the ecdysteroid biosynthesis genes are largely unknown. Here we show that a novel zinc finger transcription factor Ouija board (Ouib) is essential for activating the expression of one ecdysteroid biosynthesis gene, spookier, in the ecdysteroid producing cells. Ouib is the first invertebrate transcription factor that is predominantly expressed in the steroidogenic organs and essential for development via inducing expression of the steroidogenic gene. In addition, this is the first report showing the catalytic step-specific control of steroid hormone biosynthesis through transcriptional regulation.

Pervasive variation of transcription factor orthologs contributes to regulatory network evolution

Differences in transcriptional regulatory networks underlie much of the phenotypic variation observed across organisms. Changes to cis-regulatory elements are widely believed to be the predominant means by which regulatory networks evolve, yet examples of regulatory network divergence due to transcription factor (TF) variation have also been observed. To systematically ascertain the extent to which TFs contribute to regulatory divergence, we analyzed the evolution of the largest class of metazoan TFs, Cys2-His2 zinc finger (C2H2-ZF) TFs, across 12 Drosophila species spanning ~45 million years of evolution. Remarkably, we uncovered that a significant fraction of all C2H2-ZF 1-to-1 orthologs in flies exhibit variations that can affect their DNA-binding specificities. In addition to loss and recruitment of C2H2-ZF domains, we found diverging DNA-contacting residues in ~44% of domains shared between D. melanogaster and the other fly species. These diverging DNA-contacting residues, found in ~70% of the D. melanogaster C2H2-ZF genes in our analysis and corresponding to ~26% of all annotated D. melanogaster TFs, show evidence of functional constraint: they tend to be conserved across phylogenetic clades and evolve slower than other diverging residues. These same variations were rarely found as polymorphisms within a population of D. melanogaster flies, indicating their rapid fixation. The predicted specificities of these dynamic domains gradually change across phylogenetic distances, suggesting stepwise evolutionary trajectories for TF divergence. Further, whereas proteins with conserved C2H2-ZF domains are enriched in developmental functions, those with varying domains exhibit no functional enrichments. Our work suggests that a subset of highly dynamic and largely unstudied TFs are a likely source of regulatory variation in Drosophila and other metazoans.

Mechanisms and proteins involved in long-distance interactions

Due to advances in genome-wide technologies, consistent distant interactions within chromosomes of higher eukaryotes have been revealed. In particular, it has been shown that enhancers can specifically and directly interact with promoters by looping out intervening sequences, which can be up to several hundred kilobases long. This review is focused on transcription factors that are supposed to be involved in long-range interactions. Available data are in agreement with the model that several known transcription factors and insulator proteins belong to an abundant but poorly studied class of proteins that are responsible for chromosomal architecture.

Distinct mechanisms of transcriptional pausing orchestrated by GAGA factor and M1BP, a novel transcription factor

Thousands of genes in Drosophila have Pol II paused in the promoter proximal region. Almost half of these genes are associated with either GAGA factor (GAF) or a newly discovered factor we call M1BP. Although both factors dictate the association of Pol II at their target promoters, they are nearly mutually exclusive on the genome and mediate different mechanisms of regulation. High-resolution mapping of Pol II using permanganate-ChIP-seq indicates that pausing on M1BP genes is transient and could involve the +1 nucleosome. In contrast, pausing on GAF genes is much stronger and largely independent of nucleosomes. Distinct regulatory mechanisms are reflected by transcriptional plasticity: M1BP genes are constitutively expressed throughout development while GAF genes exhibit much greater developmental specificity. M1BP binds a core promoter element called Motif 1. Motif 1 potentially directs a distinct transcriptional mechanism from the canonical TATA box, which does not correlate with paused Pol II on the genomic scale. In contrast to M1BP and GAF genes, a significant portion of TATA box genes appear to be controlled at preinitiation complex formation.

Structural analysis and dimerization profile of the SCAN domain of the pluripotency factor Zfp206

Zfp206 (also named as Zscan10) belongs to the subfamily of C₂H₂ zinc finger transcription factors, which is characterized by the N-terminal SCAN domain. The SCAN domain mediates self-association and association between the members of SCAN family transcription factors, but the structural basis and selectivity determinants for complex formation is unknown. Zfp206 is important for maintaining the pluripotency of embryonic stem cells presumably by combinatorial assembly of itself or other SCAN family members on enhancer regions. To gain insights into the folding topology and selectivity determinants for SCAN dimerization, we solved the 1.85 Å crystal structure of the SCAN domain of Zfp206. In vitro binding studies using a panel of 20 SCAN proteins indicate that the SCAN domain Zfp206 can selectively associate with other members of SCAN family transcription factors. Deletion mutations showed that the N-terminal helix 1 is critical for heterodimerization. Double mutations and multiple mutations based on the Zfp206SCAN–Zfp110SCAN model suggested that domain swapped topology is a possible preference for Zfp206SCAN–Zfp110SCAN heterodimer. Together, we demonstrate that the Zfp206SCAN constitutes a protein module that enables C₂H₂ transcription factor dimerization in a highly selective manner using a domain-swapped interface architecture and identify novel partners for Zfp206 during embryonal development.

The Drosophila zinc finger protein trade embargo is required for double strand break formation in meiosis

Homologous recombination in meiosis is initiated by the programmed induction of double strand breaks (DSBs). Although the Drosophila Spo11 ortholog Mei-W68 is required for the induction of DSBs during meiotic prophase, only one other protein (Mei-P22) has been shown to be required for Mei-W68 to exert this function. We show here that the chromatin-associated protein Trade Embargo (Trem), a C2H2 zinc finger protein, is required to localize Mei-P22 to discrete foci on meiotic chromosomes, and thus to promote the formation of DSBs, making Trem the earliest known function in the process of DSB formation in Drosophila oocytes. We speculate that Trem may act by either directing the binding of Mei-P22 to preferred sites of DSB formation or by altering chromatin structure in a manner that allows Mei-P22 to form foci.

Function and Evolution of C2H2 Zinc Finger Arrays

Krüppel-type or C2H2 zinc fingers represent a dominant DNA-binding motif in eukaryotic transcription factor (TF) proteins. In Krüppel-type (KZNF) TFs, KZNF motifs are arranged in arrays of three to as many as 40 tandem units, which cooperate to define the unique DNA recognition properties of the protein. Each finger contains four amino acids located at specific positions, which are brought into direct contact with adjacent nucleotides in the DNA sequence as the KZNF array winds around the major groove of the alpha helix. This arrangement creates an intimate and potentially predictable relationship between the amino acid sequence of KZNF arrays and the nucleotide sequence of target binding sites. The large number of possible combinations and arrangements of modular KZNF motifs, and the increasing lengths of KZNF arrays in vertebrate species, has created huge repertoires of functionally unique TF proteins. The properties of this versatile DNA-binding motif have been exploited independently many times over the course of evolution, through attachment to effector motifs that confer activating, repressing or other activities to the proteins. Once created, some of these novel inventions have expanded in specific evolutionary clades, creating large families of TFs that are lineage- or species-unique. This chapter reviews the properties and their remarkable evolutionary history of eukaryotic KZNF TF proteins, with special focus on large families that dominate the TF landscapes in different metazoan species.

Mago Nashi, Tsunagi/Y14, and Ranshi form a complex that influences oocyte differentiation in Drosophila melanogaster

During Drosophila melanogaster oogenesis, a germline stem cell divides forming a cyst of 16 interconnected cells. One cell enters the oogenic pathway, and the remaining 15 differentiate as nurse cells. Although directed transport and localization of oocyte differentiation factors within the single cell are indispensible for selection, maintenance, and differentiation of the oocyte, the mechanisms regulating these events are poorly understood. Mago Nashi and Tsunagi/Y14, core components of the exon junction complex (a multiprotein complex assembled on spliced RNAs), are essential for restricting oocyte fate to a single cell and for localization of oskar mRNA. Here we provide evidence that Mago Nashi and Tsunagi/Y14 form an oogenic complex with Ranshi, a protein with a zinc finger-associated domain and zinc finger domains. Genetic analyses of ranshi reveal that (1) 16-cell cysts are formed, (2) two cells retain synaptonemal complexes, (3) all cells have endoreplicated DNA (as observed in nurse cells), and (4) oocyte-specific cytoplasmic markers accumulate and persist within a single cell but are not localized within the posterior pole of the presumptive oocyte. Our results indicate that Ranshi interacts with the exon junction complex to localize components essential for oocyte differentiation within the posterior pole of the presumptive oocyte.

Apprehending multicellularity: regulatory networks, genomics, and evolution

The genomic revolution has provided the first glimpses of the architecture of regulatory networks. Combined with evolutionary information, the "network view" of life processes leads to remarkable insights into how biological systems have been shaped by various forces. This understanding is critical because biological systems, including regulatory networks, are not products of engineering but of historical contingencies. In this light, we attempt a synthetic overview of the natural history of regulatory networks operating in the development and differentiation of multicellular organisms. We first introduce regulatory networks and their organizational principles as can be deduced using ideas from the graph theory. We then discuss findings from comparative genomics to illustrate the effects of lineage-specific expansions, gene-loss, and nonprotein-coding DNA on the architecture of networks. We consider the interaction between expansions of transcription factors, and cis regulatory and more general chromatin state stabilizing elements in the emergence of morphological complexity. Finally, we consider a case study of the Notch subnetwork, which is present throughout Metazoa, to examine how such a regulatory system has been pieced together in evolution from new innovations and pre-existing components that were originally functionally distinct.

Structure-function analysis of the THAP zinc finger of THAP1, a large C2CH DNA-binding module linked to Rb/E2F pathways

THAP1, the founding member of a previously uncharacterized large family of cellular proteins (THAP proteins), is a sequence-specific DNA-binding factor that has recently been shown to regulate cell proliferation through modulation of pRb/E2F cell cycle target genes. THAP1 shares its DNA-binding THAP zinc finger domain with Drosophila P element transposase, zebrafish E2F6, and several nematode proteins interacting genetically with the retinoblastoma protein pRb. In this study, we report the three-dimensional structure and structure-function relationships of the THAP zinc finger of human THAP1. Deletion mutagenesis and multidimensional NMR spectroscopy revealed that the THAP domain of THAP1 is an atypical zinc finger of approximately 80 residues, distinguished by the presence between the C2CH zinc coordinating residues of a short antiparallel beta-sheet interspersed by a long loop-helix-loop insertion. Alanine scanning mutagenesis of this loop-helix-loop motif resulted in the identification of a number of critical residues for DNA recognition. NMR chemical shift perturbation analysis was used to further characterize the residues involved in DNA binding. The combination of the mutagenesis and NMR data allowed the mapping of the DNA binding interface of the THAP zinc finger to a highly positively charged area harboring multiple lysine and arginine residues. Together, these data represent the first structure-function analysis of a functional THAP domain, with demonstrated sequence-specific DNA binding activity. They also provide a structural framework for understanding DNA recognition by this atypical zinc finger, which defines a novel family of cellular factors linked to cell proliferation and pRb/E2F cell cycle pathways in humans, fish, and nematodes.

Lineage-specific expansion of the zinc finger associated domain ZAD

The zinc finger associated domain (ZAD), present in almost 100 distinct proteins, characterizes the largest subgroup of C2H2 zinc finger proteins in Drosophila melanogaster and was initially found to be encoded by arthropod genomes only. Here, we report that the ZAD was also present in the last common ancestor of arthropods and vertebrates, and that vertebrate genomes contain a single conserved gene that codes for a ZAD-like peptide. Comparison of the ZAD proteomes of several arthropod species revealed an extensive and species-specific expansion of ZAD-coding genes in higher holometabolous insects, and shows that only few ZAD-coding genes with essential functions in Drosophila melanogaster are conserved. Furthermore, at least 50% of the ZAD-coding genes of Drosophila melanogaster are expressed in the female germline, suggesting a function in oocyte development and/or a requirement during early embryogenesis. Since the majority of the essential ZAD coding genes of Drosophila melanogaster were not conserved during arthropod or at least during insect evolution, we propose that the LSE of ZAD-coding genes shown here may provide the raw material for the evolution of new functions that allow organisms to pursue novel evolutionary paths.

Solution Structure of the THAP Domain from Caenorhabditis elegans C-terminal Binding Protein (CtBP)

The THAP (Thanatos-associated protein) domain is a recently discovered zinc-binding domain found in proteins involved in transcriptional regulation, cell-cycle control, apoptosis and chromatin modification. It contains a single zinc atom ligated by cysteine and histidine residues within a Cys-X(2-4)-Cys-X(35-53)-Cys-X(2)-His consensus. We have determined the NMR solution structure of the THAP domain from Caenorhabditis elegans C-terminal binding protein (CtBP) and show that it adopts a fold containing a treble clef motif, bearing similarity to the zinc finger-associated domain (ZAD) from Drosophila Grauzone. The CtBP THAP domain contains a large, positively charged surface patch and we demonstrate that this domain can bind to double-stranded DNA in an electrophoretic mobility-shift assay. These data, together with existing reports, indicate that THAP domains might exhibit a functional diversity similar to that observed for classical and GATA-type zinc fingers.

Crystal structures of the Mnk2 kinase domain reveal an inhibitory conformation and a zinc binding site

Human mitogen-activated protein kinases (MAPK)-interacting kinases 1 and 2 (Mnk1 and Mnk2) target the translational machinery by phosphorylation of the eukaryotic initiation factor 4E (eIF4E). Here, we present the 2.1 A crystal structure of a nonphosphorylated Mnk2 fragment that encompasses the kinase domain. The results show Mnk-specific features such as a zinc binding motif and an atypical open conformation of the activation segment. In addition, the ATP binding pocket contains an Asp-Phe-Asp (DFD) in place of the canonical magnesium binding Asp-Phe-Gly (DFG) motif. The phenylalanine of this motif sticks into the ATP binding pocket and blocks ATP binding as observed with inhibitor bound and, thus, inactive p38 kinase. Replacement of the DFD by the canonical DFG motif affects the conformation of Mnk2, but not ATP binding and kinase activity. The results suggest that the ATP binding pocket and the activation segment of Mnk2 require conformational switches to provide kinase activity.

A major bristle QTL from a selected population of Drosophila uncovers the zinc-finger transcription factor poils-au-dos, a repressor of achaete-scute

Traditional screens aiming at identifying genes regulating development have relied on mutagenesis. Here, we describe a new gene involved in bristle development, identified through the use of natural variation and selection. Drosophila melanogaster bears a pattern of 11 macrochaetes per heminotum. From a population initially sampled in Marrakech, a strain was selected for an increased number of thoracic macrochaetes. Using recombination and single nucleotide polymorphisms, the factor responsible was mapped to a single locus on the third chromosome, poils au dos, that encodes a zinc-finger-ZAD protein. The original, as well as new, presumed null, alleles of poils au dos, is associated with ectopic achaete-scute expression that results in the additional bristles. This suggests a possible role for Poils au dos as a repressor of achaete and scute. Ectopic expression appears to be independent of the activity of known cis-regulatory enhancer sequences at the achaete-scute complex that mediate activation at specific sites on the notum. The target sequences for Poils au dos activity were mapped to a 14 kb region around scute. In addition, we show that pad interacts synergistically with the repressor hairy and with Dpp signaling in posterior and anterior regions of the notum, respectively.