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
First published on 6th September 2024
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.
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.
For the flow-layer, PDMS (Sylgard 184, Dow Corning, United States) mixed at a ratio of 20:
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
:
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
:
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.
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.
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).
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.
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).
![]() | ||
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 − (PQ − PC) | (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.
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).
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†).
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.
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:
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.
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).
![]() | ||
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%.
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.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4lc00325j |
This journal is © The Royal Society of Chemistry 2024 |