Next Article in Journal
Semi-Supervised Pipeline for Autonomous Annotation of SARS-CoV-2 Genomes
Previous Article in Journal
The Impact of COVID-19 on Mortality in Spain: Monitoring Excess Mortality (MoMo) and the Surveillance of Confirmed COVID-19 Deaths
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 12
10.3390/v13122425
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

Virome Characterization in Commercial Bovine Serum Batches—A Potentially Needed Testing Strategy for Biological Products

by
Willian P. Paim
1,2,
Mayara F. Maggioli
1,
Shollie M. Falkenberg
3,
Akhilesh Ramachandran
4,
Matheus N. Weber
5,
Cláudio W. Canal
2 and
Fernando V. Bauermann
1,*
1
Department of Veterinary Pathobiology, College of Veterinary Medicine, Oklahoma State University (OSU), Stillwater, OK 74078, USA
2
Laboratório de Virologia, Faculdade de Veterinária, Universidade Federal do Rio Grande do Sul (UFRGS), Rio Grande do Sul, Porto Alegre 90040-060, Brazil
3
Ruminant Diseases and Immunology Research, Agricultural Research Service, USA Department of Agriculture, Ames, IA 50010, USA
4
Molecular Diagnostic Laboratory, Oklahoma Animal Disease Diagnostic Laboratory, Oklahoma State University (OSU), Stillwater, OK 74078, USA
5
Laboratório de Microbiologia Molecular, Instituto de Ciências da Saúde, Universidade Feevale, Rio Grande do Sul, Novo Hamburgo 93510-235, Brazil
*
Author to whom correspondence should be addressed.
Viruses 2021, 13(12), 2425; https://doi.org/10.3390/v13122425
Submission received: 1 November 2021 / Revised: 21 November 2021 / Accepted: 30 November 2021 / Published: 3 December 2021
(This article belongs to the Topic Veterinary Infectious Diseases)
Browse Figures
Versions Notes

Abstract

:
Bovine serum has been widely used as a universal supplement in culture media and other applications, including the manufacture of biological products and the production of synthetic meat. Currently, commercial bovine serum is tested for possible viral contaminants following regional guidelines. Regulatory agencies’ established tests focused on detecting selected animal origin viruses and are based on virus isolation, immunofluorescence, and hemadsorption assays. However, these tests may fail to detect new or emerging viruses in biological products. High-throughput sequencing is a powerful option since no prior knowledge of the viral targets is required. In the present study, we evaluate the virome of seven commercial batches of bovine serum from Mexico (one batch), New Zealand (two batches), and the United States (four batches) using a specific preparation and enrichment method for pooled samples and sequencing using an Illumina platform. A variety of circular replicase-encoding single-stranded (CRESS) DNA families (Genomoviridae, Circoviridae, and Smacoviridae) was identified. Additionally, CrAssphage, a recently discovered group of bacteriophage correlated with fecal contamination, was identified in 85% of the tested batches. Furthermore, sequences representing viral families with single-stranded DNA (Parvoviridae), double-stranded DNA (Polyomaviridae and Adenoviridae), single-stranded RNA (Flaviviridae, Picornaviridae, and Retroviridae), and double-stranded RNA (Reoviridae) were identified. These results support that high-throughput sequencing associated with viral enrichment is a robust tool and should be considered an additional layer of safety when testing pooled biologicals to detect viral contaminants overlooked by the current testing protocols.

1. Introduction

Bovine serum has been used as a universal supplement in cell culture applications [1]. The cell culture uses include in-vitro diagnostics, vaccine production, gene therapy, stem cell therapy, and, more recently, the production of synthetic meat [2,3,4]. Different categories of bovine serum are commercially available, including Fetal Bovine Serum (FBS), Newborn Calf Serum (NCS), Bovine Calf Serum (BCS), and Donor Bovine Serum (DBS).
Bovine origin products used to manufacture biologicals are currently tested for viral contaminants following regional guidelines. For instance, in the United States (USA), the test for adventitious viruses in cattle serum is performed in compliance with the United States Department of Agriculture (USDA) regulations for biologics products as specified in the Code of Federal Regulations (CFR9), Title 9, Part 113.53. The USDA regulations require the testing of bovine serum for Bovine Adenovirus (BAV), Bovine Herpesvirus (BoHV), Bovine Parvovirus (BPV), Bovine Respiratory Syncytial Virus (BRSV), Bovine Viral Diarrhea Virus (BVDV), Bluetongue Virus (BTV), Parainfluenza Virus 3 (PI3), Rabies Virus (RABV), and Vesicular Stomatitis Virus (VSV). The screening tests are mainly conducted using traditional virology assays, including virus isolation, immunofluorescence, and hemadsorption. Similarly, in Europe, the use of bovine serum for the manufacture of human medicinal products with a biological origin is regulated by the European Medicines Agency (EMA/CHMP/BWP/457920/2012 Rev.1). The veterinary products in Europe are regulated by the EMA-CVMP-743-00-Rev.2.
The EMA guidelines mostly rely on virus isolation in cell culture. Additionally, manufacturers are encouraged to test for specific viruses using nucleic acid detection assays and investigate possible emerging viruses. However, specific guidelines and requirements toward this end are unclear. The lack of a systematic and global strategy to unbiasedly monitor for emerging viral threats creates opportunities for viral entry and dissemination within and among countries, with grave implications to animal and human health.
There is an extensive international trade of bovine serum [5]. While part of the serum batches is sterilized by gamma irradiation, the procedure is not required. Commercialized bovine serum is typically pooled serum, frequently composed of hundreds to thousands of individual samples. In the final product label, the described product origin is based on the country where processing and packaging occurred (World Trade Organization–Rules of Origin, Article 3 {b}). The country processing the serum may not necessarily correspond to where the serum was first collected, or the origin of the animals used to collect the serum.
Sites collecting and processing bovine serum for commercialization may also handle products of various animal species, allowing the chance for contamination among products with viruses from unrelated animal species. This is especially concerning regarding emerging or high-impact animal disease dissemination. Another relevant concern is the presence of human and zoonotic viruses contaminating bovine serum. During collection and processing, humans are involved in the process, and the testing requirements do not consider the accidental or intentional risk for contaminating bovine serum batches with human viruses.
Second-generation sequencing technologies, including sequencing by synthesis and pyrosequencing platforms, provide the ability to detect known and unknown viruses in pooled samples [6,7,8,9]. These technologies associated with the enrichment of viral nucleic acid potentially increase the sequencing sensitivity. For instance, an evaluation of the virome in FBS samples used the Klenow enzyme for DNA amplification and characterized sequences originating from eukaryotic viruses [6]. On the other hand, multiple displacement amplification (MDA) was used to identify specific viruses and evaluate the test sensitivity in biological samples [10,11]. In the present study, we characterized the viral population in non-gamma irradiated bovine serum batches from Mexico, the USA, and New Zealand. The serum batches were commercially available in the USA and were subjected to a combination of steps previously described to enrich viral nucleic acid in the samples before sequencing [6,7,9,10,11,12]. Results demonstrated the presence of sequences related to several eukaryotic viruses from a wide range of animal species and a great diversity of phages, including members of the recently discovered group of cross-assembly phage (CrAssphage).

2. Materials and Methods

2.1. Samples of the Commercial Bovine Serum Batches

Seven batches of non-irradiated commercial bovine serum were purchased from three different manufacturers. One FBS from Mexico (FBS1 Mex), two NCS batches from New Zealand (NCS2 NZ and NCS3 NZ), and four batches from the United States, including one FBS (FBS2 USA), two NCS (NCS1 USA and NCS4 USA), and one BCS (BCS1 USA). The batches of commercial bovine serum were stored at −80 °C until processing. In addition, a non-related sample composed of peripheral blood mononuclear cells (PBMC) from 15 two-week-old calves was added. The sample collection was approved by the Institutional Animal Care and Use Committee (IACUC) protocol number IACUC-20-14-STW.

2.2. Sample Preparation and Sequencing

A total of 50 mL of each commercial batch of bovine serum was centrifuged at 2000× g for 30 min at 4 °C. Prior to centrifugation, the calf PBMC pooled sample was submitted to a freeze and thaw cycle and then diluted to 50 mL in PBS. The supernatants were filtered using a 0.45 μm pore size membrane (Millipore, Darmstadt, Germany). Subsequently, samples were ultracentrifuged on a 25% sucrose cushion at 100,000× g for 3 h at 4 °C in a Sorvall AH629 rotor (Thermo Fisher Scientific, Waltham, MA, USA). The pellet containing the viral particles was incubated for 1.5 h with DNase and RNase enzymes (Thermo Fisher Scientific) as previously described [13]. Viral RNA was isolated using TRIzol™ LS Reagent (Thermo Fisher Scientific), whereas the DNA extraction employed a phenol-chloroform protocol [14].
The viral DNA was enriched through multiple displacement amplification (MDA), using a GenomePlex® Complete Whole Genome Amplification (WGA) Kit (Sigma-Aldrich, St. Louis, MO, USA). Following the manufacturer’s recommendations, the viral RNA was reverse-transcribed and further enriched using the SeqPlex RNA Amplification Kit (Sigma-Aldrich). The enriched DNA products were pooled in equimolar amounts and purified using the PureLink™ Quick Gel Extraction and PCR Purification Combo Kit (Thermo Fisher Scientific). The quality and quantity of the DNA were assessed through spectrophotometry and fluorometry performed with NanoDrop™ (Thermo Fisher Scientific) and Qubit™ 4 (Thermo Fisher Scientific). Using the Nextera XT DNA Library Preparation Kit, DNA libraries were prepared with 1 ng of purified DNA. Sequencing was conducted using the Illumina iSeq 100 System using the Illumina i1 reagent kit (2 × 150 paired-end reads).

2.3. Bioinformatic Analysis

The quality of the sequences generated was evaluated using the FastQC tool. High-quality sequences (Phred score ≥ Q30) were trimmed on the Illumina iSeq 100 System at the time the FASTQ was generated. The data was de novo assembled on BaseSpace Cloud (Illumina) with the metaSPAdes genome assembler (Version 3.0). All assemblies were confirmed by mapping reads to contigs as previously described [7] and using Geneious Prime software (Version 2020.2.1). The assembled contigs were examined in search of similarities to known sequences with BLASTN software using MegaBLAST [15]. Sequences with E-values of ≤10−5 were classified as likely originating from a eukaryotic virus, bacteria, phages, eukaryote, or unknown, based on the taxonomic origin of the sequence with the best E-value as previously described [16]. Open reading frames (ORF) predictions and genome annotations of the near-full-length genomes were performed using Geneious Prime software. The alignment of contigs and reads to a reference genome was used to generate the near near-full-length genomes using Geneious Prime software (Version 2020.2.1). Gene comparisons were performed with BLASTN (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 25 January 2021). Representative sequences of viruses belonging to the families Parvoviridae and Polyomaviridae were obtained from GenBank and aligned with the sequences identified in the present study using MUSCLE software [17]. Phylogenetic trees were constructed using MEGA6 [18], applying the Neighbor-Joining (NJ) method through the p-distance model for Parvoviridae analysis and the Maximum-Likelihood (ML) method, General Time Reversible (GTR) model to Polyomaviridae family. All analyses were conducted with 1000 bootstrap replicates.

3. Results

3.1. Overview

Approximately 11,492,156 sequence reads were generated by metagenomic sequencing in the serum batches. From the total retrieved reads, about 12.7% were classified as viral sequences, whereas about 56.8% matched the bovine genome. After assembly, a total of 639,902 contig sequences were obtained; 5.9% were viral contigs. From these, 42.9% were assigned to eukaryotic viruses and 57.1% to bacteriophages. A taxonomic classification, including the number of reads recovered for each eukaryotic virus-related contig, is provided in Table 1. Sequences representing viral families with single-stranded DNA (ssDNA) (Parvoviridae), double-stranded DNA (dsDNA) (Polyomaviridae and Adenoviridae), and circular Rep-encoding single-stranded (CRESS) DNA (Genomoviridae, Circoviridae, and Smacoviridae) were observed. Additionally, sequences belonging to four RNA viral families were identified, including single-stranded RNA (ssRNA) (Flaviviridae, Picornaviridae, and Retroviridae) and double-stranded RNA (dsRNA) (Reoviridae) (Table 1). The complete coding genomes were obtained for Bovine Parvovirus 2 (BPV-2), Bovine Parvovirus 3 (BPV-3), Bosavirus (BosaV), Bovine Hokovirus 2 (BHoV-2) (Figure 1), and Bovine Polyomavirus 1 (BPyV-1) (Figure 2). Most of the viral sequences obtained shared a high degree of identity with known animal viruses. Moreover, sequences representing 12 viral families of bacteriophages were identified (Figure 3), including sequences related to CrAssphage (Table 2). The sequencing results of the PBMC from calves revealed the presence of eukaryotic virus-related contigs belonging to the families Parvoviridae, Herpesviridae, Retroviridae, Phycodnaviridae, Polydnaviridae, and Poxviridae. Additionally, bacteriophage related-sequences from the families Myoviridae, Podoviridae, and Siphoviridae were retrieved (Supplementary Table S1).

3.2. DNA Viral Families: Parvoviridae, Polyomaviridae, and Adenoviridae

A high diversity of sequences matching members of Parvoviridae were retrieved, including sequences representing four parvovirus genera: Erythroparvovirus, Copiparvovirus, Tetraparvovirus, and Bocaparvovirus. A total of 16,156 contig sequences closely related to the family Parvoviridae were observed in the tested samples (Table 1). Among these contig sequences, 11,287 were related to BPV-3 (Ungulate erythroparvovirus 1). Two complete coding genomes of BPV-3 were recovered by mapping the contigs and reads to a reference sequence. The BPV-3 genomes were 5100 nucleotides in length, with 58,148× and 10,207× average coverage. The genomes BPV-3 NCS2/NZ and BPV-3 BCS1/USA were deposited in GenBank under the accession numbers MZ502232 and MZ502233, respectively. The Parvoviridae family phylogeny is presented in Figure 1.
A total of 4342 contigs were closely related to BPV-2 (Ungulate copiparvovirus 1) (Table 1). Five complete coding genomes belonging to the Copiparvovirus genus were obtained (BPV-2 FBS2/USA, BPV-2 NCS1/USA, BPV-2 NCS3/NZ, BPV-2 NCS4/USA, and BPV-2 BCS1/USA). The sequences were deposited in GenBank under the accession numbers MZ502227, MZ502228, MZ502229, MZ502230, and MZ502231. The genome lengths were 5334, 5613, 5559, 5257, and 5613 nucleotides, and presented an average coverage of 94 ×, 38,275×, 2150×, 23×, and 10,394×, respectively. Sequences of BosaV (Ungulate copiparvovirus 5) were also described in all NCS samples. Two complete coding genomes of Copiparvovirus were retrieved (BosaV NCS1/USA and BosaV NCS3/NZ). The genomes were deposited in GenBank under the accession numbers MZ502234 and MZ502235. The genomes were 5222 and 5371 nucleotides long and had an average coverage of 105× and 541×.
Three serum batches (NCS1/USA, NCS4/USA, and BCS1/USA) showed 106 contig sequences related to bovine hokovirus (Table 1). Notably, nucleotide BLAST demonstrated that 12 contigs presented 87.8% to 100% identity with the Bovine Hokovirus 1 (BHoV-1) strain HK1 and Brazilian BHoV-1 isolate Guarapuava (GenBank accession numbers NC038898 and MG745679), while 94 contigs were related to the Chinese Bovine Hokovirus 2 (BHoV-2) strain HK-B38, and BHoV-2 strain BS-S13 (85.7–100% identity) (GenBank accession numbers JF504698 and KU172423) (Table 1). One complete Tetraparvovirus coding genome with 5097 nucleotides in length and 376× average coverage sequence was identified in the BCS1 sample. The BHoV-2 BCS1/USA genome sequence was deposited in GenBank under the accession number MZ502236.
Twenty-six contig sequences that were closely related to the Polyomaviridae family were retrieved from serum batches FBS1, FBS2, NCS3, and NCS4. One complete coding genome of BPyV-1 with 19× average coverage of the genome and 4697 nucleotides long was retrieved from FBS2/USA (GenBank accession number MZ520135). The phylogeny is shown in Figure 2.
Four contig sequences were closely related to the Adenoviridae family (Table 1). The batch NCS3/NZ showed two contig sequences that were 103 and 147 nucleotides long and 97.1–99.3% nucleotide identity with the Bovine Adenovirus 4 (BAdV-4) strain THT/62 identified in Hungary (GenBank accession number AF036092). The other two contig sequences were detected in serum batch BCS1/USA. The contig sequences were 280 and 282 nucleotides long and were related to the Bovine Adenovirus 6 (BAdV-6) strain 671130 from the Netherlands, with 98.9–97.8% nucleotide identity (GenBank accession number JQ345700).

3.3. RNA Viruses: Flaviviridae, Picornaviridae, Reoviridae, and Retroviridae

Four contig sequences were closely related to members of the Flaviviridae family, the Pestivirus A and B species (formerly classified as Bovine Viral Diarrhea Virus 1 (BVDV-1) and 2 (BVDV-2) species) (Table 1). The serum batches NCS1 and NCS4 contained three contig sequences related to the Bovine Kobuvirus (BKV; species Aichivirus B, genus Kobuvirus, family Picornaviridae) (Table 1).
A total of three contig sequences closely related to family Reoviridae were observed in serum batches NCS2/NZ (two contigs) and NCS3/NZ (one contig). From these, a total of two contig sequences were related to segments 1 and 8 (VP1 and NSP2) of the Bovine Rotavirus A (BRV-A), while one contig (NCS3/NZ) showed 96.1% nucleotide identity to segment 3 (VP3) of the Human Rotavirus (HRV) strain RVA/Human-wt/AUS/V585/2011/G10P (GenBank accession number JX567762). In addition, one contig sequence closely related to the Retroviridae family was observed in serum batch FBS1 from Mexico. The 206 nucleotides long contig had 99.5% nucleotide identity with the gag-pol polyprotein from the Moloney Murine Leukemia Virus (MoMLV) (GenBank accession number NC001501).

3.4. Circular Rep-Encoding Single-Stranded (CRESS) DNA Viruses-Genomoviridae, Circoviridae, and Smacoviridae

A total of 28 contig sequences closely related to CRESS-DNA viruses were retrieved. Fifteen contig sequences were related to ten CRESS-DNA viruses (feces associated gemycircularvirus 14 and 16, Gopherus associated genomovirus 1, Mongoose feces-associated gemycircularvirus A and B, Sewage-associated gemycircularvirus 3 and Genomoviridae sp.) belonging to the Genomoviridae family (Table 1). Twelve contig sequences (103–482 nucleotides long) were related to nine CRESS-DNA viruses in the Circoviridae family (Cyclovirus Equ1, Po-Circo-like virus 21, Po-Circo-like virus 22, Po-Circo-like virus 41, Po-Circo-like virus S20, Po-Circo-like virus GX19, Starling circovirus, Bo-Circo-like virus CH, and Circovirus sp. isolate PoCirV VIRES YN02 C2). The Smacoviridae family showed one contig sequence that was 151 nucleotides in length and shared 88.7% nucleotide identity to porcine associated porprismacovirus 10 isolate 49-Fec25-pig detected in NZ (GenBank accession number NC030126).

3.5. Bacteriophages

Viruses from the Myoviridae, Siphoviridae, and Podoviridae families represented the majority of the sequences recovered (Figure 3). A total of 21,566 contig sequences were assembled and closely related to bacteriophages. The Podoviridae family had 3361 related contig sequences. Among these contigs, 242 were closely related to CrAssphage and found in six out of the seven tested batches. The NCS serum batches demonstrated an increased number of CrAssphage sequences (Table 2).

4. Discussion

Commercial bovine serum is widely used as a supplement for cell culture in research, biotechnology, pharmaceutical manufacturing, and the production of synthetic meat [1,2,3,4]. The extensive international trade of bovine serum and other biologicals poses a possible risk of transporting and disseminating viruses into new regions. Limited consideration is given to possible emerging viruses in the existing safety guidelines, and the detection methodology currently in place is biased toward known targets. Furthermore, the pooled nature of most biological products poses an additional limitation for pathogen identification due to the dilution factor. To overcome these challenges, the present study used metagenomics associated with a combination of viral enrichment steps to comprehensively characterize the virome of seven non-irradiated commercial bovine serum batches from Mexico, New Zealand, and the USA.
The significant diversity of the virome identified and the retrieval of several near-full length genomes demonstrated the approach’s robustness. Sequences related to vertebrate viruses from families Parvoviridae, Polyomaviridae, Adenoviridae, Flaviviridae, Picornaviridae, Reoviridae, and Retroviridae were identified. Additionally, CRESS-DNA viruses (Genomoviridae, Circoviridae, and Smacoviridae) and sequences related to the recently described CrAssphage were recognized. Although it is not possible to confirm the absence of sample contamination during processing in the laboratory, the analyses of the unrelated sample (PBMC) suggest that if contamination occurred, it was limited based on the minimal overlap of viral species identified (Supplementary Table S1).
Due to the small size, viruses from the Parvoviridae, and Polyomaviridae families, and the CRESS-DNA viruses are less susceptible to the triple 0.1 µm filtration steps usually employed for serum sterilization. As an additional safety precaution, many institutions employ the use of gamma-irradiated serum. The gamma irradiation effect on virus inactivation [19] is variable depending on the virus structure. Studies have demonstrated that the typical dose of irradiation, 25–40 kGy (2.5–4.0 Mrad), may only lead to a minor viral titer reduction for parvoviruses, polyomaviruses, and circoviruses; in some cases limited to reducing about 1 or 2 logs of the original tissue culture infectious dose (TCID50) [20,21]. Reoviruses and retroviruses also demonstrated resistance to gamma irradiation, with a titer reduction of about 3 to 4 TCID50 logs [22,23]. One limitation of metagenomics is the inability to determine whether the viral sequences were from infectious or inactivated viral particles or free nucleic acid. Despite the ultracentrifugation and nucleases treatment, a significant number of reads matching the host genome were identified and demonstrated that free nucleic acids were sequenced.
Previous studies have identified a reduced number of CRESS-DNA virus species in bovine serum [24,25]. In a virome evaluation of pooled cattle plasma, the DNA synthesis was conducted using a Klenow fragment DNA polymerase, and a total of three groups of CRESS-DNA viruses was identified [24]. The lack of specific enrichment may also have limited the identification of CRESS-DNA families in commercial bovine serum [8]. In the current study, the enrichment of viral DNA was accomplished using the MDA method. A total of 17 CRESS-DNA viral species belonging to the Genomoviridae, Circoviridae, and Smacoviridae families were identified. The enrichment of viral DNA using the MDA method may have favored the circular DNA viral species [26], and therefore, the relative abundance of reads identified may not represent the abundance of virus and/or nucleic acid in the original sample. Whereas the tested serum samples may have a variable number of contaminating viral species in each study, the significantly lower number of viral species identified in other studies may suggest that a combination of viral enrichment is essential before sequencing. Moreover, due to the dilution factor of commercial bovine serum samples, ultracentrifugation of larger amounts of the product prior to nucleic acid extraction may also have favored the increased diversity of viral species in the present study compared to direct nucleic acid extraction [11].
The quality of bovine serum is also assessed by the endotoxin concentration, which may correlate to the initial bacterial load on the collected samples and translate into the overall sanitary conditions employed during collection and processing [27]. While increased endotoxin levels may indicate poor sanitary conditions, no further insights on the origin of bacterial contamination can be retrieved. Among the bacteriophages identified in the present study, CrAssphage was identified in six out of the seven tested batches. CrAssphage has not been reported as a contaminant of bovine serum in previous studies [8,24,25]. This recently discovered bacteriophage infects Bacteroides intestinalis [28,29], and was found to be one of the most abundant human gut phages. In addition to human feces, CrAssphage has been detected in non-human primates and several mammal fecal samples [28,30,31,32,33,34,35,36]. Due to its correlation with fecal material, its presence was suggested as a potential marker for fecal contamination [30,32,37,38,39]. The presence of CrAssphage in the bovine serum lots is likely due to contamination during collection or processing and packing. However, this study cannot determine the precise source of CrAssphage contamination in the tested serum batches.
The ever-growing bovine serum market may lead to product shortage and induce the use of lower quality serum or products with unreliable or untraceable origin. Furthermore, product scarcity associated with the biased testing guidelines and the lack of specific search for viruses from other animal species, including humans, are concerning. Additionally, the current guidelines overlook testing for viruses potentially involved with the disease in humans, like the CRESS DNA viruses. The results demonstrated a high diversity of sequences related to eukaryotic viral families, including CRESS DNA viruses belonging to the Genomoviridae, Circoviridae, and Smacoviridae families. Moreover, representative sequences for 12 bacteriophage families were identified, including related sequences for CrAssphage. Considering the limitations of the current testing, our results indicate that high-throughput sequencing and viral enrichment steps can be a valuable supportive method to attest to the safety and quality of pooled biological products while supporting viral surveillance and early detection of emerging threats.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13122425/s1, Table S1: Summary of sequences retrieved in the control sample (PBMC) that matched with eukaryotic viruses and bacteriophages.

Author Contributions

F.V.B., W.P.P., S.M.F. and C.W.C. conceived and designed the research; W.P.P. and F.V.B. performed the research; W.P.P., F.V.B. and M.N.W. conducted the data analyses; F.V.B., W.P.P. and M.F.M. wrote the original paper; F.V.B., M.F.M., M.N.W., W.P.P., S.M.F., A.R. and C.W.C. thoroughly revised the manuscript; F.V.B. and C.W.C. directed the project; A.R. provided resources; F.V.B. and C.W.C. acquired the funding. All authors have read and agreed to the published version of the manuscript.

Funding

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)-Finance Code 001, Fundação de Amparo à Pesquisa do Rio Grande do Sul (FAPERGS), and Pró-reitoria de Pesquisa da Universidade Federal do Rio Grande do Sul (Propesq/UFRGS).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Animal Care and Use Committee (IACUC) protocol number IACUC-20-14-STW.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. van der Valk, J.; Bieback, K.; Buta, C.; Cochrane, B.; Dirks, W.G.; Fu, J.; Hickman, J.J.; Hohensee, C.; Kolar, R.; Liebsch, M.; et al. Fetal Bovine Serum (FBS): Past–Present–Future. ALTEX 2018, 35, 99–118. [Google Scholar] [CrossRef] [Green Version]
  2. Sykes, J.E.; Rankin, S.C. Chapter 1—Isolation in Cell Culture. In Canine and Feline Infectious Diseases; Sykes, J.E., Ed.; Saunders: Saint Louis, MO, USA, 2014; pp. 2–9. ISBN 978-1-4377-0795-3. [Google Scholar]
  3. Stephens, N.; Di Silvio, L.; Dunsford, I.; Ellis, M.; Glencross, A.; Sexton, A. Bringing cultured meat to market: Technical, socio-political, and regulatory challenges in cellular agriculture. Trends Food Sci. Technol. 2018, 78, 155–166. [Google Scholar] [CrossRef]
  4. Verma, A.; Verma, M.; Singh, A. Animal tissue culture principles and applications. Anim. Biotechnol. 2020, 269–293. [Google Scholar] [CrossRef]
  5. Murray, J.; Versteegen, R. Fetal Bovine—Geographical Origin and International Trade. Bioprocess. J. 2019, 18. [Google Scholar] [CrossRef]
  6. Baylis, S.A.; Miskey, C.; Blümel, J.; Kaiser, M.; Kapusinszky, B.; Delwart, E. Identification of a novel bovine copiparvovirus in pooled fetal bovine serum. Virus Genes 2020, 56, 522–526. [Google Scholar] [CrossRef]
  7. Paim, W.P.; Maggioli, M.F.; Weber, M.N.; Rezabek, G.; Narayanan, S.; Ramachandran, A.; Canal, C.W.; Bauermann, F.V. Virome characterization in serum of healthy show pigs raised in Oklahoma demonstrated great diversity of ssDNA viruses. Virology 2021, 556, 87–95. [Google Scholar] [CrossRef]
  8. Toohey-Kurth, K.; Sibley, S.D.; Goldberg, T.L. Metagenomic assessment of adventitious viruses in commercial bovine sera. Biologicals 2017, 47, 64–68. [Google Scholar] [CrossRef]
  9. Sadeghi, M.; Kapusinszky, B.; Yugo, D.M.; Phan, T.G.; Deng, X.; Kanevsky, I.; Opriessnig, T.; Woolums, A.R.; Hurley, D.J.; Meng, X.J.; et al. Virome of US bovine calf serum. Biologicals 2017, 46, 64–67. [Google Scholar] [CrossRef] [Green Version]
  10. Cheval, J.; Sauvage, V.; Frangeul, L.; Dacheux, L.; Guigon, G.; Dumey, N.; Pariente, K.; Rousseaux, C.; Dorange, F.; Berthet, N.; et al. Evaluation of high-throughput sequencing for identifying known and unknown viruses in biological samples. J. Clin. Microbiol. 2011, 49, 3268–3275. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Gagnieur, L.; Cheval, J.; Gratigny, M.; Hébert, C.; Muth, E.; Dumarest, M.; Eloit, M. Unbiased analysis by high throughput sequencing of the viral diversity in fetal bovine serum and trypsin used in cell culture. Biologicals 2014, 42, 145–152. [Google Scholar] [CrossRef]
  12. Victoria, J.G.; Wang, C.; Jones, M.S.; Jaing, C.; McLoughlin, K.; Gardner, S.; Delwart, E.L. Viral nucleic acids in live-attenuated vaccines: Detection of minority variants and an adventitious virus. J. Virol. 2010, 84, 6033–6040. [Google Scholar] [CrossRef] [Green Version]
  13. Thurber, R.V.; Haynes, M.; Breitbart, M.; Wegley, L.; Rohwer, F. Laboratory procedures to generate viral metagenomes. Nat. Protoc. 2009, 4, 470–483. [Google Scholar] [CrossRef]
  14. Sambrook, J.; Russel, D.W. Molecular Cloning: A Laboratory Manual, 3rd ed.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 2001; 2344p. [Google Scholar]
  15. Morgulis, A.; Coulouris, G.; Raytselis, Y.; Madden, T.L.; Agarwala, R.; Schäffer, A.A. Database indexing for production MegaBLAST searches. Bioinformatics 2008, 24, 1757–1764. [Google Scholar] [CrossRef] [PubMed]
  16. Caesar, L.; Cibulski, S.P.; Canal, C.W.; Blochtein, B.; Sattler, A.; Haag, K.L. The virome of an endangered stingless bee suffering from annual mortality in Southern Brazil. J. Gen. Virol. 2019, 100, 1153–1164. [Google Scholar] [CrossRef] [PubMed]
  17. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  18. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  19. Pecora, A.; Pérez López, J.; Jordán, M.J.; Franco, L.N.; Politzki, R.; Ruiz, V.; Alvarez, I. Analysis of irradiated Argentinean fetal bovine serum for adventitious agents. J. Vet. Diagn. Investig. 2020, 32, 892–897. [Google Scholar] [CrossRef]
  20. Gauvin, G.; Nims, R. Gamma-irradiation of serum for the inactivation of adventitious contaminants. PDA J. Pharm. Sci. Technol. 2010, 64, 432–435. [Google Scholar]
  21. Nims, R.W.; Gauvin, G.; Plavsic, M. Gamma irradiation of animal sera for inactivation of viruses and mollicutes—A review. Biologicals 2011, 39, 370–377. [Google Scholar] [CrossRef]
  22. House, C.; House, J.A.; Yedloutschnig, R.J. Inactivation of viral agents in bovine serum by gamma irradiation. Can. J. Microbiol. 1990, 36, 737–740. [Google Scholar] [CrossRef] [PubMed]
  23. Plavsic, M.; Nims, R.; Wintgens, M.; Versteegen, R. Gamma irradiation of animal serum: Validation of efficacy for pathogen reduction and assessment of impacts on serum performance. Bioprocessing 2016, 15, 1538–8786. [Google Scholar] [CrossRef]
  24. Wang, H.; Li, S.; Mahmood, A.; Yang, S.; Wang, X.; Shen, Q.; Shan, T.; Deng, X.; Li, J.; Hua, X.; et al. Plasma virome of cattle from forest region revealed diverse small circular ssDNA viral genomes. Virol. J. 2018, 15, 11. [Google Scholar] [CrossRef] [Green Version]
  25. Weber, M.N.; Cibulski, S.P.; Silveira, S.; Siqueira, F.M.; Mósena, A.C.S.; da Silva, M.S.; Olegário, J.C.; Varela, A.P.M.; Teixeira, T.F.; Bianchi, M.V.; et al. Evaluation of the serum virome in calves persistently infected with Pestivirus A, presenting or not presenting mucosal disease. Virus Genes 2018, 54, 768–778. [Google Scholar] [CrossRef]
  26. Binga, E.K.; Lasken, R.S.; Neufeld, J.D. Something from (almost) nothing: The impact of multiple displacement amplification on microbial ecology. ISME J. 2008, 2, 233–241. [Google Scholar] [CrossRef] [Green Version]
  27. Kirikae, T.; Tamura, H.; Hashizume, M.; Kirikae, F.; Uemura, Y.; Tanaka, S.; Yokochi, T.; Nakano, M. Endotoxin contamination in fetal bovine serum and its influence on tumor necrosis factor production by macrophage-like cells J774.1 cultured in the presence of the serum. Int. J. Immunopharmacol. 1997, 19, 255–262. [Google Scholar] [CrossRef]
  28. Dutilh, B.E.; Cassman, N.; McNair, K.; Sanchez, S.E.; Silva, G.G.Z.; Boling, L.; Barr, J.J.; Speth, D.R.; Seguritan, V.; Aziz, R.K.; et al. A highly abundant bacteriophage discovered in the unknown sequences of human faecal metagenomes. Nat. Commun. 2014, 5, 4498. [Google Scholar] [CrossRef] [Green Version]
  29. Shkoporov, A.N.; Khokhlova, E.V.; Fitzgerald, C.B.; Stockdale, S.R.; Draper, L.A.; Ross, R.P.; Hill, C. ΦCrAss001 represents the most abundant bacteriophage family in the human gut and infects Bacteroides intestinalis. Nat. Commun. 2018, 9, 4781. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Ahmed, W.; Payyappat, S.; Cassidy, M.; Besley, C.; Power, K. Novel crAssphage marker genes ascertain sewage pollution in a recreational lake receiving urban stormwater runoff. Water Res. 2018, 145, 769–778. [Google Scholar] [CrossRef] [PubMed]
  31. Guerin, E.; Shkoporov, A.; Stockdale, S.R.; Clooney, A.G.; Ryan, F.J.; Sutton, T.D.S.; Draper, L.A.; Gonzalez-Tortuero, E.; Ross, R.P.; Hill, C. Biology and Taxonomy of crAss-like Bacteriophages, the Most Abundant Virus in the Human Gut. Cell Host Microbe 2018, 24, 653–664.e6. [Google Scholar] [CrossRef] [Green Version]
  32. Stachler, E.; Akyon, B.; de Carvalho, N.A.; Ference, C.; Bibby, K. Correlation of crAssphage qPCR Markers with Culturable and Molecular Indicators of Human Fecal Pollution in an Impacted Urban Watershed. Environ. Sci. Technol. 2018, 52, 7505–7512. [Google Scholar] [CrossRef]
  33. Edwards, R.A.; Vega, A.A.; Norman, H.M.; Ohaeri, M.; Levi, K.; Dinsdale, E.A.; Cinek, O.; Aziz, R.K.; McNair, K.; Barr, J.J.; et al. Global phylogeography and ancient evolution of the widespread human gut virus crAssphage. Nat. Microbiol. 2019, 4, 1727–1736. [Google Scholar] [CrossRef] [Green Version]
  34. Brown, B.P.; Chopera, D.; Havyarimana, E.; Wendoh, J.; Jaumdally, S.; Nyangahu, D.D.; Gray, C.M.; Martin, D.P.; Varsani, A.; Jaspan, H.B. crAssphage genomes identified in fecal samples of an adult and infants with evidence of positive genomic selective pressure within tail protein genes. Virus Res. 2021, 292, 198219. [Google Scholar] [CrossRef]
  35. Malla, B.; Ghaju Shrestha, R.; Tandukar, S.; Sherchand, J.B.; Haramoto, E. Performance Evaluation of Human-Specific Viral Markers and Application of Pepper Mild Mottle Virus and CrAssphage to Environmental Water Samples as Fecal Pollution Markers in the Kathmandu Valley, Nepal. Food Environ. Virol. 2019, 11, 274–287. [Google Scholar] [CrossRef]
  36. Li, Y.; Gordon, E.; Shean, R.C.; Idle, A.; Deng, X.; Greninger, A.L.; Delwart, E. CrAssphage and its bacterial host in cat feces. Sci. Rep. 2021, 11, 815. [Google Scholar] [CrossRef]
  37. García-Aljaro, C.; Ballesté, E.; Muniesa, M.; Jofre, J. Determination of crAssphage in water samples and applicability for tracking human faecal pollution. Microb. Biotechnol. 2017, 10, 1775–1780. [Google Scholar] [CrossRef]
  38. Stachler, E.; Kelty, C.; Sivaganesan, M.; Li, X.; Bibby, K.; Shanks, O.C. Quantitative CrAssphage PCR Assays for Human Fecal Pollution Measurement. Environ. Sci. Technol. 2017, 51, 9146–9154. [Google Scholar] [CrossRef]
  39. Park, G.W.; Ng, T.F.F.; Freeland, A.L.; Marconi, V.C.; Boom, J.A.; Staat, M.A.; Montmayeur, A.M.; Browne, H.; Narayanan, J.; Payne, D.C.; et al. CrAssphage as a Novel Tool to Detect Human Fecal Contamination on Environmental Surfaces and Hands. Emerg. Infect. Dis. 2020, 26, 1731–1739. [Google Scholar] [CrossRef]
Viruses 13 02425 g001 550
Figure 1. Genetic characterization of members of the Parvoviridae family was detected in the present study. Representative genomes belonging to the Parvoviridae family were selected from GenBank at the National Center for Biotechnology Information (NCBI). The nucleotide phylogenetic tree was constructed using the Neighbor-Joining (NJ) method based on a p-distance model using the full coding genome. Analyses were conducted with 1000 bootstrap replicates. Complete coding genomes obtained are indicated by the black dots and the sequences were deposited in GenBank under the accession numbers MZ502227, MZ502228, MZ502229, MZ502230, MZ502231, MZ502232, MZ502233, MZ502234, MZ502235, and MZ502236.
Figure 1. Genetic characterization of members of the Parvoviridae family was detected in the present study. Representative genomes belonging to the Parvoviridae family were selected from GenBank at the National Center for Biotechnology Information (NCBI). The nucleotide phylogenetic tree was constructed using the Neighbor-Joining (NJ) method based on a p-distance model using the full coding genome. Analyses were conducted with 1000 bootstrap replicates. Complete coding genomes obtained are indicated by the black dots and the sequences were deposited in GenBank under the accession numbers MZ502227, MZ502228, MZ502229, MZ502230, MZ502231, MZ502232, MZ502233, MZ502234, MZ502235, and MZ502236.
Viruses 13 02425 g001
Viruses 13 02425 g002 550
Figure 2. Phylogenetic tree based on the coding sequence of the polyomavirus large tumor antigen gene. The MEGA6 app suite was used for the phylogeny inference using the Maximum-Likelihood (ML) method based on General Time Reversible (GTR) model. Analyses were conducted with 1000 bootstrap replicates. The polyomavirus reference sequences described in the International Committee on Taxonomy of Viruses were retrieved from the GenBank database. The retrieved bovine polyomavirus 1 (BPyV-1) sequence is represented by the black dot and it was deposited in GenBank under the accession number MZ520135.
Figure 2. Phylogenetic tree based on the coding sequence of the polyomavirus large tumor antigen gene. The MEGA6 app suite was used for the phylogeny inference using the Maximum-Likelihood (ML) method based on General Time Reversible (GTR) model. Analyses were conducted with 1000 bootstrap replicates. The polyomavirus reference sequences described in the International Committee on Taxonomy of Viruses were retrieved from the GenBank database. The retrieved bovine polyomavirus 1 (BPyV-1) sequence is represented by the black dot and it was deposited in GenBank under the accession number MZ520135.
Viruses 13 02425 g002
Viruses 13 02425 g003 550
Figure 3. Taxonomic classification of bacteriophage reads. Bar chart representation of the reads based on BLASTN, E-value < 1 × 10−5 against the GenBank non-redundant database. Sample type: fetal bovine serum (FBS), newborn calf serum (NCS), and bovine calf serum (BCS). Serum origin: Mexico (Mex), New Zealand (NZ), and the United States of America (USA).
Figure 3. Taxonomic classification of bacteriophage reads. Bar chart representation of the reads based on BLASTN, E-value < 1 × 10−5 against the GenBank non-redundant database. Sample type: fetal bovine serum (FBS), newborn calf serum (NCS), and bovine calf serum (BCS). Serum origin: Mexico (Mex), New Zealand (NZ), and the United States of America (USA).
Viruses 13 02425 g003
Table
Table 1. Summary of sequences retrieved in the commercial batches of bovine serum that matched with eukaryotic viruses.
Table 1. Summary of sequences retrieved in the commercial batches of bovine serum that matched with eukaryotic viruses.
Sample Type aNo. of Raw Reads% of Viral/Non-Viral ReadsBest BLASTN Hit bNon-Assembled Reads cAssembled ReadsContig Length Range (nt)Contig Identity dGenome Coverage e
FBS1 MEX772,7540.004/99.996Bovine parvovirus 3184181–31193.7%28.2%
Bovine polyomavirus 12120599.0%4.7%
Moloney murine leukemia virus12120699.5%3.3%
FBS2 USA1,538,3780.117/99.882Bovine parvovirus 2113843107–190497.7%100%
Bovine parvovirus 33121989.0%8.6%
Bovine polyomavirus 166021106–2250100%99.4%
Bovine viral diarrhea virus 152129–18093.0%2.5%
Bovine viral diarrhea virus 262151–22494.7%2.4%
NCS1 USA3,157,4620.453/99.546Bovine bocaparvovirus 287816122–3110100%88.2%
Bosavirus180634101–3021100%100%
Bovine hokovirus 1414120–21795.0%46.9%
Bovine hokovirus 239139598.1%45.6%
Bovine kobuvirus42151–20795.4%4.4%
Bovine parvovirus 18142096.0%12.8%
Bovine parvovirus 211,470297101–2552100%100%
Bovine parvovirus 36713113–62397.7%59.2%
Cyclovirus Equ142103–209100%17.5%
Feces associated gemycircularvirus 142119297.4%9.4%
Gopherus associated genomovirus 12118182.3%10.1%
Mongoose feces-associated gemycircularvirus b2118186.2%8.2%
Po-Circo-like virus 212124693.3%6.2%
Porcine associated porprismacovirus 102115188.7%5.7%
NCS2 NZ2,806,99435.240/64.759Bosavirus184157–411100%26.7%
Bovine parvovirus 24411281–195186.0%25.3%
Bovine parvovirus 3989,1237700100–2471100%100%
Bovine rotavirus A52113–15797.3%7.7%
Starling circovirus2118895.2%9.1%
NCS3 NZ1,082,2062.081/97.918Bovine bocaparvovirus 241245100%5.0%
Bosavirus20,458365100–2050100%100%
Bovine adenovirus 442103–14799.3%0.8%
Bovine parvovirus 22002232101–137096.6%100%
Bovine polyomavirus 1 102141–221100%8.6%
Bovine rotavirus A4129895.3%9.0%
Feces associated gemycircularvirus 14 103190–23899.6%37.7%
Feces associated gemycircularvirus 162114092.1%6.3%
Gopherus associated genomovirus 12120079.0%9.2%
Human rotavirus A4110396.1%4.0%
Mongoose feces-associated gemycircularvirus a2117680.7%8.4%
Mongoose feces-associated gemycircularvirus b 224387–57199.2%88.2%
Sewage-associated gemycircularvirus 32116298.8%7.6%
NCS4 USA1,099,3440.021/99.978Bovine hokovirus 12123388.6%4.6%
Bovine kobuvirus2120394.1%2.5%
Bovine parvovirus 221726149–167998.4%100%
Bovine polyomavirus 62115–23397.6%13.0%
Bovine serum-associated circular virus4137092.0%73.0%
Genomoviridae sp. ctca3672111996.6%5.4%
BCS1 USA1,035,01841.878/58.121Bovine adenovirus 662280–28298.9%2.9%
Bovine hokovirus 1347107–184100%15.3%
Bovine hokovirus 215,00993105–1418100%100%
Bovine parvovirus 2149,9793733100–2888100%100%
Bovine parvovirus 3268,2653569100–2574100%100%
Bo-Circo-like virus CH61149100%3.8%
Circovirus sp. PoCirV VIRES YN02 C23133491.3%31.1%
Po-Circo-like virus 21123132–25295.5%14.4%
Po-Circo-like virus 2213119288.5%4.9%
Po-Circo-like virus 413119193.6%6.5%
Po-Circo-like virus S207126296.5%38.9%
Po-Circo-like virus GX1911148090.2%12.2%
a Sample type: fetal bovine serum (FBS), newborn calf serum (NCS), and bovine calf serum (BCS). Serum origin: Mexico (Mex), New Zealand (NZ), and the United States of America (USA). b Analysis performed in BLASTN (E-value < 1 × 10−5). c Number of read sequences obtained by mapping to the virus RefSeq database (NCBI Reference Sequence Database). d Highest percentage of nucleotide identity based on the reference strain available in the virus RefSeq database (NCBI Reference Sequence Database). e Coverage of the genome based on the length of the reference strain included on the virus RefSeq database (NCBI Reference Sequence Database). The analyses included contig and reads.
Table
Table 2. Characterization of retrieved CrAssphage contig sequences.
Table 2. Characterization of retrieved CrAssphage contig sequences.
Sample Type aNumber of ContigsContig Length RangeBest BLASTN Hit ID of the Longest Contig b/Sample Type/SpeciesContig Identity c
FBS1 Mex0Not applyNot apply Not apply
FBS2 USA180Uncultured crAssphage -NC_024711/Fecal/human86.08%
NCS1 USA6976-625CrAssphage FA1-2_000172F-MK415404/Gut/human99.20%
NCS2 NZ56108–409CrAssphage YS1-2_2437-MK415410/Gut/human98.78%
NCS3 NZ10764–721CrAssphage LMMB-MT006214/Fecal/human94.59%
NCS4 USA8100–241CrAssphage sp. C0531BW4-MW067003/Fecal/human97.93%
BCS1 USA174CrAssphage YS1-2_2437-MK415410/Gut/human82.43%
a Sample type: fetal bovine serum (FBS), newborn calf serum (NCS), and bovine calf serum (BCS). Serum origin: Mexico (Mex), New Zealand (NZ), and the United States of America (USA). b Analysis performed in BLASTN (E-value < 1 × 10−5). c Nucleotide identity based on the reference strain available in the virus RefSeq database (NCBI Reference Sequence Database).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Paim, W.P.; Maggioli, M.F.; Falkenberg, S.M.; Ramachandran, A.; Weber, M.N.; Canal, C.W.; Bauermann, F.V. Virome Characterization in Commercial Bovine Serum Batches—A Potentially Needed Testing Strategy for Biological Products. Viruses 2021, 13, 2425. https://doi.org/10.3390/v13122425

AMA Style

Paim WP, Maggioli MF, Falkenberg SM, Ramachandran A, Weber MN, Canal CW, Bauermann FV. Virome Characterization in Commercial Bovine Serum Batches—A Potentially Needed Testing Strategy for Biological Products. Viruses. 2021; 13(12):2425. https://doi.org/10.3390/v13122425

Chicago/Turabian Style

Paim, Willian P., Mayara F. Maggioli, Shollie M. Falkenberg, Akhilesh Ramachandran, Matheus N. Weber, Cláudio W. Canal, and Fernando V. Bauermann. 2021. "Virome Characterization in Commercial Bovine Serum Batches—A Potentially Needed Testing Strategy for Biological Products" Viruses 13, no. 12: 2425. https://doi.org/10.3390/v13122425

APA Style

Paim, W. P., Maggioli, M. F., Falkenberg, S. M., Ramachandran, A., Weber, M. N., Canal, C. W., & Bauermann, F. V. (2021). Virome Characterization in Commercial Bovine Serum Batches—A Potentially Needed Testing Strategy for Biological Products. Viruses, 13(12), 2425. https://doi.org/10.3390/v13122425

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop