Comparative genomic analysis between newly sequenced Brucella suis Vaccine Strain S2 and the Virulent Brucella suis Strain 1330
- Dong-dong Di†1,
- Hai Jiang†3,
- Li-li Tian1,
- Jing-li Kang1,
- Wen Zhang3,
- Xin-ping Yi4,
- Feng Ye4,
- Qi Zhong4,
- Bo Ni1,
- You-yu He5,
- Lin Xia5,
- Yao Yu5,
- Bu-yun Cui3Email author,
- Xiang Mao2Email author and
- Wei-xing Fan1Email author
© The Author(s). 2016
Received: 10 April 2015
Accepted: 7 September 2016
Published: 20 September 2016
Brucellosis is a bacterial disease caused by Brucella infection. In the late fifties, Brucella suis vaccine strain S2 with reduced virulence was obtained by serial transfer of a virulent B. suis biovar 1 strain in China. It has been widely used for vaccination in China since 1971. Until now, the mechanisms underlie virulence attenuation of S2 are still unknown.
In this paper, the whole genome sequencing of S2 was carried out by Illumina Hiseq2000 sequencing method. We further performed the comparative genomic analysis to find out the differences between S2 and the virulent Brucella suis strain 1330. We found premature stops in outer membrane autotransporter omaA and eryD genes. Single mutations were found in phosphatidylcholine synthase, phosphorglucosamine mutase, pyruvate kinase and FliF, which have been reported to be related to the virulence of Brucella or other bacteria. Of the other different proteins between S2 and 1330, such as Omp2b, periplasmic sugar-binding protein, and oligopeptide ABC transporter, no definitive implications related to bacterial virulence were found, which await further investigation.
The data presented here provided the rational basis for designing Brucella vaccines that could be used in other strains.
KeywordsBrucella suis vaccine strain S2 Brucella suis virulent strain 1330 Genome sequence Comparative genomic analysis
Brucellosis is a worldwide zoonotic disease caused by the infection of Brucella species, a Gram-negative and facultative intracellular bacterium. Brucellosis can be acquired by direct contact with infected animals including cattle, sheep, goats and swine. Brucella infects approximately 500,000 humans worldwide annually according to the World Health Organization (WHO) . The infection can cause reproductive failure in food animals and even the loss of human productivity. The clinical symptoms of Brucella infection include undulant fever, abortion, asthenia, endocarditis, and encephalitis. Brucellosis in food animals could be controlled by vaccination, while human brucellosis can be treated with antibiotics.
Currently, there are ten recognized Brucella species based on their host preference and phenotypic differences: B. melitensis (goats), B. abortus (cattle), B. suis (swine), B. canis (dogs), B. ovis (sheep), B. neotomae (desert mice), B. ceti (cetaceans), B. pinnipedialis (pinnipeds), B. microti (common voles) and B. inopinata (unknown) . The genome of B. melitensis strain 16M is the first one to be sequenced , followed by B. suis strain 1330 [4, 5] and B. abortus strain 2308 and 9–941 . New Brucella genome sequences are increasing, and at least 39 Brucella strains have been sequenced. These genomic resources can be served for identifying Brucella gene variability that associated with their host preference, pathogenesis, and virulence.
B. suis strain 1330 (hereafter referred to as 1330) is a swine isolate and used as the standard reference strain for B. suis biovar 1 . According to phylogenetic analysis of 10 Brucella strains, B. suis has a broader host specificity without any identified species-specific markers .
The B. suis strain S2 is isolated from swine fetus by the Institute of Inspection for Veterinary Medicine in 1952 in China. It possesses all of the characteristics of B. suis biovar 1 and a smooth colonial morphology . After serial transfer on media for years, B. suis strain S2 is naturally attenuated to a vaccine strain designated as B. suis vaccine “strain 2” or S2 (hereafter referred to as S2). The avirulent property of S2 is stable when inoculated into susceptible animals, and it does not cause the abortion of pregnant female . Therefore, it has been widely used in China since 1971 to vaccinate sheep and goats . The vaccine is administered through drinking water, a dose of 10 billion bacteria provides the protection for 2–3 years. In the middle of 1980s, S2 is introduced to the other countries . Although studies suggest that the protection rate of S2 is less than that of B. abortus “strain 19” or S19 (hereafter referred to as S19) and B. melitensis Rev. 1 (hereafter referred to as Rev. 1) [11, 12], a large-scale animal test shows that S2 can provide a satisfactory protection rate in the field [9, 13]. Based on these encouraging results, the low virulence compared with S19 and Rev.1, the low cost of production, the safety for pregnant sheep and the feasibility of administration by the oral route in cattle, pigs, sheep and goats, S2 is used successfully in China to vaccinate all target species against brucellosis by administration through drinking water [9–12]. In mice, S2 has low residual virulence and less persistence time compared with S19 and Rev.1. It provides short-time protection against virulence B. melitensis challenge compared with the other two vaccines . However, it does not show significant protection against B. melitensis 53H38 challenge when tested in lambs .
The underlying molecular and physiological mechanisms of the low virulence of S2 are not well understood so far. In this study, we performed the paired-end sequencing combined with Sanger sequencing method to determine the complete genome sequences of S2. We further identified the different genes between S2 and 1330 by genome comparison. Several putative factors were identified as the candidate genes that might be related to the virulence attenuation of S2. Further investigation of these target genes will be important for designing better vaccines to provide efficient control for brucellosis.
Results and discussion
Genome sequence properties
Genome properties of the newly sequenced genome of B. suis strain S2 in comparison with the known genome sequence of B. suis strain 1330
B. suis strain 2
B. suis 1330
Average gene length(bp)
Single nucleotide polymorphism comparison between B. suis strain S2 and 1330
The identified non-identical proteins between attenuated strain S2 and virulent strain 1330
S2 vs. 1330
hypothetical protein (found)
IclR family transcriptional regulator (paper)
glycosyl hydrolase family protein(found)
outer membrane autotransporter (paper)
H → R
V → A
S → Y
sugar ABC transporter, periplasmic sugar-binding protein (paper)
Q → R
D → A
R → H
Y → H
hypothetical protein (found)
K → T
hypothetical protein (found)
V → A
R → Q
penicillin-binding protein 1A (unknown)
H → R
acetyl-CoA hydrolase/transferase family protein (surface)
P → S
serine protease (unknown)
S → G
R → H
dedA family protein (paper)
H → R
A → V
V → A
hypothetical protein (found)
L → R
glmM (no direct interaction with substrate)
T → A
pyk (no direct interaction with substrate)
Y → H
E → K
GMC family oxidoreductase (surface)
A → V
sensor histidine kinase BvrS (unknown)
A → D
PhoH family protein (found)
IS711, transposase orfA
luciferase family protein
short chain dehydrogenase
pyridine nucleotide-disulfide oxidoreductase family protein (found)
3-hydroxyisobutyrate dehydrogenase (paper)
sugar ABC transporter ATPase (found)
IclR family transcriptional regulator (paper)
A → T
A → E
T → A
major facilitator family transporter (found)
P → Q
oligopeptide ABC transporter, periplasmic oligopeptide-binding protein (paper)
S → G
oligopeptide ABC transporter, periplasmic oligopeptide-binding protein (paper)
S → L
A → P
peptide ABC transporter, periplasmic peptide-binding protein (found)
V → A
hypothetical protein (found)
M → T
E → K
glucuronate isomerase (interface)
N → D
aminotransferase, class III (surface)
W → R
oxidoreductase, molybdopterin-binding (unknown)
S → P
WD repeat-containing protein (unknown)
C → W
hypothetical protein (found)
P → S
M → L
branched-chain amino acid ABC transporter, ATP-binding protein
outer surface protein
branched-chain amino acid ABC transporter, periplasmic amino acid-binding protein
Putative protein-encoding genes, tRNA and rRNA in S2 genome were predicted with RAST (http://rast.nmpdr.org). The predicted protein sequences were compared with non-redundant protein database. The protein sequences were also compared with COG and CDD database, the threshold is E < 1e-2.
Identification of virulence associated differences between the attenuated strain S2 and virulent strain 1330
The pairwise comparison was carried out to identify the non-identical ORFs between S2 and 1330 (Table 3), which provided useful information on the underlying mechanisms for the attenuation of S2.
Most of the different ORFs between 1330 and S2 in chromosome I and II had only 1 bp difference in nucleotide sequence. Six ORFs had 2, 6, 8, 16, 69 and 78 bp insertion or deletion in the genes respectively (Table 3). The insertion or deletion in the gene caused either frameshifting (including premature stop) or residue insertion in the coding protein. The different ORFs in S2 that encoded the same protein sequences in other Brucella strains were labeled ‘found’ and omitted in discussions. The different ORFs labeled with ‘surface’ suggested that the different residue was located on the surface of the modeled structure, which were also omitted from the analysis (Table 3).
Known virulence associated differences
In chromosome I of S2, a 2 bp deletion caused the premature stop of an outer membrane autotransporter OmaA. The autotransporters are secreted proteins that can translocate themselves through the membranes to the cell surface or to the extracellular environment. Bandara et al. studied the function of this OmaA . In mice, the OmaA-deficient strain provides greater protection against the challenge of 1330 than that of B. suis mutant strain VTRS1. The deficient strain can induce a significantly higher level of serum IgG1 and IgG2a antibodies and can be cleared quickly from the mice. Although the OmaA activity is not necessary in the course of the acute phase of infection of B. suis, it might still be a critical virulent factor during the chronic phase of the infection. Therefore, we suspected the lack of intact OmaA protein might be the major cause for S2 virulence attenuation.
The second difference between 1330 and S2 was found in the erythritol (ery) operon. The erythritol operon contains 4 ORFs (eryA, eryB, eryC and eryD) . Erythritol metabolism by Brucella has been known to be associated with the abortion in livestock. The high concentration of erythritol in the foetal tissue of the animal (cattle, sheep, goat and pig) might be advantageous for the growth of Brucella. On the other hand, erythritol catabolism pathway might play the important physiological roles in bacteria growth. It has been suggested that erythrose 4-phosphate, the precursor for the biosynthesis of aromatic compounds, may be easily obtained from the intermediates of erythritol catabolism. Compared to 1330, a 16 bp insertion caused the premature stop of the eryD gene in S2. The studies on B. abortus S19 and S2308 showed that eryD might act as a repressor of erythritol operon . The consequence of eryC and eryD ORFs deletion in S19 has been studied, which reveals that they are not sufficient or required for virulence in a murine model. No direct correlation between erythritol metabolism and in vivo colonization was found when S19 and 2308 strains were compared with three genetically engineered strains related to ery operon . Considering its repressive control on ery operon, the lack of eryD might cause the uncontrolled erythritol catabolism. Nevertheless, its consequence on virulence attenuation needs to be further explored.
Effect of erythritol on growth of Brucella suis strain S2 and B. abortus strain S19
Phosphatidylcholine synthase gene (PCS) in chromosome 2 of S2 strain had one different residue (S32L) compared with PCS protein in 1330. In bacteria, phosphatidylcholine can be synthesized through either the phospholipid N-methylation pathway or the phosphatidylcholine synthase pathway [20, 21]. The correlation between bacterial virulence and membrane phosphatidylcholine has been reported in L. pneumophila and A. tumefaciens [22–24]. The cell envelope of the Brucella contains phosphatidylethanolamine and phosphatidylcholine. Although the phosphatidylcholine level in bacterial membrane does not affect some major virulence properties of Brucella such as invasion, intracellular traffic, and intracellular replication, it is necessary to maintain a chronic infection in mice [25, 26]. A systematic study on phosphatidylcholine synthase from S. meliloti was performed by alanine scanning mutagenesis. Protein sequence alignment of PCSs from several bacteria species showed that only five of eight residues at position 32 are serine. Changing the conserved residues close to position 32 causes about 20 % decrease in PCS activity . Further studies are needed to explore if the mutation at position 32 affects the PCS activity, which might decrease S2 virulence.
Peptidoglycan is a major component of the bacterial cell wall. Its main function is to preserve the cell integrity and resist osmotic pressure . It has been shown that peptidoglycan synthesis genes affect bacterial shape, elongation, and cell division . Therefore, disruption of peptidoglycan synthesis might affect the survival rate of the bacteria in the host. UDP-N-acetylglucosamine (UDP-GlcNAc) is one of the key building blocks in peptidoglycan biosynthesis . The second step in the formation of UDP-GlcNAc is the interconversion of glucosamine-6-phosphate to glucosamine-1-phosphate, which is catalyzed by the enzyme phosphorglucosamine mutase (PNGM). A growing body of research suggests the importance of PNGM in bacterial virulence and infectivity [31–33]. The different residue in PNGM between S2 and 1330 was located at position 250. Although the residue was close to the substrate binding site, it did not interact with the substrate. The biochemical studies are needed to fully define the function of this residue on PNGM activity and bacterial virulence.
Pyruvate kinase (PK), the key enzyme in the glycolysis, is important for controlling the concentrations of glycolytic intermediates, biosynthetic precursors and nucleoside triphosphates in the cells. Pyruvate kinase transfers the phosphoryl group of phosphoenolpyruvate to ADP to yield pyruvate and ATP. Most eukaryotic PKs can be allosterically activated by fructose 1,6-bisphosphate (FBP), while bacterial PK could be activated by either FBP or AMP/sugar monophosphates . The bacteria carrying the impaired PK activity might exhibit severe attenuation as suggested in the study in C. Neoformans . In chromosome I of S2, there was one residue different (T462A) in PK compared to that in 1330. Based on result of the Protein Databank search by blastp program, Brucella PK shared some similarities with PKs from other species, such as human, rabbit, B. Stearothermophilus, S. aureus et al. In the structure of human pyruvate kinase isozyme M2, the corresponding threonine at position 462 was close to the FBP binding site and interacted with FBP through two water molecules . The threonine residue is highly conserved in many PKs, except for the serine in E. coli. Both threonine and serine contain hydroxyl groups, which might be important for its interaction with the activator. A study demonstrates that S2 grows normally in the medium containing glucose as the only carbon source . We, therefore, suspected the mutant on PK might not be a major reason for S2 attenuation. Further biochemical studies are still necessary to clarify this conclusion.
Although Brucella spp. is considered typically as non-motile, a study demonstrates that flagellar genes are necessary for the assembly of a functional flagellum [38, 39]. B. melitensis with fliF mutant is unable to replicate in macrophages and HeLa cells, suggesting the Brucella flagellar genes might be related to virulence . Another study shows that in the early log phase of its growth in 2YT nutrient broth, B. melitensis expresses genes corresponding to the basal (MS ring) and the distal parts of the flagellar apparatus. A polar and sheathed flagellar structure has been observed under transmission electron microscopy. Although the mutations encoding different parts of the flagellum structure do not show a discernible phenotype when compared with the wild-type strain in cellular infection models, all these mutants are attenuated in mice infected via the intraperitoneal route in 4 weeks after infection . We found a mutation in fliF protein (M544L) located in MS ring. Whether the mutation decreases the virulence of S2 awaits further studies.
Other different proteins found between S2 and 1330
Two IclR regulatory genes in S2 were different from those in 1330. The IclR regulatory gene located in chromosome I lacked 8 bp that began at ATG, causing the loss of the coding protein. The other different IclR regulatory gene located in chromosome II had a 1 bp deletion, resulting in the premature stop in the coding protein. Members of the IclR transcription regulator family regulate a variety of metabolic processes in various microorganisms . In one study, B. melitensis 16M mutations in AraC, ArsR, Crp, DeoR, GntR, IclR, LysR, MerR, RpiR and TetR were constructed. GntR and LysR mutants showed the decreased virulence . Two different IclR regulatory genes found in S2 genome were correspondent to IclR5 and IclR6 in the study. Therefore, the lack of these two IclR regulatory genes might not be the major cause for S2 attenuation.
In chromosome I of S2, Omp2b had two residues different from that in 1330. Omp2b is bacterial porins present in the outer membrane of Brucella. The differences in the Omp2b variants are located in the predicted external loops of the porin . The first different residue (D79A) is highly conserved in different Brucella strains, while the second different residue (R225H) is not strictly conserved. Omp2b is present in virulent B. abortus, B. melitensis and B. suis strains that are pathogenic to humans, but is absent in non-pathogenic B. ovis . Bacterial porins have been reported to be involved in apoptosis modulation [46–48]. Expression of B. melitensis Omp2b in yeast prevents the cell from death induced by the mammalian pro-apoptotic protein Bax . No further functional studies have been reported on Brucella Omp2b.
High-affinity binding protein-dependent ABC transporter consists of a high-affinity periplasmic substrate-binding subunit, two hydrophobic membrane subunits and two additional cytoplasmic subunits. ATP hydrolysis by associated ATPase provides the energy for substrate accumulation across the inner membrane against the concentration gradient. Periplasmic substrate-binding subunits are composed of two separate but similarly folded globular domains, which are connected by a hinge region made of two or three short polypeptide segments . In S2 chromosome I, 1 bp difference in periplasmic sugar-binding protein caused one residue difference (S97Y) in protein sequence compared with that in 1330 strain. The sequence alignment with structure-solute receptor GacH of Streptomyces glaucescens (Protein Databank accession number: 3JZJ ) showed that this residue was located very close to the substrate binding site. Due to the bulky structure of tyrosine, we suspected the mutation might affect the sugar binding affinity or preference of the protein. The comparative proteome study on B. melitensis vaccine strain Rev 1 and virulent strain 16 M reveals that four sugar-binding proteins are differentially expressed in Rev 1, although these two strains have similar physical properties in biotyping tests . Therefore, it is possible that the substrate preference change might help bacteria adapt to different growth conditions.
In S2 chromosome II, 2 bp differences resulted in 2 residues difference in oligopeptide ABC transporter. The first different residue (S136G) existed in other Brucella strains, while the other different residue was only seen in S2 (P284Q). Bacteria utilize peptides as a source of amino acids, carbon, nitrogen, and/or energy. A variety of peptide uptake systems mediate the translocation of peptides across the bacterial cytoplasmic membrane. Recent studies show that OppA possesses preference in substrate peptide. E. coli OppA prefers positively charged peptides with three or four amino acids in length. However, other studies also suggest that the substrate preference of the peptide uptake systems are more complicated [53, 54]. The sequence alignment with the homologous structure of oligopeptide-binding protein from Salmonella typhimurium  showed that residue 284 was located close to the substrate binding site, but not involved in direct interaction with the substrate. Compared with other available structures of oligopeptide-binding proteins, the residue 284 of OppA protein can be glutamate (S. typhimurium), serine (E. coli), alanine (E. coli) or tyrosine (B. Pseudomallei). Meanwhile, two oligopeptide-binding proteins were encoded in 1330 genome, which shared 84 % similarity in protein sequences. The residue at 284 is aspartate in the other sequence, suggesting this residue was not conserved well in different bacteria. ATP-Binding Cassette Transporter in Brucella ovis is shown to be essential in mice . The ΔabcAB mutant strain can trigger host serologic responses similar to the WT strain and a significant cellular host response .
One residue different in S2 was found in a DedA family protein. The DedA family proteins are integral inner membrane proteins, which are present in nearly all species of bacteria. E. coli encodes 8 DedA proteins. Although each of them is nonessential, they are collectively essential . An E. coli mutant with two dedA genes deletion (ΔyghB/ΔyqjA) fails to complete cell division or grow at elevated temperatures . B. burgdorferi possesses only one dedA gene. The gene deletion results in imbalanced membrane phospholipid composition, which is required for proper cell division and envelope integrity . No functional clue was suggested for Brucella DedA so far.
We identified a difference at residue 576 between S2 and 1330 in a putative molybdopterin-binding oxidoreductase. Members of the molybdopterin oxidoreductase family include formate dehydrogenase, nitrate reductase, DMSO reductase, TMAO reductase, pyrogallol hydroxytransferase and arsenate reductase. The protein in S2 showed about 40 % similarity with formate dehydrogenase H from E. coli at residues 54–420, as suggested by the protein databank search with blastp. The residues involved in Mo binding are conserved . No significant sequence similarity was found between these two proteins around residue 576.
A 1 bp insertion caused the premature stop of 3-hydroxyisobutyrate dehydrogenase in S2. 3-hydroxyisobutyrate dehydrogenase, a key enzyme for the metabolism of valine and some keto-bodies, exists widely in bacteria, yeast, and mammalian tissue. It catalyzes the reversible conversion of 3-hydroxyisobutyrate to methylmalonate semialdehyde . Although the enzymatic properties of 3-hydroxyisobutyrate dehydrogenase have been studies in several organisms, its relevance to bacterial virulence awaits further studies.
NADH-quinone oxidoreductase subunit 1 (NuoF) in S2 had one different residue compared with that of 1330 strain. NuoF is part of the hydrophilic fragment of NADH dehydrogenase I, which represents the electron input part of the enzyme NADH: ubiquinone oxidoreductase. It couples electron transport to proton translocation across the membrane, which is partially responsible for generating the proton gradient necessary for ATP production. NADH: ubiquinone oxidoreductase serves as both a proton pump and an entry point for electrons into the respiratory chain . Currently, no bacterial virulence is associated with this gene.
One residue difference was found in lepA gene between S2 and 1330. LepA protein encoded by this gene exists in most of the bacterial genomes . Although lepA knockout strain of E. coli  and Staphylococcus aureus  are viable, the gene product might be necessary for bacteria survival under certain growth conditions .
There are other different proteins between S2 and 1330 (labeled unknown in Table 3). When we analyzed the sequences, we did not find any functional implications of the proteins in bacteria. Therefore, further studies are needed to clarify their links to virulence attenuation of S2 strain.
In this report, we sequenced the attenuated strain Brucella suis S2 genome and performed comparative genome analysis between S2 and virulent strain 1330. During the analysis, we found out 59 different ORFs between these two strains. Of these different ORFs, several have been reported to be related to Brucella virulence, such as outer membrane autotransporter, eryD. The proteins encoded by some of the different genes have been suggested to be essential for bacteria growth, such as dedA, phosphorglucosamine mutase and phosphatidylcholine synthase. These gene products might be responsible for the utilization of the materials from the environment, which are necessary for bacteria growth. We found several mutations in peptide-binding ATP transporter and sugar-binding ATP transporter, but their roles on S2 attenuation could not be ruled out without further investigation.
It has been shown that the high concentration of erythritol in the foetal tissues of cattle, sheep, goats and pigs might be beneficial for the growth of Brucella. Erythritol catabolism pathway may play important physiological roles in Brucella. The studies in S19 and S2308 showed that eryD might act as the repressor of erythritol operon . In our experiment, we found that S2 but not S19 could grow in the media supplemented with erythritol, suggesting S2 has the ability to metabolize erythritol. Therefore, further investigation on the function of the eryD gene is necessary to clarify its relevance to virulence attenuation.
Although we identified some target proteins that might be related to the virulence, most of the identified proteins have not been investigated in Brucella. Therefore, all these genes have to be examined individually to ensure their relationships with the virulence. The modification of the genes identified in our study will help to reduce residual virulence, which is critical to developing more efficient Brucella vaccine strains.
Strain and genomic DNA preparation
The S2 used in sequencing was obtained from the China Institute of Veterinary Drugs Control. Total genomic DNA was extracted using the DNeasy Blood and Tissue Kit (QIAGEN, China Ltd., China). The DNA was fragmented by nebulization, and fragments ranging from 400 to 600 bp were extracted from an agarose gel after size-fractionation, and then adaptors were added according to Illumina Sample Preparation Guide.
Genome sequencing and assembly
After sequencing with Illumina Hiseq2000, 50,941,654 paired-end (2*101 bp) reads were obtained. Base calling was performed with the software CAVASA. There were 41,608,314 paired reads passed the filter, which provided coverage of chromosome of ~1000-fold. Reads were assembled into contigs and scaffolds using SOAP de-novo (Release 1.04) . Assembled contigs were compared to the published genome sequence of B. suis strain 1330 (NC_017250.1, NC_017251.1) using BLAT v.34 . The order and the orientation of the assembled contigs were determined in accordance with the control genome. Based on the extensive similarities between the published genomes of Brucella, PCR primers were designed to link the gaps between two neighboring contigs. The leftover gaps were sequenced separately by conventional Sanger sequencing. The Illumina sequence reads for the genomes are deposited to NCBI, with the accession number of SAMN03068300.
Gene prediction and annotation
Gene prediction was established with Glimmer 2.0 using the default settings . All the predicted CDS and putative intergenic sequences were subjected to further manual inspections. Exhaustive BLAST searches with an incremental stringency against the NCBI non-redundant protein database were performed to determine the homology of the predicted coding sequences. Translational start codons were identified based on protein homology, proximity to the ribosome-binding site, relative positions to predicted signal peptide and putative promoter sequences. We also classified the putative proteins according to the COG database search results. Transfer RNAs were predicted with the tRNAscan-SE software. The repeat sequences were predicted by Trf program . The results showed 55 and 23 tandem sequences in chromosome I and II. The average copy number is 2.79 and 2.86.
All the S2 proteins predicted have been compared to KEGG bacterial gene set with KAAS (KEGG Automatic Annotation Server http://www.genome.jp/tools/kaas/) program by BBH (Bi-directional best hit) method to map the genes into bacterial pathways automatically.
Comparative genomic analysis
MUMmer 3.0 program was used to investigate the differences between the genome sequences of 1330 and S2 with the default parameters . The genomes were first automatically compared for SNPs and indels by MUMmer3.0. Then SNPs between 1330 and S2 were called by home-brew perl scripts. We classified SNPs as intergenic, frameshift and premature stop according to their positions and functions, and also counted the numbers of synonymous substitutions and non-synonymous substitutions. To identify potential genes that were different between 1330 and S2, we compared S2 genome sequences using blastp program with proteins download from NCBI RefSeq database.
S2 and S19 erythritol sensitivity test
The Trypticase-soy agar plates in the presence of 1 mg/ml erythritol [TSA (Ery)] plates were prepared . The autoclaved TSA was cooled to 47 °C, and the appropriate amount of erythritol was added and mixed. The plates were poured. S2 and S19 were grown on TSA slants at 37 °C for 24 h and resuspended in saline, and were tested by placing 20 μl of tenfold dilutions (107, 106, 105, 104, 103, 102 and 101 CFU/ml) of the bacteria resuspension on TSA and TSA (Ery) plates and incubated at 37 °C. Plates were examined for 2 to 5 days. The effect of erythritol on the growth of S2 was estimated from the size of the colonies on the test plates in comparison with the control plates (S2 and S19 inoculated on TSA plates and S19 inoculated on TSA (Ery) plates). The experiments were repeated three times.
This study was funded by the Science and Technology Basic Work Program (2012FY111000) and the National Natural Science Foundation of China (81271900). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
DDD, WXF and XM conceived and designed the experiments; DDD performed the experiments; DDD, LLT, JLK, YYH, LX, NB and YY analyzed the data; BYC, HJ, XPY, FY and QZ contributed reagents/materials/analysis tools; DDD and XM wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Pappas G, Papadimitriou P, Akritidis N, Christou L, Tsianos EV. The new global map of human brucellosis. Lancet Infect Dis. 2006;6(2):91–9.View ArticlePubMedGoogle Scholar
- O'Callaghan D, Whatmore AM. Brucella genomics as we enter the multi-genome era. Brief Funct Genomics. 2011;10(6):334–41.View ArticlePubMedGoogle Scholar
- DelVecchio VG, Kapatral V, Redkar RJ, Patra G, Mujer C, Los T, Ivanova N, Anderson I, Bhattacharyya A, Lykidis A, et al. The genome sequence of the facultative intracellular pathogen Brucella melitensis. Proc Natl Acad Sci U S A. 2002;99(1):443–8.View ArticlePubMedGoogle Scholar
- Paulsen IT, Seshadri R, Nelson KE, Eisen JA, Heidelberg JF, Read TD, Dodson RJ, Umayam L, Brinkac LM, Beanan MJ, et al. The Brucella suis genome reveals fundamental similarities between animal and plant pathogens and symbionts. Proc Natl Acad Sci U S A. 2002;99(20):13148–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Tae H, Shallom S, Settlage R, Preston D, Adams LG, Garner HR. Revised genome sequence of Brucella suis 1330. J Bacteriol. 2011;193(22):6410.View ArticlePubMedPubMed CentralGoogle Scholar
- Crasta OR, Folkerts O, Fei Z, Mane SP, Evans C, Martino-Catt S, Bricker B, Yu G, Du L, Sobral BW. Genome sequence of Brucella abortus vaccine strain S19 compared to virulent strains yields candidate virulence genes. PLoS One. 2008;3(5):e2193.View ArticlePubMedPubMed CentralGoogle Scholar
- Morgan WJ, Corbel MJ. Recommendations for the description of species and biotypes of the genus Brucella. Dev Biol Stand. 1976;31:27–37.PubMedGoogle Scholar
- Moreno E, Cloeckaert A, Moriyon I. Brucella evolution and taxonomy. Vet Microbiol. 2002;90(1–4):209–27.View ArticlePubMedGoogle Scholar
- Xin X. Orally administrable brucellosis vaccine: Brucella suis strain 2 vaccine. Vaccine. 1986;4(4):212–6.View ArticlePubMedGoogle Scholar
- Deqiu S, Donglou X, Jiming Y. Epidemiology and control of brucellosis in China. Vet Microbiol. 2002;90(1–4):165–82.View ArticlePubMedGoogle Scholar
- Verger JM, Grayon M, Zundel E, Lechopier P, Olivier-Bernardin V. Comparison of the efficacy of Brucella suis strain 2 and Brucella melitensis Rev. 1 live vaccines against a Brucella melitensis experimental infection in pregnant ewes. Vaccine. 1995;13(2):191–6.View ArticlePubMedGoogle Scholar
- Blasco JM, Marin C, De Bagues Jimenez MP, Barberan M. Efficacy of Brucella suis strain 2 vaccine against Brucella ovis in rams. Vaccine. 1993;11(13):1291–4.View ArticlePubMedGoogle Scholar
- Mustafa AA, Abusowa M. Field-oriented trial of the Chinese Brucella suis strain 2 vaccine on sheep and goats in Libya. Vet Res. 1993;24(5):422–9.PubMedGoogle Scholar
- Bosseray N, Plommet M. Brucella suis S2, brucella melitensis Rev. 1 and Brucella abortus S19 living vaccines: residual virulence and immunity induced against three Brucella species challenge strains in mice. Vaccine. 1990;8(5):462–8.View ArticlePubMedGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Bandara AB, Sriranganathan N, Schurig GG, Boyle SM. Putative outer membrane autotransporter protein influences survival of Brucella suis in BALB/c mice. Vet Microbiol. 2005;109(1–2):95–104.View ArticlePubMedGoogle Scholar
- Sangari FJ, Aguero J, Garcia-Lobo JM. The genes for erythritol catabolism are organized as an inducible operon in Brucella abortus. Microbiology. 2000;146(Pt 2):487–95.View ArticlePubMedGoogle Scholar
- Sangari FJ, Grillo MJ, De Bagues Jimenez MP, Gonzalez-Carrero MI, Garcia-Lobo JM, Blasco JM, Aguero J. The defect in the metabolism of erythritol of the Brucella abortus B19 vaccine strain is unrelated with its attenuated virulence in mice. Vaccine. 1998;16(17):1640–5.View ArticlePubMedGoogle Scholar
- Barbier T, Collard F, Zuniga-Ripa A, Moriyon I, Godard T, Becker J, Wittmann C, Van Schaftingen E, Letesson JJ. Erythritol feeds the pentose phosphate pathway via three new isomerases leading to D-erythrose-4-phosphate in Brucella. Proc Natl Acad Sci U S A. 2014;111(50):17815–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Sohlenkamp C, de Rudder KE, Rohrs V, Lopez-Lara IM, Geiger O. Cloning and characterization of the gene for phosphatidylcholine synthase. J Biol Chem. 2000;275(25):18919–25.View ArticlePubMedGoogle Scholar
- de Rudder KE, Sohlenkamp C, Geiger O. Plant-exuded choline is used for rhizobial membrane lipid biosynthesis by phosphatidylcholine synthase. J Biol Chem. 1999;274(28):20011–6.View ArticlePubMedGoogle Scholar
- Klusener S, Aktas M, Thormann KM, Wessel M, Narberhaus F. Expression and physiological relevance of Agrobacterium tumefaciens phosphatidylcholine biosynthesis genes. J Bacteriol. 2009;191(1):365–74.View ArticlePubMedGoogle Scholar
- Conover GM, Martinez-Morales F, Heidtman MI, Luo ZQ, Tang M, Chen C, Geiger O, Isberg RR. Phosphatidylcholine synthesis is required for optimal function of Legionella pneumophila virulence determinants. Cell Microbiol. 2008;10(2):514–28.PubMedGoogle Scholar
- Wessel M, Klusener S, Godeke J, Fritz C, Hacker S, Narberhaus F. Virulence of Agrobacterium tumefaciens requires phosphatidylcholine in the bacterial membrane. Mol Microbiol. 2006;62(3):906–15.View ArticlePubMedGoogle Scholar
- Comerci DJ, Altabe S, de Mendoza D, Ugalde RA. Brucella abortus synthesizes phosphatidylcholine from choline provided by the host. J Bacteriol. 2006;188(5):1929–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Conde-Alvarez R, Grillo MJ, Salcedo SP, de Miguel MJ, Fugier E, Gorvel JP, Moriyon I, Iriarte M. Synthesis of phosphatidylcholine, a typical eukaryotic phospholipid, is necessary for full virulence of the intracellular bacterial parasite Brucella abortus. Cell Microbiol. 2006;8(8):1322–35.View ArticlePubMedGoogle Scholar
- Solis-Oviedo RL, Martinez-Morales F, Geiger O, Sohlenkamp C. Functional and topological analysis of phosphatidylcholine synthase from Sinorhizobium meliloti. Biochim Biophys Acta. 2012;1821(4):573–81.View ArticlePubMedGoogle Scholar
- van Heijenoort J. Formation of the glycan chains in the synthesis of bacterial peptidoglycan. Glycobiology. 2001;11(3):25R–36R.View ArticlePubMedGoogle Scholar
- Cloud KA, Dillard JP. A lytic transglycosylase of Neisseria gonorrhoeae is involved in peptidoglycan-derived cytotoxin production. Infect Immun. 2002;70(6):2752–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Barreteau H, Kovac A, Boniface A, Sova M, Gobec S, Blanot D. Cytoplasmic steps of peptidoglycan biosynthesis. FEMS Microbiol Rev. 2008;32(2):168–207.View ArticlePubMedGoogle Scholar
- Yajima A, Takahashi Y, Shimazu K, Urano-Tashiro Y, Uchikawa Y, Karibe H, Konishi K. Contribution of phosphoglucosamine mutase to the resistance of Streptococcus gordonii DL1 to polymorphonuclear leukocyte killing. FEMS Microbiol Lett. 2009;297(2):196–202.View ArticlePubMedGoogle Scholar
- Shimazu K, Takahashi Y, Uchikawa Y, Shimazu Y, Yajima A, Takashima E, Aoba T, Konishi K. Identification of the Streptococcus gordonii glmM gene encoding phosphoglucosamine mutase and its role in bacterial cell morphology, biofilm formation, and sensitivity to antibiotics. FEMS Immunol Med Microbiol. 2008;53(2):166–77.View ArticlePubMedGoogle Scholar
- Liu XD, Duan J, Guo LH. Role of phosphoglucosamine mutase on virulence properties of Streptococcus mutans. Oral Microbiol Immunol. 2009;24(4):272–7.View ArticlePubMedGoogle Scholar
- Johnsen U, Hansen T, Schonheit P. Comparative analysis of pyruvate kinases from the hyperthermophilic archaea Archaeoglobus fulgidus, Aeropyrum pernix, and Pyrobaculum aerophilum and the hyperthermophilic bacterium Thermotoga maritima: unusual regulatory properties in hyperthermophilic archaea. J Biol Chem. 2003;278(28):25417–27.View ArticlePubMedGoogle Scholar
- Price MS, Betancourt-Quiroz M, Price JL, Toffaletti DL, Vora H, Hu G, Kronstad JW, Perfect JR. Cryptococcus neoformans requires a functional glycolytic pathway for disease but not persistence in the host. mBio. 2011;2(3):e00103–00111.View ArticlePubMedPubMed CentralGoogle Scholar
- Dombrauckas JD, Santarsiero BD, Mesecar AD. Structural basis for tumor pyruvate kinase M2 allosteric regulation and catalysis. Biochemistry. 2005;44(27):9417–29.View ArticlePubMedGoogle Scholar
- Plommet M. Minimal requirements for growth of Brucella suis and other Brucella species. Zentralblatt fur Bakteriologie: Int J Med Microbiol. 1991;275(4):436–50.View ArticleGoogle Scholar
- Sanchez DO, Zandomeni RO, Cravero S, Verdun RE, Pierrou E, Faccio P, Diaz G, Lanzavecchia S, Aguero F, Frasch AC, et al. Gene discovery through genomic sequencing of Brucella abortus. Infect Immun. 2001;69(2):865–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Halling SM. On the presence and organization of open reading frames of the nonmotile pathogen Brucella abortus similar to class II, III, and IV flagellar genes and to LcrD virulence superfamily. Microb Comp Genomics. 1998;3(1):21–9.View ArticlePubMedGoogle Scholar
- Lestrate P, Dricot A, Delrue RM, Lambert C, Martinelli V, De Bolle X, Letesson JJ, Tibor A. Attenuated signature-tagged mutagenesis mutants of Brucella melitensis identified during the acute phase of infection in mice. Infect Immun. 2003;71(12):7053–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Fretin D, Fauconnier A, Kohler S, Halling S, Leonard S, Nijskens C, Ferooz J, Lestrate P, Delrue RM, Danese I, et al. The sheathed flagellum of Brucella melitensis is involved in persistence in a murine model of infection. Cell Microbiol. 2005;7(5):687–98.View ArticlePubMedGoogle Scholar
- Molina-Henares AJ, Krell T, Eugenia Guazzaroni M, Segura A, Ramos JL. Members of the IclR family of bacterial transcriptional regulators function as activators and/or repressors. FEMS Microbiol Rev. 2006;30(2):157–86.View ArticlePubMedGoogle Scholar
- Haine V, Sinon A, Van Steen F, Rousseau S, Dozot M, Lestrate P, Lambert C, Letesson JJ, De Bolle X. Systematic targeted mutagenesis of Brucella melitensis 16M reveals a major role for GntR regulators in the control of virulence. Infect Immun. 2005;73(9):5578–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Paquet JY, Diaz MA, Genevrois S, Grayon M, Verger JM, de Bolle X, Lakey JH, Letesson JJ, Cloeckaert A. Molecular, antigenic, and functional analyses of Omp2b porin size variants of Brucella spp. J Bacteriol. 2001;183(16):4839–47.View ArticlePubMedPubMed CentralGoogle Scholar
- He Y, Xiang Z. Bioinformatics analysis of Brucella vaccines and vaccine targets using VIOLIN. Immunome Research. 2010;6 Suppl 1:S5.View ArticlePubMedPubMed CentralGoogle Scholar
- Kozjak-Pavlovic V, Dian-Lothrop EA, Meinecke M, Kepp O, Ross K, Rajalingam K, Harsman A, Hauf E, Brinkmann V, Gunther D, et al. Bacterial porin disrupts mitochondrial membrane potential and sensitizes host cells to apoptosis. PLoS Pathog. 2009;5(10):e1000629.View ArticlePubMedPubMed CentralGoogle Scholar
- Massari P, King CA, Ho AY, Wetzler LM. Neisserial PorB is translocated to the mitochondria of HeLa cells infected with Neisseria meningitidis and protects cells from apoptosis. Cell Microbiol. 2003;5(2):99–109.View ArticlePubMedGoogle Scholar
- Muller A, Gunther D, Dux F, Naumann M, Meyer TF, Rudel T. Neisserial porin (PorB) causes rapid calcium influx in target cells and induces apoptosis by the activation of cysteine proteases. EMBO J. 1999;18(2):339–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Laloux G, Deghelt M, de Barsy M, Letesson JJ, De Bolle X. Identification of the essential Brucella melitensis porin Omp2b as a suppressor of Bax-induced cell death in yeast in a genome-wide screening. PLoS One. 2010;5(10):e13274.View ArticlePubMedPubMed CentralGoogle Scholar
- Diez J, Diederichs K, Greller G, Horlacher R, Boos W, Welte W. The crystal structure of a liganded trehalose/maltose-binding protein from the hyperthermophilic Archaeon Thermococcus litoralis at 1.85 A. J Mol Biol. 2001;305(4):905–15.View ArticlePubMedGoogle Scholar
- Vahedi-Faridi A, Licht A, Bulut H, Scheffel F, Keller S, Wehmeier UF, Saenger W, Schneider E. Crystal structures of the solute receptor GacH of Streptomyces glaucescens in complex with acarbose and an acarbose homolog: comparison with the acarbose-loaded maltose-binding protein of Salmonella typhimurium. J Mol Biol. 2010;397(3):709–23.View ArticlePubMedGoogle Scholar
- Eschenbrenner M, Wagner MA, Horn TA, Kraycer JA, Mujer CV, Hagius S, Elzer P, DelVecchio VG. Comparative proteome analysis of Brucella melitensis vaccine strain Rev 1 and a virulent strain, 16M. J Bacteriol. 2002;184(18):4962–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Klepsch MM, Kovermann M, Low C, Balbach J, Permentier HP, Fusetti F, de Gier JW, Slotboom DJ, Berntsson RP. Escherichia coli peptide binding protein OppA has a preference for positively charged peptides. J Mol Biol. 2011;414(1):75–85.View ArticlePubMedGoogle Scholar
- Berntsson RP, Doeven MK, Fusetti F, Duurkens RH, Sengupta D, Marrink SJ, Thunnissen AM, Poolman B, Slotboom DJ. The structural basis for peptide selection by the transport receptor OppA. EMBO J. 2009;28(9):1332–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Tame JR, Murshudov GN, Dodson EJ, Neil TK, Dodson GG, Higgins CF, Wilkinson AJ. The structural basis of sequence-independent peptide binding by OppA protein. Science. 1994;264(5165):1578–81.View ArticlePubMedGoogle Scholar
- Silva TM, Paixao TA, Costa EA, Xavier MN, Sa JC, Moustacas VS, den Hartigh AB, Carvalho Neta AV, Oliveira SC, Tsolis R, et al. Putative ATP-binding cassette transporter is essential for Brucella ovis pathogenesis in mice. Infect Immun. 2011;79(4):1706–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Silva AP, Macedo AA, Costa LF, Turchetti AP, Bull V, Pessoa MS, Araujo MS, Nascimento EF, Martins-Filho OA, Paixao TA, et al. Brucella ovis lacking a species-specific putative ATP-binding cassette transporter is attenuated but immunogenic in rams. Vet Microbiol. 2013;167(3–4):546–53.View ArticlePubMedGoogle Scholar
- Boughner LA, Doerrler WT. Multiple deletions reveal the essentiality of the DedA membrane protein family in Escherichia coli. Microbiology. 2012;158(Pt 5):1162–71.View ArticlePubMedGoogle Scholar
- Thompkins K, Chattopadhyay B, Xiao Y, Henk MC, Doerrler WT. Temperature sensitivity and cell division defects in an Escherichia coli strain with mutations in yghB and yqjA, encoding related and conserved inner membrane proteins. J Bacteriol. 2008;190(13):4489–500.View ArticlePubMedPubMed CentralGoogle Scholar
- Liang FT, Xu Q, Sikdar R, Xiao Y, Cox JS, Doerrler WT. BB0250 of Borrelia burgdorferi is a conserved and essential inner membrane protein required for cell division. J Bacteriol. 2010;192(23):6105–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Boyington JC, Gladyshev VN, Khangulov SV, Stadtman TC, Sun PD. Crystal structure of formate dehydrogenase H: catalysis involving Mo, molybdopterin, selenocysteine, and an Fe4S4 cluster. Science. 1997;275(5304):1305–8.View ArticlePubMedGoogle Scholar
- Bannerjee D, Sanders LE, Sokatch JR. Properties of purified methylmalonate semialdehyde dehydrogenase of Pseudomonas aeruginosa. J Biol Chem. 1970;245(7):1828–35.PubMedGoogle Scholar
- Braun M, Bungert S, Friedrich T. Characterization of the overproduced NADH dehydrogenase fragment of the NADH:ubiquinone oxidoreductase (complex I) from Escherichia coli. Biochemistry. 1998;37(7):1861–7.View ArticlePubMedGoogle Scholar
- Margus T, Remm M, Tenson T. Phylogenetic distribution of translational GTPases in bacteria. BMC Genomics. 2007;8:15.View ArticlePubMedPubMed CentralGoogle Scholar
- Dibb NJ, Wolfe PB. lep operon proximal gene is not required for growth or secretion by Escherichia coli. J Bacteriol. 1986;166(1):83–7.PubMedPubMed CentralGoogle Scholar
- Colca JR, McDonald WG, Waldon DJ, Thomasco LM, Gadwood RC, Lund ET, Cavey GS, Mathews WR, Adams LD, Cecil ET, et al. Cross-linking in the living cell locates the site of action of oxazolidinone antibiotics. J Biol Chem. 2003;278(24):21972–9.View ArticlePubMedGoogle Scholar
- Bijlsma JJ, Lie ALM, Nootenboom IC, Vandenbroucke-Grauls CM, Kusters JG. Identification of loci essential for the growth of Helicobacter pylori under acidic conditions. J Infect Dis. 2000;182(5):1566–9.View ArticlePubMedGoogle Scholar
- Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, Li Y, Li S, Shan G, Kristiansen K, et al. De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010;20(2):265–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Kent WJ. BLAT--the BLAST-like alignment tool. Genome Res. 2002;12(4):656–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Delcher AL, Harmon D, Kasif S, White O, Salzberg SL. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999;27(23):4636–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27(2):573–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Delcher AL, Phillippy A, Carlton J, Salzberg SL. Fast algorithms for large-scale genome alignment and comparison. Nucleic Acids Res. 2002;30(11):2478–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Jones LM, Montgomery V, Wilson JB. Characteristics of Carbon Dioxide-Independent Cultures of Brucella Abortus Isolated from Cattle Vaccinated with Strain 19. J Infect Dis. 1965;115:312–20.View ArticlePubMedGoogle Scholar