Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 25 April 2022
Sec. Environmental Health and Exposome
This article is part of the Research Topic Tackling the Presence of Pathogens in Water and Wastewater View all 6 articles

Viral Metagenomics Reveals Widely Diverse Viral Community of Freshwater Amazonian Lake

  • 1Department of Biochemistry and Immunology, Federal University of Minas Gerais, Belo Horizonte, Brazil
  • 2Institute of Biological Sciences, Federal University of Pará, Belém, Brazil
  • 3Laboratory of Genomic and Bioinformatics, Center of Genomics and System Biology, Federal University of Pará, Belém, Brazil

Despite the importance of understanding the ecology of freshwater viruses, there are not many studies on the subject compared to marine viruses. The microbiological interactions in these environments are still poorly known, especially between bacteriophages and their host bacteria and between cyanophages and cyanobacteria. Lake Bologna, Belém, capital of the Brazilian State of Pará, is a water source that supplies the city and its metropolitan region. However, it remains unexplored regarding the contents of its virome and viral diversity composition. Therefore, this work aims to explore the taxonomic diversity of DNA viruses in this lake, especially bacteriophages and cyanophages, since they can act as transducers of resistance genes and reporters of water quality for human consumption. We used metagenomic sequencing data generated by previous studies. We analyzed it at the taxonomic level using the tools Kraken2, Bracken, and Pavian; later, the data was assembled using Genome Detective, which performs the assembly of viruses. The results observed here suggest the existence of a widely diverse viral community and established microbial phage-regulated dynamics in Lake Bolonha. This work is the first ever to describe the virome of Lake Bolonha using a metagenomic approach based on high-throughput sequencing, as it contributes to the understanding of water-related public health concerns regarding the spreading of antibiotic resistance genes and population control of native bacteria and cyanobacteria.

Introduction

Amazonia is a 10 million years old unrivaled nest of biodiversity that reigns over South America; from bird-eating spiders to emperor tamarins to pink river dolphins, biologists find a new species every other day (1). Its ecosystems are essential for biodiversity preservation, climate regulation, energy production, and food and water security. The Amazon has a vital role in controlling zoonotic diseases and vector-borne and water-borne infections (2, 3).

Despite its importance, policies, laws, agreements, funds, and practical actions focused on Amazon protection have been weakened in Brazil, encouraging deforestation and culminating in losing about 20% of the original Amazon forest cover in Brazil by 2019 (4). The association between anthropogenic action in the Amazon rainforest, eutrophication of its water bodies, climate change, and alterations in vector dynamics, human migration, genetic changes in pathogens, and the poor social and environmental conditions in many Latin-American countries serve as an opportunity for the emergence and re-emergence of human infectious diseases in Brazil and other Amazonian countries (5).

Amazonian fauna hosts a vast diversity of well-known pathogens and many other potential new or even unknown pathogens (611). This abundance of microorganisms indicates that the emergence of new infections from the forest is a constant threat to human health, particularly favorable to water-borne diseases due to anthropogenic activities nearby freshwater sources, such as ponds, rivers, basins, and lakes (1214).

Lakes close to urban areas are increasingly changing their ecosystem as human population expansion occurs and commercial, recreational, and residential uses increase (15). The eutrophication process offers particular conditions for the replication of viruses, as environments of this type seem to provide high viral activation and hypothetically control host abundance, respiration, and production (16).

Regarding lake environments, little is known about the ecology of the freshwater virus when compared to marine viruses (17). Most are bacteriophages or human and other animal viruses, but plant viruses are also identified (18). Bacteriophages can play an essential function in the aquatic ecosystem as they can contribute to the acquisition and spread of antibiotic resistance genes (ARGs) (19). Some studies have shown ARGs-carrying phages to be abundant in many environments, especially those impacted by anthropogenic activities (2024), which demonstrates that these types of viruses are relevant to local microbial ecology (20).

Like bacteriophages, another essential group of DNA viruses is cyanophages, which infects cyanobacteria and have a similar morphology (25). Being abundant in both fresh and saltwater, they play an essential role in modulating cyanobacterial populations and preserving water quality (25, 26). Also, cyanophages are abundant in aquatic environments and play a fundamental role in flowering dynamics, including growth regulation and photosynthesis of cyanobacteria (27). However, unlike bacteriophages, it has many genera of possible hosts; therefore, freshwater cyanophages can be classified according to the taxonomy of their host organisms (28).

Given the need to study the diversity of different environments, techniques have been developed, such as viral metagenomics, also known as virome. This technique allows the study of various viruses from environmental samples (29). In this way, metagenomics and next-generation sequencing (NGS) have demonstrated considerable genetic complexity and inter-species and intra-species interaction by exploring viral populations both in aquatic environments and within the human microbiome (30, 31). However, despite increasing studies using the technique, there are still significant gaps in the virome databases. It has been estimated that 1,031 viral particles are infecting bacterial populations. However, <2,200 double-chain DNA virus (dsDNA) and retrovirus genomes are deposited at the National Center for Biotechnology Information (NCBI), compared to more than 45,000 bacterial genomes (32).

In order to understand the relationships that may exist between human actions and the emergence of new diseases from water sources in the Amazon, it is crucial to comprehend the role and dynamics displayed by the present viruses inside the local community. Therefore, our objective was to identify and describe the diversity of DNA viruses through viral metagenomics analyses in Lake Bolonha, especially those that have bacteriophagic and cyanophagic behavior, thus contributing to the future handling of water-borne diseases associated with resistance to antimicrobials of public health concern.

Materials and Methods

Sample Collection

All sequencing data employed in this study were generated previously by Alves et al. (33). The water samples were collected in January of 2017 at Lake Bolonha, Belém-PA, at three different points, namely: P1 (S 01°25.530″ W 048°26.043″), upstream of the Water Treatment Plant uptake; P2 (S 01°25.530″ W 048°26.018″) in the morning-glory spillway that supplies other water treatment substations, and P3 (S 01°24.992″ W 048°25.785″) in the channel connecting the lakes Água Preta and Bolonha (Figure 1). The assessment of water quality, DNA extraction, and Total Community DNA (TC-DNA) metagenomics sequencing by the Ion ProtonTM platform are described in Alves et al. (33).

FIGURE 1
www.frontiersin.org

Figure 1. Map of the lake area and identification of sampling points (P1, P2, and P3). Alves et al. (33).

Metagenomic Dataset

The raw metagenomic data used in this work is currently available at SRA/NCBI database under the access numbers SRR8893560 (P1), SRR8893561 (P2), and SRR8893559 (P3) (Table 1). This data can also be found at www.ncbi.nlm.nih.gov/sra/PRJNA506429. It is important to note that the data employed here was generated through high-throughput metagenomic sequencing of environmental samples using the Total Community DNA approach. However, we only included in our analyses the viral portion of the metagenomic dataset for this work being the virome our aimed subject of study.

TABLE 1
www.frontiersin.org

Table 1. General information about the metagenomic dataset used in this work.

Taxonomy Analysis

The raw data was used to perform the taxonomic analysis through the Kraken2 tool (34) with the parameter “–download-library viral” to download the complete viral sequences of RefSeq and classify the reads regards its taxonomy.

Assessment of Viral Diversity

The output generated by Kraken2 was submitted to the tool Bracken (35), using abundance and diversity to generate more accurate estimations on the viruses genus and species levels. The input parameters were “${CLASSIFICATION_LVL} = ‘S' (Species)” and “input data = kraken2 output (report).” All other parameters were set as default. Later, the results were displayed with the Pavian tool (36), which allows comparing the taxonomic classifications obtained by Kraken2 and Bracken and presenting abundance estimations in several samples. Veen diagrams were generated using the web-based tool InteractiVenn (37).

Viral Metagenome Assembly

The online tool Genome Detective (38) was used in default parameters to assemble the sequencing data and classify the contigs formed into their respective taxa, identified using metaSPAdes software for single-end reads (39).

Results

Virome Assembly

The viral portion of the raw metagenomic data was assembled to prepare the data for the classification of their respective taxon per each of the freshwater samples (Table 2). The species with the highest percentage of coverage after assembly were Cyanophage KBS-S-2A (P1) with 25.37% and 69.66% identity, Cladosporium fulvum T-1 virus (P2) with 23.14% coverage and 56.90% identity, and Cyanophage KBS-S-2A (P3) with 20.65% coverage and 71.18% identity.

TABLE 2
www.frontiersin.org

Table 2. List of samples (P1, P2, and P3) and their respective number of reads obtained after the sequencing, after quality control, the number of reads associated with viruses, and the number of reads used for coverage depth.

Viral Diversity of Lake Bolonha

Taxonomic analysis performed by Kraken2 revealed that there might be more than 3.500 distinct species of viruses in Lake Bolonha. A Sankey diagram is represented summarizing taxonomic diversity at samples P1, P2, and P3, respectively, in Figure 2. The complete data from the taxonomic analysis comparison between each sample site is reported as Supplementary Tables S1S3.

FIGURE 2
www.frontiersin.org

Figure 2. Sankey diagram display of the taxonomic diversity found at samples P1 (red “viruses” tag), P2 (coral “viruses” tag), and P3 (green “viruses” tag), respectively.

The 20 species with a higher reads count in all freshwater samples are represented in Figures 35. The number of reads found for each virus on this top selection was compared with their respective number of reads on the other two samples from this study. Also, a Venn diagram showcasing the overlap between these top 20 most represented viruses for samples P1, P2, and P3, is displayed in Figure 6. A similar diagram for the complete dataset of the identified viruses is reported as Supplementary Figure S1.

FIGURE 3
www.frontiersin.org

Figure 3. Bar graph showing the distribution of reads for each top 20 virophages found in sample P1, plotted against their respective representations on samples P2 and P3.

The sample P1 presented an overall abundance of Synechococcus phage, a cyanophage (Figure 3). We obtained the highest reads for a single virus in sample P2 for Choristoneura fumiferana granulovirus (Figure 4). The sample P2 also presented a much more significant amount of the genus Pandoravirus and Mimivirus compared to other collection sites. Haemophilus phage HP1 was only observed on sample P3 (Figure 5).

FIGURE 4
www.frontiersin.org

Figure 4. Bar graph showing the distribution of reads for each top 20 virophages found in sample P2, plotted against their respective representations on samples P1 and P3.

FIGURE 5
www.frontiersin.org

Figure 5. Bar graph showing the distribution of reads for each top 20 virophages found in sample P3, plotted against their respective representations on samples P1 and P2.

Discussion

Aquatic ecological environments have a broad range of viruses that are crucial in controlling bacterial communities and regulating biogeochemical cycles (40, 41). While most of the literature on water virome has concentrated on viruses in marine waters (42, 43), previous research indicated that freshwater harbored specific viral communities distinct from other aquatic environments (44). Other studies have described the characterization of freshwater viromes in ballast water (45), sewage (46, 47), lakes (48, 49), river estuaries, where marine and freshwater mix (5052), and long rivers (53). In the present study, we identified the diversity of the viral community for three different metagenomic samples from the Lake Bolonha, obtained by the previous work of Alves et al. (33), to acquire a vivid understanding of this Amazonian lake water system.

According to Alves et al. (33), the Amazonian vegetation on its shore characterizes Lake Bolonha and the propagation of large plants under its surface, resulting in eutrophication. Also, it presents increased phosphorus and total nitrogen values in physical-chemical analysis and a high fecal coliform rate (33). In this study, we noted the significant presence of cyanophages, mainly at P1 and P3 (Figure 2). This finding contributes as an essential indicator of the interference of these phages in the environment, known for their ability to perform photosynthesis by consuming oxygen and their potential for binding to nitrogen and producing toxins (54). Nitrogen or phosphorus supplies and reduced growth rate and biomass may naturally limit freshwater ecosystems, including those involving cyanobacteria accumulation (55).

When it comes to viral abundance, taxonomic analysis reveals that a plethora of more than 3.500 distinct viruses is present at freshwater Brazilian Lake Bolonha, an area of environmental preservation. It is possible to observe that the phages with a higher abundance of reads in P1 and P3 have a lower distribution among their P2 (Figures 35), which can be explained by the proximity between sampling sites P1 and P3. The alpha diversity observed in the raw sequencing data of Table 1 of Alves et al. (33) also shows significant sample abundance in P1 and P3, which could also be related to this observation.

The diversity and abundance imbalance of essential viruses, such as bacteriophages and cyanophages, can cause significant changes in the aquatic ecosystem. Some of these phages can mediate the transduction of resistance genes between bacteria, which can provide evolutionary advantages to microorganisms and affect the water quality (33). In addition, understanding this viral community's diversity could help prevent the spread of antimicrobial resistance elements and circumvent possible future multi-resistant pathogen epidemics.

Although all samples came from the same lake and environment, the overlap of viruses best represented at each sampling site showed a vast distinction when considering its top 20 viruses (Figure 6). For instance, only Escherichia phage DE3 was at the top 20 viruses for all three samples. At the top, only shared by samples P1 and P3 were five viruses: Mycobacterium phage Crossroads; Yellowstone lake phycodnavirus 3; Agrobacterium phage Atuph07; Yellowstone lake phycodnavirus 2; and Acidovorax virus ACP17. Moreover, sample P2 had no virus shared exclusively with sample P1 or P3 at its top 20. That observation shows the diversity of viral entities present in a single environment, such as Lake Bolonha. It highlights choosing different and representative sites when studying an environmental microbiome.

FIGURE 6
www.frontiersin.org

Figure 6. Venn diagram showing the overlap between the top 20 viruses for samples P1 (in red), P2 (in green), and P3 (in blue).

The sample P1 presented an abundance of Synechococcus phage (Figure 3). This phage frequently infects cyanobacteria of the Synechococcus genus in diurnal patterns of infection due to the photosynthetic activities of its host (56). This infection affects their population dynamics by killing part of this cyanobacteria population daily, estimated between 0.005 and 30% per day (27, 57, 58). In addition, it has been described that such cyanophages play a vital role in the diversity and evolution of their host cyanobacteria (18, 32, 54, 5962).

Amongst the abundant species, the cyanophage S-RIM50 was found in P1, which could infect cyanobacteria of the genus Synechococcus (63). The abundance of cyanophage S-RIM50 has been reported in both fresh and seawater (26). Cyanophages such as the cyanophage S-RIM50 and Synechococcus phage are abundant in freshwater environments. They have been isolated from various freshwater reserves, including lakes, ponds, streams, and sewage points (27, 56, 58). They have an essential contribution to maintaining the cyanobacterial community and the preservation of water quality (25).

Synechococcus phages, present in both samples P1 (Figure 3) and P3 (Figure 5), can be associated with health problems such as multiple sclerosis. This phage expresses proteins containing consensus peptide stretches that are highly homologous to the products of 16 autoantigens related to multiple sclerosis susceptibility genes (64). Other viruses associated with multiple sclerosis, such as the Epstein-Barr virus, have also shown this behavior, and the bacteriophage Synechococcus has been identified as a new relevant contributor to this phenomenon. Its cyanobacterial host prefers a temperate climate, indicating that the ecology of this cyanophage is consistent with the overall distribution and epidemiology of multiple sclerosis (63).

Yellowstone Lake phycodnaviruses, a double-stranded DNA virus that infects algae, and Escherichia phage DE3 on the sample P1 were also observed (Figure 3). A more significant number of reads associated with Shigella phage SfIV was observed compared to samples P2 and P3 (Figure 3). It is important to mention a relationship with the possible environmental presence of the Shigella bacterial host, responsible for causing intestinal infection followed or not by fever, colic, and diarrhea with blood and mucus (65). This observation demonstrates the importance of this study in characterizing the environmental conditions as a possible source of information for public health.

The P2 collection point presented many Choristoneura fumiferana granulovirus (Figure 4), part of the Baculoviridae family (66). Halovirus HF1, responsible for infecting members of the Halobacteriaceae family (67), was also identified (Figure 4). Interestingly, a small number of reads (25 reads) associated with the species Diplodia scrobiculate RNA virus 1 is only present in P2. As a preferable host, this phage has the endophytic fungus Diplodia scrobiculate, which primarily affects the genus Pinus spp. among other conifers (68), which would be odd to find on a water sample from a tropical locality.

A considerable amount of Bacillus phage Bp8pC has also been observed on sample P2 (Figure 4), which hosts the bacteria Bacillus thuringiensis and Bacillus pumilus. Both are of economic importance because they are used in agriculture as pest control bringing little harm to humans (69). The presence of the genus Bacillus in lake water may indicate the contamination of the water environment by different types of residues coming from the watershed to the lake (70).

Exclusively on P2, a much more significant amount of the genus Pandoravirus and Mimivirus was observed compared to the other collection points (Figure 4). Previous studies suggest a potential role of Mimivirus in respiratory pathology displayed during seroconversion in patients with pulmonary pneumonia. In addition, positive serology for Mimivirus is associated with increased duration of mechanical ventilation supported breathing and intensive care unit in patients with ventilator-associated pneumonia (71). Both genera are constituted of giant viruses and have Amoeba as their typical host (7274). Pandoravirus has a size of about 1 micron and may resemble some types of bacteria. Their genome contains more than 100 distinct genes and can be twice as large as the Mimivirus genome, besides the fact that their genome is quite different compared to other known organisms (75).

Haemophilus virus HP1, a bacteriophage that infects the Haemophilus influenzae bacterium (76), was only observed at sample P3 (Figure 5). Its sampling site is located at the starting point of the channel connecting both Lakes Bolonha and Água Preta. It is crucial to note that what occurs at this site may in the future influence the environment of nearby Lake Água Preta.

Overall, these results denote the presence of a diverse viral community and suggest the existence of established regulation dynamics in the local microbial environment of Lake Bolonha, highly influenced by the bacteriophages and cyanophages that inhabit the location. The dispersion of those biological entities along the water distribution channels using Lake Bolonha as a water source and general eutrophic activity might contribute to the spread of minor genetic elements like ARGs and future unbalance the microenvironment of close by freshwater sources, such as Lake Água Preta.

Conclusions

Given the importance of Lake Bolonha as a source of drinking water supply for the metropolitan region of Belém, the elucidation of the viral diversity from this environment is relevant to provide a better understanding of how its exploration can affect it. The results observed in this work indicate a widely diverse viral community, especially bacteriophages and cyanophages. These findings also suggest the existence of established micro-environmental dynamics in Lake Bolonha, possibly regulated by such phage entities. The dispersion of those viral beings bare similarity along the course of the lake, apparently more related the deeper they are into the lake (P1, P3) and the further away they are from the water evacuation sites to other treatment substations (P2). This study is the first-ever work to describe the virome of Lake Bolonha and, as such, contributes to the understanding of water-related public health concerns regarding the spreading of antibiotic resistance genes and population control of native bacteria and cyanobacteria.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/; https://www.ncbi.nlm.nih.gov/sra/PRJNA506429.

Author Contributions

BG, KP, and RR designed the study. WN, BG, and KP compiled and curated the data, performed bioinformatic analysis, and interpreted the results. AF and RR supervised and administered the project and provided funding. WN and BG wrote the original draft and manuscript with input from KP, AA, AQ, AF, and RR. All authors critically reviewed the manuscript and approved the final version.

Funding

We thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior–CAPES and PROPESP/UFPA (Pró-Reitoria de Pesquisa e Pós-Graduação / Universidade Federal do Pará), the funding agencies FAPESPA (Fundação Amazônia de Amparo à Estudos e Pesquisas), and FAPEMIG (Fundação de Amparo à Pesquisa do Estado de Minas Gerais) for financial support on this work.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2022.869886/full#supplementary-material

Supplementary Figure 1. Venn diagram representation of overlapping viruses for Samples P1, P2, and P3 complete datasets.

Supplementary Table 1. Full report table of identified viruses from taxonomic analysis performed by Kraken2 for Lake Bolonha's metagenome Sample P1.

Supplementary Table 2. Full report table of identified viruses from taxonomic analysis performed by Kraken2 for Lake Bolonha's metagenome Sample P2.

Supplementary Table 3. Full report table of identified viruses from taxonomic analysis performed by Kraken2 for Lake Bolonha's metagenome Sample P3.

Abbreviations

ARGs, Antibiotic resistance genes; NGS, Next-Generation Sequencing; dsDNA, double-chain DNA; NCBI, National Center for Biotechnology Information.

References

1. Sandy M. Why Is the Amazon Rain Forest Disappearing? Time (2019).

Google Scholar

2. Baker JCA, Spracklen DV. Climate benefits of intact amazon forests and the biophysical consequences of disturbance, front. For Glob Chang. (2019) 2:47. doi: 10.3389/ffgc.2019.00047

CrossRef Full Text | Google Scholar

3. Metzger JP, Bustamante MMC, Ferreira J, Fernandes GW, Librán-Embid F, Pillar VD, et al. Why brazil needs its legal reserves, perspect. Ecol Conserv. (2019) 17:91–103. doi: 10.1016/j.pecon.2019.07.002

CrossRef Full Text | Google Scholar

4. INPE PRODES — Coordenação-Geral de Observação da Terra Inst. Nac. Pesqui. Espac. (2020). p. 1. Available online at: http://www.obt.inpe.br/OBT/assuntos/programas/amazonia/prodes (accessed January 27, 2022).

5. Ellwanger JH, Kulmann-Leal B, Kaminski VL, Valverde-Villegas JM, VEIGA ABGDA, Spilki FR, et al. Beyond diversity loss and climate change: Impacts of Amazon deforestation on infectious diseases and public health. An Acad Bras Cienc. (2020) 92:20191375. doi: 10.1590/0001-3765202020191375

PubMed Abstract | CrossRef Full Text | Google Scholar

6. de Cássio Veloso de Barros B, Chagas EN, Bezerra LW, Ribeiro LG, Barboza Duarte JW, Pereira D, et al. Rotavirus A in wild and domestic animals from areas with environmental degradation in the Brazilian Amazon. PLoS ONE. (2018) 13:e0209005. doi: 10.1371/journal.pone.0209005

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Silva MBda, Portela JM, Li W, Jackson M, Gonzalez-Juarrero M, Hidalgo AS, et al. Evidence of zoonotic leprosy in Pará, Brazilian Amazon, and risks associated with human contact or consumption of armadillos. PLoS Negl Trop Dis. (2018) 12:e0006532. doi: 10.1371/journal.pntd.0006532

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Farikoski IO, Medeiros LS, Carvalho YK, Ashford DA, Figueiredo EES, Fernandes DVGS, et al. The urban and rural capybaras (Hydrochoerus hydrochaeris) as reservoir of Salmonella in the western Amazon, Brazil. Pesqui Vet Bras. (2019) 39:66–9. doi: 10.1590/1678-5150-pvb-5761

CrossRef Full Text | Google Scholar

9. Franco Filho LC, Barata RR, Cardoso JF, de V., Massafra JM, da S. Lemos P, Casseb LMN, et al. Metagenomic analysis of samples from three bat species collected in the amazon rain forest. Microbiol Resour Announc. (2019) 8:e01422–18. doi: 10.1128/MRA.01422-18

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Fernandes J, Guterres A, De Oliveira RC, Chamberlain J, Lewandowski K, Teixeira BR, et al. Xapuri virus, a novel mammarenavirus: Natural reassortment and increased diversity between New World viruses. Emerg Microbes Infect. (2018) 7:120. doi: 10.1038/s41426-018-0119-9

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Medeiros DBA, Vasconcelos PFC. Is the brazilian diverse environment is a crib for the emergence and maintenance of exotic arboviruses? An Acad Bras Cienc. (2019) 91:e20190407. doi: 10.1590/0001-3765201920190407

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Confalonieri UEC, Margonari C, Quintão AF. Environmental change and the dynamics of parasitic diseases in the Amazon. Acta Trop. (2014) 129:33–41. doi: 10.1016/j.actatropica.2013.09.013

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Fonseca PAM, Hacon Sde S, Reis VL, Costa D, Brown IF. O uso de dados de satélite para estudar a relação entre chuva e doenças diarreicas em uma bacia na amazônia sul-ocidental. Cienc Saude Coletiva. (2016) 21:731–42. doi: 10.1590/1413-81232015213.20162015

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Levy K, Smith SM, Carlton EJ. Climate change impacts on waterborne diseases: moving toward designing interventions. Curr Environ Heal Reports. (2018) 5:272–82. doi: 10.1007/s40572-018-0199-7

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Putt AE, MacIsaac EA, Herunter HE, Cooper AB, Selbie DT. Eutrophication forcings on a peri-urban lake ecosystem: context for integrated watershed to airshed management. PLoS ONE. (2019) 14:e0219241. doi: 10.1371/journal.pone.0219241

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Cook A, Rushton S, Allan JA, Baxter A. An evaluation of techniques to control problem bird species on landfill sites. Environ Manage. (2008) 41:834–43. doi: 10.1007/s00267-008-9077-7

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Palermo CN, Fulthorpe RR, Saati R, Short SM. Metagenomic analysis of virus diversity and relative abundance in a eutrophic freshwater harbour. Viruses. (2019) 11:792. doi: 10.3390/v11090792

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Mohiuddin M, Schellhorn HE. Spatial and temporal dynamics of virus occurrence in two freshwater lakes captured through metagenomic analysis. Front Microbiol. (2015) 6:960. doi: 10.3389/fmicb.2015.00960

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Koskella B, Brockhurst MA. Bacteria-phage coevolution as a driver of ecological and evolutionary processes in microbial communities. FEMS Microbiol Rev. (2014) 38:916–31. doi: 10.1111/1574-6976.12072

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Calero-Cáceres W, Melgarejo A, Colomer-Lluch M, Stoll C, Lucena F, Jofre J, et al. Sludge as a potential important source of antibiotic resistance genes in both the bacterial and bacteriophage fractions. Environ Sci Technol. (2014) 48:7602–11. doi: 10.1021/es501851s

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Colomer-Lluch M, Calero-Cáceres W, Jebri S, Hmaied F, Muniesa M, Jofre J. Antibiotic resistance genes in bacterial and bacteriophage fractions of Tunisian and Spanish wastewaters as markers to compare the antibiotic resistance patterns in each population. Environ Int. (2014) 73:167–75. doi: 10.1016/j.envint.2014.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Lekunberri I, Villagrasa M, Balcázar JL, Borrego CM. Contribution of bacteriophage and plasmid DNA to the mobilization of antibiotic resistance genes in a river receiving treated wastewater discharges. Sci Total Environ. (2017) 601–2:206–9. doi: 10.1016/j.scitotenv.2017.05.174

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Marti E, Variatza E, Balcázar JL. Bacteriophages as a reservoir of extended-spectrum β-lactamase and fluoroquinolone resistance genes in the environment. Clin Microbiol Infect. (2014) 20:0456–9. doi: 10.1111/1469-0691.12446

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Subirats J, Sànchez-Melsió A, Borrego CM, Balcázar JL, Simonet P. Metagenomic analysis reveals that bacteriophages are reservoirs of antibiotic resistance genes. Int J Antimicrob Agents. (2016) 48:163–7. doi: 10.1016/j.ijantimicag.2016.04.028

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Xia H, Li T, Deng F, Hu Z. Freshwater cyanophages. Virol Sin. (2013) 28:253–9. doi: 10.1007/s12250-013-3370-1

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Shane MS. Distribution of blue-green algal viruses in various types of natural waters. Water Res. (1971) 5:711–6. doi: 10.1016/0043-1354(71)90093-5

CrossRef Full Text | Google Scholar

27. Jiang X, Ha C, Lee S, Kwon J, Cho H, Gorham T, et al. Characterization of cyanophages in lake erie: Interaction mechanisms and structural damage of toxic cyanobacteria. Toxins. (2019) 11:444. doi: 10.3390/toxins11080444

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Singh P. Structural, functional and molecular basis of cyanophage-cyanobacterial interactions and its significance. AFRICAN J Biotechnol. (2012) 11:2591–608. doi: 10.5897/AJB10.790

CrossRef Full Text | Google Scholar

29. Wommack KE, Nasko DJ, Chopyk J, Sakowski EG. Counts and sequences, observations that continue to change our understanding of viruses in nature. J Microbiol. (2015) 53:181–92. doi: 10.1007/s12275-015-5068-6

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Aggarwala V, Liang G, Bushman FD. Viral communities of the human gut: Metagenomic analysis of composition and dynamics. Mob DNA. (2017) 8:1–10. doi: 10.1186/s13100-017-0095-y

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Coutinho FH, Silveira CB, Gregoracci GB, Thompson CC, Edwards RA, Brussaard CPD, et al. Marine viruses discovered via metagenomics shed light on viral strategies throughout the oceans. Nat Commun. (2017) 8:1–12. doi: 10.1038/s41467-016-0009-6

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Paez-Espino D, Eloe-Fadrosh EA, Pavlopoulos GA, Thomas AD, Huntemann M, Mikhailova N, et al. Uncovering Earth's virome. Nature. (2016) 536:425–30. doi: 10.1038/nature19094

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Alves J, Dias L, Mateus J, Marques J, Graças D, Ramos R, et al. Resistome in lake bolonha, brazilian amazon: identification of genes related to resistance to broad-spectrum antibiotics. Front Microbiol. (2020) 11:67. doi: 10.3389/fmicb.2020.00067

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol. (2019) 20:257. doi: 10.1186/s13059-019-1891-0

PubMed Abstract | CrossRef Full Text

35. Lu J, Breitwieser FP, Thielen P, Salzberg SL. Bracken: estimating species abundance in metagenomics data. PeerJ Comput Sci. (2017) 2017:e104. doi: 10.7717/peerj-cs.104

CrossRef Full Text

36. Breitwieser FP, Salzberg SL. Pavian: interactive analysis of metagenomics data for microbiome studies and pathogen identification. Bioinformatics. (2020) 36:1303–4. doi: 10.1093/bioinformatics/btz715

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Heberle H, Meirelles VG, Silva FRda, Telles GP, Minghim R. InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams. BMC Bioinformatics. (2015) 16:1–7. doi: 10.1186/s12859-015-0611-3

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Vilsker M, Moosa Y, Nooij S, Fonseca V, Ghysens Y, Dumon K, et al. Genome detective: an automated system for virus identification from high-throughput sequencing data. Bioinformatics. (2019) 35:871–3. doi: 10.1093/bioinformatics/bty695

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. MetaSPAdes: a new versatile metagenomic assembler. Genome Res. (2017) 27:824–34. doi: 10.1101/gr.213959.116

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Fischer UR, Velimirov B. High control of bacterial production by viruses in a eutrophic oxbow lake. Aquat Microb Ecol. (2002) 27:1–12. doi: 10.3354/ame027001

CrossRef Full Text | Google Scholar

41. Jover LF, Effler TC, Buchan A, Wilhelm SW, Weitz JS. The elemental composition of virus particles: implications for marine biogeochemical cycles. Nat Rev Microbiol. (2014) 12:519–28. doi: 10.1038/nrmicro3289

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Suttle CA. Environmental microbiology: viral diversity on the global stage. Nat Microbiol. (2016) 1:16205. doi: 10.1038/nmicrobiol.2016.205

PubMed Abstract | CrossRef Full Text

43. Gregory AC, Zayed AA, Conceição-Neto N, Temperton B, Bolduc B, Alberti A, et al. Marine DNA viral macro- and microdiversity from pole to pole. Cell. (2019) 177:1109–23.e14. doi: 10.1016/j.cell.2019.03.040

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Roux S, Enault F, Robin A, Ravet V, Personnic S, Theil S, et al. Assessing the diversity and specificity of two freshwater viral communities through metagenomics. PLoS ONE. (2012) 7:e33641. doi: 10.1371/journal.pone.0033641

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Kim Y, Aw TG, Teal TK, Rose JB. Metagenomic Investigation of viral communities in ballast water. Environ Sci Technol. (2015) 49:8396–407. doi: 10.1021/acs.est.5b01633

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Fernandez-Cassi X, Timoneda N, Martínez-Puchol S, Rusiñol M, Rodriguez-Manzano J, Figuerola N, et al. Metagenomics for the study of viruses in urban sewage as a tool for public health surveillance. Sci Total Environ. (2018) 618:870–80. doi: 10.1016/j.scitotenv.2017.08.249

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Martínez-Puchol S, Rusiñol M, Fernández-Cassi X, Timoneda N, Itarte M, Andrés C, et al. Characterisation of the sewage virome: comparison of NGS tools and occurrence of significant pathogens. Sci Total Environ. (2020) 713:136604. doi: 10.1016/j.scitotenv.2020.136604

PubMed Abstract | CrossRef Full Text | Google Scholar

48. López-Bueno J., Tamames D., Velázquez A., Moya A., Quesada A., Alcamí. High diversity of the viral community from an Antarctic lake. Science. (2009) 326:858–61. doi: 10.1126/science.1179287

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Ge X, Wu Y, Wang M, Wang J, Wu L, Yang X, et al. Viral metagenomics analysis of planktonic viruses in East Lake, Wuhan, China. Virol Sin. (2013) 28:280–90. doi: 10.1007/s12250-013-3365-y

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Cai L, Zhang R, He Y, Feng X, Jiao N. Metagenomic analysis of Virioplankton of the subtropical Jiulong river estuary, China. Viruses. (2016) 8:35. doi: 10.3390/v8020035

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Wolf YI, Silas S, Wang Y, Wu S, Bocek M, Kazlauskas D, et al. Doubling of the known set of RNA viruses by metagenomic analysis of an aquatic virome. Nat Microbiol. (2020) 5:1262–70. doi: 10.1038/s41564-020-0755-4

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Wu S, Zhou L, Zhou Y, Wang H, Xiao J, Yan S, et al. Diverse and unique viruses discovered in the surface water of the East China Sea. BMC Genomics. (2020) 21:441. doi: 10.1186/s12864-020-06861-y

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Lu J, Yang S, Zhang X, Tang X, Zhang J, Wang X, et al. Metagenomic analysis of viral community in the yangtze river expands known eukaryotic and prokaryotic virus diversity in freshwater. Virol Sin. (2022) 37:60–9. doi: 10.1016/j.virs.2022.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Jaskulska J. Mankiewicz-Boczek Cyanophages specific to cyanobacteria from the genus Microcystis. Ecohydrol Hydrobiol. (2020) 20:83–90. doi: 10.1016/j.ecohyd.2019.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Cheng K, Frenken T, Brussaard CPD, Van De Waal DB. Cyanophage propagation in the freshwater cyanobacterium phormidiumis constrained by phosphorus limitation and enhanced by elevated pCO2. Front Microbiol. (2019) 10:617. doi: 10.3389/fmicb.2019.00617

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Liu R, Liu Y, Chen Y, Zhan Y, Zeng Q. Cyanobacterial viruses exhibit diurnal rhythms during infection. Proc Natl Acad Sci U S A. (2019) 116:14077–82. doi: 10.1073/pnas.1819689116

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Bragg L, Tyson GW. Metagenomics using next-generation sequencing. Methods Mol Biol. (2014) 1096:183–201. doi: 10.1007/978-1-62703-712-9_15

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Toyama D, Kishi LT, Santos-Júnior CD, Soares-Costa A, Souza de Oliveira TC, de Miranda FP, et al. Metagenomics analysis of microorganisms in freshwater lakes of the Amazon Basin. Genome Announc. (2016) 4:1440–56. doi: 10.1128/genomeA.01440-16

PubMed Abstract | CrossRef Full Text

59. Handelsman J. Metagenomics: application of genomics to uncultured microorganisms, microbiol. Mol Biol Rev. (2004) 68:669–85. doi: 10.1128/MMBR.68.4.669-685.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Lorenz P, Eck J. Metagenomics and industrial applications. Nat Rev Microbiol. (2005) 3:510–6. doi: 10.1038/nrmicro1161

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Tundisi JG. Água no século 21: enfrentando a escassez (2003).

Google Scholar

62. Fedida A, Lindell D. Two synechococcus genes, two different effects on cyanophage infection. Viruses. (2017) 9:136. doi: 10.3390/v9060136

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Stoddard LI, Martiny JBH, Marston MF. Selection and characterization of cyanophage resistance in marine synechococcus strains. Appl Environ Microbiol. (2007) 73:5516–22. doi: 10.1128/AEM.00356-07

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Carter CJ. Epstein-Barr and other viral mimicry of autoantigens, myelin and vitamin D-related proteins and of EIF2B, the cause of vanishing white matter disease: massive mimicry of multiple sclerosis relevant proteins by the Synechococcus phage. Immunopharmacol Immunotoxicol. (2012) 34:21–35. doi: 10.3109/08923973.2011.572262

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Martinez-Becerra FJ, Chen X, Dickenson NE, Choudhari SP, Harrison K, Clements JD, et al. Characterization of a novel fusion protein from IpaB and IpaD of Shigella spp and its potential as a Pan-Shigella vaccine. Infect Immun. (2013) 81:4470–7. doi: 10.1128/IAI.00859-13

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Clem RJ, Passarelli AL. Baculoviruses: sophisticated pathogens of insects. PLoS Pathog. (2013) 9:e1003729. doi: 10.1371/journal.ppat.1003729

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Ventosa. Halophilic Microorganisms. Berlin; Heidelberg: Springer (2004). doi: 10.1007/978-3-662-07656-9

CrossRef Full Text | Google Scholar

68. Bihon W, Slippers B, Burgess T, Wingfield MJ, Wingfield BD. Diplodia scrobiculata found in the southern hemisphere. For Pathol. (2011) 41:175–81. doi: 10.1111/j.1439-0329.2010.00649.x

CrossRef Full Text | Google Scholar

69. Stotzky G. Persistence and biological activity in soil of insecticidal proteins from bacillus thuringiensis and of bacterial DNA bound on clays and humic acids. J Environ Qual. (2000) 29:691–705. doi: 10.2134/jeq2000.00472425002900030003x

CrossRef Full Text | Google Scholar

70. Shafi S, Kamili AN, Shah MA, Bandh SA, Dar R. Dynamics of bacterial class Bacilli in the deepest valley lake of Kashmir-the Manasbal Lake. Microb Pathog. (2017) 104:78–83. doi: 10.1016/j.micpath.2017.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Vanspauwen MJ, Franssen FME, Raoult D, Wouters EFM, Bruggeman CA, Linssen CFM. Infections with mimivirus in patients with chronic obstructive pulmonary disease. Respir Med. (2012) 106:1690–4. doi: 10.1016/j.rmed.2012.08.019

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Legendre M, Fabre E, Poirot O, Jeudy S, Lartigue A, Alempic JM, et al. Diversity and evolution of the emerging pandoraviridae family. Nat Commun. (2018) 9:1–12. doi: 10.1038/s41467-017-02088-w

PubMed Abstract | CrossRef Full Text | Google Scholar

73. La Scola B, Audic S, Robert C, Jungang L, De Lamballerie X, Drancourt M, et al. A giant virus in amoebae. Science. (2003) 299:2033. doi: 10.1126/science.1081867

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Raoult D, Audic S, Robert C, Abergel C, Renesto P, Ogata H, et al. The 1.2-megabase genome sequence of mimivirus. Science. (2004) 306:1344–50. doi: 10.1126/science.1101485

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Abergel C, Legendre M, Claverie JM. The rapidly expanding universe of giant viruses: mimivirus, pandoravirus, pithovirus and mollivirus. FEMS Microbiol Rev. (2015) 39:779–96. doi: 10.1093/femsre/fuv037

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Hakimi JM, Scocca JJ. Purification and characterization of the integrase from the Haemophilus influenzae bacteriophage HP1; identification of a four-stranded intermediate and the order of strand exchange. Mol Microbiol. (1996) 21:147–58. doi: 10.1046/j.1365-2958.1996.6311351.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Amazon, Lake Bolonha, freshwater, viruses, bacteriophages, cyanophages, virome, metagenomics

Citation: Nogueira WG, Gois BVA, Pinheiro KdC, Aragão AdO, Queiroz ALC, Silva ALd, Folador AC and Ramos RTJ (2022) Viral Metagenomics Reveals Widely Diverse Viral Community of Freshwater Amazonian Lake. Front. Public Health 10:869886. doi: 10.3389/fpubh.2022.869886

Received: 05 February 2022; Accepted: 24 March 2022;
Published: 25 April 2022.

Edited by:

Maria Ines Zanoli Sato, Companhia Ambiental do Estado de São Paulo (CETESB), Brazil

Reviewed by:

Niti B. Jadeja, Ashoka Trust for Research in Ecology and the Environment (ATREE), India
Mahmuda Yasmin, University of Dhaka, Bangladesh

Copyright © 2022 Nogueira, Gois, Pinheiro, Aragão, Queiroz, Silva, Folador and Ramos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Rommel Thiago Jucá Ramos, rommelthiago@gmail.com

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.