Genomic instability of I elements of Drosophila melanogaster in absence of dysgenic crosses

The effects of transpositions of functional I retrotransposons depend on the conditions and dose of parental exposure

PURPOSE

Transposable elements (TEs) cause destabilization of animal genomes. I retrotransposons of Drosophila melanogaster, as well as human LINE1 retrotransposons, are sources of intra- and interindividual diversity and responses to the action of internal and external factors. The aim of this study was to investigate the response to irradiation for the offspring of Drosophila melanogaster with the increased activity of inherited functional I elements.

MATERIALS AND METHODS

The material used was dysgenic Drosophila females with active I retrotransposons obtained as a result of crossing irradiated/non-irradiated parents of a certain genotype. Non-dysgenic females (without functional I elements) were used as controls. The effects of different conditions (irradiation of both parents simultaneously or separately) and doses (1-100 Gy) of parental irradiation have been assessed by analyzing SF-sterility, DNA damage and lifespan. The presence of full-size I retrotransposons was determined by PCR analysis.

RESULTS

The maternal exposure and exposure of both parents are efficient in contrast with paternal exposure. Irradiation of mothers reduces the reproductive potential and viability of their female offspring which undergo high activity of functional I retrotransposons. Though I retrotranspositions negatively affect the female gonads, irradiation of the paternal line can increase the lifespan of SF-sterile females. Radiation stress in the range of 1-100 Gy increases DNA fragmentation in both somatic and germ cells of the ovaries with high I-retrotransposition.

CONCLUSIONS

These results allow for the specificity of the radiation-induced behavior of I retrotransposons and their role in survival under conditions of strong radiation stress.

Rapid Low-Cost Assembly of the Drosophila melanogaster Reference Genome Using Low-Coverage, Long-Read Sequencing

Accurate and comprehensive characterization of genetic variation is essential for deciphering the genetic basis of diseases and other phenotypes. A vast amount of genetic variation stems from large-scale sequence changes arising from the duplication, deletion, inversion, and translocation of sequences. In the past 10 years, high-throughput short reads have greatly expanded our ability to assay sequence variation due to single nucleotide polymorphisms. However, a recent de novo assembly of a second Drosophila melanogaster reference genome has revealed that short read genotyping methods miss hundreds of structural variants, including those affecting phenotypes. While genomes assembled using high-coverage long reads can achieve high levels of contiguity and completeness, concerns about cost, errors, and low yield have limited widespread adoption of such sequencing approaches. Here we resequenced the reference strain of D. melanogaster (ISO1) on a single Oxford Nanopore MinION flow cell run for 24 hr. Using only reads longer than 1 kb or with at least 30x coverage, we assembled a highly contiguous de novo genome. The addition of inexpensive paired reads and subsequent scaffolding using an optical map technology achieved an assembly with completeness and contiguity comparable to the D. melanogaster reference assembly. Comparison of our assembly to the reference assembly of ISO1 uncovered a number of structural variants (SVs), including novel LTR transposable element insertions and duplications affecting genes with developmental, behavioral, and metabolic functions. Collectively, these SVs provide a snapshot of the dynamics of genome evolution. Furthermore, our assembly and comparison to the D. melanogaster reference genome demonstrates that high-quality de novo assembly of reference genomes and comprehensive variant discovery using such assemblies are now possible by a single lab for under $1,000 (USD).

Natural variation of piRNA expression affects immunity to transposable elements

In the Drosophila germline, transposable elements (TEs) are silenced by PIWI-interacting RNA (piRNA) that originate from distinct genomic regions termed piRNA clusters and are processed by PIWI-subfamily Argonaute proteins. Here, we explore the variation in the ability to restrain an alien TE in different Drosophila strains. The I-element is a retrotransposon involved in the phenomenon of I-R hybrid dysgenesis in Drosophila melanogaster. Genomes of R strains do not contain active I-elements, but harbour remnants of ancestral I-related elements. The permissivity to I-element activity of R females, called reactivity, varies considerably in natural R populations, indicating the existence of a strong natural polymorphism in defense systems targeting transposons. To reveal the nature of such polymorphisms, we compared ovarian small RNAs between R strains with low and high reactivity and show that reactivity negatively correlates with the ancestral I-element-specific piRNA content. Analysis of piRNA clusters containing remnants of I-elements shows increased expression of the piRNA precursors and enrichment by the Heterochromatin Protein 1 homolog, Rhino, in weak R strains, which is in accordance with stronger piRNA expression by these regions. To explore the nature of the differences in piRNA production, we focused on two R strains, weak and strong, and showed that the efficiency of maternal inheritance of piRNAs as well as the I-element copy number are very similar in both strains. At the same time, germline and somatic uni-strand piRNA clusters generate more piRNAs in strains with low reactivity, suggesting the relationship between the efficiency of primary piRNA production and variable response to TE invasions. The strength of adaptive genome defense is likely driven by naturally occurring polymorphisms in the rapidly evolving piRNA pathway proteins. We hypothesize that hyper-efficient piRNA production is contributing to elimination of a telomeric retrotransposon HeT-A, which we have observed in one particular transposon-resistant R strain.

Transposon activity in the germline is suppressed by the PIWI-interacting RNA (piRNA) pathway. The resistance of natural Drosophila strains to transposon invasion varies considerably, but the nature of this variability is unknown. We discovered that natural variation in the efficiency of primary piRNA production in the germline causes dramatic differences in the susceptibility to expansion of a newly invaded transposon. A high level of piRNA production in the germline is achieved by increased expression of piRNA precursors. In one of the most transposon-resistant strains, increased content of primary piRNA is observed in both the germline and ovarian somatic cells. We suggest that polymorphisms in piRNA pathway factors are responsible for increased piRNA production. piRNA pathway proteins have been shown to be evolving rapidly under selective pressure. Our data are the first to describe a phenotype that might be caused by this kind of polymorphism. We also demonstrate a likely explanation as to why an overly active piRNA pathway can cause more harm than good in Drosophila: Highly efficient piRNA processing leads to elimination of domesticated telomeric retrotransposons essential for telomere elongation, an effect which has been observed in a natural strain that is extremely resistant to transposon invasion.

Unique transposon landscapes are pervasive across Drosophila melanogaster genomes

To understand how transposon landscapes (TLs) vary across animal genomes, we describe a new method called the Transposon Insertion and Depletion AnaLyzer (TIDAL) and a database of >300 TLs in Drosophila melanogaster (TIDAL-Fly). Our analysis reveals pervasive TL diversity across cell lines and fly strains, even for identically named sub-strains from different laboratories such as the ISO1 strain used for the reference genome sequence. On average, >500 novel insertions exist in every lab strain, inbred strains of the Drosophila Genetic Reference Panel (DGRP), and fly isolates in the Drosophila Genome Nexus (DGN). A minority (<25%) of transposon families comprise the majority (>70%) of TL diversity across fly strains. A sharp contrast between insertion and depletion patterns indicates that many transposons are unique to the ISO1 reference genome sequence. Although TL diversity from fly strains reaches asymptotic limits with increasing sequencing depth, rampant TL diversity causes unsaturated detection of TLs in pools of flies. Finally, we show novel transposon insertions negatively correlate with Piwi-interacting RNA (piRNA) levels for most transposon families, except for the highly-abundant roo retrotransposon. Our study provides a useful resource for Drosophila geneticists to understand how transposons create extensive genomic diversity in fly cell lines and strains.

P elements and the determinants of hybrid dysgenesis have different dynamics of propagation in Drosophila melanogaster populations

Intraspecific hybrid dysgenesis (HD) appears after some strains of D. melanogaster are crossed. The predominant idea is that the movement of transposable P elements causes HD. It is believed that P elements appeared in the D. melanogaster genome in the middle of the last century by horizontal transfer, simultaneously with the appearance of HD determinants. A subsequent simultaneous expansion of HD determinants and P elements occurred. We analyzed the current distribution of HD determinants in natural populations of D. melanogaster and found no evidence of their further spread. However, full-sized P elements were identified in the genomes of all analyzed natural D. melanogaster strains independent of their cytotypes. Thus, the expansion of P elements does not correlate with the expansion of HD determinants. We found that the ovaries of dysgenic females did not contain germ cells despite the equal number of primordial germ cells in early stages in dysgenic and non-dysgenic embryos. We propose that HD does not result from DNA damage caused by P element transposition, but it would be the disruption in the regulation of dysgenic ovarian formation that causes the dysgenic phenotypes.

Comparison of methods for the determination of the transposition rate of mobile elements

There is wide-spread interest in understanding the rate of transposable element movement within populations and between species. A recent study using interprecific crosses between D. buzzatii and D. koepferae indicated that transposition rates in hybrids may be quite high. However, we suggest caution should be taken in this interpretation since AFLP methods to detect transposition events may lead to overestimated rate estimates. Comparative analyses of genome instability received by different methods suggest that transposition rates can be higher in intraspecific crosses compared to interspecific crosses.

Specific activation of an I-like element in Drosophila interspecific hybrids

The non-long terminal repeat (LTR) retrotransposon I, which belongs to the I superfamily of non-LTR retrotransposons, is well known in Drosophila because it transposes at a high frequency in the female germline cells in I-R hybrid dysgenic crosses of Drosophila melanogaster. Here, we report the occurrence and the upregulation of an I-like element in the hybrids of two sister species belonging to the repleta group of the genus Drosophila, D. mojavensis, and D. arizonae. These two species display variable degrees of pre- and postzygotic isolation, depending on the geographic origin of the strains. We took advantage of these features to explore the transposable element (TE) dynamics in interspecific crosses. We fully characterized the copies of this TE family in the D. mojavensis genome and identified at least one complete copy. We showed that this element is transcriptionally active in the ovaries and testes of both species and in their hybrids. Moreover, we showed that this element is upregulated in hybrid males, which could be associated with the male-sterile phenotype.

Euchromatic transposon insertions trigger production of novel Pi- and endo-siRNAs at the target sites in the drosophila germline

The control of transposable element (TE) activity in germ cells provides genome integrity over generations. A distinct small RNA–mediated pathway utilizing Piwi-interacting RNAs (piRNAs) suppresses TE expression in gonads of metazoans. In the fly, primary piRNAs derive from so-called piRNA clusters, which are enriched in damaged repeated sequences. These piRNAs launch a cycle of TE and piRNA cluster transcript cleavages resulting in the amplification of piRNA and TE silencing. Using genome-wide comparison of TE insertions and ovarian small RNA libraries from two Drosophila strains, we found that individual TEs inserted into euchromatic loci form novel dual-stranded piRNA clusters. Formation of the piRNA-generating loci by active individual TEs provides a more potent silencing response to the TE expansion. Like all piRNA clusters, individual TEs are also capable of triggering the production of endogenous small interfering (endo-si) RNAs. Small RNA production by individual TEs spreads into the flanking genomic regions including coding cellular genes. We show that formation of TE-associated small RNA clusters can down-regulate expression of nearby genes in ovaries. Integration of TEs into the 3′ untranslated region of actively transcribed genes induces piRNA production towards the 3′-end of transcripts, causing the appearance of genic piRNA clusters, a phenomenon that has been reported in different organisms. These data suggest a significant role of TE-associated small RNAs in the evolution of regulatory networks in the germline.

Silencing of transposable elements (TEs) in germ cells depends on a distinct class of small RNAs, Piwi-interacting RNAs (piRNAs). TE repression is provided by piRNAs derived from large heterochromatic loci enriched in fragmented TE copies, so-called piRNA clusters. According to the current model, individual TEs and their transcripts are considered merely as targets of cluster-derived primary piRNAs, which exert post-transcriptional and transcriptional silencing in Drosophila. In our work, we show that natural individual transposons become piRNA-generating loci themselves. We came to this conclusion by comparing the ovarian small RNAs and TE insertion sites of two Drosophila strains, which showed that euchromatic target sites of strain-specific TEs generate a number of novel strain-specific piRNAs. This mechanism allows production of additional small RNAs that target active TEs and provide more potent transposon suppression in the germline. Moreover, small RNA production by individual TEs spreads into the flanking genomic regions, which affects the expression of adjacent coding genes and microRNA genes. These data underline the role of individual TEs in a silencing response and explore a new level of TE impact on the gene regulatory networks in the germline.

Quantitatively increased somatic transposition of transposable elements in Drosophila strains compromised for RNAi

In Drosophila melanogaster, small RNAs homologous to transposable elements (TEs) are of two types: piRNA (piwi-interacting RNA) with size 23-29nt and siRNA (small interfering RNA) with size 19-22nt. The siRNA pathway is suggested to silence TE activities in somatic tissues based on TE expression profiles, but direct evidence of transposition is lacking. Here we developed an efficient FISH (fluorescence in Situ hybridization) based method for polytene chromosomes from larval salivary glands to reveal new TE insertions. Analysis of the LTR-retrotransposon 297 and the non-LTR retroposon DOC shows that in the argonaut 2 (Ago2) and Dicer 2 (Dcr2) mutant strains, new transposition events are much more frequent than in heterozygous strains or wild type strains. The data demonstrate that the siRNA pathway represses TE transposition in somatic cells. Nevertheless, we found that loss of one functional copy of Ago2 or Dcr2 increases somatic transpositions of the elements at a lower level depending on the genetic background, suggesting a quantitative role for RNAi core components on mutation frequency.

Drosophila interspecific hybrids phenocopy piRNA-pathway mutants

Hybrids of two Drosophila species show transposable element derepression and piRNA pathway malfunction, revealing adaptive evolution of piRNA pathway components.

The Piwi-interacting RNA (piRNA) pathway defends the germline of animals from the deleterious activity of selfish transposable elements (TEs) through small-RNA mediated silencing. Adaptation to novel invasive TEs is proposed to occur by incorporating their sequences into the piRNA pool that females produce and deposit into their eggs, which then propagates immunity against specific TEs to future generations. In support of this model, the F1 offspring of crosses between strains of the same Drosophila species sometimes suffer from germline derepression of paternally inherited TE families, caused by a failure of the maternal strain to produce the piRNAs necessary for their regulation. However, many protein components of the Drosophila piRNA pathway exhibit signatures of positive selection, suggesting that they also contribute to the evolution of host genome defense. Here we investigate piRNA pathway function and TE regulation in the F1 hybrids of interspecific crosses between D. melanogaster and D. simulans and compare them with intraspecific control crosses of D. melanogaster. We confirm previous reports showing that intraspecific crosses are characterized by derepression of paternally inherited TE families that are rare or absent from the maternal genome and piRNA pool, consistent with the role of maternally deposited piRNAs in shaping TE silencing. In contrast to the intraspecific cross, we discover that interspecific hybrids are characterized by widespread derepression of both maternally and paternally inherited TE families. Furthermore, the pattern of derepression of TE families in interspecific hybrids cannot be attributed to their paucity or absence from the piRNA pool of the maternal species. Rather, we demonstrate that interspecific hybrids closely resemble piRNA effector-protein mutants in both TE misregulation and aberrant piRNA production. We suggest that TE derepression in interspecific hybrids largely reflects adaptive divergence of piRNA pathway genes rather than species-specific differences in TE-derived piRNAs.

Eukaryotic genomes contain large quantities of transposable elements (TEs), short self-replicating DNA sequences that can move within the genome. The selfish replication of TEs has potentially drastic consequences for the host, such as disruption of gene function, induction of sterility, and initiation or exacerbation of some cancers. Like the adaptive immune system that defends our bodies against pathogens, the Piwi-interacting RNA (piRNA) pathway defends animal genomes against the harmful effects of TEs. Fundamental to piRNA-mediated defense is the production of small noncoding RNAs that act like antibodies to target replicating TEs for destruction by piRNA-effector proteins. piRNAs are expected to diverge rapidly between species in response to genome infection by increasingly disparate TEs. Here, we tested this hypothesis by examining how differences in piRNAs between two species of fruit fly relate to TE “immunity” in their hybrid offspring. Because piRNAs are maternally deposited, we expected excessive replication of paternal TEs in hybrids. Surprisingly, we observe increased activity of both maternal and paternal TEs, together with defects in piRNA production that are reminiscent of piRNA effector-protein mutants. Our observations reveal that piRNA effector-proteins do not function properly in hybrids, and we propose that adaptive evolution among piRNA effector-proteins contributes to host genome defense and leads to the functional incompatibilities that we observe in hybrids.

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.