Improvements in viral gene annotation using large language models and soft alignments.

BACKGROUND: The annotation of protein sequences in public databases has long posed a challenge in molecular biology. This issue is particularly acute for viral proteins, which demonstrate limited homology to known proteins when using alignment, k-mer, or profile-based homology search approaches. A novel methodology employing Large Language Models (LLMs) addresses this methodological challenge by annotating protein sequences based on embeddings. RESULTS: Central to our contribution is the soft alignment algorithm, drawing from traditional protein alignment but leveraging embedding similarity at the amino acid level to bypass the need for conventional scoring matrices. This method not only surpasses pooled embedding-based models in efficiency but also in interpretability, enabling users to easily trace homologous amino acids and delve deeper into the alignments. Far from being a black box, our approach provides transparent, BLAST-like alignment visualizations, combining traditional biological research with AI advancements to elevate protein annotation through embedding-based analysis while ensuring interpretability. Tests using the Virus Orthologous Groups and ViralZone protein databases indicated that the novel soft alignment approach recognized and annotated sequences that both blastp and pooling-based methods, which are commonly used for sequence annotation, failed to detect. CONCLUSION: The embeddings approach shows the great potential of LLMs for enhancing protein sequence annotation, especially in viral genomics. These findings present a promising avenue for more efficient and accurate protein function inference in molecular biology.

Ubiquitous, B12-dependent virioplankton utilizing ribonucleotide-triphosphate reductase demonstrate interseasonal dynamics and associate with a diverse range of bacterial hosts in the pelagic ocean.

Through infection and lysis of their coexisting bacterial hosts, viruses impact the biogeochemical cycles sustaining globally significant pelagic oceanic ecosystems. Currently, little is known of the ecological interactions between lytic viruses and their bacterial hosts underlying these biogeochemical impacts at ecosystem scales. This study focused on populations of lytic viruses carrying the B12-dependent Class II monomeric ribonucleotide reductase (RNR) gene, ribonucleotide-triphosphate reductase (Class II RTPR), documenting seasonal changes in pelagic virioplankton and bacterioplankton using amplicon sequences of Class II RTPR and the 16S rRNA gene, respectively. Amplicon sequence libraries were analyzed using compositional data analysis tools that account for the compositional nature of these data. Both virio- and bacterioplankton communities responded to environmental changes typically seen across seasonal cycles as well as shorter term upwelling-downwelling events. Defining Class II RTPR-carrying viral populations according to major phylogenetic clades proved a more robust means of exploring virioplankton ecology than operational taxonomic units defined by percent sequence homology. Virioplankton Class II RTPR populations showed positive associations with a broad phylogenetic diversity of bacterioplankton including dominant taxa within pelagic oceanic ecosystems such as Prochlorococcus and SAR11. Temporal changes in Class II RTPR virioplankton, occurring as both free viruses and within infected cells, indicated possible viral-host pairs undergoing sustained infection and lysis cycles throughout the seasonal study. Phylogenetic relationships inferred from Class II RTPR sequences mirrored ecological patterns in virio- and bacterioplankton populations demonstrating possible genome to phenome associations for an essential viral replication gene.

Spontaneously Produced Lysogenic Phages Are an Important Component of the Soybean Bradyrhizobium Mobilome.

The ability of Bradyrhizobium spp. to nodulate and fix atmospheric nitrogen in soybean root nodules is critical to meeting humanity's nutritional needs. The intricacies of soybean bradyrhizobia-plant interactions have been studied extensively; however, bradyrhizobial ecology as influenced by phages has received somewhat less attention, even though these interactions may significantly impact soybean yield. In batch culture, four soybean bradyrhizobia strains, Bradyrhizobium japonicum S06B (S06B-Bj), B. japonicum S10J (S10J-Bj), Bradyrhizobium diazoefficiens USDA 122 (USDA 122-Bd), and Bradyrhizobium elkanii USDA 76T (USDA 76-Be), spontaneously (without apparent exogenous chemical or physical induction) produced tailed phages throughout the growth cycle; for three strains, phage concentrations exceeded cell numbers by ~3-fold after 48 h of incubation. Phage terminase large-subunit protein phylogeny revealed possible differences in phage packaging and replication mechanisms. Bioinformatic analyses predicted multiple prophage regions within each soybean bradyrhizobia genome, preventing accurate identification of spontaneously produced prophage (SPP) genomes. A DNA sequencing and mapping approach accurately delineated the boundaries of four SPP genomes within three of the soybean bradyrhizobia chromosomes and suggested that the SPPs were capable of transduction. In addition to the phages, S06B-Bj and USDA 76-Be contained three to four times more insertion sequences (IS) and large, conjugable, broad host range plasmids, both of which are known drivers of horizontal gene transfer (HGT) in soybean bradyrhizobia. These factors indicate that SPP along with IS and plasmids participate in HGT, drive bradyrhizobia evolution, and play an outsized role in bradyrhizobia ecology. IMPORTANCE Previous studies have shown that IS and plasmids mediate HGT of symbiotic nodulation (nod) genes in soybean bradyrhizobia; however, these events require close cell-to-cell contact, which could be limited in soil environments. Bacteriophage-assisted gene transduction through spontaneously produced prophages provides a stable means of HGT not limited by the constraints of proximal cell-to-cell contact. These phage-mediated HGT events may shape soybean bradyrhizobia population ecology, with concomitant impacts on soybean agriculture.

Novel Viral DNA Polymerases From Metagenomes Suggest Genomic Sources of Strand-Displacing Biochemical Phenotypes.

Viruses are the most abundant and diverse biological entities on the planet and constitute a significant proportion of Earth's genetic diversity. Most of this diversity is not represented by isolated viral-host systems and has only been observed through sequencing of viral metagenomes (viromes) from environmental samples. Viromes provide snapshots of viral genetic potential, and a wealth of information on viral community ecology. These data also provide opportunities for exploring the biochemistry of novel viral enzymes. The in vitro biochemical characteristics of novel viral DNA polymerases were explored, testing hypothesized differences in polymerase biochemistry according to protein sequence phylogeny. Forty-eight viral DNA Polymerase I (PolA) proteins from estuarine viromes, hot spring metagenomes, and reference viruses, encompassing a broad representation of currently known diversity, were synthesized, expressed, and purified. Novel functionality was shown in multiple PolAs. Intriguingly, some of the estuarine viral polymerases demonstrated moderate to strong innate DNA strand displacement activity at high enzyme concentration. Strand-displacing polymerases have important technological applications where isothermal reactions are desirable. Bioinformatic investigation of genes neighboring these strand displacing polymerases found associations with SNF2 helicase-associated proteins. The specific function of SNF2 family enzymes is unknown for prokaryotes and viruses. In eukaryotes, SNF2 enzymes have chromatin remodeling functions but do not separate nucleic acid strands. This suggests the strand separation function may be fulfilled by the DNA polymerase for viruses carrying SNF2 helicase-associated proteins. Biochemical data elucidated from this study expands understanding of the biology and ecological behavior of unknown viruses. Moreover, given the numerous biotechnological applications of viral DNA polymerases, novel viral polymerases discovered within viromes may be a rich source of biological material for further in vitro DNA amplification advancements.

Towards an integrative view of virus phenotypes.

Understanding how phenotypes emerge from genotypes is a foundational goal in biology. As challenging as this task is when considering cellular life, it is further complicated in the case of viruses. During replication, a virus as a discrete entity (the virion) disappears and manifests itself as a metabolic amalgam between the virus and the host (the virocell). Identifying traits that unambiguously constitute a virus's phenotype is straightforward for the virion, less so for the virocell. Here, we present a framework for categorizing virus phenotypes that encompasses both virion and virocell stages and considers functional and performance traits of viruses in the context of fitness. Such an integrated view of virus phenotype is necessary for comprehensive interpretation of viral genome sequences and will advance our understanding of viral evolution and ecology.

CRISPR Spacers Indicate Preferential Matching of Specific Virioplankton Genes.

Viral infection exerts selection pressure on marine microbes, as virus-induced cell lysis causes 20 to 50% of cell mortality, resulting in fluxes of biomass into oceanic dissolved organic matter. Archaeal and bacterial populations can defend against viral infection using the clustered regularly interspaced short palindromic repeat (CRISPR)-associated (Cas) system, which relies on specific matching between a spacer sequence and a viral gene. If a CRISPR spacer match to any gene within a viral genome is equally effective in preventing lysis, no viral genes should be preferentially matched by CRISPR spacers. However, if there are differences in effectiveness, certain viral genes may demonstrate a greater frequency of CRISPR spacer matches. Indeed, homology search analyses of bacterioplankton CRISPR spacer sequences against virioplankton sequences revealed preferential matching of replication proteins, nucleic acid binding proteins, and viral structural proteins. Positive selection pressure for effective viral defense is one parsimonious explanation for these observations. CRISPR spacers from virioplankton metagenomes preferentially matched methyltransferase and phage integrase genes within virioplankton sequences. These virioplankton CRISPR spacers may assist infected host cells in defending against competing phage. Analyses also revealed that half of the spacer-matched viral genes were unknown, some genes matched several spacers, and some spacers matched multiple genes, a many-to-many relationship. Thus, CRISPR spacer matching may be an evolutionary algorithm, agnostically identifying those genes under stringent selection pressure for sustaining viral infection and lysis. Investigating this subset of viral genes could reveal those genetic mechanisms essential to virus-host interactions and provide new technologies for optimizing CRISPR defense in beneficial microbes.IMPORTANCE The CRISPR-Cas system is one means by which bacterial and archaeal populations defend against viral infection which causes 20 to 50% of cell mortality in the ocean. We tested the hypothesis that certain viral genes are preferentially targeted for the initial attack of the CRISPR-Cas system on a viral genome. Using CASC, a pipeline for CRISPR spacer discovery, and metagenome data from oceanic microbes and viruses, we found a clear subset of viral genes with high match frequencies to CRISPR spacers. Moreover, we observed a many-to-many relationship of spacers and viral genes. These high-match viral genes were involved in nucleotide metabolism, DNA methylation, and viral structure. It is possible that CRISPR spacer matching is an evolutionary algorithm pointing to those viral genes most important to sustaining infection and lysis. Studying these genes may advance the understanding of virus-host interactions in nature and provide new technologies for leveraging CRISPR-Cas systems in beneficial microbes.

Family A DNA Polymerase Phylogeny Uncovers Diversity and Replication Gene Organization in the Virioplankton.

Shotgun metagenomics, which allows for broad sampling of viral diversity, has uncovered genes that are widely distributed among virioplankton populations and show linkages to important biological features of unknown viruses. Over 25% of known dsDNA phage carry the DNA polymerase I (polA) gene, making it one of the most widely distributed phage genes. Because of its pivotal role in DNA replication, this enzyme is linked to phage lifecycle characteristics. Previous research has suggested that a single amino acid substitution might be predictive of viral lifestyle. In this study Chesapeake Bay virioplankton were sampled by shotgun metagenomic sequencing (using long and short read technologies). More polA sequences were predicted from this single viral metagenome (virome) than from 86 globally distributed virome libraries (ca. 2,100, and 1,200, respectively). The PolA peptides predicted from the Chesapeake Bay virome clustered with 69% of PolA peptides from global viromes; thus, remarkably the Chesapeake Bay virome captured the majority of known PolA peptide diversity in viruses. This deeply sequenced virome also expanded the diversity of PolA sequences, increasing the number of PolA clusters by 44%. Contigs containing polA sequences were also used to examine relationships between phylogenetic clades of PolA and other genes within unknown viral populations. Phylogenic analysis revealed five distinct groups of phages distinguished by the amino acids at their 762 (Escherichia coli IAI39 numbering) positions and replication genes. DNA polymerase I sequences from Tyr762 and Phe762 groups were most often neighbored by ring-shaped superfamily IV helicases and ribonucleotide reductases (RNRs). The Leu762 groups had non-ring shaped helicases from superfamily II and were further distinguished by an additional helicase gene from superfamily I and the lack of any identifiable RNR genes. Moreover, we found that the inclusion of ribonucleotide reductase associated with PolA helped to further differentiate phage diversity, chiefly within lytic podovirus populations. Altogether, these data show that DNA Polymerase I is a useful marker for observing the diversity and composition of the virioplankton and may be a driving factor in the divergence of phage replication components.

Aquatic hazard, bioaccumulation and screening risk assessment for ammonium 2,3,3,3-tetrafluoro-2-(heptafluoropropoxy)-propanoate.

The fluoropolymer manufacturing industry is moving to alternative polymerization processing aid technologies with more favorable toxicological and environmental profiles as part of a commitment to curtail the use of long-chain perfluoroalkyl acids (PFAAs). To facilitate the environmental product stewardship assessment and premanufacture notification (PMN) process for a candidate replacement chemical, we conducted acute and chronic aquatic toxicity tests to evaluate the toxicity of ammonium 2,3,3,3-tetrafluoro-2-(heptafluoropropoxy)-propanoate (C6HF11O3.H3N) or the acid form of the substance to the cladoceran, Daphnia magna, the green alga, Pseudokirchneriella subcapitata, and a number of freshwater fish species including the rainbow trout, Oncorhynchus mykiss, In addition, testing with the common carp, Cyprinus carpio, was conducted to determine the bioconcentration potential of the acid form of the compound. Based on the relevant criteria in current regulatory frameworks, the results of the aquatic toxicity and bioconcentration studies indicate the substance is of low concern for aquatic hazard and bioconcentration in aquatic organisms. Evaluation of environmental monitoring data in conjunction with the predicted no effect concentration (PNEC) based on the available data suggest low risk to aquatic organisms.

Aquatic hazard, bioaccumulation and screening risk assessment for 6:2 fluorotelomer sulfonate.

This study assessed the aquatic toxicity and bioaccumulation potential of 6:2 fluorotelomer sulfonate (6:2 FTSA). Acute and chronic aquatic hazard endpoints indicate 6:2 FTSA is not classified for aquatic hazard according to GHS or European CLP legislation. The aqueous bioconcentration factors for 6:2 FTSA were <40 and the dietary assimilation efficiency, growth corrected half-life and dietary biomagnification factor (BMF) were 0.435, 23.1d and 0.295, respectively. These data indicate that 6:2 FTSA is not bioaccumulative in aquatic organisms. Comparison of PNECs with the reported surface water concentrations (non-spill situations) suggests low risk to aquatic organisms from 6:2 FTSA. Future studies are needed to elucidate the biotic and abiotic fate of commercial AFFF surfactants in the environment.

Intra- and interlaboratory reliability of a cryopreserved trout hepatocyte assay for the prediction of chemical bioaccumulation potential.

Measured rates of intrinsic clearance determined using cryopreserved trout hepatocytes can be extrapolated to the whole animal as a means of improving modeled bioaccumulation predictions for fish. To date, however, the intra- and interlaboratory reliability of this procedure has not been determined. In the present study, three laboratories determined in vitro intrinsic clearance of six reference compounds (benzo[a]pyrene, 4-nonylphenol, di-tert-butyl phenol, fenthion, methoxychlor and o-terphenyl) by conducting substrate depletion experiments with cryopreserved trout hepatocytes from a single source. O-terphenyl was excluded from the final analysis due to nonfirst-order depletion kinetics and significant loss from denatured controls. For the other five compounds, intralaboratory variability (% CV) in measured in vitro intrinsic clearance values ranged from 4.1 to 30%, while interlaboratory variability ranged from 27 to 61%. Predicted bioconcentration factors based on in vitro clearance values exhibited a reduced level of interlaboratory variability (5.3-38% CV). The results of this study demonstrate that cryopreserved trout hepatocytes can be used to reliably obtain in vitro intrinsic clearance of xenobiotics, which provides support for the application of this in vitro method in a weight-of-evidence approach to chemical bioaccumulation assessment.

Comparative acute freshwater hazard assessment and preliminary PNEC development for eight fluorinated acids.

Short-term 48, 72 and 96-h aquatic toxicity tests were conducted to evaluate the acute toxicity of eight fluorinated acids to the cladoceran, Daphnia magna, the green alga, Pseudokirchneriella subcapitata, and the rainbow trout, Oncorhynchus mykiss or the fathead minnow, Pimephales promelas. The eight fluorinated acids studied were tridecafluorohexyl ethanoic acid (6:2 FTCA), heptadecafluorooctyl ethanoic acid (8:2 FTCA), 2H-dodecafluoro-2-octenoic acid (6:2 FTUCA), 2H-hexadecafluoro-2-decenoic acid (8:2 FTUCA), 2H,2H,3H,3H-undecafluoro octanoic acid (5:3 acid), 2H,2H,3H,3H-pentadecafluoro decanoic acid (7:3 acid), n-perfluoropentanoic acid (PFPeA) and n-perfluorodecanoic acid (PFDA). The results of the acute toxicity tests conducted during this study suggest that the polyfluorinated acids, 8:2 FTCA, 8:2 FTUCA, 6:2 FTCA, 6:2 FTUCA, 7:3 acid and 5:3 acid, and the perfluorinated acids PFPeA and PFDA, are generally of low to medium concern based on evaluation of their acute freshwater toxicity (EC/LC50s typically between 1 and >100 mg L(-1)) using the USEPA TSCA aquatic toxicity evaluation paradigm. For the polyfluorinated acids, aquatic toxicity generally decreased as the number of fluorinated carbons decreased and as the overall carbon chain length decreased from 12 to 8. Acute aquatic toxicity of the 5 and 10 carbon perfluorocarboxylic acids (EC/LC50s between 10.6 and >100 mg L(-1)) was greater or similar to that of the 6-9 carbon perfluorocarboxylic acids (EC/LC50s>96.5 mg L(-1)). This study also provides the first report of the acute aquatic toxicity of the 5:3 acid (EC/LC50s of 22.5 to >103 mg L(-1)) which demonstrated less aquatic toxicity than the 7:3 acid (EC/LC50s of 0.4-32 mg L(-1)). The cladoceran, D. magna and the green alga, P. subcapitata had generally similar EC50 values for a given substance while fish were typically equally or less sensitive with the exception that PFPeA was most toxic to fish. Predicted no-effect concentrations (PNECs) were estimated using approaches consistent with REACH guidance and when compared with available environmental concentrations, these PNECs suggest that the fluorinated acids tested pose little risk for aquatic organisms.