The good, the bad and the ugly of transposable elements annotation tools.

Transposable elements are repetitive and mobile DNA segments that can be found in virtually all organisms investigated to date. Their complex structure and variable nature are particularly challenging from the genomic annotation point of view. Many softwares have been developed to automate and facilitate TEs annotation at the genomic level, but they are highly heterogeneous regarding documentation, usability and methods. In this review, we revisited the existing software for TE genomic annotation, concentrating on the most often used ones, the methodologies they apply, and usability. Building on the state of the art of TE annotation software we propose best practices and highlight the strengths and weaknesses from the available solutions.

Systematics and spatio-temporal evolutionary patterns of the flavopilosa group of Drosophila (Diptera, Drosophilidae).

The Drosophila flavopilosa group comprises morphologically cryptic species that are ecologically restricted to feeding, breeding and ovipositing on flowers of Cestrum and Sessea (Solanaceae). Previous studies confirmed the monophyly of the group and the success of DNA barcoding in identifying a subset of its species, but several others remain yet to be evaluated. Furthemore, the taxonomy of the group remains incomplete, with only nine of the 17 species assigned to subgroups. Here, we accessed the phylogenetic relationships and spatio-temporal evolutionary patterns of the flavopilosa group based on a mitochondrial and two nuclear genes, providing the first molecular support to the subdivision of the group and suggesting a new taxonomic scheme for its species. Barcoding proved to be an effective tool, as all species were reciprocally monophyletic and different analyses of species delimitation yielded congruent results. The close relationship of D. flavopilosa with D. cestri and D. cordeiroi was strongly supported, suggesting that the latter should be placed in the flavopilosa subgroup together with the first. Furthermore, D. mariaehelenae was positioned as sister to D. incompta, supporting its inclusion in the nesiota subgroup. Despite new taxonomic assignments, the synapomorphic status of the diagnostic characters proposed for both subgroups was supported. Based on them, each of the remaining species were placed into one of both subgroups. Divergence time estimates suggest that their diversification coincided with the divergence of Sessea and Cestrum, providing an interesting case of coevolution.

Horizontal Transposon Transfer and Their Ecological Drivers: The Case of Flower-breeding Drosophila.

Understanding the mechanisms that shape the architecture, diversity, and adaptations of genomes and their ecological and genetic interfaces is of utmost importance to understand biological evolution. Transposable elements (TEs) play an important role in genome evolution, due to their ability to transpose within and between genomes, providing sites of nonallelic recombination. Here we investigate patterns and processes of TE-driven genome evolution associated with niche diversification. Specifically, we compared TE content, TE landscapes, and frequency of horizontal transposon transfers (HTTs) across genomes of flower-breeding Drosophila (FBD) with different levels of specialization on flowers. Further, we investigated whether niche breadth and ecological and geographical overlaps are associated with a potential for HTT rates. Landscape analysis evidenced a general phylogenetic pattern, in which species of the D. bromeliae group presented L-shaped curves, indicating recent transposition bursts, whereas D. lutzii showed a bimodal pattern. The great frequency of highly similar sequences recovered for all FBD suggests that these species probably experienced similar ecological pressures and evolutionary histories that contributed to the diversification of their mobilomes. Likewise, the richness of TEs superfamilies also appears to be associated with ecological traits. Furthermore, the two more widespread species, the specialist D. incompta and the generalist D. lutzii, presented the highest frequency of HTT events. Our analyses also revealed that HTT opportunities are positively influenced by abiotic niche overlap but are not associated with phylogenetic relationships or niche breadth. This suggests the existence of intermediate vectors promoting HTTs between species that do not necessarily present overlapping biotic niches.

Multiple long-range host shifts of major Wolbachia supergroups infecting arthropods.

Wolbachia is a genus of intracellular bacterial endosymbionts found in 20-66% of all insect species and a range of other invertebrates. It is classified as a single species, Wolbachia pipientis, divided into supergroups A to U, with supergroups A and B infecting arthropods exclusively. Wolbachia is transmitted mainly via vertical transmission through female oocytes, but can also be transmitted across different taxa by host shift (HS): the direct transmission of Wolbachia cells between organisms without involving vertically transmitted gametic cells. To assess the HS contribution, we recovered 50 orthologous genes from over 1000 Wolbachia genomes, reconstructed their phylogeny and calculated gene similarity. Of 15 supergroup A Wolbachia lineages, 10 have similarities ranging from 95 to 99.9%, while their hosts' similarities are around 60 to 80%. For supergroup B, four out of eight lineages, which infect diverse and distantly-related organisms such as Acari, Hemiptera and Diptera, showed similarities from 93 to 97%. These results show that Wolbachia genomes have a much higher similarity when compared to their hosts' genes, which is a major indicator of HS. Our comparative genomic analysis suggests that, at least for supergroups A and B, HS is more frequent than expected, occurring even between distantly-related species.

Oxidative and radiation stress induces transposable element transcription in Drosophila melanogaster.

It has been shown that stressors are capable of activating transposable elements (TEs). Currently, there is a hypothesis that stress activation of TEs may be involved in adaptive evolution, favouring the increase in genetic variability when the population is under adverse conditions. However, TE activation under stress is still poorly understood. In the present study, we estimated the fraction of differentially expressed TEs (DETEs) under ionizing radiation (144, 360 and 864 Gy) and oxidative stress (dioxin, formaldehyde and toluene) treatments. The stress intensity of each treatment was estimated by measuring the number of differentially expressed genes, and we show that several TEs families are activated by stress whereas others are repressed. The proportion of DETEs was positively related to stress intensity. However, even under the strongest stress, only a small fraction of TE families were activated (9.28%) and 17.72% were repressed. Considering all treatments together, the activated proportion was 19.83%. Nevertheless, as several TEs are incomplete or degenerated, only 10.55% of D. melanogaster mobilome is, at same time, activated by the stressors and able to transpose or at least code a protein. Thus, our study points out that although stress activates TEs, it is not a generalized activation process, and for some families, the stress induces repression.

Transcriptomic and Proteomic Tools in the Study of Hg Toxicity: What Is Missing?

Mercury is a hazardous substance that has unique neurodevelopmental toxic effects in humans. However, the precise sequence of molecular events that culminate in Hg-induced neuropathology is still unknown. Though the omics studies have been generating an enormous amount of new data about Hg toxicity, our ability to interpret such a large quantity of information is still limited. In this opinion article, we will reinforce the necessity of new high throughput and accurate analytical proteomic methodologies, especially, thiol and selenol-proteome. Overall, we posit that improvements in thiol- and selenol-proteomic analyses will be pivotal in identifying the primary cellular targets of Hg. However, a better understanding of the complex cascades and molecular pathways involved in its toxicity will require extensive complementary studies in more complex systems.

Cyclophosphamide in Drosophila promotes genes and transposable elements differential expression and mitochondrial dysfunction.

Cyclophosphamide (CPA) is an alkylating agent used for cancer chemotherapy, organ transplantation, and autoimmune disease treatment. Here, mRNA sequencing and high-resolution respirometry were performed to evaluate the alterations of Drosophila melanogaster gene expression fed with CPA under acute (0.1 mg/mL, for 24 h) and chronic (0.05 mg/mL, for 35 days) treatments. Differential expression analysis was performed using Cufflinks-Cuffdiff, DESeq2, and edgeR software. CPA affected genes are involved in several biological functions, including stress response and immune-related pathways, oxi-reduction and apoptotic processes, and cuticle and vitelline membrane formation. In particular, this is the first report of CPA-induced mitochondrial dysfunction caused by the downregulation of genes involved with mitochondria constituents. CPA treatment also changed the transcription pattern of transposable elements (TEs) from the gypsy and copia superfamilies. The results presented here provided evidence of CPA mitochondrial toxicity mechanisms and that CPA can modify TEs transcription in Drosophila flies.

The Catenulida flatworm can express genes from its microbiome or from the DNA it ingests.

Stenostomum are tiny planarians of the phylum Platyhelminthes that reproduce asexually. We transfected these worms using plasmids containing a gfp reporter gene. Here we show that they can express genes present in plasmids carried by bacteria and those that are encoded by naked DNA, such as plasmids or PCR fragments, transfected by electroporation; they can also express genes taken up during feeding. The microbiome associated with worm maintenance was evaluated, and the results indicated that when a plasmid is maintained in the microbiome, gfp gene expression is stable. When genes originate from naked DNA or bacteria not maintained in the microbiome, GFP expression is transient. Therefore, changes in the microbiome can modify the ability of worms to express foreign genes. In stable GFP-expressing worms, NSG showed that the gfp gene was maintained in the plasmid and was not integrated into the chromosome. These results suggest that, at least for some organisms such as flatworms, the expression of genes provided by the microbiome or the environment can be considered among the potential sources of phenotypic plasticity, which can have implications for evolvability.

The mobilome of Drosophila incompta, a flower-breeding species: comparison of transposable element landscapes among generalist and specialist flies.

The Drosophila genus is one of the main model organisms in evolutionary studies, including those investigating the role of transposable elements (TE) in genomic evolution both at the nucleotide and chromosome levels. D. incompta is a species with restricted ecology, using Cestrum (Solanaceae) flowers as unique sources for oviposition, feeding and development. In the present study, we deeply characterise the D. incompta mobilome and generate a curated dataset. A total of 277 elements were identified, corresponding to approximately 14% of the genome, and 164 of these elements are new, of which 32.62% are putatively autonomous and 8.9% are transcriptionally active in adult flies. The restricted ecology does not seem to influence the dynamics of TE in this fly, since the proportion and diversity of TEs in its genome are similar to that of other Drosophila species. This result is reinforced by the absence of a clear pattern when comparing the TE landscape between generalist and specialist flies. Using 32 available Drosophila genomes-24 ecologically generalist species and 8 specialist species-no difference was found between their TE landscape patterns. However, differences were found between species of the Sophophora and Drosophila subgenus, indicating there are lineage-specific factors shaping TE landscapes.

A re-annotation of the Anopheles darlingi mobilome.

The mobilome, portion of the genome composed of transposable elements (TEs), of Anopheles darlingi was described together with the genome of this species. Here, this mobilome was revised using similarity and de novo search approaches. A total of 5.6% of the A. darlingi genome is derived of TEs. Class I gypsy and copia were the most abundant superfamilies, corresponding to 22.36% of the mobilome. Non-LTR elements of the R1 and Jockey superfamilies account for 11% of the TEs. Among Class II TEs, the mariner superfamily is the most abundant (16.01%). Approximately 87% of the A. darlingi mobilome consist of short, truncated and/or degenerated copies of TEs. Only three retrotransposons, two belonging to gypsy and one to copia superfamilies, are putatively active elements. Only one Class II element, belonging to the mariner superfamily, is putatively active, having 12 copies in the genome. The TE landscape of A. darlingi is formed primarily by degenerated elements and, therefore, somewhat stable. Future applications of TE-based vectors for genetic transformation of A. darlingi should take into consideration mariner and piggyBac transposons, because full length and putatively active copies of these elements are present in its genome.

The somatic mobilization of transposable element mariner-Mos1 during the Drosophila lifespan and its biological consequences.

Transposable elements (TEs) are mobile DNA sequences on genomes. Some elements are able to transpose in somatic cells, a process known as somatic transposition (ST), which has been associated with detrimental biological effects. The mariner-Mos1 element of Drosophila promotes transposition in somatic and germline cells and is an excellent model for studies related to the biological consequence of somatic excision (SE). In this work, we used temperature stress to induce increasing transposition of mariner-Mos1 during different stages of the development of D. simulans, aiming to quantify SE during lifespan. Furthermore, strains of D. melanogaster exhibiting differential expression of mariner-Mos1 were employed for estimating some biological consequences of mariner mobilization. It is shown that SE of mariner-Mos1 was not constant during development; the larval phase had the highest rates while the pupal stage exhibited lower rates, and in the embryonic stage, no difference was detected. SE can be detrimental, as suggested by correlation in SE level and reduction in behavioral activities and embryonic viability. This study showed that mariner-Mos1 SE accumulates during the Drosophila life cycle, and can be involved in detrimental effects.

Drosophila relics hobo and hobo-MITEs transposons as raw material for new regulatory networks.

Hypermutable strains of Drosophila simulans have been studied for 20 years. Several mutants were isolated and characterized, some of which had phenotypes associated with alteration in development; for example, showing ectopic legs with eyes being expressed in place of antennae. The causal agent of this hypermutability is a non-autonomous hobo-related sequence (hoboVA). Around 100 mobilizable copies of this element are present in the D. simulans genome, and these are likely mobilized by the autonomous and canonical hobo element. We have shown that hoboVA has transcription factor binding sites for the developmental genes, hunchback and even-skipped, and that this transposon is expressed in embryos, following the patterns of these genes. We suggest that hobo and hobo-related elements can be material for the emergence of new regulatory networks.

A Mos1 transposase in vivo assay to screen new HIV-1 integrase inhibitors.

The integrase and transposase enzymes of retrovirus and transposons, respectively, share the catalytic DDE domain. In vitro assays showed that inhibitors of HIV-1 integrase generally inhibit the mariner Mos1 transposase. Using a Drosophila strain in which the mobilisation of the mariner element can be quantified by mosaic eyes, we showed that flies maintained in medium containing 210 µM to 4 mM of raltegravir, or 1 or 2 mM of dolutegravir, which are HIV-1 integrase inhibitor used in AIDS treatment, have 23-33% less somatic mobilisation in mosaic eyes when treated with raltegravir and 28-32% when treated with dolutegravir. The gene expression of the mariner transposase gene, estimated by qPCR, is similar among treated and control flies. The results suggest that in vivo assays using Drosophila can be used as a primary screening of inhibitory drugs for transposase and retroviral integrase. The advantages of this assay are that it is easy, quick, cheap and is an in vivo test, meaning that the tested substance has to have been taken in by cells and has arrived at the target site, which is not the case when in vitro assays are applied.

Drosophila relics hobo and hobo-MITEs transposons as raw material for new regulatory networks

Abstract Hypermutable strains of Drosophila simulans have been studied for 20 years. Several mutants were isolated and characterized, some of which had phenotypes associated with alteration in development; for example, showing ectopic legs with eyes being expressed in place of antennae. The causal agent of this hypermutability is a non-autonomous hobo-related sequence (hoboVA). Around 100 mobilizable copies of this element are present in the D. simulans genome, and these are likely mobilized by the autonomous and canonical hobo element. We have shown that hoboVA has transcription factor binding sites for the developmental genes, hunchback and even-skipped, and that this transposon is expressed in embryos, following the patterns of these genes. We suggest that hobo and hobo-related elements can be material for the emergence of new regulatory networks.

Somatizing the transposons action.

The somatic mobilization of transposable elements is more common than previously thought. In this review we discuss how the intensity and the biologic consequences of somatic mobilization are dependent on the transposable elements landscapes of each genome, and on the "momentum" of each particular TE with respect to the mechanisms that control its transposition and the possibility to escape this control. Additionally, the biologic consequences of somatic mobilization vary among organisms that show an early separation between the germline and somatic cells and those organisms that do not exhibit this separation or that reproduce asexually. In the former, somatic transposition can be involved in phenotypic plasticity, detrimental conditions such as disease, or processes such as aging. For the organisms without separation between the germ and soma, somatic mobilization can be a source of genetic variability.

The rearranged mitochondrial genome of Leptopilina boulardi (Hymenoptera: Figitidae), a parasitoid wasp of Drosophila.

The partial mitochondrial genome sequence of Leptopilina boulardi (Hymenoptera: Figitidae) was characterized. Illumina sequencing was used yielding 35,999,679 reads, from which 102,482 were utilized in the assembly. The length of the sequenced region of this partial mitochondrial genome is 15,417 bp, consisting of 13 protein-coding, two rRNA, and 21tRNA genes (the trnaM failed to be sequenced) and a partial A+T-rich region. All protein-coding genes start with ATN codons. Eleven protein-coding genes presented TAA stop codons, whereas ND6 and COII that presented TA, and T nucleotides, respectively. The gene pattern revealed extensive rearrangements compared to the typical pattern generally observed in insects. These rearrangements involve two protein-coding and two ribosomal genes, along with the 16 tRNA genes. This gene order is different from the pattern described for Ibalia leucospoides (Ibaliidae, Cynipoidea), suggesting that this particular gene order can be variable among Cynipoidea superfamily members. A maximum likelihood phylogenetic analysis of the main groups of Apocrita was performed using amino acid sequence of 13 protein-coding genes, showing monophyly for the Cynipoidea superfamily within the Hymenoptera phylogeny.

Unraveling the evolutionary scenario of the hobo element in populations of Drosophila melanogaster and D. simulans in South America using the TPE repeats as markers.

Transposable elements (TEs) are nucleotide sequences found in most studied genomes. These elements are highly diversified and have a large variation in nucleotide structure and mechanisms of transposition. hobo is a member of class II, belonging to hAT superfamily, described inDrosophila melanogaster, and it presents in its Open Reading Frame, a repetitive region encoding the amino acids threonine-proline-glutamic acid (TPE), which shows variability in the number of repeats in some regions of the world. Due to this variability some evolutionary scenarios of the hobo element are discussed, such as the scenario of the invasion of hobo element in populations ofD. melanogaster. In the present study, we investigated 22 DNA sequences of D. melanogaster and seven sequences ofD. simulans, both from South America, to check the number of repetitions of TPE, in order to clarify the evolutionary scenario of thehobo element in these populations. Our results showed a monomorphism in populations of both species in South America, with only three TPE repeats. Hence, we discuss and propose an evolutionary scenario of the invasion of the hobo element in populations of D. melanogaster and D. simulans.

Molecular and sensory mechanisms to mitigate sunlight-induced DNA damage in treefrog tadpoles.

The increased incidence of solar ultraviolet B (UVB) radiation has been proposed as an environmental stressor, which may help to explain the enigmatic decline of amphibian populations worldwide. Despite growing knowledge regarding the UV-induced biological effects in several amphibian models, little is known about the efficacy of DNA repair pathways. In addition, little attention has been given to the interplay between these molecular mechanisms with other physiological strategies that avoid the damage induced by sunlight. Here, DNA lesions induced by environmental doses of solar UVB and UVA radiation were detected in genomic DNA samples of treefrog tadpoles (Hypsiboas pulchellus) and their DNA repair activity was evaluated. These data were complemented by monitoring the induction of apoptosis in blood cells and tadpole survival. Furthermore, the tadpoles' ability to perceive and escape from UV wavelengths was evaluated as an additional strategy of photoprotection. The results show that tadpoles are very sensitive to UVB light, which could be explained by the slow DNA repair rates for both cyclobutane pyrimidine dimers (CPDs) and pyrimidine (6,4) pyrimidone photoproducts (6,4PPs). However, they were resistant to UVA, probably as a result of the activation of photolyases during UVA irradiation. Surprisingly, a sensory mechanism that triggers their escape from UVB and UVA light avoids the generation of DNA damage and helps to maintain the genomic integrity. This work demonstrates the genotoxic impact of both UVB and UVA radiation on tadpoles and emphasizes the importance of the interplay between molecular and sensory mechanisms to minimize the damage caused by sunlight.

BuT2 is a member of the third major group of hAT transposons and is involved in horizontal transfer events in the genus Drosophila.

The hAT superfamily comprises a large and diverse array of DNA transposons found in all supergroups of eukaryotes. Here we characterized the Drosophila buzzatii BuT2 element and found that it harbors a five-exon gene encoding a 643-aa putatively functional transposase. A phylogeny built with 85 hAT transposases yielded, in addition to the two major groups already described, Ac and Buster, a third one comprising 20 sequences that includes BuT2, Tip100, hAT-4_BM, and RP-hAT1. This third group is here named Tip. In addition, we studied the phylogenetic distribution and evolution of BuT2 by in silico searches and molecular approaches. Our data revealed BuT2 was, most often, vertically transmitted during the evolution of genus Drosophila being lost independently in several species. Nevertheless, we propose the occurrence of three horizontal transfer events to explain its distribution and conservation among species. Another aspect of BuT2 evolution and life cycle is the presence of short related sequences, which contain similar 5' and 3' regions, including the terminal inverted repeats. These sequences that can be considered as miniature inverted repeat transposable elements probably originated by internal deletion of complete copies and show evidences of recent mobilization.

hobo-brothers elements and their time and place for horizontal transfer.

Transposable elements (TEs) are ubiquitous components of nearly all genomes studied. These elements are highly variable in copy number, molecular structure and transposition strategies. They can move within and between genomes, thus increasing their copy numbers and avoiding being eliminated by stochastic and deterministic processes. hobo is a class II element isolated from Drosophila melanogaster. Previous phylogenetic analyses have shown that the canonical hobo element from D. melanogaster has a sister group formed by sequences found in D. willistoni (called howilli2) and D. mojavensis (called homo1). In the present study, we investigated 36 Drosophilidae species for sequences similar to howilli2 and homo1 using degenerate primers. Additionally, in silico searches were performed in 21 available Drosophila genomes. The obtained sequences formed a monophyletic sister group with the canonical hobo element; we termed these sequences 'hobo-brothers' elements. These elements showed a patch distribution and incongruities with the TE and host species phylogenies, suggesting possible cases of horizontal transfer (HT). Species that possess hobo-brothers sequences are from the New World, mainly Neotropical areas. In addition, the estimated divergence of the sequences found showed that these elements are or were recently active; the large number of HT events observed suggests that these elements could be experiencing an expansion process in Neotropical genomes. A comparison of these results with the literature is discussed with regard to the importance of the time and location of horizontal transposon transfer events.