Complete sequence and detailed analysis of the first indigenous plasmid from Xanthomonas oryzae pv. oryzicola

Phenotypic and Genomic Characterization of Pseudomonas wuhanensis sp. nov., a Novel Species with Promising Features as a Potential Plant Growth-Promoting and Biocontrol Agent

Plant growth-promoting rhizobacterial strain FP607T was isolated from the rhizosphere of beets in Wuhan, China. Strain FP607T exhibited significant antagonism toward several phytopathogenic bacteria, indicating that FP607T may produce antimicrobial metabolites and has a stronger biocontrol efficacy against plant pathogens. Growth-promoting tests showed that FP607T produced indole-3-acetic acid (IAA), NH3, and ferritin. The genome sequence of strain FP607T was 6,590,972 bp long with 59.0% G + C content. The optimum temperature range was 25-30 °C, and the optimum pH was 7. The cells of strain FP607T were Gram-negative, short, and rod-shaped, with polar flagella. The colonies on the King's B (KB) agar plates were light yellow, smooth, and circular, with regular edges. A phylogenetic analysis of the 16S rRNA sequence and a multilocus sequence analysis (MLSA) showed that strain FP607T was most closely related to the type of strain Pseudomonas farris SWRI79T. Based on a polyphasic taxonomic approach, strain FP607T was identified as a novel species within the genus Pseudomonas, for which the name Pseudomonas wuhanensis sp. nov. was proposed. The type of strain used was FP607T (JCM 35688, CGMCC 27743, and ACCC 62446).

Live tracking of a plant pathogen outbreak reveals rapid and successive, multidecade plasmid reduction

Quickly understanding the genomic changes that lead to pathogen emergence is necessary to launch mitigation efforts and reduce harm. In this study, we tracked in real time a 2022 bacterial plant disease outbreak in U.S. geraniums (Pelargonium × hortorum) caused by Xhp2022, a novel lineage of Xanthomonas hortorum. Genomes from 31 Xhp2022 isolates from seven states showed limited chromosomal variation and all contained a single plasmid (p93). Time tree and single nucleotide polymorphism whole-genome analysis estimated that Xhp2022 emerged within the last decade. The phylogenomic analysis determined that p93 resulted from the cointegration of three plasmids (p31, p45, and p66) found sporadically across isolates from previous outbreaks. Although p93 had a 49 kb nucleotide reduction, it retained putative fitness genes, which became predominant in the 2022 outbreak. Overall, we demonstrated, through rapid whole-genome sequencing and analysis, a recent, traceable event of genome reduction for niche adaptation typically observed over millennia in obligate and fastidious pathogens.IMPORTANCEThe geranium industry, valued at $4 million annually, faces an ongoing Xanthomonas hortorum pv. pelargonii (Xhp) pathogen outbreak. To track and describe the outbreak, we compared the genome structure across historical and globally distributed isolates. Our research revealed Xhp population has not had chromosome rearrangements since 1974 and has three distinct plasmids. In 2012, we found all three plasmids in individual Xhp isolates. However, in 2022, the three plasmids co-integrated into one plasmid named p93. p93 retained putative fitness genes but lost extraneous genomic material. Our findings show that the 2022 strain group of the bacterial plant pathogen Xanthomonas hortorum underwent a plasmid reduction. We also observed several Xanthomonas species from different years, hosts, and continents have similar plasmids to p93, possibly due to shared agricultural settings. We noticed parallels between genome efficiency and reduction that we see across millennia with obligate parasites with increased niche specificity.

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.

Cloning and Molecular Characterization of the Recombinant CVB4E2 Immunogenic Viral Protein (rVP1), as a Potential Subunit Protein for Vaccine and Immunodiagnostic Reagent Candidate

The aim of the present study was, first, to clone the VP1 gene of the human coxsackievirus B4 strain E2 (CVB4E2) in the prokaryotic pUC19 plasmid expression vector then to compare it with the structural capsid proteins of the same strain using bioinformatic tools. PCR colony amplification followed through a restriction digestion analysis and sequencing process which affirmed the success of the cloning process. SDS-PAGE and Western Blotting were used to characterize the purified recombinant viral protein expressed in bacteria cells. The BLASTN tool revealed that the nucleotide sequence of the recombinant VP1 (rVP1) expressed by pUC19 highly matched the target nucleotide sequence of the diabetogenic CVB4E2 strain. Secondary structure and three-dimension structure prediction suggested that rVP1, such as wild-type VP1, is chiefly composed of random coils and a high percentage of exposed amino acids. Linear B-cell epitope prediction showed that several antigenic epitopes are likely present in rVP1 and CVB4E2 VP1 capsid protein. Additionally, phosphorylation site prediction revealed that both proteins may affect the signal transduction of host cells and can be involved in virus virulence. The present work highlights the usefulness of cloning and bioinformatics characterizations for gene investigation. Furthermore, the collected data are helpful for future experimental research related to the development of immunodiagnostic reagents and subunit vaccines based on the expression of immunogenic viral capsid proteins.

The Complete Genome Resource of Xanthomonas oryzae pv. oryzae CIX2779 Includes the First Sequence of a Plasmid for an African Representative of This Rice Pathogen

The bacterial plant pathogen Xanthomonas oryzae pv. oryzae is responsible for the foliar rice bacterial blight disease. Genetically contrasted, continent-specific, sublineages of this species can cause important damages to rice production both in Asia and Africa. We report on the genome of the CIX2779 strain of this pathogen, previously named NAI1 and originating from Niger. Oxford Nanopore long reads assembly and Illumina short reads polishing produced a genome sequence composed of a 4,725,792-bp circular chromosome and a 39,798-bp-long circular plasmid designated pCIX2779_1. The chromosome structure and base-level sequence are highly related to reference strains of African X. oryzae pv. oryzae and encode identical transcription activator-like effectors for virulence. Importantly, our in silico analysis strongly indicates that pCIX2779_1 is a genuine conjugative plasmid, the first indigenous one sequenced from an African strain of the X. oryzae species. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Horizontal transfer of the multidrug resistance plasmid RP4 inhibits ammonia nitrogen removal dominated by ammonia-oxidizing bacteria

Antibiotic resistance genes (ARGs) have become an important public health concern. Particularly, although several ARGs have been identified in wastewater treatment plants (WWTPs), very few studies have characterized their impacts on reactor performance. Therefore, our study sought to investigate the effect of a representative conjugative transfer plasmid (RP4) encoding multidrug resistance genes on ammonia oxidation. To achieve this, we established sequencing batch reactors (SBRs) and a conjugation model with E. coli donor strains carrying the RP4 plasmid and a typical ammonia-oxidating (AOB) bacterial strain (Nitrosomonas europaea ATCC 25978) as a recipient to investigate the effect of conjugative transfer of plasmid RP4 on AOB. Our findings demonstrated that the RP4 plasmid carried by the donor strains could be transferred to AOB in the SBR and to Nitrosomonas europaea ATCC 25978. In SBR treated with donor strains carrying the RP4 plasmid, ammonia removal efficiency continuously decreased to 71%. Once the RP4 plasmid entered N. europaea ATCC 25978 in the conjugation model, ammonia removal was significantly inhibited and nitrite generation was decreased. Furthermore, the expression of several functional genes related to ammonia oxidation in AOB was suppressed following the transfer of the RP4 plasmid, including amoA, amoC, hao, nirK, and norB. In contrast, the cytL gene encoding cytochrome P460 was upregulated. These results demonstrated the ecological risk of ARGs in WWTPs, and therefore measures must be taken to avoid their transfer.

A network approach to decipher the dynamics of Lysobacteraceae plasmid gene sharing

Plasmids provide an efficient vehicle for gene sharing among bacterial populations, playing a key role in bacterial evolution. Network approaches are particularly suitable to represent multipartite relationships and are useful tools to characterize plasmid-mediated gene sharing events. The Lysobacteraceae bacterial family gathers plant commensal, plant pathogenic and opportunistic human pathogens for which plasmid mediated adaptation was reported. We searched for homologues of plasmid gene sequences from this family in all the diversity of available bacterial genome sequences and built a network of plasmid gene sharing from the results. While plasmid genes are openly shared between the bacteria of the Lysobacteraceae family, taxonomy strongly defined the boundaries of these exchanges, that only barely reached other families. Most inferred plasmid gene sharing events involved a few genes only, and evidence of full plasmid transfers were restricted to taxonomically close taxon. We detected multiple plasmid-chromosome gene transfers, among which the otherwise known sharing of a heavy metal resistance transposon. In the network, bacterial lifestyles shaped sub-structures of isolates colonizing specific ecological niches and harboring specific types of resistance genes. Genes associated to pathogenicity or antibiotic and metal resistance were among those that most importantly structured the network, highlighting the imprints of human-mediated selective pressure on pathogenic populations. A massive sequencing effort on environmental Lysobacteraceae is therefore required to refine our understanding on how this reservoir fuels the emergence and the spread of genes amongst this family and its potential impact on plant, animal and human health.

Diversity of transcription activator-like effectors and pathogenicity in strains of Xanthomonas oryzae pv. oryzicola from Yunnan

The present study aimed to evaluate transcriptional activator-like effector (TALE) genes in 86 Xanthomonas oryzae pv. oryzicola strains collected from 8 rice-growing regions in Yunnan, and to examine the relationship between TALE genotypes and virulence in 6 differential rice lines. Besides, the geographical areas, distribution of these genotypes were studied in detail. Genetic diversity was analyzed through the number and size of putative TALE genes based on TALE gene avrXa3 as a probe. We found that X. oryzae pv. oryzicola strains consist of variable number (13-27) of avrXa3-hybridizing fragments (putative TALE genes). Test strains were classified into 8 genotypes (G1-G8) with major genotypes G3 and G7 widely distributed in Yunnan. Pathogenicity of X. oryzae pv. oryzicola was evaluated by inoculating 6 differential rice lines with a single resistance gene into 9 pathotypes clusters (I-IX), the dominant Genotypes G3 and G7 consist of pathotypes I, II, and IV. Furthermore, we also detected the known TALE target genes expression in susceptible rice cultivar (cv. nipponbare) after inoculating 8 genotypes-representative X. oryzae pv. oryzicola strain. Correlation between the numbers of putative TALE genes of X. oryzae pv. oryzicola and relevant target genes in nipponbare confirmed up-regulation. Altogether, this study has given insights into the population structure of X. oryzae pv. oryzicola that may inform strategies to control BLS in rice.

A comparative genomic analysis of Xanthomonas arboricola pv. juglandis strains reveal hallmarks of mobile genetic elements in the adaptation and accelerated evolution of virulence

Xanthomonas arboricola pv. juglandis (Xaj) is the most significant aboveground walnut bacterial pathogen. Disease management uses copper-based pesticides which induce pathogen resistance. We examined the genetic repertoire associated with adaptation and virulence evolution in Xaj. Comparative genomics of 32 Xaj strains reveal the possible acquisition and propagation of virulence factors via insertion sequences (IS). Fine-scale annotation revealed a Tn3 transposon (TnXaj417) encoding copper resistance genes acquired by horizontal gene transfer and associated with adaptation and tolerance to metal-based pesticides commonly used to manage pathogens in orchard ecosystems. Phylogenomic analysis reveals IS involvement in acquisition and diversification of type III effector proteins ranging from two to eight in non-pathogenic strains, 16 to 20 in pathogenic strains, besides six other putative effectors with a reduced identity degree found mostly among pathogenic strains. Yersiniabactin, xopK, xopAI, and antibiotic resistance genes are also located near ISs or inside genomic islands and structures resembling composite transposons.

Genomics-Enabled Novel Insight Into the Pathovar-Specific Population Structure of the Bacterial Leaf Streak Pathogen Xanthomonas translucens in Small Grain Cereals

The Gram-negative bacterium Xanthomonas translucens infects a wide range of gramineous plants with a notable impact on small grain cereals. However, genomics-informed intra-species population structure and virulence repertories of the pathogen have rarely been investigated. In this study, the complete genome sequences of seven X. translucens strains representing an entire set of genetic diversity of two pathovars X. translucens pv. undulosa and X. translucens pv. translucens is provided and compared with those of seven publicly available complete genomes of the pathogen. Organization of the 25 type III secretion system genes in all the 14 X. translucens strains was exactly the same, while TAL effector genes localized singly or in clusters across four loci in X. translucens pv. translucens and five to six loci in X. translucens pv. undulosa. Beside two previously unreported endogenous plasmids in X. translucens pv. undulosa, and variations in repeat variable diresidue (RVD) of the 14 strains, tal1a of X. translucens pv. translucens strain XtKm8 encode the new RVDs HE and YI which have not previously been reported in xanthomonads. Further, a number of truncated tal genes were predicted among the 14 genomes lacking conserved BamHI site at N-terminus and SphI site at C-terminus. Our data have doubled the number of complete genomes of X. translucens clarifying the population structure and genomics of the pathogen to pave the way in the small grain cereals industry for disease resistance breeding in the 21st century's agriculture.

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

Horizontal gene transfer is of major evolutionary importance as it allows for the redistribution of phenotypically important genes among lineages. Such genes with essential functions include those involved in resistance to antimicrobial compounds and virulence factors in pathogenic bacteria. Understanding gene turnover at microevolutionary scales is critical to assess the pace of this evolutionary process. Here, we characterized and quantified gene turnover for the epidemic lineage of a bacterial plant pathogen of major agricultural importance worldwide. Relying on a dense geographic sampling spanning 39 years of evolution, we estimated both the dynamics of single nucleotide polymorphism accumulation and gene content turnover. We identified extensive gene content variation among lineages even at the smallest phylogenetic and geographic scales. Gene turnover rate exceeded nucleotide substitution rate by three orders of magnitude. Accessory genes were found preferentially located on plasmids, but we identified a highly plastic chromosomal region hosting ecologically important genes such as transcription activator-like effectors. Whereas most changes in the gene content are probably transient, the rapid spread of a mobile element conferring resistance to copper compounds widely used for the management of plant bacterial pathogens illustrates how some accessory genes can become ubiquitous within a population over short timeframes.

Secrete or perish: The role of secretion systems in Xanthomonas biology

Bacteria of the Xanthomonas genus are mainly phytopathogens of a large variety of crops of economic importance worldwide. Xanthomonas spp. rely on an arsenal of protein effectors, toxins and adhesins to adapt to the environment, compete with other microorganisms and colonize plant hosts, often causing disease. These protein effectors are mainly delivered to their targets by the action of bacterial secretion systems, dedicated multiprotein complexes that translocate proteins to the extracellular environment or directly into eukaryotic and prokaryotic cells. Type I to type VI secretion systems have been identified in Xanthomonas genomes. Recent studies have unravelled the diverse roles played by the distinct types of secretion systems in adaptation and virulence in xanthomonads, unveiling new aspects of their biology. In addition, genome sequence information from a wide range of Xanthomonas species and pathovars have become available recently, uncovering a heterogeneous distribution of the distinct families of secretion systems within the genus. In this review, we describe the architecture and mode of action of bacterial type I to type VI secretion systems and the distribution and functions associated with these important nanoweapons within the Xanthomonas genus.

HpaXpm, a novel harpin of Xanthomonas phaseoli pv. manihotis, acts as an elicitor with high thermal stability, reduces disease, and promotes plant growth

BACKGROUND

Harpins are proteins secreted by the type III secretion system of Gram-negative bacteria during pathogen-plant interactions that can act as elicitors, stimulating defense and plant growth in many types of non-host plants. Harpin-treated plants have higher resistance, quality and yields and, therefore, harpin proteins may potentially have many valuable agricultural applications. Harpins are characterized by high thermal stability at 100 °C. However, it is unknown whether harpins are still active at temperatures above 100 °C or whether different temperatures affect the activity of the harpin protein in different ways. The mechanism responsible for the heat stability of harpins is also unknown.

RESULTS

We identified a novel harpin, HpaXpm, from the cassava blight bacteria Xanthomonas phaseoli pv. manihotis HNHK. The predicted secondary structure and 3-D structure indicated that the HpaXpm protein has two β-strand domains and two major α-helical domains located at the N- and C-terminal regions, respectively. A phylogenetic tree generated using the maximum likelihood method grouped HpaXpm in clade I of the Hpa1 group along with harpins produced by other Xanthomonas spp. (i.e., HpaG-Xag, HpaG-Xcm, Hpa1-Xac, and Hpa1Xm). Phenotypic assays showed that HpaXpm induced the hypersensitive response (HR), defense responses, and growth promotion in non-host plants more effectively than Hp1Xoo (X. oryzae pv. oryzae). Quantitative real-time PCR analysis indicated that HpaXpm proteins subjected to heat treatments at 100 °C, 150 °C, or 200 °C were still able to stimulate the expression of function-related genes (i.e., the HR marker genes Hin1 and Hsr203J, the defense-related gene NPR1, and the plant growth enhancement-related gene NtEXP6); however, the ability of heat-treated HpaXpm to induce HR was different at different temperatures.

CONCLUSIONS

These findings add a new member to the harpin family. HpaXpm is heat-stable up to 200 °C and is able to stimulate powerful beneficial biological functions that could potentially be more valuable for agricultural applications than those stimulated by Hpa1Xoo. We hypothesize that the extreme heat resistance of HpaXpm is because the structure of harpin is very stable and, therefore, the HpaXpm structure is less affected by temperature.

Type VI secretion system is not required for virulence on rice but for inter-bacterial competition in Xanthomonas oryzae pv. oryzicola

The type VI secretion system (T6SS), a multifunctional protein secretion device, plays very important roles in bacterial killing and/or virulence to eukaryotic cells. Although T6SS genes have been found in many Xanthomonas species, the biological function of T6SSs has not been elucidated in most xanthomonads. In this study, we identified two phylogenetically distinct T6SS clusters, T6SS1 and T6SS2, in a newly sequenced Chinese strain GX01 of Xanthomonas oryzea pv. oryzicola (Xoc) which causes bacterial leaf streak (BLS) of rice (Oryza sativa L.). Mutational assays demonstrated that T6SS1 and T6SS2 are not required for the virulence of Xoc GX01 on rice. Nevertheless, we found that T6SS2, but not T6SS1, played an important role in bacterial killing. Transcription and secretion analysis revealed that hcp2 gene is actively expressed and that Hcp2 protein is secreted via T6SS. Moreover, several candidate T6SS effectors were predicted by bioinformatics analysis that might play a role in the antibacterial activity of Xoc. This is the first report to investigate the type VI secretion system in Xanthomonas oryzae. We speculate that Xoc T6SS2 might play an important role in inter-bacterial competition, allowing this plant pathogen to gain niche advantage by killing other bacteria.

Transposons and pathogenicity in Xanthomonas: acquisition of murein lytic transglycosylases by TnXax1 enhances Xanthomonas citri subsp. citri 306 virulence and fitness

Xanthomonas citri subsp. citri 306 (XccA) is the causal agent of type A citrus canker (CC), one of the most significant citriculture diseases. Murein lytic transglycosylases (LT), potentially involved in XccA pathogenicity, are enzymes responsible for peptidoglycan structure assembly, remodeling and degradation. They directly impact cell wall expansion during bacterial growth, septum division allowing cell separation, cell wall remodeling allowing flagellar assembly, bacterial conjugation, muropeptide recycling, and secretion system assembly, in particular the Type 3 Secretion System involved in bacterial virulence, which play a fundamental role in XccA pathogenicity. Information about the XccA LT arsenal is patchy: little is known about family diversity, their exact role or their connection to virulence in this bacterium. Among the LTs with possible involvement in virulence, two paralogue open reading frames (ORFs) (one on the chromosome and one in plasmid pXAC64) are passenger genes of the Tn3 family transposon TnXax1, known to play a significant role in the evolution and emergence of pathogenicity in Xanthomonadales and to carry a variety of virulence determinants. This study addresses LT diversity in the XccA genome and examines the role of plasmid and chromosomal TnXax1 LT passenger genes using site-directed deletion mutagenesis and functional characterization. We identified 13 XccA LTs: 12 belong to families 1A, 1B, 1C, 1D (two copies), 1F, 1G, 3A, 3B (two copies), 5A, 6A and one which is non-categorized. The non-categorized LT is exclusive to the Xanthomonas genus and related to the 3B family but contains an additional domain linked to carbohydrate metabolism. The categorized LTs are probably involved in cell wall remodeling to allow insertion of type 3, 4 and 6 secretion systems, flagellum assembly, division and recycling of cell wall and degradation and control of peptidoglycan production. The TnXax1 passenger LT genes (3B family) are not essential to XccA or for CC development but are implicated in peptidoglycan metabolism, directly impacting bacterial fitness and CC symptom enhancement in susceptible hosts (e.g., Citrus sinensis). This underlines the role of TnXax1 as a virulence and pathogenicity-propagating agent in XccA and suggests that LT acquisition by horizontal gene transfer mediated by TnXax1 may improve bacterial fitness, conferring adaptive advantages to the plant-pathogen interaction process.

Genome-Wide Screening for Novel Candidate Virulence Related Response Regulator Genes in Xanthomonas oryzae pv. oryzicola

Two-component regulatory system (TCS), a major type of cellular signal transduction system, is widely used by bacteria to adapt to different conditions and to colonize certain ecological niches in response to environmental stimuli. TCSs are of distinct functional diversity, genetic diversity, and species specificity (pathovar specificity, even strain specificity) across bacterial groups. Although TCSs have been demonstrated to be crucial to the virulence of Xanthomonas, only a few researches have been reported about the studies of TCSs in Xanthomonas oryzae pathovar oryzicola (hereafter Xoc), the pathogen of rice bacterial streak disease. In the genome of Xoc strain GX01, it has been annotated 110 TCSs genes encoding 54 response regulators (RRs), 36 orthodox histidine kinase (HKs) and 20 hybrid histidine kinase (HyHKs). To evaluate the involvement of TCSs in the stress adaptation and virulence of Xoc, we mutated 50 annotated RR genes in Xoc GX01 by homologous vector integration mutagenesis and assessed their phenotypes in given conditions and tested their virulence on host rice. 17 RR genes were identified to be likely involved in virulence of Xoc, of which 10 RR genes are novel virulence genes in Xanthomonas, including three novel virulence genes for bacteria. Of the novel candidate virulence genes, some of which may be involved in the general stress adaptation, exopolysaccharide production, extracellular protease secretion and swarming motility of Xoc. Our results will facilitate further studies on revealing the biological functions of TCS genes in this phytopathogenic bacterium.

Identification of new protein-coding genes with a potential role in the virulence of the plant pathogen Xanthomonas euvesicatoria

Background

Bacteria of the genus Xanthomonas are economically important plant pathogens. Pathogenicity of Xanthomonas spp. depends on the type III-secretion system and additional virulence determinants. The number of sequenced Xanthomonas genomes increases rapidly, however, accurate annotation of these genomes is difficult, because it relies on gene prediction programs. In this study, we used a mass-spectrometry (MS)-based approach to identify the proteome of Xanthomonas euvesicatoria (Xe) strain 85–10 also known as X. campestris pv. vesicatoria, a well-studied member of plant-pathogenic Xanthomonadaceae.

Results

Using different culture conditions, MS-datasets were searched against a six-frame-translated genome database of Xe. In total, we identified 2588 proteins covering 55% of the Xe genome, including 764 hitherto hypothetical proteins. Our proteogenomic approach identified 30 new protein-coding genes and allowed correction of the N-termini of 50 protein-coding genes. For five novel and two N-terminally corrected genes the corresponding proteins were confirmed by immunoblot. Furthermore, our data indicate that two putative type VI-secretion systems encoded in Xe play no role in bacterial virulence which was experimentally confirmed.

Conclusions

The discovery and re-annotation of numerous genes in the genome of Xe shows that also a well-annotated genome can be improved. Additionally, our proteogenomic analyses validates “hypothetical” proteins and will improve annotation of Xanthomonadaceae genomes, providing a solid basis for further studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4041-7) contains supplementary material, which is available to authorized users.

Complete sequence and comparative genomic analysis of eight native Pseudomonas syringae plasmids belonging to the pPT23A family

BACKGROUND

The pPT23A family of plasmids appears to be indigenous to the plant pathogen Pseudomonas syringae and these plasmids are widely distributed and widely transferred among pathovars of P. syringae and related species. pPT23A-family plasmids (PFPs) are sources of accessory genes for their hosts that can include genes important for virulence and epiphytic colonization of plant leaf surfaces. The occurrence of repeated sequences including duplicated insertion sequences on PFPs has made obtaining closed plasmid genome sequences difficult. Therefore, our objective was to obtain complete genome sequences from PFPs from divergent P. syringae pathovars and also from strains of P. syringae pv. syringae isolated from different hosts.

RESULTS

The eight plasmids sequenced ranged in length from 61.6 to 73.8 kb and encoded from 65 to 83 annotated orfs. Virulence genes including type III secretion system effectors were encoded on two plasmids, and one of these, pPt0893-29 from P. syringae pv. tabaci, encoded a wide variety of putative virulence determinants. The PFPs from P. syringae pv. syringae mostly encoded genes of importance to ecological fitness including the rulAB determinant conferring tolerance to ultraviolet radiation. Heavy metal resistance genes encoding resistance to copper and arsenic were also present in a few plasmids. The discovery of part of the chromosomal genomic island GI6 from P. syringae pv. syringae B728a in two PFPs from two P. syringae pv. syringae hosts is further evidence of past intergenetic transfers between plasmid and chromosomal DNA. Phylogenetic analyses also revealed new subgroups of the pPT23A plasmid family and confirmed that plasmid phylogeny is incongruent with P. syringae pathovar or host of isolation. In addition, conserved genes among seven sequenced plasmids within the same phylogenetic group were limited to plasmid-specific functions including maintenance and transfer functions.

CONCLUSIONS

Our sequence analysis further revealed that PFPs from P. syringae encode suites of accessory genes that are selected at species (universal distribution), pathovar (interpathovar distribution), and population levels (intrapathovar distribution). The conservation of type IV secretion systems encoding conjugation functions also presumably contributes to the distribution of these plasmids within P. syringae populations.

Adaptation of genetically monomorphic bacteria: evolution of copper resistance through multiple horizontal gene transfers of complex and versatile mobile genetic elements

Copper-based antimicrobial compounds are widely used to control plant bacterial pathogens. Pathogens have adapted in response to this selective pressure. Xanthomonas citri pv. citri, a major citrus pathogen causing Asiatic citrus canker, was first reported to carry plasmid-encoded copper resistance in Argentina. This phenotype was conferred by the copLAB gene system. The emergence of resistant strains has since been reported in Réunion and Martinique. Using microsatellite-based genotyping and copLAB PCR, we demonstrated that the genetic structure of the copper-resistant strains from these three regions was made up of two distant clusters and varied for the detection of copLAB amplicons. In order to investigate this pattern more closely, we sequenced six copper-resistant X. citri pv. citri strains from Argentina, Martinique and Réunion, together with reference copper-resistant Xanthomonas and Stenotrophomonas strains using long-read sequencing technology. Genes involved in copper resistance were found to be strain dependent with the novel identification in X. citri pv. citri of copABCD and a cus heavy metal efflux resistance-nodulation-division system. The genes providing the adaptive trait were part of a mobile genetic element similar to Tn3-like transposons and included in a conjugative plasmid. This indicates the system's great versatility. The mining of all available bacterial genomes suggested that, within the bacterial community, the spread of copper resistance associated with mobile elements and their plasmid environments was primarily restricted to the Xanthomonadaceae family.