Orthogroup-Based Comparative Analysis of Prophage Gene Content in Candidatus Liberibacter Asiaticus Supports a Predominantly Conserved Global Repertoire with Limited Accessory Variation

Huanglongbing, a destructive citrus disease of global importance that is also present in Saudi Arabia, is associated with Candidatus Liberibacter asiaticus (CLas) and remains a major threat to citrus production. Although previous studies have documented sequence variation and prophage polymorphism in CLas, broader comparisons of prophage-associated gene content remain limited. In particular, comparative orthogroup analysis of prophage gene-content conservation across geographically structured CLas populations has rarely been explored. In this study, we analyzed 42 CLas prophage genomes from Saudi Arabia and other geographic regions using a comparative orthogroup framework. OrthoFinder assigned 99.1% of predicted proteins (1825 of 1841) to 64 orthogroups, with only 16 genes remaining unassigned. A small number of rare orthogroups restricted to only a few genomes were identified, and no orthogroup was detected in all genomes. Presence–absence analyses supported a predominantly conserved prophage gene repertoire together with a small accessory component, while also indicating that apparent absences should be interpreted in light of mixed assembly status and prophage-region completeness. Saudi Arabian genomes were distributed within the broader global framework and exhibited generally similar gene-content profiles rather than a deeply separated lineage. Functional interpretation of representative orthogroups identified conserved prophage-associated genes related to replication, helicase activity, and phage packaging, whereas variable orthogroups were primarily associated with hypothetical or accessory prophage-related functions. Overall, these results are consistent with a model in which CLas prophage diversification is associated more with sequence-level variation and localized structural differences than with extensive gain or loss of prophage genes. These findings further refine current understanding of CLas genome evolution and highlight conserved prophage-associated targets that may support molecular diagnostics and epidemiological surveillance.