• <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
    日本-黄色视频高清免费观看| 国产精品久久久久久精品电影小说 | 日韩欧美一区视频在线观看 | 日本vs欧美在线观看视频 | 极品教师在线视频| 免费大片黄手机在线观看| 亚洲精品久久久久久婷婷小说| 国产成人精品婷婷| 狂野欧美白嫩少妇大欣赏| 免费少妇av软件| 在线天堂最新版资源| 亚洲av成人精品一二三区| 视频中文字幕在线观看| 成年免费大片在线观看| 99国产精品免费福利视频| 日日啪夜夜撸| 免费大片18禁| 美女高潮的动态| 97超碰精品成人国产| 欧美xxⅹ黑人| 日韩中文字幕视频在线看片 | 1000部很黄的大片| 在线观看三级黄色| 2021少妇久久久久久久久久久| 狂野欧美白嫩少妇大欣赏| 美女视频免费永久观看网站| 欧美日本视频| 黄色日韩在线| 十分钟在线观看高清视频www | 精品一区二区三区视频在线| 久热这里只有精品99| av在线蜜桃| 十分钟在线观看高清视频www | 大香蕉97超碰在线| 人妻 亚洲 视频| 亚洲av欧美aⅴ国产| 亚洲自偷自拍三级| 日韩 亚洲 欧美在线| 青春草视频在线免费观看| 欧美日韩一区二区视频在线观看视频在线| 国产成人精品婷婷| 亚洲欧美日韩另类电影网站 | 亚洲国产av新网站| 王馨瑶露胸无遮挡在线观看| 高清毛片免费看| 久久久久国产网址| 亚洲国产精品专区欧美| 一本久久精品| 大又大粗又爽又黄少妇毛片口| 日本-黄色视频高清免费观看| 七月丁香在线播放| 简卡轻食公司| 人妻一区二区av| 国产午夜精品久久久久久一区二区三区| 青春草亚洲视频在线观看| 少妇 在线观看| av视频免费观看在线观看| 国产乱人视频| 精品少妇久久久久久888优播| 国产精品久久久久久精品古装| 亚洲欧美成人精品一区二区| av在线观看视频网站免费| 亚洲国产精品一区三区| 亚洲国产成人一精品久久久| 国产熟女欧美一区二区| 一本色道久久久久久精品综合| 黄片无遮挡物在线观看| 午夜日本视频在线| 一级毛片电影观看| 人妻系列 视频| 精品一品国产午夜福利视频| 亚洲国产高清在线一区二区三| 草草在线视频免费看| 校园人妻丝袜中文字幕| 妹子高潮喷水视频| 国产成人精品福利久久| 26uuu在线亚洲综合色| 欧美日韩视频高清一区二区三区二| 99热全是精品| 99国产精品免费福利视频| 汤姆久久久久久久影院中文字幕| 2018国产大陆天天弄谢| 99热6这里只有精品| 美女中出高潮动态图| 亚洲电影在线观看av| 搡女人真爽免费视频火全软件| 91在线精品国自产拍蜜月| 人妻系列 视频| 亚洲欧美成人精品一区二区| 国产黄频视频在线观看| 久久99热这里只频精品6学生| 久久精品国产亚洲网站| 国产成人精品福利久久| 亚洲精品色激情综合| 国产老妇伦熟女老妇高清| 欧美97在线视频| 在线观看av片永久免费下载| 人人妻人人添人人爽欧美一区卜 | 成年美女黄网站色视频大全免费 | 菩萨蛮人人尽说江南好唐韦庄| 国模一区二区三区四区视频| 久久国产精品男人的天堂亚洲 | a级一级毛片免费在线观看| 亚洲精品乱码久久久久久按摩| 国精品久久久久久国模美| 青春草国产在线视频| 亚洲,欧美,日韩| 日韩中字成人| 热99国产精品久久久久久7| 日韩电影二区| 亚洲欧美日韩无卡精品| 一区二区三区乱码不卡18| 久久久久久久亚洲中文字幕| 中文字幕av成人在线电影| 我要看日韩黄色一级片| 蜜桃久久精品国产亚洲av| 久久久久久久久久人人人人人人| 一边亲一边摸免费视频| 国产无遮挡羞羞视频在线观看| 亚洲最大成人中文| av免费在线看不卡| 丝袜脚勾引网站| 日韩在线高清观看一区二区三区| 国产伦精品一区二区三区视频9| 六月丁香七月| 在线观看免费视频网站a站| 这个男人来自地球电影免费观看 | 777米奇影视久久| 欧美精品亚洲一区二区| 欧美激情极品国产一区二区三区 | 久久久久久久亚洲中文字幕| 日产精品乱码卡一卡2卡三| av国产精品久久久久影院| 九九爱精品视频在线观看| 精品国产露脸久久av麻豆| 一级片'在线观看视频| 亚洲精品第二区| 国产探花极品一区二区| 国产亚洲91精品色在线| 人妻少妇偷人精品九色| 国产伦理片在线播放av一区| 亚洲综合色惰| 麻豆精品久久久久久蜜桃| 在线精品无人区一区二区三 | 欧美97在线视频| 婷婷色麻豆天堂久久| 久久精品久久精品一区二区三区| 国产免费又黄又爽又色| 国产一区二区三区av在线| 美女福利国产在线 | 久久精品久久久久久噜噜老黄| 成人无遮挡网站| 日本欧美视频一区| 黄片无遮挡物在线观看| 2021少妇久久久久久久久久久| 国产在线男女| www.色视频.com| 国产淫片久久久久久久久| 人妻一区二区av| 国产精品秋霞免费鲁丝片| 老司机影院成人| 欧美高清成人免费视频www| av国产免费在线观看| 一区二区三区乱码不卡18| 久久这里有精品视频免费| 街头女战士在线观看网站| 黄色怎么调成土黄色| 丝瓜视频免费看黄片| 欧美xxⅹ黑人| 欧美zozozo另类| 久久久久久久亚洲中文字幕| 蜜臀久久99精品久久宅男| 少妇猛男粗大的猛烈进出视频| 99热这里只有是精品50| 街头女战士在线观看网站| 精品久久久久久久末码| 久热这里只有精品99| 一区在线观看完整版| 亚洲欧美清纯卡通| 免费看光身美女| 丰满人妻一区二区三区视频av| 欧美xxxx性猛交bbbb| 亚洲国产日韩一区二区| 精品久久久久久久久av| 狠狠精品人妻久久久久久综合| 亚洲国产高清在线一区二区三| 国产高清国产精品国产三级 | 中国三级夫妇交换| 免费在线观看成人毛片| 国产男女超爽视频在线观看| 丰满乱子伦码专区| 欧美日韩亚洲高清精品| 日日摸夜夜添夜夜添av毛片| 最黄视频免费看| 久久韩国三级中文字幕| 晚上一个人看的免费电影| 在线天堂最新版资源| 国产高清国产精品国产三级 | 久久久国产一区二区| 干丝袜人妻中文字幕| 高清午夜精品一区二区三区| 在线观看国产h片| 日日啪夜夜爽| 在线观看免费日韩欧美大片 | 色视频在线一区二区三区| 免费黄频网站在线观看国产| 美女主播在线视频| 亚洲欧美一区二区三区国产| 久久国产精品大桥未久av | 在线免费观看不下载黄p国产| 精品少妇黑人巨大在线播放| 七月丁香在线播放| 91精品国产九色| 男女免费视频国产| 欧美成人a在线观看| 秋霞伦理黄片| 亚洲怡红院男人天堂| 少妇人妻精品综合一区二区| 人体艺术视频欧美日本| 欧美性感艳星| 久久这里有精品视频免费| 精品国产乱码久久久久久小说| 亚洲国产色片| 看十八女毛片水多多多| a级毛片免费高清观看在线播放| 久久久久久人妻| 日韩人妻高清精品专区| 国模一区二区三区四区视频| 欧美高清成人免费视频www| 99热网站在线观看| 亚洲国产成人一精品久久久| 日日啪夜夜撸| 男人和女人高潮做爰伦理| 成人高潮视频无遮挡免费网站| 一边亲一边摸免费视频| 久久久久久久久久久丰满| 婷婷色麻豆天堂久久| 一级毛片我不卡| 人妻一区二区av| 丰满人妻一区二区三区视频av| 女的被弄到高潮叫床怎么办| 亚洲精品日韩在线中文字幕| 国产精品久久久久久精品电影小说 | 国产欧美另类精品又又久久亚洲欧美| 成年av动漫网址| 天天躁夜夜躁狠狠久久av| 亚洲成人手机| 18+在线观看网站| 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 日韩人妻高清精品专区| 在现免费观看毛片| 看十八女毛片水多多多| 国产精品一区二区在线不卡| 中文字幕久久专区| 亚洲精品乱久久久久久| 高清毛片免费看| 只有这里有精品99| 欧美xxxx性猛交bbbb| 最近中文字幕2019免费版| 免费大片18禁| 亚洲精品一二三| 国产一区亚洲一区在线观看| 伊人久久国产一区二区| 在线观看一区二区三区| 下体分泌物呈黄色| 国产有黄有色有爽视频| 日本爱情动作片www.在线观看| 亚洲丝袜综合中文字幕| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品一区蜜桃| 免费观看a级毛片全部| 亚洲人与动物交配视频| 天堂中文最新版在线下载| 高清不卡的av网站| 中国美白少妇内射xxxbb| 精品午夜福利在线看| 欧美极品一区二区三区四区| 精品人妻视频免费看| 欧美精品亚洲一区二区| 你懂的网址亚洲精品在线观看| 91在线精品国自产拍蜜月| 久久久欧美国产精品| 日韩中文字幕视频在线看片 | 亚洲av福利一区| 亚洲精品乱码久久久v下载方式| 这个男人来自地球电影免费观看 | 校园人妻丝袜中文字幕| 熟女av电影| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲av天美| 免费人妻精品一区二区三区视频| 欧美xxxx性猛交bbbb| 欧美xxⅹ黑人| 国产免费一区二区三区四区乱码| 看十八女毛片水多多多| 啦啦啦视频在线资源免费观看| 熟女av电影| 色哟哟·www| 一级爰片在线观看| 18禁裸乳无遮挡动漫免费视频| 精品99又大又爽又粗少妇毛片| 亚洲精品,欧美精品| 国产精品99久久99久久久不卡 | 少妇高潮的动态图| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 成人特级av手机在线观看| 精品久久久噜噜| 亚洲综合色惰| 亚洲美女搞黄在线观看| 男女下面进入的视频免费午夜| 美女高潮的动态| 亚洲精品亚洲一区二区| 少妇 在线观看| 亚洲精品aⅴ在线观看| 国产在视频线精品| 国产有黄有色有爽视频| 插阴视频在线观看视频| 熟女电影av网| 亚洲欧美清纯卡通| 亚洲成人av在线免费| 91精品国产国语对白视频| 国产伦理片在线播放av一区| 丰满迷人的少妇在线观看| 国产在线男女| 国产视频首页在线观看| 亚洲精品一二三| 国产精品三级大全| 99久久精品一区二区三区| 亚洲国产欧美在线一区| 久久国产精品男人的天堂亚洲 | 在线观看一区二区三区激情| 免费不卡的大黄色大毛片视频在线观看| 九九爱精品视频在线观看| 免费观看a级毛片全部| 欧美丝袜亚洲另类| 一区二区三区精品91| 免费黄色在线免费观看| 2022亚洲国产成人精品| 天堂俺去俺来也www色官网| 久久久久久九九精品二区国产| 免费观看的影片在线观看| 欧美日韩综合久久久久久| 男男h啪啪无遮挡| 日本黄色日本黄色录像| 亚洲aⅴ乱码一区二区在线播放| 一个人看视频在线观看www免费| 亚洲av电影在线观看一区二区三区| 成人午夜精彩视频在线观看| 91久久精品国产一区二区三区| 一级毛片aaaaaa免费看小| 一本色道久久久久久精品综合| 伊人久久国产一区二区| 十分钟在线观看高清视频www | 久热这里只有精品99| 免费不卡的大黄色大毛片视频在线观看| 新久久久久国产一级毛片| 亚洲伊人久久精品综合| 免费黄频网站在线观看国产| 亚洲精品色激情综合| 国产色爽女视频免费观看| 一级毛片 在线播放| 国产乱人偷精品视频| 久久99热这里只有精品18| 久久精品国产a三级三级三级| 一区二区av电影网| videos熟女内射| 国产精品一及| 中文字幕人妻熟人妻熟丝袜美| 精品一区二区免费观看| 久久久a久久爽久久v久久| 一级毛片 在线播放| 中文字幕人妻熟人妻熟丝袜美| 亚洲综合色惰| 老女人水多毛片| 国产精品国产三级国产av玫瑰| 在线观看人妻少妇| 国产一区亚洲一区在线观看| 国产成人精品一,二区| 精品午夜福利在线看| 夜夜爽夜夜爽视频| 欧美另类一区| 91精品一卡2卡3卡4卡| 搡老乐熟女国产| 日韩在线高清观看一区二区三区| 国产精品精品国产色婷婷| 中国美白少妇内射xxxbb| 亚洲国产精品一区三区| 人人妻人人爽人人添夜夜欢视频 | 国模一区二区三区四区视频| 日本免费在线观看一区| 2021少妇久久久久久久久久久| 亚洲欧美精品自产自拍| 国产精品99久久99久久久不卡 | .国产精品久久| 六月丁香七月| 精品熟女少妇av免费看| 精品一区二区三区视频在线| av网站免费在线观看视频| 最近2019中文字幕mv第一页| 国产精品福利在线免费观看| 国产精品蜜桃在线观看| 26uuu在线亚洲综合色| 亚洲av不卡在线观看| 视频区图区小说| 精品亚洲乱码少妇综合久久| 丰满人妻一区二区三区视频av| 一个人免费看片子| 有码 亚洲区| 亚洲av中文字字幕乱码综合| 免费大片18禁| 看十八女毛片水多多多| 国产成人a∨麻豆精品| 性高湖久久久久久久久免费观看| 18禁裸乳无遮挡免费网站照片| 天堂中文最新版在线下载| 国产老妇伦熟女老妇高清| 午夜福利高清视频| 美女高潮的动态| 日本色播在线视频| 亚洲不卡免费看| 少妇人妻久久综合中文| 久久久久久久亚洲中文字幕| 一边亲一边摸免费视频| 日日摸夜夜添夜夜添av毛片| 啦啦啦视频在线资源免费观看| 夫妻性生交免费视频一级片| 中文字幕精品免费在线观看视频 | 国产精品久久久久久精品电影小说 | av卡一久久| 夜夜看夜夜爽夜夜摸| 亚洲国产精品成人久久小说| 国产黄频视频在线观看| 亚洲图色成人| 一级片'在线观看视频| 免费看不卡的av| 久久99热这里只频精品6学生| 亚洲成色77777| 亚洲四区av| 男女啪啪激烈高潮av片| 美女视频免费永久观看网站| 日日撸夜夜添| 成人午夜精彩视频在线观看| 欧美变态另类bdsm刘玥| 九色成人免费人妻av| 亚洲欧美成人精品一区二区| 女人十人毛片免费观看3o分钟| 国产男女超爽视频在线观看| 久久久久精品性色| 国产黄频视频在线观看| 亚洲国产精品国产精品| 麻豆国产97在线/欧美| 亚洲性久久影院| 国产成人a∨麻豆精品| 亚洲电影在线观看av| 欧美精品亚洲一区二区| 亚洲精品日韩av片在线观看| 国产大屁股一区二区在线视频| 国产男女内射视频| 一个人看视频在线观看www免费| 嫩草影院入口| 成人毛片a级毛片在线播放| 韩国高清视频一区二区三区| 搡女人真爽免费视频火全软件| 国产精品麻豆人妻色哟哟久久| a级毛片免费高清观看在线播放| 日本黄色日本黄色录像| 我要看黄色一级片免费的| 亚洲激情五月婷婷啪啪| 精品99又大又爽又粗少妇毛片| 亚洲精品乱久久久久久| 亚洲欧美成人精品一区二区| 欧美成人a在线观看| 国产精品一区二区性色av| 三级国产精品片| 免费播放大片免费观看视频在线观看| 夜夜爽夜夜爽视频| 综合色丁香网| 一级毛片黄色毛片免费观看视频| 少妇熟女欧美另类| 九九久久精品国产亚洲av麻豆| 美女中出高潮动态图| 青青草视频在线视频观看| av天堂中文字幕网| 51国产日韩欧美| 3wmmmm亚洲av在线观看| 在线观看一区二区三区激情| 18禁在线播放成人免费| 九九在线视频观看精品| 欧美xxxx性猛交bbbb| 婷婷色综合大香蕉| 久久国产精品男人的天堂亚洲 | www.av在线官网国产| 久久国内精品自在自线图片| 午夜福利在线在线| 国产在线一区二区三区精| 精品久久国产蜜桃| 大香蕉久久网| 下体分泌物呈黄色| 精品国产乱码久久久久久小说| 精品久久久噜噜| 青青草视频在线视频观看| 欧美日韩综合久久久久久| 51国产日韩欧美| 久久久久精品久久久久真实原创| 成人一区二区视频在线观看| 亚洲,一卡二卡三卡| 在线看a的网站| 成人二区视频| 大码成人一级视频| 亚洲在久久综合| 中文字幕精品免费在线观看视频 | 中文在线观看免费www的网站| 成人特级av手机在线观看| 女人久久www免费人成看片| 国产片特级美女逼逼视频| 亚洲美女黄色视频免费看| 最近最新中文字幕大全电影3| 我要看日韩黄色一级片| 欧美区成人在线视频| 99热全是精品| 日本vs欧美在线观看视频 | 只有这里有精品99| 老师上课跳d突然被开到最大视频| 成人影院久久| 99视频精品全部免费 在线| 久久久久性生活片| 色吧在线观看| 人人妻人人爽人人添夜夜欢视频 | 国产视频内射| 亚洲一级一片aⅴ在线观看| 国产伦理片在线播放av一区| 日本一二三区视频观看| 亚洲美女黄色视频免费看| tube8黄色片| 在线观看美女被高潮喷水网站| 肉色欧美久久久久久久蜜桃| 成人二区视频| 一级片'在线观看视频| 免费大片黄手机在线观看| av女优亚洲男人天堂| 欧美bdsm另类| 久久久精品免费免费高清| 日日摸夜夜添夜夜爱| 亚洲精品aⅴ在线观看| 纯流量卡能插随身wifi吗| 麻豆国产97在线/欧美| 国产白丝娇喘喷水9色精品| 国产黄色视频一区二区在线观看| 日本wwww免费看| kizo精华| 高清欧美精品videossex| 少妇人妻一区二区三区视频| 中文字幕久久专区| 偷拍熟女少妇极品色| 九草在线视频观看| 性色avwww在线观看| 国产亚洲5aaaaa淫片| av在线老鸭窝| 免费av中文字幕在线| 国产成人一区二区在线| 91精品伊人久久大香线蕉| 亚洲精品乱码久久久久久按摩| 国产乱人偷精品视频| 偷拍熟女少妇极品色| 亚洲欧美一区二区三区黑人 | 国产男女超爽视频在线观看| a级毛色黄片| 亚洲精品色激情综合| 九草在线视频观看| 97精品久久久久久久久久精品| 亚洲av.av天堂| 在线精品无人区一区二区三 | 亚洲成人中文字幕在线播放| 免费人成在线观看视频色| 天天躁日日操中文字幕| 一区二区三区精品91| 成人美女网站在线观看视频| 日日摸夜夜添夜夜爱| 国产黄色免费在线视频| 亚洲国产毛片av蜜桃av| 亚洲色图综合在线观看| 91午夜精品亚洲一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 综合色丁香网| 国产爱豆传媒在线观看| 极品少妇高潮喷水抽搐| 天美传媒精品一区二区| 久久热精品热| 精品亚洲成国产av| 亚洲色图综合在线观看| 亚洲精品日韩在线中文字幕| tube8黄色片| 欧美zozozo另类| 这个男人来自地球电影免费观看 | av黄色大香蕉| 黄色怎么调成土黄色| 天天躁夜夜躁狠狠久久av| 精品亚洲成国产av| 亚洲,一卡二卡三卡| 久久久午夜欧美精品| 久久国产精品男人的天堂亚洲 | 日韩av不卡免费在线播放| 乱码一卡2卡4卡精品| 最近中文字幕高清免费大全6| 日韩av不卡免费在线播放| 最近最新中文字幕大全电影3| 亚洲色图综合在线观看| 这个男人来自地球电影免费观看 | 日日啪夜夜爽| 校园人妻丝袜中文字幕| 日韩av不卡免费在线播放| tube8黄色片| 国产成人a∨麻豆精品|