NCBI viral genomes resource

Reduced bonobo MHC class I diversity predicts a reduced viral peptide binding ability compared to chimpanzees

BACKGROUND

The highly polymorphic genes of the major histocompatibility complex (MHC) class I are involved in defense against viruses and other intracellular pathogens. Although several studies found reduced MHC class I diversity in bonobos in comparison to the closely related chimpanzee, it is unclear if this lower diversity also influences the functional ability of MHC class I molecules in bonobos. Here, we use a bioinformatic approach to analyze the viral peptide binding ability of all published bonobo MHC class I molecules (n = 58) in comparison to all published chimpanzee MHC class I molecules (n = 161) for the class I loci A, B, C and A-like.

RESULTS

We examined the peptide binding ability of all 219 different MHC class I molecules to 5,788,712 peptides derived from 1432 different primate viruses and analyzed the percentage of bound peptides and the overlap of the peptide binding repertoires of the two species. We conducted multiple levels of analysis on the "species"-, "population"- and "individual"-level to account for the characterization of MHC variation in a larger number of chimpanzees and their broader geographic distribution. We found a lower percentage of bound peptides in bonobos at the B locus in the "population"-level comparison and at the B and C loci in the "individual"-level comparison. Furthermore, we found evidence of a limited peptide binding repertoire in bonobos by tree-based visualization of functional clustering of MHC molecules, as well as an analysis of peptides bound by both species.

CONCLUSION

Our results suggest a reduced MHC class I viral peptide binding ability at the B and C loci in bonobos compared to chimpanzees. The effects of this finding on the immune defense against viruses in wild living bonobos are unclear. However, special caution is needed to prevent introduction and spread of new viruses to bonobos, as their defensive ability to cope with new viruses could be limited compared to chimpanzees.

Metagenomic Sequencing of HIV-1 in the Blood and Female Genital Tract Reveals Little Quasispecies Diversity during Acute Infection

Due to error-prone replication, HIV-1 generates a diverse population of viruses within a chronically infected individual. When HIV-1 is transmitted to a new individual, one or a few viruses establish the new infection, leading to a genetic bottleneck in the virus population. Understanding the timing and nature of this bottleneck may provide insight into HIV-1 vaccine design and other preventative strategies. We examined the HIV-1 population in three women enrolled in a unique prospective cohort in South Africa who were followed closely during the earliest stages of HIV-1 infection. We found very little HIV-1 diversity in the blood and female genital tract during the first 2 weeks after virus was detected in the bloodstream. These results are compatible with a very early HIV-1 population bottleneck, suggesting the need to study the HIV-1 population in the female genital tract before virus is detectable in the bloodstream.

Heterosexual transmission of human immunodeficiency virus type 1 (HIV-1) is associated with a significant bottleneck in the viral quasispecies population, yet the timing of that bottleneck is poorly understood. We characterized HIV-1 diversity in the blood and female genital tract (FGT) within 2 weeks after detection of infection in three women enrolled in a unique prospective cohort in South Africa. We assembled full-length HIV-1 genomes from matched cervicovaginal lavage (CVL) samples and plasma. Deep sequencing allowed us to identify intrahost single-nucleotide variants (iSNVs) and to characterize within-sample HIV-1 diversity. Our results demonstrated very little HIV-1 diversity in the FGT and plasma by the time viremia was detectable. Within each subject, the consensus HIV-1 sequences were identical in plasma and CVL fluid. No iSNV was present at >6% frequency. One subject had 77 low-frequency iSNVs across both CVL fluid and plasma, another subject had 14 iSNVs in only CVL fluid from the earliest time point, and the third subject had no iSNVs in CVL fluid or plasma. Overall, the small amount of diversity that we detected was greater in the FGT than in plasma and declined over the first 2 weeks after viremia was detectable, compatible with a very early HIV-1 transmission bottleneck. To our knowledge, our study represents the earliest genomic analysis of HIV-1 in the FGT after transmission. Further, the use of metagenomic sequencing allowed us to characterize other organisms in the FGT, including commensal bacteria and sexually transmitted infections, highlighting the utility of the method to sequence both HIV-1 and its metagenomic environment.

IMPORTANCE Due to error-prone replication, HIV-1 generates a diverse population of viruses within a chronically infected individual. When HIV-1 is transmitted to a new individual, one or a few viruses establish the new infection, leading to a genetic bottleneck in the virus population. Understanding the timing and nature of this bottleneck may provide insight into HIV-1 vaccine design and other preventative strategies. We examined the HIV-1 population in three women enrolled in a unique prospective cohort in South Africa who were followed closely during the earliest stages of HIV-1 infection. We found very little HIV-1 diversity in the blood and female genital tract during the first 2 weeks after virus was detected in the bloodstream. These results are compatible with a very early HIV-1 population bottleneck, suggesting the need to study the HIV-1 population in the female genital tract before virus is detectable in the bloodstream.

Minimum Information about an Uncultivated Virus Genome (MIUViG)

We present an extension of the Minimum Information about any (x) Sequence (MIxS) standard for reporting sequences of uncultivated virus genomes. Minimum Information about an Uncultivated Virus Genome (MIUViG) standards were developed within the Genomic Standards Consortium framework and include virus origin, genome quality, genome annotation, taxonomic classification, biogeographic distribution and in silico host prediction. Community-wide adoption of MIUViG standards, which complement the Minimum Information about a Single Amplified Genome (MISAG) and Metagenome-Assembled Genome (MIMAG) standards for uncultivated bacteria and archaea, will improve the reporting of uncultivated virus genomes in public databases. In turn, this should enable more robust comparative studies and a systematic exploration of the global virosphere.

Capturing sequence diversity in metagenomes with comprehensive and scalable probe design

Metagenomic sequencing has the potential to transform microbial detection and characterization, but new tools are needed to improve its sensitivity. Here we present CATCH, a computational method to enhance nucleic acid capture for enrichment of diverse microbial taxa. CATCH designs optimal probe sets, with a specified number of oligonucleotides, that achieve full coverage of, and scale well with, known sequence diversity. We focus on applying CATCH to capture viral genomes in complex metagenomic samples. We design, synthesize, and validate multiple probe sets, including one that targets the whole genomes of the 356 viral species known to infect humans. Capture with these probe sets enriches unique viral content on average 18-fold, allowing us to assemble genomes that could not be recovered without enrichment, and accurately preserves within-sample diversity. We also use these probe sets to recover genomes from the 2018 Lassa fever outbreak in Nigeria and to improve detection of uncharacterized viral infections in human and mosquito samples. The results demonstrate that CATCH enables more sensitive and cost-effective metagenomic sequencing.

A longitudinal systems immunologic investigation of acute Zika virus infection in an individual infected while traveling to Caracas, Venezuela

Zika virus (ZIKV) is an emerging mosquito-borne flavivirus linked to devastating neurologic diseases. Immune responses to flaviviruses may be pathogenic or protective. Our understanding of human immune responses to ZIKV in vivo remains limited. Therefore, we performed a longitudinal molecular and phenotypic characterization of innate and adaptive immune responses during an acute ZIKV infection. We found that innate immune transcriptional and genomic responses were both cell type- and time-dependent. While interferon stimulated gene induction was common to all innate immune cells, the upregulation of important inflammatory cytokine genes was primarily limited to monocyte subsets. Additionally, genomic analysis revealed substantial chromatin remodeling at sites containing cell-type specific transcription factor binding motifs that may explain the observed changes in gene expression. In this dengue virus-experienced individual, adaptive immune responses were rapidly mobilized with T cell transcriptional activity and ZIKV neutralizing antibody responses peaking 6 days after the onset of symptoms. Collectively this study characterizes the development and resolution of an in vivo human immune response to acute ZIKV infection in an individual with pre-existing flavivirus immunity.

Zika virus (ZIKV) is an emerging flaviviral infection that causes significant clinical disease. It is estimated that approximately one half of the world’s population is at risk for ZIKV infection. There are only a limited number of studies describing the human immune response to ZIKV infection. Carlin et al. combined conventional and genomic approaches to longitudinally analyze the innate and adaptive immune responses to acute ZIKV infection and its resolution in a person who was infected while traveling in Venezuela during the 2016 ZIKV epidemic year. Genome-wide sequencing in individual cell types revealed that although many populations respond to interferon stimulation, only specific cell populations within peripheral blood mononuclear cells upregulate important inflammatory cytokine gene expression. Additionally, analysis of open chromatin using ATAC-seq suggests that chromatin remodeling at sites containing cell-type specific transcription factor binding motifs may help us understand changes in gene expression. Consistent with previous reports, this individual with prior exposure to dengue virus (DENV), rapidly developed neutralizing anti-ZIKV responses that were cross-reactive with multiple DENV serotypes. Collectively this study combines traditional and genomic approaches to characterize the cell-type specific development of an in vivo human immune response to acute ZIKV infection.

Detecting life signatures with RNA sequence similarity measures

The RNA World is currently the most plausible hypothesis for explaining the origins of life on Earth. The supporting body of evidence is growing and it comes from multiple areas, including astrobiology, chemistry, biology, mathematics, and, in particular, from computer simulations. Such methods frequently assume the existence of a hypothetical species on Earth, around three billion years ago, with a base sequence probably dissimilar from any in known genomes. However, it is often hard to verify whether or not a hypothetical sequence has the characteristics of biological sequences, and is thus likely to be functional. The primary objective of the presented research was to verify the possibility of building a computational 'life probe' for determining whether a given genetic sequence is biological, and assessing the sensitivity of such probes to the signatures of life present in known biological sequences. We have proposed decision algorithms based on the normalized compression distance (NCD) and Levenshtein distance (LD). We have validated the proposed method in the context of the RNA World hypothesis using short genetic sequences shorter than the error threshold value (i.e., 100 nucleotides). We have demonstrated that both measures can be successfully used to construct life probes that are significantly better than a random decision procedure, while varying from each other when it comes to detailed characteristics. We also observed that fragments of sequences related to replication have better discriminatory power than sequences having other molecular functions. In a broader context, this shows that the signatures of life in short RNA samples can be effectively detected using relatively simple means.

Quito's virome: Metagenomic analysis of viral diversity in urban streams of Ecuador's capital city

In Quito, the microbiological contamination of surface water represents a public health problem, mainly due to the lack of sewage treatment from urban wastewater. Contaminated water contributes to the transmission of many enteric pathogens through direct consumption, agricultural and recreational use. Among the different pathogens present in urban discharges, viruses play an important role on disease, being causes of gastroenteritis, hepatitis, meningitis, respiratory infections, among others. This study analyzes the presence of viruses in highly impacted surface waters of urban rivers using next-generation sequencing techniques. Three representative locations of urban rivers, receiving the main discharges from Quito sewerage system, were selected. Water samples of 500 mL were concentrated by skimmed-milk flocculation method and the viral nucleic acid was extracted and processed for high throughput sequencing using Illumina MiSeq. The results yielded very relevant data of circulating viruses in the capital of Ecuador. A total of 29 viral families were obtained, of which 26 species were associated with infections in humans. Among the 26 species identified, several were related to gastroenteritis: Human Mastadenovirus F, Bufavirus, Sapporovirus, Norwalk virus and Mamastrovirus 1. Also detected were: Gammapapillomavirus associated with skin infections, Polyomavirus 1 related to cases of kidney damage, Parechovirus A described as cause of neonatal sepsis with neurological affectations and Hepatovirus A, the etiologic agent of Hepatitis A. Other emergent viruses identified, of which its pathogenicity remains to be fully clarified, were: Bocavirus, Circovirus, Aichi Virus and Cosavirus. The wide diversity of species detected through metagenomics gives us key information about the public health risks present in the urban rivers of Quito. In addition, this study describes for the first time the presence of important infectious agents not previously reported in Ecuador and with very little reports in Latin America.

Genomic Signatures Among Acanthamoeba polyphaga Entoorganisms Unveil Evidence of Coevolution

The definition of a genomic signature (GS) is "the total net response to selective pressure". Recent isolation and sequencing of naturally occurring organisms, hereby named entoorganisms, within Acanthamoeba polyphaga, raised the hypothesis of a common genomic signature despite their diverse and unrelated evolutionary origin. Widely accepted and implemented tests for GS detection are oligonucleotide relative frequencies (OnRF) and relative codon usage (RCU) surveys. A common pattern and strong correlations were unveiled from OnRFs among A. polyphaga's Mimivirus and virophage Sputnik. RCU showed a common A-T bias at third codon position. We expanded tests to the amoebal mitochondrial genome and amoeba-resistant bacteria, achieving strikingly coherent results to the aforementioned viral analyses. The GSs in these entoorganisms of diverse evolutionary origin are coevolutionarily conserved within an intracellular environment that provides sanctuary for species of ecological and biomedical relevance.

KrakenUniq: confident and fast metagenomics classification using unique k-mer counts

False-positive identifications are a significant problem in metagenomics classification. We present KrakenUniq, a novel metagenomics classifier that combines the fast k-mer-based classification of Kraken with an efficient algorithm for assessing the coverage of unique k-mers found in each species in a dataset. On various test datasets, KrakenUniq gives better recall and precision than other methods and effectively classifies and distinguishes pathogens with low abundance from false positives in infectious disease samples. By using the probabilistic cardinality estimator HyperLogLog, KrakenUniq runs as fast as Kraken and requires little additional memory. KrakenUniq is freely available at https://github.com/fbreitwieser/krakenuniq .

Discovery of highly divergent lineages of plant-associated astro-like viruses sheds light on the emergence of potyviruses

RNA viruses are believed to have originated from a common ancestor, but how this ancestral genome evolved into the large variety of genomic architectures and viral proteomes we see today remains largely unknown. Tackling this question is hindered by the lack of universally conserved proteins other than the RNA-dependent RNA polymerase (RdRp) as well as a limited RNA virus sampling. The latter is still heavily biased towards relatively few viral lineages from a non-representative collection of hosts, which complicates studies aiming to reveal possible trajectories during the evolution of RNA virus genomes that are favored over others. We report the discovery of 11 highly divergent lineages of viruses with genomic architectures that resemble those of the astroviruses. These genomes were initially identified through a sequence homology search in more than 6600 plant transcriptome projects from the Sequence Read Archive (SRA) using astrovirus representatives as query. Seed-based viral genome assembly of unprocessed SRA data for several dozens of the most promising hits resulted in two viral genome sequences with full-length coding regions, nine partial genomes and a much larger number of short sequence fragments. Genomic and phylogenetic characterization of the 11 discovered viruses, which we coined plastroviruses (plant-associated astro-like viruses), showed that they are related to both astro- and potyviruses and allowed us to identify divergent Serine protease, RdRp and viral capsid domains encoded in the plastrovirus genome. Interestingly, some of the plastroviruses shared different features with potyviruses including the replacement of the catalytic Ser by a Cys residue in the protease active site. These results suggest that plastroviruses may have reached different points on an evolutionary trajectory from astro-like to poty-like genomes. A model how potyviruses might have emerged from (pl)astro-like ancestors in a multi-step process is discussed.

A planarian nidovirus expands the limits of RNA genome size

RNA viruses are the only known RNA-protein (RNP) entities capable of autonomous replication (albeit within a permissive host environment). A 33.5 kilobase (kb) nidovirus has been considered close to the upper size limit for such entities; conversely, the minimal cellular DNA genome is in the 100–300 kb range. This large difference presents a daunting gap for the transition from primordial RNP to contemporary DNA-RNP-based life. Whether or not RNA viruses represent transitional steps towards DNA-based life, studies of larger RNA viruses advance our understanding of the size constraints on RNP entities and the role of genome size in virus adaptation. For example, emergence of the largest previously known RNA genomes (20–34 kb in positive-stranded nidoviruses, including coronaviruses) is associated with the acquisition of a proofreading exoribonuclease (ExoN) encoded in the open reading frame 1b (ORF1b) in a monophyletic subset of nidoviruses. However, apparent constraints on the size of ORF1b, which encodes this and other key replicative enzymes, have been hypothesized to limit further expansion of these viral RNA genomes. Here, we characterize a novel nidovirus (planarian secretory cell nidovirus; PSCNV) whose disproportionately large ORF1b-like region including unannotated domains, and overall 41.1-kb genome, substantially extend the presumed limits on RNA genome size. This genome encodes a predicted 13,556-aa polyprotein in an unconventional single ORF, yet retains canonical nidoviral genome organization and expression, as well as key replicative domains. These domains may include functionally relevant substitutions rarely or never before observed in highly conserved sites of RdRp, NiRAN, ExoN and 3CLpro. Our evolutionary analysis suggests that PSCNV diverged early from multi-ORF nidoviruses, and acquired additional genes, including those typical of large DNA viruses or hosts, e.g. Ankyrin and Fibronectin type II, which might modulate virus-host interactions. PSCNV's greatly expanded genome, proteomic complexity, and unique features–impressive in themselves–attest to the likelihood of still-larger RNA genomes awaiting discovery.

RNA viruses are the only known RNA-protein (RNP) entities capable of autonomous replication. The upper genome size for such entities was assumed to be <35 kb; conversely, the minimal cellular DNA genome is in the 100–300 kilobase (kb) range. This large difference presents a daunting gap for the proposed evolution of contemporary DNA-RNP-based life from primordial RNP entities. Here, we describe a nidovirus from planarians, named planarian secretory cell nidovirus (PSCNV), whose 41.1 kb genome is 23% larger than any riboviral genome yet discovered. This increase is nearly equivalent in size to the entire poliovirus genome, and it equips PSCNV with an unprecedented extra coding capacity to adapt. PSCNV has broken apparent constraints on the size of the genomic subregion that encodes core replication machinery in other nidoviruses, including coronaviruses, and has acquired genes not previously observed in RNA viruses. This virus challenges and advances our understanding of the limits to RNA genome size.

A diversity-generating retroelement encoded by a globally ubiquitous Bacteroides phage

BACKGROUND

Diversity-generating retroelements (DGRs) are genetic cassettes that selectively mutate target genes to produce hypervariable proteins. First characterized in Bordetella bacteriophage BPP-1, the DGR creates a hypervariable phage tail fiber that enables host tropism switching. Subsequent surveys for DGRs conclude that the majority identified to date are bacterial or archaeal in origin. This work examines bacteriophage and bacterial genomes for novel phage-encoded DGRs.

RESULTS

This survey discovered 92 DGRs that were only found in phages exhibiting a temperate lifestyle. The majority of phage-encoded DGRs were identified as prophages in bacterial hosts from the phyla Bacteroidetes, Proteobacteria, and Firmicutes. Sequence reads from these previously unidentified prophages were present in viral metagenomes (viromes), indicating these prophages can produce functional viruses. Five phages possessed hypervariable proteins with structural similarity to the tail fiber of BPP-1, whereas the functions of the remaining DGR target proteins were unknown. A novel temperate phage that harbors a DGR cassette targeting a protein of unknown function was induced from Bacteroides dorei. This phage, here named Bacteroides dorei Hankyphage, lysogenizes 13 different Bacteroides species and was present in 34% and 21% of whole-community metagenomes and human-associated viromes, respectively.

CONCLUSIONS

Here, the number of known DGR-containing phages is increased from four to 92. All of these phages exhibit a temperate lifestyle, including a cosmopolitan human-associated phage. Targeted hypervariation by temperate phages may be a ubiquitous mechanism underlying phage-bacteria interaction in the human microbiome.

Overlapping genes and the proteins they encode differ significantly in their sequence composition from non-overlapping genes

Overlapping genes represent a fascinating evolutionary puzzle, since they encode two functionally unrelated proteins from the same DNA sequence. They originate by a mechanism of overprinting, in which point mutations in an existing frame allow the expression (the "birth") of a completely new protein from a second frame. In viruses, in which overlapping genes are abundant, these new proteins often play a critical role in infection, yet they are frequently overlooked during genome annotation. This results in erroneous interpretation of mutational studies and in a significant waste of resources. Therefore, overlapping genes need to be correctly detected, especially since they are now thought to be abundant also in eukaryotes. Developing better detection methods and conducting systematic evolutionary studies require a large, reliable benchmark dataset of known cases. We thus assembled a high-quality dataset of 80 viral overlapping genes whose expression is experimentally proven. Many of them were not present in databases. We found that overall, overlapping genes differ significantly from non-overlapping genes in their nucleotide and amino acid composition. In particular, the proteins they encode are enriched in high-degeneracy amino acids and depleted in low-degeneracy ones, which may alleviate the evolutionary constraints acting on overlapping genes. Principal component analysis revealed that the vast majority of overlapping genes follow a similar composition bias, despite their heterogeneity in length and function. Six proven mammalian overlapping genes also followed this bias. We propose that this apparently near-universal composition bias may either favour the birth of overlapping genes, or/and result from selection pressure acting on them.

Considerations for Optimization of High-Throughput Sequencing Bioinformatics Pipelines for Virus Detection

High-throughput sequencing (HTS) has demonstrated capabilities for broad virus detection based upon discovery of known and novel viruses in a variety of samples, including clinical, environmental, and biological. An important goal for HTS applications in biologics is to establish parameter settings that can afford adequate sensitivity at an acceptable computational cost (computation time, computer memory, storage, expense or/and efficiency), at critical steps in the bioinformatics pipeline, including initial data quality assessment, trimming/cleaning, and assembly (to reduce data volume and increase likelihood of appropriate sequence identification). Additionally, the quality and reliability of the results depend on the availability of a complete and curated viral database for obtaining accurate results; selection of sequence alignment programs and their configuration, that retains specificity for broad virus detection with reduced false-positive signals; removal of host sequences without loss of endogenous viral sequences of interest; and use of a meaningful reporting format, which can retain critical information of the analysis for presentation of readily interpretable data and actionable results. Furthermore, after alignment, both automated and manual evaluation may be needed to verify the results and help assign a potential risk level to residual, unmapped reads. We hope that the collective considerations discussed in this paper aid toward optimization of data analysis pipelines for virus detection by HTS.

Genome-wide diversity in temporal and regional populations of the betabaculovirus Erinnyis ello granulovirus (ErelGV)

Background

Erinnyis ello granulovirus (ErelGV) is a betabaculovirus infecting caterpillars of the sphingid moth E. ello ello (cassava hornworm), an important pest of cassava crops (Manihot esculenta). In this study, the genome of seven field isolates of the virus ErelGV were deep sequenced and their inter- and intrapopulational sequence diversity were analyzed.

Results

No events of gene gain/loss or translocations were observed, and indels were mainly found within highly repetitive regions (direct repeats, drs). A naturally occurring isolate from Northern Brazil (Acre State, an Amazonian region) has shown to be the most diverse population, with a unique pattern of polymorphisms. Overall, non-synonymous substitutions were found all over the seven genomes, with no specific gathering of mutations on hotspot regions. Independently of their sizes, some ORFs have shown higher levels of non-synonymous changes than others. Non-core genes of known functions and structural genes were among the most diverse ones; and as expected, core genes were the least variable genes. We observed remarkable differences on diversity of paralogous genes, as in multiple copies of p10, fgf, and pep. Another important contrast on sequence diversity was found on genes encoding complex subunits and/or involved in the same biological processes, as late expression factors (lefs) and per os infectivity factors (pifs). Interestingly, several polymorphisms in coding regions lie on sequences encoding specific protein domains.

Conclusions

By comparing and integrating information about inter- and intrapopulational diversity of viral isolates, we provide a detailed description on how evolution operates on field isolates of a betabaculovirus. Our results revealed that 35–41% of the SNPs of ErelGV lead to amino acid changes (non-synonymous substitutions). Some genes, especially non-core genes of unknown functions, tend to accumulate more mutations, while core genes evolve slowly and are more conserved. Additional studies would be necessary to understand the actual effects of such gene variations on viral infection and fitness.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5070-6) contains supplementary material, which is available to authorized users.

High-Resolution Temporal and Spatial Patterns of Virome in Wastewater Treatment Systems

Wastewater treatment plants (WWTPs) are considered reservoirs of viruses, but the diversity and dynamic changes of viruses are not well understood. In this study, we recovered 8478 metagenomic viral contigs (mVCs; >5 kb) from two WWTPs (Shatin, 2806; Shek Wu Hui, 5672) in Hong Kong. Approximately 60% of the mVCs were poorly covered (<35% of genes in identified mVCs) by the current NCBI and IMG/VR viral databases. The temporal profile of the newly identified mVCs among 98 Shatin AS samples collected monthly (for approximately 9 years) revealed the presence of periodic dynamics at an interval of approximately one year (341 days). The spatial distribution pattern of the virome in the wastewater treatment systems showed that shared viral clusters (viral populations categorized based on shared gene content and network analysis) can be globally found among similar samples of wastewater treatment systems, indicating the presence of core viral communities among geographically isolated wastewater treatment systems. These results not only supplemented the current virome database of engineered systems but also, to some extent, expanded the understanding of long-term cyclical development and spatial distributions of viral communities in wastewater treatment systems.

Mimiviridae: An Expanding Family of Highly Diverse Large dsDNA Viruses Infecting a Wide Phylogenetic Range of Aquatic Eukaryotes

Since 1998, when Jim van Etten’s team initiated its characterization, Paramecium bursaria Chlorella virus 1 (PBCV-1) had been the largest known DNA virus, both in terms of particle size and genome complexity. In 2003, the Acanthamoeba-infecting Mimivirus unexpectedly superseded PBCV-1, opening the era of giant viruses, i.e., with virions large enough to be visible by light microscopy and genomes encoding more proteins than many bacteria. During the following 15 years, the isolation of many Mimivirus relatives has made Mimiviridae one of the largest and most diverse families of eukaryotic viruses, most of which have been isolated from aquatic environments. Metagenomic studies of various ecosystems (including soils) suggest that many more remain to be isolated. As Mimiviridae members are found to infect an increasing range of phytoplankton species, their taxonomic position compared to the traditional Phycodnaviridae (i.e., etymologically “algal viruses”) became a source of confusion in the literature. Following a quick historical review of the key discoveries that established the Mimiviridae family, we describe its current taxonomic structure and propose a set of operational criteria to help in the classification of future isolates.

Soil Viruses Are Underexplored Players in Ecosystem Carbon Processing

Rapidly thawing permafrost harbors ∼30 to 50% of global soil carbon, and the fate of this carbon remains unknown. Microorganisms will play a central role in its fate, and their viruses could modulate that impact via induced mortality and metabolic controls. Because of the challenges of recovering viruses from soils, little is known about soil viruses or their role(s) in microbial biogeochemical cycling. Here, we describe 53 viral populations (viral operational taxonomic units [vOTUs]) recovered from seven quantitatively derived (i.e., not multiple-displacement-amplified) viral-particle metagenomes (viromes) along a permafrost thaw gradient at the Stordalen Mire field site in northern Sweden. Only 15% of these vOTUs had genetic similarity to publicly available viruses in the RefSeq database, and ∼30% of the genes could be annotated, supporting the concept of soils as reservoirs of substantial undescribed viral genetic diversity. The vOTUs exhibited distinct ecology, with different distributions along the thaw gradient habitats, and a shift from soil-virus-like assemblages in the dry palsas to aquatic-virus-like assemblages in the inundated fen. Seventeen vOTUs were linked to microbial hosts (in silico), implicating viruses in infecting abundant microbial lineages from Acidobacteria, Verrucomicrobia, and Deltaproteobacteria, including those encoding key biogeochemical functions such as organic matter degradation. Thirty auxiliary metabolic genes (AMGs) were identified and suggested virus-mediated modulation of central carbon metabolism, soil organic matter degradation, polysaccharide binding, and regulation of sporulation. Together, these findings suggest that these soil viruses have distinct ecology, impact host-mediated biogeochemistry, and likely impact ecosystem function in the rapidly changing Arctic. IMPORTANCE This work is part of a 10-year project to examine thawing permafrost peatlands and is the first virome-particle-based approach to characterize viruses in these systems. This method yielded >2-fold-more viral populations (vOTUs) per gigabase of metagenome than vOTUs derived from bulk-soil metagenomes from the same site (J. B. Emerson, S. Roux, J. R. Brum, B. Bolduc, et al., Nat Microbiol 3:870-880, 2018, https://doi.org/10.1038/s41564-018-0190-y). We compared the ecology of the recovered vOTUs along a permafrost thaw gradient and found (i) habitat specificity, (ii) a shift in viral community identity from soil-like to aquatic-like viruses, (iii) infection of dominant microbial hosts, and (iv) carriage of host metabolic genes. These vOTUs can impact ecosystem carbon processing via top-down (inferred from lysing dominant microbial hosts) and bottom-up (inferred from carriage of auxiliary metabolic genes) controls. This work serves as a foundation which future studies can build upon to increase our understanding of the soil virosphere and how viruses affect soil ecosystem services.

Tracing mother-infant transmission of bacteriophages by means of a novel analytical tool for shotgun metagenomic datasets: METAnnotatorX

BACKGROUND

Despite the relevance of viral populations, our knowledge of (bacterio) phage populations, i.e., the phageome, suffers from the absence of a "gold standard" protocol for viral DNA extraction with associated in silico sequence processing analyses. To overcome this apparent hiatus, we present here a comprehensive performance evaluation of various protocols and propose an optimized pipeline that covers DNA extraction, sequencing, and bioinformatic analysis of phageome data.

RESULTS

Five widely used protocols for viral DNA extraction from fecal samples were tested for their performance in removal of non-viral DNA. Moreover, we developed a novel bioinformatic platform, METAnnotatorX, for metagenomic dataset analysis. This in silico tool facilitates a range of read- and assembly-based analyses, including taxonomic profiling using an iterative multi-database pipeline, classification of contigs at genus and species level, as well as functional characterizations of reads and assembled data. Performances of METAnnotatorX were assessed through investigation of seven mother-newborn pairs, leading to the identification of shared phage genotypes, of which two were genomically decoded and characterized. METAnnotatorX was furthermore employed to evaluate a protocol for the identification of contaminant non-viral DNA in sequenced datasets and was exploited to determine the amount of metagenomic data needed for robust evaluation of human adult-derived (fecal) phageomes.

CONCLUSIONS

Results obtained in this study demonstrate that a comprehensive pipeline for analysis of phageomes will be pivotal for future explorations of the ecology of phages in the gut environment as well as for understanding their impact on the physiology and bacterial community kinetics as players of dysbiosis and homeostasis in the gut microbiota.

Predicting MicroRNA Mediated Gene Regulation between Human and Viruses

MicroRNAs (miRNAs) mediate various biological processes by actively fine-tuning gene expression at the post-transcriptional level. With the identification of numerous human and viral miRNAs, growing evidence has indicated a common role of miRNAs in mediating the interactions between humans and viruses. However, there is only limited information about Cross-Kingdom miRNA target sites from studies. To facilitate an extensive investigation on the interplay among the gene regulatory networks of humans and viruses, we designed a prediction pipeline, mirTarP, that is suitable for miRNA target screening on the genome scale. By applying mirTarP, we constructed the database mirTar, which is a comprehensive miRNA target repository of bidirectional interspecies regulation between viruses and humans. To provide convenient downloading for users from both the molecular biology field and medical field, mirTar classifies viruses according to “ICTV viral category” and the “medical microbiology classification” on the web page. The mirTar database and mirTarP tool are freely available online.

WIsH: who is the host? Predicting prokaryotic hosts from metagenomic phage contigs

Summary

WIsH predicts prokaryotic hosts of phages from their genomic sequences. It achieves 63% mean accuracy when predicting the host genus among 20 genera for 3 kbp-long phage contigs. Over the best current tool, WisH shows much improved accuracy on phage sequences of a few kbp length and runs hundreds of times faster, making it suited for metagenomics studies.

Availability and implementation

OpenMP-parallelized GPL-licensed C ++ code available at https://github.com/soedinglab/wish.

Supplementary information

Supplementary data are available at Bioinformatics online.

Multiscale Analysis of Independent Alzheimer's Cohorts Finds Disruption of Molecular, Genetic, and Clinical Networks by Human Herpesvirus

Investigators have long suspected that pathogenic microbes might contribute to the onset and progression of Alzheimer's disease (AD) although definitive evidence has not been presented. Whether such findings represent a causal contribution, or reflect opportunistic passengers of neurodegeneration, is also difficult to resolve. We constructed multiscale networks of the late-onset AD-associated virome, integrating genomic, transcriptomic, proteomic, and histopathological data across four brain regions from human post-mortem tissue. We observed increased human herpesvirus 6A (HHV-6A) and human herpesvirus 7 (HHV-7) from subjects with AD compared with controls. These results were replicated in two additional, independent and geographically dispersed cohorts. We observed regulatory relationships linking viral abundance and modulators of APP metabolism, including induction of APBB2, APPBP2, BIN1, BACE1, CLU, PICALM, and PSEN1 by HHV-6A. This study elucidates networks linking molecular, clinical, and neuropathological features with viral activity and is consistent with viral activity constituting a general feature of AD.

Viromic Analysis of Wastewater Input to a River Catchment Reveals a Diverse Assemblage of RNA Viruses

Detection of viruses in the environment is heavily dependent on PCR-based approaches that require reference sequences for primer design. While this strategy can accurately detect known viruses, it will not find novel genotypes or emerging and invasive viral species. In this study, we investigated the use of viromics, i.e., high-throughput sequencing of the biosphere's viral fraction, to detect human-/animal-pathogenic RNA viruses in the Conwy river catchment area in Wales, United Kingdom. Using a combination of filtering and nuclease treatment, we extracted the viral fraction from wastewater and estuarine river water and sediment, followed by high-throughput RNA sequencing (RNA-Seq) analysis on the Illumina HiSeq platform, for the discovery of RNA virus genomes. We found a higher richness of RNA viruses in wastewater samples than in river water and sediment, and we assembled a complete norovirus genotype GI.2 genome from wastewater effluent, which was not contemporaneously detected by conventional reverse transcription-quantitative PCR (qRT-PCR). The simultaneous presence of diverse rotavirus signatures in wastewater indicated the potential for zoonotic infections in the area and suggested runoff from pig farms as a possible origin of these viruses. Our results show that viromics can be an important tool in the discovery of pathogenic viruses in the environment and can be used to inform and optimize reference-based detection methods provided appropriate and rigorous controls are included. IMPORTANCE Enteric viruses cause gastrointestinal illness and are commonly transmitted through the fecal-oral route. When wastewater is released into river systems, these viruses can contaminate the environment. Our results show that we can use viromics to find the range of potentially pathogenic viruses that are present in the environment and identify prevalent genotypes. The ultimate goal is to trace the fate of these pathogenic viruses from origin to the point where they are a threat to human health, informing reference-based detection methods and water quality management.

Viromic Analysis of Wastewater Input to a River Catchment Reveals a Diverse Assemblage of RNA Viruses

Detection of viruses in the environment is heavily dependent on PCR-based approaches that require reference sequences for primer design. While this strategy can accurately detect known viruses, it will not find novel genotypes or emerging and invasive viral species. In this study, we investigated the use of viromics, i.e., high-throughput sequencing of the biosphere's viral fraction, to detect human-/animal-pathogenic RNA viruses in the Conwy river catchment area in Wales, United Kingdom. Using a combination of filtering and nuclease treatment, we extracted the viral fraction from wastewater and estuarine river water and sediment, followed by high-throughput RNA sequencing (RNA-Seq) analysis on the Illumina HiSeq platform, for the discovery of RNA virus genomes. We found a higher richness of RNA viruses in wastewater samples than in river water and sediment, and we assembled a complete norovirus genotype GI.2 genome from wastewater effluent, which was not contemporaneously detected by conventional reverse transcription-quantitative PCR (qRT-PCR). The simultaneous presence of diverse rotavirus signatures in wastewater indicated the potential for zoonotic infections in the area and suggested runoff from pig farms as a possible origin of these viruses. Our results show that viromics can be an important tool in the discovery of pathogenic viruses in the environment and can be used to inform and optimize reference-based detection methods provided appropriate and rigorous controls are included. IMPORTANCE Enteric viruses cause gastrointestinal illness and are commonly transmitted through the fecal-oral route. When wastewater is released into river systems, these viruses can contaminate the environment. Our results show that we can use viromics to find the range of potentially pathogenic viruses that are present in the environment and identify prevalent genotypes. The ultimate goal is to trace the fate of these pathogenic viruses from origin to the point where they are a threat to human health, informing reference-based detection methods and water quality management.

A comprehensive and quantitative exploration of thousands of viral genomes

The complete assembly of viral genomes from metagenomic datasets (short genomic sequences gathered from environmental samples) has proven to be challenging, so there are significant blind spots when we view viral genomes through the lens of metagenomics. One approach to overcoming this problem is to leverage the thousands of complete viral genomes that are publicly available. Here we describe our efforts to assemble a comprehensive resource that provides a quantitative snapshot of viral genomic trends – such as gene density, noncoding percentage, and abundances of functional gene categories – across thousands of viral genomes. We have also developed a coarse-grained method for visualizing viral genome organization for hundreds of genomes at once, and have explored the extent of the overlap between bacterial and bacteriophage gene pools. Existing viral classification systems were developed prior to the sequencing era, so we present our analysis in a way that allows us to assess the utility of the different classification systems for capturing genomic trends.

Complete genome sequence of peanut virus C, a putative novel ilarvirus

We determined the complete genome sequence of a putative novel ilarvirus, tentatively named "peanut virus C" (PVC), identified in peanut (Arachis hypogaea). The three segmented genomic RNA molecules of PVC were 3474 (RNA1), 2925 (RNA2), and 2160 (RNA3) nucleotides in length, with five predicted open reading frames containing conserved domains and motifs that are typical features of ilarviruses. The three genomic RNAs shared nucleotide sequence similarity (74% identity and 93% query coverage for RNA1, 75% identity and 85% query coverage for RNA2, and 72% identity and 70% query coverage for RNA3) with the most closely related ilarvirus, parietaria mottle virus. These results suggest that PVC is a novel member of the genus Ilarvirus in the family Bromoviridae.

Reproducible protocols for metagenomic analysis of human faecal phageomes

BACKGROUND

Recent studies have demonstrated that the human gut is populated by complex, highly individual and stable communities of viruses, the majority of which are bacteriophages. While disease-specific alterations in the gut phageome have been observed in IBD, AIDS and acute malnutrition, the human gut phageome remains poorly characterised. One important obstacle in metagenomic studies of the human gut phageome is a high level of discrepancy between results obtained by different research groups. This is often due to the use of different protocols for enriching virus-like particles, nucleic acid purification and sequencing. The goal of the present study is to develop a relatively simple, reproducible and cost-efficient protocol for the extraction of viral nucleic acids from human faecal samples, suitable for high-throughput studies. We also analyse the effect of certain potential confounding factors, such as storage conditions, repeated freeze-thaw cycles, and operator bias on the resultant phageome profile. Additionally, spiking of faecal samples with an exogenous phage standard was employed to quantitatively analyse phageomes following metagenomic sequencing. Comparative analysis of phageome profiles to bacteriome profiles was also performed following 16S rRNA amplicon sequencing.

RESULTS

Faecal phageome profiles exhibit an overall greater individual specificity when compared to bacteriome profiles. The phageome and bacteriome both exhibited moderate change when stored at + 4 °C or room temperature. Phageome profiles were less impacted by multiple freeze-thaw cycles than bacteriome profiles, but there was a greater chance for operator effect in phageome processing. The successful spiking of faecal samples with exogenous bacteriophage demonstrated large variations in the total viral load between individual samples.

CONCLUSIONS

The faecal phageome sequencing protocol developed in this study provides a valuable additional view of the human gut microbiota that is complementary to 16S amplicon sequencing and/or metagenomic sequencing of total faecal DNA. The protocol was optimised for several confounding factors that are encountered while processing faecal samples, to reduce discrepancies observed within and between research groups studying the human gut phageome. Rapid storage, limited freeze-thaw cycling and spiking of faecal samples with an exogenous phage standard are recommended for optimum results.

The End of a 60-year Riddle: Identification and Genomic Characterization of an Iridovirus, the Causative Agent of White Fat Cell Disease in Zooplankton

The planktonic freshwater crustacean of the genus Daphnia are a model system for biomedical research and, in particular, invertebrate-parasite interactions. Up until now, no virus has been characterized for this system. Here we report the discovery of an iridovirus as the causative agent of White Fat Cell Disease (WFCD) in Daphnia WFCD is a highly virulent disease of Daphnia that can easily be cultured under laboratory conditions. Although it has been studied from sites across Eurasia for more than 60 years, its causative agent had not been described, nor had an iridovirus been connected to WFCD before now. Here we find that an iridovirus-the Daphnia iridescent virus 1 (DIV-1)-is the causative agent of WFCD. DIV-1 has a genome sequence of about 288 kbp, with 39% G+C content and encodes 367 predicted open reading frames. DIV-1 clusters together with other invertebrate iridoviruses but has by far the largest genome among all sequenced iridoviruses. Comparative genomics reveal that DIV-1 has apparently recently lost a substantial number of unique genes but has also gained genes by horizontal gene transfer from its crustacean host. DIV-1 represents the first invertebrate iridovirus that encodes proteins to purportedly cap RNA, and it contains unique genes for a DnaJ-like protein, a membrane glycoprotein and protein of the immunoglobulin superfamily, which may mediate host-pathogen interactions and pathogenicity. Our findings end a 60-year search for the causative agent of WFCD and add to our knowledge of iridovirus genomics and invertebrate-virus interactions.

Analysis of IAV Replication and Co-infection Dynamics by a Versatile RNA Viral Genome Labeling Method

Genome delivery to the proper cellular compartment for transcription and replication is a primary goal of viruses. However, methods for analyzing viral genome localization and differentiating genomes with high identity are lacking, making it difficult to investigate entry-related processes and co-examine heterogeneous RNA viral populations. Here, we present an RNA labeling approach for single-cell analysis of RNA viral replication and co-infection dynamics in situ, which uses the versatility of padlock probes. We applied this method to identify influenza A virus (IAV) infections in cells and lung tissue with single-nucleotide specificity and to classify entry and replication stages by gene segment localization. Extending the classification strategy to co-infections of IAVs with single-nucleotide variations, we found that the dependence on intracellular trafficking places a time restriction on secondary co-infections necessary for genome reassortment. Altogether, these data demonstrate how RNA viral genome labeling can help dissect entry and co-infections.

Metagenomics for the study of viruses in urban sewage as a tool for public health surveillance

The application of next-generation sequencing (NGS) techniques for the identification of viruses present in urban sewage has not been fully explored. This is partially due to a lack of reliable and sensitive protocols for studying viral diversity and to the highly complex analysis required for NGS data processing. One important step towards this goal is finding methods that can efficiently concentrate viruses from sewage samples. Here the application of a virus concentration method based on skimmed milk organic flocculation (SMF) using 10L of sewage collected in different seasons enabled the detection of many viruses. However, some viruses, such as human adenoviruses, could not always be detected using metagenomics, even when quantitative PCR (qPCR) assessments were positive. A targeted metagenomic assay for adenoviruses was conducted and 59.41% of the obtained reads were assigned to murine adenoviruses. However, up to 20 different human adenoviruses (HAdV) were detected by this targeted assay being the most abundant HAdV-41 (29.24%) and HAdV-51 (1.63%). To improve metagenomics' sensitivity, two different protocols for virus concentration were comparatively analysed: an ultracentrifugation protocol and a lower-volume SMF protocol. The sewage virome contained 41 viral families, including pathogenic viral species from families Caliciviridae, Adenoviridae, Astroviridae, Picornaviridae, Polyomaviridae, Papillomaviridae and Hepeviridae. The contribution of urine to sewage metavirome seems to be restricted to a few specific DNA viral families, including the polyomavirus and papillomavirus species. In experimental infections with sewage in a rhesus macaque model, infective human hepatitis E and JC polyomavirus were identified. Urban raw sewage consists of the excreta of thousands of inhabitants; therefore, it is a representative sample for epidemiological surveillance purposes. The knowledge of the metavirome is of significance to public health, highlighting the presence of viral strains that are circulating within a population while acting as a complex matrix for viral discovery.

A Reference Viral Database (RVDB) To Enhance Bioinformatics Analysis of High-Throughput Sequencing for Novel Virus Detection

To facilitate bioinformatics analysis of high-throughput sequencing (HTS) data for the detection of both known and novel viruses, we have developed a new reference viral database (RVDB) that provides a broad representation of different virus species from eukaryotes by including all viral, virus-like, and virus-related sequences (excluding bacteriophages), regardless of their size. In particular, RVDB contains endogenous nonretroviral elements, endogenous retroviruses, and retrotransposons. Sequences were clustered to reduce redundancy while retaining high viral sequence diversity. A particularly useful feature of RVDB is the reduction of cellular sequences, which can enhance the run efficiency of large transcriptomic and genomic data analysis and increase the specificity of virus detection.

Detection of distantly related viruses by high-throughput sequencing (HTS) is bioinformatically challenging because of the lack of a public database containing all viral sequences, without abundant nonviral sequences, which can extend runtime and obscure viral hits. Our reference viral database (RVDB) includes all viral, virus-related, and virus-like nucleotide sequences (excluding bacterial viruses), regardless of length, and with overall reduced cellular sequences. Semantic selection criteria (SEM-I) were used to select viral sequences from GenBank, resulting in a first-generation viral database (VDB). This database was manually and computationally reviewed, resulting in refined, semantic selection criteria (SEM-R), which were applied to a new download of updated GenBank sequences to create a second-generation VDB. Viral entries in the latter were clustered at 98% by CD-HIT-EST to reduce redundancy while retaining high viral sequence diversity. The viral identity of the clustered representative sequences (creps) was confirmed by BLAST searches in NCBI databases and HMMER searches in PFAM and DFAM databases. The resulting RVDB contained a broad representation of viral families, sequence diversity, and a reduced cellular content; it includes full-length and partial sequences and endogenous nonretroviral elements, endogenous retroviruses, and retrotransposons. Testing of RVDBv10.2, with an in-house HTS transcriptomic data set indicated a significantly faster run for virus detection than interrogating the entirety of the NCBI nonredundant nucleotide database, which contains all viral sequences but also nonviral sequences. RVDB is publically available for facilitating HTS analysis, particularly for novel virus detection. It is meant to be updated on a regular basis to include new viral sequences added to GenBank.

IMPORTANCE To facilitate bioinformatics analysis of high-throughput sequencing (HTS) data for the detection of both known and novel viruses, we have developed a new reference viral database (RVDB) that provides a broad representation of different virus species from eukaryotes by including all viral, virus-like, and virus-related sequences (excluding bacteriophages), regardless of their size. In particular, RVDB contains endogenous nonretroviral elements, endogenous retroviruses, and retrotransposons. Sequences were clustered to reduce redundancy while retaining high viral sequence diversity. A particularly useful feature of RVDB is the reduction of cellular sequences, which can enhance the run efficiency of large transcriptomic and genomic data analysis and increase the specificity of virus detection.

Searching for human oncoviruses: Histories, challenges, and opportunities

Oncoviruses contribute significantly to cancer burden. A century of tumor virological studies have led to the discovery of seven well-accepted human oncoviruses, cumulatively responsible for approximately 15% of human cancer cases. Virus-caused cancers are largely preventable through vaccination. Identifying additional oncoviruses and virus-caused tumors will advance cancer prevention and precision medicine, benefiting affected individuals, and society as a whole. The historic success of finding human oncoviruses has provided a unique lesson for directing new research efforts in the post-sequencing era. Combing the experiences from these pioneer studies with emerging high-throughput techniques will certainly accelerate new discovery and advance our knowledge of the remaining human oncoviruses.

Determinants of translation speed are randomly distributed across transcripts resulting in a universal scaling of protein synthesis times

Ribosome profiling experiments have found greater than 100-fold variation in ribosome density along mRNA transcripts, indicating that individual codon elongation rates can vary to a similar degree. This wide range of elongation times, coupled with differences in codon usage between transcripts, suggests that the average codon translation-rate per gene can vary widely. Yet, ribosome run-off experiments have found that the average codon translation rate for different groups of transcripts in mouse stem cells is constant at 5.6 AA/s. How these seemingly contradictory results can be reconciled is the focus of this study. Here, we combine knowledge of the molecular factors shown to influence translation speed with genomic information from Escherichia coli, Saccharomyces cerevisiae and Homo sapiens to simulate the synthesis of cytosolic proteins in these organisms. The model recapitulates a near constant average translation rate, which we demonstrate arises because the molecular determinants of translation speed are distributed nearly randomly amongst most of the transcripts. Consequently, codon translation rates are also randomly distributed and fast-translating segments of a transcript are likely to be offset by equally probable slow-translating segments, resulting in similar average elongation rates for most transcripts. We also show that the codon usage bias does not significantly affect the near random distribution of codon translation rates because only about 10% of the total transcripts in an organism have high codon usage bias while the rest have little to no bias. Analysis of Ribo-Seq data and an in vivo fluorescent assay supports these conclusions.

ReprDB and panDB: minimalist databases with maximal microbial representation

BACKGROUND

Profiling of shotgun metagenomic samples is hindered by a lack of unified microbial reference genome databases that (i) assemble genomic information from all open access microbial genomes, (ii) have relatively small sizes, and (iii) are compatible to various metagenomic read mapping tools. Moreover, computational tools to rapidly compile and update such databases to accommodate the rapid increase in new reference genomes do not exist. As a result, database-guided analyses often fail to profile a substantial fraction of metagenomic shotgun sequencing reads from complex microbiomes.

RESULTS

We report pipelines that efficiently traverse all open access microbial genomes and assemble non-redundant genomic information. The pipelines result in two species-resolution microbial reference databases of relatively small sizes: reprDB, which assembles microbial representative or reference genomes, and panDB, for which we developed a novel iterative alignment algorithm to identify and assemble non-redundant genomic regions in multiple sequenced strains. With the databases, we managed to assign taxonomic labels and genome positions to the majority of metagenomic reads from human skin and gut microbiomes, demonstrating a significant improvement over a previous database-guided analysis on the same datasets.

CONCLUSIONS

reprDB and panDB leverage the rapid increases in the number of open access microbial genomes to more fully profile metagenomic samples. Additionally, the databases exclude redundant sequence information to avoid inflated storage or memory space and indexing or analyzing time. Finally, the novel iterative alignment algorithm significantly increases efficiency in pan-genome identification and can be useful in comparative genomic analyses.

Know Your Enemy: Successful Bioinformatic Approaches to Predict Functional RNA Structures in Viral RNAs

Structured RNA elements may control virus replication, transcription and translation, and their distinct features are being exploited by novel antiviral strategies. Viral RNA elements continue to be discovered using combinations of experimental and computational analyses. However, the wealth of sequence data, notably from deep viral RNA sequencing, viromes, and metagenomes, necessitates computational approaches being used as an essential discovery tool. In this review, we describe practical approaches being used to discover functional RNA elements in viral genomes. In addition to success stories in new and emerging viruses, these approaches have revealed some surprising new features of well-studied viruses e.g., human immunodeficiency virus, hepatitis C virus, influenza, and dengue viruses. Some notable discoveries were facilitated by new comparative analyses of diverse viral genome alignments. Importantly, comparative approaches for finding RNA elements embedded in coding and non-coding regions differ. With the exponential growth of computer power we have progressed from stem-loop prediction on single sequences to cutting edge 3D prediction, and from command line to user friendly web interfaces. Despite these advances, many powerful, user friendly prediction tools and resources are underutilized by the virology community.

Bacteriophage Taxonomy: An Evolving Discipline

While taxonomy is an often-unappreciated branch of science it serves very important roles. Bacteriophage taxonomy has evolved from a mainly morphology-based discipline, characterized by the work of David Bradley and Hans-Wolfgang Ackermann, to the holistic approach that is taken today. The Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) takes a comprehensive approach to classifying prokaryote viruses measuring overall DNA and protein identity and phylogeny before making decisions about the taxonomic position of a new virus. The huge number of complete genomes being deposited with 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.

The unequivocal preponderance of biocomputation in clinical virology

Bioinformatics and computer based data simulation and modeling are captivating biological research, delivering great results already and promising to deliver more. As biological research is a complex, intricate, diverse field, any available support is gladly taken. With recent outbreaks and epidemics, pathogens are a constant threat to the global economy and security. Virus related plagues are somehow the most difficult to handle. Biocomputation has provided appreciable help in resolving clinical virology related issues. This review, for the first time, surveys the current status of the role of computation in virus related research. Advances made in the fields of clinical virology, antiviral drug design, viral immunology and viral oncology, through input from biocomputation, have been discussed. The amount of progress made and the software platforms available are consolidated in this review. The limitations of computation based methods are presented. Finally, the challenges facing the future of biocomputation in clinical virology are speculated upon.

Biocomputation in clinical virology.

A metagenomic assessment of viral contamination on fresh parsley plants irrigated with fecally tainted river water

Microbial food-borne diseases are still frequently reported despite the implementation of microbial quality legislation to improve food safety. Among all the microbial agents, viruses are the most important causative agents of food-borne outbreaks. The development and application of a new generation of sequencing techniques to test for viral contaminants in fresh produce is an unexplored field that allows for the study of the viral populations that might be transmitted by the fecal-oral route through the consumption of contaminated food. To advance this promising field, parsley was planted and grown under controlled conditions and irrigated using contaminated river water. Viruses polluting the irrigation water and the parsley leaves were studied by using metagenomics. To address possible contamination due to sample manipulation, library preparation, and other sources, parsley plants irrigated with nutritive solution were used as a negative control. In parallel, viruses present in the river water used for plant irrigation were analyzed using the same methodology. It was possible to assign viral taxons from 2.4 to 74.88% of the total reads sequenced depending on the sample. Most of the viral reads detected in the river water were related to the plant viral families Tymoviridae (66.13%) and Virgaviridae (14.45%) and the phage viral families Myoviridae (5.70%), Siphoviridae (5.06%), and Microviridae (2.89%). Less than 1% of the viral reads were related to viral families that infect humans, including members of the Adenoviridae, Reoviridae, Picornaviridae and Astroviridae families. On the surface of the parsley plants, most of the viral reads that were detected were assigned to the Dicistroviridae family (41.52%). Sequences related to important viral pathogens, such as the hepatitis E virus, several picornaviruses from species A and B as well as human sapoviruses and GIV noroviruses were detected. The high diversity of viral sequences found in the parsley plants suggests that irrigation on fecally-tainted food may have a role in the transmission of a wide diversity of viral families. This finding reinforces the idea that the best way to avoid food-borne viral diseases is to introduce good field irrigation and production practices. New strains have been identified that are related to the Picornaviridae and distantly related to the Hepeviridae family. However, the detection of a viral genome alone does not necessarily indicate there is a risk of infection or disease development. Thus, further investigation is crucial for correlating the detection of viral metagenomes in samples with the risk of infection. There is also an urgent need to develop new methods to improve the sensitivity of current Next Generation Sequencing (NGS) techniques in the food safety area.

Influence of Scale on Biomass Growth and Nutrient Removal in an Algal-Bacterial Leachate Treatment System

Data collected from experiments conducted at a flask scale are regularly used as input data for life cycle assessments and techno-economic analyses for predicting the potential productivities of large-scale commercial facilities. This study measures and compares nitrogen removal and biomass growth rates in treatment systems that utilize an algae-bacteria consortium to remediate landfill leachate at three scales: small (0.25 L), medium (100 L), and large (1000 L). The medium- and large-scale vessels were run for 52 consecutive weeks as semibatch reactors under variable environmental conditions. The small-scale experiments were conducted in flasks as batch experiments under controlled environmental conditions. Kolomogov-Smirnov statistical tests, which compare the distributions of entire data sets, were used to determine if the ammonia removal, total nitrogen removal, and biomass growth rates at each scale were statistically different. Results from the Kolmogov-Smirnov comparison indicate that there is a significant difference between all rates determined in the large-scale vessels compared to those in the small-scale vessels. These results suggest that small-scale experiments may not be appropriate as input data in predictive analyses of full scale algal processes. The accumulation of nitrite and nitrate within the reactor, observed midway through the experimental process, is attributed to high relative abundances of ammonia- and nitrite-oxidizing bacteria, identified via metagenomic analysis.

Coral-associated viral communities show high levels of diversity and host auxiliary functions

Stony corals (Scleractinia) are marine invertebrates that form the foundation and framework upon which tropical reefs are built. The coral animal associates with a diverse microbiome comprised of dinoflagellate algae and other protists, bacteria, archaea, fungi and viruses. Using a metagenomics approach, we analysed the DNA and RNA viral assemblages of seven coral species from the central Great Barrier Reef (GBR), demonstrating that tailed bacteriophages of the Caudovirales dominate across all species examined, and ssDNA viruses, notably the Microviridae, are also prevalent. Most sequences with matches to eukaryotic viruses were assigned to six viral families, including four Nucleocytoplasmic Large DNA Viruses (NCLDVs) families: Iridoviridae, Phycodnaviridae, Mimiviridae, and Poxviridae, as well as Retroviridae and Polydnaviridae. Contrary to previous findings, Herpesvirales were rare in these GBR corals. Sequences of a ssRNA virus with similarities to the dinornavirus, Heterocapsa circularisquama ssRNA virus of the Alvernaviridae that infects free-living dinoflagellates, were observed in three coral species. We also detected viruses previously undescribed from the coral holobiont, including a virus that targets fungi associated with the coral species Acropora tenuis. Functional analysis of the assembled contigs indicated a high prevalence of latency-associated genes in the coral-associated viral assemblages, several host-derived auxiliary metabolic genes (AMGs) for photosynthesis (psbA, psbD genes encoding the photosystem II D1 and D2 proteins respectively), as well as potential nematocyst toxins and antioxidants (genes encoding green fluorescent-like chromoprotein). This study expands the currently limited knowledge on coral-associated viruses by characterising viral composition and function across seven GBR coral species.

Novel Microbial Diversity and Functional Potential in the Marine Mammal Oral Microbiome

The vast majority of bacterial diversity lies within phylum-level lineages called "candidate phyla," which lack isolated representatives and are poorly understood. These bacteria are surprisingly abundant in the oral cavity of marine mammals. We employed a genome-resolved metagenomic approach to recover and characterize genomes and functional potential from microbes in the oral gingival sulcus of two bottlenose dolphins (Tursiops truncatus). We detected organisms from 24 known bacterial phyla and one archaeal phylum. We also recovered genomes from two deep-branching, previously uncharacterized phylum-level lineages (here named "Candidatus Delphibacteria" and "Candidatus Fertabacteria"). The Delphibacteria lineage is found in both managed and wild dolphins; its metabolic profile suggests a capacity for denitrification and a possible role in dolphin health. We uncovered a rich diversity of predicted Cas9 proteins, including the two longest predicted Cas9 proteins to date. Notably, we identified the first type II CRISPR-Cas systems encoded by members of the Candidate Phyla Radiation. Using their spacer sequences, we subsequently identified and assembled a complete Saccharibacteria phage genome. These findings underscore the immense microbial diversity and functional potential that await discovery in previously unexplored environments.

Proofreading-Deficient Coronaviruses Adapt for Increased Fitness over Long-Term Passage without Reversion of Exoribonuclease-Inactivating Mutations

The coronavirus (CoV) RNA genome is the largest among the single-stranded positive-sense RNA viruses. CoVs encode a proofreading 3'-to-5' exoribonuclease within nonstructural protein 14 (nsp14-ExoN) that is responsible for CoV high-fidelity replication. Alanine substitution of ExoN catalytic residues [ExoN(-)] in severe acute respiratory syndrome-associated coronavirus (SARS-CoV) and murine hepatitis virus (MHV) disrupts ExoN activity, yielding viable mutant viruses with defective replication, up to 20-fold-decreased fidelity, and increased susceptibility to nucleoside analogues. To test the stability of the ExoN(-) genotype and phenotype, we passaged MHV-ExoN(-) 250 times in cultured cells (P250), in parallel with wild-type MHV (WT-MHV). Compared to MHV-ExoN(-) P3, MHV-ExoN(-) P250 demonstrated enhanced replication and increased competitive fitness without reversion at the ExoN(-) active site. Furthermore, MHV-ExoN(-) P250 was less susceptible than MHV-ExoN(-) P3 to multiple nucleoside analogues, suggesting that MHV-ExoN(-) was under selection for increased replication fidelity. We subsequently identified novel amino acid changes within the RNA-dependent RNA polymerase and nsp14 of MHV-ExoN(-) P250 that partially accounted for the reduced susceptibility to nucleoside analogues. Our results suggest that increased replication fidelity is selected in ExoN(-) CoVs and that there may be a significant barrier to ExoN(-) reversion. These results also support the hypothesis that high-fidelity replication is linked to CoV fitness and indicate that multiple replicase proteins could compensate for ExoN functions during replication.IMPORTANCE Uniquely among RNA viruses, CoVs encode a proofreading exoribonuclease (ExoN) in nsp14 that mediates high-fidelity RNA genome replication. Proofreading-deficient CoVs with disrupted ExoN activity [ExoN(-)] either are nonviable or have significant defects in replication, RNA synthesis, fidelity, fitness, and virulence. In this study, we showed that ExoN(-) murine hepatitis virus can adapt during long-term passage for increased replication and fitness without reverting the ExoN-inactivating mutations. Passage-adapted ExoN(-) mutants also demonstrate increasing resistance to nucleoside analogues that is explained only partially by secondary mutations in nsp12 and nsp14. These data suggest that enhanced resistance to nucleoside analogues is mediated by the interplay of multiple replicase proteins and support the proposed link between CoV fidelity and fitness.

Do Viruses Exchange Genes across Superkingdoms of Life?

Viruses can be classified into archaeoviruses, bacterioviruses, and eukaryoviruses according to the taxonomy of the infected host. The host-constrained perception of viruses implies preference of genetic exchange between viruses and cellular organisms of their host superkingdoms and viral origins from host cells either via escape or reduction. However, viruses frequently establish non-lytic interactions with organisms and endogenize into the genomes of bacterial endosymbionts that reside in eukaryotic cells. Such interactions create opportunities for genetic exchange between viruses and organisms of non-host superkingdoms. Here, we take an atypical approach to revisit virus-cell interactions by first identifying protein fold structures in the proteomes of archaeoviruses, bacterioviruses, and eukaryoviruses and second by tracing their spread in the proteomes of superkingdoms Archaea, Bacteria, and Eukarya. The exercise quantified protein structural homologies between viruses and organisms of their host and non-host superkingdoms and revealed likely candidates for virus-to-cell and cell-to-virus gene transfers. Unexpected lifestyle-driven genetic affiliations between bacterioviruses and Eukarya and eukaryoviruses and Bacteria were also predicted in addition to a large cohort of protein folds that were universally shared by viral and cellular proteomes and virus-specific protein folds not detected in cellular proteomes. These protein folds provide unique insights into viral origins and evolution that are generally difficult to recover with traditional sequence alignment-dependent evolutionary analyses owing to the fast mutation rates of viral gene sequences.

Identification of sapovirus GV.2, astrovirus VA3 and novel anelloviruses in serum from patients with acute hepatitis of unknown aetiology

Hepatitis is a general term meaning inflammation of the liver, which can be caused by a variety of viruses. However, a substantial number of cases remain with unknown aetiology. We analysed the serum of patients with clinical signs of hepatitis using a metagenomics approach to characterize their viral species composition. Four pools of patients with hepatitis without identified aetiological agents were evaluated. Additionally, one pool of patients with hepatitis E (HEV) and pools of healthy volunteers were included as controls. A high diversity of anelloviruses, including novel sequences, was found in pools from patients with hepatitis of unknown aetiology. Moreover, viruses recently associated with gastroenteritis as sapovirus GV.2 and astrovirus VA3 were also detected only in those pools. Besides, most of the HEV genome was recovered from the HEV pool. Finally, GB virus C and human endogenous retrovirus were found in the HEV and healthy pools. Our study provides an overview of the virome in serum from hepatitis patients suggesting a potential role of these viruses not previously described in cases of hepatitis. However, further epidemiologic studies are necessary to confirm their contribution to the development of hepatitis.

Protein-Protein Interactions in Virus-Host Systems

To study virus–host protein interactions, knowledge about viral and host protein architectures and repertoires, their particular evolutionary mechanisms, and information on relevant sources of biological data is essential. The purpose of this review article is to provide a thorough overview about these aspects. Protein domains are basic units defining protein interactions, and the uniqueness of viral domain repertoires, their mode of evolution, and their roles during viral infection make viruses interesting models of study. Mutations at protein interfaces can reduce or increase their binding affinities by changing protein electrostatics and structural properties. During the course of a viral infection, both pathogen and cellular proteins are constantly competing for binding partners. Endogenous interfaces mediating intraspecific interactions—viral–viral or host–host interactions—are constantly targeted and inhibited by exogenous interfaces mediating viral–host interactions. From a biomedical perspective, blocking such interactions is the main mechanism underlying antiviral therapies. Some proteins are able to bind multiple partners, and their modes of interaction define how fast these “hub proteins” evolve. “Party hubs” have multiple interfaces; they establish simultaneous/stable (domain–domain) interactions, and tend to evolve slowly. On the other hand, “date hubs” have few interfaces; they establish transient/weak (domain–motif) interactions by means of short linear peptides (15 or fewer residues), and can evolve faster. Viral infections are mediated by several protein–protein interactions (PPIs), which can be represented as networks (protein interaction networks, PINs), with proteins being depicted as nodes, and their interactions as edges. It has been suggested that viral proteins tend to establish interactions with more central and highly connected host proteins. In an evolutionary arms race, viral and host proteins are constantly changing their interface residues, either to evade or to optimize their binding capabilities. Apart from gaining and losing interactions via rewiring mechanisms, virus–host PINs also evolve via gene duplication (paralogy); conservation (orthology); horizontal gene transfer (HGT) (xenology); and molecular mimicry (convergence). The last sections of this review focus on PPI experimental approaches and their limitations, and provide an overview of sources of biomolecular data for studying virus–host protein interactions.