Phenotypic antimicrobial susceptibility testing with deep learning

School of Electrical, Computer and Energy Engineering, Arizona State University, Tempe,. Arizona 85287 ... Consequently, healthcare providers often fa...
2 downloads 8 Views 1MB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Phenotypic antimicrobial susceptibility testing with deep learning video microscopy Hui Yu, Wenwen Jing, Rafael Iriya, Yunze Yang, Karan Syal, Manni Mo, Thomas E Grys, Shelley E Haydel, Shaopeng Wang, and Nongjian Tao Anal. Chem., Just Accepted Manuscript • DOI: 10.1021/acs.analchem.8b01128 • Publication Date (Web): 20 Apr 2018 Downloaded from http://pubs.acs.org on April 20, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Phenotypic antimicrobial susceptibility testing with deep learning video microscopy Hui Yu1,2†, Wenwen Jing2†, Rafael Iriya2,3, Yunze Yang2, Karan Syal2, Manni Mo2,3, Thomas E. Grys5, Shelley E. Haydel6,7, Shaopeng Wang2,3, and Nongjian Tao2-4* 1

Institute for Personalized Medicine, School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, 200030, China.

2

Biodesign Center for Biosensors and Bioelectronics, Arizona State University, Tempe, Arizona 85287, USA. 3

State Key Laboratory of Analytical Chemistry for Life Science, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China. 4

School of Electrical, Computer and Energy Engineering, Arizona State University, Tempe, Arizona 85287, USA. 5

Department of Laboratory Medicine and Pathology, Mayo Clinic, Phoenix, Arizona 85054, USA. 6

Biodesign Center for Immunotherapy, Vaccines, and Virotherapy, Arizona State University, Tempe, Arizona 85287, USA. 7

School of Life Sciences, Arizona State University, Tempe, Arizona 85287, USA.

*To whom correspondence should be addressed: [email protected] † These authors contribute equally to the project.

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract: Timely determination of antimicrobial susceptibility for a bacterial infection enables precision prescription, shortens treatment time and helps minimize the spread of antibiotic resistant infections. Current antimicrobial susceptibility testing (AST) methods often take several days and thus impede these clinical and health benefits. Here, we present an AST method by imaging freely moving bacterial cells in urine in real time and analyzing the videos with a deep learning algorithm. The deep learning algorithm determines if an antibiotic inhibits a bacterial cell by learning multiple phenotypic features of the cell without the need for defining and quantifying each feature. We apply the method to urinary tract infection, a common infection that affects millions of people, to determine the minimum inhibitory concentration (MIC) of pathogens from both bacteria spiked urine and clinical infected urine samples for different antibiotics within 30 minutes, and validate the results with the gold standard broth macrodilution method. The deep learning video microscopy-based AST holds great potential to contribute to the solution of increasing drug-resistant infections.

Antimicrobial resistance or emergence of “superbugs” has become a global health epidemic.1-3 Acceleration of this epidemic in recent years is primarily caused by the widespread overuse and misuse of antibiotics, prompting bacteria to evolve and develop resistance.4 To address this threat, it is critical to accurately prescribe effective antibiotics for the patient, which necessitates timely antimicrobial susceptibility testing (AST). Current AST technologies, including disk diffusion and broth dilution methods, often take several days to complete.5-7 Consequently, healthcare providers often face a dilemma: delaying treatment or prescribing potentially ineffective or broad-range empiric therapy while awaiting AST results. A rapid AST technology would help identify antimicrobial susceptibility at the earliest stage of infection, and allow

ACS Paragon Plus Environment

Page 2 of 23

Page 3 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

healthcare providers to prescribe narrow-spectrum antibiotic treatment, thus reducing patient mortality and spread of antimicrobial resistance.8 Innovative AST technologies have been pursued using either genotypic or phenotypic approaches.9 The former detects genes responsible for conferring drug resistance,10-12 which is powerful, but requires prior knowledge of the genes, detects only the potential of antibiotic resistance, and cannot differentiate viable and non-viable bacterial cells. The latter detects if a bacterium can be effectively inhibited or killed by an antibiotic by measuring its phenotypic features

using various

detection

techniques.13-24

These techniques

typically require

immobilization of bacteria on a sensor surface, in a gel, or in sophisticated microfluidic channels for imaging and detection, which raises practical difficulties in testing clinical samples. Furthermore, each of them typically measures one phenotypic feature only, limiting their capability in infections by different pathogens. Optical microscopy13,14 is especially attractive due to its capability in imaging multiple phenotypic features of individual single cell, including cell size, morphology, motion, and division. However, defining and quantifying these features with the traditional image processing method is challenging because a cell can grow in size, change in shape, divide over time, rotate, move around in the solution, and move in and out of the microscopic field of view. These challenges are further highlighted when considering that most optical images are 2D representations of 3D bacterial cells that rotate and move in solution.25 Here, we describe an AST technology that images single, non-immobilized bacterial cells and analyzes multiple phenotypic features and responses of the cells automatically with a deep learning (DL) algorithm (Figure 1). DL is an exciting new area of artificial intelligence using large neural networks, and has been used for cell segmentation and classification based on static

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

images,26,27 but applying it to AST with live videos of bacteria as input data is non-trivial, and has not yet been demonstrated. Our approach maximizes the speed and accuracy of AST by learning multiple phenotypic features at the pixel level without having to define and then quantify each of them. Its self-learning capability allows improvement of AST accuracy over time as the number of analyzed samples increases. To demonstrate this DL video microscopyenabled AST (DLVM-AST), we focused on Escherichia coli (E. coli), a bacterial pathogen that is the most common cause of urinary tract infections (UTI), and five relevant antibiotics for treating UTI: polymyxin B (PMB), streptomycin, ciprofloxacin, aztreonam, and ampicillin. These antibiotics kill or inhibit E. coli via different mechanisms, resulting in different cell phenotypic changes, such as motion, morphology and division changes. We evaluated the capability of DLVM-AST for automatically identifying and analyzing antibiotic-mediated inhibition of bacterial cells using E. coli as an example and determining the minimum inhibitory concentrations (MIC). DLVM-AST results were also compared to results obtained by the traditional imaging processing algorithm13,14 and the gold standard broth macrodilution (BMD) method. Finally, clinical isolates were tested to evaluate the translation of this technology for clinical use. Experimental Section Materials Unfiltered human urine samples (Lot#: BRH1041997) and E. coli (ATCC 43888; Biosafety Level 1 organism that does not produce either Shiga-like I or II toxins and lacks the genes for these toxins) were purchased from BioreclamationIVT Co. and Fisher Scientific, respectively. Antibiotics, including polymyxin B (PMB), ampicillin, streptomycin, ciprofloxacin, and

ACS Paragon Plus Environment

Page 4 of 23

Page 5 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

aztreonam, and all other reagents were purchased from Sigma-Aldrich. The antibiotic powders were stored in dark at -2-8°C.

Figure 1 Schematics of deep learning video microscopy-based antimicrobial susceptibility testing (DLVM-AST) method. Urine samples were mixed with antibiotics at different concentrations and imaged with a microscope in a microfluidic channel without immobilization. Videos of bacteria were recorded over time, and compressed into static images containing single cell features. A deep learning algorithm was used to determine the minimum inhibitory concentration (MIC) value from the sub-videos.

Antibiotic preparations Stock solutions of PMB, ampicillin, and streptomycin at concentrations of 200 µg/mL were prepared by directly dissolving the antibiotics in ultrapure water. Ciprofloxacin and aztreonam were first dissolved in 0.1 M HCl (1:60, m/V) and dimethylformamide: methanol solution (1: 1, V/V), respectively, and then diluted in ultrapure water to obtain stock concentrations of 200 µg/mL. These antibiotic stock solutions were stored in dark at -80°C. Before AST, the antibiotic stock solutions were thawed to room temperature and diluted in ultrapure water to various

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

concentrations for AST, following guidelines recommended by the Clinical and Laboratory Standards Institute (CLSI)28. Growth and preparation of E. coli Frozen E. coli strains were thawed, and 50 µL of which were cultured in 5 mL of Luria-Bertani (LB) medium (Per liter: 10 g peptone 140, 5 g yeast extract, and 5 g sodium chloride) at 37°C and 150 rpm for 16 hours. Saturated cultures in the volume of 20 µL were diluted into 5 mL of fresh LB medium, and growth continued at 37°C with 150 rpm for 1 hour to attain a logarithmic phase of growth. Bacterial cells were collected by centrifugation at 450g for 5 min and suspended in urine to a concentration of 2 × 107 cells/mL29. This concentration was determined by measuring the extinction coefficients for E.coli from the Optical Density (OD600) reading taken with a spectrophototometer (NanoDrop™ 2000/2000c Spectrophotometers, Thermo Scientific). The calibration factor for bacterial cell cultures estimation was 8 × 108 cells/mL per OD600 unit. Before use, the bacteria-spiked urine samples were filtered using a 5 µm syringe filter (EMD Millipore) to remove large particles. Fabrication and structure of the microfluidic chip We used a microfluidic chip with a channel volume less than 100 nL to generate a stable microenvironment for the bacterial cells (Figure S4, Supporting Information). The microfluidic chips were fabricated by multilayer soft lithography,30-32 including pneumatic control and fluidic layers made of PDMS (RTV 615, the ratio of A/B is 5: 1) and PDMS (RTV 615, the ratio of A/B is 10: 1, Momentive Specialty Chemicals), respectively. The fluidic layer included a detection channel (0.5 cm long, 200 µm wide and 25 µm high) and inlet and outlet channels. The control and fluidic layers were aligned by thermopolymerization reaction and bonded on a glass slide with oxygen plasma. The mold of the control layer was made of negative photoresist (SU8-2025,

ACS Paragon Plus Environment

Page 6 of 23

Page 7 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Microchem), and the mold of the fluidic layer was made from a positive photoresist (AZ-50XT, AZ Electronic Materials USA Corp.). The microfluidic chip has six parallel detection channels, which allowed AST detection with different concentration of antibiotics simultaneously. The fluids were kept inside the detection channel by closing the valves during video recording. AST with video microscopy The microfluidic chip was placed on an inverted microscope (Olympus IX-81) with a 40× phase contrast objective lens, and imaged with a CCD camera (Pike-032B, Allied Vision Technologies, Newburyport, MA). A 200-µL bacterial suspension (2×107 cells/mL) was mixed with an equal volume of antibiotic solution for each antibiotic concentration or an equal volume of water as a control experiment. These mixed solutions were injected into different microfluidic channels simultaneously. After the microfluidic channels were fulfilled with these mixed solutions, two microfluidic valves of each detection channel were closed simultaneously to generate a stable microenvironment. Videos of the bacterial cells were recorded at 100 frame per second (fps) immediately (0 min) and after every 30 min. Each video lasted for 30 seconds. The raw images were batch-converted to 16-bit tiff format using a Matlab program, and pre-processed to minimize background artifacts (Figure S1, Supporting Information) before being processed with the DL model. To include enough cells for reliable results, we integrated bacterial cells from three videos recorded from different experiments. Each experiment was repeated independently for three times. Deep learning model The DL model was implemented with the TensorFlowTM (Google Inc.), an open-source software library for Machine Intelligence. The training dataset included 1000 static single cell images, with 500 “inhibited”, and 500 “uninhibited” (control) cells for each antibiotic tested. Each model

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

was trained iteratively for 1000 iterations. In each iteration, 50 static images were randomly selected from the training dataset and input into the model. We examined the accuracy of the trained model using 300 static images (150 labeled “inhibited” and 150 labeled “uninhibited”) upon exposure to each antibiotic. All the computations were performed with a desktop computer (Intel ® Core ™ i7-4790 CPU @3.60GHz). AST with broth macrodilution (BMD) method For comparison, AST was also performed by the BMD method (CLSI gold standard28). The adjusted inoculum E.coli suspension is diluted in Cation-adjusted Mueller Hinton Broth (CAMHB, Sigma-Aldrich) and the concentration of E.coli cell cultures are adjusted to 1 × 106 cells/mL based on UV-Vis spectrophotometer (NanoDrop 2000, Thermo Fisher) readings at OD600. Within 15 minutes after the inoculum has been prepared, 1 mL of the adjusted inoculum is added to tubes containing 1 mL of antibiotics in two-fold dilution series or only broth (control group) and mix. This results in a 1:2 dilution of each antibiotics and inoculum concentration. After inoculation, each tube contains approximately 5 × 105 cells/mL. After incubation at 37 °C for 16 hours, the MIC values can be read as the lowest concentration without visible growth. This test was performed in triplicate. Clinical sample preparation and testing. This study was approved by the Arizona State University’s Institutional Review Board. All urine clinical samples were stored at 4℃ and detected within two hours after collection. The raw urine samples were filtered by a 5 microns filter membrane, which removes large substances such as epithelial cells and blood cells. The filtration removes less than 5% bacterial cells as confirmed by cell counting under microscope.

ACS Paragon Plus Environment

Page 8 of 23

Page 9 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Results and Discussion The workflow of DLVM-AST consists of 1) imaging single bacterial cells in urine samples with a phase contrast microscope before, during, and after exposure to each antibiotic at different concentrations, 2) compressing the videos into static images while preserving essential phenotypic features, 3) feeding the static images into a DL model (pre-trained with thousands of images), and 4) obtaining antimicrobial susceptibility and MIC for the bacterial strain (Figure 1). We describe below each of these steps. Imaging single bacterial cells without immobilization We imaged bacterial cells spiked in a urine sample directly in a microfluidic chip (see Material and Methods for details) without immobilizing them onto a sensor surface or in gels. This simplified sample preparation also allowed the bacterial cells to move freely in urine solution (e.g., swimming and tumbling)33,34, thus capturing phenotypic features that are not trackable with immobilized bacterial cells. We observed that the bacterial cells frequently moved in and out of the microscope view and focus (Movie S1, Supporting Information), and following each of them over time proved difficult with the conventional image processing method. DLVM-AST overcame this difficulty because it did not rely on tracking a specific feature of a bacterial cell. A time averaging approach was adopted to remove the static background in recorded video, which greatly improved the signal noise ratio for post analysis (Figure S1, Supporting Information). Compressing the videos into static images without losing key phenotypic features In principle, the raw bacterial videos could be used as inputs for a DL model, but the computational expense is impractical. In fact, few DL applications could use videos as input data directly even with high performance super-computers.35-37 To overcome this difficulty, we developed a method to compress the raw videos into static images without losing key phenotypic

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

features. This approach imitates human vision, which pre-processes raw images in the vision system to reduce complexity before passing them to the brain. The key phenotypic features in the video include cell division, motion, and morphology (Figure S2, Supporting Information). A bacterial cell grows and divides, so cell division is a useful phenotypic feature to indicate if the cell is killed or its growth is effectively inhibited by an antibiotic. Cell motion and morphology may also change when exposing the cell to antibiotics, thus serving additional phenotypic features for AST. We found that the motion and morphology changes of E. coli were different for different antibiotics. For example, while PMB decreased bacterial motion, aztreonam caused the bacterial cells to elongate. These observations underscore the value of tracking both the motion and morphology as phenotypic features in addition to cell division, particularly for slowly dividing bacterial strains.

We compressed the bacterial videos while preserving the essential phenotypic features described above using the strategy shown in Figure 2. It compresses each video (duration of 1 s) into two sets of static images, capturing the morphology and motion of a single cell, respectively. The image containing the morphological feature is a snapshot of each individual bacterial cell (inset e in Figure2). In contrast, the image containing the bacterial movement is the superposition of the binarized frames in the video, which represents the motion of the cell as a trace in the binary image (inset g in Figure2). The cell division feature is determined by the number of individual bacterial cells after clustering (inset c in Figure2). The two sets of static images are merged into a single set of images (inset h in Figure2) as the input data for the DL model. The model automatically learns and determines antimicrobial susceptibility from the input data at pixel level without specifically extracting high-level features. This strategy is different from the traditional

ACS Paragon Plus Environment

Page 10 of 23

Page 11 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

cell imaging analysis, which defines and then quantifies each of the phenotypic features at single cell level, such as size, perimeter length, and speed.

Figure 2 Compressing the bacterial video into computationally tractable dataset without losing key features. The raw video (a) is first converted into a binary video (b) by comparing the intensity value of each pixel with a threshold value (±1000 a.u.). To obtain a compressed image for tracking of the morphological feature, the binary images are then segmented in each of the binary images to create a set of single cell image sequence (c). Next, one binary image (d) that has the maximum number of pixels is chosen from the single cell image sequence, and the corresponding image in the raw video (e) is identified and selected to provide the morphological feature of the cell. To obtain a compressed image that keeps the motion feature, all frames of the binary video (b) were summed up over time, resulting in a trace image (f), which is then segmented into single cell track images (g). The compressed images for morphological and motion features are merged into a single compressed image (h). DL model and training The DL algorithm uses a convolutional neural network model38, which includes two hidden convolutional layers, two subsampling layers, a fully connected layer, and an output layer (Figs. 3 and S3). For each input image, the model produces an output of “1” if the cell is inhibited (or killed) by the antibiotic, or “0” if it is uninhibited. The model then determines the total number of uninhibited bacterial cells (NDL) over time for each antibiotic concentration (C), and produces

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

inhibition curves (NDL vs. C) from which the MIC value is determined, describing the minimal concentration of an antibiotic that inhibits the bacterial strain.

Figure 3 DL algorithm for antimicrobial susceptibility testing. (a) The DL workflow for determining MIC values. (b) Preparation of the DL training dataset. (c) Accuracy of the DL model in recognizing the effect of five antibiotics on each bacterial cell. The DL model produces an output for each input image in the data, and the accuracy was determined as the percentage that DL output coincided with the original label in (b). (d) The improved accuracy of the DL model with increasing training data size for AST with ampicillin. The error bar is the standard deviation of accuracy in 30 individual training runs.

Prior to AST, we trained the DL model to learn to differentiate antibiotic effects on bacterial cells with the input data (images). A practical difficulty in many DL applications is the preparation of a large training dataset to train the model.39 We overcame this difficulty by recording two sets of videos, each containing more than 100 bacterial cells. One video was for bacterial samples not exposed to antibiotics, and the second one was for samples treated with concentrated antibiotics (4X higher than the MIC value) for 6 hours to ensure inhibition of all the bacterial cells. The bacterial cells in the former video were assumed to be and labeled as “uninhibited”, and those in the latter as “inhibited” (Figure 3b) by the antibiotics. Both videos

ACS Paragon Plus Environment

Page 12 of 23

Page 13 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

were segmented into multiple sub-videos (1 s duration). We then compressed the sub-videos into single cell static images using the procedure described earlier and introduced them into the DL model for training. Figure 3c shows the accuracy of DL models for the five antibiotics. The training accuracies for the antibiotics are 97.5% for PMB, 98.5% for streptomycin, 98.0% for ciprofloxacin, 82.7% for aztreonam, and 82.3% for ampicillin, respectively. The training for each antibiotic lasted about 15 minutes using an office desktop computer (see Methods for details). A unique advantage of the DLVM-AST method is the improving accuracy with increasing size of training data sets. For example, Figure 3d shows an accuracy of 73% ± 3.1% for ampicillin if training the DL model with 500 single cell images. This accuracy increases to 87% ± 1.2% after testing additional 1500 datasets. This improving capability is particularly attractive if the DLVM-AST method is widely adopted and vast amount of validated test data become available to train and improve the DL model. Rapid AST with DL After DL training, we performed DLVM-AST for E. coli against five different antibiotics and validated the results with the gold standard BMD method. We recorded 30 s-interval videos of bacterial cells in the presence of antibiotics at a given concentration every 30 minutes over 3 hours, and then repeated these assays with different antibiotic concentrations. Figure 4a shows typical snapshots of E. coli (bright or dark rods) captured at 0, 30 and 180 minutes after exposing the sample to PMB at different concentrations, where red dashed circles marked bacterial cells were determined by the DL model to be ‘uninhibited’. We then segmented the videos into 30 sub-videos (1 s duration) and compressed each sub-video using the procedure described earlier, and introduced the single cell static images into the pre-trained DL model to determine the

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

number of uninhibited cells (NDL) in each sub-video. An inhibition curve was then plotted for each PMB concentrations, as the change of bacteria number (NDL) (normalized by the initial bacteria number in the microscopic field of view) over assay time (Figure 4b). If we define the empirical MIC as the antibiotic concentration at which less than 75% of the bacterial cells are uninhibited (dark dashed line in Figure4d) at 30 minutes after antibiotic treatment, the readout of the MIC in Figure4d is at 2 µg/mL. This MIC agrees with that obtained with the gold standard BMD method performed overnight (Table 1).

ACS Paragon Plus Environment

Page 14 of 23

Page 15 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Figure 4 AST results with PMB on E. coli. (a) Typical microscopic images recorded in the presence of PMB at different concentrations and after different treatment times. Bacterial cells appear as dark and bright spots in the images, while those cells marked with red dashed circles were recognized by the DL model to be phenotypically uninhibited upon exposure to PMB. (b) Inhibition curve plotted as the change of uninhibited bacterial number determined by the DL versus treatment time at different PMB concentrations. (c) Inhibition curve plotted as the change of total bacterial number determined by the division feature without DL as in Figure3 versus treatment time at different PMB concentrations. (d) Inhibition curves of DLVM-AST (b) after PMB exposure for 30 minutes. The bacterial number is normalized by the initial bacteria numbers in the

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

microscopic field of view. Scale bar in the figure is 5 µm. Error bars are standard deviation of bacterial numbers in 30 sub-videos.

To further validate the DLVM-AST method, we counted the total number of bacteria cells (Ntotal) from the videos (Figure 2), and plotted the change of Ntotal (mean ± std) vs. time for PMB at various concentrations (Figure 4c). The plots show that Ntotal increased over time (reflecting the growth and division of the bacterial cells) at low concentration PMB, but changed little over time (indicating effective inhibition) at concentrated PMB. We thus determined a MIC value of ~2 µg/mL (red curve) at 180 minutes as the minimal concentration that inhibits the increase of total bacterial number. This MIC value is also consistent with that obtained by DLVM-AST at 30 minutes. Using a similar procedure, we performed DLVM-AST with streptomycin, ciprofloxacin, aztreonam, and ampicillin on E. coli, each test repeated three times on different days. We plotted the MIC curves in Figure5 for all antibiotics determined from videos recorded 30 minutes after treatment, similar to that of Figure4d. The MIC results determined from the curves are summarized in Table 1, showing that the MIC values obtained by DLVM-AST are consistent with those by the gold standard BMD method for all the antibiotics. We also show the MIC values obtained by AST without DL as in Figure4c in 3 hours in Table 1.

ACS Paragon Plus Environment

Page 16 of 23

Page 17 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Figure 5 Inhibition curves obtained with DLVM-AST for four antibiotics on E. coli. (a) streptomycin, (b) ciprofloxacin, (c) aztreonam and (d) ampicillin. The MIC value is determined by the minimum concentration that less than 75% of the bacteria (dashed lines) remain uninhibited as determined by DLVM-AST.

Table 1 The MIC values (µg/mL) and assay time (h) determined by the three methods for all five antibiotics, with each test performed in triplicate. Antibiotics

MIC in 30 min (µ µg/mL) DLVM-AST

MIC in 3 hours (µ µg/mL) AST Without DL*

BMD**

Polymyxin B

2

1-2

2

Streptomycin

4

4

4

Ciprofloxacin

0.03-0.06

0.015-0.03

0.03-0.06

Aztreonam

0.12

0.12

0.12

ACS Paragon Plus Environment

MIC (µ µg/mL)

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Ampicillin

2

Page 18 of 23

2

2

*AST without DL: AST determined by division feature only (total cell number). **BMD: Broth Macrodilution.

Clinical testing To evaluate the effectiveness of the proposed approach, we then performed DLVM-AST with clinically collected urine samples. Mid-stream urine samples from 10 patients with UTI symptoms were obtained at the Mayo Clinic Hospital, in Phoenix, AZ, USA. Before AST, we first evaluated bacterial infection in each sample by automatic cell counting in the phase contrast microscopic images. Among the 10 samples, 1 case was confirmed positive for bacterial infection, and the other 9 cases were confirmed negative since few bacterial cells were observed. This evaluation corresponds well to laboratory results from isolation culture performed separately, and the pathogen in the positive case was identified to be E. coli, the most common pathogens in UTI. We then performed AST to determine the MIC values of five antibiotics against the pathogens in the urine samples, by DLVM-AST within 30 minutes, and standard BMD methods in12 hours. Figure6 shows the inhibition curves determined by DLVM-AST at 30 minutes after antibiotics treatments as described above. Table 2 shows the MIC values and susceptibility (according to CLSI QC range35) as determined by 30 minutes DLVM-AST and BMD methods. The DLVMAST determined the same MIC values as the BMD methods for antibiotics including ciprofloxacin, aztreonam, and ampicillin. For PMB and streptomycin, the MIC differs a little between the two methods, but indicated same susceptibility results. An important reason that the present DLVM-AST method is faster than the BMD and traditional microscopy-based AST methods based on cell counting without DL is its inclusion of multiple phenotypic features and analysis of the features at the pixel-level. Additional phenotypic features

ACS Paragon Plus Environment

Page 19 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

could also be included in the cell counting method based on the traditional image processing and classification techniques without using DL.14 However, defining and quantifying various phenotypic features, such as swimming, tumbling, filament formation, swelling and morphology changes, are difficult. DLVM-AST offers a universal solution to learn one or a combination of features from the videos without specifically defining and quantifying each specific feature. We note that it takes effort to train the DLVM-AST model, but once trained, the model can be used anytime afterwards and quickly (~5 minutes) predict MIC values from the videos. In fact, we pre-trained the DL model, and applied it to perform AST 4 months later to produce data shown in Figs. 4, 5, 6, and Tables 1, 2. Since the current work mainly focused on E. coli, the most common pathogens in UTI, validation with other pathogens and more clinical samples are needed to further apply the approach in the clinic. We also note that the present work focuses on videos, but the DL model is not limited to phenotypic features captured in the videos, and could be expanded to include biochemical features, such as adenosine triphosphate (ATP) consumption, proteins and nucleic acids9,20,23, to further improve its specificity and sensitivity and shorten the assay time.

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 23

Figure6 DLVM-AST results of the five antibiotics against the pathogens in clinical sample, including (a) polymyxin B, (b) streptomyxin, (c) ciprofloxacin, (d) aztreonam, (e) ampicillin. The results are obtained 30 minutes after antibiotics treatment. Table 2 The susceptibility and MIC values (µg/mL) determination reports of clinical urine samples (caused by E.coli) for all five antibiotics, with each test performed in triplicate. Antimicrobials CLSI QC range35* BMD** DLVM-AST in 30 (mg/mL) overnight minutes S I R MIC Susceptibi MIC Susceptibili Result lity*** Results ty*** (µg/mL) (µg/mL) Polymyxin B 0.25 S 0.5 S 4 ≤2 ≥8 Streptomycin 32/ >64 >64 — — — Ciprofloxacin 16 R 16 R 2 ≤1 ≥4 Aztreonam Ampicillin

≤4 ≤8

8 16

≥16 ≥32

≤0.125 16

S I

≤0.125 16

S I

*

CLSI susceptibility interpretation for Enterobacteriaceae. BMD: Broth Macrodilution. *** S: Susceptible. I: Intermediate. R: Resistant. **

Conclusion We have developed a deep learning video microscopy-based AST method (DLVM-AST). The video microscopy images single bacterial cells in urine without immobilization, which simplifies

ACS Paragon Plus Environment

Page 21 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

sample preparation, and allows capturing of the motion of each bacterial cell in free solution (urine), along with other phenotypic features, such as cell morphology and division. The deep learning model learns to differentiate bacterial cells inhibited by an antibiotic from those uninhibited cells from the videos at pixel level without defining, extracting and quantifying each specific phenotypic feature, thus overcoming the difficulty of the traditional imaging processing methods for analyzing complex features of bacterial cells. Using DLVM-AST we have rapidly determined antimicrobial susceptibilities of E. coli, the most common cause of UTI, against five different commonly used antibiotics. Our results show that DLVM-AST can accurately determine the minimum inhibitory concentrations (MIC) to inhibit E. coli within 30 minutes, which is faster than the gold standard broth macrodilution method and traditional microscopybased (cell counting) methods. Supporting Information Figure S1 Pre-processing the bacterial videos to improve contrast and remove noise in the images. Figure S2 Phenotypic features in the recorded bacteria videos in the presence and absence of antibiotics. FigureS3 Convolutional neural network model. Figure S4 Experimental setup for recording bacterial videos. Movie S1. Microscopic videos of freely moving bacteria cells spiked in urine sample.

References (1) Hancock, R. E. The end of an era?, Nature Reviews Drug Discovery 2007, 6, 28-28. (2) Neu, H. C. The crisis in antibiotic resistance, Science 1992, 257, 1064-1074. (3) Rossolini, G. M.; Arena, F.; Pecile, P.; Pollini, S. Update on the antibiotic resistance crisis, Current Opinion in Pharmacology 2014, 18, 56-60. (4) O’Neill, J. Tackling drug-resistant infections globally: final report and recommendations, The review on antimicrobial resistance 2016. (5) Dalgaard, P.; Ross, T.; Kamperman, L.; Neumeyer, K.; McMeekin, T. A. Estimation of bacterial growth rates from turbidimetric and viable count data, International journal of food microbiology 1994, 23, 391-404. (6) Reller, L. B.; Weinstein, M.; Jorgensen, J. H.; Ferraro, M. J. Antimicrobial susceptibility testing: a review of general principles and contemporary practices, Clinical infectious diseases 2009, 49, 1749-1755. (7) Wiegand, I.; Hilpert, K.; Hancock, R. E. Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances, Nature protocols 2008, 3, 163-175.

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(8) O’Neill, J. Rapid diagnostics: stopping unnecessary use of antibiotics, Review on Antimicrobial Resistance 2015. (9) Davenport, M.; Mach, K. E.; Shortliffe, L. M. D.; Banaei, N.; Wang, T.-H.; Liao, J. C. New and developing diagnostic technologies for urinary tract infections, Nat Rev Urol 2017, 14, 296-310. (10) Bergeron, M. G.; Ouellette, M. Preventing antibiotic resistance through rapid genotypic identification of bacteria and of their antibiotic resistance genes in the clinical microbiology laboratory, Journal of clinical microbiology 1998, 36, 2169-2172. (11) Dutka-Malen, S.; Evers, S.; Courvalin, P. Detection of glycopeptide resistance genotypes and identification to the species level of clinically relevant enterococci by PCR, Journal of clinical microbiology 1995, 33, 24-27. (12) Palmer, A. C.; Kishony, R. Understanding, predicting and manipulating the genotypic evolution of antibiotic resistance, Nature Reviews Genetics 2013, 14, 243-248. (13) Choi, J.; Jeong, H. Y.; Lee, G. Y.; Han, S.; Han, S.; Jin, B.; Lim, T.; Kim, S.; Kim, D. Y.; Kim, H. C.; Kim, E.C.; Song, S. H.; Kim, T. S.; Kwon, S. Direct, rapid antimicrobial susceptibility test from positive blood cultures based on microscopic imaging analysis, Scientific Reports 2017, 7, 1148. (14) Choi, J.; Yoo, J.; Lee, M.; Kim, E.-G.; Lee, J. S.; Lee, S.; Joo, S.; Song, S. H.; Kim, E.-C.; Lee, J. C. A rapid antimicrobial susceptibility test based on single-cell morphological analysis, Science translational medicine 2014, 6, 267ra174-267ra174. (15) Longo, G.; Alonso-Sarduy, L.; Rio, L. M.; Bizzini, A.; Trampuz, A.; Notz, J.; Dietler, G.; Kasas, S. Rapid detection of bacterial resistance to antibiotics using AFM cantilevers as nanomechanical sensors, Nature nanotechnology 2013, 8, 522-526. (16) Lissandrello, C.; Inci, F.; Francom, M.; Paul, M.; Demirci, U.; Ekinci, K. Nanomechanical motion of Escherichia coli adhered to a surface, Applied physics letters 2014, 105, 113701. (17) Syal, K.; Iriya, R.; Yang, Y.; Yu, H.; Wang, S.; Haydel, S. E.; Chen, H.-Y.; Tao, N. Antimicrobial susceptibility test with plasmonic imaging and tracking of single bacterial motions on nanometer scale, ACS nano 2015, 10, 845852. (18) Kinnunen, P.; Sinn, I.; McNaughton, B. H.; Newton, D. W.; Burns, M. A.; Kopelman, R. Monitoring the growth and drug susceptibility of individual bacteria using asynchronous magnetic bead rotation sensors, Biosensors and Bioelectronics 2011, 26, 2751-2755. (19) Sinn, I.; Albertson, T.; Kinnunen, P.; Breslauer, D. N.; McNaughton, B. H.; Burns, M. A.; Kopelman, R. Asynchronous magnetic bead rotation microviscometer for rapid, sensitive, and label-free studies of bacterial growth and drug sensitivity, Analytical chemistry 2012, 84, 5250-5256. (20) Liu, T.; Lu, Y.; Gau, V.; Liao, J. C.; Wong, P. K. Rapid Antimicrobial Susceptibility Testing with Electrokinetics Enhanced Biosensors for Diagnosis of Acute Bacterial Infections, Annals of Biomedical Engineering 2014, 42, 2314-2321. (21) Mann, T. S.; Mikkelsen, S. R. Antibiotic susceptibility testing at a screen-printed carbon electrode array, Analytical chemistry 2008, 80, 843-848. (22) Ertl, P.; Robello, E.; Battaglini, F.; Mikkelsen, S. R. Rapid Antibiotic Susceptibility Testing via Electrochemical Measurement of Ferricyanide Reduction by Escherichia c oli and Clostridium s porogenes, Analytical chemistry 2000, 72, 4957-4964. (23) Altobelli, E.; Mohan, R.; Mach, K. E.; Sin, M. L. Y.; Anikst, V.; Buscarini, M.; Wong, P. K.; Gau, V.; Banaei, N.; Liao, J. C. Integrated Biosensor Assay for Rapid Uropathogen Identification and Phenotypic Antimicrobial Susceptibility Testing, European Urology Focus 2016. (24) Baltekin, Ö.; Boucharin, A.; Tano, E.; Andersson, D. I.; Elf, J. Antibiotic susceptibility testing in less than 30 min using direct single-cell imaging, Proceedings of the National Academy of Sciences 2017, 114, 9170-9175. (25) Frymier, P. D.; Ford, R. M.; Berg, H. C.; Cummings, P. T. Three-dimensional tracking of motile bacteria near a solid planar surface, Proceedings of the National Academy of Sciences 1995, 92, 6195-6199. (26) Van Valen, D. A.; Kudo, T.; Lane, K. M.; Macklin, D. N.; Quach, N. T.; DeFelice, M. M.; Maayan, I.; Tanouchi, Y.; Ashley, E. A.; Covert, M. W. Deep Learning Automates the Quantitative Analysis of Individual Cells in Live-Cell Imaging Experiments, PLOS Computational Biology 2016, 12, e1005177. (27) Chen, C. L.; Mahjoubfar, A.; Tai, L.-C.; Blaby, I. K.; Huang, A.; Niazi, K. R.; Jalali, B. Deep Learning in Label-free Cell Classification, Scientific Reports 2016, 6, 21471. (28) Jean B., P.; Miller, L. A.; Franklin R., C.; Nicolau, D. P.; Bradford, P. A.; Powell, M.; Eliopoulos, G. M.; Swenson, J. M.; Hindler, J. A.; Traczewski, M. M.; Jenkins, S. G.; Turnidge, J. D.; Lewis, J. S.; Weistein, M. P.; Limbago, B.; Zimmer, B. L. Methods for Dilution Antimicrobial Susceptibility Tests for Bacteria That Grow Aerobically; Approved Standard-Tenth Edition, 2015. (29) Ying, S.-Y. Generation of cDNA Libraries, Methods in Molecular Biology 2003, 221.

ACS Paragon Plus Environment

Page 22 of 23

Page 23 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

(30) Lee, C.-C.; Sui, G.; Elizarov, A.; Shu, C. J.; Shin, Y.-S.; Dooley, A. N.; Huang, J.; Daridon, A.; Wyatt, P.; Stout, D.; Kolb, H. C.; Witte, O. N.; Satyamurthy, N.; Heath, J. R.; Phelps, M. E.; Quake, S. R.; Tseng, H.-R. Multistep Synthesis of a Radiolabeled Imaging Probe Using Integrated Microfluidics, Science 2005, 310, 1793. (31) Jing, W.; Jiang, X.; Zhao, W.; Liu, S.; Cheng, X.; Sui, G. Microfluidic platform for direct capture and analysis of airborne Mycobacterium tuberculosis, Analytical chemistry 2014, 86, 5815-5821. (32) Jing, W.; Zhao, W.; Liu, S.; Li, L.; Tsai, C.-T.; Fan, X.; Wu, W.; Li, J.; Yang, X.; Sui, G. Microfluidic device for efficient airborne bacteria capture and enrichment, Analytical chemistry 2013, 85, 5255-5262. (33) Lauga, E.; DiLuzio, W. R.; Whitesides, G. M.; Stone, H. A. Swimming in circles: motion of bacteria near solid boundaries, Biophysical journal 2006, 90, 400-412. (34) Sokolov, A.; Aranson, I. S. Physical properties of collective motion in suspensions of bacteria, Physical review letters 2012, 109, 248109. (35) Yue-Hei Ng, J.; Hausknecht, M.; Vijayanarasimhan, S.; Vinyals, O.; Monga, R.; Toderici, G. In Proceedings of the IEEE conference on computer vision and pattern recognition, 2015, pp 4694-4702. (36) Tran, D.; Bourdev, L.; Fergus, R.; Torresani, L.; Paluri, M. In Proceedings of the IEEE International Conference on Computer Vision, 2015, pp 4489-4497. (37) Simonyan, K.; Zisserman, A. In Advances in neural information processing systems, 2014, pp 568-576. (38) Krizhevsky, A.; Sutskever, I.; Hinton, G. E. In Advances in neural information processing systems, 2012, pp 1097-1105. (39) Chen, X.-W.; Lin, X. Big data deep learning: challenges and perspectives, IEEE Access 2014, 2, 514-525.

Acknowledgments: Financial support from the Moore Foundation is acknowledged. H.Y. W.J. and N.T. wrote the paper. N.T. conceived and supervised the project. W.J. and M.M. performed the AST experiments, and H.Y. and R.I. developed the DL method and analyzed the data. Y.Y., K.S., S.E.H., T.E.G., S.W. provided technical support on the project. All authors edited the paper. The authors declare no competing financial interests. Table of Content

ACS Paragon Plus Environment