Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Rapid identification of bacterial isolates using microfluidic adaptive channels and multiplexed fluorescence microscopy

Stelios Chatzimichail *ab, Piers Turner ab, Conor Feehily c, Alison Farrar ab, Derrick Crook cd, Monique Andersson cd, Sarah Oakley cd, Lucinda Barrett cd, Hafez El Sayyed ab, Jingwen Kyropoulos a, Christoffer Nellåker e, Nicole Stoesser cd and Achillefs N. Kapanidis *ab
aDepartment of Physics, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK. E-mail: stelios.chatzimichail@physics.ox.ac.uk; kapanidis@physics.ox.ac.uk
bKavli Institute for Nanoscience Discovery, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
cNuffield Department of Medicine, University of Oxford, John Radcliffe Hospital, Oxford, OX3 9DU, UK
dDepartment of Microbiology and Infectious Diseases, Oxford University Hospitals NHS Foundation Trust, Oxford, OX3 9DU, UK
eNuffield Department of Women's & Reproductive Health, University of Oxford, Big Data Institute, Oxford, OX3 7LF, UK

Received 13th April 2024 , Accepted 4th September 2024

First published on 6th September 2024


Abstract

We demonstrate the rapid capture, enrichment, and identification of bacterial pathogens using Adaptive Channel Bacterial Capture (ACBC) devices. Using controlled tuning of device backpressure in polydimethylsiloxane (PDMS) devices, we enable the controlled formation of capture regions capable of trapping bacteria from low cell density samples with near 100% capture efficiency. The technical demands to prepare such devices are much lower compared to conventional methods for bacterial trapping and can be achieved with simple benchtop fabrication methods. We demonstrate the capture and identification of seven species of bacteria with bacterial concentrations lower than 1000 cells per mL, including common Gram-negative and Gram-positive pathogens such as Escherichia coli and Staphylococcus aureus. We further demonstrate that species identification of the trapped bacteria can be undertaken in the order of one-hour using multiplexed 16S rRNA-FISH with identification accuracies of 70–98% with unsupervised classification methods across 7 species of bacteria. Finally, by using the bacterial capture capabilities of the ACBC chip with an ultra-rapid antimicrobial susceptibility testing method employing fluorescence imaging and convolutional neural network (CNN) classification, we demonstrate that we can use the ACBC chip as an imaging flow cytometer that can predict the antibiotic susceptibility of E. coli cells after identification.


Introduction

Antibiotics are critical to the treatment of bacterial infections, but their widespread use has contributed to the emergence of bacterial strains that can survive antibiotic treatment. These antibiotic-resistant bacteria pose a serious public health threat, causing more than 1.27 million deaths per year worldwide,1 and threaten a return to the pre-antibiotic era.

One of the means of controlling antibiotic resistance is by improving the accuracy and turnaround of diagnostic tests that detect the presence of bacteria and define their susceptibility to bacterial treatment directly from clinical samples. Culture-based tests for bacterial identification and antimicrobial resistance (AMR) already exist, but typically require 1–2 days to complete from the time the specimen is obtained from the patient. During that time, and especially for severe cases of infection, healthcare professionals prescribe broad-spectrum antibiotics based mainly on epidemiological data, when these are available. To address this, there is a clear need for rapid diagnostic assays that can be applied directly to clinical samples (i.e., without the need to first prepare clinical isolates, i.e., pure, monomicrobial cultures) to determine the type of infecting bacteria and enable accurate and rapid antibiotic susceptibility testing to guide patient treatment and antimicrobial stewardship.

Advanced bacterial identification methods include MALDI-TOF MS based methods which employ proteomic profile matching against known databases to rapidly identify bacteria from as little as a single-colony of material. Present state-of-the-art MALDI-TOF protocols such as the BD Rapid Sepsityper® only take 10 minutes to process a positive blood sample. However the key bottleneck with respect to identification time still remains the dependency on initial culture which has a time-to-positivity (TTP) of approximately 24 hours.2 Protocols aiming to optimize for TTP have still been found to require approximately 4–5 hours before MALDI-TOF testing, suggesting that methods capable of bypassing the bacterial culture step altogether could drastically reduce turnaround times.3 Furthermore, it should be noted that systems which rely on ensemble measurement of cells, often face challenges when dealing with polymicrobial infections as often only one of the species present in the infection may be identified.4

A promising approach for rapid bacterial identification and analysis from clinical samples involves microscopy techniques. While detection of bacterial pathogens using microscopy on clinical isolates is relatively simple, detection in complex clinical samples is much more challenging, since, in addition to the need to distinguish bacteria from the sample matrix, some clinical samples have low bacterial concentrations in the absence of pre-culture steps (e.g., 1–10 colony-forming units (CFUs) per mL for blood).5,6 To enable the study and quantification of bacteria by microscopy, it is thus important that the target bacteria are enriched and spatially concentrated so as to allow rapid identification and interrogation of the pathogen.4

To achieve bacterial enrichment, researchers have previously used a combination of size-selection modules such as nano-porous monoliths,7 hydrodynamic traps,8 and nanostructured channels,9 as well as bacterial surface affinity agents, such as apolipoprotein-H and mannose-binding lectin.10,11 Further distinguishing features, such as small differences in bacterial length and morphology, have also been leveraged in combination with viscoelastic focusing to enrich bacterial sub-populations.12 Most affinity-based isolation techniques, however, lend themselves more to genotypic methods (such as qPCR) as opposed to microscopy-based methods, mainly due to the use of beads which can preclude microscopic observation due to substantial loss of spatial information and loss of population heterogeneity information from the sample.

In contrast, microfluidic platforms have been increasingly used to enable bacterial enrichment compatible with microscopy. These platforms typically employ size-selection modules in the form of hydrodynamic traps in order to capture or slow down bacterial cells so as to render them quasi-static for microscopic investigation.8,13,14 For instance, Stratz et al.13 made use of hydrodynamic traps with 1 μm gaps to capture E. coli for single-cell studies. In their setup the authors employed PDMS elastomeric valves to create enclosed chambers following trapping and lysis of the bacteria and subsequent capture of their contents by a pulldown assay. However the above setup, while optimized for the formation of single-cell arrays of E. coli for proteomic analysis, was not necessarily tuned for high capture-efficiency which is necessary when processing low cell-density samples.

Adaptive channels are another way to render bacteria static for microscopic imaging and often comprise a microchannel underlayer where the sample flows through and a control overlayer which compresses and expands so as to tune the dimensions of the sample flow-channel. Such adaptive channels have previously been employed in tandem with nanostructured channels15 to enable the quasi-static immobilization of bacteria for mechanical phenotyping as well as to introduce bacteria to ‘mother-machine’ like nanostructures to enable growth-rate monitoring of bacteria under different conditions.16 Okumus et al.14 used adaptive channels in an imaging flow cytometer fashion by flowing E. coli cells to a field-of-view, slowing them down by semi-actuating the control layer, and finally imaging cells by fully closing the channel and rendering cells quasi-static. Cells were finally released, and this approach was repeated iteratively to exchange with new cells to facilitate the counting of low-abundance proteins in cells in a high-throughput manner.

Adaptive channels circumvent a common issue associated with conventional bacterial size-based capture platforms which are the sub-micron feature dimensions required to achieve bacterial capture. Sub-micron channel regions result in very high hydrodynamic resistances which, given the ∼5 bar burst pressure of conventional polydimethylsiloxane (PDMS)–glass devices, result in operational flow rates of only a few μLs per minute, and in turn drastically reduce the maximum sample processing speed.17 Although this limitation does not affect bacterial capture significantly (as parallel channels can be incorporated in microfluidic architectures to increase the speed at which the biofluid sample can be processed), it does limit downstream cell assays requiring multiple reagent infusions, e.g., sequential fluorescence in situ hybridization (FISH) techniques such as multiplexed error-robust fluorescence in situ hybridization (MERFISH).18 Identifying species using FISH by targetting hypervariable regions of the 16S-rRNA is a common rapid technique to ascertain the identity of bacterial species down to the strain level. In recent years, 16S-rRNA FISH techniques have also become more scalable via barcoding strategies (akin to those in MERFISH), an example of which is high phylogenetic resolution microbiome mapping by fluorescence in situ hybridization (HiPRFISH)19 which enables the assay of panels containing thousand of bacterial species.

Barcoded FISH techniques are particularly attractive for the fluorescence-based identification of messenger RNA (mRNA) in spatial gene expression studies and of ribosomal RNA (rRNA) in bacterial identification studies. These techniques address the limitation of target scalability as the panel of identifiable targets may scale combinatorially with respect to the number of reagent infusion rounds performed. Shi et al.19 and Kandavalli et al.20 overcame the aforementioned sample processing speed issue by performing a spectral analogue of these techniques by using multiple complementary imager DNA strands bearing different fluorophores. For the latter study, this enabled the dye-conjugated oligonucleotide strands attached to the 16S ribosomal subunit of collected bacteria to form a spectral barcode and thus enabled combinatorial-FISH of up to 4 different imager probes leading to the identification of 7 species while minimizing the number of infusion reagent rounds. Scaling of the bacterial assay panel beyond 10 pathogens, however, would nonetheless necessitate more infusion reagents and would therefore slow down bacterial identification on this platform as a result of long reagent exchange times.20

Here, we achieve bacterial enrichment and identification by multiplexed 16S rRNA-FISH in a streamlined pipeline (Fig. 1) by employing an adaptive channel to capture bacteria. This is achieved by robotic control and tuning of the channel backpressure using an in-line pressure sensor and an active feedback control system. Our approach enables the controlled formation of a region within the flow-channel with highly reduced dimensions which we refer to as the ‘capture region’. The capture region acts as a bacterial filtration element that can be formed or removed on-demand. Once bacteria are captured, the adaptive channel can be restored back to dimensions that enable high flow-rate infusion of reagents, making such platforms ideal for sequential assays. We prepare our devices using 3D-printed moulds, thus circumventing the need for complex and expensive fabrication techniques that would otherwise be needed to achieve sub-micron features. Finally, we demonstrate the identification of a panel of 7 species by evaluating simulated samples containing spikes of fixed isolates in aqueous buffer. We achieve this in the order of one-hour by using our device in tandem with multiplexed 16S rRNA-FISH. Using relative fluorescence intensity measurements for bacterial species classification, we show that our technique achieves high accuracy (in most cases, >90%) in species identification. In addition to the technique being rapid and scalable, we demonstrate that by identifying bacteria at the single-cell level, we can also discern the presence of mixed infections. Our work paves the way for using such fluidic devices for bacterial enrichment and identification in complex clinical specimens.


image file: d4lc00325j-f1.tif
Fig. 1 Illustration of the workflow for the bacterial capture and identification process from bacterial samples. Using the ACBC device, bacterial isolate sample is flown and captured using dynamically formed capture regions. The immobilised bacteria are then stained and assayed using multiplexed 16S rRNA FISH. The single-cell intensity of the imager probes is binarized so as to form a barcode which is then compared against a species dictionary to yield the identity of the species.

Materials and methods

Preparation of ACBC devices

Flow and control channel moulds of microfluidic devices were designed using Autodesk Fusion360 and printed with an LCD stereolithography 3D printer (Anycubic Mono 4k, Shenzhen, China) using Anycubic ECO plant-based clear resin. The moulds were washed with IPA, further UV cured for 15 minutes and placed in an oven at 70 °C for two days to eliminate uncured resin residues that have otherwise been found to inhibit PDMS crosslinking.

For the flow-layer, PDMS (Sylgard 184, Dow Corning, United States) mixed at a ratio of 20[thin space (1/6-em)]:[thin space (1/6-em)]1 was spin-coated (Ossila spin-coater, Ossila Ltd, The Netherlands) at 500 RPM on the flow-layer mould. For the control layer, PDMS was mixed at a ratio of 5[thin space (1/6-em)]:[thin space (1/6-em)]1 on a 3D printed mould until 2 mm the height of PDMS reached 2 mm. The two moulds were then placed in an oven at 70 °C for 35 minutes. The off-ratio bonding technique was then employed. Briefly, the moulds were removed from the oven and allowed to cool down. The cured control layer PDMS negatives were then cut out, and inlets were punched using a biopsy punch (1 mm) and aligned to the flow-layer devices under a microscope. PDMS mixed at a ratio of 5[thin space (1/6-em)]:[thin space (1/6-em)]1 was then poured around the control layer to form a device with a final thickness of 2 mm and cured in an oven at 70 °C for 1 hour. The PDMS device was finally detached from its mould, inlets were punched using a biopsy punch (1 mm) and finally bonded to glass after the two parts were treated with air plasma. Devices were placed in an oven at 70 °C overnight to further reinforce bonding. The devices were finally interfaced using steel pins (gauge: 5G) fitted to silicone tubing (Tygon, ID: 250 μm ID).

For medium pressure interconnects, the interconnection was reinforced by using epoxy (Araldite fast-cure) at the steel–PDMS interface, allowing the epoxy to harden for 6 hours and finally an additional baking step at 70 °C for 2 hours.

Flow instrumentation

The fluidic setup consisted of a syringe pump (Elite 11, Harvard Apparatus, USA) operated in constant pressure mode used to infuse the sample as well as 8 custom-made Arduino controlled syringe pumps used to deliver assay reagents via a Labsmith 8-port selector valve (AV801). A 3-way selector valve (Labsmith, AV201) is connected to the main chip line and selects whether the sample containing syringe pump or the assay reagent valve are engaged for infusion to the chip. Prior to the ACBC device inlet, an in-line pressure sensor (LS-uPS0800) is placed to measure the ACBC device backpressure Pchip. For adaptive channel control, a variable constant air pressure source was connected to an in-line pressure sensor (LS-uPS0800) and interfaced to the control layer of the ACBC device.

Bacterial strains and preparation of fixed bacterial isolates

Bacterial strains used, with the exception of E. coli MG1655, originate from blood culture isolates processed for diagnostic purposes and stored by the Microbiology Laboratory of the Oxford University Hospitals NHS Foundation Trust, Oxford, UK. Individual colonies from these banked clinical isolates of E. coli, P. aeruginosa, K. pneumoniae (Gram-negative) and E. faecalis, S. pneumoniae, S. agalactiae, and S. aureus (Gram-positive) were grown in brain heart infusion (BHI) broth supplemented with 2% yeast extract overnight at 37 °C and subsequently diluted in fresh broth for further growth at 37 °C until an OD600 of 0.2 was reached. A 1 ml volume of this culture was then centrifuged at 10[thin space (1/6-em)]000 × g for 5 minutes, the supernatant was discarded, and the pelleted bacterial cells were then fixed with 1 ml of a PBS solution containing 2% (v/v) paraformaldehyde for 20 min. Cells were then re-centrifuged and the pellet resuspended and permeabilized with 1 ml 1[thin space (1/6-em)]:[thin space (1/6-em)]1 ethanol/PBS solution and stored at −20 °C until use.

Preparation of bacterial samples for capture studies

Bacterial capture studies for E. coli and S. aureus were conducted with fixed cells of strains ECOL1 and SAUR1 respectively (ESI) suspended in PBS. The densities of the cell suspensions were determined using a Neubauer haemocytometer and the suspensions of the strains were serially diluted in PBS buffer to obtain the desired cell densities for capture studies. For exemplar ACBC capture videos (ESI: Videos S1 and S2), the same protocol was followed, except the cells were also stained with Nile red membrane stain (15 minute incubation, 1 μg mL−1 in PBS) for cell circumference visualization.

Bacterial capture using ACBC device

For capture efficiency studies, devices were first surface-treated with chitosan (0.015% w/w in 0.1 M acetic acid, 200 μL) for 10 minutes and then washed with 500 μL of 1× PBS. Fixed cells were then flown, visualized using their native autofluorescence arising from the growth in rich medium (473 nm excitation), and enumerated.

The ACBC control layer, filled with Milli-Q water, was then actuated at 2 bar, after which bacteria containing sample was flown through the device at a flow rate of 40 μL min−1. The backpressure of the system was monitored throughout. The formation of a capture region was observed visually as backpressure built up within the device and the syringe pump controlled via custom software was then instructed to maintain constant pressure so as to not perturb the geometry of the capture-region. The collection of bacteria was then allowed to proceed for 20 minutes, or until a sufficient amount of target was deemed to have been collected.

Imaging

Fluorescence images were captured using a wide-field Nanoimager microscope (ONI, Oxford, UK) equipped with a Hamamatsu Flash4 v3 sCMOS camera. Samples were imaged in highly inclined thin illumination (HILO) mode using a 100× oil-immersion objective. The laser illumination was focused at 2° below TIR which was at an angle of 51.5° with respect to the normal. For multiplexed FISH assays, images were acquired at an exposure time of 33 ms using a 640 nm excitation laser at a power of 0.78 kW cm−2. For tracking of the bacteria during flow-through experiments, the sample was illuminated for an exposure time of 500 ms using a 473 nm laser at a power of 1.09 kW cm−2.

Bacterial segmentation

For bacterial segmentation we employed Cellpose,21 an instance segmentation model, to segment cells from epifluorescence images of E. coli (MG1655) labelled with EUB338-Cy3 probes excited with 532 nm light. To improve segmentation performance, we trained a custom Cellpose model on our microscopy data for 100 epochs using the standard Cellpose hyperparameters. Cellpose segmentation and segmentation curation/editing was carried out in Napari-BacSeg, which is a custom Napari plugin that was built for segmenting and analysing images of bacteria. Napari-BacSeg is available from the Napari Hub.

Multiplexed 16S rRNA-FISH

Fixed bacterial isolates were centrifuged at 7000g for 5 minutes, re-suspended in PBS and further permeabilized using a solution containing lysozyme (20 mg mL−1) and lysostaphin (0.1 mg mL−1) dissolved in TEG buffer (25 mM Tris, 10 mM EDTA, 50 mM glucose, pH = 8) for 20 minutes at 37 °C. The cells were then pelleted at 7000g and resuspended in 200 μL of PBS.

The bacteria were flown for 20-minutes through chitosan-surface treated ACBC chips and captured and immobilized at capture-region forming conditions (PQ = 2 bar, Pchip = 1 bar). The permeabilized bacteria were then encoded by flowing encoding hybridization buffer (1× Denhardt's reagent, 10% dextran sulfate, 2× SSC, 30% formamide, 1 mg mL−1 yeast t-RNA, 0.04% SDS, EP ssDNA probes) for 15 minutes (flow rate: 50 μL min−1, Vinfused = 50 μL) and then washed using encoding wash buffer (2× SSC, 30% formamide) for 1-minute (flow rate: 100 μL min−1, Vinfused = 100 μL). For Gram-negative species, 1 μM EP ssDNA probes were used, whereas for Gram-positive species 5 μM of EP ssDNA probes were used. The encoding probe sequences can be found in the ESI (Table S1).

Following the encoding process, the sample was then sequentially interrogated using imager probe (I1–I6) hybridization solutions (1× Denhardt's reagent, 10% dextran sulfate, 2× SSC, 10% formamide, 1 mg mL−1 yeast t-RNA, 0.04% SDS, 100 nM IP-Cy5 ssDNA probe) for 1 minute (flow rate: 50 μL min−1, Vinfused = 50 μL), followed by washing of weakly bound probes with imager wash buffer (40% formamide, 0.2× SSC) by infusion of the wash buffer (flow rate: 100 μL min−1, Vinfused = 100 μL). Following fluorescence imaging of the imager probe, the fields of view of interest were subsequently photobleached prior to the introduction of new imager probes. The imager probe sequences can be found in the ESI (Table S2). For validation of bacterial presence the EUB338-Cy5 probe (Table S2) was used (2× SSC, 10% formamide, 1 mg mL−1 yeast t-RNA, 0.04% SDS, 100 nM EUB probe).

Ultra rapid AST via ACBC imaging flow cytometry

For the AST studies, MG1655 E. coli cells were grown in LB broth until an OD600 of 0.2. The cells were then treated with either ciprofloxacin (0.5 mg L−1) or rifampicin (100 mg L−1) for 30 minutes. A 1 ml volume of this culture was then centrifuged at 7000 × g for 5 minutes, the supernatant was discarded, and the pelleted bacterial cells were then fixed with 1 ml of a PBS solution containing 2% (v/v) paraformaldehyde for 20 min. Cells were then re-centrifuged and the pellet resuspended in 1 mL PBS solution containing Nile red (1 μg mL−1) and DAPI (1 μg mL−1) for 5 minutes. Cells were re-centrifuged and the pellet resuspended in 1 mL of PBS solution. For untreated cells, the antibiotic treatment step was omitted.

The stained cells were diluted to a concentration of 1 × 105 cells per mL and flown through ACBC chips at capture conditions (PQ = 2 bar, Pchip = 1 bar) and imaged at 33 ms exposure using alternating pulses of the 405 nm laser (5.1 kW cm−2) for DAPI visualization and 532 nm laser (16.5 kW cm−2) for Nile red visualization while operating in HILO mode.

For phenotype prediction, images of cells were segmented and exported to create standardised 64 × 64 images of each cell zero-filled outside the segmentation boundary. Pre-trained convolutional neural networks (CNN) based on the DenseNet121 architecture for ciprofloxacin treated and rifampicin treated cells were then employed to predict the phenotype of the imaged cells. For more detailed explanations as to the training of these CNNs, training code and datasets as well as the pre-trained models, the interested readers can find these in ref. 22.

Results and discussion

ACBC chip design & operational principles

To address the limitations of previous approaches (especially the low operational flow-rates associated with sub-micron flow channel dimensions) while still being able to achieve bacterial enrichment, we reasoned that active formation and removal of bacterial-capture regions would facilitate both enrichment and rapid infusion of reagents after the initial capture stage of the identification pipeline.

We also designed our Adaptive Channel Bacterial Capture (ACBC) chip to accommodate the enrichment of samples containing a range of bacterial concentrations (102–104 cells per mL, with 102 cells per mL reflecting patient specimens with relatively low CFU counts). We thus designed a device that captured bacteria with near 100% efficiency and immobilized bacteria to a small number of fields of view. To achieve that, we designed a bilayer PDMS microfluidic device comprising of a top control layer and a bottom fluidic layer, akin to a conventional top-down pneumatic valve (Fig. S1).23,24 To achieve a pre-defined location of the capture region that was invariable to user alignment errors of the control layer, we utilized an LCD stereolithography 3D printer so as to produce channels with both parabolic and rectangular cross-sections present on the same mould (Fig. S2). In doing so, the capture region (Fig. 2) was defined physically by the intersection of the channels designed with rectangular and parabolic cross-sections rather than the edge of the control channel itself, which would be prone to miniscule placement errors. We found the minimal channel width we were able to produce consistently using this approach as measured by bright field microscopy to be 34.6 ± 1.6 μm (n = 3 moulds).


image file: d4lc00325j-f2.tif
Fig. 2 Device operation under different conditions of applied pressure PQ and lift pressure PL. Each set of conditions is depicted by a schematic as well as a position-map that showcases the end-position of E. coli cells flown from a PBS spike sample containing 7 × 102E. coli per mL across 2 runs. Dashed lines serve as a guide to the eye for the channel dimensions (white) and capture-region (green). A) The device at rest has its original design dimensions, no capture region is formed and thus no bacteria are trapped. B) The control layer is actuated at a pressure above the closing pressure PC. The channel dimensions collapse and no sample flow is observed. C) When a sufficient lift pressure is applied, flow through the device is enabled and a capture-region with sub-micron dimensions is formed allowing the trapping of bacteria. The end-position of flowing cells, as can be seen in the corresponding position-map, is constrained within this region. The device's chitosan coating enables the immobilization of cells within this narrow passage. D) At high device backpressure values of PL the capture region can be expanded before it finally leaks across the two sides of the channel. The experimentally obtained images used to construct these maps can be found in the ESI. Scale bars correspond to 10 μm.

To perform bacterial captures, the pneumatically actuated control layer acts on a parabolic flow-channel resulting in an initially fully-sealed channel (Fig. 2B). Deviating from conventional pneumatic valve principles, and by infusing sample through the device (Fig. 2C) at pressures exceeding the closing pressure of the valve, the compressed fluidic channel eventually opens due to deformation of PDMS from the backpressure build-up inside the flow channel. The opening of the channel is controlled using a closed feedback-loop that maintains the channel backpressure at the desired value. When we flow bacterial cells to the capture region of our devices, they are hydrodynamically trapped at the capture-region and finally immobilized locally on the polycationic chitosan coating of the glass substrate. Chitosan has previously been shown to be a good substrate for both the immobilization of bacteria as well as its good biocompatibility, having been shown to pose no significant effects on cell fitness.25

In our experiments, we used fixed cells instead of live cells to enhance the consistency of cell immobilization on the chitosan layer following trapping. This approach was necessary as certain strains of species such as K. pneumoniae were challenging to immobilize in their live form. We postulate this to be due to factors such as the presence of active flagella, which enabled these cells to maintain their mobility despite being partially attached to the chitosan surface. Cell fixation ensured a more reliable and reproducible immobilization process across different bacterial species and strains thereof. Under these conditions, immobilisation to the chitosan surface was stable in our system, allowing the ceiling to be retracted back to full channel scale without release of the captured cells. By restoring the channel dimensions post-capture, we circumvent the issue that common sub-micron structured devices have with regards to operational flow rates and thus enable rapid delivery of FISH reagents up to 50 μL per minute, drastically reducing reagent delivery times.

When attempting to use top-down valves as sieving elements in elastomer systems, it is important to note that the geometry of the microfluidic channel can alter as a function of channel backpressure due to the elastic nature of PDMS (Video S1).26 It is therefore not advised to operate syringe drivers under constant volumetric flow-rate when trying to achieve a static channel geometry, as would be necessary for the prolonged filtration of bacteria. Expansion of the channels due to pressure build-up would have the effect of releasing the captured bacteria over time. To achieve a quasi-static geometry, the syringe pumps used here operate under constant pressure control over the 20-minute infusion of the sample. Videos of PDMS deformation in the absence and presence of pressure control can be found in the electronic ESI (Videos S1 and S2).

To form a channel with suitable dimensions for bacterial filtration it is important to pressurize the control channel at a pressure PQ which is above the closing pressure of the valve (PC). The PC value can be experimentally determined by pressurizing the control layer while the flow channel is at rest and visually examining the point at which bright-field contrast between the flow channel boundary and the PDMS walls (Fig. S3b) vanishes. We found PC to be 1.3 ± 0.2 bar (n = 5 devices) for control layers fabricated using our method. This value exceeds devices made with soft lithography, which we attribute to the larger depth of our channels, as well as their partial deviations from semi-circularity. Compera et al.27 have shown that closing pressures of 0.25–1.3 bar were observed when using channel depths of 100 μm from 3D printed moulds which were re-flown or anti-aliased respectively, the latter of which are like the ones employed here.18

Using fluorescence measurements under no net-flow (Pchip = 0 bar) we estimated the channel's resulting height at different values of PQ (Fig. S3a) and observed a linear relationship between applied pressure PQ and the resulting channel height, h, with a regression coefficient of 16.6 μm bar−1 (r2 = 0.98).

By infusing the flow channel with a syringe pump able to generate sufficient torque to exceed the difference of the control pressure PQ and the device closing pressure PC, an overall lifting pressure can be generated that raises the collapsed channel. We calculate this lift pressure, PL, experimentally from our device's in-line pressure sensors using the following equation:

 
PL = Pchip − (PQPC)(1)

The flow-channel pressure, Pchip, is calculated live and used to estimate the height of the collapsed channel of the capture region. We estimated the channel's resulting height using fluorescence measurements at actuated control-layer conditions (PQ = 2 bar) for different values of Pchip by applying flow-induced pressure to the flow layer (Fig. S3b) and observed a linear dependence of resulting channel height, h, and the applied pressure. A regression coefficient of +8.6 μm bar−1 (r2 = 0.98) was determined which was lower than the magnitude of the membrane displacement coefficient observed when flow was at rest (−16.6 μm bar−1). Given that PDMS exhibits linear elastic behaviour for strains up to 20%,28 which would not be exceeded with a membrane displacement of 20 μm (membrane thickness: 200 μm), we presume that this difference is associated with the direction of pressure as applied by the flow-channel to the PDMS membrane. Indeed, when we pressurize the flow-channel, in addition to changes in channel height we also observe receding of the membrane along the flow-path direction. At our operational conditions (PQ = 2 bar, Pchip = 1 bar, PL = 0.3 bar) the stress along the flow-channel dimensions results in the planar projection of the capture region to be parabolic in shape with maximum displacement in the order of 25 μm as can be seen below (Fig. 2C). This suggests that stress is distributed both upwards but also longitudinally to the flow-path direction and could explain the difference in the regression coefficients observed in the two scenarios.

Under the above operational conditions, h was estimated to be 2.3 μm at the edge of the capture region. This estimate is larger than the typical diameter of bacteria, however in practice we observe near 100% capture efficiency (ηeff) of E. coli bacteria at these conditions. Fluorescence measurements as used here for microfluidic height estimation have high uncertainties at larger channel heights.29 This occurs because the emitted light is not fully collected throughout the entire depth of the channel, leading to an underestimation of the intensity corresponding to that height. Conversely, when the channel is in its collapsed state with a small height, all the emitted light from the excited molecules is captured. Consequently, when the adaptive channel is pried open from its closed state, it produces an overestimated height value as the resting state of the channel was used as the reference.

We found that operating at a pressure PQ of 2 bar or higher (which is ≫PC) is more practical, since operating at PQ = PC led to more pronounced effect of backpressure variations on the final geometry of the capture region. Further, a PQ of 2 bar did not damage the PDMS membrane over prolonged periods of actuation. The boundary of the capture region can be visually observed using bright field microscopy as well as by the autofluorescence of PDMS when excited at short wavelengths (405–473 nm) in highly inclined thin illumination (HILO) mode (Fig. S4). This is a useful feature as it enables placement tuning of the capture region as desired depending on the experiment at hand. For instance, we found that for samples with high bacterial densities (107 bacteria per ml) that the narrow region of our device would result in very dense packing of bacteria over time, in which case the capture region area could be expanded towards the wider section of the device by allowing the backpressure of the device to build up further (Fig. 2D). This in turn allowed for much more uniform spreading of the cells on the chitosan surface, simplifying bacterial imaging of these samples.

Bacterial capture and recovery

To assess the device's performance at capturing bacteria at low cell densities, we performed capture experiments with two species, E. coli and S. aureus, spiked at known concentrations in phosphate buffer saline (PBS) and subsequently flowing the cell suspension through our device over the course of 20-minute runs.

The progress of the capturing process was monitored using the autofluorescence of E. coli and S. aureus (λexc = 473 nm). Bacteria grown in rich media tend to exhibit high green autofluorescence due to the secreted flavins in the external media which is a useful feature for cell-counting purposes.30 Monitoring the trajectory of flowing bacteria using this method enabled both visualization of the trajectory of incoming cells and counting of the number of cells able to pass beyond the capture region.

Using combinations of PQ and Pchip such that a capture region is formed but not over-expanded (PQ = 2 bar, Pchip = 0.8–1.1 bar), we achieved ∼100% capture efficiency (see below) of flowing cells over the course of 20-minute flow-through experiments (Fig. 3b).


image file: d4lc00325j-f3.tif
Fig. 3 a) Autofluorescence timelapse of MG1655 E. coli cells collected in the ACBC capture region from a 7 × 102E. coli per mL sample over 20 minutes. Newly arriving cells display higher autofluorescence intensity relative to cells imaged over the course of the timelapse due to photobleaching of the innate cell fluorescence over prolonged exposure. Scale bar is 15 μm. b) Number of cells captured over the course of a 20-minute filtration for a sample with a cell count of 7 × 102E. coli per mL (n = 4) and 5 × 103E. coli per mL (n = 3). On average, 27 and 191 E. coli cells were captured by the end of 20-minute course runs respectively. The dashed black lines denote the theoretical maximum capture from samples at these concentrations. c) Exemplar pressure trace obtained during capture runs. With the adaptive ceiling actuated, sample is infused through the flow-layer resulting in an initial increase of Pchip. When capture region formation is observed (i), the syringe pump is instructed to maintain this value of backpressure. Initial overshooting to the target pressure can be seen before stabilization of the backpressure to the desired value (ii).

The maximum flow-rate through the device while producing backpressures able to sustain a capture region was found to be 2.1 ± 0.3 μL min−1. For dilute samples of E. coli, containing 7 × 102 cells per mL, 20-minute flow-through runs were found to yield 26.8 ± 3.1 captured cells (n = 4, ηeff: 96.4 ± 6.3%). At this processed sample volume (42 μL) and sample concentration, stochastic sampling effects are not expected to be and therefore the expectation value of 29 cells per run at these conditions is a reasonable estimate. The above number of cells captured is within error from the average value of 29 cells expected to be captured at 100% capture efficiency (p = 0.19 > 0.05). Furthermore, no fluorescence trajectories of cells bypassing the capture region were observed during flow-through acquisitions. These findings therefore suggest that there are only negligible losses of E. coli cells through the device and that the capture region is capable of achieving capture of E. coli at near 100% capture efficiency.

For the capture of S. aureus cells (the smallest pathogen in our panel), and at the lowest concentration tested (7 × 102 cells per mL) we observed a capture efficiency of ηeff = 71.4 ± 7.6% under the same experimental conditions. This is lower than the capture efficiency of E. coli and can be attributed to the smaller diameter of S. aureus cells, enabling a proportion of the population to bypass hydrodynamic trapping by the ACBC chip. Indeed, during flow-through runs of S. aureus cells we made note of cells that were able to escape the capture region. At the higher cell density tested (5 × 102S. aureus cells per mL), we observed higher capture efficiencies (ηeff = 90.3 ± 7.6%), something which we attribute to the keystone effect becoming more prevalent at higher cell densities.31

When the adaptive ceilings of the device were not actuated (e.g., conditions as in Fig. 2A), we did not observe the capture of any bacteria within the capture region observation field suggesting that the chitosan coating alone, used to immobilize bacteria post-capture, is not sufficient to isolate an appreciable number of cells at this cell density. Fig. 2 denotes the end position of E. coli bacteria captured across a series of runs for fully sealed capture regions (n = 2) and expanded capture regions (n = 2). As seen in the fully sealed case (Fig. 2C), the end-position of bacteria is localized within a narrow region and bacteria are attached uniformly within this area. When higher Pchip pressure is applied, in addition to the height increase of the membrane, the boundaries of the capture region are further displaced downstream (16 μm displacement across the boundary), effectively expanding the capture region. Notably, when no chitosan coating was used, the end positions of bacteria were observed predominantly at the boundaries of the capture region (ESI: Videos 2 and S3).

It was important to ensure that the ACBC device was able to capture a wide range of bacterial species to ensure that a variety of bacterial pathogens could be identified using this platform. Hydrodynamic trap formats act as low-pass filters and therefore the capture region of the ACBC device could potentially be unable to trap flowing bacteria of smaller dimensions. To verify the capabilities of the device to capture smaller diameter bacterial species we performed capture runs of a S. aureus (d = 0.5–1.5 μm) strain.32 We found that S. aureus cells are captured efficiently (ηeff = 71.4 ± 7.6%), suggesting that the ACBC chip is capable of capturing smaller diameter bacteria as well. A capture timelapse for S. aureus captures can be found in the ESI (Video S4).

Following microscopic investigation, trapped bacteria could be released from the chitosan coated surface of our device by flowing alkaline lysis buffer (0.2 M NaOH, 1% SDS) through the device. This could potentially enable further downstream analysis of the sampled cells to be carried out, e.g., for sequencing methods (Video S5).

Multiplexed 16S rRNA-FISH for bacterial identification

To demonstrate that the device can support processing and further analysis of the captured bacteria, we developed methods to identify a panel of seven pathogenic bacteria commonly associated with human infections (Fig. 4). These species contribute to many cases of serious human infection, including bloodstream and joint infections, and meningitis.
image file: d4lc00325j-f4.tif
Fig. 4 a) Schematic of the multiplexed 16S rRNA assay. The rRNA of the captured bacteria is hybridized with encoding probes containing two landing pads. Six rounds of hybridizations with fluorescently tagged oligos then reveal the identity of the bacterial species. b) Confusion matrix of this binary classification method showing classification rates for all 7 species assayed (averaged across 4 strains of each species) from bacterial isolates. The inconclusive class corresponds to any barcode raised that does not correspond to an entry in the species dictionary.

Specifically, we employed encoding oligonucleotides containing a 30-nucleotide (nt) targetting sequence, flanked by two 25-nt “coding sequences”. To design these sequences, 16S rRNA targetting sequences were first identified from literature studies that deployed these probes in human samples (Table S1); this was done to ensure there was no cross-interference between the probe and common host sequences present in clinical sample matrices. These targetting sequences were aligned against strains found in the NCBI's database of 16S ribosomal RNA sequences for bacteria and archaea33 and extended to 30-nt to ensure all sequences exhibited similar melting temperatures and more resiliency towards stringent probe washes. The targetting sequences were finally flanked on both sides by 25-nt coding sequences obtained from the Elledge group's database of orthogonal oligos34 while ensuring the assembled oligo showed no continuous alignment beyond 14 base pairs to either the 16S targetting sequences nor the genome of the panel species.

To perform the assay on captured bacteria, we performed an initial 20 min hybridization step containing all encoding probes, and then washed all weakly bound probes. The coding sequences were subsequently hybridized and imaged using complementary fluorescent Cy5-imager probes which were flown sequentially through the ACBC chip (Fig. 4a). The performance of each infusion round was assessed by monitoring the integrated single-cell fluorescence intensity resulting from binding of the imager probes, and the identity of the species was decoded from a species dictionary after six imaging rounds. Beyond the multiplexed 16S rRNA assay that facilitates bacterial species identification, we also employ a pan-bacterial probe, EUB338-Cy5 (Fig. 5a). This probe serves a dual purpose: it detects the presence of bacteria not specifically targeted in the probe panel, and aids in segmenting and distinguishing bacteria from matrix components.


image file: d4lc00325j-f5.tif
Fig. 5 a) Exemplar sequential hybridization images of an E. coli bacterial strain assayed on the ACBC device. We note that a small number of bacteria detach from the chitosan coating of the device over the course of the sequential hybridization infusions. Scale bar is 15 μm. b) Exemplar single-cell bacteria images and their relative single-cell fluorescence intensities (normalised to brightest probe) after assay with the multiplexed 16S-rRNA imager probes.

To assess the performance of our assay in identifying bacterial species, we employed a classifier which first normalized the single-cell intensity signal of the 6 imager probes based on the intensity of the brightest probe, and identified the two imager probes exhibiting the highest single-cell intensity. All other probe signals were regarded as null (i.e., a “0” entry). If either of the two highest single-cell signals obtained were below a threshold value of 2 standard deviations above non-specific background, then they were also regarded as a negative (i.e., a “0” entry), and the species was subsequently classified as inconclusive. After converting the 6-imager signals to a binary six-entry vector barcode, this barcode was cross-referenced against a dictionary collated from the assigned bacterial barcodes termed here as the “species dictionary”. In addition, if a barcode was raised that was not present in this dictionary, then the cell was designated to the ‘inconclusive’ class suggesting that no identification could be made. Otherwise, one of the 7 species classes was assigned.

To verify the robustness of the assay, the classification performance was assessed across 4 strains of each of the 7 species tested. Our unsupervised approach was able to achieve high classification accuracies at the single-cell level, which we have categorized into three classes. The first class, which we refer to as ‘very high accuracy species’, includes E. coli and Klebsiella pneumoniae, both with a classification rate of 98%, and Enterococcus faecalis with rates of 96%. The second class, termed ‘high accuracy species’, consists of Staphylococcus aureus and Pseudomonas aeruginosa, achieving classification rates of 89% and 87% respectively. The third and final class, the ‘moderate accuracy species’, is represented by Streptococcus agalactiae and Streptococcus pneumoniae, which had a classification rates of 81% and 70% respectively (Fig. 4b). A detailed breakdown of per-strain classification rates can be found in the ESI (Table S3).

This variation in classification accuracy of the assay across the three classes of bacteria can be attributed, at least in part, to the differential susceptibility of bacterial populations, particularly Gram-positive species, to lysozyme treatment. Lysozyme treatment is necessary in order for FISH probes to access the ribosomal content. However, within a given bacterial population, individual bacteria may exhibit varying degrees of susceptibility to lysozyme treatment over the 20-minute treatment period (Fig. S5b). This variability in permeabilization can impact the accessibility of the probes to the ribosomal content of the more lysozyme resistant cells, which in turn can affect the classification accuracy. Specifically, species that are less susceptible to lysozyme treatment, and therefore less permeabilized, are more likely to fall into the ‘moderate’ and ‘high’ accuracy classes, with only a single species achieving a “very high” classification rate.

Permeabilization variability was not only limited to species-specific responses but also manifested on a strain-by-strain basis. For instance, some strains of S. pneumoniae displayed enhanced resilience to lysozyme treatment (SPN3, Table S3). Despite this variability, it is noteworthy that the assay consistently demonstrated a robust capacity to accurately classify the bacterial species across different strains. While higher lysozyme concentrations could enable both faster accessibility and larger fraction of permeabilized ‘high’ and ‘moderate’ accuracy class bacteria, we observed that the remaining bacteria in this panel can lyse at higher lysozyme conditions. As a result, and considering the excellent performance of the classification assay across all species, this set of conditions was chosen for its ability to universally allow classification across all species in the present panel.

When evaluated against the panel of 7 species, the Multiplexed 16S rRNA-FISH assay shows both high sensitivity and specificity across the 7 pathogens, notably: E. coli (sensitivity: 0.98, specificity: 1.00), K. pneumoniae (sensitivity: 0.98, specificity: 1.00), P. aeruginosa (sensitivity: 0.91, specificity: 0.99), E. faecalis (sensitivity: 0.97, specificity: 1.00), S. pneumoniae (sensitivity: 0.79, specificity: 0.98), S. agalactiae (sensitivity: 0.87, specificity: 0.99), and S. aureus (sensitivity: 0.94, specificity: 0.99).

To determine the minimum number of bacterial cells that would need to be assayed before calling the presence of one of the 7 pathogens within a sample, we calculated the true positive (TP), true negative (TN), false positive (FP), and false negative (FN) rates of the multiplexed 16S-rRNA FISH assay per species from the confusion matrix in Fig. 4b and calculated the assay's sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV). The calculation of these parameters is described in the ESI section and the tabulated characteristics of the assay can be found in Table S4.

Using the positive predictive value (PPV) of the assay, we determined the minimal number of positive single-cell identifications that would be needed to assert the presence of a species with 99.5% confidence (Table S4). Then, using the number of positive calls needed, we determined the number of cells that would need to be isolated by the ACBC chip to guarantee an amount of sample that yields sufficient statistical power for a confident identification. Using the number of positive calls needed and based on the TP rate of each species, we estimated the number of cells that would need to be isolated for each species to have a 99.5% confidence that a positive species assignment can be made. For each species these were: n = 2 cells for E. coli and K. pneumoniae, n = 3 cells for E. faecalis, n = 4 cells for S. aureus, n = 5 cells for P. aeruginosa, n = 7 cells for S. agalactiae and n = 9 cells for S. pneumoniae.

Of course, when guiding a diagnosis, a decision would not be made on the basis of single-cell classification alone. This is due to the inherent variability among individual cells within a population, which can lead to discrepancies when extrapolating single-cell data to a broader context. Furthermore, the FP rate of the assay for some of the species could lead to the minimum number of positive calls to be exceeded leading to confounding interpretations. To confidently ascertain the identity of the pathogens present in a sample, the species call counts determined in a sample measurement are placed in an 8-entry vector of observed counts (7 species + inconclusive class) and Fisher's exact test is employed to compare the observed counts against the expected assay performance outcome (Fig. 4b) that would arise from a monomicrobial infection of each of the 7 species in our panel. When the observed counts align with the expected assay performance above a threshold of p > 0.005, the sample is categorized as ‘very likely’ to be containing a monomicrobial infection of the pathogen in question. This approach therefore accounts for the low cross-reactivity observed between species.

An advantage of having single-cell resolution for identification is in the ability to detect the presence of multiple pathogens in an infection. The ability to differentiate between species in mixed infections is particularly crucial in clinical settings, where accurate identification can guide targeted antibiotic treatment strategies. Mixed infections often require a more nuanced approach to treatment compared to mono-bacterial infections, as they may involve pathogens with different antibiotic resistance profiles. To determine the presence of a mixed infection in the sample, if more than one species surpass the threshold for positive calls in the observed counts vector, and simultaneously fails the ‘very likely’ Fisher's exact test outcome against monomicrobial expectations, then a composite presence in the sample is suggested and the sample is classified as a mixed infection.

To demonstrate that we can distinguish mixed infections using our multiplexed 16S rRNA-FISH assay such as in the case of the two Gram-negative rods K. pneumoniae and P. aeruginosa, we stained the former with a membrane stain, wheat germ agglutinin (WGA-AF488), and immobilized the two species to a chitosan-coated surface mixed at a 1[thin space (1/6-em)]:[thin space (1/6-em)]1 cell density ratio. After immobilization and after noting the position and identity of the bacteria, the species were treated with permeabilization solution and hybridized with the encoding probe solution. The identity of the two species was finally verified using the sequential imager probe hybridizations. These two species would have been hard to distinguish by purely morphological features or Gram-stain differentiation (Fig. 6b); the precise identification and differentiation of pathogens within a mixed infection, as enabled by our assay, can inform the selection of more effective, personalized antibiotic regimens.


image file: d4lc00325j-f6.tif
Fig. 6 a) Schematic showcasing a mixed infection scenario. b) A field of view showcasing a mixed infection scenario by the Gram-negative rods K. pneumoniae and P. aeruginosa. True-positive imager probes for the former (I1 and I6) and the latter (I2 and I4) selectively bind to their respective species target suggesting minimal cross-binding of the encoding probes. Negative probes (I3 and I5) show minimal fluorescence signal above background upon excitation. The pan-bacterial kingdom probe, EUB338, is employed at the end of the assay to stain all bacteria present in the field of view. Scale bar is 5 μm.

Ultra-rapid antibiotic susceptibility testing by ACBC enabled imaging flow cytometry

A crucial advantage of microscopy methods for rapid bacterial diagnostics is their ability to preserve the integrity of cells, thus enabling antibiotic susceptibility tests to be conducted at the single-cell level. In previous work by our group, we demonstrated that ultra-rapid antibiotic susceptibility testing (AST) can be performed for E. coli by fluorescence imaging of cellular landmarks in antibiotic-treated cells. These landmarks include the spatial profiles of the nucleoid, membrane, and ribosomes.22,35 This technique offers a significant advantage over more conventional cytological profiling techniques. For instance, rapid AST methods that rely on observing bacterial growth rates under antibiotic conditions can face a bottleneck when the bacteria under investigation have intrinsically slow doubling times.36 By contrast, this technique requires only a 30-minute antibiotic treatment, followed by fixation and staining of the cells with DAPI (for nucleoid staining) and Nile red (for membrane staining).

Here, we demonstrate that the ACBC device is compatible with this ultra-rapid AST method and can function as an imaging flow cytometer for cell classification (Fig. 7). By virtue of the capture region's reduced dimensions, E. coli cells are captured in the same focal plane, thus allowing for continuous imaging of incoming cells and classification of cells as responsive or unresponsive to the antibiotic (ESI: Video S6).


image file: d4lc00325j-f7.tif
Fig. 7 Ultra-rapid AST using the ACBC device. A) MG1655 E. coli are treated with an antibiotic of choice for 30 minutes and are then fixed and stained with Nile red and DAPI. The cells are then captured using the ACBC chip, imaged and classified as susceptible or resistant using a convolutional neural network (CNN). B) Antibiotic response phenotypes of E. coli. i) Untreated/resistant phenotype, ii) ciprofloxacin sensitive phenotype (compacted nucleoid, longer length) and iii) rifampicin sensitive phenotype (expanded nucleoid). Scale bar is 3 μm. C) Capture timelapse of ciprofloxacin-treated MG1655 E. coli cells (0.5 mg L−1, 1 × 105 CFU mL−1) at early timepoints (a)–(c) of an ultra-rapid AST experiment. Cells were stained with Nile red (membrane, red) and DAPI (nucleoid, green). The dashed green line denotes the capture region boundary. Scale bar is 10 μm. The imaging flow cytometry timelapse can be found in the ESI (Video S6). D) and E) Confusion matrices for the phenotype classifier. The treatment conditions (CIP = ciprofloxacin, RIF = rifampicin) are shown on the columns and the model's prediction is shown on the rows.

We performed AST tests on E. coli MG1655 using pre-trained classifiers (classifiers can be found in ref. 22) for 2 antibiotics; ciprofloxacin and rifampicin. In our experiments, both antibiotics were used at concentrations above the minimum inhibitory concentration (MIC) of the MG1655 strain. Ciprofloxacin produces a compaction of the nucleoid towards the cell centre due to topoisomerase IV inhibition (Fig. 7Bii).37 Conversely, rifampicin, produces nucleoid decompaction by inhibiting transcription initiation by RNA polymerase (Fig. 7Biii).22

For ciprofloxacin-treated E. coli cells (concentration used: 0.5 mg L−1, MIC = 0.012 mg L−1), our classifier achieved an accuracy of 89.4% in identifying cells as responsive to the treatment (Ncells = 269) (Fig. 7D). 89.4% of treated cells were correctly classified as ciprofloxacin-sensitive, while 10.6% were misclassified as false-negatives. For untreated cells, 99.88% were correctly identified, with a false-positive rate of 0.12%.

Similarly, for rifampicin-treated E. coli (concentration used: 100 mg L−1, MIC = 8 mg L−1), the classifier demonstrated a higher accuracy rate of 96.6% for detecting response to the treatment (Ncells = 58), (Fig. 7E). As can be seen in confusion matrix, 96.6% of rifampicin-treated cells were accurately classified as responsive, and only 3.45% were misclassified as false-negative. For untreated cells, the classifier maintained an accuracy of 99.12%, with a false-positive rate of 0.88%.

Conclusions

In this study we have demonstrated the development of a microfluidic platform capable of capturing and identifying common human bacterial pathogenic species using adaptive channels capable of forming capture regions for bacteria. By carefully controlling the channel backpressure, we show that sustained bacterial trapping can be achieved. The actuated nature of the system enables the device to act either in capture mode or interrogation mode, circumventing a major issue that common hydrodynamic trapping-based devices have limits infusion of reagents at high flow-rates. In the capture mode, our device can isolate micron-sized objects from solution, albeit at low flow-rates (2.1 μL min−1); in interrogation mode, however, reagents can be flown at much faster rates (up to 50 μL min−1 for FISH reagents, up to 100 μL min−1 for aqueous washes) and therefore enable the rapid exchange of reagents as would be needed with current state of the art spatial transcriptomic techniques.

Future work in the design of the ACBC chip will focus on assessing and increasing clinical sample processing speeds, which will rely on introducing parallel elements in the device design, as well as tuning the hydrophilicity of the capture-region channel that forms upon channel collapse. Further refinements in the architecture of future devices will also be investigated to mitigate the need for surface chemistry modifications required to achieve bacterial-immobilization post-capture. Such approach will enable the microscopic investigation of bacteria otherwise not as prone to immobilization in the poly-cationic surfaces presently used (chitosan) in the ACBC device. We expect these modifications to further facilitate the passage of the aqueous sample matrix, achieve higher overall flow rates under the same channel backpressure conditions and enable universal capture and assay of bacteria found in clinical samples.

Using the ACBC device, we have shown that we can capture bacterial cells from simulated samples in which patient isolates are spiked in known concentrations into buffers. We have also shown that the species of these captured bacteria can be reliably ascertained using a multiplexed 16S-rRNA FISH molecular barcoding method that achieves high classification rates for seven species in an unsupervised manner. In addition to mono-bacterial infections, our assay's capacity to identify and differentiate multiple pathogens within a single sample can provide a more comprehensive understanding of polymicrobial infections. This could potentially reveal interactions between different bacterial species within the infection, offering further insights into the pathogenesis and progression of the disease. Ultimately, these detailed insights could contribute to the development of more effective treatment approaches for patients.

Achieving devices capable of identification applicable for clinical settings requires strong performance in multiple domains such as the amount of sample processed, processing speed, high capture efficiency, high specificity and the ability to perform phenotyping in addition to identification of species. Bacterial identification using multiplex 16S-rRNA FISH in tandem with capture using adaptive channels poses significant advantages in these areas. Compared to the state-of-the-art methods summarized in Table S5, our device showcases a lower limit of detection (7 × 102 CFU mL−1) compared to similar microscopy based methods,13,20 rapid processing time (≈60 minutes), ability to work with low volume samples (42 μL).

Additionally, the ACBC device's ability to capture and identify multiple bacterial species (7) in a single run, combined with its very high capture efficiencies and adaptability for multiplexed fluorescence in situ hybridization (FISH), highlights its potential for integration into clinical diagnostic workflows. These improvements collectively emphasize the novelty and utility of our approach in providing a rapid, sensitive, and scalable solution for bacterial pathogen identification and characterization.

By taking advantage of the ACBC chip's ability to capture bacteria in the same focal-plane, we were able to employ the device as an imaging flow cytometer. Using this feature, we demonstrated the ACBC device's compatibility with ultra-rapid AST imaging methods,22,35 effectively combining the high accuracy in detecting species identity of the pathogen from very dilute samples with the ability to perform antibiotic susceptibility after the identity of the pathogen was established. The ability to perform these tests at very low concentration of samples represents a significant advancement compared to other state-of-the art identification and AST measurement platforms (Table S5). Presently, we have only demonstrated antibiotic response phenotypes for E. coli only. Future work in this direction will aim to develop imaging-based assays to characterize antibiotic response phenotypes across a wider range of pathogens.

While a limited number of bacterial isolates were evaluated, the performance of this combined capture and identification approach yielded high bacterial species classification rates relevant to the typical CFUs of these samples. We aim to further scale the panel of pathogens that can be currently assayed although we do envisage that potential difficulties will eventually arise due to species homology overlaps in the 16S hypervariable regions as more species are included in the assay panel used. However, we remain confident that by adopting more elaborate barcoding schemes and also by introducing error correction methods in the species dictionary employed, the identification of a wider range of bacteria using this method will be realized. In addition, the inclusion of a universal bacterial probe in principle enables the detection of bacteria not present in the panel, and this result in itself could also be relevant in empiric antibiotic prescribing decisions.

Ethics

Ethical approval for the use of clinical samples and isolates processed by the John Radcliffe Hospital microbiology laboratory in the development of diagnostic assays was granted by the UK's Health Research Authority (London – Queen Square Research Ethics Committee [REC reference: 17/LO/1420]).

Data availability

Code used to analyse results and to control fluidic hardware and is available at GitHub at: https://github.com/KapanidisLab/ACBC_Chip. The data that support the findings of this study are available at Zenodo at: https://doi.org/10.5281/zenodo.13287507.

Author contributions

S. C. designed the ACBC device, carried out experiments, analysed data and wrote the manuscript. P. T. developed software tools used for segmentation, analysis and classification of bacteria. C. F. developed protocols for bacterial isolate culture and provided the strains used in this study. A. F. trained ribosome classifier models used for bacterial segmentation. H. E. S. and J. K. performed preliminary experiments and provided reagents. S. C., A. N. K., N. S. and C. N. conceived the study and interpreted results. All authors reviewed and approved the final manuscript.

Conflicts of interest

The work was carried out using a wide-field microscope from Oxford Nanoimaging, a company in which A. N. K. is a co-founder and shareholder. A. N. K. received no payment for this work from Oxford Nanoimaging, and Oxford Nanoimaging was not involved in any aspect of this work.

Acknowledgements

This work was supported by an Oxford Martin School (by the establishment of the Oxford Martin School Programme on Antimicrobial Resistance Testing; to ANK, NS, CN), by Wellcome Trust grant 110164/Z/15/Z (to ANK) and by UK Biotechnology and Biological Sciences Research Council grants BB/N018656/1 and BB/S008896/1 (to ANK) and by Oxford's EPSRC/BBSRC Impact Acceleration Account EP/X525777/1 (award 0012379, to ANK, NS, and CN). The research was additionally supported by the National Institute for Health Research (NIHR) Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance (NIHR200915) at the University of Oxford in partnership with United Kingdom Health Security Agency (UKHSA) and by the NIHR Oxford Biomedical Research Centre. This work was also supported by a John Fell Fund grant to NS (Grant application: 0008776). NS is an NIHR Oxford BRC Senior Research Fellow. Figures presented contain elements created with https://Biorender.com.

References

  1. C. J. Murray, et al. , Lancet, 2022, 399, 629–655 CrossRef CAS PubMed .
  2. P. Puerta-Alcalde, C. Cardozo, M. Suárez-Lledó, O. Rodríguez-Núñez, L. Morata, C. Fehér, F. Marco, A. Del Río, J. A. Martínez, J. Mensa, M. Rovira, J. Esteve, A. Soriano and C. Garcia-Vidal, Clin. Microbiol. Infect., 2019, 25, 447–453 CrossRef CAS PubMed .
  3. Y. W. Tsai, T. C. Lin, H. Y. Chou, H. Y. Hung, C. K. Tan, L. C. Wu, I. J. Feng and Y. L. Shiue, Diagnostics, 2021, 11(8), 1514 CrossRef CAS PubMed .
  4. L. Ponderand, P. Pavese, D. Maubon, E. Giraudon, T. Girard, C. Landelle, M. Maurin and Y. Caspar, Ann. Clin. Microbiol. Antimicrob., 2020, 19, 60 CrossRef CAS PubMed .
  5. T. H. Huang, Y. L. Tzeng and R. M. Dickson, Cytometry, Part A, 2018, 93, 639–648 CrossRef CAS PubMed .
  6. T. K. Price, T. Dune, E. E. Hilt, K. J. Thomas-White, S. Kliethermes, C. Brincat, L. Brubaker, A. J. Wolfe, E. R. Mueller and P. C. Schreckenberger, J. Clin. Microbiol., 2016, 54, 1216–1222 CrossRef CAS PubMed .
  7. J. Mai, V. V. Abhyankar, M. E. Piccini, J. P. Olano, R. Willson and A. V. Hatch, Biosens. Bioelectron., 2014, 54, 435–441 CrossRef CAS PubMed .
  8. G. Pitruzzello, S. Thorpe, S. Johnson, A. Evans, H. Gadêlha and T. F. Krauss, Lab Chip, 2019, 19, 1417–1426 RSC .
  9. Ö. Baltekin, A. Boucharin, E. Tano, D. I. Andersson and J. Elf, Proc. Natl. Acad. Sci. U. S. A., 2017, 114, 9170–9175 CrossRef PubMed .
  10. S. H. Jung, Y. K. Hahn, S. Oh, S. Kwon, E. Um, S. Choi and J. H. Kang, Small, 2018, 14(34), 1801731 CrossRef PubMed .
  11. M. R. Vutukuru, D. K. Sharma, M. S. Ragavendar, S. Schmolke, Y. Huang, W. Gumbrecht and N. Mitra, J. Microbiol. Methods, 2016, 131, 105–109 CrossRef CAS PubMed .
  12. P. Liu, H. Liu, L. Semenec, D. Yuan, S. Yan, A. K. Cain and M. Li, Microsyst. Nanoeng., 2022, 8, 7 CrossRef CAS PubMed .
  13. S. Stratz, K. Eyer, F. Kurth and P. S. Dittrich, Anal. Chem., 2014, 86, 12375–12381 CrossRef CAS PubMed .
  14. B. Okumus, C. J. Baker, J. C. Arias-Castro, G. C. Lai, E. Leoncini, S. Bakshi, S. Luro, D. Landgraf and J. Paulsson, Nat. Protoc., 2018, 13, 170–194 CrossRef CAS PubMed .
  15. J. D. Baker, D. T. Kysela, J. Zhou, S. M. Madren, A. S. Wilkens, Y. V. Brun and S. C. Jacobson, Anal. Chem., 2016, 88, 8476–8483 CrossRef CAS PubMed .
  16. H. Li, P. Torab, K. E. Mach, C. Surrette, M. R. England, D. W. Craft, N. J. Thomas, J. C. Liao, C. Puleo and P. K. Wong, Proc. Natl. Acad. Sci. U. S. A., 2019, 116, 10270–10279 CrossRef CAS PubMed .
  17. A. Borók, K. Laboda and A. Bonyár, Biosensors, 2021, 11, 292 CrossRef PubMed .
  18. K. H. Chen, A. N. Boettiger, J. R. Moffitt, S. Wang and X. Zhuang, Science, 2015, 348(6233), 412 CrossRef CAS PubMed .
  19. H. Shi, Q. Shi, B. Grodner, J. S. Lenz, W. R. Zipfel, I. L. Brito and I. De Vlaminck, Nature, 2020, 588, 676–681 CrossRef CAS PubMed .
  20. V. Kandavalli, P. Karempudi, J. Larsson and J. Elf, Nat. Commun., 2022, 13, 6215 CrossRef CAS PubMed .
  21. C. Stringer, T. Wang, M. Michaelos and M. Pachitariu, Nat. Methods, 2021, 18, 100–106 CrossRef CAS PubMed .
  22. A. Zagajewski, P. Turner, C. Feehily, H. El Sayyed, M. Andersson, L. Barrett, S. Oakley, M. Stracy, D. Crook, C. Nellåker, N. Stoesser and A. N. Kapanidis, Commun. Biol., 2023, 6, 1164 CrossRef CAS PubMed .
  23. T. Thorsen, J. S. Maerkl and R. S. Quake, Science, 2002, 298, 580–584 CrossRef CAS PubMed .
  24. T. Hansson and S. Sjölander, Valve especially for fluid handling bodies with microchannels, US Pat., US005593130A, 1997 Search PubMed .
  25. J. Tréguier, L. Bugnicourt, G. Gay, M. Diallo, S. T. Islam, A. Toro, L. David, O. Théodoly, G. Sudre and T. Mignot, mBio, 2019, e01375-19 CrossRef PubMed .
  26. L. Hunter, J. Gala de Pablo, A. C. Stammers, N. H. Thomson, S. D. Evans and J. Shim, SN Appl. Sci., 2020, 2, 1501 CrossRef CAS .
  27. N. Compera, S. Atwell, J. Wirth, B. Wolfrum and M. Meier, Lab Chip, 2021, 21, 2986–2996 RSC .
  28. O. Akogwu, D. Kwabi, S. Midturi, M. Eleruja, B. Babatope and W. O. Soboyejo, Mater. Sci. Eng., B, 2010, 170, 32–40 CrossRef CAS .
  29. D. Valesco Anez, C. Hadji, S. Carreras, É. Lorenceau and C. Picard, Microfluid. Nanofluid., 2021, 7 DOI:10.1007/s10404-020-02408-5 .
  30. I. Mihalcescu, M. Van-Melle Gateau, B. Chelli, C. Pinel and J.-L. Ravanat, Phys. Biol., 2015, 12, 10 Search PubMed .
  31. L. Ceriotti, N. F. De Rooij and E. Verpoorte, Anal. Chem., 2002, 74, 639–647 CrossRef CAS PubMed .
  32. L. G. Harris, S. J. Foster, R. G. Richards, P. Lambert, D. Stickler and A. Eley, Eur. Cells Mater., 2002, 4, 39–60 CrossRef CAS PubMed .
  33. E. W. Sayers, E. E. Bolton, J. R. Brister, K. Canese, J. Chan, D. C. Comeau, R. Connor, K. Funk, C. Kelly, S. Kim, T. Madej, A. Marchler-Bauer, C. Lanczycki, S. Lathrop, Z. Lu, F. Thibaud-Nissen, T. Murphy, L. Phan, Y. Skripchenko, T. Tse, J. Wang, R. Williams, B. W. Trawick, K. D. Pruitt and S. T. Sherry, Nucleic Acids Res., 2022, 50, D20–D26 CrossRef CAS PubMed .
  34. Q. Xu, M. R. Schlabach, G. J. Hannon, S. J. Elledge and B. Watson, Proc. Natl. Acad. Sci. U. S. A., 2009, 106, 2289 CrossRef CAS PubMed .
  35. A. Farrar, P. Turner, H. El Sayyed, C. Feehily, S. Chatzimichail, D. Crook, M. Andersson, S. Oakley, L. Barrett, C. Nellåker, N. Stoesser and A. Kapanidis, medRxiv, 2024, preprint,  DOI:10.1101/2024.06.18.24309111.
  36. S. Srivastava, S. P. Van Rijn, A. M. A. Wessels, J. W. C. Alffenaar and T. Gumbo, Antimicrob. Agents Chemother., 2016, 60, 3193–3195 CrossRef CAS PubMed .
  37. E. Vickridge, C. Planchenault, C. Cockram, I. G. Junceda and O. Espéli, Nat. Commun., 2017, 8, 14618 CrossRef PubMed .

Footnote

Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4lc00325j

This journal is © The Royal Society of Chemistry 2024
Click here to see how this site uses Cookies. View our privacy policy here.