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

    Intracellular Strain Evaluation-Based Oocyte Enucleation and Its Application in Robotic Cloning

    2023-11-14 07:44:44MingZhuSunYoWeiLiuMoShengCuiQiLiZhoXingFeiZhoYiDiZhngJingJingHungGuiZhngLuXinZho
    Engineering 2023年5期

    Ming-Zhu Sun, Yo-Wei Liu, Mo-Sheng Cui*, Qi-Li Zho, Xing-Fei Zho, Yi-Di Zhng,Jing-Jing Hung Gui-Zhng Lu, Xin Zho,*

    a Institute of Robotics and Automatic Information System, Tianjin Key Laboratory of Intelligent Robotics, Nankai University, Tianjin 300071, China

    b Tianjin Animal Science and Veterinary Research Institute, Tianjin 300201, China

    c Institute of Intelligence Technology and Robotic Systems, Shenzhen Research Institute of Nankai University, Shenzhen 518063, China

    Keywords:Oocyte enucleation Robotic cloning Intracellular strain Intracellular damage Optical flow

    ARTICLEINFO Since the first cloned sheep was produced in 1996,cloning has attracted considerable attention because of its great potential in animal breeding.Somatic cell nuclear transfer (SCNT)is widely used for creating clones.However,SCNT is very complicated to manipulate and inevitably causes intracellular damage during manipulation.Typically,only less than 1%of reconstructed embryos develop into live cloned animals.This low success rate is considered to be the major limitation in the extensive application of cloning techniques.In this study,we proposed an intracellular strain evaluation-based oocyte enucleation method to reduce potential intracellular damage in SCNT.We first calculated the intracellular strain based on the intracellular velocity field and then used the intracellular strain as a criterion to improve the enucleation operation.We then developed a robotic batch SCNT system to apply this micromanipulation method to animal cloning.Experimental results showed that we increased the blastocyst rate from 10.0% to 20.8%,and we successfully produced 17 cloned piglets by robotic SCNT for the first time.The success rate of cloning was significantly increased compared to that of traditional methods(2.50%vs 0.73%on average).In addition to the cloning technique, the intracellular strain evaluation-based enucleation method is expected to be applicable to other biological operations and for establishing a universal cell manipulation protocol to reduce intracellular damage.

    1.Introduction

    Cloning has attracted great attention since Dolly,the first cloned sheep, was produced in 1996 [1] because of its great potential in animal breeding [2–7].Somatic cell nuclear transfer (SCNT) is a technique for creating clones that involves multiple manipulation procedures, such as oocyte holding; oocyte rotation, penetration,and enucleation; somatic cell aspiration; and injection [8].However, the efficiency of SCNT is disappointingly low.In Ref.[9],researchers investigated more than 30 experiments on porcine SCNT from 2000 to 2006.The success rate of cloning(number of live offspring per number of embryos transferred)varied;with the average success rate only being 0.73%.For example,in Ref.[10],two sets of donor cells were used to produce 308 nuclear transfer embryos,which were surgically transferred to three surrogate gilts.Five female piglets were delivered by caesarian from one surrogate,resulting in a cloning success rate of 1.6%.In Ref.[11],2918 nuclear transfer embryos were transferred to 20 recipient gilts, of which three were carried to term,resulting in six live births with a success rate of only 0.21%.Despite recent technological progress, the efficiency has remained largely unchanged since the first cloned animal was delivered [12].This low success rate has significantly hindered the extensive application of cloning techniques.

    Many improvements in SCNT have been made by choosing different cell types [13], changing the cell cycle [14], using different activation and fusion patterns [15,16], transferring different numbers of cloned embryos[17], and so forth.However, relatively less attention has been paid to improving manipulation skills.SCNT is a type of complicated cell surgery,and the level of manipulation skill has a significant influence on cell viability [18].Thus, improving SCNT manipulation would effectively increase the success rate.

    Oocyte enucleation, in which genetic materials, including the nucleus and polar body, are extracted from the oocyte using an injection micropipette, is a key step in SCNT.Fig.1 illustrates the oocyte enucleation process.During enucleation, the cytoplasm is torn rapidly by pressure,and the oocyte loses part of its cytoplasm,leading to huge damage to the oocyte.In previous studies, we designed the path of the injection micropipette inside the oocyte,analyzed the cytoplasm loss for complete removal of the nucleus,and realized enucleation volume control based on the adaptive slide mode [19,20].These methods reduce cytoplasmic loss and increase developmental potential.However, oocyte enucleation is a complicated dynamic process, and its effect on oocyte viability remains unclear.

    During oocyte enucleation, the external force is exerted on the oocyte by the injection micropipette and injector, in which huge deformations of the oocyte are generated, resulting in high intracellular stress.Several force sensors for cell manipulation have recently been reported [21–23], but what the force sensors measured was the stress exerting on the cell.Because intracellular stress cannot be directly measured, we observed the intracellular linear strain instead, which is related to the change in the length of the deformed object[24].Some physical markers,including fluorescent microbeads [25] and magnetic beads [26], have been applied to intracellular strain measurements.We also proposed a simple and label-free measurement method based on computer vision [27,28].

    Several studies have demonstrated that mechanical strain can affect cell viability.Scott confirmed that apoptosis can occur in response to short-term and high-strain mechanical loading [29];Gladman also revealed that mechanical injury (20% tensile strain)led to significant neuronal cell death [30].In previous studies, we analyzed the dynamic process of oocyte penetration by evaluating intracellular strain.The experimental results indicate that cellular developmental potential is strongly negatively correlated to intracellular strain [27,28].Therefore, intracellular strain bridges the cell manipulation and developmental potential.It is important to reduce the intracellular strain during oocyte enucleation to reduce potential intracellular damage.

    In this study, we reduced the large intracellular strain by improving enucleation initialization and the dynamic enucleation process.We used a pneumatic injector instead of a conventional hydraulic injector, which is commonly used for manual enucleation.Oocyte enucleation by a hydraulic injector is very fast, but the movement of the cytoplasm is unstable due to oil inertia,leading to a large intracellular strain.As for the pneumatic injector,extra pressure, which is called the balance pressure, is required to balance the capillary force.The initial balance pressure affects the initial force exerted on the oocyte and further affects intracellular strain during enucleation.We evaluated the statistical data on intracellular strain during the entire enucleation process and obtained the proper enucleation initialization parameter.

    The dynamic process of cytoplasmic extraction is the most important step in oocyte enucleation.During the process,the force exerted inside the oocyte changes dynamically, leading to large intracellular strain.However, directly performing online enucleation control based on intracellular strain is unrealistic because of the large errors during strain detection.The dynamic strain inside the cell is affected by the change in the pressure of the injector as well as the movement of the cytoplasm.In this paper, we analyzed the relationship between intracellular strain and movement of the cytoplasm in the injection micropipette according to the evaluation results of intracellular strain.We then realized enucleation control via trajectory planning of the cytoplasmic movement to reduce the maximum intracellular strain during this dynamic process.

    Finally, we applied the proposed enucleation method to the entire SCNT process and implemented robotic batch SCNT.The experimental results indicated that the proposed oocyte enucleation method reduced the number of large intracellular strains.The maximum strain decreased 60.8% (from 0.79 to 0.31)compared to the manual operation.We then performed SCNT of 113 oocytes using the proposed enucleation method and manual operation.The blastocyst rate, which is the last stage of in vitro development of reconstructed embryos and denoting successful cloning, increased from 10.0% to 20.8% when using the proposed method.Finally, we performed SCNT on 525 oocytes using a robotic batch SCNT system.The average success rate and survival rate were 99.1% and 97.1% respectively.A total of 510 embryos were transferred to six recipients.Two recipients became pregnant and delivered ten and seven piglets, respectively.Compared to manual operation, the success rate of cloning increased significantly from an average of 0.73% to 2.50%.We greatly increased the developmental potential by improving the enucleation operation and reducing intracellular strain, achieving robotic cloning results superior to those of manual operation for the first time.

    2.Materials and methods

    2.1.Intracellular strain calculation in oocyte enucleation

    2.1.1.Visual detection in oocyte enucleation

    Oocyte enucleation can be divided into two steps according to cytoplasmic movement.

    Step 1:Before cytoplasmic movement.The initial balance pressure was set by adjusting the volume of culture medium in the injection micropipette.After micropipette penetration, the pressure was reduced gradually to generate aspiration force until the cytoplasm began to move in the micropipette.

    Step 2:Cytoplasmic movement.The cytoplasm was moved into the micropipette from the oocyte by continuously controlling the pressure, as shown in Fig.1(b).

    When enough cytoplasm was aspirated into the injection micropipette, which ensures that the genetic materials have been extracted, the micropipette withdrew from the oocytes separating the cytoplasm inside and outside the micropipette.Simultaneously,the pressure was increased to stop the cytoplasm from moving inside the micropipette, as shown in Fig.1(c).

    Fig.1.Oocyte enucleation.(a) Before enucleation, the oocyte is immobilized using the holding micropipette and penetrated by the injection micropipette; (b) during enucleation,the cytoplasm is extracted into the micropipette;(c)after enucleation,the injection micropipette withdraws from the oocyte,while the cytoplasm is still moving inside the micropipette.

    In this study, a series of image processing methods were combined to automatically detect the positions of the oocyte, tips of the holding micropipette and injection micropipette, and movement of the cytoplasm in the injection micropipette.

    First, we determined the position of holding micropipette region by image scanning and detected the injection micropipette tip outside the oocyte by using a template-matching algorithm.As shown in Fig.2(a),we set the position of the injection micropipette tip to the middle of the slope of the micropipette tip to facilitate subsequent strain analysis.Second, we obtained the region of interest (ROI) of the oocyte based on the position of the holding micropipette tip.To identify the cell region,we detected the oocyte edge using the Canny edge detection algorithm and found the oocyte contour after open operation of mathematical morphology,as shown in Figs.2(b)–(d).After the injection micropipette had penetrated the oocyte, we estimated its position according to the distance moved by the micropipette, as shown in Fig.2(e).

    In step 1 of oocyte enucleation, the cytoplasm started to move and entered the micropipette after an aspiration force was exerted on the oocyte for a duration.We used the frame difference method to determine whether the cytoplasm had entered the micropipette.Fig.3(a) shows the ROI for cytoplasmic movement detection.The difference image was obtained by subtracting two consecutive frames.The difference image was then transformed into a binary image using the Otsu binarization algorithm, as shown in Fig.3(b).Step 1 was completed when the area of the foreground was greater than the threshold, which was set to 200 pixels in the experiment.

    In step 2 of oocyte enucleation, the cytoplasm rapidly moved from the oocyte into the micropipette.We detected the position of the cytoplasmic interface in the micropipette in real time using the Lucas–Kanade optical flow, which is a type of sparse optical flow method [31].As shown in Fig.3(c), 100 initial tracking points were selected on the rightmost side of the foreground according to the frame difference result at the end of step 1.The average position of the points on the x-axis was calculated as the initial position of the cytoplasmic interface.Lucas–Kanade optical flow was then applied to cytoplasmic interface tracking,and the average position of the tracking points was calculated frame by frame to update the position of the cytoplasmic interface,as shown in Fig.3(d).When the cytoplasmic interface moved to the set position of the injection micropipette, that is, when sufficient cytoplasm has been aspirated, step 2 was considered complete.Video S1 in Appendix A shows the tracking results of the cytoplasmic interface.

    Fig.2.Visual detection of the micropipettes and the oocyte.(a) Micropipette detection by using image scanning and a template-matching algorithm; (b) region of interest (ROI) of oocyte image; (c) edge detection by the Canny detection algorithm and morphology open operator;(d)oocyte contour by contour detection algorithm; (e) injection micropipette positioning according to the distance moved after the micropipette penetration.

    Fig.3.Visual detection of cytoplasmic movement in oocyte enucleation.(a)ROI for cytoplasmic movement detection; (b) difference image binarization by the Otsu algorithm,step 1 was completed when the area of the foreground was greater than the threshold;(c)tracking point selection according to the frame difference result at the end of step 1;(d)cytoplasmic interface tracking based on Lucas–Kanade optical flow algorithm.

    2.1.2.Intracellular strain field calculation based on velocity field

    We used the dense optical flow algorithm proposed by Farneback [32] to detect intracellular motion and calculate the intracellular strain field.As shown in Fig.4(a), we considered a point as an intracellular point every ten pixels in the intracellular region according to the size of the cytoplasmic particles.The velocity field and oocyte deformation of all intracellular points were calculated by applying the optical flow algorithm to two consecutive frames, as shown in Figs.4(b) and (c).The average error of optical flow detection using forward–backward tracking was 0.67 pixels,which ensured the accuracy of the intracellular strain field calculation.Further,an additional experiment on stretching elastic rubber was conducted to verify the optical flow algorithm.The error was less than 1.00 pixel on average, indicating that the optical flow method can be applied to large deformed objects.(as shown in Fig.S1 and Video S2 in Appendix A).

    The intracellular strain fields ε of the oocytes were obtained using Eq.(1):

    where L is the initial distance between two adjacent points and ΔL is the distance variation after deformation.During oocyte enucleation, the cytoplasm was transferred into the micropipette and placed perpendicular to the slope of the injection micropipette tip.The deformation of the oocyte as well as the intracellular strain occurred mainly in this direction.We calculated the strain field diagonally at a 45° angle because the micropipette tip is usually opened at a 45° angle for SCNT.

    where L0is the initial distance between two adjacent intracellular points in one frame;(x(i),y(i))and(x(i+1),y(i+1))are the coordinates of the adjacent points in the next frame,which are calculated using the velocity field of the intracellular points.Fig.4(d) shows the intracellular strain field.The intracellular structure is stretched when the strain is greater than zero and compressed when the strain is less than zero.To reduce the large intracellular strain in oocyte enucleation, we selected a 3 × 3 intracellular-point region with the largest strain near the micropipette tip as the strain ROI for further analysis, as shown in the red box in Fig.4(d).

    Fig.4.Calculation of intracellular strain field based on velocity field.(a) Intracellular points (the red points) in the intracellular region; (b) intracellular velocity field (in yellow) in two consecutive frames, the lines with an arrow show the magnitude and direction of the velocities; (c) oocyte deformation in two consecutive frames;(d)intracellular strain field,red and blue represent the high and low values.A 3×3 intracellular point region with the largest strain near the micropipette tip was selected as the strain ROI (the red box) for further analysis.

    2.2.Oocyte enucleation based on intracellular strain evaluation

    In cell manipulations, a large intracellular strain means a large intracellular deformation; for example, when the values of strain distribution are greater than 0.5, the intracellular structure is stretched by more than a factor of 0.5.Moreover, when cytoplasm is extracted from the oocyte, a high proportion of the negative strain means that oscillation is generated in the oocyte and dominates the enucleation process.Large intracellular deformations and oscillations have an adverse effect on cellular developmental potential.In this study,we analyzed the relationship between enucleation and distribution of the intracellular strain.We then set the initial balance pressure and realized trajectory planning of cytoplasmic movement in the dynamic enucleation process based on intracellular strain evaluation to avoid large intracellular strain and reduce the negative strain.

    2.2.1.Enucleation initialization based on intracellular strain evaluation

    The initial balance pressure was set before the injection micropipette penetrated the oocytes.Our previous study demonstrated that the balance pressure decreases as more culture media are aspirated into the micropipette[33].Therefore,we provided different initial balance pressures for enucleation by aspirating different volumes of the culture media into the micropipette.As shown in Fig.5(a), five typical balance pressures (0.9, 0.8, 0.7, 0.6, and 0.5 psi) (1 psi = 6.895 kPa) were selected as the initial balance pressures, and Fig.5(b) shows the corresponding positions of the gas–liquid interface (GLI).

    Fifty oocytes collected from the same batch were randomly divided into five groups and used for intracellular strain evaluation.The pressure in the injection micropipette continuously decreased from the initial pressure required for enucleation.To ensure consistency,the micropipette was placed in the same position inside the oocyte, and a constant change rate of pressure was maintained.Specifically, during oocyte penetration, the injection micropipette moved horizontally at a speed of 50 m·s-1to 3/4 of the oocyte diameter to penetrate the oocyte; the injection micropipette then retracted by 1/2 the oocyte diameter to a suitable position for enucleation.In step 1 of oocyte enucleation, the pressure change rate was set to 0.5 psi·s-1.Because the pressure decreased rapidly, the cytoplasm was extracted as soon as it was possible to do so.In step 2, the pressure change rate was reduced from 0.5 to 0.1 psi·s-1linearly at the maximum acceleration,which decreased the speed of cytoplasmic movement and reduced cytoplasmic loss.

    The intracellular strain fields were calculated frame by frame.Fig.5(c) shows typical strain fields, with the largest intracellular strain occurring during enucleation, indicating that the large values of intracellular strain decrease, and the range shrinks as the initial balance pressure decreases during enucleation.

    We obtained a normalized histogram of all intracellular strain values in the strain ROI for each initial balance pressure.Each figure in Fig.5(d) shows the strain values of 10 enucleation experiments during the entire enucleation process,including steps 1 and 2.Figs.5(e) and (f) show the percentages of the large intracellular strain distributed in [0.5, 1] and the percentages of the negative intracellular strain with different initial balance pressures, respectively.When the initial balance pressure was high,most strain values were distributed in the positive half-plane of the histogram,which demonstrates that the intracellular structure was under stretch and there was enough driving force for enucleation.However,a large intracellular strain region([0.5,1])was also observed in the histogram, indicating that a large intracellular deformation had been generated.With the decrease in the initial pressure, the large strain region, as well as the large intracellular deformation, decreased, whereas the proportion of the negative strain increased, which means that the oscillation generated in the oocyte and the aspiration force were not large enough to facilitate smooth enucleation.Considering both large strain and negative strain, we set the initial balance pressure to 0.7 psi to avoid large intracellular strain and reduce the negative strain proportion.The intracellular stretch dominated the enucleation process, and the operation could be accomplished with smaller intracellular strain and less oscillation in the oocyte, which may lead to less potential intracellular damage.

    2.2.2.Trajectory planning of cytoplasmic interface in injection micropipette

    In step 2 of oocyte enucleation, the cytoplasm moves from the oocyte into the micropipette.We divided the enucleation process into multiple periods according to the acceleration variation in the cytoplasmic interface and analyzed the strain distributions in different periods (Fig.S2 in Appendix A).The large intracellular strain is positively correlated with the change in acceleration:the more the acceleration changes, the larger the strain.

    Fig.5.Intracellular strain with different initial balance pressures.(a) Five typical balance pressures (0.9, 0.8, 0.7, 0.6, and 0.5 psi) by aspirating different volume of culture media; (b) the corresponding positions of gas–liquid interface (GLI) of the initial balance pressures; (c) intracellular strain fields in the enucleation process with different initial balance pressures; (d) normalized histograms of strain values in strain ROI with different initial balance pressures; (e) percentages of the large intracellular strain(strain values distributed in [0.5, 1]) with different initial balance pressures; (f) percentages of negative intracellular strain with different initial balance pressures.1 psi = 6.895 kPa.

    We further selected an enucleation process with an initial balance pressure of 0.7 psi and analyzed the relationship between the intracellular strain and cytoplasmic interface movement.First,we obtained the maximum strain value for each time interval in a given time interval and calculated the standard deviation of the acceleration of the cytoplasmic interface.Second, line fitting was performed during cytoplasmic movement.Because there was a response time from cytoplasmic surface movement to intracellular strain generation, we should analyze the cytoplasmic surface movement corresponding to the intracellular strain with a specific time delay.The optimal delay time was selected based on the goodness of line fitting.Table 1 lists the line-fitting results for different time intervals.A clear linear relationship exists between the maximum strain and the standard deviation of the acceleration of the cytoplasmic interface.The delay time was almost fixed for different time intervals.The fitting result is better when the time interval is larger owing to the strain measurement error.We obtained similar linear fitting results with a time interval of 0.4 s for different initial balance pressures (0.9, 0.8, 0.7, and 0.6 psi),as shown in Fig.6.

    To reduce the large intracellular strain, the variation in the acceleration of the cytoplasmic interface should be limited.To this end,the acceleration of the cytoplasmic interface was set to a constant atat time t in step 2 of oocyte enucleation.We assumed that the velocity and acceleration of the cytoplasmic surface were a0and v0at the end of step 1, which were detected online in the experiment.The acceleration in step 2 increased from the initial acceleration a0to atin time T0:

    Table 1Line-fitting of the maximum intracellular strain and standard deviation of the acceleration of the cytoplasmic interface with different time intervals.

    Fig.6.Linear fitting results with a time interval of 0.4 s for different initial balance pressures.(a) 0.9 psi; (b) 0.8 psi; (c) 0.7 psi; (d) 0.6 psi.STD: standard deviation.

    where T is the time of cytoplasmic extraction.The trajectories of the velocity v and displacement S of the cytoplasmic interface can be expressed as follows:

    In these experiments, the movement of the cytoplasmic interface was detected in real time using the Lucas–Kanade optical flow and controlled based on an adaptive slide mode controller [19].Similar to manual operation, the velocity of the cytoplasm surface was not zero at the end of trajectory planning.Once the cytoplasm interface moved to the set position, the micropipette left the oocytes to stop oocyte enucleation, and the pressure was simultaneously increased.

    2.3.Robotic batch SCNT

    2.3.1.Robotic SCNT platform

    Robotic SCNT was implemented using a self-developed NKMR601 micromanipulation system,as shown in Fig.7.It consisted of the following parts: a standard inverted microscope (Olympus,BX-51, Japan) as the basic platform of the system; pair of inhouse developed motorized X–Y–Z micromanipulators for positioning the injection micropipette and holding micropipette; an in-house-developed motorized X–Y stage for positioning the petri dish; an in-house-developed pneumatic microinjector for supplying negative aspiration pressure and positive injection pressure;in-house-developed multi-axis controller for motor control; a charge-coupled device (CCD) camera (Panasonic, W-V-460, Japan)for visual detection; and host computer for microscopic image processing, pressure data acquisition, and motion control.The micromanipulators had a travel range of 50 mm × 50 mm ×50 mm, repeatability of 1 μm, and maximum speed of 1 mm·s-1.The X–Y stage has a travel range of 100 mm×100 mm,repeatability of 1 μm, and maximum speed of 2 mm·s-1.The pneumatic microinjector has a negative pressure range from -3 to 0 kPa, a resolution of 10 Pa for cell holding, a positive pressure range of 0–200 kPa,and a resolution of 10 Pa for cell injection.The injection valve was connected to a motor to allow the injection pressure to be controlled in real time.All images were captured using ×4 and×10 objective lenses (Olympus, Japan) with a resolution of 768 × 576 pixels and a frame rate of 20 Hz.In the experiment,the injection micropipette and holding micropipette were assembled on the right and left micromanipulators and connected to the injection channel and holding channel of the micro-injector,respectively.

    Fig.7.Self-developed NK-MR601 micromanipulation robotic system for SCNT.(a)Main parts of the micromanipulation system;(b)the Robotic SCNT system;(c)the platform and the micromanipulators.CCD: charge-coupled device.

    2.3.2.Robotic batch SCNT procedure

    In manual SCNT, the operator finds oocytes and donor cells using the lower objective lens and performs nuclear transfer using the higher objective lens.In the entire SCNT process, the operator converts and focuses the objective lens four times and switches to the different operation areas six times, which is not an efficient method of performing batch SCNT.In robotic SCNT, we designed a global map and obtained a global field of view of the batch oocyte so that batch SCNT could be implemented in an optimized sequence.First, we set up the positions of the five operation areas by human–computer interaction and built a global map in the lower objective lens, as shown in Fig.8(a).We then obtained the global field of view of the oocyte area by microscopic image stitching and automatically detected the global positions of the oocytes.Finally, the operational order of the oocytes was planned to achieve the shortest path without collision with the micropipettes,as shown in Fig.8(b).

    Fig.8.Global map and the global field of view in robotic batch SCNT.(a)Operation areas and switching order in batch SCNT, we set up the positions of the five operation areas by human–computer interaction and built the global map;(b) global field of view of the oocytes, the figure shows the detection results and the operational order of the oocytes.

    Furthermore,we divided the robotic SCNT procedure into multiple steps (Fig.S3 in Appendix A) and utilized the oocyte enucleation method based on intracellular strain evaluation in robotic SCNT.First, the initial balance pressure was adjusted to 0.7 psi by aspirating the culture media into the injection micropipette before the micropipette penetrated the oocyte.The pressure was then reduced until the cytoplasm began to move.Second,the cytoplasm, nucleus, and polar body were extracted from the oocyte according to the trajectory planning of the cytoplasmic interface in the injection micropipette.After a sufficient amount of cytoplasm was aspirated, the injection micropipette was withdrawn from the oocyte horizontally at a speed of 100 μm·s-1.

    3.Results

    3.1.Oocyte enucleation control results

    According to the enucleation experiment,the time of cytoplasm extraction T was set to 4.0 s,and T0was set to 0.2 s in Eqs.(3)–(5).The moving distance of the cytoplasm interface S was set to 150 pixels,considering the cytoplasm loss was set to 8%of the volume of the oocyte for complete removal of the nucleus [20].The blue line in Fig.9 shows the desired trajectory of the cytoplasmic interface.In this experiment,the initial velocity and acceleration of the cytoplasmic surface were 20 pixels·s-1and 8 pixels·s-2.The acceleration constant atin Eq.(5) was calculated as 8.8 pixels·s-2, and the corresponding standard deviation of the acceleration was 0.1.The movement of the cytoplasmic interface was detected in real time by Lucas–Kanade optical flow with 12 ms per frame.The average detection error was 0.61 pixels.The trajectory of the cytoplasmic surface was tracked using an adaptive slide mode controller[19].The red line in Fig.9 indicates the tracking results.

    Fig.9.Trajectory and tracking result of the cytoplasmic interface.

    Table 2 shows a comparison of the enucleation methods.Videos S3–S5 in Appendix A show the enucleation processes and strain distributions by manual operation, robotic enucleation method[19], and proposed enucleation method, respectively.The experimental results demonstrate that the enucleation process is more stable when the proposed method is used.The maximum intracellular strain of the proposed method was the smallest.Specifically,the maximum strain was reduced by 60.8% (from 0.79 to 0.31)compared to that in the typical manual operation with an initial balance pressure of 0.9 psi.

    3.2.Robotic batch SCNT results

    There are two main parts of robotic batch SCNT.First,the global map was built under the lower objective lens (×4), and then theoocytes were individually manipulated under the higher objective lens (×10), as shown in Video S6 in Appendix A.The process for each SCNT is as follows.Fig.10 shows the manipulation results.

    Table 2Enucleation method comparison.

    (1) Move the holding micropipette and the injection micropipette to the field of view, and then move the oocyte to the center of the field of view after oocyte relocalization.Fig.10(a) shows the localization results of the micropipettes and oocyte.

    (2)Hold the oocyte using a holding micropipette and rotate the oocyte along the Y-axis and Z-axis until the polar body points to the desired orientation.We derived the minimum force for oocyte rotation to ensure minimum deformation of the oocyte [34], and Fig.10(b) shows the oocyte holding and rotation results (polar body at 4 o’clock).

    (3)Adjust the initial balance pressure to 0.7 psi for oocyte enucleation,penetrate the oocyte using the injection micropipette,and move the micropipette tip to the same position in the oocyte;Fig.10(c) shows that the oocyte has been penetrated.

    (4) Reduce the balance pressure rapidly after the GLI is stable and until the cytoplasm starts to move, as shown in Fig.10(d).

    (5) Adjust the balance pressure and extract the cytoplasm, as well as the nucleus and the polar body,from the oocyte according to the trajectory planning result of the cytoplasmic interface in the injection micropipette, as shown in Fig.10(e).

    (6) Withdraw the injection micropipette from the oocyte.Fig.10(f) shows the end of enucleation.

    (7) Exclude the genetic materials from the injection micropipette, as shown in Fig.10(g).

    (8) Select a proper donor cell, and aspirate the donor cell into the injection micropipette.Fig.10(h) shows the donor cell in the micropipette.

    (9) Inject the donor cell into the oocyte.Fig.10(i) shows the donor cell injection result.

    Fig.10.Robotic SCNT procedure.(a) Micropipettes and oocyte localization; (b) oocyte holding and rotation to the desired orientation (polar body at 4 o’clock); (c) oocyte penetration after enucleation initialization; (d)balance pressure reducing, until the cytoplasm starts to move; (e) the cytoplasm,as well as the nucleus and the polar body,was extracted from the oocyte; (f) end of enucleation; (g) genetic material exclusion; (h) donor cell aspiration; (i) donor cell injection.

    In robotic batch SCNT, we need not convert the objective lens and halve the number of objective lens focusing using the robotic system based on the global field of view.Taking a group of 20 oocytes as an example, we spent approximately 21 s building the global map, but we reduced the objective lens conversion and objective lens focusing 40 and 20 folds, respectively.For each SCNT,it took approximately 86 s to manipulate one oocyte.Table 3 lists the time spent on the key steps.It took 1741 s to manipulate a group of 20 oocytes, which means that the operation speed was approximately 87 s per cell.

    To verify the effect of the enucleation method on the cellular developmental potential,we performed batch SCNT on 113 oocytes using the proposed method and manual operation.In the proposed method, the initial balance pressure was set to 0.7 psi.The cytoplasm, nucleus, and polar body were extracted from the oocyte according to the trajectory planning results of the cytoplasmic interface in the injection micropipette.The maximum intracellular strain produced during this process was 0.31.In manual operation,the initial balance pressure was set at 0.9 psi to ensure the GLI in the field of view, which is convenient for the operator.The maximum intracellular strain during the manual enucleation process was 0.79.

    The operated oocytes, which were restructured embryos, were activated and cultured for seven days.It is 11 out of 53 embryos that developed into blastocysts when the proposed method was applied; in contrast, 6 out of 60 embryos developed into blastocysts via the manual operation.The blastocyst rate increased from 10.0% to 20.8% when using the proposed method, demonstrating that large intracellular strain would cause damage to cell development and affect cellular developmental potential.This is a great achievement because blastocysts are the last stage of the in vitro development of reconstructed embryos and are a symbol of successful cloning.

    3.3.Pig cloning by robotic SCNT

    We performed SCNT from six groups of 525 oocytes using the robotic system, in which 520 oocytes were operated successfully and activated for transfer, and 510 oocytes developed into embryos.Table 4 shows the results of the robotic SCNT and those of pig cloning.The average success rate and survival rate were 99.1% and 97.1%, respectively.Subsequently, 510 embryos were transferred to six recipients.Two recipients were pregnant(Fig.S4 in Appendix A).Ten piglets were delivered by one recipient(No.206)on day 103 of gestation,consisting of seven healthy piglets and three weak piglets.Three weak piglets died 2–3 days later.The average birth weight of the piglets was 1.32 kg (range 1.25–1.46 kg).Seven piglets were delivered by another recipient (No.14304) on day 110 of gestation, consisting of six healthy piglets and one weak piglet.The weak piglet died 29 days later.The average birth weight was 1.62 kg (range 1.51–1.85 kg).We increased the success rate of pig cloning from 0.73% for manual operation to an average to 2.5% (13/510).Fig.11 shows 13 healthy piglets at 90 and 93 days of age.

    Table 3Operation time of robotic SCNT.

    We performed microsatellite analysis of the genomic DNA from each sample (Table S1 in Appendix A).The comparison results indicated that 13 of the nuclear transfer piglets were derived from the donor cell line B301 with 100% identity for all eight microsatellite markers.All nuclear transfer piglets were significantly different from the surrogate mothers (No.206 and No.14304).See the Materials and Methods section in Appendix A for more details.

    4.Discussion

    In this study, we analyzed the process of oocyte enucleation in detail and developed an intracellular strain evaluation-based oocyte enucleation method.Because intracellular strain bridges cell manipulation and cellular developmental potential, we improved the enucleation operation to reduce large intracellular strain, as well as the potential intracellular damage during oocyte enucleation.We evaluated the statistical data of intracellular strain during the entire process of enucleation and obtained a reasonable enucleation initialization parameter, the initial balance pressure for enucleation.Further, we analyzed the relationship between the intracellular strain and cytoplasmic movement in the dynamic process of enucleation and designed the trajectory of the cytoplasmic movement in the injection micropipette according to the evaluation of the intracellular strain.The experimental results indicate that the proposed method markedly reduces large intracellular strains.The maximum strain decreased 60.8% (from 0.79 to 0.31), compared with the manual operation with initial balance pressure of 0.9 psi.

    Fig.11.Thirteen healthy piglets produced by robotic SCNT.

    Table 4Robotic SCNT results and cloned pig results.

    Furthermore, we evaluated the developmental potential of the oocytes after enucleation.Reconstructed embryos developed into blastocysts after nuclear transfer.We increased the blastocyst rate from 10.0% to 20.8% using the proposed method, which is a great achievement because blastocyst is the symbol of successful cloning.

    Finally,we applied the oocyte enucleation method to SCNT and implemented a robotic-batch SCNT system.There 17 cloned piglets were successfully produced,and 13 clones developed to term with a cloning success rate of 2.5%,achieving operation results superior to those of manual operation for the first time.Robotic SCNT ensures consistent manipulation and operation results,from which biochemistry operation environments and conditions in the cloning process can benefit the most.

    The higher developmental potential implies that the proposed method reduces intracellular damage during the manipulation process.Several studies have evaluated intracellular damage at the cellular level.For example, cell mechanical damage caused by microinjection can be analyzed by counting the number of ruptured bonds [35] or evaluating the intracellular damage caused by chilling injury based on fluorescent staining [36].However, the mechanism of the intracellular damage caused by micromanipulation remains unclear.Our study provides an experimental platform to explore the mechanism of intracellular damage and an approach to reduce damage.

    From the perspective of the objects to be operated, we introduced intracellular strain into micromanipulation and provided an approach to achieve operation results superior to those of manual operation.This will promote cloning technology applications and many other biological operations, such as intracytoplasmic sperm injection and polar-body genome transfer.Robotic systems for biological manipulations, such as surgical robots, also cause biological damage to tissues and organs.However, it is difficult to evaluate the damage due to the complexity of the organisms.In this study, we explored the connection between robotic manipulations, intracellular strain,and potential intracellular damage at the cellular level,which also opens a door to exploring the damage during manipulations of surgical robots.

    In the future, we will further improve robotic SCNT.However,oocyte rotation and donor cell aspiration take a long time, as shown in Table 3.We will design new oocyte rotation and cell aspiration methods based on kinematic modeling and fluid modeling to improve operation efficiency.In addition, we will analyze the geometric and mechanical properties of oocytes and expand the robotic SCNT to other animals, such as mice and sheep.

    Acknowledgments

    This research was jointly supported by the National Key Research and Development Program of China (2018YFB1304905),the National Natural Science Foundation of China (NSFC)(62027812,U1813210,62003174,and 61903201),and China Postdoctoral Science Foundation (2020M680865).We thank Professor Li-Ning Sun and Professor Jin-Chang Liu for helpful discussions and suggestions.

    Compliance with ethic guidelines

    Ming-Zhu Sun, Yao-Wei Liu, Mao-Sheng Cui, Qi-Li Zhao, Xiang-Fei Zhao,Yi-Di Zhang,Jing-Jing Huang,Gui-Zhang Lu,and Xin Zhao declare that they have no conflict of interest or financial conflicts to disclose.

    Appendix A.Supplementary data

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

    国内毛片毛片毛片毛片毛片| 久久人妻av系列| 欧美日本视频| 国产三级在线视频| a级毛片a级免费在线| 久久久久国内视频| 国产成人啪精品午夜网站| 91九色精品人成在线观看| 我的老师免费观看完整版| 亚洲激情在线av| 欧美在线一区亚洲| 亚洲av电影不卡..在线观看| 一个人看的www免费观看视频| 亚洲乱码一区二区免费版| 男女床上黄色一级片免费看| 欧美色视频一区免费| 人妻久久中文字幕网| 国产伦在线观看视频一区| 日韩人妻高清精品专区| 欧美性猛交黑人性爽| 亚洲欧美日韩高清专用| 99精品在免费线老司机午夜| 精品电影一区二区在线| a级毛片a级免费在线| 成人性生交大片免费视频hd| 午夜视频国产福利| 午夜免费成人在线视频| 国内毛片毛片毛片毛片毛片| 色av中文字幕| 在线观看66精品国产| av欧美777| 久久6这里有精品| 欧美中文日本在线观看视频| 国产亚洲精品一区二区www| 伊人久久大香线蕉亚洲五| 哪里可以看免费的av片| 亚洲人与动物交配视频| 午夜日韩欧美国产| 色在线成人网| 亚洲av一区综合| 免费在线观看成人毛片| 国产亚洲精品一区二区www| 人人妻,人人澡人人爽秒播| 欧美性猛交黑人性爽| 黄色片一级片一级黄色片| 美女被艹到高潮喷水动态| 日韩欧美 国产精品| 五月玫瑰六月丁香| 日韩有码中文字幕| 中文字幕高清在线视频| 日韩欧美国产一区二区入口| 我要搜黄色片| 国产成人欧美在线观看| 欧美精品啪啪一区二区三区| 国产不卡一卡二| 日日干狠狠操夜夜爽| 首页视频小说图片口味搜索| 国产精品一区二区三区四区免费观看 | 宅男免费午夜| 别揉我奶头~嗯~啊~动态视频| 啦啦啦观看免费观看视频高清| 欧美三级亚洲精品| 国产av不卡久久| 久久久国产成人免费| 亚洲欧美日韩高清在线视频| 99精品欧美一区二区三区四区| 一边摸一边抽搐一进一小说| 午夜精品一区二区三区免费看| 亚洲 欧美 日韩 在线 免费| 美女高潮喷水抽搐中文字幕| 可以在线观看的亚洲视频| 国产精品一区二区免费欧美| 亚洲专区中文字幕在线| 香蕉丝袜av| 国产69精品久久久久777片| 欧美日韩国产亚洲二区| 美女免费视频网站| 成人无遮挡网站| 级片在线观看| 日韩免费av在线播放| 亚洲国产欧美人成| av片东京热男人的天堂| 成人永久免费在线观看视频| 久久99热这里只有精品18| 免费看十八禁软件| 99久久精品一区二区三区| 大型黄色视频在线免费观看| 亚洲性夜色夜夜综合| 97超级碰碰碰精品色视频在线观看| 久久精品国产综合久久久| 午夜福利18| 性色av乱码一区二区三区2| 特大巨黑吊av在线直播| 99久久九九国产精品国产免费| 国产精品香港三级国产av潘金莲| 久久亚洲真实| 久久久久久久亚洲中文字幕 | 国产国拍精品亚洲av在线观看 | 亚洲人成网站在线播放欧美日韩| 深夜精品福利| 非洲黑人性xxxx精品又粗又长| 韩国av一区二区三区四区| 日本与韩国留学比较| 日本免费a在线| 亚洲avbb在线观看| av在线蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 国产精品久久久久久久久免 | 日韩欧美 国产精品| 亚洲无线在线观看| 精品久久久久久成人av| 久久精品91蜜桃| 哪里可以看免费的av片| 中出人妻视频一区二区| 亚洲人成网站在线播放欧美日韩| 淫秽高清视频在线观看| 蜜桃久久精品国产亚洲av| 亚洲专区国产一区二区| 国产精品一区二区三区四区久久| bbb黄色大片| 好看av亚洲va欧美ⅴa在| 亚洲精品成人久久久久久| 国产精品 国内视频| 中文在线观看免费www的网站| 身体一侧抽搐| 我要搜黄色片| 99热这里只有是精品50| 精品人妻1区二区| 国产黄色小视频在线观看| 99热这里只有精品一区| 在线天堂最新版资源| 亚洲,欧美精品.| 欧美日韩一级在线毛片| 91九色精品人成在线观看| 中文字幕人成人乱码亚洲影| 别揉我奶头~嗯~啊~动态视频| 久久精品国产综合久久久| eeuss影院久久| 国产精品久久视频播放| 波多野结衣巨乳人妻| 美女大奶头视频| 天天躁日日操中文字幕| 99精品在免费线老司机午夜| 欧美bdsm另类| 亚洲成人久久爱视频| 特级一级黄色大片| 日韩欧美三级三区| 国产单亲对白刺激| 99视频精品全部免费 在线| 国产av不卡久久| 成年免费大片在线观看| 动漫黄色视频在线观看| 免费观看人在逋| 欧美中文综合在线视频| 伊人久久大香线蕉亚洲五| 一夜夜www| www.熟女人妻精品国产| 此物有八面人人有两片| 一边摸一边抽搐一进一小说| 国产亚洲精品一区二区www| 亚洲美女视频黄频| 99热6这里只有精品| 天天一区二区日本电影三级| 国产一区二区在线av高清观看| 国产视频内射| 亚洲 欧美 日韩 在线 免费| 亚洲av第一区精品v没综合| 国产av一区在线观看免费| 我的老师免费观看完整版| 免费看美女性在线毛片视频| 国产精品一区二区免费欧美| 日韩 欧美 亚洲 中文字幕| 十八禁人妻一区二区| 小蜜桃在线观看免费完整版高清| 亚洲国产精品sss在线观看| 18禁裸乳无遮挡免费网站照片| 99久久精品一区二区三区| 亚洲成人中文字幕在线播放| 99精品久久久久人妻精品| 好看av亚洲va欧美ⅴa在| www.色视频.com| 成人av在线播放网站| 国产伦一二天堂av在线观看| 黄色女人牲交| 亚洲一区高清亚洲精品| 亚洲国产高清在线一区二区三| 久久精品综合一区二区三区| 尤物成人国产欧美一区二区三区| 精品人妻1区二区| 色精品久久人妻99蜜桃| 窝窝影院91人妻| 97人妻精品一区二区三区麻豆| www.www免费av| 少妇的逼水好多| 欧美bdsm另类| 18禁黄网站禁片午夜丰满| 婷婷亚洲欧美| 嫩草影院精品99| 色尼玛亚洲综合影院| ponron亚洲| 日韩av在线大香蕉| 亚洲国产欧美网| 精品国产美女av久久久久小说| 少妇人妻精品综合一区二区 | 少妇人妻一区二区三区视频| 在线a可以看的网站| 国产国拍精品亚洲av在线观看 | bbb黄色大片| 亚洲在线观看片| 全区人妻精品视频| 国产一区二区激情短视频| 51国产日韩欧美| 9191精品国产免费久久| 无遮挡黄片免费观看| av视频在线观看入口| 色尼玛亚洲综合影院| 麻豆久久精品国产亚洲av| 中文字幕熟女人妻在线| 91麻豆av在线| 女同久久另类99精品国产91| 久久这里只有精品中国| 青草久久国产| 少妇的逼水好多| 丁香欧美五月| www.www免费av| 中文字幕久久专区| 精品久久久久久,| www.www免费av| 中文在线观看免费www的网站| 国产高潮美女av| 偷拍熟女少妇极品色| 老司机福利观看| 亚洲av二区三区四区| 国产高清激情床上av| 亚洲精华国产精华精| 国产在线精品亚洲第一网站| 男人舔奶头视频| 18禁在线播放成人免费| 亚洲人成伊人成综合网2020| 久久精品亚洲精品国产色婷小说| 亚洲成人精品中文字幕电影| 亚洲内射少妇av| 午夜精品久久久久久毛片777| 精品久久久久久久久久久久久| 一个人看的www免费观看视频| 国产真人三级小视频在线观看| 一个人免费在线观看电影| 免费看日本二区| 国产探花极品一区二区| 亚洲成人中文字幕在线播放| av视频在线观看入口| 午夜福利成人在线免费观看| 又黄又爽又免费观看的视频| 国内揄拍国产精品人妻在线| 禁无遮挡网站| av天堂中文字幕网| 国产精品女同一区二区软件 | 亚洲第一欧美日韩一区二区三区| 精品人妻1区二区| 一本一本综合久久| 亚洲成人免费电影在线观看| 91在线精品国自产拍蜜月 | 成人午夜高清在线视频| 国产亚洲精品一区二区www| 欧美黄色片欧美黄色片| 精品久久久久久久久久久久久| 久久国产乱子伦精品免费另类| 99精品久久久久人妻精品| 麻豆成人av在线观看| 香蕉丝袜av| 成人欧美大片| 超碰av人人做人人爽久久 | 91麻豆av在线| 国产午夜精品久久久久久一区二区三区 | 尤物成人国产欧美一区二区三区| 男人舔女人下体高潮全视频| 特大巨黑吊av在线直播| 日韩欧美精品v在线| 亚洲人成电影免费在线| 国内精品美女久久久久久| 亚洲精华国产精华精| 久久精品亚洲精品国产色婷小说| 国产精品自产拍在线观看55亚洲| 美女被艹到高潮喷水动态| 国产精品,欧美在线| 成熟少妇高潮喷水视频| 美女大奶头视频| 一级黄片播放器| 亚洲中文字幕一区二区三区有码在线看| 色播亚洲综合网| netflix在线观看网站| 美女黄网站色视频| 在线播放国产精品三级| 午夜福利欧美成人| 精品午夜福利视频在线观看一区| 国产成人系列免费观看| 久久精品国产99精品国产亚洲性色| 国产高清视频在线播放一区| 久久久色成人| 亚洲成人免费电影在线观看| 午夜老司机福利剧场| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美免费精品| 91在线观看av| 欧美性猛交╳xxx乱大交人| 丁香六月欧美| 亚洲精品国产精品久久久不卡| 最近最新中文字幕大全电影3| 亚洲片人在线观看| 国产成人啪精品午夜网站| 91av网一区二区| aaaaa片日本免费| 久久久久久九九精品二区国产| 桃红色精品国产亚洲av| 一级黄色大片毛片| 午夜老司机福利剧场| 欧美不卡视频在线免费观看| 久久久久久久精品吃奶| 成人国产综合亚洲| 亚洲av熟女| 欧美+日韩+精品| 国产精品爽爽va在线观看网站| 精品国内亚洲2022精品成人| 亚洲熟妇熟女久久| 国产精品野战在线观看| 久久久久久久午夜电影| 亚洲av中文字字幕乱码综合| 亚洲精品国产精品久久久不卡| 精品一区二区三区av网在线观看| 99热只有精品国产| 久久久精品大字幕| 国产国拍精品亚洲av在线观看 | 九九久久精品国产亚洲av麻豆| 亚洲精品456在线播放app | 欧美一区二区精品小视频在线| 给我免费播放毛片高清在线观看| 国内精品美女久久久久久| 黄片大片在线免费观看| 熟女电影av网| 美女cb高潮喷水在线观看| av天堂中文字幕网| a在线观看视频网站| 一个人看的www免费观看视频| 岛国在线免费视频观看| 久久久久国产精品人妻aⅴ院| 性色av乱码一区二区三区2| 国产精品亚洲美女久久久| 久久久精品欧美日韩精品| 看黄色毛片网站| 国产一区二区在线av高清观看| 欧美最黄视频在线播放免费| 18禁国产床啪视频网站| 男女那种视频在线观看| 深夜精品福利| 99热精品在线国产| 色综合站精品国产| 18+在线观看网站| 午夜福利免费观看在线| 99久久99久久久精品蜜桃| 国产亚洲欧美在线一区二区| 麻豆一二三区av精品| 国产免费av片在线观看野外av| 国产精品久久电影中文字幕| 免费av毛片视频| 日本在线视频免费播放| 美女黄网站色视频| 人妻丰满熟妇av一区二区三区| 日韩欧美三级三区| 免费在线观看日本一区| 欧美日韩中文字幕国产精品一区二区三区| 99久久99久久久精品蜜桃| 舔av片在线| 亚洲国产日韩欧美精品在线观看 | 变态另类成人亚洲欧美熟女| 性色avwww在线观看| 手机成人av网站| 色播亚洲综合网| 精品福利观看| 午夜免费激情av| 国产精品久久久久久久电影 | 亚洲乱码一区二区免费版| 日本黄色片子视频| 日韩 欧美 亚洲 中文字幕| 欧美另类亚洲清纯唯美| 午夜免费观看网址| 国产成人a区在线观看| 中文字幕人妻丝袜一区二区| 美女被艹到高潮喷水动态| 久久久久久国产a免费观看| 精品久久久久久久人妻蜜臀av| 亚洲国产欧美人成| 国产熟女xx| 观看美女的网站| 欧美绝顶高潮抽搐喷水| 夜夜爽天天搞| 日韩中文字幕欧美一区二区| 欧美日韩一级在线毛片| 国产三级在线视频| 精品99又大又爽又粗少妇毛片 | 丁香欧美五月| 日韩欧美在线二视频| 搡老岳熟女国产| 亚洲avbb在线观看| 日韩欧美三级三区| 天堂av国产一区二区熟女人妻| 琪琪午夜伦伦电影理论片6080| 天堂av国产一区二区熟女人妻| 一个人免费在线观看的高清视频| 两性午夜刺激爽爽歪歪视频在线观看| 老鸭窝网址在线观看| 国产视频一区二区在线看| 国产精华一区二区三区| 国内少妇人妻偷人精品xxx网站| 美女 人体艺术 gogo| 人人妻人人澡欧美一区二区| 老司机福利观看| 久久婷婷人人爽人人干人人爱| 天堂影院成人在线观看| 有码 亚洲区| 9191精品国产免费久久| 69av精品久久久久久| 国产精品1区2区在线观看.| 十八禁人妻一区二区| 两个人的视频大全免费| 亚洲人成伊人成综合网2020| 国产一区二区在线观看日韩 | 国产一区在线观看成人免费| 天堂影院成人在线观看| 国产一区二区在线av高清观看| 国产亚洲精品综合一区在线观看| av视频在线观看入口| 在线a可以看的网站| 成年人黄色毛片网站| 大型黄色视频在线免费观看| 成年女人毛片免费观看观看9| 国产高清激情床上av| xxx96com| 久久久久久大精品| 熟妇人妻久久中文字幕3abv| 91在线精品国自产拍蜜月 | 精品久久久久久,| 久久精品国产99精品国产亚洲性色| 少妇丰满av| 午夜久久久久精精品| 一二三四社区在线视频社区8| 亚洲中文日韩欧美视频| 亚洲欧美日韩卡通动漫| 动漫黄色视频在线观看| 亚洲人成伊人成综合网2020| 一a级毛片在线观看| 国产aⅴ精品一区二区三区波| 十八禁网站免费在线| 精品久久久久久久人妻蜜臀av| 超碰av人人做人人爽久久 | 国产黄色小视频在线观看| 一个人看的www免费观看视频| 禁无遮挡网站| 婷婷精品国产亚洲av在线| 成熟少妇高潮喷水视频| 亚洲精品色激情综合| 午夜免费男女啪啪视频观看 | 12—13女人毛片做爰片一| 97超级碰碰碰精品色视频在线观看| 日韩 欧美 亚洲 中文字幕| 观看免费一级毛片| 久久久久久大精品| 欧美不卡视频在线免费观看| 好男人在线观看高清免费视频| 亚洲真实伦在线观看| 99热这里只有精品一区| 老汉色av国产亚洲站长工具| 无限看片的www在线观看| 午夜福利在线观看吧| 国内精品美女久久久久久| 久久久久性生活片| 国产蜜桃级精品一区二区三区| 人妻丰满熟妇av一区二区三区| 听说在线观看完整版免费高清| 成人高潮视频无遮挡免费网站| 嫁个100分男人电影在线观看| 午夜福利高清视频| 人妻久久中文字幕网| 亚洲av美国av| 国产伦一二天堂av在线观看| 99riav亚洲国产免费| 嫩草影视91久久| 99久国产av精品| 美女黄网站色视频| 韩国av一区二区三区四区| 老鸭窝网址在线观看| 国产高潮美女av| 老司机午夜福利在线观看视频| 91麻豆精品激情在线观看国产| 可以在线观看毛片的网站| 性色av乱码一区二区三区2| 女同久久另类99精品国产91| 国产中年淑女户外野战色| avwww免费| 亚洲在线观看片| 免费看美女性在线毛片视频| 中文亚洲av片在线观看爽| 很黄的视频免费| 热99在线观看视频| av专区在线播放| 狠狠狠狠99中文字幕| 日韩中文字幕欧美一区二区| 桃色一区二区三区在线观看| 免费电影在线观看免费观看| 中文资源天堂在线| 国模一区二区三区四区视频| 久久久久亚洲av毛片大全| 99热这里只有精品一区| 女生性感内裤真人,穿戴方法视频| 少妇熟女aⅴ在线视频| 久久久色成人| 国产一区在线观看成人免费| 免费在线观看日本一区| 成人午夜高清在线视频| 国产高潮美女av| 色av中文字幕| 女生性感内裤真人,穿戴方法视频| 国内精品一区二区在线观看| 午夜精品在线福利| av国产免费在线观看| 日韩大尺度精品在线看网址| 岛国在线免费视频观看| 每晚都被弄得嗷嗷叫到高潮| 99久国产av精品| 久久精品影院6| 在线十欧美十亚洲十日本专区| 免费观看的影片在线观看| 男女午夜视频在线观看| 久久人妻av系列| 亚洲av一区综合| 国产极品精品免费视频能看的| 性欧美人与动物交配| 欧美黑人欧美精品刺激| 国产探花极品一区二区| 精品一区二区三区视频在线 | 精品国内亚洲2022精品成人| 久久6这里有精品| 老熟妇仑乱视频hdxx| 校园春色视频在线观看| 国产黄a三级三级三级人| 99久久九九国产精品国产免费| 欧美日韩福利视频一区二区| 午夜a级毛片| 国产精品野战在线观看| 色精品久久人妻99蜜桃| av中文乱码字幕在线| 真实男女啪啪啪动态图| 日本免费a在线| 三级国产精品欧美在线观看| 99久国产av精品| 国产黄a三级三级三级人| 亚洲色图av天堂| 在线免费观看不下载黄p国产 | eeuss影院久久| 国产av麻豆久久久久久久| 啪啪无遮挡十八禁网站| 国产精品影院久久| 久久欧美精品欧美久久欧美| 尤物成人国产欧美一区二区三区| 欧美性猛交╳xxx乱大交人| 久久精品91蜜桃| 99久久综合精品五月天人人| 国产一区二区在线av高清观看| h日本视频在线播放| 国产中年淑女户外野战色| 宅男免费午夜| 亚洲欧美日韩东京热| 在线观看免费午夜福利视频| 久久99热这里只有精品18| 韩国av一区二区三区四区| 给我免费播放毛片高清在线观看| 2021天堂中文幕一二区在线观| 制服人妻中文乱码| а√天堂www在线а√下载| 日韩免费av在线播放| 99久久成人亚洲精品观看| 日韩欧美免费精品| 国产高潮美女av| 亚洲第一电影网av| 免费看美女性在线毛片视频| 两个人的视频大全免费| 国产精品永久免费网站| 免费大片18禁| 国产精品综合久久久久久久免费| 国产av一区在线观看免费| 免费观看人在逋| 老熟妇仑乱视频hdxx| av黄色大香蕉| 国产黄色小视频在线观看| 日本免费一区二区三区高清不卡| 国产在视频线在精品| 一个人免费在线观看电影| 久久久久久九九精品二区国产| 久久久久精品国产欧美久久久| 一个人免费在线观看电影| 国产乱人视频| 国产高清三级在线| 两个人视频免费观看高清| 在线观看av片永久免费下载| 久久国产精品人妻蜜桃| 欧美激情在线99| 亚洲真实伦在线观看| 一级毛片女人18水好多| 18美女黄网站色大片免费观看| 亚洲自拍偷在线| av专区在线播放| 男人的好看免费观看在线视频| 午夜福利欧美成人| 韩国av一区二区三区四区| 国产淫片久久久久久久久 | 天天添夜夜摸|