Genetic Diversity of SARS-CoV-2 in Kazakhstan from 2020 to 2022

Coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, has had major social and economic consequences worldwide. Whole genome sequencing (WGS) is essential for genomic monitoring, enabling tracking of viral evolution, detection of emerging variants, and identification of introductions and transmission chains to inform timely public health responses. Here, we compile and harmonize SARS-CoV-2 genomic data generated by multiple laboratories across Kazakhstan together with publicly available sequences to provide a national overview of genomic dynamics across successive epidemic waves from 2020 to 2022. We analyzed 4462 genomes deposited in GISAID (including 340 generated in this study), of which 3299 passed Nextclade quality filters, and summarized lineage turnover across major phases (pre-VOC, Alpha, Delta, Omicron BA.1/BA.2, Omicron BA.4/BA.5, and a later recombinant-dominant period). Sequencing intensity varied markedly over time (0.60‱ of confirmed cases during Delta vs. 11.57‱ during the Omicron BA.5 wave), suggesting that lineage diversity and persistence may be underestimated. Pre-VOC circulation included ≥12 Pango lineages with evidence of multiple introductions and sustained local transmission, including a Kazakhstan-restricted B.4.1 lineage that emerged in Nur-Sultan/Astana and disappeared after April 2020. The Tengizchevroil oilfield outbreak comprised B.1.1 viruses with phylogenetic support for ≥three independent introductions. Alpha and Omicron waves were characterized by repeated introductions and heterogeneous origins, whereas Delta was dominated by AY.122 with an additional distinct AY.122 cluster; a notable BF.7 local transmission event was observed during BA.5. We also highlight locally enriched non-lineage-defining mutations. Overall, recurrent importations and variable local amplification shaped SARS-CoV-2 dynamics in Kazakhstan, while interpretation is constrained by strongly time-skewed sequencing.