Ribosomal RNA (rRNA) sequences from 33 globally distributed mosquito species for improved metagenomics and species identification

  1. Cassandra Koh  Is a corresponding author
  2. Lionel Frangeul
  3. Hervé Blanc
  4. Carine Ngoagouni
  5. Sébastien Boyer
  6. Philippe Dussart
  7. Nina Grau
  8. Romain Girod
  9. Jean-Bernard Duchemin
  10. Maria-Carla Saleh  Is a corresponding author
  1. Institut Pasteur, France
  2. Institut Pasteur de Bangui, Central African Republic
  3. Institut Pasteur du Cambodge, Cambodia
  4. Institut Pasteur de Madagascar, Madagascar
  5. Institut Pasteur de la Guyane, French Guiana

Abstract

Total RNA sequencing (RNA-seq) is an important tool in the study of mosquitoes and the RNA viruses they vector as it allows assessment of both host and viral RNA in specimens. However, there are two main constraints. First, as with many other species, abundant mosquito ribosomal RNA (rRNA) serves as the predominant template from which sequences are generated, meaning that the desired host and viral templates are sequenced far less. Second, mosquito specimens captured in the field must be correctly identified, in some cases to the sub-species level. Here, we generate mosquito ribosomal RNA (rRNA) datasets which will substantially mitigate both of these problems. We describe a strategy to assemble novel rRNA sequences from mosquito specimens and produce an unprecedented dataset of 234 full-length 28S and 18S rRNA sequences of 33 medically important species from countries with known histories of mosquito-borne virus circulation (Cambodia, the Central African Republic, Madagascar, and French Guiana). These sequences will allow both physical and computational removal of rRNA from specimens during RNAseq protocols. We also assess the utility of rRNA sequences for molecular taxonomy and compare phylogenies constructed using rRNA sequences versus those created using the gold standard for molecular species identification of specimens-the mitochondrial cytochrome c oxidase I (COI) gene. We find that rRNA- and COI-derived phylogenetic trees are incongruent and that 28S and concatenated 28S+18S rRNA phylogenies reflect evolutionary relationships that are more aligned with contemporary mosquito systematics. This significant expansion to the current rRNA reference library for mosquitoes will improve mosquito RNA-seq metagenomics by permitting the optimization of species-specific rRNA depletion protocols for a broader range of species and streamlining species identification by rRNA sequence and phylogenetics.

Data availability

Multiple sequence alignment files are included as source data files. All sequences generated in this study have been deposited in GenBank under the accession numbers OM350214-OM350327 for 18S rRNA sequences, OM542339-OM542460 for 28S rRNA sequences, and OM630610-OM630715 for COI sequences.

Article and author information

Author details

  1. Cassandra Koh

    Viruses and RNA Interference Unit, Institut Pasteur, Paris, France
    For correspondence
    cassandra.koh@pasteur.fr
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2466-6731
  2. Lionel Frangeul

    Viruses and RNA Interference Unit, Institut Pasteur, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  3. Hervé Blanc

    Viruses and RNA Interference Unit, Institut Pasteur, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  4. Carine Ngoagouni

    Medical Entomology Laboratory, Institut Pasteur de Bangui, Bangui, Central African Republic
    Competing interests
    The authors declare that no competing interests exist.
  5. Sébastien Boyer

    Medical and Veterinary Entomology Unit, Institut Pasteur du Cambodge, Phnom Penh, Cambodia
    Competing interests
    The authors declare that no competing interests exist.
  6. Philippe Dussart

    Virology Unit, Institut Pasteur du Cambodge, Phnom Penh, Cambodia
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1931-3037
  7. Nina Grau

    Medical Entomology Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar
    Competing interests
    The authors declare that no competing interests exist.
  8. Romain Girod

    Medical Entomology Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar
    Competing interests
    The authors declare that no competing interests exist.
  9. Jean-Bernard Duchemin

    Vectopôle Amazonien Emile Abonnenc, Institut Pasteur de la Guyane, Cayenne, French Guiana
    Competing interests
    The authors declare that no competing interests exist.
  10. Maria-Carla Saleh

    Viruses and RNA Interference Unit, Institut Pasteur, Paris, France
    For correspondence
    carla.saleh@pasteur.fr
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8593-4117

Funding

Defense Advanced Research Projects Agency (Cooperative Agreement HR001118S0017)

  • Maria-Carla Saleh

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

Copyright

© 2023, Koh 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

  • 1,454
    views
  • 207
    downloads
  • 3
    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. Cassandra Koh
  2. Lionel Frangeul
  3. Hervé Blanc
  4. Carine Ngoagouni
  5. Sébastien Boyer
  6. Philippe Dussart
  7. Nina Grau
  8. Romain Girod
  9. Jean-Bernard Duchemin
  10. Maria-Carla Saleh
(2023)
Ribosomal RNA (rRNA) sequences from 33 globally distributed mosquito species for improved metagenomics and species identification
eLife 12:e82762.
https://doi.org/10.7554/eLife.82762

Share this article

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

Further reading

    1. Biochemistry and Chemical Biology
    2. Microbiology and Infectious Disease
    Stephanie M Stuteley, Ghader Bashiri
    Insight

    In the bacterium M. smegmatis, an enzyme called MftG allows the cofactor mycofactocin to transfer electrons released during ethanol metabolism to the electron transport chain.

    1. Microbiology and Infectious Disease
    Chenghao Jia, Chenghu Huang ... Min Yue
    Research Article

    Bacterial regional demonstration after global dissemination is an essential pathway for selecting distinct finesses. However, the evolution of the resistome during the transition to endemicity remains unaddressed. Using the most comprehensive whole-genome sequencing dataset of Salmonella enterica serovar Gallinarum (S. Gallinarum) collected from 15 countries, including 45 newly recovered samples from two related local regions, we established the relationship among avian-specific pathogen genetic profiles and localization patterns. Initially, we revealed the international transmission and evolutionary history of S. Gallinarum to recent endemicity through phylogenetic analysis conducted using a spatiotemporal Bayesian framework. Our findings indicate that the independent acquisition of the resistome via the mobilome, primarily through plasmids and transposons, shapes a unique antimicrobial resistance profile among different lineages. Notably, the mobilome-resistome combination among distinct lineages exhibits a geographical-specific manner, further supporting a localized endemic mobilome-driven process. Collectively, this study elucidates resistome adaptation in the endemic transition of an avian-specific pathogen, likely driven by the localized farming style, and provides valuable insights for targeted interventions.