Stability of the Virome in Lab- and Field-Collected Aedes albopictus Mosquitoes across Different Developmental Stages and Possible Core Viruses in the Publicly Available Virome Data of Aedes Mosquitoes

Our study revealed that the virome was very stable across all developmental stages of both lab-derived and field-collected Aedes albopictus. The data representing the core virome in lab A. albopictus proved the vertical transmission route of these viruses, forming a “vertically transmitted core virome.” Field mosquitoes also contained this stable vertically transmitted core virome as well as additional viruses, which probably represented “environment-derived core virome” and which therefore were less stable over time and geography. By further screening publicly available SRA viral metagenomic data sets from mosquitoes belonging to the genus Aedes, some of the identified core ISVs were shown to be present in the majority of SRAs, such as Phasi Charoen-like phasivirus and Guadeloupe mosquito virus. How these core ISVs influence the biology of the mosquito host and arbovirus infection and evolution deserves to be further explored. ABSTRACT Aedes mosquitoes can efficiently transmit many pathogenic arboviruses, placing a great burden on public health worldwide. In addition, they also carry a number of insect-specific viruses (ISVs), and it was recently suggested that some of these ISVs might form a stable species-specific “core virome” in mosquito populations. However, little is known about such a core virome in laboratory colonies and if it is present across different developmental stages. In this study, we compared the viromes in eggs, larvae, pupae, and adults of Aedes albopictus mosquitoes collected from a lab colony and compared each to the virome of different developmental stages collected in the field. The virome in lab-derived A. albopictus was very stable across all stages, consistent with a vertical transmission route of these viruses, and formed a possible “vertically transmitted core virome.” The different stages of field-collected A. albopictus mosquitoes also contained this stable vertically transmitted core virome, as well as another set of viruses (e.g., viruses distantly related to Guadeloupe mosquito virus, Hubei virga-like virus 2, and Sarawak virus) shared by mosquitoes across different stages, which might represent an “environment-derived core virome.” To further study this core set of ISVs, we screened 48 publicly available SRA viral metagenomic data sets of mosquitoes belonging to the genus Aedes, showing that some of the identified ISVs were identified in the majority of SRAs and providing further evidence supporting the core-virome concept. IMPORTANCE Our study revealed that the virome was very stable across all developmental stages of both lab-derived and field-collected Aedes albopictus. The data representing the core virome in lab A. albopictus proved the vertical transmission route of these viruses, forming a “vertically transmitted core virome.” Field mosquitoes also contained this stable vertically transmitted core virome as well as additional viruses, which probably represented “environment-derived core virome” and which therefore were less stable over time and geography. By further screening publicly available SRA viral metagenomic data sets from mosquitoes belonging to the genus Aedes, some of the identified core ISVs were shown to be present in the majority of SRAs, such as Phasi Charoen-like phasivirus and Guadeloupe mosquito virus. How these core ISVs influence the biology of the mosquito host and arbovirus infection and evolution deserves to be further explored.

[1]  R. Irizarry ggplot2 , 2019, Introduction to Data Science.

[2]  A. Blomström,et al.  Viromics Reveal a Number of Novel RNA Viruses in Swedish Mosquitoes , 2019, Viruses.

[3]  Kwe Claude Yinda,et al.  Stable distinct core eukaryotic viromes in different mosquito species from Guadeloupe, using single mosquito viral metagenomics , 2019, Microbiome.

[4]  M. Shi,et al.  Meta-Transcriptomic Comparison of the RNA Viromes of the Mosquito Vectors Culex pipiens and Culex torrentium in Northern Europe , 2019, bioRxiv.

[5]  A. Fernández-Sesma,et al.  Aedes aegypti (Aag2)-derived clonal mosquito cell lines reveal the impact of pre-existing persistent infection with the insect-specific bunyavirus Phasi Charoen-like virus on arbovirus replication , 2019, bioRxiv.

[6]  R. Jarman,et al.  Metagenomic Analysis Reveals Three Novel and Prevalent Mosquito Viruses from a Single Pool of Aedes vexans nipponii Collected in the Republic of Korea , 2019, Viruses.

[7]  J. Ribeiro,et al.  Aedes aegypti(Aag2)-derived clonal mosquito cell lines reveal impact of pre-existing persistent co-infection with the insect-specific Bunyavirus Phasi Charoen-like virus on arbovirus replication , 2019, Access Microbiology.

[8]  Ahmad Feroz Changing Epidemiology of Dengue , 2018, Journal of Shaheed Suhrawardy Medical College.

[9]  D. Brackney,et al.  Generation of axenic Aedes aegypti demonstrate live bacteria are not required for mosquito development , 2018, Nature Communications.

[10]  H. Liu,et al.  Metagenomic Analysis of Flaviviridae in Mosquito Viromes Isolated From Yunnan Province in China Reveals Genes From Dengue and Zika Viruses , 2018, Front. Cell. Infect. Microbiol..

[11]  Alexandre Souvorov,et al.  SKESA: strategic k-mer extension for scrupulous assemblies , 2018, Genome Biology.

[12]  Jinquan Cheng,et al.  Discovery and high prevalence of Phasi Charoen-like virus in field-captured Aedes aegypti in South China. , 2018, Virology.

[13]  C. Barker,et al.  Virome of > 12 thousand Culex mosquitoes from throughout California. , 2018, Virology.

[14]  J. Connor,et al.  Dual Insect specific virus infection limits Arbovirus replication in Aedes mosquito cells. , 2018, Virology.

[15]  C. Valiente Moro,et al.  The mosquito holobiont: fresh insight into mosquito-microbiota interactions , 2018, Microbiome.

[16]  L. Krause,et al.  Mapping the virome in wild-caught Aedes aegypti from Cairns and Bangkok , 2018, Scientific Reports.

[17]  Chenyan Shi,et al.  Comparative Metagenomic Profiling of Viromes Associated with Four Common Mosquito Species in China , 2018, Virologica Sinica.

[18]  A. Zohaib,et al.  Metagenomic Virome Analysis of Culex Mosquitoes from Kenya and China , 2018, Viruses.

[19]  M. Shi,et al.  High-Resolution Metatranscriptomics Reveals the Ecological Dynamics of Mosquito-Associated RNA Viruses in Western Australia , 2017, Journal of Virology.

[20]  Thomas K. F. Wong,et al.  ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates , 2017, Nature Methods.

[21]  Alexander Lex,et al.  UpSetR: an R package for the visualization of intersecting sets and their properties , 2017, bioRxiv.

[22]  Steven M. Lakin,et al.  West African Anopheles gambiae mosquitoes harbor a taxonomically diverse virome including new insect-specific flaviviruses, mononegaviruses, and totiviruses. , 2016, Virology.

[23]  Roland Eils,et al.  Complex heatmaps reveal patterns and correlations in multidimensional genomic data , 2016, Bioinform..

[24]  Mark R. Brown,et al.  Gut bacteria differentially affect egg production in the anautogenous mosquito Aedes aegypti and facultatively autogenous mosquito Aedes atropalpus (Diptera: Culicidae) , 2016, Parasites & Vectors.

[25]  Xing Li,et al.  Characterizing a large outbreak of dengue fever in Guangdong Province, China , 2016, Infectious Diseases of Poverty.

[26]  Vishwesh P. Mokashi,et al.  Bioinformatic Characterization of Mosquito Viromes within the Eastern United States and Puerto Rico: Discovery of Novel Viruses , 2016, Evolutionary bioinformatics online.

[27]  R. Tesh,et al.  Insect-Specific Virus Discovery: Significance for the Arbovirus Community , 2015, Viruses.

[28]  H. Lanz-Mendoza,et al.  Bacterial Exposure at the Larval Stage Induced Sexual Immune Dimorphism and Priming in Adult Aedes aegypti Mosquitoes , 2015, PloS one.

[29]  Xiaomin Hu,et al.  A Metagenomic Survey of Viral Abundance and Diversity in Mosquitoes from Hubei Province , 2015, PloS one.

[30]  Andrew J Tatem,et al.  The changing epidemiology of dengue in China, 1990-2014: a descriptive analysis of 25 years of nationwide surveillance data , 2015, BMC Medicine.

[31]  H. Guzmán,et al.  Insect-specific viruses detected in laboratory mosquito colonies and their potential implications for experiments evaluating arbovirus vector competence. , 2015, The American journal of tropical medicine and hygiene.

[32]  M. Shi,et al.  Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses , 2015, eLife.

[33]  M. Oshaghi,et al.  Malpighian tubules are important determinants of Pseudomonas transstadial transmission and longtime persistence in Anopheles stephensi , 2015, Parasites & Vectors.

[34]  E. Walker,et al.  Elizabethkingia anophelis: Molecular Manipulation and Interactions with Mosquito Hosts , 2015, Applied and Environmental Microbiology.

[35]  R. Christen,et al.  Composition of Anopheles coluzzii and Anopheles gambiae microbiota from larval to adult stages. , 2014, Infection, genetics and evolution : journal of molecular epidemiology and evolutionary genetics in infectious diseases.

[36]  Chao Xie,et al.  Fast and sensitive protein alignment using DIAMOND , 2014, Nature Methods.

[37]  A. von Haeseler,et al.  IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies , 2014, Molecular biology and evolution.

[38]  Guofa Zhou,et al.  Urbanization Increases Aedes albopictus Larval Habitats and Accelerates Mosquito Development and Survivorship , 2014, PLoS neglected tropical diseases.

[39]  M. Sanchez Leon,et al.  Virome Analysis of Amblyomma americanum, Dermacentor variabilis, and Ixodes scapularis Ticks Reveals Novel Highly Divergent Vertebrate and Invertebrate Viruses , 2014, Journal of Virology.

[40]  Mark R. Brown,et al.  Mosquitoes rely on their gut microbiota for development , 2014, Molecular ecology.

[41]  J. Bruenn,et al.  Discovery and Evolution of Bunyavirids in Arctic Phantom Midges and Ancient Bunyavirid-Like Sequences in Insect Genomes , 2014, Journal of Virology.

[42]  M. Shi,et al.  A tick-borne segmented RNA virus contains genome segments derived from unsegmented viral ancestors , 2014, Proceedings of the National Academy of Sciences.

[43]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[44]  W. Lipkin,et al.  Genome characterization of Long Island tick rhabdovirus, a new virus identified in Amblyomma americanum ticks , 2014, Virology Journal.

[45]  A. James,et al.  The invasive mosquito species Aedes albopictus: current knowledge and future perspectives. , 2013, Trends in parasitology.

[46]  Susan Holmes,et al.  phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data , 2013, PloS one.

[47]  Sergey I. Nikolenko,et al.  SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing , 2012, J. Comput. Biol..

[48]  F. Simard,et al.  Notes on the blood-feeding behavior of Aedes albopictus (Diptera: Culicidae) in Cameroon , 2012, Parasites & Vectors.

[49]  Xiao-Guang Chen,et al.  A local outbreak of dengue caused by an imported case in Dongguan China , 2012, BMC Public Health.

[50]  Adam M. Phillippy,et al.  Interactive metagenomic visualization in a Web browser , 2011, BMC Bioinformatics.

[51]  H. Delatte,et al.  Blood-feeding behavior of Aedes albopictus, a vector of Chikungunya on La Réunion. , 2010, Vector borne and zoonotic diseases.

[52]  C. Paupy,et al.  Aedes albopictus, an arbovirus vector: from the darkness to the light. , 2009, Microbes and infection.

[53]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[54]  Toni Gabaldón,et al.  trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses , 2009, Bioinform..

[55]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[56]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[57]  L. Kramer,et al.  Bacteria of the genus Asaia stably associate with Anopheles stephensi, an Asian malarial mosquito vector , 2007, Proceedings of the National Academy of Sciences.

[58]  K. Katoh,et al.  MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. , 2002, Nucleic acids research.

[59]  Z. Adelman,et al.  Sindbis virus‐induced silencing of dengue viruses in mosquitoes , 2001, Insect molecular biology.

[60]  W. Hawley,et al.  Aedes albopictus (Diptera: Culicidae): a statistical model of the role of temperature, photoperiod, and geography in the induction of egg diapause. , 1994, Journal of medical entomology.

[61]  M. N. Anwar,et al.  A Viral Metagenomic Analysis Reveals Rich Viral Abundance and Diversity in Mosquitoes from Pig Farms. , 2019, Transboundary and emerging diseases.

[62]  K. Lerdthusnee,et al.  Meconial Peritrophic Membranes and the Fate of Midgut Bacteria During Mosquito (Diptera: Culicidae) Metamorphosis , 2001, Journal of medical entomology.