Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSphere
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Applied and Environmental Science

Metagenomic Analysis of the Virome of Mosquito Excreta

Ana L. Ramírez, Agathe M. G. Colmant, David Warrilow, Bixing Huang, Alyssa T. Pyke, Jamie L. McMahon, Dagmar B. Meyer, Rikki M. A. Graham, Amy V. Jennison, Scott A. Ritchie, Andrew F. van den Hurk
Garret Suen, Editor
Ana L. Ramírez
aCollege of Public Health, Medical and Veterinary Sciences, James Cook University, Smithfield, Queensland, Australia
bAustralian Institute of Tropical Health and Medicine, Smithfield, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ana L. Ramírez
Agathe M. G. Colmant
cSchool of Chemistry and Molecular Biosciences, The University of Queensland, St. Lucia, Queensland, Australia
dAustralian Infectious Diseases Research Centre, The University of Queensland, St. Lucia, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David Warrilow
ePublic Health Virology, Forensic and Scientific Services, Department of Health, Coopers Plains, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bixing Huang
ePublic Health Virology, Forensic and Scientific Services, Department of Health, Coopers Plains, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alyssa T. Pyke
ePublic Health Virology, Forensic and Scientific Services, Department of Health, Coopers Plains, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Alyssa T. Pyke
Jamie L. McMahon
ePublic Health Virology, Forensic and Scientific Services, Department of Health, Coopers Plains, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dagmar B. Meyer
aCollege of Public Health, Medical and Veterinary Sciences, James Cook University, Smithfield, Queensland, Australia
bAustralian Institute of Tropical Health and Medicine, Smithfield, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rikki M. A. Graham
fPublic Health Microbiology, Forensic and Scientific Services, Department of Health, Coopers Plains, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Amy V. Jennison
fPublic Health Microbiology, Forensic and Scientific Services, Department of Health, Coopers Plains, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Amy V. Jennison
Scott A. Ritchie
aCollege of Public Health, Medical and Veterinary Sciences, James Cook University, Smithfield, Queensland, Australia
bAustralian Institute of Tropical Health and Medicine, Smithfield, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrew F. van den Hurk
ePublic Health Virology, Forensic and Scientific Services, Department of Health, Coopers Plains, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Garret Suen
University of Wisconsin–Madison
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSphere.00587-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Traditional screening for arboviruses in mosquitoes requires a priori knowledge and the utilization of appropriate assays for their detection. Mosquitoes can also provide other valuable information, including unexpected or novel arboviruses, nonarboviral pathogens ingested from hosts they feed on, and their own genetic material. Metagenomic analysis using next-generation sequencing (NGS) is a rapidly advancing technology that allows us to potentially obtain all this information from a mosquito sample without any prior knowledge of virus, host, or vector. Moreover, it has been recently demonstrated that pathogens, including arboviruses and parasites, can be detected in mosquito excreta by molecular methods. In this study, we investigated whether RNA viruses could be detected in mosquito excreta by NGS. Excreta samples were collected from Aedes vigilax and Culex annulirostris experimentally exposed to either Ross River or West Nile viruses and from field mosquitoes collected across Queensland, Australia. Total RNA was extracted from the excreta samples, reverse transcribed to cDNA, and sequenced using the Illumina NextSeq 500 platform. Bioinformatic analyses from the generated reads demonstrate that mosquito excreta provide sufficient RNA for NGS, allowing the assembly of near-full-length viral genomes. We detected Australian Anopheles totivirus, Wuhan insect virus 33, and Hubei odonate virus 5 and identified seven potentially novel viruses closely related to members of the order Picornavirales (2/7) and to previously described, but unclassified, RNA viruses (5/7). Our results suggest that metagenomic analysis of mosquito excreta has great potential for virus discovery and for unbiased arbovirus surveillance in the near future.

IMPORTANCE When a mosquito feeds on a host, it ingests not only its blood meal but also an assortment of microorganisms that are present in the blood, thus acting as an environmental sampler. By using specific tests, it is possible to detect arthropod-borne viruses (arboviruses) like dengue and West Nile viruses in mosquito excreta. Here, we explored the use of next-generation sequencing (NGS) for unbiased detection of RNA viruses present in excreta from experimentally infected and field-collected mosquitoes. We have demonstrated that mosquito excreta provide a suitable template for NGS and that it is possible to recover and assemble near-full-length genomes of both arboviruses and insect-borne viruses, including potentially novel ones. These results importantly show the direct practicality of the use of mosquito excreta for NGS, which in the future could be used for virus discovery, environmental virome sampling, and arbovirus surveillance.

INTRODUCTION

In effect, female mosquitoes act as environmental samplers (“biological syringes”) that feed on the blood of a variety of vertebrate hosts (1). Mosquitoes harbor arthropod-borne viruses (arboviruses), which are capable of replicating within and being transmitted by mosquito vectors to vertebrate hosts, and insect-specific viruses (ISV) (2), as well as nonarboviruses that do not replicate in mosquitoes but might be present in hosts they feed upon (3). Traditionally, molecular assays widely utilized in arbovirus surveillance programs screen only for known and characterized endemic and enzootic viruses, including broadly reactive pan-flavivirus and pan-alphavirus assays. It is likely that many other viruses, regardless of pathogenicity, may remain undetected. Metagenomic analysis using next-generation sequencing (NGS) has proven to be useful even when a finite amount of sample is available and has allowed the unbiased identification of viruses, mosquito species, and endosymbionts, such as Wolbachia, from a single mosquito in a single reaction (4). Viral metagenomics is also versatile and has been successfully used in Australia for the identification of multiple arboviruses, including novel rhabdoviruses, bunyaviruses (5), and mesoniviruses (6) from field-collected mosquitoes.

Vector-enabled metagenomics could also be used as a tool to monitor human and animal diseases (1), an application often referred to as xenosurveillance (3). Xenosurveillance offers an alternative to directly sampling hosts, which involves a time-consuming process, that for humans, requires individual informed consent or, for animals, necessitates animal ethics approval prior to commencement. Alternatively, by proxy, xenosurveillance using mosquitoes has been successfully used to detect circulating H5N1 influenza virus (7), Epstein-Barr virus, canine distemper virus (3), human herpesvirus, human papillomaviruses, anelloviruses and circoviruses, among others (8). This methodology has also been used to study other pathogens, such as filarial parasites (9) or apicomplexans (10). In Sri Lanka, xenosurveillance has been successfully used to map areas with persistent circulation of Wuchereria bancrofti after mass drug administration programs (11).

Although metagenomic approaches are becoming more accessible, the costs associated with processing many samples can be prohibitive, especially in low-resource settings (12). However, analyzing mosquito excreta has the advantage of reducing the number of samples that need to be processed (1 or 2 samples versus several pools of mosquitoes in each trap). Metagenomic analysis of mosquito excreta could potentially allow unbiased identification of circulating viruses (pathogenic and nonpathogenic) from a given locale in a single reaction from a single sample.

Recently, arboviruses such as the dengue viruses (DENVs), Ross River virus (RRV), and West Nile virus (WNV) have been detected in mosquito excreta by reverse-transcription real-time PCR (RT-rtPCR) by us and others (13, 14). Additionally, hepatitis B virus, which does not replicate in the mosquito, has been detected in mosquito excreta by RT-PCR and Southern blotting up to 72 h after the ingestion of an infectious blood meal (15). A system to collect mosquito excreta in the field for the detection of circulating arboviruses has been developed, with mosquito excreta samples having been found to be positive for WNV, RRV, and Murray Valley encephalitis virus (MVEV) in Australia (16).

In this study, we evaluated the application of NGS for unbiased detection of RNA viruses, first in excreta harvested from experimentally infected mosquitoes and, later, in excreta of field-collected mosquitoes sampled from different locations in the state of Queensland, Australia. The sampling locations in Queensland have a history of transmission of medically important arboviruses, such as RRV, WNV, and Barmah Forest virus (17, 18), as well as evidence of a diverse range of insect-specific viruses (19, 20). This is the first reported study to investigate mosquito excreta as a viable and practical sample type for NGS-based metagenomics.

RESULTS

Laboratory studies.Excreta samples were collected from groups of 5 Culex annulirostris or Aedes vigilax mosquitoes exposed to either WNV (Kunjin subtype [WNVKUN]) or RRV, by swabbing a Parafilm disc placed in the bottom of each container (14). RNA extracted from all six excreta samples collected from groups of experimentally infected mosquitoes were positive for either WNVKUN or RRV by RT-rtPCR, with mean (± standard error of the mean [SEM]) cycle threshold (CT) values of 26.3 ± 0.4 and 26.8 ± 1.0, respectively. Based on these results, six libraries were subsequently prepared and next-generation sequenced using the Illumina NextSeq platform (21). The mean (± SEM) numbers of raw reads obtained from the libraries from mosquitoes exposed to WNVKUN or RRV were 5, 877,227 ± 137,555 and 9,020,173 ± 515,578, respectively. Preliminary DIAMOND/MEGAN analysis accurately classified reads obtained from experimentally infected mosquitoes as either WNVKUN or RRV. Results from subsequent assembly demonstrated that excreta collected from small groups (≤5) of mosquitoes experimentally infected with arboviruses provide sufficient template for NGS, allowing the assembly of near-full-length viral genomes (Fig. 1).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Sequence depth and coverage of virus genomes. Each sequence was obtained from excreta from a group of five experimentally infected mosquitoes exposed to either WNVKUN (A to C) or RRV (D to F). Sequences were assembled to a reference sequence (GenBank accession numbers KX394395.1 and GQ433359.1 for WNVKUN and RRV, respectively). Sequence depth (y axis) and coverage (x axis) are shown.

Field studies.To determine if viral genomes could be sequenced from excreta sampled from mosquitoes collected using the methods of Meyer et al. (16), studies were conducted in Far North Queensland (FNQ) and South East Queensland (SEQ) in March and April 2018 (Fig. S1). Excreta were sampled from 19,751 mosquitoes representing at least 25 species collected in 24 traps deployed in FNQ and from 2,531 mosquitoes representing at least 18 species collected in 22 SEQ traps (Data Set S1). Viral RNA was extracted, and RRV and BFV were detected using specific RT-rtPCR assays (22, 23). RRV RNA was detected in 2 of 46 (4%) excreta samples, which had been collected from White Rock and Cattana wetlands (FNQ), with CT values of 36.9 and 36.0, respectively (Table S1). BFV RNA was detected in a single sample from Cattana wetlands (2%; CT = 36.3), which coincidentally was also positive for RRV. A total of 703 female mosquitoes from the two traps that contained the two RRV- and/or BFV-positive excreta samples mentioned above were sorted into 42 pools of ≤50 mosquitoes and screened for RRV and BFV using RT-rtPCR. Five (12%) pools were positive for RRV, and five (12%) pools were positive for BFV (Table S2).

TABLE S1

Description of excreta samples from field-collected mosquitoes. Download Table S1, DOCX file, 0.02 MB.
Copyright © 2020 Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S2

Mosquitoes from White Rock (FNQ8) and Cattana Wetlands (FNQ21), far North Queensland, Australia, processed for virus detection from traps that yielded excreta samples positive for Ross River and Barmah Forest viruses. Download Table S2, DOCX file, 0.02 MB.
Copyright © 2020 Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S1

Study sites in Queensland, Australia, where mosquitoes and mosquito excreta were collected. (A) Far North Queensland (FNQ); (B) South East Queensland (SEQ). Download FIG S1, EPS file, 1.0 MB.
Copyright © 2020 Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

DATA SET S1

Mosquito species identified in this study. Download Data Set S1, XLSX file, 0.02 MB.
Copyright © 2020 Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

A total of 47 libraries (including a fetal calf serum [FCS] negative control), corresponding to excreta samples collected from 11 locations, were sequenced using the Illumina NextSeq platform, with a mean (± SEM) of 12,776,515 ± 565,467 raw reads per sample (Table S1). Upon preliminary DIAMOND/MEGAN analyses using a threshold of 1,000 assigned reads, no sequences from field-collected samples were taxonomically classified as known pathogenic arboviruses, including those that were positive by RT-rtPCR for RRV and BFV. However, other RNA virus genomes were detected by NGS in 22 of the 46 excreta samples sequenced, with some samples containing up to 3 different viruses (Fig. 2).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Virus genomes detected by NGS of excreta from field-collected mosquitoes. The virus names shown are from the taxonomic annotation by DIAMOND and MEGAN. These sequences were used for initial assembly and to calculate the percentage of nucleotide identity of the final consensus sequence. An asterisk indicates that only partial sequences were obtained. AATV, Australian Anopheles totivirus (NC_035674.1); HBOV5, Hubei odonate virus 5 (NC_033206.1); HiPV, Himetobi P virus (AB017037.1); HBPLV41, Hubei picorna-like virus 41 (NC_033238.1); HBTMV2, Hubei tetragnatha maxillosa virus 2 (KX883694.1); WIV33, Wuhan insect virus 33 (NC_033722.1); ArIFV, Armigeres iflavirus (LC310707.1); Daeseongdong virus 2 (NC_028489.1); DCV, Drosophila C virus (NC_001834.1); HBAV1, Hubei arthropod virus 1 (KX883297.1); HBOV7, Hubei odonate virus 7 (NC_033232.1); HBPLV61, Hubei picorna-like virus 61 (NC_033003.1); ZJMV1, Zhejiang mosquito virus 1 (NC_033716.1). Only libraries with identified virus genomes are listed.

The 13 viruses identified in this study were found to be related to members of the order Picornavirales and previously described unclassified RNA viruses. For three of the viruses detected (Himetobi P virus [HiPV], Hubei tetragnatha maxillosa virus 2 [HBTMV2], and Hubei picorna-like virus 41 [HBPLV41]) only partial sequences were obtained; for the rest of the viruses, we were able to assemble near-full-length genomes. Of these, Australian Anopheles totivirus (AATV), Wuhan insect virus 33 (WIV33), and Hubei odonate virus 5 (HBOV5) showed >90% amino acid (aa) identity over the RNA-dependent RNA polymerase (RdRp) to published sequences in GenBank, indicating that they correspond to strains of these viruses (19, 24). Among them, the four HBOV5 sequences were > 98% identical to each other.

We identified seven potentially novel virus species (Table 1) and constructed phylogenetic trees based on their RdRp sequences. Warroolaba Creek virus 2 (WRCV2) is related to Drosophila C virus (DCV), a cripavirus belonging to the family Dicistroviridae (Fig. 3). The sequence shared ∼84% aa identity with DCV. Sequences from three samples collected in SEQ contained Perrin Park virus (PPKV), which is phylogenetically similar to Armigeres iflavirus (ArIFV; 71% aa identity), an iflavirus from the family Iflaviridae, which was first isolated from Armigeres sp. mosquitoes in the Philippines (25). The three PPKV sequences generated here were >99% identical to each other.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Novel RNA virus sequences identified in this study

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Phylogenetic relationships of viruses related to the order Picornavirales and other unclassified RNA viruses discovered in field mosquito excreta. A multiple-sequence alignment of the RNA-dependent RNA polymerase amino acid sequences was used to create a maximum-likelihood phylogeny using 100 bootstrap replicates; an asterisk indicates node support of >70% bootstrap support. The tree was midpoint rooted. The potential novel viruses discovered in this study are shown in bold and color coded: green triangles for samples collected in South East Queensland and blue circles for samples collected in Far North Queensland. Unclassified RNA viruses described by Shi et al. (24) are shown with gray squares. Corresponding GenBank accession numbers for compared virus sequences are provided in parentheses.

We also obtained sequences closely related to unclassified RNA viruses identified as a part of a large-scale invertebrate virosphere survey conducted on samples from China by Shi and colleagues (24) (Fig. 3). A potentially novel virus (Warroolaba Creek virus 1 [WRCV1]) from a sample from SEQ was related to Hubei odonate virus 7 (HBOV7), sharing 36% aa identity in the RdRp. Five samples from FNQ contained Redbank virus (REBV) sequences related to Hubei arthropod virus 1 (HBAV1; 59% aa identity), with >99% similarity between them. Four samples from FNQ and a single sample from SEQ contained Old Port virus (OPTV) sequences related to Zhejiang mosquito virus 1 (ZJMV1; 85% aa identity) with >97% similarity with each other. Sequences from one sample collected (Warroolaba Creek virus 3 [WRCV3]) in SEQ were phylogenetically similar to Hubei picorna-like virus 61 (HBPLV61; 83% aa similarity), which had been previously identified in mosquitoes.

Finally, three samples from FNQ contained Smithfield permutotetra-like virus (SmPLV) sequences related to both Culex Daeseongdong-like virus and Daeseongdong virus 2 (DV2) (Fig. 4). The sequences shared ∼84% aa identity with both unclassified RNA viruses, which are themselves highly similar to each other (>99%) and were identified in Culex mosquitoes from Korea and California, respectively (26, 27). Our three SmPLV sequences were >99% identical to each other.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Phylogenetic relationships of unclassified RNA viruses discovered in field mosquito excreta. A multiple-sequence alignment of the RNA-dependent RNA polymerase amino acid sequences was used to create a maximum-likelihood phylogeny using 100 bootstrap replicates; an asterisk indicates node support of >70% bootstrap support. The tree was midpoint rooted. The potential novel viruses discovered in this study are shown in bold and color coded with blue circles for samples collected in Far North Queensland. Unclassified RNA viruses described by Shi et al. (24) are shown with gray squares. Corresponding GenBank accession numbers for compared virus sequences are provided in parentheses.

The FNQ collections appeared to have greater virus diversity, with full-length sequences from 6 different viruses identified only from FNQ (AATV, HBO5V, HiPV, WIV33, SmPLV, and REBV), compared with 4 identified only from SEQ (PPKV, WRCV1, WRCV2, and WRCV3). Only OPTV was detected in both FNQ and SEQ samples.

Evaluation of rRNA in excreta samples.Since an rRNA depletion kit was not used in our study, we evaluated the extent of rRNA presence in excreta samples where a virus was found by NGS in both laboratory-infected and field-collected mosquitoes. We did a BLAST search of a subset of reads (n = 1,000) against SILVA small-subunit (SSU) and large-subunit (LSU) rRNA sequences databases (28). The amount of SSU RNA and LSU rRNA varied from 12% to 75% of reads and from 19% to 93% of reads, respectively, suggesting a large amount of rRNA present in the samples (Table S3).

TABLE S3

rRNA sequences in mosquito excreta samples determined by BLAST against SILVA SSU and LSU databases. Download Table S3, DOCX file, 0.02 MB.
Copyright © 2020 Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

DISCUSSION

Over the past decade, unbiased metagenomic analysis using NGS has become a valuable tool and revolutionized virus discovery and surveillance (29, 30). However, as with all detection-based technologies, these techniques are prone to limitations associated with sample type, integrity, availability, and amount. The current study aimed to circumvent some of these issues pertinent in surveillance strategies by coupling the less biased sampling afforded by mosquito feeding behaviors, which is further simplified by collection of excreta, with unbiased virome detection made possible via NGS. Our laboratory findings demonstrate that excreta from experimentally infected mosquitoes provide a sufficient template for sequencing and assembly of near-full-length virus genomes. However, we were unsuccessful at sequencing RRV or BFV from the two field-collected samples that were positive by RT-rtPCR. A possible explanation for this is the likelihood that NGS is not as sensitive as RT-rtPCR for detection of viruses with low titers (31) and those present in complex samples (32). In our study, samples from experimentally infected mosquitoes contained a larger amount of starting template (as evidenced by lower CT values) and were likely subjected to less sample degradation than field samples. Since the viral titers to which mosquitoes were exposed to during laboratory infection may have been higher than those encountered by mosquitoes in the wild, future investigations should include determination of the limit of detection for this method by exposing mosquitoes to different viral titers.

Despite this, we could detect arbovirus sequences in samples with CT values of ≤28.5. Among the libraries that were positive by RT-rtPCR for RRV and/or BFV but negative by NGS, all had CT values of ≥36.0, suggesting that the failure to detect these viruses may have resulted from low abundance of viral sequences. Further, our results indicate a large amount of rRNA present in mosquito excreta, with an average of more than 50% of sequence reads matching both SSU and LSU rRNA, which could potentially impact the sensitivity of virus detection. In the future, it is likely that by adding an rRNA depletion step, sensitivity can be increased. Additionally, viral enrichment protocols such as NetoVIR (33), which have been used for single-mosquito viral metagenomics (34), could also be used to increase the sensitivity for virome analysis. However, it is important to note that we were able to assemble near-full-length genomes from excreta samples where up to 95% of sequence reads corresponded to rRNA. With increased application of NGS, it is likely that new, improved protocols will be devised, increasing the efficiency and sensitivity of this sequencing platform. Coinciding with NGS technological advancements, mosquito excreta could be utilized as a valuable alternative sample for routine arbovirus surveillance, enabling the unbiased detection of arboviruses of public health importance from a single sample, thereby increasing the surveillance scope. Furthermore, this approach together with NGS metagenomic analysis of mosquito saliva deposited on honey-soaked FTA cards (35) could potentially provide a wider and more comprehensive overview of pathogenic and nonpathogenic microbiota circulating in given locales.

By performing sequencing of field-collected mosquito excreta samples, we were able to show evidence of the circulation of 13 insect-borne viruses, of which one (AATV) had been previously identified in mosquitoes and four (PPKV, SmPLV, WRCV3, and OPTV) were related to viruses previously identified in mosquitoes. Some of these viruses have been isolated and/or characterized. In the case of AATV, although the virus has not been isolated in C6/36 mosquito cells, evidence of its replication in mosquitoes has been observed by small-RNA analysis (19). ArIFV, to which PPKV is related, causes pronounced cytopathic effect in C6/C6 cells; however, its pathogenicity in mosquitoes still needs to be evaluated (25). Some iflaviruses, such as slow bee paralysis virus and infectious flacherie virus, can cause lethal infections in their host (36, 37). Viruses related to SmPLV, WRCV3, and OPTV that have been previously identified in mosquitoes remain unclassified and have been characterized only in silico (24, 26, 27). Clearly, more work needs to be done on the novel viruses identified in this study with respect to their phenotypic characterization.

Although these insect-borne viruses are unlikely to be associated with disease in vertebrates, they can potentially affect the vector competence of mosquitoes for pathogenic viruses, as demonstrated previously for a number of insect-specific viruses (38, 39). Because of this, elucidating the mosquito virome is critical for understanding the role a mosquito species plays in arbovirus transmission cycles and potential control strategies. Further, novel insect-specific viruses have been shown to be potentially useful as biological control agents or as platforms for vaccine and diagnostic development (2), and therefore, expanded detection of these agents could have important downstream applications.

As of 2020, the costs associated with NGS and the time and bioinformatics skills required to analyze the results from thousands to tens of thousands of mosquitoes can be prohibitive, especially in low-resource settings (12). We have shown that mosquito excreta can be used as a preliminary sample for virus discovery in field populations of mosquitoes. Using mosquito excreta has the advantage of reducing costs by requiring sequencing of only one sample from a trap, instead of multiple pools of mosquitoes, which in traps containing as many as 10,000 mosquitoes could result in 200 to 400 pools depending on pool size. Based on the results obtained from excreta, the mosquitoes could be used for subsequent sequencing or in attempts to detect virus for subsequent isolation. Indeed, we detected RRV and BFV by RT-rtPCR in mosquito pools prepared from traps that yielded RT-rtPCR-positive excreta. Although these viruses were detected using RT-rtPCR, this result highlights the utility of the excreta sampling approach in identifying viral RNA in the first instance, which, in the future, could more suitably be applied to NGS, as the sensitivity of NGS increases.

A limitation of our field study is that it is almost impossible to attribute the excreta deposited on the polycarbonate substrate to a particular insect. Although mosquitoes comprise the majority of the collections, traps used to capture mosquitoes also attract nontarget insects (40), which could feed on the honey and excrete on the substrate. This can be reflected by the detection of HBOV5, which is associated with dragonflies and damselflies, and HBTMV2, which is associated with spiders (24). To overcome this technical limitation, the mosquitoes could first be sorted and transferred to clean containers in the laboratory from which excreta would be obtained and sequenced to confirm the origin of each virus. Additionally, based on results obtained from sequencing mosquito excreta, virus detection and subsequent isolation could be attempted from the mosquito homogenates for viruses that can be grown in vitro. Use of traps that are highly selective for mosquitoes, such as CO2-baited passive box traps that do not deploy lights that attract other insects, is recommended (16, 41).

It is evident that NGS technologies have many applications for the study of vectors and the pathogens they transmit (42). Here, we have demonstrated that metagenomic analysis of mosquito excreta can be used in the context of arbovirus surveillance for virus discovery, as well as for unbiased virome sampling, particularly as costs decrease and technologies become more accessible.

MATERIALS AND METHODS

Excreta collection from experimentally infected mosquitoes.Field-collected C. annulirostris and laboratory-reared A. vigilax mosquitoes were exposed to defibrinated sheep blood (Institute of Medical and Veterinary Science, Adelaide, Australia) containing either WNV (Kunjin subtype [WNVKUN]; 107.3±0.3 50% tissue culture infective doses [TCID50]/ml) or RRV (108.1±0.1 TCID50/ml) via the hanging drop method (43) and by using a Hemotek membrane feeder (Discovery Workshops, Accrington, Lancashire, UK), respectively. For each virus, three groups of five fully engorged females were placed in 200-ml polypropylene containers modified with an insect screen floor for excreta collection (14). Four to 7 days after mosquitoes fed on the infectious blood meal, a single excreta sample was collected from each container using a cotton swab moistened with growth medium (GM) (Opti-MEM [Gibco, Invitrogen Corporation, Grand Island, NY] supplemented with 3% fetal bovine serum [In Vitro Technologies, Nobel Park North, VIC, Australia], antibiotics and antimycotics) as described by Ramírez et al. (14). The swab was then placed in a 2-ml free-standing tube containing 1 ml GM and stored at −80°C. All mosquito exposures were conducted under biosafety level 2 (BSL-2) conditions, and mosquitoes were maintained at 28°C, 75% relative humidity (RH), and a 12:12 (light:dark) photoperiod within an environmental growth cabinet.

Excreta collection from field mosquitoes.Mosquito excreta samples were collected in March and April 2018 in SEQ and FNQ, Australia, from different sites that encompassed a variety of vertebrate host species, including flying foxes, wading birds, livestock, and macropods (Fig. S1). Adult mosquitoes were collected using Centers for Disease Control and Prevention (CDC) model 512 light traps (John W. Hock Company, Gainesville, FL) fitted with collection containers modified from the design described in reference 16. Collection containers housed two filter paper cards soaked with blue-dyed honey and a removable polycarbonate substrate for excreta collection. All traps were baited with 1 kg of dry ice as a source of CO2; traps deployed in FNQ were also supplemented with 1-octen-3-ol to increase capture of mosquitoes (44). Traps were operated for 14 h overnight, before being transported to the laboratory and placed in a humidified box, where mosquitoes were allowed to feed on the honey-soaked filter paper cards for an additional 24 h to increase the amount of excreta produced. Excreta were collected from the polycarbonate insert using a moistened swab as described above and stored at −80°C. The mosquitoes from each trap were cold-anesthetized and placed in 70-ml vials before being stored at −80°C. To avoid cross-contamination, the traps and polycarbonate inserts were handled with gloves during trap assembly and retrieval, and gloves were changed after each sample was collected. After collection, the polycarbonate inserts were soaked in 1% bleach and then rinsed and wiped with 70% ethanol, while the pots and mesh inserts were wiped with 1% bleach followed by 70% ethanol.

Virus assays.Thawed excreta samples were shaken to release the material from the swab using a TissueLyser II (Qiagen, Hilden, Germany) for 3 min at 26 Hz and centrifuged at 14,000 × g for 30 s to remove droplets from the lid (14). RNA from excreta from experimentally infected mosquitoes was extracted from the supernatant using the QIAamp viral RNA minikit (Qiagen, Hilden, Germany); RNA from excreta from field-collected mosquitoes was extracted using the One-For-All nucleic acid kit (Qiagen, Hilden, Germany). Extractions were carried out with and without carrier (to avoid interference in downstream NGS), according to the manufacturer’s instructions. RNA extracted with carrier was immediately tested by RT-rtPCR, and RNA extracted without carrier was stored at −80°C for sequencing.

Excreta samples from experimentally infected mosquitoes were screened for either WNVKUN (23) or RRV (22) by RT-rtPCR. Field-collected excreta samples were tested for RRV and BFV (23), as these two alphaviruses are endemic in the sampling locations. In addition, mosquitoes from traps where arboviral RNA was detected in excreta by RT-rtPCR were identified morphologically by an experienced medical entomologist (A.F.V.D.H) and placed into pools of ≤50 mosquitoes by species before being homogenized in 2 ml of GM as described above; viral RNA was extracted using the QIAamp viral RNA minikit with carrier following the manufacturer’s recommendations. Viral RNA was detected in excreta samples and mosquito pools using RT-rtPCR assays specific for RRV, WNVKUN, and BFV on a Rotor-Gene 600 real-time thermocycler (Qiagen, Australia). Every run included synthetic primer and probe controls, a positive extraction control (bovine viral diarrheal virus [BVDV]), a negative extraction control, and a no-template control (molecular-grade water). The results were evaluated qualitatively: for any sample, a threshold cycle number (CT) of >40 indicated that no RNA was detected (45).

Sequencing.Genomic DNA (gDNA) was removed from total RNA aliquots extracted without carrier RNA using the Heat & Run gDNA removal kit (ArcticZymes, Tromso, Norway). RNA was reverse transcribed, and first-strand cDNA was synthesized using the NEB Protoscript II first-strand cDNA synthesis kit (New England BioLabs, Ipswich MA) according to the manufacturer’s instructions, followed by second-strand cDNA synthesis using NEB second-strand synthesis enzyme buffer (New England BioLabs, Ipswich, MA) and second-strand DNA enzymes (DNA polymerase I [Escherichia coli], 10 U; RNase H, 0.35 U; and E. coli DNA ligase, 1.25 U; New England BioLabs, Ipswich, MA). The newly synthesized DNA was purified by ethanol precipitation. DNA libraries were prepared using the Nextera XT DNA library preparation kit (Illumina) and Illumina Nextera XT index kit. The resulting libraries were analyzed, and DNA sizing and quantification were performed using a 2200 TapeStation (Agilent Technologies). A fetal calf serum (FCS) library was prepared as described above from FCS RNA as a negative control. Libraries were diluted to 1 nM, pooled, denatured and diluted to a final concentration of 1.2 pM. Paired-end sequencing was performed using the NextSeq platform (Illumina) using a NextSeq 500 Mid Output V2 kit (Illumina) (21).

Sequence analysis and phylogenetics.Sequence reads were demultiplexed and adapters were removed using bcl2fastq version 2.20 (http://sapac.support.illumina.com/downloads/bcl2fastq-conversion-software-v2-20.html). An initial search of the raw sequences was conducted using DIAMOND (BLASTx) (46) against the NCBI-nr viral protein reference sequence database (downloaded on 20 March 2019). Taxonomic binning of the reads was performed using MEGAN CE version 6.15.2 (47) using the naive LCA algorithm (minimum score = 75.0; maximum expected = 0.1; top percent = 10.0; minimum support = 10). Based on these results, RNA virus sequences with at least 1,000 reads assigned as close relatives were selected as references for assembly. Viral sequences were assembled using Geneious Prime version 2019.0.4 by either de novo assembly or manually mapping the reads to a reference sequence obtained from GenBank. For de novo assembly, low-quality (Q < 30) and short (<100-nt) reads were trimmed using BBDuk before being assembled using SPAdes (48) with default parameters. Assembled sets of overlapping DNA sequences (contigs) were then mapped to a reference sequence, or the longest contigs containing an open reading frame were compared against the NCBI-nr database using BLAST (49) and used as a reference sequence for further assembly. Alternatively, raw reads were mapped against a reference sequence using default settings, and the consensus sequence of this assembly was used for subsequent assembly (19). This process was repeated until the length of the consensus sequence did not increase anymore. Reads from experimentally infected mosquitoes or from field-collected mosquitoes with a positive arbovirus RT-rtPCR result were mapped against a reference sequence.

Assembled sequences with > 90% amino acid identity with existing viruses over the RNA-dependent RNA polymerase (RdRp) were considered strains of these viruses and assumed to phylogenetically group with them. Because of this, only complete sequences with less than 90% amino acid identity to the reference sequence were included in phylogenetic analyses. The translated contigs were aligned with protein sequences obtained from GenBank using the results from BLAST and previously published phylogenetic trees of the related viruses (24, 25). Multiple protein alignments were done using MAFFT v7.388 (50) and trimmed using TrimAl (51). The optimal evolutionary model was selected using the Akaike information criterion in SMS (52). Maximum-likelihood phylogenies were generated using the Le Gascuel (LG) amino acid substitution model with 100 bootstrap replicates using PhyML v 3.3 (53).

Evaluation of rRNA contamination in excreta samples.Since an rRNA depletion kit was not used in our study, we evaluated the extent of rRNA contamination in our excreta samples. For this, we estimated the amount of rRNA on a subset of reads (n = 1,000) from all the excreta samples shown to contain a virus by NGS (from both laboratory-infected and field-collected mosquitoes) by performing a BLAST search (49) against SILVA small subunit (SSU) and large subunit (LSU) rRNA sequences databases (https://www.arb-silva.de/) (28) using Geneious Prime version 2020.2.2. Finally, for each sample, the percentage of sequences that got a hit was calculated for each database.

Data availability.Raw sequence reads generated from this study have been deposited in the NCBI Sequence Read Archive (SRA) repository under project number PRJNA631724. All virus genome sequences generated in this study have been deposited in GenBank under accession numbers MN784056 to MN784081.

ACKNOWLEDGMENTS

This work was supported by a National Health and Medical Research Council-funded “Improving Health Outcomes in the Tropical North: A multidisciplinary collaboration (HOT NORTH)” grant (number 1131932). A.L.R. was supported by a doctoral scholarship granted by the Australian Institute of Tropical Health and Medicine, James Cook University QLD, Australia. S.A.R. was supported, in part, by National Health and Medical Research Council research fellowship 1044698. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

We thank Noel Giebel and Barbara Williams for letting us collect mosquitoes in their property. We also thank Michael Townsend and Kyran Staunton for their assistance in the field and Jean Barcelon, Carol Kestler, Glen Hewitson, Ian Mackay, and Peter Moore for their assistance and advice with various parts of the study.

FOOTNOTES

    • Received June 29, 2020.
    • Accepted August 27, 2020.
  • Copyright © 2020 Ramírez et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Brinkmann A,
    2. Nitsche A,
    3. Kohl C
    . 2016. Viral metagenomics on blood-feeding arthropods as a tool for human disease surveillance. Int J Mol Sci 17:1743. doi:10.3390/ijms17101743.
    OpenUrlCrossRef
  2. 2.↵
    1. Bolling BG,
    2. Weaver SC,
    3. Tesh RB,
    4. Vasilakis N
    . 2015. Insect-specific virus discovery: significance for the arbovirus community. Viruses 7:4911–4928. doi:10.3390/v7092851.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Grubaugh ND,
    2. Sharma S,
    3. Krajacich BJ,
    4. Fakoli LSI,
    5. Bolay FK,
    6. Diclaro JWI,
    7. Johnson WE,
    8. Ebel GD,
    9. Foy BD,
    10. Brackney DE
    . 2015. Xenosurveillance: a novel mosquito-based approach for examining the human-pathogen landscape. PLoS Negl Trop Dis 9:e0003628. doi:10.1371/journal.pntd.0003628.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Hall-Mendelin S,
    2. Allcock R,
    3. Kresoje N,
    4. van den Hurk AF,
    5. Warrilow D
    . 2013. Detection of arboviruses and other micro-organisms in experimentally infected mosquitoes using massively parallel sequencing. PLoS One 8:e58026. doi:10.1371/journal.pone.0058026.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Coffey LL,
    2. Page BL,
    3. Greninger AL,
    4. Herring BL,
    5. Russell RC,
    6. Doggett SL,
    7. Haniotis J,
    8. Wang C,
    9. Deng X,
    10. Delwart EL
    . 2014. Enhanced arbovirus surveillance with deep sequencing: identification of novel rhabdoviruses and bunyaviruses in Australian mosquitoes. Virology 448:146–158. doi:10.1016/j.virol.2013.09.026.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Warrilow D,
    2. Watterson D,
    3. Hall RA,
    4. Davis SS,
    5. Weir R,
    6. Kurucz N,
    7. Whelan P,
    8. Allcock R,
    9. Hall-Mendelin S,
    10. O'Brien CA,
    11. Hobson-Peters J
    . 2014. A new species of mesonivirus from the Northern Territory, Australia. PLoS One 9:e91103. doi:10.1371/journal.pone.0091103.
    OpenUrlCrossRef
  7. 7.↵
    1. Barbazan P,
    2. Thitithanyanont A,
    3. Misse D,
    4. Dubot A,
    5. Bosc P,
    6. Luangsri N,
    7. Gonzalez J-P,
    8. Kittayapong P
    . 2008. Detection of H5N1 avian influenza virus from mosquitoes collected in an infected poultry farm in Thailand. Vector Borne Zoonotic Dis 8:105–109. doi:10.1089/vbz.2007.0142.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    1. Ng TFF,
    2. Willner DL,
    3. Lim YW,
    4. Schmieder R,
    5. Chau B,
    6. Nilsson C,
    7. Anthony S,
    8. Ruan Y,
    9. Rohwer F,
    10. Breitbart M
    . 2011. Broad surveys of DNA viral diversity obtained through viral metagenomics of mosquitoes. PLoS One 6:e20579. doi:10.1371/journal.pone.0020579.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Bockarie MJ
    . 2007. Molecular xenomonitoring of lymphatic filariasis. Am J Trop Med Hyg 77:591–592. doi:10.4269/ajtmh.2007.77.591.
    OpenUrlFREE Full Text
  10. 10.↵
    1. Fernandez de Marco M,
    2. Brugman VA,
    3. Hernandez-Triana LM,
    4. Thorne L,
    5. Phipps LP,
    6. Nikolova NI,
    7. Fooks AR,
    8. Johnson N
    . 2016. Detection of Theileria orientalis in mosquito blood meals in the United Kingdom. Vet Parasitol 229:31–36. doi:10.1016/j.vetpar.2016.09.012.
    OpenUrlCrossRef
  11. 11.↵
    1. Rao RU,
    2. Samarasekera SD,
    3. Nagodavithana KC,
    4. Punchihewa MW,
    5. Dassanayaka TDM,
    6. Gamini PKD,
    7. Ford E,
    8. Ranasinghe USB,
    9. Henderson RH,
    10. Weil GJ
    . 2016. Programmatic use of molecular xenomonitoring at the level of evaluation units to assess persistence of lymphatic filariasis in Sri Lanka. PLoS Negl Trop Dis 10:e0004722. doi:10.1371/journal.pntd.0004722.
    OpenUrlCrossRef
  12. 12.↵
    1. Souf S
    . 2016. Recent advances in diagnostic testing for viral infections. Biosci Horiz 9:hzw010. doi:10.1093/biohorizons/hzw010.
    OpenUrlCrossRef
  13. 13.↵
    1. Fontaine A,
    2. Jiolle D,
    3. Moltini-Conclois I,
    4. Lequime S,
    5. Lambrechts L
    . 2016. Excretion of dengue virus RNA by Aedes aegypti allows non-destructive monitoring of viral dissemination in individual mosquitoes. Sci Rep 6:24885. doi:10.1038/srep24885.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Ramírez AL,
    2. Hall-Mendelin S,
    3. Doggett SL,
    4. Hewitson GR,
    5. McMahon JL,
    6. Ritchie SA,
    7. van den Hurk AF
    . 2018. Mosquito excreta: a sample type with many potential applications for the investigation of Ross River virus and West Nile virus ecology. PLoS Negl Trop Dis 12:e0006771. doi:10.1371/journal.pntd.0006771.
    OpenUrlCrossRef
  15. 15.↵
    1. Blow JA,
    2. Turell MJ,
    3. Walker ED,
    4. Silverman AL
    . 2002. Post-bloodmeal diuretic shedding of hepatitis B virus by mosquitoes (Diptera: Culicidae). J Med Entomol 39:605–612. doi:10.1603/0022-2585-39.4.605.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Meyer DB,
    2. Ramírez AL,
    3. van den Hurk AF,
    4. Kurucz N,
    5. Ritchie SA
    . 2019. Development and field evaluation of a system to collect mosquito excreta for the detection of arboviruses. J Med Entomol 56:1116–1121. doi:10.1093/jme/tjz031.
    OpenUrlCrossRef
  17. 17.↵
    1. Jansen CC,
    2. Prow NA,
    3. Webb CE,
    4. Hall RA,
    5. Pyke AT,
    6. Harrower BJ,
    7. Pritchard IL,
    8. Zborowski P,
    9. Ritchie SA,
    10. Russell RC,
    11. Van Den Hurk AF
    . 2009. Arboviruses isolated from mosquitoes collected from urban and peri-urban areas of eastern Australia. J Am Mosq Control Assoc 25:272–278. doi:10.2987/09-5908.1.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Roche SE,
    2. Wicks R,
    3. Garner MG,
    4. East IJ,
    5. Paskin R,
    6. Moloney BJ,
    7. Carr M,
    8. Kirkland P
    . 2013. Descriptive overview of the 2011 epidemic of arboviral disease in horses in Australia. Aust Vet J 91:5–13. doi:10.1111/avj.12018.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Colmant AMG,
    2. Etebari K,
    3. Webb CE,
    4. Ritchie SA,
    5. Jansen CC,
    6. van den Hurk AF,
    7. Bielefeldt-Ohmann H,
    8. Hobson-Peters J,
    9. Asgari S,
    10. Hall RA
    . 2017. Discovery of new orbiviruses and totivirus from Anopheles mosquitoes in Eastern Australia. Arch Virol 162:3529–3534. doi:10.1007/s00705-017-3515-x.
    OpenUrlCrossRef
  20. 20.↵
    1. Colmant AMG,
    2. Hobson-Peters J,
    3. Bielefeldt-Ohmann H,
    4. van den Hurk AF,
    5. Hall-Mendelin S,
    6. Chow WK,
    7. Johansen CA,
    8. Fros J,
    9. Simmonds P,
    10. Watterson D,
    11. Cazier C,
    12. Etebari K,
    13. Asgari S,
    14. Schulz BL,
    15. Beebe N,
    16. Vet LJ,
    17. Piyasena TBH,
    18. Nguyen H-D,
    19. Barnard RT,
    20. Hall RA
    . 2017. A new clade of insect-specific flaviviruses from Australian Anopheles mosquitoes displays species-specific host restriction. mSphere 2:e00262-17. doi:10.1128/mSphere.00262-17.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Huang B,
    2. Jennison A,
    3. Whiley D,
    4. McMahon J,
    5. Hewitson G,
    6. Graham R,
    7. De Jong A,
    8. Warrilow D
    . 2019. Illumina sequencing of clinical samples for virus detection in a public health laboratory. Sci Rep 9:5409. doi:10.1038/s41598-019-41830-w.
    OpenUrlCrossRef
  22. 22.↵
    1. Hall RA,
    2. Prow NA,
    3. Pyke AT
    . 2010. Molecular diagnostics for Ross River virus. In Liu D (ed), Molecular detection of human viral pathogens. CRC Press, Boca Raton, FL.
  23. 23.↵
    1. van den Hurk AF,
    2. Hall-Mendelin S,
    3. Townsend M,
    4. Kurucz N,
    5. Edwards J,
    6. Ehlers G,
    7. Rodwell C,
    8. Moore FA,
    9. McMahon JL,
    10. Northill JA,
    11. Simmons RJ,
    12. Cortis G,
    13. Melville L,
    14. Whelan PI,
    15. Ritchie SA
    . 2014. Applications of a sugar-based surveillance system to track arboviruses in wild mosquito populations. Vector Borne Zoonotic Dis 14:66–73. doi:10.1089/vbz.2013.1373.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Shi M,
    2. Lin X-D,
    3. Tian J-H,
    4. Chen L-J,
    5. Chen X,
    6. Li C-X,
    7. Qin X-C,
    8. Li J,
    9. Cao J-P,
    10. Eden J-S,
    11. Buchmann J,
    12. Wang W,
    13. Xu J,
    14. Holmes EC,
    15. Zhang Y-Z
    . 2016. Redefining the invertebrate RNA virosphere. Nature 540:539–543. doi:10.1038/nature20167.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Kobayashi D,
    2. Isawa H,
    3. Fujita R,
    4. Murota K,
    5. Itokawa K,
    6. Higa Y,
    7. Katayama Y,
    8. Sasaki T,
    9. Mizutani T,
    10. Iwanaga S,
    11. Ohta N,
    12. Garcia-Bertuso A,
    13. Sawabe K
    . 2017. Isolation and characterization of a new iflavirus from Armigeres spp. mosquitoes in the Philippines. J Gen Virol 98:2876–2881. doi:10.1099/jgv.0.000929.
    OpenUrlCrossRef
  26. 26.↵
    1. Hang J,
    2. Klein TA,
    3. Kim H-C,
    4. Yang Y,
    5. Jima DD,
    6. Richardson JH,
    7. Jarman RG
    . 2016. Genome sequences of five arboviruses in field-captured mosquitoes in a unique rural environment of South Korea. Genome Announc 4:e01644-15. doi:10.1128/genomeA.01644-15.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Sadeghi M,
    2. Altan E,
    3. Deng X,
    4. Barker CM,
    5. Fang Y,
    6. Coffey LL,
    7. Delwart E
    . 2018. Virome of >12 thousand Culex mosquitoes from throughout California. Virology 523:74–88. doi:10.1016/j.virol.2018.07.029.
    OpenUrlCrossRef
  28. 28.↵
    1. Quast C,
    2. Pruesse E,
    3. Yilmaz P,
    4. Gerken J,
    5. Schweer T,
    6. Yarza P,
    7. Peplies J,
    8. Glöckner FO
    . 2013. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41:D590–D596. doi:10.1093/nar/gks1219.
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.↵
    1. Forbes JD,
    2. Knox NC,
    3. Ronholm J,
    4. Pagotto F,
    5. Reimer A
    . 2017. Metagenomics: the next culture-independent game changer. Front Microbiol 8:1069. doi:10.3389/fmicb.2017.01069.
    OpenUrlCrossRef
  30. 30.↵
    1. Zhang Y-Z,
    2. Chen Y-M,
    3. Wang W,
    4. Qin X-C,
    5. Holmes EC
    . 2019. Expanding the RNA virosphere by unbiased metagenomics. Annu Rev Virol 6:119–139. doi:10.1146/annurev-virology-092818-015851.
    OpenUrlCrossRef
  31. 31.↵
    1. Wylie KM,
    2. Mihindukulasuriya KA,
    3. Sodergren E,
    4. Weinstock GM,
    5. Storch GA
    . 2012. Sequence analysis of the human virome in febrile and afebrile children. PLoS One 7:e27735. doi:10.1371/journal.pone.0027735.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Batovska J,
    2. Mee PT,
    3. Lynch SE,
    4. Sawbridge TI,
    5. Rodoni BC
    . 2019. Sensitivity and specificity of metatranscriptomics as an arbovirus surveillance tool. Sci Rep 9:19398. doi:10.1038/s41598-019-55741-3.
    OpenUrlCrossRef
  33. 33.↵
    1. Conceição-Neto N,
    2. Zeller M,
    3. Lefrère H,
    4. De Bruyn P,
    5. Beller L,
    6. Deboutte W,
    7. Yinda CK,
    8. Lavigne R,
    9. Maes P,
    10. Van Ranst M,
    11. Heylen E,
    12. Matthijnssens J
    . 2015. Modular approach to customise sample preparation procedures for viral metagenomics: a reproducible protocol for virome analysis. Sci Rep 5:16532. doi:10.1038/srep16532.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Shi C,
    2. Beller L,
    3. Deboutte W,
    4. Yinda KC,
    5. Delang L,
    6. Vega-Rúa A,
    7. Failloux A-B,
    8. Matthijnssens J
    . 2019. Stable distinct core eukaryotic viromes in different mosquito species from Guadeloupe, using single mosquito viral metagenomics. Microbiome 7:121–121. doi:10.1186/s40168-019-0734-2.
    OpenUrlCrossRef
  35. 35.↵
    1. Birnberg L,
    2. Temmam S,
    3. Aranda C,
    4. Correa-Fiz F,
    5. Talavera S,
    6. Bigot T,
    7. Eloit M,
    8. Busquets N
    . 2020. Viromics on honey-baited FTA cards as a new tool for the detection of circulating viruses in mosquitoes. Viruses 12:274. doi:10.3390/v12030274.
    OpenUrlCrossRef
  36. 36.↵
    1. Bailey L,
    2. Woods RD
    . 1974. Three previously undescribed viruses from the honey bee. J Gen Virol 25:175–186. doi:10.1099/0022-1317-25-2-175.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Isawa H,
    2. Asano S,
    3. Sahara K,
    4. Iizuka T,
    5. Bando H
    . 1998. Analysis of genetic information of an insect picorna-like virus, infectious flacherie virus of silkworm: evidence for evolutionary relationships among insect, mammalian and plant picorna(-like) viruses. Arch Virol 143:127–143. doi:10.1007/s007050050273.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Vasilakis N,
    2. Tesh RB
    . 2015. Insect-specific viruses and their potential impact on arbovirus transmission. Curr Opin Virol 15:69–74. doi:10.1016/j.coviro.2015.08.007.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Hall RA,
    2. Bielefeldt-Ohmann H,
    3. McLean BJ,
    4. O'Brien CA,
    5. Colmant AMG,
    6. Piyasena TBH,
    7. Harrison JJ,
    8. Newton ND,
    9. Barnard RT,
    10. Prow NA,
    11. Deerain JM,
    12. Mah MGKY,
    13. Hobson-Peters J
    . 2016. Commensal viruses of mosquitoes: host restriction, transmission, and interaction with arboviral pathogens. Evol Bioinform Online 12:35–44. doi:10.4137/EBO.S40740.
    OpenUrlCrossRef
  40. 40.↵
    1. Li C-X,
    2. Smith ML,
    3. Fulcher A,
    4. Kaufman PE,
    5. Zhao T-Y,
    6. Xue R-D
    . 2015. Field evaluation of three new mosquito light traps against two standard light traps to collect mosquitoes (Diptera: Culicidae) and non-target insects in Northeast Florida. Fla Entomol 98:114–117. doi:10.1653/024.098.0118.
    OpenUrlCrossRef
  41. 41.↵
    1. Ritchie SA,
    2. Cortis G,
    3. Paton C,
    4. Townsend M,
    5. Shroyer D,
    6. Zborowski P,
    7. Hall-Mendelin S,
    8. van den Hurk AF
    . 2013. A simple non-powered passive trap for the collection of mosquitoes for arbovirus surveillance. J Med Entomol 50:185–194. doi:10.1603/me12112.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Rinker DC,
    2. Pitts RJ,
    3. Zwiebel LJ
    . 2016. Disease vectors in the era of next generation sequencing. Genome Biol 17:95. doi:10.1186/s13059-016-0966-4.
    OpenUrlCrossRef
  43. 43.↵
    1. Goddard LB,
    2. Roth AE,
    3. Reisen WK,
    4. Scott TW
    . 2002. Vector competence of California mosquitoes for West Nile virus. Emerg Infect Dis 8:1385–1391. doi:10.3201/eid0812.020536.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. Ritchie SA,
    2. Kline DL
    . 1995. Comparison of CDC and EVS light traps baited with carbon dioxide and octenol for trapping mosquitoes in Brisbane, Queensland (Diptera: Culicidae). Aust J Entomol 34:215–218. doi:10.1111/j.1440-6055.1995.tb01322.x.
    OpenUrlCrossRef
  45. 45.↵
    1. Pyke AT,
    2. Smith IL,
    3. van den Hurk AF,
    4. Northill JA,
    5. Chuan TF,
    6. Westacott AJ,
    7. Smith GA
    . 2004. Detection of Australasian Flavivirus encephalitic viruses using rapid fluorogenic TaqMan RT-PCR assays. J Virol Methods 117:161–167. doi:10.1016/j.jviromet.2004.01.007.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    1. Buchfink B,
    2. Xie C,
    3. Huson DH
    . 2015. Fast and sensitive protein alignment using DIAMOND. Nat Methods 12:59–60. doi:10.1038/nmeth.3176.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Huson DH,
    2. Beier S,
    3. Flade I,
    4. Gorska A,
    5. El-Hadidi M,
    6. Mitra S,
    7. Ruscheweyh H-J,
    8. Tappu R
    . 2016. MEGAN community edition—interactive exploration and analysis of large-scale microbiome sequencing data. PLoS Comput Biol 12:e1004957. doi:10.1371/journal.pcbi.1004957.
    OpenUrlCrossRef
  48. 48.↵
    1. Bankevich A,
    2. Nurk S,
    3. Antipov D,
    4. Gurevich AA,
    5. Dvorkin M,
    6. Kulikov AS,
    7. Lesin VM,
    8. Nikolenko SI,
    9. Pham S,
    10. Prjibelski AD,
    11. Pyshkin AV,
    12. Sirotkin AV,
    13. Vyahhi N,
    14. Tesler G,
    15. Alekseyev MA,
    16. Pevzner PA
    . 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455–477. doi:10.1089/cmb.2012.0021.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Johnson M,
    2. Zaretskaya I,
    3. Raytselis Y,
    4. Merezhuk Y,
    5. McGinnis S,
    6. Madden TL
    . 2008. NCBI BLAST: a better web interface. Nucleic Acids Res 36:W5–W9. doi:10.1093/nar/gkn201.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    1. Katoh K,
    2. Standley DM
    . 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30:772–780. doi:10.1093/molbev/mst010.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Capella-Gutierrez S,
    2. Silla-Martinez JM,
    3. Gabaldon T
    . 2009. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25:1972–1973. doi:10.1093/bioinformatics/btp348.
    OpenUrlCrossRefPubMedWeb of Science
  52. 52.↵
    1. Lefort V,
    2. Longueville J-E,
    3. Gascuel O
    . 2017. SMS: Smart Model Selection in PhyML. Mol Biol Evol 34:2422–2424. doi:10.1093/molbev/msx149.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Guindon S,
    2. Dufayard J-F,
    3. Lefort V,
    4. Anisimova M,
    5. Hordijk W,
    6. Gascuel O
    . 2010. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 59:307–321. doi:10.1093/sysbio/syq010.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Metagenomic Analysis of the Virome of Mosquito Excreta
Ana L. Ramírez, Agathe M. G. Colmant, David Warrilow, Bixing Huang, Alyssa T. Pyke, Jamie L. McMahon, Dagmar B. Meyer, Rikki M. A. Graham, Amy V. Jennison, Scott A. Ritchie, Andrew F. van den Hurk
mSphere Sep 2020, 5 (5) e00587-20; DOI: 10.1128/mSphere.00587-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSphere article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Metagenomic Analysis of the Virome of Mosquito Excreta
(Your Name) has forwarded a page to you from mSphere
(Your Name) thought you would be interested in this article in mSphere.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Metagenomic Analysis of the Virome of Mosquito Excreta
Ana L. Ramírez, Agathe M. G. Colmant, David Warrilow, Bixing Huang, Alyssa T. Pyke, Jamie L. McMahon, Dagmar B. Meyer, Rikki M. A. Graham, Amy V. Jennison, Scott A. Ritchie, Andrew F. van den Hurk
mSphere Sep 2020, 5 (5) e00587-20; DOI: 10.1128/mSphere.00587-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

excreta
metagenomics
mosquito
next-generation sequencing
virome

Related Articles

Cited By...

About

  • About mSphere
  • Board of Editors
  • Policies
  • For Reviewers
  • For the Media
  • Embargo Policy
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSphereJ

@ASMicrobiology

       

 

Website feedback

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5042