Multiple Mariner transposons in flatworms and hydras are related to those of insects

maT and mosquito transposons in cnidarians: evolutionary history and intraspecific differences

Transposable elements exert a significant effect on the size and structure of eukaryotic genomes. Tc1/mariner superfamily elements represent the widely distributed and highly variable group of DNA transposons. Tc1/mariner elements include TLE/DD34-38E, MLE/DD34D, maT/DD37D, Visitor/DD41D, Guest/DD39D, mosquito/DD37E, and L18/DD37E families, all of which are well or less scarcely studied. However, more detailed research into the patterns of prevalence and diversity of Tc1/mariner transposons enables one to better understand the coevolution of the TEs and the eukaryotic genomes. We performed a detailed analysis of the maT/DD37D family in Cnidaria. The study of 77 genomic assemblies demonstrated that maT transposons are found in a limited number of cnidarian species belonging to classes Cubozoa (1 species), Hydrozoa (3 species) и Scyphozoa (5 species) only. The identified TEs were classified into 5 clades, with the representatives from Pelagiidae (class Scyphozoa) forming a separate clade of maT transposons, which has never been described previously. The potentially functional copies of maT transposons were identified in the hydrae. The phylogenetic analysis and the studies of distribution among the taxons and the evolutionary dynamics of the elements suggest that maT transposons of the cnidarians are the descendants of several independent invasion events occurring at different periods of time. We also established that the TEs of mosquito/DD37E family are found in Hydridae (class Hydrozoa) only. A comparison of maT and mosquito prevalence in two genomic assemblies of Hydra viridissima revealed obvious differences, thus demonstrating that each individual organism might carry a unique mobilome pattern. The results of the presented research make us better understand the diversity and evolution of Tc1/mariner transposons and their effect on the eukaryotic genomes.

Insights on mauritiana-like Elements Diversity in Mayetiola destructor and M. hordei (Diptera: Cecidomyiidae)

Mariner-like elements (MLEs) are class II transposons belonging to the Tc1-mariner family, that have successfully invaded many insect genomes. In the current study, the availability of the Hessian fly Mayetiola destructor genome has enabled us to perform in silico analysis of MLEs using as query the previously described mariner element (Desmar1) belonging to mauritiana subfamily. Eighteen mauritiana-like elements were detected and were clustered into three main groups named Desmar1-like, MauCons1 and MauCons2. Subsequently, in vitro analysis was carried out to investigate mauritiana-like elements in M. destructor as well as in Mayetiola hordei using primers designed from TIRs of the previously identified MLEs. PCR amplifications were successful and a total of 12 and 17 mauritiana-like elements were discovered in M. destructor and M. hordei, respectively. Sequence analyses of mauritiana-like elements obtained in silico and in vitro have showed that MauCons1 and MauCons2 elements share low similarity with Desmar1 ranging from 50% to 55% suggesting different groups under mauritiana subfamily have invaded the genomes of M. destructor and M. hordei. These groups are likely inherited by vertical transmission that subsequently underwent different evolutionary histories. This work describes new mauritiana-like elements in M. destructor that are distinct from the previouslydiscovered Desmar1 and provides the first evidence of MLEs belonging to mauritiana subfamily in M. hordei.

The evolutionary history of mariner elements in stalk-eyed flies reveals the horizontal transfer of transposons from insects into the genome of the cnidarian Hydra vulgaris

The stalk-eyed flies (Diopsidae, Diptera) are a family of approximately 100 species of calypterate dipterans, characterised by extended head capsules. Species within the family have previously been shown to possess six subfamilies of mariner transposons, with nucleotide substitution patterns suggesting that at least two subfamilies are currently active. The vertumnana subfamily has been shown to have been involved in a horizontal transfer event involving Diopsidae and a second dipteran family in the Tephritidae. Presented here are cloned and sequenced mariner elements from three further diopsid species, in addition to a bioinformatic analysis of mariner elements identified in transcriptomic and genomic data from the genus Teleopsis. The newly identified mariner elements predominantly fall into previously recognised subfamilies, however the publicly available Teleopsis data also revealed a novel subfamily. Three of the seven identified subfamilies are shown to have undergone horizontal transfer, two of which appear to involve diopsid donor species. One recipient group of a diopsid mariner is the Bactrocera genus of tephritid flies, the transfer of which was previously proposed in an earlier study of diopsid mariner elements. The second horizontal transfer, of the mauritiana subfamily, can be traced from the Teleopsis genus to the cnidarian Hydra vulgaris. The mauritiana elements are shown to be active in the recipient H. vulgaris and transposase expression is observed in all body tissues examined in both species. The increased diversity of diopsid mariner elements points to a minimum of four subfamilies being present in the ancestral genome. Both vertical inheritance and stochastic loss of TEs have subsequently occurred within the diopsid radiation. The TE complement of H. vulgaris contains at least two mariner subfamilies of insect origin. Despite the phylogenetic distance between donor and recipient species, both subfamilies are shown to be active and proliferating within H. vulgaris.

Phylogenetic analysis of the Tc1/mariner superfamily reveals the unexplored diversity of pogo-like elements

Background

Tc1/mariner transposons are widespread DNA transposable elements (TEs) that have made important contributions to the evolution of host genomic complexity in metazoans. However, the evolution and diversity of the Tc1/mariner superfamily remains poorly understood. Following recent developments in genome sequencing and the availability of a wealth of new genomes, Tc1/mariner TEs have been identified in many new taxa across the eukaryotic tree of life. To date, the majority of studies focussing on Tc1/mariner elements have considered only a single host lineage or just a small number of host lineages. Thus, much remains to be learnt about the evolution of Tc1/mariner TEs by performing analyses that consider elements that originate from across host diversity.

Results

We mined the non-redundant database of NCBI using BLASTp searches, with transposase sequences from a diverse set of reference Tc1/mariner elements as queries. A total of 5158 Tc1/mariner elements were retrieved and used to reconstruct evolutionary relationships within the superfamily. The resulting phylogeny is well resolved and includes several new groups of Tc1/mariner elements. In particular, we identify a new family of plant-genome restricted Tc1/mariner elements, which we call PlantMar. We also show that the pogo family is much larger and more diverse than previously appreciated, and we review evidence for a potential revision of its status to become a separate superfamily.

Conclusions

Our study provides an overview of Tc1-mariner phylogeny and summarises the impressive diversity of Tc1-mariner TEs among sequenced eukaryotes. Tc1/mariner TEs are successful in a wide range of eukaryotes, especially unikonts (the taxonomic supergroup containing Amoebozoa, Opisthokonta, Breviatea, and Apusomonadida). In particular, ecdysozoa, and especially arthropods, emerge as important hosts for Tc1/mariner elements (except the PlantMar family). Meanwhile, the pogo family, which is by far the largest Tc1/mariner family, also includes many elements from fungal and chordate genomes. Moreover, there is evidence of the repeated exaptation of pogo elements in vertebrates, including humans, in addition to the well-known example of CENP-B. Collectively, our findings provide a considerable advancement in understanding of Tc1/mariner elements, and more generally they suggest that much work remains to improve understanding of the diversity and evolution of DNA TEs.

Diversity of <i>mariner</i>-like elements in Orthoptera

Mariner-like elements (MLEs) are among the most widespread DNA transposable elements in eukaryotes. Insects were the first organisms in which MLEs were identified, however the diversity of MLEs in the insect order Orthoptera has not yet been addressed. In the present study, we explore the diversity of MLEs elements in 16 species of Orthoptera belonging to three infraorders, Acridoidea (Caelifera), Grylloidea (Ensifera), and Tettigoniidea (Ensifera) by combining data mined from computational analysis of sequenced degenerative PCR MLE amplicons and available Orthoptera genomic scaffolds. In total, 75 MLE lineages (Ortmar) were identified in all the studied genomes. Automatic phylogeny-based classification suggested that the current known variability of MLEs can be assigned to seven statistically well-supported phylogenetic clusters (I–VII), and the identified Orthoptera lineages were distributed among all of them. The majority of the lineages (36 out of 75) belong to cluster I; 20 belong to cluster VI; and seven, six, four, one and one lineages belong to clusters II, IV, VII, III, and V, respectively. Two of the clusters (II and IV) were composed of a single Orthoptera MLE lineage each (Ortmar37 and Ortmar45, respectively) which were distributed in the vast majority of the studied Orthoptera genomes. Finally, for 16 Orthoptera MLE lineages, horizontal transfer from the distantly related taxa belonging to other insect orders may have occurred. We believe that our study can serve as a basis for future researches on the diversity, distribution, and evolution of MLEs in species of other taxa that are still lacking the sequenced genomes.

Genetic exchange in eukaryotes through horizontal transfer: connected by the mobilome

Background

All living species contain genetic information that was once shared by their common ancestor. DNA is being inherited through generations by vertical transmission (VT) from parents to offspring and from ancestor to descendant species. This process was considered the sole pathway by which biological entities exchange inheritable information. However, Horizontal Transfer (HT), the exchange of genetic information by other means than parents to offspring, was discovered in prokaryotes along with strong evidence showing that it is a very important process by which prokaryotes acquire new genes.

Main body

For some time now, it has been a scientific consensus that HT events were rare and non-relevant for evolution of eukaryotic species, but there is growing evidence supporting that HT is an important and frequent phenomenon in eukaryotes as well.

Conclusion

Here, we will discuss the latest findings regarding HT among eukaryotes, mainly HT of transposons (HTT), establishing HTT once and for all as an important phenomenon that should be taken into consideration to fully understand eukaryotes genome evolution. In addition, we will discuss the latest development methods to detect such events in a broader scale and highlight the new approaches which should be pursued by researchers to fill the knowledge gaps regarding HTT among eukaryotes.

Social evolution. Genomic signatures of evolutionary transitions from solitary to group living

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

An evaluation of the ecological relationship between Drosophila species and their parasitoid wasps as an opportunity for horizontal transposon transfer

Evidences of horizontal transfer, the exchange of genetic material between reproductively isolated species, have accumulated over the last decades, including for multicellular eukaryotic organisms. However, the mechanisms and ecological relationships that promote such phenomenon is still poorly known. Host-parasite interaction is one type of relationship usually pointed in the literature that could potentially increase the probability of the horizontal transfer between species, because the species involved in such relationships are generally in close contact. Transposable elements, which are well-known genomic parasites, are DNA entities that tend to be involved in horizontal transfer due to their ability to mobilize between different genomic locations. Using Drosophila species and their parasitoid wasps as a host-parasite model, we evaluated the hypothesis that horizontal transposon transfers (HTTs) are more frequent in this set of species than in species that do not exhibit a close ecological and phylogenetic relationship. For this purpose, we sequenced two sets of species using a metagenomic and single-species genomic sampling approach through next-generation DNA sequencing. The first set was composed of five generalist Drosophila (D. maculifrons, D. bandeirantorum, D. polymorpha, D. mercatorum and D. willistoni) species and their associated parasitoid wasps, whereas the second set was composed of D. incompta, which is a flower specialist species, and its parasitoid wasp. We did not find strong evidence of HTT in the two sets of Drosophila and wasp parasites. However, at least five cases of HTT were observed between the generalist and specialist Drosophila species. Moreover, we detected an HT event involving a Wolbachia lineage between generalist and specialist species, indicating that these endosymbiotic bacteria could play a role as HTT vectors. In summary, our results do not support the hypothesis of prevalent HTT between species with a host-parasite relationship, at least for the studied wasp-Drosophila pairs. Moreover, it suggests that other mechanisms or parasites are involved in promoting HTT between Drosophila species as the Wolbachia endosymbiotic bacteria.

Trematode Himasthla elongata mariner element (Hemar): structure and applications

We cloned and analyzed Hemar1-the full-length mariner of Himasthla elongata. Hemar1 amount and distribution in the genome is typical for the transposable elements. Hemar1 closest relatives found in databases are the mariner-like element (MLE) of Girardia tigrina with 88% similarity in the most conserved transposase domain and Cemar1 of Caenorhabditis elegans with the most similar inverted terminal repeats. Hydra's (Cnidaria) MLE are the next in similarity to Hemar1. We checked whether sequences similar to Hemar1 exist in intermediate and definitive hosts of the parasitic trematode and did not find obvious similarity. This fact, together with the data of Hemar1 evolutionary position, argues against recent MLE-mediated horizontal transfer in this parasite-host model. Our results demonstrate that H. elongata generates genomic variability in asexual parthenogenetic generations within the snail. Transposon insertional display based on full-length sequence showed that Hemar1 could be located in the regions involved in generating clonal diversity in rediae and cercariae, that is, trematode parthenitae.

Horizontal transposon transfer in eukarya: detection, bias, and perspectives

The genetic similarity observed among species is normally attributed to the existence of a common ancestor. However, a growing body of evidence suggests that the exchange of genetic material is not limited to the transfer from parent to offspring but can also occur through horizontal transfer (HT). Transposable elements (TEs) are DNA fragments with an innate propensity for HT; they are mobile and possess parasitic characteristics that allow them to exist and proliferate within host genomes. However, horizontal transposon transfer (HTT) is not easily detected, primarily because the complex TE life cycle can generate phylogenetic patterns similar to those expected for HTT events. The increasingly large number of new genome projects, in all branches of life, has provided an unprecedented opportunity to evaluate the TE content and HTT events in these species, although a standardized method of HTT detection is required before trends in the HTT rates can be evaluated in a wide range of eukaryotic taxa and predictions about these events can be made. Thus, we propose a straightforward hypothesis test that can be used by TE specialists and nonspecialists alike to discriminate between HTT events and natural TE life cycle patterns. We also discuss several plausible explanations and predictions for the distribution and frequency of HTT and for the inherent biases of HTT detection. Finally, we discuss some of the methodological concerns for HTT detection that may result in the underestimation and overestimation of HTT rates during eukaryotic genome evolution.

The rhizome of life: what about metazoa?

The increase in huge number of genomic sequences in recent years has contributed to various genetic events such as horizontal gene transfer (HGT), gene duplication and hybridization of species. Among them HGT has played an important role in the genome evolution and was believed to occur only in Bacterial and Archaeal genomes. As a result, genomes were found to be chimeric and the evolution of life was represented in different forms such as forests, networks and species evolution was described more like a rhizome, rather than a tree. However, in the last few years, HGT has also been evidenced in other group such as metazoa (for example in root-knot nematodes, bdelloid rotifers and mammals). In addition to HGT, other genetic events such as transfer by retrotransposons and hybridization between more closely related lineages are also well established. Therefore, in the light of such genetic events, whether the evolution of metazoa exists in the form of a tree, network or rhizome is highly questionable and needs to be determined. In the current review, we will focus on the role of HGT, retrotransposons and hybridization in the metazoan evolution.

The evolutionary history of mariner-like elements in Neotropical drosophilids

The evolutionary history of mariner-like elements (MLEs) in 49 mainly Neotropical drosophilid species is described. So far, the investigations about the distribution of MLEs were performed mainly using hybridization assays with the Mos1 element (the first mariner active element described) in a widely range of drosophilid species and these sequences were found principally in species that arose in Afrotropical and Sino-Indian regions. Our analysis in mainly Neotropical drosophilid species shows that twenty-three species presented MLEs from three different subfamilies in their genomes: eighteen species had MLEs from subfamily mellifera, fifteen from subfamily mauritiana and three from subfamily irritans. Eleven of these species exhibited elements from more than one subfamily in their genome. In two subfamilies, the analyzed coding region was uninterrupted and contained conserved catalytic motifs. This suggests that these sequences were probably derived from active elements. The species with these putative active elements are Drosophila mediopunctata and D. busckii for the mauritiana subfamily, and D. paramediostriata for the mellifera subfamily. The phylogenetic analysis of MLE, shows a complex evolutionary pattern, exhibiting vertical transfer, stochastic loss and putative events of horizontal transmission occurring between different Drosophilidae species, and even those belonging to more distantly related taxa such as Bactrocera tryoni (Tephritidae family), Sphyracephala europaea (Diopsoidea superfamily) and Buenoa sp. (Hemiptera order). Moreover, our data show that the distribution of MLEs is not restricted to Afrotropical and Sino-Indian species. Conversely, these TEs are also widely distributed in drosophilid species arisen in the Neotropical region.

Next-generation sequencing reveals recent horizontal transfer of a DNA transposon between divergent mosquitoes

Horizontal transfer of genetic material between complex organisms often involves transposable elements (TEs). For example, a DNA transposon mariner has been shown to undergo horizontal transfer between different orders of insects and between different phyla of animals. Here we report the discovery and characterization of an ITmD37D transposon, MJ1, in Anopheles sinensis. We show that some MJ1 elements in Aedes aegypti and An. sinensis contain intact open reading frames and share nearly 99% nucleotide identity over the entire transposon, which is unexpectedly high given that these two genera had diverged 145–200 million years ago. Chromosomal hybridization and TE-display showed that MJ1 copy number is low in An. sinensis. Among 24 mosquito species surveyed, MJ1 is only found in Ae. aegypti and the hyrcanus group of anopheline mosquitoes to which An. sinensis belongs. Phylogenetic analysis is consistent with horizontal transfer and provides the basis for inference of its timing and direction. Although report of horizontal transfer of DNA transposons between higher eukaryotes is accumulating, our analysis is one of a small number of cases in which horizontal transfer of nearly identical TEs among highly divergent species has been thoroughly investigated and strongly supported. Horizontal transfer involving mosquitoes is of particular interest because there are ongoing investigations of the possibility of spreading pathogen-resistant genes into mosquito populations to control malaria and other infectious diseases. The initial indication of horizontal transfer of MJ1 came from comparisons between a 0.4x coverage An. sinensis 454 sequence database and available TEs in mosquito genomes. Therefore we have shown that it is feasible to use low coverage sequencing to systematically uncover horizontal transfer events. Expanding such efforts across a wide range of species will generate novel insights into the relative frequency of horizontal transfer of different TEs and provide the evolutionary context of these lateral transfer events.

Occurrence and abundance of a mariner-like element in freshwater and terrestrial planarians (Platyhelminthes, Tricladida) from southern Brazil

Transposable elements are DNA sequences present in all the large phylogenetic groups, both capable of changing position within the genome and constituting a significant part of eukaryotic genomes. The mariner family of transposons is one of the few which occurs in a wide variety of taxonomic groups, including freshwater planarians. Nevertheless, so far only five planarian species have been reported to carry mariner-like elements (MLEs), although several different species have been investigated. Regarding the number of copies of MLEs, Girardia tigrina is the only planarian species in which this has been evaluated, with an estimation of 8,000 copies of the element per haploid genome. Preliminary results obtained in our laboratory demonstrated that MLE is found in a large number of different species of planarians, including terrestrial. With this in mind, the aim was to evaluate the occurrence and estimate the number of MLE copies in different planarian species collected in south Brazil. Twenty-eight individuals from 15 planarian species were analyzed. By using PCR and the hybridization of nucleic acids, it was found that MLE was present in all the analyzed species, the number of copies being high, probably over 10³ per haploid genome.

GalEa retrotransposons from galatheid squat lobsters (Decapoda, Anomura) define a new clade of Ty1/copia-like elements restricted to aquatic species

Crustacean species have not been examined in great detail for their transposable elements content. Here we focus on galatheid crabs, which are one of the most diverse and widespread taxonomic groups of Decapoda. Ty1/copia retrotransposons are a diverse and taxonomically dispersed group. Using degenerate primers, we isolated several DNA fragments that show homology with Ty1/copia retroelements reverse transcriptase gene. We named the corresponding elements from which they originated GalEa1 to GalEa3 and analyzed one of them further by isolating various clones containing segments of GalEa1. This is the first LTR retrotransposon described in crustacean genome. Nucleotide sequencing of the clones revealed that GalEa1 has LTRs (124 bp) and that the internal sequence (4,421 bp) includes a single large ORF containing gag and pol regions. Further screening identified highly related elements in six of the nine galatheid species studied. By performing BLAST searches on genome databases, we could also identify GalEa-like elements in some fishes and Urochordata genomes. These elements define a new clade of Ty1/copia retrotransposons that differs from all other Ty1/copia elements and that seems to be restricted to aquatic species.

Molecular characterization and phylogenetic position of a new mariner-like element in the coastal crab, Pachygrapsus marmoratus

Mariner-like elements (MLEs) are class-II transposable elements that move within the genome of their hosts by means of a DNA-mediated "cut and paste" mechanism. MLEs have been identified in several organisms, from most of the phyla. Nevertheless, only a few of the sequences characterized contain an intact open reading frame. Investigation of the genome of a coastal crab, Pachygrapsus marmoratus, has identified nine Pacmmar elements, two of which have an open reading frame encoding a putatively functional transposase. Nucleic acid analyses and comparison with the previous data showed that the GC contents of MLEs derived from coastal organisms such as P. marmoratus are significantly higher than those of terrestrial MLEs and significantly lower than those of hydrothermal ones. Furthermore, molecular phylogeny analyses have shown that Pacmmar elements constitute a new lineage of the irritans subfamily within the mariner family.

Species sympatry and horizontal transfers of Mariner transposons in marine crustacean genomes

Mariner-like elements (MLEs) have been widely detected in terrestrial species. The first complete MLE isolated from a marine invertebrate was detected in the genome of the hydrothermal crab Bythograea thermydron by Halaimia-Toumi et al. [Halaimia-Toumi, N., Casse, N., Demattei, M.V., Renault, S., Pradier, E., Bigot, Y., Laulier, M., 2004. The GC-rich transposon Bytmar1 from the deep-sea hydrothermal crab, Bythograea thermydron, may encode three transposase isoforms from a single ORF. J. Mol. Evol. 59, 747-760] and called Bytmar1. Here, we report the isolation of three new Bytmar1 relatives from the genomes of one hydrothermal amphipod Ventiella sulfuris (Vensmar1) and two coastal crustacea, Maia brachydactila (Maibmar1) and Cancer pagurus (Canpmar1). Like Bytmar1, these MLEs have an unusually high GC content, a high CpG ratio, and a low TpA ratio. Their consensus sequence encodes a transposase that is preceded by an N-flag, as in Bytmar1, which could be a marine feature. Only one of the 19 clones obtained, Vensmar1.3, encoded for a full-length transposase. The phylogenetic analyses revealed that all these Bytmar1-related elements can be differentiated into two clusters, corresponding to the coastal or hydrothermal origin of their hosts. They also confirmed that the irritans sub-family comprises at least four lineages that seem to depend on the taxonomical position and habitat of their hosts. Finally, we observed that elements coding for two potentially complete transposases exhibiting 99.5% similarity, Bytmar1.11 and Vensmar1.3, were present in the genome of two distantly related hydrothermal crustacea, one Amphipod and one Decapod. The hypothesis of horizontal transfers is discussed in the light of the sequence similarities observed.

The GC-rich transposon Bytmar1 from the deep-sea hydrothermal crab, Bythograea thermydron, may encode three transposase isoforms from a single ORF

Mariner-like elements (MLEs) are classII transposons with highly conserved sequence properties and are widespread in the genome of animal species living in continental environments. We describe here the first full-length MLE found in the genome of a marine crustacean species, the deep-sea hydrothermal crab Bythograea thermydron (Crustacea), named Bytmar1. A comparison of its sequence features with those of the MLEs contained in the genomes of continental species reveals several distinctive characteristics. First, Bytmar1 elements contains an ORF that may encode three transposase isoforms 349, 379, and 398 amino acids (aa) in long. The two biggest proteins are due to the presence of a 30- and 49-aa flag, respectively, at the N-terminal end of the 349-aa cardinal MLE transposase. Their GC contents are also significantly higher than those found in continental MLEs. This feature is mainly due to codon usage in the transposase ORF and directly interferes with the curvature propensities of the Bytmar1 nucleic acid sequence. Such an elevated GC content may interfere with the ability of Bytmar 1 to form an excision complex and, in consequence, with its efficiency to transpose. Finally, the origin of these characteristics and their possible consequences on transposition efficiency are discussed.

The Sinbad retrotransposon from the genome of the human blood fluke, Schistosoma mansoni, and the distribution of related Pao-like elements

Background

Of the major families of long terminal repeat (LTR) retrotransposons, the Pao/BEL family is probably the least well studied. It is becoming apparent that numerous LTR retrotransposons and other mobile genetic elements have colonized the genome of the human blood fluke, Schistosoma mansoni.

Results

A proviral form of Sinbad, a new LTR retrotransposon, was identified in the genome of S. mansoni. Phylogenetic analysis indicated that Sinbad belongs to one of five discreet subfamilies of Pao/BEL like elements. BLAST searches of whole genomes and EST databases indicated that members of this clade occurred in species of the Insecta, Nematoda, Echinodermata and Chordata, as well as Platyhelminthes, but were absent from all plants, fungi and lower eukaryotes examined. Among the deuterostomes examined, only aquatic species harbored these types of elements. All four species of nematode examined were positive for Sinbad sequences, although among insect and vertebrate genomes, some were positive and some negative. The full length, consensus Sinbad retrotransposon was 6,287 bp long and was flanked at its 5'- and 3'-ends by identical LTRs of 386 bp. Sinbad displayed a triple Cys-His RNA binding motif characteristic of Gag of Pao/BEL-like elements, followed by the enzymatic domains of protease, reverse transcriptase (RT), RNAseH, and integrase, in that order. A phylogenetic tree of deduced RT sequences from 26 elements revealed that Sinbad was most closely related to an unnamed element from the zebrafish Danio rerio and to Saci-1, also from S. mansoni. It was also closely related to Pao from Bombyx mori and to Ninja of Drosophila simulans. Sinbad was only distantly related to the other schistosome LTR retrotransposons Boudicca, Gulliver, Saci-2, Saci-3, and Fugitive, which are gypsy-like. Southern hybridization and bioinformatics analyses indicated that there were about 50 copies of Sinbad in the S. mansoni genome. The presence of ESTs representing transcripts of Sinbad in numerous developmental stages of S. mansoni along with the identical 5'- and 3'-LTR sequences suggests that Sinbad is an active retrotransposon.

Conclusion

Sinbad is a Pao/BEL type retrotransposon from the genome of S. mansoni. The Pao/BEL group appears to be comprised of at least five discrete subfamilies, which tend to cluster with host species phylogeny. Pao/BEL type elements appear to have colonized only the genomes of the Animalia. The distribution of these elements in the Ecdysozoa, Deuterostomia, and Lophotrochozoa is discontinuous, suggesting horizontal transmission and/or efficient elimination of Pao-like mobile genetic elements from some genomes.

A potentially functional mariner transposable element in the protist Trichomonas vaginalis

Mariner transposable elements encoding a D,D34D motif-bearing transposase are characterized by their pervasiveness among, and exclusivity to, animal phyla. To date, several hundred sequences have been obtained from taxa ranging from cnidarians to humans, only two of which are known to be functional. Related transposons have been identified in plants and fungi, but their absence among protists is noticeable. Here, we identify and characterize Tvmar1, the first representative of the mariner family to be found in a species of protist, the human parasite Trichomonas vaginalis. This is the first D,D34D element to be found outside the animal kingdom, and its inclusion in the mariner family is supported by both structural and phylogenetic analyses. Remarkably, Tvmar1 has all the hallmarks of a functional element and has recently expanded to several hundred copies in the genome of T. vaginalis. Our results show that a new potentially active mariner has been found that belongs to a distinct mariner lineage and has successfully invaded a nonanimal, single-celled organism. The considerable genetic distance between Tvmar1 and other mariners may have valuable implications for the design of new, high-efficiency vectors to be used in transfection studies in protists.

Medfly transposable elements: diversity, evolution, genomic impact and possible applications

The medfly genome has been shown to contain a rich assortment of transposable elements from the mariner, Tc1, hAT and gypsy/Ty3 families. These elements display different levels of diversity, abundance and distribution in the genome. The presence of actively transposing elements in the medfly genome is revealed by hybrid dysgenesis phenomena, insertion site polymorphisms and other genetic instabilities. The medfly has been a target of transformation studies involving the exogenous elements Minos, Hermes and piggyBac from three families. The presence of active endogenous homologous elements can have important implications for the stability of such transgenic lines. The potential applications of endogenous elements for medfly population analysis and control are discussed.

Boudicca, a retrovirus-like long terminal repeat retrotransposon from the genome of the human blood fluke Schistosoma mansoni

The genome of Schistosoma mansoni contains a proviral form of a retrovirus-like long terminal repeat (LTR) retrotransposon, designated BOUDICCA: Sequence and structural characterization of the new mobile genetic element, which was found in bacterial artificial chromosomes prepared from S. mansoni genomic DNA, revealed the presence of three putative open reading frames (ORFs) bounded by direct LTRs of 328 bp in length. ORF1 encoded a retrovirus-like major homology region and a Cys/His box motif, also present in Gag polyproteins of related retrotransposons and retroviruses. ORF2 encoded enzymatic domains and motifs characteristic of a retrovirus-like polyprotein, including aspartic protease, reverse transcriptase, RNase H, and integrase, in that order, a domain order similar to that of the gypsy/Ty3 retrotransposons. An additional ORF at the 3' end of the retrotransposon may encode an envelope protein. Phylogenetic comparison based on the reverse transcriptase domain of ORF2 confirmed that Boudicca was a gypsy-like retrotransposon and showed that it was most closely related to CsRn1 from the Oriental liver fluke Clonorchis sinensis and to kabuki from Bombyx mori. Bioinformatics approaches together with Southern hybridization analysis of genomic DNA of S. mansoni and the screening of a bacterial artificial chromosome library representing approximately 8-fold coverage of the S. mansoni genome revealed that numerous copies of Boudicca were interspersed throughout the schistosome genome. By reverse transcription-PCR, mRNA transcripts were detected in the sporocyst, cercaria, and adult developmental stages of S. mansoni, indicating that Boudicca is actively transcribed in this trematode.

Neutral evolution of ten types of mariner transposons in the genomes of Caenorhabditis elegans and Caenorhabditis briggsae

Ten types of mariner transposable elements (232 individual sequences) are present in the completed genomic DNA sequence of Caenorhabditis elegans and the partial sequence of Caenorhabditis briggsae. We analyze these replicated instances of mariner evolution and find that elements of a type have evolved within their genomes under no selection on their transposase genes. Seven of the ten reconstructed ancestral mariners carry defective transposase genes. Selection has acted during the divergence of some ancestral elements. The neutrally-evolving mariners are used to analyze the pattern of molecular evolution in Caenorhabditis. There is a significant mutational bias against transversions and significant variation in rates of change across sites. Deletions accumulate at a rate of 0.034 events/bp per substitution/site, with an average size of 166 bp (173 gaps observed). Deletions appear to obliterate preexisting deletions over time, creating larger gaps. Insertions accumulate at a rate of 0.019 events/bp per substitution/site, with an average size of 151 bp (61 events). Although the rate of deletion is lower than most estimates in other species, the large size of deletions causes rapid elimination of neutral DNA: a mariner's "half-life" (the time by which half an element's sequence should have been deleted) is approximately 0.1 subsitutions/site. This high rate of DNA deletion may explain the compact nature of the nematode genome.

Bmmar6, a second mori subfamily mariner transposon from the silkworm moth Bombyx mori

A second member of the divergent mori subfamily of mariner transposons, Bmmar6, is described from the silkworm moth Bombyx mori genome. A confident consensus sequence for Bmmar6 was obtained from a single genomic copy, 17 EST sequences, and the direct sequencing of a 'family' sequence from an amplification of all full-length genomic copies. Bmmar6 is most similar to Bmmar1 in the mori subfamily, which now also includes several fly and nematode transposons. These might be viewed as a discrete family of transposons within the IS630-Tc1-mariner superfamily with a distinctive D,D37D catalytic motif, and another small divergent D,D41D clade is recognized as their sister group of transposons.

Molecular evolutionary analyses of mariners and other transposable elements in fire ants (Hymenoptera: Formicidae)

Screens of a library of genomic DNA made during a recent study of the fire ant Solenopsis invicta revealed the presence of three distinct types of transposable elements (TEs). Two of the recovered sequences showed a high similarity to long-terminal repeat (LTR) retrotransposons, while the third showed a high homology to mariner elements. To investigate the distribution and relationships of mariners in related ants, we PCR-amplified these elements from additional Solenopsis species. Phylogenetic analyses showed that they form a single group within the mauritiana subfamily that is part of a larger clade derived from hymenopteran species. We also present partial sequence data for the two LTR-retrotransposons and describe their phylogenetic affinities.

Developmental signaling in Hydra: what does it take to build a "simple" animal?

Developmental processes in multicellular animals depend on an array of signal transduction pathways. Studies of model organisms have identified a number of such pathways and dissected them in detail. However, these model organisms are all bilaterians. Investigations of the roles of signal transduction pathways in the early-diverging metazoan Hydra have revealed that a number of the well-known developmental signaling pathways were already in place in the last common ancestor of Hydra and bilaterians. In addition to these shared pathways, it appears that developmental processes in Hydra make use of pathways involving a variety of peptides. Such pathways have not yet been identified as developmental regulators in more recently diverged animals. In this review I will summarize work to date on developmental signaling pathways in Hydra and discuss the future directions in which such work will need to proceed to realize the potential that lies in this simple animal.

Reverse transcriptase activity and untranslated region sharing of a new RTE-like, non-long terminal repeat retrotransposon from the human blood fluke, Schistosoma japonicum

A new RTE-like, non-long terminal repeat retrotransposon, termed SjR2, from the human blood fluke, Schistosoma japonicum, is described. SjR2 is approximately 3.9 kb in length and is constituted of a single open reading frame encoding a polyprotein with apurinic/apyrimidinic endonuclease and reverse transcriptase domains. The open reading frame is bounded by 5'- and 3'-terminal untranslated regions and, at its 3'-terminus, SjR2 bears a short (TGAC)(3) repeat. Phylogenetic analyses based on conserved domains of reverse transcriptase or endonuclease revealed that SjR2 belonged to the RTE clade of non-long terminal repeat retrotransposons. Further, SjR2 was homologous, but probably not orthologous, to SR2 from the African blood fluke, Schistosoma mansoni; this RTE-like family of non-long terminal repeat retrotransposons appears to have arisen before the divergence of the extant schistosome species. Hybridisation analyses indicated that approximately 10,000 copies of SjR2 were dispersed throughout the S. japonicum chromosomes, accounting for up to 14% of the nuclear genome. Messenger RNAs encoding the reverse transcriptase and endonuclease domains of SjR2 were detected in several developmental stages of the schistosome, indicating that the retrotransposon was actively replicating within the genome of the parasite. Exploration of the coding and non-coding regions of SjR2 revealed two notable characteristics. First, the recombinant reverse transcriptase domain of SjR2 expressed in insect cells primed reverse transcription of SjR2 mRNA in vitro. By contrast, recombinant SjR2-endonuclease did not appear to cleave schistosome or plasmid DNA. Second, the 5'-untranslated region of SjR2 was >80% identical to the 3'-untranslated region of a schistosome heat shock protein-70 gene (hsp-70) in the antisense orientation, indicating that SjR2-like elements were probably inserted into the non-coding regions of ancestral S. japonicum HSP-70, probably after the species diverged from S. mansoni.

Comparative and functional studies of Drosophila species invasion by the gypsy endogenous retrovirus

Gypsy is an endogenous retrovirus of Drosophila melanogaster. Phylogenetic studies suggest that occasional horizontal transfer events of gypsy occur between Drosophila species. gypsy possesses infective properties associated with the products of the envelope gene that might be at the origin of these interspecies transfers. We report here the existence of DNA sequences putatively encoding full-length Env proteins in the genomes of Drosophila species other than D. melanogaster, suggesting that potentially infective gypsy copies able to spread between sexually isolated species can occur. The ability of gypsy to invade the genome of a new species is conditioned by its capacity to be expressed in the naive genome. The genetic basis for the regulation of gypsy activity in D. melanogaster is now well known, and it has been assigned to an X-linked gene called flamenco. We established an experimental simulation of the invasion of the D. melanogaster genome by gypsy elements derived from other Drosophila species, which demonstrates that these non- D. melanogaster gypsy elements escape the repression exerted by the D. melanogaster flamenco gene.

The genome of the Queensland fruit fly Bactrocera tryoni contains multiple representatives of the mariner family of transposable elements

Representatives of five distinct types of transposable elements of the mariner family were detected in the genomes of the Queensland fruit fly Bactrocera tryoni and its sibling species Bactrocera neohumeralis by phylogenetic analysis of transposase gene fragments. Three mariner types were also found in an additional tephritid, Bactrocera jarvisi. Using genomic library screening and inverse PCR, full-length elements representing the mellifera subfamily (B. tryoni.mar1) and the irritans subfamily (B. tryoni.mar2) were isolated from the B. tryoni genome. Nucleotide consensus sequences for each type were derived from multiple defective copies. Predicted transposase sequences share approximately 23% amino acid identity. B. tryoni.mar1 elements have an estimated copy number of about 900 in the B. tryoni genome, whereas B. tryoni.mar2 element types appear to be present in low copy number.

Mutator transposase is widespread in the grasses

Although the Mutator (Mu) system is well characterized in maize (Zea mays), very little is known about this highly mutagenic system of transposons in other grasses. Mutator is regulated by the MuDR class of elements, which encodes two genes, one of which, mudrA, has similarity to a number of bacterial transposases. Experiments in our laboratory, as well as database searches, demonstrate that mudrA sequences are ubiquitous and diverse in the grasses. In several species it is clear that multiple paralogous elements can be present in a single genome. In some species such as wheat (Triticum aestivum) and rice (Oryza sativa), mudrA-similar sequences are represented in cDNA databases, suggesting the presence of active Mu transposon systems in these species. Further, in rice and in sorghum, mudrA-like genes are flanked by long terminal inverted repeats, as well as the short host sequence direct repeats diagnostic of insertion. Thus, there is ample evidence that systems related to Mu in maize are at least potentially active in a wide variety of grasses. However, the mudrB gene, though important for Mu activity in maize, is not necessarily a component of Mu elements in other grasses.

A current perspective on insect gene transformation

The genetic transformation of non-drosophilid insects is now possible with several systems, with germ-line transformation reported in published and unpublished accounts for about 12 species using four different transposon vectors. For some of these species, transformation can now be considered routine. Other vector systems include viruses and bacterial symbionts that have demonstrated utility in species and applications requiring transient expression, and for some, the potential exists for genomic integration. Many of these findings are quite recent, presenting a dramatic turning point in our ability to study and manipulate agriculturally and medically important insects. This review discusses these findings from the perspective of all the contributions that has made this technology a reality, the research that has yet to be done for its safe and efficient use in a broader range of species, and an overview of the available methodology to effectively utilize these systems.

Transposable elements in sexual and ancient asexual taxa

Sexual reproduction allows deleterious transposable elements to proliferate in populations, whereas the loss of sex, by preventing their spread, has been predicted eventually to result in a population free of such elements [Hickey, D. A. (1982) Genetics 101, 519-531]. We tested this expectation by screening representatives of a majority of animal phyla for LINE-like and gypsy-like reverse transcriptases and mariner/Tc1-like transposases. All species tested positive for reverse transcriptases except rotifers of the class Bdelloidea, the largest eukaryotic taxon in which males, hermaphrodites, and meiosis are unknown and for which ancient asexuality is supported by molecular genetic evidence. Mariner-like transposases are distributed sporadically among species and are present in bdelloid rotifers. The remarkable lack of LINE-like and gypsy-like retrotransposons in bdelloids and their ubiquitous presence in other taxa support the view that eukaryotic retrotransposons are sexually transmitted nuclear parasites and that bdelloid rotifers evolved asexually.

Poised for contagion: evolutionary origins of the infectious abilities of invertebrate retroviruses

Phylogenetic analyses suggest that long-terminal repeat (LTR) bearing retrotransposable elements can acquire additional open-reading frames that can enable them to mediate infection. Whereas this process is best documented in the origin of the vertebrate retroviruses and their acquisition of an envelope (env) gene, similar independent events may have occurred in insects, nematodes, and plants. The origins of env-like genes are unclear, and are often masked by the antiquity of the original acquisitions and by their rapid rate of evolution. In this report, we present evidence that in three other possible transitions of LTR retrotransposons to retroviruses, an envelope-like gene was acquired from a viral source. First, the gypsy and related LTR retrotransposable elements (the insect errantiviruses) have acquired their envelope-like gene from a class of insect baculoviruses (double-stranded DNA viruses with no RNA stage). Second, the Cer retroviruses in the Caenorhabditis elegans genome acquired their envelope gene from a Phleboviral (single ambisense-stranded RNA viruses) source. Third, the Tas retroviral envelope (Ascaris lumricoides) may have been obtained from Herpesviridae (double-stranded DNA viruses, no RNA stage). These represent the only cases in which the env gene of a retrovirus has been traced back to its original source. This has implications for the evolutionary history of retroviruses as well as for the potential ability of all LTR-retrotransposable elements to become infectious agents.

NeSL-1, an ancient lineage of site-specific non-LTR retrotransposons from Caenorhabditis elegans

Phylogenetic analyses of non-LTR retrotransposons suggest that all elements can be divided into 11 lineages. The 3 oldest lineages show target site specificity for unique locations in the genome and encode an endonuclease with an active site similar to certain restriction enzymes. The more "modern" non-LTR lineages possess an apurinic endonuclease-like domain and generally lack site specificity. The genome sequence of Caenorhabditis elegans reveals the presence of a non-LTR retrotransposon that resembles the older elements, in that it contains a single open reading frame with a carboxyl-terminal restriction-like endonuclease domain. Located near the N-terminal end of the ORF is a cysteine protease domain not found in any other non-LTR element. The N2 strain of C. elegans appears to contain only one full-length and several 5' truncated copies of this element. The elements specifically insert in the Spliced leader-1 genes; hence the element has been named NeSL-1 (Nematode Spliced Leader-1). Phylogenetic analysis confirms that NeSL-1 branches very early in the non-LTR lineage and that it represents a 12th lineage of non-LTR elements. The target specificity of NeSL-1 for the spliced leader exons and the similarity of its structure to that of R2 elements leads to a simple model for its expression and retrotransposition.

Explosive invasion of plant mitochondria by a group I intron

Group I introns are mobile, self-splicing genetic elements found principally in organellar genomes and nuclear rRNA genes. The only group I intron known from mitochondrial genomes of vascular plants is located in the cox1 gene of Peperomia, where it is thought to have been recently acquired by lateral transfer from a fungal donor. Southern-blot surveys of 335 diverse genera of land plants now show that this intron is in fact widespread among angiosperm cox1 genes, but with an exceptionally patchy phylogenetic distribution. Four lines of evidence-the intron's highly disjunct distribution, many incongruencies between intron and organismal phylogenies, and two sources of evidence from exonic coconversion tracts-lead us to conclude that the 48 angiosperm genera found to contain this cox1 intron acquired it by 32 separate horizontal transfer events. Extrapolating to the over 13,500 genera of angiosperms, we estimate that this intron has invaded cox1 genes by cross-species horizontal transfer over 1,000 times during angiosperm evolution. This massive wave of lateral transfers is of entirely recent occurrence, perhaps triggered by some key shift in the intron's invasiveness within angiosperms.

Human meiotic recombination products revealed by sequencing a hotspot for homologous strand exchange in multiple HNPP deletion patients

The HNPP (hereditary neuropathy with liability to pressure palsies) deletion and CMT1A (Charcot-Marie-Tooth disease type 1A) duplication are the reciprocal products of homologous recombination events between misaligned flanking CMT1A-REP repeats on chromosome 17p11. 2-p12. A 1.7-kb hotspot for homologous recombination was previously identified wherein the relative risk of an exchange event is 50 times higher than in the surrounding 98.7% identical sequence shared by the CMT1A-REPs. To refine the region of exchange further, we designed a PCR strategy to amplify the recombinant CMT1A-REP from HNPP patients as well as the proximal and distal CMT1A-REPs from control individuals. By comparing the sequences across recombinant CMT1A-REPs to that of the proximal and distal CMT1A-REPs, the exchange was mapped to a 557-bp region within the previously identified 1.7-kb hotspot in 21 of 23 unrelated HNPP deletion patients. Two patients had recombined sequences suggesting an exchange event closer to the mariner-like element previously identified near the hotspot. Five individuals also had interspersed patches of proximal or distal repeat specific DNA sequence indicating potential gene conversion during the exchange of genetic material. Our studies provide a direct observation of human meiotic recombination products. These results are consistent with the hypothesis that minimum efficient processing segments, which have been characterized in Escherichia coli, yeast, and cultured mammalian cells, may be required for efficient homologous meiotic recombination in humans.

Factors affecting transposition of the Himar1 mariner transposon in vitro

Mariner family transposable elements are widespread in animals, but their regulation is poorly understood, partly because only two are known to be functional. These are particular copies of the Dmmar1 element from Drosophila mauritiana, for example, Mos1, and the consensus sequence of the Himar1 element from the horn fly, Haematobia irritans. An in vitro transposition system was refined to investigate several parameters that influence the transposition of Himar1. Transposition products accumulated linearly over a period of 6 hr. Transposition frequency increased with temperature and was dependent on Mg2+ concentration. Transposition frequency peaked over a narrow range of transposase concentration. The decline at higher concentrations, a phenomenon observed in vivo with Mos1, supports the suggestion that mariners may be regulated in part by "overproduction inhibition." Transposition frequency decreased exponentially with increasing transposon size and was affected by the sequence of the flanking DNA of the donor site. A noticeable bias in target site usage suggests a preference for insertion into bent or bendable DNA sequences rather than any specific nucleotide sequences beyond the TA target site.

Molecular evolution of an ancient mariner transposon, Hsmar1, in the human genome

A confident consensus sequence for Hsmar1, the first mariner transposon recognized in the human genome, was generated using three genomic and 15 cDNA sequences. It is thought to represent the ancestrally active copy that invaded an early primate genome. The consensus is 1287 base pairs (bp) long, has 30 bp perfect inverted terminal repeats (ITRs), and encodes a 343 amino acid (aa) mariner transposase. Each copy has diverged from the consensus largely independently of the others and mostly neutrally, and most are now defective. They differ from the consensus by an average of 7.8% in DNA sequence and 7.5 indels per kilobase, both of which values indicate that the copies were formed about 50 Myr ago. On average, only 20% of the 73 surmised CpG hypermutable sites in the consensus remain. A remarkable exception to this loss of functionality is revealed by a set of ten cDNA clones derived from a particular genomic copy that has diverged only 2.4% from the consensus, retained 54% of its hypermutable CpG pairs, and which has a full-length transposase open reading frame. The complete sequence of one of these cDNAs (NIB1543) indicates that the transposase gene of this copy may have been conserved because it is spliced to a human cellular gene encoding a SET domain protein. A specific PCR assay was used to reveal the presence of Hsmar1 copies in all primates examined representing all major lineages, but not in close relatives of primates. PCR fragments cloned and sequenced from a representative sample of primates confirmed that Hsmar1 copies are present in all major lineages, and also revealed another cecropia subfamily mariner in prosimians only, and a third highly divergent mariner present in the greater slow loris Nycticebus coucang. There are about 200 copies of Hsmar1 in the human genome, as well as +/-2400 copies of a derived 80 bp paired ITR structure and +/-4600 copies of solo ITRs. Thus, this transposon had a considerable insertional mutagenic effect on past primate genomes.

Molecular evolution of the second ancient human mariner transposon, Hsmar2, illustrates patterns of neutral evolution in the human genome lineage

A consensus sequence for the second ancient mariner identified in the human genome, Hsmar2, was constructed by majority rule from full-length and partial sequences of 44 of the +/-1000 copies in the genome. This 1300 base pair (bp) consensus has 31 bp imperfect terminal repeats (ITRs) and encodes a 351 amino acid (aa) mariner transposase. The sequence of this transposase has allowed classification of Hsmar2 as a basal lineage of the irritans subfamily of mariners, sharing at most 38% aa identity with other members of the subfamily. The individual copies in the human genome are all highly mutated from the consensus, having suffered numerous small and some large insertions and deletions (indels), including many insertions of S and J subfamily Alu elements. The copies differ, on average, from the consensus by 11.6%, have suffered 11.8 indels per kilobase (kb), and only 3.7% of the 30 hypermutable CpG dinucleotide pairs in the consensus remain intact. This level of divergence indicates that the ancestrally active Hsmar2 element represented by the consensus was present in the human genome lineage about 80 million years (Myr) ago. Each copy has apparently evolved since then largely independently of the others, and with little constraint on its transposase coding capacity. This pattern of molecular evolution fits the current model for mariner transposon evolution. These copies provide multiple independent datasets for evaluating the pattern of neutral evolution in the human genome, for example, they confirm that most indels are very short and that deletions are twice as common as insertions.

The genomes of most animals have multiple members of the Tc1 family of transposable elements

A PCR assay was employed to detect sequence homologous to the transposase gene of the Tc1 family of transposable elements in a wide variety of animals. Amplification products of the appropriate size were obtained from most insects (92 of 108 examined; 85%), most other invertebrates (33 of 43; 77%), and many vertebrates (18 of 36; 50%). Sequencing a sample of cloned PCR products from eight insects, one hydra, and two frogs revealed that each had multiple distinct members of the family in their genomes. In the most extreme case, the horn fly Haematobia irritans yielded evidence of seventeen distinct types of Tc1 family elements. Most of the sequences obtained indicate that the elements are within the range of variation already known from fungi, nematodes, flies, fish and frogs. Some, however, had novel length variants or divergent sequences, indicating that they represent new subfamilies of these transposons. These results indicate that this family of transposons is extremely common in animal genomes, with multiple representatives in most genomes.