Impact of Lateral Transfers on the Genomes of Lepidoptera

A novel and diverse family of filamentous DNA viruses associated with parasitic wasps

Large dsDNA viruses from the Naldaviricetes class are currently composed of four viral families infecting insects and/or crustaceans. Since the 1970s, particles described as filamentous viruses (FVs) have been observed by electronic microscopy in several species of Hymenoptera parasitoids but until recently, no genomic data was available. This study provides the first comparative morphological and genomic analysis of these FVs. We analyzed the genomes of seven FVs, six of which were newly obtained, to gain a better understanding of their evolutionary history. We show that these FVs share all genomic features of the Naldaviricetes while encoding five specific core genes that distinguish them from their closest relatives, the Hytrosaviruses. By mining public databases, we show that FVs preferentially infect Hymenoptera with parasitoid lifestyle and that these viruses have been repeatedly integrated into the genome of many insects, particularly Hymenoptera parasitoids, overall suggesting a long-standing specialization of these viruses to parasitic wasps. Finally, we propose a taxonomical revision of the class Naldaviricetes in which FVs related to the Leptopilina boulardi FV constitute a fifth family. We propose to name this new family, Filamentoviridae.

Identification and analysis of putative tRNA genes in baculovirus genomes

Transfer RNAs (tRNAs) genes are both coded for and arranged along some viral genomes representing the entire virosphere and seem to play different biological functions during infection, other than transferring the correct amino acid to a growing peptide chain. Baculovirus genome description and annotation has focused mostly on protein-coding genes, microRNA, and homologous regions. Here we carried out a large-scale in silico search for putative tRNA genes in baculovirus genomes. Ninety-six of 257 baculovirus genomes analyzed was found to contain at least one putative tRNA gene. We found great diversity in primary and secondary structure, in location within the genome, in intron presence and size, and in anti-codon identity. In some cases, genes of tRNA-containing genomes were found to have a bias for the codons specified by the tRNAs present in such genomes. Moreover, analysis revealed that most of the putative tRNA genes possessed conserved motifs for tRNA type 2 promoters, including the A-box and B-box motifs with few mismatches from the eukaryotic canonical motifs. From publicly available small RNA deep sequencing datasets of baculovirus-infected insect cells, we found evidence that a putative Autographa californica multiple nucleopolyhedrovirus Gln-tRNA gene was transcribed and modified with the addition of the non-templated 3'-CCA tail found at the end of all tRNAs. Further research is needed to determine the expression and functionality of these viral tRNAs.

Genome analysis of Psilogramma increta granulovirus and its intrapopulation diversity

The complete genome of Psilogramma increta granulovirus (PsinGV), isolated from P. increta (Lepidoptera: Sphingidae), was ultra-deep sequenced with a Novaseq PE150 platform and de novo assembled and annotated. The PsinGV genome is a circular double-stranded DNA, 103,721 bp in length, with a G+C content of 33.0%, the third lowest G+C content in present sequenced baculoviruses. It encodes 123 putative open reading frames, including 38 baculovirus core genes, 42 lepidopteran baculovirus conserved genes, 38 betabaculovirus conserved genes, and 5 genes unique to PsinGV. Meanwhile, 3 homologous repeated regions with the core sequence TTGCAA and 3 direct repeated sequences were identified within the PsinGV genome. Kimura two-parameters distance analysis confirmed that Psilogramma increta granulovirus is a representative of a prospective new species of the genus Betabaculovirus. Phylogenetic analysis based on the baculovirus core genes showed that PsinGV is closely related to Choristoneura fumiferana granulovirus, Clostera anastomosis granulovirus-B, and Erinnyis ello granulovirus. These four species therefore share a common ancestor residing in the Betabaculovirus genus. The genome of the PsinGV isolate contained two p10 copies: p10 and p10-2. Phylogenetic reconstruction of P10 suggests a transfer event of the p10-2 gene from an alphabaculovirus to the aforementioned common ancestor. Analysis of genomic diversity showed that 203 intrahost variants, including 182 single nucleotide variants and 21 short insertions/deletions, are present within the PsinGV isolate. Meanwhile, allele frequency indicated that the isolate contains three major genotypes, implying the archived isolate consists of several P. increta carcasses killed by PsinGV with different genotypes.

Investigating bracovirus chromosomal integration and inheritance in lepidopteran host and nontarget species

Bracoviruses (BVs) are domesticated viruses found in braconid parasitoid wasp genomes. They are composed of domesticated genes from a nudivrius, coding viral particles in which wasp DNA circles are packaged. BVs are viewed as possible vectors of horizontal transfer of genetic material (HT) from wasp to their hosts because they are injected, together with wasp eggs, by female wasps into their host larvae, and because they undergo massive chromosomal integration in multiple host tissues. Here, we show that chromosomal integrations of the Cotesia typhae BV (CtBV) persist up to the adult stage in individuals of its natural host, Sesamia nonagrioides, that survived parasitism. However, while reproducing host adults can bear an average of nearly two CtBV integrations per haploid genome, we were unable to retrieve any of these integrations in 500 of their offspring using Illumina sequencing. This suggests either that host gametes are less targeted by CtBVs than somatic cells or that gametes bearing BV integrations are nonfunctional. We further show that CtBV can massively integrate into the chromosomes of other lepidopteran species that are not normally targeted by the wasp in the wild, including one which is divergent by at least 100 million years from the natural host. Cell entry and chromosomal integration of BVs are thus unlikely to be major factors shaping wasp host range. Together, our results shed new light on the conditions under which BV-mediated wasp-to-host HT may occur and provide information that may be helpful to evaluate the potential risks of uncontrolled HT associated with the use of parasitoid wasps as biocontrol agents.

Horizontal gene transfer from genetically modified plants - Regulatory considerations

Gene technology regulators receive applications seeking permission for the environmental release of genetically modified (GM) plants, many of which possess beneficial traits such as improved production, enhanced nutrition and resistance to drought, pests and diseases. The regulators must assess the risks to human and animal health and to the environment from releasing these GM plants. One such consideration, of many, is the likelihood and potential consequence of the introduced or modified DNA being transferred to other organisms, including people. While such gene transfer is most likely to occur to sexually compatible relatives (vertical gene transfer), horizontal gene transfer (HGT), which is the acquisition of genetic material that has not been inherited from a parent, is also a possibility considered during these assessments. Advances in HGT detection, aided by next generation sequencing, have demonstrated that HGT occurrence may have been previously underestimated. In this review, we provide updated evidence on the likelihood, factors and the barriers for the introduced or modified DNA in GM plants to be horizontally transferred into a variety of recipients. We present the legislation and frameworks the Australian Gene Technology Regulator adheres to with respect to the consideration of risks posed by HGT. Such a perspective may generally be applicable to regulators in other jurisdictions as well as to commercial and research organisations who develop GM plants.

Diversity and Functional Roles of the Gut Microbiota in Lepidopteran Insects

Lepidopteran insects are one of the most widespread and speciose lineages on Earth, with many common pests and beneficial insect species. The evolutionary success of their diversification depends on the essential functions of gut microorganisms. This diverse gut microbiota of lepidopteran insects provides benefits in nutrition and reproductive regulation and plays an important role in the defence against pathogens, enhancing host immune homeostasis. In addition, gut symbionts have shown promising applications in the development of novel tools for biological control, biodegradation of waste, and blocking the transmission of insect-borne diseases. Even though most microbial symbionts are unculturable, the rapidly expanding catalogue of microbial genomes and the application of modern genetic techniques offer a viable alternative for studying these microbes. Here, we discuss the gut structure and microbial diversity of lepidopteran insects, as well as advances in the understanding of symbiotic relationships and interactions between hosts and symbionts. Furthermore, we provide an overview of the function of the gut microbiota, including in host nutrition and metabolism, immune defence, and potential mechanisms of detoxification. Due to the relevance of lepidopteran pests in agricultural production, it can be expected that the research on the interactions between lepidopteran insects and their gut microbiota will be used for biological pest control and protection of beneficial insects in the future.

Horizontally transmitted parasitoid killing factor shapes insect defense to parasitoids

Interkingdom competition occurs between hymenopteran parasitoids and insect viruses sharing the same insect hosts. It has been assumed that parasitoid larvae die with the death of the infected host or as result of competition for host resources. Here we describe a gene family, parasitoid killing factor (pkf), that encodes proteins toxic to parasitoids of the Microgastrinae group and determines parasitism success. Pkfs are found in several entomopathogenic DNA virus families and in some lepidopteran genomes. We provide evidence of equivalent and specific toxicity against endoparasites for PKFs found in entomopoxvirus, ascovirus, baculovirus, and Lepidoptera through a mechanism that elicits apoptosis in the cells of susceptible parasitoids. This highlights the evolutionary arms race between parasitoids, viruses, and their insect hosts.

Discovery of Hyperactive Antifreeze Protein from Phylogenetically Distant Beetles Questions Its Evolutionary Origin

Beetle hyperactive antifreeze protein (AFP) has a unique ability to maintain a supercooling state of its body fluids, however, less is known about its origination. Here, we found that a popular stag beetle Dorcus hopei binodulosus (Dhb) synthesizes at least 6 isoforms of hyperactive AFP (DhbAFP). Cold-acclimated Dhb larvae tolerated -5 °C chilled storage for 24 h and fully recovered after warming, suggesting that DhbAFP facilitates overwintering of this beetle. A DhbAFP isoform (~10 kDa) appeared to consist of 6-8 tandem repeats of a 12-residue consensus sequence (TCTxSxNCxxAx), which exhibited 3 °C of high freezing point depression and the ability of binding to an entire surface of a single ice crystal. Significantly, these properties as well as DNA sequences including the untranslated region, signal peptide region, and an AFP-encoding region of Dhb are highly similar to those identified for a known hyperactive AFP (TmAFP) from the beetle Tenebrio molitor (Tm). Progenitor of Dhb and Tm was branched off approximately 300 million years ago, so no known evolution mechanism hardly explains the retainment of the DNA sequence for such a lo-ng divergence period. Existence of unrevealed gene transfer mechanism will be hypothesized between these two phylogenetically distant beetles to acquire this type of hyperactive AFP.

Detecting bracoviral orthologs distribution in five tsetse fly species and the housefly genomes

Objective

Mutualism between endogenous viruses and eukaryotes is still poorly understood. Several endogenous double-stranded polydnaviruses, bracoviruses, homologous to those present in parasitic braconid wasp genomes were detected in the tsetse fly (Glossina morsitans morsitans). This is peculiar since tsetse flies do not share a reproductive lifestyle similar to wasps, but deliver fully developed larvae that pupate within minutes of exiting their mothers. The objective of this study is to investigate genomic distribution of bracoviral sequences in five tsetse fly species and the housefly, and examine its value as a potential vector control strategy target point. We use comparative genomics to determine the presence, distribution across Glossina species genomes, and evolutionary relationships of bracoviruses of five tsetse fly species and the housefly.

Results

We report on homologous bracoviruses in multiple Dipteran genomes. Phylogenetic reconstruction using within-species concatenated bracoviral orthologs shows great congruence with previously reconstructed insect species phylogenies. Our findings suggest that bracoviruses present in Diptera originate from a single integration event of the viral genome that occurred in an ancestor insect before the evolutionary radiation of different insect orders.

Functional characterization of the horizontally transferred 4,5-DOPA extradiol dioxygenase gene in the domestic silkworm, Bombyx mori

4,5-DOPA dioxygenase (DODA) is a crucial enzyme in the biosynthetic pathway of betalain. Previous studies have shown that DODA is present in plants, fungi and bacteria. Using updated data, here we show that DODA genes (BmDODA) in the domestic silkworm (Bombyx mori) and other lepidopteran insects are most likely to be horizontally transferred from fungi. A synteny analysis indicated that BmDODA1 is orthologous to other lepidopteran DODAs and that BmDODA2 is a paralogous gene. To explore the function of DODA in Lepidoptera, we first examined the expression patterns of BmDODA1. BmDODA1 showed high transcriptional and translational levels in the midgut and head. Then, we exogenously expressed the BmDODA1 gene, detected 4,5-DOPA ring-cleaving activity and calculated the kinetic parameters of the recombinant BmDODA1. We found that the transcription levels of BmDODA1 were significantly induced by the pathogens Bacillus bombyseptieus and Escherichia coli. Thus, the horizontal transfer of the BmDODA gene in the silkworm may be involved in dopa metabolism and contribute to antimicrobial activity in this species. Our results provide a documented example of functional horizontal gene transfer (HGT) between fungi and animals and expand our knowledge of HGT amongst eukaryotes.

Diversification of Transposable Elements in Arthropods and Its Impact on Genome Evolution

Transposable elements (TEs) are ubiquitous in arthropods. However, analyses of large-scale and long-term coevolution between TEs and host genomes remain scarce in arthropods. Here, we choose 14 representative Arthropoda species from eight orders spanning more than 500 million years of evolution. By developing an unbiased TE annotation pipeline, we obtained 87 to 2266 TE reference sequences in a species, which is a considerable improvement compared to the reference TEs previously annotated in Repbase. We find that TE loads are diversified among species and were previously underestimated. The highly species- and time-specific expansions and contractions, and intraspecific sequence diversification are the leading driver of long terminal repeat (LTR) dynamics in Lepidoptera. Terminal inverted repeats (TIRs) proliferated substantially in five species with large genomes. A phylogenetic comparison reveals that the loads of multiple TE subfamilies are positively correlated with genome sizes. We also identified a few horizontally transferred TE candidates across nine species. In addition, we set up the Arthropod Transposable Elements database (ArTEdb) to provide TE references and annotations. Collectively, our results provide high-quality TE references and uncover that TE loads and expansion histories vary greatly among arthropods, which implies that TEs are an important driving force shaping the evolution of genomes through gain and loss.

Co-option of microbial associates by insects and their impact on plant-folivore interactions

Plants possess a suite of traits that make them challenging to consume by insect herbivores. Plant tissues are recalcitrant, have low levels of protein, and may be well defended by chemicals. Insects use diverse strategies for overcoming these barriers, including co-opting metabolic activities from microbial associates. In this review, we discuss the co-option of bacteria and fungi in the herbivore gut. We particularly focus upon chewing, folivorous insects (Coleoptera and Lepidoptera) and discuss the impacts of microbial co-option on herbivore performance and plant responses. We suggest that there are two components to microbial co-option: fixed and plastic relationships. Fixed relationships are involved in integral dietary functions and can be performed by microbial enzymes co-opted into the genome or by stably transferred associates. In contrast, the majority of gut symbionts appear to be looser and perform more facultative, context-dependent functions. This more plastic, variable co-option of bacteria likely produces a greater number of insect phenotypes, which interact differently with plant hosts. By altering plant detection of herbivory or mediating insect interactions with plant defensive compounds, microbes can effectively improve herbivore performance in real time within and between generations.

Global survey of mobile DNA horizontal transfer in arthropods reveals Lepidoptera as a prime hotspot

More than any other genome components, Transposable Elements (TEs) have the capacity to move across species barriers through Horizontal Transfer (HT), with substantial evolutionary consequences. Previous large-scale surveys, based on full-genomes comparisons, have revealed the transposition mode as an important predictor of HT rates variation across TE superfamilies. However, host biology could represent another major explanatory factor, one that needs to be investigated through extensive taxonomic sampling. Here we test this hypothesis using a field collection of 460 arthropod species from Tahiti and surrounding islands. Through targeted massive parallel sequencing, we uncover patterns of HT in three widely-distributed TE superfamilies with contrasted modes of transposition. In line with earlier findings, the DNA transposons under study (TC1-Mariner) were found to transfer horizontally at the highest frequency, closely followed by the LTR superfamily (Copia), in contrast with the non-LTR superfamily (Jockey), that mostly diversifies through vertical inheritance and persists longer within genomes. Strikingly, across all superfamilies, we observe a marked excess of HTs in Lepidoptera, an insect order that also commonly hosts baculoviruses, known for their ability to transport host TEs. These results turn the spotlight on baculoviruses as major potential vectors of TEs in arthropods, and further emphasize the importance of non-vertical TE inheritance in genome evolution.

Transposable elements are chunks of DNA that can produce copies of themselves. New copies usually insert in the genome of their carrier but are occasionally subject to horizontal transmission between organisms, sometimes belonging to evolutionarily-distant lineages. Previous surveys have established that the probability of such events is largely conditioned by the transposition mechanism. For example, elements with an RNA intermediate tend to be less frequently involved in horizontal transfers. Here we investigate host taxa as another potential explanatory factor of variation in horizontal transfer rates. Using targeted sequencing in hundreds of insects and other arthropod species collected in South Pacific islands, we found that butterflies and moths (Lepidoptera) show an abnormally elevated rate of horizontal transfers. Previous studies have established that Lepidoptera are also commonly attacked by baculoviruses, large viruses that can transport host DNA. Taken together, these findings point to baculoviruses as a major suspect for transposable elements transfers across arthropod species.

Cotesia congregata Bracovirus Circles Encoding PTP and Ankyrin Genes Integrate into the DNA of Parasitized Manduca sexta Hemocytes

Polydnaviruses (PDVs) are essential for the parasitism success of tens of thousands of species of parasitoid wasps. PDVs are present in wasp genomes as proviruses, which serve as the template for the production of double-stranded circular viral DNA carrying virulence genes that are injected into lepidopteran hosts. PDV circles do not contain genes coding for particle production, thereby impeding viral replication in caterpillar hosts during parasitism. Here, we investigated the fate of PDV circles of Cotesia congregata bracovirus during parasitism of the tobacco hornworm, Manduca sexta, by the wasp Cotesia congregata Sequences sharing similarities with host integration motifs (HIMs) of Microplitis demolitor bracovirus (MdBV) circles involved in integration into DNA could be identified in 12 CcBV circles, which encode PTP and VANK gene families involved in host immune disruption. A PCR approach performed on a subset of these circles indicated that they persisted in parasitized M. sexta hemocytes as linear forms, possibly integrated in host DNA. Furthermore, by using a primer extension capture method based on these HIMs and high-throughput sequencing, we could show that 8 out of 9 circles tested were integrated in M. sexta hemocyte genomic DNA and that integration had occurred specifically using the HIM, indicating that an HIM-mediated specific mechanism was involved in their integration. Investigation of BV circle insertion sites at the genome scale revealed that certain genomic regions appeared to be enriched in BV insertions, but no specific M. sexta target site could be identified.IMPORTANCE The identification of a specific and efficient integration mechanism shared by several bracovirus species opens the question of its role in braconid parasitoid wasp parasitism success. Indeed, results obtained here show massive integration of bracovirus DNA in somatic immune cells at each parasitism event of a caterpillar host. Given that bracoviruses do not replicate in infected cells, integration of viral sequences in host DNA might allow the production of PTP and VANK virulence proteins within newly dividing cells of caterpillar hosts that continue to develop during parasitism. Furthermore, this integration process could serve as a basis to understand how PDVs mediate the recently identified gene flux between parasitoid wasps and Lepidoptera and the frequency of these horizontal transfer events in nature.

Horizontal acquisition of transposable elements and viral sequences: patterns and consequences

It is becoming clear that most eukaryotic transposable elements (TEs) owe their evolutionary success in part to horizontal transfer events, which enable them to invade new species. Recent large-scale studies are beginning to unravel the mechanisms and ecological factors underlying this mode of transmission. Viruses are increasingly recognized as vectors in the process but also as a direct source of genetic material horizontally acquired by eukaryotic organisms. Because TEs and endogenous viruses are major catalysts of variation and innovation in genomes, we argue that horizontal inheritance has had a more profound impact in eukaryotic evolution than is commonly appreciated. To support this proposal, we compile a list of examples, including some previously unrecognized, whereby new host functions and phenotypes can be directly attributed to horizontally acquired TE or viral sequences. We predict that the number of examples will rapidly grow in the future as the prevalence of horizontal transfer in the life cycle of TEs becomes even more apparent, firmly establishing this form of non-Mendelian inheritance as a consequential facet of eukaryotic evolution.