US20180292352A1 - Methods and systems for identifying a particle using dielectrophoresis - Google Patents
Methods and systems for identifying a particle using dielectrophoresis Download PDFInfo
- Publication number
- US20180292352A1 US20180292352A1 US16/007,658 US201816007658A US2018292352A1 US 20180292352 A1 US20180292352 A1 US 20180292352A1 US 201816007658 A US201816007658 A US 201816007658A US 2018292352 A1 US2018292352 A1 US 2018292352A1
- Authority
- US
- United States
- Prior art keywords
- mhz
- frequency
- microelectrode array
- particles
- sweep rate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 239000002245 particle Substances 0.000 title claims abstract description 135
- 238000004720 dielectrophoresis Methods 0.000 title description 62
- 238000000034 method Methods 0.000 title description 30
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 238000004891 communication Methods 0.000 claims abstract description 10
- 210000003743 erythrocyte Anatomy 0.000 claims description 88
- 238000010191 image analysis Methods 0.000 claims 5
- 230000003068 static effect Effects 0.000 description 103
- 230000004044 response Effects 0.000 description 91
- 239000011324 bead Substances 0.000 description 78
- 239000004793 Polystyrene Substances 0.000 description 59
- 229920002223 polystyrene Polymers 0.000 description 57
- 230000010287 polarization Effects 0.000 description 37
- 230000006399 behavior Effects 0.000 description 33
- 238000005259 measurement Methods 0.000 description 26
- 210000004027 cell Anatomy 0.000 description 24
- 238000002474 experimental method Methods 0.000 description 24
- 210000004369 blood Anatomy 0.000 description 23
- 239000008280 blood Substances 0.000 description 23
- 238000012360 testing method Methods 0.000 description 18
- 230000001052 transient effect Effects 0.000 description 17
- 230000005684 electric field Effects 0.000 description 16
- 239000000243 solution Substances 0.000 description 16
- 238000004458 analytical method Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 9
- 238000010408 sweeping Methods 0.000 description 9
- 230000008859 change Effects 0.000 description 8
- 239000006185 dispersion Substances 0.000 description 7
- 230000001965 increasing effect Effects 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 7
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 6
- 238000013459 approach Methods 0.000 description 6
- 230000000007 visual effect Effects 0.000 description 6
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 5
- 230000001419 dependent effect Effects 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 4
- 229910052802 copper Inorganic materials 0.000 description 4
- 239000010949 copper Substances 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 239000008121 dextrose Substances 0.000 description 4
- 230000002441 reversible effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000013480 data collection Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000011521 glass Substances 0.000 description 3
- 239000012528 membrane Substances 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 239000011780 sodium chloride Substances 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000032258 transport Effects 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 210000000601 blood cell Anatomy 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012856 packing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 239000000725 suspension Substances 0.000 description 2
- 238000010865 video microscopy Methods 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 240000004808 Saccharomyces cerevisiae Species 0.000 description 1
- 241000700605 Viruses Species 0.000 description 1
- 150000001450 anions Chemical class 0.000 description 1
- 239000000427 antigen Substances 0.000 description 1
- 102000036639 antigens Human genes 0.000 description 1
- 108091007433 antigens Proteins 0.000 description 1
- WQZGKKKJIJFFOK-VFUOTHLCSA-N beta-D-glucose Chemical compound OC[C@H]1O[C@@H](O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-VFUOTHLCSA-N 0.000 description 1
- 239000012620 biological material Substances 0.000 description 1
- 239000012472 biological sample Substances 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 150000001768 cations Chemical class 0.000 description 1
- 230000003833 cell viability Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 210000003850 cellular structure Anatomy 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 210000000805 cytoplasm Anatomy 0.000 description 1
- 210000000172 cytosol Anatomy 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 239000004205 dimethyl polysiloxane Substances 0.000 description 1
- 235000013870 dimethyl polysiloxane Nutrition 0.000 description 1
- 230000003292 diminished effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000037427 ion transport Effects 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 210000000265 leukocyte Anatomy 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000001404 mediated effect Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 238000001000 micrograph Methods 0.000 description 1
- CXQXSVUQTKDNFP-UHFFFAOYSA-N octamethyltrisiloxane Chemical compound C[Si](C)(C)O[Si](C)(C)O[Si](C)(C)C CXQXSVUQTKDNFP-UHFFFAOYSA-N 0.000 description 1
- 210000003463 organelle Anatomy 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000004987 plasma desorption mass spectroscopy Methods 0.000 description 1
- 229920000435 poly(dimethylsiloxane) Polymers 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000012798 spherical particle Substances 0.000 description 1
- 210000000130 stem cell Anatomy 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/26—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating electrochemical variables; by using electrolysis or electrophoresis
- G01N27/416—Systems
- G01N27/447—Systems using electrophoresis
- G01N27/44756—Apparatus specially adapted therefor
- G01N27/44791—Microapparatus
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B03—SEPARATION OF SOLID MATERIALS USING LIQUIDS OR USING PNEUMATIC TABLES OR JIGS; MAGNETIC OR ELECTROSTATIC SEPARATION OF SOLID MATERIALS FROM SOLID MATERIALS OR FLUIDS; SEPARATION BY HIGH-VOLTAGE ELECTRIC FIELDS
- B03C—MAGNETIC OR ELECTROSTATIC SEPARATION OF SOLID MATERIALS FROM SOLID MATERIALS OR FLUIDS; SEPARATION BY HIGH-VOLTAGE ELECTRIC FIELDS
- B03C5/00—Separating dispersed particles from liquids by electrostatic effect
- B03C5/005—Dielectrophoresis, i.e. dielectric particles migrating towards the region of highest field strength
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/1031—Investigating individual particles by measuring electrical or magnetic effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/483—Physical analysis of biological material
- G01N33/487—Physical analysis of biological material of liquid biological material
- G01N33/48707—Physical analysis of biological material of liquid biological material by electrical means
-
- G01N2015/0073—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/01—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials specially adapted for biological cells, e.g. blood cells
- G01N2015/012—Red blood cells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1006—Investigating individual particles for cytology
Definitions
- the present invention was conceived while performing work under CBET 0644538, CBET 1041338, and IIP 1340126, each of which has been awarded by the National Science Foundation. The government has certain rights in the invention.
- the present invention relates to identification of particles based on dielectrophoretic responses.
- a system for identifying a particle includes a microfluidic device; a microelectrode array including a plurality of electrodes, the microelectrode array disposed within the microfluidic device; a plurality of particles suspended in a solution and delivered to the microelectrode array using the microfluidic device; a signal generator operatively coupled to the microelectrode array; a particle detector adjacent to the microelectrode array; and a controller in operative communication with the signal generator and the particle detector.
- the controller is configured to apply an oscillating voltage signal to the microelectrode array between a low frequency and a high frequency at a sweep rate, wherein the sweep rate is no more than a maximum sweep rate, and determine a distribution of the plurality of particles relative to the microelectrode array at a plurality of frequency levels between the low frequency and the high frequency.
- a method of identifying a particle includes the steps of: placing a plurality of particles adjacent a microelectrode array, the microelectrode array including a plurality of electrodes; applying an oscillating voltage signal to the microelectrode array, the oscillating voltage signal varying between a low frequency and a high frequency at a sweep rate, wherein the sweep rate is no more than a maximum sweep rate; and determining a distribution of the plurality of particles relative to the microelectrode array at a plurality of frequency levels between the low frequency and the high frequency.
- FIG. 1 shows (a) Dielectric relaxation mechanism for PS beads showing cases when i) particle polarization occurs at a static frequency, ii) ⁇ MW is shorter than the slow frequency sweep rate ( ⁇ FS ) allowing the bead interface time to polarize in response to the non-uniform AC field, and iii) ⁇ MW is longer than the ⁇ FS for fast frequency sweep rates and the bead interface does not have time to fully polarize. (b) Schematic of the quadrapole electrodes micro patterned onto a glass slide, and (c) microdevice with PDMS fluidic layer bonded above the quadrapole electrodes silver-epoxied to copper leads.
- FIG. 2 shows (a) nDEP behavior of 6 ⁇ m PS beads suspended in E-pure H 2 O 2.5 ⁇ 10 ⁇ 4 S/m and 250V pp /cm 0.0063, 0.056 and 0.17 MHz/s sweep rates from 0.010 MHz to 1.0 MHz.
- PS bead assembly at slower frequency sweep rates track static frequency responses while 0.056 MHz/s illustrates transitional behavior and frequency sweeps above 0.17 MHz/s substantially lag the true static frequency DEP responses.
- FIG. 3 shows (a) 6 ⁇ m PS beads nDEP intensity profiles for 0.00080, 0.0063, and 0.056 MHz/s and static steady state (SS) measurements (black diamonds). Intensity analysis captures bead assembly to the quadrapole center with transient and SS regions. The slowest frequency sweep rate of 0.00080 MHz/s best predicts the static DEP responses.
- FIG. 4 shows (a) nDEP behavior of RBCs suspended in 0.1 S/m dextrose buffer and 250V/cm at 0.00080 MHz/s, 0.0063 MHz/s and 0.056 MHz/s sweep rates from 0.010 MHz to 0.50 MHz.
- FIG. 5 shows an image comparison of the nDEP and pDEP behavior of A+ red blood cells suspended in 0.10 S/m dextrose solution and 1000V pp /cm at 0.00080, 0.0016, and 0.0028 MHz/s sweep rates from 0.010 MHz to 1.0 MHz.
- (1) Denotes red blood cells nDEP behavior and (2) denotes red blood cell pDEP behavior.
- the red blood cells' DEP behavior at slower frequency sweep rates correlates well with the static frequency response (top row).
- FIG. 6 illustrates (a) nDEP and (b) pDEP intensity profiles for 0.00080, 0.0016, and 0.0028 MHz/s and static steady state measurements (solid circle).
- the intensity (arbitrary units) analysis captures the RBCs assembly toward the electrode center (nDEP) and near the electrodes (pDEP).
- nDEP electrode center
- pDEP near the electrodes
- 0.00080 MHz/s is the slowest
- 0.0016 MHz/s is the fastest sweep rate to best predict RBCs' static DEP response.
- FIG. 7 illustrates (a) RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (middle and bottom row).
- R1 and R2 arc two different repeats completed for this measurement.
- Each test was completed with A+ blood in 0.10 S/m at 1000V pp /cm.
- FIG. 8 illustrates (a) RBC static images (top row) compared to RBCs response using 0.0025 MHz/s (middle and bottom row).
- R1 and R2 are two different repeats completed for this measurement.
- Each test was completed with A+ blood in 0.10 S/m at 1000V pp /cm.
- FIG. 9 illustrates (a) RBC static images (top row) compared to RBCs response using 0.0026 MHz/s (middle and bottom row).
- R1 and R2 are two different repeats completed for this measurement.
- Each test was completed with A+ blood in 0.10 S/m at 1000V pp /cm.
- FIG. 10 illustrates (a) RBC static images (top row) compared to the response of RBCs using a sweep rate of 0.0028 MHz/s (middle and bottom row).
- R1 and R2 are two different repeats completed for this measurement.
- Panels (b) nDEP and (c) pDEP show intensity profiles for 0.0028 MHz/s with RBC static steady state measurements (circles). The images show poor agreement and intensity profiles illustrate the lack of reproducibility of agreement with static measurements.
- Each test was completed with A+ blood in 0.10 S/m at 1000V pp /cm.
- FIG. 11 illustrates (a) A ⁇ RBC static images (top row) compared to A ⁇ RBCs response using 0.0024 MHz/s (bottom row). (b) nDEP and (c) pDEP intensity profiles for 0.0024 MHz/s with RBC static steady state measurements (circles). The images and intensity profiles show good agreement with static measurements. Each test was completed in 0.10 S/m at 1000V pp /cm.
- FIG. 12 illustrates (a) O+ RBC static images (top row) compared to O+ RBCs response using 0.0024 MHz/s (bottom row). (b) nDEP and (c) pDEP intensity profiles for 0.0024 MHz/s with RBC static steady state measurements (circles). The images and intensity profiles show good agreement at 0.70 MHz and fair agreement at 0.80 MHz. Each test was completed in 0.10 S/m at 1000V pp /cm.
- FIG. 13 illustrates (a) B+ RBC static images (top row) compared to OB+ RBCs response using 0.0024 MHz/s (bottom row). (b) nDEP and (c) pDEP intensity profiles for 0.0024 MHz/s with RBC static steady state measurements (circles). The images and intensity profiles show good agreement over the tested frequency range 0.60-0.84MHz. Each test was completed in 0.10 S/m at 1000V pp /cm.
- FIG. 14 shows a plot of sweep rate (MHz/s) as a function of conductivity of a solution.
- the plot points are a sweep rate at which the results for a given conductivity were no longer accurate.
- the plotted curve defines a threshold sweep rate.
- sweep rates above the threshold sweep rate for a given concentration are too fast, and any sweep rate below the curve may be used accurately and reliably. This aids in determining the fastest sweep rates that may be used for a solution with a given concentration in order to decrease the overall time needed for the procedure.
- FIG. 15 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000V pp /cm. This illustrates the visual similarity between the static response and response at this sweep rate at the solution conductivity of 0.10 S/m.
- FIG. 16 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.50 S/m at 1000V pp /cm. This illustrates the visual similarity between the static response and response at this sweep rate at the solution conductivity of 0.50 S/m.
- FIG. 17 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 1.00 S/m at 1000V pp /cm. This illustrates the visual similarity between the static response and response at this sweep rate at the solution conductivity of 1.0 S/m.
- FIG. 18 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.25 S/m at 1000V pp /cm. This plot generally shows the 0.0027 MHz/s sweep rates near the threshold sweep rate for this concentration match the static frequency data well, but the 0.0028 MHz/s sweep rate begins to falter.
- FIG. 19 shows RBC static images (top row) compared to RBCs response using 0.0027 MHz/s and 0.0028 Mhz/s (middle and bottom row). Each test was completed with A+ blood in 0.25 S/m at 1000V pp /cm. This illustrates the visual similarity between the static response and the response at a sweep rate near the threshold sweep rate for the given conductivity of the solution.
- FIG. 20 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.5 S/m at 1000V pp /cm. This plot generally shows the 0.0029 MHz/s sweep rates near the threshold sweep rate for this concentration match the static frequency data well, but the 0.0030 MHz/s sweep rate begins to falter.
- FIG. 21 shows RBC static images (top row) compared to RBCs response using 0.0029 MHz/s and 0.0030 Mhz/s (middle and bottom row). Each test was completed with A+ blood in 0.50 S/m at 1000 V pp /cm. This illustrates the visual similarity between the static response and the response at a sweep rate near the threshold sweep rate for the given conductivity of the solution.
- FIG. 22 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 1.0 S/m at 1000V pp /cm. This plot generally shows the 0.0030 MHz/s sweep rates near the threshold sweep rate for this concentration match the static frequency data well, but the 0.0030 MHz/s sweep rate begins to falter.
- FIG. 23 shows RBC static images (top row) compared to RBCs response using 0.0030 MHz/s and 0.0031 MHz/s (middle and bottom row). Each test was completed with A+ blood in 1.00 S/m at 1000 V pp /cm. This illustrates the visual similarity between the static response and the response at a sweep rate near the threshold sweep rate for the given conductivity of the solution.
- FIG. 24 shows RBC static images (top row) compared to RBCs response using a sweep rate of ⁇ 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000 V pp /cm.
- FIG. 25 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.10 S/m at 1000 V pp /cm using the reverse sweep method as illustrated by the images in FIG. 24 .
- FIG. 26 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000 V pp /cm.
- FIG. 27 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.10 S/m at 1000 V pp /cm using the reverse sweep method as illustrated by the images in FIG. 26 .
- Alternating current (AC) dielectrophoresis (DEP) experiments for biological particles in microdevices have typically been applied at fixed frequencies. Reconstructing the DEP response curve from static frequency experiments is laborious, but is important for ascertaining differences in dielectric properties of biological particles.
- the disclosed systems and methods employ the novel concept of sweeping the frequency as a function of time to rapidly determine the DEP response curve from fewer experiments.
- Homogeneous 6.08 ⁇ m polystyrene (PS) beads were initially used as a model system to determine whether sweeping the frequency would be a viable method for generating DEP responses and then to identify an optimal sweep rate.
- the particles that are analyzed may include polystyrene beads (e.g. for testing purposes) or cells such as blood cells; in particular, different subtypes of red blood cells (e.g. A+, A ⁇ , O+, O ⁇ , etc.) may be distinguished based on the surface charge differences of the red blood cell subtypes (e.g. due to differing antigens on the cell surfaces).
- the particles may range in size from about 1-50 ⁇ m and should be detectable (e.g. through optical or electrical means) using the particle detector (e.g. an imaging system).
- a system may include a microfluidic device (which in some embodiments may include an enclosed microfluidic chamber) having a microelectrode array disposed in the fluid path of the microfluidic device.
- particles are delivered to the vicinity of the microelectrode array using the microfluidic device prior to data collection.
- data collection is performed in a “batch-wise” manner, i.e. a group of particles is delivered to the microelectrode array and fluid movement is then stopped before data collection begins so that particle movements that are observed are due to dielectrophoresis.
- the microelectrode array in one embodiment is a quadrapole arrangement as shown, for example, in FIGS. 1, 2, 4, and 5 .
- four electrodes are arranged in an “X” with a gap in the center ( FIG. 5 ).
- Other possible arrangements of electrodes include interdigitated electrodes, V-shaped electrodes, circular electrodes, and T-shaped electrodes.
- the electrodes are arranged so that oppositely charged electrodes are not parallel to one another, as this would create uniform fields whereas other, non-parallel geometries create non-uniform electric fields.
- the electrodes are arranged so that they create a spatially non-uniform field.
- the electrodes are attached to the bottom of the microfluidic device and the particles that are delivered to the device are initially distributed in the vicinity of the electrodes in a random arrangement (e.g. see upper left panel of FIG. 2 a ) before any electrical signal is applied.
- the microelectrode array may be made by depositing metal strips onto a glass slide with a cover having a microfluidic channel being bonded on top of the glass slide ( FIG. 1 c ).
- opposing pairs of electrodes may be electrically coupled using copper wires, as shown, with the leads (i.e. ground and “hot” AC signal) of a signal generator being connected to the copper wires.
- pairs of electrodes may be electrically coupled as shown in FIG. 1 c so that, upon stimulation, each electrode is 90° out of phase from the others to create a traveling wave signal.
- a signal generator is used to deliver an oscillating voltage to the electrodes.
- the voltage is applied at a peak-to-peak amplitude of 0.1V pp , 1V pp , 10V pp , 100V pp , or other suitable amplitude.
- the oscillating voltage is applied at frequencies of at least about 0.001 MHz, at least about 0.005 MHz, at least about 0.01 MHz, at least about 0.05 MHz, at least about 0.1 MHz, at least about 0.5 MHz, or at least about 1.0 MHz.
- the oscillating voltage is applied at frequencies of no more than about 10.0 MHz, no more than about 5.0 MHz, no more than about 2.0 MHz, no more than about 1.0 MHz, or no more than about 0.5 MHz.
- the frequency of the oscillating voltage is varied, for example from a low frequency to a high frequency, in order to collect data at a variety of different frequencies, a process referred to as “sweeping” the frequency.
- sweeping the frequency of the oscillating voltage
- comparable results are obtained when the frequency is swept from a high frequency to a low frequency. Sweeping the oscillating voltage using a continuously varying frequency permits a relatively large amount of data to be gathered in a short period of time.
- the rate at which the frequency sweep i.e. the “sweep rate,” may vary from 0.00001 MHz/s to 0.1 MHz/s.
- the continuously varying frequency may be approximated by a series of discrete, step-wise changes in frequency with an increment ranging from about 10 nHz to about 10 Hz.
- the optimum sweep rate may depend on conditions such as the conductivity of the solution in which the particles are suspended. The present inventors have found that when the frequency is swept above a certain maximum sweep rate the frequency changes too quickly, such that the particles do not have sufficient time to respond to the voltage signal at a given frequency before the signal changes to the next frequency. If the oscillating voltage signal is varied too quickly, i.e. above the maximum sweep rate, the observed particle movements and distributions will be inaccurate and could lead to an inconclusive or erroneous particle identification.
- the maximum sweep rate is no more than about 0.003 MHz/s, no more than about 0.0029 MHz/s, no more than about 0.0028 MHz/s, no more than about 0.0027 MHz/s, no more than about 0.0026 MHz/s, no more than about 0.0025 MHz/s, no more than about 0.0020 MHz/s, no more than about 0.0015 MHz/s, no more than about 0.0010 MHz/s, no more than about 0.0008 MHz/s, or no more than about 0.0005 MHz/s.
- a minimum sweep rate of at least about 0.00005 MHz/s, at least about 0.0001 MHz/s, at least about 0.00015 MHz/s, at least about 0.0002 MHz/s, at least about 0.0004 MHz/s, at least about 0.0005 MHz/s, at least about 0.00075 MHz/s, or at least about 0.0010 MHz/s may be used.
- Particle detection may be carried out with systems which are capable of identifying the spatial distributions of the particles with sufficient temporal (e.g. operating at 0.1-10 Hz) and spatial (e.g. capable of resolving 0.1 ⁇ m ⁇ 0.1 ⁇ m areas) resolution.
- images are collected at regular intervals (e.g. at video rates of 30 frames/sec or at other, slower rates such as 1 image or frame/sec) while sweeping the oscillating voltage.
- the images may be processed (for example several sequential video-rate frames may be averaged together) and the images or subregions thereof may be analyzed to characterize particle distribution and behavior at one or several frequencies.
- the analyses may include one or more of determining the particles' intensity profiles, transient responses, and velocities. Analyses may be conducted at one or more discrete locations within the images including along one or more lines, e.g. lines running between electrode tips. Analysis patterns for a sample which includes an unknown particle may be compared to patterns generated under equivalent conditions using known particles to determine the identity of the unknown particle.
- regions and patterns on the substrate are selected for analysis based on the areas that are expected to have the greatest change in electric field patterns and hence the greatest change in particle distribution at different frequencies, so as to provide the most information for distinguishing between particle types. Additional methods for analyzing DEP behavior of particles are disclosed in Salmanzadeh et al. (2012), Rozitsky et al. (2013), and An et al. (2014), each of which is incorporated herein by reference.
- a system for carrying out embodiments of the invention may include a controller for carrying out or more of the procedures disclosed herein.
- the controller may be in operative communication with one or more of the signal generator (which in turn is in communication with the microelectrode array) and the particle detector (e.g. imaging system).
- the controller may be a part of or in communication with a computer system.
- the computer system may be part of an existing computer system (e.g. on a smartphone, desktop computer, on-board computer, etc.) or may be implemented as a separate, standalone unit that is in local or remote communication with other components.
- the computer system(s) may be in wired or wireless communication with other systems through a combination of local and global networks including the Internet.
- Each computer system may include one or more input device, output device, storage medium, and processor (e.g. a microprocessor).
- Input devices may include a microphone, a keyboard, a computer mouse, a touch pad, a touch screen, a digital tablet, a track ball, and the like.
- Output devices include a cathode-ray tube (CRT) computer monitor, an LCD or LED computer monitor, touch screen, speaker, and the like.
- CTR cathode-ray tube
- the computer system may be organized into various modules including an acquisition module and an output module along with the controller, where the controller is in communication with the acquisition module and the output module.
- the various modules for acquiring and processing data and for returning a result may be implemented by a single computer system or the modules may be implemented by several computer systems which are in either local or remote communication with one another.
- Storage media include various types of local or remote memory devices such as a hard disk, RAM, flash memory, and other magnetic, optical, physical, or electronic memory devices.
- the processor may be any known computer processor for performing calculations and directing other functions for performing input, output, calculation, and display of data in accordance with the disclosed methods.
- implementation of the disclosed invention includes generating sets of instructions and data that are stored on one or more of the storage media and operated on by a controller, where the controller may be configured to implement various embodiments of the disclosed invention.
- the system may be in the form of a portable or handheld device which may be self-contained, including the microfluidic device, microelectrode array, the controller, and additional components such as a power supply and input/output capabilities.
- DEP enables phenotypically similar biological cells to be discriminated based on dielectric properties including the conductivity and permittivity of the membrane, cytoplasm, and other structurally relevant organelles.
- Cell components and structure contribute to a cell's signature dielectric dispersion.
- a particle's complex permittivity is frequency dependent and characterized by dielectric dispersion regions (y, ⁇ , and ⁇ , where w ⁇ ⁇ w ⁇ ⁇ w y ) specific to an applied frequency.
- Polarized particles can exhibit either positive dielectrophoresis (pDEP) or negative dielectrophoresis (nDEP) as a consequence of the frequency-dependent polarizability of the particle in the surrounding medium. Particles that exhibit pDEP move to high electric field regions and particles that exhibit nDEP move to low electric field regions. This motion up and down electric field gradients is described by the Clausius-Mossotti factor for spherical particles.
- Polarization is not an instantaneous event; charge transport into the interface takes a few microseconds in response to the electric field.
- Maxwell-Wagner dielectric relaxation is a physical phenomenon related to the transport delay of cation and anion alignment in and around the interface of the dielectric particle.
- particle polarization is driven by this conductive polarization.
- charges do not have enough time to move into and around the interface double layer, so particles experience polarization lag time as a result of the rapidly modulating field and do not reach maximum polarization.
- Maxwell-Wagner dielectric relaxation is characterized by a time constant, ⁇ MW , which is unique to each particle or cell due to the time constant's dependence on the cell dielectric properties.
- the time required for a particle to reach maximum polarization is given by Eq. (3) (see Morgan et al. (2003), p. 27; see also Grosse et al. (2010) and Mittal et al. (2008), each of which is incorporated herein by reference):
- ⁇ TW ( ⁇ p + 2 ⁇ ⁇ m ) ⁇ ⁇ 0 ⁇ p ⁇ 2 ⁇ ⁇ ⁇ m .
- Typical relaxation times for particle polarization vary from pico- to microseconds (see Morgan et al. (2003), p. 27; see also Grosse et al. (2010) and Mittal et al. (2008)), and the calculated ⁇ MW for polystyrene (PS) beads in our Epure H 2 O medium at 2.5 ⁇ 10 ⁇ 4 S/m is 3.5 las.
- a single AC cycle is on the order of 0.01 to 2 ⁇ s; the time delay in ion transport within a static frequency field of 0.010 to 2.0 MHz is such that 2 to 350 AC cycles must be completed before the particle experiences full polarization.
- FIG. 1 a highlights the Maxwell-Wagner particle polarization at the interface under static frequency as well as slow and fast frequency sweep rates.
- the particle experiences a constant frequency field such that the relaxation time is not a factor and the particle fully polarizes.
- a particle in a field with a slowly changing frequency sweep has a relaxation time, ⁇ ⁇ FS , that is less than ⁇ MW and thus the particle interface fully polarizes.
- a particle in a fast frequency sweep has a relaxation time, ⁇ ⁇ FS , that is larger than ⁇ MW and the particle interface does not have time to fully polarize in the field.
- PS beads are lossy dielectric particles treated as homogeneous spheres and are thus an idealized particle to examine new techniques, devices, or approaches to dielectrophoretic characterizations.
- Our system is easily able to discern pDEP and nDEP transitional behavior and adaptable to new frequency sweep techniques.
- microfluidic and dielectrophoretic (DEP) technologies enable a wide variety of particle polarizations with nonuniform electric fields on microchips to achieve particle manipulation, concentration, separations, and property-based identification.
- Particles can include bioparticles (e.g. DNA, viruses, or proteins) as well as cells (e.g. blood cells, cancer cells, stem cells, and yeast).
- bioparticles e.g. DNA, viruses, or proteins
- cells e.g. blood cells, cancer cells, stem cells, and yeast.
- the advantages to coupling DEP with microfluidics are small sample size (on the order of microliters), rapid analysis (approximately minutes to achieve results), minimal sample preparation, and minimal waste production.
- DEP experiments are completed at static, fixed frequencies such that maximum particle polarization can be achieved and measured.
- An optimal frequency sweep rate can be estimated by calculating the Maxwell-Wagner dielectric relaxation time for the particle/cell of interest, provided the cell's permittivity and conductivity is known. The frequency sweep rate chosen for the DEP study should then remain at frequencies below the inverse dielectric relaxation time (1/ ⁇ MW ) for 5-45 s (longer times spent below the threshold give better DEP predictions).
- the microdevice shown in FIG. 1 c was fabricated according to previously published microfabrication techniques (Grom et al. (2006), incorporated herein by reference), with the 100 ⁇ m wide electrodes spaced 200 ⁇ m apart aligned at 90° along the bottom of a 70 pm deep by 1000 ⁇ m wide microfluidic chamber as shown in FIG. 1 b.
- Polystyrene beads Cat No. PP-60-10, Spherotech, Lake Forest, Ill., USA
- 6.08 ⁇ m in diameter were centrifuged at 1300 min ⁇ 1 for 5 mins to separate the beads from the liquid.
- PS beads were resuspended in Epure H 2 O (18 M ⁇ or 2.5 ⁇ 10 ⁇ 4 S/m) at a 1:10 (bead to water) volumetric dilution ratio and vortexed.
- Microdevice was pre-rinsed with Epure H 2 O and Alconox precision cleaner (Cat No. 1104, Alconox Inc, White Plains, N.Y., USA) to prevent bead adhesion.
- PS beadEpure H 2 O suspension was pumped to the microchamber using a syringe. Time was allowed for inlet and outlet pressures to equalize and flow to stop.
- the function generator (Agilent 33250A, Agilent, Santa Clara, Calif., USA) was connected via copper leads to produce a 10V pp AC sine wave with frequencies ranging from 0.010-2.0 MHz at specific frequency sweep rates 0.00080, 0.0011, 0.0030, 0.0063, 0.013, 0.021, 0.028, 0.042, 0.056, 0.083, and 0.17 MHz/s. Frequency sweeps linearly increased the applied frequency as a function of time.
- Video recordings of PS beads DEP behaviors were analyzed with ImageJ (NIH, Bethesda, Md.) using intensity, transient slope, and velocity measurements. Since PS beads only exhibit nDEP over the frequency range of interest, intensity data acquisition from images was completed by drawing a rectangular box at the device center, I CTR , and background, I BK measured in a location with no PS beads present (See FIG. 2 a ). ImageJ Z Project function was used to average the pixel intensities in the specified boxed region.
- the transient response of the PS beads was extracted from the steady-state response via signal processing in which the delay and rise time were quantified.
- the PS bead delay time, t d was characterized as the time required for the intensity response to reach 50% of the final intensity response for the first time.
- the rise time, t T was determined as the time needed for the intensity response to reach 100% of the final intensity response for the first time (Ogata et al. (1978), pp. 517-518, incorporated herein by reference). This allowed the transient response to be segmented and a linear trend line was fit between t d and t T where t d ⁇ t T .
- a comparison of the transient slope for frequency sweep rates and static frequency measurements is given in FIG. 3 c .
- PS bead velocities were determined from the original video by tracking the x-, y-pixel position of individual PS beads from 0-50 s. PS bead located within 5 ⁇ m of electrode tips were selected to control for similar electric field gradients. This procedure was repeated for at least 10 beads in each specific frequency sweep rate and static experimental video.
- blood of the appropriate type e.g. O+, A+, etc.
- blood of the appropriate type e.g. O+, A+, etc.
- the packed RBCs were removed, then resuspended at 1:75 v:v in 0.10 S/m isotonic dextrose buffer doped with 0.75% BSA (Cat No. A7906, Sigma Aldrich, St. Louis, Mo., USA) to prevent cell/device adhesion.
- Video microscopy at 25 ⁇ and 1 fps was obtained with a Zeiss Axiovert Inverted Light Microscope (Zeiss, Germany). The video images were analyzed as described herein for the PS beads.
- nDEP responses were quantified via intensity analysis as described herein for all sweeps and all static frequency experiments.
- FIG. 2 b illustrates the frequency- (and time-) dependent intensity for the 0.0063 MHz/s sweep rate images shown in FIG. 2 a .
- This quantification of the PS bead nDEP response was correlated to total bead packing via the calibration shown in the inset.
- the 188-bead count at the center deviates slightly from the initial, field off, bead count of 245 because PS beads also move down the electric field gradient to regions outside of the image field of view.
- FIG. 3 a Normalized intensities, Eq. (4), were compiled in FIG. 3 a for SS (i.e. 45 seconds) static frequency nDEP responses and 0.00080, 0.0063, 0.056 MHz/s frequency sweep rate nDEP responses.
- the time for sweep responses to achieve the true nDEP static response decreases as the sweep rate decreases.
- FIG. 3 a inset shows that the slowest 0.00080 MHz/s sweep rate more quickly aligns closely with the static frequency responses.
- the confidence intervals around the transient 0.0063 MHz/s sweeps are smaller than for 0.17 MHz/s over much of the frequency range indicating greater reproducibility at slower sweep rates.
- Faster sweep rates either do not reach SS or have a lag before reaching SS (compare to FIG. 2 a ) suggesting the bead interface does not fully polarized and thus displays attenuated nDEP motion.
- the transient behavior was quantified for all static frequencies and frequency sweeps via a transient slope analysis as compiled in FIG. 3 c .
- Four static frequency measurements 0.010, 0.60, 1.0 and 2.0 MHz are shown compared to 0.00080, 0.0063, 0.028, 0.056, and 0.17 MHz/s frequency sweep rates.
- Static frequency transient slopes range between 0.023-0.095 and are within the 95% (p ⁇ 0.05) confidence intervals of 0.00080, 0.0063, and 0.028 MHz/s frequency sweep transient slopes.
- These slower sweep rates and 0.056 MHz/s differ at p ⁇ 0.001 from the fastest sweep rate of 0.17 MHz/s, which is also significantly different at p ⁇ 0.001 from the static measurements (except for 1.0 ⁇ 10 4 Hz with p ⁇ 0.01).
- the second timescale of interest is the frequency change per time or frequency sweep rate, which determines how many consecutive cycles a particle experiences a specific frequency.
- the PS beads experience a specific frequency for a large number of cycles and thus the beads have time to polarize because the timescale of the frequency change is slower than the dielectric relaxation time.
- a particle must experience a single frequency during the sweep for a minimum of 3.5 ⁇ s for maximum interfacial polarization to be achieved.
- the particle which its current DEP force has to overcome inertia and Stokes drag to achieve observable particle motion down the electric field gradient. At static frequencies, it takes roughly 5 s for maximum velocity to be attained (see FIG.
- a frequency sweep approach can be utilized to attain accurate DEP behavior of PS beads, provided the sweep rate is slower than conductivity mediated interfacial polarization timescale. This result is reliable over frequency ranges where particle polarization is dominated by the conduction of free charges from the media.
- the charges are moving around the PS beads through the particle-liquid interface inducing a dipole, which causes PS bead movement down the electric field gradient to the electrode center.
- the rate of movement of the charges varies which varies the rate of the dipole being induced, observed as dielectric relaxation.
- Each sweep rate has a unique dielectric relaxation time and our results are consistent with Maxwell-Wagner interfacial polarization theory. 0.00080 MHz/s is the optimal sweep rate necessary to predict the true DEP behavior of PS beads because it allows for full or partial (when the frequency is above 0.28 MHz) polarization.
- the slowest frequency sweep of 0.00080 MHz/s accurately predicts the static frequency behavior and is highly reproducible, with a very narrow 95% confidence interval range ( FIG. 4 c ).
- the fastest sweep rate of 0.056 MHz/s does not predict the static behavior of the human RBCs and is much less reproducible, as evidenced by the large 95% confidence interval in FIG. 4 c .
- FIGS. 5-13 show data obtained from applying a sweeping oscillating voltage signal of 1000Vpp/cm to human red blood cells in a solution with conductivity of 0.10 S/m.
- FIGS. 5-10 show data obtained from A+ RBCs while FIG. 11 shows data from A ⁇ RBCs, FIG. 12 shows data from O+ RBCs, and FIG. 13 shows data from B+ RBCs.
- FIG. 5 shows images of A+RBCs distributed in the vicinity of quadrapole electrodes during application of an oscillating voltage having frequencies ranging from 0.01 MHz to 1.0 MHz applied statically or at sweep rates of 0.00080 MHz/s, 0.0016 MHz/s, or 0.0028 MHz/s.
- the notation “(1)” denotes RBCs demonstrating nDEP behavior and “(2)” denotes RBCs demonstrating pDEP behavior.
- the RBC's DEP behavior at slower frequency sweep rates (0.00080 MHz/s and 0.0016 MHz/s) correlates well with the static frequency response (top row), indicating that these are below the maximum sweep rate for these conditions.
- FIG. 6 compares the nDEP (panel (a)) and the nDEP (panel (b)) intensity profiles of the particles for the statically applied oscillating voltage as well as at the various sweep rates tested across the range of frequencies.
- Panel (a) of each of FIGS. 7-13 shows images of particle distributions at oscillating voltages of 0.70 MHz and 0.80 MHz collected either with static application of the oscillating voltage (“0 MHz”) or while sweeping at the indicated sweep rates.
- Panels (b) and (c) of each of FIGS. 7-13 show nDEP (panels (b)) and pDEP (panels (c)) intensity profiles throughout the frequency range relative to intensity profiles obtained with statically-applied oscillating voltages at 0.7 MHz and 0.8 MHz.
- FIGS. 7-10 include results from two different experimental runs, R1 and R2.
- Comparison of the images of particle distributions at different sweep rates relative to distributions obtained with statically applied oscillating voltages provides an indication of whether or not the sweep rate is too fast, based on whether the swept images match those obtained with statically applied oscillating voltages at the same frequency.
- the data in this example shows the applicability of the disclosed methods and in particular the similarity in maximum sweep rates for various blood types.
- the data also shows the dependence of the maximum sweep rate on conductivity of the solution in which the particles (RBCs) are suspended.
- Example 3 Conditions for Example 3 arc the same as described above for Examples 1 and 2 except where otherwise stated.
- FIGS. 14-23 show the effect of changes in conductivity on the maximum sweep rate using red blood cells.
- biological material and in particular cells such as red blood cells
- varying combinations of NaCl and dextrose were combined to achieve the stated levels of conductivity of 0.10 S/m, 0.25 S/m, 0.50 S/m, and 1.0 S/m while maintaining the solution at approximately isotonic levels for human red blood cells; increasing the amount of NaCl increases the conductivity and proportionately less dextrose is used as NaCl is increased in order to maintain a relatively constant tonicity (see An et al. 2014, incorporated herein by reference).
- increasing conductance from 0.10 S/m to 1.0 S/m has the effect of increasing the maximum sweep rate that can be used from less than 0.0026 MHz/s at 0.1 S/m to less than 0.0031 MHz/s at 1.0 S/m.
- using a sweep rate below the maximum level generates particle distributions at various frequencies that are equivalent to distributions that are obtained with the application of an oscillating voltage at a static frequency.
- Increasing the conductance permits the use of a faster sweep rate, which in turn permits data to be collected at a faster rate.
- Example 4 Conditions for Example 4 are the same as described above for Examples 1-3 except where otherwise stated.
- Example 4 demonstrate that particle DEP behavior is independent of the ‘direction’ of sweeping, i.e. sweeping the oscillating voltage signal from a high frequency to a low frequency generates equivalent results as when the oscillating voltage is swept from a low frequency to a high frequency.
- FIG. 24 shows RBC static images (top row) compared to RBCs response using a sweep rate of ⁇ 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000 V pp /cm.
- FIG. 25 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.10 S/m at 1000 V pp /cm using the reverse sweep method as illustrated by the images in FIG. 24 .
- FIG. 26 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000 V pp /cm.
- FIG. 27 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.10 S/m at 1000 V pp /cm using the reverse sweep method as illustrated by the images in FIG. 26 .
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Dispersion Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Biophysics (AREA)
- Hematology (AREA)
- Urology & Nephrology (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
A system for identifying a particle. The system includes a microfluidic device; a microelectrode array including a plurality of electrodes, the microelectrode array disposed within the microfluidic device; a plurality of particles suspended in a solution and delivered to the micro-electrode array using the microfluidic device; a signal generator operatively coupled to the microelectrode array; a particle detector adjacent to the microelectrode array; and a controller in operative communication with the signal generator and the particle detector. The controller is configured to apply an oscillating voltage signal to the microelectrode array between a low frequency and a high frequency at a sweep rate, wherein the sweep rate is no more than a maximum sweep rate, and determine a distribution of the plurality of particles relative to the microelectrode array at a plurality of frequency levels between the low frequency and the high frequency.
Description
- This application claims priority to co-pending U.S. Provisional Patent Application No. 61/887,178 filed on Oct. 4, 2013, the entire content of which is incorporated herein by reference.
- The present invention was conceived while performing work under CBET 0644538, CBET 1041338, and IIP 1340126, each of which has been awarded by the National Science Foundation. The government has certain rights in the invention.
- The present invention relates to identification of particles based on dielectrophoretic responses.
- In one embodiment, a system for identifying a particle includes a microfluidic device; a microelectrode array including a plurality of electrodes, the microelectrode array disposed within the microfluidic device; a plurality of particles suspended in a solution and delivered to the microelectrode array using the microfluidic device; a signal generator operatively coupled to the microelectrode array; a particle detector adjacent to the microelectrode array; and a controller in operative communication with the signal generator and the particle detector. The controller is configured to apply an oscillating voltage signal to the microelectrode array between a low frequency and a high frequency at a sweep rate, wherein the sweep rate is no more than a maximum sweep rate, and determine a distribution of the plurality of particles relative to the microelectrode array at a plurality of frequency levels between the low frequency and the high frequency.
- In another embodiment, a method of identifying a particle. The method includes the steps of: placing a plurality of particles adjacent a microelectrode array, the microelectrode array including a plurality of electrodes; applying an oscillating voltage signal to the microelectrode array, the oscillating voltage signal varying between a low frequency and a high frequency at a sweep rate, wherein the sweep rate is no more than a maximum sweep rate; and determining a distribution of the plurality of particles relative to the microelectrode array at a plurality of frequency levels between the low frequency and the high frequency.
- Other features and aspects of the invention will become apparent by consideration of the following detailed description and accompanying drawings.
-
FIG. 1 shows (a) Dielectric relaxation mechanism for PS beads showing cases when i) particle polarization occurs at a static frequency, ii) τMW is shorter than the slow frequency sweep rate (τΔFS) allowing the bead interface time to polarize in response to the non-uniform AC field, and iii) τMW is longer than the τΔFS for fast frequency sweep rates and the bead interface does not have time to fully polarize. (b) Schematic of the quadrapole electrodes micro patterned onto a glass slide, and (c) microdevice with PDMS fluidic layer bonded above the quadrapole electrodes silver-epoxied to copper leads. -
FIG. 2 shows (a) nDEP behavior of 6 μm PS beads suspended in E-pure H2O 2.5×10−4 S/m and 250Vpp/cm 0.0063, 0.056 and 0.17 MHz/s sweep rates from 0.010 MHz to 1.0 MHz. (b) Raw intensity (arbitrary units) profile of PS beads in the center nDEP region (boxes shown at 0.20 MHz) at 0.0063 MHz/s sweep rate. Inset is a calibration of intensity per bead. (c) Clausius-Mossotti factor for the PS beads from 0.010 MHz to 2.0 MHz at three conductivities of 2.5×10−4, 1.0×10−3, and 1.0 S/m. PS bead assembly at slower frequency sweep rates track static frequency responses while 0.056 MHz/s illustrates transitional behavior and frequency sweeps above 0.17 MHz/s substantially lag the true static frequency DEP responses. -
FIG. 3 shows (a) 6 μm PS beads nDEP intensity profiles for 0.00080, 0.0063, and 0.056 MHz/s and static steady state (SS) measurements (black diamonds). Intensity analysis captures bead assembly to the quadrapole center with transient and SS regions. The slowest frequency sweep rate of 0.00080 MHz/s best predicts the static DEP responses. (b) Bead assembly intensity (arbitrary units) profiles for 0.0063 (n=8) and 0.17 MHz/s (n=7) with 95% confidence upper and lower limits shown as dashed lines. (c) Transient slope comparison for static frequencies (0 MHz/s) as well as frequency sweeps. (d) Comparison of static frequency and frequency sweep PS bead velocities from 0 to 50 s. 0.00080 MHz/s results are consistently similar to the static frequency results. -
FIG. 4 shows (a) nDEP behavior of RBCs suspended in 0.1 S/m dextrose buffer and 250V/cm at 0.00080 MHz/s, 0.0063 MHz/s and 0.056 MHz/s sweep rates from 0.010 MHz to 0.50 MHz. (b) RBCs nDEP intensity (arbitrary units) profiles for 0.00080, 0.0063, and 0.056 MHz/s and static measurements. (c) 0.00080 and 0.056 MHz/s RBC assembly intensity (arbitrary units) profiles n=8, with 95% confidence interval upper and lower limits shown as dashed lines. -
FIG. 5 shows an image comparison of the nDEP and pDEP behavior of A+ red blood cells suspended in 0.10 S/m dextrose solution and 1000Vpp/cm at 0.00080, 0.0016, and 0.0028 MHz/s sweep rates from 0.010 MHz to 1.0 MHz. (1) Denotes red blood cells nDEP behavior and (2) denotes red blood cell pDEP behavior. The red blood cells' DEP behavior at slower frequency sweep rates correlates well with the static frequency response (top row). -
FIG. 6 illustrates (a) nDEP and (b) pDEP intensity profiles for 0.00080, 0.0016, and 0.0028 MHz/s and static steady state measurements (solid circle). The intensity (arbitrary units) analysis captures the RBCs assembly toward the electrode center (nDEP) and near the electrodes (pDEP). 0.00080 MHz/s is the slowest and 0.0016 MHz/s is the fastest sweep rate to best predict RBCs' static DEP response. -
FIG. 7 illustrates (a) RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (middle and bottom row). R1 and R2 arc two different repeats completed for this measurement. (b) nDEP and (c) pDEP intensity (arbitrary units) profiles for 0.0024 MHz/s with RBC static steady state measurements (circles). The images and intensity profiles show the threshold frequency sweep rate at which agreement with static measurements is acceptable, but begins to falter. Sweep rates less than 0.0024 MHz/s agree well. Each test was completed with A+ blood in 0.10 S/m at 1000Vpp/cm. -
FIG. 8 illustrates (a) RBC static images (top row) compared to RBCs response using 0.0025 MHz/s (middle and bottom row). R1 and R2 are two different repeats completed for this measurement. (b) nDEP and (c) pDEP intensity profiles for 0.0025 MHz/s with RBC static steady state measurements (circles). The images and intensity profiles show poor agreement with static measurements. Each test was completed with A+ blood in 0.10 S/m at 1000Vpp/cm. -
FIG. 9 illustrates (a) RBC static images (top row) compared to RBCs response using 0.0026 MHz/s (middle and bottom row). R1 and R2 are two different repeats completed for this measurement. (b) nDEP and (c) pDEP intensity profiles for 0.0026 MHz/s with RBC static steady state measurements (circles). The images show poor agreement and intensity profiles show good agreement with static measurements. Each test was completed with A+ blood in 0.10 S/m at 1000Vpp/cm. -
FIG. 10 illustrates (a) RBC static images (top row) compared to the response of RBCs using a sweep rate of 0.0028 MHz/s (middle and bottom row). R1 and R2 are two different repeats completed for this measurement. Panels (b) nDEP and (c) pDEP show intensity profiles for 0.0028 MHz/s with RBC static steady state measurements (circles). The images show poor agreement and intensity profiles illustrate the lack of reproducibility of agreement with static measurements. Each test was completed with A+ blood in 0.10 S/m at 1000Vpp/cm. -
FIG. 11 illustrates (a) A− RBC static images (top row) compared to A− RBCs response using 0.0024 MHz/s (bottom row). (b) nDEP and (c) pDEP intensity profiles for 0.0024 MHz/s with RBC static steady state measurements (circles). The images and intensity profiles show good agreement with static measurements. Each test was completed in 0.10 S/m at 1000Vpp/cm. -
FIG. 12 illustrates (a) O+ RBC static images (top row) compared to O+ RBCs response using 0.0024 MHz/s (bottom row). (b) nDEP and (c) pDEP intensity profiles for 0.0024 MHz/s with RBC static steady state measurements (circles). The images and intensity profiles show good agreement at 0.70 MHz and fair agreement at 0.80 MHz. Each test was completed in 0.10 S/m at 1000Vpp/cm. -
FIG. 13 illustrates (a) B+ RBC static images (top row) compared to OB+ RBCs response using 0.0024 MHz/s (bottom row). (b) nDEP and (c) pDEP intensity profiles for 0.0024 MHz/s with RBC static steady state measurements (circles). The images and intensity profiles show good agreement over the tested frequency range 0.60-0.84MHz. Each test was completed in 0.10 S/m at 1000Vpp/cm. -
FIG. 14 shows a plot of sweep rate (MHz/s) as a function of conductivity of a solution. The plot points are a sweep rate at which the results for a given conductivity were no longer accurate. The plotted curve defines a threshold sweep rate. Thus, sweep rates above the threshold sweep rate for a given concentration are too fast, and any sweep rate below the curve may be used accurately and reliably. This aids in determining the fastest sweep rates that may be used for a solution with a given concentration in order to decrease the overall time needed for the procedure. -
FIG. 15 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000Vpp/cm. This illustrates the visual similarity between the static response and response at this sweep rate at the solution conductivity of 0.10 S/m. -
FIG. 16 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.50 S/m at 1000Vpp/cm. This illustrates the visual similarity between the static response and response at this sweep rate at the solution conductivity of 0.50 S/m. -
FIG. 17 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 1.00 S/m at 1000Vpp/cm. This illustrates the visual similarity between the static response and response at this sweep rate at the solution conductivity of 1.0 S/m. -
FIG. 18 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.25 S/m at 1000Vpp/cm. This plot generally shows the 0.0027 MHz/s sweep rates near the threshold sweep rate for this concentration match the static frequency data well, but the 0.0028 MHz/s sweep rate begins to falter. -
FIG. 19 shows RBC static images (top row) compared to RBCs response using 0.0027 MHz/s and 0.0028 Mhz/s (middle and bottom row). Each test was completed with A+ blood in 0.25 S/m at 1000Vpp/cm. This illustrates the visual similarity between the static response and the response at a sweep rate near the threshold sweep rate for the given conductivity of the solution. -
FIG. 20 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.5 S/m at 1000Vpp/cm. This plot generally shows the 0.0029 MHz/s sweep rates near the threshold sweep rate for this concentration match the static frequency data well, but the 0.0030 MHz/s sweep rate begins to falter. -
FIG. 21 shows RBC static images (top row) compared to RBCs response using 0.0029 MHz/s and 0.0030 Mhz/s (middle and bottom row). Each test was completed with A+ blood in 0.50 S/m at 1000 Vpp/cm. This illustrates the visual similarity between the static response and the response at a sweep rate near the threshold sweep rate for the given conductivity of the solution. -
FIG. 22 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 1.0 S/m at 1000Vpp/cm. This plot generally shows the 0.0030 MHz/s sweep rates near the threshold sweep rate for this concentration match the static frequency data well, but the 0.0030 MHz/s sweep rate begins to falter. -
FIG. 23 shows RBC static images (top row) compared to RBCs response using 0.0030 MHz/s and 0.0031 MHz/s (middle and bottom row). Each test was completed with A+ blood in 1.00 S/m at 1000 Vpp/cm. This illustrates the visual similarity between the static response and the response at a sweep rate near the threshold sweep rate for the given conductivity of the solution. -
FIG. 24 shows RBC static images (top row) compared to RBCs response using a sweep rate of −0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000 Vpp/cm. -
FIG. 25 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.10 S/m at 1000 Vpp/cm using the reverse sweep method as illustrated by the images inFIG. 24 . -
FIG. 26 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000 Vpp/cm. -
FIG. 27 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.10 S/m at 1000 Vpp/cm using the reverse sweep method as illustrated by the images inFIG. 26 . - Before any embodiments of the invention are explained in detail, it is to be understood that the invention is not limited in its application to the details of construction and the arrangement of components set forth in the following description or illustrated in the following drawings. The invention is capable of other embodiments and of being practiced or of being carried out in various ways. Also, it is to be understood that the phraseology and terminology used herein is for the purpose of description and should not be regarded as limiting.
- Alternating current (AC) dielectrophoresis (DEP) experiments for biological particles in microdevices have typically been applied at fixed frequencies. Reconstructing the DEP response curve from static frequency experiments is laborious, but is important for ascertaining differences in dielectric properties of biological particles. The disclosed systems and methods, on the other hand, employ the novel concept of sweeping the frequency as a function of time to rapidly determine the DEP response curve from fewer experiments. Homogeneous 6.08 μm polystyrene (PS) beads were initially used as a model system to determine whether sweeping the frequency would be a viable method for generating DEP responses and then to identify an optimal sweep rate. Subsequent experiments were performed using the sweep rate approach with −7 μm red blood cells (RBC) to verify that this approach would also work with biological samples. A Au/Ti quadrapole electrode microfluidic device was used to separately subject particles and cells to 10Vpp AC electric fields at frequencies ranging from 0.010-2.0 MHz over sweep rates from 0.00080 to 0.17 MHz/s. PS beads exhibited negative DEP assembly over the frequencies explored, likely due to Maxwell-Wagner interfacial polarizations. Results demonstrate that frequency sweep rates must be slower than particle polarization timescales; in some embodiments, sweep rates near 0.00080 MHz/s yielded DEP behaviors very consistent with static frequency DEP responses for both PS beads and RBCs, although higher sweep rates may also be employed.
- Accordingly, disclosed herein are systems and methods for identifying a particle using dielectrophoresis (DEP). Embodiments of the methods and systems disclosed herein may be used to distinguish between different types of particles based on differences in dielectric properties of the particles. In various embodiments, the particles that are analyzed may include polystyrene beads (e.g. for testing purposes) or cells such as blood cells; in particular, different subtypes of red blood cells (e.g. A+, A−, O+, O−, etc.) may be distinguished based on the surface charge differences of the red blood cell subtypes (e.g. due to differing antigens on the cell surfaces). In general, the particles may range in size from about 1-50 μm and should be detectable (e.g. through optical or electrical means) using the particle detector (e.g. an imaging system).
- A system according to embodiments of the invention may include a microfluidic device (which in some embodiments may include an enclosed microfluidic chamber) having a microelectrode array disposed in the fluid path of the microfluidic device. In various embodiments, particles are delivered to the vicinity of the microelectrode array using the microfluidic device prior to data collection. In certain embodiments, data collection is performed in a “batch-wise” manner, i.e. a group of particles is delivered to the microelectrode array and fluid movement is then stopped before data collection begins so that particle movements that are observed are due to dielectrophoresis.
- The microelectrode array in one embodiment is a quadrapole arrangement as shown, for example, in
FIGS. 1, 2, 4, and 5 . In this arrangement four electrodes are arranged in an “X” with a gap in the center (FIG. 5 ). Other possible arrangements of electrodes include interdigitated electrodes, V-shaped electrodes, circular electrodes, and T-shaped electrodes. In various embodiments, the electrodes are arranged so that oppositely charged electrodes are not parallel to one another, as this would create uniform fields whereas other, non-parallel geometries create non-uniform electric fields. In general, the electrodes are arranged so that they create a spatially non-uniform field. The electrodes are attached to the bottom of the microfluidic device and the particles that are delivered to the device are initially distributed in the vicinity of the electrodes in a random arrangement (e.g. see upper left panel ofFIG. 2a ) before any electrical signal is applied. The microelectrode array may be made by depositing metal strips onto a glass slide with a cover having a microfluidic channel being bonded on top of the glass slide (FIG. 1c ). As shown inFIG. 1 c, opposing pairs of electrodes may be electrically coupled using copper wires, as shown, with the leads (i.e. ground and “hot” AC signal) of a signal generator being connected to the copper wires. In various embodiments, pairs of electrodes may be electrically coupled as shown inFIG. 1c so that, upon stimulation, each electrode is 90° out of phase from the others to create a traveling wave signal. - Once particles have been delivered to the microelectrode array, a signal generator is used to deliver an oscillating voltage to the electrodes. In various embodiments, the voltage is applied at a peak-to-peak amplitude of 0.1Vpp, 1Vpp, 10Vpp, 100Vpp, or other suitable amplitude. In various embodiments, the oscillating voltage is applied at frequencies of at least about 0.001 MHz, at least about 0.005 MHz, at least about 0.01 MHz, at least about 0.05 MHz, at least about 0.1 MHz, at least about 0.5 MHz, or at least about 1.0 MHz. In other embodiments, the oscillating voltage is applied at frequencies of no more than about 10.0 MHz, no more than about 5.0 MHz, no more than about 2.0 MHz, no more than about 1.0 MHz, or no more than about 0.5 MHz.
- In particular embodiments, the frequency of the oscillating voltage is varied, for example from a low frequency to a high frequency, in order to collect data at a variety of different frequencies, a process referred to as “sweeping” the frequency. In various embodiments, comparable results are obtained when the frequency is swept from a high frequency to a low frequency. Sweeping the oscillating voltage using a continuously varying frequency permits a relatively large amount of data to be gathered in a short period of time. The rate at which the frequency sweep, i.e. the “sweep rate,” may vary from 0.00001 MHz/s to 0.1 MHz/s. In certain embodiments, the continuously varying frequency may be approximated by a series of discrete, step-wise changes in frequency with an increment ranging from about 10 nHz to about 10 Hz. As discussed herein, the optimum sweep rate may depend on conditions such as the conductivity of the solution in which the particles are suspended. The present inventors have found that when the frequency is swept above a certain maximum sweep rate the frequency changes too quickly, such that the particles do not have sufficient time to respond to the voltage signal at a given frequency before the signal changes to the next frequency. If the oscillating voltage signal is varied too quickly, i.e. above the maximum sweep rate, the observed particle movements and distributions will be inaccurate and could lead to an inconclusive or erroneous particle identification. Thus, in certain embodiments the maximum sweep rate is no more than about 0.003 MHz/s, no more than about 0.0029 MHz/s, no more than about 0.0028 MHz/s, no more than about 0.0027 MHz/s, no more than about 0.0026 MHz/s, no more than about 0.0025 MHz/s, no more than about 0.0020 MHz/s, no more than about 0.0015 MHz/s, no more than about 0.0010 MHz/s, no more than about 0.0008 MHz/s, or no more than about 0.0005 MHz/s. In various embodiments, a minimum sweep rate of at least about 0.00005 MHz/s, at least about 0.0001 MHz/s, at least about 0.00015 MHz/s, at least about 0.0002 MHz/s, at least about 0.0004 MHz/s, at least about 0.0005 MHz/s, at least about 0.00075 MHz/s, or at least about 0.0010 MHz/s may be used.
- While the oscillating voltage is being applied to the microelectrode array at varying frequencies, data may be collected to determine the spatial distribution of the particles within the microfluidic device, particularly the particles in the vicinity of the electrodes, as a function of time. Particle detection may be carried out with systems which are capable of identifying the spatial distributions of the particles with sufficient temporal (e.g. operating at 0.1-10 Hz) and spatial (e.g. capable of resolving 0.1 μm×0.1 μm areas) resolution.
- In some embodiments, images are collected at regular intervals (e.g. at video rates of 30 frames/sec or at other, slower rates such as 1 image or frame/sec) while sweeping the oscillating voltage. The images may be processed (for example several sequential video-rate frames may be averaged together) and the images or subregions thereof may be analyzed to characterize particle distribution and behavior at one or several frequencies. The analyses may include one or more of determining the particles' intensity profiles, transient responses, and velocities. Analyses may be conducted at one or more discrete locations within the images including along one or more lines, e.g. lines running between electrode tips. Analysis patterns for a sample which includes an unknown particle may be compared to patterns generated under equivalent conditions using known particles to determine the identity of the unknown particle. In certain embodiments in which only a few (e.g. 2-3) types of particles need to be distinguished, it may be sufficient to analyze only two subregions of the particle distribution in order to reliable distinguish the particle types from one another. In other embodiments in which a larger number of particle types are possible it may be necessary to analyze more subregions in order to reliably distinguish among particle types. For example, in order to distinguish among the eight red blood cell subtypes (A+, A−, B+, B−, O+, O−, AB+, and AB−) it may be necessary to analyze at least 4 different subregions of the particle distribution, and greater accuracy would be achieved by including more of the particle distribution in the analysis. In general, regions and patterns on the substrate are selected for analysis based on the areas that are expected to have the greatest change in electric field patterns and hence the greatest change in particle distribution at different frequencies, so as to provide the most information for distinguishing between particle types. Additional methods for analyzing DEP behavior of particles are disclosed in Salmanzadeh et al. (2012), Rozitsky et al. (2013), and An et al. (2014), each of which is incorporated herein by reference.
- A system for carrying out embodiments of the invention may include a controller for carrying out or more of the procedures disclosed herein. The controller may be in operative communication with one or more of the signal generator (which in turn is in communication with the microelectrode array) and the particle detector (e.g. imaging system). The controller may be a part of or in communication with a computer system. The computer system may be part of an existing computer system (e.g. on a smartphone, desktop computer, on-board computer, etc.) or may be implemented as a separate, standalone unit that is in local or remote communication with other components. The computer system(s) may be in wired or wireless communication with other systems through a combination of local and global networks including the Internet. Each computer system may include one or more input device, output device, storage medium, and processor (e.g. a microprocessor). Input devices may include a microphone, a keyboard, a computer mouse, a touch pad, a touch screen, a digital tablet, a track ball, and the like. Output devices include a cathode-ray tube (CRT) computer monitor, an LCD or LED computer monitor, touch screen, speaker, and the like.
- The computer system may be organized into various modules including an acquisition module and an output module along with the controller, where the controller is in communication with the acquisition module and the output module. The various modules for acquiring and processing data and for returning a result may be implemented by a single computer system or the modules may be implemented by several computer systems which are in either local or remote communication with one another.
- Storage media include various types of local or remote memory devices such as a hard disk, RAM, flash memory, and other magnetic, optical, physical, or electronic memory devices. The processor may be any known computer processor for performing calculations and directing other functions for performing input, output, calculation, and display of data in accordance with the disclosed methods. In various embodiments, implementation of the disclosed invention includes generating sets of instructions and data that are stored on one or more of the storage media and operated on by a controller, where the controller may be configured to implement various embodiments of the disclosed invention.
- Some embodiments the system may be in the form of a portable or handheld device which may be self-contained, including the microfluidic device, microelectrode array, the controller, and additional components such as a power supply and input/output capabilities.
- Without being limited as to theory, DEP enables phenotypically similar biological cells to be discriminated based on dielectric properties including the conductivity and permittivity of the membrane, cytoplasm, and other structurally relevant organelles. Cell components and structure contribute to a cell's signature dielectric dispersion. A particle's complex permittivity is frequency dependent and characterized by dielectric dispersion regions (y, β, and α, where wα<wβ<wy) specific to an applied frequency. Certain embodiments of this work to illustrate sweep rates uses frequencies in the range of 0.010 to 2.0 MHz in the β-dispersion region because the Clausius-Mossotti factor, which governs sign and polarization strength, for polystyrene beads is nearly constant over this range. Maxwell-Wagner theory describes the polarization mechanism of particles in the β-dispersion region as interfacial polarization where moving charges build around the interface of a charged or charge-neutral particle and exchange ions with the suspending medium. Interfacial particle polarization creates an induced dipole moment such that the particle experiences disproportionate forces in each half cycle of the alternating current (AC) field resulting in net particle movement.
- Polarized particles can exhibit either positive dielectrophoresis (pDEP) or negative dielectrophoresis (nDEP) as a consequence of the frequency-dependent polarizability of the particle in the surrounding medium. Particles that exhibit pDEP move to high electric field regions and particles that exhibit nDEP move to low electric field regions. This motion up and down electric field gradients is described by the Clausius-Mossotti factor for spherical particles.
-
- where {tilde over (ϵ)}i is the complex permittivity of the particle (i=p) and of the medium (i=m), which are both functions of conductivity (σ), permittivity (ϵ), and angular frequency (ω).
- Polarization is not an instantaneous event; charge transport into the interface takes a few microseconds in response to the electric field. Maxwell-Wagner dielectric relaxation is a physical phenomenon related to the transport delay of cation and anion alignment in and around the interface of the dielectric particle. At lower frequencies (<˜10 MHz), particle polarization is driven by this conductive polarization. At higher AC frequencies, charges do not have enough time to move into and around the interface double layer, so particles experience polarization lag time as a result of the rapidly modulating field and do not reach maximum polarization.
- Maxwell-Wagner dielectric relaxation is characterized by a time constant, τMW, which is unique to each particle or cell due to the time constant's dependence on the cell dielectric properties. The time required for a particle to reach maximum polarization is given by Eq. (3) (see Morgan et al. (2003), p. 27; see also Grosse et al. (2010) and Mittal et al. (2008), each of which is incorporated herein by reference):
-
- Typical relaxation times for particle polarization vary from pico- to microseconds (see Morgan et al. (2003), p. 27; see also Grosse et al. (2010) and Mittal et al. (2008)), and the calculated τMW for polystyrene (PS) beads in our Epure H2O medium at 2.5×10−4 S/m is 3.5 las. Thus, a single AC cycle is on the order of 0.01 to 2μs; the time delay in ion transport within a static frequency field of 0.010 to 2.0 MHz is such that 2 to 350 AC cycles must be completed before the particle experiences full polarization.
- The Maxwell-Wagner dielectric timescale for charge transport into and around the interface becomes important when the frequency is swept, i.e. changes as a function of time.
FIG. 1a highlights the Maxwell-Wagner particle polarization at the interface under static frequency as well as slow and fast frequency sweep rates. At a static frequency in the β-dispersion region, the particle experiences a constant frequency field such that the relaxation time is not a factor and the particle fully polarizes. A particle in a field with a slowly changing frequency sweep has a relaxation time, τΔFS, that is less than τMW and thus the particle interface fully polarizes. Tn contrast, a particle in a fast frequency sweep has a relaxation time, τΔFS, that is larger than τMW and the particle interface does not have time to fully polarize in the field. PS beads are lossy dielectric particles treated as homogeneous spheres and are thus an idealized particle to examine new techniques, devices, or approaches to dielectrophoretic characterizations. Our system is easily able to discern pDEP and nDEP transitional behavior and adaptable to new frequency sweep techniques. The homogeneous spherical DEP polarization model for PS beads (ϵ=2.5 and σ=9.4×10−5 S/m) suspended in Epure H2O displays only nDEP behavior over 0.010 to 2.0 MHz. - Thus, microfluidic and dielectrophoretic (DEP) technologies enable a wide variety of particle polarizations with nonuniform electric fields on microchips to achieve particle manipulation, concentration, separations, and property-based identification. Particles can include bioparticles (e.g. DNA, viruses, or proteins) as well as cells (e.g. blood cells, cancer cells, stem cells, and yeast). The advantages to coupling DEP with microfluidics are small sample size (on the order of microliters), rapid analysis (approximately minutes to achieve results), minimal sample preparation, and minimal waste production. Traditionally, DEP experiments are completed at static, fixed frequencies such that maximum particle polarization can be achieved and measured. Multiple experiments are conducted, each at discrete frequencies over the range of interest to stitch together DEP response spectra; this is a labor-intensive approach. Further disadvantages are that extended field exposure times at fixed frequencies can change particle properties or cell viability. As disclosed herein, it is demonstrated that frequency can be swept with time in the β-dispersion region thus enabling interrogation of cells at multiple frequencies within a short time period. The benefits of using a frequency sweep technique are that nearly continuous DEP response curves, when coupled with automated response analysis, can be compiled in near real time and the number of experiments needed to obtain particle DEP spectra are greatly reduced.
- Traditional DEP measurements are completed at single static frequencies in order to compile frequency by frequency, the DEP spectrum for a particle or cell system. This method is laborious and, as disclosed herein, requires time for particles to fully polarize for accurate observed DEP responses. The present disclosure describes the use of frequency sweeps as a means to more efficiently interrogate multiple frequencies in a single experimental run and systematically compared the responses to the nDEP response at fixed frequencies between 0.010 and 2.0 MHz. It was observed that frequency sweep rates influence the DEP response of PS beads and RBCs and further, the permissible frequency sweep rate is particle or cell dependent. The underlying mechanism appears to be the same. At slower sweep rates, particles have more time to polarize in the electric field and therefore a more accurate and reproducible DEP spectrum can be obtained. At faster frequency sweep rates, the particles are unable to achieve maximum interfacial polarization because of the dielectric relaxation time scale so the observed DEP response does not match the true DEP behavior of the particle.
- For polystyrene beads at frequency sweep rates below 0.0063 MHz/s, responses correlate closely with dielectric responses of particles subjected to a static frequency potential. In the PS bead system, 0.056 MHz/s is the transitional sweep rate where the particle dielectric relaxation is approximately the same order of magnitude as the shifts in frequency within the sweep. Dielectric responses continue to track the static frequency responses, although reproducibility is diminished. However as this sweep rate is increased further, conductivity dominated interfacial polarizations cannot be established and the PS bead frequency sweep data does not coincide with static frequency measurements.
- For full utility in DEP experiments, this frequency sweep rate methodology must be translatable to cell systems. Results illustrated that only 0.00080 MHz/s accurately predicted the static frequency DEP responses of human RBCs. Red blood cells are substantially more morphologically and dielectrically complex than polystyrene beads. Calculation of the dielectric relaxation time, taking into account only the membrane permittivity and conductivity of 4.4 and 10−7 S/m, respectively (Gascoyne et al. (2004), incorporated herein by reference) yields a dielectric relaxation time ˜4.6 μs roughly corresponding to 0.21 MHz. This relaxation time is larger than the PS bead relaxation time of 3.5 μs, so the optimal frequency sweep rate for red blood cells would be slower than that for PS beads. This result suggests that for each new cell system of interest it is imperative to determine the optimal frequency sweep rate for accurately and reproducibly interrogating the behavior of that cell. This work outlines a systematic technique to make comparisons between frequency sweep rate and static frequency shown. For all cell systems, sweep rates that are too fast will not allow the cell adequate time to polarize and will result in inaccurate and less reproducible DEP responses. An optimal frequency sweep rate can be estimated by calculating the Maxwell-Wagner dielectric relaxation time for the particle/cell of interest, provided the cell's permittivity and conductivity is known. The frequency sweep rate chosen for the DEP study should then remain at frequencies below the inverse dielectric relaxation time (1/τMW) for 5-45 s (longer times spent below the threshold give better DEP predictions).
- Since the cell's permittivity and conductivity are determined from the frequency dependent DEP spectrum, this presents a cyclical situation. However, this work has demonstrated that frequency sweep rates slower than 0.00080 MHz/s can yield accurate DEP response of PS beads as well as RBCs. This sweep rate may therefore be translatable to other cell systems. In addition, at higher frequencies where the polarization mechanism is more heavily influenced by charge permittivity effects through the membrane and cell cytosol, it is possible that slow frequency sweep rates can still accurately capture DEP response spectra. Lastly, this frequency sweep rate technique will enable researchers to obtain accurate and continuous DEP response spectra in shorter experiment times.
- The following non-limiting Examples are intended to be purely illustrative, and show specific experiments that were carried out in accordance with embodiments of the invention.
- In this Example, dielectrophoretic responses of PS beads (model system) were quantified at both static frequencies and frequency sweeps at rates ranging from 0.00080 to 0.17 MHz/s over the β-dispersion frequency range of 0.010-2.0 MHz. PS bead motion in the electric field was imaged with video microscopy and analyzed using three techniques: intensity profiles, transient response, and particle velocities. Data shows that frequency sweep rates impact particle polarization due to dielectric relaxation limitations. This frequency sweep technique was further extended in this Example to negatively charged biconcave red blood cells (RBCs), which are an important cellular system for medical disease diagnostics.
- The microdevice shown in
FIG. 1c was fabricated according to previously published microfabrication techniques (Grom et al. (2006), incorporated herein by reference), with the 100 μm wide electrodes spaced 200 μm apart aligned at 90° along the bottom of a 70 pm deep by 1000 μm wide microfluidic chamber as shown inFIG. 1 b. Polystyrene beads (Cat No. PP-60-10, Spherotech, Lake Forest, Ill., USA), 6.08 μm in diameter were centrifuged at 1300 min−1 for 5 mins to separate the beads from the liquid. The PS beads were resuspended in Epure H2O (18 MΩ or 2.5×10−4 S/m) at a 1:10 (bead to water) volumetric dilution ratio and vortexed. Microdevice was pre-rinsed with Epure H2O and Alconox precision cleaner (Cat No. 1104, Alconox Inc, White Plains, N.Y., USA) to prevent bead adhesion. PS beadEpure H2O suspension was pumped to the microchamber using a syringe. Time was allowed for inlet and outlet pressures to equalize and flow to stop. The function generator (Agilent 33250A, Agilent, Santa Clara, Calif., USA) was connected via copper leads to produce a 10Vpp AC sine wave with frequencies ranging from 0.010-2.0 MHz at specific frequency sweep rates 0.00080, 0.0011, 0.0030, 0.0063, 0.013, 0.021, 0.028, 0.042, 0.056, 0.083, and 0.17 MHz/s. Frequency sweeps linearly increased the applied frequency as a function of time. Greater than five (n>5) static frequency experiments were completed at each frequency 0.010, 0.020, 0.030, 0.040, 0.050, 0.20, 0.40, 0.60, 0.80, 1.0, 1.2, 1.4, 1.6, 1.8, and 2.0 MHz by applying 10Vpp for 90 s. These DEP static frequency responses were compared to each frequency sweep rate DEP responses. For the static and frequency sweep experiments, the PS bead concentration was between 238-263 beads in the t=0 field of view. Video recordings of experiments were taken at 30 fps at 640×480 pixels/image using LabSmith SVM Synchronized Video Microscope with a 10× objective (LabSmith, Livermore, Calif., USA). - Video recordings of PS beads DEP behaviors were analyzed with ImageJ (NIH, Bethesda, Md.) using intensity, transient slope, and velocity measurements. Since PS beads only exhibit nDEP over the frequency range of interest, intensity data acquisition from images was completed by drawing a rectangular box at the device center, ICTR, and background, IBK measured in a location with no PS beads present (See
FIG. 2a ). ImageJ Z Project function was used to average the pixel intensities in the specified boxed region. The initial background, IBK(t=0) and center intensity, ICTR(t=0) were subtracted from the center and background intensity at each time, ICTR(t) and IBK(t), and then a normalized intensity was calculated by dividing by the maximum intensity experienced by the PS beads, (Eq. (4)): -
- This normalized intensity tracked the real-time magnitude of the PS bead DEP response, which had two distinct regions: transient where beads moved with nDEP toward the center, and steady-state (SS) where beads achieved tight packing at the device center. These two responses were analyzed separately via transient slope and particle velocity.
- The transient response of the PS beads was extracted from the steady-state response via signal processing in which the delay and rise time were quantified. The PS bead delay time, td, was characterized as the time required for the intensity response to reach 50% of the final intensity response for the first time. The rise time, tT, was determined as the time needed for the intensity response to reach 100% of the final intensity response for the first time (Ogata et al. (1978), pp. 517-518, incorporated herein by reference). This allowed the transient response to be segmented and a linear trend line was fit between td and tT where td<tT. A comparison of the transient slope for frequency sweep rates and static frequency measurements is given in
FIG. 3c . PS bead velocities were determined from the original video by tracking the x-, y-pixel position of individual PS beads from 0-50 s. PS bead located within 5 μm of electrode tips were selected to control for similar electric field gradients. This procedure was repeated for at least 10 beads in each specific frequency sweep rate and static experimental video. - For experiments involving human RBCs, blood of the appropriate type (e.g. O+, A+, etc.) was obtained from a single donor and centrifuged at 1400 rpm for 5 mins to separate the packed RBCs from the plasma and leukocytes. The packed RBCs were removed, then resuspended at 1:75 v:v in 0.10 S/m isotonic dextrose buffer doped with 0.75% BSA (Cat No. A7906, Sigma Aldrich, St. Louis, Mo., USA) to prevent cell/device adhesion. This RBC suspension was syringe-pumped to the microchamber, with time being allowed for flow to stop after pumping before the 10Vpp signal was applied over 0.010-0.50 MHz (range reduced to avoid pDEP behavior) at frequency sweep rates of 0.00080, 0.0063 and 0.056 MHz/s (n=7). RBC static frequency experiments were completed at 0.010, 0.10, 0.25 and 0.50 MHz at 10Vpp for 90 s (n=7). Video microscopy at 25× and 1 fps was obtained with a Zeiss Axiovert Inverted Light Microscope (Zeiss, Germany). The video images were analyzed as described herein for the PS beads.
- Frequency sweep rates ranging from 0.00080 to 0.17 MHz/s were explored to see if the nDEP response of PS beads would vary and/or correspond to static frequency measurements. The frequency range was chosen for the relatively consistent Clausius-Mossotti factor, Re(fCM) for a homogeneous lossy polystyrene sphere of 0.26 to 0.48 (see
FIG. 2c ) over the frequency range of 0.010 to 2.0 MHz. Static frequency experiments were completed at fixed values in this same frequency range.FIG. 2a shows still images from both static frequency experiments and the frequency sweeps at 0.20, 0.60, and 1.0 MHz. For static frequencies, the response 45 seconds after field application is shown while for frequency sweeps of 0.0063, 0.056, and 0.17 MHz/s, the image is shown at the time stamp when the specified frequency is reached. The electrodes are visible as black shadows in the images and the PS beads assemble due to nDEP forces at the central electric field gradient minima Data was examined to determine the sweep rate that most closely approximated the static frequency response. Frequency sweeps 0.00080 and 0.0063 MHz/s (shown) tracked static frequency, or true, DEP responses while the slightly faster sweep of 0.056 MHz/s begins to lag the true DEP responses and at 0.17 MHz/s and faster, particles were unable to achieve sufficient polarization to respond sufficiently in the electric field. - nDEP responses were quantified via intensity analysis as described herein for all sweeps and all static frequency experiments.
FIG. 2b illustrates the frequency- (and time-) dependent intensity for the 0.0063 MHz/s sweep rate images shown inFIG. 2a . This quantification of the PS bead nDEP response was correlated to total bead packing via the calibration shown in the inset. The 188-bead count at the center deviates slightly from the initial, field off, bead count of 245 because PS beads also move down the electric field gradient to regions outside of the image field of view. - Normalized intensities, Eq. (4), were compiled in
FIG. 3a for SS (i.e. 45 seconds) static frequency nDEP responses and 0.00080, 0.0063, 0.056 MHz/s frequency sweep rate nDEP responses. The time for sweep responses to achieve the true nDEP static response decreases as the sweep rate decreases. Frequency sweep rates 0.00080 and 0.0063 MHz/s are within the 95% confidence intervals (n=7) of the static steady-state (SS) responses.FIG. 3a inset shows that the slowest 0.00080 MHz/s sweep rate more quickly aligns closely with the static frequency responses.FIG. 3b compares average 0.0063 MHz/s (n=8) to 0.17 MHz/s (n=7) with the dashed lines signifying the upper and lower limits of the 95% confidence intervals for IDEP. The confidence intervals around the transient 0.0063 MHz/s sweeps are smaller than for 0.17 MHz/s over much of the frequency range indicating greater reproducibility at slower sweep rates. Faster sweep rates either do not reach SS or have a lag before reaching SS (compare toFIG. 2a ) suggesting the bead interface does not fully polarized and thus displays attenuated nDEP motion. - The transient behavior was quantified for all static frequencies and frequency sweeps via a transient slope analysis as compiled in
FIG. 3c . Four static frequency measurements 0.010, 0.60, 1.0 and 2.0 MHz are shown compared to 0.00080, 0.0063, 0.028, 0.056, and 0.17 MHz/s frequency sweep rates. Static frequency transient slopes range between 0.023-0.095 and are within the 95% (p<0.05) confidence intervals of 0.00080, 0.0063, and 0.028 MHz/s frequency sweep transient slopes. These slower sweep rates and 0.056 MHz/s differ at p<0.001 from the fastest sweep rate of 0.17 MHz/s, which is also significantly different at p<0.001 from the static measurements (except for 1.0×104 Hz with p<0.01). - Individual bead velocities were compiled for static as well as frequency sweeps in
FIG. 3d . PS bead velocity corroborates the intensity profile and the slope analysis that 0.00080 MHz/s frequency sweep rate closely tracks the bead velocity at static frequencies. 0.056 MHz/s gives good estimations of static frequency bead velocity at times greater than 20 s. Based on intensity, transient slope, and velocity analysis, the slow frequency sweep rate of 0.00080 MHz/s is most consistent with static frequency DEP responses. - There is an observable inverse relationship between the frequency sweep rate and particle polarization, where slower sweep rates result in comparable particle polarization characteristics to static frequency responses. Dielectric relaxation is the driving force of this relationship; the calculated dielectric relaxation time Eq. (3) for PS beads in E-pure H2O at 2.5×10−4 S/m is 3.5 ps, which corresponds to 0.28 MHz. There are two timescales that influence this behavior: the frequency itself and the change in frequency per time. The Maxwell-Wagner, conductivity-driven interfacial polarization mechanism occurs below ˜0.28 MHz; above this frequency threshold the interfacial polarization of the PS beads gradually decreases and the particle permittivity increasingly influences the DEP force. The experimental frequencies tested were within the range dominated by Maxwell-Wagner polarization such that maximum particle interfacial polarization was possible.
- The second timescale of interest is the frequency change per time or frequency sweep rate, which determines how many consecutive cycles a particle experiences a specific frequency. At slower sweep rates, the PS beads experience a specific frequency for a large number of cycles and thus the beads have time to polarize because the timescale of the frequency change is slower than the dielectric relaxation time. A particle must experience a single frequency during the sweep for a minimum of 3.5 μs for maximum interfacial polarization to be achieved. Upon polarization, the particle, which its current DEP force has to overcome inertia and Stokes drag to achieve observable particle motion down the electric field gradient. At static frequencies, it takes roughly 5 s for maximum velocity to be attained (see
FIG. 3d , AC field applied at t=5 s) and as much as 45 s for final SS at the field gradient minima to be reached. As the sweep rate increases, the dielectric relaxation time and the rate of change of the frequency approach the same order of magnitude. Results suggest that 0.056 MHz/s is a transitional sweep rate because the DEP behavior roughly corresponds to the static behavior of the PS beads. With further increases in frequency sweep rates, the timescale for frequency change surpasses the dielectric relaxation timescale such that particles are unable to fully polarize resulting in an attenuated DEP response as shown with data inFIGS. 2, 3 a, and 3 b.FIG. 3b also demonstrates that the transient behavior of the PS beads is more reproducible at slower frequency sweep rates, which can be attributed to the interfacial polarization timescale of the beads. Implications of the intensity, slope, and velocity analysis compared with static frequencies are that slow frequency sweep rates accurately predict the DEP response of PS beads because the changes in frequency are slower than the characteristic Maxwell-Wagner dielectric relaxation. - Thus, a frequency sweep approach can be utilized to attain accurate DEP behavior of PS beads, provided the sweep rate is slower than conductivity mediated interfacial polarization timescale. This result is reliable over frequency ranges where particle polarization is dominated by the conduction of free charges from the media. The charges are moving around the PS beads through the particle-liquid interface inducing a dipole, which causes PS bead movement down the electric field gradient to the electrode center. At different sweep rates the rate of movement of the charges varies which varies the rate of the dipole being induced, observed as dielectric relaxation. Each sweep rate has a unique dielectric relaxation time and our results are consistent with Maxwell-Wagner interfacial polarization theory. 0.00080 MHz/s is the optimal sweep rate necessary to predict the true DEP behavior of PS beads because it allows for full or partial (when the frequency is above 0.28 MHz) polarization.
- Given that the sweep methodology yielded accurate DEP responses for the ideal system of PS beads, the same methodology and frequency sweep rates were explored with human RBCs. The three most successful PS bead frequency sweep rates were reproduced with human red blood cells: 0.00080 MHz/s, 0.0063 MHz/s and 0.056 MHz/s. Static frequency experiments were also performed at 0.010 MHz, 0.10 MHz, 0.25 MHz and 0.50 MHz. Seen in
FIG. 4a are 25× microscope images taken of the t=45 s final static frequency frames aligned above the sweep time points that correspond to those four static frequencies. Qualitatively, the only sweep rate that accurately matches the static frequency behavior of the human RBCs is 0.00080 MHz/s. This behavior was further verified by the same intensity analysis as for PS beads. InFIG. 4b , the scaled intensity is plotted for 0.00080, 0.0063 and 0.056 MHz/s experiments (n=8) as compared to the static frequency intensities. After the initial 10 s transition for the red blood cells to polarize and overcome drag, the slowest frequency sweep of 0.00080 MHz/s accurately predicts the static frequency behavior and is highly reproducible, with a very narrow 95% confidence interval range (FIG. 4c ). The fastest sweep rate of 0.056 MHz/s does not predict the static behavior of the human RBCs and is much less reproducible, as evidenced by the large 95% confidence interval inFIG. 4c . From these experiments, we conclude that the optimal frequency sweep for determining the accurate DEP behavior of RBCs is 0.00080 MHz/s. Due to the complex dielectric properties of cells, it is necessary to carefully compare frequency sweep rates with static frequency behaviors to ascertain optimal frequency sweep rates that accurately interrogate the cell of interest. - Conditions for Example 2 are the same as described above for Example 1 except where otherwise stated.
FIGS. 5-13 show data obtained from applying a sweeping oscillating voltage signal of 1000Vpp/cm to human red blood cells in a solution with conductivity of 0.10 S/m.FIGS. 5-10 show data obtained from A+ RBCs whileFIG. 11 shows data from A− RBCs,FIG. 12 shows data from O+ RBCs, andFIG. 13 shows data from B+ RBCs. -
FIG. 5 shows images of A+RBCs distributed in the vicinity of quadrapole electrodes during application of an oscillating voltage having frequencies ranging from 0.01 MHz to 1.0 MHz applied statically or at sweep rates of 0.00080 MHz/s, 0.0016 MHz/s, or 0.0028 MHz/s. The notation “(1)” denotes RBCs demonstrating nDEP behavior and “(2)” denotes RBCs demonstrating pDEP behavior. The RBC's DEP behavior at slower frequency sweep rates (0.00080 MHz/s and 0.0016 MHz/s) correlates well with the static frequency response (top row), indicating that these are below the maximum sweep rate for these conditions.FIG. 6 compares the nDEP (panel (a)) and the nDEP (panel (b)) intensity profiles of the particles for the statically applied oscillating voltage as well as at the various sweep rates tested across the range of frequencies. - Panel (a) of each of
FIGS. 7-13 shows images of particle distributions at oscillating voltages of 0.70 MHz and 0.80 MHz collected either with static application of the oscillating voltage (“0 MHz”) or while sweeping at the indicated sweep rates. Panels (b) and (c) of each ofFIGS. 7-13 show nDEP (panels (b)) and pDEP (panels (c)) intensity profiles throughout the frequency range relative to intensity profiles obtained with statically-applied oscillating voltages at 0.7 MHz and 0.8 MHz.FIGS. 7-10 include results from two different experimental runs, R1 and R2. - Comparison of the images of particle distributions at different sweep rates relative to distributions obtained with statically applied oscillating voltages provides an indication of whether or not the sweep rate is too fast, based on whether the swept images match those obtained with statically applied oscillating voltages at the same frequency.
- The data in this example shows the applicability of the disclosed methods and in particular the similarity in maximum sweep rates for various blood types. The data also shows the dependence of the maximum sweep rate on conductivity of the solution in which the particles (RBCs) are suspended.
- Conditions for Example 3 arc the same as described above for Examples 1 and 2 except where otherwise stated.
- One factor which may affect the maximum sweep rate is the conductivity of the solution in which the particles are suspended. Accordingly, experiments were carried out to determine the extent to which conductivity impacts the maximum sweep rate.
-
FIGS. 14-23 show the effect of changes in conductivity on the maximum sweep rate using red blood cells. When using biological material, and in particular cells such as red blood cells, it is important when varying the conductivity of the solution to maintain the overall tonicity of the solution within a limited range. For the experiments inFIGS. 14-23 , varying combinations of NaCl and dextrose were combined to achieve the stated levels of conductivity of 0.10 S/m, 0.25 S/m, 0.50 S/m, and 1.0 S/m while maintaining the solution at approximately isotonic levels for human red blood cells; increasing the amount of NaCl increases the conductivity and proportionately less dextrose is used as NaCl is increased in order to maintain a relatively constant tonicity (see An et al. 2014, incorporated herein by reference). - As shown in
FIG. 14 , increasing conductance from 0.10 S/m to 1.0 S/m has the effect of increasing the maximum sweep rate that can be used from less than 0.0026 MHz/s at 0.1 S/m to less than 0.0031 MHz/s at 1.0 S/m. As seen inFIGS. 15-23 , using a sweep rate below the maximum level generates particle distributions at various frequencies that are equivalent to distributions that are obtained with the application of an oscillating voltage at a static frequency. Increasing the conductance permits the use of a faster sweep rate, which in turn permits data to be collected at a faster rate. - Conditions for Example 4 are the same as described above for Examples 1-3 except where otherwise stated.
- The experiments of Example 4 demonstrate that particle DEP behavior is independent of the ‘direction’ of sweeping, i.e. sweeping the oscillating voltage signal from a high frequency to a low frequency generates equivalent results as when the oscillating voltage is swept from a low frequency to a high frequency.
-
FIG. 24 shows RBC static images (top row) compared to RBCs response using a sweep rate of −0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000 Vpp/cm.FIG. 25 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.10 S/m at 1000 Vpp/cm using the reverse sweep method as illustrated by the images inFIG. 24 . -
FIG. 26 shows RBC static images (top row) compared to RBCs response using 0.0024 MHz/s (bottom row). Each test was completed with A+ blood in 0.10 S/m at 1000 Vpp/cm.FIG. 27 shows pDEP and nDEP plots of scaled intensity versus frequency for A+ blood in 0.10 S/m at 1000 Vpp/cm using the reverse sweep method as illustrated by the images inFIG. 26 . - Each of the following references is incorporated herein by reference in its entirety:
- A. Salmanzadeh, L. Romero, H. Shafiee, R. C. Gallo-Villanueva, M. A. Stremler, S. D. Cramer and R. V. Davalos, Lab on a Chip 12, 182 (2012).
- L. Rozitsky, A. Fine, D. Dado, S. Nussbaum-Ben-Shaul, S. Levenberg, G. Yossifon, Biomed Microdevices 15, 859 (2013).
- R. An, D. O. Wipf, A. R. Minerick, Biomicrofluidics 8,
issue 2, article 021803 (2014). - H. Xie, R. Tewari, H. Fukushima, J. Narendra, C. L. Heldt, J. King, A. Minerick, J. Vis. Exp. 88, e51696 (2014).
- F. Grom, J. Kentsch, T. Muller, T. Schnelle, and M. Stelzle, Electrophoresis 27, 1386 (2006).
- C. Grosse and A. V. Delgado, Current Opinion in Colloid & Interface Sciences 15, 145 (2010).
- H. Morgan and N.G. Green, R. Pethig, AC Electrokinetics: colloids and nanoparticles (Research Studies Press Limited, Philadelphia, 2003).
- M. Mittal, P. P. Lele, E. W. Kaler, and E. M. Furst, Journal of Chemical Physics 129, 065413 (2008).
- K. Ogata, System Dynamics (Prentice-Hall Inc, Englewood Cliffs, 1978).
- P. Gascoyne, J. Satayavivad, and M. Ruchirawat, Acta Tropica 89, 357 (2004)
- Various features of the invention are set forth in the following claims.
Claims (15)
1. A system for identifying a plurality of particles suspended in a solution, the system comprising:
a microfluidic device to receive the solution including the plurality of particles;
a microelectrode array comprising a plurality of electrodes, the microelectrode array disposed within the microfluidic device and in vicinity of the received solution, the microelectrode array including a first plurality of electrodes to provide a first charge and a second plurality of electrodes to provide a second charge, the second charge being different than the first charge, and the second plurality of electrodes being in a nonparallel relationship with the first plurality of electrodes;
a signal generator operatively coupled to the microelectrode array;
a particle detector adjacent to the microelectrode array; and
a controller in operative communication with the signal generator and the particle detector, the controller being configured to
apply an oscillating voltage signal to the microelectrode array at a plurality of frequency levels varying continuously between a low frequency and a high frequency, the plurality of frequency levels being applied at a sweep rate, wherein the sweep rate is no more than a maximum sweep rate and is no less than a minimum sweep rate, the applying of the oscillating voltage signal to the microelectrode array resulting in a spatially non-uniform field, and
determine a distribution of the plurality of particles relative to the microelectrode array at the plurality of frequency levels varying continuously between the low frequency and the high frequency.
2. The system of claim 1 , wherein the solution has a conductivity and wherein the maximum sweep rate is a function of the conductivity.
3. The system of claim 2 , wherein the conductivity is 0.10 S/m and the maximum sweep rate is less than 0.0026 MHz/s.
4. The system of claim 2 , wherein the conductivity is 1.0 S/m and the maximum sweep rate is less than 0.0031 MHz/s.
5. The system of claim 1 , wherein the low frequency is 0.01 MHz and the high frequency is 2.0 MHz.
6. The system of claim 1 , wherein the particle detector comprises an image detector and an image analysis system and wherein the controller, to determine the distribution of the plurality of particles relative to the microelectrode array, is further configured to collect an image of the microelectrode array using the image detector and determine a spatially resolvable concentration of the plurality of particles relative to the microelectrode array using the image analysis system at each of the plurality of frequency levels.
7. The system of claim 1 , wherein, to determine the distribution of the plurality of particles relative to the microelectrode array, the controller is further configured to use an image analysis system to determine a first spatial distribution of the plurality of particles and a second spatial distribution of the plurality of particles at a second location at each of the plurality of frequency levels.
8. The system of claim 7 , wherein, to determine the first spatial distribution of the plurality of particles at a first location, the controller is further configured to use the image analysis system to determine an intensity of the plurality of particles at the first location.
9. The system of claim 7 , wherein, to determine the second spatial distribution of the plurality of particles at the second location, the controller is further configured to use the image analysis system to determine an intensity of the plurality of particles at the second location.
10. The system of claim 1 , wherein the microelectrode array is a quadrapole microelectrode array.
11. The system of claim 1 , wherein the plurality of particles comprise red blood cells.
12. The system of claim 1 , wherein the sweep rate is no less than the minimum sweep rate.
13. The system of claim 12 , wherein the minimum sweep rate is 0.0008 MHz/s.
14-30. (canceled)
31. The system of claim 1 , wherein the system comprises a handheld device, the handheld device including the microfluidic device, the microelectrode array, the signal generator, the particle detector, and the controller.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US16/007,658 US20180292352A1 (en) | 2013-10-04 | 2018-06-13 | Methods and systems for identifying a particle using dielectrophoresis |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361887178P | 2013-10-04 | 2013-10-04 | |
PCT/US2014/059332 WO2015051372A1 (en) | 2013-10-04 | 2014-10-06 | Methods and systems for identifying a particle using dielectrophoresis |
US201615027014A | 2016-04-04 | 2016-04-04 | |
US16/007,658 US20180292352A1 (en) | 2013-10-04 | 2018-06-13 | Methods and systems for identifying a particle using dielectrophoresis |
Related Parent Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/027,014 Continuation US10012613B2 (en) | 2013-10-04 | 2014-10-06 | Methods and systems for identifying a particle using dielectrophoresis |
PCT/US2014/059332 Continuation WO2015051372A1 (en) | 2013-10-04 | 2014-10-06 | Methods and systems for identifying a particle using dielectrophoresis |
Publications (1)
Publication Number | Publication Date |
---|---|
US20180292352A1 true US20180292352A1 (en) | 2018-10-11 |
Family
ID=52779226
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/027,014 Active 2034-12-04 US10012613B2 (en) | 2013-10-04 | 2014-10-06 | Methods and systems for identifying a particle using dielectrophoresis |
US16/007,658 Abandoned US20180292352A1 (en) | 2013-10-04 | 2018-06-13 | Methods and systems for identifying a particle using dielectrophoresis |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/027,014 Active 2034-12-04 US10012613B2 (en) | 2013-10-04 | 2014-10-06 | Methods and systems for identifying a particle using dielectrophoresis |
Country Status (4)
Country | Link |
---|---|
US (2) | US10012613B2 (en) |
EP (1) | EP3052237A4 (en) |
CA (1) | CA2925112A1 (en) |
WO (1) | WO2015051372A1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10444184B2 (en) * | 2015-12-28 | 2019-10-15 | International Business Machines Corporation | Operation of diagnostic devices involving microchannels and electrodes |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5122805A (en) * | 1991-02-06 | 1992-06-16 | Radian Corporation | Radio acoustic sounding system for remotely determining atmospheric temperature profiles |
GB9615775D0 (en) | 1996-07-26 | 1996-09-04 | British Tech Group | Apparatus and method for characterising particles using dielectrophoresis |
US6294063B1 (en) * | 1999-02-12 | 2001-09-25 | Board Of Regents, The University Of Texas System | Method and apparatus for programmable fluidic processing |
US6537433B1 (en) | 2000-03-10 | 2003-03-25 | Applera Corporation | Methods and apparatus for the location and concentration of polar analytes using an alternating electric field |
US7063777B2 (en) * | 2002-12-12 | 2006-06-20 | Aura Biosystems Inc. | Dielectrophoretic particle profiling system and method |
WO2004059290A2 (en) | 2002-12-20 | 2004-07-15 | Board Of Regents, The University Of Texas System | Methods and apparatus for electrosmear analysis |
GB0516783D0 (en) | 2005-08-16 | 2005-09-21 | Univ Surrey | Micro-electrode device for dielectrophoretic characterisation of particles |
EP2259045A1 (en) * | 2009-06-05 | 2010-12-08 | Koninklijke Philips Electronics N.V. | Multi-frequency impedance method and apparatus for discriminating and counting particles expressing a specific marker |
WO2013116311A1 (en) | 2012-01-31 | 2013-08-08 | The Penn State Research Foundation | Microfluidic manipulation and sorting of particles using tunable standing surface acoustic wave |
-
2014
- 2014-10-06 US US15/027,014 patent/US10012613B2/en active Active
- 2014-10-06 WO PCT/US2014/059332 patent/WO2015051372A1/en active Application Filing
- 2014-10-06 CA CA2925112A patent/CA2925112A1/en not_active Abandoned
- 2014-10-06 EP EP14850207.3A patent/EP3052237A4/en not_active Withdrawn
-
2018
- 2018-06-13 US US16/007,658 patent/US20180292352A1/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
EP3052237A1 (en) | 2016-08-10 |
WO2015051372A1 (en) | 2015-04-09 |
US20160238558A1 (en) | 2016-08-18 |
US10012613B2 (en) | 2018-07-03 |
CA2925112A1 (en) | 2015-04-09 |
EP3052237A4 (en) | 2017-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6264815B1 (en) | Apparatus and method for testing using dielectrophoresis | |
Sun et al. | Simultaneous on-chip DC dielectrophoretic cell separation and quantitative separation performance characterization | |
Antunez-Vela et al. | Simultaneous determination of linear and nonlinear electrophoretic mobilities of cells and microparticles | |
Vahey et al. | High-throughput cell and particle characterization using isodielectric separation | |
US7063777B2 (en) | Dielectrophoretic particle profiling system and method | |
US8524063B2 (en) | Micro-electrode device for dielectrophoretic characterisation of particles | |
Zhang et al. | Characterization of single-cell biophysical properties and cell type classification using dielectrophoresis model reduction method | |
US20190032103A1 (en) | Analysis device and separation device | |
Boymelgreen et al. | Observing electrokinetic Janus particle–channel wall interaction using microparticle image velocimetry | |
Adams et al. | Frequency sweep rate dependence on the dielectrophoretic response of polystyrene beads and red blood cells | |
Cao et al. | Measuring nanoparticle polarizability using fluorescence microscopy | |
US8262886B2 (en) | Apparatus for analyzing characteristics of particulate with dielectrophoresis of particulate by applying angle-modulated wave and method for the same | |
KR100738071B1 (en) | A dielectrophoresis apparatus disposed of means for concentration gradient generation, method for separating a material and method for screening a suitable conditions for separating a material | |
Buyong et al. | Design, fabrication and characterization of dielectrophoretic microelectrode array for particle capture | |
Oh et al. | Selective manipulation of biomolecules with insulator-based dielectrophoretic tweezers | |
US20180292352A1 (en) | Methods and systems for identifying a particle using dielectrophoresis | |
Gao et al. | Insights into surface charge of single particles at the orifice of a nanopipette | |
Gielen et al. | Dielectric cell response in highly conductive buffers | |
Fernández-Mateo et al. | Concentration–polarization electroosmosis near insulating constrictions within microfluidic channels | |
Pethig et al. | Cell physiometry tools based on dielectrophoresis | |
Abd Samad et al. | Dielectrophoresis velocities response on tapered electrode profile: Simulation and experimental | |
Gimsa | Characterization of particles and biological cells by AC electrokinetics | |
Adams | Methods and systems for identifying a particle using dielectrophoresis | |
JP2006510020A (en) | Electrophoretic particle profiling system and method | |
Tsukahara et al. | Dielectrophoresis of microbioparticles in water with planar and capillary quadrupole electrodes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |