International Journal of Genomics
Volume 2018 (2018), Article ID 1062716, 10 pages
https://doi.org/10.1155/2018/1062716

Research Article

Insights into the Genome Sequence of Chromobacterium amazonense Isolated from a Tropical Freshwater Lake

Alexandre Bueno Santos,¹ Patrícia Silva Costa,¹ Anderson Oliveira do Carmo,¹ Gabriel da Rocha Fernandes,² Larissa Lopes Silva Scholte,² Jeronimo Ruiz,² Evanguedes Kalapothakis,¹ Edmar Chartone-Souza,¹ and Andréa Maria Amaral Nascimento¹

¹Departamento de Biologia Geral, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil
²Centro de Pesquisas René Rachou, FIOCRUZ, Belo Horizonte, MG, Brazil

Correspondence should be addressed to Andréa Maria Amaral Nascimento; amaral@ufmg.br

Received 5 December 2017; Revised 30 March 2018; Accepted 19 April 2018; Published 20 May 2018

Academic Editor: Wilfred van IJcken

Copyright © 2018 Alexandre Bueno Santos et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Members of the genus Chromobacterium have been isolated from geographically diverse ecosystems and exhibit considerable metabolic flexibility, as well as biotechnological and pathogenic properties in some species. This study reports the draft assembly and detailed sequence analysis of Chromobacterium amazonense strain 56AF. The de novo-assembled genome is 4,556,707 bp in size and contains 4294 protein-coding and 95 RNA genes, including 88 tRNA, six rRNA, and one tmRNA operon. A repertoire of genes implicated in virulence, for example, hemolysin, hemolytic enterotoxins, colicin V, lytic proteins, and Nudix hydrolases, is present. The genome also contains a collection of genes of biotechnological interest, including esterases, lipase, auxins, chitinases, phytoene synthase and phytoene desaturase, polyhydroxyalkanoates, violacein, plastocyanin/azurin, and detoxifying compounds. Importantly, unlike other Chromobacterium species, the 56AF genome contains genes for pore-forming toxin alpha-hemolysin, a type IV secretion system, among others. The analysis of the C. amazonense strain 56AF genome reveals the versatility, adaptability, and biotechnological potential of this bacterium. This study provides molecular information that may pave the way for further comparative genomics and functional studies involving Chromobacterium-related isolates and improves our understanding of the global genomic diversity of Chromobacterium species.

1. Introduction

The genus Chromobacterium, belonging to the class Betaproteobacteria and family Chromobacteriaceae (formerly Neisseriaceae) [1], was originally proposed in 1881 by Bergonzini [2]. Until 2007, the only validated species within the genus Chromobacterium was C. violaceum, the type species of the genus. The current taxonomy of Chromobacterium consists of 11 recognized species: C. violaceum [2], C. subtsugae [3], C. aquaticum [4], C. haemolyticum [5], C. piscinae, C. pseudoviolaceum [6], C. vaccinii [7], C. amazonense [8], C. rhizoryzae [9], C. alkanivorans [10], and C. sphagni [11]. It is evident that, since 2007, there has been a rapid taxonomic expansion of the Chromobacterium genus.

Species in this genus have been collected from geographically diverse ecosystems. The wide spread of Chromobacterium throughout a variety of environments, such as soil, water, and plant from tropical and subtropical regions, is due to its considerable metabolic flexibility [12–16]. Chromobacterium has attracted considerable biotechnological interest for its potential as a biocontrol agent, environmental detoxification, and bioprospecting, in addition to industrial and pharmacological uses. Among its important characteristics, there are the production of chitinase (with fungicide, insecticide, and nematicide activities), polyhydroxyalkanoates (biodegradable plastics), violacein (with anticarcinogenic and antimicrobial activities), and cyanide biogenesis associated with gold recovery production [17]. These attributes were confirmed by the complete sequencing of the first strain of the type species C. violaceum, ATCC 12472, and by further genomic studies of C. vaccinii, C. subtsugae, C. haemolyticum, C. piscinae, C. aquaticum, and C. pseudoviolaceum, which showed alternative pathways for energy generation, a large number of open reading frames (ORFs) for transport-related proteins, and systems for stress adaptation, indicating the possible versatility and adaptability of this bacterium. Moreover, in spite of the fact that Chromobacterium is mainly considered a free-living microorganism, it is recognized as an important opportunistic pathogen and occasionally leads to lethal infections in mammals [18, 19].

Despite the potential importance of Chromobacterium in industry, agriculture, and medicine, few genomic analyses have been performed to gain insights into the biotechnological applications and pathogenicity of Chromobacterium species. Here, we present high-quality draft genome sequence of C. amazonense strain 56AF, isolated from a tropical freshwater lake. This organism was selected for sequencing due to its high resistance to β-lactam and as part of an ongoing effort to investigate the bacterial diversity on protected area (undisturbed environment). Additionally, a comparative genomic analysis of Chromobacterium species sequences available in public databases was performed to gain insight into the core and unique genes. To our knowledge, this is the first reported genome sequencing of a C. amazonense isolate.

2. Materials and Methods

2.1. Bacterial Isolate and DNA Extraction

Chromobacterium sp. strain 56AF was originally recovered in 2005 from the tropical freshwater Lake Dom Helvécio located in the Rio Doce State Park (Atlantic Rain biome), Minas Gerais, Brazil [12]. This has been a RAMSAR site since 2010 (http://www.ramsar.org) in recognition of its importance for the global conservation of biological diversity. Briefly, Chromobacterium sp. strain 56 AF was isolated on 25%-strength nutrient agar (Difco Laboratories, USA) at 25°C. More details were presented in a previous study of our group [12].

For the extraction of genomic DNA (gDNA), strain 56AF was grown in nutrient broth medium (Difco Laboratories, USA) at 25°C with shaking at 150 rpm for 24 h. The gDNA was extracted using the PureLink Genomic DNA kit (Thermo Fisher Scientific, USA) according to the manufacturer’s instructions. gDNA quantification was performed with a Qubit® fluorometer (Thermo Fisher Scientific).

2.2. Phylogenetic and Phylogenomic Analyses

After the description of new Chromobacterium species, the phylogenetic and taxonomic position of Chromobacterium sp. strain 56AF was determined using 16S rRNA gene sequences (1474 bp) from 21 Neisseriaceae and Chromobacteriaceae sequences retrieved from the GenBank database (http://www.ncbi.nlm.nih.gov/). Nucleic acids were aligned using MAFFT V7 [20] with iterative refinement using the G-INS-i strategy. To optimize the dataset for phylogenetic analysis, gap-rich columns were removed from the alignment using TrimAl (Gappyout option) (version 1.3) [21] available in Phylemon [22]. The best fit model for the multiple sequence alignment was estimated using ProtTest (version 3) [23]. Altogether, 12 different evolutionary models (Blosum62, CpREV, Dayhoff, DCMut, JTT, LG, MtArt, MtMam, MtREV, RtREV, VT, and WAG) were tested with the +I, +G, and +F parameters. The evolutionary model that best fit the data was determined according to the Akaike information criterion (AIC), and support values for each node were estimated using the approximate likelihood-ratio test (aLRT).

Digital DNA-DNA hybridization (dDDH) values between the genome sequences of Chromobacterium sp. strain 56AF and the available Chromobacterium type strains (GenBank assembly accessions GCA_000711885.1, GCA_000007705.1, GCA_000971335.1, GCA_001676875.1, GCA_001953795.1, GCA_001953775.1, and GCA_001855565.1) were estimated using the online analysis tool Genome-Genome-Distance Calculator 2.1 (GGDC) [24]. The recommended distance formula 2 was taken into account to interpret the results, as it is robust against the use of incomplete genome sequences.

2.3. Genome Sequencing and Assembly

The genome of strain 56AF was prepared using the Nextera DNA Library Preparation Kit (Illumina Inc., USA) and sequenced using a paired-end approach on the Illumina MiSeq platform. The average insert length was ~550 bp. Reads were trimmed and filtered with Trimmomatic software (version 0.32) [25], and sequences with a Phred Q score < 22, length < 35 bp, and ambiguous bases were discarded.

Trimmed reads were assembled de novo with Velvet software (version 1.2.10) [26]. Several assemblies were computed using k-mer values from 31 to 119. Assembly results were compared, and the best assembly (k-mer 95) was selected based on the N50 and the lowest number of contigs. This was further aligned with the Mauve Multiple Genome Alignment program against the complete genome of C. violaceum ATCC 12472. Possible misassemblies were examined with QUAST software (version 3.0) [27] (Table S1). The genome coverage was estimated by aligning the raw reads against the assembly with BWA [28] and Samtools [29], revealing a sequence coverage of 16,482x.

2.4. Accession Number

This whole-genome shotgun project has been deposited at DDBJ/ENA/GenBank under the accession number MTBD00000000. The version described in this paper is the first version, MTBD01000000.

2.5. Genome Annotation

ORFs were predicted using the RAST server databases. Additional functional annotations were performed with the SEED database and Kyoto Encyclopedia of Genes and Genomes (KEGG). Additionally, gene discovery was performed with Prodigal software [30]. The identified ORFs were checked against prot2003-2014.fa and cog2003-2014 from the National Center for Biotechnology Information nonredundant sequence database (NCBI-Nr) through a Blastp search for the assignment of Clusters of Orthologous Group (COG) functions. tRNA and rRNA genes were identified using the Aragorn and RNAmmer software programs, respectively. The coding density of the genome was determined with Artemis software (version 16.0.0) [31] based on GenBank files generated by the RAST server. A circular genome map was generated using the BLAST Ring Image Generator (BRIG) program [32]. Phage discovery was conducted using the PHAST tool.

2.6. Comparative Genomic Analysis

The C. amazonense strain 56AF genome sequence was compared to the genomes of C. violaceum ATCC 12472 [17], C. haemolyticum DSM 19808 [33], C. vaccinii 21-1 [34], C. piscinae ND17 [35], C. aquaticum CC-SEYA-1 [36], C. subtsugae MWU 2920 [37], C. pseudoviolaceum LMG3953 [38], and C. amazonense CBMAI 310 (unpublished, GenBank assembly accession: GCA_001855565.1). The predicted proteomes were used with OrthoFinder [39] to identify orthologous genes.

3. Results and Discussion

The phylogenetic position of Chromobacterium sp. strain 56AF, previously described by our group [12], was examined using the complete 16S rRNA gene sequences from 21 Neisseriaceae and Chromobacteriaceae sequences retrieved from the GenBank database. The resulting phylogeny shows the evolutionary relationships among these sequences and revealed that Chromobacterium sp. 56AF is closely related (94% of identity) to C. amazonense CBMAI 310, isolated from water samples from the Rio Negro, in the Amazon, Brazil (Figure 1). The pairwise alignment indicated that the percent identity between these two sequences is over 99.5%. In addition, the phylogenetic tree also demonstrated that the Chromobacterium genus forms a well-supported monophyletic clade.

Figure 1: Maximum likelihood phylogenetic tree showing the evolutionary position of C. amazonense strain 56AF in relation to other Chromobacterium species. Neisseriaceae sequences were included as the outgroup. The phylogeny was reconstructed based on 22 partial 16S rRNA gene sequences and 1470 sites. Support values for each node were estimated using the approximate likelihood-ratio test (aLRT). Boxes of different colors highlight distinct families.

The C. amazonense CBMAI 310 genome is 81.6% homologous with that of Chromobacterium sp. strain 56AF, whereas the other reference genomes exhibited dDDH values ranging from 23.10 to 42.9% when compared with the Chromobacterium sp. strain 56AF genome. Thus, based on the 16S rRNA gene phylogeny and dDDH data, Chromobacterium sp. strain 56AF has been renamed as C. amazonense strain 56AF.

3.1. General Genomic Features

The C. amazonense strain 56AF genome is 4,556,707 bp long and distributed across 141 contigs. A map of the chromosome and the features of the genome is shown in Figure 2 and Table 1, respectively. The N50 is 169,722 bp, and the GC content is 61.95%, consistent with the values reported for Chromobacterium species, which range from 62.23% to 64.83%. A total of 85.8% of the final assembly was annotated as coding regions. Out of 4294 predicted coding sequences (CDSs) with an average length of 911 bp, 4205 were protein-coding genes, and 95 were RNA genes, including 88 tRNAs, six rRNAs, and one tmRNA. Among the protein-coding genes, 3298 (76.8%) were matched to functions present in the database, whereas the remaining genes were annotated as hypothetical proteins. The classification of genes into functional categories generated by RAST is shown in Table 2. Prophage sequences identified by PHAST constitute 5.8% of the bacterial chromosome, including seven intact, two incomplete, and one questionable prophage (Table S2 and Figure S1).

Figure 2: (a) Circular map of the draft genome of C. amazonense strain 56AF. From outer to inner circle: genes on the forward and reverse strands (predicted coding sequences colored), tRNA genes, rRNA genes, GC content, and GC skew. The circular genome map was generated using BRIG. (b) BLAST comparison of the draft genome of C. amazonense strain 56AF against the genomes of Chromobacterium species. The innermost ring depicts C. violaceum ATCC12472, followed by the query sequences of C. amazonense strain 56AF, C. amazonense, C. piscinae, C. vaccinii, C. subtsugae, C. haemolyticum, C. pseudoviolaceum, and C. aquaticum.

Table 1: Comparison of the genomic properties of Chromobacterium spp.

Chromobacterium spp.
Size (Mb)
Contig
GC (%)
Gene
rRNA
tRNA
Other RNA
Reference

C. amazonense 56AF
4,56
141
62
4,294
8
89
1
This study

C. violaceum ATCC 12472
4,75
1
64.8
4,438
25
98
5
[32]

C. haemolyticum DSM 19808
5,03
67
62.8
4,595
13
57
4
[33]

C. piscinae ND 17
4,09
223
62.6
3,916
4
81
1
[35]

C. subtsugae MWU 2920
4,67
152
64.9
4,313

77
4
[37]

C. vaccinii MWU 205
4,97
152
64.44
4463
6
74
4
[34]

C. aquaticum CC-SEYA-1
4.77
1,63
64.8
4492

33
4
[36]

C. pseudoviolaceum LMG 3953
4.63
326
64.7
4404

33
5
[12]

Table 2: SEED subsystems distribution of the Chromobacterium amazonense strain 56AF genome based on MG-RAST annotation.

Subsystem category
Subsystem feature counts

Sulfur metabolism
32

Phosphorus metabolism
47

Carbohydrates
242

Amino acids and derivatives
417

Fatty acids, lipid, and isoprenoids
103

Metabolism secondary
4

Nitrogen metabolism
13

Metabolism of aromatic compounds
24

Cofactors, vitamins, prosthetic groups, and pigments
242

Nucleosides and nucleotides
90

Dormancy and sporulation
1

Respiration
133

Stress response
134

DNA metabolism
99

Membrane transport
156

Regulation and cell signaling
88

Motility and chemotaxis
176

Cell division and cell cycle
35

Protein metabolism
298

RNA metabolism
135

Iron acquisition and metabolism
34

Phages, prophages, transposable elements and plasmids
65

Miscellaneous
45

Photosynthesis
0

Potassium metabolism
20

Virulence, disease, and defense
88

Cell wall and capsule
166

Among the predicted genes present in the genome of C. amazonense strain 56AF, we focused on virulence factors, toxins, and secretion systems, which might be of particular importance in terms of pathogenicity, survival, and adaptation to the environment.

3.2. Insights from the Genome Sequence

3.2.1. Virulence Factors

A number of potential Chromobacterium virulence factors have been previously identified, including different types of secretion systems and toxins that are likely associated with C. violaceum and C. haemolyticum fatal infections in humans [17, 34]. Consistent with these reports, the draft genome of C. amazonense strain 56AF includes CDSs with significant similarity to virulence factors (Table S3) that were previously described for strains belonging to the Chromobacterium genus among others, as discussed in the following sections.

3.2.2. Secretion Systems

There are five well-studied double-membrane-spanning secretion systems (T1SS, T2SS, T3SS, T4SS, and T6SS) present in a wide variety of pathogenic and nonpathogenic gram-negative bacteria. These systems secrete products with diverse functions, such as those involved in adhesion, pathogenicity, adaptation, and survival [40]. Here, we identified a number of particularly interesting proteins that are part of the T1SS, T2SS, T3SS, T4SS, and T6SS machinery.

In gram-negative bacteria, T1SS is composed of three membrane proteins, two in the inner membrane belonging to the ABC transporter and membrane fusion protein (MFP) families and a third porin-like protein in the outer membrane [41], known as LapB, LapC, and TolC, respectively. The products secreted by T1SS are diverse in size and function and frequently associated with nutrient acquisition (lipases and proteases) and virulence, for example, colicin V, which was found in the C. amazonense strain 56AF genome.

T2SS consists of a large multiprotein machinery with 12–15 components known as general secretion pathway (Gsp) proteins [40]. Analysis of the C. amazonense strain 56AF genome revealed the presence of all of these proteins with the exception of GspC, which may be missing from the genome assembly. An additional T2SS protein, GspN, was identified; however, its function and location are unclear [42].

T3SS is found in a large number of gram-negative pathogens and symbionts [43]. This secretion system is composed of structural, translocator, and effector proteins. Similar to C. violaceum ATCC 12472, C. amazonense strain 56AF possesses the Chromobacterium pathogenicity islands (Cpi) 1, 1a, and 2, whereas C. haemolyticum T124 lacks the Cpi-2-like second T3SS. These data suggest the potential virulence of C. amazonense strain 56AF.

T4SS is found in both gram-negative and gram-positive bacteria and in some archaea and allows the translocation of DNA and proteins into target cells, playing an important role in the pathogenesis of a wide range of bacteria [44]. Among the secretion systems, T4SS is unique for its ability to mediate the conjugation of plasmid DNA [40]. It consists of 12 proteins (VirB1–VirB11 and VirD4), of which only VirB3 and VirB7 were missing from the C. amazonense strain 56AF genome assembly. This is the first report of T4SS in the Chromobacterium genus, which is considered to be critical for pathogenicity.

Finally, T6SS is broadly distributed among Proteobacteria. It translocates toxic effector proteins into eukaryotic and prokaryotic cells and has an important role in pathogenesis, bacterial communication, and interactions with the environment [45–47]. T6SS consists of a 13-protein core and a set of conserved accessory proteins [48]. In the genome of C. amazonense strain 56AF, 7 out of the 13 T6SS-associated proteins (ImpA, ImpB, ImpC, ImpF, ImpG, ImpH, and ImpJ) were found, in addition to the T4SS homolog IcmF. It is also known that T6SS secretes hemolysin and the effector protein VrgG, both of which are present in the C. amazonense strain 56AF genome.

3.2.3. Toxins

Toxins can be identified by the presence of a variety of proteins, such as hemolysins, hemolytic enterotoxins, colicin V, lytic proteins, and Nudix hydrolases. These toxins are associated with pathogenic bacteria [49–52] and are present in the genome of C. amazonense strain 56AF.

Important virulence factors for host cell invasion were found, for example, CDSs encoding a hydrolase cell-wall associated protein, metalloproteases, lypolitic proteins, collagenases, and the proteolytic enzymes KpSS and KpSC. The last two are conserved proteins involved in capsular polysaccharides assembly systems [53], suggesting the potential of C. amazonense strain 56AF to synthesize capsular polysaccharides, which can act as important virulence factors in many pathogenic bacteria [53]. Moreover, genes involved in flagellar biosynthesis were identified. Flagella and swimming motility are present in all the genera of the family Chromobacteriaceae [1] and are considered to be important contributors to host colonization for most pathogens.

Unlike other Chromobacterium species, the genome of C. amazonense strain 56AF contains the gene for pore-forming toxin alpha-hemolysin (hlyA), which shares homology with the gene in Vibrio cholerae. It should be noted that hemolysin is a well-known virulence factor, and hemolytic activity has been detected in other Chromobacterium species, such as C. violaceum, C. haemolyticum, C. aquaticum, and C. subtsugae [5, 17, 54]. Thus, some Chromobacterium species have the potential to produce different cytolytic toxins that may contribute to their pathogenicity.

A set of genes with potential roles in multiple antibiotic resistance were identified in the C. amazonense strain 56AF genome and included CDSs encoding the RND efflux transporter, tripartite multidrug resistance system, Bcr/CflA family drug resistance transporter, multiple antibiotic resistance MarC protein, and resistance to β-lactams. Some efflux pumps not only confer clinically relevant resistance to antibiotics but also have a role in bacterial pathogenicity and may be beneficial for bacterial survival [55]. Regarding β-lactams, C. amazonense strain 56AF genome presented CDSs related to β-lactamases, including extended spectrum β-lactamases Ambler class C β-lactamase, penicillinases, metallo-β-lactamase, and metal-dependent hydrolases of the β -lactamase superfamily I PhnP protein. It is worth mentioning that they hydrolyze almost all clinically used β-lactams including carbapenems, thereby representing a therapeutic challenge. These results were supported by antimicrobial susceptibility testing obtained by Lima-Bittencourt et al. [12] for Chromobacterium sp. strain 56AF, which found resistance to ampicillin, ampicillin/clavulanic acid, and cefotaxime (third-generation cephalosporin), widely used for the treatment of bacterial infections. Other Chromobacterium species also exhibit β-lactamases, such as C. piscinae and C. haemolyticum [56]. Other antibiotic resistance genes encode resistance to bicyclomycin, fosfomycin, mitomycin, polymyxin (PmrJ, ArnC, and ArnT proteins), and fosmidomycin. Importantly, the genome contains genes related to resistance to several toxic compounds, including tellurite, cobalt, zinc-cadmium (CzcA), and fusaric acid, a mycotoxin produced by several Fusarium species that cause diseases in economically important crops [57].

Additionally, the annotation of the C. amazonense strain 56AF revealed its potential to produce esterases, lipase, auxins, chitinases, phytoene synthase and phytoene desaturase, polyhydroxyalkanoates, glutathione S-transferases, and violacein, suggesting that this strain might have potential for biotechnological purposes. The genome contains predicted toxin-antitoxin genes, in particular those of higAB and doc/phd, involved in plasmid stabilization [58], as well as the antitoxin YgiT (MqsA), which helps to mediate the bacterial general stress response [59], although the gene encoding MqsR toxin was absent. Finally, plastocyanin/azurin genes were also found, suggesting that this strain may participate in denitrifying processes.

3.3. Comparison of Eight Chromobacterium Genomes

Based on orthologous groups (OGs) analysis, the eight Chromobacterium genomes examined contain a total of 7220 OGs, with a conserved core genome of 2467 OGs present in all species (Figure 3). A large fraction of the flexible genes (65.8%), that is, those present in one or more but not all of the genomes, was mapped, among which 2069 were species-specific and 2684 were shared noncore. These findings are consistent with the great metabolic flexibility of Chromobacterium, conferring to their members a possible selective advantage to colonize and to survive in diverse ecosystems. It should be noted that C. violaceum contained the largest number of OGs shared with other Chromobacterium genomes. Comparisons of C. amazonense strain 56AF with the seven other genomes revealed a subset of OGs (301) that were unique to C. amazonense strain 56AF, and these may contribute to the species-specific features of C. amazonense. The unique set of C. amazonense strain 56AF genes was associated with the COG functional categories carbohydrate transport and metabolism (G); cell cycle control, cell division, and chromosome partitioning (D); amino acid transport and metabolism (E); nucleotide transport and metabolism (F); replication, recombination, and repair (L); defense mechanisms (V); and cell motility (N). Considering that seven out of the eight genomes analyzed here are draft sequences (all except for C. violaceum), it is likely that individual genes are missing from these assemblies.

Figure 3: Distribution of orthologous groups (OGs) among the genomes of C. amazonense, C. piscinae, C. pseudoviolaceum, C. violaceum, C. amazonense strain 56AF, C. haemolyticum, C. aquaticum, C. subtsugae, and C. vaccinii. (a) The bars represent the number of unique or shared OGs among Chromobacterium genomes marked with the filled circle. (b) The heat map shows the number of shared OGs between different Chromobacterium species. The dendrogram is based on Euclidean distance of the shared OG matrix.

4. Conclusions

C. amazonense strain 56AF is the first strain of the genus known to contain predicted genes associated with five of the six double-membrane-spanning secretion systems (types 1–4 and 6) present in gram-negative bacteria, of which T4SS has not been reported in any other sequenced Chromobacterium. Additionally, the presence of several predicted genes encoding potential multidrug resistance and detoxification compounds suggests that this isolate harbors a genomic repertoire capable of dealing with environmental toxins, which may be beneficial to its adaptation and survival. Moreover, C. amazonense strain 56AF contains a number of relevant virulence-related genes, suggesting its potential pathogenicity. In summary, our findings confirm the versatility, adaptability, and biotechnological potential of this bacterium. Moreover, our results provide important insights into the genomics of Chromobacterium species, leading to a better understanding of the biology of this genus. Further functional studies will contribute with new information about its pathogenic potential.

Abbreviations

ORFs: | Open reading frames
gDNA: | Genomic DNA
dDDH: | Digital DNA-DNA hybridization
KEGG: | Kyoto Encyclopedia of Genes and Genomes
COGs: | Clusters of Orthologous Groups
OGs: | Orthologous groups
TSS: | Type secretion system.

Conflicts of Interest

The authors declare no conflict of interests.

Acknowledgments

This work and publication were supported by the grants from the Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG, APQ 00801/12 and APQ-00463-13), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).

Supplementary Materials

Table S1: misassembly report based on the alignment of C. amazonense strain 56AF, extracted from QUAST. Table S2: phages present in the genome of C. amazonense strain 56AF. Table S3: virulence factors identified in the genome of C. amazonense strain 56AF according to the virulence factor database of pathogenic bacteria. Figure S1: phages present in the genome of C. amazonense strain 56AF. Circular view extracted from PHAST. (Supplementary Materials)

References

M. Adeolu and R. S. Gupta, “Phylogenomics and molecular signatures for the order Neisseriales: proposal for division of the order Neisseriales into the emended family Neisseriaceae and Chromobacteriaceae fam. nov.,” Antonie Van Leeuwenhoek, vol. 104, no. 1, pp. 1–24, 2013.

C. Bergonzini, “Um nuevo bacterio colorato,” Annuario della Societa dei Naturalisti in Moderna, vol. 2, pp. 149–158, 1881.

P. A. W. Martin, D. Gundersen-Rindal, M. Blackburn, and J. Buyer, “Chromobacterium subtsugae sp. nov., a betaproteobacterium toxic to Colorado potato beetle and other insect pests,” International Journal of Systematic and Evolutionary Microbiology, vol. 57, no. 5, pp. 993–999, 2007.

C. C. Young, A. B. Arun, W. A. Lai et al., “Chromobacterium aquaticum sp. nov., isolated from spring water samples,” International Journal of Systematic and Evolutionary Microbiology, vol. 58, no. 4, pp. 877–880, 2008.

X. Y. Han, F. S. Han, and J. Segal, “Chromobacterium haemolyticum sp. nov., a strongly haemolytic species,” International Journal of Systematic and Evolutionary Microbiology, vol. 58, no. 6, pp. 1398–1403, 2008.

P. Kämpfer, H. J. Busse, and H. C. Scholz, “Chromobacterium piscinae sp. nov. and Chromobacterium pseudoviolaceum sp. nov., from environmental samples,” International Journal of Systematic and Evolutionary Microbiology, vol. 59, no. 10, pp. 2486–2490, 2009.

S. D. Soby, S. R. Gadagkar, C. Contreras, and F. L. Caruso, “Chromobacterium vaccinii sp. nov., isolated from native and cultivated cranberry (Vaccinium macrocarpon Ait.) bogs and irrigation ponds,” International Journal of Systematic and Evolutionary Microbiology, vol. 63, Part 5, pp. 1840–1846, 2013.

C. B. A. Menezes, M. F. Tonin, D. B. A. Corrêa et al., “Chromobacterium amazonense sp. nov. isolated from water samples from the Rio Negro, Amazon, Brazil,” Antonie Van Leeuwenhoek, vol. 107, no. 4, pp. 1057–1063, 2015.

S. Zhou, X. Guo, H. Wang et al., “Chromobacterium rhizoryzae sp. nov., isolated from rice roots,” International Journal of Systematic and Evolutionary Microbiology, vol. 66, no. 10, pp. 3890–3896, 2016.

A. Bajaj, A. Kumar, S. Yadav et al., “Isolation and characterization of a novel Gram-negative bacterium Chromobacterium alkanivorans sp. nov., strain IITR-71T degrading halogenated alkanes,” International Journal of Systematic and Evolutionary Microbiology, vol. 66, no. 12, pp. 5228–5235, 2016.

M. B. Blackburn, R. R. Farrar Jr, M. E. Sparks, D. Kuhar, A. Mitchell, and D. E. Gundersen-Rindal, “Chromobacterium sphagni sp. nov., an insecticidal bacterium isolated from Sphagnum bogs,” International Journal of Systematic and Evolutionary Microbiology, vol. 67, no. 9, pp. 3417–3422, 2017.

C. I. Lima-Bittencourt, S. Astolfi-Filho, E. Chartone-Souza, F. R. Santos, and A. M. A. Nascimento, “Analysis of Chromobacterium sp. natural isolates from different Brazilian ecosystems,” BMC Microbiology, vol. 7, no. 1, p. 58, 2007.

C. I. Lima-Bittencourt, P. S. Costa, F. A. R. Barbosa, E. Chartone-Souza, and A. M. A. Nascimento, “Characterization of a Chromobacterium haemolyticum population from a natural tropical lake,” Letters in Applied Microbiology, vol. 52, no. 6, pp. 642–650, 2011.

E. S. Barreto, A. R. Torres, M. R. Barreto, A. T. R. Vasconcelos, S. Astolfi-Filho, and M. Hungria, “Diversity in antifungal activity of strains of Chromobacterium violaceum from the Brazilian Amazon,” Journal of Industrial Microbiology and Biotechnology, vol. 35, no. 7, pp. 783–790, 2008.

L. T. Dall'Agnol, R. N. Martins, A. C. R. Vallinoto, and K. T. S. Ribeiro, “Diversity of Chromobacterium violaceum isolates from aquatic environments of state of Pará, Brazilian Amazon,” Memórias do Instituto Oswaldo Cruz, vol. 103, no. 7, pp. 678–682, 2008.

K. Ponnusamy, S. Jose, I. Savarimuthu, G. P. Michael, and M. Redenbach, “Genetic diversity study of Chromobacterium violaceum isolated from Kolli Hills by amplified ribosomal DNA restriction analysis (ARDRA) and random amplified polymorphic DNA (RAPD),” Letters in Applied Microbiology, vol. 53, no. 3, pp. 341–349, 2011.

Brazilian National Genome Project Consortium, A. T. R. de Vasconcelos, D. F. de Almeida et al., “The complete genome sequence of Chromobacterium violaceum reveals remarkable and exploitable bacterial adaptability,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 20, pp. 11660–11665, 2011.

A. Chattopadhyay, V. Kumar, N. Bhat, and P. L. N. G. Rao, “Chromobacterium violaceum infection: a rare but frequently fatal disease,” Journal of Pediatric Surgery, vol. 37, no. 1, pp. 108–110, 2002.

M. Okada, R. Inokuchi, K. Shinohara et al., “Chromobacterium haemolyticum-induced bacteremia in a healthy young man,” BMC Infectious Diseases, vol. 13, no. 1, p. 406, 2013.

K. Katoh and D. M. Standley, “MAFFT multiple sequence alignment software version 7: improvements in performance and usability,” Molecular Biology and Evolution, vol. 30, no. 4, pp. 772–780, 2013.

S. Capella-Gutiérrez, J. M. Silla-Martínez, and T. Gabaldón, “trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses,” Bioinformatics, vol. 25, no. 15, pp. 1972-1973, 2009.

R. Sánchez, F. Serra, J. Tárraga et al., “Phylemon 2.0: a suite of web-tools for molecular evolution, phylogenetics, phylogenomics and hypotheses testing,” Nucleic Acids Research, vol. 39, pp. W470–W474, 2011.

D. Darriba, G. L. Taboada, R. Doallo, and D. Posada, “ProtTest 3: fast selection of best-fit models of protein evolution,” Bioinformatics, vol. 27, no. 8, pp. 1164-1165, 2011.

J. P. Meier-Kolthoff, A. F. Auch, H. P. Klenk, and M. Göker, “Genome sequence-based species delimitation with confidence intervals and improved distance functions,” BMC Bioinformatics, vol. 14, no. 1, p. 60, 2013.

A. M. Bolger, M. Lohse, and B. Usadel, “Trimmomatic: a flexible trimmer for Illumina sequence data,” Bioinformatics, vol. 30, no. 15, pp. 2114–2120, 2014.

D. R. Zerbino and E. Birney, “Velvet: algorithms for de novo short read assembly using de Bruijn graphs,” Genome Research, vol. 18, no. 5, pp. 821–829, 2008.

A. Gurevich, V. Saveliev, N. Vyahhi, and G. Tesler, “QUAST: quality assessment tool for genome assemblies,” Bioinformatics, vol. 29, no. 8, pp. 1072–1075, 2013.

H. Li and R. Durbin, “Fast and accurate long-read alignment with Burrows-Wheeler transform,” Bioinformatics, vol. 26, no. 5, pp. 589–595, 2010.

H. Li, B. Handsaker, A. Wysoker et al., “The sequence alignment/map format and SAMtools,” Bioinformatics, vol. 25, no. 16, pp. 2078-2079, 2009.

D. Hyatt, G.-L. Chen, P. F. LoCascio, M. L. Land, F. W. Larimer, and L. J. Hauser, “Prodigal: prokaryotic gene recognition and translation initiation site identification,” BMC Bioinformatics, vol. 11, no. 1, pp. 119–130, 2010.

K. Rutherford, J. Parkhill, J. Crook et al., “Artemis: sequence visualization and annotation,” Bioinformatics, vol. 16, no. 10, pp. 944-945, 2000.

N.-F. Alikhan, N. K. Petty, N. L. Ben Zakour, and S. A. Beatson, “BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons,” BMC Genomics, vol. 12, no. 1, p. 402, 2011.

T. Miki and N. Okada, “Draft genome sequence of Chromobacterium haemolyticum causing human bacteremia infection in Japan,” Genome Announcements, vol. 2, no. 6, pp. e01047–e01014, 2014.

K. Vöing, A. Harrison, and S. D. Soby, “Draft genome sequence of Chromobacterium vaccinii, a potential biocontrol agent against mosquito (Aedes aegypti) larvae,” Genome Announcements, vol. 3, no. 3, pp. e00477–e00415, 2015.

K.-G. Chan and N. Y. M. Yunos, “Whole-genome sequencing analysis of Chromobacterium piscinae strain ND17, a quorum-sensing bacterium,” Genome Announcements, vol. 4, no. 2, pp. e00081–e00016, 2016.

S. D. Soby, “Draft genome sequence of Chromobacterium aquaticum CCSEYA-1, a nonpigmented member of thegenus Chromobacterium,” Genome Announcements, vol. 5, no. 12, pp. e01661–e01616, 2017.

K. Vöing, A. Harrison, and S. D. Soby, “Draft genome sequences of three Chromobacterium subtsugae isolates from wild and cultivated cranberry bogs in southeastern Massachusetts,” Genome Announcements, vol. 3, no. 5, pp. e00998–e00915, 2015.

S. D. Soby, “Draft genome sequence of Chromobacterium pseudoviolaceum LMG 3953T, an enigmatic member of the genus Chromobacterium,” Genome Announcements, vol. 5, no. 12, pp. e01632–e01616, 2017.

D. M. Emms and S. Kelly, “OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy,” Genome Biology, vol. 16, no. 1, p. 157, 2015.

T. R. D. Costa, C. Felisberto-Rodrigues, A. Meir et al., “Secretion systems in Gram-negative bacteria: structural and mechanistic insights,” Nature Reviews Microbiology, vol. 13, no. 6, pp. 343–359, 2015.

S. Thomas, I. B. Holland, and L. Schmitt, “The type 1 secretion pathway –the hemolysin system and beyond,” Biochimica et Biophysica Acta (BBA) - Molecular Cell Research, vol. 1843, no. 8, pp. 1629–1641, 2014.

K. V. Korotkov, M. Sandkvist, and W. G. J. Hol, “The type II secretion system: biogenesis, molecular architecture and mechanism,” Nature Reviews Microbiology, vol. 10, no. 5, pp. 336–351, 2012.

B. Coburn, I. Sekirov, and B. B. Finlay, “Type III secretion systems and disease,” Clinical Microbiology Reviews, vol. 20, no. 4, pp. 535–549, 2007.

C. E. Alvarez-Martinez and P. J. Christie, “Biological diversity of prokaryotic type IV secretion systems,” Microbiology and Molecular Biology Reviews, vol. 73, no. 4, pp. 775–808, 2009.

G. Bönemann, A. Pietrosiuk, A. Diemand, H. Zentgraf, and A. Mogk, “Remodelling of VipA/VipB tubules by ClpV-mediated threading is crucial for type VI protein secretion,” EMBO Journal, vol. 28, no. 4, pp. 315–325, 2009.

B. T. Ho, T. G. Dong, and J. J. Mekalanos, “A view to a kill: the bacterial type VI secretion system,” Cell Host & Microbe, vol. 15, no. 1, pp. 9–21, 2014.

A. Zoued, Y. R. Brunet, E. Durand et al., “Architecture and assembly of the type VI secretion system,” Biochimica et Biophysica Acta (BBA) - Molecular Cell Research, vol. 1843, no. 8, pp. 1664–1673, 2014.

F. Boyer, G. Fichant, J. Berthod, Y. Vandenbrouck, and I. Attree, “Dissecting the bacterial type VI secretion system by a genome wide in silico analysis: what can be learned from available microbial genomic resources?” BMC Genomics, vol. 10, no. 1, p. 104, 2009.

J. A. Vázquez-Boland, M. Kuhn, P. Berche et al., “Listeria pathogenesis and molecular virulence determinants,” Clinical Microbiology Reviews, vol. 14, no. 3, pp. 584–640, 2001.

P. H. Edelstein, B. Hu, T. Shinzato, M. A. C. Edelstein, W. Xu, and M. J. Bessman, “Legionella pneumophila NudA is a Nudix hydrolase and virulence factor,” Infection and Immunity, vol. 73, no. 10, pp. 6567–6576, 2005.

E. Ortega, H. Abriouel, R. Lucas, and A. Gálvez, “Multiple roles of Staphylococcus aureus enterotoxins: pathogenicity, superantigenic activity, and correlation to antibiotic resistance,” Toxins, vol. 2, no. 8, pp. 2117–2131, 2010.

N. Jeßberger, R. Dietrich, S. Bock, A. Didier, and E. Märtlbauer, “Bacillus cereus enterotoxins act as major virulence factors and exhibit distinct cytotoxicity to different human cell lines,” Toxicon, vol. 77, pp. 49–57, 2014.

L. M. Willis and C. Whitfield, “KpsC and KpsS are retaining 3-deoxy-d-manno-oct-2-ulosonic acid (Kdo) transferases involved in synthesis of bacterial capsules,” Proceedings of the National Academy of Sciences of the United States of America, vol. 110, no. 51, pp. 20753–20758, 2013.

P. D. Rekha, C. C. Young, and A. B. Arun, “Identification of N-acyl-L-homoserine lactones produced by non-pigmented Chromobacterium aquaticum CC-SEYA-1T and pigmented Chromobacterium subtsugae PRAA4-1T,” Biotech, vol. 1, no. 4, pp. 239–245, 2011.

L. J. V. Piddock, “Multidrug-resistance efflux pumps - not just for resistance,” Nature Reviews Microbiology, vol. 4, no. 8, pp. 629–636, 2006.

D. D. Gudeta, V. Bortolaia, A. Jayol, L. Poirel, P. Nordmann, and L. Guardabassi, “Chromobacterium spp. harbour ambler class A β-lactamases showing high identity with KPC,” Journal of Antimicrobial Chemotherapy, vol. 71, no. 6, pp. 1493–1496, 2016.

C. W. Bacon, J. K. Porter, W. P. Norred, and J. F. Leslie, “Production of fusaric acid by Fusarium species,” Applied and Environmental Microbiology, vol. 62, no. 11, pp. 4039–4043, 1996.

K. Gerdes, S. K. Christensen, and A. Løbner-Olesen, “Prokaryotic toxin–antitoxin stress response loci,” Nature Reviews Microbiology, vol. 3, no. 5, pp. 371–382, 2005.

X. Wang, Y. Kim, S. H. Hong et al., “Antitoxin MqsA helps mediate the bacterial general stress response,” Nature Chemical Biology, vol. 7, no. 6, pp. 359–366, 2011.