Deep vertebrate roots for mammalian zinc finger transcription factor subfamilies

Identification of ZNF384 as a regulator of epigenome in leukemia

Leukemia is a complex hematologic cancer driven by genetic and epigenetic changes that impact gene expression. Understanding these molecular mechanisms is essential for improving leukemia diagnosis and prognosis. This study examines the role of the zinc finger protein ZNF384 in the epigenome and its influence on gene regulation in leukemia. We analyzed next-generation sequencing data from The Encyclopedia of DNA Elements (ENCODE), integrating datasets such as chromatin immunoprecipitation sequencing (ChIP-seq) of ZNF384 and regulatory histone marks, RNA sequencing (RNA-seq), and Hi-C data from K562 and GM12878 cells. Additionally, we used RNA-seq from K562 ZNF384 knock-down (KD) cells generated via CRISPR interference (CRISPRi) to validate our findings. This enabled us to explore the chromatin interaction patterns of ZNF384 and its regulatory impact. Our results demonstrate that ZNF384 associates with promoters and enhancers in K562 and GM12878 cells, facilitating increased transcription levels. We also found ZNF384 enriched at topologically associating domain (TAD) boundaries and chromatin loops, suggesting a role in three-dimensional (3D) chromatin organization. Furthermore, we identified a significant binding of ZNF384 at SINE-Alu elements in both K562 and GM12878 cells. In summary, this study highlights the regulatory role of ZNF384 in the leukemia epigenome and its impact on gene expression. Understanding the oncogenic implications of ZNF384 may improve leukemia diagnosis and prognosis.

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.

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.

Gene expression networks regulated by human personality

Genome-wide association studies of human personality have been carried out, but transcription of the whole genome has not been studied in relation to personality in humans. We collected genome-wide expression profiles of adults to characterize the regulation of expression and function in genes related to human personality. We devised an innovative multi-omic approach to network analysis to identify the key control elements and interactions in multi-modular networks. We identified sets of transcribed genes that were co-expressed in specific brain regions with genes known to be associated with personality. Then we identified the minimum networks for the co-localized genes using bioinformatic resources. Subjects were 459 adults from the Young Finns Study who completed the Temperament and Character Inventory and provided peripheral blood for genomic and transcriptomic analysis. We identified an extrinsic network of 45 regulatory genes from seed genes in brain regions involved in self-regulation of emotional reactivity to extracellular stimuli (e.g., self-regulation of anxiety) and an intrinsic network of 43 regulatory genes from seed genes in brain regions involved in self-regulation of interpretations of meaning (e.g., production of concepts and language). We discovered that interactions between the two networks were coordinated by a control hub of 3 miRNAs and 3 protein-coding genes shared by both. Interactions of the control hub with proteins and ncRNAs identified more than 100 genes that overlap directly with known personality-related genes and more than another 4000 genes that interact indirectly. We conclude that the six-gene hub is the crux of an integrative network that orchestrates information-transfer throughout a multi-modular system of over 4000 genes enriched in liquid-liquid-phase-separation (LLPS)-related RNAs, diverse transcription factors, and hominid-specific miRNAs and lncRNAs. Gene expression networks associated with human personality regulate neuronal plasticity, epigenesis, and adaptive functioning by the interactions of salience and meaning in self-awareness.

Rapid adaptation of cellular metabolic rate to the MicroRNA complements of mammals and its relevance to the evolution of endothermy

The metabolic efficiency of mammalian cells depends on the attenuation of intrinsic translation noise by microRNAs. We devised a metric of cellular metabolic rate (cMR), rMR/Mexp optimally fit to the number of microRNA families (mirFam), that is robust to variation in mass and sensitive to body temperature (Tb), consistent with the heat dissipation limit theory of Speakman and Król (2010). Using mirFam as predictor, an Ornstein-Uhlenbeck process of stabilizing selection, with an adaptive shift at the divergence of Boreoeutheria, accounted for 95% of the variation in cMR across mammals. Branchwise rates of evolution of cMR, mirFam and Tb concurrently increased 6- to 7-fold at the divergence of Boreoeutheria, independent of mass. Cellular MR variation across placental mammals was also predicted by the sum of model conserved microRNA-target interactions, revealing an unexpected degree of integration of the microRNA-target apparatus into the energy economy of the mammalian cell.

A zinc finger transcription factor enables social behaviors while controlling transposable elements and immune response in prefrontal cortex

The neurobiological origins of social behaviors are incompletely understood. Here we utilized synthetic biology approaches to reprogram the function of ZFP189, a transcription factor whose expression and function in rodent prefrontal cortex was previously demonstrated to be protective against stress-induced social deficits. We created novel synthetic ZFP189 transcription factors including ZFP189VPR, which activates the transcription of target genes and therefore exerts opposite functional control from the endogenous, transcriptionally repressive ZFP189WT. Following viral delivery of these synthetic ZFP189 transcription factors to mouse prefrontal cortex, we observe that ZFP189-mediated transcriptional control promotes mature dendritic spine morphology on transduced pyramidal neurons. Interestingly, inversion of ZFP189-mediated transcription in this brain area, achieved by viral delivery of synthetic ZFP189VPR, precipitates social behavioral deficits in terms of social interaction, motivation, and the cognition necessary for the maintenance of social hierarchy, without other observable behavioral deficits. RNA sequencing of virally manipulated prefrontal cortex tissues reveals that ZFP189 transcription factors of opposing regulatory function (ZFP189WT versus ZFP189VPR) have opposite influence on the expression of genetic transposable elements as well as genes that participate in adaptive immune functions. Collectively, this work reveals that ZFP189 function in the prefrontal cortex coordinates structural and transcriptional neuroadaptations necessary for complex social behaviors while regulating transposable element-rich regions of DNA and the expression of immune-related genes. Given the evidence for a co-evolution of social behavior and the brain immune response, we posit that ZFP189 may have evolved to augment brain transposon-associated immune function as a way of enhancing an animal's capacity for functioning in social groups.

Primate-specific ZNF808 is essential for pancreatic development in humans

Identifying genes linked to extreme phenotypes in humans has the potential to highlight biological processes not shared with all other mammals. Here, we report the identification of homozygous loss-of-function variants in the primate-specific gene ZNF808 as a cause of pancreatic agenesis. ZNF808 is a member of the KRAB zinc finger protein family, a large and rapidly evolving group of epigenetic silencers which target transposable elements. We show that loss of ZNF808 in vitro results in aberrant activation of regulatory potential contained in the primate-specific transposable elements it represses during early pancreas development. This leads to inappropriate specification of cell fate with induction of genes associated with liver identity. Our results highlight the essential role of ZNF808 in pancreatic development in humans and the contribution of primate-specific regions of the human genome to congenital developmental disease.

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.

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

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 zinc finger transcription factor tunes social behaviors by controlling transposable elements and immune response in prefrontal cortex

The neurobiological origins of social behaviors are incompletely understood. Here we utilized synthetic biology approaches to reprogram the function of ZFP189, a transcription factor whose expression and function in the rodent prefrontal cortex was previously determined to be protective against stress-induced social deficits. We created novel synthetic ZFP189 transcription factors including ZFP189VPR, which activates the transcription of target genes and therefore exerts opposite functional control from the endogenous, transcriptionally repressive ZFP189WT. Upon viral delivery of these synthetic ZFP189 transcription factors to mouse prefrontal cortex, we observe that ZFP189-mediated transcriptional control promotes mature dendritic spine morphology on transduced pyramidal neurons. Interestingly, dysregulation of ZFP189-mediated transcription in this brain area, achieved by delivery of synthetic ZFP189VPR, precipitates social behavioral deficits in terms of social interaction, motivation, and the cognition necessary for the maintenance of social hierarchy, without other observable behavioral deficits. By performing RNA sequencing in virally manipulated prefrontal cortex tissues, we discover that ZFP189 transcription factors of opposing regulatory function have opposite influence on the expression of genetic transposable elements as well as genes that participate in immune functions. Collectively, this work reveals that ZFP189 function in the prefrontal cortex coordinates structural and transcriptional neuroadaptations necessary for social behaviors by binding transposable element-rich regions of DNA to regulate immune-related genes. Given the evidence for a co-evolution of social behavior and the brain immune response, we posit that ZFP189 may have evolved to augment brain transposon-associated immune function as a way of enhancing an animal's capacity for functioning in social groups.

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.

Genome-wide identification of copy number variation and association with fat deposition in thin and fat-tailed sheep breeds

Copy number variants (CNVs) are a type of genetic polymorphism which contribute to phenotypic variation in several species, including livestock. In this study, we used genomic data of 192 animals from 3 Iranian sheep breeds including 96 Baluchi sheep and 47 Lori-Bakhtiari sheep as fat-tailed breeds and 47 Zel sheep as thin-tailed sheep breed genotyped with Illumina OvineSNP50K Beadchip arrays. Also, for association test, 70 samples of Valle del Belice sheep were added to the association test as thin-tailed sheep breed. PennCNV and CNVRuler software were, respectively, used to study the copy number variation and genomic association analyses. We detected 573 and 242 CNVs in the fat and thin tailed breeds, respectively. In terms of CNV regions (CNVRs), these represented 328 and 187 CNVRs that were within or overlapping with 790 known Ovine genes. The CNVRs covered approximately 73.85 Mb of the sheep genome with average length 146.88 kb, and corresponded to 2.6% of the autosomal genome sequence. Five CNVRs were randomly chosen for validation, of which 4 were experimentally confirmed using Real time qPCR. Functional enrichment analysis showed that genes harbouring CNVs in thin-tailed sheep were involved in the adaptive immune response, regulation of reactive oxygen species biosynthetic process and response to starvation. In fat-tailed breeds these genes were involved in cellular protein modification process, regulation of heart rate, intestinal absorption, olfactory receptor activity and ATP binding. Association test identified one copy gained CNVR on chromosomes 6 harbouring two protein-coding genes HGFAC and LRPAP1. Our findings provide information about genomic structural changes and their association to the interested traits including fat deposition and environmental compatibility in sheep.

Active endogenous retroviral elements in human pluripotent stem cells play a role in regulating host gene expression

Human endogenous retroviruses, also called LTR elements, can be bound by transcription factors and marked by different histone modifications in different biological contexts. Recently, individual LTR or certain subclasses of LTRs such as LTR7/HERVH and LTR5_Hs/HERVK families have been identified as cis-regulatory elements. However, there are still many LTR elements with unknown functions. Here, we dissected the landscape of histone modifications and regulatory map of LTRs by integrating 98 ChIP-seq data in human embryonic stem cells (ESCs), and annotated the active LTRs enriching enhancer/promoter-related histone marks. Notably, we found that MER57E3 functionally acted as proximal regulatory element to activate respective ZNF gene. Additionally, HERVK transcript could mainly function in nucleus to activate the adjacent genes. Since LTR5_Hs/LTR5 was bound by many early embryo-specific transcription factors, we further investigated the expression dynamics in different pluripotent states. LTR5_Hs/LTR5/HERVK exhibited higher expression level in naïve ESCs and extended pluripotent stem cells (EPSCs). Functionally, the LTR5_Hs/LTR5 with high activity could serve as a distal enhancer to regulate the host genes. Ultimately, our study not only provides a comprehensive regulatory map of LTRs in human ESCs, but also explores the regulatory models of MER57E3 and LTR5_Hs/LTR5 in host genome.

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.

KRAB zinc finger protein ZNF676 controls the transcriptional influence of LTR12-related endogenous retrovirus sequences

BACKGROUND

Transposable element-embedded regulatory sequences (TEeRS) and their KRAB-containing zinc finger protein (KZFP) controllers are increasingly recognized as modulators of gene expression. We aim to characterize the contribution of this system to gene regulation in early human development and germ cells.

RESULTS

Here, after studying genes driven by the long terminal repeat (LTR) of endogenous retroviruses, we identify the ape-restricted ZNF676 as the sequence-specific repressor of a subset of contemporary LTR12 integrants responsible for a large fraction of transpochimeric gene transcripts (TcGTs) generated during human early embryogenesis. We go on to reveal that the binding of this KZFP correlates with the epigenetic marking of these TEeRS in the germline, and is crucial to the control of genes involved in ciliogenesis/flagellogenesis, a biological process that dates back to the last common ancestor of eukaryotes.

CONCLUSION

These results illustrate how KZFPs and their TE targets contribute to the evolutionary turnover of transcription networks and participate in the transgenerational inheritance of epigenetic traits.

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.

Transposable elements shape the evolution of mammalian development

Transposable elements (TEs) promote genetic innovation but also threaten genome stability. Despite multiple layers of host defence, TEs actively shape mammalian-specific developmental processes, particularly during pre-implantation and extra-embryonic development and at the maternal-fetal interface. Here, we review how TEs influence mammalian genomes both directly by providing the raw material for genetic change and indirectly via co-evolving TE-binding Krüppel-associated box zinc finger proteins (KRAB-ZFPs). Throughout mammalian evolution, individual activities of ancient TEs were co-opted to enable invasive placentation that characterizes live-born mammals. By contrast, the widespread activity of evolutionarily young TEs may reflect an ongoing co-evolution that continues to impact mammalian development.

ZNF507 affects TGF-β signaling via TGFBR1 and MAP3K8 activation in the progression of prostate cancer to an aggressive state

Background

The progression of prostate cancer (PC) to the highly aggressive metastatic castration-resistant prostate cancer (mCRPC) or neuroendocrine prostate cancer (NEPC) is a fatal condition and the underlying molecular mechanisms are poorly understood. Here, we identified the novel transcriptional factor ZNF507 as a key mediator in the progression of PC to an aggressive state.

Methods

We analyzed ZNF507 expression in the data from various human PC database and high-grade PC patient samples. By establishment of ZNF507 knockdown and overexpression human PC cell lines, we assessed in vitro PC phenotype changes including cell proliferation, survival, migration and invasion. By performing microarray with ZNF507 knockdown PC cells, we profiled the gene clusters affected by ZNF507 knockdown. Moreover, ZNF507 regulated key signal was evaluated by dual-luciferase reporter and chromatin immunoprecipitation (ChIP) assays. Finally, we performed xenograft and in vivo metastasis assay to confirm the effect of ZNF507 knockdown in PC cells.

Results

We found that ZNF507 expression was increased, particularly in the highly graded PC. ZNF507 was also found to be associated with metastatic PC of a high grade. Loss- or gain-of-function–based analysis revealed that ZNF507 promotes the growth, survival, proliferation, and metastatic properties of PC (e.g., epithelial-mesenchymal transition) by upregulating TGF-β signaling. Profiling of gene clusters affected by ZNF507 knockdown revealed that ZNF507 positively regulated the transcription of TGFBR1, MAP3K8, and FURIN, which in turn promoted the progression of PC to highly metastatic and aggressive state.

Conclusions

Our findings suggest that ZNF507 is a novel key regulator of TGF-β signaling in the progression of malignant PC and could be a promising target for studying the development of advanced metastatic PCs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-02094-3.

Baby Genomics: Tracing the Evolutionary Changes That Gave Rise to Placentation

It has long been challenging to uncover the molecular mechanisms behind striking morphological innovations such as mammalian pregnancy. We studied the power of a robust comparative orthology pipeline based on gene synteny to address such problems. We inferred orthology relations between human genes and genes from each of 43 other vertebrate genomes, resulting in ∼18,000 orthologous pairs for each genome comparison. By identifying genes that first appear coincident with origin of the placental mammals, we hypothesized that we would define a subset of the genome enriched for genes that played a role in placental evolution. We thus pinpointed orthologs that appeared before and after the divergence of eutherian mammals from marsupials. Reinforcing previous work, we found instead that much of the genetic toolkit of mammalian pregnancy evolved through the repurposing of preexisting genes to new roles. These genes acquired regulatory controls for their novel roles from a group of regulatory genes, many of which did in fact originate at the appearance of the eutherians. Thus, orthologs appearing at the origin of the eutherians are enriched in functions such as transcriptional regulation by Krüppel-associated box-zinc-finger proteins, innate immune responses, keratinization, and the melanoma-associated antigen protein class. Because the cellular mechanisms of invasive placentae are similar to those of metastatic cancers, we then used our orthology inferences to explore the association between placenta invasion and cancer metastasis. Again echoing previous work, we find that genes that are phylogenetically older are more likely to be implicated in cancer development.

Adaptive venom evolution and toxicity in octopods is driven by extensive novel gene formation, expansion, and loss

BACKGROUND

Cephalopods represent a rich system for investigating the genetic basis underlying organismal novelties. This diverse group of specialized predators has evolved many adaptations including proteinaceous venom. Of particular interest is the blue-ringed octopus genus (Hapalochlaena), which are the only octopods known to store large quantities of the potent neurotoxin, tetrodotoxin, within their tissues and venom gland.

FINDINGS

To reveal genomic correlates of organismal novelties, we conducted a comparative study of 3 octopod genomes, including the Southern blue-ringed octopus (Hapalochlaena maculosa). We present the genome of this species and reveal highly dynamic evolutionary patterns at both non-coding and coding organizational levels. Gene family expansions previously reported in Octopus bimaculoides (e.g., zinc finger and cadherins, both associated with neural functions), as well as formation of novel gene families, dominate the genomic landscape in all octopods. Examination of tissue-specific genes in the posterior salivary gland revealed that expression was dominated by serine proteases in non-tetrodotoxin-bearing octopods, while this family was a minor component in H. maculosa. Moreover, voltage-gated sodium channels in H. maculosa contain a resistance mutation found in pufferfish and garter snakes, which is exclusive to the genus. Analysis of the posterior salivary gland microbiome revealed a diverse array of bacterial species, including genera that can produce tetrodotoxin, suggestive of a possible production source.

CONCLUSIONS

We present the first tetrodotoxin-bearing octopod genome H. maculosa, which displays lineage-specific adaptations to tetrodotoxin acquisition. This genome, along with other recently published cephalopod genomes, represents a valuable resource from which future work could advance our understanding of the evolution of genomic novelty in this family.

Multiple Isolated Transcription Factors Act as Switches and Contribute to Species Uniqueness

Mammals have variable numbers (1300-2000) of transcription factors (TFs), but the reasons for this large variation are unclear. To investigate general TF patterns, we de novo identified 156,906 TFs from 96 mammalian species. We identified more than 500 human isolated TFs that are rarely reported in human TF-to-TF networks. Mutations in the genes of these TFs were less lethal than those of connected TFs. Consequently, these isolated TFs are more tolerant of changes and have become unique during speciation. They may also serve as a source of variation for TF evolution. Reconciliation of TF-family phylogenetic trees with a mammalian species tree revealed an average of 37.8% TF gains and 15.0% TF losses over 177 million years, which implies that isolated TFs are pervasive in mammals. Compared with non-TF interacting genes, TF-interacting genes have unique TF profiles and have higher expression levels in mice than in humans. Different expression levels of the same TF-interacting gene contribute to species-specific phenotypes. Formation and loss of isolated TFs enabling unique TF profiles may provide variable switches that adjust divergent expression profiles of target genes to generate species-specific phenotypes, thereby making species unique.

Identification of an Individualized Prognostic Signature Based on the RWSR Model in Early-Stage Bladder Carcinoma

Bladder cancer (BLCA) is the fourth common cancer among males in the United States, which is also the fourth leading cause of cancer-related death in old males. BLCA has a high recurrence rate, with over 50% of patients which has at least one recurrence within five years. Due to the complexity of the molecular mechanisms and heterogeneous cancer feature, BLCA clinicians find it hard to make an efficient management decision as they lack reliable assessment of mortality risk. Meanwhile, there is currently no screening suitable prognostic signature or method recommended for early detection, which is significantly important to early-stage detection and prognosis. In this study, a novel model, named the risk-weighted sparse regression (RWSR) model, is constructed to identify a robust signature for patients of early-stage BLCA. The 17-gene signature is generated and then validated as an independent prognostic factor in BLCA cohorts from GSE13507 and TCGA_BLCA datasets. Meanwhile, a risk score model is developed and validated among the 17-gene signature. The risk score is also considered an independent factor for prognosis prediction, which is confirmed through prognosis analysis. The Kaplan-Meier with the log-rank test is used to assess survival difference. Furthermore, the predictive capacity of the signature is proved through stratification analysis. Finally, an effective patient classification is completed by a combination of the 17-gene signature and stage information, which is for better survival prediction and treatment decisions. Besides, 11 genes in the signature, such as coiled-coil domain containing 73 (CCDC73) and protein kinase, DNA-activated, and catalytic subunit (PRKDC), are proved to be prognosis marker genes or strongly associated with prognosis and progress of other types of cancer in published literature already. As a result, this paper would more accurately predict a patient's prognosis and improve surveillance in the clinical setting, which may provide a quantitative and reliable decision-making basis for the treatment plan.

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

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

DUF3669, a "domain of unknown function" within ZNF746 and ZNF777, oligomerizes and contributes to transcriptional repression

Background

ZNF746 and ZNF777 belong to a subset of the large Krüppel-associated box (KRAB) zinc finger (ZNF) transcription factor family. They contain, like four other members in human, an additional conserved domain, the “domain of unknown function 3669” (DUF3669). Previous work on members of this subfamily suggested involvement in transcriptional regulation and aberrant ZNF746 overexpression leads to neuronal cell death in Parkinson’s disease.

Results

Here we demonstrate that N-terminal protein segments of the ZNF746a major isoform and ZNF777 act in concert to exert moderate transcriptional repression activities. Full potency depended on the intact configuration consisting of DUF3669, a variant KRAB domain and adjacent sequences. While DUF3669 contributes an intrinsic weak inhibitory activity, the isolated KRAB-AB domains did not repress. Importantly, DUF3669 provides a novel protein-protein interaction interface and mediates direct physical interaction between the members of the subfamily in oligomers. The ZNF746 protein segment encoded by exons 5 and 6 boosted repressor potency, potentially due to the presence of an acceptor lysine for sumoylation at K189. Repressor activity of the potent canonical ZNF10 KRAB domain was not augmented by heterologous transfer of DUF3669, pointing to the importance of context for DUF3669’s impact on transcription. Neither ZNF746a nor ZNF777 protein segments stably associated with TRIM28 within cells. Isoform ZNF746b that contains, unlike the major isoform, a full-length KRAB-A subdomain, displayed substantially increased repressor potency. This increase is due to canonical mechanisms known for KRAB domains since it did not take place in HAP1 knockout models of TRIM28 and SETDB1. A glycine to glutamic acid replacement that complies with a bona fide conserved “MLE” sequence within KRAB-A led to a further strong gain in repressor potency to levels comparable to those of the canonical ZNF10 KRAB domain. Each gain of repressive activity was accompanied by an enhanced interaction with TRIM28 protein.

Conclusion

DUF3669 adds a protein-protein interaction surface to a subgroup of KRAB-ZNF proteins within an N-terminal configuration with variant KRAB and adjacent sequences likely regulated by sumoylation. DUF3669 contributes to transcriptional repression strength and its homo- and hetero-oligomerization characteristics probably extended the regulatory repertoire of KRAB-ZNF transcription factors during amniote evolution.

The interactome of KRAB zinc finger proteins reveals the evolutionary history of their functional diversification

Krüppel-associated box (KRAB)-containing zinc finger proteins (KZFPs) are encoded in the hundreds by the genomes of higher vertebrates, and many act with the heterochromatin-inducing KAP1 as repressors of transposable elements (TEs) during early embryogenesis. Yet, their widespread expression in adult tissues and enrichment at other genetic loci indicate additional roles. Here, we characterized the protein interactome of 101 of the ~350 human KZFPs. Consistent with their targeting of TEs, most KZFPs conserved up to placental mammals essentially recruit KAP1 and associated effectors. In contrast, a subset of more ancient KZFPs rather interacts with factors related to functions such as genome architecture or RNA processing. Nevertheless, KZFPs from coelacanth, our most distant KZFP-encoding relative, bind the cognate KAP1. These results support a hypothetical model whereby KZFPs first emerged as TE-controlling repressors, were continuously renewed by turnover of their hosts' TE loads, and occasionally produced derivatives that escaped this evolutionary flushing by development and exaptation of novel functions.

The KRAB-zinc-finger protein ZFP708 mediates epigenetic repression at RMER19B retrotransposons

Global epigenetic reprogramming is vital to purge germ cell-specific epigenetic features to establish the totipotent state of the embryo. This process transpires to be carefully regulated and is not an undirected, radical erasure of parental epigenomes. The TRIM28 complex has been shown to be crucial in embryonic epigenetic reprogramming by regionally opposing DNA demethylation to preserve vital parental information to be inherited from germline to soma. Yet the DNA-binding factors guiding this complex to specific targets are largely unknown. Here, we uncover and characterize a novel, maternally expressed, TRIM28-interacting KRAB zinc-finger protein: ZFP708. It recruits the repressive TRIM28 complex to RMER19B retrotransposons to evoke regional heterochromatin formation. ZFP708 binding to these hitherto unknown TRIM28 targets is DNA methylation and H3K9me3 independent. ZFP708 mutant mice are viable and fertile, yet embryos fail to inherit and maintain DNA methylation at ZFP708 target sites. This can result in activation of RMER19B-adjacent genes, while ectopic expression of ZFP708 results in transcriptional repression. Finally, we describe the evolutionary conservation of ZFP708 in mice and rats, which is linked to the conserved presence of the targeted RMER19B retrotransposons in these species.

Summary: Analysis of the function and targets of a maternal KRAB-zinc-finger protein, ZFP708, found to specifically mediate maintenance of DNA methylation at a subset of LTR retrotransposons during embryonic epigenetic reprogramming.

Similarity regression predicts evolution of transcription factor sequence specificity

Transcription factor (TF) binding specificities (motifs) are essential for the analysis of gene regulation. Accurate prediction of TF motifs is critical, because it is infeasible to assay all TFs in all sequenced eukaryotic genomes. There is ongoing controversy regarding the degree of motif diversification among related species that is, in part, because of uncertainty in motif prediction methods. Here we describe similarity regression, a significantly improved method for predicting motifs, which we use to update and expand the Cis-BP database. Similarity regression inherently quantifies TF motif evolution, and shows that previous claims of near-complete conservation of motifs between human and Drosophila are inflated, with nearly half of the motifs in each species absent from the other, largely due to extensive divergence in C2H2 zinc finger proteins. We conclude that diversification in DNA-binding motifs is pervasive, and present a new tool and updated resource to study TF diversity and gene regulation across eukaryotes.

Hominoid-Specific Transposable Elements and KZFPs Facilitate Human Embryonic Genome Activation and Control Transcription in Naive Human ESCs

Expansion of transposable elements (TEs) coincides with evolutionary shifts in gene expression. TEs frequently harbor binding sites for transcriptional regulators, thus enabling coordinated genome-wide activation of species- and context-specific gene expression programs, but such regulation must be balanced against their genotoxic potential. Here, we show that Krüppel-associated box (KRAB)-containing zinc finger proteins (KZFPs) control the timely and pleiotropic activation of TE-derived transcriptional cis regulators during early embryogenesis. Evolutionarily recent SVA, HERVK, and HERVH TE subgroups contribute significantly to chromatin opening during human embryonic genome activation and are KLF-stimulated enhancers in naive human embryonic stem cells (hESCs). KZFPs of corresponding evolutionary ages are simultaneously induced and repress the transcriptional activity of these TEs. Finally, the same KZFP-controlled TE-based enhancers later serve as developmental and tissue-specific enhancers. Thus, by controlling the transcriptional impact of TEs during embryogenesis, KZFPs facilitate their genome-wide incorporation into transcriptional networks, thereby contributing to human genome regulation.

•

KLFs foster EGA by activating enhancers embedded in young TEs (TEENhancers)

•

TEENhancers confer a degree of species specificity to early genome activation

•

TEENhancers stimulate the expression of KZFPs responsible for their repression

•

These KZFPs in turn facilitate TEENhancers’ exaptation as tissue-specific regulators

ZFP30 promotes adipogenesis through the KAP1-mediated activation of a retrotransposon-derived Pparg2 enhancer

Krüppel-associated box zinc finger proteins (KZFPs) constitute the largest family of mammalian transcription factors, but most remain completely uncharacterized. While initially proposed to primarily repress transposable elements, recent reports have revealed that KFZPs contribute to a wide variety of other biological processes. Using murine and human in vitro and in vivo models, we demonstrate here that one poorly studied KZFP, ZFP30, promotes adipogenesis by directly targeting and activating a retrotransposon-derived Pparg2 enhancer. Through mechanistic studies, we further show that ZFP30 recruits the co-regulator KRAB-associated protein 1 (KAP1), which, surprisingly, acts as a ZFP30 co-activator in this adipogenic context. Our findings provide an understanding of both adipogenic and KZFP-KAP1 complex-mediated gene regulation, showing that the KZFP-KAP1 axis can also function in a non-repressive manner.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Endogenous Retroviruses Function as Gene Expression Regulatory Elements During Mammalian Pre-implantation Embryo Development

Pre-implantation embryo development encompasses several key developmental events, especially the activation of zygotic genome activation (ZGA)-related genes. Endogenous retroviruses (ERVs), which are regarded as “deleterious genomic parasites”, were previously considered to be “junk DNA”. However, it is now known that ERVs, with limited conservatism across species, mediate conserved developmental processes (e.g., ZGA). Transcriptional activation of ERVs occurs during the transition from maternal control to zygotic genome control, signifying ZGA. ERVs are versatile participants in rewiring gene expression networks during epigenetic reprogramming. Particularly, a subtle balance exists between ERV activation and ERV repression in host–virus interplay, which leads to stage-specific ERV expression during pre-implantation embryo development. A large portion of somatic cell nuclear transfer (SCNT) embryos display developmental arrest and ZGA failure during pre-implantation embryo development. Furthermore, because of the close relationship between ERV activation and ZGA, exploring the regulatory mechanism underlying ERV activation may also shed more light on the enigma of SCNT embryo development in model animals.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation

BACKGROUND

Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes.

FINDINGS

Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated.

CONCLUSION

The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

ZNF445 is a primary regulator of genomic imprinting

Genomic imprinting is an epigenetic process regulated by germline-derived DNA methylation, causing parental origin-specific monoallelic gene expression. Zinc finger protein 57 (ZFP57) is critical for maintenance of this epigenetic memory during post-fertilization reprogramming, yet incomplete penetrance of ZFP57 mutations in humans and mice suggests additional effectors. We reveal that ZNF445/ZFP445, which we trace to the origins of imprinting, binds imprinting control regions (ICRs) in mice and humans. In mice, ZFP445 and ZFP57 act together, maintaining all but one ICR in vivo, whereas earlier embryonic expression of ZNF445 and its intolerance to loss-of-function mutations indicate greater importance in the maintenance of human imprints.

Complete Disruption of Autism-Susceptibility Genes by Gene Editing Predominantly Reduces Functional Connectivity of Isogenic Human Neurons

Autism spectrum disorder (ASD) is phenotypically and genetically heterogeneous. We present a CRISPR gene editing strategy to insert a protein tag and premature termination sites creating an induced pluripotent stem cell (iPSC) knockout resource for functional studies of ten ASD-relevant genes (AFF2/FMR2, ANOS1, ASTN2, ATRX, CACNA1C, CHD8, DLGAP2, KCNQ2, SCN2A, TENM1). Neurogenin 2 (NGN2)-directed induction of iPSCs allowed production of excitatory neurons, and mutant proteins were not detectable. RNA sequencing revealed convergence of several neuronal networks. Using both patch-clamp and multi-electrode array approaches, the electrophysiological deficits measured were distinct for different mutations. However, they culminated in a consistent reduction in synaptic activity, including reduced spontaneous excitatory postsynaptic current frequencies in AFF2/FMR2-, ASTN2-, ATRX-, KCNQ2-, and SCN2A-null neurons. Despite ASD susceptibility genes belonging to different gene ontologies, isogenic stem cell resources can reveal common functional phenotypes, such as reduced functional connectivity.

•

iPSC knockout resource for functional studies of ten ASD-risk genes

•

Disruption of common transcriptional networks associated with neurons and synapses

•

Reduced synaptic activity commonly observed for functionally diverse ASD-risk genes

The Cephalopod Large Brain Enigma: Are Conserved Mechanisms of Stem Cell Expansion the Key?

Within the clade of mollusks, cephalopods have developed an unusually large and complex nervous system. The increased complexity of the cephalopod centralized "brain" parallels an amazing amount of complex behaviors that culminate in one order, the octopods. The mechanisms that enable evolution of expanded brains in invertebrates remain enigmatic. While expression mapping of known molecular pathways demonstrated the conservation of major neurogenesis pathways and revealed neurogenic territories, it did not explain why cephalopods could massively increase their brain size compared to other mollusks. Such an increase is reminiscent of the expansion of the cerebral cortex in mammalians, which have enlarged their number and variety of neurogenic stem cells. We hypothesize that similar mechanisms might be at play in cephalopods and that focusing on the stem cell biology of cephalopod neurogenesis and genetic innovations might be smarter strategies to uncover the mechanism that has driven cephalopod brain expansion.

Two critical positions in zinc finger domains are heavily mutated in three human cancer types

A major goal of cancer genomics is to identify somatic mutations that play a role in tumor initiation or progression. Somatic mutations within transcription factors are of particular interest, as gene expression dysregulation is widespread in cancers. The substantial gene expression variation evident across tumors suggests that numerous regulatory factors are likely to be involved and that somatic mutations within them may not occur at high frequencies across patient cohorts, thereby complicating efforts to uncover which ones are cancer-relevant. Here we analyze somatic mutations within the largest family of human transcription factors, namely those that bind DNA via Cys2His2 zinc finger domains. Specifically, to hone in on important mutations within these genes, we aggregated somatic mutations across all of them by their positions within Cys2His2 zinc finger domains. Remarkably, we found that for three classes of cancers profiled by The Cancer Genome Atlas (TCGA)-Uterine Corpus Endometrial Carcinoma, Colon and Rectal Adenocarcinomas, and Skin Cutaneous Melanoma-two specific, functionally important positions within zinc finger domains are mutated significantly more often than expected by chance, with alterations in 18%, 10% and 43% of tumors, respectively. Numerous zinc finger genes are affected, with those containing Krüppel-associated box (KRAB) repressor domains preferentially targeted by these mutations. Further, the genes with these mutations also have high overall missense mutation rates, are expressed at levels comparable to those of known cancer genes, and together have biological process annotations that are consistent with roles in cancers. Altogether, we introduce evidence broadly implicating mutations within a diverse set of zinc finger proteins as relevant for cancer, and propose that they contribute to the widespread transcriptional dysregulation observed in cancer cells.

Interrogating the Functions of PRDM9 Domains in Meiosis

Homologous recombination is required for proper segregation of homologous chromosomes during meiosis. It occurs predominantly at recombination hotspots that are defined by the DNA binding specificity of the PRDM9 protein. PRDM9 contains three conserved domains typically involved in regulation of transcription; yet, the role of PRDM9 in gene expression control is not clear. Here, we analyze the germline transcriptome of Prdm9^−/− male mice in comparison to Prdm9^+/+ males and find no apparent differences in the mRNA and miRNA profiles. We further explore the role of PRDM9 in meiosis by analyzing the effect of the KRAB, SSXRD, and post-SET zinc finger deletions in a cell culture expression system and the KRAB domain deletion in mice. We found that although the post-SET zinc finger and the KRAB domains are not essential for the methyltransferase activity of PRDM9 in cell culture, the KRAB domain mutant mice show only residual PRDM9 methyltransferase activity and undergo meiotic arrest. In aggregate, our data indicate that domains typically involved in regulation of gene expression do not serve that role in PRDM9, but are likely involved in setting the proper chromatin environment for initiation and completion of homologous recombination.

ZSCAN5B and primate-specific paralogs bind RNA polymerase III genes and extra-TFIIIC (ETC) sites to modulate mitotic progression

Mammalian genomes contain hundreds of genes transcribed by RNA Polymerase III (Pol III), encoding noncoding RNAs and especially the tRNAs specialized to carry specific amino acids to the ribosome for protein synthesis. In addition to this well-known function, tRNAs and their genes (tDNAs) serve a variety of other critical cellular functions. For example, tRNAs and other Pol III transcripts can be cleaved to yield small RNAs with potent regulatory activities. Furthermore, from yeast to mammals, active tDNAs and related “extra-TFIIIC” (ETC) loci provide the DNA scaffolds for the most ancient known mechanism of three-dimensional chromatin architecture. Here we identify the ZSCAN5 TF family - including mammalian ZSCAN5B and its primate-specific paralogs - as proteins that occupy mammalian Pol III promoters and ETC sites. We show that ZSCAN5B binds with high specificity to a conserved subset of Pol III genes in human and mouse. Furthermore, primate-specific ZSCAN5A and ZSCAN5D also bind Pol III genes, although ZSCAN5D preferentially localizes to MIR SINE- and LINE2-associated ETC sites. ZSCAN5 genes are expressed in proliferating cell populations and are cell-cycle regulated, and siRNA knockdown experiments suggested a cooperative role in regulation of mitotic progression. Consistent with this prediction, ZSCAN5A knockdown led to increasing numbers of cells in mitosis and the appearance of cells. Together, these data implicate the role of ZSCAN5 genes in regulation of Pol III genes and nearby Pol II loci, ultimately influencing cell cycle progression and differentiation in a variety of tissues.

Individual retrotransposon integrants are differentially controlled by KZFP/KAP1-dependent histone methylation, DNA methylation and TET-mediated hydroxymethylation in naïve embryonic stem cells

Background

The KZFP/KAP1 (KRAB zinc finger proteins/KRAB-associated protein 1) system plays a central role in repressing transposable elements (TEs) and maintaining parent-of-origin DNA methylation at imprinting control regions (ICRs) during the wave of genome-wide reprogramming that precedes implantation. In naïve murine embryonic stem cells (mESCs), the genome is maintained highly hypomethylated by a combination of TET-mediated active demethylation and lack of de novo methylation, yet KAP1 is tethered by sequence-specific KZFPs to ICRs and TEs where it recruits histone and DNA methyltransferases to impose heterochromatin formation and DNA methylation.

Results

Here, upon removing either KAP1 or the cognate KZFP, we observed rapid TET2-dependent accumulation of 5hmC at both ICRs and TEs. In the absence of the KZFP/KAP1 complex, ICRs lost heterochromatic histone marks and underwent both active and passive DNA demethylation. For KAP1-bound TEs, 5mC hydroxylation correlated with transcriptional reactivation. Using RNA-seq, we further compared the expression profiles of TEs upon Kap1 removal in wild-type, Dnmt and Tet triple knockout mESCs. While we found that KAP1 represents the main effector of TEs repression in all three settings, we could additionally identify specific groups of TEs further controlled by DNA methylation. Furthermore, we observed that in the absence of TET proteins, activation upon Kap1 depletion was blunted for some TE integrants and increased for others.

Conclusions

Our results indicate that the KZFP/KAP1 complex maintains heterochromatin and DNA methylation at ICRs and TEs in naïve embryonic stem cells partly by protecting these loci from TET-mediated demethylation. Our study further unveils an unsuspected level of complexity in the transcriptional control of the endovirome by demonstrating often integrant-specific differential influences of histone-based heterochromatin modifications, DNA methylation and 5mC oxidation in regulating TEs expression.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0177-1) contains supplementary material, which is available to authorized users.

Mammalian transposable elements and their impacts on genome evolution

Transposable elements (TEs) are genetic elements with the ability to mobilize and replicate themselves in a genome. Mammalian genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. As a result, TEs have had significant impacts on mammalian evolution. Here we summarize the current understanding of TE content in mammal genomes and find that, with a few exceptions, most fall within a predictable range of observations. First, one third to one half of the genome is derived from TEs. Second, most mammalian genomes are dominated by LINE and SINE retrotransposons, more limited LTR retrotransposons, and minimal DNA transposon accumulation. Third, most mammal genome contains at least one family of actively accumulating retrotransposon. Finally, horizontal transfer of TEs among lineages is rare. TE exaptation events are being recognized with increasing frequency. Despite these beneficial aspects of TE content and activity, the majority of TE insertions are neutral or deleterious. To limit the deleterious effects of TE proliferation, the genome has evolved several defense mechanisms that act at the epigenetic, transcriptional, and post-transcriptional levels. The interaction between TEs and these defense mechanisms has led to an evolutionary arms race where TEs are suppressed, evolve to escape suppression, then are suppressed again as the defense mechanisms undergo compensatory change. The result is complex and constantly evolving interactions between TEs and host genomes.

KRAB zinc finger proteins

Krüppel-associated box domain zinc finger proteins (KRAB-ZFPs) are the largest family of transcriptional regulators in higher vertebrates. Characterized by an N-terminal KRAB domain and a C-terminal array of DNA-binding zinc fingers, they participate, together with their co-factor KAP1 (also known as TRIM28), in repression of sequences derived from transposable elements (TEs). Until recently, KRAB-ZFP/KAP1-mediated repression of TEs was thought to lead to irreversible silencing, and the evolutionary selection of KRAB-ZFPs was considered to be just the host component of an arms race against TEs. However, recent advances indicate that KRAB-ZFPs and their TE targets also partner up to establish species-specific regulatory networks. Here, we provide an overview of the KRAB-ZFP gene family, highlighting how its evolutionary history is linked to that of TEs, and how KRAB-ZFPs influence multiple aspects of development and physiology.

The Role of KRAB-ZFPs in Transposable Element Repression and Mammalian Evolution

Kruppel-associated box zinc-finger proteins (KRAB-ZFPs) make up the largest family of transcription factors in humans. These proteins emerged in the last common ancestor of coelacanth and tetrapods, and have expanded and diversified in the mammalian lineage. Although their mechanism of transcriptional repression has been well studied for over a decade, the DNA-binding activities and the biological functions of these proteins have been largely unexplored. Recent large-scale ChIP-seq studies and loss-of-function experiments have revealed that KRAB-ZFPs play a major role in the recognition and transcriptional silencing of transposable elements (TEs), consistent with an 'arms race model' of KRAB-ZFP evolution against invading TEs. However, this model is insufficient to explain the evolution of many KRAB-ZFPs that appear to domesticate TEs for novel host functions. We highlight some of the mammalian regulatory innovations driven by specific KRAB-ZFPs, including genomic imprinting, meiotic recombination hotspot choice, and placental growth.

New factors in mammalian DNA repair-the chromatin connection

In response to DNA damage mammalian cells activate a complex network of stress response pathways collectively termed DNA damage response (DDR). DDR involves a temporary arrest of the cell cycle to allow for the repair of the damage. DDR also attenuates gene expression by silencing global transcription and translation. Main function of DDR is, however, to prevent the fixation of debilitating changes to DNA by activation of various DNA repair pathways. Proper execution of DDR requires careful coordination between these interdependent cellular responses. Deregulation of some aspects of DDR orchestration is potentially pathological and could lead to various undesired outcomes such as DNA translocations, cellular transformation or acute cell death. It is thus critical to understand the regulation of DDR in cells especially in the light of a strong linkage between the DDR impairment and the occurrence of common human diseases such as cancer. In this review we focus on recent advances in understanding of mammalian DNA repair regulation and a on the function of PAXX/c9orf142 and ZNF281 proteins that recently had been discovered to play a role in that process. We focus on regulation of double-strand DNA break (DSB) repair via the non-homologous end joining pathway, as unrepaired DSBs are the primary cause of pathological cellular states after DNA damage. Interestingly these new factors operate at the level of chromatin, which reinforces a notion of a central role of chromatin structure in the regulation of cellular DDR regulation.

C2H2 Zinc Finger Proteins: The Largest but Poorly Explored Family of Higher Eukaryotic Transcription Factors

The emergence of whole-genome assays has initiated numerous genome-wide studies of transcription factor localizations at genomic regulatory elements (enhancers, promoters, silencers, and insulators), as well as facilitated the uncovering of some of the key principles of chromosomal organization. However, the proteins involved in the formation and maintenance of the chromosomal architecture and the organization of regulatory domains remain insufficiently studied. This review attempts to collate the available data on the abundant but still poorly understood family of proteins with clusters of the C2H2 zinc finger domains. One of the best known proteins of this family is a well conserved protein known as CTCF, which plays a key role in the establishment of the chromosomal architecture in vertebrates. The distinctive features of C2H2 zinc finger proteins include strong and specific binding to a long and unique DNA recognition target sequence and rapid expansion within various animal taxa during evolution. The reviewed data support a proposed model according to which many of the C2H2 proteins have functions that are similar to those of the CTCF in the organization of the chromatin architecture.