Altered Distribution of RNA Polymerase Lacking the Omega Subunit within the Prophages along the Escherichia coli K-12 Genome

The 91-amino-acid-residue small-subunit omega (the rpoZ gene product) of Escherichia coli RNA polymerase plays a structural role in the formation of RNA polymerase (RNAP) as a chaperone in folding the largest subunit (β′, of 1,407 residues in length), but except for binding of the stringent signal ppGpp, little is known of its role in the control of RNAP function. After analysis of genomewide distribution of wild-type and RpoZ-defective RNAP by the ChIP-chip method, we found alteration of the RpoZ-defective RNAP inside open reading frames, in particular, of the genes within prophages. For a set of the genes that exhibited altered occupancy of the RpoZ-defective RNAP, transcription was found to be altered as observed by qRT-PCR assay. All the observations here described indicate the involvement of RpoZ in recognition of some of the prophage genes. This study advances understanding of not only the regulatory role of omega subunit in the functions of RNAP but also the regulatory interplay between prophages and the host E. coli for adjustment of cellular physiology to a variety of environments in nature.

the regulatory interplay between prophages and the host E. coli for adjustment of cellular physiology to a variety of environments in nature. KEYWORDS Escherichia coli, RNA polymerase, omega subunit, prophage, transcription regulation R NA polymerase (RNAP) is the key enzyme of transcription. In Escherichia coli, the RNAP core enzyme is composed of four subunits, ␣ (RpoA), ␤ (RpoB), ␤= (RpoC), and (RpoZ), in the stoichiometry of ␣ 2 ␤␤= (reviewed in references 1 and 2. The core enzyme is assembled sequentially in the order 2␣ Ͼ ␣ 2 Ͼ ␣ 2 ␤ Ͼ ␣ 2 ␤␤= (3,4). In this assembly pathway, interacts with ␤=, forming ␤= intermediate (5), which then binds to the preformed ␣ 2 ␤ complex, leading to formation of the core enzyme. A model of chaperone function was proposed for 91-residue-long (RpoZ) in supporting the folding of the longest polypeptide, ␤= (RpoC), of 1,407 residues in size. RNAP purified from the rpoZ-defective mutant is associated with GroL, indicating the participation of GroL chaperone in place of in the process of RNAP formation (6). In agreement with this assembly mechanism of RNAP, the rpoZ gene is not an essential gene and the mutant lacking rpoZ is viable (7). In the core enzyme, RpoZ is present in nearstoichiometric amounts with respect to other subunits (8). The crystal structures of RNAP from Thermus aquaticus (9) and E. coli (10) confirm that RpoZ is one of the RNAP subunits.
Structure-function relationships have been extensively characterized for RpoA, RpoB, and RpoC (11,12), but except for the structural role in assembly of RNAP, the functional role played by RpoZ remains unsolved (13). When E. coli RNAP is associated with a dominant negative variant of RpoZ subunit, the resulting RNAP is defective in initiation of transcription, although preinitiated RNAP-RNA complex can elongate transcription (14). When RpoZ is tethered to DNA-binding protein, it is able to activate transcription through protein-protein contact between RNAP and the RpoZ segment (15). These findings indicate the direct molecular contact between RpoZ and other subunits of RNAP. Several lines of observation indicated the involvement of RpoZ in the functional control of RNAP (14,16). A functional link between the RpoZ subunit and the stringent response has been elucidated: first, in vitro transcription by purified RNAP without RpoZ subunit was found to be insensitive to ppGpp (17), and the RpoZ-less RNAP regains its sensitivity to (p)ppGpp upon the external addition of the RpoZ subunit (18) or the protein DksA, a collaborative player in the stringent response (19). Direct binding of ppGpp at a site near the rifampin binding site on the RpoB subunit was indicated by genetic and cross-linking studies (20)(21)(22), while the recent crystal structure of RNAP-ppGpp complex and mutational studies on E. coli RNAP showed that ppGpp binds at the interface of RpoC and RpoZ subunit (23)(24)(25). The level of RpoZ subunit influences DNA relaxation in E. coli (26), implying difference in DNA-binding properties between wild-type RNAP and RpoZ-defective RNAP. Microarray analysis of transcriptome in the absence of RpoZ indicated alteration in transcription of a set of genes, including the relA gene encoding ppGpp synthetase (7). These observations altogether indicate that the function of RNAP is controlled by the associated RpoZ subunit. Besides the ppGpp-binding site on the RpoZ subunit (ppGpp site 1), ppGpp binding was recently identified at an interface between RNAP and DksA (ppGpp site 2) (27). In the stringent control, a small transcription factor, DksA, participates in conjunction with ppGpp (19).
In this study, we made attempts to find the regulatory role of RpoZ subunit in the function of RNAP and transcription. As a shortcut approach to get insights into the overall functional role of RpoZ subunit in the transcription of the E. coli genome, we have performed, in this study, ChIP-chip (chromatin immunoprecipitation with microarray technology) analysis of the distribution of RNAP along the genome in the presence and absence of RpoZ subunit. Results of the ChIP-chip analysis indicated that the distribution of RNAP lacking RpoZ in growing E. coli K-12 cells is altered compared with that of wild-type RNAP mostly in the middle of the open reading frame (ORF) of a was observed for about 30 positions by setting a cutoff level of 5 (occupancy relative to that of wild-type RNAP) (Fig. 1). The genes that exhibited the high-level distribution of RpoZ-defective RNAP are located within some cryptic prophage regions. The test DNA obtained by ChIP-chip screening often binds more than two probes of the tiling array, forming a single peak, and thus, the total amount of RNAP binding of each single peak was calculated by combining the fluorescence intensity hybridized for each array probe within a single and the same peak. Among the total of 30 peaks of high-level distribution of RpoZ-defective RNAP detected at the cutoff level of 5.0, the highest-level distribution of RpoZ-defective RNAP was identified inside the ORF of the CP4-44 prophage flu gene encoding the Ag43 autotransporter (Fig. 1). The relative binding level of RpoZ-defective RNAP for this site was 36.6 compared with wild-type RNAP (see Table S2 for details). Almost half (47%) of the high-level distribution of the mutant RNAP belonged to the genes that are organized in specific prophages, including CP4-6, e14, Rac, CP4-44, CP4-57, and KpLE2 ( Fig. 1 shows the location of these prophages; for details, see Table S1). The increased distribution of RpoZ-defective RNAP was detected in multiple genes for e14 (5 genes), Rac (3 genes), CP4-6 (2 genes), and KpLE2 (2 genes) (Fig. 2). The set of genes that showed high-level distribution of RpoZ-lacking RNAP included ymfN (putative transcription factor [TF]), ymfK (putative repressor), ycfK (unknown protein), intE (predicted integrase), and ymfM (unknown protein) in e14 prophage; ydaY (pseudogene), stfR (predicted tail fiber protein), and ydaV (putative replication protein) in Rac prophage; ykfC (conserved protein) and ykfI (YkfI-YafW T-AT toxin) in CP4-6 prophage; and yjhH (predicted lyase) and yjhG (Dxylonate dehydrogenase) in KpLE2 prophage (Fig. 2 shows the locations of these genes in each prophage). Even though little is known about the functions of prophageencoded proteins, these cryptic phages contribute to cell physiology such as cell growth, resistance to antibiotics, stress responses, and biofilm formation (34), implying that at least some of the prophage-encoded proteins are expressed in E. coli supposedly under specific conditions. The genes showing the high-level distribution of RpoZ-lacking RNAP are not clustered but scattered along each prophage (Fig. 2). Including these genes, transcription organization of the prophage genes is not known yet (for details, see Discussion).

FIG 1
ChIP-chip analysis of wild-type (WT) and RpoZ-defective mutant RNAP. The ChIP-chip analysis was performed for E. coli K-12 wild-type BW25113 and its rpoZ-deleted mutant JW3624 under the standard procedure (78,79). The ratio of the binding level of rpoZ mutant to that of wild type was plotted along the E. coli K-12 genome. The peaks located within open reading frames are shown in orange, while the peaks located inside the spacer are shown in green. The list of genes with high-level distribution of RpoZ-lacking RNAP is given in Table 1, while the locations of these genes along prophages are shown in Fig. 2.
Among the top 30 genes that exhibited increased distribution of the RpoZ-deleted RNAP, about half are carried by the host E. coli K-12 genome (Table 1) and are located within a group of transporter genes such as the genes encoding LsrAC transporter for cell-cell communication signal AI-2 and ModABC transporter for molybdate. Signal transduction apparatus is also included, such as RcsBC two-component system (TCS) phosphotransferase and NtrBC TCS sensor. It is noteworthy that the cell surface receptor for phage N4 is also included in this group. Some of these genes could be considered the research targets for detailed analysis of the role of RpoZ in the functional control of RNAP.
The ChIP-chip pattern also indicated decreased distribution of the RpoZ-deleted RNAP inside the ORF of specific genes, again mostly within some prophages (Table S3). Marked reduction in the distribution of mutant RNAP was observed for prophage CP4-6 ( Fig. 3). The lowest difference, of 0.03, was observed inside the ORF of the CP4-6 prophage yagN and CPZ-55 prophage yffN genes (Table S3). Noteworthy is that the genes exhibiting the decreased distribution of the RpoZ-defective RNAP in the prophages are located within the so-called promoter islands containing promoter-like sequences (35). Although the level of constitutive promoters recognized by RpoD holoenzyme alone in the absence of additional supporting factors is low in the prophage regions (see below), the level of promoter islands is high in some prophage regions (35). One possibility is that these promoter-like sequences within the promoter islands might be differentially recognized by RNAP with and without RpoZ subunit (see below).
Transcription of the genes showing altered distribution of the RpoZ-defective RNAP. The altered distribution of the RpoZ-defective RNAP in some specific genes might be correlated with altered transcription, pausing, and/or attenuation-termination. To search for possible relationships between altered distribution of RpoZ-defective RNAP and altered transcription, in particular within prophages, attempts were then made to directly measure transcripts from the genes that showed altered distribution of RpoZdefective RNAP. From the set of genes that showed high-level distribution of RpoZdefective RNAP (Table 1), we selected 17 genes, i.e., 5 from the prophages and 12 from  Table 1). Locations of these genes within each prophage are shown by orange arrows along the genetic map of the respective prophage. High-level occupancy was detected in particular within e14 and Rac prophages.
Functional Role of RNA Polymerase Omega Subunit the host genome (Table 2). Using reverse transcription-quantitative PCR (RT-qPCR), transcripts of these genes were measured for these genes. As expected, the rpoZ mRNA was not detectable in the rpoZ-defective mutant strain. Among the total of 17 genes examined, the level of transcript was more than 4-fold higher for 7 genes (ydbA, ybdJ, yfiQ [CP4-57], paaE, ydaY [Rec], ybcH, and yhcD) and the transcript level was more than 2-fold higher for 6 genes (ymfN [e14], basS, yaiP, ybhA, ycfK, and ymfK [e14]) ( Table 3). This finding indicates that the high-level distribution of RpoZ-defective RNAP correlates with their high-level transcription.
In order to identify possible influence of the decreased binding of RpoZ-defective RNAP on transcription, we then analyzed the levels of transcripts for seven representative genes (Table 4) from the list of decreased binding of the mutant RNAP (Table 3). After RT-qPCR, transcripts were found to markedly decrease for yagM, yffL, yagN, yjfJ, yagL, and intF (Table 4). This finding indicates that the decreased distribution of RpoZ-defective RNAP correlates, to a certain extent, with their decreased transcription. Taken together, we predicted that the altered distribution of RpoZ-defective RNAP reflects the altered transcription of particular genes by the RpoZ-defective RNAP. Possible mechanisms of how only a specific set of genes is differently transcribed by the RpoZ-defective RNAP remain to be solved in future. In this aspect, the identification of contact partners, which collaborate with RNAP in different manners in the presence and absence of RpoZ, could be important.
Distribution of the constitutive promoters inside prophages. The prophage genes are generally silent and are not expressed under the steady state of cell growth, but some genes are induced under specific conditions, leading to influence on cell Putative protein a ChIP-chip assay was performed for E. coli K-12 wild type and its rpoZ-defective mutant. After tiling array analysis, the levels of RNAP binding were compared between two strains. The ratio between the binding of RpoZ-defective RNAP and that of wild-type RNAP was calculated for the entire probe along the genome. These values are listed in decreasing order in the Level column. When the gene is located inside prophages, the name of the prophage is given. OM, outer membrane; CoA, coenzyme A.
growth, resistance to antibiotics, stress responses, and biofilm formation (for instance, see reference 34. For instance, the minor sigma factor FecI is encoded by the KpLE2 prophage and regulates only the divergently transcribed fecABCDE operon in the same prophage. This small regulon is employed by E. coli for utilization of the ferric citrate transport system (36). The promoters recognized by E. coli RNAP should be fewer in prophages, because in the case of E. coli phages, only early genes are transcribed by the host RNAP, but afterward, the modified host RNAP by phage gene products (in the case of T-even and  Table 3). Locations of these genes within each prophage are shown by blue arrows along the genetic map of the respective prophage. Decreased-level occupancy was detected in particular within CP4-6 prophage. lambdoid phages) or the phage-encoded RNAP (in the case of T7-type phages) is responsible for the transcription of late genes (37). Except for a limited number of characterized genes such as the fec regulon in prophage KpLE2, however, almost nothing is known of the locations of promoters inside the prophages. We then analyzed  . When the target gene is located inside a prophage, the name of the prophage is given. b The level of mRNA was determined for E. coli K-12 mutants lacking the rpoZ gene by using the qRT-PCR method. The genes were selected from the list showing the decreased distribution of RpoZ-defective RNAP (Table 3 and Fig. 3). As references, two genes (a and b) were selected from the list of genes with increased distribution of RpoZ-defective RNAP ( Table 2 and Fig. 2).
the distribution of promoters in each prophage using the genomic SELEX (systematic evolution of ligands by exponential enrichment) screening system. RNAP holoenzyme was reconstituted from the sigma-free core enzyme and 4-fold excess of purified RpoD sigma. The constitutive promoters were predicted based on the binding sites of this reconstituted RpoD holoenzyme alone in the absence of other DNA-binding proteins (38). A maximum total of 669 constitutive promoters were identified on the E. coli K-12 genome (38, 39) (see Fig. S2A for the distribution of RpoD constitutive promoters). Only a small number of the constitutive promoters were identified on prophages CP4-6, e14, and Rac (Fig. 4), supposedly each contributing to the transcription of the argF gene encoding ornithine carbamoyltransferase (CP4-6), the lit gene encoding protease for cleavage of EF-Tu in collaboration with Gol protein (e14), and the trkG gene encoding K ϩ transporter (Rac), respectively. E. coli K-12 contains two genes, argF and argI, both encoding the ornithine carbamoyltransferase involved in the synthesis of L-citrulline from carbamoyl phosphate and L-ornithine along the pathway of arginine biosynthesis. Thus, the argF gene is considered to be integrated into the prophage CP4-6, together with its promoter, after duplication or transposition of the original argI gene (40). Likewise, the trkG gene encoding a K ϩ transporter is closely related to the trkH gene in Functional Role of RNA Polymerase Omega Subunit the genome of E. coli K-12 (41). Both TrkG and TrkH are active as low-affinity transporters of K ϩ and function in conjunction with TrkA, a membrane binding protein.
Thus, the trkG gene must have been inserted, together with its promoter, into the prophage Rac. The lit gene product blocks late gene expression of phage T4, leading to phage exclusion (42). This inhibitory activity depends on Gol protein of T4 gene 23, together functioning as peptidase for cleavage of EF-Tu (43). The putative promoter of the lit gene is, however, activated only after mutation (42). Taken together, it is unlikely that prophages contain many promoters recognized by E. coli K-12 RNAP. Distribution of H-NS silencer along the prophages. To avoid deleterious effects of the expression of foreign genes, including prophages, E. coli carries the gene silencing system, in which H-NS (44,45) and Rho (46)(47)(48) are known to participate as sentries. Even though these two silencing players have apparently similar functions with respect to the gene silencing, the mechanism is different between H-NS and Rho. H-NS is one of the nucleoid-associated proteins with both an architectural role in genome folding and a global regulatory role of transcription (49). The altered transcription of a set of prophage genes by RpoZ-defective RNAP might be due to differences in the interaction of RpoZ-defective RNAP with prophage genomeassociated H-NS silencer. We then analyzed the location of primary binding sites of H-NS along the genome of E. coli K-12 by using the Genomic SELEX screening system. A total of 987 binding sites were identified (38) (see Fig. S2B for the H-NS binding sites along the entire genome), of which a small number of H-NS binding sites were identified within specific prophages (Fig. 5). Noteworthy is that the strong binding sites of H-NS are almost absent for both e14 and Rac prophages, which showed the high-level distribution of RpoZ-defective RNAP. This finding implies that the binding of RpoZ-defective RNAP increased in the genes inside some prophages without strong binding sites for H-NS silencer.
Altered phenotypes of the RNAP lacking RpoZ detected by PM analysis. To see whether the differential occupancy of genes by RpoZ-less RNAP influences the phenotype of the cell, we studied the growth phenotype of the RpoZ-deleted strain with various carbon substrates and its tolerance to different environmental stresses like osmotic stress, pH, and antibiotic, using phenotype microarray (PM). PM technology allows us to find altered functions of genes by testing mutants for a large number of phenotypes simultaneously (50). An altered gene selection pattern of RpoZ-defective RNAP leads to the different phenotypes under various conditions tested. We then subjected both wild-type and rpoZ-defective E. coli strains to phenotype microarray (PM) assay.
The most prominent effect was seen on chemicals targeting the cell membrane, as the rpoZ-defective strain was remarkably sensitive toward sodium iodonitrotetrazolium (NT) violet (Fig. 6). NT violet is widely used for measurement of cellular redox activity or the respiratory activity in bacteria (50), indicating reduced transcription of the genes involved in respiratory metabolism for the rpoZ mutant. On the other hand, the growth of this mutant strain was significantly better than that of the wild type in the presence of gallic acid (PM19, A05-09) and phenethicillin (PM19, F01-04). Gallic acid, a type of phenolic acid, is known as a hydrated natural product of tannin and is commonly used in the pharmaceutical industry to produce a psychedelic alkaloid. Against bacteria, gallic acid is a toxic agent, a mutagen, and a modulator of amyloid formation (51,52). These actions of gallic acid are considered to be attributable to changes in membrane properties such as charge and hydrophobicity. Phenethicillin is a semisynthetic acidresistant penicillin, which is a methyl analog of phenoxymethyl penicillin. Penicillin group antibiotics bind to the penicillin-binding proteins (PBPs) and inhibit the crosslinking of peptidoglycan chains by PBPs (53), ultimately leading to weakening of the bacterial cell wall. Besides the antibiotics targeting the membrane, the rpoZ-defective mutant was sensitive to the ribosome-targeting oxytetracycline, a broad-spectrum tetracycline (PM20, F05-09) (Fig. 6). E. coli gains the resistance to tetracyclines through interference with its binding to ribosomes by ribosome protection proteins such as TetM (54). Resistance to tetracyclines is also mediated through decreased import or enhanced export of the drugs through membranes. In fact, the resistance by R factors is mediated through decrease of the intracellular level of tetracyclines (55).
The observed resistance of the rpoZ-defective mutant to these antibiotics detected by PM assay was then confirmed using individual liquid cultures in the presence of different concentrations of tetracycline (Fig. S3A) or penicillin (Fig. S3B). In the absence of drugs, the growth was retarded for the rpoZ mutant after late stationary phase. In the presence of these drugs, however, the growth retardation apparently disappeared. After prolonged culture, the growth of the rpoZ mutant was comparable to that of the wild type (Fig. S3) and continued longer than that of the wild type as observed by the PM assay (Fig. 6) and the tetrazolium reduction assay (Fig. S4). These observations altogether imply that the phenotypic resistance of the rpoZ mutant to these antibiotics could be attributable to the physiological modulation through some metabolic shifts due to altered transcription of as-yet-unspecified genes by the mutant RNAP. A close correlation has been established between the antibiotic sensitivity and the bacterial metabolism (56). Under starvation conditions after prolonged culture, the susceptibility to antibiotics decreases, leading to display of the phenotypic resistance. The altered distribution of RpoZ-defective RNAP in prophage regions and the altered sensitivity of the rpoZ-defective mutant to antibiotics both agree with the finding that the cryptic Functional Role of RNA Polymerase Omega Subunit prophages influence the sensitivity to a variety of environments, including the sensitivity to antibiotics (34).

DISCUSSION
E. coli RNAP functionally differentiates through two steps of protein-protein interaction: first, seven species of the sigma subunit, and second, more than 300 species of transcription factors (TFs) (57). TFs interact with one of the RNAP subunits for function (58)(59)(60). As a subunit of RNAP core enzyme, RpoZ might be involved in interaction with TFs. We then attempted to look into a broader perspective of the regulatory role played by RpoZ. Here, we identified the distribution of RNAP with and without RpoZ across the E. coli genome through ChIP-chip assay and found that the RNAP lacking RpoZ showed altered distribution within the prophage regions along the E. coli genome, supporting the prediction that RpoZ is involved in control of the regulatory function of RNAP.
A functional link between the RpoZ subunit and the stringent control has been suggested through its direct interaction with a nucleotide effector, ppGpp (17,18,24). Genetic and cross-linking studies of E. coli RNAP indicated that the direct binding of ppGpp to RpoB affects the catalytic activity of RNA polymerization (20,21), while the mutational studies of E. coli RNAP and the recent crystal structure of RNAP-ppGpp complex showed that ppGpp binds at the interface of RpoB, RpoC, and RpoZ subunit (23)(24)(25). Besides this ppGpp-binding site 1, ppGpp-binding site 2 was identified at an interface between RNAP and DksA (27). In the absence of RpoZ, marked alteration was not identified for the set of hitherto-identified genes under stringent control such as rRNA genes; the ppGpp-binding site 2 might play a major role in the stringent control.
Horizontal transfer of foreign genes is a major contributor to the evolution of prokaryotic genomes (61). The well-characterized model bacterium E. coli K-12 carries a set of 10 cryptic prophages and a short prophage segment (32, 61). Overall, a total of approximately FIG 6 PM analysis of wild-type E. coli K-12 and its rpoZ-defective mutant. PM assay was performed for the combination of E. coli wild-type BW25113 and its rpoZ-deleted mutant JW3624 according to the standard procedure (86,87). A typical pattern of the cell growth of the whole set of 20 PM plates is shown. The growth curve for wild type is shown in green, while that of the rpoZ mutant is shown in red. Differences in the growth rate were observed under several culture conditions as indicated. 304 genes, including predicted pseudogenes, exist in these prophages, and based on the similarity to phage genes, the functions have been predicted, but without experimental confirmation, for some genes in each prophage (see Table S1 in the supplemental material). One critical problem is whether these horizontally transferred genes are expressed and what the roles of these genes are in prophage survival inside the host E. coli K-12. Defective prophages Rac (62), e14 (63), DLP-12 (64), and Qin (65) are believed to carry some functional genes. Some of the prophage genes are beneficial to E. coli, including the genes encoding toxins and antibiotic resistance components for survival under various stressful natural conditions and for stable persistence in host animals. Sometimes, however, prophages kill E. coli through their induction. For instance, Rac repressor and integrase retain functional activity as conjugational transfer induces gene expression from the prophage and causes excision (66). In agreement with the expression of some functional genes from the Rac prophage, Rac is lethal to the host when its genes are expressed, resulting in the inhibition of cell division (67).
To avoid deleterious effects of prophages, E. coli carries the gene silencing system, in which H-NS (44,45,68) and Rho (46,47) are known to participate as sentries. H-NS is known as a nucleoid-associated global silencer to prevent transcription. The silencing function of H-NS is interfered with by global regulators such as LeuO (69). We predicted that the level of H-NS binding is related to the alteration of transcription in the presence and absence of RpoZ. On the other hand, Rho acts as a transcription terminator. The altered transcription of a set of prophage genes by RpoZ-defective RNAP might also be due to the difference in the interaction of RpoZ-defective RNAP with termination factor Rho. Direct interaction of RNAP with Rho or Rho-Nus protein complexes has been suggested (3,70). Involvement of Rho in the maintenance of some prophages has been suggested because the absence of the rho gene induces excision of defective prophages (47).
In the total of approximately 304 prophage genes, the genes for one RNAP sigma factor (FecI) and 14 TFs exist (Table S1), of which most are considered to control expression of the prophage genes. We have determined that the total number of constitutive promoters that are recognized by RNAP RpoD holoenzyme alone in the absence of supporting TFs is 492 to 669 (38). Along this line, the numbers of constitutive promoters recognized by the minor sigma factors were estimated to be 129 to 179 (RpoS), 101 to 142 (RpoH), 34 to 42 (RpoF), and 77 to 106 (RpoE) (71). In contrast, FecI is unique because it regulates only the divergently transcribed fecABCDE operon that encodes the ferric citrate transport system, suggesting that the regulatory target is still fixed within the KpLE2 prophage (36).
Among the total of 14 TFs, the regulatory target and function have been experimentally examined only for CP4-57-encoded AlpA (intA regulator) (72), DPL-12encoded AppY (acid phosphatase regulator) (73), and CP4-6-encoded XynR (xylonate catabolism regulator) (74). Regulatory targets of these TFs seem to be fixed to the genes inside prophages. For instance, AlpA of CP4-57 regulates only the intA gene for excision of the CP4-57 gene (72). The polysaccharide xylan, one representative renewable plant hemicellulose biopolymer, consists of D-xylose. The main pathway for utilization of D-xylose by E. coli K-12 depends on the xylFGH and xylAB operons. XylR is a TF that activates D-xylose import (xylFGH) and catabolism (xylAB) genes (75). The expression level of XylR is controlled by DicF small RNA (sRNA) encoded by Qin prophage (76). E. coli lacks the oxidation pathway of xylose, but once the oxidized product xylonate is provided by coexisting microorganisms in nature, E. coli is able to catabolize xylonate with the use of CP4-6-encoded YagEF enzymes. XynR (renamed YagI) on CP4-6 prophage regulates only the adjacent yagEF genes on the same CP4-6 prophage (74), indicating that XynR is a rare single-target TF and its regulatory target is still fixed on the CP4-6 prophage genes. Thus, E. coli gained the system for utilizing plant-derived xylose from the prophages. In utilization of plant-derived xylan, host TF (XylR) and prophage TFs (XynR and DicF) collaborate, thereby contributing to the stable maintenance of prophages CP4-6 and Qin. Likewise, DLP-12-encoded AppY is induced under anaerobic conditions and, in collaboration with ArcA, plays a role in induction of the expression of the hydrogenase 1 operon (hyaABCDEF) (77). ArcA is the response regulator of the quinone-dependent ArcAB two-component signal transduction system to respond to the change in respiratory conditions. Prophage-encoded AppY and host ArcA collaborate for antirepression against the repressor IscR, the iron-sulfur cluster [2Fe-2S] regulator. This is another example of the novel mode of prophage-host interaction, in which the prophage-encoded TF collaborates with the host TF so as to modulate the spectrum of regulation targets. Except for these three TFs, the regulatory function is not known for the other 13 prophage-encoded TFs.
We conclude that not only structurally but functionally as well, RpoZ is crucial for transcription as it guides the transcription machinery to express the genes necessary for viability under various environmental stresses. Our results indicate a future need to explore the role of RpoZ in alteration of the RNAP functions, including the interaction with as-yet-unidentified TFs, including H-NS and Rho.

MATERIALS AND METHODS
Bacterial strains. Escherichia coli strains BW25113 (parent strain) and JW3624 lacking the rpoZ gene (30) were used for all experiments. The strains were cultured in M9 medium (Difco) containing 0.4% glucose at 37°C. Cell growth was monitored by measuring the optical density at 600 nm (OD 600 ).
ChIP-chip analysis. The ChIP-chip assay was carried out as previously described (78,79) with a few modifications. BW25113 and JW3624 strains were grown in M9 glucose medium at 37°C to an optical density at 600 nm (OD 600 ) of 0.2 and then incubated in M9 medium containing formaldehyde (final concentration of 1%) for 30 min for cross-linking between proteins and genomic DNA. After termination of the cross-linking reaction by the addition of glycine, cells were collected, washed, and lysed by addition of lysozyme. The lysates were sonicated using a digital Sonifier (Branson) to fragment the genome. After centrifugation, the supernatant of the whole-cell extract was subjected to immunoprecipitation with the anti-RpoA antibody (Neoclone)-coated protein A-Dynal Dynabeads (Invitrogen). RNAP-conjugated DNAs in immunoprecipitation fractions were digested by pronase (Roche), and then the free DNA fragments were purified using a QIAquick PCR purification kit (Qiagen). Recovered DNA fragments were amplified by PCR using a pair of random primers. The amplified DNA fragments from the wild-type strain were labeled with Cy3, while another sample from the rpoZ mutant was labeled with Cy5. The labeled DNAs were mixed and hybridized to a 43,450-feature E. coli tiling array (Oxford Gene Technology) (80,81). After hybridization of samples to the DNA tiling array, the Cy5/Cy3 ratio was measured and the peaks of scanned patterns were plotted against the positions of DNA probes along the E. coli W3110 genome.
Genomic SELEX screening. The genomic SELEX screening was carried out as previously described (82,83). A mixture of DNA fragments of the E. coli K-12 W3110 genome was prepared after sonication of purified genome DNA and cloned into a multicopy plasmid, pBR322. In each SELEX screening, the DNA mixture was regenerated by PCR. For SELEX screening, 5 pmol of the mixture of DNA fragments and 10 pmol reconstituted RNAP RpoD holoenzyme or purified H-NS were mixed in a binding buffer (10 mM Tris-HCl, pH 7.8, at 4°C, 3 mM magnesium acetate, 150 mM NaCl, and 1.25 mg/ml bovine serum albumin). RNAP RpoD holoenzyme was prepared by mixing the purified sigma-free core enzyme and a 4-fold molar excess of the overexpressed and purified RpoD sigma (38) while H-NS was purified from overexpressed cells (84). The sequences of DNA fragments obtained by the genomic SELEX screening were identified by a SELEX-chip method as described previously (82,83). SELEX-chip data were submitted to the transcription factor profiling of Escherichia coli (TEC) at the National Institute of Genetics, Mishima, Japan (https://shigen.nig.ac.jp/ecoli/tec/top/; accession code RpoZ_ChIP).
RT-qPCR. Total RNA was prepared from E. coli cell as previously described (85). E. coli was grown in M9 glucose medium to an OD 600 of 0.2 at 37°C with shaking. Cells were harvested, and total RNA was prepared using hot phenol. The concentration of total RNA was determined by measuring the absorbance at 260 nm, and its purity was checked by agarose gel electrophoresis. Next, total RNAs were transcribed to cDNA with random primers using the Primer Script first-strand cDNA synthesis kit (TaKaRa Bio), and quantitative PCR (qPCR) was conducted using SYBR green PCR master mix (Applied Biosystems) as previously described (38,79). The primers used are described in Table S3 in the supplemental material. The cDNA templates were serially diluted to 2-fold and used in the qPCR assays. The levels of the 16S rRNA gene were used for normalization of data. The relative expression levels were quantified using the threshold cycle method presented by PE Applied Biosystems (PerkinElmer).
PM. Phenotype microarray (PM), a high-throughput technology for simultaneous testing of a large number of cellular phenotypes, was employed according to the manufacturer's instructions (BioLog) (86,87). In this study, PM was used for screening of the phenotypic differences between the wild-type strain and the rpoZ-deleted mutant strain. Growth difference was monitored by measuring the color intensity of oxidation of tetrazolium violet by NADH.