Skip to main content

Multiple horizontal transfer events of a DNA transposon into turtles, fishes, and a frog

Abstract

Horizontal transfer of transposable elements (HTT) has been reported across many species and the impact of such events on genome structure and function has been well described. However, few studies have focused on reptilian genomes, especially HTT events in Testudines (turtles). Here, as a consequence of investigating the repetitive content of Malaclemys terrapin terrapin (Diamondback turtle) we found a high similarity DNA transposon, annotated in RepBase as hAT-6_XT, shared between other turtle species, ray-finned fishes, and a frog. hAT-6_XT was notably absent in reptilian taxa closely related to turtles, such as crocodiles and birds. Successful invasion of DNA transposons into new genomes requires the conservation of specific residues in the encoded transposase, and through structural analysis, these residues were identified indicating some retention of functional transposition activity. We document six recent independent HTT events of a DNA transposon in turtles, which are known to have a low genomic evolutionary rate and ancient repeats.

Summary of species

Malaclemys terrapin terrapin (Diamondback turtle).

Malaclemys terrapin pileata (Mississippi diamondback terrapin turtle).

Trachemys scripta elegans (Red-eared slider turtle).

Chrysemys picta bellii (Western painted turtle).

Dermatemys mawii (Hickatee turtle).

Sternotherus odoratus (Common musk turtle).

Mesoclemmys tuberculata (Tuberculate Toad-headed turtle).

Etheostoma spectabile (Orangethroat darter fish).

Thalassophryne amazonica (Prehistoric monster fish).

Scophthalmus maximus (Turbot fish).

Syngnathus acus (Greater pipefish).

Scleropages formosus (Asian Arowana fish).

Xenopus tropicalis (Western clawed frog).

Introduction

A large proportion of eukaryotic genomes is composed of mobile repetitive sequences known as transposable elements (TEs). TEs are divided into two classes based on their mode of mobilisation. Retrotransposons (Class I) copy themselves using an RNA intermediate which is reverse-transcribed into cDNA and is integrated back into the host genome, whereas DNA transposons (Class II) move via a cut-and-paste mechanism mediated by an encoded transposase [1, 2]. One of the largest and most widespread superfamilies of DNA transposons are hAT elements. While there is sequence variation between families of hAT transposons, they are defined by ∽ 8 bp terminal target site duplication (TSDs) and ∽ 15 bp terminal inverted repeats (TIRs) [3]. The hAT transposase is a multidomain protein and the only hAT transposase crystal structure is from the Hermes transposase of the house fly (Musca domestica), which showed the presence of conserved residues and domains required for TIR recognition and transposition [4, 5].

Testudines (turtles) are a group of reptiles found in diverse ecological settings ranging from terrestrial, marine, and freshwater environments and are a sister group to Archosauria (birds and crocodilians) that diverged during the Permian-Triassic period approximately 257.4 Mya [6, 7]. The repetitive content of turtle genomes has not been studied extensively, but current studies indicate that TEs make up approximately 30% of the genome in various turtles and are dominated by LINEs and DNA transposons [8,9,10]. Turtle genomes generally do not show significant variation in size and appear to evolve slowly compared to other reptiles which makes analysis of repeats an area of interest [11]. Until now, horizontal transfer of DNA transposable elements (HTT) has not been documented in or between turtles, fishes, and a frog.

Horizontal transfer is the process by which genetic material is obtained from non-parental genomes/sources, as opposed to vertical transfer which is from parent to offspring [12, 13]. TEs, in particular, are widely spread through horizontal transfer in eukaryotes and can persist within the invaded genome [14]. HTT has been documented in several species, for example, SPIN (DNA TEs) and BovB elements (non-LTR retrotransposons) have colonised many squamate reptiles [15,16,17]. However, in both studies, evidence for HTT was notably absent in turtles. In Zhang et al. (2020), HTT of both retrotransposons and DNA transposons between turtles, fishes, and lizards were inferred, however, they did not specifically detect hAT-6_XT (first curated in Xenopus tropicalis by Kapitonov & Jurka, 2006) HTT events between species, HTT events between turtles or carry out an evolutionary analysis [18, 19]. In this study, we outline a rare horizontal transfer of a hAT-6_XT DNA TE between/into six turtles, fishes, and a frog.

Results and discussion

HTT amongst species of turtles, fishes, and a frog

Horizontal transposon transfer has been widely reported in the species discussed in this study, with new reports emerging describing HTT of both DNA and RNA transposable elements in ray-finned fish, amphibians, and reptiles [18, 20]. HTT between turtles, especially of DNA transposons, has been proposed to be a rare event based on our current understanding of turtle genome evolution [18, 20, 21]. We screened the genomes of approximately 100 species ranging from fishes, reptiles, mammals, birds, and insects to find hAT-6_XTs that share striking homology (Additional file 1; Table S1). We have identified the first case of horizontal transfer of a hAT-6_XT DNA transposon, first annotated in Xenopus tropicals, into turtles, ray-finned fishes, and a frog.

We show that hAT-6_XT TEs are distinct from other hAT and DNA TEs (Fig. 1). In addition, the high sequence similarity of these TEs from distant species, absence in closely related species, and discordant topologies of the species and hAT-6_XT phylogenies (Fig. 2; Additional file 2: Figure S1) make a compelling case for the horizontal transfer of this DNA TE. The first detected transfer of hAT-6_XT was between Xenopus tropicalis (Western clawed frog) and M.t. terrapin as it satisfied the previous criteria for HTT. A further 12 hAT-6_XTs were found across turtles, a frog, and fish. In addition, given the high sequence identity of hAT-6_XT between pairs of species (Fig. 2), hAT-6_XT may have repeatedly transferred between aquatic animals over time. While we could not determine the direction of transfer, we speculate that transfer is mediated by a parasitic donor.

Fig. 1
figure 1

Phylogeny of horizontally transferred hAT-6_XT TEs with a sample of other hAT TEs from Repbase. Tree constructed using IQTree 2 (1000 bootstraps) based on MAFFT protein alignments trimmed using Clipkit and formatted in iTOL. Support values under 65 are displayed at nodes. Turtle species are in blue, frogs are in pink, and fishes are in green

Fig. 2
figure 2

Sequence homology and phylogeny of hAT-6_XT relative to host species phylogeny. (A) Host species phylogeny (topology only) of hAT-6_XT host species. The time of divergence between some host species is shown in MYA. (B) Observed hAT-6_XT topology constructed using IQTree (1000 bootstraps) based on MAFFT protein alignments trimmed using Clipkit and plotted in iTOL. (C) BLASTN pairwise identity matrix of hAT-6_XT representative sequence alignments. Darker shading indicates a higher pairwise identity between two hAT-6_XTs. Turtle species are in blue, frogs are in pink, and fishes are in green

The structure of hAT-6_XT transposases indicates activity

To support our case for HTT of hAT-6_XT, we determined the predicted structure and the presence of functional domains required for mobility in the corresponding encoded transposases. The conservation of the TIRs/TSDs in full-length hAT-6_XT transposons is expected and required for HTT. In contrast, the absence of those features in fragmented hAT-6_XT transposons indicates degradation in the genome. We find conservation of such features in most examined hAT-6_XTs (Additional file 1; Table S2). In addition, the functional annotation of hAT-6_XT transposases shows the presence of the essential DDD/E catalytic triad required for transposition and for HTT (Additional file 2: Figure S2-3) [4]. These findings suggest that the hAT-6_XTs have the necessary sequence features and motifs to be active now and thus increase the likelihood that they were recently horizontally transferred.

hAT-6_XT expansion and divergence

As all hAT-6_XT from turtle genomes are remarkably similar, we selected two closely related species, M.t. terrapin and T.s. elegans, for downstream investigation. To understand the evolution of each hAT-6_XT transfer in the host species, we investigated the coverage of each element. The coverage and divergence plots of hAT-6_XT in turtles (Additonal file 2: Figure S4) show very low divergence which may be a product of the slow genome evolution [18] of turtles rather than recent HTT events [18]. When combined with patchy phylogenetic distribution in turtles/other species, and how unexpectedly similar hAT-6_XTs are across all species, the evidence points towards HTT in turtles. In comparison to the other species examined in this study, turtles also contained more full-length copies of hAT-6_XT. We were able to rule out horizontal transfer of hAT-6_XT into the most recent common ancestor, prior to the divergence of Trachemys scripta elegans (red-eared slider) and M.t. terrapin ∽ 16 MYA, and its subsequent preservation in the genomes of T.s. elegans ​and M.t. terrapin​. To determine if HTT of hAT-6_XT occurred independently into turtle species we used a presence/absence analysis to determine if insertions in each species were present or absent from homologous genomic regions (Additional file 2: Figure S5-S10). This analysis showed independent HTT for M.t. terrapin, ​T.s. elegans, C.p. bellii, D. mawii, and S. odoratus, but not for M.t. pilaeta (subspecies of M.t. terrapin). We did not observe full-length hAT-6_XT in the C.p. bellii genome assembly, and this could be the result of degradation of hAT-6_XT following HTT, or incompleteness of the assembly (Additional file 2: Figure S28). M.t. terrapin and M.t. pileata share eight hAT-6_XT insertions indicating HTT into the common ancestor of M.t. terrapin with a subsequent sub-species expansion of hAT-6_XT in M.t. pileata (Additional file 2: Figure S9).

In E. spectabile, where hAT-6_XT divergence is much greater, we also observed a highly amplified region of ∽ 800 bp from hAT-6_XT (Additional file 2: Figure S4). Upon further investigation, we found that this region is a hAT-6_XT derived non-autonomous DNA transposon (Additional file 2: Figure S11). By looking at the Kimura-based divergence of both the hAT-6_XT and hAT-6_XT derived non-autonomous DNA transposon from E. spectabile (hAT-6N1_XT_ESp), we see that copy number increase of hAT-6N1_XT_ESp occurred at the same time as copy number increase of hAT-6_XT (Additional file 2: Figure S12). We also observed two additional instances of hAT-6_XT copy number increase in E. spectabile, however, these two peaks were also observed in all instances where hAT-6_XT was present, which could indicate reactivation of hAT-6_XT rather than recent HTT (Additional file 2: Figure S13-S25). However, taken together with other evidence of HTT, we hypothesise that there may have been repeated HTT into E. spectabile, as opposed to a scenario where the more diverged copies identified in the alignments are degraded versions of a similar hAT.

In FigureS26 (Additional File 2), we show the frequency of hAT-6_XT insertion and divergence relative to other vertically inherited DNA TEs and retrotransposons in turtles. In both M.t. terrapin (Fig. 6.1) and T.s. elegans (Fig. 6.2), where there has been independent HTT, hAT-6_XT is present in low quantities compared to other elements but shows less than 1% divergence, further supporting the case for HTT. When comparing hAT-6_XT to hAT-3_MTT/TSE, a DNA TE shared by both M.t. terrapin and T.s. elegans, we see that the hAT-6_XT Kimura divergence and abundance are more similar than expected between the two turtles, despite the divergence between the two species estimated to be 14.5–15.6 MYA [22]. In addition, hAT-3s from both turtles are nearly identical in sequence but show a higher Kimura divergence over time which is consistent with an ancestral repeat. We also show that DNA TEs and retrotransposable elements have similar patterns of copy number increase and divergence, likely indicating the presence of population bottlenecks leading to fixation of insertions [23]. It is important to note that other vertically inherited TEs show peaks consistent with reactivation over time in the turtle genomes, which is not observed with hAT-6_XT even though hAT-6_XTs in turtles appear to have two peaks (marked with arrows) (Additional file 2: Figure S26). When investigating the sequences corresponding to the second peak between 34 and 52% divergence for hAT-6_XT, manual curation revealed that the alignments showed no resemblance to hAT-6_XT, but did resemble other hATs. We observe a similar pattern in hAT-like elements in snakes suggesting they are remnants of other degraded and ancient hATs (Additional file 2; Figure S27).

Aquatic environments and parasite-host relationships may facilitate HTT

The high percentage identity and low divergence of hAT-6_XT between species are likely the result of independent transfers from an unknown donor(s). The species in this study are largely found in semi-aquatic or aquatic environments but are mostly geographically distant. This indicates possible donor(s) that is/are likely ubiquitously aquatic and/or parasitic in nature. This finding also supports the frequent and recent transfer of transposons in aquatic environments [18, 20, 24]. Horizontal transfer is known to be facilitated by parasite-host relationships with previous findings showing HTT of a retrotransposon from parasitic nematodes [25, 26].

The presence of parasites in aquatic ecosystems involves complex life cycles and may contribute to HTT. The larval stages of certain parasites go through freshwater fish as intermediate hosts and reach turtles as the final hosts, therefore a similar mechanism may happen for frogs and tadpoles [27]. Protozoans and some Metazoans, particularly leeches, exhibit generalist feeding behaviour and can transmit pathogens like trypanosomes [28, 29]. Host specificity of parasites is variable, with some parasites reaching dead-end hosts due to dietary or environmental factors [29, 30]. Parasites in aquatic systems adeptly utilise paratenic hosts, where development is paused until the intermediate host is consumed, thus advancing the parasites through the food chain [29,30,31]. In addition, there is a geographical overlap between the turtles M.t. terrapin, T.s. elegans, C.p. belli, and S. odoratus and E. spectabile (Orangethroat darter) suggesting hAT-6_XT HTT into these turtles may have occurred from E. spectabile into turtles through a shared parasitic vector, such as darter fish parasites, as all five species have overlapping geographical distributions in North America [32]. However, we cannot rule out the possibility of HTT between turtle species. Finally, as there is no geographical overlap between X. tropicalis and the other fishes or turtles, more genomic data from geographically/ecologically overlapping species are required to determine possible donors (Additional file 1; Table S3). As a whole, our study documenting the horizontal transfer of hAT-6_XT DNA transposons among turtles, fishes, and a frog sheds light on the interplay of genetic elements across diverse aquatic species, and provides insights into their genome evolution.

Conclusions

Overall our study expands our knowledge of HTT in aquatic species and especially the evolution of HTT repeats in the slow-evolving genomes of turtles. We have documented new, recent horizontal transfer events between/into turtles, ray-finned fishes, and a frog, showing that HTT may be more common than expected in turtles. Our findings support the notion that HTT is a common occurrence in ray-finned fishes and suggest that aquatic environments may facilitate a large number of HTT events. The direction of HTT of hAT-6_XT into M.t. terrapin, T.s. elegans, C.p. belli, S. odoratus turtles is possibly from E. spectabile as the species share habitat and overlap geographically, but we cannot rule out transfer from an unknown donor into all five species. While the donors for HTT are unknown, our results and others from the literature suggest the existence of a cryptic aquatic network of horizontal transfer that is widely distributed given the geographic distances between the HTT recipients.

Methods

Identification and classification of horizontally transferred DNA transposon candidates

We performed an ab initio repeat annotation for the genome of ​M.t. terrapin using the Comprehensive ab initio Repeat Pipeline (CARP) [33]. We identified one DNA transposon in M.t. terrapin which was originally curated in X. tropicalis (Western clawed frog) as hAT-6_XT (Additional file 1; Table S2). We thus renamed the DNA TE in M.t. terrapin as hAT-6_XT_MTT. To find more similar sequences to hAT-6_XT_MTT, we performed an extensive local alignment using hAT-6_XT_MTT as a query against the following taxonomic groups: crocodilians, birds, frogs and toads, snakes, turtles, and fish using sensitive BLASTN 2.7.1 + parameters [34] (​word-size: 7, match/mismatch score: 4, -5​) against RefSeq Representative Genomes [35] (Additional file 1; Table S1). A cutoff of 75% identity and 90% coverage to hAT-6_XT_MTT was used to find full-length transposon sequences (​-blastn -e-value 1e-10​). We did global alignments of each sequence back to the query species using MAFFT v7.450 using NCBI coordinates to detect TIRs/TSDs [35, 36]. We classified sequences as full-length transposons or fragments based on the presence of TIRs and TSDs. Representative sequences were selected based on the previous characteristics when multiple copies were found in the respective genomes. Through this process, 13 sequences similar to hAT-6_XT_MTT were found. Open reading frames (ORFs) were identified using GENSCAN and a pairwise identity matrix (PIM) for all 13 sequences was made using BLASTN [37].

Construction of repeat phylogenies

We aligned the above 13 horizontally transferred hAT-6_XT repeats and representative hAT sequences from RepBase using MAFFT v7.450 [35, 36, 38] (​FFT-NS-1 model​). To select conserved regions for phylogenetic analyses we processed the multiple alignments using Clipkit [39], allowing small final blocks, gap positions within the final blocks and relaxed flanking positions (smart-gap). Two independent tree-building tools were used: IQTree 2 (JTT + G4; 1000 bootstraps​) and Fasttree 2.1 (​JTT + CAT model​) [40, 41]. All tree files were visualized and edited using iTOL [42].

Construction of species phylogeny

We used TimeTree to construct a species tree of amphibians, bony fishes, and reptiles (https://timetree.org/). The tree was visualised on the Interactive Tree of Life (iTOL) [42, 43]. In the case where a species of interest was not available on TimeTree, we substituted a species from the same clade as a proxy.

Protein structural analysis of hAT-6_XTs

We determined the structural features of three of the hAT-6_XT DNA transposons we discovered with a protein sequence exceeding 600 amino acids - hAT-6_XT_TSE (T.s. elegans), hAT-6_XT_SFo (S. formosus), and hAT-6_XT_SAc (S. acus) - using AlphaFold [40]. AlphaFold output was visualized using PYMOL (The PyMOL Molecular Graphics System, Version 2.0 Schrödinger, LLC). DDD/E and RW amino acid residues in transposase sequences were identified using M-COFFEE with alignment to the Hermes transposase (Additional file 2; Figure S3) [44].

Divergence and genome coverage of horizontally transferred DNA transposons

We identified orthologous full-length hAT-6_XT sequences from the genomes we searched (see above) using reciprocal BLAST searches [20, 45]. A custom database was made using the relevant turtle, fishes, and frog genomes as stated above to find the reciprocal best hit using hAT-6_XT_MTT as a nucleotide query from the genome of interest. TE-Aid (https://github.com/clemgoub/TE-Aid) was used to align a sample of full-length hAT-6_XT DNA transposon sequences to each other from the genome it was curated from to visualise the frequency of the complete sequence and fragments.

Kimura distances

We calculated the Kimura 2-parameter distance for all hAT-6_XT TEs using RepeatMasker version 4.1.5 and the script calcDivergenceFromAlign.pl to obtain the relative age distribution of the TEs in the genome (http://www.repeatmasker.org/RMDownload.html). We also calculated the Kimura 2-parameter distance for other repeats in the Testudine genomes (T.s. elegans and M.t. terrapin) to determine the age of hAT-6_XT compared to other repeats. We excluded short simple repeat alignments from the analysis.

Presence/absence test

To determine whether HTT of hAT-6_XT in all the turtle species resulted from a single ancestral event or independent transfers, we conducted a presence/absence test. Using BLASTN 2.7.1 + parameters [34] (blastn-short), all hAT-6_XT intervals were retrieved from the turtle genomes using the representative hAT-6_XT from that genome as a query. The blast output was cleaned to remove hits < 100 bp and converted to bed format with blast2bed.sh (https://github.com/nterhoeven/blast2bed). A genome index file was generated and the cleaned intervals were extended by 1500 bp on each side using BEDTools and SAMtools (slopBed, faidx) [46, 47]. Fasta formatted sequences were extracted for the extended bed intervals using BEDTools (getfasta). The extended hAT-6_XT sequences from M.t. terrapin were aligned to each turtle with Gepard and displayed as a dot plot [47, 48].

Data availability

The dataset(s) supporting the conclusions of this article are included within the article (and its additional.

file(s)).

References

  1. Jurka J, Kapitonov VV, Kohany O, Jurka MV. Repetitive sequences in complex genomes: structure and evolution. Annu Rev Genomics Hum Genet. 2007;8:241–59.

    Article  CAS  PubMed  Google Scholar 

  2. Feschotte C, Pritham EJ. DNA transposons and the evolution of eukaryotic genomes. Annu Rev Genet. 2007;41:331–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Atkinson PW. hAT Transposable Elements. Microbiol Spectr [Internet]. 2015;3. https://doi.org/10.1128/microbiolspec.MDNA3-0054-2014.

  4. Hickman AB, Ewis HE, Li X, Knapp JA, Laver T, Doss A-L, et al. Structural basis of hAT transposon end recognition by Hermes, an octameric DNA transposase from Musca domestica. Cell. 2014;158:353–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Hickman AB, Voth AR, Ewis H, Li X, Craig NL, Dyda F. Structural insights into the mechanism of double strand break formation by Hermes, a hAT family eukaryotic DNA transposase. Nucleic Acids Res. 2018;46:10286–301.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Chiari Y, Cahais V, Galtier N, Delsuc F. Phylogenomic analyses support the position of turtles as the sister group of birds and crocodiles (Archosauria). BMC Biol. 2012;10:65.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Wang Z, Pascual-Anaya J, Zadissa A, Li W, Niimura Y, Huang Z, et al. The draft genomes of soft-shell turtle and green sea turtle yield insights into the development and evolution of the turtle-specific body plan. Nat Genet. 2013;45:701–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Boissinot S, Bourgeois Y, Manthey JD, Ruggiero RP. The mobilome of reptiles: evolution, structure, and function. Cytogenet Genome Res. 2019;157:21–33.

    Article  PubMed  Google Scholar 

  9. Chalopin D, Naville M, Plard F, Galiana D, Volff J-N. Comparative analysis of transposable elements highlights mobilome diversity and evolution in vertebrates. Genome Biol Evol. 2015;7:567–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Tollis M, DeNardo DF, Cornelius JA, Dolby GA, Edwards T, Henen BT et al. The Agassiz’s desert tortoise genome provides a resource for the conservation of a threatened species [Internet]. PLOS ONE. 2017. p. e0177708. https://doi.org/10.1371/journal.pone.0177708.

  11. Shaffer HB, Bradley Shaffer H, Minx P, Warren DE, Shedlock AM, Thomson RC et al. The western painted turtle genome, a model for the evolution of extreme physiological adaptations in a slowly evolving lineage [Internet]. Genome Biology. 2013. p. R28. https://doi.org/10.1186/gb-2013-14-3-r28.

  12. Burmeister AR. Horizontal gene transfer. Evol Med Public Health. 2015;2015:193–4.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Danchin EGJ. Lateral gene transfer in eukaryotes: tip of the iceberg or of the ice cube? BMC Biol. 2016. p. 101.

  14. Peccoud J, Cordaux R, Gilbert C. Analyzing Horizontal Transfer of Transposable Elements on a Large Scale: Challenges and Prospects. Bioessays [Internet]. 2018;40. https://doi.org/10.1002/bies.201700177.

  15. Ivancevic AM, Kortschak RD, Bertozzi T, Adelson DL. Horizontal transfer of BovB and L1 retrotransposons in eukaryotes. Genome Biol. 2018;19:85.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Gilbert C, Hernandez SS, Flores-Benabib J, Smith EN, Feschotte C. Rampant horizontal transfer of SPIN transposons in squamate reptiles. Mol Biol Evol. 2012;29:503–15.

    Article  CAS  PubMed  Google Scholar 

  17. Walsh AM, Kortschak RD, Gardner MG, Bertozzi T, Adelson DL. Widespread horizontal transfer of retrotransposons. Proc Natl Acad Sci U S A. 2013;110:1012–6.

    Article  CAS  PubMed  Google Scholar 

  18. Zhang H-H, Peccoud J, Xu M-R-X, Zhang X-G, Gilbert C. Horizontal transfer and evolution of transposable elements in vertebrates. Nat Commun. 2020;11:1362.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Hellsten U, Harland RM, Gilchrist MJ, Hendrix D, Jurka J, Kapitonov V, et al. The genome of the western clawed Frog Xenopus tropicalis. Science. 2010;328:633.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Galbraith JD, Ludington AJ, Suh A, Sanders KL, Adelson DL. New environment, new invaders - repeated horizontal transfer of LINEs to sea snakes [Internet]. https://doi.org/10.1101/2020.02.27.968685.

  21. Kordis D. Transposable elements in reptilian and avian (sauropsida) genomes. Cytogenet Genome Res. 2009;127:94–111.

    Article  CAS  PubMed  Google Scholar 

  22. Thomson RC, Spinks PQ, Shaffer HB. A global phylogeny of turtles reveals a burst of climate-associated diversification on continental margins. Proc Natl Acad Sci U S A [Internet]. 2021;118. https://doi.org/10.1073/pnas.2012215118.

  23. Engels WR, Johnson-Schlitz DM, Eggleston WB, Sved J. High-frequency P element loss in Drosophila is homolog dependent. Cell. 1990;62:515–25.

    Article  CAS  PubMed  Google Scholar 

  24. Metzger MJ, Paynter AN, Siddall ME, Goff SP. Horizontal transfer of retrotransposons between bivalves and other aquatic species of multiple phyla. Proc Natl Acad Sci U S A. 2018;115:E4227–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Dunemann SM, Wasmuth JD. Horizontal transfer of a retrotransposon between parasitic nematodes and the common shrew. Mob DNA. 2019;10:24.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Suh A, Witt CC, Menger J, Sadanandan KR, Podsiadlowski L, Gerth M, et al. Ancient horizontal transfers of retrotransposons between birds and ancestors of human pathogenic nematodes. Nat Commun. 2016;7:11396.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Palumbo E, Cassano MJ, Alcalde L, Diaz JI. Seasonal variation of Hedruris dratini (Nematoda) parasitizing Hydromedusa tectifera (Chelidae), with focus on host’s torpor state. BMC Zoology [Internet]. 2021 [cited 2023 Dec 4];6. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10127356/.

  28. Khan RA. Host-Parasite Interactions in Some Fish Species. J Parasitol Res [Internet]. 2012 [cited 2024 Jan 22];2012. https://doi.org/10.1155/2012/237280.

  29. Arfuso F, Gaglio G, Ferrara MC, Abbate F, Giannetto S, Brianti E. First record of infestation by nasal leeches, Limnatis Nilotica (Hirudinida, Praobdellidae), from cattle in Italy. J Vet Med Sci. 2019;81:1419.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Turner WC, Kamath PL, van Heerden H, Huang Y-H, Barandongo ZR, Bruce SA et al. The roles of environmental variation and parasite survival in virulence–transmission relationships. Royal Society Open Science [Internet]. 2021 [cited 2023 Dec 4];8. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8170194/.

  31. Eberhard ML, Yabsley MJ, Zirimwabagabo H, Bishop H, Cleveland CA, Maerz JC, et al. Possible role of Fish and Frogs as paratenic hosts of Dracunculus medinensis, Chad. Emerg Infect Dis. 2016;22:1428.

    Article  PubMed  PubMed Central  Google Scholar 

  32. McDonough JM, Gleason LN. Histopathology in the rainbow darter, Etheostoma caeruleum, resulting from infections with the acanthocephalans, Pomphorhynchus bulbocolli and Acanthocephalus dirus. J Parasitol. 1981;67:403–9.

    Article  CAS  PubMed  Google Scholar 

  33. Zeng L, Kortschak RD, Raison JM, Bertozzi T, Adelson DL. Superior ab initio identification, annotation and characterisation of TEs and segmental duplications from genome assemblies. PLoS ONE. 2018;13:e0193588.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Zhang Z, Schwartz S, Wagner L, Miller W. A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000;7:203–14.

    Article  CAS  PubMed  Google Scholar 

  35. O’Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016;44:D733–45.

    Article  PubMed  Google Scholar 

  36. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Burge C, Karlin S. Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997;268:78–94.

    Article  CAS  PubMed  Google Scholar 

  38. Bao W, Kojima KK, Kohany O. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mob DNA. 2015;6:11.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26:1641–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Letunic I, Bork P. Interactive tree of life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 2019;47:W256–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Hedges SB, Dudley J, Kumar S. TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics. 2006;22:2971–2.

    Article  CAS  PubMed  Google Scholar 

  44. Di Tommaso P, Moretti S, Xenarios I, Orobitg M, Montanyola A, Chang J-M, et al. T-Coffee: a web server for the multiple sequence alignment of protein and RNA sequences using structural information and homology extension. Nucleic Acids Res. 2011;39:W13–7.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46:W296–303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Krumsiek J, Arnold R, Rattei T. Gepard: a rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics. 2007;23:1026–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Diane Barton for discussions about parasites and Terry Bertozzi for critical reading and helpful suggestions. We would also like to thank our lab for their ongoing support and input while putting together this manuscript.

Funding

This research was funded by the University of Adelaide.

Author information

Authors and Affiliations

Authors

Contributions

N.T.H, J.D.G, and D.L.A designed research; J.D.G provided scripts; N.T.H and D.L.A wrote the paper with input from J.D.G. All authors reviewed the manuscript.

Corresponding author

Correspondence to David L. Adelson.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hassan, N.T., Galbraith, J.D. & Adelson, D.L. Multiple horizontal transfer events of a DNA transposon into turtles, fishes, and a frog. Mobile DNA 15, 7 (2024). https://doi.org/10.1186/s13100-024-00318-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13100-024-00318-9

Keywords