key: cord-003166-k3jxvzfi
authors: Noh, Ji Yeong; Jeong, Dae Gwin; Yoon, Sun-Woo; Kim, Ji Hyung; Choi, Yong Gun; Kang, Shien-Young; Kim, Hye Kwon
title: Isolation and characterization of novel bat paramyxovirus B16-40 potentially belonging to the proposed genus Shaanvirus
date: 2018-08-22
journal: Sci Rep
DOI: 10.1038/s41598-018-30319-7
sha: 
doc_id: 3166
cord_uid: k3jxvzfi

The bat paramyxovirus B16-40 was first isolated in Korea in this study. Using the isolated virus, we could obtain not only genomic information, but also several biological characteristics of the virus. In the phylogenetic analysis, the virus was found to belong to the recently proposed genus Shaanvirus. Through sequence analyses and in vitro testing, the isolated virus was also found to have haemagglutinin-neuraminidase (HN) protein as one of the structural proteins. When mouse antiserum was generated against the isolated virus and tested, it was cross-reactive to human parainfluenza virus 1 in an indirect immunofluorescence assay but could not cross-neutralize human parainfluenza virus 1. In addition, the bat paramyxovirus B16-40 was not infectious in the mouse model. Collectively, this study provided basic information on further classification of the bat paramyxovirus B16-40 and related viruses in the proposed genus Shaanvirus.

Bats are considered a reservoir of severe emerging infectious diseases. Severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), Nipah virus, Hendra virus, and Ebola virus are all thought to be bat-borne viruses 1, 2 .

Notably, bats also host major mammalian paramyxoviruses from the family Paramyxoviridae, order Monone-gavirales 3, 4 . While Henipaviruses (Nipah and Hendra viruses) in South East Asia and Australia are associated with fruit bats 5 , other paramyxoviruses have been detected not only in fruit bats but in insectivorous bats worldwide [6] [7] [8] [9] . A potential pathway for Nipah virus transmission from bats to humans was found to be associated with a human-bat interface, specifically date palm sap shared by bats and humans 10 . In addition, serological evidence of possible human infection with a bat-originated paramyxovirus, Tioman virus 11 , reinforces the epidemiological role of bats in the emergence of pathogens such as paramyxoviruses in humans.

In addition to these bat paramyxoviruses with zoonotic potential, other new paramyxoviruses have been reported. These include several new mammalian paramyxoviruses such as Beilong virus and J virus, which remain unassigned under the family Paramyxoviridae 12 . Recent bat-associated paramyxoviruses were proposed to be grouped in a separate phylogenetic clade within a potentially separate genus such as Shaanvirus 13 which was distantly related to Jeilongvirus 14 . In addition, novel strains of bat paramyxoviruses in diverse genera have been reported continuously [15] [16] [17] . Based on the recent papers, bat paramyxoviruses found worldwide to date have belonged to the genera Rubulavirus, Morbillivirus, Henipavirus and the unclassified proposed genera Shaanvirus. Expanded classifications for grouping newly identified viruses in bats can be accomplished by further studying the biological characteristics of novel paramyxoviruses as well as genome characterization 18 .

In this study, active surveillance was performed to reveal paramyxoviruses circulating in Korean bats. A total of 232 bat samples were collected at 48 sites in natural bat habitats and tested for the possible existence of paramyxoviruses.

Identification of novel bat paramyxovirus B16-40. In this study, 232 bat samples were collected at 48 sites from natural bat habitats and tested for the possible existence of paramyxoviruses (Supplementary Table S1 ). Based on the RT-semi-nested PCR using the consensus paramyxovirus primers targeting the RdRp region 19 , five samples were positive and confirmed to be bat paramyxoviruses by sequencing. B16-6 was from a urine sample of Miniopterus schreibersii at BT cave in Hapcheon in March 2016, and B16-40 was from a feces sample of Miniopterus schreibersii at G cave in Danyang in April 2016. B16-148 and B16-154 were from feces samples of Murina ussuriensis at OJ cave in Youngwol and at S cave in Pyungchang in July 2016, respectively. B16-230 was from a urine sample of Miniopterus schreibersii at BLR cave in Seogwipo in December 2016. Among these, one paramyxovirus (B16-40) was successfully isolated.

The isolated virus, bat paramyxovirus B16-40, was isolated from the Korean bat Miniopterus schreibersii. The virus was successfully cultivated after three blind passages in the MARC-145 cell line (Fig. 1) , which is a clone of the MA-104 cell line derived from kidney epithelial cells of an African green monkey 20 . Viral particles were observed in a pellet of infected MARC-145 cells under a transmission electron microscope (Fig. 1B) . In terms of viral growth, infectious viruses were increased in the culture supernatant and infected cells from 24 hours post inoculation (hpi) and were consistently released for 144 hpi, as measured using tissue culture infectious dose 50 (TCID 50 ) method (Fig. 1C ).

Infectivity in C57BL/6 mice. When 6-week-old female C57BL/6 mice were inoculated via the intranasal and intragastric routes with 10 5 Tissue Culture Infectious Dose (TCID) 50 /mL or 10 2.5 TCID 50 /mL of bat paramyxovirus B16-40, we found no evidence of infection, i.e. no viral shedding, histopathological findings, or seroconversion.

Characterization of the bat paramyxovirus B16-40. Using high throughput sequencing (HiSeq. 2000 sequencing system based on the transcriptome de novo sequencing platform), 23,336 contigs totaling 8,010,635 base pairs (bps) with an average length of 343 bps were obtained. In the taxonomy annotation with MG-RAST, 12,522 contigs were annotated. While most of the annotated sequences associated with eukaryotes, 13 sequences were associated with family Paramyxoviridae. The 17,771 bases of viral genomic sequences were obtained by RT-PCR using primers designed based on the Paramyxoviridae-associated contigs and RACE sequencing.

The genomic sequence of the virus encodes 8 Open Reading Frames (ORFs) ( Fig. 2A) , which is consistent with bat paramyxovirus strain Ms-ParaV/Anhui2011 (KC154054) 13 . In addition, maximum-likelihood trees showed that the bat paramyxovirus B16-40 strain was closely related to the bat paramyxovirus Ms-ParaV/Anhui2011 (KC154054) and BtMI-ParaV/QH2013 (KJ641657) strains rather than Beilong virus, J-virus, and Tailam virus (Figs 2B and S1). Even though the HN amino acids sequences were similar to those from viruses in the proposed genera Shaanvirus, it was also related to that of Sendai virus and human parainfluenza virus 1, which belong to a different genus, Respirovirus showing (Table 1) .

Based on the partial RNA-dependent RNA polymerase (RdRp) gene, the bat paramyxovirus isolates and other Korean bat paramyxoviruses belong to the genus Shaanvirus, with 68.8-79.0% nucleotide similarities (Fig. 2C ). This suggests a single, closely related group of paramyxoviruses circulating in Korea. Consistent with the genetic findings, the bat paramyxovirus B16-40 was found to have haemagglutinin and neuraminidase activities in the haemadsorption assay (0.5% chicken red blood cells in PBS (pH 7.4)) and neuraminidase assay, respectively (Fig. 3 ).

In the preliminary studies of mouse antisera (BALB/ cA mice) against human respiratory syncytial virus A (KBPV-VR-41), human mumps virus (KBPV-VR-51) and human parainfluenza virus 1 (KBPV-VR-44), the antisera against human parainfluenza virus 1 (KBPV-VR-44) was only cross-reactive with the bat paramyxovirus B16-40 as determined using the immunofluorescence assay. Therefore, further cross-reactivity and cross-neutralizing tests between human parainfluenza virus 1 (KBPV-VR-44) and the bat paramyxovirus B16-40 were performed.

As shown in Table 2 , partial cross-reactivities were observed between the two viruses in the indirect immunofluorescence assay, even though they belong to different genera. Among the six mouse sera against the bat paramyxovirus B16-40, two sera were cross-reactive to human parainfluenza virus 1 (KBPV-VR-44) with 10 and 40 as the end-point titers for the fluorescent signals. In addition, among three pooled sera (two mouse sera each) against human parainfluenza virus 1 (KBPV-VR-44), one pooled serum was cross-reactive to the bat paramyxovirus B16-40 with 40 as the end-point titer for the fluorescent signal ( Table 2 , Fig. 4 ).

While the mouse antisera neutralized the homologous virus with 20-160 serum neutralizing titer, they did not neutralize the heterologous virus (Table 2 ).

The bat paramyxovirus B16-40 in this study was similar to viruses in the recently proposed genera Shaanvirus and Jeilongvirus. The genus Jeilongvirus includes Beilong virus, which was first isolated from human kidney mesangial cells 21 , and J virus, first isolated from mice 22 . These two viruses were only recently proposed to form Jeilongvirus as a distinct genus 14 . In addition, the newly discovered bat paramyxoviruses, Ms-ParaV/ Anhui2011 and BtMl-ParaV/QH2013, were recently proposed to form the additional genus Shaanvirus due to their genomic length and sequence identity differences 13 . However, their taxonomy has not yet been confirmed by the International Committee on Taxonomy of Viruses (ICTV, https://talk.ictvonline.org). The proposed genera Shaanvirus and Jeilongvirus have two novel genes putatively encoding SH (small hydrophobic) and TM (transmembrane) proteins 12, 13 . Two viral genes were also located between the fusion and haemagglutinin-neuraminidase genes of bat paramyxovirus B16-40. Further, their nucleotide and amino acids Reference viruses belonging to Paramyxoviridae were included and aligned by MAFFT. The phylogenetic tree was generated by the maximum-likelihood method with 1,000 replicates of bootstrap sampling and the Kimura 2-parameter model using MEGA 6. Bat paramyxovirus B16-40 isolated in this study is indicated with the blue box. (C) Phylogenetic analysis based on partial RNA-dependent RNA polymerase (RdRp) nucleotide sequences of bat paramyxoviruses detected in the Korean bats and other reference viruses belonging to Paramyxoviridae. The phylogenetic tree was generated by the maximum-likelihood method with 1,000 replicates of bootstrap sampling and the Kimura 2-parameter model using MEGA 6. Korean bat paramyxoviruses are indicated by blue boxes.

sequences were similar to the SH and TM genes identified from bat paramyxovirus strain Ms-ParaV/Anhui 2011. Therefore, the isolated bat paramyxovirus B16-40 might belong to the proposed genus Shaanvirus rather than Jeilongvirus.

The bat paramyxovirus B16-40 in this study was predicted to have a haemagglutinin-neuraminidase (HN) protein as an attachment glycoprotein in the sequence analysis. Consistent with the genetic findings, the bat paramyxovirus B16-40 was found to have haemagglutinin and neuraminidase activities in the haemadsorption assay (0.5% chicken red blood cells in PBS (pH 7.4)) and neuraminidase assay, respectively ( Fig. 2A,B) . Notably, even though the HN amino acids sequences were similar to those from viruses in the proposed genera Shaanvirus, it was also related to that of Sendai virus and human parainfluenza virus 1, which belong to a different genus, Respirovirus (Table 1) . Maximum-likelihood tree analysis based on the full HN amino acid sequences also showed that the bat paramyxovirus B16-40 was related to viruses in the genus Respirovirus ( Supplementary  Fig. S1 ). Therefore, the B16-40 HN amino acid sequence was compared to that of human parainfluenza virus 1 (KBPV-VR-44), and six regions that have at least six conserved amino acids were found. These regions were located on the sialidase superfamily domain and were structurally gathered together ( Supplementary Fig. S2 ). Among the seven neuraminidase active sites (R, D, E, R, R, Y, and E 23 ), two of these (E at the third site and Y) were located within the conserved amino acid sequences. These bioinformatics findings might indicate that the HN proteins of bat paramyxovirus B16-40 and human parainfluenza virus 1 (KBPV-VR-44) are genetically related, even though they are in different genera. Further study on their antigenic and evolutionary relationship should be followed up.

Notably, in a previous study, bat serum (diluted 1/40) from African wild bat, Micropteropus pusillus was reactive with human parainfluenza virus 1 in an indirect immunofluorescence assay 3 . There may be two possibilities to explain this. First, the African bat might have been infected with a human parainfluenza virus 1-related virus in the genus Respirovirus that has not yet been found. However, as there has been no evidence of bat paramyxoviruses related to the genus Respirovirus in recent studies 13, 15, 16 , the first possibility seems low. Second, the bat might have been infected by other bat paramyxoviruses which were cross-reactive with human parainfluenza virus 1.

In fact, in this study, when mouse antisera were made and tested against bat paramyxovirus B16-40 and human parainfluenza virus 1 (KBPV-VR-44), the two viruses were partially cross-reactive to each other in an indirect immunofluorescence assay. This result provided an evidence that there could be cross-reactivity between viruses belonging to different genera. Therefore, the bat paramyxovirus B16-40-related viruses circulating in bats might be one of the reasons why the wild bat sera was reactive with human parainfluenza virus 1 in the previous study. However, we could not find the cross-neutralization between bat paramyxovirus B16-40 and human parainfluenza virus 1. Based on the recent studies, cross-neutralization was found between the viruses in the same genus of paramyxoviruses: Sendai virus and human parainfluenza virus 1, which both belong to the genus Respirovirus 24 , and human mumps virus and African bat mumps virus, which both belong to the genus Rubulavirus 25 . Further, the engaged proteins in the cross-reactivity between bat paramyxovirus B16-40 and human parainfluenza virus 1 had not been revealed. Therefore, although bat paramyxovirus B16-40 in this study was cross-reactive with human parainfluenza virus 1 in indirect immunofluorescence assay, it is difficult to say that the bat paramyxovirus B16-40 was antigenically related with human parainfluenza virus 1.

Since a number of paramyxoviruses have been newly found around world, their classification has been difficult due to the limited criteria, such as sequence alone 18 . Additional input from the biological context may be helpful. In this study, a bat paramyxovirus B16-40 belonging to the proposed genus Shaanvirus was first isolated. Through virus isolation, we could obtain not only nucleotide sequence information, but also several biological characteristics of the virus. Its putative HN protein was demonstrated to have haemagglutinin and neuraminidase activities. Mouse antisera against the bat paramyxovirus B16-40 was also cross-reactive to human parainfluenza virus 1 in the indirect immunofluorescence assay, but could not cross-neutralize human parainfluenza virus 1. Additionally, the bat paramyxovirus B16-40 was not infectious in 6-week-old female C57BL/6 mice. Therefore, this study provided basic information on further classification of the bat paramyxovirus B16-40 and related viruses in the proposed genus Shaanvirus. Samples. In 2016, a total of 232 samples were collected at 48 sites in natural bat habitats. Most samples were collected as fresh guano under the colony of bats. In some cases, bats were captured by a fine mesh net to collect fresh feces as well as oropharyngeal and urine swabs, and the animals were then immediately released. The collected samples were placed into transport medium (Noble Bio. Co. Ltd) in a 10% suspension and were transported to the lab (BSL2) for further analysis. The major bat species at the collection sites were determined based on morphology and on previous data from bats' roosting sites 26 Table S1 ).

RNA extracted from bat samples was tested by RT-semi-nested PCR using consensus paramyxovirus primers targeting the RdRp region 19 . The nucleotide sequences of the positive samples were determined by Cosmogenetech Co. Ltd. Five of the 232 bat samples were confirmed positive to the bat paramyxoviruses by sequencing. Virus isolation was performed from the RT-semi-nested PCR-positive samples including feces, urine, and oropharyngeal swabs from bats. Each bat sample was diluted 10-fold with fresh Dulbecco Modified Eagle Medium (DMEM) and filtered with a 0.2-μm filter before being inoculated into the MARC 145 cell line, a clone from the MA-104 cell line from kidney epithelial cells of an African green monkey 20 ; blind passaging was performed up to three times. The cytopathic effect (CPE) was confirmed in cells infected with one sample (B16-40) of five PCR-positive samples. Next, a monolayer of MARC 145 was infected with bat paramyxovirus (B16-40) at a multiplicity of infection (MOI) of 0.1. Two days post-infection, the infected cells were harvested and pelleted. The infected cell pellet was fixed with 2.5% glutaraldehyde in PBS (pH 7.4) for 2 hours and sent for transmission electron microscopy (TEM) to the Advanced Analysis Center in the Korean Institute of Science and Technology in Korea. Briefly, after washing with PBS the cells were post-fixed in 2% osmium tetroxide and dehydrated in an ethanol series (70%, 80%, 85%, 90%, 95%, and 100% (two times)). Then, the sample was transited with 100% propylene oxide and embedded in Epoxy resin. The embedded sample was sectioned into 70-80 nm sections with an ultramicrotome (Ultra Cut C; Leica), stained with 2% uranyl acetate and lead citrate on a copper grid, and examined using TEM mode on a CryoTecnai F20 (FEI Ltd., USA). culture supernatant and freeze-thawed cell lysate were centrifuged at 3,000 � g for 10 minutes to remove cell pellets and then used for viral titration.

The infectious viruses were quantified from the prepared samples using the tissue culture infectious dose 50(TCID 50 ) method. Briefly, the prepared samples were 10-fold serially diluted and inoculated to monolayers of MARC-145 cells in 4 wells. The highest dilution showing 50% of viral infection was obtained by observing cytopathic effects with a light inverted microscope (IX71S1F-3, OLYMPUS CORPORATION). The TCID 50 was calculated using the Spearman-Karber method.

Genomic sequencing and analyses. A virus isolate showing cytopathic effects (CPEs) was filtered through a 0.2 μm filter and ultra-centrifuged at 250,000 � g for 1 hour. Each pellet was suspended in 500 µL of 1x digestion buffer (Turbo DNA Free Kit; Ambion) and treated with Turbo DNase; the suspension was incubated at 37 °C for 30 min. The RNA from the suspension was isolated using the Trizol LS reagent according to the manufacturer's instructions (Ambion). The extracted RNA was submitted to Macrogen (Seoul, Korea) for high throughput sequencing in a HiSeq. 2000 sequencing system based on the transcriptome de novo sequencing platform. The obtained viral contigs were analysed and annotated by the MG-RAST server 27 . The cut-off for the annotation was 10 �5 maximum e-value, 60% minimum percentage identity, and 15 for minimum alignment length. The viral genomic sequences were obtained by RT-PCR using primers based on the Paramyxoviridae-associated contigs obtained by high throughput sequencing. Additionally, 3� and 5� end sequences were obtained by RACE sequencing based on PCR using the adapter-oligo dTVN primer from Bionics Co. Ltd.

For human parainfluenza virus 1 (KBPV-VR-44), the 1,992 bases of the HN gene were obtained by RT-nested PCR following the previously published method 28 . The obtained genomic RNA sequences and RT-PCR-positive amplicons were further analysed with related sequences in GenBank using MAFFT 29 and MEGA version 6 30 .

Genomic sequence data generated in this study have been deposited in GenBank under accession number MG230624. The full HN sequence of human parainfluenza virus 1 (KBPV-VR-44) has also been deposited in GenBank under accession number MG255129. The HN amino acid sequences of human parainfluenza virus 1 (KBPV-VR-44) and bat paramyxovirus were further compared with the BioEdit program 31 . Putative HN protein structures were predicted using the Phyre2 website 32 and visualized with the UCSF Chimera 1.11.2 program 33 .

Haemadsorption and neuraminidase assays. To investigate the activity of the HN protein, the haemadsorption assay was used. Monolayer cultures of MARC 145 cells in 96-well plates were inoculated with the bat paramyxovirus B16-40 at an MOI of 0.1. After 48 hours, the cells were gently washed with PBS and incubated for 1 hour with 0.5% chicken red blood cells (RBCs) in PBS (pH 7.4) at room temperature. Then, the cells were washed with PBS twice to remove unbound RBCs and observed under a microscope. In addition, the neuraminidase assay was performed with the viral isolate using the NA-Fluor Influenza Neuraminidase Assay Kit (Thermo) according to the manufacturer's protocol. The virus was diluted 2-fold with the NA-Fluor 2� assay buffer in a black 96-well plate and incubated for 1 hour at 37 °C with the NA-Fluor substrate. After adding the stop solution, the plate was read at the excitation and emission wavelengths of 355 nm and 460 nm, respectively, using a PerkinElmer VICTOR2 fluorometer.

A total of 44 6-week-old female C57BL/6 mice maintained under specific-pathogen-free (SPF) conditions were divided into two groups of negative controls (n = 4 each) and four inoculation groups (n = 9 each). The inoculation groups were inoculated with 400 µL of 10 5 TCID 50 or 10 2.5 TCID 50 of the bat paramyxovirus B16-40 via intra-gastric administration routes, and 30 µl of 10 5 TCID 50 or 10 2.5 TCID 50 of the bat paramyxovirus B16-40 via the intranasal administration routes. Each inoculation was conducted twice with 1 week in between. Weight loss was monitored daily for up to 2 weeks. The lungs, liver, brain, and intestine were harvested at 7 or 14 days post-challenge, fixed with 4% paraformaldehyde in PBS, pH = 7.2, and submitted to BioLead Inc., Korea, for histopathological examination. Additionally, RNA was extracted from each tissue and from oral swabs and fecal samples collected daily for 7 days post-challenge. The extracted RNA was tested using RT-semi-nested PCR with the consensus paramyxovirus primers 19 . In addition, real time PCR was performed with newly designed primers and probes targeting regions of the membrane and fusion protein. Briefly, the designed primers were as follows:

forward primer: PVM-F5�-CCCAGGAGTATGGTTATCAAGTGAGG-3�; reverse primer: PVM-R 5�-TCCA TTGGGCTCTCTTTGTTTGC-3�; Taqman probe: PVM-P 5�-FAM-CCCATCCCAGACCAGCCACCA GACCC-TAMRA-3�

Real time RT-PCR was performed as follows: Reverse transcription at 45 °C for 10 min, followed by PCR �95 °C for 5 min, cycling 40 times at 95 °C for 10 sec, and 60 °C for 30 sec. Using the 10-fold diluted virus (bat paramyxovirus B16-40), standard curves were performed for every reaction.

Generation of mouse antisera against bat paramyxovirus B16-40 and human paramyxoviruses. We next observed preliminary cross-reactivity in indirect immunofluorescence assay between several human viruses belonging to the family Paramyxoviridae and the bat paramyxovirus B16-40. 6-week-old female BALB/cA mice maintained under SPF conditions were intramuscularly inoculated twice, with 2 weeks in between, with a 1:1 ratio mixture of AddaVax adjuvant and live virus stock (bat paramyxovirus B16-40, human respiratory syncytial virus A (KBPV-VR-41), human mumps virus (KBPV-VR-51), and human parainfluenza virus 1 (KBPV-VR-44), respectively, from the Korea Bank for Pathogenic Viruses. Two weeks after the last immunization, antisera against the virus were obtained by cardiac puncture for further studies related to cross-reactivity in indirect immunofluorescence assay. Mouse antisera against the bat paramyxovirus B16-40 (10 5.75 TCID 50 /mL) and human parainfluenza virus 1 (KBPV-VR-44, 10 5.75 TCID 50 /mL) were further compared through serum titration using the indirect immunofluorescence assay and serum neutralization test.

The washed cells were fixed with 100% ethanol for 1 hour at �20 °C and rehydrated with PBS (pH 7.4) for 10 min at room temperature. Primary antibodies used in the experiments were the mouse antisera against bat paramyxovirus B16-40 and human parainfluenza virus 1 (KBPV-VR-44). The sera were diluted with 50% culture supernatants from the cell lines in fresh DMEM; the corresponding mock-infected cells were treated in the same way as the target virus culture. Then, the sera were serially diluted 2-fold from an initial dilution of 1:10 and added to the fixed infected or mock-infected cells and incubated for 1 hour at room temperature. The plates were washed three times with PBS containing 0

with the same serum diluent and was added to the plates and incubated for 1 hour at room temperature. After three washes with PBS-T, the plates were viewed with a fluorescence microscope

They were then diluted 2-fold from an initial dilution of 1:10 in 96-well plates and mixed with 200 TCID 50 of each virus. The plates were then incubated at 37 °C for 1 hour. The virus-serum mixtures for the bat paramyxovirus B16-40 and human parainfluenza virus 1 (KBPV-VR-44) were added to monolayer cultures of MARC 145 and LLC-MK2 cells, respectively, and incubated at 37 °C for 1 hour. The MARC 145 cells and LLC-MK2 cells were then washed once with PBS

Bats as reservoirs of severe emerging infectious diseases

Detection of Severe Acute Respiratory Syndrome-Like, Middle East Respiratory Syndrome-Like Bat Coronaviruses and Group H Rotavirus in Faeces of Korean Bats

Bats host major mammalian paramyxoviruses

Taxonomy of the order Mononegavirales: update 2017

Hendra and Nipah infection: emerging paramyxoviruses

Isolation of multiple novel paramyxoviruses from pteropid bat urine

Novel paramyxoviruses in free-ranging European bats

Novel paramyxoviruses in bats from sub-Saharan Africa

Identification of novel paramyxoviruses in insectivorous bats of the Southwest Indian Ocean

Nipah virus transmission from bats to humans associated with drinking traditional liquor made from date palm sap

Serological evidence of possible human infection with Tioman virus, a newly described paramyxovirus of bat origin

The complete genome sequence of J virus reveals a unique genome structure in the family Paramyxoviridae

Deciphering the bat virome catalog to better understand the ecological diversity of bat viruses and the bat origin of emerging infectious diseases

The immune evasion function of J and Beilong virus V proteins is distinct from that of other paramyxoviruses, consistent with their inclusion in the proposed genus Jeilongvirus

Coronavirus and paramyxovirus in bats from Northwest Italy

Circulation of Alphacoronavirus, Betacoronavirus and Paramyxovirus in Hipposideros bat species in Zimbabwe. Infection

Novel Alphacoronaviruses and paramyxoviruses Cocirculate with type 1 and severe acute respiratory system (SARS)-related Betacoronaviruses in Synanthropic bats of Luxembourg

Problems of classification in the family Paramyxoviridae

Sensitive and broadly reactive reverse transcription-PCR assays to detect novel paramyxoviruses

Enhanced replication of porcine reproductive and respiratory syndrome (PRRS) virus in a homogeneous subpopulation of MA-104 cell line

Beilong virus, a novel paramyxovirus with the largest genome of non-segmented negative-stranded RNA viruses

A new mouse paramyxovirus (J virus)

Structural studies of the parainfluenza virus 5 hemagglutinin-neuraminidase tetramer in complex with its receptor, sialyllactose

Safety and immunogenicity of an intranasal Sendai virus-based human parainfluenza virus type 1 vaccine in 3-to 6-year-old children

Cross-neutralization between human and African bat mumps viruses

Sounds of the Bats in Korea

The metagenomics RAST server-a public resource for the automatic phylogenetic and functional analysis of metagenomes

A study of genetic variability of human parainfluenza virus type 1 in Croatia

MAFFT multiple sequence alignment software version 7: improvements in performance and usability

MEGA6: molecular evolutionary genetics analysis version 6.0. Molecular biology and evolution

In Nucleic acids symposium series

The Phyre2 web portal for protein modeling, prediction and analysis

UCSF Chimera-a visualization system for exploratory research and analysis

Data availability. The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request. 

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-30319-7.Competing Interests: The authors declare no competing interests.Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.