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

    Molecular detection of microbial communities associated with Microcystis vs Synechococcus dominated waters in Tianjin, China*

    2018-08-02 02:50:42HUODa霍達(dá)CHENYouxin陳友信LIUPing劉萍LIYiming李一鳴QIAOZhiyi喬之怡LIRenhui李仁輝
    Journal of Oceanology and Limnology 2018年4期
    關(guān)鍵詞:劉萍

    HUO Da (霍達(dá)) , CHEN Youxin (陳友信) , , LIU Ping (劉萍) , LI Yiming (李一鳴) QIAO Zhiyi (喬之怡) , LI Renhui (李仁輝) ,

    1 Department of Fisheries Sciences, Tianjin Agricultural University, National Demonstration Center for Experimental Aquaecology and Aquaculture Education, Tianjin 300384, China

    2 Key Laboratory of Algal Biology, State Key Laboratory of Freshwater Ecology and Biotechnology of China, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China

    3 Department of Basic Sciences, Tianjin Agricultural University, Tianjin 300384, China

    4 University of Chinese Academy of Sciences, Beijing 100049, China

    Abstract Harmful cyanobacterial blooms cause many ecological disasters worldwide. During the development of cyanobacterial blooms, the diversity and domination of cyanobacterial taxa are of a particular concern. In this study, the microbial community structure within a water system, such as in Yuqiao Reservoir and Haihe River in Tianjin City, China, was compared by using next-generation sequencing. A total of 5 001 operational taxonomic units were obtained and clustered from filtered 16S rDNA V3–V4 region sequences.The cyanobacterial and microbial structures greatly differed in these two water areas. Microcystis was dominant in Yuqiao, whereas Synechococcus was dominant in Haihe. Proteobacteria species were dominant among all detected samples. The relative abundances of Bacteroidetes and Planctomycetes were higher in Yuqiao Reservoir than in Haihe River, whereas Firmicutes and Verrucomicrobia were relatively abundant in Haihe River. Further analyses indicated that the domination of both cyanobacteria was strongly related to several environmental factors, such as total nitrogen, total phosphorus, and dissolved oxygen, re flecting the role of trophic states in shaping the dominance of cyanobacterial taxa. The present study provided the example for Microcystis and Synechococcus dominance along a cyanobacterial bloom in north China. Applying highthroughput sequencing could offer a wide field of vision in analyzing microbial community structures.

    Keyword: community structure; cyanobacteria; Yuqiao Reservoir; Haihe River; dominating patterns

    1 INTRODUCTION

    Massive proliferation of harmful cyanobacteria has become a growing problem in many freshwater ecosystems worldwide (Harke et al., 2016).Cyanobacteria blooms cause many environmental problems due to the increase in eutrophication level.Various genera/species of bloom-forming cyanobacteria produce toxins, further increasing the risk for aquatic organisms and humans (Fleming et al., 2002; Lone et al., 2015). Many studies have focused on the genetic diversity of cyanobacterial populations during blooms to determine the factors regulating the formation of water blooms (Zhu et al., 2012; Li et al., 2016; Kong et al., 2017). With the development of next-generation DNA sequencing techniques, several studies have recently investigated the microbial community structure during water blooms (Chen et al., 2011;Bertos-Fortis et al., 2016; Parulekar et al., 2017).China suffers from the heaviest cyanobacterial blooms, and considerable research focused on the water blooms in southern China (Zhang et al., 2016;Liu et al., 2017). However, rising temperature and carbon dioxide concentration in global and regional scales have aggravated the threat from cyanobacterial blooms (Paerl and Huisman, 2008; O’Neil et al.,2012; Verspagen et al., 2014; Visser et al., 2016), and such a consequence will lead to the extension of the water bloom range to the large area of northern China.

    Complicated relationships exist between aquatic bacteria and cyanobacterial species during cyanobacterial blooms in an aquatic ecosystem. From a general view, cyanobacteria offer a living habitat for heterotrophic bacteria; bacterial supersession could produce the source of nitrogen, phosphorus, sulfur,vitamin and carbon that cyanobacteria could utilize(Cai et al., 2013; Lee et al., 2014; Klawonn et al.,2015; Stuart et al., 2016; Xie et al., 2016). Some epiphytic bacteria are involved in complicated cyanobacterial life activities, which would contribute to cyanobacterial nutrition absorption, colony formation, and toxin degradation (Brunberg, 1999;Maruyama et al., 2003; Shi et al., 2010, 2012; Shen et al., 2011). Different cyanobacterial taxa with different cellular/colonial surface structures and physiobiochemical features will form differently attached microbial communities (Parveen et al.,2013). The most frequently occurring blooms are dominated by Microcystis species (Wirsing et al.,1998; Jav?rek et al., 2015), and various studies have focused on the microbial community structures in Microcystis dominant blooms. Nevertheless,super ficially Microcystis blooms are composed of a large number of picocyanobacterial cells, which are difficult to detect by traditional microscopic observation. Ye et al. (2011) used DGGE and realtime PCR methods and found that the two most dominant cyanobacteria in Taihu Lake are Microcystis and Synechococcus, which both have signi ficant succession during the development of water blooms.16S rDNA sequencing revealed that picocyanobacteria(mainly Synechococcus-like) are also important cyanobacterial members in Chaohu Lake (Cai and Kong, 2013). Synechococcus was found to be dominant in the eutrophic Dongqian Lake by using the targeted deep sequencing of cyanobacteria (Jiang et al., 2017). Therefore, picocyanobacterial abundance in cyanobacterial blooms might be neglected or largely underestimated when molecular detecting approaches are not used. Factors that drive the dominance and succession of different cyanobacteria during a cyanobacterial bloom or different blooms are still unknown; however, the characterization of different cyanobacterial blooms using traditional microscopic or modern molecular methods is an initial step for elucidating these driving forces.Previous studies related to picocyanobacterial dominance in cyanobacterial blooms were mainly recorded in southern and eastern China.

    In this study, we selected the Yuqiao-Haihe water system, a reservoir and river in Tianjin City, North China. This water system has annually suffered from huge cyanobacteria blooms. We used high-throughput sequencing of 16S rDNA to investigate the aquatic bacterial community structure along the water course from Yuqiao Reservoir to Haihe River during a bloom in 2016. The aquatic microbial community structures and environmental conditions in the reservoir and river systems were compared. We aimed to determine the controlling factors with regard to the distribution of Microcystis and Synechococcus-like picocyanobacteria.

    2 MATERIAL AND METHOD

    2.1 Study area and sampling strategy

    Yuqiao Reservoir and Haihe River are both located in Tianjin City, North China. Yuqiao Reservoir is at the upper part of Haihe River flowing through Tianjin City. The reservoir has a surface area of 86.8 km2with a storage capacity of 1.56×109m3, average depth of 4.74 m, and a length of 30 km. Yuqiao Reservoir is a drinking water source of Tianjin City. Haihe River was also suffering from cyanobacterial blooms due to the flowing water from Yuqiao, and the relative position between Yuqiao and Haihe is shown in Fig.1.

    Sampling was conducted on September 14, 2016,at 14 sites in the Yuqiao and Haihe area of Tianjin.Samples were collected during an active bloom. The coordinates of the sampling sites are shown in supplementary material S1. In brief, 2 L of water sample in each site was collected at 0.5 m below the surface. All water samples were subsequently divided into three subsamples: for the analysis of water quality parameters, for the filtration to high-throughput sequencing analysis, and for the counting of Microcystis cells (data shown in supplementary material S2). 500 mL of water samples was filtered through a 0.22-μm filter (Millipore, USA).

    2.2 Analysis of physicochemical parameters

    Fig.1 Schematic diagram of sampling sites

    Water temperature (WT), pH, dissolved oxygen(DO), and salinity (SAL) were measured using YSI multiparameter detector (YSI, USA). Total nitrogen(TN), total phosphorus (TP), chemical oxygen demand (COD), ammonia nitrogen (NH4+), and nitrite(NO2ˉ) were measured according to the standard methods (Wu et al., 2006).

    2.3 DNA extraction and library construction

    DNA was extracted from each filtered membrane by using the PowerWater DNA isolation kit following the manufacturer’s manual (Mo Bio Laboratories,Inc., Carlsbad, CA). The concentration of DNA was evaluated using a Qubit 2.0 fluorometer (Life Technologies Japan Ltd., Tokyo, Japan). For the ampli fication of the V3–V4 domain of bacterial 16S rDNA, primers 341F (5′-CCCTACACGACGCTCTTCCGATCTG (barcode) CCTACGGGNGGCWGCAG-3′) and 805R (5′-GACTGGAGTTCCTTGGCACCCGAGAATTCCAGACTACHVGGGTATCTAATCC-3′) (Herlemann et al., 2011) were used, both of which contained an Illumina adapter region for sequencing on the Illumina MiSeqTMplatform. In the forward primer, a 6-bp barcode was used to enable sample multiplexing. PCR was performed twice. The first-round PCR contained 10–20 ng DNA template,1× Taq Master Mix (TaKaRa, Tokyo, Japan), and 2 pmole of both bar primer and reverse primer.Reaction conditions consisted of an initial 94°C for 3 min, followed by 5 cycles of 94°C for 30 s, 45°C for 20 s, and 65°C for 30 s; 20 cycles of 94°C for 20 s,55°C for 20 s, 72°C for 30 s, and a final extension of 72°C for 5 min. Illumina Bridge PCR primer, which contained 20 ng DNA template, 1× Taq Master Mix(TaKaRa, Tokyo, Japan), and 2 pmol of each primer were added during the second-round PCR. Reaction conditions consisted of an initial 95°C for 30 s,followed by 5 cycles of 95°C for 15 s, 55°C for 15 s,and 72°C for 30 s with the final extension of 72°C for 5 min. All samples were analyzed on the T100TM Thermal Cycler (Bio-Rad, NY, USA). Amplicons were puri fied using the QIAquick 96 kit (Qiagen) and quanti fied using the Qubit2.0 fluorometer (Life Technologies Japan Ltd., Tokyo, Japan). Gel puri fication was performed using QIAquick gel extraction kit (Qiagen). All operations were performed according to the respective manufacturer’s instructions. Finally, 10 ng of each PCR product was mixed for sequencing, and the final concentration of PCR reaction was 20 pmole.

    2.4 Data analyses

    Raw Illumina FASTQ files were demultiplexed,quality filtered, and analyzed using Cutadapt (Martin,2011), PEAR (Zhang et al., 2014), Prinseq (Schmieder and Edwards, 2011), QIIME (Caporaso et al., 2010),and mothur (Schloss et al., 2009). Operational taxonomic units (OTUs) were assigned using the QIIME implementation of USEARCH (Edgar, 2010)with a threshold of 97% pairwise identity, andrepresentative sequences from each OTU were selected for taxonomic assignment. OTUs were classi fied taxonomically by using Mothur method,and singleton OTU tag was removed prior to analysis.Multiple indices (observed OTUs, Chao 1 richness,Shannon’s index, and Simpson’s index) were calculated using the QIIME package to represent diversity.

    Table 1 Environmental parameters of the samples

    2.5 Statistic

    All statistical analyses were performed in the R environment using VEGAN (Dixon, 2003), grid(Wegener et al., 2009), and ggplot2 (Wickham, 2016)packages. Results were presented as mean values ±standard deviation. PCA was calculated by euclidian distance. A detrended correspondence analysis was performed before further analyses. The unimodally distributed model was used for canonical correlation analysis (CCA). For the principal component analysis(PCA) and CCA, logarithmic results were used;logarithmic result=log(otu#+1). Signi ficance test(Student’s t-test) was completed using SPSS (version 22.0).

    3 RESULT

    3.1 Environmental parameters

    The characteristics of the surface water samples from Yuqiao Reservoir and Haihe River are summarized in Table 1. The weather during the sampling process was sunny. WT was not signi ficantly different and was in the range of 20.08°C to 25.01°C.DO in Haihe River was higher than that in Yuqiao Reservoir. In most of the sampling sites, TN in Yuqiao Reservoir was higher than that in Haihe River.However, center of the reservoir, had TN of 0.80 mg/L,which lower than in other sampling sites in Yuqiao.Nitrogen and phosphorus were spatially in fluctuating state. COD in Yuqiao Reservoir (8.17±0.69 mg/L)was signi ficantly higher than that in Haihe River(2.60±1.18 mg/L). The SAL of water in Yuqiao Reservoir was more stable than that in Haihe River.

    3.2 Amplicon of 16S rRNA gene sequencing

    After quality control, the OTU singleton was removed. A total of 5 001 OTUs were obtained and clustered from 46 9708 filtered 16S V3–V4 region sequences. The largest number of reads was obtained from sample YQ_SF, whereas the highest value of OTUs was observed at HH_AMB. Compared with Yuqiao Reservoir, Haihe River showed high diversity in all kinds of diversity indices we selected (OTU number, Chao 1, good coverage, Shannon, Simpson)during our research period.

    All samples were analyzed using a subsample for alpha diversity analysis. The minimum number of OTU counts (27 525) was assigned as e value. As summarized in Table 2, most sampling sites in Haihe River had greater richness than those in Yuqiao Reservoir. With the increasing sequence depths, the rarefaction curves had become flat. No higher diversity value was observed, and our subsample amount was sufficient to re flect the diversity for each sample. OTU number curves, Chao 1 index, andShannon index showed high diversity in Haihe River.

    Table 2 Multiple diversity indices of the 14 samples

    Fig.2 Rarefaction curves for a dissimilarity of 3% from the 14 samples

    3.3 Microbial community composition

    The script assign_taxonomy.py in QIIME software was used to assign the taxonomy information into different phylogenetic bacterial taxa. This method was called “mothur”. The relative abundance of bacterial community at phylum level is shown in Fig.3. Proteobacteria were the most abundant phylum in all experimental sites accounting for 53.3% of total bacterial sequences, followed by Bacteroidetes(10.4%), Firmicutes (7.1%), Actinobacteria (7.1%),Verrucomicrobia (6.9%), and cyanobacteria (5.7%).The relative abundances of both Bacteroidetes and Planctomycetes in Yuqiao Reservoir were higher than those in Haihe River, However, the relative abundances of Firmicutes and Verrucomicrobia in Haihe River were higher than those in Yuqiao Reservoir. The cyanobacterial genera with relative abundance greater than 0.001 were analyzed to further understand the differences in cyanobacterial genera between the two water areas (Fig.4). Microcystis had the biggest abundance in Yuqiao Reservoir; however,Synechococcus was dominant among all cyanobacterial genera in Haihe River. By analyzing the relative abundances of cyanobacterial taxa, two signi ficantly different cyanobacterial communities were found in the two water areas ( P <0.01), with Microcystis taking an absolutely large part in cyanobacterial organisms in Yuqiao and Synechococcus dominating but other cyanobacterial taxa co-existing in Haihe River (Fig.4)

    3.4 Principal component analysis for relative abundance of OTUs in two water area

    PCA analysis revealed that all samples in Yuqiao Reservoir were clearly clustered in a group, indicating that the samples in Yuqiao Reservoir had high similarity in community characters. By contrast, the samples in Haihe River were separated into different diverse zones. Results indicated that the microbial community composition in Yuqiao Reservoir was more consistent than that in Haihe River (Fig.5)

    3.5 Correlation between microbial community composition with environmental variables

    Fig.4 Cyanobacterial communities and eukaryotic algae with chloroplast based on the abundance of the samples

    Fig.5 Principal component analysis for relative abundance of OTUs in two water areas, and all OTU numbers were added one and then performed a logarithmic transformation

    Fig.6 CCA analysis for the relationship between environmental factors and the relative abundance of bacteria

    We used CCA to detect the correlations between OTU abundance of 14 sampling sites with 9 selected environmental variables (TN, TP, ratio of TN and TP,ammonia, nitrite, WT, DO, SAL, and pH). As presented in Fig.6, the water samples in Yuqiao Reservoir and Haihe River were clearly divided.Among all the environmental factors, the investigated communities were signi ficantly distributed against the gradient of DO, SAL, TN, and TP. Seems to be mainly affected by these factors. In Yuqiao Reservoir,TP and TN were positive factors to most sample sites,and DO, SAL, NH4+, and pH were negative factors.YQ_CR had weaker relation to these factors than other Yuqiao sites. Interestingly, to the species in Haihe River, the distribution in CCA plot was divided in to three different groups. They have both positive correlation to DO, SAL, and NH4+but at different levels. The HH_TTE, HH_AMB, HH_TJS, and HH_BYP were weakly affected by the listed environmental factors. HH_SZLB showed close relationship with DO, SAL, and pH. HH_JBB and HH_ZYHSR were affected by the different formations of nitrogen.

    4 DISCUSSION

    To elucidate the mechanism by which cyanobacterial blooms form and factors control their formation,various studies were conducted on the relationship between cyanobacterial biomass, cyanotoxin production, nutrients, water residence time, and circulation (Vincent et al., 1993; Palíková et al., 2007;Harvey, 2012; Romo et al., 2013; Sinang et al., 2013).Speci fically, for Microcystis-dominated blooms, an increasing number of studies were focused on spatiotemporal dynamics of Microcystis populations based on molecular genetic markers (Briand et al.,2009). These studies revealed that many Microcystis genotypes coexists in a certain bloom, and a succession of dominant genotypes occurs along the development of the bloom (Bozarth et al., 2010). Further studies based on molecular biological methods showed that picocyanobacterial Synechococcus group could largely co-occur and even dominate in super ficially Microcystis-based blooms in southern Chinese lakes(Li et al., 2014). In this study, we provided the occurrence of both Microcystis- and Synechococcusdominated blooms along the reservoir-river water system in the north of China, the latter one representing the case with higher Synechococcus abundance than the Microcystis one. Studies are rare providing information on cyanobacterial blooms in a large water body in northern China. Our study is the first study aiming at the different cyanobacterial domination patterns during a bloom and their relationship with aquatic bacterial communities in northern China. Such an integrated water system as Yuqiao Reservoir and Haihe River offer a natural experiment model to understand microbial community structures under different cyanobacterial domination patterns.

    In the present study, Illumina-based 16S rRNA gene sequencing and subsequent analyses revealed the signi ficant differences in cyanobacterial composition between Yuqiao Reservoir and Haihe Ricer ( P <0.01), with more Synechococcus abundance in Haihe River. Therefore, Yuqiao and Haihe were spatially divided into Microcystis- and Synechococcusdominated patterns, respectively, during the bloom.The higher eutrophication in Yuqiao Reservoir was suggested to cause the annual the Microcystisdominated blooms in recent years, and the nutrient levels as shown in Table 1 also indicated the trophic state in the reservoir. The dominance of Synechococcus posed a very interesting finding. The dominance of Synechococcus has been well known in both marine and freshwater environments (Li, 1998; Takasu et al.,2015). In marine systems, Synechococcus favors mesotrophic conditions (Grob et al., 2007; Post et al.,2011). However, the occurrence of Synechococcus in freshwater was found in all trophic states, and Synechococcus even dominated in oligotrophic lakes in Germany (Ruber et al., 2016) and eutrophic lakes in southern China (Cai et al., 2013). Compared with Yuqiao Reservoir, a lower trophic state in Haihe River was implicated as the main potential factor to cause the dominance of Synechococcus, in agreement with the finding of Ruber et al. (2016).

    Further analyses showed that microbial community structures in Yuqiao Reservoir and Haihe River largely differed, represented as higher Bacteroidetes and Planctomycetes in Yuqiao Reservoir and higher Firmicutes and Verrucomicrobia in Haihe River.Bacteroidetes was the most representative group in clone libraries of four bloom-forming lakes (Eiler and Bertilsson, 2004). Hu et al. (2013) observed that Bacteroidetes had an obvious correlation with the dominant Microcystis-based bloom in Lake Erhai. Shi et al. (2010) found that, by the semi-nested PCR and DGGE methods, Serratia (Proteobacteria),Pseudomonas (Proteobacteria), and Flexibacteraceae(Bacteroidetes) were the dominant bacterial groups in the M. flos- aquae, M. wesenbergii, and M. aeruginosa blooms in Lake Taihu, respectively, consistent with our analyses on the relative abundance of the microbial community in Yuqiao samples. Thus, these findings and the result obtained in this study may implicate the role of Bacteroidetes in the formation of Microcystisdominated blooms. With regard to the association of Planctomycetes with Microcystis blooms, Cai et al.(2013) indicated that Planctomycetes tightly attached to carpet-like mucilaginous cyanobacterial blooms might play a key role in the carbon and nutrient cycling. Further, they speculated that the ecological role of Planctomycetes from the Lake Taihu blooms was linked to the degradation of sulfated polysaccharides produced by cyanobacteria. This finding in Lake Taihu is meaningful to the case in Yuqiao Reservoir, and the higher abundance of Planctomycetes in Yuqiao Reservoir may signify the degradation of the Microcystis bloom. Soon the same bloom flowing into Haihe River would likely be degraded by the combined effect of this microbial community attached to the Microcystis colonies in Yuqiao Reservoir and environmental conditions in Haihe River. Consequently, Synechococcus dominance would occur in Haihe River.

    However, very few studies have been previously conducted to reveal the association of bacterial community with dominance of Synechococcus-like picocyanobacteria in freshwaters, and more samples from Synechococcus-dominated freshwater environments are needed to detect their associated microbial community structures in future studies.

    An important ecological method for understanding microbial community structures is to determine the possible environmental factors and their contribution to the variation of bacterial communities. Several environmental factors, including WT (Moll et al.,1999), DO (Stadler and Love, 2016), nitrogen (Robles et al., 2015), and phosphorus (Carvalho et al., 2007),were reported to be effective factors on the distribution of aquatic bacterial communities. The CCA analyses for bacterial communities in the two waters revealed different patterns because they were signi ficantly related to different environmental factors. Among them, DO, TN, and TP were strongly related to the microbial community structure in the Yuqiao and Haihe.

    In general, the interesting results obtained in the present study were largely due to the fast development of molecular biological technologies, especially for the availability of high-throughput sequencing method. By using the advanced approach,Synechococcus detection was done using a large amount, which is unlikely to be done under the traditional microscopic examination. It is expected that more studies in future will reveal more cases with picocyanobacterial dominance in all trophic states in freshwaters and possible linkage between this dominance and microbial community structures.

    5 CONCLUSION

    By using next-generation sequencing, the Microcysti s- and Synechococcus-dominated patterns were found in two water areas within a water system along a cyanobacterial bloom in Tianjin City. Both dominations were strongly related to several environmental factors, such as TN, TP, and DO,re flecting the role of trophic states. The Microcystisdominated pattern was shown to have a relatively high abundance of Bacteroidetes and Planctomycetes,whereas the Synechococcus dominance was possibly related to the high abundance of Firmicutes and Verrucomicrobia. The application of high-throughput sequencing could offer a broader field of vision in analyzing microbial community structures.

    6 ACKNOWLEDGMENT

    We would like to thank ZHENG Tao, JI Dongdong,QIN Yiming and LI Na at Tianjin Agricultural University for their work to sample collecting.Especially thank to Ms. CAO Qi for her support to data arrangement.

    猜你喜歡
    劉萍
    刻薄
    愛著你的愛
    劉萍作品
    杭州藝尚小鎮(zhèn),數(shù)字時尚向未來!
    單堯堯、胡豐凡、劉萍、葉小輝作品
    做大自然的“翻譯官”
    ——淺談寫生的意義
    丈夫保險不保險
    輪回
    輪回
    上海故事(2014年12期)2014-12-02 16:12:20
    檔案記載:為抗日英勇就義的徐水縣長劉萍
    檔案天地(2014年6期)2014-06-13 17:35:40
    久久久久亚洲av毛片大全| 欧美成人a在线观看| 国模一区二区三区四区视频| 国产精品伦人一区二区| 91狼人影院| 精品人妻偷拍中文字幕| 九色成人免费人妻av| www.www免费av| 免费黄网站久久成人精品 | 男女之事视频高清在线观看| 两人在一起打扑克的视频| 国产精品野战在线观看| 成人av一区二区三区在线看| 丰满人妻一区二区三区视频av| 一个人观看的视频www高清免费观看| 桃色一区二区三区在线观看| 麻豆一二三区av精品| 我的老师免费观看完整版| 我的女老师完整版在线观看| 久久久久国内视频| 级片在线观看| 熟女电影av网| 欧美黄色片欧美黄色片| 午夜福利视频1000在线观看| 十八禁国产超污无遮挡网站| 精品国产亚洲在线| 深夜精品福利| 久久午夜亚洲精品久久| 国产高清三级在线| 9191精品国产免费久久| 国产成+人综合+亚洲专区| 亚洲经典国产精华液单 | 毛片一级片免费看久久久久 | 国产免费av片在线观看野外av| 黄片小视频在线播放| 欧美一区二区亚洲| 99国产综合亚洲精品| 国产一区二区亚洲精品在线观看| 最近视频中文字幕2019在线8| 亚洲aⅴ乱码一区二区在线播放| 欧美乱妇无乱码| 欧美丝袜亚洲另类 | 国产在视频线在精品| 亚洲av中文字字幕乱码综合| 亚洲人成网站在线播| 少妇人妻精品综合一区二区 | 国产精品一区二区免费欧美| 99热精品在线国产| 亚洲成人久久性| 久久人人爽人人爽人人片va | 国产av麻豆久久久久久久| 又黄又爽又刺激的免费视频.| 好看av亚洲va欧美ⅴa在| 九九在线视频观看精品| av天堂在线播放| 欧美午夜高清在线| 精品一区二区免费观看| 婷婷精品国产亚洲av| 亚洲经典国产精华液单 | 亚州av有码| 露出奶头的视频| 国模一区二区三区四区视频| 国产中年淑女户外野战色| 国产 一区 欧美 日韩| 91九色精品人成在线观看| 熟女电影av网| 精品久久久久久,| www.999成人在线观看| av视频在线观看入口| 精品人妻熟女av久视频| 直男gayav资源| 国产精品不卡视频一区二区 | 成年人黄色毛片网站| 天堂网av新在线| 一进一出抽搐gif免费好疼| 国产三级黄色录像| 欧美色欧美亚洲另类二区| 日韩欧美在线乱码| av国产免费在线观看| 欧美zozozo另类| 亚洲avbb在线观看| 中文亚洲av片在线观看爽| 欧美黑人欧美精品刺激| 狂野欧美白嫩少妇大欣赏| 午夜日韩欧美国产| 精华霜和精华液先用哪个| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产亚洲av天美| 国产老妇女一区| 亚洲,欧美精品.| 日韩亚洲欧美综合| 91在线观看av| 亚洲成人中文字幕在线播放| 亚洲欧美日韩卡通动漫| 久久国产精品人妻蜜桃| 久久热精品热| 性色av乱码一区二区三区2| 校园春色视频在线观看| 欧美区成人在线视频| 亚洲av免费在线观看| 99在线视频只有这里精品首页| 国产欧美日韩精品一区二区| 色在线成人网| 丝袜美腿在线中文| 成人美女网站在线观看视频| 婷婷精品国产亚洲av在线| 午夜亚洲福利在线播放| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区三区视频了| 国产一区二区激情短视频| 99久久精品国产亚洲精品| 欧美绝顶高潮抽搐喷水| 国产精品一及| 嫁个100分男人电影在线观看| 日韩高清综合在线| 成人国产一区最新在线观看| 99久久99久久久精品蜜桃| 午夜两性在线视频| 88av欧美| 亚洲精品久久国产高清桃花| 国产单亲对白刺激| 亚洲美女搞黄在线观看 | 国产精品久久久久久精品电影| 欧美性感艳星| av女优亚洲男人天堂| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美日韩东京热| 午夜福利在线在线| 波多野结衣高清作品| 最近最新免费中文字幕在线| 欧美黄色淫秽网站| 美女被艹到高潮喷水动态| 国产精品久久久久久久电影| 男女下面进入的视频免费午夜| 一级a爱片免费观看的视频| 欧美一区二区国产精品久久精品| 国产精品久久久久久人妻精品电影| 老司机福利观看| 小说图片视频综合网站| 90打野战视频偷拍视频| 欧洲精品卡2卡3卡4卡5卡区| 欧美一区二区精品小视频在线| 深爱激情五月婷婷| 国产精品伦人一区二区| 国产午夜精品论理片| 日韩成人在线观看一区二区三区| 91av网一区二区| 亚洲男人的天堂狠狠| 99热这里只有是精品在线观看 | 亚洲精品亚洲一区二区| 亚洲av美国av| 亚洲av日韩精品久久久久久密| 免费大片18禁| 日本一二三区视频观看| 欧洲精品卡2卡3卡4卡5卡区| 特级一级黄色大片| 亚洲精品在线观看二区| 99在线人妻在线中文字幕| 97热精品久久久久久| 午夜福利成人在线免费观看| 欧美日韩国产亚洲二区| 成人国产一区最新在线观看| 国产成人福利小说| 成年人黄色毛片网站| 国产又黄又爽又无遮挡在线| 午夜激情福利司机影院| 亚洲人与动物交配视频| 真人做人爱边吃奶动态| 国产高清视频在线播放一区| 天堂√8在线中文| 国产视频内射| 夜夜看夜夜爽夜夜摸| 亚洲av电影在线进入| 欧美日韩亚洲国产一区二区在线观看| 12—13女人毛片做爰片一| 少妇裸体淫交视频免费看高清| 国产精品永久免费网站| av天堂在线播放| 日本 欧美在线| 中文字幕久久专区| 久久精品久久久久久噜噜老黄 | 九色国产91popny在线| 国产精品影院久久| 十八禁人妻一区二区| 人人妻人人看人人澡| 中文字幕人成人乱码亚洲影| 国产人妻一区二区三区在| 国产精品乱码一区二三区的特点| 国产成人福利小说| 久久伊人香网站| 国产乱人视频| 日韩欧美精品免费久久 | 69av精品久久久久久| 国内精品久久久久精免费| 久久国产乱子伦精品免费另类| 久久久久九九精品影院| 老司机深夜福利视频在线观看| 精品乱码久久久久久99久播| 亚洲第一欧美日韩一区二区三区| 身体一侧抽搐| 欧美黑人巨大hd| 免费一级毛片在线播放高清视频| 男女做爰动态图高潮gif福利片| 能在线免费观看的黄片| 麻豆成人av在线观看| 淫秽高清视频在线观看| 亚洲人成网站高清观看| 男女下面进入的视频免费午夜| 简卡轻食公司| 小说图片视频综合网站| 亚洲狠狠婷婷综合久久图片| 午夜久久久久精精品| 国产av在哪里看| a在线观看视频网站| 日日摸夜夜添夜夜添av毛片 | 毛片女人毛片| 日韩国内少妇激情av| 亚州av有码| 免费看日本二区| 国产精品一区二区性色av| 日本a在线网址| 亚洲三级黄色毛片| 精品一区二区免费观看| 国产成人欧美在线观看| 在线观看av片永久免费下载| 嫁个100分男人电影在线观看| 色5月婷婷丁香| 日本精品一区二区三区蜜桃| 内地一区二区视频在线| 日本与韩国留学比较| 亚洲欧美清纯卡通| 欧美不卡视频在线免费观看| 欧美性猛交╳xxx乱大交人| 真实男女啪啪啪动态图| 91av网一区二区| 国产真实乱freesex| 少妇高潮的动态图| 一级a爱片免费观看的视频| 日韩亚洲欧美综合| 三级男女做爰猛烈吃奶摸视频| 国产国拍精品亚洲av在线观看| 18禁黄网站禁片午夜丰满| 五月玫瑰六月丁香| 国产黄色小视频在线观看| 美女 人体艺术 gogo| 精品人妻一区二区三区麻豆 | av天堂在线播放| 久久久精品欧美日韩精品| 国产精品亚洲av一区麻豆| 少妇熟女aⅴ在线视频| 少妇丰满av| 在线观看午夜福利视频| 久久人妻av系列| 国产免费一级a男人的天堂| 精品久久国产蜜桃| 中文字幕精品亚洲无线码一区| 国产精品不卡视频一区二区 | 日韩欧美在线二视频| 午夜福利在线在线| 精品人妻熟女av久视频| 999久久久精品免费观看国产| 久久久久久久精品吃奶| 欧美绝顶高潮抽搐喷水| 国产一区二区三区在线臀色熟女| 窝窝影院91人妻| 他把我摸到了高潮在线观看| 校园春色视频在线观看| 91av网一区二区| bbb黄色大片| 一区二区三区激情视频| 最近最新免费中文字幕在线| 欧美+日韩+精品| 国产三级黄色录像| 精品久久久久久久久久久久久| 亚洲国产精品sss在线观看| 熟女人妻精品中文字幕| 1000部很黄的大片| 亚洲成av人片免费观看| 欧美日韩综合久久久久久 | 亚洲五月婷婷丁香| 日本成人三级电影网站| 日本 欧美在线| 日日摸夜夜添夜夜添av毛片 | 日韩 亚洲 欧美在线| 全区人妻精品视频| 一边摸一边抽搐一进一小说| 伦理电影大哥的女人| 久久久精品欧美日韩精品| 看免费av毛片| 亚洲激情在线av| 色av中文字幕| 精品久久久久久久久亚洲 | 99久久99久久久精品蜜桃| 亚洲国产欧洲综合997久久,| 亚洲人成伊人成综合网2020| 夜夜躁狠狠躁天天躁| 99久国产av精品| 欧美日韩综合久久久久久 | 又爽又黄a免费视频| 黄色日韩在线| 嫩草影院精品99| 久久久久久久亚洲中文字幕 | 99热6这里只有精品| 在线观看66精品国产| 国产精品女同一区二区软件 | 91麻豆av在线| 午夜精品一区二区三区免费看| 亚洲成人久久性| 国产精品女同一区二区软件 | 怎么达到女性高潮| 99国产极品粉嫩在线观看| а√天堂www在线а√下载| 日韩欧美在线二视频| 悠悠久久av| 日本免费a在线| 一本精品99久久精品77| 日本免费一区二区三区高清不卡| 免费大片18禁| 免费黄网站久久成人精品 | 18禁黄网站禁片午夜丰满| 激情在线观看视频在线高清| 99在线人妻在线中文字幕| 成年女人毛片免费观看观看9| 亚洲电影在线观看av| 日日干狠狠操夜夜爽| 一本综合久久免费| 国产精品精品国产色婷婷| 麻豆一二三区av精品| 99在线人妻在线中文字幕| 99热这里只有是精品50| 一二三四社区在线视频社区8| 中文在线观看免费www的网站| 国产高清激情床上av| 成熟少妇高潮喷水视频| 91av网一区二区| 美女高潮的动态| 99国产极品粉嫩在线观看| 狂野欧美白嫩少妇大欣赏| 一区二区三区高清视频在线| 亚洲成人久久爱视频| 国产高潮美女av| 好男人在线观看高清免费视频| 亚洲 欧美 日韩 在线 免费| 欧美成狂野欧美在线观看| 成人特级黄色片久久久久久久| 一边摸一边抽搐一进一小说| 全区人妻精品视频| 欧美乱妇无乱码| 精品国产亚洲在线| 最近在线观看免费完整版| 不卡一级毛片| 欧美日韩福利视频一区二区| 日韩精品青青久久久久久| 国产成人aa在线观看| 国产精品综合久久久久久久免费| 欧美xxxx性猛交bbbb| 久久热精品热| aaaaa片日本免费| 亚洲精品在线观看二区| 欧美3d第一页| 中文字幕人成人乱码亚洲影| 日本三级黄在线观看| 亚洲美女视频黄频| 午夜久久久久精精品| 欧美午夜高清在线| 成人美女网站在线观看视频| 国产黄a三级三级三级人| 欧美潮喷喷水| or卡值多少钱| 欧美绝顶高潮抽搐喷水| 亚洲黑人精品在线| 精品不卡国产一区二区三区| 亚洲成av人片免费观看| 国内少妇人妻偷人精品xxx网站| 精品熟女少妇八av免费久了| 欧美又色又爽又黄视频| 国产成人欧美在线观看| av欧美777| 色尼玛亚洲综合影院| 丰满的人妻完整版| 精品久久久久久久久久免费视频| 国产成人福利小说| 国产精品不卡视频一区二区 | 三级男女做爰猛烈吃奶摸视频| 精品久久久久久久末码| 精品一区二区三区视频在线观看免费| 99精品久久久久人妻精品| 欧美一区二区精品小视频在线| 可以在线观看的亚洲视频| 嫩草影视91久久| 老司机午夜十八禁免费视频| 91av网一区二区| 校园春色视频在线观看| 91麻豆精品激情在线观看国产| 国产色爽女视频免费观看| 国产成年人精品一区二区| 日韩大尺度精品在线看网址| a级毛片a级免费在线| 久久久久免费精品人妻一区二区| av视频在线观看入口| 白带黄色成豆腐渣| 美女cb高潮喷水在线观看| av在线蜜桃| 一进一出抽搐动态| 久久国产乱子伦精品免费另类| 久久草成人影院| avwww免费| 国内久久婷婷六月综合欲色啪| 国产私拍福利视频在线观看| 啦啦啦韩国在线观看视频| 色综合站精品国产| 国产探花极品一区二区| 国产精品久久久久久人妻精品电影| av福利片在线观看| 1024手机看黄色片| 日韩人妻高清精品专区| 亚洲成a人片在线一区二区| 亚洲国产精品999在线| 校园春色视频在线观看| 日韩欧美三级三区| 真人做人爱边吃奶动态| 亚洲av一区综合| 日本熟妇午夜| 性色avwww在线观看| 亚洲一区二区三区色噜噜| 国产精品美女特级片免费视频播放器| 午夜福利18| 精品99又大又爽又粗少妇毛片 | 桃红色精品国产亚洲av| 久久热精品热| 亚洲精品456在线播放app | 午夜福利在线观看吧| 亚洲 欧美 日韩 在线 免费| 国产成人av教育| 国产午夜精品久久久久久一区二区三区 | 18禁黄网站禁片免费观看直播| 亚洲人成网站在线播| 日韩欧美一区二区三区在线观看| 脱女人内裤的视频| 免费看美女性在线毛片视频| 国产精品1区2区在线观看.| 青草久久国产| 午夜福利在线在线| 国产精品久久电影中文字幕| 天堂动漫精品| 日日摸夜夜添夜夜添av毛片 | 国产久久久一区二区三区| 在线免费观看的www视频| 国产av麻豆久久久久久久| 欧美性感艳星| 欧美最新免费一区二区三区 | 婷婷亚洲欧美| 亚洲成av人片免费观看| 国产精品影院久久| 亚洲久久久久久中文字幕| 美女黄网站色视频| 在线国产一区二区在线| 日韩av在线大香蕉| 精品一区二区三区人妻视频| 成人三级黄色视频| 亚洲不卡免费看| 男女床上黄色一级片免费看| 18+在线观看网站| 91在线精品国自产拍蜜月| 99视频精品全部免费 在线| 成人av在线播放网站| 亚洲美女黄片视频| 日韩欧美国产在线观看| a级毛片免费高清观看在线播放| 日本熟妇午夜| 日本一二三区视频观看| 99国产极品粉嫩在线观看| 特大巨黑吊av在线直播| 深夜精品福利| 91狼人影院| 国产精品久久视频播放| 日韩免费av在线播放| 少妇人妻精品综合一区二区 | 国产精品一区二区性色av| 国产精品自产拍在线观看55亚洲| 又黄又爽又刺激的免费视频.| 久久精品国产自在天天线| 亚洲一区高清亚洲精品| 亚洲国产高清在线一区二区三| 亚洲最大成人中文| 99在线人妻在线中文字幕| 直男gayav资源| 午夜福利欧美成人| 性欧美人与动物交配| 国产男靠女视频免费网站| 日韩有码中文字幕| 欧美一区二区精品小视频在线| 国产欧美日韩一区二区精品| 久久精品91蜜桃| 午夜福利视频1000在线观看| 小说图片视频综合网站| 性插视频无遮挡在线免费观看| 亚洲色图av天堂| 老鸭窝网址在线观看| 伦理电影大哥的女人| 免费人成视频x8x8入口观看| 日韩精品中文字幕看吧| 岛国在线免费视频观看| 国产毛片a区久久久久| 一区二区三区四区激情视频 | 亚洲自拍偷在线| 亚洲精品456在线播放app | 精品乱码久久久久久99久播| 欧美高清性xxxxhd video| 精华霜和精华液先用哪个| 无人区码免费观看不卡| 日韩欧美 国产精品| 国产美女午夜福利| 亚洲精品在线观看二区| 免费看日本二区| 久久精品夜夜夜夜夜久久蜜豆| 观看免费一级毛片| 女生性感内裤真人,穿戴方法视频| 真人做人爱边吃奶动态| 精品久久国产蜜桃| 18+在线观看网站| 制服丝袜大香蕉在线| 人妻制服诱惑在线中文字幕| 性色avwww在线观看| 三级国产精品欧美在线观看| 在线免费观看不下载黄p国产 | 欧美bdsm另类| 国产久久久一区二区三区| 国内揄拍国产精品人妻在线| 国产综合懂色| 亚洲熟妇中文字幕五十中出| 国产探花在线观看一区二区| 亚洲最大成人av| 99国产精品一区二区三区| 毛片一级片免费看久久久久 | 亚洲经典国产精华液单 | 亚洲经典国产精华液单 | 中文亚洲av片在线观看爽| 禁无遮挡网站| 欧美一区二区精品小视频在线| 国产亚洲欧美98| 亚洲片人在线观看| 日本一本二区三区精品| 色综合站精品国产| 熟妇人妻久久中文字幕3abv| 欧美成人a在线观看| 国产午夜福利久久久久久| 国产成人av教育| 18禁黄网站禁片免费观看直播| 又爽又黄a免费视频| 国产成人福利小说| 国产精品99久久久久久久久| 99久久无色码亚洲精品果冻| 精品一区二区三区人妻视频| 成年人黄色毛片网站| 每晚都被弄得嗷嗷叫到高潮| 国产真实乱freesex| 久久久久亚洲av毛片大全| 欧美黄色淫秽网站| www日本黄色视频网| 精品人妻一区二区三区麻豆 | 又爽又黄a免费视频| 国产精品一及| 久久久精品大字幕| 国产一区二区激情短视频| 久久久久久久久中文| 99国产综合亚洲精品| 国产亚洲精品综合一区在线观看| 综合色av麻豆| 国产午夜精品论理片| 欧美日韩乱码在线| 国产毛片a区久久久久| 俺也久久电影网| 中文亚洲av片在线观看爽| 成人无遮挡网站| 午夜福利在线观看免费完整高清在 | 嫩草影院新地址| av福利片在线观看| 欧美成狂野欧美在线观看| 国产爱豆传媒在线观看| 日本黄色片子视频| 草草在线视频免费看| 亚洲人成网站高清观看| 欧美日韩中文字幕国产精品一区二区三区| 色5月婷婷丁香| 久久亚洲真实| 国产一区二区三区视频了| 在线免费观看不下载黄p国产 | 九九久久精品国产亚洲av麻豆| 久久国产精品人妻蜜桃| 国产一区二区亚洲精品在线观看| 国产成+人综合+亚洲专区| 毛片一级片免费看久久久久 | 精品国产亚洲在线| 午夜福利成人在线免费观看| 欧美一区二区国产精品久久精品| 久久香蕉精品热| 婷婷丁香在线五月| 90打野战视频偷拍视频| 动漫黄色视频在线观看| 亚洲av电影不卡..在线观看| 舔av片在线| 麻豆久久精品国产亚洲av| 久久精品国产亚洲av天美| 超碰av人人做人人爽久久| 国产高清视频在线播放一区| 级片在线观看| 欧美成人性av电影在线观看| 日韩欧美国产一区二区入口| 亚洲av成人精品一区久久| 99久久无色码亚洲精品果冻| 欧美日本亚洲视频在线播放| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 啪啪无遮挡十八禁网站| 亚洲av二区三区四区|