ABSTRACT
Background SARS-CoV-2 is continuously spreading worldwide at an unprecedented scale and evolved into seven clades according to GISAID where four (G, GH, GR and GV) are globally prevalent in 2020. These major predominant clades of SARS-CoV- 2 are continuously increasing COVID-19 cases worldwide; however, after an early rise in 2020, the death-case ratio has been decreasing to a plateau. G clade viruses contain four co- occurring mutations in their genome (C241T+C3037T+C14408T: RdRp.P323L+A23403G:spike.D614G). GR, GH, and GV strains are defined by the presence of these four mutations in addition to the clade-featured mutation in GGG28881- 28883AAC:N. RG203-204KR, G25563T:ORF3a.Q57H, and C22227T:spike.A222V+C28932T-N.A220V+G29645T, respectively. The research works are broadly focused on the spike protein mutations that have direct roles in receptor binding, antigenicity, thus viral transmission and replication fitness. However, mutations in other proteins might also have effects on viral pathogenicity and transmissibility. How the clade- featured mutations are linked with viral evolution in this pandemic through gearing their fitness and virulence is the main question of this study.
Methodology We thus proposed a hypothetical model, combining a statistical and structural bioinformatics approach, endeavors to explain this infection paradox by describing the epistatic effects of the clade-featured co-occurring mutations on viral fitness and virulence.
Results and Discussion The G and GR/GV clade strains represent a significant positive and negative association, respectively, with the death-case ratio (incidence rate ratio or IRR = 1.03, p <0.001 and IRR= 0.99/0.97, p < 0.001), whereas GH clade strains showed no association with the Docking analysis showed the higher infectiousness of a spike mutant through more favorable binding of G614 with the elastase-2. RdRp mutation p.P323L significantly increased genome-wide mutations (p<0.0001) since more expandable RdRp (mutant)-NSP8 interaction may accelerate replication. Superior RNA stability and structural variation at NSP3:C241T might impact upon protein or RNA interactions. Another silent 5’UTR:C241T mutation might affect translational efficiency and viral packaging. These G- featured co-occurring mutations might increase the viral load, alter immune responses in host and hence can modulate intra-host genomic plasticity. An additional viroporin ORF3a:p.Q57H mutation, forming GH-clade, prevents ion permeability by cysteine (C81)- histidine (H57) inter-transmembrane-domain interaction mediated tighter constriction of the channel pore and possibly reduces viral release and immune response. GR strains, four G clade mutations and N:p.RG203-204KR, would have stabilized RNA interaction by more flexible and hypo-phosphorylated SR-rich region. GV strains seemingly gained the evolutionary advantage of superspreading event through confounder factors; nevertheless, N:p.A220V might affect RNA binding.
Conclusion These hypotheses need further retrospective and prospective studies to understand detailed molecular and evolutionary events featuring the fitness and virulence of SARS-CoV-2.
Highlights
We speculated an association of particular SARS-CoV-2 clade with death rate.
The polymerase mutant virus can speed up replication that corresponds to higher mutations.
The impact on viral epistasis by evolving mutations in SARS-CoV-2.
How the virus changes its genotype and circulate with other types given the overall dynamics of the epidemics?
Human intervention seems to work well to control the viral virulence. This hygiene practice will control the overall severity of the pandemic situation as recommended by the WHO. Our work has given the same message but explain with the dominant co-occurring mutations.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
We get no funds to carry out this research work
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:
Our University (Jashore University of Science and Technology) has provided ethical permission to do the work. However, this is a dry lab work.
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
Data Availability
All the sequence data were taken from the GISAID (https://www.gisaid.org/) and RCSB PDB (https://www.rcsb.org/) as mentioned in the methodology section. We provide all the necessary information such as accession numbers, date-based data source for helping readers and reviewers to check the authencity of the work.