Genomic variance of Open Reading Frames (ORFs) and Spike protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2).

BACKGROUND: The outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused severe pneumonia at December 2019. Since then, it has been wildly spread from Wuhan, China, to Asia, European, and United States to become the pandemic worldwide. Now coronavirus disease 2019 were globally diagnosed over 3 084 740 cases with mortality of 212 561 toll. Current reports variants are found in SARS-CoV-2, majoring in functional ribonucleic acid (RNA) to transcribe into structural proteins as transmembrane spike (S) glycoprotein and the nucleocapsid (N) protein holds the virus RNA genome; the envelope (E) and membrane (M) alone with spike protein form viral envelope. The nonstructural RNA genome includes ORF1ab, ORF3, ORF6, 7a, 8, and ORF10 with highly conserved information for genome synthesis and replication in ORF1ab.
METHODS: We apply genomic alignment analysis to observe SARS-CoV-2 sequences from GenBank (http://www.ncbi.nim.nih.gov/genebank/): MN 908947 (China, C1); MN985325 (United States: WA, UW); MN996527 (China, C2); MT007544 (Australia: Victoria, A1); MT027064 (United States: CA, UC); MT039890 (South Korea, K1); MT066175 (Taiwan, T1); MT066176 (Taiwan, T2); LC528232 (Japan, J1); and LC528233 (Japan, J2) and Global Initiative on Sharing All Influenza Data database (https://www.gisaid.org). We adopt Multiple Sequence Alignments web from Clustalw (https://www.genome.jp/tools-bin/clustalw) and Geneious web (https://www.geneious.com.
RESULTS: We analyze database by genome alignment search for nonstructural ORFs and structural E, M, N, and S proteins. Mutations in ORF1ab, ORF3, and ORF6 are observed; specific variants in spike region are detected.
CONCLUSION: We perform genomic analysis and comparative multiple sequence of SARS-CoV-2. Large scaling sequence alignments trace to localize and catch different mutant strains in United possibly to transmit severe deadly threat to humans. Studies about the biological symptom of SARS-CoV-2 in clinic animal and humans will be applied and manipulated to find mechanisms and shield the light for understanding the origin of pandemic crisis.

1. INTRODUCTION

The outbreak of severe acute respiratory syndrome coronavirus 2 (n class="Species">SARS-CoV-2) has caused severe pneumonia at December 2019.[1] Since then, it has been wildly spread from Wuhan, China, to Asia, European, and United States to become pandemic worldwide.[2] Severe cases beginning from Huanan Seafood Wholesale market in China which confirmed human pneumonia with the infection of a novel coronavirus (2019-nCoV),[3] and named as SARS-CoV-2 by International Committee on Taxonomy of Viruses.[4,5] Now coronavirus disease 2019 were globally diagnosed over 3 084 740 cases with mortality of 212 516 toll.[6]

Current reports single nucleotide variants are found in many patients with n class="Species">SARS-CoV-2, which belongs to beta-coronavirus species. SARS-CoV-2 contains functional genomic ribonucleic acid (RNA) to transcribe into structural proteins as transmembrane spike (S) glycoprotein for mediating the virus to entry the host cell by utilizing host’s cellular angiotensin-converting enzyme 2 (ACE2), and the nucleocapsid (N) protein holds the major nuclear viral RNA genome; the envelope (E) and membrane (M) alone with spike protein form viral envelope.[7] The nonstructural RNA genome including ORF1ab, ORF3, ORF6, 7a, 8, and ORF10 contains highly conserved information for genome RNA synthesis and replication in ORF1ab and unclear-verified function in other ORF proteins.[8]

The transmission mechanisms with the start of SARS-CoV attaches host cell n class="Gene">membrane receptor and then induce the membrane endocytosis to entry host cells. ORF1 of virus genome leads its replication and synthesize the subgenomic RNAs afterward. Meanwhile, N protein and new genomic RNA assemble to form helical nucleocapsids with M protein inserted in endoplasmic reticulum (ER) and anchored Golgi of host cells. E and M proteins then begin to trigger budding processes. S together with helical N on membrane-bound ER triggers the translation-required viral structure proteins and transport to Golgi. During the final cycle, virions are released by exocytosis to finish the life cycle and replication of the virus.[9]

Previous SARS-CoV-1 in 2003 traene">nsmits possibly through Bat aene">nd Civet as its intn class="Gene">ermediate hosts, and finally to human with the symptoms of severe respiratory impacts in a 10% mortality rate. However, Wuhan SARS-CoV-2 is suspected to be transmitted from bat (RaTG13) to pangolin as intermediate hosts before transmitted to humans by some unknown mechanisms with symptoms of severe respiratory impacts with highest mortality now.[10] The genomic sequence of RaTG13 cited the 96% similarity with Wuhan coronavirus.[11] Although intermediate host is not clear at present, genomic sequence comparison obviously points out spike receptor-binding domain (RBD) of Wuhan SARS-CoV-2 with the similarity in 90% homolog of pangolin. Thus, the possibility that pangolin might contribute the spike protein region to cross-transmitted to RATG13 forms a new recombinant mutant Wuhan SARS-CoV-2 to transmit onto human finally.[12]

The S protein of SARS-CoV-1 aene">nd n class="Species">SARS-CoV-2 responsible for viral entry mediates the binding to host cell membrane of ACE2 through its RBD.[13] The surface S spike protein of SARS-CoV comprises two components (S1 and S2). The S protein of SARS-CoV-2 binds to the host receptor ACE2 through its S1 subunit, which contains RBD, and follows by fusing the viral and host membranes through the S2 subunit, which contains the fusion peptide primed by host protease.

Major six ORFs exist in SARS-CoV-2. n class="Gene">ORF1ab occupies the two-thirds length of the whole genome and subgenome RNA to play roles in viral pathogenesis excluding its replication function as well as involving in cellular signaling and modification of cellular gene expression.[14]

There is no clue for aene">ntiviral thn class="Gene">erapy and treatment for SARS-CoV-2 at present. Further study approaches the molecular genomic variants for selection and packaging is critical for developing antiviral strategies. We will verify and compare various SARS-CoV-2 sequences from different countries by analyzing the possible genomic networks of disease from its origin to evolution, providing the moving development of strategy against the worldwide SARS-CoV-2 pandemic threat.

2. METHODS

2.1. Sequence resource

Studies focus on evolutionary and phylogenetic analysis have applied in disease progression for Wuhan lung pneumonia treatmeene">nt. Hn class="Gene">erein, we apply genomic analysis to observe SARS-CoV-2 sequences from GenBank (http://www.ncbi.nim.nih.gov/genebank/): MN 908947 (China, C1); MN985325 (United States: WA, UW); MN996527 (China, C2); MT007544 (Australia: Victoria, A1); MT027064 (United States: CA, UC); MT039890 (South Korea, K1); MT066175 (Taiwan, T1); MT066176 (Taiwan, T2); MT192759 (Taiwan, T3); MT198652 (Spain, SP); LC528232 (Japan, J1); LC528233 (Japan, J2); MT093571 (Sweden, SW); MT066156 (Italy, IT); and MT050493 (India, In) for genomic sequence alignment analysis.

2.2. Method applied

Multiple Sequence Alignment by Clustalw (https://www.genome.jp/tools-bin/clustalw) web is applied as our alignment tool. Phylogenetic analysis platform performs at Geene">neious website (https://www.geene">neious.com).

3. RESULTS

3.1. ORF1ab

ORF1ab joins 16 proteins togethn class="Gene">er to perform viral genomic replication and synthesis. From the data analysis, it reveals eight mutations from a different country: During this long 6796 amino acids protein, we observe eight mutations located in different regions from various countries; position T609I mutation in California/United States sequence, G818S in Sweden and India, M902I in Korea, F3071Y in Spain, S3120L China, L3606X in Italy and L3606F in Japan, F4321L in Sweden and India, and T6891M in Korea.

3.2. ORF3a

ORF3a functions as accessory protein to help new viral synthesis and escape from the host cell. We find four position mutations; n class="Mutation">M128L in Korea, K136X in Spain, G196V in Spain, and G251V in Italy, Korea, and Sweden.

3.3. ORF6, ORF7a, ORF8, and ORF10

There are no mutations in n class="Gene">ORF6, ORF7a, and ORF10, but we do find one mutation in ORF8 located at L84S from Spain, India, and China.

3.4. E protein

E protein has a short and hydrophilic N-terminus coene">nsisting of 7-12 amino acids, followed by a large hydrophobic traene">nsn class="Gene">membrane domain of 25 amino acids, and ends with a long, hydrophilic C-carboxyl terminus (C-terminal), which comprises the majority of the E protein. Analyzing of E protein alignment, we find one amino acid mutation at L37H from Korea.

3.5. M and N protein

The M protein abundantly defines the shape of the viral envelope. N protein functions primarily to bind to RNA genome of SARS-CoV, making up the nucleocapsid.[15] Although N is most involved in processes viral geene">nome sigene">naling, it is also involved RNA replicatioene">n cycle with host cellular respoene">nse to n class="Disease">viral infection. Although many differences between SARS-CoV-1 and SARS-CoV-2 within in M and N protein, there is no variant observed in M protein but we find a point mutation S197L from Spain.

3.6. S protein

S protein mediates the attachment of SARS-CoV-1 to the host cell surface receptors aene">nd subsequeene">ntly fuse them to facilitate viral eene">ntry into the host cell.[15] The expressioene">n of S protein at the cell n class="Gene">membrane can mediate cell-cell fusion. This formation offers a strategy to spread the virus between cells to subvert function of virus-neutralizing antibodies mechanisms, which play major controlling of protein interaction. By analysis of S protein, we find four mutations from 10 countries; S221W in Korea, S247R in Australia, F737C in Sweden, and A870V in India (Figs. 3–6).

Fig. 3

Genomic analysis of ORF6, ORF7a, ORF8, and ORF10 protein amino acid sequence. There are not any mutations in ORF6, ORF7a, and ORF10, but we find one mutation in ORF8 located at L84S in Taiwan, United States, Spain, India, and China.

Fig. 6

Genomic analysis of S protein amino acid sequence. During analysis of S protein, we find four mutations from 10 countries; S221W in Korea, S247R in Australia, F737C in Sweden, and A870V in India.

4. DISCUSSION

4.1. Point mutation

Six ORFs in SARS-CoV-2 fuene">nctioene">n variously. n class="Gene">ORF1ab joins 16 proteins together to perform viral genomic replication and synthesis. Our first finding reveals eight mutations in different countries. Eight mutation in different regions from various countries are; position T609I mutation in California/United States sequence, G818S in Sweden and India, M902I in Korea, F3071Y in Spain, S3120L China, L3606X in Italy and L3606F in Japan, F4321L in Sweden and India, and T6891M in Korea. No direct evidence proves if each mutant will enhance or decrease viral RNA polymerase and replication (Fig. 1).

Fig. 1

Genomic analysis of ORF1ab protein amino acid sequence. We detect eight mutations in different regions from various countries, T609I mutation in United States, G818S in Sweden and India, M902I in Korea, F3071Y in Spain, S3120L in China, L3606X in Italy and L3606F in Japan, F4321L in Sweden and India, and T6891M in Korea.

Genomic analysis of ORF1ab protein amino acid sequeene">nce. We detect eight mutatioene">ns in diffn class="Gene">erent regions from various countries, T609I mutation in United States, G818S in Sweden and India, M902I in Korea, F3071Y in Spain, S3120L in China, L3606X in Italy and L3606F in Japan, F4321L in Sweden and India, and T6891M in Korea.

ORF3a fuene">nctioene">ns as accessory protein to help new viral syene">nthesis aene">nd escape from the host cell. We find four positioene">n mutatioene">n; n class="Mutation">M128L in Korea, K136X in Spain, G196V in Spain, and G251V in Italy, Korea, and Sweden (Fig. 2). We do not observe any mutations in ORF6, ORF7a, and ORF10 proteins, but we find one mutation in ORF8, which located at L84S from Spain, India, and China. No inclusion can explain the mutations happened at present (Fig. 3).

Fig. 2

Genomic analysis of ORF3a protein amino acid sequence. We find four position mutations; M128L in Korea, K136X in Spain, G196V in Spain, and G251V in Italy, Korea, and Sweden.

Genomic analysis of ORF3a protein amino acid sequeene">nce. We find four positioene">n mutatioene">ns; n class="Mutation">M128L in Korea, K136X in Spain, G196V in Spain, and G251V in Italy, Korea, and Sweden.

Genomic analysis of ORF6, n class="Gene">ORF7a, ORF8, and ORF10 protein amino acid sequence. There are not any mutations in ORF6, ORF7a, and ORF10, but we find one mutation in ORF8 located at L84S in Taiwan, United States, Spain, India, and China.

In comparison of 10 strains from differeene">nt couene">ntries, oene">ne mutatioene">n of E protein is obsn class="Gene">erved at L37H in Korea (Fig. 4). Inside the envelope, there is the nucleocapsid, which is formed from multiple copies of the nucleocapsid (N) protein, which are bound to the positive-sense single stranded RNA genome in a continuous beads-on-a-string type conformation.[16] The lipid bilayer envelope, membrane proteins, and nucleocapsid protect the virus when it is outside the host cell.[17]

Fig. 4

Genomic analysis of E protein amino acid sequence. We found one amino acid mutation at position 37th L37H as “H” from South Korea comparing the “L” from other nine sequences. Yellow line indicates the difference in 10 sequence alignment.

Genomic analysis of E protein amino acid sequence. We found one amino acid mutation at position 37th L37H as “H” from South Korea comparing the “L” from othn class="Gene">er nine sequences. Yellow line indicates the difference in 10 sequence alignment.

Although the N protein holds the viral RNA, and M protein joins with E and S proteins together to create the viral eene">nvelope for protectioene">n wheene">n it is outside the host cell, we do not find point mutatioene">n of M protein.

We do find a point mutation S197L of N protein in Spain. The binding of M to N stability the nucleocapsid (N protein-RNA complex), as well as the intn class="Gene">ernal core of virions, and, ultimately, promotes completion of viral assembly.[18] No evidence demonstrates if S197L will abolish function of N protein (Fig. 5).

Fig. 5

. Genomic analysis of M and N protein amino acid sequences. We do not observe any mutation in 10 sequences of M protein region but detect one mutation in Spain at S197L of N protein.

. Genomic analysis of M and N protein amino acid sequences. We do not observe aene">ny mutatioene">n in 10 sequeene">nces of M protein regioene">n but detect oene">ne mutatioene">n in Spain at n class="Mutation">S197L of N protein.

By analysis of S protein, we find four mutations from 10 countries; S221W in Korea, n class="Mutation">S247R in Australia, F737C in Sweden, and A870V in India (Fig. 6). Report[19] mentioned a single amino acid reversion (L294Q) in the S protein is sufficient to abrogate the phenotype and grows well at and below 32oC.

Genomic analysis of S protein amino acid sequence. During analysis of S protein, we find four mutations from 10 countries; S221W in Korea, n class="Mutation">S247R in Australia, F737C in Sweden, and A870V in India.

4.2. Large scaling alignment of spike protein mutations and phylogenetic analysis

Although SARS-CoV-1 aene">nd n class="Species">SARS-CoV-2 share the sequence similarity with 80% homolog. After performing the alignment, they reveal their 75% similarity in spike protein. The S protein mediates viral entry into host cells by first binding to a host receptor through the RBD in the S1 subunit and then fusing the viral and host membranes through the S2 subunit priming by host cell proteases.[20-23] Unraveling which cellular factors are used by SARS-CoV-2 for entry might provide insights into viral transmission and reveal therapeutic targets. SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV) RBDs recognize different receptors. SARS-CoV recognizes ACE2 as its receptor, whereas MERS-CoV recognizes dipeptidyl peptidase 4 as its receptor.[14,24] Since SARS-CoV-2 recognizes ACE2 as its host receptor binding to viral S protein.[25] Therefore, it is critical to define the RBD in SARS-CoV-2 S protein as the most likely target for the mechanism of virus attachment such as new developing inhibitors, neutralizing antibodies, and vaccines.

Authors from the group of Tai et al[26] demonstrate by characterizing of n class="Species">SARS-CoV-2 RBD to display a multiple sequence alignment of RBDs of SARS-CoV-2, SARS-CoV, and MERS-CoV spike (S) proteins.

They identified the RBD in SARS-CoV-2 S protein aene">nd fouene">nd that the RBD protein bouene">nd stroene">ngly to n class="Species">human and bat ACE2 receptors. SARS-CoV-2 RBD displayed significantly higher binding affinity to ACE2 receptor than SARS-CoV RBD. Subsequently, SARS-CoV RBD-specific antibodies could cross-react with SARS-CoV-2 RBD protein. Meanwhile, SARS-CoV RBD-induced antisera could cross-neutralize SARS-CoV-2 which suggested the potentials to develop SARS-CoV RBD-based vaccines for prevention of SARS-CoV-2 and SARS-CoV infection.[26]

Hoffmann group mentions SARS-CoV-1 aene">nd n class="Species">SARS-CoV-2 share 76% amino acid identity in spike protein region. By the amino acid alignment, they observe the receptor-binding motif of SARS-CoV-1 corresponding to the sequences of bat-associated beta-coronavirus S proteins. Demonstration of high or low similarity by taking advantage of ACE2 as cellular receptor reveals SARS-CoV-2 possesses crucial amino acid residues for ACE2 binding.

They also find similarity signal to points out between SARS-CoV-2 aene">nd n class="Species">SARS-CoV-1 during transmitting host cells stage and then identify a potential target for antiviral intervention. Inspecting conserved amino acids within ACE2 domain, Hoffmann group perform SARS-CoV-2 to transmit cell entry depends on ACE2 and transmembrane serine protease 2 two proteins and is blocked by applied clinically proven protease inhibitor.[27,28]

By deep and large scaling analysis of spike protein from maene">ny couene">ntries, we do have variaene">nts fouene">nd in US case including specimeene">n from east coast Uene">nited States. We do find variaene">nts in Uene">nited States comparing with China origin (Fig. 7). Mutaene">nt-1 expresses a “G” amino acid at 614 instead of China “D” (n class="Mutation">D614G). Mutant-2 strain displays the position at 614 same as China strain with “D” but other mutations found in different regions (Fig. 8A). Mutant 2-2 with same position of 614 “D” but only display one mutation same as China pointed as QIS60546 strain (Fig. 8B). Studies suggest various viral strains originally spread from China to Europe which one strain should be deadly mutations as observed and then they spread to New York finally. The other milder strains also spread to west coast in United States from China.[29] Since this report cites SARS-CoV2 acquired mutations capable of substantially changing its pathogenicity. Will this observation be matched with our finding that three variants found in New York become more severe transmitted to humans than west coast in the United States?

Fig. 7

Spike protein reveals variants in the world. We find many variants in spike protein by alignment and phylogenetic analysis.

Fig. 8

A, Spike protein in China sequences exhibit a conserved amino acid. We found a conserved amino acid “D” at position 614 of Spike protein in most China sequences. B, Analysis indicates three variants of spike protein in the United States. We observe three variants in the analysis of United States sequences; mutant-1 found with different amino acid “G” at position 614, mutant 2-1 with same “D” at position 614 same as China but various variants at other regions. Mutant 2-2 same as 2-1 at 614 but same as China in one region as QIS60546 indicated. (I) US case. (II) Phylogenetic analysis to map three mutants in United States and China.

Spike protein reveals variants in the world. We find many variants in n class="Gene">spike protein by alignment and phylogenetic analysis.

A, Spike protein in China sequeene">nces exhibit a coene">nsn class="Gene">erved amino acid. We found a conserved amino acid “D” at position 614 of Spike protein in most China sequences. B, Analysis indicates three variants of spike protein in the United States. We observe three variants in the analysis of United States sequences; mutant-1 found with different amino acid “G” at position 614, mutant 2-1 with same “D” at position 614 same as China but various variants at other regions. Mutant 2-2 same as 2-1 at 614 but same as China in one region as QIS60546 indicated. (I) US case. (II) Phylogenetic analysis to map three mutants in United States and China.

Limitedly in the study, we perform our study eithn class="Gene">er data mining by alignment and phylogenetic analysis from public domains such as Global Initiative on Sharing All Influenza Data and National Center for Biotechnology Information. There will be interesting to demonstrate biological approaches with specimens in hands to observe the correlation from clinical to lab analysis directly.

In conclusion, we analyze database by genome alignment search for nonstructural ORFs and structural E, M, N, and S proteins. Large scaling performaene">nce to catch diffn class="Gene">erent mutant strains in American possibly induce severe deadly threat to humans. More studies about the biological symptom of SARS-CoV-2 in clinic animal and humans will manipulate and shield the light for understanding the origin of pandemic crisis.

ACKNOWLEDGMENTS

This research was funded by Taipei Veteraene">ns Geene">nn class="Gene">eral Hospital (grant number V107E-002-2, V108D46-004-MY2-1, V108E-006-4, 108E-006-5, and 109VACS-003).