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

DNA-binding affinity and specificity determine the phenotypic diversity in BCL11B-related disorders

BCL11B is a Cys2-His2 zinc-finger (C2H2-ZnF) domain-containing, DNA-binding, transcription factor with established roles in the development of various organs and tissues, primarily the immune and nervous systems. BCL11B germline variants have been associated with a variety of developmental syndromes. However, genotype-phenotype correlations along with pathophysiologic mechanisms of selected variants mostly remain elusive. To dissect these, we performed genotype-phenotype correlations of 92 affected individuals harboring a pathogenic or likely pathogenic BCL11B variant, followed by immune phenotyping, analysis of chromatin immunoprecipitation DNA-sequencing data, dual-luciferase reporter assays, and molecular modeling. These integrative analyses enabled us to define three clinical subtypes of BCL11B-related disorders. It is likely that gene-disruptive BCL11B variants and missense variants affecting zinc-binding cysteine and histidine residues cause mild to moderate neurodevelopmental delay with increased propensity for behavioral and dental anomalies, allergies and asthma, and reduced type 2 innate lymphoid cells. Missense variants within C2H2-ZnF DNA-contacting α helices cause highly variable clinical presentations ranging from multisystem anomalies with demise in the first years of life to late-onset, hyperkinetic movement disorder with poor fine motor skills. Those not in direct DNA contact cause a milder phenotype through reduced, target-specific transcriptional activity. However, missense variants affecting C2H2-ZnFs, DNA binding, and "specificity residues" impair BCL11B transcriptional activity in a target-specific, dominant-negative manner along with aberrant regulation of alternative DNA targets, resulting in more severe and unpredictable clinical outcomes. Taken together, we suggest that the phenotypic severity and variability is largely dependent on the DNA-binding affinity and specificity of altered BCL11B proteins.

A transcription network underlies the dual genomic coordination of mitochondrial biogenesis

Mitochondrial biogenesis requires the expression of genes encoded by both the nuclear and mitochondrial genomes. However, aside from a handful transcription factors regulating specific subsets of mitochondrial genes, the overall architecture of the transcriptional control of mitochondrial biogenesis remains to be elucidated. The mechanisms coordinating these two genomes are largely unknown. We performed a targeted RNAi screen in developing eyes with reduced mitochondrial DNA content, anticipating a synergistic disruption of tissue development due to impaired mitochondrial biogenesis and mitochondrial DNA (mtDNA) deficiency. Among 638 transcription factors annotated in the Drosophila genome, 77 were identified as potential regulators of mitochondrial biogenesis. Utilizing published ChIP-seq data of positive hits, we constructed a regulatory network revealing the logic of the transcription regulation of mitochondrial biogenesis. Multiple transcription factors in core layers had extensive connections, collectively governing the expression of nearly all mitochondrial genes, whereas factors sitting on the top layer may respond to cellular cues to modulate mitochondrial biogenesis through the underlying network. CG1603, a core component of the network, was found to be indispensable for the expression of most nuclear mitochondrial genes, including those required for mtDNA maintenance and gene expression, thus coordinating nuclear genome and mtDNA activities in mitochondrial biogenesis. Additional genetic analyses validated YL-1, a transcription factor upstream of CG1603 in the network, as a regulator controlling CG1603 expression and mitochondrial biogenesis.

Activation of recombinases at specific DNA loci by zinc-finger domain insertions

Recombinases have several potential advantages as genome editing tools compared to nucleases and other editing enzymes, but the process of engineering them to efficiently recombine predetermined DNA targets demands considerable investment of time and labor. Here we sought to harness zinc-finger DNA-binding domains (ZFDs) to program recombinase binding by developing fusions, in which ZFDs are inserted into recombinase coding sequences. By screening libraries of hybrid proteins, we optimized the insertion site, linker length, spacing and ZFD orientation and generated Cre-type recombinases that remain dormant unless the insertionally fused ZFD binds its target site placed in the vicinity of the recombinase binding site. The developed fusion improved targeted editing efficiencies of recombinases by four-fold and abolished measurable off-target activity in mammalian cells. The ZFD-dependent activity is transferable to a recombinase with relaxed specificity, providing the means for developing fully programmable recombinases. Our engineered recombinases provide improved genome editing tools with increased precision and efficiency.

A transcription network underlies the dual genomic coordination of mitochondrial biogenesis

Mitochondrial biogenesis requires the expression of genes encoded by both the nuclear and mitochondrial genomes. However, aside from a handful transcriptional factors regulating specific subsets of mitochondrial genes, the overall architecture of the transcriptional control of mitochondrial biogenesis remains to be elucidated. The mechanisms coordinating these two genomes are largely unknown. We performed a targeted RNAi screen in developing eyes with reduced mitochondrial DNA content, anticipating a synergistic disruption of tissue development due to impaired mitochondrial biogenesis and mtDNA deficiency. Among 638 transcription factors annotated in Drosophila genome, 77 were identified as potential regulators of mitochondrial biogenesis. Utilizing published ChIP-seq data of positive hits, we constructed a regulatory network revealing the logic of the transcription regulation of mitochondrial biogenesis. Multiple transcription factors in core layers had extensive connections, collectively governing the expression of nearly all mitochondrial genes, whereas factors sitting on the top layer may respond to cellular cues to modulate mitochondrial biogenesis through the underlying network. CG1603, a core component of the network, was found to be indispensable for the expression of most nuclear mitochondrial genes, including those required for mtDNA maintenance and gene expression, thus coordinating nuclear genome and mtDNA activities in mitochondrial biogenies. Additional genetics analyses validated YL-1, a transcription factor upstream of CG1603 in the network, as a regulator controlling CG1603 expression and mitochondrial biogenesis.

Improving prediction performance of general protein language model by domain-adaptive pretraining on DNA-binding protein

DNA-protein interactions exert the fundamental structure of many pivotal biological processes, such as DNA replication, transcription, and gene regulation. However, accurate and efficient computational methods for identifying these interactions are still lacking. In this study, we propose a method ESM-DBP through refining the DNA-binding protein sequence repertory and domain-adaptive pretraining based the general protein language model. Our method considers the lacking exploration of general language model for DNA-binding protein domain-specific knowledge, so we screen out 170,264 DNA-binding protein sequences to construct the domain-adaptive language model. Experimental results on four downstream tasks show that ESM-DBP provides a better feature representation of DNA-binding protein compared to the original language model, resulting in improved prediction performance and outperforming the state-of-the-art methods. Moreover, ESM-DBP can still perform well even for those sequences with only a few homologous sequences. ChIP-seq on two predicted cases further support the validity of the proposed method.

Modulating Mitochondrial DNA Heteroplasmy with Mitochondrially Targeted Endonucleases

Mitochondria, mainly known as energy factories of eukaryotic cells, also exert several additional signaling and metabolic functions and are today recognized as major cellular biosynthetic and signaling hubs. Mitochondria possess their own genome (mitochondrial DNA-mtDNA), that encodes proteins essential for oxidative phosphorylation, and mutations in it are an important contributor to human disease. The mtDNA mutations often exist in heteroplasmic conditions, with both healthy and mutant versions of the mtDNA residing in patients' cells and the level of mutant mtDNA may vary between different tissues and organs and affect the clinical outcome of the disease. Thus, shifting the ratio between healthy and mutant mtDNA in patients' cells provides an intriguing therapeutic option for mtDNA diseases. In this review we describe current strategies for modulating mitochondrial heteroplasmy levels with engineered endonucleases including mitochondrially targeted TALENs and Zinc finger nucleases (ZFNs) and discuss their therapeutic potential. These gene therapy tools could in the future provide therapeutic help both for patients with mitochondrial disease as well as in preventing the transfer of pathogenic mtDNA mutations from a mother to her offspring.

Geometric deep learning of protein-DNA binding specificity

Predicting protein-DNA binding specificity is a challenging yet essential task for understanding gene regulation. Protein-DNA complexes usually exhibit binding to a selected DNA target site, whereas a protein binds, with varying degrees of binding specificity, to a wide range of DNA sequences. This information is not directly accessible in a single structure. Here, to access this information, we present Deep Predictor of Binding Specificity (DeepPBS), a geometric deep-learning model designed to predict binding specificity from protein-DNA structure. DeepPBS can be applied to experimental or predicted structures. Interpretable protein heavy atom importance scores for interface residues can be extracted. When aggregated at the protein residue level, these scores are validated through mutagenesis experiments. Applied to designed proteins targeting specific DNA sequences, DeepPBS was demonstrated to predict experimentally measured binding specificity. DeepPBS offers a foundation for machine-aided studies that advance our understanding of molecular interactions and guide experimental designs and synthetic biology.

Updated understanding of the protein-DNA recognition code used by C2H2 zinc finger proteins

C2H2 zinc-finger (ZF) proteins form the largest family of DNA-binding transcription factors coded by mammalian genomes. In a typical DNA-binding ZF module, there are twelve residues (numbered from -1 to -12) between the last zinc-coordinating cysteine and the first zinc-coordinating histidine. The established C2H2-ZF "recognition code" suggests that residues at positions -1, -4, and -7 recognize the 5', central, and 3' bases of a DNA base-pair triplet, respectively. Structural studies have highlighted that additional residues at positions -5 and -8 also play roles in specific DNA recognition. The presence of bulky and either charged or polar residues at these five positions determines specificity for given DNA bases: guanine is recognized by arginine, lysine, or histidine; adenine by asparagine or glutamine; thymine or 5-methylcytosine by glutamate; and unmodified cytosine by aspartate. This review discusses recent structural characterizations of C2H2-ZFs that add to our understanding of the principles underlying the C2H2-ZF recognition code.

C2H2 proteins: Evolutionary aspects of domain architecture and diversification

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

New Drosophila promoter-associated architectural protein Mzfp1 interacts with CP190 and is required for housekeeping gene expression and insulator activity

In Drosophila, a group of zinc finger architectural proteins recruits the CP190 protein to the chromatin, an interaction that is essential for the functional activity of promoters and insulators. In this study, we describe a new architectural C2H2 protein called Madf and Zinc-Finger Protein 1 (Mzfp1) that interacts with CP190. Mzfp1 has an unusual structure that includes six C2H2 domains organized in a C-terminal cluster and two tandem MADF domains. Mzfp1 predominantly binds to housekeeping gene promoters located in both euchromatin and heterochromatin genome regions. In vivo mutagenesis studies showed that Mzfp1 is an essential protein, and both MADF domains and the CP190 interaction region are required for its functional activity. The C2H2 cluster is sufficient for the specific binding of Mzfp1 to regulatory elements, while the second MADF domain is required for Mzfp1 recruitment to heterochromatin. Mzfp1 binds to the proximal part of the Fub boundary that separates regulatory domains of the Ubx and abd-A genes in the Bithorax complex. Mzfp1 participates in Fub functions in cooperation with the architectural proteins Pita and Su(Hw). Thus, Mzfp1 is a new architectural C2H2 protein involved in the organization of active promoters and insulators in Drosophila.

Relish-mediated C2H2 zinc finger protein IMZF restores Drosophila immune homeostasis via inhibiting the transcription of Imd/Tak1

The dysregulation of intensity and duration in innate immunity can result in detrimental effects on the body, emphasizing the crucial need for precise regulation. However, the intricate and accurate nature of innate immunity implies the existence of numerous undiscovered innate immunomodulators, particularly transcription factors. In this study, we have identified a Drosophila C2H2 zinc finger protein CG18262, named Immune-mediated Zinc Finger protein (IMZF), capable of suppressing immune responses of Imd pathway. Mechanistically, IMZF serves as a transcription factor that represses the expression of Imd and Tak1. Intriguingly, our findings also reveal that Relish, an NF-κB transcription factor, positively regulates the expression of IMZF, consequently inhibiting the activation of Imd and Tak1 to prevent an exaggerated immune response. Additionally, we have elucidated the pivotal role played by the Relish-IMZF-Imd/Tak1 axis in restoring immune homeostasis of Drosophila Imd pathway. In summary, our findings not only unveil a novel C2H2 zinc finger immunoregulatory transcription factor, IMZF, along with its specific mechanism of immune regulation, but also shed light on the dual functionality of Relish in different stages of the immune response by modulating distinct effectors. This discovery provides new insights and enlightenment into the complex regulation of Drosophila innate immunity.

SETDB1 deletion causes DNA demethylation and upregulation of multiple zinc-finger genes

BACKGROUND

SETDB1 (SET domain bifurcated-1) is a histone H3-lysine 9 (H3K9)-specific methyltransferase that mediates heterochromatin formation and repression of target genes. Despite the assumed functional link between DNA methylation and SETDB1-mediated H3K9 trimethylations, several studies have shown that SETDB1 operates autonomously of DNA methylation in a region- and cell-specific manner. This study analyzes SETDB1-null HAP1 cells through a linked methylome and transcriptome analysis, intending to explore genes controlled by SETDB1-involved DNA methylation.

METHODS AND RESULTS

We investigated SETDB1-mediated regulation of DNA methylation and gene transcription in human HAP1 cells using reduced-representation bisulfite sequencing (RRBS) and RNA sequencing. While two-thirds of differentially methylated CpGs (DMCs) in genic regions were hypomethylated in SETDB1-null cells, we detected a plethora of C2H2-type zinc-finger protein genes (C2H2-ZFP, 223 of 749) among the DMC-associated genes. Most C2H2-ZFPs with DMCs in their promoters were found hypomethylated in SETDB1-KO cells, while other non-ZFP genes with promoter DMCs were not. These C2H2-ZFPs with DMCs in their promoters were significantly upregulated in SETDB1-KO cells. Similarly, C2H2-ZFP genes were upregulated in SETDB1-null 293T cells, suggesting that SETDB1's function in ZFP gene repression is widespread. There are several C2H2-ZFP gene clusters on chromosome 19, which were selectively hypomethylated in SETDB1-KO cells.

CONCLUSIONS

SETDB1 collectively and specifically represses a substantial fraction of the C2H2-ZFP gene family. Through the en-bloc silencing of a set of ZFP genes, SETDB1 may help establish a panel of ZFP proteins that are expressed cell-type specifically and thereby can serve as signature proteins for cellular identity.

Keep Fingers on the CpG Islands

The post-genomic era has ushered in the extensive application of epigenetic editing tools, allowing for precise alterations of gene expression. The use of reprogrammable editors that carry transcriptional corepressors has significant potential for long-term epigenetic silencing for the treatment of human diseases. The ideal scenario involves precise targeting of a specific genomic location by a DNA-binding domain, ensuring there are no off-target effects and that the process yields no genetic remnants aside from specific epigenetic modifications (i.e., DNA methylation). A notable example is a recent study on the mouse Pcsk9 gene, crucial for cholesterol regulation and expressed in hepatocytes, which identified synthetic zinc-finger (ZF) proteins as the most effective DNA-binding editors for silencing Pcsk9 efficiently, specifically, and persistently. This discussion focuses on enhancing the specificity of ZF-array DNA binding by optimizing interactions between specific amino acids and DNA bases across three promoters containing CpG islands.

DeepPBS: Geometric deep learning for interpretable prediction of protein-DNA binding specificity

Predicting specificity in protein-DNA interactions is a challenging yet essential task for understanding gene regulation. Here, we present Deep Predictor of Binding Specificity (DeepPBS), a geometric deep-learning model designed to predict binding specificity across protein families based on protein-DNA structures. The DeepPBS architecture allows investigation of different family-specific recognition patterns. DeepPBS can be applied to predicted structures, and can aid in the modeling of protein-DNA complexes. DeepPBS is interpretable and can be used to calculate protein heavy atom-level importance scores, demonstrated as a case-study on p53-DNA interface. When aggregated at the protein residue level, these scores conform well with alanine scanning mutagenesis experimental data. The inference time for DeepPBS is sufficiently fast for analyzing simulation trajectories, as demonstrated on a molecular-dynamics simulation of a Drosophila Hox-DNA tertiary complex with its cofactor. DeepPBS and its corresponding data resources offer a foundation for machine-aided protein-DNA interaction studies, guiding experimental choices and complex design, as well as advancing our understanding of molecular interactions.

Interpretable neural architecture search and transfer learning for understanding CRISPR-Cas9 off-target enzymatic reactions

Finely tuned enzymatic pathways control cellular processes, and their dysregulation can lead to disease. Developing predictive and interpretable models for these pathways is challenging because of the complexity of the pathways and of the cellular and genomic contexts. Here we introduce Elektrum, a deep learning framework that addresses these challenges with data-driven and biophysically interpretable models for determining the kinetics of biochemical systems. First, it uses in vitro kinetic assays to rapidly hypothesize an ensemble of high-quality kinetically interpretable neural networks (KINNs) that predict reaction rates. It then employs a transfer learning step, where the KINNs are inserted as intermediary layers into deeper convolutional neural networks, fine-tuning the predictions for reaction-dependent in vivo outcomes. We apply Elektrum to predict CRISPR-Cas9 off-target editing probabilities and demonstrate that Elektrum achieves improved performance, regularizes neural network architectures and maintains physical interpretability.

Interpretable neural architecture search and transfer learning for understanding CRISPR/Cas9 off-target enzymatic reactions

Finely-tuned enzymatic pathways control cellular processes, and their dysregulation can lead to disease. Creating predictive and interpretable models for these pathways is challenging because of the complexity of the pathways and of the cellular and genomic contexts. Here we introduce Elektrum, a deep learning framework which addresses these challenges with data-driven and biophysically interpretable models for determining the kinetics of biochemical systems. First, it uses in vitro kinetic assays to rapidly hypothesize an ensemble of high-quality Kinetically Interpretable Neural Networks (KINNs) that predict reaction rates. It then employs a novel transfer learning step, where the KINNs are inserted as intermediary layers into deeper convolutional neural networks, fine-tuning the predictions for reaction-dependent in vivo outcomes. Elektrum makes effective use of the limited, but clean in vitro data and the complex, yet plentiful in vivo data that captures cellular context. We apply Elektrum to predict CRISPR-Cas9 off-target editing probabilities and demonstrate that Elektrum achieves state-of-the-art performance, regularizes neural network architectures, and maintains physical interpretability.

Structures of CTCF-DNA complexes including all 11 zinc fingers

The CCCTC-binding factor (CTCF) binds tens of thousands of enhancers and promoters on mammalian chromosomes by means of its 11 tandem zinc finger (ZF) DNA-binding domain. In addition to the 12-15-bp CORE sequence, some of the CTCF binding sites contain 5' upstream and/or 3' downstream motifs. Here, we describe two structures for overlapping portions of human CTCF, respectively, including ZF1-ZF7 and ZF3-ZF11 in complex with DNA that incorporates the CORE sequence together with either 3' downstream or 5' upstream motifs. Like conventional tandem ZF array proteins, ZF1-ZF7 follow the right-handed twist of the DNA, with each finger occupying and recognizing one triplet of three base pairs in the DNA major groove. ZF8 plays a unique role, acting as a spacer across the DNA minor groove and positioning ZF9-ZF11 to make cross-strand contacts with DNA. We ascribe the difference between the two subgroups of ZF1-ZF7 and ZF8-ZF11 to residues at the two positions -6 and -5 within each finger, with small residues for ZF1-ZF7 and bulkier and polar/charged residues for ZF8-ZF11. ZF8 is also uniquely rich in basic amino acids, which allows salt bridges to DNA phosphates in the minor groove. Highly specific arginine-guanine and glutamine-adenine interactions, used to recognize G:C or A:T base pairs at conventional base-interacting positions of ZFs, also apply to the cross-strand interactions adopted by ZF9-ZF11. The differences between ZF1-ZF7 and ZF8-ZF11 can be rationalized structurally and may contribute to recognition of high-affinity CTCF binding sites.

A universal deep-learning model for zinc finger design enables transcription factor reprogramming

Cys2His2 zinc finger (ZF) domains engineered to bind specific target sequences in the genome provide an effective strategy for programmable regulation of gene expression, with many potential therapeutic applications. However, the structurally intricate engagement of ZF domains with DNA has made their design challenging. Here we describe the screening of 49 billion protein-DNA interactions and the development of a deep-learning model, ZFDesign, that solves ZF design for any genomic target. ZFDesign is a modern machine learning method that models global and target-specific differences induced by a range of library environments and specifically takes into account compatibility of neighboring fingers using a novel hierarchical transformer architecture. We demonstrate the versatility of designed ZFs as nucleases as well as activators and repressors by seamless reprogramming of human transcription factors. These factors could be used to upregulate an allele of haploinsufficiency, downregulate a gain-of-function mutation or test the consequence of regulation of a single gene as opposed to the many genes that a transcription factor would normally influence.

THE1B may have no role in human pregnancy due to ZNF430-mediated silencing

THE1-family retrovirus invaded the primate genome more than 40 million years ago. Dunn-Fletcher et al. reported one THE1B element upstream of CRH gene alters gestation length by upregulating corticotropin-releasing hormone expression in transgenic mice and concluded it has the same role in human as well. However, no promoter or enhancer mark has been detected around this CRH-proximal element in any human tissue or cell, so probably some anti-viral factor exists in primates to prevents it from wreaking havoc. Here I report two paralogous zinc finger genes, ZNF430 and ZNF100, that emerged during the simian lineage to specifically silence THE1B and THE1A, respectively. Contact residue changes in one finger confers each ZNF the unique ability to preferentially repress one THE1 sub-family over the other. The reported THE1B element contains an intact ZNF430 binding site, thus under the repression of ZNF430 in most tissues including placenta, it is questionable whether or not this retrovirus has any role in human pregnancy. Overall, this analysis highlights the need to study human retroviruses' functions in suitable model system.

Exploration of the regulatory relationship between KRAB-Zfp clusters and their target transposable elements via a gene editing strategy at the cluster specific linker-associated sequences by CRISPR-Cas9

BACKGROUND

Krüppel Associated Box-containing Zinc Finger Proteins (KRAB-ZFPs), representing the largest superfamily of transcription factors in mammals, are predicted to primarily target and repress transposable elements (TEs). It is challenging to dissect the distinct functions of these transcription regulators due to their sequence similarity and diversity, and also the complicated repetitiveness of their targeting TE sequences.

RESULTS

Mouse KRAB-Zfps are mainly organized into clusters genomewide. In this study, we revealed that the intra-cluster members had a close evolutionary relationship, and a similar preference for zinc finger (ZnF) usage. KRAB-Zfps were expressed in a cell type- or tissue type specific manner and they tended to be actively transcribed together with other cluster members. Further sequence analyses pointed out the linker sequences in between ZnFs were conserved, and meanwhile had distinct cluster specificity. Based on these unique characteristics of KRAB-Zfp clusters, sgRNAs were designed to edit cluster-specific linkers to abolish the functions of the targeted cluster(s). Using mouse embryonic stem cells (mESC) as a model, we screened and obtained a series of sgRNAs targeting various highly expressed KRAB-Zfp clusters. The effectiveness of sgRNAs were verified in a reporter assay exclusively developed for multi-target sgRNAs and further confirmed by PCR-based analyses. Using mESC cell lines inducibly expressing Cas9 and these sgRNAs, we found that editing different KRAB-Zfp clusters resulted in the transcriptional changes of distinct categories of TEs.

CONCLUSIONS

Collectively, the intrinsic sequence correlations of intra-cluster KRAB-Zfp members discovered in this study suggest that the conserved cluster specific linkers played crucial roles in diversifying the tandem ZnF array and the related target specificity of KRAB-Zfps during clusters' evolution. On this basis, an effective CRISPR-Cas9 based approach against the linker sequences is developed and verified for rapidly editing KRAB-Zfp clusters to identify the regulatory correlation between the cluster members and their potential TE targets.

Function and Evolution of C1-2i Subclass of C2H2-Type Zinc Finger Transcription Factors in POPLAR

C2H2 zinc finger (C2H2-ZF) transcription factors participate in various aspects of normal plant growth regulation and stress responses. C1-2i C2H2-ZFs are a special subclass of conserved proteins that contain two ZnF-C2H2 domains. Some C1-2i C2H2-ZFs in Arabidopsis (ZAT) are involved in stress resistance and other functions. However, there is limited information on C1-2i C2H2-ZFs in Populus trichocarpa (PtriZATs). To analyze the function and evolution of C1-2i C2H2-ZFs, eleven PtriZATs were identified in P. trichocarpa, which can be classified into two subgroups. The protein structure, conserved ZnF-C2H2 domains and QALGGH motifs, showed high conservation during the evolution of PtriZATs in P. trichocarpa. The spacing between two ZnF-C2H2 domains, chromosomal locations and cis-elements implied the original proteins and function of PtriZATs. Furthermore, the gene expression of different tissues and stress treatment showed the functional differentiation of PtriZATs subgroups and their stress response function. The analysis of C1-2i C2H2-ZFs in different Populus species and plants implied their evolution and differentiation, especially in terms of stress resistance. Cis-elements and expression pattern analysis of interaction proteins implied the function of PtriZATs through binding with stress-related genes, which are involved in gene regulation by via epigenetic modification through histone regulation, DNA methylation, ubiquitination, etc. Our results for the origin and evolution of PtriZATs will contribute to understanding the functional differentiation of C1-2i C2H2-ZFs in P. trichocarpa. The interaction and expression results will lay a foundation for the further functional investigation of their roles and biological processes in Populus.

Learning probabilistic protein-DNA recognition codes from DNA-binding specificities using structural mappings

Knowledge of how proteins interact with DNA is essential for understanding gene regulation. Although DNA-binding specificities for thousands of transcription factors (TFs) have been determined, the specific amino acid-base interactions comprising their structural interfaces are largely unknown. This lack of resolution hampers attempts to leverage these data in order to predict specificities for uncharacterized TFs or TFs mutated in disease. Here we introduce recognition code learning via automated mapping of protein-DNA structural interfaces (rCLAMPS), a probabilistic approach that uses DNA-binding specificities for TFs from the same structural family to simultaneously infer both which nucleotide positions are contacted by particular amino acids within the TF as well as a recognition code that relates each base-contacting amino acid to nucleotide preferences at the DNA positions it contacts. We apply rCLAMPS to homeodomains, the second largest family of TFs in metazoans and show that it learns a highly effective recognition code that can predict de novo DNA-binding specificities for TFs. Furthermore, we show that the inferred amino acid-nucleotide contacts reveal whether and how nucleotide preferences at individual binding site positions are altered by mutations within TFs. Our approach is an important step toward automatically uncovering the determinants of protein-DNA specificity from large compendia of DNA-binding specificities and inferring the altered functionalities of TFs mutated in disease.

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.

Auxin-driven ecophysiological diversification of leaves in domesticated tomato

Heterobaric leaves have bundle sheath extensions (BSEs) that compartmentalize the parenchyma, whereas homobaric leaves do not. The presence of BSEs affects leaf hydraulics and photosynthetic rate. The tomato (Solanum lycopersicum) obscuravenosa (obv) mutant lacks BSEs. Here, we identify the obv gene and the causative mutation, a nonsynonymous amino acid change that disrupts a C2H2 zinc finger motif in a putative transcription factor. This mutation exists as a polymorphism in the natural range of wild tomatoes but has increased in frequency in domesticated tomatoes, suggesting that the latter diversified into heterobaric and homobaric leaf types. The obv mutant displays reduced vein density, leaf hydraulic conductance and photosynthetic assimilation rate. We show that these and other pleiotropic effects on plant development, including changes in leaf insertion angle, leaf margin serration, minor vein density, and fruit shape, are controlled by OBV via changes in auxin signaling. Loss of function of the transcriptional regulator AUXIN RESPONSE FACTOR 4 (ARF4) also results in defective BSE development, revealing an additional component of a genetic module controlling aspects of leaf development important for ecological adaptation and subject to breeding selection.

Single base-pair resolution analysis of DNA binding motif with MoMotif reveals an oncogenic function of CTCF zinc-finger 1 mutation

Defining the impact of missense mutations on the recognition of DNA motifs is highly dependent on bioinformatic tools that define DNA binding elements. However, classical motif analysis tools remain limited in their capacity to identify subtle changes in complex binding motifs between distinct conditions. To overcome this limitation, we developed a new tool, MoMotif, that facilitates a sensitive identification, at the single base-pair resolution, of complex, or subtle, alterations to core binding motifs, discerned from ChIP-seq data. We employed MoMotif to define the previously uncharacterized recognition motif of CTCF zinc-finger 1 (ZF1), and to further define the impact of CTCF ZF1 mutation on its association with chromatin. Mutations of CTCF ZF1 are exclusive to breast cancer and are associated with metastasis and therapeutic resistance, but the underlying mechanisms are unclear. Using MoMotif, we identified an extension of the CTCF core binding motif, necessitating a functional ZF1 to bind appropriately. Using a combination of ChIP-Seq and RNA-Seq, we discover that the inability to bind this extended motif drives an altered transcriptional program associated with the oncogenic phenotypes observed clinically. Our study demonstrates that MoMotif is a powerful new tool for comparative ChIP-seq analysis and characterising DNA-protein contacts.

A complex epigenome-splicing crosstalk governs epithelial-to-mesenchymal transition in metastasis and brain development

Epithelial-to-mesenchymal transition (EMT) renders epithelial cells migratory properties. While epigenetic and splicing changes have been implicated in EMT, the mechanisms governing their crosstalk remain poorly understood. Here we discovered that a C2H2 zinc finger protein, ZNF827, is strongly induced during various contexts of EMT, including in brain development and breast cancer metastasis, and is required for the molecular and phenotypic changes underlying EMT in these processes. Mechanistically, ZNF827 mediated these responses by orchestrating a large-scale remodelling of the splicing landscape by recruiting HDAC1 for epigenetic modulation of distinct genomic loci, thereby slowing RNA polymerase II progression and altering the splicing of genes encoding key EMT regulators in cis. Our findings reveal an unprecedented complexity of crosstalk between epigenetic landscape and splicing programme in governing EMT and identify ZNF827 as a master regulator coupling these processes during EMT in brain development and breast cancer metastasis.

Oncogenic inspiration for programmable activators

The discovery of regulatory domains has been limited to the investigation of transcription factors and homologous protein sequences. In this issue of Cell Genomics, motivated by an oncogenic protein fusion, Tak et al.¹ direct the regulatory potential of a nontraditional effector domain to novel genomic loci with fusions to programmable DNA-binding domains.

Neural network modeling of differential binding between wild-type and mutant CTCF reveals putative binding preferences for zinc fingers 1-2

BACKGROUND

Many transcription factors (TFs), such as multi zinc-finger (ZF) TFs, have multiple DNA binding domains (DBDs), and deciphering the DNA binding motifs of individual DBDs is a major challenge. One example of such a TF is CCCTC-binding factor (CTCF), a TF with eleven ZFs that plays a variety of roles in transcriptional regulation, most notably anchoring DNA loops. Previous studies found that CTCF ZFs 3-7 bind CTCF's core motif and ZFs 9-11 bind a specific upstream motif, but the motifs of ZFs 1-2 have yet to be identified.

RESULTS

We developed a new approach to identifying the binding motifs of individual DBDs of a TF through analyzing chromatin immunoprecipitation sequencing (ChIP-seq) experiments in which a single DBD is mutated: we train a deep convolutional neural network to predict whether wild-type TF binding sites are preserved in the mutant TF dataset and interpret the model. We applied this approach to mouse CTCF ChIP-seq data and identified the known binding preferences of CTCF ZFs 3-11 as well as a putative GAG binding motif for ZF 1. We analyzed other CTCF datasets to provide additional evidence that ZF 1 is associated with binding at the motif we identified, and we found that the presence of the motif for ZF 1 is associated with CTCF ChIP-seq peak strength.

CONCLUSIONS

Our approach can be applied to any TF for which in vivo binding data from both the wild-type and mutated versions of the TF are available, and our findings provide new potential insights binding preferences of CTCF's DBDs.

Structures and biological functions of zinc finger proteins and their roles in hepatocellular carcinoma

Zinc finger proteins are transcription factors with the finger domain, which plays a significant role in gene regulation. As the largest family of transcription factors in the human genome, zinc finger (ZNF) proteins are characterized by their different DNA binding motifs, such as C2H2 and Gag knuckle. Different kinds of zinc finger motifs exhibit a wide variety of biological functions. Zinc finger proteins have been reported in various diseases, especially in several cancers. Hepatocellular carcinoma (HCC) is the third leading cause of cancer-associated death worldwide, especially in China. Most of HCC patients have suffered from hepatitis B virus (HBV) and hepatitis C virus (HCV) injection for a long time. Although the surgical operation of HCC has been extremely developed, the prognosis of HCC is still very poor, and the underlying mechanisms in HCC tumorigenesis are still not completely understood. Here, we summarize multiple functions and recent research of zinc finger proteins in HCC tumorigenesis and progression. We also discuss the significance of zinc finger proteins in HCC diagnosis and prognostic evaluation.

OUP accepted manuscript

In oligodendrocytes of the vertebrate central nervous system a complex network of transcriptional regulators is required to ensure correct and timely myelination of neuronal axons. Here we identify Zfp276, the only mammalian ZAD-domain containing zinc finger protein, as a transcriptional regulator of oligodendrocyte differentiation and central myelination downstream of Sox10. In the central nervous system, Zfp276 is exclusively expressed in mature oligodendrocytes. Oligodendroglial deletion of Zfp276 led to strongly reduced expression of myelin genes in the early postnatal mouse spinal cord. Retroviral overexpression of Zfp276 in cultured oligodendrocyte precursor cells induced precocious expression of maturation markers and myelin genes, further supporting its role in oligodendroglial differentiation. On the molecular level, Zfp276 directly binds to and represses Sox10-dependent gene regulatory regions of immaturity factors and functionally interacts with the transcriptional repressor Zeb2 to enable fast transition of oligodendrocytes to the myelinating stage.

Zinc Fingers 10 and 11 of Miz-1 undergo conformational exchange to achieve specific DNA binding

Miz-1 (ZBTB17) is a poly-zinc finger BTB/POZ transcription factor with 12 consecutive C2H2 zinc fingers (ZFs) that binds transcriptional start sites (TSSs) to regulate the expression of genes involved in cell development and proliferation. As of now, it is not known which of the 12 consecutive ZFs are responsible for the recognition of the 24 base pair consensus sequence found at these TSSs. Evidence suggests ZFs 7-12 plays this role. We provide validation for this and describe the structural and dynamical characterization of unprecedented conformational exchange in the linker between ZFs 10 and 11. This conformational exchange uncouples ZFs 7-10 from 11 and 12 and promotes a scanning-recognition mechanism through which the two segments cooperate to bind two sub-sites at both ends of the consensus. We further show that this can result in the coiling of TSSs as part of Miz-1's mechanism of transcriptional transactivation.

Overexpression of AtDREB1 and BcZAT12 genes confers drought tolerance by reducing oxidative stress in double transgenic tomato (Solanum lycopersicum L.)

Double transgenic tomato developed by AtDREB1A and BcZAT12 genes pyramiding showed significant drought tolerance by reducing oxidative stress with enhanced yield. Although a large number of efforts have been made by different researchers to develop abiotic stress tolerance tomato for improving yield using single gene, however, no reports are available which targets AtDREB1 and BcZAT12 genes together. Hence, in the present study, double transgenic plants were developed using AtDREB1 and BcZAT12 genes to improve yield potential with better drought tolerance. Double transgenic (DZ1-DZ5) tomato lines showed enhanced drought tolerance than their counterpart non-transgenic and single transgenic plants at 0, 07, 14, and 21 days of water deficit, respectively. Double transgenic plants showed increased activity of antioxidant enzymes, like catalase (CAT), superoxide dismutase (SOD), glutathione reductase (GR), ascorbate peroxidase (APX), dehydroascorbate reductase (DHAR), monodehydroascorbate reductase (MDHAR) and guaiacol peroxidase (POD), and accumulation of non-enzymatic antioxidants like ascorbic acid, glutathione as compared to non-transgenic and single transgenic. Additionally, the transcript analysis of antioxidant enzymes revealed the increased level of gene expression in double transgenic tomato lines. Developed double-transgenic tomato plants co-over-expressing both genes exhibited more enzymatic and non-enzymatic anti-oxidative activities as compared to the non-transgenic and single transgenic control, respectively. This is the preliminary report in tomato, which forms the basis for a multigene transgenic approach to cope with drought stress.

Overexpression of AtDREB1 and BcZAT12 genes confers drought tolerance by reducing oxidative stress in double transgenic tomato (Solanum lycopersicum L.)

Double transgenic tomato developed by AtDREB1A and BcZAT12 genes pyramiding showed significant drought tolerance by reducing oxidative stress with enhanced yield. Although a large number of efforts have been made by different researchers to develop abiotic stress tolerance tomato for improving yield using single gene, however, no reports are available which targets AtDREB1 and BcZAT12 genes together. Hence, in the present study, double transgenic plants were developed using AtDREB1 and BcZAT12 genes to improve yield potential with better drought tolerance. Double transgenic (DZ1-DZ5) tomato lines showed enhanced drought tolerance than their counterpart non-transgenic and single transgenic plants at 0, 07, 14, and 21 days of water deficit, respectively. Double transgenic plants showed increased activity of antioxidant enzymes, like catalase (CAT), superoxide dismutase (SOD), glutathione reductase (GR), ascorbate peroxidase (APX), dehydroascorbate reductase (DHAR), monodehydroascorbate reductase (MDHAR) and guaiacol peroxidase (POD), and accumulation of non-enzymatic antioxidants like ascorbic acid, glutathione as compared to non-transgenic and single transgenic. Additionally, the transcript analysis of antioxidant enzymes revealed the increased level of gene expression in double transgenic tomato lines. Developed double-transgenic tomato plants co-over-expressing both genes exhibited more enzymatic and non-enzymatic anti-oxidative activities as compared to the non-transgenic and single transgenic control, respectively. This is the preliminary report in tomato, which forms the basis for a multigene transgenic approach to cope with drought stress.

Points of View on the Tools for Genome/Gene Editing

Theoretically, a DNA sequence-specific recognition protein that can distinguish a DNA sequence equal to or more than 16 bp could be unique to mammalian genomes. Long-sequence-specific nucleases, such as naturally occurring Homing Endonucleases and artificially engineered ZFN, TALEN, and Cas9-sgRNA, have been developed and widely applied in genome editing. In contrast to other counterparts, which recognize DNA target sites by the protein moieties themselves, Cas9 uses a single-guide RNA (sgRNA) as a template for DNA target recognition. Due to the simplicity in designing and synthesizing a sgRNA for a target site, Cas9-sgRNA has become the most current tool for genome editing. Moreover, the RNA-guided DNA recognition activity of Cas9-sgRNA is independent of both of the nuclease activities of it on the complementary strand by the HNH domain and the non-complementary strand by the RuvC domain, and HNH nuclease activity null mutant (H840A) and RuvC nuclease activity null mutant (D10A) were identified. In accompaniment with the sgRNA, Cas9, Cas9(D10A), Cas9(H840A), and Cas9(D10A, H840A) can be used to achieve double strand breakage, complementary strand breakage, non-complementary strand breakage, and no breakage on-target site, respectively. Based on such unique characteristics, many engineered enzyme activities, such as DNA methylation, histone methylation, histone acetylation, cytidine deamination, adenine deamination, and primer-directed mutation, could be introduced within or around the target site. In order to prevent off-targeting by the lasting expression of Cas9 derivatives, a lot of transient expression methods, including the direct delivery of Cas9-sgRNA riboprotein, were developed. The issue of biosafety is indispensable in in vivo applications; Cas9-sgRNA packaged into virus-like particles or extracellular vesicles have been designed and some in vivo therapeutic trials have been reported.

On the prediction of DNA-binding preferences of C2H2-ZF domains using structural models: application on human CTCF

Cis2-His2 zinc finger (C2H2-ZF) proteins are the largest family of transcription factors in human and higher metazoans. To date, the DNA-binding preferences of many members of this family remain unknown. We have developed a computational method to predict their DNA-binding preferences. We have computed theoretical position weight matrices (PWMs) of proteins composed by C2H2-ZF domains, with the only requirement of an input structure. We have predicted more than two-third of a single zinc-finger domain binding site for about 70% variants of Zif268, a classical member of this family. We have successfully matched between 60 and 90% of the binding-site motif of examples of proteins composed by three C2H2-ZF domains in JASPAR, a standard database of PWMs. The tests are used as a proof of the capacity to scan a DNA fragment and find the potential binding sites of transcription-factors formed by C2H2-ZF domains. As an example, we have tested the approach to predict the DNA-binding preferences of the human chromatin binding factor CTCF. We offer a server to model the structure of a zinc-finger protein and predict its PWM.

Application of Genome Editing in Tomato Breeding: Mechanisms, Advances, and Prospects

Plants regularly face the changing climatic conditions that cause biotic and abiotic stress responses. The abiotic stresses are the primary constraints affecting crop yield and nutritional quality in many crop plants. The advances in genome sequencing and high-throughput approaches have enabled the researchers to use genome editing tools for the functional characterization of many genes useful for crop improvement. The present review focuses on the genome editing tools for improving many traits such as disease resistance, abiotic stress tolerance, yield, quality, and nutritional aspects of tomato. Many candidate genes conferring tolerance to abiotic stresses such as heat, cold, drought, and salinity stress have been successfully manipulated by gene modification and editing techniques such as RNA interference, insertional mutagenesis, and clustered regularly interspaced short palindromic repeat (CRISPR/Cas9). In this regard, the genome editing tools such as CRISPR/Cas9, which is a fast and efficient technology that can be exploited to explore the genetic resources for the improvement of tomato and other crop plants in terms of stress tolerance and nutritional quality. The review presents examples of gene editing responsible for conferring both biotic and abiotic stresses in tomato simultaneously. The literature on using this powerful technology to improve fruit quality, yield, and nutritional aspects in tomato is highlighted. Finally, the prospects and challenges of genome editing, public and political acceptance in tomato are discussed.

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

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

High-Throughput Affinity Measurements of Transcription Factor and DNA Mutations Reveal Affinity and Specificity Determinants

Transcription factors (TFs) bind regulatory DNA to control gene expression, and mutations to either TFs or DNA can alter binding affinities to rewire regulatory networks and drive phenotypic variation. While studies have profiled energetic effects of DNA mutations extensively, we lack similar information for TF variants. Here, we present STAMMP (simultaneous transcription factor affinity measurements via microfluidic protein arrays), a high-throughput microfluidic platform enabling quantitative characterization of hundreds of TF variants simultaneously. Measured affinities for ∼210 mutants of a model yeast TF (Pho4) interacting with 9 oligonucleotides (>1,800 K_ds) reveal that many combinations of mutations to poorly conserved TF residues and nucleotides flanking the core binding site alter but preserve physiological binding, providing a mechanism by which combinations of mutations in cis and trans could modulate TF binding to tune occupancies during evolution. Moreover, biochemical double-mutant cycles across the TF-DNA interface reveal molecular mechanisms driving recognition, linking sequence to function. A record of this paper's Transparent Peer Review process is included in the Supplemental Information.

The geometric influence on the Cys2His2 zinc finger domain and functional plasticity

The Cys₂His₂ zinc finger is the most common DNA-binding domain expanding in metazoans since the fungi human split. A proposed catalyst for this expansion is an arms race to silence transposable elements yet it remains poorly understood how this domain is able to evolve the required specificities. Likewise, models of its DNA binding specificity remain error prone due to a lack of understanding of how adjacent fingers influence each other's binding specificity. Here, we use a synthetic approach to exhaustively investigate binding geometry, one of the dominant influences on adjacent finger function. By screening over 28 billion protein–DNA interactions in various geometric contexts we find the plasticity of the most common natural geometry enables more functional amino acid combinations across all targets. Further, residues that define this geometry are enriched in genomes where zinc fingers are prevalent and specificity transitions would be limited in alternative geometries. Finally, these results demonstrate an exhaustive synthetic screen can produce an accurate model of domain function while providing mechanistic insight that may have assisted in the domains expansion.

Balance between asymmetry and abundance in multi-domain DNA-binding proteins may regulate the kinetics of their binding to DNA

DNA sequences are often recognized by multi-domain proteins that may have higher affinity and specificity than single-domain proteins. However, the higher affinity to DNA might be coupled with slower recognition kinetics. In this study, we address this balance between stability and kinetics for multi-domain Cys₂His₂- (C₂H₂-) type zinc-finger (ZF) proteins. These proteins are the most prevalent DNA-binding domain in eukaryotes and C₂H₂ type zinc-finger proteins (C₂H₂-ZFPs) constitute nearly one-half of all known and predicted transcription factors in human. Extensive contact with DNA via tandem ZF domains confers high stability on the sequence-specific complexes. However, this can limit target search efficiency, especially for low abundance ZFPs. Earlier, we found that asymmetrical distribution of electrostatic charge among the three ZF domains of the low abundance transcription factor Egr-1 facilitates its DNA search process. Here, on a diverse set of 273 human C₂H₂-ZFP comprised of 3–15 tandem ZF domains, we find that, in many cases, electrostatic charge and binding specificity are asymmetrically distributed among the ZF domains so that neighbouring domains have different DNA-binding properties. For proteins containing 3–6 ZF domains, we show that the low abundance proteins possess a higher degree of non-specific asymmetry and vice versa. Our findings suggest that where the electrostatics of tandem ZF domains are similar (i.e., symmetrical), the ZFPs are more abundant to optimize their DNA search efficiency. This study reveals new insights into the fundamental determinants of recognition by C₂H₂-ZFPs of their DNA binding sites in the cellular landscape. The importance of electrostatic asymmetry with respect to binding site recognition by C₂H₂-ZFPs suggests the possibility that it may also be important in other ZFP systems and reveals a new design feature for zinc finger engineering.

Optimal recognition of proteins to DNA is governed by various factors among them the thermodynamics, kinetics and specificity of the protein-DNA complex. Multi-domain DNA-binding proteins are expected to have higher affinity and specificity due to the extensive interface they form with DNA. However, larger interface may result with higher friction when these proteins scan the DNA for the target site via the sliding mechanism. A way to overcome this drawback is to have asymmetry in the protein so that the interface with DNA is smaller. Alternatively, higher abundance can also increase the search speed. Here, using computational analysis of large data set of multi-domain zinc finger DNA-binding proteins, we report a trade-off between asymmetry and abundance.

Crystal structures of REF6 and its complex with DNA reveal diverse recognition mechanisms

Relative of Early Flowing 6 (REF6) is a DNA-sequence-specific H3K27me3/2 demethylase that contains four zinc finger (ZnF) domains and targets several thousand genes in Arabidopsis thaliana. The ZnF domains are essential for binding target genes, but the structural basis remains unclear. Here, we determined crystal structures of the ZnF domains and REF6-DNA complex, revealing a unique REF6-family-specific half-cross-braced ZnF (RCZ) domain and two C2H2-type ZnFs. DNA-binding induces a profound conformational change in the hinge region of REF6. Each REF6 recognizes six bases and DNA methylation reduces the binding affinity. Both the acidic region and basic region are important for the self-association of REF6. The REF6 DNA-binding affinity is determined by the sequence-dependent conformations of DNA and also the cooperativity in different target motifs. The conformational plasticity enables REF6 to function as a global transcriptional regulator that directly binds to many diverse genes, revealing the structural basis for the epigenetic modification recognition.

C2H2 Zinc Finger Proteins: Master Regulators of Abiotic Stress Responses in Plants

Abiotic stresses such as drought and salinity are major environmental factors that limit crop yields. Unraveling the molecular mechanisms underlying abiotic stress resistance is crucial for improving crop performance and increasing productivity under adverse environmental conditions. Zinc finger proteins, comprising one of the largest transcription factor families, are known for their finger-like structure and their ability to bind Zn2+. Zinc finger proteins are categorized into nine subfamilies based on their conserved Cys and His motifs, including the Cys2/His2-type (C2H2), C3H, C3HC4, C2HC5, C4HC3, C2HC, C4, C6, and C8 subfamilies. Over the past two decades, much progress has been made in understanding the roles of C2H2 zinc finger proteins in plant growth, development, and stress signal transduction. In this review, we focus on recent progress in elucidating the structures, functions, and classifications of plant C2H2 zinc finger proteins and their roles in abiotic stress responses.

Sharing DNA-binding information across structurally similar proteins enables accurate specificity determination

We are now in an era where protein-DNA interactions have been experimentally assayed for thousands of DNA-binding proteins. In order to infer DNA-binding specificities from these data, numerous sophisticated computational methods have been developed. These approaches typically infer DNA-binding specificities by considering interactions for each protein independently, ignoring related and potentially valuable interaction information across other proteins that bind DNA via the same structural domain. Here we introduce a framework for inferring DNA-binding specificities by considering protein-DNA interactions for entire groups of structurally similar proteins simultaneously. We devise both constrained optimization and label propagation algorithms for this task, each balancing observations at the individual protein level against dataset-wide consistency of interaction preferences. We test our approaches on two large, independent Cys2His2 zinc finger protein-DNA interaction datasets. We demonstrate that jointly inferring specificities within each dataset individually dramatically improves accuracy, leading to increased agreement both between these two datasets and with a fixed external standard. Overall, our results suggest that sharing protein-DNA interaction information across structurally similar proteins is a powerful means to enable accurate inference of DNA-binding specificities.

Emerging Genome Engineering Tools in Crop Research and Breeding

Recent advances in genome engineering are revolutionizing crop research and plant breeding. The ability to make specific modifications to a plant's genetic material creates opportunities for rapid development of elite cultivars with desired traits. The plant genome can be altered in several ways, including targeted introduction of nucleotide changes, deleting DNA segments, introducing exogenous DNA fragments and epigenetic modifications. Targeted changes are mediated by sequence specific nucleases (SSNs), such as zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and CRISPR (clustered regularly interspersed short palindromic repeats)-Cas (CRISPR associated protein) systems. Recent advances in engineering chimeric Cas nucleases fused to base editing enzymes permit for even greater precision in base editing and control over gene expression. In addition to gene editing technologies, improvement in delivery systems of exogenous DNA into plant cells have increased the rate of successful gene editing events. Regeneration of fertile plants containing the desired edits remains challenging; however, manipulation of embryogenesis-related genes such as BABY BOOM (BBM) has been shown to facilitate regeneration through tissue culture, often a major hurdle in recalcitrant cultivars. Epigenome reprogramming for improved crop performance is another possibility for future breeders, with recent studies on MutS HOMOLOG 1 (MSH1) demonstrating epigenetic-dependent hybrid vigor in several crops. While these technologies offer plant breeders new tools in creating high yielding, better adapted crop varieties, constantly evolving government policy regarding the cultivation of plants containing transgenes may impede the widespread adoption of some of these techniques. This chapter summarizes advances in genome editing tools and discusses the future of these techniques for crop improvement.

Molecular characterization and expression analysis reveal the roles of Cys2/His2 zinc-finger transcription factors during flower development of Brassica rapa subsp. chinensis

Conserved motif, gene structure, expression and interaction analysis of C2H2-ZFPs in Brassica rapa, and identified types of genes may play essential roles in flower development, and BrZFP38 was proved to function in flower development by affecting pollen formation. Flower development plays a central role in determining the reproduction of higher plants, and Cys₂/His₂ zinc-finger proteins (C2H2-ZFPs) widely participate in the transcriptional regulation of flower development. C2H2-ZFPs with various structures are the most widespread DNA-binding transcription factors in plants. In this study, conserved protein motif and gene structures were analyzed to investigate systematically the molecular features of Brassica rapa C2H2-ZFP genes. Expression of B. rapa C2H2-ZFPs in multiple tissues showed that more than half of the family members with different types ZFs were expressed in flowers. The specific expression profiles of these C2H2-ZFPs in different B. rapa floral bud stages were further evaluated to identify their potential roles in flower development. Interaction networks were constructed in B. rapa based on the orthology of flower-related C2H2-ZFP genes in Arabidopsis. The putative cis-regulatory elements in the promoter regions of these C2H2-ZFP genes were thoroughly analyzed to elucidate their transcriptional regulation. Results showed that the orthologs of known-function flower-related C2H2-ZFP genes were conserved and differentiated in B. rapa. A C2H2-ZFP was proved to function in B. rapa flower development. Our study provides a systematic investigation of the molecular characteristics and expression profiles of C2H2-ZFPs in B. rapa and promotes further work in function and transcriptional regulation of flower development.

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

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

AKT-mediated phosphorylation enhances protein stability and transcription activity of ZNF322A to promote lung cancer progression

ZNF322A is an oncogenic zinc-finger transcription factor. Our published results show that ZNF322A positively regulates transcription of alpha-adducin (ADD1) and cyclin D1 (CCND1) to promote tumorgenicity of lung cancer. However, the upstream regulatory mechanisms of ZNF322A protein function remain elusive. Here, we demonstrate that AKT could phosphorylate ZNF322A by in vitro kinase assay and cell-based mass spectrometry analysis. Overexpression of AKT promoted ZNF322A protein stability and transcriptional activity, whereas these effects were inhibited by knockdown of AKT or treating with AKT inhibitor. We studied AKT-mediated phosphorylation sites, viz. Thr-150, Ser-224, Thr-234, and Thr-262. ZNF322A phosphorylation at Thr-262 by AKT promoted ZNF322A protein stability thus increased ADD1 promoter activity. Interestingly, phosphorylation at Thr-150, Ser-224, and Thr-234 enhanced transcription activity without affecting protein stability of ZNF322A. Chromatin immunoprecipitation and DNA affinity precipitation assays showed that ZNF322A phosphorylation defective mutants Thr-150A, Ser-224A, and Thr-234A attenuated chromatin binding and DNA binding affinity to ADD1 and CCND1 promoters compared with wild-type ZNF322A. Furthermore, AKT-mediated Thr-150, Ser-224, Thr-234, and Thr-262 phosphorylation promoted lung cancer cell growth and metastasis in vitro and in vivo. Clinically, expression of phosphorylated ZNF322A (p-ZNF) correlated with actively phosphorylated AKT (p-AKT) in tumor specimens from 150 lung cancer patients. Multivariate Cox regression analysis indicated that combined p-AKT and p-ZNF expression profile was an independent factor to predict the clinical outcome in lung cancer patients. Our results reveal a new mechanism of AKT signaling in promoting ZNF322A protein stability and transcriptional activity in lung cancer cell, xenograft, and clinical models.

Role for first zinc finger of WT1 in DNA sequence specificity: Denys-Drash syndrome-associated WT1 mutant in ZF1 enhances affinity for a subset of WT1 binding sites

Wilms tumor protein (WT1) is a Cys2-His2 zinc-finger transcription factor vital for embryonic development of the genitourinary system. The protein contains a C-terminal DNA binding domain with four tandem zinc-fingers (ZF1-4). An alternative splicing of Wt1 can add three additional amino acids-lysine (K), threonine (T) and serine (S)-between ZF3 and ZF4. In the -KTS isoform, ZF2-4 determine the sequence-specificity of DNA binding, whereas the function of ZF1 remains elusive. Three X-ray structures are described here for wild-type -KTS isoform ZF1-4 in complex with its cognate DNA sequence. We observed four unique ZF1 conformations. First, like ZF2-4, ZF1 can be positioned continuously in the DNA major groove forming a 'near-cognate' complex. Second, while ZF2-4 make base-specific interactions with one DNA molecule, ZF1 can interact with a second DNA molecule (or, presumably, two regions of the same DNA molecule). Third, ZF1 can intercalate at the joint of two tail-to-head DNA molecules. If such intercalation occurs on a continuous DNA molecule, it would kink the DNA at the ZF1 binding site. Fourth, two ZF1 units can dimerize. Furthermore, we examined a Denys-Drash syndrome-associated ZF1 mutation (methionine at position 342 is replaced by arginine). This mutation enhances WT1 affinity for a guanine base. X-ray crystallography of the mutant in complex with its preferred sequence revealed the interactions responsible for this affinity change. These results provide insight into the mechanisms of action of WT1, and clarify the fact that ZF1 plays a role in determining sequence specificity of this critical transcription factor.