Next Article in Journal
Dendritic Cell Targeting of Bovine Viral Diarrhea Virus E2 Protein Expressed by Lactobacillus casei Effectively Induces Antigen-Specific Immune Responses via Oral Vaccination
Previous Article in Journal
Susceptibility of Chickens to Porcine Deltacoronavirus Infection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Picobirnaviruses in Respiratory and Alimentary Tracts of Cattle and Monkeys with Large Intra- and Inter-Host Diversity

1
State Key Laboratory of Emerging Infectious Diseases, The University of Hong Kong, Hong Kong, China
2
Department of Microbiology, The University of Hong Kong, Hong Kong, China
3
Carol Yu Centre for Infection, The University of Hong Kong, Hong Kong, China
4
Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The University of Hong Kong, Hong Kong, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2019, 11(6), 574; https://doi.org/10.3390/v11060574
Submission received: 21 May 2019 / Revised: 17 June 2019 / Accepted: 21 June 2019 / Published: 23 June 2019
(This article belongs to the Section Animal Viruses)

Abstract

:
Picobirnaviruses (PBVs) are mostly found in animal alimentary samples. In this study, among 576 respiratory specimens from 476 mammals and 100 chickens, genogroup I PBVs were detected in three cattle and three monkeys, and a genogroup II PBV-positive sample was collected from one cattle specimen. More than one PBV sequence type was observed in two and one genogroup I PBV-positive samples from cattle and monkeys, respectively. Twenty-four complete/near-complete segments 2 (nine from respiratory and 15 from alimentary samples) from the cattle and monkey genogroup I PBVs and one complete segment 2 from the cattle genogroup II PBV were sequenced. Similar to other studies, the cattle PBVs also showed a high diversity. In contrast, the monkey PBVs observed in this study were clustered into three distinct clades. Within each clade, all the sequences showed >99% amino acid identities. This unique phenomenon is probably due to the fact that monkeys in our locality reside in separated troops with minimal inter-troop contact.

1. Introduction

Picobirnaviruses (PBVs) are small non-enveloped bi-segmented double-stranded RNA viruses. Since their first discovery in alimentary samples of humans and rats in 1988 [1,2], PBVs have been reported in other terrestrial and marine mammals, birds, invertebrates and environmental water samples [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]. Recently, we have described the first discovery and a high diversity of PBVs in California sea lions from an oceanarium in Hong Kong and dromedary camels from the Middle East [18,20,21]. We also found that multiple strains of PBV were present in the majority of PBV-positive samples in different kinds of animals [18].
The genome of a PBV consists of two segments: segment 1 and segment 2. Segment 1 contains the capsid gene and another open reading frame which encodes a putative protein of an unknown function, whereas segment 2 contains the RNA-dependent RNA polymerase (RdRp) gene [22]. By conducting sequence and phylogenetic analyses of RdRp genes, PBVs were classified into genogroups I and II [22,23]. Recently, it was reported that a third genogroup of PBV, genogroup III, was obtained from invertebrate hosts and differs largely from the PBVs of genogroups I and II [15,24]. In our recent study in which serial samples from the same animal were collected, we observed that PBV probably evolves through mechanisms similar to other segmented RNA viruses, such as rotaviruses (double-stranded RNA virus) [18]. Rotaviruses evolve through reassortment of their genomes, resulting in enormous diversity [25,26]. This also contributes to the emergence of new rotavirus strains, resulting in disease outbreaks [25,26]. However, in a recent study on PBV detected in a marmot of the Qinghai-Tibet plateau, it was found that segmentation-associated motifs were identified in both bi-segmented and unsegmented PBVs. Unsegmented PBVs may undergo genome segmentation at the motif region and convert to segmented PBVs, which represents a unique segmentation mechanism [9].
Traditionally, PBVs were found in alimentary samples of animals. However, in two recent studies, PBVs were amplified from >10% and 0.97% of respiratory samples collected from pigs and a human, respectively [27,28]. Such discoveries have improved our understandings on the tissue tropism of PBVs in the corresponding hosts. We hypothesize that PBVs may also reside in the respiratory tracts of animals other than pigs and humans. To test this hypothesis, we carried out a molecular epidemiology study of PBVs in respiratory samples collected from mammals and poultry. Genome segment 2 sequencing, comparative genomics, and phylogenetic analysis of the PBVs detected in the respiratory and alimentary samples of cattle and monkeys were also performed.

2. Materials and Methods

2.1. Ethics Statement

The study was approved by the Committee on the Use of Live Animals for Teaching and Research, The University of Hong Kong (CULATR Ref. No.: 2284-10 and 3330-14; Date of approval: 23 March 2011 and 17 April 2014).

2.2. Mammal and Poultry Surveillance and Sample Collection

All specimens of bats, monkeys, dogs, and cats were collected with the assistance of the Department of Agriculture, Fisheries and Conservation (AFCD), Hong Kong Special Administrative Region (HKSAR); those of cattle and chickens were collected with the assistance of the Department of Food, Environmental and Hygiene (FEHD), HKSAR, from various locations in HKSAR; and those of horses were collected with the assistance of the Hong Kong Jockey Club. All respiratory specimens (nasopharyngeal swabs in cattle, dogs, cats, and horses; throat swabs in monkeys; mouth swabs in bats; and tracheal swabs in chickens) were collected using procedures described previously [29]. A total of 576 samples collected over an 83-month period (March 2006 to January 2013) in Hong Kong from 157 bats, 52 monkeys, 50 dogs, 50 cats, 116 horses, 51 cattle, and 100 chickens have been tested (Table 1).

2.3. RNA Extraction

Viral RNA was extracted from the respiratory samples using the EZ1 Virus Mini Kit v2.0 according to the manufacturer’s instructions (Qiagen, Hilden, Germany). The RNA was eluted in 60 μL of AVE buffer (Qiagen, Hilden, Germany) and used as the template for RT-PCR.

2.4. RT-PCR for PBVs and DNA Sequencing

PBV screening was performed by PCR amplification of a 205-bp (nucleotide position 967–1171) and a 207-bp fragment (nucleotide position 786–992) of the RdRp gene of genogroup I and genogroup II PBVs, respectively, using conserved primers (genogroup I: 5′-CAAARTTYGACCARCACTT-3′ and 5′-TCRTCDGCRTTGGTACCACC-3′; genogroup II: 5′-WTGGATGTTTCCGATGTC-3′ and 5′-TGYGCATCCATYTTMGTGGTGTCTC-3′) designed based on multiple alignments of the nucleotide sequences of available RdRp genes of PBVs [30]. Reverse transcription was performed using the SuperScript III kit (Invitrogen, San Diego, CA, USA) and the reaction mixture (10 μL) contained RNA, first-strand buffer (50 mM Tris-HCl pH 8.3, 75 mM KCl, 3 mM MgCl2), 5 mM dithiothreitol, 50 ng random hexamers, 500 μM of each deoxyribonucleotide triphosphate and 100 U Superscript III reverse transcriptase. The mixtures were incubated at 25 °C for 5 min, followed by 50 °C for 60 min and 70 °C for 15 min. The PCR mixture (25 μL) contained cDNA, PCR buffer (10 mM Tris-HCl pH 8.3, 50 mM KCl, 2 mM MgCl2), 200 μM of each dNTP, 10 μM of each primer, and 1.0 U Taq polymerase (Applied Biosystems, Foster City, CA, USA). The mixtures were amplified in 60 cycles of 94 °C for 1 min, 50 °C for 1 min, and 72 °C for 1 min, and a final extension at 72 °C for 10 min in an automated thermal cycler (Applied Biosystems, Foster City, CA, USA). Standard precautions were taken to avoid PCR contamination and no false-positive was observed in the negative controls.
All PCR products were gel-purified using the QIAquick gel extraction kit according to the manufacturer’s instructions (Qiagen, Hilden, Germany). Both strands of the PCR products were sequenced twice with a Prism 3700 DNA analyzer (Applied Biosystems, Foster City, CA, USA), using the two PCR primers. When multiple nucleotide peaks were observed in the sequencing results, which suggested more than one type of PBV in the sample, the corresponding purified PCR product was cloned into the pCR-II-TOPO TA vector (Invitrogen, San Diego, CA, USA) according to the manufacturer’s instructions. Both strands of 10 clones for each sample were sequenced, using primers 5′-TAATACGACTCACTATAGGG-3′ and 5′-CGGCTCGTATGTTGTGTGGA-3′, which bind to the nucleotide position 160–179 and 406–425 of the vector, respectively. The sequences of the PCR products or the clones were compared with known sequences of the RdRp of PBVs in the GenBank database.

2.5. Complete/Near-Complete Segment 2 Sequencing of Genogroups I and II PBVs

Five and four complete segment 2 samples of cattle and monkey genogroup I PBVs, and one complete segment 2 of cattle a genogroup II PBV detected in the respiratory specimens were amplified and sequenced respectively using published strategies for double-stranded RNA viruses [31]. In addition, six and nine complete/near-complete segment 2 samples of the cattle and monkey PBVs we previously detected in the alimentary specimens were also amplified and sequenced using the same strategy [31]. RNA extracted from the original respiratory and alimentary specimens, using the EZ1 virus mini kit (Qiagen, Hilden, Germany), was used as the template for RT-PCR and sequencing. An adaptor primer, with a 3′ NH2 blocking group, was ligated to the 3′ termini of the viral RNA and subjected to reverse transcription using a complementary primer. After RNA hydrolysis, reannealing, and end-filling, single-primer amplification of viral genomic segments was performed using a complementary primer and genome specific primers using a strategy described previously [32]. The 5′ and 3′ ends of the viral genomes were confirmed by the rapid amplification of cDNA ends using the 5′/3′ RACE kit (Clontech, Mountain View, CA, USA). The PCR products were gel purified and cloned into the pCR-Blunt II-TOPO vector by employing the Zero Blunt TOPO PCR cloning kit (Invitrogen, San Diego, CA, USA). Clones were sequenced using an ABI Prism 3700 DNA analyzer (Applied Biosystems, Foster City, CA, USA). Sequences were assembled and manually edited to produce the final sequences of the segment 2 genome.

2.6. Genome Segment 2 Analysis

The nucleotide sequences of genome segment 2 and the deduced amino acid sequences of the open reading frames (ORFs) were compared to those of other PBVs. Phylogenetic analysis was performed by the neighbor-joining method, using the Jukes–Cantor substitution model with gamma distributed rate variation and bootstrap values calculated from 1000 trees [30]. A clade is defined as two or more PBVs sharing >99% amino acid identities in their complete or near-complete RdRp sequences. Secondary structure prediction in the 5′UTR was performed using RNAfold [33].

2.7. Estimation of Synonymous and Non-Synonymous Substitution Rates

The number of synonymous substitutions per synonymous site, Ks, and the number of non-synonymous substitutions per non-synonymous site, Ka, for the RdRp between monkey PBVs belonging to clade M-R1 were calculated using the Nei-Gojobori method (Jukes-Cantor) in MEGA6 [30]. The Ks and Ka for the other two clades of monkey PBVs were not calculated because the number of sequences was small.

2.8. Nucleotide Sequence Accession Numbers

The genome segment 2 sequences of all cattle and monkey PBVs have been lodged within the GenBank sequence database under accession no. KY120170-KY120194.

3. Results

3.1. Detection of PBVs in Respiratory Specimens of Animals

A total of 576 respiratory specimens from 476 mammals and 100 chickens were obtained (Table 1). RT-PCR for a 205-bp fragment in the RdRp of a genogroup I PBV was positive in respiratory samples collected from three cattle and three monkeys, respectively (Table 1). Sequencing and phylogenetic analysis showed that the three positive samples from cattle (C345N, C346N, and C361N) contained three (C345N-1, C345N-2, and C345N-3), two (C346N-1 and C346N-2), and one (C361N) sequence types, respectively (Figure 1A). There were 63 (31.3%) to 73 (35.8%) nucleotide differences among C345N-1, C345N-2, and C345N-3 and 62 (30.5%) nucleotide differences between C346N-1 and C346N-2. For the three positive samples from monkeys (11T, 21T, and 35T), two sequence types (11T-1 and 11T-2) were observed in one sample (11T) and one sequence type was observed in the other two samples (21T and 35T) (Figure 1A). There were 64 (31.5%) nucleotide differences between 11T-1 and 11T-2. The PBV sequences from cattle and monkeys were widely distributed in the whole phylogenetic tree with a high diversity (Figure 1A). No genogroup I PBV was amplified from respiratory samples collected from 50 dogs, 50 cats, 157 bats, 116 horses, and 100 chickens.
RT-PCR for a 207-bp fragment in the RdRp of a genogroup II PBV was positive in one respiratory sample collected from cattle (Table 1). Sequencing and phylogenetic analysis showed that the sequence possessed 46–89 (22.8–50.6%) nucleotide differences to other genogroup II PBV sequences, being most closely related to macaque PBV (MG010908.1) (Figure 1B). No genogroup II PBV was amplified from respiratory samples collected from the 52 monkeys, 50 dogs, 50 cats, 157 bats, 116 horses, and 100 chickens. No co-infection with both genogroups of PBVs in a single sample was observed.

3.2. Complete/Near-Complete Segment 2 Sequence Analysis of Cattle and Monkey PBVs

Twenty complete and four near-complete segment 2 samples (nine from respiratory and 15 from alimentary samples) from the cattle and monkey genogroup I PBVs were successfully sequenced and assembled (Table 2, Figure 2). For the three positive alimentary samples from cattle (C343R, C369R, and C374R), each of them contained two sequence types (C343R-1 and C343R-2, C369R-1 and C369R-2, and C374R-1 and C374R-2, respectively). For the six positive alimentary samples from monkeys (1R, 2R, 12R, 13R, 14R, and 15R), two sequence types (1R-1 and 1R-2, 13R-1 and 13R-2, and 15R-1 and 15R-2, respectively) were observed in three samples (1R, 13R, and 15R) and one sequence type was observed in the other three samples (2R, 12R, and 14R) [18]. Sequencing of the other eight segment 2 samples was not successful because of low viral loads and/or failure to amplify part of the segment. These 24 sequenced segment 2 samples ranged from 1613 to 1736 bases in length, with overall G+C contents of 41.1 to 47.7%. The 5′ non-coding regions (27 to 74 bases) are AU-rich (G+C contents 14.8 to 35.4%) with five conserved bases, GUAAA, located at the 5′ end. A highly stable stem loop structure observed in other known PBVs was observed in 17 of 20 complete segment 2 sequences as a result of the pairing of 5′-GUAAA-3′ and 5′-UUUAC-3′ sequences in the 5′ non-coding region [34]. The 3′ non-coding regions (15 to 50 bases) had G+C contents that ranged from 40.0 to 59.5%, ending with four conserved bases (CUGC) in most of the genomes. Excluding the one detected in monkey PBV 1R-1 with incomplete RdRp, all the other 23 sequenced segment 2 samples possessed one long ORF (1563 to 1686 bp) encoding the putative RdRp of 520 to 561 amino acids. These putative RdRp sequences shared 43.0 to 80.0% amino acid identities with those of other genogroup I PBV strains, being most closely related to dromedary picobirnavirus (GenBank number BBC20590.1), human PBV HuPBV-E-CDC16 (Genbank number KJ663816), human PBV GPBV6C1 (Genbank number AB517731), human PBV VS10 (Genbank number GU968924), and fox PBV F5-1 (Genbank number KC692366). All segment 2 samples, except the one detected in cattle (C374R-1) with the first motif being D-E-D instead of D-T/S-D, possessed three conserved motifs (D-T/S-D, SG-T, and GDD) commonly found in the RdRp sequences of other dsRNA viruses [6]. Conserved cysteine and proline residues present in other genogroup I PBVs were also observed in all 24 segment 2 sequences [6].
One complete segment 2 from the cattle genogroup II PBV was sequenced and assembled (Table 2, Figure 2). This segment 2 was 1622 bases in length with a G+C content of 40.4%. The 5′ non-coding region (54 bases) was AU-rich (G+C content 16.7%) with the five conserved bases, GUAAA, located at the 5′ end. Similar to genogroup I PBVs, a highly stable stem loop structure was also observed in the 5′ non-coding region as a result of the pairing of 5′-GUAAA-3′ and 5′-UUUAC-3′ sequences [34]. The 3′ non-coding region (23 bases) had a G+C content of 52.2%, ending with CUCC. This segment 2 possessed one long ORF (1545 bp) encoding the putative RdRp of 514 amino acids. This putative RdRp shared 40.7 to 78.0% amino acid identities with those of other genogroup II PBV strains, being most closely related to macaque PBV (Genbank number AVD54059.1). The three conserved motifs (D-T/S-D, SG-T, and GDD) commonly found in the RdRp sequences of other dsRNA viruses were also observed [34].

3.3. Phylogenetic Analysis of Cattle and Monkey PBVs’ Complete/Near-Complete Segment 2 Sequences

In the phylogenetic tree constructed using the RdRp gene, the 13 segment 2 sequences of the monkey PBVs were clustered into three (M-R1 to M-R3) distinct clades (Figure 2). Within each clade, all the sequences showed >99% amino acid identities, suggesting that they belong to the same clone. The alimentary samples of three monkeys (1R, 13R, and 15R) possessed sequences that belonged to two different clades (M-R2 and M-R3 for sample 1R; M-R1 and M-R3 for samples 13R and 15R), and the throat swab of one monkey (11T) possessed sequences that belonged to two different clades (M-R1 and M-R2). As for the sequences from the cattle samples, except for two sequences (C346N-1 and C374R-2), no distinct clades were observed, suggesting that there was no dominant clone in these alimentary and respiratory samples (Figure 2).

3.4. Estimation of Substitution Rates

The Ka, Ks, and Ka/Ks of the RdRp in monkey PBVs belonging to clade M-R1 are 0.002, 0.049 and 0.041, respectively.

4. Discussion

We discovered PBVs from respiratory samples of cattle and monkeys. The original discovery of PBV from alimentary samples raised the question of whether PBV is associated with gastroenteritis [1,2]. The recent findings on the presence of genogroup I and II PBV in respiratory tract swab specimens of pigs, genogroup I PBV in bronchoalveolar lavage of humans, and non-genogroup I/II PBV in nasal swabs of young dairy cattle with cattle respiratory disease further prompted us to rethink whether these viruses could be respiratory tract pathogens [27,28,35]. Moreover, the high abundance of PBV in respiratory samples of pigs is also in line with its high abundance in their alimentary samples [27,36]. In this study, we further provided evidence that PBVs can be detected in respiratory samples of cattle and monkeys in Hong Kong. Although PBV was not detected in respiratory samples of dogs, cats, bats, and chickens in the present study, this may partly be due to the relatively low abundance of PBV in these animals, as our recent study on mammalian PBVs and our unpublished data on poultries also did not reveal a high abundance of PBV in their alimentary samples. Notably, similar to results in previous reports [37,38], genogroup I PBV sequences are more commonly found than genogroup II PBV sequences in the present study. This may be due to a genuine higher abundance of genogroup I PBV or the use of primers better designed for the amplification of genogroup I than genogroup II PBV as a result of the availability of many more genogroup I PBV sequences in GenBank. In the present and our previous studies [18,20], we found that the topologies of the phylogenetic trees constructed based on the partial RdRp sequences (around 200-bp) and complete RdRp sequences of PBVs were completely concordant (i.e., correct classification of PBVs into genogroup I and genogroup II). This implies that sequencing partial RdRp PCR products is already sufficient for classifying genogroup I and genogroup II PBVs.
Multiple strains of PBV were present in the respiratory samples of monkeys and cattle. Up until early 2019, although PBVs have been found in different animal species [2,3,7,8,9,10,13,17,18,19,20,21,39,40,41,42,43,44,45], more than one strain of PBV has only been described in alimentary samples collected from humans [5,38], pigs [36,39], chickens [14], monkeys [17], and buffalos [43]. In our recent study, multiple strains of PBV were observed in the same alimentary samples of cattle, monkeys, horses, pigs, and California sea lions [18]. In the present study, among the three monkeys and four cattle positive for PBVs in their respiratory samples, more than one strain of PBV was present in two cattle and one monkey, as confirmed by RdRp gene sequencing. This phenomenon of multiple strains of PBV in one respiratory sample has been reported in pigs, from which both genogroups I and II PBVs were observed in the same sample [27].
The cattle PBVs observed from the respiratory and alimentary samples showed a high diversity. In our recent study, it was found that the PBV sequences from different kinds of mammals showed a high diversity and phylogenetic analysis revealed that they were widely distributed over the phylogenetic tree [18]. Although genogroup I PBVs have been found in different kinds of animals [2,3,7,8,9,10,13,17,18,19,20,21,39,40,41,42,43,44,45], genogroup II PBVs have only been reported in humans [5,46], monkeys [47], pigs [36], a cattle calf [44], and a cat [19]. In this study, both genogroups I and II PBVs were present in cattle, with the genogroup I cattle PBVs showing a high diversity. Except for the two cattle (C346N-1 and C374R-2) of which the RdRp genes of their PBVs showed 93.6% amino acid identity, all the other cattle PBVs had <80% amino acid identities with each other. In our recent study, we proposed that PBVs evolved through genome reassortment and high mutation rates, which may partly account for the high diversities observed in PBVs from different kinds of animals [18]. Furthermore, a recent study reported that PBVs may infect prokaryotes [47], suggesting that the interplay of PBVs and prokaryotes may also account for the evolution and genetic diversity of PBVs.
In contrast to the high diversities observed in PBVs from cattle, horses, pigs, and California sea lions from the present and our previous study, three distinct PBV clones were found in the monkeys in this study. In Hong Kong, most monkeys in the wild are hybrids of the rhesus macaque and long-tailed macaque. The total number of macaques in Hong Kong is about 1800 to 2000. They live as 30 geographically separated troops. The respiratory and alimentary samples collected in this study were from monkeys that belonged to two separated troops 1.6 km apart. For clades M-R1 and M-R2, it can be observed that the same PBV clone was present in both the respiratory and alimentary samples (Figure 2). This implies that the same PBV can probably infect both the respiratory and gastrointestinal tracts of monkeys. This unique phenomenon observed in these PBVs in monkeys in this study has also given us the opportunity to perform molecular evolution studies for protein-coding genes in monkey PBVs of the M-R1 clade, which showed that the Ka/Ks ratio was <0.05, implying that this virus was stably evolving in monkeys. Further studies should be performed to examine if there is any association between PBV and diarrhea in monkeys, as well as their diversity and evolutionary mechanism.

Author Contributions

P.C.Y.W., S.K.P.L., and K.-Y.Y. conceived and designed the study and contributed reagents; J.L.L.T. coordinated the study; R.B., Y.T., A.Y.P.W., K.S.M.L., C.S.F.L., and R.Y.Y.F. performed the laboratory work; P.C.Y.W., J.L.L.T., and R.B. performed data analyses; P.C.Y.W. and J.L.L.T. wrote the manuscript. All authors read and approved the manuscript.

Funding

This work is partly supported by the Consultancy Service for Enhancing Laboratory Surveillance of Emerging Infectious Diseases and Research Capability on Antimicrobial Resistance for Department of Health, HKSAR Government, and the Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, the Ministry of Education of China.

Acknowledgments

We thank members of the AFCD and FEHD, HKSAR, for facilitation of and support with sample collection. We also thank Christopher M. Riggs and members of the Department of Veterinary Clinical Services, Hong Kong Jockey Club, for assistance in the collection of horse samples.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Pereira, H.G.; Fialho, A.M.; Flewett, T.H.; Teixeira, J.M.; Andrade, Z.P. Novel viruses in human faeces. Lancet 1988, 2, 103–104. [Google Scholar] [CrossRef]
  2. Pereira, H.G.; Flewett, T.H.; Candeias, J.A.; Barth, O.M. A virus with a bisegmented double-stranded RNA genome in rat (oryzomys nigripes) intestines. J. Gen. Virol. 1988, 69 Pt 11, 2749–2754. [Google Scholar] [CrossRef]
  3. Bodewes, R.; van der Giessen, J.; Haagmans, B.L.; Osterhaus, A.D.; Smits, S.L. Identification of multiple novel viruses, including a parvovirus and a hepevirus, in feces of red foxes. J. Virol. 2013, 87, 7758–7764. [Google Scholar] [CrossRef] [PubMed]
  4. Fregolente, M.C.; de Castro-Dias, E.; Martins, S.S.; Spilki, F.R.; Allegretti, S.M.; Gatti, M.S. Molecular characterization of picobirnaviruses from new hosts. Virus Res. 2009, 143, 134–136. [Google Scholar] [CrossRef] [PubMed]
  5. Ganesh, B.; Nagashima, S.; Ghosh, S.; Nataraju, S.M.; Rajendran, K.; Manna, B.; Ramamurthy, T.; Niyogi, S.K.; Kanungo, S.; Sur, D.; et al. Detection and molecular characterization of multiple strains of picobirnavirus causing mixed infection in a diarrhoeic child: Emergence of prototype genogroup II-like strain in kolkata, india. Int. J. Mol. Epidemiol. Genet. 2011, 2, 61–72. [Google Scholar] [PubMed]
  6. Ghosh, S.; Kobayashi, N.; Nagashima, S.; Naik, T.N. Molecular characterization of full-length genomic segment 2 of a bovine picobirnavirus (PBV) strain: Evidence for high genetic diversity with genogroup I PBVs. J. Gen. Virol. 2009, 90, 2519–2524. [Google Scholar] [CrossRef] [PubMed]
  7. Gillman, L.; Sanchez, A.M.; Arbiza, J. Picobirnavirus in captive animals from uruguay: Identification of new hosts. Intervirology 2013, 56, 46–49. [Google Scholar] [CrossRef] [PubMed]
  8. Kluge, M.; Campos, F.S.; Tavares, M.; de Amorim, D.B.; Valdez, F.P.; Giongo, A.; Roehe, P.M.; Franco, A.C. Metagenomic survey of viral diversity obtained from feces of subantarctic and south american fur seals. PLoS ONE 2016, 11, e0151921. [Google Scholar] [CrossRef]
  9. Luo, X.L.; Lu, S.; Jin, D.; Yang, J.; Wu, S.S.; Xu, J. Marmota himalayana in the qinghai-tibetan plateau as a special host for bi-segmented and unsegmented picobirnaviruses. Emerg. Microbes Infect. 2018, 7, 20. [Google Scholar] [CrossRef]
  10. Malik, Y.S.; Chandrashekar, K.M.; Sharma, K.; Haq, A.A.; Vaid, N.; Chakravarti, S.; Batra, M.; Singh, R.; Pandey, A.B. Picobirnavirus detection in bovine and buffalo calves from foothills of himalaya and central india. Trop. Anim. Health Prod. 2011, 43, 1475–1478. [Google Scholar] [CrossRef]
  11. Martinez, L.C.; Masachessi, G.; Carruyo, G.; Ferreyra, L.J.; Barril, P.A.; Isa, M.B.; Giordano, M.O.; Ludert, J.E.; Nates, S.V. Picobirnavirus causes persistent infection in pigs. Infect. Genet. Evol. J. Mol. Epidemiol. Evol. Genet. Infect. Dis. 2010, 10, 984–988. [Google Scholar] [CrossRef] [PubMed]
  12. Ng, T.F.; Mesquita, J.R.; Nascimento, M.S.; Kondov, N.O.; Wong, W.; Reuter, G.; Knowles, N.J.; Vega, E.; Esona, M.D.; Deng, X.; et al. Feline fecal virome reveals novel and prevalent enteric viruses. Vet. Microbiol. 2014, 171, 102–111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Pankovics, P.; Boros, A.; Nemes, C.; Kapusinszky, B.; Delwart, E.; Reuter, G. Molecular characterization of a novel picobirnavirus in a chicken. Arch. Virol. 2018, 163, 3455–3458. [Google Scholar] [CrossRef] [PubMed]
  14. Ribeiro Silva, R.; Bezerra, D.A.; Kaiano, J.H.; Oliveira Dde, S.; Silvestre, R.V.; Gabbay, Y.B.; Ganesh, B.; Mascarenhas, J.D. Genogroup I avian picobirnavirus detected in brazilian broiler chickens: A molecular epidemiology study. J. Gen. Virol. 2014, 95, 117–122. [Google Scholar] [CrossRef] [PubMed]
  15. Shi, M.; Lin, X.D.; Tian, J.H.; Chen, L.J.; Chen, X.; Li, C.X.; Qin, X.C.; Li, J.; Cao, J.P.; Eden, J.S.; et al. Redefining the invertebrate rna virosphere. Nature 2016, 540, 539–543. [Google Scholar] [CrossRef] [PubMed]
  16. Symonds, E.M.; Griffin, D.W.; Breitbart, M. Eukaryotic viruses in wastewater samples from the United States. Appl. Environ. Microbiol. 2009, 75, 1402–1409. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, Y.; Banyai, K.; Tu, X.; Jiang, B. Simian genogroup I picobirnaviruses: Prevalence, genetic diversity, and zoonotic potential. J. Clin. Microbiol. 2012, 50, 2779–2782. [Google Scholar] [CrossRef]
  18. Woo, P.C.; Teng, J.L.; Bai, R.; Wong, A.Y.; Martelli, P.; Hui, S.W.; Tsang, A.K.; Lau, C.C.; Ahmed, S.S.; Yip, C.C.; et al. High diversity of genogroup i picobirnaviruses in mammals. Front. Microbiol. 2016, 7, 1886. [Google Scholar] [CrossRef] [PubMed]
  19. Zhang, W.; Li, L.; Deng, X.; Kapusinszky, B.; Pesavento, P.A.; Delwart, E. Faecal virome of cats in an animal shelter. J. Gen. Virol. 2014, 95, 2553–2564. [Google Scholar] [CrossRef]
  20. Woo, P.C.; Lau, S.K.; Bai, R.; Teng, J.L.; Lee, P.; Martelli, P.; Hui, S.W.; Yuen, K.Y. Complete genome sequence of a novel picobirnavirus, otarine picobirnavirus, discovered in california sea lions. J. Virol. 2012, 86, 6377–6378. [Google Scholar] [CrossRef]
  21. Woo, P.C.; Lau, S.K.; Teng, J.L.; Tsang, A.K.; Joseph, M.; Wong, E.Y.; Tang, Y.; Sivakumar, S.; Bai, R.; Wernery, R.; et al. Metagenomic analysis of viromes of dromedary camel fecal samples reveals large number and high diversity of circoviruses and picobirnaviruses. Virology 2014, 471–473, 117–125. [Google Scholar] [CrossRef] [PubMed]
  22. Wakuda, M.; Pongsuwanna, Y.; Taniguchi, K. Complete nucleotide sequences of two RNA segments of human picobirnavirus. J. Virol. Methods 2005, 126, 165–169. [Google Scholar] [CrossRef] [PubMed]
  23. Rosen, B.I.; Fang, Z.Y.; Glass, R.I.; Monroe, S.S. Cloning of human picobirnavirus genomic segments and development of an rt-pcr detection assay. Virology 2000, 277, 316–329. [Google Scholar] [CrossRef] [PubMed]
  24. Delmas, B.; Attoui, H.; Ghosh, S.; Malik, Y.S.; Mundt, E.; Vakharia, V.N.; ICTV Report Consortium. ICTV virus taxonomy profile: Picobirnaviridae. J. Gen. Virol. 2019, 100, 133–134. [Google Scholar] [CrossRef] [PubMed]
  25. Kirkwood, C.D. Genetic and antigenic diversity of human rotaviruses: Potential impact on vaccination programs. J. Infect. Dis. 2010, 202 (Suppl. 1), S43–S48. [Google Scholar] [CrossRef]
  26. Komoto, S.; Apondi, E.W.; Shah, M.; Odoyo, E.; Nyangao, J.; Tomita, M.; Wakuda, M.; Maeno, Y.; Shirato, H.; Tsuji, T.; et al. Whole genomic analysis of human G12P[6] and G12P[8] rotavirus strains that have emerged in kenya: Identification of porcine-like NSP4 genes. Infect. Genet. Evol. 2014, 27, 277–293. [Google Scholar] [CrossRef] [PubMed]
  27. Smits, S.L.; Poon, L.L.; van Leeuwen, M.; Lau, P.N.; Perera, H.K.; Peiris, J.S.; Simon, J.H.; Osterhaus, A.D. Genogroup I and II picobirnaviruses in respiratory tracts of pigs. Emerg. Infect. Dis. 2011, 17, 2328–2330. [Google Scholar] [CrossRef] [PubMed]
  28. Smits, S.L.; van Leeuwen, M.; Schapendonk, C.M.; Schurch, A.C.; Bodewes, R.; Haagmans, B.L.; Osterhaus, A.D. Picobirnaviruses in the human respiratory tract. Emerg. Infect. Dis. 2012, 18, 1539–1540. [Google Scholar] [CrossRef] [PubMed]
  29. Lau, S.K.; Woo, P.C.; Li, K.S.; Huang, Y.; Tsoi, H.W.; Wong, B.H.; Wong, S.S.; Leung, S.Y.; Chan, K.H.; Yuen, K.Y. Severe acute respiratory syndrome coronavirus-like virus in chinese horseshoe bats. Proc. Natl. Acad. Sci. USA 2005, 102, 14040–14045. [Google Scholar] [CrossRef]
  30. Kumar, S.; Nei, M.; Dudley, J.; Tamura, K. Mega: A biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief. Bioinform. 2008, 9, 299–306. [Google Scholar] [CrossRef]
  31. Attoui, H.; Billoir, F.; Cantaloube, J.F.; Biagini, P.; de Micco, P.; de Lamballerie, X. Strategies for the sequence determination of viral dsrna genomes. J. Virol. Methods 2000, 89, 147–158. [Google Scholar] [CrossRef]
  32. Lambden, P.R.; Cooke, S.J.; Caul, E.O.; Clarke, I.N. Cloning of noncultivatable human rotavirus by single primer amplification. J. Virol. 1992, 66, 1817–1822. [Google Scholar] [PubMed]
  33. Hofacker, I.L. Vienna RNA secondary structure server. Nucleic Acids Res. 2003, 31, 3429–3431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Nates, S.V.; Gatti, M.S.V.; Ludert, J.E. The picobirnavirus: An integrated view on its biology, epidemiology and pathogenic potential. Future Virol. 2011, 6, 223–235. [Google Scholar] [CrossRef]
  35. Ng, T.F.; Kondov, N.O.; Deng, X.; Van Eenennaam, A.; Neibergs, H.L.; Delwart, E. A metagenomics and case-control study to identify viruses associated with bovine respiratory disease. J. Virol. 2015, 89, 5340–5349. [Google Scholar] [CrossRef] [PubMed]
  36. Chen, M.; Sun, H.; Lan, D.; Hua, X.; Cui, L.; Yuan, C.; Yang, Z. Molecular detection of genogroup I and II picobirnaviruses in pigs in china. Virus Genes 2014, 48, 553–556. [Google Scholar] [CrossRef] [PubMed]
  37. Ganesh, B.; Nataraju, S.M.; Rajendran, K.; Ramamurthy, T.; Kanungo, S.; Manna, B.; Nagashima, S.; Sur, D.; Kobayashi, N.; Krishnan, T. Detection of closely related picobirnaviruses among diarrhoeic children in kolkata: Evidence of zoonoses? Infect. Genet. Evol. 2010, 10, 511–516. [Google Scholar] [CrossRef] [PubMed]
  38. Van Leeuwen, M.; Williams, M.M.; Koraka, P.; Simon, J.H.; Smits, S.L.; Osterhaus, A.D. Human picobirnaviruses identified by molecular screening of diarrhea samples. J. Clin. Microbiol. 2010, 48, 1787–1794. [Google Scholar] [CrossRef]
  39. Banyai, K.; Martella, V.; Bogdan, A.; Forgach, P.; Jakab, F.; Meleg, E.; Biro, H.; Melegh, B.; Szucs, G. Genogroup I picobirnaviruses in pigs: Evidence for genetic diversity and relatedness to human strains. J. Gen. Virol. 2008, 89, 534–539. [Google Scholar] [CrossRef]
  40. Duraisamy, R.; Akiana, J.; Davoust, B.; Mediannikov, O.; Michelle, C.; Robert, C.; Parra, H.J.; Raoult, D.; Biagini, P.; Desnues, C. Detection of novel rna viruses from free-living gorillas, republic of the congo: Genetic diversity of picobirnaviruses. Virus Genes 2018, 54, 256–271. [Google Scholar] [CrossRef]
  41. Ghosh, S.; Shiokawa, K.; Aung, M.S.; Malik, Y.S.; Kobayashi, N. High detection rates of picobirnaviruses in free roaming rats (Rattus spp.): Molecular characterization of complete gene segment-2. Infect. Genet. Evol. 2018, 65, 131–135. [Google Scholar] [CrossRef] [PubMed]
  42. Kuhar, U.; Jamnikar-Ciglenecki, U. High detection rate and high genetic diversity of genogroup I picobirnaviruses from roe deer. Infect. Genet. Evol. J. Mol. Epidemiol. Evol. Genet. Infect. Dis. 2019, 73, 210–213. [Google Scholar] [CrossRef] [PubMed]
  43. Malik, Y.S.; Kumar, N.; Sharma, K.; Sharma, A.K.; Sircar, S.; Jeena, L.M.; Singh, N.K.; Mondal, A.; Joardar, S.N.; Balasubramanian, G. Molecular characterization of a genetically diverse bubaline picobirnavirus strain, india. Thai J. Vet. Med. 2013, 43, 609–613. [Google Scholar]
  44. Malik, Y.S.; Sharma, A.K.; Kumar, N.; Sharma, K.; Ganesh, B.; Kobayashi, N. Identification and characterisation of a novel genogroup ii picobirnavirus in a calf in india. Vet. Rec. 2014, 174, 278. [Google Scholar] [CrossRef] [PubMed]
  45. Malik, Y.S.; Sircar, S.; Dhama, K.; Singh, R.; Ghosh, S.; Banyai, K.; Vlasova, A.N.; Nadia, T.; Singh, R.K. Molecular epidemiology and characterization of picobirnaviruses in small ruminant populations in india. Infect. Genet. Evol. J. Mol. Epidemiol. Evol. Genet. Infect. Dis. 2018, 63, 39–42. [Google Scholar] [CrossRef] [PubMed]
  46. Bhattacharya, R.; Sahoo, G.C.; Nayak, M.K.; Rajendran, K.; Dutta, P.; Mitra, U.; Bhattacharya, M.K.; Naik, T.N.; Bhattacharya, S.K.; Krishnan, T. Detection of genogroup i and ii human picobirnaviruses showing small genomic rna profile causing acute watery diarrhoea among children in kolkata, india. Infect. Genet. Evol. 2007, 7, 229–238. [Google Scholar] [CrossRef] [PubMed]
  47. Krishnamurthy, S.R.; Wang, D. Extensive conservation of prokaryotic ribosomal binding sites in known and novel picobirnaviruses. Virology 2018, 516, 108–114. [Google Scholar] [CrossRef]
Figure 1. Phylogenetic analysis of a partial RdRp gene of genogroup I and II picobirnaviruses (PBVs) discovered in the present study. The trees were constructed using the neighbor-joining method. Bootstrap values below 70% are not shown. The scale bar indicates the number of nucleotide substitutions per site. All the accession numbers are given as cited in GenBank. (a) One hundred and ninety-one nucleotide positions were included in the analysis and rooted with genogroup II human strains, shown in italics. All PBV strains detected in the same host species were highlighted in the same color, with the highlights in bold representing the PBVs discovered in respiratory specimens. If more than one sequence type was found in the same sample, each sequence type was numbered in the order of identification (e.g., Monkey/11T-1 and Monkey/11T-2 indicated that there were two sequence types found in the same monkey sample 11T). (b) One hundred and ninety-two nucleotide positions were included in the analysis and rooted with genogroup I human PBV strains. The genogroup II PBV characterized in this study is shown in bold.
Figure 1. Phylogenetic analysis of a partial RdRp gene of genogroup I and II picobirnaviruses (PBVs) discovered in the present study. The trees were constructed using the neighbor-joining method. Bootstrap values below 70% are not shown. The scale bar indicates the number of nucleotide substitutions per site. All the accession numbers are given as cited in GenBank. (a) One hundred and ninety-one nucleotide positions were included in the analysis and rooted with genogroup II human strains, shown in italics. All PBV strains detected in the same host species were highlighted in the same color, with the highlights in bold representing the PBVs discovered in respiratory specimens. If more than one sequence type was found in the same sample, each sequence type was numbered in the order of identification (e.g., Monkey/11T-1 and Monkey/11T-2 indicated that there were two sequence types found in the same monkey sample 11T). (b) One hundred and ninety-two nucleotide positions were included in the analysis and rooted with genogroup I human PBV strains. The genogroup II PBV characterized in this study is shown in bold.
Viruses 11 00574 g001
Figure 2. Phylogenetic analysis of RdRp proteins of cattle and monkey picobirnaviruses (PBVs). The trees were constructed using the neighbor-joining method. Four hundred and seventy-five amino acid positions were included in the analysis. The scale bar indicates 0.05 substitutions per site. Bootstrap values below 70% are not shown. All PBV strains detected in the same host were highlighted with the solid dots/circles in the same color, with the highlights in bold representing the PBVs discovered in respiratory specimens. If more than one sequence type was found in the same sample, each sequence type was numbered in the order of identification (e.g., monkey/11T-1 and monkey/11T-2 indicated that there were two sequence types found in the same monkey sample 11T). The clade of each monkey PBV and their sampling locations (EN1: Lion Rock Site 1 and KS2: Kam Shan Site 2) are also indicated. All the accession numbers are given as cited in GenBank. GI and GII represent genogroup I and genogroup II PBVs, respectively.
Figure 2. Phylogenetic analysis of RdRp proteins of cattle and monkey picobirnaviruses (PBVs). The trees were constructed using the neighbor-joining method. Four hundred and seventy-five amino acid positions were included in the analysis. The scale bar indicates 0.05 substitutions per site. Bootstrap values below 70% are not shown. All PBV strains detected in the same host were highlighted with the solid dots/circles in the same color, with the highlights in bold representing the PBVs discovered in respiratory specimens. If more than one sequence type was found in the same sample, each sequence type was numbered in the order of identification (e.g., monkey/11T-1 and monkey/11T-2 indicated that there were two sequence types found in the same monkey sample 11T). The clade of each monkey PBV and their sampling locations (EN1: Lion Rock Site 1 and KS2: Kam Shan Site 2) are also indicated. All the accession numbers are given as cited in GenBank. GI and GII represent genogroup I and genogroup II PBVs, respectively.
Viruses 11 00574 g002
Table 1. Prevalence of picobirnaviruses (PBVs) in respiratory samples of animals by RT-PCR.
Table 1. Prevalence of picobirnaviruses (PBVs) in respiratory samples of animals by RT-PCR.
Animal TypeSample TypesNo. of Specimens ObtainedNo. (%) of Specimens Positive for Genogroup I PBVsNo. (%) of Specimens Positive for Genogroup II PBVs
Mammals
CattleNasopharyngeal swab513 (5.9%)1 (1.9%)
Monkey 1Throat swab523 (5.8%)0 (0%)
DogNasopharyngeal swab500 (0%)0 (0%)
CatNasopharyngeal swab500 (0%)0 (0%)
Bat 2Mouth swab1570 (0%)0 (0%)
Horse 3Nasopharyngeal swab1160 (0%)0 (0%)
Poultry
ChickenTracheal swab1000 (0%)0 (0%)
1 Most monkeys found in our locality are considered hybrids of Macaca mulatta and Macaca fascicularis. 2 Species included Hipposideros pomona (n = 103), Rhinolophus affinis (n = 14), and Rhinolophus sinicus (n = 40). 3 Thoroughbred racehorses stabled at the Hong Kong Jockey Club.
Table 2. Comparison of genomic features among monkey and cattle picobirnaviruses (PBVs) and other PBV strains with a complete/near complete segment 2 sequence available.
Table 2. Comparison of genomic features among monkey and cattle picobirnaviruses (PBVs) and other PBV strains with a complete/near complete segment 2 sequence available.
PBV Genome Features
Size (bp)G+C Content (%)ORF Features5′ UTR Features3′ UTR Features
ORFLocation (nt)Length (nt)Length (aa)FrameLength (nt)G+C Content (%)Beginning BasesLength (nt)G+C Content (%)Ending Bases
Reference genogroup I PBVs segment 2 sequences
Human 1 174546.36RdRp94-16981605534+19321.51GTAAA4753.19CTGC
Human 2 169643.63RdRp58-16501593530+15724.56GTAAA4650.00CTGC
Reference genogroup II PBV segment 2 sequence
Human 3 167443.91RdRp88-16411554517+18716.09GTAAA3357.58CTC
Genogroup I PBVs detected in this study
PBVs (Clade)GenBank accession no.
Monkey/11T-1 (M-R1;)KY120190167842.55RdRp28-16321605534+12714.81GTAAA4652.17CTGC
Monkey/11T-2 (M-R2)KY120191168141.05RdRp47-16361590529+24626.09GTAAA4551.11CTGC
Monkey/35T (M-R1)KY120193167842.61RdRp28-16321605534+12714.81GTAAA4652.17CTGC
Monkey/21T (M-R1)KY120192167842.49RdRp28-16321605534+12714.81GTAAA4652.17CTGC
Monkey/1R-1 4 (M-R3)KY120194161342.72RdRp2-15641563520NANANANA4953.06CTGC
Monkey/1R-2 (M-R2)KY120182168141.52RdRp47-16361590529+24626.09GTAAA4551.11CTGC
Monkey/2R (M-R1)KY120183167842.37RdRp28-16321605534+12714.81GTAAA4650.00CTGC
Monkey/12R (M-R1)KY120184167842.61RdRp28-16321605534+12714.81GTAAA4652.17CTGC
Monkey/13R-1 (M-R3)KY120185173441.64RdRp75-16851611536+37425.68GTAAA4953.06CTGC
Monkey/13R-2 (M-R1)KY120186167842.55RdRp28-16321605534+12714.81GTAAA4652.17CTGC
Monkey/14R (M-R1)KY120187167842.55RdRp28-16321605534+12714.81GTAAA4652.17CTGC
Monkey/15R-1 (M-R1)KY120188167842.31RdRp28-16321605534+12714.81GTAAA4650.00CTGC
Monkey/15R-2 (M-R3)KY120189173441.64RdRp75-16851611536+37425.68GTAAA4953.06CTGC
Cattle/C345N-2KY120170167143.93RdRp49-16291581526+14825.00GTAAA4252.38CATC
Cattle/C345N-3KY120171169442.15RdRp44-16571614537+24323.26GTAAA3759.46CTGC
Cattle/C346N-1KY120172170147.68RdRp49-16561608535+14829.17GTAAA4544.44CTGC
Cattle/C346N-2KY120173169644.99RdRp52-16471596531+15127.45GTAAA4946.94CTGC
Cattle/C361NKY120174169644.58RdRp49-16531605534+14835.42GTAAA4344.19CTGC
Cattle/C343R-1KY120175169244.03RdRp50-16421593530+24920.41GTAAA5058.00CTGC
Cattle/C343R-2KY120176169545.25RdRp48-16521605534+34725.53GTAAA4344.19CTGC
Cattle/C369R-1KY120177173642.45RdRp36-17211686561+33522.86GTAAA1540.00CTTC
Cattle/C369R-2 5KY120179165944.91RdRp24-16161593530NA2330.43NA4355.81CTGC
Cattle/C374R-1 5KY120180165144.16RdRp8-16061599532NA728.57NA4555.56CTGC
Cattle/C374R-2 5KY120181166947.15RdRp24-16251602533NA2334.78NA4447.73CTGC
Genogroup II PBVs detected in this study
Cattle/C372NKY120178162240.44RdRp55-15991545514+15416.67%GTAAA2352.17CTCC
1 The accession number of human PBV segment 2 is NC_007027; 2 the accession number of human PBV segment 2 is AF246939; 3 the accession number of human PBV segment 2 is AF246940; 4 the partial segment 2 with a near complete RdRp sequence; 5 the partial segment 2 with complete RdRp sequences; NA, not available; T, throat swab; R, rectal swab; N, nasal swab.

Share and Cite

MDPI and ACS Style

Woo, P.C.Y.; Teng, J.L.L.; Bai, R.; Tang, Y.; Wong, A.Y.P.; Li, K.S.M.; Lam, C.S.F.; Fan, R.Y.Y.; Lau, S.K.P.; Yuen, K.-Y. Novel Picobirnaviruses in Respiratory and Alimentary Tracts of Cattle and Monkeys with Large Intra- and Inter-Host Diversity. Viruses 2019, 11, 574. https://doi.org/10.3390/v11060574

AMA Style

Woo PCY, Teng JLL, Bai R, Tang Y, Wong AYP, Li KSM, Lam CSF, Fan RYY, Lau SKP, Yuen K-Y. Novel Picobirnaviruses in Respiratory and Alimentary Tracts of Cattle and Monkeys with Large Intra- and Inter-Host Diversity. Viruses. 2019; 11(6):574. https://doi.org/10.3390/v11060574

Chicago/Turabian Style

Woo, Patrick C. Y., Jade L. L. Teng, Ru Bai, Ying Tang, Annette Y. P. Wong, Kenneth S. M. Li, Carol S. F. Lam, Rachel Y. Y. Fan, Susanna K. P. Lau, and Kwok-Yung Yuen. 2019. "Novel Picobirnaviruses in Respiratory and Alimentary Tracts of Cattle and Monkeys with Large Intra- and Inter-Host Diversity" Viruses 11, no. 6: 574. https://doi.org/10.3390/v11060574

APA Style

Woo, P. C. Y., Teng, J. L. L., Bai, R., Tang, Y., Wong, A. Y. P., Li, K. S. M., Lam, C. S. F., Fan, R. Y. Y., Lau, S. K. P., & Yuen, K. -Y. (2019). Novel Picobirnaviruses in Respiratory and Alimentary Tracts of Cattle and Monkeys with Large Intra- and Inter-Host Diversity. Viruses, 11(6), 574. https://doi.org/10.3390/v11060574

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