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

    南黃海首次暴發(fā)“巨囊”生態(tài)型球形棕囊藻藻華*

    2022-09-21 02:37:12張清春王錦秀宋敏杰王云峰孔凡洲耿慧霞于仁成
    海洋與湖沼 2022年5期
    關(guān)鍵詞:囊體生態(tài)型黃海

    張清春 向 玲 王錦秀 劉 揚 宋敏杰 王云峰 孔凡洲 耿慧霞 顏 天 于仁成

    南黃海首次暴發(fā)“巨囊”生態(tài)型球形棕囊藻藻華*

    張清春1, 2, 4向 玲1, 3王錦秀1, 3劉 揚1, 5宋敏杰1, 3王云峰1, 2, 4孔凡洲1, 2, 4耿慧霞1, 2, 4顏 天1, 2, 4于仁成1, 2, 3, 4①

    (1. 中國科學(xué)院海洋生態(tài)與環(huán)境科學(xué)重點實驗室(中國科學(xué)院海洋研究所) 山東青島 266071; 2. 青島海洋科學(xué)與技術(shù)國家試點實驗室海洋生態(tài)與環(huán)境科學(xué)功能實驗室 山東青島 266071; 3. 中國科學(xué)院大學(xué) 北京 100049; 4. 中國科學(xué)院海洋大科學(xué)中心 山東青島 266071; 5. 青島科技大學(xué) 山東青島 266045)

    球形棕囊藻()常在熱帶至溫帶近海海域形成有害藻華, 該種具有種內(nèi)遺傳多態(tài)性, 且有囊體大小、特征色素組成等性狀分化。在我國南海形成有害藻華的球形棕囊藻可形成“巨型囊體”, 以19′-丁酰氧基巖藻黃素(But-fuco)為特征色素, 是一種獨特的“巨囊”生態(tài)型。2021年11月底至12月初, 南黃海青島沿岸暴發(fā)大規(guī)模球形棕囊藻藻華。藻華發(fā)生期間, 對青島沿岸3個站位的水文和化學(xué)要素進行了觀測, 分析了球形棕囊藻囊體的數(shù)量、直徑和色素組成, 并應(yīng)用一種高分辨率種下分子標記——葉綠體S-27基因間隔區(qū)分析了其遺傳特征。結(jié)果表明, 藻華發(fā)生期間青島沿岸表層海水溫度較低(12~14 °C), 海域營養(yǎng)鹽組成具有高溶解有機氮、低溶解無機氮的特征; 球形棕囊藻囊體豐度超過20個/L, 最大囊體直徑為18 mm, 以But-fuco為特征色素,S-27序列分析表明其與南?!熬弈摇鄙鷳B(tài)型球形棕囊藻具有相同的遺傳特征。南黃海首次暴發(fā)的球形棕囊藻藻華是由“巨囊”生態(tài)型形成, 該藻華可能對海洋生態(tài)系統(tǒng)、水產(chǎn)養(yǎng)殖業(yè)發(fā)展和核電設(shè)施運行等構(gòu)成威脅, 亟待開展藻華成因與監(jiān)測預(yù)警對策研究。

    球形棕囊藻; “巨囊”生態(tài)型; 有害藻華; 南黃海

    2021年11月底12月初, 南黃海青島沿岸突發(fā)球形棕囊藻()形成的有害藻華(harmful algal bloom, HAB)。大量球形棕囊藻囊體聚集于青島市黃島至膠南沿岸海域, 這是南黃海海域首次暴發(fā)球形棕囊藻藻華。球形棕囊藻廣泛分布于熱帶至溫帶海洋, 能產(chǎn)生二甲基丙酸鹽, 對全球氣候變化具有重要影響(Charlson, 1987; Wassmann, 2005; Whipple, 2005)。球形棕囊藻具有復(fù)雜的異形生活史, 含有游離單細胞和囊體群體兩種迥異的細胞形態(tài)(齊雨藻等, 2001), 其中囊體是由黏多糖包裹成千上萬個藻細胞形成(Van Rijssel, 1997)。球形棕囊藻常以囊體形式形成藻華, 在全球許多近海海域均有該藻藻華記錄(Qi, 2004; Smith, 2014b)。在球形棕囊藻藻華衰退過程中, 死亡藻體分解消耗大量溶解氧, 囊體分解形成泡沫物質(zhì), 影響海洋生態(tài)系統(tǒng)和旅游景觀, 亦會產(chǎn)生溶血毒素等有毒有害化學(xué)物質(zhì), 危害海水養(yǎng)殖業(yè)發(fā)展; 藻華期間形成的大量囊體易堵塞冷卻水過濾網(wǎng), 影響核電等沿海大型設(shè)施的安全運行(齊雨藻等, 2001; Yu, 2017)。

    1997年, 球形棕囊藻首次在我國東南沿海形成有害藻華(陳菊芳等, 1999), 之后在廣東、海南和廣西近海海域連年暴發(fā), 成為我國東南沿海最為突出的藻華災(zāi)害(李亞男等, 2012; 沈萍萍等, 2018)。在渤海也多次暴發(fā)該藻藻華(曲凌云等, 2008; 屠建波等, 2011)。在本次青島沿岸藻華之前, 未見有東海和南黃海發(fā)生球形棕囊藻藻華的報道。

    球形棕囊藻具有囊體大小和特征色素組成等性狀分化。包括我國近海在內(nèi)的東亞至南亞一帶沿岸海域, 球形棕囊藻囊體常達到厘米級, 稱為“巨型囊體”(Giant colony) (齊雨藻等, 2001; Smith, 2014b)。與之相比, 全球其他海域的球形棕囊藻囊體多為毫米級, 通常不超過3 mm (Rousseau, 1990; Smith, 2014b)。除了囊體大小, 該藻的特征色素組成差異明顯(Vaulot, 1994; Zapata, 2004)。巖藻黃素(fucoxanthin, Fuco)的兩種衍生物19′-己酰氧基巖藻黃素(19′-hexanoyloxyfucoxanthin, Hex-fuco)和19′-丁酰氧基巖藻黃素(19′-butanoyloxyfucoxanthin, But-fuco)為棕囊藻的特征色素(Zapata, 2000, 2004)。具體到球形棕囊藻, 部分株系含有上述兩種特征色素, 但也有株系僅含有其中一種特征色素或兩種特征色素都不含有(Vaulot, 1994; Antajan, 2004; Zapata, 2004; Seoane, 2009; 王錦秀等, 2019)。基于核糖體18S rDNA、28S rDNA和ITS等常見分子標記的研究發(fā)現(xiàn), 球形棕囊藻具有種下遺傳多樣性(Lange, 2002; Chen, 2003; Medlin, 2007), 但是現(xiàn)有分子標記過于保守或高變, 與性狀特征分化或地理分布之間無法關(guān)聯(lián)。近年的研究發(fā)現(xiàn), 葉綠體和線粒體基因組中的高變區(qū)能較好地分辨球形棕囊藻不同地理株系, 與囊體大小、特征色素組成等也有很好的對應(yīng)關(guān)系, 可以用作球形棕囊藻種下遺傳分化的分子標記(Song, 2020, 2021; Zhang, 2021; Niu, 2022)。

    基于球形棕囊藻葉綠體S-27基因間隔區(qū)和線粒體8基因高分辨率分子標記, 可以將全球不同海域20余株球形棕囊藻分為4個基因型, 其中Ⅰ型和Ⅱ型廣布于世界各海域, 而Ⅲ型株系僅分布于大西洋北海海域, Ⅳ型僅分布于我國南海海域。前3個基因型的球形棕囊藻僅形成直徑≤3 mm小囊體, 含有Hex-fuco和But-fuco兩種特征色素, 但Ⅳ型可形成“巨型囊體”, 僅含But-fuco一種特征色素(王錦秀等, 2019; Zhang, 2021, 2022; Niu, 2022)。綜合其遺傳特征、囊體大小和特征色素組成, Ⅳ型球形棕囊藻被定義為一個獨特的“巨囊”生態(tài)型( Zhang, 2022)。2015~2021年間, 南海多起球形棕囊藻藻華期間采集的囊體均以But-fuco為特征色素, 且都屬于Ⅳ型, 表明“巨囊”生態(tài)型引發(fā)了南海球形棕囊藻藻華(Wang, 2022; Zhang, 2022)。

    目前對引發(fā)青島沿岸藻華的球形棕囊藻色素組成和遺傳特征都不了解。在球形棕囊藻藻華發(fā)生期間,本文在青島沿岸采集了球形棕囊藻囊體, 對囊體大小、色素組成以及遺傳信息進行了分析, 初步研判了此次球形棕囊藻藻華原因種的生物學(xué)特征, 以期為解析南黃海球形棕囊藻發(fā)生機制和有效防控藻華危害提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 樣品采集

    2021年12月6日, 分別在青島市黃島區(qū)沿岸的唐島灣(TDW)、星光島(XGD)和古鎮(zhèn)口(GZK)海域設(shè)置采樣點(圖1)。在每個采樣點采集表層海水1~2 L, 倒入白瓷盤中, 以刻度尺測量囊體直徑, 最小可測量囊體直徑約0.5 mm。每個采樣點挑取單個或多個囊體, 測量并記錄直徑后, 以滅菌海水清洗去除囊體表面附著物, 分別在低負壓下過濾到孔徑為0.4 μm聚碳酸酯膜(Millipore, 美國) 和孔徑約為0.68 μm GF/F玻璃纖維濾膜(Whatman, 英國) 上, 迅速放入液氮中保存, 帶回實驗室于-80 °C冰箱中存放, 用于囊體分子和色素特征分析。取80 mL海水樣品用經(jīng)450 °C灼燒的GF/F玻璃纖維濾膜過濾, 裝入100 mL聚碳酸酯瓶中, 帶回實驗室存放于-40 °C冰箱, 用于營養(yǎng)鹽分析?,F(xiàn)場海水溫度、鹽度、pH以及溶解氧(DO)由梅特勒多參數(shù)測試儀(Mettler Toledo, Columbus, 美國)測定。

    1.2 營養(yǎng)鹽樣品分析

    1.3 色素樣品分析

    棕囊藻囊體色素提取和分析參照Zapata等(2000)和Kong等(2012)的方法。將囊體色素濾膜剪碎后, 加入1 400 μL 95%甲醇水溶液進行色素提取, 同時加入100 μL 8′-apo-β,ψ-carotaldehyde (750 μg/L)作為內(nèi)標, 充分攪拌后冰浴超聲處理5 min。以孔徑0.22 μm有機相針頭濾器對提取液進行過濾, 吸取800 μL濾液至1.5 mL棕色樣品瓶中, 加入160 μL超純水, 震蕩混勻, 立即進行高效液相色譜(HPLC)色素分析。

    圖1 我國渤海和黃海(a)和青島沿岸采樣站位地圖(b)

    注: TDW: 唐島灣站; XGD: 星光島站; GZK: 古鎮(zhèn)口站

    1.4 DNA提取、PCR擴增和測序

    采用十六烷基三甲基溴化銨(Hexadecyl trimethyl ammonium Bromide, CTAB)提取法提取囊體過濾膜樣基因組DNA (Winnepenninckx, 1993), 提取的DNA溶解于TE溶液中, 保存于-20 °C用于PCR擴增。利用前期研發(fā)的高分辨分子標志擴增球形棕囊藻葉綠體S-27基因間隔區(qū), 引物為S-27-F (5′-GAAGCATACGCACCAAGC-3′) 和S-27-3R (5′-ATTGTTGAGTTTACGACGGATA-3′) (Zhang, 2021), PCR體系和擴增條件分別參照Zhang等(2021)。用1%瓊脂糖凝膠電泳檢測PCR擴增產(chǎn)物, 然后送測序公司進行雙端測序, 測序獲得序列用Invitrogen 11.0 ContigExpress (ThermoFisher, 美國)進行組裝。

    1.5 系統(tǒng)進化分析

    從GenBank下載全球不同海域的20株球形棕囊藻葉綠體S-27基因間隔區(qū)序列(表1), 加上本研究獲得的6條序列, 并以南極棕囊藻() CCMP1374和約氏棕囊藻() CCMP2496為外群進行系統(tǒng)進化分析。使用ClustalW對序列進行比對; 使用MEGA10.0程序中Find Best DNA/Protein Models (ML)功能選擇最適合替換模型T92, 構(gòu)建最大似然進化樹(maximum- likelihood tree), bootstrap值為1 000 (Kumar, 2018);使用Modelfinder程序選擇貝葉斯樹(Bayesian tree)最適合替換模型為HKY+F+G4 (Kalyaanamoorthy, 2017), 使用MrBayes3.2.6構(gòu)建貝葉斯樹, 四條MCMC鏈運行2 000 000代, 每1 000代進行一次取樣(Ronquist, 2012)。

    2 研究結(jié)果

    2.1 環(huán)境參數(shù)

    2.2 球形棕囊藻囊體數(shù)量與大小

    在唐島灣采集1 L表層海水中, 觀測到直徑介于0.5~18 mm的球形棕囊藻囊體27個; 星光島采集2 L表層海水中直徑介于1~12 mm的球形棕囊藻囊體37個; 古鎮(zhèn)口1 L表層海水中直徑介于4~8 mm球形棕囊藻囊體3個(圖2)。將3個站位的囊體信息綜合分析表明, 現(xiàn)場海水中直徑在5~10 mm的囊體占優(yōu)勢(55%), 其次是直徑1~5 mm的囊體(30%), 直徑大于10 mm囊體數(shù)量較低, 最大囊體直徑可達18 mm (圖3)。

    表1 本研究中球形棕囊藻和南極棕囊藻的信息

    Tab.1 Information of Phaeocystis globosa and P. Antarctica strains for this study

    注: /表示株系不成囊或偶爾形成直徑≤3 mm小囊體; —表示具體信息不詳

    表2 2021年12月青島沿岸三個調(diào)查站位的環(huán)境因子

    Tab.2 Environmental parameters measured at the three sampling sites in the coastal water of Qingdao in December 2021

    2.3 球形棕囊藻色素組成

    藻華現(xiàn)場采集9個直徑介于4~18 mm的球形棕囊藻囊體進行色素分析, 結(jié)果如圖4所示。檢測到的主要色素包括脫植基葉綠素(chlide)、甲基脫植基葉綠素(Me-chlide)、19′-丁酰氧基巖藻黃素(But-fuco)、巖藻黃素(fuco)、硅甲藻黃素(diad)、硅藻黃素(diat)、葉綠素(chl)、α-胡蘿卜素(βε-Car)和β-胡蘿卜素(ββ-Car)。所有囊體樣品中都未檢測到Hex-fuco, 僅檢測到But-fuco。結(jié)果表明形成藻華的球形棕囊藻以But-fuco為特征色素。

    參照王錦秀等(2019)前期建立的囊體直徑-囊體細胞數(shù)方程=0.00762.120 4, 根據(jù)囊體直徑換算得出每個囊體樣品中含有藻細胞數(shù), 計算出囊體中單細胞主要色素含量(表3)。fuco和chl是球形棕囊藻囊體細胞中的主要色素組分, 含量分別為0.120~0.279 pg/cell和0.076~0.303 pg/cell; 其次為diad, 含量為0.009~0.092 pg/cell; 細胞特征色素But-fuco含量為0.005~0.019 pg/cell。

    圖2 球形棕囊藻藻華囊體的現(xiàn)場照片(a, b)和顯微照片(c)

    圖3 2021年12月青島沿岸球形棕囊藻藻華期間囊體直徑分布狀況

    2.4 球形棕囊藻種下多態(tài)性分析

    采自三個站位的6個球形棕囊藻囊體的葉綠體S-27基因間隔區(qū)序列完全一致。根據(jù)構(gòu)建的系統(tǒng)進化樹可以看出, 所有球形棕囊藻聚集在一起, 遠離南極棕囊藻和約氏棕囊藻。球形棕囊藻可分為4個獨立簇, 其中第1和2簇聚集在一起, 第3和4簇聚集一起。第1簇包括中國南海北部灣株系MEL43、MEL44、MEL68和MEL69, 南海珠江口株系MEL45、東海泉州海域株系MEL47, 越南海域株系MEL78、泰國灣株系CCMP1524以及大西洋加勒比海, 墨西哥灣和英吉利海峽株系CCMP628、CCMP629和RCC055, 這11個株系內(nèi)部還有遺傳分化。第2簇包括分離自東海霞浦海域西株系MEL71以及大西洋北美海域株系CCMP2754。第3簇包括分離自大西洋北海海域株系RCC678、RCC736、RCC1736和Pg-G(A)。第4簇包括分離自南海北部灣和大亞灣株系MEL42、MEL70和TIO393, 以及本研究中從青島沿岸采集的6個囊體(圖5)。

    3 討論

    3.1 青島沿岸球形棕囊藻藻華是由“巨囊”生態(tài)型形成

    本研究首次報道了南黃海海域暴發(fā)的球形棕囊藻藻華。球形棕囊藻主要分布于熱帶至溫帶大洋和近海海域, 可引發(fā)有害藻華。20世紀80~90年代, 球形棕囊藻藻華主要出現(xiàn)在北大西洋的北海及其鄰近海域(Rousseau, 1990; Lancelot, 2009)。20世紀末, 太平洋西部的東亞和東南亞海域頻發(fā)球形棕囊藻藻華, 位于印度洋的阿拉伯海等也有零星出現(xiàn)(Madhupratap, 2000; Qi, 2004; Hai, 2010; Smith, 2014b)。1997年10月, 中國東南沿海首次暴發(fā)大規(guī)模球形棕囊藻藻華, 隨后, 廣東、海南和廣西沿海亦頻頻發(fā)生。2004年6月至2006年10月, 渤海沿海多次發(fā)生球形棕囊藻藻華(曲凌云等, 2008; 屠建波等, 2011)。在本次青島沿岸發(fā)生球形棕囊藻藻華之前, 未曾有南黃海發(fā)生球形棕囊藻藻華的報道。

    圖4 球形棕囊藻囊體色素色譜圖(440 nm)

    表3 2021年12月青島沿岸三個調(diào)查站位球形棕囊藻囊體細胞中的主要色素含量*

    Tab.3 Content of pigments in colony cells collected from the three sampling sites in the coastal water of Qingdao in December 2021

    注:*因缺少childe標準品, 未計算細胞內(nèi)childe色素含量

    圖5 基于rbcS-rpl27序列建立的球形棕囊藻系統(tǒng)進化樹

    Figure 5 Phylogenetic tree established forbased on theS-27 intergenic spacer

    注: 圖中數(shù)值斜線左邊是最大似然bootstrap值, 右邊是貝葉斯分析值, 圖中加粗的為本研究中從青島沿岸采集的6個球形棕囊藻囊體

    前期研究發(fā)現(xiàn), 球形棕囊藻具有種下遺傳多態(tài)性, 且存在囊體大小和特征色素等性狀分化(Vaulot, 1994; Zapata, 2004; Medlin, 2007; Smith, 2014b)。基于核糖體18S rDNA、28S rDNA、ITS和L等常見分子標記研究表明, 球形棕囊藻具有復(fù)雜的種下遺傳多態(tài)性(Medlin, 2007)。在大西洋北海海域, 球形棕囊藻囊體一般小于3 mm, 含Hex-fuco和But-fuco兩種特征色素。而在東南亞沿海, 球形棕囊藻藻華中常出現(xiàn)厘米級大小的“巨型囊體”, 且僅含有But-fuco一種特征色素。隨著研究深入, 在我國近海也陸續(xù)發(fā)現(xiàn)了僅形成小囊體、以Hex-fuco和But-fuco為特征色素的球形棕囊藻株系(王錦秀等, 2019)。盡管球形棕囊藻囊體形成受到光照、溫度、營養(yǎng)鹽、硅藻, 以及浮游動物攝食等多種因素影響(Tang, 2003; Sazhin, 2007; Wang, 2010, 2011), 但上述因素并非決定球形棕囊藻大小的根本因素。棕囊藻特征色素含量會受到光照等環(huán)境因素影響(Buma, 1991; van Leeuwe, 1998, 2014; Schlüter, 2000; DiTullio, 2007), 但其特征色素組成與環(huán)境因素無關(guān)。因此, 球形棕囊藻囊體大小和特征色素等性狀應(yīng)當具有其遺傳基礎(chǔ)。早期研究中采用的分子標記因過于保守或高變, 無法合理解釋球形棕囊藻的生物地理分布狀況, 與其性狀特征分化也沒有清晰的關(guān)聯(lián)。

    生態(tài)型(ecotype)是指同一物種的不同種群長期生活在不同環(huán)境中, 由于環(huán)境選擇和隔離產(chǎn)生趨異適應(yīng), 成為遺傳上有差異、適應(yīng)不同生態(tài)環(huán)境的類群, 相當于分類學(xué)上的亞種(Turesson, 1922; Begon, 2005)?;谌~綠體S-27基因間隔區(qū)和線粒體8基因分子標記研究都表明, 分離于南海的部分球形棕囊藻藻株能形成厘米級“巨型囊體”, 以But-fuco作為特征色素, 且其在基因型屬于Ⅳ型。在南海海域2015~2021年期間的調(diào)查發(fā)現(xiàn), 該海域存在包括Ⅰ、Ⅱ和Ⅳ型之內(nèi)多種基因型球形棕囊藻游離細胞, 但藻華期間采集的球形棕囊藻囊體都屬于Ⅳ型(Zhang, 2021, 2022; Niu, 2022)。室內(nèi)模擬實驗也發(fā)現(xiàn), Ⅳ型球形棕囊藻比Ⅰ和Ⅱ型具有更高的溶血活性和生物毒性。因此, “巨型囊體”球形棕囊藻與其他球形棕囊藻相比具有明顯的性狀和遺傳差異, 應(yīng)屬于一類獨特的生態(tài)型——“巨囊”生態(tài)型(Zhang, 2022)。本研究發(fā)現(xiàn), 在青島沿岸形成藻華的球形棕囊藻可以形成厘米級“巨型囊體”, 囊體葉綠體S-27基因間隔區(qū)與“巨囊”生態(tài)型球形棕囊藻完全一致, 均屬于Ⅳ型, 且以But-fuco作為特征色素。上述結(jié)果表明此次青島沿岸球形棕囊藻藻華是由“巨囊”生態(tài)型形成。

    3.2 南黃海球形棕囊藻藻華的監(jiān)測和研究展望

    球形棕囊藻藻華暴發(fā)具有規(guī)模大、時間長等特征。如2016~2017年北部灣球形棕囊藻藻華影響海域覆蓋了從廣西防城港至廣東雷州半島的近岸海域, 持續(xù)時間從2016年12月至次年3月(賀成等, 2019); 1997年我國東南沿海首次暴發(fā)球形棕囊藻藻華, 面積超過3 000 km2, 持續(xù)6個月(陳菊芳等, 1999)。此外, 球形棕囊藻藻華往往會在暴發(fā)海域多年持續(xù)發(fā)生, 很難消除。20世紀80~90年代, 球形棕囊藻藻華在北大西洋的北海及鄰近海域幾乎連年出現(xiàn), 1997年以來, 我國南海海域也連續(xù)多年暴發(fā), 2011年至今, 北部灣幾乎每年冬季都會發(fā)生(沈萍萍等, 2021)。南黃海海域的球形棕囊藻藻華是否會連年出現(xiàn), 值得密切關(guān)注。

    球形棕囊藻藻華是典型的有害藻華, 對海洋生態(tài)系統(tǒng)、水產(chǎn)養(yǎng)殖業(yè)以及沿岸大型設(shè)施運行等有不利影響(沈萍萍等, 2021)。南黃海球形棕囊藻藻華的潛在危害同樣值得密切關(guān)注。球形棕囊藻藻華暴發(fā)時會產(chǎn)生二甲基硫化物和硫丙酸等有毒物質(zhì), 常造成海洋生物死亡。藻華消亡時, 高生物量囊體細胞分解消耗大量氧氣, 造成水體缺氧, 威脅海洋生物生存(沈萍萍等, 2018)。1997年我國東南沿海暴發(fā)的大規(guī)模藻華造成養(yǎng)殖生物大量死亡, 僅廣東拓林灣的漁業(yè)直接損失超過6 000萬元(陳菊芳等, 1999)。前人研究發(fā)現(xiàn), 球形棕囊藻南海株系可產(chǎn)生溶血毒素, 是造成大量海洋生物死亡的原因之一(何家菀等, 1999)。本文近期研究也發(fā)現(xiàn), “巨囊”生態(tài)型球形棕囊藻(Ⅳ型)比Ⅰ和Ⅱ型具有更高的溶血活性, 對輪蟲和鹵蟲等浮游生物的生物毒性更強(Zhang, 2022)。南黃海的“巨囊”生態(tài)型球形棕囊藻有可能對該海域的水產(chǎn)養(yǎng)殖活動造成損害。球形棕囊藻具有膠質(zhì)囊體, 容易堵塞冷卻系統(tǒng)過濾網(wǎng), 威脅沿海地區(qū)核電和熱點設(shè)施的安全運行,南海北部灣海域的球形棕囊藻藻華對核電運行的影響是備受關(guān)注的熱點問題(曹西華等, 2017)。近年來, 我國沿海大力發(fā)展核電和火電工程, 南黃海沿岸已建立和正在建立多個核電和火電廠, 依賴大量海水來維持冷源系統(tǒng)正常運行。南黃海海域“巨囊“生態(tài)型球形棕囊藻藻華對核電和火電冷源系統(tǒng)安全運行的影響值得高度關(guān)注。

    球形棕囊藻藻華在渤海和南海都有多次形成大規(guī)模藻華的記錄, 但在黃、東海海域很少發(fā)生, 其原因有待深入研究。本研究發(fā)現(xiàn), 此次南黃海的球形棕囊藻藻華是由“巨囊”生態(tài)型形成。以往渤海和北黃海球形棕囊藻藻華期間, 也有“巨型囊體”出現(xiàn)的記錄(李偉才等, 2006; 曲凌云等, 2008), 但由于缺乏現(xiàn)場樣品和分離的球形棕囊藻培養(yǎng)系, 無法進行進一步驗證?;?021年對渤海、黃海和東海浮游植物群樣品的分析, 我們發(fā)現(xiàn)渤海也存在Ⅳ型球形棕囊藻, 而黃海、東海海域的球形棕囊藻則以其他基因型為主。這表明渤海海域以往的球形棕囊藻藻華極有可能也是由“巨囊”生態(tài)型形成。但是, 此次藻華過程中的“巨囊”生態(tài)型球形棕囊藻從何而來, 我國南方、北方海域的“巨囊”生態(tài)型球形棕囊藻是否存在種群隔離和生理生態(tài)特征的差異等, 目前并不清楚。其次, 這次南黃海球形棕囊藻藻華也具有一些獨有的特征: 藻華期海水溫度僅為12~14 °C, 遠低于南海海域藻華發(fā)生時的海水溫度(水溫20 °C左右) (徐寧等, 2003; 賀成等, 2019)。此外, 以往研究認為高濃度硝態(tài)氮是球形棕囊藻大量繁殖和形成藻華的重要物質(zhì)基礎(chǔ)(Lancelot, 2007; Wang, 2021), 但在南黃海球形棕囊藻藻華期間, 海水中DON濃度(21~23 μmol/L)遠高于DIN濃度(5.0~9.5 μmol/L), 高濃度DON是否有利于此次球形棕囊藻藻華形成仍有待闡明。鑒于球形棕囊藻藻華對海洋生態(tài)系統(tǒng)、水產(chǎn)養(yǎng)殖業(yè)、沿岸核電和火電等設(shè)施的潛在危險, 有必要開展針對性的研究, 深入揭示我國北方海域球形棕囊藻藻華的發(fā)生機制, 在此基礎(chǔ)上進一步強化對北方海域球形棕囊藻藻華的監(jiān)測預(yù)警。

    本次南黃海球形棕囊藻藻華已被證實是由“巨囊”生態(tài)型形成。在海水中, 球形棕囊藻可以游離單細胞和囊狀群體兩種形態(tài)存在(齊雨藻等, 2001), 囊體出現(xiàn)常被視為藻華形成的標志。在藻華形成前, 球形棕囊藻以游離單細胞形式存活于海水中。球形棕囊藻游離單細胞大小僅為3~8 μm, 通過傳統(tǒng)的顯微鏡觀察手段難以實現(xiàn)對游離細胞的鑒定和監(jiān)測?!熬弈摇鄙鷳B(tài)型球形棕囊藻以But-fuco作為特征色素, 在北部灣藻華的研究中, Wang等(2022)將But-fuco含量作為診斷色素, 成功應(yīng)用于球形棕囊藻藻華的過程研究和早期診斷。根據(jù)“巨囊”生態(tài)型球形棕囊藻的遺傳序列信息, Niu等(2022)開發(fā)了以線粒體8基因為靶區(qū)的特異性qPCR定量方法, 檢測限可達到30 cells/L, 成功應(yīng)用于北部灣海域球形棕囊藻藻華監(jiān)測和研究中。對青島沿岸采集的球形棕囊藻線粒體8基因擴增測序結(jié)果表明, 其序列信息與南?!熬弈摇鄙鷳B(tài)型序列一致, 這表明前期研發(fā)的“巨囊”生態(tài)型球形棕囊藻特異性qPCR方法可以用于北方海域球形棕囊藻的檢測。這些研究手段為進一步開展北方海域球形棕囊藻的監(jiān)測和研究工作奠定了基礎(chǔ)。

    4 結(jié)論

    2021年11月底至12月初, 青島沿岸暴發(fā)球形棕囊藻藻華, 藻華中產(chǎn)生大量囊體主要集聚于青島黃島至膠南沿岸。藻華發(fā)生時, 青島沿岸表層海水溫度低, 僅為12~14 °C, 呈現(xiàn)高DON、低DIN營養(yǎng)鹽組成特征。海水中囊體豐度可超過20個/L, 最大直徑超過1 cm。藻華中形成的囊體與南?!熬弈摇鄙鷳B(tài)型球形棕囊藻一樣, 均以But-fuco為特征色素; 葉綠體S-27序列完全一致, 與南?!熬弈摇鄙鷳B(tài)型球形棕囊藻同屬于基因Ⅳ型?!熬弈摇鄙鷳B(tài)型球形棕囊藻藻華可能會對南黃海生態(tài)系統(tǒng)、水產(chǎn)養(yǎng)殖業(yè)和沿岸核電、熱電冷卻系統(tǒng)構(gòu)成威脅, 亟需加強我國北方海域球形棕囊藻藻華形成機制研究, 進一步提高監(jiān)測預(yù)警能力, 防控藻華危害。

    致謝 感謝參與現(xiàn)場調(diào)查和采樣全體科研人員, 以及煙臺大學(xué)江濤教授提供部分現(xiàn)場調(diào)查信息。

    王錦秀, 孔凡洲, 陳振帆, 等, 2019. 六株球形棕囊藻的色素組成特征研究[J]. 海洋與湖沼, 50(3): 611-620.

    曲凌云, 呂頌輝, 高春蕾, 等, 2008. 棕囊藻渤海株核糖體18S rDNA和ITS基因結(jié)構(gòu)序列分析[J]. 海洋科學(xué)進展, 26(2): 200-206.

    齊雨藻, 沈萍萍, 王艷, 2001. 棕囊藻屬()的分類與生活史(綜述)[J]. 熱帶亞熱帶植物學(xué)報, 9(2): 174-184.

    李亞男, 沈萍萍, 黃良民, 等, 2012. 棕囊藻的分類及系統(tǒng)進化研究進展[J]. 生態(tài)學(xué)雜志, 31(3): 745-754.

    李偉才, 孫軍, 宋書群, 等, 2006. 煙臺港和鄰近錨地及其入境船舶壓艙水中的浮游植物[J]. 海洋湖沼通報(4): 70-77.

    何家菀, 施之新, 張銀華, 等, 1999. 一種棕囊藻的形態(tài)特征與毒素分析[J]. 海洋與湖沼, 30(2): 172-179.

    沈萍萍, 齊雨藻, 2021. 棕囊藻屬()的種類多樣性及地理分布特征研究進展[J]. 海洋與湖沼, 52(1): 1-15.

    沈萍萍, 齊雨藻, 歐林堅, 2018. 中國沿海球形棕囊藻()的分類、分布及其藻華[J]. 海洋科學(xué), 42(10): 146-162.

    陳菊芳, 徐寧, 江天久, 等, 1999. 中國赤潮新記錄種——球形棕囊藻()[J]. 暨南大學(xué)學(xué)報(自然科學(xué)版), 20(3): 124-129.

    賀成, 宋書群, 李才文, 2019. 廣西北部灣海域球形棕囊藻囊體時空分布及其影響因素[J]. 海洋與湖沼, 50(3): 630-643.

    徐寧, 齊雨藻, 陳菊芳, 等, 2003. 球形棕囊藻(Scherffel)赤潮成因分析[J]. 環(huán)境科學(xué)學(xué)報, 23(1): 113-118.

    曹西華, 俞志明, 邱麗霞, 2017. 改性黏土法消除球形棕囊藻赤潮的現(xiàn)場實驗與效果評估[J]. 海洋與湖沼, 48(4): 753-759.

    屠建波, 張秋豐, 徐玉山, 等, 2011. 渤海灣天津近岸海域首次棕囊藻赤潮初探[J]. 海洋通報, 30(3): 334-337.

    ANTAJAN E, CHRéTIENNOT-DINET M J, LEBLANC C,, 2004. 19′-hexanoyloxyfucoxanthin may not be the appropriate pigment to trace occurrence and fate of: the case ofin Belgian coastal waters [J]. Journal of Sea Research, 52(3): 165-177.

    BEGON M, TOWNSEND C R, HARPER J L, 2005. Ecology: From Individuals to Ecosystems [M]. 4th ed. Oxford: Wiley-Blackwell.

    BUMA A G J, BANO N, VELDHUIS M J W,, 1991. Comparison of the pigmentation of two strains of the prymnesiophytesp. [J]. Netherlands Journal of Sea Research, 27(2): 173-182.

    CHARLSON R J, LOVELOCK J E, ANDREAE M O,, 1987. Oceanic phytoplankton, atmospheric sulphur, cloud albedo and climate [J]. Nature, 326(6114): 655-661.

    CHEN Y Q, SHAO P, WANG N,, 2003. Molecular identification of bloom-forming speciesglobosa (Prymnesiophyta) and its dispersal based on rDNA ITS sequence analysis [J]. Acta Oceanologica Sinica, 22(2): 243-253.

    DITULLIO G R, GARCIA N, RISEMAN S F,, 2007. Effects of iron concentration on pigment composition ingrown at low irradiance [J]. Biogeochemistry, 83(1/2/3): 71-81.

    HAI D N, LAM N N, DIPPNER J W, 2010. Development ofblooms in the upwelling waters of the South Central coast of Viet Nam [J]. Journal of Marine Systems, 83(3/4): 253-261.

    KALYAANAMOORTHY S, MINH B Q, WONG T K F,, 2017. ModelFinder: fast model selection for accurate phylogenetic estimates [J]. Nature Methods, 14(6): 587-589.

    KONG F Z, YU R C, ZHANG Q C,, 2012. Pigment characterization for the 2011 bloom in Qinhuangdao implicated “brown tide” events in China [J]. Chinese Journal of Oceanology and Limnology, 30(3): 361-370.

    KUMAR S, STECHER G, LI M,, 2018. MEGA X: molecular evolutionary genetics analysis across computing platforms [J]. Molecular Biology and Evolution, 35(6): 1547-1549.

    LANCELOT C, GYPENS N, BILLEN G,, 2007. Testing an integrated river-ocean mathematical tool for linking marine eutrophication to land use: the-dominated Belgian coastal zone (Southern North Sea) over the past 50 years [J]. Journal of Marine Systems, 64(1/2/3/4): 216-228.

    LANCELOT C, ROUSSEAU V, GYPENS N, 2009. Ecologically based indicators fordisturbance in eutrophied Belgian coastal waters (Southern North Sea) based on field observations and ecological modelling [J]. Journal of Sea Research, 61(1/2): 44-49.

    LANGE M, CHEN Y Q, MEDLIN L K, 2002. Molecular genetic delineation ofspecies (Prymnesiophyceae) using coding and non-coding regions of nuclear and plastid genomes [J]. European Journal of Phycology, 37(1): 77-92.

    MADHUPRATAP M, SAWANT S, GAUNS M, 2000. A first report on a bloom of the marine prymnesiophycean,from the Arabian Sea [J]. Oceanologica Acta, 23(1): 83-90.

    MEDLIN L, ZINGONE A, 2007. A taxonomic review of the genus[J]. Biogeochemistry, 83(1/2/3): 3-18.

    NIU Z, LIU C, ZHANG Q C,, 2022. Development of sensitive genotype-specific quantitative polymerase chain reaction methods for detection ofin the South China Sea [J]. Limnology and Oceanography: Methods, 20(3): 131-145.

    QI Y Z, CHEN J F, WANG Z H,, 2004. Some observations on harmful algal bloom (HAB) events along the coast of Guangdong, southern China in 1998 [J]. Hydrobiologia, 512(1/2/3): 209-214.

    RONQUIST F, TESLENKO M, VAN DER MARK P,, 2012. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space [J]. Systematic Biology, 61(3): 539-542.

    ROUSSEAU V, MATHOT S, LANCELOT C, 1990. Calculating carbon biomass ofsp. from microscopic observations [J]. Marine Biology, 107(2): 305-314.

    SAZHIN A F, ARTIGAS L F, NEJSTGAARD J C,, 2007. The colonization of twospecies () by pennate diatoms and other protists: a significant contribution to colony biomass [J]. Biogeochemistry, 83(1/2/3): 137-145.

    SCHLüTER L, M?HLENBERG F, HAVSKUM H,, 2000. The use of phytoplankton pigments for identifying and quantifying phytoplankton groups in coastal areas: testing the influence of light and nutrients on pigment/chlorophyllratios [J]. Marine Ecology Progress Series, 192: 49-63.

    SEOANE S, ZAPATA M, ORIVE E, 2009. Growth rates and pigment patterns of haptophytes isolated from estuarine waters [J]. Journal of Sea Research, 62(4): 286-294.

    SMITH D R, ARRIGO K R, ALDERKAMP A C,, 2014a. Massive difference in synonymous substitution rates among mitochondrial, plastid, and nuclear genes ofalgae [J]. Molecular Phylogenetics and Evolution, 71: 36-40.

    SMITH JR W O, LIU X, TANG K W,, 2014b. Giantism and its role in the harmful algal bloom species[J]. Deep Sea Research Part II: Topical Studies in Oceanography, 101: 95-106.

    SONG H Y, CHEN Y, GIBSON K,, 2021. High genetic diversity of the harmful algal bloom speciesrevealed using the molecular marker[J]. Harmful Algae, 107: 102065.

    SONG H Y, LIU F, LI Z L,, 2020. Development of a high-resolution molecular marker for trackinggenetic diversity through comparative analysis of chloroplast genomes [J]. Harmful Algae, 99: 101911.

    TANG K W, 2003. Grazing and colony size development in(Prymnesiophyceae): the role of a chemical signal [J]. Journal of Plankton Research, 25(7): 831-842.

    TURESSON G, 1922. The genotypical response of the plant species to the habitat [J]. Hereditas, 3(3): 211-350.

    VAN LEEUWE M A, STEFELS J, 1998. Effects of iron and light stress on the biochemical composition of Antarcticsp. (Prymnesiophyceae). II. Pigment composition [J]. Journal of Phycology, 34(3): 496-503.

    VAN LEEUWE M A, VISSER R J W, STEFELS J, 2014. The pigment composition of(Haptophyceae) under various conditions of light, temperature, salinity, and iron [J]. Journal of Phycology, 50(6): 1070-1080.

    VAN RIJSSEL M, HAMM C, GIESKES W, 1997.(Prymnesiophyceae) colonies: hollow structures built with small amounts of polysaccharides [J]. European Journal of Phycology, 32(2): 185-192.

    VAULOT D, BIRRIEN J L, MARIE D,, 1994. Morphology, ploidy, pigment composition, and genome size of cultured strains of(Prymnesiophyceae) [J]. Journal of Phycology, 30(6): 1022-1035.

    WANG J X, KONG F Z, GENG H X,, 2022. Pigment characterization of the giant-colony-forming haptophytein the Beibu Gulf reveals blooms of different origins [J]. Applied and Environmental Microbiology, 88(4): e0165421.

    WANG X D, SONG H Y, WANG Y,, 2021. Research on the biology and ecology of the harmful algal bloom speciesin China: progresses in the last 20 years [J]. Harmful Algae, 107: 102057.

    WANG X D, TANG K W, WANG Y,, 2010. Temperature effects on growth, colony development and carbon partitioning in threespecies [J]. Aquatic Biology, 9(3): 239-249.

    WANG X D, WANG Y, SMITH W O, 2011. The role of nitrogen on the growth and colony development of(Prymnesiophyceae) [J]. European Journal of Phycology, 46(3): 305-314.

    WASSMANN P, RATKOVA T, REIGSTAD M, 2005. The contribution of single and colonial cells ofto spring and summer blooms in the north-eastern North Atlantic [J]. Harmful Algae, 4(5): 823-840.

    WHIPPLE S J, PATTEN B C, VERITY P G, 2005. Life cycle of the marine alga: a conceptual model to summarize literature and guide research [J]. Journal of Marine Systems, 57(1/2): 83-110.

    WINNEPENNINCKX B, BACKELJAU T, DEWACHTER R, 1993. Extraction of high molecular weight DNA from molluscs [J]. Trends in Genetics, 9(12): 407.

    YU Z M, SONG X X, CAO X H,, 2017. Mitigation of harmful algal blooms using modified clays: theory, mechanisms, and applications [J]. Harmful Algae, 69: 48-64.

    ZAPATA M, JEFFREY S W, WRIGHT S W,, 2004. Photosynthetic pigments in 37 species (65 strains) of Haptophyta: implications for oceanography and chemotaxonomy [J]. Marine Ecology Progress Series, 270: 83-102.

    ZAPATA M, RODRíGUEZ F, GARRIDO J L, 2000. Separation of chlorophylls and carotenoids from marine phytoplankton: a new HPLC method using a reversed phase C8column and pyridine-containing mobile phases [J]. Marine Ecology Progress Series, 195: 29-45.

    ZHANG Q C, NIU Z, WANG J X,, 2021. Development of high-resolution chloroplast markers for intraspecific phylogeographic studies of[J]. Journal of Oceanology and Limnology, 39(2): 508-524.

    ZHANG Q C, LIU C, WANG J X,, 2022. Intense blooms ofin the South China Sea are caused by a unique “giant-colony” ecotype [J]. Harmful Algae, 114: 102227.

    FIRST REPORT OFBLOOM CAUSED BY “GIANT-COLONY” ECOTYPE IN THE SOUTH YELLOW SEA

    ZHANG Qing-Chun1, 2, 4, XIANG Ling1, 3, WANG Jin-Xiu1, 3, LIU Yang1, 5, SONG Min-Jie1, 3, WANG Yun-Feng1, 2, 4, KONG Fan-Zhou1, 2, 4, GENG Hui-Xia1, 2, 4, YAN Tian1, 2, 4, YU Ren-Cheng1, 2, 3, 4

    (1.CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. Laboratory for Marine Ecology and Environmental Science, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266071, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China; 4.Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China;5. Qingdao University of Science & technology, Qingdao 266045, China)

    , an important causative species of harmful algal blooms worldwide, exhibits varying morphological and physiological features and high genetic diversity. The bloom-formingin the South China Sea (SCS) is proposed to be a unique “giant-colony” ecotype, which features giant colonies and possesses 19′-butanoyloxyfucoxanthin (But-fuco) as its diagnostic pigment, and is genetically distinct from other strains of...blooms occur frequently in the SCS, as well as several times in the Bohai Sea and the North Yellow Sea, but never in the South Yellow Sea (SYS). A.bloom broke out in the coastal waters of Qingdao from late November to early December, 2021. During the bloom, it was found that the surface temperature ranged 12~14 °C, andexhibited high DON (dissolved organic nitrogen) content and low DIN (dissolved inorganic nitrogen) content in the Qingdao coastal waters. The colony abundance was more than 20 colonies/L with the maximum diameter of 18 mm. The sampled colonies from the Qingdao coastal waters contained a unique diagnostic pigment But-fuco, and belonged to the same genetic clade of “giant-colony” ecotype of.distributed in the SCS. This study firstly reported that thebloom in the SYS caused by the “giant-colony” ecotype. The “giant-colony”.bloom may threaten local plankton ecosystem, aquaculture, and coastal nuclear power plant. Therefore, it is our key task to detect and monitor.blooms, and study the mechanism of its occurrence in the northern waters of China.

    ; “giant-colony” ecotype; harmful algal bloom (HAB); South Yellow Sea

    Q14

    10.11693/hyhz20220100017

    *青島海洋科學(xué)與技術(shù)試點國家實驗室山東省專項經(jīng)費, 2021QNLM040001-2; 國家自然科學(xué)基金委面上項目42076140號, U20A004號, 41776127號; 國家科技基礎(chǔ)資源調(diào)查項目, 2018FY100200號。張清春, 副研究員, E-mail: qczhang@qdio.ac.cn

    于仁成, 博士生導(dǎo)師, 研究員, E-mail: rcyu@qdio.ac.cn

    2022-01-21,

    2022-02-09

    猜你喜歡
    囊體生態(tài)型黃海
    多元策略的生態(tài)型游客中心
    氮磷營養(yǎng)變化對欽州灣球形棕囊藻藻體生長的影響
    兩種充氣碟形囊體結(jié)構(gòu)力學(xué)性能比較研究*
    東方濕地 黃海明珠
    華人時刊(2020年21期)2021-01-14 01:32:28
    Agg Hab生態(tài)型聚合棲地
    囊體材料高溫撕裂強度測試方法研究
    黃海簡介
    番茄有機生態(tài)型無土栽培
    三角恒等變換去哪兒了
    微生物內(nèi)核 生態(tài)型農(nóng)資
    99久久精品国产亚洲精品| 好男人电影高清在线观看| av国产免费在线观看| 变态另类丝袜制服| 久久九九热精品免费| 又爽又黄无遮挡网站| 男女那种视频在线观看| 国产亚洲av高清不卡| 最好的美女福利视频网| 黄色视频不卡| 国产熟女午夜一区二区三区| 久久久久久久久免费视频了| 亚洲中文字幕日韩| 国产一区在线观看成人免费| 成在线人永久免费视频| 可以在线观看的亚洲视频| 日韩欧美在线二视频| 国产伦在线观看视频一区| 99久久久亚洲精品蜜臀av| 日韩 欧美 亚洲 中文字幕| 九九热线精品视视频播放| 嫁个100分男人电影在线观看| 久久久久国产一级毛片高清牌| 真人做人爱边吃奶动态| 国产爱豆传媒在线观看 | 黄色片一级片一级黄色片| 天堂影院成人在线观看| 国产精品98久久久久久宅男小说| 亚洲国产看品久久| 国产高清有码在线观看视频 | 中文亚洲av片在线观看爽| 亚洲成人中文字幕在线播放| 亚洲七黄色美女视频| 午夜激情福利司机影院| 亚洲av片天天在线观看| 亚洲av成人不卡在线观看播放网| 长腿黑丝高跟| 久久久国产成人免费| 午夜久久久久精精品| 久久久久国产精品人妻aⅴ院| 亚洲精品在线观看二区| 男女做爰动态图高潮gif福利片| 午夜亚洲福利在线播放| 精品欧美一区二区三区在线| 国产精品亚洲美女久久久| 国产爱豆传媒在线观看 | 好看av亚洲va欧美ⅴa在| 日韩中文字幕欧美一区二区| 亚洲,欧美精品.| 99久久无色码亚洲精品果冻| 成人18禁在线播放| 一二三四在线观看免费中文在| 亚洲成人中文字幕在线播放| 免费人成视频x8x8入口观看| 亚洲国产精品久久男人天堂| 亚洲精品一卡2卡三卡4卡5卡| av国产免费在线观看| 久久人妻福利社区极品人妻图片| 国产69精品久久久久777片 | av中文乱码字幕在线| 麻豆国产97在线/欧美 | 一区二区三区激情视频| 亚洲av美国av| 久久久国产成人免费| 国产成人av激情在线播放| 麻豆成人av在线观看| 草草在线视频免费看| 美女扒开内裤让男人捅视频| 国内少妇人妻偷人精品xxx网站 | 夜夜夜夜夜久久久久| 国产v大片淫在线免费观看| www国产在线视频色| 好男人在线观看高清免费视频| 国产一区二区在线av高清观看| 欧美三级亚洲精品| 成人精品一区二区免费| www.熟女人妻精品国产| 国产区一区二久久| 亚洲五月婷婷丁香| а√天堂www在线а√下载| 变态另类丝袜制服| 久久亚洲精品不卡| 91字幕亚洲| 国产精品免费视频内射| 久久久久久久久中文| 不卡av一区二区三区| 日韩成人在线观看一区二区三区| 国产成人精品无人区| 亚洲欧美精品综合久久99| 国产免费男女视频| 中文字幕人成人乱码亚洲影| 色精品久久人妻99蜜桃| 91老司机精品| 天天躁夜夜躁狠狠躁躁| 最近视频中文字幕2019在线8| 无人区码免费观看不卡| 国产午夜精品论理片| 色在线成人网| 欧美黑人巨大hd| svipshipincom国产片| 国产乱人伦免费视频| 国内揄拍国产精品人妻在线| 国模一区二区三区四区视频 | 在线a可以看的网站| 欧美午夜高清在线| 欧美乱妇无乱码| 怎么达到女性高潮| 嫩草影视91久久| 露出奶头的视频| 亚洲人成伊人成综合网2020| 99国产精品一区二区三区| 搡老岳熟女国产| а√天堂www在线а√下载| 午夜福利成人在线免费观看| 久久香蕉精品热| 岛国视频午夜一区免费看| 亚洲五月婷婷丁香| 成人精品一区二区免费| 精华霜和精华液先用哪个| 久久精品国产亚洲av高清一级| 精品久久蜜臀av无| 久热爱精品视频在线9| 一进一出好大好爽视频| 制服诱惑二区| 午夜日韩欧美国产| 啦啦啦免费观看视频1| 色播亚洲综合网| 丰满人妻熟妇乱又伦精品不卡| 亚洲专区国产一区二区| 亚洲国产欧洲综合997久久,| 精品一区二区三区四区五区乱码| 欧美日韩瑟瑟在线播放| 免费av毛片视频| 午夜福利视频1000在线观看| 成人一区二区视频在线观看| 国产三级在线视频| 久久亚洲真实| 中文字幕久久专区| 最好的美女福利视频网| 亚洲精华国产精华精| 男女那种视频在线观看| 欧美又色又爽又黄视频| 人成视频在线观看免费观看| 啦啦啦观看免费观看视频高清| 老熟妇仑乱视频hdxx| 嫩草影院精品99| 久久精品国产99精品国产亚洲性色| 两性午夜刺激爽爽歪歪视频在线观看 | 99久久精品热视频| 一级作爱视频免费观看| 一本大道久久a久久精品| 亚洲精品久久成人aⅴ小说| 中亚洲国语对白在线视频| 国产精品一区二区免费欧美| 人人妻,人人澡人人爽秒播| 午夜亚洲福利在线播放| 国产精品一区二区免费欧美| 别揉我奶头~嗯~啊~动态视频| 色综合婷婷激情| av中文乱码字幕在线| 两个人视频免费观看高清| 国产久久久一区二区三区| 一级黄色大片毛片| videosex国产| 丰满人妻熟妇乱又伦精品不卡| 色哟哟哟哟哟哟| 九九热线精品视视频播放| 久99久视频精品免费| 人人妻人人看人人澡| av国产免费在线观看| 日韩高清综合在线| 国产成人精品久久二区二区免费| 亚洲美女视频黄频| 久久九九热精品免费| 精品高清国产在线一区| 欧美在线黄色| 麻豆国产av国片精品| 9191精品国产免费久久| 在线观看一区二区三区| 毛片女人毛片| 不卡av一区二区三区| 久久天堂一区二区三区四区| 免费在线观看亚洲国产| 日本成人三级电影网站| 最近视频中文字幕2019在线8| 欧美极品一区二区三区四区| 国产高清视频在线播放一区| 免费人成视频x8x8入口观看| 99国产精品99久久久久| 丰满的人妻完整版| 97超级碰碰碰精品色视频在线观看| 亚洲午夜精品一区,二区,三区| 黄片小视频在线播放| 一进一出抽搐动态| 国产精品一区二区精品视频观看| 久久精品综合一区二区三区| 国产久久久一区二区三区| 国产爱豆传媒在线观看 | 老汉色∧v一级毛片| 国产精品自产拍在线观看55亚洲| 黑人欧美特级aaaaaa片| 亚洲国产看品久久| aaaaa片日本免费| 在线a可以看的网站| a级毛片a级免费在线| 精品久久久久久成人av| 国产三级中文精品| 久久精品成人免费网站| 欧美 亚洲 国产 日韩一| 88av欧美| √禁漫天堂资源中文www| 午夜福利高清视频| 午夜福利在线观看吧| 一级毛片高清免费大全| 久久99热这里只有精品18| 国产麻豆成人av免费视频| 手机成人av网站| 国产一区二区在线av高清观看| 亚洲精品国产一区二区精华液| 黄色 视频免费看| 99国产精品一区二区蜜桃av| 国产亚洲精品综合一区在线观看 | 成人欧美大片| 久久99热这里只有精品18| 免费在线观看亚洲国产| 一边摸一边做爽爽视频免费| 91av网站免费观看| 怎么达到女性高潮| 日本 av在线| 一进一出抽搐动态| 999久久久精品免费观看国产| 我要搜黄色片| 中文字幕久久专区| 精品国产美女av久久久久小说| av片东京热男人的天堂| 一进一出抽搐gif免费好疼| 亚洲黑人精品在线| 亚洲中文日韩欧美视频| 一二三四社区在线视频社区8| 特大巨黑吊av在线直播| 亚洲精品一区av在线观看| 男人的好看免费观看在线视频 | 在线观看免费日韩欧美大片| 激情在线观看视频在线高清| 日韩欧美免费精品| 一级毛片高清免费大全| 十八禁网站免费在线| 亚洲av电影不卡..在线观看| 男女午夜视频在线观看| 好男人在线观看高清免费视频| 手机成人av网站| 国产精品爽爽va在线观看网站| 黄色丝袜av网址大全| 亚洲真实伦在线观看| 男人的好看免费观看在线视频 | 高清毛片免费观看视频网站| 欧美不卡视频在线免费观看 | 欧美日韩精品网址| 俄罗斯特黄特色一大片| 可以在线观看毛片的网站| 婷婷六月久久综合丁香| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美精品v在线| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美精品综合一区二区三区| 亚洲一码二码三码区别大吗| 又黄又粗又硬又大视频| 亚洲性夜色夜夜综合| 精品国产乱子伦一区二区三区| 日本 av在线| 亚洲免费av在线视频| 老司机靠b影院| 99久久无色码亚洲精品果冻| 亚洲一区中文字幕在线| 国产伦人伦偷精品视频| 久久 成人 亚洲| 精品人妻1区二区| 亚洲欧美一区二区三区黑人| 国产亚洲精品综合一区在线观看 | 亚洲av成人一区二区三| 他把我摸到了高潮在线观看| 女警被强在线播放| 精品国产亚洲在线| 国产伦人伦偷精品视频| 欧美久久黑人一区二区| 国产真人三级小视频在线观看| 99热6这里只有精品| 在线观看www视频免费| 69av精品久久久久久| 亚洲国产欧美网| 99精品欧美一区二区三区四区| www.精华液| 久久国产乱子伦精品免费另类| 一级作爱视频免费观看| 老司机靠b影院| 这个男人来自地球电影免费观看| 老司机在亚洲福利影院| 两个人视频免费观看高清| 日韩欧美在线乱码| 国产精品爽爽va在线观看网站| 亚洲欧美日韩高清在线视频| 亚洲18禁久久av| 国产片内射在线| 国产成人系列免费观看| 老司机午夜十八禁免费视频| 欧美+亚洲+日韩+国产| 好男人在线观看高清免费视频| 久久国产精品影院| 夜夜爽天天搞| 午夜激情av网站| 免费观看人在逋| 黄色毛片三级朝国网站| 亚洲熟妇熟女久久| 国产探花在线观看一区二区| 哪里可以看免费的av片| 久久久精品大字幕| 深夜精品福利| 日本一本二区三区精品| 99热只有精品国产| 伦理电影免费视频| 99久久精品热视频| 级片在线观看| 激情在线观看视频在线高清| 国产伦人伦偷精品视频| 韩国av一区二区三区四区| 亚洲av成人av| 99精品在免费线老司机午夜| 亚洲成人国产一区在线观看| 日韩欧美在线二视频| 啦啦啦免费观看视频1| 国产69精品久久久久777片 | 亚洲欧美一区二区三区黑人| 国产亚洲精品一区二区www| 免费av毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 在线观看一区二区三区| 国产成人啪精品午夜网站| 亚洲国产精品999在线| 在线看三级毛片| 久久精品人妻少妇| 黄色片一级片一级黄色片| 免费在线观看影片大全网站| 好男人在线观看高清免费视频| 成人特级黄色片久久久久久久| 热99re8久久精品国产| 精品久久久久久久久久久久久| 亚洲精品在线美女| 天天添夜夜摸| 女生性感内裤真人,穿戴方法视频| 免费人成视频x8x8入口观看| 色av中文字幕| 亚洲 欧美 日韩 在线 免费| 亚洲av日韩精品久久久久久密| 麻豆国产av国片精品| 久热爱精品视频在线9| 国内毛片毛片毛片毛片毛片| 国产人伦9x9x在线观看| 国产爱豆传媒在线观看 | 国产精品一区二区免费欧美| 亚洲欧美精品综合一区二区三区| 国产av在哪里看| 精品久久久久久久末码| 男男h啪啪无遮挡| 在线视频色国产色| 丝袜人妻中文字幕| 91在线观看av| 国产私拍福利视频在线观看| 欧美人与性动交α欧美精品济南到| 波多野结衣巨乳人妻| 啦啦啦免费观看视频1| 成年人黄色毛片网站| 日韩欧美在线乱码| 国产主播在线观看一区二区| 禁无遮挡网站| 亚洲午夜理论影院| 久久天堂一区二区三区四区| av片东京热男人的天堂| 精品福利观看| 国产一区二区在线av高清观看| 老司机福利观看| 又粗又爽又猛毛片免费看| 亚洲黑人精品在线| 久久国产精品人妻蜜桃| 午夜影院日韩av| 高清毛片免费观看视频网站| 9191精品国产免费久久| 亚洲,欧美精品.| 成年女人毛片免费观看观看9| 妹子高潮喷水视频| 99国产综合亚洲精品| xxx96com| 日韩成人在线观看一区二区三区| 久久性视频一级片| 9191精品国产免费久久| 深夜精品福利| 欧美3d第一页| www日本在线高清视频| 亚洲国产欧美人成| 欧美成人免费av一区二区三区| 国产精品影院久久| 桃红色精品国产亚洲av| 久久久久久久久免费视频了| 男人舔女人下体高潮全视频| 国产视频内射| 在线a可以看的网站| 久久中文字幕一级| 日日摸夜夜添夜夜添小说| 亚洲激情在线av| 色av中文字幕| 黑人巨大精品欧美一区二区mp4| 亚洲精品粉嫩美女一区| 在线看三级毛片| 悠悠久久av| 男女之事视频高清在线观看| 黄色丝袜av网址大全| 免费看a级黄色片| 午夜福利视频1000在线观看| 国产1区2区3区精品| av国产免费在线观看| 国产真实乱freesex| 欧美激情久久久久久爽电影| 国产蜜桃级精品一区二区三区| 麻豆成人午夜福利视频| 看片在线看免费视频| 久久久久久久久久黄片| 午夜免费激情av| 日本一本二区三区精品| 两个人的视频大全免费| 天天添夜夜摸| 老司机午夜福利在线观看视频| 天天一区二区日本电影三级| 色噜噜av男人的天堂激情| 禁无遮挡网站| 午夜日韩欧美国产| 久久精品综合一区二区三区| 午夜福利成人在线免费观看| 老汉色∧v一级毛片| 色播亚洲综合网| 久久久久免费精品人妻一区二区| 亚洲成人中文字幕在线播放| 亚洲精华国产精华精| 好看av亚洲va欧美ⅴa在| 欧美黄色片欧美黄色片| 亚洲avbb在线观看| 小说图片视频综合网站| 久久精品国产亚洲av香蕉五月| 免费看日本二区| 日韩成人在线观看一区二区三区| 母亲3免费完整高清在线观看| 看黄色毛片网站| 久久精品国产清高在天天线| 狂野欧美激情性xxxx| 午夜视频精品福利| www.www免费av| 亚洲午夜精品一区,二区,三区| 免费看美女性在线毛片视频| 一级毛片女人18水好多| 欧美午夜高清在线| 熟女少妇亚洲综合色aaa.| 丁香欧美五月| 国产97色在线日韩免费| 亚洲 国产 在线| 麻豆成人午夜福利视频| 久久精品91蜜桃| 这个男人来自地球电影免费观看| 老熟妇乱子伦视频在线观看| 叶爱在线成人免费视频播放| 波多野结衣巨乳人妻| 变态另类丝袜制服| 91老司机精品| 国产精品,欧美在线| 黄色丝袜av网址大全| 久久精品aⅴ一区二区三区四区| 国产精品久久电影中文字幕| 波多野结衣高清作品| 看黄色毛片网站| 欧美zozozo另类| 在线观看www视频免费| 精品无人区乱码1区二区| 一级a爱片免费观看的视频| 亚洲天堂国产精品一区在线| 露出奶头的视频| 老熟妇乱子伦视频在线观看| 欧美人与性动交α欧美精品济南到| 最好的美女福利视频网| 伦理电影免费视频| 91成年电影在线观看| 日韩国内少妇激情av| www.www免费av| 一区二区三区激情视频| 两性夫妻黄色片| 一区二区三区高清视频在线| 婷婷丁香在线五月| 一二三四社区在线视频社区8| 久久久久免费精品人妻一区二区| 一边摸一边抽搐一进一小说| 好男人电影高清在线观看| 成人亚洲精品av一区二区| 国产高清视频在线观看网站| 好看av亚洲va欧美ⅴa在| 欧美成人一区二区免费高清观看 | 亚洲精品粉嫩美女一区| 欧美日韩瑟瑟在线播放| 久9热在线精品视频| 久久人人精品亚洲av| 国产高清videossex| 视频区欧美日本亚洲| 成年版毛片免费区| 99在线视频只有这里精品首页| 亚洲欧美日韩无卡精品| 精品第一国产精品| 91av网站免费观看| 动漫黄色视频在线观看| 最新在线观看一区二区三区| 国产91精品成人一区二区三区| 久久久久久久久中文| av超薄肉色丝袜交足视频| 精品电影一区二区在线| 国产激情久久老熟女| 免费高清视频大片| 99国产精品一区二区三区| 成年人黄色毛片网站| 日韩成人在线观看一区二区三区| 亚洲在线自拍视频| 亚洲欧美一区二区三区黑人| 三级国产精品欧美在线观看 | 欧美不卡视频在线免费观看 | 国产一区二区在线观看日韩 | 操出白浆在线播放| 亚洲国产看品久久| 久久 成人 亚洲| 男女视频在线观看网站免费 | 美女扒开内裤让男人捅视频| 亚洲精品国产精品久久久不卡| 又黄又爽又免费观看的视频| 国产欧美日韩一区二区精品| 国产一区二区在线观看日韩 | 国内精品久久久久久久电影| 久久精品国产清高在天天线| 好看av亚洲va欧美ⅴa在| 亚洲五月婷婷丁香| 亚洲国产精品999在线| 50天的宝宝边吃奶边哭怎么回事| а√天堂www在线а√下载| 人妻夜夜爽99麻豆av| 色av中文字幕| 国产亚洲欧美在线一区二区| 亚洲欧美精品综合一区二区三区| 精品一区二区三区四区五区乱码| 国产91精品成人一区二区三区| 久久久久免费精品人妻一区二区| 成人手机av| 美女大奶头视频| 脱女人内裤的视频| 老鸭窝网址在线观看| 成人国产一区最新在线观看| svipshipincom国产片| 亚洲人成77777在线视频| 亚洲人成网站在线播放欧美日韩| 欧美丝袜亚洲另类 | 精品久久久久久久末码| 亚洲精品在线观看二区| 最近最新中文字幕大全电影3| 日韩欧美免费精品| 亚洲成人久久性| 日韩欧美在线二视频| 精品久久久久久久毛片微露脸| 久久国产精品影院| 成人国语在线视频| 99在线人妻在线中文字幕| 老司机在亚洲福利影院| 久久久久性生活片| 亚洲精品色激情综合| 日韩精品中文字幕看吧| 国产精品美女特级片免费视频播放器 | 久久精品国产清高在天天线| 一本大道久久a久久精品| 99在线视频只有这里精品首页| 一区福利在线观看| 国产成人精品久久二区二区免费| 黄片大片在线免费观看| 无遮挡黄片免费观看| 免费在线观看日本一区| 亚洲五月天丁香| 99国产极品粉嫩在线观看| 观看免费一级毛片| а√天堂www在线а√下载| 大型黄色视频在线免费观看| 黄色丝袜av网址大全| 一个人免费在线观看电影 | 国产真人三级小视频在线观看| 日韩欧美国产在线观看| 琪琪午夜伦伦电影理论片6080| 欧美 亚洲 国产 日韩一| 1024视频免费在线观看| 久久久久久人人人人人| 欧美3d第一页| 国产精品一及| 亚洲av中文字字幕乱码综合| 亚洲av第一区精品v没综合| 中文亚洲av片在线观看爽| 美女 人体艺术 gogo| 99re在线观看精品视频| 久久人妻av系列| 淫秽高清视频在线观看| 免费无遮挡裸体视频| 欧美成狂野欧美在线观看| 午夜免费成人在线视频| 亚洲精品国产一区二区精华液| 国产精品一区二区三区四区久久| 亚洲av片天天在线观看| 老鸭窝网址在线观看| 国产高清视频在线播放一区| 中出人妻视频一区二区| 亚洲,欧美精品.|