The genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest.

BACKGROUND: Although native to North America, the invasion of the aphid-like grape phylloxera Daktulosphaira vitifoliae across the globe altered the course of grape cultivation. For the past 150 years, viticulture relied on grafting-resistant North American Vitis species as rootstocks, thereby limiting genetic stocks tolerant to other stressors such as pathogens and climate change. Limited understanding of the insect genetics resulted in successive outbreaks across the globe when rootstocks failed. Here we report the 294-Mb genome of D. vitifoliae as a basic tool to understand host plant manipulation, nutritional endosymbiosis, and enhance global viticulture. RESULTS: Using a combination of genome, RNA, and population resequencing, we found grape phylloxera showed high duplication rates since its common ancestor with aphids, but similarity in most metabolic genes, despite lacking obligate nutritional symbioses and feeding from parenchyma. Similarly, no enrichment occurred in development genes in relation to viviparity. However, phylloxera evolved > 2700 unique genes that resemble putative effectors and are active during feeding. Population sequencing revealed the global invasion began from the upper Mississippi River in North America, spread to Europe and from there to the rest of the world. CONCLUSIONS: The grape phylloxera genome reveals genetic architecture relative to the evolution of nutritional endosymbiosis, viviparity, and herbivory. The extraordinary expansion in effector genes also suggests novel adaptations to plant feeding and how insects induce complex plant phenotypes, for instance galls. Finally, our understanding of the origin of this invasive species and its genome provide genetics resources to alleviate rootstock bottlenecks restricting the advancement of viticulture.

Functional genomics of parental care of insects.

Parental care was likely the first step most lineages made towards sociality. However, the molecular mechanisms that generate parental care are not broadly characterized. Insects are important as an evolutionary independent group from classic models of parental care, such as, house mice. They provide an opportunity to test the generality of our understanding. With this review, I survey the functional genomics of parental care of insects, summarize several recent advances in the broader framework for studying and understanding parental care, and finish with suggested priorities for further research. Although there are too few studies to draw definitive conclusions, I argue that natural selection appears to be rewiring existing gene networks to produce parental care, that the epigenetic mechanisms influencing parental care are not well understood, and, as an interesting early consensus, that genes strongly associated with carer/offspring interactions appear biased towards proteins that are secreted. I summarize the studies that have functionally validate candidate genes and highlight the increasing need to perform this work. I finish with arguments for both conceptual and practical changes moving forward. I argue that future work can increase the use of predictive frameworks, broaden its definition of conservation of mechanism to gene networks rather than single genes, and increase the use of more established comparative methods. I further highlight the practical considerations of standardizing analyses and reporting, increasing the sampling of both carers and offspring, better characterizing gene regulatory networks, better characterizing taxonomically restricted genes and any consistent role they have underpinning parental care, and using factorial designs to disentangle the influence of multiple variables on the expression of parental care.

Distributed physiology and the molecular basis of social life in eusocial insects.

The traditional focus of physiological and functional genomic research is on molecular processes that play out within a single multicellular organism. In the colonial (eusocial) insects such as ants, bees, and termites, molecular and behavioral responses of interacting nestmates are tightly linked, and key physiological processes are regulated at the scale of the colony. Such colony-level physiological processes regulate nestmate physiology in a distributed fashion, through various social communication mechanisms. As a result of physiological decentralization over evolutionary time, organismal mechanisms, for example related to pheromone detection, hormone signaling, and neural signaling pathways, are deployed in novel contexts to influence nestmate and colony traits. Here we explore how functional genomic, physiological, and behavioral studies can benefit from considering the traits of eusocial insects in this light. We highlight functional genomic work exploring how nestmate-level and colony-level traits arise and are influenced by interactions among physiologically-specialized nestmates of various developmental stages. We also consider similarities and differences between nestmate-level (organismal) and colony-level (superorganismal) physiological processes, and make specific hypotheses regarding the physiology of eusocial taxa. Integrating theoretical models of distributed systems with empirical functional genomics approaches will be useful in addressing fundamental questions related to the evolution of eusociality and collective behavior in natural systems.

Correlated variation and population differentiation in satellite DNA abundance among lines of Drosophila melanogaster.

Tandemly repeating satellite DNA elements in heterochromatin occupy a substantial portion of many eukaryotic genomes. Although often characterized as genomic parasites deleterious to the host, they also can be crucial for essential processes such as chromosome segregation. Adding to their interest, satellite DNA elements evolve at high rates; among Drosophila, closely related species often differ drastically in both the types and abundances of satellite repeats. However, due to technical challenges, the evolutionary mechanisms driving this rapid turnover remain unclear. Here we characterize natural variation in simple-sequence repeats of 2-10 bp from inbred Drosophila melanogaster lines derived from multiple populations, using a method we developed called k-Seek that analyzes unassembled Illumina sequence reads. In addition to quantifying all previously described satellite repeats, we identified many novel repeats of low to medium abundance. Many of the repeats show population differentiation, including two that are present in only some populations. Interestingly, the population structure inferred from overall satellite quantities does not recapitulate the expected population relationships based on the demographic history of D. melanogaster. We also find that some satellites of similar sequence composition are correlated across lines, revealing concerted evolution. Moreover, correlated satellites tend to be interspersed with each other, further suggesting that concerted change is partially driven by higher order structure. Surprisingly, we identified negative correlations among some satellites, suggesting antagonistic interactions. Our study demonstrates that current genome assemblies vastly underestimate the complexity, abundance, and variation of highly repetitive satellite DNA and presents approaches to understand their rapid evolutionary divergence.

The piRNA Pathway Guards the Germline Genome Against Transposable Elements.

Transposable elements (TEs) have the capacity to replicate and insert into new genomic locations. This contributs significantly to evolution of genomes, but can also result in DNA breaks and illegitimate recombination, and therefore poses a significant threat to genomic integrity. Excess damage to the germ cell genome results in sterility. A specific RNA silencing pathway, termed the piRNA pathway operates in germ cells of animals to control TE activity. At the core of the piRNA pathway is a ribonucleoprotein complex consisting of a small RNA, called piRNA, and a protein from the PIWI subfamily of Argonaute nucleases. The piRNA pathway relies on the specificity provided by the piRNA sequence to recognize complementary TE targets, while effector functions are provided by the PIWI protein. PIWI-piRNA complexes silence TEs both at the transcriptional level - by attracting repressive chromatin modifications to genomic targets - and at the posttranscriptional level - by cleaving TE transcripts in the cytoplasm. Impairment of the piRNA pathway leads to overexpression of TEs, significantly compromised genome structure and, invariably, germ cell death and sterility.The piRNA pathway is best understood in the fruit fly, Drosophila melanogaster, and in mouse. This Chapter gives an overview of current knowledge on piRNA biogenesis, and mechanistic details of both transcriptional and posttranscriptional TE silencing by the piRNA pathway. It further focuses on the importance of post-translational modifications and subcellular localization of the piRNA machinery. Finally, it provides a brief description of analogous pathways in other systems.

Ecological genomics of host shifts in Drosophila mojavensis.

Advances in next-generation sequencing technologies have liberated our dependency on model laboratory species for answering genomic and transcriptomic level questions. These new techniques have dramatically expanded our breadth of study organisms and have allowed the analysis of species from diverse ecological environments. One such species is the cactophilic Drosophila mojavensis that inhabits the deserts of western North America. These insects feed and develop in the necrotic cacti, feeding largely on the microflora of the necrotic plant tissues. Drosophila mojavensis is composed of four geographically and ecologically separated populations. Each population (Baja California peninsula, mainland Sonoran Desert, Mojave Desert and Santa Catalina Island) utilizes the necrotic tissues of distinct cactus species. The differences in the nutritional and chemical composition of the necroses include a set of toxic compounds to which resident population must adapt. These ecological differences have facilitated many of the life history, behavior, physiological and genetic differences between the cactus host populations. Genomic resources have allowed investigators to examine the genomic and transcriptional level changes associated with the local adaptation of the four D. mojavensis populations, thereby providing further understanding of the genetic mechanism of adaptation and its role in the divergence of ecologically distinct populations.

The genomics of an adaptive radiation: insights across the Heliconius speciation continuum.

Fueled by new technologies that allow rapid and inexpensive assessment of fine scale individual genomic variation, researchers are making transformational discoveries at the interface between genomes and biological complexity. Here we review genomic research in Heliconius butterflies - a radiation characterized by extraordinary phenotypic diversity in warningly colored wing patterns and composed of a continuum of taxa across the stages of speciation. These characteristics, coupled with a 50-year legacy of ecological and behavioral research, offer exceptional prospects for genomic studies into the nature of adaptive differences and the formation of new species. Research in Heliconius provides clear connections between genotype, phenotype, and fitness of wing color patterns shown to underlie adaptation and speciation. This research is challenging our perceptions about how speciation occurs in the presence of gene flow and the role of hybridization in generating adaptive novelty. With the release of the first Heliconius genome assembly, emerging genomic studies are painting a dynamic picture of the evolving species boundary. As the field of speciation genomics moves beyond describing patterns, towards a more integrated understanding of the process of speciation, groups such as Heliconius, where there is a clear speciation continuum and the traits underlying adaptation and speciation are known, will provide a roadmap for identifying variation crucial in the origins of biodiversity.

Sex pheromone desaturase functioning in a primitive Ostrinia moth is cryptically conserved in congeners' genomes.

(E)-11- and (Z)-11-tetradecenyl acetate are the most common female sex pheromone components in Ostrinia moths. The Δ11-desaturase expressed in the pheromone gland (PG) of female moths is a key enzyme that introduces a double bond into pheromone molecules. A single Δ11-desaturase of Ostrinia nubilalis, OnubZ/E11, has been shown to produce an ∼7:3 mixture of (E)-11- and (Z)-11-tetradecenoate from the substrate tetradecanoate. In contrast, the sex pheromone of Ostrinia latipennis, a primitive species of Ostrinia, is (E)-11-tetradecenol. This pheromone is unique in that it is not acetylated, and includes no Z isomer. In the present study, through the cloning and functional analysis of a PG-specific Δ11-desaturase in O. latipennis, we showed that the absence of the Z isomer in the pheromone is attributable to the strict product specificity of the Δ11-desaturase in this species, LATPG1. Phylogenetic analysis revealed that LATPG1 was not closely related to OnubZ/E11. Rather, it was closely related to retroposon-linked cryptic Δ11-desaturases (ezi-Δ11) found in the genomes of O. nubilalis and Ostrinia furnacalis. Taken together, the results showed that an unusual Δ11-desaturase is functionally expressed in O. latipennis, although the genes encoding this enzyme appear to be cryptic in congeners.

Genome-wide identification of neuropeptides and their receptor genes in Bemisia tabaci and their transcript accumulation change in response to temperature stresses.

Insect neuropeptides play an important role in regulating physiological functions such as growth, development, behavior and reproduction. We identified temperature-sensitive neuropeptides and receptor genes of the cotton whitefly, Bemisia tabaci. We identified 38 neuropeptide precursor genes and 35 neuropeptide receptors and constructed a phylogenetic tree using additional data from other insects. As temperature adaptability enables B. tabaci to colonize a diversity of habitats, we performed quantitative polymerase chain reaction with two temperature stresses (low = 4 °C and high = 40 °C) to screen for temperature-sensitive neuropeptides. We found many neuropeptides and receptors that may be involved in the temperature adaptability of B. tabaci. This study is the first to identify B. tabaci neuropeptides and their receptors, and it will help to reveal the roles of neuropeptides in temperature adaptation of B. tabaci.

Multiple SET methyltransferases are required to maintain normal heterochromatin domains in the genome of Drosophila melanogaster.

Methylation of histone H3 lysine 9 (H3K9) is a key feature of silent chromatin and plays an important role in stabilizing the interaction of heterochromatin protein 1 (HP1) with chromatin. Genomes of metazoans such as the fruit fly Drosophila melanogaster generally encode three types of H3K9-specific SET domain methyltransferases that contribute to chromatin homeostasis during the life cycle of the organism. SU(VAR)3-9, dG9a, and dSETDB1 all function in the generation of wild-type H3K9 methylation levels in the Drosophila genome. Two of these enzymes, dSETDB1 and SU(VAR)3-9, govern heterochromatin formation in distinct but overlapping patterns across the genome. H3K9 methylation in the small, heterochromatic fourth chromosome of D. melanogaster is governed mainly by dSETDB1, whereas dSETDB1 and SU(VAR)3-9 function in concert to methylate H3K9 in the pericentric heterochromatin of all chromosomes, with dG9a having little impact in these domains, as shown by monitoring position effect variegation. To understand how these distinct heterochromatin compartments may be differentiated, we examined the developmental timing of dSETDB1 function using a knockdown strategy. dSETDB1 acts to maintain heterochromatin during metamorphosis, at a later stage in development than the reported action of SU(VAR)3-9. Surprisingly, depletion of both of these enzymes has less deleterious effect than depletion of one. These results imply that dSETDB1 acts as a heterochromatin maintenance factor that may be required for the persistence of earlier developmental events normally governed by SU(VAR)3-9. In addition, the genetic interactions between dSETDB1 and Su(var)3-9 mutations emphasize the importance of maintaining the activities of these histone methyltransferases in balance for normal genome function.

Organization and evolution of heterochromatin in malaria mosquitoes.

The African malaria mosquito Anopheles gambiae was the first disease vector chosen for genome sequencing. Although its genome assembly has been facilitated by physical mapping, large gaps still pose a serious problem for accurate annotation and genome analysis. The majority of the gaps are located in regions of pericentromeric and intercalary heterochromatin. Genomic analysis has identified protein-coding genes and various classes of repetitive elements in the Anopheles heterochromatin. Molecular and cytogenetic studies have demonstrated that heterochromatin is a structurally heterogeneous and rapidly evolving part of the malaria mosquito genome.

[Molecular combing in studies of the genome organization and DNA replication].

Molecular combing (MC) yields preparations where individual DNA molecules are uniformly stretched and are parallel to each other. Fluorescence in situ hybridization on such preparations allows an exact mapping of DNA sequences, and pulsed inclusion of halogenated deoxyuridine analogs and their detection using fluorochrome-conjugated antibodies makes it possible to visualize replication. The MC technique was adapted for studying DNA replication in isolated Drosophila melanogaster organs, and it was checked whether a mutation of the Suppressor of UnderReplication (SuUR) gene directly affected the replication fork rate.

[Intercalary heterochromatin in the genome of Drosophila].

The modern concept of intercalary heterochromatin as polytene chromosome regions exhibiting a number of specific characteristics is formulated. DNA constituting these regions is replicated late in the S period; therefore, some strands of polytene chromosomes are underrepresented; i.e., they are underreplicated. Late-replicating regions account for about 7% of the genome; genes are located there in clusters of as many as 40. In general, the gene density in the clusters is substantially lower than in the main part of the genome. Late-replicating regions have an inactivating capacity: genes incorporated into these regions as parts of transposons are inactivated with a higher probability. These regions contain a specific protein SUUR affecting the rate of replication completion.

Large gene family expansion and variable selective pressures for cathepsin B in aphids.

Aphids exclusively feed on plant phloem sap that contains much sugar and some nonessential amino acids but is poor in lipids and proteins. Conventionally, it has been believed that aphids substantially have no intestinal digestion of proteins. However, we here report an unexpected finding that cysteine protease genes of the family cathepsin B are massively amplified in the lineage of aphids and that many of the protease genes exhibit gut-specific overexpression. By making use of expressed sequence tag data, sequenced cDNAs, and genomic trace sequences of the pea aphid Acyrthosiphon pisum, we identified a total of 28 cathepsin B-like gene copies in the genome of A. pisum. Phylogenetic analyses of all the cathepsin B genes in aphids revealed that genic expansion has continuously proceeded with basal, intermediary, and recent duplications. Estimation of molecular evolutionary rates indicated that major alterations of the rates often occurred after duplications. For example, a gene copy ("348") was shown to be slow evolving and close to genes of other insects like Drosophila melanogaster, whereas the other gene copies appeared to have evolved faster with higher ratios of nonsynonymous to synonymous substitutions. We identified a number of gene copies (16 in A. pisum) that contained a replacement at the site required for catalytic activity of the protease. Among these, 2 copies were pseudogenes, whereas the remaining copies were structurally intact and possibly acquired new functions. For example, a cluster of such gene copies ("1674") has been subjected to positive selection. Quantitative reverse transcriptase-polymerase chain reaction analyses revealed that the more conserved gene copy ("348") showed a constitutive expression, whereas 5 other forms ("84," "16," "16D," "1874," and "2744") were preferentially expressed in the gut of A. pisum. Putative biological roles of the diversified cathepsin B-like gene copies in aphids are discussed in relation to their nutritional physiology specialized for plant sap feeding lifestyle.

[On the phylogenetic position of Hexapoda within the Pancrustacea].

Contemporary views on the phylogeny of arthropods are at odds with the traditional system, which recognizes four independent arthropod classes: Chelicerata, Crustacea, Myriapoda and Insecta. There is compelling evidence that insects in fact comprise a monophyletic lineage with Crustacea within a larger clade of Pancrustacea (=Tetraconata). Which crustacean group is the closest living relative of insects remains an open question. Recent phylogenetic analyses based on multiple genes suggest their sistership with "lower" crustaceans, the Branchiopoda. This relationship was often impeached to be caused by the long branch attraction artifact. We analyzed concatenated data on 77 ribosomal proteins, elongation factor 1 alpha (EF1A), initiation factor 5 alpha (alF5A) and other selected nuclear and mitochondrial proteins. Nuclear protein data supports the monophyly of Hexapoda, the clade uniting entognath and ectognath insects. Hexapoda and Branchiopoda comprise a monophyletic lineage in most analyses. Maxillopoda occupies the sister position to the Hexapoda + Branchiopoda. "Higher" crustaceans, the Malacostraca, in most reconstructions comprise a more basal lineage withinthe Pancrustacea. Molecular synapomorphies in low homoplastic regions are found for the clades Hexapoda Branchiopoda + Maxillopoda and the monophyletic Malacostraca containing the Phyllocarida. Therefore, the sistership of Hexapoda and Branchiopoda and their position within Entomostraca may in fact represent bona fide phylogenetic relationships.

Genomic Basis of Circannual Rhythm in the European Corn Borer Moth.

Synchronizing the annual timing of physiological, morphological, and behavioral transitions with seasons enables survival in temperate environments [1]. The capacity to adjust life history timing and track local seasonal cycles can facilitate geographic expansion [2], adaptation [3], and tolerance [4-6] during rapid environmental change. Understanding the proximate causes of variation in seasonal timing improves prediction of future response and persistence [7, 8]. However, relatively little is known about the molecular basis generating this diversity [9], particularly in Lepidoptera, a group with many species in decline [10, 11]. In insects, the stress-tolerant physiological state of diapause enables coping with seasonal challenges [1, 12-15]. Seasonal changes in photoperiod and temperature are used to synchronize diapause with winter, and timing of diapause transitions varies widely within and among species [9, 16]. Changes in spring diapause termination in the European corn borer moth (Ostrinia nubilalis) have allowed populations to respond to shorter winters and emerge ∼3 weeks earlier in the year [17]. Multiple whole-genome approaches suggest two circadian clock genes, period (per) and pigment-dispersing factor receptor (Pdfr), underlie this polymorphism. Per and Pdfr are within interacting quantitative trait loci (QTL) and differ in allele frequency among individuals that end diapause early or late, with alleles maintained in high linkage disequilibrium. Our results provide testable hypotheses about the physiological role of circadian clock genes in the circannual timer. We predict these gene candidates will be targets of selection for future adaptation under continued global climate change [18].

Genome-Wide Association Analysis of Anoxia Tolerance in Drosophila melanogaster.

As the genetic bases to variation in anoxia tolerance are poorly understood, we used the Drosophila Genetics Reference Panel (DGRP) to conduct a genome-wide association study (GWAS) of anoxia tolerance in adult and larval Drosophila melanogaster Survival ranged from 0-100% in adults exposed to 6 h of anoxia and from 20-98% for larvae exposed to 1 h of anoxia. Anoxia tolerance had a broad-sense heritability of 0.552 in adults and 0.433 in larvae. Larval and adult phenotypes were weakly correlated but the anoxia tolerance of adult males and females were strongly correlated. The GWA identified 180 SNPs in adults and 32 SNPs in larvae associated with anoxia tolerance. Gene ontology enrichment analysis indicated that many of the 119 polymorphic genes associated with adult anoxia-tolerance were associated with ionic transport or immune function. In contrast, the 22 polymorphic genes associated with larval anoxia-tolerance were mostly associated with regulation of transcription and DNA replication. RNAi of mapped genes generally supported the hypothesis that disruption of these genes reduces anoxia tolerance. For two ion transport genes, we tested predicted directional and sex-specific effects of SNP alleles on adult anoxia tolerance and found strong support in one case but not the other. Correlating our phenotype to prior DGRP studies suggests that genes affecting anoxia tolerance also influence stress-resistance, immune function and ionic balance. Overall, our results provide evidence for multiple new potential genetic influences on anoxia tolerance and provide additional support for important roles of ion balance and immune processes in determining variation in anoxia tolerance.

Overexpression of Methoprene-tolerant, a Drosophila melanogaster gene that is critical for juvenile hormone action and insecticide resistance.

The Methoprene-tolerant (Met) gene of Drosophila melanogaster is involved in both juvenile hormone (JH) action and resistance to JH insecticides, such as methoprene. Although the consequences of Met mutations on development and methoprene resistance are known, no studies have examined Met+ overexpression. Met+ was overexpressed in transgenic lines with various promoters that drive overexpression to different levels. Flies expressing either genomic or cDNA Met+ transgenes showed higher susceptibility to both the morphogenetic and toxic effects of methoprene, consistent with the hormone-binding property of MET. Both the sensitive period and lethal period were the same as seen for non-overexpressing Met+ flies. However, continual exposure of high-overexpressing Met+ larvae to borderline-toxic or higher methoprene doses advanced the sensitive period from prepupae to first instar and the lethal period from pharate adults to larvae and early pupae. When expression of transgenic UAS-Met+ was driven to high levels by either an actin-GAL4 or tubulin-GAL4 promoter, larvae showed high mortality in the absence of methoprene, indicating that high MET titer is lethal, perhaps resulting from expression in an inappropriate tissue. Adults overexpressing Met+ did not show enhanced oogenesis, ruling out MET as a limiting factor for this hormone-driven physiology.

Genome-wide analysis of ionotropic receptor gene repertoire in Lepidoptera with an emphasis on its functions of Helicoverpa armigera.

The functions of the Ionotropic Receptor (IR) family have been well studied in Drosophila melanogaster, but only limited information is available in Lepidoptera. Here, we conducted a large-scale genome-wide analysis of the IR gene repertoire in 13 moths and 16 butterflies. Combining a homology-based approach and manual efforts, totally 996 IR candidates are identified including 31 pseudogenes and 825 full-length sequences, representing the most current comprehensive annotation in lepidopteran species. The phylogeny, expression and sequence characteristics classify Lepidoptera IRs into three sub-families: antennal IRs (A-IRs), divergent IRs (D-IRs) and Lepidoptera-specific IRs (LS-IRs), which is distinct from the case of Drosophila IRs. In comparison to LS-IRs and D-IRs, A-IRs members share a higher degree of protein identity and are distinguished into 16 orthologous groups in the phylogeny, showing conservation of gene structure. Analysis of selective forces on 27 orthologous groups reveals that these lepidopteran IRs have evolved under strong purifying selection (dN/dS≪1). Most notably, lineage-specific gene duplications that contribute primarily to gene number variations across Lepidoptera not only exist in D-IRs, but are present in the two other sub-families including members of IR41a, 76b, 87a, 100a and 100b. Expression profiling analysis reveals that over 80% (21/26) of Helicoverpa armigera A-IRs are expressed more highly in antennae of adults or larvae than other tissues, consistent with its proposed function in olfaction. However, some are also detected in taste organs like proboscises and legs. These results suggest that some A-IRs in H. armigera likely bear a dual function with their involvement in olfaction and gustation. Results from mating experiments show that two HarmIRs (IR1.2 and IR75d) expression is significantly up-regulated in antennae of mated female moths. However, no expression difference is observed between unmated female and male adults, suggesting an association with female host-searching behaviors. Our current study has greatly extended the IR gene repertoire resource in Lepidoptera, and more importantly, identifies potential IR candidates for olfactory, gustatory and oviposition behaviors in the cotton bollworm.

[Evolution from retrotransposons to retroviruses: origin of the env gene].

In genome of Drosophila melanogaster, various families of retrotransposons with different combination of functional domens and mechanisms of transposition are present. However only retrotransposons of gypsy family are retroviruses related to errantiviruses. Other families seemingly appeared as intermediate forms of retroviruses evolution. Despite the fact that the question on origin of retroviruses remains unclear, now the hypothesis of their origin from retrotransoposons can be considered the most consistent. Infectious properties of errantiviruses are linked to the presence of the third open reading frame (the env gene). Acquisition of the env gene conversed retrotransposons into retroviruses. So, origin of this gene is of special interest. Homologues of the env gene of errantiviruses are discovered in genomes of D. melanogaster, as well as in baculoviruses and in bacteria Wolbachia pipientis, the endosymbiont of Drosophila. It was shown that homologue of the env gene come to Wolbachia genome from Drosophila genome by horizontal transfer of the gypsy group retrotransposon. Thus, Wolbachia was not a donor of the env gene for errantiviruses. Seemingly, errantiviruses captured the baculoviral homologue of the env gene (f). However origin of the f gene is not clear. At the same time the env gene homologue in D. melanogaster genome exist (Iris). It must not be ruled out that the Iris gene was the source of the env gene of errantiviruses and baculoviruses.