Abstract
There were five epidemic waves of coronavirus disease 2019 in Japan between 2020 and 2021. It remains unclear how the domestic waves arose and abated. To better understand this, we analyzed the pangenomic sequences of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and characterized the molecular epidemiological features of the five epidemic waves in Japan. In this study, we performed deep sequencing to determine the pangenomic SARS-CoV-2 sequences of 1,286 samples collected in two cities far from each other, Tokyo Metropolis and Nagoya. Then, the spatiotemporal genetic changes of the obtained sequences were compared with the sequences available in the Global Initiative on Sharing All Influenza Data (GISAID) database. A total of 873 genotypes carrying different sets of mutations were identified in the five epidemic waves. Phylogenetic analysis demonstrated that sharp displacements of lineages and genotypes occurred between consecutive waves over the 2 years. In addition, a wide variety of genotypes were observed in the early half of each wave, whereas a few genotypes were detected across Japan during an entire wave. Phylogenetically, putative descendant genotypes observed late in each wave displayed regional clustering and evolution in Japan. The genetic diversity of SARS-CoV-2 displayed uneven dynamics during each epidemic wave in Japan. Our findings provide an important molecular epidemiological basis to aid in controlling future SARS-CoV-2 epidemics.