key: cord-0821235-0duo77v8
authors: Lee, Seung-Ho; Kim, Kijin; Kim, Jongwoo; No, Jin Sun; Park, Kyungmin; Budhathoki, Shailesh; Lee, Seung Ho; Lee, Jingyeong; Cho, Seung Hye; Cho, Seungchan; Lee, Geum-Young; Hwang, Jusun; Kim, Heung-Chul; Klein, Terry A.; Uhm, Chang-Sub; Kim, Won-Keun; Song, Jin-Won
title: Discovery and Genetic Characterization of Novel Paramyxoviruses Related to the Genus Henipavirus in Crocidura Species in the Republic of Korea
date: 2021-10-07
journal: Viruses
DOI: 10.3390/v13102020
sha: 0a29b8c1bd1874dc6458e15a60e233d7eebcef69
doc_id: 821235
cord_uid: 0duo77v8

Paramyxoviruses, negative-sense single-stranded RNA viruses, pose a critical threat to human public health. Currently, 78 species, 17 genera, and 4 subfamilies of paramyxoviruses are harbored by multiple natural reservoirs, including rodents, bats, birds, reptiles, and fish. Henipaviruses are critical zoonotic pathogens that cause severe acute respiratory distress and neurological diseases in humans. Using reverse transcription-polymerase chain reaction, 115 Crocidura species individuals were examined for the prevalence of paramyxovirus infections. Paramyxovirus RNA was observed in 26 (22.6%) shrews collected at five trapping sites, Republic of Korea. Herein, we report two genetically distinct novel paramyxoviruses (genus: Henipavirus): Gamak virus (GAKV) and Daeryong virus (DARV) isolated from C. lasiura and C. shantungensis, respectively. Two GAKVs and one DARV were nearly completely sequenced using next-generation sequencing. GAKV and DARV contain six genes (3′-N-P-M-F-G-L-5′) with genome sizes of 18,460 nucleotides and 19,471 nucleotides, respectively. The phylogenetic inference demonstrated that GAKV and DARV form independent genetic lineages of Henipavirus in Crocidura species. GAKV-infected human lung epithelial cells elicited the induction of type I/III interferons, interferon-stimulated genes, and proinflammatory cytokines. In conclusion, this study contributes further understandings of the molecular prevalence, genetic characteristics and diversity, and zoonotic potential of novel paramyxoviruses in shrews.

and 'Daeryong virus' from the first discovered location of origin (Mount Gamak and Daeryong, respectively). GAKV, DARV, and rodent-borne henipavirus (MojV) appear to be distinct from previously identified bat-borne henipaviruses (Hendra virus, NiV, Cedar virus, and Kumasi virus), suggesting that they may be novel species within the genus Henipavirus. Infection with GAKV elicited the expression of type I/III interferons (IFNs), IFN-stimulated genes (ISGs), and proinflammatory cytokines in A549 cells. Therefore, this report demonstrates the molecular prevalence, genetic characteristics and diversity, and virus-host interactions of paramyxoviruses in shrews, ROK.

The US Forces Korea (USFK) approved the animal trapping procedure by USFK Regulation 40-1 "Prevention, Surveillance, and Treatment of Hemorrhagic Fever with Renal Syndrome". The trapping, experiments, and handling of animals were conducted under a protocol approved by the Korea University Institutional Animal Care and Use Committee (KUIACUC, #2016-0049).

Small mammals were captured during 2017-2018 using Sherman traps (8 cm × 9 cm × 23 cm; H. B. Sherman, Tallahassee, FL, USA). A total of 115 individuals belonging to Crocidura species were collected from the previous study [30] . Shrews were identified using morphological criteria and polymerase chain reaction (PCR), when required. Serum, brain, lung, spleen, kidney, and liver tissues were collected aseptically and frozen at −80 • C until use. Total RNA was extracted from the kidney tissues of shrews using TRI Reagent Solution (AMBION Inc., Austin, TX, USA) according to the manufacturer's instructions.

cDNA was generated from the RNA extracted from shrew-borne paramyxovirusinfected cells and kidney tissues using SISPA method was carried out according to a protocol described previously [30] .

The cDNA library was prepared using the TruSeq Nano DNA LT Sample Preparation Kit (Illumina, San Diego, CA, USA), and the amplicon was size-selected, A-tailed, ligated with indexes and adaptors, and enriched. The library was verified for quality and size with Agilent 2100 Bioanalyzer (Agilent High Sensitivity DNA kit, Agilent Technologies) and sequenced with MiSeq reagent V2 (Illumina) using MiSeq benchtop sequencer (Illumina) with 2 × 150 bp.

Libraries were prepared by using the NEBNext Ultra II Directional RNA-Seq Kit (New England Biolabs, Ipswich, MA, USA). Additionally, mRNA was isolated using the Poly(A) RNA Selection Kit (Lexogen, Vienna, Austria). The mRNAs were synthesized for cDNA and sheared according to the manufacturer's instructions. The fragments were size-selected, A-tailed, ligated with indexes and adaptors, and enriched by PCR. The libraries were evaluated for the mean fragment size by the Agilent 2100 bioanalyzer (DNA High Sensitivity Kit). Quantification was performed by using a library quantification kit and StepOne Real-Time PCR System (Life Technologies). Libraries were conducted as paired-end 100 sequencing by the HiSeq X10 system (Illumina). 

The reads were trimmed with Trimmomatic (v.0.36) to remove adapter sequences for the metagenomic approach [44] . Aligned reads against the host genomic sequences were excluded using Bowtie2 (v.2.2.6) [45] . The complete mitochondrial sequence of the species on the National Center for Biotechnology Information (NCBI) RefSeq was utilized as a host reference [46] . The remaining reads were filtered for quality by FaQCs (v0.11.5), and de novo assembly was conducted by using SPAdes (v3.11.1) [47, 48] . The assembled contigs were subsequently analyzed in the NCBI RefSeq database consisting of complete viral genomes (updated in June 2019) by nucleotide-Basic Local Alignment Search Tool (BLASTn) (v2.6.0). In addition, for the reference mapping strategy, adaptor and index sequences of reads were trimmed, and low-quality sequences were filtered using the CLC Genomics Workbench v.7.5.2 (CLC Bio, Cambridge, MA, USA). The genome sequences of Hendra virus (HeV), NiV, and MojV were used in the reference mapping method. The consensus genomic sequences of GAKV and DARV were determined by combining viral contigs extracted from CLC analysis and de novo assembly.

RNA was used for cDNA synthesis using a high-capacity RNA-to-cDNA kit (Applied Biosystems, Foster City, CA, USA). First and nested PCRs were conducted in a 25 µL reaction mixture containing 2.5 U of Ex Taq DNA polymerase (TaKaRa BIO Inc., Shiga, Japan), 1.5 µg of cDNA, and 10 pM of pan-paramyxoviral primers [49] . PCR condition and purification were carried out according to a protocol described previously [30] . Sequencing was performed in both directions of each PCR product using a BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) on an automated sequencer (ABI 3730XL DNA Analyzer; Applied Biosystems).

The 3 and 5 terminal genome sequences of paramyxovirus were determined by using a SMARTer ® RACE 5 /3 Kit (Takara Bio), according to the manufacturer's specifications. PCR products were purified by the LaboPass PCR Purification Kit (Cosmo Genetech). Sequencing was conducted in both directions of each PCR product using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) on an automated sequencer (Applied Biosystems).

DNA was extracted from liver or lung tissues using a High Pure PCR template preparation kit (Roche, Basel, Switzerland). The cytochrome b gene of mtDNA was identified by universal primers [50] .

The viral genomic sequences were aligned with MUSCLE algorithm in MEGA v.7.0 [51] and trimmed using the Lasergene program v.5 (DNASTAR, Madison, WI, USA). Phylogenetic analyses were conducted by the maximum likelihood method according to the best-fit substitution model (GTR+G+I) using MEGA 7.0. Support for the topologies was evaluated by bootstrapping for 1000 iterations. Additionally, the Markov chain Monte Carlo (MCMC) method and the BEAST package (v.1.10.4) as the Bayesian inference method were implemented to infer Bayesian phylogenetic trees [52] . The MCMC chain length was established to 100 million states by sampling every 50,000 states. The results of parameters produced sufficient sample sizes (ESS > 200). Maximum clade credibility trees were generated using TreeAnnotator (v.1.10.4) and FigTree (v.1.4.0).

To compare different phylogenetic links matching between paramyxoviruses and their hosts, a tanglegram algorithm was performed. The auxiliary lines in the center connect the phylogenetic trees. The complexity between dendrograms of phylogenies was diminished to a maximum before complete reconciliation analysis. The method was implemented using the R package "dendextend" [53] .

Vero E6 cells (ATCC, #DR-L2785) and human lung adenocarcinoma cells (A549) (ATCC, #CCL-185) were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA). Vero E6 and A549 cells were grown in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal bovine serum (FBS), 1 mM sodium pyruvate, 2 mL of L-glutamine, and 50 mg/mL gentamicin at 37 • C in a 5% carbon dioxide (CO 2 ) incubator.

Tissue homogenates of kidneys from GAKV-positive samples were inoculated on Vero E6 cells. The inoculum was removed after one and a half hours of adsorption. The culture was maintained with 5.5 mL of DMEM containing 5% FBS at 37 • C in a 5% CO 2 incubator. Cytopathic effects were inspected daily using an inverted microscope.

Vero E6 cells were prepared at a density of 1.5 × 10 6 cells per well onto 6-well plates. The cells were washed twice with phosphate-buffered saline (PBS) and inoculated serially with 10-fold diluted viruses. A 1:1 medium-melting-point agarose and overlay medium mix was added on the cells after 90 m of incubation at 37 • C with constant shaking. The agarose overlay was discarded at day 6. The plaques were stained with 0.1% crystal violet in 10% formaldehyde for visualization.

GAKV-inoculated Vero E6 cells were collected at 5 d postinfection and fixed with 2% paraformaldehyde and 2.5% glutaraldehyde in 0.1 M PBS (pH 7.4). Ultrathin sections were placed on 400-mesh square copper electron microscopy grids (Electron Microscopy Sciences) and observed by using a transmission electron microscope (Model H-7650; Hitachi, Japan).

A total of 1 × 10 6 cells were infected with GAKV at a multiplicity of infection of 0.02. The samples were then collected at 2, 6, 12, 18, and 24 h postinfection (hpi). The mRNA genes of type I/III IFNs, ISGs, and proinflammatory cytokines were examined by using the primers (Table S1 ).

Statistical analyses were conducted as shown in figures using the GraphPad Prism v.5.00 for Windows (GraphPad Software, San Diego, CA, USA; www.graphpad.com (accessed on 1 September 2021)).

A total of 115 shrews (94 C. lasiura and 21 C. shantungensis) were captured in 11 regions of the ROK from 2017 to 2018. All samples were obtained from Gangwon (Cheorwon, Chuncheon, Hwacheon, Hongcheon, and Pyeongchang), Gyeonggi (Osan, Paju, Pocheon, Suwon, and Yeoncheon), and Gyeongsangnam (Changnyeong) Provinces ( Figure 1 ). To detect paramyxovirus RNA, we performed RT-PCR using consensus primers targeting the partial sequence of L gene of paramyxoviruses. In total, 20 (21.3%) C. lasiura and six (28.6%) C. shantungensis were positive for Gamak and Daeryong henipaviruses, respectively (Table 1 Viruses 2021, 13, 2020 6 of 16 and Table S2 ). GAKV was discovered in animals from five collection sites (Paju, Yeoncheon in Gyeonggi Province, Cheorwon, Chuncheon in Gangwon Province, and Changnyeong in Gyeongsangnam Province), whereas DARV was detected in animals from one collection site (Chuncheon in Gangwon Province) (Tables S3 and S4 ). The geographic prevalence of GAKV followed 11/40 (27.5%) in Gangwon Province, 9/52 (17.3%) in Gyeonggi Province, and 1/2 (50.0%) in Gyeongsangnam Province, while DARV was detected only in Gangwon Province. A total of 16 of 58 (27.6%) males and 5 of 36 (13.9%) females were detected for GAKV RNA. Subadult (<10.0 g) and adult (10-20 g) C. lasiura harbored GAKV infection with 11/54 (20.4%) and 10/40 (25.0%), respectively. GAKV-positive shrews were found in spring, summer, and autumn. A total of one of nine males (11.1%) and four of 12 females (33.3%) infected DARV. Three subadult C. shatungensis (<5.0 g) were positive for DARV, whereas two adult C. shatungensis (≥5.0 g) were infected with the virus. DARV-infected shrews were found only in autumn.

Paju, Pocheon, Suwon, and Yeoncheon), and Gyeongsangnam (Changnyeong) Provinces ( Figure 1 ). To detect paramyxovirus RNA, we performed RT-PCR using consensus primers targeting the partial sequence of L gene of paramyxoviruses. In total, 20 (21.3%) C. lasiura and six (28.6%) C. shantungensis were positive for Gamak and Daeryong henipaviruses, respectively (Table 1 and Table S2 ). GAKV was discovered in animals from five collection sites (Paju, Yeoncheon in Gyeonggi Province, Cheorwon, Chuncheon in Gangwon Province, and Changnyeong in Gyeongsangnam Province), whereas DARV was detected in animals from one collection site (Chuncheon in Gangwon Province) (Table S3 and S4). The geographic prevalence of GAKV followed 11/40 (27.5%) in Gangwon Province, 9/52 (17.3%) in Gyeonggi Province, and 1/2 (50.0%) in Gyeongsangnam Province, while DARV was detected only in Gangwon Province. A total of 16 of 58 (27.6%) males and 5 of 36 (13.9%) females were detected for GAKV RNA. Subadult (˂10.0 g) and adult (10-20 g) C. lasiura harbored GAKV infection with 11/54 (20.4%) and 10/40 (25.0%), respectively. GAKV-positive shrews were found in spring, summer, and autumn. A total of one of nine males (11.1%) and four of 12 females (33.3%) infected DARV. Three subadult C. shatungensis (˂5.0 g) were positive for DARV, whereas two adult C. shatungensis (≥5.0 g) were infected with the virus. DARV-infected shrews were found only in autumn. 

GAKV was isolated from the kidney tissues of Cl17-32 using a cell culture-based method. The first isolate of GAKV was confirmed by passaging twice for 14 days postinoculation. GAKV particles were observed using a transmission electron microscope ( Figure 2A ). In addition, the number of infectious GAKV particles was 1 × 10 6 PFU/mL, as quantified using the plaque assay ( Figure 2B ).

We performed host rRNA depletion on SISPA-based MiSeq of the kidney tissue of Cl17-46, resulting in two contigs (602 and 17,892 nt in length) with significant similarities in the genomic sequence of paramyxoviruses. The NGS of Cl17-46 generated 245,846 reads, and the number of viral genome sequences was 3608 (Table S5) . To obtain whole-genome sequences of shrew-borne paramyxovirus, RNA-Seq (HiSeq) was performed on samples from the kidney tissues of Cl17-46 and Cs17-65 and on cells infected with the paramyxovirus isolated from Cl17-32. The genomic sequences of henipaviruses were identified across 23 contigs (408-18,382 nt in length) with significant similarities. The number of viral genome reads was 14,779,024, 6497, and 50,987 from Cl17-32, Cl17-46, and Cs17-65, respectively (Table S6) . Nearly complete genome sequences of GAKV and DARV were obtained except for both 3 and 5 terminal ends. The genomic sequences of GAKV and DARV have been deposited in GenBank (accession number: MZ574407-MZ574409). 

We performed host rRNA depletion on SISPA-based MiSeq of the kidney tissue of Cl17-46, resulting in two contigs (602 and 17,892 nt in length) with significant similarities in the genomic sequence of paramyxoviruses. The NGS of Cl17-46 generated 245,846 reads, and the number of viral genome sequences was 3608 (Table S5) . To obtain wholegenome sequences of shrew-borne paramyxovirus, RNA-Seq (HiSeq) was performed on samples from the kidney tissues of Cl17-46 and Cs17-65 and on cells infected with the paramyxovirus isolated from Cl17-32. The genomic sequences of henipaviruses were identified across 23 contigs (408-18,382 nt in length) with significant similarities. The number of viral genome reads was 14,779,024, 6497, and 50,987 from Cl17-32, Cl17-46, and Cs17-65, respectively (Table S6) . Nearly complete genome sequences of GAKV and DARV were obtained except for both 3′ and 5′ terminal ends. The genomic sequences of GAKV and DARV have been deposited in GenBank (accession number: MZ574407-MZ574409). The N, M, F, G, and L genes encode a single protein, whereas the P gene encodes multiple accessory proteins: the viral phosphoprotein, C protein, and V/W protein. GAKV and DARV possessed a putative RNA editing site (TTAAAAAAGGCA) of the P gene, a conserved motif sequence (YTAAAARRGGCA) in the members of the genus Henipavirus. The 3′ leader and 5′ trailer sequences were a length of 52 nt and 25 nt, respectively. Tables S7 and S8 show the start, stop, and intergenic region (IGR) sequences of GAKV and DARV in shrews. 

The phylogenetic analysis of all 26 currently recognized Orthoparamyxovirus species elucidated that GAKV and DARV belong to the genus Henipavirus and are closely related to the henipaviruses (Figure 4 ). In addition, the amino acid sequences (N, P, M, F, G, and L) of GAKV and DARV showed the phylogenetic shape comparable to the viral RNA genome sequences among major paramyxoviruses ( Figure S1 ). The genetic cluster of 

The phylogenetic analysis of all 26 currently recognized Orthoparamyxovirus species elucidated that GAKV and DARV belong to the genus Henipavirus and are closely related to the henipaviruses (Figure 4 ). In addition, the amino acid sequences (N, P, M, F, G, and L) of GAKV and DARV showed the phylogenetic shape comparable to the viral RNA genome sequences among major paramyxoviruses ( Figure S1 ). The genetic cluster of GAKV and DARV shared a common ancestor with MojV with an amino acid, similarity to the L gene of 62.9% and 74.1%, respectively (Table S9) .

The partial L genomic sequences of GAKV and DARV strains (13,790-14,215 nt and 14,849-15,311 nt, respectively) belong to the genus Henipavirus, subfamily Orthoparamyxovirinae ( Figure S2) . Consistently, the phylogenies of GAKV and DARV were closely related to those of MojV, HeV, NiV, Cedar virus (CedV), and Kumasi virus (KV).

GAKV and DARV shared a common ancestor with MojV with an amino acid, similarity to the L gene of 62.9% and 74.1%, respectively (Table S9) .

The partial L genomic sequences of GAKV and DARV strains (13,790-14,215 nt and 14,849-15,311 nt, respectively) belong to the genus Henipavirus, subfamily Orthoparamyxovirinae ( Figure S2) . Consistently, the phylogenies of GAKV and DARV were closely related to those of MojV, HeV, NiV, Cedar virus (CedV), and Kumasi virus (KV). 

Cophylogeny mapping demonstrated the segregation of paramyxoviruses according to the subfamily of their reservoir hosts, using consensus trees based on the L protein amino acid sequences ( 

Cophylogeny mapping demonstrated the segregation of paramyxoviruses according to the subfamily of their reservoir hosts, using consensus trees based on the L protein amino acid sequences ( Figure 5 ). The phylogenetic positions of GAKV, DARV, and shrew paramyxovirus (Shrew PV) mirrored the phylogenetic relationships of the shrew species. In contrast, the phylogenetic positions of viruses from rodents such as Bank Vole virus (BaVV) and Pohorje Myodes paramyxovirus 1 (PMPV-1) from Myodes glareolus, Mount Mabu Lophuromys virus 1 (MMLV-1) and Mount Mabu Lophuromys virus 2 (MMLV-2) from Lophuromys machangui, and PAPV-1 and PAPV-2 from Apodemus agrarius, were not directly matched, as per the tanglegram between the virus and host species.

To examine the infectivity and expression of innate antiviral genes in human lung epithelial cells, GAKV infected A549 cells for 2, 6, 12, 18, and 24 h. The replication of GAKV occurred and sequentially increased at 2, 6, 12, 18, and 24 hpi ( Figure 6A ). GAKV infection upregulated the mRNA copy numbers of type I/III IFN and ISGs, including Ifn-β, Ifnl1/Il-29, IFN-stimulated gene product 15 (ISG-15), Ifit2/Isg54, and Ifit1/Isg56 at 12, 18, and 24 hpi (Figure 6B-F) . The augmented induction of innate antiviral genes, radical S-adenosyl methionine domain containing 2 (Rsad2/Viperin), 2 , 5 -oligoadenylate synthetase 1 (OAS1), Ddx58/Rig-I, and Ifih1/Mda5, was observed at 12, 18, and 24 hpi ( Figure 6G-J) . The Il-6 mRNA gene copies were rapidly increased at 2, 6, 12, 18, and 24 hpi ( Figure 6K ). In summary, these data demonstrated that GAKV infection elicited the induction of innate antiviral genes in human lung epithelial cells. 

To examine the infectivity and expression of innate antiviral genes in human lung epithelial cells, GAKV infected A549 cells for 2, 6, 12, 18, and 24 h. The replication of GAKV occurred and sequentially increased at 2, 6, 12, 18, and 24 hpi ( Figure 6A ). GAKV infection upregulated the mRNA copy numbers of type I/III IFN and ISGs, including Ifnβ, Ifnl1/Il-29, IFN-stimulated gene product 15 (ISG-15), Ifit2/Isg54, and Ifit1/Isg56 at 12, 18, and 24 hpi ( Figure 6B-F) . The augmented induction of innate antiviral genes, radical S-adenosyl methionine domain containing 2 (Rsad2/Viperin), 2′, 5′-oligoadenylate synthetase 1 (OAS1), Ddx58/Rig-I, and Ifih1/Mda5, was observed at 12, 18, and 24 hpi ( Figure 6G-J) . The Il-6 mRNA gene copies were rapidly increased at 2, 6, 12, 18, and 24 hpi ( Figure 6K ). In summary, these data demonstrated that GAKV infection elicited the induction of innate antiviral genes in human lung epithelial cells. 

We characterized two novel paramyxoviruses, GAKV and DARV, identified in Crocidura species, ROK. GAKV and DARV showed a high positive rate: GAKV was found in five regions, while DARV was detected in one region. The genomic structure (3′-N-P-M-F-G-L-5′) of these viruses is representative of the genus of henipaviruses, including HeV, NiV, MojV, CedV, and KV. The phylogenies of GAKV and DARV shaped independent 

We characterized two novel paramyxoviruses, GAKV and DARV, identified in Crocidura species, ROK. GAKV and DARV showed a high positive rate: GAKV was found in five regions, while DARV was detected in one region. The genomic structure (3 -N-P-M-F-G-L-5 ) of these viruses is representative of the genus of henipaviruses, including HeV, NiV, MojV, CedV, and KV. The phylogenies of GAKV and DARV shaped independent genetic lineages corresponding to the paramyxovirus species distinctive criterion (the L gene amino acid distance of <0.82 for the Orthoparamyxovirinae) [26] . These observations suggest that GAKV and DARV are genetically distinct novel shrew-borne paramyxoviruses within the genus Henipavirus in the family Paramyxoviridae.

Shrews and rodents constitute the two largest groups of mammalian species [9] . They have been characterized as reservoir hosts of emerging zoonotic viruses that cause human diseases [54, 55] . Recently, many studies have been conducted on animal viromes in wildlife [56] [57] [58] . Using viral metagenomics techniques, influenza virus and herpesvirus were discovered in house shrews, urban rats, and bats [58] . In addition, arenavirus, astrovirus, hantavirus, hepatitis E virus, and ZIKV were detected in shrews and rodents [33, 55, 58] . Recently, we discovered that striped field mice harbored PAPV-1 and -2 [30] . In addition, BaVV and PMPV-1 were found in the bank vole (Myodes glareolus), and MMLV-1 and MMLV-2 were carried by Rungwe brush-furred rats (Lophuromys machangui) [29, 59] . In this study, we identified novel paramyxoviruses in samples from C. lasiura and C. shantungensis individuals.

Cophylogenetic relationships of zoonotic viruses and natural reservoir hosts play a critical role in understanding the evolutionary process for coadaptation, spillover, and host sharing [60] [61] [62] . Multiple natural reservoirs harbor Orthohantavirus species, and in most cases, it is a specific relationship of "one virus-one host" [63] . The congruent relationships between whole groups of hantaviruses and (sub)families of hosts proposed a view of the long-term coevolution of these viruses with their hosts [64, 65] . However, recent studies, based on cophylogenetic reconciliation and estimation of evolutionary rates and divergence times, demonstrate that local host-specific adaptation and preferential host switching account for the phylogenetic similarities between viruses and their mammalian hosts [66, 67] . Our cophylogenetic analysis indicates that shrew-borne henipaviruses and their host species have coadapted together differently from the phylogenies of rodent-borne paramyxoviruses. Thus, continuous and large-scale investigations will provide insights into the evolutionary history and dynamics of paramyxoviruses, harbored by shrews and rodents, in the family Orthoparamyxoviridae.

Triggering of innate antiviral genes modulates viral infectious diseases in humans and animals [68, 69] . The expression of antiviral genes, Mx1, Rsad2/Viperin, Isg15, and OAS1, was highly induced during early NiV and HeV infections in ferrets [70] . These robust innate antiviral responses were insufficient to prevent viral dissemination and tissue damage in vivo [71] . The isolate of GAKV from C. lasiura was evaluated for the innate antiviral response in human lung epithelial cells. We found that GAKV enabled replication and elicited the rapid induction of type I/III IFNs, ISGs, and pro-inflammatory cytokines in A549 cells. These results demonstrated that the novel paramyxovirus may infect and induce pro-inflammatory responses via the human respiratory tract. The characterization and effect of innate antiviral responses to shrew-borne henipaviruses in humans await for further studies.

In conclusion, this study reports two generically distinct novel paramyxoviruses, GAKV in C. lasiura and DARV in C. shantungensis. Phylogenetic inference and genomic characterization demonstrated that these viruses belong to henipaviruses within the family Paramyxoviridae. Infection with GAKV showed the upregulation of multiple human innate antiviral genes in vitro. Therefore, these results facilitate understanding of the molecular epidemiology, genetic characteristics and diversity, and zoonotic potential of shrew-borne paramyxoviruses, ROK. These observations raise awareness and caution among virologists and physicians about the henipa-like paramyxoviruses in shrews.

The following are available online at https://www.mdpi.com/article/10 .3390/v13102020/s1, Figure S1 : Phylogenetic trees using the coding sequences (CDSs) of N, P, M, F, G, and L genes of Gamak virus (GAKV) and Daeryong virus (DARV); Figure S2 : Phylogenetic analyses of the paramyxoviruses using partial genome sequences of L gene; Table S1 : Sequences of primers for quantitative reverse transcription-polymerase chain reaction; Table S2 : Characteristics of Crocidura species infected with Gamak virus (GAKV) and Daeryong virus (DARV) and the nucleotide position of GAKV and DARV RNA obtained in this study; Table S3 : Molecular prevalence of Gamak virus (GAKV) by region, sex, weight, and season in C. lasiura, the Republic of Korea from 2017 to 2018; Table S4 : Molecular prevalence of Daeryong virus (DARV) by region, sex, weight, and season in C. shantungensis, the Republic of Korea from 2017 to 2018; Table S5 : Genome coverage of consensus sequences for Gamak virus using next-generation sequencing (SISPA-based MiSeq) method; Table S6 : Genome coverage of consensus sequences for Gamak and Daeryong viruses using nextgeneration sequencing (RNA-Seq based HiSeq) method; Table S7 : Sequences of intergenic regions (IGRs) and transcriptional stop and start signals of Gamak virus; Table S8 : Sequences of intergenic regions (IGRs) and transcriptional stop and start signals of Daeryong virus; Table S9 : Amino acid similarities of GAKV and DARV with other viruses in Paramyxoviridae. 

Data Availability Statement: All the data generated for this publication have been included in the current manuscript.

Risk factors for human disease emergence

Host range and emerging and reemerging pathogens

Global shifts in mammalian population trends reveal key predictors of virus spillover risk

Spillover and pandemic properties of zoonotic viruses with high host plasticity

Assessing the global threat from Zika virus

WHO Zika Causality Working Group. Zika Virus Infection as a Cause of Congenital Brain Abnormalities and Guillain-Barré Syndrome: Systematic Review

New methodologies for the estimation of population vulnerability to diseases: A case study of Lassa fever and Ebola in Nigeria and Sierra Leone

Genomic Epidemiology and Active Surveillance to Investigate Outbreaks of Hantaviruses

Zoonotic spillover infections with Borna disease virus 1 leading to fatal human encephalitis, 1999-2019: An epidemiological investigation

Isolation and characterization of a novel arenavirus harbored by Rodents and Shrews in Zhejiang province

Discovery of a Highly Divergent Coronavirus in the Asian House Shrew from China Illuminates the Origin of the Alphacoronaviruses

Seewis virus, a genetically distinct hantavirus in the Eurasian common shrew (Sorex araneus)

Newfound Hantavirus in Chinese Mole Shrew

Characterization of Imjin Virus, a Newly Isolated Hantavirus from the Ussuri White-Toothed Shrew (Crocidura lasiura)

Divergent ancestral lineages of newfound hantaviruses harbored by phylogenetically related crocidurine shrew species in Korea

Identification of novel and diverse rotaviruses in rodents and insectivores, and evidence of cross-species transmission into humans

Discovery of a highly divergent hepadnavirus in shrews from China

Applications of Next-Generation Sequencing Technologies to Diagnostic Virology

Viral metagenomics as an emerging and powerful tool in veterinary medicine

New viruses in veterinary medicine, detected by metagenomic approaches

Next-generation sequencing technology in clinical virology

Evaluation of viral genome assembly and diversity estimation in deep metagenomes

Viral Diversity of House Mice

Zoonotic Potential of Emerging Paramyxoviruses: Knowns and Unknowns

Bats host major mammalian paramyxoviruses

Discovery and genome characterization of three new Jeilongviruses, a lineage of paramyxoviruses characterized by their unique membrane proteins

Novel Paju Apodemus paramyxovirus 1 and 2, harbored by Apodemus agrarius in the Republic of Korea

Genome Sequence of Ruloma Virus, a Novel Paramyxovirus Clustering Basally to Members of the

Molecular characterization of feline paramyxovirus in Japanese cat populations

Molecular Detection and Genetic Characterization of Novel RNA Viruses in Wild and Synanthropic Rodents and Shrews in Kenya

Isolation and characterization of novel bat paramyxovirus B16-40 potentially belonging to the proposed genus

The Epidemiological Characteristics of the Korean Bat Paramyxovirus between

A review of Nipah and Hendra viruses with an historical aside

Nipah virus: Transmission of a zoonotic paramyxovirus

The emergence of Nipah virus, a highly pathogenic paramyxovirus

Emerging trends of Nipah virus: A review

Molecular epidemiology of paramyxoviruses in Zambian wild rodents and shrews

Novel Henipa-like Virus, Mojiang Paramyxovirus, in Rats

Idiosyncratic Mòjiāng virus attachment glycoprotein directs a host-cell entry pathway distinct from genetically related henipaviruses

Functional Analysis of the Fusion and Attachment Glycoproteins of Mojiang Henipavirus

Trimmomatic: A flexible trimmer for Illumina sequence data

Fast gapped-read alignment with Bowtie 2

NCBI reference sequences (RefSeq): A curated non-redundant sequence database of genomes, transcripts and proteins

SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing

Rapid evaluation and quality control of next generation sequencing data with FaQCs

Sensitive and Broadly Reactive Reverse Transcription-PCR Assays To Detect Novel Paramyxoviruses

Evolution of the cytochrome b gene of mammals

MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets

Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10

dendextend: An R package for visualizing, adjusting and comparing trees of hierarchical clustering

Rodent-borne diseases and their risks for public health

Detection and characterization of three zoonotic viruses in wild rodents and shrews from Shenzhen city

Viral gut metagenomics of sympatric wild and domestic canids, and monitoring of viruses: Insights from an endangered wolf population

Marine viruses discovered via metagenomics shed light on viral strategies throughout the oceans

Detection of Virus-Related Sequences Associated with Potential Etiologies of Hepatitis in Liver Tissue Samples From Rats, Mice, Shrews, and Bats

Genetic characterization of bank vole virus (BaVV), a new paramyxovirus isolated from kidneys of bank voles in Russia

A Cophylogenetic Perspective of RNA-Virus Evolution

Phylogenomic Characterization of Lopma Virus and Praja Virus, Two Novel Rodent-Borne Arteriviruses

Host and viral traits predict zoonotic spillover from mammals

Hantavirus Coevolution with Their Rodent Hosts

Reconstructing the evolutionary origins and phylogeography of hantaviruses

Evolution of hantaviruses: Co-speciation with reservoir hosts for more than 100MYR

Host switch during evolution of a genetically distinct hantavirus in the American shrew mole (Neurotrichus gibbsii)

Hantavirus Evolution in Relation to Its Rodent and Insectivore Hosts: No Evidence for Codivergence

Cytokines and respiratory syncytial virus infection

Hendra and Nipah viruses: Different and dangerous

Host gene expression profiles in ferrets infected with genetically distinct henipavirus strains

Zoonotic potential of a novel bat morbillivirus. bioRxiv 2021

We thank Su-Am Kim and Hyebin Jeong for trapping the small mammals and preparing the figures, respectively.

The authors declare no competing financial interests.