# Complete Genome Assemblies of All Xanthomonas translucens Pathotype Strains Reveal Three Genetically Distinct Clades

Florian Goettelmann, Veronica Roman-Reyna, Sébastien Cunnac, Jonathan Jacobs, Claude Bragard, Bruno Studer, Ralf Koebnik, Roland Kölliker

<figure><img src="image_1.png" /></figure>

## To cite this version:

Florian Goettelmann, Veronica Roman-Reyna, Sébastien Cunnac, Jonathan Jacobs, Claude Bragard, et al.. Complete Genome Assemblies of All Xanthomonas translucens Pathotype Strains Reveal Three Genetically Distinct Clades. Frontiers in Microbiology, 2022, 12, pp.817815. 10.3389/fmicb.2021.817815. [https://hal.science/hal-03626316](https://hal.science/hal-03626316)

HAL Id: hal-03626316

https://hal.science/hal-03626316

Submitted on 12 Apr 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for
updates

# Complete Genome Assemblies of All *Xanthomonas translucens* Pathotype Strains Reveal Three Genetically Distinct Clades

Florian Goettelmann¹, Veronica Roman-Reyna²,³, Sébastien Cunnac⁴, Jonathan M. Jacobs²,³, Claude Bragard⁵, Bruno Studer¹, Ralf Koebnik⁴ and Roland Kölliker¹*

¹ Molecular Plant Breeding, Institute of Agricultural Sciences, ETH Zürich, Zurich, Switzerland, ² Department of Plant Pathology, The Ohio State University, Columbus, OH, United States, ³ Infectious Diseases Institute, The Ohio State University, Columbus, OH, United States, ⁴ Plant Health Institute of Montpellier, University of Montpellier, CIRAD, INRAE, Institut Agro, IRD, Montpellier, France, ⁵ Earth and Life Institute, UCLouvain, Louvain-la-Neuve, Belgium

## OPEN ACCESS

**Edited by:**
Feng Gao,
Tianjin University, China

**Reviewed by:**
Meng Yuan,
Huazhong Agricultural University,
China

Toni A. Chapman,
New South Wales Department
of Primary Industries, Australia

*Correspondence:
Roland Kölliker
roland.koelliker@usys.ethz.ch

**Specialty section:**
This article was submitted to Evolutionary and Genomic Microbiology, a section of the journal Frontiers in Microbiology

Received: 18 November 2021

Accepted: 29 December 2021

Published: 02 March 2022

**Citation:**
Goettelmann F, Roman-Reyna V, Cunnac S, Jacobs JM, Bragard C, Studer B, Koebnik R and Kölliker R (2022) Complete Genome Assemblies of All Xanthomonas translucens Pathotype Strains Reveal Three Genetically Distinct Clades. *Front. Microbiol.* 12:817815.
doi: 10.3389/fmicb.2021.817815

The *Xanthomonas translucens* species comprises phytopathogenic bacteria that can cause serious damage to cereals and to forage grasses. So far, the genomic resources for *X. translucens* were limited, which hindered further understanding of the host-pathogen interactions at the molecular level and the development of disease-resistant cultivars. To this end, we complemented the available complete genome sequence of the *X. translucens* pv. *translucens* pathotype strain DSM 18974 by sequencing the genomes of all the other 10 *X. translucens* pathotype strains using PacBio long-read technology and assembled complete genome sequences. Phylogeny based on average nucleotide identity (ANI) revealed three distinct clades within the species, which we propose to classify as clades Xt-I, Xt-II, and Xt-III. In addition to 2,181 core *X. translucens* genes, a total of 190, 588, and 168 genes were found to be exclusive to each clade, respectively. Moreover, 29 non-transcription activator-like effector (TALE) and 21 TALE type III effector classes were found, and clade- or strain-specific effectors were identified. Further investigation of these genes could help to identify genes that are critically involved in pathogenicity and/or host adaptation, setting the grounds for the development of new resistant cultivars.

**Keywords:** Xanthomonas translucens, complete genomes, phylogeny, host adaptation, comparative genomics, virulence factors

## INTRODUCTION

*Xanthomonas translucens* is a species of Gram-negative phytopathogenic bacteria that causes serious damage in gramineous plants. To date, a total of 11 pathovars have been defined based on their host range, each represented by one strain referred to as the pathotype strain and deposited in the appropriate bacterial strain collections. Two major groups are historically distinguished within the species: the “translucens” group and the “graminis” group. The translucens group consists of pathovars *cerealis*, *hordei*, *secalis*, *translucens*, and *undulosa*, which cause leaf streak and black chaff in economically important cereals such as wheat (*Triticum* spp.) and barley
(Hordeum vulgare) (Sapkota et al., 2020). The graminis group consists of pathovars arrhenatheri, graminis, phlei, phleipratensis, and poae, which cause bacterial wilt in forage grasses (Egli et al., 1975; Egli and Schmidt, 1982; Stead, 1989). Although recent genome sequence data loosely supported this phenotypic classification (Peng et al., 2016; Hersemann et al., 2017; Langlois et al., 2017; Shah et al., 2021), these groups are not clearly defined based on the genetic relationships between the strains. Moreover, genomic studies have previously suggested that *pv. cerealis* could be genetically distinct from the other “translucens” group pathovars (Bragard et al., 1997; Rademaker et al., 2006; Peng et al., 2016; Langlois et al., 2017). Furthermore, *Xanthomonas* strains causing dieback in pistachio (*Pistacia vera* L.) were found to belong to the *X. translucens* species, a species so far thought to infect only monocotyledonous plants (Giblot-Ducray et al., 2009). These strains were classified as *pv. pistaciae*, in which two groups were distinguished by rep-PCR and *gyrB*-based phylogeny and are referred to as groups A and B (Facelli et al., 2005; Marefat et al., 2006; Giblot-Ducray et al., 2009). The group A strain was found to be closely related to pathovars *translucens*, *secalis*, and *undulosa*, while the group B strain was closest to *pv. cerealis*, suggesting that the ability to infect pistachio may have evolved two separate times in *X. translucens*.

*X. translucens* pathovars were initially classified solely based on their host range as pathovars of the *X. campestris* species. However, a classification based on DNA-DNA hybridization was proposed by Vauterin et al. (1995), who amended these pathovars of *X. campestris* to the species level, forming the *X. translucens* species. Nonetheless, though these pathovars are now considered pathovars of the *X. translucens* species, the definition of these pathovars is still based on their host range as before. Phylogeny is crucial to better understand the evolutionary history of each pathovar, in order to determine how each strain has adapted to infect its hosts. This knowledge is key in defining targets for breeding resistant cultivars. The current definition of pathovars of *X. translucens* might not properly reflect these processes of host adaptation, and a better definition based on the genetic relationships between the strains would improve the comprehension of the host-pathogen relationships. Recent advances in sequencing technology now allow for the rapid sequencing of complete genomes of bacterial strains, which provide a better basis for phylogenetic analyses. At the time of this study, complete genome sequences were publicly available for *pv. cerealis* strain NXtc01 (Shah et al., 2019); *pv. translucens* strains DSM 18974 (Jaenicke et al., 2016), XtKm7, XtKm8, XtKm9, XtKm33, XtKm34 (Shah et al., 2021); and *pv. undulosa* strains XT4699 (Peng et al., 2016), ICMP 11055 (Falahi Charkhabi et al., 2017), LW16, P3 (Peng et al., 2019), XtFa1, XtLr8, XtKm12, and XtKm15 (Shah et al., 2021). Though the number of available complete genome sequences is continuously growing, these sequences represent only three pathovars to date, with a complete genome sequence being available for only one pathotype strain.

High-quality genome sequences also help in the identification of genes that are directly linked to the pathogen's virulence and host range. Indeed, successful infection by pathogenic bacteria is mediated by a set of virulence factors, including degradative

enzymes and effector proteins. These are often host specific, and the repertoire of virulence factors that a bacterial strain possesses defines which hosts it is able to infect. In *Xanthomonas* species, virulence factors generally depend on the type II and type III secretion systems (T2SS and T3SS, respectively). These secretion systems allow the bacteria to export virulence factors in order to enable or facilitate their proliferation and survival in the host by targeting specific host components (Büttner and Bonas, 2010; Alvarez-Martinez et al., 2021). Identifying which virulence factors are involved in the pathogen's virulence and which host components are the targets of these virulence factors is crucial, as these host components can then be the focus of resistance breeding.

The T3SS is encoded by the “hypersensitive reaction and pathogenicity” (*hrp*) gene cluster (Bonas et al., 1991). Some *hrp* genes that are highly conserved are referred to as “*hrp*-conserved” (*hrc*) genes (Bogdanove et al., 1996). Moreover, some genes in the cluster that are involved in, but not necessary to, the host-plant interaction are called “*hrp*-associated” (*hpa*) genes. In the *X. translucens* species, the core *hrp* cluster consists of 23 genes, with 8 *hrp* genes, 11 *hrc* genes, and 4 *hpa* genes (Wichmann et al., 2013; Pesce et al., 2017).

This secretion system injects effector proteins into the host cell. In *Xanthomonas*, these effectors are generally called “*Xanthomonas* outer proteins” (Xop), with 53 classes from XopA to XopBA (White et al., 2009). Other effectors are named according to their avirulence characteristics, causing hypersensitive response in the host, such as AvrBs1 to AvrBs3. These effectors are key virulence factors, as, when translocated into the host cell, they are able to target the different pathways of the host, allowing the pathogen, for example, to acquire nutrients or to evade or suppress host defenses. A specific type of effectors secreted by the T3SS are the “transcription activator-like” effectors (TALEs). Their amino acid sequences contain highly conserved repetitive sequences of ~34 amino acids, with only the 12th and 13th residues being hypervariable and referred to as the “repeat variable di-residue” (RVD). The RVD array of each TALE allows it to bind to specific nucleotide sequences in the host DNA, thus activating the expression of the neighboring genes to the pathogen’s advantage (Streubel et al., 2017; Wang et al., 2017).

Due to their repetitive sequences, TALE genes are difficult to assemble using short-read sequencing technology. Complete genome sequences based on either long-read sequencing, a very high coverage of short-read sequencing, or a mixture of both, are thus necessary to properly identify these effectors. Such high-quality genome sequences allowed the identification of eight and five TALE genes in *pv. translucens* strains DSM 18974 and UPB886, respectively (Jaenicke et al., 2016; Roman-Reyna et al., 2020), eight and seven in *pv. undulosa* strains XT4699 and ICMP 11055, respectively (Peng et al., 2016; Falahi Charkhabi et al., 2017), and two in *pv. cerealis* strains CFBP 2541 and NXtc01 (Pesce et al., 2015; Shah et al., 2019). However, of these, only four were functionally characterized and found to play a role in virulence to date (Falahi Charkhabi et al., 2017; Peng et al., 2019; Shah et al., 2019). No TALE has yet been identified in other pathovars of *X. translucens*.
Similar to the T3SS, the T2SS is responsible for the export of virulence factors, most of which are cell wall-degradative enzymes, into the host apoplasm (Jha et al., 2005). In *Xanthomonas*, two types of clusters encoding the T2SS can be found. The *xps* cluster is conserved across *Xanthomonas* species, while the *xcs* cluster is only found in some species such as *X. citri* and *X. campestris* (Szczesny et al., 2010).

In addition to the T2SS and T3SS, the type IV and type VI secretion systems (T4SS and T6SS, respectively) also secrete proteins that may affect the virulence of the pathogen. However, while the T2SS and the T3SS target host components, the T4SS and the T6SS are involved in the defense against microbial predators such as amoeba, as well as in the competition with other microorganisms (Büttner and Bonas, 2010; Alvarez-Martinez et al., 2021). The T4SS is evolutionarily related to bacterial conjugation systems and is involved in the competition with other bacteria by injecting protein effectors or protein-DNA complexes into their cells (Sgro et al., 2019). The T6SS is related to the tail of bacteriophages and, similarly to the T4SS, is able to inject effector proteins into prokaryotic as well as eukaryotic cells (Bayer-Santos et al., 2019). In the Xanthomonadales order, three subtypes of T6SS have been found: subtypes 1, 3, and 4. Moreover, subtype 3 is further subdivided into subgroups 3*, 3**, and 3***. Although not directly related to virulence, the T4SS and T6SS could be key elements in the proliferation and survival of the bacteria on the host plant (Souza et al., 2015; Choi et al., 2020).

In this study, in order to complement the available complete genome sequence of the pv. *translucens* pathotype strain DSM 18974 (Jaenicke et al., 2016), we sequenced the whole genomes of all the other 10 pathotype strains of the *X. translucens* species, as well as a representative strain of the *X. translucens* pv. *pistaciae* group B, to produce high-quality genome sequences. Using these, we built a phylogeny of these strains to clarify the taxonomy of the *X. translucens* species. We then scrutinized the genomes for major virulence features of these strains to identify genes that might be important for pathogenicity and in defining their host range.

# MATERIALS AND METHODS

## Bacterial Strains, Growth Conditions, and DNA Extraction

The relevant data for all *X. translucens* strains used in this study are listed in Table 1. The genome sequence of *X. translucens* pv. *translucens* strain DSM 18974 was retrieved from the National Center for Biotechnology Information (NCBI) GenBank database (accession number LT604072). Strains LMG 726, LMG 727, LMG 728, LMG 730, LMG 843, and UPB458 were grown at 28°C on YDC agar medium (2% dextrose, 1% yeast extract, 2% CaCO3, 1.5% agar) for 48 h. Bacteria were then dissolved in 10 ml washing buffer (50 mM TRIS-HCl pH 8.0, 50 mM Ethylenediaminetetraacetic acid (EDTA) pH 8.0, 150 mM NaCl). The genomic DNA was then extracted with the NucleoSpin® Microbial DNA kit (Macherey Nagel, Duren, Germany), according to the manufacturer's recommendations. Strains CFBP 2055, CFBP 2539, CFBP 2541, and CFBP 8304 were grown at 28°C on PSA medium (0.5% peptone, 2% sucrose,

1.5% agar) for 24 h. Bacteria were then resuspended in 10 mM MgCl₂ and diluted to an optical density at 600 nm of 1.0. Cells from 2 ml were harvested by centrifugation and washed once with 10 mM MgCl₂, and genomic DNA was isolated using QIAGEN Genomic tip 100/G (QIAGEN, Hilden, Germany) according to the manufacturer's instructions. The genomic DNA from strain ICMP 16317 was extracted following a standard phenol/chloroform method (Booher et al., 2015).

## Sequencing, Genome Assembly, and Annotation

The library preparation and DNA sequencing of strains LMG 726, LMG 727, LMG 728, LMG 730, LMG 843, and UPB458 were done at the Functional Genomics Center Zurich. For these strains, as well as for ICMP 16317, libraries were prepared and multiplexed with the PacBio SMRTbell® Express Template Prep Kit 2.0 (PacBio, Menlo Park, CA, United States) according to the published protocol<fnref n="1" />. Multiplex library preparation in pools of eight strains, including strains CFBP 2055, CFBP 2539, CFBP 2541, and CFBP 8304, and simultaneous sequencing of eight strains on one SMRTCell was conducted at the GENTYANE genotyping platform (INRA Clermont-Ferrand, France). All strains were sequenced with the PacBio Sequel technology.

Genomic sequences were then de novo assembled with Flye 2.7 for strains CFBP 2055, CFBP 2539, CFBP 2541, and CFBP 8304; Flye 2.8.1 for LMG 726, LMG 728, LMG 730, and UPB458; Flye 2.9 for LMG 727, and LMG 843; and HGAP 4 for ICMP 16317 (Chin et al., 2013; Kolmogorov et al., 2019). Assemblies produced with Flye were done using the “--plasmids --iterations 2” parameters. Genomes were functionally annotated using the NCBI Prokaryotic Genome Annotation Pipeline (PGAP) or with Prokka 1.13 (Seemann, 2014).

## Phylogeny and Comparative Analysis

Average nucleotide identity (ANI) was calculated using Pyani 0.2.11 with default parameters (Pritchard et al., 2016), and a phylogeny dendrogram was constructed using Ward's hierarchical clustering method in R version 4.1.0. The gene content of the strains used in this study was compared with Roary 3.7.0 using the Prokka annotation with default parameters (Page et al., 2015). To identify genomic rearrangements and conserved genomic regions within the three clades that we defined, we compared the genomic structure of the complete genome sequences used in this study with Mauve v20150226 (Darling, 2004).

The presence of T2SS, T3SS, and T4SS was determined by tBLASTn using the amino acid sequences of the main components of each cluster retrieved from UniProt as query (Supplementary Dataset 1) and the genomic sequences of each strain as subject, with a 0.01 e-value threshold. The presence of each secretion system was then validated by looking for the presence of the cluster in the PGAP annotation at the predicted locus. The presence of T6SS was determined using SecReT6 3.0

<footnote n="1">¹https://www.pacb.com/wp-content/uploads/Procedure-Checklist-Preparing-Multiplexed-Microbial-Libraries-Using-SMRTbell-Express-Template-Prep-Kit-2.0.pdf</footnote>
<table><caption>**TABLE 1** | X. translucens pathotype strains used in this study and characteristics of the obtained genome assemblies.</caption><thead><tr><th>Pathovar</th><th>Strain</th><th>Country</th><th>Isolated from</th><th>Origin<sup>a</sup></th><th>Coverage</th><th>Assembly length</th><th>Plasmid length</th><th>Accession</th></tr></thead><tbody><tr><td>arrhenatheri</td><td>LMG 727</td><td>Switzerland</td><td>Arrhenatherum elatius</td><td>BCCM</td><td>103</td><td>4,843,101</td><td>NA</td><td>CP086333</td></tr><tr><td>cerealis</td><td>CFBP 2541</td><td>United States</td><td>Bromus inermis</td><td>CFBP</td><td>404</td><td>4,504,942</td><td>NA</td><td>CP074364</td></tr><tr><td>graminis</td><td>LMG 726</td><td>Switzerland</td><td>Dactylis glomerata</td><td>BCCM</td><td>81</td><td>4,678,781</td><td>NA</td><td>CP076254</td></tr><tr><td>hordei</td><td>UPB458<sup>b</sup></td><td>India</td><td>Hordeum vulgare</td><td>BCCM</td><td>94</td><td>4,679,124</td><td>NA</td><td>CP076249</td></tr><tr><td>phlei</td><td>LMG 730</td><td>Norway</td><td>Phleum pratense</td><td>BCCM</td><td>87</td><td>4,569,024</td><td>NA</td><td>CP076251</td></tr><tr><td>phleipratensis</td><td>LMG 843</td><td>United States</td><td>Phleum pratense</td><td>BCCM</td><td>121</td><td>4,902,099</td><td>NA</td><td>CP086332</td></tr><tr><td>pistaciae (group A)</td><td>CFBP 8304</td><td>Australia</td><td>Pistacia vera</td><td>CFBP</td><td>301</td><td>4,599,174</td><td>NA</td><td>CP074365</td></tr><tr><td>pistaciae (group B)</td><td>ICMP 16317</td><td>Australia</td><td>Pistacia vera</td><td>ICMP</td><td>161</td><td>4,386,175</td><td>NA</td><td>CP083804</td></tr><tr><td>poae</td><td>LMG 728</td><td>Switzerland</td><td>Poa trivialis</td><td>BCCM</td><td>86</td><td>4,792,655</td><td>NA</td><td>CP076250</td></tr><tr><td>secalis</td><td>CFBP 2539</td><td>Canada</td><td>Secale cereale</td><td>CFBP</td><td>186</td><td>4,565,955</td><td>NA</td><td>CP074363</td></tr><tr><td>translucens</td><td>DSM 18947<sup>c</sup></td><td>United States</td><td>Hordeum vulgare</td><td>DSMZ</td><td>223</td><td>4,715,357</td><td>NA</td><td>LT604072</td></tr><tr><td>undulosa</td><td>CFBP 2055</td><td>Canada</td><td>Triticum turgidum subsp. durum</td><td>CFBP</td><td>379</td><td>4,653,288</td><td>46,036</td><td>CP074361-CP074362</td></tr></tbody></table>

with default settings (Li et al., 2015). The T3SS gene cluster was analyzed further by comparing the sequences of the PGAP-annotated genes using Clinker 0.0.21 with default parameters (Gilchrist and Chooi, 2021). In strains CFBP 8304 and LMG 726, the *xopF* gene found in the cluster was incorrectly annotated with PGAP, but was correctly annotated with Prokka, and was thus manually corrected in the PGAP annotation for the cluster comparison. The same was true for the *xopM* gene in strains UPB458, CFBP 8304, and LMG 728.

The presence of type III effectors was determined by BLASTp using the amino acid sequences of the effectors retrieved from http://xanthomonas.org/ as query (**Supplementary Dataset 2**) and the amino acid sequences of the genes annotated with Prokka as subject, with a 0.01 *e*-value threshold. Only hits that had > 30% identity over > 70% query sequence length were retained. To validate the presence of each effector, the amino acid sequences of the selected genes were then extracted and used as a query for a BLASTp against the type III effectors' sequences as subject. Hits that had > 30% identity over > 70% query sequence length were then considered to be putative type III effectors, as discussed in this article. The presence of TALEs, their RVD sequence, and their classification were determined using AnnoTALE 1.5 (Grau et al., 2016).

# RESULTS

## Genome Assembly

For all the sequenced strains, sequence coverages between 81- and 404-fold were obtained. This allowed to assemble complete genome sequences, consisting of one single circular chromosome for all strains (**Table 1**). Genome sizes ranged from 4,386,175 bp for ICMP 1317 to 4,902,099 for LMG 843. In strain CFBP 2055, an additional circular contig of 46,036 bp was assembled. A comparison by BLAST to the NCBI non-redundant nucleotide sequences database showed homology to plasmid sequences. This

contig was thus considered to represent a plasmid of pv. *undulosa* strain CFBP 2055.

## Phylogeny and Comparative Analysis

Phylogeny based on ANI revealed that though all strains shared > 95% ANI, three distinct groups could be observed, with ANI values above 97% (**Figure 1**). The first group consisted of strains CFBP 2055, CFBP 2539, UPB458, DSM 18794, and CFBP 8304; the second consisted of strains ICMP 16317 and CFBP 2541, and the third consisted of strains LMG 730, LMG 843, LMG 728, LMG 727, and LMG 726. These three groups will thereafter be referred to as clades Xt-I, Xt-II, and Xt-III, respectively. These three clades could also be observed when including all the publicly available complete genome sequences of *X. translucens* strains, with ANI values above 97% (**Supplementary Figure 1**). Interestingly, the two strains regarded as pv. *pistaciae* were genetically distinct and were found in two separate clades.

A comparison of the genomic structure within each clade showed 34–48 locally collinear blocks (LCB) in clade Xt-I (**Figure 2A**), 11 LCB in clade Xt-II (**Figure 2B**), and 175–312 LCB in clade Xt-III (**Figure 2C**). These LCB correspond to genomic regions that are conserved between the compared strains, showing no rearrangement. Thus, the low number of LCB in clade Xt-II showed that there are very few genomic rearrangements between the two strains of the clade, while there are more rearrangements in clade Xt-I and the most in clade Xt-III.

The number of genes found in each strain with the Prokka annotation ranged from 3,735 genes in ICMP 16317 to 4,262 genes in LMG 726. The pangenome of the 12 *X. translucens* strains used in this study consisted of 9,772 genes, while the core genome consisted of 2,181 genes (**Figure 3**). A total of 190 genes were exclusive to clade Xt-I, 588 genes to clade Xt-II, and 168 genes to clade Xt-III. Moreover, a total of 3,681 genes were

<footnote><sup>a</sup>DSMZ, German Collection of Microorganisms and Cell Cultures; CFBP, Collection of Plant Pathogenic Bacteria; ICMP, International Collection of Microorganisms from Plants; BCCM, Belgian Coordinated Collections of Microorganisms; UPB, Collection of plant pathogenic bacteria, Earth and Life Institute, UCLouvain.</footnote>

<footnote><sup>b</sup>UPB458 is the name of the strain in the UPB collection from where it was obtained and corresponds to LMG 737 in the BCCM collection.</footnote>

<footnote><sup>c</sup>Jaenicke et al. (2016).</footnote>
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Average nucleotide identity (ANI) of *X. translucens* pathotype strains and ANI-based phylogeny constructed with Ward's hierarchical clustering method. Distance in the dendrogram represents ANI dissimilarity between nodes. Orange: clade Xt-I, blue: clade Xt-II, pink: clade Xt-III. ANI is depicted as a gradient from blue (< 96%) to white (98%) to red (100%).</figcaption></figure>

exclusive to one strain only, ranging from 48 genes in strain CFBP 2539 to 966 genes in LMG 726.

In all strains, an *xps* T2SS and a T3SS were identified, while no *xcs* T2SS was identified in any strain (Figure 4A). A T6SS-i3*** was identified in all strains of clades Xt-I and Xt-II, while a T6SS-i4 was identified in strains CFBP 2055, CFBP 2539, UPB458, ICMP 16317, CFBP 2541, and LMG 843. A T4SS was identified in strains DSM 18974, CFBP 8304, LMG 730, LMG 728, and LMG 727 which all lacked a T6SS-i4. The strain LMG 726 was the only one that did not possess a T4SS or a T6SS.

A direct comparison of the *hrp* cluster showed that all strains share the same genetic organization of the cluster (Figure 5). Moreover, most of the main components of the T3SS, from *hrcC* to *hrpD*, as well as *hpaH* and *xopF* were very conserved, with > 80% identity across the species. However, the *hrpE* structural component was the most variable, with as little as 60% identity within clades Xt-I and Xt-III (Figure 6). In strains LMG 728 and LMG 726, an additional gene was found between *hrpX* and *hrcT*, in opposite directions in each strain, but showing 100% identity. This gene showed a high identity with insertion sequence (IS) 5 family transposases found in other Xanthomonas species by BLASTx against the NCBI non-redundant protein sequence database. Moreover, in strains LMG 727 and LMG 726, one additional gene was predicted between *hpaB* and *hrpG*, with 61% identity between the two strains. However, no similarity to a known gene was found by BLASTx. Additionally, *hpaC* was found to be truncated in DSM 18974 due to an early stop codon, but still showed between 71 and 86% identity with the other strains (Figure 6).

The T3SS-associated components upstream and downstream of the core *hrp* cluster (*hpaT*, *hpa1*, *xopM*, and *hgiB*) showed more variability. The putative translocon *hpaT* was very conserved

between the two clade Xt-II strains and showed between 63 and 89% identity in clade Xt-I (Figure 6). It was much more variable in clade Xt-III, where it showed between 39 and 67% identity. The second putative translocon component, *hpa1*, was very conserved among the two clade Xt-II strains (98% identity) and among the UPB458, DSM 18974, and CFBP 8304 strains in clade Xt-I. It was, however, more dissimilar between these three strains and the other two of the clade, with 45-47% identity. Within the Xt-III clade, there was a high variability, with 50-79% identity between strains. Moreover, one additional gene was predicted between *hpaT* and *hpa1* in UPB458, but no similarity to any known gene was found by BLASTx.

In clade Xt-I, the XopM effector was only present in UPB458, DSM 18974, and CFBP 8304 and was 99% identical between the three strains. The effector was also found in LMG 728, where it was 92-93% identical to the ones found in clade Xt-I. In the other four strains of clade Xt-III, the gene was disrupted by early stop codons. The *hgiB* gene was also very variable, with as low as 43% identity within clade Xt-I, and 37% identity within clade Xt-III, but with 99% identity between the two clade Xt-II strains. In LMG 728, it was found to be truncated. Furthermore, an additional gene was found between *xopM* and *hgiB* in LMG 726, which showed 100% identity to an IS4 family transposase by BLASTx.

Together with XopF and XopM, a total of 29 putative non-TALE type III effector classes were found to be present in the 12 *X. translucens* strains used in this study, ranging from 21 effectors in LMG 726 to 31 in CFBP 2541 and DSM 18974 (Figure 4B). Among the predicted type III effectors, AvrBs2, as well as effectors of classes XopC, XopF, XopK, XopL, XopN, XopP, XopQ, XopR, XopV, XopX, XopZ, and XopAM were found in all strains. The effectors of classes XopB, XopG, XopE, XopAA, and XopAF were conserved in all strains of clades Xt-I and Xt-II
<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Pairwise comparisons of the genomic structure within the three *X. translucens* clades with progressive Mauve (Darling, 2004). (A) Clade Xt-I, (B) clade Xt-II, (C) clade Xt-III. Colors represent conserved genomic regions (locally collinear blocks, LCBs), i.e., regions with no rearrangement across all the compared genome sequences. Lines between strains link LCBs that are orthologous between two genome sequences. LCBs found on the bottom part represent regions that are in reverse orientation compared to the reference. In each comparison, the sequence on top is used as reference.</figcaption></figure>

but present in only some strains of clade Xt-III. XopJ was found only in clade Xt-I, while XopAD was found only in clades Xt-II and Xt-III. Interestingly, the *xopM* gene located downstream of the *hrp* cluster was not found in this analysis, the XopM effector identified in LMG 728 being a different one. However, this XopM was found to be present in some strains by visual inspection of the gene clusters (Figure 5).

A total of 21 TALE classes were identified in the 12 strains used in this study, with up to 8 in DSM 18974 (Table 2). Clade Xt-I was the clade with the most TALEs. However, the TALE repertoires of the clade were very diverse, with only TalDA being present in all strains of the clade. Interestingly, TalDA was also found in ICMP 16317. Clade Xt-II and Xt-III strains had a much smaller set of TALEs, ranging from zero to three, with most of them being

exclusive to one strain, except for TalIT, which was found in both LMG 843 and LMG 727. Additionally, a potential pseudo-TALE was identified in LMG 727, with only two repeats. However, the second repeat is only 19 amino acids long and many stop codons were found in the N-terminus. LMG 726 was the only strain in which no TALEs were identified.

DISCUSSION

In this study, we generated high-quality complete genome sequences of all pathotype strains of *X. translucens* and make them available as a community resource for in-depth comparative genome analyses within one of the most important pathogenic
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Gene content of the 12 *X. translucens* pathotype strains used in this study. Vertical bars show the number of genes exclusive to the strains marked in the matrix below. The first group represents the core genome of the species. The second, third, and fourth groups represent genes that are exclusive to the three clades, as indicated by their respective color. Other groups represent genes that are exclusive to one strain. Horizontal bars show the total number of genes in each strain. Orange: clade Xt-I, blue: clade Xt-II, pink: clade Xt-III.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | (A) Presence of the xps and xcs type II secretion systems (T2SS), type III secretion system (T3SS), type IV secretion system (T4SS), subtype 3**** type VI secretion system (T6SS-i3**), and subtype 4 type VI secretion system (T6SS-i4) gene clusters in *X. translucens* pathotype strains. (B) Presence of putative type III effectors. Green circles indicate presence, and white circles indicate absence. Numbers denote the number of putative effectors of each class. Orange: clade Xt-I, blue: clade Xt-II, pink: clade Xt-III.</figcaption></figure>

bacterial species. These are the first complete genome sequences for strains of the pathovars *arrhenatheri*, *graminis*, *hordei*, *phlei*, *phleipratensis*, *pistaciae*, *poae*, and *secalis*. These resources complement the complete genome of the pv. *translucens* pathotype strain DSM 18974 (Jaenicke et al., 2016), as well as the already available genome sequences of pv. *cerealis*, *translucens*, and *undulosa* strains.

Phylogeny based on ANI revealed that three genetically distinct groups can be identified, with members of each clade being less than 96% identical to members of the two other clades. This is in contrast with the usual distinction of only two groups, the “translucens” group and the “graminis” group but is in line with the previous studies suggesting that pv. *cerealis* could be genetically distinct from the other pathovars (Peng et al., 2016;
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Pairwise comparisons of the genetic content of the type III secretion system cluster in the *X. translucens* pathotype strains obtained with Clinker (Gilchrist and Chooi, 2021). The phylogenetic tree is based on ANI as shown in **Figure 1**, orange: Clade Xt-I, blue: clade Xt-II, pink: clade Xt-III. Colors of the arrows represent groups of similar genes. Clusters for strains ICMP 16317, CFBP 2541, and LMG 727 are found in opposite direction and were reversed in the figure.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**FIGURE 6** | Pairwise comparisons of the nucleotide sequences of the most variable components of the type III secretion system across all *X. translucens* pathotype strains. For *xopM*, only UPB458, DSM 18974, CFBP 8304, and LMG 728 were considered, as they were the only strains where a complete gene was identified.</figcaption></figure>
<table><caption>**TABLE 2** | Classes and repeat variable diresidue (RVD) sequences of transcription activator-like effectors (TALEs) found in the *X. translucens* pathotype strains.</caption><thead><tr><th>TALE class<sup>a</sup></th><th>Strain</th><th>RVD sequence</th><th>Chromosome position (location:strand)</th></tr></thead><tbody><tr><td rowspan="5">TaICT</td><td>CFBP 2055</td><td>HN-HD-HD-HD-NI-NI-NI-HN-HD-HD-NN-NN-NI-NN-HD</td><td>1,324,426-1,327,651:1</td></tr><tr><td>CFBP 2539</td><td>NN-HD-HD-HD-NI-NI-NI-HN-HD-HD-NN-NN-NI-NN-HD</td><td>1,325,514-1,328,739:1</td></tr><tr><td>DSM 18974</td><td>NN-HD-HD-HD-NI-NI-NI-HD-HD-NN-NN-NI-NN-HD</td><td>671,336-674,561:1</td></tr><tr><td>UPB458</td><td>NN-HD-HD-NI-NI-NI-HD-HD-NN-NN-NI-NN-HD</td><td>3,936,481-3,939,712:-1</td></tr><tr><td>DSM 18974</td><td>NG-HD-HD-HN-NG-NI-HG-HG-HD-ND-NN-NN-NI-NH-QD</td><td>3,565,727-3,568,955:-1</td></tr><tr><td rowspan="2">TaICV</td><td>DSM 18974</td><td>NG-NN-HD-HD-NN-NI-HG-HD-ND-HG-NI-NN-HD</td><td>3,562,396-3,565,408:-1</td></tr><tr><td>UPB458</td><td>NG-NN-HD-HD-NN-NI-HG-HD-ND-HG-NI-NN-HD</td><td>1,087,442-1,090,460:1</td></tr><tr><td>TaICW</td><td>DSM 18974</td><td>NN-NI-HN-HD-NI-NH-NG-HN-HD-HD-NI-QD</td><td>660,311-663,329:1</td></tr><tr><td rowspan="3">TaICX</td><td>CFBP 8304</td><td>NN-HD-NG-NN-HN-KG-NI-HD-NI-NK-HD-HD-HD-NI-HN-NH-HD-QD</td><td>674,110-677,491:1</td></tr><tr><td>DSM 18974</td><td>NN-HD-NG-NI-HN-KG-NI-HD-NI-NH-NG-NN-HD-HD-NI-NN-NI-HD-QD</td><td>663,647-667,295:1</td></tr><tr><td>UPB458</td><td>NN-QD-NG-NN-HN-KG-NI-HD-NI-NH-NG-HN-HD-NI-NN-HD</td><td>3,939,840-3,943,281:-1</td></tr><tr><td>TaICY</td><td>DSM 18974</td><td>NI-NG-HN-NN-HD-NG-ND-NK-QD-NH-QD</td><td>668,858-671,207:1</td></tr><tr><td rowspan="3">TaICZ</td><td>CFBP 2055</td><td>NH-NN-HD-NN-HD-NH-HD-YK-NG-NH-Y*-HD-NN-NI-NG-QD</td><td>1,916,521-1,919,857:-1</td></tr><tr><td>CFBP 2539</td><td>NH-NN-HD-NN-HD-NH-HD-YK-NG-NH-Y*-HD-NN-NI-NG-QD</td><td>1,917,023-1,920,359:-1</td></tr><tr><td>DSM 18974</td><td>NH-NN-HD-NN-HD-NH-HD-YK-NG-NH-Y*-HD-NN-NI-NG-QD</td><td>1,995,271-1,998,607:-1</td></tr><tr><td rowspan="6">TaIDA</td><td>CFBP 2055</td><td>HD-YD-NI-NG-NG-NN-YK-NG-HD-NG-NG-ND-NG-QD-NH-HD</td><td>593,184-596,364:1</td></tr><tr><td>CFBP 2539</td><td>HD-YD-NI-NG-NG-NN-YK-NG-HD-NG-NG-ND-NG-QD-NH-HD</td><td>594,296-597,476:1</td></tr><tr><td>CFBP 8304</td><td>HD-YD-NI-NG-NG-NN-YK-NG-HD-NG-NG-ND-NK-QD-NH-QD</td><td>653,260-656,587:1</td></tr><tr><td>DSM 18974</td><td>NN-HD-NG-NG-NG-NN-YK-NG-HD-NG-NG-ND-NG-HD-NH-HD</td><td>627,243-630,573:1</td></tr><tr><td>ICMP 16317</td><td>HD-YD-NI-NG-NG-NN-YK-NG-HD-NG-NG-ND-NG-QD-NH-HD</td><td>586,137-589,323:1</td></tr><tr><td>UPB458</td><td>NN-HD-NG-NG-NG-NN-YK-NG-HD-NG-NG-ND-NG-HD-NH-HD</td><td>4,001,560-4,004,884:-1</td></tr><tr><td rowspan="3">TaIDC</td><td>CFBP 2055</td><td>NN-NG-HD-HD-HD-KG-NN-Y*-NG-HD-HD-QD-HN</td><td>1,321,073-1,324,088:1</td></tr><tr><td>CFBP 2055</td><td>NN-HD-NG-NN-HN-KG-NI-HD-NI-NN-HD-HN-HD-HD-NI-HN-HD-QD</td><td>612,552-616,092:1</td></tr><tr><td>CFBP 2539</td><td>NN-HD-NG-NN-HN-KG-NI-HD-NI-NN-HD-HN-HD-HD-NI-HN-HD-QD</td><td>613,664-617,204:1</td></tr><tr><td rowspan="2">TaIDE</td><td>CFBP 2539</td><td>NN-HD-NG-NN-HN-KG-NI-HD-NI-NN-HD-HD-NN-NI-HN-HD</td><td>1,321,735-1,325,176:1</td></tr><tr><td>CFBP 2539</td><td>NN-HD-NG-NN-HN-HN-NI-NI-NH-NN-HD-NH-NH-HD</td><td>1,717,212-1,720,551:1</td></tr><tr><td rowspan="2">TaIDF</td><td>CFBP 2539</td><td>NN-HD-NG-NN-HN-HN-NI-NI-NH-NN-HD-NH-NH-HD</td><td>1,716,713-1,720,052:1</td></tr><tr><td>CFBP 2539</td><td>NN-HD-NG-NN-HN-HN-NI-NI-NH-NN-HD-NH-NH-HD</td><td>1,711,799-1,715,336:1</td></tr><tr><td rowspan="2">TaIDO</td><td>CFBP 2539</td><td>HD-HN-HN-HD-NH-NH-HG-HD-ND-NN-Y*-NG-HD-NI-NH-NG-HD-HN</td><td>1,712,298-1,715,835:1</td></tr><tr><td>CFBP 2541</td><td>HD-HN-HN-HD-NH-NH-HG-HD-KG-NN-Y*-NG-HD-NI-NH-NG-HD-HN</td><td>2,549,947-2,553,289:1</td></tr><tr><td>TaDP</td><td>CFBP 2541</td><td>NS-KI-NI-HD-NK-GI-HD-NK-HD-NN-HD-NK</td><td>590,471-593,369:1</td></tr><tr><td>TaJU</td><td>LMG 843</td><td>NS-NG-HD-HD-NN</td><td>1,670,023-1,673,728:-1</td></tr><tr><td rowspan="2">TaJS</td><td>LMG 727</td><td>NN-NG-NN-NG-HD-NK-NI-DD-NK-HD-NG-NN-NI-NG-NN-HD-HD-QD</td><td>3,139,389-3,143,094:1</td></tr><tr><td>LMG 843</td><td>NN-NG-NN-NG-HD-NK-NI-DD-NK-DD-NG-NN-NI-NG-NN-HD-HD-QD</td><td>2,289,805-2,292,265:1</td></tr><tr><td>TaIY</td><td>UPB 458</td><td>NI-NG-HN-NK-HD-NH-HN-HD-HD-HD-HD-QD</td><td>1,816,261-1,819,174:-1</td></tr><tr><td>TaJB</td><td>CFBP 8304</td><td>NH-NN-HE-NK-HD-NK-HD-YK-NG-NH-Y*-HE-NI-NI-NG-QD</td><td>1,986,200-1,989,539:-1</td></tr><tr><td>TaJT</td><td>LMG 730</td><td>NN-NN-HK-HD-HN-HD-NN</td><td>2,152,104-2,154,627:1</td></tr><tr><td>TaJV</td><td>LMG 728</td><td>NN-NI-HD-HD-HD-HD-KT-NG-NN-NN-KT-NI-HN-HD-NG-NG-NN-HD-NK</td><td>2,784,954-2,788,575:-1</td></tr><tr><td>TaJW</td><td>LMG 728</td><td>HK-HD-HN-NI-NG-HD-HN-NI-NG-HD-NG-NN-HN-HD-NG-NN-NI-HD-NG-NN-HD-HD-QD</td><td>2,798,245-2,802,361:-1</td></tr><tr><td>TaJX</td><td>LMG 728</td><td>NK-NG-NI-HD-NG-NN-NG-HD-NK-N*-NK-HD-NN-HD-NG-NN-NI-NG-HD-HD-NN-HD-HD-NN-HD-HD-QD</td><td>1,903,396-1,907,938:-1</td></tr><tr><td>Pseudo TALE</td><td>LMG 727</td><td>NI-NN</td><td>2,236,495-2,238,688:1</td></tr></tbody></table>

Langlois et al., 2017; Shah et al., 2019). Based on this phylogeny, we propose to classify these groups as clade Xt-I, containing pathovars *hordei*, *translucens*, *undulosa*, and *secalis*; clade Xt-II, containing *pv. cerealis*; and clade Xt-III, containing pathovars *arrhenatheri*, *graminis*, *phlei*, *phleipratensis*, and *poae*.

Moreover, the two strains of *pv. pistaciae* were grouped in two different clades, with the group A strain found in clade Xt-I and the group B strain found in clade Xt-II, the two strains being only

95.4% identical. This confirmed the previous phylogeny based on *gyrB* sequences where the group A strain was more closely related to pathovars *translucens*, *secalis*, and *undulosa*, while the group B strain was closest to *pv. cerealis* (Giblot-Ducray et al., 2009). There is clear evidence that these strains are not directly related, although they share the same host.

These results raise some limitations of the classical pathovar classification used in *Xanthomonas*. Indeed, the pathotype strains

<footnote><sup>a</sup>TALE classes based on AnnoTALE (Grau et al., 2016).</footnote>

<footnote>*Thirteenth residue is missing.</footnote>
of pathovars *secalis* and *undulosa* are 99.8% identical and share similar genomic organization and virulence features. As they also have a similar host range, this could lead to them being considered as the same taxonomic entity. On the other hand, the two groups of pv. *pistaciae* are very different genetically and considering them under the same pathovar could hinder the better understanding of their respective biology. Despite these limitations, the pathovar classification was surprisingly robust for the other pathovars. Nonetheless, this work suggests that this classification, currently based solely on pathogenicity tests, should be rethought to better reflect the genetic relationships between pathovars and their evolutionary history.

Within the three clades identified by the ANI-based phylogeny, we found 190, 588, and 168 genes that were specific to clade Xt-I, Xt-II, and Xt-III, respectively, as well as 48–966 genes that were strain specific. However, these numbers are probably biased by the small number of strains included in the comparison and including more strains would result in a smaller number of clade- and strain-specific genes. Nonetheless, these constitute a valuable list of genes that could shape the host range of each clade and/or strain.

All strains had an *xps* T2SS and no *xcs* T2SS, as well as a T3SS with a similar genetic organization, although a few genes of the *hrp* cluster were variable in sequence between strains. Interestingly, strains that possessed a T4SS did not possess a T6SS-i4 and vice versa. As these two secretion systems are both involved in antimicrobial activity, it is possible that their function is redundant in *X. translucens* and only one of them is required to play this role. However, the presence or absence of these two secretion systems does not necessarily reflect the phylogenetic relationships between the strains investigated in this study. Additionally, a T6SS-i3*** was found in all strains of clades Xt-I and Xt-II. However, it was previously hypothesized that this subgroup of T6SS could be non-functional due to the lack of a PAAR and could be complemented by the presence of another subtype of T6SS (Bayer-Santos et al., 2019). Interestingly, no T4SS or T6SS was identified in strain LMG 726, which confirms the previous research that showed that many pv. *graminis* strains lack a T6SS, with only strains Xtg2, Xtg9, Xtg10, and NCPPB3709 found to harbor one (Hersemann et al., 2017).

We have identified a total of 13 type III non-TALE effector classes that could constitute a core set of effectors in *X. translucens*. Moreover, XopAD was not found in the strains of clade Xt-I in our analysis but was previously identified in strain DSM 18974 in different analyses, indicating that it could also be part of the *X. translucens* core set of type III effectors (Peng et al., 2016; Koebnik et al., 2021; Shah et al., 2021). Five additional effector classes constitute a core set of effectors in clades Xt-I and Xt-II, as well as XopJ, which was specific to clade Xt-I. Some effector classes were specific to one strain, such as AvrBs1 in pv. *cerealis* strain CFBP 2541, or XopAH in pv. *hordei* strain UPB458. Furthermore, XopAJ and XopAL class effectors were found in both pv. *translucens* and *hordei*. These effectors could have a role in the host specificity of each clade and/or pathovar. However, no type III effector has yet been functionally characterized in *X. translucens*, and the effectors identified in this study will need to be

functionally validated to confirm their role in the pathogenicity of each pathovar.

High-quality genome assemblies allowed for the identification of the first TALE to be reported in grass-infecting *X. translucens* strains, as well as in pathovars *secalis*, *hordei*, and *pistaciae*. No TALE was found in pv. *graminis* in our analysis, confirming previous research where no TALE was identified in draft genomes of the pathovar (Wichmann et al., 2013). The TALEs identified in pv. *undulosa* correspond to those previously identified in other strains of this pathovar, except for Xt4699-Tal3, for which no similar TALE was found, as was the case in strain ICMP 11055 (Peng et al., 2016; Falahi Charkhabi et al., 2017). The TALE of class TalDC that we identified has the same RVD sequence as Xt4699-Tal8 and a similar sequence to ICMP 11055-Tal4b. Additionally, the TALE of class TalDD identified in CFBP 2055 and CFBP 2539 has a similar RVD sequence to ICMP 11055-Tal2. Furthermore, the TALEs identified in pv. *cerealis* strain CFBP 2541 correspond to the TALE previously identified in strain Nxtc01, with the TALE of class TalDP having a similar sequence to Nxtc01-Tal1 (Shah et al., 2019). As these four TALEs have previously been shown to have a role in virulence, the corresponding effectors we identified in this study could play a similar role in their respective strain (Falahi Charkhabi et al., 2017; Peng et al., 2019; Shah et al., 2019). However, these are the only TALEs that have been functionally characterized in *X. translucens* to date. Nonetheless, as most of the TALEs identified in this study are clade or strain specific, they could be essential components of host adaptation. Indeed, TALEs activate the transcription of plant genes by binding to their promoter region. As these sequences vary between plant species, TALEs must be adapted to target the specific sequence of that promoter in their host and thus reflect processes of co-adaptation between the bacteria and its host (Jacques et al., 2016). However, the lack of TALEs in pv. *graminis* indicates that other mechanisms play a role in *X. translucens* host speciation.

The pv. *pistaciae* group A pathotype strain shared a very similar genomic organization, T3SS cluster organization, and gene identity, as well as type III effector and TALE repertoires with the pathotype strains of pv. *translucens* and *hordei*. The same was true for the group B pv. *pistaciae* strain and the pv. *cerealis* pathotype strain, which were even more closely related. This could indicate that a very small set of genes is responsible for their ability to infect pistachio, as a previous study showed that pv. *translucens* strain DAR 35705 was unable to cause symptoms in pistachio and other Anacardiaceae, while the pv. *pistaciae* strains of both groups were pathogenic in Anacardiaceae and Poaceae (Marefat et al., 2006). However, no type III effector or TALE exclusive to either strain of pv. *pistaciae* has been identified in this study.

Although it had the biggest set of total genes, the pv. *graminis* strain LMG 726 was found to lack many virulence features such as the T4SS, the T6SS, as well as TALEs and had the smallest set of type III effectors. This is surprising, as this pathovar is known to be the most widespread and virulent grass-infecting *X. translucens* pathovar, with the largest host range in clade Xt-III. Indeed, while pv. *phlei* and pv.
phleipratensis, pv. arrhenatheri, and pv. poae are restricted to the genera Phleum, Arrhenatherum, and Poa, respectively, pv. graminis can infect many grasses from genera such as Agrostis, Alopecurus, Dactylis, Deschampsia, Festuca, Lolium, Phalaris, Phleum, Poa, and Trisetum (Egli et al., 1975; Egli and Schmidt, 1982). Nonetheless, these results go along previous research that showed that strains of pv. graminis also lack a flagellum and have a distinct type IV pilus compared to other X. translucens pathovars (Hersemann et al., 2017). It was hypothesized that since pv. graminis is usually spread by mowing tools, the lack of a flagellum might not hinder its ability to spread in the plant, as flagellar motility is mostly necessary for the bacteria to reach points of entry in the plant. Similarly, it might encounter less competition with epiphytic microorganisms, and the presence of a T4SS and/or a T6SS might not be crucial for its survival. Additionally, as these features can act as elicitors of plant defense, their absence in pv. graminis could help to evade such defense mechanisms—an intriguing hypothesis that, however, has not yet been tested.

## CONCLUSION

In conclusion, our study substantially increased the number of complete genome sequences available for X. translucens, providing high-quality genomic resources for all the pathovars of the species. These sequences constitute a valuable basis for future studies investigating the phylogenetic relationships between X. translucens pathovars and other Xanthomonas species, as well as their key genetic features. The virulence features of X. translucens that we identified will help to better understand the biology of each clade and/or pathovar, and processes of adaptation to their respective hosts. The inclusion of the high-quality genome sequences of the additional strains of each pathovar in further comparative genomics studies will help in refining this list, laying the foundations for the development of new resistant cereal crop and forage grass cultivars.

## DATA AVAILABILITY STATEMENT

The genome sequences generated in this study can be found on https://www.ncbi.nlm.nih.gov with accession numbers CP074361-CP074365, CP076249-CP076251, CP076254, CP083804, CP086332-CP086333, and LT604072

## AUTHOR CONTRIBUTIONS

FG and RoK: design of the study with the help of RaK, CB, JJ, VR-R, and BS. FG, VR-R, RaK, and SC: genome

## REFERENCES

Alvarez-Martinez, C. E., Sgro, G. G., Araujo, G. G., Paiva, M. R. N., Matsuyama, B. Y., Guzzo, C. R., et al. (2021). Secrete or perish: the role of secretion systems in Xanthomonas biology. *Comput. Struct. Biotechnol. J.* 19, 279–302. doi: 10.1016/j.csbj.2020.12.020

sequencing and assemblies. FG: data analysis with assistance from VR-R. FG, RoK, RaK, JJ, VR-R, CB, SC, and BS: writing and reviewing. All authors contributed to the article and approved the submitted version.

## FUNDING

Funding was provided by the Swiss National Science Foundation (Grant No. IZCOZO_177062). This research was also supported by the United States Department of Agriculture (USDA) National Institute of Food and Agriculture (NIFA) (Award No. 2018-67013-28490) through the Joint National Science Foundation-NIFA Plant Biotic Interactions Program and the USDA NIFA Food and Agriculture Cyberinformatics and Tools grant program (Grant No. 2021-67021-34343) to JJ. Work in RaK's laboratory was supported by a grant from the French National Research Agency (Grant No. CROpTAL ANR-14-CE19-0002).

## ACKNOWLEDGMENTS

Library preparation and sequencing were done by the Functional Genomics Center Zurich (Zurich, Switzerland; https://fgcz.ch/) and the GENTYANE genotyping platform (INRA Clermont-Ferrand, France; https://gentyane.clermont.inrae.fr). Technical assistance by Ildiko Katalin Nagy (Enviroinvest Corp. Pecs, Hungary) is greatly appreciated. RaK and RoK are grateful for an international mobility EXPLORE grant from the Montpellier University of Excellence (MUSE). This article is based upon work from COST Action CA16107 EuroXanth, supported by COST (European Cooperation in Science and Technology).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb. 2021.817815/full#supplementary-material

**Supplementary Figure 1** | Average nucleotide identity (ANI) of all available complete genome sequences of Xanthomonas translucens strains and ANI-based phylogeny constructed with Ward's hierarchical clustering method. Distance in the dendrogram represents ANI dissimilarity between nodes. Orange: clade Xt-I, blue: clade Xt-II, pink: clade Xt-III. ANI is depicted as a gradient from blue (<96 %) to white (98%) to red (100%).

**Supplementary Dataset 1** | Amino acid sequences of all secretion system components used to predict the presence of these secretion systems by tBLASTn.

**Supplementary Dataset 2** | Amino acid sequences of all type III effectors (T3Es) used to predict the presence of these effectors by BLASTp.

Bayer-Santos, E., Ceseti, L., de, M., Farah, C. S., and Alvarez-Martinez, C. E. (2019). Distribution, function and regulation of type 6 secretion systems of Xanthomonadales. *Front. Microbiol.* 10:1635. doi: 10.3389/fmicb.2019. 01635

Bogdanove, A. J., Beer, S. V., Bonas, U., Boucher, C. A., Collmer, A., Coplin, D. L., et al. (1996). Unified nomenclature for broadly conserved *hrp* genes
of phytopathogenic bacteria. *Mol. Microbiol.* 20, 681–683. doi: 10.1046/j.1365-2958.1996.5731077.x

Bonas, U., Schulte, R., Fenselau, S., Mainsavage, G. V., Staskawicz, B. J., and Stall, R. E. (1991). Isolation of a gene cluster from *Xanthomonas campestris* pv. vesicatoria that determines pathogenicity and the hypersensitive response on and the hypersensitive response on pepper and tomato. *Mol. Plant Microbe Interact.* 4, 81–88.

Booher, N. J., Carpenter, S. C. D., Sebra, R. P., Wang, L., Salzberg, S. L., Leach, J. E., et al. (2015). Single molecule real-time sequencing of *Xanthomonas oryzae* genomes reveals a dynamic structure and complex TAL (transcription activator-like) effector gene relationships. *Microb. Genom.* 1:e000032. doi: 10.1099/mgen. 0.000032

Bragard, C., Singer, E., Alizadeh, A., Vauterin, L., Maraite, H., and Swings, J. (1997). *Xanthomonas translucens* from small grains: diversity and phytopathological relevance. *Phytopathology*® 87, 1111–1117. doi: 10.1094/PHYTO.1997.87.11. 1111

Büttner, D., and Bonas, U. (2010). Regulation and secretion of Xanthomonas virulence factors. *FEMS Microbiol. Rev.* 34, 107–133. doi: 10.1111/j.1574-6976. 2009.00192.x

Chin, C.-S., Alexander, D. H., Marks, P., Klammer, A. A., Drake, J., Heiner, C., et al. (2013). Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. *Nat. Methods* 10, 563–569. doi: 10.1038/nmeth.2474

Choi, Y., Kim, N., Mannaa, M., Kim, H., Park, J., Jung, H., et al. (2020). Characterization of type VI secretion system in *Xanthomonas oryzae* pv. oryzae and its role in virulence to rice. *Plant. Pathol. J.* 36, 289–296. doi: 10.5423/PPJ. NT.02.2020.0026

Darling, A. C. E. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. *Genome Res.* 14, 1394–1403. doi: 10.1101/gr. 2289704

Egli, T., Goto, M., and Schmidt, D. (1975). Bacterial wilt, a new forage grass disease. *J. Phytopathol.* 82, 111–121. doi: 10.1111/j.1439-0434.1975. tb02825.x

Egli, T., and Schmidt, D. (1982). Pathogenic variation among the causal agents of bacterial wilt of forage grasses. *J. Phytopathol.* 104, 138–150. doi: 10.1111/ j.1439-0434.1982.tb00520.x

Facelli, E., Taylor, C., Scott, E., Fegan, M., Huys, G., Noble, R. D., et al. (2005). Identification of the causal agent of pistachio dieback in Australia. *Eur. J. Plant Pathol.* 112, 155–165. doi: 10.1007/s10658-005-3120-9

Falahi Charkhabi, N., Booher, N. J., Peng, Z., Wang, L., Rahimian, H., Shams-Bakhsh, M., et al. (2017). Complete genome sequencing and targeted mutagenesis reveal virulence contributions of Tal2 and Tal4b of *Xanthomonas translucens* pv. *undulosa* ICMP11055 in bacterial leaf streak of wheat. *Front. Microbiol.* 8:1488. doi: 10.3389/fmicb.2017.01488

Giblot-Ducray, D., Marefat, A., Gillings, M. R., Parkinson, N. M., Bowman, J. P., Ophel-Keller, K., et al. (2009). Proposal of *Xanthomonas translucens* pv. *pistaciae* pv. nov., pathogenic to pistachio (*Pistacia vera*). *Syst. Appl. Microbiol.* 32, 549–557. doi: 10.1016/j.syapm.2009.08.001

Gilchrist, C. L. M., and Chooi, Y.-H. (2021). Clinker & clustermap.js: automatic generation of gene cluster comparison figures. *Bioinformatics* 37, 2473–2475. doi: 10.1093/bioinformatics/btab007

Grau, J., Reschke, M., Erkes, A., Streubel, J., Morgan, R. D., Wilson, G. G., et al. (2016). AnnoTALE: bioinformatics tools for identification, annotation and nomenclature of TALEs from *Xanthomonas* genomic sequences. *Sci. Rep.* 6:21077. doi: 10.1038/srep21077

Hersemann, L., Wibberg, D., Blom, J., Goesmann, A., Widmer, F., Vorhölter, F.-J., et al. (2017). Comparative genomics of host adaptive traits in *Xanthomonas translucens* pv. *graminis*. *BMC Genomics* 18:35. doi: 10.1186/s12864-016- 3422-7

Jacques, M.-A., Arlat, M., Boulanger, A., Boureau, T., Carrère, S., Cesbron, S., et al. (2016). Using ecology, physiology, and genomics to understand host specificity in *Xanthomonas*. *Annu. Rev. Phytopathol.* 54, 163–187. doi: 10.1146/annurev- phyto-080615-100147

Jaenicke, S., Bunk, B., Wibberg, D., Spröer, C., Hersemann, L., Blom, J., et al. (2016). Complete genome sequence of the barley pathogen *Xanthomonas translucens* pv. *translucens* DSM 18974 T (ATCC 19319 T). *Genome Announc.* 4:e01334–16. doi: 10.1128/genomeA.01334–16

Jha, G., Rajeshwari, R., and Sonti, R. V. (2005). Bacterial type two secretion system secreted proteins: double-edged swords for plant pathogens. *Mol. Plant Microbe Interact*®. 18, 891–898. doi: 10.1094/MPMI-18-0891

Koebnik, R., Rurokiene, D., Bragard, C., Chang, C., Saux, M. F.-L., Kölliker, R., et al. (2021). The complete genome sequence of *Xanthomonas theicola*, the causal agent of canker on tea plants, reveals novel secretion systems in clade-1 xanthomonads. *Phytopathology*® 111, 611–616. doi: 10.1094/PHYTO-07-20- 0273-SC

Kolmogorov, M., Yuan, J., Lin, Y., and Pevzner, P. A. (2019). Assembly of long, error-prone reads using repeat graphs. *Nat. Biotechnol.* 37, 540–546. doi: 10. 1038/s41587-019-0072-8

Langlois, P. A., Snelling, J., Hamilton, J. P., Bragard, C., Koebnik, R., Verdier, V., et al. (2017). Characterization of the *Xanthomonas translucens* complex using draft genomes, comparative genomics, phylogenetic analysis, and diagnostic lamp assays. *Phytopathology*® 107, 519–527. doi: 10.1094/PHYTO-08-16- 0286-R

Li, J., Yao, Y., Xu, H. H., Hao, L., Deng, Z., Rajakumar, K., et al. (2015). SecReT6: a web-based resource for type VI secretion systems found in bacteria: genomics update. *Environ. Microbiol.* 17, 2196–2202. doi: 10.1111/1462-2920.12794

Marefat, A., Scott, E. S., Ophel-Keller, K., and Sedgley, M. (2006). Genetic, phenotypic and pathogenic diversity among xanthomonads isolated from pistachio (*Pistacia vera*) in Australia. *Plant Pathol.* 55, 639–649. doi: 10.1111/ j.1365-3059.2006.01437.x

Page, A. J., Cummins, C. A., Hunt, M., Wong, V. K., Reuter, S., Holden, M. T. G., et al. (2015). Roary: rapid large-scale prokaryote pan genome analysis. *Bioinformatics* 31, 3691–3693. doi: 10.1093/bioinformatics/btv421

Peng, Z., Hu, Y., Xie, J., Potnis, N., Akhunova, A., Jones, J., et al. (2016). Long read and single molecule DNA sequencing simplifies genome assembly and TAL effector gene analysis of *Xanthomonas translucens*. *BMC Genomics* 17:21. doi: 10.1186/s12864-015-2348-9

Peng, Z., Hu, Y., Zhang, J., Huguet-Tapia, J. C., Block, A. K., Park, S., et al. (2019). *Xanthomonas translucens* commandeers the host rate-limiting step in ABA biosynthesis for disease susceptibility. *Proc. Natl. Acad. Sci. U.S.A* 116, 20938–20946. doi: 10.1073/pnas.1911660116

Pesce, C., Bolot, S., Cunnac, S., Portier, P., Fischer-Le Saux, M., Jacques, M.-A., et al. (2015). High-quality draft genome sequence of the *Xanthomonas translucens* pv. *cerealis* pathotype strain cfbp 2541. *Genome Announc.* 3:e01574-14. doi: 10.1128/genomeA.01574-14

Pesce, C., Jacobs, J. M., Berthelot, E., Perret, M., Vancheva, T., Bragard, C., et al. (2017). Comparative genomics identifies a novel conserved protein, hpaT, in proteobacterial type III secretion systems that do not possess the putative translocon protein hrpF. *Front. Microbiol.* 8:1177. doi: 10.3389/fmicb.2017. 01177

Pritchard, L., Glover, R. H., Humphris, S., Elphinstone, J. G., and Toth, I. K. (2016). Genomics and taxonomy in diagnostics for food security: soft-rotting enterobacterial plant pathogens. *Anal. Methods* 8, 12–24. doi: 10.1039/ C5AY02550H

Rademaker, J. L. W., Norman, D. J., Forster, R. L., Louws, F. J., Schultz, M. H., and de Bruijn, F. J. (2006). Classification and identification of *Xanthomonas translucens* isolates, including those pathogenic to ornamental asparagus. *Phytopathology*® 96, 876–884. doi: 10.1094/PHYTO-96-0876

Roman-Reyna, V., Luna, E. K., Pesce, C., Vancheva, T., Chang, C., Ziegle, J., et al. (2020). Genome resource of barley bacterial blight and leaf streak pathogen *Xanthomonas translucens* pv. *translucens* strain UPB886. *Plant Dis.* 104, 13–15. doi: 10.1094/PDIS-05-19-1103-A

Sapkota, S., Mergoum, M., and Liu, Z. (2020). The translucens group of *Xanthomonas translucens*: complicated and important pathogens causing bacterial leaf streak on cereals. *Mol. Plant Pathol.* 21, 291–302. doi: 10.1111/ mpp.12909

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. *Bioinformatics* 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

Sgro, G. G., Oka, G. U., Souza, D. P., Cenens, W., Baver-Santos, E., Matsuyama, B. Y., et al. (2019). Bacteria-killing type IV secretion systems. *Front. Microbiol.* 10:1078. doi: 10.3389/fmicb.2019.01078

Shah, S. M. A., Haq, F., Ma, W., Xu, X., Wang, S., Xu, Z., et al. (2019). Tal1NXtc01 in *Xanthomonas translucens* pv. *cerealis* contributes to virulence in bacterial leaf streak of wheat. *Front. Microbiol.* 10:2040. doi: 10.3389/fmicb.2019. 02040

Shah, S. M. A., Khojasteh, M., Wang, Q., Taghavi, S. M., Xu, Z., Khodaygan, P., et al. (2021). Genomics-enabled novel insight into the pathovar-specific population structure of the bacterial leaf streak pathogen *Xanthomonas translucens* in small grain cereals. *Front. Microbiol.* 12:674952. doi: 10.3389/fmicb.2021.674952
Souza, D. P., Oka, G. U., Alvarez-Martinez, C. E., Bisson-Filho, A. W., Dunger, G., Hobeika, L., et al. (2015). Bacterial killing via a type IV secretion system. *Nat. Commun.* 6:6453. doi: 10.1038/ncomms7453

Stead, D. E. (1989). Grouping of Xanthomonas campestris pathovars of cereals and grasses by fatty acid profiling 1. *EPPO Bull.* 19, 57–68.

Streubel, J., Baum, H., Grau, J., Stuttman, J., and Boch, J. (2017). Dissection of TALE-dependent gene activation reveals that they induce transcription cooperatively and in both orientations. *PLoS One* 12:e0173580. doi: 10.1371/journal.pone.0173580

Szczesny, R., Jordan, M., Schramm, C., Schulz, S., Cogez, V., Bonas, U., et al. (2010). Functional characterization of the Xcs and Xps type II secretion systems from the plant pathogenic bacterium Xanthomonas campestris pv vesicatoria. *New Phytol.* 187, 983–1002. doi: 10.1111/j.1469-8137.2010.03312.x

Vauterin, L., Hoste, B., Kersters, K., and Swings, J. (1995). Reclassification of *Xanthomonas*. *Int. J. Syst. Bacteriol.* **45**, 472–489. doi: 10.1099/00207713-45-3-472

Wang, L., Rinaldi, F. C., Singh, P., Doyle, E. L., Dubrow, Z. E., Tran, T. T., et al. (2017). TAL effectors drive transcription bidirectionally in plants. *Mol. Plant* 10, 285–296. doi: 10.1016/j.molp.2016.12.002

White, F. F., Potnis, N., Jones, J. B., and Koebnik, R. (2009). The type III effectors of Xanthomonas. *Mol. Plant Pathol.* 10, 749–766. doi: 10.1111/j.1364-3703.2009.00590.x

Wichmann, F., Vorhölter, F.-J., Hersemann, L., Widmer, F., Blom, J., Niehaus, K., et al. (2013). The noncanonical type III secretion system of Xanthomonas translucens pv. graminis is essential for forage grass infection. *Mol. Plant Pathol.* 14, 576–588. doi: 10.1111/mpp.12030

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2022 Goettelmann, Roman-Reyna, Cunnac, Jacobs, Bragard, Studer, Koebnik and Kölliker. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.