Variation in synonymous evolutionary rates in the SARS-CoV-2 genome

Qianru Sun; Jinfeng Zeng; Kang Tang; Haoyu Long; Chi Zhang; Jie Zhang; Jing Tang; Yuting Xin; Jialu Zheng; Litao Sun; Siyang Liu; Xiangjun Du

doi:10.3389/fmicb.2023.1136386

Variation in synonymous evolutionary rates in the SARS-CoV-2 genome

Front Microbiol. 2023 Mar 9:14:1136386. doi: 10.3389/fmicb.2023.1136386. eCollection 2023.

Authors

Qianru Sun^{1

2}, Jinfeng Zeng^{1

2}, Kang Tang^{1

2}, Haoyu Long^{1

2}, Chi Zhang^{1

2}, Jie Zhang^{1

2}, Jing Tang^{1

2}, Yuting Xin^{1

2}, Jialu Zheng^{1

2}, Litao Sun^{1

2}, Siyang Liu^{1

2}, Xiangjun Du^{1

2

3}

Affiliations

¹ School of Public Health (Shenzhen), Shenzhen Campus of Sun Yat-sen University, Shenzhen, China.
² School of Public Health (Shenzhen), Sun Yat-sen University, Guangzhou, China.
³ Key Laboratory of Tropical Disease Control, Ministry of Education, Sun Yat-sen University, Guangzhou, China.

Abstract

Introduction: Coronavirus disease 2019 is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Influential variants and mutants of this virus continue to emerge, and more effective virus-related information is urgently required for identifying and predicting new mutants. According to earlier reports, synonymous substitutions were considered phenotypically silent; thus, such mutations were frequently ignored in studies of viral mutations because they did not directly cause amino acid changes. However, recent studies have shown that synonymous substitutions are not completely silent, and their patterns and potential functional correlations should thus be delineated for better control of the pandemic.

Methods: In this study, we estimated the synonymous evolutionary rate (SER) across the SARS-CoV-2 genome and used it to infer the relationship between the viral RNA and host protein. We also assessed the patterns of characteristic mutations found in different viral lineages.

Results: We found that the SER varies across the genome and that the variation is primarily influenced by codon-related factors. Moreover, the conserved motifs identified based on the SER were found to be related to host RNA transport and regulation. Importantly, the majority of the existing fixed-characteristic mutations for five important virus lineages (Alpha, Beta, Gamma, Delta, and Omicron) were significantly enriched in partially constrained regions.

Discussion: Taken together, our results provide unique information on the evolutionary and functional dynamics of SARS-CoV-2 based on synonymous mutations and offer potentially useful information for better control of the SARS-CoV-2 pandemic.

Keywords: SARS-CoV-2; binding motif; codon usage; dominant variants; synonymous evolutionary rate.

Grants and funding

This study was supported by the National Key Research and Development Projects of the Ministry of Science and Technology of China under grant 2021YFC2301300; the Shenzhen Science and Technology Program under grant KQTD20180411143323605, JSGG20200225152008136, and GXWD20201231165807008; the Guangdong Frontier and Key Tech Innovation Program under grants 2019B020228001, 2019B111103001, 2021A111112007, and 2022B1111020006; and the Natural Science Foundation of Guangdong Province, China, under grant 2021A1515011592.