Abstract
Highly transmissible or immuno-evasive SARS-CoV-2 variants have intermittently emerged and outcompeted previously circulating strains, resulting in repeated COVID-19 surges, reinfections, and breakthrough infections in vaccinated individuals. With over 5 million SARS-CoV-2 genomes sequenced globally over the last 2 years, there is unprecedented data to decipher how competitive viral evolution results in the emergence of fitter SARS-CoV-2 variants. Much attention has been directed to studying how specific mutations in the Spike protein impact its binding to the ACE2 receptor or viral neutralization by antibodies, but there is limited knowledge of genomic signatures shared primarily by dominant variants. Here we introduce a methodology to quantify the genome-wide distinctiveness of polynucleotide fragments of various lengths (3- to 240-mers) that constitute SARS-CoV-2 lineage genomes. Compared to standard phylogenetic distance metrics and overall mutational load, the quantification of distinctive 9-mer polynucleotides provides a higher resolution of separation between variants of concern (Reference = 89, IQR: 65-108; Alpha = 166, IQR: 150-182; Beta 130, IQR: 113-147; Gamma = 165, IQR: 152-180; Delta = 234, IQR: 216-253; and Omicron = 294, IQR: 287-315). The similar scoring of the Alpha and Gamma variants by our methodology is consistent with these strains emerging at approximately the same time and circulating in distinct geographical regions as dominant strains. Furthermore, evaluation of genomic distinctiveness for 1,363 lineages annotated in GISAID highlights that polynucleotide diversity has increased over time (R2 = 0.37) and that VOCs show high distinctiveness compared to non-VOC contemporary lineages. To facilitate similar real-time assessments on the competitive fitness potential of future variants, we are launching a freely accessible resource for infusing pandemic preparedness with genomic inference (“GENI” – https://academia.nferx.com/GENI). This study demonstrates the value of characterizing new SARS-CoV-2 variants by their genome-wide polynucleotide distinctiveness and emphasizes the need to go beyond a narrow set of mutations at known functionally salient sites.
Competing Interest Statement
KM, MJMN, BR, PJL, PG, PA, TF, SS, RS, TJH and VS are employees of nference and have financial interests in the company. nference collaborates with bio-pharmaceutical companies, health systems, and academic medical centers on data science initiatives unrelated to this study. These collaborations had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. KM, MJMN, BR, and VS are named inventors on a provisional patent application related to this study.
Funding Statement
This study was funded by nference (https://nference.com) and did not receive any external 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
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
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
↵+ Joint first authors
Data Availability
All data produced are available online at https://academia.nferx.com/GENI