Abstract
The exceptionally widespread outbreak of human monkeypox, an emerging zoonosis caused by the monkeypox virus (MPXV), with more than 69,000 confirmed cases in 100 non-endemic countries since 2022, is a major public health concern. Codon usage patterns reflect genetic variation and adaptation to new hosts and ecological niches. However, detailed analyses of codon usage bias in MPXV based on large-scale genomic data, especially for strains responsible for the 2022 outbreak, are lacking. In this study, we analyzed codon usage in MPXV and its relationship with host adaptation. We confirmed the ongoing outbreak of MPXVs belonging to the West Africa (WA) lineage by principal component analysis based on their codon usage patterns. The 2022 outbreak strains had a relatively low codon usage bias. Codon usage of MPXVs was shaped by mutation and natural selection; however, different from past strains, codon usage in the 2022 outbreak strains was predominantly determined by mutation pressure. Additionally, as revealed by the codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) analyses, the codon usage patterns of MPXVs were also affected by their hosts. In particular, the 2022 outbreak strains showed slightly but significantly greater adaptation to many primates, including humans, and were subjected to stronger selection pressure induced by hosts. Our results suggest that MPXVs contributing to the 2022 outbreak have unique evolutionary features, emphasizing the importance of sustained monitoring of their transmission and evolution.