Multiparameter functional diversity of human C2H2 zinc finger proteins

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

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

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

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

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

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

Structures and biological functions of zinc finger proteins and their roles in hepatocellular carcinoma

Zinc finger proteins are transcription factors with the finger domain, which plays a significant role in gene regulation. As the largest family of transcription factors in the human genome, zinc finger (ZNF) proteins are characterized by their different DNA binding motifs, such as C2H2 and Gag knuckle. Different kinds of zinc finger motifs exhibit a wide variety of biological functions. Zinc finger proteins have been reported in various diseases, especially in several cancers. Hepatocellular carcinoma (HCC) is the third leading cause of cancer-associated death worldwide, especially in China. Most of HCC patients have suffered from hepatitis B virus (HBV) and hepatitis C virus (HCV) injection for a long time. Although the surgical operation of HCC has been extremely developed, the prognosis of HCC is still very poor, and the underlying mechanisms in HCC tumorigenesis are still not completely understood. Here, we summarize multiple functions and recent research of zinc finger proteins in HCC tumorigenesis and progression. We also discuss the significance of zinc finger proteins in HCC diagnosis and prognostic evaluation.

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

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

ZFTA Translocations Constitute Ependymoma Chromatin Remodeling and Transcription Factors

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

KLF10 integrates circadian timing and sugar signaling to coordinate hepatic metabolism

The mammalian circadian timing system and metabolism are highly interconnected, and disruption of this coupling is associated with negative health outcomes. Krüppel-like factors (KLFs) are transcription factors that govern metabolic homeostasis in various organs. Many KLFs show a circadian expression in the liver. Here, we show that the loss of the clock-controlled KLF10 in hepatocytes results in extensive reprogramming of the mouse liver circadian transcriptome, which in turn alters the temporal coordination of pathways associated with energy metabolism. We also show that glucose and fructose induce Klf10, which helps mitigate glucose intolerance and hepatic steatosis in mice challenged with a sugar beverage. Functional genomics further reveal that KLF10 target genes are primarily involved in central carbon metabolism. Together, these findings show that in the liver KLF10 integrates circadian timing and sugar metabolism-related signaling, and serves as a transcriptional brake that protects against the deleterious effects of increased sugar consumption.

Loss of ZNF677 Expression Is an Independent Predictor for Distant Metastasis in Middle Eastern Papillary Thyroid Carcinoma Patients

Thyroid cancer incidence has increased in recent decades. Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. Approximately 30% of PTC patients develop recurrence or distant metastasis and tend to have poor prognosis. Therefore, the identification of targetable biomarkers in this subset of patients is of great importance. Accumulating evidence indicates that zinc finger protein 677 (ZNF677), which belongs to the zinc finger protein family, is an important effector during the progression of multiple malignancies. However, its role in Middle Eastern PTC patients has not been fully illustrated. Here, we uncovered the molecular mechanism and the clinical impact of ZNF677 expression in a large cohort of more than 1200 Middle Eastern PTC and 15 metastatic tissues. We demonstrated that ZNF677 is frequently downregulated in primary PTC (13.6%, 168/1235) and showed that complete loss of expression of ZNF677 is significantly associated with aggressive clinico-pathological markers such as extrathyroidal extension (p = 0.0008) and distant metastases (p < 0.0001). We also found a significantly higher incidence of ZNF677 loss in primary tumors with distant metastases (33.3%; p < 0.0001) as well as in distant metastatic tissues (46.7%; p = 0.0002) compared to the overall cohort (13.6%). More importantly, PTC with loss of ZNF677 expression showed significantly lower metastasis-free survival (p = 0.0090). Interestingly, on multivariate logistic regression analysis, ZNF677 loss was an independent predictor of distant metastasis in PTC (Odds ratio = 2.60, 95% Confidence interval = 1.20–5.62, p = 0.0155). In addition, we found a significant association between ZNF677 loss and phospho-AKT expression (p < 0.0001). Our functional molecular results suggest that ZNF677 acts as a tumor suppressor, mediating its effect by inhibiting AKT phosphorylation. Taken together, our results highlight the pivotal role played by ZNF677 during carcinogenesis and metastasis formation in Middle Eastern PTC patients.

Functional characterization of RebL1 highlights the evolutionary conservation of oncogenic activities of the RBBP4/7 orthologue in Tetrahymena thermophila

Retinoblastoma-binding proteins 4 and 7 (RBBP4 and RBBP7) are two highly homologous human histone chaperones. They function in epigenetic regulation as subunits of multiple chromatin-related complexes and have been implicated in numerous cancers. Due to their overlapping functions, our understanding of RBBP4 and 7, particularly outside of Opisthokonts, has remained limited. Here, we report that in the ciliate protozoan Tetrahymena thermophila a single orthologue of human RBBP4 and 7 proteins, RebL1, physically interacts with histone H4 and functions in multiple epigenetic regulatory pathways. Functional proteomics identified conserved functional links for Tetrahymena RebL1 protein as well as human RBBP4 and 7. We found that putative subunits of multiple chromatin-related complexes including CAF1, Hat1, Rpd3, and MuvB, co-purified with RebL1 during Tetrahymena growth and conjugation. Iterative proteomics analyses revealed that the cell cycle regulatory MuvB-complex in Tetrahymena is composed of at least five subunits including evolutionarily conserved Lin54, Lin9 and RebL1 proteins. Genome-wide analyses indicated that RebL1 and Lin54 (Anqa1) bind within genic and intergenic regions. Moreover, Anqa1 targets primarily promoter regions suggesting a role for Tetrahymena MuvB in transcription regulation. RebL1 depletion inhibited cellular growth and reduced the expression levels of Anqa1 and Lin9. Consistent with observations in glioblastoma tumors, RebL1 depletion suppressed DNA repair protein Rad51 in Tetrahymena, thus underscoring the evolutionarily conserved functions of RBBP4/7 proteins. Our results suggest the essentiality of RebL1 functions in multiple epigenetic regulatory complexes in which it impacts transcription regulation and cellular viability.

ZNF37A promotes tumor metastasis through transcriptional control of THSD4/TGF-β axis in colorectal cancer

Poorly differentiated colorectal cancer (CRC) is characterized by aggressive invasion and stromal fibroblast activation, which results in rapid progression and poor therapeutic consequences. However, the regulatory mechanism involved remains unclear. Here, we showed that ZNF37A, a member of KRAB-ZFP family, was upregulated in poorly differentiated CRCs and associated with tumor metastasis. ZNF37A enhanced the metastatic potential of multiple CRC cell lines and promoted distant metastasis in an orthotopic CRC model. Further investigation attributed the ZNF37A-exacerbated metastasis to increased extracellular TGF-β and the consequent activation of cancer-associated fibroblasts (CAFs) in tumor microenvironment (TME). Mechanistically, ZNF37A formed a complex with KAP1 and bound to the promoter of THSD4, a TME modulator, to suppress its transcription, which is required for ZNF37A-mediated TGF-β activation and CRC metastasis. Collectively, our study indicates that ZNF37A promotes TGF-β signaling in CRC cells and activates CAFs by transcriptionally repressing THSD4 to drive CRC metastasis, implicating ZNF37A as a potential biomarker for CRC differentiation and progression.

ZNF91 deletion in human embryonic stem cells leads to ectopic activation of SVA retrotransposons and up-regulation of KRAB zinc finger gene clusters

Transposable element (TE) invasions have shaped vertebrate genomes over the course of evolution. They have contributed an extra layer of species-specific gene regulation by providing novel transcription factor binding sites. In humans, SINE-VNTR-Alu (SVA) elements are one of three still active TE families; approximately 2800 SVA insertions exist in the human genome, half of which are human-specific. TEs are often silenced by KRAB zinc finger (KZNF) proteins recruiting corepressor proteins that establish a repressive chromatin state. A number of KZNFs have been reported to bind SVAs, but their individual contribution to repressing SVAs and their roles in suppressing SVA-mediated gene-regulatory effects remains elusive. We analyzed the genome-wide binding profile for ZNF91 in human cells and found that ZNF91 interacts with the VNTR region of SVAs. Through CRISPR-Cas9-mediated deletion of ZNF91 in human embryonic stem cells, we established that loss of ZNF91 results in increased transcriptional activity of SVAs. In contrast, SVA activation was not observed upon genetic deletion of the ZNF611 gene encoding another strong SVA interactor. Epigenetic profiling confirmed the loss of SVA repression in the absence of ZNF91 and revealed that mainly evolutionary young SVAs gain gene activation-associated epigenetic modifications. Genes close to activated SVAs showed a mild up-regulation, indicating SVAs adopt properties of cis-regulatory elements in the absence of repression. Notably, genome-wide derepression of SVAs elicited the communal up-regulation of KZNFs that reside in KZNF clusters. This phenomenon may provide new insights into the potential mechanisms used by the host genome to sense and counteract TE invasions.

Conserved paradoxical relationships among the evolutionary, structural and expressional features of KRAB zinc-finger proteins reveal their special functional characteristics

Background

One striking feature of the large KRAB domain-containing zinc finger protein (KZFP) family is its rapid evolution, leading to hundreds of member genes with various origination time in a certain mammalian genome. However, a comprehensive genome-wide and across-taxa analysis of the structural and expressional features of KZFPs with different origination time is lacking. This type of analysis will provide valuable clues about the functional characteristics of this special family.

Results

In this study, we found several conserved paradoxical phenomena about this issue. 1) Ordinary young domains/proteins tend to be disordered, but most of KRAB domains are completely structured in 64 representative species across the superclass of Sarcopterygii and most of KZFPs are also highly structured, indicating their rigid and unique structural and functional characteristics; as exceptions, old-zinc-finger-containing KZFPs have relatively disordered KRAB domains and linker regions, contributing to diverse interacting partners and functions. 2) In general, young or highly structured proteins tend to be spatiotemporal specific and have low abundance. However, by integrated analysis of 29 RNA-seq datasets, including 725 samples across early embryonic development, embryonic stem cell differentiation, embryonic and adult organs, tissues in 7 mammals, we found that KZFPs tend to express ubiquitously with medium abundance regardless of evolutionary age and structural disorder degree, indicating the wide functional requirements of KZFPs in various states. 3) Clustering and correlation analysis reveal that there are differential expression patterns across different spatiotemporal states, suggesting the specific-high-expression KZFPs may play important roles in the corresponding states. In particular, part of young-zinc-finger-containing KZFPs are highly expressed in early embryonic development and ESCs differentiation into endoderm or mesoderm. Co-expression analysis revealed that young-zinc-finger-containing KZFPs are significantly enriched in five co-expression modules. Among them, one module, including 13 young-zinc-finger-containing KZFPs, showed an ‘early-high and late-low’ expression pattern. Further functional analysis revealed that they may function in early embryonic development and ESC differentiation via participating in cell cycle related processes.

Conclusions

This study shows the conserved and special structural, expressional features of KZFPs, providing new clues about their functional characteristics and potential causes of their rapid evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00346-w.

CTCF As an Example of DNA-Binding Transcription Factors Containing Clusters of C2H2-Type Zinc Fingers

In mammals, most of the boundaries of topologically associating domains and all well-studied insulators are rich in binding sites for the CTCF protein. According to existing experimental data, CTCF is a key factor in the organization of the architecture of mammalian chromosomes. A characteristic feature of the CTCF is that the central part of the protein contains a cluster consisting of eleven domains of C2H2-type zinc fingers, five of which specifically bind to a long DNA sequence conserved in most animals. The class of transcription factors that carry a cluster of C2H2-type zinc fingers consisting of five or more domains (C2H2 proteins) is widely represented in all groups of animals. The functions of most C2H2 proteins still remain unknown. This review presents data on the structure and possible functions of these proteins, using the example of the vertebrate CTCF protein and several well- characterized C2H2 proteins in Drosophila and mammals.

Structure-function relationships explain CTCF zinc finger mutation phenotypes in cancer

CCCTC-binding factor (CTCF) plays fundamental roles in transcriptional regulation and chromatin architecture maintenance. CTCF is also a tumour suppressor frequently mutated in cancer, however, the structural and functional impact of mutations have not been examined. We performed molecular and structural characterisation of five cancer-specific CTCF missense zinc finger (ZF) mutations occurring within key intra- and inter-ZF residues. Functional characterisation of CTCF ZF mutations revealed a complete (L309P, R339W, R377H) or intermediate (R339Q) abrogation as well as an enhancement (G420D) of the anti-proliferative effects of CTCF. DNA binding at select sites was disrupted and transcriptional regulatory activities abrogated. Molecular docking and molecular dynamics confirmed that mutations in residues specifically contacting DNA bases or backbone exhibited loss of DNA binding. However, R339Q and G420D were stabilised by the formation of new primary DNA bonds, contributing to gain-of-function. Our data confirm that a spectrum of loss-, change- and gain-of-function impacts on CTCF zinc fingers are observed in cell growth regulation and gene regulatory activities. Hence, diverse cellular phenotypes of mutant CTCF are clearly explained by examining structure-function relationships.

High-Throughput Discovery and Characterization of Human Transcriptional Effectors

Thousands of proteins localize to the nucleus; however, it remains unclear which contain transcriptional effectors. Here, we develop HT-recruit, a pooled assay where protein libraries are recruited to a reporter, and their transcriptional effects are measured by sequencing. Using this approach, we measure gene silencing and activation for thousands of domains. We find a relationship between repressor function and evolutionary age for the KRAB domains, discover that Homeodomain repressor strength is collinear with Hox genetic organization, and identify activities for several domains of unknown function. Deep mutational scanning of the CRISPRi KRAB maps the co-repressor binding surface and identifies substitutions that improve stability/silencing. By tiling 238 proteins, we find repressors as short as ten amino acids. Finally, we report new activator domains, including a divergent KRAB. These results provide a resource of 600 human proteins containing effectors and demonstrate a scalable strategy for assigning functions to protein domains.

Exome Sequencing Reveals a Phenotype Modifying Variant in ZNF528 in Primary Osteoporosis With a COL1A2 Deletion

We studied a family with severe primary osteoporosis carrying a heterozygous p.Arg8Phefs*14 deletion in COL1A2, leading to haploinsufficiency. Three affected individuals carried the mutation and presented nearly identical spinal fractures but lacked other typical features of either osteogenesis imperfecta or Ehlers-Danlos syndrome. Although mutations leading to haploinsufficiency in COL1A2 are rare, mutations in COL1A1 that lead to less protein typically result in a milder phenotype. We hypothesized that other genetic factors may contribute to the severe phenotype in this family. We performed whole-exome sequencing in five family members and identified in all three affected individuals a rare nonsense variant (c.1282C > T/p.Arg428*, rs150257846) in ZNF528. We studied the effect of the variant using qPCR and Western blot and its subcellular localization with immunofluorescence. Our results indicate production of a truncated ZNF528 protein that locates in the cell nucleus as per the wild-type protein. ChIP and RNA sequencing analyses on ZNF528 and ZNF528-c.1282C > T indicated that ZNF528 binding sites are linked to pathways and genes regulating bone morphology. Compared with the wild type, ZNF528-c.1282C > T showed a global shift in genomic binding profile and pathway enrichment, possibly contributing to the pathophysiology of primary osteoporosis. We identified five putative target genes for ZNF528 and showed that the expression of these genes is altered in patient cells. In conclusion, the variant leads to expression of truncated ZNF528 and a global change of its genomic occupancy, which in turn may lead to altered expression of target genes. ZNF528 is a novel candidate gene for bone disorders and may function as a transcriptional regulator in pathways affecting bone morphology and contribute to the phenotype of primary osteoporosis in this family together with the COL1A2 deletion. © 2020 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Role of CCCH-Type Zinc Finger Proteins in Human Adenovirus Infections

The zinc finger proteins make up a significant part of the proteome and perform a huge variety of functions in the cell. The CCCH-type zinc finger proteins have gained attention due to their unusual ability to interact with RNA and thereby control different steps of RNA metabolism. Since virus infections interfere with RNA metabolism, dynamic changes in the CCCH-type zinc finger proteins and virus replication are expected to happen. In the present review, we will discuss how three CCCH-type zinc finger proteins, ZC3H11A, MKRN1, and U2AF1, interfere with human adenovirus replication. We will summarize the functions of these three cellular proteins and focus on their potential pro- or anti-viral activities during a lytic human adenovirus infection.

An efficient KRAB domain for CRISPRi applications in human cells

Clustered regularly interspaced short palindromic repeat interference (CRISPRi), based on the fusion of inactive Cas9 (dCas9) to the Krüppel-associated box (KRAB) repressor, is a powerful platform for silencing gene expression. However, it suffers from incomplete silencing of target genes. We assayed 57 KRAB domains for their repressive potency and identified the ZIM3 KRAB domain as an exceptionally potent repressor. We establish that ZIM3 KRAB-dCas9 fusion silences gene expression more efficiently than existing platforms.

Silencing and Transcriptional Regulation of Endogenous Retroviruses: An Overview

Almost half of the human genome is made up of transposable elements (TEs), and about 8% consists of endogenous retroviruses (ERVs). ERVs are remnants of ancient exogenous retrovirus infections of the germ line. Most TEs are inactive and not detrimental to the host. They are tightly regulated to ensure genomic stability of the host and avoid deregulation of nearby gene loci. Histone-based posttranslational modifications such as H3K9 trimethylation are one of the main silencing mechanisms. Trim28 is one of the identified master regulators of silencing, which recruits most prominently the H3K9 methyltransferase Setdb1, among other factors. Sumoylation and ATP-dependent chromatin remodeling factors seem to contribute to proper localization of Trim28 to ERV sequences and promote Trim28 interaction with Setdb1. Additionally, DNA methylation as well as RNA-mediated targeting of TEs such as piRNA-based silencing play important roles in ERV regulation. Despite the involvement of ERV overexpression in several cancer types, autoimmune diseases, and viral pathologies, ERVs are now also appreciated for their potential positive role in evolution. ERVs can provide new regulatory gene elements or novel binding sites for transcription factors, and ERV gene products can even be repurposed for the benefit of the host.

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.

The UCSC repeat browser allows discovery and visualization of evolutionary conflict across repeat families

Background

Nearly half the human genome consists of repeat elements, most of which are retrotransposons, and many of which play important biological roles. However repeat elements pose several unique challenges to current bioinformatic analyses and visualization tools, as short repeat sequences can map to multiple genomic loci resulting in their misclassification and misinterpretation. In fact, sequence data mapping to repeat elements are often discarded from analysis pipelines. Therefore, there is a continued need for standardized tools and techniques to interpret genomic data of repeats.

Results

We present the UCSC Repeat Browser, which consists of a complete set of human repeat reference sequences derived from annotations made by the commonly used program RepeatMasker. The UCSC Repeat Browser also provides an alignment from the human genome to these references, uses it to map the standard human genome annotation tracks, and presents all of them as a comprehensive interface to facilitate work with repetitive elements. It also provides processed tracks of multiple publicly available datasets of particular interest to the repeat community, including ChIP-seq datasets for KRAB Zinc Finger Proteins (KZNFs) – a family of proteins known to bind and repress certain classes of repeats. We used the UCSC Repeat Browser in combination with these datasets, as well as RepeatMasker annotations in several non-human primates, to trace the independent trajectories of species-specific evolutionary battles between LINE 1 retroelements and their repressors. Furthermore, we document at https://repeatbrowser.ucsc.edu how researchers can map their own human genome annotations to these reference repeat sequences.

Conclusions

The UCSC Repeat Browser allows easy and intuitive visualization of genomic data on consensus repeat elements, circumventing the problem of multi-mapping, in which sequencing reads of repeat elements map to multiple locations on the human genome. By developing a reference consensus, multiple datasets and annotation tracks can easily be overlaid to reveal complex evolutionary histories of repeats in a single interactive window. Specifically, we use this approach to retrace the history of several primate specific LINE-1 families across apes, and discover several species-specific routes of evolution that correlate with the emergence and binding of KZNFs.

Jump around: transposons in and out of the laboratory

Since Barbara McClintock’s groundbreaking discovery of mobile DNA sequences some 70 years ago, transposable elements have come to be recognized as important mutagenic agents impacting genome composition, genome evolution, and human health. Transposable elements are a major constituent of prokaryotic and eukaryotic genomes, and the transposition mechanisms enabling transposon proliferation over evolutionary time remain engaging topics for study, suggesting complex interactions with the host, both antagonistic and mutualistic. The impact of transposition is profound, as over 100 human heritable diseases have been attributed to transposon insertions. Transposition can be highly mutagenic, perturbing genome integrity and gene expression in a wide range of organisms. This mutagenic potential has been exploited in the laboratory, where transposons have long been utilized for phenotypic screening and the generation of defined mutant libraries. More recently, barcoding applications and methods for RNA-directed transposition are being used towards new phenotypic screens and studies relevant for gene therapy. Thus, transposable elements are significant in affecting biology both in vivo and in the laboratory, and this review will survey advances in understanding the biological role of transposons and relevant laboratory applications of these powerful molecular tools.

Dual threshold optimization and network inference reveal convergent evidence from TF binding locations and TF perturbation responses

A high-confidence map of the direct, functional targets of each transcription factor (TF) requires convergent evidence from independent sources. Two significant sources of evidence are TF binding locations and the transcriptional responses to direct TF perturbations. Systematic data sets of both types exist for yeast and human, but they rarely converge on a common set of direct, functional targets for a TF. Even the few genes that are both bound and responsive may not be direct functional targets. Our analysis shows that when there are many nonfunctional binding sites and many indirect targets, nonfunctional sites are expected to occur in the cis-regulatory DNA of indirect targets by chance. To address this problem, we introduce dual threshold optimization (DTO), a new method for setting significance thresholds on binding and perturbation-response data, and show that it improves convergence. It also enables comparison of binding data to perturbation-response data that have been processed by network inference algorithms, which further improves convergence. The combination of dual threshold optimization and network inference greatly expands the high-confidence TF network map in both yeast and human. Next, we analyze a comprehensive new data set measuring the transcriptional response shortly after inducing overexpression of a yeast TF. We also present a new yeast binding location data set obtained by transposon calling cards and compare it to recent ChIP-exo data. These new data sets improve convergence and expand the high-confidence network synergistically.

Widespread correlation of KRAB zinc finger protein binding with brain-developmental gene expression patterns

The large family of KRAB zinc finger (KZNF) genes are transcription factors implicated in recognizing and repressing repetitive sequences such as transposable elements (TEs) in our genome. Through successive waves of retrotransposition-mediated insertions, various classes of TEs have invaded mammalian genomes at multiple timepoints throughout evolution. Even though most of the TE classes in our genome lost the capability to retrotranspose millions of years ago, it remains elusive why the KZNFs that evolved to repress them are still retained in our genome. One hypothesis is that KZNFs become repurposed for other regulatory roles. Here, we find evidence that evolutionary changes in KZNFs provide them not only with the ability to repress TEs, but also to bind to gene promoters independent of TEs. Using KZNF binding site data in conjunction with gene expression values from the Allen Brain Atlas, we show that KZNFs have the ability to regulate gene expression in the human brain in a region-specific manner. Our analysis shows that the expression of KZNFs shows correlation with the expression of their target genes, suggesting that KZNFs have a direct influence on gene expression in the developing human brain. The extent of this regulation and the impact it has on primate brain evolution are still to be determined, but our results imply that KZNFs have become widely integrated into neuronal gene regulatory networks. Our analysis predicts that gene expression networks have been repeatedly innovated throughout primate evolution, continuously gaining new layers of gene regulation mediated by both TEs and KZNFs in our genome. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Tumor edge-to-core transition promotes malignancy in primary-to-recurrent glioblastoma progression in a PLAGL1/CD109-mediated mechanism

BACKGROUND

Glioblastoma remains highly lethal due to its inevitable recurrence. Most of this recurrence is found locally, indicating that postsurgical tumor-initiating cells (TICs) accumulate at the tumor edge. These edge-TICs then generate local recurrence harboring new core lesions. Here, we investigated the clinical significance of the edge-to-core (E-to-C) signature generating glioblastoma recurrence and sought to identify its central mediators.

METHODS

First, we examined the association of E-to-C-related expression changes to patient outcome in matched primary and recurrent samples (n = 37). Specifically, we tested whether the combined decrease of the edge-TIC marker PROM1 (CD133) with the increase of the core-TIC marker CD109, representing E-to-C transition during the primary-to-recurrence progression, indicates poorer patient outcome. We then investigated the specific molecular mediators that trigger tumor recurrence driven by the E-to-C progression. Subsequently, the functional and translational significance of the identified molecule was validated with our patient-derived edge-TIC models in vitro and in vivo.

RESULTS

Patients exhibiting the CD133low/CD109high signature upon recurrence representing E-to-C transition displayed a strong association with poorer progression-free survival and overall survival among all tested patients. Differential gene expression identified that PLAGL1 was tightly correlated with the core TIC marker CD109 and was linked to shorter patient survival. Experimentally, forced PLAGL1 overexpression enhanced, while its knockdown reduced, glioblastoma edge-derived tumor growth in vivo and subsequent mouse survival, suggesting its essential role in the E-to-C-mediated glioblastoma progression.

CONCLUSIONS

E-to-C axis represents an ongoing lethal process in primary glioblastoma contributing to its recurrence, partly in a PLAGL1/CD109-mediated mechanism.

The Arms Race Between KRAB–Zinc Finger Proteins and Endogenous Retroelements and Its Impact on Mammals

Nearly half of the human genome consists of endogenous retroelements (EREs) and their genetic remnants, a small fraction of which carry the potential to propagate in the host genome, posing a threat to genome integrity and cell/organismal survival. The largest family of transcription factors in tetrapods, the Krüppel-associated box domain zinc finger proteins (KRAB-ZFPs), binds to specific EREs and represses their transcription. Since their first appearance over 400 million years ago, KRAB-ZFPs have undergone dramatic expansion and diversification in mammals, correlating with the invasions of new EREs. In this article we review our current understanding of the structure, function, and evolution of KRAB-ZFPs and discuss growing evidence that the arms race between KRAB-ZFPs and the EREs they target is a major driving force for the evolution of new traits in mammals, often accompanied by domestication of EREs themselves.

Small Drosophila zinc finger C2H2 protein with an N-terminal zinc finger-associated domain demonstrates the architecture functions

Recently, the concept has arisen that a special class of architectural proteins exists, which are responsible not only for global chromosome architecture but also for the local regulation of enhancer-promoter interactions. Here, we describe a new architectural protein, with a total size of only 375 aa, which contains an N-terminal zinc finger-associated domain (ZAD) and a cluster of five zinc finger C2H2 domains at the C-terminus. This new protein, named ZAD and Architectural Function 1 protein (ZAF1 protein), is weakly and ubiquitously expressed, with the highest expression levels observed in oocytes and embryos. The cluster of C2H2 domains recognizes a specific 15-bp consensus site, located predominantly in promoters, near transcription start sites. The expression of ZAF1 by a tissue-specific promoter led to the complete blocking of the eye enhancer when clusters of ZAF1 binding sites flanked the eye enhancer in transgenic lines, suggesting that the loop formed by the ZAF1 protein leads to insulation. The ZAF1 protein also supported long-range interactions between the yeast GAL4 activator and the white promoter in transgenic Drosophila lines. A mutant protein lacking the ZAD failed to block the eye enhancer or to support distance interactions in transgenic lines. Taken together, these results suggest that ZAF1 is a minimal architectural protein that can be used to create a convenient model for studying the mechanisms of distance interactions.

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.

Host-transposon interactions: conflict, cooperation, and cooption

Transposable elements (TEs) are mobile DNA sequences that colonize genomes and threaten genome integrity. As a result, several mechanisms appear to have emerged during eukaryotic evolution to suppress TE activity. However, TEs are ubiquitous and account for a prominent fraction of most eukaryotic genomes. We argue that the evolutionary success of TEs cannot be explained solely by evasion from host control mechanisms. Rather, some TEs have evolved commensal and even mutualistic strategies that mitigate the cost of their propagation. These coevolutionary processes promote the emergence of complex cellular activities, which in turn pave the way for cooption of TE sequences for organismal function.

LINE-1 Evasion of Epigenetic Repression in Humans

Epigenetic silencing defends against LINE-1 (L1) retrotransposition in mammalian cells. However, the mechanisms that repress young L1 families and how L1 escapes to cause somatic genome mosaicism in the brain remain unclear. Here we report that a conserved Yin Yang 1 (YY1) transcription factor binding site mediates L1 promoter DNA methylation in pluripotent and differentiated cells. By analyzing 24 hippocampal neurons with three distinct single-cell genomic approaches, we characterized and validated a somatic L1 insertion bearing a 3' transduction. The source (donor) L1 for this insertion was slightly 5' truncated, lacked the YY1 binding site, and was highly mobile when tested in vitro. Locus-specific bisulfite sequencing revealed that the donor L1 and other young L1s with mutated YY1 binding sites were hypomethylated in embryonic stem cells, during neurodifferentiation, and in liver and brain tissue. These results explain how L1 can evade repression and retrotranspose in the human body.

Zscan5b Deficiency Impairs DNA Damage Response and Causes Chromosomal Aberrations during Mitosis

Zygotic genome activation (ZGA) begins after fertilization and is essential for establishing pluripotency and genome stability. However, it is unclear how ZGA genes prevent mitotic errors. Here we show that knockout of the ZGA gene Zscan5b, which encodes a SCAN domain with C2H2 zinc fingers, causes a high incidence of chromosomal abnormalities in embryonic stem cells (ESCs), and leads to the development of early-stage cancers. After irradiation, Zscan5b-deficient ESCs displayed significantly increased levels of γ-H2AX despite increased expression of the DNA repair genes Rad51l3 and Bard. Re-expression of Zscan5b reduced γ-H2AX content, implying a role for Zscan5b in DNA damage repair processes. A co-immunoprecipitation analysis showed that Zscan5b bound to the linker histone H1, suggesting that Zscan5b may protect chromosomal architecture. Our report demonstrates that the ZGA gene Zscan5b is involved in genomic integrity and acts to promote DNA damage repair and regulate chromatin dynamics during mitosis.

•

Deficiency of zygotic genome activation gene Zscan5b causes chromosomal anomalies

•

Zscan5b binds to linker histone H1 and protects chromosomal structure

•

Irradiated Zscan5b-deficient ESCs show significantly increased DNA stress markers

•

Zscan5b-deficient ESCs develop small choriocarcinomas and embryonal carcinomas

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.

ZNF341 controls STAT3 expression and thereby immunocompetence

Signal transducer and activator of transcription 3 (STAT3) is a central regulator of immune homeostasis. STAT3 levels are strictly controlled, and STAT3 impairment contributes to several diseases including the monogenic autosomal-dominant hyper-immunoglobulin E (IgE) syndrome (AD-HIES). We investigated patients of four consanguineous families with an autosomal-recessive disorder resembling the phenotype of AD-HIES, with symptoms of immunodeficiency, recurrent infections, skeletal abnormalities, and elevated IgE. Patients presented with reduced STAT3 expression and diminished T helper 17 cell numbers, in absence of STAT3 mutations. We identified two distinct homozygous nonsense mutations in ZNF341, which encodes a zinc finger transcription factor. Wild-type ZNF341 bound to and activated the STAT3 promoter, whereas the mutant variants showed impaired transcriptional activation, partly due to nuclear translocation failure. In summary, nonsense mutations in ZNF341 account for the STAT3-like phenotype in four autosomal-recessive kindreds. Thus, ZNF341 is a previously unrecognized regulator of immune homeostasis.

Genome-wide CRISPR Screens in T Helper Cells Reveal Pervasive Crosstalk between Activation and Differentiation

T helper type 2 (Th2) cells are important regulators of mammalian adaptive immunity and have relevance for infection, autoimmunity, and tumor immunology. Using a newly developed, genome-wide retroviral CRISPR knockout (KO) library, combined with RNA-seq, ATAC-seq, and ChIP-seq, we have dissected the regulatory circuitry governing activation and differentiation of these cells. Our experiments distinguish cell activation versus differentiation in a quantitative framework. We demonstrate that these two processes are tightly coupled and are jointly controlled by many transcription factors, metabolic genes, and cytokine/receptor pairs. There are only a small number of genes regulating differentiation without any role in activation. By combining biochemical and genetic data, we provide an atlas for Th2 differentiation, validating known regulators and identifying factors, such as Pparg and Bhlhe40, as part of the core regulatory network governing Th2 helper cell fates.

The Human Transcription Factors

Transcription factors (TFs) recognize specific DNA sequences to control chromatin and transcription, forming a complex system that guides expression of the genome. Despite keen interest in understanding how TFs control gene expression, it remains challenging to determine how the precise genomic binding sites of TFs are specified and how TF binding ultimately relates to regulation of transcription. This review considers how TFs are identified and functionally characterized, principally through the lens of a catalog of over 1,600 likely human TFs and binding motifs for two-thirds of them. Major classes of human TFs differ markedly in their evolutionary trajectories and expression patterns, underscoring distinct functions. TFs likewise underlie many different aspects of human physiology, disease, and variation, highlighting the importance of continued effort to understand TF-mediated gene regulation.

Cistrome Data Browser: expanded datasets and new tools for gene regulatory analysis

The Cistrome Data Browser (DB) is a resource of human and mouse cis-regulatory information derived from ChIP-seq, DNase-seq and ATAC-seq chromatin profiling assays, which map the genome-wide locations of transcription factor binding sites, histone post-translational modifications and regions of chromatin accessible to endonuclease activity. Currently, the Cistrome DB contains approximately 47,000 human and mouse samples with about 24,000 newly collected datasets compared to the previous release two years ago. Furthermore, the Cistrome DB has a new Toolkit module with several features that allow users to better utilize the large-scale ChIP-seq, DNase-seq, and ATAC-seq data. First, users can query the factors which are likely to regulate a specific gene of interest. Second, the Cistrome DB Toolkit facilitates searches for factor binding, histone modifications, and chromatin accessibility in any given genomic interval shorter than 2Mb. Third, the Toolkit can determine the most similar ChIP-seq, DNase-seq, and ATAC-seq samples in terms of genomic interval overlaps with user-provided genomic interval sets. The Cistrome DB is a user-friendly, up-to-date, and well maintained resource, and the new tools will greatly benefit the biomedical research community. The database is freely available at http://cistrome.org/db, and the Toolkit is at http://dbtoolkit.cistrome.org.

Non-redundant functions of EMT transcription factors

Epithelial-mesenchymal transition (EMT) is a crucial embryonic programme that is executed by various EMT transcription factors (EMT-TFs) and is aberrantly activated in cancer and other diseases. However, the causal role of EMT and EMT-TFs in different disease processes, especially cancer and metastasis, continues to be debated. In this Review, we identify and describe specific, non-redundant functions of the different EMT-TFs and discuss the reasons that may underlie disputes about EMT in cancer.

Sequences encoding C2H2 zinc fingers inhibit polyadenylation and mRNA export in human cells

The large C2H2-Zinc Finger (C2H2-ZNF) gene family has rapidly expanded in primates through gene duplication. There is consequently considerable sequence homology between family members at both the nucleotide and amino acid level, allowing for coordinated regulation and shared functions. Here we show that multiple C2H2-ZNF mRNAs experience differential polyadenylation resulting in populations with short and long poly(A) tails. Furthermore, a significant proportion of C2H2-ZNF mRNAs are retained in the nucleus. Intriguingly, both short poly(A) tails and nuclear retention can be specified by the repeated elements that encode zinc finger motifs. These Zinc finger Coding Regions (ZCRs) appear to restrict polyadenylation of nascent RNAs and at the same time impede their export. However, the polyadenylation process is not necessary for nuclear retention of ZNF mRNAs. We propose that inefficient polyadenylation and export may allow C2H2-ZNF mRNAs to moonlight as non-coding RNAs or to be stored for later use.

Sequence Determinants for Nuclear Retention and Cytoplasmic Export of mRNAs and lncRNAs

Eukaryotes are divided into two major compartments: the nucleus where RNA is synthesized and processed, and the cytoplasm, where mRNA is translated into proteins. Although many different RNAs are made, only a subset is allowed access to the cytoplasm, primarily RNAs involved in protein synthesis (mRNA, tRNA, and rRNA). In contrast, nuclear retained transcripts are mostly long non-coding RNAs (lncRNAs) whose role in cell physiology has been a source of much investigation in the past few years. In addition, it is likely that many non-functional RNAs, which arise by spurious transcription and misprocessing of functional RNAs, are also retained in the nucleus and degraded. In this review, the main sequence features that dictate whether any particular mRNA or lncRNA is a substrate for retention in the nucleus, or export to the cytoplasm, are discussed. Although nuclear export is promoted by RNA-splicing due to the fact that the spliceosome can help recruit export factors to the mature RNA, nuclear export does not require splicing. Indeed, most stable unspliced transcripts are well exported and associate with these same export factors in a splicing-independent manner. In contrast, nuclear retention is promoted by specialized cis-elements found in certain RNAs. This new understanding of the determinants of nuclear retention and cytoplasmic export provides a deeper understanding of how information flow is regulated in eukaryotic cells. Ultimately these processes promote the evolution of complexity in eukaryotes by shaping the genomic content through constructive neutral evolution.

FEZF1 is an Independent Predictive Factor for Recurrence and Promotes Cell Proliferation and Migration in Cervical Cancer

The Fez family zinc finger protein 1 (FEZF1), a critical transcription factor in nervous system development, has been implicated in cancer progression recently. However, its clinical significance remains unknown. By analyzing gene expression data of eight most common cancer types from The Cancer Genome Atlas (TCGA), we found that FEZF1 prominently associated with the recurrence-free survival of cervical cancer patients (P<0.001) and was an independent diagnostic factor for cervical cancer recurrence (P=0.002). Moreover, FEZF1 expression was significantly higher in the tumor samples from cervical cancer patients with relapse in TCGA(P=0.015). By RNA interference, we knocked down FEZF1 and found that cell proliferation, growth and migration were significantly decreased in C33A and SiHa cells. Meanwhile, FEZF1 knockdown also attenuated the growth of C33A cells in nude mice. In contrast, expression of FEZF1 promoted cell proliferation, growth and migration in HeLa cells. Using chromatin immunoprecipitation (ChIP) assay, we revealed that FEZF1 could bind to multiple key genes in the Wnt signaling pathway in HeLa cells. Furthermore, analysis of the levels of β-catenin protein, the core component of the Wnt pathway, and downstream effector genes of the pathway showed that FEZF1 could activate the Wnt pathway. Together, these results suggest that FEZF1 promotes cell proliferation and migration possibly by acting as a transcriptional activator of the Wnt signaling pathway in cervical cancer, and also provide a valuable molecular predictive marker for cervical cancer recurrence.

KAP1 regulates endogenous retroviruses in adult human cells and contributes to innate immune control

Endogenous retroviruses (ERVs) have accumulated in vertebrate genomes and contribute to the complexity of gene regulation. KAP1 represses ERVs during development by its recruitment to their repetitive sequences through KRAB zinc-finger proteins (KZNFs), but little is known about the regulation of ERVs in adult tissues. We observed that KAP1 repression of HERVK14C was conserved in differentiated human cells and performed KAP1 knockout to obtain an overview of KAP1 function. Our results show that KAP1 represses ERVs (including HERV-T and HERV-S) and ZNF genes, both of which overlap with KAP1 binding sites and H3K9me3 in multiple cell types. Furthermore, this pathway is functionally conserved in adult human peripheral blood mononuclear cells. Cytosine methylation that acts on KAP1 regulated loci is necessary to prevent an interferon response, and KAP1-depletion leads to activation of some interferon-stimulated genes. Finally, loss of KAP1 leads to a decrease in H3K9me3 enrichment at ERVs and ZNF genes and an RNA-sensing response mediated through MAVS signaling. These data indicate that the KAP1-KZNF pathway contributes to genome stability and innate immune control in adult human cells.

Human cardiac cis-regulatory elements, their cognate transcription factors, and regulatory DNA sequence variants

Cis-regulatory elements (CRE), short DNA sequences through which transcription factors (TFs) exert regulatory control on gene expression, are postulated to be the major sites of causal sequence variation underlying the genetics of complex traits and diseases. We present integrative analyses, combining high-throughput genomic and epigenomic data with sequence-based computations, to identify the causal transcriptional components in a given tissue. We use data on adult human hearts to demonstrate that (1) sequence-based predictions detect numerous, active, tissue-specific CREs missed by experimental observations, (2) learned sequence features identify the cognate TFs, (3) CRE variants are specifically associated with cardiac gene expression, and (4) a significant fraction of the heritability of exemplar cardiac traits (QT interval, blood pressure, pulse rate) is attributable to these variants. This general systems approach can thus identify candidate causal variants and the components of gene regulatory networks (GRN) to enable understanding of the mechanisms of complex disorders on a tissue- or cell-type basis.

Thalidomide promotes degradation of SALL4, a transcription factor implicated in Duane Radial Ray syndrome

In historical attempts to treat morning sickness, use of the drug thalidomide led to the birth of thousands of children with severe birth defects. Despite their teratogenicity, thalidomide and related IMiD drugs are now a mainstay of cancer treatment; however, the molecular basis underlying the pleiotropic biology and characteristic birth defects remains unknown. Here we show that IMiDs disrupt a broad transcriptional network through induced degradation of several C₂H₂ zinc finger transcription factors, including SALL4, a member of the spalt-like family of developmental transcription factors. Strikingly, heterozygous loss of function mutations in SALL4 result in a human developmental condition that phenocopies thalidomide-induced birth defects such as absence of thumbs, phocomelia, defects in ear and eye development, and congenital heart disease. We find that thalidomide induces degradation of SALL4 exclusively in humans, primates, and rabbits, but not in rodents or fish, providing a mechanistic link for the species-specific pathogenesis of thalidomide syndrome.

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.

The HUSH complex cooperates with TRIM28 to repress young retrotransposons and new genes

Retrotransposons encompass half of the human genome and contribute to the formation of heterochromatin, which provides nuclear structure and regulates gene expression. Here, we asked if the human silencing hub (HUSH) complex is necessary to silence retrotransposons and whether it collaborates with TRIM28 and the chromatin remodeler ATRX at specific genomic loci. We show that the HUSH complex contributes to de novo repression and DNA methylation of an SVA retrotransposon reporter. By using naïve versus primed mouse pluripotent stem cells, we reveal a critical role for the HUSH complex in naïve cells, implicating it in programming epigenetic marks in development. Although the HUSH component FAM208A binds to endogenous retroviruses (ERVs) and long interspersed element-1s (LINE-1s or L1s), it is mainly required to repress evolutionarily young L1s (mouse-specific lineages <5 million years old). TRIM28, in contrast, is necessary to repress both ERVs and young L1s. Genes co-repressed by TRIM28 and FAM208A are evolutionarily young, or exhibit tissue-specific expression, are enriched in young L1s, and display evidence for regulation through LTR promoters. Finally, we demonstrate that the HUSH complex is also required to repress L1 elements in human cells. Overall, these data indicate that the HUSH complex and TRIM28 co-repress young retrotransposons and new genes rewired by retrotransposon noncoding DNA.

Silencing of transposable elements may not be a major driver of regulatory evolution in primate iPSCs

Transposable elements (TEs) comprise almost half of primate genomes and their aberrant regulation can result in deleterious effects. In pluripotent stem cells, rapidly evolving KRAB-ZNF genes target TEs for silencing by H3K9me3. To investigate the evolution of TE silencing, we performed H3K9me3 ChIP-seq experiments in induced pluripotent stem cells from 10 human and 7 chimpanzee individuals. We identified four million orthologous TEs and found the SVA and ERV families to be marked most frequently by H3K9me3. We found little evidence of inter-species differences in TE silencing, with as many as 82% of putatively silenced TEs marked at similar levels in humans and chimpanzees. TEs that are preferentially silenced in one species are a similar age to those silenced in both species and are not more likely to be associated with expression divergence of nearby orthologous genes. Our data suggest limited species-specificity of TE silencing across 6 million years of primate evolution.

Zinc-finger protein 471 suppresses gastric cancer through transcriptionally repressing downstream oncogenic PLS3 and TFAP2A

Zinc-finger protein 471 (ZNF471) was preferentially methylated in gastric cancer using promoter methylation array. The role of ZNF471 in human cancer is unclear. Here we elucidated the functional significance, molecular mechanisms and clinical impact of ZNF471 in gastric cancer. ZNF471 mRNA was silenced in 15 out of 16 gastric cancer cell lines due to promoter hypermethylation. Significantly higher ZNF471 promoter methylation was also observed in primary gastric cancers compared to their adjacent normal tissues (P < 0.001). ZNF471 promoter CpG-site hypermethylation correlated with poor survival of gastric cancer patients (n = 120, P = 0.001). Ectopic expression of ZNF471 in gastric cancer cell lines (AGS, BGC823, and MKN74) significantly suppressed cell proliferation, migration, and invasion, while it induced apoptosis in vitro and inhibited xenograft tumorigenesis in nude mice. Transcription factor AP-2 Alpha (TFAP2A) and plastin3 (PLS3) were two crucial downstream targets of ZNF471 demonstrated by bioinformatics modeling and ChIP-PCR assays. ZNF471 directly bound to the promoter of TFAP2A and PLS3 and transcriptionally inhibited their expression. TFAP2A and PLS3 showed oncogenic functions in gastric cancer cell lines. Moreover, ZNF471 recruited KAP1 to the promoter of the target genes, thereby inducing H3K9me3 enrichment for transcriptional repression and inhibition of oncogenic TFAP2A and PLS3. In conclusion, ZNF471 acts as a tumor suppressor in gastric cancer by transcriptionally inhibiting downstream targets TFAP2A and PLS3. KAP1 is a co-repressor of ZNF471 at the promoter of the target genes. The promoter CpG-site methylation is an independent prognostic factor for overall survival of gastric cancer patients.