Abstract
Citrus tristeza virus (CTV) is an important pathogen threatening the global citrus industry, but its evolution and transmission mechanism in wild citrus has not been clarified. Most of the existing studies are based on CTV-specific gene fragments, lacking genome-wide analysis. There is especially a lack of understanding of CTV transmission dynamics in wild citrus, which needs further investigation. In this study, wild citrus samples from three provinces of China were collected, virus genome data were obtained by high-throughput sequencing (HTS) technology and combined with public database data, and Bayesian phylogeographic inference was used to analyze virus composition characteristics in wild citrus, as well as the population genetic structure, temporal dynamic evolution, and spatial transmission mode of CTV. The results showed that Yunnan wild citrus samples contained the most abundant virus components, including CTV, Citrus Exocortis Viroid (CEVd), Citrus associated Ampelovirus 1 (CaAV-1), and Citrus Virus B (CiVB), while Jiangxi and Hunan samples only contained CTV and CEVd, with all samples showing mixed infection. Phylogenetic analysis showed that nine wild citrus CTV isolates were scattered in different evolutionary clades, and only 9.27% of genetic variation existed between the populations, while 90.72% of genetic variation existed within the populations, indicating little effect of geographic isolation on gene flow. The time to the most recent common ancestor (tMRCA) of CTV was estimated at 1360 CE, with subsequent divergence into two lineages, with population size stabilizing after a rapid increase in 1980-1990. Asia has been identified as the central source of CTV's global spread, with key migration events including Asia to North America (1746), Asia to Oceania (1829), and Asia to South America (1965), coinciding with global maritime trade and the expansion of the citrus industry.