# Phage susceptibility testing and infectious titer determination through wide-field lensless monitoring of phage plaque growth

Prisca Perlemoine, Pierre R. Marcoux, Emmanuel Picard, Emmanuel Hadji, Marc Zelsmann, Grégoire Mugnier, Aurélie Marchet, Grégory Resch, Larry O'Connell, Eric Lacot

► To cite this version:

Prisca Perlemoine, Pierre R. Marcoux, Emmanuel Picard, Emmanuel Hadji, Marc Zelsmann, et al. Phage susceptibility testing and infectious titer determination through wide-field lensless monitoring of phage plaque growth. PLoS ONE, 2021, 16 (3), pp.e0248917. 10.1371/journal.pone.0248917. cea-03188740

HAL Id: cea-03188740

https://cea.hal.science/cea-03188740

Submitted on 2 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Phage susceptibility testing and infectious titer determination through wide-field lensless monitoring of phage plaque growth

Prisca Perlemoine<sup>1*</sup>, Pierre R. Marcoux<sup>1*</sup>, Emmanuel Picard<sup>2</sup>, Emmanuel Hadji<sup>2</sup>, Marc Zelsmann<sup>3</sup>, Grégoire Mugnier<sup>3</sup>, Aurélie Marchet<sup>4</sup>, Grégory Resch<sup>4</sup>, Larry O'Connell<sup>5</sup>, Eric Lacot<sup>6</sup>

1 Department of Microtechnologies for Biology and Health, LETI, CEA, University Grenoble Alpes, Grenoble, France, 2 SINAPS, PHELIQS, DEPHY, IRIG, DRF, CEA, University Grenoble Alpes, Grenoble, France, 3 LTM-Micro and Nanotechnologies for Health, CNRS, CEA, University Grenoble Alpes, Grenoble, France, 4 Department of Fundamental Microbiology, University of Lausanne, Lausanne, Switzerland, 5 SyMMES, IRIG, DRF, CEA, University Grenoble Alpes, Grenoble, France, 6 Laboratoire Interdisciplinaire de Physique, CNRS UMR 5588, University Grenoble Alpes, St Martin d'Hères, France

* pierre.marcoux@cea.fr (PRM); prisca.perlemoine@gmail.com (PP)

<figure><img src="image_2.png" /></figure>

## Abstract

The growing number of drug-resistant bacterial infections worldwide is driving renewed interest in phage therapy. Based on the use of a personalized cocktail composed of highly specific bacterial viruses, this therapy relies on a range of tests on agar media to determine the most active phage on a given bacterial target (phage susceptibility testing), or to isolate new lytic phages from an environmental sample (enrichment of phage banks). However, these culture-based techniques are still solely interpreted through direct visual detection of plaques. The main objective of this work is to investigate computer-assisted methods in order to ease and accelerate diagnosis in phage therapy but also to study phage plaque growth kinetics. For this purpose, we designed a custom wide-field lensless imaging device, which allows continuous monitoring over a very large area sensor (3.3 cm²). Here we report bacterial susceptibility to *Staphylococcus aureus* phage in 3 hr and estimation of infectious titer in 8 hr 20 min. These are much shorter time-to-results than the 12 to 24 hours traditionally needed, since naked eye observation and counting of phage plaques is still the most widely used technique for susceptibility testing prior to phage therapy. Moreover, the continuous monitoring of the samples enables the study of plaque growth kinetics, which enables a deeper understanding of the interaction between phage and bacteria. Finally, thanks to the 4.3 µm resolution, we detect phage-resistant bacterial microcolonies of *Klebsiella pneumoniae* inside the boundaries of phage plaques and thus show that our prototype is also a suitable device to track phage resistance. Lensless imaging is therefore an all-in-one method that could easily be implemented in cost-effective and compact devices in phage laboratories to help with phage therapy diagnosis.

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Perlemoine P, Marcoux PR, Picard E, Hadji E, Zelsmann M, Mugnier G, et al. (2021) Phage susceptibility testing and infectious titer determination through wide-field lensless monitoring of phage plaque growth. PLoS ONE 16(3): e0248917. https://doi.org/10.1371/journal.pone.0248917

**Editor:** Adelaide Almeida, Universidade de Aveiro, PORTUGAL

**Received:** October 5, 2020

**Accepted:** March 8, 2021

**Published:** March 23, 2021

**Copyright:** © 2021 Perlemoine et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files. The minimal dataset used to reach conclusions drawn in the manuscript is included in S2 Dataset. More over S1 Dataset and S1 Files provide images and image processing algorithm to allow the reader to reproduce the processing used in the manuscript. Python codes cannot be included due to intellectual property reason.
**Funding:** The author(s) received no specific funding for this work.

**Competing interests:** The authors have declared that no competing interests exist.

Introduction

The increasing number of drug-resistant bacterial infections is an emerging global health crisis that is driving a growing demand for phage therapy [1]. First implemented in 1919 [2], phage therapy is based on the use of highly specific bacterial viruses called bacteriophages or phages. Through a lytic life cycle, some phages are able to replicate within the cytoplasm of bacterial cells before being released through lysis of their host. In recent years, promising clinical studies of phage therapy were conducted to treat infection of burn wounds [3], urinary tract infections [4], and chronic otitis [5] caused by antibiotic-resistant bacteria. Currently, increasing evidence is leading to a consensus in the scientific community regarding synergism between phages and antibiotics [6, 7].

Regarding the high specificity of phages to particular bacterial species and even particular strains thereof, phage therapy is envisioned as part of a shift to a new paradigm of personalized medicine wherein personalized phage cocktails are administered combining only the most effective phages for a given bacterial isolate [8]. Accordingly, this approach requires access to large phage libraries (or phage banks) that must be maintained and constantly replenished with novel lytic phages targeting otherwise phage-resistant clones. Moreover, evolutionary selection of bacterial mutants resistant to the phages composing the cocktail should be taken into account, as much as possible, in the process of developing personalized phage cocktails and carefully monitored during treatment to avoid therapeutic failures [9]. Therefore, the success of personalized phage medicine largely relies on our capacity to quickly identify and isolate new lytic phages against a diverse range of pathogens [10] as well as to monitor the development of phage-resistance.

The classical method for phage discovery and isolation is to mix liquid cultures of clinical bacterial isolates with samples expected to contain phages, [11] e.g. human feces [12], chopped chicken intestines [13], raw sewage [14–16], farmyard slurries [17] and river waters [18]. After a suitable period and temperature of incubation, remaining host cells and debris are removed through centrifugation/filtration and the aqueous phase (i.e. the so-called lysate) is assayed for the presence of phages using the double agar overlay method [8]. At this stage, the lysate may contain several different phages that have been proliferated in, especially if samples contained a complex and diverse flora. For instance, Primrose *et al*. observed as many as seven different plaque morphologies in a double agar overlay assay using a single *Escherichia coli* HfrH strain as a host, and river samples collected downstream from the effluent discharge point of a sewage works as a source of novel phages [18]. To the best of our knowledge, despite the laborious nature of this process, direct visual evaluation of plaques by the experimenter is still the sole technique routinely used in laboratories to discriminate between phages. Usually this time-consuming visual inspection of phage plaque is only the first step of the phage characterization procedure and is combined with other evaluations, such as time-kill kinetics assay through kinetic optical density measurements [19]. Moreover, genome sequencing is performed on candidate viral strains for therapeutic use to exclude temperate phages and detect the presence of undesirable genes encoding antimicrobial resistance or virulence factors [20]. Accordingly, more automated approaches are highly desired, including tools that could quickly detect incipient plaques, accurately enumerate phage titer, study and discriminate between their morphologies and growth kinetics, and simultaneously detect phage-resistance selection. Indeed, such a tool would greatly help enrichment and screening of existing phage banks. While Koch [21] reported one of the first studies of plaque growth kinetics through direct visual measurement, other groups reported attempts to automate this assay through time-lapse photography [22–24]. These early attempts allowed the study of only up to six plaques simultaneously. In this
work, we evaluate lensless imaging for the high-throughput study of plaque growth kinetics as well as the detection of phage-resistance selection.

Lensless imaging is inspired by the in-line holography technique of D. Gabor [25]. It consists of imaging a sample placed directly on top of an imaging sensor without the use of any intervening optical objectives. Lensless imaging presents several advantages. Firstly, there are no optical aberrations caused by the use of lenses. Moreover, because the sample is close to the sensor, the magnification is equal to one. The resolution and field of view (FoV) of the image is then solely limited by pixel pitch, FoV, and size of the sensor. Therefore, a wider field of view is accessible for monitoring of the sample, compared to that of conventional optical microscopy [26, 27]. Moreover, since the pixel pitch of modern sensors can be as low as a few microns, lensless imaging can resolve structures of a few dozens of microns, such as bacterial microcolonies [28]. Finally, compared to conventional optical microscopes, lensless prototypes are compact and cost-efficient devices that can easily fit into standard laboratory incubators.

The lensless imaging approach has already been implemented to detect and count viral plaques and study eukaryote cells death dynamics in plaques using eukaryotic viruses with a 24-mm² FoV [29]. However, the authors do not provide results on the plaque growth kinetics. We therefore suggest using a similar approach but with a wider 3.3-cm² FoV CMOS sensor to study phage plaques. Here we report computer-assisted detection and counting of phage plaques as well as simultaneous measurement of the growth kinetics of nineteen plaques from the same sample. Finally, we propose a lensless imaging technique to monitor phage-resistance selection through imaging of the interior of phage plaques.

# Materials and methods

## Bacteria and bacteriophages

For the study of plaque count and growth rate, the methicillin-susceptible *S. aureus* (MSSA) strain Laus102 -isolated from a healthy patient at the University Hospital of Lausanne-and the lytic phage vB_SauM_2002 -isolated from a wastewater plant in Vidy, Lausanne, Switzerland-were used. Phage vB_SauM_2002 belongs to the Myoviridae family and harbors 99.92% nucleotide sequence identity over 91% of its ca. 140Kb genome with phage Sb-1 (GenBank: HQ163896.1).

For the study of resistance detection, the clinical isolate *K. pneumoniae* strain R405 TMP-8 and the lytic phage vB_Kpn_5055 -isolated from the wastewater plant in Vongy, Thonon-les-Bains, France-were used. The closest relative of vB_Kpn_5055 is the unclassified Przondovirus phage Kund-ULIP47 (GenBank: MK380015.1) with 98.40% nucleotide sequence identity over 97% of its ca. 41Kb genome.

Genomic phage DNA was purified following the traditional phenol/chloroform extraction protocol of the Center for Phage Technology of the Texas A&M University [30]. Purified DNA was sent to Eurofins Genomics Germany GmbH (Ebersberg, Germany) for sequencing. Phage genomic DNA fragment libraries were prepared with an optimized protocol and standard Illumina adapter sequences. Sequencing was performed with Illumina technology, NovaSeq 6000 (read mode 2 x 150bp).

For both bacteria and prior to the experiments, a loopful of 20% glycerol stock kept at -80°C was plated on a fresh Luria Bertani Broth (LB) agar plate and incubated aerobically at 37°C. Sixteen hours later, a single colony was resuspended in 5 mL fresh LB and incubated aerobically for 16 hr at 37°C with agitation (200rpm).

Phage solutions were prepared through simple phage amplification in liquid. Briefly, Bacterial cultures (Laus102 or R405 TMP-8) were grown at an OD₆₀₀nm of 0.3 before the corresponding phage were added at a multiplicity of infection (MOI) of 1 (vB_SauM_2002) or 0.01
(vB_Kpn_5055). The mixture was incubated aerobically at 37°C under 200 rpm agitation for 16 hr. After a centrifugation at 8000g for 15 min, the supernatant was filtered at 0.45 µm and titered using a soft-agar assay. Briefly, 250 µL of an overnight bacterial culture (Laus102 or R405 TMP-8) was mixed with serial dilutions of the corresponding phage solution and 5mL of molten soft-agar. The mixtures were poured in petri dishes and left at room temperature to allow the agar to solidify. Petri dishes were placed at 37°C for 16 hr and phage titers were determined by plaque counting. Phage solutions were stored at 4°C until further use. For the phage plaque counting experiment, 100µL of phage vB_SauM_2002 suspension was adjusted to 1.2×10³ PFU/ml and 250 µL of host bacteria S.aureus Laus102 culture are added to 5 mL of molten soft-agar. The mixture is smoothly shaken for homogenization and 2 mL are poured into a 35 mm diameter cell culture dish (Falcon® ref.353001). Finally, 1 mL of silicone oil AR 20 (Aldrich Chemistry ref. 10836-100ML) is added on top of the solidified soft agar layer. The silicon oil allows gas exchange to the growth medium while preventing it from drying.

To study the phage/bacteria interaction we choose the soft-agar approach over a culture in liquid broth since the former is a more faithful model of spatially structured environments in which most environmental bacteria reside (e.g. biofilms, soil, or plant and animal tissues) [31].

## Image acquisition and processing

**Custom lensless prototype.** The lensless imaging system is shown in Fig 1. The apparatus is composed of a detection system and an illumination module. The detection system consists of a 22.3 × 14.9 mm² complementary metal oxide semiconductor (CMOS) advanced photo system type-C (APS-C) sensor, repurposed from a Canon 1200D DSLR. The sensor–along with its accompanying supply and readout circuitry–is incorporated into a custom housing as shown in Fig 1. The sensor consists of an array of 5344 × 3516 pixels of 4.3 µm pitch. The Petri dish, prepared as indicated above, is placed directly on the sensor and illuminated from above by a light source. The light source consists of a monochromatic green LED (518 nm) coupled into a 200 µm core diameter multimode optical fiber (Thorlabs M72L02). Images are acquired using the built-in Canon software and converted from Canon native raw format (.CR2) to.tif

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Customized lensless prototype. (a) Drawing and (b) picture of the lensless configuration used to measure phage plaque growth. The prototype is small enough to fit in a standard laboratory incubator (here a Heraterm 18L).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0248917.g001
<figure><img src="image_5.png" /><figcaption>Fig 2. Algorithm structure to count phage plaques.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0248917.g002

files using Ddraw software. The images are then processed using two different algorithms; the first processes the entire image area to detect plaques, while the second processes only a cropped sub-image of each plaque to compute the growth rate.

**Computer-assisted plaque detection and counting.** The plaque detection is performed using an in-house algorithm that uses a computer-assisted method to count the phage lysis plaques. Images are processed in order to isolate the shape of the plaque from the background, and then are labelled to yield the number of detected phage plaques (see Fig 2).

**Plaque growth kinetics.** A different algorithm is used to study the growth kinetics of phage plaques using the Fiji software package [32] (see S1 File). Each plaque is isolated and processed separately (see S1 Dataset for examples of phage plaque image stacks). First, a Gaussian blur filter with two-pixel radius is applied in order to remove high spatial frequency noise. Then the image is thresholded and binarized using the built-in Otsu algorithm (see Fig 3) [33]. The plaque area is finally retrieved using the 'Analyze particles' module (see S2 Dataset). The
<figure><img src="image_6.png" /><figcaption>Fig 3. Image processing for plaque growth kinetics study. (a) Unprocessed phage plaque image (b) After a Gaussian blur filtering (c) Binarized final image.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0248917.g003

plaque radius is computed using:

$$r = \sqrt{\frac{A}{\pi}}$$

The time-derivative of the phage plaque radius ($\frac{dr}{dt}$, expressed in µm/hr) was calculated in order to retrieve the radial velocity of phage plaques, referred to hereafter as the ‘plaque growth rate’. In the case of a plaque whose radius grows linearly with time, the plaque growth rate is determined with the slope of curve, i.e. a constant. Linear regression was performed using Origin from OriginLab.

# Results and discussion

## Plaque detection and counting

An image of a single-layer soft agar assay recorded after 21 hr 44 min of incubation at 37°C is shown in Fig 4 (A) along with the corresponding size distribution of plaques at this time point (b) and the number of plaques detected over time (c). A movie of the acquisition is also

<figure><img src="image_7.png" /><figcaption>Fig 4. Acquisition of single-layer agar assay through lensless imaging. (a) Image of a single-layer soft agar assay after 21 hr 44 minutes of incubation at 37°C, with the dark background corresponding to the area of bacterial growth and bright disks to lysis plaques. The orange arrows indicate plaques that were not excluded from the growth kinetics study. (b) Corresponding plaque radius histograms and (c) cumulative graph of plaques detected as a function of time.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0248917.g004
<figure><img src="image_8.png" /><figcaption>**Fig 5.** Evolution of the mean gray level (red curve) of bacterial lawn and its variation (black curve). The gray level variation is the derivative of the mean gray level. It highlights the evolution of the growth rate of the bacteria inside the lawn. It is maximal after 2 hr of incubation.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0248917.g005

available as a (S1 Movie). The dark background on image (a) corresponds to the bacterial lawn while bright disks correspond to plaques. It is important to note that the background monotonically darkens during the first 10 hr of incubation as bacteria proliferate within the lawn, and then stabilizes at a constant gray level of around 1000 A.U. for the following 14 hr (Fig 5, red curve). Moreover, as indicated by the gray level variation maximum of -22.5 hr⁻¹ (Fig 5, black curve), the bacteria within the lawn reached a maximum of exponential growth rate after 2 hr of incubation.

The first plaque was detected only 4 hr after the start of incubation (Fig 4C) indicating that the phage was indeed lysing the selected bacterial strain. After 8 hr 20 min of incubation, a total of 23 different plaques were counted within the FoV and this number did not further increase after the 10 hr mark, likely indicating that all plaques within the considered FoV did appear. Knowing the sensor and culture dish area (3.3 and 9.6 cm², respectively), the volume of soft-agar in the dish (2 mL) and the volume of phage suspension introduced in the culture medium (100 µL), we estimated a phage titer of 1.7×10³ PFU/mL. This titer is only 30% different from the titer estimated with standard double layer agar assay. In other words, in our settings the continuous lensless monitoring of a single-layer soft agar assay can yield phage susceptibility results within 4 hr and phage titer in solution after only 8 hr 20 min. In contrast, the technique routinely used in laboratory requires around at least 12 hr of incubation and visual inspection of the agar plate to obtain similar data.

Of note, neither plaques at the edge of the FoV nor confluent plaques were considered, leaving us with 19 in individual plaques. A great variability in the plaque radius of these 19 plaques was observed at time point 21 hr 44 min, ranging from 147 µm to 546 µm with a mean of 325 ± 140 µm.

## Study of plaque growth kinetics

When isolating each plaque within the FoV (from Fig 4), we could detect and study the growth rate of the first plaque already after 3 hr 30 min of incubation (Fig 6). Of note, the plaques first appearing at 3 hr 30 min hr grow to a radius of >400 µm by the end of the experiment (at 21 hr 44 min) whereas the last-appearing plaques appearing at 10 hr would not exceed a radius of
<figure><img src="image_9.png" /><figcaption>Fig 6. Growth of plaque radius over time. (a) Plaque expansion on plates expressed as µm radius as a function of time for the nineteen plaques considered. Larger plaques were those that started to grow earlier than smaller plaques. (b) Representative examples of linear fits of the curves obtained in the first 10 hr of incubation for plaques of various final radii. The slopes of the fitted curves correspond to the growth rate of the respective plaques.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0248917.g006

150 µm by the end of the image acquisition process. This likely reflects delayed adsorption of some phage vB_SauM_2002 on cells of its bacterial host S. aureus Laus102, with larger plaques resulting from phage particles that adsorbed early and smaller plaques resulting from phage particles that adsorbed later during the experiment [34].

We separated the plaque expansion curves into three distinct phases (Fig 6A).

**Phase I.** During the first 10 hr of incubation (phase I), the plaque growth, reflected by the increase in radius over time), was linear (Fig 6A). Therefore, we derived a constant growth rate for each plaque during phase I using the slopes of the linear fits (Fig 6B). As observed with a representative selection of five different plaques of final radii ranging between 175 to 520 µm (Fig 6B), the growth rates of the larger plaques were higher than the ones of the smaller plaques, highlighting that the plaque growth rates correlated with the final radius size of the plaques.

Linear growth of plaques has been reported previously [22–24], and was expected from theoretical models suggesting linear plaque growth rate $\frac{dr}{dt}$ dependent only on the phage diffusion rate D and the latency time L [21] following the equation:

$$ \frac{dr}{dt} = 10 \cdot \sqrt{\frac{D}{L}} $$

However, this overly simplistic model did not take into account the differing plaque growth rates in relation to plaque size and the surrounding bacterial density. More complex models demonstrated that a higher initial bacterial density inside the top-agar layer resulted in smaller
<figure><img src="image_10.png" /><figcaption>**Fig 7.** Study of plaque growth rates over time. Plaque growth rates vs time. (a) At 15 hr of incubation, the mean growth rate stabilized for all plaques independent of the size reached by the end of the experiment. (b) The mean growth rates during phase III were determined for each group.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0248917.g007

plaques as well as lower plaque growth rates [35]. Moreover, as explained above, the phage adsorption rate contributes to differences in the final size of a plaque. Taken together, these parameters induce a variance in the timing of the first infection event for each individual phage. As observed in phase I, it might be assumed that the first infection event that produces a phage progeny happened with a significant delay in smaller plaques compared to larger plaques. Indeed, the smallest plaques were first observed on plates approximately 5 hr after the first large plaque was observed. According to this significant delay, a higher initial bacterial density surrounding smaller plaques might be assumed. Therefore, our experiments would support the previously observed correlation between a higher initial bacterial density around small plaques associated to a lower plaque growth rate [35]. The associated lower growth rate of late-appearing plaques would be further supported by the differences in plaque growth rates observed in phase I (Fig 6B).

**Phase II.** During phase II, starting from 10 hr after the start of incubation, the plaque growth was no longer linear. Therefore, we analysed the growth rate of each plaque during phase II by computing the time derivative of the plaque radii (Fig 6A). For clarity, the plaques were classified into five different groups according to the size of their final radii. While plaque growth rate decreases for all groups, the plaque growth rate still correlates with the plaque size, being higher in larger plaques (Fig 7) as observed during phase I (see above). Therefore, we observed that despite a decrease in growth rate, larger plaques always grew faster than smaller plaques in both phase I and phase II.

In our experiments, the observed decreases in plaque growth rate during phase II could be correlated with the maturation of the bacterial lawn. Indeed, starting at 12 hr (i.e. at the beginning of phase II), the variation of the gray level of the bacterial lawn area was close to $0 \text{ hr}^{-1}$ (Fig 5). This indicates that the amount of light transmitted through the bacterial lawn stopped decreasing as the bacteria ceased to divide. Therefore, our results fully support the previously reported observation that the plaque growth rate decreases as the host density increases [36].
In addition, hindered diffusion of phage particles in a dense bacterial lawn could contribute to the decrease in plaque growth rate during phase II. Indeed, the growth medium is not an homogeneous medium but rather a hydrogel (soft agar) which includes many bacterial micro-colonies which may act as diffusional barriers [37].

Therefore, our results are in agreement with a general decrease in plaque growth rates in phase II resulting from impeded diffusion of the phages due to the high density of the bacterial lawn.

**Phase III.** Finally, after 15 hr of incubation (phase III), the plaque growth rates stabilized (Fig 6B). However-as in phase I and II-the plaque growth rates somehow correlated with the final plaque size as the largest size group (500 to 600 µm) had a mean growth rate of 7.15 ± 0.2 µm/hr while the smallest size group (100 to 200 µm) had a mean plaque growth rate of 2.1 ± 0.03 µm/hr. However, and interestingly, the three intermediary size group displayed similar plaque growth rates at around 4 µm/hr.

Therefore, our results suggest that the bacterial lawn reached full confluence at around 15 hr and that, despite consistently more rapid growth rates for larger plaques, the growth rates stabilized in such a homogeneous bacterial lawn.

## Imaging of phage-resistant microcolonies within plaques

Phage-resistant bacterial clones selected under phage pressure can develop as microcolonies within plaques [38, 39] and be imaged through lensless imaging (Fig 8). Indeed, the resolution of the sensor being limited to the pixel pitch (here 4.3 µm), we were able to image elements with sizes in the same order of magnitude. For instance, we imaged resistant microcolonies of

<figure><img src="image_11.png" /><figcaption>**Fig 8.** Lensless imaging of *K. pneumoniae* bacterial microcolonies developing within plaques after 23 hr of incubation. The red arrows indicate microcolonies within plaques.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0248917.g008
K. pneumoniae with a diameter of 60 µm (Fig 8, red arrows). A movie of the acquisition is also available as a (S2 Movie).

These microcolonies probably represent *K. pneumoniae* bacteria that are resistant to phage vB_Kpn_5055, but this would have to be confirmed through isolation and reculturing of these colonies followed by a susceptibility test. However this result highlights that lensless imaging can be used as a preliminary tool to detect possible phage resistance to a bacterial strain. Indeed the wide 3.3 cm² FoV coupled to the 4.3 µm pixel pitch enable the study of the evolution of millimeter-scale plaques while still permitting a resolution which is high enough to simultaneously visualize elements on the scale of around a few dozen micrometers.

## Conclusion

In this work, we report the use of a custom lensless device, over a 3.3 cm² area sensor, for the continuous monitoring of phage lysis plaques. Leveraging wide-field lensless imaging, we performed a computer-assisted detection and counting of plaques that allowed detection of bacterial susceptibility to phages within 3 hr and accurate estimation of infectious titer in only 8 hr 20 min.

Moreover, by investigating the growth rates of 19 isolated plaques we confirmed the previously observed correlation between bacterial density and phage diffusion in the soft-agar layer. However, lensless imaging allowed us to identify three very distinct phases during the growth of a plaque independent of the size of the final plaque. Moreover, we demonstrated that plaques harboring the smaller sizes at the end of the experiment were always associated with a delayed onset of growth in phase I as well as lower growth rates during the three phases.

In addition, we showed that lensless imaging could be a powerful tool to screen the emergence of phage-resistant bacteria by imaging bacterial microcolonies within plaques. Additional experiments are needed to evaluate whether this technique could lead to an estimation of the phage-resistance frequency.

Future work will also focus on the development of an algorithm allowing morphological classification of plaques (and therefore phages) according to their plaque growth rates but also their morphotypes (though we did not discuss this aspect here). Moreover, we think that lensless imaging could reveal a very important approach to study phage-antibiotic synergy (PAS). Indeed, it has for instance been reported that the inclusion of an antibiotic in the agar layer leads to significantly larger plaques [6], making visible plaques that would otherwise be invisible to the naked eye [40]. Finally, a setup based on multiple image sensors working simultaneously is currently under development, which will increase the field of view and consequently the amount of data acquired at the same time.

## Supporting information

**S1 Movie.** Lensless acquisition of single-layer soft agar assay of *S. aureus* strain Laus102 with lytic phage vB_SauM_2002.
(AVI)

**S2 Movie.** Lensless acquisition of single-layer soft agar assay of *K. pneumoniae* strain R405 TMP-8 and the lytic phage vB_Kpn_5055.
(AVI)

**S1 Dataset.** Phage plaque image stack from vB_SauM_2002 on *S. aureus* Laus102 acquisition.
(ZIP)
**S2 Dataset. Measured phage plaque areas of phage vB_SauM_2002 on S. aureus Laus102 strain prior to processing.**
(CSV)

**S1 File. FIJI macro used to measure phage plaque kinetics.**
(ZIP)

**Acknowledgments**

We thank Patrice Nordmann for the kind gift of the *K. pneumoniae* strain R405 TMP-8.

**Author Contributions**

**Conceptualization:** Prisca Perlemoine, Pierre R. Marcoux, Emmanuel Picard, Marc Zelsmann.

**Investigation:** Prisca Perlemoine, Aurélie Marchet.

**Methodology:** Prisca Perlemoine, Emmanuel Picard, Marc Zelsmann.

**Resources:** Aurélie Marchet, Grégory Resch.

**Software:** Grégoire Mugnier.

**Supervision:** Emmanuel Hadji, Eric Lacot.

**Writing – original draft:** Prisca Perlemoine.

**Writing – review & editing:** Pierre R. Marcoux, Grégory Resch, Larry O'Connell.

**References**

1. Abedon ST, Kuhl SJ, Blasdel BG, Kutter EM. Phage treatment of human infections. Bacteriophage. 2011 Mar; 1(2):66–85. https://doi.org/10.4161/bact.1.2.15845 PMID: 22334863

2. Ravat F, Jault P, Gabard J. Bactériophages et phagothérapie: utilisation de virus naturels pour traiter les infections bactériennes. Ann Burns Fire Disasters. 2015 Mar 31; 28(1):13–20. PMID: 26668557

3. Jault P, Leclerc T, Jennes S, Pirnay JP, Que Y-A, Resch G, et al. Efficacy and tolerability of a cocktail of bacteriophages to treat burn wounds infected by Pseudomonas aeruginosa (PhagoBurn): a randomised, controlled, double-blind phase 1/2 trial. The Lancet Infectious Diseases. 2019 Jan; 19(1):35–45. https://doi.org/10.1016/S1473-3099(18)30482-1 PMID: 30292481

4. Ujmajuridze A, Chanishvili N, Goderdzishvili M, Leitner L, Mehnert U, Chkhotua A, et al. Adapted Bacteriophages for Treating Urinary Tract Infections. Front Microbiol. 2018; 9:1832. https://doi.org/10.3389/fmicb.2018.01832 PMID: 30131795

5. Wright A, Hawkins CH, Änggård EE, Harper DR. A controlled clinical trial of a therapeutic bacteriophage preparation in chronic otitis due to antibiotic-resistant *Pseudomonas aeruginosa*; a preliminary report of efficacy. Clinical Otolaryngology. 2009 Aug; 34(4):349–57. https://doi.org/10.1111/j.1749-4486.2009.01973.x PMID: 19673983

6. Comeau AM, Tétart F, Trojet SN, Prère M-F, Krisch HM. Phage-Antibiotic Synergy (PAS): β-Lactam and Quinolone Antibiotics Stimulate Virulent Phage Growth. Fox D, editor. PLoS ONE. 2007 Aug 29; 2 (8):e799. https://doi.org/10.1371/journal.pone.0000799 PMID: 17726529

7. Oechslin F, Piccardi P, Mancini S, Gabard J, Moreillon P, Entenza JM, et al. Synergistic interaction between phage therapy and antibiotics clears *Pseudomonas aeruginosa* infection in endocarditis and reduces virulence. INFDIS. 2016 Dec 22;jiw632.

8. Pirnay J-P, Verbeken G, Ceyssens P-J, Huys I, De Vos D, Ameloot C, et al. The Magistral Phage. Viruses. 2018 Feb 6; 10(2):64. https://doi.org/10.3390/v10020064 PMID: 29415431

9. Oechslin F. Resistance Development to Bacteriophages Occurring during Bacteriophage Therapy. Viruses. 2018 Jun 30; 10(7):351. https://doi.org/10.3390/v10070351 PMID: 29966329

10. Mattila S, Ruotsalainen P, Jalasvuori M. On-Demand Isolation of Bacteriophages Against Drug-Resistant Bacteria for Personalized Phage Therapy. Front Microbiol [Internet]. 2015 Nov 13 [cited 2020 Apr
13];6. Available from: http://journal.frontiersin.org/Article/10.3389/fmicb.2015.01271/abstract PMID: 26617601

11. Seeley ND, Primrose SB. A Review: The isolation of bacteriophages from the environment. Journal of Applied Bacteriology. 1982 Aug; 53(1):1-17. https://doi.org/10.1111/j.1365-2672.1982.tb04729.x PMID: 6757237

12. Castro-Mejía JL, Muhammed MK, Kot W, Neve H, Franz CMAP, Hansen LH, et al. Optimizing protocols for extraction of bacteriophages prior to metagenomic analyses of phage communities in the human gut. Microbiome. 2015 Dec; 3(1):64. https://doi.org/10.1186/s40168-015-0131-4 PMID: 26577924

13. Carvalho C, Susano M, Fernandes E, Santos S, Gannon B, Nicolau A, et al. Method for bacteriophage isolation against target *Campylobacter* strains. Letters in Applied Microbiology. 2010 Feb; 50(2):192-7. https://doi.org/10.1111/j.1472-765X.2009.02774.x PMID: 20002571

14. Dias RS, Eller MR, Duarte VS, Pereira ÂL, Silva CC, Mantovani HC, et al. Use of phages against antibiotic-resistant Staphylococcus aureus isolated from bovine mastitis1. Journal of Animal Science. 2013 Aug 1; 91(8):3930-9. https://doi.org/10.2527/jas.2012-5884 PMID: 23893999

15. Goyal SM, Zerda KS, Gerba CP. Concentration of coliphages from large volumes of water and wastewater. Appl Environ Microbiol. 1980 Jan; 39(1):85-91. https://doi.org/10.1128/AEM.39.1.85-91.1980 PMID: 7356323

16. Leshkasheli L, Kutateladze M, Balarjishvili N, Bolkvadze D, Save J, Oechslin F, et al. Efficacy of newly isolated and highly potent bacteriophages in a mouse model of extensively drug-resistant Acinetobacter baumannii bacteraemia. Journal of Global Antimicrobial Resistance. 2019 Dec 1; 19:255-61. https://doi.org/10.1016/j.jgar.2019.05.005 PMID: 31100499

17. O'Flaherty S, Ross RP, Flynn J, Meaney WJ, Fitzgerald GF, Coffey A. Isolation and characterization of two anti-staphylococcal bacteriophages specific for pathogenic Staphylococcus aureus associated with bovine infections. Lett Appl Microbiol. 2005 Dec; 41(6):482-6. https://doi.org/10.1111/j.1472-765X.2005.01781.x PMID: 16305674

18. Primrose SB, Day M. Rapid Concentration of Bacteriophages from Aquatic Habitats. Journal of Applied Bacteriology. 1977 Jun; 42(3):417-21. https://doi.org/10.1111/j.1365-2672.1977.tb00709.x PMID: 407206

19. Ferry T, Kolenda C, Batailler C, Gustave C-A, Lustig S, Malatray M, et al. Phage Therapy as Adjuvant to Conservative Surgery and Antibiotics to Salvage Patients With Relapsing S. aureus Prosthetic Knee Infection. Front Med [Internet]. 2020 [cited 2020 Dec 15];7. Available from: https://www.frontiersin.org/articles/10.3389/fmed.2020.570572/full PMID: 33304911

20. Philipson CW, Voegtly LJ, Lueder MR, Long KA, Rice GK, Frey KG, et al. Characterizing Phage Genomes for Therapeutic Applications. Viruses [Internet]. 2018 Apr 10 [cited 2020 Dec 15]; 10(4). Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5923482/ https://doi.org/10.3390/v10040188 PMID: 29642590

21. Koch AL. The growth of viral plaques during the enlargement phase. Journal of Theoretical Biology. 1964 May; 6(3):413-31. https://doi.org/10.1016/0022-5193(64)90056-6 PMID: 5875213

22. Kaplan DA, Naumovski L, Rothschild B, John Collier R. Appendix: A model of plaque formation. Gene. 1981 Apr; 13(3):221-5. https://doi.org/10.1016/0378-1119(81)90027-5 PMID: 6790348

23. Lee Y, Yin J. Imaging the propagation of viruses. Biotechnology and Bioengineering. 1996; 52(3):438-42. https://doi.org/10.1002/(SICI)1097-0290(19961105)52:3<438::AID-BIT11>3.0.CO;2-F PMID: 18629915

24. Lee Y, Yin J. Detection of evolving viruses. Nat Biotechnol. 1996 Apr; 14(4):491-3. https://doi.org/10.1038/nbt0496-491 PMID: 9630926

25. Gabor D. A New Microscopic Principle. Nature. 1948 May; 161(4098):777-8. https://doi.org/10.1038/161777a0 PMID: 18860291

26. Zhang Y, Alexander M, Yang S, Bian Y, Botvinick E, Lakey JRT, et al. High-Throughput Screening of Encapsulated Islets Using Wide-Field Lens-Free On-Chip Imaging. ACS Photonics. 2018 Jun 20; 5(6):2081-6.

27. Ozcan A, Demirci U. Ultra wide-field lens-free monitoring of cells on-chip. Lab Chip. 2008; 8(1):98-106. https://doi.org/10.1039/b713695a PMID: 18094767

28. Maeda Y, Dobashi H, Sugiyama Y, Saeki T, Lim T, Harada M, et al. Colony fingerprint for discrimination of microbial species based on lensless imaging of microcolonies. Lee B-L, editor. PLoS ONE. 2017 Apr 3; 12(4):e0174723. https://doi.org/10.1371/journal.pone.0174723 PMID: 28369067

29. Han C, Yang C. Viral plaque analysis on a wide field-of-view, time-lapse, on-chip imaging platform. Analyst. 2014; 139(15):3727-34. https://doi.org/10.1039/c3an02323k PMID: 24611157
30. Protocol for Phage DNA Extraction with Phenol:Chloroform [Internet]. [cited 2020 Dec 16]. Available from: https://cpt.tamu.edu/wordpress/wp-content/uploads/2018/09/Phage-DNA-Extraction-by-PhenolChloroform-Protocol.pdf

31. Abedon ST, Yin J. Bacteriophage Plaques: Theory and Analysis. In: Clokie MRJ, Kropinski AM, editors. Bacteriophages [Internet]. Totowa, NJ: Humana Press; 2009 [cited 2020 Mar 27]. p. 161–74. (Walker JM, editor. Methods in Molecular Biology; vol. 501). Available from: http://link.springer.com/10.1007/978-1-60327-164-6_17

32. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nature Methods. 2012 Jul; 9(7):676–82. https://doi.org/10.1038/nmeth.2019 PMID: 22743772

33. Otsu N. A Threshold Selection Method from Gray-Level Histograms. IEEE Trans Syst, Man, Cybern. 1979 Jan; 9(1):62–6.

34. Adams MH. Bacteriophages [Internet]. Interscience Publishers; 1959. 592 p. Available from: https://archive.org/details/bacteriophages00adam/mode/2up

35. Mitarai N, Brown S, Sneppen K. Population Dynamics of Phage and Bacteria in Spatially Structured Habitats Using Phage λ and Escherichia coli. Silhavy TJ, editor. J Bacteriol. 2016 Jun 15; 198 (12):1783–93. https://doi.org/10.1128/JB.00965-15 PMID: 27068593

36. Ortega-Cejas V, Fort J, Méndez V, Campos D. Approximate solution to the speed of spreading viruses. Phys Rev E. 2004 Mar 24; 69(3):031909. https://doi.org/10.1103/PhysRevE.69.031909 PMID: 15089324

37. Yin J, McCaskill JS. Replication of viruses in a growing plaque: a reaction-diffusion model. Biophysical Journal. 1992 Jun; 61(6):1540–9. https://doi.org/10.1016/S0006-3495(92)81958-6 PMID: 1617137

38. Pouillot F, Chomton M, Blois H, Courroux C, Noelig J, Bidet P, et al. Efficacy of Bacteriophage Therapy in Experimental Sepsis and Meningitis Caused by a Clone O25b:H4-ST131 Escherichia coli Strain Producing CTX-M-15. Antimicrob Agents Chemother. 2012 Jul; 56(7):3568–75. https://doi.org/10.1128/AAC.06330-11 PMID: 22491690

39. Postic B, Finland M. Observations on bacteriophage typing of Pseudomonas aeruginosa. J Clin Invest. 1961 Nov 1; 40(11):2064–75. https://doi.org/10.1172/JCI104432 PMID: 14488221

40. Łoś JM, Golec P, Węgrzyn G, Węgrzyn A, Łoś M. Simple Method for Plating Escherichia coli Bacteriophages Forming Very Small Plaques or No Plaques under Standard Conditions. AEM. 2008 Aug 15; 74 (16):5113–20. https://doi.org/10.1128/AEM.00306-08 PMID: 18586961