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.     

Immunopathogenesis of coronavirus infections: implications for SARS

At the end of 2002, the first cases of severe acute respiratory syndrome (SARS) were reported, and in the following year, SARS resulted in considerable mortality and morbidity worldwide. SARS is caused by a novel species of coronavirus (SARS-CoV) and is the most severe coronavirus-mediated human disease that has been described so far. On the basis of similarities with other coronavirus infections, SARS might, in part, be immune mediated. As discussed in this Review, studies of animals that are infected with other coronaviruses indicate that excessive and sometimes dysregulated responses by macrophages and other pro-inflammatory cells might be particularly important in the pathogenesis of disease that is caused by infection with these viruses. It is hoped that lessons from such studies will help us to understand more about the pathogenesis of SARS in humans and to prevent or control outbreaks of SARS in the future.     

Evidence of the recombinant origin of a bat severe acute respiratory syndrome (SARS)-like coronavirus and its implications on the direct ancestor of SARS coronavirus

Bats have been identified as the natural reservoir of severe acute respiratory syndrome (SARS)-like and SARS coronaviruses (SLCoV and SCoV). However, previous studies suggested that none of the currently sampled bat SLCoVs is the descendant of the direct ancestor of SCoV, based on their relatively distant phylogenetic relationship. In this study, evidence of the recombinant origin of the genome of a bat SLCoV is demonstrated. We identified a potential recombination breakpoint immediately after the consensus intergenic sequence between open reading frame 1 and the S coding region, suggesting the replication intermediates may participate in the recombination event, as previously speculated for other CoVs. Phylogenetic analysis of its parental regions suggests the presence of an uncharacterized SLCoV lineage that is phylogenetically closer to SCoVs than any of the currently sampled bat SLCoVs. Using various Bayesian molecular-clock models, interspecies transfer of this SLCoV lineage from bats to the amplifying host (e.g., civets) was estimated to have happened a median of 4.08 years before the SARS outbreak. Based on this relatively short window period, we speculate that this uncharacterized SLCoV lineage may contain the direct ancestor of SCoV. This study sheds light on the possible host bat species of the direct ancestor of SCoV, providing valuable information on the scope and focus of surveillance for the origin of SCoV.     

Heterozygosity and mutation rate: evidence for an interaction and its implications

If natural selection chose where new mutations occur it might well favour placing them near existing polymorphisms, thereby avoiding disruption of areas that work while adding novelty to regions where variation is tolerated or even beneficial. Such a system could operate if heterozygous sites are recognised and ‘repaired’ during the initial stages of crossing over. Such repairs involve an extra round of DNA replication, providing an opportunity for further mutations, thereby raising the local mutation rate. If so, the changes in heterozygosity that occur when populations grow or shrink could feed back to modulate both the rate and the distribution of mutations. Here, I review evidence from isozymes, microsatellites and single nucleotide polymorphisms that this potential is realised in real populations. I then consider the likely implications, focusing particularly on how these processes might affect microsatellites, concluding that heterozygosity does impact on the rate and distribution of mutations.     

人类基因组中突变对紧邻碱基的依赖性与重组率

减数分裂重组通过基因转变、碱基替换等方式影响基因组进化。紧邻碱基对突变偏好性有很强的影响,但该“紧邻碱基效应”如何随重组率变化有待深入研究。本文提出基于条件互信息(Conditional mutual information)量化突变对紧邻碱基依赖性的方法,并利用人类SNP等相关数据,分析重组率如何影响突变对紧邻碱基的依赖性。结果表明:在全基因组水平上, SNP位点上的突变对紧邻碱基的依赖性(即平均条件互信息)随着重组率的增加而增加;具体而言,当SNPs两侧碱基为A/G、C/G或C/T时,随着重组率的增加突变偏向性增强,但两侧碱基为A/A或T/T时,重组率对SNP突变偏向性产生抑制作用;另外,重组率越高,外显子与基因间区SNP的突变偏好性越强;而内含子区域SNP的突变偏好受到高重组率的抑制。结果有助于深入理解减数分裂重组如何影响基因组进化。     

Variation in recombination rate across the genome: evidence and implications

Recent data from humans and other species provide convincing evidence of variation in recombination rate in different genomic regions. Comparison of physical and genetic maps reveals variation on a scale of megabases, with substantial differences between sexes. Recombination is often suppressed near centromeres and elevated near telomeres, but neither of these observations is true for all chromosomes. In humans, patterns of linkage disequilibrium and experimental measures of recombination from sperm-typing reveal dramatic hotspots of recombination on a scale of kilobases. Genome-wide variation in the amount of crossing-over may be due to variation in the density of hotspots, the intensity of hotspots, or both. Theoretical models of selection and linkage predict that genetic variation will be reduced in regions of low recombination, and this prediction is supported by data from several species. Heterogeneity in rates of crossing-over provides both an opportunity and a challenge for identifying disease genes: as associations occur in blocks, genomic regions containing disease loci may be identified with relatively few markers, yet identifying the causal mutations is unlikely to be achieved through associations alone.     

Viral evolution and the emergence of SARS coronavirus

The recent appearance of severe acute respiratory syndrome coronavirus (SARS-CoV) highlights the continual threat to human health posed by emerging viruses. However, the central processes in the evolution of emerging viruses are unclear, particularly the selection pressures faced by viruses in new host species. We outline some of the key evolutionary genetic aspects of viral emergence. We emphasize that, although the high mutation rates of RNA viruses provide them with great adaptability and explain why they are the main cause of emerging diseases, their limited genome size means that they are also subject to major evolutionary constraints. Understanding the mechanistic basis of these constraints, particularly the roles played by epistasis and pleiotropy, is likely to be central in explaining why some RNA viruses are more able than others to cross species boundaries. Viral genetic factors have also been implicated in the emergence of SARS-CoV, with the suggestion that this virus is a recombinant between mammalian and avian coronaviruses. We show, however, that the phylogenetic patterns cited as evidence for recombination are more probably caused by a variation in substitution rate among lineages and that recombination is unlikely to explain the appearance of SARS in humans.     

Conformational reorganization of the SARS coronavirus spike following receptor binding: implications for membrane fusion

The SARS coronavirus (SARS-CoV) spike is the largest known viral spike molecule, and shares a similar function with all class 1 viral fusion proteins. Previous structural studies of membrane fusion proteins have largely used crystallography of static molecular fragments, in isolation of their transmembrane domains. In this study we have produced purified, irradiated SARS-CoV virions that retain their morphology, and are fusogenic in cell culture. We used cryo-electron microscopy and image processing to investigate conformational changes that occur in the entire spike of intact virions when they bind to the viral receptor, angiotensin-converting enzyme 2 (ACE2). We have shown that ACE2 binding results in structural changes that appear to be the initial step in viral membrane fusion, and precisely localized the receptor-binding and fusion core domains within the entire spike. Furthermore, our results show that receptor binding and subsequent membrane fusion are distinct steps, and that each spike can bind up to three ACE2 molecules. The SARS-CoV spike provides an ideal model system to study receptor binding and membrane fusion in the native state, employing cryo-electron microscopy and single-particle image analysis.     

Architecture of the SARS coronavirus prefusion spike

The emergence in 2003 of a new coronavirus (CoV) responsible for the atypical pneumonia termed severe acute respiratory syndrome (SARS) was a stark reminder that hitherto unknown viruses have the potential to cross species barriers to become new human pathogens. Here we describe the SARS-CoV 'spike' structure determined by single-particle cryo-EM, along with the docked atomic structures of the receptor-binding domain and prefusion core.     

A point mutation within the replicase gene differentially affects coronavirus genome versus minigenome replication

During the construction of the transmissible gastroenteritis virus (TGEV) full-length cDNA clone, a point mutation at position 637 that was present in the defective minigenome DI-C was maintained as a genetic marker. Sequence analysis of the recovered viruses showed a reversion at this position to the original virus sequence. The effect of point mutations at nucleotide 637 was analyzed by reverse genetics using a TGEV full-length cDNA clone and cDNAs from TGEV-derived minigenomes. The replacement of nucleotide 637 of TGEV genome by a T, as in the DI-C sequence, or an A severely affected virus recovery from the cDNA, yielding mutant viruses with low titers and small plaques compared to those of the wild type. In contrast, T or A at position 637 was required for minigenome rescue in trans by the helper virus. No relationship between these observations and RNA secondary-structure predictions was found, indicating that mutations at nucleotide 637 most likely had an effect at the protein level. Nucleotide 637 occupies the second codon position at amino acid 108 of the pp1a polyprotein. This position is predicted to map in the N-terminal polyprotein papain-like proteinase (PLP-1) cleavage site at the p9/p87 junction. Replacement of G-637 by A, which causes a drastic amino acid change (Gly to Asp) at position 108, affected PLP-1-mediated cleavage in vitro. A correlation was found between predicted cleaving and noncleaving mutations and efficient virus rescue from cDNA and minigenome amplification, respectively.     

Identifying inhibitors of the SARS coronavirus proteinase

The Severe Acute Respiratory Syndrome (SARS) is a serious respiratory illness that has recently been reported in parts of Asia and Canada. In this study, we use molecular dynamics (MD) simulations and docking techniques to screen 29 approved and experimental drugs against the theoretical model of the SARS CoV proteinase as well as the experimental structure of the transmissible gastroenteritis virus (TGEV) proteinase. Our predictions indicate that existing HIV-1 protease inhibitors, L-700,417 for instance, have high binding affinities and may provide good starting points for designing SARS CoV proteinase inhibitors.     

RNA recognition and cleavage by the SARS coronavirus endoribonuclease

The emerging disease SARS is caused by a novel coronavirus that encodes several unusual RNA-processing enzymes, including non-structural protein 15 (Nsp15), a hexameric endoribonuclease that preferentially cleaves at uridine residues. How Nsp15 recognizes and cleaves RNA is not well understood and is the subject of this study. Based on the analysis of RNA products separated by denaturing gel electrophoresis, Nsp15 has been reported to cleave both 5' and 3' of the uridine. We used several RNAs, including some with nucleotide analogs, and mass spectrometry to determine that Nsp15 cleaves only 3' of the recognition uridylate, with some cleavage 3' of cytidylate. A highly conserved RNA structure in the 3' non-translated region of the SARS virus was cleaved preferentially at one of the unpaired uridylate bases, demonstrating that both RNA structure and base-pairing can affect cleavage by Nsp15. Several modified RNAs that are not cleaved by Nsp15 can bind Nsp15 as competitive inhibitors. The RNA binding affinity of Nsp15 increased with the content of uridylate in substrate RNA and the co-factor Mn(2+). The hexameric form of Nsp15 was found to bind RNA in solution. A two-dimensional crystal of Nsp15 in complex with RNA showed that at least two RNA molecules could be bound per hexamer. Furthermore, an 8.3 A structure of Nsp15 was developed using cyroelectron microscopy, allowing us to generate a model of the Nsp15-RNA complex.     

Immunogenicity difference between the SARS coronavirus and the bat SARS-like coronavirus spike (S) proteins

SARS-like coronavirus (SL-CoV) in bats have a similar genomic organization to the human SARS-CoV. Their cognate gene products are highly conserved with the exception of the N-terminal region of the S proteins, which have only 63-64% sequence identity. The N-terminal region of coronavirus S protein is responsible for virus-receptor interaction. In this study, the immunogenicity of the SL-CoV S protein (S_SL) was studied and compared with that of SARS-CoV (S_SARS). DNA immunization in mice with S_SL elicited a high titer of antibodies against HIV-pseudotyped S_SL. The sera had low cross-reactivity, but no neutralization activity, for the HIV-pseudotyped S_SARS. Studies using wild bat sera revealed that it is highly likely that the immunodominant epitopes overlap with the major neutralizing sites of the SL-CoV S protein. These results demonstrated that SL-CoV and SARS-CoV shared only a limited number of immunogenic epitopes in their S proteins and the major neutralization epitopes are substantially different. This work provides useful information for future development of differential serologic diagnosis and vaccines for coronaviruses with different S protein sequences.     

mRNA cap-1 methyltransferase in the SARS genome

The 3D jury system has predicted the methyltransferase fold for the nsp13 protein of the SARS coronavirus. Based on the conservation of a characteristic tetrad of residues, the mRNA cap-1 methyltransferase function has been assigned to this protein, which has potential implications for antiviral therapy.     

A three-stemmed mRNA pseudoknot in the SARS coronavirus frameshift signal

A wide range of RNA viruses use programmed −1 ribosomal frameshifting for the production of viral fusion proteins. Inspection of the overlap regions between ORF1a and ORF1b of the SARS-CoV genome revealed that, similar to all coronaviruses, a programmed −1 ribosomal frameshift could be used by the virus to produce a fusion protein. Computational analyses of the frameshift signal predicted the presence of an mRNA pseudoknot containing three double-stranded RNA stem structures rather than two. Phylogenetic analyses showed the conservation of potential three-stemmed pseudoknots in the frameshift signals of all other coronaviruses in the GenBank database. Though the presence of the three-stemmed structure is supported by nuclease mapping and two-dimensional nuclear magnetic resonance studies, our findings suggest that interactions between the stem structures may result in local distortions in the A-form RNA. These distortions are particularly evident in the vicinity of predicted A-bulges in stems 2 and 3. In vitro and in vivo frameshifting assays showed that the SARS-CoV frameshift signal is functionally similar to other viral frameshift signals: it promotes efficient frameshifting in all of the standard assay systems, and it is sensitive to a drug and a genetic mutation that are known to affect frameshifting efficiency of a yeast virus. Mutagenesis studies reveal that both the specific sequences and structures of stems 2 and 3 are important for efficient frameshifting. We have identified a new RNA structural motif that is capable of promoting efficient programmed ribosomal frameshifting. The high degree of conservation of three-stemmed mRNA pseudoknot structures among the coronaviruses suggests that this presents a novel target for antiviral therapeutics.     

Nucleomorph genome diversity and its phylogenetic implications in cryptomonad algae

The relationship between phylogeny and nucleomorph genome size was examined in 16 strains of cryptomonad algae using pulsed‐field gel electrophoresis, Southern hybridization and phylogenetic analyses. Our results suggest that all cryptomonads examined in this study contain three nucleomorph chromosomes and their total genome size ranges from 495 to 750 kb. In addition, we estimated the plastid genome size of the respective organisms. The plastid genomes of photosynthetic strains were approximately 120–160 kb in size, whereas the non‐photosynthetic Cryptomonas paramecium NIES715 possesses a genome of approximately 70 kb. Phylogenetic analysis of the nuclear small subunit ribosomal DNA (SSU rDNA) gene showed that nucleomorph genome size varies considerably within closely related strains. This result indicates that the reduction of nucleomorph genomes is a rapid phenomenon that occurred multiple times independently during cryptomonad evolution. The nucleomorph genome sizes of Cryptomonas rostratiformis NIES277 appeared to be approximately 495 kb. This is smaller than that of Guillardia theta CCMP327, which until now was thought to have the smallest known nucleomorph genome size among photosynthetic cryptomonads.     

Enzymatic activity of the SARS coronavirus main proteinase dimer

The enzymatic activity of the SARS coronavirus main proteinase dimer was characterized by a sensitive, quantitative assay. The new, fluorogenic substrate, (Ala-Arg-Leu-Gln-NH)(2)-Rhodamine, contained a severe acute respiratory syndrome coronavirus (SARS CoV) main proteinase consensus cleavage sequence and Rhodamine 110, one of the most detectable compounds known, as the reporter group. The gene for the enzyme was cloned in the absence of purification tags, expressed in Escherichia coli and the enzyme purified. Enzyme activity from the SARS CoV main proteinase dimer could readily be detected at low pM concentrations. The enzyme exhibited a high K(m), and is unusually sensitive to ionic strength and reducing agents.