Authors: Mohammed Nouaman, Alexis Darras, Christian Wagner, Steffen M. Recktenwald
Categories: Regular Articles
Source: Biomicrofluidics
Doi: 10.1063/5.0197208
The ability to change shape is essential for the proper functioning of red blood cells (RBCs) within the microvasculature. The shape of RBCs significantly influences blood flow and has been employed in microfluidic lab-on-a-chip devices, serving as a diagnostic biomarker for specific pathologies and enabling the assessment of RBC deformability. While external flow conditions, such as the vessel size and the flow velocity, are known to impact microscale RBC flow, our comprehensive understanding of how their shape-adapting ability is influenced by channel confinement in biomedical applications remains incomplete. This study explores the impact of various rectangular and square channels, each with different confinement and aspect ratios, on the in vitro RBC flow behavior and characteristic shapes. We demonstrate that rectangular microchannels, with a height similar to the RBC diameter in combination with a confinement ratio exceeding 0.9, are required to generate distinctive well-defined croissant and slipper-like RBC shapes. These shapes are characterized by their equilibrium positions in the channel cross section, and we observe a strong elongation of both stable shapes in response to the shear rate across the different channels. Less confined channel configurations lead to the emergence of unstable other shape types that display rich shape dynamics. Our work establishes an experimental framework to understand the influence of channel size on the single-cell flow behavior of RBCs, providing valuable insights for the design of biomicrofluidic single-cell analysis applications.
Red blood cells (RBCs) are the main cellular constituent of blood and are vital in facilitating gas exchange between blood and tissues within the microcirculation. The properties of RBCs, such as their deformability and shape, have a profound impact on blood flow.^1^ At rest, human RBCs have a disk-like shape with a diameter of 8μm and a thickness of 2μm. Due to their high deformability, healthy RBCs can dynamically adapt their shape according to external flow conditions, such as the flow rate, vessel confinement, and the rheological properties of the surrounding fluid.^2–4^ In the microvascular network where vessel diameters are comparable to RBC size, RBCs flow in a single-file arrangement, and various RBC shapes have been observed.^5–7^ Consequently, the study of single-cell RBC flow has been conducted experimentally in microfluidic devices,^8–15^ as well as through numerical simulations,^16–21^ to understand the impact of intrinsic cell properties and external flow conditions on cell shape. Understanding single-cell RBC flow and deformation has contributed to the advancement and development of lab-on-a-chip technologies for RBC deformability analysis of storage lesions in transfusion medicine^22–24^ and for evaluation of biomedical RBC properties in health and disease.^25–27^ Nevertheless, fundamental knowledge of how external flow conditions such as channel confinement modify the RBC shape in in vitro microscale flows is still missing.
In microscale single-cell flow, RBCs display a wide variety of stable and dynamic shapes depending on the biophysical cell properties, channel confinement, flow velocity, and the properties of the surrounding medium. Various dynamical states, including snaking, tumbling, swinging, and tank-treading cells, have been reported under steady flow conditions.^1,28,29^ Notably, in strongly confined rectangular channels with channel dimensions (height and width) similar to the RBC size, this RBC shape complexity consolidates into two dominant RBC the croissant and the slipper shape.^12,19,24^ The symmetric croissant-like shape [Fig. 1(a)] predominantly appears at low flow velocities, while the asymmetric slipper shape [Fig. 1(b)] emerges for higher cell velocities. Guckenberger et al.^19^ investigated these two dominant RBC shapes in a rectangular microfluidic channel with a width of 12μm and a height of 10μm through a combination of microfluidic experiments and numerical simulations. They introduced the concept of the so-called RBC shape phase diagram, which shows the fraction of these stable croissant and slipper shapes, as well as a class called “others” [Fig. 1(c)] that were not uniquely identifiable, as a function of the cell velocity v.^19^ In their rectangular channel, they observed that the croissant shape dominated at cell velocities v<5mm/s, while roughly 70%-80% of the RBCs exhibited a slipper shape above v>5mm/s. They also reported a strong flow-shape coupling of the highly deformable RBCs and showed that the cell’s equilibrium position in the cross section and along the channel width is inherently coupled to its shape. Hence, croissant-like RBCs preferentially flow at a central position in the channel centerline, whereas slippers flow closer to the channel sidewalls (see green marker in Fig. 1). Furthermore, the authors demonstrated that the emergence of RBC shapes is influenced not only by system parameters, such as flow velocity or channel size, but also by initial conditions, including the initial shape of the RBC and its position within the channel cross section at the onset of the microfluidic channel.^19^
FIG. 1. Representative RBC shapes in confined flows. (a) Centered croissant-like shape, (b) off-centered slipper-like shape, and (c) other shape in a microfluidic channel with W=10μm and H=8μm. Flow is from left to right. Black dashed lines correspond to the channel centerline across the channel width W. Green dots indicate the cell’s center of mass with respect to the y-direction. (Scale bars, 5μm.)
Recently, microfluidic investigations have unveiled the sensitivity of the stable croissant and slipper shapes in relation to biophysical cell properties such as RBC age and reduced membrane deformability.^15^ Utilizing RBCs from healthy donors and employing density gradient centrifugation methods, it has been demonstrated that the proportion of stable asymmetric, off-centered slipper-like cells diminishes with increasing age, while aged cells exhibit an augmented prevalence of stable symmetric croissants along the microchannel centerline. Consequently, analysis of the RBC shape phase diagram facilitates the discrimination of distinct RBC sub-populations, notably revealing variations in cell age.
Expanding upon the shape phase diagram proposed by Guckenberger et al.,^19^ additional RBC shapes in microcapillary flow, such as sphero-echinocytes, acanthocytes, and complementary pathological croissant and slipper shapes, have been introduced under pathological conditions to complement the phase diagram.^30,31^ Integrating the analysis of the cell’s equilibrium position, these cell shape classification approaches have been recently applied to assess changes in the microcapillary flow behavior of neuroacanthocytosis syndrome and COVID-19 patients and patients undergoing dialysis,^24,30–32^ demonstrating its potential as a biomarker and functional diagnostic tool for specific pathologies and to evaluate the quality of stored blood.
In addition to experiments conducted in rectangular microchannels, numerical simulations revealed a diverse range of steady and dynamic RBC shapes in round microcapillaries. Fedosov et al.^16^ constructed a comprehensive shape phase diagram encompassing stationery parachutes, swinging slippers, tumbling, and snaking discocytes as a function of the channel confinement and the shear rate in the system with a circular cross section. Note that parachute-shaped RBCs, which tend to become more asymmetric in rectangular channels, are often referred to as croissants in microfluidic experiments with square or rectangular microchannels. Similar to experimental observations,^19^ parachute-like RBCs in the round capillaries flow in the tube center, while slippers are displaced further from the channel center.^16^ The conditions of occurrence of these shapes and the transition between the dynamic shapes critically depend on the RBC properties and the viscosity contrast between RBC cytosol and blood plasma.^21^
While experimental and numerical progress has advanced our understanding of how intrinsic cell properties, such as their deformability, cytosol viscosity, viscoelasticity of the membrane, and age, affect the RBC shape and microscale flow behavior,^15,33,34^ detailed experimental investigations of how slight changes in the channel confinement affect the RBC capillary flow and the shape phase diagram remain scarce. Presently, poly(dimethylsiloxane) (PDMS) microchannels fabricated using soft lithography techniques serve as a common platform to investigate the microscale flow of biological systems.^35,36^ These techniques generally offer highly reproducible results down to the nanoscale, with low shrinkage during cure.^37^ However, inconsistencies in PDMS casting, curing, releasing, or bonding can introduce variations in the channel dimensions between different molds or chips during fabrication. Additionally, changes in the channel cross section may occur due to the deformation of flexible PDMS microchannels under pressure-driven flow.^38–42^ Therefore, understanding the effect of channel confinement on the different RBC shapes, employed for in vitro microvascular flow assessment, is paramount. This becomes particularly crucial in the context of assessing microscale flow under physiologically relevant conditions, where subtle variations in channel confinement can significantly influence the observed RBC shapes and, consequently, the accuracy of the RBC shape phase diagram.
In this study, we examine the single-cell flow of RBCs through six distinct microfluidic channels. Specifically, we explore the effect of how small differences in the order of a few micrometers affect the RBC phase diagram, the cells’ equilibrium positions, their elongation, and the occurrence of unstable cell shapes. Our study aims to conduct a comprehensive experimental investigation to explore the sensitivity of these approaches, particularly how the RBC shape phase diagrams (PDs) are influenced by the precise channel dimensions. This understanding is crucial for elucidating RBC flow dynamics in more intricate systems, such as those found in pathological conditions. For this, we use high-speed imaging of microscale RBC flow and systematically vary the channel width and height, resulting in different aspect ratios and confinement ratios of the channel cross section. Our findings highlight that the presence of stable slipper-like shapes requires a non-square rectangular cross section, while croissant-like cells also disappear gradually in square channels. Intriguingly, both stable croissant and slipper shapes disappear when the channel dimensions exceed 10μm, leading to the emergence of highly dynamical other RBC shapes. Additionally, we show how the elongation of croissants and slippers, quantified by the elongation index EI, scales with the shear rate in the microchannels, highly relevant for biomicrofluidic technologies aimed at measuring RBC deformability.^27^
Blood was collected with informed consent from three healthy voluntary donors (age 28–51 years) through needle prick and it was subsequently suspended in phosphate-buffered saline solution (Gibco PBS, Fisher Scientific, Schwerte, Germany). Following collection, samples were centrifuged for 5 min at 3000 × g. The sedimented RBCs were washed three times with PBS and final RBC samples were adjusted at a hematocrit concentration of 1%Ht in a PBS solution that contained 1g/l bovine serum albumin (BSA, Sigma-Aldrich, Taufkirchen, Germany). The PBS/BSA mixture is a Newtonian fluid with a shear viscosity of roughly η=1.2mPas.^13,42^
Blood withdrawal, sample preparation, and microfluidic experiments were performed according to the guidelines of the Declaration of Helsinki and approved by the ethics committee of the “Ärztekammer des Saarlandes” (permission number 51/18).
RBC suspensions were pumped through six distinct microfluidic chips using a high-precision pressure device (OB1-MK3, Elveflow, Paris, France) to apply constant pressure drops ranging between p=50-1000mbar. The microfluidic chips were designed with microchannel having widths of either W≈10μm or W≈15μm in combination with three different heights H≈8,10,15μm (Table I and Fig. 2). The length of all microfluidic channels is L=40mm. Channel dimensions were determined by a customized MATLAB algorithm [9.14.0.2206163 (R2023a), The MathWorks, Natick, MA] that detects the channel borders based on microscopic images obtained using brightfield microscopy. Channel dimension data were averaged from different microfluidic chips as well as different positions along the channel flow direction within a chip.
FIG. 2. Schematic representation of the used microfluidic channels. We employ different channels with a width of (a) 10μm and (b) 15μm in combination with different channel heights. Blue arrows indicate the optical access in the microfluidic setup. The red circle illustrates the size of an undeformed discocyte-shaped RBC with a diameter of 8μm relative to the used channels.
We define the channel aspect ratio AR=W/H and the confinement ratio χ=DRBC/Dh with the diameter of a discocyte-shaped RBC at rest DRBC≈8μm and the hydraulic diameter of the rectangular channel Dh=2WH/(W+H). The microfluidic chip fabrication followed standard soft lithography techniques using polydimethylsiloxane (PDMS, RTV 615A/B, Momentive Performance Materials, Waterford, NY, USA).^37^ Subsequently, the chip was bonded to a glass slide using a plasma cleaner (PDC-32G, Harrick Plasma, Ithaca, NY, USA). The inlet and the outlet of the microfluidic chips were connected with rigid medical-grade polyethylene tubing ( 0.86mm inner diameter, Scientific Commodities, Lake Havasu City, AZ, USA) to the sample and waste containers, respectively.
The microfluidic device was mounted on an inverted microscope (Eclipse TE2000-S, Nikon, Melville, NY, USA), featuring red LED illumination (SOLIS-415C, Thorlabs Inc., Newton, NJ, USA), a 60× air objective (Plan Fluor, Nikon, Melville, NY, USA) with a numerical aperture NA=1.25. RBC flow was recorded in the middle of the microchannels at L/2 using a high-speed camera (Fastec HiSpec 2G, FASTEC Imaging, San Diego, CA, USA). All microfluidic experiments were performed at 22°C.
Image sequences were post-processed using a customized MATLAB algorithm. For each single RBC, we determined the center of mass of each cell in the projection plane (Fig. 1), length a in the flow direction, and diameter b along the channel width by identifying a bounding box around the RBC shape. The cell’s elongation index is calculated as EI=(a−b)/(a+b). To determine individual cell velocities, we tracked the cell position throughout the image sequence within the field of view. A frame rate of up to 400 frames per second was used to record image sequences of RBC passing the field of view. RBC shapes in flow were classified manually following the criteria established by Guckenberger et al.^19^ Since we did not observe inter-individual variations in the results, data were averaged between the three healthy donors. Data analysis was performed on an average of 6714 cells per donor (between 2653 and 10 635 cells).
Our investigation centers on exploring the impact of channel confinement on the RBC shape in microchannels, focusing on the three dominant RBC shape classes, namely, croissants, slippers, and others, similar to the previously established phase diagram RBC shapes.^12,19,24,31^ Based on the applied pressure drop and the channel cross section, the resulting cell velocities in the microchannels are in the range of v=0.2-35mm/s (Fig. 3), with the lower range being similar to the physiological flow velocities in the microvascular network.^3,4^ Additionally, we calculated the nondimensional capillary number,
where η=1.2mPas is the shear viscosity of the surrounding fluid,^13,42^
γ˙=6v/Dh the wall shear rate in the channel,^43^
a=DRBC/2 the discoid radius of the RBC, and Gs the membrane shear elastic modulus.^44^ We use Gs=4μN/m, in agreement with previous studies.^45,46^ In the investigated cell velocity regime, we find Ca≈0.28-17 (see top x axes in Fig. 3).
FIG. 3. Shape phase diagrams (PDs) of RBCs in microfluidic (a) W=10μm and H=8μm, (b) W=10μm and H=10μm, (c) W=10μm and H=15μm, (d) W=15μm and H=8μm, (e) W=15μm and H=10μm, and (f) W=15μm and H=15μm. Fraction of croissant-like, slipper-like, and other RBC shapes as a function of the cell velocity (bottom x axes) and capillary number Ca (top x axes) in the different microfluidic channels.
In the RBC shape phase diagrams, we observe a prominent proportion of stable croissant and slipper-shaped RBCs for the smallest channel cross section [Fig. 3(a)] consistent with studies of RBC flow in similarly confined channels.^12,19^ The croissant fraction reaches a peak value of roughly 50% at a velocity of v≈1mm/s. As the cell velocity increases, the fraction of croissant-shaped RBC continuously decreases. Simultaneously, the amount of slipper-shaped RBCs increases above v>3mm/s, eventually reaching a plateau value at 70%-75% above 5mm/s. Keeping the channel width fixed and increasing the channel height to H=10μm results in a higher fraction of croissants while slipper-shaped RBCs disappear [Fig. 3(b)]. However, with a further increase of the channel height to H=15μm, the croissant fraction decreases again, and most RBCs display other shape types [Fig. 3(c)].
In the channels with a larger width of W=15μm and a height of H=8μm, RBC also exhibit both croissant and slipper-shaped RBCs [Fig. 3(d)] and a qualitatively similar phase diagram than for the smallest channel cross section [compare Fig. 3(a)]. However, the fraction of both stable shapes notably reduces in the wider channel. The prominent croissant peak reaches merely 20% at v≈1mm/s, while the slipper plateau saturates at 30%-40% above 5mm/s. Increasing the channel height results in a successive suppression of slipper-like cells [Fig. 3(e)]. At W = 15 µm and H = 15 µm, the occurrence of slipper-like RBCs is ultimately suppressed completely, whereas most cells exhibit other shapes at all investigated velocities with a few croissant-like RBC remaining [Fig. 3(f)].
Taken together, our observations emphasize that an oblong rectangular cross section coupled with a shallow channel ( H<10μm) is a prerequisite for the emergence of slippers [see Figs. 3(a) and 3(d)]. Modestly increasing the channel height by 2μm to H=10μm results in a significant reduction of slipper fraction [see Fig. 3(e)]. Consequently, stable slipper-like shapes are primarily found in strongly confined, rectangular channels with a height smaller than 10μm and a confinement ratio exceeding χ≈0.9. It is noteworthy that these slippers are categorically different than the tank-treading slipper-shaped RBCs reported in numerical simulation within round capillaries.^16,21^ In the context of round capillaries, slippers are essentially absent at χ>0.7 due to the cylindrical channel geometry and mainly appear at low confinement ratios compared to their parachute-like counterpart.
In square channels, we do not observe slipper-shaped cells [see Figs. 3(b) and 3(f)], but an increase in symmetric croissant-like RBCs when the channel dimensions do not exceed 10μm, corresponding to a confinement ratio of χ=0.75 [see Fig. 3(b)]. At χ=0.52 [see Fig. 3(f)], the fraction of croissants is drastically reduced, which is similar to observations in numerical simulations, which showed that a transition between stationary parachutes and other dynamical shapes occurs at χ≈0.5-0.55 in square microchannels.^20^
Note that our microfluidic experiments are restricted to optically accessing RBC flow through one channel side (see blue arrows in Fig. 2). As a result, cells will be classified based on their projection in the optical x- y-plane in the channel (see Fig. 2). Therefore, despite having two channels with the same dimensions swapped by 90° [Figs. 3(c) and 3(e)], we observe different shape phase diagrams. For instance, slippers detected in the channel with W=15μm and H=10μm flow at off-centered position along channel width W in y-direction in the x- y-plane. In the channel with W=10μm and H=15μm, these cells would flow at off-centered positions along the channel height H in z-direction and their projections on the x- y-observation plane do not exhibit the characteristic slipper shape. Thus, such cells would be classified as others, similar to the so-called sheared croissant class, observed in previous studies.^12^ This effect emphasizes the importance of selecting not only the appropriate channel dimensions in terms of channel height and width but also their relative aspect ratio with respect to the optical access of the channel. We did not use channels with dimensions smaller than 8μm in this study because decreasing the channel size further results in the formation of tightly squeezed, symmetric bullet-like shapes that do not transition into other asymmetric or off-centered slipper-like shapes, as demonstrated for microcapillaries with diameter 4.7-6.6μm corresponding to χ=1.7-1.2.^8,9^
The RBC shape is intrinsically linked to its equilibrium position in the microchannel. Based on the 2D projection of the cells in flow and the optical axes of the used setup (see blue arrows in Fig. 2), we evaluate the RBC equilibrium y-position along the channel width W. It is assessed using the probability density distributions (pdf) of the absolute value of the normalized y-coordinate |y/W| as a function of the velocity (Fig. 4).
FIG. 4. Equilibrium cell position across the channel width in the different (a) W=10μm and H=8μm, (b) W=10μm and H=10μm, (c) W=10μm and H=15μm, (d) W=15μm and H=8μm, (e) W=15μm and H=10μm, and (f) W=15μm and H=15μm. Probability density distributions (pdfs) of the absolute values of the RBC y-position normalized by the channel width |y/W|. Data are shown for five velocities for the different geometries based on the data of Fig. 3.
For the smallest channel cross section with W=10μm and H=8μm [Fig. 4(a)], symmetric croissants flow at a centered position |y/W|≈0 at low velocities. Increasing the cell velocity results in a reduction of the peak at the central position and an off-centered peak emerges at |y/W|≈0.22. For the square channel [ W=H=10μm, Fig. 4(b)], most cells also flow along the channel central axis |y/W|≈0 even at high velocities exceeding 10mm/s. This pronounced centered peak in the pdfs persists when keeping the channel width fixed and increasing the height further to H=15μm, yet the pdf broadens to more off-centered cell positions at the same time [Fig. 4(c)].
In the rectangular channel with W=15μm and H=8μm [Fig. 4(d)], most cells flow at a central position |y/W|≈0 at low velocities. At velocities above 5mm/s, the pdfs exhibit only an off-centered peak at |y/W|≈0.22. Increasing the height to H=10μm [Fig. 4(e)] results in the occurrence of two pronounced peaks in the distributions at v>8mm/s, with one at the centerline and an off-centered one. Upon further increase of the channel height to 15μm [Fig. 4(f)], we still observe a central peak at low velocities up to v=20mm/s. However, RBCs do not flow at a preferred position across the channel width at higher velocities, as indicated by the broad distribution at v>25mm/s.
Our analysis of the cell’s y-position distribution aligns with the observed phase diagrams. When a significant number of slipper-shaped RBCs is observed [Figs. 3(a), 3(d), and 3(e)], a pronounced off-centered peak at |y/W|≈0.22 appears in the pdfs [Figs. 4(a), 4(d), and 4(e)]. This observation is in good agreement with previous work.^15^ Interestingly, while we find only one off-centered peak at v>5mm/s in the channel with W=15μm and H=8μm [Fig. 4(d)], two peaks at |y/W|≈0 and |y/W|≈0.22 appear for the slightly deeper channel with W=15μm and H=10μm [Fig. 4(e)]. Such pdfs with two pronounced peaks at elevated velocities have been reported before^15,31^ and are also found in the most confined channel [Fig. 4(a)] due to the emergence of both off-centered slippers and central flowing croissants and other shapes. Note that both channels [Figs. 4(a) and 4(e)] have similar channel aspect ratios AR≈1.4 (see Table I). However, the singularly peaked distribution in the channel with W=15μm and H=8μm [Fig. 4(d)] at v>5mm/s suggests that the significant fraction of other-shaped cells [Fig. 3(d)] also flow at an off-centered position, similar to the concurrently appearing slippers. These results highlight that both the RBC shape and the cell’s position in the channel cross section provide valuable information for evaluating RBC flow properties relevant to multiple clinical applications.^24^
In micro-confined conditions, the RBC shape critically depends on the flow rate in the channel.^11,47^ Here, we study how the shape of stable croissant and slipper-like RBCs is affected by the flow conditions in the different microchannels [Fig. 5(a)]. To characterize this effect, we assess the cell’s elongation index EI as a function of the shear rate in the channel.
FIG. 5. Elongation of stable croissant-like and slipper-like RBC shapes in the different channels. (a) Representative images of croissants and slippers in different channels and at various shear rates γ˙. (Scale bars, 5μm.) (b) and (c) Elongation index EI=(a−b)/(a+b) in the different microfluidic channels for croissant-like and slipper-like RBCs, respectively. Data are shown as a function of the shear rate (bottom x axes) and capillary number Ca (top x axes). The length a of the RBC in the flow direction and the diameter b along the channel width are identified with a bounding box around the RBC shape, schematically shown in the top right image in (a). Error bars correspond to averaging EI data of cells at the same applied pressure drop.
At low shear rates, croissant-shaped RBCs exhibit a broad shape with a shorter length in flow direction than across the channel width, hence, EI<0 [Fig. 5(b)]. As the shear rate increases up to γ˙≈5×103s−1, there is a substantial rise of the elongation index, independent of the channel cross section. Concurrently, the increase in cell length in the flow direction leads to a transition to positive elongation indices at γ˙≈5×103s−1. Beyond 8×103s−1, cell elongation seems to saturate, eventually reaching a plateau at EI≈0.1.
For slipper-like RBCs, we also observed a gradual increase in the elongation index with the channel shear rate [Fig. 5(b)]. However, due to the initially elongated slipper shape, EI is always positive. Between 1×103 and 3×103s−1, slipper-like cells exhibit an increase in EI, followed by a saturation at EI≈0.3-0.4 beyond 5×103s−1.
Assessing the elongation and stretching of RBCs by the fluid shear stress in microcapillaries has emerged as an important technology measuring RBC deformability.^25,27^ Previous work revealed differences in the dynamic deformation behavior between control RBCs and artificially stiffened RBCs, as well as RBCs of diabetes patients in microcapillary flow.^15,48,49^ In our study, the observed increase of cell elongation for stable croissants is in good agreement with previous investigations on single RBCs.^11,50^ This is attributed to the two “tails” of the croissant-shaped RBC that seem to become longer and pointed [see the top row in Fig. 5(a)] as the velocity and shear stress increase. While the cell elongation of croissants does not seem to saturate completely even at high shear rates [Fig. 5(b)], EI of slipper-shaped RBCs clearly plateaus above γ˙≈5×103s−1, similar to the previously reported high-velocity limit of cell elongation in confined microcapillaries.^9^ Future work will be required to evaluate how the observed dynamic deformation behavior of both stable RBC configurations is connected with the intrinsic cell properties, such as cytosol viscosity and membrane elasticity.
In our experimental investigations, we primarily focus on stable RBC shapes in microcapillary flow. In an effort to explain the suppression of these stable shapes and the occurrence of other shapes (see Fig. 3), we assess temporal shape changes of the RBC shape within the region of interest with a length of roughly 300μm in the middle of the microfluidic chip under steady flow conditions. We consider an RBC to have an unstable shape when it rotates, tumbles, or shows any other dynamic shape transitions that lead to fluctuations Δy=y(t)−y¯ of the cell’s temporal y-position y(t) from its mean y-position y¯ larger than 5% during passage within the region of interest.
Croissant and slipper-shaped RBCs exhibit stable shape configurations that do not change significantly while passing the microfluidic channel [two top rows in Fig. 6(a)]. For the other cell shape, we find a stable category [third row in Fig. 6(a)], as well as an unstable category that exhibits rich temporal dynamics [bottom row in Fig. 6(a)]. While the cell position along the channel width of the three stable shape classes does not change significantly during flow, we observe strong fluctuations of the cell’s center of mass position for the category of unstable other cell shapes [representatively shown in Fig. 6(b)]. The fraction of unstable other shapes is larger in the microfluidic channels with a large cross section, i.e., 10×15, 15×10, and 15×15μm2 [Fig. 6(c)].
FIG. 6. Stable and unstable other RBC shapes. (a) Representative image sequences of a croissant, a slipper, a stable other, and an unstable other cell shape in a channel with W=10μm and H=8μm. (Scale bars, 5μm.) Time intervals between consecutive images are displayed in the second image of each sequence. (b) Deviations Δy=y(t)−y¯ in the lateral movement over time of the cell’s center of mass position of the four representative cells shown in (a). (c) Fraction of unstable other cells with respect to the total number of other shapes for the different channels used in this study.
Although our experimental field of view only covers approximately 300μm along the channel flow direction, recent cell-tracking measurements demonstrated that the croissant and slipper shapes are indeed stable.^13^ Once the cell achieves its shape after entering the microfluidic channel, it does not change under steady flow conditions. In our study, we examine the single RBC flow in the middle of the microfluidic chip at L/2=20mm. At this position, cells already reached their equilibrium y-position and final shape, and transient effects induced by the fluid inlet can be neglected.^51^ Note that the length over which tank-treading slipper-shaped RBCs in such confined channels oscillate and periodically change their y-position is usually much longer than the region of interest used in this study.^13^ Hence, tank-treading slippers are classified as a stable RBC configuration.
In contrast, the dynamic behavior of unsteady other cell shapes does not allow the RBC to reach a steady y-position [see representatively Fig. 6(b)]. These cells predominantly appear in less confined channels with χ≤0.65 [see Fig. 6(c); 10×15, 15×10, and 15×15μm2], which is in good agreement with numerical simulations that predict a transition from steady to dynamic shapes upon increasing the channel dimensions.^16,20,21^ In such large channels, RBCs can flow at off-centered streamlines and tumble and rotate more easily driven by the parabolic velocity profile than in the strongly confined channels. While we only classify other shapes as stable and unstable, previous work has revealed a plethora of dynamic RBC states, including tumbling, rolling, swinging, snaking, and tank-treading motions.^28,29,52^ Such dynamic RBC morphologies have received much attention as they can affect blood shear thinning behavior under microcirculatory flow conditions.^1,53^ Since approximately 20% of others exhibit unstable shapes in the channels that generate the highest fraction of stable croissants and slippers, a future research objective will be to investigate the exact nature of these shapes. For example, if tumbling trilobes or multilobes can be reliably detected within the field of measurements, this class could be integrated into the shape phase diagram. Under pathological conditions when RBC deformability or membrane mechanical properties are impaired, and various unstable states emerge, such refined shape phase diagrams represent a central role in evaluating RBC flow behavior and could be used in diagnostic applications.
In this study, we perform an experimental characterization of microscale RBC flow behavior, elucidating consistent RBC shape state diagrams and illustrating the complexity of RBC shapes across various confined microchannels. Based on the assessment of RBC shapes in microcapillary flow, our findings emphasize the significance of strongly confined channels, specifically with a height below 10μm and χ≈0.9 along with a rectangular cross section, to generate the characteristic croissant peak and slipper plateau at low and high velocities, respectively. Despite the potential advantages of increased throughput associated with larger channel dimensions, our results show that the number of unsteady other cell shapes that are not uniquely identifiable increases when enlarging the channel dimensions. The implications of our study are pivotal for future microfluidic designs employing RBC shape classification approaches for diagnostic purposes for specific pathologies or as quality assessments of stored blood units. In combination with recent studies highlighting the influence of intrinsic cell properties on microscale blood flow, our work underscores the critical role of confinement effects and shear rate as external conditions impacting the RBC shape, advancing our understanding of in vitro single RBC flow.
This work was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) Project No. 349558021 (WA 1336/13-1 and RE 5025/1-2). M.N. acknowledges support and funding from the “Deutsch-Französische-Hochschule” (DFH) DFDK CDFA-01-14 “Living fluids.” A.D. acknowledges funding from the Young Investigator Grant of Saarland University.
The authors have no conflicts to disclose.
Ethics approval for experiments reported in the submitted manuscript on animal or human subjects was granted. The study involving human blood was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of the Aerztekammer des Saarlandes (approval number 51/18). Written informed consent was obtained from all subjects involved in the study.
Mohammed Nouaman: Conceptualization (equal); Formal analysis (lead); Investigation (lead); Methodology (lead); Writing – review & editing (equal). Alexis Darras: Conceptualization (equal); Formal analysis (equal); Methodology (equal); Supervision (equal); Writing – review & editing (equal). Christian Wagner: Funding acquisition (lead); Resources (lead); Supervision (equal); Writing – review & editing (equal). Steffen M. Recktenwald: Conceptualization (equal); Formal analysis (equal); Funding acquisition (equal); Methodology (equal); Supervision (equal); Visualization (lead); Writing – original draft (lead); Writing – review & editing (equal).
The data that support the findings of this study are available from the corresponding author upon reasonable request.
The data that support the findings of this study are available from the corresponding author upon reasonable request.