PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2017; 12(3): e0172295.
Published online 2017 March 20. doi:  10.1371/journal.pone.0172295
PMCID: PMC5358740

Comparative genomics of Clavibacter michiganensis subspecies, pathogens of important agricultural crops

Shihui Yang, Editor

Abstract

Subspecies of Clavibacter michiganensis are important phytobacterial pathogens causing devastating diseases in several agricultural crops. The genome organizations of these pathogens are poorly understood. Here, the complete genomes of 5 subspecies (C. michiganensis subsp. michiganensis, Cmi; C. michiganensis subsp. sepedonicus, Cms; C. michiganensis subsp. nebraskensis, Cmn; C. michiganensis subsp. insidiosus, Cmi and C. michiganensis subsp. capsici, Cmc) were analyzed. This study assessed the taxonomic position of the subspecies based on 16S rRNA and genome-based DNA homology and concludes that there is ample evidence to elevate some of the subspecies to species-level. Comparative genomics analysis indicated distinct genomic features evident on the DNA structural atlases and annotation features. Based on orthologous gene analysis, about 2300 CDSs are shared across all the subspecies; and Cms showed the highest number of subspecies-specific CDS, most of which are mobile elements suggesting that Cms could be more prone to translocation of foreign genes. Cms and Cmi had the highest number of pseudogenes, an indication of potential degenerating genomes. The stress response factors that may be involved in cold/heat shock, detoxification, oxidative stress, osmoregulation, and carbon utilization are outlined. For example, the wco-cluster encoding for extracellular polysaccharide II is highly conserved while the sucrose-6-phosphate hydrolase that catalyzes the hydrolysis of sucrose-6-phosphate yielding glucose-6-phosphate and fructose is highly divergent. A unique second form of the enzyme is only present in Cmn NCPPB 2581. Also, twenty-eight plasmid-borne CDSs in the other subspecies were found to have homologues in the chromosomal genome of Cmn which is known not to carry plasmids. These CDSs include pathogenesis-related factors such as Endocellulases E1 and Beta-glucosidase. The results presented here provide an insight of the functional organization of the genomes of five core C. michiganensis subspecies, enabling a better understanding of these phytobacteria.

Introduction

Members of the species Clavibacter michiganensis (Smith 1910) are gram-positive bacteria belonging to the family Microbacteriaceae, and consist of five core subspecies. The cells are (i) rods of coryneform morphology, (ii) having B2γ-type cell wall peptidoglycan with the diaminobutyric acid MK-9 as the predominant menaquinone, (iii) phosphatidyglycerol and diphosphatidyglycerol as the basic polar lipids, and (iv) a high GC content of 72–74 mol% [1,2]. All of the subspecies are plant pathogens of important agricultural crops (attacking members of the Solanaceae, Poaceae, and Leguminosae). Given the high level economic threat that they can cause, four of these subspecies are categorized as quarantine phytosanitary organisms [3]. They cause diseases of tomato (C. michiganensis subsp. michiganensis =, Cmm), potato (C. michiganensis subsp. sepedonicus =, Cms), alfalfa (C. michiganensis subsp. insidiosus =, Cmi), corn (C. michiganensis subsp. nebraskensis =, Cmn) and pepper (C. michiganensis subsp. capsici =, Cmc). Latent systemic infections of the xylem can be caused by all subspecies. All subspecies have been reported to invade seeds, and seems to poorly survive in soil [4,5]. In addition, they may have an epiphytic saprobic mode [6].

The genome organizations of subspecies of Clavibacter michiganensis are poorly understood. Next-generation technologies have revolutionized genome sequencing and as such the number of bacterial genomes available for analysis is expanding rapidly [7,8], leading to the generation of complete chromosomal and plasmid genomes of representatives strains of five subspecies (Cmm, Cms, Cmi, Cmn and Cmc) of C. michiganensis. Detailed analyses of the genomes of Cmm and Cms identified new sets of pathogenicity-related genes [9,10]. In Cmm and Cms, plasmid-borne virulence factors have been implicated in disease induction while chromosomally encoded genes are involved in successful host colonization [11]. In Cmm, a 129-kb low G+C region ((chp/tomA) near the origin of replication was considered essential for pathogenicity [10]. For example, individual genes found in this region, such as serine proteases, are necessary for effective colonization of tomato [10]. The serine protease-encoding pat-1 gene and cellulase-encoding celA gene in Cmm are directly implicated in pathogenicity [12]. An intact orthologue occurs in Cms. However, celB, a second cellulase gene, on the genome of both subspecies, is deactivated by a nonsense mutation in Cms [9]. It is unclear whether similar or novel regions exist in the genomes of Cmn, Cmi and Cmc. The complete chromosomal genome sequences of Cmn strain NCPPB 2581 (K.H. Gartemann, GenBank accession # HE614873), Cmi strain R1-1 [13] and Cmc strain PF008 [14] were published. The Cmi genome carries 3 plasmids while that of Cmc has two plasmids which might possess similar virulence factors. The genome of Cmn 2581 is not known to carry plasmids (Gartemann, per. Comm). Plasmids are reported not to be required for the pathogenicity of Cmn since most strains isolated do not carry a plasmid [11,15]. As such, it is suggested that the virulence mechanisms might be different from those reported for Cmm or Cms [16]. Since the genome of Cmn 2581 does not carry any plasmids, it can be hypothesized that the disease-inducing virulence factors are also chromosomally encoded alongside genes involved in successful host colonization. However, Clavibacter michiganensis subspecies harboring plasmid-borne disease-inducing virulence factors on the chromosome is yet to be reported.

The goals of this study were (i) to assess the taxonomic position of the subspecies based on 16S rRNA and genome-based DNA-DNAhomology; ii) to perform a comprehensive comparison of genomes of Cmm, Cms, Cmc, Cmi and Cmn using DNA structural and annotation features; (iii) to identify some of the genes involved in survival capacity and carbon utilization; and (iv) to assess whether some of the disease-inducing plasmid-borne virulence factors are present on the chromosomal genome of Cmn strain NCPPB 2581. Analyses of DNA structural features of complete genomes can pinpoint genomic regions that are sites of certain genes and elements involved in significant biological processes. Analyzing genome sequences can confer a wide range of new knowledge [17,18] useful in highlighting species and subspecies diversity that would not be otherwise possible [19].These will enable a better understanding of the host-specificity and pathogenicity of the subspecies of C. michiganensis and identify evolutionary genomic events associated with subspeciation [9]. The results presented here suggest that most of the subspecies could be distinct species. Comparative genomics revealed that the wco-cluster involved in extracellular polysaccharide II production is conserved within the subspecies while the sucrose-6-phosphate hydrolase is not; and outlined genes that may be implicated in stress responses. Finally, the data also show that some plasmid-borne genes in Cmm, Cms, Cmi and Cmc are chromosomally encoded in Cmn, known to not carry plasmids.

Materials and methods

Genome downloads and annotation

Whole-genome data of the five C. michiganensis subspecies were downloaded from GenBank [20] at NCBI, www.ncbi.nlm.nih.gov/genome/browser. NCBI GenBank International Nucleotide Sequence Database Collaboration (INSDC) or Whole-genome-sequence (WGS) numbers was used, respectively, to download each genome in the NCBI GenBank format using the getgbk.pl script as implemented in CMG-Biotools [19]. Genome sequences were extracted from GenBank files and saved in FASTA format using the saco_convert script [21]. The complete genomes of the 5 subspecies were submitted to the RAST web-based annotation system [22] and PATRIC [23] followed by manual curation.

Basic characterization of genomes

16S rRNA and gyrB-recA-rpoB phylogenies of Clavibacter michiganensis subspecies were implemented in MEGA7 [24] using neighbor-joining method with Kimura 2-parameter and Jukes-Cantor models respectively. Branch robustness was evaluated using 1000 bootstrap replicates. genome-sequence-based digital DNA-DNA hybridization (dDDH; [25]) and MUMmer-based average nucleotide identity (ANIm;[26]) were employed to assess the taxonomic position of strains relative to the closest taxon, Rathayibacter tritici NCPPB 1953 (GenBank # CP015515). The dDDH values were calculated using the genome-to-genome distance calculator (GGDC) Version 2.1 (http://ggdc.dsmz.de; [25]). ANIm similarity values were computed as described by Kurtz et al. [26] and implemented in JSpecies [27].

Genome comparison and analysis

The structural DNA atlases were generated from complete genomes as implemented in CMG-Biotools [19,28] to show the average and standard deviation of percent AT, GC skew, global repeats, intrinsic curvature and stacking energy. Each of the parameters are computed independently through a pipeline and outputted in a circular plot, an atlas [17].

Proteome comparisons

The comparison of proteomes was implemented using PATRIC web service [23]and CMG-Biotool [19]. PATRIC was executed using default parameters. For CMG-Biotool, a blastmatrix was generated using an XML formatted input file created by makebmdest [19]. A pairwise proteome comparison using BLAST [29] was used to generate a BLAST matrix. Protein sequences were compared to each other. Two sequences are similar and collected in the same ‘‘protein family’ if the BLAST hit had at least 50% identical matches in the alignment and the length of the alignment is 50% of the longest gene in the comparison. For the comparison of two genomes, single linkage is used to build protein families. Paralogs within a proteome are also evaluated and outputted at the bottom row of the matrix. Also, the Protein Family Sorter tool of PATRIC [23] was used to examine the distribution of specific gene families, known as FIGFams, across the different genomes. Analysis of orthologous clusters was also performed using the FastOrtho (http://enews.patricbrc.org/fastortho/), a faster reimplementation of OrthoMCL [30] with default parameters (e-vlaue of 1e-5 and inflation value of 1.5).

Results

Verification of strain identity and genomic relationship

Since Cmi strain R1-1 and Cmm strain NCPPB 382 were not type strains, their identities were verified. The 16S rDNA extracted of Cmi R1-1 and NCPPB 382 genomes were compared to those of their corresponding type strains by BLAST and phylogenetic analysis. Nucleotide BLAST searches (http://blast.ncbi.nlm.nih.gov/Blast.cgi) of the GenBank database showed that the 16S rDNA sequences of both strains exhibited more than 99% nucleotide identities to their respective type strains, LMG 3663T (U09761) and DSM 46364T (X77435). Seventeen 16S rDNA sequences from different subspecies and closely related genera (Rathayibacter and Leifsonia) were selected to infer a phylogenetic tree that showed strains Cmi R1-1 and Cmm NCPPB 382 clustered perfectly with their respective type strains (S1 Fig).

Genome similarity analysis using dDDH and ANIm showed values ranging from 39.1 to 60% and 90.75–95.25% respectively (Table 1). All the dDDH values are below (70%) the proposed cut-off species boundary. Highest dDDH homology (60%) was between Cmi and Cmn and the lowest was was Cms and Cmc. Similar trend was observed for ANIm values (cut-off = 95%) with the exception of Cmi-Cmn value that was 95.2% (Table 1). A well-supported gyrB-recA-rpoB phylogeny (Fig 1) of C. michiganensis subspecies is in agreement with dDDH and ANIm results.

Table 1
Chromosomal genome similarities between Clavibacter michiganensis subspecies based on genome-to-genome digital DNA-DNA Hybridization (dDDH; lower diagonal) and MuMmer-based Average Nucleotide Identity (ANIm; upper diagonal).
Fig 1
gyrB-recA-rpoB phylogeny of core Clavibacter michiganensis subspecies inferred using the neighbor-joining method and conducted in MEGA7[24].

Summary statistics and general features of the genomes

The basic statistics and general features of the five C. michiganensis subspecies genomes are shown in Table 2. The genome sizes ranged from 3.06 (Cmn) to 3.41 Mb (Cmi). All the subspecies genomes possess 2 or 3 plasmids except Cmn that has no plasmid. High G+C content (72.42–73.19%) is characteristic of actinomycetes. The number of protein-coding genes with function is between 2201 (Cmn) and 2341 (Cmi) with 18 to 114 pseudogenes (Table 2). The atlases of these C. michiganensis subspecies visually represent structural properties of the genomic DNA molecule such as intrinsic curvature, stacking energy, position preference, and inverted and direct repeats. S1 Fig. shows DNA structures of the five genomes including the locations of rRNA operons, inverted and direct repeats as well as strongly curved regions (high stacking energy and DNA intrinsic curvature) having genes that might be involved in a functionally specific DNA structure. The genome atlas of Cms ATCC 33113 exhibits, at least, 31 inverted repeat regions while Cmi R1-1 has about 14 (S2 Fig). The other subspecies lacked any visible inverted repeats. All the subspecies have three or more direct repeats, with Cms ATCC 33113 having the highest number of direct repeats (S2 Fig).

Table 2
General genome features of the five core Clavibacter michiganensis subspecies.

Comparative genomics

Comparison of the functional categories among the five subspecies of C. michiganensis shows that the highest number of CDSs that are involved in carbohydrate metabolism, while none is involved in photosynthesis (Fig 2). Pairwise proteome comparisons using the BLAST matrix [19] between the genomes showed similarity ranging from 66.1% to 74.6% (S3 Fig). The genome of Cmm NCPPB 382 has the highest similarity (75.3%) to that of Cmn NCPPB 2581 (S3 Fig). Cms exhibited the lowest proteome similarities (66.4–68.8%) with the other C. michiganensis subspecies (S3 Fig). To identify conserved and subspecies-specific CDSs, pan-genome analyses including orthologous group classification and orthologous relationship were performed. Orthologous relationships were determined using the FastOrtho method. All the CDSs of the five subspecies were clustered into 3,155 orthologous groups with 2,274 conserved groups. The number of conserved protein-coding sequences is relatively similar across the different subspecies (Table 2). Cmn NCPPB 2581 has the lowest number of subspecies-specific CDSs while Cms has the highest number (Table 2). PATRIC proteome comparison tool was used to compare the genomes of the five C. michiganensis subspecies. An overview of the conserved (blue arrow) and specific (brown arrow) genomic regions are given in Fig 3. Also, some of the plasmid-borne CDSs showed homologies to chromosomal genome of Cmn NCPPB 2581 known to not carry plasmids (Fig 3; square bracket). Twenty-eight plasmid-borne CDSs are present in the chromosomal genome of Cmn which include pathogenesis-related factors such as Endocellulases E1 and Beta-glucosidase (Fig 4). At least 75 CDSs related to stress response were identified in the genomes after analysis and comparison (S1 Table). These include CDSs involved in oxidative and osmotic stresses, cold and heat shock, and resistance to antibiotics and toxic compounds.

Fig 2
Comparison of functional categories among five subspecies of Clavibacter michiganensis.
Fig 3
Proteomes comparison of the genomes of the five Clavibacter michiganensis subspecies.
Fig 4
The chromosome of Clavibacter michiganensis subsp. nebraskensis NCPPB 2581 (Cmn) contains protein families with homologies in at least one of the plasmids of the other subspecies (Cmm, Cms, Cmc and Cmi).

Discussion

This study compared, for the first time, the complete genomes of five C. michiganensis subspecies. Of the 5 strains analysed, two (R1-1 and NCPPB 382) are not type strains. However, based on 16S rDNA BLAST and phylogenetic these strains were confirmed to belong to the same taxonomic positions as their corresponding type strains. Genome comparisons of the subspecies based on dDDH and ANIm showed values that are significantly below the cut-off threshold for species delineation, suggesting a higher taxonomic position (species-level) for these bacteria. A formal taxonomic study will provide a better insight.

Comparative genomic analysis of genomes showed that Cmn has the smallest genome, resulting in the fewest number of protein-coding genes, suggesting that colonizing and living in corn leaf tissues requires relatively few genes. Proteome comparison revealed that the Cms has the lowest similarity to the other C. michiganensis subspecies, suggesting that Cms is a more divergent probably linked to its soil niche, a more complex environment. Also, the Cms ATCC 33113 genome showed highest number of direct repeats most of which are mobile elements constituting most of the subspecies-specific protein-coding genes. Direct repeats play a significant role in the diversification of Helicobacter pylori DNA [31,32]. Wide-ranging repetitive DNA could facilitate the plasticity of a prokaryotic genome [31], suggesting that the genome of Cms ATCC 33113 could be more prone to translocation of foreign genes than the other subspecies. Also, the genomes of Cms and Cmi had the highest number of non-functional pseudogenes which might reduce the coding capacity of these strains, suggesting possible degeneration of the genome [9]. This process is often associated with new niche adaptation by a bacterial species, making certain gene expendables [9,10].

The genomic DNA atlases also revealed differences in intrinsic curvatures. High curvature and stacking energy regions, for example, in Cmm NCPPB 382 (S1 Fig; brown arrows) indicate strongly curved regions that might be involved in specific biological function. Curved DNA portions seems to have highly expressed genes that are modulated by histone-like proteins [19]. The rRNA operons are associated with regions of high curvature, average stacking energy and low position preference in all the chromosomal genomes of the subspecies. DNA curvature plays a significant role in several biologically vital processes, including recombination [33], DNA replication [34], and positioning of nucleosome [35].

Comparisons of functional categories among the genomes of the five subspecies showed that the number genes implicated in carbohydrate metabolism and transport (Fig 2, category AA) were highest compared to the other categories within each genome. This suggests that carbohydrate metabolism is a key factor to the survival of these subspecies and could be involved in plant-pathogenic interaction. For example, in planta, genes within the wco-cluster involved in sugar metabolism were up-regulated in Cmm in late infection stages suggesting potential involvement in pathogenicity [36]. The functions of genes in this cluster include chitinases, putative glycosyltransferase, glycoamylases and GumJ proteins. In a tomato plant study [36] in Cmm, the CMM_0824 locus encoding for glycosyltransferase (wcoF) showed highest up-regulated value. Also, the GumJ protein contributes to the formation of biofilm and cells adhesion to host surfaces [3740]. Genome-wide comparison of the five subspecies showed that the genes within this cluster are generally conserved (90–99%). Seventeen CDSs were identified in all the genomes. Three CDSs identified on the Cmm genome as wcoA, wcoB and wcoP had low homology in the other genomes. wcoA, a chitinase, in Cmm had only 67.3% similarity to a potential homologous gene in Cmc. A hypothetical protein (wcoB) in Cmm had low similarities to CDSs in Cms (85%), Cmc (66.4%) and Cmn (68.7%). A transcriptional regulator of the MarR family (wcoP) present in Cmm showed about 86.8% and 36.6% in Cmn and Cmc, respectively. Given their up-regulation in planta, these genes may play an important role in utilizing plant derived nutrients.

Sucrose is a naturally abundant carbohydrate found in several plants and plant parts (Reid and Abratt, 2005). A CDS encoding for sucrose phosphate synthases associated with sucrose biosynthesis was identified in all the subspecies and showed about 93.0% homology to locus CMM_0494 found in Cmm. Two CDSs associated with sucrose catabolism were identified but only one is present in all the subspecies. Sucrose phosphorylase (S1S5 Datasets), an important enzyme that converts sucrose to D-fructose and alpha-D-glucose-1-phosphate, is present in all the subspecies with about 87.8% homology to locus CMM_2523 found in Cmm. However, a sucrose-6-phosphate hydrolase (EC 3.2.1.26) found in Cmm (CMM_2780) was identified only in Cmn and Cms (CMS_0938) with a low homology of 36.4% and 49.8% respectively. A second CDS encoding another form of sucrose-6-phosphate hydrolase (EC 3.2.1.B3) is present only in Cmn, the pathogen of corn. Corn possesses a very active sucrose-6-phosphate biosynthetic system. Cytoplasmic sucrose-6-phosphate hydrolase catalyzes the hydrolysis of sucrose-6-phosphate yielding glucose-6-phosphate and fructose [41]. It is unclear why this high divergence among the subspecies especially its absence in Cmi and Cmc. It is possible that alternate pathways exit in Cmi and Cmc. In Streptococcus mutans, Tao et al. [42] indicated that other sugar transport including sucrose is done through the MSM (multiple sugar metabolism) systems.

The survival of bacteria in a given environment depends on the ability to respond to changes in oxidative stress. At least 22 CDSs involved in oxidative stress response were identified in all the subspecies (S1 Table). The CDSs coding for catalases (EC 1.11.1.6), superoxide dismutases (EC 1.15.1.1), and ferroxidases (EC 1.16.3.1) are conserved among the five genomes of the subspecies with homology of about 99%. Other CDSs found in all the subspecies include iron-binding ferritin-like antioxidant protein and alkyl hydroperoxidase reductase subunit C-like protein. In addition, all the Clavibacter subspecies genomes encode glutathione peroxidase (EC 1.11.1.9). Also, a CDS encoding for redox-sensitive activator (SoxR), an oxidative stress response protein; furB, a zinc uptake regulation protein (ZUR), and a transcriptional regulator of the FUR family are present in all the genomes studied.

All Clavibacter subspecies encode 5 CDSs involved in biosynthesis of mycothiol, an unusual thiol compound found in the Actinobacteria with important antioxidant and detoxification functions[43]. A CDS, mshA encodes N-acetylglucosamine transferase involved in the formation of GlcNAc-Ins; mshB encodes for deacetylase; mshC (ligase) catalyses the ligation of GlcN-Ins with a cysteine [44] followed by the acetylation of Cys-GlcN-Ins to form mycothiol. This acetylation process is catalysed by mshD, acetyltransferase [45]. The fifth CDS is the mycothiol S-conjungate amidase, Mca. Mca is involved in the cleavage of the amide bond of mycothiol S-conjugates of specific xenobiotics and alkylating agents producing mercapturic acid and GlcN-Ins excreted from the cell [43]. While Mca had a homology level of 96% among the subspecies, lower homology values (89.3–91.0%) were observed for genes involved in mycothiol biosynthesis.

In Actinobacteria, mycothiol biosynthesis is also implicated in arsenate resistance [46], a process that involves chemically reducing the toxic arsenate. The reduction of the product arseno-mycothiol is catalysed by mycoredoxin (EC 1.20.4.3) to mycothiol-mycoredoxin disulfide and arsenite followed by the formation of mycothione by a second mycothiol that recycles mycoredoxin. In the genomes of Clavibacter subspecies, CDSs linked to arsenic resistance are chromosomally and plasmid encoded except for Cmn 2581 where it is in the chromosome only. Two CDSs, arsB encoding arsenic efflux pump protein and arsC2 encoding arsenate-mycothiol transferase (EC 2.8.4.2) are present in all the chromosomes of the subspecies with high homology. Also, three CDSs encoding the arsenic transcriptional repressor (arsR) are present in the chromosome of all the genomes. In addition, one CDS of arsR is carried in the plasmids of all the subspecies except the Cmn which has no plasmid. The lack of plasmid in Cmn can suggest a low tolerance to arsenic. In the Staphylococcus [47,48]- or E. coli R773 or R46 [49,50]the plasmid-borne operons confer considerably high level of arsenic resistance than the chromosomal operon.

In addition to arsenic tolerance, the survival of bacteria in their respective ecological niches is dependent on their resistance to antibiotics and toxic compounds including metals such as selenium and copper. Bacteria have developed effective homeostasis and resistance systems in order to maintain the required functional amounts of these metals while detoxifying excesses. These complicated processes involve acquisition, sequestration, and efflux of metal ions [51]. Selenium occurs naturally in the Earth’s crust; and at low concentration it is essential for living organisms [52]. Under aerobic conditions, this trace element exists as selenite and selenate, and at high levels these salts can be toxic and mutagenic to bacteria [52,53]. High selenite-resistant bacterial strains like Ralstonia metallidurans CH34 possess the dedA gene that regulatesselenite uptake [52,53]. Three dedA genes encoding the putative selenite transport protein (DedA) including various polyols permease components of the ABC transporters are present in each of the Clavibacter subspecies, suggesting that members of the species C. michiganensis can detoxify environmental selenite/selenite.

In addition, Copper, an essential trace and redox-active element, serves as a cofactor for several enzymes. In aerobic cells, excess Cu metal ion can produce cytotoxic reactive oxygen species capable of damaging DNA, lipids and proteins[51]. A CDS that is chromosomally encoding Copper-translocating P-type ATPase (copA; EC 3.6.3.4), repressor CsoR of the copZA operon, and Copper (I) chaperone CopZ; two CDSs each encoding for Copper resistance protein CopC and conserved membrane protein in copper uptake (YcnI) are present in all the subspecies. In addition, all the genomes have one CDS encoding for copD (a Copper resistance protein) except for Cmc PF008 that has two CDSs encoding for CopD. It might be interesting to elucidate why Cmc PF008 has more than one copy of the copD. Other stress response factors found in all the genomes include sigma factors (RsbW, RsbV, SigB, RsbU), Hfl operon encoding the GTP-binding protein, bacterial hemoglobin-like protein (HbO). Each subspecies has a CDS for HbO.

Cold- and heat-shock responses enable bacteria to survive changes in environmental temperature [54]. The cold shock response is governed by the expression of RNA chaperones and ribosomal factors. Two cold-shock protein (cspA and cspC) genes were identified in each of the Clavibacter michiganensis subspecies. In Escherichia coli, cspC, reported previously to be a regulator of rpoS [55], is expressed at 37°C and involved in cell division [56,57]. Bacterial responses to heat shock include heat shock proteins (HSPs) that are encoded by transcriptional up-regulation of genes. Genome of all the subspecies have a dnaK gene that encodes for heat-shock protein GrpE, chaperone proteins DnaJ and DnaK, a transcriptional repressor of the dnaK operon (hspR), hrcA, a heat-inducible repressor of transcription, and other genes (e.g. smpB, encoding HSPs). With the exception of the chaperone grpE, cold- and heat-shock response proteins are conserved (homology of 96–99%) among the subspecies. Cmi and Cmn exhibited a 98% homology with the protein GrpE while both had only a 91% similarity to the other C. michiganensis subspecies.

This study assessed the taxonomic position of the subspecies based on 16S rRNA and genome-based DNA homology and concludes that there is ample evidence to perform a detailed analysis to elevate some of the subspecies to species-level. In addition, a detailed comparative genomics of the genomes of the subspecies indicated distinct genomic features evident on the DNA structural atlases and annotation features. Orthologous gene analysis revealed that the about 2300 CDSs are conserved across all the subspecies; and Cms showed the highest number of subspecies-specific CDS, most of which are mobile elements, suggesting that Cms could be more prone to translocation of foreign genes. In addition, Cms Cmi had the highest number of pseudogenes, an indication of potential degenerating genomes. This study also summarized some of the genetic factors encoded in these subspecies to survive under different stress conditions. The study outlined some of the stress response factors that may be involved in cold/heat shock, detoxification, oxidative stress, osmo-regulation and carbon utilization. In carbon utilization, the wco cluster encoding for extracellular polysaccharide II is highly conserved while the sucrose-6-phosphate hydrolase that catalyzes the hydrolysis of sucrose-6-phosphate yielding glucose-6-phosphate and fructose is highly diverged. It will be intriguing to elucidate why this gene is absent in Cmc and Cmi. The results presented here provide an insight of the functional organization of the genomes of five C. michiganensis subspecies and as such a better understanding of these phytobacteria.

Supporting information

S1 Dataset

Annotated features of Clavibacter michiganesis subsp. michiganensis strain NCPPB 382 generated by PATRIC [23].

(XLSX)

S2 Dataset

Annotated features of Clavibacter michiganesis subsp. sepedonicus strain ATCC 33113 generated by PATRIC [23].

(XLSX)

S3 Dataset

Annotated features of Clavibacter michiganesis subsp. nebraskensis strain NCPPB 2581 generated using PATRIC [23].

(XLSX)

S4 Dataset

Annotation features of Clavibacter michiganesis subsp. insidiosus strain R1-1 generated using PATRIC [23].

(XLSX)

S5 Dataset

Annotated features of Clavibacter michiganesis subsp. capsici strain PF008 generated by PATRIC [23].

(XLSX)

S1 Fig

16S rDNA phylogenetic tree of Clavibacter michiganensis strains inferred using the neighbor-joining method1 implemented in MEGA7 [24].

The optimal tree with the sum of branch length = 0.11186829 is shown. The values next to the branches are percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates). Bootstrap values greater than 50% are shown. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Kimura 2-parameter method. The analysis involved 17 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 1,250 positions in the final dataset. Taxa in bold are strains used in genome comparison that are not type strains and clustered perfectly with the corresponding type strains. The sequence accession numbers of the taxa are given in parentheses.

(PPTX)

S2 Fig

DNA structures of complete genomes of the five Clavibacter michiganensis subspecies based on genomic atlases.

Data of DNA, RNA and gene annotation are from the published GenBank entries. Each lane of the circular representation of the chromosome shows a different DNA feature. From innermost circle: size of genome (axis), percent AT (red = high AT), GC skew (blue = most G’s; orange arrows), inverted and direct repeats (color = repeats), position preference, stacking energy and intrinsic curvature. Dark brown arrows highlight areas of the genome with significantly different DNA structures than the remaining of the genome. Blue arrows shows the locations of rRNA operons as annotated in the GenBank file. Genome atlas was generated using CMG-Biotools [19] which calculates a numerical value for each nucleotide and saved in a file that is read by GeneWiz software. See “Materials and Methods” for details.

(PPTX)

S3 Fig

BLAST matrix of proteomes of five strains of Clavibacter michiganensis subspecies based on all against protein comparison to define homologs.

A hit is considered significant if 50%/50% (identity/length coverage) requirement between-proteomes is met. Paralogs (internal homology) are proteins within a genome matching the same 50–50 rule. Cms, C. michiganesis subsp. sepedonicus; Cmn, C. m. subsp. nebraskensis; Cmm, C. m. subsp. michiganensis; Cmi, C. m. subsp. insidiosus; C. m. subsp. capsici.

(PPTX)

S1 Table

Number of protein-coding sequences related to stress response survival capacity of the five subspecies of Clavibacter michiganensis.

(PPTX)

Acknowledgments

This study was funded by Agriculture and Agri-Food Canada through projects J-000409, AAFC/GRDI# J-000011 and J-000985 and by Manitoba Corn Growers Association project # CRADA AGR-10755. I am indebted to S. Miller and K. Seifert for reviewing the draft manuscript.

Funding Statement

This study was funded by the Agriculture and Agri-Food Canada through projects J-000409, AAFC/GRDI# J-000011 and J-000985 and by the Manitoba Corn Growers Association project # CRADA AGR-10755.

This paper was supported by the following grant(s):

http://dx.doi.org/10.13039/501100000040Agriculture and Agri-Food Canada J-000409 and J-000985 to James T. Tambong.
Manitoba Corn Growers Association CRADA AGR-10755 to James T. Tambong.

Data Availability

Data Availability

All relevant data are within the paper and its Supporting Information files.

References

1. Agarkova IV, Lambrecht PA, Vidaver AK (2011) Genetic diversity and population structure of Clavibacter michiganensis subsp. nebraskensis. Can J Microbiol 57: 366–374. doi: 10.1139/W11-016 [PubMed]
2. Davis MJ, Gillaspie AG, Vidaver AK, Harris R (1984) Clavibacter: a new genus containing some phytopathogenic coryneform bacteria, including Clavibacter xyli subsp. xyli sp. nov., and Clavibacter xyli subsp. cynodontis subsp. nov., pathogens that cause ratoon stunting of sugarcane and Bermuda grass stunting disease. Int J Syst Bacteriol 34: 107–117.
3. Jacques MA, Durand K, Orgeur G, Balidas S, Fricot C, Bonneau S, et al. (2012) Phylogenetic analysis and polyphasic characterization of Clavibacter michiganensis strains isolated from tomato seeds reveal that nonpathogenic strains are distinct from C. michiganensis subsp. michiganensis. Appl Environ Microbiol 78: 8388–8402. doi: 10.1128/AEM.02158-12 [PMC free article] [PubMed]
4. Brumbley SM, Petrasovits LA, Hermann SR, Young AJ, Croft BJ (2006) Recent advances in the molecular biology of Leifsonia xyli subsp. xyli, causal organism of ratoon stunting disease. Aust Plant Pathol 35681–35689.
5. Park YH, Suzuki K, Yim DG, Lee KC, Kim E, Yoon J, et al. (1993) Suprageneric classification of peptidoglycan group B actinomycetes by nucleotide sequencing of 5S ribosomal RNA. Antonie Van Leeuwenhoek 64: 307–313. [PubMed]
6. de Souza ML, Newcombe D, Alvey S, Crowley DE, Hay A, Sadowsky MJ, et al. (1998) Molecular basis of a bacterial consortium: interspecies catabolism of atrazine. Appl Environ Microbiol 64: 178–184. [PMC free article] [PubMed]
7. Achtman M, Wagner M (2008) Microbial diversity and the genetic nature of microbial species. Nat Rev Microbiol 6: 431–440. doi: 10.1038/nrmicro1872 [PubMed]
8. Jolley KA, Bliss CM, Bennett JS, Bratcher HB, Brehony C, Colle FM, et al. (2012) Ribosomal multilocus sequence typing: universal characterization of bacteria from domain to strain. Microbiology 158: 1005–1015. doi: 10.1099/mic.0.055459-0 [PMC free article] [PubMed]
9. Bentley SD, Corton C, Brown SE, Barron A, Clark L, Doggett J, et al. (2008) Genome of the actinomycete plant pathogen Clavibacter michiganensis subsp. sepedonicus suggests recent niche adaptation. J Bacteriol 190: 2150–2160. doi: 10.1128/JB.01598-07 [PMC free article] [PubMed]
10. Gartemann KH, Abt B, Bekel T, Burger A, Engemann J, Flugel M, et al. (2008) The genome sequence of the tomato-pathogenic actinomycete Clavibacter michiganensis subsp. michiganensis NCPPB382 reveals a large island involved in pathogenicity. J Bacteriol 190: 2138–2149. doi: 10.1128/JB.01595-07 [PMC free article] [PubMed]
11. Ahmad A, Mbofung GY, Acharya J, Schmidt CL, Robertson AE (2015) Characterization and Comparison of Clavibacter michiganensis subsp. nebraskensis strains recovered from Epiphytic and symptomatic infections of maize in Iowa. PLoS One 10: e0143553 doi: 10.1371/journal.pone.0143553 [PMC free article] [PubMed]
12. Gartemann KH, Kirchner O, Engemann J, Grafen I, Eichenlaub R, Burger A (2003) Clavibacter michiganensis subsp. michiganensis: first steps in the understanding of virulence of a Gram-positive phytopathogenic bacterium. J Biotechnol 106: 179–191. [PubMed]
13. Lu Y, Samac DA, Glazebrook J, Ishimaru CA (2015) Complete Genome Sequence of Clavibacter michiganensis subsp. insidiosus R1-1 Using PacBio Single-Molecule Real-Time Technology. Genome Announc 3. [PMC free article] [PubMed]
14. Bae C, Oh EJ, Lee HB, Kim BY, Oh CS (2015) Complete genome sequence of the cellulase-producing bacterium Clavibacter michiganensis PF008. J Biotechnol 214: 103–104. doi: 10.1016/j.jbiotec.2015.09.021 [PubMed]
15. Vidaver AK, Gross DC, Wysong DS, Doupnik JR (1981) Diversity of Corynebacterium nebraskense strains causing Goss’s bacterial wilt and blight of corn. Plant Disease 65: 480–483.
16. Eichenlaub R, Gartemann K.-H., Burger A (2006) Clavibacter michiganensis, a group of Gram-positive phytopathogenic bacteria In: Gnanamanickam SS, editor. Plant-associated bacteria. Dordrecht, The Netherlands: Springer; pp. 385–422.
17. Fleischmann RD, Adams MD, White O, Clayton RA, Kirkness EF, Kerlavage AR, et al. (1995) Whole-genome random sequencing and assembly of Haemophilus influenzae Rd. Science 269: 496–512. [PubMed]
18. Fraser CM, Gocayne JD, White O, Adams MD, Clayton RA, Flieschmann RD, et al. (1995) The minimal gene complement of Mycoplasma genitalium. Science 270: 397–403. [PubMed]
19. Vesth T, Lagesen K, Acar O, Ussery D (2013) CMG-biotools, a free workbench for basic comparative microbial genomics. PLoS One 8: e60120 doi: 10.1371/journal.pone.0060120 [PMC free article] [PubMed]
20. Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW (2011) GenBank. Nucleic Acids Res 39: D32–37. doi: 10.1093/nar/gkq1079 [PMC free article] [PubMed]
21. Jensen LJ, Knudsen S (1999) Automatic discovery of regulatory patterns in promoter regions based on whole cell expression data and functional annotation. Bioinformatics: 326–333. [PubMed]
22. Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. (2008) The RAST Server: rapid annotations using subsystems technology. BMC Genomics 9: 75 doi: 10.1186/1471-2164-9-75 [PMC free article] [PubMed]
23. Wattam AR, Abraham D, Dalay O, Disz TL, Driscoll T, Gabbard JL, et al. (2014) PATRIC, the bacterial bioinformatics database and analysis resource. Nucleic Acids Res 42: D581–591. doi: 10.1093/nar/gkt1099 [PMC free article] [PubMed]
24. Kumar S, Stecher G, Tamura K (2016) MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol 33: 1870–1874. doi: 10.1093/molbev/msw054 [PubMed]
25. Meier-Kolthoff JP, Auch AF, Klenk HP, Goker M (2013) Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinformatics 14: 60 doi: 10.1186/1471-2105-14-60 [PMC free article] [PubMed]
26. Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, et al. (2004) Versatile and open software for comparing large genomes. Genome Biol 5: R12 doi: 10.1186/gb-2004-5-2-r12 [PMC free article] [PubMed]
27. Richter M, Rossello-Mora R (2009) Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A 106: 19126–19131. doi: 10.1073/pnas.0906412106 [PubMed]
28. Tambong JT, Xu R, Adam Z, Cott M, Rose K, Reid LM, et al. (2015) Draft genome sequence of Clavibacter michiganensis subsp. nebraskensis strain DOAB 397, Isolated from an Infected field corn plant in Manitoba, Canada. Genome Announc 3 (4):e00768–15. doi: 10.1128/genomeA.00768-15 [PMC free article] [PubMed]
29. Binnewies TT, Hallin PF, Staerfeldt HH, Ussery DW (2005) Genome Update: proteome comparisons. Microbiology 151: 1–4. doi: 10.1099/mic.0.27760-0 [PubMed]
30. Li L, Stoeckert CJ Jr., Roos DS (2003) OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res 13: 2178–2189. doi: 10.1101/gr.1224503 [PubMed]
31. Aras RA, Kang J, Tschumi AI, Harasaki Y, Blaser MJ (2003) Extensive repetitive DNA facilitates prokaryotic genome plasticity. Proc Natl Acad Sci U S A 100: 13579–13584. doi: 10.1073/pnas.1735481100 [PubMed]
32. Lara-Ramirez EE, Segura-Cabrera A, Guo X, Yu G, Garcia-Perez CA, Rodriguez-Perez MA (2011) New implications on genomic adaptation derived from the Helicobacter pylori genome comparison. PLoS One 6: e17300 doi: 10.1371/journal.pone.0017300 [PMC free article] [PubMed]
33. Mazin A, Milot E, Devoret R, Chartrand P (1994) KIN17, a mouse nuclear protein, binds to bent DNA fragments that are found at illegitimate recombination junctions in mammalian cells. Mol Gen Genet 244: 435–438. [PubMed]
34. Ueguchi C, Kakeda M, Yamada H, Mizuno T (1994) An analogue of the DnaJ molecular chaperone in Escherichia coli. Proc Natl Acad Sci U S A 91: 1054–1058. [PubMed]
35. Kiyama R, Trifonov EN (2002) What positions nucleosomes?—A model. FEBS Lett 523: 7–11. [PubMed]
36. Flugel M, Becker A, Gartemann KH, Eichenlaub R (2012) Analysis of the interaction of Clavibacter michiganensis subsp. michiganensis with its host plant tomato by genome-wide expression profiling. J Biotechnol 160: 42–54. doi: 10.1016/j.jbiotec.2012.01.023 [PubMed]
37. Feil H, Feil WS, Lindow SE (2007) Contribution of fimbrial and afimbrial adhesins of Xylella fastidiosa to attachment to surfaces and virulence to grape. Phytopathology 97: 318–324. doi: 10.1094/PHYTO-97-3-0318 [PubMed]
38. Guilhabert MR, Kirkpatrick BC (2005) Identification of Xylella fastidiosa antivirulence genes: hemagglutinin adhesins contribute a biofilm maturation to X. fastidiosa and colonization and attenuate virulence. Mol Plant Microbe Interact 18: 856–868. doi: 10.1094/MPMI-18-0856 [PubMed]
39. Killiny N, Almeida RP (2009) Host structural carbohydrate induces vector transmission of a bacterial plant pathogen. Proc Natl Acad Sci U S A 106: 22416–22420. doi: 10.1073/pnas.0908562106 [PubMed]
40. Meng Y, Li Y, Galvani CD, Hao G, Turner JN, Burr TJ, et al. (2005) Upstream migration of Xylella fastidiosa via pilus-driven twitching motility. J Bacteriol 187: 5560–5567. doi: 10.1128/JB.187.16.5560-5567.2005 [PMC free article] [PubMed]
41. Wang B, Kuramitsu HK (2003) Control of enzyme IIscr and sucrose-6-phosphate hydrolase activities in Streptococcus mutans by transcriptional repressor ScrR binding to the cis-active determinants of the scr regulon. J Bacteriol 185: 5791–5799. doi: 10.1128/JB.185.19.5791-5799.2003 [PMC free article] [PubMed]
42. Tao L, Sutcliffe IC, Russell RR, Ferretti JJ (1993) Transport of sugars, including sucrose, by the msm transport system of Streptococcus mutans. J Dent Res 72: 1386–1390. doi: 10.1177/00220345930720100701 [PubMed]
43. Buchmeier N, Fahey RC (2006) The mshA gene encoding the glycosyltransferase of mycothiol biosynthesis is essential in Mycobacterium tuberculosis Erdman. FEMS Microbiol Lett 264: 74–79. doi: 10.1111/j.1574-6968.2006.00441.x [PubMed]
44. Sareen D, Steffek M, Newton GL, Fahey RC (2002) ATP-dependent L-cysteine:1D-myo-inosityl 2-amino-2-deoxy-alpha-D-glucopyranoside ligase, mycothiol biosynthesis enzyme MshC, is related to class I cysteinyl-tRNA synthetases. Biochemistry 41: 6885–6890. [PubMed]
45. Koledin T, Newton GL, Fahey RC (2002) Identification of the mycothiol synthase gene (mshD) encoding the acetyltransferase producing mycothiol in actinomycetes. Arch Microbiol 178: 331–337. doi: 10.1007/s00203-002-0462-y [PubMed]
46. Ordonez E, Van Belle K, Roos G, De Galan S, Letek M, Gil JA, et al. (2009) Arsenate reductase, mycothiol, and mycoredoxin concert thiol/disulfide exchange. J Biol Chem 284: 15107–15116. doi: 10.1074/jbc.M900877200 [PMC free article] [PubMed]
47. Ji G, Silver S (1992) Regulation and expression of the arsenic resistance operon from Staphylococcus aureus plasmid pI258. J Bacteriol 174: 3684–3694. [PMC free article] [PubMed]
48. Rosenstein R, Peschel A, Wieland B, Gotz F (1992) Expression and regulation of the antimonite, arsenite, and arsenate resistance operon of Staphylococcus xylosus plasmid pSX267. J Bacteriol 174: 3676–3683. [PMC free article] [PubMed]
49. Hedges RW, Baumberg S (1973) Resistance to arsenic compounds conferred by a plasmid transmissible between strains of Escherichia coli. J Bacteriol 115: 459–460. [PMC free article] [PubMed]
50. Mobley HL, Silver S, Porter FD, Rosen BP (1984) Homology among arsenate resistance determinants of R factors in Escherichia coli. Antimicrob Agents Chemother 25: 157–161. [PMC free article] [PubMed]
51. Nawapan S, Charoenlap N, Charoenwuttitam A, Saenkham P, Mongkolsuk S, Vattanaviboon P. (2009) Functional and expression analyses of the cop operon, required for copper resistance in Agrobacterium tumefaciens. J Bacteriol 191: 5159–5168. doi: 10.1128/JB.00384-09 [PMC free article] [PubMed]
52. Yao L, Du Q, Yao H, Chen X, Zhang Z, Xu S (2015) Roles of oxidative stress and endoplasmic reticulum stress in selenium deficiency-induced apoptosis in chicken liver. Biometals 28: 255–265. doi: 10.1007/s10534-014-9819-3 [PubMed]
53. Ledgham F, Quest B, Vallaeys T, Mergeay M, Coves J (2005) A probable link between the DedA protein and resistance to selenite. Res Microbiol 156: 367–374. doi: 10.1016/j.resmic.2004.11.003 [PubMed]
54. Ron EZ (2012) Bacterial Stress Response; Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thomson F, editors. Berlin: Springer.
55. Phadtare S, Inouye M (2001) Role of CspC and CspE in regulation of expression of RpoS and UspA, the stress response proteins in Escherichia coli. J Bacteriol 183: 1205–1214. doi: 10.1128/JB.183.4.1205-1214.2001 [PMC free article] [PubMed]
56. Hu KH, Liu E, Dean K, Gingras M, DeGraff W, Trun NJ (1996) Overproduction of three genes leads to camphor resistance and chromosome condensation in Escherichia coli. Genetics 143: 1521–1532. [PubMed]
57. Rath D, Jawali N (2006) Loss of expression of cspC, a cold shock family gene, confers a gain of fitness in Escherichia coli K-12 strains. J Bacteriol 188: 6780–6785. doi: 10.1128/JB.00471-06 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science