DOI: 10.6084/m9.figshare.29577704
title: <b>SARS-CoV-2 Spike Mutation Clusters Identified via Large-Scale Statistical Linkage Analysis,</b><i>Insights into Epistatic Interactions and Evolutionary Pathways</i>

description: This dataset presents the results of a large-scale mutation co-occurrence analysis of the SARS-CoV-2 spike protein, based on <b>158,342</b> high-quality genomes collected between 2022 and 2024. Mutation profiles were extracted using standard bioinformatics pipelines (MAFFT, iVar, SAMtools, BCFtools) and validated using the GVATLAS platform. A total of <b>3,823,995</b> unique mutation pairs were tested for non-random co-occurrence using the Chi-Square test of independence . Of these, <b>47,068 </b>pairs showed statistically significant co-occurrence (<b>p &lt; 0.05</b>), representing approximately <b>1.23%</b> of all pairs , indicating that while most mutations occur independently, a small but meaningful subset co-evolves under selective pressure.The most significant pair was <b>S:T478R </b>& <b>S:T547I</b> , with a p-value of approximately<b> 0.00</b> , showing strong linkage. Other top co-occurring pairs include S:A1015S & S:L24S , suggesting potential functional or adaptive relationships.These findings highlight potential mutation clusters that may influence immune evasion , transmissibility , and antigenic drift , and align with similar patterns observed in other RNA viruses such as RSV , MERS-CoV , and influenza .We have since expanded this dataset to include <b>6.8 million accessions</b>, enabling an even more comprehensive view of mutation dynamics across time and geographies. This large-scale extension allows deeper insights into the evolution of SARS-CoV-2 and provides a robust resource for future research on viral adaptation.<br><b>SARS-CoV-2 Spike Mutation Co-Occurrence Summary</b><br><b>Based on 5,525,457 accessions</b><br><b>Total mutations observed: 109825409</b>Unique mutations: <b>7,806</b><br><b>PAIRWISE ANALYSIS</b>-----------------Total mutation pairs analyzed: <b>540,368</b>Significant pairs (p &lt; 0.05): <b>249,836</b><br><b>TOP MUTATIONS (by frequency)</b>-------------------------------S:D614G 5,490,427S:T478K 4,383,840S:P681H 3,213,690S:N501Y 3,166,083S:H655Y 2,973,639S:N679K 2,962,488S:N969K 2,957,849S:Q954H 2,945,931S:G142D 2,935,728S:D796Y 2,922,508<br><b>TOP CO-OCCURRING MUTATION PAIRS (by count)</b>------------------------------------------------------------S:D614G | S:T478K | 4,370,242S:D614G | S:P681H | 3,204,312S:D614G | S:N501Y | 3,156,208S:D614G | S:H655Y | 2,967,845S:D614G | S:N679K | 2,954,582S:D614G | S:N969K | 2,947,503S:N969K | S:Q954H | 2,942,025S:N501Y | S:P681H | 2,939,596S:N679K | S:N969K | 2,939,375S:D614G | S:Q954H | 2,935,645<br><b>STRONGEST CO-OCCURRING MUTATION PAIRS (by Cramér's V)</b>------------------------------------------------------------S:N969K | S:Q954H | 0.702040S:S373P | S:S375F | 0.701917S:L981F | S:N856K | 0.701013S:L981F | S:T547K | 0.700331S:A570D | S:D1118H | 0.700198S:L452W | S:N481K | 0.699967S:A570V | S:P621S | 0.699760S:N856K | S:T547K | 0.698985S:A570V | S:E554K | 0.697428S:A570D | S:S982A | 0.696944-------------------------------------------------<b>Interpretation:</b><b>Over 249,836 significant mutation pairs </b>indicate that while many mutations occur independently, a large number show <b>non-random co-evolution patterns </b>.Many top pairs involve mutations located in key domains like:<b>Receptor Binding Domain (RBD) </b>– e.g., N501Y, T478K<b>Furin cleavage site </b>– e.g., P681H<b>NTD and other antigenic sites </b>– e.g., G142D, S371LThis supports the hypothesis that <b>selective immune pressure </b>and <b>viral adaptation </b>drive the emergence of specific mutation combinations.<b>For inquiries regarding genomic surveillance and collaborative research, please contact:</b><br>TahirHB@gotmail.com <br><br>