Next Article in Journal
Hepatitis C Virus and Hepatitis B Virus Co-Infection
Previous Article in Journal
Infection of a Lepidopteran Cell Line with Deformed Wing Virus
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 7
10.3390/v12070740
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
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pteropine Orthoreovirus in an Angolan Soft-Furred Fruit Bat (Lissonycteris angolensis) in Uganda Dramatically Expands the Global Distribution of an Emerging Bat-Borne Respiratory Virus

by
Andrew J. Bennett
1 and
Tony L. Goldberg
1,2,*
1
Department of Pathobiological Sciences, University of Wisconsin-Madison, Madison, WI 53706, USA
2
Global Health Institute, University of Wisconsin-Madison, Madison, WI 53706, USA
*
Author to whom correspondence should be addressed.
Viruses 2020, 12(7), 740; https://doi.org/10.3390/v12070740
Submission received: 2 June 2020 / Revised: 2 July 2020 / Accepted: 7 July 2020 / Published: 9 July 2020
(This article belongs to the Section Animal Viruses)
Browse Figures
Versions Notes

Abstract

:
Pteropine orthoreovirus (PRV; Reoviridae: Spinareovirinae) is an emerging bat-borne zoonotic virus that causes influenza-like illness (ILI). PRV has thus far been found only in Australia and Asia, where diverse old-world fruit bats (Pteropodidae) serve as hosts. In this study, we report the discovery of PRV in Africa, in an Angolan soft-furred fruit bat (Lissonycteris angolensis ruwenzorii) from Bundibugyo District, Uganda. Metagenomic characterization of a rectal swab yielded 10 dsRNA genome segments, revealing this virus to cluster within the known diversity of PRV variants detected in bats and humans in Southeast Asia. Phylogeographic analyses revealed a correlation between geographic distance and genetic divergence of PRVs globally, which suggests a geographic continuum of PRV diversity spanning Southeast Asia to sub-Saharan Africa. The discovery of PRV in an African bat dramatically expands the geographic range of this zoonotic virus and warrants further surveillance for PRVs outside of Southeast Asia.

1. Introduction

Pteropine orthoreovirus (PRV; Reoviridae: Spinareovirinae) is a bat-borne zoonosis that causes upper respiratory tract infection (URTI), ranging from mild respiratory symptoms to influenza-like illness (ILI), in Southeast Asia [1,2,3,4]. PRV is a fusogenic orthoreovirus with 10 dsRNA segments (S1-4, M1-3, L1-3) [5]. PRV is only distantly related to other mammalian orthoreoviruses (MRV), and is nested within a clade of fusogenic viruses that cause significant pathology in reptiles and birds [5]. First isolated in 1968 from a flying fox (Pteropus poliocephalus) in Australia [6], variants of PRV have since been detected in diverse pteropodid bats across Asia [7,8,9,10] and in residents of, and travelers to, Southeast Asia [1,11,12,13]. Here, we describe the discovery of PRV in a bat from Equatorial Africa and discuss the implications of a dramatically expanded PRV geographic distribution.
The reoviruses (respiratory enteric orphan viruses; family Reoviridae) contain non-enveloped, segmented dsRNA viruses that cause distinct pathologies in diverse hosts [14]. The orthoreovirus genus is split into fusogenic and nonfusogenic viruses, determined by their capacity to form multinucleated syncytia mediated by uniquely simple fusion-associated small transmembrane (FAST) proteins for which there are no known homologs elsewhere in nature [15]. For nearly 30 years, the original variant of PRV, Nelson Bay virus (PRV1NB), remained a singular curiosity as the lone mammalian fusogenic orthoreovirus [16].
Since 1999, 15 PRVs (Table 1) were discovered as infecting pteropid bats in China, Indonesia, Malaysia, and the Philippines [7,8,9,10,17], and direct evidence of spillover to humans has been reported in Malaysia and Indonesia [1,11,12,13,18,19]. In 2006, Melaka virus (PRV3M), and shortly thereafter Kampar virus (PRV4K), were the first PRVs isolated from patients in Malaysia presenting to clinics with high fever, acute respiratory illness, vomiting, and diarrhea [18,19]. In each instance, the index case had had recent but indirect contact with bats, and in each instance human-to-human transmission occurred [18,19]. Subsequently, acute human PRV infections were identified in travelers who had visited caves in Bali, Indonesia [11,12]; in a Malaysian man with no clear bat exposure [20]; and in 34 out of 200 patients presenting with URTI to a Malaysian health clinic in the summer of 2012 [1]. Serological surveys suggest that zoonotic PRV infections may be common in Southeast Asia. For example, a survey of people in Tioman Island, Malaysia, detected antibodies in 13% of residents (n = 109; virus neutralization test) [18], and a survey of 272 hospital patients in Vietnam found that 4% had anti-PRV IgG antibodies by ELISA [21].
In June and July of 2017, we sampled Angolan soft-furred fruit bats (Lissonycteris angolensis ruwenzorii) in Kasitu Subcounty, Bundibugyo District, Western Region, Uganda, which was the site of a 2007–2008 outbreak of ebolavirus disease [22]. We sampled bats from a culvert beneath a road spanning a mountain stream that serves as a water source for nearby human communities. We then used metagenomic methods to survey for viruses that may have been shed into the stream. In one rectal swab, we identified the nearly complete genome of the first African bat PRV, which we provisionally named Kasama virus (PRV16K) (named for Kasitu Subcounty, Uganda, and the Konjo word for water, amaghetse).
The genome sequence of PRV16K is highly similar to that of Asian PRV variants and PRV16K clusters within the known diversity of Asian PRVs, despite its geographic distance from Asia (see below). In Southeast Asia, PRV variants appear to transmit readily between species of pteropid bats [23]. We hypothesize that this promiscuity could result in a multispecies geographic continuum of PRV diversity in pteropid bats from Southeast Asia to Equatorial Africa. In support of this hypothesis, we present evidence of isolation-by-distance (IBD) for PRV. Given the apparent frequency of PRV transmission from bats to humans in Southeast Asia, PRV may represent an unappreciated cause of ILI across Africa and Asia.

2. Materials and Methods

2.1. Collection of Samples

The bats were collected by mist-net (Avinet, Portland, ME, USA) as they exited their roost at dusk. The bats were roosting as a multispecies assemblage of ~200 Angolan soft-furred fruit bats (2 nights of double-observer roost exit counts) and an unknown number of Sundevall’s leaf-nosed bats (Hipposideros caffer). The roost was a stone and cement culvert beneath the Itojo-Sempaya Old Road in Bundibugyo District, Uganda (0°50’18.5”N and 30°10’22.1”E). The bats were held in individual cloth drawstring bags from the time of collection until sampling, they were provisioned with a sugar/water solution (30% m/v) following sample collection, and they were released on site.
Rectal swabs were collected from 11 Angolan soft-furred fruit bats using sterile rayon-polyester-tipped swabs, preserved in TRI Reagent (Zymo Research, Irvine, CA, USA), and frozen at −20 °C within 3 h of collection. All the samples were exported with permission of the Uganda Wildlife Authority and the Uganda National Council for Science and Technology, shipped in accordance with International Air Transport Association (IATA) regulations, and imported under permit number 2017-07-103 of the Centers for Disease Control and Prevention (Atlanta, GA, USA).

2.2. Metagenomics

Total RNA was extracted from the bat rectal swab samples using the Direct-zol RNA MicroPrep kit (Zymo Research, Irvine, CA, USA). Libraries were prepared for high-throughput sequencing (HTS) as previously described [24]. Briefly, sequencing libraries were prepared from 1 ng of cDNA per sample using the Nextera XT DNA (Illumina, San Diego, CA, USA) sample kit and sequenced using the Illumina MiSeq platform (Illumina, San Diego, CA, USA). HTS reads were quality-trimmed, and contiguous sequences (contigs) were assembled de novo using CLC Genomics Workbench version 11.1 (CLC Bio, Aarhus, Denmark) as previously described [25]. The assembled contigs were compared at the nucleotide and deduced protein sequence level to virus-derived and bat mitochondrial gene sequence in GenBank using BLASTn and BLASTx homology searching algorithms, respectively [26].

2.3. Phylogenetics of PRV16K Segments and Asia Pacific PRVs

PRV has 10 dsRNA genome segments (S1-3, M1-3, L1-3), nine of which contain a single open reading frame (ORF) (S2-3, M1-3, L1-3) and one (S1) with 3 overlapping ORFs [27]. The ORFs of the 10 segments of PRV16K were aligned (GenBank Accession numbers MT505315–MT505324) with homologous ORFs of the 15 previously identified strains of PRV [28] (Table 1; Genbank Accession numbers in Supplemental Table S1) using the Prank algorithm implemented in TranslatorX. The phylogenies were inferred using PhyML with “Bayesian Information Criterion” (BIC) smart model selection, “general-time reversible model with incorporation of rate of variation across site and proportion of invariable sites” (GTR+Γ+I) [29] and these are visualized in FigTree v.1.4.3 with midpoint rooting. The pairwise p-distances and standard errors (estimated from 1000 bootstraps) were calculated for all available PRV segments using MEGA7 [30].

2.4. Detection of Positive Selection on Cell-Attachment Protein (σC)

Sixteen sequences (PRV1NB–PRV16K; Table 1) of the S1 ORF encoding the cell-attachment protein (σC) were aligned as previously described [31]. The ratio of nonsynonymous (Ka) to synonymous (Ks) nucleotide substitution rates was estimated and a Ka/Ks annotated tree was generated using the Ka/Ks Calculation Tool of the University of Bergen Computational Biology Unit [32].

2.5. Correlation of Geographic Distance with Patristic Distance of Segments S1-4

To examine the patterns of the geographic spread of PRVs, we measured the pairwise geographic distance (minimizing distance traveled over water) between the site of sample collection/reported sites of origin for 15 PRVs (excluding PRV12I, which was isolated from a bat imported to Italy) (Supplemental Table S1) using ArcGIS Online (ESRI, Redlands, CA, USA). We computed the patristic distances between PRV variants from the phylogenetic trees of segments S1-S4 using PATRISTIC [33]. Mantel tests of the matrix correlation [34] were then performed to compare matrices of pairwise geographic distance to matrices of pairwise patristic distance for S1-S4 using XLSTAT [35]. The confidence intervals and the significance of Mantel’s r (rM) were computed using 10,000 Monte Carlo simulations (alpha = 0.05) and Mantel correlograms were generated using the Spatial Data Calculator [36] (20 intervals of 1000 km).

3. Results

3.1. Sequencing and PRV Genome Assembly

We collected rectal swabs from 11 bats, representing an estimated 5.5% of the approximately 200 Angolan soft-furred fruit bats in the roost. Metagenomic sequencing revealed the presence of PRV in one of these bats (prevalence 9.1%). In the infected bat, the average depth of coverage across the 10 PRV segments was 82.15 ± 32.7, with 1.2% of quality-trimmed reads mapping to a PRV segment. Gaps were filled via rt-PCR and Sanger sequencing. The coding-complete sequence of 9 out of 10 PRV segments was recovered. Based on the sequence data, PRV16K has the typical PRV genome organization including a tricistronic S1 and single ORFs on the remaining segments. For L3, 706 nt of the ORF encoding the major inner capsid protein (λ1) (nt 2544–3249 of L1) was lost to low coverage and could not be recovered by rt-PCR. This gap in the sequence alignment of L3 was excluded for all subsequent analyses. The mitochondrial gene sequences of cytochrome b (cytb) and cytochrome c oxidase subunit 1 (cox1) from the infected bat were also assembled as described above. These sequences (GenBank accession numbers MT708578 and MT707261) confirmed the bat to be an Angolan soft-furred bat of the Ruwenzorii subspecies (Lissonycteris angolensis ruwenzorii) (percent identity: cytb 99.56; cox1 100%).

3.2. Evolutionary Relationships among PRV Segments

Nucleotide-level analyses revealed different patterns for the 10 PRV segments (Figure 1). PRV S1 demonstrates the greatest diversity globally with an average pairwise distance of 35% ± 12.7%. The tricistronic S1 segment encodes the cell-attachment protein (σC), as well as two nonstructural proteins (p10: fusion-associated small transmembrane (FAST) protein, and p17 (unknown function in PRV)) [37]. Much of the global PRV S1 nucleotide diversity is found in the σC ORF (overall mean p-distance 40.6% ± 0.7%) as compared to P10 (24.1.x% ± 1.4%) and P17 (30.4% ±1.5%). PRV16K form a clade with bat variants PRV1NB, PRV12I, and PRV15G, but clusters with different human and bat variants across the remaining segments (Figure 2).
PRV16K S1 is most closely related to S1 of PRV12I, from a bat imported to Italy from Indonesia [9] (p-distance = 25.3% ± 1.1%). For the remaining segments, p-distance is remarkably low between PRV16K and known variants from humans (PRV3M, L3 (6.6% ± 0.5%); PRV4k, S2 (6.4% ± 0.4%), M2 (5.3% ± 0.5%), M3 (7.6% ± 0.5%); PRV8B, S3 (9.0% ± 0.6%); PRV10M, L2 (13.2% ± 0.6%)) and from bats (PRV2P, L1 (8.7% ± 0.5%); PRV11C, M1 (4.8% ± 0.4%); PRV12I, S4 (4.8% ± 0.4%)) of Southeast Asia (Figure 1).

3.3. Geographic Relationships among PRV Segments S1-4

The location of origin of 14 PRV variants were mapped with the color of points indicating p-distance relative to PRV16K of S1-4 (Figure 3). Although p-distances generally increase with geographic distance for S1-S4, the low S1 pairwise distance between PRV16K and PRV1NB is an outlier for the segment. The global Mantel test of the correlation between geographic and patristic distance of PRV S1-4 indicated significant correlation for S1, S2, and S4 but not for S3, and correlation was strongest for S1 and S4 (Table 2). The mantel correlograms (Figure 4) show some distance decay of Mantel’s r (rM), but genetic distance significantly declines with geographic distance out to the greatest pairwise geographic distance (19,420 km) for S2 and S4. The correlogram for S1 shows a a steep decline in rM from 2000 to 5000 km, which suggests that forces other than geographic distance, such as diversifying selection, may be driving geographic clustering of S1 genetic diversity. The global Mantel test for S3 was not significant, but the mantel correlogram for S3 shows significant correlation at short and intermediate distance intervals (0–4000 km, 5000–9000 km; Figure 4).

3.4. Diversifying Selection on Cell-Attachment Protein (σC)

We detected strong positive selection (Ka/Ks = 1.231) at node 15, where the two major clades of PRV S1 diverged. The purifying selection predominates (Ka/Ks << 1) on all other branches aside from the more moderate positive selection (Ka/Ks = 1.007) on branch 2 of Node 7 (Supplementary Figure S1).

4. Discussion

PRV was first discovered in 1968, with the isolation of PRV1NB from a grey-headed flying fox (Pteropus poliocephalus) [6]. PRV1NB was a novelty among mammalian orthoreoviruses, as it was capable of forming multinucleated syncytia via a FAST protein akin to the avian reoviruses (ARV) [5]. PRV2P, isolated from a small flying fox (Pteropus hypomelanus) in Malaysia in 1999, expanded the range of PRV in bats [10], but the recognition of PRV as a cause of respiratory illness in humans did not come until the discovery of PRV3M and PRV4K associated with URTI in Malaysia [18,19]. This realization led to further investigations showing that PRV may be transmitted regularly from bats to humans [1,20], placing it in a rare class of zoonotic bat viruses that transmit frequently from bats to humans [38]. The discovery of a closely related PRV strain in a Ugandan bat expands the geographic distribution of the virus to a new continent. The high similarity of PRV16K to Southeast Asian variants, and the strong and significant correlation between the geographic and genetic distances of PRV segments, evince that PRV16K is not the extreme outlier its far-flung origin would suggest, but perhaps a distant point on a geographic continuum of PRV diversity.
PRV16K, discovered in a rectal swab from an apparently healthy Angolan soft-furred fruit bat in Bundibugyo District, Uganda, is closely related to PRVs isolated from humans and bats in Asia Pacific (Figure 1). The ORFs of 8 out of 10 segments (S2-4, M1-3, L1, L2) all share ≥91% nucleotide identity with a known PRV variant. S1 identity is lower, but globally S1 appears to be less constrained than the remaining segments across PRVs. That no segments of PRV16K appear as an outlier to an Asian clade of PRVs (Figure 2) suggests that the virus did not co-diverge with African fruit bats from their ancestral Asian lineages in the Miocene [39], but rather spread between African and Asian bats more recently. The Angolan soft-furred bats are split into five subspecies. PRV16K was identified in the Ruwenzorii subspecies, which has been identified across sub-Saharan Africa [40].
The genetic and geographic distances between PRVs correlated significantly for segments S1, S2, and S4 by the global Mantel test (Table 2). Significant Mantel correlation supports an isolation-by-distance (IBD) hypothesis [41], wherein genetic distance between PRV variants increases with geographic distance across a distribution of host species [42]. For S2 and S4, Mantel correlograms (Figure 4) showed a relatively smooth decrease in rM, with significant correlation out to the greatest distance interval. This provides evidence that S2 and S4 diversity, across the geographic distribution of PRV, is not significantly structured by factors such as host species. In other words, PRV appears to have a broad host tropism [3]. We note that our analyses assume that human PRV infections were acquired locally, or, for imported cases, at the reported site of infection. Furthermore, PRV in bats may also be oversampled at some sites, which could cause other geographic biases. The strong association between PRV genetic and geographic distance that we documented may therefore change in magnitude as new viral variants are described from other locations.
Globally, PRV segments show different evolutionary histories (Figure 3), probably due to frequent inter-species transmission and genome reassortment during PRV co-infection [23]. A study of bat roosts in the Philippines found evidence of cross-species transmission of PRV isolates, neutralizing antibodies in the sera of infected and uninfected bats, and one bat apparently co-infected with two genetically divergent PRV isolates [7]. PRV has thus far been isolated from cave- and canopy-roosting frugivorous Pteropodid bats [43]. Cave-roosting bats show a higher propensity for interspecies virus transmission than do canopy roosting bats, due to the co-roosting of multiple species [44]. Although the mechanism by which PRV transmits between bat species across the range of the virus remains incompletely known, these studies suggested that co-roosting behavior is a major driver of such transmission. The transmission of PRV between related bat species should occur with relative ease under permissive ecological conditions, given the aforementioned zoonotic potential of PRV and further evidence of infection of wild crab-eating macaques (Macaca fascicularis) in Singapore [3] and Thailand [45]. We anticipate that future surveys of pteropids bats of south Asia, the Middle East, and the Horn of Africa may further expand the range of PRV and fill the gaps in the geographic continuum of genetic distance that we have documented.
The S1 Mantel correlogram (Figure 4) shows some evidence of genetic clustering, due to the steep decline in rM in the intervals from 2000 to 5000 km. We hypothesized that diversifying selection may have significantly impacted the evolutionary history of S1. We found evidence of strong positive selection on σC, and in particular the ratio of nonsynonymous to synonymous mutations showed strong positive selection acting on the node separating the two S1 clades (Supplemental Figure S1). σC is necessary for cell attachment and may determine cell/host tropism, although the host ligand(s) for PRV σC is not certain [37]. The structural and functional homolog in the more thoroughly studied ARV (σ1) induces high levels of strain specific anti-ARV neutralizing antibodies in birds [46], and so the strong positive selection on σC which orders the evolutionary history of S1 in PRV may be due to immune selection.
Although we do not yet have evidence that PRV16K is zoonotic, the high similarity of many PRV16K segments to PRVs isolated from humans with ILI suggests that further investigation is warranted. PRV infection dynamics are poorly understood, and the modes of transmission between bats, and to spillover hosts, are largely unknown, but the detection of virus RNA in bat excreta indicates that fecal–oral transmission could play a role. That the bat population in which PRV16K was discovered was associated with a human water source in Bundibugyo District, Uganda, is of concern. Reoviruses are notably stable in the environment due to their non-enveloped, “double-shelled” capsid [47]. Although waterborne transmission of PRV has not been reported, infectious mammalian reoviruses are frequently isolated from surface water [48,49]. Additionally, the original zoonotic variants of PRV (PRV3M, PRV4K) were acquired via indirect contact with bats [18,19], implicating environmental transmission of the viruses. We also note that Bundibugyo District was the site of an ebolavirus disease outbreak in 2007–2008 [22], suggesting a high risk for bat-associated zoonoses in this area.
Surveillance for influenza and ILI in many sub-Saharan African countries is poor, although the disease burden of ILI is significant [50,51,52,53]. Even when active influenza surveillance is undertaken, differential diagnosis requires multiplex diagnostics to accurately identify the diverse etiologies of ILI [54]. It is therefore entirely conceivable that PRV is a cause of URTI in Africa that has gone undetected. Furthermore, given the pan-equatorial African distribution of the Angolan soft-furred fruit bat and the apparent ease of transmission of PRV between bat species reported in Southeast Asia, PRV may be broadly distributed across the African continent. The discovery of PRV in Africa dramatically expands the range of this bat-borne zoonosis and warrants further investigation to assess whether PRV is actively transmitted in sub-Saharan Africa, as it is in Southeast Asia.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/12/7/740/s1, Figure S1: Ka/Ks annotated tree of cell-attachment protein ORF of Pteropine orthoreovirus Segment 1, Table S1: Accession numbers and locations of PRVs used in this study.

Author Contributions

Conceptualization, A.J.B. and T.L.G.; methodology, A.J.B and T.L.G.; formal analysis, A.J.B; investigation, A.J.B.; resources, A.J.B. and T.L.G.; data curation, A.J.B.; writing—original draft preparation, A.J.B.; writing—review and editing, A.J.B. and T.L.G.; visualization, A.J.B.; supervision, T.L.G.; funding acquisition, A.J.B. and T.L.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded through a grant from the University of Wisconsin-Madison Global Health Institute (A.J.B.), and by the University of Wisconsin-Madison through the John D. MacArthur Professorship Chair (T.L.G.).

Acknowledgments

We would like to thank David Hyeroba (Makerere University) for his invaluable intellectual and logistic contributions to this project (in memoriam). We also thank Jessica Carag (University of Wisconsin-Madison) for assistance with sample collection and Erica Moulton (University of Wisconsin-Madison) for assistance with the preparation of figures for this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Voon, K.; Tan, Y.F.; Leong, P.P.; Teng, C.L.; Gunnasekaran, R.; Ujang, K.; Chua, K.B.; Wang, L.-F. Pteropine orthoreovirus infection among out-patients with acute upper respiratory tract infection in Malaysia. J. Med. Virol. 2015, 87, 2149–2153. [Google Scholar] [CrossRef] [PubMed]
  2. Uehara, A.; Tan, C.W.; Mani, S.; Chua, K.B.; Leo, Y.S.; Anderson, D.E.; Wang, L.-F. Serological evidence of human infection by bat orthoreovirus in Singapore. J. Med. Virol. 2018, 91, 707–710. [Google Scholar] [CrossRef] [PubMed]
  3. Tan, C.W.; Wittwer, K.; Lim, X.F.; Uehara, A.; Mani, S.; Wang, L.-F.; Anderson, D.E. Serological evidence and experimental infection of cynomolgus macaques with pteropine orthoreovirus reveal monkeys as potential hosts for transmission to humans. Emerg. Microbes Infect. 2019, 8, 787–795. [Google Scholar] [CrossRef] [PubMed]
  4. Tan, Y.F.; Teng, C.L.; Chua, K.B.; Voon, K. Pteropine orthoreovirus: An important emerging virus causing infectious disease in the tropics? J. Infect. Dev. Ctries. 2017, 11, 215–219. [Google Scholar] [CrossRef] [Green Version]
  5. Day, J.M. The diversity of the orthoreoviruses: Molecular taxonomy and phylogentic divides. Infect. Genet. Evol. 2009, 9, 390–400. [Google Scholar] [CrossRef]
  6. Gard, G.P.; Marshall, I.D. Nelson Bay virus. Arch. Virol. 1973, 43, 34–42. [Google Scholar] [CrossRef]
  7. Taniguchi, S.; Maeda, K.; Horimoto, T.; Masangkay, J.S.; Puentespina, R.; Alvarez, J.; Eres, E.; Cosico, E.; Nagata, N.; Egawa, K.; et al. First isolation and characterization of pteropine orthoreoviruses in fruit bats in the Philippines. Arch. Virol. 2017, 162, 1529–1539. [Google Scholar] [CrossRef]
  8. Hu, T.-S.; Qiu, W.; He, B.; Zhang, Y.; Yu, J.; Liang, X.; Zhang, W.; Chen, G.; Zhang, Y.; Wang, Y.; et al. Characterization of a novel orthoreovirus isolated from fruit bat, China. BMC Microbiol. 2014, 14, 293. [Google Scholar] [CrossRef] [Green Version]
  9. Lorusso, A.; Teodori, L.; Leone, A.; Marcacci, M.; Mangone, I.; Orsini, M.; Dondona, A.C.; Cammà, C.; Monaco, F.; Savini, G. A new member of the Pteropine Orthoreovirus species isolated from fruit bats imported to Italy. Infect. Genet. Evol. 2015, 30, 55–58. [Google Scholar] [CrossRef] [Green Version]
  10. Pritchard, L.I.; Chua, K.B.; Cummins, D.; Hyatt, A.; Crameri, G.; Eaton, B.T.; Wang, L.-F. Pulau virus; A new member of the Nelson Bay orthoreovirus species isolated from fruit bats in Malaysia. Arch. Virol. 2005, 151, 229–239. [Google Scholar] [CrossRef]
  11. Wong, A.H.; Cheng, P.K.; Lai, M.Y.; Leung, P.C.; Wong, K.K.; Lee, W.; Lim, W.W. Virulence Potential of Fusogenic Orthoreoviruses. Emerg. Infect. Dis. 2012, 18, 944–948. [Google Scholar] [CrossRef] [PubMed]
  12. Cheng, P.; Lau, C.S.; Lai, A.; Ho, E.; Leung, P.; Chan, F.; Wong, A.; Lim, W. A novel reovirus isolated from a patient with acute respiratory disease. J. Clin. Virol. 2009, 45, 79–80. [Google Scholar] [CrossRef] [PubMed]
  13. Singh, H.; Yoshikawa, T.; Kobayashi, T.; Fukushi, S.; Tani, H.; Taniguchi, S.; Fukuma, A.; Yang, M.; Sugamata, M.; Shimojima, M.; et al. Rapid whole genome sequencing of Miyazaki-Bali/2007 Pteropine orthoreovirus by modified rolling circular amplification with adaptor ligation—Next generation sequencing. Sci. Rep. 2015, 5, 16517. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Reoviridae—DsRNA Viruses—DsRNA Viruses. 2011. Available online: https://talk.ictvonline.org/ictv-reports/ictv_9th_report/dsrna-viruses-2011/w/dsrna_viruses/188/reoviridae (accessed on 22 April 2020).
  15. Shmulevitz, M.; Duncan, R. A new class of fusion-associated small transmembrane (FAST) proteins encoded by the non-enveloped fusogenic reoviruses. EMBO J. 2000, 19, 902–912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Duncan, R.; Murphy, F.A.; Mirkovic, R.R. Characterization of a Novel Syncytium-Inducing Baboon Reovirus. Virology 1995, 212, 752–756. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Du, L.; Lu, Z.; Fan, Y.; Meng, K.; Jiang, Y.; Zhu, Y.; Wang, S.; Gu, W.; Zou, X.; Tu, C. Xi River virus, a new bat reovirus isolated in southern China. Arch. Virol. 2010, 155, 1295–1299. [Google Scholar] [CrossRef]
  18. Chua, K.B.; Crameri, G.; Hyatt, A.; Yu, M.; Tompang, M.R.; Rosli, J.; McEachern, J.; Crameri, S.; Kumarasamy, V.; Eaton, B.T.; et al. A previously unknown reovirus of bat origin is associated with an acute respiratory disease in humans. Proc. Natl. Acad. Sci. USA 2007, 104, 11424–11429. [Google Scholar] [CrossRef] [Green Version]
  19. Chua, K.B.; Voon, K.; Crameri, G.; Tan, H.S.; Rosli, J.; McEachern, J.A.; Suluraju, S.; Yu, M.; Wang, L.-F. Identification and Characterization of a New Orthoreovirus from Patients with Acute Respiratory Infections. PLoS ONE 2008, 3, e3803. [Google Scholar] [CrossRef]
  20. Chua, K.B.; Voon, K.; Yu, M.; Keniscope, C.; Rasid, K.A.; Wang, L.-F. Investigation of a Potential Zoonotic Transmission of Orthoreovirus Associated with Acute Influenza-Like Illness in an Adult Patient. PLoS ONE 2011, 6, e25434. [Google Scholar] [CrossRef] [Green Version]
  21. Singh, H.; Shimojima, M.; Ngoc, T.C.; Huy, N.V.Q.; Chuong, T.X.; le van, A.; Saijo, M.; Yang, M.; Sugamata, M. Serological evidence of human infection with Pteropine orthoreovirus in Central Vietnam. J. Med Virol. 2015, 87, 2145–2148. [Google Scholar] [CrossRef] [Green Version]
  22. Wamala, J.F.; Lukwago, L.; Malimbo, M.; Nguku, P.; Yoti, Z.; Musenero, M.; Amone, J.; Mbabazi, W.; Nanyunja, M.; Zaramba, S.; et al. Ebola Hemorrhagic Fever Associated with Novel Virus Strain, Uganda, 2007–2008. Emerg. Infect. Dis. 2010, 16, 1087–1092. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Voon, K.; Chua, K.B.; Yu, M.; Crameri, G.; Barr, J.; Malik, Y.; Wang, L.-F. Evolutionary relationship of the L- and M-class genome segments of bat-borne fusogenic orthoreoviruses in Malaysia and Australia. J. Gen. Virol. 2011, 92, 2930–2936. [Google Scholar] [CrossRef] [PubMed]
  24. Bennett, A.J.; Sibley, S.D.; Lauck, M.; Weny, G.; Hyeroba, D.; Tumukunde, A.; Friedrich, T.C.; O’Connor, D.H.; Johnson, C.A.; Rothman, J.M.; et al. Naturally Circulating Hepatitis A Virus in Olive Baboons, Uganda. Emerg. Infect. Dis. 2016, 22, 1308–1310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Toohey-Kurth, K.; Sibley, S.; Goldberg, T.L. Metagenomic assessment of adventitious viruses in commercial bovine sera. Biologicals 2017, 47, 64–68. [Google Scholar] [CrossRef]
  26. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  27. Schiff, L.A.; Nibert, M.L.; Tyler, K.L. Orthoreoviruses and their replication. In Fields Virology; Lippincott/Williams & Wilkins: Philadelphia, PA, USA, 2007; pp. 1853–1915. [Google Scholar]
  28. Wang, L.-F.; Anderson, D.E. Viruses in bats and potential spillover to animals and humans. Curr. Opin. Virol. 2019, 34, 79–89. [Google Scholar] [CrossRef]
  29. Lefort, V.; Longueville, J.-E.; Gascuel, O. SMS: Smart Model Selection in PhyML. Mol. Biol. Evol. 2017, 34, 2422–2424. [Google Scholar] [CrossRef] [Green Version]
  30. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  31. Löytynoja, A. Phylogeny-aware alignment with PRANK. In Advanced Structural Safety Studies; Springer Science and Business Media LLC: Berlin/Heidelberg, Germany, 2013; Volume 1079, pp. 155–170. [Google Scholar]
  32. Computational Biology Unit—Ka/Ks Calculation Tool. Available online: http://services.cbu.uib.no/tools/kaks (accessed on 20 May 2020).
  33. Fourment, M.; Gibbs, M.J. PATRISTIC: A program for calculating patristic distances and graphically comparing the components of genetic change. BMC Evol. Biol. 2006, 6, 1. [Google Scholar] [CrossRef] [Green Version]
  34. Bailey, K.; Sokal, R.R.; Rohlf, F. Biometry: The Principles and Practice of Statistics in Biological Research. J. Am. Stat. Assoc. 1982, 77, 946. [Google Scholar] [CrossRef]
  35. XLSTAT Statistical and Data Analysis Solution; Addinsoft: Long Island, NY, USA, 2019.
  36. Remm, K.; Kelviste, T. An online calculator for spatial data and its applications. Comput. Ecol. Softw. 2014, 4, 22–34. [Google Scholar]
  37. Kawagishi, T.; Kanai, Y.; Tani, H.; Shimojima, M.; Saijo, M.; Matsuura, Y.; Kobayashi, T. Reverse Genetics for Fusogenic Bat-Borne Orthoreovirus Associated with Acute Respiratory Tract Infections in Humans: Role of Outer Capsid Protein σC in Viral Replication and Pathogenesis. PLoS Pathog. 2016, 12, e1005455. [Google Scholar] [CrossRef] [Green Version]
  38. Plowright, R.; Eby, P.; Hudson, P.J.; Smith, I.L.; Westcott, D.; Bryden, W.L.; Middleton, D.; Reid, P.A.; McFarlane, R.; Martin, G.; et al. Ecological dynamics of emerging bat virus spillover. Proc. R. Soc. Biol. Sci. 2015, 282, 20142124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Almeida, F.C.; Giannini, N.P.; Simmons, N.B. The evolutionary history of the African fruit bats (Chiroptera: Pteropodidae). Acta Chiropterol. 2016, 18, 73–90. [Google Scholar] [CrossRef]
  40. African Chiroptera Report 2019. Available online: https://www.academia.edu/41472741/AFRICAN_CHIROPTERA_REPORT_2019 (accessed on 1 June 2020).
  41. Wright, S. Isolation by Distance. Genetics 1943, 28, 114–138. [Google Scholar] [PubMed]
  42. Tao, Y.; Shi, M.; Chommanard, C.; Queen, K.; Zhang, J.; Markotter, W.; Kuzmin, I.V.; Holmes, E.C.; Tong, S. Surveillance of Bat Coronaviruses in Kenya Identifies Relatives of Human Coronaviruses NL63 and 229E and Their Recombination History. J. Virol. 2017, 91, e01953-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wilson, D.E.; Reeder, D.M. Mammal Species of the World: A Taxonomic and Geographic Reference; JHU Press: Baltimore, MD, USA, 2005; ISBN 978-0-8018-8221-0. [Google Scholar]
  44. Willoughby, A.R.; Phelps, K.L.; Olival, K.J. PREDICT Consortium A Comparative Analysis of Viral Richness and Viral Sharing in Cave-Roosting Bats. Diversity 2017, 9, 35. [Google Scholar] [CrossRef] [Green Version]
  45. Kosoltanapiwat, N.; Reamtong, O.; Okabayashi, T.; Ampawong, S.; Rungruengkitkun, A.; Thiangtrongjit, T.; Thippornchai, N.; Leaungwutiwong, P.; Mahittikorn, A.; Mori, H.; et al. Mass spectrometry-based identification and whole-genome characterisation of the first pteropine orthoreovirus isolated from monkey faeces in Thailand. BMC Microbiol. 2018, 18, 135. [Google Scholar] [CrossRef]
  46. Wickramasinghe, R.; Meanger, J.; Enriquez, C.E.; Wilcox, G.E. Avian Reovirus Proteins Associated with Neutralization of Virus Infectivity. Virology 1993, 194, 688–696. [Google Scholar] [CrossRef]
  47. Payne, S. Family Reoviridae. In Viruses; Elsevier BV: Amsterdam, The Netherlands, 2017; pp. 219–226. [Google Scholar]
  48. Spinner, M.L.; di Giovanni, G.D. Detection and Identification of Mammalian Reoviruses in Surface Water by Combined Cell Culture and Reverse Transcription-PCR. Appl. Environ. Microbiol. 2001, 67, 3016–3020. [Google Scholar] [CrossRef] [Green Version]
  49. Lim, M.C.Y.; Wang, Y.-F.; Huang, S.-W.; Yang, J.-Y.; Wang, J.-R. High Incidence of Mammalian Orthoreovirus Identified by Environmental Surveillance in Taiwan. PLoS ONE 2015, 10, e0142745. [Google Scholar] [CrossRef] [PubMed]
  50. McMorrow, M.L.; Emukule, G.O.; Njuguna, H.N.; Bigogo, G.; Montgomery, J.M.; Nyawanda, B.; Audi, A.; Breiman, R.F.; Katz, M.A.; Cosmas, L.; et al. The Unrecognized Burden of Influenza in Young Kenyan Children, 2008–2012. PLoS ONE 2015, 10, e0138272. [Google Scholar] [CrossRef] [Green Version]
  51. Nair, H.; Brooks, W.A.; Katz, M.; Roca, A.; Berkley, J.A.; Madhi, S.A.; Simmerman, J.M.; Gordon, A.; Sato, M.; Howie, S.; et al. Global burden of respiratory infections due to seasonal influenza in young children: A systematic review and meta-analysis. Lancet 2011, 378, 1917–1930. [Google Scholar] [CrossRef] [Green Version]
  52. Katz, M.A.; Schoub, B.D.; Heraud, J.-M.; Breiman, R.F.; Njenga, M.K.; Widdowson, M.-A. Influenza in Africa: Uncovering the Epidemiology of a Long-Overlooked Disease. J. Infect. Dis. 2012, 206, S1–S4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Schoub, B. Surveillance and management of influenza on the African continent. Expert Rev. Respir. Med. 2010, 4, 167–169. [Google Scholar] [CrossRef]
  54. Razuri, H.; Romero, C.; Tinoco, Y.; Guezala, M.C.; Ortiz, E.; Silva, M.; Reaves, E.; Williams, M.; Laguna-Torres, V.A.; Halsey, E.S.; et al. Population-based active surveillance cohort studies for influenza: Lessons from Peru. Bull. World Heal. Organ. 2012, 90, 318–320. [Google Scholar] [CrossRef] [Green Version]
Viruses 12 00740 g001 550
Figure 1. Heatmap of p-distance for the segments of 15 previously isolated PRVs relative to PRV16K. The pairwise p-distance was estimated using MEGA7 (1000 bootstrap replicates). The GenBank accession numbers for PRV isolates are available in Supplemental Table S1. The missing segments of PRV isolates are indicated by black and white hashed cells. The outlined cells indicate the lowest p-distance to PRV16K for each segment. The heatmap color key shows the range of the p-distance values estimated.
Figure 1. Heatmap of p-distance for the segments of 15 previously isolated PRVs relative to PRV16K. The pairwise p-distance was estimated using MEGA7 (1000 bootstrap replicates). The GenBank accession numbers for PRV isolates are available in Supplemental Table S1. The missing segments of PRV isolates are indicated by black and white hashed cells. The outlined cells indicate the lowest p-distance to PRV16K for each segment. The heatmap color key shows the range of the p-distance values estimated.
Viruses 12 00740 g001
Viruses 12 00740 g002 550
Figure 2. Maximum likelihood trees of the 10 PRV genome segments. The trees were constructed from codon-based alignments of open reading frames from each segment (GenBank Accession numbers in Supplemental Table S1), with molecular evolution model selection (GTR+ +I) by PhyML smart model selection (BIC). Circles on nodes indicate >75% confidence based on 1000 bootstrap replicates. The scale bar indicates nucleotide substitutions per site.
Figure 2. Maximum likelihood trees of the 10 PRV genome segments. The trees were constructed from codon-based alignments of open reading frames from each segment (GenBank Accession numbers in Supplemental Table S1), with molecular evolution model selection (GTR+ +I) by PhyML smart model selection (BIC). Circles on nodes indicate >75% confidence based on 1000 bootstrap replicates. The scale bar indicates nucleotide substitutions per site.
Viruses 12 00740 g002
Viruses 12 00740 g003 550
Figure 3. Phylogeographic analysis of PRV segments S1-S4. The p-distances between PRV16K segments S1-4 and segments S1-S4 of 14 PRV isolates (MEGA7, 1000 bootstrap replicates) are indicated by hexagons indicating the locations of origin of PRV isolates, with hexagons colored to indicate p-distance values (scale bars indicate ranges for each segment). The numbers indicate PRV isolate: (1) PRV1NB, (2) PRV2P, (3) PRV3M, (4) PRV4K, (5) PRV5HK, (6) PRV6XR, (7) PRV7S, (8) PRV8B, (9) PRV9HK, (10) PRV10M, (11) PRV11C, (12) PRV12I, (13) PRV13P, (14) PRV14P, (15) PRV15G. The GenBank accession numbers are in Supplemental Table S1. The sequences for S2 and S4 were not available for (6) PRV6XR (indicated by a black hexagon).
Figure 3. Phylogeographic analysis of PRV segments S1-S4. The p-distances between PRV16K segments S1-4 and segments S1-S4 of 14 PRV isolates (MEGA7, 1000 bootstrap replicates) are indicated by hexagons indicating the locations of origin of PRV isolates, with hexagons colored to indicate p-distance values (scale bars indicate ranges for each segment). The numbers indicate PRV isolate: (1) PRV1NB, (2) PRV2P, (3) PRV3M, (4) PRV4K, (5) PRV5HK, (6) PRV6XR, (7) PRV7S, (8) PRV8B, (9) PRV9HK, (10) PRV10M, (11) PRV11C, (12) PRV12I, (13) PRV13P, (14) PRV14P, (15) PRV15G. The GenBank accession numbers are in Supplemental Table S1. The sequences for S2 and S4 were not available for (6) PRV6XR (indicated by a black hexagon).
Viruses 12 00740 g003
Viruses 12 00740 g004 550
Figure 4. Mantel correlograms for segments S1-4 of 16 PRV isolates. The mantel correlograms show the relationship between geographic distance (km) and the correlation between pairwise genetic and geographic distances (Mantel’s r, or rM) across 20 distance intervals (1000 km). The black lines show Mantel’s r computed for each interval, while the red lines show the results of 10,000 random permutations (randomizing matrix associations). For each interval, Mantel’s r values (black lines) outside of the limits of the results of 10,000 random permutations (red lines) are considered significant (alpha = 0.5).
Figure 4. Mantel correlograms for segments S1-4 of 16 PRV isolates. The mantel correlograms show the relationship between geographic distance (km) and the correlation between pairwise genetic and geographic distances (Mantel’s r, or rM) across 20 distance intervals (1000 km). The black lines show Mantel’s r computed for each interval, while the red lines show the results of 10,000 random permutations (randomizing matrix associations). For each interval, Mantel’s r values (black lines) outside of the limits of the results of 10,000 random permutations (red lines) are considered significant (alpha = 0.5).
Viruses 12 00740 g004
Table
Table 1. Pteropine orthoreoviruses detected in bats and humans.
Table 1. Pteropine orthoreoviruses detected in bats and humans.
IsolateCommon NameYearHostCountry of Origin
PRV1NBNelson Bay virus1968Bat (Pteropus policephalus)Australia
PRV2PPulau virus1999Bat (Pteropus hypomelanus)Malaysia
PRV3MMelaka virus2006HumanMalaysia
PRV4KKampar virus2006HumanMalaysia
PRV5HKHK23629/072007HumanIndonesia
PRV6XRXi River virus2006Bat (Rousettus leschenaultii)China
PRV7SSikamat virus2010HumanMalaysia
PRV8BHK46886/092009HumanIndonesia
PRV9HKHK50842/102010HumanIndonesia
PRV10MMiyazaki-Bali 20072007HumanIndonesia
PRV11CCangyuan virus2012Bat (Rousettus leschenaultii)China
PRV12IIndonesia/20102010Bat (Pteropus vampyrus)Indonesia
PRV13PSamal-242013Bat (Eonycteris spelaea)Philippines
PRV14PTalikud-802013Bat (Rousettus amplexicaudatus)Philippines
PRV15GGarut-692017Bat (Pteropus vampyrus)Indonesia
PRV16K*Kasama virus2017Bat (Lissonycteris angolensis ruwenzorii)Uganda
* identified in this study.
Table
Table 2. Mantel test showing the correlation of geographic and genetic distances for PRV S1, S2, and S4.
Table 2. Mantel test showing the correlation of geographic and genetic distances for PRV S1, S2, and S4.
PRV Genome SegmentVirus Protein(s)Mantel’s r (rM)p-Value (Two-Tailed)
[99% CI]
S1P10 (FAST),
P17 (unknown),
σC (cell attachment protein)		0.462<0.0001 *
S2σ1 (major inner-capsid protein)0.2290.031 [0.027–0.036] *
S3σNS (nonstructural replication protein)0.1320.179 [0.169–0.189]
S4σ2 (major outer-capsid protein)0.373<0.0001 *
* significant to α = 0.05.

Share and Cite

MDPI and ACS Style

Bennett, A.J.; Goldberg, T.L. Pteropine Orthoreovirus in an Angolan Soft-Furred Fruit Bat (Lissonycteris angolensis) in Uganda Dramatically Expands the Global Distribution of an Emerging Bat-Borne Respiratory Virus. Viruses 2020, 12, 740. https://doi.org/10.3390/v12070740

AMA Style

Bennett AJ, Goldberg TL. Pteropine Orthoreovirus in an Angolan Soft-Furred Fruit Bat (Lissonycteris angolensis) in Uganda Dramatically Expands the Global Distribution of an Emerging Bat-Borne Respiratory Virus. Viruses. 2020; 12(7):740. https://doi.org/10.3390/v12070740

Chicago/Turabian Style

Bennett, Andrew J., and Tony L. Goldberg. 2020. "Pteropine Orthoreovirus in an Angolan Soft-Furred Fruit Bat (Lissonycteris angolensis) in Uganda Dramatically Expands the Global Distribution of an Emerging Bat-Borne Respiratory Virus" Viruses 12, no. 7: 740. https://doi.org/10.3390/v12070740

APA Style

Bennett, A. J., & Goldberg, T. L. (2020). Pteropine Orthoreovirus in an Angolan Soft-Furred Fruit Bat (Lissonycteris angolensis) in Uganda Dramatically Expands the Global Distribution of an Emerging Bat-Borne Respiratory Virus. Viruses, 12(7), 740. https://doi.org/10.3390/v12070740

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