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

HAL
open science

# Microfluidic dose-response platform to track the dynamics of drug response in single mycobacterial cells

Maxime Mistretta, Nicolas Gangneux, Giulia Manina

► To cite this version:

Maxime Mistretta, Nicolas Gangneux, Giulia Manina. Microfluidic dose-response platform to track the dynamics of drug response in single mycobacterial cells. Scientific Reports, 2022, 12 (1), pp.19578. 10.1038/s41598-022-24175-9. hal-04613629

HAL Id: hal-04613629

https://hal.science/hal-04613629

Submitted on 16 Jun 2024

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# Microfluidic dose-response platform to track the dynamics of drug response in single mycobacterial cells

Maxime Mistretta, Nicolas Gangneux, Giulia Manina

► To cite this version:

Maxime Mistretta, Nicolas Gangneux, Giulia Manina. Microfluidic dose-response platform to track the dynamics of drug response in single mycobacterial cells. Scientific Reports, 2022, 12 (1), pp.19578. 10.1038/s41598-022-24175-9. hal-04613629

HAL Id: hal-04613629

https://hal.science/hal-04613629v1

Submitted on 16 Jun 2024

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

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

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

OPEN

# Microfluidic dose-response platform to track the dynamics of drug response in single mycobacterial cells

Maxime Mistretta, Nicolas Gangneux & Giulia Manina*

Preclinical analysis of drug efficacy is critical for drug development. However, conventional bulk-cell assays statically assess the mean population behavior, lacking resolution on drug-escaping cells. Inaccurate estimation of efficacy can lead to overestimation of compounds, whose efficacy will not be confirmed in the clinic, or lead to rejection of valuable candidates. Time-lapse microfluidic microscopy is a powerful approach to characterize drugs at high spatiotemporal resolution, but hard to apply on a large scale. Here we report the development of a microfluidic platform based on a pneumatic operating principle, which is scalable and compatible with long-term live-cell imaging and with simultaneous analysis of different drug concentrations. We tested the platform with mycobacterial cells, including the tubercular pathogen, providing the first proof of concept of a single-cell dose-response assay. This dynamic in-vitro model will prove useful to probe the fate of drug-stressed cells, providing improved predictions of drug efficacy in the clinic.

Persistence of bacterial pathogens results in long-term host colonization, which can be silent or symptomatic¹,². Persistent infections are also typically recalcitrant to treatment, due to diverse bacterial drug-escape mechanisms, with consequences for morbidity and mortality³⁻⁵. Not surprisingly, pathogens responsible for persistent infections have been often implicated in drug resistance, and populate the priority list of antibiotic-resistant pathogens⁶. Drug resistance results from spontaneous mutations, which are primarily associated with drug-target modification, drug inactivation or decreased membrane permeability, causing an increase in the minimum inhibitory concentration (MIC); conversely, drug persistence is a transient phenotypic change that is not genetically inherited, and does not cause MIC change at the population level⁴. Drug persistence is insidious and can result from several host and pathogen factors, such as: drug absorption, distribution, metabolism and elimination⁷⁻⁹; drug-refractory barriers built by bacterial populations, growing in the form of cellular aggregates or organized communities³,¹⁰; and phenotypic variation, an inherent feature of any clonal cell population that increases during adverse conditions. This phenotypic diversification is likely to enhance the adaptive potential of the bacterial population, contributing to drug-escape mechanisms and ultimately altering the efficacy of drugs¹¹⁻¹⁷. Drug persisters can be thought of as phenotypic variants within a clonal population that endure drug cidality, despite being genetically indistinguishable from drug-susceptible bacteria. Drug persistence is commonly regarded as a specific subset of drug tolerance, leading to a slowdown in the rate of drug-dependent killing, with no change in the MIC, often but not always associated with decreased growth rate⁴,¹³,¹⁸,¹⁹. A few studies have also begun to demonstrate a relationship between drug tolerance and resistance²⁰⁻²⁵, implying that hampering drug tolerance mechanisms might prevent the onset of drug resistance.

*Mycobacterium tuberculosis* is the epitome of a persistent pathogen, and tuberculosis treatment is exceptionally protracted, due to both pharmacokinetic (PK) and pharmacodynamic (PD) limitations⁹,²⁶. In particular, the tubercular pathogen occupies host niches that are refractory to drugs, where it exhibits marked growth and metabolic variation, which is associated with drug persistence⁷,⁸,¹¹,¹⁵. Although non-growing bacilli are enriched with drug persisters¹⁻⁵, phenotypic drug-tolerance is a more complex phenomenon, contingent on variation in different cellular processes, and reported also in replicating cells³,¹¹,¹³,²⁷. Understanding of this phenomenon is still incomplete, mostly due to indirect investigative approaches, and this undermines preclinical drug discovery and optimization. Conventional dose-response relationship studies of drug efficacy depend on the readout used to determine the impact of drugs on pathogen fitness²⁶,²⁸. These approaches primarily rely on the determination

Microbial Individuality and Infection Laboratory, Institut Pasteur, Université de Paris, 75015 Paris, France. * email: giulia.manina@pasteur.fr
<figure><img src="image_2.png" /></figure>

OPEN

# Microfluidic dose-response platform to track the dynamics of drug response in single mycobacterial cells

Maxime Mistretta, Nicolas Gangneux & Giulia Manina*

Preclinical analysis of drug efficacy is critical for drug development. However, conventional bulk-cell assays statically assess the mean population behavior, lacking resolution on drug-escaping cells. Inaccurate estimation of efficacy can lead to overestimation of compounds, whose efficacy will not be confirmed in the clinic, or lead to rejection of valuable candidates. Time-lapse microfluidic microscopy is a powerful approach to characterize drugs at high spatiotemporal resolution, but hard to apply on a large scale. Here we report the development of a microfluidic platform based on a pneumatic operating principle, which is scalable and compatible with long-term live-cell imaging and with simultaneous analysis of different drug concentrations. We tested the platform with mycobacterial cells, including the tubercular pathogen, providing the first proof of concept of a single-cell dose-response assay. This dynamic in-vitro model will prove useful to probe the fate of drug-stressed cells, providing improved predictions of drug efficacy in the clinic.

Persistence of bacterial pathogens results in long-term host colonization, which can be silent or symptomatic¹,². Persistent infections are also typically recalcitrant to treatment, due to diverse bacterial drug-escape mechanisms, with consequences for morbidity and mortality³⁻⁵. Not surprisingly, pathogens responsible for persistent infections have been often implicated in drug resistance, and populate the priority list of antibiotic-resistant pathogens⁶. Drug resistance results from spontaneous mutations, which are primarily associated with drug-target modification, drug inactivation or decreased membrane permeability, causing an increase in the minimum inhibitory concentration (MIC); conversely, drug persistence is a transient phenotypic change that is not genetically inherited, and does not cause MIC change at the population level⁴. Drug persistence is insidious and can result from several host and pathogen factors, such as: drug absorption, distribution, metabolism and elimination⁷⁻⁹; drug-refractory barriers built by bacterial populations, growing in the form of cellular aggregates or organized communities³,¹⁰; and phenotypic variation, an inherent feature of any clonal cell population that increases during adverse conditions. This phenotypic diversification is likely to enhance the adaptive potential of the bacterial population, contributing to drug-escape mechanisms and ultimately altering the efficacy of drugs¹¹⁻¹⁷. Drug persisters can be thought of as phenotypic variants within a clonal population that endure drug cidality, despite being genetically indistinguishable from drug-susceptible bacteria. Drug persistence is commonly regarded as a specific subset of drug tolerance, leading to a slowdown in the rate of drug-dependent killing, with no change in the MIC, often but not always associated with decreased growth rate⁴,¹³,¹⁸,¹⁹. A few studies have also begun to demonstrate a relationship between drug tolerance and resistance²⁰⁻²⁵, implying that hampering drug tolerance mechanisms might prevent the onset of drug resistance.

*Mycobacterium tuberculosis* is the epitome of a persistent pathogen, and tuberculosis treatment is exceptionally protracted, due to both pharmacokinetic (PK) and pharmacodynamic (PD) limitations⁹,²⁶. In particular, the tubercular pathogen occupies host niches that are refractory to drugs, where it exhibits marked growth and metabolic variation, which is associated with drug persistence⁷,⁸,¹¹,¹⁵. Although non-growing bacilli are enriched with drug persisters¹⁻⁵, phenotypic drug-tolerance is a more complex phenomenon, contingent on variation in different cellular processes, and reported also in replicating cells³,¹¹,¹³,²⁷. Understanding of this phenomenon is still incomplete, mostly due to indirect investigative approaches, and this undermines preclinical drug discovery and optimization. Conventional dose-response relationship studies of drug efficacy depend on the readout used to determine the impact of drugs on pathogen fitness²⁶,²⁸. These approaches primarily rely on the determination

Microbial Individuality and Infection Laboratory, Institut Pasteur, Université de Paris, 75015 Paris, France. *email: giulia.manina@pasteur.fr
of colony-forming units and of other viability parameters under static environmental conditions, inherently over-
looking minor subpopulations involved in persistence<sup>29,30</sup>. In an effort to capture drug persistence, a few bulk-cell
assays have been introduced<sup>31</sup>. For instance, drug-dependent killing kinetics that show a biphasic pattern with
a reduction in the minimum duration of killing, and growth resumption of colonies following disk-diffusion
test are both suggestive of persistence. Although valuable, these approaches are limited and retrospective, espe-
cially for slow-growing pathogens such as M. tuberculosis, where the long time required for colony formation
is likely to bring inferential and inaccurate results. As a consequence, preclinical in-vitro characterization of
anti-tubercular drugs has seldom been a good predictor of sterilizing activity in humans, especially against drug-
persistent subpopulations<sup>26,32</sup>. Given the complexity of bacterial drug-evasion mechanisms and their implication
in therapeutic failures, preclinical study of drug efficacy requires the integration of cutting-edge approaches,
which examine the PK-PD properties at higher spatial and temporal resolution, providing direct insights into
the single-cell responses to drugs. Improved understanding of the dynamics of formation of drug-recalcitrant
subpopulations would support decision making on the ranking of clinically relevant drugs, and provide more
robust information to guide the transition from preclinical evaluation to clinical implementation<sup>33</sup>.

Over the past two decades, nano- and micro-technologies, such as lithography, molding and device packaging, have proven their utility in biotechnology and biomedical fields³⁴,³⁵. The development of microfluidics-based microsystems in conjunction with time-resolved fluorescence microscopy of live cells has made it possible to quantitatively access the complexity and dynamics of microbial populations under tight environmental control with unprecedented accuracy³⁶⁻⁴⁵. The transition from bulk-scale analysis of non-adherent cell suspensions to two-dimensional tracking of individual cells over time at the micrometric scale required enclosure systems compatible with cell size, growth, viability, changing environment and imaging. To this aim, several microfluidic systems based on different operating principles have been introduced, typically made of micropatterned polydimethylsiloxane (PDMS) and glass⁴⁶,⁴⁷. The most popular examples are the mother machine and its derivative microfluidic chemostats³⁸,⁴⁸. They consist of a series of narrow PDMS microchannels, whereby single progenitor cells give rise to successive generations that can be monitored until they flow into perpendicular channels, which function as both nutrient source and exit. Other PDMS-glass devices consist of growth microchannels or microchambers, of diverse shape and size, which are interconnected in fluidic networks⁴²,⁴⁵,⁴⁹,⁵⁰. While these systems are relatively functional for large adherent cells, bacterial cells find room to move, leading to three-dimensional growth and making cell tracking unstable over time. To overcome these limitations, other systems were developed to force monolayer growth, either by embedding bacteria within hydrocolloid layers, or by applying gentle mechanical compression, and providing nutrients from PDMS microchannels, which are separated from the bacterial monolayer by semipermeable membranes¹³,³⁶,³⁷. Although these systems provide superior resolution and more stable two-dimensional growth over time, the growth areas cannot be easily multiplied nor can they be sealed independently, making these systems useless for multi-condition applications. Another way to confine cells is based on the use of pneumatic valves, originally created to regulate fluids, by opening and closing multilayer networks of microchannels⁵¹. Although easy to control, valve-based microsystems operate at very high pressures, which are appropriate for high-resolution imaging of intracellular compartments, but are not compatible with long-term cell viability³⁹. On the other hand, valve arrays controlled by pneumatic channels were also used to create large dome-shaped microchambers suitable either for imaging large adherent cells or for enzymatic assays on single bacterial lysates⁵²,⁵³. A separate category of microfluidic devices, compatible not only with cell survival and molecular analysis but also with high throughput assays, is the one based on hydrolipidic droplets⁵⁴,⁵⁵. Droplet-based microfluidics are being used to trap different cell types, which grow in the form of aggregates. This growth configuration impairs single-cell resolution, and requires the acquisition of several image stacks, causing phototoxicity and challenging image-analysis workflows. In summary, although several microsystems have been introduced to probe single cells, there is yet no system that stably traps very small microorganisms under separate environmental conditions, preserving their viability during prolonged periods of observation by microscopy.

Here, we report the development of a prototype microfluidic system consisting of multiple micropatterned PDMS layers bonded to glass, where an array of independent culture chambers is combined with a classical dilution tree⁵⁶. The peculiarity of this platform lies in its fundamental unit, which is a microchamber that traps cells through an unconventional hydro-pneumatic mechanism, allowing them to grow as a monolayer, receiving nutrients by active diffusion from the periphery of the cell-growth area. Our system is compatible with time-resolved imaging of single bacterial cells and microcolonies on a single focal plane, thus reducing phototoxicity. Furthermore, thanks to the multiplication of its fundamental unit and the integration of a double-input dilution tree, we achieve the formation of a concentration gradient, which simultaneously reaches separate chambers. The whole system is driven by a multi-channel pressure controller, allowing constant adjustment of either pressure or flow in the functional layers of the microfluidic device. With this platform, we carried out long-term imaging of both fast- and slow-growing mycobacterial cells, notably establishing a procedure compatible with biosafety level 3 (BSL3) containment. In this work, we also demonstrate the use of our system for dynamic dose-response relationship studies, by simulating PK parameters in vitro. As a proof of concept, we analyzed the relationship between concentration and effect of a relevant antitubercular drug over time, at the subpopulation and single-cell level. We found that moxifloxacin impacts mycobacterial growth in a dose- and time-dependent manner, and that clonal cells exhibit heterogeneous responses at the same drug concentration. Furthermore, while inhibition of the target occurred at higher drug concentrations, for concentrations close to the MIC the target was upregulated. We speculate that this may trigger adaptive mechanisms, leading to persistence in a fraction of the population.

We expect our microsystem to ultimately provide superior single-cell PK-PD data, which will better correlate with in-vivo efficacy. This will allow for improved predictions of drug activity against persistent subpopulations, which are challenging to eradicate, thus enhancing the control of infections.
of colony-forming units and of other viability parameters under static environmental conditions, inherently over-
looking minor subpopulations involved in persistence<sup>29,30</sup>. In an effort to capture drug persistence, a few bulk-cell
assays have been introduced<sup>31</sup>. For instance, drug-dependent killing kinetics that show a biphasic pattern with
a reduction in the minimum duration of killing, and growth resumption of colonies following disk-diffusion
test are both suggestive of persistence. Although valuable, these approaches are limited and retrospective, espe-
cially for slow-growing pathogens such as M. tuberculosis, where the long time required for colony formation
is likely to bring inferential and inaccurate results. As a consequence, preclinical in-vitro characterization of
anti-tubercular drugs has seldom been a good predictor of sterilizing activity in humans, especially against drug-
persistent subpopulations<sup>26,32</sup>. Given the complexity of bacterial drug-evasion mechanisms and their implication
in therapeutic failures, preclinical study of drug efficacy requires the integration of cutting-edge approaches,
which examine the PK-PD properties at higher spatial and temporal resolution, providing direct insights into
the single-cell responses to drugs. Improved understanding of the dynamics of formation of drug-recalcitrant
subpopulations would support decision making on the ranking of clinically relevant drugs, and provide more
robust information to guide the transition from preclinical evaluation to clinical implementation<sup>33</sup>.

Over the past two decades, nano- and micro-technologies, such as lithography, molding and device packaging, have proven their utility in biotechnology and biomedical fields³⁴,³⁵. The development of microfluidics-based microsystems in conjunction with time-resolved fluorescence microscopy of live cells has made it possible to quantitatively access the complexity and dynamics of microbial populations under tight environmental control with unprecedented accuracy³⁶⁻⁴⁵. The transition from bulk-scale analysis of non-adherent cell suspensions to two-dimensional tracking of individual cells over time at the micrometric scale required enclosure systems compatible with cell size, growth, viability, changing environment and imaging. To this aim, several microfluidic systems based on different operating principles have been introduced, typically made of micropatterned polydimethylsiloxane (PDMS) and glass⁴⁶,⁴⁷. The most popular examples are the mother machine and its derivative microfluidic chemostats³⁸,⁴⁸. They consist of a series of narrow PDMS microchannels, whereby single progenitor cells give rise to successive generations that can be monitored until they flow into perpendicular channels, which function as both nutrient source and exit. Other PDMS-glass devices consist of growth microchannels or microchambers, of diverse shape and size, which are interconnected in fluidic networks⁴²,⁴⁵,⁴⁹,⁵⁰. While these systems are relatively functional for large adherent cells, bacterial cells find room to move, leading to three-dimensional growth and making cell tracking unstable over time. To overcome these limitations, other systems were developed to force monolayer growth, either by embedding bacteria within hydrocolloid layers, or by applying gentle mechanical compression, and providing nutrients from PDMS microchannels, which are separated from the bacterial monolayer by semipermeable membranes¹³,³⁶,³⁷. Although these systems provide superior resolution and more stable two-dimensional growth over time, the growth areas cannot be easily multiplied nor can they be sealed independently, making these systems useless for multi-condition applications. Another way to confine cells is based on the use of pneumatic valves, originally created to regulate fluids, by opening and closing multilayer networks of microchannels⁵¹. Although easy to control, valve-based microsystems operate at very high pressures, which are appropriate for high-resolution imaging of intracellular compartments, but are not compatible with long-term cell viability³⁹. On the other hand, valve arrays controlled by pneumatic channels were also used to create large dome-shaped microchambers suitable either for imaging large adherent cells or for enzymatic assays on single bacterial lysates⁵²,⁵³. A separate category of microfluidic devices, compatible not only with cell survival and molecular analysis but also with high throughput assays, is the one based on hydrolipidic droplets⁵⁴,⁵⁵. Droplet-based microfluidics are being used to trap different cell types, which grow in the form of aggregates. This growth configuration impairs single-cell resolution, and requires the acquisition of several image stacks, causing phototoxicity and challenging image-analysis workflows. In summary, although several microsystems have been introduced to probe single cells, there is yet no system that stably traps very small microorganisms under separate environmental conditions, preserving their viability during prolonged periods of observation by microscopy.

Here, we report the development of a prototype microfluidic system consisting of multiple micropatterned PDMS layers bonded to glass, where an array of independent culture chambers is combined with a classical dilution tree⁵⁶. The peculiarity of this platform lies in its fundamental unit, which is a microchamber that traps cells through an unconventional hydro-pneumatic mechanism, allowing them to grow as a monolayer, receiving nutrients by active diffusion from the periphery of the cell-growth area. Our system is compatible with time-resolved imaging of single bacterial cells and microcolonies on a single focal plane, thus reducing phototoxicity. Furthermore, thanks to the multiplication of its fundamental unit and the integration of a double-input dilution tree, we achieve the formation of a concentration gradient, which simultaneously reaches separate chambers. The whole system is driven by a multi-channel pressure controller, allowing constant adjustment of either pressure or flow in the functional layers of the microfluidic device. With this platform, we carried out long-term imaging of both fast- and slow-growing mycobacterial cells, notably establishing a procedure compatible with biosafety level 3 (BSL3) containment. In this work, we also demonstrate the use of our system for dynamic dose-response relationship studies, by simulating PK parameters in vitro. As a proof of concept, we analyzed the relationship between concentration and effect of a relevant antitubercular drug over time, at the subpopulation and single-cell level. We found that moxifloxacin impacts mycobacterial growth in a dose- and time-dependent manner, and that clonal cells exhibit heterogeneous responses at the same drug concentration. Furthermore, while inhibition of the target occurred at higher drug concentrations, for concentrations close to the MIC the target was upregulated. We speculate that this may trigger adaptive mechanisms, leading to persistence in a fraction of the population.

We expect our microsystem to ultimately provide superior single-cell PK-PD data, which will better correlate
with in-vivo efficacy. This will allow for improved predictions of drug activity against persistent subpopulations,
which are challenging to eradicate, thus enhancing the control of infections.
## Results

**Construction of a microfluidic culture chamber for time-lapse imaging of bacterial cells.** We conceived an unconventional microfluidic cell-culture chamber (Fig. 1a), inspired by micromechanical push-down valves⁵¹, but surpassing their standards of fabrication. In particular, we designed a hydro-pneumatic microfluidic valve composed of two main PDMS layers, covalently bonded to a glass coverslip. Both the upper and the lower layer, referred to as control layer (CL) and flow layer (FL), respectively, are micropatterned with circular chambers of 1-mm diameter that are superimposed. On each layer the chambers are crossed by 300-µm wide microchannels, which are perpendicular to each other (Fig. 1a). The total height of the CL is 8 mm and the height of its micropatterned structures is 200 µm. The total height of the FL is 50 µm and the height of its micropatterned structures is about 30 µm, with a flexible PDMS top layer of about 20 µm, also referred to as membrane, whose controlled movement constitutes the operating principle of our culture chamber (Fig. 1b).

The top and bottom layers serve different but complementary functions. On the one hand, the CL is filled with water and closed at one end. By using a flow controller to fine-tune the pressure inside the CL, the PDMS membrane is lowered by pneumatic inflation until it contacts the underlying glass coverslip. Lowering of the membrane creates a contact zone between PDMS and glass, suitable for microbial cell containment and two-dimensional growth (Fig. 1b). Using water instead of air to actuate the CL improves the imaging quality, since the refractive indexes of water and PDMS are similar. On the other hand, the FL is continuously perfused with growth medium, providing nutrients to the cells trapped between the PDMS membrane and the glass coverslip. Importantly, the CL is composed of a standard (1:5) ratio between PDMS crosslinker and pre-polymer, as opposed to the FL, composed of a (1:20) ratio, in order to bond the two layers by migration of the crosslinking agent⁵⁷. In addition, we found that treating the FL with oxygen plasma improves the penetration of liquid culture medium (see also "Methods"), due to increased hydrophilicity⁵⁸.

We fabricated a first microsystem prototype, by arraying eight culture chambers in a row, fed by a common flow and controlled by a single CL, laid out perpendicular to the FL (Supplementary Fig. 1a). To prevent accumulation of debris in the FL and clogging of the system, we also integrated PDMS pillars of 100- and 20-µm diameter at the inlet port of the FL. To test the membrane motion, we perfused the FL with a solution of fluorescein isothiocyanate (FITC) at a constant rate of 150 µL/h, measured at the outlet, bringing the membrane to a convex position relative to the coverslip. Next, we applied a gradient of pressures to the CL, and acquired sequential images of the chamber at increasing pressure steps of one mbar per image (Fig. 1c). By measuring the fluorescence profile along the membrane diameter, we found that movement of the membrane from convex, to flat, to concave, is associated with a decrease in fluorescence (Fig. 1d; Supplementary Video 1). Starting from the point of maximum convexity of the membrane, an increment of 5 mbar was sufficient to flatten the membrane, and a further increment of 20 mbar generated a contact zone of approximately 250 µm in diameter (Fig. 1d; Supplementary Fig. 1b). This contact area is sufficiently large to visualize the growth of several bacterial microcolonies, and leaves sufficient room for lateral flow of growth medium. The unusual aspect ratio between diameter and height of the lower microchamber (33:1), over three fold above the manufacturing standards typically used to prevent collapse of microstructures⁵¹, enabled us to use actuation pressures in the CL more than one order of magnitude lower than those used in conventional micromechanical valves⁵⁹. Our ability to work at low pressure in the CL was likely to be beneficial to cellular viability.

Next, we used the FL to load a suspension of *Mycobacterium smegmatis* cells, a non-pathogenic tuberculosis model that duplicates in vitro approximately every 3 h. To identify the best culture conditions, we tested different combinations of pressure and flow rate in the CL and in the FL, and measured the doubling time of microcolonies within the two-dimensional growth area (Fig. 1e). We found that several combinations are compatible with bacterial survival and normal doubling time. However, a drastic decrease of flow rate in the FL impaired bacterial growth. This was likely due to decreased nutrients influx and possibly to excessive pressure exerted on the cells from the CL. Based on these results, we were able to choose the operating conditions of the device that are most compatible with mycobacterial growth, namely, high flow rate in the FL (150 µL/h) and moderate pressure in the CL (30 mbar). Under these conditions, we not only measured similar generation times between separate microfluidic culture chambers (Fig. 1f), implying good reproducibility across the system, but also found that mycobacterial growth was comparable to that obtained with our former single-condition microfluidic device, which is based on a different working principle and operates in the absence of shear stress¹³,⁴⁴. In addition, we calculated that the growth rate of microcolonies was independent of their location in the chamber (Fig. 1g; Supplementary Fig. 1b), suggesting homogeneous growth throughout the two-dimensional growth area. Overall, these data imply that bacteria do not experience apparent stress inside the microchambers, and that shear stress is also likely to be negligible in our system.

In conclusion, we constructed a modular microfluidic culture chamber based on a hydro-pneumatic mechanism. We also proved that this device can trap live bacterial cells, and can be used to monitor their two-dimensional growth.

**Evolution of the microfluidic culture chamber into a multi-condition platform.** To increase the versatility of our microfluidic culture chamber, we arrayed 20 of them, generating five separate groups of four microchambers (Fig. 2a, Supplementary Data 1), and we connected each group to the outlet branch of a classic microfluidic dilution tree⁵⁶,⁶⁰,⁶¹. In this configuration, the CL consists of a straight microchannel that widens in correspondence of each microchamber, ensuring the same pneumatic mechanism of the single functional unit (Fig. 1). The resulting five-condition platform is supervised by a computer-driven multichannel flow-controller. The first channel regulates a tripartite tubing, by applying pressure into three separate bottles, containing either culture medium alone or the highest concentration of the molecule to be tested dissolved in culture medium (Fig. 2a; Supplementary Fig. 2a). In addition, a bidirectional switch valve is placed downstream of two bottles to
## Results

**Construction of a microfluidic culture chamber for time-lapse imaging of bacterial cells.** We conceived an unconventional microfluidic cell-culture chamber (Fig. 1a), inspired by micromechanical push-down valves⁵¹, but surpassing their standards of fabrication. In particular, we designed a hydro-pneumatic microfluidic valve composed of two main PDMS layers, covalently bonded to a glass coverslip. Both the upper and the lower layer, referred to as control layer (CL) and flow layer (FL), respectively, are micropatterned with circular chambers of 1-mm diameter that are superimposed. On each layer the chambers are crossed by 300-µm wide microchannels, which are perpendicular to each other (Fig. 1a). The total height of the CL is 8 mm and the height of its micropatterned structures is 200 µm. The total height of the FL is 50 µm and the height of its micropatterned structures is about 30 µm, with a flexible PDMS top layer of about 20 µm, also referred to as membrane, whose controlled movement constitutes the operating principle of our culture chamber (Fig. 1b).

The top and bottom layers serve different but complementary functions. On the one hand, the CL is filled with water and closed at one end. By using a flow controller to fine-tune the pressure inside the CL, the PDMS membrane is lowered by pneumatic inflation until it contacts the underlying glass coverslip. Lowering of the membrane creates a contact zone between PDMS and glass, suitable for microbial cell containment and two-dimensional growth (Fig. 1b). Using water instead of air to actuate the CL improves the imaging quality, since the refractive indexes of water and PDMS are similar. On the other hand, the FL is continuously perfused with growth medium, providing nutrients to the cells trapped between the PDMS membrane and the glass coverslip. Importantly, the CL is composed of a standard (1:5) ratio between PDMS crosslinker and pre-polymer, as opposed to the FL, composed of a (1:20) ratio, in order to bond the two layers by migration of the crosslinking agent⁵⁷. In addition, we found that treating the FL with oxygen plasma improves the penetration of liquid culture medium (see also "Methods"), due to increased hydrophilicity⁵⁸.

We fabricated a first microsystem prototype, by arraying eight culture chambers in a row, fed by a common flow and controlled by a single CL, laid out perpendicular to the FL (Supplementary Fig. 1a). To prevent accumulation of debris in the FL and clogging of the system, we also integrated PDMS pillars of 100- and 20-µm diameter at the inlet port of the FL. To test the membrane motion, we perfused the FL with a solution of fluorescein isothiocyanate (FITC) at a constant rate of 150 µL/h, measured at the outlet, bringing the membrane to a convex position relative to the coverslip. Next, we applied a gradient of pressures to the CL, and acquired sequential images of the chamber at increasing pressure steps of one mbar per image (Fig. 1c). By measuring the fluorescence profile along the membrane diameter, we found that movement of the membrane from convex, to flat, to concave, is associated with a decrease in fluorescence (Fig. 1d; Supplementary Video 1). Starting from the point of maximum convexity of the membrane, an increment of 5 mbar was sufficient to flatten the membrane, and a further increment of 20 mbar generated a contact zone of approximately 250 µm in diameter (Fig. 1d; Supplementary Fig. 1b). This contact area is sufficiently large to visualize the growth of several bacterial microcolonies, and leaves sufficient room for lateral flow of growth medium. The unusual aspect ratio between diameter and height of the lower microchamber (33:1), over three fold above the manufacturing standards typically used to prevent collapse of microstructures⁵¹, enabled us to use actuation pressures in the CL more than one order of magnitude lower than those used in conventional micromechanical valves⁵⁹. Our ability to work at low pressure in the CL was likely to be beneficial to cellular viability.

Next, we used the FL to load a suspension of *Mycobacterium smegmatis* cells, a non-pathogenic tuberculosis model that duplicates in vitro approximately every 3 h. To identify the best culture conditions, we tested different combinations of pressure and flow rate in the CL and in the FL, and measured the doubling time of microcolonies within the two-dimensional growth area (Fig. 1e). We found that several combinations are compatible with bacterial survival and normal doubling time. However, a drastic decrease of flow rate in the FL impaired bacterial growth. This was likely due to decreased nutrients influx and possibly to excessive pressure exerted on the cells from the CL. Based on these results, we were able to choose the operating conditions of the device that are most compatible with mycobacterial growth, namely, high flow rate in the FL (150 µL/h) and moderate pressure in the CL (30 mbar). Under these conditions, we not only measured similar generation times between separate microfluidic culture chambers (Fig. 1f), implying good reproducibility across the system, but also found that mycobacterial growth was comparable to that obtained with our former single-condition microfluidic device, which is based on a different working principle and operates in the absence of shear stress¹³,⁴⁴. In addition, we calculated that the growth rate of microcolonies was independent of their location in the chamber (Fig. 1g; Supplementary Fig. 1b), suggesting homogeneous growth throughout the two-dimensional growth area. Overall, these data imply that bacteria do not experience apparent stress inside the microchambers, and that shear stress is also likely to be negligible in our system.

In conclusion, we constructed a modular microfluidic culture chamber based on a hydro-pneumatic mechanism. We also proved that this device can trap live bacterial cells, and can be used to monitor their two-dimensional growth.

**Evolution of the microfluidic culture chamber into a multi-condition platform.** To increase the versatility of our microfluidic culture chamber, we arrayed 20 of them, generating five separate groups of four microchambers (Fig. 2a, Supplementary Data 1), and we connected each group to the outlet branch of a classic microfluidic dilution tree⁵⁶,⁶⁰,⁶¹. In this configuration, the CL consists of a straight microchannel that widens in correspondence of each microchamber, ensuring the same pneumatic mechanism of the single functional unit (Fig. 1). The resulting five-condition platform is supervised by a computer-driven multichannel flow-controller. The first channel regulates a tripartite tubing, by applying pressure into three separate bottles, containing either culture medium alone or the highest concentration of the molecule to be tested dissolved in culture medium (Fig. 2a; Supplementary Fig. 2a). In addition, a bidirectional switch valve is placed downstream of two bottles to
<figure><img src="image_4.png" /><figcaption>b</figcaption></figure>

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

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

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

<figure><img src="image_8.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_13.png" /></figure>
<figure><img src="image_3.png" /><figcaption>b</figcaption></figure>

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

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

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

<figure><img src="image_7.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_12.png" /></figure>
►**Figure 1.** Operating principle and validation of the microfluidic culture chamber. (a,b) Schematics of the microfluidic culture chamber. Axonometric projection (a) and cross section (b) of the different layers constituting the chamber. From top to bottom are the control layer (CL), the membrane (M) integrated on the top of the flow layer (FL), and the glass coverslip (G). Magenta rods represent bacteria being trapped between M and G as the pressure in the CL increases. Bacterial growth is ensured by constant flow of the culture medium in the FL, and is monitored through an inverted microscope. Numbers represent the dimensions of the device in µm. (c) Representative fluorescence image series taken from the bottom of the microfluidic chamber as the pressure in the CL increases, while a FITC solution (100 µM) is injected at a constant flow rate of 150 µL/h. Numbers represent incremental pressure steps of one mbar per image in the CL, starting from 55 mbar set in the flow controller. Scale bar 285 µm. (d) Fluorescence profile along the diameter of the microfluidic chamber at a constant flow rate of 150 µL/h in the FL and increasing the pressure in the CL by one mbar at a time, as in C. Fluorescence is expressed as mean ± SEM over 1690 sections (N=16 chambers, from two independent experiments). CL and FL are color-coded as in A and B. (e) Bubble plot of *M. smegmatis* doubling time as a function of the flow rate in the FL and pressure in the CL (N=24 microcolonies, from two independent experiments). (f) Comparison of *M. smegmatis* microcolony generation time in eight culture chambers from the same device, compared to our control (CT) device⁴⁴. Black lines indicate mean ± SD (33 ≤ N ≤ 39, from three independent experiments). No significant difference is found between microcolonies growing in independent chambers, by one-way ANOVA. (g) Pearson correlation between *M. smegmatis* microcolony grow rate as a function of its position in the culture chamber. The absence of correlation suggests that bacteria behave similarly at different locations within the two-dimensional growth area.

switch between conditions. The second channel regulates a water supply for the actuation of the CL. The third channel is connected to a sealed waste receptacle, and used to maintain a stable flow rate inside the FL, in association with a flowmeter placed downstream of the outlet port (Fig. 2a).

We tested whether the five-condition platform is suitable for monolayer growth of both fast-growing *M. smegmatis* and slow-growing *Mycobacterium tuberculosis*. To this aim, we established a specific sequence of actions to prepare, assemble and load the different components of our platform, complying with the standards for BSL3 containment (Supplementary Fig. 2a-c). The FL was loaded with the bacterial suspension from a secondary outlet tubing. Following a short period of incubation, the control layer was actuated to block bacterial cells between the PDMS membrane and the coverslip (Supplementary Fig. 2b). After checking for leaks, the reservoir and waste bottles were connected to the device and to the flow controller, and the whole platform was mounted on the microscope stage (Fig. 2a, Supplementary Figs. 2c, 3a,b). By carrying out time lapse microscopy (Fig. 2b,c; Supplementary Fig. 3b), we not only confirmed that both mycobacterial species grew at standard growth rate within different microchambers of the five-condition platform, but also proved its stability for long-term imaging (Fig. 2d,e).

**Generation of different gradient patterns to mimic the dynamics of drug exposure.** To generate a gradient of concentrations in our five-condition platform (Fig. 2a), we inject two different solutions from the inlet ports, reaching two independent serpentine-shaped microchannels. The latter bifurcate into an outer serpentine, where the solutions remain at their original concentration, and an inner serpentine that acts as a mixer, where the incoming solutions are blended. Each of the three serpentines bifurcates, originating two outer serpentines and two mixers, which further bifurcate and give rise to the last dilution level with two outer serpentines and three inner mixers. The flow in the serpentines is the same at each dilution level, while the flow in the connecting sections at each bifurcation is greater towards the outside of the system and smaller towards the inside⁶⁰. The concentration gradient resulting from this network structure is linear. From the last dilution level, the solutions flow into five separate clusters (C1-C5) of four microchambers each. Finally, each group of chambers is connected to the same outlet port, aiming to maximize flow stability in the whole microsystem. We illustrate the operation of our platform by sequential mixing of dyes (Fig. 3a).

We quantified the gradient by injecting a solution of FITC dissolved into 7H9 medium into the first inlet port, and non-fluorescent 7H9 medium into the second inlet port (Supplementary Fig. 4a). At steady state, the measured gradient minimally differed from the theoretical expectations⁶⁰, in that the FITC solution was undiluted in C1, diluted to about 72% in C2, 48% in C3, 25% in C4, and was absent in C5 (Fig. 3b; Supplementary Fig. 4a). Aiming to mimic a typical PK profile (Supplementary Fig. 4b), we also generated pulse-shaped FITC gradients (Fig. 3c), whose concentrations at peak differed from those obtained with static gradients. In particular, the FITC solution at the time of maximum concentration was undiluted in C1, diluted to about 82% in C2, 54% in C3, 21% in C4, and was absent in C5. Static and dynamic gradients (Fig. 3b,c) were produced by alternating the perfusion of non-fluorescent and fluorescent solutions for different time intervals, by means of a bidirectional switch valve. Interestingly, we could obtain pulsing PK-mimetic gradients only when capacitors were integrated upstream of the inlet ports (Fig. 2a; Supplementary Fig. 2c). In summary, by integrating a dilution-tree network with our microfluidic culture chambers, we created a microfluidic platform suitable to titrate molecules over space and time.

**Moxifloxacin acts in a dose- and time-dependent manner against *M. smegmatis* microcolonies.** To test the operation of our platform, we examined the effect of different drug concentrations on the growth rate of *M. smegmatis*. We focused on moxifloxacin, a fluoroquinolone drug targeting the type II topoisomerase DNA gyrase, and critical for antitubercular therapy against multidrug-resistant strains³³. To monitor the effect on the target, we generated a red fluorescent reporter of GyrA, by fusing the *mCherry* gene in frame
►**Figure 1.** Operating principle and validation of the microfluidic culture chamber. (a,b) Schematics of the microfluidic culture chamber. Axonometric projection (a) and cross section (b) of the different layers constituting the chamber. From top to bottom are the control layer (CL), the membrane (M) integrated on the top of the flow layer (FL), and the glass coverslip (G). Magenta rods represent bacteria being trapped between M and G as the pressure in the CL increases. Bacterial growth is ensured by constant flow of the culture medium in the FL, and is monitored through an inverted microscope. Numbers represent the dimensions of the device in µm. (c) Representative fluorescence image series taken from the bottom of the microfluidic chamber as the pressure in the CL increases, while a FITC solution (100 µM) is injected at a constant flow rate of 150 µL/h. Numbers represent incremental pressure steps of one mbar per image in the CL, starting from 55 mbar set in the flow controller. Scale bar 285 µm. (d) Fluorescence profile along the diameter of the microfluidic chamber at a constant flow rate of 150 µL/h in the FL and increasing the pressure in the CL by one mbar at a time, as in C. Fluorescence is expressed as mean ± SEM over 1690 sections (N=16 chambers, from two independent experiments). CL and FL are color-coded as in A and B. (e) Bubble plot of *M. smegmatis* doubling time as a function of the flow rate in the FL and pressure in the CL (N=24 microcolonies, from two independent experiments). (f) Comparison of *M. smegmatis* microcolony generation time in eight culture chambers from the same device, compared to our control (CT) device⁴⁴. Black lines indicate mean ± SD (33 ≤ N ≤ 39, from three independent experiments). No significant difference is found between microcolonies growing in independent chambers, by one-way ANOVA. (g) Pearson correlation between *M. smegmatis* microcolony grow rate as a function of its position in the culture chamber. The absence of correlation suggests that bacteria behave similarly at different locations within the two-dimensional growth area.

switch between conditions. The second channel regulates a water supply for the actuation of the CL. The third channel is connected to a sealed waste receptacle, and used to maintain a stable flow rate inside the FL, in association with a flowmeter placed downstream of the outlet port (Fig. 2a).

We tested whether the five-condition platform is suitable for monolayer growth of both fast-growing *M. smegmatis* and slow-growing *Mycobacterium tuberculosis*. To this aim, we established a specific sequence of actions to prepare, assemble and load the different components of our platform, complying with the standards for BSL3 containment (Supplementary Fig. 2a-c). The FL was loaded with the bacterial suspension from a secondary outlet tubing. Following a short period of incubation, the control layer was actuated to block bacterial cells between the PDMS membrane and the coverslip (Supplementary Fig. 2b). After checking for leaks, the reservoir and waste bottles were connected to the device and to the flow controller, and the whole platform was mounted on the microscope stage (Fig. 2a, Supplementary Figs. 2c, 3a,b). By carrying out time lapse microscopy (Fig. 2b,c; Supplementary Fig. 3b), we not only confirmed that both mycobacterial species grew at standard growth rate within different microchambers of the five-condition platform, but also proved its stability for long-term imaging (Fig. 2d,e).

**Generation of different gradient patterns to mimic the dynamics of drug exposure.** To generate a gradient of concentrations in our five-condition platform (Fig. 2a), we inject two different solutions from the inlet ports, reaching two independent serpentine-shaped microchannels. The latter bifurcate into an outer serpentine, where the solutions remain at their original concentration, and an inner serpentine that acts as a mixer, where the incoming solutions are blended. Each of the three serpentines bifurcates, originating two outer serpentines and two mixers, which further bifurcate and give rise to the last dilution level with two outer serpentines and three inner mixers. The flow in the serpentines is the same at each dilution level, while the flow in the connecting sections at each bifurcation is greater towards the outside of the system and smaller towards the inside⁶⁰. The concentration gradient resulting from this network structure is linear. From the last dilution level, the solutions flow into five separate clusters (C1-C5) of four microchambers each. Finally, each group of chambers is connected to the same outlet port, aiming to maximize flow stability in the whole microsystem. We illustrate the operation of our platform by sequential mixing of dyes (Fig. 3a).

We quantified the gradient by injecting a solution of FITC dissolved into 7H9 medium into the first inlet port, and non-fluorescent 7H9 medium into the second inlet port (Supplementary Fig. 4a). At steady state, the measured gradient minimally differed from the theoretical expectations⁶⁰, in that the FITC solution was undiluted in C1, diluted to about 72% in C2, 48% in C3, 25% in C4, and was absent in C5 (Fig. 3b; Supplementary Fig. 4a). Aiming to mimic a typical PK profile (Supplementary Fig. 4b), we also generated pulse-shaped FITC gradients (Fig. 3c), whose concentrations at peak differed from those obtained with static gradients. In particular, the FITC solution at the time of maximum concentration was undiluted in C1, diluted to about 82% in C2, 54% in C3, 21% in C4, and was absent in C5. Static and dynamic gradients (Fig. 3b,c) were produced by alternating the perfusion of non-fluorescent and fluorescent solutions for different time intervals, by means of a bidirectional switch valve. Interestingly, we could obtain pulsing PK-mimetic gradients only when capacitors were integrated upstream of the inlet ports (Fig. 2a; Supplementary Fig. 2c). In summary, by integrating a dilution-tree network with our microfluidic culture chambers, we created a microfluidic platform suitable to titrate molecules over space and time.

**Moxifloxacin acts in a dose- and time-dependent manner against *M. smegmatis* microcolonies.** To test the operation of our platform, we examined the effect of different drug concentrations on the growth rate of *M. smegmatis*. We focused on moxifloxacin, a fluoroquinolone drug targeting the type II topoisomerase DNA gyrase, and critical for antitubercular therapy against multidrug-resistant strains³³. To monitor the effect on the target, we generated a red fluorescent reporter of GyrA, by fusing the *mCherry* gene in frame
<figure><img src="image_14.png" /><figcaption>**Figure 2.** Five-condition platform set-up and validation. (a) Schematic of the five-condition microfluidic device connected to a multi-channel microfluidic flow controller (Fluigent). The first channel is connected to three independent reservoirs, in turn connected to a software-controlled bidirectional position valve to switch between conditions. The reservoirs are connected to two inlet ports (c,d) of the FL (yellow pattern), and two capacitors (gray cylinders) are integrated upstream of the inlet ports. A software-driven flowmeter is integrated downstream the outlet port of the FL (e), which flows into a sealed trash. The outlet tubing is bifurcated via a Y-connector, for bacteria injection. The second channel is connected to the water reservoir, which actuates the membrane from the CL (blue pattern). Filters (0.2 µm) are used at the entry of reservoir and waste bottles to prevent contamination of both fluidics and environment. Time-lapse imaging is carried out from the glass bottom of the device, using a motorized wide-field inverted microscope. (b,c) Representative phase-contrast time-lapse image series of exponentially growing *M. smegmatis* (b) and *M. tuberculosis* (c) inside the five-condition device. Cells were imaged at 20-min (b) or 3-h intervals (c), and numbers represent hours. Scale bars 10 µm. (d,e) Microcolony growth rates of *M. smegmatis* (d) and *M. tuberculosis* (e) growing inside different chambers of the five-condition device, and compared to our CT device⁴⁴. Black lines indicate mean ± SD (15 ≤ N ≤ 60, from at least two independent experiments per strain). No significant difference is found by one-way ANOVA followed by Kruskal-Wallis multiple comparisons test.</figcaption></figure>

with the *gyrA* chromosomal copy (Supplementary Fig. 5a). We seeded log-phase cells of *M. smegmatis* GyrA-mCherry reporter inside the 5-condition platform (Fig. 2a). First, we grew the bacteria under constant perfusion of fresh 7H9 growth medium for 6 h (about two generation times), next we injected moxifloxacin from one inlet
<figure><img src="image_13.png" /><figcaption>**Figure 2.** Five-condition platform set-up and validation. (a) Schematic of the five-condition microfluidic device connected to a multi-channel microfluidic flow controller (Fluigent). The first channel is connected to three independent reservoirs, in turn connected to a software-controlled bidirectional position valve to switch between conditions. The reservoirs are connected to two inlet ports (c,d) of the FL (yellow pattern), and two capacitors (gray cylinders) are integrated upstream of the inlet ports. A software-driven flowmeter is integrated downstream the outlet port of the FL (e), which flows into a sealed trash. The outlet tubing is bifurcated via a Y-connector, for bacteria injection. The second channel is connected to the water reservoir, which actuates the membrane from the CL (blue pattern). Filters (0.2 µm) are used at the entry of reservoir and waste bottles to prevent contamination of both fluidics and environment. Time-lapse imaging is carried out from the glass bottom of the device, using a motorized wide-field inverted microscope. (b,c) Representative phase-contrast time-lapse image series of exponentially growing *M. smegmatis* (b) and *M. tuberculosis* (c) inside the five-condition device. Cells were imaged at 20-min (b) or 3-h intervals (c), and numbers represent hours. Scale bars 10 µm. (d,e) Microcolony growth rates of *M. smegmatis* (d) and *M. tuberculosis* (e) growing inside different chambers of the five-condition device, and compared to our CT device⁴⁴. Black lines indicate mean ± SD (15 ≤ N ≤ 60, from at least two independent experiments per strain). No significant difference is found by one-way ANOVA followed by Kruskal-Wallis multiple comparisons test.</figcaption></figure>

with the *gyrA* chromosomal copy (Supplementary Fig. 5a). We seeded log-phase cells of *M. smegmatis* GyrA-mCherry reporter inside the 5-condition platform (Fig. 2a). First, we grew the bacteria under constant perfusion of fresh 7H9 growth medium for 6 h (about two generation times), next we injected moxifloxacin from one inlet
<figure><img src="image_15.png" /><figcaption>**Figure 3.** Characterization of gradient formation in the five-condition platform. (a) Picture of the five-condition device filled with dyes from two separate inlet ports (IN1 and IN2) of the FL, in the absence of CL. Blue and yellow dyes form a color gradient in the dilution tree, via serpentine mixers (M), reach five independent groups of four chambers (C1–C5), and flow into a common outlet (OUT). Scale bar 10 mm. (b,c) Formation of gradients through progressive mixing of a FITC solution (100 µM) injected into IN1 with a non-fluorescent solution injected into IN2. Mean fluorescence is measured from each culture chamber. Steady (b) and pulsing (c) gradient patterns of different duration are shown. Filled lines and error bars indicate mean ± SEM (N=2 or 4). Colors indicate different groups of chambers, color coded as in (a).</figcaption></figure>

port to generate a pulsing gradient, and finally we switched back to drug-free medium (Fig. 4a; Supplementary Video 2). To mimic the PK profile of moxifloxacin⁶², we empirically formed the gradient, producing a positively skewed injection pattern (Fig. 4b; Supplementary Fig. 4b), and scaling its duration to the generation time of *M. smegmatis*. At one end of the device we left a group of chambers without drug to monitor the normal growth, and at the opposite end we tested the highest drug concentration (250 ng/mL; 5X-MIC), reaching the peak of the pulse about one and a half hour after injection (Fig. 4b). The concentrations of moxifloxacin present in the three middle groups of microchambers were estimated from the dilution factors measured by fluorescence (Fig. 3c). We also verified that moxifloxacin was not significantly adsorbed by PDMS, by measuring its inlet and outlet concentration via spectrofluorimetry over time (Supplementary Fig. 4c,d). We found that at the beginning of infusion, the outgoing concentration is lower than the incoming concentration by about 15%, suggesting an initial adsorption of moxifloxacin by PDMS. However, in less than an hour, the incoming and outgoing drug levels become comparable. In addition, when fresh medium was injected in lieu of medium containing moxifloxacin, the presence of drug quickly became undetectable, suggesting that PDMS does not retain the drug.

We implemented a semi-automated microcolony analysis workflow (Supplementary Fig. 5b, Supplementary Data 2–7), enabling us to measure size and fluorescence of several microcolonies per condition, for the entire
<figure><img src="image_14.png" /><figcaption>**Figure 3.** Characterization of gradient formation in the five-condition platform. (a) Picture of the five-condition device filled with dyes from two separate inlet ports (IN1 and IN2) of the FL, in the absence of CL. Blue and yellow dyes form a color gradient in the dilution tree, via serpentine mixers (M), reach five independent groups of four chambers (C1–C5), and flow into a common outlet (OUT). Scale bar 10 mm. (b,c) Formation of gradients through progressive mixing of a FITC solution (100 µM) injected into IN1 with a non-fluorescent solution injected into IN2. Mean fluorescence is measured from each culture chamber. Steady (b) and pulsing (c) gradient patterns of different duration are shown. Filled lines and error bars indicate mean ± SEM (N=2 or 4). Colors indicate different groups of chambers, color coded as in (a).</figcaption></figure>

port to generate a pulsing gradient, and finally we switched back to drug-free medium (Fig. 4a; Supplementary Video 2). To mimic the PK profile of moxifloxacin⁶², we empirically formed the gradient, producing a positively skewed injection pattern (Fig. 4b; Supplementary Fig. 4b), and scaling its duration to the generation time of *M. smegmatis*. At one end of the device we left a group of chambers without drug to monitor the normal growth, and at the opposite end we tested the highest drug concentration (250 ng/mL; 5X-MIC), reaching the peak of the pulse about one and a half hour after injection (Fig. 4b). The concentrations of moxifloxacin present in the three middle groups of microchambers were estimated from the dilution factors measured by fluorescence (Fig. 3c). We also verified that moxifloxacin was not significantly adsorbed by PDMS, by measuring its inlet and outlet concentration via spectrofluorimetry over time (Supplementary Fig. 4c,d). We found that at the beginning of infusion, the outgoing concentration is lower than the incoming concentration by about 15%, suggesting an initial adsorption of moxifloxacin by PDMS. However, in less than an hour, the incoming and outgoing drug levels become comparable. In addition, when fresh medium was injected in lieu of medium containing moxifloxacin, the presence of drug quickly became undetectable, suggesting that PDMS does not retain the drug.

We implemented a semi-automated microcolony analysis workflow (Supplementary Fig. 5b, Supplementary Data 2–7), enabling us to measure size and fluorescence of several microcolonies per condition, for the entire
<figure><img src="image_16.png" /><figcaption>**Figure 4.** Effect of a pulsing moxifloxacin gradient on *M. smegmatis* GyrA-mCherry reporter. (a,b) Time-lapse microscopy of exponentially growing *M. smegmatis* seeded into the five-condition device and stressed with a pulsing gradient of moxifloxacin. Time-lapse image series representative of each chamber group (a). Phase-contrast and fluorescence images were acquired every 20 min and numbers represent hours. Scale bar 10 µm. Circles indicate the absence or the presence of different drug concentrations, color-coded in (b). Concentrations shown in the inset relative to the MIC (dashed line) refer to the peak concentration reached during moxifloxacin exposure, according to the dilution factors measured in Fig. 3c. (c,d) Normalized microcolony area (c) and GyrA-mCherry fluorescence (d). Space between dotted lines indicates drug exposure period. Lines and error bars are color-coded as in (b) and represent mean ± SEM (10≤N≤42 microcolonies, from three independent experiments). Significance by 2-way ANOVA followed by Dunnett's multiple comparisons test: ***P=0.0005; ****P<0.0001. (e) Microcolony growth rates measured for each condition color-coded as in (b), during different experimental phases, before (pre-drug), during (drug) and after moxifloxacin exposure (recovery). Lines and error bars represent mean ± SD (10≤N≤42). Significance by one-way ANOVA followed by Dunnett's multiple comparisons test: *P=0.01; **P=0.001; ***P=0.0003; ****P<0.0001. (f) Pearson r correlation between different PK parameters and averaged microcolonies growth rates (10≤N≤42), measured over the entire duration of the experiment, for each condition color-coded as in (b). Left y-axis: ratio of area under the concentration-time curve to MIC (AUC/MIC), and time above MIC (T>MIC); right y-axis: ratio of peak concentration to MIC (Peak/MIC).</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 4.** Effect of a pulsing moxifloxacin gradient on *M. smegmatis* GyrA-mCherry reporter. (a,b) Time-lapse microscopy of exponentially growing *M. smegmatis* seeded into the five-condition device and stressed with a pulsing gradient of moxifloxacin. Time-lapse image series representative of each chamber group (a). Phase-contrast and fluorescence images were acquired every 20 min and numbers represent hours. Scale bar 10 µm. Circles indicate the absence or the presence of different drug concentrations, color-coded in (b). Concentrations shown in the inset relative to the MIC (dashed line) refer to the peak concentration reached during moxifloxacin exposure, according to the dilution factors measured in Fig. 3c. (c,d) Normalized microcolony area (c) and GyrA-mCherry fluorescence (d). Space between dotted lines indicates drug exposure period. Lines and error bars are color-coded as in (b) and represent mean ± SEM (10≤*N*≤42 microcolonies, from three independent experiments). Significance by 2-way ANOVA followed by Dunnett's multiple comparisons test: ***P=0.0005; ****P<0.0001. (e) Microcolony growth rates measured for each condition color-coded as in (b), during different experimental phases, before (pre-drug), during (drug) and after moxifloxacin exposure (recovery). Lines and error bars represent mean ± SD (10≤*N*≤42). Significance by one-way ANOVA followed by Dunnett's multiple comparisons test: *P*=0.01; **P*=0.001; ***P=0.0003; ****P<0.0001. (f) Pearson *r* correlation between different PK parameters and averaged microcolonies growth rates (10≤*N*≤42), measured over the entire duration of the experiment, for each condition color-coded as in (b). Left y-axis: ratio of area under the concentration-time curve to MIC (AUC/MIC), and time above MIC (T>MIC); right y-axis: ratio of peak concentration to MIC (Peak/MIC).</figcaption></figure>
duration of the experiment (Fig. 4c,d). While in the absence of moxifloxacin bacteria exhibited exponential growth, we measured a significant and progressive decrease of microcolony size at increasing concentrations of drug, and a rapid growth arrest at the highest concentration (Fig. 4c). Compared to untreated cells, we observed significant target inhibition at the highest concentration of moxifloxacin, and target induction at lower concentrations (Fig. 4d), with implications for the induction of DNA damage response mechanisms⁶³ and possible consequences for the onset of mutations. During the first 6 h of drug exposure, we did not measure significant changes in growth rate among microcolonies growing in the absence or presence of different drug concentrations, whereas the growth rate started to slow down significantly during the last 6 h of drug exposure and continued to decrease during the recovery phase in the absence of drug (Fig. 4e). In addition, we observed marked heterogeneity in growth rate among clonal microcolonies exposed to the same drug concentration (Fig. 4e). Finally, we checked whether there was a correlation between growth rate and the PK parameters simulated in vitro (Fig. 4b; Supplementary Fig. 4b). Interestingly, we measured an inverse relationship between the growth rate and both dose (R² = 89.8%) and time above the MIC (R² = 90.8%). In contrast, we found no significant correlation between the growth rate and the area under the curve (AUC) to MIC ratio (Fig. 4f). To conclude, we applied our platform to analyze PD properties and PD-PK relationships at high cellular resolution in space and time, and found that moxifloxacin exerts its effect on *M. smegmatis* microcolonies in a dose- and time-dependent manner.

**Response heterogeneity among single *M. tuberculosis* cells subjected to a moxifloxacin gradient.** To test the stability of our platform for long-term imaging, we examined the effect of different concentrations of moxifloxacin against *M. tuberculosis* by multi-phasic time-lapse microscopy. We created the red fluorescent reporter strain of GyrA in *M. tuberculosis* (Supplementary Fig. 6a), and further labeled it with a constitutively expressed green fluorescent protein (GFP<sub>cyt</sub>). While GyrA-mCherry provided information about the moxifloxacin target, the cytosolic marker helped us to identify individual bacilli and to recognize their loss of viability or lysis from gradual or abrupt fluorescence decrease, respectively. We first seeded log-phase GFP<sub>cyt</sub>-GyrA-mCherry bacilli inside the 5-condition platform, and grew them in fresh medium for almost 4 days (about four generation times). Next, we injected the maximum moxifloxacin concentration (250 ng/mL; 5X-MIC) from one inlet port to generate a stable gradient for 2 days, leaving the opposite end of the device devoid of drug. Finally, we switched back to drug-free medium for 5 days to monitor the recovery phase (Supplementary Fig. 6b; Supplementary Video 3). Also in this case, the concentrations of moxifloxacin present in the three middle groups of microchambers were estimated from the dilution factors previously measured by fluorescence (Fig. 3b).

To increase the analytical resolution, we quantified individual bacilli and their fate over time. In the absence of drug, the number of bacilli grew exponentially until they saturated the field of view, leading to a decrease of the division rate and making cell counting impractical (Fig. 5a). At increasing concentrations of moxifloxacin, *M. tuberculosis* cell number plateaued, and the rate of division decreased until it stopped (Fig. 5a). Interestingly, cells subjected to lower drug concentrations exhibited significant increase in size, consistent with induction of DNA damage response¹³,⁶³. Conversely, bacilli stressed with the highest concentration of moxifloxacin (250 ng/mL; 5X-MIC) underwent abrupt growth arrest without significant changes in their size (Fig. 5b). Overall, these results imply that small differences in drug concentrations, close to the MIC, could have important implications for the induction of adaptive stress responses in individual cells, and for the onset of drug resistance. As drug concentration increased, we also observed a progressive decrease in GFP<sub>cyt</sub> fluorescence, concurrent with an increase in GFP<sub>cyt</sub>-dim bacilli, suggestive of gradual loss of cell viability (Fig. 5c-e). In contrast, the presence of residual GFP<sub>cyt</sub> fluorescence in intact bacilli that do not grow may be attributed to a persistent-like state, which could eventually revert to active growth under appropriate conditions and in due time. Thus, in this experimental set-up, we regarded as survivors both cells that resumed growth after drug removal and those that retained weak GFP<sub>cyt</sub> signal in the absence of growth. The fraction of these viable cells decreased considerably starting from a dose of moxifloxacin equal to 2.4X-MIC, when GFP<sub>cyt</sub>-dim cells were prevailing but only a few cells lysed (Fig. 5f). We measured a half maximal inhibitory concentration (IC50) just below the MIC value, while the maximum drug efficacy was achieved at 2.4X-MIC (Fig. 5g). This was also consistent with the effect on the target, in that only cells treated with the lowest concentration (1.25X-MIC) exhibited an increase of GyrA-mCherry, consistent with induction of DNA damage response by fluoroquinolones¹³,⁶³. In contrast, GyrA-mCherry expression decreased at higher concentrations (Fig. 5h; Supplementary Fig. 6c), implying dose-dependent inhibition.

Here, we report the first proof of concept of a dynamic dose-response microfluidic assay with a slow-growing and clinically relevant pathogen. We proved the functionality of the platform for long-term microscopy in BSL3, and the ability to visualize diverse single-cell responses according to small changes in drug concentration. Understanding of these dynamics is crucial for predicting adaptive responses of the pathogen in vivo under fluctuating drug concentrations. In conclusion, our platform has the potential to improve preclinical characterization of drug efficacy and the understanding of the mechanism of action. This approach can complement conventional bulk assays, offering direct information on the dynamics of killing and survival of drug-evading cells, responsible for the loss of efficacy of antimicrobial agents.

## Discussion

The ability of some pathogenic bacteria to tolerate prolonged exposure to antimicrobials can be conducive to antimicrobial resistance²⁰⁻²⁵, and poses major limitations to the management of infections¹,⁶⁴. Although this phenomenon has been known for a long time³,⁴, direct identification and analysis of drug-persistent subpopulations are troublesome, hindering preclinical identification of drugs effective against them. As opposed to bulk-cell studies, spatiotemporal analysis of single bacterial cells cultured inside microfluidic devices offers insights into the dynamics of cellular responses to antimicrobial treatment, and the formation of drug-persistent cells¹¹,¹³,²⁰,⁴⁰,⁴³.
duration of the experiment (Fig. 4c,d). While in the absence of moxifloxacin bacteria exhibited exponential growth, we measured a significant and progressive decrease of microcolony size at increasing concentrations of drug, and a rapid growth arrest at the highest concentration (Fig. 4c). Compared to untreated cells, we observed significant target inhibition at the highest concentration of moxifloxacin, and target induction at lower concentrations (Fig. 4d), with implications for the induction of DNA damage response mechanisms⁶³ and possible consequences for the onset of mutations. During the first 6 h of drug exposure, we did not measure significant changes in growth rate among microcolonies growing in the absence or presence of different drug concentrations, whereas the growth rate started to slow down significantly during the last 6 h of drug exposure and continued to decrease during the recovery phase in the absence of drug (Fig. 4e). In addition, we observed marked heterogeneity in growth rate among clonal microcolonies exposed to the same drug concentration (Fig. 4e). Finally, we checked whether there was a correlation between growth rate and the PK parameters simulated in vitro (Fig. 4b; Supplementary Fig. 4b). Interestingly, we measured an inverse relationship between the growth rate and both dose (R² = 89.8%) and time above the MIC (R² = 90.8%). In contrast, we found no significant correlation between the growth rate and the area under the curve (AUC) to MIC ratio (Fig. 4f). To conclude, we applied our platform to analyze PD properties and PD-PK relationships at high cellular resolution in space and time, and found that moxifloxacin exerts its effect on *M. smegmatis* microcolonies in a dose- and time-dependent manner.

**Response heterogeneity among single *M. tuberculosis* cells subjected to a moxifloxacin gradient.** To test the stability of our platform for long-term imaging, we examined the effect of different concentrations of moxifloxacin against *M. tuberculosis* by multi-phasic time-lapse microscopy. We created the red fluorescent reporter strain of GyrA in *M. tuberculosis* (Supplementary Fig. 6a), and further labeled it with a constitutively expressed green fluorescent protein (GFP<sub>cyt</sub>). While GyrA-mCherry provided information about the moxifloxacin target, the cytosolic marker helped us to identify individual bacilli and to recognize their loss of viability or lysis from gradual or abrupt fluorescence decrease, respectively. We first seeded log-phase GFP<sub>cyt</sub>-GyrA-mCherry bacilli inside the 5-condition platform, and grew them in fresh medium for almost 4 days (about four generation times). Next, we injected the maximum moxifloxacin concentration (250 ng/mL; 5X-MIC) from one inlet port to generate a stable gradient for 2 days, leaving the opposite end of the device devoid of drug. Finally, we switched back to drug-free medium for 5 days to monitor the recovery phase (Supplementary Fig. 6b; Supplementary Video 3). Also in this case, the concentrations of moxifloxacin present in the three middle groups of microchambers were estimated from the dilution factors previously measured by fluorescence (Fig. 3b).

To increase the analytical resolution, we quantified individual bacilli and their fate over time. In the absence of drug, the number of bacilli grew exponentially until they saturated the field of view, leading to a decrease of the division rate and making cell counting impractical (Fig. 5a). At increasing concentrations of moxifloxacin, *M. tuberculosis* cell number plateaued, and the rate of division decreased until it stopped (Fig. 5a). Interestingly, cells subjected to lower drug concentrations exhibited significant increase in size, consistent with induction of DNA damage response¹³,⁶³. Conversely, bacilli stressed with the highest concentration of moxifloxacin (250 ng/mL; 5X-MIC) underwent abrupt growth arrest without significant changes in their size (Fig. 5b). Overall, these results imply that small differences in drug concentrations, close to the MIC, could have important implications for the induction of adaptive stress responses in individual cells, and for the onset of drug resistance. As drug concentration increased, we also observed a progressive decrease in GFP<sub>cyt</sub> fluorescence, concurrent with an increase in GFP<sub>cyt</sub>-dim bacilli, suggestive of gradual loss of cell viability (Fig. 5c-e). In contrast, the presence of residual GFP<sub>cyt</sub> fluorescence in intact bacilli that do not grow may be attributed to a persistent-like state, which could eventually revert to active growth under appropriate conditions and in due time. Thus, in this experimental set-up, we regarded as survivors both cells that resumed growth after drug removal and those that retained weak GFP<sub>cyt</sub> signal in the absence of growth. The fraction of these viable cells decreased considerably starting from a dose of moxifloxacin equal to 2.4X-MIC, when GFP<sub>cyt</sub>-dim cells were prevailing but only a few cells lysed (Fig. 5f). We measured a half maximal inhibitory concentration (IC50) just below the MIC value, while the maximum drug efficacy was achieved at 2.4X-MIC (Fig. 5g). This was also consistent with the effect on the target, in that only cells treated with the lowest concentration (1.25X-MIC) exhibited an increase of GyrA-mCherry, consistent with induction of DNA damage response by fluoroquinolones¹³,⁶³. In contrast, GyrA-mCherry expression decreased at higher concentrations (Fig. 5h; Supplementary Fig. 6c), implying dose-dependent inhibition.

Here, we report the first proof of concept of a dynamic dose-response microfluidic assay with a slow-growing and clinically relevant pathogen. We proved the functionality of the platform for long-term microscopy in BSL3, and the ability to visualize diverse single-cell responses according to small changes in drug concentration. Understanding of these dynamics is crucial for predicting adaptive responses of the pathogen in vivo under fluctuating drug concentrations. In conclusion, our platform has the potential to improve preclinical characterization of drug efficacy and the understanding of the mechanism of action. This approach can complement conventional bulk assays, offering direct information on the dynamics of killing and survival of drug-evading cells, responsible for the loss of efficacy of antimicrobial agents.

## Discussion

The ability of some pathogenic bacteria to tolerate prolonged exposure to antimicrobials can be conducive to antimicrobial resistance²⁰⁻²⁵, and poses major limitations to the management of infections¹,⁶⁴. Although this phenomenon has been known for a long time³,⁴, direct identification and analysis of drug-persistent subpopulations are troublesome, hindering preclinical identification of drugs effective against them. As opposed to bulk-cell studies, spatiotemporal analysis of single bacterial cells cultured inside microfluidic devices offers insights into the dynamics of cellular responses to antimicrobial treatment, and the formation of drug-persistent cells¹¹,¹³,²⁰,⁴⁰,⁴³.
<figure><img src="image_17.png" /></figure>

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

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

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

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

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

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

<figure><img src="image_24.png" /><figcaption>**Figure 5.** Effect of a static moxifloxacin gradient on *M. tuberculosis* GFP<sub>cyt</sub>-GyrA-mCherry reporter. (a) Assessment of *M. tuberculosis* GFP<sub>cyt</sub>-GyrA-mCherry growth during multi-phasic time-lapse microscopy. Number of bacilli (left panel) and cell division rates (right panel) in different chamber groups of the five-condition device before, during and after moxifloxacin exposure (dotted lines). Lines and error bars represent mean ± SEM (279 ≤ *N* ≤ 842 bacilli, from two independent experiments). Colors refer to the concentration of moxifloxacin used in each chambers group relative to the MIC (50 ng/mL) according to the dilution factors measured in Fig. 3b, as denoted in the inset. In the absence of drug (gray lines), cell count stops earlier due to excessive growth. (b,e,h) Single-cell size (b), GFP<sub>cyt</sub> (e) and GyrA-mCherry (h) fluorescence at different experimental time points just before (Pre) and just after (Post) treatment, and at the end of the recovery period (End). Black lines indicate mean ± SD (253 ≤ *N* ≤ 308 bacilli, from two independent experiments). Significance by one-way ANOVA followed by Kruskal-Wallis multiple comparisons test: \*P = 0.05; \*\*P < 0.005; \*\*\*P < 0.0008; \*\*\*\*P < 0.0001. (c,d) Fractions of subpopulations over time displaying bright GFP<sub>cyt</sub> (c) or dim GFP<sub>cyt</sub> (d) fluorescence. Space between dotted lines represents the period of drug exposure. Lines and error bars represent mean ± SEM (222 ≤ *N* ≤ 1327 bacilli, from two independent experiments). Significance by two-way ANOVA followed by Tukey's multiple comparisons test: \*\*P = 0.006; \*\*\*P < 0.0001. (f) Fractions of subpopulation fates at different drug concentrations (295 ≤ *N* ≤ 1327, from two independent experiments). (g) Cumulative fraction of survivors (regrowing and intact GFP<sub>cyt</sub>-bright) at the end of the recovery period as a function of moxifloxacin concentration. Symbols and error bars represent mean ± 95% CI (295 ≤ *N* ≤ 1327, from two independent experiments). Black line represents fitting of the data with a four-parameter dose-response curve: R² = 0.88; d.f. 98. The IC50 value is shown (dashed lines).</figcaption></figure>

Despite the potential of this approach, the challenge remains of how to scale single-condition devices to multi-condition platforms, while preserving device stability, user control and long-term cell viability.

Here we report the construction of a modular microfluidic cell-culture chamber that operates on a hydro-pneumatic mechanism, and is characterized by a particularly high aspect ratio. By multiplying this functional unit and combining it with a double-entry microfluidic dilution system, we derived a 5-condition microfluidic platform. The latter is composed of two inlets followed by a series of serpentine channels, which sequentially achieve titration of the incoming solutions across five independent groups of chambers. Time-resolved observation of these chambers by wide-field inverted microscopy allowed us to directly capture growth, mortality and survival rates of single mycobacterial cells before, during and after drug exposure. We tested the operation of our platform using both non-pathogenic and pathogenic mycobacterial species, which are particularly small, below 0.5 µm in thickness and between 1 and 6 µm in length, and have generation times ranging from 3 to 24 h. Mycobacterial cells were cultured within a two-dimensional growth area, generated by gentle adhesion of a thin
<figure><img src="image_16.png" /></figure>

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

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

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

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

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

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

<figure><img src="image_23.png" /><figcaption>**Figure 5.** Effect of a static moxifloxacin gradient on *M. tuberculosis* GFP<sub>cyt</sub>-GyrA-mCherry reporter. (a) Assessment of *M. tuberculosis* GFP<sub>cyt</sub>-GyrA-mCherry growth during multi-phasic time-lapse microscopy. Number of bacilli (left panel) and cell division rates (right panel) in different chamber groups of the five-condition device before, during and after moxifloxacin exposure (dotted lines). Lines and error bars represent mean ± SEM (279≤N≤842 bacilli, from two independent experiments). Colors refer to the concentration of moxifloxacin used in each chambers group relative to the MIC (50 ng/mL) according to the dilution factors measured in Fig. 3b, as denoted in the inset. In the absence of drug (gray lines), cell count stops earlier due to excessive growth. (b,e,h) Single-cell size (b), GFP<sub>cyt</sub> (e) and GyrA-mCherry (h) fluorescence at different experimental time points just before (Pre) and just after (Post) treatment, and at the end of the recovery period (End). Black lines indicate mean ± SD (253≤N≤308 bacilli, from two independent experiments). Significance by one-way ANOVA followed by Kruskal-Wallis multiple comparisons test: \*P=0.05; \*\*P<0.005; \*\*\*P<0.0008; \*\*\*\*P<0.0001. (c,d) Fractions of subpopulations over time displaying bright GFP<sub>cyt</sub> (c) or dim GFP<sub>cyt</sub> (d) fluorescence. Space between dotted lines represents the period of drug exposure. Lines and error bars represent mean ± SEM (222≤N≤1327 bacilli, from two independent experiments). Significance by two-way ANOVA followed by Tukey's multiple comparisons test: \*\*P=0.006; \*\*\*P<0.0001. (f) Fractions of subpopulation fates at different drug concentrations (295≤N≤1327, from two independent experiments). (g) Cumulative fraction of survivors (regrowing and intact GFP<sub>cyt</sub>-bright) at the end of the recovery period as a function of moxifloxacin concentration. Symbols and error bars represent mean ± 95% CI (295≤N≤1327, from two independent experiments). Black line represents fitting of the data with a four-parameter dose-response curve: R²=0.88; d.f. 98. The IC50 value is shown (dashed lines).</figcaption></figure>

Despite the potential of this approach, the challenge remains of how to scale single-condition devices to multi-condition platforms, while preserving device stability, user control and long-term cell viability.

Here we report the construction of a modular microfluidic cell-culture chamber that operates on a hydro-pneumatic mechanism, and is characterized by a particularly high aspect ratio. By multiplying this functional unit and combining it with a double-entry microfluidic dilution system, we derived a 5-condition microfluidic platform. The latter is composed of two inlets followed by a series of serpentine channels, which sequentially achieve titration of the incoming solutions across five independent groups of chambers. Time-resolved observation of these chambers by wide-field inverted microscopy allowed us to directly capture growth, mortality and survival rates of single mycobacterial cells before, during and after drug exposure. We tested the operation of our platform using both non-pathogenic and pathogenic mycobacterial species, which are particularly small, below 0.5 µm in thickness and between 1 and 6 µm in length, and have generation times ranging from 3 to 24 h. Mycobacterial cells were cultured within a two-dimensional growth area, generated by gentle adhesion of a thin
PDMS membrane on a coverslip. We probed single-cell growth relative to drug gradients, which were generated under diverse patterns of perfusion. Overall, this work reports the first proof of concept of a single-cell dose-response assay for slow-growing bacterial cells, and demonstrates the ability of our platform to investigate PK-PD relationships at high spatial and temporal resolution. As opposed to bulk-cell methods that indirectly infer the presence of drug-persistent cells, our platform allows for direct detection of diverse single-cell fates over time, providing a technological advance in drug discovery.

Microfluidic devices used to trap, culture and image living cells are mostly made of micropatterned PDMS, which can be activated or functionalized on its surface, to generate irreversible chemical bonds with compatible materials, such as glass or other polymers⁴⁶. Typically, surface activation is obtained inside a vacuum chamber, by ionization of oxygen via an electric field that oscillates as a function of time. This so-called plasma treatment removes impurities and generates functional silanol groups, which make the surface hydrophilic, and covalently bind the same groups generated on other substrates⁶⁵. To seal our platforms, we use this plasma treatment followed by manual alignment, in the absence of pressure from the control layer, aiming to prevent sticking of the PDMS membrane to the underlying glass. However, the current manufacturing process of our 5-condition platform has a success rate around 80%, primarily due to manual alignment errors between the different layers that constitute the device, and to patchy adhesion of the PDMS membrane to the glass surface in some growth microchambers. These limitations in the manufacturing of the platform can impair the stability of the system, and the delivery of nutrients to cells as well as their survival. Many variables are involved in the process of surface oxidation of PDMS, thereby producing inter-device variation, which is hard to predict. By fine-tuning plasma dose, exposure time and pressure, the surface oxidation of PDMS can lead to local vitrification⁶⁶. Vitrification increases the stiffness of PDMS and leads to the formation of nanometer-sized wrinkles, which could not only prevent unintended bonding, but also improve the circulation of the nutrients within the 2D-growth area. Although PDMS vitrification holds promise for device optimization, evaluation of the conditions, functioning of the platform and cell viability remain to be tested. In order to make the manufacturing process increasingly reproducible, PDMS could also be blended with amphiphilic or hydrophilic additives⁶⁷. The resulting hybrid materials exhibit more-stable properties over time, and are less susceptible to non-specific adsorption of hydrophobic molecules, possibly improving the circulation of drug candidates. Overall, PDMS is an excellent material for prototyping, but it has a number of drawbacks, including adsorption, liquid evaporation, and low-throughput manufacturing and packaging processes³⁴.

Aiming to improve the reproducibility and throughput of the fabrication process, the use of alternative materials and more efficient replication techniques can be considered. In particular, transparent amorphous thermoplastics, such as polystyrene, polymethylmethacrylate, cyclic olefin copolymer, or sulfone polymers, can be thermoformed by hot embossing against a metal mold heated to high temperatures using a press; by injection molding inside a mold⁶⁸,⁶⁹; or by additive stereolithography⁷⁰. To replace the flexible membrane of our device, we could instead turn to transparent thermoplastic elastomers, copolymers of plastic and rubber that can be stretched without breaking and returning to their original position, such as styrenic block copolymers; Dryflex-T; or Hytrel. Irrespective of the manufacturing process, the modularity and inherent versatility of our platform make it suitable to structural modifications and upgrades. For instance, it could be scaled up by increasing the number of inlet ports, branches of the dilution tree and microchambers. However, any modification from the original design will require verifying device stability and controllability, gradient formation and homogeneous cell growth across microchambers from scratch.

Since our platform was designed to trap and monitor the behavior of very small bacterial cells, using a flexible PDMS membrane that exerts low pressure on cells, allows gases to permeate and nutrients to circulate, and limits shear stress, it will most likely be compatible with a broad spectrum of microbial cells, eukaryotic cells, organoids, and may be adaptable to more complex multicellular models⁷¹,⁷². Additionally, the current platform allows us to combine experimental conditions in a virtually unlimited manner. For instance, nutrients, cell stimuli, flow dynamics, atmosphere and pressure can be combined and modified as a function of time, aiming to simulate different environmental micro-niches. In particular, to study the effect of drugs it is important to reproduce fundamental PK parameters⁷². Although at present our system does not include components of the immune response, we can create diverse drug exposure patterns, and quantify the dynamics of response of individual cells under clinically relevant infusion patterns. To this aim, target-controlled infusion software systems could be implemented, in conjunction with multiple inlet ports or multi-directional valves, aiming to fine-tune each fluidic input according to in-silico models inferred from human PK parameters⁷³⁻⁷⁵. Spatiotemporal drug gradients could also be achieved by modulated oscillation of the inlet tubing of microfluidic devices. This would generate harmonic flow patterns that simulate blood and lymphatic vasculatures, which can be sites of dissemination and invasion for pathogenic bacteria⁷⁶. Drug stimulation could also be combined with customized thermal gradients, thanks to the implementation of convective heat exchangers that fine-tune the local temperature of the microfluidic device⁷⁷. Temperature fluctuations would make it possible to simulate febrile illness, for instance, and to investigate the bacterial response to combined stresses.

In-vitro models enabling accurate PK-PD analysis with high predictive potential for the in-vivo phase are essential but limited. Chief among them is the so-called hollow fiber, a well-known dynamic macrosystem used for preclinical assessment of drug activity and selection of dosing strategies, also for antitubercular drugs⁷⁸,⁷⁹. The hollow fiber is a cylindrical bioreactor controlled by peristaltic pumps, where the flow of nutrients and drugs reaches bacterial cells by diffusion through a series of semi-permeable capillary fibers. Growth inhibition and recovery during fluctuating drug exposure, which mimics typical human PK profiles, are monitored by regular sampling and assessment of colony-forming units⁷⁹⁻⁸¹. As a consequence, the hollow-fiber system is inherently limited by the ability of colonies to resume growth on plate, and cannot detect cells that are non-growing but metabolically active and pose a risk for relapse¹¹. Unlike the hollow fiber, our microsystem can determine PD parameters such as division and killing rates at the single-cell and microcolony level, and correlate them with
PDMS membrane on a coverslip. We probed single-cell growth relative to drug gradients, which were generated under diverse patterns of perfusion. Overall, this work reports the first proof of concept of a single-cell dose-response assay for slow-growing bacterial cells, and demonstrates the ability of our platform to investigate PK-PD relationships at high spatial and temporal resolution. As opposed to bulk-cell methods that indirectly infer the presence of drug-persistent cells, our platform allows for direct detection of diverse single-cell fates over time, providing a technological advance in drug discovery.

Microfluidic devices used to trap, culture and image living cells are mostly made of micropatterned PDMS, which can be activated or functionalized on its surface, to generate irreversible chemical bonds with compatible materials, such as glass or other polymers⁴⁶. Typically, surface activation is obtained inside a vacuum chamber, by ionization of oxygen via an electric field that oscillates as a function of time. This so-called plasma treatment removes impurities and generates functional silanol groups, which make the surface hydrophilic, and covalently bind the same groups generated on other substrates⁶⁵. To seal our platforms, we use this plasma treatment followed by manual alignment, in the absence of pressure from the control layer, aiming to prevent sticking of the PDMS membrane to the underlying glass. However, the current manufacturing process of our 5-condition platform has a success rate around 80%, primarily due to manual alignment errors between the different layers that constitute the device, and to patchy adhesion of the PDMS membrane to the glass surface in some growth microchambers. These limitations in the manufacturing of the platform can impair the stability of the system, and the delivery of nutrients to cells as well as their survival. Many variables are involved in the process of surface oxidation of PDMS, thereby producing inter-device variation, which is hard to predict. By fine-tuning plasma dose, exposure time and pressure, the surface oxidation of PDMS can lead to local vitrification⁶⁶. Vitrification increases the stiffness of PDMS and leads to the formation of nanometer-sized wrinkles, which could not only prevent unintended bonding, but also improve the circulation of the nutrients within the 2D-growth area. Although PDMS vitrification holds promise for device optimization, evaluation of the conditions, functioning of the platform and cell viability remain to be tested. In order to make the manufacturing process increasingly reproducible, PDMS could also be blended with amphiphilic or hydrophilic additives⁶⁷. The resulting hybrid materials exhibit more-stable properties over time, and are less susceptible to non-specific adsorption of hydrophobic molecules, possibly improving the circulation of drug candidates. Overall, PDMS is an excellent material for prototyping, but it has a number of drawbacks, including adsorption, liquid evaporation, and low-throughput manufacturing and packaging processes³⁴.

Aiming to improve the reproducibility and throughput of the fabrication process, the use of alternative materials and more efficient replication techniques can be considered. In particular, transparent amorphous thermoplastics, such as polystyrene, polymethylmethacrylate, cyclic olefin copolymer, or sulfone polymers, can be thermoformed by hot embossing against a metal mold heated to high temperatures using a press; by injection molding inside a mold⁶⁸,⁶⁹; or by additive stereolithography⁷⁰. To replace the flexible membrane of our device, we could instead turn to transparent thermoplastic elastomers, copolymers of plastic and rubber that can be stretched without breaking and returning to their original position, such as styrenic block copolymers; Dryflex-T; or Hytrel. Irrespective of the manufacturing process, the modularity and inherent versatility of our platform make it suitable to structural modifications and upgrades. For instance, it could be scaled up by increasing the number of inlet ports, branches of the dilution tree and microchambers. However, any modification from the original design will require verifying device stability and controllability, gradient formation and homogeneous cell growth across microchambers from scratch.

Since our platform was designed to trap and monitor the behavior of very small bacterial cells, using a flexible PDMS membrane that exerts low pressure on cells, allows gases to permeate and nutrients to circulate, and limits shear stress, it will most likely be compatible with a broad spectrum of microbial cells, eukaryotic cells, organoids, and may be adaptable to more complex multicellular models⁷¹,⁷². Additionally, the current platform allows us to combine experimental conditions in a virtually unlimited manner. For instance, nutrients, cell stimuli, flow dynamics, atmosphere and pressure can be combined and modified as a function of time, aiming to simulate different environmental micro-niches. In particular, to study the effect of drugs it is important to reproduce fundamental PK parameters⁷². Although at present our system does not include components of the immune response, we can create diverse drug exposure patterns, and quantify the dynamics of response of individual cells under clinically relevant infusion patterns. To this aim, target-controlled infusion software systems could be implemented, in conjunction with multiple inlet ports or multi-directional valves, aiming to fine-tune each fluidic input according to in-silico models inferred from human PK parameters⁷³⁻⁷⁵. Spatiotemporal drug gradients could also be achieved by modulated oscillation of the inlet tubing of microfluidic devices. This would generate harmonic flow patterns that simulate blood and lymphatic vasculatures, which can be sites of dissemination and invasion for pathogenic bacteria⁷⁶. Drug stimulation could also be combined with customized thermal gradients, thanks to the implementation of convective heat exchangers that fine-tune the local temperature of the microfluidic device⁷⁷. Temperature fluctuations would make it possible to simulate febrile illness, for instance, and to investigate the bacterial response to combined stresses.

In-vitro models enabling accurate PK-PD analysis with high predictive potential for the in-vivo phase are essential but limited. Chief among them is the so-called hollow fiber, a well-known dynamic macrosystem used for preclinical assessment of drug activity and selection of dosing strategies, also for antitubercular drugs⁷⁸,⁷⁹. The hollow fiber is a cylindrical bioreactor controlled by peristaltic pumps, where the flow of nutrients and drugs reaches bacterial cells by diffusion through a series of semi-permeable capillary fibers. Growth inhibition and recovery during fluctuating drug exposure, which mimics typical human PK profiles, are monitored by regular sampling and assessment of colony-forming units⁷⁹⁻⁸¹. As a consequence, the hollow-fiber system is inherently limited by the ability of colonies to resume growth on plate, and cannot detect cells that are non-growing but metabolically active and pose a risk for relapse¹¹. Unlike the hollow fiber, our microsystem can determine PD parameters such as division and killing rates at the single-cell and microcolony level, and correlate them with
critical PK parameters, also providing real-time insights into bactericidal versus bacteriostatic activity and risk of relapse for individual cells. Additionally, several drug doses can be tested simultaneously, surpassing the throughput of standard microfluidic systems<sup>46,47</sup>. Finally, using appropriate fluorescent reporters, we can also monitor the dynamics of induction and repression of the drug target as a function of drug dose and exposure time. The analysis presented here reveals not only heterogeneous responses to the same drug concentrations in individual cells, but also how small changes in drug concentration can impact the population dynamics, conceivably triggering survival mechanisms that might be conducive to drug resistance<sup>20-25</sup>.

Because effective antimicrobials are inexorably running out, posing one of the biggest challenges of the 21st century<sup>6</sup>, we urgently need new tools to help us better understand how drugs work and how best to use them. Our platform will prove useful as a dynamic single-cell PK-PD in-vitro model, to clarify the activity of antimicrobials and how certain bacterial cells escape treatment. Furthermore, the integration of single-cell PD datasets with clinical ones, through mathematical modeling and simulations<sup>29,82</sup>, might also help to better predict in-vivo efficacy. We ultimately hope that our work will contribute to the establishment of more robust drug-development pipelines, aiming to prevent disease relapse and to improve cure rates in patients.

## Methods

**Bacterial strains and growth conditions.** Plasmids cloning, sequencing and amplification steps were carried out in *Escherichia coli* TOP10 strains, grown at 37 °C under shaking conditions in LB medium containing the appropriate selective antibiotic: 100 µg/mL ampicillin; 50 µg/mL kanamycin; 100 µg/mL hygromycin. All mycobacterial strains were cultured at 37 °C under shaking conditions in Middlebrook 7H9 broth enriched with 0.5% BSA, 0.5% glycerol, 0.2% glucose, 0.085% NaCl, and 0.05% Tween-80. Middlebrook 7H10 agar plates were supplemented with 10% OADC and 0.5% glycerol. Mycobacterial transformants were selected on 20 µg/mL kanamycin, or a combination of 15 µg/mL kanamycin and 50 µg/mL hygromycin. For the construction of GyrA-mCherry fluorescent reporters, second recombination events were selected on Middlebrook 7H10 containing 5% sucrose for *M. smegmatis* or 2.5% sucrose for *M. tuberculosis*.

Bacterial stocks were prepared by adding 15% glycerol into late exponential-phase cultures (OD<sub>600nm</sub> 1.0) obtained from single colonies, frozen at -80 °C, and used only once to start primary cultures (1:100 dilution). For final assays, primary cultures were grown up to mid-log phase (OD<sub>600nm</sub> 0.5-0.8) and used to inoculate final secondary cultures (OD<sub>600nm</sub> 0.025 dilution) until mid-log phase. For strains carrying chromosomal integrative plasmids, the selective antibiotic was maintained in both primary and secondary cultures, but removed during microscopy experiments.

**Strains construction.** Translational fluorescent reporter strains were generated using a stable mCherry variant<sup>83</sup>. An oligopeptide linker was used to fuse GyrA and mCherry, so as not to affect the folding of the native protein. Plasmids used for cloning were checked by sequencing and final constructs were confirmed by restriction enzyme profiling. Stable GFP variant under the control of a constitutive promoter was inserted into the L5 phage integration site<sup>13</sup>. Mycobacteria were transformed by electroporation.

Fluorescent reporters of *M. smegmatis* *gyrA* (*MSMEG_0006*) and *M. tuberculosis* *gyrA* (*rv0006*) were obtained by two events of homologous recombination, using the suicide vector pJG1100<sup>13</sup>. In particular, the *mCherry* open reading frame was inserted into the *gyrA* chromosomal locus just upstream the stop codon of *gyrA* via an oligonucleotide sequence, generating C terminus translational reporters. Amplicons of about 600 bp, upstream and downstream the chromosomal insertion site, were ligated to mCherry using compatible restriction sites, and cloned into pJG1100. After electroporation, the first recombination event was selected on antibiotics, while the second recombination event was selected on sucrose, and both events were confirmed by PCR analysis of genomic DNA (Supplementary Figs. 5a, 6a). The final GyrA-mCherry reporters were named GMS5 for *M. smegmatis* and GMT18 for *M. tuberculosis*. GMT18 was further transformed with a plasmid constitutively expressing GFP, and the final double fluorescent reporter was named GMT18-GFP (Supplementary Table 1).

**MIC assessment by resazurin assay.** Secondary cultures were diluted to OD<sub>600nm</sub> 0.005 using pre-warmed Middlebrook 7H9 medium without detergent, and 100 µL of bacterial suspension were dispensed into a 96-well plate, except the well with the highest concentration of drug to be tested that was filled with 200 µL of cell suspension. The drug was added to this well and then two-fold serially diluted. Plates were incubated at 37 °C for 24 h for *M. smegmatis* and 7 days for *M. tuberculosis*, then 10 µL of 0.01% resazurin were added to each well and incubated for 24 h at 37 °C, before visual quantification. The lowest drug concentration causing cell mortality, assessed as blue color, accounted for the MIC.

**Microfabrication.** The fabrication of the microfluidic devices consisted of the following steps: mask design and printing; wafers fabrication; soft lithography; PDMS layers assembly and PDMS-glass bonding. Masks were designed using AutoCAD 2017, and low-resolution photomasks were printed at Selba S.A. (Supplementary Data 1).

For each device, two micropatterned silicon wafers were created, one for the flow layer (FL-W) and the second one for the control layer (CL-W). Before coating, standard silicon wafers were dehydrated at 200 °C on a hot plate for 30 min. The FL-W was produced by spin-coating SU8-2025 photoresist onto a wafer at 3000 rpm for 30 s, to generate a 30-µm thick layer, and soft-backed at 65 °C for 1 min, followed by incubation at 95 °C for 5.5 min. The photoresist was crosslinked under 155 mJ/cm² UV light for 11.5 s, and the wafer backed at 65 °C for 1 min, followed by incubation at 95 °C for 5 min. The wafer was developed by immersion into propylene glycol methyl ether acetate (PGMEA) for 4 min and washed with isopropanol. The CL-W was generated by spin-coating the SU8-2100 photoresist onto a standard silicon wafer at 1500 rpm for 30 s, to generate a 200-µm thick layer, and
critical PK parameters, also providing real-time insights into bactericidal versus bacteriostatic activity and risk of relapse for individual cells. Additionally, several drug doses can be tested simultaneously, surpassing the throughput of standard microfluidic systems<sup>46,47</sup>. Finally, using appropriate fluorescent reporters, we can also monitor the dynamics of induction and repression of the drug target as a function of drug dose and exposure time. The analysis presented here reveals not only heterogeneous responses to the same drug concentrations in individual cells, but also how small changes in drug concentration can impact the population dynamics, conceivably triggering survival mechanisms that might be conducive to drug resistance<sup>20-25</sup>.

Because effective antimicrobials are inexorably running out, posing one of the biggest challenges of the 21st century<sup>6</sup>, we urgently need new tools to help us better understand how drugs work and how best to use them. Our platform will prove useful as a dynamic single-cell PK-PD in-vitro model, to clarify the activity of antimicrobials and how certain bacterial cells escape treatment. Furthermore, the integration of single-cell PD datasets with clinical ones, through mathematical modeling and simulations<sup>29,82</sup>, might also help to better predict in-vivo efficacy. We ultimately hope that our work will contribute to the establishment of more robust drug-development pipelines, aiming to prevent disease relapse and to improve cure rates in patients.

## Methods

**Bacterial strains and growth conditions.** Plasmids cloning, sequencing and amplification steps were carried out in *Escherichia coli* TOP10 strains, grown at 37 °C under shaking conditions in LB medium containing the appropriate selective antibiotic: 100 µg/mL ampicillin; 50 µg/mL kanamycin; 100 µg/mL hygromycin. All mycobacterial strains were cultured at 37 °C under shaking conditions in Middlebrook 7H9 broth enriched with 0.5% BSA, 0.5% glycerol, 0.2% glucose, 0.085% NaCl, and 0.05% Tween-80. Middlebrook 7H10 agar plates were supplemented with 10% OADC and 0.5% glycerol. Mycobacterial transformants were selected on 20 µg/mL kanamycin, or a combination of 15 µg/mL kanamycin and 50 µg/mL hygromycin. For the construction of GyrA-mCherry fluorescent reporters, second recombination events were selected on Middlebrook 7H10 containing 5% sucrose for *M. smegmatis* or 2.5% sucrose for *M. tuberculosis*.

Bacterial stocks were prepared by adding 15% glycerol into late exponential-phase cultures (OD<sub>600nm</sub> 1.0) obtained from single colonies, frozen at -80 °C, and used only once to start primary cultures (1:100 dilution). For final assays, primary cultures were grown up to mid-log phase (OD<sub>600nm</sub> 0.5-0.8) and used to inoculate final secondary cultures (OD<sub>600nm</sub> 0.025 dilution) until mid-log phase. For strains carrying chromosomal integrative plasmids, the selective antibiotic was maintained in both primary and secondary cultures, but removed during microscopy experiments.

**Strains construction.** Translational fluorescent reporter strains were generated using a stable mCherry variant<sup>83</sup>. An oligopeptide linker was used to fuse GyrA and mCherry, so as not to affect the folding of the native protein. Plasmids used for cloning were checked by sequencing and final constructs were confirmed by restriction enzyme profiling. Stable GFP variant under the control of a constitutive promoter was inserted into the L5 phage integration site<sup>13</sup>. Mycobacteria were transformed by electroporation.

Fluorescent reporters of *M. smegmatis* *gyrA* (*MSMEG_0006*) and *M. tuberculosis* *gyrA* (*rv0006*) were obtained by two events of homologous recombination, using the suicide vector pJG1100<sup>13</sup>. In particular, the *mCherry* open reading frame was inserted into the *gyrA* chromosomal locus just upstream the stop codon of *gyrA* via an oligonucleotide sequence, generating C terminus translational reporters. Amplicons of about 600 bp, upstream and downstream the chromosomal insertion site, were ligated to mCherry using compatible restriction sites, and cloned into pJG1100. After electroporation, the first recombination event was selected on antibiotics, while the second recombination event was selected on sucrose, and both events were confirmed by PCR analysis of genomic DNA (Supplementary Figs. 5a, 6a). The final GyrA-mCherry reporters were named GMS5 for *M. smegmatis* and GMT18 for *M. tuberculosis*. GMT18 was further transformed with a plasmid constitutively expressing GFP, and the final double fluorescent reporter was named GMT18-GFP (Supplementary Table 1).

**MIC assessment by resazurin assay.** Secondary cultures were diluted to OD<sub>600nm</sub> 0.005 using pre-warmed Middlebrook 7H9 medium without detergent, and 100 µL of bacterial suspension were dispensed into a 96-well plate, except the well with the highest concentration of drug to be tested that was filled with 200 µL of cell suspension. The drug was added to this well and then two-fold serially diluted. Plates were incubated at 37 °C for 24 h for *M. smegmatis* and 7 days for *M. tuberculosis*, then 10 µL of 0.01% resazurin were added to each well and incubated for 24 h at 37 °C, before visual quantification. The lowest drug concentration causing cell mortality, assessed as blue color, accounted for the MIC.

**Microfabrication.** The fabrication of the microfluidic devices consisted of the following steps: mask design and printing; wafers fabrication; soft lithography; PDMS layers assembly and PDMS-glass bonding. Masks were designed using AutoCAD 2017, and low-resolution photomasks were printed at Selba S.A. (Supplementary Data 1).

For each device, two micropatterned silicon wafers were created, one for the flow layer (FL-W) and the second one for the control layer (CL-W). Before coating, standard silicon wafers were dehydrated at 200 °C on a hot plate for 30 min. The FL-W was produced by spin-coating SU8-2025 photoresist onto a wafer at 3000 rpm for 30 s, to generate a 30-µm thick layer, and soft-backed at 65 °C for 1 min, followed by incubation at 95 °C for 5.5 min. The photoresist was crosslinked under 155 mJ/cm² UV light for 11.5 s, and the wafer backed at 65 °C for 1 min, followed by incubation at 95 °C for 5 min. The wafer was developed by immersion into propylene glycol methyl ether acetate (PGMEA) for 4 min and washed with isopropanol. The CL-W was generated by spin-coating the SU8-2100 photoresist onto a standard silicon wafer at 1500 rpm for 30 s, to generate a 200-µm thick layer, and
soft-backed at 65 °C for 5 min, followed by incubation at 95 °C for 20 min. The photoresist was crosslinked under 240 mJ/cm² UV light for 21.8 s, and the wafer backed at 65 °C for 5 min, followed by incubation at 95 °C for 10 min. The wafer was developed by immersion into PGMEA for 15 min and washed with isopropanol. The FL-W and CL-W were hard-backed at 180 °C for 2 h, and silanized with a drop of Trichloro (1H,1H,2H,2H-perfluorooctyl) silane under a vacuum chamber overnight. The height of the microchannels was confirmed using a profilometer DektakXT (Bruker).

The FL-W and CL-W were used as molds for PDMS soft lithography. To promote bonding between the two layers, we used a different amount of cross-linking agent, which diffuses from the layer with the highest concentration to the layer with the lowest concentration. The FL-W was covered with a degassed mixture of 20 g of silicone elastomer SYLGARD 184 (Dow Corning) and 1 g of cross-linking agent (20:1 ratio), spin-coated at 1500 rpm for 60 s to obtain a 50-µm thick layer, incubated at room temperature for 20 min on a flat surface, and baked at 80 °C for 18 min. The CL-W was placed into a square Petri dish and covered with a degassed mixture of 40 g of silicone elastomer and 8 g of cross-linking agent (5:1 ratio), further degassed for 30 min in a vacuum chamber, baked at 60 °C for 30 min, and the patterned PDMS area was cut with a blade and gently separated from the mold.

The inlet ports were generated using a biopsy punch (1.0 mm OD, Harris Uni-Core). The PDMS CL was manually aligned with the FL, using the circular chambers in the two layers as references. Metallic connectors (0.8/1.2 mm ID/OD, Phymep) were inserted into the inlet port of the top element. Finally, the merged FL and CL were immersed in a degassed mixture of 30 g of silicone elastomer and 3 g of cross-linking agent (10:1 ratio), placed in a vacuum chamber for 30 min until degassing was completed, and baked at 80 °C overnight. The whole PDMS assembly was gently separated from the FL-W, the metallic connectors were unplugged from the inlet and outlet ports of the CL, and the inlet and outlet ports of the FL were generated with a biopsy punch (1.0 mm OD).

The PDMS assembly was bonded with a glass coverslip by oxygen plasma treatment. For the single-condition device we used a 24 × 50 × 0.175 mm coverslip, and for the 5-condition platform we used a 40 × 60 × 0.175 mm coverslip. Coverslips were cleaned with isopropanol prior to plasma treatment. Both the coverslip and the PDMS assembly, with the microchannels facing upward, were placed inside the plasma chamber. Initially the vacuum was created until 0.1 mbar, then oxygen was injected until a pressure of 0.25 mbar was reached, and the treatment was maintained for 30 s for the single-condition device, and for 60 s for the 5-condition device. Upon plasma activation, the PDMS assembly was placed on the coverslip and incubated at 80 °C overnight to complete the bonding. For decontamination, the device was exposed to ozone/UV treatment for 15 min inside a UVO-Cleaner (Jelight) and stored in a sterile 20-mm petri dish.

**Characterization of PDMS membrane deformation.** The deformation of the PDMS membrane was evaluated in the single-condition device, using a 100 µM solution of fluorescein (FITC) perfused into the FL at 50 mbar, and applying a pressure range from 55 to 80 mbar in the CL with incremental steps of 1 mbar. The FL inlet was connected to the FITC reservoir and the FL outlet to a waste bottle. The CL inlet was connected to a water reservoir, and the CL outlet to a short tubing, which was closed with a plier upon water injection. The pressures in the CL and FL were fine-tuned via the MFCS-EZ flow controller (Fluigent), which is driven by the MAESFLO software. FITC fluorescence (Ex 475/28, Em 525/48 nm), 10% T and 10 ms exposure time was acquired from each microchamber at 30 s intervals, using an inverted DeltaVision Elite Microscope (Leica) equipped with an APOCHROMAT 20x/1.0 air objective (ZEISS) and a high-speed sCMOS camera, 2560 × 2160 pixels, pixel size 6.5 × 6.5 µm, 15-bit, spectral range of 370–1100 nm. Image stacks were analysed on ImageJ 1.53a, tracing a line along the diameter of each microchamber and acquiring the fluorescence intensity over time with the Plot Profile function. The background mean fluorescence intensity was extracted from areas of the device without liquid circulation and subtracted to the microchamber measurements. Fluorescence decrease in the microchambers was considered to be proportional to increasing pressure in the CL.

**FITC gradients generation and quantification.** Three 125 mL bottles were prepared, two containing 7H9 medium and the third containing a 100 µM FITC solution in 7H9 medium, and were connected to a bidirectional M-switch valve (Fluigent) with 10 inlet and 1 outlet port. One medium and the FITC bottle were connected to two separate inlet ports of the valve, and the outlet port of the valve was connected to the left inlet of the device. The right inlet port of the device was connected to the third bottle of medium. The outlet port of the device was connected to a flowmeter (Fluigent) and then to an empty waste bottle. The switch between conditions was achieved via the ESS Control software (Fluigent). A Tygon tubing was used to connect the CL inlet to a water bottle, and the CL outlet was connected to a short tubing. The three bottles were driven by the first channel (C1) of the flow controller (Fig. 2a). The M-switch was set on medium perfusion mode. The water bottle was driven by the second channel (C2) and the waste by the third channel (C3). The pressure was set at 100 mbar on C1 and C2 and, once the CL and FL were filled with liquids, their outlet tubing was closed with pliers and the pressure was increased to 400 mbar on C1. After 20 min of incubation to remove bubbles, the plier closing the FL outlet was removed. The pressure was set at 100 mbar in C1, 30 mbar in C2, and 20 mbar in C3. The MFCS-EZ flow controller (Fluigent) has a pressure sensor resolution of 0.03% full scale; a pressure stability < 0.1% CV; a response time down to 40 ms; and a settling time down to 100 ms. C1 (FL inlet) and C3 (FL outlet) range from 3 to 1000 mbar, and C2 (CL inlet) ranges from 1.035 to 345 mbar. Based on the measured flow rate in the FL, the pressure was fine-tuned in C1 to reach a flow rate of 150 µL/h. To verify the formation of different gradient patterns, FITC fluorescence was acquired from each microchamber over time. To generate static gradients (Fig. 3b), 7H9 medium and FITC solution were injected alternately for targeted periods of time by switching the position on the M-switch valve, in the absence of capacitors. To generate pulsing (PK-like) gradients (Fig. 3c), two capacitors were integrated just upstream the inlet ports of the FL. Capacitors were 1-cm Tygon tubing (2 mm ID,
soft-backed at 65 °C for 5 min, followed by incubation at 95 °C for 20 min. The photoresist was crosslinked under 240 mJ/cm² UV light for 21.8 s, and the wafer backed at 65 °C for 5 min, followed by incubation at 95 °C for 10 min. The wafer was developed by immersion into PGMEA for 15 min and washed with isopropanol. The FL-W and CL-W were hard-backed at 180 °C for 2 h, and silanized with a drop of Trichloro (1H,1H,2H,2H-perfluorooctyl) silane under a vacuum chamber overnight. The height of the microchannels was confirmed using a profilometer DektakXT (Bruker).

The FL-W and CL-W were used as molds for PDMS soft lithography. To promote bonding between the two layers, we used a different amount of cross-linking agent, which diffuses from the layer with the highest concentration to the layer with the lowest concentration. The FL-W was covered with a degassed mixture of 20 g of silicone elastomer SYLGARD 184 (Dow Corning) and 1 g of cross-linking agent (20:1 ratio), spin-coated at 1500 rpm for 60 s to obtain a 50-µm thick layer, incubated at room temperature for 20 min on a flat surface, and baked at 80 °C for 18 min. The CL-W was placed into a square Petri dish and covered with a degassed mixture of 40 g of silicone elastomer and 8 g of cross-linking agent (5:1 ratio), further degassed for 30 min in a vacuum chamber, baked at 60 °C for 30 min, and the patterned PDMS area was cut with a blade and gently separated from the mold.

The inlet ports were generated using a biopsy punch (1.0 mm OD, Harris Uni-Core). The PDMS CL was manually aligned with the FL, using the circular chambers in the two layers as references. Metallic connectors (0.8/1.2 mm ID/OD, Phymep) were inserted into the inlet port of the top element. Finally, the merged FL and CL were immersed in a degassed mixture of 30 g of silicone elastomer and 3 g of cross-linking agent (10:1 ratio), placed in a vacuum chamber for 30 min until degassing was completed, and baked at 80 °C overnight. The whole PDMS assembly was gently separated from the FL-W, the metallic connectors were unplugged from the inlet and outlet ports of the CL, and the inlet and outlet ports of the FL were generated with a biopsy punch (1.0 mm OD).

The PDMS assembly was bonded with a glass coverslip by oxygen plasma treatment. For the single-condition device we used a 24 × 50 × 0.175 mm coverslip, and for the 5-condition platform we used a 40 × 60 × 0.175 mm coverslip. Coverslips were cleaned with isopropanol prior to plasma treatment. Both the coverslip and the PDMS assembly, with the microchannels facing upward, were placed inside the plasma chamber. Initially the vacuum was created until 0.1 mbar, then oxygen was injected until a pressure of 0.25 mbar was reached, and the treatment was maintained for 30 s for the single-condition device, and for 60 s for the 5-condition device. Upon plasma activation, the PDMS assembly was placed on the coverslip and incubated at 80 °C overnight to complete the bonding. For decontamination, the device was exposed to ozone/UV treatment for 15 min inside a UVO-Cleaner (Jelight) and stored in a sterile 20-mm petri dish.

**Characterization of PDMS membrane deformation.** The deformation of the PDMS membrane was evaluated in the single-condition device, using a 100 µM solution of fluorescein (FITC) perfused into the FL at 50 mbar, and applying a pressure range from 55 to 80 mbar in the CL with incremental steps of 1 mbar. The FL inlet was connected to the FITC reservoir and the FL outlet to a waste bottle. The CL inlet was connected to a water reservoir, and the CL outlet to a short tubing, which was closed with a plier upon water injection. The pressures in the CL and FL were fine-tuned via the MFCS-EZ flow controller (Fluigent), which is driven by the MAESFLO software. FITC fluorescence (Ex 475/28, Em 525/48 nm), 10% T and 10 ms exposure time was acquired from each microchamber at 30 s intervals, using an inverted DeltaVision Elite Microscope (Leica) equipped with an APOCHROMAT 20x/1.0 air objective (ZEISS) and a high-speed sCMOS camera, 2560 × 2160 pixels, pixel size 6.5 × 6.5 µm, 15-bit, spectral range of 370–1100 nm. Image stacks were analysed on ImageJ 1.53a, tracing a line along the diameter of each microchamber and acquiring the fluorescence intensity over time with the Plot Profile function. The background mean fluorescence intensity was extracted from areas of the device without liquid circulation and subtracted to the microchamber measurements. Fluorescence decrease in the microchambers was considered to be proportional to increasing pressure in the CL.

**FITC gradients generation and quantification.** Three 125 mL bottles were prepared, two containing 7H9 medium and the third containing a 100 µM FITC solution in 7H9 medium, and were connected to a bidirectional M-switch valve (Fluigent) with 10 inlet and 1 outlet port. One medium and the FITC bottle were connected to two separate inlet ports of the valve, and the outlet port of the valve was connected to the left inlet of the device. The right inlet port of the device was connected to the third bottle of medium. The outlet port of the device was connected to a flowmeter (Fluigent) and then to an empty waste bottle. The switch between conditions was achieved via the ESS Control software (Fluigent). A Tygon tubing was used to connect the CL inlet to a water bottle, and the CL outlet was connected to a short tubing. The three bottles were driven by the first channel (C1) of the flow controller (Fig. 2a). The M-switch was set on medium perfusion mode. The water bottle was driven by the second channel (C2) and the waste by the third channel (C3). The pressure was set at 100 mbar on C1 and C2 and, once the CL and FL were filled with liquids, their outlet tubing was closed with pliers and the pressure was increased to 400 mbar on C1. After 20 min of incubation to remove bubbles, the plier closing the FL outlet was removed. The pressure was set at 100 mbar in C1, 30 mbar in C2, and 20 mbar in C3. The MFCS-EZ flow controller (Fluigent) has a pressure sensor resolution of 0.03% full scale; a pressure stability < 0.1% CV; a response time down to 40 ms; and a settling time down to 100 ms. C1 (FL inlet) and C3 (FL outlet) range from 3 to 1000 mbar, and C2 (CL inlet) ranges from 1.035 to 345 mbar. Based on the measured flow rate in the FL, the pressure was fine-tuned in C1 to reach a flow rate of 150 µL/h. To verify the formation of different gradient patterns, FITC fluorescence was acquired from each microchamber over time. To generate static gradients (Fig. 3b), 7H9 medium and FITC solution were injected alternately for targeted periods of time by switching the position on the M-switch valve, in the absence of capacitors. To generate pulsing (PK-like) gradients (Fig. 3c), two capacitors were integrated just upstream the inlet ports of the FL. Capacitors were 1-cm Tygon tubing (2 mm ID,
BioRad), having Luer connectors on each extremity. The pressure in the FL was fine-tuned to 120 µL/h on 7H9 medium mode. To initiate a single pulse, the M-switch valve was moved to the FITC position for 12 min and switched back to 7H9 medium position for 12 h. To generate sequential pulses, the M-switch valve was moved to FITC position for 12 min, switched back to 7H9 medium position for 6 h, and then the sequence was repeated. Fluorescence imaging was carried out as described above, and pixel intensity was measured on regions of interest (ROI) of equal size and subtracting the background intensity from areas of the device without liquid circulation.

**Device assembly, cell loading and fluids management.** To make the PDMS-glass devices suitable for cell loading and tracking by time lapse microscopy, different tubing, interfaces and fittings were prepared and autoclaved prior to final assembly (Supplementary Fig. 2).

Inlet ports of the FL were used for medium and drug perfusion from the reservoir bottles, which were closed with pressure lids (Fluigent), and the flangeless fittings (IDEX Health & Science) of these lids were connected to 50-cm Teflon tubing. In turn, the inlet ports were connected to 5-cm Tygon tubing (Saint-Gobain), connected to a metallic fitting at the other end. Then Teflon (TFE) and Tygon tubing were joined with a piece of silastic and sealed with PDMS, in the absence or presence of a capacitor. The outlet port of the FL, used both to seed the cell suspension and to drain waste fluid into the waste bottle, was connected to a tripartite tubing as follows. A PVDF Y-connector (Buerkle) was connected to two 10-cm Tygon tubing and one 50-cm Tygon tubing (1.2 mm ID, Saint-Gobain). The Tygon tubing used for cell loading was connected to a female Luer-to-barb connector and to a Luer lock (Cole Parmer). The other Tygon tubing was plugged into the device via a metallic fitting. The longer tubing was connected to a 250 mL waste bottle closed with a pressure lid (Fluigent) via flangeless fittings. A flowmeter (Fluigent) was also integrated between the outlet tubing and the waste bottle, to maintain the flow rate constant.

The inlet port of the CL was connected to a 50-cm Tygon tubing, having a metallic fitting at one end and a flangeless fitting at the other end, plugged into the water reservoir. The outlet port of the CL was connected to a 5-cm Tygon via metallic fitting and closed with a plier.

To load bacterial cells, the outlet tubing was closed with a plier. Exponentially growing cells were disaggregated through a pre-equilibrated 5-µm filter (Merck), and about 400 µL of cell suspension were gently injected into the ad hoc Tygon tubing by hand, using a 1-mL syringe. When loading pathogenic *M. tuberculosis* cells, 0.2 µm filters were also added at the extremities of the inlet tubing, to prevent the pathogen from coming out of the device during the loading procedure. Next, the device was incubated at room temperature for 10 min. The entire procedure was performed under a class II biosafety cabinet, and the waste bottle was filled up for one fifth with a 2.6% sodium hypochlorite solution. At the end of the incubation period the syringe was unplugged, replaced by the Luer lock, and the junction was sprayed with mycobactericidal disinfectant. Finally, the device was connected to the reservoir and waste bottles and mounted on the microscope stage (Fig. 2a; Supplementary Fig. 2).

To generate the drug gradient, one sterile reservoir bottle was filled with 7H9 medium, the second one with 7H9 medium supplemented with 250 ng/mL of moxifloxacin (5X-MIC). These bottles were connected to the first and second port of an M-switch valve, in turn connected to the inlet port 'c' of the device. The third reservoir bottle was filled with 7H9 medium, and connected to the inlet port 'd' of the device (Fig. 2a). Fluids were controlled from the multi-channel ESS Control software via pressure caps: C1 drove the device inlet of the FL, C2 drove the CL, and C3 drove the device outlet of the FL. After cell seeding, the CL was filled with water at 100 mbar and closed with a plier. The FL was perfused with 7H9 medium at 200 mbar for 20 min and the outlet closed with a plier for other 20 min to degas the pipe. Next, in the single condition device, C1 was set at 60 mbar and C2 was set at 30 mbar, targeting a final flowrate of about 100 µL/h. In the 5-condition platform, C1 was set at 60 mbar and both C2 and C3 were set at 30 mbar, targeting a final flowrate of about 150 µL/h. To inject the drug, the M-switch valve was moved to the second position, and the pressure was set at 60 mbar on C1, 30 mbar on C2, and 20 mbar on C3. At the end of drug injection, the M-switch valve was moved to the initial position.

**Spectrofluorimetric assessment of moxifloxacin adsorption by PDMS.** To measure the concentration of moxifloxacin in 7H9 medium before and after flowing through the microfluidic device, we used a spectrofluorimetric method, which was previously reported to assess the concentration of moxifloxacin in various biological fluids⁸⁴. 7H9 medium containing moxifloxacin at a concentration of 500 ng/mL was used both for the calibration curves and to perfuse inside the device after a 1:2 dilution in 7H9. The device was loaded from one inlet port via C1 of the flow controller, while the second inlet port was closed. The outflow of medium was regulated via a flowmeter, from which the outlet tubing was forked to either a waste container or a collection microtube, in turn connected to C2 of the flow controller. The medium was perfused at a rate of 200 µL/h. The input samples were collected directly from the external reservoir, while to collect the output samples, the flow was redirected from the waste container to the collection microtube using a binder clip. To flush out the drug, the external reservoir was replaced with a clean one containing only 7H9 medium. All samples were stored at 4 °C until the end of collection. Both calibration curves and samples reading were performed at the same time, using 100 µL of sample mixed with 100 µL of phosphoric acid-sodium phosphate buffer pH 8.2 (0.1 M). Fluorescence was quantified in 96-well plates with black bottom and transparent lid, using a plate reader (Ex 287, Em 465 nm), measurement time 100 ms and excitation bandwidth 12 nm.

**Time-lapse imaging.** The assembled device was mounted on the stage of our inverted DeltaVision Elite Microscope, equipped with a UPLFLN100XO2/PH3/1.30 oil objective (Olympus), using a custom-made metal support (Supplementary Fig. 3). Exposure conditions for *M. smegmatis*: phase contrast 100% T, 150 ms; mCherry (Ex 575/25, Em 625/45 nm) 50% T, 150 ms. Exposure conditions for *M. tuberculosis*: phase contrast
BioRad), having Luer connectors on each extremity. The pressure in the FL was fine-tuned to 120 µL/h on 7H9 medium mode. To initiate a single pulse, the M-switch valve was moved to the FITC position for 12 min and switched back to 7H9 medium position for 12 h. To generate sequential pulses, the M-switch valve was moved to FITC position for 12 min, switched back to 7H9 medium position for 6 h, and then the sequence was repeated. Fluorescence imaging was carried out as described above, and pixel intensity was measured on regions of interest (ROI) of equal size and subtracting the background intensity from areas of the device without liquid circulation.

**Device assembly, cell loading and fluids management.** To make the PDMS-glass devices suitable for cell loading and tracking by time lapse microscopy, different tubing, interfaces and fittings were prepared and autoclaved prior to final assembly (Supplementary Fig. 2).

Inlet ports of the FL were used for medium and drug perfusion from the reservoir bottles, which were closed with pressure lids (Fluigent), and the flangeless fittings (IDEX Health & Science) of these lids were connected to 50-cm Teflon tubing. In turn, the inlet ports were connected to 5-cm Tygon tubing (Saint-Gobain), connected to a metallic fitting at the other end. Then Teflon (TFE) and Tygon tubing were joined with a piece of silastic and sealed with PDMS, in the absence or presence of a capacitor. The outlet port of the FL, used both to seed the cell suspension and to drain waste fluid into the waste bottle, was connected to a tripartite tubing as follows. A PVDF Y-connector (Buerkle) was connected to two 10-cm Tygon tubing and one 50-cm Tygon tubing (1.2 mm ID, Saint-Gobain). The Tygon tubing used for cell loading was connected to a female Luer-to-barb connector and to a Luer lock (Cole Parmer). The other Tygon tubing was plugged into the device via a metallic fitting. The longer tubing was connected to a 250 mL waste bottle closed with a pressure lid (Fluigent) via flangeless fittings. A flowmeter (Fluigent) was also integrated between the outlet tubing and the waste bottle, to maintain the flow rate constant.

The inlet port of the CL was connected to a 50-cm Tygon tubing, having a metallic fitting at one end and a flangeless fitting at the other end, plugged into the water reservoir. The outlet port of the CL was connected to a 5-cm Tygon via metallic fitting and closed with a plier.

To load bacterial cells, the outlet tubing was closed with a plier. Exponentially growing cells were disaggregated through a pre-equilibrated 5-µm filter (Merck), and about 400 µL of cell suspension were gently injected into the ad hoc Tygon tubing by hand, using a 1-mL syringe. When loading pathogenic *M. tuberculosis* cells, 0.2 µm filters were also added at the extremities of the inlet tubing, to prevent the pathogen from coming out of the device during the loading procedure. Next, the device was incubated at room temperature for 10 min. The entire procedure was performed under a class II biosafety cabinet, and the waste bottle was filled up for one fifth with a 2.6% sodium hypochlorite solution. At the end of the incubation period the syringe was unplugged, replaced by the Luer lock, and the junction was sprayed with mycobactericidal disinfectant. Finally, the device was connected to the reservoir and waste bottles and mounted on the microscope stage (Fig. 2a; Supplementary Fig. 2).

To generate the drug gradient, one sterile reservoir bottle was filled with 7H9 medium, the second one with 7H9 medium supplemented with 250 ng/mL of moxifloxacin (5X-MIC). These bottles were connected to the first and second port of an M-switch valve, in turn connected to the inlet port 'c' of the device. The third reservoir bottle was filled with 7H9 medium, and connected to the inlet port 'd' of the device (Fig. 2a). Fluids were controlled from the multi-channel ESS Control software via pressure caps: C1 drove the device inlet of the FL, C2 drove the CL, and C3 drove the device outlet of the FL. After cell seeding, the CL was filled with water at 100 mbar and closed with a plier. The FL was perfused with 7H9 medium at 200 mbar for 20 min and the outlet closed with a plier for other 20 min to degas the pipe. Next, in the single condition device, C1 was set at 60 mbar and C2 was set at 30 mbar, targeting a final flowrate of about 100 µL/h. In the 5-condition platform, C1 was set at 60 mbar and both C2 and C3 were set at 30 mbar, targeting a final flowrate of about 150 µL/h. To inject the drug, the M-switch valve was moved to the second position, and the pressure was set at 60 mbar on C1, 30 mbar on C2, and 20 mbar on C3. At the end of drug injection, the M-switch valve was moved to the initial position.

**Spectrofluorimetric assessment of moxifloxacin adsorption by PDMS.** To measure the concentration of moxifloxacin in 7H9 medium before and after flowing through the microfluidic device, we used a spectrofluorimetric method, which was previously reported to assess the concentration of moxifloxacin in various biological fluids⁸⁴. 7H9 medium containing moxifloxacin at a concentration of 500 ng/mL was used both for the calibration curves and to perfuse inside the device after a 1:2 dilution in 7H9. The device was loaded from one inlet port via C1 of the flow controller, while the second inlet port was closed. The outflow of medium was regulated via a flowmeter, from which the outlet tubing was forked to either a waste container or a collection microtube, in turn connected to C2 of the flow controller. The medium was perfused at a rate of 200 µL/h. The input samples were collected directly from the external reservoir, while to collect the output samples, the flow was redirected from the waste container to the collection microtube using a binder clip. To flush out the drug, the external reservoir was replaced with a clean one containing only 7H9 medium. All samples were stored at 4 °C until the end of collection. Both calibration curves and samples reading were performed at the same time, using 100 µL of sample mixed with 100 µL of phosphoric acid-sodium phosphate buffer pH 8.2 (0.1 M). Fluorescence was quantified in 96-well plates with black bottom and transparent lid, using a plate reader (Ex 287, Em 465 nm), measurement time 100 ms and excitation bandwidth 12 nm.

**Time-lapse imaging.** The assembled device was mounted on the stage of our inverted DeltaVision Elite Microscope, equipped with a UPLFLN100XO2/PH3/1.30 oil objective (Olympus), using a custom-made metal support (Supplementary Fig. 3). Exposure conditions for *M. smegmatis*: phase contrast 100% T, 150 ms; mCherry (Ex 575/25, Em 625/45 nm) 50% T, 150 ms. Exposure conditions for *M. tuberculosis*: phase contrast
100% T, 150 ms; FITC (Ex 475/28, Em 525/48 nm) 32% T, 10 ms; mCherry (Ex 575/25, Em 625/45 nm) 50% T, 250 ms. In each microchamber between 5 and 30 x, y, z coordinates were recorded, each one containing from 3 to 10 bacilli. The image acquisition rate was set at intervals of 20 or 30 min for *M. smegmatis* and of 3 h for *M. tuberculosis*. At the end of the experiment, the entire platform was washed with 2.6% sodium hypochlorite, the device was separated from tubing and discarded into a sharp waste container. The electronic M-switch valve and the flowmeter were extensively washed with 2.6% sodium hypochlorite, 70% ethanol and sterile water for 1 h. Tubing, fluidic interfaces and fittings were placed into plastic bags, sprayed with mycobactericidal disinfectant and autoclaved. Each experiment was performed at least twice independently.

**Microcolony analysis.** The growth of individual microcolonies from the single condition device was assessed by manual segmentation over time. The analysis of microcolonies from the 5-condition device was carried out in a semi-automated manner (Supplementary Fig. 5b; Supplementary Data 2–7). Image stacks were separated by channel using a customized macro (Ch-Extraction.ijm). To generate binary masks, a specific fluorescence channel was used: mCherry for *M. smegmatis* (mCherry-Mask.ijm), and FITC for *M. tuberculosis* (GFP-Mask.ijm). Images were first smoothed and then the Auto threshold plugin was used to distinguish the foreground from the background pixels. The Huang method<sup>85</sup> was used for the mCherry channel, and the Intermodes method<sup>86</sup> was used for the FITC channel. The binary masks were used to measure area and mean fluorescence from the ROI (ROI-Measure.ijm), and the inverse of the binary masks were used to measure the background fluorescence (BG-Measure.ijm), to be subtracted from the corresponding ROI. Data were saved as .csv files, which were merged and processed, using a customized R script (csv-Management.R). The microcolony coordinates were also recorded from the image metadata files, in order to determine the position relative to the center of the microchamber.

**Single-cell analysis.** *Mycobacterium tuberculosis* single-cell area and mean fluorescence were measured from independent ROI using the Selection Brush Tool and the ROI Manager macro at given time points. In each image, the background mean fluorescence was also measured and subtracted from the fluorescence of individual cells, which were segmented in the corresponding image.

*Mycobacterium tuberculosis* single-cell counting over time was performed using the Cell Counter plugin, and scoring the total cell count 24 h before drug exposure; dividing cells (alive); lysed cells (dead); GFP-dim cells (likely dead); cells that resumed growth after drug exposure (alive); and non-growing GFP-positive cells (likely viable). Quantified subpopulations were expressed as a function of time and as fractions.

**Calculation of parameters.** The growth rate of individual microcolonies was measured by nonlinear regression analysis of the area increase as function of time, fitting an exponential growth equation: $Y = Y_0 \cdot \exp(k \cdot t)$, where $Y_0$ is the initial microcolony size (µm²) and $k$ is the rate constant (time⁻¹). The doubling time is $\ln2/k$.

The AUC and associated parameters were measured based on the highest concentration of moxifloxacin injected, according to the dilution factors assessed from the formation of FITC gradients, and using the MIC concentration (50 ng/mL) as the baseline.

Dividing and lysed cells were counted starting from the number of cells 24 h before drug exposure, and progressively adding the number of divisions and subtracting the number of lysis events, obtaining the total count throughout the experiment. Next, the number of division events was divided by the cell count at each time point and by the image acquisition time frequency, obtaining the division rate per hour.

The relative IC50 was computed from the fraction of viable cells at different moxifloxacin doses, by interpolating a four-parameter dose-response curve according to the variable-slope model. We considered viable those cells that were non-growing but GFP positive and those that resumed growth.

**Statistical analysis.** Plots and statistical analyses were generated using Prism 9 (GraphPad). One-way ANOVA and multiple comparisons tests were used to assess significant differences within datasets containing more than two groups. Linear relationship between *x*, *y* datasets was determined by computing the Pearson correlation coefficient. Plots include datasets derived from at least two independent experimental replicates. *N* values, *P* values and specific statistical tests are indicated in each legend.

**Data availability**

Requests for data, materials and methods should be addressed to and will be fulfilled by the corresponding author. Plasmids and bacterial strains generated in this study are available from the corresponding author with a completed Materials Transfer Agreement.

Received: 31 May 2022; Accepted: 11 November 2022

Published online: 15 November 2022

**References**

1. Grant, S. S. & Hung, D. T. Persistent bacterial infections, antibiotic tolerance, and the oxidative stress response. *Virulence* **4**, 273–283 (2014).

2. Bakkeren, E. et al. *Salmonella* persisters promote the spread of antibiotic resistance plasmids in the gut. *Nature* **573**, 276–280 (2019).

3. Van den Bergh, B., Fauvart, M. & Michiels, J. Formation, physiology, ecology, evolution and clinical importance of bacterial persisters. *FEMS Microbiol. Rev.* **41**, 219–251 (2017).

4. Balaban, N. Q. et al. Definitions and guidelines for research on antibiotic persistence. *Nat. Rev. Microbiol.* **17**, 441–448 (2019).

5. Kaldalu, N. et al. In vitro studies of persister cells. *Microbiol. Mol. Biol. Rev.* **84**, e00070 (2020).
100% T, 150 ms; FITC (Ex 475/28, Em 525/48 nm) 32% T, 10 ms; mCherry (Ex 575/25, Em 625/45 nm) 50% T, 250 ms. In each microchamber between 5 and 30 x, y, z coordinates were recorded, each one containing from 3 to 10 bacilli. The image acquisition rate was set at intervals of 20 or 30 min for *M. smegmatis* and of 3 h for *M. tuberculosis*. At the end of the experiment, the entire platform was washed with 2.6% sodium hypochlorite, the device was separated from tubing and discarded into a sharp waste container. The electronic M-switch valve and the flowmeter were extensively washed with 2.6% sodium hypochlorite, 70% ethanol and sterile water for 1 h. Tubing, fluidic interfaces and fittings were placed into plastic bags, sprayed with mycobactericidal disinfectant and autoclaved. Each experiment was performed at least twice independently.

**Microcolony analysis.** The growth of individual microcolonies from the single condition device was assessed by manual segmentation over time. The analysis of microcolonies from the 5-condition device was carried out in a semi-automated manner (Supplementary Fig. 5b; Supplementary Data 2–7). Image stacks were separated by channel using a customized macro (Ch-Extraction.ijm). To generate binary masks, a specific fluorescence channel was used: mCherry for *M. smegmatis* (mCherry-Mask.ijm), and FITC for *M. tuberculosis* (GFP-Mask.ijm). Images were first smoothed and then the Auto threshold plugin was used to distinguish the foreground from the background pixels. The Huang method<sup>85</sup> was used for the mCherry channel, and the Intermodes method<sup>86</sup> was used for the FITC channel. The binary masks were used to measure area and mean fluorescence from the ROI (ROI-Measure.ijm), and the inverse of the binary masks were used to measure the background fluorescence (BG-Measure.ijm), to be subtracted from the corresponding ROI. Data were saved as .csv files, which were merged and processed, using a customized R script (csv-Management.R). The microcolony coordinates were also recorded from the image metadata files, in order to determine the position relative to the center of the microchamber.

**Single-cell analysis.** *Mycobacterium tuberculosis* single-cell area and mean fluorescence were measured from independent ROI using the Selection Brush Tool and the ROI Manager macro at given time points. In each image, the background mean fluorescence was also measured and subtracted from the fluorescence of individual cells, which were segmented in the corresponding image.

*Mycobacterium tuberculosis* single-cell counting over time was performed using the Cell Counter plugin, and scoring the total cell count 24 h before drug exposure; dividing cells (alive); lysed cells (dead); GFP-dim cells (likely dead); cells that resumed growth after drug exposure (alive); and non-growing GFP-positive cells (likely viable). Quantified subpopulations were expressed as a function of time and as fractions.

**Calculation of parameters.** The growth rate of individual microcolonies was measured by nonlinear regression analysis of the area increase as function of time, fitting an exponential growth equation: $Y = Y_0 \* \exp(k \* t)$, where $Y_0$ is the initial microcolony size (µm²) and $k$ is the rate constant (time⁻¹). The doubling time is $\ln2/k$.

The AUC and associated parameters were measured based on the highest concentration of moxifloxacin injected, according to the dilution factors assessed from the formation of FITC gradients, and using the MIC concentration (50 ng/mL) as the baseline.

Dividing and lysed cells were counted starting from the number of cells 24 h before drug exposure, and progressively adding the number of divisions and subtracting the number of lysis events, obtaining the total count throughout the experiment. Next, the number of division events was divided by the cell count at each time point and by the image acquisition time frequency, obtaining the division rate per hour.

The relative IC50 was computed from the fraction of viable cells at different moxifloxacin doses, by interpolating a four-parameter dose-response curve according to the variable-slope model. We considered viable those cells that were non-growing but GFP positive and those that resumed growth.

**Statistical analysis.** Plots and statistical analyses were generated using Prism 9 (GraphPad). One-way ANOVA and multiple comparisons tests were used to assess significant differences within datasets containing more than two groups. Linear relationship between *x*, *y* datasets was determined by computing the Pearson correlation coefficient. Plots include datasets derived from at least two independent experimental replicates. *N* values, *P* values and specific statistical tests are indicated in each legend.

**Data availability**

Requests for data, materials and methods should be addressed to and will be fulfilled by the corresponding author. Plasmids and bacterial strains generated in this study are available from the corresponding author with a completed Materials Transfer Agreement.

Received: 31 May 2022; Accepted: 11 November 2022

Published online: 15 November 2022

**References**

1. Grant, S. S. & Hung, D. T. Persistent bacterial infections, antibiotic tolerance, and the oxidative stress response. *Virulence* **4**, 273–283 (2014).

2. Bakkeren, E. et al. *Salmonella* persisters promote the spread of antibiotic resistance plasmids in the gut. *Nature* **573**, 276–280 (2019).

3. Van den Bergh, B., Fauvart, M. & Michiels, J. Formation, physiology, ecology, evolution and clinical importance of bacterial persisters. *FEMS Microbiol. Rev.* **41**, 219–251 (2017).

4. Balaban, N. Q. et al. Definitions and guidelines for research on antibiotic persistence. *Nat. Rev. Microbiol.* **17**, 441–448 (2019).

5. Kaldalu, N. et al. In vitro studies of persister cells. *Microbiol. Mol. Biol. Rev.* **84**, e00070 (2020).
6. Tacconelli, E. et al. Discovery, research, and development of new antibiotics: The WHO priority list of antibiotic-resistant bacteria and tuberculosis. *Lancet Infect. Dis.* **18**, 318–327 (2018).

7. Sarathy, J. P. et al. Extreme drug tolerance of *Mycobacterium tuberculosis* in caseum. *Antimicrob. Agents Chemother.* **62**, 25 (2018).

8. Ordonez, A. A. et al. Dynamic imaging in patients with tuberculosis reveals heterogeneous drug exposures in pulmonary lesions. *Nat. Med.* **26**, 529–534 (2020).

9. Cicchese, J. M., Dartois, V., Kirschner, D. E. & Linderman, J. J. Both pharmacokinetic variability and granuloma heterogeneity impact the ability of the first-line antibiotics to sterilize tuberculosis granulomas. *Front. Pharmacol.* **11**, 25 (2020).

10. Jo, J., Price-Whelan, A. & Dietrich, L. E. P. Gradients and consequences of heterogeneity in biofilms. *Nat. Rev. Microbiol.* **25**, 25 (2022).

11. Dhar, N., McKinney, J. & Manina, G. Phenotypic heterogeneity in *Mycobacterium tuberculosis*. *Microbiol. Spectr.* **4**, TBTB2-0021–2016 (2016).

12. Larrouy-Maumus, G. et al. Cell-envelope remodeling as a determinant of phenotypic antibacterial tolerance in *Mycobacterium tuberculosis*. *ACS Infect. Dis.* **2**, 352–360 (2016).

13. Manina, G., Griego, A., Singh, L. K., McKinney, J. D. & Dhar, N. Preexisting variation in DNA damage response predicts the fate of single mycobacteria under stress. *EMBO J.* **38**, e101876 (2019).

14. Levien, E., Kondev, J. & Amir, A. The interplay of phenotypic variability and fitness in finite microbial populations. *J. R. Soc. Interface* **17**, 20190827 (2020).

15. Mishra, R., Yadav, V., Guha, M. & Singh, A. Heterogeneous host-pathogen encounters coordinate antibiotic resilience in *Mycobacterium tuberculosis*. *Trends Microbiol.* **29**, 606–620 (2021).

16. Eisenreich, W., Rudel, T., Heesemann, J. & Goebel, W. Persistence of intracellular bacterial pathogens—With a focus on the metabolic perspective. *Front. Cell. Infect. Microbiol.* **10**, 25 (2021).

17. Baquero, F. & Levin, B. R. Proximate and ultimate causes of the bactericidal action of antibiotics. *Nat. Rev. Microbiol.* **19**, 123–132 (2020).

18. Lopatkin, A. J. et al. Bacterial metabolic state more accurately predicts antibiotic lethality than growth rate. *Nat. Microbiol.* **4**, 2109–2117 (2019).

19. Moreno-Gámez, S. et al. Wide lag time distributions break a trade-off between reproduction and survival in bacteria. *PNAS* **117**, 18729–18736 (2020).

20. Levin-Reisman, I. et al. Antibiotic tolerance facilitates the evolution of resistance. *Science* **355**, 826–830 (2017).

21. Sebastian, J. et al. De novo emergence of genetically resistant mutants of *Mycobacterium tuberculosis* from the persistence phase cells formed against antituberculosis drugs in vitro. *Antimicrob. Agents Chemother.* **61**, 25 (2017).

22. Levin-Reisman, I., Brauner, A., Ronin, I. & Balaban, N. Q. Epistasis between antibiotic tolerance, persistence, and resistance mutations. *PNAS* **116**, 14734–14739 (2019).

23. Windels, E. M. et al. Bacterial persistence promotes the evolution of antibiotic resistance by increasing survival and mutation rates. *ISME J.* **13**, 1239–1251 (2019).

24. Sulaiman, J. E. & Lam, H. Proteomic investigation of tolerant *Escherichia coli* populations from cyclic antibiotic treatment. *J. Prot. Res.* **19**, 900–913 (2020).

25. Liu, J., Gefen, O., Ronin, I., Bar-Meir, M. & Balaban, N. Q. Effect of tolerance on the evolution of antibiotic resistance under drug combinations. *Science* **367**, 200–204 (2020).

26. Davies, G. Pharmacologic considerations in use and development of antituberculosis drugs. *Cold Spring Harb. Perspect. Med.* **5**, a021170–a021170 (2014).

27. Sakatos, A. et al. Posttranslational modification of a histone-like protein regulates phenotypic resistance to isoniazid in mycobacteria. *Sci. Adv.* **4**, eaao1478 (2018).

28. Chua, H. C. et al. Combatting the rising tide of antimicrobial resistance: Pharmacokinetic/pharmacodynamic dosing strategies for maximal precision. *Int. J. Antimicrob. Agents* **57**, 106269 (2021).

29. Clewe, O., Faraj, A., Hu, Y., Coates, A. R. M. & Simonsson, U. S. H. A model-based analysis identifies differences in phenotypic resistance between in vitro and in vivo: Implications for translational medicine within tuberculosis. *J. Pharmacokinet. Pharmacodyn.* **47**, 421–430 (2020).

30. Fleischmann, S., Robben, C., Alter, T., Rossmanith, P. & Mester, P. How to evaluate non-growing cells-current strategies for determining antimicrobial resistance of VBNC bacteria. *Antibiotics (Basel)* **10**, 115 (2021).

31. Huemer, M., Mairpady Shambat, S., Brugger, S. D. & Zinkernagel, A. S. Antibiotic resistance and persistence—implications for human health and treatment perspectives. *EMBO Rep.* **21**, 25 (2020).

32. Oh, S., Trifonov, L., Yadav, V. D., Barry, C. E. & Boshoff, H. I. Tuberculosis drug discovery: A decade of hit assessment for defined targets. *Front. Cell. Infect. Microbiol.* **11**, 25 (2021).

33. Tiberi, S. et al. Tuberculosis: Progress and advances in development of new drugs, treatment regimens, and host-directed therapies. *Lancet Infect. Dis.* **18**, e183–e198 (2018).

34. Sackmann, E. K., Fulton, A. L. & Beebe, D. J. The present and future role of microfluidics in biomedical research. *Nature* **507**, 181–189 (2014).

35. Duncombe, T. A., Tentori, A. M. & Herr, A. E. Microfluidics: Reframing biological enquiry. *Nat. Rev. Mol. Cell Biol.* **16**, 554–567 (2015).

36. Golchin, S. A., Stratford, J., Curry, R. J. & McFadden, J. A microfluidic system for long-term time-lapse microscopy studies of mycobacteria. *Tuberculosis* **92**, 489–496 (2012).

37. Wakamoto, Y. et al. Dynamic persistence of antibiotic-stressed mycobacteria. *Science* **339**, 91–95 (2013).

38. Long, Z. et al. Measuring bacterial adaptation dynamics at the single-cell level using a microfluidic chemostat and time-lapse fluorescence microscopy. *Analyst* **139**, 5254–5262 (2014).

39. Okumus, B. et al. Mechanical slowing-down of cytoplasmic diffusion allows in vivo counting of proteins in individual cells. *Nat. Commun.* **7**, 25 (2016).

40. Bamford, R. A. et al. Investigating the physiology of viable but non-culturable bacteria by microfluidics and time-lapse microscopy. *BMC Biol.* **15**, 25 (2017).

41. Kaganovitch, E. et al. Microbial single-cell analysis in picoliter-sized batch cultivation chambers. *New Biotech.* **47**, 50–59 (2018).

42. Sanfilippo, J. E. et al. Microfluidic-based transcriptomics reveal force-independent bacterial rheosensing. *Nat. Microbiol.* **4**, 1274–1281 (2019).

43. Bakshi, S. et al. Tracking bacterial lineages in complex and dynamic environments with applications for growth control and persistence. *Nat. Microbiol.* **6**, 783–791 (2021).

44. Manina, G. & Dhar, N. Single-cell analysis of mycobacteria using microfluidics and time-lapse microscopy. *Methods Mol. Biol.* **2314**, 205–229 (2021).

45. Wang, H., Conover, G. M., Han, S.-I., Sacchettini, J. C. & Han, A. Development of single-cell-level microfluidic technology for long-term growth visualization of living cultures of *Mycobacterium smegmatis*. *Microsys. Nanoeng.* **7**, 25 (2021).

46. Potvin-Trottier, L., Luro, S. & Paulsson, J. Microfluidics and single-cell microscopy to study stochastic processes in bacteria. *Curr. Opin. Microbiol.* **43**, 186–192 (2018).

47. Scheler, O., Postek, W. & Garstecki, P. Recent developments of microfluidics as a tool for biotechnology and microbiology. *Curr. Opin. Biotech.* **55**, 60–67 (2019).
6. Tacconelli, E. et al. Discovery, research, and development of new antibiotics: The WHO priority list of antibiotic-resistant bacteria and tuberculosis. *Lancet Infect. Dis.* **18**, 318–327 (2018).

7. Sarathy, J. P. et al. Extreme drug tolerance of *Mycobacterium tuberculosis* in caseum. *Antimicrob. Agents Chemother.* **62**, 25 (2018).

8. Ordonez, A. A. et al. Dynamic imaging in patients with tuberculosis reveals heterogeneous drug exposures in pulmonary lesions. *Nat. Med.* **26**, 529–534 (2020).

9. Cicchese, J. M., Dartois, V., Kirschner, D. E. & Linderman, J. J. Both pharmacokinetic variability and granuloma heterogeneity impact the ability of the first-line antibiotics to sterilize tuberculosis granulomas. *Front. Pharmacol.* **11**, 25 (2020).

10. Jo, J., Price-Whelan, A. & Dietrich, L. E. P. Gradients and consequences of heterogeneity in biofilms. *Nat. Rev. Microbiol.* **25**, 25 (2022).

11. Dhar, N., McKinney, J. & Manina, G. Phenotypic heterogeneity in *Mycobacterium tuberculosis*. *Microbiol. Spectr.* **4**, TBTB2-0021–2016 (2016).

12. Larrouy-Maumus, G. et al. Cell-envelope remodeling as a determinant of phenotypic antibacterial tolerance in *Mycobacterium tuberculosis*. *ACS Infect. Dis.* **2**, 352–360 (2016).

13. Manina, G., Griego, A., Singh, L. K., McKinney, J. D. & Dhar, N. Preexisting variation in DNA damage response predicts the fate of single mycobacteria under stress. *EMBO J.* **38**, e101876 (2019).

14. Levien, E., Kondev, J. & Amir, A. The interplay of phenotypic variability and fitness in finite microbial populations. *J. R. Soc. Interface* **17**, 20190827 (2020).

15. Mishra, R., Yadav, V., Guha, M. & Singh, A. Heterogeneous host-pathogen encounters coordinate antibiotic resilience in *Mycobacterium tuberculosis*. *Trends Microbiol.* **29**, 606–620 (2021).

16. Eisenreich, W., Rudel, T., Heesemann, J. & Goebel, W. Persistence of intracellular bacterial pathogens—With a focus on the metabolic perspective. *Front. Cell. Infect. Microbiol.* **10**, 25 (2021).

17. Baquero, F. & Levin, B. R. Proximate and ultimate causes of the bactericidal action of antibiotics. *Nat. Rev. Microbiol.* **19**, 123–132 (2020).

18. Lopatkin, A. J. et al. Bacterial metabolic state more accurately predicts antibiotic lethality than growth rate. *Nat. Microbiol.* **4**, 2109–2117 (2019).

19. Moreno-Gámez, S. et al. Wide lag time distributions break a trade-off between reproduction and survival in bacteria. *PNAS* **117**, 18729–18736 (2020).

20. Levin-Reisman, I. et al. Antibiotic tolerance facilitates the evolution of resistance. *Science* **355**, 826–830 (2017).

21. Sebastian, J. et al. De novo emergence of genetically resistant mutants of *Mycobacterium tuberculosis* from the persistence phase cells formed against antituberculosis drugs in vitro. *Antimicrob. Agents Chemother.* **61**, 25 (2017).

22. Levin-Reisman, I., Brauner, A., Ronin, I. & Balaban, N. Q. Epistasis between antibiotic tolerance, persistence, and resistance mutations. *PNAS* **116**, 14734–14739 (2019).

23. Windels, E. M. et al. Bacterial persistence promotes the evolution of antibiotic resistance by increasing survival and mutation rates. *ISME J.* **13**, 1239–1251 (2019).

24. Sulaiman, J. E. & Lam, H. Proteomic investigation of tolerant *Escherichia coli* populations from cyclic antibiotic treatment. *J. Prot. Res.* **19**, 900–913 (2020).

25. Liu, J., Gefen, O., Ronin, I., Bar-Meir, M. & Balaban, N. Q. Effect of tolerance on the evolution of antibiotic resistance under drug combinations. *Science* **367**, 200–204 (2020).

26. Davies, G. Pharmacologic considerations in use and development of antituberculosis drugs. *Cold Spring Harb. Perspect. Med.* **5**, a021170–a021170 (2014).

27. Sakatos, A. et al. Posttranslational modification of a histone-like protein regulates phenotypic resistance to isoniazid in mycobacteria. *Sci. Adv.* **4**, eaao1478 (2018).

28. Chua, H. C. et al. Combatting the rising tide of antimicrobial resistance: Pharmacokinetic/pharmacodynamic dosing strategies for maximal precision. *Int. J. Antimicrob. Agents* **57**, 106269 (2021).

29. Clewe, O., Faraj, A., Hu, Y., Coates, A. R. M. & Simonsson, U. S. H. A model-based analysis identifies differences in phenotypic resistance between in vitro and in vivo: Implications for translational medicine within tuberculosis. *J. Pharmacokinet. Pharmacodyn.* **47**, 421–430 (2020).

30. Fleischmann, S., Robben, C., Alter, T., Rossmanith, P. & Mester, P. How to evaluate non-growing cells-current strategies for determining antimicrobial resistance of VBNC bacteria. *Antibiotics (Basel)* **10**, 115 (2021).

31. Huemer, M., Mairpady Shambat, S., Brugger, S. D. & Zinkernagel, A. S. Antibiotic resistance and persistence—implications for human health and treatment perspectives. *EMBO Rep.* **21**, 25 (2020).

32. Oh, S., Trifonov, L., Yadav, V. D., Barry, C. E. & Boshoff, H. I. Tuberculosis drug discovery: A decade of hit assessment for defined targets. *Front. Cell. Infect. Microbiol.* **11**, 25 (2021).

33. Tiberi, S. et al. Tuberculosis: Progress and advances in development of new drugs, treatment regimens, and host-directed therapies. *Lancet Infect. Dis.* **18**, e183–e198 (2018).

34. Sackmann, E. K., Fulton, A. L. & Beebe, D. J. The present and future role of microfluidics in biomedical research. *Nature* **507**, 181–189 (2014).

35. Duncombe, T. A., Tentori, A. M. & Herr, A. E. Microfluidics: Reframing biological enquiry. *Nat. Rev. Mol. Cell Biol.* **16**, 554–567 (2015).

36. Golchin, S. A., Stratford, J., Curry, R. J. & McFadden, J. A microfluidic system for long-term time-lapse microscopy studies of mycobacteria. *Tuberculosis* **92**, 489–496 (2012).

37. Wakamoto, Y. et al. Dynamic persistence of antibiotic-stressed mycobacteria. *Science* **339**, 91–95 (2013).

38. Long, Z. et al. Measuring bacterial adaptation dynamics at the single-cell level using a microfluidic chemostat and time-lapse fluorescence microscopy. *Analyst* **139**, 5254–5262 (2014).

39. Okumus, B. et al. Mechanical slowing-down of cytoplasmic diffusion allows in vivo counting of proteins in individual cells. *Nat. Commun.* **7**, 25 (2016).

40. Bamford, R. A. et al. Investigating the physiology of viable but non-culturable bacteria by microfluidics and time-lapse microscopy. *BMC Biol.* **15**, 25 (2017).

41. Kaganovitch, E. et al. Microbial single-cell analysis in picoliter-sized batch cultivation chambers. *New Biotech.* **47**, 50–59 (2018).

42. Sanfilippo, J. E. et al. Microfluidic-based transcriptomics reveal force-independent bacterial rhösensing. *Nat. Microbiol.* **4**, 1274–1281 (2019).

43. Bakshi, S. et al. Tracking bacterial lineages in complex and dynamic environments with applications for growth control and persistence. *Nat. Microbiol.* **6**, 783–791 (2021).

44. Manina, G. & Dhar, N. Single-cell analysis of mycobacteria using microfluidics and time-lapse microscopy. *Methods Mol. Biol.* **2314**, 205–229 (2021).

45. Wang, H., Conover, G. M., Han, S.-I., Sacchettini, J. C. & Han, A. Development of single-cell-level microfluidic technology for long-term growth visualization of living cultures of *Mycobacterium smegmatis*. *Microsys. Nanoeng.* **7**, 25 (2021).

46. Potvin-Trottier, L., Luro, S. & Paulsson, J. Microfluidics and single-cell microscopy to study stochastic processes in bacteria. *Curr. Opin. Microbiol.* **43**, 186–192 (2018).

47. Scheler, O., Postek, W. & Garstecki, P. Recent developments of microfluidics as a tool for biotechnology and microbiology. *Curr. Opin. Biotech.* **55**, 60–67 (2019).
48. Wang, P. et al. Robust growth of *Escherichia coli*. *Curr. Biol.* **20**, 1099–1103 (2010).

49. Changou, C. A. et al. Arginine starvation-associated typical cellular death involves mitochondrial dysfunction, nuclear DNA leakage, and chromatin autophagy. *PNAS* **111**, 14147–14152 (2014).

50. Hansen, A. S., Hao, N. & O'Shea, E. K. High-throughput microfluidics to control and measure signaling dynamics in single yeast cells. *Nat. Protoc.* **10**, 1181–1197 (2015).

51. Melin, J. & Quake, S. R. Microfluidic large-scale integration: The evolution of design rules for biological automation. *Annu. Rev. Biophys. Biomolec. Struct.* **36**, 213–231 (2007).

52. Kim, J. et al. A programmable microfluidic cell array for combinatorial drug screening. *Lab Chip.* **12**, 1813 (2012).

53. Stratz, S., Eyer, K., Kurth, F. & Dittrich, P. S. On-chip enzyme quantification of single escherichia coli bacteria by immunoassay-based analysis. *Analyt. Chem.* **86**, 12375–12381 (2014).

54. Lareau, C. A. et al. Droplet-based combinatorial indexing for massive-scale single-cell chromatin accessibility. *Nat. Biotech.* **37**, 916–924 (2019).

55. Tomasi, R. F. X., Sart, S., Champetier, T. & Baroud, C. N. Individual control and quantification of 3D spheroids in a high-density microfluidic droplet array. *Cell Rep.* **31**, 107670 (2020).

56. Wang, X., Liu, Z. & Pang, Y. Concentration gradient generation methods based on microfluidic systems. *RCS Adv.* **7**, 29966–29984 (2017).

57. Unger, M. A., Chou, H. P., Thorsen, T., Scherer, A. & Quake, S. R. Monolithic microfabricated valves and pumps by multilayer soft lithography. *Science* **288**, 113–116 (2000).

58. Tan, S. H., Nguyen, N. T., Chua, Y. C. & Kang, T. G. Oxygen plasma treatment for reducing hydrophobicity of a sealed polydimethylsiloxane microchannel. *Biomicrofluidics* **4**, 32204 (2010).

59. Thorsen, T., Maerkl, S. J. & Quake, S. R. Microfluidic large-scale integration. *Science* **298**, 580–584 (2002).

60. Jeon, N. L. et al. Generation of solution and surface gradients using microfluidic systems. *Langmuir* **16**, 8311–8316 (2000).

61. Lin, F. et al. Generation of dynamic temporal and spatial concentration gradients using microfluidic devices. *Lab Chip.* **4**, 164 (2004).

62. Óbrink-Hansen, K. et al. Moxifloxacin pharmacokinetic profile and efficacy evaluation in empiric treatment of community-acquired pneumonia. *Antimicrob. Agents Chemother.* **59**, 2398–2404 (2015).

63. Iacobino, A., Piccaro, G., Pardini, M., Fattorini, L. & Giannoni, F. Moxifloxacin activates the SOS response in *Mycobacterium tuberculosis* in a dose- and time-dependent manner. *Microorganisms* **9**, 255 (2021).

64. Meylan, S., Andrews, I. W. & Collins, J. J. Targeting antibiotic tolerance, pathogen by pathogen. *Cell* **172**, 1228–1238 (2018).

65. Borók, A., Laboda, K. & Bonyár, A. PDMS bonding technologies for microfluidic applications: A review. *Biosensors* **11**, 292 (2021).

66. Nania, M., Matar, O. K. & Cabral, J. T. Frontal vitrification of PDMS using air plasma and consequences for surface wrinkling. *Soft Matter* **11**, 3067–3075 (2015).

67. Gökaltun, A., Kang, Y. B., Yarmush, M. L., Usta, O. B. & Asatekin, A. Simple surface modification of poly(dimethylsiloxane) via surface segregating smart polymers for biomicrofluidics. *Sci. Rep.* **9**, 25 (2019).

68. Kristiansen, P. M., Karpik, A., Werder, J., Guilherme, M. & Grob, M. Thermoplastic microfluidics. *Methods Mol. Biol.* **2373**, 39–55 (2022).

69. Scott, S. M. & Ali, Z. Fabrication methods for microfluidic devices: An overview. *Micromachines* **12**, 319 (2021).

70. Nielsen, A. V., Beauchamp, M. J., Nordin, G. P. & Woolley, A. T. 3D printed microfluidics. *Annu. Rev. Analyt. Chem.* **13**, 45–65 (2020).

71. Rossi, G., Manfrin, A. & Lutolf, M. P. Progress and potential in organoid research. *Nat. Rev. Gen.* **19**, 671–687 (2018).

72. Jalili-Firoozinezhad, S., Miranda, C. C. & Cabral, J. M. S. Modeling the human body on microfluidic chips. *Trends Biotech.* **39**, 838–852 (2021).

73. Woodruff, K. & Maerkl, S. J. Microfluidic module for real-time generation of complex multimolecule temporal concentration profiles. *Anal. Chem.* **90**, 696–701 (2017).

74. Komen, J. et al. Controlled pharmacokinetic anti-cancer drug concentration profiles lead to growth inhibition of colorectal cancer cells in a microfluidic device. *Lab Chip.* **20**, 3167–3178 (2020).

75. Zhang, T. et al. A microfluidic concentration gradient maker with tunable concentration profiles by changing feed flow rate ratios. *Micromachines* **11**, 284 (2020).

76. Thurgood, P. et al. Tunable harmonic flow patterns in microfluidic systems through simple tube oscillation. *Small* **16**, e2003612 (2020).

77. Zhu, J. Y. et al. Reconfigurable, self-sufficient convective heat exchanger for temperature control of microfluidic systems. *Anal. Chem.* **91**, 15784–15790 (2019).

78. Ernest, J. P. et al. Development of new tuberculosis drugs: Translation to regimen composition for drug-sensitive and multidrug-resistant tuberculosis. *Annu. Rev. Pharmacol. Toxicol.* **61**, 495–516 (2021).

79. Gumbo, T. et al. Selection of a moxifloxacin dose that suppresses drug resistance in *Mycobacterium tuberculosis*, by use of an in vitro pharmacodynamic infection model and mathematical modeling. *J. Infect. Dis.* **190**, 1642–1651 (2004).

80. Gumbo, T., Pasipanodya, J. G., Nuermberger, E., Romero, K. & Hanna, D. Correlations between the hollow fiber model of tuberculosis and therapeutic events in tuberculosis patients: Learn and confirm. *Clin. Infect. Dis.* **61**, S18–S24 (2015).

81. Kloprogge, F., Hammond, R., Kipper, K., Gillespie, S. H. & Della Pasqua, O. Mimicking in-vivo exposures to drug combinations in-vitro: Anti-tuberculosis drugs in lung lesions and the hollow fiber model of infection. *Sci. Rep.* **9**, 25 (2019).

82. Martinez, A. & Abel zur Wiesch, P. Estimating treatment prolongation for persistent infections. *Pathog. Dis.* **76**, 25 (2018).

83. Ocaña, J. A., Barragán, F. J. & Callejón, M. Spectrofluorimetric determination of moxifloxacin in tablets, human urine and serum. *Analyst.* **125**, 2322–2325 (2000).

84. Huang, L.-K. & Wang, M.-J.J. Image thresholding by minimizing the measures of fuzziness. *Pattern Recogn.* **28**, 41–51 (1995).

85. Prewitt, J. M. S. & Mendelsohn, M. L. The analysis of cell images*. *Ann. N. Y. Acad. Sci.* **128**, 1035–1053 (2006).

86. Hailey, D. W., Davis, T. N. & Muller, E. G. Fluorescence resonance energy transfer using color variants of green fluorescent protein. *Meth. Enzymol.* **351**, 34–49 (2002).

## Acknowledgements

We thank Pietro Slavich (LPTHE) for reading of the manuscript and critical discussions. We thank Théo Guillerm for his contribution in analyzing the gradient formation with dyes. We thank Gaël Millot (Bioinformatics and Biostatistics Hub of the Institut Pasteur), for support in writing the R script for managing .csv files. We thank Samy Gobaa for providing us access to the Biomaterials and Microfluidics core facility of the Institut Pasteur. This work was supported by the following Grants to G.M.: French Medical Research Foundation Grant ING20160435202; French National Research Agency Grant ANR-17-CE11-0007-01; French National Research Agency Grant ANR-10-LABX-62-IBEID; IMI 2 Joint Undertaking Grant 853989, receiving support from the European Union's Horizon 2020 research and innovation programme and EFPIA and Global Alliance for TB Drug Development non-profit organization, Bill & Melinda Gates Foundation, University of Dundee; Institut Pasteur core funding.
48. Wang, P. et al. Robust growth of *Escherichia coli*. *Curr. Biol.* **20**, 1099–1103 (2010).

49. Changou, C. A. et al. Arginine starvation-associated typical cellular death involves mitochondrial dysfunction, nuclear DNA leakage, and chromatin autophagy. *PNAS* **111**, 14147–14152 (2014).

50. Hansen, A. S., Hao, N. & O'Shea, E. K. High-throughput microfluidics to control and measure signaling dynamics in single yeast cells. *Nat. Protoc.* **10**, 1181–1197 (2015).

51. Melin, J. & Quake, S. R. Microfluidic large-scale integration: The evolution of design rules for biological automation. *Annu. Rev. Biophys. Biomolec. Struct.* **36**, 213–231 (2007).

52. Kim, J. et al. A programmable microfluidic cell array for combinatorial drug screening. *Lab Chip.* **12**, 1813 (2012).

53. Stratz, S., Eyer, K., Kurth, F. & Dittrich, P. S. On-chip enzyme quantification of single escherichia coli bacteria by immunoassay-based analysis. *Analyt. Chem.* **86**, 12375–12381 (2014).

54. Lareau, C. A. et al. Droplet-based combinatorial indexing for massive-scale single-cell chromatin accessibility. *Nat. Biotech.* **37**, 916–924 (2019).

55. Tomasi, R. F. X., Sart, S., Champetier, T. & Baroud, C. N. Individual control and quantification of 3D spheroids in a high-density microfluidic droplet array. *Cell Rep.* **31**, 107670 (2020).

56. Wang, X., Liu, Z. & Pang, Y. Concentration gradient generation methods based on microfluidic systems. *RCS Adv.* **7**, 29966–29984 (2017).

57. Unger, M. A., Chou, H. P., Thorsen, T., Scherer, A. & Quake, S. R. Monolithic microfabricated valves and pumps by multilayer soft lithography. *Science* **288**, 113–116 (2000).

58. Tan, S. H., Nguyen, N. T., Chua, Y. C. & Kang, T. G. Oxygen plasma treatment for reducing hydrophobicity of a sealed polydimethylsiloxane microchannel. *Biomicrofluidics* **4**, 32204 (2010).

59. Thorsen, T., Maerkl, S. J. & Quake, S. R. Microfluidic large-scale integration. *Science* **298**, 580–584 (2002).

60. Jeon, N. L. et al. Generation of solution and surface gradients using microfluidic systems. *Langmuir* **16**, 8311–8316 (2000).

61. Lin, F. et al. Generation of dynamic temporal and spatial concentration gradients using microfluidic devices. *Lab Chip.* **4**, 164 (2004).

62. Óbrink-Hansen, K. et al. Moxifloxacin pharmacokinetic profile and efficacy evaluation in empiric treatment of community-acquired pneumonia. *Antimicrob. Agents Chemother.* **59**, 2398–2404 (2015).

63. Iacobino, A., Piccaro, G., Pardini, M., Fattorini, L. & Giannoni, F. Moxifloxacin activates the SOS response in *Mycobacterium tuberculosis* in a dose- and time-dependent manner. *Microorganisms* **9**, 255 (2021).

64. Meylan, S., Andrews, I. W. & Collins, J. J. Targeting antibiotic tolerance, pathogen by pathogen. *Cell* **172**, 1228–1238 (2018).

65. Borók, A., Laboda, K. & Bonyár, A. PDMS bonding technologies for microfluidic applications: A review. *Biosensors* **11**, 292 (2021).

66. Nania, M., Matar, O. K. & Cabral, J. T. Frontal vitrification of PDMS using air plasma and consequences for surface wrinkling. *Soft Matter* **11**, 3067–3075 (2015).

67. Gökaltun, A., Kang, Y. B., Yarmush, M. L., Usta, O. B. & Asatekin, A. Simple surface modification of poly(dimethylsiloxane) via surface segregating smart polymers for biomicrofluidics. *Sci. Rep.* **9**, 25 (2019).

68. Kristiansen, P. M., Karpik, A., Werder, J., Guilherme, M. & Grob, M. Thermoplastic microfluidics. *Methods Mol. Biol.* **2373**, 39–55 (2022).

69. Scott, S. M. & Ali, Z. Fabrication methods for microfluidic devices: An overview. *Micromachines* **12**, 319 (2021).

70. Nielsen, A. V., Beauchamp, M. J., Nordin, G. P. & Woolley, A. T. 3D printed microfluidics. *Annu. Rev. Analyt. Chem.* **13**, 45–65 (2020).

71. Rossi, G., Manfrin, A. & Lutolf, M. P. Progress and potential in organoid research. *Nat. Rev. Gen.* **19**, 671–687 (2018).

72. Jalili-Firoozinezhad, S., Miranda, C. C. & Cabral, J. M. S. Modeling the human body on microfluidic chips. *Trends Biotech.* **39**, 838–852 (2021).

73. Woodruff, K. & Maerkl, S. J. Microfluidic module for real-time generation of complex multimolecule temporal concentration profiles. *Anal. Chem.* **90**, 696–701 (2017).

74. Komen, J. et al. Controlled pharmacokinetic anti-cancer drug concentration profiles lead to growth inhibition of colorectal cancer cells in a microfluidic device. *Lab Chip.* **20**, 3167–3178 (2020).

75. Zhang, T. et al. A microfluidic concentration gradient maker with tunable concentration profiles by changing feed flow rate ratios. *Micromachines* **11**, 284 (2020).

76. Thurgood, P. et al. Tunable harmonic flow patterns in microfluidic systems through simple tube oscillation. *Small* **16**, e2003612 (2020).

77. Zhu, J. Y. et al. Reconfigurable, self-sufficient convective heat exchanger for temperature control of microfluidic systems. *Anal. Chem.* **91**, 15784–15790 (2019).

78. Ernest, J. P. et al. Development of new tuberculosis drugs: Translation to regimen composition for drug-sensitive and multidrug-resistant tuberculosis. *Annu. Rev. Pharmacol. Toxicol.* **61**, 495–516 (2021).

79. Gumbo, T. et al. Selection of a moxifloxacin dose that suppresses drug resistance in *Mycobacterium tuberculosis*, by use of an in vitro pharmacodynamic infection model and mathematical modeling. *J. Infect. Dis.* **190**, 1642–1651 (2004).

80. Gumbo, T., Pasipanodya, J. G., Nuermberger, E., Romero, K. & Hanna, D. Correlations between the hollow fiber model of tuberculosis and therapeutic events in tuberculosis patients: Learn and confirm. *Clin. Infect. Dis.* **61**, S18–S24 (2015).

81. Kloprogge, F., Hammond, R., Kipper, K., Gillespie, S. H. & Della Pasqua, O. Mimicking in-vivo exposures to drug combinations in-vitro: Anti-tuberculosis drugs in lung lesions and the hollow fiber model of infection. *Sci. Rep.* **9**, 25 (2019).

82. Martinez, A. & Abel zur Wiesch, P. Estimating treatment prolongation for persistent infections. *Pathog. Dis.* **76**, 25 (2018).

83. Ocaña, J. A., Barragán, F. J. & Callejón, M. Spectrofluorimetric determination of moxifloxacin in tablets, human urine and serum. *Analyst.* **125**, 2322–2325 (2000).

84. Huang, L.-K. & Wang, M.-J.J. Image thresholding by minimizing the measures of fuzziness. *Pattern Recogn.* **28**, 41–51 (1995).

85. Prewitt, J. M. S. & Mendelsohn, M. L. The analysis of cell images*. *Ann. N. Y. Acad. Sci.* **128**, 1035–1053 (2006).

86. Hailey, D. W., Davis, T. N. & Muller, E. G. Fluorescence resonance energy transfer using color variants of green fluorescent protein. *Meth. Enzymol.* **351**, 34–49 (2002).

## Acknowledgements

We thank Pietro Slavich (LPTHE) for reading of the manuscript and critical discussions. We thank Théo Guillerm for his contribution in analyzing the gradient formation with dyes. We thank Gaël Millot (Bioinformatics and Biostatistics Hub of the Institut Pasteur), for support in writing the R script for managing .csv files. We thank Samy Gobaa for providing us access to the Biomaterials and Microfluidics core facility of the Institut Pasteur. This work was supported by the following Grants to G.M.: French Medical Research Foundation Grant ING20160435202; French National Research Agency Grant ANR-17-CE11-0007-01; French National Research Agency Grant ANR-10-LABX-62-IBEID; IMI 2 Joint Undertaking Grant 853989, receiving support from the European Union's Horizon 2020 research and innovation programme and EFPIA and Global Alliance for TB Drug Development non-profit organization, Bill & Melinda Gates Foundation, University of Dundee; Institut Pasteur core funding.
**Author contributions**

Conceptualization: G.M. and M.M.; methodology: M.M. and G.M.; investigation: M.M., N.G. and G.M.; resources: M.M. and G.M.; data curation: M.M.; formal analysis: M.M., N.G. and G.M.; visualization: M.M. and G.M.; writing—original draft: G.M.; writing—review and editing: M.M., N.G. and G.M.; supervision: G.M.; project administration: G.M.; funding acquisition: G.M.

**Competing interests**

G.M. and M.M. are designated as inventors in the pending international patent application WO 2020/229629 filed by the Institut Pasteur. This patent application includes the microfluidic device for use in single cell analysis described in the manuscript.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-24175-9.

**Correspondence** and requests for materials should be addressed to G.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022
**Author contributions**

Conceptualization: G.M. and M.M.; methodology: M.M. and G.M.; investigation: M.M., N.G. and G.M.; resources: M.M. and G.M.; data curation: M.M.; formal analysis: M.M., N.G. and G.M.; visualization: M.M. and G.M.; writing—original draft: G.M.; writing—review and editing: M.M., N.G. and G.M.; supervision: G.M.; project administration: G.M.; funding acquisition: G.M.

**Competing interests**

G.M. and M.M. are designated as inventors in the pending international patent application WO 2020/229629 filed by the Institut Pasteur. This patent application includes the microfluidic device for use in single cell analysis described in the manuscript.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-24175-9.

Correspondence and requests for materials should be addressed to G.M.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022