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

    Integration of Droplet Microfluidic Tools for Single-cell Functional Metagenomics: An Engineering Head Start

    2021-02-24 03:06:12DavidConchousoAmaniAlMaabadiHayedehBehzadMohammedAlarawiMasahitoHosokawaYoheiNishikawaHarukoTakeyamaKatsuhikoMinetaTakashiGojobori
    Genomics,Proteomics & Bioinformatics 2021年3期

    David Conchouso, Amani Al-Ma’abadi, Hayedeh Behzad, Mohammed Alarawi, Masahito Hosokawa, Yohei Nishikawa, Haruko Takeyama,7,Katsuhiko Mineta, Takashi Gojobori,

    1Department of Industrial and Mechanical Engineering, Universidad de las Américas Puebla, Puebla 72810, Mexico

    2Computational Bioscience Research Center, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia

    3Biological and Environmental Sciences and Engineering Division,King Abdullah University of Science and Technology,Thuwal 23955-6900, Saudi Arabia

    4Research Organization for Nano & Life Innovation, Waseda University, Tokyo 162-0041, Japan

    5Department of Life Science and Medical Bioscience, Waseda University, Tokyo 162-8480, Japan

    6Institute for Advanced Research of Biosystem Dynamics, Waseda Research Institute for Science and Engineering, Waseda University,Tokyo 169-8555, Japan

    7Computational Bio Big-Data Open Innovation Laboratory, AIST-Waseda University, Tokyo 169–0072, Japan

    8Computer, Electrical, and Mathematical Sciences and Engineering Division, King Abdullah University of Science and Technology,Thuwal 23955-6900, Saudi Arabia

    Abstract Droplet microfluidic techniques have shown promising outcome to study single cells at high throughput.However, their adoption in laboratories studying “-omics” sciences is still irrelevant due to the complex and multidisciplinary nature of the field. To facilitate their use, here we provide engineering details and organized protocols for integrating three droplet-based microfluidic technologies into the metagenomic pipeline to enable functional screening of bioproducts at high throughput. First, a device encapsulating single cells in droplets at a rate of~250 Hz is described considering droplet size and cell growth. Then, we expand on previously reported fluorescence-activated droplet sorting systems to integrate the use of 4 independent fluorescence-exciting lasers (i.e., 405, 488, 561, and 637 nm) in a single platform to make it compatible with different fluorescence-emitting biosensors.For this sorter,both hardware and software are provided and optimized for effortlessly sorting droplets at 60 Hz.Then,a passive droplet merger is also integrated into our pipeline to enable adding new reagents to already-made droplets at a rate of 200 Hz.Finally,we provide an optimized recipe for manufacturing these chips using silicon dry-etching tools. Because of the overall integration and the technical details presented here,our approach allows biologists to quickly use microfluidic technologies and achieve both single-cell resolution and high-throughput capability (>50,000 cells/day)for mining and bioprospecting metagenomic data.

    KEYWORDS Droplet microfluidics;Single cell;Metagenomics;Droplet sorter;Biotechnology

    Introduction

    Droplet microfluidic technologies are revolutionary tools enabling chemical and biological discoveries [1–5]. In biology, droplet microfluidics has an extraordinary impact on developing novel assays performed at the granularity of a single cell and high throughput [6–12]. For example, it is now possible to carry out millions of single-cell studies in microfluidic droplets acting as microreactors of only~ 1–1000 pl each [13]. In these droplets, one can encapsulate single cells [6,11], insert controlled amounts of reagents [7,12–15], add unique barcodes to each droplet[7,8,16,17],and deliver biosensors to detect cell characteristics at unprecedented precision and high-throughput capability[6,7,11,18].

    Other characteristics that are key for the integration of these technologies into biological research include accurate detection of biological products[5,11,15,18–20],rapid heat transfer and cycling [21], low-cost [22,23], as well as simple integration with automation, data processing tools[6,8,18,24],and portable electronic devices[25].

    These features create new opportunities for developing ambitious protocols to mine and discover novel genes [11]and entire genomes of environmental microorganisms[26,27].The potential of droplet microfluidics is massive,as the pool of environmental microbes contains trillions of unique life forms [28,29], each containing fundamental pieces of information in their genes that are the key to understanding the mechanisms of evolution and discovering new biochemical products. These protocols can also study cell-to-cell heterogeneity between cells belonging to the same microbial population [27,30] or the same tissue or organ[8].Consequently,droplet microfluidics can open the doors to numerous unexplored areas.

    Traditional methods cannot practically handle single-cell studies,as 99%of the microorganisms cannot be cultured in the laboratory[29]and their throughput is low compared to the number of single-cell experiments needed. Some techniques, such as 16S rRNA gene sequencing, can obtain phylogenetic relations of uncultured microorganisms [31]but cannot mine genes for bioprospecting. Other metagenomic techniques enable gene mining from environmental microbiomes by cloning exogenous DNA in surrogate host cells to perform sequence-based and function-driven analyses in a well-controlled host system [29]. Nevertheless,these methods still require many manual operations,making them impractical.

    Droplet microfluidics is the only technology capable of providing the high-throughput analysis required from sample collection to sequencing.Moreover,this technology also offers seamless integration with automation[20,32,33]and digital processing tools that are vital for analyzing the Big Data produced by the gigantic pool of living microorganisms and next-generation sequencing(NGS)techniques[29,34,35].

    Despite the evidence of successful microfluidic protocols for studying single cells at unprecedented throughput and resolution [6–8,11,12], only a few laboratories working on“Lab on a Chip”use them.Some of the reasons behind this poor implementation are:1)the limited access of biologists to cleanroom facilities and the equipment needed to make microfluidic chips and conduct experiments;2)the complex and multidisciplinary nature of microfluidics (e.g., microbiology,biosensing,NGS,bioinformatics,fluid mechanics,electronics, controls, and nanofabrication); 3) the lack of standardized protocols and methodologies; and 4) the significant budget required to start a single-cell microfluidics facility.This study will focus on solving points 2)and 3)by providing critical information and engineering insights to facilitate biologists using droplet microfluidics.

    Our group studies the metagenomic diversity of the Red Sea because of its unique conditions such as high salinity and high temperature [28,30,34–37]. During this process,we assembled the set of droplet microfluidic tools, here described, to perform single-cell studies directly from metagenomic libraries of the Red Sea. Including the automatic screening of these libraries based on phenotypic responses (i.e., functional metagenomics), we selectively carried out NGS and identified new species and biological products.This study presents the integration of these tools and provides technical engineering details to implement these tools in biological laboratories. The chosen methods include the main functions needed for single-cell metagenomic studies (Figure 1): 1) single-cell encapsulation in microfluidic droplets using droplet generators with two aqueous inlets(Figure 1A)or one aqueous inlet(Figure 1B)[6,11], 2) fluorescence-based droplet sorting (Figure 1C)[6], and 3) passive droplet merging, which enables the addition of precise amounts of reagents to the droplets(Figure 1D) [12].

    Figure 1 Droplet microfluidic techniques reviewed and integrated in this studyA.Droplet generation device for single-cell encapsulation with two inlets for introducing two different sample fluids.B.Droplet generation chip with only one inlet for the sample fluid. Both droplet generators can produce droplets at rates of 250–350 droplets/s and effectively encapsulate single cells at concentration of 0.1 cells/droplet. C. Automatic droplet sorting chip based on fluorescence. The high voltage signal applied to the terminals creates a dielectrophoretic effect, resulting in droplet sorting. The electrodes conducting the signal are colored in red and black. This picture also shows the interrogation window,where fluorescence inspection is taking place by shining each droplet with a laser and capturing the emitted light with a PMT sensor.D.Passive droplet merging device capable of merging droplets at 200 droplets/s.Note.All the channels in the insets colored in blue are filled with the oil phase, whereas channels colored in orange carry the sample fluids’ aqueous phases. The designs of these devices were adapted from [6,12]. Panels A–C were adapted by permission from Mazutis and colleagues [6]. PMT, photomultiplier tube; HVA, high voltage amplifier.

    First, a guide for manufacturing these devices is provided. This guide includes recipes for producing longlasting molds on silicon using dry etching techniques.Then,single-cell encapsulation is discussed to identify critical problems arising from this process, and a guidance is provided to solve them.Later, we present technical details for building a fully automated droplet sorter using fluorescence signals, electronic controls, lasers, and optics. We also provide the program developed for this procedure,insights about using the electronic instruments, and a summary of the suggested lenses and mirrors for different fluorescenceemitting biosensors working at distinctive excitation and emission wavelengths. Additionally, a passive droplet merger is described, and details of its operation are examined. Finally, we provide a discussion on the overall integration of these methods toward functional metagenomic studies.

    This study is an engineering head start for any laboratory willing to incorporate droplet microfluidics into their research, because it gives a step-by-step guide to implement and integrate droplet microfluidic techniques into the metagenomic pipeline.

    Method

    Materials, equipment, and consumables used in this study

    Table 1shows a list of the materials and equipment used in this study, including reagents, optical components, devicefabrication materials, and microfluidic and electronic equipment. We also report on model numbers, names, and locations of the vendors to quickly obtain the exact products.A picture of the facility designed to integrate the three microfluidic technologies, here described, is shown in Figure S1.

    Table 1 Materials and equipment used in this study

    Table 2 Process parameters for DRIE etching the master molds on silicon

    (Continued)

    The MFCS series pressure controllers from Fluigent are used to drive the fluids through the microfluidic chips.These controllers have a pressure range from 0 to 700 kPa and were managed through the company’s proprietary software MAESFLO. The flow rates were measured using flow-rate sensors from the same company. Although these sensors are usually calibrated for water, it is possible to calibrate them to work with other fluids using a proportionality parameter.

    Fabrication method for manufacturing silicon master molds for soft lithography

    Various manufacturing techniques for microfluidic devices were developed using techniques borrowed from the microelectro-mechanical system (MEMS) industry. Examples include lasers and conventional milling machines[32,33,38,39], hot embossing [40], injection molding [41],digital light-projecting [42], 3D printing [43,44], capillary assembly [45], poly(ethylene)-glycol-diacrylate (PEGDA)casting [46],and paper-based patterning amongst others.Soft lithography[6,47,48]is the most popular technique for fabricating microfluidic chips. This process involves creating high-quality master molds for casting polydimethylsiloxane (PDMS) replicas. Although soft lithography was initially designed as a rapid prototyping technique for testing microfluidic designs,several protocols and commercially available products consider using PDMS chips in single-cell studies because of the following advantages:PDMS devices are affordable,the replicas inherit the high resolution of master molds, PDMS is biocompatible, and there are several available surface-crafting techniques to render the wetting properties[49–52].

    A disadvantage of using PDMS is the low reusability of the devices in biological applications. One cannot use the same chip twice due to cross-contamination between experiments.Therefore,many PDMS chips are required daily,and consequently,the master molds must last for numerous casting processes. In this situation, the master mold’s lifetime is critical. Here, we describe the most common challenges associated with the fabrication of these molds by two methods, SU-8 lithography and bulk micromachining in silicon.

    Fabrication challenges and benefits of SU-8 master molds

    SU-8 lithography is an excellent technique for making master molds for soft lithography,because SU-8 photoresist can be deposited at different thicknesses(e.g.,0.5–200 μm)and can also be patterned at high-aspect ratios (usually> 10:1).This method also allows for creating multi-height structures with a smooth surface, making it a versatile technology for creating high-quality masters[6,53].

    Although SU-8 presents excellent advantages for rapid prototyping and low-cost manufacturing of masters, depositing thick SU-8 layers (> 50 μm) can become challenging for someone unfamiliar with photolithography due to the high viscosity of the photoresist. An inexperienced person may face problems with thick SU-8 layers, such as trapped bubbles,nonuniform spin-coating,poor adhesion to the wafer, incorrect UV exposure times, internal stresses[54,55], and poor mask contact resulting in low-feature resolution due to diffraction[56].

    SU-8 masters can be used for many PDMS casting procedures, demonstrating excellent lifetime usage. However,sometimes these molds can show signs of wear, and the structures can be peeled off after some casting cycles due to any of the fabrication challenges mentioned above.

    Fabrication of silicon-based master molds using dry etching tools

    Another type of master molds,with a long-lasting lifetime,is made by bulk micromachining directly on a silicon wafer using deep reactive ion etching (DRIE) and cleanroom facilities [57,58]. Although this process requires specialized equipment and not all the laboratories have access to it,the fabrication of this type of master molds becomes less variable and more robust than SU-8.

    These molds can be used indefinitely to make PDMS chips, showing exceptional lifetime usage as features are made in bulk silicon.The DRIE process is not trivial,but we provide a recipe for Oxford Instruments’ DRIE tool with inductively coupled plasma (ICP) as a starting point for anybody who wants to manufacture these devices for the first time.

    The process begins with a new p-type {100} silicon wafer,in which an ECI-3027 photoresist mask is deposited and patterned using standard lithography. We selectively etch the silicon and create high-aspect-ratio features directly onto the silicon substrate(Figure 2).The process details are described below.

    Figure 2 Fabrication of silicon-based master molds and PDMS chipsA.SEM image of the master mold for a droplet generator device.The channel features are well defined with smooth surfaces and straight walls.B.Casting of PDMS over the master mold.We prepared a mixture of the polymeric base and the curing agent(10:1)as suggested by the Sylgard?184 Silicone-based Elastomeric Kit. This mixture was degassed in a vacuum chamber and then poured on top of the silicon master. Once the PDMS mixture is fully crosslinked, the inlets/outlets are punched, and the resulting PDMS replica is bonded to a glass slide using O2 plasma treatment. C. The inset presents an example of a finalized chip. SEM, scanning electron microscopy; PDMS, polydimethylsiloxane.

    Mask printing

    This process was performed using a μPG-101 from Heidelberg Instruments.The mask designs were modified from previous studies[6,12],and they were printed in their light field configuration for positive photoresist(i.e.,opposite to SU-8 masking). A copy of the mask files for photolithography is provided as Datasets S1–S3.

    Wafer pretreatment(optional)

    The wafer was pretreated to improve the adhesion of the photoresist by removing the existing natural oxide layer.This treatment was achieved by submerging the wafer in a 20-s buffered oxide etch (BOE) bath, rinsing with deionized (DI) water for 10 min, and drying on a hot plate. Warning, BOE contains hydrofluoric (HF) acid,and improper handling may result in serious injuries or death.A second possible pretreatment is to deposit a thin film of the hexamethyldisilizane (HMDS) adhesion promoter; however, this process requires specialized equipment [59].

    Standard photolithography

    A 4-μm layer of ECI-3027 positive photoresist was spun,prebaked, exposed, and developed to act as a mask to etch the silicon wafer selectively. Since this is a standard procedure, the resulting lithography was uniform and welldefined.

    DRIE

    This process consists of a series of etching and deposition(i.e., passivation) cycles that alternate continuously to produce high-aspect structures on silicon. A summary of the process parameters is shown inTable 2. The initial conditions for the process were set to achieve a high vacuum(e.g., 7.5 × 10?9Torr) and reach a stable temperature of?20 oC. The number of cycle repetitions depends on the feature density in our wafer. If the feature density is low,there is a larger surface area to etch away,and therefore,the etching rate is slow (e.g.,~0.328 μm/cycle for the droplet generator chip). Conversely, if the feature density is high,there is less silicon to remove, and the etching rate is fast(~ 0.4 μm/cycle for the droplet sorter chip). In any case, a feature with a height of 50 μm takes from~120–150 cycles.Figure 2A shows a scanning electron microscopy (SEM)image of the silicon mold.

    Photoresist strip

    The remaining ECI-3027 photoresist was then stripped away using acetone.

    Oxygen plasma clean In the case of any remaining photoresist left in the silicon wafer,30-s O2plasma clean is recommended.This process also helps to activate the surface in preparation for the next step.

    Pretreatment for PDMS casting

    The masters were treated with trichloro (1H,1H,2H,2Hperfluorooctyl) silane to facilitate demolding. A drop of 50–100 μl was poured into a container and placed next to the master in a closed chamber.A vacuum was then applied for 4 h to evaporate and deposit the silane on the master’s substrate. Caution, silane is toxic, and this step must be performed under a fume hood.

    PDMS casting (Figure 2B)and bonding

    After the master was completed, one can follow the standard soft lithography procedures [6,47,48].

    Fabricated devices are shown in Figure 2C. To store these chips, we recommend covering all surfaces with frosted Scotch tape until their first use to avoid capturing dust and debris on the chip’s surface.Finally,the channels’walls were rendered to favor the wetting properties of the carrier fluid. Since fluorocarbon (FC) oils such as FC-40 and HFE-7500 are commonly used in single-cell experiments,one must render the PDMS channels hydrophobic.A couple of methods to achieve this surface modification include flushing the devices with Aquapel water-repellent solution(30 s)[6,7]or with a 1%trichloro(1H,1H,2H,2Hperfluorooctyl) silane solution in HFE-7500 oil (v/v)[60,61]. In both cases, the devices must be rinsed with the oil and blow-dried immediately after treatment [11,12,62].This rinsing step is of utmost importance for Aquapel treatment because this solution can quickly react with oxygen and form particles that clog the devices.

    Single-cell encapsulation in droplet microfluidics

    Single-cell encapsulation is the basis of this research. Millions of cells can be individually analyzed in a pico-liter(pl)droplet that acts like a micro-Eppendorf tube where only the appropriate reagents and a single cell are present. Two critical factors are relevant for encapsulating only one cell per droplet:generation of droplets with stable and uniform size and a constant cell concentration in the sample fluid[6,11].Here,we describe some of the possible solutions to address these challenges.

    Microfluidic droplet generation results from the interaction between two immiscible fluids forced to flow together inside a microchannel.Droplet formation then occurs when the interfacial tension and the fluid’s viscous force are balanced [32]. Although there are different channel configurations for droplet generation [5], flow-focusing generators are frequently preferred for single-cell encapsulation [6–8,11,63]. In this configuration, the droplet size heavily depends on the flow rate ratio of two liquids.The higher the oil-water flow rate ratio is, the smaller the produced droplets [64]. The droplet’s volume is proportional to the cubic of the characteristic length;therefore,one must carefully choose a droplet size and keep it constant throughout the experiments. For example, a droplet of 100 μm in diameter has a volume of 524 pl,whereas a 40-μm droplet has only 33.5 pl. Therefore, a slight difference in droplet size can directly affect the encapsulation efficiency.

    Additionally, if droplets have different sizes, we may compromise the droplets’operational handling in the other two processes (i.e., droplet sorting or merging). Therefore,the size of the droplets must be carefully designed. In the case of our experiments, the droplet size was 45 μm in diameter.

    The stability of the droplets is addressed by modifying the surface tension between the liquids.Usually,this can be achieved using surfactants such as 2%–5% (w/w) Pico-Surf1 [11] (Sphere Fluidics Limited, Cambridge, UK) or 2%(w/w)deprotonated Krytox [7] surfactant in FC oils.

    Cell concentration can be measured optically using a bright-field microscope and a hemocytometer or electronically with a cell counter. The cell concentration must be 0.1 cells/droplet to avoid encapsulating more than one cell per droplet [11,65]. Working with living cells in nutrientrich solutions [66,67], such asE.coli-based metagenomic libraries in lysogeny broth (LB) medium, brings other complications to cell encapsulation efficiency.Cell division can occur in 20–60 min [68], and therefore, the cell concentration may not remain constant at 0.1 cells/droplet during the encapsulation process. Cell division can be delayed by lowering the temperature of the tubing in an ice bath. If cell concentration varies due to gravity, adding OptiPrep density gradient medium from Sigma-Aldrich can reduce this problem.

    Fluorescence-activated droplet sorting

    One of the biggest challenges of working simultaneously with millions of droplets (i.e., experiments) is to quickly screen positive droplets from negative ones. This task requires automation and electronics to properly handle the throughput needed[69–71].Here,we describe six essential elements to sort the droplets selectively: detection signal,control system, actuation mechanism, connection diagram and flow of information, optics and fluorescent excitation using lasers, and software development[71].

    Detection signal

    In metagenomics, fluorescent biosensors are commonly used to detect a phenotypic response in cells. These biosensors can be chosen based on particular enzymatic functions;thus,they can help us narrow down the specific cells of interest with a high probability of obtaining novel findings.In this regard,fluorescent biosensors are vital elements that enable the functional detection of targeted droplets.

    Control system

    A control system is a decision-making unit enabling time synchronization between detecting a positive droplet and switching on the actuation mechanism [71]. The system’s response time is a critical design parameter to achieve efficient sorting.The overall response time must consider the reaction time of the electronic instruments, and the time required by the actuation mechanism. Those two times should then be synchronized with the speed of the moving droplets.

    Actuation mechanism

    Several actuation mechanisms can be used to manipulate droplets in microfluidic channels. These mechanisms include magnetic field manipulation,pneumatic control, dielectrophoretic action, acoustic wave interactions, and thermally actuated mechanisms [71]. Among these techniques, dielectrophoretic action was chosen in this study,because it offers the fastest actuation time (e.g., sorting up to 30 kHz [69]) and is the best candidate for single-cell screening at high throughput [6,69,72]. Dielectrophoretic sorting consists of selectively manipulating the droplets by activating a high-voltage alternating current (AC) electric field in a localized area of a microfluidic channel.This field induces a dipole moment to slightly push a positive droplet to a different channel using a“Y-junction” [6,69,72].

    Connection diagram and flow of information

    We show the equipment and connections to successfully implement the dielectrophoretic droplet sorter in this subsection.The design of the sorter was obtained from Mazutis and colleagues[6].However,we focus on expanding some technical details because biologists are generally unfamiliar with these electronic components.Figure 3shows the flow of information from one component to the other. The process starts by detecting a fluorescence signal and ends when the high-voltage AC signal is triggered at the microfluidic chip,sorting the positive droplets out.The processing time for executing these electronic commands across the different instruments was measured to be 360 μs. The effective sorting rate achieved with this system was limited by the speed of the passing droplets but not by the electronic instrumentation.

    Optics and fluorescent excitation using lasers

    The implementation of four laser sources (i.e., 405, 488,561, and 637 nm) for exciting different fluorescence-emitting light biosensors using the same setup is presented inFigure 4. By selecting mirrors that reflect specific wavelengths and transmit the others(Table 3),one can combine the light path of four or more lasers,as shown in Figure 4A.Alternatively,these lasers can be mounted on an OBIS LX/LS Laser box,and then the lasers’light can be guided to an inverted microscope by optical fibers.

    Once excitation light arrives at the inverted microscope,it hits a filter box and then goes to the sorting chip through the objective lens. The emitted light then bounces back across the filter box and arrives at the high-speed camera and the photomultiplier tube (PMT) ports at a ratio of 90%:10%(Figure 4B).The filter box is at the center of the microscope and selectively filters and directs the excitation and emitting lights to the right locations.The filter box is a set of components composed of an excitation filter, an emission filter, and a dichroic mirror. Each of these components must be chosen depending on the excitation and emitting wavelengths of the biosensor. Fortunately, there are filter box sets that are compatible with two or more laser wavelengths. For example, the Laser Quad Band Set from Chroma (Catalog No. TRF89901-EMv2-ET-405/488/561/640nm) can be used with the four proposed laser wavelengths.

    Figure 3 Flow of information through the automatic droplet sorting chipWhen the PMT sensor detects fluorescent light,it sends a“true”signal to the FPGA board through the DAQ board.The FPGA logic board then decides to send back a triggering signal to the waveform generator (externally triggerred) also through the DAQ board. On this command, the waveform generator creates a train of 600 pulses at a frequency of 30 kHz, with an amplitude of 500 mVpp and a DC offset of 250 mV. Simultaneously, the output of the waveform generator is amplified at a ratio of 1:1000 by the HVA and sent to the microfluidic chip terminals.The reaction time for the electric signals was measured to be 360 μs. FPGA, field programable gate array; DAQ, data acquisition; DC, direct current.

    Figure 4 Lasers’ path and optical setup for fluorescence detectionA. Laser-light combiner. Dichroic mirrors are used to align the laser beams in a single path guided to the epi-fluorescent illumination port of the microscope. The mirrors are selected to pass the light wavelengths to their left and reflect the wavelength of the corresponding laser. Using this configuration,one can excite the droplets with one or more laser sources.B.Optical setup and light’s path through the inverted microscope.The excitation light travels from the input port to the sample chip,and the emitting light goes from the chip to the high-speed camera and the PMT ports.A key element reflecting and filtering the light is the filter box placed at the inverted microscope center.

    Connecting a high-speed camera and a PMT sensor at the same time is challenging.On the one hand,the high-speed camera requires an intense illumination backlight to record high-speed video. On the other hand, the PMT sensor is extremely sensitive, and it gets quickly saturated even with ambient light. For this reason, we use a beam splitter of a 90%:10%ratio to send most of the light to the camera port.Additionally,we installed a backlight filter and a PMTsingleband filter to allow only a narrow wavelength band to pass through each sensing device.These two filters are required to have independent wavelength bands, and their selection depends on the fluorescent biosensor of interest.For example,a biosensor emitting 517 nm light can use an ET510/20m as a PMT filter and an ET665lp as a backlight filter.

    Since the four laser sources are installed on the same platform,they can be used simultaneously or independently to excite the passing droplets. However, the sorting device has only a single “Y-junction”, and therefore, it can only sort droplets in two different groups:fluorescence-emitting droplets and non-fluorescence-emitting droplets. If more than two groups are required,one must repeat the process as many times as different wavelengths exist.

    Software development

    The use of field programmable gate arrays (FPGAs) has enabled the design of mid-to-high-complexity systems in a programable and manageable way [73].Since droplet sorting requires fast response times (μs range), we decided to implement the sorting controls on LabVIEW FPGA hardware(Figure 3).The code running on this device has a clock cycle of only 25 ns, thus allowing our entire program to execute in only 66 μs.This program comprises three subroutines and a user interface as depicted in Figures S2–S5.Figure S2 shows the subroutine for configuring and initializing the PMT, and Figure S3 illustrates the subroutine for reading the PMT sensor and sending a triggering signal to the function generator. Likewise, Figure S4 presents a subroutine for simulating the detection of a positive droplet,and finally,Figure S5 shows the overall user interface.

    The connection between the microfluidic system and the FPGA hardware was achieved using a data acquisition(DAQ) board (Figure 3).Table 4describes all the wiring required to run the FPGA program and implement the sorting system.

    Table 3 Optical set up

    Table 4 Electrical connections required to run the FPGA program

    The system also requires a high voltage amplifier(HVA)to produce a large electric field capable of deviating a positive droplet to a different output utilizing dielectrophoresis.

    One of the most crucial factors for successfully implementing this sorting device is to measure the time it takes for a droplet to go from the laser interrogation window to the beginning of the splitting channels (Figure 1C). This time is measured with a high-speed camera, and then it is used to synchronize the controls with the droplets’ flow effectively. A useful tool to evaluate droplets’ speed in a channel from a high-speed video was reported by Basu and colleagues [24].

    Passive droplet merger for the addition of other reagents

    The third droplet microfluidic function,reviewed here,is a passive droplet merger developed by Hosokawa and colleagues [12]. This function enables the addition of precise amounts of reagents to each droplet for further processing.For example,droplet merging has been used to add multiple displacement amplification (MDA) reagents to cell lysate droplets to amplify single-cell genomes [12]. Likewise,others have merged droplets with barcoding DNA sequences[66]to tag and then identify each experiment.In general,droplet merging is a useful tool to enable multiplestep reactions in droplet microfluidics at high throughput.

    There are active[7,65,74–77]and passive[12,14,78–80]droplet merging devices. Regardless of the merging approach, 1:1 droplet synchronization is one of the most challenging aspects. Our group has published the experimental details of the device described in this section [12].When this system runs efficiently, it can yield a droplet merging success rate of 95% at a rate as high as 12,000 droplets/min. In this merging device, partially stabilized droplets are merged with newly formed droplets of a particular reagent (e.g., lysis buffer, MDA reagents, and barcoding).Once the droplets are merged at the zigzag channel,an additional surfactant is added to stabilize the newly formed droplets.

    Droplet merging often implies reinjecting already-made droplets into a merging device.If one is not careful with this reinjection process,droplets can be poorly packed together and unevenly spaced, leading to synchronization problems when merging in a 1:1 ratio.

    For a successful reinjection,the droplet suspension must settle first until most of the droplets aggregate at the top of the container. Then these closely packed droplets can be carefully retrieved using a piece of tubing and a syringe pump with a constant withdrawing rate (e.g., 10 μl/min).Once the tubing is loaded with the droplets,one can detach the tubing with the droplet suspension from the syringe pump and connect it directly to the merging chip.

    Results and discussion

    In this section,we discuss the overall integration of the three microfluidic techniques presented here and show a practical way for achieving the high-throughput requirements of single-cell functional metagenomics studies.

    In functional metagenomics [81,82], characterization of metagenomic fragments is essential for finding applicable values as genetic resources. In practice, functional metagenomics deals with several experimental processes, such as extraction of microbial DNA from environmental samples, construction of an expression library with the microbial genetic material,and a screening method for detecting the desired phenotypic responses. The first two processes are obtained in large quantities(e.g.,millions of samples at a time), and they are discussed elsewhere [83]. However,the third process is more challenging because it involves searching for different phenotypes in each cell in a library of millions. Therefore, traditional techniques in plates cannot practically handle the number of experiments required.The required screening throughput can only be achieved using droplet-based microfluidic methods. For these reasons, in this study, we decided to integrate the three droplet microfluidic techniques to enable a practical screening of functional metagenomic libraries [6,12] in a single platform.This integration can overcome the obstacles that were not adequately addressed previously.

    For example,the protocol by Mazutis et al.[6]lacks the capability of adding precise amounts of regents to alreadyformed droplets as these droplets are closed micro-reaction chambers with no extra inputs. Similarly, the protocol by Hosokawa et al. [12] massively overcomes this limitation with a droplet merger,but unfortunately,their throughput is tremendously hindered by manually sorting the droplets with a pipette.

    Our combinatory approach of the three droplet microfluidic tools is described inFigure 5and in Movie S1.The output of the first technology is the encapsulation of millions of single cells in microdroplets (Figure 5A). This process is carried out at a cell concentration of 0.1 cells/droplet,and our device generates droplets at 250–350 droplets/s. Therefore, a hundred thousand cells can be encapsulated in approximately 1 h.

    With the second technology, one can automatically sort out the droplets that show positive fluorescence (i.e.,functional metagenomic information) from the negative ones (Figure 5B). Our sorting system offers a quick electrical response from fluorescence detection to the triggering signal, as only 360 μs is needed per event. However, the limiting factor in the sorting process is the speed at which the droplets can move through the detection zone without compromising their integrity. Therefore, our automatic droplet sorting device achieved a maximum sorting rate of 60 droplet/s.Sorting half of the droplet population produced(i.e.,~500,000 droplets or~50,000 droplets with single cells)takes roughly 2.5 h.

    Finally, the third chip enables us to combine those positive-fluorescence droplets with other reagents for further processing or amplification(Figure 5C). The droplet merging device shows a success rate of 95% at a high throughput of 200 droplets/s. Therefore, merging the positive-fluorescence droplets (e.g., sorted out from the previous device)with lysis buffer or amplification reagents can take less than 0.5 h.

    Figure 5 Droplet microfluidic capabilities enabling functional metagenomic studiesThis figure summarizes the three droplet microfluidic tools integrated in this study.A.Process diagram fully describing single-cell encapsulation.A critical step of this procedure is to measure the cell concentration before the encapsulation to guarantee that only one cell is encapsulated per droplet.B.Protocol for droplet sorting.This protocol includes selecting the right optical components and identifying possible errors on the hardware and software.Measuring the time from droplet detection to the sorting junction is fundamental for synchronizing the control system with the device’s droplet flow. C. Merging process of two droplets to add new reagents to already-made droplets.The most crucial step in this process is to adjust the inlet flow rates to obtain a 1:1 droplet synchronization.

    Once the metagenomic material is obtained from the positive droplets, the next steps consist of performing sequencing and bioinformatics.

    This protocol allows the finding of functional genes out of vast metagenomic libraries. Our group is currently working on a case study that focuses on the metagenomic application of the methods presented in this publication.In this follow-up work,we apply the techniques presented here to search for lipolytic enzymes in Red Sea environmental samples. This metagenomic investigation lead to identifying three novel lipase genes. These genes showed thermostability and high activity based on alignment tests and functional assay results.

    Conclusion

    A recipe for bulk micromachining high-quality and longlasting master molds on silicon is provided for the first time.Although our procedure requires specialized DRIE equipment, it shows a more robust solution to mold manufacturing than traditional master mold techniques such as SU-8 lithography because these molds can be used indefinitely(e.g.,>100 times)without showing signs of wear.

    Three microfluidic droplet techniques for single-cell studies are also described:single-cell encapsulation,fluorescence-based automatic droplet sorting, and droplet merging. For encapsulation, cell concentration is the most crucial parameter,and therefore it must be measured before starting any experiment.A concentration of 0.1 cells/droplet provides the best results to obtain single-cell encapsulation.Our droplet encapsulation system generates droplets at rates of~250–350 droplets/s and effectively encapsulates approximately 25 droplets/s.Automatic droplet sorting is then described entirely, providing guides to optics, instrumentation, and controls. The electronics response time from fluorescence detection to the activation of triggering signal takes only 360 μs per event. However, controlled delays are programmed as we can only move the droplets at~60 droplets/s through the channel without tearing them apart. Our setup is robust and can work with four wavelengths of fluorescent excitation (i.e., 405, 488, 561, and 637 nm), making it compatible with several fluorescenceemitting biosensors. A summary of fluorochromes compatible with this setup is also provided as a quick selection guide.We present a passive droplet merging device that can achieve efficiencies of 95% and at~200 droplets/s. This merger facilitates the addition of precise amounts of reagents for various purposes,such as starting or terminating a reaction and adding new reagents.

    Finally,we discuss the potential of these three microfluidic technologies in mining the immense pool of metagenomic resources for bioprospecting.

    CRediT author statement

    David Conchouso:Methodology, Investigation, Writing -original draft, Writing - review & editing.Amani Al-Ma’abadi:Resources.Hayedeh Behzad:Resources.Mohammed Alarawi:Resources.Masahito Hosokawa:Methodology.Yohei Nishikawa:Methodology.HarukoTakeyama:Methodology.Katsuhiko Mineta:Methodology, Supervision, Writing - review & editing.Takashi Gojobori:Conceptualization, Methodology, Supervision,Writing - review & editing. All authors have read and approved the final manuscript.

    Competing interests

    The authors have declared no competing interests.

    Acknowledgments

    The work was supported by the grants from King Abdullah University of Science and Technology (KAUST), Saudi Arabia(Grant Nos.BAS/1/1059/01/01,URF/1/1976/03/01,URF/1/1976-17-01, URF/1/1976-20-01, and FCS/1/3326-01-01). Finally, we thank the nanofabrication and the microfluidics Core Laboratories at KAUST for their support in realizing these microfluidic technologies.

    Supplementary material

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

    ORCID

    0000-0002-9788-0977 (David Conchouso)

    0000-0001-8895-5423 (Amani Al-Ma’abadi)

    0000-0001-7216-8337 (Hayedeh Behzad)

    0000-0001-8986-8209 (Mohammed Alarawi)

    0000-0002-7070-4690 (Masahito Hosokawa)

    0000-0001-5443-7230 (Yohei Nishikawa)

    0000-0002-2058-8185 (Haruko Takeyama)

    0000-0002-4727-045X (Katsuhiko Mineta)

    0000-0001-7850-1743 (Takashi Gojobori)

    国产综合懂色| 成人美女网站在线观看视频| 色尼玛亚洲综合影院| 欧美色欧美亚洲另类二区| 18禁在线播放成人免费| 精品一区二区三区视频在线观看免费| 好男人电影高清在线观看| 国产av在哪里看| 精品一区二区三区视频在线| 亚洲aⅴ乱码一区二区在线播放| 国产成人福利小说| 男女床上黄色一级片免费看| 久久这里只有精品中国| 国产蜜桃级精品一区二区三区| 婷婷六月久久综合丁香| 97超视频在线观看视频| 长腿黑丝高跟| 国产一区二区三区在线臀色熟女| 亚洲av免费在线观看| 国产真实伦视频高清在线观看 | 一个人免费在线观看的高清视频| 精品午夜福利在线看| 午夜福利在线在线| 动漫黄色视频在线观看| 欧美另类亚洲清纯唯美| 免费在线观看日本一区| 一个人看视频在线观看www免费| .国产精品久久| 国产精品98久久久久久宅男小说| 91字幕亚洲| 亚洲精品日韩av片在线观看| 国产一区二区在线av高清观看| 深爱激情五月婷婷| 男人的好看免费观看在线视频| 能在线免费观看的黄片| 亚洲avbb在线观看| av专区在线播放| 国产伦在线观看视频一区| 国产精品一区二区三区四区免费观看 | 淫妇啪啪啪对白视频| 男女那种视频在线观看| 久久久精品大字幕| 国产激情偷乱视频一区二区| 国产伦精品一区二区三区四那| 亚洲天堂国产精品一区在线| 麻豆成人av在线观看| 日韩欧美精品v在线| 亚洲欧美日韩高清在线视频| 黄色配什么色好看| 97碰自拍视频| 国产一区二区亚洲精品在线观看| 久久天躁狠狠躁夜夜2o2o| 如何舔出高潮| 中文资源天堂在线| 亚洲真实伦在线观看| 99热只有精品国产| 久久久久九九精品影院| 久久国产精品人妻蜜桃| 亚洲一区二区三区色噜噜| 亚洲av.av天堂| 成人午夜高清在线视频| 亚洲不卡免费看| 特级一级黄色大片| 少妇被粗大猛烈的视频| 国产精品久久视频播放| 国产成人a区在线观看| 亚洲经典国产精华液单 | 一区二区三区免费毛片| 国产精品永久免费网站| 熟女电影av网| 黄片小视频在线播放| 最近最新免费中文字幕在线| 乱码一卡2卡4卡精品| 免费无遮挡裸体视频| 91av网一区二区| 国产精品女同一区二区软件 | 看片在线看免费视频| 不卡一级毛片| АⅤ资源中文在线天堂| 少妇人妻精品综合一区二区 | 中文字幕人成人乱码亚洲影| 日本 欧美在线| 欧美极品一区二区三区四区| 变态另类成人亚洲欧美熟女| 精品乱码久久久久久99久播| 五月伊人婷婷丁香| 亚洲七黄色美女视频| 久久精品国产99精品国产亚洲性色| 国产亚洲精品av在线| 亚洲 国产 在线| 91麻豆精品激情在线观看国产| 欧美性猛交黑人性爽| 精品人妻熟女av久视频| 国产熟女xx| 亚洲精品一区av在线观看| 久久久精品欧美日韩精品| 日韩 亚洲 欧美在线| 国产探花在线观看一区二区| 97碰自拍视频| 观看美女的网站| 国产视频内射| 一本综合久久免费| 成年女人毛片免费观看观看9| 日韩 亚洲 欧美在线| 亚洲av中文字字幕乱码综合| 成人三级黄色视频| 天天躁日日操中文字幕| 国产在线精品亚洲第一网站| 中文字幕av成人在线电影| 亚洲av电影在线进入| 国产美女午夜福利| 老司机深夜福利视频在线观看| 精品人妻视频免费看| 亚洲七黄色美女视频| 99热这里只有精品一区| 三级国产精品欧美在线观看| 我要看日韩黄色一级片| 亚洲avbb在线观看| 国产极品精品免费视频能看的| 国产国拍精品亚洲av在线观看| 亚洲人成电影免费在线| 国产精品99久久久久久久久| 高清日韩中文字幕在线| 18禁黄网站禁片午夜丰满| 日韩欧美在线乱码| 国产伦精品一区二区三区视频9| 中亚洲国语对白在线视频| www.熟女人妻精品国产| 一进一出抽搐gif免费好疼| 成年女人永久免费观看视频| 久久草成人影院| 久久久久国产精品人妻aⅴ院| 亚洲美女视频黄频| 熟女人妻精品中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 色噜噜av男人的天堂激情| 在线观看午夜福利视频| 国产激情偷乱视频一区二区| 男人舔女人下体高潮全视频| 午夜精品在线福利| 欧美绝顶高潮抽搐喷水| 久久中文看片网| 在现免费观看毛片| 一本综合久久免费| 在线免费观看的www视频| 欧美色视频一区免费| 9191精品国产免费久久| 婷婷六月久久综合丁香| 别揉我奶头~嗯~啊~动态视频| 91在线观看av| 久久亚洲精品不卡| 亚洲人成电影免费在线| 看十八女毛片水多多多| 欧美日韩中文字幕国产精品一区二区三区| 免费黄网站久久成人精品 | 岛国在线免费视频观看| 精品一区二区三区av网在线观看| 自拍偷自拍亚洲精品老妇| 91久久精品国产一区二区成人| av黄色大香蕉| 成人国产一区最新在线观看| 久久天躁狠狠躁夜夜2o2o| 天堂网av新在线| 午夜福利在线在线| 日韩欧美精品免费久久 | 一级毛片久久久久久久久女| 中文字幕人妻熟人妻熟丝袜美| 一进一出抽搐gif免费好疼| 深爱激情五月婷婷| 老女人水多毛片| 日韩免费av在线播放| 美女xxoo啪啪120秒动态图 | 亚洲最大成人av| 国产单亲对白刺激| 12—13女人毛片做爰片一| 舔av片在线| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美 国产精品| 99国产综合亚洲精品| 久久久久性生活片| 成人国产综合亚洲| 国产黄a三级三级三级人| 久久99热6这里只有精品| 黄色日韩在线| 丰满乱子伦码专区| 亚洲第一电影网av| av天堂在线播放| 男女视频在线观看网站免费| 精品久久久久久久人妻蜜臀av| 搡老岳熟女国产| 免费观看的影片在线观看| 波多野结衣巨乳人妻| 亚洲国产精品999在线| 18美女黄网站色大片免费观看| 国产欧美日韩精品一区二区| 热99在线观看视频| 51国产日韩欧美| 亚洲欧美日韩东京热| 日韩欧美三级三区| 成人亚洲精品av一区二区| 三级国产精品欧美在线观看| 男女视频在线观看网站免费| 日日摸夜夜添夜夜添av毛片 | 久久久久久久精品吃奶| av专区在线播放| 婷婷色综合大香蕉| 日韩 亚洲 欧美在线| 黄色一级大片看看| 午夜福利免费观看在线| 亚洲人成网站高清观看| 在线观看一区二区三区| 国产视频内射| 两个人的视频大全免费| 老鸭窝网址在线观看| 亚洲av中文字字幕乱码综合| 97超视频在线观看视频| 亚洲avbb在线观看| 人人妻人人澡欧美一区二区| 国产精品嫩草影院av在线观看 | 免费av不卡在线播放| 久久精品综合一区二区三区| 天堂影院成人在线观看| 成人三级黄色视频| 国产中年淑女户外野战色| 久99久视频精品免费| 少妇的逼水好多| av欧美777| 久久久国产成人免费| 特级一级黄色大片| 欧美日韩中文字幕国产精品一区二区三区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本 欧美在线| 成人精品一区二区免费| 男女下面进入的视频免费午夜| 一本一本综合久久| 久久国产精品影院| 日日干狠狠操夜夜爽| 日日夜夜操网爽| 国产一区二区在线观看日韩| 亚洲人成伊人成综合网2020| 999久久久精品免费观看国产| 免费大片18禁| 窝窝影院91人妻| 热99在线观看视频| 成人性生交大片免费视频hd| 小说图片视频综合网站| 亚洲av熟女| 欧美最黄视频在线播放免费| 亚洲综合色惰| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久九九精品二区国产| 欧美激情在线99| 精品久久久久久久人妻蜜臀av| 免费av毛片视频| 国产日本99.免费观看| 麻豆成人av在线观看| 深爱激情五月婷婷| 国产欧美日韩精品一区二区| 男女床上黄色一级片免费看| av黄色大香蕉| av在线观看视频网站免费| 欧美性猛交╳xxx乱大交人| 精品久久久久久久久久久久久| 国产精品亚洲一级av第二区| 99在线人妻在线中文字幕| 成人亚洲精品av一区二区| 成人国产一区最新在线观看| 99在线视频只有这里精品首页| 国产乱人伦免费视频| 最近视频中文字幕2019在线8| 老女人水多毛片| 婷婷色综合大香蕉| 国产成人欧美在线观看| 国产三级中文精品| 99在线视频只有这里精品首页| 国产精品98久久久久久宅男小说| 女人十人毛片免费观看3o分钟| 天堂av国产一区二区熟女人妻| 国产蜜桃级精品一区二区三区| 美女高潮的动态| 色哟哟哟哟哟哟| 最近最新免费中文字幕在线| 欧美不卡视频在线免费观看| 国产高清有码在线观看视频| 国产午夜精品论理片| 亚洲熟妇中文字幕五十中出| 久久久久久久亚洲中文字幕 | 亚洲在线自拍视频| 免费av毛片视频| 精品午夜福利视频在线观看一区| 最好的美女福利视频网| 在现免费观看毛片| 国产一级毛片七仙女欲春2| 欧美在线一区亚洲| 一个人看视频在线观看www免费| a级毛片免费高清观看在线播放| 国产在线精品亚洲第一网站| 黄片小视频在线播放| 一卡2卡三卡四卡精品乱码亚洲| av在线观看视频网站免费| 国产真实乱freesex| 麻豆成人av在线观看| av视频在线观看入口| 看黄色毛片网站| 亚洲av成人av| 啦啦啦观看免费观看视频高清| 国产主播在线观看一区二区| 亚洲久久久久久中文字幕| av在线观看视频网站免费| 别揉我奶头~嗯~啊~动态视频| 国产亚洲精品久久久久久毛片| 亚洲七黄色美女视频| 在线观看舔阴道视频| 女生性感内裤真人,穿戴方法视频| 九色国产91popny在线| 色噜噜av男人的天堂激情| 首页视频小说图片口味搜索| 亚洲自拍偷在线| 日韩 亚洲 欧美在线| 国产精品一区二区性色av| 性色avwww在线观看| 亚洲色图av天堂| 日日干狠狠操夜夜爽| 一本综合久久免费| 我要搜黄色片| 内射极品少妇av片p| 丰满乱子伦码专区| 在线国产一区二区在线| 欧美日韩中文字幕国产精品一区二区三区| 午夜视频国产福利| 琪琪午夜伦伦电影理论片6080| 又黄又爽又刺激的免费视频.| 最新中文字幕久久久久| 精品日产1卡2卡| 精品久久久久久,| 国产日本99.免费观看| 大型黄色视频在线免费观看| 51国产日韩欧美| 91久久精品电影网| 国产成人aa在线观看| 国产免费一级a男人的天堂| 禁无遮挡网站| 成人特级黄色片久久久久久久| 久久天躁狠狠躁夜夜2o2o| 少妇高潮的动态图| 91久久精品国产一区二区成人| 欧美精品啪啪一区二区三区| 久久久久久久精品吃奶| 国产白丝娇喘喷水9色精品| 一夜夜www| 久久久久免费精品人妻一区二区| 国产精品自产拍在线观看55亚洲| 免费av观看视频| 亚洲av日韩精品久久久久久密| 欧美潮喷喷水| 久久久久久大精品| 日韩欧美在线乱码| 国内毛片毛片毛片毛片毛片| 日韩国内少妇激情av| 欧美性猛交╳xxx乱大交人| 日韩欧美在线二视频| 精品人妻偷拍中文字幕| 午夜福利免费观看在线| 免费高清视频大片| 免费在线观看成人毛片| 成年女人看的毛片在线观看| 亚洲欧美日韩卡通动漫| 小蜜桃在线观看免费完整版高清| 久久久久久久久大av| 精品午夜福利视频在线观看一区| 嫩草影视91久久| 99riav亚洲国产免费| 色5月婷婷丁香| 有码 亚洲区| 又爽又黄a免费视频| 深夜a级毛片| 三级毛片av免费| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品影视一区二区三区av| 亚洲va日本ⅴa欧美va伊人久久| 乱人视频在线观看| 国产色婷婷99| 久久国产乱子免费精品| 99精品久久久久人妻精品| 亚洲中文日韩欧美视频| 成人国产一区最新在线观看| 国产精品久久视频播放| 欧美日本视频| 国产精品久久久久久精品电影| 欧美色欧美亚洲另类二区| 国产午夜福利久久久久久| 欧美激情在线99| 丁香六月欧美| 精品一区二区免费观看| 亚洲一区高清亚洲精品| 午夜福利在线观看免费完整高清在 | 熟女电影av网| 精品免费久久久久久久清纯| 日本在线视频免费播放| 亚洲精品在线观看二区| 18美女黄网站色大片免费观看| 日本一本二区三区精品| av在线蜜桃| 欧美一区二区亚洲| 久久久久性生活片| 亚洲欧美日韩高清在线视频| 国产精品综合久久久久久久免费| 久久久久久国产a免费观看| 国产色婷婷99| 亚洲av电影在线进入| 久久久久免费精品人妻一区二区| 岛国在线免费视频观看| 小说图片视频综合网站| 精品一区二区三区视频在线| 91午夜精品亚洲一区二区三区 | 少妇的逼水好多| 亚洲性夜色夜夜综合| 亚洲一区二区三区不卡视频| 欧美在线一区亚洲| 夜夜看夜夜爽夜夜摸| 毛片女人毛片| 嫩草影视91久久| 亚洲成a人片在线一区二区| 国产野战对白在线观看| 免费黄网站久久成人精品 | 波多野结衣巨乳人妻| 乱人视频在线观看| 激情在线观看视频在线高清| 桃红色精品国产亚洲av| 两个人的视频大全免费| 欧美色欧美亚洲另类二区| 少妇人妻一区二区三区视频| 久久伊人香网站| 91麻豆精品激情在线观看国产| 淫妇啪啪啪对白视频| 亚洲国产欧美人成| 欧美绝顶高潮抽搐喷水| 18美女黄网站色大片免费观看| 动漫黄色视频在线观看| 国产淫片久久久久久久久 | 给我免费播放毛片高清在线观看| 亚洲成人久久性| 能在线免费观看的黄片| 午夜两性在线视频| 免费搜索国产男女视频| 看黄色毛片网站| 婷婷六月久久综合丁香| 日本撒尿小便嘘嘘汇集6| 国产精品亚洲一级av第二区| 深爱激情五月婷婷| 国产精品99久久久久久久久| 又黄又爽又免费观看的视频| 最后的刺客免费高清国语| 国产免费一级a男人的天堂| 一级黄片播放器| av国产免费在线观看| 国产成人啪精品午夜网站| 国产精品免费一区二区三区在线| 国产激情偷乱视频一区二区| 亚洲av第一区精品v没综合| 成人国产综合亚洲| 亚洲最大成人手机在线| 精品一区二区三区人妻视频| 久久午夜福利片| av天堂中文字幕网| 啪啪无遮挡十八禁网站| 亚洲最大成人中文| 国产精品爽爽va在线观看网站| 精品一区二区三区av网在线观看| 给我免费播放毛片高清在线观看| 一卡2卡三卡四卡精品乱码亚洲| 香蕉av资源在线| 免费观看精品视频网站| 免费大片18禁| 欧美区成人在线视频| 午夜精品久久久久久毛片777| 人人妻人人澡欧美一区二区| 国产精品一及| 国产激情偷乱视频一区二区| 亚洲七黄色美女视频| 在线观看免费视频日本深夜| 成人性生交大片免费视频hd| 18禁黄网站禁片免费观看直播| 日韩欧美国产在线观看| 90打野战视频偷拍视频| 免费无遮挡裸体视频| 毛片女人毛片| 村上凉子中文字幕在线| 中文字幕人妻熟人妻熟丝袜美| АⅤ资源中文在线天堂| 内射极品少妇av片p| 最后的刺客免费高清国语| 在线观看av片永久免费下载| 老司机福利观看| 美女高潮喷水抽搐中文字幕| 午夜日韩欧美国产| 99国产精品一区二区蜜桃av| av在线蜜桃| 亚洲一区高清亚洲精品| 国产91精品成人一区二区三区| 成人一区二区视频在线观看| 69人妻影院| 国产老妇女一区| 狂野欧美白嫩少妇大欣赏| 国产不卡一卡二| 夜夜躁狠狠躁天天躁| or卡值多少钱| 在线观看66精品国产| 久久久国产成人免费| 亚洲欧美日韩无卡精品| 国产精品一区二区免费欧美| 深爱激情五月婷婷| 熟女电影av网| 两人在一起打扑克的视频| 丰满的人妻完整版| 午夜日韩欧美国产| 中亚洲国语对白在线视频| 国产爱豆传媒在线观看| 99久国产av精品| 宅男免费午夜| 久久久久免费精品人妻一区二区| www.www免费av| 久久热精品热| 内地一区二区视频在线| 亚洲国产欧美人成| 国产白丝娇喘喷水9色精品| 十八禁国产超污无遮挡网站| 免费大片18禁| 亚洲第一区二区三区不卡| x7x7x7水蜜桃| 久久婷婷人人爽人人干人人爱| 国产探花在线观看一区二区| 成人午夜高清在线视频| 91久久精品电影网| 亚洲人成网站在线播| 波多野结衣高清作品| 两性午夜刺激爽爽歪歪视频在线观看| 97超视频在线观看视频| 国产精品亚洲美女久久久| 一本久久中文字幕| 村上凉子中文字幕在线| 久久久久久久亚洲中文字幕 | 亚洲,欧美,日韩| 99riav亚洲国产免费| 免费看a级黄色片| 国产高清视频在线观看网站| 国产亚洲精品av在线| 99riav亚洲国产免费| 亚洲一区二区三区色噜噜| 别揉我奶头~嗯~啊~动态视频| 亚洲av成人不卡在线观看播放网| 免费在线观看亚洲国产| 深夜a级毛片| 99热这里只有是精品在线观看 | 国产淫片久久久久久久久 | 嫩草影院精品99| 久久久久久久久久成人| 国产黄色小视频在线观看| 激情在线观看视频在线高清| 日本 欧美在线| 国产精品久久久久久久电影| 亚洲欧美日韩高清专用| 成人毛片a级毛片在线播放| 精品人妻熟女av久视频| 亚洲真实伦在线观看| 精品午夜福利视频在线观看一区| 国产主播在线观看一区二区| 国产大屁股一区二区在线视频| 久久精品影院6| 午夜日韩欧美国产| а√天堂www在线а√下载| 婷婷精品国产亚洲av| 亚洲人成伊人成综合网2020| 丰满人妻一区二区三区视频av| 18+在线观看网站| 嫁个100分男人电影在线观看| 免费看a级黄色片| 亚洲人成网站高清观看| 欧美成人性av电影在线观看| 日韩人妻高清精品专区| 精品久久久久久久久久久久久| 俄罗斯特黄特色一大片| 日本黄大片高清| 极品教师在线视频| 久久久久久久午夜电影| 在线观看免费视频日本深夜| 国产单亲对白刺激| 有码 亚洲区| 在线a可以看的网站| 国产精品亚洲av一区麻豆| 99在线人妻在线中文字幕| 婷婷精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| 高清日韩中文字幕在线| 老司机深夜福利视频在线观看| 99精品久久久久人妻精品| 久久精品久久久久久噜噜老黄 | 国产一级毛片七仙女欲春2| 国内精品久久久久久久电影| 午夜福利欧美成人| 久久久国产成人精品二区| 免费无遮挡裸体视频| 久久性视频一级片| 精品人妻1区二区| 欧美另类亚洲清纯唯美| 国产精品久久久久久精品电影| 国产精品一区二区免费欧美| 久久精品夜夜夜夜夜久久蜜豆| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲人成网站在线播放欧美日韩| 天天躁日日操中文字幕| 一区二区三区高清视频在线| 精品人妻熟女av久视频| avwww免费| 夜夜躁狠狠躁天天躁|