Abstract
The devastating SARS-CoV2 pandemic is worsening with relapsing surges, emerging mutants, and increasing mortalities. Despite enormous efforts, it is not clear how SARS-CoV2 adapts and evolves in a clonal background. Laboratory research is hindered by high biosafety demands. However, the rapid sequence availability opened doors for bioinformatics. Using different bioinformatics programs, we investigated 6305 sequences for clonality, expressions strategies, and evolutionary dynamics. Results showed high nucleotide identity of 99.9% among SARS-CoV2 indicating clonal evolution and genome. High sequence identity and phylogenetic tree concordance were obtained with isolates from different regions. In any given tree topology, similar to 50% of isolates in a country formed country-specific sub-clusters. However, abundances of subtle overexpression strategies were found including transversions, signature-sequences and slippery-structures. Five different short tracks dominated with identical location patterns in all genomes where Slippery-4 AAGAA was the most abundant. Interestingly, transversion and transition substitutions mostly affected the same amino acid residues implying compensatory changes. To ensure these strategies were independent of sequence clonality, we simultaneously examined sequence homology indicators; tandem-repeats, restriction-site, and 3'UTR, 5' UTR-caps and stem-loop locations in addition to stringent alignment parameters for 100% identity which all confirmed stability. Nevertheless, two rare events; a rearrangement in two SARS-CoV2 isolates against betacoronavirus ancestor and a polymorphism in S gene, were detected. Thus, we report on abundance of transversions, slippery sequences, and ON/OFF molecular structures, implying adaptive expressions had occurred, despite clonal evolution and genome stability. Furthermore, functional validation of the point mutations would provide insights into mechanisms of SARS-CoV2 virulence and adaptation.