key: cord-0883555-gmtqj6bt
authors: Zeng, Jian; Yang, Zijun; Guo, Wentao; Wang, Xiaochun; Yang, Shixing; Shen, Quan; Wang, Hao; Zhang, Wen
title: Identification and genome characterization of novel parechovirus sequences from Hipposideros armiger in China
date: 2022-05-15
journal: Virol J
DOI: 10.1186/s12985-022-01806-1
sha: dcb938fb57c5d0108acfb2691f1c191ac74f3edb
doc_id: 883555
cord_uid: gmtqj6bt

BACKGROUND: Bats were identified as a natural reservoir of emerging and re-emerging infectious pathogens threatening human health and life. METHODS: This study collected 21 fecal samples of Hipposideros armiger in Mengla County of Xishuangbanna Prefecture Yunnan Province to combine one pool for viral metagenomic sequencing. RESULTS: Two nearly complete genomes of parechoviruses, BPeV11 and BPeV20, were sequenced. Genome analysis revealed that BPeV11 and BPeV20 follow a 3-3-4 genome layout: 5′ UTR-VP0-VP3-VP1-2A-2B-2C-3A-3B-3C-3D-3′ UTR. The prevalence of BPev11 and BPev20 by Nested-PCR showed that 1 of 21 fecal samples was positive. Based on amino acid identity comparison and phylogenetic analysis of P1, 2C, and 3D, BPeV11 and BPeV20 were closely related to but distinct from FPeVs. CONCLUSION: It was probably proposed to be a novel species in the genus Parechovirus of the family Picornaviridae. The isolation of BPev11 and BPev20 from H. armiger in China is the first complete genome of parechovirus isolations from bat feces of the genus Hipposideros. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12985-022-01806-1.

Parechovirus C [9] (formerly Sebokele virus 1) contains only one genotype. Parechovirus D [10] (ferret parechovirus), Parechovirus E [11] (falcon parechovirus), and Parechovirus F [12] (gecko parechovirus) are the same as Parechovirus C. The knowledge of picornaviruses host range, geographical distribution and genome organization has recently exploded due to the use of highthroughput sequencing and the identification of novel picornaviruses from various species [10] .

Bats, the only flying mammal and account for more than 20% of the subsistent mammals, were recently identified as a natural reservoir of emerging and re-emerging infectious pathogens [13] , many of which could spillover into animal and human populations, such as severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), coronaviral disease-19, Nipah virus, Hendra virus, and Ebola virus [14, 15] . Hence, investigating viruses in bats is critical for improved control and prevention of large epidemics.

In this study, two nearly complete genomes of parechoviruses were sequenced and analyzed from the fecal samples of Hipposideros armiger in Mengla County of Xishuangbanna Prefecture Yunan Province by metagenomic analysis.

During 2017, a total of 21 fresh fecal samples were collected from wild H. armiger. The study was conducted in Mengla County of Xishuangbanna Prefecture Yunnan province, which is connected with Laos in the east and south and faces Myanmar across lancang River in the west, as shown in Fig. 1 . All samples were collected with disposable materials, shipped on dry ice and stored at − 80 °C for further study. The collected fecal samples were mixed into a group, suspended in 600 µl of Dulbecco's phosphate-buffered saline (DPBS), and then vigorously vortex oscillation for 5 min. The 500 µl supernatants were then collected from each pool after centrifugation (5 min, 15,000 g, 4 °C).

The 500 µl supernatant was purified through a 0.45-μm filter (Millipore) to remove eukaryotic and bacterial cell-sized particles. The filtrate was treated for 60 min at 37 °C with a DNases mixture (Turbo DNase from Ambion, Baseline-ZERO from Epicentre), benzonase (Novagen), and RNase (Fermentas) to digest unprotected nucleic acid. Nucleic acids (total DNA and RNA) were then extracted using a QIAamp Viral RNA Mini Kit (QIAGEN) following the manufacturer's instructions. Total nucleic acids were subjected to RT reactions with SuperScript III reverse transcriptase (Invitrogen), following second-strand cDNA synthesis with Large (Klenow) fragment (NEB). The library was then constructed using Nextera XT DNA Sample Preparation Kit (Illumina) and sequenced using the MiSeq Illumina platform with 250 base pair-ends with dual barcoding.

Paired-end reads of 250 bp generated by MiSeq sequencing were debarcoded using vendor software from Illumina. An in-house analysis pipeline running on a 32-node Linux cluster was used to process the data. Reads were considered duplicates if bases 5 to 55 were identical and only one random copy of duplicates was kept. Clonal reads were removed, and low-sequencingquality tails were trimmed using Phred. Adaptors were trimmed using VecScreen with the default parameters, which uses NCBI BLASTn with specific parameters designed for adapter removal. The cleaned reads were then compared to an in-house non-virus non-redundant (NVNR) protein database to remove false-positive viral hits using DIAMOND BLASTx search with default parameters [16] . The NVNR database was compiled using non-viral protein sequences extracted from an NCBI nr fasta file (based on annotation taxonomy, excluding the virus kingdom). Then, taxonomic classification for DIA-MOND results was parsed using MEGAN to perform the LCA-assignment algorithm according to default parameters. Gene assembly, prediction, and annotation were completed with Geneious software [17] .

Nested PCR was performed using rTaq DNA Polymerase (Takara) to amplify complete DNA or RNA and determine whether exist the target viruses. The specific primer sequences are shown in Table 1 .

The predicted potential proteins were aligned with their corresponding homologs of reference viruses using the MUSCLE multiple sequence alignment program with default settings [18] . The RdRp is the only conservedsequence domain across all RNA viruses and was used for phylogenetic inference [19] . All phylogenetic analysis was performed based on a Bayesian method implemented in MrBayes version 3.2.7 [20, 21] . In the MrBayes analyses, we used two simultaneous runs of Markov chain Monte Carlo sampling, and the runs were terminated upon convergence (standard deviation of the split frequencies < 0.01) [22] . The visualization and beautification of the phylogenetic trees were achieved by Figtree version 1.4.4 (available from http:// tree. bio. ed. ac. uk/ softw are/ figtr ee/).

All protein prediction was conducted by Geneious prime version 2019.2.3 [19] . The conserved domains were determined using the NCBI conserved domain search in combination with the Pfam conserved domain search [23, 24] . The cleavage sites of the BPev11 and BPev20 polyproteins were predicted by sequence alignment comparisons to the polyproteins of other viruses within the genus Parechovirus (See Additional file 1) using Geneious prime version 2019.2.3. The same approach was implemented in order to predict the polyprotein cleavage sites of three other closely related fish picornaviruses (Wenling bighead beaked sandfish picornavirus, Guangdong spotted longbarbel catfish picornavirus, West African lungfish picornavirus, and clownfish picornavirus) that had not previously been annotated [12, 25] . Pairwise genetic comparisons of the aa sequences of the P1, 2C, and 3D regions of BPev11 and BPev20 polyprotein were compared to those parechoviruses ( Fig. 2B ) using the Sequence Demarcation Tool v1.2 [26] , with the MUSCLE alignment option implemented.

The 21 fecal samples of H. armiger in the Yunnan Province were combined into one pool for viral metagenomic sequencing. The Illumina MiSeq outputted a total of 581,026 of 250-base pair-end reads and then removed duplicate, low-quality, and cellular reads. After de novo assembly, we acquired two contigs that be similar to parechovirus. Then, the two novel parechovirus sequences, BatPeVyn11 (BPeV11) and BatPeVyn20 (BPeV20) were identified by reads mapping. Nested-PCR was employed to investigate the prevalence of BPev11 and BPev20 using a primer pair.

Amplicons of approximately 316 and 330 base pairs were obtained from 2 of 21 fecal samples (sample 11 and 20), respectively. The primer sequences are shown in Table 1 .

The nearly complete genomic sequence of BPeV11 is 7143 nucleotides (nt) in length, which is within the range found for other picornaviruses (6938-9035 nt). A single ORF (nt 522-7142) was predicted to encode a polyprotein precursor of 2207 aa, which can be artificially divided into three parts: P1, P2, and P3. Similarly, the nearly complete 7096-nt-long RNA genome of BPeV20 has a single ORF (nt 320 -7048), which encodes a polyprotein precursor of 2243 aa ( Fig. 2A ). BPev11 and BPev20 both follow a 3-3-4 genome layout: 5′ UTR-VP0-VP3-VP1-2A-2B-2C-3A-3B-3C-3D-3′ UTR ( Fig. 2A) . The Predicted cleavage sites for genes of BPev11 and BPev20 are shown in Fig. 2A and Table 1 , based on a Muscle-alignment with fourteen previously annotated parechoviruses showing cleavage sites. In the polyprotein of BPeV11 and BPeV20, five conserved domains could be identified by an NCBI combined with a Pfam conserved domain search ( Fig. 2A) . The specific sites and lengths are listed in Table 2 .

According to BLASTx search, the ORF sequence of BPev11 and BPev20 shared the 50.44% and 50.47% identity at aa level with that of Ferret parechovirus (FPeV) isolate MpPeV1 (GenBank no. NC_034453) collected from Mustela putorius furo, respectively [3] . The 2C region of the BPev11 exhibited the greatest (60%) aa identity to FPeV, while its P1 and 3D displayed 46 and 56% identity to FPeV, respectively (Fig. 2B ). For BPev20, the 3D region showed the highest (59%) aa identity to FPeV, while its P1 and 2C displayed 46 and 56% identity to FPeV, respectively (Fig. 2B ). BPev11 and BPev20 were most closely to the genus Parechovirus of picornaviruses. Therefore, the representative members in the genus Parechovirus and other representative genera and species in Picornaviridae were selected as reference strains for the phylogenetic analysis. The phylogenetic relationships between BPev11 and BPev20 and the representative picornaviruses based on aa sequences of the different picornavirus coding regions (P1and 3CD) are shown in Fig. 3 . Phylogenetic analysis of the P1 confirmed that BPev11 and BPev20 formed a monophyletic branch with the members within Parechoviruses E (two Ferret parechoviruses and an unclassified parechovirus), while GPeV in Parechovirus F were also involved in this monophyletic clade based on the 2C and 3CD phylogenetic trees. However, BPev11 and BPev20 have distant relatedness to the viruses in same branch. According to the species demarcation criteria proposed 

With the rise of viral metagenomics analysis, new batborne viruses have been continuously discovered around the world, including influenza A virus [27] , Phlebovirus [28] , and Banyangvirus [29] . From the wide geographical locations of the various bat viruses detected so far, it is almost certain that we will see more and more disease outbreaks caused by bat viruses [14] . Bats harbour a more significant proportion of zoonotic viruses per host species than other mammalian orders [30, 31] . In some cases, outbreaks have been linked to bat roosting or foraging in close proximity to human settlements [32] [33] [34] . The driving factors for the increasing spillover events of bat viruses, particularly bat CoVs, are complex and most likely a combination of bat habitat disruption through climate change, increased urbanization pressure from humans, wildlife trade and animal markets [35, 36] . The SARS CoV outbreak in China, which caused more than 8000 cases of severe respiratory disease in humans resulting in 10% Mortality [37] , was linked to Rhinolophus sp. bats and the wildlife trade [38] . Up to today, the virus members from families of Rhabdoviridae, Orthomyxoviridae, Paramyxoviridae, Coronaviridae, Togaviridae, Flaviviridae, Bunyaviridae, Reoviridae, Arenaviridae, Herpesviridae, Picornaviridae, Hepesviridae and Adenoviridae, have been isolated from different bat species [39] , but no complete genome of Parechovirus has been found, only a 1080 nt parechovirus contig was recovered in a previous study [40] . As a natural reservoir, the bats deserve more studies to prevent the outbreaks of diseases caused by viruses. This study collected 21 fecal samples of H. armiger from Mengla County of Xishuangbanna Prefecture, Yunnan province, to process metagenomic analysis, then acquired two nearly complete sequences BPev11 and BPev20. The isolation of BPev11 and BPev20 from H. armiger in China is the first complete genome of parechovirus isolations from bat feces of the genus Hipposideros. Nested-PCR was employed to investigate the prevalence of BPev11 and BPev20 using a primer pair. The positive rates of BPev11 and BPev20 both were 1/21. BPev11 showed an electrophoresis band in sample 11, while BPev20 showed an electrophoresis band in sample 20 (See Additional file 2). Parechovirus is a genus of RNA viruses with a poly(A) tail. Regrettably, after a series of trials, such as a 5′ Rapid Amplification for cDNA End (RACE) PCR and 3' RACE, the 5'UTR and 3'UTR of BPev11 and BPev20 are not complete. Sequence identity matrices for P1, 2C, and 3D regions (Fig. 2B ) and phylogenetic tree for P1 and 3D regions of BPev11 and BPev20 suggest BPev11 and BPev20 may be a novel species of genus Parechovirus. Although BPev11 and BPev20 are more distantly related to the viruses of Parechovirus A and Parechovirus B, further studying the host range restriction and pathogenicity is necessary.

This study found two nearly complete genomes of parechoviruses, BPeV11 and BPeV20, in 21 fecal samples of H. armiger collected in Mengla County of Xishuangbanna Prefecture Yunnan Province, China. The prevalence of BPev11 and BPev20 by Nested-PCR showed that 1 of 21 fecal samples was positive. Based on amino acid identity comparison and phylogenetic analysis of P1, 2C, and 3D, BPeV11 and BPeV20 were closely related to but distinct from FPeVs. It was probably proposed to be a novel species in the genus Parechovirus of the family Picornaviridae. Until now, this is the first time that complete genome of parechovirus has been found in bat feces. However, further research is needed to uncover BPeV11 and BPeV20 in the pathogenic mechanism in animals and humans. 

ICTV virus taxonomy profile: picornaviridae

Complete genome sequences of the first parechoviruses a associated with sporadic pediatric acute gastroenteritis in Russia

Enteroviruses and parechoviruses: echoviruses, coxsackieviruses, and others

Pediatric parechovirus infections

Enterovirus and parechovirus infection in children: a brief overview

A new picornavirus isolated from bank voles (Clethrionomys glareolus)

Zoonotic Ljungan virus associated with central nervous system malformations in terminated pregnancy

Parechovirus a pathogenesis and the enigma of genotype A-3

Human parechovirus: an increasingly recognized cause of sepsis-like illness in young infants

Ljungan/Sebokele-like picornavirus in birds of prey, common kestrel (Falco tinnunculus) and red-footed falcon (F. vespertinus)

The evolutionary history of vertebrate RNA viruses

Bat and virus

Viruses in bats and potential spillover to animals and humans

A pneumonia outbreak associated with a new coronavirus of probable bat origin

Fast and sensitive protein alignment using DIAMOND

Potential and active functions in the gut microbiota of a healthy human cohort

MUSCLE: multiple sequence alignment with high accuracy and high throughput

Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data

MRBAYES: Bayesian inference of phylogenetic trees

MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space

Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses

The Pfam protein families database in 2019

CDD/SPARCLE: the conserved domain database in 2020

Characterization of a novel picornavirus isolated from moribund aquacultured clownfish

SDT: a virus classification tool based on pairwise sequence alignment and identity calculation

A distinct lineage of influenza A virus from bats

Malsoor virus, a novel bat phlebovirus, is closely related to severe fever with thrombocytopenia syndrome virus and heartland virus

Zwiesel bat banyangvirus, a potentially zoonotic • fast, convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year • At BMC, research is always in progress. Learn more biomedcentral.com/submissions Ready to submit your research Ready to submit your research ? Choose BMC

Huaiyangshan banyangvirus (Formerly known as SFTS)-like banyangvirus in Northern bats from

Host and viral traits predict zoonotic spillover from mammals

A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special?

A metagenomic viral discovery approach identifies potential zoonotic and novel mammalian viruses in Neoromicia bats within South Africa

Middle east respiratory syndrome coronavirus in bats, Saudi Arabia

Identification of a novel coronavirus in bats

Global trends in emerging infectious diseases

Impacts of biodiversity on the emergence and transmission of infectious diseases

Severe acute respiratory syndrome (SARS)-paradigm of an emerging viral infection

Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor

Bats: important reservoir hosts of emerging viruses

A viral metagenomic survey identifies known and novel mammalian viruses in bats from Saudi Arabia

Publisher's Note

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

We are grateful for the valuable support from National Key Research and Development Programs of China for Virome in Important Wildlife, Jiangsu Provincial Key Research and Development Projects. 

The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s12985-022-01806-1.Additional file 1. BPev11 and BPev20 MAFFT-alignment with other parechovirus polyproteins.Additional file 2. The result of Nested PCR about BPev11 and BPev20.