Adaptive evolution in zinc finger transcription factors

TvZNF1 is a C2H2 zinc finger protein of Trichomonas vaginalis

The zinc fingers proteins (ZNF) are the largest family of DNA binding proteins and can act as transcriptional factors in eukaryotes. ZNF are implicated in activation in response to environmental stimulus by biometals such as Zn²⁺. Many of these proteins have the classical C₂H₂ zinc finger motifs (C₂H₂-ZNFm) of approximately 30 amino acids, where a Zn²⁺ ion is coordinated by two cysteine and two histidine residues. Trichomonas vaginalis is a protozoan parasite than responds to environmental changes including Zn²⁺. Until now has not been described any ZNF that could be involved in the regulation of genic expression of T. vaginalis. Here, we characterized in silico and experimentally an annoted ZNF (TvZNF1) from T. vaginalis and isolated the gene, tvznf1 encoding it. TvZNF1 have eight C₂H₂-ZNFm with residues that maybe involved in the structural stability of DNA binding motifs. In this work we confirmed the Zn²⁺ upregulation expression of tvznf1 gene. Recombinant TvZNF1 was able to bind to specific DNA sequences according to EMSA assay. Additionally, we demonstrated that recombinant TvZNF1 bind to MRE signature in vitro, which strongly suggests its role in transcriptional regulation, similar to the one observed for mammalian MTF-1. This result suggested a conserved mechanism of genic regulation mediated by ZNFs in T. vaginalis.

KRAB zinc finger proteins

Krüppel-associated box domain zinc finger proteins (KRAB-ZFPs) are the largest family of transcriptional regulators in higher vertebrates. Characterized by an N-terminal KRAB domain and a C-terminal array of DNA-binding zinc fingers, they participate, together with their co-factor KAP1 (also known as TRIM28), in repression of sequences derived from transposable elements (TEs). Until recently, KRAB-ZFP/KAP1-mediated repression of TEs was thought to lead to irreversible silencing, and the evolutionary selection of KRAB-ZFPs was considered to be just the host component of an arms race against TEs. However, recent advances indicate that KRAB-ZFPs and their TE targets also partner up to establish species-specific regulatory networks. Here, we provide an overview of the KRAB-ZFP gene family, highlighting how its evolutionary history is linked to that of TEs, and how KRAB-ZFPs influence multiple aspects of development and physiology.

The Role of KRAB-ZFPs in Transposable Element Repression and Mammalian Evolution

Kruppel-associated box zinc-finger proteins (KRAB-ZFPs) make up the largest family of transcription factors in humans. These proteins emerged in the last common ancestor of coelacanth and tetrapods, and have expanded and diversified in the mammalian lineage. Although their mechanism of transcriptional repression has been well studied for over a decade, the DNA-binding activities and the biological functions of these proteins have been largely unexplored. Recent large-scale ChIP-seq studies and loss-of-function experiments have revealed that KRAB-ZFPs play a major role in the recognition and transcriptional silencing of transposable elements (TEs), consistent with an 'arms race model' of KRAB-ZFP evolution against invading TEs. However, this model is insufficient to explain the evolution of many KRAB-ZFPs that appear to domesticate TEs for novel host functions. We highlight some of the mammalian regulatory innovations driven by specific KRAB-ZFPs, including genomic imprinting, meiotic recombination hotspot choice, and placental growth.

Deciphering the DNA code for the function of the Drosophila polydactyl zinc finger protein Suppressor of Hairy-wing

Polydactyl zinc finger (ZF) proteins have prominent roles in gene regulation and often execute multiple regulatory functions. To understand how these proteins perform varied regulation, we studiedDrosophila Suppressor of Hairy-wing [Su(Hw)], an exemplar multifunctional polydactyl ZF protein. We identified separation-of-function (SOF) alleles that encode proteins disrupted in a single ZF that retain one of the Su(Hw) regulatory activities. Through extended in vitro analyses of the Su(Hw) ZF domain, we show that clusters of ZFs bind individual modules within a compound DNA consensus sequence. Through in vivo analysis of SOF mutants, we find that Su(Hw) genomic sites separate into sequence subclasses comprised of combinations of modules, with subclasses enriched for different chromatin features. These data suggest a Su(Hw) code, wherein DNA binding dictates its cofactor recruitment and regulatory output. We propose that similar DNA codes might be used to confer multiple regulatory functions of other polydactyl ZF proteins.

Non-base-contacting residues enable kaleidoscopic evolution of metazoan C2H2 zinc finger DNA binding

Background

The C2H2 zinc finger (C2H2-ZF) is the most numerous protein domain in many metazoans, but is not as frequent or diverse in other eukaryotes. The biochemical and evolutionary mechanisms that underlie the diversity of this DNA-binding domain exclusively in metazoans are, however, mostly unknown.

Results

Here, we show that the C2H2-ZF expansion in metazoans is facilitated by contribution of non-base-contacting residues to DNA binding energy, allowing base-contacting specificity residues to mutate without catastrophic loss of DNA binding. In contrast, C2H2-ZF DNA binding in fungi, plants, and other lineages is constrained by reliance on base-contacting residues for DNA-binding functionality. Reconstructions indicate that virtually every DNA triplet was recognized by at least one C2H2-ZF domain in the common progenitor of placental mammals, but that extant C2H2-ZF domains typically bind different sequences from these ancestral domains, with changes facilitated by non-base-contacting residues.

Conclusions

Our results suggest that the evolution of C2H2-ZFs in metazoans was expedited by the interaction of non-base-contacting residues with the DNA backbone. We term this phenomenon “kaleidoscopic evolution,” to reflect the diversity of both binding motifs and binding motif transitions and the facilitation of their diversification.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1287-y) contains supplementary material, which is available to authorized users.

Multilayered Control of Alternative Splicing Regulatory Networks by Transcription Factors

Networks of coordinated alternative splicing (AS) events play critical roles in development and disease. However, a comprehensive knowledge of the factors that regulate these networks is lacking. We describe a high-throughput system for systematically linking trans-acting factors to endogenous RNA regulatory events. Using this system, we identify hundreds of factors associated with diverse regulatory layers that positively or negatively control AS events linked to cell fate. Remarkably, more than one-third of the regulators are transcription factors. Further analyses of the zinc finger protein Zfp871 and BTB/POZ domain transcription factor Nacc1, which regulate neural and stem cell AS programs, respectively, reveal roles in controlling the expression of specific splicing regulators. Surprisingly, these proteins also appear to regulate target AS programs via binding RNA. Our results thus uncover a large "missing cache" of splicing regulators among annotated transcription factors, some of which dually regulate AS through direct and indirect mechanisms.

Modelling the evolution of transcription factor binding preferences in complex eukaryotes

Transcription factors (TFs) exert their regulatory action by binding to DNA with specific sequence preferences. However, different TFs can partially share their binding sequences due to their common evolutionary origin. This "redundancy" of binding defines a way of organizing TFs in "motif families" by grouping TFs with similar binding preferences. Since these ultimately define the TF target genes, the motif family organization entails information about the structure of transcriptional regulation as it has been shaped by evolution. Focusing on the human TF repertoire, we show that a one-parameter evolutionary model of the Birth-Death-Innovation type can explain the TF empirical repartition in motif families, and allows to highlight the relevant evolutionary forces at the origin of this organization. Moreover, the model allows to pinpoint few deviations from the neutral scenario it assumes: three over-expanded families (including HOX and FOX genes), a set of "singleton" TFs for which duplication seems to be selected against, and a higher-than-average rate of diversification of the binding preferences of TFs with a Zinc Finger DNA binding domain. Finally, a comparison of the TF motif family organization in different eukaryotic species suggests an increase of redundancy of binding with organism complexity.

Proteins Recognizing DNA: Structural Uniqueness and Versatility of DNA-Binding Domains in Stem Cell Transcription Factors

Proteins in the form of transcription factors (TFs) bind to specific DNA sites that regulate cell growth, differentiation, and cell development. The interactions between proteins and DNA are important toward maintaining and expressing genetic information. Without knowing TFs structures and DNA-binding properties, it is difficult to completely understand the mechanisms by which genetic information is transferred between DNA and proteins. The increasing availability of structural data on protein-DNA complexes and recognition mechanisms provides deeper insights into the nature of protein-DNA interactions and therefore, allows their manipulation. TFs utilize different mechanisms to recognize their cognate DNA (direct and indirect readouts). In this review, we focus on these recognition mechanisms as well as on the analysis of the DNA-binding domains of stem cell TFs, discussing the relative role of various amino acids toward facilitating such interactions. Unveiling such mechanisms will improve our understanding of the molecular pathways through which TFs are involved in repressing and activating gene expression.

The evolution of duplicate gene expression in mammalian organs

Gene duplications generate genomic raw material that allows the emergence of novel functions, likely facilitating adaptive evolutionary innovations. However, global assessments of the functional and evolutionary relevance of duplicate genes in mammals were until recently limited by the lack of appropriate comparative data. Here, we report a large-scale study of the expression evolution of DNA-based functional gene duplicates in three major mammalian lineages (placental mammals, marsupials, egg-laying monotremes) and birds, on the basis of RNA sequencing (RNA-seq) data from nine species and eight organs. We observe dynamic changes in tissue expression preference of paralogs with different duplication ages, suggesting differential contribution of paralogs to specific organ functions during vertebrate evolution. Specifically, we show that paralogs that emerged in the common ancestor of bony vertebrates are enriched for genes with brain-specific expression and provide evidence for differential forces underlying the preferential emergence of young testis- and liver-specific expressed genes. Further analyses uncovered that the overall spatial expression profiles of gene families tend to be conserved, with several exceptions of pronounced tissue specificity shifts among lineage-specific gene family expansions. Finally, we trace new lineage-specific genes that may have contributed to the specific biology of mammalian organs, including the little-studied placenta. Overall, our study provides novel and taxonomically broad evidence for the differential contribution of duplicate genes to tissue-specific transcriptomes and for their importance for the phenotypic evolution of vertebrates.

Reconstruction and evolutionary history of eutherian chromosomes

Whole-genome assemblies of 19 placental mammals and two outgroup species were used to reconstruct the order and orientation of syntenic fragments in chromosomes of the eutherian ancestor and six other descendant ancestors leading to human. For ancestral chromosome reconstructions, we developed an algorithm (DESCHRAMBLER) that probabilistically determines the adjacencies of syntenic fragments using chromosome-scale and fragmented genome assemblies. The reconstructed chromosomes of the eutherian, boreoeutherian, and euarchontoglires ancestor each included >80% of the entire length of the human genome, whereas reconstructed chromosomes of the most recent common ancestor of simians, catarrhini, great apes, and humans and chimpanzees included >90% of human genome sequence. These high-coverage reconstructions permitted reliable identification of chromosomal rearrangements over ∼105 My of eutherian evolution. Orangutan was found to have eight chromosomes that were completely conserved in homologous sequence order and orientation with the eutherian ancestor, the largest number for any species. Ruminant artiodactyls had the highest frequency of intrachromosomal rearrangements, and interchromosomal rearrangements dominated in murid rodents. A total of 162 chromosomal breakpoints in evolution of the eutherian ancestral genome to the human genome were identified; however, the rate of rearrangements was significantly lower (0.80/My) during the first ∼60 My of eutherian evolution, then increased to greater than 2.0/My along the five primate lineages studied. Our results significantly expand knowledge of eutherian genome evolution and will facilitate greater understanding of the role of chromosome rearrangements in adaptation, speciation, and the etiology of inherited and spontaneously occurring diseases.

The C2H2-ZF transcription factor Zfp335 recognizes two consensus motifs using separate zinc finger arrays

Here, Han et al. show that transcription factor Zfp335 binds DNA and drives transcription via recognition of two distinct consensus motifs by separate ZF clusters and identify the specific motif interaction disrupted by the mutation R1092W. This study presents Zfp335 as a model for understanding how C2H2-ZF TFs may use multiple recognition motifs to control gene expression.

The complexities of DNA recognition by transcription factors (TFs) with multiple Cys2–His2 zinc fingers (C2H2-ZFs) remain poorly studied. We previously reported a mutation (R1092W) in the C2H2-ZF TF Zfp335 that led to selective loss of binding at a subset of targets, although the basis for this effect was unclear. We show that Zfp335 binds DNA and drives transcription via recognition of two distinct consensus motifs by separate ZF clusters and identify the specific motif interaction disrupted by R1092W. Our work presents Zfp335 as a model for understanding how C2H2-ZF TFs may use multiple recognition motifs to control gene expression.

Methylomics of breast cancer: Seeking epimarkers in peripheral blood of young subjects

Critical roles of epigenomic alterations in the pathogenesis of breast cancer have recently seized great attentions toward finding epimarkers in either non-invasive or semi-non-invasive samples as well as peripheral blood. In this way, methylated DNA immunoprecipitation microarray (MeDIP-chip) was performed on DNA samples isolated from white blood cells of 30 breast cancer patients compared to 30 healthy controls. A total of 1799 differentially methylated regions were identified including SLC6A3, Rab40C, ZNF584, and FOXD3 whose significant methylation differences were confirmed in breast cancer patients through quantitative real-time polymerase chain reaction. Hypermethylation of APC, HDAC1, and GSK1 genes has been previously reported in more than one study on tissue samples of breast cancer. Methylation of those aforementioned genes in white blood cells of our young patients not only relies on their importance in breast cancer pathogenesis but also may highlight their potential as early epimarkers that makes further assessments necessary in large cohort studies.

The mouse genome displays highly dynamic populations of KRAB-zinc finger protein genes and related genetic units

KRAB-containing poly-zinc finger proteins (KZFPs) constitute the largest family of transcription factors encoded by mammalian genomes, and growing evidence indicates that they fulfill functions critical to both embryonic development and maintenance of adult homeostasis. KZFP genes underwent broad and independent waves of expansion in many higher vertebrates lineages, yet comprehensive studies of members harbored by a given species are scarce. Here we present a thorough analysis of KZFP genes and related units in the murine genome. We first identified about twice as many elements than previously annotated as either KZFP genes or pseudogenes, notably by assigning to this family an entity formerly considered as a large group of Satellite repeats. We then could delineate an organization in clusters distributed throughout the genome, with signs of recombination, translocation, duplication and seeding of new sites by retrotransposition of KZFP genes and related genetic units (KZFP/rGUs). Moreover, we harvested evidence indicating that closely related paralogs had evolved through both drifting and shifting of sequences encoding for zinc finger arrays. Finally, we could demonstrate that the KAP1-SETDB1 repressor complex tames the expression of KZFP/rGUs within clusters, yet that the primary targets of this regulation are not the KZFP/rGUs themselves but enhancers contained in neighboring endogenous retroelements and that, underneath, KZFPs conserve highly individualized patterns of expression.

Discovery of a novel oocyte-specific Krüppel-associated box domain-containing zinc finger protein required for early embryogenesis in cattle

Zinc finger (ZNF) transcription factors interact with DNA through zinc finger motifs and play important roles in a variety of cellular functions including cell growth, proliferation, development, apoptosis, and intracellular signal transduction. One-third of ZNF proteins in metazoans contain a highly conserved N-terminal motif known as the Krüppel-associated box (KRAB) domain, which acts as a potent, DNA-binding dependent transcriptional repression module. Analysis of RNA-Seq data generated from a bovine oocyte cDNA library identified a novel transcript, which encodes a KRAB-containing ZNF transcription factor (named ZNFO). Characterization of ZNFO mRNA expression revealed that it is exclusively expressed in bovine oocytes and early embryos. A GFP reporter assay demonstrated that ZNFO protein localizes specifically to the nucleus, supporting its role in transcriptional regulation. To test the role of ZNFO in early embryonic development, zygotes were generated by in vitro maturation and fertilization of oocytes, and injected with small interfering RNA (siRNA) designed to knockdown ZNFO. Cleavage rates were not affected by ZNFO siRNA injection. However, embryonic development to 8- to 16-cell stage and blastocyst stage was significantly reduced relative to the uninjected and negative control siRNA-injected embryos. Further, interaction of ZNFO with the highly conserved co-factor, KRAB-associated protein-1 (KAP1), was demonstrated, and evidence supporting transcriptional repression by ZNFO was demonstrated using a GAL4-luciferase reporter system. Results of described studies demonstrate that ZNFO is a maternally-derived oocyte-specific nuclear factor required for early embryonic development in cattle, presumably functioning by repressing transcription.

Genome-wide Regulatory Roles of the C2H2-type Zinc Finger Protein ZNF764 on the Glucocorticoid Receptor

The C2H2-type zinc finger protein ZNF764 acts as an enhancer for several steroid hormone receptors, and haploinsufficiency of this gene may be responsible for tissue resistance to multiple steroid hormones including glucocorticoids observed in a patient with 16p11.2 microdeletion. We examined genome-wide regulatory actions of ZNF764 on the glucocorticoid receptor (GR) in HeLa cells as a model system. ZNF764- and GR-binding sites demonstrated similar distribution in various genomic features. They positioned predominantly around 50-500 kbs from the transcription start sites of their nearby genes, and were closely localized with each other, overlapping in ~37% of them. ZNF764 demonstrated differential on/off effects on GR-binding and subsequent mRNA expression: some genes were highly dependent on the presence/absence of ZNF764, but others were not. Pathway analysis revealed that these 3 gene groups were involved in distinct cellular activities. ZNF764 physically interacted with GR at ligand-binding domain through its KRAB domain, and both its physical interaction to GR and zinc finger domain appear to be required for ZNF764 to regulate GR transcriptional activity. Thus, ZNF764 is a cofactor directing GR transcriptional activity toward specific biologic pathways by changing GR binding and transcriptional activity on the glucocorticoid-responsive genes.

The Complex Role of the ZNF224 Transcription Factor in Cancer

ZNF224 is a member of the Kruppel-associated box zinc finger proteins (KRAB-ZFPs) family. It was originally identified as a transcriptional repressor involved in gene-specific silencing through the recruitment of the corepressor KAP1, chromatin-modifying activities, and the arginine methyltransferase PRMT5 on the promoter of its target genes. Recent findings indicate that ZNF224 can behave both as a tumor suppressor or an oncogene in different human cancers. The transcriptional regulatory properties of ZNF224 in these systems appear to be complex and influenced by specific sets of interactors. ZNF224 can also act as a transcription cofactor for other DNA-binding proteins. A role for ZNF224 in transcriptional activation has also emerged. Here, we review the state of the literature supporting both roles of ZNF224 in cancer. We also examine the functional activity of ZNF224 as a transcription factor and the influence of protein partners on its dual behavior. Increasing information on the mechanism through which ZNF224 can operate could lead to the identification of agents capable of modulating ZNF224 function, thus potentially paving the way to new therapeutic strategies for treatment of cancer.

Independent evolution of genomic characters during major metazoan transitions

Metazoan evolution encompasses a vast evolutionary time scale spanning over 600 million years. Our ability to infer ancestral metazoan characters, both morphological and functional, is limited by our understanding of the nature and evolutionary dynamics of the underlying regulatory networks. Increasing coverage of metazoan genomes enables us to identify the evolutionary changes of the relevant genomic characters such as the loss or gain of coding sequences, gene duplications, micro- and macro-synteny, and non-coding element evolution in different lineages. In this review we describe recent advances in our understanding of ancestral metazoan coding and non-coding features, as deduced from genomic comparisons. Some genomic changes such as innovations in gene and linkage content occur at different rates across metazoan clades, suggesting some level of independence among genomic characters. While their contribution to biological innovation remains largely unclear, we review recent literature about certain genomic changes that do correlate with changes to specific developmental pathways and metazoan innovations. In particular, we discuss the origins of the recently described pharyngeal cluster which is conserved across deuterostome genomes, and highlight different genomic features that have contributed to the evolution of this group. We also assess our current capacity to infer ancestral metazoan states from gene models and comparative genomics tools and elaborate on the future directions of metazoan comparative genomics relevant to evo-devo studies.

Multiparameter functional diversity of human C2H2 zinc finger proteins

C2H2 zinc finger proteins represent the largest and most enigmatic class of human transcription factors. Their C2H2-ZF arrays are highly variable, indicating that most will have unique DNA binding motifs. However, most of the binding motifs have not been directly determined. In addition, little is known about whether or how these proteins regulate transcription. Most of the ∼700 human C2H2-ZF proteins also contain at least one KRAB, SCAN, BTB, or SET domain, suggesting that they may have common interacting partners and/or effector functions. Here, we report a multifaceted functional analysis of 131 human C2H2-ZF proteins, encompassing DNA binding sites, interacting proteins, and transcriptional response to genetic perturbation. We confirm the expected diversity in DNA binding motifs and genomic binding sites, and provide motif models for 78 previously uncharacterized C2H2-ZF proteins, most of which are unique. Surprisingly, the diversity in protein-protein interactions is nearly as high as diversity in DNA binding motifs: Most C2H2-ZF proteins interact with a unique spectrum of co-activators and co-repressors. Thus, multiparameter diversification likely underlies the evolutionary success of this large class of human proteins.

The pioneering role of PRDM9 indel mutations in tarsier evolution

PRDM9 is currently the sole speciation gene found in vertebrates causing hybrid sterility probably due to incompatible alleles. Its role in defining the double strand break loci during the meiotic prophase I is crucial for proper chromosome segregation. Therefore, the rapid turnover of the loci determining zinc finger array seems to be causative for incompatibilities. We here investigated the zinc finger domain-containing exon of PRDM9 in 23 tarsiers. Tarsiers, the most basal extant haplorhine primates, exhibit two frameshifting indels at the 5'-end of the array. The first mutation event interrupts the reading frame and function while the second compensates both. The fixation of this allele variant in tarsiers led to hypothesize that de- and reactivation of the zinc finger domain drove the speciation in early haplorhine or tarsiiform primates. Moreover, the high allelic diversity within Tarsius points to multiple effects of genetic drift reflecting their phylogeographic history since the Miocene.

The Transcriptional Repressive Activity of KRAB Zinc Finger Proteins Does Not Correlate with Their Ability to Recruit TRIM28

KRAB domain Zinc finger proteins are one of the most abundant families of transcriptional regulators in higher vertebrates. The prevailing view is that KRAB domain proteins function as potent transcriptional repressors by recruiting TRIM28 and promoting heterochromatin spreading. However, the extent to which all KRAB domain proteins are TRIM28-dependent transcriptional repressors is currently unclear. Our studies on mouse ZFP568 revealed that TRIM28 recruitment by KRAB domain proteins is not sufficient to warrant transcriptional repressive activity. By using luciferase reporter assays and yeast two-hybrid experiments, we tested the ability of ZFP568 and other mouse KRAB domain proteins to repress transcription and bind TRIM28. We found that some mouse KRAB domain proteins are poor transcriptional repressors despite their ability to recruit TRIM28, while others showed strong KRAB-dependent transcriptional repression, but no TRIM28 binding. Together, our results show that the transcriptional repressive activity of KRAB-ZNF proteins does not correlate with their ability to recruit TRIM28, and provide evidence that KRAB domains can regulate transcription in a TRIM28-independent fashion. Our findings challenge the current understanding of the molecular mechanisms used by KRAB domain proteins to control gene expression and highlight that a high percentage of KRAB domain proteins in the mouse genome differ from the consensus KRAB sequence at amino acid residues that are critical for TRIM28 binding and/or repressive activity.

Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates

Since their discovery, a growing body of evidence has emerged demonstrating that transposable elements are important drivers of species diversity. These mobile elements exhibit a great variety in structure, size and mechanisms of transposition, making them important putative actors in organism evolution. The vertebrates represent a highly diverse and successful lineage that has adapted to a wide range of different environments. These animals also possess a rich repertoire of transposable elements, with highly diverse content between lineages and even between species. Here, we review how transposable elements are driving genomic diversity and lineage-specific innovation within vertebrates. We discuss the large differences in TE content between different vertebrate groups and then go on to look at how they affect organisms at a variety of levels: from the structure of chromosomes to their involvement in the regulation of gene expression, as well as in the formation and evolution of non-coding RNAs and protein-coding genes. In the process of doing this, we highlight how transposable elements have been involved in the evolution of some of the key innovations observed within the vertebrate lineage, driving the group's diversity and success.

Transposable Elements and Their KRAB-ZFP Controllers Regulate Gene Expression in Adult Tissues

KRAB-containing zinc finger proteins (KRAB-ZFPs) are early embryonic controllers of transposable elements (TEs), which they repress with their cofactor KAP1 through histone and DNA methylation, a process thought to result in irreversible silencing. Using a target-centered functional screen, we matched murine TEs with their cognate KRAB-ZFP. We found the paralogs ZFP932 and Gm15446 to bind overlapping but distinguishable subsets of ERVK (endogenous retrovirus K), repress these elements in embryonic stem cells, and regulate secondarily the expression of neighboring genes. Most importantly, we uncovered that these KRAB-ZFPs and KAP1 control TEs in adult tissues, in cell culture and in vivo, where they partner up to modulate cellular genes. Therefore, TEs and KRAB-ZFPs establish transcriptional networks that likely regulate not only development but also many physiological events. Given the high degree of species specificity of TEs and KRAB-ZFPs, these results have important implications for understanding the biology of higher vertebrates, including humans.

Diverse and pervasive subcellular distributions for both coding and long noncoding RNAs

Here, Wilk et al. examined ∼8000 mRNA transcripts throughout Drosophila embryogenesis. They found that almost all RNAs, both coding and noncoding RNAs, are subcellularly localized at some stage of development, thus providing an important resource for functional gene analysis.

In a previous analysis of 2300 mRNAs via whole-mount fluorescent in situ hybridization in cellularizing Drosophila embryos, we found that 70% of the transcripts exhibited some form of subcellular localization. To see whether this prevalence is unique to early Drosophila embryos, we examined ∼8000 transcripts over the full course of embryogenesis and ∼800 transcripts in late third instar larval tissues. The numbers and varieties of new subcellular localization patterns are both striking and revealing. In the much larger cells of the third instar larva, virtually all transcripts observed showed subcellular localization in at least one tissue. We also examined the prevalence and variety of localization mechanisms for >100 long noncoding RNAs. All of these were also found to be expressed and subcellularly localized. Thus, subcellular RNA localization appears to be the norm rather than the exception for both coding and noncoding RNAs. These results, which have been annotated and made available on a recompiled database, provide a rich and unique resource for functional gene analyses, some examples of which are provided.

ATRX binds to atypical chromatin domains at the 3' exons of zinc finger genes to preserve H3K9me3 enrichment

ATRX is a SWI/SNF chromatin remodeler proposed to govern genomic stability through the regulation of repetitive sequences, such as rDNA, retrotransposons, and pericentromeric and telomeric repeats. However, few direct ATRX target genes have been identified and high-throughput genomic approaches are currently lacking for ATRX. Here we present a comprehensive ChIP-sequencing study of ATRX in multiple human cell lines, in which we identify the 3' exons of zinc finger genes (ZNFs) as a new class of ATRX targets. These 3' exonic regions encode the zinc finger motifs, which can range from 1-40 copies per ZNF gene and share large stretches of sequence similarity. These regions often contain an atypical chromatin signature: they are transcriptionally active, contain high levels of H3K36me3, and are paradoxically enriched in H3K9me3. We find that these ZNF 3' exons are co-occupied by SETDB1, TRIM28, and ZNF274, which form a complex with ATRX. CRISPR/Cas9-mediated loss-of-function studies demonstrate (i) a reduction of H3K9me3 at the ZNF 3' exons in the absence of ATRX and ZNF274 and, (ii) H3K9me3 levels at atypical chromatin regions are particularly sensitive to ATRX loss compared to other H3K9me3-occupied regions. As a consequence of ATRX or ZNF274 depletion, cells with reduced levels of H3K9me3 show increased levels of DNA damage, suggesting that ATRX binds to the 3' exons of ZNFs to maintain their genomic stability through preservation of H3K9me3.

Charactering the ZFAND3 gene mapped in the sex-determining locus in hybrid tilapia (Oreochromis spp.)

Zinc finger AN1-type domain 3 (ZFAND3) is essential for spermatogenesis in mice. However, its function in teleosts remains unclear. In this study, we characterized the ZFAND3 gene (termed as OsZFAND3) in an important food fish, tilapia. The OsZFAND3 cDNA sequence is 1,050 bp in length, containing an ORF of 615 bp, which encodes a putative peptide of 204 amino acid residues. Quantitative real-time PCR revealed that the OsZFAND3 transcripts were exclusively expressed in the testis and ovary. In situ hybridization showed that the high expression of OsZFAND3 transcripts was predominantly localized in the spermatocyte and spermatid. These results suggest that OsZFAND3 is involved in male germ cell maturation. Three single nucleotide polymorphisms (SNPs) were detected in the introns of OsZFAND3. The OsZFAND3 gene was mapped in the sex-determining locus on linkage group 1 (LG1). The three SNPs in the OsZFAND3 gene were strictly associated with sex phenotype, suggesting that the OsZFAND3 gene is tightly linked to the sex-determining locus. Our study provides new insights into the functions of the OsZFAND3 gene in tilapia and a foundation for further detailed analysis of the OsZFAND3 gene in sex determination and differentiation.

Combgap contributes to recruitment of Polycomb group proteins in Drosophila

Polycomb group (PcG) proteins are responsible for maintaining the silenced transcriptional state of many developmentally regulated genes. PcG proteins are organized into multiprotein complexes that are recruited to DNA via cis-acting elements known as "Polycomb response elements" (PREs). In Drosophila, PREs consist of binding sites for many different DNA-binding proteins, some known and others unknown. Identification of these DNA-binding proteins is crucial to understanding the mechanism of PcG recruitment to PREs. We report here the identification of Combgap (Cg), a sequence-specific DNA-binding protein that is involved in recruitment of PcG proteins. Cg can bind directly to PREs via GTGT motifs and colocalizes with the PcG proteins Pleiohomeotic (Pho) and Polyhomeotic (Ph) at the majority of PREs in the genome. In addition, Cg colocalizes with Ph at a number of targets independent of Pho. Loss of Cg leads to decreased recruitment of Ph at only a subset of sites; some of these sites are binding sites for other Polycomb repressive complex 1 (PRC1) components, others are not. Our data suggest that Cg can recruit Ph in the absence of PRC1 and illustrate the diversity and redundancy of PcG protein recruitment mechanisms.

Effect of the linkers between the zinc fingers in zinc finger protein 809 on gene silencing and nuclear localization

Zinc finger protein 809 (ZFP809) belongs to the Kruppel-associated box-containing zinc finger protein (KRAB-ZFP) family and functions in repressing the expression of Moloney murine leukemia virus (MoMLV). ZFP809 binds to the primer-binding site (PBS)located downstream of the MoMLV-long terminal repeat (LTR) and induces epigenetic modifications at integration sites, such as repressive histone modifications and de novo DNA methylation. KRAB-ZFPs contain consensus TGEKP linkers between C2H2 zinc fingers. The phosphorylation of threonine residues within linkers leads to the inactivation of zinc finger binding to target sequences. ZFP809 also contains consensus linkers between zinc fingers. However, the function of ZFP809 linkers remains unknown. In the present study, we constructed ZFP809 proteins containing mutated linkers and examined their ability to silence transgene expression driven by MLV, binding ability to MLV PBS, and cellular localization. The results of the present study revealed that the linkers affected the ability of ZFP809 to silence transgene expression. Furthermore, this effect could be partly attributed to changes in the localization of ZFP809 proteins containing mutated linkers. Further characterization of ZFP809 linkers is required for understanding the functions and features of KRAB-ZFP-containing linkers.

An epigenetic switch ensures transposon repression upon dynamic loss of DNA methylation in embryonic stem cells

DNA methylation is extensively remodeled during mammalian gametogenesis and embryogenesis. Most transposons become hypomethylated, raising the question of their regulation in the absence of DNA methylation. To reproduce a rapid and extensive demethylation, we subjected mouse ES cells to chemically defined hypomethylating culture conditions. Surprisingly, we observed two phases of transposon regulation. After an initial burst of de-repression, various transposon families were efficiently re-silenced. This was accompanied by a reconfiguration of the repressive chromatin landscape: while H3K9me3 was stable, H3K9me2 globally disappeared and H3K27me3 accumulated at transposons. Interestingly, we observed that H3K9me3 and H3K27me3 occupy different transposon families or different territories within the same family, defining three functional categories of adaptive chromatin responses to DNA methylation loss. Our work highlights that H3K9me3 and, most importantly, polycomb-mediated H3K27me3 chromatin pathways can secure the control of a large spectrum of transposons in periods of intense DNA methylation change, ensuring longstanding genome stability.

Transposable Elements, Polydactyl Proteins, and the Genesis of Human-Specific Transcription Networks

Transposable elements (TEs) may account for up to two-thirds of the human genome, and as genomic threats they are subjected to epigenetic control mechanisms engaged from the earliest stages of embryonic development. We previously determined that an important component of this process is the sequence-specific recognition of TEs by KRAB (Krüppel-associated box)-containing zinc-finger proteins (KRAB-ZFPs), a large family of tetrapod-restricted transcription factors that act by recruiting inducers of heterochromatin formation and DNA methylation. We further showed that KRAB-ZFPs and their cofactor KAP1 exert a marked influence on the transcription dynamics of embryonic stem cells via their docking of repressor complexes at TE-contained regulatory sequences. It is generally held that, beyond this early embryonic period, TEs become permanently silenced, and that the evolutionary selection of KRAB-ZFPs and other TE controllers is the result of a simple evolutionary arms race between the host and these genetics invaders. Here, I discuss recent evidence that invalidates this dual assumption and instead suggests that KRAB-ZFPs are the instruments of a massive enterprise of TE domestication, whereby transposon-based regulatory sequences and their cellular ligands establish species-specific transcription regulation networks that influence multiple aspects of human development and physiology.

A Large-Scale Functional Screen to Identify Epigenetic Repressors of Retrotransposon Expression

Deposition of epigenetic marks is an important layer of the transcriptional control of retrotransposons, especially during early embryogenesis. Krüppel-associated box domain zinc finger proteins (KRAB-ZFPs) are one of the largest families of transcription factors, and collectively partake in this process by tethering to thousands of retroelement-containing genomic loci their cofactor KAP1, which acts as a scaffold for a heterochromatin-inducing machinery. However, while the sequence-specific DNA binding potential of the poly-zinc finger-containing KRAB-ZFPs is recognized, very few members of the family have been assigned specific targets. In this chapter, we describe a large-scale functional screen to identify the retroelements bound by individual murine KRAB-ZFPs. Our method is based on the automated transfection of a library of mouse KRAB-ZFP-containing vectors into 293T cells modified to express GFP from a PGK promoter harboring in its immediate vicinity a KAP1-recruiting retroelement-derived sequence. Analysis is then performed by plate reader and flow cytometry fluorescence readout. Such large-scale DNA-centered functional approach can not only help to identify the trans-acting factors responsible for silencing retrotransposons, but also serve as a model for dissecting the transcriptional networks influenced by retroelement-derived cis-acting sequences.

The evolution of gene expression and binding specificity of the largest transcription factor family in primates

The KRAB-containing zinc finger (KRAB-ZF) proteins represent the largest family of transcription factors (TFs) in humans, yet for the great majority, their function and specific genomic target remain unknown. However, it has been shown that a large fraction of these genes arose from segmental duplications, and that they have expanded in gene and zinc finger number throughout vertebrate evolution. To determine whether this expansion is linked to selective pressures acting on different domains, we have manually curated all KRAB-ZF genes present in the human genome together with their orthologous genes in three closely related species and assessed the evolutionary forces acting at the sequence level as well as on their expression profiles. We provide evidence that KRAB-ZFs can be separated into two categories according to the polymorphism present in their DNA-contacting residues. Those carrying a nonsynonymous single nucleotide polymorphism (SNP) in their DNA-contacting amino acids exhibit significantly reduced expression in all tissues, have emerged in a recent lineage, and seem to be less strongly constrained evolutionarily than those without such a polymorphism. This work provides evidence for a link between age of the TF, as well as polymorphism in their DNA-contacting residues and expression levels-both of which may be jointly affected by selection.

Spotting the enemy within: Targeted silencing of foreign DNA in mammalian genomes by the Krüppel-associated box zinc finger protein family

Tandem C2H2-type zinc finger proteins (ZFPs) constitute the largest transcription factor family in animals. Tandem-ZFPs bind DNA in a sequence-specific manner through arrays of multiple zinc finger domains that allow high flexibility and specificity in target recognition. In tetrapods, a large proportion of tandem-ZFPs contain Krüppel-associated-box (KRAB) repression domains, which are able to induce epigenetic silencing through the KAP1 corepressor. The KRAB-ZFP family continuously amplified in tetrapods through segmental gene duplications, often accompanied by deletions, duplications, and mutations of the zinc finger domains. As a result, tetrapod genomes contain unique sets of KRAB-ZFP genes, consisting of ancient and recently evolved family members. Although several hundred human and mouse KRAB-ZFPs have been identified or predicted, the biological functions of most KRAB-ZFP family members have gone unexplored. Furthermore, the evolutionary forces driving the extraordinary KRAB-ZFP expansion and diversification have remained mysterious for decades. In this review, we highlight recent studies that associate KRAB-ZFPs with the repression of parasitic DNA elements in the mammalian germ line and discuss the hypothesis that the KRAB-ZFP family primarily evolved as an adaptive genomic surveillance system against foreign DNA. Finally, we comment on the computational, genetic, and biochemical challenges of studying KRAB-ZFPs and attempt to predict how these challenges may be soon overcome.

Functional Domains of ZFP809 Essential for Nuclear Localization and Gene Silencing

Zinc finger protein 809 (ZFP809) is a member of the Kruppel-associated box-containing zinc finger protein (KRAB-ZFP) family, and is highly expressed in mouse immature cells. ZFP809 is known to inhibit the expression of transduced genes driven by Moloney murine leukemia virus (MoMLV)-typed retroviral vectors by binding to the primer binding site (PBS) located downstream of the MLV-long terminal repeat (LTR) of the vectors and recruiting protein complexes that introduce epigenetic silencing marks such as histone modifications and DNA methylation at the MLV-LTR. However, it remains undetermined what domains of ZFP809 among the KRAB domain at N-terminus and the seven zinc fingers are critical for gene silencing. In this study, we assessed subcellular localization, gene silencing ability, and binding ability to the PBS of a series of truncated and mutated ZFP809 proteins. We revealed the essential role of the KRAB A box for all functions assessed, together with the accessory roles of a subset of zinc fingers. Our data also suggest that interaction between KAP1 and the KRAB A box of ZFP809 is critical in KAP1-dependent control of gene silencing for ZFP809 targets.

The developmental control of transposable elements and the evolution of higher species

Transposable elements (TEs) account for at least 50% of the human genome. They constitute essential motors of evolution through their ability to modify genomic architecture, mutate genes and regulate gene expression. Accordingly, TEs are subject to tight epigenetic control during the earliest phases of embryonic development via histone and DNA methylation. Key to this process is recognition by sequence-specific RNA- and protein-based repressors. Collectively, these mediators are responsible for silencing a very broad range of TEs in an evolutionarily dynamic fashion. As a consequence, mobile elements and their controllers exert a marked influence on transcriptional networks in embryonic stem cells and a variety of adult tissues. The emerging picture is not that of a simple arms race but rather of a massive and sophisticated enterprise of TE domestication for the evolutionary benefit of the host.

A structural approach reveals how neighbouring C2H2 zinc fingers influence DNA binding specificity

Development of an accurate protein–DNA recognition code that can predict DNA specificity from protein sequence is a central problem in biology. C₂H₂ zinc fingers constitute by far the largest family of DNA binding domains and their binding specificity has been studied intensively. However, despite decades of research, accurate prediction of DNA specificity remains elusive. A major obstacle is thought to be the inability of current methods to account for the influence of neighbouring domains. Here we show that this problem can be addressed using a structural approach: we build structural models for all C₂H₂-ZF–DNA complexes with known binding motifs and find six distinct binding modes. Each mode changes the orientation of specificity residues with respect to the DNA, thereby modulating base preference. Most importantly, the structural analysis shows that residues at the domain interface strongly and predictably influence the binding mode, and hence specificity. Accounting for predicted binding mode significantly improves prediction accuracy of predicted motifs. This new insight into the fundamental behaviour of C₂H₂-ZFs has implications for both improving the prediction of natural zinc finger-binding sites, and for prioritizing further experiments to complete the code. It also provides a new design feature for zinc finger engineering.

Human Genes Encoding Transcription Factors and Chromatin-Modifying Proteins Have Low Levels of Promoter Polymorphism: A Study of 1000 Genomes Project Data

The expression level of each gene is controlled by its regulatory regions, which determine the precise regulation in a tissue-specific manner, according to the developmental stage of the body and the necessity of a response to external stimuli. Nucleotide substitutions in regulatory gene regions may modify the affinity of transcription factors to their specific DNA binding sites, affecting the transcription rates of genes. In our previous research, we found that genes controlling the sensory perception of smell and genes involved in antigen processing and presentation were overrepresented significantly among genes with high SNP contents in their promoter regions. The goal of our study was to reveal functional features of human genes containing extremely small numbers of SNPs in promoter regions. Two functional groups were found to be overrepresented among genes whose promoters did not contain SNPs: (1) genes involved in gene-specific transcription and (2) genes controlling chromatin organization. We revealed that the 5′-regulatory regions of genes encoding transcription factors and chromatin-modifying proteins were characterized by reduced genetic variability. One important exception from this rule refers to genes encoding transcription factors with zinc-coordinating DNA-binding domains (DBDs), which underwent extensive expansion in vertebrates, particularly, in primate evolution. Hence, we obtained new evidence for evolutionary forces shaping variability in 5′-regulatory regions of genes.

Transcriptomic changes in Ca²⁺-depleted cells: Role of elevated intracellular [Na⁺]/[K⁺] ratio

Previously, we reported that Ca(2+) depletion increased permeability of the plasma membrane for Na(+). This study examined the relative impact of [Na(+)]i/[K(+)]i-mediated signaling on transcriptomic changes in cultured vascular smooth muscle cells from rat aorta (VSMC) subjected to Ca(2+)-depletion by extra-(EGTA) and intracellular (BAPTA-AM) Ca(2+) chelators. Na(+),K(+)-ATPase inhibition in K(+)-free medium during 3 h led to elevation of [Na(+)]i and attenuation of [K(+)]i by ∼7- and 10-fold, whereas Ca(2+)-depletion resulted in alteration of these parameters by ∼3- and 2-fold, respectively. Augmented VSMC permeability for Na(+) and elevation of the [Na(+)]i/[K(+)]i ratio was triggered by addition to Ca(2+)-free medium 50 μM EGTA and was not affected by 10 μM BAPTA-AM. Na(+),K(+)-ATPase inhibition and Ca(2+)-depletion changed expression of 3677 and 4610 mRNA transcripts, respectively. We found highly significant (p<10(-12)) positive (R(2)>0.51) correlation between levels of expression of 2071 transcripts whose expression was affected by both stimuli. Among genes whose expression in Ca(2+)-depleted cells was augmented by more than 7-fold we noted cyclic AMP-dependent transcription factor Atf3, early growth response protein Egr1 and nuclear receptor subfamily 4, group A member Nr4a1. Dissipation of transmembrane gradients of monovalent cations in high-K(+), low-Na(+)-medium abolished the increments of the [Na(+)]i/[K(+)]i ratio as well as the augmented expression of these genes triggered by incubation of VSMC in EGTA containing medium. Thus, our results demonstrate, for the first time, that robust transcriptomic changes triggered by Ca(2+)-depletion in the presence of extracellular Ca(2+)-chelators are at least partially mediated by elevation of the [Na(+)]i/[K(+)]i ratio and activation of Ca(2+)i-independent, [Na(+)]i/[K(+)]i-mediated mechanism of excitation-transcription coupling. These results shad a new light on analysis of data obtained in cells subjected to long-term exposure to Ca(2+) chelators.

Retrotransposons shape species-specific embryonic stem cell gene expression

Over half of our genome is composed of retrotransposons, which are mobile elements that can readily amplify their copy number by replicating through an RNA intermediate. Most of these elements are no longer mobile but still contain regulatory sequences that can serve as promoters, enhancers or repressors for cellular genes. Despite dominating our genetic content, little is known about the precise functions of retrotransposons, which include both endogenous retroviruses (ERVs) and non-LTR elements like long interspersed nuclear element 1 (LINE-1). However, a few recent cutting-edge publications have illustrated how retrotransposons shape species-specific stem cell gene expression by two opposing mechanisms, involving their recruitment of stem cell-enriched transcription factors (TFs): firstly, they can activate expression of genes linked to naïve pluripotency, and secondly, they can induce repression of proximal genes. The paradox that different retrotransposons are active or silent in embryonic stem cells (ESCs) can be explained by differences between retrotransposon families, between individual copies within the same family, and between subpopulations of ESCs. Since they have coevolved with their host genomes, some of them have been co-opted to perform species-specific beneficial functions, while others have been implicated in genetic disease. In this review, we will discuss retrotransposon functions in ESCs, focusing on recent mechanistic advances of how HERV-H has been adopted to preserve human naïve pluripotency and how particular LINE-1, SVA and ERV family members recruit species-specific transcriptional repressors. This review highlights the fine balance between activation and repression of retrotransposons that exists to harness their ability to drive evolution, while minimizing the risk they pose to genome integrity.

Inflammation and Nutritional Science for Programs/Policies and Interpretation of Research Evidence (INSPIRE)

An increasing recognition has emerged of the complexities of the global health agenda—specifically, the collision of infections and noncommunicable diseases and the dual burden of over- and undernutrition. Of particular practical concern are both 1) the need for a better understanding of the bidirectional relations between nutritional status and the development and function of the immune and inflammatory response and 2) the specific impact of the inflammatory response on the selection, use, and interpretation of nutrient biomarkers. The goal of the Inflammation and Nutritional Science for Programs/Policies and Interpretation of Research Evidence (INSPIRE) is to provide guidance for those users represented by the global food and nutrition enterprise. These include researchers (bench and clinical), clinicians providing care/treatment, those developing and evaluating programs/interventions at scale, and those responsible for generating evidence-based policy. The INSPIRE process included convening 5 thematic working groups (WGs) charged with developing summary reports around the following issues: 1) basic overview of the interactions between nutrition, immune function, and the inflammatory response; 2) examination of the evidence regarding the impact of nutrition on immune function and inflammation; 3) evaluation of the impact of inflammation and clinical conditions (acute and chronic) on nutrition; 4) examination of existing and potential new approaches to account for the impact of inflammation on biomarker interpretation and use; and 5) the presentation of new approaches to the study of these relations. Each WG was tasked with synthesizing a summary of the evidence for each of these topics and delineating the remaining gaps in our knowledge. This review consists of a summary of the INSPIRE workshop and the WG deliberations.

The KRAB zinc finger protein ZFP809 is required to initiate epigenetic silencing of endogenous retroviruses

Endogenous retroviruses (ERVs) are epigenetically silenced during development, yet the cellular factors recognizing ERVs in a sequence-specific manner remain elusive. Wolf et al. find that ZFP809 initiates the silencing of ERVs in a sequence-specific manner via recruitment of heterochromatin-inducing complexes. ERV reactivation is accompanied by an epigenetic shift from repressive to active histone modifications. ZFP809 is required to initiate ERV silencing during embryonic development but becomes largely dispensable in somatic tissues.

Retroviruses have been invading mammalian germlines for millions of years, accumulating in the form of endogenous retroviruses (ERVs) that account for nearly one-tenth of the mouse and human genomes. ERVs are epigenetically silenced during development, yet the cellular factors recognizing ERVs in a sequence-specific manner remain elusive. Here we demonstrate that ZFP809, a member of the Krüppel-associated box zinc finger protein (KRAB-ZFP) family, initiates the silencing of ERVs in a sequence-specific manner via recruitment of heterochromatin-inducing complexes. ZFP809 knockout mice display highly elevated levels of ZFP809-targeted ERVs in somatic tissues. ERV reactivation is accompanied by an epigenetic shift from repressive to active histone modifications but only slight destabilization of DNA methylation. Importantly, using conditional alleles and rescue experiments, we demonstrate that ZFP809 is required to initiate ERV silencing during embryonic development but becomes largely dispensable in somatic tissues. Finally, we show that the DNA-binding specificity of ZFP809 is evolutionarily conserved in the Muroidea superfamily of rodents and predates the endogenization of retroviruses presently targeted by ZFP809 in Mus musculus. In sum, these data provide compelling evidence that ZFP809 evolved to recognize foreign DNA and establish histone modification-based epigenetic silencing of ERVs.

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.

C2H2 zinc finger proteins greatly expand the human regulatory lexicon

Cys2-His2 zinc finger (C2H2-ZF) proteins represent the largest class of putative human transcription factors. However, for most C2H2-ZF proteins it is unknown whether they even bind DNA or, if they do, to which sequences. Here, by combining data from a modified bacterial one-hybrid system with protein-binding microarray and chromatin immunoprecipitation analyses, we show that natural C2H2-ZFs encoded in the human genome bind DNA both in vitro and in vivo, and we infer the DNA recognition code using DNA-binding data for thousands of natural C2H2-ZF domains. In vivo binding data are generally consistent with our recognition code and indicate that C2H2-ZF proteins recognize more motifs than all other human transcription factors combined. We provide direct evidence that most KRAB-containing C2H2-ZF proteins bind specific endogenous retroelements (EREs), ranging from currently active to ancient families. The majority of C2H2-ZF proteins, including KRAB proteins, also show widespread binding to regulatory regions, indicating that the human genome contains an extensive and largely unstudied adaptive C2H2-ZF regulatory network that targets a diverse range of genes and pathways.

The Krüppel-associated box repressor domain induces reversible and irreversible regulation of endogenous mouse genes by mediating different chromatin states

The Krüppel-associated box (KRAB) domain is a transcription repression module from the largest family of transcriptional regulators encoded by higher vertebrates. We developed a drug-controllable regulation system based on an artificial KRAB-containing repressor (tTS) that targets the endogenous Hprt gene to explore the regulatory mechanism and molecular basis of KRAB-containing regulators within the context of an endogenous gene in vivo. We show that KRAB can mediate irreversible and reversible regulation of endogenous genes in mouse that is dependent on embryonic developmental stage. KRAB-induced stable DNA methylation within the KRAB binding region during the early embryonic stage, resulting in irreversible gene repression. In later stages, KRAB mainly induced de-acetylation and methylation of histone, resulting in reversible gene repression. Thus, we have characterized the KRAB-mediated regulation system within the context of an endogenous gene and multiple spatiotemporal ranges, thereby providing a basis for identifying the function of KRAB-containing regulators and aiding development of novel KRAB-based gene regulation tools in vivo.

A systematic survey of the Cys2His2 zinc finger DNA-binding landscape

Cys2His2 zinc fingers (C2H2-ZFs) comprise the largest class of metazoan DNA-binding domains. Despite this domain's well-defined DNA-recognition interface, and its successful use in the design of chimeric proteins capable of targeting genomic regions of interest, much remains unknown about its DNA-binding landscape. To help bridge this gap in fundamental knowledge and to provide a resource for design-oriented applications, we screened large synthetic protein libraries to select binding C2H2-ZF domains for each possible three base pair target. The resulting data consist of >160 000 unique domain-DNA interactions and comprise the most comprehensive investigation of C2H2-ZF DNA-binding interactions to date. An integrated analysis of these independent screens yielded DNA-binding profiles for tens of thousands of domains and led to the successful design and prediction of C2H2-ZF DNA-binding specificities. Computational analyses uncovered important aspects of C2H2-ZF domain-DNA interactions, including the roles of within-finger context and domain position on base recognition. We observed the existence of numerous distinct binding strategies for each possible three base pair target and an apparent balance between affinity and specificity of binding. In sum, our comprehensive data help elucidate the complex binding landscape of C2H2-ZF domains and provide a foundation for efforts to determine, predict and engineer their DNA-binding specificities.

The TCF C-clamp DNA binding domain expands the Wnt transcriptome via alternative target recognition

LEF/TCFs direct the final step in Wnt/β-catenin signalling by recruiting β-catenin to genes for activation of transcription. Ancient, non-vertebrate TCFs contain two DNA binding domains, a High Mobility Group box for recognition of the Wnt Response Element (WRE; 5'-CTTTGWWS-3') and the C-clamp domain for recognition of the GC-rich Helper motif (5'-RCCGCC-3'). Two vertebrate TCFs (TCF-1/TCF7 and TCF-4/TCF7L2) use the C-clamp as an alternatively spliced domain to regulate cell-cycle progression, but how the C-clamp influences TCF binding and activity genome-wide is not known. Here, we used a doxycycline inducible system with ChIP-seq to assess how the C-clamp influences human TCF1 binding genome-wide. Metabolic pulse-labeling of nascent RNA with 4'Thiouridine was used with RNA-seq to connect binding to the Wnt transcriptome. We find that the C-clamp enables targeting to a greater number of gene loci for stronger occupancy and transcription regulation. The C-clamp uses Helper sites concurrently with WREs for gene targeting, but it also targets TCF1 to sites that do not have readily identifiable canonical WREs. The coupled ChIP-seq/4'Thiouridine-seq analysis identified new Wnt target genes, including additional regulators of cell proliferation. Thus, C-clamp containing isoforms of TCFs are potent transcriptional regulators with an expanded transcriptome directed by C-clamp-Helper site interactions.

KAP1 promotes proliferation and metastatic progression of breast cancer cells

KAP1 (TRIM28) is a transcriptional regulator in embryonic development that controls stem cell self-renewal, chromatin organization, and the DNA damage response, acting as an essential corepressor for KRAB family zinc finger proteins (KRAB-ZNF). To gain insight into the function of this large gene family, we developed an antibody that recognizes the conserved zinc fingers linker region (ZnFL) in multiple KRAB-ZNF. Here, we report that the expression of many KRAB-ZNF along with active SUMOlyated KAP1 is elevated widely in human breast cancers. KAP1 silencing in breast cancer cells reduced proliferation and inhibited the growth and metastasis of tumor xenografts. Conversely, KAP1 overexpression stimulated cell proliferation and tumor growth. In cells where KAP1 was silenced, we identified multiple downregulated genes linked to tumor progression and metastasis, including EREG/epiregulin, PTGS2/COX2, MMP1, MMP2, and CD44, along with downregulation of multiple KRAB-ZNF proteins. KAP1-dependent stabilization of KRAB-ZNF required direct interactions with KAP1. Together, our results show that KAP1-mediated stimulation of multiple KRAB-ZNF contributes to the growth and metastasis of breast cancer.