Previous Article | Next Article ![]()
Journal of Bacteriology, April 2005, p. 2715-2726, Vol. 187, No. 8
0021-9193/05/$08.00+0 doi:10.1128/JB.187.8.2715-2726.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Bacterial Diseases of Livestock Unit, National Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Ames, Iowa,1 Department of Microbiology and Biomedical Genomics Center, University of Minnesota, Minneapolis, Minnesota2
Received 8 August 2004/ Accepted 14 January 2005
|
|
|---|
|
|
|---|
Analysis of 16S rRNA sequences places Brucella spp. as members of the alpha-2 Proteobacteria (31). The genus Brucella has six recognized species, all of which exhibit distinct host preferences (25, 26). The high degree of similarity among the brucellae (1, 3, 13, 33) lends support to the proposal that the classical species of Brucella are actually strains of Brucella melitensis (40). However, this view conflicts with the hypothesized evolutionary isolation of these classical species due to their intracellular existence and host preference (29). Common host-pathogen associations among the classical Brucella species are as follows: B. abortus, cattle; B. suis, swine; B. melitensis, goats; B. ovis, sheep; B. canis, dogs; B. neotomae, desert wood rats. Although these host-pathogen associations represent the norm in nature, cross-species infections do occur. Recently, brucellae have also been isolated from marine mammals (11, 34). Brucellae may be more widespread than previously recognized.
Brucellae can be rapidly identified by PCR assays, such as those based on the insertion sequence IS711 (6, 7). Though the sequences of brucellae are very similar, biovars of some of the classical species were differentiated by DNA sequence determination of several outer membrane proteins (OMPs) (8, 12, 41). Strains of B. abortus biovar 1 were distinguishable by analysis of a multilocus variable nucleotide tandem repeat (4).
Pulsed-field gel electrophoresis (PFGE) maps of the classical Brucella spp. genomes are composed of two circular chromosomes of approximately 2.1 and 1.2 Mbp (22, 27, 28), with the exception of B. suis biovar 3, which has a single chromosome of 3.1 Mbp. PFGE studies revealed other differences, including a 640-kb inversion in the small chromosome of B. abortus 544 and a deletion in the small chromosome of B. ovis. The two chromosomes of brucellae differ in important ways (33). The origin of replication of the large chromosome (Chr I) is typical of bacterial chromosomes, while that of the small chromosome (Chr II) is plasmid like. Further, most of the essential genes are located on Chr I. The G+C content of the two chromosomes is nearly identical, consistent with the assertion that the assimilation and stabilization of a plasmid was an ancient event (33) in brucellae.
The genome sequences of B. melitensis and B. suis have been determined (10, 33). Comparative analyses revealed both that the two genomes are extremely similar and that they have many similarities to both bacterial plant and animal pathogens and symbionts (33, 38). The sequence identity for most open reading frames (ORFs) was 99% or higher. Nevertheless, unique fragments were reported to exist between these two genomes (33). Prior to sequencing the B. abortus genome, a large number of short sequences were available in GenBank. Many of these sequences were derived from analyses of plasmids estimated to cover 20% of the genome from a random shotgun library of B. abortus S2308 (36). In this study, we present the completed B. abortus genome sequence and compare it to the genomes of B. melitensis and B. suis. Taken together, the genome sequences of these classical Brucella species provide a firm foundation for further research into the genetic bases for host preference, pathogenesis, virulence, and biotype differences.
|
|
|---|
Purification of genomic DNA. Total genomic DNA was purified from B. abortus strain 9-941 by a modification of a previously described method (17). Bacteria were harvested from agar plates in saline and killed by the addition of two volumes of methanol. Approximately 1010 bacteria were pelleted, washed in TE buffer (10 mM Tris, 1 mM EDTA; pH 8.0), and treated with 0.5% Zwittergent 3-14 in TE buffered with citrate at pH 4.0 for 1 h at 50°C. The treated bacteria were washed in TE, lysed in a solution containing 4% sarcosine, 0.5% sodium dodecyl sulfate, 125 mg of proteinase K/ml, 10 mM EDTA, and 20 mM Tris (pH 8.0) for 20 min at 65°C, and the lysate was treated with RNase A. The DNA was precipitated in ethanol, removed by spooling, and resuspended in DNAzol (catalog no. 10503-027; Life Technologies, Grand Island, N.Y.). The DNA was precipitated with ethanol a second time, dissolved in 8 mM sodium hydroxide, and adjusted to pH 7.4 with 10 mM HEPES for storage at 4°C.
Library construction. A random 2-kb insert library of B. abortus 9-941 was constructed by shearing whole genomic DNA using a nebulizer and compressed nitrogen according to protocols developed by Bruce Roe's laboratory and posted at The University of Oklahoma's Advanced Center for Genome Technology website (http://www.genome.ou.edu). The sheared DNA fragments were separated by gel electrophoresis, and fragments of 2 to 3 kb were excised from the gel and purified. The ends of the purified fragments were polished by the addition of nucleotides and Klenow fragment (New England BioLabs, Beverly, Mass.) and ligated to a SmaI-restricted calf intestinal alkaline phosphatase-treated pUC18 vector for cloning by electroporation into Escherichia coli. The library, which consisted of >90% recombinant clones, was used to construct a culture collection for sequence determination.
Sequence determination. Plasmid DNA was extracted using the QIAprep 96 Turbo kit (QIAGEN, Santa Clarita, Calif.), quantitated using PicoGreen (dsDNA quantitation kit; Molecular Probes, Eugene, Oreg.), and labeled (DyeDeoxy Terminator cycle sequencing kit; ABI automated DNA sequencing chemistry guide, ABI, Foster City, Calif.) for sequencing in the presence of dimethyl sulfoxide. The sequence was determined (ABI Prism 3700 DNA analyzer) and assembled using Phred/Phrap/Consed software obtained from the University of Washington Genome Center (http://www.genome.washington.edu/UWGC/) and the MacVector 7.0 DNA analysis package (Accelrys Inc., San Diego, Calif.). Contigs were linked and gaps were filled by predicting linkages based on putative colinearity of sequences with B. suis. Linkages were confirmed by amplifying and sequencing genomic DNA from B. abortus 9-941. The genome sequence was derived from 37,718 plasmids, and coverage was 10-fold. Mean Phrap quality score was 86. The confirmed mean read length was 819 bp.
Annotation. Artemis (releases 4 and 5; The Sanger Centre [www.sanger.ac.uk/software/ACT/]) (35) was used to identify putative genes by determining which of the ORFs with 50 or more amino acids (aa) encoded homologs in GenBank searches using BLASTP (2). B. abortus ORF annotations were modeled after those of the homologous ORFs in GenBank, especially those of the B. suis genome. B. abortus ORFs that were truncated due to premature stops or had frameshifts compared to homologs from B. suis or other entries in GenBank were designated pseudogenes. Pseudogenes of B. suis were obtained from the GenBank accession numbers for the two chromosomes of B. suis 1330. For B. melitensis, pseudogenes were identified by visual inspection of BLASTP searches of GenBank with Artemis and by cross-comparisons of the three genomes' DNA sequences using Act version 2 (The Sanger Centre).
Identification of differential ORFs. ORFs that differed in length among the Brucella genomes due to frameshifts or premature stops were labeled as differential ORFs. ORFs were not categorized as differential if their lengths differed solely due to selection of alternative start codons during the annotation processes.
Transposon search. ORFs from The Institute for Genomic Research transposon role category database of the comprehensive microbial resource website (http://www.tigr.org/tigr-scripts/CMR2/CMRHomePage.spl) for B. suis were used to identify shared transposon-related sequences among B. melitensis and B. abortus by use of MacVector 7.0.
Unique genes. Putative genes were designated as unique if no homolog was identified by aligning ORFs from the three Brucella genomes with each other (BLASTP version 2.2.6, -e 0.01, -F F) and to GenBank. Products of unique genes and ORFs were denoted as hypothetical proteins.
SNPs.
Single nucleotide polymorphisms
(SNPs) were identified using pairwise chromosome alignment data
generated with the MUMmer 3.0 run-mummer3 script
(24). No minimum
separation distance restrictions were placed on neighboring SNPs. SNP
totals were used as a measure of genetic distance for the
neighbor-joining tree [Ps =
(
SNP count/1,000)] construction with MEGA2
(23). Mesorhizobium
loti (GenBank accession numbers
NC_002678,
pA:NC_002679,
and
pB:NC_002682)
was used as an outgroup for rooting the
tree.
Nucleotide sequence accession numbers. The GenBank accession numbers for B. abortus 9-941 Chr I and Chr II are AE017223 and AE017224, respectively. The genome accession numbers used for B. suis 1330 were AE014291 and AE014292, and for B. melitensis M16 they were NC_003317 and NC_003318.
Raw
data are accessible at http://patric.vbi.vt.edu/. Note that a later
annotation of the B. melitensis genome sequence is available
at the website
http://serine.urbm.fundp.ac.be/
seqbruce/GENOMES/.
|
|
|---|
Many of the annotation differences among the genomes are related to small ORFs. The number and annotation of short hypothetical ORFs of 100 aa or less are similar between B. suis and B. abortus. The annotated genome of B. melitensis has fewer short ORFs than B. abortus and B. suis. A total of 551 ORFs of less than 100 aa are annotated in B. abortus, while in B. melitensis there are only 304. The disparity in the number of small ORFs is larger when those less than 50 aa are considered. While 161 of the 551 short ORFs in B. abortus are less than 50 aa, only 11 of the 304 short ORFs in B. melitensis are. Functional assays such as microarrays and proteomic studies will be necessary to identify genes and their products.
Occurrence of fragments in B. abortus not shared between B. suis and B. melitensis. Sequences larger than 100 bp that were previously reported to be unique in either B. suis or B. melitensis (see Table 1 in reference 33) were aligned to the B. abortus genome to determine their presence or absence (Tables 1 and 2). Many of these sequences had been shown to be related to mobile genetic elements, while others had not. Table 1 lists loci of fragments related to phages, transposable elements, or plasmids. The loci of the remaining fragments are in Table 2. Many of the fragments that were unique between B. suis and B. melitensis were found in B. abortus. The genome of B. abortus shared more fragments with B. suis and B. melitensis than B. suis and B. melitensis did with each other. B. abortus shared more fragments with B. melitensis than B. suis. This agrees with other analyses that showed B. abortus and B. melitensis being more closely related than B. abortus and B. suis (22, 27, 28).
|
View this table: [in a new window] |
TABLE 1. Occurrence
in B. abortus of fragments unique to either B. suis
or B. melitensis that have homology to phage,
transposable elements, or plasmids
|
|
View this table: [in a new window] |
TABLE 2. Occurrence
in B. abortus of regions unique between B. suis or
B. melitensis without homology to phage, transposable
elements, or plasmids
|
Regions containing sequences specific to B. abortus. The loci BruAb1_0072 and BruAb2_0168 (Fig. 1A and C) have sequences specific to B. abortus (Fig. 1B and D) relative to B. suis and B. melitensis and contain sequences that are repeated. There are eight copies of a 250-bp sequence in the 2-kb region and two copies of a 500-bp sequence in the 4-kb region. These sequences are not homologous.
![]() View larger version (45K): [in a new window] |
FIG. 1. Comparison
of region containing B. abortus-specific BruAb1_0072
and BruAb2_0168 to homologs in B. suis and B.
melitensis. (A) Representation of ORFs from B.
suis (top line), B. abortus (middle line), and B.
melitensis (bottom line) in the region with homology to
BruAb1_0072. Arrows show direction and extent of the ORF, and
hatching shows the ORF annotation. Like patterns of hatching indicate
homologs. (B) Pustell (MacVector 7.2) comparisons of
BruAb1_0072 and contiguous upstream ORF to same (1);
BruAb1_0072 and contiguous upstream to ORF BR0072 (2);
BruAb1_0072 and contiguous upstream ORF to BMEI1873 and
BMEI1872 (3); and BR0072 to BMEI1873 and BMEI1872 (4). The window size
for Pustell comparisons was 20 and minimum score is 90, with a hash
value of 6 and jump of 1. (C) Representation of ORFs from
B. suis (top line), B. abortus (middle line), and
B. melitensis (bottom line) in the region containing homology
to the locus BruAb2_0168. Arrows show directions of ORFs, and
hatching shows annotation of ORFs. Like patterns of hatching indicate
homologs. (D) Pustell comparisons of regions represented in
panel C: BruAb2_0168 to same (1); BruAb2_0168 to BR0173
(2); BruAb2_0168 to BMEII1070 and BMEII1069 (3); and BR0173 to
BMEII1070 and BMEII1069 (4). The window size for Pustell comparisons
was 20 and minimum score is 90, with a hash value of 6 and jump of
1.
|
The 4-kb region is on Chr II and encodes an autotransporter in B. abortus and B. suis, BruAb2_0168 (1,983 aa) and BRA0173 (1,113 aa), respectively (Fig. 1C). There were several possible start codons in the 6,062-bp ORF in B. abortus with homology to BRA0173. The start codon selected generating BruAb2_0168 was near the beginning of the large ORF and near a putative ribosome binding site. Experimental studies will be necessary to establish if this is an authentic start codon. In B. melitensis, a frameshift relative to the other genomes results in two relatively short ORFs, BMEII1069 (488 aa) and BMEII1070 (114 aa) (Fig. 1C). As neither of these ORFs encodes the domains of autotransporters (19), they appear to be pseudogenes. The B. abortus sequence has two direct 500-bp repeats separated by approximately 2,750 bp of sequence; the 2,750-bp sequence was B. abortus specific (Fig. 1D, graphs 1 to 3). There is a single copy of the B. abortus 500-bp repeat in B. suis (Fig. 1D, graph 2), while in B. melitensis there is only a partial copy (180 bp) (Fig. 1D, graph 3). The amino termini of BruAb2_0168 and BMEII1070 are more similar to each other than either is to the amino termini of BRA1073 (Fig. 1D, graphs 2 to 4).
The two regions encoded OMPs that were previously suggested as potentially affecting the pathogenicity or host preference of the brucellae (30).
Inversions. Sequences surrounding and at the site of the large inversion previously described in Chr II of B. abortus 544 relative to the other Brucella (27) were analyzed in B. abortus. The inversion disrupted B. abortus homologs of ORFs BRA1003 and BRA0235 of B. suis, resulting in four pseudogenes (BruAb2_0230, BruAb2_0231, BruAb2_0943, and BruAb2_0944). The ORFs BRA1003 and BRA0235 encode a putative GAF/GGDEF prokaryotic signaling domain protein and a hypothetical protein, respectively. The sites disrupted in the large inversion relative to B. melitensis are between ORFs BMEII0292 and BMEII0293 and within BMEII1009, a homolog of BRA1003. In B. abortus, a short distance downstream of the large inversion there is an indel of 838 bp. This affected a locus with homology to B. suis ORFs BRA1004 and BRA1005 and resulted in the pseudogene BruAb2_0945 in B. abortus. The finding of a single copy of the sequence 5'-CCA-GCA-CCG-CCT-GC-3' (bp 949172 to bp 949185) in B. abortus and two copies in both B. suis and B. melitensis separated by 810 bp is consistent with the indel in B. abortus arising from either homologous recombination or slipped-strand mispairing during replication. The inversion site and the 838-bp indel were described recently in B. abortus 2308 (37). The large inversion in the small chromosome is not found in all biovars of B. abortus (27). Though it was found in biovars 1, 2, 3, and 4 by PFGE, we detected it by PCR (15) only in biovars found in the United States, biovars 1, 2, and 4 (1).
Two small inversions were found in B. abortus. An inversion of 2,185 bp that was unique to B. abortus is near a 780-bp indel unique to B. suis (Table 2). This inversion occurs in a homolog of the B. suis proline dipeptidase BR1062, creating pseudogenes BruAb1_1065 and BruAb1_1067. A second smaller inversion of 2,150 bp was found in both B. abortus and B. melitensis, disrupting B. suis homologs BRA0485 and BRA0487 in these genomes. These ORFs encode a putative protein and a glycosyl transferase family 25 protein, respectively. The glycosyl transferase family 25 protein could possibly affect lipopolysaccharide structure.
Polymorphic regions encoding outer membrane proteins. Two regions in the Brucella genomes encoding homologs of OMPs predicted to be virulence associated in Brucella (30, 33) were found to have greater sequence variation than that calculated for the genomes as a whole. One of the regions encoded a putative bacterial immunoglobulin-like protein with a group 1 domain (PFAM protein family PF02369) common to bacterial surface proteins invasins and adhesins. The sequence variation affected the sizes of the homologs among the three genomes and shifted the ORF in the carboxy end in B. abortus. The OMP in B. suis, BR2009, is 500 aa, while the OMPs in B. melitensis and B. abortus, BMEI0063 and BruAb1_1984, respectively, are less than 400 aa. The sequence variation resulted in the truncation of the amino terminus of BMEI0063 relative to BR2009 and the carboxy terminus of BruAb1_1984 relative to BR2009. BR2009 and BMEI0063 have proline-rich regions in their carboxy ends. In the proline-rich stretches, 21 of 25 aa are proline in B. suis and 20 of 25 aa are proline in B. melitensis. Due to a frameshift and sequence differences in the carboxy end of BruAb1_1984 from B. abortus, there is a leucine- and histidine-rich region rather than a proline one. Only a few proline- and leucine-rich regions were found by BLASTP, and these are in disparate proteins. These adhesins may affect pathogenicity and host preference (30, 33).
The second variable region encodes an autotransporter in B. suis, BR2013, and putative pseudogenes in B. melitensis and B. abortus. This region had two in-frame stop codons in B. abortus relative to B. suis and is annotated as a pseudogene (BruAb1_1988). In B. melitensis, the homolog of BR2013 (BMEI0058) appeared to be a pseudogene also, due to an in-frame stop codon. Among the three genomes, only the B. suis locus encoded all the functional domains of an autotransporter, which may represent a unique virulence factor of B. suis.
Variable sizes of ORFs. All ORFs that varied in size among the three genomes were compiled (Tables 3 and 4). If the differences in sizes of ORFs resulted solely from selection of an alternative start codon in the annotation process, the ORFs were not labeled as being variable. The genome of B. suis was used as the reference for determining which ORFs were variable, because the annotation of its ORFs was more similar with annotation of protein homologs identified by BLASTP searches than the annotation of B. melitensis. There were almost as many variable ORFs on the large chromosomes as there were on the smaller chromosomes, even though the large chromosome has approximately twice the number of ORFs as the small chromosome. As the larger chromosome has many of the genes that encode core metabolic functions of the bacterium (33), mutations in the large chromosome may be selected against or lethal. Furthermore, the small chromosome, which appears to be a stabilized megaplasmid, may have more genes that were acquired horizontally. These genes may not be essential and, thus, might not be under positive selective pressure.
|
View this table: [in a new window] |
TABLE 3. Variable
ORFs among B. suis, B. abortus, and B.
melitensis: Chr I (BR)
|
|
View this table: [in a new window] |
TABLE 4. Variable
ORFs among B. suis, B. abortus, and B.
melitensis: Chr II (BRA)
|
Unique ORFs. Several Chr I ORFs of B. abortus had no homologs in GenBank. These ORFs are referred to here as unique ORFs and are associated with regions containing homologs of phage or insertion sequences. Unique BruAb1_1085 is near a homolog of a site-specific integrase, phage family protein. Unique BruAb1_1088 flanks a resolvase family protein. Unique BruAb1_1833 is within Tn2020 (18). Though B. suis has a homolog of this ORF, it was not annotated. Most of the unique ORFs occur in a 20-kb phage-associated fragment shared only by B. melitensis (BMEI1674-BMEI1703) and B. abortus (BruAb1_0274-BruAb1_0242). BruAb1_0272 and BMEI1774 were annotated at homologous loci but on opposing strands. BruAb1_0272 was annotated on the opposing strand because it appeared to occur in an operon. BruAb1_0246 and BruAb1_0263 in B. abortus have no protein homologs in B. melitensis. Though a number of ORFs had homology to phages, the function of the encoded peptides was rarely known. Thus, the contribution of the phage- and plasmid-related regions to brucella metabolism and infectious cycle are unknown.
Metabolic capabilities. Two genomic fragments of 7,738 and 2,653 bp identified as B. suis specific relative to B. melitensis on Chr II by Paulsen and colleagues (33) correlated with the ability of B. suis but not B. melitensis to oxidize ornithine, citrulline, arginine, and lysine and were aligned with the genome of B. abortus. Although the 2,653-bp fragment is present in B. abortus, the 7,738-bp fragment is not. Like B. melitensis, B. abortus also does not oxidize these compounds (1). This supports that the 7,738-bp fragment plays a vital role in these reactions.
B. abortus Chr I has two urease clusters, as described for B. suis (33). B. abortus ORFs BruAb1_0267-BruAb1_0273 and BruAb1_1356-BruAb1_1363 are homologs of B. suis urease cluster 1 ORFs BR0267-BR0273 and urease cluster 2 ORFs BR1356-BR1362, respectively. While there were no pseudogenes in cluster 1 of B. abortus, ureE and ureA are pseudogenes in cluster 1 in B. suis and B. melitensis. While there are no pseudogenes in cluster 2 of B. suis and B. melitensis, ureE is a pseudogene and ureD has a 6-bp insert in the urease cluster 2 in B. abortus. The sequence differences among the clusters correlate with differences in the rate of urea hydrolysis among the three bacteria. Urea is hydrolyzed in 1 to 2 h by B. abortus, compared to 0 to 30 min by B. suis (1). The rate of hydrolysis of urease in B. melitensis is variable, suggesting that at least another locus influences hydrolysis of urea. B. suis infects the urinary tract, while B. abortus and B. melitensis do not. The ability to quickly hydrolyze urea by B. suis may aid in its infection of, and excretion from, the urinary tract and subsequent spread in swine herds.
Bru-RS elements. The number and orientation relative to ORFs of the small, palindromic Bru-RS1 and Bru-RS2 elements were determined. Twenty-two whole or partial Bru-RS1 elements (14) were identified on Ch I and 18 were identified on Chr II. Short ORFs comprised largely of Bru-RS1 or Bru-RS2 elements were not annotated in B. abortus as ORFs but are in B. suis. The Bru-RS elements were not clustered, and their orientation relative to ORFs appeared to be random. A copy of the Bru-RS1 element in B. abortus was identified in one of the homologs of proline racemase (BruAb1_0363) and in a probable transcription regulator (BruAb1_1398). There were nine whole or partial Bru-RS2 elements on Chr I and five on Chr II. None of these elements occurs within genes of known function. The elements, which are just over 100 bp, could affect gene expression.
Transposable elements. There were only minor differences in the distribution and presence of transposase-related ORFs among the genomes, and those differences were associated with IS711, also known as IS6501 (16, 32). It was known from previous studies (6, 16, 32) that sequences of IS711 elements are not identical and, though the genomes have copies at the same loci, they have at least one copy at a unique locus also. The unique insertion locus of IS711 in B. suis is on Chr I, whereas the unique insertion loci of IS711 in B. melitensis and B. abortus are on Chr II (Tables 1 and 2). As described above, one of the common insertion site copies of IS711 is truncated along with ORF C of an ISBm1 copy in B. abortus relative to B. suis (BRA0551-BRA0559) and B. melitensis.
An analysis of the Brucella sequences suggests that the only mobile transposable element in Brucella is IS711. The B. abortus 9-941 genome has the same number of copies of IS711 as that found in B. abortus 544 (7). There is an additional copy in B. abortus 2308 (7), a biovar 1 strain that is commonly used as a vaccine challenge strain. The rough vaccine strain B. abortus RB51 has one more copy than its parental strain, S2308 (7, 39). Other brucellae have many more copies of IS711 than the sequenced genomes. In Southern blot analyses, B. ovis was estimated to have at least 30 copies (17), and Southern blot analyses have shown that the marine isolates have even more copies than that (5). This insertion sequence has not been identified in other bacterial genera and is most closely related to IS427 from the phylogenetically related bacterium Agrobacterium tumefaciens (32).
Sequence of field isolate versus vaccine challenge strain. The sequence of B. abortus 9-941 was analyzed for loci detected by suppressive subtractive hybridization (SSH) in B. melitensis 16 M but not B. abortus S2308 (37). Homologs of B. melitensis 16 M loci BMEI0888, BMEI0943, BMEI0971, BMEI1055, BMEI1331, and BMEI919 missing in B. abortus S2308 were also missing in B. abortus 9-941. Several homologs not identified in B. abortus S2308 in initial SSH studies were identified in the sequence of B. abortus 9-941, namely BMEI0943, BMEI0971, BMEI1055, BMEII25, BMEI1331, BMEI1380, and BME1919.
Evolutionary and genomic analyses. SNPs were identified for the shared sequences of the genomes of B. abortus, B. suis, and B. melitensis (Table 5) by using the Mummer whole genome comparison tool. This analysis identified 7,208 SNP mutations between B. abortus and B. suis genomic sequences, 6,342 SNP mutations between B. abortus and B. melitensis, and 7,844 SNP mutations between B. suis and B. melitensis. The mean for the genomes was 1 SNP for approximately 463 nucleotides.
|
View this table: [in a new window] |
TABLE 5. SNPs
among the three brucella genomes
|
![]() View larger version (14K): [in a new window] |
FIG. 2. Evolutionary
relationships among the sequenced brucellae, as shown in a
neighbor-joining dendrogram generated from an SNP-based genetic
distance data tree and rooted using the related rhizobial M.
loti genome as the outgroup. The scale bar represents 1,000
SNPs.
|
Mention of trade names or commercial products herein is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»