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

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

Retrotransposons encompass half of the human genome and contribute to the formation of heterochromatin, which provides nuclear structure and regulates gene expression. Here, we asked if the human silencing hub (HUSH) complex is necessary to silence retrotransposons and whether it collaborates with TRIM28 and the chromatin remodeler ATRX at specific genomic loci. We show that the HUSH complex contributes to de novo repression and DNA methylation of an SVA retrotransposon reporter. By using naıve versus primed mouse pluripotent stem cells, we reveal a critical role for the HUSH complex in naıve cells, implicating it in programming epigenetic marks in development. Although the HUSH component FAM208A binds to endogenous retroviruses (ERVs) and long interspersed element-1s (LINE-1s or L1s), it is mainly required to repress evolutionarily young L1s (mouse-specific lineages <5 million years old). TRIM28, in contrast, is necessary to repress both ERVs and young L1s. Genes co-repressed by TRIM28 and FAM208A are evolutionarily young, or exhibit tissue-specific expression, are enriched in young L1s, and display evidence for regulation through LTR promoters. Finally, we demonstrate that the HUSH complex is also required to repress L1 elements in human cells. Overall, these data indicate that the HUSH complex and TRIM28 co-repress young retrotransposons and new genes rewired by retrotransposon noncoding DNA.
[Supplemental material is available for this article.] Although <2% of DNA sequence in the human genome codes for proteins, the vast majority plays an enigmatic role and has thus been referred to as genomic dark matter (Diederichs et al. 2016). However, this extra DNA serves a purpose: First, it contains regulatory elements that control when and where genes are expressed, a role only now being realized (Sanjana et al. 2016). Second, it is involved in building heterochromatin, for example at the nuclear periphery (Lemaitre and Bickmore 2015). Little is understood about how heterochromatin is formed, but its content is dominated by retrotransposons, which contribute to its establishment from plants to mammals (Lippman et al. 2004;Matsui et al. 2010). Retrotransposons replicate through an RNA intermediate, which has allowed them and their regulatory sequences to accumulate and coevolve with their hosts (Robbez-Masson and Rowe 2015; Thompson et al. 2016).
The human silencing hub (HUSH) complex, composed of FAM208A (also known as TASOR), MPHOSPH8 (also known as mpp8), and PPHLN1 (periphilin 1) is recruited to genomic loci rich in H3K9me3 (Brummelkamp and van Steensel 2015;Tchasovnikarova et al. 2015;Timms et al. 2016) and interacts with SETDB1 and MORC2 (Tchasovnikarova et al. 2017). The HUSH complex mediates position-effect variegation at reporter constructs that are integrated into silent chromatin (Tchasovnikarova et al. 2015), and depletion of HUSH components reduces H3K9me3 and alters transcription (Timms et al. 2016). It is un-known whether the HUSH complex is required for the repression of retrotransposons.
We asked here if the HUSH complex, like TRIM28, is necessary for retrotransposon repression and whether it cooperates with TRIM28 and ATRX at specific genomic loci.

Results
The HUSH complex contributes to repression of an SVA retrotransposon reporter We used a human retrotransposon reporter (Fig. 1A) that is repressed through ZNF91 binding to an SVA-type D variable number tandem repeat (SVA VNTR) sequence (Jacobs et al. 2014). We found that the SVA reporter was repressed (3.3×) in POU5F1-positive human embryonal NTERA-2 cells, which naturally express ZNF91, but not in 293T cells (Fig. 1B). Reporter repression could be engineered in POU5F1-expressing mouse embryonic stem cells (mESCs) by transfection of the cognate ZNF91, but not of a control KZNF, ZNF93 as expected (Supplemental Fig. S1A;Jacobs et al. 2014), and even in 293T cells (Supplemental Fig. S1B). We used 293T cells for subsequent reporter assays because of their amenity to genetic manipulation.
We validated that reporter repression (6.7×) (Fig. 1C) was TRIM28-dependent by using TRIM28 knockout and complemented 293T cells. Having set up this system, we depleted the HUSH complex using shRNA before introduction of the reporter gene (Fig. 1D). All three HUSH components contributed to SVA repression ( Fig. 1E), with the most striking effect observed for MPHOSPH8 depletion (14× derepression). Although these results were obtained with single hairpins, we validated the role of Mphosph8 with an independent hairpin in mESCs (Supplemental Fig. S1C). Considering that the HUSH complex is thought to function similarly to CBX proteins (Brummelkamp and van Steensel 2015), we depleted all three CBX family members, which revealed that they too all contributed to repression in 293Ts and HeLa cells (Supplemental Fig. S1D). Finally, we identified a role for H3F3A (4× derepression) and its ATRX-DAXX chaperone complex in establishing SVA repression ( De novo DNA methylation of the SVA reporter depends on the HUSH complex We reasoned that the HUSH complex may be targeted to the repressed SVA reporter through the chromodomain of MPHOSPH8 that interacts with H3K9me3 (Kokura et al. 2010); indeed, we found H3K9me3 to be enriched (2.4×) on the repressed reporter (Supplemental Fig. S1E). Of note, plasmids are chromatinized and subject to H3K9me3 (Barde et al. 2009). We also detected cognate ZNF-specific de novo DNA methylation (23.5%) (Fig. 1F), as observed before for retroviral reporters (Wolf and Goff 2009;Rowe et al. 2013a). DNA methylation was not necessary for repression but contributed to it, since we detected . Western blot for TRIM28 using PCNA as a loading control (insets from the same blot).
(D) 293T cells were transduced with shRNA vectors against epigenetic factors and puromycin-selected before reporter assays. (E) Following the assay in D, data for ZNF91 were normalized to ZNF93 (left). Unpaired t-tests were used to assess differences between the shControl (transduced with the same vector lacking a hairpin that was puromycin-selected in parallel) and the target shRNA. Experiments were repeated at least three times for each candidate, and representative data are shown. Two-tailed unpaired t-tests were done: ( * ) P = <0.05; ( * * ) P = <0.01; ( * * * ) P = <0.001. qRT-PCR was used to assess the knockdown efficiency of each mRNA ( 4.5× SVA repression in Dnmt knockout mESCs, compared to 9× repression in wild-type mESCs (Supplemental Fig. S1F). This suggested HUSH may be necessary for handover of H3K9me3 to DNA methylation and in support of this, depletion of HUSH components coincided with a decrease (5.8×) in de novo DNA methylation at the SV40 promoter along with a reduction (up to 2.7×) in SVA reporter repression ( Fig. 1G; Supplemental Fig. S1G). In sum, these data on the nonintegrated reporter suggest that HUSH may be required for the maintenance of retrotransposon repression.
The HUSH complex is critical for endogenous retrotransposon repression in naıve pluripotent cells We reasoned that the HUSH complex may exert its greatest impact in naïve pluripotent cells, in which chromatin awaits stable epigenetic programming (Ying et al. 2008;Ficz et al. 2013). We depleted epigenetic factors in mESCs ( Fig.  2A) and compared mixed population to naïve cultures, the latter in which we verified enhanced NANOG expression (Fig.  2B). Depletion of Mphosph8, as well as Atrx and Trim28, was sufficient to reactivate retrotransposons in two different strains of serum-cultured mESCs (Supplemental Fig. S2A). However, although Trim28 or Atrx depletion mainly affected intracisternal A-particle (IAP) elements, Mphosph8 depletion mainly affected L1s. Parallel culture of J1 ESCs in serum versus 2i conditions consistently led to more pronounced reactivation in the naïve cells ( TRIM28 and FAM208A co-repress a set of protein-coding genes We sought to identify genomic sites where TRIM28 collaborates with FAM208A, SETDB1, and ATRX and therefore performed mRNA-sequencing (  Table  S4). These three groups of TRIM28-repressed genes are likely direct rather than indirect targets, because up to 81% overlapped TRIM28 peaks (TRIM28-FAM208A and TRIM28-SETDB1 gene sets) and up to 77% overlapped H3K9me3 (all three gene sets), compared to randomly-selected gene groups ("Random") ( One representative experiment of three is shown. Atrx was not examined in the third experiment, excluding it from statistical analyses. Two-tailed paired t-tests were done for 2i + LIF samples. (D) Western blot for IAP GAG p73 using a rabbit IAP GAG antibody or PCNA as control in 2i + LIF J1 ESCs. The antibody detects p73 as well as GAG-POL and GAG cleavage products, including p41, representing partially processed GAG. Samples were re-run on a second gel and reblotted for L1 ORF1 protein (40 kDa    groups were enriched in protein-coding genes (Fig. 3D), and follow-up of the TRIM28-FAM208A genes revealed them to be functionally related and involved in developmental pathways (Fig.  3E), unlike random genes. Of note, we verified that up-regulation was detectable across the length of the transcripts within this group (Supplemental Fig. S3C).
TRIM28 and FAM208A exert nonredundant roles in repressing young L1 elements TRIM28 regulates genes through binding repeats (Rowe et al. 2013b;Hummel et al. 2017). We therefore asked if these sets of genes were enriched in repeats and divided repeats by size ( Fig.  3F; Supplemental Fig. S3D). TRIM28-SETDB1 genes were enriched for all repeats here (9× for LINEs >5 kb), whereas TRIM28-ATRX genes were enriched for full-length LTRs (4×); most interestingly, TRIM28-FAM208A genes were enriched for LINE elements of any size but particularly those >5 kb (7×) (Fig. 3F; Supplemental Fig.  S3D) and for satellites (2×). In fact, 81.72% of TRIM28-FAM208A genes contained an L1 or satellite or both within 20 kb (only TRIM28-regulated L1s were included) , and L1s were often present in arrays (Fig. 3G). We verified that TRIM28-FAM208A genes contained significantly more TRIM28-L1s than random genes (e.g., 2.3× for L1Md_T) (Supplemental Fig. S3E). Because TRIM28-FAM208A genes were most enriched in LINEs >5 kb ( Fig. 3F; Supplemental Fig. S3D), we selected all of these (within 20 kb), ordered them by their evolutionary age (Sookdeo et al. 2013), and found 98% to be <5 million years old and absent from rat genomes. We also ordered them by mean divergence (Supplemental Fig. S3F) and found them all to be <10% diverged from consensus sequences and the L1Md_A integrants to be the youngest by this method (mean divergence 1.53%).
In a complimentary approach, we mapped mRNA-sequencing reads to Repbase, which showed repeats to be overexpressed in all four treatment groups ( Fig. 3I; Supplemental Table S5). Scoring the top five repeats derepressed in each treatment group showed that IAP elements are co-repressed by TRIM28, SETDB1, and ATRX, whereas FAM208A mainly represses young L1s (from the TF and GF families) (Supplemental Fig. S3G). In total, 25 families of repeats were co-repressed by TRIM28 and FAM208A (Fig.  3J). This included 13 L1 families that we classified by age (Sookdeo et al. 2013) and found 77% to be <3 million years old (Fig. 3K). The rest of the TRIM28-FAM208A co-repressed repeats fell into the ERV or satellite classes with satellites most highly derepressed (14×) (Fig. 3L,M). Overall, these data suggest that TRIM28 and FAM208A co-repress young L1s. Of note, we did not assess polymorphic or de novo L1 insertions.

TRIM28-FAM208A coregulated genes are enriched in tissue-specific and new genes
LTRs provide genetic material (promoters, enhancers, and first exons) to create new genes or new expression patterns (Franke et al. 2017), and active L1s can create new genes through retroposing cellular mRNAs (Carelli et al. 2016). We asked if TRIM28-FAM208A genes (Supplemental Table S6) were enriched in new genes, (which we define here as mouse-specific) and tissue-specific genes. We first found 41% of genes were not conserved across placental mammals and were mouse-specific, compared to 14% of random genes (Fig.  4A, left); focusing only on the protein-coding genes and their last common ancestor as a measure of their evolutionary age also revealed TRIM28-FAM208A genes to be enriched in mouse-specific genes (11.8% compared to 0.9% of all genes in the mouse genome) (Fig. 4A, right). Forty percent of TRIM28-FAM208A genes exhibited tissue-specific expression patterns (Fig. 4B), whereas 44% had an unknown expression pattern and this group was enriched for new genes (Fig. 4B). We verified that example loci of these new or tissue-specific genes were associated with arrays of young L1s or ERVs and epigenetic regulation ( Fig. 4C; Supplemental Fig. S4). Finally, we observed that LTRs 3 kb upstream of TRIM28-FAM208A genes were biased to reside in a sense orientation (69%) suggesting they may function as promoters (Fig. 4D).

FAM208A binds primarily to ERVs and L1 elements
We asked if FAM208A mainly regulates young L1s because it binds selectively to young L1s. We addressed this by performing ChIPseq using an antibody recognizing mouse FAM208A and mapped reads to Repbase. FAM208A binds a range of retrotransposons, primarily ERVs and L1s ( Fig. 5A; Supplemental Table S7), coating their 3 ′ halves (Supplemental Fig. S5A). Repbase L1s bound by FAM208A (22 families with an enrichment of more than 4×) were inactive families lacking full-length copies (Sookdeo et al. 2013) mostly older than 13 million years, because they were present before the mouse-rat split (Fig. 5B). We therefore also mapped reads to the genome to see if we could detect FAM208A binding to young L1s. We found 1045 FAM208A peaks ( Fig. 5C; Supplemental Table S8), which clustered together (75% of ChIP-seq peaks were within 50 kb of another FAM208A peak) (Fig. 5D), suggesting FAM208A binding spreads with 34% of FAM208A peaks overlapping H3K9me3 (Fig. 5E; Rowe et al. 2013b). We could detect FAM208A binding to young L1s: Of 194 peaks targeting L1s, 6% of them were young L1s, whereas the rest were inactive L1s ( Fig.  5F; for an example of a bound young L1, see Supplemental Fig.  S5B). FAM208A may spread to young L1s mainly from its tethering to inactive L1s and ERVs. In support of this, we found that 61% and 87% of TRIM28-FAM208A genes (that are enriched for young L1s) (Fig. 3H) contained a FAM208A-bound L1 or an ERVK within 20 kb, respectively (Fig. 5G).

Discussion
A key question has been whether the HUSH complex participates in retrotransposon repression and whether it collaborates with TRIM28 and its cofactors. Here, we show that the HUSH complex and TRIM28 exert nonredundant roles at evolutionarily young L1s of <5 million years old in naïve pluripotent cells. These young elements are likely prone to derepression because they bind TRIM28 and FAM208A only weakly. In the case of TRIM28, weak binding to young retrotransposons is known to result from them escaping KZNF recognition due to sequence divergence (Jacobs et al. 2014). In contrast, TRIM28 and FAM208A are strongly associated with ERVs, which recruit dense heterochromatin and at which only TRIM28 but not HUSH is required for repression. We found FAM208A binding spreads through chromatin, consistent with the ability of HUSH to mediate position-effect variegation (Tchasovnikarova et al. 2015). Our data reveal that TRIM28 and FAM208A coregulate new genes and tissue-specific genes, which are associated with young L1s. This suggests that regions of the genome that are enriched in young L1s and leaky heterochromatin/ part euchromatin may be hot spots for the evolution of new genes and new regulation of existing genes. Such genes may hijack incomplete epigenetic repression to gain tissue-specific expression (for a summary model, see Fig. 7). TRIM28 and FAM208A cooption into these gene regulatory networks, therefore, appears to be a byproduct of their leaky regulation of young retrotransposons.
In summary, this work illustrates the complexity of epigenetic repression pathways that have evolved to regulate cellular genes through retrotransposon sequences. Future work focused on specific retrotransposon integrants will enable us to understand how these elements have been coerced to roles in gene regulation and chromatin organization in mammalian cells.

Cell culture
Human embryonal NTERA-2 cells (NT2/D1, a kind gift from Peter Andrews, University of Sheffield) were grown in Dulbecco's Modified Eagle's Medium (DMEM, Gibco) high glucose with 2 mM L-Glutamine, 10% fetal calf serum (FCS), and 1% Penicillin/ Streptomycin (P/S). They were split in half by cell scraping. 293T cells were grown in standard DMEM + 10% FCS and P/S. J1 ESCs (129S4/SvJae) or their derived triple knockout (TKO) cells, which are knockout for Dnmt3a, Dnmt3b, and Dnmt1 (from M. Okano) were used where stated, or as an independent mESC line, ES3 ESCs (Rowe et al. 2010). All mESCs were cultured on gelatin-coated  Table Browser). For the TRIM28-FAM208A genes, the nonconserved ones were verified to be mouse-specific using the Ensembl GeneTree. Fisher's exact test one-sided P-value = 2.204 × 10 −5 . (Right) Only protein-coding TRIM28-FAM208A genes (n = 68) were selected and their Last Common Ancestor extracted from the Ensembl database (version 90) using R version 3.3.1, compared to all genes in the mouse genome. Fisher's exact tests on 2 × 2 tables were significant for Mus musculus (P = 1.58 × 10 −7 ).  Figure S4. (D) 3 kb regions were identified upstream of each TRIM28-FAM208A coregulated gene or the random genes, and the orientation of all LTRs in these regions was assessed using the UCSC Table Browser. In the TRIM28-FAM208A group, LTRs were shown to be biased to be in a sense orientation (P = 0.005692, Fisher's exact one-sided test).
plates (0.2%), either in standard media as previously described (Rowe et al. 2013a), or in media containing two small-molecule inhibitors of MEK and GSK3 plus LIF (2i + LIF), as described (Ying et al. 2008;Ficz et al. 2013). Cells were split every other day using accutase (Gibco, A11105-01) and cultured in 2i media for 7 d before transduction.

Luciferase assays
Dual luciferase assays were performed at ratios detailed elsewhere (Jacobs et al. 2014). 293T cells were plated at 5 × 10 4 cells per well in a 24-well plate and transfected with 200 ng of KZNF plasmid, 20 ng of luciferase reporter plasmid, and 2 ng of pRT-TK_Renilla luciferase-encoding control plasmid using 1.5 µL of FuGENE 6 (Promega) per well in triplicate wells (or for ESCs, we used 500 ng KZNF, 50 ng luciferase reporter, and 5 ng Renilla luciferase plasmid in 12-well plates). Forty-eight hours post-transfection, cells were lysed, and luciferase was measured using the Dual Luciferase assay kit (Promega, E1910), an opti-plate, and a GloMax 96 microplate Luminometer (Promega) using the Dual Glow program. Raw luciferase values were normalized to Renilla luciferase values to control transfection efficiency.

RNA extraction and quantification
Total RNA was extracted using an RNeasy mini kit (Qiagen), treated with DNase (Ambion AM1907), and 500 ng was reverse transcribed using SuperScript II (Thermo Fisher Scientific) and random primers following the manufacturer's instructions. Samples were run on an ABI 7500 Real Time PCR System (Applied Biosystems) using SYBR green Fast PCR mastermix (Life Technologies). CT values for test genes were normalized to B2M and GAPDH for human genes or Cox6a1 and Gapdh for mouse genes using the −ΔΔC t method to calculate fold change. See Supplemental Table S1 for primers.

Western blotting
293T cells were washed with PBS and lysed in 2X Sample Buffer (Invitrogen) supplemented with 5% β-mercaptoethanol (Sigma). ESCs were washed 2× in ice-cold PBS and resuspended in radioimmunoprecipitation buffer with protease inhibitors added (Roche 11836170001) and quantified using the BCA Assay (Pierce). Ten micrograms was loaded per well on mini Tris/glycine gels, run in Tris/glycine/SDS buffer and mini-PROTEAN tanks (Bio-Rad) or on precast Bis-Tris gels (4%-12% or 8% gels), run in MOPs SDS buffer, followed by wet transfers to PVDF membranes. Primary antibodies are in Supplemental Table  S2. Secondary antibodies were horseradish peroxidase-conjugated (GE healthcare), and membranes were developed using ECL kits (Amersham).

mRNA-sequencing
Mouse 2i + LIF cultured J1 ESCs treated with different shRNA vectors were used for mRNA-sequencing. Cells were cultured for 5 d

A B
C E F G D Figure 5. FAM208A binds primarily to ERVs and L1 elements. (A) Reads from TI and FAM208A IP samples were mapped to rodent Repbase. Duplicates were averaged and RPKM ratios calculated between TIs and IPs. Repeats were selected giving ≥fourfold enrichment in the IPs. The ERV2 class includes ERVK elements from ETN and IAP families. (B) All L1 elements from A are displayed here with family name and age. (C ) After mapping reads to mm10, 1045 peaks were identified (present in both duplicates and not in the TIs). (D) Peaks from C were sorted into those that clustered by their presence within 50 kb of a second FAM208A peak. (E) Intersection of FAM208A peaks with H3K9me3 peaks (either any overlap or 80% overlap considered). (F ) FAM208A peaks from C overlap with young versus inactive L1s. (G) TRIM28-FAM208A genes were assessed for the percentage that contain either an inactive L1 (from the families L1_Rod, L1MB7, L1_Mur1, L1_Mur2, L1_Mur3, Lx8, Lx9, and Lx10) or an ERVK within 20 kb.
following puromycin selection before RNA extraction, and samples from three independent experiments were used. See Supplemental Methods for further details.

Repbase analysis
Reads were mapped to rodent Repbase (https://www.girinst.org/ downloads/) and the latest release downloaded (Repbase 20.06 used here). The SAMtools v.1.19 idxstats utility (Li et al. 2009) was used to extract the number of mapped reads per repeat, that were inputed into the R package DESeq2 (https://bioconductor. org/packages/3.2/bioc/html/DESeq2.html) to identify differentially expressed repeats between samples depleted of epigenetic modifiers and controls, as previously described (Love et al. 2014). P-values were adjusted for multiple testing with the Benjamini-Hochberg false discovery rate (FDR) procedure. Chromatin Immunoprecipitation (ChIP) 293T cells were harvested using trypsin, whereas NTERA-2 cells and 2i + LIF grown mESCs were harvested using accutase, and chromatin was cross-linked, quenched, and prepared as described (Rowe et al. 2013b), except that sonication was performed on a Bioruptor (Diagenode). See Supplemental Methods for further details.

DNA methylation analysis
DNA was purified with a DNeasy Blood and Tissue Kit (Qiagen), and 1 µg DNA was used for bisulfite conversion using an EpiTect Bisulfite Kit (Qiagen). Four microliters converted DNA was amplified by PCR using primers in Supplemental Table S1, which are from Rowe et al. (2013a) or were designed using http://urogene. org/methprimer/, and PCR products were cloned using the TOPO TA-Cloning Kit (Thermo Fisher Scientific), and at least 10 colonies were sent for sequencing using the T7P primer. Results were analyzed using the QUMA online tool (http://quma.cdb. riken.jp) from the Riken Institute.

Statistical analysis
All data in the figures are presented as the standard deviation (where there are three or more samples) or by standard error of the mean (SEM) and assessed by unpaired or paired two-tailed Student t-tests (see figure legends for details). A P-value of <0.05 was considered statistically significant ( * * * P < 0.001; * * P < 0.01; * P < 0.05).

Data access
mRNA-sequencing and ChIP-sequencing data and processed files from this study have been submitted to the NCBI Gene Expression Omnibus (GEO; http://www. ncbi.nlm.nih.gov/geo/) under accession number GSE107840, and processed files are also included in Supplemental Material. FAM208A also binds to ERVs and is known to interact with H3K9me3 through MPHOSPH8. FAM208A binding spreads through chromatin and overlaps H3K9me3, suggesting HUSH uses H3K9me3 as a platform on which to spread. TRIM28 is required to repress ERVs but FAM208A is not, likely because it is redundant at these sites of dense H3K9me3. Young L1s, in contrast, reside in "leaky heterochromatin" or part euchromatin, which exhibits weak TRIM28 and FAM208A binding and low levels of H3K9me3 and DNA methylation. Both TRIM28 and FAM208A exert nonredundant roles at young L1s. These sites are also rich in new and tissue-specific genes and are flanked by upstream sense LTRs. This suggests that genes may hijack repeats and incomplete epigenetic repression to rewire their expression patterns.