Domesticated retroviral GAG gene in Drosophila: new functions for an old gene

Association of a rapidly selected 4.3kb transposon-containing structural variation with a P450-based resistance to pyrethroids in the African malaria vector Anopheles funestus

Deciphering the evolutionary forces controlling insecticide resistance in malaria vectors remains a prerequisite to designing molecular tools to detect and assess resistance impact on control tools. Here, we demonstrate that a 4.3kb transposon-containing structural variation is associated with pyrethroid resistance in central/eastern African populations of the malaria vector Anopheles funestus. In this study, we analysed Pooled template sequencing data and direct sequencing to identify an insertion of 4.3kb containing a putative retro-transposon in the intergenic region of two P450s CYP6P5-CYP6P9b in mosquitoes of the malaria vector Anopheles funestus from Uganda. We then designed a PCR assay to track its spread temporally and regionally and decipher its role in insecticide resistance. The insertion originates in or near Uganda in East Africa, where it is fixed and has spread to high frequencies in the Central African nation of Cameroon but is still at low frequency in West Africa and absent in Southern Africa. A marked and rapid selection was observed with the 4.3kb-SV frequency increasing from 3% in 2014 to 98% in 2021 in Cameroon. A strong association was established between this SV and pyrethroid resistance in field populations and is reducing pyrethroid-only nets' efficacy. Genetic crosses and qRT-PCR revealed that this SV enhances the expression of CYP6P9a/b but not CYP6P5. Within this structural variant (SV), we identified putative binding sites for transcription factors associated with the regulation of detoxification genes. An inverse correlation was observed between the 4.3kb SV and malaria parasite infection, indicating that mosquitoes lacking the 4.3kb SV were more frequently infected compared to those possessing it. Our findings highlight the underexplored role and rapid spread of SVs in the evolution of insecticide resistance and provide additional tools for molecular surveillance of insecticide resistance.

The Lifespan of D. melanogaster Depends on the Function of the Gagr Gene, a Domesticated gag Gene of Drosophila LTR Retrotransposons

(1) Background: The Gagr gene in Drosophila melanogaster's genome originated from the molecular domestication of retrotransposons and retroviruses' gag gene. In all Drosophila species, the Gagr protein homologs exhibit a conserved structure, indicative of a vital role. Previous studies have suggested a potential link between the Gagr gene function and stress responses. (2) Methods: We compared flies with Gagr gene knockdown in all tissues to control flies in physiological tests and RNA-sequencing experiments. (3) Results: Flies with the Gagr gene knockdown exhibited shorter lifespans compared to control flies. Transcriptome analysis revealed that Gagr knockdown flies showed elevated transcription levels of immune response genes. We used ammonium persulfate, a potent stress inducer, to elicit a stress response. In control flies, ammonium persulfate activated the Toll, JAK/STAT, and JNK/MAPK signaling pathways. In contrast, flies with the Gagr gene knockdown displayed reduced expression of stress response genes. Gene ontology enrichment analysis identified categories of genes upregulated under ammonium persulfate stress in control flies but not in Gagr knockdown flies. These genes are involved in developmental control, morphogenesis, and central nervous system function. (4) Conclusion: Our findings indicate the significance of the Gagr gene in maintaining immune response and homeostasis.

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.

Domesticated LTR-Retrotransposon gag-Related Gene (Gagr) as a Member of the Stress Response Network in Drosophila

The most important sources of new components of genomes are transposable elements, which can occupy more than half of the nucleotide sequence of the genome in higher eukaryotes. Among the mobile components of a genome, a special place is occupied by retroelements, which are similar to retroviruses in terms of their mechanisms of integration into a host genome. The process of positive selection of certain sequences of transposable elements and retroviruses in a host genome is commonly called molecular domestication. There are many examples of evolutionary adaptations of gag (retroviral capsid) sequences as new regulatory sequences of different genes in mammals, where domesticated gag genes take part in placenta functioning and embryogenesis, regulation of apoptosis, hematopoiesis, and metabolism. The only gag-related gene has been found in the Drosophila genome—Gagr. According to the large-scale transcriptomic and proteomic analysis data, the Gagr gene in D. melanogaster is a component of the protein complex involved in the stress response. In this work, we consider the evolutionary processes that led to the formation of a new function of the domesticated gag gene and its adaptation to participation in the stress response. We discuss the possible functional role of the Gagr as part of the complex with its partners in Drosophila, and the pathway of evolution of proteins of the complex in eukaryotes to determine the benefit of the domesticated retroelement gag gene.

Unearthing LTR Retrotransposon gag Genes Co-opted in the Deep Evolution of Eukaryotes

LTR retrotransposons comprise a major component of the genomes of eukaryotes. On occasion, retrotransposon genes can be recruited by their hosts for diverse functions, a process formally referred to as co-option. However, a comprehensive picture of LTR retrotransposon gag gene co-option in eukaryotes is still lacking, with several documented cases exclusively involving Ty3/Gypsy retrotransposons in animals. Here, we use a phylogenomic approach to systemically unearth co-option of retrotransposon gag genes above the family level of taxonomy in 2,011 eukaryotes, namely co-option occurring during the deep evolution of eukaryotes. We identify a total of 14 independent gag gene co-option events across more than 740 eukaryote families, eight of which have not been reported previously. Among these retrotransposon gag gene co-option events, nine, four, and one involve gag genes of Ty3/Gypsy, Ty1/Copia, and Bel-Pao retrotransposons, respectively. Seven, four, and three co-option events occurred in animals, plants, and fungi, respectively. Interestingly, two co-option events took place in the early evolution of angiosperms. Both selective pressure and gene expression analyses further support that these co-opted gag genes might perform diverse cellular functions in their hosts, and several co-opted gag genes might be subject to positive selection. Taken together, our results provide a comprehensive picture of LTR retrotransposon gag gene co-option events that occurred during the deep evolution of eukaryotes and suggest paucity of LTR retrotransposon gag gene co-option during the deep evolution of eukaryotes.

Nanopore RNA Sequencing Revealed Long Non-Coding and LTR Retrotransposon-Related RNAs Expressed at Early Stages of Triticale SEED Development

The intergenic space of plant genomes encodes many functionally important yet unexplored RNAs. The genomic loci encoding these RNAs are often considered “junk”, DNA as they are frequently associated with repeat-rich regions of the genome. The latter makes the annotations of these loci and the assembly of the corresponding transcripts using short RNAseq reads particularly challenging. Here, using long-read Nanopore direct RNA sequencing, we aimed to identify these “junk” RNA molecules, including long non-coding RNAs (lncRNAs) and transposon-derived transcripts expressed during early stages (10 days post anthesis) of seed development of triticale (AABBRR, 2n = 6x = 42), an interspecific hybrid between wheat and rye. Altogether, we found 796 lncRNAs and 20 LTR retrotransposon-related transcripts (RTE-RNAs) expressed at this stage, with most of them being previously unannotated and located in the intergenic as well as intronic regions. Sequence analysis of the lncRNAs provide evidence for the frequent exonization of Class I (retrotransposons) and class II (DNA transposons) transposon sequences and suggest direct influence of “junk” DNA on the structure and origin of lncRNAs. We show that the expression patterns of lncRNAs and RTE-related transcripts have high stage specificity. In turn, almost half of the lncRNAs located in Genomes A and B have the highest expression levels at 10–30 days post anthesis in wheat. Detailed analysis of the protein-coding potential of the RTE-RNAs showed that 75% of them carry open reading frames (ORFs) for a diverse set of GAG proteins, the main component of virus-like particles of LTR retrotransposons. We further experimentally demonstrated that some RTE-RNAs originate from autonomous LTR retrotransposons with ongoing transposition activity during early stages of triticale seed development. Overall, our results provide a framework for further exploration of the newly discovered lncRNAs and RTE-RNAs in functional and genome-wide association studies in triticale and wheat. Our study also demonstrates that Nanopore direct RNA sequencing is an indispensable tool for the elucidation of lncRNA and retrotransposon transcripts.

Domesticated gag Gene of Drosophila LTR Retrotransposons Is Involved in Response to Oxidative Stress

Drosophila melanogaster is one of the most extensively used genetic model organisms for studying LTR retrotransposons that are represented by various groups in its genome. However, the phenomenon of molecular domestication of LTR retrotransposons has been insufficiently studied in Drosophila, as well as in other invertebrates. The present work is devoted to studying the role of the domesticated gag gene, Gagr, in the Drosophila genome. The Gagr gene has been shown to be involved in the response to stress caused by exposure to ammonium persulfate, but not in the stress response to oligomycin A, zeomycin, and cadmium chloride. Ammonium persulfate tissue specifically activates the expression of Gagr in the tissues of the carcass, but not in the gut. We found that the Gagr gene promoter contains one binding motif for the transcription factor kayak, a component of the JNK signaling pathway, and two binding motifs for the transcription factor Stat92E, a component of the Jak-STAT signaling pathway. Remarkably, Gagr orthologs contain the second binding motif for Stat92E only in D. melanogaster, D. simulans and D. sechellia, whereas in D. yakuba and D. erecta, Gagr orthologs contain a single motif, and there are no binding sites for Stat92E in the promoters of Gagr orthologs in D. ananassae and in species outside the melanogaster group. The data obtained indicate the formation of the protective function of the Gagr gene during evolution.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

Symbiosis: Viruses as Intimate Partners

Viruses must establish an intimate relationship with their hosts and vectors in order to infect, replicate, and disseminate; hence, viruses can be considered as symbionts with their hosts. Symbiotic relationships encompass different lifestyles, including antagonistic (or pathogenic, the most well-studied lifestyle for viruses), commensal (probably the most common lifestyle), and mutualistic (important beneficial partners). Symbiotic relationships can shape the evolution of the partners in a holobiont, and placing viruses in this context provides an important framework for understanding virus-host relationships and virus ecology. Although antagonistic relationships are thought to lead to coevolution, this is not always clear in virus-host interactions, and impacts on evolution may be complex. Commensalism implies a hitchhiking role for viruses-selfish elements just along for the ride. Mutualistic relationships have been described in detail in the past decade, and they reveal how important viruses are in considering host ecology. Ultimately, symbiosis can lead to symbiogenesis, or speciation through fusion, and the presence of large amounts of viral sequence in the genomes of everything from bacteria to humans, including some important functional genes, illustrates the significance of viral symbiogenesis in the evolution of all life on Earth.

Mechanisms of LTR-Retroelement Transposition: Lessons from Drosophila melanogaster

Long terminal repeat (LTR) retrotransposons occupy a special place among all mobile genetic element families. The structure of LTR retrotransposons that have three open reading frames is identical to DNA forms of retroviruses that are integrated into the host genome. Several lines of evidence suggest that LTR retrotransposons share a common ancestry with retroviruses and thus are highly relevant to understanding mechanisms of transposition. Drosophila melanogaster is an exceptionally convenient model for studying the mechanisms of retrotransposon movement because many such elements in its genome are transpositionally active. Moreover, two LTRretrotransposons of D. melanogaster, gypsy and ZAM, have been found to have infectious properties and have been classified as errantiviruses. Despite numerous studies focusing on retroviral integration process, there is still no clear understanding of integration specificity in a target site. Most LTR retrotransposons non-specifically integrate into a target site. Site-specificity of integration at vertebrate retroviruses is rather relative. At the same time, sequence-specific integration is the exclusive property of errantiviruses and their derivatives with two open reading frames. The possible basis for the errantivirus integration specificity is discussed in the present review.

Close ecological relationship among species facilitated horizontal transfer of retrotransposons

Background

Horizontal transfer (HT) of genetic materials is increasingly being found in both animals and plants and mainly concerns transposable elements (TEs). Many crustaceans have big genome sizes and are thus likely to harbor high TE contents. Their habitat might offer them ample opportunities to exchange genetic materials with organisms that are ecologically close but taxonomically distant to them.

Results

In this study, we analyzed the transcriptome of Pacific white shrimp (Litopenaeus vannamei), an important economic crustacean, to explore traces of HT events. From a collection of newly assembled transcripts, we identified 395 high reliable TE transcripts, most of which were retrotransposon transcripts. One hundred fifty-seven of those transcripts showed highest similarity to sequences from non-arthropod organisms, including ray-finned fishes, mollusks and putative parasites. In total, 16 already known L. vannamei TE families are likely to be involved in horizontal transfer events. Phylogenetic analyses of 10 L. vannamei TE families and their homologues (protein sequences) revealed that L. vannamei TE families were generally more close to sequences from aquatic species. Furthermore, TEs from other aquatic species also tend to group together, although they are often distantly related in taxonomy. Sequences from parasites and microorganisms were also widely present, indicating their possible important roles in HT events. Expression profile analyses of transcripts in two NCBI BioProjects revealed that transcripts involved in HT events are likely to play important roles in antiviral immunity. More specifically, those transcripts might act as inhibitors of antiviral immunity.

Conclusions

Close ecological relationship, especially predation, might greatly facilitate HT events among aquatic species. This could be achieved through exchange of parasites and microorganisms, or through direct DNA flow. The occurrence of HT events may be largely incidental, but the effects could be beneficial for recipients.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0767-0) contains supplementary material, which is available to authorized users.

Transposable Elements versus the Fungal Genome: Impact on Whole-Genome Architecture and Transcriptional Profiles

Transposable elements (TEs) are exceptional contributors to eukaryotic genome diversity. Their ubiquitous presence impacts the genomes of nearly all species and mediates genome evolution by causing mutations and chromosomal rearrangements and by modulating gene expression. We performed an exhaustive analysis of the TE content in 18 fungal genomes, including strains of the same species and species of the same genera. Our results depicted a scenario of exceptional variability, with species having 0.02 to 29.8% of their genome consisting of transposable elements. A detailed analysis performed on two strains of Pleurotus ostreatus uncovered a genome that is populated mainly by Class I elements, especially LTR-retrotransposons amplified in recent bursts from 0 to 2 million years (My) ago. The preferential accumulation of TEs in clusters led to the presence of genomic regions that lacked intra- and inter-specific conservation. In addition, we investigated the effect of TE insertions on the expression of their nearby upstream and downstream genes. Our results showed that an important number of genes under TE influence are significantly repressed, with stronger repression when genes are localized within transposon clusters. Our transcriptional analysis performed in four additional fungal models revealed that this TE-mediated silencing was present only in species with active cytosine methylation machinery. We hypothesize that this phenomenon is related to epigenetic defense mechanisms that are aimed to suppress TE expression and control their proliferation.

Transposable elements (TEs) are enigmatic genetic units that have played important roles in the evolution of eukaryotic genomes. Since their discovery in the 1950s, they have gained increasing attention and are known today as active genome modelers in multiple species. Although these elements have been widely studied in plants, much less is known about their occurrence and impact on the fungal kingdom. Using a diverse set of basidiomycete and ascomycete fungi, we quantified and characterized a huge diversity of DNA and RNA transposable elements, and we identified species that had 0.02 to 29.8% of their genomes occupied by transposable elements. In addition, using our basidiomycete model Pleurotus ostreatus, we demonstrated how TE insertions produced detrimental effects on the expression of upstream and downstream genes, which were downregulated compared with the control groups. This silencing mechanism was present in the basidiomycetes tested but exhibited a patchy distribution in ascomycetes, and might be related to specific genome defense mechanisms that control transposon proliferation. This finding reveals the broader impact of transposable elements in fungi. In addition to their importance as long-term evolutionary forces, they play major roles in the more dynamic transcriptome regulation of certain species.

In between: gypsy in Drosophila melanogaster reveals new insights into endogenous retrovirus evolution

Retroviruses are RNA viruses that are able to synthesize a DNA copy of their genome and insert it into a chromosome of the host cell. Sequencing of different eukaryote genomes has revealed the presence of many such endogenous retroviral sequences. The mechanisms by which these retroviral sequences have colonized the genome are still unknown, and the endogenous retrovirus gypsy of Drosophila melanogaster is a powerful experimental model for deciphering this process in vivo. Gypsy is expressed in a layer of somatic cells, and then transferred into the oocyte by an unknown mechanism. This critical step is the start of the endogenization process. Moreover gypsy has been shown to have infectious properties, probably due to its envelope gene acquired from a baculovirus. Recently we have also shown that gypsy maternal transmission is reduced in the presence of the endosymbiotic bacterium Wolbachia. These studies demonstrate that gypsy is a unique and powerful model for understanding the endogenization of retroviruses.