key: cord-0923644-o912ihsl authors: Lu, Congyu; Peng, Yousong title: Computational Viromics: Applications of the Computational Biology in Viromics Studies date: 2021-05-31 journal: Virol Sin DOI: 10.1007/s12250-021-00395-7 sha: 2d637de0ca2ae9a2af81738ad47d7664fa477521 doc_id: 923644 cord_uid: o912ihsl nan Viruses are a kind of biological entities which rely on host cells for survival. Depending on the genetic materials and replication mode, they can be grouped into double-stranded DNA (dsDNA), single-stranded DNA (ssDNA), doublestranded RNA (dsRNA), positive-sense single-stranded RNA (?ssRNA), negative-sense single-stranded RNA (-ssRNA), ssRNA reverse transcriptase viruses (ssRNA-RT) and dsDNA reverse transcriptase viruses (dsDNA-RT) (Walker et al. 2020) . Viruses can infect most kinds of biological entities, including viruses, bacteria, archaea and eukaryote (La Scola et al. 2008; Fermin, 2018) . They have a great impact on the earth by shaping bacterial population dynamics and balancing the global ecosystem (Suttle, 2007) . For humans, viruses, on the one hand, can cause high human morbidity and mortality and serious economic loss (Baud et al. 2020 ), on the other hand, they can promote and maintain the healthy balance of the gut microbiome (Seo and Kweon, 2019) . Besides, some phages can be applied as the therapy of bacterial infections, especially for the bacterial strains resistant to multiple antibiotics (Altamirano and Barr, 2019) . The viromics studies based on the high-throughput sequencing technology have become increasingly popular in recent years, and novel viruses are being discovered at an unprecedented pace (Gregory et al. 2019) . For example, the Tara Oceans Project recently identified 195,728 viral populations which were more than 10 times as many as the known global ocean DNA virome (Gregory et al. 2019) . However, several challenges exist in analyzing the sequencing data from viromics studies. Firstly, it is difficult to identify all viral nucleotide sequences from the nucleotide sequences that mixed with the sequences of other species and the possible pollutions (Roux et al. 2015a; Ren et al. 2017; Fang et al. 2019; Kieft et al. 2020) ; secondly, the annotation of viral nucleotide sequences is still challenging, especially for those with remote or no homology with the known viruses (Roux et al. 2015b; McNair et al. 2019; Zhang et al. 2019a) ; thirdly, the taxonomic assignment of novel viruses is difficult due to a lack of a unified classification system for viruses (Low et al. 2019) ; fourthly, rapid functional characterization of a large number of newly discovered viruses such as identifying the viral hosts is extremely difficult to achieve by using traditional experimental methods (Jofre and Muniesa, 2020) . According to the above analysis, an emerging area of computational viromics which is defined as using the computational methods to solve the problems in viromics studies was proposed in the present study. It includes but not limited to the following aspects: The first step of viromics studies is to identify the viral nucleotide sequences from the metagenomic sequencing data which often contain lots of DNA sequences of host cells such as bacteria and human (Edwards and Rohwer, 2005) . Both the eukaryotic and prokaryotic viruses are usually identified using approaches based on the homology of marker genes or the genomic sequences with known viral nucleotide sequences, such as ViromeScan (Rampelli et al. 2016) . Unfortunately, viruses lack the universal marker genes like the 16s ribosomal RNA (rRNA) in other species. Therefore, the marker genes used for viral sequence identification should be carefully selected to ensure the coverage of viruses. For example, the RNAdependent RNA polymerases (RdRp) can be taken as the marker gene for the RNA viruses (Wolf et al. 2018) . Nevertheless, the homology-based methods show a significant limitation when they are used in identifying the viruses with large diversification from the known viruses (Gregory et al. 2019) . To overcome the limitation, the sequence homology independent methods have been developed Fang et al. 2019) . For example, the VirFinder identified viral sequences based on the k-mer frequencies in the nucleotide sequences ). In addition, the methods such as VIBRANT combining homology-based methods and homology-independent methods were also further developed to improve the viral sequence identification (Kieft et al. 2020 ). The annotation of viral genomes is important for the further characterization of viruses, and the identification of the genes in viral genomes is most crucial to genome annotation. Currently, there are a large number of methods that can be used for gene prediction (Hyatt et al. 2010; McNair et al. 2019; Zhang et al. 2019a) . Although most of them are not designed for viruses, they should be suitable for the gene predictions of viruses since the viruses can use the machinery of their host cells for transcription and translation (Stern-Ginossar et al. 2019). The characteristics of viral genes including the compact gene structure, no or few introns, overlapped or co-transcribed genes, can be incorporated to optimize the gene prediction methods, such as the PHANOTATE and Vgas (McNair et al. 2019; Zhang et al. 2019a) . Besides, the combination of several different tools may improve gene predictions (Zhang et al. 2019a ). Moreover, the emerging metaproteomics which is defined as large-scale identification and quantification of proteins from microbial communities may help the identification of genes in viromics studies to a great extent (Brum et al. 2016 ). Most of the newly-discovered viruses lack biological features and cannot be classified by the current classification system proposed by the International Committee on Taxonomy of Viruses (ICTV) (Walker et al. 2020) . The usage of viral sequences in virus taxonomy is valid and is supported by the ICTV (Simmonds et al. 2017) . Currently, the homology-based methods have been used for taxonomic assignment of the virome. However, they can only be suitable for a small proportion of viruses with sequence homology to those with known taxonomy (Gregory et al. 2020) . A comprehensive classification of the whole viral sequence space remains challenging due to the lack of universal marker genes in viruses. The gene content-based methods such as vConTACT and GRAViTy proposed in recent studies can accurately classify viruses and may provide novel frameworks for a unified classification of all viruses (Eloe-Fadrosh, 2019). In the era of viromics, molecular evolution on the virome level can provide a global view on the origin, diversity and evolution of viruses (Wolf et al. 2018; Gregory et al. 2019) . Due to the lack of universal marker genes and the large diversity of viral nucleotide sequences, evolutionary analysis was often conducted on a group of viruses which share one or more marker genes (Wolf et al. 2018; Low et al. 2019) . For example, Wolf et al. analyzed the origins and evolution of the RNA virome using the RNA-dependent RNA polymerase (RdRp) (Wolf et al. 2018) . The results obtained from the study determined the evolutionary relationship among the double-stranded RNA, positivestranded RNA and negative-stranded RNA viruses, and revealed the extensive gene module exchange among diverse viruses and the horizontal virus transfer between the distantly related hosts (Wolf et al. 2018) . The identification of viral hosts is essential for characterizing viruses, as viruses must rely on host cells for survival. Host predictions of eukaryotic viruses have been usually conducted based on viral sequences alone, such as those for influenza viruses and coronaviruses (Xu et al. 2017; Tian, 2020) ; while those of prokaryotic viruses have been usually conducted based on the similarity of sequence features or sequences between viruses and hosts. At present, two kinds of computational methods have been developed to predict prokaryotic virus hosts based on genomic sequences (Edwards et al. 2016; Ahlgren et al. 2017; Galiez et al. 2017; Lu et al. 2021 ). The first kind of methods rely on the sequence similarity search between the query viruses and the candidate host genomes since viruses and their hosts may share the same genes and/or short nucleotide sequences such as the spacer sequences used in CRISPR systems (Edwards et al. 2016) . This kind of method can usually predict viral hosts with high accuracy, especially the CRISPR-spacer-based method (Edwards et al. 2016 ). However, they can be only used for a small proportion of viruses since only some viruses have sequence similarities with their hosts (Edwards et al. 2016) . Another kind of methods can predict the viral hosts based on the sequence composition similarity between viruses and their hosts, such as the Prokaryotic virus Host Predictor (PHP) (Lu et al. 2021) , VirHostMatcher and WIsH (Galiez et al. 2017) . Although the latter kind of method predicts viral hosts with lower accuracy than the former, they can be used for any prokaryotic viruses. Viruses can change their hosts or spill over into other species after genetic mutations or recombinations (Letko et al. 2020) , which poses a great challenge for the prediction of viral hosts. For better prediction of viral hosts, it is important to understand the host specificity of viruses determined by the interactions between viruses and their hosts. The complex interactions between virome and hosts are difficult to be resolved in viromics studies (Seo and Kweon, 2019) . Besides, novel viruses are being discovered at an unprecedented pace (Gregory et al. 2019) . So, the high-throughput experimental techniques and computational methods are in urgent need to analyze the interactions between viruses and their hosts (Lasso et al. 2019 ; Lian et al. 2021; Zhang et al. 2020a Zhang et al. , 2020b . The proteinprotein interaction (PPI) prediction methods can help the identification of the interactions between virome and hosts to a great extent (Lasso et al. 2019; Lian et al. 2021) . For example, Lasso et al. used the structural information to develop a computational framework to predict the PPIs between 1,001 human-infecting viruses and human, and they obtained a series of new findings about human-virus interactions such as the shared and unique machinery employed across human-infecting viruses and the previously unappreciated cellular circuits that act on humaninfecting viruses (Lasso et al. 2019 ). The virus isolation and cultivation are the basis for further studies of the virus. Culturomics is defined as a systematic method to find the optimum culture conditions such as the culture medium and the incubation temperature for microbial cultivation (Greub, 2012) . Many achievements in the field of bacteria culturomics have been obtained (Lagier et al. 2018) . For example, Oberhardt et al. integrated known medium databases and a novel prediction tool into a platform that predicts the culture medium given an organism's 16S rRNA sequence (Oberhardt et al. 2015) . Moreover, this platform can also predict culture media for new organisms using a transitivity property and a phylogeny-based collaborative filtering method. Similar to the work by Oberhardt et al., it is possible to predict the cell line or tissue that can be used in virus cultivation based on the similarity of genomic sequences and the predicted PPIs between viruses and hosts. The virome has a significant impact on human health. Previous studies have shown that the virome is associated with multiple diseases. However, the detailed mechanism is still unknown due to the complex interactions between the virome and their hosts (Clooney et al. 2019) . Computational methods are needed to identify the viruses and their roles in causing human diseases. For example, Zhu et al. developed a metagenomic data analysis pipeline, Micro-Pro, to analyze the association between the microbes in the human body and complex diseases . The virome is also closely related to the early warnings of newly emerging viruses. The Global Virome Project (GVP) has estimated that there are 631,000-827,000 unknown viruses with the potential of infecting humans (Carroll et al. 2018) . Recent studies have developed machinelearning methods to identify the human-infecting virome based on sequence features (Zhang et al. 2019b ). More efforts are needed to validate their usage in applications. Taken together, this perspective provides an overall view of computational viromics which includes the identification, annotation and taxonomic assignment of viral genomics sequences, phenotype prediction of viruses, evolution of viromes, virus-host interactions, virus culturomics, association of the virome and human health, and so on ( Table 1) . The computational viromics is still in the beginning stage. Much more computational methods and experimental efforts are needed to characterize the virome and its interactions with the hosts and environments considering the huge diversity of the global virome. Alignmentfree oligonucleotide frequency dissimilarity measure improves prediction of hosts from metagenomically-derived viral sequences Phage therapy in the postantibiotic era Real estimates of mortality following Covid-19 infection Illuminating structural proteins in viral ''dark matter'' with metaproteomics The global virome project Whole-virome analysis sheds light on viral dark matter in inflammatory bowel disease Viral metagenomics Computational approaches to predict bacteriophage-host relationships Towards a genome-based virus taxonomy Ppr-meta: A tool for identifying phages and plasmids from metagenomic fragments using deep learning Host range, host-virus interactions, and virus transmission WIsH: Who is the host? Predicting prokaryotic hosts from metagenomic phage contigs Marine DNA viral macro-and microdiversity from pole to pole The gut virome database reveals age-dependent patterns of virome diversity in the human gut Culturomics: a new approach to study the human microbiome Prodigal: Prokaryotic gene recognition and translation initiation site identification Bacteriophage isolation and characterization: phages of escherichia coli. In horizontal gene transfer Vibrant: Automated recovery, annotation and curation of microbial viruses, and evaluation of viral community function from genomic sequences The virophage as a unique parasite of the giant mimivirus Culturing the human microbiota and culturomics A structure-informed atlas of human-virus interactions Bat-borne virus diversity, spillover and emergence Current status and future perspectives of computational studies on human-virus proteinprotein interactions Evaluation of a concatenated protein phylogeny for classification of tailed double-stranded DNA viruses belonging to the order caudovirales Prokaryotic virus host predictor: a gaussian model for host prediction of prokaryotic viruses in metagenomics Phanotate: A novel approach to gene identification in phage genomes Harnessing the landscape of microbial culture media to predict new organism-media pairings Viromescan: a new tool for metagenomic viral community profiling Virfinder: A novel k-mer based tool for identifying viral sequences from assembled metagenomic data Virsorter: mining viral signal from microbial genomic data Viral dark matter and virus-host interactions resolved from publicly available microbial genomes Virome-host interactions in intestinal health and disease Consensus statement: Virus taxonomy in the age of metagenomics Translational control in virus-infected cells Marine viruses-major players in the global ecosystem The potential intermediate hosts for SARS-CoV-2 Changes to virus taxonomy and the statutes ratified by the international committee on taxonomy of viruses Origins and evolution of the global rna virome Predicting the host of influenza viruses based on the word vector Vgas: a viral genome annotation system Rapid identification of human-infecting viruses Phage protein receptors have multiple interaction partners and high expressions Prediction of the receptorome for the human-infecting virome Micropro: Using metagenomic unmapped reads to provide insights into human microbiota and disease associations Computational Biology in Viromics Studies Acknowledgements This work was supported by the National Key Plan for Scientific Research and Development of China (2016YFD0500300) and Hunan Provincial Natural Science Foundation of China (2020JJ3006).