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

    Transcriptional responses to starvation of pathogenic Vibrio harveyi strain DY1*

    2020-03-19 12:31:44LIUXiaodanGAOXiaojianCHENNanZHANGYingyingLIXixiZHANGYueZHANGXiaojun
    Journal of Oceanology and Limnology 2020年2期

    LIU Xiaodan, GAO Xiaojian, CHEN Nan, ZHANG Yingying, LI Xixi, ZHANG Yue,ZHANG Xiaojun

    College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China

    Abstract Vibrio harveyi is a pathogen of various aquatic organisms that has been recently associated with massive mortality episodes in the aquaculture industry. Recurrent outbreaks of vibriosis are closely correlated with the capacity of this bacterial species to survive long-term starvation conditions. To study the regulation mechanism of gene expression at the transcriptional level in V. harveyi under starvation conditions, the transcriptomic response prof iles were determined of the Portunus trituberculatus pathogen V. harveyi strain DY1 under normal conditions and after four weeks of starvation. A total of 4 679 and 4 661 genes were expressed in the non-starved and starved cells, respectively. The signif icantly diff erentially expressed genes (DEGs) between non-starved and starved groups were identif ied, in which 255 genes were up-regulated and 411 genes were down-regulated. GO analysis and KEGG enrichment analysis were used to analyze the DEGs and revealed the involvement of these DEGs in many pathways, including ABC transporters, f lagellum assembly, and fatty acid metabolism. Several DEGs were randomly selected and their expression levels were conf irmed by quantitative real-time PCR (qRT-PCR). This is the f irst comprehensive transcriptomic analysis of starvation eff ects in V. harveyi. Our f indings will facilitate future study on stress adaptation and survival mechanisms of V. harveyi.

    Keyword: Vibrio harveyi; starvation stress; transcriptome sequencing; diff erentially expressed genes;adaptation and survival mechanisms

    1 INTRODUCTION

    Vibrioharveyiis a facultatively anaerobic,bioluminescent, and Gram-negative rod-shaped bacterium (Johnson and Shunk, 1936). Bacterial infection causes serious vibriosis in marine f ish and invertebrates (Austin and Zhang, 2006), particularly in juvenile populations (Pujalte et al., 2003), leading to massive deaths and major economic losses in aquaculture.V.harveyihas been isolated from multiple sources including common snook(Centropomusundecimalis) in USA (Kraxberger-Beatty et al., 1990), silvery black porgy(Acanthopagruscuvieri) and cultured brown spotted grouper (Epinepelustauvina) in Kuwait (Saeed,1995), and sunf ish (Molamola) in Spain (Hispano et al., 1997).V.harveyi-associated diseases are major constraint on the production of marine invertebrates,particularly in South America and Asia.V.harveyiinfections have also been reported in common dentex(DentexdentexL.) (Company et al., 1999), salmonids(SalmosalarL.) (Zhang and Austin, 2000), seahorse(Hippocampuskuda) (Alcaide et al., 2010), farmed sole (Soleasenegalensis) (Zorrilla et al., 2003),summer f lounder (Paralichthysdentatus)(Gauger et al., 2006), and gilthead sea bream (Sparusaurata)(Haldar et al., 2010). In addition,V.harveyihas been isolated from several crustacean larvae, including black tiger prawn (Penaeusmonodon) (Lavilla-Pitogo et al., 1990), kuruma prawns (Penaeusjaponicas)(Liu et al., 1996), rock lobster (Jasusverreauxi)(Diggles et al., 2000), and oysters (Ortigosa et al.,1994). Infected animals became inappetent and developed inf lammation and necrotic subdermal cysts in various species of f ishes, such as crustaceans and mollusks (Austin, 2010).

    Bacteria could enter viable but non-culturable(VBNC) or normal dormant state by a variety of environmental stresses (Wolf and Oliver, 1992;Biosca et al., 1996; Jiang and Chai, 1996). Moreover,nutrient def iciency is the most common environmental stress in natural ecosystems for microorganisms. It is reported that some species could survive for a long time in response to starvation condition such asV.harveyi,V.cholerae,V.f ischeri,V.anguillarum,V.campbellii,V.mimicus,V.proteolyticus,V.vulnif icus(Jiang and Chai, 1996; McDougald et al.,1998). Like other bacterial species,V.harveyimay be chronically exposed to nutrient-poor waters and its ability to survive in water for extended periods of time is an adaptive strategy against stress. Imaging showed that starved cells alter their morphology from rod to spherule cells (Novitsky and Morita, 1976; Sun et al., 2016). Starved cells will adjust their gene expression to adapt to environmental stress, with potential changes to transcriptional regulators, heat shock proteins, virulence regulators, and genes related to metabolism. However, little is known about gene expression prof iles under starvation stress.

    Hence, we obtained whole transcriptomic prof iles ofV.harveyiDY1 strain under starvation stress by high-throughput sequencing and determined the gene expression changes compared to non-starved cells.This is the f irst comprehensive transcriptomic analysis of starvation eff ects in pathogenicV.harveyi, and the results greatly strengthen our understanding of the molecular pathogenic and transmission mechanisms ofV.harveyistrain DY1.

    2 MATERIAL AND METHOD

    2.1 Bacterial strain

    TheV.harveyiDY1 strain was f irst isolated from cultured megalopa of swimming crab (Portunus trituberculatus) in 2011 by Zhang et al. (2014) in Jiangsu, China. DY1 strain was recovered from -80°C as glycerol stocks and incubated with LB broth(Hopebio Qingdao, China) at 28°C overnight with shaking.

    2.2 Starvation stress and cell enumeration

    Overnight cultures were diluted (1:100) with fresh LB broth and grown to the mid-log phase at 28°C.Cultures were collected by centrifugation at 12 000×gfor 10 min. The pellets were washed twice with sterile artif icial seawater to rinse nutrients off and suspended in the starvation regimes. Then the cells in triplicate were placed in Erlenmeyer f lasks containing 100 mL sterile artif icial seawater to give a f inal concentration of 108CFU/mL. Cells were incubated in the dark at 25°C and monitored for a period of four weeks according to the previous method. Cells in starvation regimes were sampled at 1, 7, 14, 21, and 28 days and counted by plate count method (Sun et al., 2016). The cell number was counted and converted to base-10 logarithms to f it normal distribution model.

    2.3 Morphology analysis

    The morphology changes ofV.harveyiDY1 strain after four weeks of starvation was observed by SEM as previously described (Arias et al., 2012). The nonstarved and starved cells were centrifuged at 10 000×gfor 10 min at 4°C. The pellets were washed with PBS,and f ixed with glutaraldehyde (0.25 g/L). A graded ethanol series was used to dehydrate the samples,then coated with gold palladium alloy in an Electron Microscopy Science (EMS 550X), and examined by Zeiss EVO 50 (Zeiss, Germany).

    2.4 Library construction and Illumina sequencing

    Cells starved for four weeks and the non-starved cells were subjected to total RNA isolation and cDNA preparations. Three independent trials were performed for sequencing. Total RNA was extracted from cells using the EasyPure RNA kit (TransGen Biotech,Beijing, China). rRNA was removed after total RNA was collected from prokaryocyte. Fragmentation buff er was added for splitting mRNA to short fragments. Taking these short fragments as templates,random hexamer-primer were used to synthesize the f irst-strand cDNA and the second-strand cDNA was synthesized using buff er, RNase H, dNTPs, DNA polymerase I. Short fragments were purif ied with QiaQuick PCR extraction kit and resolved with elution buff er for end reparation and adding poly (A).After that, the short fragments were connected with sequencing adapters and the UNG enzyme was used to degrade the second-strand cDNA, and the product was purif ied by MiniElute PCR purif ication kit before PCR amplif ication. At last, the library was sequencedusing Illumina HiSeq2000. Images generated by sequencers were converted by base calling into nucleotide sequences, which were called raw data or raw reads and are stored in FASTQ format. Raw reads produced from sequencing machines contain dirty reads that contain adapters, unknown or low quality bases. After removing reads with adaptors and low quality reads, clean reads were obtained. The clean reads were mapped to reference genome and genes sequences respectively using SOAP2 (Li et al., 2009).

    Table 1 Primers used for the detection of DEGs by qRT-PCR

    2.5 Gene annotation and diff erence analysis

    After data f iltering, we compared the clean reads with the reference genome and sequence which has been deposited at DDBJ/EMBL/GenBank under the accession numbers CP009467.2 and CP009468.1.Gene coverage is the percentage of a gene covered by reads. This value equals to ratio of the number of bases in a gene covered by unique mapping reads to number of total bases in that gene. The threshold of theP-value in multiple tests was determined by setting FDR (false discovery rate) at 0.000 1. The FPKM (fragments per kb per Million reads) method was applied to calculate gene expression (Mortazavi et al., 2008). To screen the DEGs in DY1 cells potentially related to starvation stress, statistical and analysis of gene expression was applied using fold change >2 andP≤0.001 as standards.

    2.6 Functional annotation and GO/KEGG enrichment analysis

    For gene ontology analysis, the alignment results were parsed for assigning GO terms by using Blast2GO software. Gene sequences ofV.harveyiDY1 strain were aligned against Gene Ontology database. The calculatedP-value went through Bonferroni Correction, taking correctedP≤0.05 as a threshold. GO terms fulf illing this condition are def ined as signif icantly enriched GO terms in DEGs.

    KEGG pathway enrichment analysis identif ies signif icantly enriched metabolic pathways or signal transduction pathways in DEGs comparing with the whole genome background. The calculating formula is the same as that in GO analysis.

    2.7 Quantitative real-time PCR analysis

    RNA were extracted using an EasyPure RNA kit(TransGen Biotech, Beijing, China) and reversetranscribed by TransScript One-Step gDNA Removal and cDNA Synthesis Supermix (TransGen Biotech,Beijing, China) according to the manufacturer’s instructions. The cDNA was synthesized using anchored oligo (dT) 18 primer and incubated for 15 min at 42°C. Real-time PCR was performed with 10 μL SYBR?SuperMix, 1 μL diluted cDNA, 0.4 μL of forward primer, 0.4 μL of reverse primer and 8.2 μL of nuclease-free water. Reactions were incubated at 95°C for 10 min, followed by 40 cycles of 95°C for 10 s, 60°C for 30 s. Each sample was run in triplicate for analysis. The relative gene expression were calculated using 2-ΔΔCtmethod and16S rRNA gene was chosen as an internal control for normalization.Specif ic primers were designed according to the corresponding sequences ofV.harveyi(Table 1).

    2.8 Statistics analysis

    The data are presented as the mean±SD (n=3).One-Way Analysis of variance (one-way ANOVA)was used to evaluate the diff erential expression by SPSS (20.0).

    Fig.1 SEM micrographs of V. harveyi

    Table 2 Summary of the reads and annotated genes in nonstarved group and starved group

    3 RESULT

    3.1 Morphology analysis

    Vibrioharveyiwas found decreased in size and changed in shape from a rod to a sphere after four weeks of starvation shown in SEM. There was also a signif icant decrease in size of cells after starvation compared to the initial population (Fig.1).

    3.2 De novo assembly and gene annotation

    In this study, RNA-Sequencing was used to study the adaptation ofV.harveyiunder four-week starvation using Illumina HisSeqTM2000 system.After f iltering the raw data, 14 422 512 and 8 662 506 clean reads were obtained in control and starvation stress group (Table 2). Proportions of clean reads mapped back to genome and genes can provide an overall assessment of the sequencing. The reads were aligned toV.harveyireference genome using the SOAPaligner/soap2 software and mismatches no more than 5 bases were allowed in the alignment.There were 12 375 624 (85.81%) and 6 846 716(79.04%) clean reads were mapped to reference genome in control and starvation stress group (Table 2). There were 5 605 402 (38.87%) and 3 477 390(40.14%) reads and ratio of perfect matches in total mapped reads in the control and starvation stress group (Table 2). There were 5 337 genes annotated in the transcriptome analysis.V.harveyistrain DY1 revealed 4 679 and 4 661 annotated genes before and after starvation stress, respectively.

    Fig.2 The number of DEGs in control and starvation stress groups

    3.3 Gene expression diff erence analysis

    A total of 666 diff erentially expressed genes(DEGs) were identif ied in the samples from the fourweek starved group compared to the samples from the untreated group, including 411 down-regulated genes and 255 up-regulated genes (Fig.2). Thus, about 14.2% of the total number of genes showed signif icantly altered expression levels after starvation stress for four weeks. The up-regulated DEGs included 20 genes with more than 10-fold expression change between the starved group and the untreated group. Similarly, of the down-regulated DEGs, there were four genes with more than 10-fold expression change.

    Fig.3 GO analysis of diff erential expression genes between the initial and the four-week starved V. harveyi

    3.4 GO functional classif ication of DEGs

    GO analyses were performed to classify gene functions ofV.harveyiunder starvation stress. GO analysis of DEGs inV.harveyiafter starvation revealed that in the biological process catalogue, the most abundant category of the diff erentially expressed genes control metabolism. There were 136 upregulated genes assigned metabolic process, such as hydroxyethylthiazole kinase and peptidase. There were 170 downregulated genes assigned metabolic process including hemolysin, oxidoreductase,6-phosphogluconate dehydrogenase, et al. In cellular component, the top three categories were cell (55 upregulated, 88 down regulated), cell part (53 upregulated, 86 down regulated) and membrane process (51 upregulated, 64 down regulated). In the molecular function catalogue, the most abundant terms of the diff erential expression gene were catalytic activity, followed by binding and transporter activity process, primarily including cell division protein ZapB, f lagellar motor switch protein FliM, sulfur transfer protein TusE and transcriptional regulator(Fig.3).

    3.5 KEGG enrichment analysis of DEGs

    Expression of some genes plays a signif icant role in protection against environmental stress. KEGG enrichment analysis was used to identify the most enriched pathways inV.harveyiafter starvation stress. In KEGG enrichment, the DEGs were assigned into 120 pathways. The most abundant pathways were metabolic pathways, microbial metabolism in diverse environments and Biosynthesis of secondary metabolites. Flagellar assembly pathway and Butirosin and neomycin biosynthesis pathway were the most representative pathways related to starvation stress (Fig.4). Moreover, ABC transporters pathway,f lagellar assembly pathway, fatty acid metabolism pathway, fatty acid degradation pathway, biosynthesis of antibiotics pathway and beta-Lactam resistance pathway also involved many diff erentially expressed genes.

    3.6 Verif ication of the DEGs by qRT-PCR

    To validate the sequencing results, eight diff erentially expressed genes (danK,f imT,hcp,impA,aphB,rnfA,fur,htpG) ofV.harveyiwere randomly selected for qRT-PCR. Among these genes,four genes (danK,f imT,impA,aphB) are up-regulated genes and the other four (hcp,rnfA,fur,htpG) are down-regulated genes. The cell samples that subjected for qRT-PCR are the same as those used for transcriptome sequencing. The expression prof ile of these DEGs are similar to those obtained in sequencing results (Fig.5).

    Fig.4 KEGG pathway enrichment analysis of the DEGs between the initial and the four-week starved V. harveyi

    4 DISCUSSION

    High mortality rates in the early stages of the marine f ish and invertebrates caused byV.harveyiis one of the most important reasons for economic losses during the production period. However, the survival strategies ofV.harveyiin the aquatic environment are not well understood. Microorganisms must adapt to environmental stress to grow and survive (Poindexter,1981; Kjellerberg et al., 1983; Kunttu et al., 2009).Previous studies showed that the most widespread environmental stress for microbes in a natural environment is a lack of nutrients (Wai et al., 1999;Vatsos et al., 2003; Suzina et al., 2004; Montánchez et al., 2014; Kaberdin et al., 2015; Parada et al., 2016;Montánchez et al., 2019). To study the physiological response and the long-term survival mechanism upon starvation stress is urgent and important.

    In this study, we compared the size of the starved cells with cells from the control group, and observed that the starved cells were obviously shorter (Fig.1).This implied thatV.harveyilikely triggered specif ic regulatory mechanisms to adjust its shape and size.This result was previously reported (Sun et al., 2016).Parada et al. (2016) found that morphological change was not directly related to the entry of cells into the VBNC state, although they observed a gradual decrease in the size ofV.harveyicells during incubation. Montánchez et al. (2014) found that incubation in cold seawater for 12 h did not cause any signif icant morphological changes inV.harveyi,suggesting thatV.harveyilikely elicits specif ic adaptation mechanisms maintaining its culturability under stress conditions.

    Fig.5 Verif ication of the diff erentially expressed genes of V.harveyi before and after starvation stress by RNASeq and qRT-PCR

    The reads were mapped to the reference genome.The gene coverage was above 70% for both groups.The RNA-seq data can be used to assess the variation in expression of virulence-related and metabolismrelated factors inV.harveyi, which suggest the strategies used by these bacteria to survive under starvation conditions. There were 255 up-regulated genes and 411 down-regulated genes under starvation stress, and these genes very likely act in the response to starvation stress (Fig.2). The analysis shown in Fig.3 revealed 35 GO categories relevant to starvation stress ofV.harveyi. The main terms associated with starvation stress were binding and transporter activity,metabolic process, and cell and catalytic activity.Previous proteomic study revealed that the level of membrane proteins participating in cellular transport,maintenance of cell structure, and bioenergetics processes remained unchanged during starvation at low temperature, suggesting thatV.harveyimight need these proteins for long-term survival or for the resuscitation process after dormancy (Parada et al.,2016).

    The analysis of KEGG enrichment may help us to understand the molecular mechanisms utilized byV.harveyifor long-term survival under starvation stress. The most abundant KEGG pathways included antibiotic biosynthesis, ABC transporters, butanoate metabolism pathway, starch and sucrose metabolism,and beta lactam resistance pathways, which are processes that have previously been related to starvation stress (Higgins, 2001; Kim et al., 2013;Svensson et al., 2014). Genes related to nutrient transport, such as phospholipids, oligopeptides, and simple sugars, were down-regulated after four weeks starvation. However, genes related to nutrient transfer,such as the phospholipid transfer protein gene and the monosaccharide transporter gene, were signif icantly up-regulated under starvation stress, representing an active metabolism pathway inV.harveyi. This change may be supplementation sparked by the compensatory response, and changes in the expression of genes controlling biosynthesis of lipids and molecular transport likely aff ect the composition and properties of theV.harveyicell envelope, an obvious adaptation to stress. We suspect thatV.harveyicomprehensively utilized multiple strategies and adaptation mechanisms to sustain key physiological functions under nutritional stress. At least in part,V.harveyicompensated for the reduced expression of biosynthetic genes by the upregulation of transporter genes controlling the uptake of amino acids (e.g. amino acid ABC transporters). Previous research found that limitation of nutrients leads to the signif icant downregulation of genes controlling central carbon metabolism,biosynthesis of lipids, amino acids, and nucleotides(Montánchez et al., 2014).

    To better understand the survival strategies ofV.harveyiunder starvation stress, the fatty acid metabolism pathway was analyzed. The genes related with fatty acid synthesis were signif icantly downregulated and the expression level of fatty acid degradation related genes were signif icantly increased,producing a large number of acetyl coenzyme A for TCA cycle. The observation of the up-regulation of fatty acid degradation correlates to previous work by Kaberdin et al. (2015). Moreover, the f ine-tuning of these metabolic pathways appears to be attained through the action of small regulatory RNAs known for their essential roles in post-transcriptional control of gene expression (Kaberdin and Bl?si, 2006).

    Virulence genes played a crucial role in environmental adaptation and virulence. Previous study analyzedV.harveyiadaptation in sea water microcosms at elevated temperature and found that elevated temperature also aff ected regulation ofV.harveyigenes controlling its virulence and ancillary mechanisms (i.e. production and secretion of virulence factors, biof ilm formation and motility)(Montánchez et al., 2019). The sequencing analysis revealed downregulated expression levels of virulence regulator genetoxRand f lagellin A. These virulence genes regulate a considerable number of genes involved in environmental adaptation and virulence.Other down-regulated genes include f lagellinassociated genes, such asf lgB,f lgC,f lgD,f lgF,f lgG,f lgH,f lgI,f lgK,f lgL,f liM, andf liO. Down-regulation of f lagellin-associated genes will directly reduce f lagellar synthesis ofV.harveyi, resulting in reduced motility and adhesive strength, which may be critical during the initials steps of infection. Thus, the pathogenicity may decrease under starvation stress.The decreased expression of f lagellin-associated genes suggests that under starvation conditions,V.harveyireduced energy requirements by reducing the synthesis of non-essential structures.

    Under starvation stress, many genes were enriched in metabolic pathways involved in the synthesis of macrolides and ketolide antibiotics (Fig.4). The sequencing revealed up-regulated expression levels of many genes related to the synthesis of macrolides and ketolides and beta lactam resistance genes,includingTolC(outer membrane protein),acrA(membrane fusion protein, multidrug effl ux system),andopp(oligopeptide transport system). These resistance-related genes may signif icantly improve the survival rate ofV.harveyi.

    5 CONCLUSION

    This study presented the transcriptomic response prof iles of thePortunustrituberculatuspathogenV.harveyistrain DY1 under starvation conditions, to elucidate the changes inV.harveygene expression due to starvation in seawater.V.harveyimay survive under starvation conditions by regulating the expression of virulence and metabolism-related genes. The identif ied genes may be important highvalue drug targets, suggesting new ways to eff ectively control clinical infection ofV.harveyi. Our f indings should facilitate future study on stress adaptation and environmental survival mechanisms ofV.harveyi.

    6 DATA AVAILABILITY STATEMENT

    All sequence data that support the findings of this study have been deposited in the NCBI Short Read.The sequence read archive (SRA) accession number:PRJNA507871.

    在线免费观看不下载黄p国产| 久久97久久精品| 黄色一级大片看看| 欧美日韩一区二区视频在线观看视频在线| 只有这里有精品99| 在线观看av片永久免费下载| h日本视频在线播放| 黑丝袜美女国产一区| 在线天堂最新版资源| a级毛片在线看网站| 在线观看免费日韩欧美大片 | 久久女婷五月综合色啪小说| 高清黄色对白视频在线免费看 | 三级国产精品片| 老熟女久久久| 欧美一级a爱片免费观看看| 亚洲精品国产av蜜桃| 黄色欧美视频在线观看| 国产精品人妻久久久久久| 日日摸夜夜添夜夜添av毛片| 不卡视频在线观看欧美| 女人久久www免费人成看片| 国产男女超爽视频在线观看| 美女脱内裤让男人舔精品视频| 亚洲av欧美aⅴ国产| 99re6热这里在线精品视频| 蜜桃在线观看..| 亚洲成色77777| 欧美区成人在线视频| av在线观看视频网站免费| 国产欧美日韩综合在线一区二区 | 大片免费播放器 马上看| 人妻少妇偷人精品九色| 18禁动态无遮挡网站| 十八禁高潮呻吟视频 | 国产一区二区在线观看日韩| 欧美日韩在线观看h| 久久99精品国语久久久| 日韩一本色道免费dvd| 国产成人精品无人区| 免费av不卡在线播放| 亚洲无线观看免费| 国产欧美日韩精品一区二区| 三上悠亚av全集在线观看 | av.在线天堂| 黄色日韩在线| 国产在线视频一区二区| 欧美高清成人免费视频www| 亚洲三级黄色毛片| 十八禁高潮呻吟视频 | 日产精品乱码卡一卡2卡三| 亚洲av日韩在线播放| 波野结衣二区三区在线| 亚洲国产精品999| 少妇精品久久久久久久| 大片免费播放器 马上看| av免费在线看不卡| 亚洲av欧美aⅴ国产| 日韩熟女老妇一区二区性免费视频| 日韩在线高清观看一区二区三区| 日韩电影二区| 日本91视频免费播放| 国产91av在线免费观看| 妹子高潮喷水视频| 秋霞在线观看毛片| 精品一区二区三区视频在线| 亚洲综合精品二区| 亚洲国产最新在线播放| 嫩草影院新地址| 久久午夜综合久久蜜桃| 色哟哟·www| 在线观看三级黄色| 九九久久精品国产亚洲av麻豆| 少妇人妻一区二区三区视频| 美女大奶头黄色视频| 日本黄大片高清| 高清午夜精品一区二区三区| av福利片在线观看| 欧美日本中文国产一区发布| 女人久久www免费人成看片| 久久人妻熟女aⅴ| 免费看光身美女| 欧美激情国产日韩精品一区| 欧美xxⅹ黑人| 永久免费av网站大全| 在线观看免费日韩欧美大片 | 亚洲伊人久久精品综合| 精品久久国产蜜桃| 国产精品女同一区二区软件| 啦啦啦视频在线资源免费观看| 人妻人人澡人人爽人人| 2021少妇久久久久久久久久久| 男人和女人高潮做爰伦理| 亚洲精品日本国产第一区| 一本一本综合久久| 2021少妇久久久久久久久久久| 男人和女人高潮做爰伦理| 国产淫片久久久久久久久| 婷婷色综合www| 色婷婷av一区二区三区视频| 精华霜和精华液先用哪个| 最新中文字幕久久久久| 欧美日韩精品成人综合77777| 欧美日韩在线观看h| 国产精品无大码| 国产精品无大码| 日本av免费视频播放| 只有这里有精品99| 国产精品国产三级专区第一集| 亚洲精品视频女| 久久免费观看电影| 婷婷色麻豆天堂久久| 日韩av在线免费看完整版不卡| 国内少妇人妻偷人精品xxx网站| 在线看a的网站| 天堂中文最新版在线下载| 精品国产露脸久久av麻豆| 你懂的网址亚洲精品在线观看| 精品国产一区二区三区久久久樱花| 免费av不卡在线播放| 人人妻人人澡人人爽人人夜夜| 国产成人精品福利久久| 国产白丝娇喘喷水9色精品| 久久久国产精品麻豆| 黄色日韩在线| 99视频精品全部免费 在线| 熟女电影av网| 国产免费一级a男人的天堂| 永久免费av网站大全| 成人毛片a级毛片在线播放| 男的添女的下面高潮视频| 欧美精品一区二区大全| 两个人的视频大全免费| 寂寞人妻少妇视频99o| 成年美女黄网站色视频大全免费 | 免费观看av网站的网址| 久久99蜜桃精品久久| 丰满人妻一区二区三区视频av| 观看美女的网站| 下体分泌物呈黄色| 国产精品麻豆人妻色哟哟久久| 精品人妻熟女毛片av久久网站| 永久免费av网站大全| 少妇人妻 视频| 久久ye,这里只有精品| 午夜91福利影院| 国产日韩欧美视频二区| 国产亚洲av片在线观看秒播厂| 伦精品一区二区三区| 男女免费视频国产| 少妇人妻 视频| 搡老乐熟女国产| 久久久久精品性色| 久久精品国产鲁丝片午夜精品| 亚洲久久久国产精品| 黄色日韩在线| 午夜福利在线观看免费完整高清在| 亚洲国产精品一区二区三区在线| 国产成人精品无人区| 国产精品秋霞免费鲁丝片| 少妇人妻 视频| av天堂久久9| 中文字幕精品免费在线观看视频 | 少妇被粗大猛烈的视频| 多毛熟女@视频| 天美传媒精品一区二区| 六月丁香七月| 简卡轻食公司| 欧美日韩亚洲高清精品| 久久久a久久爽久久v久久| 免费人成在线观看视频色| 亚洲精华国产精华液的使用体验| 看免费成人av毛片| 夜夜看夜夜爽夜夜摸| 亚洲精品亚洲一区二区| 国产熟女午夜一区二区三区 | 极品人妻少妇av视频| 国产国拍精品亚洲av在线观看| 18禁在线播放成人免费| 久久午夜福利片| 欧美日韩av久久| 晚上一个人看的免费电影| 97超碰精品成人国产| 免费观看av网站的网址| 99久久精品热视频| 成人影院久久| 高清欧美精品videossex| 欧美3d第一页| 人妻 亚洲 视频| av福利片在线观看| 观看av在线不卡| 又大又黄又爽视频免费| 少妇人妻久久综合中文| xxx大片免费视频| 搡老乐熟女国产| 亚洲色图综合在线观看| 精品国产国语对白av| 亚洲欧美清纯卡通| 亚洲国产精品成人久久小说| 免费在线观看成人毛片| 国产成人一区二区在线| 精华霜和精华液先用哪个| 91在线精品国自产拍蜜月| 热re99久久精品国产66热6| 欧美精品人与动牲交sv欧美| 亚洲怡红院男人天堂| 亚洲成人手机| 高清在线视频一区二区三区| av.在线天堂| 九色成人免费人妻av| 麻豆成人午夜福利视频| 在线观看一区二区三区激情| 国产白丝娇喘喷水9色精品| 99热网站在线观看| 国产精品久久久久成人av| 欧美一级a爱片免费观看看| 国产av精品麻豆| 菩萨蛮人人尽说江南好唐韦庄| 亚洲婷婷狠狠爱综合网| 亚洲国产成人一精品久久久| 伦理电影大哥的女人| 九九久久精品国产亚洲av麻豆| 午夜福利视频精品| 狠狠精品人妻久久久久久综合| 久久精品熟女亚洲av麻豆精品| 国产成人freesex在线| 国产永久视频网站| 日韩中字成人| 国产亚洲午夜精品一区二区久久| 国精品久久久久久国模美| 国产黄色视频一区二区在线观看| 91成人精品电影| 高清视频免费观看一区二区| 黑丝袜美女国产一区| 午夜av观看不卡| 中国三级夫妇交换| 久久综合国产亚洲精品| 久久国内精品自在自线图片| 热re99久久精品国产66热6| 国产欧美日韩精品一区二区| 熟女电影av网| 国产伦精品一区二区三区视频9| 精品一区在线观看国产| 久久午夜福利片| 国产乱来视频区| 免费黄网站久久成人精品| 日韩在线高清观看一区二区三区| 久久人妻熟女aⅴ| 成人国产麻豆网| 久久影院123| 国产精品久久久久成人av| 亚洲国产欧美日韩在线播放 | 99热网站在线观看| 美女xxoo啪啪120秒动态图| 日韩亚洲欧美综合| 亚洲人成网站在线播| 成人午夜精彩视频在线观看| 亚洲中文av在线| 亚洲一区二区三区欧美精品| 日日撸夜夜添| av在线播放精品| 亚洲欧美一区二区三区国产| 亚洲不卡免费看| 国产淫语在线视频| 伊人久久精品亚洲午夜| av网站免费在线观看视频| 国产黄色免费在线视频| 久久久久视频综合| 国产综合精华液| 国产精品久久久久成人av| 亚洲精品中文字幕在线视频 | 日本黄色日本黄色录像| 熟女电影av网| 亚洲欧美成人综合另类久久久| 色哟哟·www| 91久久精品国产一区二区成人| 亚洲欧美中文字幕日韩二区| 国产亚洲一区二区精品| 人人妻人人看人人澡| 日本-黄色视频高清免费观看| 亚洲精品日韩av片在线观看| 久久综合国产亚洲精品| 久久精品国产亚洲av天美| 五月开心婷婷网| 日本猛色少妇xxxxx猛交久久| 女的被弄到高潮叫床怎么办| 欧美xxⅹ黑人| 久久综合国产亚洲精品| 精品亚洲成a人片在线观看| 伊人亚洲综合成人网| 国产高清国产精品国产三级| 超碰97精品在线观看| 亚洲中文av在线| 国产一区亚洲一区在线观看| 国产精品熟女久久久久浪| 午夜影院在线不卡| 国产探花极品一区二区| 视频区图区小说| 久久青草综合色| 蜜桃久久精品国产亚洲av| 大香蕉久久网| 亚洲精品,欧美精品| 视频中文字幕在线观看| 纯流量卡能插随身wifi吗| 亚洲欧美一区二区三区黑人 | 国产免费一区二区三区四区乱码| 国产成人午夜福利电影在线观看| 搡女人真爽免费视频火全软件| 水蜜桃什么品种好| 麻豆成人av视频| 一级毛片 在线播放| 免费在线观看成人毛片| 亚洲无线观看免费| 国产熟女午夜一区二区三区 | 97超视频在线观看视频| 欧美日韩国产mv在线观看视频| 极品少妇高潮喷水抽搐| 精品少妇黑人巨大在线播放| 久久久久精品久久久久真实原创| 色婷婷久久久亚洲欧美| 国产又色又爽无遮挡免| 亚洲美女视频黄频| 久久精品夜色国产| 久久久久国产网址| 黄片无遮挡物在线观看| 国产免费视频播放在线视频| 亚洲无线观看免费| 日韩熟女老妇一区二区性免费视频| 亚洲精品,欧美精品| 国产成人a∨麻豆精品| 中文字幕久久专区| 街头女战士在线观看网站| 大片免费播放器 马上看| 亚洲国产精品国产精品| 性高湖久久久久久久久免费观看| 菩萨蛮人人尽说江南好唐韦庄| 日本欧美国产在线视频| a级一级毛片免费在线观看| av播播在线观看一区| 国产免费一区二区三区四区乱码| 久久精品熟女亚洲av麻豆精品| 制服丝袜香蕉在线| √禁漫天堂资源中文www| 国产精品熟女久久久久浪| 80岁老熟妇乱子伦牲交| 亚洲精品乱码久久久久久按摩| 久久人人爽av亚洲精品天堂| 又爽又黄a免费视频| 日韩av在线免费看完整版不卡| 亚洲精品久久午夜乱码| 18禁动态无遮挡网站| 永久网站在线| 人人澡人人妻人| 久久久久久人妻| 亚洲精品日韩在线中文字幕| 91在线精品国自产拍蜜月| 蜜桃久久精品国产亚洲av| 黄色一级大片看看| 美女国产视频在线观看| 夜夜骑夜夜射夜夜干| 日韩视频在线欧美| 男女边摸边吃奶| 偷拍熟女少妇极品色| 晚上一个人看的免费电影| 亚洲经典国产精华液单| 亚洲av男天堂| 久久久久网色| 少妇的逼水好多| 日本av免费视频播放| 国精品久久久久久国模美| 亚洲国产精品国产精品| 久久人妻熟女aⅴ| 亚洲精品国产色婷婷电影| 妹子高潮喷水视频| 你懂的网址亚洲精品在线观看| 97在线视频观看| 中文字幕亚洲精品专区| 高清午夜精品一区二区三区| 男女国产视频网站| 亚洲不卡免费看| 成人毛片60女人毛片免费| 亚洲精品国产色婷婷电影| 如何舔出高潮| 插逼视频在线观看| 国产精品欧美亚洲77777| 久久人人爽人人片av| 亚洲无线观看免费| 亚洲欧洲国产日韩| av不卡在线播放| kizo精华| 久久精品久久久久久噜噜老黄| 婷婷色综合www| 久久久欧美国产精品| 3wmmmm亚洲av在线观看| 夜夜骑夜夜射夜夜干| 赤兔流量卡办理| 国产一区二区在线观看日韩| 春色校园在线视频观看| 欧美变态另类bdsm刘玥| 国产精品无大码| 一级毛片我不卡| av视频免费观看在线观看| 婷婷色综合大香蕉| 我要看黄色一级片免费的| 少妇被粗大的猛进出69影院 | 91久久精品电影网| 久久久精品94久久精品| 亚洲精品日本国产第一区| 丁香六月天网| 午夜日本视频在线| av在线播放精品| 一个人免费看片子| 免费高清在线观看视频在线观看| 91精品国产国语对白视频| 国产色婷婷99| 如日韩欧美国产精品一区二区三区 | 精品人妻一区二区三区麻豆| 久久久久久久久久久久大奶| 嫩草影院新地址| 国产高清三级在线| 国产精品国产三级专区第一集| 国精品久久久久久国模美| 少妇精品久久久久久久| 日日啪夜夜撸| 亚洲高清免费不卡视频| 久久婷婷青草| 少妇丰满av| 下体分泌物呈黄色| 国产欧美亚洲国产| 欧美精品国产亚洲| 国产无遮挡羞羞视频在线观看| 波野结衣二区三区在线| 成人18禁高潮啪啪吃奶动态图 | 成人影院久久| 蜜桃久久精品国产亚洲av| 高清黄色对白视频在线免费看 | 日韩视频在线欧美| 自拍偷自拍亚洲精品老妇| 精品少妇内射三级| 在线观看免费视频网站a站| 深夜a级毛片| 一本久久精品| 午夜免费男女啪啪视频观看| 日韩中字成人| 色婷婷久久久亚洲欧美| 午夜福利,免费看| 亚洲国产av新网站| 欧美精品亚洲一区二区| 黄色日韩在线| av在线播放精品| av视频免费观看在线观看| 国产成人精品婷婷| 亚洲三级黄色毛片| 国产国拍精品亚洲av在线观看| 国产日韩欧美在线精品| av免费在线看不卡| 91精品国产九色| 蜜臀久久99精品久久宅男| 亚洲国产色片| 在线精品无人区一区二区三| 一级毛片电影观看| 国产色婷婷99| 美女国产视频在线观看| 另类精品久久| 交换朋友夫妻互换小说| 插阴视频在线观看视频| 欧美日韩精品成人综合77777| 色视频在线一区二区三区| 国产精品国产av在线观看| 精品午夜福利在线看| 老熟女久久久| 国产极品粉嫩免费观看在线 | 久久国产亚洲av麻豆专区| 中国美白少妇内射xxxbb| 啦啦啦在线观看免费高清www| 久久精品国产鲁丝片午夜精品| 日本色播在线视频| 国产乱来视频区| 久久热精品热| 国产老妇伦熟女老妇高清| 高清午夜精品一区二区三区| 亚洲成人手机| 大香蕉久久网| 久久青草综合色| 久久久久久久久久成人| 国产精品一区二区在线不卡| 久久狼人影院| 极品人妻少妇av视频| 啦啦啦啦在线视频资源| 两个人的视频大全免费| 亚洲成人手机| 内射极品少妇av片p| 有码 亚洲区| 中文在线观看免费www的网站| 少妇人妻一区二区三区视频| 国产 精品1| 最近中文字幕2019免费版| 综合色丁香网| 少妇的逼水好多| 久久久久视频综合| 熟妇人妻不卡中文字幕| 日本与韩国留学比较| 欧美日韩av久久| 在线天堂最新版资源| av在线播放精品| kizo精华| av一本久久久久| 91精品一卡2卡3卡4卡| 欧美少妇被猛烈插入视频| 国产淫语在线视频| 国产精品人妻久久久影院| 十八禁高潮呻吟视频 | 亚洲av综合色区一区| 3wmmmm亚洲av在线观看| 免费看日本二区| 男人舔奶头视频| 男人爽女人下面视频在线观看| 黄色一级大片看看| 久久久久久久精品精品| 高清在线视频一区二区三区| 日韩制服骚丝袜av| 成人二区视频| 观看av在线不卡| 91精品国产九色| 久久久久久久久大av| 精品一区二区三卡| 精品午夜福利在线看| 精品久久久噜噜| 亚洲av福利一区| 在线观看一区二区三区激情| 国产精品人妻久久久影院| 一级爰片在线观看| 亚洲精品色激情综合| 国产精品人妻久久久久久| 18禁在线播放成人免费| 伦精品一区二区三区| 亚洲,一卡二卡三卡| 亚洲久久久国产精品| av国产精品久久久久影院| 看免费成人av毛片| 精品卡一卡二卡四卡免费| 欧美区成人在线视频| 中文字幕精品免费在线观看视频 | 亚洲av综合色区一区| 国产黄频视频在线观看| 国产精品久久久久久精品古装| 亚洲美女视频黄频| 夜夜爽夜夜爽视频| 男女边吃奶边做爰视频| 国产av一区二区精品久久| 91久久精品国产一区二区三区| 国产午夜精品一二区理论片| av又黄又爽大尺度在线免费看| 在线播放无遮挡| av在线播放精品| 九草在线视频观看| 亚洲国产最新在线播放| 精品国产露脸久久av麻豆| 老司机影院成人| 久久久午夜欧美精品| 高清毛片免费看| 国产欧美另类精品又又久久亚洲欧美| 夜夜爽夜夜爽视频| 青春草亚洲视频在线观看| 欧美+日韩+精品| 久久韩国三级中文字幕| 欧美日韩一区二区视频在线观看视频在线| 免费黄色在线免费观看| 亚洲天堂av无毛| 国产深夜福利视频在线观看| 一本—道久久a久久精品蜜桃钙片| 国产色婷婷99| 噜噜噜噜噜久久久久久91| 免费大片黄手机在线观看| 99久久综合免费| 永久网站在线| 国产精品嫩草影院av在线观看| 日韩人妻高清精品专区| 美女视频免费永久观看网站| 亚洲伊人久久精品综合| 男人添女人高潮全过程视频| 伦理电影免费视频| 伦理电影大哥的女人| 色视频在线一区二区三区| 少妇高潮的动态图| 精品少妇内射三级| 亚洲欧美精品自产自拍| 久久精品夜色国产| 美女主播在线视频| 在线观看人妻少妇| 亚洲欧洲精品一区二区精品久久久 | 一区在线观看完整版| 国产精品伦人一区二区| 中国三级夫妇交换| 国产伦理片在线播放av一区| 最近中文字幕高清免费大全6| 看十八女毛片水多多多| 男人和女人高潮做爰伦理| 成人毛片a级毛片在线播放| 韩国av在线不卡| 免费播放大片免费观看视频在线观看| 乱码一卡2卡4卡精品| 成人特级av手机在线观看| 成人亚洲精品一区在线观看| 一本一本综合久久| 只有这里有精品99| 亚洲精品一区蜜桃| 一边亲一边摸免费视频| 日日摸夜夜添夜夜添av毛片| 9色porny在线观看| 成人免费观看视频高清| 国产精品99久久99久久久不卡 | 黑人巨大精品欧美一区二区蜜桃 | 超碰97精品在线观看| 日韩制服骚丝袜av|