Abstract
COVID-19 vaccination resistance has become a major challenge to prevent global SARS-CoV-2 transmission. Here we report that the vaccination coverage rate is inversely correlated to the mutation frequency of the full genome (R2=0.878) and spike gene (R2=0.829) of SARS-CoV-2 delta variants in 16 countries, suggesting that full vaccination against COVID-19, with other mitigation strategies, is critical to suppress emergent mutations. Neutrality analysis of DH and Zeng’s E tests suggested that directional selection was the major driving force of delta variant evolution. To eliminate the homogenous effects (population expansion, selective sweep etc.), the synonymous (Dsyn) and nonsynonymous (Dnonsyn) polymorphisms of the delta variant spike gene were estimated with Tajima’s D statistic. Both D ratio (Dnonsyn/Dsyn) and ΔD (Dsyn-Dnonsyn) have positive correlation with the full vaccination rate (R2= 0.723 and 0.505, respectively) in 19 countries, indicating that purifying selection pressure of SARS-CoV-2 spike gene increased as the vaccination coverage rate increased. Taken together, our data suggests that vaccination plays an important role in the purifying selection force of spike protein of SARS-CoV-2 delta variants.
Competing Interest Statement
The authors have declared no competing interest.
Clinical Trial
Not applicable
Funding Statement
No funding
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
No patients were involved and no IRB was required in this study.
All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Footnotes
This paper is dedicated to Professor Austin Hughes (1949-2015), an evolution biologist at University of South Carolina, for his lifetime achievement in neural theory and adaptive molecular evolution.
The Tajima D forecasting in the previous version was replaced by the synonymous (Dsyn) and nonsynonymous (Dnonsyn) polymorphisms of the delta variant spike gene. Both D ratio (Dnonsyn/Dsyn) and deltaD (Dsyn-Dnonsyn) have positive correlation with the full vaccination rate (R2= 0.723 and 0.505, respectively) in 19 countries, indicating that purifying selection pressure of SARS-CoV-2 spike gene increased as the vaccination coverage rate increased.
Data Availability
1.The World Health organization. Tracking SARS-CoV-2 variants. 2. Global Initiative on Sharing All Influenza Data. 3. Coronavirus (COVID-19) Vaccinations. Our World in Data.
https://www.who.int/en/activities/tracking-SARS-CoV-2-variants/