{'fraction_non_alphanumeric': 0.03203253391270438, 'fraction_numerical': 0.015403470655413498, 'mean_word_length': 5.169788623474824, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Elimination of Foreign Sequences in Eukaryotic Viral Reference Genomes Improves the Accuracy of Virome Analysis

ABSTRACT Widespread in public databases, foreign contaminant sequences pose a substantial obstacle in genomic analyses. Such contamination in viral genome databases is also notorious but more complicated and often causes questionable results in various applications, particularly in virome-based virus detection. Here, we conducted comprehensive screening and identification of the foreign sequences hidden in the largest eukaryotic viral genome collections of GenBank and UniProt using a scrutiny pipeline, which enables us to rigorously detect those problematic viral sequences (PVSs) with origins in hosts, vectors, and laboratory components. As a result, a total of 766 nucleotide PVSs and 276 amino acid PVSs with lengths up to 6,605 bp were determined, which were widely distributed in 39 families with many involving highly public health-concerning viruses, such as hepatitis C virus, Crimean-Congo hemorrhagic fever virus, and filovirus. The majority of these PVSs are genomic fragments of hosts including humans and bacteria. However, they cannot simply be regarded as foreign contaminants, since parts of them are results of natural occurrence or artificial engineering of viruses. Nevertheless, they severely disturb such sequence-based analyses as genome annotation, taxonomic assignment, and virome profiling. Therefore, we provide a clean version of the eukaryotic viral reference data set by the removal of these PVSs, which allows more accurate virome analysis with less time consumed than with other comprehensive databases. IMPORTANCE High-throughput sequencing-based viromics highly depends on reference databases, but foreign contamination is widespread in public databases and often leads to confusing and even wrong conclusions in genomic analysis and viromic profiling. To address this issue, we systematically detected and identified the contamination in the largest viral sequence collections of GenBank and UniProt based on a stringent scrutiny pipeline. We found hundreds of PVSs that are related to hosts, vectors, and laboratory components. By the removal of them, the resulting data set greatly improves the accuracy and efficiency of eukaryotic virome profiling. These results refresh our knowledge of the type and origin of PVSs and also have warning implications for viromic analysis. Viromic practitioners should be aware of these problems caused by PVSs and need to realize that a careful review of bioinformatic results is necessary for a reliable conclusion.

screening of other databases (6)(7)(8)(9)(10)24), we concluded that PVSs are related to hosts, vectors, laboratory components, and misclassification. Thus, we established a consecutive scrutiny pipeline, which is composed of five parts and managed using Snakemake (see Fig. S1 in the supplemental material). Because we aim at diagnosis of viral diseases and discovery of eukaryotic viruses, the first preliminary filtration step removes sequences of viruses infecting bacteria, archaea, fungi, or microorganisms or those shorter than 200 bp. The second step is host genome scrutiny, enabling us to detect sequences with origins of 28 mammalian, avian, and arthropod species, covering humans, domestic animals, and natural hosts and vectors, which are the predominant targets in virome analyses. The following step is to detect sequences derived from backbones or nonviral functional cassettes of vectors using two rounds of searches. Moreover, we noted that some viral reference sequences have incorrect lineage definitions, especially at high taxonomic levels, which is an important source of misannotation in virome analyses. Thus, we employed an exhaustive all-against-all comparison to identify misclassified sequences. Finally, we close the pipeline with multiple checks to detect LCD sequences using viral metagenomic annotation of 15 raw data sets that originated from different viromic studies spanning humans and other animals throughout the world. In order to avoid sacrificing the database's representativeness, we used very strict criteria in these steps to prevent sequences from being misclassified as PVSs. For example, a blast search is considered positive in a certain step only if the alignment achieves an E value of #1e250 and is longer than 500 (see Materials and Methods).
Host sequences are predominant in PVSs. The viral division (gbvrl) of GenBank is the largest resource of eukaryotic viral nt sequences and has been widely used in virological research, even to construct other specialized subdatabases (25,26), from which the Viral Genome Resources are derived to serve as a set of high-quality curated viral reference genomes and their validated genomic neighbors but lacking the full spectrum of viral diversity (27). As of 4 March 2021, gbvrl and the Viral Genome Resources have archived 3,316,373 and 288,226 nt sequences, respectively. They overlapped by 263,895 sequences; hence, we added the remaining 24,331 sequences of the Viral Genome Resources into gbvrl, which resulted in a preliminary data set (PDS) of 3,340,704 sequences. The PDS was subjected to PVS detection using the consecutive scrutiny pipeline. The preliminary filtration removed 91,549 sequences that contribute little to eukaryotic virome profiling of animals. After four rounds of scrutiny, we identified and deleted 766 nt PVSs from PDS (Data Set S1).
Regarding viral aa reference sequences, we retrieved all sequences under the Taxonomy of Viruses in UniProtKB (version 2021_03). UniProtKB is mainly based on the translation of genome sequences submitted to the International Nucleotide Sequences Database Collaboration (INSDC) source databases and also supplemented by genomes sequenced and/or annotated by other academic groups, making it the most comprehensive set of protein sequences (28). The aa sequence scrutiny was summarized in Data Set S1. In general, UniProt aa sequences have fewer PVSs than GenBank nt sequences, in that translation itself is a recognized validation method of viral genomes. Furthermore, foreign insertion often occurs as a flanking sequence in the untranslated region at the terminus of nt sequence. Finally, a total of 267 PVSs were detected, with most being counterparts in nt scrutiny (Data Set S1); hence, they will not be discussed below.
The causation of PVS occurrence: natural versus artificial. Contamination is very common in nonviral public databases and is caused by the foreign DNA of other species and even species from other kingdoms (6-10). However, we found that these PVSs cannot simply be considered foreign contaminants but rather should be ascribed to natural (na), intentionally artificial (ia), and unintentionally artificial (ua) occurrences.
(i) Natural occurrence. Some na-PVSs are naturally acquired by viruses in the process of replication, which is essential for certain viruses to gain new features. Bovine viral diarrhea virus (BVDV) is a worldwide-distributed pathogen and can cause severe consequences to cattle and sheep (29). Almost all PVSs within the family Flaviviridae are inserts of bovine hybrid ribosomal S27a and ubiquitin sequences into the BVDV genomes ( Fig. 2A). The in-frame insertion of the host sequence into the NS3 gene is essential for the virus to gain cytopathogenicity in cell culture (30). Hepatitis E virus (HEV) is hard to culture using cell systems; the integration of a short piece of human S17 ribosomal protein fragment into the hypervariable region of the HEV genome enables some variants to grow in HepG2/C3A cells (Data Set S1) (31). Besides host sequences, genomic fragments of other viral families can also integrate into some viral genomes, particularly during coinfection of multiple viruses. For some large DNA viruses, viral DNA replicates within the cellular nucleus or cytoplasm, providing an opportunity for a retrovirus to be integrated into a viral genome. Thus, avian retrovirus was shown to be integrated into the genome of Marek's disease virus, an avian herpesvirus (Data Set S1) (32). We also detected reticuloendotheliosis virus sequences of various lengths, even nearly full length, integrated into the genomes of some fowlpox viruses (Fig. 2B), which likely enhanced the pathogenicity of the viruses (33,34). Interfamily recombination can also occur in RNA viruses. A betacoronavirus detected in bats contained a unique gene at the 39 end of its genome that most likely originated from the p10 gene of a bat orthoreovirus, a gene that can induce the formation of cell syncytia (35). Elimination of Viral Foreign Sequences mSystems (ii) Intentionally artificial occurrence. Some viral genomes are intentionally engineered to contain foreign sequences that might derive from nonviral artifacts or viruses of other families, by which these engineered viruses provide important tools to study viral infection, deliver foreign proteins, or even combat diseases. We found that most vector PVSs (87.2%), a few misclassification PVSs (n = 3), and no host PVSs are intentional artifacts (Data Set S1). Among ia-PVSs with vector origins, green fluorescent proteins are very common (41.5%) ( Fig. 3A and Data Set S1), and elements like neomycin phosphotransferase, mCherry, and firefly luciferase were also detected (Data Set S1). The three misclassified ia-PVSs are all associated with avian paramyxovirus (Data Set S1). These are actually artificial recombinants designed to serve as vaccine vectors to combat avian influenza (Fig. 3B) (36).
(iii) Unintentionally artificial occurrence. The ua-PVSs are technically contaminating sequences but are unintentionally introduced into viral genomes. They are widely distributed in host, vector, laboratory component, and misclassification PVSs (Data Set S1). The host-origin ua-PVSs can be full-length sequences, e.g., a 399-bp-long human mRNA was erroneously defined as hepatitis C virus (Fig. 4A). Misassembly of HTS reads will result in chimeric ua-PVSs at the termini of a sequence, e.g., a 1,636-bp-long human sorting nexin 10 fragment was misassembled into the 39 terminus of segment M of a Crimean-Congo hemorrhagic fever orthonairovirus (CCHFV) (Fig. 4B). As to ua-PVSs of vector origin, we found two short stealth virus sequences that are actually vector backbones (Data Set S1). Through multiple checks using viral metagenomes from different hosts, we found some viral reference sequences sharing .99% nt identities Elimination of Viral Foreign Sequences mSystems with viromic contigs of different host species. Viruses harbored by different host species are usually distinct from each other due to the genetic adaptation to a specific host species. If a virus is detected in multiple host species, it should be particularly noted whether it is the result of cross-species transmission or simply foreign contamination. Further verification showed that these reference sequences are all nonviral but are genomic fragments of bacteria (Data Set S1). For example, a bluetongue virus sequence (AY397620) frequently encountered in our viral metagenomic analyses is a Mycoplasma bovis chromosomal sequence (Data Set S1). Cross-family misclassification can occur between eukaryotic viral families, even between eukaryotic and prokaryotic viral families. Three sequences wrapping circovirusfeatured rep and cap genes should be classified into the family Circoviridae but are defined as dependoparvoviruses within the family Parvoviridae (Fig. S2). A 558-bp-long sapovirus sequence (AB212270) defined within the family Caliciviridae actually originated from a Salmonella phage (Data Set S1). If a viral sequence is highly novel with very low similarity to known references, it would be likely misclassified at the family level. A 4,047-bp-long sequence recovered from a bird metagenome was defined as Parvoviridae sp., but it had very few blastn hits in the nt database and several blastx hits against major capsid proteins of microviruses (Fig. S3A). Profile comparison showed that one of its proteins is perfectly matched to the capsid of microviruses with a probability of 100% (Fig. S3B). Accordingly, it should be classified as a bacteriophage rather than a parvovirus.
Although these na-and ia-PVSs are essentially different from the foreign contaminants in other databases, they can, together with ua-PVSs, result in misleading conclusions in Elimination of Viral Foreign Sequences mSystems applications like virus detection, virome analysis, and taxonomic assignment. Therefore, we deleted or trimmed these PVSs from the PDS to build a clean version of the viral reference data set. However, the resulting data set is still redundant, with a high level of identical sequences. Thus, we removed the redundancy at 99% identity and 90% coverage, which downsized the nt and aa data sets for ;6 and ;3 times, respectively (Fig. S1). Functional enhancement of the data set. We built this data set primarily for virome-based virus detection. However, the HTS-based viral metagenome is prone to be contaminated by LCD viral sequences and vector-derived viral functional cassettes (11,12,16,37), which are also recorded in GenBank. This issue cannot be addressed by removing them from the reference database because those sequences would still be annotated by their genetically close relatives. Thus, we added these risk sequences into the data set by labeling them "LCD" or "Vector," which can provide a warning; if a query shows extremely high similarity to them, it should be concerning whether the sample is contaminated by exogenous LCD sequences or vectors (12). Besides, attenuated viruses are widely used in human and animal vaccinations to combat infectious diseases. It is critical to distinguish them from field strains in clinical diagnosis. Therefore, we enhanced the data set function by the addition of labeled LCD (n = 155), viral functional cassette (n = 79), and vaccine (n = 40) sequences into the nt branch (Data Set S2). Vaccine sequences added here covered 15 attenuated viruses commonly used in humans and animals against mumps, equine infectious anemia, porcine epidemic diarrhea, etc. (Data Set S2). By such enhancement, the final eukaryotic viral reference data set (EVRD) was achieved, with the nt and aa sequences archived in the Elimination of Viral Foreign Sequences mSystems EVRD-nt and EVRD-aa branches, respectively. EVRD-nt has 558,638 sequences with an average length of 2,943 bp covering 117 families, while EVRD-aa catalogs 1,256,089 sequences from 115 families with an average length of 371 aa. At the time of manuscript preparation, we have updated the data set into version 2.0 with addition of 34,037 and 221,426 new sequences into EVRD-nt and EVRD-aa, respectively. These new sequences covered 111 families with the majority being retroviruses, picornaviruses, and reoviruses. EVRD improves the accuracy and efficiency of virome analysis. The performance of EVRD in virome analysis was evaluated by comparison of its accuracy, coverage, and efficiency with GenBank (for nt) and UniProt (for aa) virus branches and RVDB (v21.0). RVDB is a reference viral database that provides a broad representation of different virus species from eukaryotes by including all viral, virus-like, and virus-related sequences (26). The detailed comparisons were summarized in Text S1. Briefly, using nine viral metagenomes of pigs, bats, and humans, the blastn annotation revealed that the majority (88.1%) of virus-like reads (VLRs) were coannotated by the three databases (Fig. 5), but EVRD massively reduced false annotations to calicivirus, BVDV, pneumovirus, etc. (Fig. 5). Annotation using contigs at the nt and aa levels also showed that EVRD prevented misannotation to calicivirus, reovirus, parvovirus, etc. (Fig. S4). These results indicated that EVRD does not sacrifice the detection spectrum of eukaryotic viruses but rather greatly improves the specificity and accuracy of viromic annotation (Text S1). Notably, we detected cocirculation of field and vaccine strains of porcine reproductive and respiratory syndrome viruses in viral metagenome AH (Text S1 and Fig. S5), which should be especially concerning, since new viruses could be generated through recombination between field viruses and vaccine strains. Viromic annotation is quite consuming of time and computing resources. A small-scale reference database can save the analytic time and minimize the computing resources. With an entry-level platform, analyses of reads or contigs at nt or aa levels using EVRD were 1.8 to 3.3 and 1.9 to 3.2 times faster than using GenBank/UniProt and RVDB, respectively (Fig. S6), indicating that EVRD is more time-saving and easier to configure.

DISCUSSION
A high-quality database is critical in bioinformatic analyses in biology, biotechnology, and medicine. To maintain the quality of the database, the National Center for Biotechnology Information applies several protocols to screen submitted sequences for contamination, such as using the program VecScreen and cross-checking using BLAST search (7,38). However, such measures are insufficient to provide adequate quality assurance. The original submitters should carefully check the quality of the sequence(s) before they submit it. As to viral sequences, many aspects can be considered to improve their quality. First, the viral genome should be thoroughly annotated to indicate all foreign sequences that occur naturally or are introduced intentionally, so that users can understand that they are related to evolutionary processes or biomedical use, as well as take note of their heterogeneity. Second, although HTS-based viral metagenomics has greatly promoted the rate of new virus discovery, there is an inevitable issue in which viral genomes can be misassembled to include a foreign insert(s) and a misplaced gene(s) (39). Thus, those new sequences should be validated by multiple means, such as PCR amplification and Sanger sequencing. Third, the emerging profile-based and machine learning-based virome annotation algorithms allow us to explore those remote viruses that show very limited similarities to known references (40,41), but it would be difficult to accurately assign their taxonomic lineages. In such cases, classification based on multiple criteria such as virus hallmark gene annotation and genomic organization is helpful. Virus hallmark genes encode very conserved proteins that are usually shared by broad ranges of viruses, which can be used as indicators in high-rank taxonomic classification (42), while viruses within a low-rank taxon have common genomic organization, which helps to determine taxa below family level. Last but not least, users should be aware of these problems caused by PVSs FIG 5 Comparison of the virus-like read numbers in nine viromic data sets annotated by blastn searches against EVRD-nt (highlighted in orange), GenBank, and RVDB-nt. Viral families are divided into categories of "Shared," "EVRD," and "Other," corresponding to families that are coannotated by the three reference databases, not annotated by EVRD in certain data sets, and annotated by one or two reference databases in certain data sets, respectively.
Elimination of Viral Foreign Sequences mSystems when using databases and need to realize that a careful review of bioinformatic results is necessary for a reliable conclusion. In other databases, those foreign sequences are considered contaminating ones, among which human sequences contribute a substantial part (6-10). It is not surprising that host sequences are the predominant PVSs because host cells provide the obligatory venue for viruses to replicate. However, unlike organisms in other kingdoms, viruses are characterized by unique and diverse replication strategies, with some capable of exchanging genomic fragments with their hosts (43), and they are also easily genetically modified to contain foreign DNA for various biomedical purposes. Thus, these host PVSs cannot be arbitrarily considered foreign contaminants, since some of them are actually natural products or intentional artifacts, which explains why the foreign sequences in viral databases are much more complicated and more difficult to differentiate than those in databases of other organisms (7). However, no matter how they originate, these PVSs can cause a variety of problems in downstream studies, particularly for evolutionary analysis and genome annotation.
EVRD is maintained in FASTA format with informative identifiers. It can be easily localized and configured according to any similarity-based sequence analytical software, such as BLAST suite for sequence search or MMSeqs2 for sequence cluster. EVRD provides a valuable resource for, but not limited to, these application scenarios below. Accurate determination of causative agents is a priority in clinical diagnosis of EIDs. EVRD can be widely used in virome-based viral disease diagnosis since it avoids those false positives caused by misannotation of PVSs. EVRD can also improve the taxonomic classification of viral sequences in the assessment of virus diversity (39). In such analysis, viral contigs need to be clustered with reference sequences, but these PVSs, specifically the cross-family misclassified ones, will disturb the boundary of virus clusters and even result in incorrect taxonomic classification. In addition, multiple sequence alignments (MSAs) are prone to being corrupted by PVSs. The refined EVRD sequences can help build high-quality MSAs that are the basis of profiles of clustered sequences, thus favoring the exploration of remote viruses.
Besides utilizing a high-quality reference database, other measures can be considered to improve the reliability of viromics-based virus detection. According to different purposes, we should choose a reasonable bioinformatic pipeline. Annotation using reads provides richer information than using contigs and hence is helpful to capture ultralow-abundance viruses in virus detection (44,45). But sequence completeness is a priority in viral ecology, so assembly is preferentially performed before annotation (39). The criterion used to determine a viral sequence has a substantial impact on virome annotation (see Text S1 in the supplemental material). Regarding read, the criterion is mainly based on E value, but the alignment length is also an important factor to help increase the confidence level of annotation. Besides E value and length, the requirement of a minimum gene number has been widely considered in contig annotation (39). The quality of assemblies should be seriously considered in contig annotation. There are many means to improve the assembly quality, such as choosing suitable software (46), employing a rational sample treatment protocol (47), and reducing the bias induced by random amplification (48). A classification of host and other microorganism reads prior to de novo assembly could help reduce chimeric contigs (47). Besides, a final check provides an additional guarantee for high-quality annotation (49), e.g., host contamination should be checked and eliminated as much as possible. Prokaryotic host contamination can be determined using CheckV, but a different strategy is needed to deal with eukaryotic host contamination (49,50).
In conclusion, we detected these PVSs hidden in the largest viral nt and aa databases using a systematic scrutiny pipeline, followed by careful identification of their types and origins, which not only refreshes our knowledge of the occurrence of PVSs but also provides a clean version of the eukaryotic viral reference data set. The scrutiny pipeline is designed for eukaryotic viruses but is not suitable for prokaryotic ones. Bacteriophages, specifically lysogenic ones, frequently exchange their genomic elements with their hosts, making it very difficult to differentiate foreign contamination from genomic integration. But it is feasible to identify host contamination in phage genomes by creating and querying a distinctive bacterial gene set. Besides, it is important to note that we excluded those large DNA viruses infecting eukaryotic microorganisms, due to their extraordinarily large and complicated genomes and the lack of evidence that they cause diseases in vertebrates (51)(52)(53). Though we have deleted hundreds of PVSs of vertebrate-infecting large DNA viruses from families like Herpesviridae and Poxviridae, there are still some ambiguous sequences that can be treated as host PVSs if using loose criteria. Thus, annotations to these viruses using EVRD should still be verified with caution. Additionally, these tagged warning sequences in EVRD are very useful, but they are partial and represent only the sequences we have searched so far. We will keep the database updated with new advances in this regard.

MATERIALS AND METHODS
Scrutiny pipeline for nucleotide sequences. (i) Preliminary filtration. We first generated the taxonomic lineages of all sequences and then removed those lineages infecting bacteria, archaea, fungi, and eukaryotic microorganisms using the relationship of virus and host recorded in the ViralZone database (54). In addition, as there are a large number of sequences that cannot be assigned to a complete lineage, we searched their definition using keywords and removed the sequences related to prokaryotic and environmental viruses and metagenomes, such as bacteriophage/phage, environmental and uncultured viruses, and ameba viruses. The gbvrl division also deposits a large number of sequences with lengths of #200 bp, which are highly similar to the longer ones and contribute little to genomic analysis; hence, they were also removed.
(ii) Host genome scrutiny. In this part, fragments of host genomes in the remaining sequences of the PDS were detected. Genomic assemblies of human (n = 1), pig (n = 1), bat (n = 7), rodent (n = 2), arthropod (n = 11), cattle (n = 1), dog (n = 1), cat (n = 1), sheep (n = 1), chicken (n = 1), and mallard (n = 1) were used to perform a blastn search against these sequences with a maximum of 1,000 subjects to show alignments (length of $150 and identity of $85%). Retroviruses can infect almost all vertebrates, resulting in thousands of loci of retroviral sequences in vertebrate genomes (55). Here, we did not challenge the known ambiguity of retroviruses, and hence, hits to retroviruses were not considered. The aligned sequences of the subject were extracted and subjected to blastn search against the nt database to further validate their identities. The top 100 hits of each sequence were kept, and if $80% of these hits were annotated as nonviral, the aligned sequence was considered problematic. The original sequence was removed from the PDS if its problematic part comprised $80% of its length or was trimmed by deleting the problematic parts. Such a threshold was also applied to the following treatments. Such treatment was iterated until no host genomic fragments were found.
(iii) Vector sequence scrutiny. To detect PVSs derived from backbones or functional cassettes of vectors, the UniVec database and sequences of $1,000 bp under the GenBank taxonomy of vectors (uid: 29278) were downloaded. As vectors have many functional cassettes originating from viruses, such as simian virus 40 (SV40) and cytomegalovirus (CMV) promoters and retroviral gag and pol elements, these vector-originating PVSs in the PDS were carefully detected and examined using the following procedure to prevent any erroneous deletions of genuine viral sequences. We generated a nonviral protein core (NVPC) that consists of nonviral expression elements (n = 13,287) originating in vectors (see Text S1 in the supplemental material). Sequences in the PDS were subjected to a blastx search against the NVPC using Diamond with those showing $99% similarity over an alignment of $60 aa to subjects being pruned. In addition, UniVec was used to identify adapters, linkers, and primers that are often used to clone sequences. The remaining sequences in the PDS were further scrutinized using the procedure introduced under "(ii) Host genome scrutiny" with the same criteria. Briefly, these downloaded vector sequences were used as queries to search for possible subjects in the PDS using blastn. Hits in the PDS were further validated by blastn searches against the nt database. After the removal of those vectororiginating sequences, the rest of the PDS was reexamined until no vector sequences existed.
(iv) Misclassification scrutiny. Erroneous taxonomic annotation of viral sequences was detected by all-against-all blastn searches with a maximum of 1,000 subjects to show alignments. We found that there are a large number of sequences with correct taxonomic annotation showing intrafamily crossspecies/genus blastn hits, such as Betacoronavirus/Gammacoronavirus within the family Coronaviridae, Tetraparvovirus/Protoparvovirus of the family Parvoviridae, and Circovirus/Cyclovirus within the family Circoviridae, which were likely ascribed to high similarity between species/genera. Hence, we inspected annotation at the family level. Here, we defined that a blastn hit is significant if its E value is #1e250 and length is $500. If the proportion of alignments that were generated by a query against subjects of a different family to all alignments of the query was $80%, the query was considered to be possibly misclassified and was further subjected to genomic organization identification, in which if the genomic organization of the query was not typical of the features its defined taxonomic lineage should have had, the query was truly misclassified and was removed from the PDS. During treatment, we noted that some sequences had a few alignments (usually #10) that showed #80% similarities with subjects of different families, and we kept their original annotations since there are not enough references in GenBank to determine their true taxonomic lineages.

Elimination of Viral Foreign Sequences mSystems
(v) LCD sequence detection. Previous studies showed that some viral contaminant sequences are highly prevalent in cross-host HTS-based viromic data, which might be linked to biological or synthetic products (11,12,16). To examine whether cross-host sequences exist in the database, the remaining sequences in the PDS were subjected to a cross-check of viral metagenomes. A total of 15 viromic raw data sets covering humans, bats, ticks, rodents, bovine, pigs, and avians were downloaded from the SRA and respectively de novo assembled. Contigs of $1,000 bp were subjected to blastn search against the PDS with a maximum of 1,000 subjects to show alignments. If a subject was matched by contigs from viromic data sets of two or more different hosts with alignment of $150 bp and identity of $80%, it was classified as a suspicious sequence and further validated by blastn search against the nt database. If a suspicious sequence was annotated to nonviral species by blastn searches against the nt database, it was considered a truly exogenous contaminating sequence and removed from the PDS. However, if a suspicious sequence was still annotated to virus and shared 99% nt identities with viromic contigs of two or more different hosts, it was considered a truly viral sequence but one probably originating from laboratory-component-derived viral sequence contamination and hence was retained in the PDS but was tagged as LCD. The remaining suspicious sequences were passed and kept in the PDS.
Scrutiny pipeline for viral protein sequences. The protein sequences were retrieved from the UniProt virus division. We first checked their representativeness. In case there were any coding regions not annotated by the original submitters, all proteins of the PDS nt sequences prior to filtration were de novo predicted using Prodigal v2.6.3 with meta mode. Proteins of $50 aa were blastp searched against the UniProt viral division (E value of #1e210 and sequence identity of $90), and the results revealed that the UniProt viral division has high representativeness with 99.6% consistency with the prediction of GenBank viruses. Therefore, we simply used UniProt sequences to screen PVSs as described for scrutiny of nt sequences with minor modification (Text S1). Although these nt and aa PVSs were independently detected, we compared the counterparts of aa PVSs in their nt sequences with these nt PVSs to further confirm that these sequences are problematic.
EVRD finalization. After the above scrutiny, the sequences in the PDS were still very redundant, and hence, a deredundance procedure was applied to downsize the PDS. Clustering of viral nt and aa sequences was performed using MMseq2 (56) with a sequence similarity threshold of 0.99 and 90% coverage of the short sequence. The LCD sequences reported previously (12,16) and identified here were tagged using "LCD" and added into the PDS as risk sequences. To distinguish viral functional cassettes from true virus sequences, the sequences corresponding to the regulatory classes of promoter, terminator, and enhancer and/or with notes containing the word "virus" were extracted from vectors and subjected to blastn search against the nonredundant PDS. The vector sequences verified to be viral were dereplicated and also added into the PDS with tags of "Vector." In addition, we collected vaccine strains commonly used in humans and other animals by searching publications or by personal communication (Data Set S2). These vaccine nt sequences were also tagged using "Vaccine" and added into the PDS.
Performance evaluation of EVRD. We evaluated the performance of EVRD by comparison with gbvrl (for nt), UniProt (for aa), and RVDB at the read and contig levels. A total of nine viral metagenomic raw data sets (Table S1) were subjected to host genome removal using Bowtie2 (v2.4.1) with sensitive mode and then taxonomically classified using Kraken2 (v2.0.9-beta) to remove bacterial, archaeal, and fungal reads. First, the unassigned reads were directly blastn (E value of #10e25 and length of $120) and blastx (E value of #10e25 and length of $40) searched against these databases. Then, they were de novo assembled using megahit (v1.2.9). Contigs of $1,000 bp were retained for blastn (v2.10.0) and Diamond blastx (v0.9.35) search against the nt and aa reference databases, respectively. The blastn hit of a contig to a subject with one alignment with an E value of #10e210 and a length of $450 or two or more alignments with an E value of #10e25 and a length of $150 was considered positive, and the blastx hit to a subject was recognized as positive if it had one alignment with an E value of #10e210 and a length of $150 or two or more alignments with an E value of #10e25 and a length of $50. The positive reads and contigs were further verified by blastn/x search against nt/nr databases (18). All blast searches were performed using 12 Â 86_64 central processing units (CPUs) of an Inter Xeon Gold 2.660-GHz processor. To detect if a data set contained warning sequences, we defined a rigorous cutoff, i.e., if a sequence hit a tagged subject with an identity of $99% and coverage of the query of $90% by blastn search against EVRD-nt, it was considered to be a warning sequence.
Data availability. All data used here were downloaded from relevant databases. The accession numbers of viral metagenomic raw data and genome identifiers (IDs) are summarized in Fig. S1 and Text S1. The key intermediate data (NVPC) and essential codes are available from http://github.com/BH-Lab/ EVRD. The code for viromic analysis using EVRD has been published elsewhere (18). EVRD reported here is scheduled to annually update, and the current version is v2.0, which is freely accessible at http://cvri .caas.cn/kxyj/yjfx/bfbd/EVRD/index.htm.

SUPPLEMENTAL MATERIAL
Supplemental material is available online only. TEXT S1, DOCX file, 0.02 MB.