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

    Recent advances in droplet microfluidics for microbiology

    2022-06-20 06:18:48ZiyiHeHoWuXinghuYnWuLiu
    Chinese Chemical Letters 2022年4期

    Ziyi He,Ho Wu,b,Xinghu Yn,Wu Liu

    a State Key Laboratory of Agricultural Microbiology,College of Veterinary Medicine,Huazhong Agricultural University,Wuhan,430070,China

    b Interdisciplinary Sciences Institute,Huazhong Agricultural University,Wuhan 430070,China

    c State Key Laboratory of Agricultural Microbiology,College of Animal Sciences and Technology,Huazhong Agricultural University,Wuhan 430070,China

    d Hubei Hongshan Laboratory,Wuhan 430070,China

    e The Cooperative Innovation Center of Sustainable Pig Production,Wuhan 430070,China

    f Hubei Provincial Engineering Laboratory for Pig Precision Feeding and Feed Safety Technology,Wuhan 430070,China

    g Institute of Pharmaceutical Analysis,School of Pharmaceutical Sciences,Cheeloo College of Medicine,Shandong University,Jinan 250012,China

    ABSTRACT Advances in microbiology rely on innovations in technology.Droplet microfluidics,as a versatile and powerful technique that allows high-throughput generation and manipulation of subnanoliter volume droplets,has become an indispensable tool shifting experimental paradigms in microbiology.Droplet microfluidics has opened new avenues to various microbiological research,from resolving single-cell heterogeneity to investigating spatiotemporal dynamics of microbial communities,from precise quantitation of microbiota to systematic decipherment of microbial interactions,and from isolating rare and uncultured microbes to improving genetic engineered strains.In this review,we present recent advances of droplet microfluidics in various fields of microbiology:i)microbial cultivation,ii)microorganism detection and characterization,iii)antibiotic susceptibility testing,iv)microbial interactions,v)microbial biotechnology.We also provide our perspectives on the challenges and future directions for droplet microfluidic-based microbiology research.

    Keywords:Droplet microfluidics Microorganisms Antibiotic susceptibility testing Microbial interactions Microbial biotechnology

    1.Introduction

    Microorganisms,including bacteria,viruses,fungi and algae,live in almost every habitat on the planet,from deep-sea vents to the highest clouds to the digestive tracts of humans[1].They are so small,diverse,highly adaptable and highly mutable,that advanced technologies are essential to identify,characterize,isolate,culture,engineer and utilize microorganisms,and to uncover complex interactions in microbial ecosystems and/or between microorganisms and their hosts[2,3].Droplet microfluidics,as one subcategory of microfluidics that generates and manipulates subnanoliter volume droplets inside microdevices,has become an indispensable tool changing the paradigm of experimentation in microbiology[4-6].Droplet microfluidics provides multiple unique benefits for microbiological research.Firstly,by encapsulating microorganisms in microdroplets,a small quantity of cells,or even single cells can be confined in isolated environments[7];interspecies competition and biases due to growth rate differences are eliminated,which facilitate recovery of rare or slow-growing microorganisms from complex ecosystems[8];fast accumulation of metabolites in droplets enables easy activation of concentration-dependent processes such as quorum sensing[9].Secondly,droplet microfluidics allows rapid droplet production with a speed up to 20,000 per second and massive droplet analysis,opening up the possibility of ultrahigh-throughput identification,screening and strain improvement of microorganisms[5,10-12].Thirdly,in droplet microfluidics,configurable channel designs and incorporated control modules engender precise droplet manipulations,including fusion,mixing,splitting,long-term incubation and sorting[13];these allow,for example,introducing multiple detection reagents or stimuli to microorganisms[14,15],creating well-controlled and variable culture environments[16],or selectively collecting improved microbial strains[17].

    Fig.1.Applications of droplet microfluidics in various fields of microbiology.

    The above features collectively make droplet microfluidics a robust and valuable platform for microbiology studies.Combining with advanced sequencing technologies[18]or dynamic fluorescence imaging[19],droplet microfluidics helps bring microbiologists closer to the ultimate goal of microbial ecology:observing and understanding who does what,when,where and why[3].In the present review,we focus on the advances in droplet microfluidics for microbiology within the past five years.We start the review with a brief introduction of droplet generation strategies in microfluidics,followed by introducing the advances and applications of droplet microfluidics in various fields of microbiology:i)microbial cultivation,ii)microorganism detection and characterization,iii)antibiotic susceptibility testing(AST),iv)microbial interactions,v)microbial biotechnology(Fig.1).Finally,we will discuss the challenges and future perspectives of droplet microfluidics for microbiology research,inspiring the collaborations between microbiologists and microfluidic researchers to push the boundaries of this field beyond what we now consider possible.

    2.Droplet generation in microfluidics

    Most droplet generation approaches in microfluidics could be categorized into flow-based strategy and array-based strategy.In flow-based strategy,droplets are formed through the extrusion and shearing of two mutually immiscible phases in a microchannel[4].According to different channel configurations,the flow-based strategy can be classified as co-flowing,T-junction,flow-focusing and step emulsification methods[20,21].Flow-focusing is the most common geometry,where the dispersed phase(aqueous phase)and continuous phases(oil phase)flow co-axially through a contraction region and generate an elongation filament that eventually breaks into droplets[7].Droplet volumes are controlled by the flow rate ratios and channel dimensions[22].To realize on-demand droplet generation,active modules,such as electrical,magnetic,thermal,or mechanical controls,can be incorporated into flow-based droplet microfluidics,offering additional handles to regulate droplet volumes and generation rates[23].In array-based strategy,aqueous droplets are spatially organized in arrays of microwells,microchambers or hydrophilic micropatterns,which are isolated from each other by the oil phase[24,25].Droplet volumes are mainly controlled by the dimensions of microwells,microchambers or micropatterns.Array-based strategy allows droplets to be indexed by their spatial layout,which will be beneficial to identify,monitor,and address individual droplets[26,27].Both droplet generation strategies have been successfully utilized in microbiology research.By using a microorganism suspension as the aqueous phase,microorganisms can be encapsulated in droplets,and the distribution of the cell number per droplet follows the Poisson distribution[28].

    3.Microbial cultivation

    Microbial cultivation is the foundational method in microbiology,to grow and diagnose microorganisms,to isolate and obtain pure cultures,and to investigate effects of environmental factors[29].In traditional culture methods,such as agar plating and nutrient broth,it takes a long time for microbes to grow and multiply to be detectable,and rare or slow-growing species in the microbial community can be omitted because of interspecies competitions and growth rate differences[30].Droplet-based microbial cultivation can resolve these issues by isolating and culturing a small quantity of cells or single cells in microdroplets,which enables real-time monitoring of microbial growth by microscopy,prevents interspecies competition,and allows precise control and flexible regulation of microbial growth conditions.

    Droplet microfluidics represents one of the most powerful techniques to isolate and cultivate single microbial cells in confined environments[31].The growth of microbes within individual droplets can be tracked by recording cell numbers under bright field or performing fluorescence-based quantification over time[26,32].Dynamic microbial monitoring is valuable for investigating the effects of different nutrient sources as well as stress conditions,such as acid stress[32,33],heavy metal exposition[34,35],antibiotics[36]or electric field[37],on microorganisms.Juskováet al.presented a microfluidic system named“Basicles”,which generated giant unilamellar lipid vesicles(GUVs)for single-bacteria encapsulation and immobilized the GUVs by hydrodynamic traps for real-time monitoring of microbial growth and production[36].Escherichia coli(E.coli)could be cultivated in GUVs for up to 15 h,and the growth of bacteria per GUV was monitored based on the increase of sfGFP(super-folder green fluorescent protein)fluorescence.A significant advantage of this system was that membranepermeable reagents could passively diffuse across the liposomal membrane[38],achieving controlled feeding and stimulation of the encapsulated bacteria.Therefore,riboflavin biosynthesis behaviors ofE.coliwere assessed under different conditions with inducers or antibiotics.Compared to optical-based growth assays,molecular techniques,such as quantitative PCR(qPCR)and gene sequencing,can enable more comprehensive and systematic growth measurements of hundreds of gut bacterial species in droplets.Villaet al.developed a MicDrop platform to isolate and culture individual gut bacteria in droplets from human gut microbiota(Fig.2a)[39].Since a large proportion of gut microbiome are obligate anaerobes,microbial encapsulation and incubation in droplets were performed in an anaerobic chamber to ensure high cell viability.To track microbial growth,bacteria were recovered from droplets at different culture time points,followed by DNA sequencing and total quantification of the 16S rRNA gene,providing absolute levels of each sequence variant across all droplets.MicDrop enabled to cultivate 2.8 times more bacterial taxa than standard bulk culture.The high-throughput and quantitative growth measurements allowed characterizing carbohydrate utilizations of gut bacterial taxa from distinct human stool samples.

    Droplet-based microbial cultivation facilitates isolation of rare or slow-growing microorganisms from complex community.The Du group described the microfluidic streak plate(MSP)for highthroughput single-microbe separation and cultivation in nanoliter droplets(Fig.2b)[40].Droplets generated by a microfluidic device were streaked onto petri dishes prefilled with carrier oil,forming sessile-droplet arrays.Microbes could be cultured for up to five months at room temperature.Rapid recovery of microbes after growth could be easily achieved by using a toothpick or an automatic droplet picker[41],allowing scale-up cultivation and downstream characterization.Results from 16S rRNA sequencing demonstrated that MSP method enabled detection of the richest species diversity with better coverage of rare species compared to the agar plate method.The MSP platform provided a powerful tool in isolating the rare species from complex environmental samples,including soil[42],termite gut microbiome[43]and deep-sea sediments[41].Wattersonet al.utilized droplet microfluidics in an anaerobic chamber to isolate and cultivate gut microbiota[44].A unique feature of this system was the incorporation of a droplet sorting unit,which allowed isolating slow growing organisms based on the optical density(OD)of droplets and further amplified the abundance of rare taxa.This system offered an increase of taxonomic richness by 15% to 410% over conventional plates depending on the culture medium,and the droplet sorting enriched six separate Alistipes populations,which were either not detected or in very low abundance in plate-based cultivations.

    Fig.2.Microbial cultivation in droplets.(a)MicDrop platform for isolating and culturing gut bacteria in droplets;distinct colony morphologies could be observed across droplets of gut microbial community;scale bars:(i)100 μm,(ii)50 μm.Reproduced with permission[39].Creative Commons Attribution 4.0 International License.(b)Microfluidic streak plate(MSP)for single-microbe separation and cultivation in nanoliter droplets;scale bars:(i)1 cm,(ii)1 mm.Reproduced with permission[40].Copyright 2016,American Society for Microbiology.(c)Comparisons of bacterial biofilm formations in microfluidic flow cell and dynamic sessile-droplet habitats(DSHs);scale bars:100 μm.Reproduced with permission[45].Copyright 2018,Wiley-VCH.

    Apart from single microorganism culture,droplet microfluidics also enables controllable cultivation of bacterial biofilms.The Du group introduced a dynamic sessile-droplet habitat(DSH)platform,containing an array of hydrophilic patterns in different sizes at the downstream of a T-junction droplet generator(Fig.2c)[45].By flowing nanoliter plugs containing bacteria through the patterns,an array of discrete sessile droplets with inoculated bacteria was produced,which could form spatially confined biofilms and culture over a long term by flowing media plugs through the patterns subsequently.The DSH method could avoid channel clogging and flow disruption that always appeared in microfluidic flow cellbased method,providing more consistent hydrodynamic environments for biofilm growth and stimulation.This platform allowed lifespan monitoring ofPseudomonas aeruginosabiofilm formation by fluorescence,and the roles of flagella and pili in biofilm development were investigated.

    4.Microorganism detection and characterization

    Detection and characterization of microorganisms are pivotal from fundamental microbiology to infectious disease diagnosis to food industry[46,47].Attributed to the capabilities of rapid droplet production and massive droplet analysis,droplet microfluidics paves the way to ultrahigh-throughput and highly parallelized microbial analysis.Recently,various detection techniques have been utilized to analyze microorganisms in droplets,including optical-,electrical-,mass spectrometry-,and molecular biology-based methods[10,48-51].In this section,we will discuss the recent advances in droplet microfluidic-based microorganism detection and characterization,mainly focusing on the optical detection and molecular biology methods including droplet digital PCR,sequencing-based characterization and Cas13a-based detection.

    4.1.Optical detection

    Optical detection of droplets is one of the most frequently used techniques[52],which is straightforward and enables real-time analysis.Hsiehet al.developed resazurin-amplified picoarray detection(RAPiD)for simple and precise counting of viable bacteria[53].Using vacuum-assisted sample loading and oil-driven digitization,single bacteria were enclosed in picochambers with resazurin,an indicator dye that could be reduced by viable cells to strongly fluorescent resorufin.By counting the number of chambers exhibiting fluorescence signals,viable bacteria inE.coliandStaphylococcus aureus(S.aureus)samples could be quantified with high precision.Viable bacteria enumeration could also be achieved by using microdroplet turbidity imaging-based digital counting,which was label-free and allowed signal acquisition using a smartphone camera,benefiting point-of-care applications[54].

    Despite being highly convenient for microbial analysis,fluorescent detection relies on the labeling of fluorescent tags or introduction of fluorescent substrates,which is only applicable to some specific target molecules and may affect microbial metabolism.In contrast,single-cell Raman spectrum(SCRS),as a label-free and non-invasive technique for metabolite detection in a given cell based on the vibrational signals from molecules,can provide intrinsic biochemical“fingerprints”for individual cells[55,56].The Xu group coupled the SCRS with droplet microfluidics for culture-free and label-free microbial characterization and screening.They developed Raman-activated gravity-driven single-cell encapsulation and sequencing(RAGE-Seq)for phenome–genome profiling of single bacterial cells[57].Individual cells were phenotypically screenedviaSCRS at first,and the targeted cells were respectively encapsulated in picoliter droplets,which were then exported in a precisely indexed,“one-cell-one-tube”manner.Genomic DNA from the sorted cell was individually amplified in droplet and subjected to whole genome sequencing.The RAGEseq system allowed to link metabolic phenotype with high-quality genotype at the one-cell resolution.However,the throughput of the RAGE-seq system was 2 cells per minute,which was relatively low.To achieve higher throughput,they further presented positive dielectrophoresis-based Raman-activated droplet sorting(pDEP-RADS),in which fast-moving single microbes were effi-ciently trapped by pDEP field,profiledviaSCRS,encapsulated in droplets for DEP-based sorting,and collected for sequencing-based identification(Fig.3a)[58].The pDEP-RADS was capable of processing hundreds of cells per minute,allowing culture-free and label-free screening of enzyme functionsin vivo.Two previously unknown diacylglycerol acyltransferases(DGATs)variants were discovered from an engineered yeast library,which were undetectableviaconventional fluorescence-activated cell sorting.

    4.2.Droplet digital PCR

    Droplet digital PCR(ddPCR)is the first commercially significant application of droplet microfluidics[4].As a powerful technology allowing absolute nucleic acid quantitation,ddPCR has been widely adopted in the detection of pathogenic microorganisms[59].ddPCR involves partitioning PCR solution into tens of thousands of nanoliter droplets,with the requirement that each droplet contains either“one”or“zero”target nucleic acid(Fig.3b)[60].Droplets are then amplified individually,and the concentration of target nucleic acid can be obtained by counting the number of positive(fluorescent)droplets at the end point and calculating based on Poisson statistics[61].ddPCR is independent of external standard calibrations and exhibits higher sensitivity and precision compared to quantitative PCR(qPCR).Recently,it has been successfully applied to quantify various bacterial species in diverse samples[62-65],such asLactococcus lactisgroup during milk fermentation[66],Listeria monocytogenesrecovered from biofilms[67],total coliforms andE.coliin urban water[68],andHelicobacter pyloriin stomach and stool samples[69].Besides specific bacteria detection,ddPCR could be combined with gene sequencing to provide a quantitative framework measuring absolute abundances of individual bacterial taxa in complex communities[70-72].ddPCR was utilized to measure the total amount of 16S rRNA gene copies,which represented the total microbial load in the sample,and 16S rRNA gene sequencing offered the relative abundance of each taxon;the absolute abundances of taxa were calculated by multiplying the total load with the relative abundance.Absolute instead of relative microbial abundances could provide more accurate analyses of changes in taxa among different samples.This framework has been applied to determine the differential effects of diet on each taxon in murine lumenal and mucosal samples along the gastrointestinal(GI)tract[71].ddPCR is also advantageous for early diagnosis of viral infection[73-77].For example,Donget al.reported the rapid detection of fowl adenovirus serotypes 4/10(FAdV-4/10)[78].Compared to conventional PCR detection,ddPCR exhibited 1000-fold higher sensitivity.This assay may serve as a tool for detection of FAdV-4 contamination in live-attenuated vaccines.Jianget al.introduced the design and implementation of an automatic,integrated instrument with a built-in power supply for continuousflow ddPCR,as demonstrated by absolute quantitation of hepatitis B virus(HBV)in serum samples[79].The Shen group developed a portable self-partitioning SlipChip(sp-SlipChip)microfluidic device for quantification of BK virus(BKV)DNA[80].A major feature of this method is a“chain-of-pearls”microchannel,with which individual droplets for PCR reactions are generated in a robust and reproducible manner through a slip-induced self-partitioning mechanism(Fig.3c).As a result,this system allows for“sample-todigital-result”detection,making it ideal for point-of-care analysis.Earlier,the group also utilized this droplet formation strategy for digital loop-mediated isothermal amplification(LAMP)quantification of human papillomavirus(HPV)DNA[81].

    In response to the outbreak of SARS-CoV-2,real-time reversetranscriptase PCR(RT-PCR)has been considered as the gold standard for clinical diagnosis.However,because of the low viral load collected by throat swabs,a considerable false-negative rate has been recognized as a key challenge to minimize the risk of viral transmission based on a false sense of security[82].Thus,Suoet al.proposed an optimized ddPCR assay as a complementary to reduce the false negatives[83].Based on the results with clinical samples from 77 patients,ddPCR exhibited not only a significantly enhancement in sensitivity,but also a much lower negative likelihood ratio.The overall accuracy was improved from 47% for RT-PCR to 95% for ddPCR.Similarly,Yuet al.parallelly applied RT-PCR and ddPCR tests to 323 samples from 76 SARS-CoV-2-infected patients[84].ddPCR was found more accurate for low-viral-load samples.Besides,several new droplet manipulation strategies have been utilized for point-of-care diagnosis of SARS-CoV-2,including the commercialized thin-film-transistor digital-microfluidics[85].Recently,to avoid the use of bulky lab-bound equipment and thus accelerate the workflow for SARS-CoV-2 detection,Chenet al.used a head-flattened pipette tip with an elliptical cross-section to manually produce monodispersed droplets ~150 μm to 350 μm in size for ddPCR[86].Moreover,using a smartphone as the detector,a limit of detection(LOD)of 3.8 copies per 20 μL could be achieved,comparable to commercial PCR platforms.

    Fig.3.Microorganism detection and characterization in droplets.(a)Positive dielectrophoresis–based Raman-activated droplet sorting;scale bar:200 μm.Reproduced,with permission[58].Creative Commons Attribution 4.0 International License.(b)Schematic diagram of droplet digital PCR.Reproduced,with permission[60].Creative Commons Attribution 4.0 International License.(c)Slip-induced self-partitioning for droplet generation and ddPCR quantification of BK virus DNA in different concentrations;scale bar:2 mm.Reproduced with permission[80].Copyright 2021,Elsevier.

    4.3.Sequencing-based characterization

    In recent years,the revolutionary improvements of sequencing technology allow high-throughput microbial identification and genotyping from a large,mixed communities without prior cultivation[46,87,88].Conventional sequencing techniques,such as metagenomic sequencing and 16S rRNA sequencing,require homogenization of input material,which will result in the loss of the heterogeneity at the single-cell level as well as the spatial microbiome information[89].Combining droplet microfluidics with sequencing techniques can circumvent the above limitations,by confining and indexing individual cells or spatially organized cell clusters in thousands of independent microdroplets[90,91].Lanet al.described single-cell genome sequencing(SiC-seq)using droplet microfluidics(Fig.4)[14].The principal strategy of SiC-seq is to isolate single cells,extract and purified genomic DNA from the cellular matter,and label all DNA sequences from the same genome with a unique barcode sequence for single-cell indexing.The whole workflow was assisted by droplet microfluidics as it is capable of processing millions of compartmentalized picoliter reactions and introducing reagentsviadroplet merging.Single cells were first encapsulated in agarose microgels,which enabled confined genomic extraction and purification from individual cells through a series of bulk enzymatic and detergent lysis steps.By re-encapsulating the microgels into droplets containing tagmentation reagents and introducing PCR materials and barcode reagents sequentially,genomic DNAs were fragmented and attached with barcodes,which were ready for subsequent sequencing.The sequencing data were grouped according to the barcodes,providing a library of singlecell genomes for downstream analysis.The SiC-seq demonstrated an ultra-high-throughput of>50,000 cells per run in a few hours,which was highly promising to uncover the genetic heterogeneity in diverse cell populations.Using SiC-seq,the distributions of antibiotic-resistance genes,virulence factors,and transduction potential in microbial communities from environmental samples were characterized.Droplet-based barcoding strategy has also been utilized in single-cell RNA sequencing(scRNA-seq),which has recently become a gold standard method for molecularly defining cell states and phenotypes[92].However,conventional scRNA-seq is optimized for mammalian cells,which is ineffective when directly applied to microbial cells that contain less mRNA and have cell walls.To resolve this issue,Liuet al.developed isogenic colony sequencing(ICO-seq),in which single yeast cells were encapsulated in microgels and cultured to yield isogenic yeast colonies of tens to hundreds of cells before subjected to RNA sequencing[93].The culture process amplified the amount of RNA available for deep sequencing and reduced errors.Cell walls were removedviabulk enzymatic treatment of the microgels,making the mRNA more accessible.Released mRNA was captured on barcoded beads,and followed by reverse transcription,amplification and sequencing,transcriptome-wide screening of thousands of isogenic yeast colonies could be achieved in a single experiment.

    To connect microbial function to phylogeny in single cells from a community,the epicPCR(Emulsion,Paired Isolation and Concatenation PCR)method was developed[94].Single cells enclosed in polyacrylamide microgels were lysed and re-encapsulated in droplets with fused primers.The key step of the epicPCR is the in-droplet fusion of functional target genes with 16S rRNA gene(Fig.5a),which presents in almost all bacteria and contains species-specific signature sequences useful for bacteria identification.Sequencing the fused amplicons provided each read containing 16S rRNA and functional genes simultaneously,allowing phylogenetic identification of the functional taxa.The epicPCR provides a versatile technique to disclose the linked phylogenetic and functional information from millions of uncultured single cells in parallel.It has been applied to discover new putative sulfate reducers in a microbial community from a freshwater lake[94]and capture virus–host interactions[95].To characterize spatial organization of a microbiome at micrometer-scale resolution,Shethet al.presented the metagenomic plot sampling by sequencing(MaPSseq)[96].In MaPS-seq,intact microbiome samples were first immobilized in a gel matrix and cryofractured into particles.After cell lysis,the particles contained genomic DNA immobilized in the original arrangement,preserving local spatial information.Similar to the aforementioned methods,individual particles were coencapsulated with barcoded gel beads in droplets for 16S rRNA amplification,and microbiome in the corresponding particles could be obtained by deep sequencing.MaPS-seq is a promising technique for studying microbial biogeography in complex habitats,and it has been used to survey microbiome in three regions of the mouse GI tract,revealing heterogeneous microbial distributions with positive and negative co-associations between specific taxa(Fig.5b).

    Fig.4.Single-cell genome sequencing(SiC-seq)using droplet microfluidics;scale bar:100 μm.Reproduced with permission[14].Copyright 2017,Springer Nature.

    Apart from droplet-based sequencing,droplet microfluidics can also serve as a robust platform for whole-genome amplification[97].Multiple displacement amplification(MDA)is a widely used technique for amplification of minute amounts of DNA samples to a reasonable quantity before whole genome sequencing[98,99].Random bias of amplification could be reduced through compartmentation[100].Shiet al.integrated a microbe enrichment microfluidic device with droplet-based MDA for the preparation of airway microbiome DNA from sputa for shotgun metagenomic sequencing[101].A microfluidic device with four cascades of 20 selection units was developed to deplete the host cells and enrich particles with diameters less than 5 μm,which were suitable for most bacteria and virus.Genomic DNA was extracted and then amplified with MDA in uniformed droplets,following by the shotgun metagenome sequencing.This method presented greater microbiome complexity and more genome information at a strain level comparing with direct sequencing,with higher sensitivity for identifying prophages and DNA viruses.

    Fig.5.Microbial characterization in droplets based on sequencing technique.(a)epicPCR for phylogenetic identification of the functional taxa.Reproduced with permission[94].Creative Commons Attribution 4.0 International License.(b)Metagenomic plot sampling by sequencing(MaPS-seq)for characterizing spatial organization of a microbiome at micrometer-scale resolution.Reproduced with permission[96].Copyright 2019,Springer Nature.

    Fig.6.CARMEN platform combined with Cas13a-based detection for highly scalable,multiplexed detection of pathogens;scale bar = 200 μm.Reproduced with permission[104].Creative Commons Attribution 4.0 International License.

    4.4.Cas13a-based detection

    Rapid,sensitive,high-throughput and point-of-care detection of pathogenic microorganisms is critical for disease diagnosis and surveillance.CRISPR-Cas13a is a recently discovered RNA-guided RNase,which can be programmed with a guide CRISPR RNA(cr-RNA)to target a single strand RNA of interest[102].Once binding to its RNA target,activated Cas13a exhibits nonspecific endonuclease activity that can cleave nearby nontargeted RNAs.This mechanism motivated the development of SHERLOCK(specific high-sensitivity enzymatic reporter unLOCKing)technique,which enables rapid detection of nucleic acids with attomolar sensitivity based on nucleic acid amplification and Cas13a-mediated collateral cleavage of a reporter RNA[103].Combining Cas13abased detection with droplet microfluidics opens a new possibility for pathogen detection.An outstanding work was presented by the Blainey group,which developed the CARMEN(combinatorial arrayed reactions for multiplexed evaluation of nucleic acids)platform and combined it with SHERLOCK for highly scalable,multiplexed detection of pathogens(Fig.6)[104].First,Cas13adetection mixes(containing different specific crRNA)or amplified samples were emulsified by droplet microfluidics to form nanoliter droplets,and each of the individual droplets carried a color code.Droplets were then pooled together and loaded into CARMEN chip,which contained more than 150,000 microwells and each well accommodates exactly two droplets from the pool at random.Based on the color codes of the droplets,pairwise combinations of droplets containing specific Cas13a-detection mixes and droplets containing specific samples could be determined and addressed.Detection reactions were initiated by applying an electric field to merge the droplet pairs in all the wells.When Cas13a finds a target,reporter RNA was cleaved,generating fluorescence in the corresponding droplet.The CARMEN-Cas13a system was able to test more than 4,500 crRNA-target pairs on a single array,with the reagent cost per test decreased by more than 300-fold.It was able to simultaneously test dozens of samples for all 169 human-associated viruses,including the novel coronavirus SARSCoV-2.This highly scalable and multiplexed assay will greatly benefit early disease diagnosis and promote comprehensive infection surveillance,helping to improve patient care and public health.Another exciting work was introduced by Tianet al.,who combined droplet digital concept with Cas13a-based detection to realize single-molecule RNA diagnostics[105].Cas13a assay mix was partitioned in more than 20,000 picoliter droplets,and droplet involving target RNA exhibited fluorescence ater merely a 60-min isothermal incubation at 37 °C,allowing absolute quantitation of target RNA in samples based on Poisson distribution.This method enabled ultrasensitive detection of uropathogen 16S rRNAs and diagnosis of SARS-CoV-2 infection from clinical samples with attomolar sensitivity.The unique advantage of the digital Cas13a assay is the elimination of reverse transcription and nucleic acid amplification,which avoids the introduction of errors due to sample loss,cross-contamination,and amplification bias.For further improvement of this method,more work needs to be done to achieve multiplexed target detection and increase the throughput.

    5.Antibiotic susceptibility testing

    The inappropriate use and massive consumption of antibiotics have caused antibiotic resistance of pathogens,which has become a global health crisis in recent decades[106].AST is a widely used approach for evaluating the susceptibility of microorganisms to antibiotics and identifying possible antibiotic resistance[107,108].It provides essential information to determine suitable antibiotics for pathogen infection treatment.Droplet microfluidics offers a powerful and promising platform for rapid,sensitive and highthroughput AST[109-112].It enables AST down to a single-cell level,uncovering heteroresistance in an isogenic population[113].Kaushiket al.presented dropFAST platform,which confined single bacteria in picoliter droplets with antibiotics and the viable cell indicator resazurin[114].After 1 h incubation,fluorescence of droplets was measured for quantitation of growing bacteria and assessment of antibiotic susceptibility.Introduction of resazurin into the droplets at the very beginning may increase the background fluorescence of droplets.To increase the contrast in fluorescence between droplets containing resistant cells and susceptible cells,Lyuet al.developed a droplet AST system in which resazurin was injected into droplets after incubating the cells with antibiotics for 1 h(Fig.7a)[115].This system allowed to detect a resistant sub-population that comprised as low as 10-6of the entire population.To further enhance the testing throughput,Posteket al.developed a microfluidic platform that could produce a series of emulsions each containing densely packed droplets(called as tankers)and separated from each other by a third immiscible phase(Fig.7b)[116].Single bacteria were encapsulated in droplets with antibiotics.Each tanker contained different concentrations of antibiotic and the same concentration of bacteria,and it was spatially addressable by its position in the train.This system enabled single-cell AST with hundreds of replications in a range of antibiotic concentrations within a single experimental run.To realize multiplexed AST,Kanget al.presented a platform consisting of four integrated microdroplet arrays,each hosting over 8000 docking sites[117].Four combinations of antibiotics/pathogens could be screened simultaneously,providing antibiotic sensitivity assessments in 15-30 min.

    Combinatorial drug treatment is an appealing strategy to overcoming antibiotic resistance[118].However,the discovery of new combinations has been slowed due to the high complexity,cost,and compound consumption of conventional screening methods[119].Recently,the Blainey group introduced a high-throughput combinatorial drug screening platform that automatically constructs drug combinations from nanoliter-scale droplets[120].The design and working principle of this platform is the same as the aforementioned CARMEN chip.Pools of droplets containing bacteria with different drugs were loaded into the microwell array,with each microwell capturing two droplets at random.The pairwise droplets were also indexed by the color codes,and after merging,bacteria could be incubated with differently drug combinations and their growth were assayed by constitutive GFP fluorescence.This system realized the prediction of synergy between more than 4,000 investigational and approved drugs and a panel of 10 antibiotics againstE.coli,and discovered a range of drugs not previously indicated for infectious disease that synergized with antibiotics.

    To shorten the turnaround time for AST,several rapid and portable AST systems have been developed[121,122].Avesaret al.proposed a system that combined the resazurin assay with stationary nanoliter droplet arrays(SNDA)containing lyophilized antibiotics within each well,allowing rapid and scalable phenotypic AST within less than 5.5 h(Fig.7c)[123].Sample loading and bacteria confinement were achieved simply by pipette injection of bacterial suspension plug and oil plug sequentially,increasing its range of applications,such as for resource-poor settings.Yiet al.described an integrated microfluidic SlipChip device for direct AST of bloodstream pathogens without the need of overnight subculture to obtain isolated colonies[124].Bacteria were extracted and enriched from positive blood cultures by dielectrophoresis,and SlipChip strategy enabled parallel inoculation of the extracted bacteria into nanoliter-scale broth droplets to perform multiplexed ASTs simultaneously.Entropy-based image analysis was used for the characterization of bacterial susceptibility patterns,eliminating the requirement for single-cell morphological analysis and fluorescence labeling.The whole on-chip AST process took 3-8 h,and it has been applied to evaluate the susceptibility patterns ofE.coliandS.aureusagainst several broad-spectrum antibiotics.

    Droplet microfluidics have also been utilized to evaluate antiviral agents.Mashaghiet al.proposed a viral fusion assay for screening of candidate compounds[125].By using highly quantitative fluorescent probes,the fusion kinetics of influenza A virions with target liposomes in microdroplets were obtained with sub-second temporal resolution.Drug candidates on viral fusion dynamics could be characterized for information on the mechanism of action.Since viral fusion is an essential step for enveloped viruses to enter their host cells,insights into this process should facilitate the development of antiviral therapeutics.Furthermore,to determine the impact of latency reversing agents(LRAs)on HIV reactivation,a microfluidic device was developed for high-throughput encapsulation and direct characterization of HIV infected cells[126].HIV transcriptionally reactivated cells in response to LRAs,which produce unspliced RNA and multiply spliced RNA,were identified through intra-droplet PCR amplification and then sorted out for downstream quantification of genomic DNA or mRNA.Further analysis of single CD4+T cells obtained from individuals receiving antiretroviral therapy found that LRAs enhance transcription from active cells in some cases while increase the number of transcriptionally active cells in others,demonstrating the importance of single-cell analysis for deeper mechanistic understanding of HIV reactivation strategies.Liuet al.characterized molecular inhibitors of poliovirus by tracking thousands of individual virus–cell interactions in microdroplets[127,128].Single-cell analysis has the ability to reveal inhibitor-specific signatures of antiviral action.

    Fig.7.Antibiotic susceptibility testing in droplets.(a)Single-cell AST in droplet microfluidics for quantitative phenotyping of heteroresistance;scale bar = 200 μm.Reproduced with permission[115].Copyright 2018,Elsevier.(b)Generation of a series of emulsions each containing densely packed droplets for high-throughput single-cell AST;scale bar = 400 μm.Reproduced with permission[116].Copyright 2018,The Royal Society of Chemistry.(c)Stationary nanoliter droplet arrays(SNDA)combined the resazurin assay for rapid and scalable phenotypic AST.Reproduced permission[123].Copyright 2017,the National Academy of Sciences of the United States of America.

    6.Microbial interactions

    Microbial interactions are the major drivers of microbial community dynamics and functions[129-131].These interactions occur by the environmental recognition followed by transference of molecular and genetic information.Interactions patterns could be positive,negative,or negligible[129].Detecting and investigating various types of interactions in microbial ecosystems remain challenging due to the complex,highly diverse and highly dynamic nature of these communities[132].Droplet microfluidics facilitates overcoming this challenge by encapsulating and co-culturing mixtures of microbes in droplets,which allow high-throughput characterization of interactions in highly controlled environments and conditions[133-135].Tanet al.co-cultivated sub-communities from a fecal sample in droplets,supporting ecological interactions between a reduced number of cells[136].After a 1-week coculture,22 individual droplets with strong bacterial co-growth were isolated for whole genome sequencing to characterize the bacterial sub-community.Mixed cultures of microalgae in millions of parallel droplets enabled investigations into the interactions that might affect the biomass accumulation,offering a high-throughput approach to detect transgressive overyielding,which occurs when the biomass production of a community of algae species exceeds that of the most productive constitutive species[137].

    Droplet-based high-throughput isolation allows screening of probiotics with antimicrobial activity.Terekhovet al.developed an ultrahigh-throughput microfluidic platform,which co-encapsulated microbiota with GFP-labeled target pathogen in tens of thousands water-in-oil-in-water double emulsions(Fig.8a)[138,139].After co-cultivation,calcein violet AM was added to MDE for viable bacteria indication.Unlike water-in-oil droplets whose external phase is oil,double emulsions are in aqueous carrier fluid,making them compatible with the flow cytometry.Droplets demonstrating anti-pathogen features(high initial target load,low GFP level,and high metabolic activity)were isolated by fluorescenceactivated cell sorting and characterized by sequencing and liquid chromatography-mass spectrometry,providing genotype,phenotype and secretome descriptions of the encapsulated microorganisms.This platform has been successfully utilized to profile the anti-pathogen properties of an entire oral microbial community of the Siberian bear[139].An antibiotic amicoumacin A produced by aBacillus pumilusstrain was identified,which was responsible for growth inhibiting ofS.aureus.

    Gene transformation is the genetic alteration of a recipient bacterium by directly uptake of exogenous genetic material released from donor bacteria[140,141].It is one important pathway leading to antibiotic resistance of pathogens[142].Lamet al.developed a droplet microfluidic system for encapsulation of individual pneumococcal cell pairs in femtoliter droplets[143,144].Gene transformation between cell combinations was characterized using wholegenome sequencing,and the results indicated that both multiple coincident transfers and transfers of large blocks or tracts were displayed after short-term cell–cell competence interactions.

    Fig.8.Microbial interactions in droplets.(a)Ultrahigh-throughput microfluidic platform for droplet-based screening and isolation of probiotics with antimicrobial activity;scale bar = 50 μm.Reproduced with permission[139].Creative Commons Attribution 4.0 International License.(b)Microbial interaction network inference in microdroplets(MINI-Drop)to identify the interactions within microbial communities across different environments;scale bar = 40 μm.Reproduced with permission[145].Copyright 2019,Elsevier.(c)Droplet printing enabling to position interacting microbes in specific patterns in sub-millimetre arrays;scale bar = 200 μm.Reproduced with permission[148].Creative Commons Attribution 4.0 International License.

    Besides pairwise microbial interactions,higher-order interactions,which are challenging to be controlled and analyzed by conventional methods,can be deciphered using microfluidic droplets in a high-throughput manner.Hsuet al.presented microbial interaction network inference in microdroplets(MINI-Drop)to identify the interactions within microbial communities across different environments(Fig.8b)[145].Two or three bacterial strains were mixed and randomly encapsulated in droplets,enabling parallel culturing of many sub-communities.The highlight of this work was coupling fluorescence microscopy to an automated computational method,which allowed rapid determination of the absolute abundance of each strain in hundreds to thousands of droplets per condition,inferring the type and strength of microbial interactions after co-cultivation.The MINI-Drop could accurately depict pairwise and higher-order interactions in synthetic consortia with an inferred interaction network.The complex web of interactions linking antibiotics and different species was thus elucidated in a synthetic consortium.The Blainey group introduced the kChip,a platform similar to their CARMAN chip,performing rapid,massively parallel,bottom-up construction and screening of synthetic microbial communities[146].Different strains(in different environmental conditions)were distributed into nanoliter droplets and randomly self-assembled in kChip.The kChip contained a high-density array of microwells,each of which was designed to precisely group k droplets(k = 1,2,…,7,19).The grouped droplets were subsequently merged,allowing for rapid functional profiling of fluorescently labeled and unlabeled strains across libraries of conditions(e.g.,antibiotics,natural products,carbon sources)with flexible temporal resolution.The kChip was successfully applied to screen ~100,000 multispecies communities comprising up to 19 soil isolates,and specific compositions of bacteria that promoted the growth of a model plant symbiontHerbaspirillum frisingensewere identified.It is a promising technique for the discovery of microbial consorftia with functional properties,such as pathogen suppression and recalcitrant substrate degradation.

    In microbial community,the spatial arrangement of strains and species is considered critical for their ecology[147].Kumaret al.presented a high-resolution droplet printing method that enabled interacting microbes to be positioned in specific patterns in sub-millimetre arrays(Fig.8c)[148].Using piezo droplet printers,droplets containing bacteria and agarose were deposited into 2D patterns by line-by-line printing in an oil/phospholipid solution.Stabilized structures were formed after gelation and then transferred to medium for long-term cultivation.Results revealed that spatial structure impacted the outcome of both interference and exploitative competition.By printing different potential barriers between a toxin-producing strain and a susceptible strain,the ecological refuge effect was investigated,demonstrating that cells closest to the toxin producer mopped up the toxin and protected their clonemates.This method realizes precisely patterning and manipulating bacterial interactions at the relevant spatial scales,which is highly valuable for establishing a causal link between the spatial organization of bacterial strains and ecological outcomes.

    Fig.9.Microbial biotechnology in droplets.(a)High-throughput optimization of recombinant protein production in microfluidic gel droplets;scale bar = 50 μm.Reproduced with permission[158].Creative Commons Attribution 4.0 International License.(b)Droplet-based ESI-MS for microbial metabolite characterization.Reproduced with permission[159].Copyright 2020,American Chemical Society.(c)Microfluidic automated plasmid library enrichment for biosynthetic gene cluster discovery;scale bar = 100 μm.Reproduced with permission[164].Creative Commons Attribution 4.0 International License.

    7.Microbial biotechnology

    Microbial biotechnology deals with the manipulation of genetic engineering of microorganisms or their components to produce useful products for the human welfare[149].It is critical for advances in food industry,agriculture and pharmaceuticals[150-152].Droplet microfluidics provides an efficient platform for strain improvement,metabolite screening and biosynthetic gene cluster discovery in microbial biotechnology[17,153,154].Jianet al.reported a microbial microdroplet culture system(MMC),which performed automated microbial cultivation,growth monitoring and adaptive evolution in up to 200 replicate droplets[155].An 18-day process of adaptive evolution of methanol-essentialE.colistrain was conducted in MMC,obtaining two strains exhibiting higher growth rates compared with the parent strain.By co-encapsulating screening microbes with a sensor strain or metabolite capturing particles,the microbial production of metabolites can be evaluated in droplets for strain improvement.Hernandez-Valdeset al.developed a growth-based sensor strain to detect essential amino acids[156].Encapsulating potential amino acid producing strains with the sensor strain in droplets enabled to identify strains that naturally secreted amino acids.This method isolated three mutatedLactococcus lactisIPLA838 strains with 5~10-fold increased amino acid-secretion compared to the wild type.Similarly,Mahleret al.implemented a high-throughput screening of antibiotic producers from natural microbial communities by pico-injecting bacterial reporter strains(E.coliandBacillus subtilis,as representatives of Gram-negative and Gram-positive bacteria)into all droplets[157].Napiorkowskaet al.co-encapsulated single microbes with proteincapturing microparticles in agarose gel droplets(Fig.9a)[158].During microbial cultivation,the secreted proteins with hexahistidine sequences(His6-tag)could be continuously captured by the microparticles.The agarose gel beads were then recovered into an aqueous buffer,treated for protein crosslinking and labeling,and sorted by FACS.The highly productive variants could be obtained by isolating the most fluorescent fraction of agarose gel beads.This system permitted the analysis of>106 genotypes per day and was applied to optimize cutinase secretion inKomagataella phaffii,isolating a strain with 5.7-fold improvement.

    Combing droplet microfluidics with mass spectrometry allows label-free characterization of microbial products.Mahleret al.integrated droplet microfluidic device with electrospray mass spectrometry(ESI-MS)(Fig.9b)[159,160].A few of microbes were cultivated in droplets,which were then electrosprayed from an integrated steel emitter and analyzed by MS for the detection of accumulated metabolites.This system enabled analysis of the conversion of D-glucose to L-lysine byCorynebacterium glutamicum.Currently,this method only provided qualitative results of lysine production.Limited by the attenuated detection sensitivity due to the ion suppression by buffer salt,it encapsulated multiple microbes in each droplet for metabolite accumulation.Future works need to be done to realize metabolite quantitation and improve sensitivity,allowing us to disclose metabolism heterogeneity of single microbes in a complex community.

    Microbial communities provide a rich source of bioactive molecules such as enzymes,but their experimental identification is challenging.Coupling metagenomic library construction with droplet-based screening can promote the discovery of biosynthetic gene clusters from nature[159-161].Maet al.constructed a metagenomic library hosted inE.colifrom domestic running water samples[162].SingleE.colicells were then encapsulated in water-in-oil-in-water double emulsion droplets,along with the enzyme substrate.After cell lysis,the released enzyme could catalyze the substrate to generate a fluorescent product,and the droplets were sorted by FACS to enrich the positive clones.This platform enabled to identify a novel esterase with high catalytic efficiency and distinct evolutionary origin from other lipolytic enzymes.Similar system was also presented by Tauzinet al.,who constructed a metagenomic library from mucus-associated human gut microbiome and screened for host-glycan degrading enzymes[163].Results revealed novel pathways of metabolization of human glycans by gut bacteria,in particular pathobionts.To offer a more generalizable way to identify sequences of interest,Xuet al.developed a screening system which merged the colony droplets with PCR reagent droplets(Fig.9c)[164].After PCR,positive droplets were sorted directly on chip by dielectrophoretic force.DNA was then recovered from droplets for downstream sequencing to identify the target gene clusters.This system was successfully applied to isolate and sequencing type I polyketide synthase gene clusters from an Antarctic soil metagenome.

    8.Challenges and future perspectives

    Droplet microfluidics,as a versatile and powerful technique that allows high-throughput and precise manipulation of microorganisms,has expanded the scope and context of microbiology.The technical innovations of droplet microfluidics have opened new possibilities to various fields of microbiology,from defining single-cell phenotypic and genetic heterogeneity to investigating spatiotemporal dynamics of microbial communities,from precise quantitation of microbiota to systematic decipherment of microbial interactions,from isolating rare and uncultured microbes to improving genetic engineered strains.Attributed to its unique advantages,droplet microfluidics is being expected to serve as basic tools for microbiologists to explore the world of microorganisms.However,there are still challenges encountered in droplet microfluidics before its extensive adoption in mainstream microbiology.

    There are several issues related to the inherent properties of droplets.Firstly,droplet stability during droplet manipulation and microbial cultivation is highly important to maintain the isolation of microorganisms;despite having a wide variety of surfactants[165,166],it is not easy to balance the high droplet stability with low microbial cytotoxicity for microbial culture.The second issue is the limited capability of substance exchange in droplets,which throws up difficulties in long-term cultivation,especially for metabolically active microorganisms.Encapsulating microbes in giant unilamellar lipid vesicles[36,38]or hydrogel microcapsules[167-169]can solve this issue to some extent,but these methods will lead to crosstalk among compartments,and individually regulating compositions of each droplet is still hard to achieve.Thirdly,for microorganisms,the oxygen requirement varies greatly among species.Currently,the growth of aerobic microbes in droplets are mainly supported by the dissolved oxygen in oil phase(fluorinated oil such as FC-40 and HFE-7500),while anaerobic microbes in droplets require to be generated,manipulated and cultured in anaerobic chamber.Precise control and flexible regulation of oxygen level in droplets is still challenging to realize in droplet microfluidic systems.Fourthly,recent studies revealed that small molecules,particularly lipophilic molecules,can partly dissolve in oil phase and exchange among droplets[170];although some research takes advantage of this phenomenon for investigation of microbial interactions among droplets[171],it may often lead to undesirable cross-contamination,affecting the outcomes of droplet-based microbial study.Finally,but not least,the liquid environment of droplets may be not suitable for some microorganisms that need to grow on solid medium or require airwater interface[172].

    There are also some challenges toward the widespread application of microfluidics in microbiology.The first is the lack of standardized devices and procedures.A variety of droplet microfluidic platforms have been developed for microbiological studies in various fields.Owing to the complexity of experimental conditions,comparisons among systems could be rather difficult,inevitably limiting the adoption of microfluidic methods by a broader community.Currently,ddPCR has made a good example of standardization and practical transformation of novel microfluidic technologies,which will inspire further development of the microfluidic industry.Secondly,most of current droplet microfluidic systems require tedious device fabrication as well as complex fluid operations,which are challenging for non-experts to learn and use.To address this issue,several off-chip droplet systems have been recently developed,such as cross-interface emulsification(XiE)[173,174],spinning-based emulsification[175-177],vertical step emulsification[178,179]and vibrating sharp-tip capillary[180],which use commercially available“off-the-shelf”components for system setup and rely less on complicated chip manufacture and experienced operators.This issue can also be solved by exploring further automation of liquid handling protocols,improving the accessibility of droplet microfluidics to microbiologists.

    Recently,droplet microfluidic technology has become one of the most powerful and promising tools for microbiological studies.It has been extensively utilized in various fields and is still in rapid progression.There are fascinating opportunities for droplet microfluidics in exploration of novel or uncultured microbial resources,diagnosis and surveillance of pathogen infections,discovery of antimicrobial drugs or probiotics,elucidation of complex microbial interactions,and directed evolution of microbial communities.For further development,one important direction is the automation and simplification of the whole droplet microfluidic system,providing standardized equipment and routine protocols for diverse aspects of microbiological applications.Another direction is the improvements of downstream analysis techniques.Currently,molecular techniques,such as PCR,sequencing and CRISPRCas13,have been commendably integrated with droplet microfluidics to acquire genome or transcriptome information of microorganisms in droplets.However,proteome or metabolome information of microorganisms in droplets is less documented.Labelfree techniques,including single-cell Raman spectrum and mass spectrometry,are highly promising for microbial proteomic and metabolomic analysis;but their integration with droplet microfluidics still needs further optimization,to improve the detection sensitivity as well as the processing throughput.We believe that indepth collaborations between multidisciplinary researchers(microbiologists,chemists,engineers and clinicians)will promote further development and maturation of droplet microfluidics,continually introduce new paradigms for microbiology,and contribute to frontier microbial ecology.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    This work was supported by National Natural Science Foundation of China(Nos.82173774,31925037,22104041),the Fundamental Research Funds for the Central Universities(Nos.2662021DKQD001,2662021JC001),and the Cheeloo Scholar Program of Shandong University(to W.Liu).

    在线 av 中文字幕| 中文字幕精品免费在线观看视频 | 2021少妇久久久久久久久久久| 国产精品一区二区在线观看99| 丝袜在线中文字幕| 精品人妻熟女av久视频| 最近中文字幕高清免费大全6| 国产精品一区二区在线观看99| 国产免费一区二区三区四区乱码| 一区二区三区四区激情视频| 国产精品国产三级专区第一集| 乱码一卡2卡4卡精品| 国产中年淑女户外野战色| 久久99热这里只频精品6学生| √禁漫天堂资源中文www| 哪个播放器可以免费观看大片| 亚洲精品视频女| 制服丝袜香蕉在线| 五月开心婷婷网| 岛国毛片在线播放| 国产成人精品福利久久| 大又大粗又爽又黄少妇毛片口| 欧美+日韩+精品| 九九爱精品视频在线观看| 亚洲综合精品二区| 22中文网久久字幕| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久精品电影小说| 亚洲四区av| 一区二区三区乱码不卡18| 国产美女午夜福利| 最新的欧美精品一区二区| 熟女电影av网| 黑丝袜美女国产一区| 一个人看视频在线观看www免费| 成年女人在线观看亚洲视频| 日本黄色日本黄色录像| 欧美日韩在线观看h| 国产视频首页在线观看| 久久99热6这里只有精品| 男女无遮挡免费网站观看| 成人亚洲欧美一区二区av| 欧美日韩精品成人综合77777| 久久综合国产亚洲精品| 亚洲国产欧美日韩在线播放 | 男人添女人高潮全过程视频| 人妻夜夜爽99麻豆av| 永久网站在线| 亚洲天堂av无毛| 亚洲一级一片aⅴ在线观看| 18禁在线播放成人免费| 国产av国产精品国产| 国产精品久久久久成人av| 最近的中文字幕免费完整| 亚洲精品中文字幕在线视频 | 少妇的逼水好多| 亚洲av二区三区四区| 欧美区成人在线视频| 国产精品.久久久| 日韩av在线免费看完整版不卡| 最近最新中文字幕免费大全7| 久久ye,这里只有精品| 精品少妇内射三级| 毛片一级片免费看久久久久| 乱系列少妇在线播放| 欧美 日韩 精品 国产| 亚洲第一av免费看| 亚洲在久久综合| 成年女人在线观看亚洲视频| 免费黄频网站在线观看国产| 国产精品久久久久久久久免| 国产乱来视频区| 成人二区视频| 成年av动漫网址| 晚上一个人看的免费电影| 久久久久久久国产电影| 丝瓜视频免费看黄片| 中文字幕人妻熟人妻熟丝袜美| 国产一区有黄有色的免费视频| 成人美女网站在线观看视频| 午夜影院在线不卡| a级片在线免费高清观看视频| 精品国产一区二区三区久久久樱花| 亚洲怡红院男人天堂| 美女福利国产在线| 国产成人精品一,二区| 欧美日韩一区二区视频在线观看视频在线| 国产精品三级大全| 男人狂女人下面高潮的视频| 久久久久久久久久久免费av| 国产欧美另类精品又又久久亚洲欧美| 国产 精品1| 亚洲第一av免费看| 久久久久久久久久成人| 亚洲av国产av综合av卡| 99国产精品免费福利视频| 国产伦理片在线播放av一区| a级一级毛片免费在线观看| 人人妻人人看人人澡| 国产成人一区二区在线| 成年人免费黄色播放视频 | 在线观看av片永久免费下载| 日本午夜av视频| 免费播放大片免费观看视频在线观看| 日韩中字成人| 国产精品国产三级国产av玫瑰| 日本猛色少妇xxxxx猛交久久| 毛片一级片免费看久久久久| 青春草视频在线免费观看| 国产女主播在线喷水免费视频网站| 国产成人精品婷婷| 久久久久精品性色| 欧美激情国产日韩精品一区| 最黄视频免费看| 亚洲成人手机| 又爽又黄a免费视频| 国产欧美日韩一区二区三区在线 | 精品视频人人做人人爽| 国产爽快片一区二区三区| 亚洲av在线观看美女高潮| 久久久国产欧美日韩av| 人妻 亚洲 视频| 亚洲欧美日韩东京热| 精品人妻熟女毛片av久久网站| h日本视频在线播放| 伊人久久精品亚洲午夜| 80岁老熟妇乱子伦牲交| 少妇的逼水好多| 国产av精品麻豆| 人妻制服诱惑在线中文字幕| 久久99精品国语久久久| 久久97久久精品| av女优亚洲男人天堂| 国产精品麻豆人妻色哟哟久久| av女优亚洲男人天堂| 老女人水多毛片| 中文字幕人妻丝袜制服| √禁漫天堂资源中文www| 99热全是精品| av天堂中文字幕网| 日韩av免费高清视频| 国产美女午夜福利| 看十八女毛片水多多多| 少妇的逼好多水| 少妇 在线观看| 99国产精品免费福利视频| 又大又黄又爽视频免费| 熟女人妻精品中文字幕| 国产精品一区二区三区四区免费观看| 国产男人的电影天堂91| 香蕉精品网在线| 亚洲经典国产精华液单| xxx大片免费视频| 人妻少妇偷人精品九色| 久久精品国产亚洲av涩爱| 亚洲精品日韩av片在线观看| 国产老妇伦熟女老妇高清| 两个人免费观看高清视频 | 乱人伦中国视频| 天堂中文最新版在线下载| a级毛片免费高清观看在线播放| 久久精品久久久久久久性| 欧美日本中文国产一区发布| 又粗又硬又长又爽又黄的视频| 亚洲精品成人av观看孕妇| 国产在线视频一区二区| 2018国产大陆天天弄谢| 婷婷色av中文字幕| 男女免费视频国产| 91久久精品国产一区二区成人| 精品亚洲成国产av| av播播在线观看一区| 美女内射精品一级片tv| 国产成人午夜福利电影在线观看| 亚洲伊人久久精品综合| 国产日韩欧美视频二区| 黄色配什么色好看| 97超视频在线观看视频| 国产中年淑女户外野战色| 精品酒店卫生间| 各种免费的搞黄视频| 久久精品夜色国产| 伦理电影大哥的女人| xxx大片免费视频| 色视频www国产| 国产精品一区二区在线观看99| 一级片'在线观看视频| 国产亚洲一区二区精品| 国产av精品麻豆| 久久久精品94久久精品| 久久精品国产鲁丝片午夜精品| 美女大奶头黄色视频| 久久精品久久精品一区二区三区| 人妻 亚洲 视频| 精品国产乱码久久久久久小说| 美女cb高潮喷水在线观看| 欧美变态另类bdsm刘玥| 色吧在线观看| 99热全是精品| 男女免费视频国产| 精品少妇内射三级| 伦精品一区二区三区| 亚洲综合色惰| 免费看av在线观看网站| 欧美精品一区二区大全| 精品少妇内射三级| av视频免费观看在线观看| 精品久久久久久久久亚洲| 在现免费观看毛片| 伦理电影大哥的女人| 免费久久久久久久精品成人欧美视频 | 色5月婷婷丁香| 久久精品久久久久久噜噜老黄| 亚洲av免费高清在线观看| 欧美人与善性xxx| 七月丁香在线播放| 亚洲精品456在线播放app| av.在线天堂| 天堂8中文在线网| 精品亚洲成a人片在线观看| 亚洲av欧美aⅴ国产| 精品午夜福利在线看| 亚洲欧洲国产日韩| 伊人久久精品亚洲午夜| 大陆偷拍与自拍| 久久99蜜桃精品久久| 爱豆传媒免费全集在线观看| 亚洲,一卡二卡三卡| 晚上一个人看的免费电影| 国产欧美日韩精品一区二区| 久久久国产精品麻豆| 亚洲精品日韩av片在线观看| 久久人人爽人人爽人人片va| videossex国产| 久热久热在线精品观看| 中国美白少妇内射xxxbb| 国产永久视频网站| 免费久久久久久久精品成人欧美视频 | 女人久久www免费人成看片| 色哟哟·www| 亚洲精品乱码久久久v下载方式| 亚洲激情五月婷婷啪啪| 国模一区二区三区四区视频| 亚洲无线观看免费| 国国产精品蜜臀av免费| 97在线人人人人妻| 天堂中文最新版在线下载| 最新中文字幕久久久久| 成人综合一区亚洲| 尾随美女入室| 免费人妻精品一区二区三区视频| 国产精品一区二区在线不卡| 一本一本综合久久| 日本91视频免费播放| 国产av码专区亚洲av| 日韩中字成人| 精品少妇内射三级| 黄色欧美视频在线观看| 22中文网久久字幕| 91久久精品国产一区二区三区| 水蜜桃什么品种好| 亚洲精品乱码久久久v下载方式| 亚洲精品一二三| 插逼视频在线观看| 一区在线观看完整版| 一级毛片久久久久久久久女| 日韩 亚洲 欧美在线| 国产欧美日韩精品一区二区| 边亲边吃奶的免费视频| 亚洲精品色激情综合| 日韩欧美精品免费久久| 高清在线视频一区二区三区| av国产精品久久久久影院| 高清毛片免费看| 中国美白少妇内射xxxbb| 人人妻人人澡人人看| 97精品久久久久久久久久精品| 夫妻性生交免费视频一级片| 啦啦啦中文免费视频观看日本| 日韩精品免费视频一区二区三区 | 日韩一区二区视频免费看| 亚洲av成人精品一二三区| 一本大道久久a久久精品| 日日啪夜夜爽| 欧美日韩综合久久久久久| 日韩av在线免费看完整版不卡| 七月丁香在线播放| 大片免费播放器 马上看| 国产精品久久久久久久久免| 亚洲欧美日韩东京热| 亚洲欧洲国产日韩| 亚洲婷婷狠狠爱综合网| 伊人久久国产一区二区| 中文字幕人妻丝袜制服| 亚洲一级一片aⅴ在线观看| 国产精品一区www在线观看| 国产欧美日韩精品一区二区| 国产精品偷伦视频观看了| 免费大片18禁| 一级av片app| 伦理电影大哥的女人| 国产一区亚洲一区在线观看| 国产黄色视频一区二区在线观看| 赤兔流量卡办理| 18禁裸乳无遮挡动漫免费视频| 狠狠精品人妻久久久久久综合| 亚洲精品,欧美精品| a 毛片基地| 国产熟女欧美一区二区| 日韩欧美精品免费久久| 国产男女超爽视频在线观看| 卡戴珊不雅视频在线播放| 国产精品不卡视频一区二区| 成人影院久久| 亚洲欧美成人精品一区二区| 成人特级av手机在线观看| 一级爰片在线观看| 色吧在线观看| 国产色爽女视频免费观看| 99热国产这里只有精品6| 三级国产精品片| 永久免费av网站大全| 能在线免费看毛片的网站| 黄色配什么色好看| 中文字幕av电影在线播放| a级片在线免费高清观看视频| 免费久久久久久久精品成人欧美视频 | 日韩伦理黄色片| 国产视频内射| 国产极品天堂在线| 少妇熟女欧美另类| 亚洲精品亚洲一区二区| 国产中年淑女户外野战色| 亚洲av欧美aⅴ国产| 亚洲av免费高清在线观看| 国产av精品麻豆| 美女福利国产在线| 久久女婷五月综合色啪小说| 18禁在线无遮挡免费观看视频| 国产探花极品一区二区| 久久青草综合色| 中国美白少妇内射xxxbb| 久久久精品94久久精品| 日韩欧美一区视频在线观看 | 欧美老熟妇乱子伦牲交| 成人影院久久| 午夜福利,免费看| 成人影院久久| 色5月婷婷丁香| 我要看日韩黄色一级片| 国产伦理片在线播放av一区| av天堂久久9| 日韩欧美一区视频在线观看 | 久久精品国产亚洲av涩爱| 9色porny在线观看| 人人妻人人爽人人添夜夜欢视频 | 国产色婷婷99| 国产男女内射视频| h视频一区二区三区| xxx大片免费视频| 亚洲国产毛片av蜜桃av| 亚州av有码| 狠狠精品人妻久久久久久综合| 免费观看av网站的网址| 在线看a的网站| 全区人妻精品视频| 少妇人妻一区二区三区视频| 热re99久久国产66热| 国产精品麻豆人妻色哟哟久久| av网站免费在线观看视频| 18禁在线播放成人免费| 又黄又爽又刺激的免费视频.| 亚洲精品乱码久久久久久按摩| 美女xxoo啪啪120秒动态图| 亚洲精品一区蜜桃| 亚洲av欧美aⅴ国产| 高清不卡的av网站| 黄色欧美视频在线观看| 亚洲美女黄色视频免费看| 日韩欧美精品免费久久| 亚洲欧美成人精品一区二区| 99久久精品热视频| 秋霞在线观看毛片| 中文字幕人妻丝袜制服| 精品少妇久久久久久888优播| 少妇精品久久久久久久| 高清av免费在线| 成年女人在线观看亚洲视频| 一边亲一边摸免费视频| 国产又色又爽无遮挡免| 成人漫画全彩无遮挡| 人妻一区二区av| 欧美日韩视频高清一区二区三区二| 少妇被粗大的猛进出69影院 | 亚洲va在线va天堂va国产| 免费大片黄手机在线观看| 国产综合精华液| 80岁老熟妇乱子伦牲交| 2018国产大陆天天弄谢| 亚洲国产毛片av蜜桃av| 亚州av有码| 亚洲av在线观看美女高潮| 国产极品粉嫩免费观看在线 | 久久6这里有精品| 国产伦理片在线播放av一区| 偷拍熟女少妇极品色| 日韩在线高清观看一区二区三区| a级毛片免费高清观看在线播放| 男女无遮挡免费网站观看| 久久久国产一区二区| 精品久久久噜噜| 成人国产av品久久久| 亚洲欧美中文字幕日韩二区| 亚洲中文av在线| 日本av手机在线免费观看| 少妇被粗大的猛进出69影院 | 乱人伦中国视频| 精品久久久久久电影网| 国产国拍精品亚洲av在线观看| 又爽又黄a免费视频| 久久久久久久久大av| 亚洲精品一区蜜桃| 99久久中文字幕三级久久日本| 亚洲av不卡在线观看| 黄片无遮挡物在线观看| 亚洲精品国产av蜜桃| 日韩不卡一区二区三区视频在线| 老熟女久久久| 久久久久久久久大av| 国产乱人偷精品视频| 少妇丰满av| 中文字幕亚洲精品专区| 免费播放大片免费观看视频在线观看| 久久精品国产亚洲av天美| 日韩不卡一区二区三区视频在线| 三上悠亚av全集在线观看 | 香蕉精品网在线| 26uuu在线亚洲综合色| 精品亚洲乱码少妇综合久久| 久久久国产一区二区| 人人妻人人看人人澡| 黑人猛操日本美女一级片| 欧美一级a爱片免费观看看| 亚洲久久久国产精品| 国产精品不卡视频一区二区| 国产一区二区在线观看日韩| 少妇猛男粗大的猛烈进出视频| 国产一区有黄有色的免费视频| 又大又黄又爽视频免费| 一级黄片播放器| 中文字幕人妻丝袜制服| 国产日韩一区二区三区精品不卡 | 国精品久久久久久国模美| 国产成人freesex在线| 国产欧美另类精品又又久久亚洲欧美| 亚洲人成网站在线播| 蜜臀久久99精品久久宅男| 久久久久精品性色| 精品熟女少妇av免费看| 蜜臀久久99精品久久宅男| 国产老妇伦熟女老妇高清| 一级毛片黄色毛片免费观看视频| 久久97久久精品| 色5月婷婷丁香| 亚洲精品国产色婷婷电影| 成人午夜精彩视频在线观看| 丰满乱子伦码专区| 不卡视频在线观看欧美| 一级毛片我不卡| 永久免费av网站大全| 99热这里只有是精品在线观看| 国产一区二区在线观看日韩| 国产一区有黄有色的免费视频| 欧美另类一区| 国产欧美亚洲国产| 日本欧美国产在线视频| 欧美亚洲 丝袜 人妻 在线| 久久久久久久国产电影| 亚洲国产欧美在线一区| 在线播放无遮挡| 久久97久久精品| 高清午夜精品一区二区三区| 99热6这里只有精品| 观看av在线不卡| 极品少妇高潮喷水抽搐| a级毛片免费高清观看在线播放| 亚洲在久久综合| 日韩欧美精品免费久久| 亚洲精品乱码久久久久久按摩| 久热久热在线精品观看| 亚州av有码| 日韩在线高清观看一区二区三区| 国产精品久久久久久av不卡| 天天操日日干夜夜撸| 韩国高清视频一区二区三区| 色婷婷av一区二区三区视频| 国产在线男女| 黄色毛片三级朝国网站 | 少妇 在线观看| 看非洲黑人一级黄片| 一级爰片在线观看| 人人澡人人妻人| 人体艺术视频欧美日本| 在线免费观看不下载黄p国产| 免费观看a级毛片全部| 极品少妇高潮喷水抽搐| 狂野欧美激情性xxxx在线观看| 少妇高潮的动态图| 美女福利国产在线| 亚洲经典国产精华液单| 亚洲高清免费不卡视频| 黄色一级大片看看| 精品少妇久久久久久888优播| 80岁老熟妇乱子伦牲交| 亚洲欧洲日产国产| 日本欧美视频一区| 十分钟在线观看高清视频www | 插逼视频在线观看| 最后的刺客免费高清国语| 老司机影院毛片| 人妻人人澡人人爽人人| .国产精品久久| videossex国产| 人体艺术视频欧美日本| 人妻人人澡人人爽人人| 中文字幕人妻丝袜制服| 国产精品熟女久久久久浪| 国产免费福利视频在线观看| 亚洲av免费高清在线观看| 三级国产精品欧美在线观看| 亚洲精品日本国产第一区| 亚洲精品aⅴ在线观看| 久久久久网色| 国产精品一区二区性色av| 男人爽女人下面视频在线观看| 交换朋友夫妻互换小说| 久久久久网色| 亚洲性久久影院| 国产极品天堂在线| 国产美女午夜福利| 91久久精品电影网| 色视频www国产| 日韩av在线免费看完整版不卡| 国内精品宾馆在线| 欧美人与善性xxx| 亚洲精品久久午夜乱码| 中国美白少妇内射xxxbb| 9色porny在线观看| 国产一区二区三区综合在线观看 | 国产精品一区二区三区四区免费观看| 日日爽夜夜爽网站| h日本视频在线播放| 午夜老司机福利剧场| 亚洲电影在线观看av| 国产精品嫩草影院av在线观看| 国产精品一区二区性色av| 精品人妻熟女av久视频| 这个男人来自地球电影免费观看 | 久久6这里有精品| 久久女婷五月综合色啪小说| 亚洲成人手机| 国产精品久久久久久久电影| 热99国产精品久久久久久7| 蜜臀久久99精品久久宅男| 91成人精品电影| 视频中文字幕在线观看| 午夜福利网站1000一区二区三区| 熟妇人妻不卡中文字幕| 噜噜噜噜噜久久久久久91| 亚洲av二区三区四区| 少妇熟女欧美另类| 一本—道久久a久久精品蜜桃钙片| 中文资源天堂在线| 免费观看的影片在线观看| 国产精品无大码| 青春草视频在线免费观看| 中文字幕av电影在线播放| 男女啪啪激烈高潮av片| av在线app专区| 少妇裸体淫交视频免费看高清| 三级国产精品欧美在线观看| 国产爽快片一区二区三区| 国产亚洲最大av| 你懂的网址亚洲精品在线观看| 一级av片app| 午夜久久久在线观看| 精品人妻熟女毛片av久久网站| 国产精品久久久久久精品电影小说| 亚洲欧洲国产日韩| 亚洲欧美成人精品一区二区| 大香蕉97超碰在线| 少妇丰满av| 日韩在线高清观看一区二区三区| 久久国内精品自在自线图片| 亚洲天堂av无毛| 久热这里只有精品99| 亚洲婷婷狠狠爱综合网| 国产精品成人在线| 大又大粗又爽又黄少妇毛片口| 国产精品久久久久久av不卡| 亚洲真实伦在线观看| 黄色欧美视频在线观看| 亚洲成人手机| 国产午夜精品久久久久久一区二区三区| 日日撸夜夜添| 免费黄频网站在线观看国产| 黑人猛操日本美女一级片| 久久鲁丝午夜福利片| 免费观看的影片在线观看| 亚洲一区二区三区欧美精品| 99热这里只有精品一区| 伦理电影免费视频| 亚洲欧美精品自产自拍| 自线自在国产av| 久久综合国产亚洲精品|