A comprehensive catalog of human KRAB-associated zinc finger genes: insights into the evolutionary history of a large family of transcriptional repressors

From the cauldron of conflict: Endogenous gene regulation by piRNA and other modes of adaptation enabled by selfish transposable elements

Transposable elements (TEs) provide a prime example of genetic conflict because they can proliferate in genomes and populations even if they harm the host. However, numerous studies have shown that TEs, though typically harmful, can also provide fuel for adaptation. This is because they code functional sequences that can be useful for the host in which they reside. In this review, I summarize the "how" and "why" of adaptation enabled by the genetic conflict between TEs and hosts. In addition, focusing on mechanisms of TE control by small piwi-interacting RNAs (piRNAs), I highlight an indirect form of adaptation enabled by conflict. In this case, mechanisms of host defense that regulate TEs have been redeployed for endogenous gene regulation. I propose that the genetic conflict released by meiosis in early eukaryotes may have been important because, among other reasons, it spurred evolutionary innovation on multiple interwoven trajectories - on the part of hosts and also embedded genetic parasites. This form of evolution may function as a complexity generating engine that was a critical player in eukaryotic evolution.

Cellular zinc status alters chromatin accessibility and binding of p53 to DNA

Zn2+ is an essential metal required by approximately 850 human transcription factors. How these proteins acquire their essential Zn2+ cofactor and whether they are sensitive to changes in the labile Zn2+ pool in cells remain open questions. Using ATAC-seq to profile regions of accessible chromatin coupled with transcription factor enrichment analysis, we examined how increases and decreases in the labile zinc pool affect chromatin accessibility and transcription factor enrichment. We found 685 transcription factor motifs were differentially enriched, corresponding to 507 unique transcription factors. The pattern of perturbation and the types of transcription factors were notably different at promoters versus intergenic regions, with zinc-finger transcription factors strongly enriched in intergenic regions in elevated Zn2+ To test whether ATAC-seq and transcription factor enrichment analysis predictions correlate with changes in transcription factor binding, we used ChIP-qPCR to profile six p53 binding sites. We found that for five of the six targets, p53 binding correlates with the local accessibility determined by ATAC-seq. These results demonstrate that changes in labile zinc alter chromatin accessibility and transcription factor binding to DNA.

DiPRO1 distinctly reprograms muscle and mesenchymal cancer cells

We have recently identified the uncharacterized ZNF555 protein as a component of a productive complex involved in the morbid function of the 4qA locus in facioscapulohumeral dystrophy. Subsequently named DiPRO1 (Death, Differentiation, and PROliferation related PROtein 1), our study provides substantial evidence of its role in the differentiation and proliferation of human myoblasts. DiPRO1 operates through the regulatory binding regions of SIX1, a master regulator of myogenesis. Its relevance extends to mesenchymal tumors, such as rhabdomyosarcoma (RMS) and Ewing sarcoma, where DiPRO1 acts as a repressor via the epigenetic regulators TIF1B and UHRF1, maintaining methylation of cis-regulatory elements and gene promoters. Loss of DiPRO1 mimics the host defense response to virus, awakening retrotransposable repeats and the ZNF/KZFP gene family. This enables the eradication of cancer cells, reprogramming the cellular decision balance towards inflammation and/or apoptosis by controlling TNF-α via NF-kappaB signaling. Finally, our results highlight the vulnerability of mesenchymal cancer tumors to si/shDiPRO1-based nanomedicines, positioning DiPRO1 as a potential therapeutic target.

Highly efficient and specific regulation of gene expression using enhanced CRISPR-Cas12f system

The recently developed CRISPR activator (CRISPRa) system uses a CRISPR-Cas effector-based transcriptional activator to effectively control the expression of target genes without causing DNA damage. However, existing CRISPRa systems based on Cas9/Cas12a necessitate improvement in terms of efficacy and accuracy due to limitations associated with the CRISPR-Cas module itself. To overcome these limitations and effectively and accurately regulate gene expression, we developed an efficient CRISPRa system based on the small CRISPR-Cas effector Candidatus Woesearchaeota Cas12f (CWCas12f). By engineering the CRISPR-Cas module, linking activation domains, and using various combinations of linkers and nuclear localization signal sequences, the optimized eCWCas12f-VPR system enabled effective and target-specific regulation of gene expression compared with that using the existing CRISPRa system. The eCWCas12f-VPR system developed in this study has substantial potential for controlling the transcription of endogenous genes in living organisms and serves as a foundation for future gene therapy and biological research.

A novel gain-of-function phosphorylation site modulates PTPN22 inhibition of TCR signaling

Protein tyrosine phosphatase nonreceptor type 22 (PTPN22) is encoded by a major autoimmunity gene and is a known inhibitor of T cell receptor (TCR) signaling and drug target for cancer immunotherapy. However, little is known about PTPN22 posttranslational regulation. Here, we characterize a phosphorylation site at Ser325 situated C terminal to the catalytic domain of PTPN22 and its roles in altering protein function. In human T cells, Ser325 is phosphorylated by glycogen synthase kinase-3 (GSK3) following TCR stimulation, which promotes its TCR-inhibitory activity. Signaling through the major TCR-dependent pathway under PTPN22 control was enhanced by CRISPR/Cas9-mediated suppression of Ser325 phosphorylation and inhibited by mimicking it via glutamic acid substitution. Global phospho-mass spectrometry showed Ser325 phosphorylation state alters downstream transcriptional activity through enrichment of Swi3p, Rsc8p, and Moira domain binding proteins, and next-generation sequencing revealed it differentially regulates the expression of chemokines and T cell activation pathways. Moreover, in vitro kinetic data suggest the modulation of activity depends on a cellular context. Finally, we begin to address the structural and mechanistic basis for the influence of Ser325 phosphorylation on the protein's properties by deuterium exchange mass spectrometry and NMR spectroscopy. In conclusion, this study explores the function of a novel phosphorylation site of PTPN22 that is involved in complex regulation of TCR signaling and provides details that might inform the future development of allosteric modulators of PTPN22.

Transposable elements regulate thymus development and function

Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/β. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

Zinc finger Protein207 orchestrates glioma migration through regulation of epithelial-mesenchymal transition

BACKGROUND

Glioma represents the predominant primary malignant brain tumor. For several years, molecular profiling has been instrumental in the management and therapeutic stratification of glioma, providing a deeper understanding of its biological complexity. Accumulating evidence unveils the putative involvement of zinc finger proteins (ZNFs) in cancer. This study aimed to elucidate the role and significance of ZNF207 in glioma.

METHODS

Utilizing online data such as The Cancer Genome Atlas (TCGA), the Chinese Glioma Genome Atlas (CGGA), the Genotype-Tissue Expression (GTEx) project, the Clinical Proteomic Tumor Analysis Consortium (CPTAC), and the Human Protein Atlas (HPA) databases, in conjunction with bioinformatics methodologies including GO, KEGG, GSEA, CIBERSORT immune cell infiltration estimation, and protein-protein interaction (PPI) analysis, enabled a comprehensive exploration of ZNF207's involvement in gliomagenesis. Immunohistochemistry and RT-PCR techniques were employed to validate the expression level of ZNF207 in glioma samples. Subsequently, the biological effects of ZNF207 on glioma cells were explored through in vitro assays.

RESULTS

Our results demonstrate elevated expression of ZNF207 in gliomas, correlating with unfavorable patient outcomes. Stratification analyses were used to delineate the prognostic efficacy of ZNF207 in glioma with different clinicopathological characteristics. Immunocorrelation analysis revealed a significant association between ZNF207 expression and the infiltration levels of T helper cells, macrophages, and natural killer (NK) cells. Utilizing ZNF207 expression and clinical features, we constructed an OS prediction model and displayed well discrimination with a C-index of 0.861. Moreover, the strategic silencing of ZNF207 attenuated glioma cell advancement, evidenced by diminished cellular proliferation, weakened cell tumorigenesis, augmented apoptotic activity, and curtailed migratory capacity alongside the inhibition of the epithelial-mesenchymal transition (EMT) pathway.

CONCLUSIONS

ZNF207 may identify as a prospective biomarker and therapeutic candidate for glioma prevention, providing valuable insights into understanding glioma pathogenesis and treatment strategies.

Genomic profiling of six human somatic histone H1 variants denotes that H1X accumulates at recently incorporated transposable elements

Histone H1, a vital component in chromatin structure, binds to linker DNA and regulates nuclear processes. We have investigated the distribution of histone H1 variants in a breast cancer cell line using ChIP-Seq. Two major groups of variants are identified: H1.2, H1.3, H1.5 and H1.0 are abundant in low GC regions (B compartment), while H1.4 and H1X preferentially localize in high GC regions (A compartment). Examining their abundance within transposable elements (TEs) reveals that H1X and H1.4 are enriched in recently-incorporated TEs (SVA and SINE-Alu), while H1.0/H1.2/H1.3/H1.5 are more abundant in older elements. Notably, H1X is particularly enriched in SVA families, while H1.4 shows the highest abundance in young AluY elements. Although low GC variants are generally enriched in LINE, LTR and DNA repeats, H1X and H1.4 are also abundant in a subset of recent LINE-L1 and LTR repeats. H1X enrichment at SVA and Alu is consistent across multiple cell lines. Further, H1X depletion leads to TE derepression, suggesting its role in maintaining TE repression. Overall, this study provides novel insights into the differential distribution of histone H1 variants among repetitive elements, highlighting the potential involvement of H1X in repressing TEs recently incorporated within the human genome.

Copy number variation of the ZNF679 gene in cattle and its association analysis with growth traits

Copy number variation (CNV) is an important member of genetic structural variation that exists widely in animal genomes and is between 50 bp and several Mb in length and widely used in research's of animal genetics and breeding. ZNF679 is an important transcription factor, which has been found association with diseases in the human genome many times. This gene has also been found to be associated with cattle growth traits in previous re-sequencing studies. We tested the CNVs of the ZNF679 gene in 809 individuals from 7 Chinese cattle breeds and tested the association between the CNVs and growth traits in 552 individuals from 5 breeds. The results demonstrated the correlation the correlation between the CNVs of the ZNF679 gene and some Chinese cattle (QC cattle and XN cattle) growth traits. To sum up, this study indicated that ZNF679-CNVs can be used as a candidate gene for molecular genetic marker-assisted selection breeding for cattle growth traits to contribute to the development of genetic improvement of Chinese cattle.

PhylteR: Efficient Identification of Outlier Sequences in Phylogenomic Datasets

In phylogenomics, incongruences between gene trees, resulting from both artifactual and biological reasons, can decrease the signal-to-noise ratio and complicate species tree inference. The amount of data handled today in classical phylogenomic analyses precludes manual error detection and removal. However, a simple and efficient way to automate the identification of outliers from a collection of gene trees is still missing. Here, we present PhylteR, a method that allows rapid and accurate detection of outlier sequences in phylogenomic datasets, i.e. species from individual gene trees that do not follow the general trend. PhylteR relies on DISTATIS, an extension of multidimensional scaling to 3 dimensions to compare multiple distance matrices at once. In PhylteR, these distance matrices extracted from individual gene phylogenies represent evolutionary distances between species according to each gene. On simulated datasets, we show that PhylteR identifies outliers with more sensitivity and precision than a comparable existing method. We also show that PhylteR is not sensitive to ILS-induced incongruences, which is a desirable feature. On a biological dataset of 14,463 genes for 53 species previously assembled for Carnivora phylogenomics, we show (i) that PhylteR identifies as outliers sequences that can be considered as such by other means, and (ii) that the removal of these sequences improves the concordance between the gene trees and the species tree. Thanks to the generation of numerous graphical outputs, PhylteR also allows for the rapid and easy visual characterization of the dataset at hand, thus aiding in the precise identification of errors. PhylteR is distributed as an R package on CRAN and as containerized versions (docker and singularity).

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.

ZNF714 Supports Pro-Oncogenic Features in Lung Cancer Cells

Despite the ongoing progress in diagnosis and treatments, cancer remains a threat to more than one-third of the human population. The emerging data indicate that many Krüppel-associated box zinc finger proteins (KRAB-ZNF) belonging to a large gene family may be involved in carcinogenesis. Our previous study identified Zinc Finger Protein 714 (ZNF714), a KRAB-ZNF gene of unknown function, as being commonly overexpressed in many tumors, pointing to its hypothetical oncogenic role. Here, we harnessed The Cancer Genome Atlas (TCGA)-centered databases and performed functional studies with transcriptomic and methylomic profiling to explore ZNF714 function in cancer. Our pan-cancer analyses confirmed frequent ZNF714 overexpression in multiple tumors, possibly due to regional amplification, promoter hypomethylation, and Nuclear Transcription Factor Y Subunit Beta (NFYB) signaling. We also showed that ZNF714 expression correlates with tumor immunosuppressive features. The in vitro studies indicated that ZNF714 expression positively associates with proliferation, migration, and invasion. The transcriptomic analysis of ZNF714 knocked-down cells demonstrated deregulation of cell adhesion, migration, proliferation, apoptosis, and differentiation. Importantly, we provided evidence that ZNF714 negatively regulates the expression of several known TSGs indirectly via promoter methylation. However, as ZNF714 did not show nuclear localization in our research model, the regulatory mechanisms exerted by ZNF714 require further investigation. In conclusion, our results reveal, for the first time, that ZNF714 may support pro-oncogenic features in lung cancer 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.

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.

Transposable Elements: Epigenetic Silencing Mechanisms or Modulating Tools for Vertebrate Adaptations? Two Sides of the Same Coin

Transposable elements constitute one of the main components of eukaryotic genomes. In vertebrates, they differ in content, typology, and family diversity and played a crucial role in the evolution of this taxon. However, due to their transposition ability, TEs can be responsible for genome instability, and thus silencing mechanisms were evolved to allow the coexistence between TEs and eukaryotic host-coding genes. Several papers are highlighting in TEs the presence of regulatory elements involved in regulating nearby genes in a tissue-specific fashion. This suggests that TEs are not sequences merely to silence; rather, they can be domesticated for the regulation of host-coding gene expression, permitting species adaptation and resilience as well as ensuring human health. This review presents the main silencing mechanisms acting in vertebrates and the importance of exploiting these mechanisms for TE control to rewire gene expression networks, challenging the general view of TEs as threatening elements.

Systematic characterization of regulatory variants of blood pressure genes

High blood pressure (BP) is the major risk factor for cardiovascular disease. Genome-wide association studies have identified genetic variants for BP, but functional insights into causality and related molecular mechanisms lag behind. We functionally characterize 4,608 genetic variants in linkage with 135 BP loci in vascular smooth muscle cells and cardiomyocytes by massively parallel reporter assays. High densities of regulatory variants at BP loci (i.e., ULK4, MAP4, CFDP1, PDE5A) indicate that multiple variants drive genetic association. Regulatory variants are enriched in repeats, alter cardiovascular-related transcription factor motifs, and spatially converge with genes controlling specific cardiovascular pathways. Using heuristic scoring, we define likely causal variants, and CRISPR prime editing finally determines causal variants for KCNK9, SFXN2, and PCGF6, which are candidates for developing high BP. Our systems-level approach provides a catalog of functionally relevant variants and their genomic architecture in two trait-relevant cell lines for a better understanding of BP gene regulation.

ZFP57 promotes ovarian cancer progression by transcriptionally regulating BRCA1 and managing G1 checkpoint

Ovarian cancer (OC) which is one of the frequently-occurring gynecologic malignant tumors, endangers the health of women. The zinc finger protein 57 (ZFP57) plays crucial functions during the progression of cancer and is reported as a prognostic and therapeutic candidate in a variety of cancer. However, the biological function as well as the underlying mechanism of ZFP57 during OC progression remains unknown. Here, ZFP57 expression was found prominently increased in OC tissues and correlated with the prognosis of OC patients. Knock down of ZFP57 in OC cells inhibited the cell proliferation and migration, and also arrested the cells at G1 phase as well as accelerated the apoptosis. Additionally, ZFP57 transcriptionally regulated BRCA1 expression in OC, indicating that ZFP57 may affect BRCA1 mediated G1 checkpoint to regulate the cell cycle of OC cells and further influence the progression of OC. Taken together, our present study discovered a novel function of ZFP57 in OC, suggesting that ZFP57 could be potentially treated as a prognostic biomarker and therapeutic target for OC patients.

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.

Acetylation-Mimic Mutation of TRIM28-Lys304 to Gln Attenuates the Interaction with KRAB-Zinc-Finger Proteins and Affects Gene Expression in Leukemic K562 Cells

TRIM28/KAP1/TIF1β is a crucial epigenetic modifier. Genetic ablation of trim28 is embryonic lethal, although RNAi-mediated knockdown in somatic cells yields viable cells. Reduction in TRIM28 abundance at the cellular or organismal level results in polyphenism. Posttranslational modifications such as phosphorylation and sumoylation have been shown to regulate TRIM28 activity. Moreover, several lysine residues of TRIM28 are subject to acetylation, but how acetylation of TRIM28 affects its functions remains poorly understood. Here, we report that, compared with wild-type TRIM28, the acetylation-mimic mutant TRIM28-K304Q has an altered interaction with Krüppel-associated box zinc-finger proteins (KRAB-ZNFs). The TRIM28-K304Q knock-in cells were created in K562 erythroleukemia cells by CRISPR-Cas9 (Clustered regularly interspaced short palindromic repeats/CRISPR-associated protein nuclease 9) gene editing method. Transcriptome analysis revealed that TRIM28-K304Q and TRIM28 knockout K562 cells had similar global gene expression profiles, yet the profiles differed considerably from wild-type K562 cells. The expression levels of embryonic-related globin gene and a platelet cell marker integrin-beta 3 were increased in TRIM28-K304Q mutant cells, indicating the induction of differentiation. In addition to the differentiation-related genes, many zinc-finger-proteins genes and imprinting genes were activated in TRIM28-K304Q cells; they were inhibited by wild-type TRIM28 via binding with KRAB-ZNFs. These results suggest that acetylation/deacetylation of K304 in TRIM28 constitutes a switch for regulating its interaction with KRAB-ZNFs and alters the gene regulation as demonstrated by the acetylation mimic TRIM28-K304Q.

A Dual-Function "TRE-Lox" System for Genetic Deletion or Reversible, Titratable, and Near-Complete Downregulation of Cathepsin D

Commonly employed methods for reversibly disrupting gene expression, such as those based on RNAi or CRISPRi, are rarely capable of achieving >80-90% downregulation, making them unsuitable for targeting genes that require more complete disruption to elicit a phenotype. Genetic deletion, on the other hand, while enabling complete disruption of target genes, often produces undesirable irreversible consequences such as cytotoxicity or cell death. Here we describe the design, development, and detailed characterization of a dual-function "TRE-Lox" system for effecting either (a) doxycycline (Dox)-mediated downregulation or (b) genetic deletion of a target gene-the lysosomal aspartyl protease cathepsin D (CatD)-based on targeted insertion of a tetracycline-response element (TRE) and two LoxP sites into the 5' end of the endogenous CatD gene (CTSD). Using an optimized reverse-tetracycline transrepressor (rtTR) variant fused with the Krüppel-associated box (KRAB) domain, we show that CatD expression can be disrupted by as much as 98% in mouse embryonic fibroblasts (MEFs). This system is highly sensitive to Dox (IC₅₀ = 1.46 ng/mL) and results in rapid (t_1/2 = 0.57 d) and titratable downregulation of CatD. Notably, even near-total disruption of CatD expression was completely reversed by withdrawal of Dox. As expected, transient expression of Cre recombinase results in complete deletion of the CTSD gene. The dual functionality of this novel system will facilitate future studies of the involvement of CatD in various diseases, particularly those attributable to partial loss of CatD function. In addition, the TRE-Lox approach should be applicable to the regulation of other target genes requiring more complete disruption than can be achieved by traditional methods.

Classification of domains in predicted structures of the human proteome

Recent advances in protein structure prediction have generated accurate structures of previously uncharacterized human proteins. Identifying domains in these predicted structures and classifying them into an evolutionary hierarchy can reveal biological insights. Here, we describe the detection and classification of domains from the human proteome. Our classification indicates that only 62% of residues are located in globular domains. We further classify these globular domains and observe that the majority (65%) can be classified among known folds by sequence, with a smaller fraction (33%) requiring structural data to refine the domain boundaries and/or to support their homology. A relatively small number (966 domains) cannot be confidently assigned using our automatic pipelines, thus demanding manual inspection. We classify 47,576 domains, of which only 23% have been included in experimental structures. A portion (6.3%) of these classified globular domains lack sequence-based annotation in InterPro. A quarter (23%) have not been structurally modeled by homology, and they contain 2,540 known disease-causing single amino acid variations whose pathogenesis can now be inferred using AF models. A comparison of classified domains from a series of model organisms revealed expansions of several immune response-related domains in humans and a depletion of olfactory receptors. Finally, we use this classification to expand well-known protein families of biological significance. These classifications are presented on the ECOD website (http://prodata.swmed.edu/ecod/index_human.php).

DNA methylation-mediated low expression of ZNF582 promotes the proliferation, migration, and invasion of clear cell renal cell carcinoma

OBJECTIVE

The methylation of DNA promoter region mediates the low expression of many tumor suppressor genes and plays an essential part in cancer progression. We investigated methylation and expression of ZNF582 in clear cell renal cell carcinoma (ccRCC), and to study the function of ZNF582 in ccRCC cells.

METHODS

Methylation data and mRNA expression data of TCGA-KIRC were obtained from TCGA database to screen methylation-driven genes. Survival analysis and gene set enrichment analysis (GSEA) were done for the target gene. The methylation degree and mRNA level of ZNF582 in ccRCC cell line were detected by methylation-specific PCR (MSP) and qRT-PCR, respectively. Effects of overexpression of ZNF582 on ccRCC cells were assessed via CCK-8, flow cytometry, wound healing, Transwell, and cell adhesion assays.

RESULTS

Eighteen methylation-driven genes were identified via bioinformatics methods. Among them, ZNF582 was noticeably hypermethylated and lowly expressed in tumor tissue, and ZNF582 methylation and expression levels were pronouncedly associated with prognosis and clinical stage. MSP also displayed that the ZNF582 DNA promoter region was hypermethylated in ccRCC cells, and the mRNA expression of ZNF582 was dramatically elevated after demethylation. In vitro cell experiments disclosed that overexpression of ZNF582 markedly hindered cell proliferation, invasion, migration, and fostered cell apoptosis and adhesion of ccRCC.

CONCLUSION

ZNF582 was hypermethylated in ccRCC, which mediated its low level. Overexpression of ZNF582 inhibited tumor cell proliferation, migration and invasion. This study generates novel ideas for ccRCC diagnosis and treatment.

Poly zinc finger protein ZFP14 suppresses lymphomagenesis and abnormal inflammatory response via the HOXA gene cluster

Poly zinc finger proteins (ZFP) that contain a KRAB (Krüppel-associated box) domain represent the largest class of transcription factors in higher organisms, but their roles in development and pathogenesis are largely undefined. ZFP14 (also known as ZNF531) contains thirteen zinc fingers and is highly conserved across species. Notably, we found that ZFP14 is frequently down-regulated in a multitude of human cancers, which correlates with poor prognosis of patients. Since ZFP14 has never been characterized, we generated a Zfp14-deficient mouse model to investigate the role of ZFP14 in development and pathogenesis. We showed that the mice deficient in Zfp14 had a short lifespan and were prone to diffuse large B-cell lymphoma (DLBCL), hyperplasia in multiple organs, systemic chronic inflammation, liver steatosis, and pancreatitis. Additionally, several pro-inflammatory cytokines, including IL-1β, IL18, and TNFα, were highly expressed in inflamed Zfp14-/- mice spleens, livers, kidneys and lungs. To determine the underlying mechanism, RNA-seq analysis was performed and showed that the loss of ZFP14 led to increased expression of inflammatory and tumor-promoting genes. Out of the various tumor-promoting genes upregulated by ZFP14 loss, the HOXA gene cluster, which is known to promote lymphomagenesis and conserved between mouse and human, is consistently induced by loss of ZFP14. Moreover, we showed that the HOXA gene expression was inversely correlated with that of ZFP14 in human cancer patients and higher HOXA1 expression was correlated with poor patient prognosis. Together, we postulate that ZFP14 suppresses lymphomagenesis and abnormal inflammatory response by maintaining proper expression of the HOXA gene cluster.

KRAB family is involved in network shifts in response to osmotic stress in camels

A feature of the camel is its tolerance to osmotic stress. However, few studies of osmotic stress in vivo or comparative analyses between different tissues of the camel have been performed. Here, we report the roles of Krüppel-associated box domain containing zinc-finger repressor proteins (KRAB-ZFPs) in transcriptional networks under osmotic stress in camels by analyzing transcriptomes of four different tissues under various osmotic conditions. We found that 273 of 278 KRAB-ZFPs were expressed in our data set, being involved in all of the 65 identified networks and exhibiting their extensive functional diversity. We also found that 110 KRAB-ZFPs were hub genes involved in more than half of the networks. We demonstrated that the osmotic stress response is involved in network shifts and that KRAB-ZFPs mediate this process. Finally, we presented the diverse mechanisms of osmotic stress responses in different tissues. These results revealed the genetic architecture of systematic physiological response in vivo to osmotic stress in camels. Our work will lead to new directions for studying the mechanism of osmotic stress response in anti-arid mammals.

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

Background

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

Results

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

Conclusions

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

Widespread expression of the ancient HERV-K (HML-2) provirus group in normal human tissues

Human endogenous retrovirus (HERV) transcripts are known to be highly expressed in cancers, yet their activity in nondiseased tissue is largely unknown. Using the GTEx RNA-seq dataset from normal tissue sampled at autopsy, we characterized individual expression of the recent HERV-K (HML-2) provirus group across 13,000 different samples of 54 different tissues from 948 individuals. HML-2 transcripts could be identified in every tissue sampled and were elevated in the cerebellum, pituitary, testis, and thyroid. A total of 37 different individual proviruses were expressed in 1 or more tissues, representing all 3 LTR5 subgroups. Nine proviruses were identified as having long terminal repeat (LTR)-driven transcription, 7 of which belonged to the most recent LTR5HS subgroup. Proviruses of different subgroups displayed a bias in tissue expression, which may be associated with differences in transcription factor binding sites in their LTRs. Provirus expression was greater in evolutionarily older proviruses with an earliest shared ancestor of gorilla or older. HML-2 expression was significantly affected by biological sex in 1 tissue, while age and timing of death (Hardy score) had little effect. Proviruses containing intact gag, pro, and env open reading frames (ORFs) were expressed in the dataset, with almost every tissue measured potentially expressing at least 1 intact ORF (gag).

Human endogenous retrovirus (HERV) transcripts are known to be highly expressed in cancers, but what is their activity in normal tissue? This study uses unique patterns of HERV-K RNA expression in the large GEx dataset from non-diseased tissue sites to provide new insights into both the coevolution of HERV-K with our primate ancestors and their current role in human biology.

Mapping the interaction between Trim28 and the KRAB domain at the center of Trim28 silencing of endogenous retroviruses

Transcription of endogenous retroviral elements are tightly regulated during development by members of the KRAB‐containing zinc finger proteins (KRAB‐ZFPs) and the co‐repressor Trim28 (also known as Kap‐1 or Tif1β). KRAB‐ZFPs form the largest family of transcription regulators in mammals and initiate transcriptional silencing by tethering Trim28 to a target locus. Subsequently, Trim28 recruits chromatin modifying effectors resulting in the formation of heterochromatin. In the present study, we identify surface exposed residues on the central six turns of the Trim28 coiled‐coil region forming the binding interface for the KRAB domain. Using AlphaFold2 (AF2) we provide high confidence models of the interface between Trim28 and the KRAB domain and identified leucine 301 on each chain of the Trim28 monomer to act as a pin extending into a hydrophobic pocket on the KRAB domain surface. Site directed mutations in the Trim28‐KRAB binding interface abolished binding to the KRAB domain. Our work provides a detailed understanding of the specific interactions between the KRAB domain and the Trim28 coiled‐coil and how this interaction may be regulated during silencing events.

A resource for integrated genomic analysis of the human liver

In this study, we generated whole-transcriptome RNA-Seq from n = 192 genotyped liver samples and used these data with existing data from the GTEx Project (RNA-Seq) and previous liver eQTL (microarray) studies to create an enhanced transcriptomic sequence resource in the human liver. Analyses of genotype-expression associations show pronounced enrichment of associations with genes of drug response. The associations are primarily consistent across the two RNA-Seq datasets, with some modest variation, indicating the importance of obtaining multiple datasets to produce a robust resource. We further used an empirical Bayesian model to compare eQTL patterns in liver and an additional 20 GTEx tissues, finding that MHC genes, and especially class II genes, are enriched for liver-specific eQTL patterns. To illustrate the utility of the resource to augment GWAS analysis with small sample sizes, we developed a novel meta-analysis technique to combine several liver eQTL data sources. We also illustrate its application using a transcriptome-enhanced re-analysis of a study of neutropenia in pancreatic cancer patients. The associations of genotype with liver expression, including splice variation and its genetic associations, are made available in a searchable genome browser.

Targeted demethylation at ZNF154 promotor upregulates ZNF154 expression and inhibits the proliferation and migration of Esophageal Squamous Carcinoma cells

Zinc finger protein 154 (ZNF154) is hypermethylated at the promoter in many epithelial-derived solid tumors. However, its methylation status and function in esophageal squamous carcinoma (ESCC) are poorly understood. We found that the ZNF154 promoter is hypermethylated in ESCC and portends poor prognosis. In addition, ZNF154 functions as a tumor suppressor gene (TSG) in ESCC, and is downregulated by promoter hypermethylation. We established a targeted demethylation strategy based on CRISPR/dCas9 technology and found that the hypermethylation of ZNF154 promoter repressed ZNF154 induction, which in turn promoted the proliferation and migration of ESCC cells in vitro and in vivo. Finally, high-throughput CUT&Tag analysis, GEPIA software and qPCR were used to revealed the role of ZNF154 as a transcription factor to upregulate the expression of ESCC-associated tumor suppressor genes. Taken together, hypermethylation of the ZNF154 promoter plays an important role in the development of ESCC, and epigenetic editing is a promising tool for inhibiting ESCC cells with aberrant DNA methylation.

19q13.12 KRAB zinc-finger protein ZNF383 represses p53 signaling pathway by interacting with p53

As one of the most important tumor suppressors, the activity of p53 is precisely regulated. However, the mechanism of p53 regulation is still being elucidated and new regulatory molecules for p53 have also been frequently identified. Our previous works revealed that two members of the KRAB zinc-finger protein (KZFP) family Apak and PISA, which are located on human 19q13.12, participated in the regulation of p53 signaling pathway. KZFPs genes are mainly amplified via tandem in situ duplication during evolution, which indicates that similar sequences and functions may be conserved in evolutionarily and physically close KZFPs. Here, we revealed that ZNF383, another member of the KZFPs mapped at 19q13.12, could inhibit p53-mediated apoptosis and the activation of IFN-β pathway by decreasing the H3K36me2 level at p53's binding sites and the attenuating the binding of p53 to its target genes. We further explored the effect of other KZFPs clustered on 19q13.12 on p53, and found that 85% of these KZFPs exerted p53-repressive activity. Intriguingly, an acidic amino acid-enriched sequence, the SAcL motif in the zinc-finger domains of these KZFPs, was found to be critical for p53 binding. Taken together, our findings revealed the KZFPs cluster located at 19q13.12 not only was involved in p53 regulation but also exhibited different features in the selective regulation of p53 and functional mechanisms, and for the first time, confirmed a motif in KZFPs that mediates the interaction of KZFPs and p53. These results would enrich the knowledge on the role, sequence features, and functional mechanisms of the KZFP family in p53 regulation.

Along the Bos taurus genome, uncover candidate imprinting control regions

Background

In mammals, Imprinting Control Regions (ICRs) regulate a subset of genes in a parent-of-origin-specific manner. In both human and mouse, previous studies identified a set of CpG-rich motifs occurring as clusters in ICRs and germline Differentially Methylated Regions (gDMRs). These motifs consist of the ZFP57 binding site (ZFBS) overlapping a subset of MLL binding units known as MLL morphemes. MLL or MLL1 (Mixed Lineage Leukemia 1) is a relatively large multidomain protein that plays a central role in the regulation of transcription. The structures of both MLL1 and MLL2 include a domain (MT) that binds CpG-rich DNA and a conserved domain (SET) that methylates lysine 4 in histone H3 producing H3K4me3 marks in chromatin.

Results

Since genomic imprinting impacts many developmental and key physiological processes, we followed a previous bioinformatics strategy to pinpoint ICR positions in the Bos taurus genome. Initial genome-wide analyses involved finding the positions of ZFP57 binding sites, and the CpG-rich motifs (ZFBS-morph overlaps) along cattle chromosomal DNA. By creating plots displaying the density of ZFBS-morph overlaps, we removed background noise and thus improved signal detection. With the density-plots, we could view the positions of peaks locating known and candidate ICRs in cattle DNA. Our evaluations revealed the correspondence of peaks in plots to reported known and inferred ICRs/DMRs in cattle. Beside peaks pinpointing such ICRs, the density-plots also revealed additional peaks. Since evaluations validated the robustness of our approach, we inferred that the additional peaks may correspond to candidate ICRs for imprinted gene expression.

Conclusion

Our bioinformatics strategy offers the first genome-wide approach for systematically localizing candidate ICRs. Furthermore, we have tailored our datasets for upload onto the UCSC genome browser so that researchers could find known and candidate ICRs with respect to a wide variety of annotations at all scales: from the positions of Single Nucleotide Polymorphisms (SNPs), to positions of genes, transcripts, and repeated DNA elements. Furthermore, the UCSC genome browser offers tools to produce enlarged views: to uncover the genes in the vicinity of candidate ICRs and thus discover potential imprinted genes for experimental validations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08694-3.

ZNF582 promoter methylation predicts cervical cancer radiosensitivity and ZNF582 protein overexpression reduces radiosensitivity by cell cycle arrest in S phase

This study aimed to investigate the relationship between ZNF582 promoter methylation (ZNF582^m) level and radiosensitivity of cervical cancer and its biological basis. This was a prospective multicenter clinical study, comprised of two independent cohorts of locally advanced cervical cancer patients. Exfoliated cervical cells were collected at 0, 24, 30, 36, 48, and 64 Gy to test ZNF582^m levels. Radiotherapy response evaluated according to RECIST Version. RT-PCR, WT were used to detect the mRNA and protein expression levels; MTT, flow cytometry were used to detect the cell viability and cell cycle, respectively. While clone formation and subcutaneous tumorigenesis in nude mice were used to detect the growth of HeLa cells with/without ZNF582 overexpression. In the first cohort, 22 cases achieved complete remission (CR) or partial response (PR), and the other 28 cases exhibited stable disease (SD). Radiotherapy reduced ZNF582^m levels among all patients. Initial lever of ZNF582^m was significantly higher in the Responder (CR + PR) group than in SD group. Also, patients with higher initial lever ZNF582^m were more sensitive toward radiotherapy than ZNF582^m-low patients. The second cohort confirmed above results. The amplitude of ZNF582^m levels were related to radiotherapeutic response; some patients of ZNF582^m-low showed a transient increase in ZNF582^m, and present greater radiosensitivity than other ZNF582^m-low patients. In vitro, ZNF582 protein overexpression promoted cell cycle arrest in S phase. These results suggested that higher ZNF582^m levels predicted greater radiosensitivity in clinical cervical cancer cases. Overexpressed ZNF582 conferred radioresistance by cell cycle arrest in vitro.

Biallelic variants in ZNF142 lead to a syndromic neurodevelopmental disorder

Biallelic variants of the gene encoding for the zinc-finger protein 142 (ZNF142) have recently been associated with intellectual disability (ID), speech impairment, seizures, and movement disorders in nine individuals from five families. In this study, we obtained phenotype and genotype information of 26 further individuals from 16 families. Among the 27 different ZNF142 variants identified in the total of 35 individuals only four were missense. Missense variants may give a milder phenotype by changing the local structure of ZF motifs as suggested by protein modelling; but this correlation should be validated in larger cohorts and pathogenicity of the missense variants should be investigated with functional studies. Clinical features of the 35 individuals suggest that biallelic ZNF142 variants lead to a syndromic neurodevelopmental disorder with mild to moderate ID, varying degrees of delay in language and gross motor development, early onset seizures, hypotonia, behavioral features, movement disorders, and facial dysmorphism. The differences in symptom frequencies observed in the unpublished individuals compared to those of published, and recognition of previously underemphasized facial features are likely to be due to the small sizes of the previous cohorts, which underlines the importance of larger cohorts for the phenotype descriptions of rare genetic disorders. This article is protected by copyright. All rights reserved.

Establishment of H3K9-methylated heterochromatin and its functions in tissue differentiation and maintenance

Heterochromatin is characterized by dimethylated or trimethylated histone H3 Lys9 (H3K9me2 or H3K9me3, respectively) and is found at transposable elements, satellite repeats and genes, where it ensures their transcriptional silencing. The histone methyltransferases (HMTs) that methylate H3K9 — in mammals Suppressor of variegation 3–9 homologue 1 (SUV39H1), SUV39H2, SET domain bifurcated 1 (SETDB1), SETDB2, G9A and G9A-like protein (GLP) — and the ‘readers’ of H3K9me2 or H3K9me3 are highly conserved and show considerable redundancy. Despite their redundancy, genetic ablation or mistargeting of an individual H3K9 methyltransferase can correlate with impaired cell differentiation, loss of tissue identity, premature aging and/or cancer. In this Review, we discuss recent advances in understanding the roles of the known H3K9-specific HMTs in ensuring transcriptional homeostasis during tissue differentiation in mammals. We examine the effects of H3K9-methylation-dependent gene repression in haematopoiesis, muscle differentiation and neurogenesis in mammals, and compare them with mechanistic insights obtained from the study of model organisms, notably Caenorhabditis elegans and Drosophila melanogaster. In all these organisms, H3K9-specific HMTs have both unique and redundant roles that ensure the maintenance of tissue integrity by restricting the binding of transcription factors to lineage-specific promoters and enhancer elements.

Histone H3 Lys9 (H3K9)-methylated heterochromatin ensures transcriptional silencing of repetitive elements and genes, and its deregulation leads to impaired cell and tissue identity, premature aging and cancer. Recent studies in mammals clarified the roles H3K9-specific histone methyltransferases in ensuring transcriptional homeostasis during tissue differentiation.

Prognostic Relevance of ZNF844 and Chr 19p13.2 KRAB-Zinc Finger Proteins in Clear Cell Renal Carcinoma

BACKGROUND/AIM

Clear-cell renal cell carcinoma (ccRCC) is the most common and aggressive form of all urological cancers, with poor prognosis and high mortality. Despite growing evidence of involvement in carcinogenesis, the role of KRAB-ZFP in ccRCC has not been fully explored. KRAB Zinc finger proteins (KRAB-ZFPs) are the largest family of mammalian transcription regulators. They are differentially expressed in various tissues during cellular development and phenotypic differentiation.

MATERIALS AND METHODS

In this study, the levels of transcripts of ccRCC from The Cancer Genome Atlas (TCGA) dataset were used to identify prognostic biomarkers in this disease.

RESULTS

Using bioinformatics techniques, we demonstrate that approximately 60% of KRAB zinc finger proteins located on chromosome 19p13.2 are differentially expressed, with all but two being down-regulated in ccRCC. Moreover, ZNF844, a paralog of ZNF433, was the most down-regulated across all histological grades and pathological stages (p<0.001). In addition, the decrease in ZNF844 expression was associated with poor patient survival (HR=0.41; 95% CI=0.3-0.56; p<0.0001). Gene Set Enrichment Analysis of genes inversely co-expressed with ZNF844 revealed that enriched pathways were consistently related to immune and translation processes (p<0.05, FDR <0.05). Lastly, ZNF844 expression showed moderate, inverse correlation to Helper T-cell (CD4 or Th1) subtype 1 (R=-0.558, p=5.15×10^-39) infiltration and with the exhausted T-cell phenotype (R=-0.37; p=4.1×10^-21).

CONCLUSION

Down-regulation of KRAB-ZFPs at 19p13.2 may represent a signature for ccRCC. Moreover, ZNF844 is a prognostic marker for ccRCC and may serve as a putative immune-related tumor suppressor gene.

Molecular profile of pure squamous cell carcinoma of the bladder identifies major roles for OSMR and YAP signalling

Pure squamous cell carcinoma (SCC) is the most common pure variant form of bladder cancer, found in 2–5% of cases. It often presents late and is unresponsive to cisplatin‐based chemotherapy. The molecular features of these tumours have not been elucidated in detail. We carried out whole‐exome sequencing (WES), copy number, and transcriptome analysis of bladder SCC. Muscle‐invasive bladder cancer (MIBC) samples with no evidence of squamous differentiation (non‐SD) were used for comparison. To assess commonality of features with urothelial carcinoma with SD, we examined data from SD samples in The Cancer Genome Atlas (TCGA) study of MIBC. TP53 was the most commonly mutated gene in SCC (64%) followed by FAT1 (45%). Copy number analysis revealed complex changes in SCC, many differing from those in samples with SD. Gain of 5p and 7p was the most common feature, and focal regions on 5p included OSMR and RICTOR. In addition to 9p deletions, we found some samples with focal gain of 9p24 containing CD274 (PD‐L1). Loss of 4q35 containing FAT1 was found in many samples such that all but one sample analysed by WES had FAT1 mutation or deletion. Expression features included upregulation of oncostatin M receptor (OSMR), metalloproteinases, metallothioneins, keratinisation genes, extracellular matrix components, inflammatory response genes, stem cell markers, and immune response modulators. Exploration of differentially expressed transcription factors identified BNC1 and TFAP2A, a gene repressed by PPARG, as the most upregulated factors. Known urothelial differentiation factors were downregulated along with 72 Kruppel‐associated (KRAB) domain‐containing zinc finger family protein (KZFP) genes. Novel therapies are urgently needed for these tumours. In addition to upregulated expression of EGFR, which has been suggested as a therapeutic target in basal/squamous bladder cancer, we identified expression signatures that indicate upregulated OSMR and YAP/TAZ signalling. Preclinical evaluation of the effects of inhibition of these pathways alone or in combination is merited.

KRAB-ZFPs and cancer stem cells identity

Studies on carcinogenesis continue to provide new information about different disease-related processes. Among others, much research has focused on the involvement of cancer stem cells (CSCs) in tumor initiation and progression. Studying the similarities and differences between CSCs and physiological stem cells (SCs) allows for a better understanding of cancer biology. Recently, it was shown that stem cell identity is partially governed by the Krϋppel-associated box domain zinc finger proteins (KRAB-ZFPs), the biggest family of transcription regulators. Several KRAB-ZFP factors exert a known effect in tumor cells, acting as tumor suppressor genes (TSGs) or oncogenes, yet their role in CSCs is still poorly characterized. Here, we review recent studies regarding the influence of KRAB-ZFPs and their cofactor protein TRIM28 on CSCs phenotype, stemness features, migration and invasion potential, metastasis, and expression of parental markers.

Host and Viral Zinc-Finger Proteins in COVID-19

An unprecedented effort to tackle the ongoing COVID-19 pandemic has characterized the activity of the global scientific community over the last two years. Hundreds of published studies have focused on the comprehension of the immune response to the virus and on the definition of the functional role of SARS-CoV-2 proteins. Proteins containing zinc fingers, both belonging to SARS-CoV-2 or to the host, play critical roles in COVID-19 participating in antiviral defenses and regulation of viral life cycle. Differentially expressed zinc finger proteins and their distinct activities could thus be important in determining the severity of the disease and represent important targets for drug development. Therefore, we here review the mechanisms of action of host and viral zinc finger proteins in COVID-19 as a contribution to the comprehension of the disease and also highlight strategies for therapeutic developments.

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.

A Multi-Omics Network of a Seven-Gene Prognostic Signature for Non-Small Cell Lung Cancer

There is an unmet clinical need to identify patients with early-stage non-small cell lung cancer (NSCLC) who are likely to develop recurrence and to predict their therapeutic responses. Our previous study developed a qRT-PCR-based seven-gene microfluidic assay to predict the recurrence risk and the clinical benefits of chemotherapy. This study showed it was feasible to apply this seven-gene panel in RNA sequencing profiles of The Cancer Genome Atlas (TCGA) NSCLC patients (n = 923) in randomly partitioned feasibility-training and validation sets (p < 0.05, Kaplan-Meier analysis). Using Boolean implication networks, DNA copy number variation-mediated transcriptional regulatory network of the seven-gene signature was identified in multiple NSCLC cohorts (n = 371). The multi-omics network genes, including PD-L1, were significantly correlated with immune infiltration and drug response to 10 commonly used drugs for treating NSCLC. ZNF71 protein expression was positively correlated with epithelial markers and was negatively correlated with mesenchymal markers in NSCLC cell lines in Western blots. PI3K was identified as a relevant pathway of proliferation networks involving ZNF71 and its isoforms formulated with CRISPR-Cas9 and RNA interference (RNAi) profiles. Based on the gene expression of the multi-omics network, repositioning drugs were identified for NSCLC treatment.

Loss of zinc-finger protein 212 leads to Purkinje cell death and locomotive abnormalities with phospholipase D3 downregulation

Although Krüppel-associated box domain-containing zinc-finger proteins (K-ZNFs) may be associated with sophisticated gene regulation in higher organisms, the physiological functions of most K-ZNFs remain unknown. The Zfp212 protein was highly conserved in mammals and abundant in the brain; it was mainly expressed in the cerebellum (Cb). Zfp212 (mouse homolog of human ZNF212) knockout (Zfp212-KO) mice showed a reduction in survival rate compared to wild-type mice after 20 months of age. GABAergic Purkinje cell degeneration in the Cb and aberrant locomotion were observed in adult Zfp212-KO mice. To identify genes related to the ataxia-like phenotype of Zfp212-KO mice, 39 ataxia-associated genes in the Cb were monitored. Substantial alterations in the expression of ataxin 10, protein phosphatase 2 regulatory subunit beta, protein kinase C gamma, and phospholipase D3 (Pld3) were observed. Among them, Pld3 alone was tightly regulated by Flag-tagged ZNF212 overexpression or Zfp212 knockdown in the HT22 cell line. The Cyclic Amplification and Selection of Targets assay identified the TATTTC sequence as a recognition motif of ZNF212, and these motifs occurred in both human and mouse PLD3 gene promoters. Adeno-associated virus-mediated introduction of human ZNF212 into the Cb of 3-week-old Zfp212-KO mice prevented Purkinje cell death and motor behavioral deficits. We confirmed the reduction of Zfp212 and Pld3 in the Cb of an alcohol-induced cerebellar degeneration mouse model, suggesting that the ZNF212–PLD3 relationship is important for Purkinje cell survival.

ZNF224 Protein: Multifaceted Functions Based on Its Molecular Partners

The transcription factor ZNF224 is a Kruppel-like zinc finger protein that consists of 707 amino acids and contains 19 tandemly repeated C₂H₂ zinc finger domains that mediate DNA binding and protein-protein interactions. ZNF224 was originally identified as a transcriptional repressor of genes involved in energy metabolism, and it was demonstrated that ZNF224-mediated transcriptional repression needs the interaction of its KRAB repressor domain with the co-repressor KAP1 and its zinc finger domains 1-3 with the arginine methyltransferase PRMT5. Furthermore, the protein ZNF255 was identified as an alternative isoform of ZNF224 that possesses different domain compositions mediating distinctive functional interactions. Subsequent studies showed that ZNF224 is a multifunctional protein able to exert different transcriptional activities depending on the cell context and the variety of its molecular partners. Indeed, it has been shown that ZNF224 can act as a repressor, an activator and a cofactor for other DNA-binding transcription factors in different human cancers. Here, we provide a brief overview of the current knowledge on the multifaceted interactions of ZNF224 and the resulting different roles of this protein in various cellular contexts.

Low Levels of TRIM28-Interacting KRAB-ZNF Genes Associate with Cancer Stemness and Predict Poor Prognosis of Kidney Renal Clear Cell Carcinoma Patients

Simple Summary

This is the first report investigating the involvement of TRIM28-interacting KRAB-ZNFs in kidney cancer progression. We demonstrate a significant negative association between KRAB-ZNFs and cancer stemness followed by an attenuated immune-suppressive response and reveal the prognostic role for several KRAB-ZNFs. Our findings may help better understand the molecular basis of kidney cancer and ultimately pave the way to more appropriate prognostic tools and novel therapeutic strategies directly eradicating the dedifferentiated compartment of the tumor.

Abstract

Krüppel-associated box zinc finger (KRAB-ZNF) proteins are known to regulate diverse biological processes, such as embryonic development, tissue-specific gene expression, and cancer progression. However, their involvement in the regulation of cancer stemness-like phenotype acquisition and maintenance is scarcely explored across solid tumor types, and to date, there are no data for kidney renal clear cell cancer (KIRC). We have harnessed The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database transcriptomic data and used several bioinformatic tools (i.e., GEPIA2, GSCALite, TISIDB, GSEA, CIBERSORT) to verify the relation between the expression and genomic alterations in KRAB-ZNFs and kidney cancer, focusing primarily on tumor dedifferentiation status and antitumor immune response. Our results demonstrate a significant negative correlation between KRAB-ZNFs and kidney cancer dedifferentiation status followed by an attenuated immune-suppressive response. The transcriptomic profiles of high KRAB-ZNF-expressing kidney tumors are significantly enriched with stem cell markers and show a depletion of several inflammatory pathways known for favoring cancer stemness. Moreover, we show for the first time the prognostic role for several KRAB-ZNFs in kidney cancer. Our results provide new insight into the role of selected KRAB-ZNF proteins in kidney cancer development. We believe that our findings may help better understand the molecular basis of KIRC.

Developing ZNF Gene Signatures Predicting Radiosensitivity of Patients with Breast Cancer

Adjuvant radiotherapy is one of the main treatment methods for breast cancer, but its clinical benefit depends largely on the characteristics of the patient. This study aimed to explore the relationship between the expression of zinc finger (ZNF) gene family proteins and the radiosensitivity of breast cancer patients. Clinical and gene expression data on a total of 976 breast cancer samples were obtained from The Cancer Genome Atlas (TCGA) database. ZNF gene expression was dichotomized into groups with a higher or lower level than the median level of expression. Univariate and multivariate Cox regression analyses were used to evaluate the relationship between ZNF gene expression levels and radiosensitivity. The Molecular Taxonomy Data of the International Federation of Breast Cancer (METABRIC) database was used for validation. The results revealed that 4 ZNF genes were possible radiosensitivity markers. High expression of ZNF644 and low expression levels of the other 3 genes (ZNF341, ZNF541, and ZNF653) were related to the radiosensitivity of breast cancer. Hierarchical cluster, Cox, and CoxBoost analysis based on these 4 ZNF genes indicated that patients with a favorable 4-gene signature had better overall survival on radiotherapy. Thus, this 4-gene signature may have value for selecting those patients most likely to benefit from radiotherapy. ZNF gene clusters could act as radiosensitivity signatures for breast cancer patients and may be involved in determining the radiosensitivity of cancer.

Transposable elements and their KZFP controllers are drivers of transcriptional innovation in the developing human brain

Transposable elements (TEs) account for more than 50% of the human genome and many have been co-opted throughout evolution to provide regulatory functions for gene expression networks. Several lines of evidence suggest that these networks are fine-tuned by the largest family of TE controllers, the KRAB-containing zinc finger proteins (KZFPs). One tissue permissive for TE transcriptional activation (termed "transposcription") is the adult human brain, however comprehensive studies on the extent of this process and its potential contribution to human brain development are lacking. To elucidate the spatiotemporal transposcriptome of the developing human brain, we have analyzed two independent RNA-seq data sets encompassing 16 brain regions from eight weeks postconception into adulthood. We reveal a distinct KZFP:TE transcriptional profile defining the late prenatal to early postnatal transition, and the spatiotemporal and cell type-specific activation of TE-derived alternative promoters driving the expression of neurogenesis-associated genes. Long-read sequencing confirmed these TE-driven isoforms as significant contributors to neurogenic transcripts. We also show experimentally that a co-opted antisense L2 element drives temporal protein relocalization away from the endoplasmic reticulum, suggestive of novel TE dependent protein function in primate evolution. This work highlights the widespread dynamic nature of the spatiotemporal KZFP:TE transcriptome and its importance throughout TE mediated genome innovation and neurotypical human brain development. To facilitate interactive exploration of these spatiotemporal gene and TE expression dynamics, we provide the "Brain TExplorer" web application freely accessible for the community.

Gv1, a Zinc Finger Gene Controlling Endogenous MLV Expression

The genomes of inbred mice harbor around 50 endogenous murine leukemia virus (MLV) loci, although the specific complement varies greatly between strains. The Gv1 locus is known to control the transcription of endogenous MLVs and to be the dominant determinant of cell-surface presentation of MLV envelope, the G_IX antigen. Here, we identify a single Krüppel-associated box zinc finger protein (ZFP) gene, Zfp998, as Gv1 and show it to be necessary and sufficient to determine the GIX+ phenotype. By long-read sequencing of bacterial artificial chromosome clones from 129 mice, the prototypic GIX+ strain, we reveal the source of sufficiency and deficiency as splice-acceptor variations and highlight the varying origins of the chromosomal region encompassing Gv1. Zfp998 becomes the second identified ZFP gene responsible for epigenetic suppression of endogenous MLVs in mice and further highlights the prominent role of this gene family in control of endogenous retroviruses.

Disruption of RING and PHD Domains of TRIM28 Evokes Differentiation in Human iPSCs

TRIM28, a multi-domain protein, is crucial in the development of mouse embryos and the maintenance of embryonic stem cells’ (ESC) self-renewal potential. As the epigenetic factor modulating chromatin structure, TRIM28 regulates the expression of numerous genes and is associated with progression and poor prognosis in many types of cancer. Because of many similarities between highly dedifferentiated cancer cells and normal pluripotent stem cells, we applied human induced pluripotent stem cells (hiPSC) as a model for stemness studies. For the first time in hiPSC, we analyzed the function of individual TRIM28 domains. Here we demonstrate the essential role of a really interesting new gene (RING) domain and plant homeodomain (PHD) in regulating pluripotency maintenance and self-renewal capacity of hiPSC. Our data indicate that mutation within the RING or PHD domain leads to the loss of stem cell phenotypes and downregulation of the FGF signaling. Moreover, impairment of RING or PHD domain results in decreased proliferation and impedes embryoid body formation. In opposition to previous data indicating the impact of phosphorylation on TRIM28 function, our data suggest that TRIM28 phosphorylation does not significantly affect the pluripotency and self-renewal maintenance of hiPSC. Of note, iPSC with disrupted RING and PHD functions display downregulation of genes associated with tumor metastasis, which are considered important targets in cancer treatment. Our data suggest the potential use of RING and PHD domains of TRIM28 as targets in cancer therapy.

A DNA Regulatory Element Haplotype at Zinc Finger Genes Is Associated with Host Resilience to Small Ruminant Lentivirus in Two Sheep Populations

Small ruminant lentivirus (SRLV) causes Maedi-Visna or Ovine Progressive Pneumonia in sheep and creates insidious livestock production losses. This retrovirus is closely related to human immunodeficiency virus and currently has no vaccines or cure. Genetic marker assisted selection for sheep disease resiliency presents an attractive management solution. Previously, we identified a region containing a cluster of zinc finger genes that had association with ovine SRLV proviral concentration. Trait-association analysis validated a small insertion/deletion variant near ZNF389 (rs397514112) in multiple sheep breeds. In the current study, 543 sheep from two distinct populations were genotyped at 34 additional variants for fine mapping of the regulatory elements within this locus. Variants were selected based on ChIP-seq annotation data from sheep alveolar macrophages that defined active cis-regulatory elements predicted to influence zinc finger gene expression. We present a haplotype block of variants within regulatory elements that have improved associations and larger effect sizes (up to 4.7-fold genotypic difference in proviral concentration) than the previously validated ZNF389 deletion marker. Hypotheses for the underlying causal mutation or mutations are presented based on changes to in silico transcription factor binding sites. These variants offer alternative markers for selective breeding and are targets for future functional mutation assays.