Broad surveys of DNA viral diversity obtained through viral metagenomics of mosquitoes

Vector-Enabled Metagenomic (VEM) Surveys Using Whiteflies (Aleyrodidae) Reveal Novel Begomovirus Species in the New and Old Worlds

Whitefly-transmitted viruses belonging to the genus Begomovirus (family Geminiviridae) represent a substantial threat to agricultural food production. The rapid evolutionary potential of these single-stranded DNA viruses combined with the polyphagous feeding behavior of their whitefly vector (Bemisia tabaci) can lead to the emergence of damaging viral strains. Therefore, it is crucial to characterize begomoviruses circulating in different regions and crops globally. This study utilized vector-enabled metagenomics (VEM) coupled with high-throughput sequencing to survey begomoviruses directly from whiteflies collected in various locations (California (USA), Guatemala, Israel, Puerto Rico, and Spain). Begomoviruses were detected in all locations, with the highest diversity identified in Guatemala where up to seven different species were identified in a single field. Both bipartite and monopartite viruses were detected, including seven new begomovirus species from Guatemala, Puerto Rico, and Spain. This begomovirus survey extends the known diversity of these highly damaging plant viruses. However, the new genomes described here and in the recent literature appear to reflect the outcome of interactions between closely-related species, often resulting from recombination, instead of unique, highly divergent species.

Genome Sequences of Poaceae-Associated Gemycircularviruses from the Pacific Ocean Island of Tonga

We sampled and analyzed 43 Poaceae plants from the Pacific Ocean island of Tonga for the presence of circular DNA viruses. From these samples, we recovered three gemycircularvirus genomes, which share >99% identity, from Brachiaria deflexa (n = 2) and sugarcane (n = 1). These genomes share <61% genome-wide identity with other gemycircularvirus sequences in public databases.

A novel gemycircularvirus from experimental rats

Recently, gemycircularviruses have been found in humans and various species of animals. Here, a novel gemycircularvirus named Ch-zjrat-01 from blood samples of experimental rats was characterized. The novel gemycircularvirus encodes two major proteins, including a capsid protein (Cap) and a replication-associated protein (Rep). Phylogenetic analysis based on the amino acid sequence of Rep indicated that Ch-zjrat-01 clusters with two gemycircularviruses discovered from bird (KF371635) and mosquito (HQ335086), sharing 48.7 and 49.4 % sequence identities with them, respectively.

Novel circular single-stranded DNA viruses identified in marine invertebrates reveal high sequence diversity and consistent predicted intrinsic disorder patterns within putative structural proteins

Viral metagenomics has recently revealed the ubiquitous and diverse nature of single-stranded DNA (ssDNA) viruses that encode a conserved replication initiator protein (Rep) in the marine environment. Although eukaryotic circular Rep-encoding ssDNA (CRESS-DNA) viruses were originally thought to only infect plants and vertebrates, recent studies have identified these viruses in a number of invertebrates. To further explore CRESS-DNA viruses in the marine environment, this study surveyed CRESS-DNA viruses in various marine invertebrate species. A total of 27 novel CRESS-DNA genomes, with Reps that share less than 60.1% identity with previously reported viruses, were recovered from 21 invertebrate species, mainly crustaceans. Phylogenetic analysis based on the Rep revealed a novel clade of CRESS-DNA viruses that included approximately one third of the marine invertebrate associated viruses identified here and whose members may represent a novel family. Investigation of putative capsid proteins (Cap) encoded within the eukaryotic CRESS-DNA viral genomes from this study and those in GenBank demonstrated conserved patterns of predicted intrinsically disordered regions (IDRs), which can be used to complement similarity-based searches to identify divergent structural proteins within novel genomes. Overall, this study expands our knowledge of CRESS-DNA viruses associated with invertebrates and explores a new tool to evaluate divergent structural proteins encoded by these viruses.

A metagenomic survey of viral abundance and diversity in mosquitoes from Hubei province

Mosquitoes as one of the most common but important vectors have the potential to transmit or acquire a lot of viruses through biting, however viral flora in mosquitoes and its impact on mosquito-borne disease transmission has not been well investigated and evaluated. In this study, the metagenomic techniquehas been successfully employed in analyzing the abundance and diversity of viral community in three mosquito samples from Hubei, China. Among 92,304 reads produced through a run with 454 GS FLX system, 39% have high similarities with viral sequences belonging to identified bacterial, fungal, animal, plant and insect viruses, and 0.02% were classed into unidentified viral sequences, demonstrating high abundance and diversity of viruses in mosquitoes. Furthermore, two novel viruses in subfamily Densovirinae and family Dicistroviridae were identified, and six torque tenosus virus1 in family Anelloviridae, three porcine parvoviruses in subfamily Parvovirinae and a Culex tritaeniorhynchus rhabdovirus in Family Rhabdoviridae were preliminarily characterized. The viral metagenomic analysis offered us a deep insight into the viral population of mosquito which played an important role in viral initiative or passive transmission and evolution during the process.

Small circular single stranded DNA viral genomes in unexplained cases of human encephalitis, diarrhea, and in untreated sewage

Viruses with small circular ssDNA genomes encoding a replication initiator protein can infect a wide range of eukaryotic organisms ranging from mammals to fungi. The genomes of two such viruses, a cyclovirus (CyCV-SL) and gemycircularvirus (GemyCV-SL) were detected by deep sequencing of the cerebrospinal fluids of Sri Lankan patients with unexplained encephalitis. One and three out of 201 CSF samples (1.5%) from unexplained encephalitis patients tested by PCR were CyCV-SL and GemyCV-SL DNA positive respectively. Nucleotide similarity searches of pre-existing metagenomics datasets revealed closely related genomes in feces from unexplained cases of diarrhea from Nicaragua and Brazil and in untreated sewage from Nepal. Whether the tropism of the cyclovirus and gemycircularvirus reported here include humans or other cellular sources in or on the human body remains to be determined.

Factors responsible for the emergence of arboviruses; strategies, challenges and limitations for their control

Slave trading of Africans to the Americas, during the 16th to the 19th century was responsible for the first recorded emergence in the New World of two arthropod-borne viruses (arboviruses), yellow fever virus and dengue virus. Many other arboviruses have since emerged from their sylvatic reservoirs and dispersed globally due to evolving factors that include anthropological behaviour, commercial transportation and land-remediation. Here, we outline some characteristics of these highly divergent arboviruses, including the variety of life cycles they have developed and the mechanisms by which they have adapted to evolving changes in habitat and host availability. We cite recent examples of virus emergence that exemplify how arboviruses have exploited the consequences of the modern human lifestyle. Using our current understanding of these viruses, we also attempt to demonstrate some of the limitations encountered in developing control strategies to reduce the impact of future emerging arbovirus diseases. Finally, we present recommendations for development by an international panel of experts reporting directly to World Health Organization, with the intention of providing internationally acceptable guidelines for improving emerging arbovirus disease control strategies. Success in these aims should alleviate the suffering and costs encountered during recent decades when arboviruses have emerged from their sylvatic environment.

RNA shotgun metagenomic sequencing of northern California (USA) mosquitoes uncovers viruses, bacteria, and fungi

Mosquitoes, most often recognized for the microbial agents of disease they may carry, harbor diverse microbial communities that include viruses, bacteria, and fungi, collectively called the microbiota. The composition of the microbiota can directly and indirectly affect disease transmission through microbial interactions that could be revealed by its characterization in natural populations of mosquitoes. Furthermore, the use of shotgun metagenomic sequencing (SMS) approaches could allow the discovery of unknown members of the microbiota. In this study, we use RNA SMS to characterize the microbiota of seven individual mosquitoes (species include Culex pipiens, Culiseta incidens, and Ochlerotatus sierrensis) collected from a variety of habitats in California, USA. Sequencing was performed on the Illumina HiSeq platform and the resulting sequences were quality-checked and assembled into contigs using the A5 pipeline. Sequences related to single stranded RNA viruses of the Bunyaviridae and Rhabdoviridae were uncovered, along with an unclassified genus of double-stranded RNA viruses. Phylogenetic analysis finds that in all three cases, the closest relatives of the identified viral sequences are other mosquito-associated viruses, suggesting widespread host-group specificity among disparate viral taxa. Interestingly, we identified a Narnavirus of fungi, also reported elsewhere in mosquitoes, that potentially demonstrates a nested host-parasite association between virus, fungi, and mosquito. Sequences related to 8 bacterial families and 13 fungal families were found across the seven samples. Bacillus and Escherichia/Shigella were identified in all samples and Wolbachia was identified in all Cx. pipiens samples, while no single fungal genus was found in more than two samples. This study exemplifies the utility of RNA SMS in the characterization of the natural microbiota of mosquitoes and, in particular, the value of identifying all microbes associated with a specific host.

Metagenomic profile of the viral communities in Rhipicephalus spp. ticks from Yunnan, China

Besides mosquitoes, ticks are regarded as the primary source of vector-borne infectious diseases. Indeed, a wide variety of severe infectious human diseases, including those involving viruses, are transmitted by ticks in many parts of the world. To date, there are no published reports on the use of next-generation sequencing for studying viral diversity in ticks or discovering new viruses in these arthropods from China. Here, Ion-torrent sequencing was used to investigate the presence of viruses in three Rhipicephalus spp. tick pools (NY-11, NY-13, and MM-13) collected from the Menglian district of Yunnan, China. The sequencing run resulted in 3,641,088, 3,106,733, and 3,871,851 reads in each tick pool after trimming. Reads and assembled contiguous sequences (contigs) were subject to basic local alignment search tool analysis against the GenBank database. Large numbers of reads and contigs related to known viral sequences corresponding to a broad range of viral families were identified. Some of the sequences originated from viruses that have not been described previously in ticks. Our findings will facilitate better understanding of the tick virome, and add to our current knowledge of disease-causing viruses in ticks living under natural conditions.

Xenosurveillance: a novel mosquito-based approach for examining the human-pathogen landscape

BACKGROUND

Globally, regions at the highest risk for emerging infectious diseases are often the ones with the fewest resources. As a result, implementing sustainable infectious disease surveillance systems in these regions is challenging. The cost of these programs and difficulties associated with collecting, storing and transporting relevant samples have hindered them in the regions where they are most needed. Therefore, we tested the sensitivity and feasibility of a novel surveillance technique called xenosurveillance. This approach utilizes the host feeding preferences and behaviors of Anopheles gambiae, which are highly anthropophilic and rest indoors after feeding, to sample viruses in human beings. We hypothesized that mosquito bloodmeals could be used to detect vertebrate viral pathogens within realistic field collection timeframes and clinically relevant concentrations.

METHODOLOGY/PRINCIPAL FINDINGS

To validate this approach, we examined variables influencing virus detection such as the duration between mosquito blood feeding and mosquito processing, the pathogen nucleic acid stability in the mosquito gut and the pathogen load present in the host's blood at the time of bloodmeal ingestion using our laboratory model. Our findings revealed that viral nucleic acids, at clinically relevant concentrations, could be detected from engorged mosquitoes for up to 24 hours post feeding by qRT-PCR. Subsequently, we tested this approach in the field by examining blood from engorged mosquitoes from two field sites in Liberia. Using next-generation sequencing and PCR we were able to detect the genetic signatures of multiple viral pathogens including Epstein-Barr virus and canine distemper virus.

CONCLUSIONS/SIGNIFICANCE

Together, these data demonstrate the feasibility of xenosurveillance and in doing so validated a simple and non-invasive surveillance tool that could be used to complement current biosurveillance efforts.

Identification of novel Bromus- and Trifolium-associated circular DNA viruses

The genomes of a large number of highly diverse novel circular DNA viruses from a wide range of sources have been characterised in recent years, including circular single-stranded DNA (ssDNA) viruses that share similarities with plant-infecting ssDNA viruses of the family Geminiviridae. Here, we describe six novel circular DNA viral genomes that encode replication-associated (Rep) proteins that are most closely related to those of either geminiviruses or gemycircularviruses (a new group of ssDNA viruses that are closely related to geminiviruses). Four possible viral genomes were recovered from Bromus hordeaceus sampled in New Zealand, and two were recovered from B. hordeaceus and Trifolium resupinatum sampled in France. Two of the viral genomes from New Zealand (one from the North Island and one from the South Island each) share >99 % sequence identity, and two genomes recovered from B. hordeaceus and T. resupinatum sampled in France share 74 % identity. All of the viral genomes that were recovered were found to have a major open reading frame on both their complementary and virion-sense strands, one of which likely encodes a Rep and the other a capsid protein. Although future infectivity studies are needed to identify the host range of these viruses, this is the first report of circular DNA viruses associated with grasses in New Zealand.

An overview of malaria transmission from the perspective of Amazon Anopheles vectors

In the Americas, areas with a high risk of malaria transmission are mainly located in the Amazon Forest, which extends across nine countries. One keystone step to understanding the Plasmodium life cycle in Anopheles species from the Amazon Region is to obtain experimentally infected mosquito vectors. Several attempts to colonise Anopheles species have been conducted, but with only short-lived success or no success at all. In this review, we review the literature on malaria transmission from the perspective of its Amazon vectors. Currently, it is possible to develop experimental Plasmodium vivax infection of the colonised and field-captured vectors in laboratories located close to Amazonian endemic areas. We are also reviewing studies related to the immune response to P. vivax infection of Anopheles aquasalis, a coastal mosquito species. Finally, we discuss the importance of the modulation of Plasmodium infection by the vector microbiota and also consider the anopheline genomes. The establishment of experimental mosquito infections with Plasmodium falciparum, Plasmodium yoelii and Plasmodium berghei parasites that could provide interesting models for studying malaria in the Amazonian scenario is important. Understanding the molecular mechanisms involved in the development of the parasites in New World vectors is crucial in order to better determine the interaction process and vectorial competence.

Characterisation of a diverse range of circular replication-associated protein encoding DNA viruses recovered from a sewage treatment oxidation pond

Our knowledge of circular replication-associated protein encoding single-stranded (CRESS) DNA virus diversity has increased dramatically in recent years, largely due to advances in high-throughput sequencing technologies. These viruses are apparently major virome components in most terrestrial and aquatic environments and it is therefore of interest to determine their diversity at the interfaces between these environments. Treated sewage water is a particularly interesting interface between terrestrial and aquatic viromes in that it is directly pumped into waterways and is likely to contain virus populations that have been strongly impacted by humans. We used a combination of high-throughput sequencing, full genome PCR amplification, cloning and Sanger sequencing to investigate the diversity of CRESS DNA viruses present in a sewage oxidation pond. Using this approach, we recovered 50 putatively complete novel CRESS viral genomes (it remains possible that some are components of multipartite viral genomes) and 11 putatively sub-genome-length circular DNA molecules which may be either defective genomes or components of multipartite genomes. Thirteen of the genomes have bidirectional genome organisations and share similar conserved replication-associated protein (Rep) motifs to those of the gemycircularviruses: a group that in turn is most closely related to the geminiviruses. The remaining 37 viral genomes share very low degrees of Rep similarity to those of all other known CRESS DNA viruses. This number of highly divergent CRESS DNA virus genomes within a single sewage treatment pond further reinforces the notion that there likely exist hundreds of completely unknown genus/family level CRESS DNA virus groupings.

Bionomics of temperate and tropical Culicoides midges: knowledge gaps and consequences for transmission of Culicoides-borne viruses

Culicoides midges are abundant hematophagous flies that vector arboviruses of veterinary and medical importance. Dramatic changes in the epidemiology of Culicoides-borne arboviruses have occurred since 1998, including the emergence of exotic viruses in northern temperate regions, increases in global disease incidence, and enhanced virus diversity in tropical zones. Drivers may include changes in climate, land use, trade, and animal husbandry. New Culicoides species and new wild reservoir hosts have been implicated in transmission, highlighting the dynamic nature of pathogen-vector-host interactions. Focusing on potential vector species worldwide and key elements of vectorial capacity, we review the sensitivity of Culicoides life cycles to abiotic and biotic factors. We consider implications for designing control measures and understanding impacts of environmental change in different ecological contexts. Critical geographical, biological, and taxonomic knowledge gaps are prioritized. Recent developments in genomics and mathematical modeling may enhance ecological understanding of these complex arbovirus systems.

Characterization of a novel circo-like virus in Aedes vexans mosquitoes from Germany: evidence for a new genus within the family Circoviridae

Over recent decades, metagenomic studies have expanded the number of newly described, often unclassified, viruses within the family Circoviridae. Using broad-spectrum circovirus and cyclovirus PCRs, we characterized a novel circo-like virus in Aedes vexans mosquitoes from Germany whose main putative ORFs shared very low amino acid identity with those of previously characterized circoviruses and cycloviruses. Phylogenetic and genetic distance analysis revealed that this new virus species defined, together with previously described mosquito- and bat faeces-derived circo-like viruses, a different genus, tentatively called Krikovirus, within the family Circoviridae. We further demonstrated that viruses of the putative genus Krikovirus all shared a genomic organization that was unique among the family Circoviridae. Further investigations are needed to determine the host range, tissue tropism and transmission route(s). This report increases the current knowledge of the genetic diversity and evolution of the members of the family Circoviridae.

Species-specific viromes in the ancestral holobiont Hydra

Recent evidence showing host specificity of colonizing bacteria supports the view that multicellular organisms are holobionts comprised of the macroscopic host in synergistic interdependence with a heterogeneous and host-specific microbial community. Whereas host-bacteria interactions have been extensively investigated, comparatively little is known about host-virus interactions and viral contribution to the holobiont. We sought to determine the viral communities associating with different Hydra species, whether these viral communities were altered with environmental stress, and whether these viruses affect the Hydra-associated holobiont. Here we show that each species of Hydra harbors a diverse host-associated virome. Primary viral families associated with Hydra are Myoviridae, Siphoviridae, Inoviridae, and Herpesviridae. Most Hydra-associated viruses are bacteriophages, a reflection of their involvement in the holobiont. Changes in environmental conditions alter the associated virome, increase viral diversity, and affect the metabolism of the holobiont. The specificity and dynamics of the virome point to potential viral involvement in regulating microbial associations in the Hydra holobiont. While viruses are generally regarded as pathogenic agents, our study suggests an evolutionary conserved ability of viruses to function as holobiont regulators and, therefore, constitutes an emerging paradigm shift in host-microbe interactions.

Anopheles gambiae densovirus (AgDNV) has negligible effects on adult survival and transcriptome of its mosquito host

Mosquito densoviruses (DNVs) are candidate agents for paratransgenic control of malaria and other vector-borne diseases. Unlike other mosquito DNVs, the Anopheles gambiae DNV (AgDNV) is non-pathogenic to larval mosquitoes. However, the cost of infection upon adults and the molecular mechanisms underpinning infection in the mosquito host are unknown. Using life table analysis, we show that AgDNV infection has minimal effects on An. gambiae survival (no significant effect in 2 replicates and a slight 2 day survival decrease in the third replicate). Using microarrays, we show that AgDNV has very minimal effect on the adult mosquito transcriptome, with only 4-15 genes differentially regulated depending on the statistical criteria imposed. The minimal impact upon global transcription provides some mechanistic understanding of lack of virus pathogenicity, suggesting a long co-evolutionary history that has shifted towards avirulence. From an applied standpoint, lack of strong induced fitness costs makes AgDNV an attractive agent for paratransgenic malaria control.

Metagenomic shotgun sequencing of a Bunyavirus in wild-caught Aedes aegypti from Thailand informs the evolutionary and genomic history of the Phleboviruses

Arthropod-borne viruses significantly impact human health. They span multiple families, all of which include viruses not known to cause disease. Characterizing these representatives could provide insights into the origins of their disease-causing counterparts. Field-caught Aedes aegypti mosquitoes from Nakhon Nayok, Thailand, underwent metagenomic shotgun sequencing to reveal a Bunyavirus closely related to Phasi Charoen (PhaV) virus, isolated in 2009 from Ae. aegypti near Bangkok. Phylogenetic analysis of this virus suggests it is basal to the Phlebovirus genus thus making it ideally positioned phylogenetically for understanding the evolution of these clinically important viruses. Genomic analysis finds that a gene necessary for virulence in vertebrates, but not essential for viral replication in arthropods, is missing. The sequencing of this phylogenetically-notable and genomically-unique Phlebovirus from wild mosquitoes exemplifies the utility and efficacy of metagenomic shotgun sequencing for virus characterization in arthropod vectors of human diseases.

Metagenomics, paratransgenesis and the Anopheles microbiome: a portrait of the geographical distribution of the anopheline microbiota based on a meta-analysis of reported taxa

Anophelines harbour a diverse microbial consortium that may represent an extended gene pool for the host. The proposed effects of the insect microbiota span physiological, metabolic and immune processes. Here we synthesise how current metagenomic tools combined with classical culture-dependent techniques provide new insights in the elucidation of the role of the Anopheles-associated microbiota. Many proposed malaria control strategies have been based upon the immunomodulating effects that the bacterial components of the microbiota appear to exert and their ability to express anti-Plasmodium peptides. The number of identified bacterial taxa has increased in the current "omics" era and the available data are mostly scattered or in "tables" that are difficult to exploit. Published microbiota reports for multiple anopheline species were compiled in an Excel® spreadsheet. We then filtered the microbiota data using a continent-oriented criterion and generated a visual correlation showing the exclusive and shared bacterial genera among four continents. The data suggested the existence of a core group of bacteria associated in a stable manner with their anopheline hosts. However, the lack of data from Neotropical vectors may reduce the possibility of defining the core microbiota and understanding the mosquito-bacteria interactive consortium.

Viral metagenomics on animals as a tool for the detection of zoonoses prior to human infection?

Many human viral infections have a zoonotic, i.e., wild or domestic animal, origin. Several zoonotic viruses are transmitted to humans directly via contact with an animal or indirectly via exposure to the urine or feces of infected animals or the bite of a bloodsucking arthropod. If a virus is able to adapt and replicate in its new human host, human-to-human transmissions may occur, possibly resulting in an epidemic, such as the A/H1N1 flu pandemic in 2009. Thus, predicting emerging zoonotic infections is an important challenge for public health officials in the coming decades. The recent development of viral metagenomics, i.e., the characterization of the complete viral diversity isolated from an organism or an environment using high-throughput sequencing technologies, is promising for the surveillance of such diseases and can be accomplished by analyzing the viromes of selected animals and arthropods that are closely in contact with humans. In this review, we summarize our current knowledge of viral diversity within such animals (in particular blood-feeding arthropods, wildlife and domestic animals) using metagenomics and present its possible future application for the surveillance of zoonotic and arboviral diseases.

Testing the Hamilton-Zuk hypothesis: past, present, and future

Hamilton and Zuk proposed a good-genes model of sexual selection in which genetic variation can be maintained when females prefer ornaments that indicate resistance to parasites. When trait expression depends on a male's resistance, the co-adaptive cycles between host resistance and parasite virulence provide a mechanism in which genetic variation for fitness is continually renewed. The model made predictions at both the intraspecific and interspecific levels. In the three decades since its publication, these predictions have been theoretically examined in models of varying complexity, and empirically tested across many vertebrate and invertebrate taxa. Despite such prolonged interest, however, it has turned out to be extremely difficult to empirically demonstrate the process described, in part because we have not been able to test the underlying mechanisms that would unequivocally identify how parasites act as mediators of sexual selection. Here, we discuss how the use of high-throughput sequencing datasets available from modern genomic approaches might improve our ability to test this model. We expect that important contributions will come through the ability to identify and quantify the suite of parasites likely to influence the evolution of hosts' resistance, to confidently reconstruct phylogenies of both host and parasite taxa, and, perhaps most exciting, to detect generational cycles of heritable variants in populations of hosts and parasites. Integrative approaches, building on systems undergoing parasite-mediated selection with genomic resources already available, will be particularly useful in moving toward robust tests of this hypothesis. We finish by presenting case studies of well-studied host-parasite relationships that represent promising avenues for future research.

Diverse small circular single-stranded DNA viruses identified in a freshwater pond on the McMurdo Ice Shelf (Antarctica)

Antarctica has some of the harshest environmental conditions for existence of life on Earth. In this pilot study we recovered eight diverse circular single-stranded DNA (ssDNA) viral genome sequences (1904-3120 nts) from benthic mats dominated by filamentous cyanobacteria in a freshwater pond on the McMurdo Ice Shelf sampled in 1988. All genomes contain two to three major open reading frames (ORFs) that are uni- or bi-directionally transcribed and all have an ORF encoding a replication-associated protein (Rep). In one genome, the second ORF has similarity to a capsid protein (CP) of Nepavirus which is most closely related to geminiviruses. Additionally, all genomes have two intergenic regions that contain putative stem loop structures, six genomes have NANTATTAC as the nonanucleotide motif, while one has CCTTATTAC, and another has a non-canonical stem loop. In the large intergenic region, we identified iterative sequences flanking the putative stem-loop elements which are a hallmark of most circular ssDNA viruses encoding rolling circle replication (RCR) initiators of the HUH endonuclease superfamily. The Reps encoded by ssDNA viral genomes recovered in this study shared <38% pairwise identity to all other Reps of known ssDNA viruses. A previous study on Lake Limnopolar (Livingston Island, South Shetland Islands), using next-generation sequencing identified circular ssDNA viruses and their putative Reps share <35% pairwise identity to those from the viral genomes removed in this study. It is evident from our pilot study that the global diversity of ssDNA viruses is grossly underestimated and there is limited knowledge on ssDNA viruses in Antarctica.

Metavir 2: new tools for viral metagenome comparison and assembled virome analysis

Background

Metagenomics, based on culture-independent sequencing, is a well-fitted approach to provide insights into the composition, structure and dynamics of environmental viral communities. Following recent advances in sequencing technologies, new challenges arise for existing bioinformatic tools dedicated to viral metagenome (i.e. virome) analysis as (i) the number of viromes is rapidly growing and (ii) large genomic fragments can now be obtained by assembling the huge amount of sequence data generated for each metagenome.

Results

To face these challenges, a new version of Metavir was developed. First, all Metavir tools have been adapted to support comparative analysis of viromes in order to improve the analysis of multiple datasets. In addition to the sequence comparison previously provided, viromes can now be compared through their k-mer frequencies, their taxonomic compositions, recruitment plots and phylogenetic trees containing sequences from different datasets. Second, a new section has been specifically designed to handle assembled viromes made of thousands of large genomic fragments (i.e. contigs). This section includes an annotation pipeline for uploaded viral contigs (gene prediction, similarity search against reference viral genomes and protein domains) and an extensive comparison between contigs and reference genomes. Contigs and their annotations can be explored on the website through specifically developed dynamic genomic maps and interactive networks.

Conclusions

The new features of Metavir 2 allow users to explore and analyze viromes composed of raw reads or assembled fragments through a set of adapted tools and a user-friendly interface.

Whole genomic characterization of a Korean human parechovirus type 1 (HPeV1) identifies recombination events

A human parechovirus (HPeV) strain CAU10-NN was detected from a stool sample of a 2-year-old healthy female infant in South Korea using a metagenomic approach. The CAU10-NN virus was isolated using a cell culture system and its whole genome was analyzed. The RNA genome of the CAU10-NN strain consists of 7,348 nucleotides (nt), excluding a poly(A) tail. A large open reading frame of 6,540 nt that encodes a putative polyprotein precursor of 2,180 amino acids is flanked by a 5'-untranslated region (UTR) of 708 nt and 3'-UTR of 88 nt followed by a poly(A) tail. Phylogenetic analysis shows that the CAU10-NN strain belongs to HPeV1. SimPlot and Bootscan analyses reveal that the virus genome is composed of regions related to corresponding genomic regions of other HPeVs. Recombination analysis indicates that the CAU10-NN strain might be a product of more than one genomic recombination event that occurred among HPeV1, HPeV3, and HPeV4 strains.

RNA viral metagenome of whiteflies leads to the discovery and characterization of a whitefly-transmitted carlavirus in North America

Whiteflies from the Bemisia tabaci species complex have the ability to transmit a large number of plant viruses and are some of the most detrimental pests in agriculture. Although whiteflies are known to transmit both DNA and RNA viruses, most of the diversity has been recorded for the former, specifically for the Begomovirus genus. This study investigated the total diversity of DNA and RNA viruses found in whiteflies collected from a single site in Florida to evaluate if there are additional, previously undetected viral types within the B. tabaci vector. Metagenomic analysis of viral DNA extracted from the whiteflies only resulted in the detection of begomoviruses. In contrast, whiteflies contained sequences similar to RNA viruses from divergent groups, with a diversity that extends beyond currently described viruses. The metagenomic analysis of whiteflies also led to the first report of a whitefly-transmitted RNA virus similar to Cowpea mild mottle virus (CpMMV Florida) (genus Carlavirus) in North America. Further investigation resulted in the detection of CpMMV Florida in native and cultivated plants growing near the original field site of whitefly collection and determination of its experimental host range. Analysis of complete CpMMV Florida genomes recovered from whiteflies and plants suggests that the current classification criteria for carlaviruses need to be reevaluated. Overall, metagenomic analysis supports that DNA plant viruses carried by B. tabaci are dominated by begomoviruses, whereas significantly less is known about RNA viruses present in this damaging insect vector.

New insights into mycoviruses and exploration for the biological control of crop fungal diseases

Mycoviruses are viruses that infect fungi. A growing number of novel mycoviruses have expanded our knowledge of virology, particularly in taxonomy, ecology, and evolution. Recent progress in the study of mycoviruses has comprehensively improved our understanding of the properties of mycoviruses and has strengthened our confidence to explore hypovirulence-associated mycoviruses that control crop diseases. In this review, the advantages of using hypovirulence-associated mycoviruses to control crop diseases are discussed, and, as an example, the potential for Sclerotinia sclerotiorum hypovirulence-associated DNA virus 1 (SsHADV-1) to control the stem rot of rapeseed (Brassica napus) is also introduced. Fungal vegetative incompatibility is likely to be the key factor that limits the wide utilization of mycoviruses to control crop diseases; however, there are suggested strategies for resolving this problem.

Assessment of viral community functional potential from viral metagenomes may be hampered by contamination with cellular sequences

Although the importance of viruses in natural ecosystems is widely acknowledged, the functional potential of viral communities is yet to be determined. Viral genomes are traditionally believed to carry only those genes that are directly pertinent to the viral life cycle, though this view was challenged by the discovery of metabolism genes in several phage genomes. Metagenomic approaches extended these analyses to a community scale, and several studies concluded that microbial and viral communities encompass similar functional potentials. However, these conclusions could originate from the presence of cellular DNA within viral metagenomes. We developed a computational method to estimate the proportion and origin of cellular sequences in a set of 67 published viromes. A quarter of the datasets were found to contain a substantial amount of sequences originating from cellular genomes. When considering only viromes with no cellular DNA detected, the functional potential of viral and microbial communities was found to be fundamentally different-a conclusion more consistent with the actual picture drawn from known viruses. Yet a significant number of cellular metabolism genes was still retrieved in these viromes, suggesting that the presence of auxiliary genes involved in various metabolic pathways within viral genomes is a general trend in the virosphere.

Functional Immunomics of the Squash Bug, Anasa tristis (De Geer) (Heteroptera: Coreidae)

The Squash bug, Anasa tristis (De Geer), is a major piercing/sucking pest of cucurbits, causing extensive damage to plants and fruits, and transmitting phytopathogens. No genomic resources to facilitate field and laboratory studies of this pest were available; therefore the first de novo exome for this destructive pest was assembled. RNA was extracted from insects challenged with bacterial and fungal immunoelicitors, insects fed on different cucurbit species, and insects from all life stages from egg to adult. All treatments and replicates were separately barcoded for subsequent analyses, then pooled for sequencing in a single lane using the Illumina HiSeq2000 platform. Over 211 million 100-base tags generated in this manner were trimmed, filtered, and cleaned, then assembled into a de novo reference transcriptome using the Broad Institute Trinity assembly algorithm. The assembly was annotated using NCBIx NR, BLAST2GO, KEGG and other databases. Of the >130,000 total assemblies 37,327 were annotated identifying the sequences of candidate gene silencing targets from immune, endocrine, reproductive, cuticle, and other physiological systems. Expression profiling of the adult immune response was accomplished by aligning the 100-base tags from each biological replicate from each treatment and controls to the annotated reference assembly of the A. tristis transcriptome.

Novel virus discovery and genome reconstruction from field RNA samples reveals highly divergent viruses in dipteran hosts

We investigated whether small RNA (sRNA) sequenced from field-collected mosquitoes and chironomids (Diptera) can be used as a proxy signature of viral prevalence within a range of species and viral groups, using sRNAs sequenced from wild-caught specimens, to inform total RNA deep sequencing of samples of particular interest. Using this strategy, we sequenced from adult Anopheles maculipennis s.l. mosquitoes the apparently nearly complete genome of one previously undescribed virus related to chronic bee paralysis virus, and, from a pool of Ochlerotatus caspius and Oc. detritus mosquitoes, a nearly complete entomobirnavirus genome. We also reconstructed long sequences (1503-6557 nt) related to at least nine other viruses. Crucially, several of the sequences detected were reconstructed from host organisms highly divergent from those in which related viruses have been previously isolated or discovered. It is clear that viral transmission and maintenance cycles in nature are likely to be significantly more complex and taxonomically diverse than previously expected.

The possible impact of persistent virus infection on the function of the RNAi machinery in insects: a hypothesis

RNAi experiments in insects are characterized by great variability in efficiency; for instance beetles and locusts are very amenable to dsRNA-mediated gene silencing, while other insect groups, most notably lepidopterans, are more refractory to RNAi. Several factors can be forwarded that could affect the efficiency of RNAi, such as the composition and function of the intracellular RNAi machinery, the mechanism of dsRNA uptake, the presence of dsRNA- and siRNA-degrading enzymes and non-specific activation of the innate immune response. In this essay, we investigate the evidence whether persistent infection with RNA viruses could be a major factor that affects the response to exogenous dsRNA in insects. The occurrence of RNA viruses in different insect groups will be discussed, as well as several mechanisms by which viruses could interfere with the process of RNAi. Finally, the impact of RNA virus infection on the design of dsRNA-based insect control strategies will be considered.

Enhanced arbovirus surveillance with deep sequencing: Identification of novel rhabdoviruses and bunyaviruses in Australian mosquitoes

Viral metagenomics characterizes known and identifies unknown viruses based on sequence similarities to any previously sequenced viral genomes. A metagenomics approach was used to identify virus sequences in Australian mosquitoes causing cytopathic effects in inoculated mammalian cell cultures. Sequence comparisons revealed strains of Liao Ning virus (Reovirus, Seadornavirus), previously detected only in China, livestock-infecting Stretch Lagoon virus (Reovirus, Orbivirus), two novel dimarhabdoviruses, named Beaumont and North Creek viruses, and two novel orthobunyaviruses, named Murrumbidgee and Salt Ash viruses. The novel virus proteomes diverged by ≥ 50% relative to their closest previously genetically characterized viral relatives. Deep sequencing also generated genomes of Warrego and Wallal viruses, orbiviruses linked to kangaroo blindness, whose genomes had not been fully characterized. This study highlights viral metagenomics in concert with traditional arbovirus surveillance to characterize known and new arboviruses in field-collected mosquitoes. Follow-up epidemiological studies are required to determine whether the novel viruses infect humans.

Identification and characterisation of a highly divergent geminivirus: evolutionary and taxonomic implications

During a large scale "non a priori" survey in 2010 of South African plant-infecting single stranded DNA viruses, a highly divergent geminivirus genome was isolated from a wild spurge, Euphorbia caput-medusae. In addition to being infectious in E. caput-medusae, the cloned viral genome was also infectious in tomato and Nicotiana benthamiana. The virus, named Euphorbia caput-medusae latent virus (EcmLV) due to the absence of infection symptoms displayed by its natural host, caused severe symptoms in both tomato and N. benthamiana. The genome organisation of EcmLV is unique amongst geminiviruses and it likely expresses at least two proteins without any detectable homologues within public sequence databases. Although clearly a geminivirus, EcmLV is so divergent that we propose its placement within a new genus that we have tentatively named Capulavirus. Using a set of highly divergent geminiviruses genomes, it is apparent that recombination has likely been a primary process in the genus-level diversification of geminiviruses. It is also demonstrated how this insight, taken together with phylogenetic analyses of predicted coat protein and replication associated protein (Rep) amino acid sequences indicate that the most recent common ancestor of the geminiviruses was likely a dicot-infecting virus that, like modern day mastreviruses and becurtoviruses, expressed its Rep from a spliced complementary strand transcript.

Two novel circo-like viruses detected in human feces: complete genome sequencing and electron microscopy analysis

The application of viral metagenomic techniques and a series of PCRs in a human fecal sample enabled the detection of two novel circular unisense DNA viral genomes with 92% nucleotide similarity. The viruses were tentatively named circo-like virus-Brazil (CLV-BR) strains hs1 and hs2 and have genome lengths of 2526 and 2533 nucleotides, respectively. Four major open reading frames (ORFs) were identified in each of the genomes, and differences between the two genomes were primarily observed in ORF 2. Only ORF 3 showed significant amino acid similarities to a putative rolling circle replication initiator protein (Rep), although with low identity (36%). Our phylogenetic analysis, based on the Rep protein, demonstrated that the CLV-BRs do not cluster with members of the Circoviridae, Nanoviridae or Geminiviridae families and are more closely related to circo-like genomes previously identified in reclaimed water and feces of a wild rodent and of a bat. The CLV-BRs are members of a putative new family of circular Rep-encoding ssDNA viruses. Electron microscopy revealed icosahedral (~23 nm) structures, likely reflecting the novel viruses, and rod-shaped viral particles (~65-460 × 21 × 10 nm in length, diameter, and axial canal, respectively). Circo-like viruses have been detected in stool samples from humans and other mammals (bats, rodents, chimpanzees and bovines), cerebrospinal fluid and sera from humans, as well as samples from many other sources, e.g., insects, meat and the environment. Further studies are needed to classify all novel circular DNA viruses and elucidate their hosts, pathogenicity and evolutionary history.

Evaluation of rapid and simple techniques for the enrichment of viruses prior to metagenomic virus discovery

•

The effect of simple virus enrichment methods were tested on a metagenomics dataset.

•

Centrifugation, filtration or nuclease-treatment was evaluated.

•

A multi-step enrichment method increased the proportion of virus sequences.

•

This evaluation guides researchers in their choice of enrichment methodology.

The discovery of new or divergent viruses using metagenomics and high-throughput sequencing has become more commonplace. The preparation of a sample is known to have an effect on the representation of virus sequences within the metagenomic dataset yet comparatively little attention has been given to this. Physical enrichment techniques are often applied to samples to increase the number of viral sequences and therefore enhance the probability of detection. With the exception of virus ecology studies, there is a paucity of information available to researchers on the type of sample preparation required for a viral metagenomic study that seeks to identify an aetiological virus in an animal or human diagnostic sample. A review of published virus discovery studies revealed the most commonly used enrichment methods, that were usually quick and simple to implement, namely low-speed centrifugation, filtration, nuclease-treatment (or combinations of these) which have been routinely used but often without justification. These were applied to a simple and well-characterised artificial sample composed of bacterial and human cells, as well as DNA (adenovirus) and RNA viruses (influenza A and human enterovirus), being either non-enveloped capsid or enveloped viruses. The effect of the enrichment method was assessed by both quantitative real-time PCR and metagenomic analysis that incorporated an amplification step. Reductions in the absolute quantities of bacteria and human cells were observed for each method as determined by qPCR, but the relative abundance of viral sequences in the metagenomic dataset remained largely unchanged. A 3-step method of centrifugation, filtration and nuclease-treatment showed the greatest increase in the proportion of viral sequences. This study provides a starting point for the selection of a purification method in future virus discovery studies, and highlights the need for more data to validate the effect of enrichment methods on different sample types, amplification, bioinformatics approaches and sequencing platforms. This study also highlights the potential risks that may attend selection of a virus enrichment method without any consideration for the sample type being investigated.

Metagenomic identification of a nodavirus and a circular ssDNA virus in semi-purified viral nucleic acids from the hepatopancreas of healthy Farfantepenaeus duorarum shrimp

Fisheries and aquaculture are impacted sporadically by newly emerged viral diseases. In most cases, searches for a causative pathogen only occur after a serious disease has emerged. As random shotgun sequencing (metagenomics) offers opportunities to identify novel viruses preemptively, the method was tested on nucleic acids extracted from the hepatopancreas of 12 healthy northern pink shrimp Farfantepenaeus duorarum captured from the Gulf of Mexico. Among the sequences, a nodavirus (Farfantepenaeus duorarum nodavirus, FdNV) and a virus with similarities to circoviruses and cycloviruses that possess circular single-stranded DNA (ssDNA) genomes, were identified. The FdNV genome sequence was most closely related phylogenetically to nodaviruses causing white tail disease in Macrobrachium rosenbergii and muscle necrosis disease in Litopenaeus vannamei. While the circular ssDNA virus represents the third to be detected in association with a marine invertebrate, transmission trials are needed to confirm its infectivity for F. duorarum. This study highlights the potential for using metagenomic approaches in fisheries and aquaculture industries to identify new potential pathogens in asymptomatic marine invertebrates, uncharacterized pathogens causing a new disease, or multiple pathogens associated with disease syndromes.

Novel myco-like DNA viruses discovered in the faecal matter of various animals

A wide variety of novel single-stranded DNA (ssDNA) viruses have been found in faecal matter of chimpanzees, cows, rodents, bats, badgers, foxes and pigs over the last few years. Using a combination of rolling circle amplification coupled with restriction enzyme digests based approach as well as a next generation sequencing informed approach, we have recovered fourteen full genomes of ssDNA viruses which exhibit genomic features described for members of the recently proposed gemycircularvirus group from a wide variety of mammal and bird faecal samples across New Zealand. The fourteen novel ssDNA viruses (2122-2290nt) encode two major proteins, a replication associated protein (Rep) and a capsid protein (Cp) which are bi-directionally transcribed. Interestingly, the Rep of these novel viruses are similar to gemycircularviruses detected in insects, cassava leaves, and badger faecal matter, the novel viruses share sequence similarities with the mycovirus sclerotinia sclerotiorum hypovirulence-associated DNA virus 1 (SsHADV-1) and Rep-like sequences found in fungal genomes. Pairwise sequence similarities between the 14 novel genomes with other related viral isolates (gemycircularviruses) indicated that they share greater than 55.8% genome-wide identity. Additionally, they share between 55% and 59% pairwise identity with putative novel ssDNA virus genomes recently isolated from sewage baminivirus, niminivirus and nephavirus. Based on the similarities to SsHADV-1 and Rep-like sequences found in fungal genomes, these novel gemycircularviruses may infect fungi.

Metagenomic analysis of the ferret fecal viral flora

Ferrets are widely used as a small animal model for a number of viral infections, including influenza A virus and SARS coronavirus. To further analyze the microbiological status of ferrets, their fecal viral flora was studied using a metagenomics approach. Novel viruses from the families Picorna-, Papilloma-, and Anelloviridae as well as known viruses from the families Astro-, Corona-, Parvo-, and Hepeviridae were identified in different ferret cohorts. Ferret kobu- and hepatitis E virus were mainly present in human household ferrets, whereas coronaviruses were found both in household as well as farm ferrets. Our studies illuminate the viral diversity found in ferrets and provide tools to prescreen for newly identified viruses that potentially could influence disease outcome of experimental virus infections in ferrets.

Discovery of Sclerotinia sclerotiorum Hypovirulence-Associated Virus-1 in Urban River Sediments of Heathcote and Styx Rivers in Christchurch City, New Zealand

In samples of benthic and bank river sediments of two urban rivers in Christchurch city (New Zealand), we identified and recovered isolates of Sclerotinia sclerotiorum hypovirulence-associated virus-1 (SsHADV-1), a fungus-infecting circular single-stranded DNA virus. This is the first report of SsHADV-1 outside of China and in environmental samples.

Identification of an astrovirus commonly infecting laboratory mice in the US and Japan

Mice (Mus musculus) are the most commonly used laboratory animals. Viral metagenomics on tissues of immunodeficient mice revealed sequences of a novel mammalian astrovirus. Using PCR, we screened mice from 4 breeders, 4 pharmaceutical companies, 14 research institutes and 30 universities in the US and Japan. Mice from one US breeder tested positive while none from Japanese breeders were positive for MuAstV. Mice in over half of the universities (19/30), institutes (7/14) and pharmaceutical animal facilities (2/4) investigated revealed the presence of MuAstV. Nine mice strains tested positive including both immunodeficient strains (NSG, NOD-SCID, NSG-3GS, C57BL6-Timp-3 (-/-), and uPA-NOG) and immunocompetent strains (B6J, ICR, Bash2, BALB/c). Our data indicates that MuAstV has a wide geographical, institutional and host strain distribution. Comparison of the MuAstV RdRp sequences showed numerous mutations indicating ongoing viral divergence in different facilities. This study demonstrates the need for metagenomic screening of laboratory animals to identify adventitious infections that may affect experimental outcomes.

Safety analysis of a Russian phage cocktail: from metagenomic analysis to oral application in healthy human subjects

Phage therapy has a long tradition in Eastern Europe, where preparations are comprised of complex phage cocktails whose compositions have not been described. We investigated the composition of a phage cocktail from the Russian pharmaceutical company Microgen targeting Escherichia coli/Proteus infections. Electron microscopy identified six phage types, with numerically T7-like phages dominating over T4-like phages. A metagenomic approach using taxonomical classification, reference mapping and de novo assembly identified 18 distinct phage types, including 7 genera of Podoviridae, 2 established and 2 proposed genera of Myoviridae, and 2 genera of Siphoviridae. De novo assembly yielded 7 contigs greater than 30 kb, including a 147-kb Myovirus genome and a 42-kb genome of a potentially new phage. Bioinformatic analysis did not reveal undesired genes and a small human volunteer trial did not associate adverse effects with oral phage exposure.

Evidence of the megavirome in humans

BACKGROUND

Megavirales is a proposed new virus order composed of Mimivirus, Marseillevirus and closely related viruses, as well as members of the families Poxviridae, Iridoviridae, Ascoviridae, Phycodnaviridae and Asfarviridae. The Megavirales virome, which we refer to as the megavirome, has been largely neglected until now because of the use of technical procedures that have jeopardized the discovery of giant viruses, particularly the use of filters with pore sizes in the 0.2-0.45-μm range. Concurrently, there has been accumulating evidence supporting the role of Mimivirus, discovered while investigating a pneumonia outbreak using amoebal coculture, as a causative agent in pneumonia.

OBJECTIVES

In this paper, we describe the detection of sequences related to Mimivirus and Marseillevirus in the gut microbiota from a young Senegalese man. We also searched for sequences related to Megavirales in human metagenomes publicly available in sequence databases.

RESULTS

We serendipitously detected Mimivirus- and Marseillevirus-like sequences while using a new metagenomic approach targeting bacterial DNA that subsequently led to the isolation of a new member of the family Marseilleviridae, named Senegalvirus, from human stools. This discovery demonstrates the possibility of the presence of giant viruses of amoebae in humans. In addition, we detected sequences related to Megavirales members in several human metagenomes, which adds to previous findings by several groups.

CONCLUSIONS

Overall, we present convergent evidence of the presence of mimiviruses and marseilleviruses in humans. Our findings suggest that we should re-evaluate the human megavirome and investigate the prevalence, diversity and potential pathogenicity of giant viruses in humans.

Metagenomic analysis of viral communities in (hado)pelagic sediments

In this study, we analyzed viral metagenomes (viromes) in the sedimentary habitats of three geographically and geologically distinct (hado)pelagic environments in the northwest Pacific; the Izu-Ogasawara Trench (water depth = 9,760 m) (OG), the Challenger Deep in the Mariana Trench (10,325 m) (MA), and the forearc basin off the Shimokita Peninsula (1,181 m) (SH). Virus abundance ranged from 10(6) to 10(11) viruses/cm(3) of sediments (down to 30 cm below the seafloor [cmbsf]). We recovered viral DNA assemblages (viromes) from the (hado)pelagic sediment samples and obtained a total of 37,458, 39,882, and 70,882 sequence reads by 454 GS FLX Titanium pyrosequencing from the virome libraries of the OG, MA, and SH (hado)pelagic sediments, respectively. Only 24-30% of the sequence reads from each virome library exhibited significant similarities to the sequences deposited in the public nr protein database (E-value <10(-3) in BLAST). Among the sequences identified as potential viral genes based on the BLAST search, 95-99% of the sequence reads in each library were related to genes from single-stranded DNA (ssDNA) viral families, including Microviridae, Circoviridae, and Geminiviridae. A relatively high abundance of sequences related to the genetic markers (major capsid protein [VP1] and replication protein [Rep]) of two ssDNA viral groups were also detected in these libraries, thereby revealing a high genotypic diversity of their viruses (833 genotypes for VP1 and 2,551 genotypes for Rep). A majority of the viral genes predicted from each library were classified into three ssDNA viral protein categories: Rep, VP1, and minor capsid protein. The deep-sea sedimentary viromes were distinct from the viromes obtained from the oceanic and fresh waters and marine eukaryotes, and thus, deep-sea sediments harbor novel viromes, including previously unidentified ssDNA viruses.

Detection of arboviruses and other micro-organisms in experimentally infected mosquitoes using massively parallel sequencing

Human disease incidence attributed to arbovirus infection is increasing throughout the world, with effective control interventions limited by issues of sustainability, insecticide resistance and the lack of effective vaccines. Several promising control strategies are currently under development, such as the release of mosquitoes trans-infected with virus-blocking Wolbachia bacteria. Implementation of any control program is dependent on effective virus surveillance and a thorough understanding of virus-vector interactions. Massively parallel sequencing has enormous potential for providing comprehensive genomic information that can be used to assess many aspects of arbovirus ecology, as well as to evaluate novel control strategies. To demonstrate proof-of-principle, we analyzed Aedes aegypti or Aedes albopictus experimentally infected with dengue, yellow fever or chikungunya viruses. Random amplification was used to prepare sufficient template for sequencing on the Personal Genome Machine. Viral sequences were present in all infected mosquitoes. In addition, in most cases, we were also able to identify the mosquito species and mosquito micro-organisms, including the bacterial endosymbiont Wolbachia. Importantly, naturally occurring Wolbachia strains could be differentiated from strains that had been trans-infected into the mosquito. The method allowed us to assemble near full-length viral genomes and detect other micro-organisms without prior sequence knowledge, in a single reaction. This is a step toward the application of massively parallel sequencing as an arbovirus surveillance tool. It has the potential to provide insight into virus transmission dynamics, and has applicability to the post-release monitoring of Wolbachia in mosquito populations.