Cenote-Taker 2 democratizes virus discovery and sequence annotation

Opportunities and barriers in omics-based biomarker discovery for steatotic liver diseases

The rising prevalence of liver diseases related to obesity and excessive use of alcohol is fuelling an increasing demand for accurate biomarkers aimed at community screening, diagnosis of steatohepatitis and significant fibrosis, monitoring, prognostication and prediction of treatment efficacy. Breakthroughs in omics methodologies and the power of bioinformatics have created an excellent opportunity to apply technological advances to clinical needs, for instance in the development of precision biomarkers for personalised medicine. Via omics technologies, biological processes from the genes to circulating protein, as well as the microbiome - including bacteria, viruses and fungi, can be investigated on an axis. However, there are important barriers to omics-based biomarker discovery and validation, including the use of semi-quantitative measurements from untargeted platforms, which may exhibit high analytical, inter- and intra-individual variance. Standardising methods and the need to validate them across diverse populations presents a challenge, partly due to disease complexity and the dynamic nature of biomarker expression at different disease stages. Lack of validity causes lost opportunities when studies fail to provide the knowledge needed for regulatory approvals, all of which contributes to a delayed translation of these discoveries into clinical practice. While no omics-based biomarkers have matured to clinical implementation, the extent of data generated has enabled the hypothesis-free discovery of a plethora of candidate biomarkers that warrant further validation. To explore the many opportunities of omics technologies, hepatologists need detailed knowledge of commonalities and differences between the various omics layers, and both the barriers to and advantages of these approaches.

Genomes of cressdnaviruses, microviruses, a gyrovirus, and a caudovirus identified in the fecal samples of a mallard and double-crested cormorant

Mallards and double-crested cormorants have a broad distribution across North America. In the fecal sample from two individual mallard and double-crested cormorant, we determined the genomes of a caudovirus, microviruses (n = 6), cressdnaviruses (n = 35), and a gyrovirus (chicken anemia virus, CAV). Here, we report double-crested cormorant as a CAV host.

Diverse Circular DNA Viral Communities in Blood, Oral, and Fecal Samples of Captive Lemurs

Few studies have addressed viral diversity in lemurs despite their unique evolutionary history on the island of Madagascar and high risk of extinction. Further, while a large number of studies on animal viromes focus on fecal samples, understanding viral diversity across multiple sample types and seasons can reveal complex viral community structures within and across species. Groups of captive lemurs at the Duke Lemur Center (Durham, NC, USA), a conservation and research center, provide an opportunity to build foundational knowledge on lemur-associated viromes. We sampled individuals from seven lemur species, i.e., collared lemur (Eulemur collaris), crowned lemur (Eulemur coronatus), blue-eyed black lemur (Eulemur flavifrons), ring-tailed lemur (Lemur catta), Coquerel's sifaka (Propithecus coquereli), black-and-white ruffed lemur (Varecia variegata variegata), and red ruffed lemur (Varecia rubra), across two lemur families (Lemuridae, Indriidae). Fecal, blood, and saliva samples were collected from Coquerel's sifaka and black-and-white ruffed lemur individuals across two sampling seasons to diversify virome biogeography and temporal sampling. Using viral metagenomic workflows, the complete genomes of anelloviruses (n = 4), cressdnaviruses (n = 47), caudoviruses (n = 15), inoviruses (n = 34), and microviruses (n = 537) were determined from lemur blood, feces, and saliva. Many virus genomes, especially bacteriophages, identified in this study were present across multiple lemur species. Overall, the work presented here uses a viral metagenomics approach to investigate viral communities inhabiting the blood, oral cavity, and feces of healthy captive lemurs.

Metagenome-assembled microvirus and cressdnavirus genomes from fecal samples of house mice (Mus musculus)

House mice, Mus musculus, are highly adapted to anthropogenic spaces. Fecal samples were collected from house mice entering primate enclosure areas at the Duke Lemur Center (Durham, NC, USA). We identified 14 cressdnavirus and 59 microvirus genomes in these mouse feces.

Longitudinal alterations in the urinary virome of kidney transplant recipients are influenced by BK viremia and patient sex

Little is known about the urinary virome and how it interacts with the host, particularly in renal transplant diseases. Using metagenomic sequencing, we characterized the urinary virome of 23 kidney transplant recipients longitudinally (11 BKV+ patients and 12 BKV- patients). We applied linear mixed effects models, PERMANOVA, k-means clustering, and MaAsLin2 algorithms to determine virome signatures associated with post-transplant time, BK viremia status, and patient sex. We found that the richness and alpha diversity of urinary virome were significantly different in renal transplant recipients with BKV+ over time in comparison to BKV- (richness P = 0.012, alpha P < 0.0001). Female BKV- patients had significantly higher virome richness than males (P = 0.0063). Virome beta diversity was significantly different between patients by BKV status (P < 0.001). Additionally, we identified underlying interactions between patient sex and BKV status, in terms of virome beta diversity (P = 0.008). BK polyomavirus infections were primarily of subtypes IA, IB1, and IB2. The non-BK dominant samples clustered into six urinary virome community states. BKV- samples had more anelloviruses than BKV+ samples though this difference was not statistically significant. Lastly, we identified specific viruses, associated with BKV+ and time in our samples. Our results indicate that dynamic alterations in the urinary virome over the post-transplant period in kidney transplant recipients can be shaped by BK viremia and patient sex. These findings advance our fundamental understanding of the urinary virome and support a new line of investigation in renal disease and transplantation.

IMPORTANCE

The urinary microbiome is increasingly implicated in renal health and disease. While most research focuses on bacteria communities of the microbiome, factors that influence the urinary virome are not understood. Here, we investigated the urinary virome of 23 adult kidney transplant recipients longitudinally over 14 weeks post-transplant. We show that alterations in the urinary virome are associated with kidney transplant recipients with BK polyomavirus viremia that can lead to BK nephropathy and allograft rejection. By modeling the temporal dynamics post-transplant, we delineated specific profiles of the urinary virome associated with patient sex and urinary community states. These findings reveal fundamental aspects of the urinary virome that can be leveraged to better manage kidney diseases.

Sequencing-based analysis of microbiomes

Microbiomes occupy a range of niches and, in addition to having diverse compositions, they have varied functional roles that have an impact on agriculture, environmental sciences, and human health and disease. The study of microbiomes has been facilitated by recent technological and analytical advances, such as cheaper and higher-throughput DNA and RNA sequencing, improved long-read sequencing and innovative computational analysis methods. These advances are providing a deeper understanding of microbiomes at the genomic, transcriptional and translational level, generating insights into their function and composition at resolutions beyond the species level.

Identification and classification of the genomes of novel microviruses in poultry slaughterhouse

Microviridae is a family of phages with circular ssDNA genomes and they are widely found in various environments and organisms. In this study, virome techniques were employed to explore potential members of Microviridae in a poultry slaughterhouse, leading to the identification of 98 novel and complete microvirus genomes. Using a similarity clustering network classification approach, these viruses were found to belong to at least 6 new subfamilies within Microviridae and 3 higher-level taxonomic units. Genome size, GC content and genome structure of these new taxa showed evident regularities, validating the rationality of our classification method. Our method can divide microviruses into about 45 additional detailed clusters, which may serve as a new standard for classifying Microviridae members. Furthermore, by addressing the scarcity of host information for microviruses, the current study significantly broadened their host range and discovered over 20 possible new hosts, including important pathogenic bacteria such as Helicobacter pylori and Vibrio cholerae, as well as different taxa demonstrated different host specificities. The findings of this study effectively expand the diversity of the Microviridae family, providing new insights for their classification and identification. Additionally, it offers a novel perspective for monitoring and controlling pathogenic microorganisms in poultry slaughterhouse environments.

Two Lineages of Papillomaviruses Identified from Caracals (Caracal caracal) in South Africa

Papillomaviruses (PV) infect epithelial cells and can cause hyperplastic or neoplastic lesions. In felids, most described PVs are from domestic cats (Felis catus; n = 7 types), with one type identified in each of the five wild felid species studied to date (Panthera uncia, Puma concolor, Leopardus wiedii, Panthera leo persica and Lynx rufus). PVs from domestic cats are highly diverse and are currently classified into three genera (Lambdapapillomavirus, Dyothetapapillomavirus, and Taupapillomavirus), whereas those from wild felids, although diverse, are all classified into the Lambdapapillomavirus genus. In this study, we used a metagenomic approach to identify ten novel PV genomes from rectal swabs of five deceased caracals (Caracal caracal) living in the greater Cape Town area, South Africa. These are the first PVs to be described from caracals, and represent six new PV types, i.e., Caracal caracal papillomavirus (CcarPV) 1-6. These CcarPV fall into two phylogenetically distinct genera: Lambdapapillomavirus, and Treisetapapillomavirus. Two or more PV types were identified in a single individual for three of the five caracals, and four caracals shared at least one of the same PV types with another caracal. This study broadens our understanding of wild felid PVs and provides evidence that there may be several wild felid PV lineages.

Gut virome and microbiome dynamics before and after SARS-CoV-2 infection in women living with HIV and their infants

Microbiome perturbations can have long-term effects on health. The dynamics of the gut microbiome and virome in women living with HIV (WLHIV) and their newborn infants is poorly understood. Here, we performed metagenomic sequencing analyses on longitudinal stool samples including 23 mothers (13 WLHIV, 10 HIV-negative) and 12 infants that experienced SARS-CoV-2 infection with mild disease, as well as 40 mothers (18 WLHIV, 22 HIV-negative) and 60 infants that remained SARS-CoV-2 seronegative throughout the study follow-up. Regardless of HIV or SARS-CoV-2 status, maternal bacterial and viral profiles were distinct from infants. Using linear mixed effects models, we showed that while the microbiome alpha diversity trajectory was not significantly different between SARS-CoV-2 seropositive and seronegative women. However, seropositive women's positive trajectory while uninfected was abruptly reversed after SARS-CoV-2 infection (p = 0.015). However, gut virome signatures of women were not associated with SARS-CoV-2. Alterations in infant microbiome and virome diversities were generally not impacted by SARS-CoV-2 but were rather driven by development. We did not find statistically significant interactions between HIV and SARS-CoV-2 on the gut microbiome and virome. Overall, our study provides insights into the complex interplay between maternal and infant bacterial microbiome, virome, and the influence of SARS-CoV-2 and HIV status.

Two Novel Geminiviruses Identified in Bees (Apis mellifera and Nomia sp.)

Members of the Geminviridae family are circular single-stranded DNA plant-infecting viruses, some of which impact global food production. Geminiviruses are vectored by sap-feeding insects such as leafhoppers, treehoppers, aphids, and whiteflies. Additionally, geminivirus sequences have also been identified in other insects such as dragonflies, mosquitoes, and stingless bees. As part of a viral metagenomics study on honeybees and solitary bees (Nomia sp.), two geminivirus genomes were identified. These represent a novel citlodavirus (from honeybees collected from Westmoreland, Jamaica) and a mastrevirus-like genome (from a solitary bee collected from Tempe, Arizona, USA). The novel honeybee-derived citlodavirus genome shares ~61 to 69% genome-wide nucleotide pairwise identity with other citlodavirus genome sequences and is most closely related to the passion fruit chlorotic mottle virus identified in Brazil. Whereas the novel solitary bee-derived mastrevirus-like genome shares ~55 to 61% genome-wide nucleotide identity with other mastreviruses and is most closely related to tobacco yellow dwarf virus identified in Australia, based on pairwise identity scores of the full genome, replication-associated protein, and capsid protein sequences. Previously, two geminiviruses in the Begomovirus genus were identified in samples of stingless bee (Trigona spp.) samples. Here, we identify viruses that represent two new species of geminiviruses from a honeybee and a solitary bee, which continues to demonstrate that plant pollinators can be utilized for the identification of plant-infecting DNA viruses in ecosystems.

The virome of the kitome: small circular virus-like genomes in laboratory reagents

In the course of studying the virome of protozoan parasites, we identified small circular genomes resembling viruses, which turned out to be contaminants from an RNA purification kit. We report their sequences here so others can detect possible contamination in their samples by aligning them to these targets.

Phage-specific immunity impairs efficacy of bacteriophage targeting Vancomycin Resistant Enterococcus in a murine model

Bacteriophage therapy is a promising approach to address antimicrobial infections though questions remain regarding the impact of the immune response on clinical effectiveness. Here, we develop a mouse model to assess phage treatment using a cocktail of five phages from the Myoviridae and Siphoviridae families that target Vancomycin-Resistant Enterococcus gut colonization. Phage treatment significantly reduces fecal bacterial loads of Vancomycin-Resistant Enterococcus. We also characterize immune responses elicited following administration of the phage cocktail. While minimal innate responses are observed after phage administration, two rounds of treatment induces phage-specific neutralizing antibodies and accelerate phage clearance from tissues. Interestingly, the myophages in our cocktail induce a more robust neutralizing antibody response than the siphophages. This anti-phage immunity reduces the effectiveness of the phage cocktail in our murine model. Collectively, this study shows phage-specific immune responses may be an important consideration in the development of phage cocktails for therapeutic use.

Correlation between the gut microbiome and neurodegenerative diseases: a review of metagenomics evidence

A growing body of evidence suggests that the gut microbiota contributes to the development of neurodegenerative diseases via the microbiota-gut-brain axis. As a contributing factor, microbiota dysbiosis always occurs in pathological changes of neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. High-throughput sequencing technology has helped to reveal that the bidirectional communication between the central nervous system and the enteric nervous system is facilitated by the microbiota's diverse microorganisms, and for both neuroimmune and neuroendocrine systems. Here, we summarize the bioinformatics analysis and wet-biology validation for the gut metagenomics in neurodegenerative diseases, with an emphasis on multi-omics studies and the gut virome. The pathogen-associated signaling biomarkers for identifying brain disorders and potential therapeutic targets are also elucidated. Finally, we discuss the role of diet, prebiotics, probiotics, postbiotics and exercise interventions in remodeling the microbiome and reducing the symptoms of neurodegenerative diseases.

Widespread Horizontal Gene Transfer Among Animal Viruses

The initial objective of this study was to shed light on the evolution of small DNA tumor viruses by analyzing de novo assemblies of publicly available deep sequencing datasets. The survey generated a searchable database of contig snapshots representing more than 100,000 Sequence Read Archive records. Using modern structure-aware search tools, we iteratively broadened the search to include an increasingly wide range of other virus families. The analysis revealed a surprisingly diverse range of chimeras involving different virus groups. In some instances, genes resembling known DNA-replication modules or known virion protein operons were paired with unrecognizable sequences that structural predictions suggest may represent previously unknown replicases and novel virion architectures. Discrete clades of an emerging group called adintoviruses were discovered in datasets representing humans and other primates. As a proof of concept, we show that the contig database is also useful for discovering RNA viruses and candidate archaeal phages. The ancillary searches revealed additional examples of chimerization between different virus groups. The observations support a gene-centric taxonomic framework that should be useful for future virus-hunting efforts.

New duck papillomavirus type identified in a mallard in Missouri, USA

Papillomaviruses are small circular DNA viruses that infect epithelial and mucosal cells and have co-evolved with their hosts. Some papillomaviruses in mammals are well studied (especially those associated with disease). However, there is limited information on papillomaviruses associated with avian hosts. From a cloacal swab sample of a mallard (Anas platyrhynchos) sampled in Missouri, USA (6 Jan 2023), we identified a papillomavirus (7839 nt) that shares ~68% genome-wide nucleotide sequence identity with Anas platyrhynchos papillomavirus 1 (AplaPV1) from a mallard sampled in Newfoundland (Canada) and ~40% with AplaPV2 from a mallard sampled in Minnesota (USA) with mesenchymal dermal tumors. The papillomavirus we identified shares 73.6% nucleotide sequence identity in the L1 gene with that of AplaPV1 and thus represents a new AplaPV type (AplaPV3). The genome sequence of AplaPV3 shares >97% identity with three partial PV genome sequences (1316, 1997, and 4241 nt) identified in a mallard in India, indicating that that virus was also AplaPV3.

Benchmarking informatics approaches for virus discovery: caution is needed when combining in silico identification methods

Understanding the ecological impacts of viruses on natural and engineered ecosystems relies on the accurate identification of viral sequences from community sequencing data. To maximize viral recovery from metagenomes, researchers frequently combine viral identification tools. However, the effectiveness of this strategy is unknown. Here, we benchmarked combinations of six widely used informatics tools for viral identification and analysis (VirSorter, VirSorter2, VIBRANT, DeepVirFinder, CheckV, and Kaiju), called "rulesets." Rulesets were tested against mock metagenomes composed of taxonomically diverse sequence types and diverse aquatic metagenomes to assess the effects of the degree of viral enrichment and habitat on tool performance. We found that six rulesets achieved equivalent accuracy [Matthews Correlation Coefficient (MCC) = 0.77, Padj ≥ 0.05]. Each contained VirSorter2, and five used our "tuning removal" rule designed to remove non-viral contamination. While DeepVirFinder, VIBRANT, and VirSorter were each found once in these high-accuracy rulesets, they were not found in combination with each other: combining tools does not lead to optimal performance. Our validation suggests that the MCC plateau at 0.77 is partly caused by inaccurate labeling within reference sequence databases. In aquatic metagenomes, our highest MCC ruleset identified more viral sequences in virus-enriched (44%-46%) than in cellular metagenomes (7%-19%). While improved algorithms may lead to more accurate viral identification tools, this should be done in tandem with careful curation of sequence databases. We recommend using the VirSorter2 ruleset and our empirically derived tuning removal rule. Our analysis provides insight into methods for in silico viral identification and will enable more robust viral identification from metagenomic data sets.

IMPORTANCE

The identification of viruses from environmental metagenomes using informatics tools has offered critical insights in microbial ecology. However, it remains difficult for researchers to know which tools optimize viral recovery for their specific study. In an attempt to recover more viruses, studies are increasingly combining the outputs from multiple tools without validating this approach. After benchmarking combinations of six viral identification tools against mock metagenomes and environmental samples, we found that these tools should only be combined cautiously. Two to four tool combinations maximized viral recovery and minimized non-viral contamination compared with either the single-tool or the five- to six-tool ones. By providing a rigorous overview of the behavior of in silico viral identification strategies and a pipeline to replicate our process, our findings guide the use of existing viral identification tools and offer a blueprint for feature engineering of new tools that will lead to higher-confidence viral discovery in microbiome studies.

Genomic surveillance of Clostridioides difficile transmission and virulence in a healthcare setting

Clostridioides difficile infection (CDI) is a major cause of healthcare-associated diarrhea, despite the widespread implementation of contact precautions for patients with CDI. Here, we investigate strain contamination in a hospital setting and the genomic determinants of disease outcomes. Across two wards over 6 months, we selectively cultured C. difficile from patients (n = 384) and their environments. Whole-genome sequencing (WGS) of 146 isolates revealed that most C. difficile isolates were from clade 1 (131/146, 89.7%), while only one isolate of the hypervirulent ST1 was recovered. Of culture-positive admissions (n = 79), 19 (24%) patients were colonized with toxigenic C. difficile on admission to the hospital. We defined 25 strain networks at ≤2 core gene single nucleotide polymorphisms; two of these networks contain strains from different patients. Strain networks were temporally linked (P < 0.0001). To understand the genomic correlates of the disease, we conducted WGS on an additional cohort of C. difficile (n = 102 isolates) from the same hospital and confirmed that clade 1 isolates are responsible for most CDI cases. We found that while toxigenic C. difficile isolates are associated with the presence of cdtR, nontoxigenic isolates have an increased abundance of prophages. Our pangenomic analysis of clade 1 isolates suggests that while toxin genes (tcdABER and cdtR) were associated with CDI symptoms, they are dispensable for patient colonization. These data indicate that toxigenic and nontoxigenic C. difficile contamination persist in a hospital setting and highlight further investigation into how accessory genomic repertoires contribute to C. difficile colonization and disease.

IMPORTANCE

Clostridioides difficile infection remains a leading cause of hospital-associated diarrhea, despite increased antibiotic stewardship and transmission prevention strategies. This suggests a changing genomic landscape of C. difficile. Our study provides insight into the nature of prevalent C. difficile strains in a hospital setting and transmission patterns among carriers. Longitudinal sampling of surfaces and patient stool revealed that both toxigenic and nontoxigenic strains of C. difficile clade 1 dominate these two wards. Moreover, quantification of transmission in carriers of these clade 1 isolates underscores the need to revisit infection prevention measures in this patient group. We identified unique genetic signatures associated with virulence in this clade. Our data highlight the complexities of preventing transmission of this pathogen in a hospital setting and the need to investigate the mechanisms of in vivo persistence and virulence of prevalent lineages in the host gut microbiome.

Complete genome sequence of a novel papillomavirus in Siamese fighting fish (Betta splendens) from Trinidad and Tobago

Here, we announce the complete genome of a previously undescribed papillomavirus from a betta fish, Betta splendens. The genome is 5,671 bp with a GC content of 38.2%. Variants were detected in public databases. This genome is most similar to papillomaviruses that infect sea bass (52.9 % nucleotide identity).

Forty-nine metagenomic-assembled genomes from an aquatic virome expand Caudoviricetes by 45 potential new families and the newly uncovered Gossevirus of Bamfordvirae

Twenty complete genomes (29-63 kb) and 29 genomes with an estimated completeness of over 90 % (30-90 kb) were identified for novel dsDNA viruses in the Yangshan Harbor metavirome. These newly discovered viruses contribute to the expansion of viral taxonomy by introducing 46 potential new families. Except for one virus, all others belong to the class Caudoviricetes. The exception is a novel member of the recently characterized viral group known as Gossevirus. Fifteen viruses were predicted to be temperate. The predicted hosts for the viruses appear to be involved in various aspects of the nitrogen cycle, including nitrogen fixation, oxidation and denitrification. Two viruses were identified to have a host of Flavobacterium and Tepidimonas fonticaldi, respectively, by matching CRISPR spacers with viral protospacers. Our findings provide an overview for characterizing and identifying specific viruses from Yangshan Harbor. The Gossevirus-like virus uncovered emphasizes the need for further comprehensive isolation and investigation of polinton-like viruses.

Five Species of Wild Freshwater Sport Fish in Wisconsin, USA, Reveal Highly Diverse Viromes

Studies of marine fish have revealed distant relatives of viruses important to global fish and animal health, but few such studies exist for freshwater fish. To investigate whether freshwater fish also host such viruses, we characterized the viromes of five wild species of freshwater fish in Wisconsin, USA: bluegill (Lepomis macrochirus), brown trout (Salmo trutta), lake sturgeon (Acipenser fulvescens), northern pike (Esox lucius), and walleye (Sander vitreus). We analyzed 103 blood serum samples collected during a state-wide survey from 2016 to 2020 and used a metagenomic approach for virus detection to identify known and previously uncharacterized virus sequences. We then characterized viruses phylogenetically and quantified prevalence, richness, and relative abundance for each virus. Within these viromes, we identified 19 viruses from 11 viral families: Amnoonviridae, Circoviridae, Coronaviridae, Hepadnaviridae, Peribunyaviridae, Picobirnaviridae, Picornaviridae, Matonaviridae, Narnaviridae, Nudnaviridae, and Spinareoviridae, 17 of which were previously undescribed. Among these viruses was the first fish-associated coronavirus from the Gammacoronavirus genus, which was present in 11/15 (73%) of S. vitreus. These results demonstrate that, similar to marine fish, freshwater fish also harbor diverse relatives of viruses important to the health of fish and other animals, although it currently remains unknown what effect, if any, the viruses we identified may have on fish health.

Large language models improve annotation of prokaryotic viral proteins

Viral genomes are poorly annotated in metagenomic samples, representing an obstacle to understanding viral diversity and function. Current annotation approaches rely on alignment-based sequence homology methods, which are limited by the paucity of characterized viral proteins and divergence among viral sequences. Here we show that protein language models can capture prokaryotic viral protein function, enabling new portions of viral sequence space to be assigned biologically meaningful labels. When applied to global ocean virome data, our classifier expanded the annotated fraction of viral protein families by 29%. Among previously unannotated sequences, we highlight the identification of an integrase defining a mobile element in marine picocyanobacteria and a capsid protein that anchors globally widespread viral elements. Furthermore, improved high-level functional annotation provides a means to characterize similarities in genomic organization among diverse viral sequences. Protein language models thus enhance remote homology detection of viral proteins, serving as a useful complement to existing approaches.

Hiding in plain sight: The discovery of complete genomes of 11 hypothetical spindle-shaped viruses that putatively infect mesophilic ammonia-oxidizing archaea

The genome of a putative Nitrosopumilaceae virus with a hypothetical spindle-shaped particle morphology was identified in the Yangshan Harbour metavirome from the East China Sea through protein similarity comparison and structure analysis. This discovery was accompanied by a set of 10 geographically dispersed close relatives found in the environmental virus datasets from typical locations of ammonia-oxidizing archaeon distribution. Its host prediction was supported by iPHoP prediction and protein sequence similarity. The structure of the predicted major capsid protein, together with the overall N-glycosylation site, the transmembrane helices prediction, the hydrophilicity profile, and the docking simulation of the major capsid proteins, indicate that these viruses resemble spindle-shaped viruses. It suggests a similarly assembled structure and, consequently, a possibly spindle-shaped morphology of these newly discovered archaeal viruses.

Host species and geography impact bee-associated RNA virus communities with evidence for isolation by distance in viral populations

Virus symbionts are important mediators of ecosystem function, yet we know little of their diversity and ecology in natural populations. The alarming decline of pollinating insects in many regions of the globe, especially the European honey bee, Apis mellifera, has been driven in part by worldwide transmission of virus pathogens. Previous work has examined the transmission of known honey bee virus pathogens to wild bee populations, but only a handful of studies have investigated the native viromes associated with wild bees, limiting epidemiological predictors associated with viral pathogenesis. Further, variation among different bee species might have important consequences in the acquisition and maintenance of bee-associated virome diversity. We utilized comparative metatranscriptomics to develop a baseline description of the RNA viromes associated with wild bee pollinators and to document viral diversity, community composition, and structure. Our sampling includes five wild-caught, native bee species that vary in social behavior as well as managed honey bees. We describe 26 putatively new RNA virus species based on RNA-dependent RNA polymerase phylogeny and show that each sampled bee species was associated with a specific virus community composition, even among sympatric populations of distinct host species. From 17 samples of a single host species, we recovered a single virus species despite over 600 km of distance between host populations and found strong evidence for isolation by distance in associated viral populations. Our work adds to the small number of studies examining viral prevalence and community composition in wild bees.

Exploring the Archaeal Virosphere by Metagenomics

During the past decade, environmental research has demonstrated that archaea are abundant and widespread in nature and play important ecological roles at a global scale. Currently, however, the majority of archaeal lineages cannot be cultivated under laboratory conditions and are known exclusively or nearly exclusively through metagenomics. A similar trend extends to the archaeal virosphere, where isolated representatives are available for a handful of model archaeal virus-host systems. Viral metagenomics provides an alternative way to circumvent the limitations of culture-based virus discovery and offers insight into the diversity, distribution, and environmental impact of uncultured archaeal viruses. Presently, metagenomics approaches have been successfully applied to explore the viromes associated with various lineages of extremophilic and mesophilic archaea, including Asgard archaea (Asgardarchaeota), ANME-1 archaea (Methanophagales), thaumarchaea (Nitrososphaeria), altiarchaea (Altiarchaeota), and marine group II archaea (Poseidoniales). Here, we provide an overview of methods widely used in archaeal virus metagenomics, covering metavirome preparation, genome annotation, phylogenetic and phylogenomic analyses, and archaeal host assignment. We hope that this summary will contribute to further exploration and characterization of the enigmatic archaeal virome lurking in diverse environments.

Interspecies Papillomavirus Type Infection and a Novel Papillomavirus Type in Red Ruffed Lemurs (Varecia rubra)

The Papillomaviridae are a family of vertebrate-infecting viruses of oncogenic potential generally thought to be host species- and tissue-specific. Despite their phylogenetic relatedness to humans, there is a scarcity of data on papillomaviruses (PVs) in speciose non-human primate lineages, particularly the lemuriform primates. Varecia variegata (black-and-white ruffed lemurs) and Varecia rubra (red ruffed lemurs), two closely related species comprising the Varecia genus, are critically endangered with large global captive populations. Varecia variegata papillomavirus (VavPV) types -1 and -2, the first PVs in lemurs with a fully identified genome, were previously characterized from captive V. variegata saliva. To build upon this discovery, saliva samples were collected from captive V. rubra with the following aims: (1) to identify PVs shared between V. variegata and V. rubra and (2) to characterize novel PVs in V. rubra to better understand PV diversity in the lemuriform primates. Three complete PV genomes were determined from V. rubra samples. Two of these PV genomes share 98% L1 nucleotide identity with VavPV2, denoting interspecies infection of V. rubra by VavPV2. This work represents the first reported case of interspecies PV infection amongst the strepsirrhine primates. The third PV genome shares <68% L1 nucleotide identity with that of all PVs. Thus, it represents a new PV species and has been named Varecia rubra papillomavirus 1 (VarPV1). VavPV1, VavPV2, and VarPV1 form a new clade within the Papillomaviridae family, likely representing a novel genus. Future work diversifying sample collection (i.e., lemur host species from multiple genera, sample type, geographic location, and wild populations) is likely to uncover a world of diverse lemur PVs.

Comparative genomics reveals the correlations of stress response genes and bacteriophages in developing antibiotic resistance of Staphylococcus saprophyticus

IMPORTANCE

Staphylococcus saprophyticus is the second most common bacteria associated with urinary tract infections (UTIs) in women. The antimicrobial treatment regimen for uncomplicated UTI is normally nitrofurantoin, trimethoprim-sulfamethoxazole (TMP-SMX), or a fluoroquinolone without routine susceptibility testing of S. saprophyticus recovered from urine specimens. However, TMP-SMX-resistant S. saprophyticus has been detected recently in UTI patients, as well as in our cohort. Herein, we investigated the understudied resistance patterns of this pathogenic species by linking genomic antibiotic resistance gene (ARG) content to susceptibility phenotypes. We describe ARG associations with known and novel SCCmec configurations as well as phage elements in S. saprophyticus, which may serve as intervention or diagnostic targets to limit resistance transmission. Our analyses yielded a comprehensive database of phenotypic data associated with the ARG sequence in clinical S. saprophyticus isolates, which will be crucial for resistance surveillance and prediction to enable precise diagnosis and effective treatment of S. saprophyticus UTIs.

Inoviruses

Tropini and Burckhardt introduce inoviruses and discuss how they are unique amongst bacteriophages.

Genomic surveillance of Clostridioides difficile transmission and virulence in a healthcare setting

Clostridioides difficile infection (CDI) is a major cause of healthcare-associated diarrhea, despite the widespread implementation of contact precautions for patients with CDI. Here, we investigate strain contamination in a hospital setting and genomic determinants of disease outcomes. Across two wards over six months, we selectively cultured C. difficile from patients (n=384) and their environments. Whole-genome sequencing (WGS) of 146 isolates revealed that most C. difficile isolates were from clade 1 (131/146, 89.7%), while only one isolate of the hypervirulent ST1 was recovered. Of culture-positive admissions (n=79), 19 (24%) of patients were colonized with toxigenic C. difficile on admission to the hospital. We defined 25 strain networks at ≤ 2 core gene SNPs; 2 of these networks contain strains from different patients. Strain networks were temporally linked (p<0.0001). To understand genomic correlates of disease, we conducted WGS on an additional cohort of C. difficile (n=102 isolates) from the same hospital and confirmed that clade 1 isolates are responsible for most CDI cases. We found that while toxigenic C. difficile isolates are associated with the presence of cdtR , nontoxigenic isolates have an increased abundance of prophages. Our pangenomic analysis of clade 1 isolates suggests that while toxin genes ( tcdABER and cdtR ) were associated with CDI symptoms, they are dispensable for patient colonization. These data indicate toxigenic and nontoxigenic C. difficile contamination persists in a hospital setting and highlight further investigation into how accessory genomic repertoires contribute to C. difficile colonization and disease.

Brief Report: The Virome of Bladder Tumors Arising in People Living With HIV

BACKGROUND

People living with HIV (PLWH) have elevated risk for developing virus-related cancers. Bladder cancer risk is not increased in PLWH but is elevated among immunosuppressed solid organ transplant recipients (SOTRs). BK polyomavirus and, to a lesser extent, other viruses have been detected in bladder cancers from SOTRs.

OBJECTIVE

To characterize the virome of bladder tumors in PLWH.

DESIGN

Retrospective case series.

METHODS

We sequenced DNA and RNA from archived formalin-fixed bladder tumors from PLWH. Nonhuman reads were assembled and matched to a database of known viruses.

RESULTS

Fifteen bladder tumors from PLWH (13 carcinomas, 2 benign tumors) were evaluated. Fourteen tumors were in men, and the median age at diagnosis was 59 years (median CD4 count 460 cells/mm3). All but 1 tumor yielded both sufficient DNA and RNA. One bladder cancer, arising in a 52-year-old man with a CD4 count of 271 cells/mm3, manifested diverse Alphatorquevirus DNA and RNA sequences. A second cancer arising in a 58-year-old male former smoker (CD4 count of 227 cells/mm3) also showed Alphatorquevirus and Gammatorquevirus DNA sequences. Neither tumor exhibited viral integration.

CONCLUSIONS

Alphatorqueviruses and Gammatorqueviruses are anelloviruses, which have also been detected in bladder cancers from SOTRs, but anelloviruses are common infections, and detection may simply reflect increased abundance in the setting of immunosuppression. The lack of detection of BK polyomavirus among bladder tumors from PLWH parallels the lower level of bladder cancer risk seen in PLWH compared with SOTRs, indirectly supporting a role for BK polyomavirus in causing the excess risk in SOTRs.

Phage-antibiotic combinations against multidrug-resistant Pseudomonas aeruginosa in in vitro static and dynamic biofilm models

Biofilm-producing Pseudomonas aeruginosa infections pose a severe threat to public health and are responsible for high morbidity and mortality. Phage-antibiotic combinations (PACs) are a promising strategy for combatting multidrug-resistant (MDR), extensively drug-resistant (XDR), and difficult-to-treat P. aeruginosa infections. Ten MDR/XDR P. aeruginosa strains and five P. aeruginosa-specific phages were genetically characterized and evaluated based upon their antibiotic susceptibilities and phage sensitivities. Two selected strains, AR351 (XDR) and I0003-1 (MDR), were treated singly and in combination with either a broad-spectrum or narrow-spectrum phage, phage EM-T3762627-2_AH (EM), or 14207, respectively, and bactericidal antibiotics of five classes in biofilm time-kill analyses. Synergy and/or bactericidal activity was demonstrated with all PACs against one or both drug-resistant P. aeruginosa strains (average reduction: -Δ3.32 log10 CFU/cm2). Slightly improved ciprofloxacin susceptibility was observed in both strains after exposure to phages (EM and 14207) in combination with ciprofloxacin and colistin. Based on phage cocktail optimization with four phages (EM, 14207, E20050-C (EC), and 109), we identified several effective phage-antibiotic cocktails for further analysis in a 4-day pharmacokinetic/pharmacodynamic in vitro biofilm model. Three-phage cocktail, EM + EC + 109, in combination with ciprofloxacin demonstrated the greatest biofilm reduction against AR351 (-Δ4.70 log10 CFU/cm2 from baseline). Of remarkable interest, the addition of phage 109 prevented phage resistance development to EM and EC in the biofilm model. PACs can demonstrate synergy and offer enhanced eradication of biofilm against drug-resistant P. aeruginosa while preventing the emergence of resistance.

Ciprofloxacin in combination with bacteriophage cocktails against multi-drug resistant Pseudomonas aeruginosa in ex vivo simulated endocardial vegetation models

Pseudomonas aeruginosa-associated infective endocarditis represents difficult-to-treat, deep-seated infections. Phage-antibiotic combinations have shown to eradicate multi-drug resistant (MDR) P. aeruginosa, limit the development of phage resistance, and restore antibiotic sensitivity. The objective of this study was to evaluate the activity of phage-ciprofloxacin (CIP) combinations in 4-day ex vivo simulated endocardial vegetation (SEV) models against drug-resistant P. aeruginosa isolates. Two P. aeruginosa isolates, extensively drug-resistant AR351 and MDR I0003-1, were selected for their drug resistance and sensitivity to phage. Three phages [LL-5504721-AH (LL), E2005-C (EC), and 109] and CIP were evaluated alone and in combination for their activity and influence on drug and phage resistance using 24-h time-kill analysis. The three-phage cocktail (q24h) in combination with CIP (400 mg q12h) was then tested in dynamic 4-day ex vivo SEV models, with reduction of log10 CFU/mL compared using ANOVA with Bonferroni analysis. Compared to other combinations, CIP-LL-EC-109 demonstrated synergistic and bactericidal activity from starting CFU/g against AR351 and I0003-1 (-Δ5.65 and 6.60 log10 CFU/g, respectively; P < 0.001). Additionally, CIP-LL-EC-109 mitigated phage resistance, while all other therapies had a high degree of resistance to >1 phages, and all phage-containing regimens prevented CIP mean inhibitory concentration increases compared to CIP alone for both AR351 and I0003-1 at 96 h.

Genomic Analysis and Taxonomic Characterization of Seven Bacteriophage Genomes Metagenomic-Assembled from the Dishui Lake

Viruses in aquatic ecosystems exhibit remarkable abundance and diversity. However, scattered studies have been conducted to mine uncultured viruses and identify them taxonomically in lake water. Here, whole genomes (29-173 kbp) of seven uncultured dsDNA bacteriophages were discovered in Dishui Lake, the largest artificial lake in Shanghai. We analyzed their genomic signatures and found a series of viral auxiliary metabolic genes closely associated with protein synthesis and host metabolism. Dishui Lake phages shared more genes with uncultivated environmental viruses than with reference viruses based on the gene-sharing network classification. Phylogeny of proteomes and comparative genomics delineated three new genera within two known viral families of Kyanoviridae and Autographiviridae, and four new families in Caudoviricetes for these seven novel phages. Their potential hosts appeared to be from the dominant bacterial phyla in Dishui Lake. Altogether, our study provides initial insights into the composition and diversity of bacteriophage communities in Dishui Lake, contributing valuable knowledge to the ongoing research on the roles played by viruses in freshwater ecosystems.

Phage-bacteria dynamics during the first years of life revealed by trans-kingdom marker gene analysis

Humans are colonized with commensal bacteria soon after birth, and, while this colonization is affected by lifestyle and other factors, bacterial colonization proceeds through well-studied phases. However, less is known about phage communities in early human development due to small study sizes, inability to leverage large databases, and lack of appropriate bioinformatics tools. In this study, whole genome shotgun sequencing data from the TEDDY study, composed of 12,262 longitudinal samples from 887 children in 4 countries, is reanalyzed to assess phage and bacterial dynamics simultaneously. Reads from these samples were mapped to marker genes from both bacteria and a new database of tens of thousands of phage taxa from human microbiomes. We uncover that each child is colonized by hundreds of different phages during the early years, and phages are more transitory than bacteria. Participants' samples continually harbor new phage species over time whereas the diversification of bacterial species begins to saturate. Phage data improves the ability for machine learning models to discriminate samples by country. Finally, while phage populations were individual-specific, striking patterns arose from the larger dataset, showing clear trends of ecological succession amongst phages, which correlated well with putative host bacteria. Improved understanding of phage-bacterial relationships may reveal new means by which to shape and modulate the microbiome and its constituents to improve health and reduce disease, particularly in vulnerable populations where antibiotic use and/or other more drastic measures may not be advised.

Identification of a novel polyomavirus in wild Sonoran Desert rodents of the family Heteromyidae

Rodents are the largest and most diverse group of mammals. Covering a wide range of structural and functional adaptations, rodents successfully occupy virtually every terrestrial habitat, and they are often found in close association with humans, domestic animals, and wildlife. Although a significant amount of research has focused on rodents' prominence as known reservoirs of zoonotic viruses, there has been less emphasis on the viral ecology of rodents in general. Here, we utilized a viral metagenomics approach to investigate polyomaviruses in wild rodents from the Baja California peninsula, Mexico, using fecal samples. We identified a novel polyomavirus in fecal samples from two rodent species, a spiny pocket mouse (Chaetodipus spinatus) and a Dulzura kangaroo rat (Dipodomys simulans). These two polyomaviruses represent a new species in the genus Betapolyomavirus. Sequences of this polyomavirus cluster phylogenetically with those of other rodent polyomaviruses and two other non-rodent polyomaviruses (WU and KI) that have been identified in the human respiratory tract. Through our continued work on seven species of rodents, we endeavor to explore the viral diversity associated with wild rodents on the Baja California peninsula and expand on current knowledge of rodent viral ecology and evolution.

Characterization of Diverse Anelloviruses, Cressdnaviruses, and Bacteriophages in the Human Oral DNA Virome from North Carolina (USA)

The diversity of viruses identified from the various niches of the human oral cavity-from saliva to dental plaques to the surface of the tongue-has accelerated in the age of metagenomics. This rapid expansion demonstrates that our understanding of oral viral diversity is incomplete, with only a few studies utilizing passive drool collection in conjunction with metagenomic sequencing methods. For this pilot study, we obtained 14 samples from healthy staff members working at the Duke Lemur Center (Durham, NC, USA) to determine the viral diversity that can be identified in passive drool samples from humans. The complete genomes of 3 anelloviruses, 9 cressdnaviruses, 4 Caudoviricetes large bacteriophages, 29 microviruses, and 19 inoviruses were identified in this study using high-throughput sequencing and viral metagenomic workflows. The results presented here expand our understanding of the vertebrate-infecting and microbe-infecting viral diversity of the human oral virome in North Carolina (USA).

Identification and characterization of a polyomavirus in the thornback skate (Raja clavata)

Members of the family Polyomaviridae have a circular double-stranded DNA genome that have been identified in various hosts ranging from mammals to arachnids. Here we report the identification and analysis of a complete genome sequence of a novel polyomavirus, Raja clavata polyomavirus (RcPyV1), from a cartilaginous fish, the thornback skate (Raja clavata). The genome sequence was determined using a metagenomics approach with an aim to provide baseline viral data in cartilaginous fish in different ecosystems. The RcPyV1 genome (4,195 nucleotides) had typical organization of polyomavirus, including early antigens (small T; Large T) encoded on one strand and late viral proteins (VP1; VP2) on the complementary strand. Maximum-likelihood phylogenetic analysis of the large T-antigen revealed that RcPyV1 clusters with a polyomavirus obtained from another cartilaginous fish, the guitarfish polyomavirus 1 (GfPyV1). These two share ~ 56% pairwise identity in LT and VP1 protein sequences. These analyses support the hypothesis that cartilaginous fishes have a specific lineage of polyomaviruses.

Viruses of Freshwater Mussels during Mass Mortality Events in Oregon and Washington, USA

Freshwater mussels (Unionida) are globally imperiled, in part due to largely unexplained mass mortality events (MMEs). While recent studies have begun to investigate the possibility that mussel MMEs in the Eastern USA may be caused by infectious diseases, mussels in the Western USA have received relatively little attention in this regard. We conducted a two-year epidemiologic investigation of the role of viruses in ongoing MMEs of the Western pearlshell (Margaritifera falcata) and the Western ridged mussel (Gonidea angulata) in the Chehalis River and Columbia River watersheds in the Western USA. We characterized viromes of mussel hemolymph from 5 locations in 2018 and 2020 using metagenomic methods and identified 557 viruses based on assembled contiguous sequences, most of which are novel. We also characterized the distribution and diversity of a previously identified mussel Gammarhabdovirus related to pathogenic finfish viruses. Overall, we found few consistent associations between viruses and mussel health status. Variation in mussel viromes was most strongly driven by location, with little influence from date, species, or health status, though these variables together only explained ~1/3 of variation in virome composition. Our results demonstrate that Western freshwater mussels host remarkably diverse viromes, but no single virus or combination of viruses appears to be associated with morbidity or mortality during MMEs. Our findings have implications for the conservation of imperiled freshwater mussels, including efforts to enhance natural populations through captive propagation.

Phages are unrecognized players in the ecology of the oral pathogen Porphyromonas gingivalis

BACKGROUND

Porphyromonas gingivalis (hereafter "Pg") is an oral pathogen that has been hypothesized to act as a keystone driver of inflammation and periodontal disease. Although Pg is most readily recovered from individuals with actively progressing periodontal disease, healthy individuals and those with stable non-progressing disease are also colonized by Pg. Insights into the factors shaping the striking strain-level variation in Pg, and its variable associations with disease, are needed to achieve a more mechanistic understanding of periodontal disease and its progression. One of the key forces often shaping strain-level diversity in microbial communities is infection of bacteria by their viral (phage) predators and symbionts. Surprisingly, although Pg has been the subject of study for over 40 years, essentially nothing is known of its phages, and the prevailing paradigm is that phages are not important in the ecology of Pg.

RESULTS

Here we systematically addressed the question of whether Pg are infected by phages-and we found that they are. We found that prophages are common in Pg, they are genomically diverse, and they encode genes that have the potential to alter Pg physiology and interactions. We found that phages represent unrecognized targets of the prevalent CRISPR-Cas defense systems in Pg, and that Pg strains encode numerous additional mechanistically diverse candidate anti-phage defense systems. We also found that phages and candidate anti-phage defense system elements together are major contributors to strain-level diversity and the species pangenome of this oral pathogen. Finally, we demonstrate that prophages harbored by a model Pg strain are active in culture, producing extracellular viral particles in broth cultures.

CONCLUSION

This work definitively establishes that phages are a major unrecognized force shaping the ecology and intra-species strain-level diversity of the well-studied oral pathogen Pg. The foundational phage sequence datasets and model systems that we establish here add to the rich context of all that is already known about Pg, and point to numerous avenues of future inquiry that promise to shed new light on fundamental features of phage impacts on human health and disease broadly. Video Abstract.

Diversity, Relationship, and Distribution of Virophages and Large Algal Viruses in Global Ocean Viromes

Virophages are a group of small double-stranded DNA viruses that replicate and proliferate with the help of the viral factory of large host viruses. They are widely distributed in aquatic environments but are more abundant in freshwater ecosystems. Here, we mined the Global Ocean Viromes 2.0 (GOV 2.0) dataset for the diversity, distribution, and association of virophages and their potential host large viruses in marine environments. We identified 94 virophage sequences (>5 kbp in length), of which eight were complete genomes. The MCP phylogenetic tree showed that the GOV virophages were widely distributed on the global virophage tree but relatively clustered on three major branches. The gene-sharing network divided GOV virophages into 21 outliers, 2 overlaps, and 14 viral clusters, of which 4 consisted of only the GOV virophages. We also identified 45 large virus sequences, 8 of which were >100 kbp in length and possibly involved in cell-virus-virophage (C-V-v) trisome relationships. The potential eukaryotic hosts of these eight large viruses and the eight virophages with their complete genomes identified are likely to be algae, based on comparative genomic analysis. Both homologous gene and codon usage analyses support a possible interaction between a virophage (GOVv18) and a large algal virus (GOVLV1). These results indicate that diverse and novel virophages and large viruses are widespread in global marine environments, suggesting their important roles and the presence of complicated unknown C-V-v relationships in marine ecosystems.

Roving methyltransferases generate a mosaic epigenetic landscape and influence evolution in Bacteroides fragilis group

Three types of DNA methyl modifications have been detected in bacterial genomes, and mechanistic studies have demonstrated roles for DNA methylation in physiological functions ranging from phage defense to transcriptional control of virulence and host-pathogen interactions. Despite the ubiquity of methyltransferases and the immense variety of possible methylation patterns, epigenomic diversity remains unexplored for most bacterial species. Members of the Bacteroides fragilis group (BFG) reside in the human gastrointestinal tract as key players in symbiotic communities but also can establish anaerobic infections that are increasingly multi-drug resistant. In this work, we utilize long-read sequencing technologies to perform pangenomic (n = 383) and panepigenomic (n = 268) analysis of clinical BFG isolates cultured from infections seen at the NIH Clinical Center over four decades. Our analysis reveals that single BFG species harbor hundreds of DNA methylation motifs, with most individual motif combinations occurring uniquely in single isolates, implying immense unsampled methylation diversity within BFG epigenomes. Mining of BFG genomes identified more than 6000 methyltransferase genes, approximately 1000 of which were associated with intact prophages. Network analysis revealed substantial gene flow among disparate phage genomes, implying a role for genetic exchange between BFG phages as one of the ultimate sources driving BFG epigenome diversity.

Using archived and biocollection samples towards deciphering the DNA virus diversity associated with rodent species in the families cricetidae and heteromyidae

Rodentia is the most speciose order of mammals, and they are known to harbor a wide range of viruses. Although there has been significant research on zoonotic viruses in rodents, research on the diversity of other viruses has been limited, especially for rodents in the families Cricetidae and Heteromyidae. In fecal and liver samples of nine species of rodents, we identify 346 distinct circular DNA viral genomes. Of these, a large portion are circular, single-stranded DNA viruses in the families Anelloviridae (n = 3), Circoviridae (n = 5), Genomoviridae (n = 7), Microviridae (n = 297), Naryaviridae (n = 4), Vilyaviridae (n = 15) and in the phylum Cressdnaviricota (n = 13) that cannot be assigned established families. We also identified two large bacteriophages of 36 and 50 kb that are part of the class Caudoviricetes. Some of these viruses are clearly those that infect rodents, however, most of these likely infect various organisms associated with rodents, their environment or their diet.

Large language models improve annotation of viral proteins

Viral sequences are poorly annotated in environmental samples, a major roadblock to understanding how viruses influence microbial community structure. Current annotation approaches rely on alignment-based sequence ho-mology methods, which are limited by available viral sequences and sequence divergence in viral proteins. Here, we show that protein language model representations capture viral protein function beyond the limits of remote sequence homology by targeting two axes of viral sequence annotation: systematic labeling of protein families and function identification for biologic discovery. Protein language model representations capture protein functional properties specific to viruses and expand the annotated fraction of ocean virome viral protein sequences by 37%. Among unannotated viral protein families, we identify a novel DNA editing protein family that defines a new mobile element in marine picocyanobacteria. Protein language models thus significantly enhance remote homology detection of viral proteins and can be utilized to enable new biological discovery across diverse functional categories.

Evidence for virus-mediated oncogenesis in bladder cancers arising in solid organ transplant recipients

A small percentage of bladder cancers in the general population have been found to harbor DNA viruses. In contrast, up to 25% of tumors of solid organ transplant recipients, who are at an increased risk of developing bladder cancer and have an overall poorer outcomes, harbor BK polyomavirus (BKPyV). To better understand the biology of the tumors and the mechanisms of carcinogenesis from potential oncoviruses, we performed whole genome and transcriptome sequencing on bladder cancer specimens from 43 transplant patients. Nearly half of the tumors from this patient population contained viral sequences. The most common were from BKPyV (N=9, 21%), JC polyomavirus (N=7, 16%), carcinogenic human papillomaviruses (N=3, 7%), and torque teno viruses (N=5, 12%). Immunohistochemistry revealed variable Large T antigen expression in BKPyV-positive tumors ranging from 100% positive staining of tumor tissue to less than 1%. In most cases of BKPyV-positive tumors, the viral genome appeared to be clonally integrated into the host chromosome consistent with microhomology-mediated end joining and coincided with focal amplifications of the tumor genome similar to other virus-mediated cancers. Significant changes in host gene expression consistent with the functions of BKPyV Large T antigen were also observed in these tumors. Lastly, we identified four mutation signatures in our cases, with those attributable to APOBEC3 and SBS5 being the most abundant. Mutation signatures associated with an antiviral drug, ganciclovir, and aristolochic acid, a nephrotoxic compound found in some herbal medicines, were also observed. The results suggest multiple pathways to carcinogenesis in solid organ transplant recipients with a large fraction being virus-associated.

Cervicovaginal Human Papillomavirus Genomes, Microbiota Composition and Cytokine Concentrations in South African Adolescents

The interaction between cervicovaginal virome, bacteriome and genital inflammation has not been extensively investigated. We assessed the vaginal DNA virome from 33 South African adolescents (15-19 years old) using shotgun DNA sequencing of purified virions. We present analyses of eukaryote-infecting DNA viruses, with a focus on human papillomavirus (HPV) genomes and relate these to the vaginal bacterial microbiota (assessed by 16S rRNA gene sequencing) and cytokines (assessed by Luminex). The DNA virome included single-stranded (Anelloviridae, Genomoviridae) and double-stranded DNA viruses (Adenoviridae, Alloherpesviridae, Herpesviridae, Marseilleviridae, Mimiviridae, Polyomaviridae, Poxviridae). We identified 110 unique, complete HPV genomes within two genera (Alphapapillomavirus and Gammapapillomavirus) representing 40 HPV types and 12 species. Of the 40 HPV types identified, 35 showed positive co-infection patterns with at least one other type, mainly HPV-16. HPV-35, a high-risk genotype currently not targeted by available vaccines, was the most prevalent HPV type identified in this cohort. Bacterial taxa commonly associated with bacterial vaginosis also correlated with the presence of HPV. Bacterial vaginosis, rather than HPV, was associated with increased genital inflammation. This study lays the foundation for future work characterizing the vaginal virome and its role in women's health.

Diverse DNA virus genomes identified in fecal samples of Mexican free-tailed bats (Tadarida brasiliensis) captured in Chiricahua Mountains of southeast Arizona (USA)

Bats (order Chiroptera) are some of the most abundant mammals on earth and their species ecology strongly influences zoonotic potential. While substantial research has been conducted on bat-associated viruses, particularly on those that can cause disease in humans and/or livestock, globally, limited research has focused on endemic bats in the USA. The southwest region of the US is of particular interest because of its high diversity of bat species. We identified 39 single-stranded DNA virus genomes in the feces of Mexican free-tailed bats (Tadarida brasiliensis) sampled in the Rucker Canyon (Chiricahua Mountains) of southeast Arizona (USA). Twenty-eight of these belong to the virus families Circoviridae (n = 6), Genomoviridae (n = 17), and Microviridae (n = 5). Eleven viruses cluster with other unclassified cressdnaviruses. Most of the viruses identified represent new species. Further research on identification of novel bat-associated cressdnaviruses and microviruses is needed to provide greater insights regarding their co-evolution and ecology relative to bats.

What Lies Beneath? Taking the Plunge into the Murky Waters of Phage Biology

The sequence revolution revealed that bacteria-infecting viruses, known as phages, are Earth's most abundant biological entities. Phages have far-reaching impacts on the form and function of microbial communities and play a fundamental role in ecological processes. However, even well into the sequencing revolution, we have only just begun to explore the murky waters around the phage biology iceberg. Many viral reads cannot be assigned to a culturable isolate, and reference databases are biased toward more easily collectible samples, which likely distorts our conclusions. This minireview points out alternatives to mapping reads to reference databases and highlights innovative bioinformatic and experimental approaches that can help us overcome some of the challenges in phage research and better decipher the impact of phages on microbial communities. Moving beyond the identification of novel phages, we highlight phage metabolomics as an important influencer of bacterial host cell physiology and hope to inspire the reader to consider the effects of phages on host metabolism and ecosystems at large. We encourage researchers to report unassigned/unknown sequencing reads and contigs and to continue developing alternative methods to investigate phages within sequence data.

Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and Casposons

Profile hidden Markov models (HMMs) are a powerful way of modeling biological sequence diversity and constitute a very sensitive approach to detecting divergent sequences. Here, we report the development of protocols for the rational design of profile HMMs. These methods were implemented on TABAJARA, a program that can be used to either detect all biological sequences of a group or discriminate specific groups of sequences. By calculating position-specific information scores along a multiple sequence alignment, TABAJARA automatically identifies the most informative sequence motifs and uses them to construct profile HMMs. As a proof-of-principle, we applied TABAJARA to generate profile HMMs for the detection and classification of two viral groups presenting different evolutionary rates: bacteriophages of the Microviridae family and viruses of the Flavivirus genus. We obtained conserved models for the generic detection of any Microviridae or Flavivirus sequence, and profile HMMs that can specifically discriminate Microviridae subfamilies or Flavivirus species. In another application, we constructed Cas1 endonuclease-derived profile HMMs that can discriminate CRISPRs and casposons, two evolutionarily related transposable elements. We believe that the protocols described here, and implemented on TABAJARA, constitute a generic toolbox for generating profile HMMs for the highly sensitive and specific detection of sequence classes.

When Plaquing Is Not Possible: Computational Methods for Detecting Induced Phages

High-throughput sequencing of microbial communities has uncovered a large, diverse population of phages. Frequently, phages found are integrated into their bacterial host genome. Distinguishing between phages in their integrated (lysogenic) and unintegrated (lytic) stage can provide insight into how phages shape bacterial communities. Here we present the Prophage Induction Estimator (PIE) to identify induced phages in genomic and metagenomic sequences. PIE takes raw sequencing reads and phage sequence predictions, performs read quality control, read assembly, and calculation of phage and non-phage sequence abundance and completeness. The distribution of abundances for non-phage sequences is used to predict induced phages with statistical confidence. In silico tests were conducted to benchmark this tool finding that PIE can detect induction events as well as phages with a relatively small burst size (10×). We then examined isolate genome sequencing data as well as a mock community and urinary metagenome data sets and found instances of induced phages in all three data sets. The flexibility of this software enables users to easily include phage predictions from their preferred tool of choice or phage sequences of interest. Thus, genomic and metagenomic sequencing now not only provides a means for discovering and identifying phage sequences but also the detection of induced prophages.

Evaluation of computational phage detection tools for metagenomic datasets

Introduction

As new computational tools for detecting phage in metagenomes are being rapidly developed, a critical need has emerged to develop systematic benchmarks.

Methods

In this study, we surveyed 19 metagenomic phage detection tools, 9 of which could be installed and run at scale. Those 9 tools were assessed on several benchmark challenges. Fragmented reference genomes are used to assess the effects of fragment length, low viral content, phage taxonomy, robustness to eukaryotic contamination, and computational resource usage. Simulated metagenomes are used to assess the effects of sequencing and assembly quality on the tool performances. Finally, real human gut metagenomes and viromes are used to assess the differences and similarities in the phage communities predicted by the tools.

Results

We find that the various tools yield strikingly different results. Generally, tools that use a homology approach (VirSorter, MARVEL, viralVerify, VIBRANT, and VirSorter2) demonstrate low false positive rates and robustness to eukaryotic contamination. Conversely, tools that use a sequence composition approach (VirFinder, DeepVirFinder, Seeker), and MetaPhinder, have higher sensitivity, including to phages with less representation in reference databases. These differences led to widely differing predicted phage communities in human gut metagenomes, with nearly 80% of contigs being marked as phage by at least one tool and a maximum overlap of 38.8% between any two tools. While the results were more consistent among the tools on viromes, the differences in results were still significant, with a maximum overlap of 60.65%. Discussion: Importantly, the benchmark datasets developed in this study are publicly available and reusable to enable the future comparability of new tools developed.

Novel polyomaviruses identified in fecal samples from four carnivore species

Polyomaviruses are oncogenic viruses that are generally thought to have co-evolved with their hosts. While primate and rodent polyomaviruses are increasingly well-studied, less is known about polyomaviruses that infect other mammals. In an effort to gain insight into polyomaviruses associated with carnivores, we surveyed fecal samples collected in the USA from bobcats (Lynx rufus), pumas (Puma concolor), Canada lynxes (Lynx canadensis), and grizzly bears (Ursus arctos). Using a viral metagenomic approach, we identified six novel polyomavirus genomes. Surprisingly, four of the six genomes showed a phylogenetic relationship to polyomaviruses found in prey animals. These included a putative rabbit polyomavirus from a bobcat fecal sample and two possible deer-trophic polyomaviruses from Canada lynx feces. One polyomavirus found in a grizzly bear sample was found to be phylogenetically distant from previously identified polyomaviruses. Further analysis of the grizzly bear fecal sample showed that it contained anelloviruses that are known to infect pigs, suggesting that the bear might have preyed on a wild or domestic pig. Interestingly, a polyomavirus genome identified in a puma fecal sample was found to be closely related both to raccoon polyomavirus 1 and to Lyon-IARC polyomavirus, the latter of which was originally identified in human saliva and skin swab specimens but has since been found in samples from domestic cats (Felis catus).