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

SQuIRE reveals locus-specific regulation of interspersed repeat expression

Transposable elements (TEs) are interspersed repeat sequences that make up much of the human genome. Their expression has been implicated in development and disease. However, TE-derived RNA-seq reads are difficult to quantify. Past approaches have excluded these reads or aggregated RNA expression to subfamilies shared by similar TE copies, sacrificing quantitative accuracy or the genomic context necessary to understand the basis of TE transcription. As a result, the effects of TEs on gene expression and associated phenotypes are not well understood. Here, we present Software for Quantifying Interspersed Repeat Expression (SQuIRE), the first RNA-seq analysis pipeline that provides a quantitative and locus-specific picture of TE expression (https://github.com/wyang17/SQuIRE). SQuIRE is an accurate and user-friendly tool that can be used for a variety of species. We applied SQuIRE to RNA-seq from normal mouse tissues and a Drosophila model of amyotrophic lateral sclerosis. In both model organisms, we recapitulated previously reported TE subfamily expression levels and revealed locus-specific TE expression. We also identified differences in TE transcription patterns relating to transcript type, gene expression and RNA splicing that would be lost with other approaches using subfamily-level analyses. Altogether, our findings illustrate the importance of studying TE transcription with locus-level resolution.

Protecting and Diversifying the Germline

Gametogenesis represents the most dramatic cellular differentiation pathways in both female and male flies. At the genome level, meiosis ensures that diploid germ cells become haploid gametes. At the epigenome level, extensive changes are required to turn on and shut off gene expression in a precise spatiotemporally controlled manner. Research applying conventional molecular genetics and cell biology, in combination with rapidly advancing genomic tools have helped us to investigate (1) how germ cells maintain lineage specificity throughout their adult reproductive lifetime; (2) what molecular mechanisms ensure proper oogenesis and spermatogenesis, as well as protect genome integrity of the germline; (3) how signaling pathways contribute to germline-soma communication; and (4) if such communication is important. In this chapter, we highlight recent discoveries that have improved our understanding of these questions. On the other hand, restarting a new life cycle upon fertilization is a unique challenge faced by gametes, raising questions that involve intergenerational and transgenerational epigenetic inheritance. Therefore, we also discuss new developments that link changes during gametogenesis to early embryonic development-a rapidly growing field that promises to bring more understanding to some fundamental questions regarding metazoan development.

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.

Evolution and dynamics of small RNA response to a retroelement invasion in Drosophila

Although small RNAs efficiently control transposition activity of most transposons in the host genome, such an immune system is not always applicable against a new transposon's invasions. Here, we explored a possibility to introduce potentially mobile copy of the Penelope retroelement previously implicated in hybrid dysgenesis syndrome in Drosophila virilis into the genomes of two distant Drosophila species. The consequences of such introduction were monitored at different phases after experimental colonization as well as in D. virilis species, which is apparently in the process of ongoing Penelope invasion. We investigated the expression of Penelope and biogenesis of Penelope-derived small RNAs in D. virilis and D. melanogaster strains originally lacking active copies of this element after experimental Penelope invasion. These strains were transformed by constructs containing intact Penelope copies. We show that immediately after transformation, which imitates the first stage of retroelement invasion, Penelope undergoes transposition predominantly in somatic tissues, and may produce siRNAs that are apparently unable to completely silence its activity. However, at the later stages of colonization Penelope copies may jump into one of the piRNA-clusters, which results in production of homologous piRNAs that are maternally deposited and can silence euchromatic transcriptionally active copies of Penelope in trans and, hence, prevent further amplification of the invader in the host genome. Intact Penelope copies and different classes of Penelope-derived small RNAs were found in most geographical strains of D. virilis collected throughout the world. Importantly, all strains of this species containing full-length Penelope tested do not produce gonadal sterility in dysgenic crosses and, hence, exhibit neutral cytotype. To understand whether RNA interference mechanism able to target Penelope operates in related species of the virilis group, we correlated the presence of full-length and potentially active Penelope with the occurrence of piRNAs homologous to this transposable element in the ovaries of species comprising the group. It was demonstrated that Penelope-derived piRNAs are present in all virilis group species containing full-length but transcriptionally silent copies of this element that probably represent the remnants of its previous invasions taking place in the course of the virilis species divergent evolution.

Expression of Drosophila virilis retroelements and role of small RNAs in their intrastrain transposition

Transposition of two retroelements (Ulysses and Penelope) mobilized in the course of hybrid dysgenesis in Drosophila virilis has been investigated by in situ hybridization on polytene chromosomes in two D. virilis strains of different cytotypes routinely used to get dysgenic progeny. The analysis has been repeatedly performed over the last two decades, and has revealed transpositions of Penelope in one of the strains, while, in the other strain, the LTR-containing element Ulysses was found to be transpositionally active. The gypsy retroelement, which has been previously shown to be transpositionally inactive in D. virilis strains, was also included in the analysis. Whole mount is situ hybridization with the ovaries revealed different subcellular distribution of the transposable elements transcripts in the strains studied. Ulysses transpositions occur only in the strain where antisense piRNAs homologous to this TE are virtually absent and the ping-pong amplification loop apparently does not take place. On the other hand small RNAs homologous to Penelope found in the other strain, belong predominantly to the siRNA category (21nt), and consist of sense and antisense species observed in approximately equal proportion. The number of Penelope copies in the latter strain has significantly increased during the last decades, probably because Penelope-derived siRNAs are not maternally inherited, while the low level of Penelope-piRNAs, which are faithfully transmitted from mother to the embryo, is not sufficient to silence this element completely. Therefore, we speculate that intrastrain transposition of the three retroelements studied is controlled predominantly at the post-transcriptional level.

Evolution of a distinct genomic domain in Drosophila: comparative analysis of the dot chromosome in Drosophila melanogaster and Drosophila virilis

The distal arm of the fourth ("dot") chromosome of Drosophila melanogaster is unusual in that it exhibits an amalgamation of heterochromatic properties (e.g., dense packaging, late replication) and euchromatic properties (e.g., gene density similar to euchromatic domains, replication during polytenization). To examine the evolution of this unusual domain, we undertook a comparative study by generating high-quality sequence data and manually curating gene models for the dot chromosome of D. virilis (Tucson strain 15010-1051.88). Our analysis shows that the dot chromosomes of D. melanogaster and D. virilis have higher repeat density, larger gene size, lower codon bias, and a higher rate of gene rearrangement compared to a reference euchromatic domain. Analysis of eight "wanderer" genes (present in a euchromatic chromosome arm in one species and on the dot chromosome in the other) shows that their characteristics are similar to other genes in the same domain, which suggests that these characteristics are features of the domain and are not required for these genes to function. Comparison of this strain of D. virilis with the strain sequenced by the Drosophila 12 Genomes Consortium (Tucson strain 15010-1051.87) indicates that most genes on the dot are under weak purifying selection. Collectively, despite the heterochromatin-like properties of this domain, genes on the dot evolve to maintain function while being responsive to changes in their local environment.

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.

Drosophila viruses and the study of antiviral host-defense

The fruit fly Drosophila melanogaster is a powerful model to study host-pathogen interactions. Most studies so far have focused on extracellular pathogens such as bacteria and fungi. More recently, viruses have come to the front, and RNA interference was shown to play a critical role in the control of viral infections in drosophila. We review here our current knowledge on drosophila viruses. A diverse set of RNA viruses belonging to several families (Rhabdoviridae, Dicistroviridae, Birnaviridae, Reoviridae, Errantiviridae) has been reported in D. melanogaster. By contrast, no DNA virus has been recovered up to now. The drosophila viruses represent powerful tools to study virus-cell interactions in vivo. Analysis of the literature however reveals that for many of them, important gaps exist in our understanding of their replication cycle, genome organization, morphology or pathogenesis. The data obtained in the past few years on antiviral defense mechanisms in drosophila, which point to evolutionary conserved pathways, highlight the potential of the D. melanogaster model to study antiviral innate immunity and to better understand the complex interaction between arthropod-borne viruses and their insect vectors.

Gypsy endogenous retrovirus maintains potential infectivity in several species of Drosophilids

Background

Sequences homologous to the gypsy retroelement from Drosophila melanogaster are widely distributed among drosophilids. The structure of gypsy includes an open reading frame resembling the retroviral gene env, which is responsible for the infectious properties of retroviruses.

Results

In this study we report molecular and phylogeny analysis of the complete env gene from ten species of the obscura group of the genus Drosophila and one species from the genus Scaptomyza.

Conclusion

The results indicate that in most cases env sequences could produce a functional Env protein and therefore maintain the infectious capability of gypsy in these species.

Multiple invasions of Errantivirus in the genus Drosophila

Aiming to contribute to the knowledge of the evolutionary history of Errantivirus, a phylogenetic analysis of the env gene sequences of Errantivirus gypsy, gtwin, gypsy2, gypsy3, gypsy4 and gypsy6 was carried out in 33 Drosophilidae species. Most sequences were obtained from in silico searches in the Drosophila genomes. The complex evolutionary pattern reported by other authors for the gypsy retroelement was also observed in the present study, including vertical transmission, ancestral polymorphism, stochastic loss and horizontal transfer. Moreover, the elements gypsy2, gypsy3, gypsy4 and gypsy6 were shown to have followed an evolutionary model that is similar to gypsy. Fifteen new possible cases of horizontal transfer were suggested. The infectious potential of these elements may help elucidate the evolutionary scenario described in the present study.

The LTR retrotransposon micropia in the cardini group of Drosophila (Diptera: Drosophilidae): a possible case of horizontal transfer

The presence of the micropia retroelement from the Ty1-copia family of LTR retroelements was investigated in three species of the Drosophila cardini group. Southern blot analysis suggested the existence of at least four micropia copies in the genomes of D. cardinoides, D. neocardini and D. polymorpha populations. The high sequence similarity between dhMiF2 and Dm11 clones (micropia retroelements isolated from D. hydei and D. melanogaster, respectively) with micropia sequences amplified from D. cardini group genome supports the hypothesis that this retroelement plays an active role in horizontal transfer events between D. hydei and the D. cardini group.

Polymorphism of canonical and noncanonical gypsy sequences in different species of Drosophila melanogaster subgroup: possible evolutionary relations

Mobile genetic elements constitute a substantial part of eukaryotic genome and play an important role in its organization and functioning. Co-evolution of retrotransposons and their hosts resulted in the establishment of control systems employing mechanisms of RNA interference that seem to be impossible to evade. However, "active" copies of endogenous retrovirus gypsy escape cellular control in some cases, while its evolutionary elder "inactive" variants do not. To clarify the evolutionary relationship between "active" and "inactive" gypsy we combined two approaches: the analysis of gypsy sequences, isolated from G32 Drosophila melanogaster strain and from different Drosophila species of the melanogaster subgroup, as well as the study of databases, available on the Internet. No signs of "intermediate" (between "active" and "inactive") gypsy form were found in GenBank, and four full-size G32 gypsy copies demonstrated a convergence that presumably involves gene conversion. No "active" gypsy were revealed among PCR generated gypsy ORF3 sequences from the various Drosophila species indicating that "active" gypsy appeared in some population of D. melanogaster and then started to spread out. Analysis of sequences flanking gypsy variants in G32 revealed their predominantly heterochromatic location. Discrepancy between the structure of actual gypsy sites in G32 and corresponding sequences in database might indicate significant inter-strain heterochromatin diversity.

A novel repeat-associated small interfering RNA-mediated silencing pathway downregulates complementary sense gypsy transcripts in somatic cells of the Drosophila ovary

Replication of the gypsy endogenous retrovirus involves contamination of the female germ line by adjacent somatic tissues. This is prevented by flam, an as-yet-uncloned heterochromatic pericentromeric locus, at the level of transcript accumulation in these somatic ovarian tissues. We tested the effect of a presumptive RNA silencing mechanism on the accumulation of RNAs produced by constructs containing various gypsy sequences and report that the efficiency of silencing is indeed correlated with the amount of complementary RNAs, 25 to 30 nucleotides in length, in the ovary. For instance, while these RNAs were found to display a three- to fivefold excess of the antisense strands, only the transcripts that contain the complementary sense gypsy sequences could be repressed, indicating that they are targeted at the RNA, not DNA, level. Their size and asymmetry in strand polarity are typical of the novel repeat-associated small interfering RNA (rasiRNA)-mediated pathway, recently suspected to prevent the deleterious expression of selfish DNA specifically in the germ line. Unlike microRNAs (but like rasiRNAs and, surprisingly, siRNAs as well), gypsy rasiRNAs are modified at the 3' end. The rasiRNA-associated protein Piwi (but not Aub) is required for gypsy silencing, whereas Dicer-2 (which makes siRNAs) is not. In contrast, piwi, aub, and flam do not appear to affect somatic siRNA-mediated silencing. The amount of gypsy rasiRNAs is genetically determined by the flam locus in a provirus copy number-independent manner and is triggered in the somatic tissues by some pericentromeric provirus(es), which are thereby able to protect the germ line from retroviral invasion.

Evolutionary pattern of the gtwin retrotransposon in the Drosophila melanogaster subgroup

The gtwin retrotransposon was recently discovered in the Drosophila melanogaster genome and it is evolutionarily closer to gypsy endogenous retrovirus. This study has identified gtwin homologous sequences in the genome of D. simulans, D. sechellia, D. erecta and D. yakuba by performing homology searches against the public genome database of Drosophila species. The phylogenetic analyses of the gtwin env gene sequences of these species have shown some incongruities with the host species phylogeny, suggesting some horizontal transfer events for this retroelement. Moreover, we reported the existence of DNA sequences putatively encoding full-length Env proteins in the genomes of Drosophila species other than D. melanogaster. The results suggest that the gtwin element may be an infectious retrovirus able to invade the genome of new species, supporting the gtwin evolutionary picture shown in this work.

A genome-wide screening of BEL-Pao like retrotransposons in Anopheles gambiae by the LTR_STRUC program

The advanced status of assembly of the nematoceran Anopheles gambiae genomic sequence allowed us to perform a wide genome analysis to looking at the presence of Long Terminal Repeats (LTRs) in the range of 10 kb by means of the LTR_STRUC tool. More than three hundred sequences were retrieved and 210 were treated as putative complete retrotransposons that were individually analysed with respect to known retrotransposons of A. gambiae and D. melanogaster. The results show that the vast majority of the retrotransposons analysed belong to the Ty3/gypsy class and only 8% to the Ty1/copia class. In addition, phylogenetic analysis allowed us to characterize in more detail the relationship of a large BEL-Pao lineage in which a single family was shown to harbour an additional env gene.

Endogenous retroviruses and animal reproduction

Endogenous retroviruses (ERV), as part of the host genetic heritage, are transmissible to the next generation in a Mendelian way. Their abundance in animal genomes and their expression primarily detected in germ cells, embryonic tissues and cancer cell lines, raised the question of their biological significance. This article reviews the possible role of ERVs in the physiology and diseases of animal reproduction, from Drosophila to human. In males, there is no trivial involvement of ERVs in a physiological process. Conversely, a spermatogenesis defect was associated in the human male with HERV-K expression and HERV15-induced chromosomal alteration, leading to cancer and infertility, respectively. In females, the study of insect ERVs (IERV) pointed out the overlap between genetics and virology with the genetic-dependent regulation of ZAM and the non-infectious and infectious life cycles of gypsy. The pattern of ERVs expression in rodent, ovine and human females suggest a hormone-dependent mechanism consistent with the mammalian oestrus cycle regulation. The differentiation of the mammary epithelium and breast tumorigenesis involving the mouse mammary tumour viruses (MMTV) illustrate the intimate connection between endogenous and exogenous retroviruses. Last, as a major site of ERVs transcription, placenta contributed to our understanding of ERVs modulation of neighbouring gene expression. As an interface, i.e. a site of conflicts and exchanges, placenta should resist infection and protect the foetus against the maternal immune system. Retroviral envelopes could theoretically provide such features due to receptor interference, immunosuppression and fusion properties, as shown by the HERV-W envelope involved in the syncytiotrophoblast formation. We conclude with an insight on the evolutionary and epigenetic consequences of the relationships of ERV guests with their animal hosts.

Evidence for a piwi-dependent RNA silencing of the gypsy endogenous retrovirus by the Drosophila melanogaster flamenco gene

In Drosophila melanogaster, the endogenous retrovirus gypsy is repressed by the functional alleles (restrictive) of an as-yet-uncloned heterochromatic gene called flamenco. Using gypsy-lacZ transcriptional fusions, we show here that this repression takes place not only in the follicle cells of restrictive ovaries, as was previously observed, but also in restrictive larval female gonads. Analyses of the role of gypsy cis-regulatory sequences in the control of gypsy expression are also presented. They rule out the hypothesis that gypsy would contain a single binding region for a putative Flamenco repressor. Indeed, the ovarian expression of a chimeric yp3-lacZ construct was shown to become sensitive to the Flamenco regulation when any of three different 5'-UTR gypsy sequences (ranging from 59 to 647 nucleotides) was incorporated into the heterologous yp3-lacZ transcript. The piwi mutation, which is known to affect RNA-mediated homology-dependent transgene silencing, was also shown to impede the repression of gypsy in restrictive female gonads. Finally, a RNA-silencing model is also supported by the finding in ovaries of short RNAs (25-27 nucleotides long) homologous to sequences from within the gypsy 5'-UTR.

Complex evolution of gypsy in Drosophilid species

In an endeavor to contribute to the comprehension of the evolution of transposable elements (TEs) in the genome of host species, we investigated the phylogenetic relationships of sequences homologous to the retrotransposon gypsy of Drosophila melanogaster in 19 species of Drosophila, in Scaptodrosophila latifasciaeformis, and in Zaprionus indianus. This phylogenetic study was based on approximately 500 base pairs of the env gene. Our analyses showed considerable discrepancy between the phylogeny of gypsy elements and the relationship of their host species, and they allow us to infer a complex evolutionary pattern that could include ancestral polymorphism, vertical transmission, and several cases of horizontal transmission.

Comparative sequence analysis and predictions for the envelope glycoproteins of insect endogenous retroviruses

Insect endogenous retroviruses (IERVs) are present in the genome of several species. Previous studies have shown a relationship between the envelope glycoproteins (Envs) and fusion proteins (FPs) of several baculoviruses. We used this sequence similarity to predict fusion domains in the Envs of IERVs. We suggest that FPs and Envs share several specific sequence and structural motifs with other RNA viruses in the viral transmembrane protein superfamily.

Drosophila germline invasion by the endogenous retrovirus gypsy: involvement of the viral env gene

The endogenous retrovirus gypsy is expressed at high levels in mutant flamenco female flies. Gypsy viral particles extracted from such flies can infect naive flamenco individuals raised in the presence of these extracts mixed into their food. This results in the integration of new proviruses into the germline genome. These proviruses can then increase their copy number by (1) expression in the flamenco female somatic cells, (2) transfer into the oocyte and (3) integration into the genome of the progeny. Surprisingly, unlike the infection observed in the feeding experiments, this strategy of endogenous proviral multiplication does not seem to involve the expression of the viral env gene.