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

    浮游植物種群對海洋酸化和光照強(qiáng)度變化的響應(yīng):以長江口南毗鄰海域?yàn)槔?/h1>
    2023-11-15 11:12:58李童童馮媛媛王建才王雅婻宋書群白有成
    海洋科學(xué)進(jìn)展 2023年4期
    關(guān)鍵詞:硅藻高光拷貝數(shù)

    李童童,馮媛媛,王建才,蔡 婷,王雅婻,宋書群,白有成

    (1. 天津科技大學(xué) 海洋與環(huán)境學(xué)院,天津 300457;2. 上海交通大學(xué) 海洋學(xué)院,上海 200030;3. 上海市極地前沿科學(xué)研究基地,上海 200030;4. 中國科學(xué)院 海洋研究所,山東 青島 266071;5. 自然資源部 第二海洋研究所,浙江 杭州 310012)

    海洋是重要的無機(jī)碳儲(chǔ)庫,隨著CO2含量逐漸升高,海水碳酸鹽體系受到影響,海水的pH隨之降低,導(dǎo)致海洋酸化(Zeebe, 2012)。海洋浮游植物作為海洋中最主要的初級生產(chǎn)者,通過光合作用將CO2轉(zhuǎn)化為有機(jī)碳,并通過食物網(wǎng)傳輸及碳沉降等過程使部分有機(jī)碳被轉(zhuǎn)移并埋藏至深海中(Ducklow et al, 2001)。由于當(dāng)今海水中游離的CO2濃度低于大部分浮游植物光合作用所需的飽和CO2濃度,大部分浮游植物種群進(jìn)化產(chǎn)生了碳濃縮機(jī)制(Carbon Concentration Mechanism, CCM)(Badger et al, 1992),可以同時(shí)使用海水中的CO2和HCO3?作為進(jìn)行光合作用的無機(jī)碳源,升高核酮糖-1, 5-二磷酸羧化酶/加氧酶(Rubisco)附近的CO2濃度,以維持細(xì)胞較高的光合作用效率(Falkowski et al, 2013)。由于不同浮游植物種群的CCM效率不同,且受海水CO2濃度的影響(Wu et al, 2010),CO2濃度升高對某些浮游植物物種光合作用的促進(jìn)效率可能會(huì)高于其他物種,對較大粒徑的浮游植物影響尤其明顯(Wu et al, 2014),從而導(dǎo)致海洋浮游植物的群落結(jié)構(gòu)發(fā)生變化(Moroney et al,1999)。

    光照強(qiáng)度也是影響浮游植物光合作用的重要環(huán)境因素(Singh et al, 2015)。浮游植物需要足夠的光照強(qiáng)度來進(jìn)行高效的光合作用,光照強(qiáng)度對浮游植物的色素復(fù)合體合成以及光合固碳作用具有顯著影響(Hihara et al, 1998)。過高光照強(qiáng)度通常會(huì)限制光合色素的生產(chǎn)從而降低浮游植物光合作用的效率(陶宗婭等, 1999),較低光照強(qiáng)度可刺激浮游植物葉綠素的生產(chǎn),但同時(shí)導(dǎo)致光合作用速率降低(Beardall et al, 1976);不同浮游植物種群的光合作用最適宜光照強(qiáng)度有所不同(Thompson, 2010),因此光照強(qiáng)度升高會(huì)促進(jìn)某些未到達(dá)光飽和點(diǎn)的浮游植物種群的光合作用,但可能抑制一些光飽和點(diǎn)比較低的浮游植物種群的光合固碳,影響種間競爭和演替(Anderson et al, 1997)。光照強(qiáng)度變化也會(huì)影響浮游植物細(xì)胞合成脂質(zhì),細(xì)胞在較高光照強(qiáng)度下更易累積飽和脂肪酸,在低光下則更多累積不飽和脂肪酸(Khotimchenko et al, 2005)。

    光照強(qiáng)度往往與氣候變化背景下快速變化的其他環(huán)境因素同時(shí)作用于該海域的初級生產(chǎn)者并產(chǎn)生交互效應(yīng),影響生物固碳及相關(guān)元素的生物地球化學(xué)循環(huán)過程(Finkel et al, 2006)。研究表明,南大洋2株硅藻短擬脆桿藻(Fragilariopsis curta)和魏氏齒盤藻(Odontella weisflogii)對CO2分壓(pCO2)的變化很敏感,較高的光照強(qiáng)度可促進(jìn)這2株硅藻在低pCO2下顆粒有機(jī)碳的生產(chǎn)(Heiden et al,2016)。在浮游植物群落的層面,海洋酸化與光照強(qiáng)度變化協(xié)同降低了南海的浮游植物固碳作用(Gao et al, 2012),并可改變浮游植物種群組成(Domingues et al, 2014)。

    東海是我國的邊緣海,特別是長江口南毗鄰海域的生產(chǎn)力較高,在海洋碳循環(huán)過程中有重要的作用(Tseng et al, 2011)。受長江沖淡水的影響,長江口海域理化環(huán)境比較復(fù)雜,長江徑流帶來的大量懸浮泥沙及海流擾動(dòng)產(chǎn)生的再懸浮泥沙導(dǎo)致該海域表層形成透明度3 m及以下的渾濁帶(Fang et al, 2006)。長江口毗鄰海域沿鹽度梯度可分為近河口光限制區(qū)、過渡帶光和磷酸鹽限制區(qū)、遠(yuǎn)河口的氮鹽限制區(qū)(Filardo et al, 1985)。其中光限制區(qū)雖然具有較高含量的營養(yǎng)鹽水平,但渾濁帶的消光作用導(dǎo)致該區(qū)域的初級生產(chǎn)力潛在受到光限制(何文珊等, 2001)。我們以長江口南毗鄰海域浮游植物群落為研究對象,設(shè)置不同pCO2和光照強(qiáng)度條件,對長江口海域浮游植物群落進(jìn)行現(xiàn)場模擬培養(yǎng)實(shí)驗(yàn),研究海洋酸化與光照強(qiáng)度變化對浮游植物群落組成及相關(guān)生物地球化學(xué)參數(shù)的交互效應(yīng),以期分析和預(yù)測該海域浮游植物群落對海洋酸化和光照環(huán)境變化的響應(yīng)趨勢。

    1 材料和方法

    1.1 采樣和培養(yǎng)實(shí)驗(yàn)設(shè)置

    2020年春季在“向陽紅18”科學(xué)考察船上現(xiàn)場采樣并進(jìn)行了為期5 d的現(xiàn)場船基受控培養(yǎng)實(shí)驗(yàn)(Hutchins et al, 2003)。2020-05-20于站位S01-1(122°42′00″E,29°59′54″N)采集表層海水水樣,采樣站位水深35 m,表層海水水溫:20.64 ℃,鹽度:29.73;營養(yǎng)鹽濃度:1.15 μmol·L?1、7.84 μmol·L?1、0.42 μmol·L?1、-Si 20.33 μmol·L?1(圖1)。

    圖1 采樣站位()示意圖Fig. 1 Schamatic diagram of sampling station ()

    將采集的海水水樣經(jīng)孔徑200 μm的篩絹過濾,分裝至12個(gè)4.5 L酸洗后的潔凈透明聚碳酸酯培養(yǎng)瓶中,并進(jìn)行營養(yǎng)鹽加富(添加濃度為μmol·L?1、1 μmol·L?1、μmol·L?1)。通過CO2加富器(CE-100型,中國武漢瑞華有限公司生產(chǎn))調(diào)節(jié)CO2和空氣混合比例,并使用鼓泡曝氣的方法(Mcgraw et al, 2010)控制培養(yǎng)瓶中的pCO2分別為模擬現(xiàn)代空氣中的41 Pa和預(yù)測世紀(jì)末的101 Pa(Caldeira et al, 2005)。我們通過測定,采樣站位表層海水的光照強(qiáng)度為海表光照強(qiáng)度的30%~50%。參照Feng等的工作(Feng et al, 2010),為了分別模擬表層海水及由于混合層變淺導(dǎo)致的混合層內(nèi)光照強(qiáng)度升高,通過在培養(yǎng)槽上覆蓋中型密度網(wǎng)控制培養(yǎng)瓶內(nèi)光照強(qiáng)度為50%和75%海表光照強(qiáng)度,共設(shè)置4個(gè)實(shí)驗(yàn)處理:①對照組(41 Pa,50%海表光照強(qiáng)度);②酸化組(101 Pa,50%海表光照強(qiáng)度);③高光組(41 Pa,75% 海表光照強(qiáng)度);④交互組(101 Pa,75% 海表光照強(qiáng)度)。每組3個(gè)平行樣,培養(yǎng)水槽中的水浴溫度使用冷暖水機(jī)(HC-1000BH型,中國海利有限公司生產(chǎn))控制為(20 ± 1) ℃。培養(yǎng)期間每天午間采用光照數(shù)據(jù)采集器(LI-1500型,美國LICOR有限公司生產(chǎn))監(jiān)測培養(yǎng)槽內(nèi)光照強(qiáng)度(圖2)。每天上午08:00?10:00對各培養(yǎng)瓶總?cè)~綠素(Chla)質(zhì)量濃度測定,隔天采樣測定粒徑分級Chla質(zhì)量濃度(粒徑范圍分別為0.2~ <2.0 μm,2.0~ <20.0 μm,≥20.0 μm),在最終取樣日(第5天)收集用于測定顆粒有機(jī)磷(POP)、顆粒有機(jī)碳(POC)、顆粒有機(jī)氮(PON)、生物硅(BSi)濃度的樣品、微微型浮游植物樣品、16S/18S rRNA擴(kuò)增子測序樣品以及浮游植物群落顯微鏡檢查樣品。

    圖2 培養(yǎng)過程中光照強(qiáng)度Fig. 2 Light intensity during the course of incubation

    1.2 海水碳酸鹽化學(xué)測定

    每天上午08:00?10:00測定船基受控培養(yǎng)樣品的pH(Mcgraw et al, 2010)和總堿度(Dickson et al,2007)。采用pH計(jì)(SevenCompactTM S210K型,瑞士Mettler Toledo有限公司生產(chǎn))測定培養(yǎng)體系pH(Mcgraw et al, 2010)。海水總堿度采用總堿度自動(dòng)滴定儀(G20S型,瑞士Mettler Toledo有限公司生產(chǎn))測定(Dickson et al, 2007)。海水碳酸鹽化學(xué)各參數(shù)根據(jù)樣品的pH和總堿度使用CO2SYS軟件(Lewis et al, 1998)計(jì)算獲取(Dickson et al, 1987)。

    1.3 浮游植物種群組成分析

    使用體積分?jǐn)?shù)40%的甲醛溶液(購自國藥平臺(tái))對采集后的浮游植物群落樣品進(jìn)行固定,用量為水樣體積的4%。固定后常溫保存帶回實(shí)驗(yàn)室。2020-06-02在天津科技大學(xué)海洋與環(huán)境學(xué)院實(shí)驗(yàn)室采用Uterm?hl計(jì)數(shù)法檢測(Uterm?hl, 1958)。將樣品沉降24 h后放置在倒置顯微鏡(AE2000型,廈門麥克奧迪有限公司生產(chǎn))下進(jìn)行定種和細(xì)胞計(jì)數(shù)。

    1.4 16s/18s rRNA高通量測序

    測序樣品過濾到孔徑0.6 μm聚碳酸酯濾膜(購自美國Millipore有限公司),保存于液氮罐(?180℃)帶回實(shí)驗(yàn)室后保存在超低溫(?180 ℃)冰箱中。2021-03-24委托天津諾禾致源生物信息科技有限公司采用CTAB/SDS(十六烷基三甲基溴化銨/十二烷基磺酸鈉)方法提取樣品基因組的DNA。利用PCR儀(T100型,美國Bio-rad有限公司生產(chǎn))擴(kuò)增出不同區(qū)域的16S/18S rDNA基因,產(chǎn)物經(jīng)檢測后使用GeneJET試劑盒(購自美國Thermo Scientific有限公司)進(jìn)行純化,再使用TruSeq?DNA PCR-Free Sample Preparation Kit試劑盒(購自美國Illumina有限公司)構(gòu)建文庫,經(jīng)Qubit(美國Thermo Scientific有限公司生產(chǎn))定量和文庫檢測合格后使用NovaSeq6000進(jìn)行上機(jī)測序。對原數(shù)據(jù)進(jìn)行拼接、過濾后得到有效數(shù)據(jù)(Mago et al, 2011),利用Uparse算法(UPARSE v7.0.1001)(Haas et al, 2011)對所有樣品的有效數(shù)據(jù)進(jìn)行聚類,以97%的一致性將序列聚類成基因序列分類單元(Operational Taxonomic Units,OTUs),得到對應(yīng)的物種信息和基于物種的豐度分布情況(Edgar, 2013)。

    1.5 微微型浮游植物種群組成分析

    將含有浮游植物樣品的海水經(jīng)孔徑20 μm篩絹后取1.8 mL,使用200 μL多聚甲醛溶液(體積分?jǐn)?shù)為4%,購自國藥平臺(tái))進(jìn)行固定后保存在液氮罐(?180 ℃)中,2021-09-23在中國科學(xué)院海洋研究所實(shí)驗(yàn)室使用流式細(xì)胞儀(Accuri C6型,美國Dickinson and Company有限公司生產(chǎn))進(jìn)行測定(Jiang et al, 2017)。通過流式細(xì)胞儀發(fā)出的488 nm激光照射浮游植物樣品后反饋的不同信號(hào)來劃分不同類群:微微型真核藻類的FSC-H信號(hào)小于作為內(nèi)標(biāo)的2 μm綠色熒光微球,但是FSC-H和FL3-H信號(hào)大于原核自養(yǎng)生物;聚球藻含有藻紅蛋白,所以FL2-H信號(hào)明顯高于其他類群;原綠球藻細(xì)胞最小,葉綠素含量比異養(yǎng)細(xì)菌的高,可根據(jù)FSC-H和FL3-H信號(hào)強(qiáng)度進(jìn)行區(qū)分(Jiang et al, 2017)。

    1.6 浮游植物群落Chl a質(zhì)量濃度測定

    樣品采用微壓(<0.04 MPa)過濾樣品至孔徑0.2 μm GF/F膜(購自美國Whatman有限公司生產(chǎn))上用于測定樣品的總Chla質(zhì)量濃度,粒徑分級Chla的樣品依次通過孔徑為20 μm、2 μm的聚碳酸酯濾膜(購自美國Millipore有限公司)和GF/F膜過濾采樣,于?20 ℃保存。2020-06-02在天津科技大學(xué)海洋與環(huán)境學(xué)院實(shí)驗(yàn)室進(jìn)行分析,采用體積分?jǐn)?shù)為90%的丙酮溶液萃取后,用Turner熒光儀(美國Turner Designs有限公司生產(chǎn))酸化模塊(Welschmeyer, 1994)測定細(xì)胞Chla質(zhì)量濃度。

    1.7 浮游植物群落元素組成測定

    POC和PON樣品過濾至灼燒(450 ℃,4 h)過的GF/F膜上常溫保存,2021-03-30在自然資源部第二海洋研究所烘干后采用元素分析儀(ECS4010型,意大利Costech有限公司生產(chǎn))測定(Hutchins et al, 1998)。POP樣品經(jīng)灼燒過的GF/F膜過濾后轉(zhuǎn)入20 mL同樣灼燒過的玻璃瓶中,加入2 mL濃度為0.017 mol·L?1的MgSO4溶液,于60 ℃烘干后常溫保存。樣品于2020-06-14在天津科技大學(xué)海洋與環(huán)境學(xué)院實(shí)驗(yàn)室使用紫外可見分光光度計(jì)(UV-2550型,日本島津有限公司生產(chǎn))采用磷鉬藍(lán)顯色法測定(Solórzano et al, 1980)。BSi樣品用0.6 μm孔徑聚碳酸酯濾膜(購自美國Millipore有限公司)過濾后裝入10 mL離心管中常溫保存,于測定POP樣品的第2天放入60 ℃的烘箱內(nèi)烘干后用1 mol·L?1的HCl溶液消解后采用硅鉬藍(lán)顯色法測定樣品中BSi的含量(Brzezinski et al, 1995)。

    1.8 浮游植物群落沉降速率測定

    船基培養(yǎng)實(shí)驗(yàn)的第0天和第5天在船上采用同源性采樣法(SETCOL)測定浮游植物群落的沉降速率(Bienfang, 1981)。根據(jù)Bienfang的沉降公式(Bienfang, 1981)計(jì)算浮游植物群落藻細(xì)胞的沉降速率。

    1.9 數(shù)據(jù)分析

    單一環(huán)境因素酸化(A)、高光(L)及2個(gè)環(huán)境因素共同變化(A+L)下對某生理生化參數(shù)的觀測影響(observed effect, OE)依照該處理組與對照組的變化百分比計(jì)算而得(正值為正向升高效應(yīng),負(fù)值為負(fù)向降低效應(yīng)),二者的交互效應(yīng)(multiplicative effect, ME)采用Folt等方法(Feng et al, 2020)計(jì)算。當(dāng)OEA和OEL均為正或均為負(fù)時(shí),在|OEA+L| > |MEA+L|條件下,2種環(huán)境因素為協(xié)同性交互效應(yīng),在|OEA+L| < |MEA+L|條件下,2種環(huán)境因素為拮抗性交互效應(yīng);當(dāng)OEA和OEL為一正一負(fù)時(shí),2種環(huán)境因素為拮抗性交互效應(yīng)。采用雙因素方差分析(ANOVA)檢驗(yàn)變量之間的交互作用(α=0.05水平),組間顯著性差異采用t-test分析。

    2 結(jié) 果

    2.1 碳酸鹽化學(xué)參數(shù)

    在培養(yǎng)期間,不同pCO2處理組之間的pH存在顯著差異。最終采樣時(shí)各處理組海水碳酸鹽化學(xué)各參數(shù)見表1。

    表1 培養(yǎng)第5天時(shí)培養(yǎng)瓶中碳酸鹽化學(xué)參數(shù)Table 1 Chemical parameters of carbonate in the four experimental treatment groups on the fifth day incubation

    2.2 Chl a質(zhì)量濃度

    對照組和交互組的總Chla質(zhì)量濃度均在培養(yǎng)第2天達(dá)到最大值(7.84 ± 1.00) μg·L?1,酸化組在第3天達(dá)到(9.21 ± 2.87) μg·L?1,高光組最大值(9.48 ± 1.47) μg·L?1出現(xiàn)在第4天(圖3a)。與其他處理組相比,高光組中的總Chla質(zhì)量濃度在第4天和第5天均為各組中最高。

    最終采樣日(第5天)粒徑分級Chla質(zhì)量濃度結(jié)果表明,在浮游植物種群中粒徑范圍2.0~ <20.0 μm的微型浮游植物占比最高,對照組為(61 ± 6)%、酸化組為(60 ± 2)%、高光組為(52 ±3)%、交互組為(64 ± 6)%,均超過50%。酸化主要使粒徑≥20.0 μm的小型浮游植物Chla質(zhì)量濃度占比升高,隨光照強(qiáng)度升高小型和微微型(0.2~ <2.0 μm)浮游植物Chla質(zhì)量濃度占比升高。小型浮游植物Chla質(zhì)量濃度在對照組最低為(30 ± 9)%,酸化組(34 ± 2)%較對照組上升了4%,而高光組(37 ±6)%上升了7%,交互組(31 ± 5)%上升了1%;微微型浮游植物Chla質(zhì)量濃度占比最低,與對照組的(10 ± 3)%相比,高光組的(11 ± 4)%升高2%,酸化組的(7 ± 1)%和交互組的(5 ± 1)%分別降低3%和4%(圖3b)。

    2.3 浮游植物種群組成

    酸化顯著提升了微微型浮游植物的細(xì)胞豐度,對聚球藻(Synechococcus)的影響尤其明顯。流式細(xì)胞儀的檢測結(jié)果表明酸化組和交互組分別是對照組[(426 ± 177)個(gè)·L?1]的32.72倍和31.58倍(p<0.05,表2),而高光組[(864 ± 46)個(gè)·L?1]較對照組升高10.28%(p<0.05);真核微微型浮游植物只在交互組顯著高于對照組(p<0.05)。在各處理中未有原綠球藻(Prochlorococcus)信號(hào)檢出。

    表2 基于流式細(xì)胞儀檢測的培養(yǎng)第5天4個(gè)處理組的微微型浮游植物細(xì)胞豐度(個(gè)·L?1)Table 2 Abundances of microtype phytoplankton cells based on flowcytometry on the fifth day incubation (cells·L?1)

    顯微鏡檢結(jié)果表明硅藻中海鏈藻(Thalassiosiraspp.)、甲藻中卡爾藻(Karlodiniumsp.)、原甲藻(Prorocentrumsp.)在各組均為優(yōu)勢種,海鏈藻在對照組、酸化組、高光組和交互組中的占比分別為42%、59%、66%和56%。對照組中海鏈藻豐度為5.60×105個(gè)·L?1,酸化組、高光組和交互組分別比對照組降低48%、34%和30%;卡爾藻在對照組中豐度為2.66×105個(gè)·L?1,酸化組和交互組分別比對照組升高5%和11%,高光組比對照組降低53%;原甲藻在酸化組、高光組和交互組中的豐度分別比對照組2.95×104個(gè)·L?1降低76%、91%和94%(圖4)。篩鏈藻(Coscinosira polychorda)、柔弱擬菱形藻(Pseudo-nitzschia delicatissima)和冠蓋藻(Stephanopyxisspp.)分別在酸化組、高光組和交互組中占據(jù)一定優(yōu)勢。由于有效樣品數(shù)量不足,最終數(shù)據(jù)間沒有體現(xiàn)出顯著性差異。

    圖4 培養(yǎng)第 5 天顯微鏡檢的 4 個(gè)處理組浮游植物優(yōu)勢物種豐度和各處理組的物種組成百分比Fig. 4 Microscopic dominant phytoplankton species composition in the batch culture and the percentage of the phytoplankton species composition in the fourexperimental treatment groups on the fifth day incubation

    高通量測序結(jié)果與顯微鏡檢結(jié)果相似,酸化組、高光組和交互組中硅藻的基因拷貝數(shù)均高于對照組,在高光組尤為顯著,較對照組提升了231%(圖5a);金藻的基因拷貝數(shù)在酸化組最高,較對照組提升了3464%,而高光組和交互組比對照組分別提升了237%和1128%(圖5b);藍(lán)藻的基因拷貝數(shù)最高值也出現(xiàn)在交互組,酸化組、高光組和交互組分別比對照組增長156%、78%和359%(圖5c);與對照組相比,甲藻的基因拷貝數(shù)在其他3個(gè)處理組均有降低,在酸化組降低程度最高(68%,圖5d);浮游植物群落中隱藻的基因拷貝數(shù)的變化和硅藻相似,在高光組的基因拷貝數(shù)最高(圖5e),由于隱藻基因拷貝數(shù)在對照組、酸化組、高光組和交互組中組成百分比均僅為0.2%及以下,故不體現(xiàn)在各種群組成百分比中(圖5f)。

    圖5 培養(yǎng)第 5 天高通量測序的 4 個(gè)處理組基因拷貝數(shù)和各種群組成百分比Fig. 5 Gene copies and the percentage of different populations in the four experimental treatment groups on the fifth day incubation based on high-throughput sequencing

    對照組中共發(fā)現(xiàn)浮游植物5門85種。與顯微鏡檢結(jié)果相似,硅藻中海鏈藻屬(Thalassiosira)基因拷貝數(shù)最高,甲藻中卡爾藻屬(Karlodinium)、金藻中近囊胞藻屬(Paraphysomonas)、藍(lán)藻中聚球藻屬(Synechococcus)基因拷貝數(shù)較高。酸化組記錄5門103種,其中硅藻中海鏈藻屬基因拷貝數(shù)最高,甲藻中別什萊藻屬(Biecheleria)和原甲藻屬(Prorocentrun)基因拷貝數(shù)相近,對照組中的占優(yōu)的卡爾藻屬基因拷貝數(shù)比對照組降低95%;金藻近囊胞藻屬的基因拷貝數(shù)則為對照組的40倍;藍(lán)藻中聚球藻屬基因拷貝數(shù)較高,比對照組升高151%,聚球藻屬的拷貝數(shù)變化趨勢與流式細(xì)胞儀檢測結(jié)果相吻合(表2)。高光組中記錄了5門92種,其中硅藻海鏈藻屬基因拷貝數(shù)最高,是對照組的2倍(圖4e),這個(gè)結(jié)果同顯微鏡檢結(jié)果相似;甲藻中卡爾藻屬比對照組降低52%;金藻近囊胞藻屬比對照組升高232%;藍(lán)藻優(yōu)勢屬聚球藻屬比對照組升高72%。交互組記錄5門97種,其中硅藻海鏈藻屬高于對照組45%(圖4e),與顯微鏡檢結(jié)果相印證;甲藻中別什萊藻屬的基因拷貝數(shù)為對照組的18倍,其次為華鰭藻屬(Sinophysis),其基因拷貝數(shù)為對照組的167倍之多;金藻中近囊胞藻屬的拷貝數(shù)為對照組10倍;與流式細(xì)胞儀結(jié)果相似(表2),藍(lán)藻中聚球藻屬基因拷貝數(shù)最高,且在4個(gè)處理組中交互組基因拷貝數(shù)最高,次高出現(xiàn)在酸化組(圖6)。

    圖6 培養(yǎng)第5天4個(gè)處理組的浮游植物優(yōu)勢物種組成熱圖(16s/18s rRNA擴(kuò)增子測序基因拷貝數(shù))Fig. 6 Heat map of main phytoplankton dominant species composition in the batch culture in the four experimental treatment groups on the fifth day incubation (16S/18S rRNA amplicon sequencing copy numbers)

    2.4 元素組成和沉降速率

    高光組和交互組的POC、PON及BSi濃度高于對照組的及酸化組的,最高值均出現(xiàn)在交互組。但在各處理組間POC、PON、POP及BSi濃度及元素比由于數(shù)據(jù)間方差較大,沒有顯示出統(tǒng)計(jì)學(xué)顯著差異(p>0.05, 表3)。

    表3 培養(yǎng)第5天4個(gè)處理組中各元素濃度與元素摩爾比Table 3 Elemental concentrations and elemental molar ratios in the four experimental treatments groups on the fifth day of the incubation

    培養(yǎng)實(shí)驗(yàn)表明,酸化條件對浮游植物種群的沉降有顯著作用(圖7)。培養(yǎng)第5天酸化組的沉降速率為(1.47 ± 0.13) m·d?1,顯著高于其他3個(gè)處理組(p<0.05),是對照組的3倍;交互組的沉降速率為高光組的2倍。

    圖7 培養(yǎng)第5天4個(gè)處理組的浮游植物種群沉降速率(誤差棒代表平行樣間標(biāo)準(zhǔn)偏差,n=3)Fig. 7 Sinking rates of phytoplankton composition in batch culture of four experimental treatment groups on the fifth day incubation(the error bars represent standard deviations of triplicate samples, n=3)

    2.5 雙因素交互作用

    計(jì)算結(jié)果表明酸化和高光照強(qiáng)度對POC濃度、POP濃度具有協(xié)同性交互效應(yīng),而對總Chla質(zhì)量濃度、小型、微型和微微型浮游植物Chla質(zhì)量濃度、PON濃度、BSi濃度、C/N比值、N/P比值、C/P比值、C/Si比值、沉降速率、硅藻以及甲藻豐度產(chǎn)生拮抗性交互效應(yīng)(表4)。

    3 討論

    我們的研究結(jié)果表明海洋酸化程度與光照強(qiáng)度變化對浮游植物生長及種群組成產(chǎn)生了不同的影響,2種環(huán)境要素的變化可產(chǎn)生交互效應(yīng)。在350~1200 μmol·photons·m?2·s?1光照強(qiáng)度下,光照強(qiáng)度升高促進(jìn)浮游植物的生長及生物量積累,如中心綱硅藻的豐度在高光條件下顯著升高。pCO2為101 Pa的酸化更有利于聚球藻的生長,可顯著降低微微型及微型浮游植物的Chla質(zhì)量濃度組成百分比,可提升小型浮游植物Chla質(zhì)量濃度占比,浮游植物種群粒徑增大導(dǎo)致沉降速率升高。酸化和高光均可降低甲藻的優(yōu)勢度,而在酸化低光條件下更適合金藻和藍(lán)藻的生長。

    長江口海域因?yàn)槭艿介L江沖淡水、臺(tái)灣暖流和江浙沿岸流等多個(gè)復(fù)雜水動(dòng)力條件的共同作用,形成了特有的海洋理化環(huán)境(馬永存, 2013),長江徑流每年向長江口水域匯入大量泥沙,疊加潮流的影響,在長江的咸淡水交匯區(qū)出現(xiàn)長江口渾濁帶(周華君, 1994),從而降低海表浮游植物受到的有效輻射光照強(qiáng)度,因此光照強(qiáng)度成為該海域浮游植物的典型限制因素(Domingues et al, 2014)。為了提高光子捕獲效率,浮游植物單位細(xì)胞的Chla含量通常在高光條件下有所降低(Agustí, 1991)。本調(diào)查站位位于長江口以南毗鄰海域的渾濁帶中,同樣受到復(fù)雜海洋理化環(huán)境影響。在本實(shí)驗(yàn)高光條件下浮游植物總Chla質(zhì)量濃度顯著高于低光條件的,這也印證了浮游植物生長的潛在光限制。影響浮游植物吸收光照的最主要因素是色素(梁橋等, 2020),不同浮游植物種群的色素組成不同會(huì)導(dǎo)致對太陽光吸收、利用的特性不同,例如硅藻和甲藻雖然僅有類胡蘿卜素類的輔助色素組成不同,但它們的吸收特性仍表現(xiàn)出差異(梁橋等, 2020)。這導(dǎo)致不同的光照強(qiáng)度下浮游植物種群組成也有所不同。在本研究中,硅藻的基因拷貝數(shù)在高光條件下最高,海鏈藻、角毛藻(Chaetoceros)和骨條藻(Skeletonema)的基因拷貝數(shù)均較對照組顯著升高。另外,我們還發(fā)現(xiàn)在高光條件下細(xì)胞粒徑大的浮游植物的占比增加,印證了大粒徑(>10 μm)浮游植物的光能利用效率更高的結(jié)論(Hashimoto et al, 2002)。

    酸化對浮游植物種群組成產(chǎn)生的影響與不同浮游植物種群的碳濃縮機(jī)制(CCMs)效率有關(guān)。在本研究中,酸化條件下金藻和聚球藻的基因拷貝數(shù)顯著升高,這是因?yàn)檫@些種群具有較低的CCM效率(Wei et al, 2021),pCO2升高可緩解其光合作用的潛在碳限制進(jìn)而促進(jìn)其生長(Low-décarie et al,2014)。而像硅藻這類具有高效CCMs的種群對pCO2上升的響應(yīng)并不顯著(Reinfelder, 2011; Wu et al,2014)。我們的實(shí)驗(yàn)結(jié)果表明在酸化以及交互的條件下,硅藻基因拷貝數(shù)的增長明顯低于高光條件。但已有研究表明,pCO2升高更有利于體積較大的硅藻的生長,例如Tortell和Trimborn對柔弱角毛藻(Chaetoceros debilis)和偽菱形藻(Pseudo-nitzschia)在不同pCO2下(16 Pa、40 Pa和101 Pa)的培養(yǎng)實(shí)驗(yàn)結(jié)果表明體積更大的中心綱硅藻柔弱角毛藻(Chaetoceros debilis)的生長速率受到pCO2增長的刺激更顯著(Tortell et al, 2008; Trimborn et al, 2013)。因此,我們觀測到酸化組中小型浮游植物生物量占比也顯著高于對照組。

    浮游植物細(xì)胞的沉降速率與浮游植物的碳輸出和碳沉降一定程度上存在相關(guān)性。本研究的結(jié)果表明酸化和高光對浮游植物群落的沉降速率具有相反的影響。與41 Pa相比,pCO2為101 Pa的酸化條件顯著升高了浮游植物的沉降速率,這主要是由于酸化引起的種群組成變化導(dǎo)致的。在酸化條件下浮游植物種群粒徑增大,根據(jù)斯托克斯公式(Stokes' Law)(Stokes, 1850),其沉降速率也會(huì)升高。此外,酸化也可導(dǎo)致浮游植物細(xì)胞內(nèi)生化組成發(fā)生改變(Bienfang et al, 1982),使其相對于海水的密度升高,另外單個(gè)浮游植物細(xì)胞的大小、細(xì)胞是否成鏈也會(huì)影響其沉降速率(Guo et al, 2016),本研究顯微鏡檢結(jié)果表明酸化處理組中海鏈藻較少成鏈,這也可能是導(dǎo)致沉降加快的因素。500~1200 μmol·photons·m?2·s?1的高光條件則降低了浮游植物的沉降速率。一些之前的研究也表明浮游植物的沉降速率在較低光照強(qiáng)度下會(huì)有所增加(Riebesell, 1989),Bienfang發(fā)現(xiàn)亞熱帶海域的浮游植物細(xì)胞在夜晚時(shí)的沉降速率是白天光照充足時(shí)的2倍(Bienfang, 1985)。Steele和Yentschd發(fā)現(xiàn)藻細(xì)胞處于活躍分裂狀態(tài)時(shí)的沉降速率比衰老狀態(tài)時(shí)減緩50%(Steele et al, 1960)。而較高的光照強(qiáng)度也提供給浮游植物足夠的能量來調(diào)整自身的沉降能力(Riebesell, 1989)。在酸化和高光的耦合作用下浮游植物的沉降速率處于升高的狀態(tài),但升高的程度較單酸化條件有所緩慢,2個(gè)條件對浮游植物的沉降速率表現(xiàn)出拮抗作用。十幾年來,長江口區(qū)域底層缺氧現(xiàn)象可引發(fā)的海洋相關(guān)生態(tài)災(zāi)難引起了越來越多的關(guān)注(Zhu et al, 2011)。而除水文動(dòng)力學(xué)條件影響,表層水體浮游植物水華后顆粒物會(huì)大量向深層沉降,有機(jī)物的降解過程的耗氧為該海域水體缺氧的一個(gè)重要產(chǎn)生因素。因此表層水體中浮游植物沉降速率的變化對該海域的碳輸出產(chǎn)生影響,也對其缺氧現(xiàn)象的產(chǎn)生有重要意義。我們的研究結(jié)果一定程度上佐證了酸化引起的浮游植物沉降速率升高可促進(jìn)長江口海域的底層缺氧的發(fā)生(Zhu et al, 2011)。

    4 結(jié)論

    2020年5月在長江口南毗鄰海域的渾濁帶區(qū)域?qū)Ω∮沃参锓N群現(xiàn)場采樣并進(jìn)行船基受控培養(yǎng),通過測定培養(yǎng)后浮游植物群落的不同化學(xué)參數(shù)并結(jié)合傳統(tǒng)顯微鏡檢和高通量測序兩種方法來鑒定浮游植物的種類、統(tǒng)計(jì)浮游植物的豐度,探究了海洋酸化和光照強(qiáng)度變化以及二者之間的交互作用對長江口以南毗鄰光限制區(qū)域的浮游植物生理指標(biāo)和群落組成變化的影響。

    我們的研究結(jié)果表明海水酸化程度與光照強(qiáng)度變化對浮游植物生長及種群組成產(chǎn)生了不同的影響,2種環(huán)境要素的變化可產(chǎn)生交互效應(yīng)。測定獲得的15個(gè)參數(shù)充分體現(xiàn)了浮游植物群落的生物量、化學(xué)元素組成和群落結(jié)構(gòu)組成等方面隨環(huán)境變化做出的響應(yīng)。500~1200 μmol·photons·m?2·s?1的高光更促進(jìn)浮游植物的生長及生物量積累,尤其是中心綱硅藻的豐度;pCO2為101 Pa的酸化更有利于聚球藻的生長,同時(shí)顯著提高了小型浮游植物Chla質(zhì)量濃度組成百分比并導(dǎo)致浮游植物種群粒徑增大;酸化和高光均降低了甲藻的優(yōu)勢度,而酸化條件下更適于金藻和藍(lán)藻的生長。浮游植物種群組成的改變將對該海域的碳沉降及缺氧區(qū)的形成產(chǎn)生深遠(yuǎn)影響。以上研究結(jié)果有助于我們進(jìn)一步理解和預(yù)測海洋酸化和光照強(qiáng)度變化對近岸低光渾濁帶的長江口南毗鄰海域的浮游植物的耦合作用及其生物地球化學(xué)效應(yīng)。已有的研究結(jié)果在時(shí)空上還存在較大局限性,今后的研究應(yīng)加強(qiáng)更多季節(jié)、更多典型站位研究,結(jié)合長期觀測數(shù)據(jù)進(jìn)一步拓展完善我們的趨勢性和機(jī)制性認(rèn)知。此外,有研究表明溫度升高和營養(yǎng)鹽濃度水平的變化(江志兵等, 2014; Liu et al, 2020)也可引起長江口海域浮游植物種群組成的顯著變化,對能量向上傳遞及漁業(yè)資源產(chǎn)生影響。今后的研究還要進(jìn)一步結(jié)合多重環(huán)境變化因素深入探究海洋環(huán)境快速變化對相應(yīng)海域的生態(tài)、生物地球化學(xué)效應(yīng)及其對氣候的反饋機(jī)制。

    致謝:本研究的數(shù)據(jù)及樣品采集得到國家自然科學(xué)基金委員會(huì)共享航次計(jì)劃項(xiàng)目(項(xiàng)目批準(zhǔn)號(hào):41676160和41606128) 的資助。 該航次(航次編號(hào):NORC2020-02)由“向陽紅18”號(hào)科考船實(shí)施,在此一并致謝。

    猜你喜歡
    硅藻高光拷貝數(shù)
    硅藻——海里的“生物寶石”
    大自然探索(2023年7期)2023-11-14 13:07:36
    硅藻實(shí)驗(yàn)室“打工記”
    見證彼此高光
    睿士(2023年8期)2023-08-23 13:07:06
    線粒體DNA拷貝數(shù)變異機(jī)制及疾病預(yù)測價(jià)值分析
    想做聯(lián)名爆款?超級飛俠、熊出沒、奧特曼等“高光”IP等你來聊
    玩具世界(2021年5期)2021-03-08 08:42:02
    硅藻指數(shù)在水環(huán)境監(jiān)測與評價(jià)中的應(yīng)用分析
    胎兒染色體組拷貝數(shù)變異與產(chǎn)前超聲異常的相關(guān)性分析
    硅藻泥涂料的現(xiàn)狀和發(fā)展趨勢
    上海建材(2018年1期)2018-04-18 12:15:16
    釋放出來自HDR的高光力量 Optoma奧圖碼UHZ880
    DNA序列拷貝數(shù)變化決定黃瓜性別

    免费看a级黄色片| 日本爱情动作片www.在线观看| 国语自产精品视频在线第100页| 欧美成人免费av一区二区三区| 中文字幕熟女人妻在线| 小说图片视频综合网站| 性欧美人与动物交配| 国产精品一区二区三区四区久久| 黄色一级大片看看| 久久九九热精品免费| 久久精品夜夜夜夜夜久久蜜豆| 成年免费大片在线观看| 亚洲欧美日韩东京热| 性插视频无遮挡在线免费观看| 悠悠久久av| 久久这里有精品视频免费| 久久九九热精品免费| 97人妻精品一区二区三区麻豆| 波野结衣二区三区在线| 一个人观看的视频www高清免费观看| 国产乱人视频| 国产v大片淫在线免费观看| 国产一区二区激情短视频| 亚洲欧美精品自产自拍| 毛片女人毛片| 亚洲精品日韩在线中文字幕 | 久久亚洲国产成人精品v| 26uuu在线亚洲综合色| 午夜激情福利司机影院| 久久人人精品亚洲av| av黄色大香蕉| 长腿黑丝高跟| 高清午夜精品一区二区三区 | 精品无人区乱码1区二区| 免费看a级黄色片| 国产伦理片在线播放av一区 | 日韩欧美三级三区| 欧美人与善性xxx| 美女 人体艺术 gogo| 成年女人永久免费观看视频| av.在线天堂| 欧美潮喷喷水| 午夜精品在线福利| 亚洲欧洲国产日韩| 九草在线视频观看| 欧美成人精品欧美一级黄| 一级毛片久久久久久久久女| 大型黄色视频在线免费观看| 国产探花在线观看一区二区| 三级毛片av免费| 国产精品人妻久久久影院| 99久久九九国产精品国产免费| 国产亚洲欧美98| 青青草视频在线视频观看| 久久欧美精品欧美久久欧美| 中文字幕av在线有码专区| 热99re8久久精品国产| a级毛片a级免费在线| 夜夜看夜夜爽夜夜摸| 国产色爽女视频免费观看| 蜜臀久久99精品久久宅男| 啦啦啦韩国在线观看视频| 一本久久中文字幕| 精品午夜福利在线看| 91久久精品电影网| 欧美另类亚洲清纯唯美| 亚洲无线在线观看| 三级国产精品欧美在线观看| 如何舔出高潮| 偷拍熟女少妇极品色| 亚洲电影在线观看av| 色尼玛亚洲综合影院| 国产精品一区二区三区四区久久| 人人妻人人澡人人爽人人夜夜 | 伊人久久精品亚洲午夜| 国产麻豆成人av免费视频| 校园人妻丝袜中文字幕| 亚洲av中文av极速乱| 精品久久久噜噜| 午夜激情欧美在线| 激情 狠狠 欧美| 插逼视频在线观看| 国产麻豆成人av免费视频| 亚洲av一区综合| 悠悠久久av| 久久久久免费精品人妻一区二区| 久久精品国产亚洲av天美| 看非洲黑人一级黄片| 亚洲性久久影院| 色视频www国产| 噜噜噜噜噜久久久久久91| 国产精品.久久久| 欧美3d第一页| 一级二级三级毛片免费看| 国产乱人视频| 日韩国内少妇激情av| 99热全是精品| 天美传媒精品一区二区| 国产女主播在线喷水免费视频网站 | 97超视频在线观看视频| 成人美女网站在线观看视频| 色综合站精品国产| 国产乱人视频| 久久久久国产网址| av黄色大香蕉| 国内精品宾馆在线| 美女国产视频在线观看| 成人特级黄色片久久久久久久| 老师上课跳d突然被开到最大视频| 好男人在线观看高清免费视频| 国产精品一区二区三区四区久久| 国内少妇人妻偷人精品xxx网站| 九九爱精品视频在线观看| 久久久久久久午夜电影| 亚洲国产精品sss在线观看| 久久久久国产网址| 欧美一级a爱片免费观看看| 嘟嘟电影网在线观看| 一区二区三区免费毛片| 又粗又爽又猛毛片免费看| 国产av麻豆久久久久久久| 久99久视频精品免费| or卡值多少钱| 日本一本二区三区精品| 最近中文字幕高清免费大全6| 久久99蜜桃精品久久| 白带黄色成豆腐渣| 蜜桃久久精品国产亚洲av| 成人特级av手机在线观看| 亚洲乱码一区二区免费版| 黄色配什么色好看| 女人被狂操c到高潮| 亚洲av熟女| 男女边吃奶边做爰视频| 亚洲国产精品sss在线观看| 你懂的网址亚洲精品在线观看 | 免费电影在线观看免费观看| 夜夜看夜夜爽夜夜摸| 国产精品乱码一区二三区的特点| 99久久人妻综合| 只有这里有精品99| 国产欧美日韩精品一区二区| 亚洲aⅴ乱码一区二区在线播放| 国产色爽女视频免费观看| 色哟哟·www| 亚洲精品日韩av片在线观看| 91av网一区二区| 国产高清激情床上av| 国产私拍福利视频在线观看| 一边亲一边摸免费视频| 亚洲av一区综合| 校园人妻丝袜中文字幕| 69av精品久久久久久| 免费在线观看成人毛片| 国产久久久一区二区三区| 亚洲成av人片在线播放无| 伊人久久精品亚洲午夜| 欧美+日韩+精品| 亚洲人成网站在线观看播放| 亚洲av中文字字幕乱码综合| 啦啦啦韩国在线观看视频| 内射极品少妇av片p| 亚洲av中文av极速乱| 国产av不卡久久| 床上黄色一级片| 国产视频首页在线观看| 少妇裸体淫交视频免费看高清| 欧美zozozo另类| 禁无遮挡网站| 国产成人a区在线观看| 亚洲成人av在线免费| 精华霜和精华液先用哪个| 亚洲欧洲国产日韩| 久久中文看片网| 黑人高潮一二区| 国产成人福利小说| 黄色视频,在线免费观看| 欧美性感艳星| 我要搜黄色片| 成人午夜高清在线视频| 午夜亚洲福利在线播放| 久久国产乱子免费精品| 午夜福利高清视频| 欧美又色又爽又黄视频| 午夜福利在线在线| 免费观看a级毛片全部| 麻豆成人午夜福利视频| 国产精品一区www在线观看| 在线免费观看不下载黄p国产| 久久精品国产亚洲网站| 在线免费十八禁| 国产亚洲欧美98| 国产久久久一区二区三区| 国产在视频线在精品| 十八禁国产超污无遮挡网站| 国产欧美日韩精品一区二区| 国产亚洲精品av在线| 特级一级黄色大片| 亚洲国产日韩欧美精品在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 不卡视频在线观看欧美| 免费av不卡在线播放| 99久久九九国产精品国产免费| 看十八女毛片水多多多| 亚洲成人久久性| 99热只有精品国产| 欧美性感艳星| 亚洲欧美日韩卡通动漫| 18禁裸乳无遮挡免费网站照片| 国产精品嫩草影院av在线观看| 久久精品夜色国产| 亚洲久久久久久中文字幕| 真实男女啪啪啪动态图| 免费av观看视频| 大香蕉久久网| 亚洲成人久久性| 哪个播放器可以免费观看大片| 26uuu在线亚洲综合色| 亚洲av.av天堂| 色综合站精品国产| 亚洲欧美日韩卡通动漫| 好男人在线观看高清免费视频| 久久草成人影院| 又爽又黄a免费视频| 可以在线观看毛片的网站| 国产伦理片在线播放av一区 | 色综合色国产| av.在线天堂| 亚洲精品亚洲一区二区| 最新中文字幕久久久久| 亚洲国产欧洲综合997久久,| 最近视频中文字幕2019在线8| 欧美成人a在线观看| 国产精品久久久久久精品电影小说 | 欧美另类亚洲清纯唯美| 免费观看精品视频网站| 不卡视频在线观看欧美| 亚洲欧美精品自产自拍| 黄片无遮挡物在线观看| 国产精品国产三级国产av玫瑰| 日韩强制内射视频| 久久久久网色| 精品人妻一区二区三区麻豆| 高清午夜精品一区二区三区 | 91久久精品国产一区二区成人| 最后的刺客免费高清国语| 免费看av在线观看网站| 国产在线精品亚洲第一网站| 成人午夜精彩视频在线观看| 国产成人freesex在线| 永久网站在线| 美女大奶头视频| 亚洲自偷自拍三级| 超碰av人人做人人爽久久| 亚洲婷婷狠狠爱综合网| 欧美xxxx性猛交bbbb| 直男gayav资源| 国产精品电影一区二区三区| 亚洲无线观看免费| 欧美区成人在线视频| 两个人的视频大全免费| 欧美bdsm另类| 一级毛片久久久久久久久女| 一本一本综合久久| 99热这里只有精品一区| 成人一区二区视频在线观看| 色综合色国产| 久久99精品国语久久久| 久久人人爽人人爽人人片va| 日日干狠狠操夜夜爽| 亚洲欧美日韩高清专用| 色尼玛亚洲综合影院| 国产伦理片在线播放av一区 | 国产伦精品一区二区三区视频9| a级毛片免费高清观看在线播放| 亚洲精品久久国产高清桃花| 又爽又黄无遮挡网站| 大又大粗又爽又黄少妇毛片口| 一个人看的www免费观看视频| 日韩人妻高清精品专区| 菩萨蛮人人尽说江南好唐韦庄 | 国内少妇人妻偷人精品xxx网站| 日韩制服骚丝袜av| 亚洲av一区综合| 国产美女午夜福利| 亚洲不卡免费看| 18+在线观看网站| av.在线天堂| 婷婷六月久久综合丁香| 青春草亚洲视频在线观看| 插阴视频在线观看视频| 亚洲在久久综合| 日本成人三级电影网站| 18禁黄网站禁片免费观看直播| 国产成人福利小说| 亚洲精品久久国产高清桃花| 中出人妻视频一区二区| 中文字幕av在线有码专区| 成人特级黄色片久久久久久久| 欧美日韩精品成人综合77777| 最近视频中文字幕2019在线8| 能在线免费观看的黄片| 国产高清三级在线| 看十八女毛片水多多多| 日本五十路高清| 国产成人福利小说| 精品欧美国产一区二区三| 亚洲精品乱码久久久v下载方式| 久久久久免费精品人妻一区二区| 欧美+亚洲+日韩+国产| 黑人高潮一二区| 亚洲婷婷狠狠爱综合网| 女的被弄到高潮叫床怎么办| 色5月婷婷丁香| 国内揄拍国产精品人妻在线| 91久久精品国产一区二区成人| 欧美3d第一页| 中文资源天堂在线| 日韩欧美 国产精品| 深夜精品福利| 男人的好看免费观看在线视频| 精品久久久久久成人av| 色综合色国产| 哪个播放器可以免费观看大片| 观看美女的网站| 非洲黑人性xxxx精品又粗又长| 欧美一级a爱片免费观看看| 国产精品一区二区性色av| 一区福利在线观看| 日本av手机在线免费观看| 免费看日本二区| 国产精品国产高清国产av| 成人国产麻豆网| 六月丁香七月| 最近最新中文字幕大全电影3| 日本成人三级电影网站| 国产成人91sexporn| 中国美女看黄片| 精品人妻一区二区三区麻豆| 国内精品美女久久久久久| 亚洲精品粉嫩美女一区| 狂野欧美白嫩少妇大欣赏| 蜜桃久久精品国产亚洲av| 搞女人的毛片| 91久久精品电影网| 一区福利在线观看| 午夜老司机福利剧场| 中文字幕久久专区| 91午夜精品亚洲一区二区三区| 青春草亚洲视频在线观看| 老司机福利观看| 成人高潮视频无遮挡免费网站| 亚洲美女视频黄频| 成人无遮挡网站| 十八禁国产超污无遮挡网站| 人妻夜夜爽99麻豆av| 国产一区二区激情短视频| 免费看av在线观看网站| 日韩,欧美,国产一区二区三区 | 欧美+亚洲+日韩+国产| 乱人视频在线观看| 国产69精品久久久久777片| 岛国毛片在线播放| 热99在线观看视频| 在线观看午夜福利视频| 国产成人精品一,二区 | 国产精华一区二区三区| 亚洲丝袜综合中文字幕| 国产三级在线视频| 亚洲av第一区精品v没综合| 在线播放国产精品三级| 成人特级黄色片久久久久久久| 床上黄色一级片| 色视频www国产| 白带黄色成豆腐渣| 波多野结衣高清作品| 狂野欧美白嫩少妇大欣赏| 伦精品一区二区三区| 精品一区二区免费观看| 亚洲欧美日韩高清专用| 91精品一卡2卡3卡4卡| 九九久久精品国产亚洲av麻豆| 麻豆成人午夜福利视频| 赤兔流量卡办理| 深夜a级毛片| 天天躁日日操中文字幕| 十八禁国产超污无遮挡网站| 在线观看av片永久免费下载| 国产精品精品国产色婷婷| 精品日产1卡2卡| 日本黄色视频三级网站网址| 久久久久久国产a免费观看| 国产一级毛片在线| 小说图片视频综合网站| 国产成人福利小说| 女人被狂操c到高潮| 国产老妇女一区| 亚洲熟妇中文字幕五十中出| 精品无人区乱码1区二区| 久久久久国产网址| 丰满人妻一区二区三区视频av| 最近2019中文字幕mv第一页| 久久精品国产亚洲av涩爱 | 成人鲁丝片一二三区免费| 国产精品久久久久久久电影| 九九热线精品视视频播放| 午夜福利成人在线免费观看| 精品不卡国产一区二区三区| 亚洲内射少妇av| 极品教师在线视频| 亚洲人成网站在线播放欧美日韩| 色综合站精品国产| 精品日产1卡2卡| 成年女人看的毛片在线观看| 亚洲精品国产成人久久av| 免费看av在线观看网站| 日本成人三级电影网站| 国产一区二区亚洲精品在线观看| 成人鲁丝片一二三区免费| 人人妻人人澡人人爽人人夜夜 | 日韩欧美精品免费久久| 免费观看精品视频网站| 国产欧美日韩精品一区二区| av天堂在线播放| 美女内射精品一级片tv| 精品国内亚洲2022精品成人| 国产午夜精品久久久久久一区二区三区| 亚洲五月天丁香| 日韩一区二区视频免费看| 欧美日韩在线观看h| 六月丁香七月| av专区在线播放| 18禁黄网站禁片免费观看直播| 好男人视频免费观看在线| 精品久久久久久久久久免费视频| 亚洲精品自拍成人| 午夜亚洲福利在线播放| 观看免费一级毛片| 亚洲精品亚洲一区二区| 一区二区三区四区激情视频 | 在线观看一区二区三区| 日本欧美国产在线视频| 91精品一卡2卡3卡4卡| 日韩制服骚丝袜av| 三级毛片av免费| 国产探花极品一区二区| 久久中文看片网| 成人鲁丝片一二三区免费| 国产成人精品久久久久久| 亚洲av第一区精品v没综合| 亚洲在线观看片| 美女大奶头视频| 性色avwww在线观看| 国产午夜精品一二区理论片| av天堂在线播放| 亚洲va在线va天堂va国产| 两性午夜刺激爽爽歪歪视频在线观看| 国产亚洲5aaaaa淫片| 男女视频在线观看网站免费| a级一级毛片免费在线观看| 美女大奶头视频| 久久国内精品自在自线图片| 尾随美女入室| 欧美日韩在线观看h| 日韩亚洲欧美综合| 久久久久久久久久成人| 99热只有精品国产| 成熟少妇高潮喷水视频| 成人毛片60女人毛片免费| 岛国在线免费视频观看| 可以在线观看毛片的网站| 亚洲成人中文字幕在线播放| 欧美色视频一区免费| 欧美日韩综合久久久久久| 岛国在线免费视频观看| 国产精品电影一区二区三区| 国产成人91sexporn| 欧美3d第一页| 大型黄色视频在线免费观看| 一级毛片我不卡| 成人高潮视频无遮挡免费网站| 色哟哟·www| 久久久久久久久久久丰满| 男人舔奶头视频| 日本一二三区视频观看| 国产午夜精品论理片| 亚洲av成人av| 高清在线视频一区二区三区 | 天堂网av新在线| 国内揄拍国产精品人妻在线| 日本欧美国产在线视频| 中文资源天堂在线| 国产一区二区激情短视频| 老女人水多毛片| 久久99热这里只有精品18| 国产单亲对白刺激| 日本免费一区二区三区高清不卡| av视频在线观看入口| 欧美变态另类bdsm刘玥| 亚洲国产欧美在线一区| 校园春色视频在线观看| 亚洲成人精品中文字幕电影| 亚洲最大成人手机在线| 波野结衣二区三区在线| 国产高清有码在线观看视频| 国产av麻豆久久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 少妇人妻一区二区三区视频| 精品一区二区三区视频在线| 婷婷亚洲欧美| 美女xxoo啪啪120秒动态图| 赤兔流量卡办理| 国产亚洲91精品色在线| 男女那种视频在线观看| 看非洲黑人一级黄片| 日本撒尿小便嘘嘘汇集6| 国产高潮美女av| 97热精品久久久久久| 国产伦一二天堂av在线观看| 变态另类丝袜制服| 高清在线视频一区二区三区 | 91狼人影院| 国产成人福利小说| 观看免费一级毛片| 18禁黄网站禁片免费观看直播| 桃色一区二区三区在线观看| 国产乱人视频| 国产美女午夜福利| 国产午夜福利久久久久久| 性色avwww在线观看| 国产精品一区二区三区四区免费观看| 国产精品爽爽va在线观看网站| 老女人水多毛片| 色噜噜av男人的天堂激情| 少妇熟女欧美另类| 午夜福利在线观看免费完整高清在 | 久久久久九九精品影院| av天堂中文字幕网| 禁无遮挡网站| 国产成人精品久久久久久| 男人和女人高潮做爰伦理| 欧美精品一区二区大全| 午夜爱爱视频在线播放| 成人特级黄色片久久久久久久| 最后的刺客免费高清国语| 国产免费男女视频| 国产精品女同一区二区软件| 丝袜喷水一区| 国产成人91sexporn| 国内久久婷婷六月综合欲色啪| 男人舔女人下体高潮全视频| 久久精品国产99精品国产亚洲性色| 日韩精品青青久久久久久| 一本一本综合久久| 欧美区成人在线视频| 狂野欧美激情性xxxx在线观看| 老师上课跳d突然被开到最大视频| 自拍偷自拍亚洲精品老妇| 欧美xxxx性猛交bbbb| 亚洲国产高清在线一区二区三| 一个人观看的视频www高清免费观看| 天天一区二区日本电影三级| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华液的使用体验 | 一卡2卡三卡四卡精品乱码亚洲| ponron亚洲| 高清毛片免费观看视频网站| 亚洲中文字幕日韩| 高清毛片免费观看视频网站| 中文亚洲av片在线观看爽| 欧美成人一区二区免费高清观看| .国产精品久久| 六月丁香七月| 日本成人三级电影网站| 春色校园在线视频观看| 亚洲人成网站在线播放欧美日韩| 久久午夜福利片| 亚洲,欧美,日韩| 亚洲av一区综合| 国产精品.久久久| 97热精品久久久久久| 成年女人看的毛片在线观看| 欧美色视频一区免费| 精品一区二区三区人妻视频| 18禁在线无遮挡免费观看视频| av福利片在线观看| 一本久久精品| 免费观看在线日韩| 国产亚洲av嫩草精品影院| 午夜福利在线观看免费完整高清在 | 久久久色成人| 国产一级毛片七仙女欲春2| 日日干狠狠操夜夜爽| 久久鲁丝午夜福利片| 综合色丁香网| 老熟妇乱子伦视频在线观看| 卡戴珊不雅视频在线播放| 国产精品一区二区三区四区免费观看| 亚洲国产精品成人综合色| 日本熟妇午夜| 日本免费a在线| .国产精品久久| 高清毛片免费观看视频网站| 亚洲电影在线观看av| 午夜精品在线福利| 国产女主播在线喷水免费视频网站 | 蜜桃亚洲精品一区二区三区| 精华霜和精华液先用哪个| 日本在线视频免费播放| 免费电影在线观看免费观看| 亚洲不卡免费看| 亚洲精品色激情综合| 精品久久久久久久久久免费视频| 九色成人免费人妻av| 国产午夜精品一二区理论片|