Molecular evolution and genomic characteristics of genotype H1 of measles virus

Measles is one of the most infectious diseases of humans. It is caused by the measles virus (MeV) and can lead to serious illness, lifelong complications, and even death. Whole-genome sequencing (WGS) is now available to study molecular epidemiology and identify MeV transmission pathways. In the present study, WGS of 23 MeV strains of genotype H1, collected in Mainland China between 2006 and 2018, were generated and compared to 31 WGSs from the public domain to analyze genomic characteristics, evolutionary rates and date of emergence of H1 genotype. The noncoding region between M and F protein genes (M/F NCR) was the most variable region throughout the genome. Although the nucleotide substitution rate of H1 WGS was around 0.75 × 10^-3 substitution per site per year, the M/F NCR had an evolutionary rate three times higher, with 2.44 × 10^-3 substitution per site per year. Phylogenetic analysis identified three distinct genetic groups. The Time of the Most Recent Common Ancestor (TMRCA) of H1 genotype was estimated at approximately 1988, while the first genetic group appeared around 1995 followed by two other genetic groups in 1999-2002. Bayesian skyline plot showed that the genetic diversity of the H1 genotype remained stable even though the number of MeV cases decreased 50 times between 2014 (52 628) and 2020 (993). The current coronavirus disease 2019 (COVID-19) pandemic might have some effect on the measles epidemic and further studies will be necessary to assess the genetic diversity of the H1 genotype in a post-COVID area.