Epidemiology of TMUV
To investigate the phylogeny trend of TMUV, 40 strains were obtained from GenBank, and were used as representative strains of the complete genome of TMUV obtained from mainland China from 2010 to 2019. All 40 sequences were aligned using ClustalW, which was applied in MEGA 6. The best substitution model was selected by jModelTest according to the Bayesian information criterion (BIC) score. Furthermore, to study the epidemiological trends of TMUV, a Bayesian skyline plot (BSP) was established. Additionally, to ensure the feasibility of the results, two independent runs of the chain were performed and compared to ensure adequate sampling; a chain length comprised 10,000,000 generations with sampling every 1,000 generations with 10% burn-in removed.