• <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)資
    久久ye,这里只有精品| 乱码一卡2卡4卡精品| 久久久国产一区二区| 免费少妇av软件| 水蜜桃什么品种好| 狂野欧美激情性xxxx在线观看| 丝袜脚勾引网站| 国产av精品麻豆| 精品亚洲成国产av| 少妇猛男粗大的猛烈进出视频| 日韩成人av中文字幕在线观看| 亚洲精品亚洲一区二区| 久久99一区二区三区| 欧美三级亚洲精品| 色94色欧美一区二区| 免费大片黄手机在线观看| 国产69精品久久久久777片| 国产av国产精品国产| 黄色配什么色好看| av线在线观看网站| 午夜激情av网站| 久久人人爽av亚洲精品天堂| 日韩av免费高清视频| 免费看光身美女| 亚洲av不卡在线观看| 97精品久久久久久久久久精品| 中文字幕精品免费在线观看视频 | 亚洲av不卡在线观看| av不卡在线播放| 久久精品久久精品一区二区三区| 国产精品免费大片| 交换朋友夫妻互换小说| 精品人妻熟女毛片av久久网站| 99精国产麻豆久久婷婷| 亚洲精品日韩av片在线观看| 有码 亚洲区| 欧美日韩视频高清一区二区三区二| 我要看黄色一级片免费的| 妹子高潮喷水视频| 亚洲第一av免费看| 99精国产麻豆久久婷婷| 国产欧美日韩综合在线一区二区| 搡老乐熟女国产| 99九九在线精品视频| 日韩一区二区三区影片| 国产无遮挡羞羞视频在线观看| 高清在线视频一区二区三区| 国产av一区二区精品久久| 国产探花极品一区二区| 看十八女毛片水多多多| 免费高清在线观看日韩| 色94色欧美一区二区| 一级毛片黄色毛片免费观看视频| 亚洲丝袜综合中文字幕| 久久99热6这里只有精品| 蜜臀久久99精品久久宅男| 午夜福利网站1000一区二区三区| 男人操女人黄网站| 久久99蜜桃精品久久| 欧美日韩精品成人综合77777| 蜜臀久久99精品久久宅男| av免费观看日本| 欧美精品一区二区免费开放| 七月丁香在线播放| 亚洲av电影在线观看一区二区三区| av国产久精品久网站免费入址| 水蜜桃什么品种好| 97在线人人人人妻| 性高湖久久久久久久久免费观看| 七月丁香在线播放| 国产精品99久久99久久久不卡 | 欧美日韩视频精品一区| 丝袜美足系列| 国产精品麻豆人妻色哟哟久久| 天堂中文最新版在线下载| 一本—道久久a久久精品蜜桃钙片| videossex国产| av在线app专区| 熟女电影av网| 男女边吃奶边做爰视频| 妹子高潮喷水视频| 夫妻性生交免费视频一级片| 啦啦啦在线观看免费高清www| 久久久久网色| 亚洲性久久影院| 91国产中文字幕| 久久精品久久精品一区二区三区| 91成人精品电影| 日本-黄色视频高清免费观看| 国产欧美另类精品又又久久亚洲欧美| 欧美xxxx性猛交bbbb| 另类亚洲欧美激情| 国产一区二区在线观看日韩| 一级毛片我不卡| 飞空精品影院首页| 亚洲少妇的诱惑av| 一级二级三级毛片免费看| 亚洲精品色激情综合| 五月玫瑰六月丁香| 男女边吃奶边做爰视频| 久久青草综合色| av国产久精品久网站免费入址| 在现免费观看毛片| 晚上一个人看的免费电影| 一级毛片aaaaaa免费看小| 亚洲人成77777在线视频| 最新的欧美精品一区二区| 99久久中文字幕三级久久日本| 好男人视频免费观看在线| 乱码一卡2卡4卡精品| 高清午夜精品一区二区三区| 日韩大片免费观看网站| 亚洲欧美一区二区三区国产| 精品人妻熟女毛片av久久网站| 免费播放大片免费观看视频在线观看| 大香蕉97超碰在线| 涩涩av久久男人的天堂| 国产淫语在线视频| 黄色欧美视频在线观看| 精品国产国语对白av| 春色校园在线视频观看| 丝袜脚勾引网站| 亚洲成人一二三区av| 狠狠精品人妻久久久久久综合| 在线播放无遮挡| 午夜福利在线观看免费完整高清在| 一区在线观看完整版| 激情五月婷婷亚洲| 人妻人人澡人人爽人人| 啦啦啦啦在线视频资源| 精品酒店卫生间| 国产成人免费观看mmmm| 99热全是精品| 亚洲国产日韩一区二区| 亚洲国产欧美在线一区| 超碰97精品在线观看| 人人妻人人澡人人爽人人夜夜| 最近最新中文字幕免费大全7| 下体分泌物呈黄色| 蜜桃在线观看..| 街头女战士在线观看网站| 中国国产av一级| 国产精品久久久久久精品电影小说| 99热这里只有是精品在线观看| 日韩欧美精品免费久久| 免费看av在线观看网站| 婷婷成人精品国产| 日韩精品免费视频一区二区三区 | 在线天堂最新版资源| h视频一区二区三区| 欧美激情国产日韩精品一区| 欧美精品人与动牲交sv欧美| 各种免费的搞黄视频| 人人妻人人添人人爽欧美一区卜| 国产 精品1| 日韩欧美一区视频在线观看| 午夜福利网站1000一区二区三区| 一本色道久久久久久精品综合| 亚洲欧美日韩卡通动漫| 精品少妇内射三级| 久久久久久久大尺度免费视频| 80岁老熟妇乱子伦牲交| 亚洲欧美一区二区三区黑人 | 亚洲精品色激情综合| 日本爱情动作片www.在线观看| 午夜激情福利司机影院| 免费人成在线观看视频色| 男女免费视频国产| 99热网站在线观看| 国产男人的电影天堂91| 美女内射精品一级片tv| 少妇熟女欧美另类| 日本91视频免费播放| 亚洲中文av在线| 日韩中字成人| 啦啦啦视频在线资源免费观看| 人妻系列 视频| 91精品三级在线观看| 制服诱惑二区| 色视频在线一区二区三区| 亚洲成人一二三区av| 亚洲,欧美,日韩| 亚洲av欧美aⅴ国产| 色婷婷av一区二区三区视频| 黄色毛片三级朝国网站| 亚洲精品一区蜜桃| 一级毛片黄色毛片免费观看视频| 亚洲av国产av综合av卡| 多毛熟女@视频| 日本欧美国产在线视频| 久久久国产欧美日韩av| 99视频精品全部免费 在线| 欧美精品一区二区免费开放| 国产男女内射视频| 我要看黄色一级片免费的| 精品亚洲成a人片在线观看| 久久久久久久久久久免费av| 婷婷色麻豆天堂久久| 精品午夜福利在线看| 丰满迷人的少妇在线观看| 男人爽女人下面视频在线观看| 中文字幕免费在线视频6| 在线精品无人区一区二区三| 男女边摸边吃奶| 免费av中文字幕在线| av不卡在线播放| 国产色婷婷99| 国产成人一区二区在线| 熟妇人妻不卡中文字幕| 欧美97在线视频| 伦精品一区二区三区| 美女大奶头黄色视频| 国产极品天堂在线| 777米奇影视久久| 亚洲精品国产av成人精品| 中文字幕制服av| 免费大片黄手机在线观看| 只有这里有精品99| 久久ye,这里只有精品| 秋霞在线观看毛片| 日韩在线高清观看一区二区三区| 亚洲精品日韩在线中文字幕| av一本久久久久| 飞空精品影院首页| freevideosex欧美| 2021少妇久久久久久久久久久| 亚洲精华国产精华液的使用体验| 亚洲综合色网址| 99热这里只有是精品在线观看| 国产亚洲av片在线观看秒播厂| 在现免费观看毛片| 99热这里只有是精品在线观看| 日日啪夜夜爽| 亚洲av成人精品一区久久| 久久亚洲国产成人精品v| 久久久精品区二区三区| 欧美精品高潮呻吟av久久| 亚洲成人av在线免费| 免费人妻精品一区二区三区视频| 一边亲一边摸免费视频| 亚洲国产精品成人久久小说| 国产在线一区二区三区精| 久久毛片免费看一区二区三区| 日韩av在线免费看完整版不卡| 国产一区二区在线观看日韩| 国产av精品麻豆| 国产精品偷伦视频观看了| 精品午夜福利在线看| 日韩亚洲欧美综合| 中文字幕久久专区| 美女大奶头黄色视频| 18禁在线播放成人免费| 熟妇人妻不卡中文字幕| 亚洲第一区二区三区不卡| 久久久久久久精品精品| 赤兔流量卡办理| 最近最新中文字幕免费大全7| 国产精品久久久久久精品古装| 两个人的视频大全免费| 夜夜骑夜夜射夜夜干| 亚洲久久久国产精品| 午夜久久久在线观看| 99久久人妻综合| 在线观看免费日韩欧美大片 | 香蕉精品网在线| 天天躁夜夜躁狠狠久久av| 午夜免费鲁丝| 国产午夜精品一二区理论片| 黑人高潮一二区| 嘟嘟电影网在线观看| 男的添女的下面高潮视频| 欧美日韩国产mv在线观看视频| 麻豆成人av视频| 一本色道久久久久久精品综合| 五月天丁香电影| 亚洲欧美成人精品一区二区| 高清av免费在线| 飞空精品影院首页| av播播在线观看一区| 国产老妇伦熟女老妇高清| 夫妻性生交免费视频一级片| 91精品伊人久久大香线蕉| 亚洲欧美成人精品一区二区| 国产精品欧美亚洲77777| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩视频精品一区| 国产高清有码在线观看视频| 国产精品不卡视频一区二区| 男女高潮啪啪啪动态图| 一区二区三区免费毛片| 亚洲国产精品999| 国产国拍精品亚洲av在线观看| 高清毛片免费看| 日本与韩国留学比较| 国产精品嫩草影院av在线观看| 久久久国产精品麻豆| 国产亚洲最大av| 亚洲欧美色中文字幕在线| 久久久久精品久久久久真实原创| 妹子高潮喷水视频| 国产精品偷伦视频观看了| 黄色视频在线播放观看不卡| av有码第一页| 91在线精品国自产拍蜜月| 欧美变态另类bdsm刘玥| 国产男女内射视频| 大又大粗又爽又黄少妇毛片口| 菩萨蛮人人尽说江南好唐韦庄| 成人二区视频| 久久久午夜欧美精品| 老司机亚洲免费影院| 99国产综合亚洲精品| 久久国产亚洲av麻豆专区| 尾随美女入室| 99九九线精品视频在线观看视频| 欧美激情 高清一区二区三区| 有码 亚洲区| 亚洲第一av免费看| 啦啦啦在线观看免费高清www| 国产亚洲精品第一综合不卡 | 久久国产亚洲av麻豆专区| 欧美+日韩+精品| 国产亚洲午夜精品一区二区久久| 国产一区有黄有色的免费视频| 日韩不卡一区二区三区视频在线| 九九爱精品视频在线观看| 日本色播在线视频| 特大巨黑吊av在线直播| 欧美日韩精品成人综合77777| 日本黄色片子视频| 黄片无遮挡物在线观看| 成人二区视频| 美女主播在线视频| 嘟嘟电影网在线观看| 中文欧美无线码| 欧美xxxx性猛交bbbb| 男女免费视频国产| av不卡在线播放| 在线天堂最新版资源| 美女福利国产在线| 国产免费现黄频在线看| 王馨瑶露胸无遮挡在线观看| 精品国产一区二区三区久久久樱花| 毛片一级片免费看久久久久| 国产免费现黄频在线看| 午夜91福利影院| .国产精品久久| 亚洲精品av麻豆狂野| 少妇的逼水好多| 免费久久久久久久精品成人欧美视频 | 国产欧美日韩综合在线一区二区| 日日摸夜夜添夜夜爱| 国产探花极品一区二区| 高清不卡的av网站| 国产白丝娇喘喷水9色精品| 黄片播放在线免费| 国产日韩欧美在线精品| 亚洲国产欧美在线一区| 精品一品国产午夜福利视频| 国国产精品蜜臀av免费| 18+在线观看网站| 狂野欧美激情性xxxx在线观看| 午夜免费男女啪啪视频观看| 中文字幕最新亚洲高清| 欧美人与善性xxx| 免费大片18禁| 欧美激情极品国产一区二区三区 | 久热久热在线精品观看| 精品久久久久久久久亚洲| 欧美成人精品欧美一级黄| 欧美日韩视频精品一区| 中文字幕免费在线视频6| 久久精品久久久久久噜噜老黄| 在线亚洲精品国产二区图片欧美 | tube8黄色片| 成人国产av品久久久| 大陆偷拍与自拍| 啦啦啦中文免费视频观看日本| 日本黄色日本黄色录像| 少妇猛男粗大的猛烈进出视频| 女人久久www免费人成看片| 亚洲,一卡二卡三卡| 99国产精品免费福利视频| 国产成人av激情在线播放 | 有码 亚洲区| 久久这里有精品视频免费| √禁漫天堂资源中文www| 80岁老熟妇乱子伦牲交| 黑人欧美特级aaaaaa片| 欧美激情极品国产一区二区三区 | 久久久久久久国产电影| 99热网站在线观看| 丁香六月天网| 精品久久久久久久久av| 成人午夜精彩视频在线观看| 9色porny在线观看| 九草在线视频观看| 一级毛片我不卡| 肉色欧美久久久久久久蜜桃| 久久久久网色| 视频在线观看一区二区三区| 亚洲av综合色区一区| 亚洲色图 男人天堂 中文字幕 | 国产亚洲精品第一综合不卡 | 日本欧美国产在线视频| 国产精品 国内视频| 国产精品久久久久久精品电影小说| 成人毛片a级毛片在线播放| 亚洲人与动物交配视频| 视频在线观看一区二区三区| 亚洲精品456在线播放app| 一区在线观看完整版| 亚洲国产欧美日韩在线播放| 日韩制服骚丝袜av| 夜夜爽夜夜爽视频| 狂野欧美激情性xxxx在线观看| 久久av网站| 观看美女的网站| 看免费成人av毛片| 精品卡一卡二卡四卡免费| 亚洲成人av在线免费| 高清欧美精品videossex| 男女免费视频国产| 麻豆乱淫一区二区| 99热国产这里只有精品6| av黄色大香蕉| 午夜激情久久久久久久| 国产精品人妻久久久久久| 午夜久久久在线观看| 亚洲国产毛片av蜜桃av| 视频区图区小说| 国产老妇伦熟女老妇高清| 成人免费观看视频高清| 久久影院123| 久久狼人影院| av福利片在线| 久久午夜福利片| 中文字幕精品免费在线观看视频 | 美女cb高潮喷水在线观看| 亚洲在久久综合| √禁漫天堂资源中文www| 国产精品麻豆人妻色哟哟久久| 国产熟女午夜一区二区三区 | 精品少妇黑人巨大在线播放| 国产视频首页在线观看| 看非洲黑人一级黄片| 国产成人aa在线观看| 九九在线视频观看精品| 久久青草综合色| 国产精品国产三级国产专区5o| 久久婷婷青草| 日日撸夜夜添| 亚洲少妇的诱惑av| 久久99精品国语久久久| 亚洲精品日本国产第一区| 欧美三级亚洲精品| 亚洲成人手机| 免费看光身美女| 在线观看国产h片| av线在线观看网站| 亚洲av成人精品一二三区| av天堂久久9| 韩国av在线不卡| 免费观看在线日韩| 国产高清有码在线观看视频| 少妇被粗大的猛进出69影院 | 91精品一卡2卡3卡4卡| 大片电影免费在线观看免费| 久久久久久久大尺度免费视频| 国产精品成人在线| 人妻少妇偷人精品九色| 人人妻人人澡人人爽人人夜夜| 制服人妻中文乱码| 免费大片黄手机在线观看| 亚洲综合精品二区| freevideosex欧美| 国产深夜福利视频在线观看| 国产男女内射视频| 亚洲精品美女久久av网站| 人妻人人澡人人爽人人| 亚洲av国产av综合av卡| 中文字幕制服av| 成人国语在线视频| 在线播放无遮挡| 久久韩国三级中文字幕| 边亲边吃奶的免费视频| 飞空精品影院首页| 国产一区亚洲一区在线观看| 久久亚洲国产成人精品v| 国产精品久久久久久av不卡| 国产成人精品无人区| 少妇 在线观看| 日韩中字成人| 国产精品麻豆人妻色哟哟久久| 日韩强制内射视频| 久久久久久久久久久免费av| 欧美日韩av久久| 一区在线观看完整版| 蜜桃国产av成人99| 欧美日韩视频精品一区| 插阴视频在线观看视频| 下体分泌物呈黄色| 在线天堂最新版资源| 好男人视频免费观看在线| 久久久国产精品麻豆| 美女福利国产在线| 波野结衣二区三区在线| 午夜福利视频精品| 人人妻人人澡人人爽人人夜夜| 美女福利国产在线| 国产免费现黄频在线看| 综合色丁香网| 亚洲欧美中文字幕日韩二区| 91久久精品电影网| 少妇人妻精品综合一区二区| 91午夜精品亚洲一区二区三区| 国产精品嫩草影院av在线观看| 精品熟女少妇av免费看| 女性被躁到高潮视频| 亚洲精品自拍成人| 啦啦啦啦在线视频资源| 免费播放大片免费观看视频在线观看| 亚洲少妇的诱惑av| 精品午夜福利在线看| 亚洲第一av免费看| 日韩中文字幕视频在线看片| 日韩电影二区| 国产男女内射视频| 少妇的逼水好多| 色视频在线一区二区三区| 国产欧美亚洲国产| 人妻人人澡人人爽人人| 熟女电影av网| 日本欧美国产在线视频| 日韩av不卡免费在线播放| 最近2019中文字幕mv第一页| 日产精品乱码卡一卡2卡三| 久久精品国产亚洲网站| 91国产中文字幕| 人妻系列 视频| 日本av免费视频播放| 免费看av在线观看网站| 亚洲色图 男人天堂 中文字幕 | 成人漫画全彩无遮挡| 亚洲国产精品999| 黑人欧美特级aaaaaa片| 国产精品蜜桃在线观看| 亚洲国产精品999| 久久99热这里只频精品6学生| 人妻制服诱惑在线中文字幕| 久久 成人 亚洲| 最近中文字幕2019免费版| 久久国内精品自在自线图片| 制服丝袜香蕉在线| 国产精品久久久久久精品电影小说| 久久久久久久久久久久大奶| 99久久精品国产国产毛片| 永久网站在线| 在现免费观看毛片| 欧美精品人与动牲交sv欧美| 免费黄色在线免费观看| 午夜福利在线观看免费完整高清在| 精品人妻偷拍中文字幕| 能在线免费看毛片的网站| 欧美日韩综合久久久久久| 我要看黄色一级片免费的| 欧美精品一区二区大全| 欧美日韩av久久| 大香蕉久久成人网| 欧美日韩亚洲高清精品| 日韩精品有码人妻一区| 老女人水多毛片| 亚洲av福利一区| 国产精品国产av在线观看| 黄色怎么调成土黄色| 伦理电影大哥的女人| 王馨瑶露胸无遮挡在线观看| 欧美日韩成人在线一区二区| 国产日韩欧美在线精品| 久久99热这里只频精品6学生| av网站免费在线观看视频| 乱人伦中国视频| 久久人人爽av亚洲精品天堂| 日韩熟女老妇一区二区性免费视频| 最近手机中文字幕大全| 久久久久久久久久久久大奶| 七月丁香在线播放| 亚洲av不卡在线观看| 最近最新中文字幕免费大全7| av在线app专区| 在线观看三级黄色| 一区二区日韩欧美中文字幕 | 国产精品 国内视频| 热re99久久精品国产66热6| 亚洲av成人精品一二三区| 精品久久久噜噜| 国产亚洲最大av| 免费大片黄手机在线观看| 一区二区三区乱码不卡18| 国产免费一级a男人的天堂| 亚洲欧美一区二区三区国产| 99热网站在线观看| 精品人妻熟女毛片av久久网站| 国产精品一区www在线观看| 日日撸夜夜添| 亚洲国产日韩一区二区| 午夜激情福利司机影院| 黑人高潮一二区| 久久久久久久久久久免费av| 国产精品一二三区在线看| 国产成人精品婷婷| 久久久久久久久久人人人人人人| 亚洲成人一二三区av| 97在线视频观看| 精品国产一区二区三区久久久樱花|