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

Biallelic loss-of-function variants in ZNF142 are associated with a robust DNA methylation signature affecting a limited number of genomic loci

Biallelic inactivating variants in ZNF142 underlie a clinically variable neurodevelopmental disorder. ZNF142 is a zinc-finger transcription factor with potential roles on chromatin organization, implying a possible association of ZNF142 loss of function with perturbed genome-wide DNA methylation (DNAm) pattern. We performed EPIC array-based methylation profiling of peripheral blood-derived DNA samples from 27 individuals with biallelic ZNF142 inactivating variants, together with 6 heterozygous carriers and 40 controls. A DNAm signature discovery pipeline was applied by using 440 controls for discovery and validation analyses, and a machine-learning model was trained to classify 8 individuals carrying ZNF142 variants of uncertain clinical significance. Analyses directed to explore the genome-wide DNAm landscape in affected individuals revealed 88 differentially methylated probes constituting the minimal informative set specific to ZNF142 loss of function. This reproducible pattern of DNAm changes involved regulatory regions of a small number of genes. The DNAm signature derived from peripheral blood allowed us to diagnose individuals carrying biallelic inactivating ZNF142 variants when applied to fibroblasts. Our findings provide evidence that biallelic loss-of-function ZNF142 variants result in a specific and robust DNAm signature. The identified DNAm pattern suggests occurrence of a methylation disturbance involving a small number of loci that appears to be shared by different cell lineages.

RETRACTED: 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.

Transposon-host arms race: a saga of genome evolution

Once considered 'junk DNA,' transposons or transposable elements (TEs) are now recognized as key drivers of genome evolution, contributing to genetic diversity, gene regulation, and species diversification. However, their ability to move within the genome poses a potential threat to genome integrity, promoting the evolution of robust host defense systems such as Krüppel-associated box (KRAB) domain-containing zinc finger proteins (KRAB-ZFPs), the human silencing hub (HUSH) complex, 4.5SH RNAs, and PIWI-interacting RNAs (piRNAs). This ongoing evolutionary arms race between TEs and host defenses continuously reshapes genome architecture and function. This review outlines various host defense mechanisms and explores the dynamic coevolution of TEs and host defenses in animals, highlighting how the defense mechanisms not only safeguard the host genomes but also drive genetic innovation through the arms race.

Transposable elements as instructors of the immune system

Transposable elements (TEs) are mobile repetitive nucleic acid sequences that have been incorporated into the genome through spontaneous integration, accounting for almost 50% of human DNA. Even though most TEs are no longer mobile today, studies have demonstrated that they have important roles in different biological processes, such as ageing, embryonic development, and cancer. TEs influence these processes through various mechanisms, including active transposition of TEs contributing to ongoing evolution, transposon transcription generating RNA or protein, and by influencing gene regulation as enhancers. However, how TEs interact with the immune system remains a largely unexplored field. In this Perspective, we describe how TEs might influence different aspects of the immune system, such as innate immune responses, T cell activation and differentiation, and tissue adaptation. Furthermore, TEs can serve as a source of neoantigens for T cells in antitumour immunity. We suggest that TE biology is an important emerging field of immunology and discuss the potential to harness the TE network therapeutically, for example, to improve immunotherapies for cancer and autoimmune and inflammatory diseases.

Multimeric transcription factor BCL11A utilizes two zinc-finger tandem arrays to bind clustered short sequence motifs

BCL11A, a transcription factor, is vital for hematopoiesis, including B and T cell maturation and the fetal-to-adult hemoglobin switch. Mutations in BCL11A are linked to neurodevelopmental disorders. BCL11A contains two DNA-binding zinc-finger arrays, low-affinity ZF2-3 and high-affinity ZF4-6, separated by a 300-amino-acid linker. ZF2-3 and ZF4-5 share 73% identity, including five out of six DNA base-interacting residues. These arrays bind similar short sequence motifs in clusters, with the linker enabling a broader binding span. Crystallographic structures of ZF4-6, in complex with oligonucleotides from the β-globin locus region, reveal DNA sequence recognition by residues Asn756 (ZF4), Lys784 and Arg787 (ZF5). A Lys784-to-Thr mutation, linked to a neurodevelopmental disorder with persistent fetal globin expression, reduces DNA binding over 10-fold but gains interaction with a variable base pair. BCL11A isoforms may form oligomers, enhancing chromatin occupancy and repressor functions by allowing multiple copies of both low- and high-affinity ZF arrays to bind DNA. These distinctive properties, apparently conserved among vertebrates, provide essential functional flexibility to this crucial regulator.

Transcriptomic characterisation of acute myeloid leukemia cell lines bearing the same t(9;11) driver mutation reveals different molecular signatures

BACKGROUND

Acute myeloid leukemia (AML) is the most common type of acute leukemia, accounting for 20% of cases in children and adolescents. Genome-wide studies have identified genes that are commonly mutated in AML, including many epigenetic regulators involved in either DNA methylation (DNMT3A, TET2, IDH1/2) or histone post-translational modifications (ASXL1, EZH2, MLL1). Several cell lines derived from AML patients are widely used in cancer research. Whether important differences in these cell lines exist remains poorly characterised.

RESULTS

Here, we used RNA sequencing (RNA-Seq) to contrast the transcriptome of four commonly used AML-derived cell lines: THP-1, NOMO-1, MOLM-13 bearing the common initiating t(9;11) translocation, and MV4.11 bearing the t(4;11) translocation. Gene set enrichment analyses and comparison of key transcription and epigenetic regulator genes revealed important differences in the transcriptome, distinguishing these AML models. Among these, we found striking differences in the expression of clusters of genes located on chromosome 19 encoding Zinc Finger (ZNF) transcriptional repressors. Low expression of many ZNF genes within these clusters is associated with poor survival in AML patients.

CONCLUSION

The present study offers a valuable resource by providing a detailed comparative characterisation of the transcriptome of cell lines within the same AML subtype used as models for leukemia research.

Young KRAB-zinc finger gene clusters are highly dynamic incubators of ERV-driven genetic heterogeneity in mice

KRAB-zinc finger proteins (KZFPs) comprise the largest family of mammalian transcription factors, rapidly evolving within and between species. Most KZFPs repress endogenous retroviruses (ERVs) and other retrotransposons, with KZFP gene numbers correlating with the ERV load across species, suggesting coevolution. How new KZFPs emerge in response to ERV invasions is currently unknown. Using a combination of long-read sequencing technologies and genome assembly, we present a first detailed comparative analysis of young KZFP gene clusters in the mouse lineage, which has undergone recent KZFP gene expansion and ERV infiltration. Detailed annotation of KZFP genes in a cluster on Mus musculus Chromosome 4 revealed parallel expansion and diversification of this locus in different mouse strains (C57BL/6J, 129S1/SvImJ and CAST/EiJ) and species (Mus spretus and Mus pahari). Our data supports a model by which new ERV integrations within young KZFP gene clusters likely promoted recombination events leading to the emergence of new KZFPs that repress them. At the same time, ERVs also increased their numbers by duplication instead of retrotransposition alone, unraveling a new mechanism for ERV enrichment at these loci.

p85β acts as a transcription cofactor and cooperates with BCLAF1 in the nucleus

p85β is a regulatory subunit of the phosphoinositide 3-kinase (PI3K). Emerging evidence suggests that p85β goes beyond its role in the PI3K and is functional in the nucleus. In this study, we discover that nuclear p85β is enriched at gene loci and regulates gene transcription and that this regulatory role contributes to the oncogenic potential of nuclear p85β. A multi-omics approach reveals the physical interaction and functional cooperativity between nuclear p85β and a transcription factor BCLAF1. We observe genome-wide co-occupancy of p85β and BCLAF1 at gene targets associated with transcriptional responses. Intriguingly, the targetome includes BCLAF1 of which transcription is activated by p85β and BCLAF1, indicating a positive autoregulation. While BCLAF1 recruits p85β to BCLAF1 loci, p85β facilitates the assembly of BCLAF1, the scaffold protein TRIM28 and the zinc finger transcription factor ZNF263, which together act in concert to activate BCLAF1 transcription. Collectively, this study provides functional evidence and mechanistic basis to support a role of nuclear p85β in modulating gene transcription.

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.

A whole-genome scan for evidence of recent positive and balancing selection in aye-ayes (Daubentonia madagascariensis) utilizing a well-fit evolutionary baseline model

The aye-aye (Daubentonia madagascariensis) is one of the 25 most endangered primate species in the world, maintaining amongst the lowest genetic diversity of any primate measured to date. Characterizing patterns of genetic variation within aye-aye populations, and the relative influences of neutral and selective processes in shaping that variation, is thus important for future conservation efforts. In this study, we performed the first whole-genome scans for recent positive and balancing selection in the species, utilizing high-coverage population genomic data from newly sequenced individuals. We generated null thresholds for our genomic scans by creating an evolutionarily appropriate baseline model that incorporates the demographic history of this aye-aye population, and identified a small number of candidate genes. Most notably, a suite of genes involved in olfaction - a key trait in these nocturnal primates - were identified as experiencing long-term balancing selection. We also conducted analyses to quantify the expected statistical power to detect positive and balancing selection in this population using site frequency spectrum-based inference methods, once accounting for the potentially confounding contributions of population history, recombination and mutation rate variation, and purifying and background selection. This work, presenting the first high-quality, genome-wide polymorphism data across the functional regions of the aye-aye genome, thus provides important insights into the landscape of episodic selective forces in this highly endangered species.

Advances of SS18-SSX fusion gene in synovial sarcoma: Emerging novel functions and therapeutic potentials

Synovial sarcoma is a rare type of soft tissue sarcoma that primarily affects adolescents and young adults, featured by aggressive behavior and a high potential for metastasis. Genetically, synovial sarcoma is defined by the fusion oncogene SS18-SSX arising from the translocation of t(X;18)(p11;q11). SS18-SSX fusion gene is the major driver of the oncogenic event in synovial sarcoma. SS18-SSX fusion protein, while not containing any DNA-binding motifs, binds to the SWI/SNF (BAF) complex, a major epigenetic regulator, leading to the disruption of gene expression which results in tumor initiation and progression. Emerging studies on the molecular mechanisms of SS18-SSX associated signaling pathway hold promise for developments in diagnosis and treatments. Advanced diagnostic methods facilitate early and precise detection of the tumor, enabling disease monitoring and prognostic improvements. Treatment of synovial sarcoma typically comprises local surgery, radiotherapy and chemotherapy, while novel managements such as immunotherapy, targeted therapies and epigenetic modifiers are explored. This review focuses on the recent studies of SS18-SSX fusion gene, epigenetic landscape, signaling pathways, diagnostic techniques, and relevant therapeutic advances, aiming to inhibit the oncogenic processes and improve outcomes for patients with synovial sarcoma.

Landscape of evolutionary arms races between transposable elements and KRAB-ZFP family

Transposable elements (TEs) are mobile parasitic sequences that have expanded within the host genome. It has been hypothesized that host organisms have expanded the Krüppel-associated box-containing zinc finger proteins (KRAB-ZFPs), which epigenetically suppress TEs, to counteract disorderly TE transpositions. This process is referred to as the evolutionary arms race. However, the extent to which this evolutionary arms race occurred across various TE families remains unclear. In the present study, we systematically explored the evolutionary arms race between TE families and human KRAB-ZFPs using public ChIP-seq data. We discovered and characterized new instances of evolutionary arms races with KRAB-ZFPs in endogenous retroviruses. Furthermore, we found that the regulatory landscape shaped by this arms race contributed to the gene regulatory networks. In summary, our results provide insight into the impact of the evolutionary arms race on TE families, the KRAB-ZFP family, and host gene regulatory networks.

Comprehensive Annotation and Expression Profiling of C2H2 Zinc Finger Transcription Factors across Chicken Tissues

As the most abundant class of transcription factors in eukaryotes, C2H2-type zinc finger proteins (C2H2-ZFPs) play critical roles in various biological processes. Despite being extensively studied in mammals, C2H2-ZFPs remain poorly characterized in birds. Recent accumulation of multi-omics data for chicken enables the genome-wide investigation of C2H2-ZFPs in birds. The purpose of this study is to reveal the genomic occurrence and evolutionary signature of chicken C2H2-ZFPs, and further depict their expression profiles across diverse chicken tissues. Here, we annotated 301 C2H2-ZFPs in chicken genome, which are associated with different effector domains, including KRAB, BTB, HOMEO, PHD, SCAN, and SET. Among them, most KRAB-ZFPs lack orthologues in mammals and tend to form clusters by duplication, supporting their fast evolution in chicken. We also annotated a unique and previously unidentified SCAN-ZFP, which is lineage-specific and highly expressed in ovary and testis. By integrating 101 RNA-seq datasets for 32 tissues, we found that most C2H2-ZFPs have tissue-specific expression. Particularly, 74 C2H2-ZFPs-including 27 KRAB-ZFPs-show blastoderm-enriched expression, indicating their association with early embryo development. Overall, this study performs comprehensive annotation and expression profiling of C2H2 ZFPs in diverse chicken tissues, which gives new insights into the evolution and potential function of C2H2-ZFPs in avian species.

Empowering brain tumor management: chimeric antigen receptor macrophage therapy

Brain tumors pose formidable challenges in oncology due to the intricate biology and the scarcity of effective treatment modalities. The emergence of immunotherapy has opened new avenues for innovative therapeutic strategies. Chimeric antigen receptor, originally investigated in T cell-based therapy, has now expanded to encompass macrophages, presenting a compelling avenue for augmenting anti-tumor immune surveillance. This emerging frontier holds promise for advancing the repertoire of therapeutic options against brain tumors, offering potential breakthroughs in combating the formidable malignancies of the central nervous system. Tumor-associated macrophages constitute a substantial portion, ranging from 30% to 50%, of the tumor tissue and exhibit tumor-promoting phenotypes within the immune-compromised microenvironment. Constructing CAR-macrophages can effectively repolarize M2-type macrophages towards an M1-type phenotype, thereby eliciting potent anti-tumor effects. CAR-macrophages can recruit T cells to the brain tumor site, thereby orchestrating a remodeling of the immune niche to effectively inhibit tumor growth. In this review, we explore the potential limitations as well as strategies for optimizing CAR-M therapy, offering insights into the future direction of this innovative therapeutic approach.

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.

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.

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.