Next Article in Journal
Effects of Block Copolymer Terminal Groups on Toughening Epoxy-Based Composites: Microstructures and Toughening Mechanisms
Next Article in Special Issue
Squeezing Droplet Formation in a Flow-Focusing Micro Cross-Junction
Previous Article in Journal
Design of a Femtosecond Laser Percussion Drilling Process for Ni-Based Superalloys Based on Machine Learning and the Genetic Algorithm
Previous Article in Special Issue
Particle Counting Methods Based on Microfluidic Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deformability-Based Isolation of Circulating Tumor Cells in Spiral Microchannels

Department of Mechanical Engineering, Shahrekord University, Shahrekord 88186-34141, Iran
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(11), 2111; https://doi.org/10.3390/mi14112111
Submission received: 29 September 2023 / Revised: 8 November 2023 / Accepted: 17 November 2023 / Published: 17 November 2023
(This article belongs to the Special Issue Microfluidic Chips and Microdevices for Biomedical Engineering)

Abstract

:
The isolation of circulating tumor cells (CTCs) and their analysis are crucial for the preliminary identification of invasive cancer. One of the effective properties that can be utilized to isolate CTCs is their deformability. In this paper, inertial-based spiral microchannels with various numbers of loops are employed to sort deformable CTCs using the finite element method (FEM) and an arbitrary Lagrangian–Eulerian (ALE) approach. The influences of cell deformability, cell size, number of loops, and channel depth on the hydrodynamic behavior of CTCs are discussed. The results demonstrate that the trajectory of cells is affected by the above factors when passing through the spiral channel. This approach can be utilized for sorting and isolating label-free deformable biological cells at large scales in clinical systems.

1. Introduction

The isolation of rare cells, including CTCs, is a complicated process because they are notably unprotected from unintended losses [1]. The isolation of single cells from heterogeneous samples is pivotal in medical and biological applications [2,3]. On the other hand, the deformability or mechanical stiffness of individual cells affects the cellular phenotype in biological assessments [4]. It has been demonstrated that the deformability of CTCs is excessively diminished in comparison with normal tissues and correlated with metastatic potential [5].
Two major criteria for isolating CTCs from blood cells are their size and deformability. The diameter of epithelial CTCs varies from 14 μm to 26 μm and the diameter of white blood cells (WBCs) changes from 8 μm to 20 μm. Various passive and active label-free techniques have been utilized to separate cells, such as hydrodynamic methods that isolate cells based on their size, filtration methods that separate cells due to their size and deformability, and inertial focusing by employing cell intrinsic characteristics, i.e., size and stiffness [6,7,8].
The influence of the inertial lift force, F L , was first considered by Di Carlo et al. [9]. They assessed the trajectory of particles in straight and curved microchannels. It has been demonstrated that the channel Reynolds number, R e c h = U m a x D h / ν , and the particle Reynolds number, R e p = R e c h D 2 / D h 2 , where U m a x = 1.5 U a v g , D h , ν , and D are the maximum fluid velocity, hydraulic diameter, kinematic viscosity, and particle diameter, respectively, describe the impact of the inertial force on the motion of microparticles [10]. Di Carlo et al. [9] examined the influence of R e p , D / D h , the density of microparticles, etc., and suggested the optimal ranges of these factors to augment the impact of F L on particle focusing. Di Carlo et al. [11] also revealed that the optimal amount of flow rate for their proposed microfluidic device is 0.9 μL/min to isolate platelets from whole blood. Ozkumur et al. [12] carried out experimental examinations to sort various kinds of CTCs, including prostate, breast, pancreas, melanoma, and lung, utilizing inertial force. Their microfluidics device involved three outlets. Red blood cells (RBCs), platelets, and other blood components exited from the first outlet. Then, CTCs labeled with magnetic beads and WBCs were sorted based on their size and deformability in a curved channel. They were then separated by applying a magnetic field and exited from the second and third outlets. Oakey et al. [13] introduced a staged microchannel consisting of straight and curved sections for 0.2 R e p 5.9 . The curved section was the same as the one utilized by Ozkumur et al. [12]. The effectiveness of their proposed device was augmented by enhancing the fluid throughput and particle concentration. Kim et al. [14] carried out a parametric study to evaluate the impact of inertial focusing in a curved microchannel on the separation of CTCs. The effective parameters included the height and width of the microchannel as well as the curvature radius. Using sheath-less spiral microchannels, they isolated human breast cancer cells from WBCs with a purity of 86.76%. Abdulla et al. [15] employed a cascade inertial focusing device consisting of two spiral microchannels and a zigzag one to separate lung and breast cancer cells from WBCs and reached cell viability > 95%.
Almost all numerical simulations have not examined the deformability of CTCs during their continuous separation because the stiffness can be considered by utilizing specific numerical techniques, such as fluid-solid interaction (FSI), for individual cells. Hence, numerical works have been limited to size-based continuous separation of CTCs. For example, Ozbey et al. [16] employed COMSOL multiphysics software 6.1 to assess the focusing treatment of four types of CTCs utilizing a curvilinear geometry. It was found that 11–22 μm CTCs can be isolated from 8–17 μm CTCs when the curvature angle is 280° and R e c h = 121. Shiriny and Bayareh [17] utilized a single-loop spiral microfluidic device to isolate breast and epithelial cervical cancer cells and reported that the separation efficiency may reach 100% for 90 R e c h 110 . Recently, Pakhira et al. [18] carried out three-dimensional simulations to isolate various kinds of CTCs, such as breast, prostate, and lung cells, using inertial-based single- and double-loop spiral microchannels for 55.9 R e c h 65 . It was revealed that the single-loop device is more efficient than the double-loop one.
Even though several experimental studies have proposed deformability-based inertial devices [19,20], a few numerical works have assessed deformability influences on the focusing behavior of CTCs in different microfluidic geometries. For instance, Quek et al. [21] utilized the immersed boundary scheme to separate deformable cells in deterministic lateral displacement (DLD) devices consisting of regular arrays of microposts. They reported three types of cell trajectories, including zigzag, laterally displaced, and dispersive. The third type was observed for large and rigid particles when their stiffness was > 500 MPa. Additionally, Mohammadali et al. [22] employed the FSI approach to predict the impact of deformability on the trajectory of WBCs and two types of CTCs in a DLD device. They demonstrated that deformability influence on the trajectory of soft cells is higher than stiff ones, especially at low amounts of R e c h . It was also revealed that I-shaped micropillars lead to larger lateral displacement of CTCs. Guzniczak et al. [23] assessed the hydrodynamic behavior of cells in a spiral channel and demonstrated that equal-sized cells have different equilibrium positions according to their deformability. It was reported that the final position of soft and stiff cells is in the vicinity of the outer and inner walls, respectively.
In the present paper, the FSI approach is employed for the first time to evaluate the impact of deformability on the trajectory of WBCs and two kinds of CTCs in spiral microchannels. The hydrodynamic behavior of soft and stiff cells of various sizes is studied when the number of loops is 0.5, 0.75, 1, and 2 (Figure 1). This approach can be utilized for sorting and isolating label-free deformable biological cells at large scales in clinical systems. The rest of the paper is organized as follows. The deformable particle theory and mathematical model are presented in Section 2. Section 3 provides boundary and initial conditions. Grid study and validation are provided in Section 4. Section 5 discusses the numerical results, and concluding remarks are presented in Section 6.

2. Theoretical Background

2.1. Fluid Flow

As Figure 1 demonstrates, spiral microfluidic devices with various numbers of loops are considered. A spherical WBC or CTC with different radiuses is located at a distance of 200 μm from the channel inlet and a distance of 250 μm from the channel walls. Since the blood is diluted for real applications, the number of RBCs per volume is reduced considerably; hence, it can be assumed to be a Newtonian fluid. The governing equations, i.e., continuity and Navier–Stokes equations, for Newtonian incompressible fluid are expressed as follows [17]:
· u f = 0
ρ f u f t + ρ f u f · u f = · p f I + μ f u f + u f T + F f
where u f and p f represent fluid velocity vector and pressure, respectively. The density and dynamic viscosity of fluid are ρ f = 1000 k g / m 3 and μ f = 0.001 P a . s , respectively. Since volume forces, such as gravitational force, are not applied in the present simulation, F f = 0 .
For rigid particles, spiral microchannels work based on the balance between the inertial lift force ( F L ) and Dean drag force ( F D ) [24,25]:
F L = ρ U m a x 2 r P 4 D h 2 C L
F D 5.4 × 10 4 π μ D e 1.63 r P
where r P and C L are the particle radius and lift coefficient, respectively. D e = R e c h D h / 2 R represents Dean number, where R is the radius of the channel curvature. Thus, F s = F L + F D for rigid cells.
For deformable particles moving in straight [8] and curved [23] microchannels, an additional lift force, named deformability-induced lift force F D L , affects their migration. Guzniczak et al. [23] examined the contribution of F D L to the migration mechanism of soft particles in a microchannel with a cross-section of 360 × 60 μm2. They concluded that the impact of F D L can be neglected when R e c h is low. It was revealed that the lateral equilibrium position of deformable cells becomes closer to the outer wall of a curved microchannel due to the influence of F D L .
In the present work, D h = 2 W H / W + H , where W and H denote the microchannel width and depth, respectively. The ratio of H / W determines the number of lateral focusing positions [8]. For continuous separation of particles, when H / W > 2 , the particle overlap possibility is diminished.

2.2. Solid Mechanics

Many microfluidic devices employ inertial forces to focus or sort cells using passive methods. The trajectory of deformable cells can be described as follows utilizing the governing equations of linear elastodynamics [24]:
ρ s 2 u s 2 t = F s + · σ s
σ s = C ε s
ε s = 1 2 [ u s T + u s + u s T ( u s ) ]
Here, u s , σ s , ε s , ρ s , F s , and C are the solid displacement field, Cauchy stress tensor, strain tensor, solid density, volumetric force per unit volume, and stiffness matrix, respectively. The stiffness matrix is defined as follows:
C = E s 1 + R s 1 2 R s 1 R s R s R s 0 0 0 R s 1 R s R s 0 0 0 R s R s 1 R s 0 0 0 0 0 0 1 2 R s 0 0 0 0 0 0 1 2 R s 0 0 0 0 0 0 1 2 R s
Here, E s and R s represent Young’s modulus and the Poisson ratio of the solid, respectively. Equations (3)–(6) can be solved when the amounts of ρ s , F s , E s , and R s are known. The density of CTCs and WBCs is assumed to be 1050 k g / m 3 . The properties of cells are presented in Table 1.

2.3. FSI

Fluid flow and solid mechanics are coupled by the FSI to capture the interface between two phases. This approach describes the impact of fluid flow on the solid boundary and the influence of solid displacement on fluid velocity:
f s = n · p f I + μ f u f + u f T
u w = u s t ,   u w = u f
where f s , n , and u w are the total force exerted on the solid boundary, the normal vector, and the change rate of u s .
The present system is solved using the FEM and ALE method by employing COMSOL Multiphysics software. The ALE approach has been verified for deformable particles [27]. This technique combines the interface between the fluid and the solid.
The governing equations are discretized by the FEM using simple functions. A Lagrangian approach and a material frame are utilized to formulate the solid mechanics. Since the solid domain has a constant number of nodes, the following relation can be expressed between f s and F s [28]:
F s = f s d v d V
Here,   d v and d V denote the scale factors for the grid elements of the spatial and material frames, respectively.

3. Boundary and Initial Conditions

The fluid flows in the spiral microchannel due to the pressure gradient between its inlet and outlet. The fluid flow is assumed to be fully developed and laminar and its parabolic velocity distribution is described as follows, where Y is the material frame coordinate along the microchannel inlet [28].
U i n = U a v g 6 W Y Y W 2
The boundary condition at the outlet is the Dirichlet condition:
p f = 0 ,   μ f u f + u f T n = 0
Additionally, the no-slip boundary condition is also used for the microchannel walls:
u f = 0
Also, the initial amounts of u f , p f , and u s are assumed to be zero for simplicity.

4. Grid Study and Validation

Due to the motion of cells, the grid is deformed during the simulations, resulting in a reduction in the grid quality when the cell moves through the microchannel. Thus, the computational domain is re-meshed each time to provide a grid quality below 0.7. Four grid resolutions are utilized to perform the grid independence test. Figure 2a illustrates the velocity distribution for a cross-sectional area located at a 500 μm distance from the inlet. Four grids with 37,695, 72,148, 162,374, and 249,860 elements are considered in this figure, indicating that the grids with 162,374 and 249,860 elements result in the same velocity profile. Thus, the grid with 162,374 elements can be selected for further simulations.
The present numerical results are compared and verified using the experimental data and numerical results of Xu et al. [28]. For numerical simulations, they utilized the FEM to evaluate the hydrodynamic trapping of deformable particles in a microfluidic device consisting of microposts. Figure 2b illustrates the variations of microparticle velocity as a function of time, indicating a good agreement between the present results and experimental data.

5. Results and Discussion

In this section, the isolating process of single cells is examined based on their size and deformability. The influences of cell stiffness, cell size, the number of loops, and microchannel depth on cell trajectory are examined.

5.1. Effect of Cell Deformability

The most important reason for the separation of CTCs is the formation of Dean flow in spiral channels, leading to the creation of a secondary flow in the channel cross-section. In this section, CTC1, CTC2, and WBCs with the same diameter of 16 μm enter the spiral microchannel along with the diluted blood flow when R e c h = 1 . Cells are subjected to lift force in the direction perpendicular to the flow and drag force in the flow direction. Due to the small size of the cells, the gravity force is neglected. As mentioned previously, deformable cells also experience an additional lift force F D L , which pushes them to an equilibrium position closer to the outer wall of the channel, while rigid particles are concentrated closer to the inner wall. For a constant cell size, deformability has a great effect on the trajectory of CTCs. Figure 3 demonstrates that the cells tend to move toward the outer wall as deformability is augmented. Softer cells are more deformed due to pressure and shear stress (Figure 4). In other words, CTCs show higher stress than WBCs due to their larger size and various mechanical characteristics. Additionally, the stress distribution indicates a core-like pattern in CTCs and WBCs that can be utilized to design inertial-based microfluidic separation devices.

5.2. Effect of Cell Size

In a spiral microchannel, the deformable cell equilibrium position is determined due to the balance between F L , F D , and F D L . The cells with different sizes experience various hydrodynamic behaviors due to the formation of Dean flow. As shown in Figure 5, larger cells are directed toward the inner channel wall. Smaller cells follow the direction of F D , i.e., the secondary flow, and move toward the outer channel wall, which is consistent with the experimental results of Guan et al. [29]. The Dean vortices in the vicinity of the walls are weak and cannot re-entrain smaller cells. This difference in the focusing behavior of large and small cells provides an opportunity for their separation. In other words, larger cells enter the Dean flow due to the influence of lift and drag forces; however, smaller ones do not remain in Dean vortices due to not overcoming the centrifugal force and move in the outer part of the microchannel.
When the length of the spiral channel and its radius are not significantly large, bending effects are ignored. Therefore, the major force that affects the cells is the shear force, which is caused by the motion of fluid on the boundary of the cell. In this case, the difference in the amount of stress applied to the cells depends on their distance from the channel wall and their distance from each other. Cells that are closer to the channel wall are subjected to more stress. For certain kinds of cells, the one with a larger diameter is in more contact with the fluid layer, resulting in more applied shear stress (Figure 6).

5.3. Effect of the Number of Loops

In this section, the effect of the number of spiral loops on the hydrodynamic behavior of the deformable cells is evaluated. CTC1 with a diameter of 16 μm along with diluted blood enter the microchannel with various number of loops when R e c h = 1 . The number of loops has a great effect on the final position of cells. By augmenting the number of loops, the length of the path that the cells travel is enhanced, reducing their velocity due to an increment in the drag force. Therefore, to choose the optimal number of loops for the isolation of CTCs, the balance between different forces acting on the cells should be considered. Figure 7 depicts the deviation of CTC1 from the channel centerline of different geometries. It is demonstrated that the amount of deviation is 15.92 μm and 16.3 μm toward the outer wall for the 0.5- and 0.75-loop spiral microchannels, respectively. Additionally, the amount of deviation for 1- and 2-loop microchannels is 1.8 and 71.81 μm toward the inner wall, respectively, indicating that the equilibrium position of cells is strongly affected by the lengthening of the microchannel. The increase in the number of loops, i.e., the increment in the radius of curvature, leads to a higher Dean number D e = R e c h D h / 2 R , which causes the cell to deviate more toward the inner wall. For the case of the 2-loop microchannel, the deviation of CTCs toward the inner wall is much higher than the 0.5-, 0.75, and 1-loop microchannels due to the larger radius of curvature and formation of stronger Dean vortices. Hence, it can be concluded that the 2-loop spiral microchannel can be selected as the optimal case for isolating single cells.
The amount of stress exerted on the cells by the fluid is illustrated in Figure 8. Since the increase in the number of loops affects the pressure and velocity of the fluid, the change in the shape of the cells in microchannels is significant. By augmenting the number of loops, the length of the microchannel increases, leading to a decrease in the fluid velocity. Thus, the centrifugal force and the magnitude of the shear stress are reduced. Figure 8 also depicts the contour of the velocity distribution at the outlet of the microchannels.

5.4. Effect of the Microchannel Depth

The Dean flow in a spiral channel leads to a shorter distance of cells to reach their equilibrium position compared to the straight channel [23]. In this section, the effect of the depth of the spiral microchannels is discussed for CTC1 with a diameter of 22 μm when R e c h = 1 . Four values of the channel depth, h = 100, 200, 300, and 500 μm, are considered. Figure 9 shows that the channel depth (channel aspect ratio) has a significant effect on the isolation of CTCs. One of the important parameters for focusing microparticles/cells in microchannels is the blockage ratio, which is defined as the ratio of the cell diameter to the channel hydraulic diameter, β = 2 r p / D h . It was known that β > 0.07 is required to focus microparticles/cells [30]. At β < 0.07 , the influence of F L on the focusing process is diminished. The amount of β is 0.132, 0.077, 0.058, and 0.044 for h = 100, 200, 300, and 500 μm, respectively, when r p = 11 μm. It should be pointed out that F D and F D L are two other significant forces that affect the trajectory of cells. In this section, the deformability of CTCs is maintained constant. Since the amount of R e c h and the channel width are kept constant, the fluid velocity is reduced by enhancing the channel depth (or D h ). Thus, the centrifugal force caused by the blood flow is decreased with h, reducing the amount of F D . In other words, there is a competition between F L and F D during the motion of CTCs in spiral microchannels. Figure 9 demonstrates that the CTC deviation is augmented with h. This is due to the smaller effect of centrifugal force, i.e.,   F D , on CTC trajectories compared to F L .

6. Conclusions

Spiral microchannels with various loops are designed and analyzed to isolate deformable CTCs utilizing an ALE and FEM numerically. The trajectory of CTCs is controlled by their deformability and size, number of loops, and channel depth. The present simulations suggest that an inertial-based microfluidic device can be employed to sort deformable cells in clinical systems. It is demonstrated that the cells move toward the outer wall of the microchannel as deformability is augmented. Additionally, larger and smaller deformable cells are directed toward the inner and outer channel walls, respectively, due to the balance between F L , F D , and F D L . The results reveal that the 2-loop spiral microchannel can be considered the optimal geometry. Furthermore, the microchannel aspect ratio affects the trajectory and deviation of deformable cells. Since diluted blood samples are employed for practical applications, the combination of micromixers [31,32,33] and microfluidic separation devices is recommended for the isolation of CTCs from blood cells. Additionally, inertial focusing of non-spherical deformable cells can be considered to examine the influence of cell shapes on their isolation mechanisms.

Author Contributions

R.M.: methodology, data curation, and software. M.B.: the preparation of the original manuscript draft, manuscript editing, project administration, and final manuscript review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adams, A.A.; Okagbare, P.I.; Feng, J.; Hupert, M.L.; Patterson, D.; Göttert, J.; McCarley, R.L.; Nikitopoulos, D.; Murphy, M.C.; Soper, S.A. Highly Efficient Circulating Tumor Cell Isolation from Whole Blood and Label-Free Enumeration Using Polymer-Based Microfluidics with an Integrated Conductivity Sensor. J. Am. Chem. Soc. 2008, 130, 8633–8641. [Google Scholar] [CrossRef] [PubMed]
  2. Shiriny, A.; Bayareh, M. On magnetophoretic separation of blood cells using Halbach array of magnets. Meccanica 2020, 55, 1903–1916. [Google Scholar] [CrossRef]
  3. Paterlini-Brechot, P.; Benali, N.L. Circulating tumor cells (CTC) detection: Clinical impact and future directions. Cancer Lett. 2007, 253, 180–204. [Google Scholar] [CrossRef] [PubMed]
  4. Guck, J.; Schinkinger, S.; Lincoln, B.; Wottawah, F.; Ebert, S.; Romeyke, M.; Lenz, D.; Erickson, H.M.; Ananthakrishnan, R.; Mitchell, D.; et al. Optical Deformability as an Inherent Cell Marker for Testing Malignant Transformation and Metastatic Competence. Biophys. J. 2005, 88, 3689–3698. [Google Scholar] [CrossRef] [PubMed]
  5. Cross, S.E.; Jin, Y.-S.; Rao, J.; Gimzewski, J.K. Nanomechanical analysis of cells from cancer patients. Nat. Nanotechnol. 2007, 2, 780–783. [Google Scholar] [CrossRef] [PubMed]
  6. Hochmuth, R.M. Micropipette aspiration of living cells. J. Biomech. 2000, 33, 15–22. [Google Scholar] [CrossRef] [PubMed]
  7. Rosenbluth, M.J.; Lam, W.A.; Fletcher, D.A. Force Microscopy of Nonadherent Cells: A Comparison of Leukemia Cell Deformability. Biophys. J. 2006, 90, 2994–3003. [Google Scholar] [CrossRef] [PubMed]
  8. Hur, S.C.; Henderson-MacLennan, N.K.; McCabe, E.R.B.; Di Carlo, D. Deformability-based cell classification and enrichment using inertial microfluidics. Lab A Chip 2011, 11, 912. [Google Scholar] [CrossRef]
  9. Di Carlo, D.; Irimia, D.; Tompkins, R.G.; Toner, M. Continuous inertial focusing, ordering, and separation of particles in microchannels. Proc. Natl. Acad. Sci. USA 2007, 104, 18892–18897. [Google Scholar] [CrossRef]
  10. Asmolov, E.S. The inertial lift on a spherical particle in a plane Poiseuille flow at large channel Reynolds number. J. Fluid Mech. 1999, 381, 63–87. [Google Scholar] [CrossRef]
  11. Di Carlo, D.; Edd, J.F.; Irimia, D.; Tompkins, R.G.; Toner, M. Equilibrium Separation and Filtration of Particles Using Differential Inertial Focusing. Anal. Chem. 2008, 80, 2204–2211. [Google Scholar] [CrossRef] [PubMed]
  12. Ozkumur, E.; Shah, A.M.; Ciciliano, J.C.; Emmink, B.L.; Miyamoto, D.T.; Brachtel, E.; Yu, M.; Chen, P.-I.; Morgan, B.; Trautwein, J.; et al. Inertial Focusing for Tumor Antigen-Dependent and -Independent Sorting of Rare Circulating Tumor Cells. Sci. Transl. Med. 2013, 5, 179ra47. [Google Scholar] [CrossRef] [PubMed]
  13. Oakey, J.; Applegate, R.W.; Arellano, E.; Carlo, D.D.; Graves, S.W.; Toner, M. Particle Focusing in Staged Inertial Microfluidic Devices for Flow Cytometry. Anal. Chem. 2010, 82, 3862–3867. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, T.H.; Yoon, H.J.; Stella, P.; Nagrath, S. Cascaded spiral microfluidic device for deterministic and high purity continuous separation of circulating tumor cells. Biomicrofluidics 2014, 8, 064117. [Google Scholar] [CrossRef] [PubMed]
  15. Abdulla, A.; Liu, W.; Gholamipour-Shirazi, A.; Sun, J.; Ding, X. High-Throughput Isolation of Circulating Tumor Cells Using Cascaded Inertial Focusing Microfluidic Channel. Anal. Chem. 2018, 90, 4397–4405. [Google Scholar] [CrossRef]
  16. Ozbey, A.; Karimzadehkhouei, M.; Kocaturk, N.M.; Bilir, S.E.; Kutlu, O.; Gozuacik, D.; Kosar, A. Inertial focusing of cancer cell lines in curvilinear microchannels. Micro Nano Eng. 2019, 2, 53–63. [Google Scholar] [CrossRef]
  17. Shiriny, A.; Bayareh, M. Inertial focusing of CTCs in a novel spiral microchannel. Chem. Eng. Sci. 2020, 229, 116102. [Google Scholar] [CrossRef]
  18. Writtick Pakhira, R.K.; Ibrahimi, K.M. Distinct separation of multiple CTCs using inertial focusing phenomena utilizing single-looped spiral microfluidic lab-on-chip. Chem. Eng. Sci. 2023, 275, 118724. [Google Scholar] [CrossRef]
  19. Liu, Z.; Chen, R.; Li, Y.; Liu, J.; Wang, P.; Xia, X.; Qin, L. Integrated Microfluidic Chip for Efficient Isolation and Deformability Analysis of Circulating Tumor Cells. Adv. Biosyst. 2018, 1800200. [Google Scholar] [CrossRef]
  20. McFaul, S.M.; Lin, B.K.; Ma, H. Cell separation based on size and deformability using microfluidic funnel ratchets. Lab A Chip 2012, 12, 2369. [Google Scholar] [CrossRef]
  21. Quek, R.; Le, D.V.; Chiam, K.-H. Separation of deformable particles in deterministic lateral displacement devices. Phys. Rev. E 2011, 83, 056301. [Google Scholar] [CrossRef] [PubMed]
  22. Mohammadali, R.; Bayareh, M.; Ahmadi Nadooshan, A. Numerical investigation on the effects of cell deformability and DLD microfluidic device geometric parameters on the isolation of circulating tumor cells. Iran. J. Chem. Chem. Eng. 2023. [Google Scholar] [CrossRef]
  23. Guzniczak, E.; Otto, O.; Whyte, G.; Willoughby, N.A.; Jimenez, M.; Bridle, H.L. Deformability-induced lift force in spiral microchannels for cell separation. Lab A Chip 2020, 20, 614–625. [Google Scholar] [CrossRef]
  24. Bower, A.F. Applied Mechanics of Solids; CRC Press: Boca Raton, FL, USA, 2009. [Google Scholar]
  25. Nathamgari, S.S.P.; Dong, B.; Zhou, F.; Kang, W.; Giraldo-Vela, J.P.; McGuire, T.; McNaughton, R.; Sun, C.; Kessler, J.A.; Espinosa, H.D. Isolating single cells in a neurosphere assay using inertial microfluidics. Lab A Chip 2015, 15, 4591–4597. [Google Scholar] [CrossRef] [PubMed]
  26. Lee, J.E.; Moon, S.; Park, T. A New Approach for Separation of Circulating Tumor Cells (CTCs) in the Microfluidics to Enhance the Capture Purity with COMSOL Simulation. Master’s Thesis, DGIST, Daegu, Republic of Korea, 2015. [Google Scholar]
  27. Ai, Y.; Mauroy, B.; Sharma, A.; Qian, S. Electrokinetic motion of a deformable particle: Dielectrophoretic effect. Electrophoresis 2011, 32, 2282–2291. [Google Scholar] [CrossRef] [PubMed]
  28. Xu, X.; Zhenyu, L.; Arye, N. Finite element simulations of hydrodynamic trapping in microfluidic particle-trap array systems. Biomicrofluidics 2013, 7, 054108. [Google Scholar] [CrossRef] [PubMed]
  29. Guan, G.; Wu, L.; Bhagat, A.A.; Li, Z.; Chen, P.C.Y.; Chao, S.; Ong, C.J.; Han, J. Spiral microchannel with rectangular and trapezoidal cross-sections for size based particle separation. Sci. Rep. 2013, 3, 1475. [Google Scholar] [CrossRef]
  30. Bhagat, A.A.S.; Kuntaegowdanahalli, S.S.; Papautsky, I. Continuous particle separation in spiral microchannels using dean flows and differential migration. Lab A Chip 2008, 8, 1906. [Google Scholar] [CrossRef]
  31. Yin, B.; Yue, W.; Muhtasim Fuad Sohan, A.S.M.; Wan, X.; Zhou, T.; Shi, L.; Qian, C.; Lin, X. Construction of a desirable hyperbolic microfluidic chip for ultrasensitive determination of PCT based on chemiluminescence. J. Mater. Chem. B 2023, 11, 1978–1986. [Google Scholar] [CrossRef]
  32. Bahrami, D.; Bayareh, M. Impacts of channel wall twisting on the mixing enhancement of a novel spiral micromixer. Chem. Pap. 2022, 76, 465–476. [Google Scholar] [CrossRef]
  33. Usefian, A.; Bayareh, M. Numerical and experimental investigation of an efficient convergent–divergent micromixer. Meccanica 2020, 55, 1025–1035. [Google Scholar] [CrossRef]
Figure 1. Schematic of spiral microchannels with various numbers of loops (all the dimensions are in micrometers).
Figure 1. Schematic of spiral microchannels with various numbers of loops (all the dimensions are in micrometers).
Micromachines 14 02111 g001
Figure 2. (a) Velocity distribution in a cross-section located at a 500 μm distance from the inlet for various grid resolutions and (b) microparticle velocity versus time: a comparison between the present numerical results and the experimental data of Xu et al. [28].
Figure 2. (a) Velocity distribution in a cross-section located at a 500 μm distance from the inlet for various grid resolutions and (b) microparticle velocity versus time: a comparison between the present numerical results and the experimental data of Xu et al. [28].
Micromachines 14 02111 g002
Figure 3. Trajectory of CTC1, CTC2, and WBCs in a 0.5-loop spiral microchannel when R e c h = 1 .
Figure 3. Trajectory of CTC1, CTC2, and WBCs in a 0.5-loop spiral microchannel when R e c h = 1 .
Micromachines 14 02111 g003
Figure 4. The stress contour inside CTC1, CTC2, and WBCs when they reach the end of the 0.5-loop spiral microchannel.
Figure 4. The stress contour inside CTC1, CTC2, and WBCs when they reach the end of the 0.5-loop spiral microchannel.
Micromachines 14 02111 g004
Figure 5. Trajectory of CTC1 of various sizes in a 0.5-loop spiral microchannel when R e c h = 1 .
Figure 5. Trajectory of CTC1 of various sizes in a 0.5-loop spiral microchannel when R e c h = 1 .
Micromachines 14 02111 g005
Figure 6. The stress contour inside CTC1 when it reaches the end of the 0.5-loop spiral microchannel when R e c h = 1 : (a) r P = 6   μ m , (b) r P = 8   μ m , and (c) r P = 12   μ m .
Figure 6. The stress contour inside CTC1 when it reaches the end of the 0.5-loop spiral microchannel when R e c h = 1 : (a) r P = 6   μ m , (b) r P = 8   μ m , and (c) r P = 12   μ m .
Micromachines 14 02111 g006
Figure 7. The deviation of CTC1 at the outlet of spiral microchannels with various numbers of loops with respect to the channel centerline when R e c h = 1 . The centerline axis at the outlet passes through (X, Y) = (0, 8900), (3962.03, 5150), (−1000, 0), and (−2000, 0) for the 0.5-, 1-, 1.5-, and 2-loop spiral microchannels, respectively. The deviation of CTC1 is calculated utilizing their position and the coordinate of the centerline axis at the microchannel outlet.
Figure 7. The deviation of CTC1 at the outlet of spiral microchannels with various numbers of loops with respect to the channel centerline when R e c h = 1 . The centerline axis at the outlet passes through (X, Y) = (0, 8900), (3962.03, 5150), (−1000, 0), and (−2000, 0) for the 0.5-, 1-, 1.5-, and 2-loop spiral microchannels, respectively. The deviation of CTC1 is calculated utilizing their position and the coordinate of the centerline axis at the microchannel outlet.
Micromachines 14 02111 g007
Figure 8. Velocity and stress contours resulting from fluid effects on the CTC1 at the end of the microchannels when R e c h = 1 : (a) the number of loops is 0.5, (b) the number of loops is 0.75, (c) the number of loops is 1, and (d) the number of loops is 2.
Figure 8. Velocity and stress contours resulting from fluid effects on the CTC1 at the end of the microchannels when R e c h = 1 : (a) the number of loops is 0.5, (b) the number of loops is 0.75, (c) the number of loops is 1, and (d) the number of loops is 2.
Micromachines 14 02111 g008
Figure 9. The trajectory of CTC1 in a 0.5-loop spiral microchannel for different amounts of channel depth when R e c h = 1 .
Figure 9. The trajectory of CTC1 in a 0.5-loop spiral microchannel for different amounts of channel depth when R e c h = 1 .
Micromachines 14 02111 g009
Table 1. Properties of CTCs and WBCs [26].
Table 1. Properties of CTCs and WBCs [26].
R s E s , kPa r P , μm
CTC10.4100011
CTC210−818
WBC10−90.16
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mohammadali, R.; Bayareh, M. Deformability-Based Isolation of Circulating Tumor Cells in Spiral Microchannels. Micromachines 2023, 14, 2111. https://doi.org/10.3390/mi14112111

AMA Style

Mohammadali R, Bayareh M. Deformability-Based Isolation of Circulating Tumor Cells in Spiral Microchannels. Micromachines. 2023; 14(11):2111. https://doi.org/10.3390/mi14112111

Chicago/Turabian Style

Mohammadali, Roya, and Morteza Bayareh. 2023. "Deformability-Based Isolation of Circulating Tumor Cells in Spiral Microchannels" Micromachines 14, no. 11: 2111. https://doi.org/10.3390/mi14112111

APA Style

Mohammadali, R., & Bayareh, M. (2023). Deformability-Based Isolation of Circulating Tumor Cells in Spiral Microchannels. Micromachines, 14(11), 2111. https://doi.org/10.3390/mi14112111

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop