De novo prediction of DNA-binding specificities for Cys2His2 zinc finger proteins

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.

Structure-based learning to predict and model protein-DNA interactions and transcription-factor co-operativity in cis-regulatory elements

Transcription factor (TF) binding is a key component of genomic regulation. There are numerous high-throughput experimental methods to characterize TF-DNA binding specificities. Their application, however, is both laborious and expensive, which makes profiling all TFs challenging. For instance, the binding preferences of ∼25% human TFs remain unknown; they neither have been determined experimentally nor inferred computationally. We introduce a structure-based learning approach to predict the binding preferences of TFs and the automated modelling of TF regulatory complexes. We show the advantage of using our approach over the classical nearest-neighbor prediction in the limits of remote homology. Starting from a TF sequence or structure, we predict binding preferences in the form of motifs that are then used to scan a DNA sequence for occurrences. The best matches are either profiled with a binding score or collected for their subsequent modeling into a higher-order regulatory complex with DNA. Co-operativity is modelled by: (i) the co-localization of TFs and (ii) the structural modeling of protein-protein interactions between TFs and with co-factors. We have applied our approach to automatically model the interferon-β enhanceosome and the pioneering complexes of OCT4, SOX2 (or SOX11) and KLF4 with a nucleosome, which are compared with the experimentally known structures.

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.

Protein G-quadruplex interactions and their effects on phase transitions and protein aggregation

The SERF family of proteins were originally discovered for their ability to accelerate amyloid formation. Znf706 is an uncharacterized protein whose N-terminus is homologous to SERF proteins. We show here that human Znf706 can promote protein aggregation and amyloid formation. Unexpectedly, Znf706 specifically interacts with stable, non-canonical nucleic acid structures known as G-quadruplexes. G-quadruplexes can affect gene regulation and suppress protein aggregation; however, it is unknown if and how these two activities are linked. We find Znf706 binds preferentially to parallel G-quadruplexes with low micromolar affinity, primarily using its N-terminus, and upon interaction, its dynamics are constrained. G-quadruplex binding suppresses Znf706's ability to promote protein aggregation. Znf706 in conjunction with G-quadruplexes therefore may play a role in regulating protein folding. RNAseq analysis shows that Znf706 depletion specifically impacts the mRNA abundance of genes that are predicted to contain high G-quadruplex density. Our studies give insight into how proteins and G-quadruplexes interact, and how these interactions affect both partners and lead to the modulation of protein aggregation and cellular mRNA levels. These observations suggest that the SERF family of proteins, in conjunction with G-quadruplexes, may have a broader role in regulating protein folding and gene expression than previously appreciated.

Exploration of Genome-Wide Recombination Rate Variation Patterns at Different Scales in Pigs

Meiotic recombination is a prevalent process in eukaryotic sexual reproduction organisms that plays key roles in genetic diversity, breed selection, and species evolution. However, the recombination events differ across breeds and even within breeds. In this study, we initially computed large-scale population recombination rates of both sexes using approximately 52 K SNP genotypes in a total of 3279 pigs from four different Chinese and Western breeds. We then constructed a high-resolution historical recombination map using approximately 16 million SNPs from a sample of unrelated individuals. Comparative analysis of porcine recombination events from different breeds and at different resolutions revealed the following observations: Firstly, the 1Mb-scale pig recombination maps of the same sex are moderately conserved among different breeds, with the similarity of recombination events between Western pigs and Chinese indigenous pigs being lower than within their respective groups. Secondly, we identified 3861 recombination hotspots in the genome and observed medium- to high-level correlation between historical recombination rates (0.542~0.683) and estimates of meiotic recombination rates. Third, we observed that recombination hotspots are significantly far from the transcription start sites of pig genes, and the silico-predicted PRDM9 zinc finger domain DNA recognition motif is significantly enriched in the regions of recombination hotspots compared to recombination coldspots, highlighting the potential role of PRDM9 in regulating recombination hotspots in pigs. Our study analyzed the variation patterns of the pig recombination map at broad and fine scales, providing a valuable reference for genomic selection breeding and laying a crucial foundation for further understanding the molecular mechanisms of pig genome recombination.

The BTB-ZF gene Bm-mamo regulates pigmentation in silkworm caterpillars

The color pattern of insects is one of the most diverse adaptive evolutionary phenotypes. However, the molecular regulation of this color pattern is not fully understood. In this study, we found that the transcription factor Bm-mamo is responsible for black dilute (bd) allele mutations in the silkworm. Bm-mamo belongs to the BTB zinc finger family and is orthologous to mamo in Drosophila melanogaster. This gene has a conserved function in gamete production in Drosophila and silkworms and has evolved a pleiotropic function in the regulation of color patterns in caterpillars. Using RNAi and clustered regularly interspaced short palindromic repeats (CRISPR) technology, we showed that Bm-mamo is a repressor of dark melanin patterns in the larval epidermis. Using in vitro binding assays and gene expression profiling in wild-type and mutant larvae, we also showed that Bm-mamo likely regulates the expression of related pigment synthesis and cuticular protein genes in a coordinated manner to mediate its role in color pattern formation. This mechanism is consistent with the dual role of this transcription factor in regulating both the structure and shape of the cuticle and the pigments that are embedded within it. This study provides new insight into the regulation of color patterns as well as into the construction of more complex epidermal features in some insects.

The N-terminal dimerization domains of human and Drosophila CTCF have similar functionality

BACKGROUND

CTCF is highly likely to be the ancestor of proteins that contain large clusters of C2H2 zinc finger domains, and its conservation is observed across most bilaterian organisms. In mammals, CTCF is the primary architectural protein involved in organizing chromosome topology and mediating enhancer-promoter interactions over long distances. In Drosophila, CTCF (dCTCF) cooperates with other architectural proteins to establish long-range interactions and chromatin boundaries. CTCFs of various organisms contain an unstructured N-terminal dimerization domain (DD) and clusters comprising eleven zinc-finger domains of the C2H2 type. The Drosophila (dCTCF) and human (hCTCF) CTCFs share sequence homology in only five C2H2 domains that specifically bind to a conserved 15 bp motif.

RESULTS

Previously, we demonstrated that CTCFs from different organisms carry unstructured N-terminal dimerization domains (DDs) that lack sequence homology. Here we used the CTCFattP(mCh) platform to introduce desired changes in the Drosophila CTCF gene and generated a series of transgenic lines expressing dCTCF with different variants of the N-terminal domain. Our findings revealed that the functionality of dCTCF is significantly affected by the deletion of the N-terminal DD. Additionally, we observed a strong impact on the binding of the dCTCF mutant to chromatin upon deletion of the DD. However, chromatin binding was restored in transgenic flies expressing a chimeric CTCF protein with the DD of hCTCF. Although the chimeric protein exhibited lower expression levels than those of the dCTCF variants, it efficiently bound to chromatin similarly to the wild type (wt) protein.

CONCLUSIONS

Our findings suggest that one of the evolutionarily conserved functions of the unstructured N-terminal dimerization domain is to recruit dCTCF to its genomic sites in vivo.

Protein G-quadruplex interactions and their effects on phase transitions and protein aggregation

The SERF family of proteins were originally discovered for their ability to accelerate amyloid formation. Znf706 is an uncharacterized protein whose N-terminus is homologous to SERF proteins. We show here that human Znf706 can promote protein aggregation and amyloid formation. Unexpectedly, Znf706 specifically interacts with stable, non-canonical nucleic acid structures known as G-quadruplexes. G-quadruplexes can affect gene regulation and suppress protein aggregation; however, it is unknown if and how these two activities are linked. We find Znf706 binds preferentially to parallel G-quadruplexes with low micromolar affinity, primarily using its N-terminus, and upon interaction, its dynamics are constrained. G-quadruplex binding suppresses Znf706's ability to promote protein aggregation. Znf706 in conjunction with G-quadruplexes therefore may play a role in regulating protein folding. RNAseq analysis shows that Znf706 depletion specifically impacts the mRNA abundance of genes that are predicted to contain high G-quadruplex density. Our studies give insight into how proteins and G-quadruplexes interact, and how these interactions affect both partners and lead to the modulation of protein aggregation and cellular mRNA levels. These observations suggest that the SERF family of proteins, in conjunction with G-quadruplexes, may have a broader role in regulating protein folding and gene expression than previously appreciated.

Natural variation in the zinc-finger-encoding exon of Prdm9 affects hybrid sterility phenotypes in mice

PRDM9-mediated reproductive isolation was first described in the progeny of Mus musculus musculus (MUS) PWD/Ph and Mus musculus domesticus (DOM) C57BL/6J inbred strains. These male F1 hybrids fail to complete chromosome synapsis and arrest meiosis at prophase I, due to incompatibilities between the Prdm9 gene and hybrid sterility locus Hstx2. We identified 14 alleles of Prdm9 in exon 12, encoding the DNA-binding domain of the PRDM9 protein in outcrossed wild mouse populations from Europe, Asia, and the Middle East, 8 of which are novel. The same allele was found in all mice bearing introgressed t-haplotypes encompassing Prdm9. We asked whether 7 novel Prdm9 alleles in MUS populations and the t-haplotype allele in 1 MUS and 3 DOM populations induce Prdm9-mediated reproductive isolation. The results show that only combinations of the dom2 allele of DOM origin and the MUS msc1 allele ensure complete infertility of intersubspecific hybrids in outcrossed wild populations and inbred mouse strains examined so far. The results further indicate that MUS mice may share the erasure of PRDM9msc1 binding motifs in populations with different Prdm9 alleles, which implies that erased PRDM9 binding motifs may be uncoupled from their corresponding Prdm9 alleles at the population level. Our data corroborate the model of Prdm9-mediated hybrid sterility beyond inbred strains of mice and suggest that sterility alleles of Prdm9 may be rare.

Host-encoded CTCF regulates human cytomegalovirus latency via chromatin looping

Human cytomegalovirus (HCMV) is a prevalent pathogen that establishes life-long latent infection in hematopoietic cells. While this infection is usually asymptomatic, immune dysregulation leads to viral reactivation, which can cause significant morbidity and mortality. However, the mechanisms underpinning reactivation remain incompletely understood. The HCMV major immediate early promoter (MIEP)/enhancer is a key factor in this process, as its transactivation from a repressed to active state helps drive viral gene transcription necessary for reactivation from latency. Numerous host transcription factors bind the MIE locus and recruit repressive chromatin modifiers, thus impeding virus reactivation. One such factor is CCCTC-binding protein (CTCF), a highly conserved host zinc finger protein that mediates chromatin conformation and nuclear architecture. However, the mechanisms by which CTCF contributes to HCMV latency were previously unexplored. Here, we confirm that CTCF binds two convergent sites within the MIE locus during latency in primary CD14+ monocytes, and following cellular differentiation, CTCF association is lost as the virus reactivates. While mutation of the MIE enhancer CTCF binding site does not impact viral lytic growth in fibroblasts, this mutant virus fails to maintain latency in myeloid cells. Furthermore, we show the two convergent CTCF binding sites allow looping to occur across the MIEP, supporting transcriptional repression during latency. Indeed, looping between the two sites diminishes during virus reactivation, concurrent with activation of MIE transcription. Taken together, our data reveal that three-dimensional chromatin looping aids in the regulation of HCMV latency and provides insight into promoter/enhancer regulation that may prove broadly applicable across biological systems.

Interaction of MLE with CLAMP zinc finger is involved in proper MSL proteins binding to chromosomes in Drosophila

The Drosophila male-specific lethal (MSL) complex binds to the male X chromosome to activate transcription. It comprises five proteins (MSL1, MSL2, MSL3, male absent on the first (MOF), and maleless (MLE)) and two long noncoding RNAs (lncRNAs; roX1 and roX2). The MLE helicase remodels the roX lncRNAs, enabling the lncRNA-mediated assembly of the Drosophila dosage compensation complex. MSL2 is expressed only in males and interacts with the N-terminal zinc finger of the transcription factor chromatin-linked adapter for MSL proteins (CLAMP), which is important for the specific recruitment of the MSL complex to the male X chromosome. Here, we found that MLE's unstructured C-terminal region interacts with the sixth zinc-finger domain of CLAMP. In vitro, 4-5 zinc fingers are critical for the specific DNA-binding of CLAMP with GA repeats, which constitute the core motif at the high affinity binding sites for MSL proteins. Deleting the CLAMP binding region in MLE decreases the association of MSL proteins with the male X chromosome and increases male lethality. These results suggest that interactions of unstructured regions in MSL2 and MLE with CLAMP zinc finger domains are important for the specific recruitment of the MSL complex to the male X chromosome.

Patterns of recombination in snakes reveal a tug-of-war between PRDM9 and promoter-like features

In some mammals, notably humans, recombination occurs almost exclusively where the protein PRDM9 binds, whereas in vertebrates lacking an intact PRDM9, such as birds and canids, recombination rates are elevated near promoter-like features. To determine whether PRDM9 directs recombination in nonmammalian vertebrates, we focused on an exemplar species with a single, intact PRDM9 ortholog, the corn snake (Pantherophis guttatus). Analyzing historical recombination rates along the genome and crossovers in pedigrees, we found evidence that PRDM9 specifies the location of recombination events, but we also detected a separable effect of promoter-like features. These findings reveal that the uses of PRDM9 and promoter-like features need not be mutually exclusive and instead reflect a tug-of-war that is more even in some species than others.

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.

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.

Overexpression of the Transcription Factor Azf1 Reveals Novel Regulatory Functions and Impacts β-Glucosidase Production in Trichoderma reesei

The fungus Trichoderma reesei is an essential producer of enzymes that degrade lignocellulosic biomass to produce value-added bioproducts. The cellulolytic system of T. reesei is controlled by several transcription factors (TFs) that efficiently regulate the production of these enzymes. Recently, a new TF named Azf1 was identified as a positive regulator of cellulase expression. Here, we investigated novel regulatory functions of Azf1 by its overexpression. In the mutant strain OEazf1, overexpression of azf1 was achieved under both repression and induction conditions. Although azf1 was more abundant in transcript and protein, overexpression of this TF did not activate transcription of the cellulase gene in the presence of the repressor glucose, suggesting that Azf1 may be subject to posttranslational regulation. In cellulose, the expression of swo, encoding the accessory protein swollenin, and the β-glucosidases cel1a, cel1b, cel3b, and cel3g increases in the early stages of cultivation. The increased production of these β-glucosidases increases the hydrolysis rate of cellobiose and sophorose, which activates carbon catabolite repression (CCR) and causes repression of cellulase genes and the regulator Xyr1 in the later stages of cultivation. Moreover, overexpression of azf1 led to increased cellulase activity in T. reesei during long-term cultivation in cellulose and sugarcane bagasse. Our results provide new insights into the mechanisms regulating Azf1 and novel genes that are important targets of this TF. This work contributes to a better understanding of the complex mechanisms regulating cellulase expression in T. reesei. It will contribute to the development of strains with higher production of these essential enzymes.

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

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

Structural Evolution of Gene Promoters Driven by Primate-Specific KRAB Zinc Finger Proteins

Krüppel-associated box (KRAB) zinc finger proteins (KZNFs) recognize and repress transposable elements (TEs); TEs are DNA elements that are capable of replicating themselves throughout our genomes with potentially harmful consequences. However, genes from this family of transcription factors have a much wider potential for genomic regulation. KZNFs have become integrated into gene-regulatory networks through the control of TEs that function as enhancers and gene promoters; some KZNFs also bind directly to gene promoters, suggesting an additional, more direct layer of KZNF co-option into gene-regulatory networks. Binding site analysis of ZNF519, ZNF441, and ZNF468 suggests the structural evolution of KZNFs to recognize TEs can result in coincidental binding to gene promoters independent of TE sequences. We show a higher rate of sequence turnover in gene promoter KZNF binding sites than neighboring regions, implying a selective pressure is being applied by the binding of a KZNF. Through CRISPR/Cas9 mediated genetic deletion of ZNF519, ZNF441, and ZNF468, we provide further evidence for genome-wide co-option of the KZNF-mediated gene-regulatory functions; KZNF knockout leads to changes in expression of KZNF-bound genes in neuronal lineages. Finally, we show that the opposite can be established upon KZNF overexpression, further strengthening the support for the role of KZNFs as bona-fide gene regulators. With no eminent role for ZNF519 in controlling its TE target, our study may provide a snapshot into the early stages of the completed co-option of a KZNF, showing the lasting, multilayered impact that retrovirus invasions and host response mechanisms can have upon the evolution of our genomes.

Unraveling the role of ZNF506 as a human PBS-pro-targeting protein for ERVP repression

Krüppel-associated box zinc finger proteins (KZFPs) function as a defense mechanism to maintain the genome stability of higher vertebrates by regulating the transcriptional activities of transposable elements (TEs). While previous studies have characterized ZFP809 as responsible for binding and repressing ERVs containing a proline tRNA primer-binding site (PBS-Pro) in mice, comparable KZFPs have not been identified in humans yet. Here, we identified ZNF506 as a PBS-Pro-binding protein in humans, which functions as a transcriptional repressor of PBS-Pro-utilizing retroviruses by recruiting heterochromatic modifications. Although they have similar functions, the low protein similarities between ZNF506, ZFP809 and KZFPs of other species suggest their independent evolution against the invasion of PBS-Pro-utilizing retroviruses into their respective ancestor genomes after species divergence. We also explored the link between ZNF506 and leukemia. Our findings suggest that ZNF506 is a unique human KZFP that can bind to PBS-Pro, highlighting the diverse evolution of KZFPs in defending against retroviral invasions. Additionally, our study provides insights into the potential role of ZNF506 in leukemia, contributing to the expanding knowledge of KZFPs' crucial function in disease and genome stability.

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.

ZBTB34 is a hepatocellular carcinoma-associated protein with a monopartite nuclear localization signal

ZBTB34 is a novel zinc finger protein with an unknown function. In this study, the gene expression and survival prognosis of ZBTB34 were analyzed across tumors based on the TCGA datasets. According to the bioinformatics analysis and qPCR results, liver hepatocellular carcinomas exhibit a high level of ZBTB34 expression. Additionally, the experiment supported the bioinformatics analysis findings that ZBTB34 expression was regulated by miR-125b-5p and that ZBTB34 affected ZBTB10, POLR1B, and AUH expression in HepG2 cells. Biological software analysis further revealed that ZBTB34 contains a monopartite nuclear localization signal (NLS). Arginine and lysine inside the putative NLS were substituted using the alanine-scanning mutagenesis method. The findings showed that the ability of ZBTB34 to enter the nucleus was abolished by the alanine substitution of the sequence 320RGGRARQKRA329 and the mutation of Lys327 and Arg328 residues. ZBTB34 was co-immunoprecipitated with importin α1, importin α3, importin α4, and importin β1, according to the results of the co-immunoprecipitation assay. In conclusion, ZBTB34 is a hepatocellular carcinoma-associated protein with a monopartite NLS. The nuclear import of ZBTB34 is mediated by importin α1, importin α3, importin α4, and importin β1. ZBTB34 performs its biological functions via a putative miR-125b-5p/ZBTB34/(ZBTB10, POLR1B, and AUH) signaling axis in HepG2 cells.

Transposable elements drive the evolution of metazoan zinc finger genes

Cys2-His2 zinc finger genes (ZNFs) form the largest family of transcription factors in metazoans. ZNF evolution is highly dynamic and characterized by the rapid expansion and contraction of numerous subfamilies across the animal phylogeny. The forces and mechanisms underlying rapid ZNF evolution remain poorly understood, but there is growing evidence that, in tetrapods, the targeting and repression of lineage-specific transposable elements (TEs) plays a critical role in the evolution of the Krüppel-associated box ZNF (KZNF) subfamily. Currently, it is unknown whether this function and coevolutionary relationship is unique to KZNFs or is a broader feature of metazoan ZNFs. Here, we present evidence that genomic conflict with TEs has been a central driver of the diversification of ZNFs in animals. Sampling from 3221 genome assemblies, we show that the copy number of retroelements correlates with that of ZNFs across at least 750 million years of metazoan evolution. Using computational predictions, we show that ZNFs preferentially bind TEs in diverse animal species. We further investigate the largest ZNF subfamily found in cyprinid fish, which is characterized by a conserved sequence we dubbed the fish N-terminal zinc finger-associated (FiNZ) domain. Zebrafish possess approximately 700 FiNZ-ZNFs, many of which are evolving adaptively under positive selection. Like mammalian KZNFs, most zebrafish FiNZ-ZNFs are expressed at the onset of zygotic genome activation, and blocking their translation using morpholinos during early embryogenesis results in derepression of transcriptionally active TEs. Together, these data suggest that ZNF diversification has been intimately connected to TE expansion throughout animal evolution.

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

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

Patterns of recombination in snakes reveal a tug of war between PRDM9 and promoter-like features

In vertebrates, there are two known mechanisms by which meiotic recombination is directed to the genome: in humans, mice, and other mammals, recombination occurs almost exclusively where the protein PRDM9 binds, while in species lacking an intact PRDM9, such as birds and canids, recombination rates are elevated near promoter-like features. To test if PRDM9 also directs recombination in non-mammalian vertebrates, we focused on an exemplar species, the corn snake (Pantherophis guttatus). Unlike birds, this species possesses a single, intact PRDM9 ortholog. By inferring historical recombination rates along the genome from patterns of linkage disequilibrium and identifying crossovers in pedigrees, we found that PRDM9 specifies the location of recombination events outside of mammals. However, we also detected an independent effect of promoter-like features on recombination, which is more pronounced on macro- than microchromosomes. Thus, our findings reveal that the uses of PRDM9 and promoter-like features are not mutually-exclusive, and instead reflect a tug of war, which varies in strength along the genome and is more lopsided in some species than others.

Crol contributes to PRE-mediated repression and Polycomb group proteins recruitment in Drosophila

The Polycomb group (PcG) proteins are fundamental epigenetic regulators that control the repressive state of target genes in multicellular organisms. One of the open questions is defining the mechanisms of PcG recruitment to chromatin. In Drosophila, the crucial role in PcG recruitment is thought to belong to DNA-binding proteins associated with Polycomb response elements (PREs). However, current data suggests that not all PRE-binding factors have been identified. Here, we report the identification of the transcription factor Crooked legs (Crol) as a novel PcG recruiter. Crol is a C2H2-type Zinc Finger protein that directly binds to poly(G)-rich DNA sequences. Mutation of Crol binding sites as well as crol CRISPR/Cas9 knockout diminish the repressive activity of PREs in transgenes. Like other PRE-DNA binding proteins, Crol co-localizes with PcG proteins inside and outside of H3K27me3 domains. Crol knockout impairs the recruitment of the PRC1 subunit Polyhomeotic and the PRE-binding protein Combgap at a subset of sites. The decreased binding of PcG proteins is accompanied by dysregulated transcription of target genes. Overall, our study identified Crol as a new important player in PcG recruitment and epigenetic regulation.

Prediction of cooperative homeodomain DNA binding sites from high-throughput-SELEX data

Homeodomain proteins constitute one of the largest families of metazoan transcription factors. Genetic studies have demonstrated that homeodomain proteins regulate many developmental processes. Yet, biochemical data reveal that most bind highly similar DNA sequences. Defining how homeodomain proteins achieve DNA binding specificity has therefore been a long-standing goal. Here, we developed a novel computational approach to predict cooperative dimeric binding of homeodomain proteins using High-Throughput (HT) SELEX data. Importantly, we found that 15 of 88 homeodomain factors form cooperative homodimer complexes on DNA sites with precise spacing requirements. Approximately one third of the paired-like homeodomain proteins cooperatively bind palindromic sequences spaced 3 bp apart, whereas other homeodomain proteins cooperatively bind sites with distinct orientation and spacing requirements. Combining structural models of a paired-like factor with our cooperativity predictions identified key amino acid differences that help differentiate between cooperative and non-cooperative factors. Finally, we confirmed predicted cooperative dimer sites in vivo using available genomic data for a subset of factors. These findings demonstrate how HT-SELEX data can be computationally mined to predict cooperativity. In addition, the binding site spacing requirements of select homeodomain proteins provide a mechanism by which seemingly similar AT-rich DNA sequences can preferentially recruit specific homeodomain factors.

On the dependent recognition of some long zinc finger proteins

The human genome contains about 800 C2H2 zinc finger proteins (ZFPs), and most of them are composed of long arrays of zinc fingers. Standard ZFP recognition model asserts longer finger arrays should recognize longer DNA-binding sites. However, recent experimental efforts to identify in vivo ZFP binding sites contradict this assumption, with many exhibiting short motifs. Here we use ZFY, CTCF, ZIM3, and ZNF343 as examples to address three closely related questions: What are the reasons that impede current motif discovery methods? What are the functions of those seemingly unused fingers and how can we improve the motif discovery algorithms based on long ZFPs' biophysical properties? Using ZFY, we employed a variety of methods and find evidence for 'dependent recognition' where downstream fingers can recognize some previously undiscovered motifs only in the presence of an intact core site. For CTCF, high-throughput measurements revealed its upstream specificity profile depends on the strength of its core. Moreover, the binding strength of the upstream site modulates CTCF's sensitivity to different epigenetic modifications within the core, providing new insight into how the previously identified intellectual disability-causing and cancer-related mutant R567W disrupts upstream recognition and deregulates the epigenetic control by CTCF. Our results establish that, because of irregular motif structures, variable spacing and dependent recognition between sub-motifs, the specificities of long ZFPs are significantly underestimated, so we developed an algorithm, ModeMap, to infer the motifs and recognition models of ZIM3 and ZNF343, which facilitates high-confidence identification of specific binding sites, including repeats-derived elements. With revised concept, technique, and algorithm, we can discover the overlooked specificities and functions of those 'extra' fingers, and therefore decipher their broader roles in human biology and diseases.

ZFP92, a KRAB domain zinc finger protein enriched in pancreatic islets, binds to B1/Alu SINE transposable elements and regulates retroelements and genes

Repressive KRAB domain-containing zinc-finger proteins (KRAB-ZFPs) are abundant in mammalian genomes and contribute both to the silencing of transposable elements (TEs) and to the regulation of developmental stage- and cell type-specific gene expression. Here we describe studies of zinc finger protein 92 (Zfp92), an X-linked KRAB-ZFP that is highly expressed in pancreatic islets of adult mice, by analyzing global Zfp92 knockout (KO) mice. Physiological, transcriptomic and genome-wide chromatin binding studies indicate that the principal function of ZFP92 in mice is to bind to and suppress the activity of B1/Alu type of SINE elements and modulate the activity of surrounding genomic entities. Deletion of Zfp92 leads to changes in expression of select LINE and LTR retroelements and genes located in the vicinity of ZFP92-bound chromatin. The absence of Zfp92 leads to altered expression of specific genes in islets, adipose and muscle that result in modest sex-specific alterations in blood glucose homeostasis, body mass and fat accumulation. In islets, Zfp92 influences blood glucose concentration in postnatal mice via transcriptional effects on Mafb, whereas in adipose and muscle, it regulates Acacb, a rate-limiting enzyme in fatty acid metabolism. In the absence of Zfp92, a novel TE-Capn11 fusion transcript is overexpressed in islets and several other tissues due to de-repression of an IAPez TE adjacent to ZFP92-bound SINE elements in intron 3 of the Capn11 gene. Together, these studies show that ZFP92 functions both to repress specific TEs and to regulate the transcription of specific genes in discrete tissues.

Structural insights into the recognition of telomeric variant repeat TTGGGG by broad-complex, tramtrack and bric-à-brac - zinc finger protein ZBTB10

Multiple proteins bind to telomeric DNA and are important for the role of telomeres in genome stability. A recent study established a broad-complex, tramtrack and bric-à-brac - zinc finger (BTB-ZF) protein, ZBTB10 (zinc finger and BTB domain-containing protein 10), as a telomeric variant repeat-binding protein at telomeres that use an alternative method for lengthening telomeres). ZBTB10 specifically interacts with the double-stranded telomeric variant repeat sequence TTGGGG by employing its tandem C2H2 zinc fingers (ZF1-2). Here, we solved the crystal structure of human ZBTB10 ZF1-2 in complex with a double-stranded DNA duplex containing the sequence TTGGGG to assess the molecular details of this interaction. Combined with calorimetric analysis, we identified the vital residues in TTGGGG recognition and determined the specific recognition mechanisms that are different from those of TZAP (telomere zinc finger-associated protein), a recently defined telomeric DNA-binding protein. Following these studies, we further identified a single amino-acid mutant (Arg767Gln) of ZBTB10 ZF1-2 that shows a preference for the telomeric DNA repeat TTAGGG sequence. We solved the cocrystal structure, providing a structural basis for telomeric DNA recognition by C2H2 ZF proteins.

The Essential Role of Prolines and Their Conformation in Allosteric Regulation of Kaiso Zinc Finger DNA-Binding Activity by the Adjacent C-Terminal Loop

Kaiso is a methyl-DNA-binding protein containing three C2H2 zinc fingers with a C-terminal extension that participates in DNA binding. The linker between the last zinc finger and the DNA-binding portion of the extension contains two prolines that are highly conserved in vertebrates and in cognate ZBTB4 and ZBTB38 proteins. Prolines provide chain rigidity and can exist in cis and trans conformations that can be switched by proline isomerases, affecting protein function. We found that substitution of the conserved proline P588, but not of P577, to alanine, negatively affected KaisoDNA-binding according to molecular dynamics simulation and in vitro DNA-binding assays. Molecular dynamics simulations of the Kaiso DNA-binding domain with P588 either substituted to alanine or switched to the cis-conformation revealed similar alterations in the H-bonding network and uncovered allosteric effects leading to structural rearrangements in the entire domain that resulted in the weakening of DNA-binding affinity. The substitution of proline with a large hydrophobic residue led to the same negative effects despite its ability to partially rescue the intrinsic DNA-binding activity of the C-terminal loop. Thus, the presence of the C-terminal extension and cis-conformation of proline residues are essential for efficient Kaiso-DNA binding, which likely involves intramolecular tension squeezing the DNA chain.

The chromatin factor ROW cooperates with BEAF-32 in regulating long-range inducible genes

Insulator proteins located at the boundaries of topological associated domains (TAD) are involved in higher-order chromatin organization and transcription regulation. However, it is still not clear how long-range contacts contribute to transcriptional regulation. Here, we show that relative-of-WOC (ROW) is essential for the long-range transcription regulation mediated by the boundary element-associated factor of 32kD (BEAF-32). We find that ROW physically interacts with heterochromatin proteins (HP1b and HP1c) and the insulator protein (BEAF-32). These proteins interact at TAD boundaries where ROW, through its AT-hook motifs, binds AT-rich sequences flanked by BEAF-32-binding sites and motifs. Knockdown of row downregulates genes that are long-range targets of BEAF-32 and bound indirectly by ROW (without binding motif). Analyses of high-throughput chromosome conformation capture (Hi-C) data reveal long-range interactions between promoters of housekeeping genes bound directly by ROW and promoters of developmental genes bound indirectly by ROW. Thus, our results show cooperation between BEAF-32 and the ROW complex, including HP1 proteins, to regulate the transcription of developmental and inducible genes through long-range interactions.

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

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

The genes within our DNA encode the essentials of our body plan and how each task in the body is achieved. However, our genome also contains many repetitive regions of DNA that do not encode functional genes. Some of these regions are genetic parasites known as transposons that try to multiply and spread around the DNA of their host. To prevent transposon DNA from interfering with the way the body operates, humans and other animals have evolved elaborate defense mechanisms to identify transposons and prevent them from multiplying.

In one such mechanism, known as the piRNA pathway, the host makes small molecules known as piRNAs that have sequences complementary to those of transposons, and act as guides to silence the transposons. The instructions to make these piRNAs are stored in the form of transposon fragments in dedicated regions of host DNA called piRNA clusters. These clusters thereby act as genetic memory, allowing the host to recognize and silence specific transposons in other locations within the host’s genome. In fruit flies, a protein called Rhino binds to piRNA clusters that are densely packed to allow piRNAs to be made. However, it remained unclear how Rhino is able to identify and bind to piRNA clusters, but not to other similarly densely packed regions of DNA.

Baumgartner et al. used a combination of genetic, genomic, and imaging approaches to study how Rhino finds its way in the fruit fly genome. They found that another protein called Kipferl interacts with Rhino and is required for Rhino to bind to nearly all piRNA clusters. Since Kipferl can by itself bind to the sequences that Rhino needs to find, the results suggest that Kipferl acts to recruit and initiate Rhino binding within densely packed piRNA clusters. Further experiments found that, in flies lacking Kipferl, Rhino binds to regions of DNA called Satellite repeats, hinting that these selfish sequences may compete for Rhino for their own benefit.

The finding that Kipferl and Rhino work together to define the memory system of the piRNA pathway strongly advances our understanding of how a sequence-specific defense system based on small RNAs can be established.

1,8-Dihydroxynaphthalene Melanin Biosynthesis in Colletotrichum fructicola Is Developmentally Regulated and Requires the Cooperative Function of Two Putative Zinc Finger Transcription Factors

In ascomycetes, 1,8-dihydroxynaphthalene (DHN) melanin plays important protective functions and its production is usually coupled with development and environmental stress responses. The regulation of melanin biosynthesis, however, remains obscure. Colletotrichum fructicola is a phytopathogen with a broad host range that produces melanized appressoria and perithecia. In this study, we annotated melanin genes in a high-quality C. fructicola genome and characterized two zinc finger transcription factors (TFs) (cmr1 and cmr2) that form a loosely organized gene cluster with several melanin biosynthesis genes. Deleting either TF abolished melanization in both mycelia and perithecia but did not affect appressoria. The deletion mutants also showed perithecial development defects. Overexpressing cmr1 in Δcmr2 strongly activated the expression of melanin biosynthesis genes including pks1, scd1, t4hr1, and thr1 and caused hyper-accumulation of charcoal to black pigment(s). On the other hand, overexpressing cmr2 in Δcmr1 activated pks1, t4hr1, and thr1, but not scd1. We conclude that proper DHN melanin accumulation in C. fructicola requires the cooperative function of two in-cluster TFs that also regulate perithecial development. The study clarifies DHN melanin regulations in C. fructicola and expands the function of melanin in-cluster TFs to sex regulation.

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.

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.

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.

Evolution of the recombination regulator PRDM9 in minke whales

Background

PRDM9 is a key regulator of meiotic recombination in most metazoans, responsible for reshuffling parental genomes. During meiosis, the PRDM9 protein recognizes and binds specific target motifs via its array of C₂H₂ zinc-fingers encoded by a rapidly evolving minisatellite. The gene coding for PRDM9 is the only speciation gene identified in vertebrates to date and shows high variation, particularly in the DNA-recognizing positions of the zinc-finger array, within and between species. Across all vertebrate genomes studied for PRDM9 evolution, only one genome lacks variability between repeat types – that of the North Pacific minke whale. This study aims to understand the evolution and diversity of Prdm9 in minke whales, which display the most unusual genome reference allele of Prdm9 so far discovered in mammals.

Results

Minke whales possess all the features characteristic of PRDM9-directed recombination, including complete KRAB, SSXRD and SET domains and a rapidly evolving array of C₂H₂-type-Zincfingers (ZnF) with evidence of rapid evolution, particularly at DNA-recognizing positions that evolve under positive diversifying selection. Seventeen novel PRDM9 variants were identified within the Antarctic minke whale species, plus a single distinct PRDM9 variant in Common minke whales – shared across North Atlantic and North Pacific minke whale subspecies boundaries.

Conclusion

The PRDM9 ZnF array evolves rapidly, in minke whales, with at least one DNA-recognizing position under positive selection. Extensive PRDM9 diversity is observed, particularly in the Antarctic in minke whales. Common minke whales shared a specific Prdm9 allele across subspecies boundaries, suggesting incomplete speciation by the mechanisms associated with PRDM9 hybrid sterility.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08305-1.

Protein-DNA/RNA interactions: Machine intelligence tools and approaches in the era of artificial intelligence and big data

With the development of artificial intelligence technologies and the availability of large amounts of biological data, computational methods for proteomics have undergone a developmental process from traditional machine learning to deep learning. This review focuses on computational approaches and tools for the prediction of protein-DNA/RNA interactions using machine intelligence techniques. We provide an overview of the development progress of computational methods and summarize the advantages and shortcomings of these methods. We further compiled applications in tasks related to the protein-DNA/RNA interactions, and pointed out possible future application trends. Moreover, biological sequence-digitizing representation strategies used in different types of computational methods are also summarized and discussed. This article is protected by copyright. All rights reserved.

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.

Cataloging Human PRDM9 Allelic Variation Using Long-Read Sequencing Reveals PRDM9 Population Specificity and Two Distinct Groupings of Related Alleles

The PRDM9 protein determines sites of meiotic recombination in humans by directing meiotic DNA double-strand breaks to specific loci. Targeting specificity is encoded by a long array of C₂H₂ zinc fingers that bind to DNA. This zinc finger array is hypervariable, and the resulting alleles each have a potentially different DNA binding preference. The assessment of PRDM9 diversity is important for understanding the complexity of human population genetics, inheritance linkage patterns, and predisposition to genetic disease. Due to the repetitive nature of the PRDM9 zinc finger array, the large-scale sequencing of human PRDM9 is challenging. We, therefore, developed a long-read sequencing strategy to infer the diploid PRDM9 zinc finger array genotype in a high-throughput manner. From an unbiased study of PRDM9 allelic diversity in 720 individuals from seven human populations, we detected 69 PRDM9 alleles. Several alleles differ in frequency among human populations, and 32 alleles had not been identified by previous studies, which were heavily biased to European populations. PRDM9 alleles are distinguished by their DNA binding site preferences and fall into two major categories related to the most common PRDM9-A and PRDM9-C alleles. We also found that it is likely that inter-conversion between allele types is rare. By mapping meiotic double-strand breaks (DSBs) in the testis, we found that small variations in PRDM9 can substantially alter the meiotic recombination landscape, demonstrating that minor PRDM9 variants may play an under-appreciated role in shaping patterns of human recombination. In summary, our data greatly expands knowledge of PRDM9 diversity in humans.

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.

A cis-acting structural variation at the ZNF558 locus controls a gene regulatory network in human brain development

The human forebrain has expanded in size and complexity compared to chimpanzees despite limited changes in protein-coding genes, suggesting that gene expression regulation is an important driver of brain evolution. Here, we identify a KRAB-ZFP transcription factor, ZNF558, that is expressed in human but not chimpanzee forebrain neural progenitor cells. ZNF558 evolved as a suppressor of LINE-1 transposons but has been co-opted to regulate a single target, the mitophagy gene SPATA18. ZNF558 plays a role in mitochondrial homeostasis, and loss-of-function experiments in cerebral organoids suggests that ZNF558 influences developmental timing during early human brain development. Expression of ZNF558 is controlled by the size of a variable number tandem repeat that is longer in chimpanzees compared to humans, and variable in the human population. Thus, this work provides mechanistic insight into how a cis-acting structural variation establishes a regulatory network that affects human brain evolution.

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.

Genome organization and evolution of a eukaryotic nicotinate co-inducible pathway

In Aspergillus nidulans a regulon including 11 hxn genes (hxnS, T, R, P, Y, Z, X, W, V, M and N) is inducible by a nicotinate metabolic derivative, repressible by ammonium and under stringent control of the nitrogen-state-sensitive GATA factor AreA and the specific transcription factor HxnR. This is the first report in a eukaryote of the genomic organization of a possibly complete pathway of nicotinate utilization. In A. nidulans the regulon is organized in three distinct clusters, this organization is variable in the Ascomycota. In some Pezizomycotina species all 11 genes map in a single cluster; in others they map in two clusters. This variable organization sheds light on cluster evolution. Instances of gene duplication followed by or simultaneous with integration in the cluster, partial or total cluster loss, and horizontal gene transfer of several genes (including an example of whole cluster re-acquisition in Aspergillus of section Flavi) were detected, together with the incorporation in some clusters of genes not found in the A. nidulans co-regulated regulon, which underlie both the plasticity and the reticulate character of metabolic cluster evolution. This study provides a comprehensive phylogeny of six members of the cluster across representatives of all Ascomycota classes.

Identification of the transcription factor MAZ as a regulator of erythropoiesis

Erythropoiesis requires a combination of ubiquitous and tissue-specific transcription factors (TFs). Here, through DNA affinity purification followed by mass spectrometry, we have identified the widely expressed protein MAZ (Myc-associated zinc finger) as a TF that binds to the promoter of the erythroid-specific human α-globin gene. Genome-wide mapping in primary human erythroid cells revealed that MAZ also occupies active promoters as well as GATA1-bound enhancer elements of key erythroid genes. Consistent with an important role during erythropoiesis, knockdown of MAZ reduces α-globin expression in K562 cells and impairs differentiation in primary human erythroid cells. Genetic variants in the MAZ locus are associated with changes in clinically important human erythroid traits. Taken together, these findings reveal the zinc-finger TF MAZ to be a previously unrecognized regulator of the erythroid differentiation program.

Genetic variation in recombination rate in the pig

Background

Meiotic recombination results in the exchange of genetic material between homologous chromosomes. Recombination rate varies between different parts of the genome, between individuals, and is influenced by genetics. In this paper, we assessed the genetic variation in recombination rate along the genome and between individuals in the pig using multilocus iterative peeling on 150,000 individuals across nine genotyped pedigrees. We used these data to estimate the heritability of recombination and perform a genome-wide association study of recombination in the pig.

Results

Our results confirmed known features of the recombination landscape of the pig genome, including differences in genetic length of chromosomes and marked sex differences. The recombination landscape was repeatable between lines, but at the same time, there were differences in average autosome-wide recombination rate between lines. The heritability of autosome-wide recombination rate was low but not zero (on average 0.07 for females and 0.05 for males). We found six genomic regions that are associated with recombination rate, among which five harbour known candidate genes involved in recombination: RNF212, SHOC1, SYCP2, MSH4 and HFM1.

Conclusions

Our results on the variation in recombination rate in the pig genome agree with those reported for other vertebrates, with a low but nonzero heritability, and the identification of a major quantitative trait locus for recombination rate that is homologous to that detected in several other species. This work also highlights the utility of using large-scale livestock data to understand biological processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00643-0.

The duality of PRDM proteins: epigenetic and structural perspectives

PRDF1 and RIZ1 homology domain containing (PRDMs) are a subfamily of Krüppel-like zinc finger proteins controlling key processes in metazoan development and in cancer. PRDMs exhibit unique dualities: (a) PR domain/ZNF arrays-their structure combines a SET-like domain known as a PR domain, typically found in methyltransferases, with a variable array of C2H2 zinc fingers (ZNF) characteristic of DNA-binding transcription factors; (b) transcriptional activators/repressors-their physiological function is context- and cell-dependent; mechanistically, some PRDMs have a PKMT activity and directly catalyze histone lysine methylation, while others are rather pseudomethyltransferases and act by recruiting transcriptional cofactors; (c) oncogenes/tumor suppressors-their pathological function depends on the specific PRDM isoform expressed during tumorigenesis. This duality is well known as the 'Yin and Yang' of PRDMs and involves a complex regulation of alternative splicing or alternative promoter usage, to generate full-length or PR-deficient isoforms with opposing functions in cancer. In conclusion, once their dualities are fully appreciated, PRDMs represent a promising class of targets in oncology by virtue of their widespread upregulation across multiple tumor types and their somatic dispensability, conferring a broad therapeutic window and limited toxic side effects. The recent discovery of a first-in-class compound able to inhibit PRDM9 activity has paved the way for the identification of further small molecular inhibitors able to counteract PRDM oncogenic activity.