|
|
|
Journal of Bacteriology, May 2002, p . 2344-2351, Vol . 184, No . 9 Whole-Genome Analysis of Genes Regulated by the Bacillus subtilis Competence Transcription Factor ComKMitsuo Ogura,1 Hirotake Yamaguchi,2 Kazuo Kobayashi,3 Naotake Ogasawara,3 Yasutaro Fujita,2 and Teruo Tanaka1* School of Marine Science and Technology, Tokai University, 3-20-1 Orido, Shimizu, Shizuoka 424-8610,1 Faculty of Engineering, Fukuyama University, Fukuyama 729-0292,2 Nara Advanced Institute of Science and Technology, 8916-5 Takayama, Ikoma, Nara 630-01, Japan3 Received 5 November 2001/ Accepted 26 January 2002
MC medium is a minimal salt-glucose medium supplemented with potassium glutamate and Casamino Acids (21) . RNA preparation, cDNA synthesis, hybridization, and microarray analysis. Total RNAs were isolated from CU741 and CU741c cells that were grown in 100 ml of synthetic MC medium and were harvested 3 h after the end of the logarithmic growth phase . The procedures used for RNA preparation have been described previously (46) . The fluorescently labeled cDNA probes used for hybridization to DNA microarrays were prepared by a two-step procedure as described previously (32) . cDNAs from CU741and CU741c were labeled with Cy3 and Cy5, respectively . DNA microarrays were prepared by the method described previously (46) and were supplied by Takara Shuzo (Shiga, Japan) . A single microarray plate used in this study contained two grids, each containing 4,055 protein-encoding genes and 39 calf thymus DNA spots used as negative controls . The grids did not contain 45 genes for which PCR amplification was not successful . The hybridization and microarray analyses were performed as described previously (32, 46) . The mean values for the Cy3 and Cy5 fluorescence intensities for each gene on the two grids were calculated after subtraction of the background values obtained by determining the average values for the intensities of the 39 calf thymus DNA spots, and ratios of gene expression were calculated (Table 1) . transferrin receptor mRNA and the complementary human transferrin receptor primer were used as a positive, internal control (46) .
Transformation and ß-galactosidase assay. Transformation and ß-galactosidase activity analyses were carried out as described previously (29) . For the lacZ fusions showing low ß-galactosidase activities, the integrity of each fusion was checked by nucleotide sequence determination . To determine transformation efficiencies, two or more separate experiments were carried out with DNA from ODF200, and the averages are shown in Table 1 .
We carried out microarray analyses three times using RNAs isolated from independent cell cultures of CU741 and the comK-deficient mutant CU714c grown simultaneously in MC medium (see Materials and Methods) . The numbers of putative ComK target candidates showing ratios of more than three varied from experiment to experiment (159, 88, and 61 candidates for experiments 1, 2, and 3, respectively) . It is known that only a small fraction of a B . subtilis cell population becomes competent (12, 14), and the competent cell fraction is different in different experiments . In the three experiments, the intensity of the comK spot on the microarray was 3.2 to 4.1 times higher in experiment 1 than in experiments 2 and 3 despite the fact that in the three experiments the levels of expression for non-ComK-dependent genes were similar (data not shown) . Furthermore, the ratios of fluorescence intensity for known ComK target genes in ComK+ and ComK- cells were much higher in experiment 1 than in experiments 2 and 3 (data not shown) . These results suggest that the cells used in experiment 1 contained a higher level of ComK, and therefore, the variation in the number of target candidates may be attributed to different contents of competent cells in the cell population . It is also possible that mRNAs for some ComK-regulated genes were not extracted efficiently in experiments 2 and 3 due to their instability, although the same method was used for each preparation . On the basis of the reproducibility of the results, we selected 32 genes that were detected in all three experiments and 25 genes that were detected in both experiment 1 and either experiment 2 or experiment 3 . These genes are identified in Table 1 according to the experiments in which they were found . We found five genes that exhibited more than threefold expression in both experiment 2 and experiment 3 but not in experiment 1 . Three of these genes, degU, ywfI, and ycgM, were, however, false positives, since expression of the lacZ fusions of these genes was not dependent on ComK (data not shown) . The ß-galactosidase activities derived from the lacZ fusions of the remaining two genes, yndF and ytnI, were too low to determine the ComK dependence (data not shown), and therefore, these genes are not shown in Table 1 . We also include four genes in Table 1 that were detected only in experiment 1 but showed ComK dependence as determined with the lacZ fusions (see below) . The data in Table 1 indicate the levels of transcription of the genes observed in experiment 1 . Some of the genes detected were grouped in operons or putative operons by reference to the JAFAN database (http://bacillus.genome.ad.jp) . It should be noted that obviously not all of the ComK-regulated genes were detected by the microarray technique used here, because (i) the RNA was isolated under one set of experimental conditions, (ii) some ComK-regulated genes may not show a threefold change in expression, and (iii) some ComK target genes could have been among the 45 genes not present on the microarray (see Materials and Methods) . Known genes in ComK regulon. The known ComK-regulated genes include the late competence genes (the comC, comE, comF, and comG operon genes), nucA-nin, recA, addBA, and comK itself (15, 16, 43, 45) . All of these ComK-dependent genes except addBA were identified, indicating that gross transcriptional changes caused by comK deficiency could be detected under our assay conditions . One of the ComK-regulated genes, recA, showed the smallest ratio (12.1-fold) within this set of genes . The smaller ratio for recA than for the late competence genes may reflect expression of the recA gene in both the growth and competence development phases . It should be noted that none of the three microarray experiments revealed addBA as a target, although the operon is known to be regulated by ComK and carries the ComK box (15, 17) . The DNA for addBA was on the microarray plate, as shown by a relatively strong fluorescence spot on the microarray (data not shown) . The reason for the discrepancy, however, is not known . Examination of ComK dependence by ß-galactosidase assay. To confirm that the newly identified genes are indeed regulated by ComK, we examined expression of several such genes by using lacZ fusions constructed by and available from the Japan and EU consortia of B . subtilis functional genomics (see Materials and Methods) . The expression profiles of the genes are shown in Fig . 1 in the order of appearance in Table 1 . We found that expression of most of the lacZ fusions tested started after the end of the exponential growth phase, which is in agreement with the regulatory pattern of ComK . Except for the thdF, yhcF, and yhcE genes (Fig . 1A, M, and N, respectively), the levels of lacZ expression were greatly reduced during the stationary phase in the comK-deficient cells, demonstrating that 13 genes were dependent on ComK for expression . The ComK dependence of the thdF, yhcF, and yhcE genes became apparent only after 2 or 3 h of growth time relative to the end of vegetative growth . It is possible that expression of these genes is also controlled by some other regulator(s) in addition to ComK, which would result in only small differences in the levels of expression between the wild-type and comK-deficient cells . Although ComK dependence was clearly demonstrated for yvrP, yneB, sbcD, and ybdK (Fig . 1E, H, J, and P, respectively), these genes were detected only in experiment 1 . This could have been due to one or more of the possible reasons described above . Some of the fusion experiments did not show a difference in ß-galactosidase activities, apparently because the activities were too low to be evaluated (Table 1) . yjaZ was detected twice as a ComK target candidate in the microarray analysis (Table 1), but there was not a significant difference in the levels of expression of yjaZ-lacZ between the wild-type and comK-deficient cells (Table 1) . One interpretation is that expression of yjaZ is autoregulated by YjaZ, since the yjaZ gene is disrupted by insertion of lacZ in the host strain (41) . Another possibility is that the structure of the mRNA of yjaZ-lacZ is different from that of yjaZ and the mRNA of yjaZ-lacZ is subject to regulation different from ComK regulation . A definitive conclusion concerning this issue awaits further study .
lacZ fusion strains were not available for several of the genes that were tentatively shown to be regulated by ComK (Table 1) . We, therefore, examined the pattern of expression of the putative ComK target genes in the cells grown in minimal medium (9) by microarray-based time course analysis and found that expression of all the candidate genes was induced during the stationary phase, although the induction levels were relatively low for ytzD, yitT, yhzC, and yfnG (data not shown) . These results are consistent with the notion that the genes are regulated by ComK, but they do not prove it . Search for ComK boxes in upstream regions of candidate genes. It has been shown that ComK-regulated genes carry the ComK box (Fig . 2) in their control regions (17) . We, therefore, searched for a ComK box on the basis of both the typical nucleotide sequence motifs and the spacing between them and found 21 genes or operons that met the criteria (Table 1 and Fig . 2) . Among the 16 genes for which ComK dependence was observed in the lacZ fusion experiments (Fig . 1), possible ComK boxes were found for 10 genes (Fig . 2A), whereas no such sequences were identified for the remaining 4 genes (yvyF, yneA, sbcD, and yhcH) . It was previously shown that the yvyF gene is transcribed from the ComK-dependent comF promoter located upstream of this gene (24), and this may explain the ComK dependence of yvyF-lacZ expression (Fig . 1D) . The ysxA, yvrP, and smf loci carry a putative ComK box (Fig . 2A), but the sequence deviation from the typical ComK box is somewhat large . Therefore, it is possible that yneA, sbcD, yhcH, ysxA, yvrP, and smf are the secondary targets of ComK . A putative ComK box was found in the upstream region of yjaZ, for which a difference in expression in ComK+ and ComK- cells was detected by the microarray analysis but not by the lacZ fusion assay (see above) .
It should also be noted that the med-comZ operon, which has a canonical ComK box in the regulatory region (29, 30), was detected, although it was detected only once in experiment 1 (8.2- and 4.6-fold for med and comZ, respectively) . Newly identified ComK target genes . (i) Genes involved in transformation. The smf disruptant exhibited significantly reduced transformation efficiency (Table 1) . As determined by the BLAST search described in the JAFAN database, the putative Smf protein has sequence similarity with Haemophilus influenzae DprA, and disruption of the dprA gene results in reduced transformation efficiency (19) . It has been suggested that DprA might be part of the cellular apparatus for DNA transport and/or recombination in H . influenzae (19) . In Helicobacter pylori, the transformation efficiency is reduced by disruption of dprA (1) . Furthermore, H . influenzae DprA is homologous to Streptococcus pneumoniae CilB, and it has been shown that disruption of cilB results in competence deficiency (3) . From these observations, it may be reasonable to assume that smf is also involved in transformation in B . subtilis. Disruption of the rest of the newly identified genes did not result in significant decreases in transformation efficiency (Table 1) . These genes may play roles in some processes other than DNA transformation, although we cannot rule out the possibility that multiple disruption of such genes leads to transformation deficiency . (ii) Gene homologues induced in the competent state in other bacteria. B . subtilis ssb (Table 1) is homologous to S . pneumoniae cilA, which has been shown to be competence inducible (3, 36) . As shown in Fig . 1B, expression of yyaF-lacZ was comK dependent . Since yyaF, rpsF, and ssb were demonstrated to constitute an operon on the basis of Northern experiments (JAFAN database), ssb expression should also be comK dependent . Also, we note that a ysxA homologue of S . pneumoniae is induced during competence development (36) . Hence, the two homologues might play similar roles in transformation in the two organisms, although their functions are not known . (iii) Genes homologous to those acting on DNA structure or cleavage. Based on sequence similarities, the genes that are thought to affect DNA structure are shown in Table 1; these genes include ssb (encoding single-stranded DNA-binding protein), ywpH (encoding a single-stranded DNA-binding protein homologue), ysxA (encoding a DNA repair protein homologue), yneB (encoding a resolvase homologue), smf (encoding a DNA processing Smf homologue), and sbcD (encoding an exonuclease SbcD homologue, which degrades palindromic nucleotide structure) (5) . These proteins might be involved in DNA processing after the DNA translocation or recombination events . In fact, disruption of smf resulted in a decrease in transformation efficiency (see above) . ysxA and yneB disruptants, however, showed normal levels of transformation efficiency (95 and 102%, respectively), while ywpH and sbcD disruptants showed slight decreases in transformation (20 and 9%, respectively) . The ssb gene is essential (JAFAN database), which precluded us from testing the effect of gene disruption . Recently, NucA has been shown to be required for DNA cleavage during transformation, but nucA deficiency does not affect transformation efficiency (34, 43) . Taking these results into account, the slight reduction in transformation efficiency observed for ywpH and sbcD disruptants might indicate that these genes are involved in transformation . (iv) Genes involved in arginine biosynthesis, transportation, and formation of cell shape. The yhcHG genes, encoding a putative ABC transporter (35), were found to be regulated by ComK (Table 1 and Fig . 1L), which might suggest that an unknown substrate(s) is transported by YhcHG in competent cells . Surprisingly, expression of the argGH and argBD-carAB-argF operon genes was greatly reduced in comK cells (Table 1) . The comK-deficient cells, however, can grow in minimal medium (data not shown), indicating that these genes are expressed at least in noncompetent cells . Therefore, it is likely that expression of the arginine biosynthetic genes is enhanced in competent cells, although the biological significance of this is not known . The Maf protein is involved in cell shape determination, and overexpression of maf results in inhibition of septum formation, leading to filamentation of the cells and loss of viability (2) . The implications of the effect of ComK on maf are not known . (v) Regulatory factors. Expression of eight regulatory genes, including four putative genes (yyaN, sacXY, glcR, sinI, rapH, ybdK, and yhcF), was found to be ComK dependent, which suggests that secondary targets of ComK may be included in the list shown in Table 1 . Although sinR (37), the sacXY genes (6), and glcR (38) are known to be involved in certain stationary-phase events, including exoprotease production and sucrose metabolism, it is difficult to identify which genes listed in Table 1 are regulated by these genes . Further microarray analysis should provide valuable information on the regulatory network starting from ComK . It is interesting that ybdK encoding the sensor kinase of the putative ybdJK two-component system (20, 22) is regulated by ComK (Table 1 and Fig . 1P), which suggests that the YbdJ response regulator might affect competence-specific gene transcription . To examine this possibility further, we performed a microarray analysis of YbdJ-overproducing cells grown in MC medium and found that transcription of several of the comK regulon genes, including comK, was reduced in the range from 0.24- to 0.45-fold compared to transcription in the wild-type cells (data not shown) . These results suggest that there is a functional relationship between ComK and the YbdJK two-component system . We note that the rapH gene encoding a putative aspartyl-phosphate phosphatase for a phosphorylated response regulator (22, 23) is regulated by ComK (Table 1 and Fig . 1O) . It is tempting to speculate that the RapH phosphatase might have some functional relationship with the phosphorylated YbdJ response regulator, since rapH and ybdJ are in the same regulon . How these two genes are functionally correlated will be investigated in the future .
This work was supported by a Grant-in-Aid for Scientific Research on Priority Areas (C) "Genome Biology" and by Grants-in-Aid for Scientific Research (B) and (C) from the Ministry of Education, Science and Sports and Culture .
What Is Anthrax?,
What Is Fermentation?,
What Is Functional Genomics?,
What Is Biofilter?,
What Is Bioremediation?,
s,
Microbiology,
e,
Microbe,
i,
Bacteria,
n,
Microorganisms,
o,
Microorganism,
r,
Cell cultures,
o,
Erwinia,
e,
Thermophile,
e,
Enterococci,
n,
Acinetobacter,
e,
Clostridia,
a,
E coli O157,
c,
Antibiotics,
a,
Bacteroides,
c,
Minimum inhibiting concentration,
r,
Escherichia coli,
n,
Microorganisms,
n,
Water treatment,
n,
Fermentations,
c,
Staphylococcus,
r,
Prokaryotes,
i,
Microorganism,
c,
Salmonella typhimurium,
n,
Salmonella,
a,
Salmonella,
a,
Staphylococcus aureus
|
© 2005
Transgalactic Ltd (manufacturer of Bioscreen C software) |
Privacy Statement | P.O. Box
1393, 00101 Helsinki, Finland,
Last modified: May 25, 2005
| ||||||