Isolation of MERS-related coronavirus from lesser bamboo bats that uses DPP4 and infects human-DPP4-transgenic mice.

While a number of human coronaviruses are believed to be originated from ancestral viruses in bats, it remains unclear if bat coronaviruses are ready to cause direct bat-to-human transmission. Here, we report the isolation of a MERS-related coronavirus, Tylonycteris-bat-CoV-HKU4, from lesser bamboo bats. Tylonycteris-bat-CoV-HKU4 replicates efficiently in human colorectal adenocarcinoma and hepatocarcinoma cells with cytopathic effects, and can utilize human-dipeptidyl-peptidase-4 and dromedary camel-dipeptidyl-peptidase-4 as the receptors for cell entry. Flow cytometry, co-immunoprecipitation and surface plasmon resonance assays show that Tylonycteris-bat-CoV-HKU4-receptor-binding-domain can bind human-dipeptidyl-peptidase-4, dromedary camel-dipeptidyl-peptidase-4, and Tylonycteris pachypus-dipeptidyl-peptidase-4. Tylonycteris-bat-CoV-HKU4 can infect human-dipeptidyl-peptidase-4-transgenic mice by intranasal inoculation with self-limiting disease. Positive virus and inflammatory changes were detected in lungs and brains of infected mice, associated with suppression of antiviral cytokines and activation of proinflammatory cytokines and chemokines. The results suggest that MERS-related bat coronaviruses may overcome species barrier by utilizing dipeptidyl-peptidase-4 and potentially emerge in humans by direct bat-to-human transmission.

Identification of a Novel Betacoronavirus (Merbecovirus) in Amur Hedgehogs from China.

While dromedaries are the immediate animal source of Middle East Respiratory Syndrome (MERS) epidemic, viruses related to MERS coronavirus (MERS-CoV) have also been found in bats as well as hedgehogs. To elucidate the evolution of MERS-CoV-related viruses and their interspecies transmission pathway, samples were collected from different mammals in China. A novel coronavirus related to MERS-CoV, Erinaceus amurensis hedgehog coronavirus HKU31 (Ea-HedCoV HKU31), was identified from two Amur hedgehogs. Genome analysis supported that Ea-HedCoV HKU31 represents a novel species under Merbecovirus, being most closely related to Erinaceus CoV from European hedgehogs in Germany, with 79.6% genome sequence identity. Compared to other members of Merbecovirus, Ea-HedCoV HKU31 possessed unique non-structural proteins and putative cleavage sites at ORF1ab. Phylogenetic analysis showed that Ea-HedCoV HKU31 and BetaCoV Erinaceus/VMC/DEU/2012 were closely related to NeoCoV and BatCoV PREDICT from African bats in the spike region, suggesting that the latter bat viruses have arisen from recombination between CoVs from hedgehogs and bats. The predicted HKU31 receptor-binding domain (RBD) possessed only one out of 12 critical amino acid residues for binding to human dipeptidyl peptidase 4 (hDPP4), the MERS-CoV receptor. The structural modeling of the HKU31-RBD-hDPP4 binding interphase compared to that of MERS-CoV and Tylonycteris bat CoV HKU4 (Ty-BatCoV HKU4) suggested that HKU31-RBD is unlikely to bind to hDPP4. Our findings support that hedgehogs are an important reservoir of Merbecovirus, with evidence of recombination with viruses from bats. Further investigations in bats, hedgehogs and related animals are warranted to understand the evolution of MERS-CoV-related viruses.

Novel Picobirnaviruses in Respiratory and Alimentary Tracts of Cattle and Monkeys with Large Intra- and Inter-Host Diversity.

Picobirnaviruses (PBVs) are mostly found in animal alimentary samples. In this study, among 576 respiratory specimens from 476 mammals and 100 chickens, genogroup I PBVs were detected in three cattle and three monkeys, and a genogroup II PBV-positive sample was collected from one cattle specimen. More than one PBV sequence type was observed in two and one genogroup I PBV-positive samples from cattle and monkeys, respectively. Twenty-four complete/near-complete segments 2 (nine from respiratory and 15 from alimentary samples) from the cattle and monkey genogroup I PBVs and one complete segment 2 from the cattle genogroup II PBV were sequenced. Similar to other studies, the cattle PBVs also showed a high diversity. In contrast, the monkey PBVs observed in this study were clustered into three distinct clades. Within each clade, all the sequences showed >99% amino acid identities. This unique phenomenon is probably due to the fact that monkeys in our locality reside in separated troops with minimal inter-troop contact.

Bats host diverse parvoviruses as possible origin of mammalian dependoparvoviruses and source for bat-swine interspecies transmission.

Compared to the enormous species diversity of bats, relatively few parvoviruses have been reported. We detected diverse and potentially novel parvoviruses from bats in Hong Kong and mainland China. Parvoviruses belonging to Amdoparvovirus, Bocaparvovirus and Dependoparvovirus were detected in alimentary, liver and spleen samples from 16 different chiropteran species of five families by PCR. Phylogenetic analysis of partial helicase sequences showed that they potentially belonged to 25 bocaparvovirus, three dependoparvovirus and one amdoparvovirus species. Nearly complete genome sequencing confirmed the existence of at least four novel bat bocaparvovirus species (Rp-BtBoV1 and Rp-BtBoV2 from Rhinolophus pusillus, Rs-BtBoV2 from Rhinolophus sinicus and Rol-BtBoV1 from Rousettus leschenaultii) and two novel bat dependoparvovirus species (Rp-BtAAV1 from Rhinolophus pusillus and Rs-BtAAV1 from Rhinolophus sinicus). Rs-BtBoV2 was closely related to Ungulate bocaparvovirus 5 with 93, 72.1 and 78.7 % amino acid identities in the NS1, NP1 and VP1/VP2 genes, respectively. The detection of bat bocaparvoviruses, including Rs-BtBoV2, closely related to porcine bocaparvoviruses, suggests recent interspecies transmission of bocaparvoviruses between bats and swine. Moreover, Rp-BtAAV1 and Rs-BtAAV1 were most closely related to human AAV1 with 48.7 and 57.5 % amino acid identities in the rep gene. The phylogenetic relationship between BtAAVs and other mammalian AAVs suggests bats as the ancestral origin of mammalian AAVs. Furthermore, parvoviruses of the same species were detected from multiple bat species or families, supporting the ability of bat parvoviruses to cross species barriers. The results extend our knowledge on the diversity of bat parvoviruses and the role of bats in parvovirus evolution and emergence in humans and animals.

High prevalence of four novel astrovirus genotype species identified from rodents in China.

Astroviruses cause gastrointestinal and neurological infections in humans and animals. Since astrovirus is genetically diverse and different astrovirus genotypes can be found in the same animal species, astrovirus is a potential zoonotic threat to humans. In this study, we screened for astroviruses in rodents from Hong Kong, Hunan and Guangxi. Astrovirus was detected in 11.9 % (67/562) of rectal swab specimens. Phylogenetic analysis of the ORF1b region, which encodes the RdRp, showed that there were four distinct clusters (clusters A, B, C and D). Whole genome sequencing was performed for 11 representative strains from each of these four clusters. The mean amino acid genetic distances (p-dist) of full-length ORF2 were >0.634 between clusters A, B, C and other known astroviruses. The p-dist between clusters A and B, A and C, and B and C were 0.371-0.375, 0.517-0.549 and 0.524-0.555, respectively. Within cluster C, the p-dist between HN-014 and GX-006 was 0.372. Since strains with p-dist of ≥0.368 in ORF2 are now considered to be of separate genotypes species, cluster A, cluster B, cluster C-HN-014 and cluster C-GX-006 can be classified as novel genotype species. Cluster D was most closely related to the rodent astrovirus previously identified in Hong Kong. Since rodents live in close proximity to humans, interspecies jumping of these novel astroviruses may represent a threat to human health.

Identification and genomic characterization of a novel rat bocavirus from brown rats in China.

Despite recent discoveries of novel animal bocaparvoviruses, current understandings on the diversity and evolution of bocaparvoviruses are still limited. We report the identification and genome characterization of a novel bocaparvovirus, rat bocaparvovirus (RBoV), in brown rats (Rattus norvegicus) in China. RBoV was detected in 11.5%, 2.4%, 16.2% and 0.3% of alimentary, respiratory, spleen and kidney samples respectively, of 636 brown rats by PCR, but not in samples of other rodent species, suggesting that brown rats are the primary reservoir of RBoV. Six RBoV genomes sequenced from three brown rats revealed the presence of three ORFs, characteristic of bocaparvoviruses. Phylogenetic analysis showed that RBoV was distantly related to other bocaparvoviruses, forming a distinct cluster within the genus, with ≤55.5% nucleotide identities to the genome of ungulate bocaparvovirus 3, supporting its classification as a novel bocaparvovirus species. RBoV possessed a putative second exon encoding the C-terminal region of NS1 and conserved RNA splicing signals, similar to human bocaparvoviruses and canine bocaparvovirus. In contrast to human, feline and canine bocaparvoviruses which demonstrates inter/intra-host viral diversity, partial VP1/VP2 sequences of 49 RBoV strains demonstrated little inter-host genetic diversity, suggesting a single genetic group. Although the pathogenicity of RBoV remains to be determined, its presence in different host tissues suggests wide tissue tropism. RBoV represents the first bocaparvovirus in rodents with genome sequenced, which extends our knowledge on the host range of bocaparvoviruses. Further studies are required to better understand the epidemiology, genetic diversity and pathogenicity of bocaparvoviruses in different rodent populations.

High Diversity of Genogroup I Picobirnaviruses in Mammals.

In a molecular epidemiology study using 791 fecal samples collected from different terrestrial and marine mammals in Hong Kong, genogroup I picobirnaviruses (PBVs) were positive by RT-PCR targeting the partial RdRp gene in specimens from five cattle, six monkeys, 17 horses, nine pigs, one rabbit, one dog, and 12 California sea lions, with 11, 9, 23, 17, 1, 1, and 15 sequence types in the positive specimens from the corresponding animals, respectively. Phylogenetic analysis showed that the PBV sequences from each kind of animal were widely distributed in the whole tree with high diversity, sharing 47.4-89.0% nucleotide identities with other genogroup I PBV strains based on the partial RdRp gene. Nine complete segment 1 (viral loads 1.7 × 104 to 5.9 × 106/ml) and 15 segment 2 (viral loads 4.1 × 103 to 1.3 × 106/ml) of otarine PBVs from fecal samples serially collected from California sea lions were sequenced. In the two phylogenetic trees constructed using ORF2 and ORF3 of segment 1, the nine segment 1 sequences were clustered into four distinct clades (C1-C4). In the tree constructed using RdRp gene of segment 2, the 15 segment 2 sequences were clustered into nine distinct clades (R1-R9). In four sea lions, PBVs were detected in two different years, with the same segment 1 clade (C3) present in two consecutive years from one sea lion and different clades present in different years from three sea lions. A high diversity of PBVs was observed in a variety of terrestrial and marine mammals. Multiple sequence types with significant differences, representing multiple strains of PBV, were present in the majority of PBV-positive samples from different kinds of animals.

Identification of Novel Rosavirus Species That Infects Diverse Rodent Species and Causes Multisystemic Dissemination in Mouse Model.

While novel picornaviruses are being discovered in rodents, their host range and pathogenicity are largely unknown. We identified two novel picornaviruses, rosavirus B from the street rat, Norway rat, and rosavirus C from five different wild rat species (chestnut spiny rat, greater bandicoot rat, Indochinese forest rat, roof rat and Coxing's white-bellied rat) in China. Analysis of 13 complete genome sequences showed that "Rosavirus B" and "Rosavirus C" represent two potentially novel picornavirus species infecting different rodents. Though being most closely related to rosavirus A, rosavirus B and C possessed distinct protease cleavage sites and variations in Yn-Xm-AUG sequence in 5'UTR and myristylation site in VP4. Anti-rosavirus B VP1 antibodies were detected in Norway rats, whereas anti-rosavirus C VP1 and neutralizing antibodies were detected in Indochinese forest rats and Coxing's white-bellied rats. While the highest prevalence was observed in Coxing's white-bellied rats by RT-PCR, the detection of rosavirus C from different rat species suggests potential interspecies transmission. Rosavirus C isolated from 3T3 cells causes multisystemic diseases in a mouse model, with high viral loads and positive viral antigen expression in organs of infected mice after oral or intracerebral inoculation. Histological examination revealed alveolar fluid exudation, interstitial infiltration, alveolar fluid exudate and wall thickening in lungs, and hepatocyte degeneration and lymphocytic/monocytic inflammatory infiltrates with giant cell formation in liver sections of sacrificed mice. Since rosavirus A2 has been detected in fecal samples of children, further studies should elucidate the pathogenicity and emergence potential of different rosaviruses.

Comparative genome and evolutionary analysis of naturally occurring Beilong virus in brown and black rats.

Recently, we reported the presence of Beilong virus in spleen and kidney samples of brown rats and black rats, suggesting that these rodents could be natural reservoirs of Beilong virus. In this study, four genomes of Beilong virus from brown rats and black rats were sequenced. Similar to the Beilong virus genome sequenced from kidney mesangial cell line culture, those of J-virus from house mouse and Tailam virus from Sikkim rats, these four genomes from naturally occurring Beilong virus also contain the eight genes (3'-N-P/V/C-M-F-SH-TM-G-L-5'). In these four genomes, the attachment glycoprotein encoded by the G gene consists of 1046 amino acids; but for the original Beilong virus genome sequenced from kidney mesangial cell line, the G CDS was predicted to be prematurely terminated at position 2205 (TGG→TAG), resulting in a 734-amino-acid truncated G protein. This phenomenon of a lack of nonsense mutation in naturally occurring Beilong viruses was confirmed by sequencing this region of 15 additional rodent samples. Phylogenetic analyses showed that the cell line and naturally occurring Beilong viruses were closely clustered, without separation into subgroups. In addition, these viruses were further clustered with J-virus and Tailam virus, with high bootstrap supports of >90%, forming a distinct group in Paramyxoviridae. Brown rats and black rats are natural reservoirs of Beilong virus. Our results also supports that the recently proposed genus, Jeilongvirus, should encompass Beilong virus, J-virus and Tailam virus as members.

A six-year descriptive epidemiological study of human coronavirus infections in hospitalized patients in Hong Kong.

We conducted a six-year epidemiological study on human coronaviruses (HCoVs) circulating in Hong Kong, using 8275 nasopharyngeal samples from patients with acute respiratory tract infections. HCoVs were detected in 77 (0.93%) of the samples by a pan-HCoV RT-PCR assay. The most frequently detected HCoV species was HCoV-OC43 (0.58%), followed by HCoV-229E (0.15%), HCoV-HKU1 (0.13%) and HCoV-NL63 (0.07%). HCoVs were detected throughout the study period (September 2008-August 2014), with the highest detection rate from September 2010 to August 2011 (22/1500, 1.47%). Different seasonal patterns of each HCoV species in Hong Kong were noted. HCoV-OC43 was predominant in the fall and winter, whereas HCoV-HKU1 showed peak activity in winter, with a few cases occurred in spring and summer. HCoV-229E mainly occurred in winter and spring, while HCoV-NL63 was predominant in summer and autumn. HCoVs most commonly infect the elderly and young children, with median age of 79.5 years (range, 22 days to 95 years). Intriguingly, the detection rate of HCoV-OC43 in the age group of > 80 years (26/2380, 1.09%) was significantly higher than that in the age group of 0-10 years (12/2529, 0.47%) (P < 0.05). These data provides new insight into the epidemiology of coronaviruses.

Severe Acute Respiratory Syndrome (SARS) Coronavirus ORF8 Protein Is Acquired from SARS-Related Coronavirus from Greater Horseshoe Bats through Recombination.

UNLABELLED: Despite the identification of horseshoe bats as the reservoir of severe acute respiratory syndrome (SARS)-related coronaviruses (SARSr-CoVs), the origin of SARS-CoV ORF8, which contains the 29-nucleotide signature deletion among human strains, remains obscure. Although two SARS-related Rhinolophus sinicus bat CoVs (SARSr-Rs-BatCoVs) previously detected in Chinese horseshoe bats (Rhinolophus sinicus) in Yunnan, RsSHC014 and Rs3367, possessed 95% genome identities to human and civet SARSr-CoVs, their ORF8 protein exhibited only 32.2 to 33% amino acid identities to that of human/civet SARSr-CoVs. To elucidate the origin of SARS-CoV ORF8, we sampled 348 bats of various species in Yunnan, among which diverse alphacoronaviruses and betacoronaviruses, including potentially novel CoVs, were identified, with some showing potential interspecies transmission. The genomes of two betacoronaviruses, SARSr-Rf-BatCoV YNLF_31C and YNLF_34C, from greater horseshoe bats (Rhinolophus ferrumequinum), possessed 93% nucleotide identities to human/civet SARSr-CoV genomes. Although these two betacoronaviruses displayed lower similarities than SARSr-Rs-BatCoV RsSHC014 and Rs3367 in S protein to civet SARSr-CoVs, their ORF8 proteins demonstrated exceptionally high (80.4 to 81.3%) amino acid identities to that of human/civet SARSr-CoVs, compared to SARSr-BatCoVs from other horseshoe bats (23.2 to 37.3%). Potential recombination events were identified around ORF8 between SARSr-Rf-BatCoVs and SARSr-Rs-BatCoVs, leading to the generation of civet SARSr-CoVs. The expression of ORF8 subgenomic mRNA suggested that the ORF8 protein may be functional in SARSr-Rf-BatCoVs. The high Ka/Ks ratio among human SARS-CoVs compared to that among SARSr-BatCoVs supported that ORF8 is under strong positive selection during animal-to-human transmission. Molecular clock analysis using ORF1ab showed that SARSr-Rf-BatCoV YNLF_31C and YNLF_34C diverged from civet/human SARSr-CoVs in approximately 1990. SARS-CoV ORF8 originated from SARSr-CoVs of greater horseshoe bats through recombination, which may be important for animal-to-human transmission. IMPORTANCE: Although horseshoe bats are the primary reservoir of SARS-related coronaviruses (SARSr-CoVs), it is still unclear how these bat viruses have evolved to cross the species barrier to infect civets and humans. Most human SARS-CoV epidemic strains contain a signature 29-nucleotide deletion in ORF8, compared to civet SARSr-CoVs, suggesting that ORF8 may be important for interspecies transmission. However, the origin of SARS-CoV ORF8 remains obscure. In particular, SARSr-Rs-BatCoVs from Chinese horseshoe bats (Rhinolophus sinicus) exhibited <40% amino acid identities to human/civet SARS-CoV in the ORF8 protein. We detected diverse alphacoronaviruses and betacoronaviruses among various bat species in Yunnan, China, including two SARSr-Rf-BatCoVs from greater horseshoe bats that possessed ORF8 proteins with exceptionally high amino acid identities to that of human/civet SARSr-CoVs. We demonstrated recombination events around ORF8 between SARSr-Rf-BatCoVs and SARSr-Rs-BatCoVs, leading to the generation of civet SARSr-CoVs. Our findings offer insight into the evolutionary origin of SARS-CoV ORF8 protein, which was likely acquired from SARSr-CoVs of greater horseshoe bats through recombination.

Discovery of a novel bottlenose dolphin coronavirus reveals a distinct species of marine mammal coronavirus in Gammacoronavirus.

While gammacoronaviruses mainly comprise infectious bronchitis virus (IBV) and its closely related bird coronaviruses (CoVs), the only mammalian gammacoronavirus was discovered from a white beluga whale (beluga whale CoV [BWCoV] SW1) in 2008. In this study, we discovered a novel gammacoronavirus from fecal samples from three Indo-Pacific bottlenose dolphins (Tursiops aduncus), which we named bottlenose dolphin CoV (BdCoV) HKU22. All the three BdCoV HKU22-positive samples were collected on the same date, suggesting a cluster of infection, with viral loads of 1 × 10(3) to 1 × 10(5) copies per ml. Clearance of virus was associated with a specific antibody response against the nucleocapsid of BdCoV HKU22. Complete genome sequencing and comparative genome analysis showed that BdCoV HKU22 and BWCoV SW1 have similar genome characteristics and structures. Their genome size is about 32,000 nucleotides, the largest among all CoVs, as a result of multiple unique open reading frames (NS5a, NS5b, NS5c, NS6, NS7, NS8, NS9, and NS10) between their membrane (M) and nucleocapsid (N) protein genes. Although comparative genome analysis showed that BdCoV HKU22 and BWCoV SW1 should belong to the same species, a major difference was observed in the proteins encoded by their spike (S) genes, which showed only 74.3 to 74.7% amino acid identities. The high ratios of the number of synonymous substitutions per synonymous site (Ks) to the number of nonsynonymous substitutions per nonsynonymous site (Ka) in multiple regions of the genome, especially the S gene (Ka/Ks ratio, 2.5), indicated that BdCoV HKU22 may be evolving rapidly, supporting a recent transmission event to the bottlenose dolphins. We propose a distinct species, Cetacean coronavirus, in Gammacoronavirus, to include BdCoV HKU22 and BWCoV SW1, whereas IBV and its closely related bird CoVs represent another species, Avian coronavirus, in Gammacoronavirus.

Genetic characterization of Betacoronavirus lineage C viruses in bats reveals marked sequence divergence in the spike protein of pipistrellus bat coronavirus HKU5 in Japanese pipistrelle: implications for the origin of the novel Middle East respiratory syndrome coronavirus.

While the novel Middle East respiratory syndrome coronavirus (MERS-CoV) is closely related to Tylonycteris bat CoV HKU4 (Ty-BatCoV HKU4) and Pipistrellus bat CoV HKU5 (Pi-BatCoV HKU5) in bats from Hong Kong, and other potential lineage C betacoronaviruses in bats from Africa, Europe, and America, its animal origin remains obscure. To better understand the role of bats in its origin, we examined the molecular epidemiology and evolution of lineage C betacoronaviruses among bats. Ty-BatCoV HKU4 and Pi-BatCoV HKU5 were detected in 29% and 25% of alimentary samples from lesser bamboo bat (Tylonycteris pachypus) and Japanese pipistrelle (Pipistrellus abramus), respectively. Sequencing of their RNA polymerase (RdRp), spike (S), and nucleocapsid (N) genes revealed that MERS-CoV is more closely related to Pi-BatCoV HKU5 in RdRp (92.1% to 92.3% amino acid [aa] identity) but is more closely related to Ty-BatCoV HKU4 in S (66.8% to 67.4% aa identity) and N (71.9% to 72.3% aa identity). Although both viruses were under purifying selection, the S of Pi-BatCoV HKU5 displayed marked sequence polymorphisms and more positively selected sites than that of Ty-BatCoV HKU4, suggesting that Pi-BatCoV HKU5 may generate variants to occupy new ecological niches along with its host in diverse habitats. Molecular clock analysis showed that they diverged from a common ancestor with MERS-CoV at least several centuries ago. Although MERS-CoV may have diverged from potential lineage C betacoronaviruses in European bats more recently, these bat viruses were unlikely to be the direct ancestor of MERS-CoV. Intensive surveillance for lineage C betaCoVs in Pipistrellus and related bats with diverse habitats and other animals in the Middle East may fill the evolutionary gap.

Complete genome sequences of novel rat noroviruses in Hong Kong.

We report two genome sequences of novel noroviruses isolated from fecal swab specimens of brown rats in Hong Kong. The complete genome is approximately 7.5 kb in length and consists of 3 overlapping open reading frames encoding ORF1 polyprotein, VP1, and VP2, respectively. Sequence analysis suggested that these noroviruses should be classified in genogroup V, but they are distinct from other known rodent noroviruses and represent a novel cluster within the genogroup.

Recent transmission of a novel alphacoronavirus, bat coronavirus HKU10, from Leschenault's rousettes to pomona leaf-nosed bats: first evidence of interspecies transmission of coronavirus between bats of different suborders.

Although coronaviruses are known to infect various animals by adapting to new hosts, interspecies transmission events are still poorly understood. During a surveillance study from 2005 to 2010, a novel alphacoronavirus, BatCoV HKU10, was detected in two very different bat species, Ro-BatCoV HKU10 in Leschenault's rousettes (Rousettus leschenaulti) (fruit bats in the suborder Megachiroptera) in Guangdong and Hi-BatCoV HKU10 in Pomona leaf-nosed bats (Hipposideros pomona) (insectivorous bats in the suborder Microchiroptera) in Hong Kong. Although infected bats appeared to be healthy, Pomona leaf-nosed bats carrying Hi-BatCoV HKU10 had lower body weights than uninfected bats. To investigate possible interspecies transmission between the two bat species, the complete genomes of two Ro-BatCoV HKU10 and six Hi-BatCoV HKU10 strains were sequenced. Genome and phylogenetic analyses showed that Ro-BatCoV HKU10 and Hi-BatCoV HKU10 represented a novel alphacoronavirus species, sharing highly similar genomes except in the genes encoding spike proteins, which had only 60.5% amino acid identities. Evolution of the spike protein was also rapid in Hi-BatCoV HKU10 strains from 2005 to 2006 but stabilized thereafter. Molecular-clock analysis dated the most recent common ancestor of all BatCoV HKU10 strains to 1959 (highest posterior density regions at 95% [HPDs], 1886 to 2002) and that of Hi-BatCoV HKU10 to 1986 (HPDs, 1956 to 2004). The data suggested recent interspecies transmission from Leschenault's rousettes to Pomona leaf-nosed bats in southern China. Notably, the rapid adaptive genetic change in BatCoV HKU10 spike protein by ~40% amino acid divergence after recent interspecies transmission was even greater than the ~20% amino acid divergence between spike proteins of severe acute respiratory syndrome-related Rhinolophus bat coronavirus (SARSr-CoV) in bats and civets. This study provided the first evidence for interspecies transmission of coronavirus between bats of different suborders.

Isolation and characterization of a novel Betacoronavirus subgroup A coronavirus, rabbit coronavirus HKU14, from domestic rabbits.

We describe the isolation and characterization of a novel Betacoronavirus subgroup A coronavirus, rabbit coronavirus HKU14 (RbCoV HKU14), from domestic rabbits. The virus was detected in 11 (8.1%) of 136 rabbit fecal samples by reverse transcriptase PCR (RT-PCR), with a viral load of up to 10(8) copies/ml. RbCoV HKU14 was able to replicate in HRT-18G and RK13 cells with cytopathic effects. Northern blotting confirmed the production of subgenomic mRNAs coding for the HE, S, NS5a, E, M, and N proteins. Subgenomic mRNA analysis revealed a transcription regulatory sequence, 5'-UCUAAAC-3'. Phylogenetic analysis showed that RbCoV HKU14 formed a distinct branch among Betacoronavirus subgroup A coronaviruses, being most closely related to but separate from the species Betacoronavirus 1. A comparison of the conserved replicase domains showed that RbCoV HKU14 possessed <90% amino acid identities to most members of Betacoronavirus 1 in ADP-ribose 1″-phosphatase (ADRP) and nidoviral uridylate-specific endoribonuclease (NendoU), indicating that RbCoV HKU14 should represent a separate species. RbCoV HKU14 also possessed genomic features distinct from those of other Betacoronavirus subgroup A coronaviruses, including a unique NS2a region with a variable number of small open reading frames (ORFs). Recombination analysis revealed possible recombination events during the evolution of RbCoV HKU14 and members of Betacoronavirus 1, which may have occurred during cross-species transmission. Molecular clock analysis using RNA-dependent RNA polymerase (RdRp) genes dated the most recent common ancestor of RbCoV HKU14 to around 2002, suggesting that this virus has emerged relatively recently. Antibody against RbCoV was detected in 20 (67%) of 30 rabbit sera tested by an N-protein-based Western blot assay, whereas neutralizing antibody was detected in 1 of these 20 rabbits.

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.

Recently, we reported the discovery of three novel coronaviruses, bulbul coronavirus HKU11, thrush coronavirus HKU12, and munia coronavirus HKU13, which were identified as representatives of a novel genus, Deltacoronavirus, in the subfamily Coronavirinae. In this territory-wide molecular epidemiology study involving 3,137 mammals and 3,298 birds, we discovered seven additional novel deltacoronaviruses in pigs and birds, which we named porcine coronavirus HKU15, white-eye coronavirus HKU16, sparrow coronavirus HKU17, magpie robin coronavirus HKU18, night heron coronavirus HKU19, wigeon coronavirus HKU20, and common moorhen coronavirus HKU21. Complete genome sequencing and comparative genome analysis showed that the avian and mammalian deltacoronaviruses have similar genome characteristics and structures. They all have relatively small genomes (25.421 to 26.674 kb), the smallest among all coronaviruses. They all have a single papain-like protease domain in the nsp3 gene; an accessory gene, NS6 open reading frame (ORF), located between the M and N genes; and a variable number of accessory genes (up to four) downstream of the N gene. Moreover, they all have the same putative transcription regulatory sequence of ACACCA. Molecular clock analysis showed that the most recent common ancestor of all coronaviruses was estimated at approximately 8100 BC, and those of Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Deltacoronavirus were at approximately 2400 BC, 3300 BC, 2800 BC, and 3000 BC, respectively. From our studies, it appears that bats and birds, the warm blooded flying vertebrates, are ideal hosts for the coronavirus gene source, bats for Alphacoronavirus and Betacoronavirus and birds for Gammacoronavirus and Deltacoronavirus, to fuel coronavirus evolution and dissemination.

Complete genome analysis of three novel picornaviruses from diverse bat species.

Although bats are important reservoirs of diverse viruses that can cause human epidemics, little is known about the presence of picornaviruses in these flying mammals. Among 1,108 bats of 18 species studied, three novel picornaviruses (groups 1, 2, and 3) were identified from alimentary specimens of 12 bats from five species and four genera. Two complete genomes, each from the three picornaviruses, were sequenced. Phylogenetic analysis showed that they fell into three distinct clusters in the Picornaviridae family, with low homologies to known picornaviruses, especially in leader and 2A proteins. Moreover, group 1 and 2 viruses are more closely related to each other than to group 3 viruses, which exhibit genome features distinct from those of the former two virus groups. In particular, the group 3 virus genome contains the shortest leader protein within Picornaviridae, a putative type I internal ribosome entry site (IRES) in the 5'-untranslated region instead of the type IV IRES found in group 1 and 2 viruses, one instead of two GXCG motifs in 2A, an L→V substitution in the DDLXQ motif in 2C helicase, and a conserved GXH motif in 3C protease. Group 1 and 2 viruses are unique among picornaviruses in having AMH instead of the GXH motif in 3C(pro). These findings suggest that the three picornaviruses belong to two novel genera in the Picornaviridae family. This report describes the discovery and complete genome analysis of three picornaviruses in bats, and their presence in diverse bat genera/species suggests the ability to cross the species barrier.