Transposable Prophage Mu Is Organized as a Stable Chromosomal Domain of E. coli
2013

R P Saha
Z Lou
L Meng
R M Harshey
Transposable Prophage Mu Is Organized as a Stable Chromosomal Domain of E. coli

PLoS Genet
91110039022013Received July 19, 2013; Accepted September 6, 2013;Editor: William F. Burkholder, Institute of Molecular and Cell Biology (IMCB), A*STAR, Singapore Funding: This work was supported by National Institutes of Health grant GM 33247 and in part by the Robert Welch Foundation grant F-1351. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * E-mail: firstname.lastname@example.com . These authors contributed equally to this work.
The E. coli chromosome is compacted by segregation into 400-500 supercoiled domains by both active and passive mechanisms, for example, transcription and DNA-protein association. We find that prophage Mu is organized as a stable domain bounded by the proximal location of Mu termini L and R, which are 37 kbp apart on the Mu genome. Formation/ maintenance of the Mu 'domain' configuration, reported by Cre-loxP recombination and 3C (chromosome conformation capture), is dependent on a strong gyrase site (SGS) at the center of Mu, the Mu L end and MuB protein, and the E. coli nucleoid proteins IHF, Fis and HU. The Mu domain was observed at two different chromosomal locations tested. By contrast, prophage l does not form an independent domain. The establishment/maintenance of the Mu domain was promoted by low-level transcription from two phage promoters, one of which was domain dependent. We propose that the domain confers transposition readiness to Mu by fostering topological requirements of the reaction and the proximity of Mu ends. The potential benefits to the host cell from a subset of proteins expressed by the prophage may in turn help its long-term stability.

Introduction

Bacterial chromatin is spatially organized and condensed ,1000 fold to fit inside a bacterial cell [1]. Referred to as a nucleoid, E. coli chromatin is organized in a series of negatively supercoiled loops [2,3], segregated by dynamic domain barriers (defined as entities that prevent the free diffusion of supercoils) and compacted by several nucleoid-associated proteins (NAPs) including HU, IHF, Fis and H-NS [4,5]. The chromosome is not randomly condensed, but rather has a ring organization with four structured macro-domains and two less-structured regions; interactions between these regions are highly restricted as determined by cytological and genetic analyses [6]. Macro-domains are thought to orchestrate chromosome movements during the cell cycle [7].

Supercoiling not only plays a vital role in compacting the chromosome, but a proper degree of supercoiling is crucial for all DNA-related processes [1]. Segregation of supercoils into topological domains protects these processes by preventing DNA breaks from relaxing the entire chromosome [2,8]. The level of DNA superhelicity is tightly controlled by the combined activities of topoisomerases and histone-like proteins; the latter not only constrain negative supercoils and generate diffusion barriers for the formation of topological domains, but are also global regulators of gene transcription [4,[9][10][11]. The critical importance of DNA supercoiling in interconnecting chromosome structure and global gene transcription was reinforced in recent evolution experiments where supercoiling was observed to be under strong selection in E. coli populations [12].

Transposable phage Mu is a temperate phage that integrates into essentially random locations on the E. coli chromosome [13][14][15]. Transposition from mini-Mu plasmids in vitro requires DNA supercoiling for formation of a high-order transpososome within which the two Mu ends are interwound and synapsed [13,16]. Supercoiling is inferred to be similarly important in vivo [17], where Mu end pairing during replicative transposition additionally requires a centrally located gyrase binding site SGS [18]. This site is the strongest such site studied, and is found only in Mu-like prophages [19]. Highly processive supercoiling by gyrase bound at the SGS has been proposed to propagate a supercoiled loop, with SGS at the apex, assisting transposase-mediated synapsis of Mu ends at the base [20].

In vitro studies found that the Mu transposase mediates an ordered interaction of three cis-acting sites -Mu L and R ends and an enhancer element E -which traps five supercoils before pairing the L and R ends in their reactive configuration [21]. The original goal of the present study was to test if the three sites interact in a similar order in vivo upon initiation of transposition. The in vitro studies employed Cre recombinase-mediated exchange at two strategically placed loxP sites to determine the topology of interactions between a given Mu site and the other two [16]. Using this same strategy in vivo, we found to our surprise that the ends were already paired in a prophage, and that the transposase was not essential for their pairing. We have investigated the basis of this pairing using both the Cre-lox system as well as the 3C crosslinking system [22,23]. We show that Mu SGS is important for Mu end pairing, that other Mu cis and trans factors and several host NAPs contribute as well, and that the MuB protein, expressed at a low level in the prophage, likely provides a NAP-like function. We discuss the implications of this work for the maintenance of large selfish DNA elements on a bacterial genome.

Results

The two ends of the 37 kbp Mu prophage genome behave as if they are paired

The Cre-loxP site-specific recombination system has simple requirements, needing only two loxP sites; neither additional cofactors nor special DNA topology is required [24]. Cre recombinase can carry out both DNA inversion and deletion equally well, depending on the relative orientation of its loxP target sites [25]. The synapsis of loxP sites occurs by random collision, therefore the frequency of recombination between these sites will indicate their spatial proximity [26]. This property of Cre is used here to estimate the distance between recombining loxP sites engineered within the E. coli chromosome. Other site-specific recombinases have been similarly used in the past [6,27].

The experimental strategy for assessing Cre recombination efficiency is diagrammed in Figure 1A, using as an example loxP sites flanking a Mu prophage. All assays in this study used the deletion reaction i.e. loxP pairs were configured in a direct orientation. Cre recombinase was provided from a plasmid, and reaction conditions optimized as described under Materials and Methods ( Figure S1). After recombination, the intervening DNA between the loxP sites will be excised, leaving one loxP site on the chromosome and the other site in the excised product; the latter will be lost during cell growth. As diagrammed in Figure 1A, the amounts of substrate and product, both chromosomal, were assessed by qPCR after amplification across the loxP sites with appropriate primers. Recombination efficiency (RE) was calculated as the ratio of the recombination product to the starting substrate as described in Materials and Methods.

The distance-dependence of loxP recombination on the E. coli chromosome was first assessed by varying the distance between a pair of directly oriented loxP sites from ,190 bp to 37 kbp engineered within the malF locus ( Figure S2A). This locus was chosen because the Mu prophage we wished to monitor in later experiments was located there. In the loxP-engineered strains, the log of loxP RE decreased linearly over distance, giving a first order decay function ( Figure 1B and inset; see also Figure S2D). loxP sites were next engineered on either side of a malF::Mu lysogen, ,70 bp outside each L and R end as shown in Figure 1A ( Figure 1C, this wild-type loxP-Mu-loxP construct in ZL524 is labeled MU throughout). The RE of the MU loxP sites (set at 1) was closest to the loxP-site pair placed 190 bp apart within malF in the non-Mu strain ( Figure 1C, malF). To control for recombination at distances similar to the length of the Mu prophage around this region of the chromosome, loxP sites were placed 37 kbp upstream (yjcF) or downstream (purH) of a loxP site in malF in the non-Mu strain (see Figure S2A). The RE of both these loxP pairs was similar ( Figure 1C), and reflected their linear distance as determined from the graph in Figure 1B. We conclude that reduction of the linear 37 kbp distance between the L and R Mu ends to a distance equivalent to 190 bp as measured by Cre recombination, is indicative of some form of 'synapsis' of the Mu ends.

The centrally located strong-gyrase-site (SGS) within prophage Mu is important for end-synapsis Central location of SGS ( Figure 1A) is obligatory for optimal replication of Mu after prophage induction [18]. Deletion of this site results in inefficient pairing of Mu ends in vivo, as judged by a low efficiency of transposase-mediated 39 nicking at the ends [28,29]. Pato and colleagues have proposed that the requirement for SGS in vivo but not in vitro (where the distance between the Mu ends is typically ,2 kbp on mini-Mu plasmid substrates), is an adaptation for aiding synapsis of reactive sites located at large distances [20]. We therefore tested the importance of both the presence and position of SGS on the RE of loxP sites flanking Mu. Deletion of the central SGS decreased loxP recombination 30-fold, while an asymmetric location of SGS to the left (L) or right (R) of center in separate strains showed 7-and 3-fold reduction, respectively, compared to wild-type ( Figure 2A; the RE values of 0.14260.086 and 0.30360.11 for the SGS (L) and SGS (R) prophage strains are not significantly different at p,0.01). Reduction of the Mu genome by 20 kbp upon introducing symmetrical 10 kbp deletions on either side of SGS (i.e. genome size of 17 kbp), still showed Mu end synapsis on the smaller Mu genome, which was disrupted by SGS deletion (Figure 2A). A similar reduction of effective distance was not observed when an SGS site was engineered at the center of E. coli DNA segments ranging from 5-37 kbp, each flanked by loxP ( Figure S2B, C). This shows that the SGS site alone is not sufficient to synapse distant loxP sites; Mu sequences are required in addition.

Because the effect of SGS is proposed to be mediated via gyrasepromoted supercoiling, the temperature sensitive gyrase allele gyrB402ts was introduced into a Muc + lysogen, which is not temperature inducible. Since gyrase temperature-sensitive mutants are reported to have decreased supercoiling even at the permissive temperature [30], loxP REs were measured at both permissive (26uC) and non-permissive (37uC) temperatures. loxP sites behaved as if they were unpaired only at the non-permissive temperature ( Figure 2B; the wild-type loxP-Muc + -loxP construct is labeled Mu).

We conclude that DNA supercoiling is important for reducing the distance between the Mu ends, that SGS plays a critical role in this process when located centrally either on a 37 kbp or a 17 kbp Mu genome, but that SGS does not similarly contribute when located within non-Mu E. coli DNA. These results support the Mu end-pairing function of SGS as deduced by the transposition/ replication results of Pato and colleagues. However, our data were derived in the absence of prophage induction i.e. presumably in

Author Summary

A majority of sequenced bacterial genomes harbor prophage sequences. Some prophages are viable, while others have decayed from accumulating mutations and genome rearrangements. Prophages, including defective ones, can contribute important biological properties such as antibiotic resistance, toxins, and serum resistance that increase the survival and ecological range of their hosts. We show in this study that the 37 kbp transposable prophage Mu exists in a unique configuration we call the 'Mu domain', where its two ends are paired, segregating the Mu sequences from those of the host chromosome. This is the largest stable chromosomal domain in E. coli mapped to date. The Mu domain configuration promotes low-level transcription from an early prophage promoter, which controls the expression of several genes, not all essential for phage growth. Some non-essential genes include DNA repair functions. We suggest that the Mu domain provides long-term survival benefits to both the prophage and the host: to the prophage in bestowing transposition-ready topological properties unique to the Mu reaction, and to the host in contributing extraneous DNA housekeeping functions.

the absence of the transpososome proposed to stabilize the synapsed ends.

SGS-mediated Mu DNA synapsis does not extend far outside the Mu ends

To test whether Mu ends define the base of the SGS-mediated Mu DNA synapsis loop, we moved the loxP-site pairs symmetrically from 5 kbp inside Mu (In-loxP) to 5-25 kbp outside Mu (Out-loxP) in separate strains ( Figure 3A). RE of the internal loxPs (In-5 kbp) was similar to wild-type MU, while that of external loxPs [Out-5 kbp (0.32160.12) and Out-10 kbp (0.1860.07)] decreased 3-5 fold over 5 and 10 kbp distances, which is a total linear distance of 47 and 57 kbp between the loxP pairs, respectively ( Figure 3B). Synapsis was no longer evident between the Out-25 kbp pair (87 kbp linear distance), as determined by .40 fold lower RE values (0.02360.011) compared to wild-type MU. We conclude that the SGS effect extends 5-10 kbp outside Mu ends into the flanking E. coli DNA; beyond this length, the RE is reflective of the linear rather than paired distance between the loxP sites (see standard graph in Figure 1B).

3C reveals an interaction between Mu prophage ends

In 3C analysis, protein-protein and protein-DNA crosslinking by formaldehyde is used to permanently capture interactions , and the deletion products of Cre recombination. The excised Mu prophage will be lost during cell growth. Amounts of chromosomal starting substrate and product were assessed by qPCR across the loxP sites in the substrate (primer pairs 1 or 2) and product (primer pair 3). (B) RE of loxP sites as a function of distance on the E. coli chromosome (see Figure S2A). RE (recombination efficiency) calculation is described under Materials and Methods. RE of the 190 bp loxP pair was arbitrarily set at 1 (this pair is designated malF in Figure 1C). Inset: log (base e) RE values were plotted against loxP distance, and fitted to a straight line equation y~mxzc, where y~{7:089x{4:083; the slope (m) is ,7 kbp. The plot is not quite linear because it includes two data points that fall within the plateau region of the graph (25 kbp and 37 kbp). A plot excluding the 37 kbp value is shown in Figure S2D. (C) RE of loxP sites placed ,70 bp outside each end of the wild-type malF::Mu prophage in ZL524 (MU) is set at 1, and compared to the smallest (190 bp) (malF, ZL582) and largest (37 kbp) (malF/yjcF, ZL592 and purH/malF, ZL594) loxP pairs within the malF locus of the parent non-Mu strain (see Figure S2A). Bottom panel, precise fold differences in RE. Error bars are standard deviation from the mean. doi:10.1371/journal.pgen.1003902.g001 between two genomic loci [22]. After appropriate restriction enzyme digestion and ligation at low DNA concentration, the suspected junctions can be probed by PCR using locus-specific primers. While this methodology has been widely used to generate DNA contact maps in eukaryotic cells [31], it is only beginning to be used in bacteria [23,32]. We applied this strategy to test the proximity of Mu L and R ends predicted from the Cre-loxP recombination assay.

To assay DNA interactions both inside and outside Mu, we used digestion at PstI and EcoRI sites, respectively, whose positions in malF::Mu are shown in Figure 4A. The PstI sites closest to the L and R ends inside the Mu genome are ,27 kbp apart, whereas the closest EcoRI sites outside the Mu genome are ,52 kbp apart (Lproximal site is ,13 kbp upstream and R-proximal site is ,2 kbp downstream of the prophage location). Products of the expected size were detected for both PstI and EcoRI joints ( Figure 4B,  arrowheads) in a crosslinking-ligation dependent manner. Their identities were confirmed by DNA sequencing. No PCR product corresponding to the joining of DNA cut at two EcoRI sites within Mu was detected, possibly due to cross-linked proteins interfering with ligation.

The ligation products obtained in Figure 4B were quantified further by qPCR and compared to similar reaction products from an isogenic DSGS strain ( Figure 4C). Crosslinking efficiency is defined as the ratio of qPCR signal from the products of the DSGS strain compared to those from its wild-type parent. Cross-linking efficiency of PstI ends within Mu was 4-fold (0.24760.11) higher in the presence of SGS compared to its absence ( Figure 4C, left), and that of EcoRI ends outside Mu was 2.5 fold (0.41460.106) higher under similar conditions ( Figure 4C, right). The folddifferences in crosslinking efficiencies versus recombination efficiencies, of wild-type MU and its DSGS derivative ( Figure 4C vs Figure 2A), are likely due to differences in methodology.

Background (non-specific) levels of signal generated by unligated but crosslinked samples (lanes 4 and 8 in Figure 4B) are shown in each panel (NL), along with similar controls for DNA around the malF locus in a strain where Mu had been excised ( Figure 4C, DMu).

For additional and independent quantitation of the crosslinked product, we subjected it to Cre-loxP recombination in vitro using a titrated amount of Cre. While the in vitro efficiency cannot be directly compared to the in vivo efficiency, at the lowest Cre concentration required to observe ,30% recombination in the wild-type crosslinked substrate, the RE of loxP sites in the DSGS substrate was 3-fold (0.35160.11) lower than the wild-type ( Figure 4D). Taken together, these results are an independent confirmation of the close spatial proximity of Mu ends in the Mu prophage and the important contribution of SGS to this arrangement. We shall henceforth refer to this apparent Mu-loop as a 'Mu domain'. purple and green dots, PstI and EcoRI sites, respectively; small arrows, primers used to amplify the DNA ligation product; red circle, paired L and R ends. (B) PCR products of ligation. Left: Primers were designed to produce a 155 bp fragment after PstI digestionligation (lane 2, arrowhead), and a 195 bp fragment after EcoRI digestion-ligation (lane 6, arrowhead); the fainter bands above the specific products in lanes 2 and 6 could not be re-amplified, hence are non-specific. The specific products were not observed in uncrosslinked (lanes 1, 5) and unligated (lanes 4, 8) samples. The band migrating at ,100 bp in these lanes is non-specific. Lane 3, 7 DNA size marker ladder. (C) Quantitation of the ligation products. The qPCR signal obtained from wild-type MU ligation was set at 1. Crosslinking efficiency is defined as the ratio of qPCR signal from the ligation product in the DSGS strain compared to that in its wild-type parent. NL is the signal obtained from the non-ligated, crosslinked product in the wild-type reactions shown in lanes 4 and 8, and DMu is a similar control in a strain where Mu has been excised from ZL524 via recombination of the flanking loxP sites. The same set of primer pairs were used for all strains in either the PstI or the EcoRI panels. MU (ZL524), DSGS (ZL562), DMu (ZL580). (D) In vitro Cre-loxP recombination of the cross-linked DNA from the indicated strains before digestion with restriction enzymes. doi:10.1371/journal.pgen.1003902.g004 Importance of cis-and trans-acting Mu transposition factors to the Mu domain Given that the SGS effect is Mu-specific, we wondered if the L and R ends of Mu are important for closing the Mu loop at its base, and if so, whether the Mu transposase (A protein), which binds to the Mu ends, is expressed in the prophage. We therefore individually deleted the L and R ends as well as the MuA and MuB genes from the prophage (MuB regulates MuA function allosterically; [13]). Deletion of the L end showed a 15-fold (0.06660.049) reduction of RE, but deletion of the R end had no significant effect ( Figure 5A). Deletion of the MuA gene had no effect, but deletion of MuB had an 8-fold effect (0.12560.064) ( Figure 5B). To confirm the B gene deletion result, MuB was supplied to this strain from a plasmid [pMuB (pJG8; [15])]; RE levels were restored to wild-type in this strain. The nonrequirement for the R end and for MuA, but the requirement for MuB, in Mu domain formation/maintenance, will be discussed later.

MuA and B genes are expressed from the early lytic Pe promoter, expected to be repressed in a prophage ( Figure S3A; [33]). To determine if MuB was expressed in the lysogen, we engineered into the prophage genome a functional EGFP-MuB fusion [34]. Low-level expression of MuB was detected in all the cells in this strain in the absence of Mu induction ( Figure S3B). 59-RACE-PCR experiments showed transcripts originating from Pe, as well as from a second site internal to MuA, which we have named Pe* ( Figure S3C). Pe* has significant homology to the sigma 70 sequence ( Figure S3D). Deletion of each of these promoters had a small effect on EGFP-MuB expression as measured by fluorescence, but deletion of both promoters eliminated expression ( Figure 5C). Note that either promoter deletion will render the strain uninducible for Mu lytic growth. A control EGFP fusion to a late Mu gene E expressed during lytic growth, showed no fluorescence in the prophage ( Figure 5C). Thus, the transcription results cannot be due to spontaneous Mu induction in a subpopulation of cells, because all cells expressed EGFP-MuB, none expressed EGFP-MuE, the Pe deletion interrupts the lytic growth program, and the Pe* deletion does the same by disrupting the MuA gene.

To test if the low-level expression of MuB from Pe and Pe* was related to the Mu domain, we monitored EGFP-MuB expression in a DSGS strain. Deletion of SGS diminished EGFP-MuB fluorescence compared to its wild-type parent. To test if this reduction was specific to either Pe or Pe* we measured MuB-EGFP fluorescence in the DSGS strain carrying separate Pe or Pe* deletions ( Figure 5C). Expression from Pe was the most impacted by deletion of SGS.

Single or double deletion of Pe and Pe* led to a reduction in the RE of loxP sites flanking the wild-type Mu strain, with the doublepromoter deletion showing a value similar to that of MuB gene deletion ( Figure 5D, compare with Figure 5B). As described above, the Pe-Pe* deletion disrupts the MuA gene as well as early transcription, thereby eliminating Mu lytic development. This allowed us to delete the Mu c gene encoding the lysogenic repressor Rep in the Pe-Pe* deletion strain, in order to test its contribution to the Mu domain (see Figure S3A). The data showed that absence of Rep caused an additional 3-fold decrease in RE (0.05360.028) ( Figure 5D).

We conclude that the L end is important but that the R end is dispensable to the Mu domain. MuB, but not MuA, is important for domain formation/maintenance, and Rep likely contributes as well. Low levels of MuB are expressed in the lysogen from both Pe and a newly identified promoter Pe*. The Mu domain configu-ration is important for transcription from Pe but not Pe*. Thus activity of Pe is domain-dependent, while that of Pe* is not.

Cellular NAPs are critical for maintenance of the Mu domain

E. coli NAPs such as H-NS, IHF, FIS, and HU are implicated in maintaining chromosomal supercoiled domains via their DNA bending and bridging properties. In contrast to the major E. coli NAPs, which were found largely scattered throughout the nucleoid, H-NS was reported to form two compact clusters per chromosome, sequestering and juxtaposing into these clusters numerous H-NS regulated DNA segments distributed throughout the chromosome; deleting H-NS led to substantial chromosome reorganization [23]. IHF, FIS and HU have in addition, specific binding sites on the Mu genome, from where they exert effects on Pe transcription, G-segment recombination, and transposition [13] ( Figure 6A).

To determine the importance of these NAPs to the Mu domain, we assayed for changes in RE of loxP sites flanking Mu in strains individually deleted for genes expressing these proteins ( Figure 6B). Absence of H-NS had no effect on the Mu domain. Absence of HU and of Fis had 12 to 30-fold effects, respectively (0.08460.054 and 0.03460.025). The strongest effect was observed in the absence of IHF, which essentially abrogated the Mu domain. That NAP deletions do not affect Cre recombination per se was controlled for by simultaneously monitoring the recombination of loxP site pairs placed outside the Mu domain in both wild-type and DNAP strains (see Materials and Methods).

To test if the effects of IHF, Fis and HU were exerted at the specific binding sites for these proteins on the Mu genome, we deleted these sites individually within the prophage. The dramatic reduction in RE seen in the IHF mutant was not observed with deletion of the IHF binding site ( Figure 6B). However, the 6-fold effect observed (0.1760.086), could be due to the negative effect on Pe transcription from deleting the IHF site [35][36][37][38]. Deletion of the HU-binding site at the Mu L end had a small effect, while deletion of the Fis-binding enhancer site sis (Dsis) had no effect. Recent experiments have identified a set of three Fis-binding sites within the promoter region of the mom gene near the R end [39]. A deletion spanning all three sites (DPmom) also had no effect. However, a combination of sis-attR or sis-Pmom-attR deletions reduced RE 5-7 fold (0.20960.091 and 0.14360.068). We conclude that IHF, HU and Fis affect the Mu domain configuration primarily via their global effects on chromosome structure.

The domain organization is unique to Mu, and is not observed for prophage l

The four structured macro-domain regions of the E. coli chromosome are called Ori, Right, Ter, and Left; the two less or non-structured regions (NS) are located on either side of Ori [7] ( Figure 7A). The malF::Mu prophage used in all the experiments thus far is located in the Ori macro-domain. To determine if the Mu domain organization is specific to its location in a macrodomain, we also tested for its presence in a lacZ::Muc + prophage located in the less-structured NS region between Ori and Right. The results were similar to those seen with the malF::Mu prophage, with similar negative effects of deletion of SGS or absence of IHF on the domain structure ( Figure 7B, compare to similar data in Figures 2A and 6B). Absence of H-NS had no effect at this location as well.

If the domain organization of Mu were designed to pre-engage Mu ends in a transposition-ready mode for lytic growth, might a similar arrangement be expected for other prophages that depend on pairing of their ends at the start of lytic growth? Prophage l is an example of an insertion element which must pair its attL and attR ends for excisive recombination from the E. coli chromosome [40]. The attB insertion site of l is in the Right macro-domain ( Figure 7A). To determine if l prophage ends were paired, loxP sites were engineered outside the 48.5 kbp l genome, as done for Mu. In contrast to Mu, the Cre RE of these sites reflected the linear distance between the l ends ( Figure 7B). Thus, domain organization does not occur for phage l, and may be specific for Mu.

Discussion

The 37 kbp Mu prophage domain we report in this study, inferred from Cre recombination and supported independently by a crosslinking assay, is the largest stable chromosomal domain in E. coli mapped to date. 'Stable' implies that the configuration is long-lived enough to be consistently detected by both genetic recombination and biochemical crosslinking. This sets the Mu domain apart from the dynamic configuration of the 400-500 supercoiled domains that condense the E. coli nucleoid [4]. While the nucleoid exhibits different degrees of compaction around the circular genome referred to as macro-domains [41,42], and distant chromosomal regions within these macro-domains cluster via H-NS [23], the specific supercoiled domain adopted by the Mu prophage is unique in that it represents a more or less permanent feature of the E. coli chromosome. The formation/maintenance of the Mu domain requires contributions from both prophage and its host, as discussed below.

Cre-loxP recombination as a reporter for chromosomal domains

The organization of the bacterial chromosome into supercoiled domains has been studied earlier using different strategies: trimethylpsoralen binding, electron microscopy, transcription of supercoiling sensitive genes, or site-specific recombination by Res and Int [2,3,43,44]. The Res system initially yielded an average domain size of 25 kbp for a non-essential region spanning ,100 kbp of the Salmonella typhimurium genome [43]. Recombination efficiency was found to decrease linearly with the distance between target res sites over the entire region analyzed. Since the topological constrains of the resolvase reaction requires the res sites to be housed within the same supercoiled domain [45], it was concluded that topological domains are dynamic, with stochastically distributed end points. Using an improved system in which the Res protein was designed to have a shorter half-life, the average domain size in Salmonella was re-calculated to be approximately 10 kbp [27]. A domain size of 10 kbp agrees with results obtained in E. coli, where measurements from electron microscopy and the spread of DNA relaxation from double strand breaks were used to estimate a domain size centered around 10 kbp, within a 2-66 kbp range [3]. When prophage l attL and attR sites, which can recombine within an Int synapse arranged by random collision [40,46], were placed all over the chromosome, Int recombination efficiencies suggested that accessibility of individual loci was not uniform in different regions of the E. coli and Salmonella chromosomes [6,44]. The macro-domain organization of the E. coli chromosome deduced from these and cytological studies is shown in Figure 7A [6].

The utility of the Cre-loxP system in probing chromosomal domains stems from the simple requirements of Cre recombination, the ease of integrating loxP sites at desired chromosomal locales, and the in vivo distance-dependence of the reaction revealed in this study. The ,7 kbp value of the slope derived for Cre recombination ( Figure 1B, inset), which is the distance at which there is a 50% probability that barriers to supercoil diffusion exist, defines a domain size in this region of the E. coli chromosome in reasonable agreement with the average 10 kbp domain estimate of Postow et al. [3]. Since the Mu genome is much larger than the average E. coli chromosomal domain, efficient Cre recombination at loxP sites placed at the extremities of Mu would be consistent with these sites being contained within a domain. We note that although loxP recombination by Cre is analogous to attL-attR recombination by l Int in not requiring negative supercoiling and in following the random collision mechanism, the distancedependence of recombination frequencies observed for the two systems in vivo cannot be strictly compared because of many differences in experimental conditions such as bidirectional (Cre) versus unidirectional (l Int) configurations of the recombining sites, differences in recombinase levels and reaction times, different methods for estimating REs (colony color, PCR, qPCR), and differences in growth media and growth conditions [6,44].

Sequestration of Mu into an independent supercoiled domain: Mu and host factors

The domain organization of prophage Mu, anchored by Mu L and R ends, was observed at two structurally different regions of the E. coli chromosome (Figures 1-4, and Figure 7). The pairing of the neighboring DNA arms was seen to extend 5-10 kbp outside Mu into the E. coli DNA (Figure 3). Both SGS and gyrase were critical to domain integrity (Figure 2). The requirement of SGS for the formation of the prophage domain is consistent with the role for SGS originally proposed in promoting Mu end synapsis for transposition [18,47]. In each case, the gyrase-mediated processive supercoiling initiated at the center of Mu may be cemented at the L and R ends by either the transpososome or by boundary proteins and the Mu repressor (see below) to establish two functionally distinct DNA domains.

The Mu L end and the MuB protein, but not the R end and the MuA transposase, were required for formation of the prophage domain ( Figure 5A, B). The non-requirement of the R end is puzzling. While the L and R ends have specific binding sites for MuA and for lysogenic repressor Rep [48], they also have an ATrich character [49]. We speculate that while the R end normally is required, other AT-rich elements can substitute in its absence. This possibility is suggested by the observation that an effect of Rend deletion on the Mu domain is only manifested when combined with deletion of AT-rich Fis-binding sites near this end ( Figure 6B). In the E. coli genome, AT-rich elements or A-tracts are overrepresented and distributed 'quasi-regularly' with a 10-12 bp periodicity throughout the genome, organized in ,100 bp long clusters [50]. Such elements have been proposed to constitute a 'structural code' for DNA compaction via NAP binding. Thus, in the absence of the R end, synapsis of Mu termini could be assisted by NAPs. Similarly, absence of the MuA transposase could be compensated for by the presence of the lysogenic repressor Rep, which shares sequence homology with the transposase and binds to Mu ends [48,51,52] (Figure 5D).

The requirement for MuB in domain organization/stability is consistent with the detectable but low level domain-dependent transcription from the early lytic promoter Pe and activity of a domain-independent promoter Pe* in the prophage state ( Figure  S3 and 5C). Deletion of the IHF site is expected to impact Pe transcription ( Figure 6). How might MuB assist the Mu domain? MuB is known to have a binding preference for AT-rich sequences [53,54] and can compete with NAPs for such sequences [15]. The distribution of EGFP-MuB fluorescence throughout the cell is noteworthy in this context ( Figure S3B). A potential prophagespecific NAP-like activity would be a novel role to be identified for this multifunctional protein, which affects target site selection, modulates transposase activity and promotes target site immunity during transposition [13]. We note that an earlier Res recombination study in Salmonella using a Bam Mu prophage, did not detect the Mu domain we report here [55]. It is likely that the lack of MuB impacted the results, including the possibility that the Bam Mu prophage used in that study carried in addition a large insertion of an Amp-Lac segment near the R end that might have disrupted the symmetry of the SGS site; additionally, Salmonella and E. coli chromosomes show differences in supercoiling [56].

The most abundant NAPs in E. coli -Fis, HU, H-NS and IHFengender at least partially overlapping functions, as absence of any one of these proteins results in rather subtle phenotypes [4,57]. The formation of a stable Mu domain, however, was essentially abrogated by the absence of IHF, and was strongly impeded by the absence of Fis or HU; lack of H-NS had no effect on domain establishment ( Figure 6). Although IHF regulates MuA and MuB levels from Pe [35][36][37][38], and Fis has been implicated in regulating Rep levels [58,59], while HU is required for Mu transposition [60,61], the magnitude of the effects of deleting the genes for these proteins was far greater than deleting their known binding sites. Thus, IHF, HU and Fis proteins appear to facilitate Mu domain formation independently of their site-specific interactions within the Mu genome relevant to transposition or to inversion of the Gsegment via site-specific recombination [13]. Rather, the process is likely assisted by their global role in nucleoid organization.

A model for the Mu prophage domain: Functional implications

The model for the Mu prophage domain that we propose ( Figure 8) incorporates our present findings with the earlier proposal of Pato and colleagues [20]. According to this model, the processivity of gyrase bound to the SGS site located at the center of the Mu genome, helps to align the left and right arms and promote synapsis of the L and R termini. End-binding proteins such as the transposase had been proposed earlier to seal the Mu loop and stabilize the synapse. In the prophage state, it is likely that the Mu repressor Rep rather than MuA is involved in end pairing, since more repressor molecules are expected to be present. This would explain why deletion of MuA had no effect on the Mu domain ( Figure 5B). However, Rep has a lower affinity for the ends compared to the MuA [48]. Since MuA must be expressed at a low level from the Pe transcript, a scenario where both proteins contribute to closing the Mu loop is also a plausible mechanism for ensuring domain stability. MuB likely serves as a NAP [15], assisting cellular NAPs such as FIS and IHF in stabilizing the Mu domain. Other essential NAPs such as SMC proteins could also be involved in maintaining the Mu domain [62] (Figure 8). The domain configuration promotes basal level transcription from Pe, ensuring domain maintenance via MuB.

Like Mu, prophage l depends on pairing of its attL and attR termini for excision prior to entry into lytic growth. Yet a domain organization was not detected for l (Figure 7). The contrasts between Mu and l in the architecture of their prophage genomes perhaps reflect the topological and mechanistic distinctions of the transposition and recombination reactions, respectively, which set forth each phage on the lytic path. A closed supercoiled domain may not offer a special advantage to l excision as the Int bound attL and attR sites find each other by random collision, even when present on unlinked DNA molecules [46]. In contrast, the Mu synapse is arranged by an ordered series of interactions between three sites -the L/R ends and the enhancer E -all of which must be present in cis, on the same DNA molecule [13]. The first Mu interactions between E and R, which subsequently engage L to generate an LER synapse, trap five DNA supercoils within the synapse [63]. The organization of this highly specific topological filter, that presages the chemical steps of Mu transposition, would be aided by the SGS-assisted formation of a self-contained supercoiled domain [16,64].

Prophages, despite being largely repressed in gene activity, are major contributors of genome diversity in some bacterial species [65,66]. Many of these prophages appear to be in a state of mutational decay. Both intact and defective prophages can contribute important biological properties to their bacterial hosts. The best example of 'fixation' of defective prophage genes are Shiga toxin genes in Shigella dysenteriae, sopE2 in Salmonella enterica, sspH and pertussis-like toxin genes in Salmonella typhi [66]. Some of the glycosyl transferase (gtr) genes of Salmonella may be another such example [67]. The genes for virulence factors and antibiotic resistance often carried by prophages, add to the fitness of their hosts, thereby ensuring long-term selfpropagation as well [68]. In prophage l, a small subset of its genes -rex, lom and bor -transcribed at a low level in the prophage, are involved in conferring on the host bacterium resistance to lytic phages or to serum [69,70]. In this context, the low level transcription we observe from Pe may also confer some advantage to the host ( Figures 5C and S3B, C). This transcript encompasses a number of 'semi-essential' (SE) genes whose functions are largely unknown ( Figure S3A). A subset of such genes could potentially benefit the host. The gem gene in the SE region has been reported to modulate host ligase and gyrase functions even in a lysogen [71], and gam encodes an orthologue of the eukaryotic protein Ku, which participates in double-strand break repair [72]. Pe*, whose activity is domain independent, could serve as a back-up promoter for the expression of MuB and other SE genes during brief periods when the Mu domain is disrupted, either stochastically or for functional reasons such as chromosome replication. We suggest that the positioning of SGS within the Mu genome and the SGS-induced structuring of the Mu domain are beneficial to both the phage and the host, in keeping the former transposition ready and in providing the latter the NAP-like MuB protein and proteins such as gem and gam that function in cellular physiology. The proposal that the Mu prophage confers a fitness advantage is supported by earlier competition assays under glucose-limiting conditions, that demonstrated a selective advantage for Mu and other prophage containing bacteria over their prophage-free counterparts [73].

Summary

This is the first report of a distinct organization of a prophage genome into a stable supercoiled 'loop' structure which we call the Mu domain. A closed-loop structure of the Mu genome, where synapsis of prophage ends is assisted by processive supercoiling by gyrase bound to a centrally located SGS site, had been proposed earlier to aid the transposition of Mu during its lytic phase of growth. The novel result we report is that the Mu domain exists even in the quiescent prophage state, and requires in addition to SGS, several Mu proteins and host NAPs for its formation/ maintenance. The Mu domain regulates an early promoter that controls the expression of several genes, many with unknown functions; known functions include DNA transposition, DNA repair and nucleoid structure maintenance. A domain structure likely benefits the prophage by holding its ends in a transpositionready configuration, and benefits the host by providing extra housekeeping functions. The latter proposition can be tested in long-term evolution experiments, where a Mu lysogen with an intact SGS site would be expected to outcompete one without an SGS site.

Materials and Methods

Strain construction

All strains used in this work were derivatives of E. coli K-12 and are listed in Table 1. Plasmids are listed in Table 2. Gene disruptions, substitutions, deletions and insertions on the chromosome were made using the phage l red-mediated homologous recombination methodology [74,75]. Position of mutations is listed in Table S1. Primer sequences are listed in Table S2. In strains with the temperature-inducible Mucts prophage, all incubation steps were at 30uC. All gene deletions with the kanamycin gene kan replaced the start codon to the stop codon of the gene to be deleted by amplifying the flanking regions of the gene using primers with 50 nt homology extensions using pKD4 as the kan template, and selecting for kanamycin resistance (50 mg/ml). For Mu A, B and c gene deletions, the 1.6 kbp kan cassette was retained at the site of deletion in order to maintain symmetry of the SGS site. Other deletions were created by a two-step procedure: First, the DNA to be deleted was replaced by a dual selection cassette -either cat-sacB (amplified from strain SIMD30; [76]) or kan-ccdB (amplified from pKD45, where the ccdB gene is under a rhamnose-inducible promoter; [77]) into the sequence to be replaced. Selection for the cassettes was on chloramphenicol (Cam) (100 mg/ml) or kanamycin. Next, the cassettes were replaced by homologous recombination with appropriate DNA to create the desired mutation, selecting on either LB plates supplemented with 6% sucrose or 0.5% rhamnose to eliminate cat-sacB or kan-ccdB, respectively. The gyrase mutant in the Muc + lysogen was made by moving the gyrB402ts allele [78] flanked by a cat-sacB cassette from ZL940 into ZL911 by P1 transduction, selecting for Cam R . loxP sites flanking l prophage were created by inserting them on either side of the l integration site in attB site prior to lysogenization with lcI857 (Cam) obtained from SIMD30. All constructs were confirmed by DNA sequencing. In the DHU strain, hupA is deleted [79]. In the DIHF strain, himA is deleted [80].

EGFP fusions at the N-termini of MuB and MuE were constructed by amplifying the corresponding genes from MP1999 and cloning them into BglII -SalI restriction enzyme sites on plasmid pEGFP-C1 (Clontech), which generated a 5amino acid intervening linker SGLRS. The fused genes were transferred back into the Mu prophage in MP1999 by the l Red recombination methodology. A rhamnose-inducible Cre expression vector was constructed by amplifying the gene for this recombinase from pBAD24-his-Cre plasmid and cloning into SalI-XbaI restriction enzyme sites within the rhaTRS locus on plasmid pRHA113 [81], to generate plasmid pRHA113-Cre, where Cre expression is driven from the rhaT promoter.

Cre recombination in vivo

For experiments using Cre expressed from the pBAD24-Cre (Ara) plasmid [82], M9 glucose minimal media were used in the in vivo recombination assay because in LB media, basal level leaky expression of Cre from this plasmid resulted in complete recombination by the time the cultures were grown up after plasmid transformation. For plasmid transformation, overnight (O/N) cell cultures in LB were diluted 1:100 into 20 ml of the same media, and grown at 30uC for 4-5 hr until OD 600 reached 0.6. They were washed thrice with ice cold 10% glycerol and brought to the final volume of 200 ml. 40 ml of the cells were electroporated (Biorad Gene pulser, 1.8 kV, 1 mm cuvette) with 90 ng of the plasmid. After recovery for an hour in 1 ml of minimal media at 30uC, a 1:200 dilution of the culture in the same medium with added ampicillin (100 mg/ml) was propagated at 30uC ( Figure S1A). Aliquots at different times of growth were tested for extent of recombination with and without inducer (1 mM arabinose), followed by DNA extraction from 1 ml of culture using the Wizard Genomic DNA purification kit from Promega ( Figure S1B). Recombination products were assayed by qPCR as described below. An optimal substrate recombination of ,25-30% for the wild-type malF::Mu substrate (ZL524) was observed either in early-to mid-log phase cultures (OD 600 0.5-0.6) with arabinose added for 20 min, or in late-log phase cultures (OD 600 1.2-1.3) without added inducer. However, in inducer- added mid-log phase cultures, there were large variations in the recombined fraction in different strains, whereas without added inducer, this fraction was reliably reproducible in late-log cultures of all strains. The latter conditions were therefore chosen for all the experiments reported in this study, except when testing recombination in the gyrase ts strain as described below. Later in the course of this study we acquired a rhamnoseinducible plasmid (Rha), where the basal-level leaky Cre expression was negligible. We confirmed that upon induction of Cre with rhamnose in late-log phase, recombination efficiencies (REs) were comparable to those obtained with uninduced Cre expression from the Ara plasmid ( Figure S1C); REs obtained with the rhamnose-induced Cre were comparable between mid-and late-log cultures ( Figure S1D).

Cre assays in the gyrase ts mutant were carried out with the Rha plasmid. Wild-type Muc+ (ZL911) and its isogenic gyrase ts mutant (ZL941) strains transformed with the pRHA113-Cre plasmid were grown in LB media (supplemented with 0.2% glucose) at 26uC until OD 600 reached 0.5, then either maintained 26uC or shifted to the non-permissive temperature 37uC for 30 min [78], followed by addition of 1 mM rhamnose for 20 min before DNA extraction.

To control for the effect of NAPs on Cre recombination per se, RE of loxP site pairs placed outside the Mu domain was monitored simultaneously in both wild-type and DNAP strains. The REs were not affected by deletion of any of the NAPs tested (Zheng Lou, Ph.D. dissertation).

Real-time qPCR

Aliquots with 50 ng of DNA, 10 ml SYBR master mix (Applied Biosystems Inc; includes dNTPs, enzyme and buffer), 0.4 ml of each primer (10 mM) and 8.2 ml of double distilled H 2 O were held for 10 min at 95uC, followed by 40 cycles of 15 sec at 95uC and 1 min at 60uC (7900HT; Applied Biosystems). Three independent biological replicates were tested, and for each biological replicate three independent technical replicates were performed. Product integrity was checked using the dissociation curve. Cycle Threshold (C t ) was read out, and the starting template amount was quantified based on the value of C t assuming exponential growth at early stages of amplification.

Recombination efficiencies were calculated based on the threshold cycle (Ct). The relative threshold cycles of each sample were calculated as, DCt~Ct  Table S2).

Primer efficiency of a primer pair (say A & B) was determined as follows: Primer A was linked to pUC19t (forward) and primer B was linked to pUC19b (reverse); the pUC primers anneal to pUC19 plasmid and amplify a common 180 bp fragment. The PCR products were purified by Qiaquick PCR purification KitH (Qiagen) and used as templates for qPCR in the following reaction: 12.5 ml SYBR mix (Qiagen), 0.75 ml each of primer A (10 mM) and primer B (10 mM), 1 ml of template (10 ng/ml) (the 180 bp fragment, as described above) and 10 ml of double distilled H 2 O. PCR cycles were as described above. Another qPCR reaction was performed using the internal primer pair pUC19t and pUC19b. The primer efficiency of primer pairs A/B was calculated as the ratio of Ct values of the PCR product obtained using primer A-pUC19t+primer B-pUC19b to that from primers pUC19t+-pUC19b.

Chromosome Conformation Capture (3C) assay

The methodology was modified from published protocols [83,84]. An O/N cell culture in Luria broth (LB) was diluted 1:1000 into 50 ml of fresh medium and grown with shaking at 30uC until OD 600 reached 0.5-0.6.

DNA crosslinking and cell lysis. 1.35 ml of formaldehyde (Fisher-Scientific, 37%) was added to the 50 ml cell culture and incubated for 20 min at room temperature with slow shaking. Crosslinking was stopped with addition of 12 ml of 2.5 M glycine for 5 min at room temperature (r.t.). Cells were centrifuged at 50006 g for 10 min at 4uC, washed twice in 10 ml ice-cold PBS pH 7.5, resuspended in 1 ml PBS and divided into 100 ml aliquots in 10 eppendorf tubes. The aliquots were centrifuged at 140006 g 4uC for 5 min and the pellets resuspended in 200 ml of 16 restriction enzyme buffer. The centrifugation step was repeated, and the pellet resuspended in 192 ml digestion buffer (16 restriction enzyme buffer, 16 complete protease inhibitor from Roche). Next, 2 ml of 35 KU/ml Ready-Lyse lysozyme (Epicenter Biotechnologies) was added to the suspension and incubated at r.t. for 20 min. Finally, 6 ml of 10% SDS solution was added, followed by O/N incubation at 37uC.

Digestion and ligation. 100 ml of the cell lysate prepared above was added to 200 ml digestion buffer (see above) and 30 ml of 20% Triton X-100. After incubation at 37uC for 1.5 hr, restriction digestion was performed thrice at 37uC as follows: 3 ml of 100 U/ml enzyme EcoRI or PstI (New England BioLabs) for 3 hr, another 3 ml of enzyme for 3 hr, and finally another 2 ml of enzyme for O/N incubation.

On the following day, digestion was stopped by addition of 50 ml of 10% SDS, incubated at 65uC for 20 min. Samples were centrifuged at 16,000 g for 5 min at r.t. and their supernatants (,390 ml) transferred to 15 ml tubes, which contained 4095 ml of pre-ligation buffer (25 ml of 20% Triton X-100, 100 ml of 256 complete protease inhibitor, 50 ml of 1M Tris-HCl pH 7.8, 50 ml of 10 mg/ml BSA, and 3645 ml of H 2 O). After incubation for 1 hr at 37uC with shaking, 500 ml of 106 T4 DNA ligase buffer and 15 ml of 2000 U/ml T4 DNA ligase (New England BioLabs) were added to a final DNA concentration of 0.8 ng/ml. The ligation reaction was incubated at 4uC for 3 days, replenishing ATP each day by addition of 50 ml of 100 mM ATP. Finally, the ligation reaction was incubated for 1 hr at r.t., followed by addition of Proteinase K mixture (210 ml of 5M NaCl, 10 ml of 0.5 M EDTA, 105 ml of 10 mg/ml Proteinase K) to stop the reaction. Crosslinks were reversed by O/N incubation at 65uC.

DNA purification. To remove RNA from the samples, 30 ml of 10 mg/ml RNase A (Promega) was added for 45 min at 37uC. DNA was extracted twice by an equal volume of phenolchloroform-isoamyl alcohol pH 6.7 (Fisher-Scientific), and with pure chloroform once. DNA was precipitated by addition of glycogen (final concentration 50 mg/ml, Affymetrix), 500 ml 3M sodium acetate, and 10 ml isopropyl alcohol, per 5 ml DNA solution. The air-dried DNA pellet was resuspended in 100 ml of 10 mM Tris-HCl pH 7.5 for qPCR quantification as described above. All qPCR results were validated by regular PCR, DNA electrophoresis, and DNA sequencing.

Cre recombination in vitro

Cre recombination was performed on genomic DNA that was crosslinked with formaldehyde and treated as described above, just prior to addition of restriction enzymes. 10 ml of Cre (1 mg/ml) was added for 3 hr at 37uC (,30% recombination in wild-type prophage DNA substrate). Cre protein was a gift from Dr. Makkuni Jayaram [85]. DNA was treated with RNAse and precipitated as described above under 'DNA purification'.

RNA Isolation and RACE

For RNA isolation, cells were grown with shaking at 30uC in 10 ml of LB until OD 600 reached 0.6. Two ml of culture (,1610 8 cells) were harvested for RNA isolation using ToTALLY RNA Kit from Ambion according to their specification. MICROBExpress Kit from Ambion was used to enrich for mRNA from 10 mg of purified total RNA by removing the 16S and 23S ribosomal RNAs (rRNAs). The final yield of enriched mRNA was ,1 mg. The quality of total RNA was checked by agarose gel electrophoresis and the RNA concentration was determined by measuring OD at 260 and 280 nm. RNA samples were stored at 280uC until use.

59-RACE (59-Rapid Amplification of cDNA Ends) was used to determine the transcriptional start sites, using the SMARTer RACE cDNA Amplification Kit (Clontech, Mountain View, CA). Gene-specific primers (GSP1-3 for Pe* and GSP4-5 for Pe promoter) were used to amplify the 59-end of isolated mRNA. The RACE PCR amplified bands were gel-purified and sequenced directly.

Fluorescence measurement of EGFP-MuB/MuE strains

The fluorescence intensity of EGFP strains was recorded using a PTI Quanta Master Model C scanning spectrofluorometer. Strains were sub-cultured by 1:50 dilution from an O/N culture into 5 ml of LB, and grown at 30uC until OD 600 reached 0.6. Three ml of the culture cell were placed in a Bio-Rad VersaFluor cuvette with a path length of 10 mm. Excitation and emission wavelengths were set at 488 nm and 507 nm, respectively. Fluoresence measurements were obtained from three independent cultures propagated on different days, each measured in triplicate. AU is arbitrary units. Each fluorescent value was derived by calculating fluorescent data of EGFP-MuB/MuE strain minus fluorescent data of WT strain without EGFP fusion.

Fluorescence microscopy

2 ml of culture prepared as described above was placed on a glass slide, and examined with Olympus BX53 microscope equipped with a GFP filter. Photographs were taken with Olympus XM10 camera and processed with Photoshop (Adobe Systems, Palo Alto, CA). Figure S1 (A) Growth curves of the wild-type loxP-malF::Mu-loxP strain ZL524 (MU) at 30uC in M9 glucose minimal media, containing one of two Cre-expressing plasmids. ZL524 contains either the arabinose-inducible Cre plasmid pBAD24-His-Cre (Ara), or the rhamnose-inducible Cre plasmid pRHA113-Cre (Rha). Growth curves are monitored without added inducer. (B) Cre recombination in ZL524 (MU) carrying the Ara plasmid was estimated at different times during the growth curve shown in (A). Recombination was assessed either in the absence (2) or presence (+) of 1 mM arabinose inducer, added for 20 min. Percentage of substrate recombined was estimated by measuring the initial (without Cre plasmid) and final substrate concentration (with Cre plasmid) using qPCR. The data are derived from three technical repeats of three biologically independent samples. (C) Strains containing either the Ara or the Rha plasmid were propagated in minimal media for ,40 hr. No inducer was added in strains with the Ara plasmid; 1 mM rhamnose was added for 20 min to strains with the Rha plasmid. Recombination efficiency (RE) was calculated as described in Materials and Methods. RE of loxP sites flanking the wild-type malF::Mu prophage is set to 1 (MU, ZL524), and compared to a pair of loxP sites separated by an equivalent 37 kbp on chromosomal DNA in malF region (malF/ yjcF, ZL592). (D) As in (C), except at ,30 hr of growth. (TIF) Figure S2 (A) Position of pairs of loxP sites at different distances around the malF locus on the E. coli chromosome. These loxP pairs were engineered after excision of Mu from malF in ZL524. Strain numbers for the indicated distances (shown in parentheses) between loxP pairs are: malF, ZL582 (190 bp); malF-lamB, ZL704 (5 kbp); malF-ubiC, ZL708 (9 kbp); malF-yjbJ, ZL810 (15 kbp); malF-yjbY, RS059 (17 kbp); malF-yjbS, ZL712 (25 kbp); malF-yjcF, ZL592 (37 kbp); purH-malF, ZL594 (37 kbp). Strains are listed in Table 1, and the exact position of loxP sites is found in Table S1. Primers used are listed in Table S2 , and Pe controls a long early transcript from ner to C, which encodes not only the transposition functions A and B, but also largely uncharacterized functions in the semi-essential (SE) region. Pm and Plys are late promoters active during lytic growth. Pe* is a new promoter identified in this study [49,86]. (B) EGFP-MuB fluorescence in strains containing Mu prophages without (WT) or with (EGFP-MuB) EGFP fused to the B. Both strains were grown at 30uC, where the prophage does not enter lytic growth. WT (MP1999), EGFP-MuB (RS033). (C) Characterization of the Mu early transcripts in a Mu lysogen. Total RNA was isolated from the uninduced strain MP1999. 59-RACE-PCR was performed on the first-strand cDNA synthesized from the total RNA using primers within MuA and MuB genes, and the products were directly sequenced to identify the 59 ends as described in Materials and Methods. Two products were initially obtained -Pe and Pe*. These were characterized separately using gene-specific primers (GSPs) placed are varying distances to confirm that the size of the product varied as predicted from the identity of the 59 terminus. GSP positions on the Mu genome are shown in the schematic below. Lanes 3 and 7 contain DNA size markers. (D) Position of Pe and Pe* with respect to ner and A gene ORFs and their homology to the E. coli sigma 70 promoter consensus sequence. Transcription start sites as determined by 59RACE are indicated by magenta coloring of the A nucleotide starts determined for both transcripts. Start of the Pe transcript matches that reported previously by S1 mapping [87]. Conserved nucleotides in both promoters are underlined. Compared to the sigma 70 consensus promoter, Pe* has the same number of conserved nucleotides as found in the Pe i.e. 5/6 at 235 and 4/6 at 210. (TIF)

Supporting Information

Figure 1 .
1Cre recombination assay and loxP recombination as a function of distance. (A) Schematic shows directly oriented loxP sites flanking prophage Mu attL and attR (L/R in text)

Figure 2 .
2Importance of SGS and DNA supercoiling to Mu end synapsis. (A) RE of loxP sites flanking Mu in SGS-deleted, SGS-displaced to left (L) or right (R) of center, and reduced-Mu genome strains. MU (ZL524), DSGS (ZL562), SGS(L) (ZL573), SGS(R) (ZL578), 17 kb Mu (RS020), 17 kb DSGS (RS025). (B) RE of loxP sites flanking a Muc+ lysogen (ZL911) and its isogenic gyrBts strain (ZL941) measured at 26uC and 37uC, using the rhamnoseinducible Cre plasmid as described under Materials and Methods. Other descriptions as in Figure 1. doi:10.1371/journal.pgen.1003902.g002

Figure 3 .
3Propagation of the Mu domain outside Mu ends. (A) Position of loxP-site pairs inside (In) or outside (Out) Mu. The exact location of the sites is given in Table S1. (B) RE of loxP-site pairs diagrammed in A. In-5 kb, Out-5 kb etc. refers to RE of pairs of symmetrically placed sites within and outside the L and R ends of Mu in different strains at the indicated distances. MU (ZL524), In-5 kb (ZL732), Out-5 kb (ZL720), Out-10 kb (ZL724), Out-25 kb (ZL728). doi:10.1371/journal.pgen.1003902.g003

Figure 4 .
4Interaction of prophage Mu ends probed by 3C methodology. (A) Experimental design (see text and Materials and Methods). Blue line, Mu DNA; black line, E. coli DNA;

Figure 5 .
5Cisand trans-acting Mu transposition factors required for Mu domain formation, and discovery of early rightward transcripts in the prophage. (A) RE of loxPs in DattL (ZL552) and DattR (ZL556) strains compared to their wild-type parent MU. (B) RE of loxPs in strains deleted for MuA (ZL536) and MuB (ZL530) genes, and MuDB strain complemented with MuB from plasmid pJG8 (pMuB). (C) Expression of EGFP-MuB in a wild-type (WT) Mu lysogen and isogenic strains carrying SGS and early promoter deletions. Excitation and emission wavelengths were 488 nm and 507 nm, respectively. Fluorescence values were subtracted from the background fluorescence of the parental Mu strain without the EGFP-MuB (MP1999), and expressed as arbitrary units (AU). Error bars are standard deviation from the mean. Strains: Wild-type Mucts prophage expressing EGFP-MuB (RS033) and its deletion derivatives DPe (RS048), DPe* (RS102), DPe DPe* (RS103), DSGS (RS088), DSGS DPe (RS106), DSGS DPe* (RS107) and EGFP-MuE (RS101). (D) RE of loxPs in Mu prophage strains carrying early promoter deletions. DPe (RS053), DPe* (RS092), DPe DPe* (RS093), DRep DPe DPe* (ZL951). In panels A, B and D, MU is ZL524. doi:10.1371/journal.pgen.1003902.g005

Figure 6 .
6Role of E. coli NAPs and their binding sites on the Mu genome, in Mu domain configuration. (A) Schematic showing binding sites for IHF, Fis and HU on the Mu genome. O/E, operator/enhancer; G, invertible G segment; sis, Fis-binding enhancer site for G inversion; Pmom, promoter for the mom gene. (B) RE of loxP sites in strains either deleted for the indicated NAPs or for their binding sites on Mu. Strains: MU (ZL524), DH-NS (ZL624), DHU (ZL634), DFis (ZL614), DIHF (ZL604), DIHF site (ZL656), DHU site (ZL652), Dsis (ZL660), DPmom (ZL670), Dsis DattR (ZL662), Dsis DPmom DattR (ZL672). doi:10.1371/journal.pgen.1003902.g006

Figure 7 .
7Mu domain at an NS chromosomal location, and probing for a domain configuration for prophage l. (A) E. coli chromosome macro-domains and position of prophages at the loci examined in this study. (B) RE of loxP sites flanking Muc+ prophage (Mu) located in lacZ, and its isogenic mutants. Mu (ZL911), DSGS (ZL921), DH-NS (ZL931), DIHF (ZL936). (C) RE of loxP sites flanking prophage l (ZL808) compared to Mu (ZL911). doi:10.1371/journal.pgen.1003902.g007

Figure 8 .
8Model of the Mu prophage domain. See text for description. The single supercoiled Mu loop shown is not intended to imply absence of branching. NAPs not tested in this study, such as SMC-like proteins, may also be involved in domain maintenance. A variation of this model was proposed earlier for replicating Mu, to account for strong MuB binding only within Mu[47]. doi:10.1371/journal.pgen.1003902.g008

(P){Ct(S) where (P) represents recombination product and (S) the substrate before recombination. The recombination efficiencies (RE) of different samples are normalized to set the recombination efficiency of the control loxP sites as 1. The relative recombination efficiency of each sample is calculated as RE~2 DCt(WT){DCt(Mutant) . The primer pair for amplifying the starting substrates were (1) 'RT attL loxP t and RT attL loxP b' for malF::Mu in MP1999, (2) 'RT LR loxP t and RT LR loxP b' for malF without Mu, (3) 'RT lacZ L t and RT lacZ L b' for lacZ::Mu, and (4) 'RT l L t and RT l L b' for l (other primers are listed in

. (B) The SGS site was engineered at the center of the 37 kbp yjcF-malF E. coli DNA segment (see A) and RE of flanking loxP sites measured. MU (ZL524), yjcF/malF (ZL592), yjcF/ malF+SGS (ZL598). (C) Effect of SGS on the RE of loxP site pairs at varying distances in E. coli DNA. The SGS site was introduced at the center of DNA flanked by loxP pairs separated by 5-37 kbp shown in A. The RE of these sites is compared in strains with (red) and without (blue) SGS. The strains without SGS are listed in A. Those with added SGS are: ZL706 (5 kbp), ZL710 (9 kbp), ZL714 (25 kbp), ZL598 (37 kbp). (D) Double log plot of RE vs distance as described for the data in Figure 1B, except that the RE value at 37 kbp is omitted. Here, y~{5:41x{1:86. (TIF) Figure S3 Genetic map of the left end of Mu and identification of a new promoter Pe*. (A) PcM and Pe are divergent promoters that control the lysogeny-lysis decision; PcM drives transcription of the lysogenic repressor gene c (Rep protein)

Table 1 .
1Strains.Strains
Genotype*

Table 1 .
1Cont.Strains
Genotype*

Table S1
S1Exact position of insertions and deletions in Mu and in E. coli. * For Mu, the numbers indicate nucleotide positions starting at 1 at the L end of Mu; for E. coli, they indicate nucleotide positions starting at 1 on the E. coli genome. (DOCX)

Table S2
S2Primers used. (DOCX)
PLOS Genetics | www.plosgenetics.org
November 2013 | Volume 9 | Issue 11 | e1003902
AcknowledgmentsWe thank M. Jayaram for helpful comments.Author ContributionsConceived and designed the experiments: RPS ZL LM RMH. Performed the experiments: RPS ZL LM. Analyzed the data: RPS ZL LM RMH. Contributed reagents/materials/analysis tools: RPS ZL LM RMH. Wrote the paper: RMH.
Bacterial chromatin. A Travers, G Muskhelishvili, Curr Opin Genet Dev. 15Travers A, Muskhelishvili G (2005) Bacterial chromatin. Curr Opin Genet Dev 15: 507-514.

Chromosomes in living Escherichia coli cells are segregated into domains of supercoiling. R R Sinden, D E Pettijohn, Proc Natl Acad Sci U S A. 78Sinden RR, Pettijohn DE (1981) Chromosomes in living Escherichia coli cells are segregated into domains of supercoiling. Proc Natl Acad Sci U S A 78: 224-228.

Topological domain structure of the Escherichia coli chromosome. L Postow, C D Hardy, J Arsuaga, N R Cozzarelli, Genes Dev. 18Postow L, Hardy CD, Arsuaga J, Cozzarelli NR (2004) Topological domain structure of the Escherichia coli chromosome. Genes Dev 18: 1766-1779.

Bacterial nucleoid-associated proteins, nucleoid structure and gene expression. S C Dillon, C J Dorman, Nat Rev Microbiol. 8Dillon SC, Dorman CJ (2010) Bacterial nucleoid-associated proteins, nucleoid structure and gene expression. Nat Rev Microbiol 8: 185-195.

Architectural organization in E. coli nucleoid. M Macvanin, S Adhya, Biochim Biophys Acta. 1819Macvanin M, Adhya S (2012) Architectural organization in E. coli nucleoid. Biochim Biophys Acta 1819: 830-835.

Macrodomain organization of the Escherichia coli chromosome. M Valens, S Penaud, M Rossignol, F Cornet, F Boccard, EMBO J. 23Valens M, Penaud S, Rossignol M, Cornet F, Boccard F (2004) Macrodomain organization of the Escherichia coli chromosome. EMBO J 23: 4330-4341.

Spatial arrangement and macrodomain organization of bacterial chromosomes. F Boccard, E Esnault, M Valens, Mol Microbiol. 57Boccard F, Esnault E, Valens M (2005) Spatial arrangement and macrodomain organization of bacterial chromosomes. Mol Microbiol 57: 9-16.

On the structure of the folded chromosome of Escherichia coli. A Worcel, E Burgi, J Mol Biol. 71Worcel A, Burgi E (1972) On the structure of the folded chromosome of Escherichia coli. J Mol Biol 71: 127-147.

A genetic selection for supercoiling mutants of Escherichia coli reveals proteins implicated in chromosome structure. C D Hardy, N R Cozzarelli, Mol Microbiol. 57Hardy CD, Cozzarelli NR (2005) A genetic selection for supercoiling mutants of Escherichia coli reveals proteins implicated in chromosome structure. Mol Microbiol 57: 1636-1652.

The architectural role of nucleoid-associated proteins in the organization of bacterial chromatin: a molecular perspective. M S Luijsterburg, M C Noom, G J Wuite, R T Dame, J Struct Biol. 156Luijsterburg MS, Noom MC, Wuite GJ, Dame RT (2006) The architectural role of nucleoid-associated proteins in the organization of bacterial chromatin: a molecular perspective. J Struct Biol 156: 262-272.

DNA supercoiling and its role in DNA decatenation and unknotting. G Witz, A Stasiak, Nucleic Acids Res. 38Witz G, Stasiak A (2010) DNA supercoiling and its role in DNA decatenation and unknotting. Nucleic Acids Res 38: 2119-2133.

Parallel genetic and phenotypic evolution of DNA superhelicity in experimental populations of Escherichia coli. E Crozat, C Winkworth, J Gaffe, P F Hallin, M A Riley, Mol Biol Evol. 27Crozat E, Winkworth C, Gaffe J, Hallin PF, Riley MA, et al. (2010) Parallel genetic and phenotypic evolution of DNA superhelicity in experimental populations of Escherichia coli. Mol Biol Evol 27: 2113-2128.

Transposition of phage Mu DNA. G Chaconas, R M Harshey, Chaconas G, Harshey RM (2002) Transposition of phage Mu DNA;

. N L Craig, R Craigie, M Gellert, A M Lambowitz, ASM PresspWashington DCCraig NL, Craigie R, Gellert M, Lambowitz AM, editors. Washington DC: ASM Press. 384-402 p.

Microarray analysis of transposition targets in Escherichia coli: the impact of transcription. D Manna, A M Breier, N P Higgins, Proc Natl Acad Sci U S A. 101Manna D, Breier AM, Higgins NP (2004) Microarray analysis of transposition targets in Escherichia coli: the impact of transcription. Proc Natl Acad Sci U S A 101: 9780-9785.

Analysis of phage Mu DNA transposition by whole-genome Escherichia coli tiling arrays reveals a complex relationship to distribution of target selection protein B, transcription and chromosome architectural elements. J Ge, Z Lou, H Cui, L Shang, R M Harshey, J Biosci. 36Ge J, Lou Z, Cui H, Shang L, Harshey RM (2011) Analysis of phage Mu DNA transposition by whole-genome Escherichia coli tiling arrays reveals a complex relationship to distribution of target selection protein B, transcription and chromosome architectural elements. J Biosci 36: 587-601.

The Mu transpososome through a topological lens. R M Harshey, M Jayaram, Crit Rev Biochem Mol Biol. 41Harshey RM, Jayaram M (2006) The Mu transpososome through a topological lens. Crit Rev Biochem Mol Biol 41: 387-405.

DNA gyrase requirements distinguish the alternate pathways of Mu transposition. T D Sokolsky, T A Baker, Mol Microbiol. 47Sokolsky TD, Baker TA (2003) DNA gyrase requirements distinguish the alternate pathways of Mu transposition. Mol Microbiol 47: 397-409.

Central location of the Mu strong gyrase binding site is obligatory for optimal rates of replicative transposition. M L Pato, Proc Natl Acad Sci U S A. 91Pato ML (1994) Central location of the Mu strong gyrase binding site is obligatory for optimal rates of replicative transposition. Proc Natl Acad Sci U S A 91: 7056-7060.

Mu-like prophage strong gyrase site sequences: analysis of properties required for promoting efficient Mu DNA replication. M Oram, M L Pato, J Bacteriol. 186Oram M, Pato ML (2004) Mu-like prophage strong gyrase site sequences: analysis of properties required for promoting efficient Mu DNA replication. J Bacteriol 186: 4575-4584.

The Mu strong gyrase-binding site promotes efficient synapsis of the prophage termini. M L Pato, M Banerjee, Mol Microbiol. 22Pato ML, Banerjee M (1996) The Mu strong gyrase-binding site promotes efficient synapsis of the prophage termini. Mol Microbiol 22: 283-292.

Path of DNA within the Mu transpososome. Transposase interactions bridging two Mu ends and the enhancer trap five DNA supercoils. S Pathania, M Jayaram, R M Harshey, Cell. 109Pathania S, Jayaram M, Harshey RM (2002) Path of DNA within the Mu transpososome. Transposase interactions bridging two Mu ends and the enhancer trap five DNA supercoils. Cell 109: 425-436.

Capturing chromosome conformation. J Dekker, K Rippe, M Dekker, N Kleckner, Science. 295Dekker J, Rippe K, Dekker M, Kleckner N (2002) Capturing chromosome conformation. Science 295: 1306-1311.

Chromosome organization by a nucleoid-associated protein in live bacteria. W Wang, G W Li, C Chen, X S Xie, X Zhuang, Science. 333Wang W, Li GW, Chen C, Xie XS, Zhuang X (2011) Chromosome organization by a nucleoid-associated protein in live bacteria. Science 333: 1445-1449.

The nature of the interaction of the P1 recombinase Cre with the recombining site loxP. R Hoess, K Abremski, N Sternberg, Cold Spring Harb Symp Quant Biol. 49Hoess R, Abremski K, Sternberg N (1984) The nature of the interaction of the P1 recombinase Cre with the recombining site loxP. Cold Spring Harb Symp Quant Biol 49: 761-768.

A structural view of tyrosine recombinase site-specific recombination. Gd ; Van Duyne, N L Craig, R Craigie, M Gellert, A M Lambowitz, ASM PresspWashington DCVan Duyne GD (2002) A structural view of tyrosine recombinase site-specific recombination; Craig NL, Craigie R, Gellert M, Lambowitz AM, editors. Washington DC: ASM Press. 93-117 p.

Studies on the properties of P1 sitespecific recombination: evidence for topologically unlinked products following recombination. K Abremski, R Hoess, N Sternberg, Cell. 32Abremski K, Hoess R, Sternberg N (1983) Studies on the properties of P1 site- specific recombination: evidence for topologically unlinked products following recombination. Cell 32: 1301-1311.

Measuring chromosome dynamics on different time scales using resolvases with varying half-lives. R A Stein, S Deng, N P Higgins, Mol Microbiol. 56Stein RA, Deng S, Higgins NP (2005) Measuring chromosome dynamics on different time scales using resolvases with varying half-lives. Mol Microbiol 56: 1049-1061.

Characterization of Mu prophage lacking the central strong gyrase binding site: localization of the block in replication. M L Pato, M Karlok, C Wall, N P Higgins, J Bacteriol. 177Pato ML, Karlok M, Wall C, Higgins NP (1995) Characterization of Mu prophage lacking the central strong gyrase binding site: localization of the block in replication. J Bacteriol 177: 5937-5942.

Replication of Mu prophages lacking the central strong gyrase site. M L Pato, Res Microbiol. 155Pato ML (2004) Replication of Mu prophages lacking the central strong gyrase site. Res Microbiol 155: 553-558.

Rates of gyrase supercoiling and transcription elongation control supercoil density in a bacterial chromosome. N Rovinskiy, A A Agbleke, O Chesnokova, Z Pang, N P Higgins, PLoS Genet. 81002845Rovinskiy N, Agbleke AA, Chesnokova O, Pang Z, Higgins NP (2012) Rates of gyrase supercoiling and transcription elongation control supercoil density in a bacterial chromosome. PLoS Genet 8: e1002845.

A decade of 3C technologies: insights into nuclear organization. E De Wit, W De Laat, Genes Dev. 26de Wit E, de Laat W (2012) A decade of 3C technologies: insights into nuclear organization. Genes Dev 26: 11-24.

The three-dimensional architecture of a bacterial genome and its alteration by genetic perturbation. M A Umbarger, E Toro, M A Wright, G J Porreca, D Bau, Mol Cell. 44Umbarger MA, Toro E, Wright MA, Porreca GJ, Bau D, et al. (2011) The three-dimensional architecture of a bacterial genome and its alteration by genetic perturbation. Mol Cell 44: 252-264.

Regulation of transcription. N Goosen, P Van De Putte, Symonds N, Toussaint A, Van de Putte P, Howe MMPhage Mu. Cold Spring HarborNY: Cold Spring Harbor LaboratoryGoosen N, van de Putte P (1987) Regulation of transcription. In: Symonds N, Toussaint A, Van de Putte P, Howe MM, editors. Phage Mu. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory. pp. 41-52.

Dynamics of a protein polymer: the assembly and disassembly pathways of the MuB transposition target complex. E C Greene, K Mizuuchi, EMBO J. 21Greene EC, Mizuuchi K (2002) Dynamics of a protein polymer: the assembly and disassembly pathways of the MuB transposition target complex. EMBO J 21: 1477-1486.

Regulation of Mu transposition. II. The Escherichia coli HimD protein positively controls two repressor promoters and the early promoter of bacteriophage Mu. N Goosen, M Van Heuvel, G F Moolenaar, P Van De Putte, Gene. 32Goosen N, van Heuvel M, Moolenaar GF, van de Putte P (1984) Regulation of Mu transposition. II. The Escherichia coli HimD protein positively controls two repressor promoters and the early promoter of bacteriophage Mu. Gene 32: 419-426.

Positive and negative regulation of the Mu operator by Mu repressor and Escherichia coli integration host factor. H M Krause, N P Higgins, J Biol Chem. 261Krause HM, Higgins NP (1986) Positive and negative regulation of the Mu operator by Mu repressor and Escherichia coli integration host factor. J Biol Chem 261: 3744-3752.

Integration host factor of Escherichia coli regulates early-and repressor transcription of bacteriophage Mu by two different mechanisms. P A Van Rijn, N Goosen, P Van De Putte, Nucleic Acids Res. 16van Rijn PA, Goosen N, van de Putte P (1988) Integration host factor of Escherichia coli regulates early-and repressor transcription of bacteriophage Mu by two different mechanisms. Nucleic Acids Res 16: 4595-4605.

Supercoiling and integration host factor change the DNA conformation and alter the flow of convergent transcription in phage Mu. N P Higgins, D A Collier, M W Kilpatrick, H M Krause, J Biol Chem. 264Higgins NP, Collier DA, Kilpatrick MW, Krause HM (1989) Supercoiling and integration host factor change the DNA conformation and alter the flow of convergent transcription in phage Mu. J Biol Chem 264: 3035-3042.

Silencing of toxic gene expression by Fis. S Karambelkar, G Swapna, V Nagaraja, Nucleic Acids Res. 40Karambelkar S, Swapna G, Nagaraja V (2012) Silencing of toxic gene expression by Fis. Nucleic Acids Res 40: 4358-4367.

. M A Azaro, A Landy, American Society for MicrobiologypWashington, D.C.Azaro MA, Landy A, editors (2002) l integrase and the l Int family. Washington, D.C.: American Society for Microbiology. 118-148 p.

Organization of the Escherichia coli chromosome into macrodomains and its possible functional implications. O Espeli, F Boccard, J Struct Biol. 156Espeli O, Boccard F (2006) Organization of the Escherichia coli chromosome into macrodomains and its possible functional implications. J Struct Biol 156: 304- 310.

Escherichia coli and its chromosome. R Reyes-Lamothe, X Wang, D Sherratt, Trends Microbiol. 16Reyes-Lamothe R, Wang X, Sherratt D (2008) Escherichia coli and its chromosome. Trends Microbiol 16: 238-245.

Surveying a supercoil domain by using the gamma delta resolution system in Salmonella typhimurium. N P Higgins, X Yang, Q Fu, J R Roth, J Bacteriol. 178Higgins NP, Yang X, Fu Q, Roth JR (1996) Surveying a supercoil domain by using the gamma delta resolution system in Salmonella typhimurium. J Bacteriol 178: 2825-2835.

Unequal access of chromosomal regions to each other in Salmonella: probing chromosome structure with phage lambda integrase-mediated longrange rearrangements. N Garcia-Russell, T G Harmon, T Q Le, N H Amaladas, R D Mathewson, Mol Microbiol. 52Garcia-Russell N, Harmon TG, Le TQ, Amaladas NH, Mathewson RD, et al. (2004) Unequal access of chromosomal regions to each other in Salmonella: probing chromosome structure with phage lambda integrase-mediated long- range rearrangements. Mol Microbiol 52: 329-344.

Rate and selectively of synapsis of res recombination sites by Tn3 resolvase. M A Watson, M R Boocock, W M Stark, J Mol Biol. 257Watson MA, Boocock MR, Stark WM (1996) Rate and selectively of synapsis of res recombination sites by Tn3 resolvase. J Mol Biol 257: 317-329.

Site-specific recombination of bacteriophage lambda. The change in topological linking number associated with exchange of DNA strands. H A Nash, T J Pollock, J Mol Biol. 170Nash HA, Pollock TJ (1983) Site-specific recombination of bacteriophage lambda. The change in topological linking number associated with exchange of DNA strands. J Mol Biol 170: 19-38.

Immunity of replicating Mu to self-integration: a novel mechanism employing MuB protein. J Ge, Z Lou, R M Harshey, Mobile DNA. 18Ge J, Lou Z, Harshey RM (2010) Immunity of replicating Mu to self-integration: a novel mechanism employing MuB protein. Mobile DNA 1: 8.

Site-specific recognition of the bacteriophage Mu ends by the Mu A protein. R Craigie, M Mizuuchi, K Mizuuchi, Cell. 39Craigie R, Mizuuchi M, Mizuuchi K (1984) Site-specific recognition of the bacteriophage Mu ends by the Mu A protein. Cell 39: 387-394.

Bacteriophage Mu genome sequence: analysis and comparison with Mu-like prophages in Haemophilus, Neisseria and Deinococcus. G J Morgan, G F Hatfull, S Casjens, R W Hendrix, J Mol Biol. 317Morgan GJ, Hatfull GF, Casjens S, Hendrix RW (2002) Bacteriophage Mu genome sequence: analysis and comparison with Mu-like prophages in Haemophilus, Neisseria and Deinococcus. J Mol Biol 317: 337-359.

A-tract clusters may facilitate DNA packaging in bacterial nucleoid. M Y Tolstorukov, K M Virnik, S Adhya, V B Zhurkin, Nucleic Acids Res. 33Tolstorukov MY, Virnik KM, Adhya S, Zhurkin VB (2005) A-tract clusters may facilitate DNA packaging in bacterial nucleoid. Nucleic Acids Res 33: 3907- 3918.

Primary structure of phage Mu transposase: homology to Mu repressor. R M Harshey, E D Getzoff, D L Baldwin, J L Miller, G Chaconas, Proc Natl Acad Sci U S A. 82Harshey RM, Getzoff ED, Baldwin DL, Miller JL, Chaconas G (1985) Primary structure of phage Mu transposase: homology to Mu repressor. Proc Natl Acad Sci U S A 82: 7676-7680.

DNA sequence of the control region of phage D108: the N-terminal amino acid sequences of repressor and transposase are similar both in phage D108 and in its relative, phage Mu. M Mizuuchi, R A Weisberg, K Mizuuchi, Nucleic Acids Res. 14Mizuuchi M, Weisberg RA, Mizuuchi K (1986) DNA sequence of the control region of phage D108: the N-terminal amino acid sequences of repressor and transposase are similar both in phage D108 and in its relative, phage Mu. Nucleic Acids Res 14: 3813-3825.

Visualizing the assembly and disassembly mechanisms of the MuB transposition targeting complex. E C Greene, K Mizuuchi, J Biol Chem. 279Greene EC, Mizuuchi K (2004) Visualizing the assembly and disassembly mechanisms of the MuB transposition targeting complex. J Biol Chem 279: 16736-16743.

Congruence of in vivo and in vitro insertion patterns in hot E. coli gene targets of transposable element Mu: opposing roles of MuB in target capture and integration. J Ge, R M Harshey, J Mol Biol. 380Ge J, Harshey RM (2008) Congruence of in vivo and in vitro insertion patterns in hot E. coli gene targets of transposable element Mu: opposing roles of MuB in target capture and integration. J Mol Biol 380: 598-607.

Transcription induces a supercoil domain barrier in bacteriophage Mu. K E Scheirer, N P Higgins, Biochimie. 83Scheirer KE, Higgins NP (2001) Transcription induces a supercoil domain barrier in bacteriophage Mu. Biochimie 83: 155-159.

Growth rate toxicity phenotypes and homeostatic supercoil control differentiate Escherichia coli from Salmonella enterica serovar Typhimurium. K Champion, N P Higgins, J Bacteriol. 189Champion K, Higgins NP (2007) Growth rate toxicity phenotypes and homeostatic supercoil control differentiate Escherichia coli from Salmonella enterica serovar Typhimurium. J Bacteriol 189: 5839-5849.

Growth phasedependent variation in protein composition of the Escherichia coli nucleoid. Ali Azam, T Iwata, A Nishimura, A Ueda, S Ishihama, A , J Bacteriol. 181Ali Azam T, Iwata A, Nishimura A, Ueda S, Ishihama A (1999) Growth phase- dependent variation in protein composition of the Escherichia coli nucleoid. J Bacteriol 181: 6361-6370.

Inhibition of bacteriophage Mu transposition by Mu repressor and Fis. C M Van Drunen, C Van Zuylen, E J Mientjes, N Goosen, P Van De Putte, Mol Microbiol. 10van Drunen CM, van Zuylen C, Mientjes EJ, Goosen N, van de Putte P (1993) Inhibition of bacteriophage Mu transposition by Mu repressor and Fis. Mol Microbiol 10: 293-298.

Involvement of Escherichia coli FIS protein in maintenance of bacteriophage Mu lysogeny by the repressor: control of early transcription and inhibition of transposition. M Betermier, I Poquet, R Alazard, M Chandler, J Bacteriol. 175Betermier M, Poquet I, Alazard R, Chandler M (1993) Involvement of Escherichia coli FIS protein in maintenance of bacteriophage Mu lysogeny by the repressor: control of early transcription and inhibition of transposition. J Bacteriol 175: 3798-3811.

A defined system for the DNA strand-transfer reaction at the initiation of bacteriophage Mu transposition: protein and DNA substrate requirements. R Craigie, D J Arndt-Jovin, K Mizuuchi, Proc Natl Acad Sci U S A. 82Craigie R, Arndt-Jovin DJ, Mizuuchi K (1985) A defined system for the DNA strand-transfer reaction at the initiation of bacteriophage Mu transposition: protein and DNA substrate requirements. Proc Natl Acad Sci U S A 82: 7570- 7574.

Site-specific HU binding in the Mu transpososome: conversion of a sequence-independent DNA-binding protein into a chemical nuclease. B D Lavoie, G Chaconas, Genes Dev. 7Lavoie BD, Chaconas G (1993) Site-specific HU binding in the Mu transpososome: conversion of a sequence-independent DNA-binding protein into a chemical nuclease. Genes Dev 7: 2510-2519.

At the heart of the chromosome: SMC proteins in action. T Hirano, Nat Rev Mol Cell Biol. 7Hirano T (2006) At the heart of the chromosome: SMC proteins in action. Nat Rev Mol Cell Biol 7: 311-322.

A unique right end-enhancer complex precedes synapsis of Mu ends: the enhancer is sequestered within the transpososome throughout transposition. S Pathania, M Jayaram, R M Harshey, EMBO J. 22Pathania S, Jayaram M, Harshey RM (2003) A unique right end-enhancer complex precedes synapsis of Mu ends: the enhancer is sequestered within the transpososome throughout transposition. EMBO J 22: 3725-3736.

Role of DNA topology in Mu transposition: mechanism of sensing the relative orientation of two DNA segments. R Craigie, K Mizuuchi, Cell. 45Craigie R, Mizuuchi K (1986) Role of DNA topology in Mu transposition: mechanism of sensing the relative orientation of two DNA segments. Cell 45: 793-800.

Prophages and bacterial genomics: what have we learned so far?. S Casjens, Mol Microbiol. 49Casjens S (2003) Prophages and bacterial genomics: what have we learned so far? Mol Microbiol 49: 277-300.

The impact of prophages on bacterial chromosomes. C Canchaya, G Fournous, H Brussow, Mol Microbiol. 53Canchaya C, Fournous G, Brussow H (2004) The impact of prophages on bacterial chromosomes. Mol Microbiol 53: 9-18.

Horizontally acquired glycosyltransferase operons drive Salmonellae lipopolysaccharide diversity. M R Davies, S E Broadbent, S R Harris, N R Thomson, M W Van Der Woude, PLoS Genet. 91003568Davies MR, Broadbent SE, Harris SR, Thomson NR, van der Woude MW (2013) Horizontally acquired glycosyltransferase operons drive Salmonellae lipopolysaccharide diversity. PLoS Genet 9: e1003568.

Mobile genetic elements: the agents of open source evolution. L S Frost, R Leplae, A O Summers, A Toussaint, Nature reviews Microbiology. 3Frost LS, Leplae R, Summers AO, Toussaint A (2005) Mobile genetic elements: the agents of open source evolution. Nature reviews Microbiology 3: 722-732.

Increased reproductive fitness of Escherichia coli lambda lysogens. L Lin, R Bitner, G Edlin, J Virol. 21Lin L, Bitner R, Edlin G (1977) Increased reproductive fitness of Escherichia coli lambda lysogens. J Virol 21: 554-559.

A bacterial virulence determinant encoded by lysogenic coliphage lambda. J J Barondess, J Beckwith, Nature. 346Barondess JJ, Beckwith J (1990) A bacterial virulence determinant encoded by lysogenic coliphage lambda. Nature 346: 871-874.

A case of lysogenic conversion: modification of cell phenotype by constitutive expression of the Mu gem operon. L Paolozzi, P Ghelardini, Res Microbiol. 143Paolozzi L, Ghelardini P (1992) A case of lysogenic conversion: modification of cell phenotype by constitutive expression of the Mu gem operon. Res Microbiol 143: 237-243.

The Gam protein of bacteriophage Mu is an orthologue of eukaryotic Ku. F Di Fagagna, G R Weller, A J Doherty, S P Jackson, EMBO Rep. 4d'Adda di Fagagna F, Weller GR, Doherty AJ, Jackson SP (2003) The Gam protein of bacteriophage Mu is an orthologue of eukaryotic Ku. EMBO Rep 4: 47-52.

Reproductive fitness of P1, P2, and Mu lysogens of Escherichia coli. G Edlin, L Lin, R Bitner, J Virol. 21Edlin G, Lin L, Bitner R (1977) Reproductive fitness of P1, P2, and Mu lysogens of Escherichia coli. J Virol 21: 560-564.

One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. K A Datsenko, B L Wanner, Proc Natl Acad Sci U S A. 97Datsenko KA, Wanner BL (2000) One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci U S A 97: 6640- 6645.

Recombineering: in vivo genetic engineering in E. coli, S. enterica, and beyond. J A Sawitzke, L C Thomason, N Costantino, M Bubunenko, S Datta, Materials and Methods Enzymol. 421Sawitzke JA, Thomason LC, Costantino N, Bubunenko M, Datta S, et al. (2007) Recombineering: in vivo genetic engineering in E. coli, S. enterica, and beyond. Materials and Methods Enzymol 421: 171-199.

Identification and analysis of recombineering functions from Gram-negative and Gram-positive bacteria and their phages. S Datta, N Costantino, X Zhou, D L Court, Proc Natl Acad Sci U S A. 105Datta S, Costantino N, Zhou X, Court DL (2008) Identification and analysis of recombineering functions from Gram-negative and Gram-positive bacteria and their phages. Proc Natl Acad Sci U S A 105: 1626-1631.

Rare codons play a positive role in the expression of the stationary phase sigma factor RpoS (sigma(S)) in Escherichia coli. T Kolmsee, R Hengge, RNA Biol. 8Kolmsee T, Hengge R (2011) Rare codons play a positive role in the expression of the stationary phase sigma factor RpoS (sigma(S)) in Escherichia coli. RNA Biol 8: 913-921.

The gyrB gene product functions in both initiation and chain polymerization of Escherichia coli chromosome replication: suppression of the initiation deficiency in gyrB-ts mutants by a class of rpoB mutations. M Filutowicz, P Jonczyk, Mol Gen Genet. 191Filutowicz M, Jonczyk P (1983) The gyrB gene product functions in both initiation and chain polymerization of Escherichia coli chromosome replication: suppression of the initiation deficiency in gyrB-ts mutants by a class of rpoB mutations. Mol Gen Genet 191: 282-287.

Complete amino-acid sequences of DNA-binding proteins HU-1 and HU-2 from Escherichia coli. B Laine, D Kmiecik, P Sautiere, G Biserte, M Cohen-Solal, Eur J Biochem. 103Laine B, Kmiecik D, Sautiere P, Biserte G, Cohen-Solal M (1980) Complete amino-acid sequences of DNA-binding proteins HU-1 and HU-2 from Escherichia coli. Eur J Biochem 103: 447-461.

Integration host factor: a protein for all reasons. D I Friedman, Cell. 55Friedman DI (1988) Integration host factor: a protein for all reasons. Cell 55: 545-554.

Toxic protein expression in Escherichia coli using a rhamnose-based tightly regulated and tunable promoter system. M J Giacalone, A M Gentile, B T Lovitt, N L Berkley, C W Gunderson, Biotechniques. 40Giacalone MJ, Gentile AM, Lovitt BT, Berkley NL, Gunderson CW, et al. (2006) Toxic protein expression in Escherichia coli using a rhamnose-based tightly regulated and tunable promoter system. Biotechniques 40: 355-364.

Tight regulation, modulation, and high-level expression by vectors containing the arabinose PBAD promoter. L M Guzman, D Belin, M J Carson, J Beckwith, J Bacteriol. 177Guzman LM, Belin D, Carson MJ, Beckwith J (1995) Tight regulation, modulation, and high-level expression by vectors containing the arabinose PBAD promoter. J Bacteriol 177: 4121-4130.

Quantitative analysis of chromosome conformation capture assays (3C-qPCR). H Hagege, P Klous, C Braem, E Splinter, J Dekker, Nature protocols. 2Hagege H, Klous P, Braem C, Splinter E, Dekker J, et al. (2007) Quantitative analysis of chromosome conformation capture assays (3C-qPCR). Nature protocols 2: 1722-1733.

High-resolution circular chromosome conformation capture assay. A Gondor, C Rougier, R Ohlsson, Nature protocols. 3Gondor A, Rougier C, Ohlsson R (2008) High-resolution circular chromosome conformation capture assay. Nature protocols 3: 303-313.

Reactions of Cre with methylphosphonate DNA: similarities and contrasts with Flp and vaccinia topoisomerase. C H Ma, A H Kachroo, A Macieszak, T Y Chen, P Guga, PloS one. 47248Ma CH, Kachroo AH, Macieszak A, Chen TY, Guga P, et al. (2009) Reactions of Cre with methylphosphonate DNA: similarities and contrasts with Flp and vaccinia topoisomerase. PloS one 4: e7248.

. N Symonds, A Toussaint, P Van De Putte, M M Howe, Phage Mu. Cold Spring Harbor, New York: Cold Spring Harbor LaboratorySymonds N, Toussaint A, Van de Putte P, Howe MM (1987) Phage Mu. Cold Spring Harbor, New York: Cold Spring Harbor Laboratory.

The early promoter of bacteriophage Mu: definition of the site of transcript initiation. H M Krause, M R Rothwell, N P Higgins, Nucleic Acids Res. 11Krause HM, Rothwell MR, Higgins NP (1983) The early promoter of bacteriophage Mu: definition of the site of transcript initiation. Nucleic Acids Res 11: 5483-5495.