Many pathogenic bacteria are recognized as species complexes and uncertainties regarding the organization of their genetic diversity are challenges for research efforts. Within Agrobacterium tumefaciens, multiple genomospecies have been identified; however, the exact species boundaries are unclear, which causes chaos in nomenclature and hampers communication. In this work, we conducted targeted genome sequencing to achieve a comprehensive and balanced taxon sampling within this complex. Our results from genome-wide sequence identity, core genome phylogeny, and gene content not only supported that those recognized genomospecies are distinct biological entities but also identified novel genomospecies. Based on the fully resolved phylogeny, we further investigated the evolution of genes critical in Agrobacterium fitness and ecology. For the type VI secretion system (T6SS) involved in interbacterial competition, multiple losses and one horizontal gene transfer (HGT) event were inferred. For the tumor-inducing plasmids (pTi) and the pTi-encoded type IV secretion system (T4SS) that determine Agrobacterium phytopathogenicity, the evolution of these accessory replicons was decoupled from the chromosomes, thus contributing to another level of complexity. Intriguingly, for both T6SS and T4SS, genes that encode the structural components are highly conserved, whereas extensive diversity exists at multiple levels (i.e., between-species, within-species, intra-genome, and intra-gene) for genes that encode effectors and associated proteins. These findings suggest that opposite modes of selection may act on components conferring different functions within a system. In conclusion, this work provides insights into the genomic diversification of these bacteria and sheds light on the modularity of their molecular evolution.

The economically important plant pathogen Xylella fastidiosa has been reported in multiple regions of the globe during the last two decades, threatening a growing list of crops and industries. Xylella fastidiosa subspecies fastidiosa causes disease in grapevines (Pierce's disease of grapevines, PD), a current problem in the United States (US), Spain, and Taiwan. We studied PD-causing subsp. fastidiosa populations and compared the genome sequences of 33 isolates found in Central Taiwan with 171 isolates from the US and two from Spain. Phylogenetic relationships, haplotype network, and genetic diversity analyses confirm that subsp. fastidiosa was recently introduced into Taiwan from the Southeast US (i.e., the PD-I lineage in Georgia based on available data). Recent core genome recombination events were detected among introduced subsp. fastidiosa isolates in Taiwan and contributed to the development of genetic diversity, particularly in the Houli District of Taichung City in Central Taiwan. Unexpectedly, despite comprehensive sampling of all regions with high PD incidences in Taiwan, the genetic diversity observed include contributions through recombination from unknown donors, suggesting that higher diversity exists in the region. Nevertheless, no recombination event was detected between X. fastidiosa subsp. fastidiosa and the endemic sister species Xylella taiwanensis. In summary, this study improved our understanding of the genetic diversity of PD-causing subsp. fastidiosa after invasion to a new region.

Phytoplasmas are insect-transmitted plant pathogens that cause substantial losses in agriculture. In addition to economic impact, phytoplasmas induce distinct disease symptoms in infected plants, thus attracting attention for research on molecular plant-microbe interactions and plant developmental processes. Due to the difficulty of establishing an axenic culture of these bacteria, culture-independent genome characterization is a crucial tool for phytoplasma research. However, phytoplasma genomes have strong nucleotide composition biases and are repetitive, which make it challenging to produce complete assemblies. In this study, we utilized Illumina and Oxford Nanopore sequencing technologies to obtain the complete genome sequence of 'Candidatus Phytoplasma luffae' strain NCHU2019 that is associated with witches' broom disease of loofah (Luffa aegyptiaca) in Taiwan. The fully assembled circular chromosome is 769 kb in size and is the first representative genome sequence of group 16SrVIII phytoplasmas. Comparative analysis with other phytoplasmas revealed that NCHU2019 has an exceptionally repetitive genome, possessing a pair of 75 kb repeats and at least 13 potential mobile units (PMUs) that account for ~25% of its chromosome. This level of genome repetitiveness is exceptional for bacteria, particularly among obligate pathogens with reduced genomes. Our genus-level analysis of PMUs demonstrated that these phytoplasma-specific mobile genetic elements can be classified into three major types that differ in gene organization and phylogenetic distribution. Notably, PMU abundance explains nearly 80% of the variance in phytoplasma genome sizes, a finding that provides a quantitative estimate for the importance of PMUs in phytoplasma genome variability. Finally, our investigation found that in addition to horizontal gene transfer, PMUs also contribute to intra-genomic duplications of effector genes, which may provide redundancy for neofunctionalization or subfunctionalization. Taken together, this work improves the taxon sampling for phytoplasma genome research and provides novel information regarding the roles of mobile genetic elements in phytoplasma evolution.