
==== Front
F1000Res
F1000Res
F1000Research
2046-1402
F1000 Research Limited London, UK

10.12688/f1000research.72896.4
Research Article
Articles
Prediction of the effects of the top 10 synonymous mutations from 26645 SARS-CoV-2 genomes of early pandemic phase
[version 4; peer review: 1 approved

Boon Wan Xin Formal Analysis Investigation Visualization Writing – Original Draft Preparation Writing – Review & Editing https://orcid.org/0000-0003-0424-5243
1
Sia Boon Zhan Investigation Methodology 1
Ng Chong Han Conceptualization Funding Acquisition Project Administration Supervision Writing – Original Draft Preparation Writing – Review & Editing https://orcid.org/0000-0003-2926-9831
a1
1 Faculty of Information Science and Technology, Multimedia University, Bukit Beruang, Melaka, 75450, Malaysia
a chng@mmu.edu.my
No competing interests were disclosed.

18 9 2024
2021
10 105311 9 2024
Copyright: © 2024 Boon WX et al.
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Background

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) had led to a global pandemic since December 2019. SARS-CoV-2 is a single-stranded RNA virus, which mutates at a higher rate. Multiple works had been done to study nonsynonymous mutations, which change protein sequences. However, there is little study on the effects of SARS-CoV-2 synonymous mutations, which may affect viral fitness. This study aims to predict the effect of synonymous mutations on the SARS-CoV-2 genome.

Methods

A total of 26645 SARS-CoV-2 genomic sequences retrieved from Global Initiative on Sharing all Influenza Data (GISAID) database were aligned using MAFFT. Then, the mutations and their respective frequency were identified. Multiple RNA secondary structures prediction tools, namely RNAfold, IPknot++ and MXfold2 were applied to predict the effect of the mutations on RNA secondary structure and their base pair probabilities was estimated using MutaRNA. Relative synonymous codon usage (RSCU) analysis was also performed to measure the codon usage bias (CUB) of SARS-CoV-2.

Results

A total of 150 synonymous mutations were identified. The synonymous mutation identified with the highest frequency is C3037U mutation in the nsp3 of ORF1a. Of these top 10 highest frequency synonymous mutations, C913U, C3037U, U16176C and C18877U mutants show pronounced changes between wild type and mutant in all 3 RNA secondary structure prediction tools, suggesting these mutations may have some biological impact on viral fitness. These four mutations show changes in base pair probabilities. All mutations except U16176C change the codon to a more preferred codon, which may result in higher translation efficiency.

Conclusion

Synonymous mutations in SARS-CoV-2 genome may affect RNA secondary structure, changing base pair probabilities and possibly resulting in a higher translation rate. However, lab experiments are required to validate the results obtained from prediction analysis.

COVID-19
SARS-CoV-2
Synonymous mutations
RNA secondary structure
Multimedia University, Malaysia IRFund 2.0MMUI/210119 This research is supported by Multimedia University, Malaysia, IRFund 2.0 (grant number MMUI/210119 awarded to Chong Han, Ng). The funder has no role in study design, data analysis, decision to publish or manuscript preparation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Revised Amendments from Version 3

The title of the paper has been revised to, “Prediction of the effects of the top 10 synonymous mutations from 26645 SARS-CoV-2 genomes of early pandemic phase”. The introduction section with information on Alpha variant has been added and some introduction section has been revised. Additional paragraphs in the discussion section on identification of synonymous mutation, RNA secondary structure prediction and codon bias usage have been included to improve the clarity of the manuscript. Some citations have been removed, added or updated accordingly.
==== Body
pmcIntroduction

In December 2019, coronavirus disease 2019 (COVID-19) cases first emerged from Wuhan, China 1 . Soon after, rapid spread of COVID-19 has resulted in a serious global outbreak. COVID-19 is an infectious and potentially lethal disease caused by a newly found coronavirus strain, known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The virus causes clinical manifestation ranging from asymptomatic to severe pneumonia and in the worst scenario, death 2 . SARS-CoV-2 seems to have a higher transmission rate 3 but lower mortality rate 2 in comparison to Middle East respiratory coronavirus (MERS-CoV) and severe acute respiratory syndrome coronavirus (SARS-CoV).

SARS-CoV-2 is a single-stranded RNA virus with a genome size of 29,903 bases. In general, RNA viruses have a higher mutation rate than DNA viruses and this allows them to evolve rapidly, escaping the host immune defence response 4 . Different SARS-CoV-2 variants with multiple synonymous and nonsynonymous mutations have been reported since the beginning of the outbreak 5 . Some variants are classified as variants of concern (VOCs) since they are associated with the change in viral pathogenicity such as, higher disease severity, higher transmission rate, lower immunity response in the host as a consequence of the mutations 6 . However, it is expected that most of these mutations in SARS-CoV-2 genome are either neutral or mildly deleterious 7 . Numerous studies have been carried out to understand the molecular mechanisms of these nonsynonymous mutations on the functions of different SARS-CoV-2 proteins 6 . For example, the Alpha variant (B.1.1.7) of SARS-CoV-2, first identified in the UK in late 2020, is characterized by several mutations, including the D614G mutation in the spike protein, which enhances its binding affinity to the ACE2 receptor 8 . This variant exhibits increased transmissibility compared to the original virus, which led to its rapid spread globally 9 . However, there are only a few studies on the synonymous mutations of SARS-CoV-2 genome 10, 11 .

Synonymous mutations are also known as silent mutations because the nucleotide mutations result in a change in the RNA sequence without altering the amino acid sequence 12 . Synonymous mutations have been suggested to have no functional consequence on the fitness of organisms and their evolution in long term 13 . However, numerous recent studies had showed that synonymous mutations may affect the folding and stability of RNA structures 14 . Interestingly a large scale study of synonymous mutations in multiple yeast genes has shown that most of synonymous mutations are not neutral, affecting the fitness of the cell 15 . For RNA viruses, even though synonymous mutations generally do not change their pathogenicity directly, some studies reveal that synonymous mutations may affect the RNA secondary structure of the virus 16 and also change the codon usage bias of the genes in the virus 17, 18 . The use of mRNA-based COVID-19 vaccines reduce the severity of the disease. However, mRNA molecule is susceptible to the degradation due to the presence of 2’ OH group in the ribose. To improve the stability of mRNA vaccine, Zhang et al. (2023) designed a novel algorithm, which optimizes the codon usage and RNA secondary structure by using synonymous codon 19 .

The synonymous mutations play some important biological roles, which may affect viral fitness and pathogenicity. However, the study of biological consequences of synonymous mutations have been largely overlooked. In this study, we identified synonymous mutations of SARS-CoV-2 genome from early pandemic phase. We predicted the effects of these synonymous mutations of the top 10 highest frequency on RNA secondary structures and codon usage bias of SARS-CoV-2 genome. These findings allow the researchers to prioritize these mutations for function analysis in the future.

Methods

Sequence retrieval

30,229 SARS-CoV-2 genomic sequences were downloaded from GISAID database (Global Initiative on Sharing All Influenza Data, RRID:SCR_018251) 20 ranging from 31 December 2019 to 22 March 2021. SARS-CoV-2 genomic sequences were filtered by setting parameters to keep only sequences with complete genome and high coverage. The sequences were further filtered to remove those sequences with higher than 0.1% “N” unresolved nucleotides and ambiguous letters. A total of 3,584 sequences were removed by applying this filter. The reference sequence of SARS-CoV-2 genome ( NC_045512.2) 21 was retrieved in fasta format from NCBI database (NCBI, RRID:SCR_006472). It is a Wuhan isolate with a complete genome which comprises of 29,903 bases.

Multiple sequence alignment

The rapid calculation available in MAFFT online server (MAFFT, version 7.467, RRID:SCR_011811) 22 was used to perform multiple sequence alignment (MSA) for 26,645 SARS-CoV-2 genomes. This option supports the alignment of more than 20,000 sequences with approximately 30,000 sites. The alignment length was kept, which means the insertions at the mutated sequences were removed, to keep the alignment length the same as the reference sequence. While other parameters were left as default.

Identification of mutations and their frequency in SARS-CoV-2 genomes

A simple Python script was written to identify the mutations in 26,645 SARS-CoV-2 genomes. To determine whether the identified mutations are synonymous or nonsynonymous, MEGA X software, version 10.2.5 build 10210330 ( MEGA Software, RRID:SCR_000667) 23 was utilized to perform the translation for inspection purposes. The presence of amino acid changes was identified by referring to the genomic position of the nucleotide mutations. Synonymous mutations with the top 10 highest frequencies were generated.

SARS-CoV-2 RNA secondary structure prediction

The RNA secondary structure of wild type and mutant sequences were predicted using RNAfold program, version 2.4.18 (Vienna RNA, RRID:SCR_008550) 24 with the incorporation of SHAPE reactivity data obtained from the study done by Manfredonia et al. (2020) 25 . The RNA secondary structure prediction was performed using a sequence length of 250 nucleotides upstream and downstream of the mutation site. Other than RNAfold, another two programs which are IPknot++ version 2.2.1 (SCR_022557) 26 , and MXFold2 (SCR_022558) 27 were also used to perform the RNA secondary structure prediction of SARS-CoV-2 wild type and mutants.

Base pair probability estimation

To predict how the mutations affect RNA local folding, base pair probability was estimated by utilizing MutaRNA, version 1.3.0 (MutaRNA, RRID:SCR_021723) 28 . MutaRNA is a web-based tool that allows prediction and visualization of the structure changes induced by a single nucleotide polymorphism (SNP) in an RNA sequence. It includes the base pair probabilities within RNA molecule of both wild type and mutant. The parameters used in MutaRNA were set as default except the window size was changed to 501nt.

Relative Synonymous Codon Usage (RSCU)

Relative synonymous codon usage (RSCU) represents the ratio of the observed frequency of codons appearing in a gene to the expected frequency under equal codon usage. RSCU is calculated using the formula:

RSCUi=Xi1n∑i=1nXi,

where X i implies the number of occurrences of codon i and n stands for the number of synonymous codons encoded for that particular amino acid.

Results and discussion

A synonymous mutation is a change in the nucleotide that does not cause any changes in the encoded amino acid. Synonymous mutations were previously considered to be less important, but they are now proven to have some effects on RNA folding, RNA stability, miRNA binding and translational efficiency 29 . Synonymous mutations may have significant effects on the adaptation, virulence, and evolution of RNA viruses 30 . Another study done also indicated that synonymous mutations have association with more than 50 human diseases such as hemophilia B, tuberculosis (TB), cystic fibrosis (CF), Alzheimer, schizophrenia, chronic hepatitis C and so on 31 . All these studies show that increasing importance has been associated with synonymous mutations over these years. Hence, it is necessary for us to study the effects of synonymous mutations of SARS-CoV-2 genome.

Identification of SARS-CoV-2 synonymous mutations

A total of 381 mutations were found in SARS-CoV-2 genomes by using python script, in which 150 of them are synonymous mutations. The distribution of these 150 synonymous mutations in 11 coding regions is shown in Figure 1. Among these mutations, ORF1a and ORF1b have a higher number of synonymous mutations at 76 and 33, respectively, which might be due to their longer sequence length. Besides that, our findings also show high C to U mutation rate in SARS-CoV-2 genome and this mutational skews are in line with multiple studies 32– 35 . The high C to U mutation rate may be driven by host APOBEC-mediated RNA editing system and overexpression of APOBEC3 protein promotes viral replication and propagation in the human colon epithelial cell line 36 . These mutational skews are necessary to be considered when deducing the selection acting on synonymous variants in SARS-CoV-2 evolution 11 . Synonymous mutations are assumed subject to a lower selective pressure than nonsynonymous mutations, presumably the purifying selection force has stronger negative impact on the frequencies of nonsynonymous mutations. Interestingly there may be some selection force on synonymous mutations shown by a few studies, suggesting that these synonymous mutations are not random and neutral, may have some biological impact on viral fitness 11, 32, 37 .

Figure 1. Distribution of SARS-CoV-2 synonymous mutations in 11 coding regions.

The synonymous mutations in SARS-CoV-2 genomes with the top 10 highest frequency obtained from the analysis of 150 synonymous mutations were listed in Table 1. Our sequence samples are obtained from December 2019 to March 2021 and this period overlapped with the peak of Alpha variant (B.1.1.7) outbreak 5 . The defining synonymous SNPs of Alpha variant include C241T, C913T, C3037T, C5986T, C14676T, C15279T and T16176C 5 , and all except C241T are reported in our study as well. As shown in Table 1, synonymous mutations with the highest frequency identified from SARS-CoV-2 genomes is C3037U mutation located in nsp3 of ORF1a, followed by C313U mutation in nsp1 of ORF1a and C9286U mutation in nsp4 of ORF1a. Mutations with higher frequency are mostly found in ORF1a and ORF1b. Although there are some overlapping ORFs in the SARS-CoV-2 genome, such as ORF1a and ORF1b, ORF3a and ORF3c 38 , the top 10 highest frequency synonymous mutations are not located in these overlapping sites. It is of great interest to find out the effect of these top 10 synonymous mutations on SARS-CoV-2 genome. However, it is important to take note that the high frequency of some mutations is not necessarily due to their positive effects. They may emerge during early stage of pandemic and are transmitted to all of their descendants, even though they have no or little effect on viral fitness 39 .

Table 1. SARS-CoV-2 synonymous mutations with the top 10 highest frequency.

ORF	Position	Nucleotide
Variation	Frequency	
1a	nsp1	313	C > U	8212	
nsp2	913	C > U	2820	
nsp3	3037	C > U	24651	
5986	C > U	2860	
nsp4	9286	C > U	3880	
1b	nsp12	14676	C > U	2859	
15279	C > U	2842	
16176	U > C	2838	
nsp14	18877	C > U	2603	
M	26735	C > U	2432	

Similar to another companion paper, which focuses on the prediction analysis of nonsynonymous mutations of SARS-CoV-2 proteins 40 , the same SARS-CoV-2 virus genome data from GISAID database ranging from 1st January 20 to 22 March 21 were used in this study. The data collection time was overlapping with the period when the frequency of alpha variant reached the highest numbers around March–May 21 34 . There are seven synonymous mutations identified as the defining mutations in the alpha variant, of which all except C241T are also reported in our study. Due to the rapid evolution of SARS-CoV-2 genome, it is beyond the scope of our study to keep track SARS-CoV-2 mutational profile and to predict the consequences of these mutations. Two independent studies reported that alpha or alpha-like SARS-CoV-2 variants are circulating among wild deer population in North America in late 2021 41, 42 . Although there is no reported case of viral spillback from deer to human transmission, we can’t simply rule out this possibility yet. Hence, our findings remain relevant despite of not using the latest genome dataset.

RNA secondary structure prediction and base pair probability estimation analysis

SARS-CoV-2 virus can form highly structured RNA elements, which may affect viral replication, discontinuous transcription and translation 43, 44 . For example, SARS-CoV-2 forms a three-stemmed pseudoknot structure to promote programmed -1 ribosomal frameshifting to increase the synthesis of the proteins required for viral replication 43, 44 . There are numerous high throughput studies on the characterization of RNA secondary structure of SARS-CoV-2 genome 25, 45– 48 . In these recent high throughput studies, the RNA secondary structures of SARS-CoV-2 genome were determined experimentally using chemical probing methods, such as SHAPE-MaP 25, 45 or proximity ligation methods, such as RIC-seq 47 , COMRADES 48 . Although these data are very useful to determine the RNA secondary structures of SARS-CoV-2 virus, there is very little study on the effect of the synonymous mutations on RNA secondary structure, which may be beneficial or deleterious to the viral fitness. Therefore, we performed RNA secondary structure prediction and base pair probability estimation analysis of these top 10 highest frequency of synonymous mutations.

To improve the outcome of the study, multiple RNA secondary structure prediction tools, namely RNAfold with SHAPE reactivity data 24 , IPknot++ 26 and MXfold2 27 were applied in our study. In addition, MutaRNA analysis tool was used to estimate the base pair potential of the wild type and mutant sequences. RNAfold with SHAPE reactivity data uses thermodynamic approach to calculate the minimum free energy for the most probable RNA secondary structure by incorporating the nucleotide reactivity data derived from the experiments. If the reactivity value is high, the nucleotide is less likely to be paired, or vice versa. SARS-CoV-2 virus can form pseudoknot structures, which promote ribosomal frameshifting 44 . However, many RNA secondary structure prediction programmes don’t predict pseudoknot structure since the calculation is computationally demanding. IPknot++ is one of the few programmes, which can predict pseudoknot structure. MXfold2 predicts RNA secondary structure using deep learning method with a large amount of training dataset.

Although different tools may produce similar results for identical RNA sequences, it's important to note that there can be variations in prediction outputs due to differences in algorithms, thermodynamic parameter settings, inclusion of pseudoknot calculation, incorporation of experimental data and the assumptions of each tool. RNAfold, IPknot++ and MXfold2 apply the nearest neighbour model, using different thermodynamic parameters. In addition, MXfold2 implements deep learning models with max margin framework 49 . Multiple experimental genome-wide mapping of RNA secondary structures studies showed that 5’ UTR of SARS-CoV-2 RNA genome forms 7 conserved stem-loop structures, and in some studies, 8 SLs, depending on the sequence length 25, 45– 48, 50 . To demonstrate the usability of prediction tools, we predicted RNA secondary structure of the sequence of 5’ UTR (1-480 nt) of SARS-CoV-2 (Extended data 1) 51 . The RNA secondary structures predicted by RNAfold with SHAPE data, IPknot++ and MXfold2 are similar, especially SL1, SL5-8 regions and they are comparable to most of the published experimental data 25, 45– 48, 50 . Both RNAfold with SHAPE data, and MXfold2 successfully predicted SL4, but IPknot++ predicted SL4 with pseudoknot structure, which has not been reported in other studies. Interestingly the result obtained from RNAfold without SHAPE data is quite different, possibly due to the missing experimental data. In addition, it has been shown that SARS-CoV-2 may adopt different RNA secondary structure conformations 7, 19, 36, 37, 39, 41 . Our study is aimed to predict if the sSNP may affect RNA secondary structure and the outcomes allow us to prioritize variants for the experiment functional studies in the future. Using multiple prediction tools may help to increase the accuracy and reliability of the prediction result. The prediction results for all 10 synonymous mutations using these 3 tools and the base pair probability estimation results are summarized in the Table 2 (✓ - changes, × - no change). The results for all 10 synonymous mutations predicted with RNAfold, IPknot++ and MXfold2 are available in Extended data 2, 3 and 4, respectively 51 . The base pair probabilities for all 10 synonymous mutations are shown as circular plots in Extended data 5 51 . The darker the edge is, the more likely the two connected bases to form base pair. Of these 10 synonymous mutations, four mutants which are all located in ORF1ab, namely C913U, C3037U, U16176C and C18877U mutants show pronounced changes between wild types and mutants in all 3 prediction tools, suggesting these synonymous mutations may have some biological impact on viral fitness. Having say that, it is also possible that other mutants with only one or two changes predicted by these analyses, may also affect RNA secondary structures, having some impact on viral fitness. It has been shown that SARS-CoV-2 virus can form elaborated RNA secondary structures at 5’ and 3’UTRs, and frameshifting element (FSE), located between the boundary of ORF1a and ORF1ab 7, 19, 36, 37, 39, 41 . The 5’ UTR of SARS-CoV-2 is important for viral mRNA stability 52 and protein translation 53 while the 3’ UTR may be involved in viral proliferation in the host cell 54 . Interestingly it has been observed that base substitution type, transitions from C to U base occurred at higher frequency in the stem region of RNA secondary structure of 5’ and 3’ UTR of SARS-CoV-2 genome, possibly due to the less detrimental effect on the structure 34 . The FSE can form pseudoknot structures, which regulate the relative protein expression of ORF1a and ORF1ab during viral infection 43, 44 .

Table 2. Summary of RNA secondary structure prediction and base pair probability estimation analysis of SARS-CoV-2 synonymous mutations.

	RNAfold (SHAPE)	IPknot++	MXFold2	MutaRNA	
C313U	×	×	×	×	
C913U	✓	✓	✓	✓	
C3037U	✓	✓	✓	✓	
C5986U	✓	×	✓	✓	
C9286U	×	✓	×	✓	
C14676U	✓	✓	×	✓	
C15279U	×	✓	×	×	
U16176C	✓	✓	✓	✓	
C18877U	✓	✓	✓	✓	
C26735U	✓	✓	×	✓	

Other than 5’ and 3’ UTRs, Huston et al. (2021) found that ORF1ab region forms extensive RNA secondary structure network 45 . Coincidentally all four mutations, C913U, C3037U, U16176C and C18877U reported in our study are located within ORF1ab. C913U mutation is found in the Nsp2, near the start codon (position 806) in ORF1a in SARS-CoV-2 genome. As shown in Figure 2, the wild type structures predicted by RNAfold and MXfold2 shares some degree of similarity around position 95–330 of 501 base long structure. C913U mutation has a pronounced effect on RNA secondary structure predicted by RNAfold. C913U mutation results in the appearance or disappearance of multiple loops, not only at the nearby mutated residue, but also at the sites further apart, suggesting this mutation may affect its long-range RNA interaction. While MXfold2 predicts that U913 mutant forms a shorter stem and a larger hairpin loop compared to C913 wild type. However, the structure predicted by IPknot++ is quite different from others, in which, C913U results in change of pseudoknot structure. Figure 2D shows that the base pair interactions of wild type RNA are changed substantially by U913 mutation. Previously it has been shown that Nsp2 protein suppresses host immune response by inhibiting the mRNA translation of interferon gene 55 . Although C913U mutation does not alter the amino acid residue of Nsp2 protein, it may be worthwhile to see if this C913U mutation plays a direct or indirect role in host immune response through Nsp2 protein. Since C913U is near to Nsp1 and Nsp2 protein boundaries, the altered RNA secondary structure may affect ribosome stalling, which, in turn, affect folding of nascent polypeptides and translation initiation. In addition, Nsp1 protein facilitates viral propagation by inhibiting host protein translation machinery 56 and promoting host mRNA degradation 57 . It will be interesting to investigate if the C913U mutation affects these functions.

Figure 2. The effect of C913U mutation on RNA secondary structure of nsp2 in ORF1a.

( A) RNA secondary structure of C913 wild type and U913 mutant predicted using RNAfold. ( B) RNA secondary structure of C913 wild type and U913 mutant predicted using IPknot++. ( C) RNA secondary structure of C913 wild type and U913 mutant predicted using MXfold2. ( D) MutaRNA circular plots of base pairing probabilities of C913 wild type and U913 mutant. The black arrow indicates the position of WT and mutated nucleotides while the red arrow indicates the starting position of the query sequence.

C3037U mutation is found in the Nsp3 in ORF1a. As shown in Figure 3, both IPknot++ and MXfold2 predict that U3037 mutant forms longer stem and smaller internal loop compared to wild type. On the contrary, RNAfold predicts that a small internal loop fuses into a bigger internal loop in U3037 mutant. MutaRNA circular plot shows that there is some minor difference in base pair probabilities between C3037 wild type and U3037 mutant. Nsp3 is a papain-like protease, which hydrolyzes several Nsp proteins, involved in viral replication 58 . Hence, we should investigate the effect of this mutation on its cleavage activity, probably through the change in transcription or translation level of Nsp3.

Figure 3. The effect of C3037U mutation on RNA secondary structure of nsp3 in ORF1a.

( A) RNA secondary structure of C3037 wild type and U3037 mutant predicted using RNAfold. ( B) RNA secondary structure of C3037 wild type and U3037 mutant predicted using IPknot++. ( C) RNA secondary structure of C3037 wild type and U3037 mutant predicted using MXfold2. ( D) MutaRNA circular plots of the base pairing probabilities of C3037 wild type and U3037 mutant. The black arrow indicates the position of WT and mutated nucleotides while the red arrow indicates the starting position of the query sequence.

U16176C mutation is located in the Nsp12, close to the boundary of Nsp12 and Nsp13 genes in ORF1b. As shown in Figure 4, U16176C mutation results in a drastic change in RNA secondary structure predicted using RNAfold. IPknot++ predicts C16176 mutant forms new pseudoknot structures, which are absent in wild type U16176. On the other hand, MXfold2 predicts C16176 mutant forms a larger multi-branched loop and a shorter stem compared to wild type. Similarly, MutaRNA result shows C16176 mutant affects base pair potential at multiple sites. Nsp12 is one of the subunits of RNA-dependent RNA polymerase (RdRp), which is required for RNA synthesis 59 . A study showed that a 1.4-kb-long SARS-CoV-2 RNA sequence (residues 15071–16451) located in the Nsp12 and Nsp13 regions is required to facilitate viral RNA packaging 60 . Since U16176C mutation may affect RNA secondary structure, it will be interesting to see if it affects viral RNA packaging. U16176C together with C14676U and C15279U have very similar number of frequencies as shown in Table 1. Interestingly IPknot++ predicted all of them result in changes in pseudoknot structure as shown in Extended data 3. We speculated that these three sSNPs may be functionally related. These mutations are located downstream of the frameshifting element (residues 13405–13488) and this element forms a pseudoknot to promote ribosomal frameshifting during viral replication 61 . It has been demonstrated that synonymous mutations affect both RNA secondary structure of the ribosomal frameshift signal and frameshifting efficiency in SARS-CoV virus 62 . Another study had shown that this ribosomal frameshifting structure in SARS-CoV-2 virus involves long-range sequence interaction of 1.5 kb 48 . It remains to be seen whether the long-range sequence interaction for ribosomal frameshifting can go beyond 1.5kb long.

Figure 4. The effect of U16176C mutation on RNA secondary structure of nsp12 in ORF1b.

( A) RNA secondary structure of U16176 wild type and C16176 mutant predicted using RNAfold. ( B) RNA secondary structure of U16176 wild type and C16176 mutant predicted using IPknot++. ( C) RNA secondary structure of U16176 wild type and C16176 mutant predicted using MXfold2. ( D) MutaRNA circular plots of the base pairing probabilities of U16176 wild type and C16176 mutant. The black arrow indicates the position of WT and mutated nucleotides while the red arrow indicates the starting position of the query sequence.

C18877U mutation is located in Nsp14 in ORF1b. As shown in Figure 5, an additional internal loop is formed in U18877 mutant predicted by RNAfold. IPknot++ predicts U18877 mutant forms extra internal loops and longer hairpin near the mutated residue and it also affects the pseudoknot structure at 2 different sites further from the mutated residue. While MXfold2 predicts U18877 mutant forms one hairpin with multiple loops instead of one hairpin as seen in wild type. The changes at multiple base pairing sites due to the U18877 mutation is also observed in MutaRNA circular plot. Nsp14 is important to maintain high fidelity during viral RNA synthesis 63 .

Figure 5. The effect of C18877U mutation on RNA secondary structure of nsp14 in ORF1b.

( A) RNA secondary structure of C18877 wild type and U18877 mutant predicted using RNAfold. ( B) RNA secondary structure of C18877 wild type and U18877 mutant predicted using IPknot++. ( C) RNA secondary structure of C18877 wild type and U18877 mutant predicted using MXfold2. ( D) MutaRNA circular plots of the base pairing probabilities of C18877 wild type and U18877 mutant. The black arrow indicates the position of WT and mutated nucleotides while the red arrow indicates the starting position of the query sequence.

RSCU analysis of SARS-CoV-2

Other than affecting RNA secondary structure, it has been shown that synonymous mutations may affect protein translation efficiency and accuracy through the formation of codon usage bias (CUB), which is non-random usage of synonymous codons, common in all species 64 . It is a phenomenon where some codons are preferred over others for a specific amino acid. SARS-CoV-2 replicates using host cell’s machinery and synthesizes its protein by utilizing host cellular components. Hence, codon usage bias may affect the replication of viruses 65 .

Relative synonymous codon usage (RSCU) is a widely used statistical approach 66 that can be used to measure codon usage bias in coding sequences. The RSCU values of SARS-CoV-2 are shown in Table 3 and the most preferred codons for each amino acid are marked in bold. Stop codons (UAA, UAG, UGA) and codons which code for an amino acid uniquely (AUG, UGG) are excluded from RSCU analysis.

Table 3. RSCU values of SARS-CoV-2 genome.

Amino Acid	Synonymous Codons	RSCU	
Ala	GCA	1.09	
GCC	0.58	
GCG	0.16	
GCU	2.17	
Arg	AGA	2.67	
AGG	0.81	
CGA	0.29	
CGC	0.58	
CGG	0.19	
CGU	1.46	
Asn	AAC	0.65	
AAU	1.35	
Asp	GAC	0.72	
GAU	1.28	
Cys	UGC	0.45	
UGU	1.55	
Gln	CAA	1.39	
CAG	0.61	
Glu	GAA	1.44	
GAG	0.56	
Gly	GGA	0.82	
GGC	0.71	
GGG	0.12	
GGU	2.34	
His	CAC	0.61	
CAU	1.39	
Ile	AUA	0.92	
AUC	0.56	
AUU	1.53	
Leu	CUA	0.66	
CUC	0.59	
CUG	0.30	
CUU	1.75	
UUA	1.63	
UUG	1.06	
Lys	AAA	1.31	
AAG	0.69	
Phe	UUC	0.59	
UUU	1.41	
Pro	CCA	1.59	
CCC	0.29	
CCG	0.17	
CCU	1.94	
Ser	AGC	0.36	
AGU	1.43	
UCA	1.67	
UCC	0.46	
UCG	0.11	
UCU	1.96	
Thr	ACA	1.64	
ACC	0.38	
ACG	0.20	
ACU	1.78	
Tyr	UAC	0.78	
UAU	1.22	
Val	GUA	0.91	
GUC	0.56	
GUG	0.58	
GUU	1.95	

Based on the RSCU values, the synonymous codons can be classified into five groups: i) codons with RSCU value equals to 1.0 are unbiased codons; ii) codons with RSCU value > 1.0 are codons preferred in a genome; iii) codons with RSCU value < 1.0 are codons less preferred in a genome; iv) codons with RSCU value > 1.6 are codons which are over-represented in a genome; v) codons with RSCU value < 1.6 are codons which are under-represented in a genome 65 . There are 15 preferred codons (RSCU value > 1.0) and 11 over-represented codons (RSCU value > 1.6) in SARS-CoV-2 genome as shown in Table 3. The preferred codons in SARS-CoV-2 genome are GCA (Ala), CGU (Arg), AAU (Asn), GAU (Asp), UGU (Cys), CAA (Gln), GAA (Glu), CAU (His), AUU (Ile), UUG (Leu), AAA (Lys), UUU (Phe), CCA (Pro), AGU (Ser) and UAU (Tyr) while the over-represented codons are GCU (Ala), AGA (Arg), GGU (Gly), CUU (Leu), UUA (Leu), CCU (Pro), UCA (Ser), UCU (Ser), ACA (Thr), ACU (Thr), and GUU (Val). The presence of the preferred and over-presented codons in a genome increases the protein synthesis rate.

Table 4 shows the RSCU analysis of the top 10 synonymous mutations. The codons in bold in the ‘codon change’ column are the codons with higher RSCU value, which means they are more preferred in SARS-CoV-2 genome. Most of the mutations change the codon to a more preferred codon as shown in Table 4. Nine of the ten synonymous mutations involve changes from C to U nucleotides and eight of them are located at the third position of codon, suggesting these changes are not random and possibly subjected to some selection pressure. In agreement with our study, the excessive changes of C to U nucleotides in SARS-CoV-2 genome has been reported in multiple studies 32– 35 . Since the preferred codons may have a better translation efficiency and accuracy compared to the nonpreferred codons 64 , it is possible that most of these mutations may increase the viral fitness. While a study show that RNA secondary structures may be functionally linked to protein translation based on the evidence obtained from experimental work 67 , it is difficult for us to establish the connection solely using in silico studies.

Table 4. RSCU analysis of the top 10 synonymous mutations of SARS-CoV-2 genome.

ORF	Mutation	Codon Change	
1a	nsp1	C313U	CUC -> CUU	
nsp2	C913U	UCC -> UCU	
nsp3	C3037U	UUC -> UUU	
C5986U	UUC -> UUU	
nsp4	C9286U	AAC -> AAU	
1b	nsp12	C14676U	CCC -> CCU	
C15279U	CAC -> CAU	
U16176C	ACU -> ACC	
nsp14	C18877U	CUA -> UUA	
M	C26735U	UAC -> UAU	

Conclusions

The effects of SARS-CoV-2 synonymous mutations in various aspects such as RNA secondary structure and codon usage bias were studied, even though they do not cause changes in amino acid residue of the protein. C913U, C3037U, U16176C and C18877U mutants show pronounced changes between wild type and mutant predicted in all 3 RNA secondary structure prediction tools, suggesting these mutations may have some biological impact on viral fitness. In addition, these mutations showed changes in base pair potential estimated by MutaRNA. All mutations except U16176C change the codon to a more preferred codon, which may result in higher translation efficiency. Due to the shortcomings of prediction tools, experimental studies, such as protein translation assays, RNA packaging assays, are needed to give a more comprehensive understanding of the biological consequences of synonymous mutations on SARS-CoV-2 virus.

Ethics and dissemination

No ethical approval is required for data analysis in this study (EA2702021).

Data and software availability

Underlying data

SARS-CoV-2 virus genome sequence data were obtained from the GISAID Database. The multiple alignment data can be assessed through FigShare.

Figshare: MSA (SARS-CoV-2). https://doi.org/10.6084/m9.figshare.20486178.v1 68

Extended data

Figshare: RNA secondary structure prediction and base pair probability estimation analysis

https://doi.org/10.6084/m9.figshare.20486166.v6 52

Extended data 1. Comparation of RNA secondary structure of SARS-CoV-2 5’ UTR (1-480 nt) predicted using RNAfold without SHAPE data, RNAfold with SHAPE data, IPknot++, MXfold2

Extended data 2. The RNA secondary structure of SARS-CoV-2 genome predicted using RNAfold.

Extended data 3. The RNA secondary structure of SARS-CoV-2 genome predicted using IPknot++.

Extended data 4. The RNA secondary structure of SARS-CoV-2 genome predicted using MXfold2.

Extended data 5. The base pair probabilities of SARS-CoV-2 genome estimated using MutaRNA

Data are available under the terms of the Creative Commons Attribution 4.0 International license (CC-BY 4.0).

Software

The python code for the identification of SARS-CoV-2 genome mutations can be assessed through GitHub.

Author contributions

CHN contributes to the concept, design, supervision of the project. WXB and SBZ contribute to the design, methodology, and data collection. WXB contributed to the analysis, and interpretation of data. All authors were involved in drafting and revising the manuscript and approved the final version.

10.5256/f1000research.162881.r304727
Reviewer response for version 3
Schlick Tamar 1Referee https://orcid.org/0000-0002-2392-2062

1 Department of Chemistry and Courant Institute of Mathematical Sciences, New York University, New York, New York, USA
11 9 2024 Copyright: © 2024 Schlick T
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Version 3recommendationapprove-with-reservations
In this work, synonymous mutations of the SARS-CoV-2 genome are explored, with the rationale that these mutations impact function through altered RNA folding, despite unaltered protein products. Specifically, the researchers find 150 synonymous mutations after performing multiple sequence alignment and examining the protein products, and mapped the distributions in the coding regions. RNA secondary structure predictions and base pair probability calculations are then presented for 4 mutations C913U, C3037U, U16176C and C18877U that show pronounced changes between wildtype and mutant structures. Different prediction tools were used to confirm the structure predictions. The value of the study, while interesting, is limited because the long RNA lengths used make the disparate predicted structures unreliable. If the Information Classification: General authors examine some mutant systems experimentally with chemical reactivity like SHAPE or DMS, the results may be more meaningful.

Is the work clearly and accurately presented and does it cite the current literature?

Yes

If applicable, is the statistical analysis and its interpretation appropriate?

Partly

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Partly

Are the conclusions drawn adequately supported by the results?

Partly

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.

10.5256/f1000research.162881.r304726
Reviewer response for version 3
Forni Diego 1Referee
1 Scientific Institute IRCCS E Medea, Bosisio Parini, Italy
27 7 2024 Copyright: © 2024 Forni D
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Version 3recommendationreject
In this manuscript, Dr Boon et al. analyzed the most common synonymous mutations of  early pandemic SARS-CoV-2 genomes. They identified substitutions that influence the RNA secondary structure surrounding these positions. The authors have already modified and updated their manuscript based on previous comments, still there are few open questions in my opinion.

First of all, it is clear that the dataset is not up to date and it will not be changed, still this should be clearly stated, starting from the title, reporting or the time period or something like "at the beginning of the pandemic''.

- Following this, also the introduction section should be focused on what we know about genomic variability and viral evolution in the time period being analyzed herein. If I understood well, most of these sequences belong to the alpha variant.

- Moreover, " variant of concern" (abbreviated VOC not VOI as stated in the text) refers to a specific lineage that carries several and distinct mutations. As it is stated in the text now it seems that a single mutation is a VOC. Please check it.

- It is expected that most mutations are found in ORF1ab, as this covers most of the genome, so counts should be normalized by ORF length (e.g. fig.1).

- Table 1 should report mutation frequency not count. 

- Authors should discuss more why different predictors give very different results in terms of rna secondary structure, and how this can affect their analysis.

- Why did the author select 250 bp flanking positions for their prediction? Which is the rationale? How does this influence the results?

-  Are there any other mutations in the 501 nucleotide regions? Could mutations have an impact on predictions?

- How did the authors handle overlapping/internal ORFs? In these situations a syn mutations could also be non-synonymous  for the other ORF. I do not think this is the case for the top 10 mutations, but still this should be explained.

- The CUB section seems to me not really linked to the rest of the analyses and could be expanded more.

Is the work clearly and accurately presented and does it cite the current literature?

Partly

If applicable, is the statistical analysis and its interpretation appropriate?

Yes

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Partly

Are the conclusions drawn adequately supported by the results?

Partly

Are sufficient details of methods and analysis provided to allow replication by others?

Partly

Reviewer Expertise:

viral genomics, viral evolution

I confirm that I have read this submission and believe that I have an appropriate level of expertise to state that I do not consider it to be of an acceptable scientific standard, for reasons outlined above.

Ng Chong Han Multimedia University, Malaysia

5 9 2024 In this manuscript, Dr Boon et al. analyzed the most common synonymous mutations of early pandemic SARS-CoV-2 genomes. They identified substitutions that influence the RNA secondary structure surrounding these positions. The authors have already modified and updated their manuscript based on previous comments, still there are few open questions in my opinion.

First of all, it is clear that the dataset is not up to date and it will not be changed, still this should be clearly stated, starting from the title, reporting or the time period or something like "at the beginning of the pandemic''.

Author response:  The title has been revised to “Prediction of the effects of the top 10 synonymous mutations from 26645 SARS-CoV-2 genomes of early pandemic phase”.

Reviewer comment: Following this, also the introduction section should be focused on what we know about genomic variability and viral evolution in the time period being analyzed herein. If I understood well, most of these sequences belong to the alpha variant.

Author response: We added a few lines about Alpha variant in introduction, and result and discussion section (Identification of SARS-CoV-2 synonymous mutations).

Reviewer comment: Moreover, " variant of concern" (abbreviated VOC not VOI as stated in the text) refers to a specific lineage that carries several and distinct mutations. As it is stated in the text now it seems that a single mutation is a VOC. Please check it.

Author response: The abbreviation has been corrected from VOI to VOC. In the text, we mentioned that different SARS-CoV-2 variants with multiple synonymous and nonsynonymous mutations have been reported since the beginning of the outbreak. Therefore, it is clear that VOCs contain multiple mutations.

Reviewer comment: It is expected that most mutations are found in ORF1ab, as this covers most of the genome, so counts should be normalized by ORF length (e.g. fig.1).

Author response: Figure 1 show the distribution of 150 different types of synonymous mutations identified from 26645 SARS-CoV-2 genomes. It is not about the count or the average of the mutations in each ORFs, thus it should not be normalized.

Reviewer comment: Table 1 should report mutation frequency not count.

Author response: Table 1 shows the mutation frequency as reported in this paper.

Reviewer comment: Authors should discuss more why different predictors give very different results in terms of rna secondary structure, and how this can affect their analysis.

Author response: In the paper, there is about two paragraphs on comparison of these RNA secondary structure tools and we highlighted their similarities and differences of these tools. As mentioned in the discussion, we expected some differences in the outcomes obtained from different RNA secondary structure prediction tools due to the differences in algorithms, thermodynamic parameter settings, inclusion of pseudoknot calculation, incorporation of experimental data and the assumptions of each tool. Furthermore, SARS-CoV-2 may adopt different RNA secondary structure conformations. However, our major focus is to predict if the sSNP may affect/change RNA secondary structure, rather than the predicted structures per se. Importantly the outcomes allow us to prioritize variants for the experiment functional studies in the future. In addition, we demonstrated the usability of prediction tools by using 5’ UTR (1-480 nt) of SARS-CoV-2, which is well-characterized in other studies in the Extended data 1. For the four sSNPs which showed pronounced changes, we also compared the differences of these structures predicted by different tools.

Reviewer comment: Why did the author select 250 bp flanking positions for their prediction? Which is the rationale? How does this influence the results?

Author response: Some sSNPs affect only RNA secondary structures of neighboring nucleotides while some sSNPs may have long-ranging effects, in addition to the changes in the neighboring nucleotides. The length of RNA sequences used should not be too short since we may miss out the long-ranging effects. On the other hand, we need to bear in mind that most of RNA secondary structure prediction tools are not suitable for long sequences analysis. We tried out different lengths of flanking sequences, including 100b, 250b, 500b and 1000b, and we found out that 250b flanking sequences produce the optimal results.

Reviewer comment:  Are there any other mutations in the 501 nucleotide regions? Could mutations have an impact on predictions?

Author response: Yes, there may be other synonymous mutations within the 501 nucleotide regions, but they are not from the top 10 synonymous mutations. We are not doing epitasis analysis of different mutations since it is out of the scope of this study.

Reviewer comment: How did the authors handle overlapping/internal ORFs? In these situations a syn mutations could also be non-synonymous for the other ORF. I do not think this is the case for the top 10 mutations, but still this should be explained.

Author response: Yes, we don’t find any synonymous mutation of one ORF, which is also a synonymous mutation for another ORF. We included one sentence about the overlapping ORFs in the result and discussion section in Identification of SARS-CoV-2 synonymous mutation. Since some mutations, including C913U, U16176C are located near the boundaries of different ORFs, we postulated their potential long ranging effects at the neighboring ORFs in the discussion section.

Reviewer comment: The CUB section seems to me not really linked to the rest of the analyses and could be expanded more.

Author response: We expanded CUB section, including additional relevant references, elaborating the results further to make it clearer to the reader.

10.5256/f1000research.137703.r208617
Reviewer response for version 2
Huber Roland 1Referee https://orcid.org/0000-0001-5093-5988

1 Bioinformatics Institute, A*STAR, Singapore
18 10 2023 Copyright: © 2023 Huber R
2023
https://creativecommons.org/licenses/by/4.0/ This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Version 2recommendationreject
I agree with previous reviewers that the analyzed sequences represent only a limited sample of variation in SARS-CoV-2, specifically from early in the pandemic. This might introduce unexpected biases in the analysis. E.g. it would be more likely to observe host adaption early on which would be consistent with more favourable codon usage.

With regard to the data that was analysed, the structure models obtained show limited consistency. The authors state that they do not expect the used tools to concur on the structures since they use different algorithms. This is concerning, as one would expect unambiguous structures to be consistent, even using different methodologies. Other tools, e.g. RNAstructure, also allow the inclusion of shape data and the prediction of pseudoknots. We are thus left with a series of diverging structure predictions and unsure what, if any, effect these specific mutations have. This is not helped by inconsistent presentation of the results. Figures 2-5 use different visualisations for the results of the 4 tools employed, which makes comparisons of the structures difficult for the reader.

The study unfortunately does nothing to associate the regions or structural elements with any type of functional or biological information. We are thus left with a study that analyses a limited set of synonymous mutations using inconsistent structure predictions and offers no additional biological insight.

Is the work clearly and accurately presented and does it cite the current literature?

Partly

If applicable, is the statistical analysis and its interpretation appropriate?

Partly

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Partly

Are the conclusions drawn adequately supported by the results?

Partly

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

Computational Biology, Structural Genomics, RNA biology, Virology

I confirm that I have read this submission and believe that I have an appropriate level of expertise to state that I do not consider it to be of an acceptable scientific standard, for reasons outlined above.

10.5256/f1000research.137703.r149467
Reviewer response for version 2
Nithin Chandran 1Referee https://orcid.org/0000-0001-8212-6093

1 Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, Warsaw, Poland
6 9 2022 Copyright: © 2022 Nithin C
2022
https://creativecommons.org/licenses/by/4.0/ This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Version 2recommendationapprove
Authors have addressed the comments. However, the data set used in the study is outdated.

Is the work clearly and accurately presented and does it cite the current literature?

No

If applicable, is the statistical analysis and its interpretation appropriate?

Partly

Are all the source data underlying the results available to ensure full reproducibility?

Partly

Is the study design appropriate and is the work technically sound?

Partly

Are the conclusions drawn adequately supported by the results?

No

Are sufficient details of methods and analysis provided to allow replication by others?

Partly

Reviewer Expertise:

RNA structure prediction; Computational Structural Biology; RNA-protein complexes.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

10.5256/f1000research.76505.r135257
Reviewer response for version 1
Nithin Chandran 1Referee https://orcid.org/0000-0001-8212-6093

1 Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, Warsaw, Poland
28 4 2022 Copyright: © 2022 Nithin C
2022
https://creativecommons.org/licenses/by/4.0/ This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Version 1recommendationreject
The manuscript titled "Prediction of the Effects of Synonymous Variants on SARS-CoV-2 Genome" analyzes 30,299 SARS-CoV-2 genomes retrieved from GISAID for the period between 31st December 2019 and 22nd March 2021. The study identifies 381 mutations in the genome, of which 150 are synonymous. The study discusses two mutations, C913U and C26735U. The dataset used in this study is outdated and represents only a short period between 31st December 2019 and 22nd March 2021. It will be necessary to update the data to include more sequences from a wider time frame.

The authors had indicated (through the editorial manager) that they are unwilling to update the dataset; the authors need to indicate the same in the title and detail the limitations of the dataset in the introduction section.

The authors have used the RNAfold program to predict the secondary structure of wild-type and mutant sequences. The authors need to use additional prediction tools and check for consensus between the predicted structure. Moreover, authors need to employ prediction tools to predict pseudoknots in the RNA secondary structure.

The secondary structure of SARS-CoV2 determined with the help of data from SHAPE-MaP, icSHAPE, and DMS-MaPseq experiments is available in the literature. 1 , 2 , 3 The authors should utilize this data to augment this study.

The study requires the inclusion of appropriate controls which quantify the effect of mutations on the stability of the RNA structure. The authors state that "The minimum free energy value of the mutant (- 146.90 kcal/mol) is slightly less negative than that of the wild type (- 147.40 kcal/mol), which makes it a less thermodynamically stable structure compared to the wild type." This claim needs to be explored in detail and supported with evidence from experimental data and literature.

The authors state that C26735U mutation induces changes in the predicted RNA secondary structure by forming an extra multibranch loop at the mutation site. The authors should compare the same with secondary systems mentioned in point 3 (See Ref 1, 2, and 3). The authors need to perform a similar analysis for C913U.

The study identifies 381 mutations in the genome, of which 150 are synonymous; however, it discusses only two. The authors need to provide the data on the remaining 148 synonymous mutations and briefly discuss some of them in the manuscript.

The authors need to expand and include relevant literature in the introduction section.

The manuscript needs to be edited for language errors.

Is the work clearly and accurately presented and does it cite the current literature?

No

If applicable, is the statistical analysis and its interpretation appropriate?

Partly

Are all the source data underlying the results available to ensure full reproducibility?

Partly

Is the study design appropriate and is the work technically sound?

Partly

Are the conclusions drawn adequately supported by the results?

No

Are sufficient details of methods and analysis provided to allow replication by others?

Partly

Reviewer Expertise:

RNA structure prediction; Computational Structural Biology; RNA-protein complexes.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to state that I do not consider it to be of an acceptable scientific standard, for reasons outlined above.

Ng Chong Han Multimedia University, Malaysia

22 8 2022 The manuscript titled "Prediction of the Effects of Synonymous Variants on SARS-CoV-2 Genome" analyzes 30,299 SARS-CoV-2 genomes retrieved from GISAID for the period between 31st December 2019 and 22nd March 2021. The study identifies 381 mutations in the genome, of which 150 are synonymous. The study discusses two mutations, C913U and C26735U. The dataset used in this study is outdated and represents only a short period between 31st December 2019 and 22nd March 2021. It will be necessary to update the data to include more sequences from a wider time frame.

The authors had indicated (through the editorial manager) that they are unwilling to update the dataset; the authors need to indicate the same in the title and detail the limitations of the dataset in the introduction section.

Response: The landscape of mutational profile of SARS-CoV-2 genome is very dynamic, changing rapidly. As of 14 August 2022, there are more than 12 million SARS-CoV-2 virus genomes sequences deposited in GISAID database. Although it is useful to get more recent dataset, it is out of the scope of our study to update the data. Both identification of the SARS-CoV-2 mutations and the prediction analysis of the biological consequences of these mutations are very time-consuming processes. Based on the timelines of emergence of different variants, the cut-off date of our data collection overlaps with the period when the frequency of the alpha variant peaked. All the synonymous mutations except C241T reported in the alpha variant have been identified in our study as well. Although the alpha variant has been replaced by other variants in human population, it is detected from wildlife animals based on a few recent genomic surveillance analysis. It is possible that these alpha variant from wildlife animal may jump back to human, hence it remains relevant to predict the biological consequences of these synonymous mutations. The title of the paper has been revised accordingly. The limitation of the dataset was included in the discussion.

The authors have used the RNAfold program to predict the secondary structure of wild-type and mutant sequences. The authors need to use additional prediction tools and check for consensus between the predicted structure. Moreover, authors need to employ prediction tools to predict pseudoknots in the RNA secondary structure.

Response: To improve the outcome of the RNA secondary structure analysis, 3 prediction tools, namely RNAfold, IPknot++ and MXfold2 are used in our study. Of these, IPknot++ predicted some synonymous mutations may affect the pseudoknot formation. The prediction results for all top 10 synonymous mutations using these 3 tools are summarized in the Table 2

The secondary structure of SARS-CoV2 determined with the help of data from SHAPE-MaP, icSHAPE, and DMS-MaPseq experiments is available in the literature.1,2,3 The authors should utilize this data to augment this study

Response: To improve the outcome of the prediction analysis, the RNA secondary structure of wild type and mutant sequences were predicted using RNAfold program with the incorporation of SHAPE reactivity data obtained from the study done by Manfredonia et al. (2020).

The study requires the inclusion of appropriate controls which quantify the effect of mutations on the stability of the RNA structure. The authors state that "The minimum free energy value of the mutant (- 146.90 kcal/mol) is slightly less negative than that of the wild type (- 147.40 kcal/mol), which makes it a less thermodynamically stable structure compared to the wild type." This claim needs to be explored in detail and supported with evidence from experimental data and literature.

Response: Instead of comparing minimum free energy between the wild type and mutant, we applied three RNA secondary structure prediction tools to predict the effect of the mutation on RNA secondary structure.

The authors state that C26735U mutation induces changes in the predicted RNA secondary structure by forming an extra multibranch loop at the mutation site. The authors should compare the same with secondary systems mentioned in point 3 (See Ref 1, 2, and 3). The authors need to perform a similar analysis for C913U.

Response: The results of the 4 mutations showing changes in all 3 prediction tools are included in Figure 2-5.  

The study identifies 381 mutations in the genome, of which 150 are synonymous; however, it discusses only two. The authors need to provide the data on the remaining 148 synonymous mutations and briefly discuss some of them in the manuscript.

Response: The RNA secondary structures prediction results using 3 prediction tools for the top 10 synonymous mutations have been included in Extended data 1-3. There were 4 mutations showing changes in all 3 prediction tools and the results are discussed in the manuscript.

The authors need to expand and include relevant literature in the introduction section.

Response: The relevant literature in the introduction section has been updated.

The manuscript needs to be edited for language errors.

Response: The language errors have been edited accordingly.

10.5256/f1000research.76505.r101828
Reviewer response for version 1
Wang Leyi 1Referee https://orcid.org/0000-0001-5813-9505

1 Veterinary Diagnostic Laboratory and Department of Veterinary Clinical Medicine, College of Veterinary Medicine, University of Illinois, Urbana, IL, USA
20 12 2021 Copyright: © 2021 Wang L
2021
https://creativecommons.org/licenses/by/4.0/ This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Version 1recommendationapprove-with-reservations
The manuscript entitled “Prediction of the Effects of Synonymous Variants on SARS-CoV-2 Genome” reports an analysis of the SARS-CoV-2 genome obtained in the GISAID database. The authors analyzed synonymous changes of the 30,229 SARS-CoV-2 sequences available online, and the effects of these changes on RNA secondary structures and base-pair probability. They identified 150 synonymous mutations in 11 coding regions. Mutations with higher frequency are mostly found in ORF1a and ORF1b and two mutations (C913U and C26735U) found with effect on the predicted RNA secondary structure are not those with higher frequency.

This study collected 30,229 sequences in the database in the period of 31 Dec 2019 to 22 March 2021. As of 17 Dec 2021, there are over 4.4 million complete genomes with high coverage. Will these sequences used in the study include all variants reported or what sequences of variants are included? A table to summarize all variants with these mutations or not is needed.

A major concern is that the sequences even filtered by setting parameters to keep only sequencing with complete genome and high coverage still contain bad sequences, how did the authors exclude those sequences?

In the abstract section, the results part only included those mutations with high frequency but did not contain results of other analyses. Instead, they included the results of RNA secondary structure analysis in the conclusion.

Is the work clearly and accurately presented and does it cite the current literature?

Partly

If applicable, is the statistical analysis and its interpretation appropriate?

Yes

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Partly

Are the conclusions drawn adequately supported by the results?

Partly

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

Clinical Virology Diagnosis

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.

Ng Chong Han Multimedia University, Malaysia

22 8 2022 The manuscript entitled “Prediction of the Effects of Synonymous Variants on SARS-CoV-2 Genome” reports an analysis of the SARS-CoV-2 genome obtained in the GISAID database. The authors analyzed synonymous changes of the 30,229 SARS-CoV-2 sequences available online, and the effects of these changes on RNA secondary structures and base-pair probability. They identified 150 synonymous mutations in 11 coding regions. Mutations with higher frequency are mostly found in ORF1a and ORF1b and two mutations (C913U and C26735U) found with effect on the predicted RNA secondary structure are not those with higher frequency.

This study collected 30,229 sequences in the database in the period of 31 Dec 2019 to 22 March 2021. As of 17 Dec 2021, there are over 4.4 million complete genomes with high coverage. Will these sequences used in the study include all variants reported or what sequences of variants are included? A table to summarize all variants with these mutations or not is needed.

Response: We downloaded the raw data on 23 rd March 2021 and completed the data analysis in late June 2021. We submitted the manuscript in late August 2021. The landscape of mutational profile of SARS-CoV-2 genome is very dynamic, changing rapidly. It is far beyond the scope of this study to include all variants reported. Based on the timelines of emergence of different variants, the cut-off date of our data collection overlaps with the period when the frequency of the alpha variant peaked. All the synonymous mutations except C241T reported in the alpha strain have been identified in our study as well.  We added new paragraph in the discussion section to discuss the impact of our findings and to explain why the study of alpha variant remains relevant now. In addition, we used 3 RNA secondary prediction analysis tool, instead of one, to improve the outcome of the analysis. The prediction results for all top 10 synonymous mutations using these 3 tools and the base pair probability estimation results are summarized in the Table 2.

A major concern is that the sequences even filtered by setting parameters to keep only sequencing with complete genome and high coverage still contain bad sequences, how did the authors exclude those sequences?

Response: In our original method, we used the high coverage filter when downloading from GISAID database, in which only entries with less 1% N and 0.05% unique amino acids mutations are included. To further reduce bad sequences, we filtered to remove those sequences with higher than 0.1% N unresolved nucleotides and ambiguous letters. A total of 3584 sequences were removed by applying this filter. The list of the synonymous mutations remains unchanged despite of the revision of the mutation frequency.

In the abstract section, the results part only included those mutations with high frequency but did not contain results of other analyses. Instead, they included the results of RNA secondary structure analysis in the conclusion.

Response: The abstract has been revised to update the results and conclusion.

10.5256/f1000research.76505.r97299
Reviewer response for version 1
Koyama Takahiko 1Referee https://orcid.org/0000-0003-1694-9061

1 IBM TJ Watson Research Center, Yorktown Heights, NY, USA
2 11 2021 Copyright: © 2021 Koyama T
2021
https://creativecommons.org/licenses/by/4.0/ This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Version 1recommendationreject
Boon et al. made an investigation of synonymous mutations of SARS-CoV-2 genomes. Authors first identified common synonymous mutations such as 3037C>T followed by secondary structure predictions using RNAfold. 

First, the number of genomes authors used is too small and too old. Currently, nearly 5 million genomes are uploaded at GISAID; however, authors used only 30,229. The set does not represent the current mutational landscape of SARS-CoV-2. Authors need to update the data. I have 3 million genomes analyzed and 18108 distinct synonymous mutations identified. The list of highest frequency in Table 1 looks completely different with the latest dataset.

Secondly, authors employed RNAfold for RNA secondary structure prediction. Although RNAfold is a standard software for the task, the accuracy of the prediction is not impressive. Authors need to explore various other tools for the secondary structure or justify the usage of RNAfold.

Is the work clearly and accurately presented and does it cite the current literature?

Partly

If applicable, is the statistical analysis and its interpretation appropriate?

Yes

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Partly

Are the conclusions drawn adequately supported by the results?

Partly

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

Genome Analysis, SARS-CoV-2, Cancer, Immunology, Stem cell

I confirm that I have read this submission and believe that I have an appropriate level of expertise to state that I do not consider it to be of an acceptable scientific standard, for reasons outlined above.

Ng Chong Han Multimedia University, Malaysia

4 11 2021 Boon et al. made an investigation of synonymous mutations of SARS-CoV-2 genomes. Authors first identified common synonymous mutations such as 3037C>T followed by secondary structure predictions using RNAfold.

First, the number of genomes authors used is too small and too old. Currently, nearly 5 million genomes are uploaded at GISAID; however, authors used only 30,229. The set does not represent the current mutational landscape of SARS-CoV-2. Authors need to update the data. I have 3 million genomes analyzed and 18108 distinct synonymous mutations identified. The list of highest frequency in Table 1 looks completely different with the latest dataset.

RE: Thank you for the comments and feedbacks. We downloaded SARS-CoV-2 virus genome data submitted from 31 December 19 to 22 March 21 on 23 March 21. There were about a total of 842,603 SARS-CoV-2 virus genomes sequences before applying filters. SARS-CoV-2 genomic sequences were filtered by setting parameters to keep only sequences with complete genome, high coverage, low coverage excluded and patient status options. The use of the low coverage excluded and patient status filters were unintentionally left out in the manuscript and they will be added later. In our study, the use of these filters allows us to analyze high quality data to identify SARS-CoV-2 mutations without the requirement of high computing power. Our data is representative of the mutational landscape of SARS-CoV-2 for this specified period. Same synonymous mutations, including C913T, C3037T, C5986T, C14676T, C15279T and T16176C were identified in alpha strain by other research groups ( https://covariants.org/variants/20I.Alpha.V1) and ( https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/959438/Technical_Briefing_VOC_SH_NJL2_SH2.pdf). The frequency of alpha strain peaked in April-May 21 and subsequently it has been quickly replaced by delta strain. Interestingly the defining synonymous mutations of delta strain ( https://covariants.org/variants/21A.Delta) are very different from those of alpha strain. The mutational landscape of SARS-CoV-2 genome is very dynamic as shown by other research groups. Bear in mind that the main objective of our paper is not to monitor the mutational changes of SARS-CoV-2 genome. Although it is useful to get more recent dataset, it is out of the scope of our study to update the data. Both identification of the SARS-CoV-2 mutations and the prediction analysis of the biological consequences of these mutations are very time-consuming processes. The first draft of the manuscript was prepared in late June 21 and submitted in late August 21. Since different SARS-CoV-2 variants may have different sets of mutations, we will gain more insight to study the effects of different mutations in different SARS-CoV-2 variants. In our follow-up study, we plan to predict and compare the effect of the mutations of different variants. These analyses may provide some plausible explanation why some variants are replaced by other variants eventually.

Secondly, authors employed RNAfold for RNA secondary structure prediction. Although RNAfold is a standard software for the task, the accuracy of the prediction is not impressive. Authors need to explore various other tools for the secondary structure or justify the usage of RNAfold.

RE: We chose RNAfold because it is the most commonly used tools to predict the RNA secondary structure. Other than RNAfold, we also considered using Mfold to perform the RNA secondary structure prediction as well. However, the Mfold server was down during the time of the study. RNAfold computes an optimal structure over the whole sequence length, which means its performance and accuracy can be affected as the sequences used get longer. However, in our case, the accuracy of the prediction is still acceptable since the sequence we used to do the structure prediction is relatively short. Besides, the prediction results we obtained from RNAfold is further supported with the prediction results obtained from MutaRNA, which predicts the structural changes induced by the mutation by estimating the base pairing probabilities. In our results, the circular plot from MutaRNA shows the changes in the base pairing probabilities near the mutation site correlates well with the RNA secondary structure predicted by RNAfold.

Ng Chong Han Multimedia University, Malaysia

21 11 2022 First, the number of genomes authors used is too small and too old. Currently, nearly 5 million genomes are uploaded at GISAID; however, authors used only 30,229. The set does not represent the current mutational landscape of SARS-CoV-2. Authors need to update the data. I have 3 million genomes analyzed and 18108 distinct synonymous mutations identified. The list of highest frequency in Table 1 looks completely different with the latest dataset.

RE: We submitted a revised version in Sep 2022. Although we didn’t use the latest dataset, we provide some justification in the discussion to explain why our data may remain relevant. We hope you will consider to review our manuscript again. The landscape of mutational profile of SARS-CoV-2 genome is very dynamic, changing rapidly. It is far beyond the scope of this study to include all variants reported. Based on the timelines of emergence of different variants, the cut-off date of our data collection overlaps with the period when the frequency of the alpha variant peaked. All the synonymous mutations except C241T reported in the alpha strain have been identified in our study as well. Although the alpha variant has been replaced by other variants in human population, it is detected from wildlife animals based on a few recent genomic surveillance analysis. It is possible that these alpha variant from wildlife animal may jump back to human, hence it remains relevant to predict the biological consequences of these synonymous mutations. The title of the paper has been revised accordingly. The limitation of the dataset was included in the discussion.

Secondly, authors employed RNAfold for RNA secondary structure prediction. Although RNAfold is a standard software for the task, the accuracy of the prediction is not impressive. Authors need to explore various other tools for the secondary structure or justify the usage of RNAfold.

RE: In the revised version, we used 3 RNA secondary prediction analysis tool, namely RNAfold, IPknot++, to improve the outcome of the analysis. The prediction results for all top 10 synonymous mutations using these 3 tools and the base pair probability estimation results are summarized in the Table 2.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.

Competing interests: No competing interests were disclosed.
==== Refs
1 Huang C Wang Y Li X : Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395 (10223 ):497–506. 10.1016/S0140-6736(20)30183-5 31986264
2 Wu D Wu T Liu Q : The SARS-CoV-2 outbreak: what we know. Int J Infect Dis. 2020;94 :44–48. 10.1016/j.ijid.2020.03.004 32171952
3 Sharma A Tiwari S Deb MK : Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2): a global pandemic and treatment strategies. Int J Antimicrob Agents. 2020;56 (2 ): 106054. 10.1016/j.ijantimicag.2020.106054 32534188
4 Sanjuán R Domingo-Calap P : Mechanisms of viral mutation. Cell Mol Life Sci. Birkhauser Verlag AG,2016;73 (23 ): 4433–4448. 10.1007/s00018-016-2299-6 27392606
5 CoVariants. Reference Source
6 Tao K Tzou PL Nouhin J : The biological and clinical significance of emerging SARS-CoV-2 variants. Nat Rev Genet. 2021;22 (12 ):757–773. 10.1038/s41576-021-00408-x 34535792
7 Carabelli AM Peacock TP Thorne LG : SARS-CoV-2 variant biology: immune escape, transmission and fitness. Nat Rev Microbiol. 2023;21 (3 ):162–177. 10.1038/s41579-022-00841-7 36653446
8 Ozono S Zhang Y Ode H : SARS-CoV-2 D614G spike mutation increases entry efficiency with enhanced ACE2-binding affinity. Nat Commun. 2021;12 (1 ): 848. 10.1038/s41467-021-21118-2 33558493
9 Plante JA Liu Y Liu J : Spike mutation D614G alters SARS-CoV-2 fitness. Nature. 2021;592 (7852 ):116–121. 10.1038/s41586-020-2895-3 33106671
10 Zhu L Wang Q Zhang W : Evidence for selection on SARS-CoV-2 RNA translation revealed by the evolutionary dynamics of mutations in UTRs and CDSs. RNA Biol. 2022;19 (1 ):866–876. 10.1080/15476286.2022.2092351 35762570
11 de Maio N Walker CR Turakhia Y : Mutation rates and selection on synonymous mutations in SARS-CoV-2. Genome Biol Evol. 2021;13 (5 ): evab087. 10.1093/gbe/evab087 33895815
12 Bin Y Wang X Zhao L : An analysis of mutational signatures of synonymous mutations across 15 cancer types. BMC Med Genet. 2019;20 (Suppl 2 ): 190. 10.1186/s12881-019-0926-4 31815613
13 Sharp PM Averof M Lloyd AT : DNA sequence evolution: the sounds of silence. Philos Trans R Soc Lond B Biol Sci. 1995;349 (1329 ):241–247. 10.1098/rstb.1995.0108 8577834
14 Chamary JV Parmley JL Hurst LD : Hearing silence: Non-neutral evolution at synonymous sites in mammals. Nat Rev Genet. 2006;7 (2 ):98–108. 10.1038/nrg1770 16418745
15 Shen X Song S Li C : Synonymous mutations in representative yeast genes are mostly strongly non-neutral. Nature. 2022;606 (7915 ):725–731. 10.1038/s41586-022-04823-w 35676473
16 Burrill CP Westesson O Schulte MB : Global RNA structure analysis of poliovirus Identifies a conserved RNA structure involved in viral replication and infectivity. J Virol. 2013;87 (21 ):11670–11683. 10.1128/JVI.01560-13 23966409
17 Mueller S Papamichail D Coleman JR : Reduction of the rate of poliovirus protein synthesis through large-scale codon deoptimization causes attenuation of viral virulence by lowering specific infectivity. J Virol. 2006;80 (19 ):9687–9696. 10.1128/JVI.00738-06 16973573
18 Lauring AS Acevedo A Cooper SB : Codon usage determines the mutational robustness, evolutionary capacity, and virulence of an RNA virus. Cell Host Microbe. 2012;12 (5 ):623–632. 10.1016/j.chom.2012.10.008 23159052
19 Zhang H Zhang L Lin A : Algorithm for optimized mRNA design improves stability and immunogenicity. Nature. 2023;621 (7978 ):396–403. 10.1038/s41586-023-06127-z 37130545
20 GISAID Initiative. [Accessed: 23-Sep-2021]. Reference Source
21 Wu F Zhao S Yu B : A new coronavirus associated with human respiratory disease in China. Nature. 2020;579 (7798 ):265–269. 10.1038/s41586-020-2008-3 32015508
22 MAFFT - a multiple sequence alignment program. [Accessed: 23-Sep-2021]. Reference Source
23 Kumar S Stecher G Li M : MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol. 2018;35 (6 ):1547–1549. 10.1093/molbev/msy096 29722887
24 Lorenz R Luntzer D Hofacker IL : SHAPE directed RNA folding. Bioinformatics. 2016;32 (1 ):145–7. 10.1093/bioinformatics/btv523 26353838
25 Manfredonia I Nithin C Ponce-Salvatierra A : Genome-wide mapping of SARS-CoV-2 RNA structures identifies therapeutically-relevant elements. Nucleic Acids Res. 2020;48 (22 ):12436–12452. 10.1093/nar/gkaa1053 33166999
26 Sato K Kato Y : Prediction of RNA secondary structure including pseudoknots for long sequences. Brief Bioinform. 2022;23 (1 ): bbab395. 10.1093/bib/bbab395 34601552
27 Sato K Akiyama M Sakakibara Y : RNA secondary structure prediction using deep learning with thermodynamic integration. Nat Commun. 2021;12 (1 ): 941. 10.1038/s41467-021-21194-4 33574226
28 Miladi M Raden M Diederichs S : MutaRNA: analysis and visualization of mutation-induced changes in RNA structure. Nucleic Acids Res. 2020;48 (W1 ):W287–W291. 10.1093/nar/gkaa331 32392303
29 Sharma Y Miladi M Dukare S : A pan-cancer analysis of synonymous mutations. Nat Commun. 2019;10 (1 ): 2569. 10.1038/s41467-019-10489-2 31189880
30 Mochizuki T Ohara R Roossinck MJ : Large-Scale Synonymous Substitutions in Cucumber Mosaic Virus RNA 3 Facilitate Amino Acid Mutations in the Coat Protein. J Virol. 2018;92 (22 ):e01007–18. 10.1128/JVI.01007-18 30185595
31 Sauna ZE Kimchi-Sarfaty C : Synonymous Mutations as a Cause of Human Genetic Disease. In: eLS.Chichester, UK: John Wiley & Sons, Ltd,2013. 10.1002/9780470015902.a0025173
32 Morales AC Rice AM Ho AT : Causes and consequences of purifying selection on SARS-CoV-2. Genome Biol Evol. 2021;13 (10 ): evab196. 10.1093/gbe/evab196 34427640
33 Simmonds P : Rampant C→U hypermutation in the genomes of SARS-CoV-2 and other coronaviruses: causes and consequences for their short- and long-term evolutionary trajectories. mSphere. 2020;5 (3 ):e00408–20. 10.1128/mSphere.00408-20 32581081
34 Dash M Meher P Kumar A : High frequency of transition to transversion ratio in the stem region of RNA secondary structure of untranslated region of SARS-CoV-2. PeerJ. 2024;12 : e16962. 10.7717/peerj.16962 38666080
35 Forni D Cagliani R Pontremoli C : The substitution spectra of coronavirus genomes. Brief Bioinform. 2022;23 (1 ): bbab382. 10.1093/bib/bbab382 34518866
36 Kim K Calabrese P Wang S : The roles of APOBEC-mediated RNA editing in SARS-CoV-2 mutations, replication and fitness. Sci Rep. 2022;12 (1 ): 14972. 10.1038/s41598-022-19067-x 36100631
37 Sun Q Zeng J Tang K : Variation in synonymous evolutionary rates in the SARS-CoV-2 genome. Front Microbiol. 2023;14 : 1136386. 10.3389/fmicb.2023.1136386 36970680
38 Finkel Y Mizrahi O Nachshon A : The coding capacity of SARS-CoV-2. Nature. 2021;589 (7840 ):125–130. 10.1038/s41586-020-2739-1 32906143
39 Lauring AS Hodcroft EB : Genetic variants of SARS-CoV-2-what do they mean? JAMA. 2021;325 (6 ):529–531. 10.1001/jama.2020.27124 33404586
40 Sia BZ Boon WX Yap YY : Prediction of the effects of the top 10 nonsynonymous variants from 30229 SARS-CoV-2 strains on their proteins [version 2; peer review: 2 approved]. F1000Res. 2022;11 :9. 10.12688/f1000research.72904.2 35707000
41 Pickering B Lung O Maguire F : Divergent SARS-CoV-2 variant emerges in white-tailed deer with deer-to-human transmission. Nat Microbiol. 2022;7 (12 ):2011–2024. 10.1038/s41564-022-01268-9 36357713
42 Marques AD Sherrill-Mix S Everett JK : Multiple introductions of SARS-CoV-2 Alpha and Delta variants into white-tailed deer in Pennsylvania. mBio. 2022;13 (5 ): e0210122. 10.1128/mbio.02101-22 36000731
43 Bhatt PR Scaiola A Loughran G : Structural basis of ribosomal frameshifting during translation of the SARS-CoV-2 RNA genome. Science. 2021;372 (6548 ):1306–1313. 10.1126/science.abf3546 34029205
44 Kelly JA Olson AN Neupane K : Structural and functional conservation of the programmed −1 ribosomal frameshift signal of SARS Coronavirus 2 (SARS-CoV-2). J Biol Chem. 2020;295 (31 ):10741–10748. 10.1074/jbc.AC120.013449 32571880
45 Huston NC Wan H Strine MS : Comprehensive in vivo secondary structure of the SARS-CoV-2 genome reveals novel regulatory motifs and mechanisms. Mol Cell. 2021;81 (3 ):584–598. e5. 10.1016/j.molcel.2020.12.041 33444546
46 Lan TCT Allan MF Malsick LE : Secondary structural ensembles of the SARS-CoV-2 RNA genome in infected cells. Nat Commun. 2022;13 (1 ): 1128. 10.1038/s41467-022-28603-2 35236847
47 Cao C Cai Z Xiao X : The architecture of the SARS-CoV-2 RNA genome inside virion. Nat Commun. 2021;12 (1 ): 3917. 10.1038/s41467-021-22785-x 34168138
48 Ziv O Price J Shalamova L : The Short- and Long-Range RNA-RNA Interactome of SARS-CoV-2. Mol Cell. 2020;80 (6 ):1067–1077. e5. 10.1016/j.molcel.2020.11.004 33259809
49 Sato K Hamada M : Recent trends in RNA informatics: a review of machine learning and deep learning for RNA secondary structure prediction and RNA drug discovery. Brief Bioinform. 2023;24 (4 ): bbad186. 10.1093/bib/bbad186 37232359
50 Yang SL DeFalco L Anderson DE : Comprehensive mapping of SARS-CoV-2 interactions in vivo reveals functional virus-host interactions. Nat Commun. 2021;12 (1 ): 5113. 10.1038/s41467-021-25357-1 34433821
51 Boon WX Ng CH : RNA secondary structure prediction and base pair probability estimation analysis.2024. Reference Source
52 Finkel Y Gluck A Nachshon A : SARS-CoV-2 uses a multipronged strategy to impede host protein synthesis. Nature. 2021;594 (7862 ):240–245. 10.1038/s41586-021-03610-3 33979833
53 Vora SM Fontana P Mao T : Targeting stem-loop 1 of the SARS-CoV-2 5' UTR to suppress viral translation and Nsp1 evasion. Proc Natl Acad Sci U S A. 2022;119 (9 ): e2117198119. 10.1073/pnas.2117198119 35149555
54 Yuan J Feng Z Wang Q : 3’UTR of SARS-CoV-2 spike gene hijack host miR-296 or miR-520h to disturb cell proliferation and cytokine signaling. Front Immunol. 2022;13 : 924667. 10.3389/fimmu.2022.924667 36238276
55 Xu Z Choi JH Dai DL : SARS-CoV-2 impairs interferon production via NSP2-induced repression of mRNA translation. Proc Natl Acad Sci U S A. 2022;119 (32 ): e2204539119. 10.1073/pnas.2204539119 35878012
56 Schubert K Karousis ED Jomaa A : SARS-CoV-2 Nsp1 binds the ribosomal mRNA channel to inhibit translation. Nat Struct Mol Biol. 2020;27 (10 ):959–966. 10.1038/s41594-020-0511-8 32908316
57 Burke JM St Clair LA Perera R : SARS-CoV-2 infection triggers widespread host mRNA decay leading to an mRNA export block. RNA. 2021;27 (11 ):1318–1329. 10.1261/rna.078923.121 34315815
58 Wolff G Limpens RWAL Zevenhoven-Dobbe JC : A molecular pore spans the double membrane of the coronavirus replication organelle. Science. 2020;369 (6509 ):1395–1398. 10.1126/science.abd3629 32763915
59 Naydenova K Muir KW Wu LF : Structure of the SARS-CoV-2 RNA-dependent RNA polymerase in the presence of favipiravir-RTP. Proc Natl Acad Sci U S A. 2021;118 (7 ): e2021946118. 10.1073/pnas.2021946118 33526596
60 Terasaki K Narayanan K Makino S : Identification of a 1.4-kb-Long sequence located in the nsp12 and nsp13 coding regions of SARS-CoV-2 genomic RNA that mediates efficient viral RNA packaging. J Virol. 2023;97 (7 ): e0065923. 10.1128/jvi.00659-23 37367225
61 Schlick T Zhu Q Jain S : Structure-altering mutations of the SARS-CoV-2 frameshifting RNA element. Biophys J. 2021;120 (6 ):1040–1053. 10.1016/j.bpj.2020.10.012 33096082
62 Plant EP Sims AC Baric RS : Altering SARS coronavirus frameshift efficiency affects genomic and subgenomic RNA production. Viruses. 2013;5 (1 ):279–94. 10.3390/v5010279 23334702
63 Moeller NH Shi K Demir Ö : Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN. Proc Natl Acad Sci U S A. 2022;119 (9 ): e2106379119. 10.1073/pnas.2106379119 35165203
64 Sun M Zhang J : Preferred synonymous codons are translated more accurately: proteomic evidence, among-species variation, and mechanistic basis. Sci Adv. 2022;8 (27 ): eabl9812. 10.1126/sciadv.abl9812 35857447
65 Wong EH Smith DK Rabadan R : Codon usage bias and the evolution of influenza a viruses. Codon usage biases of influenza virus. BMC Evol Biol. 2010;10 (1 ): 253. 10.1186/1471-2148-10-253 20723216
66 Sharp PM Tuohy TM Mosurski KR : Codon usage in yeast: cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Res. 1986;14 (13 ):5125–5143. 10.1093/nar/14.13.5125 3526280
67 Mauger DM Cabral BJ Presnyak V : mRNA structure regulates protein expression through changes in functional half-life. Proc Natl Acad Sci U S A. 2019;116 (48 ):24075–24083. 10.1073/pnas.1908052116 31712433
68 Boon WX Ng CH : MSA (SARS-CoV-2).(accessed Aug. 17, 2022). Reference Source
