- Open Access
Detection of alpha- and betacoronaviruses in rodents from Yunnan, China
Virology Journal volume 14, Article number: 98 (2017)
Rodents represent the most diverse mammals on the planet and are important reservoirs of human pathogens. Coronaviruses infect various animals, but to date, relatively few coronaviruses have been identified in rodents worldwide. The evolution and ecology of coronaviruses in rodent have not been fully investigated.
In this study, we collected 177 intestinal samples from thress species of rodents in Jianchuan County, Yunnan Province, China. Alphacoronavirus and betacoronavirus were detected in 23 rodent samples from three species, namely Apodemus chevrieri (21/98), Eothenomys fidelis (1/62), and Apodemus ilex (1/17). We further characterized the full-length genome of an alphacoronavirus from the A. chevrieri rat and named it as AcCoV-JC34. The AcCoV-JC34 genome was 27,649 nucleotides long and showed a structure similar to the HKU2 bat coronavirus. Comparing the normal transcription regulatory sequence (TRS), 3 variant TRS sequences upstream the spike (S), ORF3, and ORF8 genes were found in the genome of AcCoV-JC34. In the conserved replicase domains, AcCoV-JC34 was most closely related to Rattus norvegicus coronavirus LNRV but diverged from other alphacoronaviruses, indicating that AcCoV-JC34 and LNRV may represent a novel alphacoronavirus species. However, the S and nucleocapsid proteins showed low similarity to those of LRNV, with 66.5 and 77.4% identities, respectively. Phylogenetic analysis revealed that the S genes of AcCoV-JC34, LRNV, and HKU2 formed a distinct lineage with all known coronaviruses.
Both alphacoronaviruses and betacoronaviruses were detected in Apodemus chevrieri in the Yunnan Province of China, indicating that Apodemus chevrieri is an important host for coronavirus. Several new features were identified in the genome of an Apodemus chevrieri coronavirus. The phylogenetic distance to other coronaviruses suggests a variable origin and evolutionary route of the S genes of AcCoV-JC34, LRNV, and HKU2. These results indicate that the diversity of rodent coronaviruses is much higher than previously expected. Further surveillance and functional studies of these coronaviruses will help to better understand the importance of rodent as host for coronaviruses.
Coronaviruses (CoVs) are enveloped viruses in the Coronaviridae family that contain a positive-sense and single-stranded RNA genome of approximately 30 kilobases . CoVs consist of 4 genera and have been identified in a wide range of animals and in humans. Members of the Alphacoronavirus (α-CoV) and Betacoronavirus (β-CoV) infect mammals, and members of the Gammacoronavirus (γ-CoV) and Deltacoronavirus (δ-CoV) mainly infect avian species [2,3,4]. As important etiological agents, CoVs have been recognized in human and animals and cause upper respiratory diseases in most cases. To date, 6 human CoVs were discovered: 4 of them (HCoV-229E, NL63, OC43, and HKU1) mainly cause mild respiratory diseases, and the other 2, severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) cause severe respiratory diseases [5, 6].
The SARS-CoV outbreak boosted the discovery of novel CoVs in various animals, particularly in bats. Over 140 novel bat coronaviruses (species or genotypes) have been discovered since the SARS outbreak [7, 8]. Furthermore, there is strong evidence to show that SARS-CoV, MERS-CoV, and HCoV229E may have evolved from bat CoVs [9,10,11,12,13].
Rodents are the most diverse mammals on the planet and have been documented as important carriers of human diseases . Although murine hepatitis virus (MHV) has been used as a model to study CoV for a long time, limited information is available regarding the prevalence and diversity of rodent CoVs [15,16,17,18]. Recently, several novel α-CoVs and β-CoVs (LRNV, LAMV, LRLV, and HKU24) were identified in rodents in China and Europe [19,20,21]. These discoveries suggested that rodents may carry diverse, unrecognized CoVs . In the present study, we describe the first discovery of CoVs in 3 different rodent species in the Yunnan Province of China and report a much higher (21.4%) detection rate of CoV nucleic acid in A. chevrieri than in other rodent species studied previously (<5%) [19, 20]. In addition, this is the first report of finding α-CoV and β-CoV in the same rodent species in China.
In October 2011, for pest control and routine pathogen surveillance, 177 rodents were captured in the bush and grass near the cropland ridge in Jianchuan County of the Yunnan Province (Additional file 1: Figure S1). Animal intestines were collected and transferred to liquid nitrogen for temporary preservation and transport. Following arrival at the laboratory, the samples were stored at –80 °C until they were used for virus detection. Animal species were first identified based on morphology and further by DNA sequencing of the mitochondrial cytochrome b (CytB) gene with ready-to-use methods .
To extract viral RNA, 50 mg of intestinal tissue samples were homogenized using 1 ml PBS. Homogenates were centrifuged and RNA was extracted from 140 μL supernatant using the QIAamp Viral RNA Mini Kit (Qiagen) following the manufacturer’s instructions. Extracted RNA was used as a template for amplifying the mitochondrial CytB gene with the primers CytBF (5′-ATGATATGAAAAACCATCGTTG-3′) and CytBR (5′-TTTCCNTTTCTGGTTTACAAGAC-3′). The 1.2-kb replicons were gel purified (Promega, Madison, USA) and directly sequenced using the forward and reverse primers with a 3100 Sequencer (Applied Biosystems, Waltham, MA, USA).
Reverse transcription PCR (RT-PCR) screening of CoVs
The 440-bp RNA-dependent RNA polymerase gene (RdRp) fragment of CoVs was amplified by previously described methods using a One-Step RT-PCR (Invitrogen, San Diego, USA) . PCR target bands were gel purified and sequenced on a 3100 Sequencer (ABI, Waltham, MA, USA). Standard precautions were taken to avoid PCR contamination, and no false positives were observed in the negative controls. To determine the heterogeneity of the amplicons, we inspected the sequencing chromatograms. No overlapping multicolor peaks were found, indicating that no CoV co-infection occurred in the animals examined in this study. To confirm the PCR results, positive samples were verified by performing two independent PCRs. The CoV-positive samples were named using the rodent species name, the location (Jianchuan County), and the sample number. For example, a CoV detected in A. chevrieri (sample number 54) was named as A. chevrieri CoV JC54 (AcCoV-JC54).
Three positive rodent samples representing different CoVs (JC30, β-CoV; JC34 and JC54, α-CoV) were used to perform viral isolation in Vero E6 cells (African green monkey kidney cells, ATCC: CRL-1586).
To sequence the viral genome, 140 μL supernatant from a JC34 tissue homogenate was treated using viral metagenomics procedures and ready-to-use methods . Synthesized DNA was used to construct the sequencing library, and next-generation sequencing (NGS) was performed using an HiSeq-PE150 instrument (Illumina/Solexa). BLAST searches were performed against the CoV database, and 413,599 reads homologous to CoV sequences were found. The reads were processed using Geneious (Version 5.5.9, Biomatters Limited, Auckland, New Zealand) to assemble a near full-length CoV genome contig. Subsequently, 5′ and 3′ RACE (Takara) were performed to confirm the ends of the genome sequence using two primers (5′-CAGGACGTCTAATGCAATACCT-3′ and 5′-AACACACTGAAATCAGACCTTG-3′), which were designed based on the obtained contig sequences and primers supplied in the kits. The replicons were both end sequenced. Finally, all sequences were assembled with the CoV contig to obtain a full-length CoV genome, designated as AcCoV-JC34.
The genome sequence of AcCoV-JC34 was compared to other representative CoVs with complete genomes available. The open reading frames (ORFs), deduced amino acid sequences, and potential cleavage sites in orf1ab were predicted by ORF Finder (NCBI) and ZCURVE_CoV 2.0 . Sequence alignment and editing were performed using ClustalW (Version 2.0), BioEdit (Version 7.1.9), and Geneious (Version 5.5.9) [27, 28]. The spike (S) protein structure of AcCoV-JC34 was searched against sequences in the Protein Data Bank (PDB) and predicted using a web-based SWISS-MODEL server. The cleavage sites in the S protein were predicted by comparing amino acid sequences, combined with analysis using a web-based ProP server . Phylogenetic trees were constructed using the maximum-likelihood (ML) algorithm, with bootstrap values determined by 1000 replicates in MEGA6 and PhyML software [30, 31].
Detection of α-CoVs and β-CoVs in rodents
A total of 177 intestinal samples were obtained from rodents, including three different species. By RT-PCR, targeting a 440 base pairs (bp) partial RdRp gene sequence of CoVs, 23 rodents were identified as CoV positive, which included 21 of 98 (21.4%) A. chevrieri, 1 of 17 (5.9%) A. ilex, and 1 of 62 (1.6%) E. fidelis samples (Table 1). The obtained CytB gene sequences were deposited in GenBank under accession numbers KX964655–KX964657. The isolation of rodent CoV from VeroE6 cells was not successful.
Partial RdRp sequences were searched against the CoV database, and the results indicated that 21 out of 23 sequences were β-CoVs, whereas the remaining two were α-CoVs. The β-CoV-related sequences had high nucleotide (nt) identities of 95–99% compared to the unclassified β-CoVs detected in rodents in China, Longquan-343 (A. agrarius) and HKU24 (R. norvegicus). The A. chevrieri CoV JC34 (AcCoV-JC34) and E. fidelis CoV JC54 (EfCoV-JC54) shared the highest nt identities of 84 and 98%, respectively, with the unclassified α-CoVs detected from R. norvegicus, LRNV (GenBank no: KF294380) in China or UKRn1 (GenBank no: KU739071) in Europe (Fig. 1). LRNV has a complete genome of 28,763 bp and UKRn1 (KU739071) has a partial RdRp sequence of 630 bp. Sequenced RdRp fragments in this study were deposited in GenBank under accession numbers: KX964650–KX964654.
Genome organization and ORFs of AcCoV-JC34
One positive sample (JC34) was chosen for further sequencing to obtain the full-length genome because it showed low sequence similarity to other CoVs and appeared to be a novel CoV. By random PCR and Illumina sequencing, a near full-length genome of CoV was assembled from 413,599 reads. After sequencing 5′- and 3′-rapid amplification of cDNA end replicons, a complete genome was characterized. This virus was named rodent AcCoV-JC34 and the complete genome sequence was deposited in GenBank under accession number KX964649.
The genome size of AcCoV-JC34 was 27,649 bp and the G + C content was 40%. Similar to other α-CoVs, AcCoV-JC34 has a concise genomic organization and genes characteristic of CoV, including (from 5′ to 3′) the ORF1ab, S, envelope (E), membrane (M), and nucleocapsid (N) genes (Table 2 and Fig. 2). In addition, ORFs likely coding for accessory proteins ORF3, ORF6, ORF8, and ORF9 were also found.
A hexanucleotide transcriptional regulatory sequence (TRS) is in the 5′-leader sequence and is required for the transcription of subgenomic RNAs, which is a unique characteristic of CoVs. Similar to rat CoV LRNV, bat CoV HKU2, and human CoV NL63, a putative TRS motif, 5′-AACUAA-3′ was found upstream of each ORF except for S, ORF3, and ORF8 in the AcCoV-JC34 genome. Instead of the 5′-AACUAA-3′ motif, the S, ORF3, and ORF8 genes had a variant TRS, 5′-AACUUA-3′, 5′-UACUAA-3′, and 5′-CACUAA-3′, respectively (Table 3).
Sixteen putative nonstructural proteins (nsp1 to nsp16) coded by ORF1ab of the AcCoV-JC34 were predicted (Table 4). The overall amino acid (aa) identity between the ORF1a and ORF1b polyproteins of AcCoV-JC34 and those of LRNV were 76 and 93.5%, respectively, but <60% relative to those of the other α-CoVs. The most conserved proteins 3CLpro (nsp5), RdRp (nsp12), and Hel (nsp13) of AcCoV-JC34 possessed high aa identities, ranging from 91.9 to 96% compared to those of LRNV, but possessed low aa identities ranging from 57 to 77.9% compared to those of other known α-CoVs (Table 2). In addition, similar to the normal cleavage sites found in polyprotein of α-CoVs, 10 different cleavage sites were predicted between the nsps of AcCoV-JC34 (Table 4).
The S protein of AcCoV-JC34, consisting of 1126 amino acid residues, is predicted to be a type-I membrane glycoprotein with a signal peptide (residues 1 to 19), an extracellular region (residues 20 to 1070), a transmembrane domain (residues 1071 to 1093), and an intracellular region (residues 1094 to 1126) (Additional file 1: Figure S2). A fusion peptide (FP) and two heptad repeats (HR1 and HR2) important for membrane fusion and viral entry were located at residues 674 to 692 for FP, 753 to 840 for HR1, and 1029 to 1058 for HR2. The S protein of AcCoV-JC34 showed the highest aa similarity of 66.5% compared with rat CoV LRNV, followed by 39.2% compared with BtCoV HKU2. Proteolysis of the S protein plays a pivotal role in the activation of viral and cell membrane fusion. Two cleavage sites, one located at residue 508 at the S1/S2 interface (RRAR/AR), and the other located at residue 674 (R/S) at the S2′ position, were predicted by comparing aa sequences based on analysis with a web-based ProP program (Additional file 1: Figure S2). The S1 region of AcCoV-JC34 has an N-terminal domain (NTD) and C-terminal domain (CTD). Both the NTD and CTD showed low aa sequence identities of <25% with those of other very well characterized CoVs. One of them was responsible for receptor recognition and binding, but due to the high dissimilarity with known receptor-binding domains (RBDs), it was difficult to determine the precise location of the RBD of AcCoV-JC34.
The AcCoV-JC34 proteins ORF3, E, M, ORF6, N, ORF8, and ORF9 also had low aa identities with those of other known α-CoVs. The structural proteins E, M, and N of AcCoV-JC34 showed differences compared to homologues of known CoVs. The most conserved M protein had 46.3 to 92.3% sequence identity relative to those of α-CoVs. The N protein was most variable with only 21.3 to 77.4% sequence identity compared to those of α-CoVs at the aa-sequence level. Homologues of the ORF3, ORF6, ORF8, and ORF9 proteins could be found among some CoVs but with low similarity. Previous studies have shown that the ORF3, ORF6, and ORF9 proteins of CoVs may play different functions for the viral life cycle and pathogenesis, although more studies are needed to discern the functions of the NS proteins of AcCoV-JC34.
Phylogenetic features of rodent CoVs
The first phylogenetic tree was constructed based on 400-bp RdRp sequences. In this tree, JC54 and JC34 clustered in the α-CoVs, within rodent and shrew CoVs (Fig. 1). JC34 was distantly related to the branch formed by the closely related CoV strains JC54, UKRn1, and LNRV (Lucheng-19). The other 21 CoV sequences detected from A. chevrieri or A. ilex clustered in β-CoVs and formed an independent lineage together with HKU24 from R. norvegicus and Longquan-353 from A. agrarius, in China. The 20 sequences detected from A. chevrieri were further divided into two branches.
Using predicted protein sequences, we further analyzed the phylogenetic features of AcCoV-JC34. In the phylogenetic trees constructed based on polyprotein 1a and 1b, AcCoV-JC34 clustered in the same branch with a rat CoV LRNV. Interestingly, in the tree based on the S protein, AcCoV-JC34 clustered with a rat CoV LRNV and a bat CoV HKU2 and formed a branch that appeared distinct from α-CoVs, β-CoVs, γ-CoVs, and δ-CoVs (Fig. 3). In the trees based on other genes, AcCoV-JC34 and LRNV together formed independent branches. These results indicated that AcCoV-JC34 possessed a special evolutionary position and may have a common origin with LRNV and HKU2 for the S protein (Fig. 3).
We detected CoVs in three different rodent species (A. chevrieri, A. ilex, and E. fidelis) from the Yunnan Province of China. In this study, we found a much higher (21.4%) detection rate of CoV nucleic acids in A. chevrieri than detected previously in other rodent species (<5%) [19, 20]. In addition, both α-CoV and β-CoV were found in A. chevrieri, suggesting that A. chevrieri may play an important role as a natural CoV host. A. chevrieri is known as Chevrier’s field mouse and is a dominant and endemic species in southwest China [32, 33]. In the Yunnan Province, A. chevrieri is an important pest in agriculture and human diseases that has been identified as a natural reservoir for plague bacilli and hantavirus . The detection of both α-CoV and β-CoV in A. chevrieri with high infection rates highlighted the importance of viral surveillance in A. chevrieri in the Yunnan Province, which may be helpful for disease prevention and control.
We further characterized a full-length genome of a novel α-CoV, AcCoV-JC34, from A. chevrieri. In all five conserved replicase domains, AcCoV-JC34 was the most closely related to a R. norvegicus CoV LNRV, but diverged from other α-CoVs, indicating that AcCoV-JC34 and LNRV belong to a novel α-CoV species. To our knowledge, AcCoV-JC34 is one of the few rodent α-CoVs with a complete genome.
The genome of AcCoV-JC34 had some unique features compared to other CoVs, such as a shorter nsp5 (3CLpro) and three variant TRSs. These sequences containing the genes or elements were verified by PCR and NGS. Analysis of the aa sequence showed that the proteins encoded by AcCoV-JC34 had very low similarities to other α-CoVs. In particular, the S protein sequence had <20% sequence identity to other α-CoVs (except for LNRV and HKU2), but had a little higher identity (20.6 to 22.1%) compared to β-CoVs. In addition, the N proteins normally were conserved among each of CoV genera, but the N protein of AcCoV-JC34 only shared approximately 25% aa sequence identity with other α-CoVs (except for LNRV) (Table 2). The phylogenetic trees for ORF1a, 1b, and N showed that AcCoV-JC34 and LNRV formed a distinct branch within but at the root position of α-CoVs, suggesting that AcCoV-JC34 and LNRV may represent a special evolutionary position among α-CoVs. More interestingly, in the phylogenetic trees of S, AcCoV-JC34, LNRV, and HKU2 formed a root branch including all CoVs. These results suggested that other unknown CoVs exist in rodents in nature. Further studies should be continued to reveal the prevalence, diversity, and evolution of rodent CoVs.
All samples used in this study were from rodent intestines, suggesting a possible enteric tropism of rodent CoVs. During previous decades of research, different tissue tropisms of rodent CoVs have been observed. As the prototype of rodent CoVs, different MHV strains can infect variant tissues, and the A59 strain is primarily hepatotropic, but the JHM strain is neurotropic [15,16,17,18]. RCoV and a strain of sialodacryoadenitis virus (SDAV) both primarily infect the respiratory tract . However, the tropisms of recently identified rodent CoVs from China and Europe have not been confirmed. A CoV in lineage A of β-CoV detected in the alimentary tract samples of Norway rats, HKU24, probably has enteric tropism . Another cluster of α-CoVs (PLMg1, UKMa2, UKMa1, and UKRn1) were only detected in liver samples of Norway rats, the bank vole, the wood mouse, and the noncyclic field vole, suggesting that they are hepatotropic . Additional research identified rodent α-CoV LRNV and β-CoVs LAMV and LRLV, which came from diverse tissue types that made it difficult to predict the tissue tropism of these viruses . Nonetheless, the extensively studied rodent CoVs (MHV and RCoV) could lead to severe or mild diseases in their hosts. Further studies are needed to determine the potential pathogenicity of AcCoV-JC34 along with other recently detected rodent CoVs.
In the AcCoV-JC34 genome, a predicted ORF3 protein (214 aa) was located between the S and E genes. The ORF3 protein of AcCoV-JC34 possessed 30 to 78% aa sequence identity with the homologous proteins encoded by other α-CoVs. This protein has variant names in different CoVs and was named ORF4 protein in human coronavirus 229E , non-structural protein 3 in human coronavirus NL63 , 3c-like protein or non-structural protein 3c in ferret coronavirus, 3c protein in feline coronavirus , 3b protein in transmissible gastroenteritis virus (TGEV) , and ORF3 protein in porcine epidemic diarrhea virus (PEDV) . Normally, the ORF3 protein was considered as an accessory non-structural protein, but several studies showed that the ORF3 protein was a membrane protein related to virulence [35, 37, 38, 40]. However, with low similarities between the ORF3 of AcCoV-JC34 and previously studied proteins, more experiments are needed to understand its function.
The S protein of CoVs is responsible for receptor recognition, binding, and membrane fusion, and serves as the first key factor of host restriction by meditating viral entry. In different CoVs, the RBD can be located at the NTD or CTD in S1. For example, among the α-CoVs, CTD was characterized as RBD in HCoV NL63 (aa 476–616), 229E (aa 417–547), and PEDV [41,42,43,44], but the NTD was characterized as an RBD in the TGEV . Here, the S1 of AcCoV-JC34 shared <20% aa sequence identity with those of very well characterized α-CoVs, which made it difficult to predict whether the RBD was located in NTD or CTD and which host molecule could be the possible receptor for AcCoV-JC34. The S2 of AcCoV-JC34 showed 40 to 50% identity to that of β-CoVs. By sequence alignment and SWISS-MODEL analysis (data not shown), we deduced the precise positions of FH, HR1, and HR2. The higher similarities between S2 of AcCoV-JC34 (HKU2) and β-CoVs than that to α-CoVs suggested that the structure and functional mechanism of S2 of AcCoV-JC34 may more homologous to β-CoVs.
Emerging infectious diseases caused by CoV are mostly due to interspecies transmission from animals to humans. Previous data indicated that bats are natural reservoirs for α- and β-CoVs . A number of human CoVs, including SARS-CoV, MERS-CoV, HCoV229E, and NL63 might have originated from bats [47, 48]. Among the rodent CoVs, the receptor usage, tissue tropism, and pathogenesis of MHV have been studied in detail . However, novel CoVs, like AcCoV-JC34, HKU24, LRNV, LAMV, and LRLV are not fully understood. Identification of the receptor for these viruses could help in evaluating the potential host range and ability for interspecies transmission from rodents to other mammals. Although most of these novel rodent CoVs have been characterized with full-length or near full-length genome sequences, the lack of successfully isolating those viruses thoroughly restricts future studies. More positive samples and cell lines will facilitate viral culture in the future. In addition, more attention should be paid to the diversity of CoVs in rodent, which could help to better understand the role of rodents in the evolution and ecology of CoVs.
The results of this study revealed that diverse α-CoVs and β-CoVs are circulating in rodents in the Yunnan Province of China and highlighted the importance of rodents as a natural reservoir for CoVs. The complete genome of Ac-JC34 with new characteristics and a special S gene provided new insights into the genetics and evolution of CoVs. These findings should be useful for future genomic studies of CoVs and for further functional studies of S proteins.
- A. chevrieri :
- A. ilex :
- E. fidelis :
Middle East respiratory syndrome coronavirus
Mouse hepatitis virus
Open reading frames
Protein Data Bank
Porcine epidemic diarrhea virus
RNA-dependent RNA polymerase gene
Reverse transcription PCR
Severe acute respiratory syndrome coronavirus
Transmissible gastroenteritis virus
Transcription regulatory sequence
Masters PS, Perlman S. Chapter 28: Coronaviridae. In: Knipe DM, Howley PM, editors. Philadelphia: Lippincott Williams & Wilkins; 2013. p. 825–79.
Cavanagh D, Britton P. Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses. Family Coronaviridae. In: Andrew MQ, King MJA, Carstens EB, Lefkowitz EJ, editors. International Committee on Taxonomy of Viruses. Leiden: Academic Press Elsevier; 2012. p. 770–92.
Perlman S, Netland J. Coronaviruses post-SARS: update on replication and pathogenesis. Nat Rev Microbiol. 2009;7:439–50.
Woo PC, Lau SK, Lam CS, Lau CC, Tsang AK, et al. Discovery of seven novel Mammalian and avian coronaviruses in the genus deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavirus and avian coronaviruses as the gene source of gammacoronavirus and deltacoronavirus. J Virol. 2012;86:3995–4008.
Peiris JS, Lai ST, Poon LL, Guan Y, Yam LY, et al. Coronavirus as a possible cause of severe acute respiratory syndrome. Lancet. 2003;361:1319–25.
Zaki AM, van Boheemen S, Bestebroer TM, Osterhaus AD, Fouchier RA. Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia. N Engl J Med. 2012;367:1814–20.
Ge XY, Hu B, Shi ZL. Bat coronaviruses. In: Lin-Fa Wang CC, editor. Bats and viruses: from pathogen discovery to host genomics. NY: Wiley; 2015. p. 127–55.
Drexler JF, Corman VM, Drosten C. Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS. Antiviral Res. 2014;101:45–56.
Lau SK, Woo PC, Li KS, Huang Y, Tsoi HW, et al. Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proc Natl Acad Sci U S A. 2005;102:14040–5.
Li W, Shi Z, Yu M, Ren W, Smith C, et al. Bats are natural reservoirs of SARS-like coronaviruses. Science. 2005;310:676–9.
Ge XY, Li JL, Yang XL, Chmura AA, Zhu G, et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature. 2013;503:535–8.
Corman VM, Ithete NL, Richards LR, Schoeman MC, Preiser W, et al. Rooting the phylogenetic tree of middle East respiratory syndrome coronavirus by characterization of a conspecific virus from an African bat. J Virol. 2014;88:11297–303.
Graham RL, Donaldson EF, Baric RS. A decade after SARS: strategies for controlling emerging coronaviruses. Nat Rev Microbiol. 2013;11:836–48.
Luis AD, Hayman DTS, O'Shea TJ, Cryan PM, Gilbert AT, et al. A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special? Proc R Soc B-Biol Sci. 2013;280. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3574368/.
Lane TE, Hosking MP. The pathogenesis of murine coronavirus infection of the central nervous system. Crit Rev Immunol. 2010;30:119–30.
Weiner LP. Pathogenesis of demyelination induced by a mouse hepatitis. Arch Neurol. 1973;28:298–303.
De Albuquerque N, Baig E, Ma X, Zhang J, He W, et al. Murine hepatitis virus strain 1 produces a clinically relevant model of severe acute respiratory syndrome in A/J mice. J Virol. 2006;80:10382–94.
Funk CJ, Manzer R, Miura TA, Groshong SD, Ito Y, et al. Rat respiratory coronavirus infection: replication in airway and alveolar epithelial cells and the innate immune response. J Gen Virol. 2009;90:2956–64.
Wang W, Lin XD, Guo WP, Zhou RH, Wang MR, et al. Discovery, diversity and evolution of novel coronaviruses sampled from rodents in China. Virology. 2015;474:19–27.
Lau SK, Woo PC, Li KS, Tsang AK, Fan RY, et al. Discovery of a novel coronavirus, China Rattus coronavirus HKU24, from Norway rats supports the murine origin of Betacoronavirus 1 and has implications for the ancestor of Betacoronavirus lineage A. J Virol. 2015;89:3076–92.
Tsoleridis T, Onianwa O, Horncastle E, Dayman E, Zhu M, et al. Discovery of Novel Alphacoronaviruses in European Rodents and Shrews. Viruses. 2016;8:84.
Fabre PH, Hautier L, Dimitrov D, Douzery EJ. A glimpse on the pattern of rodent diversification: a phylogenetic approach. BMC Evol Biol. 2012;12:88.
Irwin DM, Kocher TD, Wilson AC. Evolution of the Cytochrome-B Gene of Mammals. J Mol Evol. 1991;32:128–44.
de Souza Luna LK, Heiser V, Regamey N, Panning M, Drexler JF, et al. Generic detection of coronaviruses and differentiation at the prototype strain level by reverse transcription-PCR and nonfluorescent low-density microarray. J Clin Microbiol. 2007;45:1049–52.
Ge X, Li Y, Yang X, Zhang H, Zhou P, et al. Metagenomic analysis of viruses from bat fecal samples reveals many novel viruses in insectivorous bats in China. J Virol. 2012;86:4620–30.
Gao F, Ou HY, Chen LL, Zheng WX, Zhang CT. Prediction of proteinase cleavage sites in polyproteins of coronaviruses and its applications in analyzing SARS-CoV genomes. FEBS Lett. 2003;553:451–6.
Thompson JD, Gibson TJ, Higgins DG. Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. 2002;Chapter 2:Unit 2–3.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.
Duckert P, Brunak S, Blom N. Prediction of proprotein convertase cleavage sites. Protein Eng Des Sel. 2004;17:107–12.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Criscuolo A. morePhyML: improving the phylogenetic tree space exploration with PhyML 3. Mol Phylogenet Evol. 2011;61:944–8.
Xing-Yuan M, Xian-Guo G, Wen-Ge D, Ai-Qin N, Ti-Jun Q, et al. Ectoparasites of Chevrier's field mouse, Apodemus chevrieri, in a focus of plague in southwest China. Med Vet Entomol. 2007;21:297–300.
Yue H, Fan Z, Liu S, Liu Y, Song Z, et al. A mitogenome of the Chevrier's field mouse (Apodemus chevrieri) and genetic variations inferred from the cytochrome b gene. DNA Cell Biol. 2012;31:460–9.
Liang SC, Schoeb TR, Davis JK, Simecka JW, Cassell GH, et al. Comparative severity of respiratory lesions of sialodacryoadenitis virus and Sendai virus infections in LEW and F344 rats. Vet Pathol. 1995;32:661–7.
Dijkman R, Jebbink MF, Wilbrink B, Pyrc K, Zaaijer HL, et al. Human coronavirus 229E encodes a single ORF4 protein between the spike and the envelope genes. Virol J. 2006;3:106.
Dominguez SR, Sims GE, Wentworth DE, Halpin RA, Robinson CC, et al. Genomic analysis of 16 Colorado human NL63 coronaviruses identifies a new genotype, high sequence diversity in the N-terminal domain of the spike gene and evidence of recombination. J Gen Virol. 2012;93:2387–98.
Bank-Wolf BR, Stallkamp I, Wiese S, Moritz A, Tekes G, et al. Mutations of 3c and spike protein genes correlate with the occurrence of feline infectious peritonitis. Vet Microbiol. 2014;173:177–88.
Zhang X, Hasoksuz M, Spiro D, Halpin R, Wang S, et al. Complete genomic sequences, a key residue in the spike protein and deletions in nonstructural protein 3b of US strains of the virulent and attenuated coronaviruses, transmissible gastroenteritis virus and porcine respiratory coronavirus. Virology. 2007;358:424–35.
Kim YK, Cho YY, An BH, Lim SI, Lim JA, et al. Molecular characterization of the spike and ORF3 genes of porcine epidemic diarrhea virus in the Philippines. Arch Virol. 2016;161:1323–8.
O'Connor JB, Brian DA. The major product of porcine transmissible gastroenteritis coronavirus gene 3b is an integral membrane glycoprotein of 31 kDa. Virology. 1999;256:152–61.
Lin HX, Feng Y, Wong G, Wang L, Li B, et al. Identification of residues in the receptor-binding domain (RBD) of the spike protein of human coronavirus NL63 that are critical for the RBD-ACE2 receptor interaction. J Gen Virol. 2008;89:1015–24.
Wu K, Li W, Peng G, Li F. Crystal structure of NL63 respiratory coronavirus receptor-binding domain complexed with its human receptor. Proc Natl Acad Sci U S A. 2009;106:19970–4.
Bonavia A, Zelus BD, Wentworth DE, Talbot PJ, Holmes KV. Identification of a receptor-binding domain of the spike glycoprotein of human coronavirus HCoV-229E. J Virol. 2003;77:2530–8.
Deng F, Ye G, Liu Q, Navid MT, Zhong X, et al. Identification and Comparison of Receptor Binding Characteristics of the Spike Protein of Two Porcine Epidemic Diarrhea Virus Strains. Viruses. 2016;8:55.
Reguera J, Ordono D, Santiago C, Enjuanes L, Casasnovas JM. Antigenic modules in the N-terminal S1 region of the transmissible gastroenteritis virus spike protein. J Gen Virol. 2011;92:1117–26.
Woo PC, Lau SK, Lam CS, Tsang AK, Hui SW, et al. Discovery of a novel bottlenose dolphin coronavirus reveals a distinct species of marine mammal coronavirus in gammacoronavirus. J Virol. 2014;88:1318–31.
Huynh J, Li S, Yount B, Smith A, Sturges L, et al. Evidence supporting a zoonotic origin of human coronavirus strain NL63. J Virol. 2012;86:12816–25.
Pfefferle S, Oppong S, Drexler JF, Gloza-Rausch F, Ipsen A, et al. Distant relatives of severe acute respiratory syndrome coronavirus and close relatives of human coronavirus 229E in bats, Ghana. Emerg Infect Dis. 2009;15:1377–84.
Peng G, Sun D, Rajashankar KR, Qian Z, Holmes KV, et al. Crystal structure of mouse coronavirus receptor-binding domain complexed with its murine receptor. Proc Natl Acad Sci U S A. 2011;108:10696–701.
This work was jointly funded by the National Natural Science Foundation of China (81260437, 81660558), the Scientific and Technological Basis Special Project from the Minister of Science and Technology of China (2013FY113500) (to Y-ZZ), the Yunnan Centre for Collaborative Innovation in Public Health and Disease Prevention and Control (2015YNPHXT05) (to Y-ZZ), and the China Mega-Project for Infectious Disease (2014ZX10004001-003) (to Z-LS).
Availability of data and materials
Nucleotide sequences are available in GenBank under accession numbers KX964649–KX964657. Other relevant data are included in the manuscript.
YZ-Z, WH-Y, and JH-Z collected the samples. XY-G, B-L, YZ-Z, and W-Z performed the experiments. XY-G and ZL-S analyzed the data and prepared the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the Yunnan Institute of Endemic Diseases Control and Prevention. Sampling procedures were performed following the guidelines and protocols for the Laboratory Animal Use and Care from the Chinese Centers for Disease Control and the Rules for the Implementation of Laboratory Animal Medicine (1998) from the Ministry of Health, China.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Geographical map of Jianchuan country and the sampling areas. Figure S2. Alignment and predicted domains and cleavage sites of AcCoV-JC34 spike protein. NTD, N-terminal domain; RBD, receptor-binding domain; HR, heptad repeat; TM, transmembrane anchor. The signal peptide corresponds to residues 1 to 19. The cleavage sites were indicated by arrows. (PDF 223 kb)