Single mosquito metatranscriptomics identifies vectors, emerging pathogens and reservoirs in one assay

  1. Joshua Batson
  2. Gytis Dudas
  3. Eric Haas-Stapleton
  4. Amy L Kistler  Is a corresponding author
  5. Lucy M Li
  6. Phoenix Logan
  7. Kalani Ratnasiri
  8. Hanna Retallack
  1. Chan Zuckerberg Biohub, United States
  2. Gothenburg Global Biodiversity Centre, Sweden
  3. Alameda County Mosquito Abatement District, United States
  4. Stanford University, United States
  5. University of California, San Francisco, United States

Abstract

Mosquitoes are major infectious disease-carrying vectors. Assessment of current and future risks associated with the mosquito population requires knowledge of the full repertoire of pathogens they carry, including novel viruses, as well as their blood meal sources. Unbiased metatranscriptomic sequencing of individual mosquitoes offers a straightforward, rapid and quantitative means to acquire this information. Here, we profile 148 diverse wild-caught mosquitoes collected in California and detect sequences from eukaryotes, prokaryotes, 24 known and 46 novel viral species. Importantly, sequencing individuals greatly enhanced the value of the biological information obtained. It allowed us to a) speciate host mosquito, b) compute the prevalence of each microbe and recognize a high frequency of viral co-infections, c) associate animal pathogens with specific blood meal sources, and d) apply simple co-occurrence methods to recover previously undetected components of highly prevalent segmented viruses. In the context of emerging diseases, where knowledge about vectors, pathogens, and reservoirs is lacking, the approaches described here can provide actionable information for public health surveillance and intervention decisions.

Data availability

Raw and assembled sequencing data are deposited in NCBI Bioproject PRJNA605178. Code is available on Github at https://github.com/czbiohub/california-mosquito-study. Derived data (including all contigs) and supplementary data are available on Figshare at dx.doi.org/10.6084/m9.figshare.11832999 .

The following data sets were generated
The following previously published data sets were used

Article and author information

Author details

  1. Joshua Batson

    Data Science, Chan Zuckerberg Biohub, San Francisco, United States
    Competing interests
    The authors declare that no competing interests exist.
  2. Gytis Dudas

    Gothenburg Global Biodiversity Centre, Gothenburg, Sweden
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0227-4158
  3. Eric Haas-Stapleton

    Mosquito surveillance and control, Alameda County Mosquito Abatement District, Hayward, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Amy L Kistler

    Infectious Disease Initiative, Chan Zuckerberg Biohub, San Francisco, United States
    For correspondence
    amy.kistler@czbiohub.org
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1112-719X
  5. Lucy M Li

    Data Science, Chan Zuckerberg Biohub, San Francisco, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6562-4004
  6. Phoenix Logan

    Data Sciences, Chan Zuckerberg Biohub, San Francisco, United States
    Competing interests
    The authors declare that no competing interests exist.
  7. Kalani Ratnasiri

    Program in Immunology, Stanford University, Davis, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5953-0004
  8. Hanna Retallack

    University of California, San Francisco, San Francisco, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0533-9102

Funding

Chan Zuckerberg Initiative

  • Joshua Batson

UCSF Medical Science Training Program

  • Hanna Retallack

National Science Foundation

  • Kalani Ratnasiri

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. Dominique Soldati-Favre, University of Geneva, Switzerland

Version history

  1. Received: March 22, 2021
  2. Accepted: April 9, 2021
  3. Accepted Manuscript published: April 27, 2021 (version 1)
  4. Version of Record published: May 10, 2021 (version 2)

Copyright

© 2021, Batson et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 5,228
    views
  • 715
    downloads
  • 68
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Joshua Batson
  2. Gytis Dudas
  3. Eric Haas-Stapleton
  4. Amy L Kistler
  5. Lucy M Li
  6. Phoenix Logan
  7. Kalani Ratnasiri
  8. Hanna Retallack
(2021)
Single mosquito metatranscriptomics identifies vectors, emerging pathogens and reservoirs in one assay
eLife 10:e68353.
https://doi.org/10.7554/eLife.68353

Share this article

https://doi.org/10.7554/eLife.68353

Further reading

    1. Epidemiology and Global Health
    Yuchen Zhang, Yitang Sun ... Kaixiong Ye
    Research Article

    Background:

    Circulating omega-3 and omega-6 polyunsaturated fatty acids (PUFAs) have been associated with various chronic diseases and mortality, but results are conflicting. Few studies examined the role of omega-6/omega-3 ratio in mortality.

    Methods:

    We investigated plasma omega-3 and omega-6 PUFAs and their ratio in relation to all-cause and cause-specific mortality in a large prospective cohort, the UK Biobank. Of 85,425 participants who had complete information on circulating PUFAs, 6461 died during follow-up, including 2794 from cancer and 1668 from cardiovascular disease (CVD). Associations were estimated by multivariable Cox proportional hazards regression with adjustment for relevant risk factors.

    Results:

    Risk for all three mortality outcomes increased as the ratio of omega-6/omega-3 PUFAs increased (all Ptrend <0.05). Comparing the highest to the lowest quintiles, individuals had 26% (95% CI, 15–38%) higher total mortality, 14% (95% CI, 0–31%) higher cancer mortality, and 31% (95% CI, 10–55%) higher CVD mortality. Moreover, omega-3 and omega-6 PUFAs in plasma were all inversely associated with all-cause, cancer, and CVD mortality, with omega-3 showing stronger effects.

    Conclusions:

    Using a population-based cohort in UK Biobank, our study revealed a strong association between the ratio of circulating omega-6/omega-3 PUFAs and the risk of all-cause, cancer, and CVD mortality.

    Funding:

    Research reported in this publication was supported by the National Institute of General Medical Sciences of the National Institute of Health under the award number R35GM143060 (KY). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

    1. Ecology
    2. Epidemiology and Global Health
    Aleksandra Kovacevic, David RM Smith ... Lulla Opatowski
    Research Article

    Non-pharmaceutical interventions implemented to block SARS-CoV-2 transmission in early 2020 led to global reductions in the incidence of invasive pneumococcal disease (IPD). By contrast, most European countries reported an increase in antibiotic resistance among invasive Streptococcus pneumoniae isolates from 2019 to 2020, while an increasing number of studies reported stable pneumococcal carriage prevalence over the same period. To disentangle the impacts of the COVID-19 pandemic on pneumococcal epidemiology in the community setting, we propose a mathematical model formalizing simultaneous transmission of SARS-CoV-2 and antibiotic-sensitive and -resistant strains of S. pneumoniae. To test hypotheses underlying these trends five mechanisms were built into the model and examined: (1) a population-wide reduction of antibiotic prescriptions in the community, (2) lockdown effect on pneumococcal transmission, (3) a reduced risk of developing an IPD due to the absence of common respiratory viruses, (4) community azithromycin use in COVID-19 infected individuals, (5) and a longer carriage duration of antibiotic-resistant pneumococcal strains. Among 31 possible pandemic scenarios involving mechanisms individually or in combination, model simulations surprisingly identified only two scenarios that reproduced the reported trends in the general population. They included factors (1), (3), and (4). These scenarios replicated a nearly 50% reduction in annual IPD, and an increase in antibiotic resistance from 20% to 22%, all while maintaining a relatively stable pneumococcal carriage. Exploring further, higher SARS-CoV-2 R0 values and synergistic within-host virus-bacteria interaction mechanisms could have additionally contributed to the observed antibiotic resistance increase. Our work demonstrates the utility of the mathematical modeling approach in unraveling the complex effects of the COVID-19 pandemic responses on AMR dynamics.