Wastewater surveillance uncovers regional diversity and dynamics of SARS-CoV-2 variants across nine states in the USA

Wastewater-based epidemiology (WBE) is a non-invasive and cost-effective approach for monitoring the spread of a pathogen within a community. WBE has been adopted as one of the methods to monitor the spread and population dynamics of the SARS-CoV-2 virus, but significant challenges remain in the bioinformatic analysis of WBE-derived data. Here, we have developed a new distance metric, CoVdist, and an associated analysis tool that facilitates the application of ordination analysis to WBE data and the identification of viral population changes based on nucleotide variants. We applied these new approaches to a large-scale dataset from 18 cities in nine states of the USA using wastewater collected from July 2021 to June 2022. We found that the trends in the shift between the Delta and Omicron SARS-CoV-2 lineages were largely consistent with what was seen in clinical data, but that wastewater analysis offered the added benefit of revealing significant differences in viral population dynamics at the state, city, and even neighborhood scales. We also were able to observe the early spread of variants of concern and the presence of recombinant lineages during the transitions between variants, both of which are challenging to analyze based on clinically-derived viral genomes. The methods outlined here will be beneficial for future applications of WBE to monitor SARS-CoV-2, particularly as clinical monitoring becomes less prevalent. Additionally, these approaches are generalizable, allowing them to be applied for the monitoring and analysis of future viral outbreaks.

Evolution of anelloviruses from a circovirus-like ancestor through gradual augmentation of the jelly-roll capsid protein

Anelloviruses are highly prevalent in diverse mammals, including humans, but so far have not been linked to any disease and are considered to be part of the 'healthy virome'. These viruses have small circular single-stranded DNA (ssDNA) genomes and encode several proteins with no detectable sequence similarity to proteins of other known viruses. Thus, anelloviruses are the only family of eukaryotic ssDNA viruses currently not included in the realm Monodnaviria. To gain insights into the provenance of these enigmatic viruses, we sequenced more than 250 complete genomes of anelloviruses from nasal and vaginal swab samples of Weddell seal (Leptonychotes weddellii) from Antarctica and a fecal sample of grizzly bear (Ursus arctos horribilis) from the USA and performed a comprehensive family-wide analysis of the signature anellovirus protein ORF1. Using state-of-the-art remote sequence similarity detection approaches and structural modeling with AlphaFold2, we show that ORF1 orthologs from all Anelloviridae genera adopt a jelly-roll fold typical of viral capsid proteins (CPs), establishing an evolutionary link to other eukaryotic ssDNA viruses, specifically, circoviruses. However, unlike CPs of other ssDNA viruses, ORF1 encoded by anelloviruses from different genera display remarkable variation in size, due to insertions into the jelly-roll domain. In particular, the insertion between β-strands H and I forms a projection domain predicted to face away from the capsid surface and function at the interface of virus-host interactions. Consistent with this prediction and supported by recent experimental evidence, the outermost region of the projection domain is a mutational hotspot, where rapid evolution was likely precipitated by the host immune system. Collectively, our findings further expand the known diversity of anelloviruses and explain how anellovirus ORF1 proteins likely diverged from canonical jelly-roll CPs through gradual augmentation of the projection domain. We suggest assigning Anelloviridae to a new phylum, 'Commensaviricota', and including it into the kingdom Shotokuvirae (realm Monodnaviria), alongside Cressdnaviricota and Cossaviricota.

DNA virome composition of two sympatric wild felids, bobcat (Lynx rufus) and puma (Puma concolor) in Sonora, Mexico

With viruses often having devastating effects on wildlife population fitness and wild mammals serving as pathogen reservoirs for potentially zoonotic diseases, determining the viral diversity present in wild mammals is both a conservation and One Health priority. Additionally, transmission from more abundant hosts could increase the extinction risk of threatened sympatric species. We leveraged an existing circular DNA enriched metagenomic dataset generated from bobcat (Lynx rufus, n = 9) and puma (Puma concolor, n = 13) scat samples non-invasively collected from Sonora, Mexico, to characterize fecal DNA viromes of each species and determine the extent that viruses are shared between them. Using the metaWRAP pipeline to co-assemble viral genomes for comparative metagenomic analysis, we observed diverse circular DNA viruses in both species, including circoviruses, genomoviruses, and anelloviruses. We found that differences in DNA virome composition were partly attributed to host species, although there was overlap between viruses in bobcats and pumas. Pumas exhibited greater levels of alpha diversity, possibly due to bioaccumulation of pathogens in apex predators. Shared viral taxa may reflect dietary overlap, shared environmental resources, or transmission through host interactions, although we cannot rule out species-specific host-virus coevolution for the taxa detected through co-assembly. However, our detection of integrated feline foamy virus (FFV) suggests Sonoran pumas may interact with domestic cats. Our results contribute to the growing baseline knowledge of wild felid viral diversity. Future research including samples from additional sources (e.g., prey items, tissues) may help to clarify host associations and determine the pathogenicity of detected viruses.

Leveraging an established neighbourhood-level, open access wastewater monitoring network to address public health priorities: a population-based study

BACKGROUND

Before the COVID-19 pandemic, the US opioid epidemic triggered a collaborative municipal and academic effort in Tempe, Arizona, which resulted in the world's first open access dashboard featuring neighbourhood-level trends informed by wastewater-based epidemiology (WBE). This study aimed to showcase how wastewater monitoring, once established and accepted by a community, could readily be adapted to respond to newly emerging public health priorities.

METHODS

In this population-based study in Greater Tempe, Arizona, an existing opioid monitoring WBE network was modified to track SARS-CoV-2 transmission through the analysis of 11 contiguous wastewater catchments. Flow-weighted and time-weighted 24 h composite samples of untreated wastewater were collected at each sampling location within the wastewater collection system for 3 days each week (Tuesday, Thursday, and Saturday) from April 1, 2020, to March 31, 2021 (Area 7 and Tempe St Luke's Hospital were added in July, 2020). Reverse transcription quantitative PCR targeting the E gene of SARS-CoV-2 isolated from the wastewater samples was used to determine the number of genome copies in each catchment. Newly detected clinical cases of COVID-19 by zip code within the City of Tempe, Arizona were reported daily by the Arizona Department of Health Services from May 23, 2020. Maricopa County-level new positive cases, COVID-19-related hospitalisations, deaths, and long-term care facility deaths per day are publicly available and were collected from the Maricopa County Epidemic Curve Dashboard. Viral loads of SARS-CoV-2 (genome copies per day) measured in wastewater from each catchment were aggregated at the zip code level and city level and compared with the clinically reported data using root mean square error to investigate early warning capability of WBE.

FINDINGS

Between April 1, 2020, and March 31, 2021, 1556 wastewater samples were analysed. Most locations showed two waves in viral levels peaking in June, 2020, and December, 2020-January, 2021. An additional wave of viral load was seen in catchments close to Arizona State University (Areas 6 and 7) at the beginning of the fall (autumn) semester in late August, 2020. Additionally, an early infection hotspot was detected in the Town of Guadalupe, Arizona, starting the week of May 4, 2020, that was successfully mitigated through targeted interventions. A shift in early warning potential of WBE was seen, from a leading (mean of 8·5 days [SD 2·1], June, 2020) to a lagging (-2·0 days [1·4], January, 2021) indicator compared with newly reported clinical cases.

INTERPRETATION

Lessons learned from leveraging an existing neighbourhood-level WBE reporting dashboard include: (1) community buy-in is key, (2) public data sharing is effective, and (3) sub-ZIP-code (postal code) data can help to pinpoint populations at risk, track intervention success in real time, and reveal the effect of local clinical testing capacity on WBE's early warning capability. This successful demonstration of transitioning WBE efforts from opioids to COVID-19 encourages an expansion of WBE to tackle newly emerging and re-emerging threats (eg, mpox and polio).

FUNDING

National Institutes of Health's RADx-rad initiative, National Science Foundation, Virginia G Piper Charitable Trust, J M Kaplan Fund, and The Flinn Foundation.

Molecular identification and phylogenetic characterization of A-strain isolates of maize streak virus from western Ethiopia

The A-strain of maize streak virus (MSV) causes maize streak disease (MSD), which is a major biotic threat to maize production in sub-Saharan Africa. Previous studies have described different MSV strains of economic importance from southern and eastern African countries and how eastern African regions are hubs for MSV diversification. Despite these efforts, due to a lack of extensive sampling, there is limited knowledge about the MSV-A diversity in Ethiopia. Here, field sampling of maize plants and wild grasses with visible MSD symptoms was carried out in the western Ethiopian regions of Gambela, Oromia, and Benishangul-Gumuz during the maize-growing season of 2019. The complete genomes of MSV isolates (n = 60) were cloned and sequenced by the Sanger method. We used a model-based phylogenetic approach to analyse 725 full MSV genome sequences available in the GenBank database together with newly determined genome sequences from Ethiopia to determine their subtypes and identify recombinant lineages. Of the 127 fields accessed, MSD prevalence was highest, at 96%, in the Gambela region and lowest in Oromia, at 66%. The highest mean symptom severity of 4/5 (where 5 is the highest and 1 the lowest) was observed in Gambela and Benishangul-Gumuz. Our results show that these newly determined MSV isolates belong to recombinant lineage V of the A₁ subtype, with the widest dissemination and greatest economic significance in sub-Saharan Africa and the adjacent Indian Ocean islands.

Virion-Associated Nucleic Acid-Based Metagenomics: A Decade of Advances in Molecular Characterization of Plant Viruses

Over the last decade, viral metagenomic studies have resulted in the discovery of thousands of previously unknown viruses. These studies are likely to play a pivotal role in obtaining an accurate and robust understanding of how viruses affect the stability and productivity of ecosystems. Among the metagenomics-based approaches that have been developed since the beginning of the 21st century, shotgun metagenomics applied specifically to virion-associated nucleic acids (VANA) has been used to disentangle the diversity of the viral world. We summarize herein the results of 24 VANA-based studies, focusing on plant and insect samples conducted over the last decade (2010 to 2020). Collectively, viruses from 85 different families were reliably detected in these studies, including capsidless RNA viruses that replicate in fungi, oomycetes, and plants. Finally, strengths and weaknesses of the VANA approach are summarized and perspectives of applications in detection, epidemiological surveillance, environmental monitoring, and ecology of plant viruses are provided. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A parasite outbreak in notothenioid fish in an Antarctic fjord

Climate changes can promote disease outbreaks, but their nature and potential impacts in remote areas have received little attention. In a hot spot of biodiversity on the West Antarctic Peninsula, which faces among the fastest changing climates on Earth, we captured specimens of two notothenioid fish species affected by large skin tumors at an incidence never before observed in the Southern Ocean. Molecular and histopathological analyses revealed that X-cell parasitic alveolates, members of a genus we call Notoxcellia, are the etiological agent of these tumors. Parasite-specific molecular probes showed that xenomas remained within the skin but largely outgrew host cells in the dermis. We further observed that tumors induced neovascularization in underlying tissue and detrimentally affected host growth and condition. Although many knowledge gaps persist about X-cell disease, including its mode of transmission and life cycle, these findings reveal potentially active biotic threats to vulnerable Antarctic ecosystems.

Discovery of novel fish papillomaviruses: From the Antarctic to the commercial fish market

Fish papillomaviruses form a newly discovered group broadly recognized as the Secondpapillomavirinae subfamily. This study expands the documented genomes of the fish papillomaviruses from six to 16, including one from the Antarctic emerald notothen, seven from commercial market fishes, one from data mining of sea bream sequence data, and one from a western gull cloacal swab that is likely diet derived. The genomes of secondpapillomaviruses are ∼6 kilobasepairs (kb), which is substantially smaller than the ∼8 kb of terrestrial vertebrate papillomaviruses. Each genome encodes a clear homolog of the four canonical papillomavirus genes, E1, E2, L1, and L2. In addition, we identified open reading frames (ORFs) with short linear peptide motifs reminiscent of E6/E7 oncoproteins. Fish papillomaviruses are extremely diverse and phylogenetically distant from other papillomaviruses suggesting a model in which terrestrial vertebrate-infecting papillomaviruses arose after an evolutionary bottleneck event, possibly during the water-to-land transition.

Agricultural practices drive biological loads, seasonal patterns and potential pathogens in the aerobiome of a mixed-land-use dryland

Air carries a diverse load of particulate microscopic biological matter in suspension, either aerosolized or aggregated with dust particles, the aerobiome, which is dispersed by winds from sources to sinks. The aerobiome is known to contain microbes, including pathogens, as well as debris or small-sized propagules from plants and animals, but its variability and composition has not been studied comprehensibly. To gain a dynamic insight into the aerobiome existing over a mixed-use dryland setting, we conducted a biologically comprehensive, year-long survey of its composition and dynamics for particles less than 10 μm in diameter based on quantitative analyses of DNA content coupled to genomic sequencing. Airborne biological loads were more dependent on seasonal events than on meteorological conditions and only weakly correlated with dust loads. Core aerobiome species could be understood as a mixture of high elevation (e.g. Microbacteriaceae, Micrococcaceae, Deinococci), and local plant and soil sources (e.g. Sphingomonas, Streptomyces, Acinetobacter). Despite the mixed used of the land surrounding the sampling site, taxa that contributed to high load events were largely traceable to proximal agricultural practices like cotton and livestock farming. This included not only the predominance of specific crop plant signals over those of native vegetation, but also that of their pathogens (bacterial, viral and eukaryotic). Faecal bacterial loads were also seasonally important, possibly sourced in intensive animal husbandry or manure fertilization activity, and this microbial load was enriched in tetracycline resistance genes. The presence of the native opportunistic pathogen, Coccidioides spp., by contrast, was detected only with highly sensitive techniques, and only rarely. We conclude that agricultural activity exerts a much stronger influence that the native vegetation as a mass loss factor to the land system and as an input to dryland aerobiomes, including in the dispersal of plant, animal and human pathogens and their genetic resistance characteristics.

Novel circular DNA virus identified in Opuntia discolor (Cactaceae) that codes for proteins with similarity to those of geminiviruses

Viral metagenomic studies have enabled the discovery of many unknown viruses and revealed that viral communities are much more diverse and ubiquitous than previously thought. Some viruses have multiple genome components that are encapsidated either in separate virions (multipartite viruses) or in the same virion (segmented viruses). In this study, we identify what is possibly a novel bipartite plant-associated circular single-stranded DNA virus in a wild prickly pear cactus, Opuntia discolor, that is endemic to the Chaco ecoregion in South America. Two ~1.8 kb virus-like circular DNA components were recovered, one encoding a replication-associated protein (Rep) and the other a capsid protein (CP). Both of the inferred protein sequences of the Rep and CP are homologous to those encoded by members of the family Geminiviridae. These two putatively cognate components each have a nonanucleotide sequence within a likely hairpin structure that is homologous to the origins of rolling-circle replication (RCR), found in diverse circular single-stranded DNA viruses. In addition, the two components share similar putative replication-associated iterative sequences (iterons), which in circular single-stranded DNA viruses are important for Rep binding during the initiation of RCR. Such molecular features provide support for the possible bipartite nature of this virus, which we named utkilio virus (common name of the Opuntia discolor in South America) components A and B. In the infectivity assays conducted in Nicotiana benthamiana plants, only the A component of utkilio virus, which encodes the Rep protein, was found to move and replicate systemically in N. benthamiana. This was not true for component B, for which we did not detect replication, which may have been due to this being a defective molecule or because of the model plants (N. benthamiana) used for the infection assays. Future experiments need to be conducted with other plants, including O. discolor, to understand more about the biology of these viral components.

High-throughput sequencing of SARS-CoV-2 in wastewater provides insights into circulating variants

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) likely emerged from a zoonotic spill-over event and has led to a global pandemic. The public health response has been predominantly informed by surveillance of symptomatic individuals and contact tracing, with quarantine, and other preventive measures have then been applied to mitigate further spread. Non-traditional methods of surveillance such as genomic epidemiology and wastewater-based epidemiology (WBE) have also been leveraged during this pandemic. Genomic epidemiology uses high-throughput sequencing of SARS-CoV-2 genomes to inform local and international transmission events, as well as the diversity of circulating variants. WBE uses wastewater to analyse community spread, as it is known that SARS-CoV-2 is shed through bodily excretions. Since both symptomatic and asymptomatic individuals contribute to wastewater inputs, we hypothesized that the resultant pooled sample of population-wide excreta can provide a more comprehensive picture of SARS-CoV-2 genomic diversity circulating in a community than clinical testing and sequencing alone. In this study, we analysed 91 wastewater samples from 11 states in the USA, where the majority of samples represent Maricopa County, Arizona (USA). With the objective of assessing the viral diversity at a population scale, we undertook a single-nucleotide variant (SNV) analysis on data from 52 samples with >90% SARS-CoV-2 genome coverage of sequence reads, and compared these SNVs with those detected in genomes sequenced from clinical patients. We identified 7973 SNVs, of which 548 were "novel" SNVs that had not yet been identified in the global clinical-derived data as of 17^th June 2020 (the day after our last wastewater sampling date). However, between 17^th of June 2020 and 20^th November 2020, almost half of the novel SNVs have since been detected in clinical-derived data. Using the combination of SNVs present in each sample, we identified the more probable lineages present in that sample and compared them to lineages observed in North America prior to our sampling dates. The wastewater-derived SARS-CoV-2 sequence data indicates there were more lineages circulating across the sampled communities than represented in the clinical-derived data. Principal coordinate analyses identified patterns in population structure based on genetic variation within the sequenced samples, with clear trends associated with increased diversity likely due to a higher number of infected individuals relative to the sampling dates. We demonstrate that genetic correlation analysis combined with SNVs analysis using wastewater sampling can provide a comprehensive snapshot of the SARS-CoV-2 genetic population structure circulating within a community, which might not be observed if relying solely on clinical cases.

High-throughput sequencing of SARS-CoV-2 in wastewater provides insights into circulating variants

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) likely emerged from a zoonotic spill-over event and has led to a global pandemic. The public health response has been predominantly informed by surveillance of symptomatic individuals and contact tracing, with quarantine, and other preventive measures have then been applied to mitigate further spread. Non-traditional methods of surveillance such as genomic epidemiology and wastewater-based epidemiology (WBE) have also been leveraged during this pandemic. Genomic epidemiology uses high-throughput sequencing of SARS-CoV-2 genomes to inform local and international transmission events, as well as the diversity of circulating variants. WBE uses wastewater to analyse community spread, as it is known that SARS-CoV-2 is shed through bodily excretions. Since both symptomatic and asymptomatic individuals contribute to wastewater inputs, we hypothesized that the resultant pooled sample of population-wide excreta can provide a more comprehensive picture of SARS-CoV-2 genomic diversity circulating in a community than clinical testing and sequencing alone. In this study, we analysed 91 wastewater samples from 11 states in the USA, where the majority of samples represent Maricopa County, Arizona (USA). With the objective of assessing the viral diversity at a population scale, we undertook a single-nucleotide variant (SNV) analysis on data from 52 samples with >90% SARS-CoV-2 genome coverage of sequence reads, and compared these SNVs with those detected in genomes sequenced from clinical patients. We identified 7973 SNVs, of which 548 were "novel" SNVs that had not yet been identified in the global clinical-derived data as of 17^th June 2020 (the day after our last wastewater sampling date). However, between 17^th of June 2020 and 20^th November 2020, almost half of the novel SNVs have since been detected in clinical-derived data. Using the combination of SNVs present in each sample, we identified the more probable lineages present in that sample and compared them to lineages observed in North America prior to our sampling dates. The wastewater-derived SARS-CoV-2 sequence data indicates there were more lineages circulating across the sampled communities than represented in the clinical-derived data. Principal coordinate analyses identified patterns in population structure based on genetic variation within the sequenced samples, with clear trends associated with increased diversity likely due to a higher number of infected individuals relative to the sampling dates. We demonstrate that genetic correlation analysis combined with SNVs analysis using wastewater sampling can provide a comprehensive snapshot of the SARS-CoV-2 genetic population structure circulating within a community, which might not be observed if relying solely on clinical cases.

A novel lineage of polyomaviruses identified in bark scorpions

Polyomaviruses are non-enveloped viruses with circular double-stranded DNA genomes (~4-7 kb). Initially identified in mammals, polyomaviruses have now been identified in birds and a few fish species. Although fragmentary polyomavirus-like sequences have been detected as apparent 'hitchhikers' in shotgun genomics datasets of various arthropods, the possible diversity of these viruses in invertebrates remains unclear. Scorpions are predatory arachnids that are among the oldest terrestrial animals. Using high-throughput sequencing and traditional molecular techniques we determine the genome sequences of eight novel polyomaviruses in scorpions (Centruroides sculpturatus) from the greater Phoenix area, Arizona, USA. Analysis of Centruroides transcriptomic datasets elucidated the splicing of the viral late gene array, which is more complex than that of vertebrate polyomaviruses. Phylogenetic analysis provides further evidence of co-divergence of polyomaviruses with their hosts, suggesting that at least one ancestral species of polyomaviruses was circulating amongst the primitive common ancestors of arthropods and chordates.

Genome characterization of parsley severe stunt-associated virus in Iran

Parsley severe stunt-associated virus (PSSaV) is a recently identified nanovirus first reported in Germany. During a survey for identification of nanoviruses infecting apiaceous plants in south-eastern Iran, PSSaV was identified and characterized using a combination of rolling circle amplification (RCA) and high-throughput sequencing. Parsley plant samples were collected from vegetable production farms in Kerman province. From two symptomatic samples (39Ba and 40Ba), seven PSSaV components (DNA-C, -S, -M, -R, -N, -U1 and -U2) with two phylogenetically distinct variants of DNA-R (R1 and R2) were identified. In common with the German isolate of PSSaV, no DNA-U4 component was identified. In addition, associated alphasatellite molecules were identified in samples 39Ba [n = 6] and 40Ba [n = 5]. Sequence analyses showed that concatenated component sequences of the two Iranian PSSaVs share 97.2% nucleotide identity with each other and 82% to the German isolate. The coat proteins (CPs) of the PSSaV Iranian sequences share 97.2% amino acid identity and ~ 84% identity with that of the German isolate. Sequence and phylogenetic analyses of a total of 11 recovered alphasatellites from the two samples can be classified into the genera Fabenesatellite [n = 2], Milvetsatellite [n = 1], Mivedwarsatellite [n = 2], Subclovsatellite [n = 2], Sophoyesatellite [n = 4] in the family Alphasatellitidae. Identification of PSSaV and other nanoviruses in wild and cultivated plants in Iran reveals that nanoviruses could be causing yield reduction in crops plants in this country.

Circular DNA viruses identified in short-finned pilot whale and orca tissue samples

Members of the Delphinidae family are widely distributed across the world's oceans. We used a viral metagenomic approach to identify viruses in orca (Orcinus orca) and short-finned pilot whale (Globicephala macrorhynchus) muscle, kidney, and liver samples from deceased animals. From orca tissue samples (muscle, kidney, and liver), we identified a novel polyomavirus (Polyomaviridae), three cressdnaviruses, and two genomoviruses (Genomoviridae). In the short-finned pilot whale we were able to identify one genomovirus in a kidney sample. The presence of unclassified cressdnavirus within two samples (muscle and kidney) of the same animal supports the possibility these viruses might be widespread within the animal. The orca polyomavirus identified here is the first of its species and is not closely related to the only other dolphin polyomavirus previously discovered. The identification and verification of these viruses expands the current knowledge of viruses that are associated with the Delphinidae family.

High-throughput sequencing of SARS-CoV-2 in wastewater provides insights into circulating variants

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged from a zoonotic spill-over event and has led to a global pandemic. The public health response has been predominantly informed by surveillance of symptomatic individuals and contact tracing, with quarantine, and other preventive measures have then been applied to mitigate further spread. Non-traditional methods of surveillance such as genomic epidemiology and wastewater-based epidemiology (WBE) have also been leveraged during this pandemic. Genomic epidemiology uses high-throughput sequencing of SARS-CoV-2 genomes to inform local and international transmission events, as well as the diversity of circulating variants. WBE uses wastewater to analyse community spread, as it is known that SARS-CoV-2 is shed through bodily excretions. Since both symptomatic and asymptomatic individuals contribute to wastewater inputs, we hypothesized that the resultant pooled sample of population-wide excreta can provide a more comprehensive picture of SARS-CoV-2 genomic diversity circulating in a community than clinical testing and sequencing alone. In this study, we analysed 91 wastewater samples from 11 states in the USA, where the majority of samples represent Maricopa County, Arizona (USA). With the objective of assessing the viral diversity at a population scale, we undertook a single-nucleotide variant (SNV) analysis on data from 52 samples with >90% SARS-CoV-2 genome coverage of sequence reads, and compared these SNVs with those detected in genomes sequenced from clinical patients. We identified 7973 SNVs, of which 5680 were novel SNVs that had not yet been identified in the global clinical-derived data as of 17th June 2020 (the day after our last wastewater sampling date). However, between 17th of June 2020 and 20th November 2020, almost half of the SNVs have since been detected in clinical-derived data. Using the combination of SNVs present in each sample, we identified the more probable lineages present in that sample and compared them to lineages observed in North America prior to our sampling dates. The wastewater-derived SARS-CoV-2 sequence data indicates there were more lineages circulating across the sampled communities than represented in the clinical-derived data. Principal coordinate analyses identified patterns in population structure based on genetic variation within the sequenced samples, with clear trends associated with increased diversity likely due to a higher number of infected individuals relative to the sampling dates. We demonstrate that genetic correlation analysis combined with SNVs analysis using wastewater sampling can provide a comprehensive snapshot of the SARS-CoV-2 genetic population structure circulating within a community, which might not be observed if relying solely on clinical cases.

Diverse cressdnaviruses and an anellovirus identified in the fecal samples of yellow-bellied marmots

Over that last decade, coupling multiple strand displacement approaches with high throughput sequencing have resulted in the identification of genomes of diverse groups of small circular DNA viruses. Using a similar approach but with recovery of complete genomes by PCR, we identified a diverse group of single-stranded viruses in yellow-bellied marmot (Marmota flaviventer) fecal samples. From 13 fecal samples we identified viruses in the family Genomoviridae (n = 7) and Anelloviridae (n = 1), and several others that ware part of the larger Cressdnaviricota phylum but not within established families (n = 19). There were also circular DNA molecules identified (n = 4) that appear to encode one viral-like gene and have genomes of <1545 nts. This study gives a snapshot of viruses associated with marmots based on fecal sampling.

Novel Circoviruses Detected in Feces of Sonoran Felids

Sonoran felids are threatened by drought and habitat fragmentation. Vector range expansion and anthropogenic factors such as habitat encroachment and climate change are altering viral evolutionary dynamics and exposure. However, little is known about the diversity of viruses present in these populations. Small felid populations with lower genetic diversity are likely to be most threatened with extinction by emerging diseases, as with other selective pressures, due to having less adaptive potential. We used a metagenomic approach to identify novel circoviruses, which may have a negative impact on the population viability, from confirmed bobcat (Lynx rufus) and puma (Puma concolor) scats collected in Sonora, Mexico. Given some circoviruses are known to cause disease in their hosts, such as porcine and avian circoviruses, we took a non-invasive approach using scat to identify circoviruses in free-roaming bobcats and puma. Three circovirus genomes were determined, and, based on the current species demarcation, they represent two novel species. Phylogenetic analyses reveal that one circovirus species is more closely related to rodent associated circoviruses and the other to bat associated circoviruses, sharing highest genome-wide pairwise identity of approximately 70% and 63%, respectively. At this time, it is unknown whether these scat-derived circoviruses infect felids, their prey, or another organism that might have had contact with the scat in the environment. Further studies should be conducted to elucidate the host of these viruses and assess health impacts in felids.

Transmission of the Bean-Associated Cytorhabdovirus by the Whitefly Bemisia tabaci MEAM1

The knowledge of genomic data of new plant viruses is increasing exponentially; however, some aspects of their biology, such as vectors and host range, remain mostly unknown. This information is crucial for the understanding of virus–plant interactions, control strategies, and mechanisms to prevent outbreaks. Typically, rhabdoviruses infect monocot and dicot plants and are vectored in nature by hemipteran sap-sucking insects, including aphids, leafhoppers, and planthoppers. However, several strains of a potentially whitefly-transmitted virus, papaya cytorhabdovirus, were recently described: (i) bean-associated cytorhabdovirus (BaCV) in Brazil, (ii) papaya virus E (PpVE) in Ecuador, and (iii) citrus-associated rhabdovirus (CiaRV) in China. Here, we examine the potential of the Bemisia tabaci Middle East-Asia Minor 1 (MEAM1) to transmit BaCV, its morphological and cytopathological characteristics, and assess the incidence of BaCV across bean producing areas in Brazil. Our results show that BaCV is efficiently transmitted, in experimental conditions, by B. tabaci MEAM1 to bean cultivars, and with lower efficiency to cowpea and soybean. Moreover, we detected BaCV RNA in viruliferous whiteflies but we were unable to visualize viral particles or viroplasm in the whitefly tissues. BaCV could not be singly isolated for pathogenicity tests, identification of the induced symptoms, and the transmission assay. BaCV was detected in five out of the seven states in Brazil included in our study, suggesting that it is widely distributed throughout bean producing areas in the country. This is the first report of a whitefly-transmitted rhabdovirus.

Diverse genomoviruses representing twenty-nine species identified associated with plants

Genomoviruses (family Genomoviridae) are circular single-stranded DNA viruses that have been mainly identified through metagenomics studies in a wide variety of samples from various environments. Here, we describe 98 genomes of genomoviruses found associated with members of 19 plant families from Australia, Brazil, France, South Africa and the USA. These 98 genomoviruses represent 29 species, 26 of which are new, in the genera Gemykolovirus (n = 37), Gemyduguivirus (n = 9), Gemygorvirus (n = 8), Gemykroznavirus (n = 6), Gemycircularvirus (n = 21) and Gemykibivirus (n = 17).

A Novel Divergent Geminivirus Identified in Asymptomatic New World Cactaceae Plants

Cactaceae comprise a diverse and iconic group of flowering plants which are almost exclusively indigenous to the New World. The wide variety of growth forms found amongst the cacti have led to the trafficking of many species throughout the world as ornamentals. Despite the evolution and physiological properties of these plants having been extensively studied, little research has focused on cactus-associated viral communities. While only single-stranded RNA viruses had ever been reported in cacti, here we report the discovery of cactus-infecting single-stranded DNA viruses. These viruses all apparently belong to a single divergent species of the family Geminiviridae and have been tentatively named Opuntia virus 1 (OpV1). A total of 79 apparently complete OpV1 genomes were recovered from 31 different cactus plants (belonging to 20 different cactus species from both the Cactoideae and Opuntioideae clades) and from nine cactus-feeding cochineal insects (Dactylopius sp.) sampled in the USA and Mexico. These 79 OpV1 genomes all share > 78.4% nucleotide identity with one another and < 64.9% identity with previously characterized geminiviruses. Collectively, the OpV1 genomes display evidence of frequent recombination, with some genomes displaying up to five recombinant regions. In one case, recombinant regions span ~40% of the genome. We demonstrate that an infectious clone of an OpV1 genome can replicate in Nicotiana benthamiana and Opuntia microdasys. In addition to expanding the inventory of viruses that are known to infect cacti, the OpV1 group is so distantly related to other known geminiviruses that it likely represents a new geminivirus genus. It remains to be determined whether, like its cactus hosts, its geographical distribution spans the globe.

Viruses representing two new genomovirus species identified in citrus from Tunisia

Using a high-throughput sequencing approach, we identified four genomoviruses (family Genomoviridae) associated with a sweet orange (Citrus sinensis) plant collected in Tunisia. The ssDNA genomes of these genomoviruses, which were amplified, cloned and Sanger sequenced, range in size from 2156 to 2191 nt. Three of these viruses share > 99% full-genome pairwise sequence identity and are referred to as citrus Tunisia genomovirus 1 (CTNGmV-1). The CTNGmV-1 isolates share < 62% genome-wide pairwise nucleotide sequence identity with other genomoviruses and belong to the genus Gemykolovirus. The genome of the fourth virus, which was called CTNGmV-2, shares < 68% nucleotide sequence identity with other genomoviruses and belongs to the genus Gemycircularvirus. Based on the species demarcation criteria for members of the family Genomoviridae, CTNGmV-1 and -2 would each represent a new species. Although found associated with Citrus sp. plants, it is likely that these viruses infect fungi or other organisms associated with the plants.

Virus Discovery in Desert Tortoise Fecal Samples: Novel Circular Single-Stranded DNA Viruses

The Sonoran Desert tortoise Gopherus morafkai is adapted to the desert, and plays an important ecological role in this environment. There is limited information on the viral diversity associated with tortoises (family Testudinidae), and to date no DNA virus has been identified associated with these animals. This study aimed to assess the diversity of DNA viruses associated with the Sonoran Desert tortoise by sampling their fecal matter. A viral metagenomics approach was used to identify the DNA viruses in fecal samples from wild Sonoran Desert tortoises in Arizona, USA. In total, 156 novel single-stranded DNA viruses were identified from 40 fecal samples. Those belonged to two known viral families, the Genomoviridae (n = 27) and Microviridae (n = 119). In addition, 10 genomes were recovered that belong to the unclassified group of circular-replication associated protein encoding single-stranded (CRESS) DNA virus and five circular molecules encoding viral-like proteins.

Single-Stranded DNA Viruses in Antarctic Cryoconite Holes

Antarctic cryoconite holes, or small melt-holes in the surfaces of glaciers, create habitable oases for isolated microbial communities with tightly linked microbial population structures. Viruses may influence the dynamics of polar microbial communities, but the viromes of the Antarctic cryoconite holes have yet to be characterized. We characterize single-stranded DNA (ssDNA) viruses from three cryoconite holes in the Taylor Valley, Antarctica, using metagenomics. Half of the assembled metagenomes cluster with those in the viral family Microviridae (n = 7), and the rest with unclassified circular replication associated protein (Rep)-encoding single-stranded (CRESS) DNA viruses (n = 7). An additional 18 virus-like circular molecules encoding either a Rep, a capsid protein gene, or other unidentified but viral-like open reading frames were identified. The samples from which the genomes were identified show a strong gradient in microbial diversity and abundances, and the number of viral genomes detected in each sample mirror that gradient. Additionally, one of the CRESS genomes assembled here shares ~90% genome-wide pairwise identity with a virus identified from a freshwater pond on the McMurdo Ice Shelf (Antarctica). Otherwise, the similarity of these viruses to those previously identified is relatively low. Together, these patterns are consistent with the presence of a unique regional virome present in fresh water host populations of the McMurdo Dry Valley region.

Single Stranded DNA Viruses Associated with Capybara Faeces Sampled in Brazil

Capybaras (Hydrochoerus hydrochaeris), the world's largest rodents, are distributed throughout South America. These wild herbivores are commonly found near water bodies and are well adapted to rural and urban areas. There is limited information on the viruses circulating through capybaras. This study aimed to expand the knowledge on the viral diversity associated with capybaras by sampling their faeces. Using a viral metagenomics approach, we identified diverse single-stranded DNA viruses in the capybara faeces sampled in the Distrito Federal, Brazil. A total of 148 complete genomes of viruses in the Microviridae family were identified. In addition, 14 genomoviruses (family Genomoviridae), a novel cyclovirus (family Circoviridae), and a smacovirus (family Smacoviridae) were identified. Also, 37 diverse viruses that cannot be assigned to known families and more broadly referred to as unclassified circular replication associated protein encoding single-stranded (CRESS) DNA viruses were identified. This study provides a snapshot of the viral diversity associated with capybaras that may be infectious to these animals or associated with their microbiota or diet.

Unravelling the Single-Stranded DNA Virome of the New Zealand Blackfly

Over the last decade, arthropods have been shown to harbour a rich diversity of viruses. Through viral metagenomics a large diversity of single-stranded (ss) DNA viruses have been identified. Here we examine the ssDNA virome of the hematophagous New Zealand blackfly using viral metagenomics. Our investigation reveals a plethora of novel ssDNA viral genomes, some of which cluster in the viral families Genomoviridae (n = 9), Circoviridae (n = 1), and Microviridae (n = 108), others in putative families that, at present, remain unclassified (n = 20) and one DNA molecule that only encodes a replication associated protein. Among these novel viruses, two putative multi-component virus genomes were recovered, and these are most closely related to a Tongan flying fox faeces-associated multi-component virus. Given that the only other known multi-component circular replication-associated (Rep) protein encoding single-stranded (CRESS) DNA viruses infecting plants are in the families Geminiviridae (members of the genus Begomovirus) and Nanoviridae, it appears these are likely a new multi-component virus group which may be associated with animals. This study reiterates the diversity of ssDNA viruses in nature and in particular with the New Zealand blackflies.

Diverse single-stranded DNA viruses associated with honey bees (Apis mellifera)

Honey bees (Apis mellifera) research has increased in light of their progressive global decline over the last decade and the important role they play in pollination. One expanding area of honey bee research is analysis of their microbial community including viruses. Several RNA viruses have been characterized but little is known about DNA viruses associated with bees. Here, using a metagenomics based approach, we reveal the presence of a broad range of novel single-stranded DNA viruses from the hemolymph and brain of nurse and forager (worker divisions of labour) bees belonging to two honey bees subspecies, Italian (Apis mellifera linguistica) and New World Carniolan (Apis mellifera carnica). Genomes of 100 diverse viruses were identified, designated into three groupings; genomoviruses (family Genomoviridae) (n = 4), unclassified replication associated protein encoding single-stranded DNA viruses (n = 28), and microviruses (family Microviridae; subfamily Gokushovirinae) (n = 70). Amongst the viruses identified, it appears that nurses harbour a higher diversity of these viruses comparative to the foragers. Between subspecies, the most striking outcome was the extremely high number of diverse microviruses identified in the Italian bees comparative to the New World Carniolan, likely indicating an association to the diversity of the bacterial community associated with these subspecies.

Diverse papillomaviruses identified in Weddell seals

Papillomaviridae is a diverse family of circular, double-stranded DNA (dsDNA) viruses that infect a broad range of mammalian, avian and fish hosts. While papillomaviruses have been characterized most extensively in humans, the study of non-human papillomaviruses has contributed greatly to our understanding of their pathogenicity and evolution. Using high-throughput sequencing approaches, we identified 7 novel papillomaviruses from vaginal swabs collected from 81 adult female Weddell seals (Leptonychotes weddellii) in the Ross Sea of Antarctica between 2014-2017. These seven papillomavirus genomes were amplified from seven individual seals, and six of the seven genomes represented novel species with distinct evolutionary lineages. This highlights the diversity of papillomaviruses among the relatively small number of Weddell seal samples tested. Viruses associated with large vertebrates are poorly studied in Antarctica, and this study adds information about papillomaviruses associated with Weddell seals and contributes to our understanding of the evolutionary history of papillomaviruses.

Discovery of the first maize-infecting mastrevirus in the Americas using a vector-enabled metagenomics approach

The genus Mastrevirus (family Geminiviridae) is composed of single-stranded DNA viruses that infect mono- and dicotyledonous plants and are transmitted by leafhoppers. In South America, there have been only two previous reports of mastreviruses, both identified in sweet potatoes (from Peru and Uruguay). As part of a general viral surveillance program, we used a vector-enabled metagenomics (VEM) approach and sampled leafhoppers (Dalbulus maidis) in Itumbiara (State of Goiás), Brazil. High-throughput sequencing of viral DNA purified from the leafhopper sample revealed mastrevirus-like contigs. Using a set of abutting primers, a 2746-nt circular genome was recovered. The circular genome has a typical mastrevirus genome organization and shares <63% pairwise identity with other mastrevirus isolates from around the world. Therefore, the new mastrevirus was tentatively named "maize striate mosaic virus". Seventeen maize leaf samples were collected in the same field as the leafhoppers, and ten samples were found to be positive for this mastrevirus. Furthermore, the ten genomes recovered from the maize samples share >99% pairwise identity with the one from the leafhopper. This is the first report of a maize-infecting mastrevirus in the Americas, the first identified in a non-vegetatively propagated mastrevirus host in South America, and the first mastrevirus to be identified in Brazil.

A novel geminivirus identified in tomato and cleome plants sampled in Brazil

Viruses in the family Geminiviridae have single-stranded DNA genomes encapsulated in geminate icosahedral particles. High throughput sequencing (HTS) for metagenomic approaches are being extensively used for the identification of known and novel viruses. Using a HTS approach, we identified a novel geminivirus in a tomato (Solanum lycopersicum) sample and a Cleome sp. sample collected in the midwest region of Brazil. The genomes from the two samples share 99.96% identity and ∼61-63% to genomes in the genus Capulavirus. The novel virus has been tentatively named tomato associated geminivirus 1 (TaGV1). No visual symptoms were observed in the field tomato plant or in the inoculated Nicotiana benthamiana where the virus established a systemic infection. The replication associated protein of TaGV1 is most similar to that encoded by capulaviruses (sharing 62-70% identity), whereas the CP is most similar to that of tomato pseudo curly top virus (sharing ∼31% identity). In the TaGV1 positive Cleome sp. sample, begomovirus DNA A and B components were also detected sharing 96% and 90% sequence identity to cleome leaf crumple virus DNA A and B components, respectively. Using a HTS approach, we identified TaGV1 in tomato and Cleome sp. samples and this is the first report of a geminivirus that is non-begomovirus in Brazil.

Molecular and biological characterization of corchorus mottle virus, a new begomovirus from Brazil

A begomovirus infecting Orinoco jute (Corchorus hirtus) from Brazil was characterized. Molecular analysis revealed a bipartite genomic organization, which is typical of the New World begomoviruses. Sequence analysis and phylogenetic data showed that both genomic components have the closest relationship with abutilon mosaic Brazil virus, with an identity of 87.3 % for DNA-A, indicating that this virus is a member of a new begomovirus species for which the name "Corchorus mottle virus" (CoMoV) is proposed. Sida rhombifolia plants inoculated by biolistics with an infectious clone of CoMoV showed systemic vein chlorosis, mottling and leaf deformation symptoms, while Nicotiana benthamiana and tomato plants had symptomless infection. CoMoV is the first corchorus-infecting begomovirus reported in Brazil.

First Report of Sida micrantha mosaic virus in Phaseolus vulgaris in Brazil

Snap and common beans (Phaseolus vulgaris L.) are severely affected by Bean golden mosaic virus (BGMV) infection, so far the only begomovirus reported on these crops in Brazil (1). Samples of snap and common beans colonized by the whitefly Bemisia tabaci biotype B and displaying golden mosaic, chlorotic spots, and leaf distortion were collected in three production regions in Goiás State (Goianápolis, Luziânia, and Itaberaí) between 2003 and 2007. Total DNA extracted from leaf samples was used as template in PCR assays using universal primers targeting conserved regions of the DNA-A and DNA-B genomes (3). Begomovirus-specific amplicons were observed only with DNA template from symptomatic plants. Two single amplicons were observed for both genomic segments, indicating the presence of bipartite species in all samples. Sequence analysis of four isolates (named as GO-176, GO-260, GO-354, and GO-368) obtained from common bean samples indicated identity levels of approximately 95% with the DNA-A segment of BGMV (GenBank Accession No. FJ665283). However, the complete DNA-A sequence (GenBank Accession No. HM357459.1) of the GO-060 isolate (from a symptomatic snap bean plant collected in Goianápolis) displayed 76% identity with BGMV (GenBank Accession No. FJ665283) and 95% identity with the DNA-A of a Sida micrantha mosaic virus (SimMV) isolate (GenBank Accession No. EU908733.1) reported to be infecting okra (Abelmoschus esculentus L.) and 94.8% with a SimMV isolate reported to be infecting soybean (GenBank Accession No. FJ686693) in Brazil (2). Koch's postulates were fulfilled for the isolate GO-060 by inoculating a set of soybean and bean accessions via a biolistic approach. The ratio of positive PCR amplicons per total of inoculated plants were 15 of 16 for snap bean cv. Trepador, 9 of 10 for snap bean cv. Fartura, 18 of 24 for common bean cv. Olate Pinto, and 19 of 25 for common bean cv. Carioca. The isolate was also able to infect eight of nine soybean 'Doko' plants. Sequence analysis using symptomatic leaf samples (15 days after inoculation) confirmed SimMV as the causal agent. To our knowledge, this is the first report of a SimMV isolate infecting P. vulgaris. This virus is apparently fast expanding its host range from Malvaceae to Solanaceae species and leguminous hosts after the introduction of B. tabaci biotype B (2). More extensive surveys are necessary to access the current epidemiological importance of SimMV in both snap and common beans in Brazil. References: (1) J.C. Faria and D. P. Maxwell. Phytopathology 89:262, 1999. (2) F. R. Fernandes et al. Arch. Virol. 154:1567, 2009. (3) M. R. Rojas et al. Plant Dis. 77:340, 1993.