|
|
||||||||
Short Communication |
1 State Key Laboratory of Emerging Infectious Diseases, Department of Microbiology, The University of Hong Kong, Hong Kong SAR
2 HKU-Pasteur Research Centre, Hong Kong SAR
Correspondence
J. S. M. Peiris
malik{at}hkucc.hku.hk
L. L. M. Poon
llmpoon{at}hkucc.hku.hk
| ABSTRACT |
|---|
|
|
|---|
The GenBank/EMBL/DDBJ accession numbers for the sequences reported in this paper are EU420137–EU420139.
Supplementary material is available with the online version of this paper.
| MAIN TEXT |
|---|
|
|
|---|
We previously reported the discovery of the first BtCoV (Poon et al., 2005
). In addition, we also performed a longitudinal study on Miniopterus magnater and Miniopterus pusillus bats that cohabit an abandoned cave in Hong Kong (Chu et al., 2006
). BtCoV HKU7 and bat-SARS-CoV were detected only once in our specimens, while BtCoVs 1A, 1B and HKU8 were detected repeatedly in these bats. A clear host species restriction was observed with BtCoV 1A being found in M. magnater and BtCoV 1B in M. pusillus even though these two CoVs are genetically very similar and the two bat species co-inhabit the same cave. Thus, these virus surveillance studies might allow unique opportunities to understand better the dynamics and prevalence of coronaviruses within a single geographical location. To understand better the genetic relationships of the BtCoVs that are repeatedly found in the local Miniopterus spp, representative faecal samples positive for BtCoVs 1A, 1B and HKU8 in our previous investigations were selected for this study. Genomic sequences of these three viruses were determined and analysed as described in Supplementary Material (available in JGV Online).
The genomes of BtCoVs 1A, 1B and HKU8, excluding the 3' poly(A) tails, have sizes of 28309, 28460 and 28757 nt, respectively (GenBank accession nos EU420137
[GenBank]
–EU420139
[GenBank]
). The G+C content of 1A, 1B and HKU8 is 38.2, 38.6 and 41.8 %, respectively. The genome organizations of these CoVs are typical of group 1 CoVs (Table 1
) and the amino acid sequence identities of the major viral proteins between these viruses are shown in Supplementary Table S1 (available in JGV Online). The 3' UTRs and N genes of these viruses were reported previously and will not be discussed further (Chu et al., 2006
). Interestingly, HKU8 was found to contain an extra ORF at its 3' end. This predicted ORF (designated ORF 7) overlaps with the N gene and encodes a putative protein with 248 aa residues (Table 1
).
|
These ORFs 1 occupy approximately 70 % of the genomes. ORF1 consists of ORF1a and ORF1b that encodes viral pp1a and pp1ab, respectively. Putative features responsible for the –1 ribosomal frame shifting, e.g. the slippage sequence (5'-UUUAAAC-3'), were identified in the genomes (data not shown). Based on in silico analysis, the pp1a and pp1ab proteins are predicted to be cleaved into a total of 16 non-structural proteins by virus proteases (Table 1
) (Chen et al., 2003
). The predicted 3CL protease cleavage sites in these polyproteins (Supplementary Table S2 available in JGV Online) agree with previously defined substrate specificity for other 3CL proteases (Ziebuhr et al., 2000
). All pp1ab of BtCoV 1A, 1B and HKU8 contain functional units that are typical for CoVs. For example, all these viruses contain two papain-like (PL) protease domains in nsp3, a 3C-like (3CL) protease in nsp5, CoV RNA-dependent RNA polymerase (RdRp) in nsp12, RNA helicase (Hel) in nsp13 and an mRNA cap-1 methyltransferase (MT) functional domain in nsp16. It was noted that the putative nsp1 of HKU8 has a size of 18 aa, which is significantly smaller than sizes of nsp1 predicted in other group 1 CoVs.
Similar to other CoVs, the N-terminal regions of these S genes are genetically diverse. The protein sequence identities between the putative S1 domains (the first 740 residues approximately) of HKU8 to 1A, HKU8 to 1B, and 1A to 1B are 38.5, 38.7 and 86.6 %, respectively. The corresponding sequence identities for the putative S2 domains are 56.2 (HKU8 to 1A), 55.9 (HKU8 to 1B) and 91.3 % (1A to 1B), respectively. Twenty nine, 31 and 30 potential N-glycosylation sites were predicted in 1A, 1B and HKU8 S proteins, respectively. Two putative heptad repeat (HR) regions, which are conserved in all CoVs (Bosch et al., 2003
), were identified in these S proteins (both 1A and 1B: HR1, aa 973–1122; HR2, aa 1252–1311. HKU8: HR1, aa 976–1121; HR2, aa 1261–1317).
ORF3 of 1A and 1B have the same size of 660 nt and that of HKU8 is 669 nt. In terms of amino acid sequence similarity, ORF3 is the least conserved gene amongst these three viruses (Supplementary Table S3 available in JGV Online). Signal anchors in these ORF3 proteins were predicted with the probability of 0.545–0.689 by SignalP. Multiple transmembrane motifs were also predicted in these ORF3 proteins, suggesting these proteins might be surface proteins.
The E genes of these viruses are both the same size of 225 nt. One transmembrane signal was predicted in each of these E proteins (P>0.990), with the N terminus exposed to the surface.
Signal peptides were predicted in these M proteins (P>0.933). Cleavage sites were located between aa 17 and 18 in these M proteins. The N-terminal signal peptide was not observed in other BtCoV M proteins, including BtCoV/512/2005 and HKU2. Three strong transmembrane helices and one marginal transmembrane helix were predicted with the N terminus exposed to the surface. It has been reported that the M protein of the transmissible gastroenteritis virus (TGEV) adopts two different topologies in the virus envelope with three transmembrane helices and the C terminus located inside the virus particle, or with four transmembrane helices and the C terminus exposed to the virion surface (Escors et al., 2001
). This may also be the case for M proteins of these group 1 BtCoVs.
ORF7 is present only in HKU8, and it is located downstream of the N protein. A BLAST search of the ORF7 amino acid sequence revealed no homologous gene in GenBank. Prediction of an ORF (or ORFs) downstream of the N gene was also reported in the genomes of some group 1 CoVs, including TGEV (GenBank accession no. NC_002306
[GenBank]
), feline infectious peritonitis virus (FIPV) (GenBank accession no. AY994055
[GenBank]
), BtCoV/512/2005 (Tang et al., 2006
) and BtCoV HKU2 (GenBank accession no. NC_009988
[GenBank]
). ORF7 in TGEV was shown to produce an endoplasmic reticulum associated protein, suggesting it plays a role in the association of replication complexes or in virion assembly (Tung et al., 1992
). The ORF7ab of FIPV has been suggested to affect virulence (Haijema et al., 2004
). However, there is no significant sequence similarity between these ORFs and it is likely that these viruses have acquired their unique ORFs from different sources.
The putative core sequence 5'-CUAAAC-3' in the transcription regulatory sequence (TRS), which was identified in other group 1 CoVs (Dye & Siddell, 2005
; Hiscox et al., 1995
; Lau et al., 2007
; Pyrc et al., 2004
), is present in the leader TRSs of the studied genomes. They start at position 67 in BtCoV 1A and position 68 in both BtCoV 1B and HKU8. All the core TRS sequences are highly similar to its corresponding leader TRS (Supplementary Fig. S2) except the TRSs for all ORF3 and the E gene of HKU8. It has been reported that the leader TRS can also join the body segment at an unusual consensus sequence preceding the ORFs without the presence of the core sequence (Hofmann et al., 1993
). However, we were not able to demonstrate this in these bat viruses because none of these viruses could be cultured in vitro and we have no means of obtaining sufficient viral RNA for detailed investigation. Thus, the authentic TRSs of these bat viruses, in particular those for ORF3, are yet to be confirmed.
Phylogenetic analyses of RdRp, helicase, S, E, M and N proteins from these BtCoVs were performed. Mouse hepatitis virus was included as an out-group in the analyses. In all cases, the three viruses studied formed a clade that is distinct to group 1 members isolated from other animals (Fig. 1
). The BtCoVs 1A, 1B and HKU8 viruses phylogenetically cluster with BtCoVs Shandong/977/2006 and Fujian/773/2005 detected in mainland China (Fig. 1
). Sequence analysis on partial RdRp sequences revealed that the Fujian/773/2005 found in Miniopterus schreibersii could be a HKU7-like virus (not shown). The previously published S gene and helicase gene sequences of BtCoV Shandong/977/2006 (Vijaykrishna et al., 2007
), which was detected in bats of the species Rhinolophus ferrumequinum, show a close genetic relationship to the three CoVs described in this report. The S protein amino acid sequence similarity between BtCoV Shandong/977/2006 and BtCoV 1A is 46.5 %, while BtCoV 1A to HKU8 is 46.4 % (Supplementary Table S3). A similar relationship was also observed for the helicase gene (Shandong/977/2006 to 1A: 74.2 %; Shandong/977/2006 to HKU8: 75.5 %; 1A to HKU8: 76.5 %). Therefore, CoVs related to that found in the genus Miniopterus could also be found in other genera of bats. Considering the relationship between these CoV S proteins, which is one of the major determinants of host species specificity of infection (Dveksler et al., 1991
), and the relatively low detection rates of BtCoVs HKU8 and HKU7 in Miniopterus (Chu et al., 2006
), it is possible that the two viruses are naturally harboured in species of bats other than Miniopterus bats. Overall, the above findings highlight that some BtCoVs might have the ability for interspecies transmission, an event that is relevant to the genesis of SARS CoVs in humans.
|
The limitation of direct sequencing of PCR amplicons is that the presence of co-infections may be masked by the sequence of the dominant virus. Given the high detection rates of BtCoVs 1A and 1B in our specimens and by the apparent species restriction of these two viruses to two different bat species that co-inhabit the same cave, we designed species-specific PCR assays for detecting BtCoVs 1A, 1B and HKU8 in our specimens. We applied these assays to all specimens that were positive in our CoV consensus screening RT-PCR in order to detect the presence of co-infections. These viruses were selected for screening as they were repeatedly detected in Miniopterus bats found within the same cave in Hong Kong. A total of 38 specimens collected in a period from 2004 to 2006 in Hong Kong, which were known to be positive for CoV (Chu et al., 2006
), were tested by six different species-specific RT-PCR assays (RdRp and N genes of BtCoV 1A, 1B and HKU8). The identities of the amplified products were all confirmed by sequencing. In summary, infection with 1A alone was detected in 16 specimens collected from M. magnater, infection with 1B alone was detected in 14 specimens collected from M. pusillus and infection with BtCoV HKU8 alone was detected in one specimen collected from M. magnater. There was no evidence of co-infections with the more genetically closely related viruses 1A and 1B, even though these two viruses co-habited the same cave and were more frequently detected in our specimens from their corresponding hosts. These results suggest that there is a strong host restriction for these two viruses that apparently hinders interspecies transmission between these two species of bat. Interestingly, co-infections of 1B and HKU8 were detected in seven M. pusillus specimens collected in 2004 (n=2) and 2006 (n=5) (Table 2
).
|
In this report, we describe the full genomes of three group 1 BtCoVs. Co-infection with different CoVs in their host provides chances for producing recombinants between these CoVs. In our study, there was no evidence of recombination events detected between the reported genomes of BtCoV 1A, 1B and HKU8. But the diversity of CoVs found in bats and the detection of BtCoV co-infections suggest bats might play a key role in shaping the evolution and ecology of CoV. The availability of genetic information for these viruses helps us better understand the taxonomy and evolution of CoVs in general. But more importantly, since SARS and a number of other new human diseases have emerged through interspecies transmission of viruses harboured by bats, the biological basis of host restriction of bat CoVs remains an important subject for study for global public health.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Chen, L. L., Ou, H. Y., Zhang, R. & Zhang, C. T. (2003). ZCURVE_CoV: a new system to recognize protein coding genes in coronavirus genomes, and its applications in analyzing SARS-CoV genomes. Biochem Biophys Res Commun 307, 382–388.[CrossRef][Medline]
Chu, D. K., Poon, L. L., Chan, K. H., Chen, H., Guan, Y., Yuen, K. Y. & Peiris, J. S. (2006). Coronaviruses in bent-winged bats (Miniopterus spp.). J Gen Virol 87, 2461–2466.
Dominguez, S. R., O'Shea, T. J., Oko, L. M. & Holmes, K. V. (2007). Detection of group 1 coronaviruses in bats in North America. Emerg Infect Dis 13, 1295–1300.[Medline]
Dong, B. Q., Liu, W., Fan, X. H., Vijaykrishna, D., Tang, X. C., Gao, F., Li, L. F., Li, G. J., Zhang, J. X. & other authors (2007). Detection of a novel and highly divergent coronavirus from Asian leopard cats and Chinese ferret badgers in Southern China. J Virol 81, 6920–6926.
Dveksler, G. S., Pensiero, M. N., Cardellichio, C. B., Williams, R. K., Jiang, G. S., Holmes, K. V. & Dieffenbach, C. W. (1991). Cloning of the mouse hepatitis virus (MHV) receptor: expression in human and hamster cell lines confers susceptibility to MHV. J Virol 65, 6881–6891.
Dye, C. & Siddell, S. G. (2005). Genomic RNA sequence of Feline coronavirus strain FIPV WSU-79/1146. J Gen Virol 86, 2249–2253.
Escors, D., Camafeita, E., Ortego, J., Laude, H. & Enjuanes, L. (2001). Organization of two transmissible gastroenteritis coronavirus membrane protein topologies within the virion and core. J Virol 75, 12228–12240.
Gorbalenya, A. E., Enjuanes, L., Ziebuhr, J. & Snijder, E. J. (2006). Nidovirales: evolving the largest RNA virus genome. Virus Res 117, 17–37.[CrossRef][Medline]
Haijema, B. J., Volders, H. & Rottier, P. J. (2004). Live, attenuated coronavirus vaccines through the directed deletion of group-specific genes provide protection against feline infectious peritonitis. J Virol 78, 3863–3871.
Hiscox, J. A., Mawditt, K. L., Cavanagh, D. & Britton, P. (1995). Investigation of the control of coronavirus subgenomic mRNA transcription by using T7-generated negative-sense RNA transcripts. J Virol 69, 6219–6227.
Hofmann, M. A., Chang, R. Y., Ku, S. & Brian, D. A. (1993). Leader-mRNA junction sequences are unique for each subgenomic mRNA species in the bovine coronavirus and remain so throughout persistent infection. Virology 196, 163–171.[CrossRef][Medline]
Lau, S. K., Woo, P. C., Li, K. S., Huang, Y., Tsoi, H. W., Wong, B. H., Wong, S. S., Leung, S. Y., Chan, K. H. & Yuen, K. Y. (2005). Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proc Natl Acad Sci U S A 102, 14040–14045.
Lau, S. K., Woo, P. C., Li, K. S., Huang, Y., Wang, M., Lam, C. S., Xu, H., Guo, R., Chan, K. H. & other authors (2007). Complete genome sequence of bat coronavirus HKU2 from Chinese horseshoe bats revealed a much smaller spike gene with a different evolutionary lineage from the rest of the genome. Virology 367, 428–439.[CrossRef][Medline]
Liu, P., Li, L., Millership, J. J., Kang, H., Leibowitz, J. L. & Giedroc, D. P. (2007). A U-turn motif-containing stem-loop in the coronavirus 5' untranslated region plays a functional role in replication. RNA 13, 763–780.
Makino, S., Keck, J. G., Stohlman, S. A. & Lai, M. M. (1986). High-frequency RNA recombination of murine coronaviruses. J Virol 57, 729–737.
Muller, M. A., Paweska, J. T., Leman, P. A., Drosten, C., Grywna, K., Kemp, A., Braack, L., Sonnenberg, K., Niedrig, M. & Swanepoel, R. (2007). Coronavirus antibodies in African bat species. Emerg Infect Dis 13, 1367–1370.[Medline]
Peiris, J. S., Lai, S. T., Poon, L. L., Guan, Y., Yam, L. Y., Lim, W., Nicholls, J., Yee, W. K., Yan, W. W. & other authors (2003). Coronavirus as a possible cause of severe acute respiratory syndrome. Lancet 361, 1319–1325.[CrossRef][Medline]
Poon, L. L., Chu, D. K., Chan, K. H., Wong, O. K., Ellis, T. M., Leung, Y. H., Lau, S. K., Woo, P. C., Suen, K. Y. & other authors (2005). Identification of a novel coronavirus in bats. J Virol 79, 2001–2009.
Pyrc, K., Jebbink, M. F., Berkhout, B. & van der Hoek, L. (2004). Genome structure and transcriptional regulation of human coronavirus NL63. Virol J 1, 7[CrossRef][Medline]
Tang, X. C., Zhang, J. X., Zhang, S. Y., Wang, P., Fan, X. H., Li, L. F., Li, G., Dong, B. Q., Liu, W. & other authors (2006). Prevalence and genetic diversity of coronaviruses in bats from China. J Virol 80, 7481–7490.
Tung, F. Y., Abraham, S., Sethna, M., Hung, S. L., Sethna, P., Hogue, B. G. & Brian, D. A. (1992). The 9-kDa hydrophobic protein encoded at the 3' end of the porcine transmissible gastroenteritis coronavirus genome is membrane-associated. Virology 186, 676–683.[CrossRef][Medline]
Vijaykrishna, D., Smith, G. J., Zhang, J. X., Peiris, J. S., Chen, H. & Guan, Y. (2007). Evolutionary insights into the ecology of coronaviruses. J Virol 81, 4012–4020.
Woo, P. C., Lau, S. K., Li, K. S., Poon, R. W., Wong, B. H., Tsoi, H. W., Yip, B. C., Huang, Y., Chan, K. H. & Yuen, K. Y. (2006). Molecular diversity of coronaviruses in bats. Virology 351, 180–187.[CrossRef][Medline]
Ziebuhr, J., Snijder, E. J. & Gorbalenya, A. E. (2000). Virus-encoded proteinases and proteolytic processing in the Nidovirales. J Gen Virol 81, 853–879.
Received 21 November 2007;
accepted 30 January 2008.
This article has been cited by other articles:
![]() |
P. C. Y. Woo, S. K. P. Lau, Y. Huang, and K.-Y. Yuen Coronavirus Diversity, Phylogeny and Interspecies Jumping Exp Biol Med, October 1, 2009; 234(10): 1117 - 1127. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. C. Y. Woo, S. K. P. Lau, C. S. F. Lam, K. K. Y. Lai, Y. Huang, P. Lee, G. S. M. Luk, K. C. Dyrting, K.-H. Chan, and K.-Y. Yuen Comparative Analysis of Complete Genome Sequences of Three Avian Coronaviruses Reveals a Novel Group 3c Coronavirus J. Virol., January 15, 2009; 83(2): 908 - 917. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. K. W. Chu, L. L. M. Poon, Y. Guan, and J. S. M. Peiris Novel Astroviruses in Insectivorous Bats J. Virol., September 15, 2008; 82(18): 9107 - 9114. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |