Previous Article | Next Article 
Journal of Bacteriology, January 2007, p. 591-602, Vol. 189, No. 2
0021-9193/07/$08.00+0 doi:10.1128/JB.01381-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Essential Bacterial Functions Encoded by Gene Pairs
,||
Helena B. Thomaides,1*
Ella J. Davison,1,
Lisa Burston,1,
Hazel Johnson,1
David R. Brown,1
Alison C. Hunt,2,¶
Jeffery Errington,2,
and
Lloyd Czaplewski1
Prolysis Ltd., Begbroke Science Park, Oxfordshire, United Kingdom,1
Sir William Dunn School of Pathology, University of Oxford, Oxford, United Kingdom2
Received 30 August 2006/
Accepted 27 October 2006

ABSTRACT
To address the need for new antibacterials, a number of bacterial
genomes have been systematically disrupted to identify essential
genes. Such programs have focused on the disruption of single
genes and may have missed functions encoded by gene pairs or
multiple genes. In this work, we hypothesized that we could
predict the identity of pairs of proteins within one organism
that have the same function. We identified 135 putative protein
pairs in
Bacillus subtilis and attempted to disrupt the genes
forming these, singly and then in pairs. The single gene disruptions
revealed new genes that could not be disrupted individually
and other genes required for growth in minimal medium or for
sporulation. The pairwise disruptions revealed seven pairs of
proteins that are likely to have the same function, as the presence
of one protein can compensate for the absence of the other.
Six of these pairs are essential for bacterial viability and
in four cases show a pattern of species conservation appropriate
for potential antibacterial development. This work highlights
the importance of combinatorial studies in understanding gene
duplication and identifying functional redundancy.

INTRODUCTION
Gene disruption is a useful tool to explore the function of
genes and to identify essential genes. Several publications
have described functional analysis of genomes by systematic
disruption of genes, particularly in prokaryotes. Genes can
be disrupted by a variety of methods, such as transposition
or the use of plasmids with conditional replication and homologous
recombination. In particular, in
Bacillus subtilis, genes that
encode essential functions were identified, while other nonessential
genes were assigned to functional categories (
56). This was
based on the examination of key cellular processes, such as
metabolism, cell division, DNA segregation, competence, or secretion,
in generated mutant strains (
56). Similarly, in
Haemophilus influenzae (
2),
Staphylococcus aureus (
35),
Salmonella enterica (
55), or
Helicobacter pylori (
98), a large number of genes that
are required for growth under nutrient-rich conditions were
identified.
So far, systematic gene disruption has been applied only to single genes. Yet, studies of genome structure have revealed the prevalence of sequence duplication (58, 92). Whole genes or parts of gene sequences are frequently duplicated. As a consequence, a number of proteins within the same organism have similar sequences. Experimental studies indicate that even if the level of sequence identity is high two proteins can have significant functional differences. For example, the type II topoisomerases gyrase and topoisomerase IV are 40 to 60% identical, although each has a distinct function essential for cell viability (17, 49). Nevertheless, the extent of sequence similarity between these essential proteins means that the same antibacterial inhibitors often inhibit both type II topoisomerases (49, 94).
Equally, pairs of proteins with little sequence similarity can nonetheless have significant functional similarity. In B. subtilis, the peptide deformylases Def and YkrB show similarity only across a short length (motif) but either is sufficient for carrying out a deformylase reaction essential for cell viability (40). That the function they carry out is essential was revealed when both genes were knocked out in combination, because the presence of one protein could compensate for the absence of the other (40). Indeed, both of these enzymes are inhibited by the same antibacterial inhibitor, despite the limited similarity in sequence and the redundancy in function (40). Our hypothesis is that many such pairs of genes, which we term isologous pairs, may have been missed by current functional genomics approaches. So far, this term has been applied to proteins in different organisms (72), but we now adopt it to describe proteins, and their encoding genes, within the same organism. Isologous proteins are similar to each other just as they are similar to their homologues in other species.
In this study, we have examined the prevalence of essential isologous proteins within B. subtilis by systematic disruption of gene pairs. Prioritization of the pairs of genes selected was based on the level of sequence identity between their encoded protein sequences and bioinformatic comparison with proteins from five other bacterial species. In addition, we characterized the growth of the mutants under several conditions. Our analysis confirms the existence of essential functions encoded by isologous gene pairs. Such pairs had not been recognized as targets for antibacterial discovery programs until recently (during the course of this study) (40, 42). Antibacterial targets tend to be single proteins that are conserved in a number of bacterial pathogens (broad spectrum), but not in humans. In an age where there is an escalating clinical need for new antibacterials, our findings could provide a new focus for future antibacterial discovery programs.

MATERIALS AND METHODS
Analysis of protein sequence similarity.
Protein sequences of the whole genomes of
B. subtilis 168,
Escherichia coli K-12,
S. aureus Mu50,
Streptococcus pneumoniae R6,
H. influenzae Rd, and
Pseudomonas aeruginosa PAO1 were obtained in FASTA format
from the NCBI (
www.ncbi.nlm.nih.gov/genomes/lproks.cgi). They
were converted into BLAST databases using the FORMATDB utility.
Searches were carried out within the total protein sequences
of one organism at a time, using the utility BLASTALL and preserving
at most the top four matches. Initially, the query sequences
and the sequences queried were the total proteins of a particular
organism. In later searches the query sequences were more specific,
while the sequences queried were the total protein sequences
of a different organism. Expectation value (E value) thresholds
were set at less than 10
10 for searches in which the
query and the result sequences were obtained from one organism
or at less than 10
3 for searches across different organisms.
The E value is a descriptor of the statistical probability that
the particular sequences of two proteins will be identical.
It is much more quantitative than the percent identity descriptor
often quoted and is therefore suitable for larger-scale analyses
such as these. Identical proteins have E values of 10
180,
while BLAST does not recognize alignments with E values greater
than 10
3.
Computer scripts (written in Perl) were used to extract the protein identifiers for each query sequence that had found at least one matching sequence (hit), the E value of the match, and the protein identifiers of the best hits resulting from the search. For the searches within one organism, only information on searches in which the query sequences had found one matching sequence (other than themselves) was extracted. This is essentially a modified version of the bidirectional best-hit strategy (31). For the searches across organisms, the best two hit sequences were extracted from the search. Hit sequences resulting from such searches were used as query sequences in searches of the organism from which the original query sequences had been derived (reciprocal searches). Data extracted from the results of both searches were analyzed, using Microsoft Access and Excel, to identify similar proteins. For example, if four proteins (two pairs, each from a separate organism) found each other as top hits in the four corresponding BLAST searches (one within each organism and two across organisms), then the proteins were considered similar.
Creation of pHT35.
Standard methods for cloning and PCR were used (99). The gus neo fragment was extracted from pMLK83 (53) digested with NotI. This fragment was ligated to pACYC184 (93) that had been cut with SalI and BamHI and blunted, thus disrupting the tetracycline gene and resulting in vector pHT32. Additional cloning sites were created upstream of the Pxyl promoter of pRD96 (21) by the annealing of complementary oligonucleotide sequences topL (5'-AGCTATCGATGGATCCACTAGTAAGCTTACTGCA-3') and botL (5'-TAGCTACCTA GGTGATCATTCGAATG-3') and the ligation of the resulting linker molecule to pRD96 cut with PstI. The cat Pxyl polylinker fragment from the resulting recombinant vector was extracted following ClaI and BglII digestion and blunting at both ends. This fragment was ligated to pHT32 that had been cut with HindIII and ClaI and then blunted at both ends to create pHT34. Primers 34-F (5'-ACGGGATCCTGGCGACCACACCCGTCCTGTGGG-3') and 34-R (5'-CGATGGATCCACTAGTAAGCTTACTGCAGG-3') were used in a PCR to introduce a second BamHI site into pHT34 upstream of the gusA ribosome binding site. This vector was digested with BamHI and reannealed, selecting for loss of a 323-bp fragment and resulting in pHT35 (see Table S1 in the supplemental material).
Plasmid constructs for disruption and conditional mutants.
PCR was used to amplify an
300-bp fragment from each gene of the putative pair, starting from approximately 200 bp within the gene for the generation of knockout constructs or from immediately upstream of the predicted ribosome binding site for the generation of conditional constructs. Amplification primers were designed to contain a restriction site for subsequent cloning into a relevant vector (pMUTIN4 [119] for one member of the potential pair and pHT35 for the other). Following digestion of both the insert and the vector with the relevant restriction enzymes, insert and vector were ligated. The ligation mixture was used to transform competent E. coli MC1061 (12) cells. Transformants were selected on 25 µg/ml chloramphenicol (pHT35) or 100 µg/ml ampicillin (pMUTIN4). Insert-positive constructs were identified by PCR using vector-specific primers. Plasmids isolated from insert-positive E. coli colonies (see Table S1 in the supplemental material) were sequenced and sequence data checked using Vector NTI software.
Creation of fabHA conditional mutant.
Primers 99BI-F (5'-TCGAAGCTTAAGATATAAAGGAGTC-3') and 99BI-R (5'-TTAGGGATCCTGTGTGCACCTCACCTT-3') were used to amplify the entire fabHA gene and to introduce HindIII and BamHI primers, respectively. The amplified and digested product was then ligated to similarly digested pJPR1 (J. Rawlins, unpublished). The resulting plasmid was pHT48.
Creation of polA conditional mutant.
Primers polAmis-F (5'-GAAGGAAATTGATTGACGGAACGAAAAAAATTAGTGCTTGTAG-3') and polAmis-R (5'-CAATTTCCTTCTAAGAAGCCTCAAGCTTAATTG-3') were used to introduce a mutation of the start codon of polA, from ATG to TTG, into plasmid pLI135B (see Table S1 in the supplemental material) by amplifying around the plasmid by PCR. This mutation was necessary because the original efficient translation signal seemed to give PolA product sufficient to allow growth even under repressive conditions.
Creation of polA and ypcP deletion mutants.
To construct a complete deletion of PolA, primers prepolA-F2 (5'-CGAGGGCCCAAAAGGCGAATTGTGGGGGATGCTG-3') and prepolNK-R (5'-TTTAAGCTTCCGTCATTCAATTTCCTCC-3') were used to amplify a fragment upstream of polA by PCR and to introduce a HindIII site at the downstream end for cloning. Primers postpolNK-F (5'-GTAGGATCCGAAATAAACAGAGATAGG-3') and postpolA-R (5'-CAAGCATGCGATTGAGCATTTCCTGCCTGAC-3') were used in a PCR to amplify a fragment downstream of polA and introduce a BamHI site at the upstream end for cloning. The kanamycin resistance cassette from pBEST501 (4) was extracted following digestion with HindIII and BamHI. The cassette was ligated to the two PCR fragments, which had been digested with HindIII and BamHI. Primers postpolA-R and prepolA-F2 were used to amplify the ligation product by PCR.
To obtain a C-terminal deletion of PolA, primers prepoltrn-F (5'-GACGGATCCAAAGGCGAATTGTGGGGGATGCTG-3') and prepoltrn-R (5'-ATCGGTACCTCTACTCTTCATGATAATTATCGCCAATCTGTTCAACG-3') were used to amplify a fragment within polA and introduce BamHI and KpnI sites for cloning. The PCR product, digested with BamHI and KpnI, was ligated to pBEST501 that had been cut similarly. The resulting plasmid was then digested with SphI and SalI and ligated to a fragment downstream of polA that had been amplified by PCR, using primers postpolAtrn-F (5'-GTAGCATGCGAAATAAACAGAGATAG-3') and postpoltrn-R (5'-GGGGTCGACAATTGCATTTAAAGCGAGACG-3'), and digested with SphI and SalI. The resulting plasmid was pED5-2 (see Table S1 in the supplemental material).
To construct a complete deletion of ypcP, primers preYpcPN-F (5'-CAGAAGCTTAAACAGTGACCGATTATGTCAGCG-3') and preYpcPN-R (5'-GCAAAGATCTATTATTATTCATCAAATAACTCC-3') were used to amplify a fragment upstream of ypcP and introduce cloning sites by PCR. The fragment was digested with HindIII and BglII and ligated to similarly cut pRD96. The resulting vector was digested with BamHI and ligated to a fragment downstream of ypcP, which had been amplified by PCR using postYpcPN-F (5'-GCTGGATCCTAGAGAGATCGTTTAGATCC-3') and postYpcPN-R (5'-GACGGATCCCGTTCCTTTTGCTTCTATTTTCAAATGCG-3') primers and cut with BamHI. The resulting vector was partially cut with BamHI and EcoRI. Its staggered ends were blunted and ligated. The resulting vector was pED4-2 (see Table S1 in the supplemental material).
Cloning in B. subtilis.
For single gene disruptions, deletions, and creation of conditional mutants (see Table S1 in the supplemental material), B. subtilis 168CA was made competent (59, 70) and transformed with recombinant plasmids or linear DNA fragments, created as described above. Selection was for resistance to kanamycin (5 µg/ml) in the presence of 0.5% xylose (pHT35 and pBEST501) or to erythromycin (0.5 µg/ml) in the presence of 1 mM isopropyl-ß-D-thiogalactoside (IPTG) (pMUTIN4) or resistance to chloramphenicol (5 µg/ml) (pRD96 and pJPR1). To create the C-terminal PolA deletion, loss of the original selective marker, the kanamycin cassette, was achieved by growth for several generations in the absence of antibiotic. To create double gene knockouts, a single knockout strain was made competent and transformed with chromosomal DNA from the strain with a knockout of its gene partner, selecting only for the incoming mutation. Resulting colonies were patched onto separate antibiotic plates to check for retention of the resident mutation in the presence and absence of inducer. Colonies were also diagnosed based on two PCRs; one reaction used two plasmid-specific primers and one insert amplification primer to diagnose integration of the plasmid, and the other used two primers in regions outside the insert to diagnose disruption of the target sequence.
Analysis of medium dependence.
Single knockouts, double knockouts, and conditional strains were all streaked on nutrient agar or minimal agar (as described in reference 39, but without trace elements) media in the presence or absence of appropriate antibiotics and inducers. Selective conditions and colonies were examined by eye after at least 16 h at 37°C.
Some strains were also examined in liquid media, either Oxoid Penassay broth or minimal medium (as described in reference 39, but without trace elements), by diluting cultures 1/1,500 and incubating them at 37°C in the presence (1% xylose or 1 mM IPTG, as appropriate) or absence of inducer. For measurement of reporter activity, samples were incubated with assay buffer containing lysozyme, 4-methyl-umbelliferyl ß-D-galactoside, resorufin ß-D-glucuronide, and Triton X-100 for 1 h at room temperature before measurement of fluorescence on a BMG Fluostar Galaxy instrument. To rescue growth on minimal medium, 200 µg/ml of different amino acids was added separately and in combination.

RESULTS AND DISCUSSION
Protein duplication in bacteria.
The
B. subtilis genome contains a high level of sequence duplication
that does not appear to result from duplication of large segments
of DNA (
60). To quantify the extent of protein duplication,
BLASTP was used to compare the gene products of the whole
B. subtilis genome against itself. Seven hundred two proteins that
could be reciprocally paired (351 pairs) were identified, based
on an average threshold E value for the reciprocal similarities
of

10
20. This cutoff is reasonably, but not greatly,
different from the minimum level of similarity recognized in
BLASTP searches of 10
3 (Fig.
1). Triplets and larger
groupings of similar proteins were excluded from this analysis.
Analysis of five other bacterial genomes also revealed large
numbers of duplicated proteins in each of them (Fig.
1; Table
1). Of the organisms analyzed,
H. influenzae appeared to contain
the least redundancy. However, of the 194 proteins (97 pairs)
in
H. influenzae, 130 had previously been reported to have essential
functions (
2). Such analyses confirmed that pairs of similar
proteins exist in a number of bacterial species and as such
were worth experimental investigation to find uncharacterized
essential functions.
Systematic disruption of 135 gene pairs.
Several groups have completed functional analysis programs based
on high-throughput disruption of single genes and aimed at finding
essential genes (
2,
35,
55-
56,
98). However, essential functions
protected by redundancy, so that the loss of one protein can
be compensated for by the activity of another, would not have
been detected. We therefore undertook a systematic analysis
of gene pairs in
B. subtilis.
Bioinformatic analysis suggested that extending the threshold E value to 1010 would include protein pairs whose sequence similarity is restricted to short motifs. We therefore identified a total of 523 pairs in B. subtilis, and to these a number of criteria were applied to select gene pairs for experimental evaluation. These pairs loosely fell into three categories. In the first category were pairs encoding proteins with high similarity. Seventy-one pairs were chosen based on an arbitrary cutoff value (E value of
1070, equivalent to an observed >26% sequence identity). The second category consisted of 61 pairs selected from the remaining 452 gene pairs encoding proteins connected with lower similarity scores (E values of 1010 to 1069). This category had been filtered for the proteins that were duplicated in at least one more bacterial species, on the simple basis that they were more likely to contain common motifs. Within these we did discern some proteins with characterized motifs. In addition, these 61 pairs lacked protein homologues in humans and therefore represented higher-priority antibacterial targets. Finally, we included protein pairs that were also significantly similar to other proteins in B. subtilis and so were parts of larger similarity groupings but were putative homologues of protein pairs identified in H. influenzae. The three pairs in this third category were included to ensure coverage of duplication in pathogenic bacteria. In total, 135 pairs of genes were selected for systematic disruption individually and in combination (Table 2).
Examination of the 270 genes making up these pairs indicated
that 123 had been functionally characterized to some extent
elsewhere (Table
2) (
1,
5-
11,
13-
15,
18-
19,
22,
24,
28-
30,
32-
34,
37-
38,
41-
48,
51-
52,
54,
57,
61-
69,
71,
73-
81,
83-
89,
91,
95-
97,
101-
118,
120-
128). For the majority of the pairs, only one gene
had been characterized. Where both genes had been characterized,
the paired genes appeared to have a function in common (28 pairs)
(
1,
5,
7-
9,
11,
15,
29,
32-
33,
37,
42,
44-
45,
54,
57,
65,
74,
83-
84,
86,
88,
97,
105,
108,
110,
114-
117,
123-
124,
126). Primarily,
the characterization for these pairs was either incomplete or
based solely on biochemical data. This made the need for an
in vivo study more relevant. Therefore, all of the genes selected
bioinformatically were included in our experimental analysis.
To make double disruption strains readily, two vectors were used in parallel. One was pMUTIN4 (119), previously used for the functional analysis of single genes in B. subtilis (56). The other vector, pHT35, was constructed to be pMUTIN4-like. The selective marker (neo), reporter gene (gus), and inducible promoter (Pxyl) of the new vector, pHT35, were all different from those of pMUTIN4. pHT35 was based on a nonhomologous E. coli replicon, P15A (93), to ensure that its DNA sequence was unrelated to that of pMUTIN4 (Fig. 2). This plasmid combination should be valuable in many situations in which expression of two genes needs to be manipulated simultaneously and/or measured in parallel.
Of the 270 genes analyzed, 244 were successfully disrupted (Table
2). Our results were in agreement with earlier publications,
as only two (
topA and
yqiD) of the 244 genes we disrupted had
previously been reported to be essential (
56).
Identification of new essential genes.
The remaining 26 genes could not be disrupted, despite multiple attempts (Table 2). Eleven of these genes (ctaD, gapA, mntA, yfkN, ysgA, ycgG, yloU, yddT, yomL, ytbE, and yqhY) had previously been reported to be nonessential (14, 32, 56, 88, 120). Five of these genes (yloU, yqhY, gapA, ysgA, and ycgG) had, however, been defined as "persistent" based on the strong conservation of homologues in gram-positive bacteria (31). We propose that these genes should be added to the list of essential genes in B. subtilis. Of these 11, we have prioritized mntA, ctaD, ysgA, ycgG, yloU, and yfkN for future characterization as potential antibacterial targets because they are widely conserved and absent from humans.
Other single gene effects.
The mutants of 11 nonessential genes were auxotrophic (Table 3). Eight of these genes (yjbV, thiE, pyrC, purK, gltA, tenI, hisF, and hisA) had known metabolic roles (61, 64, 82, 83, 90, 104) or were significantly similar to characterized metabolic proteins in other species (25). An additional set of 13 mutants showed defective sporulation (Table 3). Seven of these genes (nasA, oppA, yvgN, ytxN, speA, opuE, and recG) had been characterized previously at least partly (3, 80, 85, 96, 109, 113, 122), but only oppA and ytxN had been implicated directly in sporulation (85, 113). These observations provided preliminary information on the possible general functions of ysaA, ypvA, yfjR, yodR, yqiD, yjbM, and ywaC, genes that were previously without functional assignment.
Remarkably, the effects of one single knockout were suppressed
by disruption of the partner gene (Table
3). Disruption of
yodR blocked sporulation, but when a
yodR mutation was combined with
a
scoB mutation the double knockout sporulated as well as the
wild type. Whereas the double mutant had approximately 70% spores
after a 2-day incubation on nutrient agar, no spores were detected
with the
yodR mutant. Similarly, the
yodR mutant was unable
to grow in minimal medium, while the double mutant could (data
not shown). All three strains grew well in rich medium. Bioinformatic
analysis suggests that both genes encode probable coenzyme A
transferases but with different substrates. Indeed the two proteins
are reasonably similar (E value of 10
53 or 50% identity
across their length). We speculate that this phenotypic effect
may be indirect and a consequence of variation in levels of
different metabolites in these strains. Alternatively, ScoB
may be partially isologous to YodR, causing the protein substitution
to be detrimental under some growth conditions.
Auxotrophic isologous pairs.
Knockouts of the nonessential genes, comprising 114 pairs in total, were combined to produce double mutants. This resulted in 108 viable double knockouts (Table 2). One hundred one of these mutants had not been examined previously. Our results also confirmed the nonessentiality of six pairs characterized in vivo elsewhere (29, 33, 45, 65, 97). There is controversy over the essentiality of the ykuA pbpA pair (124). All of the double knockout mutants were examined for defects in growth and sporulation on both rich and minimal media.
One double knockout (ywaA ybgE) showed medium-dependent growth. For the ywaA ybgE pair, the double knockout strain showed good growth in rich medium, whereas in minimal medium the mutant showed little growth over 5 h (Fig. 3). In contrast, the single mutants grew well under both conditions, indicating that YwaA and YbgE are isologous, with a common function essential only in minimal medium. Bioinformatic analysis suggested that YwaA and YbgE were likely branched-chain amino acid aminotransferases. Indeed, when minimal medium was supplemented with both isoleucine and valine, growth of the double knockout strain was rescued, confirming an effect on branched-chain amino acid metabolism (Fig. 3). The presence of a leucine dehydrogenase (Bcd) in B. subtilis (23) may explain why the amino acid requirements we observed in vivo did not include leucine, the sole other branched-chain amino acid. During the period of this study, YwaA and YbgE were shown to catalyze an aminotransferase reaction with these amino acids as substrates (rather than products) for the synthesis of methionine (8). This is only partially consistent with our data as, in vivo, methionine did not rescue activity (data not shown).
Essential isologous pairs.
Consistent with our original hypothesis, we also identified
six pairs of genes that appeared to be essential when both genes
were inactivated (Table
4). In these instances, the genes in
the pair could be knocked out individually but not in combination.
Within these pairs, the genes were far enough apart on the chromosome
to preclude gene linkage as an explanation to these results.
To validate that these pairs were essential, for each pair one
gene was disrupted and its partner was made conditional. As
shown in Fig.
4, for all six pairs, repressing the expression
of one gene led to a decrease in growth, thus confirming that
the proteins were both essential and isologous. Their characteristics
will now be described in turn.
ThrS and ThrZ.
ThrS and ThrZ have been described previously as isologous threonyl
tRNA synthetases (
36,
86). We also find that the isologous pair
is essential for cell growth (Fig.
4) (
86). They share the strongest
level of similarity of all of the essential pairs newly identified.
Most bacteria outside the
Bacillus spp. possess a single gene
encoding threonyl tRNA synthetase, rather than two genes. Similar
threonyl tRNA synthetases exist in humans, making these proteins
less attractive targets for antibacterial discovery.
TopA and TopB.
Most bacteria have either topoisomerase I (TopA) or topoisomerase III (TopB), but B. subtilis, E. coli, and S. aureus are relatively unusual in having both proteins. In B. subtilis, one of the proteins, TopA, was considered to be essential (56). Based on our results, the proteins in B. subtilis are isologous and essential as a pair. It is likely that so are those in E. coli and S. aureus. Owing to sequence similarity to human topoisomerases, TopA and TopB are lower-priority antibacterial targets.
FabI and FabL.
The genes coding for enoyl reductases FabI and FabL were shown elsewhere during the course of this study to have an essential function in fatty acid biosynthesis, so that the absence of one protein can be compensated for by the presence of the other (42). Our results corroborate that these proteins are isologous and essential as a pair (Fig. 4; Table 4). Interestingly, depletion of both B. subtilis proteins was shown to increase expression of fabI, whereas depletion of either protein alone caused no effect on fabI expression (Fig. 5). This feedback effect confirmed that both enoyl reductases are active in vivo. They are conserved as a single enzyme in E. coli, H. influenzae, Enterococcus faecalis, and S. aureus. They are absent from humans, which makes them potential antibacterial targets, as validated by the antibacterial triclosan (42).
FabHA and FabHB.
For FabHA and FabHB, inducer dependence was seen only with
fabHA as the conditional allele (not shown). Presumably, the repressible
construct for
fabHB allows sufficient residual expression to
allow growth. FabHA and FabHB are therefore isologous, and biochemical
data (
15,
20) suggest that they are 3-oxoacyl-[acyl carrier
protein] synthases involved in fatty acid biosynthesis. They
are duplicated in
P. aeruginosa. They are conserved as single
enzymes in
E. coli,
S. aureus,
S. pneumoniae,
E. faecalis, and
H. influenzae and absent from humans, so this enzyme is a good
antibacterial target. Indeed, there are reports of antibacterial
discovery programs directed at the single enzyme in
E. coli (
20).
SpoIIIJ and YqjG.
Homologues of SpoIIIJ and YqjG are found in a wide range of other bacteria, usually as a single protein but as a pair in Streptococcus spp. Various studies with B. subtilis and E. coli suggest that the protein function is in secretion (100, 116). As no homologues exist in humans, this pair may have potential as a novel antibacterial target.
PolA and YpcP.
Among the identified isologous pairs, the pair of PolA and YpcP is the most novel and validates the rationale for this study. Previous functional genomics studies have suggested that DNA polymerase I (PolA) is essential in S. pneumoniae (27), E. coli (50), and S. enterica (55) but not in B. subtilis (56), S. aureus (35), or H. influenzae (2). We have found a possible explanation for this in that, in B. subtilis at least, polA becomes essential when ypcP is deleted.
YpcP is a small protein of previously unknown function, and its sequence similarity is restricted to the N-terminal domain of PolA, which contains a 5'-to-3' exonuclease motif shown to be functional in S. pneumoniae and E. coli (26, 129) (Fig. 6). Standard single crossover insertions disrupting ypcP (with pMUTIN4) and polA (with pHT35) could not be combined. Similarly, complete deletions of each gene, though viable singly, could not be combined. However, a disruption that truncated PolA after the end of the N-terminal exonuclease domain (after codon 334) was viable and grew normally (data not shown). This and the depletion experiment mentioned above (Fig. 4) demonstrate that this putative 5'-to-3' exonuclease domain is essential and that its function can be supplied either by the N-terminal part of PolA or by YpcP. Indeed, we have biochemical evidence, with a His tag fusion to YpcP, that confirms that exonuclease activity with a preference for double-stranded DNA is the essential function (E. Davison and H. Thomaides, unpublished data).
These data may explain the previous findings on the essentiality
of PolA in some of the other species. They also suggest strongly
that in these other organisms it is the 5'-to-3' exonuclease
domain of PolA that is essential. This domain duplication is
present in various pathogens, including
S. aureus,
Mycobacterium tuberculosis,
S. enterica, and
E. coli (Fig.
6). Experimental
evidence for some of these pathogens, such as
E. coli and
S. enterica (
50,
55), suggests that although the domain is duplicated
the single PolA remains essential. In other pathogens, only
the single PolA protein is present, and in the case of
S. pneumoniae,
at least, it is essential (
27). Importantly, this domain appears
to be completely absent in humans. PolA/YpcP therefore represents
another unexploited antibacterial target that needs to be characterized
further.
Perspective.
This study opens many avenues for future work: 11 new essential genes have been identified and need characterization, and several previously uncharacterized genes have now been shown to be needed for growth in minimal medium. There is a whole new set of sporulation genes to be investigated. Some of these genes had initially been paired together through our bioinformatic work (e.g., yurU and yurX, yloU and yqhY, and tenI and thiE). Consequently, those that were found to be essential may encode proteins similar enough as to be inhibited by the same antibacterial, as is the case for gyrase and topoisomerase IV (94). Importantly, we have confirmed the status of three essential isologous pairs in B. subtilis (ThrS/ThrZ, FabI/FabL, and SpoIIIJ/YqjG) and identified four new ones (YwaA/YbgE, TopA/TopB, FabHA/FabHB, and PolA/YpcP) (Table 4). An extension of this work could be the development of inhibitors against SpoIIIJ/YqjG and PolA/YpcP, the unexploited isologous pairs that have no homologues in eukaryotes and therefore have potential as targets for antibacterials.
The isologous proteins do not appear to have any features, such as function, level of sequence similarity, and species conservation, distinguishing them from other proteins we paired. Although whole-genome bioinformatic analyses suggest a clustering of gene duplication (16, 60), we did not observe any clustering of the isologous genes on the chromosome. Interestingly, in each of four of the isologous pairs one gene has been identified as being "persistent" within the gram-positive clade of Firmicutes (31). This bioinformatic definition extends the number of essential bacterial functions that cannot be shown experimentally (31). The def ykrB essential pair identified elsewhere also has only one persistent gene (31, 40). That only one gene in an essential pair "persists" is consistent with our observation that the extent to which an essential function is protected by redundancy (e.g., via isologous pairs) is species specific. From our analysis, the two essential pairs involved in fatty acid biosynthesis (fabHA fabHB and fabI fabL) and the auxotrophic pair (ywaA ybgE) were not persistent, but then neither were many essential single genes (Table 2) (31). Similar analysis of the remaining pairs in this work (data not shown) suggests that an experimental approach is most suited to identify essential isologous proteins, though limited by chosen experimental conditions.
Why duplication in genomes is maintained is a heavily debated issue. We speculate that this systematic study has revealed the purpose of some of the redundancy within the B. subtilis genome by identifying some of the isologous pairs. Isologous proteins directly contribute to the fitness of the organism, while perhaps offering the possibility of acquiring mutations that fine-tune or extend protein functions. The experimental evaluation of pairs in other species or of triplets or quadruplets of genes would also investigate this issue; we detected such groupings in our preliminary analysis (D. Brown, unpublished data).
Our work sets a precedent for the systematic identification of essential isologous proteins. Similar studies will enhance our understanding of the biological significance of sequence duplication.

ACKNOWLEDGMENTS
We thank the Biotechnology and Biological Sciences Research
Council and the Department of Trade and Industry, United Kingdom,
for funding of this work through the award of a LINK grant in
Applied Genomics.
We thank Celia Caulcott for her management of the LINK project. We thank Steve Ruston of Prolysis Ltd. for his support. We thank various colleagues, Neil Stokes in particular, for critical reading of the manuscript.

FOOTNOTES
* Corresponding author. Mailing address: Prolysis Ltd., Begbroke Science Park, Sandy Lane, Yarnton OX5 1PF, Oxfordshire, United Kingdom. Phone: 44 1865 854700. Fax: 44 1865 854799. E-mail:
helena.thomaides{at}prolysis.com.

Published ahead of print on 17 November 2006. 
|| Supplemental material for this article may be found at http://jb.asm.org/. 
Present address: OCTO, Department of Clinical Pharmacology, Radcliffe Infirmary, Woodstock Road, Oxford OX2 6HE, United Kingdom. 
Present address: GPAR, Building 20 GlaxoSmithKline, Greenford Road, Greenford, Middlesex UB6 0HE, United Kingdom. 
¶ Present address: Department of Medical Microbiology, Aberdeen Royal Infirmary, Foresterhill, Aberdeen AB25 2ZN, United Kingdom. 
Present address: Institute for Cell and Molecular Biosciences, The Medical School, University of Newcastle, Catherine Cookson Building, Framlington Place, Newcastle NE2 4HH, United Kingdom. 

REFERENCES
1 - Ahmed, M., L. Lyass, P. N. Markham, S. S. Taylor, N. Vazquez-Laslop, and A. A. Neyfakh. 1995. Two highly similar multidrug transporters of Bacillus subtilis whose expression is differentially regulated. J. Bacteriol. 177:3904-3910.[Abstract/Free Full Text]
2 - Akerley, B. J., E. J. Rubin, V. L. Novick, K. Amaya, N. Judson, and J. J. Mekalanos. 2002. A genome-scale analysis for identification of genes required for growth or survival of Haemophilus influenzae. Proc. Natl. Acad. Sci. USA 99:966-971.[Abstract/Free Full Text]
3 - Alonso, J. C., G. Luder, and R. H. Tailor. 1991. Characterization of Bacillus subtilis recombination pathways. J. Bacteriol. 173:3977-3980.[Abstract/Free Full Text]
4 - Autret, S., R. Nair, and J. Errington. 2001. Genetic analysis of the chromosome segregation protein Spo0J of Bacillus subtilis: evidence for separate domains involved in DNA binding and interactions with Soj protein. Mol. Microbiol. 41:743-755.[CrossRef][Medline]
5 - Ayora, S., F. Rojo, N. Ogasawara, S. Nakai, and J. C. Alonso. 1996. The Mfd protein of Bacillus subtilis 168 is involved in both transcription-coupled DNA repair and DNA recombination. J. Mol. Biol. 256:301-318.[CrossRef][Medline]
6 - Balassa, G., B. Dod, V. Jeannoda, P. Milhaud, J. Zucca, J. C. Sousa, and M. T. Silva. 1978. Pleiotropic control mutations affecting the sporulation of Bacillus subtilis. Ann. Microbiol. (Paris) 129:537-549.
7 - Baranova, N. N., A. Danchin, and A. A. Neyfakh. 1999. Mta, a global MerR-type regulator of the Bacillus subtilis multidrug-efflux transporters. Mol. Microbiol. 31:1549-1559.[CrossRef][Medline]
8 - Berger, B. J., S. English, G. Chan, and M. H. Knodel. 2003. Methionine regeneration and aminotransferases in Bacillus subtilis, Bacillus cereus, and Bacillus anthracis. J. Bacteriol. 185:2418-2431.[Abstract/Free Full Text]
9 - Bookstein, C., C. W. Edwards, N. V. Kapp, and F. M. Hulett. 1990. The Bacillus subtilis 168 alkaline phosphatase III gene: impact of a phoAIII mutation on total alkaline phosphatase synthesis. J. Bacteriol. 172:3730-3737.[Abstract/Free Full Text]
10 - Bower, S., J. B. Perkins, R. R. Yocum, C. L. Howitt, P. Rahaim, and J. Pero. 1996. Cloning, sequencing, and characterization of the Bacillus subtilis biotin biosynthetic operon. J. Bacteriol. 178:4122-4130.[Abstract/Free Full Text]
11 - Bunai, K., M. Ariga, T. Inoue, M. Nozaki, S. Ogane, H. Kakeshita, T. Nemoto, H. Nakanishi, and K. Yamane. 2004. Profiling and comprehensive expression analysis of ABC transporter solute-binding proteins of Bacillus subtilis membrane based on a proteomic approach. Electrophoresis 25:141-155.[CrossRef][Medline]
12 - Casadaban, M. J., and S. N. Cohen. 1980. Analysis of gene control signals by DNA fusion and cloning in Escherichia coli. J. Mol. Biol. 138:179-207.[CrossRef][Medline]
13 - Cha, M. H., W. M. Yong, S. M. Lee, Y. S. Lee, and I. Y. Chung. 2000. The biochemical and molecular characterization of recombinant Bacillus subtilis tripeptidase (PepT) as a zinc-dependent metalloenzyme. Mol. Cells 10:423-431.[Medline]
14 - Chambert, R., Y. Pereira, and M. F. Petit-Glatron. 2003. Purification and characterization of YfkN, a trifunctional nucleotide phosphoesterase secreted by Bacillus subtilis. J. Biochem. (Tokyo) 134:655-660.[Abstract/Free Full Text]
15 - Choi, K. H., R. J. Heath, and C. O. Rock. 2000. Beta-ketoacyl-acyl carrier protein synthase III (FabH) is a determining factor in branched-chain fatty acid biosynthesis. J. Bacteriol. 182:365-370.[Abstract/Free Full Text]
16 - Clifton, S., and B. A. Roe. 1998. The beautifully simple but intriguingly complex world of small genomes. Genome Res. 8:331-333.[Free Full Text]
17 - Corbett, K. D., A. J. Schoeffler, N. D. Thomsen, and J. M. Berger. 2005. The structural basis for substrate specificity in DNA topoisomerase IV. J. Mol. Biol. 351:545-561.[CrossRef][Medline]
18 - Cruz Ramos, H., L. Boursier, I. Moszer, F. Kunst, A. Danchin, and P. Glaser. 1995. Anaerobic transcription activation in Bacillus subtilis: identification of distinct FNR-dependent and -independent regulatory mechanisms. EMBO J. 14:5984-5994.[Medline]
19 - Dahl, M. K., D. Schmiedel, and W. Hillen. 1995. Glucose and glucose-6-phosphate interaction with Xyl repressor proteins from Bacillus spp. may contribute to regulation of xylose utilization. J. Bacteriol. 177:5467-5472.[Abstract/Free Full Text]
20 - Daines, R. A., I. Pendrak, K. Sham, G. S. Van Aller, A. K. Kostandinidis, J. T. Lonsdale, C. A. Janson, X. Qiu, M. Brandt, S. S. Khandekar, C. Silverman, and M. S. Head. 2003. First X-ray cocrystal structure of a bacterial FabH condensing enzymes and a small molecule inhibitor achieved using rational design and homology modelling. J. Med. Chem. 46:5-8.[CrossRef][Medline]
21 - Daniel, R. A., E. J. Harry, V. L. Katis, R. G. Wake, and J. Errington. 1998. Characterization of the essential cell division gene ftsL (yIID) of Bacillus subtilis and its role in the assembly of the division apparatus. Mol. Microbiol. 29:593-604.[CrossRef][Medline]
22 - Darmon, E., D. Noone, A. Masson, S. Bron, O. P. Kuipers, K. M. Devine, and J. M. van Dijl. 2002. A novel class of heat and secretion stress-responsive genes is controlled by the autoregulated CssRS two-component system of Bacillus subtilis. J. Bacteriol. 184:5661-5671.[Abstract/Free Full Text]
23 - Debarbouille, M., R. Gardan, M. Arnaud, and G. Rapoport. 1999. Role of bkdR, a transcriptional activator of the sigL-dependent isoleucine and valine degradation pathway in Bacillus subtilis. J. Bacteriol. 181:2059-2066.[Abstract/Free Full Text]
24 - Delobbe, A., H. Chalumeau, and P. Gay. 1975. Existence of two alternative pathways for fructose and sorbitol metabolism in Bacillus subtilis Marburg. Eur. J. Biochem. 51:503-510.[Medline]
25 - Delorme, C., S. D. Ehrlich, and P. Renault. 1992. Histidine biosynthesis genes in Lactococcus lactis subsp. lactis. J. Bacteriol. 174:6571-6579.[Abstract/Free Full Text]
26 - Diaz, A., M. E. Pons, S. A. Lacks, and P. Lopez. 1992. Streptococcus pneumoniae DNA polymerase I lacks 3'-to-5' exonuclease activity: localization of the 5'-to-3' exonucleolytic domain. J. Bacteriol. 174:2014-2024.[Abstract/Free Full Text]
27 - Diaz, A., S. A. Lacks, and P. Lopez. 1992. The 5' to 3' exonuclease activity of DNA polymerase I is essential for Streptococcus pneumoniae. Mol. Microbiol. 6:3009-3019.[CrossRef][Medline]
28 - Droge, M. J., R. Bos, and W. J. Quax. 2001. Paralogous gene analysis reveals a highly enantioselective 1,2-O-isopropylideneglycerol caprylate esterase of Bacillus subtilis. Eur. J. Biochem. 268:3332-3338.[Medline]
29 - Eggert, T., G. van Pouderoyen, B. W. Dijkstra, and K. E. Jaeger. 2001. Lipolytic enzymes LipA and LipB from Bacillus subtilis differ in regulation of gene expression, biochemical properties, and three-dimensional structure. FEBS Lett. 502:89-92.[CrossRef][Medline]
30 - Ellwood, D. C., and D. W. Temoest. 1969. Control of teichoic acid and teichuronic acid biosyntheses in chemostat cultures of Bacillus subtilis var. Niger. Biochem. J. 111:1-5.
31 - Fang, G., E. Rocha, and A. Danchin. 2005. How essential are nonessential genes? Mol. Biol. Evol. 22:2147-2156.[Abstract/Free Full Text]
32 - Fillinger, S., S. Boschi-Muller, S. Azza, E. Dervyn, G. Branlant, and S. Aymerich. 2000. Two glyceraldehyde-3-phosphate dehydrogenases with opposite physiological roles in a nonphotosynthetic bacterium. J. Biol. Chem. 275:14031-14037.[Abstract/Free Full Text]
33 - Fisher, S. H., and L. V. Wray, Jr. 2002. Bacillus subtilis 168 contains two differentially regulated genes encoding L-asparaginase. J. Bacteriol. 184:2148-2154.[Abstract/Free Full Text]
34 - Fitzpatrick, T. B., N. Amrhein, and P. Macheroux. 2003. Characterization of YqjM, an Old Yellow Enzyme homolog from Bacillus subtilis involved in the oxidative stress response. J. Biol. Chem. 278:19891-19897.[Abstract/Free Full Text]
35 - Forsyth, R. A., R. J. Haselbeck, K. L. Ohlsen, R. T. Yamamoto, H. Xu, J. D. Trawick, D. Wall, L. Wang, V. Brown-Driver, J. M. Froelich, P. King, M. McCarthy, C. Malone, B. Misiner, D. Robbins, Z. Tan, Z. Y. Zhu, G. Carr, D. A. Mosca, C. Zamudio, J. G. Foulkes, and J. W. Zyskind. 2002. A genome-wide strategy for the identification of essential genes in Staphylococcus aureus. Mol. Microbiol. 43:1387-1400.[CrossRef][Medline]
36 - Gendron, N., H. Putzer, and M. Grunberg-Manago. 1994. Expression of both Bacillus subtilis threonyl-tRNA synthetase genes is autogenously regulated. J. Bacteriol. 176:486-494.[Abstract/Free Full Text]
37 - Gerth, U., J. Kirstein, J. Mostertz, T. Waldminghaus, M. Miethke, H. Kock, and M. Hecker. 2004. Fine-tuning in regulation of Clp protein content in Bacillus subtilis. J. Bacteriol. 186:179-191.[Abstract/Free Full Text]
38 - Gomez, A., D. Ramon, and P. Sanz. 1994. The Bacillus subtilis lipoprotein LplA causes cell lysis when expressed in Escherichia coli. Microbiology 140:1839-1845.[Abstract/Free Full Text]
39 - Greve, B., H. Saxild, and P. Nygaard. 2000. In W. Schumann, S. D. Ehrlich, and N. Ogasawara (ed.), Functional analysis of bacterial genes: a practical manual. John Wiley & Sons Ltd., Sussex, England.
40 - Haas, M., D. Beyer, R. Gahlmann, and C. Freiberg. 2001. YkrB is the main peptide deformylase in Bacillus subtilis, a eubacterium containing two functional peptide deformylases. Microbiology 147:1783-1791.[Abstract/Free Full Text]
41 - Hartig, E., A. Hartmann, M. Schatzle, A. M. Albertini, and D. Jahn. 2006. The Bacillus subtilis nrdEF genes, encoding a class Ib ribonucleotide reductase, are essential for aerobic and anaerobic growth. Appl. Environ. Microbiol. 72:5260-5265.[Abstract/Free Full Text]
42 - Heath, R. J., N. Su, C. K. Murphy, and C. O. Rock. 2000. The enoyl-[acyl-carrier-protein] reductases FabI and FabL from Bacillus subtilis. J. Biol. Chem. 275:40128-40133.[Abstract/Free Full Text]
43 - Heaton, M. P., R. B. Johnston, and T. L. Thompson. 1988. Controlled lysis of bacterial cells utilizing mutants with defective synthesis of D-alanine. Can. J. Microbiol. 34:256-261.[Medline]
44 - Holtmann, G., E. P. Bakker, N. Uozumi, and E. Bremer. 2003. KtrAB and KtrCD: two K+ uptake systems in Bacillus subtilis and their role in adaptation to hypertonicity. J. Bacteriol. 185:1289-1298.[Abstract/Free Full Text]
45 - Homuth, G., A. Rompf, W. Schumann, and D. Jahn. 1999. Transcriptional control of Bacillus subtilis hemN and hemZ. J. Bacteriol. 181:5922-5929.[Abstract/Free Full Text]
46 - Hullo, M. F., S. Auger, E. Dassa, A. Danchin, and I. Martin-Verstraete. 2004. The metNPQ operon of Bacillus subtilis encodes an ABC permease transporting methionine sulfoxide, D- and L-methionine. Res. Microbiol. 155:80-86.[Medline]
47 - Hunt, A., J. P. Rawlins, H. B. Thomaides, and J. Errington. 2006. Functional analysis of 11 putative essential genes in Bacillus subtilis. Microbiology 152:2895-2907.[Abstract/Free Full Text]
48 - Inaoka, T., Y. Matsumura, and T. Tsuchido. 1998. Molecular cloning and nucleotide sequence of the superoxide dismutase gene and characterization of its product from Bacillus subtilis. J. Bacteriol. 180:3697-3703.[Abstract/Free Full Text]
49 - Jacoby, G. A. 2005. Mechanisms of resistance to quinolones. Clin. Infect. Dis. 41:S120-S126.[CrossRef]
50 - Joyce, C. M., and N. D. Grindley. 1984. Method for determining whether a gene of Escherichia coli is essential: application to the polA gene. J. Bacteriol. 158:636-643.[Abstract/Free Full Text]
51 - Kamionka, A., and M. K. Dahl. 2001. Bacillus subtilis contains a cyclodextrin-binding protein which is part of a putative ABC-transporter. FEMS Microbiol. Lett. 204:55-60.[CrossRef][Medline]
52 - Kanamaru, K., S. Stephenson, and M. Perego. 2002. Overexpression of the PepF oligopeptidase inhibits sporulation initiation in Bacillus subtilis. J. Bacteriol. 184:43-50.[Abstract/Free Full Text]
53 - Karow, M. L., and P. J. Piggot. 1995. Construction of gusA transcriptional fusion vectors for Bacillus subtilis and their utilization for studies of spore formation. Gene 163:69-74.[CrossRef][Medline]
54 - Kawamura, F., and R. H. Doi. 1984. Construction of a Bacillus subtilis double mutant deficient in extracellular alkaline and neutral proteases. J. Bacteriol. 160:442-444.[Abstract/Free Full Text]
55 - Knuth, K., H. Niesalla, C. Hueck, and T. M. Fuchs. 2004. Large-scale identification of essential Salmonella genes by trapping lethal insertions. Mol. Microbiol. 51:1729-1744.[CrossRef][Medline]
56 - Kobayashi, K., S. D. Ehrlich, A. Albertini, G. Amati, K. K. Andersen, M. Arnaud, K. Asai, S. Ashikaga, S. Aymerich, P. Bessieres, F. Boland, S. C. Brignell, S. Bron, K. Bunai, J. Chapuis, L. C. Christiansen, A. Danchin, M. Debarbouille, E. Dervyn, E. Deuerling, K. Devine, S. K. Devine, O. Dreesen, J. Errington, S. Fillinger, S. J. Foster, Y. Fujita, A. Galizzi, R. Gardan, C. Eschevins, T. Fukushima, K. Haga, C. R. Harwood, M. Hecker, D. Hosoya, M. F. Hullo, H. Kakeshita, D. Karamata, Y. Kasahara, F. Kawamura, K. Koga, P. Koski, R. Kuwana, D. Imamura, M. Ishimaru, S. Ishikawa, I. Ishio, D. Le Coq, A. Masson, C. Mauel, R. Meima, R. P. Mellado, A. Moir, S. Moriya, E. Nagakawa, H. Nanamiya, S. Nakai, P. Nygaard, M. Ogura, T. Ohanan, M. O'Reilly, M. O'Rourke, Z. Pragai, M. Pooley, G. Rapoport, J. P. Rawlins, L. A. Rivas, C. Rivolta, A. Sadaie, Y. Sadaie, M. Sarvas, T. Sato, H. H. Saxild, E. Scanlan, W. Schumann, J. F. Seegers, J. Sekiguchi, A. Sekowska, S. J. Seror, M. Simon, P. Stragier, R. Studer, H. Takamatsu, T. Tanaka, M. Takeuchi, H. B. Thomaides, V. Vagner, J. M. van Dijl, K. Watabe, A. Wipat, H. Yamamoto, M. Yamamoto, Y. Yamamoto, K. Yamane, K. Yata, K. Yoshida, H. Yoshikawa, U. Zuber, and N. Ogasawara. 2003. Essential Bacillus subtilis genes. Proc. Natl. Acad. Sci. USA 100:4678-4683.[Abstract/Free Full Text]
57 - Krom, B. P., R. Aardema, and J. S. Lolkema. 2001. Bacillus subtilis YxkJ is a secondary transporter of the 2-hydroxycarboxylate transporter family that transports L-malate and citrate. J. Bacteriol. 183:5862-5869.[Abstract/Free Full Text]
58 - Kunin, V., and C. A. Ouzounis. 2003. The balance of driving forces during genome evolution in prokaryotes. Genome Res. 13:1589-1594.[Abstract/Free Full Text]
59 - Kunst, F., and G. Rapoport. 1995. Salt stress is an environmental signal affecting degradative enzyme synthesis in Bacillus subtilis. J. Bacteriol. 177:2403-2407.[Abstract/Free Full Text]
60 - Kunst, F., N. Ogasawara, I. Moszer, A. M. Albertini, G. Alloni, V. Azevedo, M. G. Bertero, P. Bessieres, A. Bolotin, S. Borchert, R. Borriss, L. Boursier, A. Brans, M. Braun, S. C. Brignell, S. Bron, S. Brouillet, C. V. Bruschi, B. Caldwell, V. Capuano, N. M. Carter, S.-K. Choi, J.-J. Codani, I. F. Connerton, N. J. Cummings, R. A. Daniel, F. Denizot, K. M. Devine, A. Düsterhöft, S. D. Ehrlich, P. T. Emmerson, K. D. Entian, J. Errington, C. Fabret, E. Ferrari, D. Foulger, C. Fritz, M. Fujita, Y. Fujita, S. Fuma, A. Galizzi, N. Galleron, S.-Y. Ghim, P. Glaser, A. Goffeau, E. J. Golightly, G. Grandi, G. Guiseppi, B. J. Guy, K. Haga, J. Haiech, C. R. Harwood, A. Hénaut, H. Hilbert, S. Holsappel, S. Hosono, M.-F. Hullo, M. Itaya, L. Jones, B. Joris, D. Karamata, Y. Kasahara, M. Klaerr-Blanchard, C. Klein, Y. Kobayashi, P. Koetter, G. Koningstein, S. Krogh, M. Kumano, K. Kurita, A. Lapidus, S. Lardinois, J. Lauber, V. Lazarevic, S.-M. Lee, A. Levine, H. Liu, S. Masuda, C. Mauël, C. Médigue, N. Medina, R. P. Mellado, M. Mizuno, D. Moestl, S. Nakai, M. Noback, D. Noone, M. O'Reilly, K. Ogawa, A. Ogiwara, B. Oudega, S.-H. Park, V. Parro, T. M. Pohl, D. Portetelle, S. Porwollik, A. M. Prescott, E. Presecan, P. Pujic, B. Purnelle, G. Rapoport, M. Rey, S. Reynolds, M. Rieger, C. Rivolta, E. Rocha, B. Roche, M. Rose, Y. Sadaie, T. Sato, E. Scanlan, S. Schleich, R. Schroeter, F. Scoffone, J. Sekiguchi, A. Sekowska, S. J. Seror, P. Serror, B.-S. Shin, B. Soldo, A. Sorokin, E. Tacconi, T. Takagi, H. Takahashi, K. Takemaru, M. Takeuchi, A. Tamakoshi, T. Tanaka, P. Terpstra, A. Tognoni, V. Tosato, S. Uchiyama, M. Vandenbol, F. Vannier, A. Vassarotti, A. Viari, R. Wambutt, E. Wedler, H. Wedler, T. Weitzenegger, P. Winters, A. Wipat, H. Yamamoto, K. Yamane, K. Yasumoto, K. Yata, K. Yoshida, H.-F. Yoshikawa, E. Zumstein, H. Yoshikawa, and A. Danchin. 1997. The complete genome sequence of the gram-positive bacterium Bacillus subtilis. Nature 390:249-256.[CrossRef][Medline]
61 - Lawhorn, B. G., S. Y. Gerdes, and T. P. Begley. 2004. A genetic screen for the identification of thiamin metabolic genes. J. Biol. Chem. 279:43555-43559.[Abstract/Free Full Text]
62 - Leatherbarrow, A. J., M. A. Yazdi, J. P. Curson, and A. Moir. 1998. The gerC locus of Bacillus subtilis, required for menaquinone biosynthesis, is concerned only indirectly with spore germination. Microbiology 144:2125-2130.[Abstract/Free Full Text]
63 - Lemma, E., J. Simon, H. Schagger, and A. Kroger. 1995. Properties of the menaquinol oxidase (Qox) and of qox deletion mutants of Bacillus subtilis. Arch. Microbiol. 163:432-438.[Medline]
64 - Lerner, C. G., B. T. Stephenson, and R. L. Switzer. 1987. Structure of the Bacillus subtilis pyrimidine biosynthetic (pyr) gene cluster. J. Bacteriol. 169:2202-2206.[Abstract/Free Full Text]
65 - Lerondel, G., T. Doan, N. Zamboni, U. Sauer, and S. Aymerich. 2006. YtsJ has the major physiological role of the four paralogous malic enzyme isoforms in Bacillus subtilis. J. Bacteriol. 188:4727-4736.[Abstract/Free Full Text]
66 - Lorca, G., B. Winnen, and M. H. Saier, Jr. 2003. Identification of the L-aspartate transporter in Bacillus subtilis. J. Bacteriol. 185:3218-3222.[Abstract/Free Full Text]
67 - Malle, D., T. Itoh, W. Hashimoto, K. Murata, S. Utsumi, and B. Mikami. 2006. Overexpression, purification and preliminary X-ray analysis of pullulanase from Bacillus subtilis strain 168. Acta Crystallogr. F 62:381-384.[CrossRef]
68 - Mansilla, M. C., and D. de Mendoza. 1997. L-Cysteine biosynthesis in Bacillus subtilis: identification, sequencing, and functional characterization of the gene coding for phosphoadenylylsulfate sulfotransferase. J. Bacteriol. 179:976-981.[Abstract/Free Full Text]
69 - Marino, M., T. Hoffmann, R. Schmid, H. Mobitz, and D. Jahn. 2000. Changes in protein synthesis during the adaptation of Bacillus subtilis to anaerobic growth conditions. Microbiology 146:97-105.[Abstract/Free Full Text]
70 - Marston, A. L., H. B. Thomaides, D. H. Edwards, M. E. Sharpe, and J. Errington. 1998. Polar localisation of the MinD protein of Bacillus subtilis and its role in selection of the mid-cell division site. Genes Dev. 12:3419-3430.[Abstract/Free Full Text]
71 - Martinez, S., P. Lopez, M. Espinosa, and S. A. Lacks. 1987. Complementation of Bacillus subtilis polA mutants by DNA polymerase I from Streptococcus pneumoniae. Mol. Gen. Genet. 210:203-210.[CrossRef][Medline]
72 - Mercuri, P. S., Y. Ishii, L. Ma, G. M. Rossolini, F. Luzzaro, G. Amicosante, N. Franceschini, J. M. Frere, and M. Galleni. 2002. Clonal diversity and metallo-beta-lactamase production in clinical isolates of Stenotrophomonas maltophilia. Microb. Drug Resist. 8:193-200.[CrossRef][Medline]
73 - Miao, L. X., J. W. Cao, R. J. Liu, Y. L. Wang, and Y. H. Zeng. 2002. Cloning and sequencing of the proBA gene from the selected mutant resistant to proline analogue from Bacillus subtilis. Yi Chuan Xue Bao 29:1111-1117. (In Chinese.)[Medline]
74 - Miethke, M., M. Hecker, and U. Gerth. 2006. Involvement of Bacillus subtilis ClpE in CtsR degradation and protein quality control. J. Bacteriol. 188:4610-4619.[Abstract/Free Full Text]
75 - Mijakovic, I., D. Petranovic, B. Macek, T. Cepo, M. Mann, J. Davies, P. R. Jensen, and D. Vujaklija. 2006. Bacterial single-stranded DNA-binding proteins are phosphorylated on tyrosine. Nucleic Acids Res. 34:1588-1596.[Abstract/Free Full Text]
76 - Nakano, M. M., F. Hajarizadeh, Y. Zhu, and P. Zuber. 2001. Loss-of-function mutations in yjbD result in ClpX- and ClpP-independent competence development of Bacillus subtilis. Mol. Microbiol. 42:383-394.[CrossRef][Medline]
77 - Nakano, M. M., T. Hoffmann, Y. Zhu, and D. Jahn. 1998. Nitrogen and oxygen regulation of Bacillus subtilis nasDEF encoding NADH-dependent nitrite reductase by TnrA and ResDE. J. Bacteriol. 180:5344-5350.[Abstract/Free Full Text]
78 - Noirot-Gros, M. F., P. Soultanas, D. B. Wigley, S. D. Ehrlich, P. Noirot, and M. A. Petit. 2002. The beta-propeller protein YxaL increases the processivity of the PcrA helicase. Mol. Genet. Genomics 267:391-400.[CrossRef][Medline]
79 - Ochiai, A., M. Yamasaki, T. Itoh, B. Mikami, W. Hashimoto, and K. Murata. 2006. Crystallization and preliminary X-ray analysis of the rhamnogalacturonan lyase YesW from Bacillus subtilis strain 168, a member of polysaccharide lyase family 11. Acta Crystallogr. F 62:438-440.[CrossRef]
80 - Ogawa, K., E. Akagawa, K. Yamane, Z. W. Sun, M. LaCelle, P. Zuber, and M. M. Nakano. 1995. The nasB operon and nasA gene are required for nitrate/nitrite assimilation in Bacillus subtilis. J. Bacteriol. 177:1409-1413.[Abstract/Free Full Text]
81 - Oussenko, I. A., R. Sanchez, and D. H. Bechhofer. 2004. Bacillus subtilis YhcR, a high-molecular-weight, nonspecific endonuclease with a unique domain structure. J. Bacteriol. 186:5376-5383.[Abstract/Free Full Text]
82 - Pang, A. S., S. Nathoo, and S. L. Wong. 1991. Cloning and characterization of a pair of novel genes that regulate production of extracellular enzymes in Bacillus subtilis. J. Bacteriol. 173:46-54.[Abstract/Free Full Text]
83 - Park, J. H., K. Burns, C. Kinsland, and T. P. Begley. 2004. Characterization of two kinases involved in thiamine pyrophosphate and pyridoxal phosphate biosynthesis in Bacillus subtilis: 4-amino-5-hydroxymethyl-2methylpyrimidine kinase and pyridoxal kinase. J. Bacteriol. 186:1571-1573.[Abstract/Free Full Text]
84 - Payne, M. S., and E. N. Jackson. 1991. Use of alkaline phosphatase fusions to study protein secretion in Bacillus subtilis. J. Bacteriol. 173:2278-2282.[Abstract/Free Full Text]
85 - Perego, M., C. F. Higgins, S. R. Pearce, M. P. Gallagher, and J. A. Hoch. 1991. The oligopeptide transport system of Bacillus subtilis plays a role in the initiation of sporulation. Mol. Microbiol. 5:173-185.[CrossRef][Medline]
86 - Putzer, H., N. Gendron, and M. Grunberg-Manago. 1992. Co-ordinate expression of the two threonyl-tRNA synthetase genes in Bacillus subtilis: control by transcriptional antitermination involving a conserved regulatory sequence. EMBO J. 11:3117-3127.[Medline]
87 - Qi, Y., Y. Kobayashi, and F. M. Hulett. 1997. The pst operon of Bacillus subtilis has a phosphate-regulated promoter and is involved in phosphate transport but not in regulation of the pho regulon. J. Bacteriol. 179:2534-2539.[Abstract/Free Full Text]
88 - Que, Q., and J. D. Helmann. 2000. Manganese homeostasis in Bacillus subtilis is regulated by MntR, a bifunctional regulator related to the diphtheria toxin repressor family of proteins. Mol. Microbiol. 35:1454-1468.[CrossRef][Medline]
89 - Reents, H., R. Munch, T. Dammeyer, D. Jahn, and E. Hartig. 2006. The Fnr regulon of Bacillus subtilis. J. Bacteriol. 188:1103-1112.[Abstract/Free Full Text]
90 - Reysset, G. 1981. New class of Bacillus subtilis glutamine-requiring mutants. J. Bacteriol. 148:653-658.[Abstract/Free Full Text]
91 - Richardson, D. J., B. C. Berks, D. A. Russell, S. Spiro, and C. J. Taylor. 2001. Functional, biochemical and genetic diversity of prokaryotic nitrate reductases. Cell. Mol. Life Sci. 58:165-178.[CrossRef][Medline]
92 - Rocha, E. P. C., A. Danchin, and A. Viari. 1999. Functional and evolutionary roles of long repeats in prokaryotes. Mol. Biol. Evol. 16:1219-1230.[Abstract]
93 - Rose, R. E. 1988. The nucleotide sequence of pACYC184. Nucleic Acids Res. 16:355.[Free Full Text]
94 - Ruiz, J. 2003. Mechanisms of resistance to quinolones: target alterations, decreased accumulation and DNA gyrase protection. J. Antimicrob. Chemother. 51:1109-1117.[Abstract/Free Full Text]
95 - Ruzal, S. M., and C. Sanchez-Rivas. 2003. Effect of glutamate synthase (GOGAT) activity on Bacillus subtilis spore properties. Curr. Microbiol. 47:208-213.[CrossRef][Medline]
96 - Sakai, A., K. Katayama, T. Katsuragi, and Y. Tani. 2001. Glycolaldehyde-forming route in Bacillus subtilis in relation to vitamin B6 biosynthesis. J. Biosci. Bioeng. 91:147-152.[CrossRef][Medline]
97 - Sakamoto, J. J., M. Sasaki, and T. Tsuchido. 2001. Purification and characterization of a Bacillus subtilis 168 nuclease, YokF, involved in chromosomal DNA degradation and cell death caused by thermal shock treatments. J. Biol. Chem. 276:47046-47051.[Abstract/Free Full Text]
98 - Salama, N. R., B. Shepherd, and S. Falkow. 2004. Global transposon mutagenesis and essential gene analysis of Helicobacter pylori. J. Bacteriol. 186:7926-7935.[Abstract/Free Full Text]
99 - Sambrook, J., and D. W. Russell. 2001. Molecular cloning: a laboratory manual, 3rd ed. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY.
100 - Samuelson, J. C., M. Chen, F. Jiang, I. Moller, M. Wiedmann, A. Kuhn, G. J. Phillips, and R. E. Dalbey. 2000. YidC mediates membrane protein insertion in bacteria. Nature 406:637-641.[CrossRef][Medline]
101 - Sanchez, H., D. Kidane, M. Castillo Cozar, P. L. Graumann, and J. C. Alonso. 2006. Recruitment of Bacillus subtilis RecN to DNA double-strand breaks in the absence of DNA end processing. J. Bacteriol. 188:353-360.[Abstract/Free Full Text]
102 - Sasaki, M., and Y. Kurusu. 2004. Analysis of spontaneous base substitutions generated in mutator strains of Bacillus subtilis. FEMS Microbiol. Lett. 234:37-42.[CrossRef][Medline]
103 - Saxild, H. H., J. H. Jacobsen, and P. Nygaard. 1994. Genetic and physiological characterization of a formate-dependent 5'-phosphoribosyl-1-glycinamide transformylase activity in Bacillus subtilis. Mol. Gen. Genet. 242:415-420.[Medline]
104 - Saxild, H. H., J. H. Jacobsen, and P. Nygaard. 1995. Functional analysis of the Bacillus subtilis purT gene encoding formate-dependent glycinamide ribonucleotide transformylase. Microbiology 141:2211-2218.[Abstract/Free Full Text]
105 - Saxild, H. H., K. Brunstedt, K. I. Nielsen, H. Jarmer, and P. Nygaard. 2001. Definition of the Bacillus subtilis PurR operator using genetic and bioinformatic tools and expansion of the PurR regulon with glyA, guaC, pbuG, xpt-pbuX, yqhZ-folD, and pbuO. J. Bacteriol. 183:6175-6183.[Abstract/Free Full Text]
106 - Schonert, S., S. Seitz, H. Krafft, E. A. Feuerbaum, I. Andernach, G. Witz, and M. K. Dahl. 2006. Maltose and maltodextrin utilization by Bacillus subtilis. J. Bacteriol. 188:3911-3922.[Abstract/Free Full Text]
107 - Schultz, A. C., P. Nygaard, and H. H. Saxild. 2001. Functional analysis of 14 genes that constitute the purine catabolic pathway in Bacillus subtilis and evidence for a novel regulon controlled by the PucR transcription activator. J. Bacteriol. 183:3293-3302.[Abstract/Free Full Text]
108 - Sciochetti, S. A., P. J. Piggot, D. J. Sherratt, and G. Blakely. 1999. The ripX locus of Bacillus subtilis encodes a site-specific recombinase involved in proper chromosome partitioning. J. Bacteriol. 181:6053-6062.[Abstract/Free Full Text]
109 - Sekowska, A., P. Bertin, and A. Danchin. 1998. Characterization of the polyamine synthesis pathway in Bacillus subtilis 168. Mol. Microbiol. 29:851-858.[CrossRef][Medline]
110 - Sievers, J., B. Raether, M. Perego, and J. Errington. 2002. Characterization of the parB-like yyaA gene of Bacillus subtilis. J. Bacteriol. 184:1102-1111.[Abstract/Free Full Text]
111 - Skarlatos, P., and M. K. Dahl. 1998. The glucose kinase of Bacillus subtilis. J. Bacteriol. 180:3222-3226.[Abstract]
112 - Tadesse, S., and P. L. Graumann. 2006. Differential and dynamic localization of topoisomerases in Bacillus subtilis. J. Bacteriol. 188:3002-3011.[Abstract/Free Full Text]
113 - Takamatsu, H., T. Kodama, and K. Watabe. 1999. Assembly of the CotSA coat protein into spores requires CotS in Bacillus subtilis. FEMS Microbiol. Lett. 174:201-206.[CrossRef][Medline]
114 - Tanaka, K., K. Kobayashi, and N. Ogasawara. 2003. The Bacillus subtilis YufLM two-component system regulates the expression of the malate transporters MaeN (YufR) and YflS, and is essential for utilization of malate in minimal medium. Microbiology 149:2317-2329.[Abstract/Free Full Text]
115 - Tanner, A., L. Bowater, S. A. Fairhurst, and S. Bornemann. 2001. Oxalate decarboxylase requires manganese and dioxygen for activity. Overexpression and characterization of Bacillus subtilis YvrK and YoaN. J. Biol. Chem. 276:43627-43634.[Abstract/Free Full Text]
116 - Tjalsma, H., S. Bron, and J. M. van Dijl. 2003. Complementary impact of paralogous Oxa1-like proteins of Bacillus subtilis on post-translocational stages in protein secretion. J. Biol. Chem. 278:15622-15632.[Abstract/Free Full Text]
117 - Tran, L., X. C. Wu, and S. L. Wong. 1991. Cloning and expression of a novel protease gene encoding an extracellular neutral protease from Bacillus subtilis. J. Bacteriol. 173:6364-6372.[Abstract/Free Full Text]
118 - Tsuge, K., T. Ano, M. Hirai, Y. Nakamura, and M. Shoda. 1999. The genes degQ, pps, and lpa-8 (sfp) are responsible for conversion of Bacillus subtilis 168 to plipastatin production. Antimicrob. Agents Chemother. 43:2183-2192.[Abstract/Free Full Text]
119 - Vagner, V., E. Dervyn, and S. D. Ehrlich. 1998. A vector for systematic gene inactivation in Bacillus subtilis. Microbiology 144:3097-3104.[Abstract/Free Full Text]
120 - van der Oost, J., C. von Wachenfeld, L. Hederstedt, and M. Saraste. 1991. Bacillus subtilis cytochrome oxidase mutants: biochemical analysis and genetic evidence for two aa3-type oxidases. Mol. Microbiol. 5:2063-2072.[CrossRef][Medline]
121 - Vincent, F., G. J. Davies, and J. A. Brannigan. 2005. Structure and kinetics of a monomeric glucosamine 6-phosphate deaminase: missing link of the NagB superfamily? J. Biol. Chem. 280:19649-19655.[Abstract/Free Full Text]
122 - von Blohn, C., B. Kempf, R. M. Kappes, and E. Bremer. 1997. Osmostress response in Bacillus subtilis: characterization of a proline uptake system (OpuE) regulated by high osmolarity and the alternative transcription factor sigma B. Mol. Microbiol. 25:175-187.[CrossRef][Medline]
123 - Wei, Y., A. A. Guffanti, M. Ito, and T. A. Krulwich. 2000. Bacillus subtilis YqkI is a novel malic/Na+-lactate antiporter that enhances growth on malate at low protonmotive force. J. Biol. Chem. 275:30287-30292.[Abstract/Free Full Text]
124 - Wei, Y., T. Havasy, D. C. McPherson, and D. L. Popham. 2003. Rod shape determination by the Bacillus subtilis class B penicillin-binding proteins encoded by pbpA and pbpH. J. Bacteriol. 185:4717-4726.[Abstract/Free Full Text]
125 - Winstedt, L., K.-I. Yoshida, Y. Fujita, and C. von Wachenfeldt. 1998. Cytochrome bd biosynthesis in Bacillus subtilis: characterization of the cydABCD operon. J. Bacteriol. 180:6571-6580.[Abstract/Free Full Text]
126 - Woods, S. A., J. S. Miles, R. E. Roberts, and J. R. Guest. 1986. Structural and functional relationships between fumarase and aspartase. Nucleotide sequences of the fumarase (fumC) and aspartase (aspA) genes of Escherichia coli K12. Biochem. J. 237:547-557.[Medline]
127 - Yoshida, K.-I., D. Aoyama, I. Ishio, T. Shibayama, and Y. Fujita. 1997. Organization and transcription of the myo-inositol operon, iol, of Bacillus subtilis. J. Bacteriol. 179:4591-4598.[Abstract/Free Full Text]
128 - Yoshida, K.-I., Y. Fujita, and S. D. Ehrlich. 2000. An operon for a putative ATP-binding cassette transport system involved in acetoin utilization of Bacillus subtilis. J. Bacteriol. 182:5454-5461.[Abstract/Free Full Text]
129 - Zhang, W., and D. H. Evans. 1995. DNA strand transfer catalyzed by the 5'-3' exonuclease domain of Escherichia coli DNA polymerase I. Nucleic Acids Res. 23:4620-4627.[Abstract/Free Full Text]
Journal of Bacteriology, January 2007, p. 591-602, Vol. 189, No. 2
0021-9193/07/$08.00+0 doi:10.1128/JB.01381-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
This article has been cited by other articles:
-
Commichau, F. M., Rothe, F. M., Herzberg, C., Wagner, E., Hellwig, D., Lehnik-Habrink, M., Hammer, E., Volker, U., Stulke, J.
(2009). Novel Activities of Glycolytic Enzymes in Bacillus subtilis: INTERACTIONS WITH ESSENTIAL PROTEINS INVOLVED IN mRNA PROCESSING. Mol. Cell. Proteomics
8: 1350-1360
[Abstract]
[Full Text]
-
Beek, A. T., Keijser, B. J. F., Boorsma, A., Zakrzewska, A., Orij, R., Smits, G. J., Brul, S.
(2008). Transcriptome Analysis of Sorbic Acid-Stressed Bacillus subtilis Reveals a Nutrient Limitation Response and Indicates Plasma Membrane Remodeling. J. Bacteriol.
190: 1751-1761
[Abstract]
[Full Text]