Ecological and Evolutionary Insights into Xanthomonas citri Pathovar Diversity

Phylogenomic analysis of 343 Xanthomonas citri pv. citri strains unravels introduction history and dispersal paths

Xanthomonas citri pv. citri (Xcc) causes the devastating citrus canker disease. Xcc is known to have been introduced into Florida, USA in at least three different events in 1915, 1986 and 1995 with the first two claimed to be eradicated. It was questioned whether the Xcc introduction in 1986 has been successfully eradicated. Furthermore, it is unknown how Xcc has spread throughout the citrus groves in Florida. In this study, we investigated the population structure of Xcc to address these questions. We sequenced the whole genome of 343 Xcc strains collected from Florida groves between 1997 and 2016. Our analysis revealed two distinct clusters of Xcc. Our data strongly indicate that the claimed eradication of the 1986 Xcc introduction was not successful and Xcc strains from 1986 introduction were present in samples from at least 8 counties collected after 1994. Importantly, our data revealed that the Cluster 2 strains, which are present in all 20 citrus-producing counties sampled in Florida, originated from the Xcc introduction event in the Miami area in 1995. Our data suggest that Polk County is the epicenter of the dispersal of Cluster 2 Xcc strains, which is consistent with the fact that three major hurricanes passed through Polk County in 2004. As copper-based products have been extensively used to control citrus canker, we also investigated whether Xcc strains have developed resistance to copper. Notably, none of the 343 strains contained known copper resistance genes. Twenty randomly selected Xcc strains displayed sensitivity to copper. Overall, this study provides valuable insights into the introduction, eradication, spread, and copper resistance of Xcc in Florida.

Comparative genomics reveals the emergence of copper resistance in a non-pigmented Xanthomonas pathogen of grapevine

Xanthomonas citri pv. viticola (Xcv) is the causal agent of bacterial canker in grapevine. The pathogen is restricted to India, where it was first reported in the 1970s, and Brazil. In the present study, we report the first complete genome sequence of Xcv LMG965, which is a reference pathotype strain. We also report genome sequences of additional isolates from India and comparative genome-based studies of isolates from Brazil. Apart from revealing the monophyletic origin of the pathovar, we could also confirm a common frameshift mutation in a gene that is part of the Xanthomonadin pigment biosynthetic gene cluster in all the isolates. The comparative study also revealed multiple intrinsic copper resistance-related genes in Brazilian isolates, suggesting intense selection, possibly because of heavy and indiscriminate usage of copper as an antimicrobial agent in the orchards. There is also the association of a Tn3-like transposase in the vicinity of the copper resistance genes, indicating a potential for rapid diversification through horizontal gene transfer events. The findings, along with genomic resources, will allow for systematic genetic and functional studies of Xcv.

A novel BLUF photoreceptor modulates the Xanthomonas citri subsp. citri-host plant interaction

Plant-pathogen interaction is influenced by multiple environmental factors, including temperature and light. Recent works have shown that light modulates not only the defense response of plants but also the pathogens virulence. Xanthomonas citri subsp. citri (Xcc) is the bacterium responsible for citrus canker, an important plant disease worldwide. The Xcc genome presents four genes encoding putative photoreceptors: one bacteriophytochrome and three blue light photoreceptors, one LOV and two BLUFs (bluf1: XAC2120 and bluf2: XAC3278). The presence of two BLUFs proteins is an outstanding feature of Xcc. In this work we show that the bluf2 gene is functional. The mutant strain, XccΔbluf2, was constructed demonstrating that BLUF2 regulates swimming-type motility, adhesion to leaves, exopolysaccharide production and biofilm formation, features involved in the Xcc virulence processes. An important aspect during the plant-pathogen interaction is the oxidative response of the host and the consequent reaction of the pathogen. We observed that ROS detoxification is regulated by Xcc bluf2 gene. The phenotypes of disease in orange plants produced by WT and XccΔbluf2 strains were evaluated, observing different phenotypes. Altogether, these results show that BLUF2 negatively regulates virulence during citrus canker. This work constitutes the first report on BLUF-like receptors in plant pathogenic bacteria.

Phylogenomic Analysis Supports the Transfer of 20 Pathovars from Xanthomonas campestris into Xanthomonas euvesicatoria

The Gram-negative bacterial genus Xanthomonas includes numerous infra-specific taxa known as pathovars, which are defined primarily on host range and disease symptoms. With the advent of molecular sequence data, many pathovars have been transferred from X. campestris into other Xanthomonas species to better harmonise taxonomy and phylogeny. We performed whole-genome shotgun sequencing on pathotype strains of the following X. campestris pathovars: blepharidis, carissae, clerodendri, convolvuli, coriandri, daturae, euphorbiae, fici, heliotropii, ionidii, lawsoniae, mirabilis, obscurae, paulliniae, pennamericanum, spermacoces, uppalii, vernoniae, viegasii and zingibericola. These genomes showed more than 98% average nucleotide identity with the type-strain of X. euvesicatoria and less than 88% with the type-strain of X. campestris. We propose the transfer of these pathovars into X. euvesicatoria and present an emended species description for X. euvesicatoria.

Citrus Canker Pathogen, Its Mechanism of Infection, Eradication, and Impacts

Citrus canker is a ravaging bacterial disease threatening citrus crops. Its major types are Asiatic Canker, Cancrosis B, and Cancrosis C, caused by Xanthomonas citri pv. citri (Xcc), Xanthomonas citri pv. aurantifolii pathotype-B (XauB), and pathotype-C (XauC), respectively. The bacterium enters its host through stomata and wounds, from which it invades the intercellular spaces in the apoplast. It produces erumpent corky necrotic lesions often surrounded by a chlorotic halo on the leaves, young stems, and fruits, which causes dark spots, defoliation, reduced photosynthetic rate, rupture of leaf epidermis, dieback, and premature fruit drop in severe cases. Its main pathogenicity determinant gene is pthA, whose variants are present in all citrus canker-causing pathogens. Countries where citrus canker is not endemic adopt different methods to prevent the introduction of the pathogen into the region, eradicate the pathogen, and minimize its dissemination, whereas endemic regions require an integrated management program to control the disease. The main aim of the present manuscript is to shed light on the pathogen profile, its mechanism of infection, and fruitful strategies for disease management. Although an adequate method to completely eradicate citrus canker has not been introduced so far, many new methods are under research to abate the disease.

Genome dynamics mediated by repetitive and mobile elements in Xanthomonas citri pv. durantae

Xanthomonas is a highly evolved group of phytopathogenic bacteria infecting nearly 400 host plants having vast genomic resources available with heterogenicity in representation from different species and pathovars. Unfortunately, the wealth of data is extremely biased and restricted to a few Xanthomonas pathogens that infect economically important plants, while those reported to infect the most diverse plants remain neglected. In the present study, we report the first complete genome sequence of Xanthomonas citri pv. durantae that was reported to infect Duranta repens L. or golden dewdrop, a hedge plant of ornamental importance native to the American region. Phylogenomic analysis with its closest relatives placed it amongst X. citri pv. citri A* pathotype strains and further comparative studies revealed various large unique genomic regions of chromosomal origin. The association of integrative and conjugative elements and prophages with unique genomic regions suggests the role of mobilome in genome dynamics. A large number of IS elements and transcription activator-like effectors encoding genes on each of the four plasmids indicate the further scope of diversification in Xanthomonas .

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.

Deep Population Genomics Reveals Systematic and Parallel Evolution at a Lipopolysaccharide Biosynthetic Locus in Xanthomonas Pathogens That Infect Rice and Sugarcane

The advent of high-throughput sequencing and population genomics has enabled researchers to investigate selection pressure at hypervariable genomic loci encoding pathogen-associated molecular pattern (PAMP) molecules like lipopolysaccharide (LPS). Xanthomonas is a model and a major group of phytopathogenic bacteria that infect hosts in tissue-specific manner. Our in-depth population-based genomic investigation revealed the emergence of major lineages in two Xanthomonas pathogens that infect xylem of rice and sugarcane is associated with the acquisition and later large-scale replacement by distinct type of LPS cassettes. In the population of the rice xylem pathogen, Xanthomonas oryzae pv. oryzae (Xoo) and sugarcane pathogens Xanthomonas sacchari (Xsac) and Xanthomonas vasicola (Xvv), the BXO8 type of LPS cassette is replaced by a BXO1 type of cassette in Xoo and by Xvv type LPS cassette in Xsac and Xvv. These findings suggest a wave of parallel evolution at an LPS locus mediated by horizontal gene transfer (HGT) events during its adaptation and emergence. Aside from xylem pathogens, two closely related lineages of Xoo that infect parenchyma of rice and Leersia hexandra grass have acquired an LPS cassette from Xanthomonas pathogens that infect parenchyma of citrus, walnut, and strawberries, indicating yet another instance of parallel evolution mediated by HGT at an LPS locus. Our targeted and megapopulation-based genome dynamic studies revealed the acquisition and dominance of specific types of LPS cassettes in adaptation and success of a major group of phytopathogenic bacteria. IMPORTANCE Lipopolysaccharide (LPS) is a major microbe associated molecular pattern and hence a major immunomodulator. As a major and outer member component, it is expected that LPS is a frontline defense mechanism to deal with different host responses. Limited studies have indicated that LPS loci are also highly variable at strain and species level in plant-pathogenic bacteria, suggesting strong selection pressure from plants and associated niches. The advent of high-throughput genomics has led to the availability of a large set of genomic resources at taxonomic and population levels. This provides an exciting and important opportunity to carryout megascale targeted and population-based comparative genomic/association studies at important loci like those encoding LPS biosynthesis to understand their role in the evolution of the host, tissue specificity, and also predominant lineages. Such studies will also fill major gap in understanding host and tissue specificity in pathogenic bacteria. Our pioneering study uses the Xanthomonas group of phytopathogens that are known for their characteristic host and tissue specificity. The present deep phylogenomics of diverse Xanthomonas species and its members revealed lineage association and dominance of distinct types of LPS in accordance with their origin, host, tissue specificity, and evolutionary success.

Phylo-Taxonogenomics Supports Revision of Taxonomic Status of 20 Xanthomonas Pathovars to Xanthomonas citri

Based on phylo-taxonogenomics criteria, we present amended descriptions for 20 pathovars to Xanthomonas citri. Incidentally, 18 were first reported from India. Seven out of twenty are classified as X. axonopodis, 12 out of 20 as X. campestris, and one as X. cissicola. In this study, we have generated genome sequence data of four pathovars, and the genomes of the remaining 16 were used from the published data. Comprehensive genome-based phylogenomic and taxonogenomic analyses reveal that all these pathovars belong to X. citri and need to reconcile their taxonomic status. This proposal will aid in systematic studies of a major species and its constitutent members that infect economically important plants.

Pest categorisation of Xanthomonas citri pv. viticola

The EFSA Plant Health Panel performed a pest categorisation of Xanthomonas citri pv. viticola (Nayudu) Dye, a Gram-negative bacterium belonging to the Xanthomonadaceae family. The pathogen is a well-defined taxonomic unit and is the causal agent of the leaf spot and bacterial canker of Vitis vinifera. This bacterium is present in India and Brazil, where it affects table grape cultivation; the same pathogen is able to cause a disease on Azadirachta indica and on some weed species. Reports indicate that the bacterium is present in Thailand as well. The pathogen has never been reported from the EU territory and it is not included in EU Commission Implementing Regulation 2019/2072. The pathogen can be detected on its host plants using direct isolation, serological or PCR-based methods. Its identification is achieved using biochemical and nutritional assays, together with a multilocus sequence analysis based on seven housekeeping genes. The main pathway for the entry of the pathogen into the EU territory is plant propagation material. In the EU, there is large availability of host plants, with grapevine being one of the most important crops in Europe and more specifically in its Mediterranean areas. Since X. citri pv. viticola is only reported in tropical and subtropical areas (BSh and Aw climatic zones according to the Köppen-Geiger classification), there is uncertainty whether the climatic conditions in the EU territory are suitable for its establishment. Nevertheless, due to the great importance of grapevine for the EU agriculture, any disease outbreak may have a high-economic impact. Phytosanitary measures are available to prevent the introduction of the pathogen into the EU. X. citri pv. viticola satisfies the criteria that are within the remit of EFSA to assess for this species to be regarded as a potential Union quarantine pest.

Deep phylo-taxono genomics reveals Xylella as a variant lineage of plant associated Xanthomonas and supports their taxonomic reunification along with Stenotrophomonas and Pseudoxanthomonas

Genus Xanthomonas is a group of phytopathogens that is phylogenetically related to Xylella, Stenotrophomonas, and Pseudoxanthomonas, having diverse lifestyles. Xylella is a lethal plant pathogen with a highly reduced genome, atypical GC content and is taxonomically related to these three genera. Deep phylo-taxono genomics reveals that Xylella is a variant Xanthomonas lineage that is sandwiched between Xanthomonas clades. Comparative studies suggest the role of unique pigment and exopolysaccharide gene clusters in the emergence of Xanthomonas and Xylella clades. Pan-genome analysis identified a set of unique genes associated with sub-lineages representing plant-associated Xanthomonas clade and nosocomial origin Stenotrophomonas clade. Overall, our study reveals the importance of reconciling classical phenotypic data and genomic findings in reconstituting the taxonomic status of these four genera. SIGNIFICANCE STATEMENT: Xylella fastidiosa is a devastating pathogen of perennial dicots such as grapes, citrus, coffee, and olives. An insect vector transmits the pathogen to its specific host wherein the infection leads to complete wilting of the plants. The genome of X. fastidiosa is significantly reduced both in terms of size (2 Mb) and GC content (50%) when compared with its relatives such as Xanthomonas, Stenotrophomonas, and Pseudoxanthomonas that have higher GC content (65%) and larger genomes (5 Mb). In this study, using systematic and in-depth genome-based taxonomic and phylogenetic criteria and comparative studies, we assert the need to unify Xanthomonas with its relatives (Xylella, Stenotrophomonas and Pseudoxanthomonas). Interestingly, Xylella revealed itself as a minor variant lineage embedded within two major Xanthomonas lineages comprising member species of different hosts.

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.

Genome Sequencing and Functional Characterization of Xanthomonas cucurbitae, the Causal Agent of Bacterial Spot Disease of Cucurbits

Bacterial leaf spot disease caused by Xanthomonas cucurbitae has severely affected the pumpkin industries in the Midwestern region of United States, with the bacteria mainly infecting pumpkin leaves and fruits, and leading to significant yield losses. In this study, we utilized genomics and genetics approaches to elucidate X. cucurbitae molecular mechanisms of pathogenesis during interaction with its host. We generated the first reference-quality whole-genome sequence of the X. cucurbitae type isolate and compared with other Xanthomonas species, X. cucurbitae has a smaller genome size with fewer virulence-related genes. RNA-seq analysis of X. cucurbitae under plant-mimicking media conditions showed altered transcriptional responses, with upregulation of virulence genes and downregulation of cellular homeostasis genes. Additionally, characterization of key virulence genes using gene deletion methods revealed that both type II enzymes and type III effectors are necessary for X. cucurbitae to cause infection in the pumpkin host.

Mining genome traits that determine the different gut colonization potential of Lactobacillus and Bifidobacterium species

Although the beneficial effects of probiotics are likely to be associated with their ability to colonize the gut, little is known about the characteristics of good colonizers. In a systematic analysis of the comparative genomics, we tried to elucidate the genomic contents that account for the distinct host adaptability patterns of Lactobacillus and Bifidobacterium species. The Bifidobacterium species, with species-level phylogenetic structures affected by recombination among strains, broad mucin-foraging activity, and dietary-fibre-degrading ability, represented niche conservatism and tended to be host-adapted. The Lactobacillus species stretched across three lifestyles, namely free-living, nomadic and host-adapted, as characterized by the variations of bacterial occurrence time, guanine-cytosine (GC) content and genome size, evolution event frequency, and the presence of human-adapted bacterial genes. The numbers and activity of host-adapted factors, such as bile salt hydrolase and intestinal tissue-anchored elements, were distinctly distributed among the three lifestyles. The strains of the three lifestyles could be separated with such a collection of colonization-related genomic content (genes, genome size and GC content). Thus, our work provided valuable information for rational selection and gut engraftment prediction of probiotics. Here, we have found many interesting predictive results for bacterial gut fitness, which will be validated in vitro and in vivo.

Xanthomonas axonopodis pv. punicae depends on multiple non-TAL (Xop) T3SS effectors for its coveted growth inside the pomegranate plant through repressing the immune responses during bacterial blight development

Xanthomonas axonopodis pv. punicae (Xap), the bacterial blight pathogen of pomegranate, incurs substantial loss to yield and reduces export quality of this economically important fruit crop. During infection, the bacterium secretes six non-TAL (Xop) effectors into the pomegranate cells through a specialized type three secretion system (T3SS). Previously, we demonstrated the role of two key effectors, XopL and XopN in pathogenesis. Here, we investigate the role of rest effectors (XopC2, XopE1, XopQ and XopZ) on disease development. We generated null mutants for each individual effector and mutant bacterial suspension was infiltrated into pomegranate leaves. Compared to Xap wild, the mutant bacterial growth was reduced by 2.7-11.5 folds. The mutants produced lesser water-soaked lesions when infiltrated on leaves by 1.13-2.21 folds. Among the four effectors, XopC2 contributes highest for in planta bacterial growth and disease development. XopC2 efficiently suppressed the defense responses like callose deposition, reactive oxygen species (ROS) and the activation of immune responsive genes. Being a major contributor, we further characterize XopC2 for its subcellular localization, its protein structure and networking. XopC2 is localized to the plasma membrane of Nicotiana benthamiana like XopL and XopN. XopC2 is a 661 amino acids protein having 15 alpha and 17 beta helix. Our STRING and I-TASSER based analysis hinted that XopC2 interacts with multiple membrane localized plant proteins including transcription regulator of CCR4-NOT family, TTN of maintenance of chromosome family and serine/threonine-protein phosphatase 2A (PP2A) isoform. Based on the interaction it is predicted that XopC2 might involve in diverse functions like nuclear-transcribed mRNA catabolic process, maintenance of chromosome, hormone signaling and protein dephosphorylation activities and thereby suppress the plant immunity. Altogether, our study suggests that Xap largely depends on three non-TAL (Xop) effectors, including XopC2, XopL and XopN, to modulate pomegranate PTI for its unrestricted proliferation during bacterial blight development.

The Variability of Thymol and Carvacrol Contents Reveals the Level of Antibacterial Activity of the Essential Oils from Different Accessions of Oliveria decumbens

Oliveria decumbens (Apiaceae) is an aromatic herb traditionally employed in the Persian medicine for the treatment of infectious and gastrointestinal disorders. In the present study, we analyzed the chemical composition of essential oils obtained from different Iranian populations and evaluated their efficacy on a panel of human pathogens (Staphylococcus aureus and Escherichia coli), probiotic (Bacillus subtilis), and phytopathogens (Clavibacter michiganensis, Curtobacterium flaccumfaciens, Xanthomonas citri, and Agrobacterium tumefaciens). The gas chromatographic-mass spectrometry analysis put in evidence four main volatile constituents such as thymol (20.3–36.4%), carvacrol (18.8–33.1%), γ-terpinene (10.6–25.9%), and p-cymene (9.5–17.3%), though with significant variability from an essential oil to another. Notably, the oils from the populations sited in Nourabad Mamasani and Dehdasht showed the highest amount of the phenolic monoterpenes thymol (36.4 and 35.2%, respectively) and carvacrol (33.1 and 30.6%, respectively). The antibacterial activity of O. decumbens essential oils was assessed by the minimum inhibitory concentration (MIC) and minimum bactericidal concentration (MBC) methods, showing high activity for the samples from Nourabad Mamasani and Dehdasht populations exhibiting high level of the above phenolics. The obtained MIC and MBC values (mg/ml) were in the ranges 0.0625–2 mg/ml and 1–16 mg/ml, respectively. Noteworthy, in some cases, the antibacterial activity of O. decumbens essential oils was higher than that of chloramphenicol used as positive control. The average MBCs displayed by the O. decumbens samples showed that C. flaccumfaciens had the highest sensitivity to the essential oils. Based on these results, our work shed light on selected O. decumbens populations deserving proper breeding and cultivation strategies in order to warrantee production of bioactive essential oils to be used at pharmaceutical and agricultural level to combat several pathogens.

Phylogenomic Insights into Diversity and Evolution of Nonpathogenic Xanthomonas Strains Associated with Citrus

Xanthomonas species are primarily known as a group of phytopathogens infecting diverse plants. Recent molecular studies reveal the existence of potential novel species and strains of Xanthomonas following a nonpathogenic lifestyle. In the present study, we report whole-genome sequences of four nonpathogenic strains from citrus (NPXc). Taxonogenomics revealed the surprising diversity, as each of these three isolates were found to be potential novel species that together form a citrus-associated nonpathogenic Xanthomonas species complex (NPXc complex). Interestingly, this NPXc complex is related to another nonpathogenic species, Xanthomonas sontii, from rice (NPXr). On the other hand, the fourth NPXc isolate was found to be related to nonpathogenic isolates from walnut (NPXw); altogether, they form a potential taxonomic outlier of pathogenic Xanthomonas arboricola species. Furthermore, genomic investigation of well-characterized pathogenicity clusters in NPXc isolates revealed lifestyle-specific gene content dynamics. Primarily, genes essential for virulence (i.e., type 1 secretion system [T1SS], T2SS and its effectors, T3SS and its effectors, T4SS, T6SS, adhesins, and rpf gene cluster) and adaptation (i.e., gum, iron uptake and utilization, xanthomonadin, and two-component systems) were depicted by comparative genomics of a Xanthomonas community comprising diverse lifestyles. Overall, the present analysis confers that nonpathogenic isolates of diverse hosts phylogenomically converge and are evolving in parallel with their pathogenic counterparts. Hence, there is a need to understand the world of nonpathogenic isolates from diverse and economically important hosts. Genomic knowledge and resources of nonpathogenic strains will be invaluable in both basic and applied research of the genus XanthomonasIMPORTANCEXanthomonas citri is one of the top phytopathogenic bacteria and is the causal agent of citrus canker. Interestingly, Xanthomonas is also reported to be associated with healthy citrus plants. The advent of the genomic era enabled us to carry out a detailed evolutionary study of a Xanthomonas community associated with citrus and other plants. Our genome-based investigations have revealed hidden and extreme interstrain diversity of nonpathogenic Xanthomonas strains from citrus plants, warranting further large-scale studies. This indicates an unexplored world of Xanthomonas from healthy citrus plant species that may be coevolving as a species complex with the host, unlike the variant pathogenic species. The knowledge and genomic resources will be valuable in evolutionary studies exploring its hidden potential and management of pathogenic species.

Extensive Genomic Rearrangements along with Distinct Mobilome and TALome are Associated with Extreme Pathotypes of a Rice Pathogen

Xanthomonas oryzae pv. oryzae (Xoo) is a serious pathogen of rice which displays tremendous interstrain variation. The emergence of highly-virulent strains of Xoo is a major threat to rice cultivation. Evolutionary insights into genome dynamics of highly virulent strains as compared with the less-virulent ones are crucial for understanding the molecular basis of exceptional success of Xoo as a highly evolved plant pathogen. In the present study, we report complete genome sequence of Xoo strains with extreme-virulent pathotypes (XVPs) characterized based on their reaction toward ten resistance (Xa) genes. One strain, IXO1088, can overcome resistance mediated by all the ten resistance genes while the other strain IXO704 cannot overcome any of them. Interestingly, our investigation revealed that XVPs display dramatic variation in the genome structure with numerous rearrangements/inversions. Moreover, XVPs also possess distinct transposon content and prophage elements that may provide genomic flux required for the acquisition of novel gene cassettes and structural changes in the genome. Interestingly, analysis of transcription activator-like effector proteins, which are major virulence determinants of Xanthomonas pathogen show marked variation in the transcription activator-like effector content and DNA binding domain of tal genes. Overall, the present study indicates the possible role of mobilomes and repetitive elements in major structural and sequence alterations, which may be leading to the emergence of novel and extreme pathotypes. The knowledge and resource of XVPs will be invaluable in the further systematic understanding of evolution and management of variant pathotypes of Xoo.

Complete Genome Sequence Reveals Evolutionary and Comparative Genomic Features of Xanthomonas albilineans Causing Sugarcane Leaf Scald

Leaf scald (caused by Xanthomonas albilineans) is an important bacterial disease affecting sugarcane in most sugarcane growing countries, including China. High genetic diversity exists among strains of X. albilineans from diverse geographic regions. To highlight the genomic features associated with X. albilineans from China, we sequenced the complete genome of a representative strain (Xa-FJ1) of this pathogen using the PacBio and Illumina platforms. The complete genome of strain Xa-FJ1 consists of a circular chromosome of 3,724,581 bp and a plasmid of 31,536 bp. Average nucleotide identity analysis revealed that Xa-FJ1 was closest to five strains from the French West Indies and the USA, particularly to the strain GPE PC73 from Guadeloupe. Comparative genomic analysis between Xa-FJ1 and GPE PC73 revealed prophage integration, homologous recombination, transposable elements, and a clustered regulatory interspaced short palindromic repeats (CRISPR) system that were linked with 16 insertions/deletions (InDels). Ten and 82 specific genes were found in Xa-FJ1 and GPE PC73, respectively, and some of these genes were subjected to phage-related proteins, zona occludens toxin, and DNA methyltransferases. Our findings highlight intra-species genetic variability of the leaf scald pathogen and provide additional genomic resources to investigate its fitness and virulence.

CRISPR elements provide a new framework for the genealogy of the citrus canker pathogen Xanthomonas citri pv. citri

Background

Xanthomonads are an important clade of Gram-negative bacteria infecting a plethora of economically important host plants, including citrus. Knowledge about the pathogen’s diversity and population structure are prerequisite for epidemiological surveillance and efficient disease management. Rapidly evolving genetic loci, such as Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), are of special interest to develop new molecular typing tools.

Results

We analyzed CRISPR loci of 56 Xanthomonas citri pv. citri strains of world-wide origin, a regulated pathogen causing Asiatic citrus canker in several regions of the world. With one exception, 23 unique sequences built up the repertoire of spacers, suggesting that this set of strains originated from a common ancestor that already harbored these 23 spacers. One isolate originating from Pakistan contained a string of 14 additional, probably more recently acquired spacers indicating that this genetic lineage has or had until recently the capacity to acquire new spacers. Comparison of CRISPR arrays with previously obtained molecular typing data, such as amplified fragment length polymorphisms (AFLP), variable-number of tandem-repeats (VNTR) and genome-wide single-nucleotide polymorphisms (SNP), demonstrated that these methods reveal similar evolutionary trajectories. Notably, genome analyses allowed to generate a model for CRISPR array evolution in X. citri pv. citri, which provides a new framework for the genealogy of the citrus canker pathogen.

Conclusions

CRISPR-based typing will further improve the accuracy of the genetic identification of X. citri pv. citri outbreak strains in molecular epidemiology analyses, especially when used concomitantly with another genotyping method.

The Green Revolution shaped the population structure of the rice pathogen Xanthomonas oryzae pv. oryzae

The impact of modern agriculture on the evolutionary trajectory of plant pathogens is a central question for crop sustainability. The Green Revolution replaced traditional rice landraces with high-yielding varieties, creating a uniform selection pressure that allows measuring the effect of such intervention. In this study, we analyzed a unique historical pathogen record to assess the impact of a major resistance gene, Xa4, in the population structure of Xanthomonas oryzae pv. oryzae (Xoo) collected in the Philippines in a span of 40 years. After the deployment of Xa4 in the early 1960s, the emergence of virulent pathogen groups was associated with the increasing adoption of rice varieties carrying Xa4, which reached 80% of the total planted area. Whole genomes analysis of a representative sample suggested six major pathogen groups with distinctive signatures of selection in genes related to secretion system, cell-wall degradation, lipopolysaccharide production, and detoxification of host defense components. Association genetics also suggested that each population might evolve different mechanisms to adapt to Xa4. Interestingly, we found evidence of strong selective sweep affecting several populations in the mid-1980s, suggesting a major bottleneck that coincides with the peak of Xa4 deployment in the archipelago. Our study highlights how modern agricultural practices facilitate the adaptation of pathogens to overcome the effects of standard crop improvement efforts.

Origin and diversification of Xanthomonas citri subsp. citri pathotypes revealed by inclusive phylogenomic, dating, and biogeographic analyses

BACKGROUND

Xanthomonas citri subsp. citri pathotypes cause bacterial citrus canker, being responsible for severe agricultural losses worldwide. The A pathotype has a broad host spectrum, while A* and Aw are more restricted both in hosts and in geography. Two previous phylogenomic studies led to contrasting well-supported clades for sequenced genomes of these pathotypes. No extensive biogeographical or divergence dating analytic approaches have been so far applied to available genomes.

RESULTS

Based on a larger sampling of genomes than in previous studies (including six new genomes sequenced by our group, adding to a total of 95 genomes), phylogenomic analyses resulted in different resolutions, though overall indicating that A + AW is the most likely true clade. Our results suggest the high degree of recombination at some branches and the fast diversification of lineages are probable causes for this phylogenetic blurring effect. One of the genomes analyzed, X. campestris pv. durantae, was shown to be an A* strain; this strain has been reported to infect a plant of the family Verbenaceae, though there are no reports of any X. citri subsp. citri pathotypes infecting any plant outside the Citrus genus. Host reconstruction indicated the pathotype ancestor likely had plant hosts in the family Fabaceae, implying an ancient jump to the current Rutaceae hosts. Extensive dating analyses indicated that the origin of X. citri subsp. citri occurred more recently than the main phylogenetic splits of Citrus plants, suggesting dispersion rather than host-directed vicariance as the main driver of geographic expansion. An analysis of 120 pathogenic-related genes revealed pathotype-associated patterns of presence/absence.

CONCLUSIONS

Our results provide novel insights into the evolutionary history of X. citri subsp. citri as well as a sound phylogenetic foundation for future evolutionary and genomic studies of its pathotypes.

Variation and inheritance of the Xanthomonas raxX-raxSTAB gene cluster required for activation of XA21-mediated immunity

The rice XA21-mediated immune response is activated on recognition of the RaxX peptide produced by the bacterium Xanthomonas oryzae pv. oryzae (Xoo). The 60-residue RaxX precursor is post-translationally modified to form a sulfated tyrosine peptide that shares sequence and functional similarity with the plant sulfated tyrosine (PSY) peptide hormones. The 5-kb raxX-raxSTAB gene cluster of Xoo encodes RaxX, the RaxST tyrosylprotein sulfotransferase, and the RaxA and RaxB components of a predicted type I secretion system. To assess raxX-raxSTAB gene cluster evolution and to determine its phylogenetic distribution, we first identified rax gene homologues in other genomes. We detected the complete raxX-raxSTAB gene cluster only in Xanthomonas spp., in five distinct lineages in addition to X. oryzae. The phylogenetic distribution of the raxX-raxSTAB gene cluster is consistent with the occurrence of multiple lateral (horizontal) gene transfer events during Xanthomonas speciation. RaxX natural variants contain a restricted set of missense substitutions, as expected if selection acts to maintain peptide hormone-like function. Indeed, eight RaxX variants tested all failed to activate the XA21-mediated immune response, yet retained peptide hormone activity. Together, these observations support the hypothesis that the XA21 receptor evolved specifically to recognize Xoo RaxX.

Genomic sequence analysis reveals diversity of Australian Xanthomonas species associated with bacterial leaf spot of tomato, capsicum and chilli

Background

The genetic diversity in Australian populations of Xanthomonas species associated with bacterial leaf spot in tomato, capsicum and chilli were compared to worldwide bacterial populations. The aim of this study was to confirm the identities of these Australian Xanthomonas species and classify them in comparison to overseas isolates. Analysis of whole genome sequence allows for the investigation of bacterial population structure, pathogenicity and gene exchange, resulting in better management strategies and biosecurity.

Results

Phylogenetic analysis of the core genome alignments and SNP data grouped strains in distinct clades. Patterns observed in average nucleotide identity, pan genome structure, effector and carbohydrate active enzyme profiles reflected the whole genome phylogeny and highlight taxonomic issues in X. perforans and X. euvesicatoria. Circular sequences with similarity to previously characterised plasmids were identified, and plasmids of similar sizes were isolated. Potential false positive and false negative plasmid assemblies were discussed. Effector patterns that may influence virulence on host plant species were analysed in pathogenic and non-pathogenic xanthomonads.

Conclusions

The phylogeny presented here confirmed X. vesicatoria, X. arboricola, X. euvesicatoria and X. perforans and a clade of an uncharacterised Xanthomonas species shown to be genetically distinct from all other strains of this study. The taxonomic status of X. perforans and X. euvesicatoria as one species is discussed in relation to whole genome phylogeny and phenotypic traits. The patterns evident in enzyme and plasmid profiles indicate worldwide exchange of genetic material with the potential to introduce new virulence elements into local bacterial populations.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5600-x) contains supplementary material, which is available to authorized users.

Xanthomonas citri pv. viticola Affecting Grapevine in Brazil: Emergence of a Successful Monomorphic Pathogen

The pathovar viticola of Xanthomonas citri causes bacterial canker of grapevine. This disease was first recorded in India in 1972, and later in Brazil in 1998, where its distribution is currently restricted to the northeastern region. A multilocus sequence analysis (MLSA) based on seven housekeeping genes and a multilocus variable number of tandem repeat analysis (MLVA) with eight loci were performed in order to assess the genetic relatedness among strains from India and Brazil. Strains isolated in India from three related pathovars affecting Vitaceae species and pathogenic strains isolated from Amaranthus sp. found in bacterial canker-infected vineyards in Brazil were also included. MLSA revealed lack of diversity in all seven genes and grouped grapevine and Amaranthus strains in a monophyletic group in X. citri. The VNTR (variable number of tandem repeat) typing scheme conducted on 107 strains detected 101 haplotypes. The total number of alleles per locus ranged from 5 to 12. A minimum spanning tree (MST) showed that Brazilian strains were clearly separated from Indian strains, which showed unique alleles at three loci. The two strains isolated from symptomatic Amaranthus sp. presented unique alleles at two loci. STRUCTURE analyses revealed three groups congruent with MST and a fourth group with strains from India and Brazil. Admixture among populations were observed in all groups. MST, STRUCTURE and e-BURST analyses showed that the strains collected in 1998 belong to two distinct groups, with predicted founder genotypes from two different vineyards in the same region. This suggest that one introduction of grape planting materials contaminated with genetically distinct strains took place, which was followed by pathogen adaptation. Genome sequencing of one Brazilian strain confirmed typical attributes of pathogenic xanthomonads and allowed the design of a complementary VNTR typing scheme dedicated to X. citri pv. viticola that will allow further epidemiological survey of this genetically monomorphic pathovar.

A highly specific tool for identification of Xanthomonas vasicola pv. musacearum based on five Xvm-specific coding sequences

Xanthomonas vasicola pv. musacearum (Xvm) is a bacterial pathogen responsible for the economically important Xanthomonas wilt disease on banana and enset crops in Sub-Saharan Africa. Given that the symptoms are similar to those of other diseases, molecular diagnosis is essential to unambiguously identify this pathogen and distinguish it from closely related strains not pathogenic on these hosts. Currently, Xvm identification is based on polymerase chain reaction (PCR) with GspDm primers, targeting the gene encoding general secretory protein D. Experimental results and examination of genomic sequences revealed poor specificity of the GspDm PCR. Here, we present and validate five new Xvm-specific primers amplifying only Xvm strains.

Complete Genome Sequence Reveals Evolutionary Dynamics of an Emerging and Variant Pathovar of Xanthomonas euvesicatoria

Xanthomonas, a complex group of pathogens, infects more than 400 plants, which is expanding to new hosts causing serious diseases. Genome-based studies are transforming our understanding on diversity and relationship of host-specific members, known as pathovars. In this study, we report complete genome sequence of a novel pathovar Xanthomonas axonopodis pv. commiphorae (Xcom) from India. It causes gumming disease of Commiphora wightii, a medicinally important plant. Genome-based phylogenetic and taxonomic investigation revealed that the pathovar belongs to Xanthomonas euvesicatoria and not X. axonopodis as reported earlier. Interestingly, it is a novel host and novel geographic origin for a X. euvesicatoria pathovar. A core-genome-based phylogenetic analysis resolved the pathovar complex of this species on the basis of their hosts. Interestingly, this pathovar harbors a unique 35-kb plasmid encoding type III effectors and toxin-antitoxin gene that is absent in other X. euvesicatoria pathovars and infects tomato, pepper, rose, onion, philodendron, alfalfa, and citrus plants. The pathovar contains two TAL (transcription activator-like) genes, one on plasmid and another on genomic region with an additional pseudo TAL gene flanked by IS elements in the plasmid. Further, Xcom has acquired a novel set of lipopolysaccharide biosynthesis genes after its divergence from the closely related pathovar that infects rose and supports the role of horizontal gene transfer in hypervariation at this locus in the species. Complete genome sequence of this variant pathovar has provided novel insights into evolution of an emerging pathovar in Xanthomonas and will be valuable resource in pathogenomics of X. euvesicatoria.

Identification of an apiosyltransferase in the plant pathogen Xanthomonas pisi

The rare branched-chain sugar apiose, once thought to only be present in the plant kingdom, was found in two bacterial species: Geminicoccus roseus and Xanthomonas pisi. Glycans with apiose residues were detected in aqueous methanol-soluble fractions as well as in the insoluble pellet fraction of X. pisi. Genes encoding bacterial uridine diphosphate apiose (UDP-apiose) synthases (bUASs) were characterized in these bacterial species, but the enzyme(s) involved in the incorporation of the apiose into glycans remained unknown. In the X. pisi genome two genes flanking the XpUAS were annotated as hypothetical glycosyltransferase (GT) proteins. The first GT (here on named XpApiT) belongs to GT family 90 and has a Leloir type B fold and a putative lipopolysaccharide-modifying (LPS) domain. The second GT (here on XpXylT) belongs to GT family 2 and has a type A fold. The XpXylT and XpApiT genes were cloned and heterologously expressed in E. coli. Analysis of nucleotide sugar extracts from E. coli expressing XpXylT or XpApiT with UAS showed that recombinant XpApiT utilized UDP-apiose and XpXylT utilized UDP-xylose as substrate. Indirect activity assay (UDP-Glo) revealed that XpApiT is an apiosyltransferase (ApiT) able to specifically use UDP-apiose. Further support for the apiosyltransferase activity was demonstrated by in microbe co-expression of UAS and XpApiT in E. coli showing the utilization of UDP-apiose to generate an apioside detectable in the pellet fraction. This work provides evidence that X. pisi developed the ability to synthesize an apioside of indeterminate function; however, the evolution of the bacterial ApiT remains to be determined. From genetic and evolutionary perspectives, the apiose operon may provide a unique opportunity to examine how genomic changes reflect ecological adaptation during the divergence of a bacterial group.

Taxonomic Repositioning of Xanthomonas campestris pv. viticola (Nayudu 1972) Dye 1978 as Xanthomonas citri pv. viticola (Nayudu 1972) Dye 1978 comb. nov. and Emendation of the Description of Xanthomonas citri pv. anacardii to Include Pigmented Isolates Pathogenic to Cashew Plant

Grapevine bacterial canker, which is caused by Xanthomonas campestris pv. viticola, is one of the most important grapevine diseases in the northeastern region of Brazil. This disease causes severe damage and represents a high potential risk to the development of Brazilian viticulture. In turn, pigmented isolates pathogenic to cashew plant, making cashew fruit unfit for sale, also have been detected in Northeastern Brazil. Given that the taxonomic position of these bacteria is unclear, the multilocus sequence analysis (MLSA) technique, average nucleotide identity (ANI) values and tetranucleotide frequency correlation coefficients (TETRA) were used to analyze their phylogenetic relationship in relation to other Xanthomonas species. X. campestris pv. viticola was closely related to X. citri pv. mangiferaeindicae (repetitive-polymerase chain reaction [rep-PCR], MLSA, and ANI) and X. citri subsp. citri (MLSA and ANI). Pigmented isolates pathogenic to cashew plant were closely related to X. citri pv. anacardii (rep-PCR, MLSA, ANI, and TETRA). The results obtained in this study support the emendation of the description of X. citri pv. anacardii to include pigmented isolates of Xanthomonas pathogenic to cashew plant. In addition, the reclassification of X. campestris pv. viticola as X. citri pv. viticola comb. nov. is suggested.