ChimeraTE: a pipeline to detect chimeric transcripts derived from genes and transposable elements.

Transposable elements (TEs) produce structural variants and are considered an important source of genetic diversity. Notably, TE-gene fusion transcripts, i.e. chimeric transcripts, have been associated with adaptation in several species. However, the identification of these chimeras remains hindered due to the lack of detection tools at a transcriptome-wide scale, and to the reliance on a reference genome, even though different individuals/cells/strains have different TE insertions. Therefore, we developed ChimeraTE, a pipeline that uses paired-end RNA-seq reads to identify chimeric transcripts through two different modes. Mode 1 is the reference-guided approach that employs canonical genome alignment, and Mode 2 identifies chimeras derived from fixed or insertionally polymorphic TEs without any reference genome. We have validated both modes using RNA-seq data from four Drosophila melanogaster wild-type strains. We found ∼1.12% of all genes generating chimeric transcripts, most of them from TE-exonized sequences. Approximately ∼23% of all detected chimeras were absent from the reference genome, indicating that TEs belonging to chimeric transcripts may be recent, polymorphic insertions. ChimeraTE is the first pipeline able to automatically uncover chimeric transcripts without a reference genome, consisting of two running Modes that can be used as a tool to investigate the contribution of TEs to transcriptome plasticity.

A machine learning based framework to identify and classify long terminal repeat retrotransposons.

Transposable elements (TEs) are repetitive nucleotide sequences that make up a large portion of eukaryotic genomes. They can move and duplicate within a genome, increasing genome size and contributing to genetic diversity within and across species. Accurate identification and classification of TEs present in a genome is an important step towards understanding their effects on genes and their role in genome evolution. We introduce TE-Learner, a framework based on machine learning that automatically identifies TEs in a given genome and assigns a classification to them. We present an implementation of our framework towards LTR retrotransposons, a particular type of TEs characterized by having long terminal repeats (LTRs) at their boundaries. We evaluate the predictive performance of our framework on the well-annotated genomes of Drosophila melanogaster and Arabidopsis thaliana and we compare our results for three LTR retrotransposon superfamilies with the results of three widely used methods for TE identification or classification: RepeatMasker, Censor and LtrDigest. In contrast to these methods, TE-Learner is the first to incorporate machine learning techniques, outperforming these methods in terms of predictive performance, while able to learn models and make predictions efficiently. Moreover, we show that our method was able to identify TEs that none of the above method could find, and we investigated TE-Learner's predictions which did not correspond to an official annotation. It turns out that many of these predictions are in fact strongly homologous to a known TE.

Insights into the evolutionary dynamics of Neotropical biomes from the phylogeography and paleodistribution modeling of Bromelia balansae.

PREMISE OF THE STUDY: Historical abiotic and biotic factors have strongly affected species diversification and speciation. Although pre-Pleistocene events have been linked to the divergence of several Neotropical organisms, studies have highlighted a more prominent role of Pleistocene climatic oscillations in shaping current patterns of genetic variation of plants. METHODS: We performed phylogeographic analyses based on plastidial markers and modeled the current distribution and paleodistribution of Bromelia balansae (Bromeliaceae), an herbaceous species with a wide geographical distribution in South America, to infer the processes underlying its evolutionary history. KEY RESULTS: Combined molecular and paleodistributional modeling analyses indicated retraction during the Last Glacial Maximum followed by interglacial expansion. Populations occurring in the semideciduous Atlantic Forest and the Cerrado formed two distinct genetic clusters, which have been historically or ecologically isolated since late Pliocene to early Pleistocene. Populations located in the transition zone had higher levels of genetic diversity, as expected by the long-term climatic stability in the region detected in our ecological niche models. CONCLUSIONS: Our study adds important information on how herbaceous species have been affected by past climate in Central and Southeast Brazil, helping to disentangle the complex processes that have triggered the evolution of Neotropical biota.

Genome of Rhodnius prolixus, an insect vector of Chagas disease, reveals unique adaptations to hematophagy and parasite infection.

Rhodnius prolixus not only has served as a model organism for the study of insect physiology, but also is a major vector of Chagas disease, an illness that affects approximately seven million people worldwide. We sequenced the genome of R. prolixus, generated assembled sequences covering 95% of the genome (∼ 702 Mb), including 15,456 putative protein-coding genes, and completed comprehensive genomic analyses of this obligate blood-feeding insect. Although immune-deficiency (IMD)-mediated immune responses were observed, R. prolixus putatively lacks key components of the IMD pathway, suggesting a reorganization of the canonical immune signaling network. Although both Toll and IMD effectors controlled intestinal microbiota, neither affected Trypanosoma cruzi, the causal agent of Chagas disease, implying the existence of evasion or tolerance mechanisms. R. prolixus has experienced an extensive loss of selenoprotein genes, with its repertoire reduced to only two proteins, one of which is a selenocysteine-based glutathione peroxidase, the first found in insects. The genome contained actively transcribed, horizontally transferred genes from Wolbachia sp., which showed evidence of codon use evolution toward the insect use pattern. Comparative protein analyses revealed many lineage-specific expansions and putative gene absences in R. prolixus, including tandem expansions of genes related to chemoreception, feeding, and digestion that possibly contributed to the evolution of a blood-feeding lifestyle. The genome assembly and these associated analyses provide critical information on the physiology and evolution of this important vector species and should be instrumental for the development of innovative disease control methods.

Genome-wide analysis of transposable elements in the coffee berry borer Hypothenemus hampei (Coleoptera: Curculionidae): description of novel families.

The coffee berry borer (CBB) Hypothenemus hampei is the most limiting pest of coffee production worldwide. The CBB genome has been recently sequenced; however, information regarding the presence and characteristics of transposable elements (TEs) was not provided. Using systematic searching strategies based on both de novo and homology-based approaches, we present a library of TEs from the draft genome of CBB sequenced by the Colombian Coffee Growers Federation. The library consists of 880 sequences classified as 66% Class I (LTRs: 46%, non-LTRs: 20%) and 34% Class II (DNA transposons: 8%, Helitrons: 16% and MITEs: 10%) elements, including families of the three main LTR (Gypsy, Bel-Pao and Copia) and non-LTR (CR1, Daphne, I/Nimb, Jockey, Kiri, R1, R2 and R4) clades and DNA superfamilies (Tc1-mariner, hAT, Merlin, P, PIF-Harbinger, PiggyBac and Helitron). We propose the existence of novel families: Hypo, belonging to the LTR Gypsy superfamily; Hamp, belonging to non-LTRs; and rosa, belonging to Class II or DNA transposons. Although the rosa clade has been previously described, it was considered to be a basal subfamily of the mariner family. Based on our phylogenetic analysis, including Tc1, mariner, pogo, rosa and Lsra elements from other insects, we propose that rosa and Lsra elements are subfamilies of an independent family of Class II elements termed rosa. The annotations obtained indicate that a low percentage of the assembled CBB genome (approximately 8.2%) consists of TEs. Although these TEs display high diversity, most sequences are degenerate, with few full-length copies of LTR and DNA transposons and several complete and putatively active copies of non-LTR elements. MITEs constitute approximately 50% of the total TEs content, with a high proportion associated with DNA transposons in the Tc1-mariner superfamily.

Transposable elements in the Anopheles funestus transcriptome.

Transposable elements (TEs) are present in most of the eukaryotic genomes and their impact on genome evolution is increasingly recognized. Although there is extensive information on the TEs present in several eukaryotic genomes, less is known about the expression of these elements at the transcriptome level. Here we present a detailed analysis regarding the expression of TEs in Anopheles funestus, the second most important vector of human malaria in Africa. Several transcriptionally active TE families belonging both to Class I and II were identified and characterized. Interestingly, we have identified a full-length putative active element (including the presence of full length TIRs in the genomic sequence) belonging to the hAT superfamily, which presents active members in other insect genomes. This work contributes to a comprehensive understanding of the landscape of transposable elements in A. funestus transcriptome. Our results reveal that TEs are abundant and diverse in the mosquito and that most of the TE families found in the genome are represented in the mosquito transcriptome, a fact that could indicate activity of these elements.The vast diversity of TEs expressed in A. funestus suggests that there is ongoing amplification of several families in this organism.

Learning HMMs for nucleotide sequences from amino acid alignments.

Profile hidden Markov models (profile HMMs) are known to efficiently predict whether an amino acid (AA) sequence belongs to a specific protein family. Profile HMMs can also be used to search for protein domains in genome sequences. In this case, HMMs are typically learned from AA sequences and then used to search on the six-frame translation of nucleotide (NT) sequences. However, this approach demands additional processing of the original data and search results. Here, we propose an alternative and more direct method which converts an AA alignment into an NT one, after which an NT-based HMM is trained to be applied directly on a genome.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH).

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Analysis of HCV quasispecies dynamic under selective pressure of combined therapy.

BACKGROUND: The quasispecies composition of Hepatitis C virus (HCV) could have important implications with regard to viral persistence and response to interferon-based therapy. The complete NS5A was analyzed to evaluate whether the composition of NS5A quasispecies of HCV 1a/1b is related to responsiveness to combined interferon pegylated (PEG-IFN) and ribavirin therapy. METHODS: Viral RNA was isolated from serum samples collected before, during and after treatment from virological sustained responder (SVR), non-responder (NR) and the end-of-treatment responder patients (ETR). NS5A region was amplified, cloned and sequenced. Six hundred and ninety full-length NS5A sequences were analyzed. RESULTS: This study provides evidence that lower nucleotide diversity of the NS5A region pre-therapy is associated with viral clearance. Analysis of samples of NRs and the ETRs time points showed that genetic diversity of populations tend to decrease over time. Post-therapy population of ETRs presented higher genetic distance from baseline probably due to the bottleneck phenomenon observed for those patients in the end of treatment. The viral effective population of those patients also showed a strong decrease after therapy. Otherwise, NRs demonstrated a continuous variation or stability of effective populations and genetic diversity over time that did not seem to be related to therapy. Phylogenetic relationships concerning complete NS5A sequences obtained from patients did not demonstrate clustering associated with specific response patterns. However, distinctive clustering of pre/post-therapy sequences was observed. In addition, the evolution of quasispecies over time was subjected to purifying or relaxed purifying selection. Codons 157 (P03), 182 and 440 (P42), 62 and 404 (P44) were found to be under positive selective pressure but it failed to be related to the therapy. CONCLUSION: These results confirm the hypothesis that a relationship exists between NS5A heterogeneity and response to therapy in patients infected with chronic hepatitis C.

Losing identity: structural diversity of transposable elements belonging to different classes in the genome of Anopheles gambiae.

BACKGROUND: Transposable elements (TEs), both DNA transposons and retrotransposons, are genetic elements with the main characteristic of being able to mobilize and amplify their own representation within genomes, utilizing different mechanisms of transposition. An almost universal feature of TEs in eukaryotic genomes is their inability to transpose by themselves, mainly as the result of sequence degeneration (by either mutations or deletions). Most of the elements are thus either inactive or non-autonomous. Considering that the bulk of some eukaryotic genomes derive from TEs, they have been conceived as "TE graveyards." It has been shown that once an element has been inactivated, it progressively accumulates mutations and deletions at neutral rates until completely losing its identity or being lost from the host genome; however, it has also been shown that these "neutral sequences" might serve as raw material for domestication by host genomes. RESULTS: We have analyzed the sequence structural variations, nucleotide divergence, and pattern of insertions and deletions of several superfamilies of TEs belonging to both class I (long terminal repeats [LTRs] and non-LTRs [NLTRs]) and II in the genome of Anopheles gambiae, aiming at describing the landscape of deterioration of these elements in this particular genome. Our results describe a great diversity in patterns of deterioration, indicating lineage-specific differences including the presence of Solo-LTRs in the LTR lineage, 5'-deleted NLTRs, and several non-autonomous and MITEs in the class II families. Interestingly, we found fragments of NLTRs corresponding to the RT domain, which preserves high identity among them, suggesting a possible remaining genomic role for these domains. CONCLUSIONS: We show here that the TEs in the An. gambiae genome deteriorate in different ways according to the class to which they belong. This diversity certainly has implications not only at the host genomic level but also at the amplification dynamic and evolution of the TE families themselves.

Vertical inheritance and bursts of transposition have shaped the evolution of the BS non-LTR retrotransposon in Drosophila.

The history of transposable elements over evolutionary time can often be partially reconstructed on the basis of genome analysis. In this study, we identified and extensively characterized the NLTR BS retrotransposon in 12 sequenced Drosophila genomes, by its sequence diversity within and among genomes, its degeneration pattern and its transcriptional activity. We show that the BS element has a variable copy number and patchy distribution within the Drosophila genus, that it is at distinct stages of the evolutionary cycle in the different Drosophila species and that its evolution is characterized by vertical transmission and by bursts of transposition in certain species.

Comparative transcriptomics between Drosophila mojavensis and D. arizonae reveals transgressive gene expression and underexpression of spermatogenesis-related genes in hybrid testes.

Interspecific hybridization is a stressful condition that can lead to sterility and/or inviability through improper gene regulation in Drosophila species with a high divergence time. However, the extent of these abnormalities in hybrids of recently diverging species is not well known. Some studies have shown that in Drosophila, the mechanisms of postzygotic isolation may evolve more rapidly in males than in females and that the degree of viability and sterility is associated with the genetic distance between species. Here, we used transcriptomic comparisons between two Drosophila mojavensis subspecies and D. arizonae (repleta group, Drosophila) and identified greater differential gene expression in testes than in ovaries. We tested the hypothesis that the severity of the interspecies hybrid phenotype is associated with the degree of gene misregulation. We showed limited gene misregulation in fertile females and an increase in the amount of misregulation in males with more severe sterile phenotypes (motile vs. amotile sperm). In addition, for these hybrids, we identified candidate genes that were mostly associated with spermatogenesis dysfunction.

Plasmodium vivax circumsporozoite genotypes: a limited variation or new subspecies with major biological consequences?

BACKGROUND: Plasmodium vivax circumsporozoite variants have been identified in several geographical areas. The real implication of the genetic variation in this region of the P. vivax genome has been questioned for a long time. Although previous studies have observed significant association between VK210 and the Duffy blood group, we present here that evidences of this variation are limited to the CSP central portion. METHODS: The phylogenetic analyses were accomplished starting from the amplification of conserved domains of 18 SSU RNAr and Cyt B. The antibodies responses against the CSP peptides, MSP-1, AMA-1 and DBP were detected by ELISA, in plasma samples of individuals infected with two P. vivax CS genotypes: VK210 and P. vivax-like. RESULTS: These analyses of the two markers demonstrate high similarity among the P. vivax CS genotypes and surprisingly showed diversity equal to zero between VK210 and P. vivax-like, positioning these CS genotypes in the same clade. A high frequency IgG antibody against the N- and C-terminal regions of the P. vivax CSP was found as compared to the immune response to the R- and V- repetitive regions (p = 0.0005, Fisher's Exact test). This difference was more pronounced when the P. vivax-like variant was present in the infection (p = 0.003, Fisher's Exact test). A high frequency of antibody response against MSP-1 and AMA-1 peptides was observed for all P. vivax CS genotypes in comparison to the same frequency for DBP. CONCLUSIONS: This results target that the differences among the P. vivax CS variants are restrict to the central repeated region of the protein, mostly nucleotide variation with important serological consequences.

Genetic diversity of NS5A protein from hepatitis C virus genotype 3a and its relationship to therapy response.

BACKGROUND: The quasispecies nature of HCV may have important implications for viral persistence, pathogenicity and resistance to antiviral agents. The variability of one of the viral proteins, NS5A, is believed to be related to the response to IFN therapy, the standard treatment for infection. In this study we analyzed the quasispecies composition of NS5A protein in patients infected with HCV genotype 3a, before IFN therapy. METHODS: Viral RNA was isolated from samples of 12 patients: four sustained virological responders (SVR), four non-responders (NR), and four end-of-treatment responders (ETR). cDNA was synthesized, the NS5A region was amplified and the fragments obtained were cloned. Fifteen clones from each patient were sequenced with eight primers, generating 179 contigs. RESULTS: Higher values for substitution (either synonymous or non-synonymous) and for distance were found in the SVR group. However, the NR group showed relatively more non-synonymous mutations than the other groups, owing to the higher values of dN/dS in complete NS5A and most specific regions. Overall, NS5A protein is undergoing purifying selection, since all dN/dS ratios values are below 0.5. CONCLUSIONS: Our study provides an overview of the genetic variability of complete NS5A protein in HCV genotype 3a.

Homology-Free Detection of Transposable Elements Unveils Their Dynamics in Three Ecologically Distinct Rhodnius Species.

Transposable elements (TEs) are widely distributed repetitive sequences in the genomes across the tree of life, and represent an important source of genetic variability. Their distribution among genomes is specific to each lineage. A phenomenon associated with this feature is the sudden expansion of one or several TE families, called bursts of transposition. We previously proposed that bursts of the Mariner family (DNA transposons) contributed to the speciation of Rhodnius prolixus Stål, 1859. This hypothesis motivated us to study two additional species of the R. prolixus complex: Rhodnius montenegrensis da Rosa et al., 2012 and Rhodnius marabaensis Souza et al., 2016, together with a new, de novo annotation of the R. prolixus repeatome using unassembled short reads. Our analysis reveals that the total amount of TEs present in Rhodnius genomes (19% to 23.5%) is three to four times higher than that expected based on the original quantifications performed for the original genome description of R. prolixus. We confirm here that the repeatome of the three species is dominated by Class II elements of the superfamily Tc1-Mariner, as well as members of the LINE order (Class I). In addition to R. prolixus, we also identified a recent burst of transposition of the Mariner family in R. montenegrensis and R. marabaensis, suggesting that this phenomenon may not be exclusive to R. prolixus. Rather, we hypothesize that whilst the expansion of Mariner elements may have contributed to the diversification of the R. prolixus-R. robustus species complex, the distinct ecological characteristics of these new species did not drive the general evolutionary trajectories of these TEs.

A single P115Q mutation modulates specificity in the Corynebacterium pseudotuberculosis arginine repressor.

The arginine repressor (ArgR) regulates the expression of genes involved in arginine biosynthesis. Upon attaining a threshold concentration of arginine in the cytoplasm, the trimeric C-terminal domain of ArgR binds three arginines in a shallow surface cleft and subsequently hexamerizes forming a dimer of trimers containing six Arg co-repressor molecules which are buried at the subunit interfaces. The N-terminal domains of this complex bind to the DNA promoter thereby interrupting the transcription of the genes related to Arg biosynthesis. The crystal structures of the wild type and mutant Pro115Gln ArgR from Corynebacterium pseudotuberculosis determined at 1.7 Å demonstrate that a single amino acid substitution switches co-repressor specificity from Tyr to Arg. Molecular dynamics simulations indicate that the first step, i.e., the binding of the co-repressor, occurs in the trimeric state and that Pro115Gln ArgR preferentially binds Arg. It was also shown that, in Pro115 ArgR hexamers, the concomitant binding of sodium ions shifts selectivity to Tyr. Structural data combined with phylogenetic analyses of ArgR from C. pseudotuberculosis suggest that substitutions in the binding pocket at position 115 may alter its specificity for amino acids and that the length of the protein interdomain linker can provide further functional flexibility. These results support the existence of alternative ArgR regulatory mechanisms in this pathogenic bacterium.

The evolutionary dynamics of the Helena retrotransposon revealed by sequenced Drosophila genomes.

BACKGROUND: Several studies have shown that genomes contain a mixture of transposable elements, some of which are still active and others ancient relics that have degenerated. This is true for the non-LTR retrotransposon Helena, of which only degenerate sequences have been shown to be present in some species (Drosophila melanogaster), whereas putatively active sequences are present in others (D. simulans). Combining experimental and population analyses with the sequence analysis of the 12 Drosophila genomes, we have investigated the evolution of Helena, and propose a possible scenario for the evolution of this element. RESULTS: We show that six species of Drosophila have the Helena transposable element at different stages of its evolution. The copy number is highly variable among these species, but most of them are truncated at the 5' ends and also harbor several internal deletions and insertions suggesting that they are inactive in all species, except in D. mojavensis in which quantitative RT-PCR experiments have identified a putative active copy. CONCLUSION: Our data suggest that Helena was present in the common ancestor of the Drosophila genus, which has been vertically transmitted to the derived lineages, but that it has been lost in some of them. The wide variation in copy number and sequence degeneration in the different species suggest that the evolutionary dynamics of Helena depends on the genomic environment of the host species.

Multiple invasions of Gypsy and Micropia retroelements in genus Zaprionus and melanogaster subgroup of the genus Drosophila.

BACKGROUND: The Zaprionus genus shares evolutionary features with the melanogaster subgroup, such as space and time of origin. Although little information about the transposable element content in the Zaprionus genus had been accumulated, some of their elements appear to be more closely related with those of the melanogaster subgroup, indicating that these two groups of species were involved in horizontal transfer events during their evolution. Among these elements, the Gypsy and the Micropia retroelements were chosen for screening in seven species of the two Zaprionus subgenera, Anaprionus and Zaprionus. RESULTS: Screening allowed the identification of diverse Gypsy and Micropia retroelements only in species of the Zaprionus subgenus, showing that they are transcriptionally active in the sampled species. The sequences of each retroelement were closely related to those of the melanogaster species subgroup, and the most parsimonious hypothesis would be that 15 horizontal transfer events shaped their evolution. The Gypsy retroelement of the melanogaster subgroup probably invaded the Zaprionus genomes about 11 MYA. In contrast, the Micropia retroelement may have been introduced into the Zaprionus subgenus and the melanogaster subgroup from an unknown donor more recently (~3 MYA). CONCLUSION: Gypsy and Micropia of Zaprionus and melanogaster species share similar evolutionary patterns. The sharing of evolutionary, ecological and ethological features probably allowed these species to pass through a permissive period of transposable element invasion, explaining the proposed waves of horizontal transfers.

The protist Trichomonas vaginalis harbors multiple lineages of transcriptionally active Mutator-like elements.

BACKGROUND: For three decades the Mutator system was thought to be exclusive of plants, until the first homolog representatives were characterized in fungi and in early-diverging amoebas earlier in this decade. RESULTS: Here, we describe and characterize four families of Mutator-like elements in a new eukaryotic group, the Parabasalids. These Trichomonas vaginalis Mutator- like elements, or TvMULEs, are active in T. vaginalis and patchily distributed among 12 trichomonad species and isolates. Despite their relatively distinctive amino acid composition, the inclusion of the repeats TvMULE1, TvMULE2, TvMULE3 and TvMULE4 into the Mutator superfamily is justified by sequence, structural and phylogenetic analyses. In addition, we identified three new TvMULE-related sequences in the genome sequence of Candida albicans. While TvMULE1 is a member of the MuDR clade, predominantly from plants, the other three TvMULEs, together with the C. albicans elements, represent a new and quite distinct Mutator lineage, which we named TvCaMULEs. The finding of TvMULE1 sequence inserted into other putative repeat suggests the occurrence a novel TE family not yet described. CONCLUSION: These findings expand the taxonomic distribution and the range of functional motif of MULEs among eukaryotes. The characterization of the dynamics of TvMULEs and other transposons in this organism is of particular interest because it is atypical for an asexual species to have such an extreme level of TE activity; this genetic landscape makes an interesting case study for causes and consequences of such activity. Finally, the extreme repetitiveness of the T. vaginalis genome and the remarkable degree of sequence identity within its repeat families highlights this species as an ideal system to characterize new transposable elements.

High frequency of horizontal transfer in Jockey families (LINE order) of drosophilids.

BACKGROUND: The use of large-scale genomic analyses has resulted in an improvement of transposable element sampling and a significant increase in the number of reported HTT (horizontal transfer of transposable elements) events by expanding the sampling of transposable element sequences in general and of specific families of these elements in particular, which were previously poorly sampled. In this study, we investigated the occurrence of HTT events in a group of elements that, until recently, were uncommon among the HTT records in Drosophila - the Jockey elements, members of the LINE (long interspersed nuclear element) order of non-LTR (long terminal repeat) retrotransposons. The sequences of 111 Jockey families deposited in Repbase that met the criteria of the analysis were used to identify Jockey sequences in 48 genomes of Drosophilidae (genus Drosophila, subgenus Sophophora: melanogaster, obscura and willistoni groups; subgenus Drosophila: immigrans, melanica, repleta, robusta, virilis and grimshawi groups; subgenus Dorsilopha: busckii group; genus/subgenus Zaprionus and genus Scaptodrosophila). RESULTS: Phylogenetic analyses revealed 72 Jockey families in 41 genomes. Combined analyses revealed 15 potential HTT events between species belonging to different genera and species groups of Drosophilidae, providing evidence for the flow of genetic material favoured by the spatio-temporal sharing of these species present in the Palaeartic or Afrotropical region. CONCLUSIONS: Our results provide phylogenetic, biogeographic and temporal evidence of horizontal transfers of the Jockey elements, increase the number of rare records of HTT in specific families of LINE elements, increase the number of known occurrences of these events, and enable a broad understanding of the evolutionary dynamics of these elements and the host species.