Ascoviruses are large, enveloped DNA viruses that induce remarkable changes in cellular architecture during which the cell is partitioned into numerous vesicles for viral replication. Previous studies have shown that these vesicles arise from a process resembling apoptosis yet which differs after nuclear lysis in that mitochondria are not degraded but are modified by the virus, changing in size, shape, and motility. Moreover, infection does not provoke an obvious innate immune response. Thus, we used in vivo RNA sequencing to determine whether infection by the Spodoptera frugiperda ascovirus 1a (SfAV-1a) modified expression of host mitochondrial, cytoskeletal, and innate immunity genes. We show that transcripts from many mitochondrial genes were similar to those from uninfected controls, whereas others increased slightly during vesicle formation, including those for ATP6, ATP8 synthase, and NADH dehydrogenase subunits, supporting electron microscopy (EM) data that these organelles were conserved for virus replication. Transcripts from 58 of 106 cytoskeletal genes studied increased or decreased more than 2-fold postinfection. More than half coded for mitochondrial motor proteins. Similar increases occurred for innate immunity transcripts and their negative regulators, including those for Toll, melanization, and phagocytosis pathways. However, those for many antimicrobial peptides, such as moricin, increased more than 20-fold. In addition, transcripts for gloverin-3, spod_x_tox, Hdd23, and lebocin, also antimicrobial, increased more than 20-fold. Interestingly, a phenoloxidase inhibitor transcript increased 12-fold, apparently to interfere with melanization. SfAV-1a destroys most fat body cells by 7 days postinfection, so innate immunity gene transcripts apparently occur in remaining cells in this tissue and possibly other major tissues, namely, epidermis and tracheal matrix.IMPORTANCE Ascoviruses are large DNA viruses that infect insects, inducing a cellular pathology that resembles apoptosis but which differs by causing enormous cellular hypertrophy followed by cleavage of the cell into numerous viral vesicles for replication. Previous EM studies suggest that mitochondria are important for vesicle formation. Transcriptome analyses of Spodoptera frugiperda larvae infected with SfAV-1a showed that mitochondrial transcripts were similar to those from uninfected controls or increased slightly during vesicle formation, especially for ATP6, ATP8 synthase, and NADH dehydrogenase subunits. This pattern resembles that for chronic disease-inducing viruses, which conserve mitochondria, differing markedly from viruses causing short-term viral diseases, which degrade mitochondrial DNA. Though mitochondrial transcript increases were low, our results demonstrate that SfAV-1a alters host mitochondrial expression more than any other virus. Regarding innate immunity, although SfAV-1a destroys most fat body cells, certain immunity genes were highly upregulated (greater than 20-fold), suggesting that these transcripts may originate from other tissues.

Ascoviruses are large dsDNA viruses characterized by the extraordinary changes they induce in cellular pathogenesis and architecture whereby after nuclear lysis and extensive hypertrophy, each cell is cleaved into numerous vesicles for virion reproduction. However, the level of viral replication and transcription in vesicles compared to other host tissues remains uncertain. Therefore, we applied RNA-Sequencing to compare the temporal transcriptome of Spodoptera frugiperda ascovirus (SfAV) and Trichoplusia ni ascovirus (TnAV) at 7, 14, and 21 days post-infection (dpi). We found most transcription occurred in viral vesicles, not in initial tissues infected, a remarkably novel reproduction mechanism compared to all other viruses and most other intracellular pathogens. Specifically, the highest level of viral gene expression occurred in hemolymph, for TnAV at 7 dpi, and SfAV at 14 dpi. Moreover, we found that host immune genes were partially down-regulated in hemolymph, where most viral replication occurred in highly dense accumulations of vesicles.

Ascoviruses are double-stranded DNA (dsDNA) viruses that attack caterpillars and differ from all other viruses by inducing nuclear lysis followed by cleavage of host cells into numerous anucleate vesicles in which virus replication continues as these grow in the blood. Ascoviruses are also unusual in that most encode a caspase or caspase-like proteins. A robust cell line to study the novel molecular biology of ascovirus replication in vitro is lacking. Therefore, we used strand-specific transcriptome sequencing (RNA-Seq) to study transcription in vivo in third instars of Spodoptera frugiperda infected with the type species, Spodoptera frugiperda ascovirus1a (SfAV-1a), sampling transcripts at different time points after infection. We targeted transcription of two types of SfAV-1a genes; first, 44 core genes that occur in several ascovirus species, and second, 26 genes predicted in silico to have metabolic functions likely involved in synthesizing viral vesicle membranes. Gene cluster analysis showed differences in temporal expression of SfAV-1a genes, enabling their assignment to three temporal classes: early, late, and very late. Inhibitors of apoptosis (IAP-like proteins; ORF016, ORF025, and ORF074) were expressed early, whereas its caspase (ORF073) was expressed very late, which correlated with apoptotic events leading to viral vesicle formation. Expression analysis revealed that a Diedel gene homolog (ORF121), the only known "virokine," was highly expressed, implying that this ascovirus protein helps evade innate host immunity. Lastly, single-nucleotide resolution of RNA-Seq data revealed 15 bicistronic and tricistronic messages along the genome, an unusual occurrence for large dsDNA viruses.IMPORTANCE Unlike all other DNA viruses, ascoviruses code for an executioner caspase, apparently involved in a novel cytopathology in which viral replication induces nuclear lysis followed by cell cleavage, yielding numerous large anucleate viral vesicles that continue to produce virions. Our transcriptome analysis of genome expression in vivo by the Spodoptera frugiperda ascovirus shows that inhibitors of apoptosis are expressed first, enabling viral replication to proceed, after which the SfAV-1a caspase is synthesized, leading to viral vesicle synthesis and subsequent extensive production of progeny virions. Moreover, we detected numerous bicistronic and tricistronic mRNA messages in the ascovirus transcriptome, implying that ascoviruses use other noncanonical translational mechanisms, such as internal ribosome entry sites (IRESs). These results provide the first insights into the molecular biology of a unique coordinated gene expression pattern in which cell architecture is markedly modified, more than in any other known eukaryotic virus, to promote viral reproduction and transmission.

Abstract Background The piRNA pathway has been shown in model organisms to be involved in silencing of transposons thereby providing genome stability. In D. melanogaster the majority of piRNAs map to these sequences. The medically important mosquito species Aedes aegypti has a large genome size, a high transposon load which includes Miniature Inverted repeat Transposable Elements (MITES) and an expansion of the piRNA biogenesis genes. Studies of transgenic lines of Ae. aegypti have indicated that introduced transposons are poorly remobilized and we sought to explore the basis of this. We wished to analyze the piRNA profile of Ae. aegypti and thereby determine if it is responsible for transposon silencing in this mosquito. Results Estimated piRNA sequence diversity was comparable between Ae. aegypti and D. melanogaster, but surprisingly only 19% of mosquito piRNAs mapped to transposons compared to 51% for D. melanogaster. Ae. aegypti piRNA clusters made up a larger percentage of the total genome than those of D. melanogaster but did not contain significantly higher percentages of transposon derived sequences than other regions of the genome. Ae. aegypti contains a number of protein coding genes that may be sources of piRNA biogenesis with two, traffic jam and maelstrom, implicated in this process in model organisms. Several genes of viral origin were also targeted by piRNAs. Examination of six mosquito libraries that had previously been transformed with transposon derived sequence revealed that new piRNA sequences had been generated to the transformed sequences, suggesting that they may have stimulated a transposon inactivation mechanism. Conclusions Ae. aegypti has a large piRNA complement that maps to transposons but primarily gene sequences, including many viral-derived sequences. This, together the more uniform distribution of piRNA clusters throughout its genome, suggest that some aspects of the piRNA system differ between Ae. aegypti and D. melanogaster.

Background

Anopheles stephensi is the key vector of malaria throughout the Indian subcontinent and Middle East and an emerging model for molecular and genetic studies of mosquito-parasite interactions. The type form of the species is responsible for the majority of urban malaria transmission across its range.

Results

Here, we report the genome sequence and annotation of the Indian strain of the type form of An. stephensi. The 221 Mb genome assembly represents more than 92% of the entire genome and was produced using a combination of 454, Illumina, and PacBio sequencing. Physical mapping assigned 62% of the genome onto chromosomes, enabling chromosome-based analysis. Comparisons between An. stephensi and An. gambiae reveal that the rate of gene order reshuffling on the X chromosome was three times higher than that on the autosomes. An. stephensi has more heterochromatin in pericentric regions but less repetitive DNA in chromosome arms than An. gambiae. We also identify a number of Y-chromosome contigs and BACs. Interspersed repeats constitute 7.1% of the assembled genome while LTR retrotransposons alone comprise more than 49% of the Y contigs. RNA-seq analyses provide new insights into mosquito innate immunity, development, and sexual dimorphism.

Conclusions

The genome analysis described in this manuscript provides a resource and platform for fundamental and translational research into a major urban malaria vector. Chromosome-based investigations provide unique perspectives on Anopheles chromosome evolution. RNA-seq analysis and studies of immunity genes offer new insights into mosquito biology and mosquito-parasite interactions.

Background

Anopheles stephensi is the key vector of malaria throughout the Indian subcontinent and Middle East and an emerging model for molecular and genetic studies of mosquito-parasite interactions. The type form of the species is responsible for the majority of urban malaria transmission across its range.

Results

Here, we report the genome sequence and annotation of the Indian strain of the type form of An. stephensi. The 221 Mb genome assembly represents more than 92% of the entire genome and was produced using a combination of 454, Illumina, and PacBio sequencing. Physical mapping assigned 62% of the genome onto chromosomes, enabling chromosome-based analysis. Comparisons between An. stephensi and An. gambiae reveal that the rate of gene order reshuffling on the X chromosome was three times higher than that on the autosomes. An. stephensi has more heterochromatin in pericentric regions but less repetitive DNA in chromosome arms than An. gambiae. We also identify a number of Y-chromosome contigs and BACs. Interspersed repeats constitute 7.1% of the assembled genome while LTR retrotransposons alone comprise more than 49% of the Y contigs. RNA-seq analyses provide new insights into mosquito innate immunity, development, and sexual dimorphism.

Conclusions

The genome analysis described in this manuscript provides a resource and platform for fundamental and translational research into a major urban malaria vector. Chromosome-based investigations provide unique perspectives on Anopheles chromosome evolution. RNA-seq analysis and studies of immunity genes offer new insights into mosquito biology and mosquito-parasite interactions.

Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. To investigate the genomic basis of vectorial capacity and explore new avenues for vector control, we sequenced the genomes of 16 anopheline mosquito species from diverse locations spanning ~100 million years of evolution. Comparative analyses show faster rates of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses, relative to Drosophila. Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts.

Tsetse flies are the sole vectors of human African trypanosomiasis throughout sub-Saharan Africa. Both sexes of adult tsetse feed exclusively on blood and contribute to disease transmission. Notable differences between tsetse and other disease vectors include obligate microbial symbioses, viviparous reproduction, and lactation. Here, we describe the sequence and annotation of the 366-megabase Glossina morsitans morsitans genome. Analysis of the genome and the 12,308 predicted protein-encoding genes led to multiple discoveries, including chromosomal integrations of bacterial (Wolbachia) genome sequences, a family of lactation-specific proteins, reduced complement of host pathogen recognition proteins, and reduced olfaction/chemosensory associated genes. These genome data provide a foundation for research into trypanosomiasis prevention and yield important insights with broad implications for multiple aspects of tsetse biology.