Evaluation of Affymetrix Severe Acute Respiratory Syndrome Resequencing GeneChips in Characterization of the Genomes of Two Strains of Coronavirus Infecting Humans

Severe acute respiratory syndrome (SARS) was discovered during a recent global outbreak of atypical pneumonia. A number of immunologic and molecular studies of the clinical samples led to the conclusion that a novel coronavirus (SARS-CoV) was associated with the outbreak. Later, a SARS resequencing GeneChip was developed by Affymetrix to characterize the complete genome of SARS-CoV on a single GeneChip. The present study was carried out to evaluate the performance of SARS resequencing GeneChips. Two human SARS-CoV strains (CDC#200301157 and Urbani) were resequenced by the SARS GeneChips. Five overlapping PCR amplicons were generated for each strain and hybridized with these GeneChips. The successfully hybridized GeneChips generated nucleotide sequences of nearly complete genomes for the two SARS-CoV strains with an average call rate of 94.6%. Multiple alignments of nucleotide sequences obtained from SARS GeneChips and conventional sequencing revealed full concordance. Furthermore, the GeneChip-based analysis revealed no additional polymorphic sites. The results of this study suggest that GeneChip-based genome characterization is fast and reproducible. Thus, SARS resequencing GeneChips may be employed as an alternate tool to obtain genome sequences of SARS-CoV strains pathogenic for humans in order to further understand the transmission dynamics of these viruses.     

Moderate mutation rate in the SARS coronavirus genome and its implications

Background

The outbreak of severe acute respiratory syndrome (SARS) caused a severe global epidemic in 2003 which led to hundreds of deaths and many thousands of hospitalizations. The virus causing SARS was identified as a novel coronavirus (SARS-CoV) and multiple genomic sequences have been revealed since mid-April, 2003. After a quiet summer and fall in 2003, the newly emerged SARS cases in Asia, particularly the latest cases in China, are reinforcing a wide-spread belief that the SARS epidemic would strike back. With the understanding that SARS-CoV might be with humans for years to come, knowledge of the evolutionary mechanism of the SARS-CoV, including its mutation rate and emergence time, is fundamental to battle this deadly pathogen. To date, the speed at which the deadly virus evolved in nature and the elapsed time before it was transmitted to humans remains poorly understood.

Results

Sixteen complete genomic sequences with available clinical histories during the SARS outbreak were analyzed. After careful examination of multiple-sequence alignment, 114 single nucleotide variations were identified. To minimize the effects of sequencing errors and additional mutations during the cell culture, three strategies were applied to estimate the mutation rate by 1) using the closely related sequences as background controls; 2) adjusting the divergence time for cell culture; or 3) using the common variants only. The mutation rate in the SARS-CoV genome was estimated to be 0.80 – 2.38 × 10^-3 nucleotide substitution per site per year which is in the same order of magnitude as other RNA viruses. The non-synonymous and synonymous substitution rates were estimated to be 1.16 – 3.30 × 10^-3 and 1.67 – 4.67 × 10^-3 per site per year, respectively. The most recent common ancestor of the 16 sequences was inferred to be present as early as the spring of 2002.

Conclusions

The estimated mutation rates in the SARS-CoV using multiple strategies were not unusual among coronaviruses and moderate compared to those in other RNA viruses. All estimates of mutation rates led to the inference that the SARS-CoV could have been with humans in the spring of 2002 without causing a severe epidemic.

     

Recognition and analysis of protein-coding genes in severe acute respiratory syndrome associated coronavirus

MOTIVATION: The recent outbreak of severe acute respiratory syndrome (SARS) caused by SARS coronavirus (SARS-CoV) has necessitated an in-depth molecular understanding of the virus to identify new drug targets. The availability of complete genome sequence of several strains of SARS virus provides the possibility of identification of protein-coding genes and defining their functions. Computational approach to identify protein-coding genes and their putative functions will help in designing experimental protocols. RESULTS: In this paper, a novel analysis of SARS genome using gene prediction method GeneDecipher developed in our laboratory has been presented. Each of the 18 newly sequenced SARS-CoV genomes has been analyzed using GeneDecipher. In addition to polyprotein 1ab(1), polyprotein 1a and the four genes coding for major structural proteins spike (S), small envelope (E), membrane (M) and nucleocapsid (N), six to eight additional proteins have been predicted depending upon the strain analyzed. Their lengths range between 61 and 274 amino acids. Our method also suggests that polyprotein 1ab, polyprotein 1a, S, M and N are proteins of viral origin and others are of prokaryotic. Putative functions of all predicted protein-coding genes have been suggested using conserved peptides present in their open reading frames. AVAILABILITY: Detailed results of GeneDecipher analysis of all the 18 strains of SARS-CoV genomes are available at http://www.igib.res.in/sarsanalysis.html     

SARS—CoV N蛋白与病毒mRNA相互作用的初步研究

SARS冠状病毒(SARS-CoV)是一种新型的冠状病毒,其基因组大小约为30,000 nt,为单股正链RNA。病毒基因中的1-72个核甘酸为前导序列。核衣壳(Nucleocapsid,N)蛋白是冠状病毒的主要结构蛋白,它在病毒基因转录,翻译以及病毒颗粒包装中起重要作用。在本研究中,我们通过PCR的方法从SARS-CoV cDNA中克隆N基因,将基因克隆到大肠杆菌表达载体中,经表达纯化获得大量重组蛋白,通过亲和层析和凝胶过滤层析获得高纯度的N蛋白。同时构建前导RNA的转录模板,经体外转录得到地高辛标记的RNA。使用Northwestern分析技术,我们证实纯化的N蛋白在体外可以与RNA发生特异性的结合。N蛋白与病毒RNA的结合特性及其在病毒生活周期中的所起作用的初步研究,为下一步设计出有效的阻断病毒周期从而达到抗病毒目的的药物或疫苗奠定了基础。     

SARS相关冠状病毒及其基因组

正在全球部分地区流行的严重急性呼吸综合征(SARS),由于其传染性强、危害性大而引起了广泛关注。各国实验室密切协作,在数月时间内分离出了SARS冠状病毒(SARSCoV),测定了病毒基因组序列,并在猴体内初步再现出SARSCoV所致肺部疾病与人相似,这些工作为遏制SARS的蔓延发挥了重要作用。现就SARSCoV的鉴定、基因组及其产物的结构与功能做一综述。     

Viral discovery and sequence recovery using DNA microarrays

Because of the constant threat posed by emerging infectious diseases and the limitations of existing approaches used to identify new pathogens, there is a great demand for new technological methods for viral discovery. We describe herein a DNA microarray-based platform for novel virus identification and characterization. Central to this approach was a DNA microarray designed to detect a wide range of known viruses as well as novel members of existing viral families; this microarray contained the most highly conserved 70mer sequences from every fully sequenced reference viral genome in GenBank. During an outbreak of severe acute respiratory syndrome (SARS) in March 2003, hybridization to this microarray revealed the presence of a previously uncharacterized coronavirus in a viral isolate cultivated from a SARS patient. To further characterize this new virus, approximately 1 kb of the unknown virus genome was cloned by physically recovering viral sequences hybridized to individual array elements. Sequencing of these fragments confirmed that the virus was indeed a new member of the coronavirus family. This combination of array hybridization followed by direct viral sequence recovery should prove to be a general strategy for the rapid identification and characterization of novel viruses and emerging infectious disease.     

用SARS冠状病毒全基因组芯片杂交方法分析SARS-CoV

为从临床样品中检测和分析SARSCoV病毒打基础,并为分析SARSCoV病毒的复制和转录等机理提供一种有效方法。以SARS冠状病毒TOR2株序列作为标准设计和制备一种覆盖SARS冠状病毒全基因组的寡聚核苷酸芯片,探针长度为70nt,每相邻的探针序列重复25nt,共660条。用该芯片分析了细胞培养的SARSCoV病毒总RNA、7个SARSCoV病毒的基因克隆片段。对RNA样品用随机引物进行反转录PCR获得cDNA。对DNA用随机引物扩增和dUTPcy3标记。结果用这种芯片杂交检测SARSCoV病毒RNA可见阳性信号呈全基因组分布,并且有多处连续的阳性信号点;用正常人的白细胞RNA为对照,杂交未出现明显阳性信号。检测7个SARSCoV病毒基因克隆片段,在该片段相应的探针区段出现连续阳性信号点。这种方法可有效地检测和分析样品中SARS冠状病毒全基因组的信息。     

SARS病毒与其他冠状病毒的基因组比较

本文利用生物信息学方法比较SARS病毒和其他冠状病毒基因组。通过数据库搜索,找出与SARS病毒基因组相似的核酸或蛋白质序列,并对相似序列进行比对,分析它们的共性和差异。结果表明,SARS病毒在基因组的组织上及结构蛋白质方面与现有冠状病毒有比较大的相似性,SARS病毒基因组与冠状病毒基因组相关。但是,SARS病毒基因组还存在一些特异性序列,ORF1a和S蛋白(特别是S1)的变化以及SARS—CoV特异性的非结构蛋白可能是SARS发病机理与传染特性区别于其他冠状病毒的分子基础。在全基因组水平上进行核酸单词出现频率分析,结果表明,SARS病毒远离已知的其他冠状病毒,单独成为一类。     

Evolution and variation of the SARS-CoV genome

Knowledge of the evolution of pathogens is of great medical and biological significance to the prevention, diagnosis, and therapy of infectious diseases. In order to understand the origin and evolution of the SARS-CoV (severe acute respiratory syndrome-associated coronavirus), we collected complete genome sequences of all viruses available in GenBank, and made comparative analyses with the SARS-CoV. Genomic signature analysis demonstrates that the coronaviruses all take the TGTT as their richest tetranucleotide except the SARS-CoV. A detailed analysis of the forty-two complete SARS-CoV genome sequences revealed the existence of two distinct genotypes, and showed that these isolates could be classified into four groups. Our manual analysis of the BLASTN results demonstrates that the HE (hemagglutinin-esterase) gene exists in the SARS-CoV, and many mutations made it unfamiliar to us.     

Receptor and viral determinants of SARS-coronavirus adaptation to human ACE2

Human angiotensin-converting enzyme 2 (ACE2) is a functional receptor for SARS coronavirus (SARS-CoV). Here we identify the SARS-CoV spike (S)-protein-binding site on ACE2. We also compare S proteins of SARS-CoV isolated during the 2002-2003 SARS outbreak and during the much less severe 2003-2004 outbreak, and from palm civets, a possible source of SARS-CoV found in humans. All three S proteins bound to and utilized palm-civet ACE2 efficiently, but the latter two S proteins utilized human ACE2 markedly less efficiently than did the S protein obtained during the earlier human outbreak. The lower affinity of these S proteins could be complemented by altering specific residues within the S-protein-binding site of human ACE2 to those of civet ACE2, or by altering S-protein residues 479 and 487 to residues conserved during the 2002-2003 outbreak. Collectively, these data describe molecular interactions important to the adaptation of SARS-CoV to human cells, and provide insight into the severity of the 2002-2003 SARS epidemic.     

Molecular epidemiology of the coronavirus associated with severe acute respiratory syndrome: a review of data from the Chinese University of Hong Kong

The epidemic of the severe acute respiratory syndrome (SARS) has swept through the globe with more than 8000 reported probable cases. In Hong Kong, the hardest hit areas included our teaching hospital and the Amoy Gardens apartment complex. A novel coronavirus, SARS-coronavirus (SARS-CoV), with a single-stranded plus sense RNA genome, was promptly implicated as the causative agent and subsequently fulfilled Koch's postulates. To aid the understanding of SARS-CoV, groups of investigators rapidly sequenced viral isolates around the world. We were the third group in the world to release the complete SARS-CoV genome sequence (isolate CUHK-W1) on the world-wide web. With other isolates from patients of distinct epidemiological backgrounds, we additionally sequenced four complete (CUHK-Su10, CUHK-AG01, CUHK-AG02, CUHK-AG03) and two partial SARS-CoV genomes. The reviewed data obtained from representative patients from the hospital and community outbreaks has documented the evolution of the virus in this epidemic. Their sequence variations also revealed a remarkable epidemiological correlation. We demonstrate that sequence variations in the SARS-CoV genome can be applied as a powerful molecular tool in tracing the route of transmission, when used adjunctively with standard epidemiology.     

Bioinformatics analysis of SARS coronavirus genome polymorphism

Background  

We have compared 38 isolates of the SARS-CoV complete genome. The main goal was twofold: first, to analyze and compare nucleotide sequences and to identify positions of single nucleotide polymorphism (SNP), insertions and deletions, and second, to group them according to sequence similarity, eventually pointing to phylogeny of SARS-CoV isolates. The comparison is based on genome polymorphism such as insertions or deletions and the number and positions of SNPs.     

Molecular evolution and multilocus sequence typing of 145 strains of SARS-CoV

In this study, we have identified 876 polymorphism sites in 145 complete or partial genomes of SARS-CoV available in the NCBI GenBank. One hundred and seventy-four of these sites existed in two or more SARS-CoV genome sequences. According to the sequence polymorphism, all SARS-CoVs can be divided into three groups: (I) group 1, animal-origin viruses (such as SARS-CoV SZ1, SZ3, SZ13 and SZ16); (II) group 2, all viruses with clinical origin during first epidemic; and (III) group 3, SARS-CoV GD03T0013. According to 10 special loci, group 2 again can be divided into genotypes C and T, which can be further divided into sub-genotypes C1-C4 and T1-T4. Positive Darwinian selections were identified between any pair of these three groups. Genotype C gives neutral selection. Genotype T, however, shows negative selection. By comparing the death rates of SARS patients in the different regions, it was found that the death rate caused by the viruses of the genotype C was lower than that of the genotype T. SARS-CoVs might originate from an unknown ancestor.     

Lack of peripheral memory B cell responses in recovered patients with severe acute respiratory syndrome: a six-year follow-up study

Six years have passed since the outbreak of severe acute respiratory syndrome (SARS). Previous studies indicated that specific Abs to SARS-related coronavirus (SARS-CoV) waned over time in recovered SARS patients. It is critical to find out whether a potential anamnestic response, as seen with other viral infections, exists to protect a person from reinfection in case of another SARS outbreak. Recovered SARS patients were followed up to 6 y to estimate the longevity of specific Ab. The specific memory B cell and T cell responses to SARS-CoV Ags were measured by means of ELISPOT assay. Factors in relation to humoral and cellular immunity were investigated. Six years postinfection, specific IgG Ab to SARS-CoV became undetectable in 21 of the 23 former patients. No SARS-CoV Ag-specific memory B cell response was detected in either 23 former SARS patients or 22 close contacts of SARS patients. Memory T cell responses to a pool of SARS-CoV S peptides were identified in 14 of 23 (60.9%) recovered SARS patients, whereas there was no such specific response in either close contacts or healthy controls. Patients with more severe clinical manifestations seemed to present a higher level of Ag-specific memory T cell response. SARS-specific IgG Ab may eventually vanish and peripheral memory B cell responses are undetectable in recovered SARS patients. In contrast, specific T cell anamnestic responses can be maintained for at least 6 y. These findings have applications in preparation for the possible reemergence of SARS.     

SARS transmission pattern in Singapore reassessed by viral sequence variation analysis

Background

Epidemiological investigations of infectious disease are mainly dependent on indirect contact information and only occasionally assisted by characterization of pathogen sequence variation from clinical isolates. Direct sequence analysis of the pathogen, particularly at a population level, is generally thought to be too cumbersome, technically difficult, and expensive. We present here a novel application of mass spectrometry (MS)–based technology in characterizing viral sequence variations that overcomes these problems, and we apply it retrospectively to the severe acute respiratory syndrome (SARS) outbreak in Singapore.

Methods and Findings

The success rate of the MS-based analysis for detecting SARS coronavirus (SARS-CoV) sequence variations was determined to be 95% with 75 copies of viral RNA per reaction, which is sufficient to directly analyze both clinical and cultured samples. Analysis of 13 SARS-CoV isolates from the different stages of the Singapore outbreak identified nine sequence variations that could define the molecular relationship between them and pointed to a new, previously unidentified, primary route of introduction of SARS-CoV into the Singapore population. Our direct determination of viral sequence variation from a clinical sample also clarified an unresolved epidemiological link regarding the acquisition of SARS in a German patient. We were also able to detect heterogeneous viral sequences in primary lung tissues, suggesting a possible coevolution of quasispecies of virus within a single host.

Conclusion

This study has further demonstrated the importance of improving clinical and epidemiological studies of pathogen transmission through the use of genetic analysis and has revealed the MS-based analysis to be a sensitive and accurate method for characterizing SARS-CoV genetic variations in clinical samples. We suggest that this approach should be used routinely during outbreaks of a wide variety of agents, in order to allow the most effective control.     

Prediction of functional class of the SARS coronavirus proteins by a statistical learning method

The complete genome of severe acute respiratory syndrome coronavirus (SARS-CoV) reveals the existence of putative proteins unique to SARS-CoV. Identification of their function facilitates a mechanistic understanding of SARS infection and drug development for its treatment. The sequence of the majority of these putative proteins has no significant similarity to those of known proteins, which complicates the task of using sequence analysis tools to probe their function. Support vector machines (SVM), useful for predicting the functional class of distantly related proteins, is employed to ascribe a possible functional class to SARS-CoV proteins. Testing results indicate that SVM is able to predict the functional class of 73% of the known SARS-CoV proteins with available sequences and 67% of 18 other novel viral proteins. A combination of the sequence comparison method BLAST and SVMProt can further improve the prediction accuracy of SMVProt such that the functional class of two additional SARS-CoV proteins is correctly predicted. Our study suggests that the SARS-CoV genome possibly contains a putative voltage-gated ion channel, structural proteins, a carbon-oxygen lyase, oxidoreductases acting on the CH-OH group of donors, and an ATP-binding cassette transporter. A web version of our software, SVMProt, is accessible at http://jing.cz3.nus.edu.sg/cgi-bin/svmprot.cgi .