Next Article in Journal
Mapping Attenuation Determinants in Enterovirus-D68
Next Article in Special Issue
Depicting the RNA Virome of Hematophagous Arthropods from Belgrade, Serbia
Previous Article in Journal
Indoleamine 2,3-Dioxygenase Is Involved in Interferon Gamma’s Anti-BKPyV Activity in Renal Cells
Previous Article in Special Issue
Discovery of Two Novel Negeviruses in a Dungfly Collected from the Arctic
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 8
10.3390/v12080866
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Aedes aegypti from Amazon Basin Harbor High Diversity of Novel Viral Species

by
Geovani de Oliveira Ribeiro
1,†,
Vanessa S. Morais
2,†,
Fred Julio Costa Monteiro
3,†,
Edcelha Soares D’Athaide Ribeiro
3,
Marlisson Octavio da S Rego
3,
Raimundo Nonato Picanço Souto
4,
Fabiola Villanova
1,
Roozbeh Tahmasebi
2,
Philip Michael Hefford
5,
Xutao Deng
6,7,
Eric Delwart
6,7,
Ester Cerdeira Sabino
2,
Licia Natal Fernandes
2,‡,
Antonio Charlys da Costa
2,*,‡ and
Élcio Leal
1,*,‡
1
Institute of Biological Sciences, Federal University of Pará, Belém, Pará, 66075-000, Brazil
2
Institute of Tropical Medicine, University of São Paulo, São Paulo 05403-000, Brazil
3
Public Health Laboratory of Amapa-LACEN/AP, Health Surveillance Superintendence of Amapa, Macapa 68905-230, Amapa, Brazil
4
Arthropoda Laboratory, Federal University of Amapa, Macapa 68905-230, Amapa, Brazil
5
University Hospitals Plymouth NHS Trust, Derriford Road, Crownhill, Plymouth PL6 8DH, UK
6
Vitalant Research Institute, 270 Masonic Avenue, San Francisco, CA 94118-4417, USA
7
Department Laboratory Medicine, University of California San Francisco, San Francisco, CA 94143, USA
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
‡
These authors jointly supervised this work.
Viruses 2020, 12(8), 866; https://doi.org/10.3390/v12080866
Submission received: 25 June 2020 / Revised: 25 July 2020 / Accepted: 3 August 2020 / Published: 8 August 2020
(This article belongs to the Special Issue Evolution and Diversity of Insect Viruses)
Browse Figures
Versions Notes

Abstract

:
Viruses are the most diverse and abundant microorganisms on earth, highly adaptive to a wide range of hosts. Viral diversity within invertebrate hosts has gained notoriety in recent years in public health as several such viruses have been of medical importance. Aedes aegypti serves as a vector for several viruses that have caused epidemics within the last year throughout Brazil; including Dengue, Zika and Chikungunya. This study aimed to identify new viral agents within Aedes aegypti mosquito in a city of the Amazonian region, where it is highly endemic. Metagenomic investigation was performed on 60 mosquito pools and viral RNA sequences present in their microbiota were characterized using genomic and phylogenetic tools. In total, we identified five putative novel virus species related to the Sobemovirus genus, Iflavirus genus and Permutatetraviridae family. These findings indicate a diverse taxonomy of viruses present in the mosquito microbiota of the Amazon, the region with the greatest invertebrate diversity in the world.

1. Introduction

Next-generation sequencing (NGS) is a revolutionary tool in molecular biology research. Courtesy of NGS, great numbers of insect microbiota have been explored, allowing the discovery of novel microorganisms, especially viruses [1,2,3]. Although the collective insect microbiota harbors many human pathogenic viruses, most viruses are non-pathogenic and have no direct public health impact. Despite this, some viruses are hypothesized to influence on mosquito susceptibility to certain arboviruses [4,5,6,7,8]. Furthermore, some phylogenetic studies indicate the evolution of pathogenic viruses from insect-specific virus to that capable of dual host tropism [1,9,10]. Therefore, metagenomic mosquito surveillance will enable new insights into the diversity and evolution of arboviruses.
The Amazon forest is considered one of the largest mosquito-related viral reservoirs in the world. Climatic and geographic conditions (frequent rainfall, year-round high temperature and dense forests) favors dissemination of several species of hematophagous diptera (mosquito, sandflies and ticks) and sylvatic animals [11]. Many mosquito-borne virus of medical importance are endemic throughout the Amazon, such as Dengue virus (DENV), Yellow Fever Virus, Oropouche virus and Mayaro virus; constituting significant morbidity for its population [12].
Ae. aegypti is an exotic mosquito species introduced to the Americas via slave ships coming from Africa [13]. Globally, the vector is closely associated with human habitations, demonstrating a great adaptive capacity to varied environments. The favorable climatic conditions of the Amazon have contributed to its domestication here [14]. Today, Ae. aegypti is responsible for transmitting a number of arboviruses such as dengue, Zika and chikungunya viruses, causing several epidemics in Brazil in the last few years [15].
Macapa in North Brazil is a typical Amazon city whereby it rains from December to May and has a dry summer from June to November. A recent study demonstrated high infestation rates of female Ae. Aegypti in Macapa [16]. It is well known that Dengue virus, Zika virus and Chikungunya virus are prevalent in this mosquito population. There is significant potential for the prevalence of other not-yet-characterized viruses within Ae. Aegypti of the Amazon and as such, may pose an emergent threat to public health. Previous studies from our group have described unusual viral sequences from Flaviridae and Reoviridae families in Ae. aegypti mosquitoes in the Amazon [17,18]. Using an NGS metagenomic approach, this study aimed to identify novel viruses in Ae. aegypti mosquito captured throughout an area of high-infestation from Macapa city, North Brazil. We found five potential novel viral species, suggestive of a rich taxonomy of arboviruses yet to be discovered in this Amazon region.

2. Materials and Methods

2.1. Mosquitoes Collection

Mosquitoes (Diptera: Culicidae) were collected from city of Macapá, Amapá state, North Brazil (see Figure S1), twice a month from January to March 2017. Electric manual aspirators and entomological nets were used to collect the mosquitoes. The mosquitoes were then transported to the laboratory, euthanized with ethyl acetate and morphologically identified using the dichotomous keys of Consoli and Lourenço-de-Oliveira [19] legs and wings removed. Between one and five females were grouped in pools according to their taxonomic category, place and date of collection. In total, 60 pools of mosquitoes were stored in a −80 °C freezer.

2.2. Sample Processing and Next Generation Sequencing (NGS)

The following metagenomics deep sequencing protocol was used. Initially, each mosquito pool was homogenized in 2 mL impact-resistant tube containing lysing matrix C (MP Biomedicals, USA) added to 900 μL of Hanks’ buffered salt solution (HBSS). The homogenized sample was centrifuged at 12,000× g for 10 min and approximately 300 μL of the supernatant was then filtrated through a 0.45 μm filter (Merck Millipore, Billerica, MA, USA). Next, 100 μL of cold PEG-it Virus Precipitation Solution (System Biosciences, CA, USA) was added to the obtained filtrate, mixed and incubated at 4 °C for 24 h. After, the mixture was centrifuged at 10,000× g for 30 min at 4 °C and supernatant discarded. The pellet rich in viral particles was treated with a mix of nuclease enzymes to digest unprotected nucleic acids. Viral nucleic acids were obtained using ZR & ZR-96 Viral DNA/RNA Kit (Zymo Research, CA, USA) according to the manufacturer’s protocol. The cDNA synthesis was conducted with AMV Reverse transcription (Promega, WI, USA). A second strand of cDNA synthesis was conducted using DNA Polymerase I Large Fragment (Promega, WI, USA). Then, DNA library was performed using Nextera XT Sample Preparation Kit (Illumina, CA, USA). The library was deep-sequenced using the HiSeq 2500 Sequencer (Illumina, CA, USA) with 126 bp ends. Bioinformatic analysis was performed according to the protocol previously described by Deng et al. [20]. The singlets and contigs were analyzed via BLAST (BLASTn and BLASTx) to look for similarity to viral sequence in GenBank’s Virus.

2.3. Phylogeny and Viral Annotation

Firstly, viral sequences identified in this study were used to query against NCBI protein database using the BLASTp tool to determine the closest sequences, its taxonomic classification and similarity. Secondly, based on BLAST result, the best hit sequences were download and aligned using Mafft software online [21] and phylogenetic trees were constructed using PhyML software [22] by Maximum Likelihood approach. Branch support values were assessed using the approximate likelihood ratio test (aLRT) on a Shimodaira-Hasegawa-like test. Evolutionary models and gamma distribution were selected according to the Bayesian information criterion (BIC) implemented in the jModeltest software [23]. Thirdly, the ORFs were annotated using InterProScan [24] and CD-search web using the CCD 3030 database and e-value < 0.05 [25].

3. Results

A total of 60 pools of Ae. aegypti female was collected (each pool containing between 1 and 5 specimens of mosquitoes, see the locations of sampling in the Figure S1 and characteristics of these pools in the Table S1), of which 24 were from Central and 36 were from Marabaixo and subsequently submitted to NGS protocol. Raw data were processed and after assembly the viral sequences were identified based on similarity of BLASTX comparison against to all RefSeq database in GenBank (details of viral richness in pools of pools contained viruses describe in this study were summarized in the Figures S4–S9, Supplementary Materials). We found nine virus-like sequences in two samples (AP59 and AP60) with <90% amino acid identity to different unclassified viruses, related to Sobemo-like virus, Iflavi-like virus and Permutotetra-like virus. These sequences represent five putative novel viruses, named Aedes Sobemo-like virus, Aedes Iflavi-like virus 1, Aedes Iflavi-like virus 2, Aedes permutotetra-like virus 1 and Aedes permutotetra-like virus 2 (Table 1). All sequences generated in this study were deposited in the GenBank (GenbBank acession: MT808014-MT808054) and they are also available in a fasta formatted file (Supplementary Materials; File S1).

3.1. Sobemo-Related Virus

Two Sobemo-like virus sequences (3000 and 2768 nt) were identified in two pools, tentatively named Aedes Sobemo-like virus (ASLV). Sobemo-like virus belongs to an unclassified group distantly related to the International Committee on Taxonomy of viruses (ICTV) Sobemovirus genus and Luteoviridae family. Sobemo-like viruses are widely found in insects and have a particular genomic organization; many of these novel viruses have (bi) segmented genomes, different to Sobemovirus genus and Luteoviridae family viruses which are monopartite [1]. We characterized two complete segment 1s of ASLV, which possesses two open read-frame (ORFs) corresponding to a putative peptidase and RNA-dependent RNA-polymarase (RdRp), respectively (Figure 1a). The ASLV RdRp gene encodes a 450-aa protein with amino acid identity varying between 72% and 83% with hypothetical protein 2 of Wenzhou Sobemo-like virus 4, the closest aligned sequence in the blastx, while the capsid shares 36–54% amino acid identity with same virus (Table 1). The ASLV has a typical overlapping reading frame, −1 frameshift and a protein layout similar to that of other known sobemoviruses (Figure 1a).
In addition, eighteen viral genomes with complete coding regions similar to Guadeloupe Mosquito virus (GMV) were obtained from several analyzed pools. The viral genomes contain two segments, encoding a putative peptidase and RdRp protein on segment 1 (2400–3000 nt) (Figure 1a) and a putative capsid and one hypothetical protein encoded by segment 2 (1000–1800 nt) (data not shown). Similarly to ASLV, GMV, Wenzhou sobemo-like virus 4 and Hubei mosquito virus 2 are all currently unclassified viruses with a distant relationship to the Luteoviridae family and Sobemovirus genus [26]. Phylogenetic analysis based on segment 1 indicates that GMV Brazilian sequences are highly similar to GMV, recently detected in Guadeloupe, and Renna virus isolated from Mexico City, sharing about 100% nucleotide sequence identity, clustered into a unique clade (Figure 1b).

3.2. Iflavi-Related Virus

One contig corresponding to the helicase gene and two contigs corresponding to the capsid and RdRp genes of two putative novel viruses; designated Aedes Iflavi-like virus 1 (AILV 1) and Aedes Iflavi-like virus 2 (AILV 2), were found in two mosquitoes pools (Figure 2a). All three fragments have low amino acid identity (<56%) with Yongsan picorna-like virus 1, the best hit with blastx (Table 1). This low identity is reflected in the topology of RdRp-based phylogeny, with AILV 1 and AILV 2 grouped in a separated clade from other Iflavi-like viruses (Figure 2b). The similar topologies were observed in the helicase and capsid-based ML phylogeny (Supplementary Materials). The phylogeny of best hits on blastx and Iflavirus members show diversity of Iflavi-like viruses, with all Iflavi-like and Iflavirus members previously isolated from arthropods, mostly insects [27]. AILV strains grouped into a cluster which shares a common ancestor with other viruses originally described in mosquitoes, dismembered from other insect-viruses (Figure 2b), however, only Yongsan picorna-like virus 1, AILV 1 and AILV 2 have been found in Aedes mosquitoes (unpublished). Since we used NGS to amplify viral sequences it is possible that contigs of ALV1 were derived from distinct viral genomes, likewise contigs of ALV2 also may have been amplified from distinct genomes. Nevertheless the phylogenetic analysis showed ALV1 and ALV2 are not the same virus and they likely represent new species.

3.3. Permutotetra-Like Virus

We found three partial genomic segments from two putative novel Permutotetra-like viruses (Figure 3). A RdRp sequence (3321 nt) presented 53% of amino acid identity with Culex Daeseongdong-like virus, the most similar virus. This putative novel virus was named Aedes permutotetra-like virus 1 (APLV1) (Figure 3a). Another two capsid sequences (882 and 1219 nt) belonging to this group shared ~48% amino acid identity with the most similar virus, Sarawak virus. This putative novel virus was named Aedes permutotetra-like virus 2 (APLV2) (Figure 3b). The cluster formed by APLV-1, Culex Daeseongdong-like, Daeseongdong virus 2 and Smothfield permutotetra-like virus have been found in mosquitoes (Figure 3c) [3,28]; similarly, the clade formed by APLV-2, Culex permutotetra virus, Shinobi tetravirus and Sarawak virus have also been detected in different mosquito species (Figure 3d) [29,30,31].

4. Discussion

In this study, we analyzed 60 pools of Ae. aegypti, and, as expected, we found several highly divergent sequences, which possibly represent novel viral species. These viruses belong to the Luteo-sobemo-related virus, Iflavirus and Alphapermutotetravirus genus.
One novel Luteo-sobemo related virus was found in two Ae. aegypti samples, named Aedes sobemo-like virus (ASLV). The sobemo-like viruses are (+) ssRNA unclassified viruses distantly related to the Sobemovirus genus and Luteoviridae family (Figure 1b), which infects plants and is known to be vectored by arthropods. Although viruses belonging to Sobemovirus genus and Luteoviridae family are known plant viruses and are of monopartite genome, Sobemo-like virus members have bi-segmented genomes and have been isolated primarily from insects [1,26]. It has been speculated that this group of viruses should be proposed as a new family [32].
Iflavirus members are a new recognized family called Iflaviridae (order Picornavirales), under the Iflavirus genus. All Iflavirus members are insect-infecting viruses and they have been identified in a wide range of hosts belonging to the class Insecta, although plant-infecting Iflavirus-like virus has been reported from tomato (Solanum lycopersicum) [33]. Currently, there are fifteen species in the genus Iflavirus recognized in the last report of ICTV [34]. However, sequence identity at the amino acid level of the capsid proteins above 90% is used for species demarcation criteria for the Iflavirus genus and several tentative novel viruses have been identified showing sequence similarity to members of the genus Iflavirus and yet are classified as iflavi-like viruses. Through NGS analysis, we assembled five contigs that showed similarity to Iflavi-like viruses, here named Aedes Iflavi-like virus 1 (AILV1) and Aedes Iflavi-like virus 2 (AILV2). For this reason, there is a possibility of these contigs belong to distinct viral genomes. Our analysis shows that ALV1 and ALV2 share a common ancestor that diverge from each other and likely represent two new viral species. BLASTp searches showed that both AILV1 and AILV2 shared low sequence identity (less than 90%) with other members of Iflaviridae at the amino acid level (Table 1), indicating that both are novel species of Iflaviridae family [34]. Additionally, sequence analysis showed that the AILV1 and AILV2 shared 50% capsid amino acid sequence identity (data not shown), suggesting that they are members of the different species. The phylogenetic analysis showed that AILV1 and AILV2 form a well-supported clade, suggesting the representation of a novel clade within the Iflaviridae family. According to ICTV, “The Iflaviridae family is expanding rapidly and will likely undergo revision in the near future” and possibly new species and genus will be included in the official taxonomy.
Additionally, we also detected three partial genomes of two putative novel viruses (APLV 1 and APLV2) closely related to unclassified permutotetra-like viruses. Permutotetraviridae is a recent classified family with a single genus (Alphapermutotetravirus) and two prototype species (Euprosterna elaeasa virus and Thosea asigna virus), restricted mainly to insects in the order Lepidoptera (butterflies and moths). In recent years, a wide range of highly divergent viruses distantly related to Permutotetraviridae family has been identified [1,3,30,35]. The lack of common genomic organizations in permutotetra-like virus members and the formation of two large and well-supported clades (Figure 3) support the need to create new groups for the current unclassified viruses of this family. So, the permutotetra-like viruses (APLV1 and APLV2) found in this study may represent new species within different genus/family. Both viruses are grouped with viruses isolated only from mosquitoes, indicating a likely common origin within their respective clades (Figure 3c,d).
All novel viruses reported here share a common ancestor with other viruses originally described in mosquitoes, dismembered from other insect-viruses, suggesting a close evolution with their mosquito hosts. Recent phylogenetic studies in several RNA insect-virus families have indicated that they are ancient agents with highly distinct lineages, leading to the credence of probable co-evolution and expansion with their insect hosts [10,36,37]. The hypothesis that insect-viruses have been closely associated with their insect hosts for a long period of time is supported by studies that demonstrate vertical/transovarial transmission (TOT), whereas some become integrated into the genomes of their own arthropod hosts [38,39,40]. Another possibility regarding insect-virus evolution is of host association, whereby dual host viruses evolved from insect-specific progenitors, with many arthropod-borne viruses possibly emerging to vertebrates and plants in this way, including complete adaptation to vertebrate or plant hosts and thereby losing the need for an invertebrate host [2,41]. Importantly, none of these novel viruses are closely related to known vector-borne pathogens of humans or other mammals. The alignments used to construct maximum likelihood trees (ML) have a poor phylogenetic signal, as determined by the high proportion of star-like trees in the likelihood mapping analysis. Although the low quality of alignments has little effect on the topology of trees constructed with other methods, it may have a significant outcome to estimates mutation rates or to the measurement of divergent-time. Consequently, in the future with the identification of new viral sequences the interpretation of data may also be affected.
The diversity of arboviruses remains to be explored, especially in the Amazon, known for being a rich bioma with many viral species. For example, the Amazon region is identified as the starting point for transmission of the yellow fever virus in a recent outbreak in Brazil that killed almost 700 people between December 2016 and March 2018 and the Amazon rainforest functions as a ‘reservoir’ region for several arboviruses [42]. Other studies have been identified novel virus presents in the mosquito microbiota in this region, in the host Aedes aegypti and anophelines [17,18,43], reinforcing the idea that our current knowledge about the diversity of viruses is still very limited. Furthermore, insect-specific virus (ISV) compose the majority of mosquito virome the virus-virus interactions may affect the transmission of some viral pathogens [44,45]. From description and characterization of these viral agents, we can gain knowledge useful to some biotechnological strategies in combating epidemic viruses, such as Dengue, Zika and Chikungunya.
Recent massive metagenomics studies have expanded our knowledge about diversity of a great number of invertebrate viruses, which include many unclassified groups, inclusive of luteo-sobemo-like virus, Ifla-like virus, permutotetra-like virus, among others [1,2,3,46]. The ICTV officialization of these taxonomic proposals requires time, alongside large quantities of sequence analysis and therefore, the identification of novel viral sequence within this study greatly contributes to the correct taxonomic classification of these “virus-like” sequences. Ultimately, our results highlight the importance of identifying and characterizing novel viruses to expand our understanding of the taxonomic diversity of viral groups (families, genera and species), which is currently poor. The host range influence and vector biology for these viruses, alongside the ecological and evolutionary history of Aedes aegypti microbiota, the principle arboviral vectors of Brazil, need to be further studied.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/12/8/866/s1, File 1: Text containing all contigs in fasta format; File 2: Figure S1. Map showing the coordinates where mosquitoes were captures in Brazil; Figure S2. Phylogenetic tree of Iflavi-like virus and Iflavirus genus members based on the alignment of capsid region protein; Figure S3. Phylogenetic tree of Iflavi-like virus and Iflavirus genus members based on the alignment of helicase region protein. Figures S4–S8. Virome composition of pools of Aedes mosquitoes. Figure S9. Comparison of viral richness of pool AP59 of Aedes mosquitoes. Table S1. Characteristics of pool samples of Aedes mosquitoes from Macapa.

Author Contributions

Specimen Collection or Sample Screening, V.S.M.; F.J.C.M.; E.S.D.R.; M.O.d.S.R. and R.N.P.S.; Conceptualization, E.C.S.; S.D.; E.D. and A.C.d.C.; sequencing and contig analysis, A.C.d.C.; L.N.F.; X.D. and E.D.; Phylogenetic analysis: G.d.O.R.; F.V. and É.L.; Data Investigation: A.C.d.C.; É.L. and E.C.S.; Writing original draft: G.d.O.R. and É.L.; Manuscript editing: R.T. and P.M.H. Supervision, E.C.S.; A.C.d.C.; E.D. and É.L.; All authors have read and agreed to the published version of the manuscript. All authors critically revised the manuscript for intellectual content and approved the final version.

Funding

This study was supported by FAPESP #2016/01735-2, CNPq #400354/2016-0, #408054/2016-6 and FAPEAP/CNPq (PPSUS 250.203.0333/2016). Antonio Charlys da Costa is funded by FAPESP #2017/00021-9 and Vanessa S. Morais is supported by a scholarship provided by Fapesp (#2019/21706-5). Élcio Leal is supported by CNPq #302677/2019-4. Geovani de Oliveira Ribeiro is supported by a scholarship provided by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES).

Acknowledgments

We thank Luciano Monteiro da Silva for the administrative support. Our thanks also to the Sage Science Inc., MP Biomedicals Inc. and Zymo Research Corporation, for the donation of reagents. We also thank Pró-reitoria de pesquisa e pós-graduação of UFPA for supporting the publication costs.

Conflicts of Interest

The authors declare no competing interests.

References

  1. Shi, M.; Lin, X.D.; Tian, J.H.; Chen, L.J.; Chen, X.; Li, C.X.; Qin, X.C.; Li, J.; Cao, J.P.; Eden, J.S.; et al. Redefining the invertebrate RNA virosphere. Nature 2016, 540, 539–543. [Google Scholar] [CrossRef]
  2. Li, C.X.; Shi, M.; Tian, J.H.; Lin, X.D.; Kang, Y.J.; Chen, L.J.; Qin, X.C.; Xu, J.; Holmes, E.C.; Zhang, Y.Z. Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses. eLife 2015, 4, e05378. [Google Scholar] [CrossRef]
  3. Sadeghi, M.; Altan, E.; Deng, X.; Barker, C.M.; Fang, Y.; Coffey, L.L.; Delwart, E. Virome of > 12 thousand Culex mosquitoes from throughout California. Virology 2018, 523, 74–88. [Google Scholar] [CrossRef]
  4. Kenney, J.L.; Solberg, O.D.; Langevin, S.A.; Brault, A.C. Characterization of a novel insect-specific flavivirus from Brazil: Potential for inhibition of infection of arthropod cells with medically important flaviviruses. J. Gen. Virol. 2014, 95, 2796. [Google Scholar] [CrossRef]
  5. Nasar, F.; Erasmus, J.H.; Haddow, A.D.; Tesh, R.B.; Weaver, S.C. Eilat virus induces both homologous and heterologous interference. Virology 2015, 484, 51–58. [Google Scholar] [CrossRef] [Green Version]
  6. Öhlund, P.; Lundén, H.; Blomström, A.L. Insect-specific virus evolution and potential effects on vector competence. Virus Genes 2019, 55, 127–137. [Google Scholar]
  7. Romo, H.; Kenney, J.L.; Blitvich, B.J.; Brault, A.C. Restriction of Zika virus infection and transmission in Aedes aegypti mediated by an insect-specific flavivirus. Emerg. Microbes Infect. 2018, 7, 1–13. [Google Scholar] [CrossRef] [Green Version]
  8. Vasilakis, N.; Tesh, R.B. Insect-specific viruses and their potential impact on arbovirus transmission. Curr. Opin. Virol. 2015, 15, 69–74. [Google Scholar] [CrossRef] [Green Version]
  9. Halbach, R.; Junglen, S.; van Rij, R.P. Mosquito-specific and mosquito-borne viruses: Evolution, infection, and host defense. Curr. Opin. Insect Sci. 2017, 22, 16–27. [Google Scholar] [CrossRef]
  10. Marklewitz, M.; Zirkel, F.; Kurth, A.; Drosten, C.; Junglen, S. Evolutionary and phenotypic analysis of live virus isolates suggests arthropod origin of a pathogenic RNA virus family. Proc. Natl. Acad. Sci. USA 2015, 112, 7536–7541. [Google Scholar] [CrossRef] [Green Version]
  11. Vasconcelos, P.F.; Travassos da Rosa, A.P.; Rodrigues, S.G.; Travassos da Rosa, E.S.; Dégallier, N.; Travassos da Rosa, J.F. Inadequate management of natural ecosystem in the Brazilian Amazon region results in the emergence and reemergence of arboviruses. Cad. Saúde Pública/Ministério da Saúde Fundação Oswaldo Cruz Esc. Nac. Saúde Pública 2001, 17, 155–164. [Google Scholar] [CrossRef]
  12. Donalisio, M.R.; Freitas, A.R.R.; Zuben, A.P.B. Von Arboviruses emerging in Brazil: Challenges for clinic and implications for public health. Rev. Saude Publica 2017, 51, 30. [Google Scholar] [CrossRef]
  13. Brown, J.E.; Evans, B.R.; Zheng, W.; Obas, V.; Barrera-Martinez, L.; Egizi, A.; Zhao, H.; Caccone, A.; Powell, J.R. Human impacts have shaped historical and recent evolution in Aedes aegypti, the dengue and yellow fever mosquito. Evolution (N. Y.) 2014, 68, 514–525. [Google Scholar] [CrossRef] [Green Version]
  14. Guagliardo, S.A.J.; Lee, Y.; Pierce, A.A.; Wong, J.; Chu, Y.Y.; Morrison, A.C.; Astete, H.; Brosi, B.; Vazquez-Prokopec, G.; Scott, T.W.; et al. The genetic structure of Aedes aegypti populations is driven by boat traffic in the Peruvian Amazon. PLoS Negl. Trop. Dis. 2019, 13. [Google Scholar] [CrossRef] [Green Version]
  15. Paixão, E.S.; Teixeira, M.G.; Rodrigues, L.C. Zika, chikungunya and dengue: The causes and threats of new and re-emerging arboviral diseases. BMJ Glob. Health 2017, 3, e000530. [Google Scholar] [CrossRef] [Green Version]
  16. Monteiro, F.J.C.; Mourão, F.R.P.; Ribeiro, E.S.D.; da Rêgo, M.O.S.; da Frances, P.A.C.; Souto, R.N.P.; Façanha, M.D.S.; Tahmasebi, R.; da Costa, A.C. Prevalence of dengue, zika and chikungunya viruses in aedes (Stegomyia) aegypti (diptera: Culicidae) in a medium-sized city, Amazon, Brazil. Rev. Inst. Med. Trop. Sao Paulo 2020, 62. [Google Scholar] [CrossRef] [Green Version]
  17. Fernandes, L.N.; Coletti, T.M.; Monteiro, F.J.C.; Rego, M.O.S.; Ribeiro, E.S.D.; Ribeiro, G.O.; Marinho, R.S.S.; Komninakis, S.V.; Witkin, S.S.; Deng, X.; et al. A novel highly divergent strain of cell fusing agent virus (Cfav) in mosquitoes from the Brazilian Amazon region. Viruses 2018, 10, 666. [Google Scholar] [CrossRef] [Green Version]
  18. de Ribeiro, G.O.; Monteiro, F.J.C.; da Rego, M.O.S.; Ribeiro, E.S.D.; de Castro, D.F.; Caseiro, M.M.; Marinho, S.; dos Santos, R.; Komninakis, S.V.; Witkin, S.S. Detection of RNA-Dependent RNA Polymerase of Hubei Reo-Like Virus 7 by Next-Generation Sequencing in Aedes aegypti and Culex quinquefasciatus Mosquitoes from Brazil. Viruses 2019, 11, 147. [Google Scholar] [CrossRef] [Green Version]
  19. Consoli, R.A.G.B.; de Oliveira, R.L. Principais Mosquitos de Importância Sanitária No Brasil, 1st ed.; Editora Fiocruz: Rio de Janeiro, Brasil, 1994; ISBN 85-85676-03-5. [Google Scholar]
  20. Deng, X.; Naccache, S.N.; Ng, T.; Federman, S.; Li, L.; Chiu, C.Y.; Delwart, E.L. An ensemble strategy that significantly improves de novo assembly of microbial genomes from metagenomic next-generation sequencing data. Nucleic Acids Res. 2015, 43, e46. [Google Scholar] [CrossRef]
  21. Katoh, K.; Rozewicki, J.; Yamada, K.D. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. Brief. Bioinform. 2017. [Google Scholar] [CrossRef] [Green Version]
  22. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Posada, D. jModelTest: Phylogenetic model averaging. Mol. Biol. Evol. 2008, 25, 1253–1256. [Google Scholar] [CrossRef]
  24. Quevillon, E.; Silventoinen, V.; Pillai, S.; Harte, N.; Mulder, N.; Apweiler, R.; Lopez, R. InterProScan: Protein domains identifier. Nucleic Acids Res. 2005, 33, W116–W120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Marchler-Bauer, A.; Derbyshire, M.K.; Gonzales, N.R.; Lu, S.; Chitsaz, F.; Geer, L.Y.; Geer, R.C.; He, J.; Gwadz, M.; Hurwitz, D.I.; et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015, 43, 222–226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Shi, C.; Beller, L.; Deboutte, W.; Yinda, K.C.; Delang, L.; Vega-Rúa, A.; Failloux, A.B.; Matthijnssens, J. Stable distinct core eukaryotic viromes in different mosquito species from Guadeloupe, using single mosquito viral metagenomics. Microbiome 2019, 7, 121. [Google Scholar] [CrossRef]
  27. Oers, M.M. Genomics and biology of lflaviruses. Insect Virol. 2010, 231–250. [Google Scholar]
  28. Hang, J.; Klein, T.A.; Kim, H.-C.; Yang, Y.; Jima, D.D.; Richardson, J.H.; Jarman, R.G. Genome sequences of five arboviruses in field-captured mosquitoes in a unique rural environment of South Korea. Genome Announc. 2016, 4. [Google Scholar] [CrossRef] [Green Version]
  29. Amoa-Bosompem, M.; Kobayashi, D.; Murota, K.; Faizah, A.N.; Itokawa, K.; Fujita, R.; Osei, J.H.N.; Agbosu, E.; Pratt, D.; Kimura, S. Entomological Assessment of the Status and Risk of Mosquito-borne Arboviral Transmission in Ghana. Viruses 2020, 12, 147. [Google Scholar] [CrossRef] [Green Version]
  30. Fujita, R.; Kato, F.; Kobayashi, D.; Murota, K.; Takasaki, T.; Tajima, S.; Lim, C.K.; Saijo, M.; Isawa, H.; Sawabe, K. Persistent viruses in mosquito cultured cell line suppress multiplication of flaviviruses. Heliyon 2018, 4, e00736. [Google Scholar] [CrossRef] [Green Version]
  31. Sadeghi, M.; Popov, V.; Guzman, H.; Phan, T.G.; Vasilakis, N.; Tesh, R.; Delwart, E. Genomes of viral isolates derived from different mosquitos species. Virus Res. 2017, 242, 49–57. [Google Scholar] [CrossRef] [Green Version]
  32. Shi, C.; Zhao, L.; Zeng, W.; Atoni, E.; Hu, X.; Matthijnssens, J.; Yuan, Z.; Xia, H. The conservation of a core virome in Aedes mosquitoes across different developmental stages and continents. bioRxiv 2020. [Google Scholar] [CrossRef]
  33. Saqib, M.; Wylie, S.J.; Jones, M.G.K. Serendipitous identification of a new Iflavirus-like virus infecting tomato and its subsequent characterization. Plant Pathol. 2015, 64, 519–527. [Google Scholar] [CrossRef] [Green Version]
  34. Valles, S.M.; Chen, Y.; Firth, A.E.; Guérin, D.M.A.; Hashimoto, Y.; Herrero, S.; De Miranda, J.R.; Ryabov, E. ICTV virus taxonomy profile: Iflaviridae. J. Gen. Virol. 2017, 98, 527–528. [Google Scholar] [CrossRef] [PubMed]
  35. Dalmon, A.; Gayral, P.; Decante, D.; Klopp, C.; Bigot, D.; Thomasson, M.; AHerniou, E.; Alaux, C.; Conte, Y. Le Viruses in the invasive hornet vespa velutina. Viruses 2019, 11, 1041. [Google Scholar] [CrossRef] [Green Version]
  36. Dudas, G.; Obbard, D.J. Are arthropods at the heart of virus evolution? eLife 2015, 4, e06837. [Google Scholar] [CrossRef] [Green Version]
  37. Cook, S.; Chung, B.Y.W.; Bass, D.; Moureau, G.; Tang, S.; McAlister, E.; Culverwell, C.L.; Glücksman, E.; Wang, H.; Brown, T.D.K.; et al. Novel virus discovery and genome reconstruction from field rna samples reveals highly divergent viruses in dipteran hosts. PLoS ONE 2013, 8, e80720. [Google Scholar] [CrossRef] [Green Version]
  38. Fort, P.; Albertini, A.; Van-Hua, A.; Berthomieu, A.; Roche, S.; Delsuc, F.; Pasteur, N.; Capy, P.; Gaudin, Y.; Weill, M. Fossil rhabdoviral sequences integrated into arthropod genomes: Ontogeny, evolution, and potential functionality. Mol. Biol. Evol. 2012, 29, 381–390. [Google Scholar] [CrossRef]
  39. Bolling, B.G.; Olea-Popelka, F.J.; Eisen, L.; Moore, C.G.; Blair, C.D. Transmission dynamics of an insect-specific flavivirus in a naturally infected Culex pipiens laboratory colony and effects of co-infection on vector competence for West Nile virus. Virology 2012, 427, 90–97. [Google Scholar] [CrossRef] [Green Version]
  40. Lequime, S.; Lambrechts, L. Discovery of flavivirus-derived endogenous viral elements in Anopheles mosquito genomes supports the existence of Anopheles-associated insect-specific flaviviruses. Virus Evol. 2017, 3. [Google Scholar] [CrossRef] [Green Version]
  41. Kuno, G.; Chang, G.-J.J. Biological transmission of arboviruses: Reexamination of and new insights into components, mechanisms, and unique traits as well as their evolutionary trends. Clin. Microbiol. Rev. 2005, 18, 608–637. [Google Scholar] [CrossRef] [Green Version]
  42. Faria, N.R.; Kraemer, M.U.G.; Hill, S.C.; De Jesus, J.G.; Aguiar, R.S.; Iani, F.C.M.; Xavier, J.; Quick, J.; Du Plessis, L.; Dellicour, S.; et al. Genomic and epidemiological monitoring of yellow fever virus transmission potential. Science 2018, 361, 894–899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Scarpassa, V.M.; Debat, H.J.; Alencar, R.B.; Saraiva, J.F.; Calvo, E.; Arcà, B.; Ribeiro, J.M.C. An insight into the sialotranscriptome and virome of Amazonian anophelines. BMC Genom. 2019, 20, 166. [Google Scholar] [CrossRef] [PubMed]
  44. Blitvich, B.J.; Firth, A.E. Insect-specific flaviviruses: A systematic review of their discovery, host range, mode of transmission, superinfection exclusion potential and genomic organization. Viruses 2015, 7, 1927–1959. [Google Scholar] [CrossRef] [Green Version]
  45. Kuwata, R.; Isawa, H.; Hoshino, K.; Sasaki, T.; Kobayashi, M.; Maeda, K.; Sawabe, K. Analysis of mosquito-borne Flavivirus superinfection in Culex tritaeniorhynchus (Diptera: Culicidae) cells persistently infected with Culex Flavivirus (Flaviviridae). J. Med. Entomol. 2015, 52, 222–229. [Google Scholar] [CrossRef] [PubMed]
  46. Käfer, S.; Paraskevopoulou, S.; Zirkel, F.; Wieseke, N.; Donath, A.; Petersen, M.; Jones, T.C.; Liu, S.; Zhou, X.; Middendorf, M.; et al. Re-assessing the diversity of negative strand RNA viruses in insects. PLoS Pathog. 2019, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Viruses 12 00866 g001 550
Figure 1. Luteo-sobemo-related viruses genomes map and phylogeny. (a) Open read-frame organization of closest related virus to Aedes sobemo-like virus and Guadeloupe mosquito virus. (b) Nucleotide maximum likelihood tree based on segment 1 for Sobemo-like viruses and Luteovirus and Sobemovirus genus. Aedes sobemo-like virus and Guadeloupe mosquito virus sequenced in this study are highlighted in orange and green, respectively. The diagram in the base of the tree is the likelihood map of the nucleotide alignment of genomes of Luteo-sobemo-related viruses. The likelihood quartet mapping is a method that allows to visualize the tree-likeness of all quartets in a single graph and provide a direct measure of the phylogenetic signal in an alignment. The triangle shows the location of all quartets calculated with the alignment used to infer the ML tree. Values in the center of the triangle represent the percentage of unresolved quartet trees (star-like trees), values in the vertices represent the percentage of fully resolved trees and values in the intermediate areas (between vertices) are the percentage of conflicting trees. The analysis was performed using GTR+ gamma correction model as is implemented in the tree puzzle software v 5.3 (http://www.tree-puzzle.de).
Figure 1. Luteo-sobemo-related viruses genomes map and phylogeny. (a) Open read-frame organization of closest related virus to Aedes sobemo-like virus and Guadeloupe mosquito virus. (b) Nucleotide maximum likelihood tree based on segment 1 for Sobemo-like viruses and Luteovirus and Sobemovirus genus. Aedes sobemo-like virus and Guadeloupe mosquito virus sequenced in this study are highlighted in orange and green, respectively. The diagram in the base of the tree is the likelihood map of the nucleotide alignment of genomes of Luteo-sobemo-related viruses. The likelihood quartet mapping is a method that allows to visualize the tree-likeness of all quartets in a single graph and provide a direct measure of the phylogenetic signal in an alignment. The triangle shows the location of all quartets calculated with the alignment used to infer the ML tree. Values in the center of the triangle represent the percentage of unresolved quartet trees (star-like trees), values in the vertices represent the percentage of fully resolved trees and values in the intermediate areas (between vertices) are the percentage of conflicting trees. The analysis was performed using GTR+ gamma correction model as is implemented in the tree puzzle software v 5.3 (http://www.tree-puzzle.de).
Viruses 12 00866 g001
Viruses 12 00866 g002 550
Figure 2. Iflavi-related contigs position in Yongsan picorna-like virus genome. (a) Contigs position of Aedes Iflavi-like virus and structure-based alignment with relative closest virus, Yongsan picorna-like virus. (b) Maximum likelihood tree (ML) for RdRp protein of Aedes Iflavi-like virus (orange) with Iflavirus genus (indicated in a blue area) and Iflavi-like sequences related to Aedes Iflavi-like virus by Blastp search. Sequences from the current study are colored in red. The diagram in the base of the tree is the likelihood map of the nucleotide alignment of genomes of Iflavi-related viruses. The likelihood quartet mapping is a method that allows to visualize the tree-likeness of all quartets in a single graph and provide a direct measure of the phylogenetic signal in an alignment. The triangle shows the location of all quartets calculated with the alignment used to infer the ML tree. Values in the center of the triangle represent the percentage of unresolved quartet trees (star-like trees), values in the vertices represent the percentage of fully resolved trees and values in the intermediate areas (between vertices) are the percentage of conflicting trees. The analysis was performed using JTT model as is implemented in the tree puzzle software v 5.3 (http://www.tree-puzzle.de).
Figure 2. Iflavi-related contigs position in Yongsan picorna-like virus genome. (a) Contigs position of Aedes Iflavi-like virus and structure-based alignment with relative closest virus, Yongsan picorna-like virus. (b) Maximum likelihood tree (ML) for RdRp protein of Aedes Iflavi-like virus (orange) with Iflavirus genus (indicated in a blue area) and Iflavi-like sequences related to Aedes Iflavi-like virus by Blastp search. Sequences from the current study are colored in red. The diagram in the base of the tree is the likelihood map of the nucleotide alignment of genomes of Iflavi-related viruses. The likelihood quartet mapping is a method that allows to visualize the tree-likeness of all quartets in a single graph and provide a direct measure of the phylogenetic signal in an alignment. The triangle shows the location of all quartets calculated with the alignment used to infer the ML tree. Values in the center of the triangle represent the percentage of unresolved quartet trees (star-like trees), values in the vertices represent the percentage of fully resolved trees and values in the intermediate areas (between vertices) are the percentage of conflicting trees. The analysis was performed using JTT model as is implemented in the tree puzzle software v 5.3 (http://www.tree-puzzle.de).
Viruses 12 00866 g002
Viruses 12 00866 g003 550
Figure 3. Permutotetra-related contigs position in most related viruses genomes. (a) Schematic representation of Culex Daeseongdong-like virus, Aedes permutetra-like virus 1, (b) Sarawak virus and Aedes permutetra-like virus 2. (c) Maximum likelihood phylogenetic tree for RdRp protein and (d) maximum likelihood phylogenetic tree for capsid protein. Viruses originally found in mosquitoes are marked by orange. Sequences from this study are indicated in red. * Indicate the location in the RdRp tree of the reference viruses used to construct the capsid tree. The diagram in the base of the tree is the likelihood map. The likelihood quartet mapping is a method that allows to visualize the tree-likeness of all quartets in a single graph and provide a direct measure of the phylogenetic signal in an alignment. The triangle shows the location of all quartets calculated with the alignment used to infer the ML tree. Values in the center of the triangle represent the percentage of unresolved quartet trees (star-like trees), values in the vertices represent the percentage of fully resolved trees and values in the intermediate areas (between vertices) are the percentage of conflicting trees. The analysis was performed using JTT+gamma correction model as is implemented in the tree puzzle software v 5.3 (http://www.tree-puzzle.de).
Figure 3. Permutotetra-related contigs position in most related viruses genomes. (a) Schematic representation of Culex Daeseongdong-like virus, Aedes permutetra-like virus 1, (b) Sarawak virus and Aedes permutetra-like virus 2. (c) Maximum likelihood phylogenetic tree for RdRp protein and (d) maximum likelihood phylogenetic tree for capsid protein. Viruses originally found in mosquitoes are marked by orange. Sequences from this study are indicated in red. * Indicate the location in the RdRp tree of the reference viruses used to construct the capsid tree. The diagram in the base of the tree is the likelihood map. The likelihood quartet mapping is a method that allows to visualize the tree-likeness of all quartets in a single graph and provide a direct measure of the phylogenetic signal in an alignment. The triangle shows the location of all quartets calculated with the alignment used to infer the ML tree. Values in the center of the triangle represent the percentage of unresolved quartet trees (star-like trees), values in the vertices represent the percentage of fully resolved trees and values in the intermediate areas (between vertices) are the percentage of conflicting trees. The analysis was performed using JTT+gamma correction model as is implemented in the tree puzzle software v 5.3 (http://www.tree-puzzle.de).
Viruses 12 00866 g003
Table
Table 1. Amino acid similarity of sequences identified in metagenomics analysis from Ae. aegypti pools.
Table 1. Amino acid similarity of sequences identified in metagenomics analysis from Ae. aegypti pools.
Virus NameClosely Related Viruses 1GeneLength (nt)Cover 1Amino Acid Identity 1
Aedes Sobemo-like virus strain AP60-1 *Wenzhou Sobemo-like virus 4Peptidase171942%36%
RdRp130888%72%
Aedes Sobemo-like virus strain AP60-2 *Wenzhou Sobemo-like virus 4Peptidase276899%54%
RdRp1308100%83%
Aedes Iflavi-like virus 1Yongsan picorna-like virus 1Capsid256160%49%
Helicase1464100%46%
RdRp201697%48%
Aedes Iflavi-like virus 2Yongsan picorna-like virus 1Capsid196549%47%
RdRp128550%52%
Aedes permutotetra-like virus 1Culex Daeseongdong-like virusRdRp332192%53%
Aedes permutotetra-like virus 2 strain AP59 *Sarawak virusCapsid121991%48%
Aedes permutotetra-like virus 2 strain AP60 *Sarawak virusCapsid88294%49%
1 Based on Blastp. * Same virus (>90% amino acid identity between them).

Share and Cite

MDPI and ACS Style

Ribeiro, G.d.O.; Morais, V.S.; Monteiro, F.J.C.; Ribeiro, E.S.D.; Rego, M.O.d.S.; Souto, R.N.P.; Villanova, F.; Tahmasebi, R.; Hefford, P.M.; Deng, X.; et al. Aedes aegypti from Amazon Basin Harbor High Diversity of Novel Viral Species. Viruses 2020, 12, 866. https://doi.org/10.3390/v12080866

AMA Style

Ribeiro GdO, Morais VS, Monteiro FJC, Ribeiro ESD, Rego MOdS, Souto RNP, Villanova F, Tahmasebi R, Hefford PM, Deng X, et al. Aedes aegypti from Amazon Basin Harbor High Diversity of Novel Viral Species. Viruses. 2020; 12(8):866. https://doi.org/10.3390/v12080866

Chicago/Turabian Style

Ribeiro, Geovani de Oliveira, Vanessa S. Morais, Fred Julio Costa Monteiro, Edcelha Soares D’Athaide Ribeiro, Marlisson Octavio da S Rego, Raimundo Nonato Picanço Souto, Fabiola Villanova, Roozbeh Tahmasebi, Philip Michael Hefford, Xutao Deng, and et al. 2020. "Aedes aegypti from Amazon Basin Harbor High Diversity of Novel Viral Species" Viruses 12, no. 8: 866. https://doi.org/10.3390/v12080866

APA Style

Ribeiro, G. d. O., Morais, V. S., Monteiro, F. J. C., Ribeiro, E. S. D., Rego, M. O. d. S., Souto, R. N. P., Villanova, F., Tahmasebi, R., Hefford, P. M., Deng, X., Delwart, E., Cerdeira Sabino, E., Fernandes, L. N., da Costa, A. C., & Leal, É. (2020). Aedes aegypti from Amazon Basin Harbor High Diversity of Novel Viral Species. Viruses, 12(8), 866. https://doi.org/10.3390/v12080866

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop