Whole genome sequencing in Drosophila virilis identifies Polyphemus, a recently activated Tc1-like transposon with a possible role in hybrid dysgenesis
© Blumenstiel; licensee BioMed Central Ltd. 2014
Received: 6 October 2013
Accepted: 28 January 2014
Published: 20 February 2014
Hybrid dysgenic syndromes in Drosophila have been critical for characterizing host mechanisms of transposable element (TE) regulation. This is because a common feature of hybrid dysgenesis is germline TE mobilization that occurs when paternally inherited TEs are not matched with a maternal pool of silencing RNAs that maintain transgenerational TE control. In the face of this imbalance TEs become activated in the germline and can cause F1 sterility. The syndrome of hybrid dysgenesis in Drosophila virilis was the first to show that the mobilization of one dominant TE, the Penelope retrotransposon, may lead to the mobilization of other unrelated elements. However, it is not known how many different elements contribute and no exhaustive search has been performed to identify additional ones. To identify additional TEs that may contribute to hybrid dysgenesis in Drosophila virilis, I analyzed repeat content in genome sequences of inducer and non-inducer lines.
Here I describe Polyphemus, a novel Tc1-like DNA transposon, which is abundant in the inducer strain of D. virilis but highly degraded in the non-inducer strain. Polyphemus expression is also increased in the germline of progeny of the dysgenic cross relative to reciprocal progeny. Interestingly, like the Penelope element, it has experienced recent re-activation within the D. virilis lineage.
Here I present the results of a comprehensive search to identify additional factors that may cause hybrid dysgenesis in D. virilis. Polyphemus, a novel Tc1-like DNA transposon, has recently become re-activated in Drosophila virilis and likely contributes to the hybrid dysgenesis syndrome. It has been previously shown that the Penelope element has also been re-activated in the inducer strain. This suggests that TE co-reactivation within species may synergistically contribute to syndromes of hybrid dysgenesis.
KeywordsHybrid dysgenesis Drosophila virilis Transposable element Penelope piRNA Genome instability Epigenetics
Hybrid dysgenesis, a syndrome of sterility and increased mutation in crosses between different strains of the same species, was first shown in Drosophila melanogaster to be driven by P elements that are inherited paternally, but not maternally [1–4]. Activation of this DNA transposon subsequently leads to germline DNA damage and sterility [5–8]. In D. melanogaster, transposable element (TE) mediated syndromes of hybrid dysgenesis are also driven by the I element retrotransposon  and the hobo DNA transposon . An important syndrome of hybrid dysgenesis has also been characterized in D. virilis[11–14]. This syndrome is significant as it is accompanied by mobilization of different, unrelated transposable element families [13, 15]. Critically, even elements such as the Ulysses retrotransposon that are evenly distributed between strains become mobilized in this cross. Recent studies using genome sequencing approaches indicate that P elements also induce the mobilization of other elements in D. melanogaster. Thus, co-mobilization of TEs may be a common feature of dysgenic syndromes but the mechanism by which TE co-mobilization occurs is poorly understood. One mechanism that has been proposed is that DNA damage caused by the activation of one TE family disrupts piRNA silencing mechanisms in the germline via the DNA damage response . In this model, disrupted piRNA silencing in turn leads to activation of normally repressed and unrelated TEs. Alternatively, DNA damage arising from transposition may drive activation of TEs through other mechanisms. Additionally, it has been proposed that, like viruses, TEs may encode suppressors of RNA silencing . In this case, the expression of a suppressor of RNA silencing encoded by a single activated TE could lead to global TE de-repression. Finally, it is important to consider the possibility that multiple TE families may be more abundant within inducer strains . Thus, the activation of multiple TE families in a dysgenic syndrome may also be explained by independent mechanisms acting across each family.
To distinguish among these hypotheses, it is critical to define the landscape of TE copy number imbalance between inducer and non-inducer strains in hybrid dysgenic syndromes. Many previous studies indicate that the Penelope element is likely to be the main driver of hybrid dysgenesis in Drosophila virilis. It is the only known element with multiple, active copies in the inducer strain and for which active copies are entirely absent from the non-inducer strain [11, 17, 18]. In addition, its expression is greatly increased in the gonads of dysgenic progeny and injection of embryos with Penelope constructs can lead to increased incidence of TE mediated mutation [11, 14, 19]. However, additional studies indicate that while Penelope may be the dominant cause of sterility, other factors may also contribute. For one, some strains of D. virilis that behave as neutral strains - maternally protecting against dysgenesis but not inducing it paternally - lack piRNAs from the Penelope element in their ovaries . If Penelope is the sole cause of dysgenesis, it is difficult to explain how these strains protect against the induction of dysgenesis since mothers would be unable to provide Penelope piRNA to the next generation. Second, two additional TEs - Helena and Paris - also show high abundance of active, euchromatic copies in the inducer strain and lower abundance in the non-inducer strain. Evidence suggests these two elements also contribute to the sterility phenotype . Whether these three elements act synergistically to cause sterility is not known. It is also not known whether they jointly contribute to the mobilization of other elements such as Ulysses.
The discovery of these candidate inducer elements - Penelope, Helena, and Paris - was facilitated by the recovery of TE insertions that gave rise to visible mutations during dysgenic co-mobilization. Thus, it has not been clear whether additional elements may contribute to the dysgenic syndrome in D. virilis. Here, I present the first systematic effort to identify additional TEs that may cause hybrid dysgenesis in D. virilis.
Results and discussion
Penelope RNA-seq analysis: ANOVA results
Treat X age
Polyphemus RNA-seq analysis: ANOVA results
Treat X age
Here I describe Polyphemus, a new Tc1-like transposable element in D. virilis that may contribute to the hybrid dysgenic syndrome. Whereas highly similar copies are abundant in the inducer strain, only degraded copies are found in the non-inducer strain. Nonetheless, the lack of active copies in the non-inducer strain does not suggest an entirely new invasion of Polyphemus into D. virilis. Instead, phylogenetic analysis indicates that different lineages of Polyphemus have persisted in D. virilis for many years. Against this history, it appears that a Polyphemus variant has now become re-activated. This re-activation may have occurred via a horizontal transfer event from another member of the D. virilis group that has maintained an active Polyphemus lineage since divergence from D. virilis. Alternatively, an active lineage of Polyphemus may have continuously persisted in D. virilis. In this case, individuals within the species would have been segregating with respect to rare active copies of Polyphemus. In this scenario, strains or populations that maintain rare active copies may have functioned as reservoirs for later re-activation of Polyphemus. Interestingly, the Penelope element shows a similar pattern. Active copies of Penelope are highly abundant in the inducer strain, but this is not because Penelope is entirely new to the D. virilis lineage. The non-inducer strain 9 also possesses old, inactive copies and different variant lineages of Penelope are seen in different species of the D. virilis group. Thus, it appears that both of these elements have become re-activated from within the D. virilis group and now may jointly cause hybrid dysgenesis.
In D. melanogaster, three different element families are known to cause hybrid dysgenesis. In the most well understood P-M system, the P element invaded via horizontal from the distant D. willistoni species . However, for the I-R system and the hobo system, both I elements  and hobo elements  have remnant copies residing in the D. melanogaster genome. This suggests that dysgenic syndromes may frequently result from re-activation of TE lineages that, in contrast to the P elements, have been long-term genomic residents. Interestingly, the presence of multiple different elements, such as Paris and Helena, that are also in excess in the inducer strain of D. virilis seems to indicate that TE control has been diminished in strain 160. Perhaps this has occurred by the same mechanism that leads to TE co-mobilization. If so, Penelope might be the cause of this TE excess in strain 160, but might not be the sole proximate cause of sterility in hybrid dysgenesis. For this reason, the mechanisms that are responsible for TE co-mobilization in dysgenic syndromes may also be relevant to understanding global TE dynamics within species after one or more TEs becomes re-activated.
Expression of Polyphemus is higher in the germline of dysgenic progeny, though not to the same magnitude as Penelope. In light of this, it is important to note that the sterility syndrome is evident early in development , not at the time that gonadal expression is typically measured. For this reason, TE expression in adult females that have escaped sterility - a necessary condition for measuring germline gene expression - may not be a perfect proxy for understanding TE expression early in development. The change in Penelope expression during the aging process indicates that TE activity is likely to be dynamic in the life of a dysgenic F1. For this reason, it will be critical to determine the patterns of germline activity for all four TEs early in development. Combined with genetic approaches, this may elucidate the causal factors of sterility and TE co-mobilization in the hybrid dysgenic syndrome of Drosophila virilis.
Genome sequencing and analysis of Polyphemus
DNA was collected from wandering third instar larvae from strain 160 (the inducer strain) and strain 9 (the non-inducer strain). DNA was then sonicated and fragments between 400 to 500 bp were selected for Illumina library preparation. Each library was 100 bp, paired-end sequenced on an individual lane of a GAII using, yielding 43.7 million (strain 160) and 37.6 million (strain 9) read pairs. Since pairs are not independent samples, only single ends of each pair were selected for this analysis. Single reads were quality trimmed using the Sickle application with default settings. Subsequent to quality trimming, reads were mapped using BWA-MEM to a pre-computed PILER library of repeat sequences from D. virilis (ftp://ftp.flybase.net/genomes/aaa/transposable_elements/PILER-DF). Total read counts for each centroid were normalized to the total number of reads mapping back to the reference and the ratio of 160:9 normalized reads was determined. From this, the 1720 bp PILER centroid.25.39 was identified as highly abundant in strain 160 but not strain 9. This centroid was used in a blastn search of the D. virilis reference genome and many nearly exact copies were identified. Three of these elements were extracted with 1000 bp of flanking sequence from each side. Reciprocal pairwise blast between these three fragments identified a core sequence 1,708 bp long of near identify among these fragments. Further annotation of this sequence was performed using Geneious. Phylogenetic analysis of Tc1/mariner members was performed using MrBayes on a MUSCLE amino acid alignment until the average standard deviation of split frequencies was less than 0.01. Analysis of heterogeneity among fragments was performed using piledriver to examine sequence heterogeneity of reads mapping to the respective TE. From the piledriver output, the frequency of the most common base at each position was calculated. Phylogenetic analysis was performed by collecting all blastn fragments from the D. virilis reference genome (excluding inverted repeats) with E values better than E-5. Fragments smaller than 100 bp were removed from the blast output and the blast output anchored to the active sequence was used as an alignment. Tree searching was performed using GARLI on the CIPRES server (http://www.phylo.org/) with a GTR model, no rate heterogeneity and empirical base frequencies. Terminal branch length distributions were extracted from the resulting tree file.
RNA was collected from 0 to 2-hour-old embryos that were laid by F1 females of the dysgenic (strain 9 mothers and strain 160 fathers) and non-dysgenic (strain 160 mothers and strain 9 fathers) crosses. Embryos of F1 females were chosen to avoid the confounding effects that are presented by ovaries that may differ in somatic vs. germline tissue representation in dysgenic crosses. Dysgenic and non-dysgenic crosses were set up en masse and hundreds of F1 females were collected soon after eclosion. Hundreds of reciprocal F1 males from the same crosses were also collected, combined in equal proportions, then reallocated equally to the collected dysgenic and non-dysgenic F1 females in mating cages. This was done to ensure a sufficient egg lay from the F1 females escapers of the dysgenic cross. Such females, if only provided their dysgenic brothers, lay few eggs. Providing equivalent but mixed populations of reciprocal males to the female pools also ensures proper genetic control over paternal effects en masse since reciprocal males of the dysgenic and non-dysgenic crosses are genetically different. From these cages, eggs were collected over 0 to 2-hour egg lay durations and each collection was immediately flash frozen in liquid nitrogen. Collections were pooled into two different age classes based on the age of the mother: 12 to 16 days old (about 1 week after sexual maturity) and 19 to 21 days old.
Pooled RNA was collected and Illumina libraries were generated for single-end, 50 bp RNA-seq. To control for index effects each RNA sample was used to generate two index libraries for a total of eight libraries, each of which were run in replicate on two different lanes. Trimming and filtering was performed using the Galaxy server (https://usegalaxy.org/). Up to 16 bp were quality trimmed from the 3' end from each read and remaining reads with more than 2 bp with quality less than 20 were removed. Trimmed and filtered reads were mapped to TE sequences using CLC with mismatch, insertion and deletion scores equal to 2.3 and 3, respectively. Expression levels were measured by the number of reads that mapped each TE normalized by the total number of reads mapping to the reference genome. ANOVA was performed in R using the aov command and a model that included the effects of treatment (dysgenic vs. non-dysgenic), age, index, and an interaction between age and treatment.
Generalised time reversible
I thank Chris Harrison for generating the RNA-seq libraries. DNA library prep and sequencing was performed by the David H. Murdock Research Institute. Funding was provided by NSF MCB-1022165, an internal grant through NIH-P20GM103638, and the University of Kansas.
- Engels WR: The P-Family of transposable elemens in Drosophila . Annu Rev Genet 1983, 17: 313-344.View ArticleGoogle Scholar
- Bingham PM, Kidwell MG, Rubin GM: The molecular basis of P-M dysgenesis - the role of the P-element, a P-strain-specific transposon family. Cell 1982, 29: 995-1004. 10.1016/0092-8674(82)90463-9View ArticlePubMedGoogle Scholar
- Engels WR: Hybrid dysgenesis in Drosophila melanogaster - rules of inheritance of female sterility. Genet Res 1979, 33: 219-236. 10.1017/S001667230001836XView ArticleGoogle Scholar
- Kidwell MG, Novy JB: Hybrid dysgenesis in Drosophila melanogaster - sterility resulting from gonadal-dysgenesis in the P-M system. Genetics 1979, 92: 1127-1140.PubMed CentralPubMedGoogle Scholar
- Khurana JS, Wang J, Xu J, Koppetsch BS, Thomson TC, Nowosielska A, Li C, Zamore PD, Weng Z, Theurkauf WE: Adaptation to P element transposon invasion in Drosophila melanogaster . Cell 2011, 147: 1551-1563. 10.1016/j.cell.2011.11.042PubMed CentralView ArticlePubMedGoogle Scholar
- Margulies L: A high level of hybrid dysgenesis in Drosophila : high thermosensitivity, dependence on DNA repair, and incomplete cytotype regulation. Mol Gen Genet 1990, 220: 448-455. 10.1007/BF00391752View ArticlePubMedGoogle Scholar
- Engels WR: In Mobile DNA. Edited by Berg DE, Howe MH. Washington, DC: American Society for Microbiology; 1989.Google Scholar
- Engels WR, Preston CR: Identifying P factors in Drosophila by means of chromosome breakage hotspots. Cell 1981, 26: 421-428. 10.1016/0092-8674(81)90211-7View ArticlePubMedGoogle Scholar
- Bucheton A, Paro R, Sang HM, Pelisson A, Finnegan DJ: The molecular basis of the I-R hybrid dysgenesis syndrome in Drosophila melanogaster - identification, cloning and properties of the I-factor. Cell 1984, 38: 153-163. 10.1016/0092-8674(84)90536-1View ArticlePubMedGoogle Scholar
- Yannopoulos G, Stamatis N, Monastirioti M, Hatzopoulos P, Louis C: hobo is responsible for the induction of hybrid dysgenesis by strains of Drosophila melanogater bearing the male recombination factor 23.5MRF. Cell 1987, 49: 487-495. 10.1016/0092-8674(87)90451-XView ArticlePubMedGoogle Scholar
- Evgenev MB, Zelentsova H, Shostak N, Kozitsina M, Barskyi V, Lankenau DH, Corces VG: Penelope , a new family of transposable elements and its possible role in hybrid dysgenesis in Drosophila virilis . Proc Natl Acad Sci U S A 1997, 94: 196-201. 10.1073/pnas.94.1.196View ArticleGoogle Scholar
- Lozovskaya ER, Scheinker VS, Evgenev MB: A hybrid dysgenesis syndrome in Drosophila virilis . Genetics 1990, 126: 619-623.PubMed CentralPubMedGoogle Scholar
- Scheinker VS, Lozovskaya ER, Bishop JG, Corces VG, Evgenev MB: A long terminal repeat-containing retrotransposon in mobilized during hybrid dysgenesis in Drosophila virilis . Proc Natl Acad Sci U S A 1990, 87: 9615-9619. 10.1073/pnas.87.24.9615PubMed CentralView ArticlePubMedGoogle Scholar
- Rozhkov NV, Aravin AA, Zelentsova ES, Schostak NG, Sachidanandam R, McCombie WR, Hannon GJ, Evgen’ev MB: Small RNA-based silencing strategies for transposons in the process of invading Drosophila species. RNA 2010, 16: 1634-1645. 10.1261/rna.2217810PubMed CentralView ArticlePubMedGoogle Scholar
- Petrov DA, Schutzman JL, Hartl DL, Lozovskaya ER: Diverse transposable elements are mobilized in hybrid dysgenesis in Drosophila virilis . Proc Natl Acad Sci U S A 1995, 92: 8050-8054. 10.1073/pnas.92.17.8050PubMed CentralView ArticlePubMedGoogle Scholar
- Blumenstiel JP, Hartl DL: Evidence for maternally transmitted small interfering RNA in the repression of transposition in Drosophila virilis . Proc Natl Acad Sci U S A 2005, 102: 15965-15970. 10.1073/pnas.0508192102PubMed CentralView ArticlePubMedGoogle Scholar
- Vieira J, Vieira CP, Hartl DL, Lozovskaya ER: Factors contributing to the hybrid dysgenesis syndrome in Drosophila virilis . Genet Res 1998, 71: 109-117. 10.1017/S001667239800322XView ArticlePubMedGoogle Scholar
- Evgen’ev MB: What happens when Penelope comes? An unusual retroelement invades a host species genome exploring different strategies. Mobile Genet Elem 2013, 3: e24542. 10.4161/mge.24542View ArticleGoogle Scholar
- Pyatkov KI, Shostak NG, Zelentsova ES, Lyozin GT, Melekhin MI, Finnegan DJ, Kidwell MG, Evgen’ev MB: Penelope retroelements from Drosophila virilis are active after transformation of drosophila melanogaster . Proc Natl Acad Sci U S A 2002, 99: 16150-16155. 10.1073/pnas.252641799PubMed CentralView ArticlePubMedGoogle Scholar
- Rozhkov NV, Schostak NG, Zelentsova ES, Yushenova IA, Zatsepina OG, Evgen’ev MB: Evolution and dynamics of small RNA response to a retroelement invasion in Drosophila . Mol Biol Evol 2013, 30: 397-408. 10.1093/molbev/mss241View ArticlePubMedGoogle Scholar
- Clark AG, Eisen MB, Smith DR, Bergman CM, Oliver B, Markow TA, Kaufman TC, Kellis M, Gelbart W, Iyer VN, Pollard DA, Sackton TB, Larracuente AM, Singh ND, Abad JP, Abt DN, Adryan B, Aguade M, Akashi H, Anderson WW, Aguadro CF, Ardell DH, Arguello R, Artieri CG, Barbash DA, Barker D, Barsanti P, Batterham P, Batzoglu S, Drosphila 12 Genomes Consortium, et al.: Evolution of genes and genomes on the Drosophila phylogeny. Nature 2007, 450: 203-218. 10.1038/nature06341View ArticlePubMedGoogle Scholar
- Bosco G, Campbell P, Leiva-Neto JT, Markow TA: Analysis of Drosophila species genome size and satellite DNA content reveals significant differences among strains as well as between species. Genetics 2007, 177: 1277-1290. 10.1534/genetics.107.075069PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics 2010, 26: 589-595. 10.1093/bioinformatics/btp698PubMed CentralView ArticlePubMedGoogle Scholar
- Smith CD, Edgar RC, Yandell MD, Smith DR, Celniker SE, Myers EW, Karpen GH: Improved repeat identification and masking in dipterans. Gene 2007, 389: 1-9. 10.1016/j.gene.2006.09.011PubMed CentralView ArticlePubMedGoogle Scholar
- Merriman PJ, Grimes CD, Ambroziak J, Hackett DA, Skinner P, Simmons MJ: S elements: a family of Tc1-like transposons in the genome of Drosophila melanogaster . Genetics 1995, 141: 1425-1438.PubMed CentralPubMedGoogle Scholar
- Lyozin GT, Makarova KS, Velikodvorskaja VV, Zelentsova HS, Khechumian RR, Kidwell MG, Koonin EV, Evgen’ev MB: The structure and evolution of penelope in the virilis species group of Drosophila : an ancient lineage of retroelements. J Mol Evol 2001, 52: 445-456.PubMedGoogle Scholar
- Morales-Hojas R, Vieira CP, Vieira J: The evolutionary history of the transposable element penelope in the Drosophila virilis group of species. J Mol Evol 2006, 63: 262-273. 10.1007/s00239-005-0213-1View ArticlePubMedGoogle Scholar
- Zwickl D: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. Austin, TX: The University of Texas; 2006.Google Scholar
- Li HP, Stephan W: Inferring the demographic history and rate of adaptive substitution in Drosophila . PLoS Genet 2006, 2: 1580-1589.Google Scholar
- Daniels SB, Peterson KR, Strausbaugh LD, Kidwell MG, Chovnick A: Evidence for horizontal transmission of the P transposable element between Drosophila species. Genetics 1990, 124: 339-355.PubMed CentralPubMedGoogle Scholar
- Bucheton A, Simonelig M, Vaury C, Crozatier M: Sequences similar to the I transposable element involved in I-R hybrid dysgenesis in D. melanogaster occur in other Drosophila species. Nature 1986, 322: 650-652. 10.1038/322650a0View ArticleGoogle Scholar
- Boussy IA, Itoh M: Wanderings of hobo : a transposon in Drosophila melanogaster and its close relatives. Genetica 2004, 120: 125-136.View ArticlePubMedGoogle Scholar
- Sokolova MI, Zelentsova ES, Shostak NG, Rozhkov NV, Evgen’ev MB: Ontogenetic consequences of dysgenic crosses in Drosophila virilis . Int J Dev Biol 2013, 57: 731-739. 10.1387/ijdb.120189meView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.