Multiparameter functional diversity of human C2H2 zinc finger proteins

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.

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.

Identifying human pre-mRNA cleavage and polyadenylation factors by genome-wide CRISPR screens using a dual fluorescence readthrough reporter

Messenger RNA precursors (pre-mRNA) generally undergo 3' end processing by cleavage and polyadenylation (CPA), which is specified by a polyadenylation site (PAS) and adjacent RNA sequences and regulated by a large variety of core and auxiliary CPA factors. To date, most of the human CPA factors have been discovered through biochemical and proteomic studies. However, genetic identification of the human CPA factors has been hampered by the lack of a reliable genome-wide screening method. We describe here a dual fluorescence readthrough reporter system with a PAS inserted between two fluorescent reporters. This system enables measurement of the efficiency of 3' end processing in living cells. Using this system in combination with a human genome-wide CRISPR/Cas9 library, we conducted a screen for CPA factors. The screens identified most components of the known core CPA complexes and other known CPA factors. The screens also identified CCNK/CDK12 as a potential core CPA factor, and RPRD1B as a CPA factor that binds RNA and regulates the release of RNA polymerase II at the 3' ends of genes. Thus, this dual fluorescence reporter coupled with CRISPR/Cas9 screens reliably identifies bona fide CPA factors and provides a platform for investigating the requirements for CPA in various contexts.

Single-mitosis dissection of acute and chronic DNA mutagenesis and repair

How chronic mutational processes and punctuated bursts of DNA damage drive evolution of the cancer genome is poorly understood. Here, we demonstrate a strategy to disentangle and quantify distinct mechanisms underlying genome evolution in single cells, during single mitoses and at single-strand resolution. To distinguish between chronic (reactive oxygen species (ROS)) and acute (ultraviolet light (UV)) mutagenesis, we microfluidically separate pairs of sister cells from the first mitosis following burst UV damage. Strikingly, UV mutations manifest as sister-specific events, revealing mirror-image mutation phasing genome-wide. In contrast, ROS mutagenesis in transcribed regions is reduced strand agnostically. Successive rounds of genome replication over persisting UV damage drives multiallelic variation at CC dinucleotides. Finally, we show that mutation phasing can be resolved to single strands across the entire genome of liver tumors from F1 mice. This strategy can be broadly used to distinguish the contributions of overlapping cancer relevant mutational processes.

The kruppel-like factor (KLF) family, diseases, and physiological events

The Kruppel-Like Factor family of regulatory proteins, which has 18 members, is transcription factors. This family contains zinc finger proteins, regulates the activation and suppression of transcription, and binds to DNA, RNA, and proteins. Klfs related to the immune system are Klf1, Klf2, Klf3, Klf4, Klf6, and Klf14. Klfs related to adipose tissue development and/or glucose metabolism are Klf3, Klf7, Klf9, Klf10, Klf11, Klf14, Klf15, and Klf16. Klfs related to cancer are Klf3, Klf4, Klf5, Klf6, Klf7, Klf8, Klf9, Klf10, Klf11, Klf12, Klf13, Klf14, Klf16, and Klf17. Klfs related to the cardiovascular system are Klf4, Klf5, Klf10, Klf13, Klf14, and Klf15. Klfs related to the nervous system are Klf4, Klf7, Klf8, and Klf9. Klfs are associated with diseases such as carcinogenesis, oxidative stress, diabetes, liver fibrosis, thalassemia, and the metabolic syndrome. The aim of this review is to provide information about the relationship of Klfs with some diseases and physiological events and to guide future studies.

Zinc Finger and BTB Domain-Containing 20: A Newly Emerging Player in Pathogenesis and Development of Human Cancers

Zinc finger and BTB domain-containing 20 (ZBTB20), which was initially identified in human dendritic cells, belongs to a family of transcription factors (TFs) with an N-terminal BTB domain and one or more C-terminal DNA-binding zinc finger domains. Under physiological conditions, ZBTB20 acts as a transcriptional repressor in cellular development and differentiation, metabolism, and innate immunity. Interestingly, multiple lines of evidence from mice and human systems have revealed the importance of ZBTB20 in the pathogenesis and development of cancers. ZBTB20 is not only a hotspot of genetic variation or fusion in many types of human cancers, but also a key TF or intermediator involving in the dysregulation of cancer cells. Given the diverse functions of ZBTB20 in both health and disease, we herein summarize the structure and physiological roles of ZBTB20, with an emphasis on the latest findings on tumorigenesis and cancer progression.

Single-cell multi-omics profiling of human preimplantation embryos identifies cytoskeletal defects during embryonic arrest

Human in vitro fertilized embryos exhibit low developmental capabilities, and the mechanisms that underlie embryonic arrest remain unclear. Here using a single-cell multi-omics sequencing approach, we simultaneously analysed alterations in the transcriptome, chromatin accessibility and the DNA methylome in human embryonic arrest due to unexplained reasons. Arrested embryos displayed transcriptome disorders, including a distorted microtubule cytoskeleton, increased genomic instability and impaired glycolysis, which were coordinated with multiple epigenetic reprogramming defects. We identified Aurora A kinase (AURKA) repression as a cause of embryonic arrest. Mechanistically, arrested embryos induced through AURKA inhibition resembled the reprogramming abnormalities of natural embryonic arrest in terms of the transcriptome, the DNA methylome, chromatin accessibility and H3K4me3 modifications. Mitosis-independent sequential activation of the zygotic genome in arrested embryos showed that YY1 contributed to human major zygotic genome activation. Collectively, our study decodes the reprogramming abnormalities and mechanisms of human embryonic arrest and the key regulators of zygotic genome activation.

ZNF283, a Krüppel-associated box zinc finger protein, inhibits RNA synthesis of porcine reproductive and respiratory syndrome virus by interacting with Nsp9 and Nsp10

Porcine reproductive and respiratory syndrome virus (PRRSV) is a viral pathogen with substantial economic implications for the global swine industry. The existing vaccination strategies and antiviral drugs offer limited protection. Replication of the viral RNA genome encompasses a complex series of steps, wherein a replication complex is assembled from various components derived from both viral and cellular sources, as well as from the viral genomic RNA template. In this study, we found that ZNF283, a Krüppel-associated box (KRAB) containing zinc finger protein, was upregulated in PRRSV-infected Marc-145 cells and porcine alveolar macrophages and that ZNF283 inhibited PRRSV replication and RNA synthesis. We also found that ZNF283 interacts with the viral proteins Nsp9, an RNA-dependent RNA polymerase, and Nsp10, a helicase. The main regions involved in the interaction between ZNF283 and Nsp9 were determined to be the KRAB domain of ZNF283 and amino acids 178-449 of Nsp9. The KRAB domain of ZNF283 plays a role in facilitating Nsp10 binding. In addition, ZNF283 may have an affinity for the 3' untranslated region of PRRSV. These findings suggest that ZNF283 is an antiviral factor that inhibits PRRSV infection and extend our understanding of the interactions between KRAB-containing zinc finger proteins and viruses.

HOCOMOCO in 2024: a rebuild of the curated collection of binding models for human and mouse transcription factors

We present a major update of the HOCOMOCO collection that provides DNA binding specificity patterns of 949 human transcription factors and 720 mouse orthologs. To make this release, we performed motif discovery in peak sets that originated from 14 183 ChIP-Seq experiments and reads from 2554 HT-SELEX experiments yielding more than 400 thousand candidate motifs. The candidate motifs were annotated according to their similarity to known motifs and the hierarchy of DNA-binding domains of the respective transcription factors. Next, the motifs underwent human expert curation to stratify distinct motif subtypes and remove non-informative patterns and common artifacts. Finally, the curated subset of 100 thousand motifs was supplied to the automated benchmarking to select the best-performing motifs for each transcription factor. The resulting HOCOMOCO v12 core collection contains 1443 verified position weight matrices, including distinct subtypes of DNA binding motifs for particular transcription factors. In addition to the core collection, HOCOMOCO v12 provides motif sets optimized for the recognition of binding sites in vivo and in vitro, and for annotation of regulatory sequence variants. HOCOMOCO is available at https://hocomoco12.autosome.org and https://hocomoco.autosome.org.

Unraveling the role of miRNAs in the diagnosis, progression, and therapeutic intervention of Parkinson's disease

Parkinson's disease (PD) is a debilitating neurological disorder characterized by the impairment of the motor system, resulting in symptoms such as resting tremor, cogwheel rigidity, bradykinesia, difficulty with gait, and postural instability. The occurrence of striatal dopamine insufficiency can be attributed to a notable decline in dopaminergic neurons inside the substantia nigra pars compacta. Additionally, the development of Lewy bodies serves as a pathological hallmark of PD. While current therapy approaches for PD aim to preserve dopaminergic neurons or replenish dopamine levels in the brain, it is important to acknowledge that achieving complete remission of the condition remains elusive. MicroRNAs (miRNAs, miR) are a class of small, non-coding ribonucleic acids involved in regulating gene expression at the post-transcriptional level. The miRNAs play a crucial part in the underlying pathogenic mechanisms of several neurodegenerative illnesses, including PD. The aim of this review is to explore the role of miRNAs in regulating genes associated with the onset and progression of PD, investigate the potential of miRNAs as a diagnostic tool, assess the effectiveness of targeting specific miRNAs as an alternative therapeutic strategy to impede disease advancement, and discuss the utilization of newly developed nanoparticles for delivering miRNAs as neurodegenerative therapies.

Dimerization choice and alternative functions of ZBTB transcription factors

Zinc Finger DNA-binding domain-containing proteins are the most populous family among eukaryotic transcription factors. Among these, members of the BTB domain-containing ZBTB sub-family are mostly known for their transcriptional repressive functions. In this Viewpoint article, we explore molecular mechanisms that potentially diversify the function of ZBTB proteins based on their homo and heterodimerization, alternative splicing and post-translational modifications. We describe how the BTB domain is as much a scaffold for the assembly of co-repressors, as a domain that regulates protein stability. We highlight another mechanism that regulates ZBTB protein stability: phosphorylation in the zinc finger domain. We explore the non-transcriptional, structural roles of ZBTB proteins and highlight novel findings that describe the ability of ZBTB proteins to associate with poly adenosine ribose in the nucleus during the DNA damage response. Herein, we discuss the contribution of BTB domain scaffolds to the formation of transcriptional repressive complexes, to chromosome compartmentalization and their non-transcriptional, purely structural functions in the nucleus.

Essential transcription factors for induced neuron differentiation

Neurogenins are proneural transcription factors required to specify neuronal identity. Their overexpression in human pluripotent stem cells rapidly produces cortical-like neurons with spiking activity and, because of this, they have been widely adopted for human neuron disease models. However, we do not fully understand the key downstream regulatory effectors responsible for driving neural differentiation. Here, using inducible expression of NEUROG1 and NEUROG2, we identify transcription factors (TFs) required for directed neuronal differentiation by combining expression and chromatin accessibility analyses with a pooled in vitro CRISPR-Cas9 screen targeting all ~1900 TFs in the human genome. The loss of one of these essential TFs (ZBTB18) yields few MAP2-positive neurons. Differentiated ZBTB18-null cells have radically altered gene expression, leading to cytoskeletal defects and stunted neurites and spines. In addition to identifying key downstream TFs for neuronal differentiation, our work develops an integrative multi-omics and TFome-wide perturbation platform to rapidly characterize essential TFs for the differentiation of any human cell type.

ZNF432 stimulates PARylation and inhibits DNA resection to balance PARPi sensitivity and resistance

Zinc finger (ZNF) motifs are some of the most frequently occurring domains in the human genome. It was only recently that ZNF proteins emerged as key regulators of genome integrity in mammalian cells. In this study, we report a new role for the Krüppel-type ZNF-containing protein ZNF432 as a novel poly(ADP-ribose) (PAR) reader that regulates the DNA damage response. We show that ZNF432 is recruited to DNA lesions via DNA- and PAR-dependent mechanisms. Remarkably, ZNF432 stimulates PARP-1 activity in vitro and in cellulo. Knockdown of ZNF432 inhibits phospho-DNA-PKcs and increases RAD51 foci formation following irradiation. Moreover, purified ZNF432 preferentially binds single-stranded DNA and impairs EXO1-mediated DNA resection. Consequently, the loss of ZNF432 in a cellular system leads to resistance to PARP inhibitors while its overexpression results in sensitivity. Taken together, our results support the emerging concept that ZNF-containing proteins can modulate PARylation, which can be embodied by the pivotal role of ZNF432 to finely balance the outcome of PARPi response by regulating homologous recombination.

Take a walk on the KRAB side

Canonical Krüppel-associated box (KRAB)-containing zinc finger proteins (KZFPs) act as major repressors of transposable elements (TEs) via the KRAB-mediated recruitment of the heterochromatin scaffold KRAB-associated protein (KAP)1. KZFP genes emerged some 420 million years ago in the last common ancestor of coelacanth, lungfish, and tetrapods, and dramatically expanded to give rise to lineage-specific repertoires in contemporary species paralleling their TE load and turnover. However, the KRAB domain displays sequence and function variations that reveal repeated diversions from a linear TE-KZFP trajectory. This Review summarizes current knowledge on the evolution of KZFPs and discusses how ancestral noncanonical KZFPs endowed with variant KRAB, SCAN or DUF3669 domains have been utilized to achieve KAP1-independent functions.

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.

MAPK1 promotes the metastasis and invasion of gastric cancer as a bidirectional transcription factor

BACKGROUND

The Mitogen-activated protein kinase 1 (MAPK1) has both independent functions of phosphorylating histones as a kinase and directly binding the promoter regions of genes to regulate gene expression as a transcription factor. Previous studies have identified elevated expression of MAPK1 in human gastric cancer, which is associated with its role as a kinase, facilitating the migration and invasion of gastric cancer cells. However, how MAPK1 binds to its target genes as a transcription factor and whether it modulates related gene expressions in gastric cancer remains unclear.

RESULTS

Here, we integrated biochemical assays (protein interactions and chromatin immunoprecipitation (ChIP)), cellular analysis assays (cell proliferation and migration), RNA sequencing, ChIP sequencing, and clinical analysis to investigate the potential genomic recognition patterns of MAPK1 in a human gastric adenocarcinoma cell-line (AGS) and to uncover its regulatory effect on gastric cancer progression. We confirmed that MAPK1 promotes AGS cells invasion and migration by regulating the target genes in different directions, up-regulating seven target genes (KRT13, KRT6A, KRT81, MYH15, STARD4, SYTL4, and TMEM267) and down-regulating one gene (FGG). Among them, five genes (FGG, MYH15, STARD4, SYTL4, and TMEM267) were first associated with cancer procession, while the other three (KRT81, KRT6A, and KRT13) have previously been confirmed to be related to cancer metastasis and migration.

CONCLUSION

Our data showed that MAPK1 can bind to the promoter regions of these target genes to control their transcription as a bidirectional transcription factor, promoting AGS cell motility and invasion. Our research has expanded the understanding of the regulatory roles of MAPK1, enriched our knowledge of transcription factors, and provided novel candidates for cancer therapeutics.

Mechanism of herpesvirus UL24 protein regulating viral immune escape and virulence

Herpesviruses have evolved a series of abilities involved in the process of host infection that are conducive to virus survival and adaptation to the host, such as immune escape, latent infection, and induction of programmed cell death for sustainable infection. The herpesvirus gene UL24 encodes a highly conserved core protein that plays an important role in effective viral infection. The UL24 protein can inhibit the innate immune response of the host by acting on multiple immune signaling pathways during virus infection, and it also plays a key role in the proliferation and pathogenicity of the virus in the later stage of infection. This article reviews the mechanism by which the UL24 protein mediates herpesvirus immune escape and its effects on viral proliferation and virulence by influencing syncytial formation, DNA damage and the cell cycle. Reviewing these studies will enhance our understanding of the pathogenesis of herpesvirus infection and provide evidence for new strategies to combat against viral infection.

Integrated multi-omics analyses reveal Jorunnamycin A as a novel suppressor for muscle-invasive bladder cancer by targeting FASN and TOP1

BACKGROUND

Bladder cancer is a urological carcinoma with high incidence, among which muscle invasive bladder cancer (MIBC) is a malignant carcinoma with high mortality. There is an urgent need to develop new drugs with low toxicity and high efficiency for MIBC because existing medication has defects, such as high toxicity, poor efficacy, and side effects. Jorunnamycin A (JorA), a natural marine compound, has been found to have a high efficiency anticancer effect, but its anticancer function and mechanism on bladder cancer have not been studied.

METHODS

To examine the anticancer effect of JorA on MIBC, Cell Counting Kit 8, EdU staining, and colony formation analyses were performed. Moreover, a xenograft mouse model was used to verify the anticancer effect in vivo. To investigate the pharmacological mechanism of JorA, high-throughput quantitative proteomics, transcriptomics, RT-qPCR, western blotting, immunofluorescence staining, flow cytometry, pulldown assays, and molecular docking were performed.

RESULTS

JorA inhibited the proliferation of MIBC cells, and the IC50 of T24 and UM-UC-3 was 0.054 and 0.084 μM, respectively. JorA-induced significantly changed proteins were enriched in "cancer-related pathways" and "EGFR-related signaling pathways", which mainly manifested by inhibiting cell proliferation and promoting cell apoptosis. Specifically, JorA dampened the DNA synthesis rate, induced phosphatidylserine eversion, and inhibited cell migration. Furthermore, it was discovered that fatty acid synthase (FASN) and topoisomerase 1 (TOP1) are the JorA interaction proteins. Using DockThor software, the 3D docking structures of JorA binding to FASN and TOP1 were obtained (the binding affinities were - 8.153 and - 7.264 kcal/mol, respectively).

CONCLUSIONS

The marine compound JorA was discovered to have a specific inhibitory effect on MIBC, and its potential pharmacological mechanism was revealed for the first time. This discovery makes an important contribution to the development of new high efficiency and low toxicity drugs for bladder cancer therapy.

Genetic features and genomic targets of human KRAB-zinc finger proteins

Krüppel-associated box (KRAB) domain-containing zinc finger proteins (KZFPs) are one of the largest groups of transcription factors encoded by tetrapods, with 378 members in human alone. KZFP genes are often grouped in clusters reflecting amplification by gene and segment duplication since the gene family first emerged more than 400 million years ago. Previous work has revealed that many KZFPs recognize transposable element (TE)-embedded sequences as genomic targets, and that KZFPs facilitate the co-option of the regulatory potential of TEs for the benefit of the host. Here, we present a comprehensive survey of the genetic features and genomic targets of human KZFPs, notably completing past analyses by adding data on close to a hundred family members. General principles emerge from our study of the TE-KZFP regulatory system, which point to multipronged evolutionary mechanisms underlaid by highly complex and combinatorial modes of action with strong influences on human speciation.

An ancestral molecular response to nanomaterial particulates

The varied transcriptomic response to nanoparticles has hampered the understanding of the mechanism of action. Here, by performing a meta-analysis of a large collection of transcriptomics data from various engineered nanoparticle exposure studies, we identify common patterns of gene regulation that impact the transcriptomic response. Analysis identifies deregulation of immune functions as a prominent response across different exposure studies. Looking at the promoter regions of these genes, a set of binding sites for zinc finger transcription factors C2H2, involved in cell stress responses, protein misfolding and chromatin remodelling and immunomodulation, is identified. The model can be used to explain the outcomes of mechanism of action and is observed across a range of species indicating this is a conserved part of the innate immune system.

KLF12 overcomes anti-PD-1 resistance by reducing galectin-1 in cancer cells

BACKGROUNDS

Immune checkpoint blockade has revolutionized cancer treatment and has improved the survival of a subset of patients with cancer. However, numerous patients do not benefit from immunotherapy, and treatment resistance is a major challenge. Krüppel-like factor 12 (KLF12) is a transcriptional inhibitor whose role in tumor immunity is unclear.

METHODS

We demonstrated a relationship between KLF12 and CD8+ T cells in vivo and in vitro by flow cytometry. The role and underlying mechanism that KLF12 regulates CD8+ T cells were investigated using reverse transcription and quantitative PCR, western blot FACS, chromatin immunoprecipitation-PCR and Dual-Luciferase reporter assays, etc, and employing small interfering RNA (siRNA) and inhibitors. In vivo efficacy studies were conducted with multiple mouse tumor models, employing anti-programmed cell death protein 1 combined with KLF12 or galectin-1 (Gal-1) inhibitor.

RESULTS

Here, we found that the expression of tumor KLF12 correlates with immunotherapy resistance. KLF12 suppresses CD8+ T cells infiltration and function in vitro and in vivo. Mechanistically, KLF12 inhibits the expression of Gal-1 by binding with its promoter, thereby improving the infiltration and function of CD8+ T cells, which plays a vital role in cancer immunotherapy.

CONCLUSIONS

This work identifies a novel pathway regulating CD8+ T-cell intratumoral infiltration, and targeting the KLF12/Gal-1 axis may serve as a novel therapeutic target for patients with immunotherapy resistance.

The Novel Role of Zfp296 in Mammalian Embryonic Genome Activation as an H3K9me3 Modulator

The changes in epigenetic modifications during early embryonic development significantly impact mammalian embryonic genome activation (EGA) and are species-conserved to some degree. Here, we reanalyzed the published RNA-Seq of human, mouse, and goat early embryos and found that Zfp296 (zinc finger protein 296) expression was higher at the EGA stage than at the oocyte stage in all three species (adjusted p-value < 0.05 |log2(foldchange)| ≥ 1). Subsequently, we found that Zfp296 was conserved across human, mouse, goat, sheep, pig, and bovine embryos. In addition, we identified that ZFP296 interacts with the epigenetic regulators KDM5B, SMARCA4, DNMT1, DNMT3B, HP1β, and UHRF1. The Cys₂-His₂(C2H2) zinc finger domain TYPE2 TYPE3 domains of ZFP296 co-regulated the modification level of the trimethylation of lysine 9 on the histone H3 protein subunit (H3K9me3). According to ChIP-seq analysis, ZFP296 was also enriched in Trim28, Suv39h1, Setdb1, Kdm4a, and Ehmt2 in the mESC genome. Then, knockdown of the expression of Zfp296 at the late zygote of the mouse led to the early developmental arrest of the mouse embryos and failure resulting from a decrease in H3K9me3. Together, our results reveal that Zfp296 is an H3K9me3 modulator which is essential to the embryonic genome activation of mouse embryos.

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.

Inherited human ZNF341 deficiency

Typical hyper-IgE syndromes (HIES) are caused by autosomal-dominant-negative (DN) variants of STAT3 (Signal Transducer And Activator Of Transcription 3) or IL6ST (Interleukin 6 Cytokine Family Signal Transducer), biallelic partial loss-of-function (LOF) variants of IL6ST, or biallelic complete LOF variants of ZNF341 (Zinc Finger Protein 341). Including the two new cases described in this review, only 20 patients with autosomal-recessive (AR) ZNF341 deficiency have ever been reported. Patients with AR ZNF341 deficiency have clinical and immunological phenotypes resembling those of patients with autosomal-dominant STAT3 deficiency, but with a usually milder clinical presentation and lower NK (Natural Killer) cell counts. ZNF341-deficient cells have 50% the normal level of STAT3 in the resting state. However, as there is no clear evidence that STAT3 haploinsufficiency causes HIES, this decrease alone is probably insufficient to explain the HIES phenotype observed in the ZNF341-deficient patients. The combination of decreased basal expression level and impaired autoinduction of STAT3 observed in ZNF341-deficient lymphocytes is considered a more likely pathophysiological mechanism. We review here what is currently known about the ZNF341 gene and ZNF341 deficiency, and briefly discuss possible roles for this protein in addition to its control of STAT3 activity.

ZNF677 inhibits oral squamous cell carcinoma growth and tumor stemness by regulating FOXO3a

Oral squamous cell carcinoma (OSCC) is a common cancer with an increasing incidence worldwide. Zinc-finger proteins 677 (ZNF677) is involved in the progression and methylation of various cancers, but its role and mechanism in OSCC remain indeterminate. The expression of ZNF677 was analyzed by online database and immunohistochemistry, while the methylation level of ZNF677 was determined by the methylation-specific PCR. The role and mechanism of ZNF677 in the tumor cell growth, migration, invasion and stemness were addressed by cell counting kit-8 (CCK-8), 5-ethynyl-2'-deoxyuridine (EdU) incorporation, Transwell, wound-healing, sphere‑formation, and western blot assays. In addition, its function was also investigated in a xenografted mice model. The results showed that ZNF677 was lowly expressed in OSCC with a hypermethylation level, which predicted poor overall survival in patients with HNSC. Upregulation of ZNF677 reduced the cell viability, Edu positive cells, numbers of invasion cells, the migration ability, numbers of spheres formation and the expression of proliferation, migration and stemness related proteins in CAL-27 and SCC25 cells. Mechanically, the relative levels of p-AKT/AKT were decreased and the levels of p-FOXO3a/FOXO3a were increased in both cells overexpressed with ZNF677, which were reversed by the SC79 treatment. Moreover, interference of FOXO3a recovered the suppressive effects of ZNF677 overexpression on cell proliferation, migration, invasion and stemness of OSCC cells. Furthermore, overexpression of ZNF677 reduced the tumor volume and weight, and the relative protein level of p-AKT/AKT with an increased level of p-FOXO3a/FOXO3a, and improved pathological symptoms in vivo. Collectively, ZNF677 suppressed OSCC cells growth, migration, invasion and stemness through inhibiting AKT/FOXO3a pathway.

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.

Prdm6 drives ductus arteriosus closure by promoting ductus arteriosus smooth muscle cell identity and contractility

Based upon our demonstration that the smooth muscle cell-selective (SMC-selective) putative methyltransferase, Prdm6, interacts with myocardin-related transcription factor-A, we examined Prdm6's role in SMCs in vivo using cell type-specific knockout mouse models. Although SMC-specific depletion of Prdm6 in adult mice was well tolerated, Prdm6 depletion in Wnt1-expressing cells during development resulted in perinatal lethality and a completely penetrant patent ductus arteriosus (DA) phenotype. Lineage tracing experiments in Wnt1Cre2 Prdm6fl/fl ROSA26LacZ mice revealed normal neural crest-derived SMC investment of the outflow tract. In contrast, myography measurements on DA segments isolated from E18.5 embryos indicated that Prdm6 depletion significantly reduced DA tone and contractility. RNA-Seq analyses on DA and ascending aorta samples at E18.5 identified a DA-enriched gene program that included many SMC-selective contractile associated proteins that was downregulated by Prdm6 depletion. Chromatin immunoprecipitation-sequencing experiments in outflow tract SMCs demonstrated that 50% of the genes Prdm6 depletion altered contained Prdm6 binding sites. Finally, using several genome-wide data sets, we identified an SMC-selective enhancer within the Prdm6 third intron that exhibited allele-specific activity, providing evidence that rs17149944 may be the causal SNP for a cardiovascular disease GWAS locus identified within the human PRDM6 gene.

Sibling rivalry among the ZBTB transcription factor family: homodimers versus heterodimers

BTB domains potentially can form homo- or heterodimers. The study examines the dimerization choice of several BTB domains and finds only one heterodimer, while all tested pairs can homodimerize.

The BTB domain is an oligomerization domain found in over 300 proteins encoded in the human genome. In the family of BTB domain and zinc finger–containing (ZBTB) transcription factors, 49 members share the same protein architecture. The N-terminal BTB domain is structurally conserved among the family members and serves as the dimerization site, whereas the C-terminal zinc finger motifs mediate DNA binding. The available BTB domain structures from this family reveal a natural inclination for homodimerization. In this study, we investigated the potential for heterodimer formation in the cellular environment. We selected five BTB homodimers and four heterodimer structures. We performed cell-based binding assays with fluorescent protein–BTB domain fusions to assess dimer formation. We tested the binding of several BTB pairs, and we were able to confirm the heterodimeric physical interaction between the BTB domains of PATZ1 and PATZ2, previously reported only in an interactome mapping experiment. We also found this pair to be co-expressed in several immune system cell types. Finally, we used the available structures of BTB domain dimers and newly constructed models in extended molecular dynamics simulations (500 ns) to understand the energetic determinants of homo- and heterodimer formation. We conclude that heterodimer formation, although frequently described as less preferred than homodimers, is a possible mechanism to increase the combinatorial specificity of this transcription factor family.

Exploring the genetic space of the DNA damage response for cancer therapy through CRISPR-based screens

The concepts of synthetic lethality and viability have emerged as powerful approaches to identify vulnerabilities and resistances within the DNA damage response for the treatment of cancer. Historically, interactions between two genes have had a longstanding presence in genetics and have been identified through forward genetic screens that rely on the molecular basis of the characterized phenotypes, typically caused by mutations in single genes. While such complex genetic interactions between genes have been studied extensively in model organisms, they have only recently been prioritized as therapeutic strategies due to technological advancements in genetic screens. Here, we discuss synthetic lethal and viable interactions within the DNA damage response and present how CRISPR-based genetic screens and chemical compounds have allowed for the systematic identification and targeting of such interactions for the treatment of cancer.

A satellite DNA array barcodes chromosome 7 and regulates totipotency via ZFP819

Mammalian genomes are a battleground for genetic conflict between repetitive elements and KRAB-zinc finger proteins (KZFPs). We asked whether KZFPs can regulate cell fate by using ZFP819, which targets a satellite DNA array, ZP3AR. ZP3AR coats megabase regions of chromosome 7 encompassing genes encoding ZSCAN4, a master transcription factor of totipotency. Depleting ZFP819 in mouse embryonic stem cells (mESCs) causes them to transition to a 2-cell (2C)-like state, whereby the ZP3AR array switches from a poised to an active enhancer state. This is accompanied by a global erosion of heterochromatin roadblocks, which we link to decreased SETDB1 stability. These events result in transcription of active LINE-1 elements and impaired differentiation. In summary, ZFP819 and TRIM28 partner up to close chromatin across Zscan4, to promote exit from totipotency. We propose that satellite DNAs may control developmental fate transitions by barcoding and switching off master transcription factor genes.

Snail maintains the stem/progenitor state of skin epithelial cells and carcinomas through the autocrine effect of matricellular protein Mindin

Preservation of a small population of cancer stem cells (CSCs) within a heterogeneous carcinoma serves as a paradigm to understand how select cells in a tissue maintain their undifferentiated status. In both embryogenesis and cancer, Snail has been correlated with stemness, but the molecular underpinning of this phenomenon remains largely ill-defined. In models of cutaneous squamous cell carcinoma (cSCC), we discovered a non-epithelial-mesenchymal transition function for the transcription factor Snail in maintaining the stemness of epidermal keratinocytes. Snail-expressing cells secrete the matricellular protein Mindin, which functions in an autocrine fashion to activate a Src-STAT3 pathway to reinforce their stem/progenitor phenotype. This pathway is activated by the engagement of Mindin with the leukocyte-specific integrin, CD11b (ITGAM), which is also unexpectedly expressed by epidermal keratinocytes. Interestingly, disruption of this signaling module in human cSCC attenuates tumorigenesis, suggesting that targeting Mindin would be a promising therapeutic approach to hinder cancer recurrence.

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.

Regulation of alternative polyadenylation by the C2H2-zinc-finger protein Sp1

Alternative polyadenylation (APA) enhances gene regulatory potential by increasing the diversity of mRNA transcripts. 3' UTR shortening through APA correlates with enhanced cellular proliferation and is a widespread phenomenon in tumor cells. Here, we show that the ubiquitously expressed transcription factor Sp1 binds RNA in vivo and is a common repressor of distal poly(A) site usage. RNA sequencing identified 2,344 genes (36% of the total mapped mRNA transcripts) with lengthened 3' UTRs upon Sp1 depletion. Sp1 preferentially binds the 3' UTRs of such lengthened transcripts and inhibits cleavage at distal sites by interacting with the subunits of the core cleavage and polyadenylation (CPA) machinery. The 3' UTR lengths of Sp1 target genes in breast cancer patient RNA-seq data correlate with Sp1 expression levels, implicating Sp1-mediated APA regulation in modulating tumorigenic properties. Taken together, our findings provide insights into the mechanism for dynamic APA regulation by unraveling a previously unknown function of the DNA-binding transcription factor Sp1.

Predicting which genes will respond to transcription factor perturbations

The ability to predict which genes will respond to the perturbation of a transcription factor serves as a benchmark for our systems-level understanding of transcriptional regulatory networks. In previous work, machine learning models have been trained to predict static gene expression levels in a biological sample by using data from the same or similar samples, including data on their transcription factor binding locations, histone marks, or DNA sequence. We report on a different challenge-training machine learning models to predict which genes will respond to the perturbation of a transcription factor without using any data from the perturbed cells. We find that existing transcription factor location data (ChIP-seq) from human cells have very little detectable utility for predicting which genes will respond to perturbation of a transcription factor. Features of genes, including their preperturbation expression level and expression variation, are very useful for predicting responses to perturbation of any transcription factor. This shows that some genes are poised to respond to transcription factor perturbations and others are resistant, shedding light on why it has been so difficult to predict responses from binding locations. Certain histone marks, including H3K4me1 and H3K4me3, have some predictive power when located downstream of the transcription start site. However, the predictive power of histone marks is much less than that of gene expression level and expression variation. Sequence-based or epigenetic properties of genes strongly influence their tendency to respond to direct transcription factor perturbations, partially explaining the oft-noted difficulty of predicting responsiveness from transcription factor binding location data. These molecular features are largely reflected in and summarized by the gene's expression level and expression variation. Code is available at https://github.com/BrentLab/TFPertRespExplainer.

Toward a base-resolution panorama of the in vivo impact of cytosine methylation on transcription factor binding

Background

While methylation of CpG dinucleotides is traditionally considered antagonistic to the DNA-binding activity of most transcription factors (TFs), recent in vitro studies have revealed a more complex picture, suggesting that over a third of TFs may preferentially bind to methylated sequences. Expanding these in vitro observations to in vivo TF binding preferences is challenging since the effect of methylation of individual CpG sites cannot be easily isolated from the confounding effects of DNA accessibility and regional DNA methylation. Thus, in vivo methylation preferences of most TFs remain uncharacterized.

Results

We introduce joint accessibility-methylation-sequence (JAMS) models, which connect the strength of the binding signal observed in ChIP-seq to the DNA accessibility of the binding site, regional methylation level, DNA sequence, and base-resolution cytosine methylation. We show that JAMS models quantitatively explain TF occupancy, recapitulate cell type-specific TF binding, and have high positive predictive value for identification of TFs affected by intra-motif methylation. Analysis of 2209 ChIP-seq experiments results in high-confidence JAMS models for 260 TFs, revealing a negative association between in vivo TF occupancy and intra-motif methylation for 45% of studied TFs, as well as 16 TFs that are predicted to bind to methylated sites, including 11 novel methyl-binding TFs mostly from the multi-zinc finger family.

Conclusions

Our study substantially expands the repertoire of in vivo methyl-binding TFs, but also suggests that most TFs that prefer methylated CpGs in vitro present themselves as methylation agnostic in vivo, potentially due to the balancing effect of competition with other methyl-binding proteins.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02713-y.

Learning representations of chromatin contacts using a recurrent neural network identifies genomic drivers of conformation

Despite the availability of chromatin conformation capture experiments, discerning the relationship between the 1D genome and 3D conformation remains a challenge, which limits our understanding of their affect on gene expression and disease. We propose Hi-C-LSTM, a method that produces low-dimensional latent representations that summarize intra-chromosomal Hi-C contacts via a recurrent long short-term memory neural network model. We find that these representations contain all the information needed to recreate the observed Hi-C matrix with high accuracy, outperforming existing methods. These representations enable the identification of a variety of conformation-defining genomic elements, including nuclear compartments and conformation-related transcription factors. They furthermore enable in-silico perturbation experiments that measure the influence of cis-regulatory elements on conformation.

Host Zinc-finger CCHC-type containing protein 3 inhibits pseudorabies virus proliferation by regulating type I interferon signaling

Zinc finger CCHC-type containing protein 3 (ZCCHC3) acts as an antiviral factor that interacts with RIG-I and cGAS to modulate innate signaling against viral infections. Here, we investigated the role of porcine ZCCHC3 during pseudorabies virus (PRV) proliferation. We found that porcine ZCCHC3 plays an inhibitory role in the proliferation of PRV by regulating cellular innate immune responses. Further, overexpression of ZCCHC3 inhibited gB protein levels and viral titers, whereas knockdown of ZCCHC3 promoted viral growth. ZCCHC3 overexpression increased IFN-β expression to upregulate downstream gene expression, thus leading to the suppression of viral replication. However, PRV infection reduced the endogenous expression of ZCCHC3 in permissive cells. Importantly, PRV-encoded UL13 and UL24 proteins were identified to inhibit the expression of ZCCHC3, thus antagonizing its antiviral effect. Collectively, our data underscore the important role of ZCCHC3 against PRV infection and promote understandings of viral proteins in PRV pathogenesis.

A role for the mitotic proteins Bub3 and BuGZ in transcriptional regulation of catalase-3 expression

The spindle assembly checkpoint factors Bub3 and BuGZ play critical roles in mitotic process, but little is known about their roles in other cellular processes in eukaryotes. In aerobic organisms, transcriptional regulation of catalase genes in response to developmental or environmental stimuli is necessary for redox homeostasis. Here, we demonstrate that Bub3 and BuGZ negatively regulate cat-3 transcription in the model filamentous fungus Neurospora crassa. The absence of Bub3 caused a significant decrease in BuGZ protein levels. Our data indicate that BuGZ and Bub3 interact directly via the GLEBS domain of BuGZ. Despite loss of the interaction, the amount of BuGZ mutant protein negatively correlated with the cat-3 expression level, indicating that BuGZ amount rather than Bub3-BuGZ interaction determines cat-3 transcription level. Further experiments demonstrated that BuGZ binds directly to the cat-3 gene and responses to cat-3 overexpression induced by oxidative stresses. However, the zinc finger domains of BuGZ have no effects on DNA binding, although mutations of these highly conserved domains lead to loss of cat-3 repression. The deposition of BuGZ along cat-3 chromatin hindered the recruitment of transcription activators GCN4/CPC1 and NC2 complex, thereby preventing the assembly of the transcriptional machinery. Taken together, our results establish a mechanism for how mitotic proteins Bub3 and BuGZ functions in transcriptional regulation in a eukaryotic organism.

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.

Krüppel-like Transcription Factor 7 Is a Causal Gene in Autism Development

Background: Autism spectrum disorder (ASD) is a complex neurodevelopmental disease. To date, more than 1000 genes have been shown to be associated with ASD, and only a few of these genes account for more than 1% of autism cases. Klf7 is an important transcription factor of cell proliferation and differentiation in the nervous system, but whether klf7 is involved in autism is unclear. Methods: We first performed ChIP-seq analysis of klf7 in N2A cells, then performed behavioral tests and RNA-seq in klf7^+/− mice, and finally restored mice with adeno-associated virus (AAV)-mediated overexpression of klf7 in klf7^+/− mice. Results: Klf7 targeted genes are enriched with ASD genes, and 631 ASD risk genes are also differentially expressed in klf7^+/− mice which exhibited the core symptoms of ASD. When klf7 levels were increased in the central nervous system (CNS) in klf7^+/− adult mice, deficits in social interaction, repetitive behavior and majority of dysregulated ASD genes were rescued in the adults, suggesting transcriptional regulation. Moreover, knockdown of klf7 in human brain organoids caused dysregulation of 517 ASD risk genes, 344 of which were shared with klf7^+/− mice, including some high-confidence ASD genes. Conclusions: Our findings highlight a klf7 regulation of ASD genes and provide new insights into the pathogenesis of ASD and promising targets for further research on mechanisms and treatments.

Species-specific KRAB-ZFPs function as repressors of retroviruses by targeting PBS regions

Hosts often target the relatively conserved regions in rapidly mutating retroviruses to inhibit their replication. One of these regions is called a primer binding site (PBS), which has to be complementary to the host tRNA to initiate reverse transcription. By analyzing endogenous retroviral elements, we found that host cells use this sequence as a target in efforts to block the expression of viral elements. A specific type of zinc finger protein targets the PBS in a host genome, which not only inhibits the transcription of endogenous viruses but also inhibits the replication of exogenous retroviruses with the same PBS. Thus, our study sheds light on a strategy for searching for host restriction factors targeting retroviruses.

Eukaryotic genomes harbor sequences derived from the chromosomal integration of ancient viruses, such as endogenous retroviruses (ERVs), which comprise 8% of the human genome. Like exogenous retroviruses, ERVs retain many common functional elements, including the corresponding DNA sequences of transfer RNA (tRNA) primer binding sites (PBSs), which are utilized for reverse transcription initiation by exogenous retroviruses. Here, through a medium-scale analysis of PBS loci positioned within ERVs, coupled with chromatin immunoprecipitation sequencing (ChIP-seq) of Kruppel-associated box zinc finger proteins (KRAB-ZFPs), we identified multiple ZFPs that specifically bind to different PBS loci. Among these, we focused on PBS-Lys, which is utilized by HIV-1, and identified its specific binding proteins to be mouse ZFP961 and human ZNF417/ZNF587. We found that these proteins not only repress ERV transcription but also inhibit retrovirus integration and transcription. Disruption of these ZFPs rendered cells more susceptible to HIV-1 infection. Thus, our research provides a methodology for identifying potential host factors that target retroviruses by ERVs.

ZNF146/OZF and ZNF507 target LINE-1 sequences

Repetitive sequences including transposable elements and transposon-derived fragments account for nearly half of the human genome. While transposition-competent transposable elements must be repressed to maintain genomic stability, mutated and fragmented transposable elements comprising the bulk of repetitive sequences can also contribute to regulation of host gene expression and broader genome organization. Here, we analyzed published ChIP-seq data sets to identify proteins broadly enriched on transposable elements in the human genome. We show 2 of the proteins identified, C2H2 zinc finger-containing proteins ZNF146 (also known as OZF) and ZNF507, are targeted to distinct sites within LINE-1 ORF2 at thousands of locations in the genome. ZNF146 binding sites are found at old and young LINE-1 elements. In contrast, ZNF507 preferentially binds at young LINE-1 sequences correlated to sequence changes in LINE-1 elements at ZNF507's binding site. To gain further insight into ZNF146 and ZNF507 function, we disrupt their expression in HEK293 cells using CRISPR/Cas9 and perform RNA sequencing, finding modest gene expression changes in cells where ZNF507 has been disrupted. We further identify a physical interaction between ZNF507 and PRMT5, suggesting ZNF507 may target arginine methylation activity to LINE-1 sequences.

Screening and Validation of p38 MAPK Involved in Ovarian Development of Brachymystax lenok

Brachymystax lenok (lenok) is a rare cold-water fish native to China that is of high meat quality. Its wild population has declined sharply in recent years, and therefore, exploring the molecular mechanisms underlying the development and reproduction of lenoks for the purposes of artificial breeding and genetic improvement is necessary. The lenok comparative transcriptome was analyzed by combining single molecule, real-time, and next generation sequencing (NGS) technology. Differentially expressed genes (DEGs) were identified in five tissues (head kidney, spleen, liver, muscle, and gonad) between immature [300 days post-hatching (dph)] and mature [three years post-hatching (ph)] lenoks. In total, 234,124 and 229,008 full-length non-chimeric reads were obtained from the immature and mature sequencing data, respectively. After NGS correction, 61,405 and 59,372 non-redundant transcripts were obtained for the expression level and pathway enrichment analyses, respectively. Compared with the mature group, 719 genes with significantly increased expression and 1,727 genes with significantly decreased expression in all five tissues were found in the immature group. Furthermore, DEGs and pathways involved in the endocrine system and gonadal development were identified, and p38 mitogen-activated protein kinases (MAPKs) were identified as potentially regulating gonadal development in lenok. Inhibiting the activity of p38 MAPKs resulted in abnormal levels of gonadotropin-releasing hormone, follicle-stimulating hormone, and estradiol, and affected follicular development. The full-length transcriptome data obtained in this study may provide a valuable reference for the study of gene function, gene expression, and evolutionary relationships in B. lenok and may illustrate the basic regulatory mechanism of ovarian development in teleosts.

SARS-CoV-2 nucleocapsid protein binds host mRNAs and attenuates stress granules to impair host stress response

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) nucleocapsid (N) protein is essential for viral replication, making it a promising target for antiviral drug and vaccine development. SARS-CoV-2 infected patients exhibit an uncoordinated immune response; however, the underlying mechanistic details of this imbalance remain obscure. Here, starting from a functional proteomics workflow, we cataloged the protein-protein interactions of SARS-CoV-2 proteins, including an evolutionarily conserved specific interaction of N with the stress granule resident proteins G3BP1 and G3BP2. N localizes to stress granules and sequesters G3BPs away from their typical interaction partners, thus attenuating stress granule formation. We found that N binds directly to host mRNAs in cells, with a preference for 3' UTRs, and modulates target mRNA stability. We show that the N protein rewires the G3BP1 mRNA-binding profile and suppresses the physiological stress response of host cells, which may explain the imbalanced immune response observed in SARS-CoV-2 infected patients.

The KRAB Domain of ZNF10 Guides the Identification of Specific Amino Acids That Transform the Ancestral KRAB-A-Related Domain Present in Human PRDM9 into a Canonical Modern KRAB-A Domain

Krüppel-associated box (KRAB) zinc finger proteins are a large class of tetrapod transcription factors that usually exert transcriptional repression through recruitment of TRIM28/KAP1. The evolutionary root of modern KRAB domains (mKRAB) can be traced back to an ancestral motif (aKRAB) that occurs even in invertebrates. Here, we first stratified three subgroups of aKRAB sequences from the animal kingdom (PRDM9, SSX and coelacanth KZNF families) and defined ancestral subdomains for KRAB-A and KRAB-B. Using human ZNF10 mKRAB-AB as blueprints for function, we then identified the necessary amino acid changes that transform the inactive aKRAB-A of human PRDM9 into an mKRAB domain capable of mediating silencing and complexing TRIM28/KAP1 in human cells when employed as a hybrid with ZNF10-B. Full gain of function required replacement of residues KR by the conserved motif MLE (positionsA32-A34), which inserted an additional residue, and exchange of A9/S for F, A20/M for L, and A27/R for V. AlphaFold2 modelling documented an evolutionary conserved L-shaped body of two α-helices in all KRAB domains. It is transformed into a characteristic spatial arrangement typical for mKRAB-AB upon the amino acid replacements and in conjunction with a third helix supplied by mKRAB-B. Side-chains pointing outward from the core KRAB 3D structure may reveal a protein-protein interaction code enabling graded binding of TRIM28 to different KRAB domains. Our data provide basic insights into structure-function relationships and emulate transitions of KRAB during evolution.

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.

The Transcription Factor FEZF1, a Direct Target of EWSR1-FLI1 in Ewing Sarcoma Cells, Regulates the Expression of Neural-Specific Genes

Ewing sarcoma is a rare pediatric tumor characterized by chromosomal translocations that give rise to aberrant chimeric transcription factors (e.g., EWSR1-FLI1). EWSR1-FLI1 promotes a specific cellular transcriptional program. Therefore, the study of EWSR1-FLI1 target genes is important to identify critical pathways involved in Ewing sarcoma tumorigenesis. In this work, we focused on the transcription factors regulated by EWSR1-FLI1 in Ewing sarcoma. Transcriptomic analysis of the Ewing sarcoma cell line A673 indicated that one of the genes more strongly upregulated by EWSR1-FLI1 was FEZF1 (FEZ family zinc finger protein 1), a transcriptional repressor involved in neural cell identity. The functional characterization of FEZF1 was performed in three Ewing sarcoma cell lines (A673, SK-N-MC, SK-ES-1) through an shRNA-directed silencing approach. FEZF1 knockdown inhibited clonogenicity and cell proliferation. Finally, the analysis of the FEZF1-dependent expression profile in A673 cells showed several neural genes regulated by FEZF1 and concomitantly regulated by EWSR1-FLI1. In summary, FEZF1 is transcriptionally regulated by EWSR1-FLI1 in Ewing sarcoma cells and is involved in the regulation of neural-specific genes, which could explain the neural-like phenotype observed in several Ewing sarcoma tumors and cell lines.

ZFTA Translocations Constitute Ependymoma Chromatin Remodeling and Transcription Factors

ZFTA (C11orf95)-a gene of unknown function-partners with a variety of transcriptional coactivators in translocations that drive supratentorial ependymoma, a frequently lethal brain tumor. Understanding the function of ZFTA is key to developing therapies that inhibit these fusion proteins. Here, using a combination of transcriptomics, chromatin immunoprecipitation sequencing, and proteomics, we interrogated a series of deletion-mutant genes to identify a tripartite transformation mechanism of ZFTA-containing fusions, including: spontaneous nuclear translocation, extensive chromatin binding, and SWI/SNF, SAGA, and NuA4/Tip60 HAT chromatin modifier complex recruitment. Thereby, ZFTA tethers fusion proteins across the genome, modifying chromatin to an active state and enabling its partner transcriptional coactivators to promote promiscuous expression of a transforming transcriptome. Using mouse models, we validate further those elements of ZFTA-fusion proteins that are critical for transformation-including ZFTA zinc fingers and partner gene transactivation domains-thereby unmasking vulnerabilities for therapeutic targeting. SIGNIFICANCE: Ependymomas are hard-to-treat brain tumors driven by translocations between ZFTA and a variety of transcriptional coactivators. We dissect the transforming mechanism of these fusion proteins and identify protein domains indispensable for tumorigenesis, thereby providing insights into the molecular basis of ependymoma tumorigenesis and vulnerabilities for therapeutic targeting.This article is highlighted in the In This Issue feature, p. 2113.