Spotting the enemy within: Targeted silencing of foreign DNA in mammalian genomes by the Krüppel-associated box zinc finger protein family

Tandem C2H2-type zinc finger proteins (ZFPs) constitute the largest transcription factor family in animals. Tandem-ZFPs bind DNA in a sequence-specific manner through arrays of multiple zinc finger domains that allow high flexibility and specificity in target recognition. In tetrapods, a large proportion of tandem-ZFPs contain Krüppel-associated-box (KRAB) repression domains, which are able to induce epigenetic silencing through the KAP1 corepressor. The KRAB-ZFP family continuously amplified in tetrapods through segmental gene duplications, often accompanied by deletions, duplications, and mutations of the zinc finger domains. As a result, tetrapod genomes contain unique sets of KRAB-ZFP genes, consisting of ancient and recently evolved family members. Although several hundred human and mouse KRAB-ZFPs have been identified or predicted, the biological functions of most KRAB-ZFP family members have gone unexplored. Furthermore, the evolutionary forces driving the extraordinary KRAB-ZFP expansion and diversification have remained mysterious for decades. In this review, we highlight recent studies that associate KRAB-ZFPs with the repression of parasitic DNA elements in the mammalian germ line and discuss the hypothesis that the KRAB-ZFP family primarily evolved as an adaptive genomic surveillance system against foreign DNA. Finally, we comment on the computational, genetic, and biochemical challenges of studying KRAB-ZFPs and attempt to predict how these challenges may be soon overcome.


Introduction
Tandem C2H2-type zinc finger proteins (ZFPs) make up the single largest transcription factor family in mice and humans with approximately 600 and 700 genes, respectively [1]. The largest of several ZFP subtypes are the Krüppel-associated box (KRAB) domain-containing ZFPs, called KRAB-ZFPs in mice and KZNFs in humans (hereafter all referred to as KRAB-ZFPs), with estimates of approximately 200 and 300 genes in mice and humans, respectively [1,2].
KRAB-ZFPs contain a potent KRAB repression domain and tandem arrays of zinc fingers (ZNFs) that mediate DNA binding. What makes KRAB-ZFPs exceptional among other DNA binding transcription factors is their ability to bind to long stretches of DNA by combinatorial use of up to several dozen ZNFs that serve as modular DNA binding units. These exceptional modular DNA binding properties were co-opted for use in gene-editing applications, forming the basis of the first generation of engineered sequence-specific DNA modifying enzymes called zinc finger nucleases [3,4]. However, the natural target sites of mammalian KRAB-ZFPs are largely unknown. Importantly, some of the characterized KRAB-ZFPs are associated with metabolism, differentiation, apoptosis, and other cellular functions [2,5], but overall very few KRAB-ZFPs have been functionally investigated.
KRAB-ZFPs are also unique among transcription factor families in that a large fraction of their members have DNA binding domains that are rapidly evolving. This rapid evolution may contribute to morphological and behavioral evolution by controlling expression of developmental genes [2,6]. However, recent experimental and computational discoveries have provided compelling evidence that perhaps a large proportion of evolutionarily young KRAB-ZFPs function as part of a surveillance system that protects mammalian genomes from infectious retroviruses, their endogenous counterparts, and nonretroviral retrotransposons [7][8][9][10][11]. In this hypothesis, the species-specific amplification and diversification of mammalian KRAB-ZFPs are consequences of ancient and ongoing germ-line colonization events by mobile DNA elements. Here we will summarize recent progress that implicated KRAB-ZFPs as molecular guardians of genomic integrity and discuss the possible connections between anti-mobile DNA KRAB-ZFPs and those that have evolved to fulfill functions beyond genomic defense.

Review
Evolution and function of the KRAB-ZFP family Origins of the KRAB-ZFP family KRAB-ZFPs are believed to have evolved from the Meisetz (PRDM9) gene [12], which has KRAB and SET domains and a tandem array of C2H2 ZNFs. Meisetz homologues have been identified in sea urchins and tunicates, indicating that the ancestral KRAB domain arose before the common deuterostome ancestor of vertebrates and echinoderms at least 520 million years ago [12]. However, unlike most mammalian KRAB-ZFPs, which interact with KAP1 (also known as TRIM28 or TIF1β) and are therefore potential transcriptional repressors, Meisetz acts as an H3K4 methyltransferase through its SET domain [13,14]. Thus the ancestral KRAB domain might have acted as a transcriptional activator instead of a repressor [12,15] and changes in the KRAB domain or the evolution of new co-repressors may have facilitated KRAB-ZFPs to function as repression factors. Indeed, KRAB domains evolved rapidly [16], and the TRIM family, to which KAP1 belongs, is highly diversified in vertebrates [17]. Although it is not known when KRAB-ZFPs began recruiting KAP1 or possibly other KAP1-like co-repressors, it seems that a Meisetz-derived KRAB-ZFP lost its SET domain at some time during evolution. Functional changes in the KRAB domains and/or KAP1 may then have resulted in a novel, highly specific transcriptional repression factor that rapidly amplified and diversified throughout tetrapod evolution.
KRAB-ZFPs with a SCAN domain have been found in mammals and lizards but are absent in frog and chicken, indicating that this domain was acquired in KRAB-ZFPs around the root of the amniote branch but subsequently got lost in some species [18]. Interestingly, the SCAN domain shows striking homology to the C-terminal portion of the gag capsid protein from the Gmr1-like family of Gypsy/Ty3-like LTR retrotransposons. It was therefore hypothesized that a retrotransposon insertion into a KRAB-ZFP gene resulted in the exaptation of this domain [18]. However, the biological function of SCAN domains in mammalian KRAB-ZFPs is entirely unknown.

DNA binding and initiation of epigenetic silencing
To date, the vast majority of KRAB-ZFP research has focused on human and mouse KRAB-ZFPs. Therefore, it is important to point out that the findings from these studies may not always be applicable to tetrapod KRAB-ZFPs in general. Nevertheless, all KRAB-ZFPs contain tandem arrays of up to 36 C2H2-type ZNFs, usually encoded by a single exon at the 3' end of the gene [2,6,8]. Each ZNF directly interacts with three consecutive nucleotides and one nucleotide of the reverse-complement strand within the adjacent trinucleotide (Fig. 1). The amino acids mainly responsible for the DNA interaction, and therefore the binding specificity of ZNFs, are located at positions −1, 2, 3, and 6 of the DNA-contacting alpha helix. The looped structure of the ZNF is stabilized by a zinc ion that is characteristically contacted by two cysteine and histidine residues [19] (Fig. 1). Unlike most transcription factors that bind rather short DNA sequences, KRAB-ZFPs can use their tandem ZNF array structure to specifically target large stretches of DNA that are unlikely to be found in significant numbers in the genome by chance. On average, mouse and human KRAB-ZFPs have about eight ZNFs, thus the average KRAB-ZFP target motif is expected to have about 24 nucleotides [1]. However, it has been argued that not all ZNFs of a single KRAB-ZFP are necessarily involved Fig. 1 Model of KRAB-ZFP binding to DNA and induction of heterochromatin formation. Protein-DNA interaction between ZNFs and DNA are mainly mediated by four amino acids at positions −1, 2, 3, and 6 of the α-helix (colored circles). KAP1 is recruited through the KRAB domain and interacts with the NURD/HDAC repressor complex and histone methyltransferases (HMTs) (e.g. SETDB1), which catalyze the removal of H3K9ac and the addition of H3K9me3, respectively. HP1γ interacts with both KAP1 and H3K9me3. DNA methyltransferases (DNMTs) methylate genomic CpG sites, leading to inheritable silencing in DNA binding [20][21][22]. For example, Gli, a non-KRAB containing ZNF protein, uses only four of its five ZNFs to interact with DNA [23], and ZNF91, one of the largest known KRAB-ZFP with 36 ZNFs, mainly uses the 11 most N-terminal ZNFs to bind efficiently to its genomic targets, with some 12 C-terminal ZNFs being dispensable for DNA-recognition [8]. Furthermore, CTCF, a KRAB-less tandem-ZFP with 11 ZNFs was shown to bind to various motifs via clustering its ZNFs in several combinations [24].
Several DNA binding prediction models for tandem-ZFPs have been developed by using bacterial one-hybrid systems, empirical calculations of pairwise amino acidnucleotide interaction energies, and knowledge from the X-ray crystal structure of a three-fingered, C2H2-type ZFP (Zlf268 or Egr1) [4,9,[25][26][27]. However, even empirically based tandem-ZFP binding predictions generally rely on data gained by testing the DNA binding preferences of individual ZNFs in heterologous hybrid proteins and in a nonchromosomal context, which may not always reflect their true DNA binding specificity. Furthermore, amino acids of ZNFs other than the four "specificity residues" can influence binding specificity [9] and ZNF "context" may contribute to a given finger's preferred binding site. In addition, some ZNFs within an array may not interact with DNA and the DNA sequence flanking the binding motif can interfere with ZNF binding [24]. Moreover, certain ZNFs bind specifically to methylated DNA [28,29], indicating that epigenetic modifications can also influence DNA binding of ZFPs. Therefore, it will be very difficult if not impossible to ever reliably predict genomic tandem-ZFP binding sites without experimental testing. Nevertheless, solving the crystal structure of several large KRAB-ZFPs bound to DNA would be a technical milestone that may be necessary to help improve existing models.

Expansion and diversification of KRAB-ZFP genes
Tandem-ZFP genes are predominantly organized in genomic clusters [2,6]. For instance, about one-third of all human tandem-ZFP genes are located within six clusters on chromosome 19, the largest one containing 72 tandem-ZFP genes within a 3.5-Mb region [1]. Generally, human tandem-ZFPs that are located in the same cluster also group together phylogenetically, indicating that they result from local gene duplication events [6]. Through chromosomal translocations and other genomic rearrangements new tandem-ZFP genes can gradually disperse and act as seeds for new clusters [1,6]. Interestingly, the chromatin landscape of KRAB-ZFP clusters is distinct from the rest of the genome. A recent method (in situ Hi-C) generated a 3D map of the human genome and correlated this information with epigenetic marks. This analysis revealed that KRAB-ZFP clusters possess a unique chromatin organization, consisting of both active (e.g. H3K36me3) and repressive H3K9me3 histone modfications [44]. These findings were consistent with two previous studies: the first finding enriched levels of HP1 at the 3' ends of KRAB-ZFP genes [45], and the second demonstrating a combination of low CpG density in gene bodies together with H3K9me3 and H3K36me3 at KRAB-ZFP loci [46]. Interestingly, the human KRAB-ZFP ZNF274 is itself responsible for H3K9me3 enrichment at the ZNF regions of KRAB-ZFP genes [47]. It has been speculated that the recruitment of KAP1 and H3K9me3 to the 3' end of KRAB-ZFP genes (and perhaps more broadly the unique combination of histone marks over KRAB-ZFP clusters) protects against ectopic, non-allelic homologous recombination to some degree [47,48]. Nevertheless, tandem-ZFP gene duplications occurred frequently over evolutionary time scales and resulted in several hundreds of KRAB-ZFP genes in mammals. This rapid amplification is likely catalyzed by the repetitive ZNFs of these genes, which are prone to illegitimate recombination and replication slippage [1,49]. Thus, a fine evolutionary balance has been reached between the need for evolvability of KRAB-ZFPs, and the need to prevent loss of important KRAB-ZFPs by recombination events.
Although the KRAB domains and amino acids required for ZNF structure are generally well conserved amongst mammalian KRAB-ZFPs, positive selection at the residues that confer DNA binding specificity is common, especially between recently duplicated gene pairs [1,[49][50][51]. Thus, once a KRAB-ZFP gene duplicates, one of the pair may keep its ZNF structure to fulfill its original function, whereas the daughter gene becomes available to alter its DNA binding specificity and potentially gain new functions.
Why are there so many KRAB-ZFP genes in mammals?
Although the progenitor of the KRAB domain apparently dates to the last common deuterostome ancestor of chordates and echinoderms [12], KRAB-ZFPs are with few exceptions restricted to tetrapod vertebrates and are most abundant in mammals [1,52]. Interestingly, the KRAB domain of KRAB-ZFPs has changed significantly during tetrapod evolution. A computational analysis of vertebrate KRAB domains revealed that in chicken, lizard and frog KRAB-ZFPs, some of the amino acids that are essential for KRAB-KAP1 interaction in mammals are not conserved [15,31]. It is therefore unclear whether KRAB-ZFPs in these species can recruit KAP1. Furthermore, some of the oldest mammalian KRAB-ZFPs do not interact with KAP1 [15,31], and some function as transcriptional activators instead [53,54]. Thus, structural changes to the KRAB domain at some point during tetrapod evolution may have caused KRAB-ZFPs to recruit KAP1, establishing a new class of epigenetic repressors that subsequently rapidly amplified. However, experimental testing of interactions between nonmammalian KRAB domains and KAP1 and possibly other factors will be required before such a conclusion can be drawn with certainty. Interestingly, tandem-ZFPs with other domains than KRAB have expanded via gene duplications in insects and amphibians [55,56]. This indicates that lineage-specific tandem-ZFP amplification and diversification is not restricted to KRAB-ZFPs.
Some KRAB-ZFPs have been associated with metabolism, differentiation, apoptosis, and human diseases [2,5,57,58], but in most cases their genomic binding sites are unknown. Since the majority of KRAB-ZFPs are predicted to interact with KAP1, most KRAB-ZFPs are believed to repress transcription. Indeed, one of the first genome-wide studies of KRAB-ZFP DNA binding identified binding sites for ZNF263, a human KRAB-ZFP with a SCAN domain, near gene promoters. Importantly, ZNF263 knockdown derepressed a subset of ZNF263targeted genes [59].
The rapid amplification and diversification of KRAB-ZFPs in tetrapods and especially mammals suggest that the bulk of recently emerged KRAB-ZFPs are involved in functions specific to these animals. Although a recent analysis of transcription factor expression during human fetal development demonstrated that the KRAB-ZFP family generally displays less tissue-specific expression levels than other transcription factor families (Siebenthall, K.T., personal communication), a fraction of KRAB-ZFPs are differentially expressed in adult tissues [60]. Notably, many KRAB-ZFPs are highly expressed in evolutionarily recent tissues, such as the mammalian-specific placenta [15]. Mammals might also require a large number of KRAB-ZFPs to control mammal-specific innovations in processes such as erythropoiesis [61] and development of the adaptive immune system [62]. Since even closely related species such as higher primates differ in their KRAB-ZFP arsenal, it was suggested that KRAB-ZFPs may also contribute to human brain development [63]. Interestingly, the transcriptional activity of certain KRAB-ZFP orthologues greatly varies between human and chimpanzee brain, suggesting that KRAB-ZFPs may change their expression levels after speciation to adapt to new functions [51,63]. While the increase in the complexity of mammalian development might explain some of the KRAB-ZFP diversification, the number of KRAB-ZFPs does not correlate with brain size or the duration of embryonic development [64]. According to a recent study, opossums have nearly twice as many KRAB-ZFP genes as humans [15]. Indeed, recent findings support the hypothesis that the majority of KRAB-ZFPs function as repressors of parasitic DNA rather than as conventional gene-regulating transcription factors. Moreover, many KRAB-ZFPs might regulate genes through targeting nearby remnants of parasitic DNA that has been coopted as novel regulatory sequences. In the following section, we will briefly discuss the impact of parasitic DNA elements on mammalian evolution and review recent findings suggesting that those elements triggered KRAB-ZFP expansion and diversification.

KRAB-ZFPs are adaptive repressors of foreign DNA Retrotransposons and the need for an adaptive repression system
Retroviruses have been invading mammalian germ lines for millions of years, accumulating in the form of endogenous retroviruses (ERVs) that account for approximately 8 % of the human genome [65]. Mammalian genomes also contain many nonretroviral retrotransposons-long interspersed nuclear elements (LINE) and short interspersed nuclear elements (SINEs)-that cannot form infectious particles but amplify through retrotransposition in host cells.
Although no replication-competent human ERVs (HERVs) have been described yet, HERV fragments are associated with human lymphomas and other cancer types, and HERV-K particles were detected in human preimplantation embryos [70,[77][78][79]. The HERV-K subgroup HML2 is responsive to the HIV-1 transactivator protein (Tat) [80], and some of these elements encode functional envelope and integrase proteins. HERV-K envelope proteins can be incorporated into HIV particles [81] and may be a biomarker for HIV latency [82]. Furthermore, the long terminal repeats (LTRs) of HERV-K contain many binding sites for inflammatory transcription factors, suggesting that these ERVs contribute to the pathology of inflammatory disease [83]. Finally, several non-retroviral retrotransposons have been linked to many human diseases [84][85][86][87], and it has been speculated that LINEs decrease longevity by eroding genomic integrity [88]. These studies highlight the potentially damaging effects of uncontrolled activation of retrotransposons.
To defend their genomes against exogenous retroviruses and EREs, mammals rely on a wide range of defense mechanisms, including APOBEC proteins [89], PIWIinteracting RNAs (piRNAs) [90], nucleic acid sensors [91], and transcriptional repression [92][93][94]. EREs are transcriptionally repressed by stable epigenetic silencing that can be maintained through cell division. This silencing mechanism requires distinct and partially overlapping machinery in pluripotent and somatic tissues. In somatic tissues, EREs are repressed by DNA methylation, as revealed by mutations in DNA methyltransferases [95][96][97]. In pluripotent embryonic stem cells (ESCs), ERVs and LINEs are repressed primarily by machinery that creates repressive histone modifications, most notably H3K9me3 [98], whereas DNA methylation is largely dispensable [99,100]. However, histone modifications established in early development seem to be required to initiate or stabilize heritable DNA methylation at EREs during differentiation and development [36,101].
Although epigenetic repression of retroviral DNA in ESCs has been the subject of numerous studies, the factors that target epigenetic silencing machinery to EREs in mammals have remained elusive. One strategy developed by eukaryotes to cope with mobile DNA diversity relies on short RNAs that are expressed by the transposons themselves and allow guiding of the silencing machinery to the expressed element through base pairing [102]. These short RNAs include small interfering RNAs (siRNAs) that regulate LTR transposons in yeast [103], siRNAs that target DNA methylation at heterochromatin in plants [104], and piRNAs that guide silencing in animal germ cells [90,[105][106][107][108]. Tetrapods likely employ an additional and equally important strategy to recognize and silence EREs: genetic encoding of an army of evolutionarily selected DNA-binding transcription factors. To silence EREs with minimal off-target effects, such transcription factors need to be able to bind large DNA motifs that are unlikely to appear in the genome by chance. Furthermore, these factors need to be evolutionarily adaptable to recognize newly emerging EREs and possess a potent repression domain to stably silence these elements. The KRAB-ZFP family alone fulfills all these criteria. Indeed, while small RNAs may be the predominant way to target EREs in plants and mammalian germ cells, mammals seem to rely on KRAB-ZFPs to recognize and silence retroviruses and EREs during early embryonic development.

KRAB-ZFPs repress exogenous and endogenous retroviruses
The most compelling direct evidence that KRAB-ZFPs repress retroviruses and EREs comes from the identification of the ZFP809/KAP1 murine leukemia virus (MuLV) repression complex and from two loss-of-function studies of the KRAB-ZFP corepressors KAP1 and SETDB1 in ESCs that revealed ERV activation phenotypes.
It has been long known that a multi-component repressor complex binds to a 17-bp sequence within the proline tRNA primer binding site (PBS pro ) of integrated MuLV in murine pluripotent stem cells [109,110]. The identification of KAP1 as an integral component of the PBS pro targeting repressor complex [111] strongly implied that a KRAB-ZFP is the DNA binding factor that tethers the complex to MuLV. Indeed, ZFP809, a mouse KRAB-ZFP with no human orthologue, was subsequently identified as the recognition module that targets the PBS pro and recruits KAP1 [112]. Shortly thereafter, genetic removal of KAP1 or its interacting protein SETDB1 revealed a broad requirement for these proteins in heterochromatin formation and ERV silencing in ESCs and in viability [36,100,113]. However, knockout of KAP1 in murine embryonic fibroblasts (MEFs) does not affect ERV expression [113], consistent with the observation that KAP1 repression during early embryogenesis leads to irreversible silencing that is maintained by DNA methylation and does not persistently require KAP1 [35]. These studies thus laid the foundation supporting the KRAB-ZFP family as the likely candidate for ERV recognition and transcriptional silencing.
Direct evidence linking an individual KRAB-ZFP to ERE silencing came from genome-wide binding and genetic knockout studies of ZFP809. In ESCs, ZFP809 binds to several PBS pro containing ERVs and recruits the KAP1/ SETDB1 repressor complex to these elements. Moreover, ZFP809 knockout leads to a strong reactivation of VL30 elements with a PBS pro (VL30 Pro ) in postimplantation embryos and in most organs and tissues of adult animals [7]. Interestingly, VL30 Pro elements are inactive in preimplantation embryos and ESCs even in the absence of ZFP809, most likely because certain transcription factors are missing in these embryos/cells. Nevertheless, ZFP809 is required to initiate epigenetic silencing of these elements in ESCs to prevent ERV reactivation during differentiation. Once silenced by ZFP809 in ESCs, VL30 Pro remain transcriptionally silent in differentiated cells, even when ZFP809 is no longer present [7]. These findings support the model in which KRAB/KAP1 silencing is initiated in early embryos and heritably maintained in somatic tissues without a continuous requirement for KRAB-ZFPs or KAP1 [35].
However, KAP1 is also required for ERV silencing in neural progenitor cells, indicating that ERV repression by KRAB/KAP1 is not strictly restricted to ESCs [114].
Moreover, SETDB1 appears to be continuously required to maintain ERV silencing in some differentiated cell types, as conditional SETDB1 deletion in MEFs and B lymphocytes leads to massive reactivation of several ERV groups [7,115]. Importantly, different groups of ERVs become reactivated in SETDB1 knockout ESCs and B lymphocytes, indicating that not only the loss of repressive chromatin marks but also the presence of possibly tissue-specific transcription factors determines which ERVs become de-repressed [115].
Although ZFP809 is so far the only KRAB-ZFP whose role in ERV silencing is supported by convincing biochemical and genetic evidence, several other KRAB-ZFPs have been implicated (Table 1). For example, ZFP819 knockdown led to a significant upregulation of IAP ERVs and other EREs in murine ESCs. Although a defined target motif for ZFP819 has not been identified in these elements, overexpression of ZFP819 also inhibits expression of a luciferase reporter containing an IAP LTR fragment [116].
Several lines of evidence suggest that human ERVs (HERVs) are repressed by KRAB-ZFPs. KAP1 is enriched at Class I and II HERVs in human ESCs and although the KRAB-ZFPs that recruit KAP1 to these elements have not been identified, a 39-bp sequence was demonstrated to be critical for KAP1-dependent silencing of HERV-K elements in reporter assays [11]. Intriguingly, this sequence overlapped with the PBS of these ERVs, which is complementary to a human lysine tRNA. Moreover, many human KRAB-ZFPs interact with specific ERV classes when overexpressed as GFP-fusion proteins in 293 T cells [9]. However, none of these interactions were validated by other types of binding or functional assays. Thus it is premature to conclude that all these proteins are indeed ERV silencers.

KAP1/KRAB-ZFP repression of nonretroviral retrotransposons
The first indirect evidence that KRAB-ZFPs also repress nonretroviral EREs came from two studies that investigated the genome-wide binding patterns of KAP1 in human ESCs and primary human T lymphocytes [10,11]. These studies showed that KAP1 is bound to a defined subset of LINE-1 (L1) transposons and several groups of SINE-VNTR-Alu (SVA) elements. The identified target sequences in these elements induced epigenetic silencing of reporter genes in human ESCs. Furthermore, the murine KRAB-ZFP Gm6871 was identified as a L1 binding protein by ChIP-seq [10].
Shortly thereafter, it was shown that several humanand primate-specific SVA and L1 elements are derepressed in trans-chromosomic murine ESCs that contain a copy of human chromosome 11 [8]. This indicated that the repression factors controlling these elements in humans are absent in mice. In a subsequent screen of a selection of 14 highly expressed primate-specific KRAB-ZFPs that appeared in the catarrhine lineage 25-35 million years ago, ZNF91 and ZNF93 repressed reporter constructs containing SVA and L1PA-type retrotransposons, respectively. Overexpression of these KRAB-ZFPs in trans-chromosomic murine ESCs resulted in resilencing of their target elements. Furthermore, ChIP-seq confirmed ZNF93 binding to endogenous L1PA elements, KO genetic knockout studies, Rep reporter assays, KD knockdown by siRNA, EMSA electrophoretic mobility shift assays, n.d. not determined providing conclusive evidence that human KRAB-ZFPs bind and repress retrotransposons [8].
Notably, no KRAB-ZFPs that target DNA transposons have been identified to date. This is not surprising since DNA transposons only constitute a small fraction of mammalian genomes. Furthermore, DNA transposons are generally not active in mammals and repression of these elements by KRAB-ZFPs might therefore be unnecessary.
A potential role of KRAB-ZFPs in repressing HIV, HTLV-1, and nonviral transgenes Several years before KAP1 and KRAB-ZFPs were identified as ERE repressors, a handful of studies explored the use of artificially designed KRAB-ZFPs to restrict HIV infection. These artificial KRAB-ZFPs successfully repressed HIV transcription by binding to the proviral LTR or PBS [117][118][119], showing that KRAB-ZFPs can be "designed" to repress a virus of interest. However, this strategy has not been adapted for clinical applications.
Retrovirus and ERE repressing KRAB-ZFPs are generally believed to evolve when a species is continuously exposed to these elements [8][9][10]64]. The identification of these four human KRAB-ZFPs as repressors of lentiviruses is therefore somewhat surprising because ZNF10, ZNF175, ZNF350 and ZNF282 (unlike ZFP809, ZNF91, and ZNF93) are well conserved in mammals. In fact, ZNF282 is one of only three human KRAB-ZFPs that have orthologues in nonmammalian amniotes [15]. Although lentiviruses have been infecting primates for millions of years [124], humans have not been exposed to HIV for more than a few decades and the HIV-related simian immunodeficiency virus was estimated to be only 32,000 years old [125]. Similarly, although HTLV is believed to have infected humans for tens of thousands of years [126] and HTLV-related simian T-cell leukemia virus (STLV) strains have been found in nonhuman primates [127], HTLV clearly emerged after ZNF282. However, it cannot be excluded that these conserved KRAB-ZFPs originally evolved to repress ancient lentiviruses or lentivirus-like elements and therefore still recognize current HIV and HTLV strains.
Although there is some evidence that ZNF175 expression is correlated with HIV infection [128,129], it is unclear whether any of the KRAB-ZFPs mentioned above inhibit HIV or HTLV in vivo. The HIV and HTLV LTRs may simply contain sequences that resemble endogenous targets of these KRAB-ZFPs and binding to these viral sequences is not evolutionarily intended or of any biological consequence. Indeed, the HIV-1 LTR sequence that is necessary for ZNF10-mediated repression overlaps with NF-κB and Sp1 binding sites [123] which are commonly found in gene promoters. Furthermore, ZNF350 was identified as a tumor-suppressor gene [130][131][132][133], ZNF282 was associated with cancer progression [134,135], and ZNF175 might play a role in neuronal survival [136].
Thus, these ancient KRAB-ZFPs may have functions other than retroviral restriction. The question remains why HIV and HTLV did not mutate to evade these potential repressor proteins. Possibly, transcriptional repression by KRAB-ZFPs is advantageous for these viruses under some circumstances (e.g., establishment of latency). Therefore, HIV and other viruses might have actually evolved to be bound by certain KRAB-ZFPs.
Surprisingly, a murine KRAB-ZFP has been associated with silencing of a bacterial transgene in mice. A 0.9-kb fragment of the bacterial xanthine-guanine phosphoribosyltransferase (gpt) gene is rapidly silenced by DNA methylation and histone modifications in mice of some strains (e.g., C57BL/6) but not others (e.g., DBA/2) [137,138]. Breeding mice carrying the unmethylated transgene to mice that methylate it leads to transgene methylation, indicating that a dominant factor, present only in some mouse strains, is required for silencing. This factor, initially named strain-specific modifier 1, was later identified as the KRAB-ZFP gene 2610305D13Rik [139]. Indeed, the observation that silencing commences during implantation and that DNA methylation spreads into neighboring regions is consistent with KRAB-ZFPmediated silencing [137,140]. However, direct binding of this KRAB-ZFP to the transgene has not been shown. Furthermore, It is unclear what the genomic targets of 2610305D13Rik are or whether the gptcontaining transgenic sequence resembles an ERE or another genomic target [139].

An evolutionary arms race between KRAB-ZFPs and EREs
Nonretroviral transposable elements such as SINEs and LINEs are vertically transmitted from generation to generation, allowing the host to evolve repressive mechanisms to restrain their activity. Although these elements can replicate within the genome, their mutation rates are rather limited. In contrast, ERVs are derived from exogenous retroviruses that can be horizontally transmitted between animals. Exogenous retroviruses can evolve much faster than EREs [141] and, in the case of a germ-line colonization, may represent new genomic elements that have few or no similarities to EREs already residing in the host genome. Thus, the host is forced to quickly develop an effective repression mechanism.
Germ-line colonization by ERV-like LTR retrotransposons has been described in insects [142]. However, vertebrates and tetrapods have accumulated a much more diverse retroviral load during evolution than other animals [64,66]. What caused the immense ERV diversity in these animals? Many retroviruses that infect mammals replicate by co-opting properties of immune cells that normally operate during intercellular communication, such as antigen presentation and T-cell activation [143]. Thus, cells of the adaptive immune system, which first emerged in jawed vertebrates [144], might have served as replication sites for retroviruses and therefore favored ERV diversity. On the other hand, adaptive immunity might have exerted selective pressure on retroviruses, contributing to the rapid diversification of these parasites. Either way, it is possible that the KRAB-ZFP repression system evolved in response to the increasingly diverse burden of horizontally transmittable retroviruses in tetrapods.
To our knowledge, only human and murine KRAB-ZFPs have been analyzed in genome-wide binding studies so far. One can therefore only speculate that EREs are indeed the main drivers of KRAB-ZFP diversification in other mammals and non-mammalian tetrapods. Consistent with this idea, a computational analysis of 16 mammalian genomes revealed a striking correlation between the number of endogenous LTR elements and the number of tandem-ZFP genes (Fig. 2) [64]. Moreover, the emergence of new LTR elements correlated with tandem-ZFP gene duplication events [64]. Surprisingly, the same correlation was observed in a selection of nonmammalian vertebrates, in which only few or none of the tandem-ZFPs contain KRAB domains [64]. It was therefore speculated that these species use tandem-ZFPs with alternative repressor domains to restrict EREs [64]. Indeed, the BTB/POZ domain, found in both vertebrate and invertebrate tandem-ZFPs, interacts with HDAC co-repressor complexes [145].
Some KRAB-ZFPs may even restrict retroviral activity without the help of transcriptional corepressors. The KRAB-ZFP associated SCAN domain is not only believed to be derived from a retrotransposon but is also structurally similar to the HIV C-terminal capsid [18]. Furthermore, the ability of the SCAN domain to multimerize by a domain-swapping mechanism resembles the multimerization of capsid domains to form the capsid structure of retroviruses [146,147]. It was therefore speculated that this domain may target KRAB-ZFPs to cytoplasmic retroviral capsids, allowing sequestration of newly synthesized retroviral DNA [18]. Thus, it is imaginable that the exaptation of the SCAN domain enabled KRAB-ZFPs to restrict retroviral activity in nonmammalian tetrapods, possibly before the KRAB domain was able to interact with KAP1. The emergence of KAP1-interacting KRAB domains may then have provided an additional repression mechanism. Importantly, KRAB/KAP1 transcriptionally represses both chromosomal and nonintegrated DNA [148] and KAP1 inhibits genomic integration of HIV [149], suggesting that KAP1 can restrict retroviral replication by multiple mechanisms. Intriguingly, several tandem-ZFPs with a SCAN domain in lizards were predicted to bind Gmr1-like EREs, one of them precisely at the PBS [18]. One may therefore speculate that the SCAN domain played an important role in the evolution of ERE repressing KRAB-ZFPs. However, experimental evidence for retrotransposon repression by the SCAN domain or nonmammalian KRAB-ZFPs in general is still lacking.  Fig. 2 An evolutionary arms race between EREs and KRAB-ZFPs. Estimated number of LTR elements [64] and KRAB-ZFPs [15] in vertebrates. The phylogenetic tree is an approximate reprint of a previously published tree [64] The hypothesis that KRAB-ZFPs evolved to defend host genomes from parasitic DNA implies that many of these transcription factors bind to ERVs and other EREs. Using computational motif prediction tools, it was predicted that many human KRAB-ZFPs bind EREs that entered the human genome around the time these KRAB-ZFPs appeared [150]. Recently, about 70 randomly selected human ZFPs were epitope tagged and expressed in a human cell line to identify their genome-wide binding patterns by ChIP-seq. Of 18 KRAB-ZFPs, 16 bound to some extent to specific EREs, versus only about 10 % of non-KRAB ZFPs [9]. Thus, the majority of human KRAB-ZFPs can bind EREs, although it is not known whether they are required for ERE repression. Intriguingly, the estimated ages of most ERE-binding KRAB-ZFPs correlate with the ages of the EREs they bind [9]. However, two KRAB-ZFPs that are well conserved in mammals, ZNF382 and ZNF33A, bind to currently active LINE L1HS-like elements and hominoid-specific SVA elements, respectively [9]. Thus, these EREs are not likely the primary targets but they simply tolerate KRAB-ZFP binding, or alternatively, these KRAB-ZFPs were recently co-opted to repress EREs. In support of the latter possibility, there is a strong signature of recent positive selection at the ZNF33A locus in humans [9].
Evidence for an ongoing evolutionary arms race between host KRAB-ZFPs and transposable elements also came from the few studies that identified individual ERErepressing KRAB-ZFPs. Both Gm6871 and KAP1 bind predominantly to L1 elements that entered the mouse genome 4-7 million years ago [10]. Similarly, human KAP1 binds preferentially to L1 elements estimated to be 8-27 million years old [10]. Moreover, in human ESCs, the youngest human L1 elements that are not bound by KAP1 are expressed at higher levels than their older counterparts. Depletion of the three DNA methyltransferases strongly upregulated these young L1 elements, but older L1 families were relatively unaffected [10]. The PIWI-piRNA pathway is involved in L1 silencing in human pluripotent stem cells [106]. In a proposed model, expression of newly emerging L1 elements is silenced by the PIWI-piRNA system, which is targeted to these elements by L1-derived piRNAs. Over time, KRAB-ZFPs evolve to recognize those transposable elements and take over repression until their target EREs become too degenerated to be recognized. By that time, the accumulated mutations and deletions of the EREs had already led to their inactivation so repression is no longer required [10]. Indeed, vertebrate genomes have many tandem-ZFP pseudogenes [64], suggesting that many KRAB-ZFP genes have become obsolete after their target EREs had been inactivated by genetic drift.
ZNF91 and ZNF93 emerged in the last common ancestor of apes and Old-World monkeys and are members of a KRAB-ZFP cluster that has amplified and diversified throughout the evolution of apes and humans [49]. Both ZNF91 and ZNF93 underwent several drastic structural changes in the last common ancestor of orangutans and humans 12-18 million years ago [8]. Intriguingly, these changes were crucial for the ability of ZNF91/93 to repress retrotransposons in humans. Around the same time as the ZNF91/93 changes, the ZNF93-targeted L1 elements also changed substantially. A 129-bp sequence within L1PA subgroups is lost in evolutionarily younger L1PA elements, indicating a potential repression escape mutation that allowed these elements to be expressed [8].
Altogether these findings imply that evolutionarily young KRAB-ZFPs repress retroviruses and EREs, whereas older and well-conserved KRAB-ZFPs fulfill other functions. To gain insight into the evolutionary history of ERE repressing KRAB-ZFPs, we generated a phylogenetic tree of murine and human KRAB-ZFPs. Indeed, ZFP809 and Gm6871 are closely related to other murine KRAB-ZFPs but do not phylogenetically cluster with any human KRAB-ZFPs. On the other hand, the L1 and SVA repressors ZNF91/93 form a cluster with many other human KRAB-ZFPs but are not closely related to any mouse KRAB-ZFP (Fig. 3). In contrast, ancient and well-conserved KRAB-ZFPs such as PRDM9, ZFP/ ZNF282 and ZNF/ZFP777 are present as one-to-one orthologues in mice and humans and do not group in clusters of species-specific KRAB-ZFPs (Fig. 3). This indicates that some KRAB-ZFPs are constrained in mammalian evolution whereas others, including the ancestors of ERE repressors such as ZFP809 and ZNF91/93, are prone to frequent gene duplications and diversification. This is in agreement with the model of KRAB-ZFP amplification and diversification as a response to invading foreign DNA.
ZFP809 binds to VL30 and MmERV elements (both ERVs of the ERV1 family) that contain a PBS pro , but many of these elements contain a PBS complementary to a glycine tRNA instead (PBS gly ) [7] (Fig. 4). Thus, even closely related elements within the same ERV group can escape repression by a specific KRAB-ZFP. At the same time, ZFP809 binds weakly to several hundred genomic RLTR10 and MERVL elements, ERVs belonging to the ERVK and ERVL families, respectively (Wolf et al., unpublished data). Although target motifs similar to the PBS pro can be found at these ZFP809 binding sites (Fig. 4), no co-occupation with KAP1 and SETDB1 was observed, presumably because ZFP809 binding is not sufficient to assemble the KAP1 repressor complex at these targets ( [7] and Wolf et al. unpublished data). However, the weak binding affinity of ZFP809 to these elements also shows that imperfect binding sites for an ERV-repressing KRAB-ZFP can appear by chance in unrelated ERVs, possibly because of the general GC richness in these elements. ERV-targeting KRAB-ZFPs might have a general potential to bind weakly to other ERVs. If such a KRAB-ZFP is duplicated or becomes available after its original target ERV was inactivated over time, mutations of the ZNFs might allow stronger binding and therefore functional silencing of newly emerged ERVs.

A possible link between KRAB-ZFPs and ERE adaptation Transcriptional regulation of genes through ERE repressing KRAB-ZFPs
EREs and especially ERVs have a profound impact on patterns of mammalian gene expression. Retroviral LTRs contain strong promoter elements to ensure efficient expression of their proviral genome. When integrated near a cellular gene, transcription from these LTRs can drive the expression of that gene [151,152]. Furthermore enhancers within LTRs can influence the expression of distant cellular genes and contribute to the innovation of gene regulatory networks [68,69,[153][154][155][156].
As discussed above, KRAB-ZFPs that repress newly emerged EREs may result from a duplication of an existing ERE repressor, followed by mutations in the DNA binding domain that leads to recognition of the new targets. Alternatively, gene-targeting KRAB-ZFPs  [197] to infer a Maximum-likelihood phylogenetic tree using MEGA version 6 with default parameters [198]. All KRAB sequences are provided as Additional file 1. Exemplary proportions of the tree that contain ERE-silencing KRAB-ZFPs or KRAB-ZFPs conserved between mouse and human (e.g., PRDM9, ZNF282, and ZNF777) are shown in more detail might duplicate and change their binding specificity towards new EREs (Fig. 5). Furthermore, since EREs have been continually co-opted as gene regulatory elements, it is plausible that some KRAB-ZFPs regulate gene expression by binding to EREs (Fig. 5). Indeed, knockout of KAP1 not only de-repressed ERVs but also many genes near those ERVs [157]. Also ZFP809 knockout led to the upregulation of a handful of cellular genes near ZFP809-targeted ERVs [7].
It has been also reported that an IAP LTR that drives gene expression is silenced by a gene located in a KRAB-ZFP cluster [158]. Although this might be another example of an ERV/KRAB-ZFP that got co-opted as a gene regulator, the identity of the repressor gene is yet to be determined. Furthermore, a human-specific SVA element that integrated into the fibroblast growth factor 2 (FGF2) promoter (Greenberg et al. unpublished data) may enhance the expression of this key protein PBS gly Fig. 4 Differential ZFP809 binding to various ERVs. ZFP809 target sequences identified by ChIP-seq [7] are shown with differences from the canonical PBS pro highlighted in red. Dashed arrows indicate weak ZFP809 binding that is not sufficient to form the KAP1/SETDB1 repressor complex  [159], possibly regulated by ZNF91 [8].

ZFP57 protects genomic imprints in retrotransposed genes
One of the best-characterized KRAB-ZFPs, ZFP57, is required to maintain a subset of genomic imprints in mice [160], and mutations in human ZFP57 have been associated with transient neonatal diabetes [161]. ZFP57 binds to a methylated hexanucleotide within imprinted control regions and recruits KAP1 and SETDB1 to establish H3K9me3 [29,162]. By binding to the methylated imprinted control region, ZFP57 also protects the methylated region from the genome-wide demethylation that occurs during mammalian preimplantation development. Although the imprinted control regions bound by ZFP57 in mammals are not associated with particular EREs, five murine imprinted genes arose by retrotransposition (Mcts2, Nap1l15, U2af1-rs1, Inpp5f_v2, and Peg12), and another two are derived from retrotransposons (Rtl1 and Peg10) [163]. Furthermore, DNA methylation is considered to have evolved primarily as a defense against foreign DNA [95]. Thus, genomic imprinting itself and its maintenance by KRAB-ZFPs might originate from retrotransposon repression.

KRAB-ZFPs regulate sexually dimorphic gene expression patterns through ERV targeting
Sexually dimorphic gene expression in liver is a complex phenomenon in mice and humans. Sex-specific expression of growth hormones can induce gene expression, and genes can be repressed in a sex-specific manner. The mouse regulator of sex-limitation (Rsl) locus encodes two KRAB-ZFP genes, Rsl1 and Rsl2, which are regulated directly in the kidney by androgen or indirectly in the liver by growth hormones [164]. One Rsl-repressed gene, Cyp2d9, is a member of the large cytochrome P450 family, which participates in many metabolic processes, such as detoxification of foreign chemicals, hormone synthesis and breakdown, and cholesterol synthesis [165]. Cyp2d9 and other sexually dimorphic cytochrome P450 genes were also upregulated in KAP1 knockout liver [166], confirming a role for KRAB-ZFPs in sexually dimorphic patterns of gene expression. Interestingly, cytochrome P450 genes have been diversified in mammals through gene duplications and positive selection, similar to KRAB-ZFP genes [167]. Because of their abundance and evolutionary dynamic, these genes might be prone to ERV-mediated repression by KRAB-ZFPs. Indeed, one of the few genes that was upregulated in ZFP809 knockout mice was a cytochrome P450 gene (Cyp4f37) that contained an ancient ERV insertion near the promoter region [7].
More importantly, another target of Rsl, which encodes sex-limited protein (Slp) was reported to be controlled by an ancient ERV LTR located 2 kb upstream of this gene [168,169]. Indeed, Rsl1 binds a defined sequence within this LTR, suggesting that Rsl1-mediated control of Slp evolved from retroviral repression [170]. Intriguingly, Rsl1 is located within a cluster of recently duplicated KRAB-ZFP genes found only in the Mus lineage [171].

From meiotic recombination control to ERV-repression
In most mammals, homologous recombination during meiosis tends to occur at specific segments of the genome. Interestingly, the placement and activity of these so-called hotspots varies greatly between closely related Mus species [172], between humans and primates [173,174], and even between human individuals [175]. These hotspots are enriched in H3K4me3, a histone mark usually found at active and poised enhancers. Human recombination hotspots often contain a 13-mer sequence motif [176]. The KRAB-ZFP PRDM9 governs recombination activity in humans and mice [177][178][179] through the H3K4 trimethyltransferase activity of its SET domain [13,14]. Intriguingly, many human PRDM9 alleles have been identified [180], and accelerated evolution of the PRDM9 DNA binding domain has been reported [181][182][183]. PRDM9 was therefore suggested to bind rapidly evolving repetitive DNA elements [176,181]. Indeed, THE1A and THE1B LTR elements, members of the Mammalian apparent LTR-retrotransposons (MaLRs) family, contain a PRDM9 binding motif and are overrepresented in PRDM9-associated hotspots [183,184].
Could this indicate that PRDM9 originally evolved as an ERV repressor? Although PRDM9 predates the emergence of human PRDM9-bound THE1 LTR elements, MaLRs colonized the genomes of eutherian mammals at least 80-100 million years ago [185]. An ancient connection between LTR elements and PRDM9 in mammals can therefore not be excluded. However, the KRAB domain of PRDM9 lacks the amino acid sequences that have been identified as essential for KAP1 interaction [15,31], indicating that PRDM9 is not a part of the KRAB/KAP1 ERV repression system. Nevertheless, ERV integrations might have re-organized recombination hotspots by introducing new PRDM9 binding sites. Moreover, the mechanism of DNA binding by PRDM9-which is believed to be highly specific yet permissive at the same time [176,186]-and its ability to rapidly change DNA specificity may have been the perfect attributes to trigger expansion and evolution of ERV-repressing KRAB-ZFPs.

Challenges and future directions of KRAB-ZFP research
Despite recent progress, KRAB-ZFPs are not only one of the largest but also one of the least understood transcription factor families in mammals. In fact, many functional KRAB-ZFP genes might not even be annotated yet, whereas some predicted KRAB-ZFPs will turn out to be pseudogenes. The highly repetitive nature of KRAB-ZFP genes makes conventional annotations difficult, and estimates of their copy numbers vary as they strongly depend on the inclusion criteria used [1,2,6,15]. Especially in low-quality genomes, the real number might be underestimated, and alternative splice isoforms may additionally contribute to the diversity of expressed KRAB-ZFPs.
Moreover, the assumption that all KRAB-ZFPs are DNA binding transcription factors might be premature and it is possible that some tandem-ZFPs function outside the nucleus. Binding of C2H2-type ZNFs to RNA and proteins has been reported [187,188], suggesting that some KRAB-ZFPs do not act as DNA binding transcription factors. Cross-linking and immunoprecipitation followed by RNA sequencing (CLIP-seq) analysis of"orphan" C2H2type ZNFs not known to bind DNA or for characterized ZNFs that contain orphan ZNFs might yield novel insights into RNA biology and RNA recognition.
A key to understanding the KRAB-ZFPs that function as DNA binding transcription factors is to determine their genome-wide binding patterns. Although ZNF prediction tools are improving, they cannot, and perhaps never will, be used to reliably predict genome-wide DNA binding sites. We used several of these tools to predict a ZFP809 binding motif and compared the outcome with the experimentally determined ZFP809 binding site, the PBS pro [7,112]. Although the predicted motifs showed some similarity to the PBS pro , the fraction of overlapping nucleotides was rather small (Fig. 6a). ZFP809 and probably most KRAB-ZFPs tolerate very few mismatches for efficient binding [112]. Therefore, every single falsely predicted nucleotide drastically increases the proportion of falsely predicted binding sites in the genome. To test the accuracy of the predicted ZFP809 binding motif, we screened the mouse genome for targets resembling this motif and analyzed the 500 top-scored genomic sites for ZFP809 enrichment using published ChIP-seq data [7]. Indeed, ZFP809 was not enriched at these predicted genomic binding sites (Fig. 6b). Although the predicted motifs of some KRAB-ZFPs strikingly resemble the experimentally determined target motif [9], the case of ZFP809 highlights that KRAB-ZFP binding prediction without experimental testing remains highly unreliable.
Improved ChIP-seq protocols such as ChIP-exo and ChIP-nexus allow transcription factor binding sites to be determined at near nucleotide resolution [189,190]. These new techniques and the reduced cost of high-throughput sequencing applications will greatly facilitate identification of the genomic binding patterns of mammalian KRAB-ZFPs within the next few years. Importantly, these data will also help to improve tandem-ZFP prediction tools. However, the similarity between KRAB-ZFPs makes the generation of specific antibodies extremely challenging. Overexpression of epitope-tagged KRAB-ZFPs might partially overcome this problem [7,9]. However, this approach can be misleading because overexpression of KRAB-ZFPs can lead to the identification of binding sites Genomic PBS pro (n=149) Top-scored predicted motifs in genome (n=500)

A B
Predicted ZFP809 binding motifs Fig. 6 Comparison of predicted ZFP809 target motifs and experimentally identified target sequence of ZFP809. a The ZFP809 protein sequence was used to infer predicted target motifs using an expanded linear support vector machine (SVM) (top) or a polynominal SVM (middle) [27] and a prediction tool based on bacterial one-hybrid screens and ChIP-seq data (bottom) [9]. The canonical PBS pro sequence is shown below. b One of the predicted motifs [9] was used to screen matching genomic sites of the mouse genome using the software tool FIMO [199]. The 500 top-scored sites and the 149 genomic PBS pro were screened for ZFP809 enrichment by NGS.plot [200] and published FLAG-ZFP809 ChIP-seq data in murine embryonic carcinoma cells [7] that would not be bound by the protein expressed at normal levels [7]. It is also unknown whether ChIP-seq with overexpressed KRAB-ZFPs in a certain cell type will indicate the binding sites in other cells or tissues, since possibly cell-specific posttranslational modifications can influence the DNA binding properties of KRAB-ZFPs [191]. KRAB-ZFPs have evolved through natural section to bind precisely to sequences that needed to be transcriptionally repressed (e.g. ERVs). Nevertheless, binding to imperfect target sites is likely to occur. We speculate that such binding is kept to a minimum by selectively expressing KRAB-ZFPs at relatively low levels, ensuring that only the preferred binding sites are occupied. However, it is possible that, by regulating KRAB-ZFP expression levels, different cell types can determine the number of functionally relevant KRAB-ZFP binding sites. A powerful strategy to circumvent KRAB-ZFP overexpression in genome-wide binding assays will be tagging of endogenous KRAB-ZFP genes with epitopes using the CRISPR/Cas9 system.
Ultimately, loss-of-function studies will be required to determine the functional roles of KRAB-ZFPs. However, the repetitive and clustered characteristics of KRAB-ZFP genes-especially of evolutionarily young ones that have recently duplicated-remain a major obstacle. Conventional gene targeting by homologous recombination as well as RNAi-and CRISPR/Cas9-based approaches depend on unique sequences to ensure that only the right target is affected by knockdown or knockout. Alternatively, gene-trap collections of ESC clones might be used to generate KRAB-ZFPs knockout mice. However, mapping of gene-trap insertions at repetitive genes is problematic as well.
Thus new strategies will be needed to test the biological requirements of individual KRAB-ZFPs. We recommend targeting the rather unique regions between KRAB domains and ZNFs with CRISPR/Cas9 guide RNAs to mutate KRAB-ZFP genes. However, in the case of recently duplicated KRAB-ZFP genes, even these regions might be too similar for specific targeting by CRISPR guide-RNAS or siRNAs. Since recently duplicated KRAB-ZFP genes are usually located in close proximity, one might consider genetic depletion of several KRAB-ZFP genes within a cluster at once by inserting loxP sites in the unique flanking regions using CRISPR/ Cas9 or conventional gene targeting, followed by Cremediated recombination. Stepwise rescue experiments and ChIP-seq with epitope-tagged proteins may be used to assign the observed phenotypes in such KRAB-ZFP cluster knockout cells or animals to a single gene.
Once a larger number of KRAB-ZFPs have been characterized in detail, several important questions about the function and evolution of KRAB-ZFPs can be addressed.
For instance, little is known about the time it takes for KRAB-ZFPs to evolve to bind newly emerged sequences such as ERVs. It is also somewhat puzzling how KRAB-ZFPs can keep up with active and therefore mutating ERVs. One possibility is that KRAB-ZFPs (such as ZFP809) that bind to retroviral sequences that are at least partially conserved in various ERV groups are preferably selected as repressors of parasitic elements. Alternatively, KRAB-ZFPs might primarily repress ERVs and other retrotransposons that have already lost the ability to replicate, whereas other repression mechanisms act on newly emerged active ERE families, as previously suggested [10]. Nevertheless, insights could be gained by experimental testing how long (how many mutations) it takes for a KRAB-ZFP to start binding to a new ERE. One possibility would be to use a target that is only weakly bound by a KRAB-ZFP and test a large number of KRAB-ZFP mutants in a high-throughput screen to test how many mutations it takes to improve binding to the new target. Such an assay would also be useful to improve the target specificity of artificially engineered KRAB-ZFPs. In the long term, such optimized engineered KRAB-ZFPs might be used to suppress transposable elements and genes that cause disease in patients. Furthermore, replacing the KRAB domain with activating or other functional domains will allow us to tightly control expression of mobile DNA and regular genes. This approach might be used to induce transcription of ERE-linked genes to boost stem cell pluripotency or help differentiation into certain tissues.
The current KRAB-ZFP sets in mammals are likely a mixture of KRAB-ZFPs that are under purifying or positive selection and KRAB-ZFPs that arose by recent gene duplications and subsequent mutations but remain nonessential for the host. These KRAB-ZFPs will become pseudogenes and eventually disappear by genetic drift. Additionally, one might expect that duplicating and mutating KRAB-ZFPs sometimes reduce host fitness and are therefore rapidly removed by negative selection. How many useless or harmful KRAB-ZFPs are necessary before a beneficial one evolves? This question could be addressed by analyzing the genomes of very closely related species (e.g., wild mice and domestic mouse strains). One might even expect that there are differences in the number of KRAB-ZFPs between individuals of the same species. More importantly, germ-line mutations and segmental duplications of KRAB-ZFPs might influence human development or cause disorders and disease. Indeed, the KRAB-ZFP ZNF568 exist as three different alleles in humans, and a correlation between these alleles and the brain size of newborns was reported [192].
According to a recent study, natural occurring nonsynonymous single nucleotide polymorphisms (SNPs) at "specificity residues" of human tandem-ZFPs are rare, indicating that ZNF mutations that change the binding specificity of tandem-ZFPs are rapidly removed from the population by negative selection [193]. Furthermore, these rare SNPs generally do not correlate with altered gene expression profiles [193]. However, the highly repetitive nature of tandem-ZFP genes, especially at the ZNF coding regions, makes SNP calling extremely challenging and might have resulted in an underestimation of such polymorphisms. Moreover, conventional expression databases do not report expression levels of repetitive elements and SNP effects on ERE expression might therefore have remained undetected. Although a reanalysis of existing RNA-seq data might reveal ERE repression deficiencies associated with mutations in tandem-ZFP genes, many RNA-seq studies are still based on short (36 bp) read sequencing, which makes it problematic to assign reads to individual ERE copies and therefore hinders accurate quantification of ERE expression. Furthermore, the study excluded frameshift mutations and did not try to identify the loss or duplication of tandem-ZFPs in individuals [193]. Importantly, tandem-ZFP clusters on human chromosome 19 have been associated with unusually high copy number variation [150]. Although it will be difficult to identify events such as duplications of a single KRAB-ZFP by genome analysis, a thorough analysis of high coverage genome sequencing data and 100 bp paired-end RNAseq data might lead to the identification of physiologically relevant KRAB-ZFP polymorphisms in humans.
Little is known about how a loss-of-function mutation of a single ERE-targeting KRAB-ZFP would affect the host organism. While reactivation of a replication-competent ERV or a high-copy retrotransposon could have immediate deleterious consequences for the host, de-repression of non-autonomous ERVs or low copy transposons might only have subtle effects. In support of this idea, the drastic upregulation of a small subset of non-autonomous VL30 elements in ZFP809 knockout mice did not seem to impair their health or fitness [7]. Nevertheless, these mice were not monitored for more than two generations; deleterious effects might have emerged in later generations. Another possibility is that potentially hazardous ZFP809repressed ERVs are polymorphic among mouse strains and simply not present in the strain that was used in this study. Furthermore, ERV reactivation caused by KRAB-ZFP deletion might impair the host only under certain circumstances, such as physiological stress or during pathogen infection.
Interestingly, it seems that SVA and SVA-related LAVA elements have expanded in some primates that lack ZNF91 [8,[194][195][196], suggesting that ZNF91 prevents genomic SVA amplification. Moreover, polymorphic human KRAB-ZFPs may allow for mobilization of DNA in a certain haplogroup/population. Indeed, certain island populations carry SVA insertions [87] that may result from ZNF91 mutations. Future work on KRAB-ZFP knockout mouse models and genome-wide association studies of human KRAB-ZFP polymorphisms will reveal how a failure of the KRAB-ZFP ERE repression system impacts fitness and health of a host.

Conclusions
It is increasingly evident that transposable elements have a profound impact on mammals. Therefore, understanding the factors that keep these elements under control is of high importance for both basic and applied medical research. The recent evidence summarized in this review strongly supports the hypothesis that a large fraction of KRAB-ZFPs evolved to bind and possibly repress mobile parasitic DNA in mammals. However, since only a very small number of KRAB-ZFPs have been thoroughly investigated, the true spectrum of KRAB-ZFP functions cannot be anticipated yet. Without doubt, future research will yield exciting and unexpected insights into this enigmatic protein family.

Competing interests
The authors declare that they have no competing interests.
Authors' contributions GW wrote the manuscript with the support of DG and under the guidance of TSM. All authors read and approved the final manuscript.