- Research article
- Open Access
High-throughput sequencing of sorted expression libraries reveals inhibitors of bacterial cell division
BMC Genomics volume 19, Article number: 781 (2018)
Bacterial filamentation occurs when rod-shaped bacteria grow without dividing. To identify genetically encoded inhibitors of division that promote filamentation, we used cell sorting flow cytometry to enrich filamentous clones from an inducible expression library, and then identified the cloned DNA with high-throughput DNA sequencing. We applied the method to an expression library made from fragmented genomic DNA of uropathogenic E. coli UTI89, which undergoes extensive reversible filamentation in urinary tract infections and might encode additional regulators of division.
We identified 55 genomic regions that reproducibly caused filamentation when expressed from the plasmid vector, and then further localized the cause of filamentation in several of these to specific genes or sub-fragments. Many of the identified genomic fragments encode genes that are known to participate in cell division or its regulation, and others may play previously-unknown roles. Some of the prophage genes identified were previously implicated in cell division arrest. A number of the other fragments encoded potential short transcripts or peptides.
The results provided evidence of potential new links between cell division and distinct cellular processes including central carbon metabolism and gene regulation. Candidate regulators of the UTI-associated filamentation response or others were identified amongst the results. In addition, some genomic fragments that caused filamentation may not have evolved to control cell division, but may have applications as artificial inhibitors. Our approach offers the opportunity to carry out in depth surveys of diverse DNA libraries to identify new genes or sequences encoding the capacity to inhibit division and cause filamentation.
Bacterial cell division genes were first discovered in filamentous, temperature-sensitive (fts) mutants . Over 50 years later, it is expected that many more genes involved in the mechanism and regulation of division remain to be discovered [2, 3]. In addition to identifying these, it is of great interest to identify inhibitors of cell division, as potential anti-microbial agents or research tools [4, 5]. Filamentation is also observed in nature as a conditional response. For example, in uropathogenic Escherichia coli (UPEC), extensive, reversible filamentation occurs during the late stage of bladder cell infection in a UTI .
Bacterial cell division starts with the identification of an appropriate mid-cell site and the assembly of division proteins into a ring-shaped structure around the cell, at the inner surface of the cell envelope. This develops into a structure known as the divisome, which constricts the cell and builds new cell poles for the two daughter cells . Most of the known regulators act at the initial assembly stage, although subsequent cell constriction and closure by the mature divisome must also be coordinated with chromosome replication and segregation.
In the E. coli cell cycle, known regulators of division include the nucleoid occlusion, Min, and Ter-linkage systems . Nucleoid occlusion inhibits cell division in areas where the bacterial chromosomes (nucleoids) are located. This involves the binding of the SlmA protein to specific sites in the chromosome, where it inhibits polymerisation and assembly of the tubulin-like division protein, FtsZ . The MinCD proteins concentrate on the inner membrane towards the poles of the cell, inhibiting divisome assembly there . Overproduction of SlmA-DNA or MinCD results in filamentation [8, 10]. The Ter-linkage system provides positive regulation of division. Interactions between MatP, bound specifically to sites in the chromosome replication terminus region, and the early-stage division protein ZapB, promote divisome assembly at mid-cell . This appears to help coordinate chromosome segregation and cell division [11, 12]. It has been recognised that as yet unknown cell cycle regulators of division are also likely to be present in E. coli [11, 13].
When compared to the processes involved in housekeeping cell cycle regulation of division described above, fewer regulatory pathways that lead to arrest of cell division and filamentation under specific conditions have been characterised. The best-known case is the DNA damage (or “SOS”) response . DNA damage is sensed by the RecA protein, which then triggers auto-degradation of the SOS repressor LexA, and expression of SOS genes. One of the SOS genes encodes the FtsZ assembly inhibitor, SulA. The resulting filamentation allows time for DNA repair; when the SOS response is switched back off, the filamentous cells can resume division to generate short rods with intact chromosomes .
Another important example of conditional filamentation is seen during urinary tract infections (UTI). During a late stage of the intracellular infection cycle, UPEC grow into extensive filaments that may be hundreds of microns long and erupt from the overwhelmed host cell . UPEC filaments are thought to help avoid phagocytosis by macrophages during infection and may improve bacterial dispersal and surface adhesion . It has been suggested that UPEC filamentation in UTI occurs via a SulA-dependent pathway . However, extensive SulA-independent filamentation by UPEC in experimental infection was recently observed, which instead required a known cell division protein, DamX . DamX normally contributes to cell envelope synthesis during division , and was previously noted to cause extensive filamentation when overproduced [20, 21]. The exact roles of these proteins in UPEC need clarification, and the signalling pathway(s) that lead to filamentation during UTI remain to be identified.
Here, we report the development and utilization of a high-throughput DNA sequencing approach to identify genes, or genomic DNA fragments, that result in filamentation after induced over-expression from a plasmid-based DNA library. We applied this method to a UPEC strain capable of UTI-associated filamentation. This identified new and known inhibitors or mediators of E. coli cell division as well as potential components of the UTI-associated or other filamentation responses.
Bacterial strains, plasmids and growth conditions
E. coli BW25113 and JW0941–1  were obtained from the E. coli Genetic Stock Center (Yale University, USA). E. coli UTI89  and plasmids pBAD24  and pLau80  were kindly provided by J. Moller-Jensen. Strains were routinely cultured at 37 °C on solid or in liquid Luria Broth (LB), with 1% (w/v) NaCl.
For the cloning of individual ORFs, DNA fragments containing fsaA, ybiY, aceE, pdhR, sdaA, lipA, pptE and pflC were amplified from E. coli UTI89 genomic DNA (gDNA) with Phusion DNA polymerase (New England Biolabs - NEB). DNA oligonucleotide sequences and the integrated restriction sites are shown in Additional file 1: Table S2. The products were purified using the Isolate II PCR and gel-extraction kit (Bioline) before digestion and cloning into pBAD24, utilizing the restriction enzymes indicated in Additional file 1: Table S2. The ligation mixtures were transformed into electro-competent E. coli BW25113 cells, and clones selected on solid LB agar supplemented with ampicillin (100 μg/ml). Ampicillin resistant colonies were selected, the expected plasmid inserts were confirmed by Sanger single-read sequencing (Australian Genome Research Facility).
Construction of the shotgun UTI89 genomic DNA expression library
Purified genomic DNA from E. coli UTI89 was partially digested with FatI (New England Biolabs), a 4-bp-recognition endonuclease that includes the ATG start codon, which would increase the frequency of in-frame ORF ligation to the expression vector compared to other 4-bp-recognition endonucleases. The gDNA fragments (1–5 kb range) were gel purified and ligated to pBAD24, which had been prepared by digestion with NcoI (NEB), dephosphorylated with Antarctic Phosphatase (NEB) and verified not to self-ligate by transformation of a vector-only ligation mixture. The library ligation mixture was used to transform E. coli JW0941–1 (ΔsulA) by electroporation, and clones were selected on LB agar supplemented with 100 μg/ml Ampicillin. Colonies arising after ~ 18 h growth at 37 °C were counted and then pooled by suspension and mixing in LB + 16% (v/v) glycerol, and stored in aliquots at − 80 °C.
Cell sorting for enrichment of filamentous clones from the expression library
Mid-log cultures of E. coli JW0941–1/pBAD24 (vector only), JW0941–1/pLau80 (encoding FtsZ-YFP under PBAD control) and the UTI89 gDNA expression library were diluted to A600 = 0.05 and induced with 0.2% (w/v) L-arabinose. Incubation was continued until A600 = 0.80. Culture samples were placed on ice and then analysed and sorted with an Aria II flow cytometer (BD Biosciences), as described previously . Specifically, 500,000 events from the “filamentous” gate were initially collected (“yield sort”), and then the sample was re-sorted to obtain ~ 20,000 events of higher purity (“purity sort”). Cultures were analysed at an average of 37,000 events per second for the initial yield sort, and at 14,000 events per second for the purity sort. Purity sorted cells and 1 mL of the unsorted sample were used to inoculate 5 mL LB media containing ampicillin and 0.2% (w/v) D-glucose and grown overnight at 37 °C with shaking (~ 150 rpm). Plasmid DNA from both cultures was extracted and purified using the Isolate II Plasmid Mini Kit (Bioline).
High-throughput multiplex DNA sequencing and data analysis
Plasmid DNA (5 ng) was prepared for sequencing using the Nextera DNA Library Prep Kit (Illumina), according to the manufacturer’s instructions. After the Nextera tagmentation and PCR amplification (tagging and barcoding) procedures, the DNA fragments were quantified on an Agilent Bioanalyzer and then underwent SPRI-select magnetic bead clean-up (Beckman-Coulter) at 0.8–0.5x left and right ratios (to select DNA of ~ 150–800 bp), according to the manufacturer’s instructions. The sorted and unsorted libraries were pooled, and quantified with an Agilent Bioanalyzer. The DNA was then diluted to 20 pM, including 5% PhiX genomic adapter-ligated DNA control library (Illumina). Sequencing was performed with an Illumina V3 150-cycle paired-end (2 × 75 cycles) flow cell and MiSeq instrument, according to the manufacturer’s instructions.
Sequence reads from the primary FASTQ files were demultiplexed according to their sample barcodes and then underwent quality control analysis and filtering using the FastQC software tool (Andrews, S. 2010), version 0.11.2. ( http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ ). Sequence read mapping to the UTI89 genome (Genbank NC_007946.1)  was done using Bowtie2 . Reads were first strictly mapped (zero mismatches allowed) to the pBAD24 plasmid. Unaligned reads were then strictly mapped to the plasmid pUTI89. Finally, the unaligned reads remaining were mapped strictly to the chromosome of UTI89. Statistics of sequencing reads and the data in Table 1 were generated by using the SAMtools software package and Flagstat command . The frequency of mapped reads and the analysis of data for Fig. 2 were determined with deepTools2 software . Artemis  was used to visualize the alignments of the reference genome with the read data (BAM files) for inspecting hit regions and annotating the output from the peak detection (described below).
Genomic regions causing filamentation were identified by comparing the read depth of the filamentous-sorted sample compared to the read-depth of reference-unsorted sample with MACS (version 2.1.1) peak-detection software . The alignment data of reads matched to the UTI89 genome were loaded into MACS using the Galaxy interface environment . The paired-end data for each sample were then combined, and the filamentous-sorted data were input into MACS as the test/tag file and the reference-unsorted data as the control file. Sequence reads of 75 bp and genome length of 5.06 Mbp were specified. The P-value threshold cutoff for peak detection and the bandwidth parameter were the default values of 0.00005 and 300, respectively. The ranges for calculating regional lambda were defined as 1000, 5000 and 10,000 bp in the control file. The value d, as determined by MACS, was 200 for both replicate screens.
Expression of cloned ORFs and analysis of cell phenotypes
Single colonies of BW25113-based strains, containing the cloned ORF’s in pBAD24, were used to inoculate liquid LB medium supplemented with ampicillin (100 μg/ml), and incubated overnight at 37 °C (~ 150 rpm). Cultures were then diluted in M9 minimal-media (supplemented with 0.4% (v/v) glycerol and 0.1% (w/v) casamino acids) to A600 = 0.02 and induced with the addition of 0.2% L-arabinose, and repressed with the addition of 0.2% (w/v) glucose, respectively. Cultures at A600 = 0.8 were fixed with volume of 2% (w/v) formaldehyde (prepared from a paraformaldehyde stock) and stored at 4 °C. Cell volume distributions were determined using a Multisizer 4 Coulter counter/cytometer (Beckman-Coulter), with a 50-μm aperture tube and Isoflow electrolyte (Beckman-Coulter) and calibrated with 2-μm latex beads. For microscopy, the cells were stained with FM4–64 (12.5 μg/mL) and Hoechst 33342 (0.5 μg/mL) and were then mounted on a 1% (w/v) agarose gel pad (in phosphate-buffered saline) and then imaged in phase-contrast and epifluorescence with a Zeiss Axioplan 2 microscope, with a 1.4 NA objective lens.
Amplification and sub-cloning of fragments from identified genomic regions
The two loci, ybiY-fsaA-moeB-moeA and ybeM-tatE-lipA were first amplified from UTI89 gDNA using Phusion High-Fidelity DNA polymerase (NEB). The amplicons were purified and then partially digested with FatI (NEB). This DNA was ligated to the compatible NcoI site of pBAD24 using T4 DNA ligase. The ligation mixture was transformed by electroporation into E. coli BW25113 and clones selected on solid LB agar supplemented with 100 μg/ml ampicillin. Forty-eight ampicillin resistant colonies were picked into 96-well plates containing LB or M9 media with 100 μg/ml ampicillin and 0.2% (w/v) L-arabinose. BW25113 containing pBAD24 only was also used as a control strain. Plates were incubated for 7 h at 37 °C (~ 150 rpm) and cells adhered to 15-well poly-L-lysine-coated glass slides and examined by phase-contrast microscopy.
A high-throughput screen based on flow cytometric sorting of DNA libraries
In order to identify genes or genomic regions that contain the potential to induce filamentation, a plasmid-based expression library of UPEC genomic DNA fragments was first constructed. Genomic DNA fragments of 1–5 kb from UPEC cystitis isolate UTI89 were ligated to the vector pBAD24, to place them under control of the arabinose-inducible promoter (PBAD). Approximately 35,000 colonies arising from transformation of E. coli JW0941–1 (K-12 background, ΔsulA) were resuspended and pooled for storage. The ΔsulA host strain was selected in order to identify mediators of filamentation that are independent of the well-characterised SulA division inhibitor.
Transcription from PBAD was induced by the addition of 0.2% (w/v) L-arabinose to a mid-log phase culture of the library, for 4 generations. These conditions were chosen to allow strong induction of transcription from PBAD, such that moderate or strong inhibitors of division can function, whilst providing sufficient time for these clones to generate filaments of lengths that were amenable to purification by flow cytometry . Flow cytometry cell sorting was then used to collect cells that had developed a filamentous morphology. Control cultures of JW0941–1/pBAD24 (“short” cells; Fig. 1a) and JW0941–1/pLau80 (“filamentous” cells, caused by over-production of FtsZ-YFP; Fig. 1b) were analysed too.
Approximately 20,000 cells that displayed a filamentous phenotype were isolated from the induced expression library (Fig. 1c and d, “filamentous” gate) on the basis of the side-scatter (pulse width) parameter, and they were then transferred to LB medium including glucose—to supress transcription from PBAD. These cells were then incubated overnight at 37 °C to allow reversal (division) of filaments and growth of viable bacteria (ongoing growth of filaments was expected to be eventually lethal). Plasmid DNA from this filamentous-sorted library was then purified and subjected to high-throughput sequencing, to identify the cloned UTI89 DNA that caused inducible filamentation. A reference sample from the same culture, treated in the same manner but not sorted, was also sequenced in the same multiplexed sequencing run. The whole experiment was performed in duplicate, starting with two frozen aliquots of the same library.
The mapped sequence-read frequencies for each sample showed that the expected majority of reads mapped to the vector pBAD24, with most of the remainder matching the UTI89 chromosome or native plasmid, pUTI89 (Table 1). The reads matching pBAD24 were excluded from further analysis. A plot of the reference library’s distribution of read depth per bp showed that the library covers 100% of the UTI89 chromosome by at least 1 read per bp, while < 1% was covered by at least 30 reads per bp (reference-unsorted 2 data, Additional file 2: Figure S1A).
We then detected genomic regions that were significantly enriched in reads from the filamentous-sorted samples, compared to the respective reference-unsorted samples. A global analysis of the degree of enrichment of specific genomic regions is shown in Fig. 2a; the filamentous-sorted sample showed a strong enrichment of reads in relatively fewer 10 kb genomic windows (or bins) compared to the reference-unsorted library. For example, 20% of the mapped reads from the filamentous-sorted sample were located in ~ 90% of the bins, while the remaining 80% of mapped reads were clustered in only 10% of the bins (dotted line, Fig. 2a). In contrast, 20% of mapped reads from the reference-unsorted sample were located in ~ 30% of bins, while the remaining 80% of mapped reads were distributed throughout ~ 70% of bins. This comparison demonstrates the effectiveness of the combination of cell sorting and high-throughput sequencing of plasmid DNA in enriching and identifying genomic regions. We characterised this enrichment further by plotting the number of mapped sequencing reads in 1 kb windows, to obtain an overview of the relative frequency of mapped reads along the entire chromosome. As may be seen in Fig. 2b, strong peaks of enrichment of read counts in the filamentous-sorted sample compared to the reference-unsorted sample were evident at loci distributed throughout the chromosome.
Identifying genomic regions that encode mediators of bacterial filamentation
We identified the significantly enriched regions in the filamentous-sorted data by using the peak-detection software MACS . The input data were the combined paired-end reads aligned to the UTI89 genome, from the filamentous-sorted and reference-unsorted samples. The aligned data (deposited in the European Nucleotide Archive under accession number PRJEB23331) were visualized with Artemis  to inspect regions identified by MACS. The complete annotated output of enriched regions from each replicate, and a list of the 55 regions that were found common to both replicate screens, are provided in Additional file 3: Table S1. We ranked these regions with respect to their significance score (−log10[P-value]), and found that the number of enriched regions identified by the two replicate screens corresponded strongly, above a significance score of ~ 70 (Additional file 2: Figure S1B). Furthermore, the fraction of enriched regions specifically in common between replicates was maintained near the maximum correspondence (~ 20%), at significance scores ≥70, but reduced markedly at lower significance scores (Fig. 2c). By applying this high-stringency significance threshold (≥70), 23 regions were identified in common between both replicate screens (Table 2). These regions are therefore most strongly implicated in causing filamentation.
The list in Table 2 includes regions that contain genes with a known involvement in cell division and its regulation, such as zapB  and minCDE. Figure 2d shows the data for the minCDE locus, encoding the well-known spatial regulator of division. The strongest region of relative enrichment coincides with the minC gene, which is the FtsZ inhibitor of the system . The region of enrichment extended downstream of minCDE, including a conserved gene of unknown function (C1358), and a pseudogene C1357 (the homolog of ycgI in E. coli K-12), suggesting that these genes or fragments in this region might affect cell division too. Alternatively, they may be represented due to their linkage to minCDE in our library. Some of the strongest hits were located in prophages or prophage remnants that have previously-demonstrated roles in blocking cell division, including kil1 and kil2 (λ-prophages), dicB/dicF (Qin prophage), and the ymfN-ymfK region of the CP-933O (e14) prophage. There were also numerous genes encoding membrane-associated proteins, such as ompA and rfaG (waaG), and others involved in peptidoglycan cell wall structure, including nlpD (Additional file 3: Table S1B). We also identified a region containing damX, which plays a role in cell division  and was previously shown to inhibit cell division when overexpressed, causing a filamentous phenotype . DamX has recently been specifically implicated in the conditional filamentation seen during urinary tract infections .
Amongst the expanded list of 55 identified regions that were common to both screens at the default threshold (Additional file 3: Table S1B) were other genes with a known involvement in cell division, such as zapA , suggesting that this list also contains significant inducers of filamentation. We expect that some regions only identified with strong significance in one of the screens may also contain genuine inducers of filamentation (Additional file 3: Table S1A), but they might have been missed in one of the screens due to limited sampling of cells from the library by the flow cytometry. Other regions encoding a capacity to cause filamentation may not have been cloned in the correct position with respect to PBAD in the library and were therefore not detected in either screen.
Identification and verification of fragments that cause filamentation
To investigate whether particular open reading frames (ORFs) from the identified DNA regions were responsible for the observed filamentation, we initially selected eight candidate ORFs from six of the hit regions of differing statistical strength, including some that were only observed in one screen and others that overlap the edge of a hit region (Additional file 3: Table S1A). These ORFs were individually cloned into pBAD24: fsaA, ybiY, aceE, pdhR, sdaA, lipA, pptE and pflC. Two of the cloned ORFs produced a reproducibly elongated or filamentous phenotype after the controlled induction period, namely pptE and pdhR, which showed a mean cell volume of 3.6 μm3 and 2.8 μm3, respectively, compared to 1.4 μm3 for the strain containing the vector-only (BW25113/pBAD24) (Fig. 3).
As may be seen in Fig. 3a and b, the cell size distributions of cultures overexpressing pptE and pdhR differed. Overexpression of pdhR, encoding the pyruvate dehydrogenase operon repressor, caused a mixed population of short cells and filaments. By comparison, overexpression of pptE, encoding a putative phosphoenolpyruvate-protein phosphotransferase, caused a much higher proportion of the cells to form filaments. Both genes caused their same characteristic filamentation in the absence of the SOS response regulator recA (Fig. 3b). Strikingly, staining of the membrane (FM4–64) and DNA (Hoechst 33342) in the pptE overexpression strain revealed irregularly positioned nucleoids in the filaments, and discontinuous zones of strong FM4–64 staining that corresponded to apparent inclusions in phase-contrast images (Fig. 3c). Similar inclusions were also present during conditions expected to result in significantly lower pptE expression (0.02% L-arabinose), although these conditions did not cause significant filamentation (Additional file 4: Figure S2). These observations suggest that the inclusions may be lipid-rich, and that filamentation occurs in this strain as a secondary effect, requiring higher levels of expression.
The six tested ORFs that did not show an obvious filamentous phenotype had average cell volumes of 1.3–1.6 μm3 (Additional file 5: Figure S3A). We considered it likely that other ORFs, a combination of multiple ORFs, a sub-fragment of an ORF or an expressed intragenic region might have caused the observed enrichment in the screens. We therefore further investigated two of the hit regions in which the preliminary whole-ORF testing above did not reveal filamentation, i.e. ybeM-tatE-lipA (2186 bp) and ybiY-fsaA-moeB-moeA (3754 bp). Sub-libraries of the two PCR-amplified regions were created by cloning FatI partial-digest products into pBAD24, and 48 transformants from each were cultured individually in M9 medium with L-arabinose for 7 h. Microscopy revealed a total of 12 clones from the two sets as filamentous, and the cloned fragments from these were sequenced (Fig. 4).
Seven of the filamentous clones contained fragments from within the ybeM ORF that were out-of-frame with respect to PBAD in the library. Two overlapping fragments of 112 bp and 87 bp were found in six of these clones (Fig. 4j); interestingly, both orientations of these fragments were observed amongst the clones, suggesting that two different products could trigger filamentation in this fragment. Alternatively, the region could conceivably contain its own Ara-inducible division inhibitor. The remaining clone (D10) contained a downstream fragment of only 29 bp within ybeM (encoding an amino acid sequence MFLQRRGAH).
From the ybiY-fsaA-moeB-moeA region, four of the five sub-clones identified as filamentous contained a 161 bp fragment in either orientation from within the ybiY ORF, despite the previous observation that the full-length ybiY ORF did not cause filamentation (Additional file 5: Figure S3A). Also, a larger downstream region of 2197 bp was identified (clone A11). This region spanned across the majority of the ybiY-fsaA-moeB-moeA region (Fig. 4k), with transcription in pBAD24 expected to start from within ybiY in the reverse orientation with respect to the expected transcription of ybiY in its genomic context. These findings show that sub-ORF in-frame fragments (e.g. clone C1 from ybiY) and out-of-frame fragments (e.g. clones H8, H9 and D10 from ybeM) can induce strong filamentation, even though the full-length ORFs (ybiY and ybeM) did not (Additional file 5: Figure S3B), suggesting that these fragments inhibit cell division via a mechanism that did not evolve for such a function in the whole ORFs they are located in.
We have developed a widely applicable high-throughput screen, combining flow cytometry cell sorting with high-throughput sequencing of a plasmid-based DNA expression library, to identify genes or genomic DNA fragments that act as inducible inhibitors of bacterial cell division and lead to bacterial filamentation. In this study, the filamentous clones isolated after expression of the genomic fragment library were recovered and amplified by overnight growth under conditions that supress expression. Since unlimited growth of bacterial filaments is expected to result in lysis, the recovery period was expected to enrich for reversible inhibitors of division, which were considered likely to be more physiologically relevant. We applied the method to E. coli UTI89, and identified 55 genomic regions that were enriched in the replicate screens. The identified genomic regions included numerous genes with known roles in cell division and its regulation, confirming that the method successfully identifies genes of biological relevance. Several previously uncharacterised genes were also identified that may represent new inhibitors or mediators of division, or be part of the UTI-associated or other filamentation responses. In addition, we identified a number of short predicted transcripts or peptides that are not predicted to be expressed from their genomic location in wild-type cells but were detected as sub-ORF, out-of-frame, or intergenic fragments. One fragment of only 29 bp, encoding a potential peptide (MFLQRRGAH), acts as an inhibitor of division in our screens. These or other similarly identified sequences could represent novel peptide or small RNA inhibitors of bacterial division. It is therefore expected that expanded screens of large and diverse DNA libraries could reveal other new inhibitors with potential applications.
An expectation of using a genomic fragment library for the screens described here is that not all genes will be appropriately positioned to be expressed from the plasmid. Furthermore, the size of our library (~ 35,000 clones) may also somewhat limit the screens, and they are therefore not considered fully comprehensive. Consistent with this, we identified six out of the 12 genomic regions previously identified by Sanger sequencing of clones obtained after a similar flow-cytometry screen was performed with E. coli K-12 (DH5α) . The genomic regions in common with this study contained genes such as damX (that has been implicated in the UTI-filamentation response pathway ), rplL-rpoB, and two regions of prophage elements e14 (ymfL/ymfM) and Rac (kil). Another possible limitation relates to the yields obtained from the cell sorter flow cytometer (~ 20,000 cells per screen). Consistent with such a sampling limitation, a number of enriched regions were detected with statistical significance in only one of the two screens (Fig. 2c and Additional file 3: Table S1). Genes from some of these regions play known roles in cell division (e.g. dedD). We therefore recommend that other such regions are also considered as potential candidates for verification. These data may be directly assessed by referring to Additional file 3: Table S1A and the genomic alignments available from the European Nucleotide Archive (accession number PRJEB23331).
Given our identification of both physiological and likely non-physiological classes of genomic fragments that lead to filamentation, we suggest that further screens aiming to focus on physiologically-relevant gene identification could utilize comprehensive complementary DNA (cDNA) expression libraries, whereas others focusing on artificial peptide inhibitors could use diverse small fragment libraries (e.g. from environmental or randomized synthetic DNA sources). The cDNA screen results should still be followed up by determining whether the identified ORFs relate to a specific pathway in controlling wild-type cell division, or whether they cause filamentation via an artefactual consequence of overexpression (e.g. via substantive changes to envelope physiology that secondarily inhibits normal cell division). Both approaches would benefit from the development of the flow cytometry or other methods to provide a significantly greater yield of filamentous bacteria from the expression library.
An intriguing aspect of our screen results was the presence of genes within the regions causing inducible filamentation that have known roles in central carbon metabolism, including gluconeogenesis (zwf) (Table 2 and Additional file 3: Table S1). Glucose metabolism has previously been implicated in inhibiting cell division in Bacillus subtilis, via the glucosyltransferase, ugtP, which acts on FtsZ to delay cell division in the presence of UDP-glucose . A similar UDP-glucose-dependent mechanism has also been found in E. coli with the functional analog, opgH . Interestingly, our screens revealed a region containing a fragment of rssB and the whole downstream galU gene (Additional file 3: Table S1). GalU catalyzes the synthesis of UDP-glucose, which could conceivably stimulate the division-inhibition mediated by OpgH. However, previous overproduction of a GalU fusion protein did not cause filamentation . Preliminary hits from our screen such as the rssB-galU region need to be checked by defining the specific regions and conditions that cause the filamentation.
We chose several genomic regions for following up and verifying the results of the screens, some of which encoded potentially novel division-metabolism links. Of particular note was the identification of a genomic region in one of our replicate screens containing the pdhR and aceE genes, encoding the pyruvate dehydrogenase complex regulator and pyruvate dehydrogenase e1 subunit, respectively. We showed that over-expression of the PdhR ORF caused filamentation, independent of the SOS pathway (recA) (Fig. 3). The PdhR regulator is controlled by the concentration of pyruvate in the cell. When ample pyruvate is present, PdhR does not repress its target genes (including its own PDH operon), whereas in low-pyruvate conditions, PdhR specifically binds DNA and represses expression of numerous metabolic genes and the cell division gene cluster, dcw, which includes ftsZ . In this way, low pyruvate (nutrient poor) conditions might reduce the frequency of division, coupling it to the reduced cell growth rate. Our finding that pdhR overexpression inhibits division therefore lends support to this notion. In B. subtilis, pyruvate metabolism also plays an important role in regulating division; pyruvate dehydrogenase was found to affect the assembly of the FtsZ-ring in a pyruvate-dependent manner . It was suggested that this positively regulates FtsZ assembly and division in nutrient-rich conditions to appropriately couple division with growth rate. Therefore, this important final stage of glycolysis appears to be a nexus in the coupling of cell growth and division in these two highly divergent bacteria, but the underlying mechanisms that evolved to establish this regulation appear to be very different. A similar occurrence of convergent evolution was noted for the UDP-glucose-dependent mechanisms of division control, mentioned above .
The second whole-ORF that we discovered to encode a strong inducer of filamentation was pptE (Fig. 3). PptE is uncharacterised in UTI89, but is a putative phosphoenol-pyruvate-protein phosphotransferase, that is present in extraintestinal pathogenic E. coli (ExPEC) strains (phylogenetic group B2) including UPEC, avian pathogenic E. coli (APEC) and neonatal meningitis E. coli (NMEC). The pptE gene varies in size substantially amongst E. coli strains. A smaller remnant of pptE (0.34 kb out of the total 2.38 kb ORF) is evident across some UPEC strains , although UTI89 has the full length ORF. Overexpression of pptE led to major cell structure and morphology effects in E. coli (Fig. 3, Additional file 4: Figure S2). It is possible that pptE plays an additional uncharacterised role in UTI89 cell division or structure, but an understanding of its normal physiological role(s) awaits further characterization. The overexpression phenotypes of pptE and other genes identified here offer new leads towards defining their roles.
Some of the strongest and most frequent inducers of filamentation from our screens were regions that contained specific prophage genes, several of which have been previously studied as strong inhibitors of E. coli cell division. For example, the kil1 gene in a λ-prophage was identified in our screens, and this gene encodes a known inhibitor of cell division that acts on FtsZ assembly . Other prophage regions in the list include a large region of the Qin prophage, containing known division inhibitors dicB and dicF. The gene dicB encodes an activator of the endogenous minC division regulator/inhibitor, whereas dicF is a small RNA inhibitor of ftsZ gene expression [43, 44]. The ymfN-ymfK region of the e14 prophage was identified (Table 2), which was previously reported to contain an SOS-inducible cell division inhibitor, originally termed sfiC . Furthermore, known prophage toxins from toxin-antitoxin pairs affecting cell division, including cbtA (yeeV) and gef (hok) [46, 47], were identified. These common functions of the prophage-associated genes therefore appear to reflect a wide-spread bacteriophage strategy to block cell division during lysis, which likely aids bacteriophage proliferation or dissemination. At least some of these genes in remnant or defective prophages appear to have been co-opted by the host bacteria for regulatory purposes in stress responses [48, 49].
We applied our method to E. coli UTI89, as it might contain specific regulators of the UPEC filamentation response seen in UTIs that are not present in other strains of E. coli. Alternatively, the filamentation response in UTI may involve a more general response pathway that is activated strongly under the conditions bacteria experience during UTI, such as the SOS response that was first postulated for UTI-associated filamentation . We identified seven genomic regions in the reproducibly filamentous expression clones that did not contain direct homologs in E. coli K-12. These included regions containing prophage genes that are implicated in inhibition of cell division (Additional file 3: Table S1B), and some of these might have been co-opted by cells for regulatory purposes . At present it is not clear whether these or other more wide-spread genes are involved in the UPEC filamentation response. However, the relevance of our screen towards the goal of identifying the components that are potentially involved in the UPEC filamentation response is exemplified by our identification of damX, which was previously identified amongst a set of genes that were upregulated during the filamentation stage of infection . DamX has a known involvement in division, and its overexpression causes filamentation, particularly when the bacteria are surface bound [18, 21]. DamX was essential for the UPEC filamentation response observed during infection .
Other components of the regulatory pathway that leads to UPEC filamentation (via DamX) are still unknown. Interestingly, amongst our list of 55 UTI89 genomic regions causing inducible filamentation (Additional file 3: Table S1B), there were 16 ORFs that showed elevated expression during the filamentation phase of bladder cell infection . These genes encode the outer membrane protein, OmpA (which is also under strong positive selection in UPEC strains ) and other proteins involved in peptidoglycan cell wall structure (NlpD and PbpG). Several relatively poorly characterized genes, such as ygaV, were also detected in both studies. YgaV is a predicted transcriptional regulator that is part of a two-component regulator (YgaVP) responsive to tributyltin exposure, but is of unknown function . These genes were therefore recognized as candidates that might act as regulators in the UTI-filamentation response.
A high-throughput sequencing-based method for genome-wide identification of genes and DNA fragments that encode a capacity to induce filamentation in E. coli was successfully developed and applied. This revealed genes from several prophages, as well as loci that have known or novel roles in bacterial cell division or filamentation, and has helped refine a list of candidate genes that may be involved in the UPEC filamentation response pathway seen extensively in UTIs. Several short DNA fragments that cause filamentation were also identified and these may have potential uses as inhibitors of cell division. Our approach can be used to carry out comprehensive surveys of diverse DNA libraries, to identify new sequences encoding the capacity to arrest cell division and induce filamentation.
Uropathogenic E. coli
Urinary Tract Infection
Van De Putte P, Van D, Roersch A. The selection of mutants of Escherichia coli with impaired cell division at elevated temperature. Mutat Res. 1964;106:121–8.
Adams DW, Errington J. Bacterial cell division: assembly, maintenance and disassembly of the Z ring. Nat Rev Microbiol. 2009;7(9):642–53.
Haeusser DP, Margolin W. Splitsville: structural and functional insights into the dynamic bacterial Z ring. Nat Rev Microbiol. 2016;14(5):305–19.
Lock RL, Harry EJ. Cell-division inhibitors: new insights for future antibiotics. Nat Rev Drug Discov. 2008;7(4):324–38.
Ma S, Ma S. The development of FtsZ inhibitors as potential antibacterial agents. ChemMedChem. 2012;7(7):1161–72.
Justice SS, Hung C, Theriot JA, Fletcher DA, Anderson GG, Footer MJ, Hultgren SJ. Differentiation and developmental pathways of uropathogenic Escherichia coli in urinary tract pathogenesis. Proc Natl Acad Sci U S A. 2004;101(5):1333–8.
Mannik J, Bailey MW. Spatial coordination between chromosomes and cell division proteins in Escherichia coli. Front Microbiol. 2015;6:306.
Cho H, McManus HR, Dove SL, Bernhardt TG. Nucleoid occlusion factor SlmA is a DNA-activated FtsZ polymerization antagonist. Proc Natl Acad Sci U S A. 2011;108(9):3773–8.
Pichoff S, Lutkenhaus J. Escherichia coli division inhibitor MinCD blocks septation by preventing Z-ring formation. J Bacteriol. 2001;183(22):6630–5.
Bi E, Lutkenhaus J. Cell division inhibitors SulA and MinCD prevent formation of the FtsZ ring. J Bacteriol. 1993;175(4):1118–25.
Bailey MW, Bisicchia P, Warren BT, Sherratt DJ, Mannik J. Evidence for divisome localization mechanisms independent of the min system and SlmA in Escherichia coli. PLoS Genet. 2014;10(8):e1004504.
Espeli O, Borne R, Dupaigne P, Thiel A, Gigant E, Mercier R, Boccard F. A MatP-divisome interaction coordinates chromosome segregation with cell division in E. coli. EMBO J. 2012;31(14):3198–211.
Lutkenhaus J, Du S. E. coli cell cycle machinery. Subcell Biochem. 2017;84:27–65.
Radman M. SOS repair hypothesis: phenomenology of an inducible DNA repair which is accompanied by mutagenesis. Basic Life Sci. 1975;5A:355–67.
Justice SS, Hunstad DA, Cegelski L, Hultgren SJ. Morphological plasticity as a bacterial survival strategy. Nat Rev Microbiol. 2008;6(2):162–8.
Horvath DJ Jr, Li B, Casper T, Partida-Sanchez S, Hunstad DA, Hultgren SJ, Justice SS. Morphological plasticity promotes resistance to phagocyte killing of uropathogenic Escherichia coli. Microbes Infection. 2011;13(5):426–37.
Justice SS, Hunstad DA, Seed PC, Hultgren SJ. Filamentation by Escherichia coli subverts innate defenses during urinary tract infection. Proc Natl Acad Sci U S A. 2006;103(52):19884–9.
Khandige S, Asferg CA, Rasmussen KJ, Larsen MJ, Overgaard M, Andersen TE, Moller-Jensen J. DamX Controls Reversible Cell Morphology Switching in Uropathogenic Escherichia coli. MBio. 2016;7(4). https://doi.org/10.1128/mBio.00642-16.
Gerding MA, Liu B, Bendezu FO, Hale CA, Bernhardt TG, de Boer PA. Self-enhanced accumulation of FtsN at division sites and roles for other proteins with a SPOR domain (DamX, DedD, and RlpA) in Escherichia coli cell constriction. J Bacteriol. 2009;191(24):7383–401.
Burke C, Liu M, Britton W, Triccas JA, Thomas T, Smith AL, Allen S, Salomon R, Harry E. Harnessing single cell sorting to identify cell division genes and regulators in bacteria. PLoS One. 2013;8(4):e60964.
Lyngstadaas A, Lobner-Olesen A, Boye E. Characterization of three genes in the dam-containing operon of Escherichia coli. Mol Gen Genet. 1995;247(5):546–54.
Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, Tomita M, Wanner BL, Mori H. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006;2:2006 0008. https://doi.org/10.1038/msb4100050.
Mulvey MA, Schilling JD, Hultgren SJ. Establishment of a persistent Escherichia coli reservoir during the acute phase of a bladder infection. Infect Immun. 2001;69(7):4572–9.
Guzman LM, Belin D, Carson MJ, Beckwith J. Tight regulation, modulation, and high-level expression by vectors containing the arabinose PBAD promoter. J Bacteriol. 1995;177(14):4121–30.
Lau IF, Filipe SR, Soballe B, Okstad OA, Barre FX, Sherratt DJ. Spatial and temporal organization of replicating Escherichia coli chromosomes. Mol Microbiol. 2003;49(3):731–43.
Chen SL, Hung CS, Xu J, Reigstad CS, Magrini V, Sabo A, Blasiar D, Bieri T, Meyer RR, Ozersky P, et al. Identification of genes subject to positive selection in uropathogenic strains of Escherichia coli: a comparative genomics approach. Proc Natl Acad Sci U S A. 2006;103(15):5977–82.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. Genome project data processing S: the sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Ramirez F, Ryan DP, Gruning B, Bhardwaj V, Kilpert F, Richter AS, Heyne S, Dundar F, Manke T. DeepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 2016;44(W1):W160–5.
Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16(10):944–5.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137.
Afgan E, Baker D, van den Beek M, Blankenberg D, Bouvier D, Cech M, Chilton J, Clements D, Coraor N, Eberhard C, et al. The galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update. Nucleic Acids Res. 2016;44(W1):W3–W10.
Ebersbach G, Galli E, Moller-Jensen J, Lowe J, Gerdes K. Novel coiled-coil cell division factor ZapB stimulates Z ring assembly and cell division. Mol Microbiol. 2008;68(3):720–35.
Hu Z, Mukherjee A, Pichoff S, Lutkenhaus J. The MinC component of the division site selection system in Escherichia coli interacts with FtsZ to prevent polymerization. Proc Natl Acad Sci U S A. 1999;96(26):14819–24.
Arends SJ, Williams K, Scott RJ, Rolong S, Popham DL, Weiss DS. Discovery and characterization of three new Escherichia coli septal ring proteins that contain a SPOR domain: DamX, DedD, and RlpA. J Bacteriol. 2010;192(1):242–55.
Galli E, Gerdes K. Spatial resolution of two bacterial cell division proteins: ZapA recruits ZapB to the inner face of the Z-ring. Mol Microbiol. 2010;76(6):1514–26.
Weart RB, Lee AH, Chien AC, Haeusser DP, Hill NS, Levin PA. A metabolic sensor governing cell size in bacteria. Cell. 2007;130(2):335–47.
Hill NS, Buske PJ, Shi Y, Levin PA. A moonlighting enzyme links Escherichia coli cell size with central metabolism. PLoS Genet. 2013;9(7):e1003663.
Gohler AK, Kokpinar O, Schmidt-Heck W, Geffers R, Guthke R, Rinas U, Schuster S, Jahreis K, Kaleta C. More than just a metabolic regulator--elucidation and validation of new targets of PdhR in Escherichia coli. BMC Syst Biol. 2011;5:197.
Monahan LG, Hajduk IV, Blaber SP, Charles IG, Harry EJ. Coordinating bacterial cell division with nutrient availability: a role for glycolysis. MBio. 2014;5(3):e00935–14.
Homeier T, Semmler T, Wieler LH, Ewers C. The GimA locus of extraintestinal pathogenic E. coli: does reductive evolution correlate with habitat and pathotype? PLoS One. 2010;5(5):e10877.
Haeusser DP, Hoashi M, Weaver A, Brown N, Pan J, Sawitzke JA, Thomason LC, Court DL, Margolin W. The Kil peptide of bacteriophage lambda blocks Escherichia coli cytokinesis via ZipA-dependent inhibition of FtsZ assembly. PLoS Genet. 2014;10(3):e1004217.
Balasubramanian D, Ragunathan PT, Fei J, Vanderpool CK. A prophage-encoded small RNA controls metabolism and cell division in Escherichia coli. mSystems. 2016;1(1):e00021–15.
Murashko ON, Lin-Chao S. Escherichia coli responds to environmental changes using enolasic degradosomes and stabilized DicF sRNA to alter cellular morphology. Proc Natl Acad Sci U S A. 2017;114(38):E8025–34.
D'Ari R, Huisman O. Novel mechanism of cell division inhibition associated with the SOS response in Escherichia coli. J Bacteriol. 1983;156(1):243–50.
Heller DM, Tavag M, Hochschild A. CbtA toxin of Escherichia coli inhibits cell division and cell elongation via direct and independent interactions with FtsZ and MreB. PLoS Genet. 2017;13(9):e1007007.
Pedersen K, Gerdes K. Multiple hok genes on the chromosome of Escherichia coli. Mol Microbiol. 1999;32(5):1090–102.
Barshishat S, Elgrably-Weiss M, Edelstein J, Georg J, Govindarajan S, Haviv M, Wright PR, Hess WR, Altuvia S. OxyS small RNA induces cell cycle arrest to allow DNA damage repair. EMBO J. 2017;37:413–26.
Wang X, Kim Y, Ma Q, Hong SH, Pokusaeva K, Sturino JM, Wood TK. Cryptic prophages help bacteria cope with adverse environments. Nat Commun. 2010;1:147.
Gueune H, Durand MJ, Thouand G, DuBow MS. The ygaVP genes of Escherichia coli form a tributyltin-inducible operon. Appl Environ Microbiol. 2008;74(6):1954–8.
This study was funded by the University of Technology Sydney. Iain Duggin was supported during part of this study by an Australian Research Council (ARC) Future Fellowship (FT160100010). The funding sources had no input into the study design or analysis and interpretation of outcomes.
Availability of data and materials
Sequence data in FASTQ and BAM file formats were deposited in the European Nucleotide Archive (accession number PRJEB23331). Materials used in this study are available to researchers on request to the corresponding author.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S2. Primer sequences used for cloning the indicated ORFs into pBAD24. PCR products generated using these primers in the appropriate combinations were digested with the indicated restriction enzymes as appropriate for cloning. (XLSX 9 kb)
Figure S1. Additional analyses of DNA sequencing data. (A) Plot showing the global read coverage of the UTI89 chromosome of the UTI89 genomic library, based on the read data from the “reference-unsorted 2” sample (Table 1). The library effectively covers the complete genome of UTI89, at a depth of 1–30 reads per bp. (B) A comparison of the number of identified enriched regions in replicate screens 1 and 2 by the MACS peak detection software (see Materials and Methods). The number of identified enriched regions were ranked by their significance score (−log10[P-value]), and the significance scores plotted against raw numbers of the identified regions. The dotted line represents the threshold of ≥70 P-value as the high-stringency criterion used to generate Table 2. (PDF 176 kb)
Table S1. E. coli UTI-89 genomic regions identified in overexpression screens for mediators of filamentation. Table S1A contains the combined MACS output from the two replicate screens. Regions indicated in bold text represent those found in common between the replicates (with an overlap of 100 bp or more). Those underlined indicate the 23 high-stringency regions found in common between replicates (with an overlap of 100 bp or more) and also were above the P-value threshold of 70 (see also Table 2). The gene name(s) included were determined by annotating the region from start to end (bp) and listing all ORFs within or overlapping those coordinates. This therefore will include partial ORFs that may not directly cause filamentation; further mapping may be required in these cases. Table S1B shows only the regions identified in common between the two replicate screens (with overlap of 100 bp or greater), and identifies ORFs within these that were not found to have direct homologs in E. coli K-12. (XLSX 69 kb)
Figure S2. pptE over-expression causes major effects on E. coli cell structure and morphology. BW25113 + pBAD24/pptE was grown in M9 minimal medium and induced with L-arabinose at the indicated concentrations and then fixed and stained with Hoechst 33342 (DNA) and FM4–64 (membrane) stains and then visualised by fluorescence microscopy. (A) 0.2% L-arabinose. (B) 0.02% L-arabinose. The intracellular inclusion-like structures and FM4–64 staining abnormalities were observed at both induction concentrations, whereas significant filamentation was only observed at the higher L-arabinose concentration. (PDF 1184 kb)
Figure S3. ORFs fsaA, ybiY, aceE, lipA, sdaA, and pflC do not cause substantial filamentation when overexpressed. The ORFs were expressed from pBAD24 in BW25113 by induction with 0.2% L-arabinose in M9 medium. Cells were fixed at OD600 = 0.8. (A) Cell volume distributions of the cell populations were determined by Coulter cytometry. The pflC expression strain shows a mild cell division defect or delay. (B) BW25113 + pBAD24/ybeM and BW25113 + pBAD24/ybiY strains examined by phase-contrast microscopy. Expression of ybeM appears to cause a minor effect on cell length, compared to the control and ybiY expression strains. (PDF 1231 kb)
About this article
Cite this article
Mediati, D.G., Burke, C.M., Ansari, S. et al. High-throughput sequencing of sorted expression libraries reveals inhibitors of bacterial cell division. BMC Genomics 19, 781 (2018). https://doi.org/10.1186/s12864-018-5187-7
- Cell division
- Uropathogenic E. coli
- Genetic screen
- DNA sequencing