Time-calibrated genomic evolution of a monomorphic bacterium during its establishment as an endemic crop pathogen

Evolution and spread of Xanthomonas citri subsp. citri in the São Paulo, Brazil, citrus belt inferred from 758 novel genomes

The São Paulo state citrus belt in Brazil is a major citrus production region. Since at least 1957, citrus plantations in this region have been affected by citrus canker, an economically damaging disease caused by Xanthomonas citri subsp. citri (Xcc). For about 50 years, until 2017, a citrus canker eradication programme was carried out in this region. In this work, our aim was to investigate the effects of the eradication programme on genetic variability and evolution of Xcc. To this end, we sequenced and analysed 758 Xcc genomes sampled in the São Paulo citrus belt, together with 730 publicly available Xcc genomes from around the world. Our phylogenomic analyses show that these genomes can be grouped into seven major lineages and that in São Paulo, lineage L7 is dominant. Our time estimate for its appearance closely matches the date when citrus production expanded. L7 can be subdivided into lineages L7.1 and L7.2. In our samples, L7.2, which we estimate to have emerged around 1964, is by far the most abundant, showing that the eradication programme had little impact on strain diversification. On the other hand, oscillations in the estimated effective population size of L7.2 strains over time closely match the shifts in the eradication programme. In sum, we present a detailed view of the genomic diversity of Xcc in the world and in São Paulo, the largest such effort in terms of a number of genomes for a crop pathogen undertaken so far. The methods employed here can form the basis for active genomic surveillance of Xcc in major citrus production areas.

Genetic Signatures of Contrasted Outbreak Histories of "Candidatus Liberibacter asiaticus", the Bacterium That Causes Citrus Huanglongbing, in Three Outermost Regions of the European Union

In an era of trade globalization and climate change, crop pathogens and pests are a genuine threat to food security. The detailed characterization of emerging pathogen populations is a prerequisite for managing invasive species pathways and designing sustainable disease control strategies. Huanglongbing is the disease that causes the most damage to citrus, a crop that ranks #1 worldwide in terms of fruit production. Huanglongbing can be caused by three species of the phloem-limited alpha-proteobacterium, "Candidatus Liberibacter," which are transmitted by psyllids. Two of these bacteria are of highest concern, "Ca. Liberibacter asiaticus" and "Ca. Liberibacter africanus," and have distinct thermal optima. These pathogens are unculturable, which complicates their high-throughput genetic characterization. In the present study, we used several genotyping techniques and an extensive sample collection to characterize Ca. Liberibacter populations associated with the emergence of huanglongbing in three French outermost regions of the European Union (Guadeloupe, Martinique and Réunion). The outbreaks were primarily caused by "Ca. Liberibacter asiaticus," as "Ca. Liberibacter africanus" was only found at a single location in Réunion. We emphasize the low diversity and high genetic relatedness between samples from Guadeloupe and Martinique, which suggests the putative movement of the pathogen between the two islands and/or the independent introduction of closely related strains. These samples were markedly different from the samples from Réunion, where the higher genetic diversity revealed by tandem-repeat markers suggests that the disease was probably overlooked for years before being officially identified in 2015. We show that "Ca. Liberibacter asiaticus" occurs from sea level to an altitude of 950 m above sea level and lacks spatial structure. This suggests the pathogen's medium- to long-distance movement. We also suggest that backyard trees acted as relays for disease spread. We discuss the implications of population biology data for surveillance and management of this threatful disease.

Amylase-associated genetic pattern in Xanthomonas euvesicatoria on pepper

Bacterial leaf spot of pepper (BSP), primarily caused by Xanthomonas euvesicatoria (Xe), poses a significant challenge to pepper production worldwide. Despite its impact, the genetic diversity of this pathogen remains underexplored, which limits our understanding of its population structure. To bridge this knowledge gap, we conducted a comprehensive analysis using 103 Xe strains isolated from pepper in southwest Florida to characterize genomic and type III effector (T3E) variation in this population. Phylogenetic analysis of core genomes revealed a major distinct genetic lineage associated with amylolytic activity. This amylolytic lineage was represented in Xe strains globally. Molecular clock analysis dated the emergence of amylolytic strains in Xe to around 1972. Notably, non-amylolytic strains possessed a single base pair frameshift deletion in the ⍺-amylase gene yet retained a conserved C-terminus. GUS assay revealed the expression of two open reading frames in non-amylolytic strains, one at the N-terminus and another that starts 136 base pairs upstream of the ⍺-amylase gene. Analysis of T3Es in the Florida Xe population identified variation in 12 effectors, including two classes of mutations in avrBs2 that prevent AvrBs2 from triggering a hypersensitive response in Bs2-resistant pepper plants. Knowledge of T3E variation could be used for effector-targeted disease management. This study revealed previously undescribed population structure in this economically important pathogen.IMPORTANCEBacterial leaf spot (BSP), a significant threat to pepper production globally, is primarily caused by Xanthomonas euvesicatoria (Xe). Limited genomic data has hindered detailed studies on its population diversity. This study analyzed the whole-genome sequences of 103 Xe strains from peppers in southwest Florida, along with additional global strains, to explore the pathogen's diversity. The study revealed two major distinct genetic groups based on their amylolytic activity, the ability to break down starch, with non-amylolytic strains having a mutation in the ⍺-amylase gene. Additionally, two classes of mutations in the avrBs2 gene were found, leading to susceptibility in pepper plants with the Bs2 resistance gene, a commercially available resistance gene for BSP. These findings highlight the need to forecast the emergence of such strains, identify genetic factors for innovative disease management, and understand how this pathogen evolves and spreads.

Xanthomonas as a Model System for Studying Pathogen Emergence and Evolution

In this review, we highlight studies in which whole-genome sequencing, comparative genomics, and population genomics have provided unprecedented insights into past and ongoing pathogen evolution. These include new understandings of the adaptive evolution of secretion systems and their effectors. We focus on Xanthomonas pathosystems that have seen intensive study and improved our understanding of pathogen emergence and evolution, particularly in the context of host specialization: citrus canker, bacterial blight of rice, and bacterial spot of tomato and pepper. Across pathosystems, pathogens appear to follow a pattern of bursts of evolution and diversification that impact host adaptation. There remains a need for studies on the mechanisms of host range evolution and genetic exchange among closely related but differentially host-specialized species and to start moving beyond the study of specific strain and host cultivar pairwise interactions to thinking about these pathosystems in a community context.

Bacterial Mutation During Seasonal Epidemics

Rapidly evolving bacterial pathogens pose a unique challenge for long-term plant disease management. In this study, we investigated the types and rate of mutations in bacterial populations during seasonal disease epidemics. Two phylogenetically distinct strains of the bacterial spot pathogen, Xanthomonas perforans, were marked, released in tomato fields, and recaptured at several time points during the growing season. Genomic variations in recaptured isolates were identified by comparative analysis of their whole-genome sequences. In total, 180 unique variations (116 substitutions, 57 insertions/deletions, and 7 structural variations) were identified from 300 genomes, resulting in the overall host-associated mutation rate of ∼0.3 to 0.9/genome/week. This result serves as a benchmark for bacterial mutation during epidemics in similar pathosystems. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Population Genomics Reveals an Emerging Lineage of Xanthomonas perforans on Pepper

Xanthomonas perforans-the dominant causal agent of bacterial leaf spot of tomato-is an emerging pathogen of pepper, indicative of a potential host expansion across the southeastern United States. However, studies of the genetic diversity and evolution of X. perforans from pepper remain limited. In this study, the whole-genome sequences of 35 X. perforans strains isolated from pepper from four fields and two transplant facilities across southwest Florida between 2019 and 2021 were used to compare genomic divergence, evolution, and variation in type III secreted effectors. Phylogenetic analysis based on core genes revealed that all 35 X. perforans strains formed one genetic cluster with pepper and tomato strains from Alabama and Turkey and were closely related to strains isolated from tomato in Indiana, Mexico, and Louisiana. The in planta population growth of tomato strains isolated from Indiana, Mexico, Louisiana, and Turkey in pepper leaf mesophyll was on par with pepper X. perforans and X. euvesicatoria strains. Molecular clock analysis of the 35 Florida strains dated their emergence to approximately 2017. While strains varied in copper tolerance, all sequenced strains harbored the avrHah1 transcription activation-like effector located on a conjugative plasmid, not previously reported in Florida. Our findings suggest that there is a geographically distributed lineage of X. perforans strains on tomato that has the genetic background to cause disease on pepper. Moreover, this study clarifies potential adaptive variants of X. perforans on pepper that could help forecast the emergence of such strains and enable immediate or preemptive intervention.

Herbarium specimen sequencing allows precise dating of Xanthomonas citri pv. citri diversification history

Herbarium collections are an important source of dated, identified and preserved DNA, whose use in comparative genomics and phylogeography can shed light on the emergence and evolutionary history of plant pathogens. Here, we reconstruct 13 historical genomes of the bacterial crop pathogen Xanthomonas citri pv. citri (Xci) from infected Citrus herbarium specimens. Following authentication based on ancient DNA damage patterns, we compare them with a large set of modern genomes to estimate their phylogenetic relationships, pathogenicity-associated gene content and several evolutionary parameters. Our results indicate that Xci originated in Southern Asia ~11,500 years ago (perhaps in relation to Neolithic climate change and the development of agriculture) and diversified during the beginning of the 13th century, after Citrus diversification and before spreading to the rest of the world (probably via human-driven expansion of citriculture through early East-West trade and colonization).

Clustered Regularly Interspaced Short Palindromic Repeats in Xanthomonas citri—Witnesses to a Global Expansion of a Bacterial Pathogen over Time

Xanthomonas citri pv. citri, a Gram-negative bacterium, is the causal agent of citrus canker, a significant threat to citrus production. Understanding of global expansion of the pathogen and monitoring introduction into new regions are of interest for integrated disease management at the local and global level. Genetic diversity can be assessed using genomic approaches or information from partial gene sequences, satellite markers or clustered regularly interspaced short palindromic repeats (CRISPR). Here, we compared CRISPR loci from 355 strains of X. citri pv. citri, including a sample from ancient DNA, and generated the genealogy of the spoligotypes, i.e., the absence/presence patterns of CRISPR spacers. We identified 26 novel spoligotypes and constructed their likely evolutionary trajectory based on the whole-genome information. Moreover, we analyzed ~30 additional pathovars of X. citri and found that the oldest part of the CRISPR array was present in the ancestor of several pathovars of X. citri. This work presents a framework for further analyses of CRISPR loci and allows drawing conclusions about the global spread of the citrus canker pathogen, as exemplified by two introductions in West Africa.

Pathotyping Citrus Ornamental Relatives with Xanthomonas citri pv. citri and X. citri pv. aurantifolii Refines Our Understanding of Their Susceptibility to These Pathogens

Xanthomonas citri pv. citri (Xcc) and X. citri pv. aurantifolii (Xca) are causal agents of Citrus Bacterial Canker (CBC), a devastating disease that severely affects citrus plants. They are harmful organisms not reported in Europe or the Mediterranean Basin. Host plants are in the Rutaceae family, including the genera Citrus, Poncirus, and Fortunella, and their hybrids. In addition, other genera of ornamental interest are reported as susceptible, but results are not uniform and sometimes incongruent. We evaluated the susceptibility of 32 ornamental accessions of the Rutaceae family belonging to the genera Citrus, Fortunella, Atalantia, Clausena, Eremocitrus, Glycosmis, Microcitrus, Murraya, Casimiroa, Calodendrum, and Aegle, and three hybrids to seven strains of Xcc and Xca. Pathotyping evaluation was assessed by scoring the symptomatic reactions on detached leaves. High variability in symptoms and bacterial population was shown among the different strains in the different hosts, indicative of complex host–pathogen interactions. The results are mostly consistent with past findings, with the few discrepancies probably due to our more complete experimental approach using multiple strains of the pathogen and multiple hosts. Our work supports the need to regulate non-citrus Rutaceae plant introductions into areas, like the EU and Mediterranean, that are currently free of this economically important pathogen.

Comparative Genomic Analysis Uncovered Evolution of Pathogenicity Factors, Horizontal Gene Transfer Events, and Heavy Metal Resistance Traits in Citrus Canker Bacterium Xanthomonas citri subsp. citri

Background: Worldwide citrus production is severely threatened by Asiatic citrus canker which is caused by the proteobacterium Xanthomonas citri subsp. citri. Foliar sprays of copper-based bactericides are frequently used to control plant bacterial diseases. Despite the sequencing of many X. citri strains, the genome diversity and distribution of genes responsible for metal resistance in X. citri subsp. citri strains from orchards with different management practices in Taiwan are not well understood. Results: The genomes of three X. citri subsp. citri strains including one copper-resistant strain collected from farms with different management regimes in Taiwan were sequenced by Illumina and Nanopore sequencing and assembled into complete circular chromosomes and plasmids. CRISPR spoligotyping and phylogenomic analysis indicated that the three strains were located in the same phylogenetic lineages and shared ∼3,000 core-genes with published X. citri subsp. citri strains. These strains differed mainly in the CRISPR repeats and pathogenicity-related plasmid-borne transcription activator-like effector (TALE)-encoding pthA genes. The copper-resistant strain has a unique, large copper resistance plasmid due to an unusual ∼40 kbp inverted repeat. Each repeat contains a complete set of the gene cluster responsible for copper and heavy metal resistance. Conversely, the copper sensitive strains carry no metal resistance genes in the plasmid. Through comparative analysis, the origin and evolution of the metal resistance clusters was resolved. Conclusion: Chromosomes remained constant among three strains collected in Taiwan, but plasmids likely played an important role in maintaining pathogenicity and developing bacterial fitness in the field. The evolution of pathogenicity factors and horizontal gene transfer events were observed in the three strains. These data suggest that agricultural management practices could be a potential trigger for the evolution of citrus canker pathogens. The decrease in the number of CRISPR repeats and pthA genes might be the result of adaptation to a less stressful environment. The metal resistance genes in the copper resistant X. citri strain likely originated from the Mauritian strain not the local copper-resistant X. euvesicatoria strain. This study highlights the importance of plasmids as 'vehicles' for exchanging genetic elements between plant pathogenic bacteria and contributing to bacterial adaptation to the environment.

First historical genome of a crop bacterial pathogen from herbarium specimen: Insights into citrus canker emergence

Over the past decade, ancient genomics has been used in the study of various pathogens. In this context, herbarium specimens provide a precious source of dated and preserved DNA material, enabling a better understanding of plant disease emergences and pathogen evolutionary history. We report here the first historical genome of a crop bacterial pathogen, Xanthomonas citri pv. citri (Xci), obtained from an infected herbarium specimen dating back to 1937. Comparing the 1937 genome within a large set of modern genomes, we reconstructed their phylogenetic relationships and estimated evolutionary parameters using Bayesian tip-calibration inferences. The arrival of Xci in the South West Indian Ocean islands was dated to the 19th century, probably linked to human migrations following slavery abolishment. We also assessed the metagenomic community of the herbarium specimen, showed its authenticity using DNA damage patterns, and investigated its genomic features including functional SNPs and gene content, with a focus on virulence factors.