INTRODUCTION

Top Abstract Introduction Materials & Methods Results Discussion References

Whooping cough is an acute respiratory disease caused by the small gram-negative bacterium Bordetella pertussis. B. pertussis expresses several factors that contribute to its ability to cause disease (15, 16, 33, 43, 45). Several of these factors, including filamentous hemagglutinin, pertactin, fimbriae, and the recently described tracheal colonization factor, contribute to the interaction between the bacterium and host cells. Other virulence factors are exotoxins which impair the function of immune cells and/or are capable of causing damage to host tissues. These include pertussis toxin, adenylate cyclase toxin, dermonecrotic toxin, and tracheal cytotoxin. Expression of these virulence factors, with the exception of tracheal cytotoxin, is activated at the level of transcription by a single locus, referred to as the bvg locus (originally designated the vir locus) (3, 38, 39, 42, 44). The bvg locus encodes a two-component regulatory system consisting of a sensor protein, BvgS, and a transcriptional activator, BvgA. Although the relevant signals for regulation of the bvg locus in vivo are unknown, activity of the bvg locus is repressed when cells are grown in the presence of MgSO₄ or nicotinic acid or when they are grown at a reduced temperature in vitro (26). This bvg-mediated change in the patterns of transcription in response to environmental signals is referred to as phenotypic modulation. Under nonmodulating conditions, autophosphorylation of BvgS on a conserved histidine residue is followed by two intramolecular phosphotransfer reactions and the transfer of the phosphate moiety to a conserved aspartate residue on BvgA (41). Upon phosphorylation by BvgS, BvgA binds to cis-acting sequences in the promoter regions of the bvg-activated genes and activates transcription (7, 8, 22).

In addition to the virulence factors which are activated by the bvg locus, a second class of genes, which is repressed by the bvg locus, has been described (24). The function(s) of the bvg-repressed genes is unknown. The fact that these genes are repressed by the locus responsible for regulation of known virulence genes suggests that they play a role in the pathogenesis of the bacterium, perhaps by contributing to late stages in the infectious cycle. Alternatively, their inappropriate expression may interfere with the bacterium's ability to cause disease. Five bvg-repressed genes have been identified to date (vrg6, vrg18, vrg24, vrg53, and vrg73), and the DNA sequences of the 5' ends of these genes have been determined (4, 5). Examination of the upstream regions of these genes revealed the presence of a conserved sequence element in four of the five genes (vrg6, vrg18, vrg24, and vrg53). The exception, vrg73, does not appear to contain this element. A 6-bp linker inserted into the conserved sequence element, as well as a single base pair change within the element in one of these genes, vrg6, eliminated responsiveness to modulation and resulted in constitutive expression (4, 5). A construct in which the sequences upstream of the initiating codon in vrg6 were replaced with a constitutive B. pertussis promoter was still regulated normally, demonstrating that the cis-acting sequences required for bvg-dependent repression of this locus are located downstream of the translation start site (5). The conserved sequence element in vrg6 was shown by Southwestern analysis to be bound by a 34-kDa protein which is present in nonmodulated cells but absent in modulated cells (5). These results, taken together, have led to a proposed model of bvg repression in which the bvg-repressed genes are regulated by a repressor protein, the expression or activity of which is activated by the bvg locus. A locus required for expression of repressor activity has been identified and shown to be located immediately downstream of the bvgS gene (30). This locus has been designated bvgR.

In this report, we describe the identification and characterization of the BvgR open reading frame. The presence of several errors in the published sequence of the region downstream of bvgS resulted in a failure to recognize previously the presence of this large open reading frame in the bvg locus (3). Analysis of the corrected sequence, presented herein, and the results of complementation analyses demonstrate that the bvgR gene is located immediately downstream of bvgS and that it is transcribed convergently to bvgAS. Activation of expression of the bvgR gene by BvgA was demonstrated to be at the level of transcription.

	MATERIALS AND METHODS

Top Abstract Introduction Materials & Methods Results Discussion References

Bacterial strains, plasmids, oligonucleotides, and media. The bacterial strains and plasmids used in this study are presented in Table 1. Escherichia coli strains were grown on L agar or in L broth supplemented with antibiotics when appropriate (31). B. pertussis strains were grown on Bordet-Gengou agar (Difco) containing 1% proteose peptone (Difco) and 15% defibrinated sheep blood. Concentrations of antibiotics, unless stated otherwise, were as follows: gentamicin sulfate, 10 µg/ml; kanamycin sulfate, 10 µg/ml; nalidixic acid, 50 µg/ml; rifampin, 50 µg/ml; streptomycin sulfate, 100 µg/ml. Plasmids were transformed into E. coli DH5 (Bethesda Research Laboratories, Bethesda, Md.).

Strain and plasmid construction. pTM193 derivatives bearing B. pertussis chromosomal fragments derived from the region downstream of the bvgA and bvgS genes were constructed as follows. Oligonucleotides MCS-1 (5'-CGCCG CGGATCCATATGAGCTC TAGATC TAG TAC TAG TCGACCATGG TACCAATTGAATTCGCTAGCATGC-3') and MCS-2 (5'-CGGCATGCTAGCGAAT TCAAT TGGTACCATGG TCGAC TAG TAC TAGATC TAGAGC TCATATGGATCCGCGG-3') were annealed to generate a linker bearing the restriction enzyme sites NarI, SacII, BamHI, NdeI, SacI, XbaI, BglII, ScaI, SpeI, SalI, HincII, NcoI, KpnI, MunI, EcoRI, NheI, SphI, and NarI. This linker was cloned into the NarI site of plasmid pBBRKAN to generate plasmid pTM193. The 9.0-kb BglII-BamHI, 6.0-kb SalI, 4.0-kb XhoI, 3.5-kb BglII-XhoI, or 1.3-kb SalI-XhoI restriction fragment derived from the region downstream of the bvg locus was cloned into the multiple cloning site of plasmid pTM193 to generate plasmids pTM193:BgB, pTM193:S, pTM193:X, pTM193:BgX, and pTM193:SX, respectively. Restriction fragments bearing BglII-proximal deletions of the 3.5-kb BglII-XhoI fragment were generated by PCR with primers X1 (5'-CTGACTCGAGAATGGCCTGCCGGTCCGCCACATCGAGCAG-3'), and either Bg2 (5'-GTCAAGATCTGGCTACGAATTGGCGCGCCGCATACGCGCC-3'), Bg3 (5'-GTCAAGATCTGGCCGAAGGTCTCGGACATGGCGCACAGGC-3'), Bg4 (5'-GTCAAGATCTTGGTACGTACGCTTGCGGCGCAGTCCGCCG-3'), Bg5 (5'-GTCAAGATCTCCTTACGGATAATTGGGCGGCATCTCGCGG-3'), or Bg6 (5'-GTCAAGATCTCATGGGCGGGGCCACATGGTCGCCCAGCAG-3'). Restriction fragments bearing XhoI-proximal deletions of the 3.5-kb BglII-XhoI fragment were generated by PCR with primers Bg1 (5'-GACTAGATCTCGAGAATGGCCTGCCGGTCCGCCACATC-3') and either X2 (5'-GTCACTCGAGCGCTGGTGGTTTCGATGCGGCCGCTGAAGG-3'), X3 (5'-GTCACTCGAGGCGACCCCGCATTGATCGGTATCGCCGACC-3'), X4 (5'-GTCACTCGAGCGGCCACATGCTGGTCACCGATCTGCTCAAC-3'), or X5 (5'-GTCACTCGAGCATGATCCACTGGACCAACAGCGCTCGCAG-3'). The products of these PCRs were double-stranded DNA fragments which contained terminal BglII and XhoI restriction enzyme sites. After digestion with BglII and XhoI, these fragments were inserted into plasmid pTM193, which had been digested with BglII and SalI. Although the entire sequences of the PCR products were not determined, at least four independent PCRs were performed to generate those products and the PCR products from each reaction gave the same result, indicating that differences in activity between the truncated DNA fragments generated by PCR and the full-length wild-type fragment are consequences of the deletion of required sequences at the termini rather than results of misincorporation of nucleotides during the extension reactions.

View larger version (19K): [in this window] [in a new window]   FIG. 1.   B. pertussis TM1311, TM1312, TM1315, and TM1316. The orf1-phoA transcriptional fusion present in strains TM1312 and TM1316 and the orf2-phoA transcriptional fusion present in strains TM1311 and TM1315 are shown. The putative BvgR (ORF1 and ORF2) coding sequences are represented by black boxes. The 3' terminus of the bvgS coding sequence is represented by gray boxes. The alkaline phosphatase coding sequence is represented by open boxes. The putative directions of transcription of the native genes and transcriptional fusions are indicated by arrows. Restriction enzyme recognition sequences are indicated as follows: E, EcoRI; B, BamHI; S, ScaI; X, XhoI; N, NcoI; A, ApaI.

Bacterial conjugations. Matings between E. coli and B. pertussis strains were performed by swabbing bacteria from fresh plate cultures of each strain onto a Bordet-Gengou agar plate supplemented with 10 mM MgCl₂ but without any antibiotics. After 3 h of incubation at 37°C, bacteria were swabbed onto Bordet-Gengou agar containing the appropriate antibiotics for selection of exconjugates, and incubation was continued at 37°C. Prior to mating, B. pertussis strains were grown for 3 days and E. coli strains were grown overnight at 37°C.

Quantitative alkaline phosphatase and -galactosidase assays. Bacteria to be assayed were recovered by sterile swabs into 3.5 ml of Tris-HCl, pH 8.0, and the absorbance at 600 nm was measured. For measurement of -galactosidase activity, 0.05 ml of cell suspension was added to 1 ml of Z buffer, cells were permeabilized by the addition of 30 µl of 0.1% sodium dodecyl sulfate and 30 µl of chloroform followed by vortexing, and the assay was completed as described by Miller (31). For measurement of alkaline phosphatase, 0.5 ml of cell suspension was added to 0.5 ml of Tris-HCl (pH 8.0), the cells were permeabilized as above, and the assay was completed as described by Brickman and Beckwith (9). Units in both cases were defined by the following equation: Units = [1,000 × A₄₂₀  (1.75 × A₅₅₀)]/(T × V × A₆₀₀) where T is the incubation time, in minutes, and V is the volume of permeabilized cells added to the assay mixture, in milliliters.

Sequence analysis. DNA fragments cloned in plasmid pBS KS⁺ were sequenced by the dideoxy sequencing method with the Sequenase 7-deaza dGTP DNA sequencing kit (United States Biochemical, Cleveland, Ohio). Computer analysis of DNA and protein sequences was performed with the GCG sequence analysis software package (Genetics Computer Group Inc., Madison, Wis.) and with the MacVector sequence analysis programs (International Biotechnologies Inc., New Haven, Conn.).

	RESULTS

Top Abstract Introduction Materials & Methods Results Discussion References

Complementation studies of the bvgR locus. We identified previously a locus in B. pertussis that is required for regulation of the bvg-repressed genes (30). This locus, designated bvgR, was mapped to a position immediately downstream of bvgAS. In order to functionally define the limits of the bvgR locus, we utilized the fact that bvgR-dependent regulation of a vrg6-phoA fusion can be restored in the bvgR-deficient strain TM1126 when the intact bvgR locus is provided in trans (30). A collection of restriction fragments containing different segments of DNA from the region immediately downstream of bvgAS were cloned into plasmid pTM193 (Fig. 2). These constructs were transferred by conjugation into strain TM1126, which contains a 12-bp in-frame insertion within the bvgR locus, a transcriptional fusion of the E. coli lacZ gene to the bvg-activated gene fhaB, and a transcriptional fusion of the E. coli phoA gene to the bvg-repressed gene vrg6. Quantitative enzyme assays for alkaline phosphatase and -galactosidase activities after growth in either the presence or absence of 50 mM MgSO₄ were performed in order to evaluate the ability of the exconjugates to complement the bvgR deficiency in strain TM1126. In all of the strains assayed, expression of the fha-lacZ fusion was induced between 13- and 38-fold, indicating that the functions of BvgA and BvgS were not altered in these constructs relative to the wild-type strain. As seen previously, the activity of the vrg6-phoA fusion in the wild-type strain shows an approximately 10-fold reduction in activity when the cells are grown in the absence of MgSO₄, while the activity of the vrg6-phoA fusion in strain TM1126 is reduced less than twofold under the same conditions (Fig. 2A). Strain TM1126 bearing plasmid pTM193 alone or pTM193:S, pTM193:X, or pTM193:SX failed to restore wild-type regulation of the vrg6-phoA fusion in strain TM1126 (Fig. 2A). Introduction of plasmid pTM193:BgB or pTM193:BgX restored wild-type regulation of the reporter fusion in strain TM1126. These results indicated that the sequences required for regulated expression of BvgR were located between the BglII site in bvgS and the XhoI site 1.5 kb downstream of bvgS.

View larger version (23K): [in this window] [in a new window]   FIG. 2.   Deletion and complementation analysis of the bvgR locus. The transcriptional activities of the vrg6-phoA transcriptional fusions in wild-type (WT) strain TM1081, strain TM1126, and strain TM1126 bearing the indicated pTM193 derivatives are shown. The bvg sequences inserted into the multiple cloning site of each pTM193 derivative are indicated by solid lines. (A) Effects of pTM193 alone and pTM193 bearing large restriction fragments derived from the bvg locus on expression of the vrg6 locus. (B) Effects of pTM193 bearing truncated versions of the bvg BglII-XhoI restriction fragment on expression of the vrg6 locus. The endpoints of deletions are indicated in parentheses. Nucleotide numbers correspond to those in the published sequence of bvgAS (3). Alkaline phosphatase activities are reported relative to strain TM1081 grown in the presence of MgSO4 (36.1 units). All reported values are averages of at least six independent assays. Restriction enzyme recognition sequences are indicated as follows: E, EcoRI; B, BamHI; S, SalI; X, XhoI; Bg, BglII.

Sequencing of the bvgR locus. A previous examination of the published sequence of the region immediately downstream of the bvgAS genes revealed three open reading frames which were disrupted by all of the mutations known to eliminate repressor activity in B. pertussis (30). Although the work of Beattie et al. (5) suggested that the bvg-dependent repressor was a protein of approximately 34 kDa, none of the open reading frames downstream of bvgS were predicted to encode a protein of that size. In addition, the boundaries of the bvgR locus, as determined by deletion and complementation analysis, were mapped to positions nearly 200 bp downstream from the ends of any of the three predicted open reading frames. In order to examine the possibility that the region defined by the deletion and complementation analysis contained an open reading frame large enough to encode a protein of approximately 34 kDa, we determined the DNA sequence between the SalI site at position 4272 in the published sequence and the XhoI site at position 5991. We identified 8 nucleotide residues between the 3' end of the bvgS gene and the XhoI site downstream of bvgS which were in error in the published sequence (see Fig. 3). Sequencing was performed on both strands of DNA for the entire region specified above. Nucleotides which were not in agreement with the published sequence were sequenced at least two times on each strand. Thus, a total of at least four independent sequencing reactions were performed, and the readings for each corrected nucleotide were in agreement in all of the reactions. Examination of the corrected sequence revealed the presence of two open reading frames predicted to encode proteins of approximately 30 kDa (Fig. 1). Open reading frame 1 (ORF1) extends from position 547 to position 1422 in the corrected sequence (Fig. 3) and is predicted to encode a 32-kDa protein. Open reading frame 2 (ORF2) extends from position 1347 to position 520 in the corrected sequence and is predicted to encode a 29-kDa protein. Transcription of ORF1 and bvgAS is predicted to occur convergently. Sequences which are good matches to the E. coli ribosomal binding site and the E. coli 10 promoter element are found upstream of ORF1 at positions 474 to 479 and positions 510 to 516, respectively. ORF2 and bvgAS are predicted to be transcribed in the same direction. Sequences which are good matches to the E. coli ribosomal binding site and 10 promoter element are also found upstream of ORF2, at positions 1379 to 1373 and 1417 to 1412, respectively.

View larger version (45K): [in this window] [in a new window]   FIG. 3.   Nucleotide sequence of bvgR. The nucleotide sequence of bvgR and the 3' end of bvgS and the corresponding amino acids are shown. A potential promoter element (10) and ribosomal binding site (RBS) are indicated by boxed-in sequences. Corrections of the published sequence are indicated as follows: nucleotides that are deleted in the published sequence are indicated by asterisks, and at those positions where a nucleotide was inserted or substituted in the published sequence, the nucleotide present in the published sequence is shown above the corrected sequence.

Activation of bvgR. In an effort to determine which of the two open reading frames encodes the bvg-dependent repressor, an ApaI fragment containing a promoterless E. coli phoA gene was inserted in both orientations into the ApaI site located at position 5853 of the published bvgAS sequence (Fig. 1). This allowed construction of transcriptional fusions of the phoA gene to ORF1 and ORF2. These constructs were transferred by conjugation into strain BP947, which contains a transcriptional fusion of the E. coli lacZ gene to fha, and strain BP907, which contains an in-frame deletion of the bvgA locus. The orientation of the phoA gene with respect to the bvg locus was determined by restriction enzyme digestion of plasmid DNA prior to introduction into strains BP947 and BP907 and also by PCR analysis of chromosomal DNA derived from strains TM1311, TM1312, TM1315, and TM1316 (data not shown). The activities of the bvg-phoA and fha-lacZ fusions in each strain were determined by analysis of the exconjugates with quantitative enzyme assays for alkaline phosphatase and -galactosidase after growth at 37°C. The results of this analysis are shown in Fig. 4. The fha-lacZ fusions in strains TM1315 and TM1316 demonstrated high levels of expression when the cells were grown in the absence of modulators. This activity was reduced approximately 25-fold upon growth of the cells in the presence of 50 mM MgSO₄ or 20 mM nicotinic acid and approximately 5-fold upon growth at 25°C (data not shown). Low-level constitutive alkaline phosphatase activity was detected under all of the conditions examined for the transcriptional fusion of the phoA gene to ORF2 (orf2-phoA) (strains TM1311 and TM1315). The transcriptional fusion of the phoA gene to ORF1 demonstrated high levels of expression in the wild-type background (strain TM1316) in the absence of modulators. This expression was reduced approximately 6.5-fold when the bacteria were grown in the presence of 50 mM MgSO₄. The orf1-phoATG fusion demonstrated low-level constitutive activity in strain TM1312, which bears an in-frame deletion of the bvgA locus. We investigated the regulation of expression of the orf1-phoA fusion further by determining the level of alkaline phosphatase activity expressed in strain TM1316 after growth in the presence of 20 mM nicotinic acid or after growth at 25°C, two conditions that have also been shown to modulate the activity of bvg-activated genes. The high level of alkaline phosphatase activity seen in strain TM1316 upon growth in the absence of modulators was reduced approximately 6.5-fold when the cells were grown in the presence of 20 mM nicotinic acid and approximately 4-fold when the cells were grown at 25°C.

View larger version (20K): [in this window] [in a new window]   FIG. 4.   Transcriptional activation of bvgR. Alkaline phosphatase activities of transcriptional fusions of the E. coli gene encoding alkaline phosphatase to ORF1 and ORF2 are shown after growth in the presence and absence of modulators. Black bars, no modulator; dotted bars, 50 mM MgSO4; diagonally striped bars, 20 mM nicotinic acid; white bar, growth at 25°C. Activities are reported relative to strain TM1316 grown in the absence of MgSO4 (153.26 units). All reported values are averages of at least four independent assays.

Analysis of the bvgR sequence. The predicted protein sequence of BvgR was determined and compared to sequences in the Swiss Protein, Protein Information Resource, and GenBank translated sequence databases. A search of the databases with the BLAST program showed that BvgR has significant homology to 18 proteins or predicted proteins found in a variety of eubacterial organisms, although the functions of these proteins are unknown (Fig. 5). The rtn locus of E. coli was identified by a transposon insertion that conferred an increased rate of adaptive mutation (19). The mechanism by which the rtn locus confers this phenotype, or if disruption of the rtn open reading frame is even responsible for the phenotype, has yet to be determined. The 22-kDa antigen of Borrelia burgdorferi has been identified, but its function is unknown (28). The functions of two of these predicted proteins have only been inferred, based upon their sequence homologies to characterized proteins in other organisms. These are the NifL homolog of Synechocystis sp. and the FixL homolog of E. coli (21) (D90789). It should be noted that the homology between these predicted proteins and BvgR does not extend to the characterized examples of these proteins (NifL of Klebsiella pneumoniae and FixL of Rhizobium meliloti). The functions of the remaining 14 of these predicted proteins are unknown, although their open reading frames are within or are tightly linked to operons with defined functions. Sequence alignment by the LFASTA program initially showed that a region of 55 amino acid residues was highly conserved among 16 of the 19 proteins (Fig. 5A). Domain I of BvgR is most closely related to domain I of YAHA_ECOLI, with 31% identity and 47% similarity over the 55-amino-acid region. Alignment of the 15 protein sequences identified as having homology to BvgR within domain I revealed the presence of a second region of high sequence homology (Fig. 5B). Although this family of proteins was originally defined by homology to domain I of BvgR, BvgR itself shows very poor conservation of domain II. Within this family, we have observed that in addition to the high degree of sequence conservation within the two domains, the spacing between the two domains and the distance between domain II and the carboxy-terminal ends of the predicted proteins are also conserved. (Fig. 5C). It should be noted that we have included in this analysis only those proteins and predicted proteins for which some information is available. An additional 12 open reading frames in the Synechocystis genome, 9 open reading frames in the E. coli genome, and 2 open reading frames in the Mycobacterium tuberculosis genome that contain these two conserved sequence elements were identified, bringing the total number of potential members of this protein family to 42. 

View larger version (73K): [in this window] [in a new window]   FIG. 5.   Alignment of the predicted amino acid sequence of BvgR with a new family of proteins. BVGR_BPER, predicted protein BvgR of B. pertussis; YAHA_ECOLI, hypothetical 40.7-kDa protein encoded by the betT 3' region of E. coli (27); NIFL_SLL, NifL homolog of Synechocystis sp. (21); STYRESORF, open reading frame 4 of the Salmonella typhimurium resolvase operon (25); FIMK_KPNEU, FimK protein of K. pneumoniae (17); RTN_ECOLI, Rtn protein of E. coli (19, 20); ISTN501ORF, unidentified reading frame 2 of the ISTN501-encoded mercuric ion resistance operon of Pseudomonas aeruginosa (10-12, 14, 32); ECOMERORF, unidentified reading frame 2 of the E. coli mercuric ion resistance operon (R100) (1, 11); SFLMERORF, unidentified reading frame 2 of the Shigella flexneri mercuric ion resistance operon (Tn501) (11); SMAMERORF, unidentified reading frame 2 of the Serratia marcesens mercuric ion resistance operon (pDU1358) (18); FIXL_ECOLI, FixL homolog of E. coli (D90789); PHAX_RMEL, unidentified reading frame X of the pha operon of R. meliloti (X93358); YNTC_AZOCA, hypothetical 80.5-kDa protein encoded by the ntrC 5' region of Azorhizobium caulinodans (34); YJCC_ECOLI, hypothetical 60.8-kDa protein encoded by the ssb-soxS intergenic region (6); YHJK_ECOLI, hypothetical 73.1-kDa protein encoded by the dctA-dppF intergenic region (36); YFEA_ECOLI, hypothetical protein encoded by the gltX 5' region of E. coli (13, 46). FCCORF_WOSU, unidentified open reading frame in the fcc 5' region of Wolinella succinogenes (Y10581); DEFORF_CALO, unidentified open reading frame in the def 3' region of Calothrix sp. (29); 22AG_BBUG, 22-kDa antigen of B. burgdorferi (28). Black boxes indicate positions conserved in at least 51% of the aligned sequences.

	DISCUSSION

Top Abstract Introduction Materials & Methods Results Discussion References

The locus required for regulation of the bvg-repressed genes in B. pertussis is located immediately downstream of bvgAS (30). This locus was designated bvgR, for Bordetella virulence gene repression. Examination of the published sequence of the bvgAS locus revealed the presence of three open reading frames in the region immediately downstream of bvgAS that are predicted to be affected by all of the mutations that have been demonstrated to abolish repressor activity (30). The largest of these open reading frames is predicted to encode a protein of approximately 25 kDa. Previous results published by Beattie et al. have suggested that the repressor protein that binds to the conserved sequence element in the vrg6 gene is a protein of approximately 34 kDa (5). It is possible that rather than repressing the bvg-repressed genes directly, the product of the bvgR locus may be required for the expression or activity of the repressor protein. If, however, bvgR does encode the repressor protein, then either the repressor is smaller than 34 kDa in size or the bvgR open reading frame is larger than those open reading frames present in the published sequence. In order to begin to address these possibilities, we sought to functionally define the upstream and downstream limits of the bvgR locus. Various plasmid derivatives of plasmid pTM193 bearing sequences derived from the bvgR locus were introduced into strain TM1126, which bears an insertional mutation in the bvgR locus which abrogates bvgR function. The ability of these derivatives to provide bvgR function in trans was determined. This analysis allowed localization of the cis-acting sequences required for regulated expression of bvgR to between positions 4447 and 5639 in the published bvgAS sequence. It should be noted, however, that although this system provides an effective method for determining the boundaries of the bvgR locus, it provides only an indirect measure of the level of bvgR expression. It is possible, for example, that promoter elements contributing to bvgR expression reside upstream of position 5639, since a requirement for these sequences would not be detected by the assay utilized in this study if less than 100% of bvgR expression was sufficient to fully activate expression at the vrg6 locus. Future studies will focus on the identification and characterization of all cis-acting sequences that contribute to bvgR expression.

Having determined the functional boundaries of the bvgR locus, we next undertook the sequencing of the region between the end of the bvgS open reading frame and the XhoI site 1.5 kb downstream of this site. This analysis resulted in the identification of 8 nucleotides which were in error in the published sequence (Fig. 3). Examination of the corrected sequence revealed the presence of two large open reading frames downstream of bvgS within the boundaries determined by deletion and complementation analysis to be sufficient to provide bvgR activity in trans (Fig. 3). One open reading frame (ORF2) was located almost in its entirety within the other open reading frame (ORF1), although the two open reading frames are encoded on opposite strands of the DNA. Both open reading frames would be expected to be disrupted by all of the mutations known to abolish repressor activity in B. pertussis, and both open reading frames are large enough to encode the protein identified as the repressor protein in the studies of Beattie et al. (5). The presence of two large open reading frames immediately downstream of bvgAS, in opposite orientation to each other, suggests two models for the regulated transcription of bvgR. If bvgR is encoded by ORF2, it lies immediately downstream of bvgS oriented in the same direction as bvgA and bvgS, which would allow all three proteins to be expressed from a single transcript under the regulation of the bvgAS promoter. If, however, bvgR is encoded by ORF1, its transcription is presumably driven from its own bvgA-activated promoter. Analysis of the transcription in the region downstream of bvgS revealed that there is only a low constitutive level of transcription through ORF2. In contrast, there was a high level of transcriptional activity through ORF1 that was reduced upon growth under conditions known to modulate the expression of bvg-activated genes; this transcription was dependent upon an intact bvgA locus. These results demonstrate that ORF1, but not ORF2, encodes a bvg-activated gene and support the conclusion that bvgR is encoded by ORF1. Although it remains to be demonstrated that bvgR encodes the actual repressor of the bvg-repressed genes, the fact that the product of the bvgR locus is predicted to be a protein of 32 kDa is consistent with that conclusion. It is interesting that the bvg-activated transcription of bvgA and bvgS does not extend through the bvgR locus. This is perhaps a little surprising, given that the end of bvgR lies only 43 nucleotides downstream from bvgS. Examination of the sequence between bvgS and bvgR reveals the presence of several large potential stem-loop structures. The presence of one or more of these stem-loop structures may lead to transcriptional termination in this region. Perhaps this is required in order to ensure that overlapping transcription of bvgAS and bvgR does not interfere with expression of the products of these genes.

The predicted BvgR protein encoded by ORF1 contains a domain that is strongly conserved among a large number of predicted proteins from a variety of eubacterial species (Fig. 5). This family of predicted proteins has been identified based on its homology to domain I of BvgR. Most of the members of this family are characterized by two domains of strong sequence homology which are located at the carboxy-terminal ends of the proteins and which reside approximately 130 residues apart. Several members of this family of proteins have either domain I or domain II but not both. The assignment of BvgR to this family of proteins is based on its high degree of conservation within domain I. In this region, BvgR is most closely related to YAHA_ECOLI, with 30% identity and 46% similarity over the 55 residues. BvgR shows weak but discernible conservation in domain II, suggesting that BvgR either has lost the requirement for the function of Domain II or has adapted this domain for an alternative function. The functional role(s), if any, of domain I and domain II is only speculative at this time, since no known function has been demonstrated for any of the members of this protein family. It is interesting to note that members of this family of proteins are encoded by open reading frames which either lie within or are very closely linked to functional operons. The identification and characterization of BvgR allows, for the first time, the assignment of an activity to one of the members of this family of proteins. It will be interesting to evaluate whether the members of this family are involved in regulation of expression of the operons to which they are so closely linked.

In this study, we have identified the bvgR open reading frame. The close correlation between the predicted size of the bvgR product and the protein shown by Beattie et al. (5) to bind the putative repressor binding site in the vrg6 gene suggests that BvgR is the actual repressor protein that binds to the conserved sequence element found within the coding sequences of the bvg-repressed genes. Our results demonstrate that the expression of bvgR is activated at the level of transcription by the products of the bvgAS genes. We propose that this activation is the result of binding of phosphorylated BvgA to the bvgR promoter. Ongoing and future work is and will be focused on characterizing in detail the activation of bvgR expression and determining the mechanism by which BvgR represses the expression of its target genes.