NCBI viral genomes resource

The potential of Chlorella spp. as antiviral source against African swine fever virus through a virtual screening pipeline

African swine fever (ASF) causes high mortality in pigs and threatens global swine production. There is still a lack of therapeutics available, with two vaccines under scrutiny and no approved small-molecule drugs. Eleven (11) viral proteins were used to identify potential antivirals in in silico screening of secondary metabolites (127) from Chlorella spp. The metabolites were screened for affinity and binding selectivity. High-scoring compounds were assessed through in silico ADMET (Absorption, Distribution, Metabolism, Excretion, Toxicity) predictions, compared to structurally similar drugs, and checked for off-target docking with prepared swine receptors. Molecular dynamics (MD) simulations determined binding stability while binding energy was measured in Molecular Mechanics - Generalized Born Surface Area (MMGBSA) or Poisson-Boltzmann Surface Area (MMPBSA). Only six (6) compounds passed until MD analyses, of which five (5) were stable after 100 ns of MD runs. Of these five compounds, only three had binding affinities that were comparable to or stronger than controls. Specifically, phytosterols 24,25-dihydrolanosterol and CID 4206521 that interact with the RNA capping enzyme (pNP868R), and ergosterol which bound to the Erv-like thioreductase (pB119L). The compounds identified in this study can be used as a theoretical basis for in vitro screening to develop potent antiviral drugs against ASFV.

Innovative omics strategies in fermented fruits and vegetables: Unveiling nutritional profiles, microbial diversity, and future prospects

Fermented fruits and vegetables (FFVs) are not only rich in essential nutrients but also contain distinctive flavors, prebiotics, and metabolites. Although omics techniques have gained widespread recognition as an analytical strategy for FFVs, its application still encounters several challenges due to the intricacies of biological systems. This review systematically summarizes the advances, obstacles and prospects of genomics, transcriptomics, proteomics, metabolomics, and multi-omics strategies in FFVs. It is evident that beyond traditional applications, such as the exploration of microbial diversity, protein expression, and metabolic pathways, omics techniques exhibit innovative potential in deciphering stress response mechanisms and uncovering spoilage microorganisms. The adoption of multi-omics strategies is paramount to acquire a multidimensional network fusion, thereby mitigating the limitations of single omics strategies. Although substantial progress has been made, this review underscores the necessity for a comprehensive repository of omics data and the establishment of universal databases to ensure precision in predictions. Furthermore, multidisciplinary integration with other physical or biochemical approaches is imperative, as it enriches our comprehension of this intricate process.

Identification and epidemiological analysis of a putative novel hantavirus in Australian flying foxes

In July 2017, an investigation into the cause of neurological signs in a black flying fox (Pteropus alecto, family Pteropodidae) identified a putative novel hantavirus (Robina virus, ROBV, order Bunyavirales, family Hantaviridae, genus Mobatvirus) in its brain. Analysis of the evolutionary relationship between other hantaviruses using maximum-likelihood, a systematic Bayesian clustering approach, and a minimum spanning tree, all suggest that ROBV is most closely related to another Mobatvirus, Quezon virus, previously identified in the lung of a Philippine frugivorous bat (Rousettus amplexicaudatus, also family Pteropodidae). Subsequently, between March 2018 and October 2023, a total of 495 bats were opportunistically screened for ROBV with an experimental qRT-PCR. The total prevalence of ROBV RNA detected in Pteropus spp. was 4.2% (95% CI 2.8-6.4%). Binomial modelling identified that there was substantial evidence supporting an increase (P = 0.033) in the detection of ROBV RNA in bats in 2019 and 2020 suggesting of a possible transient epidemic. There was also moderate evidence to support the effect of season (P = 0.064), with peak detection in the cooler seasons, autumn, and winter, possibly driven by physiological and ecological factors similar to those already identified for other bat-borne viruses. This is Australia's first reported putative hantavirus and its identification could expand the southern known range of hantaviruses in Australasia.

Differentiation of the Anammox core microbiome: Unraveling the evolutionary impetus of scalable gene flow

Anaerobic ammonium oxidation bacteria (AAOB), distinguished by their unique autotrophic nitrogen metabolism, hold pivotal positions in the global nitrogen cycle and environmental biotechnologies. However, the ecophysiology and evolution of AAOB remain poorly understood, attributed to the absence of monocultures. Hence, a comprehensive elucidation of the AAOB-dominated core microbiome, anammox core, is imperative to further completing the theory of engineered nitrogen removal and ecological roles of anammox. Performing taxonomic and phylogenetic analyses on collected genome repertoires, we show here that Candidatus Brocadia and Candidatus Kuenenia possesses a more compact core than Candidatus Jettenia, which partly explains why the latter has a less common ecological presence. Evidence of gene flow is particularly striking in functions related to biosynthesis and oxygen detoxification, underscoring the evolutionary forces driving lineage and core differentiation. Furthermore, CRISPR spacer traceback of the AAOB metagenome-assembled genomes (MAGs) reveals a series of genetic traces for the concealed phages. By reconceptualizing the functional divergence of AAOB with the historical role of phages, we ultimately propose a coevolutionary framework to understand the evolutionary trajectory of anammox microecology. The discoveries provided in this study offer new insights into understanding the evolution of AAOB and the ecology of anammox.

Computational Analysis of the Accumulation of Mutations in Therapeutically Important RNA Viral Proteins During Pandemics with Special Emphasis on SARS-CoV-2

Single stranded RNA viruses are primary causative agents for pandemics, causing extensive morbidity and mortality worldwide. A pivotal question in pandemic preparedness and therapeutic intervention is what are the specific mutations which are more likely to emerge during such global health crises? This study aims to identify markers for mutations with the highest probability of emergence in these pandemics, focusing on the SARS-CoV-2 spike protein, an essential and therapeutically significant viral protein, starting from sequence information from the onset of the pandemic until July 2022. Quite consistently, we observed that emerged mutations tended to demonstrate a high genetic score, which reflects high similarity of the type of codon required for translation between an amino acid and to the mutated one. Further, this pattern is also observed in therapeutically significant proteins of other ssRNA pandemic viruses, including influenza (HA, NA), spike proteins of Ebola, envelope of Dengue and Chikungunya. We propose that the genetic score serves as an initial indicator, preceding the actual impact of the mutation on viral fitness. Finally, we developed a comprehensive computational pipeline to further explore and predict the subsequent effects of mutations on viral fitness. We believe that our pipeline can narrow down and predict future mutations in therapeutically important viral proteins during a pandemic.

Genome sequences of human anelloviruses in the Lamedtorquevirus, Memtorquevirus, and Samektorquevirus genera identified from the female genital tract

We identified and characterized seven anellovirus genome sequences in the female genital tract through virome metagenomic sequencing of cervicovaginal lavage specimens from women living with HIV in Peru. Phylogenetic and genomic analyses indicate that they belong to three newly proposed Lamedtorquevirus, Memtorquevirus, and Samektorquevirus genera in the Anelloviridae family.

Computationally designed Spike antigens induce neutralising responses against the breadth of SARS-COV-2 variants

Updates of SARS-CoV-2 vaccines are required to generate immunity in the population against constantly evolving SARS-CoV-2 variants of concerns (VOCs). Here we describe three novel in-silico designed spike-based antigens capable of inducing neutralising antibodies across a spectrum of SARS-CoV-2 VOCs. Three sets of antigens utilising pre-Delta (T2_32), and post-Gamma sequence data (T2_35 and T2_36) were designed. T2_32 elicited superior neutralising responses against VOCs compared to the Wuhan-1 spike antigen in DNA prime-boost immunisation regime in guinea pigs. Heterologous boosting with the attenuated poxvirus - Modified vaccinia Ankara expressing T2_32 induced broader neutralising immune responses in all primed animals. T2_32, T2_35 and T2_36 elicited broader neutralising capacity compared to the Omicron BA.1 spike antigen administered by mRNA immunisation in mice. These findings demonstrate the utility of structure-informed computationally derived modifications of spike-based antigens for inducing broad immune responses covering more than 2 years of evolved SARS-CoV-2 variants.

SARS-CoV-2 Genotyping Highlights the Challenges in Spike Protein Drift Independent of Other Essential Proteins

As of 2024, SARS-CoV-2 continues to propagate and drift as an endemic virus, impacting healthcare for years. The largest sequencing initiative for any species was initiated to combat the virus, tracking changes over time at a full virus base-pair resolution. The SARS-CoV-2 sequencing represents a unique opportunity to understand selective pressures and viral evolution but requires cross-disciplinary approaches from epidemiology to functional protein biology. Within this work, we integrate a two-year genotyping window with structural biology to explore the selective pressures of SARS-CoV-2 on protein insights. Although genotype and the Spike (Surface Glycoprotein) protein continue to drift, most SARS-CoV-2 proteins have had few amino acid alterations. Within Spike, the high drift rate of amino acids involved in antibody evasion also corresponds to changes within the ACE2 binding pocket that have undergone multiple changes that maintain functional binding. The genotyping suggests selective pressure for receptor specificity that could also confer changes in viral risk. Mapping of amino acid changes to the structures of the SARS-CoV-2 co-transcriptional complex (nsp7-nsp14), nsp3 (papain-like protease), and nsp5 (cysteine protease) proteins suggest they remain critical factors for drug development that will be sustainable, unlike those strategies targeting Spike.

vAMPirus: A versatile amplicon processing and analysis program for studying viruses

Amplicon sequencing is an effective and increasingly applied method for studying viral communities in the environment. Here, we present vAMPirus, a user-friendly, comprehensive, and versatile DNA and RNA virus amplicon sequence analysis program, designed to support investigators in exploring virus amplicon sequencing data and running informed, reproducible analyses. vAMPirus intakes raw virus amplicon libraries and, by default, performs nucleotide- and amino acid-based analyses to produce results such as sequence abundance information, taxonomic classifications, phylogenies and community diversity metrics. The vAMPirus analytical framework leverages 16 different opensource tools and provides optional approaches that can increase the ratio of biological signal-to-noise and thereby reveal patterns that would have otherwise been masked. Here, we validate the vAMPirus analytical framework and illustrate its implementation as a general virus amplicon sequencing workflow by recapitulating findings from two previously published double-stranded DNA virus datasets. As a case study, we also apply the program to explore the diversity and distribution of a coral reef-associated RNA virus. vAMPirus is streamlined within Nextflow, offering straightforward scalability, standardization and communication of virus lineage-specific analyses. The vAMPirus framework is designed to be adaptable; community-driven analytical standards will continue to be incorporated as the field advances. vAMPirus supports researchers in revealing patterns of virus diversity and population dynamics in nature, while promoting study reproducibility and comparability.

Metatranscriptomic Sequencing of Sheath Blight-Associated Isolates of Rhizoctonia solani Revealed Multi-Infection by Diverse Groups of RNA Viruses

Rice sheath blight, caused by the soil-borne fungus Rhizoctonia solani (teleomorph: Thanatephorus cucumeris, Basidiomycota), is one of the most devastating phytopathogenic fungal diseases and causes yield loss. Here, we report on a very high prevalence (100%) of potential virus-associated double-stranded RNA (dsRNA) elements for a collection of 39 fungal strains of R. solani from the rice sheath blight samples from at least four major rice-growing areas in the Philippines and a reference isolate from the International Rice Research Institute, showing different colony phenotypes. Their dsRNA profiles suggested the presence of multiple viral infections among these Philippine R. solani populations. Using next-generation sequencing, the viral sequences of the three representative R. solani strains (Ilo-Rs-6, Tar-Rs-3, and Tar-Rs-5) from different rice-growing areas revealed the presence of at least 36 viruses or virus-like agents, with the Tar-Rs-3 strain harboring the largest number of viruses (at least 20 in total). These mycoviruses or their candidates are believed to have single-stranded RNA or dsRNA genomes and they belong to or are associated with the orders Martellivirales, Hepelivirales, Durnavirales, Cryppavirales, Ourlivirales, and Ghabrivirales based on their coding-complete RNA-dependent RNA polymerase sequences. The complete genome sequences of two novel RNA viruses belonging to the proposed family Phlegiviridae and family Mitoviridae were determined.

Discovery and characterization of BRBV-sheep virus in nasal swabs from domestic sheep in China

Introduction

The escalating occurrence of infectious disease outbreaks in humans and animals necessitates innovative, effective, and integrated research to better comprehend their transmission and dynamics. Viral infection in livestock has led to profound economic losses globally. Pneumonia is the prevalent cause of death in sheep. However, very few studies exist regarding virus-related pathogens in sheep. Metagenomics sequencing technologies in livestock research hold significant potential to elucidate these contingencies and enhance our understanding.

Methods

Therefore, this study aims to characterize respiratory viromes in paired nasal swabs from Inner Mongolian feedlot sheep in China using metaviromic sequencing. Through deep sequencing, de novo assembly, and similarity searches using translated protein sequences, several previously uncharacterized and known viruses were identified in this study.

Results

Among these discoveries, a novel Bovine Rhinitis B Virus (BRBV) (BRBV-sheep) strain was serendipitously detected in the nasal swabs of domestic sheep (Ovis aries). To facilitate further molecular epidemiological studies, the entire genome of BRBV-sheep was also determined. Owing to the unique sequence characteristics and phylogenetic position of BRBV-sheep, genetically distinct lineages of BRBV in sheep may exist. A TaqMan-based qRT-PCR assay targeting the 3D polymerase gene was developed and used to screen 592 clinical sheep specimens. The results showed that 44.59% of the samples (264/592) were positive. These findings suggest that BRBV sheep are widespread among Inner Mongolian herds.

Conclusion

This discovery marks the initial identification of BRBV in sheep within Inner Mongolia, China. These findings contribute to our understanding of the epidemiology and genetic evolution of BRBV. Recognizing the presence of BRBV in sheep informs strategies for disease management and surveillance and the potential development of targeted interventions to control its spread.

Viral genome sequencing methods: benefits and pitfalls of current approaches

Whole genome sequencing of viruses provides high-resolution molecular insights, enhancing our understanding of viral genome function and phylogeny. Beyond fundamental research, viral sequencing is increasingly vital for pathogen surveillance, epidemiology, and clinical applications. As sequencing methods rapidly evolve, the diversity of viral genomics applications and catalogued genomes continues to expand. Advances in long-read, single molecule, real-time sequencing methodologies present opportunities to sequence contiguous, haplotype resolved viral genomes in a range of research and applied settings. Here we present an overview of nucleic acid sequencing methods and their applications in studying viral genomes. We emphasise the advantages of different viral sequencing approaches, with a particular focus on the benefits of third-generation sequencing technologies in elucidating viral evolution, transmission networks, and pathogenesis.

Characterisation of ten NS2B-NS3 proteases: Paving the way for pan-flavivirus drugs

Flaviviruses can cause severe illness in humans. Effective and safe vaccines are available for some species; however, for many flaviviruses disease prevention or specific treatments remain unavailable. The viral replication cycle depends on the proteolytic activity of the NS2B-NS3 protease, which releases functional viral proteins from a non-functional polyprotein precursor, rendering the protease a promising drug target. In this study, we characterised recombinant NS2B-NS3 proteases from ten flaviviruses including three unreported proteases from the Usutu, Kyasanur forest disease and Powassan viruses. All protease constructs comprise a covalent Gly4-Ser-Gly4 linker connecting the NS3 serine protease domain with its cofactor NS2B. We conducted a comprehensive cleavage site analysis revealing areas of high conversion. While all proteases were active in enzymatic assays, we noted a 1000-fold difference in catalytic efficiency across proteases from different flaviviruses. Two bicyclic peptide inhibitors displayed anti-pan-flaviviral protease activity with inhibition constants ranging from 10 to 1000 nM.

Limited potexvirus diversity in eastern Gulf of Mexico seagrass meadows

Turtlegrass virus X, which infects the seagrass Thalassia testudinum, is the only potexvirus known to infect marine flowering plants. We investigated potexvirus distribution in seagrasses using a degenerate reverse transcription polymerase chain reaction (RT-PCR) assay originally designed to capture potexvirus diversity in terrestrial plants. The assay, which implements Potex-5 and Potex-2RC primers, successfully amplified a 584 nt RNA-dependent RNA polymerase (RdRp) fragment from TVX-infected seagrasses. Following validation, we screened 74 opportunistically collected, apparently healthy seagrass samples for potexviruses using this RT-PCR assay. The survey examined the host species T. testudinum, Halodule wrightii, Halophila stipulacea, Syringodium filiforme, Ruppia maritima, and Zostera marina. Potexvirus PCR products were successfully generated only from T. testudinum samples and phylogenetic analysis of sequenced PCR products revealed five distinct TVX sequence variants. Although the RT-PCR assay revealed limited potexvirus diversity in seagrasses, the expanded geographic distribution of TVX shown here emphasizes the importance of future studies to investigate T. testudinum populations across its native range and understand how the observed fine-scale genetic diversity affects host-virus interactions.

Stabilized trimeric peptide immunogens of the complete HIV-1 gp41 N-heptad repeat and their use as HIV-1 vaccine candidates

Efforts to develop an HIV-1 vaccine include those focusing on conserved structural elements as the target of broadly neutralizing monoclonal antibodies. MAb D5 binds to a highly conserved hydrophobic pocket on the gp41 N-heptad repeat (NHR) coiled coil and neutralizes through prevention of viral fusion and entry. Assessment of 17-mer and 36-mer NHR peptides presenting the D5 epitope in rodent immunogenicity studies showed that the longer peptide elicited higher titers of neutralizing antibodies, suggesting that neutralizing epitopes outside of the D5 pocket may exist. Although the magnitude and breadth of neutralization elicited by NHR-targeting antigens are lower than that observed for antibodies directed to other epitopes on the envelope glycoprotein complex, it has been shown that NHR-directed antibodies are potentiated in TZM-bl cells containing the FcγRI receptor. Herein, we report the design and evaluation of covalently stabilized trimeric 51-mer peptides encompassing the complete gp41 NHR. We demonstrate that these peptide trimers function as effective antiviral entry inhibitors and retain the ability to present the D5 epitope. We further demonstrate in rodent and nonhuman primate immunization studies that our 51-mer constructs elicit a broader repertoire of neutralizing antibody and improved cross-clade neutralization of primary HIV-1 isolates relative to 17-mer and 36-mer NHR peptides in A3R5 and FcγR1-enhanced TZM-bl assays. These results demonstrate that sensitive neutralization assays can be used for structural enhancement of moderately potent neutralizing epitopes. Finally, we present expanded trimeric peptide designs which include unique low-molecular-weight scaffolds that provide versatility in our immunogen presentation strategy.

A novel approach to exploring the dark genome and its application to mapping of the vertebrate virus fossil record

BACKGROUND

Genomic regions that remain poorly understood, often referred to as the dark genome, contain a variety of functionally relevant and biologically informative features. These include endogenous viral elements (EVEs)-virus-derived sequences that can dramatically impact host biology and serve as a virus fossil record. In this study, we introduce a database-integrated genome screening (DIGS) approach to investigate the dark genome in silico, focusing on EVEs found within vertebrate genomes.

RESULTS

Using DIGS on 874 vertebrate genomes, we uncover approximately 1.1 million EVE sequences, with over 99% originating from endogenous retroviruses or transposable elements that contain EVE DNA. We show that the remaining 6038 sequences represent over a thousand distinct horizontal gene transfer events across 10 virus families, including some that have not previously been reported as EVEs. We explore the genomic and phylogenetic characteristics of non-retroviral EVEs and determine their rates of acquisition during vertebrate evolution. Our study uncovers novel virus diversity, broadens knowledge of virus distribution among vertebrate hosts, and provides new insights into the ecology and evolution of vertebrate viruses.

CONCLUSIONS

We comprehensively catalog and analyze EVEs within 874 vertebrate genomes, shedding light on the distribution, diversity, and long-term evolution of viruses and reveal their extensive impact on vertebrate genome evolution. Our results demonstrate the power of linking a relational database management system to a similarity search-based screening pipeline for in silico exploration of the dark genome.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Tying the knot: Unraveling the intricacies of the coronavirus frameshift pseudoknot

Understanding and targeting functional RNA structures towards treatment of coronavirus infection can help us to prepare for novel variants of SARS-CoV-2 (the virus causing COVID-19), and any other coronaviruses that could emerge via human-to-human transmission or potential zoonotic (inter-species) events. Leveraging the fact that all coronaviruses use a mechanism known as -1 programmed ribosomal frameshifting (-1 PRF) to replicate, we apply algorithms to predict the most energetically favourable secondary structures (each nucleotide involved in at most one pairing) that may be involved in regulating the -1 PRF event in coronaviruses, especially SARS-CoV-2. We compute previously unknown most stable structure predictions for the frameshift site of coronaviruses via hierarchical folding, a biologically motivated framework where initial non-crossing structure folds first, followed by subsequent, possibly crossing (pseudoknotted), structures. Using mutual information from 181 coronavirus sequences, in conjunction with the algorithm KnotAli, we compute secondary structure predictions for the frameshift site of different coronaviruses. We then utilize the Shapify algorithm to obtain most stable SARS-CoV-2 secondary structure predictions guided by frameshift sequence-specific and genome-wide experimental data. We build on our previous secondary structure investigation of the singular SARS-CoV-2 68 nt frameshift element sequence, by using Shapify to obtain predictions for 132 extended sequences and including covariation information. Previous investigations have not applied hierarchical folding to extended length SARS-CoV-2 frameshift sequences. By doing so, we simulate the effects of ribosome interaction with the frameshift site, providing insight to biological function. We contribute in-depth discussion to contextualize secondary structure dual-graph motifs for SARS-CoV-2, highlighting the energetic stability of the previously identified 3_8 motif alongside the known dominant 3_3 and 3_6 (native-type) -1 PRF structures. Using a combination of thermodynamic methods and sequence covariation, our novel predictions suggest function of the attenuator hairpin via previously unknown pseudoknotted base pairing. While certain initial RNA folding is consistent, other pseudoknotted base pairs form which indicate potential conformational switching between the two structures.

DengueSeq: a pan-serotype whole genome amplicon sequencing protocol for dengue virus

BACKGROUND

The increasing burden of dengue virus on public health due to more explosive and frequent outbreaks highlights the need for improved surveillance and control. Genomic surveillance of dengue virus not only provides important insights into the emergence and spread of genetically diverse serotypes and genotypes, but it is also critical to monitor the effectiveness of newly implemented control strategies. Here, we present DengueSeq, an amplicon sequencing protocol, which enables whole-genome sequencing of all four dengue virus serotypes.

RESULTS

We developed primer schemes for the four dengue virus serotypes, which can be combined into a pan-serotype approach. We validated both approaches using genetically diverse virus stocks and clinical specimens that contained a range of virus copies. High genome coverage (>95%) was achieved for all genotypes, except DENV2 (genotype VI) and DENV 4 (genotype IV) sylvatics, with similar performance of the serotype-specific and pan-serotype approaches. The limit of detection to reach 70% coverage was 10-100 RNA copies/μL for all four serotypes, which is similar to other commonly used primer schemes. DengueSeq facilitates the sequencing of samples without known serotypes, allows the detection of multiple serotypes in the same sample, and can be used with a variety of library prep kits and sequencing instruments.

CONCLUSIONS

DengueSeq was systematically evaluated with virus stocks and clinical specimens spanning the genetic diversity within each of the four dengue virus serotypes. The primer schemes can be plugged into existing amplicon sequencing workflows to facilitate the global need for expanded dengue virus genomic surveillance.

Biogeographic patterns and drivers of soil viromes

Viruses are crucial in shaping soil microbial functions and ecosystems. However, studies on soil viromes have been limited in both spatial scale and biome coverage. Here we present a comprehensive synthesis of soil virome biogeographic patterns using the Global Soil Virome dataset (GSV) wherein we analysed 1,824 soil metagenomes worldwide, uncovering 80,750 partial genomes of DNA viruses, 96.7% of which are taxonomically unassigned. The biogeography of soil viral diversity and community structure varies across different biomes. Interestingly, the diversity of viruses does not align with microbial diversity and contrasts with it by showing low diversity in forest and shrubland soils. Soil texture and moisture conditions are further corroborated as key factors affecting diversity by our predicted soil viral diversity atlas, revealing higher diversity in humid and subhumid regions. In addition, the binomial degree distribution pattern suggests a random co-occurrence pattern of soil viruses. These findings are essential for elucidating soil viral ecology and for the comprehensive incorporation of viruses into soil ecosystem models.

Coding genomes with gapped pattern graph convolutional network

MOTIVATION

Genome sequencing technologies reveal a huge amount of genomic sequences. Neural network-based methods can be prime candidates for retrieving insights from these sequences because of their applicability to large and diverse datasets. However, the highly variable lengths of genome sequences severely impair the presentation of sequences as input to the neural network. Genetic variations further complicate tasks that involve sequence comparison or alignment.

RESULTS

Inspired by the theory and applications of "spaced seeds," we propose a graph representation of genome sequences called "gapped pattern graph." These graphs can be transformed through a Graph Convolutional Network to form lower-dimensional embeddings for downstream tasks. On the basis of the gapped pattern graphs, we implemented a neural network model and demonstrated its performance on diverse tasks involving microbe and mammalian genome data. Our method consistently outperformed all the other state-of-the-art methods across various metrics on all tasks, especially for the sequences with limited homology to the training data. In addition, our model was able to identify distinct gapped pattern signatures from the sequences.

AVAILABILITY AND IMPLEMENTATION

The framework is available at https://github.com/deepomicslab/GCNFrame.

Molecular and phylogenetic characterization of the monkeypox outbreak in the South of Spain

Background and Aim

Until the May 2022 Monkeypox (MPXV) outbreak, which spread rapidly to many non-endemic countries, the virus was considered a viral zoonosis limited to some African countries. The Andalusian circuit of genomic surveillance was rapidly applied to characterize the MPXV outbreak in the South of Spain.

Methods

Whole genome sequencing was used to obtain the genomic profiles of samples collected across the south of Spain, representative of all the provinces of Andalusia. Phylogenetic analysis was used to study the relationship of the isolates and the available sequences of the 2022 outbreak.

Results

Whole genome sequencing of a total of 160 MPXV viruses from the different provinces that reported cases were obtained. Interestingly, we report the sequences of MPXV viruses obtained from two patients who died. While one of the isolates bore no noteworthy mutations that explain a potential heightened virulence, in another patient the second consecutive genome sequence, performed after the administration of tecovirimat, uncovered a mutation within the A0A7H0DN30 gene, known to be a prime target for tecovirimat in its Vaccinia counterpart. In general, a low number of mutations were observed in the sequences reported, which were very similar to the reference of the 2022 outbreak (OX044336), as expected from a DNA virus. The samples likely correspond to several introductions of the circulating MPXV viruses from the last outbreak. The virus sequenced from one of the two patients that died presented a mutation in a gene that bears potential connections to drug resistance. This mutation was absent in the initial sequencing before treatment.

Exploring the gut DNA virome in fecal immunochemical test stool samples reveals associations with lifestyle in a large population-based study

Stool samples for fecal immunochemical tests (FIT) are collected in large numbers worldwide as part of colorectal cancer screening programs. Employing FIT samples from 1034 CRCbiome participants, recruited from a Norwegian colorectal cancer screening study, we identify, annotate and characterize more than 18000 DNA viruses, using shotgun metagenome sequencing. Only six percent of them are assigned to a known taxonomic family, with Microviridae being the most prevalent viral family. Linking individual profiles to comprehensive lifestyle and demographic data shows 17/25 of the variables to be associated with the gut virome. Physical activity, smoking, and dietary fiber consumption exhibit strong and consistent associations with both diversity and relative abundance of individual viruses, as well as with enrichment for auxiliary metabolic genes. We demonstrate the suitability of FIT samples for virome analysis, opening an opportunity for large-scale studies of this enigmatic part of the gut microbiome. The diverse viral populations and their connections to the individual lifestyle uncovered herein paves the way for further exploration of the role of the gut virome in health and disease.

Discovering and exploring the hidden diversity of human gut viruses using highly enriched virome samples

Viruses are an abundant and crucial component of the human microbiome, but accurately discovering them via metagenomics is still challenging. Currently, the available viral reference genomes poorly represent the diversity in microbiome samples, and expanding such a set of viral references is difficult. As a result, many viruses are still undetectable through metagenomics even when considering the power of de novo metagenomic assembly and binning, as viruses lack universal markers. Here, we describe a novel approach to catalog new viral members of the human gut microbiome and show how the resulting resource improves metagenomic analyses. We retrieved >3,000 viral-like particles (VLP) enriched metagenomic samples (viromes), evaluated the efficiency of the enrichment in each sample to leverage the viromes of highest purity, and applied multiple analysis steps involving assembly and comparison with hundreds of thousands of metagenome-assembled genomes to discover new viral genomes. We reported over 162,000 viral sequences passing quality control from thousands of gut metagenomes and viromes. The great majority of the retrieved viral sequences (~94.4%) were of unknown origin, most had a CRISPR spacer matching host bacteria, and four of them could be detected in >50% of a set of 18,756 gut metagenomes we surveyed. We included the obtained collection of sequences in a new MetaPhlAn 4.1 release, which can quantify reads within a metagenome matching the known and newly uncovered viral diversity. Additionally, we released the viral database for further virome and metagenomic studies of the human microbiome.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 Spike (S) protein. With this approach, we first identified Candidate Adaptive Polymorphisms (CAPs) of the SARS-CoV-2 Spike protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Bioinformatic analysis of defective viral genomes in SARS-CoV-2 and its impact on population infection characteristics

DVGs (Defective Viral Genomes) are prevalent in RNA virus infections. In this investigation, we conducted an analysis of high-throughput sequencing data and observed widespread presence of DVGs in SARS-CoV-2. Comparative analysis between SARS-CoV-2 and diverse DNA viruses revealed heightened susceptibility to damage and increased sequencing sample heterogeneity within the SARS-CoV-2 genome. Whole-genome sequencing depth variability analysis exhibited a higher coefficient of variation for SARS-CoV-2, while DVG analysis indicated a significant proportion of recombination sites, signifying notable genome heterogeneity and suggesting that a large proportion of assembled virus particles contain incomplete RNA sequences. Moreover, our investigation explored the sequencing depth and DVG content differences among various strains. Our findings revealed that as the virus evolves, there is a notable increase in the proportion of intact genomes within virus particles, as evidenced by third-generation sequencing data. Specifically, the proportion of intact genome in the Omicron strain surpassed that of the Delta and Alpha strains. This observation effectively elucidates the heightened infectiousness of the Omicron strain compared to the Delta and Alpha strains. We also postulate that this improvement in completeness stems from enhanced virus assembly capacity, as the Omicron strain can promptly facilitate the binding of RNA and capsid protein, thereby reducing the exposure time of vulnerable virus RNA in the host environment and significantly mitigating its degradation. Finally, employing mathematical modeling, we simulated the impact of DVG effects under varying environmental factors on infection characteristics and population evolution. Our findings provide an explanation for the close association between symptom severity and the extent of virus invasion, as well as the substantial disparity in population infection characteristics caused by the same strain under distinct environmental conditions. This study presents a novel approach for future virus research and vaccine development.

Metagenomic analyses of 7000 to 5500 years old coprolites excavated from the Torihama shell-mound site in the Japanese archipelago

Coprolites contain various kinds of ancient DNAs derived from gut micro-organisms, viruses, and foods, which can help to determine the gut environment of ancient peoples. Their genomic information should be helpful in elucidating the interaction between hosts and microbes for thousands of years, as well as characterizing the dietary behaviors of ancient people. We performed shotgun metagenomic sequencing on four coprolites excavated from the Torihama shell-mound site in the Japanese archipelago. The coprolites were found in the layers of the Early Jomon period, corresponding stratigraphically to 7000 to 5500 years ago. After shotgun sequencing, we found that a significant number of reads showed homology with known gut microbe, viruses, and food genomes typically found in the feces of modern humans. We detected reads derived from several types of phages and their host bacteria simultaneously, suggesting the coexistence of viruses and their hosts. The food genomes provide biological evidence for the dietary behavior of the Jomon people, consistent with previous archaeological findings. These results indicate that ancient genomic analysis of coprolites is useful for understanding the gut environment and lifestyle of ancient peoples.

In silico design of a broad-spectrum multiepitope vaccine against influenza virus

Influenza remains a global health concern due to its potential to cause pandemics as a result of rapidly mutating influenza virus strains. Existing vaccines often struggle to keep up with these rapidly mutating flu viruses. Therefore, the development of a broad-spectrum peptide vaccine that can stimulate an optimal antibody response has emerged as an innovative approach to addressing the influenza threat. In this study, an immunoinformatic approach was employed to rapidly predict immunodominant epitopes from different antigens, aiming to develop an effective multiepitope influenza vaccine (MEV). The immunodominant B-cell linear epitopes of seasonal influenza strains hemagglutinin (HA) and neuraminidase (NA) were predicted using an antibody-peptide microarray, involving a human cohort including vaccinees and infected patients. On the other hand, bioinformatics tools were used to predict immunodominant cytotoxic T-cell (CTL) and helper T-cell (HTL) epitopes. Subsequently, these epitopes were evaluated by various immunoinformatic tools. Epitopes with high antigenicity, high immunogenicity, non-allergenicity, non-toxicity, as well as exemplary conservation were then connected in series with appropriate linkers and adjuvants to construct a broad-spectrum MEV. Moreover, the structural analysis revealed that the MEV candidates exhibited good stability, and the docking results demonstrated their strong affinity to Toll-like receptors 4 (TLR4). In addition, molecular dynamics simulation confirmed the stable interaction between TLR4 and MEVs. Three injections with MEVs showed a high level of B-cell and T-cell immune responses according to the immunological simulations in silico. Furthermore, in-silico cloning was performed, and the results indicated that the MEVs could be produced in considerable quantities in Escherichia coli (E. coli). Based on these findings, it is reasonable to create a broad-spectrum MEV against different subtypes of influenza A and B viruses in silico.

Bacteriophage Taxonomy: A Continually Evolving Discipline

While taxonomy is an often underappreciated branch of science, it serves very important roles. Bacteriophage taxonomy has evolved from a discipline based mainly on morphology, characterized by the work of David Bradley and Hans-Wolfgang Ackermann, to the sequence-based approach that is taken today. The Bacterial Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) takes a holistic approach to classifying prokaryote viruses by measuring overall DNA and protein similarity and phylogeny before making decisions about the taxonomic position of a new virus. The huge number of complete genomes being deposited with the National Center for Biotechnology Information (NCBI) and other public databases has resulted in a reassessment of the taxonomy of many viruses, and the future will see the introduction of new viral families and higher orders.

Bioinformatic Approaches for Comparative Analysis of Viruses

The field of viral genomic studies has experienced an unprecedented increase in data volume. New strains of known viruses are constantly being added to the GenBank database and so are completely new species with little or no resemblance to our databases of sequences. In addition to this, metagenomic techniques have the potential to further increase the number and rate of sequenced genomes. Besides, it is important to consider that viruses have a set of unique features that often break down molecular biology dogmas, e.g., the flux of information from RNA to DNA in retroviruses and the use of RNA molecules as genomes. As a result, extracting meaningful information from viral genomes remains a challenge and standard methods for comparing the unknown and our databases of characterized sequences may need adaptations. Thus, several bioinformatic approaches and tools have been created to address the challenge of analyzing viral data. This chapter offers descriptions and protocols of some of the most important bioinformatic techniques for comparative analysis of viruses. The authors also provide comments and discussion on how viruses' unique features can affect standard analyses and how to overcome some of the major sources of problems. Protocols and topics emphasize online tools (which are more accessible to users) and give the real experience of what most bioinformaticians do in day-by-day work with command-line pipelines. The topics discussed include (1) clustering related genomes, (2) whole genome multiple sequence alignments for small RNA viruses, (3) protein alignment for marker genes and species affiliation, (4) variant calling and annotation, and (5) virome analyses and pathogen identification.

With a pinch of salt: metagenomic insights into Namib Desert salt pan microbial mats and halites reveal functionally adapted and competitive communities

IMPORTANCE

The hyperarid Namib Desert is one of the oldest deserts on Earth. It contains multiple clusters of playas which are saline-rich springs surrounded by halite evaporites. Playas are of great ecological importance, and their indigenous (poly)extremophilic microorganisms are potentially involved in the precipitation of minerals such as carbonates and sulfates and have been of great biotechnological importance. While there has been a considerable amount of microbial ecology research performed on various Namib Desert edaphic microbiomes, little is known about the microbial communities inhabiting its multiple playas. In this work, we provide a comprehensive taxonomic and functional potential characterization of the microbial, including viral, communities of sediment mats and halites from two distant salt pans of the Namib Desert, contributing toward a better understanding of the ecology of this biome.

Viral diversity and dynamics and CRISPR-Cas-mediated immunity in a robust alkaliphilic cyanobacterial consortium

IMPORTANCE

Biotechnology applications utilizing the function of microbial communities have become increasingly important solutions as we strive for sustainable applications. Although viral infections are known to have a significant impact on microbial turnover and nutrient cycling, viral dynamics have remained largely overlooked in these engineered communities. Predatory perturbations to the functional stability of these microbial biotechnology applications must be investigated in order to design more robust applications. In this study, we closely examine virus-microbe dynamics in a model microbial community used in a biotechnology application. Our findings suggest that viral dynamics change significantly with environmental conditions and that microbial immunity may play an important role in maintaining functional stability. We present this study as a comprehensive template for other researchers interested in exploring predatory dynamics in engineered microbial communities.

Hepatitis B Virus and microRNAs: A Bioinformatics Approach

In recent decades, microRNAs (miRNAs) have emerged as key regulators of gene expression, and the identification of viral miRNAs (v-miRNAs) within some viruses, including hepatitis B virus (HBV), has attracted significant attention. HBV infections often progress to chronic states (CHB) and may induce fibrosis/cirrhosis and hepatocellular carcinoma (HCC). The presence of HBV can dysregulate host miRNA expression, influencing several biological pathways, such as apoptosis, innate and immune response, viral replication, and pathogenesis. Consequently, miRNAs are considered a promising biomarker for diagnostic, prognostic, and treatment response. The dynamics of miRNAs during HBV infection are multifaceted, influenced by host variability and miRNA interactions. Given the ability of miRNAs to target multiple messenger RNA (mRNA), understanding the viral-host (human) interplay is complex but essential to develop novel clinical applications. Therefore, bioinformatics can help to analyze, identify, and interpret a vast amount of miRNA data. This review explores the bioinformatics tools available for viral and host miRNA research. Moreover, we introduce a brief overview focusing on the role of miRNAs during HBV infection. In this way, this review aims to help the selection of the most appropriate bioinformatics tools based on requirements and research goals.

Putting everything in its place: using the INSDC compliant Pathogen Data Object Model to better structure genomic data submitted for public health applications

Fast, efficient public health actions require well-organized and coordinated systems that can supply timely and accurate knowledge. Public databases of pathogen genomic data, such as the International Nucleotide Sequence Database Collaboration (INSDC), have become essential tools for efficient public health decisions. However, these international resources began primarily for academic purposes, rather than for surveillance or interventions. Now, queries need to access not only the whole genomes of multiple pathogens but also make connections using robust contextual metadata to identify issues of public health relevance. Databases that over time developed a patchwork of submission formats and requirements need to be consistently organized and coordinated internationally to allow effective searches.To help resolve these issues, we propose a common pathogen data structure called the Pathogen Data Object Model (DOM) that will formalize the minimum pieces of sequence data and contextual data necessary for general public health uses, while recognizing that submitters will likely withhold a wide range of non-public contextual data. Further, we propose contributors use the Pathogen DOM for all pathogen submissions (bacterial, viral, fungal, and parasites), which will simplify data submissions and provide a consistent and transparent data structure for downstream data analyses. We also highlight how improved submission tools can support the Pathogen DOM, offering users additional easy-to-use methods to ensure this structure is followed.

Emergence of a new designated clade 16 with significant antigenic drift in hemagglutinin gene of H9N2 subtype avian influenza virus in eastern China

H9N2 avian influenza viruses (AIVs) pose an increasing threat to the poultry industry worldwide and have pandemic potential. Vaccination has been principal prevention strategy to control H9N2 in China since 1998, but vaccine effectiveness is persistently challenged by the emergence of the genetic and/or antigenic variants. Here, we analysed the genetic and antigenic characteristics of H9N2 viruses in China, including 70 HA sequences of H9N2 isolates from poultry, 7358 from online databases during 2010-2020, and 15 from the early reference strains. Bayesian analyses based on hemagglutinin (HA) gene revealed that a new designated clade16 emerged in April 2012, and was prevalent and co-circulated with clade 15 since 2013 in China. Clade 16 viruses exhibited decreased cross-reactivity with those from clade 15. Antigenic Cartography analyses showed represent strains were classified into three antigenic groups named as Group1, Group2 and Group3, and most of the strains in Group 3 (15/17, 88.2%) were from Clade 16 while most of the strains in Group2 (26/29, 89.7%) were from Clade 15. The mean distance between Group 3 and Group 2 was 4.079 (95%CI 3.605-4.554), revealing that major switches to antigenic properties were observed over the emergence of clade 16. Genetic analysis indicated that 11 coevolving amino acid substitutions primarily at antigenic sites were associated with the antigenic differences between clade 15 and clade 16. These data highlight complexities of the genetic evolution and provide a framework for the genetic basis and antigenic characterization of emerging clade 16 of H9N2 subtype avian influenza virus.

An In Silico Analysis of PCR-Based Monkeypox Virus Detection Assays: A Case Study for Ongoing Clinical Surveillance

The 2022 global Mpox outbreak swiftly introduced unforeseen diversity in the monkeypox virus (MPXV) population, resulting in numerous Clade IIb sublineages. This propagation of new MPXV mutations warrants the thorough re-investigation of previously recommended or validated primers designed to target MPXV genomes. In this study, we explored 18 PCR primer sets and examined their binding specificity against 5210 MPXV genomes, representing all the established MPXV lineages. Our results indicated that only five primer sets resulted in almost all perfect matches against the targeted MPXV lineages, and the remaining primer sets all contained 1-2 mismatches against almost all the MPXV lineages. We further investigated the mismatched primer-genome pairs and discovered that some of the primers overlapped with poorly sequenced and assembled regions of the MPXV genomes, which are consistent across multiple lineages. However, we identified 173 99% genome-wide conserved regions across all 5210 MPXV genomes, representing 30 lineages/clades with at least 80% lineage-specific consensus for future primer development and primer binding evaluation. This exercise is crucial to ensure that the current detection schemes are robust and serve as a framework for primer evaluation in clinical testing development for other infectious diseases.

In silico study identifies peptide inhibitors that negate the effect of non-synonymous mutations in major drug targets of SARS-CoV-2 variants

Since its advent in December 2019, SARS-CoV-2 has diverged into multiple variants with differing levels of virulence owing to the accumulation of mutations in its genome. The structural changes induced by non-synonymous mutations in major drug targets of the virus are known to alter the binding of potential antagonistic inhibitors. Here, we analyzed the effects of non-synonymous mutations in major targets of SARS-CoV-2 in response to potential peptide inhibitors. We screened 12 peptides reported to have anti-viral properties against RBD and 5 peptides against Mpro of SARS-CoV-2 variants using molecular docking and simulation approaches. The mutational landscape of RBD among SARS-CoV-2 variants had 21 non-synonymous mutations across 18 distinct sites. Among these, 14 mutations were present in the RBM region directly interacting with the hACE2 receptor. However, Only 3 non-synonymous mutations were observed in Mpro. We found that LCB1 - a de novo-synthesized peptide has the highest binding affinity to RBD despite non-synonymous mutations in variants and engages key residues of RBD-hACE2 interaction such as K417, E484, N487, and N501. Similarly, an antimicrobial peptide; 2JOS, was identified against Mpro with high binding affinity as it interacts with key residues in dimerization sites such as E166 and F140 crucial for viral replication. MD simulations affirm the stability of RBD-LCB1 and Mpro-2JOS complexes with an average RMSD of 1.902 and 2.476 respectively. We ascertain that LCB1 and 2JOS peptides are promising inhibitors to combat emerging variants of SARS-CoV-2 and thus warrant further investigations using in-vitro and in-vivo analysis.Communicated by Ramaswamy H. Sarma.

Open questions in the social lives of viruses

Social interactions among viruses occur whenever multiple viral genomes infect the same cells, hosts, or populations of hosts. Viral social interactions range from cooperation to conflict, occur throughout the viral world, and affect every stage of the viral lifecycle. The ubiquity of these social interactions means that they can determine the population dynamics, evolutionary trajectory, and clinical progression of viral infections. At the same time, social interactions in viruses raise new questions for evolutionary theory, providing opportunities to test and extend existing frameworks within social evolution. Many opportunities exist at this interface: Insights into the evolution of viral social interactions have immediate implications for our understanding of the fundamental biology and clinical manifestation of viral diseases. However, these opportunities are currently limited because evolutionary biologists only rarely study social evolution in viruses. Here, we bridge this gap by (1) summarizing the ways in which viruses can interact socially, including consequences for social evolution and evolvability; (2) outlining some open questions raised by viruses that could challenge concepts within social evolution theory; and (3) providing some illustrative examples, data sources, and conceptual questions, for studying the natural history of social viruses.

Endogenous Viral Elements in Ixodid Tick Genomes

The documentation of endogenous viral elements (EVEs; virus-derived genetic material integrated into the genome of a nonviral host) has offered insights into how arthropods respond to viral infection via RNA interference pathways. Small non-coding RNAs derived from EVE loci serve to direct RNAi pathways in limiting replication and infection from cognate viruses, thus benefiting the host's fitness and, potentially, vectorial capacity. Here we use informatic approaches to analyze nine available genome sequences of hard ticks (Acari: Ixodidae; Rhipicephalus sanguineus, R. microplus, R. annulatus, Ixodes ricinus, I. persulcatus, I. scapularis, Hyalomma asiaticum, Haemaphysalis longicornis, and Dermacentor silvarum) to identify endogenous viral elements and to illustrate the shared ancestry of all elements identified. Our results highlight a broad diversity of viral taxa as having given rise to 1234 identified EVEs in ticks, with Mononegavirales (specifically Rhabdoviridae) well-represented in this subset of hard ticks. Further investigation revealed extensive adintovirus integrations in several Ixodes species, the prevalence of Bunyavirales EVEs (notably not observed in mosquitoes), and the presence of several elements similar to known emerging human and veterinary pathogens. These results will inform subsequent work on current and past associations with tick species with regard to the viruses from which their "viral fossils" are derived and may serve as a reference for quality control of various tick-omics data that may suffer from misidentification of EVEs as viral genetic material.

DengueSeq: A pan-serotype whole genome amplicon sequencing protocol for dengue virus

Background

The increasing burden of dengue virus on public health due to more explosive and frequent outbreaks highlights the need for improved surveillance and control. Genomic surveillance of dengue virus not only provides important insights into the emergence and spread of genetically diverse serotypes and genotypes, but it is also critical to monitor the effectiveness of newly implemented control strategies. Here, we present DengueSeq, an amplicon sequencing protocol, which enables whole-genome sequencing of all four dengue virus serotypes.

Results

We developed primer schemes for the four dengue virus serotypes, which can be combined into a pan-serotype approach. We validated both approaches using genetically diverse virus stocks and clinical specimens that contained a range of virus copies. High genome coverage (>95%) was achieved for all genotypes, except DENV2 (genotype VI) and DENV 4 (genotype IV) sylvatics, with similar performance of the serotype-specific and pan-serotype approaches. The limit of detection to reach 70% coverage was 101-102 RNA copies/μL for all four serotypes, which is similar to other commonly used primer schemes. DengueSeq facilitates the sequencing of samples without known serotypes, allows the detection of multiple serotypes in the same sample, and can be used with a variety of library prep kits and sequencing instruments.

Conclusions

DengueSeq was systematically evaluated with virus stocks and clinical specimens spanning the genetic diversity within each of the four dengue virus serotypes. The primer schemes can be plugged into existing amplicon sequencing workflows to facilitate the global need for expanded dengue virus genomic surveillance.

RCoV19: A One-stop Hub for SARS-CoV-2 Genome Data Integration, Variant Monitoring, and Risk Pre-warning

The Resource for Coronavirus 2019 (RCoV19) is an open-access information resource dedicated to providing valuable data on the genomes, mutations, and variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In this updated implementation of RCoV19, we have made significant improvements and advancements over the previous version. Firstly, we have implemented a highly refined genome data curation model. This model now features an automated integration pipeline and optimized curation rules, enabling efficient daily updates of data in RCoV19. Secondly, we have developed a global and regional lineage evolution monitoring platform, alongside an outbreak risk pre-warning system. These additions provide a comprehensive understanding of SARS-CoV-2 evolution and transmission patterns, enabling better preparedness and response strategies. Thirdly, we have developed a powerful interactive mutation spectrum comparison module. This module allows users to compare and analyze mutation patterns, assisting in the detection of potential new lineages. Furthermore, we have incorporated a comprehensive knowledgebase on mutation effects. This knowledgebase serves as a valuable resource for retrieving information on the functional implications of specific mutations. In summary, RCoV19 serves as a vital scientific resource, providing access to valuable data, relevant information, and technical support in the global fight against COVID-19. The complete contents of RCoV19 are available to the public at https://ngdc.cncb.ac.cn/ncov/.

A conserved class of viral RNA structures regulate translation reinitiation through dynamic ribosome interactions

Certain viral RNAs encode proteins downstream of the main protein coding region, expressed through "termination-reinitiation" events, dependent on RNA structure. RNA elements located upstream of the first stop codon within these viral mRNAs bind the ribosome, preventing ribosome recycling and inducing reinitiation. We used bioinformatic methods to identify new examples of viral reinitiation-stimulating RNAs and experimentally verified their secondary structure and function. We determined the structure of a representative viral RNA-ribosome complex using cryoEM. 3D classification and variability analyses reveal that the viral RNA structure can sample a range of conformations while remaining tethered to the ribosome, which enabling the ribosome to find a reinitiation start site within a limited range of mRNA sequence. Evaluating the conserved features and constraints of this entire RNA class in the context of the cryoEM reconstruction provides insight into mechanisms enabling reinitiation, a translation regulation strategy employed by many other viral and eukaryotic systems.

Host specific adaptations of Ligilactobacillus aviarius to poultry

The genus Ligilactobacillus encompasses species adapted to vertebrate hosts and fermented food. Their genomes encode adaptations to the host lifestyle. Reports of gut microbiota from chicken and turkey gastrointestinal tract have shown a high persistence of Ligilactobacillus aviarius along the digestive system compared to other species found in the same host. However, its adaptations to poultry as a host has not yet been described. In this work, the pan-genome of Ligilactobacillus aviarius was explored to describe the functional adaptability to the gastrointestinal environment. The core genome is composed of 1179 gene clusters that are present at least in one copy that codifies to structural, ribosomal and biogenesis proteins. The rest of the identified regions were classified into three different functional clusters of orthologous groups (clusters) that codify carbohydrate metabolism, envelope biogenesis, viral defence mechanisms, and mobilome inclusions. The pan-genome of Ligilactobacillus aviarius is a closed pan-genome, frequently found in poultry and highly prevalent across chicken faecal samples. The genome of L. aviarius codifies different clusters of glycoside hydrolases and glycosyltransferases that mediate interactions with the host cells. Accessory features, such as antiviral mechanisms and prophage inclusions, variate amongst strains from different GIT sections. This information provides hints about the interaction of this species with viral particles and other bacterial species. This work highlights functional adaptability traits present in L. aviarius that make it a dominant key member of the poultry gut microbiota and enlightens the convergent ecological relation of this species to the poultry gut environment.

Genes, inflammatory response, tolerance, and resistance to virus infections in migratory birds, bats, and rodents

Normally, the host immunological response to viral infection is coordinated to restore homeostasis and protect the individual from possible tissue damage. The two major approaches are adopted by the host to deal with the pathogen: resistance or tolerance. The nature of the responses often differs between species and between individuals of the same species. Resistance includes innate and adaptive immune responses to control virus replication. Disease tolerance relies on the immune response allowing the coexistence of infections in the host with minimal or no clinical signs, while maintaining sufficient viral replication for transmission. Here, we compared the virome of bats, rodents and migratory birds and the molecular mechanisms underlying symptomatic and asymptomatic disease progression. We also explore the influence of the host physiology and environmental influences on RNA virus expression and how it impacts on the whole brain transcriptome of seemingly healthy semipalmated sandpiper (Calidris pusilla) and spotted sandpiper (Actitis macularius). Three time points throughout the year were selected to understand the importance of longitudinal surveys in the characterization of the virome. We finally revisited evidence that upstream and downstream regulation of the inflammatory response is, respectively, associated with resistance and tolerance to viral infections.

Navigating the Landscape: A Comprehensive Review of Current Virus Databases

Viruses are abundant and diverse entities that have important roles in public health, ecology, and agriculture. The identification and surveillance of viruses rely on an understanding of their genome organization, sequences, and replication strategy. Despite technological advancements in sequencing methods, our current understanding of virus diversity remains incomplete, highlighting the need to explore undiscovered viruses. Virus databases play a crucial role in providing access to sequences, annotations and other metadata, and analysis tools for studying viruses. However, there has not been a comprehensive review of virus databases in the last five years. This study aimed to fill this gap by identifying 24 active virus databases and included an extensive evaluation of their content, functionality and compliance with the FAIR principles. In this study, we thoroughly assessed the search capabilities of five database catalogs, which serve as comprehensive repositories housing a diverse array of databases and offering essential metadata. Moreover, we conducted a comprehensive review of different types of errors, encompassing taxonomy, names, missing information, sequences, sequence orientation, and chimeric sequences, with the intention of empowering users to effectively tackle these challenges. We expect this review to aid users in selecting suitable virus databases and other resources, and to help databases in error management and improve their adherence to the FAIR principles. The databases listed here represent the current knowledge of viruses and will help aid users find databases of interest based on content, functionality, and scope. The use of virus databases is integral to gaining new insights into the biology, evolution, and transmission of viruses, and developing new strategies to manage virus outbreaks and preserve global health.

Molecular Changes in the Brain of the Wintering Calidris pusilla in the Mangroves of the Amazon River Estuary

Migrant birds prepare differently to fly north for breeding in the spring and for the flight to lower latitudes during autumn, avoiding the cold and food shortages of the Northern Hemisphere's harsh winter. The molecular events associated with these fundamental stages in the life history of migrants include the differential gene expression in different tissues. Semipalmated sandpipers (Calidris pusilla) are Arctic-breeding shorebirds that migrate to the coast of South America during the non-breeding season. In a previous study, we demonstrated that between the beginning and the end of the wintering period, substantial glial changes and neurogenesis occur in the brain of C. pusilla. These changes follow the epic journey of the autumn migration when a 5-day non-stop transatlantic flight towards the coast of South America and the subsequent preparation for the long-distance flight of the spring migration takes place. Here, we tested the hypothesis that the differential gene expressions observed in the brains of individuals captured in the autumn and spring windows are consistent with the previously described cellular changes. We searched for differential gene expressions in the brain of the semipalmated sandpiper, of recently arrived birds (RA) from the autumnal migration, and that of individuals in the premigratory period (PM) in the spring. All individuals were collected in the tropical coastal of northern Brazil in the mangrove region of the Amazon River estuary. We generated a de novo neurotranscriptome for C. pusilla individuals and compared the gene expressions across libraries. To that end, we mapped an RNA-Seq that reads to the C. pusilla neurotranscriptome in four brain samples of each group and found that the differential gene expressions in newly arrived and premigratory birds were related with neurogenesis, metabolic pathways (ketone body biosynthetic and the catabolic and lipid biosynthetic processes), and glial changes (astrocyte-dopaminergic neuron signaling, astrocyte differentiation, astrocyte cell migration, and astrocyte activation involved in immune response), as well as genes related to the immune response to virus infections (Type I Interferons), inflammatory cytokines (IL-6, IL-1β, TNF, and NF-κB), NLRP3 inflammasome, anti-inflammatory cytokines (IL-10), and cell death pathways (pyroptosis- and caspase-related changes).

VAPEX: an interactive web server for the deep exploration of natural virus and phage genomes

MOTIVATION

Studying the genetic makeup of viruses and phages through genome analysis is crucial for comprehending their function in causing diseases, progressing medicine, tracing their evolutionary history, monitoring the environment, and creating innovative biotechnologies. However, accessing the necessary data can be challenging due to a lack of dedicated comparative genomic tools and viral and phage databases, which are often outdated. Moreover, many wet bench experimentalists may not have the computational proficiency required to manipulate large amounts of genomic data.

RESULTS

We have developed VAPEX (Virus And Phage EXplorer), a web server which is supported by a database and features a user-friendly web interface. This tool enables users to easily perform various genomic analysis queries on all natural viruses and phages that have been fully sequenced and are listed in the NCBI compendium. VAPEX therefore excels in producing visual depictions of fully resolved synteny maps, which is one of its key strengths. VAPEX has the ability to exhibit a vast array of orthologous gene classes simultaneously through the use of symbolic representation. Additionally, VAPEX can fully analyze user-submitted viral and phage genomes, including those that have not yet been annotated.

AVAILABILITY AND IMPLEMENTATION

VAPEX can be accessed from all current web browsers such as Chrome, Firefox, Edge, Safari, and Opera. VAPEX is freely accessible at https://archaea.i2bc.paris-saclay.fr/vapex/.

The dominance of co-circulating SARS-CoV-2 variants in wastewater

SARS-CoV-2 remains one of the biggest global health problems, which has already reached our wastewater through fecal shedding by COVID-19 patients. While the development of vaccines has mitigated the threats of the COVID-19 pandemic, the evolutionary dynamics of SARS-CoV-2 in wastewater lack monitoring and understanding. In this study, SARS-CoV-2 variants in wastewater were identified by analyzing 8511 wastewater-derived genome sequences from 9 countries from March 2020 to May 2023. The dominance of co-circulating variants was observed, namely B.1 in 2020, Alpha and Delta in 2021, then superseded by Omicron lineages in 2022 with a three-times increase. Mutations were also profiled, revealing nearly 5031 unique amino acid substitutions occurring approximately 371,591 times, some of which were associated with enhanced viral transmission and fitness. This study provided the first long-term multi-country overview of the prevalence of co-circulating SARS-CoV-2 lineages and mutations in wastewater and showed its comparison with conventional epidemiological surveillance. The results highlight the ability of wastewater-based genome monitoring to supplement clinical surveillance efforts in rapidly detecting viruses up to the strain level to keep track of their potential transmission routes and evolutionary dynamics.

Using machine learning to detect coronaviruses potentially infectious to humans

Establishing the host range for novel viruses remains a challenge. Here, we address the challenge of identifying non-human animal coronaviruses that may infect humans by creating an artificial neural network model that learns from spike protein sequences of alpha and beta coronaviruses and their binding annotation to their host receptor. The proposed method produces a human-Binding Potential (h-BiP) score that distinguishes, with high accuracy, the binding potential among coronaviruses. Three viruses, previously unknown to bind human receptors, were identified: Bat coronavirus BtCoV/133/2005 and Pipistrellus abramus bat coronavirus HKU5-related (both MERS related viruses), and Rhinolophus affinis coronavirus isolate LYRa3 (a SARS related virus). We further analyze the binding properties of BtCoV/133/2005 and LYRa3 using molecular dynamics. To test whether this model can be used for surveillance of novel coronaviruses, we re-trained the model on a set that excludes SARS-CoV-2 and all viral sequences released after the SARS-CoV-2 was published. The results predict the binding of SARS-CoV-2 with a human receptor, indicating that machine learning methods are an excellent tool for the prediction of host expansion events.

Chromosome folding and prophage activation reveal specific genomic architecture for intestinal bacteria

BACKGROUND

Bacteria and their viruses, bacteriophages, are the most abundant entities of the gut microbiota, a complex community of microorganisms associated with human health and disease. In this ecosystem, the interactions between these two key components are still largely unknown. In particular, the impact of the gut environment on bacteria and their associated prophages is yet to be deciphered.

RESULTS

To gain insight into the activity of lysogenic bacteriophages within the context of their host genomes, we performed proximity ligation-based sequencing (Hi-C) in both in vitro and in vivo conditions on the 12 bacterial strains of the OMM12 synthetic bacterial community stably associated within mice gut (gnotobiotic mouse line OMM12). High-resolution contact maps of the chromosome 3D organization of the bacterial genomes revealed a wide diversity of architectures, differences between environments, and an overall stability over time in the gut of mice. The DNA contacts pointed at 3D signatures of prophages leading to 16 of them being predicted as functional. We also identified circularization signals and observed different 3D patterns between in vitro and in vivo conditions. Concurrent virome analysis showed that 11 of these prophages produced viral particles and that OMM12 mice do not carry other intestinal viruses.

CONCLUSIONS

The precise identification by Hi-C of functional and active prophages within bacterial communities will unlock the study of interactions between bacteriophages and bacteria across conditions (healthy vs disease). Video Abstract.