{'fraction_non_alphanumeric': 0.03327071462948875, 'fraction_numerical': 0.028007638446849142, 'mean_word_length': 4.727038321330133, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 5, 'https://': 2, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 5, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Monitoring Insect Transposable Elements in Large Double-Stranded DNA Viruses Reveals Host-to-Virus and Virus-to-Virus Transposition

Abstract The mechanisms by which transposable elements (TEs) can be horizontally transferred between animals are unknown, but viruses are possible candidate vectors. Here, we surveyed the presence of host-derived TEs in viral genomes in 35 deep sequencing data sets produced from 11 host–virus systems, encompassing nine arthropod host species (five lepidopterans, two dipterans, and two crustaceans) and six different double-stranded (ds) DNA viruses (four baculoviruses and two iridoviruses). We found evidence of viral-borne TEs in 14 data sets, with frequencies of viral genomes carrying a TE ranging from 0.01% to 26.33% for baculoviruses and from 0.45% to 7.36% for iridoviruses. The analysis of viral populations separated by a single replication cycle revealed that viral-borne TEs originating from an initial host species can be retrieved after viral replication in another host species, sometimes at higher frequencies. Furthermore, we detected a strong increase in the number of integrations in a viral population for a TE absent from the hosts’ genomes, indicating that this TE has undergone intense transposition within the viral population. Finally, we provide evidence that many TEs found integrated in viral genomes (15/41) have been horizontally transferred in insects. Altogether, our results indicate that multiple large dsDNA viruses have the capacity to shuttle TEs in insects and they underline the potential of viruses to act as vectors of horizontal transfer of TEs. Furthermore, the finding that TEs can transpose between viral genomes of a viral species sets viruses as possible new niches in which TEs can persist and evolve.


Introduction
Like any other genome component, transposable elements (TEs) are vertically transmitted from one generation to the next through reproduction. But TEs can also bypass vertical transmission and cross species boundaries through a process not involving reproduction, called horizontal transfer (HT). The inference of HT of TEs (HTT) derives from the numerous observations that TE sequences from different host organisms show much lower genetic distance than expected given the divergence time of the hosts (Peccoud et al. 2018). Since the seminal report of the P element transfer between Drosophila willistoni and D. melanogaster (Daniels et al. 1990), dozens of studies have characterized HTTs involving many branches of the eukaryote tree Dotto et al. 2018). Large-scale HTT surveys in plants, insects, and vertebrates revealed recurrent transfers, which seeded a large fraction of the TE copies found today in these taxa (Bartolom e et al. 2009;El Baidouri et al. 2014;Ivancevic et al. 2018;Reiss et al. 2019;Zhang et al. 2020). Given the strong impact TEs have on genome structure and dynamics (Cordaux and Batzer 2009;Bourque et al. 2018), HTT is considered as an important process shaping the evolution of eukaryote genomes (Gilbert and Feschotte 2018).
Several important questions about HTT remain, perhaps first and foremost that of the factors facilitating these transfers. Large-scale studies have shown that HTT are more likely to occur between closely related species living in the same biogeographical realm than between more distantly related species living in different realms (Bartolom e et al. 2009;Peccoud et al. 2017). Interestingly, some host taxa such as teleost fish among vertebrates and moths and butterflies among arthropods seem to be more prone to withinphylum HTT than others, a trend that remains unexplained (Reiss et al. 2019;Zhang et al. 2020). Furthermore, many HTT events have been reported that involve parasites and their hosts, suggesting host-parasite relationships may facilitate HTT Kuraku et al. 2012;Walsh et al. 2013;Guo et al. 2014;Suh et al. 2016). However, the molecular processes underlying HTT remain largely unknown.
Several scenarios have been proposed to explain how a TE can escape from a donor organism and enter the germline of a recipient host (Silva et al. 2004;Loreto et al. 2008;Wallau et al. 2012). Experimental observations currently support two possible HTT routes. The first proposes that extracellular vesicles (EVs) could act as vectors of HTT between animals. These 50-500 nm membrane-derived vesicles are secreted by most cell types, may carry proteins, lipids and/or genetic material, and are naturally present in biological fluids (van Niel et al. 2018). EVs were shown to shuttle retrotransposons and to mediate their HT in laboratory conditions between different human cell lines or between cell culture media and mouse cell lines or embryos (Kawamura et al. 2019;Ono et al. 2019). However, the extent to which EVs may also shuttle TEs between species in natural conditions remains to be evaluated.
The second HTT route receiving some experimental support involves viruses . Early studies using low-throughput targeted approaches identified TEs integrated in the genomes of several baculoviruses, which are double-stranded (ds) DNA viruses belonging to the Baculoviridae family, that were passaged in moth cell cultures (Miller and Miller 1982;Fraser et al. 1985Fraser et al. , 1995 or whole larvae (Jehle et al. 1995(Jehle et al. , 1998. These pioneering works demonstrated that during infection, TEs can jump from the host genome to virus genomes, and that baculoviruses can receive and potentially carry a foreign genetic load from their host. More recent works using high-throughput sequencing showed that in addition to viral genomes, multiple host RNA sequences including TEs could be packaged in capsids of RNA viruses (Routh et al. 2012;Ghoshal et al. 2015;Eckwahl et al. 2016;Telesnitsky and Wolin 2016). These results further emphasized the potential role of some viruses as vectors of HTT and suggested that TEs may not have to be integrated into viral genomes to be shuttled by viruses. Using an ultradeep sequencing approach, we revisited early works on baculoviruses and characterized the whole spectrum and frequency of host TEs integrated in genomes of the Autographa californica multiple nucleopolyhedrovirus (AcMNPV) purified from infected moth larvae (Gilbert et al. 2014(Gilbert et al. , 2016. Our results revealed that a large diversity of TEs are able to jump from moth to virus genomes at each infection cycle, with an average of 4.8% of sequenced AcMNPV genomes carrying at least one host TE. Studies of TEs segregating in baculovirus populations raised a number of outstanding questions. First, only a limited number of virus-host systems have been surveyed, such that it is still unclear whether the capacity of viral populations to carry host TEs is widespread among many viruses. Second, evidence exists showing that a virus-borne TE (TCl4.7) acquired by a baculovirus (CpGV) from a given host (the false codling moth Thaumatotibia (¼Cryptophebia) leucotreta) can persist in the virus population during passages of the virus in another host (the codling moth Cydia pomonella) (Jehle et al. 1995). However, we were unable to detect any TE insertion shared between an initial AcMNPV population (called G0) replicated in the cabbage looper (Trichoplusia ni) and populations purified after ten successive infection cycles of the G0 on ten lines of the beet armyworm (Spodoptera exigua) (Gilbert et al. 2016). Thus, it is unclear whether persistence of virus-borne TEs is common among different viruses. In fact, we showed that each individual TE insertion segregates at low frequency in AcMNPV and was always purged out of the viral population within less than ten replication cycles (Gilbert et al. 2016). Third, the genome sequences of T. ni and S. exigua were not available at the time so we could not exclude that some of the AcMNPVborne TEs originated from hosts other than T. ni or S. exigua in which the virus replicated before we conducted our study. Finally, resequencing an AcMNPV population using the longread PacBio technology unveiled many full length TE copies integrated into viral genomes, suggesting that such copies have the capacity to encode the entire machinery necessary to transpose from viral genomes to other DNA molecules (Loiseau et al. 2020). Yet, more direct evidence supporting transposition of virus-borne TEs is still lacking.
Here, we monitored the presence, nature, and frequency of TEs in 35 deep-sequenced viral genomes obtained from 11 virus-host systems involving two iridovirus and four baculovirus species. The finding of moth TEs in non-AcMNPV baculoviruses and in the iridescent virus 6 suggests that the capacity to carry host TEs may be widespread among large dsDNA viruses. Importantly, our results also demonstrate that persistence of virus-borne TEs over more than one replication cycle occurs in different viruses (e.g., AcMNPV and the invertebrate iridescent virus 6 [IIV6] iridovirus) after two successive infection cycles in divergent species (e.g., a fly and a moth). Finally, we show that virus-borne TEs are able to transpose into other viral genomes during the course of an infection cycle.

Overview of Sequenced Host-Virus Systems
A total of 35 Illumina sequencing data sets were analyzed in this study (table 1). One data set was produced by sequencing genomes of the IIV6 purified from D. melanogaster S2 cells ( fig. 1). IIV6 particles purified from S2 cells were then used to infect whole D. melanogaster flies, as well as D. hydei flies and two moth species, the spotted stalk borer (Chilo partellus) and the maize corn borer (Sesamia nonagrioides) (five data sets in total). In addition, six data sets were produced by sequencing genomes of the invertebrate iridescent virus 31 (IIV31) from two crustacean species (the terrestrial isopods Armadillidium vulgare and Porcellio dilatatus) ( fig. 1). We also analyzed a population of the AcMNPV baculovirus initially replicated on T. ni (Chateigner et al. 2015) and here used to Monitoring Insect TEs in dsDNA Viruses . doi:10.1093/molbev/msab198 MBE infect the maize corn borer (one data set) ( fig. 2). Finally, we surveyed 22 baculovirus sequencing data sets produced as part of other studies (Gueli Alletti, Eigenbrod, et al. 2017;Gueli Alletti, Sauer, et al. 2017;Gueli Alletti et al. 2018;, including one data set of the Agrotis segetum granulovirus (AgseGV) purified from the turnip moth (Agrotis segetum), six data sets of the Agrotis segetum nucleopolyhedrovirus B (AgseNPV-B) purified from a cell line of the black cutworm (A. ipsilon), and 15 data sets of the Cydia pomonella granulovirus (CpGV) purified from the codling moth (Cydia pomonella) ( fig. 2). All moth TEs found integrated into viral genomes in this study are listed in supplementary table S1, Supplementary Material online, and their sequence, as well as that of all TEs included in our library, is provided in supplementary file S1, Supplementary Material online.

Moth TEs in Baculoviruses Other Than AcMNPV
We first sought to assess whether TEs can transpose into genomes of baculoviruses other than AcMNPV using bulk deep sequencing, and to characterize the origin of such TEs as well as their frequency in viral populations. In brief, we found a total of six different TEs integrated in 7 out of 22 non-AcMNPV baculovirus data sets, which confirms that baculoviruses other than AcMNPV carry TEs integrated in their genomes (Jehle et al. 1995(Jehle et al. , 1998. More specifically, we began by searching for TEs integrated in genomes of six out of ten successive passages of AgseNPV-B replicated in cells of A. ipsilon. These genomes were sequenced at depths varying from 839Â to 2,930Â. We found TEs integrated in genomes of only one passage (pp7), with a total of 90 TE-virus chimeras covering 13 different positions in the viral genome and corresponding to a frequency of 3.4% of viral genomes carrying  fig. 3; supplementary fig S1., Supplementary Material online). The three TEs were absent from the whole-genome sequence (WGS) of A. ipsilon (NCBI accession number: PNFC00000000.1), the species from which the cell line used to passage AgseNPV-B was established. The WGS of the A. ipsilon cell line is not available. Thus, the origin of the TEs found in pp7 cannot be easily explained. In any case, the piggybac and Sola TEs were identified in our assembly of nonviral reads obtained from the AgseGV/A. segetum larvae system and the nonautonomous element was previously identified in the plutellid moth Plutella xylostella (Han et al. 2016). The fact that we only retrieved these TEs in passage pp7 of the virus and not in earlier (pp0, pp1, pp3, pp5) or later (pp10) passages may be due to the stochasticity inherent to the detection of low frequency virus-borne TEs, which is likely pronounced given the relatively low depths at which these samples were sequenced. We then searched for TEs integrated in genomes of AgseGV replicated in A. segetum larvae and sequenced at 7,997Â. A total of 33 TE-virus chimeras were detected, covering 28 different positions and corresponding to a frequency of 0.3% of viral genomes carrying an insertion (table 1). All chimeras mapped to the aforementioned Sola TE (supplementary table S1, Supplementary Material online).
We also detected TE-virus chimeras in 5 of the 15 CpGV isolates replicated in C. pomonella larvae and sequenced at depths varying from 523Â to 3,809Â. and table S1, Supplementary Material online). All three TEs are present in the C. pomonella WGS, suggesting they transposed from the C. pomonella genome to that of the virus.

Retrotransposons from D. melanogaster S2 Cells Transpose into the IIV6 Iridovirus
To assess the potential of large dsDNA viruses other than baculoviruses to shuttle TEs between arthropods, we searched for TEs integrated in genomes of the IIV31 iridovirus purified from three P. dilatatus and three A. vulgare individuals (including one previously analyzed by Loiseau et al. [2020]) and in genomes of the IIV6 iridovirus after a passage on D. melanogaster S2 cells. Sequencing depths varied from 112,750Â to 220,537Â (table 1). Although all IIV31 data sets were found to be entirely devoid of TE-virus chimeras, we found thousands of low-frequency TE copies integrated in IIV6 genomes (table 1). This shows that the presence of unfixed TEs in viral populations is not a feature restricted to baculoviruses, but that it is also not systematic. FIG. 2. Overview of baculovirus infections and number and frequency of TEs found in viral genomes. (A) An AcMNPV isolate called "G0" purified from Trichoplusia ni larvae (Gilbert et al. 2016) was used to infect Sesamia nonagrioides larvae. The AcMNPV population purified from these larvae was called "G1." (B) Illumina reads produced as part of other studies for three other baculoviruses were surveyed for the presence of TEs integrated into viral genomes (Gueli Alletti, Eigenbrod, et al. 2017;Gueli Alletti, Sauer, et al. 2017;Gueli Alletti et al. 2018;. Monitoring Insect TEs in dsDNA Viruses . doi:10.1093/molbev/msab198 MBE More specifically, we found 6,044 TE-virus chimeras in IIV6 genomes purified from D. melanogaster S2 cells, corresponding to 2,831 different positions along the viral genome ( fig. 4) and 4.33% of sequenced IIV6 genomes carrying at least one TE (table 1). TE-virus junctions involved one D. melanogaster non-LTR retrotransposon and 19 LTR retrotransposons including five (MDG1_LTR, DM297_LTR, DM176, IDEFIX_LTR, and Drosophila_melanogaster: rnd-3_family-276#LTR/Gypsy) that accounted for 500 or more (up to 3,192) chimeras (supplementary table S1, Supplementary Material online).
As this is the first report of TEs segregating in an iridovirus population, we set out to independently validate the biological nature of TE-virus chimeras by comparing these chimeras with those involving human DNA added to the IIV6 DNA sample prior to constructing the sequencing library. We investigated the structure of chimeras occurring between human TEs and the IIV6 genome. These can only be technical as IIV6 has not been in contact with human cells. Given such chimeras are not generated by transposition, they are not expected to preferentially involve the extremities of TEs.

LTR Retrotransposons from D. melanogaster Whole Flies Transpose into the IIV6 Iridovirus
The results presented above demonstrate that D. melanogaster TEs can transpose from a cell line (S2 cells) to IIV6 genomes. Here, we now show that TEs can also transpose from whole insects into IIV6 genomes. We infected two batches of 80 D. melanogaster adult flies with the IIV6 population purified from S2 cells and sequenced IIV6 viral genomes purified from these infections at 170,859Â and 211,039Â (fig. 1). We found more TE-virus chimeras in the two IIV6/adult flies' batches (12,193 and 12,005 in batches 1 and 2, respectively) than in the parental IIV6/S2 cells data set (6,044) and a higher frequency of viral genomes carrying at least one TE (7.36% and 5.86% in batches 1 and 2, respectively) than in the parental IIV6/S2 cells population (4.33%) (table 1; supplementary table S1, Supplementary Material online). A total of 28 and 29 different TEs were identified in batches 1 and 2, respectively (supplementary tables S1 and S4, Supplementary Material online). All 20 TEs identified in the parental IIV6 purified from S2 cells were found to be integrated at 2,263 and 2,345 different positions in IIV6 from whole flies (batches 1 and 2, respectively) (black and green bars in D. melanogaster 1 and 2 graphs in fig. 4). In principle, these TE insertions could have two sources. First, they could have been present in the parental IIV6 purified from S2 cells and have persisted in the viral population during replication in whole flies, implying that IIV6 genomes bearing TE insertions were encapsidated again in whole flies. Second, the TE insertions could result from transposition of TE copies located in the genomes of whole flies. Given that the TE content of D. melanogaster S2 cells and whole flies genomes is likely very similar, it is difficult to assess whether and what portion of insertions of the shared TEs come from the parental IIV6 population or from whole flies.
Interestingly, the majority of TE-virus chimeras found in IIV6 purified from whole flies (6,154 and 5,706 in batches 1 and 2, respectively) map to GYPSY6_LTR and GYPS-7_Dsim, two LTR retrotransposons not detected in the parental IIV6 purified from S2 cells (orange bars in D. melanogaster 1 and 2 graphs in fig. 4; fig. 3; supplementary fig. S3 and tables S1 and S4, Supplementary Material online). Another six and seven TEs absent from IIV6 purified from S2 cells were found in batches 1 and 2, respectively, including one (DMLTR5) involved in more than 25 TE-virus chimeras in the two whole fly batches (supplementary fig. S3 and tables S1 and S4, Supplementary Material online). Given their absence in the parental IIV6/S2 cell data set, de novo integration of these TEs in IIV6 genomes is most likely due to transposition of TE copies located in the genomes of whole flies.

TEs from Other Hosts in IIV6
Here, we show that a TE from a moth host can transpose into IIV6 genomes, though at low frequency. We again used IIV6 particles purified from S2 cells to infect whole flies from another Drosophila species (D. hydei) and larvae from the C. partellus moth ( fig. 1). We also reanalyzed another IIV6 sequencing data set from Loiseau et al. (2020) generated by infecting larvae from the S. nonagrioides moth with the IIV6 particles purified from S2 cells. In Loiseau et al. (2020), we only searched for S. nonagrioides TEs integrated into IIV6 genomes and found none. Here we extended our search to a much larger diversity of TEs, including D. melanogaster TEs. Sequencing depth varied from 52,931Â to 325,206Â (table 1). Across the three species, we found a single non-D. melanogaster TE integrated in IIV6, that is, a piggybac transposon involved in seven TE-virus chimeras in the C. partellus data set (supplementary table S1, Supplementary Material online). There is one copy identical to this TE over its full length in the C. suppressalis WGS and we independently validated its presence in the C. partellus genome by polymerase chain reaction (PCR) and Sanger sequencing (piggyBac_Chilo_partellus in supplementary file S1, Supplementary Material online). Thus, this TE likely transposed from C. partellus to IIV6 during replication of the virus in the moth larvae.

Persistence of Fly TEs during IIV6 Replication
Here, we demonstrate that many D. melanogaster TEs present in the IIV6 isolate extracted from D. melanogaster S2 cells can be retrieved in this virus after replication in another host. Contrasting with the absence or paucity of D. hydei, S. nonagrioides and C. partellus TEs integrated into IIV6 genomes, Loiseau et al. . doi:10.1093/molbev/msab198 MBE we found 1,805, 1,962 and 1,151 TE-virus chimeras involving D. melanogaster TEs in IIV6 genomes purified from the three respective hosts. All D. melanogaster TEs recovered in these IIV6 populations were present in the parental IIV6 purified from S2 cells and the eight most frequent TEs present in the IIV6/S2 cells data set (>0.1%) are also the most frequent in the daughter IIV6 populations (supplementary fig. S3 and table S4, Supplementary Material online). We verified that none of these D. melanogaster TEs is present in the genomes of D. hydei, S. nonagrioides, or C. partellus using BlastN similarity searches. Thus, we conclude that D. melanogaster TEs integrated into IIV6 genomes purified from these three hosts come from the parental IIV6 population purified from S2 cells. The various chimeric reads involving D. melanogaster TEs found in IIV6 purified from D. hydei, S. nonagrioides and C. partellus correspond to 450, 58 and 22 different positions along the IIV6 genome, respectively (table 1, fig. 4). Consistent with persistence of D. melanogaster TEs in IIV6 purified from the three other hosts, 149, 31, and 12 of these positions were shared with the parental IIV6 population purified from S2 cells, respectively. This implies that many TEbearing IIV6 genomes persisted through one replication cycle and were again encapsidated in three different hosts. By contrast, 232, 27, and 10 insertions found in IIV6 purified from D. hydei, S. nonagrioides, and C. partellus were not found in the parental IIV6 purified from S2 cells. A majority (94% in D. hydei) of these insertions were supported by less than ten chimeric reads. Thus, insertions not shared with the parental IIV6 population purified from S2 cells may in fact have been present at very low frequency in this population but not sequenced. Interestingly, it is also possible that some of them were generated through virus-to-virus transposition  1). Black bars correspond to TE insertions shared by the parental IIV6 population and at least one daughter IIV6 population. These shared insertions correspond to the same TE at the exact same position in the parental IIV6 population purified from D. melanogaster S2 cells and other IIV6 populations purified from other hosts. Green bars correspond to TE insertions not shared with the parental IIV6 population. TEs involved in these unshared insertions were present in the parental IIV6 population but not at the position they were found in the daughter populations. These unshared insertions may be due to de novo virus-to-virus transposition during infection of Drosophila whole flies or moth larvae. Alternatively, they may correspond to TE insertions that were present in the parental population but were not detected by sequencing. Orange bars correspond to TEs that are absent from the parental IIV6 population. These are likely to result from de novo transposition from the host genome. Numbers of insertions are given in supplementary

Persistence of TEs Integrated in AcMNPV Genomes
To assess whether TEs integrated into AcMNPV genomes can be recovered after one replication cycle, as observed for IIV6, we infected S. nonagrioides larvae with a viral population (termed "G0") purified from T. ni larvae (Chateigner et al. 2015) and known to contain thousands of TE copies belonging to at least 29 TE families (Gilbert et al. 2016). Our search for TEs integrated into the AcMNPV population purified from S. nonagrioides (hereafter termed "G1") unveiled a total of five TEs, two of which were also present in the G0 population: Tni_contig_27 (piggybac) and Tni_contig_21 (mariner) (supplementary table S1, Supplementary Material online). A BlastN similarity search of these two TEs in the S. nonagrioides genome did not reveal any significant hit and a PCR screening using three primer sets designed to amplify three regions of these elements supported their absence from S. nonagrioides. Thus, we conclude that Tni_contig_27 and Tni_contig_21 present in AcMNPV genomes purified from S. nonagrioides were most likely carried over from the G0 population. These results confirm that some TEs integrated into AcMNPV genomes can persist over at least one replication cycle, as observed for IIV6. Interestingly, aside from the shared presence of Tni_contig_27 and Tni_contig_21, the TE landscape of the G0 and G1 populations markedly differed, as none of the 27 other TEs present in the G0 population was recovered after replication on S. nonagrioides ( fig. 5). This is in stark contrast with IIV6, for which TE contents and frequencies were relatively similar between populations purified from D. melanogaster S2 cells and D. hydei.

AcMNPV-to-AcMNPV Transposition in Moths
Here, we compare the number of TE insertions in two AcMNPV populations separated by one replication cycle and found strong evidence supporting virus-to-virus transposition. A striking feature of the G1 population purified from S. nonagrioides is that 26.33% of AcMNPV genomes carried at least one TE insertion. This high frequency was mainly explained by Tni_contig_27, as the frequency of AcMNPV genomes bearing other TEs is only 0.04%. Compared with 0.38% in the G0 population purified from T. ni larvae (Gilbert et al. 2016), the frequency of genomes carrying a Tni_contig_27 insertion underwent a 69-fold increase during replication in S. nonagrioides. Similarly, the number of different Tni_contig_27 insertions was much higher in the G1 (n ¼ 417) than in the G0 (n ¼ 17) AcMNPV population ( fig.  6). Though this difference may at first sight be explained by virus-to-virus transposition, two other hypotheses have to be considered. First, the G0 population was sequenced at a higher depth (124,221Â) than the G1 population (82,103Â), invalidating the hypothesis according to which many Tni_contig_27 insertions would have been missed in the G0 AcMNPV population. Second, fluctuations in insertion frequencies driven by genetic drift could explain the difference between the G0 and G1 populations. Under genetic drift, the frequency of an insertion has a probability of 0.5 to have increased during replication in S. nonagrioides larvae.
In such a scenario, all 417 insertions found in the G1 population would have been present but missed because of their frequency being too low. Given that 12 out of the 17 G0 insertions are also found in the G1 population, this scenario would imply that a total of 417 þ 17 -12 ¼ 422 Tni_contig_27 insertions were present in the G0 population.
Assuming independence between insertions, the probability that as many as 417 insertions out of 422 increased in frequency due to genetic drift alone is almost null (P ¼ 2.2 Â 10 À16 , exact binomial test). Though inconsistent with drift alone, such a global increase could in principle be observed if all or most Tni_contig_27 insertions were beneficial to the virus. Such a scenario, whereby several hundreds of Tni_contig_27 insertions would be adaptive irrespective of their location along the AcMNPV genome, is highly unlikely. In fact, current evidence strongly suggests that the vast majority of TE insertions are deleterious to viruses as they very rarely reach high frequencies during viral replication in natural hosts (Gilbert et al. 2016;Gilbert and Cordaux 2017). This is best illustrated by the rare occurrence of TEs in consensus viral genomes (Sun et al. 2015;Gilbert and Cordaux 2017;Fil ee 2018). Thus, the nearly 25-fold increase in the number of Tni_contig_27 insertions observed in the G1 (from 17 to 417) is unlikely to be explained by fluctuation in insertion frequency. Instead, we contend that during replication in S. nonagrioides larvae, the Tni_contig_27 TE most likely transposed between AcMNPV genomes, and possibly within AcMNPV genomes as well.
High Frequency of a Single Tni_Contig_27 Insertion in the AcMNPV GTA Gene In this section, we show that our ability to detect virus-tovirus transposition of the Tni_contig_27 TE is likely due to the fact that a single insertion of this TE underwent a sharp increase in frequency during replication of AcMNPV in S. nonagrioides. We noticed that although virus-to-virus transposition likely explains the increase in number of Tni_contig_27 insertions in the AcMNPV G1 population, it does not account for the overall increase in AcMNPV genomes bearing a Tni_contig_27 insertion (from 0.38% in the G0 to 26.29% in the G1 population). In fact, the vast majority of TE-virus chimeras involving Tni_contig_27 in the G1 population (1,069 out of 2,245 without PCR duplicates, or 35,189 out of 37,897 including PCR duplicates) map to a single insertion site, located at position 30,918 of the WP10 AcMNPV genome, in the AcMNPV global transactivator (Ac-GTA) gene ( fig. 6). This insertion, together with 11 other Tni_contig_27 insertions (green stars in fig. 6), is present in the parental G0 population but at a much lower frequency. The frequency increase of this insertion during replication in S. nonagrioides larvae may be due to repeated integration of Tni_contig_27 at the exact same site and/or to preferential replication of viral genomes bearing this insertion. Our earlier analysis of ten AcMNPV populations obtained after infection of T. ni larvae with the G0 population did not reveal any frequency increase of this insertion (Gilbert et al. 2016). Repeated integration at position 30918 may have Loiseau et al. . doi:10.1093/molbev/msab198 MBE occurred but it is unlikely to solely explain the sharp increase in frequency of the Tni_contig_27 insertion during replication in S. nonagrioides larvae. Instead, the frequency of this insertion may have mainly increased through replication of the viral genomes bearing it. It is noteworthy that Ac-GTA is involved in transcription regulation, DNA recombination and repair, chromatin unwinding, and other functions (Rohrmann 2019). Deletion of this gene from the Bombyx mori nucleopolyhedrovirus (BmNPV) only led to mild effects on the infection outcome (i.e., delayed killing time), suggesting that it may not be essential for viral replication in some contexts (Katsuma et al. 2008). We reasoned that the increase in frequency of the Tni_contig_27 insertion in this gene may be due to low functional constraints inducing relaxation of purifying selection. If true, other mutations inactivating this gene may also be expected to have increased in frequency during replication in S. nonagrioides. However, we found no evidence of other high-frequency TE insertions or high-frequency indels in Ac-GTA ( fig. 6; supplementary fig. S4, Supplementary Material online). Thus, the strong increase in frequency of the Tni_contig_27 insertion in Ac-GTA is unlikely to be due to the absence of functional constraints acting on this gene in our experiment. The possibility that this insertion increased in frequency because of a positive effect on viral replication cannot be excluded at this stage. Another possibility worth testing in future experiments is that genomes bearing this insertion may have a higher propensity to be encapsidated compared with other genomes.

Origin of TEs Found in AcMNPV
Though our earlier study identified Tni_contig_27 and Tni_contig_21 into AcMNPV genomes purified from T. ni larvae, we were unable to trace the origin of these TEs because no WGS was available for T. ni at the time (Gilbert et al. 2016). A BlastN similarity search for these TEs against the two now available T. ni genomes (Chen et al. 2019;Talsania et al. 2019) revealed no hit. This result suggested the presence of Tni_contig_27 and Tni_contig_21 in the G0 AcMNPV population purified from T. ni larvae did not result from transposition of TEs located in the T. ni genome. Instead, the two TEs must have originated from an earlier replication cycle of AcMNPV on another host. In addition to Tni_contig_27 and Tni_contig_21, we found three S. nonagrioides DNA TEs (1 piggybac and 2 sola) involved in 5, 6, and 5 TE-AcMNPV chimeras, respectively (supplementary table S1, Supplementary Material online). Thus, TEs became integrated into AcMNPV genomes not only through virus-to-virus transposition but also through transposition of copies located in the S. nonagrioides genome. In turn, this result shows that TEs from yet another host can transpose into AcMNPV genomes  Tni_Contig_11  Tni_Contig_12  Tni_Contig_13  Tni_Contig_14  Tni_Contig_15  Tni_Contig_16  Tni_Contig_20  Tni_Contig_21  Tni_Contig_22  Tni_Contig_23  Tni_Contig_24  Tni_Contig_25  Tni_Contig_27  Tni_Contig_6  Tni_Contig_8  Tni_Contig_9 FIG. 5. Frequency of specific TE insertions in an AcMNPV population purified from Trichoplusia ni (G0) plotted against the frequency of the same TE insertions in an AcMNPV population purified from Sesamia nonagrioides (G1). The AcMNPV baculovirus was first replicated and purified from T. ni larvae (Chateigner et al. 2015;Gilbert et al. 2016). This isolate was called generation 0 (G0). The G0 was here subsequently replicated on and purified from S. nonagrioides larvae. Frequency is expressed in logarithm. The zero on the Y and X axis are for illustrative purposes as zeros cannot be shown at the log scale. Only two TEs, Tni_Contig_27 (turquoise cross symbol) and Tni_Contig_21 (gray triangle symbol), have shared insertions (i.e., located at the same position in the viral genome) between both AcMNPV populations. Only TEs inserted at least at five different positions along the AcMNPV genome are shown.
Monitoring Insect TEs in dsDNA Viruses . doi:10.1093/molbev/msab198 MBE during the course of an infection, though at lower rates than previously reported in T. ni or S. exigua (Gilbert et al. 2016).

HT of TEs Found Integrated into Viral Genomes
A previous survey showed that 8 of the 29 D. melanogaster TEs we found integrated into IIV6 genomes have undergone HT between D. melanogaster and its sister species D. simulans, which diverged less than 5 Ma (supplementary table S4 . 7). When BLAST searches were launched using these seven TEs as queries against all insect WGS available as of March 2020, we retrieved hits longer than 200 bp over at least 80% of the TE length and with >80% nucleotide identity. The synonymous distances (dS) calculated for each of these hits all fall below the 0.5% quantile of the distribution of orthologous gene dS (supplementary fig. S5, Supplementary Material online), strongly suggesting that these TEs were horizontally transferred between species, rather than vertically inherited from a common ancestor. Interestingly, all insect species in which we found these TEs are lepidopterans, perhaps

Discussion
The diversity and frequency of TEs integrated in baculovirus genomes previously assessed using ultra-deep sequencing identified 29 and 40 different TEs integrated in AcMNPV replicated in T. ni and S. exigua, respectively, with frequencies ranging from 1.1% to 14.3% of viral genomes carrying at least one host-borne TE (Gilbert et al. 2016;Loiseau et al. 2020).
Here, we show that TEs from another host (S. nonagrioides) can transpose into AcMNPV, though at lower rates than in T. ni and S. exigua. We further report measurable frequencies of viral genomes carrying TEs in another NPV (AgseNPV) replicated in A. ipsilon cells as well as in two granuloviruses (AgseGV and CpGV) replicated in A. segetum and C. pomonella larvae, respectively. We also show that large numbers of D. melanogaster TEs, as well as a few C. partellus TEs, can integrate in genomes of the IIV6 iridovirus during replication in D. melanogaster cells or in whole flies and moths. Interestingly, earlier works (Bartolom e et al. 2009) and our own inferences of HTT based on comparisons of synonymous distances between TEs and host genes show that many TEs we found integrated in viral genomes have been transmitted through HT between insect species. Altogether, these results extend our knowledge of host-virus systems for which TEs  MBE have been demonstrated to jump into viral genomes. They are in line with previous studies that support a role for large dsDNA viruses as possible vectors of HTT between eukaryotes (Miller and Miller 1982;Carstens 1987;Fraser et al. 1995;Jehle et al. 1995Jehle et al. , 1998Sun et al. 2015;Gilbert and Cordaux 2017).

Highly Contrasted TE Contents in Populations of Different Viruses
Although this study demonstrates the capacity of some viruses to receive and carry TEs, it also reveals important variation in occurrence and frequency of TE integration into viral genomes across host-virus systems. No TE integration was detected in IIV31 iridoviruses purified from two isopod species (A. vulgare and P. dilatatus), in agreement with an earlier report (Loiseau et al. 2020). Although we retrieved many D. melanogaster TEs and a few from C. partellus in IIV6, no TE from D. hydei or S. nonagrioides transposed into genomes of this virus. We found moth TEs integrated in AgseNPVs and CpGVs purified from only 1 out of 6 and 4 out of 15 viral isolates, respectively. Thus, although baculoviruses and iridoviruses have the capacity to encapsidate and shuttle TEbearing viral genomes, transposition into viral genomes does not systematically occur during infection. It may not even occur at all in some host-virus systems. We cannot exclude that some TE insertions may have been missed, especially in some relatively low sequencing-depth data sets (e.g., lower than 1,000Â for six baculoviruses), or because our TE library may not include all host TEs due to unavailable WGS (e.g., C. partellus). Yet, we did not find any TE insertion in several ultra-deeply sequenced viral genomes that had been purified from hosts with available WGS (i.e., A. vulgare, S. nonagrioides, D. hydei; table 1). The absence of TEs in these data sets strongly suggests biological underpinnings.

Factors Possibly Influencing TE Content in Viral Populations
Multiple host and virus factors are likely involved in shaping the numbers of host TEs found integrated in viral genomes. On the host side, the level of TE activity may be involved. For example, the absence of A. vulgare, S. nonagrioides and D.
hydei TEs in IIV31, AcMNPV and IIV6 genomes, respectively, could be due to the lack of currently active TEs in the genome of these three species. However, the TE landscape of A. vulgare and D. hydei (and TE expression patterns for A. vulgare) does not support this hypothesis, as a relatively large fraction of TE copies is nearly identical to their cognate family consensus sequence, which is highly suggestive of current activity (Becking et al. 2020;supplementary fig. S6, Supplementary Material online). In addition, the finding of three S. nonagrioides TEs integrated into AcMNPV genomes provides direct evidence that some TEs are indeed active in this moth. On the virus side, we are not aware of any major difference between baculovirus and iridovirus replication that could explain varying propensities of these viruses to receive and carry hosts TEs. In particular, replication involves transportation of viral DNA into the nucleus for both virus types. For baculoviruses, genome replication and assembly of nucleocapsids fully take place in the virogenic stroma within the nucleus (Rohrmann 2019). For iridoviruses, the first step of genome replication also occurs in the host nuclei, in which the viral polymerase synthesizes unit to twice unit sizes of the viral genome (Williams et al. 2005). The proximity of host and viral genomes during replication in host nuclei may facilitate hostto-virus transposition for the two virus types. However, among baculoviruses, the median frequency of viral genomes carrying a TE is significantly lower for GVs than for NPVs (Wilcoxon unpaired test, W ¼ 0; P-value ¼ 3.4 Â 10 À4 ). Complementation between genomes may be more frequent in multicapsid NPVs, such as AcMNPV and AgseNPV, which occlusion-derived virions (ODV) contain many nucleocapsids (and thus genomes), than in GVs which encapsidate only one genome per ODV and OB (Rohrmann 2019). Thus, defective TE-bearing NPV genomes may be more likely to be replicated than defective TE-bearing GVs, which may explain the trend we observe. In fact, many NPVs are known to be transmitted as multicapsid virions, whereas other large dsDNA viruses, including iridoviruses, are transmitted as monocapsid virions. This may in part explain why no A. vulgare and P. dilatatus TE and only few C. partellus TEs were found in iridovirus genomes. If true, NPVs would stand out as potentially being more efficient HTT vectors than other viruses. In turn, given that NPVs mainly infect lepidopterans (Goulson 2003;Th ez e et al. 2018), their high propensity to shuttle TEs may explain why lepidopterans are seemingly more prone to HTT than other insects, as proposed by Reiss et al. (2019). However, the trend underlined above indicating a higher frequency of virusborne TEs in NPVs than in GVs should be interpreted cautiously as it relies on a limited number of viruses. A more thorough assessment of the higher tolerance of multicapsid viruses to TEs should include single-capsid NPVs (such as the Helicoverpa zea nucleopolyhedrovirus for example). And in fact, the frequency of viral genomes carrying a TE is as high as or even higher in IIV6 purified from Drosophila S2 cells or whole flies (from 3.6% to 6%) than in NPVs, which contradicts the view that the mono-nucleocapsid nature of iridoviruses would dampen their ability to carry TEs. In this context, it is noteworthy that IIV6, which was originally sampled from the moth C. fumiferana (Fukaya and Nasu 1966), is not known to naturally infect flies. Thus, this artificial host-virus system involving a passage in cell culture and infection of whole flies by injection may have created favorable conditions to transposition and persistence of fly TEs in IIV6. In any case, our results call for future studies specifically dedicated to deciphering the relative contribution of host and virus factors involved in shaping transposition of host TEs into viral genomes, as well as TE persistence in viral populations over multiple replication rounds.

Persistence of Viral-Borne TEs over Two Successive Replication Cycles in Different Hosts
Importantly, this study provides direct evidence that TEs from a given insect species integrated into genomes of a viral population can be recovered after subsequent purification of the same virus from another species. Our previous study (Gilbert et al. 2016) showed that the number and frequency of TEs could be similar in a parental (G0) AcMNPV population Loiseau et al. . doi:10.1093/molbev/msab198 MBE replicated in one host (T. ni) and its daughter populations (G10) separated from G0 by ten successive infection cycles on another host (S. exigua). However, no TE was shared between the two populations (Gilbert et al. 2016). This result revealed a continuous dynamics of gain (via transposition) and loss (via drift and/or purifying selection) of TEs in AcMNPV populations and showed that the persistence of a given TE community in AcMNPV populations was less than ten infection cycles. Here, we found D. melanogaster TEs integrated in IIV6 genomes purified from two moths (C. partellus and S. nonagrioides) and another fly (D. hydei). We also recovered TEs present in the AcMNPV G0 data set in AcMNPV genomes purified from S. nonagrioides larvae. This demonstrates that TE-bearing viral genomes can be encapsidated, as shown earlier for the TCp3.2 and TCl4.7 TE in CpGV mutants MCp4 and MCp5, respectively (Jehle et al. 1995(Jehle et al. , 1998, and that persistence of virus-borne TEs over at least few infection cycles may be common among large dsDNA viruses. Interestingly, we observe large differences in TE diversity and number of viral genomes bearing TE insertions after one replication cycle of the same starting viral population in different hosts. For example, although the frequency of TE-bearing genomes in the IIV6 population purified from D. hydei (3.37%) is close to that of the parental IIV6 population purified from S2 cells (4.33%), it dropped down to 1.19% and 0.45% in IIV6 purified from C. partellus and S. nonagrioides, respectively. Furthermore, although many TEs originating from the genome of S2 cells were recovered in IIV6 purified from D. hydei, C. partellus and S. nonagrioides, only 2 out of 29 TEs initially found in AcMNPV purified from T. ni persisted during replication of this virus in S. nonagrioides. Such differences may be due to host-virus interactions other than hostto-virus transposition. Interstingly, the two transposoncarrying CpGV mutants identified in Jehle et al. (1995Jehle et al. ( , 1998 showed similar virulence parameters, in terms of median lethal concentration and survival time, and virus progeny production as the wild-type (wt) CpGV, but they were efficiently out-competed by the wt-CpGV, when host larvae or cultured host cells were co-infected with one of the mutants and wt-CpGV (Arends et al. 2005). It would be interesting to further evaluate the role played by the host in shaping TE content and frequency in viral population in more controlled viral replication assays.

Virus-to-Virus Transposition
Interestingly, this study also showed that most of the TE insertions present in AcMNPV genomes purified from S. nonagrioides were generated by virus-to-virus transposition. This remarkable observation shows that virus-borne TEs can transpose in new genomes (here a virus), even if the recipient species is distantly related to the host from which the TEs originated (S. nonagrioides and T. ni diverged 60 Ma [Toussaint et al. 2012]). It confirms and extends earlier findings based on long-read sequencing showing that TEs can integrate as full-length copies into AcMNPV genomes and are thus likely able to transpose from the virus genome to another DNA molecule (Loiseau et al. 2020). Our results further indicate that only some of the TEs found in a given viral population may undergo measurable virus-to-virus transposition during subsequent infection cycles. Indeed, strong evidence of virus-to-virus transposition in AcMNPV purified from S. nonagrioides was found for only one TE (Tni_contig_27) out of the 29 present in the AcMNPV G0 population purified from T. ni. It is noteworthy that this TE is the second most frequent autonomous TE in the AcMNPV G0 population, being present in 0.38% of AcMNPV genomes (Gilbert et al. 2016). This suggests that the frequency of virusto-virus transposition events involving a TE may depend on the TE frequency in the infecting viral population.
Tni_contig_27 is the only TE present in the G0 population for which a single insertion underwent a sharp increase in frequency during replication in S. nonagrioides larvae. We speculate that the increase of the Tni_contig_27 insertion located in the Ac-GTA gene may not result from transposition. Instead, the frequency of this insertion may have increased through viral replication that may be driven by positive selection and/or by a higher propensity of viral genomes bearing this insertion to be encapsidated. In turn, the large number of Tni_contig_27 copies present in the viral population could have allowed this TE to reach a high enough rate of virus-to-virus transposition to be detectable by our approach. If so, many other TEs may have transposed between viral genomes at similar rates per TE copy, but these events would have been too rare to be detected as all these TEs were hundreds to thousands of times less frequent than Tni_contig_27 in the viral population. In any case, many TEs, host and virus factors are likely to shape the likelihood for a TE to undergo virus-to-virus transposition. Altogether, our results suggest that viruses may not only serve as launching platforms for TEs to colonize naive cellular genomes (a hypothesis that still needs to be formally tested), but they may also be viewed as additional niches in which TEs could persist through time and evolve under different constraints than in cellular hosts. This is reminiscent of the molecular symbiosis hypothesis proposed by Fil ee (2018) to characterize the intricate interactions occurring between various types of mobile genetic elements and giant viruses infecting single-cell eukaryotes.

Infection of D. melanogaster S2 Cells with IIV6
The IIV6 viral strain used to infect D. melanogaster S2 cells is the one originally described in Fukaya and Nasu (1966). S2 cells were collected from two confluent T75 flasks and transferred into 50-ml tubes ( fig. 1). The cells were counted and 2 Â 10 8 cells were transferred to a fresh 50-ml tube. The cells were pelleted at 500 Â g for 10 min at room temperature. Then, 5 ml of an IIV6 suspension was added to the cell pellet at a multiplicity of infection of 0.01. Cells were kept under the hood for 1 h with gentle inversions every 10 min so that the cells stayed in suspension. After 1 h, 1.25 ml (i.e., 0.5 Â 10 8 cells) was added into each of the four T75 flasks containing 13.5 ml of complete medium. The cells were incubated at 25 C for 5 days. Cells were detached by pipetting up and down with a 10-ml pipet and transferred into a 50-ml tube. The cells Monitoring Insect TEs in dsDNA Viruses . doi:10.1093/molbev/msab198 MBE were broken down by three cycles of freezing at À20 C and thawing at 37 C. Cell debris were pelleted by centrifugation at 5,000 Â g for 30 min at 4 C. Supernatant was poured into an ULTRACLEAR tube (Beckman #361706) and underlaid with 1.5 ml of 30% (wt/wt) sucrose in 50 mM HEPES, pH 7, 0.1% BSA. Viral particles were pelleted at 35,000 Â g for 90 min at 4 C. The resulting virus pellet, showing the characteristic blue opalescence of insect iridoviruses, was resuspended in 300-400 ml 10mM Tris-HCl buffer (pH 7.2), transferred to an Eppendorf tube and debris were pelleted by centrifugation for 5 min at 3,800 Â g at 4 C. The supernatant was aliquoted, the viral titer determined (4.42 Â 10 10 PFUs/ml) and this stock suspension was stored at À80 C. Viral DNA was then extracted from an aliquot of this suspension using the QIAamp DNA Mini Kit (Qiagen). The 80 individuals were first grinded with a plastic pestle in Tris-HCl buffer (10 mM; pH 7.2) and two 5-min centrifugation steps at 500 Â g were performed to eliminate most of host cells and tissues. Then, an ultracentrifugation step on sucrose cushion was performed at 35,000 Â g for 90 min at 4 C as described above. The pellet was resuspended in 100 ml of Tris-HCl buffer and viral DNA was extracted using the QIAamp DNA Mini Kit (Qiagen).

Infections of C. partellus with IIV6
Ten fourth instar C. partellus larvae were infected with the IIV6 stock suspension purified from Drosophila S2 cells (see above and fig. 1). Larvae were pricked using a thin needle soaked in the viral suspension. Fourteen days later, the larvae presented a purple iridescence and they finally died about 4 weeks after infection. Upon host death, viral particles were filtered through cheesecloth and two centrifugation steps were performed to eliminate most of host cells and tissues. Virus purification and DNA extraction were performed as for IIV6 in Drosophila flies (see above).

Infections of A. vulgare and P. dilatatus with IIV31
The IIV31 virus used to infect A. vulgare and P. dilatatus isopods was obtained through grinding a piece of cuticle from a naturally infected A. vulgare individual collected on the campus of the University of California Riverside (same viral isolate as in Loiseau et al. [2020]). Three P. dilatatus and two A. vulgare individuals were pricked with a thin needle soaked in the virus suspension ( fig. 1). Four weeks postinfection, bluish dead isopods were individually placed and crushed in a 1.5-ml tube in a Tris-HCl buffer (10 mM; pH 7.2). Virus purification and DNA extraction were performed as for IIV6 in Drosophila flies (see above).

Infection of S. nonagrioides with AcMNPV
The AcMNPV-WP10 isolate (Chateigner et al. 2015) was used to infect ten fourth instar larvae of S. nonagrioides using the diet plug method (Sparks et al. 2008) (fig. 2). Each moth larva was fed %100,000 occlusion bodies (OBs) per 5-mm 3 diet plug. Upon host death, which occurred 2-5 days postinfection, OBs were first filtered through cheesecloth, purified twice by centrifugation (10 min at 4,500 Â g) with 0.1% sodium dodecyl sulfate, then distilled water, and finally resuspended in water. Approximately 1.5 Â 10 10 OBs were treated as described in Gilbert et al. (2014) to provide about 5 mg of high-quality dsDNA.

Infection of A. segetum Larvae with AgseGV
A previously sequenced AgseGV (isolate DA) was analyzed (Gueli Alletti, Eigenbrod, et al. 2017). Briefly, third to fourth instar A. segetum larvae were starved overnight and were fed the subsequent day with a small cube of artificial diet that contained 10 6 OBs of AgseGV (Gueli Alletti, Eigenbrod, et al. 2017). Larvae that consumed the entire piece within 12 h were transferred to virus-free diet and were checked daily for symptoms of viral infection. Larvae deceased by viral infection were collected and stored at À20 C until virus purification. OB purification and DNA extraction were performed as described previously (Gueli Alletti, Eigenbrod, et al. 2017).

Sequencing of Viral Genomes
Viral genomes were sequenced in three batches. The first batch included the AcMNPV sample purified from S. nonagrioides and the IIV6 samples purified from D. hydei and D. melanogaster flies, as well as from D. melanogaster S2 cells. One microgram of human DNA was added to the 2 mg of IIV6 purified from S2 cells DNA before sequencing to characterize artificial chimeras involving human TEs and the IIV6 genome (see below). These samples were sequenced by Novogen. A library was constructed for each sample with the NEBnext kit (average insert size was 350 bp), which was sequenced on a HiSeqX instrument in 2 Â 150 bp mode, generating 113.9-590.0 million paired-end reads. The second batch included the IIV6 samples purified from S. nonagrioides and from C. partellus, as well as the six IIV31 samples purified from three P. dilatatus and three A. vulgare individuals. These samples were sequenced by G enome Qu ebec. A library was constructed for each sample with the NEB ultra II kit (average insert size was 260 bp), which was sequenced on a HiSeqX instrument in 2 Â 150 bp paired-end mode, generating 246.1-330.9 million paired-end reads. The third batch included the 15 isolates of CpGV, the AgseGV-DA isolate, the AgseNPV-B PP2 isolate (¼ #0) and its five serial infections (#1, #3, #5, #7, and #10). Sequencing was performed on 50-100 ng purified viral DNA. Libraries were generated using a NexteraXT library preparation kit and sequencing was conducted with an Illumina NextSeq500 system (StarSEQ Ltd., Mainz, Germany) generating 0.3-10.9 million paired-end 150-bp-long reads (Gueli Alletti, Sauer, et al. 2017;Wennmann et al. 2020).

Detection of TEs Integrated in Viral Genomes
The aim of this study was to characterize as comprehensively as possible the diversity of host TEs that can become integrated into viral genomes. To this end, we searched for TEvirus junctions in sequencing reads using the approach developed by Gilbert et al. (2016). Reads were aligned separately on TE sequences and the viral genome using BlastN (-task megablast). Chimeric reads for which a portion aligned on a TE only and the other portion aligned on the viral genome only were then identified based on alignment coordinates. The R scripts used to detect and filter chimeric reads are provided in supplementary files S2 and S3, Supplementary Material online. Though artificial chimeras are certainly present in the various data sets analyzed here, several lines of evidence indicate that chimeras retained for counting and calculating the frequencies of TEs integrated into viral genomes are not artificial and result from bona fide integration that took place during replication of the virus. First, we considered that a TE was integrated into a viral genome only if we found at least three chimeric reads between a given TE and the virus (they were not required to map to the same position along the viral genome). Second, canonical transposition of a TE into a viral genome is expected to generate junction mapping in most cases at the very extremities of the TE sequences, that is, not anywhere along the TE sequence (Craig et al. 2015). Thus, if transposition into viral genomes occurred, chimeras covering the very extremities of the TE sequences should outnumber chimeras mapping internally to the TE sequences. This prediction was verified for all TEs ( fig.  3, supplementary fig. S2 and table S1, Supplementary Material online). Third, the addition of human DNA to the IIV6 DNA extracted from S2 cells as a control to estimate the presence of artificial chimeras allowed us to validate our approach by showing that contrary to Drosophila TEs, artificial chimeras involving human TEs mapped almost exclusively internally to human TE sequences (see below). Fourth, many TEs are known to duplicate a small target sequence motif upon integration called target site duplication (TSD). As in our earlier study (Gilbert et al. 2016), we were able to identify TSD motifs for several TEs found integrated into viral genomes ( fig. 3). Viral and TE sequences on which sequencing reads were aligned derive from WGSs and TE library that are fully described below. All moth TEs found integrated into viral genomes are provided in supplementary data 1, Supplementary Material online.

Viral Genomes Used for BLAST Searches
Genome sequences of four out of six viruses used in this study were retrieved from NCBI under accession number KR584663 (AgseGV-DA), KM102981 (AgseNPV-B), KM217575 (CpGV-M), and KM667940.1 (AcMNPV strain E2). For IIV31, we used the genome provided in the supplementary data of Loiseau et al. (2020). The IIV6 genome sequence available in NCBI was generated after replication on CF-124 cells of the C. fumiferana moth (Fukaya and Nasu 1966). As moths are distantly related to Drosophila flies, we reasoned that adaptation of IIV6 to S2 cells may have resulted in changes in its genome. We thus assembled the IIV6 genome de novo with the tadpole program (https://jgi.doe.gov/data-and-tools/bbtools/bbtools-user-guide/tadpole-guide/, version of December 2018, options used: "k ¼ 17"; "k ¼ 31"; "k ¼ 60"; "k ¼ 90" with "mincov ¼ 100"). The different assemblies obtained with 17, 31, 60, and 90 mers with 500Â and 1,500Â depth coverage were then fused with Geneious version 11.0.2 (https://www. geneious.com, last accessed on july 2020, options: de novo assembly, Geneious assembler, high sensitivity). The final assembly was annotated based on the available IIV6 (NC_003038.1) genome with the General Annotation Transfer Utility program (Tcherepanov et al. 2006).

TE Library
The TE library included all TE consensus sequences downloaded from Repbase in October 2017 (Bao et al. 2015), all those identified in 195 insect genomes by Peccoud et al. (2017) as well as all TEs identified integrated in the AcMNPV genome after replication of the virus onto T. ni and S. exigua moths (Gilbert et al. 2016). The library also Monitoring Insect TEs in dsDNA Viruses . doi:10.1093/molbev/msab198 MBE included all D. melanogaster TE sequences downloaded from Flybase and all A. vulgare TE consensus sequences characterized in Chebbi et al. (2019). In addition, our TE library was complemented with the TE consensus sequences constructed using RepeatModeler2 (Flynn et al. 2020) from the WGS of A. ipsilon (NCBI accession number: PNFC00000000.1), C. suppressalis (PNFC00000000.1), C. pomonella (QFTL00000000.2), D. hydei (QMEQ00000000.2), and S. nonagrioides (Muller et al. 2021). No WGS was available for A. segetum and C. partellus moths and the P. dilatatus isopod. To increase our chances to identify TEs in viral genomes, the sequencing reads obtained from A. segetum and C. partellus infections were also aligned on all lepidopteran WGS available in LepBase (as of May 17, 2018) as well as on all Noctuidae WGS and transcriptomes available in NCBI (as of October 15, 2017). In addition, for A. segetum, C. partellus, and P. dilatatus, we aligned reads on contigs that we de novo assembled using all sequencing reads not mapping onto viral genomes, that is, reads resulting from sequencing residual host and nonviral/ nonhost DNA that remained in the virus suspension after purification. Nonviral contig assembly was performed with the Tadpole assembler from the BBMap tools. Residual viral chunks were filtered out from the resulting contigs. All 248,482 TEs included in our library or characterized in partial assemblies of A. segetum, C. partellus, and P. dilatatus are provided in supplementary file S1, Supplementary Material online.

Inferring HTT between Insects
To assess whether moth TEs found integrated into viral genomes underwent HT during insect evolution, we used them as queries to perform BlastN searches (option megablast) on all 352 insect WGS available in GenBank as of March 2020. BlastN hits showing >90% nucleotide identity to the TE over >200 bp and >80% of its length were considered further. To assess whether such high levels of TE identity were due to vertical or horizontal transmission, we compared TE synonymous distances (dS) with the distribution of dS expected to occur under vertical transmission between each species pair of interest. The dS between TEs were calculated over the longest open reading frame found in the two copies involved in the best alignment for each species pair of interest. The distribution of dS expected under vertical transmission was generated for each species pair of interest by calculating dS between all single-copy and complete genes extracted from the WGS of each species using BUSCO version 4.0.2 (Simão et al. 2015;Waterhouse et al. 2018). All genes showing homology to the Insecta gene set (i.e., from 259 to 1,029 genes depending on species pairs) were used. MAFFT version 7.310 (Katoh and Standley 2013) was used to align nucleotide sequences of genes that are shared (orthologous) between species. Gene and TE alignments were analyzed in R version 3.6.1 (R Core Team 2019) with the "seqinr" package (Charif et al. 2005) to compute dS values.

Frequency of TEs in Viral Populations
The frequency of genomes carrying TE insertions was computed following the method described in Gilbert et al. (2016).
Briefly, we considered the mean number of insertions per virus genome to be the frequency of chimeras among viral reads divided by two-as one insertion yields two chimerasand multiplied by the ratio of the virus genome length over the length of a read. The latter was shortened to consider the minimum alignment length needed by BLAST to output a hit (see Gilbert et al. [2016] for details). The portion of viral genomes carrying insertions was computed as the probability that the number of inserted TEs into a viral genome exceeds zero, assuming that this number follows a Poisson distribution parameterized with the mean number of insertions per viral genome. To compute the frequencies of viral genomes carrying individual insertions in the AcMNPV-S. nonagrioides data set, another approach was used considering insertions distant by 5 bp or fewer as single insertion points and dividing the number of chimeric reads resulting from these insertions by the total number of reads covering the position (R script available in supplementary data 2, Supplementary Material online).

Supplementary Material
Supplementary data are available at Molecular Biology and Evolution online.