- Open Access
The UCSC repeat browser allows discovery and visualization of evolutionary conflict across repeat families
Mobile DNA volume 11, Article number: 13 (2020)
Nearly half the human genome consists of repeat elements, most of which are retrotransposons, and many of which play important biological roles. However repeat elements pose several unique challenges to current bioinformatic analyses and visualization tools, as short repeat sequences can map to multiple genomic loci resulting in their misclassification and misinterpretation. In fact, sequence data mapping to repeat elements are often discarded from analysis pipelines. Therefore, there is a continued need for standardized tools and techniques to interpret genomic data of repeats.
We present the UCSC Repeat Browser, which consists of a complete set of human repeat reference sequences derived from annotations made by the commonly used program RepeatMasker. The UCSC Repeat Browser also provides an alignment from the human genome to these references, uses it to map the standard human genome annotation tracks, and presents all of them as a comprehensive interface to facilitate work with repetitive elements. It also provides processed tracks of multiple publicly available datasets of particular interest to the repeat community, including ChIP-seq datasets for KRAB Zinc Finger Proteins (KZNFs) – a family of proteins known to bind and repress certain classes of repeats. We used the UCSC Repeat Browser in combination with these datasets, as well as RepeatMasker annotations in several non-human primates, to trace the independent trajectories of species-specific evolutionary battles between LINE 1 retroelements and their repressors. Furthermore, we document at https://repeatbrowser.ucsc.edu how researchers can map their own human genome annotations to these reference repeat sequences.
The UCSC Repeat Browser allows easy and intuitive visualization of genomic data on consensus repeat elements, circumventing the problem of multi-mapping, in which sequencing reads of repeat elements map to multiple locations on the human genome. By developing a reference consensus, multiple datasets and annotation tracks can easily be overlaid to reveal complex evolutionary histories of repeats in a single interactive window. Specifically, we use this approach to retrace the history of several primate specific LINE-1 families across apes, and discover several species-specific routes of evolution that correlate with the emergence and binding of KZNFs.
Transposable elements are significant drivers of eukaryotic genome evolution. In humans and other primates, transposons constitute nearly half the genome; the majority of these repeat elements are retrotransposons, although some DNA transposons are also present. Despite the high repeat content of the human genome, many genomic analyses struggle to deal with these regions as sequencing reads can often be assigned nearly equally well to multiple regions in the genome. Masking or filtering these reads is often considered a “conservative” approach in that it avoids mis-assigning the genomic location of a read, but it prevents the discovery of potentially important biological functions of repeat elements . Indeed, many repeats already have established roles in important biological processes, complex behavioral phenotypes, and disease [2,3,4,5].
One of the major challenges in proper repeat-analysis is establishing a set of standardized sequences, nomenclature and annotation sets that can be universally understood by the scientific community. The most commonly used databases and tools to study repeats are Repbase  and RepeatMasker . Repbase began as a hand-curated list in 1992 of 53 prototypic repeat sequences identified in the human genome . By 2015, it contained more than 38,000 sequences in 134 species , making curation and comprehension of each repeat family a daunting challenge. RepeatMasker is a program that screens DNA (e.g. a newly sequenced genome) for repeat elements by filtering, merging and joining human genome alignment matches based on a database of partial repeat sequences. The strategies are optimized for the different classes of repeats. For instance, some full length repeat elements such as LINE-1 elements are built after alignments for sequences of smaller subparts (e.g. 5′ and 3′ UTRs) are identified (Fig. 1a).
One option for the sequence database as the basis of these alignments is RepBase (RepBase RepeatMasker Edition). More recent versions of RepeatMasker utilize the Dfam database , an open database of hidden Markov model profiles of repeat elements, as RepBase now requires a commercial license. All sequences and RepeatMasker annotations described in the following are therefore based on Dfam.
Although a variety of tools and methods already exist to study repeats , tools to dynamically visualize genomic data utilizing existing annotation sets on repeats (e.g. protein coding regions, conservation with other sequences and the list of matches in the genome) are still needed. Generating and mapping to a consensus version of individual repeats has proven successful in illustrating novel biological features of transposon insertions, but has largely been limited to static visualizations on targeted elements of interest and specific families of these repeats [11,12,13].
Here we present the UCSC Repeat Browser, which simplifies analysis of genomic data on repeats by using consensus sequences for all human repeat element classifications within RepeatMasker. The Repeat Browser allows “lifting” of human genomic data to these consensuses where they can be analyzed in conjunction with a precomputed set of comprehensive annotations in an interactive genomic browser environment (Fig. 1). Further, we demonstrate the utility of the Repeat Browser in illustrating how a primate specific class of retrotransposons has evolved to delete portions of their own sequence that likely allow them to evade the binding of repressors.
Generating reference sequences for human repeats
We ran RepeatMasker using Dfam on the two most recent assemblies of the human genome (hg19 and hg38) and generated a list of repeat families (excluding the simple repeats) annotated on the genome. For families that had a Dfam equivalent (Table S1), we utilized the Dfam sequence as a consensus. For families that have no exact Dfam equivalent (Table S1), where the repeats are annotated by joining alignments of partial sequences, we generated our own consensuses.
To do so, we downloaded all nucleotide sequences and their annotations in the RepeatMasker annotation track on the UCSC Human Genome Browser (hg38). We observed that extremely long repeats tended to represent recombination or misannotation events and therefore removed the longest 1% of sequences in all classes. We then aligned the 50 longest remaining sequences of each class, as this produced a tractable number of sequences that allowed manual inspection of each alignment, and because insertions relative to the consensus are difficult to show on shorter sequences. For each repeat family, these fifty sequences were realigned with MUSCLE  to create a consensus sequence. Each of these consensus sequences was then stored as a “reference” in the Repeat Browser in a manner analogous to a single chromosome on the UCSC Human Genome Browser [15, 16]. Each alignment is provided as a link in a “consensus alignment” track for additional visual inspection by the user.
In addition, we added two consensus sequences manually, HERVH-full and HERVK-full, which represent full length reconstructions (internal regions + LTR) that have been shown to produce viral particles [17, 18]. We added these sequences as HERV-K and HERV-H are well studied ERVs with established biological roles, and users may wish to analyze these sequences in the context of a full ERV instead of the smaller subparts (ERV internal regions and LTRs) that RepeatMasker does not join in its final output. However the Repeat Browser also contains the Dfam consensuses for these elements (HERVH and HERVK, which correspond only to the internal regions of the ERVs as well as their corresponding LTRs) meaning that data mapped to any of these regions in the human genome corresponds to multiple consensus sequences in the Repeat Browser.
We also generated a track with existing Dfam annotations for the Dfam consensus sequences. For our custom generated consensuses, we aligned the RepeatMasker Peptide Library  (filtered only to include peptides derived from human repeats) with tBLASTn  and took the highest scoring hits to automatically annotate ORFs. We also ran Tandem Repeats Finder  on these consensus sequences to annotate tandem repeats on each consensus (e.g. variable nucleotide tandem repeat regions within repeats). The combination of these consensus sequences (Dfam and custom Repeat Browser sequences) serve as the “genome” for all repeat content and are collectively called “hg38reps”. The script to build these consensus sequences is available as part of the “buildSeqs” step in the Repeat Browser source code.
Building LiftOver chains between the genome and repeat browser
In order to map a coordinate set on the human genome to coordinates on our consensus sequences on the Repeat Browser (a process more generally known as “lifting” ), we first mapped every consensus to all annotated genomic instances (both hg19 and hg38) of the same type using BLAT  (Fig. 1a). From this process, we generated a coverage plot illustrating where the instances align to the consensus (Fig. 1b). For example, the primate-specific LINE-1 sub-family, L1PA5, shows the distribution expected for recently active LINE elements: most individual L1PA5 instances, are defective short 3′ truncations, with only a few near full length elements containing the 5′ portion. Therefore the 3′ end of the consensus is found relatively often across the human genome (Fig. 1b). Using these alignments between consensus and genome, we then generated liftOver chains that allow genomic data to be visualized on the Repeat Browser. The script to align the consensuses and build the liftOver chains files is provided as the “buildLiftOver” step in the Repeat Browser source code.
Lifting of genomic annotations
These liftOver chains allow us to map existing genome annotation sets from hg19 and hg38 to the Repeat Browser consensus sequences. We lifted several standard genomic annotations to the Repeat Browser (Table 1); for instance, lifting of human gene annotation sets (e.g. GENCODE ) allows visualization of genes which contain repeat sequence within them. For example, L1HS sequences (Fig. 2) that have been incorporated into protein coding genes tend to derive from the untranslated regions (UTRs) of the repeats, and also tend to incorporate into the UTRs of the protein coding genes. Conversely, ncRNAs tend to contain the 5’UTR of L1HS elements and can span longer portions of the consensus. The result of these procedures produces a fully annotated and interactive consensus sequence that requires minimal prior knowledge of the genomic organization of the repeat being analyzed and at the same time allows lifting of any genome annotation from either hg19 or hg38 available at online at repeatbrowser.ucsc.edu.
Mapping of existing genomic datasets
We also mapped genomic loci bound by histone-modifying enzymes from ENCODE datasets  as well as large-scale ChIP-seq collections KRAB Zinc Finger Proteins (KZNFs) [26, 27] to the Repeat Browser. KZNFs are particularly compelling factors as they engage in evolutionary “arms races” in which KZNFs evolve unique DNA binding properties to bind and repress retrotransposons [11, 28]. These retrotransposons then accumulate mutations that allow evasion of KZNF-mediated repression . In order to map this ChIP-seq data to the Repeat Browser, we first downloaded raw ChIP-seq reads from the Sequence Read Archive (SRA) , mapped them to the reference genome (hg19) using bowtie2  and called peaks using macs2  (Fig. 3a). After this standard genomic mapping and peak calling, we then took the summits of these peaks, extended them by 5 nt in both directions, and lifted them to the Repeat Browser consensus sequences. In essence, this approach leverages each repeat instance as a technical replicate, with the mapping to the consensus representing a combination of many genomic “replicates” (Fig. 3a) of DNA binding summits called on individual instances of a repeat family that individually produce a noisy set of mappings; however hundreds of them combined yield a clear overall signal, better identifying the actual binding site. We call this “summit of summits” (obtained by combining the summits on individual transposon instances into a single summit on the Repeat Browser consensus) the “meta-summit”. In order to determine these “meta-summits”, we used a peak caller (macs2) on the repeat consensus to generate a list of “meta-summits” which represent the most likely location of the DNA binding site for a specific DNA-binding factor. We then generated a track which summarizes these meta-summits for each consensus sequence allowing easy and quick determination of factors with correlated binding patterns (Fig. 3b; visualized on a full length sequence of Human Endogenous Retrovirus H, HERV-H). These meta-summits serve as a quick summary of which factors bind which consensus sequence. These meta-summits can then be investigated more deeply by examining the coverage and mappings of each individual factor which are provided as separate tracks.
Comparative analysis of L1PA elements
In order to demonstrate the power of the UCSC Repeat Browser, we studied the evolution of recent L1PA families. The L1PA lineage is a group of LINE-1 retrotransposon families specific to primates. These elements are fully autonomous, encoding proteins ORF1 (Gag in the RepeatMasker Peptide Library) and ORF2 (Pol in the RepeatMasker Peptide Library) responsible for reverse transcription and re-integration of the retrotransposon. L1PA families evolve in bursts; higher numbers (e.g. L1PA17) indicate ancient evolutionary origins, as evidenced by shared copies across species (Fig. 4a). Lower numbers indicate more recent activity and are derived from the older, higher number families (note L1PA1 is also known as L1HS, human-specific) . Although this nomenclature generally corresponds to speciation events on the phylogenetic tree of the hosts of L1PA retrotransposons, many families have had overlapping periods of activity meaning that the correspondence is not exact .
Comparison of primate repeat elements reveals a large number of gibbon specific L1PA4 elements
In order to trace the evolution of L1PAs in different species, we downloaded the complete sequences for every L1PA7 and younger L1PA family, as annotated in their UCSC Genome Browser RepeatMasker tracks, in rhesus macaque (rheMac10), gibbon (nomLeu3), orangutan (ponAbe3), chimp (panTro6), gorilla (gorGor5), bonobo (panPan2) and human (hg38). We further restricted our analysis to near full-length elements by filtering out elements less than 6000 nucleotides in length (in humans active L1 elements are ~ 6000 nt). Although most of the L1 elements in this analysis are not active, they serve as a genomic “fossil” record of once active elements that can be used to trace L1 evolution across species.
As expected, the number of elements in older families were largely similar amongst all species that shared a common ancestor when the retrotransposon was active: for instance, L1PA7, active prior to the emergence of the last common ancestor of all primates in this study, was found at a relatively constant amount in all genomes (Fig. 4b). On the other hand, human specific elements were found only (barring a few likely mis-annotations) in the human genome. In certain species (gibbon, orangutan and gorilla) instances of retrotransposon families that were active near their divergence from human, were present in much greater copy number than in human (Fig. 4b). Specifically, the number of L1PA4 elements was greater in gibbon then all other apes, while a similar pattern was seen for L1PA3 and orangutan, and L1PA2 and gorilla. These results are consistent with these primates having lineage specific expansion of these elements in a manner distinct from humans. Notably, bonobos had a markedly lower number of L1PA2 elements which may indicate stronger repression of these elements by a species-specific factor; however, the bonobo assembly was one of the older, short-read primate assemblies used in this study, and therefore the lack of L1PA2 elements may simply reflect greater difficulty in resolving these regions in the genome assembly. Note also that the UCSC RepeatMasker track for rheMac10 contains no annotated instances of L1PA5, but this simply reflects the fact that RepeatMasker taxonomy splits the L1PA5 family into L1_RS families that are rhesus-specific compared to the other primates in this study . These L1_RS instances are added to the L1PA5 count in Fig. 4b.
All apes display evidence of ZNF93 evasion in the 5’UTR of L1PA
In order to examine the selection pressures that might explain species-specific expansion and restriction of L1PA elements, we combined our primate L1PA analysis with the ChIP-seq data of KRAB Zinc Finger Proteins (KZNFs) on the Repeat Browser [26, 35]. KZNFs rapidly evolve in order to directly target retrotransposons and initiate transcriptional silencing of these elements. We previously demonstrated that a 129 bp deletion occurred and fixed in the L1PA3 subfamily (and subsequent lineages of L1PA) in order to evade repression mediated by ZNF93. In order to discover additional cases where a retrotransposon may have deleted a portion of itself to escape KZNF-mediated repression, we searched for L1 sequences with the following characteristics: 1) deletion events proximal to KZNF binding sites, and 2) increasing number of retrotransposon instances with that deletion (demonstrating increased retrotransposon activity). Comparisons of these events across primate species provides evidence for unique, species-specific mechanisms of escape.
In order to look for these signatures of L1PA families escaping repression, we used BLAT to align each individual full-length (> 6000 nt) primate L1PA of the same class instance to the human Repeat Browser consensus from the primate genomes under study. We then generated coverage tracks of these full-length elements mapped to the human consensus for each species and each L1PA family. The ZNF93-associated deletion is clearly visible as evidenced by a massive drop in coverage in the 129-bp region in human L1PA3 instances (Fig. 5a). This same drop in coverage is found in all great apes (orangutan, gorilla, bonobo, chimp, and human) confirming that this event occurred in a common ancestor. To explore how more this retroelement lineage may have evolved in more distant primates, we examined gibbon L1PA4 elements on the Repeat Browser and found no evidence of the large 129-bp deletion seen in great apes. Instead, some gibbon L1PA4 elements contain a small 20 bp deletion - at the base of the ZNF93 peak (Fig. 5b). As humans and other great apes do not contain this deletion, we believe that this mutation first occurred in gibbon L1PA4 elements after the human-gibbon divergence. Thus young, gibbon-specific L1PAs may have gained this 20 bp deletion to evade ZNF93 while young great-ape specific elements gained the 129 bp deletion we observed previously . Gibbon L1PA elements with this 20-bp deletion may even hold a selective advantage over the more drastic 129 bp L1PA3 deletions, as elegant work from the Moran lab has recently shown that the 129 bp deletion in human L1PA3 elements alters L1PA splicing in a manner that can generate defective spliced integrated retrotransposed elements (SpIREs) . The smaller deletion found in gibbons may avoid generating these intermediates. Additionally, gibbon L1PA4 elements also experience a smaller coverage drop (typically near the ZNF765 binding site (Fig. 5b). Coverage drops in this area are found predominantly in L1PA4 instances with the ZNF93 binding site already deleted, indicating that this deletion (and the presumed escape from ZNF765 control) occurred after escape from ZNF93 control (Fig. 5c).
Novel orangutan-specific deletions are visible on the UCSC repeat browser
L1PA3 elements display an increased copy number in the orangutan genome, suggesting that these elements also had a lineage specific expansion, driven by escape from KZNFs or other restriction factors. Aligning of orangutan L1PA3 elements on the Repeat Browser L1PA3 consensus displayed a clear 11 bp deletion ~ 230 bp into the 5′ UTR that is not present in human, chimp or bonobo elements (Fig. 6a). However, analysis of existing KZNF ChIP-seq data, shows no specific factor that clearly correlates with this deletion. We may simply lack ChIP-seq data for the appropriate factor (including the possibility that the KZNF driving these changes evolved specifically within the orangutan lineage) to explain the evolutionary pattern seen in these orangutan-specific elements; alternatively, this mutation might alter some other aspect of L1PA fitness (e.g. splicing). Regardless, L1PA3 elements with this deletion were highly successful in spreading throughout the orangutan genome. Furthermore, L1PA3 instances with deletions in this region also harbor the 129 bp ZNF93-associated deletion, suggesting that this 11 bp deletion occurred after orangutan L1PA3 elements escaped ZNF93 control (Fig. 6b).
No major deletions are visible in primate L1PA2 elements
Mapping of L1PA2 elements in gorilla, bonobo, chimp and human to the Repeat Browser reveals only minor changes between these relatively young elements. (Fig. 6c) Although gorilla L1PA2 elements have greatly expanded compared to other primates, no significant gorilla-specific deletions are visible in our coverage plots; therefore the spread of gorilla elements may reflect the lack of a control factor that evolved in bonobo, chimpanzees and humans, or may reflect more subtle point mutations as we recently demonstrated for L1PA escape from ZNF649 control . Curiously, all four species show minor coverage drops in the area around nucleotide 250 (Fig. 6c), a region that roughly corresponds to the deletion event observed in orangutan L1PA3 elements. Although the deletion frequencies in primate L1PA2 are relatively low compared to the 11 bp L1PA3 orangutan deletion, this overall behavior is consistent with the model that this region is under adaptive selection - possibly to escape repression from a still unknown KZNF.
The UCSC Repeat Browser provides an interactive and accessible environment to study repeat biology and side-steps the problem of mistakenly mapping reads to an incorrect locus by generating consensus representations of every repeat class, and focusing on how genome-wide datasets interact with repeat sequences independent of their genomic locus. Here we use this consensus-based approach to identify deletion events in repeats across species that suggest a model by which L1PA escape occurs differently across the phylogenetic tree of old world monkeys (Fig. 7). We provide liftOver files so researchers can map their own genome-wide data to the consensus sequences and provide step-by-step instructions at https://repeatbrowser.ucsc.edu/tutorial/.
However, several caveats should be noted about Repeat Browser-based analyses. First, they rely entirely on RepeatMasker classifications (and in turn Dfam) and therefore depend on the quality of the annotations established in these collections. Second, the Repeat Browser often uses its own consensus sequences to display genomic data, with these choices biased by length in order to ensure proper visualization, which can otherwise be problematic in regions where sequence is inserted. However, custom versions of the browser allow users to provide a custom consensus sequence. Indeed, we previously used this approach to create consensuses of L1PA3 subclasses when tracing an evolutionary arms race between ZNF93 and L1PA3 elements . However, the Repeat Browser and other consensus-based approaches risk diluting important, biologically relevant signal driven by a few instances of a repeat type that may affect the cell by virtue of their genomic location instead of their sequence. In these cases, the majority of instances in these families may generate no signal and produce an underwhelming “composite” Repeat Browser signal whereas an individual genomic locus may produce a strong, reproducible, and functionally relevant signal. Therefore, we recommend that Repeat Browser analysis be used in combination with existing genomic approaches for repeat analysis [10, 38,39,40]. Finally, the existence of the UCSC Repeat Browser as a complete “repeat genome collection” available for download should allow manipulation and utilization of repeat consensus sequences with a large set of existing, standard genomics tools, thereby enhancing the investigation of repeat sequence biology. We expect that the repeat community will make creative use of this tool beyond the workflows suggested here.
The UCSC Repeat Browser provides a fully interactive environment, analogous to the UCSC Human Genome Browser, to study repeats. We show here that this environment provides an intuitive visualization tool for analysis and hypothesis-generation. For instance, here we use the Repeat Browser to demonstrate that sequence-specific deletions in repeats potentially driven by the activity of cellular repressors occurs independently in different species. The Repeat Browser is currently available at: https://repeatbrowser.ucsc.edu.
Availability and requirements
Project name: The UCSC Repeat Browser.
Project home page: https://repeatbrowser.ucsc.edu
Operating system(s): Standard Web Browser.
Programming language: Python, bash.
License: Freely available for academic, nonprofit, and personal use.
Any restrictions to use by non-academics: Use of liftOver requires commercial license: http://genome.ucsc.edu/license
Slotkin RK. The case for not masking away repetitive DNA. Mob DNA. 2018;9:15.
Chuong EB, Elde NC, Feschotte C. Regulatory activities of transposable elements: from conflicts to benefits. Nat Rev Genet. 2016. https://doi.org/10.1038/nrg.2016.139.
Pastuzyn ED, et al. The neuronal gene arc encodes a repurposed Retrotransposon gag protein that mediates intercellular RNA transfer. Cell. 2018;173:275.
Ding Y, Berrocal A, Morita T, Longden KD, Stern DL. Natural courtship song variation caused by an intronic retroelement in an ion channel gene. Nature. 2016;536:329–32.
Tam OH, et al. Postmortem Cortex Samples Identify Distinct Molecular Subtypes of ALS: Retrotransposon Activation, Oxidative Stress, and Activated Glia. Cell Rep. 2019;29:1164–1177.e5.
Bao W, Kojima KK, Kohany O. Repbase update, a database of repetitive elements in eukaryotic genomes. Mob DNA. 2015;6:11.
Smit, A., Hubley, R. & Green, P. RepeatMasker Open-4.0. http://www.repeatmasker.org.
Jurka J, Walichiewicz J, Milosavljevic A. Prototypic sequences for human repetitive DNA. J Mol Evol. 1992;35:286–91.
Hubley R, et al. The Dfam database of repetitive DNA families. Nucleic Acids Res. 2016;44:D81–9.
Goerner-Potvin, P. & Bourque, G. Computational tools to unmask transposable elements. Nat. Rev. Genet. 1 (2018). https://doi.org/10.1038/s41576-018-0050-x.
Jacobs FMJ, et al. An evolutionary arms race between KRAB zinc-finger genes ZNF91/93 and SVA/L1 retrotransposons. Nature. 2014;516:242–5.
Sun X, et al. Transcription factor profiling reveals molecular choreography and key regulators of human retrotransposon expression. Proc Natl Acad Sci U S A. 2018;115:E5526–35.
Ito J, et al. Systematic identification and characterization of regulatory elements derived from human endogenous retroviruses. PLoS Genet. 2017;13.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Kent WJ, et al. The human genome browser at UCSC. Genome Res. 2002;12:996–1006.
Haeussler M, et al. The UCSC genome browser database: 2019 update. Nucleic Acids Res. 2019;47:D853–8.
Jern P, Sperber GO, Ahlsén G, Blomberg J. Sequence variability, gene structure, and expression of full-length human endogenous retrovirus H. J Virol. 2005;79:6325–37.
Young NL, Bieniasz PD. Reconstitution of an infectious human endogenous retrovirus. PLoS Pathog. 2007;3:0119–30.
Kohany O, Gentles AJ, Hankus L, Jurka J. Annotation, submission and screening of repetitive elements in Repbase: RepbaseSubmitter and censor. BMC Bioinformatics. 2006;7:474.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27:573–80.
Hinrichs AS, The UCSC. Genome browser database: update 2006. Nucleic Acids Res. 2006;34:D590–8.
Kent WJ. BLAT---The BLAST-Like Alignment Tool. Genome Res. 2002;12:656–64.
Harrow J, et al. GENCODE: the reference human genome annotation for the ENCODE Project. Genome Res. 2012;22:1760–74.
ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(57–74).
Imbeault M, Helleboid P-Y, Trono D. KRAB zinc-finger proteins contribute to the evolution of gene regulatory networks. Nature. 2017;543:550–4.
Najafabadi HS, et al. C2H2 zinc finger proteins greatly expand the human regulatory lexicon. Nat Biotechnol. 2015;33:555–62.
Thomas JH, Schneider S. Coevolution of retroelements and tandem zinc finger genes. Genome Res. 2011;21:1800–12.
Leinonen R, Sugawara H. Shumway, M & International Nucleotide Sequence Database Collaboration The sequence read archive. Nucleic Acids Res. 2011;39:D19–21.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9.
Zhang Y, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137.
Khan H, Smit A, Boissinot S. Molecular evolution and tempo of amplification of human LINE-1 retrotransposons since the origin of primates. Genome Res. 2006;16:78–87.
Konkel MK, Walker JA, Batzer MA. LINEs and SINEs of primate evolution. Evol Anthropol. 2010;19:236–49.
Han K, et al. Mobile DNA in Old World monkeys: A glimpse through the rhesus macaque genome. Science (80-. ). 2007;316:238–40.
Schmitges FW, et al. Multiparameter functional diversity of human C2H2 zinc finger proteins. Genome Res. 2016;26:1742–52.
Larson PA, et al. Spliced integrated retrotransposed element (SpIRE) formation in the human genome. PLoS Biol. 2018;16.
Fernandes JD, et al. KRAB Zinc Finger Proteins coordinate across evolutionary time scales to battle retroelements. bioRxiv. 2018;429563. https://doi.org/10.1101/429563.
Jeong, H. H., Yalamanchili, H. K., Guo, C., Shulman, J. M. & Liu, Z. An ultra-fast and scalable quantification pipeline for transposable elements from next generation sequencing data. in Pacific Symposium on Biocomputing 0, 168–179 (World Scientific Publishing Co. Pte Ltd, 2018).
Jin Y, Tam OH, Paniagua E, Hammell M. TEtranscripts: a package for including transposable elements in differential expression analysis of RNA-seq datasets. Bioinformatics. 2015;31:3593–9.
Kong Y, et al. Transposable element expression in tumors is associated with immune infiltration and increased antigenicity. Nat Commun. 2019;10:5228.
We thank A. Smit, R. Hubley, and A. Ewing for helpful discussions about repeat consensus choice and annotations and Dfam. We thank J. Armstrong, F. Jacobs and D. Greenberg and all members of the Haussler lab for helpful comments and discussion.
This work was supported by EMBO ALTF 292–2011 and 4U41HG002371 to MH, F32GM125388 to JDF and 1R01HG010329 to SRS. DH is an investigator of the Howard Hughes Medical Institute.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Fernandes, J.D., Zamudio-Hurtado, A., Clawson, H. et al. The UCSC repeat browser allows discovery and visualization of evolutionary conflict across repeat families. Mobile DNA 11, 13 (2020). https://doi.org/10.1186/s13100-020-00208-w
- Krab zinc finger proteins