INTRODUCTION
CRISPR-Cas system is a widespread adaptive immune system in prokaryotes including almost half of bacteria and most archaea [
1]. Short fragments of foreign DNA, or spacers, are integrated into Clustered, Regularly Interspaced Short Palindromic Repeat (CRISPR) loci on the host genome. The prokaryotes then apply the acquired spacers, which serve as memory elements, to target and cleave invading nucleic acids to defend against viruses and plasmids [
1–
5]. In the CRISPR loci, the spacers are separated by host-derived, short, repeated sequences, or repeats. Other important components of CRISPR system include CRISPR-associated (
cas) genes and a leader sequence, which are found near the CRISPR region. There are multiple kinds of
cas genes associated with a specific CRISPR [
2,
6]. During CRISPR action, the repeat-spacer array is first transcribed into a single transcript, which is further processed by Cas proteins into small RNAs (CRISPR RNA, or crRNA) [
4,
7,
8]. The mature crRNA and Cas proteins then assemble to form multicomponent CRISPR ribonucleoprotein (crRNP) complexes. The crRNP scans invading DNA by base pairing and on successful recognition, and then the target is eventually degraded by Cas nucleases [
6,
9].
According to the recent classification of CRISPR-Cas systems, there are two classes with fundamental differences in the organization of the effector module [
10]. Class 1 systems include types I, III and IV, and class 2 systems include types II, V, and VI [
1,
5]. Much of our understanding of type I system is gained from the studies of type I-E system of
E. coli K12 [
6,
11]. Recently, research has significantly advanced our knowledge about the less understood processes such as the acquisition of new spacers, for example, how the new spacers are produced and integrated into the CRISPR arrays [
5,
12–
14]. In addition, the mysteries of the structures and mechanisms of the target recognition have been largely uncovered recently [
3,
15–
19]. In the
E. coli type I-E system, there is a surveillance complex known as Cascade (CRISPR-associated complex for antiviral defense), a 405-kD complex of short crRNAs and five Cas proteins [
20]. Cascade binds to target DNA sequences and recruits Cas3 enzyme (a trans-acting nuclease-helicase, the signature Cas protein for type I system) to unwind and degrade the bound foreign DNA [
11,
16,
19]. Recently,
in vitro studies showed that linear DNA target is degraded in a unidirectional manner, and the degradation rate of negatively supercoiled DNA is 4.5 fold faster than the linear one [
11]. However, direct demonstration of DNA cleavage
in vivo has not been provided yet.
Here, we performed
in vivo experiments with the type I-E CRISPR system of
E. coli at the single-cell level. An artificial, plasmid-based CRISPR system was introduced into
E. coli to target bacteriophage lambda [
20], where we monitored CRISPR action by infecting cells with our phage strain. With our reporter system, we were able to track
in vivo lambda DNA degradation over time under the fluorescence microscope. This work provides insights on how CRISPR breaks down invading DNA.
RESULTS
Phage lysogenization is significantly reduced with the CRISPR system
For phage lambda, it was found that the
E. coli CRISPR system protects against lambda lysogenization, where the lysogenization frequency is 100-fold lower in the presence of the CRISPR compared to without [
21]. To test the efficiency of the artificial CRISPR system in our
E. coli phage DNA reporter host cell BA16 (or MG1655
seqA-yfpdam-) [
22], we utilized the same CRISPR plasmids: pWUR397A (
cas3) [
21], pWUR400 (
casABCDE) and pWUR478 (l phage targeting spacers) / pWUR477 (control spacers) [
20]. The CRISPR plasmid pWUR478 contains four anti-lambda spacers with homology to lambda genes
J, O, R, and
E. The control plasmid pWUR477 contains spacers without homology to lambda genome. To examine the function of this CRISPR system against the phage lysogenization, we performed experiments with the CRISPR or control plasmids at two different APIs (average phage input, i.e., average number of phages per cell) through a bulk lysogenization assay [
23]. Following the function of the CRISPR-Cas system, foreign phage DNA should be degraded and be unable to integrate into the host cell genome, resulting in reduced efficiency of lysogenization. As shown in Figure 1A, the lysogenization frequency of our host cells with the CRISPR or control plasmids is very low, within the order of 10
-7 – 10
-3 (the colony forming units, or cfu, and plaque forming units, or pfu, counts together with API and the calculated lysogenization frequencies are listed in Table S1). To determine the efficiency of this CRISPR-Cas system, we normalized the lysogenization frequency to that of the cells containing the control plasmid within two categories, i.e. low API (1–6) and high API (15 – 65). The lysogenization of the CRISPR targeting plasmid is only about 2%–7% as efficient compared to the control, as shown in Figure 1B. This reduced lysogenization frequency indicates that the artificial CRISPR-Cas system works at the bulk level. However, the 2%–7% efficiency in lysogenization is slightly higher than the 1% of previously reported
hns mutant strain harboring the same CRISPR plasmids at API of 1 [
21]. This could be due to our host strain MG1655 containing the
hns product, heat-stable nucleoid structuring (H-NS) protein, a global transcriptional repressor in
E. coli, or other host differences in our phage DNA reporter strain, for example, our cells with the
dam– allele lack host DNA methylation, thereby abrogating the regulation of
E. coli replication initiation by SeqA [
24–
27].
Visualizing CRISPR function in single cells
To visualize the phage DNA degradation process, we employed a reporter system for visualizing infecting phage particles on the cell surface, and phage DNA inside the cell [
23]. The infecting phage l
LZ760 is fluorescently labeled by the co-expression of gpD-mTurquoise2 and wild type gpD (the capsid decoration protein) on the phage capsid, and the packaged phage DNA is fully methylated [
28]. The host CRISPR/control strains (LZ1437/LZ1436) constitutively express a fluorescent SeqA fusion, SeqA-YFP, and the host DNA is not methylated owing to a
dam- mutation (methylation deficient) [
22]. SeqA specifically binds to fully methylated and hemi-methylated DNA, so the fully methylated DNA of phage l
LZ760 forms an YFP focus once ejected into the cell. This reporter system allows the tracking of each initial phage DNA and its first, hemi-methylated, replicated copy. However, the replicated DNAs are not always necessarily visualized as two separate fluorescent spots since the replicated DNAs can stick together and spatially overlap, or stay not far away from each other, thus appearing as one spot during imaging under the microscope.
To determine how efficiently CRISPR functions, we compared the cell fates after phage infection in both control movies (i.e., cells containing three plasmids: pWUR397A (
cas3), pWUR400 (
casABCDE) and pWUR477 (control spacers)) and CRISPR movies (i.e., cells containing three plasmids: pWUR397A (
cas3), pWUR400 (
casABCDE) and pWUR478 (l phage targeting spacers)). A successful phage infection is indicated by the appearance of a YFP focus during imaging. In the control movies (592 cells in 7 movies), 398 cells (67%) start with one YFP focus staying with one focus throughout the movies (120 minutes) (classified as group 1), 186 cells (31%) splitting into two foci (group 2), and 8 cells (2%) splitting into 3 or more probably due to recombination (group 3) [
23]. In the CRISPR movies (866 cells in 12 movies), there are 647 cells (75%) in group 1, 206 cells (24%) in group 2, and 13 cells (1%) in group 3. For the simplicity of spot tracking, we only analyzed cells with only one YFP focus throughout the movies (cells in group 1). Two typical cell fates are observed: lytic (86%, 341 out of 398 cells in the control movies; 66%, 425 out of 647 cells in the CRISPR movies) and CRISPR cells (26%, 167 out of 647 cells in the CRISPR movies). There are also rare cases of lysogeny or non-growing cells in both types of movies. The lytic cells show the accumulation of the fluorescent lytic reporter, gpD-mTurquoise2 protein inside the cell. The CRISPR cells are defined as those that have phage DNA with much shorter SeqA-YFP spot lifetime and do not lyse, instead show normal cell growth. Figure 2 shows the snapshots of representative cells over time. At 0 min, phage DNA of a fluorescently labeled phage (appears as a red focus on the cell surface) has been ejected into the cell forming a YFP focus (appears as a green focus). Over time, the fluorescence intensity of the phage DNA (green or yellow focus, yellow color shown in the overlay images is due to the overlay of green and red fluorescence) decreases. DNA foci then disappear at 110 min (the lytic cell in the control movie), 120 min (the lytic cell in the CRISPR movie) and 45 min (the CRISPR cell in the CRISPR movie). The lytic cells eventually lyse, releasing fluorescent phage particles outside the cell while the CRISPR cell divides and grows normally like the uninfected cell. The time-lapse movies are shown in Supplementary Movies S1‒ S3. Five more examples of each fate are shown in Supplementary Movies S4‒S6. As our movies only last 2 hours, some lytic cells do not lyse by the end of the movies, but the lytic reporter (gpD-mTurquoise2, appeared as red color) accumulates greatly indicating their lytic fate.
Factors associated with effective CRISPR degradation
To determine and quantify the degradation rate of DNA by CRISPR in each cell, we plot the spot intensity over time. The spot intensity in the cell population mostly maximizes within 5 minutes from the beginning of the movie indicating that phage ejection is fairly synchronized among cells, shown in Supplementary Figure S1. The small variation in the timing of maximum intensity is likely due to the variability of DNA ejection time and the time delay from setting up our microscopy movies [
23,
29]. In addition, the average maximum DNA focus intensities in lytic cells are similar to those in CRISPR cells in CRISPR movies, therefore all cells have similar levels of SeqA-YFP binding prior to CRISPR functioning independent of cell fate (Supplementary Figure S2). To simplify the analysis and eliminate the variance of the maximum intensities among cells, we normalized spot intensities at each time point by the maximum intensity of each trace. We find that the spot intensity in CRISPR cells decreases much faster than that in lytic cells, shown in Figure 3A, which indicates that CRISPR is functioning to degrade labeled phage DNA. The DNA foci become undetectable by 43.9±1.5 minutes to provide an estimated degradation rate (Figure 3B). Lytic cells also show fluorescence signal loss, but this is probably due to photobleaching, which provides a control for the effect of photobleaching under our experimental conditions. The spots in lytic cells usually disappear much later than those in CRISPR cells or the disappearance is not observed over the course of the movies (120 minutes) (Figure 3C). The spot disappearance in lytic cells is likely due to the phage DNA packaging, which would strip SeqA-YFP proteins off the phage DNA. As shown in Figure 3D, the spot disappearance time is correlated with cell lysis time with a correlation coefficient of 0.91, p-value of 0.01. When we compare the lytic cells in CRISPR movies with those in control movies, their fluorescence decreases similarly indicating that CRISPR spacers can fail in targeting the invading DNA, similarly to the control spacers (Supplementary Figure S3). Under our experimental conditions, the efficiency of CRISPR in targeting and degrading the phage DNA is around 26%, 167 CRISPR-active cells out of 647 cells with one focus throughout the movies.
We next examined various factors associated with the time of CRISPR spot disappearance corresponding to the degradation of invading phage DNA. Considering the initial maximum spot intensity, due to time variation of taking movies, both between experiments and different stages within the same movie (up to a few minutes), the maximum intensity captured for each cell trace may correspond to a partially degraded phage DNA. If this were the case, the maximum intensity would be expected to have a positive correlation with spot disappearance time. Indeed, when we plot spot maximum intensity against disappearance time, we find the correlation coefficient to be 0.97 with p-value of 0.03 (Figure 4A). Next, we examined the effect of initial cell size, as cell size (or volume, or length as proxy) influences protein concentration, here the Cas proteins. Different concentrations of Cas proteins may degrade phage DNA at different rates. From our experiments, it seems the phage degradation does not change with cell size (correlation coefficient of-0.55, with a large p-value of 0.45) as shown in Figure 4B. This indicates that the expressed Cas proteins are sufficiently concentrated to degrade phage DNA. Finally, we examined how phage DNA location inside the cell affects CRISPR action. Previous work has showed that phage lambda tends to adsorb to the cell pole and eject its DNA [
30–
32]. Although the exact mechanism underlying this polar preference has not been uncovered, cell pole is a special location where many proteins associated with the infection cycle, like ManYZ and FtsH accumulate [
30]. We then ask whether the degradation of phage DNA inside the cell has any correlation with the initial phage DNA location. After the phage is ejected inside the cell, our recent work has shown that some phage DNAs move locally while others move over the whole cell [
23]. When we correlate the phage DNA location with the efficiency of CRISPR degradation, we find that there appears to be no difference in CRISPR efficiency, whether the initial phage location is at polar/mid-cell (future pole) or non-polar positions (27%±4% at pole versus 25%±2% at non-pole, p= 0.02) (Figure 4C). Additionally, the timing of phage DNA degradation (spot disappearance time) is unaffected by phage DNA location (42.1±2.4 min at pole/mid-cell and 44.7±1.8 min at non-pole, p= 0.02) (Figure 4D). This indicates that the Cas proteins can effectively target invading DNA regardless of where the DNA is located.
DISCUSSION
In this study, we used our recently developed fluorescence technique to visualize phage DNA inside the cell upon infection by phage lambda to explore single-cell CRISPR function. With an artificial type I CRISPR system expressed in the host, the CRISPR system is able to degrade the invading phage DNA effectively with an efficiency of around 26% at the single-cell level. Consistent with the literature, the CRISPR system decreases the lysogenization frequency of the phage to be around 2%-7%, which shows the effectiveness of CRISPR system degrading the invading phage DNA. To our surprise, the fluorescent DNA foci in CRISPR cells do not divide into several smaller foci upon degradation, since the targeting plasmid pWUR478 contains four anti-lambda spacers with homology to lambda genes
J, O, R, and
E, which are fairly far apart from each other. If each spacer is targeted, the expected result would be that the phage genome is broken into multiple pieces, resulting in a few small fluorescent foci. However, our time-lapse movies show that the initial fluorescent DNA focus always stays as one before it finally disappears. Therefore, it is possible that only one anti-lambda spacer can target the phage DNA effectively one time with a synergistic effect of individual spacers and the Cas proteins degrade the DNA from that site thereafter, which may be due to the different sensitivities of the phage lambda DNA to the anti-lambda spacers [
20]. But we cannot rule out the possibility that our SeqA-YFP system holds the cut DNA together somehow such that there is only one visible spot.
Under our experimental conditions under the fluorescence microscope, the lambda DNA takes about 44 minutes to be degraded (Figure 3A and B). As a comparison, in
in vitro experiments, a 5 kb negatively supercoiled DNA and linear DNA can be degraded within a few minutes (~ 1 min and 5 min) [
11]. The much longer lambda DNA, 48.5 kb compared to this 5 kb DNA, almost 10-fold longer, might be the major cause why it takes almost 10-fold or more time to be digested. In addition, another contributor of different degrading time is the different experimental conditions of the
in vitro experiments versus our
in vivo experiments under the fluorescence microscope. The average fluorescence decrease in the lytic cells in the control and CRISPR movies (Supplementary Figure S4) is very similar, probably representing the photobleaching effect. The cell-to-cell variability of all the individual cells is probably due to the fluorescently labeled phage DNA moving in and out of focus, stochastic photobleaching of fluorescent molecules on the attached phage DNA and some unknown cell-to-cell variability. In addition to the stochastic photobleaching effect, the cell-to-cell variability of the CRISPR cells might also come from how easily the invading lambda DNA, probably the supercoiled form, is exposed for the Cas proteins to dock and degrade (Figure 3A).
From our single-cell analysis, the time to fully degrade the invading phage DNA is dependent on the initial fluorescent intensity of the phage DNA at the beginning of the movie. The lower maximum intensity of a spot in the cell might correspond to a partially degraded phage DNA due to the time to set up a movie, which is consistent with a fast initiation of CRISPR action. However, the phage DNA degradation does not depend on the initial cell size or phage DNA location. It is reasonable that Cas proteins are well mixed inside the cell and have high enough concentration to degrade the invading phage DNA with no preference of the location.
MATERIALS AND METHODS
Bacterial strains, plasmids and phages
The CRISPR/control strains (LZ1437/LZ1436) are BA16 (or MG1655
seqA-yfp dam-) [
22] harboring three plasmids pWUR397A (
cas3, AmpR) [
21], pWUR400 (
casABCDE, StrR) and pWUR478 (l phage targeting spacers,
CmR) / pWUR477 (control spacers,
CmR) [
20].
Fully methylated phage l
LZ760 was produced by inducing the lysogen MG1655(l
D-mturqoise2cI857bor::KanR)[pBR322-pLate-
D] harboring a Dam methylase overproducing plasmid pZA32-dam (with 1 mM IPTG induction) [
28]. The crude lysate was then purified from a standard CsCl purification procedure following the protocol described in [
33].
All bacterial strains, plasmids and phages are listed in Table 1.
Fluorescence microscopy and imaging
An overnight culture of the host CRISPR/control strains LZ1437/LZ1436 was diluted 1:100 in M9 minimal medium (11.3 g/l M9 salts, 1 mM MgSO4, 0.5 mg/mL thiamine HCl, 0.1% casamino acids, 100 mM CaCl2) supplemented with 0.4% maltose (M9M) at 37°C. 0.2% L-Arabinose and 0.1 mM IPTG were added into the culture to induce CRISPR plasmids for about 1.5 hours. When cells were grown to OD600 ≈ 0.3, and then concentrated and resuspended in ice-cold M9M with 0.2% L-Arabinose and 0.1 mM IPTG to OD600 ≈ 3. lLZ760 phages were added to reach a desired API (10 ml of cells+ 10 ml of phages), followed by incubation on ice for 30 min and an additional 5 min incubation at 35°C to trigger phage DNA ejection. The phage-cell mixture was diluted 1:10 in M9 with 0.2% L-Arabinose and 0.1 mM IPTG, and 1 ml of the diluted phage-cell mixture was placed on a thin 1.5% agarose slab of M9M with 0.2% L-Arabinose and 0.1 mM IPTG (~ 1 mm thick). After 1 min, a coverslip (No.1.5, Fisher Scientific) was gently laid over the mixture and the sample was imaged under the fluorescence microscope at 30°C with a cage incubator (InVivo Scientific, St. Louis, MO). Microscopy was performed on an inverted epifluorescence microscope (Ti-E, Nikon, Tokyo, Japan) using a 100× objective (Plan Fluo, NA 1.40, oil immersion) with a 2.5× TV relay lens and standard filter sets. Images were acquired using a cooled EMCCD camera (iXon3 897, Andor, Belfast, United Kingdom). Acquisition was performed using Nikon Elements software.
To localize all phages surrounding the cells, a series of 5 z-axis (vertical) images at a spacing of 400 nm were taken through the mTurquoise2/CFP channel using 200 ms exposure. To obtain more data in each time-lapse movie, cells were imaged at multiple stage positions (typically 16) in each experiment. During the time-lapse movie, the sample was imaged in phase contrast (100 ms exposure, for cell recognition), YFP (200 ms exposure, for phage DNA detection inside the cell, with a series of 5 z-axis images at a spacing of 400 nm in order to localize phage DNA foci), mTurquoise2/CFP (200 ms exposure, for phage lytic reporter) channels. The time-lapse movies were taken at a time interval of 5 min, until the cell fate was observable (2 hours).
Data analysis
All data analysis was performed in Matlab (The MathWorks, Natick, MA). Images were processed using MicrobeTracker [
34]. Briefly, cells were first outlined by the phase-contrast channel using MicrobeTracker, after which spots were recognized first automatically using SpotFinderZ, then manually corrected using SpotFinderM [
34]. Cell lineage tracking and fluorescent intensity were processed and calculated using our home-made Matlab script.
Bulk lysogenization assay
To measure the lysogenization frequency, we followed the protocol as described in [
32]. Briefly, 2 mL of the host CRISPR/control strains LZ1437/LZ1436 was grown in LBMM (LB+ 0.2% maltose+ 10 mM MgSO
4) with appropriate antibiotics for overnight and subsequently diluted 1:100 into 12 mL of LBMM and grown at 37
oC. 0.2% L-Arabinose and 0.1 mM IPTG were added into the culture to induce CRISPR plasmids for about 1.5 hours. When grown to OD
600 ≈ 0.4, cells were centrifuged (1000×
g for 10 min at 4
oC), concentrated 10x and resuspended to OD
600 ≈ 4 in pre-chilled LBMM with 0.2% L-Arabinose and 0.1 mM IPTG. 20 ml of the resuspended cells were then infected with 20 ml of phages at different concentrations by incubation for 30 min on ice. The samples were then transferred to 35
oC water bath for 5 min to allow for phage DNA ejection, followed by 100-fold dilution into pre-warmed LBGM (LB+ 0.2% glucose+ 10 mM MgSO
4) with 0.2% L-Arabinose and 0.1 mM IPTG and incubation with shaking at 265 rpm at 30
oC for 45 min. The samples were then properly diluted and plated on LB+Kn plates to allow around 100 colonies to grow.
Higher Education Press and Springer-Verlag Berlin Heidelberg