• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Dynamic Deformation Measurement of an Intact Single Cell via Microfluidic Chip with Integrated Liquid Exchange

    2023-11-14 07:43:58XuDuDiChngShingoKnekoHistkMruymHirotkSugiurMsruTsujiiNouyukiUozumiFumihitoAri
    Engineering 2023年5期

    Xu Du, Di Chng, Shingo Kneko, Histk Mruym, Hirotk Sugiur, Msru Tsujii,Nouyuki Uozumi, Fumihito Ari,,*

    a Department of Micro-Nano Mechanical Science and Engineering, Nagoya University, Nagoya 464-8603, Japan

    b Department of Mechanical Engineering, The University of Tokyo, Tokyo 113-8656, Japan

    c Department of Biomolecular Engineering, Graduate School of Engineering, Tohoku University, Sendai 980-8579, Japan

    Keywords:Microfluidic chip Mechanical properties Dynamic deformation Single cell Liquid exchange

    ARTICLEINFO This paper reports a method to measure the mechanical properties of a single cell using a microfluidic chip with integrated force sensing and a liquid exchange function.A single cell is manipulated and positioned between a pushing probe and a force sensor probe using optical tweezers.These two on-chip probes were designed to capture and deform the cells.The single cell is deformed by moving the pushing probe,which is driven by an external force.The liquid–liquid interface is formed between the probes by laminar flow to change the extracellular environment.The position of the interface is shifted by controlling the injection pressure.Two positive pressures and one negative pressure are adjusted to balance the diffusion and perturbation of the flow.The mechanical properties of a single Synechocystis sp.strain PCC 6803 were measured in different osmotic concentration environments in the microfluidic chip.The liquid exchange was achieved in approximately 0.3–0.7 s,and the dynamic deformation of a single cell was revealed simultaneously.Measurements of two Young’s modulus values under alterable osmotic concentrations and the dynamic response of a single cell in osmotic shock can be collected within 30 s.Dynamic deformations of wild-type (WT) and mutant Synechocystis cells were investigated to reveal the functional mechanism of mechanosensitive (MS) channels.This system provides a novel method for monitoring the real-time mechanical dynamics of a single intact cell in response to rapid external osmotic changes;thus,it opens up novel opportunities for characterizing the accurate physiological function of MS channels in cells.

    1.Introduction

    The mechanical properties of cells are unique and important indexes that can be used to evaluate cells’ biophysical properties.For example, it has been reported that mechanical forces play an important role in the formation, stabilization, and dissociation of focal adhesions and adherens junctions [1].Other examples include the deformability of red blood cells, which decreases with stored time [2], and the stiffness of cancer cells, which is macroscopically lower than that of normal cells [3,4].In addition, the mechanical properties of cells are critical factors in cellular division, migration, and growth [5,6].

    The mechanical properties of cells can be measured with conventional methods such as atomic force microscopy (AFM) [7–10], microplate manipulation [11,12], micropipette aspiration[13,14], and optical tweezers [2,15].AFM has high sensitivity and resolution.However, it requires fixing the cells to a substrate,which is not easily done for suspended cells in general [16].In some methods, the accuracy of microplate manipulation and micropipette aspiration is lower, which is unsuitable for the measurement of cells within a few micrometers.Moreover, the probe or pipette must be set in a three-dimensional (3D) space, which makes it difficult to capture images with a microscope[17].Optical tweezers are suitable for manipulating single cells with a pico-Newton-order force, but the applied force is too weak to deform certain cells [16].

    Robot-integrated microfluidic chips are a promising tool for characterizing the mechanical properties of a single cell.Confinement of the target cell’s position in microchannels and the integration of a series of functional components on the chip enable high throughput and highly accurate manipulations and measurements [18,19].Our previous papers describe the use of a series of microfluidic chips to measure the mechanical properties of cells individually and continuously [16,17,20,21].Two on-chip probes were designed to capture and deform the cells.The pushing probe was moved by a piezo-actuator, and the reaction force was simultaneously measured using an on-chip force sensor probe composed of a hollow folded beam.To explore cellular characteristics in detail,we use phase detection of the moiré fringe to improve the measurement accuracy of movement of the probes.

    Recently, we evaluated the response of mechanosensitive (MS)channels to changes in osmotic concentration.The mechanical properties of a model cyanobacterium, Synechocystis sp.strain PCC6803 (hereafter called Synechocystis), were measured via a robot-integrated microfluidic chip [16].The diameter of a Synechocystis cell is approximately 2 μm, which is generally too small to manipulate.Members of the MS channel family in bacteria,including Synechocystis, can sense the cell membrane tension and release excess cytoplasmic solutes into the extracellular environment to protect the cell from a decrease in the extracellular osmotic concentration (i.e., osmotic down shock) [22,23].Among the several MS channels in a Synechocystis cell, a large-conductance MS channel homolog(MscL)—the plasma membrane-located MscL(slr0875),plays a crucial role in cell volume recovery after osmotic down shock [23].In 2018, we designed an opened chip to prevent the disturbance caused by bubbles in the microfluidic channel and altered the extracellular environment using a pipette [16].This approach, which helps decrease the disturbance of the extracellular environment, enabled us to stably measure the mechanical property of a single cell.However, we could only measure the deformation of a single cell after an osmotic concentration change.Therefore, we could not measure cell deformation dynamically.In addition,an opened chip exposes cells to the external environment(and thus to the risk of contamination),and high throughput is not easily implemented.

    In this paper, we propose a closed microfluidic chip integrated with a liquid exchange function.A liquid–liquid interface forms in the measurement area by laminar flow and is moved by controlling the injection pressure of the solutions.The pressure around the measurement area is adjusted by two positive pressures and one negative pressure to reduce flow perturbations.Based on this method, we measured the mechanical properties of single Synechocystis cells in environments with different osmotic concentrations.Furthermore, we analyzed the difference in the dynamic deformation processes of wild-type(WT)and mutant Synechocystis cells.Compared with the methods reported in previous research,this microfluidic chip can achieve a faster liquid exchange with low turbulence and can measure the dynamic deformation process and Young’s modulus of a single cell within 30 s.

    2.Material and methods

    2.1.Overview of the on-chip cellular measurement system

    The cell manipulation and measurement system consists of three parts (Fig.1).The mechanical control part includes the jig,piezo-actuator, and probes in the robot-integrated microfluidic chip.Two on-chip probes were designed to capture and deform cells.The pushing probe is driven by a piezo-actuator,and the reaction force is simultaneously measured using a force sensor probe composed of a hollow folded beam.A displacement measurement method, called the sampling moiré method, is applied to enhance the sensing accuracy [17].Small holes are patterned on the tips of the two probes and on the edge close to the measurement area.These small holes are used to generate moiré fringes that contribute to amplifying the displacement of the two probes,preventing the background noise of the image data from affecting the measurement.

    The cell observation and manipulation part of the system includes optical tweezers, a microscope, and two charge-coupleddevice (CCD) cameras.A single target cell in a microchannel is trapped at the focal point of the laser and transported by controlling the position of the focal point.Two high-speed cameras are integrated with this system to accurately measure the displacement of the two probes and the liquid–liquid interface, respectively.One monochrome camera with a wide field of view is used for measuring the displacement of the probes; the other red–green–blue (RGB) camera, with a high frame rate, is used to capture the distance moved by the liquid–liquid interface surrounding the cell.

    The fluid control part of the system includes a microfluidic flow control system(MFCS-FLEX;Fluigent Corporation,France),syringe pump (KDS120; KD Scientific Inc., USA), and chip microchannels.There are three inlets in the microfluidic chip; two are used to inject solutions with different osmotic concentrations, and one is a spare microchannel(Fig.1(b)).The microfluidic flow control system is connected to three inlets, which are used to apply positive pressure.The syringe pump is connected to an outlet and used to apply negative pressure to balance the pressure around the measurement area.

    The fabrication process of the microfluidic chip is based on micro-electromechanical systems (MEMS) technology (supplementary information in Appendix A).

    2.2.Cell culture and preparation

    All the Synechocystis cells used in this study were provided by Tohoku University.Since the intracellular pressure of a Synechocystis cell is adjusted by MS channels when the osmotic concentration decreases, Synechocystis WT and MscL-defective mutant (ΔmscL) cells were used.The mutant cells had the insertion of a spectinomycin resistance gene in symscL gene (one gene encoding a putative large conductance mechanosensitive channel homolog, slr0875) [23].Both cell types were cultured in a BG11 medium at a temperature of 28 °C for 6–7 days under continuous illumination (3.5 W light emitting diode (LED) white-light source).BG11 mixed with 0.5 mol·L-1sorbitol was used as a low osmotic concentration (LOC) solution, and BG11 mixed with 1 mol·L-1sorbitol was used as a high osmotic concentration(HOC) solution [23].Prior to the experiment, the cells cultured in BG11 were centrifuged at 620g (g is gravitational acceleration)for 2 min and resuspended in HOC solution.A small amount of rhodamine B was introduced into the LOC solution to distinguish between the two solutions and make it possible to observe the liquid–liquid interface clearly.In addition, we mixed the green fluorescent dye SYTO 9 and the red fluorescent dye propidium iodide (PI) into the solutions, as these dyes make it possible to remove dead or damaged cells before an experiment (supplementary information).

    2.3.Compression model and force sensor calibration

    The compression concept is shown in Fig.2(a).A single cell or bead is squeezed by driving a piezo-actuator to move the pushing probe.The sensor probe is connected to a folded beam force sensor,which can transduce the displacement to force (Fig.2(b)).Hence,the deformation and stress of the target object can be estimated by measuring the movement of the two probes.When the pushing probe is driven, the point where the sensor probe has a displacement is defined as the deformation start position.We then successively drive the pushing probe to squeeze the cell and the sensor probe.The sensor probe and the pushing probe have displacements of δsand δp, respectively.

    Fig.1.Overview of the experimental setup and the microfluidic chip.(a) Conceptual diagram of the measurement area at the microfluidic chip; (b) fabricated robotintegrated microfluidic chip; (c) microscopic image of the measurement area; (d) photograph of the experimental setup.CCD: charge-coupled-device.

    Fig.2.The compression concept, force sensor design, and calibration.(a) Concept of the measurement process; (b) design of the force sensor part; (c) an example of the theoretical value and calibration data of the force sensor (dark blue dots show data measured using the microfluidic chip after calibration; the orange curve shows the theoretical values when the Young’s modulus of polydimethylsiloxane(PDMS)is 0.98 MPa);(d)stability of displacement of the sensor probe.D0:the original diameter of the cell;δs,δp:the displacement of the sensor probe and the pushing probe, respectively; L, w,h:the length,width, and thickness of the rectangular beam,respectively (m);σ:standard deviation (SD); T: time.

    The deformation of the cell, δc, can be calculated by the difference in the displacement of two probes:

    To measure the original diameter of the target object,the pushing probe is driven to touch the tip of the sensor probe,and the difference in the displacement with and without the cell is taken as the original diameter of the cell, D0.

    Hence, the diameter of the cell, Dc, can be represented by the following equation:

    The Hertzian model can be used to describe the relationship between the deformation of a sphere and the force for nonadhesive elastic contact [24,25].This model is widely used for measuring the Young’s modulus of cells with a small deformation[26,27].The equation can be simplified as follows:

    where F represents the reaction force(N);k is the spring constant of the force sensor(N·m-1);Ecis the Young’s modulus of the cell(Pa),and ν is the Poisson’s ratio of the cell.The Poisson’s ratio is taken to be 0.5 by assuming that the cell is composed of incompressible materials.

    The spring constant of an end-loaded thin beam with a rectangular cross-section can be estimated by the following equation[28]:

    where Ebis the Young’s modulus of the beam (Pa); L, w, and h represent the length, width, and thickness of the rectangular beam,respectively(m).In this study,Eq.(4)is used to evaluate the spring constant of the force sensor beam.

    The spring constant of the force sensor beam in the microfluidic chip was calibrated using polydimethylsiloxane(PDMS)beads.The specific production process has been reported in our previous work[16].The PDMS beads and PDMS bars were fabricated under the same conditions.Hence, the Young’s moduli of the beads and bars are similar values.The Young’s modulus of the PDMS bar was approximately 0.98 MPa, as measured by a load sensor (LVS-10GA; Kyowa Electronic Instruments Co., Ltd., Japan).The PDMS beads were then measured to calibrate the spring constant of the force sensor (Fig.2(c)).The average calibrated spring constant(mean (standard deviation (SD))) was 0.0488 (0.0160) N·m-1.

    We calculate the reaction force by measuring the displacement of the sensor probe.Therefore,the resolution of the force measurement depends on the resolution of the displacement measurement.We used the sampling moiré method to enhance the resolution of the displacement measurement[17].The principle of the sampling moiré method is described below.First, grating patterns with a constant period are etched on the tips of two probes (Fig.1(c)).A grayscale image of the two probes is obtained by using the assembled CCD camera of the microscope.Next, a down sampling process is applied for every few pixels of the image, and the moiré fringe can be exposed by means of an interpolation process.Fourier transformation is used to analyze the intensity of the moiré fringe.Hence, the displacement of the probe is determined by measuring the phase shift of the moiré fringe[16,17].The measurement noise for the displacement of the sensor probe in the experimental environment was recorded for 5 s to verify the accuracy of the microfluidic chip (Fig.2(d)).The stability of the displacement achieved by 3σ is 17.7 nm(where σ is the SD),which corresponded to a force resolution of 0.864 nN.

    2.4.Control of the liquid exchange

    An important process in cell measurement is changing the extracellular environment to simultaneously expose cells to stimulation and detection.With laminar flow, the liquid environment surrounding the cells can be exchanged quickly and simply [29].Flow switching using laminar flow is often used in cell sorting,labeling, and analysis [30–34].However, it is difficult to use laminar flow in force measurements, since the pressure variation in laminar flow control will disturb the force sensor.To accurately measure the dynamic characteristics of cells, the effects of the liquid exchange should be reduced as much as possible.Decreasing the flow velocity and the pressure required for the liquid exchange can reduce these effects.However,the time required for the liquid exchange and the amount of diffusion in the laminar flow will be increased.It is important to ensure that the measured cell is completely separated from the initial solution and immersed in the exchanged reagent.

    The liquid exchange process is shown in Fig.3.In Fig.3(a),pressure is applied to remove bubbles in the microchannels and to verify that a stable liquid–liquid interface is present in the main microchannel.Subsequently, the pressure of the LOC solution is removed to cause the measurement area to be covered by the HOC solution.Hence, the cells in the microchannel are immersed in the HOC environment.In Fig.3(c),all pressure sources are closed to reduce the flow speed.The cell is trapped by optical tweezers and moved to the area between the two probes.After trapping the cell,the cell is pinched by driving the piezo-actuator.The pinch force must be carefully controlled to prevent large deformations.Otherwise, the Young’s modulus measured in the HOC condition will not be reliable.In fact, the extracellular polysaccharides of Synechocystis are soft and sticky, so it is possible to slightly compress the cell by moving the pushing probe, and then release the pushing probe after a few seconds.The cell will adhere to the tip of the probe, which prevents it from being washed away by the stream.Alternatively, it is possible to control the two probes with an appropriate distance and place the cell upstream of the probes using optical tweezers.When the pressure of the HOC solution is slightly increased to detach the cell from the optical tweezers,the cell will flow between the probes.Because the extracellular polysaccharides are soft, the flowing cell will be captured by the two probes, and the deformation of the polysaccharides will not change the stiffness of the cell.

    After the flow speed becomes stable, the recording (t0) is started.In Fig.3(d), the target cell is squeezed by driving the piezo-actuator, and the Young’s modulus of a single cell in the HOC condition can be measured.Subsequently, the pressure of the HOC solution is gradually increased, and the pressure of the LOC solution is slightly increased to maintain the interface on the left side of the cell;this keeps the cell immersed in the HOC condition.After trapping and measuring the cell in the HOC condition,the injection pressure of the LOC solution is rapidly increased (t1)to exchange the solutions.In this step, the dynamic deformation of the target cell and the process of liquid exchange can be recorded by high-speed cameras.After a few seconds,the pushing probe is released by the actuator(t2);then,the Young’s modulus of the cell in the LOC condition is measured, as shown in Fig.3(e).

    We performed a series of experiments to determine the appropriate pressure range with low perturbation.Solutions with the same osmotic concentration were injected into both microchannels to remove the effect of osmotic shock and expose the disturbance caused by the liquid exchange (Fig.4(a)).

    Fig.3.The cell manipulation and liquid exchange process, with the curves underneath showing the pressure corresponding to each step.(a) Pressure is applied to remove bubbles in the microchannels and verify that a stable liquid–liquid interface is present in the middle of the main microchannel;(b)the cells in the microchannel are immersed in HOC solution; (c)the single cell is manipulated by optical tweezers and probes; (d)the target cell is compressed, and the Young’s modulus of the cell is measured in the HOC;(e)the Young’s modulus of the target cell is measured in the LOC;(f)the cell is released.t0:the start of the recording;t1:the time at which the solutions are exchanged;t2:the time at which the pushing probe is released;P0:the initial pressure applied by the microfluidics flow control system;v0:the initial flow speed applied by the syringe pump.

    Fig.4.Pressure control and liquid exchange in the microfluidic chip.(a) Microscopic image of the measurement area.In the perturbation test, both microchannels were injected with the LOC solution.(b)The pressure range for a small disturbance is shown with a blue outline,and the blue dots are data used to determine the boundary of the range.(c)The curve represents the dynamic deformation of the cell during the perturbation measurement(blue dots),and the figures above the curve show the measurement area.

    When the liquid exchange was performed, the pressure of the LOC increased sharply to shift the interface from the left side to the right side of the cell.However, when the injection pressure of the LOC was too low to pass the interface through the cell,liquid exchange was not achieved (for no liquid exchange, see Fig.4(b)).When the pressure of the LOC is much greater than the other pressure or is greater than 6 kPa,the sharply increased pressure applies a huge impact on the sensor probe.The dynamic response of the cell deformation collected in such a case is not reliable.Therefore,the suitable pressure is limited to a small range.

    An example of the perturbation test is shown in Fig.4(c).The curve represents the dynamic deformation of the cell, and the figures above the curve show microscopic images of the measurement area.The liquid–liquid interface can be seen in the second image;the curve corresponding to this image shows the perturbation when the liquid exchange was performed,which was approximately 30 nm.Given that there was no osmotic change during this process, the diameter change (i.e., the deformation) should be caused by the disturbance of the liquid exchange.This perturbation is small enough (which will be explained in Section 3), so this experimental condition is used in the later experiments.

    In the experiments,we used a syringe pump to control the flow speed of the main microchannel at approximately 0.1 mL·h-1(20–30 mm·s-1in the measurement area).The injection pressure of the HOC inlet was 4 kPa.When the injection pressure of the LOC inlet was increased (from 0–1 to 4 kPa) to exchange the solutions, the interface shifted by approximately 40 μm.To ensure that the measured cell was completely immersed in the exchanged solution,the movement distance of the interface was made to be larger than the range of diffusion(i.e.,diffusion distance).We measured the diffusion distance by analyzing the grayscale level around the position of the interface.It was approximately 5.3 μm (supplementary information).The movement distance of the interface was far larger than the diffusion distance;hence,the extracellular solution was completely replaced.The following experiments were conducted under these conditions.

    In addition to controlling the pressures, we reduced the disturbance of the liquid exchange by means of the structural designs of the microfluidic chip (supplementary information).

    3.Results

    When driving the pushing probe to squeeze the cell, the maximum compression force was stable at approximately 0.2 μN(yùn),which deformed the cell by 20%–30%.We calculated the Young’s modulus of the cell using the Hertzian model (Eq.(3)).Only 10% deformation was used to calculate the Young’s modulus of the cell,because the Hertzian model is only reliable for a small deformation of an elastic sphere.Although a cell is not a perfectly elastic body, cells always have very complex mechanical properties, such as viscoelasticity, poroelasticity, and plasticity [35,36].We have observed that cells cannot recover from deformation after a long period of compression.However,in our previous work,we demonstrated the feasibility of using the Hertzian model to calculate the Young’s modulus of Synechocystis with a small deformation[16].A study has also reported that there is no major difference between the fitting curves using viscoelastic and poroelastic models[37].In this experiment, the compression of a single cell took several seconds, and the damping effect of the viscoelasticity could be ignored.The maximum compression force duration (i.e., the time of the liquid exchange process)was limited to less than 5 s to prevent a change in the Young’s modulus because of compression for a long duration (Fig.5).The time points t1and t2in Fig.3 are determined under these conditions.

    The two curves in Fig.5(b) represent the displacements of the pushing and sensor probes, while the increasing direction of the vertical axis represents the displacement to the right (i.e., the direction of the pushing probe).For the first 7 s,we drove the pushing probe to squeeze the cell.As mentioned above,a cell deformation of 10%was used to calculate the Young’s modulus.Hence,the data at approximately 0–5 s of Fig.5(b) can be converted to Fig.5(c) (blue coverage area on the left).The Young’s modulus of the cell in the HOC environment can be evaluated by the fitting curve in Fig.5(c).

    The deformation ratio δc/D0at the horizontal axis label of Figs.5(c)and(e)represent the ratio of the amount of deformation to the original diameter D0.From T=7 s to 11 s,the pushing probe was stopped, and the solution environment around the cell was alternated from HOC to LOC at t1.The dynamic deformation of a single cell can be observed in these 4 s.The curve in Fig.5(d)shows that the squeezed Synechocystis cell swells during the liquid exchange and shrinks slightly after the liquid exchange—a finding that is unreported in any research.In addition,the peak of the cell deformation before and after the osmotic down shock is approximately 15% (0.33 μm).However, as shown in Section 2.4,the peak of the perturbation before and after the liquid exchange is approximately 30 nm under the same experimental conditions,which is far smaller than that of the cell deformation during the osmotic shock.Hence, the realized low-disturbance liquid exchange system is feasible for cell deformation measurement.

    Fig.5.An example of the measurement result for WT Synechocystis.(a)Schematic of the position of the two probes and the deformation of a single cell.(b)The blue dots show the displacement of the sensor probe, and the green dots show the displacement of the pushing probe.(c) Force–deformation ratio for the HOC.Blue dots show the experimental data and the red line shows the fitted Hertzian model.(d)Dynamic deformation(cell size change)ratio with respect to the time of a single cell during the liquid exchange.(e) Force–deformation ratio for the LOC.Blue dots show the experimental data, and the red line shows the fitted Hertzian model.(f) The red bar represents the average Young’s modulus(mean(SD))in the HOC solution(0.51(0.32)MPa),and the blue bar represents the average Young’s modulus(mean(SD))in the LOC solution(1.13(0.37)MPa).t′1:the time at which the extracellular environment is exchanged;EHOC:the Young’s modulus of a single cell in the HOC solution;ELOC:the Young’s modulus of a single cell in the LOC solution.

    The cell size ratio Dc/D0at the vertical axis label of Fig.5(d)represents the ratio of the cell diameter after deformation to the original diameter.After the liquid exchange, the pushing probe was released;the data in approximately 15–25 s of Fig.5(b)can be converted to Fig.5(e) (blue coverage area on the right).The Young’s modulus in the LOC environment was calculated in this step.Notably, the disturbance of the liquid exchange happened when the pushing probe was stopped (at approximately 7–11 s in Fig.5(b)), which is revealed by the data used to calculate the Young’s modulus (at approximately 0–5 and 15–25 s; the blue coverage area in Fig.5(b)).Therefore,the disturbance does not affect the calculation accuracy of the Young’s modulus.In short, this study can measure the two Young’s moduli of the same cell under different osmotic concentrations,as well as the dynamic response of a single cell in osmotic shock within 30 s.

    The average value of the Young’s modulus (mean (SD)) in the LOC solution was 1.13 (0.37) MPa (95% confidence interval (CI),0.96–1.30), and the average value of the Young’s modulus (mean(SD)) in the HOC solution was 0.51 (0.32) MPa (95% CI, 0.35–0.65)(Fig.5(f)).Twenty cells were measured in each bar.Compared with the results from our previous research, the values of the Young’s modulus in the HOC solution(BG11 mixed with 1 mol·L-1sorbitol)were very close,and the value of the Young’s modulus in LOC(BG11 mixed with 0.5 mol·L-1sorbitol)was larger than that in HOC and smaller than that in pure BG11 solution[16].These comparisons elucidated that the measurement results in this experiment were reliable.

    To explore the function of MS channels in Synechocystis cells,we measured the dynamic deformation of WT and ΔmscL Synechocystis cells.Figs.6(a)and(b)show examples of the dynamic deformation of a single WT cell and a single mutant cell during the liquid exchange.As mentioned above,the cell size ratio Dc/D0in the vertical axis label of the figure represents the ratio of the cell diameter after deformation to the original diameter.The peak D1represents the maximum cell expansion ratio before and after the liquid exchange.D2represents the deformation ratio change 3 s after the liquid exchange started, and D3represents the decrease in the deformation ratio by shrinkage, which is also the difference between D1and D2.Because Synechocystis has a spherical cell shape without obvious polarity [38], the cell volume can be estimated as the volume of a sphere.As shown in Figs.6(a) and (b),the difference in D1between the WT and ΔmscL cells was not significant.However, the cell volume of the WT cells decreased faster than that of the ΔmscL cells.This finding accounts well for the function of MscL when the cells were subjected to hypoosmotic shock.The statistical result for the apparent cell volume is shown in Fig.6(c).The error bars in Figs.6(c) and (d) show the SD.The apparent cell volume of the mutant cells (95% CI, 4.37–5.22) is slightly larger than that of the WT cells (95%CI,3.75–5.41),which is consistent with the findings in previous research[39].A comparison of three typical deformation values in the dynamic deformation curve between WT and ΔmscL is shown in Fig.6(d).Over 10 cells were measured in each bar.The deformation ratio of D1and D2in WT are smaller than in ΔmscL, and the shrinkage value D3in WT is larger than that in ΔmscL.The 95% CI of the deformation ratio data is as follows: WT D1, 95% CI, 0.07–0.10; WT D2, 95% CI,-0.04–0.04; WT D3, 95% CI, 0.05–0.12; ΔmscL D1, 95% CI, 0.07–0.13; ΔmscL D2, 95% CI, 0.01–0.04; and ΔmscL D3, 95% CI, 0.04–0.11.

    Fig.6.Dynamic deformation (cell size change) of WT and ΔmscL Synechocystis cells.Examples of (a) the dynamic deformation of a single WT cell and (b) the dynamic deformation of a single mutant cell during the liquid exchange;(c)apparent cell volume of WT(orange)and mutant(green)cells;(d)comparison of D1,D2,and D3 between WT (orange) and mutant (green) cells.

    4.Discussion

    A microfluidic chip with an integrated liquid exchange function was developed and was used to measure the dynamic deformation of a single Synechocystis cell.The cell exposed an intriguing deformation when liquid was exchanged (Fig.5(d)).The deformation first reached a peak (i.e., the swelling of a single cell from t1towhere the time to the peak depended on the speed of the liquid exchange.The speed of the liquid exchange was approximately 0.3–0.7 s (the average value of the liquid exchange time (mean(SD)) in WT was 0.49 (0.09) s, and that in ΔmscL was 0.47 (0.15)s), which was fast enough to expose the dynamic deformation of the cell.The value of this peak (i.e., the cell deformation before and after the osmotic down shock) was approximately 0.3 μm.The perturbation of the liquid exchange was approximately 0.03 μm, which could be ignored in comparison with the value of the peak.After reaching the peak, the cell slightly shrank.A study has measured cell size or volume using a non-contact method and reported that Synechocystis will rapidly swell under osmotic down shock; subsequently, the size of the cell will remain unchanged[23].

    To explain this phenomenon, we provide a hypothesis here:When the osmotic concentration decreases (i.e., osmotic down shock), water enters the cell and makes it swell; the MS channel then senses the increased membrane tension force and releases the intracellular components[40,41].However,in our experiment,the cell was compressed, which increased the membrane tension force.After osmotic down shock, the MS channel still received a response indicating a high membrane tension force;hence,it continued to release the intracellular components, causing the cell to shrink.To prove this hypothesis, we compared the difference in the dynamic deformations of WT and mutant Synechocystis cells.The difference in the peak value D1was not apparent in the statistical result.However, the recovery of the ΔmscL cells to the original cell volume was slower than that of the WT cells (as indicated by the values of D2and D3),since the ΔmscL cells could not release the intracellular components using MS channels to adapt to the osmotic down shock (Fig.6).After the liquid exchange was completed, the WT cells showed obvious shrinkage compared with the ΔmscL cells, which proved that the MS channel can sense an external force and release cytoplasm to protect the cell from bursting.

    5.Conclusions

    In this work,we proposed a system that integrates the functions of cell manipulation, flow control,and cellular force measurement for evaluating the mechanical properties of a single cell.A robotintegrated microfluidic chip was fabricated using MEMS technology and combined with a liquid exchange function to alter osmotic concentrations.A liquid–liquid interface can be formed at the measurement area by laminar flow and can be shifted by adjusting the injection pressures.The low disturbance range of the injection pressures was revealed by an experiment without an osmotic concentration change.The Hertzian model was used to calculate the Young’s modulus of Synechocystis.The deformation and reaction forces of one cell under different osmotic concentrations were measured.The average value of the Young’s modulus in the LOC solution (BG11 mixed with 0.5 mol·L-1sorbitol) was (1.13 ± 0.37)MPa, and the average value of the Young’s modulus in the HOC solution (BG11 mixed with 1 mol·L-1sorbitol) was (0.51 ± 0.32)MPa.Compared with our previous research, the Young’s modulus values obtained in this study were reasonable.In addition, the dynamic deformation of a single Synechocystis cell was revealed in this study.The squeezed cell first swelled when the osmotic concentration decreased and then slightly shrank.Moreover, the dynamic deformation processes of WT and mutant Synechocystis were compared.These results confirm that MS channels work to release cytoplasm to adapt to the extracellular environment.

    Acknowledgments

    This work was partially supported by Grants-in-Aid for Scientific Research from the Ministry of Education,Culture,Sports,Science, and Technology (18H03762 and 21H04543 to Fumihito Arai and Nobuyuki Uozumi).We thank the Chinese Scholarship Council for supporting the author Du Xu in his PhD study.

    Compliance with ethics guidelines

    Xu Du, Di Chang,Shingo Kaneko,Hisataka Maruyama, Hirotaka Sugiura, Masaru Tsujii, Nobuyuki Uozumi, and Fumihito Arai declare that they have no conflict of interest or financial conflicts to disclose.

    Appendix A.Supplementary data

    Supplementary material to this article can be found online at https://doi.org/10.1016/j.eng.2022.08.020.

    九九久久精品国产亚洲av麻豆| 80岁老熟妇乱子伦牲交| 日韩亚洲欧美综合| 国产爽快片一区二区三区| a级一级毛片免费在线观看| 免费看不卡的av| 国产伦精品一区二区三区四那| 中文资源天堂在线| 精品一区二区三卡| 肉色欧美久久久久久久蜜桃 | 亚洲欧美日韩卡通动漫| 亚洲精品成人av观看孕妇| 国产一区有黄有色的免费视频| 久久鲁丝午夜福利片| 国产永久视频网站| 日韩大片免费观看网站| 偷拍熟女少妇极品色| 亚洲av成人精品一二三区| 亚洲欧美日韩无卡精品| 在线观看国产h片| 在线 av 中文字幕| 国产成人精品一,二区| 黄色配什么色好看| 国产一区二区三区av在线| 一级毛片电影观看| 亚洲国产色片| 在线观看一区二区三区激情| 亚洲成人av在线免费| 少妇人妻一区二区三区视频| 久久人人爽av亚洲精品天堂 | 夜夜爽夜夜爽视频| 国产午夜福利久久久久久| 国内揄拍国产精品人妻在线| 久久99热这里只频精品6学生| 在线免费观看不下载黄p国产| 在线观看一区二区三区激情| 免费观看av网站的网址| 18禁裸乳无遮挡动漫免费视频 | 国产高潮美女av| 精品一区二区免费观看| 午夜福利视频精品| 亚洲美女搞黄在线观看| 女人久久www免费人成看片| 成人黄色视频免费在线看| 国产高清国产精品国产三级 | 成人欧美大片| 天天躁夜夜躁狠狠久久av| 亚洲av免费在线观看| 亚洲精品久久午夜乱码| 久久精品久久久久久久性| 少妇的逼水好多| av网站免费在线观看视频| av又黄又爽大尺度在线免费看| 大片免费播放器 马上看| 国产欧美日韩一区二区三区在线 | 国语对白做爰xxxⅹ性视频网站| 久久久久久久大尺度免费视频| 国产伦理片在线播放av一区| av播播在线观看一区| 亚洲内射少妇av| 美女cb高潮喷水在线观看| 中文字幕亚洲精品专区| 不卡视频在线观看欧美| 欧美日韩精品成人综合77777| 一区二区av电影网| 美女脱内裤让男人舔精品视频| 国产午夜精品一二区理论片| 国产在线一区二区三区精| 欧美成人午夜免费资源| 日韩亚洲欧美综合| 色哟哟·www| 国产男女内射视频| 91久久精品国产一区二区三区| 综合色丁香网| 大片免费播放器 马上看| 久久久久网色| 亚洲欧美日韩另类电影网站 | 97热精品久久久久久| 免费黄色在线免费观看| 成年av动漫网址| 成人综合一区亚洲| 日韩国内少妇激情av| 成人亚洲欧美一区二区av| 国产爱豆传媒在线观看| 高清av免费在线| 99久久精品国产国产毛片| 涩涩av久久男人的天堂| 国产精品一区二区性色av| 最近最新中文字幕免费大全7| 亚洲精品影视一区二区三区av| 午夜精品国产一区二区电影 | 99久久中文字幕三级久久日本| 亚洲国产最新在线播放| 中文在线观看免费www的网站| 在线a可以看的网站| 国国产精品蜜臀av免费| 下体分泌物呈黄色| 久久亚洲国产成人精品v| 男女边摸边吃奶| 亚洲最大成人av| 国产精品99久久99久久久不卡 | 中国国产av一级| 精品一区二区三区视频在线| 国产欧美日韩精品一区二区| 亚洲av欧美aⅴ国产| 国产av不卡久久| 在线观看国产h片| 亚洲欧美一区二区三区国产| 亚洲国产精品国产精品| 成人高潮视频无遮挡免费网站| 亚洲精品日韩av片在线观看| 69av精品久久久久久| 久久久久久久大尺度免费视频| 特级一级黄色大片| 在线观看一区二区三区| 最后的刺客免费高清国语| 国产综合懂色| 久久久亚洲精品成人影院| 激情五月婷婷亚洲| 美女视频免费永久观看网站| 内射极品少妇av片p| 久久国产乱子免费精品| 久久午夜福利片| videos熟女内射| 亚洲电影在线观看av| 亚洲精品色激情综合| 亚洲最大成人手机在线| 成人一区二区视频在线观看| 精品人妻偷拍中文字幕| 99九九线精品视频在线观看视频| 超碰97精品在线观看| 女人久久www免费人成看片| 视频区图区小说| 精品久久久精品久久久| 日韩一区二区视频免费看| 日韩av不卡免费在线播放| 成人漫画全彩无遮挡| 亚洲成人一二三区av| 亚洲欧美日韩东京热| 26uuu在线亚洲综合色| 天天躁日日操中文字幕| 天堂网av新在线| 1000部很黄的大片| av在线天堂中文字幕| 男人狂女人下面高潮的视频| .国产精品久久| 日韩大片免费观看网站| 国产精品一区www在线观看| 熟女av电影| 一级黄片播放器| 如何舔出高潮| www.av在线官网国产| 精品人妻视频免费看| 91久久精品国产一区二区成人| 免费黄网站久久成人精品| 国产精品一区二区在线观看99| 亚洲av男天堂| 国产片特级美女逼逼视频| 18+在线观看网站| 97人妻精品一区二区三区麻豆| 日日摸夜夜添夜夜爱| 亚洲av在线观看美女高潮| 日本免费在线观看一区| 亚洲精品国产成人久久av| 国产亚洲一区二区精品| 三级经典国产精品| 一级a做视频免费观看| 韩国av在线不卡| 久久久久国产网址| 又粗又硬又长又爽又黄的视频| 秋霞在线观看毛片| 久久鲁丝午夜福利片| 国产成人免费观看mmmm| 黄色视频在线播放观看不卡| 91狼人影院| 免费av毛片视频| 国产片特级美女逼逼视频| 熟女av电影| 视频区图区小说| 热99国产精品久久久久久7| 人人妻人人看人人澡| av免费观看日本| 精品久久久精品久久久| 国产精品不卡视频一区二区| 少妇人妻 视频| 禁无遮挡网站| 亚洲性久久影院| 国产伦精品一区二区三区视频9| 最近手机中文字幕大全| 欧美日本视频| 国产成人免费观看mmmm| 国产精品国产三级国产专区5o| 伊人久久精品亚洲午夜| av免费在线看不卡| 日韩 亚洲 欧美在线| 亚洲最大成人中文| 国产男人的电影天堂91| 蜜桃亚洲精品一区二区三区| 国产精品一区二区在线观看99| 高清av免费在线| 色视频www国产| 亚洲av二区三区四区| 高清av免费在线| 免费看av在线观看网站| 制服丝袜香蕉在线| 久久精品久久精品一区二区三区| 精品久久久噜噜| 日本与韩国留学比较| 欧美成人一区二区免费高清观看| 日韩一区二区视频免费看| 国产一区二区三区av在线| 性色avwww在线观看| 日韩在线高清观看一区二区三区| av免费在线看不卡| 美女视频免费永久观看网站| 亚洲国产高清在线一区二区三| 我要看日韩黄色一级片| 国产 一区 欧美 日韩| 久久精品人妻少妇| videossex国产| 亚洲丝袜综合中文字幕| 搡女人真爽免费视频火全软件| 99久国产av精品国产电影| 久久亚洲国产成人精品v| www.色视频.com| 欧美xxxx性猛交bbbb| 成人高潮视频无遮挡免费网站| 国产成年人精品一区二区| 欧美xxxx性猛交bbbb| 国产精品.久久久| 大片电影免费在线观看免费| 亚洲精品国产av蜜桃| 国产成人精品久久久久久| 视频中文字幕在线观看| 国产色婷婷99| 91在线精品国自产拍蜜月| 又黄又爽又刺激的免费视频.| 只有这里有精品99| 大香蕉97超碰在线| 国产一区二区三区综合在线观看 | 色婷婷久久久亚洲欧美| 白带黄色成豆腐渣| 亚洲最大成人av| 观看免费一级毛片| 亚洲欧美日韩无卡精品| 国产黄色视频一区二区在线观看| 网址你懂的国产日韩在线| 一边亲一边摸免费视频| 国产久久久一区二区三区| 久久久亚洲精品成人影院| 性插视频无遮挡在线免费观看| 免费人成在线观看视频色| 97人妻精品一区二区三区麻豆| 男人和女人高潮做爰伦理| 日韩在线高清观看一区二区三区| 亚洲精品乱久久久久久| 狠狠精品人妻久久久久久综合| 日本黄大片高清| 亚洲色图综合在线观看| 在线观看一区二区三区激情| 五月伊人婷婷丁香| 狂野欧美激情性bbbbbb| 九九在线视频观看精品| 成人漫画全彩无遮挡| 国产日韩欧美亚洲二区| 国产精品成人在线| 男插女下体视频免费在线播放| 亚洲人成网站在线观看播放| 99久久精品热视频| 最近的中文字幕免费完整| 一个人看视频在线观看www免费| 99热网站在线观看| 亚洲精品第二区| 国产白丝娇喘喷水9色精品| 国产精品三级大全| 亚洲精品乱码久久久v下载方式| 亚洲精品亚洲一区二区| 国语对白做爰xxxⅹ性视频网站| 久久亚洲国产成人精品v| 七月丁香在线播放| 美女视频免费永久观看网站| 久久热精品热| 高清欧美精品videossex| 日韩成人伦理影院| 欧美精品国产亚洲| 国产精品蜜桃在线观看| 婷婷色麻豆天堂久久| 嫩草影院入口| 欧美+日韩+精品| 丰满人妻一区二区三区视频av| 久久人人爽人人片av| 最近中文字幕2019免费版| 成人特级av手机在线观看| 直男gayav资源| 三级国产精品欧美在线观看| av国产久精品久网站免费入址| 日韩在线高清观看一区二区三区| 亚洲国产欧美在线一区| 中文字幕免费在线视频6| 亚洲av中文av极速乱| 最后的刺客免费高清国语| 国产成人freesex在线| av又黄又爽大尺度在线免费看| 国产欧美日韩一区二区三区在线 | 亚洲久久久久久中文字幕| 精品人妻偷拍中文字幕| 国产精品熟女久久久久浪| 成人高潮视频无遮挡免费网站| 国产亚洲91精品色在线| 中文字幕亚洲精品专区| 久久国产乱子免费精品| 三级国产精品片| 色网站视频免费| 青青草视频在线视频观看| 一级毛片aaaaaa免费看小| 亚洲伊人久久精品综合| 欧美 日韩 精品 国产| 两个人的视频大全免费| av国产久精品久网站免费入址| 看黄色毛片网站| 国产精品一及| 又爽又黄a免费视频| 成人黄色视频免费在线看| 黄色视频在线播放观看不卡| 永久免费av网站大全| 精品久久久久久久末码| 日韩,欧美,国产一区二区三区| 久久精品综合一区二区三区| 日本一本二区三区精品| 亚洲伊人久久精品综合| 99久久人妻综合| 免费观看a级毛片全部| 欧美xxxx性猛交bbbb| 大话2 男鬼变身卡| 男插女下体视频免费在线播放| 白带黄色成豆腐渣| 在线观看免费高清a一片| 尾随美女入室| 国产高清国产精品国产三级 | 国产爱豆传媒在线观看| av女优亚洲男人天堂| 国产黄片美女视频| 久久6这里有精品| 中国美白少妇内射xxxbb| 日韩av免费高清视频| 国产在线一区二区三区精| 亚洲熟女精品中文字幕| 建设人人有责人人尽责人人享有的 | 亚洲国产精品国产精品| 欧美日韩在线观看h| 老司机影院成人| 成人亚洲精品一区在线观看 | 嫩草影院精品99| 99热这里只有是精品50| 成人国产麻豆网| 欧美 日韩 精品 国产| av网站免费在线观看视频| 91精品伊人久久大香线蕉| 久久久久久久久久成人| 深爱激情五月婷婷| 国产精品久久久久久精品古装| 国产精品久久久久久久电影| 亚洲va在线va天堂va国产| 精品亚洲乱码少妇综合久久| 国产成人免费观看mmmm| 成人二区视频| 国产黄色免费在线视频| 日韩欧美精品免费久久| .国产精品久久| 亚洲av一区综合| 国产精品不卡视频一区二区| 欧美另类一区| 亚洲av中文字字幕乱码综合| 欧美另类一区| 久热这里只有精品99| av免费在线看不卡| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品成人久久久久久| 麻豆精品久久久久久蜜桃| 男女边摸边吃奶| 国内精品宾馆在线| 制服丝袜香蕉在线| 天天躁日日操中文字幕| 中国美白少妇内射xxxbb| 成人漫画全彩无遮挡| 黄色欧美视频在线观看| 国产成年人精品一区二区| av黄色大香蕉| 精品久久久久久电影网| 美女视频免费永久观看网站| 只有这里有精品99| 99热全是精品| 久久99蜜桃精品久久| 老司机影院成人| www.av在线官网国产| 亚洲欧美精品专区久久| 免费电影在线观看免费观看| 亚洲天堂av无毛| 99热国产这里只有精品6| 3wmmmm亚洲av在线观看| av在线播放精品| 成年av动漫网址| 亚洲精品亚洲一区二区| 国产淫片久久久久久久久| 天堂网av新在线| 少妇丰满av| 亚洲欧美清纯卡通| 国产精品国产av在线观看| 少妇丰满av| 三级国产精品片| 99视频精品全部免费 在线| www.av在线官网国产| 亚洲综合色惰| 韩国av在线不卡| 看十八女毛片水多多多| 亚洲怡红院男人天堂| 在现免费观看毛片| 人人妻人人澡人人爽人人夜夜| 国产av不卡久久| 最近2019中文字幕mv第一页| 美女脱内裤让男人舔精品视频| av在线天堂中文字幕| 国产色爽女视频免费观看| 婷婷色麻豆天堂久久| 国产成人a∨麻豆精品| av在线app专区| 久久99蜜桃精品久久| 亚洲精品日本国产第一区| 国产成人一区二区在线| 777米奇影视久久| 欧美高清性xxxxhd video| freevideosex欧美| 男人和女人高潮做爰伦理| 久久久久国产网址| 精品久久久久久久久av| 熟妇人妻不卡中文字幕| 亚洲一区二区三区欧美精品 | 国产精品国产av在线观看| 国产精品国产三级国产av玫瑰| 国产精品久久久久久精品古装| 亚洲av欧美aⅴ国产| 欧美xxxx黑人xx丫x性爽| 亚洲欧美成人精品一区二区| 黑人高潮一二区| 亚洲欧美日韩另类电影网站 | 少妇人妻久久综合中文| 久久久a久久爽久久v久久| 各种免费的搞黄视频| 在线观看一区二区三区激情| 97在线视频观看| 日日摸夜夜添夜夜添av毛片| 亚洲aⅴ乱码一区二区在线播放| 免费av不卡在线播放| 一本一本综合久久| 免费看a级黄色片| 亚洲精品国产av蜜桃| 成人国产av品久久久| 免费av毛片视频| 男女那种视频在线观看| 久久久精品94久久精品| 久久精品国产亚洲av涩爱| 一区二区三区四区激情视频| 九九在线视频观看精品| 少妇 在线观看| 免费观看a级毛片全部| 女的被弄到高潮叫床怎么办| 在线观看一区二区三区| 精品久久久噜噜| 听说在线观看完整版免费高清| 五月开心婷婷网| 2021少妇久久久久久久久久久| av网站免费在线观看视频| 啦啦啦啦在线视频资源| 久久99热这里只有精品18| 亚洲在线观看片| 精品一区在线观看国产| 国产欧美亚洲国产| 男人和女人高潮做爰伦理| 只有这里有精品99| 日日啪夜夜爽| 一级a做视频免费观看| 内射极品少妇av片p| 亚洲国产最新在线播放| 午夜日本视频在线| 亚洲不卡免费看| 中文精品一卡2卡3卡4更新| 欧美高清成人免费视频www| 91在线精品国自产拍蜜月| 精品国产一区二区三区久久久樱花 | 一个人观看的视频www高清免费观看| 日本-黄色视频高清免费观看| 最近的中文字幕免费完整| eeuss影院久久| 一级爰片在线观看| 嫩草影院精品99| 亚洲成色77777| 亚洲精品日本国产第一区| 中文欧美无线码| av黄色大香蕉| 久久久精品欧美日韩精品| 色综合色国产| 国产精品久久久久久精品电影小说 | 日韩制服骚丝袜av| 午夜免费观看性视频| 啦啦啦中文免费视频观看日本| 99精国产麻豆久久婷婷| 少妇熟女欧美另类| 99久久精品热视频| 国产精品国产av在线观看| 在线观看免费高清a一片| 午夜福利网站1000一区二区三区| 啦啦啦啦在线视频资源| 各种免费的搞黄视频| 国产伦理片在线播放av一区| 亚洲va在线va天堂va国产| 亚洲av国产av综合av卡| 亚洲图色成人| 久久久久久国产a免费观看| 国产精品一区二区三区四区免费观看| 日韩亚洲欧美综合| 精华霜和精华液先用哪个| 亚洲在线观看片| 校园人妻丝袜中文字幕| 日本黄大片高清| 美女视频免费永久观看网站| 777米奇影视久久| 国产探花极品一区二区| 五月天丁香电影| 亚洲精品成人久久久久久| 欧美97在线视频| 91在线精品国自产拍蜜月| 身体一侧抽搐| 久久99热6这里只有精品| www.色视频.com| 97超碰精品成人国产| 18禁在线无遮挡免费观看视频| 亚洲激情五月婷婷啪啪| 伊人久久国产一区二区| 亚洲精品国产成人久久av| av福利片在线观看| 免费黄色在线免费观看| 亚洲av男天堂| 人人妻人人爽人人添夜夜欢视频 | 大片电影免费在线观看免费| 九草在线视频观看| 精品人妻熟女av久视频| 国产精品秋霞免费鲁丝片| 久久久欧美国产精品| 亚洲国产最新在线播放| 国产精品成人在线| 亚洲在线观看片| 亚洲色图综合在线观看| 国产爱豆传媒在线观看| 99久久人妻综合| 在线观看国产h片| 你懂的网址亚洲精品在线观看| 免费看av在线观看网站| 欧美精品国产亚洲| 白带黄色成豆腐渣| 成人黄色视频免费在线看| 精品一区二区三卡| 男男h啪啪无遮挡| 久久久亚洲精品成人影院| 直男gayav资源| 欧美性猛交╳xxx乱大交人| 国产熟女欧美一区二区| 久久精品久久久久久久性| 在线观看三级黄色| 成人高潮视频无遮挡免费网站| 国产人妻一区二区三区在| 久久久久精品久久久久真实原创| 亚洲人成网站在线观看播放| 一本久久精品| 亚洲成人中文字幕在线播放| 日韩一区二区视频免费看| 国产免费一区二区三区四区乱码| 2021天堂中文幕一二区在线观| 国产成人精品一,二区| 女人十人毛片免费观看3o分钟| 日韩人妻高清精品专区| 尤物成人国产欧美一区二区三区| 久久久精品欧美日韩精品| 免费观看在线日韩| 女的被弄到高潮叫床怎么办| 一级毛片电影观看| 全区人妻精品视频| 国产v大片淫在线免费观看| 夫妻午夜视频| 欧美日韩综合久久久久久| 久久99热6这里只有精品| 欧美精品国产亚洲| 亚洲成色77777| 久久精品久久久久久久性| 国产av不卡久久| 老司机影院毛片| 成年免费大片在线观看| 一边亲一边摸免费视频| av天堂中文字幕网| 成年免费大片在线观看| 亚洲内射少妇av| 成人亚洲欧美一区二区av| 亚洲精品456在线播放app| 久久精品综合一区二区三区| 日本免费在线观看一区| 国产欧美日韩一区二区三区在线 | av天堂中文字幕网| 婷婷色av中文字幕| 一级a做视频免费观看| 九九在线视频观看精品| 欧美日韩亚洲高清精品| 亚洲最大成人av| 人妻夜夜爽99麻豆av| 国产乱人偷精品视频| 免费观看性生交大片5| 久久精品人妻少妇|