Function and Evolution of C2H2 Zinc Finger Arrays

The functions and roles of C2H2 zinc finger proteins in hepatocellular carcinoma

C2H2 zinc finger (C2H2-ZF) proteins are the majority group of human transcription factors and they have many different molecular functions through different combinations of zinc finger domains. Hepatocellular carcinoma (HCC) is one of the most prevalent malignant tumors and the main reason for cancer-related deaths worldwide. More and more findings support the abnormal expression of C2H2-ZF protein in the onset and progression of HCC. The C2H2-ZF proteins are involved in various biological functions in HCC, such as EMT, stemness maintenance, metabolic reprogramming, cell proliferation and growth, apoptosis, and genomic integrity. The study of anti-tumor drug resistance also highlights the pivotal roles of C2H2-ZF proteins at the intersection of biological functions (EMT, stemness maintenance, autophagy)and chemoresistance in HCC. The involvement of C2H2-ZF protein found recently in regulating different molecules, signal pathways and pathophysiological activities indicate these proteins as the possible therapeutic targets, and diagnostic or prognostic biomarkers for HCC.

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.

DeepZF: improved DNA-binding prediction of C2H2-zinc-finger proteins by deep transfer learning

MOTIVATION

Cys2His2 zinc-finger (C2H2-ZF) proteins are the largest class of human transcription factors and hence play central roles in gene regulation and cell function. C2H2-ZF proteins are characterized by a DNA-binding domain containing multiple ZFs. A subset of the ZFs bind diverse DNA triplets. Despite their central roles, little is known about which of their ZFs are binding and how the DNA-binding preferences are encoded in the amino acid sequence of each ZF.

RESULTS

We present DeepZF, a deep-learning-based pipeline for predicting binding ZFs and their DNA-binding preferences given only the amino acid sequence of a C2H2-ZF protein. To the best of our knowledge, we compiled the first in vivo dataset of binding and non-binding ZFs for training the first ZF-binding classifier. Our classifier, which is based on a novel protein transformer, achieved an average AUROC of 0.71. Moreover, we took advantage of both in vivo and in vitro datasets to learn the recognition code of ZF-DNA binding through transfer learning. Our newly developed model, which is the first to utilize deep learning for the task, achieved an average Pearson correlation greater than 0.94 over each of the three DNA binding positions. Together, DeepZF outperformed extant methods in the task of C2H2-ZF protein DNA-binding preferences prediction: it achieved an average Pearson correlation of 0.42 in motif similarity compared with an average correlation smaller than 0.1 achieved by extant methods. By applying established interpretability techniques, we show that DeepZF inferred biologically relevant binding principles, such as the effect of amino acid residue positions on ZF DNA-binding potential.

AVAILABILITY AND IMPLEMENTATION

DeepZF code, model, and results are available via github.com/OrensteinLab/DeepZF.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Novel ZNF414 activity characterized by integrative analysis of ChIP-exo, ATAC-seq and RNA-seq data

Transcription factor binding to DNA is a central mechanism regulating gene expression. Thus, thorough characterization of this process is essential for understanding cellular biology in both health and disease. We combined data from three sequencing-based methods to unravel the DNA binding function of the novel ZNF414 protein in cells representing two tumor types. ChIP-exo served to map protein binding sites, ATAC-seq allowed identification of open chromatin, and RNA-seq examined the transcriptome. We show that ZNF414 is a DNA-binding protein that both induces and represses gene expression. This transcriptional response has an impact on cellular processes related to proliferation and other malignancy-associated functions, such as cell migration and DNA repair. Approximately 20% of the differentially expressed genes harbored ZNF414 binding sites in their promoters in accessible chromatin, likely representing direct targets of ZNF414. De novo motif discovery revealed several putative ZNF414 binding sequences, one of which was validated using EMSA. In conclusion, this study illustrates a highly efficient integrative approach for the characterization of the DNA binding and transcriptional activity of transcription factors.

Zinc finger protein Zfp335 controls early T cell development and survival through β-selection-dependent and -independent mechanisms

T-cell development in the thymus undergoes the process of differentiation, selective proliferation, and survival from CD4⁻CD8⁻ double negative (DN) stage to CD4⁺CD8⁺ double positive (DP) stage prior to the formation of CD4⁺ helper and CD8⁺ cytolytic T cells ready for circulation. Each developmental stage is tightly regulated by sequentially operating molecular networks, of which only limited numbers of transcription regulators have been deciphered. Here, we identified Zfp335 transcription factor as a new player in the regulatory network controlling thymocyte development in mice. We demonstrate that Zfp335 intrinsically controls DN to DP transition, as T-cell-specific deficiency in Zfp335 leads to a substantial accumulation of DN3 along with reduction of DP, CD4⁺, and CD8⁺ thymocytes. This developmental blockade at DN stage results from the impaired intracellular TCRβ (iTCRβ) expression as well as increased susceptibility to apoptosis in thymocytes. Transcriptomic and ChIP-seq analyses revealed a direct regulation of transcription factors Bcl6 and Rorc by Zfp335. Importantly, enhanced expression of TCRβ and Bcl6/Rorc restores the developmental defect during DN3 to DN4 transition and improves thymocytes survival, respectively. These findings identify a critical role of Zfp335 in controlling T-cell development by maintaining iTCRβ expression-mediated β-selection and independently activating cell survival signaling.

The Honey Bee Gene Bee Antiviral Protein-1 Is a Taxonomically Restricted Antiviral Immune Gene

Insects have evolved a wide range of strategies to combat invading pathogens, including viruses. Genes that encode proteins involved in immune responses often evolve under positive selection due to their co-evolution with pathogens. Insect antiviral defense includes the RNA interference (RNAi) mechanism, which is triggered by recognition of non-self, virally produced, double-stranded RNAs. Indeed, insect RNAi genes (e.g., dicer and argonaute-2) are under high selective pressure. Honey bees (Apis mellifera) are eusocial insects that respond to viral infections via both sequence specific RNAi and a non-sequence specific dsRNA triggered pathway, which is less well-characterized. A transcriptome-level study of virus-infected and/or dsRNA-treated honey bees revealed increased expression of a novel antiviral gene, GenBank: MF116383, and in vivo experiments confirmed its antiviral function. Due to in silico annotation and sequence similarity, MF116383 was originally annotated as a probable cyclin-dependent serine/threonine-protein kinase. In this study, we confirmed that MF116383 limits virus infection, and carried out further bioinformatic and phylogenetic analyses to better characterize this important gene-which we renamed bee antiviral protein-1 (bap1). Phylogenetic analysis revealed that bap1 is taxonomically restricted to Hymenoptera and Blatella germanica (the German cockroach) and that the majority of bap1 amino acids are evolving under neutral selection. This is in-line with the results from structural prediction tools that indicate Bap1 is a highly disordered protein, which likely has relaxed structural constraints. Assessment of honey bee gene expression using a weighted gene correlation network analysis revealed that bap1 expression was highly correlated with several immune genes-most notably argonaute-2. The coexpression of bap1 and argonaute-2 was confirmed in an independent dataset that accounted for the effect of virus abundance. Together, these data demonstrate that bap1 is a taxonomically restricted, rapidly evolving antiviral immune gene. Future work will determine the role of bap1 in limiting replication of other viruses and examine the signal cascade responsible for regulating the expression of bap1 and other honey bee antiviral defense genes, including coexpressed ago-2, and determine whether the virus limiting function of bap1 acts in parallel or in tandem with RNAi.

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.

Transcription Factor KLF2 and Its Role in the Regulation of Inflammatory Processes

KLF2 is a member of the Krüppel-like transcription factor family of proteins containing highly conserved DNA-binding zinc finger domains. KLF2 participates in the differentiation and regulation of the functional activity of monocytes, T lymphocytes, adipocytes, and vascular endothelial cells. The activity of KLF2 is controlled by several regulatory systems, including the MEKK2,3/MEK5/ERK5/MEF2 MAP kinase cascade, Rho family G-proteins, histone acetyltransferases CBP and p300, and histone deacetylases HDAC4 and HDAC5. Activation of KLF2 in endothelial cells induces eNOS expression and provides vasodilatory effect. Many KLF2-dependent genes participate in the suppression of blood coagulation and aggregation of T cells and macrophages with the vascular endothelium, thereby preventing atherosclerosis progression. KLF2 can have a dual effect on the gene transcription. Thus, it induces expression of multiple genes, but suppresses transcription of NF-κB-dependent genes. Transcription factors KLF2 and NF-κB are reciprocal antagonists. KLF2 inhibits induction of NF-κB-dependent genes, whereas NF-κB downregulates KLF2 expression. KLF2-mediated inhibition of NF-κB signaling leads to the suppression of cell response to the pro-inflammatory cytokines IL-1β and TNFα and results in the attenuation of inflammatory processes.

Evolutionary Analysis of the Zinc Finger and Homeoboxes Family of Proteins Identifies Multiple Conserved Domains and a Common Early Chordate Ancestor

The Zinc Fingers and Homeoboxes (Zhx) proteins, Zhx1, Zhx2, and Zhx3, comprise a small family of proteins containing two amino-terminal C₂–H₂ zinc fingers and four or five carboxy-terminal homeodomains. These multiple homeodomains make Zhx proteins unusual because the majority of homeodomain-containing proteins contain a single homeodomain. Studies in cultured cells and mice suggest that Zhx proteins can function as positive or negative transcriptional regulators. Zhx2 regulates numerous hepatic genes, and all three Zhx proteins have been implicated in different cancers. Because Zhx proteins contain multiple predicted homeodomains, are associated with interesting physiological traits, and seem to be only present in the vertebrate lineage, we investigated the evolutionary history of this small family by comparing Zhx homologs from a wide range of chordates. This analysis indicates that the zinc finger motifs and homeodomains are highly similar among all Zhx proteins and also identifies additional Zhx-specific conserved regions, including a 13 amino acid amino-terminal motif that is nearly identical among all gnathostome Zhx proteins. We found single Zhx proteins in the sea lamprey (Petromyzon marinus) and in the nonvertebrate chordates sea squirt (Ciona intestinalis) and lancelet (Branchiostoma floridae); these Zhx proteins are most similar to gnathostome Zhx3. Based on our analyses, we propose that a duplication of the primordial Zhx gene gave rise to Zhx3 and the precursor to Zhx1 and Zhx2. A subsequent tandem duplication of this precursor generated Zhx1 and Zhx2 found in gnathostomes.

The COMET toolkit for composing customizable genetic programs in mammalian cells

Engineering mammalian cells to carry out sophisticated and customizable genetic programs requires a toolkit of multiple orthogonal and well-characterized transcription factors (TFs). To address this need, we develop the COmposable Mammalian Elements of Transcription (COMET)-an ensemble of TFs and promoters that enable the design and tuning of gene expression to an extent not, to the best of our knowledge, previously possible. COMET currently comprises 44 activating and 12 inhibitory zinc-finger TFs and 83 cognate promoters, combined in a framework that readily accommodates new parts. This system can tune gene expression over three orders of magnitude, provides chemically inducible control of TF activity, and enables single-layer Boolean logic. We also develop a mathematical model that provides mechanistic insights into COMET performance characteristics. Altogether, COMET enables the design and construction of customizable genetic programs in mammalian cells.

The Human Transcription Factors

Transcription factors (TFs) recognize specific DNA sequences to control chromatin and transcription, forming a complex system that guides expression of the genome. Despite keen interest in understanding how TFs control gene expression, it remains challenging to determine how the precise genomic binding sites of TFs are specified and how TF binding ultimately relates to regulation of transcription. This review considers how TFs are identified and functionally characterized, principally through the lens of a catalog of over 1,600 likely human TFs and binding motifs for two-thirds of them. Major classes of human TFs differ markedly in their evolutionary trajectories and expression patterns, underscoring distinct functions. TFs likewise underlie many different aspects of human physiology, disease, and variation, highlighting the importance of continued effort to understand TF-mediated gene regulation.

ZFX acts as a transcriptional activator in multiple types of human tumors by binding downstream from transcription start sites at the majority of CpG island promoters

High expression of the transcription factor ZFX is correlated with proliferation, tumorigenesis, and patient survival in multiple types of human cancers. However, the mechanism by which ZFX influences transcriptional regulation has not been determined. We performed ChIP-seq in four cancer cell lines (representing kidney, colon, prostate, and breast cancers) to identify ZFX binding sites throughout the human genome. We identified roughly 9000 ZFX binding sites and found that most of the sites are in CpG island promoters. Moreover, genes with promoters bound by ZFX are expressed at higher levels than genes with promoters not bound by ZFX. To determine if ZFX contributes to regulation of the promoters to which it is bound, we performed RNA-seq analysis after knockdown of ZFX by siRNA in prostate and breast cancer cells. Many genes with promoters bound by ZFX were down-regulated upon ZFX knockdown, supporting the hypothesis that ZFX acts as a transcriptional activator. Surprisingly, ZFX binds at +240 bp downstream from the TSS of the responsive promoters. Using Nucleosome Occupancy and Methylome Sequencing (NOMe-seq), we show that ZFX binds between the open chromatin region at the TSS and the first downstream nucleosome, suggesting that ZFX may play a critical role in promoter architecture. We have also shown that a closely related zinc finger protein ZNF711 has a similar binding pattern at CpG island promoters, but ZNF711 may play a subordinate role to ZFX. This functional characterization of ZFX provides important new insights into transcription, chromatin structure, and the regulation of the cancer transcriptome.

Individual retrotransposon integrants are differentially controlled by KZFP/KAP1-dependent histone methylation, DNA methylation and TET-mediated hydroxymethylation in naïve embryonic stem cells

Background

The KZFP/KAP1 (KRAB zinc finger proteins/KRAB-associated protein 1) system plays a central role in repressing transposable elements (TEs) and maintaining parent-of-origin DNA methylation at imprinting control regions (ICRs) during the wave of genome-wide reprogramming that precedes implantation. In naïve murine embryonic stem cells (mESCs), the genome is maintained highly hypomethylated by a combination of TET-mediated active demethylation and lack of de novo methylation, yet KAP1 is tethered by sequence-specific KZFPs to ICRs and TEs where it recruits histone and DNA methyltransferases to impose heterochromatin formation and DNA methylation.

Results

Here, upon removing either KAP1 or the cognate KZFP, we observed rapid TET2-dependent accumulation of 5hmC at both ICRs and TEs. In the absence of the KZFP/KAP1 complex, ICRs lost heterochromatic histone marks and underwent both active and passive DNA demethylation. For KAP1-bound TEs, 5mC hydroxylation correlated with transcriptional reactivation. Using RNA-seq, we further compared the expression profiles of TEs upon Kap1 removal in wild-type, Dnmt and Tet triple knockout mESCs. While we found that KAP1 represents the main effector of TEs repression in all three settings, we could additionally identify specific groups of TEs further controlled by DNA methylation. Furthermore, we observed that in the absence of TET proteins, activation upon Kap1 depletion was blunted for some TE integrants and increased for others.

Conclusions

Our results indicate that the KZFP/KAP1 complex maintains heterochromatin and DNA methylation at ICRs and TEs in naïve embryonic stem cells partly by protecting these loci from TET-mediated demethylation. Our study further unveils an unsuspected level of complexity in the transcriptional control of the endovirome by demonstrating often integrant-specific differential influences of histone-based heterochromatin modifications, DNA methylation and 5mC oxidation in regulating TEs expression.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0177-1) contains supplementary material, which is available to authorized users.

Integrated analysis of motif activity and gene expression changes of transcription factors

The ability to predict transcription factors based on sequence information in regulatory elements is a key step in systems-level investigation of transcriptional regulation. Here, we have developed a novel tool, IMAGE, for precise prediction of causal transcription factors based on transcriptome profiling and genome-wide maps of enhancer activity. High precision is obtained by combining a near-complete database of position weight matrices (PWMs), generated by compiling public databases and systematic prediction of PWMs for uncharacterized transcription factors, with a state-of-the-art method for PWM scoring and a novel machine learning strategy, based on both enhancers and promoters, to predict the contribution of motifs to transcriptional activity. We applied IMAGE to published data obtained during 3T3-L1 adipocyte differentiation and showed that IMAGE predicts causal transcriptional regulators of this process with higher confidence than existing methods. Furthermore, we generated genome-wide maps of enhancer activity and transcripts during human mesenchymal stem cell commitment and adipocyte differentiation and used IMAGE to identify positive and negative transcriptional regulators of this process. Collectively, our results demonstrate that IMAGE is a powerful and precise method for prediction of regulators of gene expression.

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

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

A sequence-specific core promoter-binding transcription factor recruits TRF2 to coordinately transcribe ribosomal protein genes

Ribosomal protein (RP) genes must be coordinately expressed for proper assembly of the ribosome yet the mechanisms that control expression of RP genes in metazoans are poorly understood. Recently, TATA-binding protein-related factor 2 (TRF2) rather than the TATA-binding protein (TBP) was found to function in transcription of RP genes in Drosophila. Unlike TBP, TRF2 lacks sequence-specific DNA binding activity, so the mechanism by which TRF2 is recruited to promoters is unclear. We show that the transcription factor M1BP, which associates with the core promoter region, activates transcription of RP genes. Moreover, M1BP directly interacts with TRF2 to recruit it to the RP gene promoter. High resolution ChIP-exo was used to analyze in vivo the association of M1BP, TRF2 and TFIID subunit, TAF1. Despite recent work suggesting that TFIID does not associate with RP genes in Drosophila, we find that TAF1 is present at RP gene promoters and that its interaction might also be directed by M1BP. Although M1BP associates with thousands of genes, its colocalization with TRF2 is largely restricted to RP genes, suggesting that this combination is key to coordinately regulating transcription of the majority of RP genes in Drosophila.

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.

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.

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.

Mechanisms of transcription factor evolution in Metazoa

Transcriptions factors (TFs) are pivotal for the regulation of virtually all cellular processes, including growth and development. Expansions of TF families are causally linked to increases in organismal complexity. Here we study the evolutionary dynamics, genetic causes and functional implications of the five largest metazoan TF families. We find that family expansions dominate across the whole metazoan tree; however, some branches experience exceptional family-specific accelerated expansions. Additionally, we find that such expansions are often predated by modular domain rearrangements, which spur the expansion of a new sub-family by separating it from the rest of the TF family in terms of protein-protein interactions. This separation allows for radical shifts in the functional spectrum of a duplicated TF. We also find functional differentiation inside TF sub-families as changes in expression specificity. Furthermore, accelerated family expansions are facilitated by repeats of sequence motifs such as C2H2 zinc fingers. We quantify whole genome duplications and single gene duplications as sources of TF family expansions, implying that some, but not all, TF duplicates are preferentially retained. We conclude that trans-regulatory changes (domain rearrangements) are instrumental for fundamental functional innovations, that cis-regulatory changes (affecting expression) accomplish wide-spread fine tuning and both jointly contribute to the functional diversification of TFs.

Identification and Characterization of Wor4, a New Transcriptional Regulator of White-Opaque Switching

The human fungal pathogen Candida albicans can switch between two cell types, "white" and "opaque," each of which is heritable through many cell divisions. Switching between these two cell types is regulated by six transcriptional regulators that form a highly interconnected circuit with multiple feedback loops. Here, we identify a seventh regulator of white-opaque switching, which we have named Wor4. We show that ectopic expression of Wor4 is sufficient to drive switching from the white to the opaque cell type, and that deletion of Wor4 blocks switching from the white to the opaque cell type. A combination of ectopic expression and deletion experiments indicates that Wor4 is positioned upstream of Wor1, and that it is formally an activator of the opaque cell type. The combination of ectopic expression and deletion phenotypes for Wor4 is unique; none of the other six white-opaque regulators show this pattern. We determined the pattern of Wor4 binding across the genome by ChIP-seq and found it is highly correlated with that of Wor1 and Wor2, indicating that Wor4 is tightly integrated into the existing white-opaque regulatory circuit. We previously proposed that white-to-opaque switching relies on the activation of a complex circuit of feedback loops that remains excited through many cell divisions. The identification of a new, central regulator of white-opaque switching supports this idea by indicating that the white-opaque switching mechanism is considerably more complex than those controlling conventional, nonheritable patterns of gene expression.

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.

KIT Mutation and Loss of 14q May Be Sufficient for the Development of Clinically Symptomatic Very Low-Risk GIST

The aim of this study was to determine the minimal set of genetic alterations required for the development of a very low risk clinically symptomatic gastro-intestinal stromal tumour within the stomach wall. We studied the genome of a very low-risk gastric gastro-intestinal stromal tumour by whole-genome sequencing, comparative genomic hybridisation and methylation profiling. The studied tumour harboured two typical genomic lesions: loss of the long arm of chromosome 14 and an activating mutation in exon 11 of KIT. Besides these genetic lesions, only two point mutations that may affect tumour progression were identified: A frame-shift deletion in RNF146 and a missense mutation in a zinc finger of ZNF407. Whilst the frameshift deletion in RNF146 seemed to be restricted to this particular tumour, a similar yet germline mutation in ZNF407 was found in a panel of 52 gastro-intestinal stromal tumours from different anatomical sites and different categories. Germline polymorphisms in the mitotic checkpoint proteins Aurora kinase A and BUB1 kinase B may have furthered tumour growth. The epigenetic profile of the tumour matches that of other KIT-mutant tumours. We have identified mutations in three genes and loss of the long arm of chromosome 14 as the so far minimal set of genetic abnormalities sufficient for the development of a very low risk clinically symptomatic gastric stromal tumour.

Mapping and analysis of Caenorhabditis elegans transcription factor sequence specificities

Caenorhabditis elegans is a powerful model for studying gene regulation, as it has a compact genome and a wealth of genomic tools. However, identification of regulatory elements has been limited, as DNA-binding motifs are known for only 71 of the estimated 763 sequence-specific transcription factors (TFs). To address this problem, we performed protein binding microarray experiments on representatives of canonical TF families in C. elegans, obtaining motifs for 129 TFs. Additionally, we predict motifs for many TFs that have DNA-binding domains similar to those already characterized, increasing coverage of binding specificities to 292 C. elegans TFs (∼40%). These data highlight the diversification of binding motifs for the nuclear hormone receptor and C2H2 zinc finger families and reveal unexpected diversity of motifs for T-box and DM families. Motif enrichment in promoters of functionally related genes is consistent with known biology and also identifies putative regulatory roles for unstudied TFs.

DOI:http://dx.doi.org/10.7554/eLife.06967.001

Many scientists use ‘model’ species—such as the fruit fly or a nematode worm called Caenorhabditis elegans—in their research because these organisms have useful features that make it easier to carry out many experiments. For example, C. elegans has a smaller genome compared to many other animals, which is useful for studying the roles of individual genes or stretches of DNA.

Transcription factors are a type of protein that can bind to specific stretches of DNA and help to switch certain genes on or off. These ‘motifs’ may be close to the gene or further away in the genome, and therefore, must stand out amongst the rest of the DNA, like lights on a landing strip. However, the motifs for only 10% of the estimated 763 transcription factors in C. elegans have been identified so far.

In this study, Narasimhan, Lambert, Yang et al. used a technique called a ‘protein binding microarray’ to identify the motifs for many more of the C. elegans transcription factors. These findings were then used to predict motifs for other transcription factors. Together, these methods increased the proportion of C. elegans transcription factors with known DNA-binding motifs from 10% to around 40%.

Now that more DNA motifs have been identified, it is possible to look for similarities and differences between them. For example, Narasimhan, Lambert, Yang et al. found that transcription factors with similar sequences can bind to very varied motifs. On the other hand, some transcription factors that are very different are able to recognize very similar motifs. The experiments also indicate that motifs found very close to genes—in sequences known as ‘promoters’—may be able to interact with many proteins to influence the activity of genes.

Narasimhan, Lambert, Yang et al.'s findings increase the number of C. elegans transcription factors with a motif, bringing the knowledge of these proteins more in line with the better-studied transcription factors of humans and fruit flies. The next challenge is to identify DNA motifs for the remaining 60% of transcription factors.

DOI:http://dx.doi.org/10.7554/eLife.06967.002

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.

Determination and inference of eukaryotic transcription factor sequence specificity

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

Deep vertebrate roots for mammalian zinc finger transcription factor subfamilies

While many vertebrate transcription factor (TF) families are conserved, the C2H2 zinc finger (ZNF) family stands out as a notable exception. In particular, novel ZNF gene types have arisen, duplicated, and diverged independently throughout evolution to yield many lineage-specific TF genes. This evolutionary dynamic not only raises many intriguing questions but also severely complicates identification of those ZNF genes that remain functionally conserved. To address this problem, we searched for vertebrate “DNA binding orthologs” by mining ZNF loci from eight sequenced genomes and then aligning the patterns of DNA-binding amino acids, or “fingerprints,” extracted from the encoded ZNF motifs. Using this approach, we found hundreds of lineage-specific genes in each species and also hundreds of orthologous groups. Most groups of orthologs displayed some degree of fingerprint divergence between species, but 174 groups showed fingerprint patterns that have been very rigidly conserved. Focusing on the dynamic KRAB-ZNF subfamily—including nearly 400 human genes thought to possess potent KRAB-mediated epigenetic silencing activities—we found only three genes conserved between mammals and nonmammalian groups. These three genes, members of an ancient familial cluster, encode an unusual KRAB domain that functions as a transcriptional activator. Evolutionary analysis confirms the ancient provenance of this activating KRAB and reveals the independent expansion of KRAB-ZNFs in every vertebrate lineage. Most human ZNF genes, from the most deeply conserved to the primate-specific genes, are highly expressed in immune and reproductive tissues, indicating that they have been enlisted to regulate evolutionarily divergent biological traits.

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.