Fake IDs? Widespread misannotation of DNA transposons as a general transcription factor
Genome Biology volume 24, Article number: 260 (2023)
Accurate annotation of genes and transposable elements (TEs) is vital for understanding genomes, but current annotation pipelines often misannotate TEs as genes. This study reveals how the general transcription factor II-I repeat domain-containing protein 2 (GTF2IRD2) erroneously annotated DNA transposons in non-mammalian species, as it contains a 3′ fused hAT transposase domain. We also demonstrate the generality of this problem by identifying misannotated TEs as genes in other vertebrate genomes. Such misannotations can lead to errors in phylogenetic analyses and wasted time for investigators. The study proposes adding a final TE-check to gene annotation pipelines to mitigate this problem.
Annotation of transposable elements (TEs) is essential for understanding genome structure and function; however, misannotation may result in erroneous classification of TEs as trans-regulatory elements such as transcription factors. TEs can introduce genetic novelty to the host genome, and an example of this is the exaptation of TEs into the human general transcription factor II-I repeat domain-containing protein 2 (GTF2IRD2). GTF2IRD2 contains a Charlie8-like element positioned at the 3′ end (C terminus) of the gene model/full protein that has retained transposable element features such as catalytically active DDE/RW amino acids (D; aspartic acid, E; glutamic acid, R; arginine, and W; tryptophan) required for transposition . The Charlie transposon (DNA transposon; hAT superfamily) is an ancient autonomous group of transposons abundant in mammalian genomes, including humans . Charlie transposons are defined by their target site duplication (TSDs) and terminally inverted repeats (TIRs). At the same time, the protein-coding sequence of the transposase may vary between Charlie elements from different species . GTF2IRD2 is found in mammals and is predicted to be in many reptiles, amphibians, and bony fishes . However, upon closer inspection, some non-mammalian GTF2IRD2 sequences appear to be hAT transposons, not transcription factors. Here, we have used structural and phylogenetic analysis to resolve the widespread misannotation of non-mammalian DNA transposons as GTF2IRD2 transcription factors. In addition, we demonstrated that the issue of misannotation is indeed widespread by finding several instances of TEs of different classes incorrectly predicted as genes in a selection of genome assemblies. We believe that this paper addresses an important issue; it has implications not only for the study of TEs but is also relevant as similar misannotations could also cause the misinterpretation of other results that depend on reliable gene annotation.
Results and discussion
While annotating hAT-6 transposons in Testudines genomes, we noticed that GTF2IRD2/2A, a human general transcription factor, was the top BLASTN result when using hAT-6 transposons as a query to search non-mammalian genomes. This was unexpected as hAT-6 has hallmarks of a functional transposase, such as TIRs and TSDs, and has the functional motifs required for transposition, such as the DDE/RW residues in the translated open reading frame (ORF) [4, 5]. NCBI’s eukaryotic gene annotation pipeline prefers to use experimental evidence when annotating genes but uses an ab initio model to predict optimal coding sequence (CDS) alignments when there is no experimental data (https://www.ncbi.nlm.nih.gov/genome/annotation_euk/process/). Specifically, Gnomon is meant to exclude gene predictions with high homology to transposable or retro-transposable elements from the final gene models; however, the eukaryotic annotation pipeline appears to lack a final TE filtering step after integrating RefSeq annotations. This may explain how genes such as GTF2IRD2/2A, which contains an integrated Charlie8-like element (a hAT-like transposase), can lead to hAT-6 being predicted to be a transcription factor in non-mammalian genomes . When searching Interpro for GTF2IRD2/2A sequences, we saw that 485 proteins were annotated as GTF2IRD2/2A proteins, but only 4 of these have been reviewed in human, cow and mouse genomes (https://www.ebi.ac.uk/interpro/entry/InterPro/IPR042224/protein/reviewed/#table, Accessed 11th November 2022). To determine if additional hAT transposons were incorrectly annotated as GTF2IRD2/2A, we examined the phylogeny of a set of protein sequences annotated as GTF2IRD2/2A or GTF2IRD2/2A-like from mammals, reptiles, bony fishes, and amphibians (Fig. 1A). Birds were excluded from analysis as the sequences annotated as GTF2IRD2/2A had no significant similarity to either mammalian GTF2IRD2/2A or to any DNA TE structures. Any sequence homology in birds was limited to the N terminus of the mammalian GTF2IRD2/2A, an indication of potential similarity to the ancestral GTF2IRD2/2A protein prior to the exaptation of Charlie8 (Additional file 1: Table S1).
There was a clear distinction between mammalian and non-mammalian GTF2IRD2/2A gene models. Mammalian GTF2IRD2/2A are correctly annotated as transcription factors as they have an N terminus ~ 400 aa long containing a GTF2I-Like repeat domain, a zinc finger binding domain, and an integrated Charlie8-like element at the C terminus . However, this is not the case with non-mammalian sequences. Multiple alignment of GTF2IRD2/2A from mammals and non-mammals with hAT-6 transposons exclusively shows high similarity alignment at the 3′ end (Fig. 1B, C). This is consistent with the position of the Charlie-like element in GTF2IRD2/2A, demonstrating that non-mammalian GTF2IRD2/2A are not likely to be transcription factors (TFs), but TEs.
Non-mammalian GTF2IRD2/2A from fish, reptiles, and frogs were manually curated to identify signature motifs of DNA transposons as they appear more closely related to hAT-6 transposons than mammalian GTF2IRD2/2A. Twenty-eight sequences were classified as autonomous hAT transposons from several different species, particularly from the genome of the Atlantic salmon (Salmo salar) (Additional file 1: Table S2). TSDs and TIRs are characteristic of hAT transposons and were found in most sequences, while sequences without them were classified as partial, non-autonomous transposons (Additional file 1: Table S3). Furthermore, DDE/RW residues required for transposition were identified for most of the newly curated hAT transposons (Additional file 2). Finally, the 5′ and 3′ TIRs were mostly conserved across the misannotated hAT transposons, which further demonstrates the degree of misannotation (Additional file 1: Table S4).
While we have demonstrated the widespread effects of misannotation in a single case of hAT-6 transposons, we wanted to determine if this was an isolated case or whether the problem extends to other TE classes and predicted gene models. To examine the potential magnitude of TE misannotation, we assessed a dataset of mammalian and non-mammalian vertebrate genome assemblies for instances of TE annotations intersecting with gene annotations from the UCSC Genome Browser . We focused specifically on instances where 100% of a TE sequence overlapped with the protein-coding region (CDS) of gene and gene predictions.
As shown in Fig. 2A, we found hundreds of instances of TEs overlapping completely with the CDS annotations of non-mammalian vertebrate genomes. Notably, overlapping events were very abundant in zebrafish and the western clawed frog genomes, both of which are used widely in a variety of biochemical and evolutionary analyses. In both the zebrafish and western clawed frog, a diverse range of TE classes were found overlapping with genes (Fig. 2B, C). The predominant class of TE were hAT DNA transposons and Ty3/Gypsy LTR retrotransposons. This observation is consistent with both the distribution and high abundance of hAT and Ty3/Gypsy elements across vertebrates [2, 7].
In contrast, in mammalian genomes, all overlapping events were correctly annotated as ‘transposase-like’ in the corresponding gene records in NCBI. This was expected as there have been significantly more concerted efforts to curate mammalian genomes with respect to both genes and repeats. The same was true for most avian genomes investigated for misannotation as these have benefited from developing interest in high-quality annotation . This level of annotation detail was notably absent in the other vertebrate genomes studied as instances where TEs overlapped with genes were not noted in the gene records.
While the number of overlapping events in vertebrates highlights how widespread misannotation can be, we further confirmed this by manually curating a sample of genes for all vertebrate genomes investigated for TE(s). Some genes predicted to be proteins or left uncharacterised are actually TEs (Additional file 1: Table S5). DNA transposable elements were often misannotated as zinc finger binding proteins. This is not surprising, as autonomous DNA TE ORFs encode DNA binding domains which can explain why some TEs were mistakenly predicted as non-TE genes. Uncharacterised loci were the predominant type of misannotation for a variety of both class 1 and 2 TEs. This finding shows a deficiency in genome annotation pipelines in combining RepeatMasker results to properly characterise unknown genes, which with minimal manual curation can be correctly annotated as TE sequences. Figure 2D illustrates a selection of TEs overlapping with gene annotations. In most cases, the TE spans the majority of the CDS, further showing that the misannotation of TEs is not limited to a particular class of TE or type of gene prediction.
We focused specifically on the CDS of gene predictions and instances where 100% of the TE overlaps with those predictions. Even with stringent parameters and a limited number of organisms, we have demonstrated that misannotation occurs frequently and we strongly encourage better curation of repeats as it affects gene annotation for a more comprehensive picture of genomes, especially of emerging fish and reptile genomes. A simple additional screening for sequence homology of gene models with TEs as the last step of gene model annotation would mitigate/eliminate this problem.
In this study, we present a case of widespread incorrect annotation of hAT DNA transposons as GTF2IRD2/2A. This has led to at least 28 new instances of hATs from reptiles, bony fishes, and amphibians that had been overlooked and could have been incorrectly used as transcription factors in other analyses. We have also demonstrated that misannotation occurs frequently for different types of both class 1 and class 2 TEs, especially across vertebrate genomes. Correct annotation is a vital step in furthering our understanding of genome evolution, and misannotation of TEs as trans-regulatory genes such as TFs affects downstream research and can confound phylogenetic analysis.
Manual curation of hAT transposons from GTF2IRD2/2A sequences
A set of GTF2IRD2/2A and GTF2IRD2/2A-like protein sequences were downloaded from NCBI. The search was limited to species belonging to the Actinopterygii, Reptilia, and Amphibia classes. To determine whether these GTF2IRD2/2A proteins were actually DNA transposons, extensive manual curation was performed to locate characteristic sequence features such as TIRs, TSDs, and ORFs. To identify hAT transposons, hAT-6_TSE (manuscript in prep) was used as a query in a BLASTP 2.7.1 + search against a set of mammal, reptile, amphibian, bony fish, and bird genomes containing GTF2IRD2/2A gene annotations . The corresponding nucleotide sequence of each top hit was extended 1000 bp in flanking regions where possible and used for manual annotation of TIRs and TSDs characteristic of hAT transposons. ORFs were searched using GENSCAN (http://hollywood.mit.edu/GENSCAN.html) and searched for DD/E and RW residues. Sequences that contained 5′ and 3′ TSDs, TIRs, and an intact ORF were classified as autonomous hAT transposons [9, 10] (Additional file 1: Table S6). The best match for each new autonomous hAT was found using Repbase (https://www.girinst.org/repbase/; accessed August 2020), and both 5′ and 3′ TIRs were aligned using MAFFT v7.310 to view conserved nucleotides [11, 12]. hATs misannotated as GTF2IRD2-like were aligned using MAFFT to mammalian GTF2IRD2/2A to confirm they had homology to GTF2IRD2/2A’s Charlie8-like domain (Additional file 3).
GTF2IRD2/2A and GTF2IRD2/2A-like protein sequences from mammals, reptiles, amphibians, and bony fishes were aligned to the Hermes and hAT-6 transposons using MAFFT v7.319 . The alignment was trimmed using CLipKit v1.3.0, and IQTree v1.6.12 was used for tree reconstruction with JTT + F + I + G4 as the best-fit model with 20 maximum likelihood trees and 1000 bootstraps [14,15,16,17].
Determining whether the misannotation of TEs is a recurring event or limited to hAT-6 s
To examine the breadth of TE misannotation in other organisms, gene and gene prediction files were downloaded from the UCSC Genome Browser (https://genome.ucsc.edu/index.html) and the corresponding RepeatMasker files for mammals and vertebrates. Using BEDTools v2.30.0, the RepeatMasker output was overlaid with the gene and gene prediction files to extract regions of overlap between TEs and genes . Stringent parameters were used to exclude only simple repeats and include TEs that overlapped 100% with protein-coding regions (CDS), which were selected for analysis (bedtools intersect -a repeats.bed -b cds.bed -wo -f 1.00 > 100_overlap_cds.bed). A sample of genes that overlapped with TEs greater than 1000 base pairs were selected to check whether any true misannotation event occurred. The best TE match for the sample of genes was found through Repbase . Assemblies used for analysis, frequency (TE:Protein-coding gene count), and coordinates of misannotation events are in Additional file 1: Table S7.
Availability of data and materials
The dataset(s) supporting the conclusions of this article are included within the article (and its additional file(s)).
General transcription factor II-I repeat domain-containing protein 2
General transcription factor II-I repeat domain-containing protein 2A
Tipney HJ, Hinsley TA, Brass A, Metcalfe K, Donnai D, Tassabehji M. Isolation and characterisation of GTF2IRD2, a novel fusion gene and member of the TFII-I family of transcription factors, deleted in Williams–Beuren syndrome. Eur J Human Genet. 2004:551–60. https://doi.org/10.1038/sj.ejhg.5201174
Arensburger P, Hice RH, Zhou L, Smith RC, Tom AC, Wright JA, et al. Phylogenetic and functional characterization of the hAT transposon superfamily. Genetics. 2011;188:45–57.
Jurka J. Repeats in genomic DNA: mining and meaning. Curr Opin Struct Biol. 1998:333–7. https://doi.org/10.1016/s0959-440x(98)80067-5
Blum M, Chang H-Y, Chuguransky S, Grego T, Kandasaamy S, Mitchell A, et al. The InterPro protein families and domains database: 20 years on. Nucleic Acids Res. 2021;49:D344–54.
Nesmelova IV, Hackett PB. DDE transposases: structural similarity and diversity. Adv Drug Deliv Rev. 2010;62:1187–95.
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The human genome browser at UCSC. Genome Res. 2002;12:996–1006.
Goodier JL. Restricting retrotransposons: a review. Mob DNA. 2016;7:16.
Eöry L, Gilbert MTP, Li C, Li B, Archibald A, Aken BL, et al. Avianbase: a community resource for bird genomics. Genome Biol. 2015;16:21.
Zhang Z, Schwartz S, Wagner L, Miller W. A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000;7:203–14.
Burge C, Karlin S. Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997:78–94. https://doi.org/10.1006/jmbi.1997.0951
Bao W, Kojima KK, Kohany O. Repbase update, a database of repetitive elements in eukaryotic genomes. Mob DNA. 2015;6:11.
Kuraku S, Zmasek CM, Nishimura O, Katoh K. aLeaves facilitates on-demand exploration of metazoan gene family trees on MAFFT sequence alignment server with enhanced interactivity. Nucleic Acids Res. 2013;41:W22–8.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Steenwyk JL, Buida TJ 3rd, Li Y, Shen X-X, Rokas A. ClipKIT: a multiple sequence alignment trimming software for accurate phylogenomic inference. PLoS Biol. 2020;18:e3001007.
Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37:1530–4.
Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35:518–22.
Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.
We would like to thank Terry Bertozzi and Christine Elsik for their expertise and helpful comments. We would also like to thank our lab members for their ongoing support throughout this study and beyond.
The review history is available as Additional file 5.
Peer review information
Tim Sands was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
This research was funded by the University of Adelaide.
Ethics approval and consent to participate
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Pairwise Identity Matrix of predicted GTF2IRD2/2A sequences in birds vs GTF2IRD2/2A in mammals. Table S2. Species from the Actinopterygii, Reptilia and Amphibia class containing hAT transposons that were incorrectly annotated as GTF2IRD2/2A. Table S3. 3′ and 5′ TIRs of TEs misannotated as GTF2IRD2/2A. Table S4. Multiple alignment of the 5′ and 3′ TIRs of Actinopterygii, Amphibia and Reptilia hAT transposons derived from sequences misannotated as GTF2IRD2/2A. Table S5. A sample of TEs from vertebrates which were incorrectly annotated as various other proteins and uncharacterised loci. Table S6. Nucleotide sequences TEs misannotated as GTF2IRD2/2A. Table S7. Sample of misannotation events. Table S8. Conserved domains of GTF2IRD2/2A compared to TEs
DDE/RW positions in new hAT TEs.
Text file of multiple alignment of mammalian GTF2IRD2/2A to predicted GTF2IRD2/2A sequences in vertebrates
Peer review history.
About this article
Cite this article
Hassan, N.T., Adelson, D.L. Fake IDs? Widespread misannotation of DNA transposons as a general transcription factor. Genome Biol 24, 260 (2023). https://doi.org/10.1186/s13059-023-03102-9