Kazuma
Yoda
a,
Yoshiyasu
Ichikawa
bc and
Masahiro
Motosuke
*bc
aDepartment of Mechanical Engineering, Graduate School of Engineering, Tokyo University of Science, Japan
bDepartment of Mechanical Engineering, Faculty of Engineering, Tokyo University of Science, Japan. E-mail: mot@rs.tus.ac.jp
cWater Frontier Research Center, Research Institute for Science and Technology, Tokyo University of Science, Japan
First published on 23rd October 2023
This paper presents the concept of a newly developed high-throughput measurement device for determining the dielectric properties of cancer cells. The proposed continuous-flow electrorotation (cROT) device can induce electrorotation (ROT) with vertical rotation using two sets of interdigitated electrodes on the top and bottom substrates to torque the cells. In the developed device, multiple rotating cells flowing in a microchannel are aligned between electrodes using dielectrophoresis. This allows for the measurement of the rotational behavior of the cells with continuous flow, resulting in a significant improvement in throughput compared to the conventional ROT devices reported previously. The dielectric properties, permittivity of the cell membrane and conductivity of the cell cytoplasm, of HeLa cells obtained by simultaneous measurements using the developed cROT device were 9.13 ± 1.02 and 0.93 ± 0.10 S m−1, respectively. Moreover, the measurement throughput was successfully increased to 2700 cells per h using the cROT technique.
The dielectric properties of cancer cells reflect cancer types and their status of drug resistance.1,3,12,13 Alazzam et al. reported the microfluidic separation of circulating tumor cells (CTCs) from blood utilizing these characteristics.19 Compared to normal blood cells, CTCs show different dielectric properties, which have different polarities of dielectrophoresis (DEP) forces in a non-uniform electric field. In addition to CTCs, exosomes and circulating tumor DNA (ctDNA) have recently been used as biomarkers, showing the possibility of separation and collection by DEP.20–22 To realize applications for advanced cancer diagnosis and monitoring using electrical properties, it is necessary to acquire the properties of cancer cells and establish their database in advance. Thus, methods for determining the electrical properties of cancer cells are in high demand.23–25
Electrorotation (ROT) has been proposed as a non-invasive technique to determine the dielectric properties of cells. ROT is a phenomenon in which cells rotate depending on their dielectric properties, and the permittivity and conductivity of individual cells can be determined from their rotational behavior.13,23 ROT can provide the dielectric properties of cancer cells without contact between the cell and the electrode; thus, it is not affected by the shape of the electrode or its contact resistance. As shown in Fig. 1(a), a rotating electric field was generated by applying four AC voltages with a phase difference of 90° using four electrodes, #1–4. Consequently, the cell placed at the center of the electrodes rotates owing to the rotating electric field. The angular velocity of the rotation varied depending on the frequency of the applied AC electric field. The type and state of the cells can be characterized by analyzing the frequency-dependent angular velocity. To determine the dielectric properties of a cell by ROT, individual cells must be placed at the center of each electrode when measuring a single cell or multiple cells simultaneously. Because negative DEP can act as a force on the cell in the direction away from the electrode, many previous studies exploited DEP with ROT to capture target cells.23,26–30 This method enabled the measurement of the stable rotation behavior of a cell because the DEP force levitates the cells and they rotate without contact with the electrode or substrate. However, some cells were not appropriately captured or the measurement was disrupted by other irruptive cells during analysis. Therefore, it is necessary to perform the capture operation for the measurement of a large number of cells.17–19 In addition, high voltage is required to provide sufficient torque for cell rotation during the capture of the cell by DEP.23,31
We developed a device that simplifies the capture of cancer cells and the measurement procedures of their dielectric properties, resulting in improved throughput for characterization. This technique is called continuous-flow ROT (cROT) because the rotational behavior of multiple cells can be measured continuously and simultaneously using a single device. The electrode locations for the cROT are schematically shown in Fig. 1(b). In the cROT, four sets of electrodes (#1–4) are arranged in the xz cross-section, as shown in Fig. 1(b), using two comb-shaped electrodes on the top and bottom of the device. This electrode structure can rotate the cells while maintaining the same torque between the electrodes because of the uniform electrode array on any plane in the xz cross-section. This structure allows the cells to flow in the y-direction in the microfluidic device while being aligned by DEP, with simultaneous measurement of their electrical properties by ROT. Additionally, cell manipulations other than pumping were not required. Therefore, the cROT has the potential to realize high-throughput measurements of the electrical properties of cancer cells. We believe that this would improve the throughput of conventional ROT devices.29–31 Moreover, cROT has another advantage over conventional ROT in terms of the simpler device configuration. Conventional ROT devices for multiple cell analysis need an electrode array with at least three layers: wiring, insulation, and electrodes.31 In contrast, the cROT device can be fabricated with minimal alignment by combining single-layer electrode substrates with the same top and bottom patterns. Hence, the cROT device is simpler than that of conventional ROT.
This study was aimed at demonstrating the validity of the cROT concept for the simultaneous measurement of dielectric properties, particularly the permittivity of the cell membrane and the conductivity of the cell cytoplasm, with improved throughput. First, the rotational behavior of the cells was obtained without flow in the device to evaluate the performance of the cROT electrode structure (Fig. 1(b)) and measure the dielectric properties. Here, we introduced a procedure to acquire the properties based on image processing and confirmed that cROT could obtain the dielectric properties of HeLa cells similarly to conventional ROT. Second, the rotational behavior of cells continuously flowing between the electrodes in a microfluidic device was measured, and it was confirmed that the throughput of the cROT was significantly higher than that of the conventional ROT.
The remainder of this manuscript is organized as follows. Section 2 presents a conceptual study of the device, including the rotation mechanism and numerical results on the design of the device geometry. Section 3 presents the experimental setup. Section 4 presents an evaluation of the measurement performance of the cROT with and without flow conditions and discusses the throughput of cROT-based measurements of the dielectric properties of flowing cancer cells. Finally, section 5 summarizes the conclusions of this study.
The electric field induced by the cROT was evaluated by numerical simulations to investigate the mechanism of cell rotation between the electrodes. COMSOL Multiphysics® (v6.0) was utilized for conducting the three-dimensional numerical simulations in this study. The calculation domain was set at 160 μm × 160 μm × 50 μm (width (x) × depth (y) × height (z)) (see Fig. S1† in the ESI† for details on the simulation of the calculation domain, governing equations, and boundary conditions). The width of each electrode on the substrate is 30 μm, and the gaps between them are 50 μm. The channel height is also 50 μm. AC voltages of 1.5 Vpp and a frequency of 200 kHz were applied to the four electrodes, which were phase-shifted by 90° with respect to each other to induce ROT.
Fig. 3(a) shows the time-series electric field in the xz cross-section acquired by numerical simulation. The contour plot indicates the magnitude of the electric field and the arrows indicate the direction. Here, we focus on the area surrounded by a red frame located between the electrodes, where a cell is to be placed, to investigate the time-dependent variation of the electric field (bottom of Fig. 3(a)). T indicates the oscillation period of the applied electric field, and the time interval t in each field corresponds to T/8. Fig. 3(b) shows the ratio of the x-component (Ex) to the z-component (Ez) of the electric field, Ex/Ez, at the cell location, which indicates the direction of the torque. To determine the value of Ex/Ez equivalent to the tangent value in the rotating frame of reference, a distance of 20 μm from the center point of the electrode was utilized. Fig. 3(a) and (b) show that a rotating electric field is generated between the four electrodes in the calculation domain. Therefore, the cells placed between the electrodes can rotate because of the torque generated by the rotating electric field.
The angular velocity Ω of the cell receiving the torque from the rotating electric field is given as follows,
(2.1) |
(2.2) |
(2.3) |
In the cROT, DEP was applied to cells simultaneously with ROT to place cancer cells between the electrodes. In DEP, when a dielectric particle is placed in an electric field, it polarizes and experiences an electrical force. DEP can be generated by applying 180°-phase-shifted AC voltages to the electrodes; thus, the same electrodes as those shown in Fig. 2 can be used. The DEP force acting on a particle FDEP, is expressed as follows:
FDEP = 2πεmR3Re[K*(ω)]∇|E|2 | (2.4) |
Fig. 3(c) shows the magnitude of the electric field and the direction of the electric field gradient when a negative DEP is induced. Here, 180°-phase-shifted AC voltages of 1.5 Vpp and a frequency of 20 kHz were applied to the electrodes of the same calculation model (Fig. S1†), to arrange the cells between the electrodes. As shown in Fig. 3(c), the electric field increases near the electrode and decreases further away from the electrode. Thus, the cells in the device were placed between the electrodes because the DEP force worked on the cells in the direction away from the electrodes, as indicated by the arrows in Fig. 3(c). Because the DEP force works even if flow occurs in the device, cells are transported from upstream to downstream, maintaining their positions between the electrodes when the flow direction is perpendicular to the xz-plane, i.e. along the y-direction. Therefore, the cROT can be achieved under these conditions.
To increase the ROT torque and DEP forces acting on the cancer cells, the electrode gap (channel height) was adjusted appropriately in the actual device. Numerical simulations were also used to investigate the effects of the device geometry on the ROT and DEP. Fig. 4 shows that the intensity of the electric field increases near the electrodes. Additionally, from eqn (2.1), the angular velocity of the cell increases with the field intensity. Since the typical diameter of the cancer cells to be measured in this study ranges from 12 to 20 μm, the electric field at the center of the electrode for a 20 μm-sized cell (the location corresponding to the dashed line shown in Fig. 4) was compared. Because cells come into contact with the top or bottom of the channel when the channel height is 20 μm, the minimum spacing of the electrode and channel was set to 30 μm. Electric field simulations were performed for electrode gaps and channel heights of 30, 40, and 50 μm. A voltage of 1.5 Vpp at a frequency of 200 kHz was applied. The calculation conditions, excluding the channel height, were the same as those used in the previous simulation. Fig. 4(b) indicates that the electric field at the floating position of the cell increases as the electrode gap decreases, and the electric field at a gap of 30 μm is approximately double that at 50 μm. When the channel height is 30 μm, the rotation speed of the cell is four times faster than that at a height of 50 μm, based on eqn (2.1). The smaller the electrode gap, the easier it is to measure the rotational movement of the cell. The magnitude of the DEP force required to align the cells between the electrodes also depends on the electric field, based on eqn (2.3). Therefore, a smaller electrode gap enables a stronger DEP force to work on the cells, and the cell is held more stably between the electrodes. From these results, both the electrode gap and channel height in this study were set to 30 μm.
Fig. 5 Observed image of rotating HeLa cells and their dielectric properties using the cROT electrode without flow. (a) Multiple HeLa cells were arranged, and they rotated between the electrodes while maintaining their positions. The measurement positions (i)–(vi) are surrounded by four electrodes (#1–4; these numbers correspond to those shown in Fig. 1(b)), creating a uniform rotating electric field. Depending on the electrodes and line positions, the rotation direction of the cell is different between odd positions ((i), (iii), (v)) and even positions ((ii), (iv), (vi)). The applied voltage signal is a superposition of 1.5 Vpp with 100 kHz for ROT and 5.0 Vpp with 20 kHz for DEP. (b) Examples of rotation analysis of each cell. The rotation speed was determined by image processing. The image processing used template matching to compare the image of frame 0 (template image) with the images of other successive frames (compared images). The comparison results were output as a similarity, confirming the periodicity corresponding to the rotational behavior. Therefore, the rotation speed of the cells was calculated from the periodicity of the similarity. (c) ROT spectrum obtained from HeLa cells. Dielectric properties were calculated by curve fitting. Since the dielectric properties are a parameter included in the CM factor, the dielectric properties were estimated by the least-squares method for the imaginary part of the CM factor calculated from the measured angular velocity. The error bars in the figure correspond to the standard error (n ≈ 10). |
Subsequently, the angular velocity of each rotating cell was determined using image processing. Template matching utilizing OpenCV was employed to analyze the rotational speed. The first frame of each cell image was employed as the template image, and the successively obtained images in other frames were used for comparison. The similarity between the template and compared images was calculated as the normalized correlation coefficient, and the time variation of the similarity in the successive images for each cell was investigated to determine the rotational periodicity. The validity of the rotation analysis method was preliminary evaluated by artificial cell images with virtual rotation, and the accuracy of the angular velocity was less than 0.25%.
The time-series similarity among cell No. 2, 5, and 7 depicted in Fig. 5(a) was obtained, as shown in Fig. 5(b). During rotation, the similarity values show periodic behavior with decreasing and increasing values; when a cell rotates once, the similarity value approaches unity. Notably, the similarity is equal to one when the complete cell image matches the template and compared images, but in the actual case, the maximum value of the similarity after a rotation is smaller than one since the appearance is not exactly the same because of the slight change in the background. The period of similarity for the rotating cells (shown in Video S2 in the ESI†) is consistent with the cycle of cell rotation. Therefore, it was possible to determine the rotational speed of each cell using a periodicity analysis based on template matching.
The image processing indicated above was applied to all the observed rotating cells, and the angular velocity of the cells was obtained for each frequency of the applied voltage as shown in Fig. 5(c). The imaginary part of the CM function was calculated using eqn (2.1)- based on the angular velocities which were acquired from image processing. Then the dielectric properties, εmem and σcyt in eqn (2.2), were estimated using the least-squares method. In this calculation, we obtained two parameters, εmem and σcyt, among the four variables in eqn (2.3) because the other two parameters, σmem and εcyt, showed little change in the frequency range normally used in ROT.13 The fixed parameters in the least-squares fitting are shown in Table 1.30 The obtained values of εmem and σcyt were 7.34 ± 0.75 (±standard error, S.E.) and 1.07 ± 0.11 (±S.E.) S m−1, respectively. The dielectric properties of HeLa cells measured by ROT have been reported as εmem = 7.40–19.78 and σcyt = 0.36–1.25, with a relatively wide range of variation.13,30,35 In our previous study,34 the dielectric properties of HeLa cells measured using a ROT device based on the concept shown in Fig. 1(a) were 10.20 ± 1.74 (±S.E.) and 1.01 ± 0.13 (±S.E.) S m−1. The variation in these measurement values between the present results and those obtained in our previous study is within the variation of the reported results; therefore, the measurement results obtained in this study are considered reasonable. Therefore, it can be said that the dielectric properties of cancer cells could be measured using the electrode structure in our developed device based on the cROT concept (shown in Fig. 1(b)).
Because multiple cells placed between the electrodes can be measured simultaneously in our device, these measurement schemes, which rely on the cROT electrode structure, are expected to improve the throughput of dielectric property measurements. Furthermore, in this experiment, cells, except for aggregated cells (approximately 10% of the observed cells), rotated between the electrodes; thus, the cell capture rate was achieved at approximately 90% of the infused cells in the device. This implies the potential of cROT to improve the low capture rate, which is an issue in conventional ROTs.29,31
Under these experimental conditions, the dielectric properties of HeLa cells were measured using the cROT. The device and applied voltage conditions are the same as in section 4.1. Fig. 6(a) shows an example of the images acquired for rotating cells while flowing at t = 0 s and 10 s at a frequency of 100 kHz. Flowing cells were automatically detected using in-house image processing software and are surrounded by white boxes, as shown in Fig. 6(a) and Video S3.† The cells flowed between electrodes at a uniform speed. Here, the detected objects (cells surrounded by white boxes in Fig. 6(a)) with a width of 20 μm or larger were determined as aggregated cells based on the size of the detected cells, and they were excluded from the analysis. Thus, more than 30 vertically rotating cells were detected in the observed area, and it was possible to analyze them simultaneously.
Fig. 6 Dielectric property measurement of HeLa cells by cROT under flow conditions. (a) Multiple HeLa cells aligned between electrodes by DEP showed rotational behavior by ROT. Since a uniform electric field was formed between the electrodes in the flow direction, the cells continued rotating along the flow direction between the electrodes. All observed cells flowed with almost the same velocity of 20 μm s−1. Since the length of the observation area is 748.8 μm, the cells to be measured are automatically replaced in ≈40 s. The applied voltage signal is a superposition of 1.5 Vpp with 800 kHz for ROT and 1.5 Vpp with 20 kHz for DEP. (b) The rotational behavior of HeLa cells measured by cROT was analyzed by image processing as in Fig. 5(b). The same periodic similarity as that shown in Fig. 5(b) is obtained under the flow conditions, and the angular velocity of the cell could be calculated. (c) ROT spectrum obtained by cROT. Dielectric properties were calculated by curve fitting using the angular velocity of the cells obtained by image processing. The error bars correspond to the standard error (n ≈ 30). |
The rotational behavior of the detected cells was analyzed using the image processing method described in section 4.1. Some detected cells could not be analyzed for angular velocity owing to low image quality or noise, and at most, 32 cells could be measured simultaneously. Fig. 6(b) shows the typical time-series similarity of the cells in this experiment. In addition, in the cROT with flow, a periodicity in similarity was observed. The angular velocities of the cells were then measured. The dielectric properties, εmem and σcyt, calculated from the ROT spectrum shown in Fig. 6(c) were 9.13 ± 1.02 (±S.E.) and 0.93 ± 0.10 (±S.E.) S m−1, respectively. The dielectric properties of the HeLa cells obtained by conventional ROT and cROT are summarized in Table 2. As described in section 4.1, this result is within a reasonable value considering the range of variation of previously reported properties.13,30,35 Therefore, it is possible to say that the cROT can properly measure the dielectric properties under continuous flow conditions. Additionally, εmem and σcyt measured in the cROT device with and without flow can be regarded as equivalent, and it is considered that the analysis for the flowing cells has no influence on the dielectric property measurements.
ROT34 | cROT (w/o flow) | cROT (w/ flow) | |
---|---|---|---|
Membrane permittivity (εmem) [−] | 10.20 ± 1.74 | 7.34 ± 0.75 | 9.13 ± 1.02 |
Cytoplasm conductivity (σcyt) [S m−1] | 1.01 ± 0.13 | 1.07 ± 0.11 | 0.93 ± 0.10 |
Finally, the measurable number of cells per unit time (cells per h) was evaluated as the throughput of dielectric property measurements. The throughput was calculated based on the number of cells that could be captured simultaneously and the measurement time. In addition, the throughput was compared with those obtained in previous studies. These values are summarized in Table 3. Here, the measurement time was the sum of the time required to obtain the ROT spectra and the time required to introduce and replace cells. In the conventional ROT, the process of capturing a new cell at the measurement location where the rotating electric field is induced must be repeated if the device can measure only a single cell.13,23,34,36 The time required to acquire the ROT spectra was approximately 3 min, and the cell replacement time was approximately 5 min; thus, the throughput was calculated to be 12–20 cells per h.
Electrodes [−] | Cells [−] | Measurement [h] | Throughput [cells per h] | |
---|---|---|---|---|
Conventional ROT13,23,34,36 | 1 | 1 | 0.05–0.08 | 12–20 |
Tsuchiya et al.31 | 36 | 5 | 0.05–0.08 | 60–100 |
Suzuki et al.37 | 225 | 70 | 0.25 | 420 |
Keim et al.30 | 39 | 39 | 0.025 | 600 |
cROT (this study) | 13 | 32 | 0.011 | 2700 |
Some studies have reported the simultaneous measurement of the dielectric properties of cells with improved throughput using an array of electrodes. Tsuchiya et al.31 fabricated a 6 × 6 array of pillar-shaped electrodes that could simultaneously measure the rotational behaviors of multiple beads. They indicated that the measurements could be conducted simultaneously at five locations. Suzuki et al.37 fabricated 225 sets of electrode arrays using microband electrodes placed above and below the microwells; the cells were captured in microwells in 10–15 min, and 70 cells were measured simultaneously. Keim et al.30 developed a DEP actuator that could independently capture and release individual cells. This technique enabled 100% cell capture at 39 measurement positions and reduced the cell replacement time to 12 s. These improvements in cell manipulation have increased the measurement throughput to 600 cells per h.
Compared to the previous ROT studies presented above, the cROT developed in this study allows the cell replacement time to be ignored by introducing a continuous flow for cell replacement and measurement. This is because the cells automatically flowed into and out of the observation area. The cells passed through the 748.8 μm length of the observation area at a flow velocity of 20 μm s−1 so that the ROT spectra could be obtained during this 40 s. For this case, the measurement throughput of cROT is calculated to be 2700 cells per h. Although the throughput itself is not as high as the impedance flow cytometry (IFC),38,39 ROT measurement under continuous flow in the proposed cROT has the potential to realize a significantly higher throughput of dielectric property measurements in cancer cells.
When the rotational behavior of cancer cells under continuous flow is measured, cells that settle and adhere to the substrate surface cannot be measured. In addition, cells that adhere to the substrate aggregate with the flowing cells, clogging the channel. In our study, the number of cells adhering to the substrate was approximately 10–20% of all flowing cells. The adhesion rate of cells can be improved by appropriate surface treatments or coatings, and the throughput can be further improved by increasing the number of measurable cells. Furthermore, the number of cells that can be observed is limited by current imaging systems, particularly by the resolution and frame rate of the camera. Further improvement in the measurement throughput is expected by capturing a wider domain with higher spatial resolution and frame rates.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3lc00301a |
This journal is © The Royal Society of Chemistry 2023 |