High-Throughput Discovery and Characterization of Human Transcriptional Effectors

Systematic identification of transcriptional activation domains from non-transcription factor proteins in plants and yeast

Transcription factors can promote gene expression through activation domains. Whole-genome screens have systematically mapped activation domains in transcription factors but not in non-transcription factor proteins (e.g., chromatin regulators and coactivators). To fill this knowledge gap, we employed the activation domain predictor PADDLE to analyze the proteomes of Arabidopsis thaliana and Saccharomyces cerevisiae. We screened 18,000 predicted activation domains from >800 non-transcription factor genes in both species, confirming that 89% of candidate proteins contain active fragments. Our work enables the annotation of hundreds of nuclear proteins as putative coactivators, many of which have never been ascribed any function in plants. Analysis of peptide sequence compositions reveals how the distribution of key amino acids dictates activity. Finally, we validated short, "universal" activation domains with comparable performance to state-of-the-art activation domains used for genome engineering. Our approach enables the genome-wide discovery and annotation of activation domains that can function across diverse eukaryotes.

CRISPR technologies for genome, epigenome and transcriptome editing

Our ability to edit genomes lags behind our capacity to sequence them, but the growing understanding of CRISPR biology and its application to genome, epigenome and transcriptome engineering is narrowing this gap. In this Review, we discuss recent developments of various CRISPR-based systems that can transiently or permanently modify the genome and the transcriptome. The discovery of further CRISPR enzymes and systems through functional metagenomics has meaningfully broadened the applicability of CRISPR-based editing. Engineered Cas variants offer diverse capabilities such as base editing, prime editing, gene insertion and gene regulation, thereby providing a panoply of tools for the scientific community. We highlight the strengths and weaknesses of current CRISPR tools, considering their efficiency, precision, specificity, reliance on cellular DNA repair mechanisms and their applications in both fundamental biology and therapeutics. Finally, we discuss ongoing clinical trials that illustrate the potential impact of CRISPR systems on human health.

Post-Transcriptional Modular Synthetic Receptors

Inspired by the power of transcriptional synthetic receptors and hoping to complement them to expand the toolbox for cell engineering, we establish LIDAR (Ligand-Induced Dimerization Activating RNA editing), a modular post-transcriptional synthetic receptor platform that harnesses RNA editing by ADAR. LIDAR is compatible with various receptor architectures in different cellular contexts, and enables the sensing of diverse ligands and the production of functional outputs. Furthermore, LIDAR can sense orthogonal signals in the same cell and produce synthetic spatial patterns, potentially enabling the programming of complex multicellular behaviors. Finally, LIDAR is compatible with compact encoding and can be delivered by synthetic mRNA. Thus, LIDAR expands the family of synthetic receptors, holding the promise to empower basic research and therapeutic applications.

The transcription factor ZNF469 regulates collagen production in liver fibrosis

Non-alcoholic fatty liver disease (NAFLD) - characterized by excess accumulation of fat in the liver - now affects one third of the world's population. As NAFLD progresses, extracellular matrix components including collagen accumulate in the liver causing tissue fibrosis, a major determinant of disease severity and mortality. To identify transcriptional regulators of fibrosis, we computationally inferred the activity of transcription factors (TFs) relevant to fibrosis by profiling the matched transcriptomes and epigenomes of 108 human liver biopsies from a deeply-characterized cohort of patients spanning the full histopathologic spectrum of NAFLD. CRISPR-based genetic knockout of the top 100 TFs identified ZNF469 as a regulator of collagen expression in primary human hepatic stellate cells (HSCs). Gain- and loss-of-function studies established that ZNF469 regulates collagen genes and genes involved in matrix homeostasis through direct binding to gene bodies and regulatory elements. By integrating multiomic large-scale profiling of human biopsies with extensive experimental validation we demonstrate that ZNF469 is a transcriptional regulator of collagen in HSCs. Overall, these data nominate ZNF469 as a previously unrecognized determinant of NAFLD-associated liver fibrosis.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Improving prime editing with an endogenous small RNA-binding protein

Prime editing enables the precise modification of genomes through reverse transcription of template sequences appended to the 3' ends of CRISPR-Cas guide RNAs1. To identify cellular determinants of prime editing, we developed scalable prime editing reporters and performed genome-scale CRISPR-interference screens. From these screens, a single factor emerged as the strongest mediator of prime editing: the small RNA-binding exonuclease protection factor La. Further investigation revealed that La promotes prime editing across approaches (PE2, PE3, PE4 and PE5), edit types (substitutions, insertions and deletions), endogenous loci and cell types but has no consistent effect on genome-editing approaches that rely on standard, unextended guide RNAs. Previous work has shown that La binds polyuridine tracts at the 3' ends of RNA polymerase III transcripts2. We found that La functionally interacts with the 3' ends of polyuridylated prime editing guide RNAs (pegRNAs). Guided by these results, we developed a prime editor protein (PE7) fused to the RNA-binding, N-terminal domain of La. This editor improved prime editing with expressed pegRNAs and engineered pegRNAs (epegRNAs), as well as with synthetic pegRNAs optimized for La binding. Together, our results provide key insights into how prime editing components interact with the cellular environment and suggest general strategies for stabilizing exogenous small RNAs therein.

Defining the mechanisms and properties of post-transcriptional regulatory disordered regions by high-throughput functional profiling

Disordered regions within RNA binding proteins are required to control mRNA decay and protein synthesis. To understand how these disordered regions modulate gene expression, we surveyed regulatory activity across the entire disordered proteome using a high-throughput functional assay. We identified hundreds of regulatory sequences within intrinsically disordered regions and demonstrate how these elements cooperate with core mRNA decay machinery to promote transcript turnover. Coupling high-throughput functional profiling with mutational scanning revealed diverse molecular features, ranging from defined motifs to overall sequence composition, underlying the regulatory effects of disordered peptides. Machine learning analysis implicated aromatic residues in particular contexts as critical determinants of repressor activity, consistent with their roles in forming protein-protein interactions with downstream effectors. Our results define the molecular principles and biochemical mechanisms that govern post-transcriptional gene regulation by disordered regions and exemplify the encoding of diverse yet specific functions in the absence of well-defined structure.

Single-molecule chromatin configurations link transcription factor binding to expression in human cells

The binding of multiple transcription factors (TFs) to genomic enhancers activates gene expression in mammalian cells. However, the molecular details that link enhancer sequence to TF binding, promoter state, and gene expression levels remain opaque. We applied single-molecule footprinting (SMF) to measure the simultaneous occupancy of TFs, nucleosomes, and components of the transcription machinery on engineered enhancer/promoter constructs with variable numbers of TF binding sites for both a synthetic and an endogenous TF. We find that activation domains enhance a TF's capacity to compete with nucleosomes for binding to DNA in a BAF-dependent manner, TF binding on nucleosome-free DNA is consistent with independent binding between TFs, and average TF occupancy linearly contributes to promoter activation rates. We also decompose TF strength into separable binding and activation terms, which can be tuned and perturbed independently. Finally, we develop thermodynamic and kinetic models that quantitatively predict both the binding microstates observed at the enhancer and subsequent time-dependent gene expression. This work provides a template for quantitative dissection of distinct contributors to gene activation, including the activity of chromatin remodelers, TF activation domains, chromatin acetylation, TF concentration, TF binding affinity, and TF binding site configuration.

A modular dCas9-based recruitment platform for combinatorial epigenome editing

Targeted epigenome editing tools allow precise manipulation and investigation of genome modifications, however they often display high context dependency and variable efficacy between target genes and cell types. While systems that simultaneously recruit multiple distinct 'effector' chromatin regulators can improve efficacy, they generally lack control over effector composition and spatial organisation. To overcome this we have created a modular combinatorial epigenome editing platform, called SSSavi. This system is an interchangeable and reconfigurable docking platform fused to dCas9 that enables simultaneous recruitment of up to four different effectors, allowing precise control of effector composition and spatial ordering. We demonstrate the activity and specificity of the SSSavi system and, by testing it against existing multi-effector targeting systems, demonstrate its comparable efficacy. Furthermore, we demonstrate the importance of the spatial ordering of the recruited effectors for effective transcriptional regulation. Together, the SSSavi system enables exploration of combinatorial effector co-recruitment to enhance manipulation of chromatin contexts previously resistant to targeted editing.

Large-scale evaluation of the ability of RNA-binding proteins to activate exon inclusion

RNA-binding proteins (RBPs) modulate alternative splicing outcomes to determine isoform expression and cellular survival. To identify RBPs that directly drive alternative exon inclusion, we developed tethered function luciferase-based splicing reporters that provide rapid, scalable and robust readouts of exon inclusion changes and used these to evaluate 718 human RBPs. We performed enhanced cross-linking immunoprecipitation, RNA sequencing and affinity purification-mass spectrometry to investigate a subset of candidates with no prior association with splicing. Integrative analysis of these assays indicates surprising roles for TRNAU1AP, SCAF8 and RTCA in the modulation of hundreds of endogenous splicing events. We also leveraged our tethering assays and top candidates to identify potent and compact exon inclusion activation domains for splicing modulation applications. Using these identified domains, we engineered programmable fusion proteins that outperform current artificial splicing factors at manipulating inclusion of reporter and endogenous exons. This tethering approach characterizes the ability of RBPs to induce exon inclusion and yields new molecular parts for programmable splicing control.

Hold out the genome: a roadmap to solving the cis-regulatory code

Gene expression is regulated by transcription factors that work together to read cis-regulatory DNA sequences. The 'cis-regulatory code' - how cells interpret DNA sequences to determine when, where and how much genes should be expressed - has proven to be exceedingly complex. Recently, advances in the scale and resolution of functional genomics assays and machine learning have enabled substantial progress towards deciphering this code. However, the cis-regulatory code will probably never be solved if models are trained only on genomic sequences; regions of homology can easily lead to overestimation of predictive performance, and our genome is too short and has insufficient sequence diversity to learn all relevant parameters. Fortunately, randomly synthesized DNA sequences enable testing a far larger sequence space than exists in our genomes, and designed DNA sequences enable targeted queries to maximally improve the models. As the same biochemical principles are used to interpret DNA regardless of its source, models trained on these synthetic data can predict genomic activity, often better than genome-trained models. Here we provide an outlook on the field, and propose a roadmap towards solving the cis-regulatory code by a combination of machine learning and massively parallel assays using synthetic DNA.

Deep mutational scanning of proteins in mammalian cells

Protein mutagenesis is essential for unveiling the molecular mechanisms underlying protein function in health, disease, and evolution. In the past decade, deep mutational scanning methods have evolved to support the functional analysis of nearly all possible single-amino acid changes in a protein of interest. While historically these methods were developed in lower organisms such as E. coli and yeast, recent technological advancements have resulted in the increased use of mammalian cells, particularly for studying proteins involved in human disease. These advancements will aid significantly in the classification and interpretation of variants of unknown significance, which are being discovered at large scale due to the current surge in the use of whole-genome sequencing in clinical contexts. Here, we explore the experimental aspects of deep mutational scanning studies in mammalian cells and report the different methods used in each step of the workflow, ultimately providing a useful guide toward the design of such studies.

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.

Compact engineered human mechanosensitive transactivation modules enable potent and versatile synthetic transcriptional control

Engineered transactivation domains (TADs) combined with programmable DNA binding platforms have revolutionized synthetic transcriptional control. Despite recent progress in programmable CRISPR-Cas-based transactivation (CRISPRa) technologies, the TADs used in these systems often contain poorly tolerated elements and/or are prohibitively large for many applications. Here, we defined and optimized minimal TADs built from human mechanosensitive transcription factors. We used these components to construct potent and compact multipartite transactivation modules (MSN, NMS and eN3x9) and to build the CRISPR-dCas9 recruited enhanced activation module (CRISPR-DREAM) platform. We found that CRISPR-DREAM was specific and robust across mammalian cell types, and efficiently stimulated transcription from diverse regulatory loci. We also showed that MSN and NMS were portable across Type I, II and V CRISPR systems, transcription activator-like effectors and zinc finger proteins. Further, as proofs of concept, we used dCas9-NMS to efficiently reprogram human fibroblasts into induced pluripotent stem cells and demonstrated that mechanosensitive transcription factor TADs are efficacious and well tolerated in therapeutically important primary human cell types. Finally, we leveraged the compact and potent features of these engineered TADs to build dual and all-in-one CRISPRa AAV systems. Altogether, these compact human TADs, fusion modules and delivery architectures should be valuable for synthetic transcriptional control in biomedical applications.

Engineering custom morpho- and chemotypes of Populus for sustainable production of biofuels, bioproducts, and biomaterials

Humans have been modifying plant traits for thousands of years, first through selection (i.e., domestication) then modern breeding, and in the last 30 years, through biotechnology. These modifications have resulted in increased yield, more efficient agronomic practices, and enhanced quality traits. Precision knowledge of gene regulation and function through high-resolution single-cell omics technologies, coupled with the ability to engineer plant genomes at the DNA sequence, chromatin accessibility, and gene expression levels, can enable engineering of complex and complementary traits at the biosystem level. Populus spp., the primary genetic model system for woody perennials, are among the fastest growing trees in temperate zones and are important for both carbon sequestration and global carbon cycling. Ample genomic and transcriptomic resources for poplar are available including emerging single-cell omics datasets. To expand use of poplar outside of valorization of woody biomass, chassis with novel morphotypes in which stem branching and tree height are modified can be fabricated thereby leading to trees with altered leaf to wood ratios. These morphotypes can then be engineered into customized chemotypes that produce high value biofuels, bioproducts, and biomaterials not only in specific organs but also in a cell-type-specific manner. For example, the recent discovery of triterpene production in poplar leaf trichomes can be exploited using cell-type specific regulatory sequences to synthesize high value terpenes such as the jet fuel precursor bisabolene specifically in the trichomes. By spatially and temporally controlling expression, not only can pools of abundant precursors be exploited but engineered molecules can be sequestered in discrete cell structures in the leaf. The structural diversity of the hemicellulose xylan is a barrier to fully utilizing lignocellulose in biomaterial production and by leveraging cell-type-specific omics data, cell wall composition can be modified in a tailored and targeted specific manner to generate poplar wood with novel chemical features that are amenable for processing or advanced manufacturing. Precision engineering poplar as a multi-purpose sustainable feedstock highlights how genome engineering can be used to re-imagine a crop species.

The H3.3 K36M oncohistone disrupts the establishment of epigenetic memory through loss of DNA methylation

Histone H3.3 is frequently mutated in cancers, with the lysine 36 to methionine mutation (K36M) being a hallmark of chondroblastomas. While it is known that H3.3K36M changes the cellular epigenetic landscape, it remains unclear how it affects the dynamics of gene expression. Here, we use a synthetic reporter to measure the effect of H3.3K36M on silencing and epigenetic memory after recruitment of KRAB: a member of the largest class of human repressors, commonly used in synthetic biology, and associated with H3K9me3. We find that H3.3K36M, which decreases H3K36 methylation, leads to a decrease in epigenetic memory and promoter methylation weeks after KRAB release. We propose a new model for establishment and maintenance of epigenetic memory, where H3K36 methylation is necessary to convert H3K9me3 domains into DNA methylation for stable epigenetic memory. Our quantitative model can inform oncogenic mechanisms and guide development of epigenetic editing tools.

Single-cell chromatin state transitions during epigenetic memory formation

Repressive chromatin modifications are thought to compact chromatin to silence transcription. However, it is unclear how chromatin structure changes during silencing and epigenetic memory formation. We measured gene expression and chromatin structure in single cells after recruitment and release of repressors at a reporter gene. Chromatin structure is heterogeneous, with open and compact conformations present in both active and silent states. Recruitment of repressors associated with epigenetic memory produces chromatin compaction across 10-20 kilobases, while reversible silencing does not cause compaction at this scale. Chromatin compaction is inherited, but changes molecularly over time from histone methylation (H3K9me3) to DNA methylation. The level of compaction at the end of silencing quantitatively predicts epigenetic memory weeks later. Similarly, chromatin compaction at the Nanog locus predicts the degree of stem-cell fate commitment. These findings suggest that the chromatin state across tens of kilobases, beyond the gene itself, is important for epigenetic memory formation.

Clusters of acidic and hydrophobic residues can predict acidic transcriptional activation domains from protein sequence

Transcription factors activate gene expression in development, homeostasis, and stress with DNA binding domains and activation domains. Although there exist excellent computational models for predicting DNA binding domains from protein sequence, models for predicting activation domains from protein sequence have lagged, particularly in metazoans. We recently developed a simple and accurate predictor of acidic activation domains on human transcription factors. Here, we show how the accuracy of this human predictor arises from the clustering of aromatic, leucine, and acidic residues, which together are necessary for acidic activation domain function. When we combine our predictor with the predictions of convolutional neural network (CNN) models trained in yeast, the intersection is more accurate than individual models, emphasizing that each approach carries orthogonal information. We synthesize these findings into a new set of activation domain predictions on human transcription factors.

Compact Programmable Control of Protein Secretion in Mammalian Cells

Synthetic biology currently holds immense potential to engineer the spatiotemporal control of intercellular signals for biomedicine. Programming behaviors using protein-based circuits has advantages over traditional gene circuits such as compact delivery and direct interactions with signaling proteins. Previously, we described a generalizable platform called RELEASE to enable the control of intercellular signaling through the proteolytic removal of ER-retention motifs compatible with pre-existing protease-based circuits. However, these tools lacked the ability to reliably program complex expression profiles and required numerous proteases, limiting delivery options. Here, we harness the recruitment and antagonistic behavior of endogenous 14-3-3 proteins to create RELEASE-NOT to turn off protein secretion in response to protease activity. By combining RELEASE and RELEASE-NOT, we establish a suite of protein-level processing and output modules called Compact RELEASE (compRELEASE). This innovation enables functions such as logic processing and analog signal filtering using a single input protease. Furthermore, we demonstrate the compactness of the post-translational design by using polycistronic single transcripts to engineer cells to control protein secretion via lentiviral integration and leverage mRNA delivery to selectively express cell surface proteins only in engineered cells harboring inducible proteases. CompRELEASE enables complex control of protein secretion and enhances the potential of synthetic protein circuits for therapeutic applications, while minimizing the overall genetic payload.

Frontiers and techniques in plant gene regulation

Understanding plant gene regulation has been a priority for generations of plant scientists. However, due to its complex nature, the regulatory code governing plant gene expression has yet to be deciphered comprehensively. Recently developed methods-often relying on next-generation sequencing technology and state-of-the-art computational approaches-have started to further our understanding of the gene regulatory logic used by plants. In this review, we discuss these methods and the insights into the regulatory code of plants that they can yield.

High-throughput functional characterization of combinations of transcriptional activators and repressors

Despite growing knowledge of the functions of individual human transcriptional effector domains, much less is understood about how multiple effector domains within the same protein combine to regulate gene expression. Here, we measure transcriptional activity for 8,400 effector domain combinations by recruiting them to reporter genes in human cells. In our assay, weak and moderate activation domains synergize to drive strong gene expression, whereas combining strong activators often results in weaker activation. In contrast, repressors combine linearly and produce full gene silencing, and repressor domains often overpower activation domains. We use this information to build a synthetic transcription factor whose function can be tuned between repression and activation independent of recruitment to target genes by using a small-molecule drug. Altogether, we outline the basic principles of how effector domains combine to regulate gene expression and demonstrate their value in building precise and flexible synthetic biology tools. A record of this paper's transparent peer review process is included in the supplemental information.

Optimization of Cas12a for multiplexed genome-scale transcriptional activation

Cas12a CRISPR technology, unlike Cas9, allows for facile multiplexing of guide RNAs from a single transcript, simplifying combinatorial perturbations. While Cas12a has been implemented for multiplexed knockout genetic screens, it has yet to be optimized for CRISPR activation (CRISPRa) screens in human cells. Here, we develop a new Cas12a-based transactivation domain (TAD) recruitment system using the ALFA nanobody and demonstrate simultaneous activation of up to four genes. We screen a genome-wide library to identify modulators of growth and MEK inhibition, and we compare these results with those obtained with open reading frame (ORF) overexpression and Cas9-based CRISPRa. We find that the activity of multiplexed arrays is largely predictable from the best-performing guide and provide criteria for selecting active guides. We anticipate that these results will greatly accelerate the exploration of gene function and combinatorial phenotypes at scale.

Pooled endogenous protein tagging and recruitment for scalable discovery of effectors for induced proximity therapeutics

The field of induced proximity therapeutics is in its ascendancy but is limited by a lack of scalable tools to systematically explore effector-target protein pairs in an unbiased manner. Here, we combined Scalable POoled Targeting with a LIgandable Tag at Endogenous Sites (SPOTLITES) for the high-throughput tagging of endogenous proteins, with generic small molecule-based protein recruitment to screen for novel proximity-based effectors. We apply this methodology in two orthogonal screens for targeted protein degradation: the first using fluorescence to monitor target protein levels directly, and the second using a cellular growth phenotype that depends on the degradation of an essential protein. Our screens revealed a multitude of potential new effector proteins for degradation and converged on members of the CTLH complex which we demonstrate potently induce degradation. Altogether, we introduce a platform for pooled induction of endogenous protein-protein interactions that can be used to expand our toolset of effector proteins for targeted protein degradation and other forms of induced proximity.

Systematic identification of transcriptional activator domains from non-transcription factor proteins in plants and yeast

Transcription factors promote gene expression via trans-regulatory activation domains. Although whole genome scale screens in model organisms (e.g. human, yeast, fly) have helped identify activation domains from transcription factors, such screens have been less extensively used to explore the occurrence of activation domains in non-transcription factor proteins, such as transcriptional coactivators, chromatin regulators and some cytosolic proteins, leaving a blind spot on what role activation domains in these proteins could play in regulating transcription. We utilized the activation domain predictor PADDLE to mine the entire proteomes of two model eukaryotes, Arabidopsis thaliana and Saccharomyces cerevisiae ( 1 ). We characterized 18,000 fragments covering predicted activation domains from >800 non-transcription factor genes in both species, and experimentally validated that 89% of proteins contained fragments capable of activating transcription in yeast. Peptides with similar sequence composition show a broad range of activities, which is explained by the arrangement of key amino acids. We also annotated hundreds of nuclear proteins with activation domains as putative coactivators; many of which have never been ascribed any function in plants. Furthermore, our library contains >250 non-nuclear proteins containing peptides with activation domain function across both eukaryotic lineages, suggesting that there are unknown biological roles of these peptides beyond transcription. Finally, we identify and validate short, 'universal' eukaryotic activation domains that activate transcription in both yeast and plants with comparable or stronger performance to state-of-the-art activation domains. Overall, our dual host screen provides a blueprint on how to systematically discover novel genetic parts for synthetic biology that function across a wide diversity of eukaryotes.

Significance Statement

Activation domains promote transcription and play a critical role in regulating gene expression. Although the mapping of activation domains from transcription factors has been carried out in previous genome-wide screens, their occurrence in non-transcription factors has been less explored. We utilize an activation domain predictor to mine the entire proteomes of Arabidopsis thaliana and Saccharomyces cerevisiae for new activation domains on non-transcription factor proteins. We validate peptides derived from >750 non-transcription factor proteins capable of activating transcription, discovering many potentially new coactivators in plants. Importantly, we identify novel genetic parts that can function across both species, representing unique synthetic biology tools.

Evolution of molecular switches for regulation of transgene expression by clinically licensed gluconate

Synthetic biology holds great promise to improve the safety and efficacy of future gene and engineered cell therapies by providing new means of endogenous or exogenous control of the embedded therapeutic programs. Here, we focused on gluconate as a clinically licensed small-molecule inducer and engineered gluconate-sensitive molecular switches to regulate transgene expression in human cell cultures and in mice. Several switch designs were assembled based on the gluconate-responsive transcriptional repressor GntR from Escherichia coli. Initially we assembled OFF- and ON-type switches by rewiring the native gluconate-dependent binding of GntR to target DNA sequences in mammalian cells. Then, we utilized the ability of GntR to dimerize in the presence of gluconate to activate gene expression from a split transcriptional activator. By means of random mutagenesis of GntR combined with phenotypic screening, we identified variants that significantly enhanced the functionality of the genetic devices, enabling the construction of robust two-input logic gates. We also demonstrated the potential utility of the synthetic switch in two in vivo settings, one employing implantation of alginate-encapsulated engineered cells and the other involving modification of host cells by DNA delivery. Then, as proof-of-concept, the gluconate-actuated genetic switch was connected to insulin secretion, and the components encoding gluconate-induced insulin production were introduced into type-1 diabetic mice as naked DNA via hydrodynamic tail vein injection. Normoglycemia was restored, thereby showcasing the suitability of oral gluconate to regulate in situ production of a therapeutic protein.

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.

Pooled endogenous protein tagging and recruitment for scalable discovery of effectors for induced proximity therapeutics

The field of induced proximity therapeutics is in its ascendancy but is limited by a lack of scalable tools to systematically explore effector-target protein pairs in an unbiased manner. Here, we combined Scalable POoled Targeting with a LIgandable Tag at Endogenous Sites (SPOTLITES) for the high-throughput tagging of endogenous proteins, with generic small molecule-based protein recruitment to screen for novel proximity-based effectors. We apply this methodology in two orthogonal screens for targeted protein degradation: the first using fluorescence to monitor target protein levels directly, and the second using a cellular growth phenotype that depends on the degradation of an essential protein. Our screens revealed a multitude of potential new effector proteins for degradation and converged on members of the CTLH complex which we demonstrate potently induce degradation. Altogether, we introduce a platform for pooled induction of endogenous protein-protein interactions that can be used to expand our toolset of effector proteins for targeted protein degradation and other forms of induced proximity.

High-throughput discovery and characterization of viral transcriptional effectors in human cells

Viruses encode transcriptional regulatory proteins critical for controlling viral and host gene expression. Given their multifunctional nature and high sequence divergence, it is unclear which viral proteins can affect transcription and which specific sequences contribute to this function. Using a high-throughput assay, we measured the transcriptional regulatory potential of over 60,000 protein tiles across ∼1,500 proteins from 11 coronaviruses and all nine human herpesviruses. We discovered hundreds of transcriptional effector domains, including a conserved repression domain in all coronavirus Spike homologs, dual activation-repression domains in viral interferon regulatory factors (VIRFs), and an activation domain in six herpesvirus homologs of the single-stranded DNA-binding protein that we show is important for viral replication and late gene expression in Kaposi's sarcoma-associated herpesvirus (KSHV). For the effector domains we identified, we investigated their mechanisms via high-throughput sequence and chemical perturbations, pinpointing sequence motifs essential for function. This work massively expands viral protein annotations, serving as a springboard for studying their biological and health implications and providing new candidates for compact gene regulation tools.

The trans-regulatory landscape of gene networks in plants

The transcriptional effector domains of transcription factors play a key role in controlling gene expression; however, their functional nature is poorly understood, hampering our ability to explore this fundamental dimension of gene regulatory networks. To map the trans-regulatory landscape in a complex eukaryote, we systematically characterized the putative transcriptional effector domains of over 400 Arabidopsis thaliana transcription factors for their capacity to modulate transcription. We demonstrate that transcriptional effector activity can be integrated into gene regulatory networks capable of elucidating the functional dynamics underlying gene expression patterns. We further show how our characterized domains can enhance genome engineering efforts and reveal how plant transcriptional activators share regulatory features conserved across distantly related eukaryotes. Our results provide a framework to systematically characterize the regulatory role of transcription factors at a genome-scale in order to understand the transcriptional wiring of biological systems.

Integrative approaches to study enhancer-promoter communication

The spatiotemporal control of gene expression in complex multicellular organisms relies on noncoding regulatory sequences such as enhancers, which activate transcription of target genes often over large genomic distances. Despite the advances in the identification and characterization of enhancers, the principles and mechanisms by which enhancers select and control their target genes remain largely unknown. Here, we review recent interdisciplinary and quantitative approaches based on emerging techniques that aim to address open questions in the field, notably how regulatory information is encoded in the DNA sequence, how this information is transferred from enhancers to promoters, and how these processes are regulated in time.

The next wave of interactomics: Mapping the SLiM-based interactions of the intrinsically disordered proteome

Short linear motifs (SLiMs) are a unique and ubiquitous class of protein interaction modules that perform key regulatory functions and drive dynamic complex formation. For decades, interactions mediated by SLiMs have accumulated through detailed low-throughput experiments. Recent methodological advances have opened this previously underexplored area of the human interactome to high-throughput protein-protein interaction discovery. In this article, we discuss that SLiM-based interactions represent a significant blind spot in the current interactomics data, introduce the key methods that are illuminating the elusive SLiM-mediated interactome of the human cell on a large scale, and discuss the implications for the field.

Interaction modules that impart specificity to disordered protein

Intrinsically disordered regions (IDRs) are especially enriched among proteins that regulate chromatin and transcription. As a result, mechanisms that influence specificity of IDR-driven interactions have emerged as exciting unresolved issues for understanding gene regulation. We review the molecular elements frequently found within IDRs that confer regulatory specificity. In particular, we summarize the differing roles of disordered low-complexity regions (LCRs) and short linear motifs (SLiMs) towards selective nuclear regulation. Examination of IDR-driven interactions highlights SLiMs as organizers of selectivity, with widespread roles in gene regulation and integration of cellular signals. Analysis of recurrent interactions between SLiMs and folded domains suggests diverse avenues for SLiMs to influence phase-separated condensates and highlights opportunities to manipulate these interactions for control of biological activity.

Large-scale mapping and mutagenesis of human transcriptional effector domains

Human gene expression is regulated by more than 2,000 transcription factors and chromatin regulators1,2. Effector domains within these proteins can activate or repress transcription. However, for many of these regulators we do not know what type of effector domains they contain, their location in the protein, their activation and repression strengths, and the sequences that are necessary for their functions. Here, we systematically measure the effector activity of more than 100,000 protein fragments tiling across most chromatin regulators and transcription factors in human cells (2,047 proteins). By testing the effect they have when recruited at reporter genes, we annotate 374 activation domains and 715 repression domains, roughly 80% of which are new and have not been previously annotated3-5. Rational mutagenesis and deletion scans across all the effector domains reveal aromatic and/or leucine residues interspersed with acidic, proline, serine and/or glutamine residues are necessary for activation domain activity. Furthermore, most repression domain sequences contain sites for small ubiquitin-like modifier (SUMO)ylation, short interaction motifs for recruiting corepressors or are structured binding domains for recruiting other repressive proteins. We discover bifunctional domains that can both activate and repress, some of which dynamically split a cell population into high- and low-expression subpopulations. Our systematic annotation and characterization of effector domains provide a rich resource for understanding the function of human transcription factors and chromatin regulators, engineering compact tools for controlling gene expression and refining predictive models of effector domain function.

Systematic discovery of recombinases for efficient integration of large DNA sequences into the human genome

Large serine recombinases (LSRs) are DNA integrases that facilitate the site-specific integration of mobile genetic elements into bacterial genomes. Only a few LSRs, such as Bxb1 and PhiC31, have been characterized to date, with limited efficiency as tools for DNA integration in human cells. In this study, we developed a computational approach to identify thousands of LSRs and their DNA attachment sites, expanding known LSR diversity by >100-fold and enabling the prediction of their insertion site specificities. We tested their recombination activity in human cells, classifying them as landing pad, genome-targeting or multi-targeting LSRs. Overall, we achieved up to seven-fold higher recombination than Bxb1 and genome integration efficiencies of 40-75% with cargo sizes over 7 kb. We also demonstrate virus-free, direct integration of plasmid or amplicon libraries for improved functional genomics applications. This systematic discovery of recombinases directly from microbial sequencing data provides a resource of over 60 LSRs experimentally characterized in human cells for large-payload genome insertion without exposed DNA double-stranded breaks.

Deciphering the multi-scale, quantitative cis-regulatory code

Uncovering the cis-regulatory code that governs when and how much each gene is transcribed in a given genome and cellular state remains a central goal of biology. Here, we discuss major layers of regulation that influence how transcriptional outputs are encoded by DNA sequence and cellular context. We first discuss how transcription factors bind specific DNA sequences in a dosage-dependent and cooperative manner and then proceed to the cofactors that facilitate transcription factor function and mediate the activity of modular cis-regulatory elements such as enhancers, silencers, and promoters. We then consider the complex and poorly understood interplay of these diverse elements within regulatory landscapes and its relationships with chromatin states and nuclear organization. We propose that a mechanistically informed, quantitative model of transcriptional regulation that integrates these multiple regulatory layers will be the key to ultimately cracking the cis-regulatory code.

Base editor scanning charts the DNMT3A activity landscape

DNA methylation is critical for regulating gene expression, necessitating its accurate placement by enzymes such as the DNA methyltransferase DNMT3A. Dysregulation of this process is known to cause aberrant development and oncogenesis, yet how DNMT3A is regulated holistically by its three domains remains challenging to study. Here, we integrate base editing with a DNA methylation reporter to perform in situ mutational scanning of DNMT3A in cells. We identify mutations throughout the protein that perturb function, including ones at an interdomain interface that block allosteric activation. Unexpectedly, we also find mutations in the PWWP domain, a histone reader, that modulate enzyme activity despite preserving histone recognition and protein stability. These effects arise from altered PWWP domain DNA affinity, which we show is a noncanonical function required for full activity in cells. Our findings highlight mechanisms of interdomain crosstalk and demonstrate a generalizable strategy to probe sequence-activity relationships of nonessential chromatin regulators.

Systematic identification and characterization of repressive domains in Drosophila transcription factors

All multicellular life relies on differential gene expression, determined by regulatory DNA elements and DNA-binding transcription factors that mediate activation and repression via cofactor recruitment. While activators have been extensively characterized, repressors are less well studied: the identities and properties of their repressive domains (RDs) are typically unknown and the specific co-repressors (CoRs) they recruit have not been determined. Here, we develop a high-throughput, next-generation sequencing-based screening method, repressive-domain (RD)-seq, to systematically identify RDs in complex DNA-fragment libraries. Screening more than 200,000 fragments covering the coding sequences of all transcription-related proteins in Drosophila melanogaster, we identify 195 RDs in known repressors and in proteins not previously associated with repression. Many RDs contain recurrent short peptide motifs, which are conserved between fly and human and are required for RD function, as demonstrated by motif mutagenesis. Moreover, we show that RDs that contain one of five distinct repressive motifs interact with and depend on different CoRs, such as Groucho, CtBP, Sin3A, or Smrter. These findings advance our understanding of repressors, their sequences, and the functional impact of sequence-altering mutations and should provide a valuable resource for further studies.

A Novel CRISPR Interference Effector Enabling Functional Gene Characterization with Synthetic Guide RNAs

While CRISPR interference (CRISPRi) systems have been widely implemented in pooled lentiviral screening, there has been limited use with synthetic guide RNAs for the complex phenotypic readouts enabled by experiments in arrayed format. Here we describe a novel deactivated Cas9 fusion protein, dCas9-SALL1-SDS3, which produces greater target gene repression than first or second generation CRISPRi systems when used with chemically modified synthetic single guide RNAs (sgRNAs), while exhibiting high target specificity. We show that dCas9-SALL1-SDS3 interacts with key members of the histone deacetylase and Swi-independent three complexes, which are the endogenous functional effectors of SALL1 and SDS3. Synthetic sgRNAs can also be used with in vitro-transcribed dCas9-SALL1-SDS3 mRNA for short-term delivery into primary cells, including human induced pluripotent stem cells and primary T cells. Finally, we used dCas9-SALL1-SDS3 for functional gene characterization of DNA damage host factors, orthogonally to small interfering RNA, demonstrating the ability of the system to be used in arrayed-format screening.

Structure and functional mapping of the KRAB-KAP1 repressor complex

Transposable elements are a genetic reservoir from which new genes and regulatory elements can emerge. However, expression of transposable elements can be pathogenic and is therefore tightly controlled. KRAB domain-containing zinc finger proteins (KRAB-ZFPs) recruit the co-repressor KRAB-associated protein 1 (KAP1/TRIM28) to regulate many transposable elements, but how KRAB-ZFPs and KAP1 interact remains unclear. Here, we report the crystal structure of the KAP1 tripartite motif (TRIM) in complex with the KRAB domain from a human KRAB-ZFP, ZNF93. Structure-guided mutations in the KAP1-KRAB binding interface abolished repressive activity in an epigenetic transcriptional silencing assay. Deposition of H3K9me3 over thousands of loci is lost genome-wide in cells expressing a KAP1 variant with mutations that abolish KRAB binding. Our work identifies and functionally validates the KRAB-KAP1 molecular interface, which is critical for a central transcriptional control axis in vertebrates. In addition, the structure-based prediction of KAP1 recruitment efficiency will enable optimization of KRABs used in CRISPRi.

Quantification of Genome Editing and Transcriptional Control Capabilities Reveals Hierarchies among Diverse CRISPR/Cas Systems in Human Cells

CRISPR/Cas technologies have revolutionized the ability to redesign genomic information and tailor endogenous gene expression. Nevertheless, the discovery and development of new CRISPR/Cas systems has resulted in a lack of clarity surrounding the relative efficacies among these technologies in human cells. This deficit makes the optimal selection of CRISPR/Cas technologies in human cells unnecessarily challenging, which in turn hampers their adoption, and thus ultimately limits their utility. Here, we designed a series of endogenous testbed systems to methodically quantify and compare the genome editing, CRISPRi, and CRISPRa capabilities among 10 different natural and engineered Cas protein variants spanning Type II and Type V CRISPR/Cas families. We show that although all Cas protein variants are capable of genome editing and transcriptional control in human cells, hierarchies exist, particularly for genome editing and CRISPRa applications, wherein Cas9 ≥ Cas12a > Cas12e/Cas12j. Our findings also highlight the utility of our modular testbed platforms to rapidly and systematically quantify the functionality of practically any natural or engineered genomic-targeting Cas protein in human cells.

Transcription factors perform a 2-step search of the nucleus

Transcription factors regulate gene expression by binding to regulatory DNA and recruiting regulatory protein complexes. The DNA-binding and protein-binding functions of transcription factors are traditionally described as independent functions performed by modular protein domains. Here, I argue that genome binding can be a 2-part process with both DNA-binding and protein-binding steps, enabling transcription factors to perform a 2-step search of the nucleus to find their appropriate binding sites in a eukaryotic genome. I support this hypothesis with new and old results in the literature, discuss how this hypothesis parsimoniously resolves outstanding problems, and present testable predictions.

Functional inference of gene regulation using single-cell multi-omics

Cells require coordinated control over gene expression when responding to environmental stimuli. Here we apply scATAC-seq and single-cell RNA sequencing (scRNA-seq) in resting and stimulated human blood cells. Collectively, we generate ~91,000 single-cell profiles, allowing us to probe the cis-regulatory landscape of the immunological response across cell types, stimuli, and time. Advancing tools to integrate multi-omics data, we develop functional inference of gene regulation (FigR), a framework to computationally pair scA-TAC-seq with scRNA-seq cells, connect distal cis-regulatory elements to genes, and infer gene-regulatory networks (GRNs) to identify candidate transcription factor (TF) regulators. Utilizing these paired multi-omics data, we define domains of regulatory chromatin (DORCs) of immune stimulation and find that cells alter chromatin accessibility and gene expression at timescales of minutes. Construction of the stimulation GRN elucidates TF activity at disease-associated DORCs. Overall, FigR enables elucidation of regulatory interactions across single-cell data, providing new opportunities to understand the function of cells within tissues.

Single-cell methods for measuring chromatin accessibility (ATAC-seq) and gene expression (RNA-seq) are rapidly evolving, but tools to integrate data and infer gene-regulatory relationships remain limited. Here we generate multi-omics data of resting and stimulated human blood cells and present a new computational framework for constructing gene-regulatory networks (GRNs). Specifically, we describe functional inference of gene regulation (FigR), a workflow to (1) pair scATAC-seq with scRNA-seq, (2) connect cis-regulatory elements to target genes, and (3) identify TF-gene relationships.

A Novel Gene Controls a New Structure: PiggyBac Transposable Element-Derived 1, Unique to Mammals, Controls Mammal-Specific Neuronal Paraspeckles

Although new genes can arrive from modes other than duplication, few examples are well characterized. Given high expression in some human brain subregions and a putative link to psychological disorders [e.g., schizophrenia (SCZ)], suggestive of brain functionality, here we characterize piggyBac transposable element-derived 1 (PGBD1). PGBD1 is nonmonotreme mammal-specific and under purifying selection, consistent with functionality. The gene body of human PGBD1 retains much of the original DNA transposon but has additionally captured SCAN and KRAB domains. Despite gene body retention, PGBD1 has lost transposition abilities, thus transposase functionality is absent. PGBD1 no longer recognizes piggyBac transposon-like inverted repeats, nonetheless PGBD1 has DNA binding activity. Genome scale analysis identifies enrichment of binding sites in and around genes involved in neuronal development, with association with both histone activating and repressing marks. We focus on one of the repressed genes, the long noncoding RNA NEAT1, also dysregulated in SCZ, the core structural RNA of paraspeckles. DNA binding assays confirm specific binding of PGBD1 both in the NEAT1 promoter and in the gene body. Depletion of PGBD1 in neuronal progenitor cells (NPCs) results in increased NEAT1/paraspeckles and differentiation. We conclude that PGBD1 has evolved core regulatory functionality for the maintenance of NPCs. As paraspeckles are a mammal-specific structure, the results presented here show a rare example of the evolution of a novel gene coupled to the evolution of a contemporaneous new structure.

Dynamic spreading of chromatin-mediated gene silencing and reactivation between neighboring genes in single cells

In mammalian cells genes that are in close proximity can be transcriptionally coupled: silencing or activating one gene can affect its neighbors. Understanding these dynamics is important for natural processes, such as heterochromatin spreading during development and aging, and when designing synthetic gene regulation circuits. Here, we systematically dissect this process in single cells by recruiting and releasing repressive chromatin regulators at dual-gene synthetic reporters, and measuring how fast gene silencing and reactivation spread as a function of intergenic distance and configuration of insulator elements. We find that silencing by KRAB, associated with histone methylation, spreads between two genes within hours, with a time delay that increases with distance. This fast KRAB-mediated spreading is not blocked by the classical cHS4 insulators. Silencing by histone deacetylase HDAC4 of the upstream gene can also facilitate background silencing of the downstream gene by PRC2, but with a days-long delay that does not change with distance. This slower silencing can sometimes be stopped by insulators. Gene reactivation of neighboring genes is also coupled, with strong promoters and insulators determining the order of reactivation. Our data can be described by a model of multi-gene regulation that builds upon previous knowledge of heterochromatin spreading, where both gene silencing and gene reactivation can act at a distance, allowing for coordinated dynamics via chromatin regulator recruitment.

Beyond the marks: reader-effectors as drivers of epigenetics and chromatin engineering

Chromatin is a system of proteins and DNA that regulates chromosome organization and gene expression in eukaryotes. Essential features that support these processes include biochemical marks on histones and DNA, 'writer' enzymes that generate or remove these marks and proteins that translate the marks into transcriptional regulation: reader-effectors. Here, we review recent studies that reveal how reader-effectors drive chromatin-mediated processes. Advances in proteomics and epigenomics have accelerated the discovery of chromatin marks and their correlation with gene states, outpacing our understanding of the corresponding reader-effectors. Therefore, we summarize the current state of knowledge and open questions about how reader-effectors impact cellular function and human disease and discuss how synthetic biology can deepen our knowledge of reader-effector activity.

Oncogenic inspiration for programmable activators

The discovery of regulatory domains has been limited to the investigation of transcription factors and homologous protein sequences. In this issue of Cell Genomics, motivated by an oncogenic protein fusion, Tak et al.¹ direct the regulatory potential of a nontraditional effector domain to novel genomic loci with fusions to programmable DNA-binding domains.

Compendium of human transcription factor effector domains

Transcription factors (TFs) regulate gene expression by binding to DNA sequences and modulating transcriptional activity through their effector domains. Despite the central role of effector domains in TF function, there is a current lack of a comprehensive resource and characterization of effector domains. Here, we provide a catalog of 924 effector domains across 594 human TFs. Using this catalog, we characterized the amino acid composition of effector domains, their conservation across species and across the human population, and their roles in human diseases. Furthermore, we provide a classification system for effector domains that constitutes a valuable resource and a blueprint for future experimental studies of TF effector domain function.

Proteome-scale identification of transcriptional activators in human cells

In this issue of Molecular Cell, Alerasool et al. (2022) present a proteome-scale functional screen to systematically uncover human proteins that can activate transcription.

Directed mutational scanning reveals a balance between acidic and hydrophobic residues in strong human activation domains

Acidic activation domains are intrinsically disordered regions of the transcription factors that bind coactivators. The intrinsic disorder and low evolutionary conservation of activation domains have made it difficult to identify the sequence features that control activity. To address this problem, we designed thousands of variants in seven acidic activation domains and measured their activities with a high-throughput assay in human cell culture. We found that strong activation domain activity requires a balance between the number of acidic residues and aromatic and leucine residues. These findings motivated a predictor of acidic activation domains that scans the human proteome for clusters of aromatic and leucine residues embedded in regions of high acidity. This predictor identifies known activation domains and accurately predicts previously unidentified ones. Our results support a flexible acidic exposure model of activation domains in which the acidic residues solubilize hydrophobic motifs so that they can interact with coactivators. A record of this paper’s transparent peer review process is included in the supplemental information.

Transcriptional activation domains are poorly conserved, intrinsically disordered regions of the transcription factors that remain difficult to predict from protein sequences. A high-throughput method reveals how strong activation domains require a balance between acidic and hydrophobic residues. This balance powers an accurate predictor of activation domains on human transcription factors.

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.