Next Article in Journal
Immune-Escape Hepatitis B Virus Mutations Associated with Viral Reactivation upon Immunosuppression
Previous Article in Journal
Interactions between Autophagy and DNA Viruses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Novel Mammalian orthorubulavirus 5 Discovered as Accidental Cell Culture Contaminant

by
Brandi J. Feehan
1,
Aleksey A. Penin
2,
Alexey N. Mukhin
3,
Deepak Kumar
1,
Anna S. Moskvina
4,
Kizkhalum M. Khametova
3,
Anton G. Yuzhakov
3,
Maria I. Musienko
4,
Alexey D. Zaberezhny
5,
Taras I. Aliper
3,
Douglas Marthaler
1,* and
Konstantin P. Alekseev
3,*
1
Veterinary Diagnostic Laboratory, College of Veterinary Medicine, Kansas State University, Manhattan, KS 66053, USA
2
Institute for Information Transmission Problems of the Russian Academy of Sciences, Moscow 127051, Russia
3
N. F. Gamaleya Federal Research Center for Epidemiology & Microbiology, Moscow 123098, Russia
4
Diagnostic and Prevention Research Institute for Human and Animal Diseases, Moscow 123098, Russia
5
Federal State Budget Scientific Institution “Fedral Scientific Centre VIEV”, Moscow 109428, Russia
*
Authors to whom correspondence should be addressed.
Viruses 2019, 11(9), 777; https://doi.org/10.3390/v11090777
Submission received: 16 July 2019 / Revised: 21 August 2019 / Accepted: 22 August 2019 / Published: 23 August 2019
(This article belongs to the Section Animal Viruses)

Abstract

:
A distinct Russian Mammalian orthorubulavirus 5 (PIV5) was detected in cell culture exhibiting cytopathic effect and hypothesized to be contaminated by a scientist with respiratory symptoms. The identification of the divergent strain indicated a lack of knowledge on the diversity of PIV5 strains and calls for surveillance of global PIV5 strains.

1. Introduction

Mammalian orthorubulavirus 5 (PIV5), formerly named parainfluenza virus 5, resides within the Paramyxoviridae family [1,2,3]. The Paramyxoviridae family contains four subfamilies, Rubulavirinae, Avulavirinae, Metaparamyxovirinae, and Orthoparamyxovirinae. PIV5 is classified with seven other viruses in the Orthorubulavirus genus, one of two genera within Rubulavirinae [1,2]. PIV5 is globally distributed and has been associated with respiratory disease in canine, cattle, swine, and lesser panda [4,5,6,7]. While the pathogenesis of PIV5 in humans is unknown, both PIV5 and contagious Mumps orthorubulavirus belong to the Orthorubulavirus genus [3]. Symptoms of Mumps orthorubulavirus in humans include fever, swollen salivary glands and possibly miscarriages or hearing loss [2,8]. We discovered a divergent Russian PIV5 strain (Moskva) in Vero cells exhibiting cytopathic effect (CPE). The cells were suspected to be accidentally contaminated by an ill laboratory technician, since the samples used for inoculation and the cells were negative for Moskva by PCR. In addition, Moskva contained a small hydrophobic (SH) protein coding sequence (CDS) present in human PIV5 strains.

2. The Study

In November 2015, an outbreak of enteric disease occurred in three-day-old piglets on a farm residing in the Republic of Buryatia, Russia (Supplemental Materials and Methods) [9]. Fecal sample supernatants were passaged on Vero cells in the hope of isolating RVB [9]. Concurrently, the scientist passaging the cell cultures exhibited mild respiratory symptoms including tussis and rhinorrhea. Cytopathic effect (syncytia formation) was observed on the fourth passage of the cells (Figure 1). The nucleic acid from the 15th passage was extracted, and NGS was performed to identify the CPE causative pathogen using an Illumina HiSeq2000 (Illumina, San Diego, CA, USA) [9]. NGS generated 9,280,254 paired reads while 8,949,411 reads resulted after trimming [9]. Kraken identified ~85% (7,602,876) paired reads as Eukaryotic and ~9% (818,499) as viral. Of viral reads, ~65% (534,325) were classified as Paramyxoviridae and ~33% (271,890) were classified as Retroviridae (retrovirus) [10]. De novo assembly of the viral reads generated a PIV5 genome (Moskva, MK593539), which upon BLASTn yielded an 87% nucleotide identity to human PIV5 strain DEN, isolated from the United Kingdom in 1980 (JQ743322). The highest BLASTn hit to a non-PIV5 strain was to a Human orthorubulavirus 2 strain (68%) from Vietnam (MH006623). Primers (forward, CCGGATCACGTGTCCTCAAA; and reverse, ACCAGGAACCGCACTAATGG) were designed to Moskva, and a PCR was developed utilizing RevertAid H Minus First Strand cDNA Synthesis Kit using recommended manufacturer protocols (Thermo Fisher Scientific, Waltham, MA, USA). The original Vero cell stock and fecal samples were negative for Moskva by PCR. Testing of the cell culture from passage four demonstrated PIV5 titers of ~103.5 TCID50/mL [11], supporting the hypothesis of human contamination of the cell culture.
The Moskva sequence was 15,218 nucleotides in length, within an incompletely sequenced 5’ untranslated region due to low coverage (<10 reads). Moskva contained a genome organization of nucleocapsid (N), V, phosphoprotein (P), matrix (M), fusion (F), SH, hemagglutinin-neuraminidase (HN), and large (L) genes [3]. The M protein is the most abundant capsid protein, and the L protein is crucial for RNA synthesis in paramyxoviruses [3]. Both the M and L CDS regions exhibited novel PIV5 truncations of 9 and 11 residues, respectively, which exemplified the diversity of Moskva. In addition, Moskva contained the putative SH gene. The presence of the SH gene varies amongst PIV5 strains, and six of seven human strains contain the SH gene. The SH CDS was present in Moskva, supporting the hypothesis of human contamination of the Vero cell culture [4]. The SH protein has been suggested to limit apoptosis via tumor necrosis factor α (TNFα) and causes a decrease in CPE in vitro [3,12]. The identified CPE and proposed function of the SH protein contradicts our observed CPE in Moskva.
Pairwise nucleotide and amino acid MAFFA alignments were generated in Geneious version 11.5.1 to calculate identities (Geneious, Newark, NJ, USA). Nucleotide percent identities of the complete genome increased the diversity of the PIV5 strains in GenBank from 95.5–100% to 85.4–100% (Table 1). Using the 80% whole genome nucleotide identity classification set by the International Committee on Taxonomy of Viruses, Moskva belongs to the Rubulavirus genus [1,2]. Moskva demonstrated an 84.9–88.8% nucleotide identity to the other PIV5 strains in GenBank regardless of the gene or host species (Table 1), highlighting the diversity of Moskva. Moskva had a 96.9% amino acid identity of RNA dependent RNA polymerase (Large protein) to the other PIV5 strains, illustrating Moskva as a member of the PIV5 species. The divergence of Moskva was reflected as a distinct lineage in a whole genome phylogenetic tree (Figure 2). The 500 bootstrapped maximum-likelihood tree was created in Geneious and edited in FigTree (Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, Scotland). Host and geographic clades were absent in the tree, which was reflected in all individual gene phylogenetic trees (Figure S1). The divergence in genomic identities and phylogenetic trees suggests undiscovered divergent PIV5 strains are circulating in geographically distinct regions.
The PIV5 genome encodes two glycoproteins, F and HN, which are essential for viral attachment, fusion, and entry. With Moskva exhibiting human origin, the F and HN amino acid sequences of Moskva were compared to the F and HN of human strains. The protein models illustrated the location and magnitude of substitutions using SWISS-MODEL (Swiss Institute of Bioinformatics, Lausanne, Switzerland; F, 4jf7.1; HN, 4wsg.1), Phyre2 (Imperial College London, London, UK; F, c2b9bA; HN, c4jf7B), and Chimera (Resource for Biocomputing, Visualization and Informatics, University of California, San Francisco, CA, USA). Consistent with Moskva’s evolutionary divergence, the F and HN proteins of Moskva had more residue substitutions (32 and 31, respectively) relative to consensus sequences of human PIV5 strains (Figure 3A). Of the 32 Moskva F protein substitutions, the majority of the substitutions (n = 25) occurred in unclassified regions of PIV while 5 and 2 substitutions occurred in domains I and III, respectively [13,14]. Domain II (including the hydrophobic loop), the fusion peptide, and heptad repeat region B (HRB) lacked residue substitutions in Moskva [13,15]. The hydrophobic loop plays a crucial role interacting with the HN protein for fusion activation, and the region is conserved in all the strains of Paramyxoviridae [13,15]. In the HN protein, the stalk was more conserved than the head, with 23 residue substitutions in the head versus only two substitutions in the middle of the stalk (Figure 3B). The stalk interacts with the F protein for fusion, and residue substitutions in the stalk near the head of the protein are less detrimental than those substitutions further from the head of the protein [16]. Alterations in the F and HN glycoproteins further represented the divergence of Moskva while supporting regions of conservation.

3. Conclusions

A novel PIV5 strain (Moskva) was reported, and Moskva is hypothesized to have resulted from accidental cross-contamination of Vero cells by a human experiencing respiratory illness. While PIV5 is a global pathogen, our genetic and phylogenetic analysis indicated greater diversity of PIV5 strains than previously reported. The branching of Moskva suggests distinct PIV5 strains are likely circulating globally undetected. Moreover, Moskva had a noteworthy number of residue substitutions in glycoproteins, as well as shorted M and L genes, indicating considerably greater PIV5 diversity than was previously known. Further research is required to understand the implications of PIV5 global presence, the existence of the SH protein, genomic substitutions, and protein truncations (M and F), especially as these topics pertain to human health.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/9/777/s1, Figure SA–H: Individual gene phylogenetic trees of PIV5 strains.

Funding

Funding of the HiSeq sequencing was supported by Russian Academy of Sciences (Projects 0053-2019-0005). Brandi Feehan was partially supported by the Veterinary Diagnostic Laboratory at Kansas State University.

Acknowledgments

The authors would like to thank Doctor Anastasia Vlasova at The Ohio State University for her critical insight.

About the Author

Brandi J. Feehan is a master’s student within the Veterinary Biomedical Sciences program under Douglas Marthaler at Kansas State University. Her primary research objective is to employ bioinformatics to elucidate viromes, resolve quasispecies, and determine relationships amongst different viral populations infecting animals.

Ethics Statement

The authors confirm that the ethical policies of the journal, as noted on the journal’s author guidelines page, have been adhered to and the appropriate ethical review committee approval has been received. The guidelines of the local Ethical and Animal Welfare Committee of the Federal State Budget Scientific Institution “Federal Scientific Centre VIEV”, (Moscow, Russia) were followed.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. International Committee on Taxonomy of Viruses (ICTV). Available online: https://talk.ictvonline.org//taxonomy/p/taxonomy-history?taxnode_id=201856468 (accessed on 30 April 2019).
  2. N.d.: Virus Pathogen Database and Analysis Resource (ViPR) Paramyxoviridae. Available online: https://www.viprbrc.org/brc/aboutPathogen.spg?decorator=paramyxo (accessed on 20 August 2019).
  3. Knipe, D.M.; Howley, P.M. Fields’ Virology; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2007. [Google Scholar]
  4. Lee, Y.N.; Lee, C. Complete genome sequence of a novel porcine parainfluenza virus 5 isolate in Korea. Arch. Virol. 2013, 158, 1765–1772. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, Y.; Li, N.; Zhang, S.; Zhang, F.; Lian, H.; Hu, R. Parainfluenza Virus 5 as Possible Cause of Severe Respiratory Disease in Calves, China. Emerg. Infect. Dis. 2015, 21, 2242–2244. [Google Scholar] [CrossRef] [PubMed]
  6. Oem, J.-K.; Kim, S.-H.; Kim, Y.-H.; Lee, M.-H.; Lee, K.-K. Molecular characteristics of canine parainfluenza viruses type 5 (CPIV-5) isolated in Korea. Can. J. Vet. Res. 2015, 79, 64–67. [Google Scholar] [PubMed]
  7. Zhai, J.-Q.; Zhai, S.-L.; Lin, T.; Liu, J.-K.; Wang, H.-X.; Li, B.; Zhang, H.; Zou, S.-Z.; Zhou, X.; Wu, M.-F.; et al. First complete genome sequence of parainfluenza virus 5 isolated from lesser panda. Arch. Virol. 2017, 162, 1413–1418. [Google Scholar] [CrossRef] [PubMed]
  8. Mayo Clinic, n.d.: Mumps Symptoms and Causes. Available online: http://www.mayoclinic.org/diseases-conditions/mumps/symptoms-causes/syc-20375361 (accessed on 1 December 2018).
  9. Alekseev, K.P.; Penin, A.A.; Mukhin, A.N.; Khametova, K.M.; Grebennikova, T.V.; Yuzhakov, A.G.; Moskvina, A.S.; Musienko, M.I.; Raev, S.A.; Mishin, A.M.; et al. Genome Characterization of a Pathogenic Porcine Rotavirus B Strain Identified in Buryat Republic, Russia in 2015. Pathogens 2018, 7, 46. [Google Scholar] [CrossRef] [PubMed]
  10. Knutson, T.P.; Velayudhan, B.T.; Marthaler, D.G. A porcine enterovirus G associated with enteric disease contains a novel papain-like cysteine protease. J. Gen. Virol. 2017, 98, 1305–1310. [Google Scholar] [CrossRef] [PubMed]
  11. Reed, L.J.; Muench, H. A simple method of estimating fifty per cent endpoints. Am. J. Epidemiol. 1938, 27, 493–497. [Google Scholar] [CrossRef]
  12. Lin, Y.; Bright, A.C.; Rothermel, T.A.; He, B. Induction of Apoptosis by Paramyxovirus Simian Virus 5 Lacking a Small Hydrophobic Gene. J. Virol. 2003, 77, 3371–3383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Bose, S.; Heath, C.M.; Shah, P.A.; Alayyoubi, M.; Jardetzky, T.S.; Lamb, R.A. Mutations in the Parainfluenza Virus 5 Fusion Protein Reveal Domains Important for Fusion Triggering and Metastability. J. Virol. 2013, 87, 13520–13531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Poor, T.A.; Song, A.S.; Welch, B.D.; Kors, C.A.; Jardetzky, T.S.; Lamb, R.A. On the Stability of Parainfluenza Virus 5 F Proteins. J. Virol. 2015, 89, 3438–3441. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Waning, D.L.; Russell, C.J.; Jardetzky, T.S.; Lamb, R.A. Activation of a paramyxovirus fusion protein is modulated by inside-out signaling from the cytoplasmic tail. Proc. Natl. Acad. Sci. USA 2004, 101, 9217–9222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Bose, S.; Welch, B.D.; Kors, C.A.; Yuan, P.; Jardetzky, T.S.; Lamb, R.A. Structure and Mutagenesis of the Parainfluenza Virus 5 Hemagglutinin-Neuraminidase Stalk Domain Reveals a Four-Helix Bundle and the Role of the Stalk in Fusion Promotion. J. Virol. 2011, 85, 12855–12866. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Cell cultures illustrating cytopathic effect (CPE) in the sequenced culture with Moskva (right) relative to the control (left).
Figure 1. Cell cultures illustrating cytopathic effect (CPE) in the sequenced culture with Moskva (right) relative to the control (left).
Viruses 11 00777 g001
Figure 2. Whole genome phylogenetic tree of Mammalian orthorubulavirus 5 (PIV5) strains. Bootstrap values greater than 70 are shown at major nodes. Strains are colored according to host with black representing human strains; red, canine; orange, swine; green, simian; teal, pangolin; blue, tiger; pink, panda; and purple, bovine. Moskva is in bold.
Figure 2. Whole genome phylogenetic tree of Mammalian orthorubulavirus 5 (PIV5) strains. Bootstrap values greater than 70 are shown at major nodes. Strains are colored according to host with black representing human strains; red, canine; orange, swine; green, simian; teal, pangolin; blue, tiger; pink, panda; and purple, bovine. Moskva is in bold.
Viruses 11 00777 g002
Figure 3. Protein models illustrating the residue substitutions in Moskva. Monomer and trimer models in uncleaved Fusion (A) and Hemagglutinin-neuraminidase (HN) (B) proteins. Regions and substitutions are colored according to the legend.
Figure 3. Protein models illustrating the residue substitutions in Moskva. Monomer and trimer models in uncleaved Fusion (A) and Hemagglutinin-neuraminidase (HN) (B) proteins. Regions and substitutions are colored according to the legend.
Viruses 11 00777 g003aViruses 11 00777 g003b
Table 1. Comparison of the nucleotide (nt) and amino acid (aa) lengths and identity for PIV5 strains.
Table 1. Comparison of the nucleotide (nt) and amino acid (aa) lengths and identity for PIV5 strains.
Gene/GenomeLengthExcluding Moskva n = 26Including Moskva n = 27Human n = 7Simian n = 2Swine n = 4Bovine, Panda, Pangolin, Tiger n = 6
ntaantaantaantaantaantaantaa
Nucleocapsid (N)153050995.8–10097.7–10086.3–10093.5–10086.8–87.594.195.187.894.587.3–87.793.9–94.787.5–87.793.7–94.5
V (V)66922295.5–10094.6–10087.3–10094.6–10087.788.894.6–96.988.897.387.7–88.395.5–96.987.4–88.295.5–96.9
Phosphoprotein (P)117739195.8–10095.9–10087.1–10094.1–10088.188.694.4–95.988.696.287.8–88.394.6–95.787.9–88.294.9–95.7
Matrix (M)1134 *377 *95.3–10096.0–10086.2–10093.7–10087.087.795.096.687.4–87.595.5–95.887.1–87.495.5–95.886.6–87.494.4–95.8
Fusion (F)1656 551 94.8–10095.9–10083.3–10091.5–10084.985.893.394.685.293.7–93.885.0–85.193.3–93.584.9–85.193.5–93.7
Hemagglutinin-neuraminidase (HN)169856595.6–10095.1–10084.8–10091.3–10085.5–85.892.894.585.6–85.792.685.6–85.891.9–92.684.9–85.892.2–92.8
Large (L)6678 2256 96.9–10098.1–10087.4–10096.9–10087.988.297.397.888.297.587.9–88.096.9–97.387.888.297.2–97.6
Entire genome15,246§ 95.5–100 85.4–100 86.286.7 86.686.7 86.2–86.4 86.4–86.5
* Moskva: 1107 nt and 368 aa. Simian strains: 1589 nt and 529 aa; Human MEL strain: 1605 nt and 534 aa. Moskva: 6735 nt and 2245 aa. § Moskva: 15,218 nt. Boldface indicates the highest percent identity to Moskva in the corresponding gene/genome.

Share and Cite

MDPI and ACS Style

Feehan, B.J.; Penin, A.A.; Mukhin, A.N.; Kumar, D.; Moskvina, A.S.; Khametova, K.M.; Yuzhakov, A.G.; Musienko, M.I.; Zaberezhny, A.D.; Aliper, T.I.; et al. Novel Mammalian orthorubulavirus 5 Discovered as Accidental Cell Culture Contaminant. Viruses 2019, 11, 777. https://doi.org/10.3390/v11090777

AMA Style

Feehan BJ, Penin AA, Mukhin AN, Kumar D, Moskvina AS, Khametova KM, Yuzhakov AG, Musienko MI, Zaberezhny AD, Aliper TI, et al. Novel Mammalian orthorubulavirus 5 Discovered as Accidental Cell Culture Contaminant. Viruses. 2019; 11(9):777. https://doi.org/10.3390/v11090777

Chicago/Turabian Style

Feehan, Brandi J., Aleksey A. Penin, Alexey N. Mukhin, Deepak Kumar, Anna S. Moskvina, Kizkhalum M. Khametova, Anton G. Yuzhakov, Maria I. Musienko, Alexey D. Zaberezhny, Taras I. Aliper, and et al. 2019. "Novel Mammalian orthorubulavirus 5 Discovered as Accidental Cell Culture Contaminant" Viruses 11, no. 9: 777. https://doi.org/10.3390/v11090777

APA Style

Feehan, B. J., Penin, A. A., Mukhin, A. N., Kumar, D., Moskvina, A. S., Khametova, K. M., Yuzhakov, A. G., Musienko, M. I., Zaberezhny, A. D., Aliper, T. I., Marthaler, D., & Alekseev, K. P. (2019). Novel Mammalian orthorubulavirus 5 Discovered as Accidental Cell Culture Contaminant. Viruses, 11(9), 777. https://doi.org/10.3390/v11090777

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