Prevalence and genetic diversity of coronaviruses in bats from China

The Widespread Evolutionary Significance of Viruses

In the last 30 years, the study of virus evolution has undergone a transformation. Originally concerned with disease and its emergence, virus evolution had not been well integrated into the general study of evolution. This chapter reviews the developments that have brought us to this new appreciation for the general significance of virus evolution to all life. We now know that viruses numerically dominate all habitats of life, especially the oceans. Theoretical developments in the 1970s regarding quasispecies, error rates, and error thresholds have yielded many practical insights into virus–host dynamics. The human diseases of HIV-1 and hepatitis C virus cannot be understood without this evolutionary framework. Yet recent developments with poliovirus demonstrate that viral fitness can be the result of a consortia, not one fittest type, a basic Darwinian concept in evolutionary biology. Darwinian principles do apply to viruses, such as with Fisher population genetics, but other features, such as reticulated and quasispecies-based evolution distinguish virus evolution from classical studies. The available phylogenetic tools have greatly aided our analysis of virus evolution, but these methods struggle to characterize the role of virus populations. Missing from many of these considerations has been the major role played by persisting viruses in stable virus evolution and disease emergence. In many cases, extreme stability is seen with persisting RNA viruses. Indeed, examples are known in which it is the persistently infected host that has better survival. We have also recently come to appreciate the vast diversity of phage (DNA viruses) of prokaryotes as a system that evolves by genetic exchanges across vast populations (Chapter 10). This has been proposed to be the “big bang” of biological evolution. In the large DNA viruses of aquatic microbes we see surprisingly large, complex and diverse viruses. With both prokaryotic and eukaryotic DNA viruses, recombination is the main engine of virus evolution, and virus host co-evolution is common, although not uniform. Viral emergence appears to be an unending phenomenon and we can currently witness a selective sweep by retroviruses that infect and become endogenized in koala bears.

Determination and application of immunodominant regions of SARS coronavirus spike and nucleocapsid proteins recognized by sera from different animal species

Knowledge of immunodominant regions in major viral antigens is important for rational design of effective vaccines and diagnostic tests. Although there have been many reports of such work done for SARS–CoV, these were mainly focused on the immune responses of humans and mice. In this study, we aim to search for and compare immunodominant regions of the spike (S) and nucleocapsid (N) proteins which are recognized by sera from different animal species, including mouse, rat, rabbit, civet, pig and horse. Twelve overlapping recombinant protein fragments were produced in Escherichia coli, six each for the S and N proteins, which covered the entire coding region of the two proteins. Using a membrane-strip based Western blot approach, the reactivity of each antigen fragment against a panel of animal sera was determined. Immunodominant regions containing linear epitopes, which reacted with sera from all the species tested, were identified for both proteins. The S3 fragment (aa 402–622) and the N4 fragment (aa 220–336) were the most immunodominant among the six S and N fragments, respectively. Antibodies raised against the S3 fragment were able to block the binding of a panel of S-specific monoclonal antibodies (mAb) to SARS–CoV in ELISA, further demonstrating the immunodominance of this region. Based on these findings, one-step competition ELISAs were established which were able to detect SARS–CoV antibodies from human and at least seven different animal species. Considering that a large number of animal species are known to be susceptible to SARS–CoV, these assays will be a useful tool to trace the origin and transmission of SARS–CoV and to minimise the risk of animal-to-human transmission.

Difference in receptor usage between severe acute respiratory syndrome (SARS) coronavirus and SARS-like coronavirus of bat origin

Severe acute respiratory syndrome (SARS) is caused by the SARS-associated coronavirus (SARS-CoV), which uses angiotensin-converting enzyme 2 (ACE2) as its receptor for cell entry. A group of SARS-like CoVs (SL-CoVs) has been identified in horseshoe bats. SL-CoVs and SARS-CoVs share identical genome organizations and high sequence identities, with the main exception of the N terminus of the spike protein (S), known to be responsible for receptor binding in CoVs. In this study, we investigated the receptor usage of the SL-CoV S by combining a human immunodeficiency virus-based pseudovirus system with cell lines expressing the ACE2 molecules of human, civet, or horseshoe bat. In addition to full-length S of SL-CoV and SARS-CoV, a series of S chimeras was constructed by inserting different sequences of the SARS-CoV S into the SL-CoV S backbone. Several important observations were made from this study. First, the SL-CoV S was unable to use any of the three ACE2 molecules as its receptor. Second, the SARS-CoV S failed to enter cells expressing the bat ACE2. Third, the chimeric S covering the previously defined receptor-binding domain gained its ability to enter cells via human ACE2, albeit with different efficiencies for different constructs. Fourth, a minimal insert region (amino acids 310 to 518) was found to be sufficient to convert the SL-CoV S from non-ACE2 binding to human ACE2 binding, indicating that the SL-CoV S is largely compatible with SARS-CoV S protein both in structure and in function. The significance of these findings in relation to virus origin, virus recombination, and host switching is discussed.

Detection of group 1 coronaviruses in bats using universal coronavirus reverse transcription polymerase chain reactions

The zoonotic transmission of SARS coronavirus from animals to humans revealed the potential impact of coronaviruses on mankind. This incident also triggered several surveillance programs to hunt for novel coronaviruses in human and wildlife populations. Using classical RT-PCR assays that target a highly conserved sequence among coronaviruses, we identified the first coronaviruses in bats. These assays and the cloning and sequencing of the PCR products are described in this chapter. Using the same approach in our subsequent studies, we further detected several novel coronaviruses in bats. These findings highlighted the fact that bats are important reservoirs for coronaviruses.

Infectious bronchitis viruses with a novel genomic organization

A number of novel infectious bronchitis viruses (IBVs) were previously identified in commercial poultry in Australia, where they caused significant economic losses. Since there has been only limited characterization of these viruses, we investigated the genomic and phenotypic differences between these novel IBVs and other, classical IBVs. The 3' 7.5 kb of the genomes of 17 Australian IBV strains were sequenced, and growth properties of 6 of the strains were compared. Comparison of sequences of the genes coding for structural and nonstructural proteins revealed the existence of two IBV genotypes: classical and novel. The genomic organization of the classical IBVs was typical of those of other group III coronaviruses: 5'-Pol-S-3a-3b-E-M-5a-5b-N-untranslated region (UTR)-3'. However, the novel IBV genotype lacked either all or most of the genes coding for nonstructural proteins at the 3' end of the genome and had a unique open reading frame, X1. The gene order was either 5'-Pol-S-X1-E-M-N-UTR-3' or 5'-Pol-S-X1-E-M-5b-N-UTR-3'. Phenotypically, novel and classical IBVs also differed; novel IBVs grew at a slower rate and reached lower titers in vitro and in vivo and were markedly less immunogenic in chicks. Although the novel IBVs induced histopathological lesions in the tracheas of infected chicks that were comparable to those induced by classical strains, they did not induce lesions in the kidneys. This study has demonstrated for the first time the existence of a naturally occurring IBV genotype devoid of some of the genes coding for nonstructural proteins and has also indicated that all of the accessory genes are dispensable for the growth of IBV and that such viruses are able to cause clinical disease and economic loss. The phylogenic differences between these novel IBVs and other avian coronaviruses suggest a reservoir host distinct from domestic poultry.

Severe acute respiratory syndrome coronavirus as an agent of emerging and reemerging infection

Before the emergence of severe acute respiratory syndrome (SARS) coronavirus (SARS-CoV) in 2003, only 12 other animal or human coronaviruses were known. The discovery of this virus was soon followed by the discovery of the civet and bat SARS-CoV and the human coronaviruses NL63 and HKU1. Surveillance of coronaviruses in many animal species has increased the number on the list of coronaviruses to at least 36. The explosive nature of the first SARS epidemic, the high mortality, its transient reemergence a year later, and economic disruptions led to a rush on research of the epidemiological, clinical, pathological, immunological, virological, and other basic scientific aspects of the virus and the disease. This research resulted in over 4,000 publications, only some of the most representative works of which could be reviewed in this article. The marked increase in the understanding of the virus and the disease within such a short time has allowed the development of diagnostic tests, animal models, antivirals, vaccines, and epidemiological and infection control measures, which could prove to be useful in randomized control trials if SARS should return. The findings that horseshoe bats are the natural reservoir for SARS-CoV-like virus and that civets are the amplification host highlight the importance of wildlife and biosecurity in farms and wet markets, which can serve as the source and amplification centers for emerging infections.

CoVDB: a comprehensive database for comparative analysis of coronavirus genes and genomes

The recent SARS epidemic has boosted interest in the discovery of novel human and animal coronaviruses. By July 2007, more than 3000 coronavirus sequence records, including 264 complete genomes, are available in GenBank. The number of coronavirus species with complete genomes available has increased from 9 in 2003 to 25 in 2007, of which six, including coronavirus HKU1, bat SARS coronavirus, group 1 bat coronavirus HKU2, groups 2c and 2d coronaviruses, were sequenced by our laboratory. To overcome the problems we encountered in the existing databases during comparative sequence analysis, we built a comprehensive database, CoVDB (http://covdb.microbiology.hku.hk), of annotated coronavirus genes and genomes. CoVDB provides a convenient platform for rapid and accurate batch sequence retrieval, the cornerstone and bottleneck for comparative gene or genome analysis. Sequences can be directly downloaded from the website in FASTA format. CoVDB also provides detailed annotation of all coronavirus sequences using a standardized nomenclature system, and overcomes the problems of duplicated and identical sequences in other databases. For complete genomes, a single representative sequence for each species is available for comparative analysis such as phylogenetic studies. With the annotated sequences in CoVDB, more specific blast search results can be generated for efficient downstream analysis.

Cytosine deamination and selection of CpG suppressed clones are the two major independent biological forces that shape codon usage bias in coronaviruses

Using the complete genome sequences of 19 coronavirus genomes, we analyzed the codon usage bias, dinucleotide relative abundance and cytosine deamination in coronavirus genomes. Of the eight codons that contain CpG, six were markedly suppressed. The mean NNU/NNC ratio of the six amino acids using either NNC or NNU as codon is 3.262, suggesting cytosine deamination. Among the 16 dinucleotides, CpG was most markedly suppressed (mean relative abundance 0.509). No correlation was observed between CpG abundance and mean NNU/NNC ratio. Among the 19 coronaviruses, CoV-HKU1 showed the most extreme codon usage bias and extremely high NNU/NNC ratio of 8.835. Cytosine deamination and selection of CpG suppressed clones by the immune system are the two major independent biochemical and biological selective forces that shape codon usage bias in coronavirus genomes. The underlying mechanism for the extreme codon usage bias, cytosine deamination and G + C content in CoV-HKU1 warrants further studies.

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

Apart from bat-SARS-CoV, we have identified a novel group 1 coronavirus, bat-CoV HKU2, in Rhinolophus sinicus (Chinese horseshoe bats). Since it has been suggested that the receptor-binding motif (RBM) of SARS-CoV may have been acquired from a group 1 coronavirus, we conducted a surveillance study and identified bat-SARS-CoV and bat-CoV HKU2 in 8.7% and 7.5% respectively of R. sinicus in Hong Kong and Guangdong. Complete genome sequencing of four strains of bat-CoV HKU2 revealed the smallest coronavirus genome (27164 nucleotides) and a unique spike protein evolutionarily distinct from the rest of the genome. This spike protein, sharing similar deletions with other group 2 coronaviruses in its C-terminus, also contained a 15-amino acid peptide homologous to a corresponding peptide within the RBM of spike protein of SARS-CoV, which was absent in other coronaviruses except bat-SARS-CoV. These suggest a common evolutionary origin in the spike protein of bat-CoV HKU2, bat-SARS-CoV, and SARS-CoV.

Review of bats and SARS

TOC Summary: The discovery of SARS-like coronaviruses in horseshoe bats highlights the possibility of future outbreaks caused by different coronaviruses of bat origin.

Bats have been identified as a natural reservoir for an increasing number of emerging zoonotic viruses, including henipaviruses and variants of rabies viruses. Recently, we and another group independently identified several horseshoe bat species (genus Rhinolophus) as the reservoir host for a large number of viruses that have a close genetic relationship with the coronavirus associated with severe acute respiratory syndrome (SARS). Our current research focused on the identification of the reservoir species for the progenitor virus of the SARS coronaviruses responsible for outbreaks during 2002–2003 and 2003–2004. In addition to SARS-like coronaviruses, many other novel bat coronaviruses, which belong to groups 1 and 2 of the 3 existing coronavirus groups, have been detected by PCR. The discovery of bat SARS-like coronaviruses and the great genetic diversity of coronaviruses in bats have shed new light on the origin and transmission of SARS coronaviruses.

Detection of a novel and highly divergent coronavirus from asian leopard cats and Chinese ferret badgers in Southern China

Since an outbreak of severe acute respiratory syndrome (SARS) was averted in 2004, many novel coronaviruses have been recognized from different species, including humans. Bats have provided the most diverse assemblages of coronaviruses, suggesting that they may be the natural reservoir. Continued virological surveillance has proven to be the best way to avert this infectious disease at the source. Here we provide the first description of a previously unidentified coronavirus lineage detected from wild Asian leopard cats (Prionailurus bengalensis) and Chinese ferret badgers (Melogale moschata) during virological surveillance in southern China. Partial genome analysis revealed a typical coronavirus genome but with a unique putative accessory gene organization. Phylogenetic analyses revealed that the envelope, membrane, and nucleoprotein structural proteins and the two conserved replicase domains, putative RNA-dependent RNA polymerase and RNA helicase, of these novel coronaviruses were most closely related to those of group 3 coronaviruses identified from birds, while the spike protein gene was most closely related to that of group 1 coronaviruses from mammals. However, these viruses always fell into an outgroup phylogenetic relationship with respect to other coronaviruses and had low amino acid similarity to all known coronavirus groups, indicating that they diverged early in the evolutionary history of coronaviruses. These results suggest that these viruses may represent a previously unrecognized evolutionary pathway, or possibly an unidentified coronavirus group. This study demonstrates the importance of systematic virological surveillance in market animals for understanding the evolution and emergence of viruses with infectious potential.

Emerging viruses: coming in on a wrinkled wing and a prayer

The role that bats have played in the emergence of several new infectious diseases has been under review. Bats have been identified as the reservoir hosts of newly emergent viruses such as Nipah virus, Hendra virus, and severe acute respiratory syndrome-like coronaviruses. This article expands on recent findings about bats and viruses and their relevance to human infections. It briefly reviews the history of chiropteran viruses and discusses their emergence in the context of geography, phylogeny, and ecology. The public health and trade impacts of several outbreaks are also discussed. Finally, we attempt to predict where, when, and why we may see the emergence of new chiropteran viruses.

Evolutionary insights into the ecology of coronaviruses

Although many novel members of the Coronaviridae have recently been recognized in different species, the ecology of coronaviruses has not been established. Our study indicates that bats harbor a much wider diversity of coronaviruses than any other animal species. Dating of different coronavirus lineages suggests that bat coronaviruses are older than those recognized in other animals and that the human severe acute respiratory syndrome (SARS) coronavirus was directly derived from viruses from wild animals in wet markets of southern China. Furthermore, the most closely related bat and SARS coronaviruses diverged in 1986, an estimated divergence time of 17 years prior to the outbreak, suggesting that there may have been transmission via an unknown intermediate host. Analysis of lineage-specific selection pressure also indicated that only SARS coronaviruses in civets and humans were under significant positive selection, also demonstrating a recent interspecies transmission. Analysis of population dynamics revealed that coronavirus populations in bats have constant population growth, while viruses from all other hosts show epidemic-like increases in population. These results indicate that diverse coronaviruses are endemic in different bat species, with repeated introductions to other animals and occasional establishment in other species. Our findings suggest that bats are likely the natural hosts for all presently known coronavirus lineages and that all coronaviruses recognized in other species were derived from viruses residing in bats. Further surveillance of bat and other animal populations is needed to fully describe the ecology and evolution of this virus family.

Severe Acute Respiratory Syndrome Coronavirus Entry as a Target of Antiviral Therapies

The identification in 2003 of a coronavirus as the aetiological agent of severe acute respiratory syndrome (SARS) intensified efforts to understand the biology of corona-viruses in general and SARS coronavirus (SARS-CoV) in particular. Rapid progress was made in describing the SARS-CoV genome, evolution and lifecycle. Identification of angiotensin-converting enzyme 2 (ACE2) as an obligate cellular receptor for SARS-CoV contributed to understanding of the SARS-CoV entry process, and helped to characterize two targets of antiviral therapeutics: the SARS-CoV spike protein and ACE2. Here we describe the role of these proteins in SARS-CoV replication and potential therapeutic strategies aimed at preventing entry of SARS-CoV into target cells.

Human coronavirus 229E encodes a single ORF4 protein between the spike and the envelope genes

Background

The genome of coronaviruses contains structural and non-structural genes, including several so-called accessory genes. All group 1b coronaviruses encode a single accessory protein between the spike and envelope genes, except for human coronavirus (HCoV) 229E. The prototype virus has a split gene, encoding the putative ORF4a and ORF4b proteins. To determine whether primary HCoV-229E isolates exhibit this unusual genome organization, we analyzed the ORF4a/b region of five current clinical isolates from The Netherlands and three early isolates collected at the Common Cold Unit (CCU) in Salisbury, UK.

Results

All Dutch isolates were identical in the ORF4a/b region at amino acid level. All CCU isolates are only 98% identical to the Dutch isolates at the nucleotide level, but more closely related to the prototype HCoV-229E (>98%). Remarkably, our analyses revealed that the laboratory adapted, prototype HCoV-229E has a 2-nucleotide deletion in the ORF4a/b region, whereas all clinical isolates carry a single ORF, 660 nt in size, encoding a single protein of 219 amino acids, which is a homologue of the ORF3 proteins encoded by HCoV-NL63 and PEDV.

Conclusion

Thus, the genome organization of the group 1b coronaviruses HCoV-NL63, PEDV and HCoV-229E is identical. It is possible that extensive culturing of the HCoV-229E laboratory strain resulted in truncation of ORF4. This may indicate that the protein is not essential in cell culture, but the highly conserved amino acid sequence of the ORF4 protein among clinical isolates suggests that the protein plays an important role in vivo.

The cytoplasmic tail of the severe acute respiratory syndrome coronavirus spike protein contains a novel endoplasmic reticulum retrieval signal that binds COPI and promotes interaction with membrane protein

Like other coronaviruses, severe acute respiratory syndrome coronavirus (SARS CoV) assembles at and buds into the lumen of the endoplasmic reticulum (ER)-Golgi intermediate compartment (ERGIC). Accumulation of the viral envelope proteins at this compartment is a prerequisite for virus assembly. Previously, we reported the identification of a dibasic motif (KxHxx) in the cytoplasmic tail of the SARS CoV spike (S) protein that was similar to a canonical dilysine ER retrieval signal. Here we demonstrate that this motif is a novel and functional ER retrieval signal which reduced the rate of traffic of the full-length S protein through the Golgi complex. The KxHxx motif also partially retained two different reporter proteins in the ERGIC region and reduced their rates of trafficking, although the motif was less potent than the canonical dilysine signal. The dibasic motif bound the coatomer complex I (COPI) in an in vitro binding assay, suggesting that ER retrieval may contribute to the accumulation of SARS CoV S protein near the virus assembly site for interaction with other viral structural proteins. In support of this, we found that the dibasic motif on the SARS S protein was required for its localization to the ERGIC/Golgi region when coexpressed with SARS membrane (M) protein. Thus, the cycling of SARS S through the ER-Golgi system may be required for its incorporation into assembling virions in the ERGIC.

Comparative analysis of twelve genomes of three novel group 2c and group 2d coronaviruses reveals unique group and subgroup features

Twelve complete genomes of three novel coronaviruses-bat coronavirus HKU4 (bat-CoV HKU4), bat-CoV HKU5 (putative group 2c), and bat-CoV HKU9 (putative group 2d)-were sequenced. Comparative genome analysis showed that the various open reading frames (ORFs) of the genomes of the three coronaviruses had significantly higher amino acid identities to those of other group 2 coronaviruses than group 1 and 3 coronaviruses. Phylogenetic trees constructed using chymotrypsin-like protease, RNA-dependent RNA polymerase, helicase, spike, and nucleocapsid all showed that the group 2a and 2b and putative group 2c and 2d coronaviruses are more closely related to each other than to group 1 and 3 coronaviruses. Unique genomic features distinguishing between these four subgroups, including the number of papain-like proteases, the presence or absence of hemagglutinin esterase, small ORFs between the membrane and nucleocapsid genes and ORFs (NS7a and NS7b), bulged stem-loop and pseudoknot structures downstream of the nucleocapsid gene, transcription regulatory sequence, and ribosomal recognition signal for the envelope gene, were also observed. This is the first time that NS7a and NS7b downstream of the nucleocapsid gene has been found in a group 2 coronavirus. The high Ka/Ks ratio of NS7a and NS7b in bat-CoV HKU9 implies that these two group 2d-specific genes are under high selective pressure and hence are rapidly evolving. The four subgroups of group 2 coronaviruses probably originated from a common ancestor. Further molecular epidemiological studies on coronaviruses in the bats of other countries, as well as in other animals, and complete genome sequencing will shed more light on coronavirus diversity and their evolutionary histories.

Characterization of White bream virus reveals a novel genetic cluster of nidoviruses

The order Nidovirales comprises viruses from the families Coronaviridae (genera Coronavirus and Torovirus), Roniviridae (genus Okavirus), and Arteriviridae (genus Arterivirus). In this study, we characterized White bream virus (WBV), a bacilliform plus-strand RNA virus isolated from fish. Analysis of the nucleotide sequence, organization, and expression of the 26.6-kb genome provided conclusive evidence for a phylogenetic relationship between WBV and nidoviruses. The polycistronic genome of WBV contains five open reading frames (ORFs), called ORF1a, -1b, -2, -3, and -4. In WBV-infected cells, three subgenomic RNAs expressing the structural proteins S, M, and N were identified. The subgenomic RNAs were revealed to share a 42-nucleotide, 5' leader sequence that is identical to the 5'-terminal genome sequence. The data suggest that a conserved nonanucleotide sequence, CA(G/A)CACUAC, located downstream of the leader and upstream of the structural protein genes acts as the core transcription-regulating sequence element in WBV. Like other nidoviruses with large genomes (>26 kb), WBV encodes in its ORF1b an extensive set of enzymes, including putative polymerase, helicase, ribose methyltransferase, exoribonuclease, and endoribonuclease activities. ORF1a encodes several membrane domains, a putative ADP-ribose 1"-phosphatase, and a chymotrypsin-like serine protease whose activity was established in this study. Comparative sequence analysis revealed that WBV represents a separate cluster of nidoviruses that significantly diverged from toroviruses and, even more, from coronaviruses, roniviruses, and arteriviruses. The study adds to the amazing diversity of nidoviruses and appeals for a more extensive characterization of nonmammalian nidoviruses to better understand the evolution of these largest known RNA viruses.