Birth, School, Work, Death, and Resurrection: The Life Stages and Dynamics of Transposable Element Proliferation

Transcript

1 G C A T T A C G genes G C A T Review Birth, School, Work, Death, and Resurrection: The Life Stages and Dynamics of Transposable Element Proliferation Justin P. Blumenstiel Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS 66049, USA; [email protected]    Received: 21 March 2019; Accepted: 23 April 2019; Published: 3 May 2019 Abstract: Transposable elements (TEs) can be maintained in sexually reproducing species even if they are harmful. However, the evolutionary strategies that TEs employ during proliferation can modulate ff erent life stages of a TE lineage, from birth to proliferation their impact. In this review, I outline the di to extinction. Through their interactions with the host, TEs can exploit diverse strategies that range from long-term coexistence to recurrent movement across species boundaries by horizontal transfer. TEs can also engage in a poorly understood phenomenon of TE resurrection, where TE lineages can apparently go extinct, only to proliferate again. By determining how this is possible, we may obtain new insights into the evolutionary dynamics of TEs and how they shape the genomes of their hosts. transposable element; horizontal transfer; arms race; LINE-1; Alu Keywords: hobo ; I element ; 1. Introduction “And he that was dead came forth, bound hand and foot with graveclothes.” John 11:44. Transposable elements (TEs) have an intimate relationship with the genomes of their hosts. Like any form of parasite they cause harm but they are also dependent on the host for fitness. However, unlike typical parasites, they are directly embedded in the genomes of their hosts. How can such parasites spread if they are harmful? Alleles that are harmful are expected to be lost, but transposable elements exist in essentially all forms of life. In eukaryotes, the persistence of TEs is explained by the fact that sexual reproduction allows TEs to spread even if their net e ff ect is a reduction in host fitness. Gamete fusion allows TEs to colonize new genomes [ ] and recombination breaks up the 1 association between progenitor copies and harmful descendant copies [ 2 , 3 ]. However, if TEs proliferate too rapidly within genomes, the consequences of their harm can indeed become too high and impede their success [ 4 ]. Transposable elements must walk a fine line between a su ffi cient rate of proliferation and one that is not so great that TEs become too burdened by the harmful e ff ects that they impose. The nature of this tension depends on the degree of intimacy with the host genome and is illuminated by considering the moment when a TE and the host genome first meet. This occurs during 5 – 9 ]). horizontal transfer, which is the first stage in the life cycle of a TE (see reviews on TE life cycles [ When a TE first invades a genome, it is a particularly fragile moment for the TE family because such events are likely to be serendipitous. For an element to be successful during the early stages of invasion, it must exploit these chance moments and avoid being lost from the population by drift. Studies show that the optimal TE strategy during horizontal transfer is to have a very high initial transposition rate [ 4 , 10 ]. This arises from the fact that the probability that a new TE becomes established is similar to the probability of fixation for a new beneficial allele. In the case of a new beneficial allele, the probability of fixation is ~2 s , where s is the beneficial selection coe ffi cient. For a transposon, the probability is the transposition rate and of establishment is ~2( s ), where u − s is the selection coe ffi cient that u genes 2019 , 10 , 336; doi:10.3390 / genes10050336 www.mdpi.com / journal / Genes

2 Genes 2019 10 , 336 2 of 14 , ff ect of each new single insertion [ , 10 , 11 ]. Establishment is achieved measures the average harmful e 4 when, on average, each individual in the population has one copy. Rather than fixation, I consider establishment to be a more appropriate term for TE families because fixation is a term that is more appropriate for alleles. Transposable elements insertions within the population are non-allelic if they reside at di ff erent locations in the genome. So, if each individual on average carries one insertion, the TE family can be considered established. A single TE insertion allele can be considered fixed if there are no non-insertion variants segregating in the population at that locus. For both a new beneficial allele and a new transposable element, the fixation (or establishment) probabilities do not depend much on the population size since the dynamics of stochastic loss by drift when the novel variant first appears are the same whether the population size is one million or one trillion. However, a transposition rate that is too high, while it will increase the probability that a TE becomes established, may also impose such a burden that the host may become extinct if the selection regime fails to limit the ever-increasing copy number. Thus, it has been shown that the optimal strategy for a transposable element is to have a high transposition rate during early invasion, 4 ]. This lower rate of movement may followed afterwards by a period with a lower rate of movement [ be enabled by host TE suppression mechanisms such as small RNA silencing. It is not apparent that selection on TE lineages would be e ffi cient enough to directly select such a tunable strategy. However, this tension reveals that optimal TE strategies will depend on the nature of the relationship with a genome. On one end of the continuum, TEs may be long term residents. On the other end, TEs may adopt a strategy of rapid invasion and movement from species to species. In the first part of this review, I discuss the nature and implications of these two strategies. Then, I consider an interesting phenomenon of TE lineages that appear to reside within genomes, go extinct, and then apparently come “back to life” many generations later. I will argue that TEs that show this pattern—I will designate them Lazarus elements—may highlight interesting aspects of TE biology and host interaction. 2. Long-Lasting Relationships Some TE lineages are long-lived residents of their host genomes. In some cases, this is because TEs have adopted a cooperative strategy with the host. For example, in Drosophila , telomere function has been assumed by TEs [ 12 ]. However, for TEs that remain parasitic with respect to the host, there may be no better example of long-term coexistence than the LINE-1 elements of mammalian genomes. LINE-1 elements are a member of the non-LTR retrotransposon class and have been residents of mammalian 13 – 15 ]. In humans, the LINE-1 element has had a genomes since early in the radiation of mammals [ 16 , profound role in shaping the genome and there are approximately 500,000 copies of this element [ ]. 17 The LINE-1 family is shared across most mammals due to continued vertical transmission since early in the mammalian radiation [ 18 , 19 ]. Vertebrates that include reptiles, amphibians and fish also share LINE-1 elements, suggesting that the LINE-1 element may have been present since before the origin of mammals [ 14 ]. Alternately, it has been proposed that LINE-1 elements entered the therian mammal ancestor (rather than the ancestor of all mammals) through horizontal transfer. This is suggested by the observation that monotremes lack LINE-1 elements and have no clear signature of their previous activity [ 20 ]. In either case, the LINE-1 lineage shows a striking level of persistence and success across mammals through ongoing vertical transmission. What has enabled this intimate relationship for millions of years within mammals? Phylogenetic analysis of LINE-1 elements within mammals has revealed a particular feature of LINE-1 persistence. In particular, phylogenetic trees of LINE-1 elements within a genome frequently have a “ladder-like” appearance [ 14 , 21 , 22 ]. This represents a scenario in which, through evolutionary time, there is typically only one or few proliferating lineages. This phylogenetic pattern has been proposed to be driven by an ongoing evolutionary arms-race with the host [ 14 ]. In particular, as mechanisms of LINE-1 control evolve on the part of the host, evolutionary innovation on the part of the TE lineage enables escape from host control. Recurrent cycles of adaptation and innovation—in both host and TE—can

3 Genes 2019 10 , 336 3 of 14 , 21 ]. This pattern may also be driven thus lead to the persistence of a single successful TE lineage [ ective population sizes that are likely more common in mammals. In very large by the smaller e ff populations, the fixation of an active and harmful TE insertion allele by drift is unlikely. However, in smaller populations, drift may allow such insertion alleles to fix. When an active copy becomes fixed at a particular locus, only decay into a non-functional state will allow the active copy to be lost from the population. Thus, fixation of an active TE insertion allele represents a critical stage in TE-host dynamics. Faced with the continued presence of the LINE-1 element over millions of years, specialized modes of host control are proposed to contribute to the evolutionary dynamics that yield the arms-race driven “ladder” phylogeny. In particular, new active LINE-1 lineages may carry key innovations that enable specialized modes of escape from repression [ 23 ]. Diverse proteins that restrict LINE-1 24 ]. Signatures of recurrent transposition include APOBEC3, MOV10, ZAP, SAMHD1 and ZNF93 [ LINE-1 adaptation that allow evasion from these restricting factors have also been found. For example, within mammals, the 5’ UTR of LINE-1 is highly dynamic [ , 25 – 27 ]. This has been proposed to be 22 driven by the ongoing evolution of KRAB zinc fingers that can evolve specificity to target particular sequences in LINE-1 for repression. In response to this, it appears that selection on the LINE-1 lineage has driven removal of particular target sequences from the 5’ UTR [28]. The ongoing persistence of one or few evolving LINE-1 lineages is likely enforced by within ff erent modes of adaptation to evolve on distinct TE lineage competition. Otherwise, we might expect di lineages, followed by successful diversification. Competition for host factors required for transposition 29 ]. Strikingly, and in contrast to mammalian systems, has been proposed to contribute to this dynamic [ the proliferation of one or few LINE-1 element lineages does not seem to apply in other vertebrates [ 14 ]. Rather, multiple lineages of LINE-1 elements have expanded and proliferated in the genomes of reptiles, amphibians and fish [ 30 – 32 ]. This represents a distinct mode of long-term coexistence within the genomes of non-mammalian species. Di ff erences in demographic history and the strength of selection ff are likely to contribute to this di erence. Compared to mammals, some non-mammalian species with greater LINE-1 diversity also show a stronger signature of selection acting to limit the fixation of TE 33 ]. This suggests that di ff erent selection regimes may contribute to the di ff erence in insertion alleles [ 34 ]). One di ff erence LINE-1 dynamics between mammalian and non-mammalian species (but see [ may arise from di ff erences in the probability of ectopic recombination between dispersed repeats [ 29 ]. Selection against ectopic recombination is an important determinant of TE dynamics and a low rate across mammalian genomes may decrease the strength of selection against insertions and allow the accumulation of repetitious sequences [ – 40 ]. In addition, if lower levels of ectopic recombination 35 allow greater TE accumulation, persisting copies that fix by drift may intensify competition for host factors. Thus, as genomic copy number increases due to reduced levels of genome-wide ectopic recombination, the magnitude of competition for host factors may increase among competing copies and lineages. This may lead to a greater tendency for a single lineage to outcompete all other lineages. For these reasons, selection on LINE-1 lineages may not simply be to evade host restriction factors. Selection to increase access to host factors that enable transposition, amidst a genome filled with many other copies, may also be critical. 3. Horizontal Transfer: Fast, Cheap and Out of Control “Based on our experience in building ground based mobile robots (legged and wheeled), we argue here for fast, cheap missions using large numbers of mass produced simple autonomous robots...” Brooks and Flynn. 1989. Fast, Cheap and Out of Control: A robot invasion of the solar system. These contrasting modes of LINE-1 evolution—the proliferation of a few lineages in mammals vs. diversification in reptiles, amphibians and fish—represent two forms of long-term co-existence. As previously indicated, long-term co-existence can also be maintained if TEs adapt a strategy of cooperation, as seen in the case of Drosophila telomeres. However, for selfish TEs that display parasitic

4 Genes 2019 10 , 336 4 of 14 , behavior with respect to the host, another strategy relies on horizontal transfer and recurrent invasion. If TEs have the capacity to invade genomes through horizontal transfer, long-term persistence may , 42 ]. If a TE family can invade a species, proliferate, be enabled by a ‘live fast, die young’ strategy [ 41 and jump to a new species, it may conceivably persist even if it is unlikely to endure within any mariner Drosophila illustrate how such a strategy is single species. Studies of the DNA transposon in – 45 possible [ mariner was discovered in D. mauritiana , a close relative of D. melanogaster . However, 43 ]. D. melanogaster species subgroup is considered “spotty” [ 46 , 47 ]. In particular, its presence within the D. melanogaster it appears in several close relatives of D. melanogaster itself. It has but is absent from mariner apparently been lost. Interestingly, an additional lineage is also found in the genomes of other members of the melanogaster species subgroup, including D. erecta , but was apparently lost from the / D. simulans clade [ 48 ]. This latter mariner D. melanogaster family also shares 97% sequence similarity with a element found in the cat flea, indicating horizontal transfer several million years ago. mariner mariner dynamics can be explained by a dynamic process of Overall, these patterns indicate that 48 ]. In contrast to mammals, it appears that horizontal recurrent horizontal transfer and extinction [ transfer is rampant in insect species. In a comprehensive analysis of the genomes of nearly 200 insect species, more than 2000 horizontal transfer events were found to have occurred within a span of about 10 million years [ 49 ]. Strikingly, the Tc1 / mariner class of DNA transposons shows the greatest frequency of horizontal transfer. This high propensity for horizontal transfer has been attributed to a lack of dependence on host factors for transposition [ ]. Tc1 / mariner cis regulatory sequences that drive 50 ffi 51 ]. Within transcription in diverse genomes may also facilitate e cient movement across species [ a single species, a TE lineage can proliferate if its transposition rate is su ffi ciently high so that it can increase at a rate faster than its removal due to negative selection. The same principle should also apply across species. If a TE can invade, by horizontal transfer, the genomes of new species at a rate faster than the within species extinction rate, the lineage will also find success. In this case, since TE success ffi cient for depends on being able to move across species, it is unlikely that natural selection will be su adaptation, on the part of a TE lineage, to a particular host genome. Rather, natural selection will favor a “generalist” strategy that enables movement in the genomes of many species. 4. Extinction Whether a TE is adapted for continued vertical transmission (as observed for LINE-1 elements) or ongoing movement across species (as perhaps observed for mariner elements), TE lineages are not guaranteed perpetual success. Rather, they can also go extinct within a species. Across mammals, LINE-1 extinction has been observed in the rhinoceros and lineages of rodents, bats, insectivores and Afrotherians [ 15 , 52 – 55 ]. Several mechanisms have been proposed to contribute to LINE-1 extinction. In one scenario, mechanisms of host suppression may be su ffi cient. It has been noted that the fate of a TE lineage depends on the balance between transposition and the rate of accumulation for degenerating mutations [ ]. If the transposition rate is lower than the rate of mutation that renders an element 56 inactive, then the TE lineage will decay. Thus, host control mechanisms that drive a significantly low transposition rate may also drive extinction by decay. Other factors are also likely to contribute to extinction. TE families may drive other TE families to extinction through direct competition for host factors. For example, LINE-1 extinction in a group of sigmodontine rodents may have been influenced by competition for host factors with an expanding endogenous retrovirus lineage [ 57 , 58 ]. Extinction may also be driven by other TE lineages through direct sequestration of TE-encoded factors that enable transposition. SINE elements, such as the Alu element, hijack LINE-1 encoded factors to favor their own increase [ ]. Thus, Alu amplification may 59 drive LINE-1 extinction through competitive saturation of LINE-1 encoded factors required for LINE-1 transposition [ 60 ]. Finally, extinction may also be enabled by a form of lineage “suicide”. In the case of DNA transposons, internally deleted copies may titrate functional transposase from fully functional copies [ 61 ]. As internally deleted copies within the genome increase, active DNA transposon lineages may lose su ffi cient access to their own encoded factors.

5 Genes 2019 , 336 5 of 14 , 10 Finally, stochastic loss and demographic factors may also contribute to lineage extinction. In populations where an active TE does not fix at any particular location, selection or drift may simply lead to the loss of every active copy in the genome. This will be most likely when the ffi ciently low, so it is likely to be enhanced by host suppression mechanisms. transposition rate is su The dynamics of stochastic loss, in many ways, are likely to be similar to loss by the mutational degeneration of active copies. How long will it take for a TE family to be lost by this mechanism? Using simulation, I have shown that total copy number within the population—rather than population size or per genome copy number—dominates the dynamics of stochastic loss assuming no individual insertion becomes fixed (Figure 1). Selection also plays a role. N = 10,000, s = 5/20,000 A 25,000 20,000 15,000 10,000 5,000 0 5 10 15 20 0 Median Number of Generations Until Loss Number of Copies Per Individual B Copies in Population = 50,000, s = 5/20,000 16,000 14,000 12,000 10,000 8,000 6,000 4,000 2,000 0 100,000 0 20,000 40,000 60,000 80,000 Median Number of Generations Until Loss Population Size C N = 10,000, Copies in Population = 50,000 35,000 30,000 25,000 20,000 15,000 10,000 5,000 0 10 8 0 2 4 6 Median Number of Generations Until Loss Selective effect = i/20,000 Figure 1. Dynamics of stochastic loss. N is the diploid population size and s is the selection coe ffi cient acting against single insertions. All simulations were performed by simple binomial sampling of insertion

6 Genes 2019 10 , 336 6 of 14 , / 2 . Sampling was iterated according to frequency in the population alleles starting at frequency of 1 N for a given number of copies. This procedure implicitly assumes there is no linkage. In addition, by assuming no actual transposition or degradation specifically, it is suitable to a scenario where the rate of transposition is equal to the rate of mutation to a non-functional state. Selection was simulated by A ffi ) Fixed population size adjusting the probability of sampling according to the selection coe cient. ( and negative selection coe ffi cient. The time until loss increases with per individual copy number. Note ) A fixed number of copies in the population, distributed among B that the rate of increase declines. ( erent population sizes. The time until loss is not a ff ected by population size. ( individuals of di ) An ff C ffi increasing selection coe cient, as expected, decreases the time until loss. 5. Resurrection Overall, the canonical life-cycle of a TE family starts with invasion followed by proliferation and eventual extinction. The duration for each of these stages may vary and be influenced by a wide variety of factors, as outlined previously. Extinction is certainly not guaranteed but there are many examples of where this appears to be the case. More striking, however, is that in some cases, extinction seems to be followed by resurrection. This is a mysterious phase of TE dynamics and worthy of investigation because it may shed light on the evolution of TE life-strategies that range between recurrent invasion and long-term coevolution. Resurrection, also known as re-invasion, occurs when an active TE lineage becomes quiescent and perhaps even extinct, and then later proliferates. Syndromes of hybrid dysgenesis were the first I-R syndrome of dysgenesis. Hybrid dysgenesis is a to reveal this phenomenon, in particular the syndrome of intraspecific sterility that occurs when active TE families transmitted paternally are absent or nearly absent from the maternal genome [ – 64 ]. In the absence of abundant maternal copies, a pool 62 65 , ]. This leads to activation of piRNAs that maintain TE repression is not provisioned to the zygote [ 66 of paternally inherited TEs and sterility. Perhaps the most well understood syndrome of hybrid P-M system. dysgenesis is the dysgenesis occurs when P elements inherited from P strain males, P-M M strain females, cause germline cell death [ 67 mated with 68 ] due to excessive transposition in the , absence of maternal element piRNAs. In the P-M system, the asymmetry in the P element abundance P P and M strains can be explained by recent horizontal transfer rather than resurrection [ 69 ]. between M laboratory strains devoid of P elements were established in the early part of the 20th century. P element invasion of natural populations via horizontal transfer occurred at a similar time, so natural P elements [ , 71 ]. In contrast, I-R dysgenesis seems to have arisen from 70 populations now carry many elements. I-R dysgenesis—observed as hatch failure in eggs laid by F1 females—occurs I resurrected I (inducer) strain males, carrying abundant non-LTR I retrotransposons, mate with R (reactive) when 63 strain females that lack active copies [ R ]. However, in contrast to the P-M system, the genomes of I ]. 72 – 75 strains are littered with degraded elements that are the fossils of a previous proliferation event [ In fact, under certain conditions, the degraded elements can contribute to the piRNA pool and I mediate repression of the newer I elements [ 76 ]. Thus, the genome retains a memory of past invasion and still retains some capacity to restrain new I elements. Two additional cases of hybrid dysgenesis reveal a similar scenario, indicating that resurrection may be a common but poorly appreciated part of the life cycle of TEs. In a third Drosophila melanogaster, hobo hobo is a DNA transposon that causes hybrid case of hybrid dysgenesis is driven by the element. hobos are mated with females that lack them [ dysgenesis when males carrying multiple active ]. 64 Studies indicate that American populations lacked active hobo elements in the 1950s [ 77 ]. Strikingly, the genomes of D. melanogaster as well as close relatives all carry degraded copies of hobo . This suggests that an active version of the element was present in an ancestor of all D. melanogaster , was lost, hobo but now has proliferated to the extent that it can cause hybrid dysgenesis. This new hobo variant also appears among close relatives that include D. simulans [ 78 – 81 ]. Finally, a similar scenario is observed in the hybrid dysgenesis syndrome of Drosophila virilis . The Penelope element likely contributes to this

7 Genes , 7 of 14 2019 10 , 336 syndrome and it also represents a case of re-invasion. Multiple degraded copies within the genome of 82 ]. D. virilis , and relatives suggest a scenario of proliferation and extinction, followed by re-invasion [ 83 As phylogenetic analysis shows (Figures 2 and 3), the mode of past and current proliferation hobo erence cases. For the varies across these di ff element, it appears proliferation occurred millions of clade. This corresponds to the years ago, prior to the divergence of D. melanogaster and the D. simulans , D. simulans and D. sechellia are phylogeny where lineages of hobo upper portion of the D. melanogaster intermingled in a complex manner. This was followed by another proliferation, perhaps only in the clade. As active and nearly identical copies currently exist in D. sechellia , D. simulans / D. melanogaster , a recent wave of re-activation appears to have occurred across all three D. sechellia D. simulans and D. melanogaster . Phylogenetic analysis suggests that the currently active species, but especially in lineage was derived from a clade of proliferating elements derived solely from D. simulans , as has been proposed [ , the previous proliferation of the I element seems more ongoing, 80 ]. In contrast to hobo genome. However, it also appears D. melanogaster but with one large wave of historical activity in the element resides in D. simulans . This supports the possibility that that the closest relative of the I and the element were introduced into D. melanogaster from both currently active variants of I hobo D. simulans . While these species do not readily produce fertile hybrids, they can in some cases. sim sech mel mel sim mel mel mel mel mel sech sim sech sech mel mel sim sech mel mel mel mel mel mel mel mel sech sim sim sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech sech hobo sech sech sech sech sech sech sech sech sim sim sech sim sim sech sim sim sim sech sech sech sim sech sech sech sech sech sim sech sim mel mel mel mel mel mel sim mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel mel 0 . 0 3 Phylogenetic analysis of hobo element fragments extracted from the genomes of D. melanogaster , Figure 2. D. simulans and D. sechellia . Alignments from BLAST output (E cuto ff − 100) were subjected to phylogenetic analysis with GARLI. Blue indicates , red indicates D. sechellia and orange D. melanogaster indicates D. simulans . The hobo phylogeny shows evidence of a previous wave in D. melanogaster , as well as new proliferation derived from a lineage residing within D. simulans or D. sechellia .

8 Genes 8 of 14 , 336 2019 , 10 sech mel mel sech mel mel sech mel mel mel sim mel mel sim mel mel mel sim mel sech sim sech mel mel sim sim sech sim sim sim sim sim sim sim sech mel sech sim sech sech sech sim sim sim sim sim mel mel sech mel mel sim sech mel sim sech sech sech sim sim sech mel sim mel mel sech mel mel sech sim sech mel sech mel mel mel sim mel mel sech mel mel mel sim sech sech sim sech mel sim sim mel sim sech sech sim sech sim sech sech sim mel mel mel sim sim sech sim sech sim sim sech sech sim sech sech sech sech sech sim sim sech sim sim sim sim sech sim sech sim sech sech sech sech sim sim sim sim sech sech mel mel sech mel sech sim sech sim sim sim sim sim sim sim sech sech sim sech sech sech sech sech sech mel mel sech sech sech sech sech sech sech sech sech sech sech sim sech sech sech sim sech sech sech sech sech sech sim sim mel sech mel mel mel sech sech mel sech mel mel mel sech mel mel sim mel mel mel sech mel mel sech mel mel mel sech mel mel sech mel mel sech mel mel mel Prior wave in sech mel mel sech mel mel mel sim mel mel sech mel D. melanogaster mel mel sech mel mel sech mel mel mel sech mel mel sech mel mel mel sech mel mel sech mel mel mel sech mel mel sech mel mel mel sech mel mel sech mel mel sech mel mel mel sech mel mel sech mel sech sim sech mel sech sech sim sim sech sech sech sech sech sim sim sim sim sim sim sim sim sim sech sim sech sim sim sim sech sim sim sech sim sim sim sech sim sech sech sim sech sech sim sim sech sim sech sech sim sim sim sim sim sim sech sech sech sech sech sim sech sech sech sim sech sech sech sim sech sim sim sech sim sech sech sech sech sech sech sech sim sim sech sech sim sech sech sech sech sech sech sim sech sech sech sech sech I element sech sech sech sech sech sim sech sech sech sech sech sech sech sech sech sech sim sech sech sim sech sech sech sech sech sech sech sech sech sech sim sech sim sim sim sim sim sech sim sim sech sech sech sech sim sech sim sech sech sech sech sech sech sech sech sech sech sech mel sech sech sech sech mel sech sech sech mel sech sech sim sim sech sech sech sim mel sim sech sech sim sech mel sech sech mel mel mel sim mel sech mel mel mel Current wave in mel sim mel mel 0 . 0 3 mel mel mel mel mel sech D. melanogaster sech sech mel mel mel sech mel mel mel 0 . 0 3 Figure 3. Phylogenetic analysis of element fragments extracted from the genomes of D. melanogaster , I D. simulans D. sechellia . Alignments from BLAST output (E cuto ff − 100) were subjected to and phylogenetic analysis with GARLI. Blue indicates D. melanogaster , red indicates D. sechellia and orange indicates . A zoom-in is provided for further detail of recent dynamics. Like hobo , D. simulans the I element shows evidence of a previous wave in D. melanogaster , as well as new proliferation derived from a lineage residing within D. simulans or D. sechellia .

9 Genes , , 336 9 of 14 2019 10 9 of 14 , x FOR PEER REVIEW 10 , Genes 2019 How does re-invasion occur? What does the time between waves of proliferation tell us about How does re-invasion occur? What does the time between waves of proliferation tell us about the likely mechanism of re-invasion? This is a critical question and I consider two possibilities. One the likely mechanism of re-invasion? This is a critical question and I consider two possibilities. possible explanation is that re-invasion arises through iterated rounds of invasion by horizontal One possible explanation is that re-invasion arises through iterated rounds of invasion by horizontal ding in a reservoir species (of any kind, not just a transfer (Figure 4A). Specifically, a TE lineage resi transfer (Figure 4A). Specifically, a TE lineage residing in a reservoir species (of any kind, not just close relative), jumps into the host. After the firs t round of horizontal transfer, proliferation is a close relative), jumps into the host. After the first round of horizontal transfer, proliferation is transfer event is followed by an additional round followed by extinction. In turn, a second horizonal followed by extinction. In turn, a second horizonal transfer event is followed by an additional round of proliferation. The reservoir source species of th e original horizontal transfer event may not be the of proliferation. The reservoir source species of the original horizontal transfer event may not be the same as the second. Under this mechanism, the time between proliferation events may be significant. same as the second. Under this mechanism, the time between proliferation events may be significant. A B Figure 4. Dynamics of transposable element (TE) resurrection. ( A ) Iterated rounds of horizontal transfer (either through hybridization or other mechanisms) between sympatric close relatives allows extinction Figure 4. A ) Iterated rounds of horizontal Dynamics of transposable element (TE) resurrection. ( B ) Population subdivisions (indicated with events to be followed by new rounds of proliferation. ( transfer (either through hybridization or other me chanisms) between sympatric close relatives allows the blue dashed line) allow a TE lineage to decrease in abundance, go extinct from one population, extinction events to be followed B by new rounds of proliferation. ( ) Population subdivisions but become resurrected through contact with an isolated refuge population. (indicated with the blue dashed line) allow a TE li neage to decrease in abundance, go extinct from A second possibility is that a TE lineage becomes quiescent, but active copies linger within the me resurrected through one population, but beco ed refuge population. contact with an isolat species (Figure 4B). In this case, the TE lineage will have experienced only apparent, but not complete A second possibility is that a TE lineage becomes quiescent, but active copies linger within the extinction. How could this occur? One scenario is that host suppression mechanisms first become ll have experienced only species (Figure 4B). In this case, the TE lineage wi apparent, but not complete su ffi cient so as to essentially cease transposition. This may occur through the fixation of host suppressor that host suppression mechanisms first become extinction. How could this occur? One scenario is alleles. For example, a TE copy may land in a piRNA cluster and drive piRNA suppression that is sufficient so as to essentially cease transposition. This may occur through the fixation of host strong enough to drive an extremely low transposition rate. Subsequently, if no TE insertion alleles suppressor alleles. For example, a TE copy may land in a piRNA cluster and drive piRNA are fixed within the population (due to selection against insertions), then all functional TE copies will suppression that is strong enough to drive an ex tremely low transposition rate. Subsequently, if no eventually be lost by a combination of drift, selection and degeneration. However, this may take a TE insertion alleles are fixed within the population (due to selection against insertions), then all very long period of time. During this lingering phase, several outcomes are possible. First, selection to functional TE copies will eventually be lost by a combination of drift, selection and degeneration. retain host suppressor alleles (such as piRNA silencing alleles) may become reduced because there ering phase, several outcomes are period of time. During this ling However, this may take a very long are few TE copies that remain a threat. However, if lingering TE copies take advantage of degraded possible. First, selection to retain host suppres sor alleles (such as piRNA silencing alleles) may suppressor alleles that have become frequent within the population, the TE lineage may return to become reduced because there are few TE copies that remain a threat. However, if lingering TE copies an active state and in turn prevent the fixation of degraded host suppressor alleles. Depending on take advantage of degraded suppressor alleles that have become frequent within the population, the the timing of these dynamics, this may simply appear to represent a case of continued proliferation. rn prevent the fixation of degraded host suppressor TE lineage may return to an active state and in tu An alternate outcome may give the illusion of resurrection. In this case, host suppression alleles may alleles. Depending on the timing of these dynamics, this may simply appear to represent a case of become degraded and these degraded variants may increase in frequency and perhaps even fix. In the the illusion of resurrection. In this case, host continued proliferation. An alternate outcome may give scenario where non-functional host suppressor alleles are neutral, perhaps due to an exceedingly e degraded variants may increase in frequency suppression alleles may become degraded and thes -functional host suppressor alleles are neutral, and perhaps even fix. In the scenario where non

10 Genes 2019 10 , 336 10 of 14 , low abundance of lingering TE copies, the time until suppressor decay will be similar to the rate of mutation to the non-functional state. This may be of the order of millions of generations and it may be unlikely that an active lineage could persist for this time. But if lingering TE copies persist, the loss of a host suppressor allele may lead to a new round of TE proliferation. This scenario may be more likely if isolated populations can function as source refugia for new waves of TE proliferation. Depending on the timing of these dynamics, this may give the appearance of TE resurrection. Even if a host suppression allele is not lost, TE resurrection may also be apparent if a novel TE variant arises during the lingering phase. Such a variant may confer resistance to the host suppression allele, thus allowing a new round of TE proliferation. This would be analogous to evolutionary rescue of the TE lineage. Finally, true resurrection may be possible if recombination between non-functional copies leads to restoration of a functional copy. These two scenarios may explain the appearance of Lazarus elements. In the first case, iterated rounds of horizontal transfer give the appearance of resurrection. In the latter case, an active lineage becomes quiescent for a duration, but lingers. After some period of time, it becomes active again. The key distinction between these two scenarios is the source for the newly proliferating lineage. Is it “from without” or “from within”? The likelihood of these two explanations depends on a large number of unknown parameters. It is apparent that horizontal transfer is quite common for transposable elements and this fact lends support to the “from without” hypothesis. In fact, at least for , Drosophila it appears that iterated rounds of horizontal transfer among close relatives might play an important hobo role in the dynamics of re-invasion. For both the I element, copies shared among close and the relatives are highly similar. Thus, rare hybridization may contribute to continued TE exchange and enable iterated bouts of re-invasion. Horizontal transfer of TEs among close relatives of D. melanogaster appears rampant [ 84 ]. A similar pattern of ongoing exchange of TEs also appears in close relatives of D. pseudoobscura that are in sympatry [ ]. Overall, the genomes of closely related species in sympatry 85 ff may function as an e ective higher-level ecosystem for TEs. As modes of TE silencing in one species decay, horizontal transfer mediated by either hybridization or other mechanisms may allow iterated rounds of re-introduction. This dynamic is likely to be shaped by the dynamics of decay of host suppressor alleles. Nonetheless, it is also apparent that systems of host suppression are in rapid flux, and selection for host suppression in any single species may be quite weak. Thus, host alleles that suppress TE movement—such as TE insertions into piRNA clusters—may decay after the threat imposed by a TE family becomes reduced, even if the TE lineage has not yet become completely extinct within the species. In this case, new proliferation events may arise “from within”. How might we test between these hypotheses? The key may be to identify the source refugia for resurrected TEs. For the “from without” hypothesis, these may be other species that live in physical proximity, such as close relatives that can hybridize. Alternately, shared parasites that enable horizontal transfer may also function as refugia [ 86 ]. However, horizontal transfer can be a very rare event, so proof of source may be extremely challenging. For the “from within” hypothesis, this may require a closer study of TE dynamics in large populations that have geographic structure. Theoretical studies may also examine whether it is plausible for a functional TE lineage to persist at low frequency until a host suppressor allele is lost, followed by re-invasion. Examples of apparent resurrection, where divergence between active and degraded copies are in the 5% to 10% range, suggest that the time would perhaps be too long. However, in the case of the I element where degraded copies may be quite young, the “from within” hypothesis might be more plausible. Altogether, distinguishing between these two possibilities will provide insight into the evolutionary strategies that TEs employ to ensure their continued presence across diverse species. Funding: Funding was provided by NSF Award 1413532 and the University of Kansas. Acknowledgments: I would like to thank two anonymous reviewers for the very helpful suggestions. Conflicts of Interest: I declare no conflict of interest.

11 Genes 2019 10 , 336 11 of 14 , References 1. 1982 , 101 , 519–531. [PubMed] Hickey, D.A. Selfish DNA: A sexually-transmitted nuclear parasite. Genetics 2. Charlesworth, B.; Sniegowski, P.; Stephan, W. The evolutionary dynamics of repetitive DNA in eukaryotes. Nature , 215–220. [CrossRef] 1994 371 , Genet. Res. , 42 , Charlesworth, B.; Charlesworth, D. The population dynamics of transposable elements. 3. 1983 1–27. [CrossRef] 4. Le Rouzic, A.; Capy, P. The first steps of transposable elements invasion: Parasitic strategy vs. genetic drift. Genetics 169 , 1033–1043. [CrossRef] 2005 , Am. Nat. Song, M.J.; Schaack, S. Evolutionary Conflict between Mobile DNA and Host Genomes. , 192 , 5. 2018 263–273. [CrossRef] [PubMed] Venner, S.; Miele, V.; Terzian, C.; Biemont, C.; Daubin, V.; Feschotte, C.; Pontier, D. Ecological networks to 6. PLoS Biol. 2017 , 15 , e2001536. [CrossRef] [PubMed] unravel the routes to horizontal transposon transfers. Venner, S.; Feschotte, C.; Biemont, C. Dynamics of transposable elements: Towards a community ecology of 7. Trends Genet. TIG 2009 , 25 , 317–323. [CrossRef] [PubMed] the genome. Kidwell, M.G.; Lisch, D.R. Perspective: Transposable elements, parasitic DNA, and genome evolution. 8. , 55 , 1–24. [CrossRef] Evolution 2001 Hartl, D.L.; Lozovskaya, E.R.; Nurminsky, D.I.; Lohe, A.R. What restricts the activity of mariner-like 9. transposable elements. 1997 , 13 , 197–201. [CrossRef] Trends Genet. TIG Kaplan, N.; Darden, T.; Langley, C.H. Evolution and extinction of transposable elements in Mendelian 10. Genetics 1985 , populations. , 459–480. [PubMed] 109 11. Groth, S.B.; Blumenstiel, J.P. Horizontal Transfer Can Drive a Greater Transposable Element Load in Large J. Hered. , 108 , 36–44. [CrossRef] Populations. 2017 12. Casacuberta, E.; Pardue, M.L. Transposon telomeres are widely distributed in the Drosophila genus: TART elements in the virilis group. 2003 , 100 , 3363–3368. [CrossRef] Proc. Natl. Acad. Sci. USA 13. Platt, R.N., 2nd; Vandewege, M.W.; Ray, D.A. Mammalian transposable elements and their impacts on genome evolution. Chromosome Res. 2018 , 26 , 25–43. [CrossRef] [PubMed] 14. Boissinot, S.; Sookdeo, A. The Evolution of LINE-1 in Vertebrates. 2016 , 8 , 3485–3507. Genome Biol. Evol. [CrossRef] [PubMed] 15. Ivancevic, A.M.; Kortschak, R.D.; Bertozzi, T.; Adelson, D.L. LINEs between Species: Evolutionary Dynamics 2016 , 8 of LINE-1 Retrotransposons across the Eukaryotic Tree of Life. Genome Biol. Evol. , 3301–3322. [CrossRef] [PubMed] 16. Boissinot, S.; Chevret, P.; Furano, A.V. L1 (LINE-1) retrotransposon evolution and amplification in recent Mol. Biol. Evol. 2000 , 17 , 915–928. [CrossRef] [PubMed] human history. 17. Lander, E.S.; Linton, L.M.; Birren, B.; Nusbaum, C.; Zody, M.C.; Baldwin, J.; Devon, K.; Dewar, K.; Doyle, M.; Nature , 409 , 860–921. FitzHugh, W.; et al. Initial sequencing and analysis of the human genome. 2001 18. Burton, F.H.; Loeb, D.D.; Voliva, C.F.; Martin, S.L.; Edgell, M.H.; Hutchison, C.A., 3rd. Conservation throughout mammalia and extensive protein-encoding capacity of the highly repeated DNA long interspersed J. Mol. Biol. 1986 , 187 , 291–304. [CrossRef] sequence one. 19. Smit, A.F.; Toth, G.; Riggs, A.D.; Jurka, J. Ancestral, mammalian-wide subfamilies of LINE-1 repetitive sequences. J. Mol. Biol. 1995 , 246 , 401–417. [CrossRef] Ivancevic, A.M.; Kortschak, R.D.; Bertozzi, T.; Adelson, D.L. Horizontal transfer of BovB and L1 20. Genome Biol. , 19 , 85. [CrossRef] retrotransposons in eukaryotes. 2018 21. Furano, A.V. The biological properties and evolutionary dynamics of mammalian LINE-1 retrotransposons. 2000 , Prog. Nucleic Acid Res. Mol. Biol. , 255–294. [PubMed] 64 22. Khan, H.; Smit, A.; Boissinot, S. Molecular evolution and tempo of amplification of human LINE-1 retrotransposons since the origin of primates. Genome Res. 2006 , 16 , 78–87. [CrossRef] 23. Boissinot, S.; Furano, A.V. Adaptive evolution in LINE-1 retrotransposons. Mol. Biol. Evol. 2001 , 18 , 2186–2194. [CrossRef] [PubMed] Goodier, J.L. Restricting retrotransposons: A review. Mob. DNA 2016 , 7 , 16. [CrossRef] [PubMed] 24. 25. Sookdeo, A.; Hepp, C.M.; McClure, M.A.; Boissinot, S. Revisiting the evolution of mouse LINE-1 in the genomic era. Mob. DNA 2013 , 4 , 3. [CrossRef] [PubMed]

12 Genes 2019 10 , 336 12 of 14 , Hayward, B.E.; Zavanelli, M.; Furano, A.V. Recombination creates novel L1 (LINE-1) elements in Rattus 26. 1997 , 146 , 641–654. [PubMed] norvegicus. Genetics 27. Wincker, P.; Jubier-Maurin, V.; Roizes, G. Unrelated sequences at the 5’ end of mouse LINE-1 repeated , 15 , 8593–8606. [CrossRef] Nucleic Acids Res. elements define two distinct subfamilies. 1987 Jacobs, F.M.; Greenberg, D.; Nguyen, N.; Haeussler, M.; Ewing, A.D.; Katzman, S.; Paten, B.; Salama, S.R.; 28. / 93 and SVA / L1 Haussler, D. An evolutionary arms race between KRAB zinc-finger genes ZNF91 2014 , 516 retrotransposons. Nature , 242–245. [CrossRef] Furano, A.V.; Duvernell, D.D.; Boissinot, S. L1 (LINE-1) retrotransposon diversity di ers dramatically 29. ff 2004 , 9–14. [CrossRef] 20 Trends Genet. TIG between mammals and fish. , 30. Novick, P.A.; Basta, H.; Floumanhaft, M.; McClure, M.A.; Boissinot, S. The evolutionary dynamics of autonomous non-LTR retrotransposons in the lizard Anolis carolinensis shows more similarity to fish than Mol. Biol. Evol. 2009 , 26 , 1811–1822. [CrossRef] mammals. 31. Hellsten, U.; Harland, R.M.; Gilchrist, M.J.; Hendrix, D.; Jurka, J.; Kapitonov, V.; Ovcharenko, I.; Putnam, N.H.; Science 2010 328 , 633–636. Shu, S.; Taher, L.; et al. The genome of the Western clawed frog Xenopus tropicalis. , [CrossRef] Duvernell, D.D.; Pryor, S.R.; Adams, S.M. Teleost fish genomes contain a diverse array of L1 retrotransposon 32. 2004 , 59 , 298–308. [CrossRef] lineages that exhibit a low copy number and high rate of turnover. J. Mol. Evol. [PubMed] Tollis, M.; Boissinot, S. Lizards and LINEs: Selection and demography a ff ect the fate of L1 retrotransposons 33. Anolis carolinensis ). Genome Biol. Evol. 2013 in the genome of the green anole ( 5 , 1754–1768. [CrossRef] , 34. ff erential E ff ect of Selection against LINE Xue, A.T.; Ruggiero, R.P.; Hickerson, M.J.; Boissinot, S. Di Retrotransposons among Vertebrates Inferred from Whole-Genome Data and Demographic Modeling. Genome Biol. Evol. , 10 , 1265–1281. [CrossRef] 2018 Petrov, D.A.; Aminetzach, Y.T.; Davis, J.C.; Bensasson, D.; Hirsh, A.E. Size matters: Non-LTR retrotransposable 35. Mol. Biol. Evol. 2003 , 20 , 880–892. [CrossRef] [PubMed] elements and ectopic recombination in Drosophila. 36. Montgomery, E.A.; Huang, S.M.; Langley, C.H.; Judd, B.H. Chromosome rearrangement by ectopic recombination in Drosophila melanogaster Genetics 1991 , 129 , 1085–1098. : Genome structure and evolution. Lee, Y.C.; Langley, C.H. Long-term and short-term evolutionary impacts of transposable elements on 37. 2012 , 192 , 1411–1432. [CrossRef] [PubMed] Drosophila. Genetics Charlesworth, B.; Langley, C.H.; Sniegowski, P.D. Transposable element distributions in Drosophila. 38. Genetics 1997 , 147 , 1993–1995. 39. Drosophila transposable elements. Annu. Rev. Charlesworth, B.; Langley, C.H. The population genetics of 1989 , 23 , 251–287. [CrossRef] Genet. Langley, C.H.; Montgomery, E.; Hudson, R.; Kaplan, N.; Charlesworth, B. On the role of unequal exchange in 40. Genet. Res. 1988 , 52 , 223–235. [CrossRef] [PubMed] the containment of transposable element copy number. Schaack, S.; Gilbert, C.; Feschotte, C. Promiscuous DNA: Horizontal transfer of transposable elements and 41. Trends Ecol. Evol. 2010 , 25 , 537–546. [CrossRef] [PubMed] why it matters for eukaryotic evolution. 42. Gilbert, C.; Feschotte, C. Horizontal acquisition of transposable elements and viral sequences: Patterns and consequences. Curr. Opin. Genet. Dev. 2018 , 49 , 15–24. [CrossRef] 43. Hartl, D.L.; Lohe, A.R.; Lozovskaya, E.R. Modern thoughts on an ancyent marinere: Function, evolution, Annu. Rev. Genet. , 31 , 337–358. [CrossRef] regulation. 1997 44. Wallau, G.L.; Capy, P.; Loreto, E.; Le Rouzic, A.; Hua-Van, A. VHICA, a New Method to Discriminate between Vertical and Horizontal Transposon Transfer: Application to the Mariner Family within Drosophila. 2016 , Mol. Biol. Evol. , 1094–1109. [CrossRef] 33 Wallau, G.L.; Capy, P.; Loreto, E.; Hua-Van, A. Genomic landscape and evolutionary dynamics of mariner 45. transposable elements within the Drosophila genus. BMC Genom. 2014 , 15 , 727. [CrossRef] 46. Capy, P.; David, J.R.; Hartl, D.L. Evolution of the transposable element mariner in the Drosophila melanogaster species group. 1992 , 86 , 37–46. [CrossRef] Genetica 47. Maruyama, K.; Hartl, D.L. Evolution of the transposable element mariner in Drosophila species. Genetics 1991 , 128 , 319–329. [PubMed] 48. Lohe, A.R.; Moriyama, E.N.; Lidholm, D.A.; Hartl, D.L. Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements. Mol. Biol. Evol. 1995 , 12 , 62–72. [CrossRef] [PubMed]

13 Genes 2019 10 , 336 13 of 14 , Peccoud, J.; Loiseau, V.; Cordaux, R.; Gilbert, C. Massive horizontal transfer of transposable elements in 49. 2017 , 114 , 4721–4726. [CrossRef] [PubMed] insects. Proc. Natl. Acad. Sci. USA ffi cient to mediate 50. Lampe, D.J.; Churchill, M.E.; Robertson, H.M. A purified mariner transposase is su , 15 transposition in vitro. 1996 EMBO J. , 5470–5479. [CrossRef] 51. Palazzo, A.; Caizzi, R.; Viggiano, L.; Marsano, R.M. Does the Promoter Constitute a Barrier in the Horizontal 2017 , 9 , 1637–1645. Transposon Transfer Process? Insight from Bari Transposons. Genome Biol. Evol. [CrossRef] [PubMed] Sookdeo, A.; Hepp, C.M.; Boissinot, S. Contrasted patterns of evolution of the LINE-1 retrotransposon in 52. 2018 , 9 , 12. [CrossRef] perissodactyls: The history of a LINE-1 extinction. Mob. DNA 53. Cantrell, M.A.; Scott, L.; Brown, C.J.; Martinez, A.R.; Wichman, H.A. Loss of LINE-1 activity in the megabats. Genetics 178 , 393–404. [CrossRef] [PubMed] 2008 , Gene 2012 , 54. Platt, R.N., 2nd; Ray, D.A. A non-LTR retroelement extinction in Spermophilus tridecemlineatus. , 47–53. [CrossRef] 500 Grahn, R.A.; Rinehart, T.A.; Cantrell, M.A.; Wichman, H.A. Extinction of LINE-1 activity coincident with a 55. Cytogenet. Genome Res. 2005 , 110 , 407–415. [CrossRef] major mammalian radiation in rodents. Nuzhdin, S.V. Sure facts, speculations, and open questions about the evolution of transposable element copy 56. 1999 , 107 , 129–137. [CrossRef] [PubMed] number. Genetica Erickson, I.K.; Cantrell, M.A.; Scott, L.; Wichman, H.A. Retrofitting the genome: L1 extinction follows 57. , 12315–12323. [CrossRef] , 85 endogenous retroviral expansion in a group of muroid rodents. 2011 J. Virol. Cantrell, M.A.; Ederer, M.M.; Erickson, I.K.; Swier, V.J.; Baker, R.J.; Wichman, H.A. MysTR: An endogenous 58. retrovirus family in mammals that is undergoing recent amplifications to unprecedented copy numbers. 2005 , 79 J. Virol. , 14698–14707. [CrossRef] 59. Genome Biol. 2011 , 12 , 236. [CrossRef] Deininger, P. Alu elements: Know the SINEs. 60. Yang, L.; Wichman, H.A. Tracing the history of LINE and SINE extinction in sigmodontine rodents. bioRxiv . [CrossRef] 2018 Hartl, D.L.; Lohe, A.R.; Lozovskaya, E.R. Regulation of the transposable element mariner. 61. 1997 , 100 , Genetica 177–184. [CrossRef] 62. Bingham, P.M.; Kidwell, M.G.; Rubin, G.M. The molecular basis of P-M dysgenesis—The role of the P-element, a P-Strain-Specific Transposon Family. Cell , 29 , 995–1004. [CrossRef] 1982 Bucheton, A.; Paro, R.; Sang, H.M.; Pelisson, A.; Finnegan, D.J. The molecular basis of the I-R hybrid 63. : Identification, cloning and properties of the I-Factor. Cell dysgenesis syndrome in Drosophila melanogaster 38 , , 153–163. [CrossRef] 1984 hobo 64. Yannopoulos, G.; Stamatis, N.; Monastirioti, M.; Hatzopoulos, P.; Louis, C. is responsible for the Drosophila melanogater bearing the male recombination factor induction of hybrid dysgenesis by strains of Cell 1987 , 49 , 487–495. [CrossRef] 23.5MRF. 65. Blumenstiel, J.P.; Hartl, D.L. Evidence for maternally transmitted small interfering RNA in the repression of Proc. Natl. Acad. Sci. USA , 102 , 15965–15970. [CrossRef] transposition in Drosophila virilis. 2005 66. Brennecke, J.; Malone, C.D.; Aravin, A.A.; Sachidanandam, R.; Stark, A.; Hannon, G.J. An Epigenetic Role for Science 2008 , Maternally Inherited piRNAs in Transposon Silencing. , 1387–1392. [CrossRef] [PubMed] 322 67. Dorogova, N.V.; Bolobolova, E.U.; Zakharenko, L.P. Cellular aspects of gonadal atrophy in Drosophila P-M hybrid dysgenesis. Dev. Biol. 2017 , 424 , 105–112. [CrossRef] [PubMed] 68. Niki, Y.; Chigusa, S.I. Developmental analysis of gonadal sterility of P-M dysgenesis of Drosophila melanogaster . Jpn. J. Genet. 1986 61 , 147–157. [CrossRef] , Daniels, S.B.; Peterson, K.R.; Strausbaugh, L.D.; Kidwell, M.G.; Chovnick, A. Evidence for horizontal 69. Genetics 1990 , 124 , 339–355. transmission of the P transposable element between Drosophila species. Crow, J.F. The genesis of dysgenesis. Genetics 1988 , 120 , 315–318. 70. 71. Engels, W.R. Invasions of P elements. Genetics 1997 , 145 , 11–15. 72. Crozatier, M.; Vaury, C.; Busseau, I.; Pelisson, A.; Bucheton, A. Structure and genomic organization of I Drosophila melanogaster . Nucleic Acids Res. elements involved in I-R hybrid dysgenesis in , 16 , 9199–9213. 1988 [CrossRef] 73. Bucheton, A.; Simonelig, M.; Vaury, C.; Crozatier, M. Sequences similar to the I transposable element involved Nature in I-R hybrid dysgenesis in occur in other Drosophila species . D. melanogaster 1986 , 322 , 650–652. [CrossRef]

14 Genes 2019 10 , 336 14 of 14 , Bucheton, A.; Vaury, C.; Chaboissier, M.C.; Abad, P.; Pelisson, A.; Simonelig, M. I elements and the Drosophila 74. 1992 , 86 , 175–190. [CrossRef] genome. Genetica Vaury, C.; Abad, P.; Pelisson, A.; Lenoir, A.; Bucheton, A. Molecular characteristics of the heterochromatic 75. I elements from a reactive strain of J. Mol. Evol. 1990 , 31 , 424–431. [CrossRef] . Drosophila melanogaster [PubMed] 76. Grentzinger, T.; Armenise, C.; Brun, C.; Mugat, B.; Serrano, V.; Pelisson, A.; Chambeyron, S. piRNA-mediated 2012 , transgenerational inheritance of an acquired trait. , 1877–1888. [CrossRef] [PubMed] Genome Res. 22 77. Periquet, G.; Hamelin, M.H.; Bigot, Y.; Lepissier, A. Geographical and historical patterns of the distribution Drosophila melanogaster populations. J. Evol. Biol. of hobo elements in , 2 , 223–229. [CrossRef] 1989 78. Drosophila melanogaster and its close relatives. Boussy, I.A.; Itoh, M. Wanderings of hobo: A transposon in 2004 120 , 125–136. [CrossRef] Genetica , 79. Ragagnin, G.T.; Bernardo, L.P.; Loreto, E.L. Unraveling the evolutionary scenario of the hobo element in Drosophila melanogaster and D. simulans in South America using the TPE repeats as markers. populations of Genet. Mol. Biol. 2016 39 , 145–150. [CrossRef] , Boussy, I.A.; Daniels, S.B. hobo transposable elements in 80. and D. simulans . Genet. Res. Drosophila melanogaster 1991 , 58 , 27–34. [CrossRef] 81. Simmons, G.M. Horizontal transfer of hobo transposable elements within the species Drosophila melanogaster Mol. Biol. Evol. , 9 , 1050–1060. [CrossRef] [PubMed] complex: Evidence from DNA sequencing. 1992 82. Lyozin, G.T.; Makarova, K.S.; Velikodvorskaja, V.V.; Zelentsova, H.S.; Khechumian, R.R.; Kidwell, M.G.; Koonin, E.V.; Evgen’ev, M.B. The structure and evolution of Penelope in the virilis species group of Drosophila: J. Mol. Evol. 2001 , 52 An ancient lineage of retroelements. , 445–456. [CrossRef] 83. Evgen’ev, M.B. What happens when Penelope comes?: An unusual retroelement invades a host species ff erent strategies. Mob. Genet. Elem. genome exploring di , 3 , 397–408. [CrossRef] [PubMed] 2013 84. Bartolome, C.; Bello, X.; Maside, X. Widespread evidence for horizontal transfer of transposable elements across Drosophila genomes. Genome Biol. 2009 , 10 , R22. [CrossRef] [PubMed] 85. Hill, T.; Betancourt, A.J. Extensive exchange of transposable elements in the group. Drosophila pseudoobscura 2018 Mob. DNA 9 , 20. [CrossRef] [PubMed] , 86. Gilbert, C.; Schaack, S.; Pace, J.K., 2nd; Brindley, P.J.; Feschotte, C. A role for host-parasite interactions in the horizontal transfer of transposons across phyla. Nature 2010 , 464 , 1347–1350. [CrossRef] © 2019 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http: // creativecommons.org / ). / by / 4.0 / licenses

Related documents