To reveal the distribution and transmission pathway of Paulownia witches’-broom (PaWB) disease, which is caused by phytoplasmas related to genetic variation, and the adaptability to the hosts and environments of the pathogenic population in different geographical regions in China, in this study, we used ten housekeeping gene fragments, including rp, fusA, secY, tuf, secA, dnaK, rpoB, pyrG, gyrB, and ipt, for multilocus sequence typing (MLST). A total of 142 PaWB phytoplasma strains were collected from 18 provinces or municipalities. The results showed that the genetic diversity was comparatively higher among the PaWB phytoplasma strains, and substantially different from that of the other 16SrI subgroup strains. The number of gene variation sites for different housekeeping genes in the PaWB phytoplasma strains ranged from 1 to 14 SNPs. Among them, rpoB (1.47%) and dnaK (1.12%) had higher genetic variation, and rp (0.20%) had the least genetic variation. The tuf and rpoB genes showed the fixation of positively selected beneficial mutations in the PaWB phytoplasma populations, and all housekeeping genes except tuf followed the neutral evolutionary model. We found an absence of recombination among PaWB phytoplasma sequence types (STs) for each housekeeping gene except dnaK, and no evidence for such recombination events for concatenated sequences of PaWB phytoplasma strains. The 22 sequence types were identified among the concatenated sequences of seven housekeeping genes (rp, fusA, secY, secA, tuf, dnaK, and rpoB) from 105 representative strains. We analyzed all 22 STs by goeBURST algorithm, forming two clonal complexes (CCs) and three singletons. Among them, ST1, as the primary founder of CC1, had the widest geographical distribution, accounting for 72.38% of all strains, with a high frequency of shared sequence type. The results of phylogenetic analysis of the concatenated sequences further revealed that the 105 strains were clustered into two representative lineages of PaWB phytoplasma, with obvious geographical differentiation. The ST1 strains of highly homogeneous lineage-1 were a widespread and predominant population in diseased areas. Lineage-2 contained strains from Jiangxi, Fujian, and Shaanxi provinces, highlighting the close genetic relatedness of the strains in these regions, which was also consistent with the results of most single-gene phylogenetic analysis of each gene. We also found that the variability in the northwest China population was higher than in other geographical populations; the range of genetic differentiation between the south of the Yangtze River population and the Huang-huai-hai Plain (or southwest China) population was relatively large. The achieved diversity and evolution data, as well as the MLST technique, are helpful for epidemiological studies and guiding PaWB disease control decisions.

The thymidylate kinase (tmk) gene is indispensable for the proliferation and survival of phytoplasma. To reveal the molecular variation and phylogeny of the tmk genes of Candidatus phytoplasma ziziphi, in this study, the tmk genes of 50 phytoplasma strains infecting different resistant and susceptible jujube cultivars from different regions in China were amplified and analyzed. Two sequence types, tmk-x and tmk-y, were identified using clone-based sequencing. The JWB phytoplasma strains were classified into three types, type-X, type-Y, and type-XY, based on the sequencing chromatograms of the tmk genes. The type-X and type-Y strains contained only tmk-x and tmk-y genes, respectively. The type-XY strain contained both tmk-x and tmk-y genes. The type-X, type-Y, and type-XY strains comprised 42%, 12%, and 46% of all the strains, respectively. The type-X and type-XY strains were identified in both susceptible and resistant jujube cultivars, while type-Y strain was only identified in susceptible cultivars. Phylogenetic analysis indicated that the tmk genes of the phytoplasmas were divided into two categories: phylo-S and phylo-M. The phylo-S tmk gene was single-copied in the genome, with an evolutionary pattern similar to the 16S rRNA gene; the phylo-M tmk gene was multi-copied, related to PMU-mediated within-genome transposition and between-genome transfer. Furthermore, the phylogenetic tree suggested that the tmk genes shuttled between the genomes of the Paulownia witches’ broom phytoplasma and JWB phytoplasma. These findings provide insights into the evolutionary and adaptive mechanisms of phytoplasmas.

Background: Paulownia witches’ broom disease (PaWB), caused by phytoplasmas, results in economic losses. Thymidine kinase (TDK) is a key enzyme in the deoxythymidine triphosphate (dTTP) synthesis pathway of phytoplasma, and dTTP serves as an essential precursor for DNA replication during cell proliferation, but its role in regulating phytoplasma proliferation and virulence remains unclear. Results: To characterize the role of TDK in phytoplasma proliferation and virulence, in this study, two PaWB phytoplasma tdk genotypes were identified via sampling in China, and tdk-1 was the most widely distributed. Phylogenetic analysis revealed that single-copied tdk exhibited higher phylogenetic resolution than 16S rRNA for phytoplasma molecular classification. In silico prediction showed that TDK functions as a thymidine kinase in these phytoplasmas. Furthermore, TDK protein expression was detected in PaWB-infected tissues. Immunofluorescence microscopy showed that TDK co-localized with phytoplasmas in sieve tube elements. Immunoelectron microscopy showed that TDK mainly localized in dividing phytoplasma cells. Absolute quantitative analysis revealed that tdk mRNA expression was positively associated with phytoplasma proliferation. Overexpressing tdk in Pseudomonas syringae pv. tomato ( Pst ) DC3000 significantly enhanced bacterial growth and increased its virulence in Arabidopsis . Conclusions: Overall, these results identify phytoplasma TDK as a critical proliferation-associated pathogenicity factor that promotes pathogen proliferation, providing a potential molecular target for controlling PaWB.