Skip to main content
Download PDF

Molecular detection and genomic characterization of Samak Micromys paramyxovirus-1 and -2 in Micromys minutus, Republic of Korea

Virology Journal volume 21, Article number: 255 (2024) Cite this article

Abstract

Background

The discovery of viruses in small mammalian populations, particularly rodents, has expanded the family Paramyxoviridae. The overlap in habitats between rodents and humans increases the risk of zoonotic events, underscoring the importance of active surveillance. Rodent species, such as Apodemus agrarius, are natural hosts for Paramyxoviridae in the Republic of Korea (ROK). However, it is unknown whether Paramyxoviridae is present in Micromys minutus, another common rodent.

Method

Here, we screened M. minutus collected from the Gangwon Province in the ROK for paramyxoviruses using nested polymerase chain reaction and confirm positive samples by next-generation metagenomic sequencing. Complete paramyxovirus genomes were further characterized by phylogenetic analysis, amino acid similarity, secondary structure, and cophylogeny.

Result

Overall, 57 of 145 (39.3%) M. minutus kidney samples tested positive for paramyxoviruses. Among them, four whole genome sequences were identified and clustered within the genus Jeilongvirus. One sequence was determined as Samak Micromys paramyxovirus 1 (SMPV-1; 19,911 nucleotides long) and three sequences as Samak Micromys paramyxovirus 2 (SMPV-2; 18,199 nucleotides long). SMPV-1 has a smaller hydrophobic gene and a longer glycoprotein gene than SMPV-2. Cophylogenetic analysis suggests that SMPV-1 evolved through co-divergence, whereas SMPV-2 was inferred to have undergone transfer events.

Conclusion

These findings highlight the prevalence of paramyxoviruses in the wild and the potential of M. minutus as a natural viral reservoir. The discovery of SMPV-1 and SMPV − 2 also reveals the genetic diversity and evolutionary history of the genus Jeilongvirus in the Paramyxoviridae.

Introduction

Paramyxoviruses pose a critical public health threat as emerging infectious diseases (EIDs) [1]. Paramyxoviruses can infect diverse hosts, including large mammals (primate, cattle, canine, and feline), small mammals (bats, shrews, and rodents) and marine species [2]. Hendra and Nipah viruses, which are members of the family Paramyxoviridae, have caused outbreaks through spillover events in horses and pigs [3]. Recent studies have reported novel paramyxoviruses of the Henipavirus and Jeilongvirus genera isolated from small mammals in the Republic of Korea (ROK) [4, 5]. A shrew-borne Henipavirus, Langya virus, was reported as a novel etiological agent of respiratory infectious diseases from clinical samples [6]. Investigation of novel paramyxoviruses as potential EIDs is necessary because of the prevalence of zoonotic events within this virus family [7].

Paramyxoviridae are enveloped, negative-sense, non-segmented RNA genomes classified into nine subfamilies, encompassing 23 genera and 153 species [8]. The virus genomes range from 14,296 to 20,148 nucleotides (nt) and encode six to ten proteins including a nucleoprotein (N), polymerase-related/nonstructural V/C proteins (P/V/C), a matrix protein (M), a fusion protein (F), a receptor binding protein (RBP), and a polymerase protein (large protein L). The receptor glycoprotein genes differ between genera according to the protein function [9].

The Respirovirus and Ferlavirus genera have a hemagglutinin-neuraminidase protein (HN), the Henipaviruses have a glycoprotein (G), and the Morbilliviruses have a hemagglutinin (H) [10]. Additional structural proteins, including the small hydrophobic protein (SH) and the transmembrane protein (TM), are limited to certain genera of paramyxoviruses. The SH gene is found in Rubulavirus, Pneumovirus, Metapneumovirus, and several members of Jeilongvirus whereas the TM gene is found only in Jeilongvirus [11].

Jeilongvirus is a recently established genus belonging to the Orthoparamyxovirinae subfamily [12, 13]. Jeilongvirus members have been found in rodents, hedgehogs, bats, and cats [9]. In addition to the TM gene, Jeilongvirus is distinct from other Paramyxoviruses because of the large size of its G protein and high variance in its C-terminal region [14]. The Jeilongvirus includes J virus from Mus musculus in Australia, Belerina virus from Erinaceus europaeus in Belgium, Ruloma virus from Lophuromys machangui in Tanzania, and Paju Apodemus paramyxovirus 1 and 2 (PAPV-1 and − 2) from Apodemus agrarius in the ROK [5, 11, 15, 16]. Recently, China reported two additional novel Jeilongvirus from Rattus tanezumi and A. agrarius [17]. Another observation on J virus, a Jeilongvirus detected in 1972, reported its ability to infect three different species of small wild mammals, specifically A. agrarius, A. peninsulae, and Microtus fortis [18].

In this study, we asked whether Paramyxoviruses were harbored by Micromys minutus in the ROK. M. minutus is known for its nest-building ability on the stems of wheat or other plants and widely distributed in temperate regions of Asia and Europe [19]. M. minutus has been reported as a reservoir of medically relevant pathogens such as Campylobacter jejuni [20]. However, the natural virus diversity of M. minutus remains unexplored. Wild small mammals are subjected to virus surveillance because of their significance as natural hosts for zoonotic viruses [21]. The increasing size of agriculture in the ROK raises the risk of viral spillover due to more frequent contact between farmers and small wild mammals [22, 23]. Screening and exploration of RNA viruses targeted within specific mammal species could maximize zoonotic virus discovery and aid in preparedness for future pandemics [24].

Materials and methods

Ethics statement

Animal trapping and all animal procedures were conducted according to protocols approved by the Institutional Animal Care and Use Committee of Hallym University (Hallym2016-37, Hallym 2018-6).

Animal trapping and organ sampling

M. minutus was captured from perilla fields in Chuncheon and Hongcheon, Gangwon Province, ROK (Supplementary Fig. 1). The small mammals were sacrificed under carbon dioxide. Individuals were identified using a field guide, and morphometric measurements were taken. A total of 145 rodents were necropsied and various tissues (serum, brain, lung, spleen, kidney, and liver) were collected aseptically. All samples were stored at -80 °C until further analysis.

Genetic identification of hosts by mitochondrial DNA (mtDNA) analysis

M. minutus species were identified based on field morphological characteristics and confirmed by sequencing of the mitochondrial cytochrome b (mtDNA cyt b) gene. Total DNA was extracted from kidney tissues using the TRIzol™ reagent. The full-length (1,140 nt) mtDNA cyt b was amplified as described in our previous study, and the amplicon was sequenced using the Sanger method [4]. The collected sequences were aligned using MUSCLE, and the results were analyzed using the IQ-TREE web server [25].

Molecular screening of paramyxovirus from M. minutus samples

Total RNA was extracted from M. minutus kidney tissues using the TRIzol™ reagent (Invitrogen, USA) according to the manufacturer’s instructions. Using random hexamers, cDNA was synthesized using a high-capacity RNA-to-cDNA kit (Applied Biosystems, Foster City, CA, USA). Paramyxovirus screening was performed using nested PCR reaction. The primers targeted the RNA-dependent RNA polymerase region of paramyxoviruses [26]. First and second PCRs were performed in a 25-µL reaction mixture containing 2.5 U of Supertherm DNA polymerase (JMR Holdings, Kent, UK), 2 µg of cDNA, 10 pM of each primer. Oligonucleotide primer sequences for the PCR were PAR-F1 (outer): 5ʹ-GAA GGI TAT TGT CAI AAR NTN TGG AC-3ʹ, PAR-F2 (inner): 5ʹ-GTT GCT TCA ATG GTT CAR GGN GAY AA-3, PAR-R (outer and inner): 5ʹ-GCT GAA GTT ACI GGI TCI CCD ATR TTN C-3ʹ. Initial denaturation was performed at 95℃ for 5 min, followed by six cycles of denaturation at 94℃ for 40 s, annealing at 37℃ for 40 s and elongation at 72℃ for 1 min, then thirty-two cycles of denaturation at 94℃ for 40 s, annealing at 42℃ for 40 s and elongation at 72℃ for 1 min (ProFlex PCR System, Life Technology, CA, USA).

PCR products were purified by MinElute® PCR purification kit (Qiagen, Hilden, Germany), and sequencing was performed in both directions of each PCR product using a BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Waltham, USA) on an automated sequencer (ABI 3730XL DNA Analyzer, Applied Biosystems). The forward and reverse sequences were assembled and manually checked using SeqMan of the DNAStar Lasergene (DNAStar Inc.). Sequence consensus was identified using the nucleotide-basic Local Alignment Search Tool (BLASTn) in the National Center for Biotechnology Information (NCBI) RefSeq database [27].

Metagenomic sequencing and de novo assembly for novel paramyxovirus identification

Libraries were prepared from total RNA using previously described protocols [5]. High-throughput sequencing was performed as paired-end 100 sequencing using HiSeq X10 (Illumina, USA). In this metagenomic approach, the reads were trimmed using Trimmomatic (v.0.36) to remove adapter sequences [28]. Additional host-sequence removal was done using Bowtie (v.2.2.6) [29]. The complete genome sequence of the species from the NCBI RefSeq database was used as the host reference [30]. The remaining reads were filtered for quality using FaQCs (v0.11.5), and de novo assembly was conducted using SPAdes (v3.15.5) [31, 32]. The assembled contigs were subsequently analyzed in the NCBI RefSeq database consisting of complete viral genomes by BLASTn (v2.6.0) [27, 30]. The contigs matched with paramyxovirus species were extracted and evaluated for coding sequences. For the reference mapping strategy, raw reads were trimmed, filtered, and mapped using CLC Genomics Workbench v.24.0.1 (CLC Bio, Cambridge, MA, USA). The Paramyxoviridae species reference was used as the mapping framework. Consensus genomic sequences of the novel virus were determined by combining viral contigs extracted from de novo assembly and reference mapping analysis, with an average depth cutoff of 10-fold.

Genomic characterization and analysis for N-glycosylation potential at glycoprotein gene

Viral genomic sequences were annotated using Geneious Prime (v. 2023.2.1). The open reading frames (ORF) extracted from the whole genome were translated and submitted to BLASTP to check for protein sequence homology. The complete genome of the virus species identified through BLASTP was subsequently collected to serve as the reference backbone for phylogenetic analysis.

For protein similarity analysis, the translated ORF of the novel virus was aligned with the coding sequence (CDS) translation of the curated viral sequence using CLUSTAL Omega 1.2.2, and the similarity percentage was evaluated based on the BLOSUM62 matrix [33, 34]. The translated glycoprotein gene of the novel paramyxovirus sequence was then analyzed using NetNGlyc-1.0, to predict the N-glycosylation sites [35]. The secondary structure prediction was performed using web-based Network Protein Sequence Analysis (NPSA), with Multivariate Linear Regression Combination (MLRC), Discrimination of Secondary Structure Class (DSC), and Profile Network from HeiDelberg (PHD) chosen as the default tools [36].

Phylogenetic and cophylogenetic analysis

Whole genome sequences were aligned and trimmed using MUSCLE 5.1 [37]. Phylogenetic trees were generated using maximum likelihood methods according to the best-fit substitution model in the IQ-TREE web server, and the tree was visualized using FigTree (v.1.4.4.; https://github.com/rambaut/figtree) software. The topologies were supported by 1,000 bootstrapping iterations.

Cophylogenetic analysis was performed using the systematic cophylogeny reconciliation tool called eMPRess [38]. This software reconciles symbiont and host trees using the duplication transfer-loss model with the assumption that a smaller cost increases the likelihood of the event [39]. The input was the Maximum Likelihood phylogenetic trees of cytochrome b from the host mitochondria and the L gene of the virus. The analysis was conducted using the following eMPRess parameters: duplication cost = 1, transfer cost = 1, and loss cost = 1. The original cost of reconciliation was lower than expected by chance (p < 0.01).

Results

Epidemiological survey of M. minutus and molecular prevalence of paramyxovirus

A total of 145 rodents were collected from two cities in Gangwon province (Hongcheon and Chuncheon) from March to April between 2016 and 2020. The rodents were categorized by year, trapping site, weight, and sex. Body weight was used to estimate rodent maturity. The collected animals comprised 81 males and 64 females. The animals were categorized as < 3 g for babies, 3–6 g for adolescents, 6–9 g for adults, and 9 g or more for the late-adults (Table 1). M. minutus was identified through mitochondrial Cytb gene sequencing (Supplementary Fig. 2). Phylogenetic analysis revealed that the animals in the present study shared a common ancestor with M. minutus from European countries.

Table 1 The epidemiological characteristics of Paramyxovirus infection in Micromys minutus
Full size table

A total of 57 (39.3%) M. minutus isolates tested positive for paramyxovirus using pan-paramyxovirus primers. The positivity rate was highest in late-adult animals, followed by adolescent and adult animals. In terms of sex, the positive result was higher in males than in females M. minutus in Gangwon Province, ROK (Table 1).

Metagenomic sequencing and genomic characterization of Samak Micromys paramyxovirus from M. Minutus

The de novo assembly of raw reads from samples Mm16-29, Mm16-32, Mm17-42, and Mm19-42 constructed long contigs of more than 18,000 base pairs. These contigs matched Jeilongvirus species with a percent identity below 85% and contained ORFs that matched the CDS of Jeilongvirus. After evaluating the coverage to generate consensus and re-mapping the raw reads to validate the sequence, we identified four complete novel paramyxovirus genome sequences from the M. minutus samples. The whole genome obtained from sample Mm17-42 had distinct characteristics compared to those obtained from samples Mm16-29, Mm16-32, and Mm19-42.

The length of the genome from Mm17-42, designated as Samak Micromys paramyxovirus (SMPV) 1, was 19,911 nt with a GC content of 40.8%. The genome structure of SMPV-1 consisted of eight genes in the order 3ʹ-N-P-M-F-SH-TM-G-L-5ʹ (Fig. 1). The whole genomes identified from Mm16-29, Mm16-32, and Mm19-42 shared similar genomic characteristics and were referred to as SMPV-2. The full-length genome of SMPV-2 was 18,199 nt with a GC content of 40.18% in Mm16-29, 40.15% in Mm16-32, and 40.19% in Mm19-42. The genome structure of SMPV-2 consisted of seven genes in the order 3ʹ-N-P-M-F-TM-G-L-5ʹ (Fig. 1 ). Each gene and their translation shared similarities with genes from other Paramyxoviridae(Fig. 2, Supplementary Tables 1 to 3). SMPV-1 and SMPV-2 had attachment glycoproteins of different sizes, although both infect the same host species. We attempted to observe their host immune-evasion potential through the proportion of amino acids and the glycosylation potential of the glycoprotein domain. Threonine, proline, and serine were the most abundant amino acids within the G domain (Supplementary Table 4). The glycoprotein of SMPV-1 contained 11.6% threonine, 9.6% serine, and 6.7% proline, whereas SMPV-2 consisted of 9.3% threonine, 8.8% proline, and 8.0% serine. N-glycosylation analysis of glycoprotein translation revealed seven potential sites on SMPV-1 and five on SMPV-2 (Supplementary Fig. 3). SMPV-1 exhibited recurring secondary structure patterns and potential glycosylation sites with PAPV-1 and Ninomys virus, whereas SMPV-2 did with PAPV-2 and Mount Mabu Lophuromys virus 1 (MMLPV-1). The start, stop, and intergenic region sequences of SMPV-1 and − 2 are listed in Table 2.

Fig. 1
figure 1

Comparative genome organization of Samak Micromys paramyxovirus 1 and 2 with other Paramyxoviruses. The genome organization of Samak Micromys parayxovirus 1 and 2 (SMPV-1 and − 2) compared with Jeilongvirus members and different genera within Paramyxoviridae family. SMPV-1 had longer whole genome sequence and G protein domain compared to SMPV-2. Abbreviations: N, nucleocapsid protein; P, phosphoprotein; M, matrix protein; F, fusion protein; SH, small hydrophobic protein; TM, transmembrane protein; G, glycoprotein; H, hemagglutinin protein; HN, hemagglutinin-neuraminidase protein; L, large protein

Full size image
Fig. 2
figure 2

Amino acid similarity of Samak Micromys paramyxovirus 1 and 2 with other species of Paramyxoviridae. The amino acid similarity percentage of each gene was based on BLOSUM62 matrix after multiple sequence alignment. Each gene was presented with a different color according to the legend. (A) Amino acid similarity percentage with SMPV-1 as the query. Ninomys virus showed high similarity with SMPV-1 (B) Amino acid similarity percentage with SMPV-2 as the query. PAPV-2 and MMLV-1 displayed high similarity with SMPV-2. Abbreviations: N, nucleocapsid protein; P, phosphoprotein; M, matrix protein; F, fusion protein; SH, small hydrophobic protein; TM, transmembrane protein; G, glycoprotein; H, hemagglutinin protein; HN, hemagglutinin-neuraminidase protein; L, large protein; SMPV-1: Samak Micromys paramyxovirus 1; SMPV-2: Samak Micromys paramyxovirus 1; PAPV-2: Paju Apodemus paramyxovirus 2; MMLV-1: Mount Mabu Lophuromys virus 1

Full size image
Table 2 Sequence of intergenic regions (IGR) and transcriptional start and stop signals of Samak Micromys paramyxovirus 1 and 2
Full size table

Phylogenetic and cophylogenetic analysis of SMPV

The complete and partial L genome sequences of these novel paramyxoviruses clustered with viruses of the Jeilongvirus genus and formed two distinct genetic lineages, SMPV-1 and SMPV-2 (Fig. 3 and Supplementary Fig. 4). The SMPV-1 cluster shared a common ancestor with Ninomys virus and Longquan Niviventer fulvescens jeilongvirus 2. The SMPV-2 cluster was closer to MMLPV-1 and PAPV-2. The SMPV-1 and − 2 clusters consisted of thirty-seven and twenty partial L-genome sequences, respectively. Cophylogenetic analysis resulted in eight co-divergence and nine transfer events (Fig. 4). According to the Maximum Parsimony Reconciliation (MPR) analysis, SMPV-1 and Ninomys viruses acquired their hosts through co-divergence, while SMPV-2 acquired its host through a transfer event. The co-divergence event was detected between Ninomys virus and SMPV-1, harbored by M. minutus from Belgium and Korea, with 100% probability of occurrence. When we traced back to the nearest common ancestor, a transfer event from PAPV-1, harbored by A. agrarius, was followed by a co-divergence event between viruses harbored by Rattus sp. (Beilong virus and Tailam virus) and Micromys sp. (Ninomys virus and SMPV-1). Another transfer event (with a probability of occurrence 50%) was detected from MMLPV-1, harbored by Lophuromys machangui, which became SMPV-2 in M. minutus. SMPV-2 is also indicated to have undergone a transfer event, resulting in PAPV-2, harbored by A. agrarius (50% probability of occurrence).

Fig. 3
figure 3

Phylogenetic tree of Samak Micromys paramyxovirus 1 and 2 within Paramyxoviridae members. The phylogenetic tree was constructed using maximum likelihood analysis by IQTREE web server, with GTR + F + I + G4 model chosen according to BIC and 1000 bootstrapping. The blue-colored label corresponds to Samak Micromys paramyxovirus 1 (SMPV-1), while the red-colored label represents Samak Micromys paramyxovirus 2 (SMPV-2). These two viruses cluster within the genus Jeilongvirus. SMPV-1 shared common ancestor with Ninomys virus, while SMPV-2 shared common ancestor with Paju Apodemus Paramyxovirus 2 and Mount Mabu Lophuromys virus 1

Full size image
Fig. 4
figure 4

Cophylogenetic analysis of Samak Micromys paramyxovirus 1 and 2 with members of Jeilongvirus and their respective host. The SMPV-1 and-2 are labelled with bold blue font. The black colored phylogenetic tree is the host tree based on mtDNA cytb sequences. The blue colored phylogenetic tree is the Jeilongvirus tree based on L gene of respective virus. The orange circle and number are the nodes for co-divergence and the probability of occurrence. The blue arrow showed the direction of the transfer event. The pink square and number are the nodes for transfer event, followed by the probability of occurrence in percent (%)

Full size image

Discussion

The virus surveillance focusing on small mammals continuously discovers novel paramyxoviruses [7]. Our study identifies two novel viruses within the M. minutus population in the Korean Peninsula. We obtained four whole-genome sequences, consisting of one SMPV-1 and three SMPV-2 sequences, using NGS. The genomes of SMPV-1 and − 2 follow the Paramyxoviridae genome organization criteria 3ʹ-N-P-M-F-(SH)-TM-G-L-5ʹ, with variations in the P gene attributed to an additional ORF. SMPV-1 has the SH gene and a larger G gene than SMPV-2. The presence of the SH gene separated Jeilongvirus into two groups. J virus, Beilongvirus, PAPV-1, and Ninomys virus harbor the SH gene, while PAPV-2, MMLPV-1, and MMLPV-2 lack it. However, the SH gene is not exclusive to Jeilongvirus, as Mumps virus (MuV) from the genus Orthorubulavirus carries this gene. Notably, SMPV-1 and − 2 exhibit different sizes and secondary structures of the G protein despite infecting the same species. The G protein is involved in the binding of viral particles to target cells [11]. The G proteins within Jeilongvirus species were varied, and the C-terminal region was enriched with proline (P), threonine (T), and serine (S) residues [14]. SMPV-1 and − 2 had dominant P/T/S residues in their glycoprotein amino acid composition and potential N-glycosylation sites. In addition, other rodent species similarly exhibit host sharing of phylogenetically distant paramyxoviruses, including Lophuromys machangui with MMLPV-1 and − 2, A. agrarius with PAPV-1 and − 2, and Myodes glareolus with Pohorje Myodes Paramyxovirus 1 (PMPV-1) and bank vole virus [5, 14, 42]. Our initial analysis of Jeilongvirus glycoprotein domain implied that the glycoprotein might elicit different host immune responses. Since co-infection with SMPV-1 and SMPV-2 was not observed in the specimens, each virus might selectively infect M. minutus individuals with distinct immune profiles. Further study characterizing the SMPV-1/-2 infected M. minutus immune response will improve our understanding of the host-virus interaction.

Mutations on RNA virus genomes are often dominated by synonymous substitutions in CDS maintaining the amino acid sequence of the protein [43]. The analysis on the CDS of SMPV-1 revealed that its genetic diversity did not translate to significant variation at the amino acid level. The amino acid sequence of SMPV-1 showed high similarity percentage with Ninomys virus, a Jeilongvirus found in M. minutus from Belgium, despite substantial phylogenetic divergence in complete and partial genome sequences [44]. RNA viruses utilize host cellular machinery for its replication, driving a coevolutionary process between the virus and the host [45]. Cophylogenetic reconstruction indicated that co-divergence events influenced the evolution of SMPV-1 and Ninomys viruses within the host population. M. minutus presumed to have low diversity owing to bottlenecking and lineage sorting during the Quaternary glacial cycles [46]. It is possible that speciation events of SMPV-1 and Ninomys virus coincided with the geographic genetic differentiation of M. minutus across Europe and Asia, with their evolution potentially constrained by the low genetic variability of the host. Observations in the Flaviviridae have shown that genetic constraints imposed by hosts and vectors pressure the viral coding region to evolve in adaptation to specific host-vector lineages [47]. Purifying selection during inter-host transmission might have eliminated mutations from intra-host pressure, leading to the fixation of advantageous mutations that efficiently spread the virus in larger populations [48]. This phenomenon might account for the observed disparity between genetic divergence and amino acid similarity in SMPV-1 and Ninomys virus, where amino acid compositions were likely conserved to preserve essential protein functions.

While SMPV-1 showed a phylogenetic relationship to the Ninomys virus, SMPV-2 shared a common ancestor with PAPV-2 and MMLV-1 hosted by other rodent species. Cophylogenetic analysis revealed that the cross-species transfer event from Lophuromys sp. preceded the speciation of SMPV-2. Since Lophuromys sp. is endemic to Africa, the transfer event likely occurred following their interaction with M. minutus populations that migrated from Europe to Asia during the Pliocene/Pleistocene.(19, 49) However, the current phylogenetic tree might not be fully comprehensive, as the transfer node only showed a 50% occurrence probability of the event; thus, there could be additional virus and host species serving as intermediaries bridging the transfer event from Lophuromys to Micromys. A transfer event detected from M. minutus to A. agrarius in the ROK antecedes PAPV-2. While direct interactions between the two host species on the Korean Peninsula might facilitate host-switching events, the potential for involvement of intermediate hosts or viruses in the evolution requires further investigation.

Compared with other rodent species, M. minutus rarely subjected to zoonosis surveillance. Being the smallest rodent in the family Muridae, M. minutus is difficult to collect [19]. Studies on other rodent and shrew genera in the ROK revealed their capacity as natural hosts of novel paramyxoviruses. Herein, SMPV detection within M. minutus population demonstrated its capacity to harbor and transmit rodent-borne paramyxoviruses. In particular, the SMPV positivity rate was higher in male compared to that in female M. minutus. This gender-specific prevalence of SMPV-1 and − 2 is similar to that of Orthohantavirus puumalaense in the bank vole of the Ural Mountains, Orthohantavirus hantanense, and Thottimvirus imjinense in small mammals of the Korean Peninsula [50,51,52]. Aggression is believed to influence the positivity rate, as a study on Orthohantavirus seoulense among R. norvegicus showed a correlation between wound severity and viral loads in male [53]. In addition to the behavioral aspect, sex dimorphism, hormones, and disparity in immune response between genders could be attributed to the higher SMPV prevalence in male specimens [54,55,56]. Since our study did not observe the wounding, behavior, or any physiological information of the animals, the transmission routes and risk factors for SMPV infection in M. minutus remain unclear.

The limitation of this study is the absence of SMPV-1 and-2 virus isolation and RACE PCR results owing to the small sample size. We were unable to provide additional information on SMPV pathogenesis or virulence factors. Moreover, the immune profile of M. minutus remains unexplored. Further investigations of M. minutus populations are required to address this issue in the ROK.

In conclusion, our study identified both SMPV-1 and − 2 from M. minutus population in the ROK through metagenomic analysis by next-generation sequencing. Although current evidence has not detected zoonotic transmission from the genus Jeilongvirus, the inclusion of SMPV-1 and − 2 sequences in paramyxoviruses is crucial for the surveillance of emerging infectious diseases. This study provides insights into the genetic diversity and evolutionary history of paramyxoviruses in small wild mammals.

Data availability

All the viral and Ctyb sequences can be accessed at Genbank with accession number listed in supplementary Table 5. The raw NGS data can be provided upon request from the corresponding author.

References

  1. Duprex WP, Dutch RE, Paramyxoviruses. Pathogenesis, vaccines, antivirals, and prototypes for pandemic preparedness. J Infect Dis. 2023;228(Supplement6):S390–7. https://doi.org/10.1093/infdis/jiad123.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Kuhn JH, Adkins S, Agwanda BR, Al Kubrusli R, Alkhovsky SV, Amarasinghe GK, et al. 2021 taxonomic update of phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders bunyavirales and Mononegavirales. Arch Virol. 2021;166(12):3513–66. https://doi.org/10.1007/s00705-021-05143-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Centers for Disease Control and Prevention (CDC). Outbreak of Hendra-Like virus–Malaysia and Singapore, 1998–1999. MMWR Morb Mortal Wkly Rep. 1999;48(13):265–9.

    Google Scholar 

  4. Lee SH, Kim K, Kim J, No JS, Park K, Budhathoki S, et al. Discovery and genetic characterization of Novel paramyxoviruses related to the Genus Henipavirus in Crocidura Species in the Republic of Korea. Viruses. 2021;13(10). https://doi.org/10.3390/v13102020.

  5. Lee SH, No JS, Kim K, Budhathoki S, Park K, Lee GY, et al. Novel Paju Apodemus paramyxovirus 1 and 2, harbored by Apodemus agrarius in the Republic of Korea. Virology. 2021;562:40–9. https://doi.org/10.1016/j.virol.2021.06.011.

    Article  CAS  PubMed  Google Scholar 

  6. Zhang XA, Li H, Jiang FC, Zhu F, Zhang YF, Chen JJ, et al. A zoonotic Henipavirus in Febrile patients in China. N Engl J Med. 2022;387(5):470–2. https://doi.org/10.1056/NEJMc2202705.

    Article  PubMed  Google Scholar 

  7. Thibault PA, Watkinson RE, Moreira-Soto A, Drexler JF, Lee B. Zoonotic potential of emerging paramyxoviruses: knowns and unknowns. Adv Virus Res. 2017;98:1–55. https://doi.org/10.1016/bs.aivir.2016.12.001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. International Committee on Taxonomy of Viruses: ICTV. Virus Taxonomy: 2023 Release. https://ictv.global/taxonomy (2023). Accesed 07 July 2024.

  9. Rima B, Balkema-Buschmann A, Dundon WG, Duprex P, Easton A, Fouchier R, et al. ICTV Virus Taxonomy Profile: Paramyxoviridae. J Gen Virol. 2019;100(12):1593–4. https://doi.org/10.1099/jgv.0.001328.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Aguilar HC, Henderson BA, Zamora JL, Johnston GP. Paramyxovirus glycoproteins and the membrane Fusion process. Curr Clin Microbiol Rep. 2016;3(3):142–54. https://doi.org/10.1007/s40588-016-0040-8.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Vanmechelen B, Vergote V, Merino M, Verbeken E, Maes P. Common occurrence of Belerina virus, a novel paramyxovirus found in Belgian hedgehogs. Sci Rep. 2020;10(1):19341. https://doi.org/10.1038/s41598-020-76419-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Audsley MD, Marsh GA, Lieu KG, Tachedjian M, Joubert DA, Wang L-F, et al. 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. J Gen Virol. 2016;97(3):581–92. https://doi.org/10.1099/jgv.0.000388.

    Article  CAS  PubMed  Google Scholar 

  13. Kuhn JH, Adkins S, Alkhovsky SV, Avšič-Županc T, Ayllón MA, Bahl J, et al. 2022 taxonomic update of phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders bunyavirales and Mononegavirales. Arch Virol. 2022;167(12):2857–906. https://doi.org/10.1007/s00705-022-05546-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Vanmechelen B, Bletsa M, Laenen L, Lopes AR, Vergote V, Beller L, et al. Discovery and genome characterization of three new jeilongviruses, a lineage of paramyxoviruses characterized by their unique membrane proteins. BMC Genomics. 2018;19(1):617. https://doi.org/10.1186/s12864-018-4995-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Jack PJ, Boyle DB, Eaton BT, Wang LF. The complete genome sequence of J virus reveals a unique genome structure in the family Paramyxoviridae. J Virol. 2005;79(16):10690–700. https://doi.org/10.1128/JVI.79.16.10690-10700.2005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Vanmechelen B, Meurs S, Zisi Z, Goüy de Bellocq J, Bletsa M, Lemey P, et al. Genome sequence of Ruloma Virus, a Novel Paramyxovirus Clustering Basally to members of the Genus Jeilongvirus. Microbiol Resour Announc. 2021;10(18). https://doi.org/10.1128/MRA.00325-21.

  17. Gan M, Hu B, Ding Q, Zhang N, Wei J, Nie T, et al. Discovery and characterization of novel jeilongviruses in wild rodents from Hubei, China. Virol J. 2024;21(1):146. https://doi.org/10.1186/s12985-024-02417-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Zhang Y, Zhang J, Wang Y, Tian F, Zhang X, Wang G, et al. Genetic diversity and expanded host range of J paramyxovirus detected in Wild Small mammals in China. Viruses. 2022;15(1). https://doi.org/10.3390/v15010049.

  19. Chen Z, Pei X, Song J, Song W, Shi Z, Onditi KO, et al. Systematics and evolutionary history of the genus Micromys (Mammalia: Rodentia: Muridae). Mammalian Biology. 2023;1–15. https://doi.org/10.1007/s42991-023-00360-9.

  20. Kim J, Guk JH, Mun SH, An JU, Kim W, Lee S, et al. The Wild mouse (Micromys minutus): Reservoir of a Novel Campylobacter jejuni strain. Front Microbiol. 2019;10:3066. https://doi.org/10.3389/fmicb.2019.03066.

    Article  PubMed  Google Scholar 

  21. Gibb R, Redding DW, Chin KQ, Donnelly CA, Blackburn TM, Newbold T, et al. Zoonotic host diversity increases in human-dominated ecosystems. Nature. 2020;584(7821):398–402. https://doi.org/10.1038/s41586-020-2562-8.

    Article  CAS  PubMed  Google Scholar 

  22. O’Shaughnessy SA, Kim M, Lee S, Kim Y, Kim H, Shekailo J. Towards smart farming solutions in the U.S. and South Korea: a comparison of the current status. Geogr Sustain. 2021;2(4):312–27. https://doi.org/10.1016/j.geosus.2021.12.002.

    Article  Google Scholar 

  23. Jones BA, Grace D, Kock R, Alonso S, Rushton J, Said MY, et al. Zoonosis emergence linked to agricultural intensification and environmental change. Proc Natl Acad Sci U S A. 2013;110(21):8399–404. https://doi.org/10.1073/pnas.1208059110.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Levinson J, Bogich TL, Olival KJ, Epstein JH, Johnson CK, Karesh W, et al. Targeting surveillance for zoonotic virus discovery. Emerg Infect Dis. 2013;19(5):743–7. https://doi.org/10.3201/eid1905.121042.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Trifinopoulos J, Nguyen L-T, von Haeseler A, Minh BQ. W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. 2016;44(W1):W232–5. https://doi.org/10.1093/nar/gkw256.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Tong S, Chern SW, Li Y, Pallansch MA, Anderson LJ. Sensitive and broadly reactive reverse transcription-PCR assays to detect novel paramyxoviruses. J Clin Microbiol. 2008;46(8):2652–8. https://doi.org/10.1128/JCM.00192-08.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10. https://doi.org/10.1016/S0022-2836(05)80360-2.

    Article  CAS  PubMed  Google Scholar 

  28. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9. https://doi.org/10.1038/nmeth.1923.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35(Database issue):D61–5. https://doi.org/10.1093/nar/gkl842.

    Article  CAS  PubMed  Google Scholar 

  31. Lo CC, Chain PS. Rapid evaluation and quality control of next generation sequencing data with FaQCs. BMC Bioinformatics. 2014;15(1):366. https://doi.org/10.1186/s12859-014-0366-2.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77. https://doi.org/10.1089/cmb.2012.0021.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Sievers F, Higgins DG. Clustal Omega for making accurate alignments of many protein sequences. Protein Sci. 2018;27(1):135–45. https://doi.org/10.1002/pro.3290.

    Article  CAS  PubMed  Google Scholar 

  34. Pearson P WR. Selecting the right similarity-scoring Matrix. Curr Protoc Bioinf. 2013;43(351–359). https://doi.org/10.1002/0471250953.bi0305s43.

  35. Gupta R, Brunak S. Prediction of glycosylation across the human proteome and the correlation to protein function. Pac Symp Biocomput. 2002:310–22.

  36. Combet C, Blanchet C, Geourjon C, Deléage G. NPS@: network protein sequence analysis. Trends Biochem Sci. 2000;25(3):147–50. https://doi.org/10.1016/s0968-0004(99)01540-6.

    Article  CAS  PubMed  Google Scholar 

  37. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7. https://doi.org/10.1093/nar/gkh340.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Santichaivekin S, Yang Q, Liu J, Mawhorter R, Jiang J, Wesley T, et al. eMPRess: a systematic cophylogeny reconciliation tool. Bioinformatics. 2020;37(16):2481–2. https://doi.org/10.1093/bioinformatics/btaa978.

    Article  CAS  Google Scholar 

  39. Ronquist F. Parsimony analysis of coevolving species associations. In: Page RDM, editor. Tangled trees: phylogeny, cospeciation and coevolution. Chicago: University of Chicago Press; 2002. pp. 22–64.

    Google Scholar 

  40. Abraham M, Arroyo-Diaz NM, Li Z, Zengel J, Sakamoto K, He B. Role of small hydrophobic protein of J paramyxovirus in virulence. J Virol. 2018;92(20). https://doi.org/10.1128/JVI.00653-18.

  41. Xu P, Li Z, Sun D, Lin Y, Wu J, Rota PA, et al. Rescue of wild-type mumps virus from a strain associated with recent outbreaks helps to define the role of the SH ORF in the pathogenesis of mumps virus. Virology. 2011;417(1):126–36. https://doi.org/10.1016/j.virol.2011.05.003.

    Article  CAS  PubMed  Google Scholar 

  42. Alkhovsky S, Butenko A, Eremyan A, Shchetinin A. Genetic characterization of bank Vole virus (BaVV), a new paramyxovirus isolated from kidneys of bank voles in Russia. Arch Virol. 2018;163(3):755–9. https://doi.org/10.1007/s00705-017-3639-z.

    Article  CAS  PubMed  Google Scholar 

  43. Simmonds P, Aiewsakun P, Katzourakis A. Prisoners of war - host adaptation and its constraints on virus evolution. Nat Rev Microbiol. 2019;17(5):321–8. https://doi.org/10.1038/s41579-018-0120-2.

    Article  CAS  PubMed  Google Scholar 

  44. Horemans M, Van Bets J, Joly Maes T, Maes P, Vanmechelen B. Discovery and genome characterization of six new orthoparamyxoviruses in small Belgian mammals. Virus Evol. 2023;9(2):vead065. https://doi.org/10.1093/ve/vead065.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Koonin EV, Gorbalenya AE. Evolution of RNA genomes: does the high mutation rate necessitate high rate of evolution of viral proteins? J Mol Evol. 1989;28(6):524–7. https://doi.org/10.1007/BF02602932.

    Article  CAS  PubMed  Google Scholar 

  46. Yasuda SP, Vogel P, Tsuchiya K, Han S-H, Lin L-K, Suzuki H. Phylogeographic patterning of mtDNA in the widely distributed harvest mouse (Micromys minutus) suggests dramatic cycles of range contraction and expansion during the mid-to late Pleistocene. Can J Zool. 2005;83(11):1411–20.

    Article  CAS  Google Scholar 

  47. Lobo FP, Mota BE, Pena SD, Azevedo V, Macedo AM, Tauch A, et al. Virus-host coevolution: common patterns of nucleotide motif usage in Flaviviridae and their hosts. PLoS ONE. 2009;4(7):e6282. https://doi.org/10.1371/journal.pone.0006282.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Duffy S, Shackelton LA, Holmes EC. Rates of evolutionary change in viruses: patterns and determinants. Nat Rev Genet. 2008;9(4):267–76. https://doi.org/10.1038/nrg2323.

    Article  CAS  PubMed  Google Scholar 

  49. Komarova VA, Kostin DS, Bryja J, Mikula O, Bryjová A, Čížková D, et al. Complex reticulate evolution of speckled brush-furred rats (Lophuromys) in the Ethiopian centre of endemism. Mol Ecol. 2021;30(10):2349–65. https://doi.org/10.1111/mec.15891.

    Article  PubMed  Google Scholar 

  50. Bernshtein AD, Apekina NS, Mikhailova TV, Myasnikov YA, Khlyap LA, Korotkov YS, et al. Dynamics of Puumala hantavirus infection in naturally infected bank voles (Clethrinomys Glareolus). Arch Virol. 1999;144(12):2415–28. https://doi.org/10.1007/s007050050654.

    Article  CAS  PubMed  Google Scholar 

  51. Klein TA, Kim HC, Chong ST, Kim JA, Lee SY, Kim WK, Nunn PV, Song JW. Hantaan virus surveillance targeting small mammals at nightmare range, a high elevation military training area, Gyeonggi Province, Republic of Korea. PLoS ONE. 2015;10(4):e0118483. https://doi.org/10.1371/journal.pone.0118483.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Lee SH, Kim WK, No JS, Kim JA, Kim JI, Gu SH, Kim HC, Klein TA, Park MS, Song JW. Dynamic circulation and Genetic Exchange of a shrew-borne Hantavirus, Imjin virus, in the Republic of Korea. Sci Rep. 2017;7:44369. https://doi.org/10.1038/srep44369.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Hinson ER, Shone SM, Zink MC, Glass GE, Klein SL. Wounding: the primary mode of Seoul virus transmission among male Norway rats. Am J Trop Med Hyg. 2004;70(3):310–7.

    Article  PubMed  Google Scholar 

  54. Brownstein DG, Gras L. Chromosome mapping of Rmp-4, a gonad-dependent gene encoding host resistance to mousepox. J Virol. 1995;69(11):6958–64. https://doi.org/10.1128/JVI.69.11.6958-6964.1995.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Geurs TL, Hill EB, Lippold DM, French AR. Sex differences in murine susceptibility to systemic viral infections. J Autoimmun. 2012;38(2):J245–53. https://doi.org/10.1016/j.jaut.2011.12.003.

    Article  CAS  PubMed  Google Scholar 

  56. Hannah MF, Bajic VB, Klein SL. Sex differences in the recognition of and innate antiviral responses to Seoul virus in Norway rats. Brain Behav Immun. 2008;22(4):503–16. https://doi.org/10.1016/j.bbi.2007.10.005.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Nayeon Jang (Hallym University) for supporting the study.

Funding

This study was provided by the Korea Institute of Marine Science & Technology Promotion (KIMST) funded by the Ministry of Oceans and Fisheries, Korea (20210466) and the Government-wide R&D to Advance Infectious Disease Prevention and Control, Republic of Korea (grant number: HG23C1623). This research was also supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) by the Ministry of Education (RS-202300249142) and by a Novo Nordisk Foundation PAD award to CBL (number: NF22SA0082041). In addition, this research was supported by the Basic Research Program through the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (2023R1A2C2006105) and the Institute for Basic Science (IBS) (IBS-R801‐D9‐A03), Republic of Korea.

Author information

Authors and Affiliations

  1. Department of Microbiology, College of Medicine, Hallym University, Chuncheon, 24252, Republic of Korea

    Augustine Natasha, Haryo Seno Pangestu, Jieun Park, Sara P. Prayitno, Bohyeon Kim, Shailesh Budhathoki & Won-Keun Kim

  2. Department of Molecular Microbiology and Center for Women Infectious Disease Research, Washington University School of Medicine in St. Louis, St. Louis, MO, USA

    Sarah E. Pye & Carolina B. López

  3. Department of Biomedical Science, College of Natural Sciences, Hallym University, Chuncheon, 24252, Republic of Korea

    Seung Hye Cho

  4. Department of Microbiology, Korea University College of Medicine, Seoul, 02841, Republic of Korea

    Kyungmin Park, Jongwoo Kim & Jin-Won Song

  5. BK21 Graduate Program, Biomedical Sciences, Korea University College of Medicine, Seoul, Republic of Korea

    Jin-Won Song

  6. College of Veterinary Medicine and Institute of Veterinary Science, Kangwon National University, Chuncheon, 24341, Republic of Korea

    Yeonsu Oh

  7. Department of Medical Genetics, College of Medicine, Hallym University, Chuncheon, 24252, Republic of Korea

    Jong Sun Lee & Jun Gyo Suh

  8. Institute of Medical Science, College of Medicine, Hallym University, Chuncheon, 24252, Republic of Korea

    Won-Keun Kim

Authors
  1. Augustine Natasha
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Sarah E. Pye
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Seung Hye Cho
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Haryo Seno Pangestu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jieun Park
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Kyungmin Park
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Sara P. Prayitno
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Bohyeon Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Jong Sun Lee
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Jongwoo Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Shailesh Budhathoki
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Yeonsu Oh
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Jin-Won Song
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Carolina B. López
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Jun Gyo Suh
    View author publications

    You can also search for this author in PubMed Google Scholar

  16. Won-Keun Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

A.N. and W-K.K. contributed to conceptualization, investigation, writing of the original draft, review, and editing. A.N. and S.E.P. contributed to data analysis and interpretation. S.H.C, H.S.P., J.P., K.P., S.P.P., B.K., J.K.,S.B., contributed equally to methodology, sample processing, data curation, and figures preparation. J.G.S. and J.S.L. contributed to specimen collection, figure preparation, and rodent data analysis. C.B.L, J.W.S, and Y.O. contributed equally to review and editing. All authors have reviewed and approved the final manuscript.

Corresponding authors

Correspondence to Jun Gyo Suh or Won-Keun Kim.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

12985_2024_2532_MOESM1_ESM.docx

Supplementary Material 1: Supplementary Figure 1. Distribution of Paramyxovirus-positive samples in the city of Chuncheon and Hongcheon, Gangwon Province. (A) The map shows the trapping sites in Chuncheon and Hongcheon, two cities located in Gangwon province. The red dot in Chuncheon represents the area where the captured M. minutus tested positive for paramyxovirus screening. Conversely, the blue dot in Hongcheon represents the area with negative screening results. (B) The picture shows a Micromys minutus specimen captured in this studySupplementary Figure 2. Mitochondrial cytb sequences result and phylogenetic tree of Micromys minutus. The phylogenetic tree was constructed using maximum likelihood analysis by IQTREE web server, with TPM3u+F+G4 model chosen according to BIC and 1000 bootstrapping. The blue colored label is mtDNA cytb sequences of the M. minutus in this study. The red colored label is the cytb sequences of M. minutus infected by Jeilongvirus from Horemans et al. study (2023). The sample from the Republic of Korea shared common ancestor with M. minutus from EuropeSupplementary Figure 3. The predicted N-glycosylation site of Samak Micromys paramyxovirus 1 and 2, in comparison with members of Paramyxoviridae. The potential N-glycosylation site marked by green bar crossing the purple threshold line. The secondary structure of the glycoprotein was added beneath the N-glycosylation graphics for additional visualization. Seven viruses are members of Jeilongvirus while Nipah virus from Henipavirus is added as control. Abbreviations: SMPV-1: Samak Micromys paramyxovirus 1; SMPV-2: Samak Micromys paramyxovirus 1; PAPV-1: Paju Apodemus paramyxovirus 1; PAPV-2: Paju Apodemus paramyxovirus 2; MMLV-1: Mount Mabu Lophuromys virus 1Supplementary Figure 4. Phylogenetic tree of SMPV-1 and -2 partial genome with other paramyxovirus. The partial sequences were generated from amplicon sequencing of the pan-paramyxovirus nested PCR product. The phylogenetic tree was constructed using maximum likelihood analysis by IQTREE web server, with TPM3u+F+G4 model chosen according to BIC and 1000 bootstrapping. The partial sequences generated two distinct clusters, in concordance with the phylogenetic tree generated from the whole genome sequences. The blue-colored label represents Samak Micromys paramyxovirus 1, and the red-colored label represents Samak Micromys paramyxovirus 2. The dot represents the samples from which the whole genome sequences originated

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Natasha, A., Pye, S.E., Cho, S.H. et al. Molecular detection and genomic characterization of Samak Micromys paramyxovirus-1 and -2 in Micromys minutus, Republic of Korea. Virol J 21, 255 (2024). https://doi.org/10.1186/s12985-024-02532-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-024-02532-6

Virology Journal

ISSN: 1743-422X