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

    縊蟶耐高鹽性能及其與生長性狀的相關(guān)性分析

    2022-11-26 05:59:10陳義華牛東紅董志國李家樂
    海洋科學(xué) 2022年10期
    關(guān)鍵詞:縊蟶耐鹽性家系

    陳義華, 牛東紅, 董志國, 李家樂

    縊蟶耐高鹽性能及其與生長性狀的相關(guān)性分析

    陳義華1, 3, 牛東紅1, 2, 董志國3, 李家樂1, 2

    (1. 上海海洋大學(xué) 上海水產(chǎn)養(yǎng)殖動(dòng)物良種培育與綠色養(yǎng)殖協(xié)同創(chuàng)新中心, 上海 201306; 2. 上海海洋大學(xué) 水產(chǎn)種質(zhì)資源發(fā)掘與利用教育部重點(diǎn)實(shí)驗(yàn)室, 上海 201306; 3. 江蘇海洋大學(xué) 江蘇省海洋生物技術(shù)重點(diǎn)實(shí)驗(yàn)室,江蘇 連云港 222005)

    縊蟶()是我國重要的水產(chǎn)經(jīng)濟(jì)貝類, 開展耐高鹽新品種選育對(duì)平衡不同地區(qū)縊蟶良種養(yǎng)殖意義重大。本研究對(duì)不同家系縊蟶幼貝進(jìn)行高鹽脅迫, 運(yùn)用統(tǒng)計(jì)學(xué)方法和相關(guān)性分析探究了縊蟶耐高鹽性能及其與生長性狀的相關(guān)性。結(jié)果顯示, 縊蟶家系間殼長、殼高、殼寬和體質(zhì)量表型值差異較大, 體質(zhì)量變異系數(shù)達(dá)15.53%。不同家系對(duì)高鹽耐受性差異較大, 家系間死亡率范圍在0~100%。脅迫實(shí)驗(yàn)共持續(xù)9 d, 各家系縊蟶首粒死亡時(shí)間在第1~9 d之間, 而且死亡率較高的家系內(nèi)縊蟶死亡時(shí)間離散程度較大, 說明不同個(gè)體間高鹽耐受性差異也較大。縊蟶開始出現(xiàn)大量死亡的時(shí)間主要在脅迫后第8 d, 第9 d部分家系死亡率高于80%, 說明縊蟶對(duì)高鹽耐受存在一個(gè)閾值, 超過該閾值會(huì)開始或大量死亡。皮爾遜和斯皮爾曼相關(guān)性分析顯示, 縊蟶存活時(shí)間和生長性狀為正相關(guān), 其中存活時(shí)間與縊蟶殼長、體質(zhì)量相關(guān)性均達(dá)到顯著水平(<0.05)??O蟶存活率與體質(zhì)量、生長速度均顯著正相關(guān)(<0.05)。整體而言, 殼長較長和體質(zhì)量較重的縊蟶個(gè)體耐高鹽性能較強(qiáng), 在高鹽環(huán)境中具有生長優(yōu)勢。本研究揭示了縊蟶的高鹽耐受性及其與生長性狀的相關(guān)性, 為縊蟶耐高鹽選育工作提供了重要的參考依據(jù)。

    縊蟶(); 耐鹽性; 生長性能; 相關(guān)性分析

    縊蟶()是雙殼綱, 竹蟶科的一種軟體動(dòng)物, 因其殼體表面有一條形似被繩索勒過的痕跡而得名[1]??O蟶具有生長快, 養(yǎng)殖周期短等優(yōu)點(diǎn), 一年就可達(dá)到商品規(guī)格[2], 在浙江、福建和江蘇等省份養(yǎng)殖廣泛, 是我國重要的水產(chǎn)經(jīng)濟(jì)物種。

    鹽度是影響海水貝類生長發(fā)育重要的環(huán)境因子之一。不同海水貝類最適生長鹽度范圍不同, 例如縊蟶在15~20鹽度下生長較快, 而波紋巴菲蛤()適宜鹽度范圍為25~30[3]。正常海水的平均鹽度為35, 由于氣候和內(nèi)陸徑流等原因, 河口及周圍地區(qū)海水鹽度變動(dòng)較大[4]。當(dāng)水體鹽度超過海洋貝類適應(yīng)范圍將會(huì)對(duì)其生長、存活造成巨大壓力[5-6]。高鹽對(duì)縊蟶生長、存活等具有一定影響, 研究發(fā)現(xiàn)鹽度30下縊蟶存活率和特定生長率均顯著低于低鹽度組[7]。在我國北方沿海和西北地區(qū)存在大量的高鹽水體, 例如山東鹽場有面積約14萬hm2的鹽度在30以上的水體[8]。而在江蘇省連云港部分養(yǎng)殖池塘, 干旱時(shí)期水體鹽度最高可達(dá)到28以上。生長性能是影響縊蟶產(chǎn)量的重要因素, 選育適合高鹽水域生長的縊蟶新品種對(duì)平衡南北區(qū)縊蟶養(yǎng)殖差異具有重要意義。

    據(jù)早期研究發(fā)現(xiàn), 縊蟶重要?dú)んw指標(biāo)和體質(zhì)量等生長性狀在高鹽下的遺傳力為中高水平, 可以通過選育來提升高鹽水域中縊蟶的生長性能[9]。當(dāng)把耐鹽、耐溫、抗病等抗逆性狀視作敏感與否時(shí), 則稱為閾性狀; 以存活時(shí)間評(píng)判時(shí), 則為數(shù)量性狀[10]。上述二者都無法從表型直接判斷其耐鹽性, 但是表型性狀較為直觀。關(guān)于貝類耐受性的研究較為多見, 例如縊蟶[11]、蝦夷扇貝()[12]的高溫耐受性, 香港牡蠣()的高鹽耐受性[13], 菲律賓蛤仔()的低氧耐受性[14]等, 但是關(guān)于貝類耐受性和生長性狀相關(guān)性的研究較為少見。本研究首次通過高鹽脅迫實(shí)驗(yàn)評(píng)估了縊蟶不同家系耐鹽性, 在此基礎(chǔ)上分析了縊蟶耐鹽性能和生長性狀的相關(guān)性。本研究結(jié)果將為縊蟶耐鹽新品種選育奠定基礎(chǔ)。

    1 材料與方法

    1.1 材料

    實(shí)驗(yàn)用對(duì)象為本實(shí)驗(yàn)于2020年10月建立的縊蟶家系, 繁育基地位于浙江省臺(tái)州市三門東航水產(chǎn)育苗科技有限公司。

    1.2 縊蟶家系的建立

    縊蟶繁育用親本殼長(5.31±0.93) cm, 體質(zhì)量(11.45±1.07) g。采用陰干加流水刺激的方法進(jìn)行催產(chǎn), 先將親本均勻置于產(chǎn)床上, 黑暗環(huán)境下陰干6~8 h, 后將產(chǎn)床直接轉(zhuǎn)移至育苗池中, 底部曝氣模擬水流。海水鹽度12, 水溫26 ℃。待親本排精、排卵時(shí), 將其拿出分別放入預(yù)先裝好同鹽度過濾海水的燒杯中。采用一雄配一雌的配對(duì)方法, 建立了30個(gè)全同胞家系。

    1.3 家系培育和池塘養(yǎng)殖

    幼蟲培育方法參照已有研究[11, 15], 并稍作調(diào)整??O蟶不同家系分別養(yǎng)殖在200 L的長方形塑料箱內(nèi), 餌料為育苗場培養(yǎng)的球等邊金藻()和牟氏角毛藻()。D形幼蟲及以前的時(shí)期投喂球等邊金藻, 投喂次數(shù)根據(jù)水體藻色而定, 每天換水量50%。后期投喂牟氏角毛藻, 每天換水量100%, 換水時(shí)直接加藻水, 視藻水濃度調(diào)整每天換水次數(shù), 期間嚴(yán)格控制家系間交叉混雜。

    經(jīng)過約2個(gè)月室內(nèi)養(yǎng)殖, 當(dāng)年12月將縊蟶全部家系轉(zhuǎn)移至育苗場自有室外池塘進(jìn)行中間培育。不同家系間以圍網(wǎng)隔離, 標(biāo)記位置并記錄家系號(hào)。2021年4月將縊蟶家系挖出, 轉(zhuǎn)移至連云港養(yǎng)殖池塘進(jìn)行養(yǎng)殖, 池塘海水周年平均鹽度24。養(yǎng)殖實(shí)驗(yàn)選用17個(gè)縊蟶家系, 初始?xì)らL(2.67±0.38) cm, 體質(zhì)量(1.48±0.26) g??O蟶家系采用圍網(wǎng)養(yǎng)殖(1.2 m× 1 m), 家系內(nèi)縊蟶個(gè)數(shù)為200~300粒。

    1.4 高鹽脅迫實(shí)驗(yàn)

    池塘養(yǎng)殖2個(gè)月后, 采集縊蟶樣本, 除其中一個(gè)家系縊蟶樣本數(shù)為45粒外, 其余16個(gè)家系縊蟶均采集60粒。縊蟶由池塘運(yùn)至實(shí)驗(yàn)室后, 在長方形塑料箱(0.45 m×0.35 m×0.15 m)中暫養(yǎng)3 d, 然后進(jìn)行高鹽條件下的脅迫實(shí)驗(yàn)。根據(jù)縊蟶半致死鹽度為35.8[16], 本實(shí)驗(yàn)設(shè)置脅迫鹽度為36。實(shí)驗(yàn)用高鹽水體用海水晶配置, 期間每天100%換水, 餌料為小球藻濃縮液(6×109~8×109ind./mL)。脅迫前期縊蟶正常攝食, 水體小球藻濃度(4×105~5×105)ind./mL。由于長期暴露于高鹽環(huán)境, 縊蟶后期攝食能力下降, 此時(shí)視攝食情況投喂。實(shí)驗(yàn)開始后, 每隔1 h觀察1次, 及時(shí)將死亡縊蟶撈出, 擦干體表海水, 記錄殼長、殼高、殼寬、體質(zhì)量等生長性狀。死亡判定標(biāo)準(zhǔn)為: 縊蟶閉殼肌失去作用, 觸摸刺激后閉殼肌無反應(yīng), 雙殼無閉合動(dòng)作。當(dāng)其中一個(gè)家系內(nèi)縊蟶全部死亡后實(shí)驗(yàn)結(jié)束。

    1.5 數(shù)據(jù)處理

    數(shù)據(jù)處理用Excel 2016完成, 箱線圖用Sigma-Plot 12.5繪制, 皮爾遜相關(guān)系數(shù)(Pearson correlation coefficient, PCC)和斯皮爾曼相關(guān)系數(shù)(Spearman's rank correlation coefficient, SCC)由SPSS 24完成。<0.05, 表示差異顯著;<0.01, 表示差異極顯著。根據(jù)之前的研究結(jié)果[17], 長期生活在高鹽環(huán)境下的縊蟶有較強(qiáng)的耐鹽性, 即縊蟶對(duì)高鹽環(huán)境有適應(yīng)能力。因此為增加結(jié)果準(zhǔn)確性, 在觀察縊蟶耐鹽性和生長速度的相關(guān)性時(shí)取縊蟶最近1個(gè)月的特定增長率。文中所用相關(guān)數(shù)據(jù)的計(jì)算公式如下:

    存活率(%)=c/m×100, (1)

    死亡率(%)=1–存活率, (2)

    特定增長率(%/d)=(ln2–ln1)/×100[18], (3)

    變異系數(shù)(%)=(標(biāo)準(zhǔn)差/平均值)×100, (4)

    式中,c為實(shí)驗(yàn)結(jié)束時(shí)縊蟶各家系存活個(gè)數(shù),m為實(shí)驗(yàn)開始時(shí)各家系縊蟶總個(gè)數(shù);2為各家系縊蟶平均殼長, 殼高, 殼寬和體質(zhì)量,1為采樣前1個(gè)月縊蟶的平均殼長, 殼高, 殼寬和體質(zhì)量;為養(yǎng)殖時(shí)間, 30 d。

    2 結(jié)果

    2.1 縊蟶家系間生長性狀生物學(xué)數(shù)據(jù)分析

    經(jīng)過2個(gè)月養(yǎng)殖, 不同家系間縊蟶殼長、殼高、殼寬和體質(zhì)量表型值出現(xiàn)較大的差異(表1)。其中, 對(duì)于體質(zhì)量和殼長性狀, 最大值和最小值間的差異在生長性狀中較為突出。本次用于脅迫實(shí)驗(yàn)的縊蟶殼長為(3.759 7±0.152 7)cm, 體質(zhì)量為(3.299 4±0.512 3) g。根據(jù)變異系數(shù), 家系間縊蟶體質(zhì)量離散程度最大, 超過15%; 殼高的離散程度最小。

    表1 縊蟶家系間生長性狀生物學(xué)數(shù)據(jù)

    2.2 高鹽脅迫下縊蟶不同家系死亡情況

    本實(shí)驗(yàn)設(shè)定其中一個(gè)家系內(nèi)縊蟶死亡率達(dá)到100%時(shí)實(shí)驗(yàn)結(jié)束, 實(shí)驗(yàn)結(jié)束時(shí)仍有一個(gè)家系死亡率為0(圖1)。圖1a為不同家系縊蟶的死亡率, 高鹽脅迫后縊蟶不同家系間死亡率由低到高差異較大, 部分家系死亡率高于80%。剔除沒有統(tǒng)計(jì)學(xué)意義的兩個(gè)家系(圖1b), 死亡率低的家系由于死亡個(gè)數(shù)較少, 死亡時(shí)間離散程度較小。死亡率較高的家系, 其縊蟶死亡時(shí)間離散程度較大, 死亡時(shí)間范圍較廣。

    圖1 高鹽脅迫下縊蟶各家系的死亡情況

    2.3 高鹽脅迫下縊蟶不同家系死亡時(shí)間

    自開始高鹽脅迫至其中一個(gè)家系縊蟶全部死亡, 脅迫歷時(shí)9 d(圖2)。實(shí)驗(yàn)期間, 不同家系縊蟶首粒死亡時(shí)間差異較大(圖2A)。結(jié)合死亡率分析表明, 縊蟶首粒死亡時(shí)間較晚的家系其整體死亡率較低, 說明該類家系內(nèi)縊蟶整體對(duì)高鹽有較強(qiáng)的耐受性。縊蟶各家系大量死亡暴發(fā)時(shí)間開始于第8 d, 第9 d死亡數(shù)量最多 (圖2b)。

    圖2 高鹽脅迫下縊蟶各家系死亡歷時(shí)

    2.4 縊蟶存活時(shí)間和生長性能相關(guān)性

    縊蟶存活時(shí)間和生長性狀的皮爾遜和斯皮爾曼相關(guān)性如表2所示。根據(jù)相關(guān)性分析理論, 結(jié)合各家系死亡數(shù)量, 對(duì)6個(gè)家系進(jìn)行統(tǒng)計(jì)分析。各家系縊蟶存活時(shí)間與殼長、殼高、殼寬、體質(zhì)量的皮爾遜和斯皮爾曼相關(guān)性除家系16外差異不大。家系12、15、16縊蟶存活時(shí)間與生長性狀的相關(guān)整體呈現(xiàn)負(fù)相關(guān), 家系13、14、17呈現(xiàn)正相關(guān)。將所有死亡縊蟶個(gè)體做一個(gè)整體統(tǒng)計(jì), 縊蟶存活時(shí)間與生長性狀的皮爾遜和斯皮爾曼相關(guān)性均是正相關(guān)關(guān)系。縊蟶存活率和殼長、體質(zhì)量的皮爾遜相關(guān)性達(dá)到顯著水平(<0.05), 而所有性狀的斯皮爾曼相關(guān)性達(dá)到了極顯著水平(<0.01)。

    表2 縊蟶存活時(shí)間和生長性狀的相關(guān)性

    注: * 代表<0.05; ** 代表<0.01。下同

    2.5 各家系縊蟶存活率和生長性能的相關(guān)性

    將各家系縊蟶存活率和生長指標(biāo)均值、增長率分別做相關(guān)性分析, 結(jié)果如表3和表4所示。各家系縊蟶存活率和生長性狀的皮爾遜相關(guān)性顯示, 存活率與殼長、殼高為負(fù)相關(guān), 與殼寬、體質(zhì)量為正相關(guān), 且與體質(zhì)量相關(guān)水平顯著(<0.05)。斯皮爾曼相關(guān)性顯示, 存活率僅與殼長為負(fù)相關(guān), 且與殼寬(<0.05)、體質(zhì)量(<0.01)相關(guān)性水平均顯著。各家系縊蟶存活率和生長性狀增長率的相關(guān)性方面, 綜合皮爾遜和斯皮爾曼相關(guān)結(jié)果, 相關(guān)性均為正相關(guān)且水平顯著(<0.05), 其中存活率與體質(zhì)量相關(guān)性達(dá)到極顯著水平(<0.01)。

    表3 各家系縊蟶存活率和生長性狀的相關(guān)性

    表4 各家系縊蟶存活率和生長性狀增長率的相關(guān)性

    3 討論

    近年來, 養(yǎng)殖縊蟶面臨品質(zhì)下降, 抗逆性差等眾多挑戰(zhàn)[19], 不同養(yǎng)殖區(qū)域水體環(huán)境的差異是導(dǎo)致縊蟶品質(zhì)不同的一個(gè)重要因素。選育縊蟶抗逆新品種對(duì)改善縊蟶養(yǎng)殖現(xiàn)狀, 消除縊蟶區(qū)域品質(zhì)和產(chǎn)量差異意義重大。類似耐溫、耐鹽等性狀選育中抗性指標(biāo)的判定是一個(gè)難點(diǎn), 即如何判斷哪個(gè)個(gè)體具有潛在抗性??O蟶的耐鹽性狀為典型的數(shù)量性狀范疇, 還沒有技術(shù)能通過殼體形態(tài)等比較直觀的外在性狀判斷。目前有關(guān)于水產(chǎn)動(dòng)物抗逆性判斷的方法大致有三種, 第一種為急性脅迫法, 主要思路是將大量受試個(gè)體直接放入極端的條件, 然后記錄死亡率, 先死亡個(gè)體定義沒有抗逆性, 后死亡或存活個(gè)體定義有抗逆性[20]; 第二種為緩慢脅迫法, 與急性脅迫法類似, 不同之處在于將直接改變生存環(huán)境變?yōu)榫徛淖兩姝h(huán)境[21]; 第三種就是一種基于心跳指標(biāo)的快速評(píng)估貝類抗性的方法[22]。但是以上三種方法在實(shí)際的育種中實(shí)現(xiàn)難度還較大, 現(xiàn)在主流育種應(yīng)用較多的還是將某一優(yōu)良性狀作為定向選育指標(biāo)的傳統(tǒng)選育方式。殼體指標(biāo)可以較直觀地反應(yīng)雙殼貝類的生長狀況[23], 因此研究高鹽養(yǎng)殖縊蟶生長性狀與耐鹽性的相關(guān)關(guān)系, 對(duì)高鹽選育親本的選擇以及高鹽養(yǎng)殖下縊蟶的生長性能評(píng)估具有非常重要的指導(dǎo)意義。

    縊蟶是廣鹽性貝類, 但是研究發(fā)現(xiàn)其濾水率、攝食率、同化率在鹽度5~30的范圍內(nèi)先上升后下降, 鹽度15~20是一個(gè)比較適宜的范圍[24]。高鹽環(huán)境會(huì)增加其滲透調(diào)節(jié)和生物代謝的耗能, 對(duì)其生長造成不利影響[7]。本研究中縊蟶首先在連云港池塘養(yǎng)殖兩個(gè)月, 該養(yǎng)殖水體相對(duì)最適養(yǎng)殖條件為高鹽環(huán)境, 對(duì)縊蟶有一定的選擇作用。選擇是育種的基礎(chǔ)[25], 表型變異是選擇的首要條件。本研究不同家系間縊蟶重要?dú)んw指標(biāo)和體質(zhì)量均出現(xiàn)不同程度的表型變異, 尤其是體質(zhì)量, 家系間表型值離散程度最大。而在貝類養(yǎng)殖中產(chǎn)量多以重量計(jì)算, 所以高鹽環(huán)境下提高縊蟶養(yǎng)殖產(chǎn)量, 以體質(zhì)量為主選性狀是可行且必要的。以死亡率判定水產(chǎn)動(dòng)物抗逆性的研究已有報(bào)道[26-27], 本研究縊蟶不同家系死亡率差異較大。實(shí)驗(yàn)結(jié)束后, 有的家系死亡率達(dá)到100%, 而有的家系卻無死亡個(gè)體出現(xiàn), 并且家系間第一只縊蟶死亡的時(shí)間差異也較大, 說明縊蟶不同家系對(duì)高鹽敏感性差異較大, 可以通過死亡率來判定縊蟶抗性家系和敏感家系, 進(jìn)而進(jìn)行家系間選育。家系內(nèi)選育是動(dòng)植物育種中遺傳改良的有效方法[28], 本研究中死亡率高的家系內(nèi)縊蟶死亡時(shí)間差別較大, 說明家系內(nèi)縊蟶不同個(gè)體間高鹽敏感性差異也較大, 可根據(jù)家系內(nèi)縊蟶個(gè)體的表現(xiàn)進(jìn)行選擇[29]。閾性狀是指它們的遺傳是由多基因決定的, 而它們的表型是非連續(xù)性的一類性狀[30]。閾性狀是一類重要的數(shù)量性狀, 動(dòng)植物包括人類在內(nèi)的抗病能力如“患病”或“正?!? “存活”或“死亡”等都可視為閾性狀。這類性狀存在一個(gè)閾值, 超過這一閾值才會(huì)顯示相應(yīng)的性狀[30]。在持續(xù)9 d的高鹽脅迫中, 縊蟶大量死亡時(shí)間爆發(fā)于第8 d。第8 d以前, 各家系縊蟶死亡緩慢, 死亡數(shù)量較少, 并且家系間首??O蟶死亡時(shí)間具有早晚差異, 說明縊蟶對(duì)高鹽耐受有一個(gè)閾值, 超過該閾值會(huì)開始或大量死亡。結(jié)合死亡情況和死亡時(shí)間來看, 不同家系的閾值不同, 不同個(gè)體的閾值也不同。

    性狀間相關(guān)性是反映物種各性狀間相關(guān)程度的重要參數(shù), 對(duì)指導(dǎo)育種早期選擇和間接選擇等策略的準(zhǔn)確制定具有重要意義[31]。間接選擇指利用一個(gè)性狀與另一個(gè)性狀的顯著相關(guān), 通過改變其中一個(gè)性狀來間接改變另一個(gè)性狀的選擇[32]。直接施加選擇的性狀叫副性狀, 要改良的目標(biāo)性狀叫主性狀。目前缺少判斷縊蟶個(gè)體耐鹽性的精確指標(biāo), 依照傳統(tǒng)的選育理論, 可進(jìn)行高鹽脅迫或高鹽養(yǎng)殖, 然后留取存活個(gè)體做親本繼續(xù)選育??紤]產(chǎn)量問題, 個(gè)體大、體質(zhì)量重的個(gè)體應(yīng)被優(yōu)先選擇當(dāng)做親本。因此, 探究縊蟶殼體指標(biāo)及體質(zhì)量與耐鹽性的相關(guān)性, 對(duì)于間接選擇至關(guān)重要。皮爾遜相關(guān)和斯皮爾曼相關(guān)是分析相關(guān)性的兩個(gè)常用相關(guān)指標(biāo)[33-34], 本研究在分析縊蟶存活時(shí)間和生長性狀相關(guān)性時(shí), 兩種指標(biāo)分析結(jié)果一致。不同家系, 其縊蟶存活時(shí)間與殼體指標(biāo)的相關(guān)性有正有負(fù), 主要原因可能是樣本量較小, 這一觀點(diǎn)可以從基于整體統(tǒng)計(jì)后的結(jié)果得到證實(shí)。通過對(duì)所有個(gè)體的數(shù)據(jù)分析表明, 縊蟶存活時(shí)間和所有生長性狀均為正相關(guān), 其中殼長、體質(zhì)量與存活時(shí)間為顯著正相關(guān)。首先, 殼長、殼高、殼寬和體質(zhì)量是數(shù)量性狀, 其數(shù)值是連續(xù)的, 存活時(shí)間也為連續(xù)的, 可以用皮爾遜相關(guān)性和斯皮爾曼相關(guān)性進(jìn)行分析。但是, 皮爾遜相關(guān)性分析是要求數(shù)值符合正態(tài)分布或者接近正態(tài)的單峰分布[35]。不同家系內(nèi)縊蟶雖然具有統(tǒng)計(jì)學(xué)意義, 但是數(shù)量相對(duì)較少, 同時(shí)存活時(shí)間跨度較大, 這可能是導(dǎo)致皮爾遜相關(guān)性和斯皮爾曼相關(guān)性結(jié)果差異的主要原因。本研究中參與統(tǒng)計(jì)的全部死亡縊蟶數(shù)量達(dá)到300粒以上, 誤差相對(duì)單個(gè)家系而言更小, 分析結(jié)果更加可靠。因此, 整體而言個(gè)體大的縊蟶具有較強(qiáng)的耐鹽性。分析縊蟶各個(gè)家系存活率和其生長性狀平均值的皮爾遜和斯皮爾曼相關(guān)性, 二者均顯示存活率和不同生長性狀的相關(guān)性有正有負(fù), 這可能是由家系內(nèi)縊蟶不同個(gè)體耐受性差異較大導(dǎo)致。但是皮爾遜和斯皮爾曼相關(guān)性顯示存活率和體質(zhì)量間為顯著正相關(guān), 說明體質(zhì)量重的家系高鹽耐受性整體較高。特定增長率可以校正分析數(shù)據(jù)的不等方差[18], 降低數(shù)據(jù)分析誤差。各家系縊蟶存活率和生長性狀增長率的皮爾遜和斯皮爾曼相關(guān)性結(jié)果一致, 都為正相關(guān), 說明生長速度快的縊蟶一般耐鹽性較強(qiáng)。張興志等[10]在香港牡蠣家系耐高鹽性篩選的研究發(fā)現(xiàn), 幼蟲期生長與存活率兩個(gè)性狀間不存在相關(guān)性, 與本結(jié)果不同, 這種差異可能與物種、生長階段、養(yǎng)殖環(huán)境等因素有關(guān)[26]。綜合分析, 殼長、體質(zhì)量、在高鹽環(huán)境中生長速度快的縊蟶個(gè)體具有較強(qiáng)的耐高鹽性能。

    4 結(jié)論

    本研究結(jié)果證明了縊蟶家系間以及個(gè)體間均存在高鹽耐受性差異, 說明了縊蟶耐鹽性是閾性狀。整體而言, 殼長較長和體質(zhì)量較重的縊蟶個(gè)體耐高鹽性能較強(qiáng), 在高鹽環(huán)境中具有生長優(yōu)勢。本研究為縊蟶耐高鹽性狀與生長性能的關(guān)系判定提供了有力證據(jù), 對(duì)進(jìn)一步開展縊蟶耐高鹽選育具有重要意義。

    [1] 常亞青. 貝類增養(yǎng)殖學(xué)[M]. 北京: 中國農(nóng)業(yè)出版社. 2007: 108-109.

    CHANG Yaqing. Science of shellfish culture and pro-pagation[M]. Beijing: China Agricultural Press. 2007: 108-109.

    [2] NIU D H, FENG B B, LIU D B, et al. Significant genetic differentiation among ten populations of the razor clamalong the coast of china revealed by a microsatellite analysis[J]. Zoological Studies, 2012, 51(3): 406-414.

    [3] 張鵬飛. 波紋巴非蛤生理生態(tài)學(xué)研究[D]. 廈門: 廈門大學(xué), 2018.

    ZHANG Pengfei. Study on physiological ecology of the undulated surf clam,[D]. Xiamen: Xiamen University, 2018.

    [4] 陳敏. 化學(xué)海洋學(xué)[M]. 北京: 海洋出版社, 2009: 1- 200.

    CHEN Min. Chemical oceanography[M]. Beijing: China Ocean Press, 2009: 1-200.

    [5] BURDON D, CALLAWAY R, ELLIOTT M, et al. Mass mortalities in bivalve populations: A review of the edible cockle[J]. Estuarine, Coastal and Shelf Science, 2014, 150: 271-280.

    [6] CARREGOSA V, FIGUEIRA E, GIL A M, et al. Tolerance ofto salinity: Osmotic and metabolic aspects[J]. Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology, 2014, 171: 36-43.

    [7] 丁紅兵, 李浩宇, 陳義華, 等. 高鹽對(duì)縊蟶生長存活、Na+/K+-ATPase活性及能量代謝相關(guān)指標(biāo)的影響[J]. 上海海洋大學(xué)學(xué)報(bào), 2021, 5: 1-17.

    DING Hongbing, LI Haoyu, CHEN Yihua, et al. Effects of high salinity on growth and survival, Na+/K+-ATPase activity and energy metabolism related indexes of razor clam[J]. Journal of Shanghai Ocean University, 2021, 5: 1-17.

    [8] 李娜, 趙玉超, 王仁杰, 等. 高鹽脅迫對(duì)凡納濱對(duì)蝦消化及免疫相關(guān)酶活力的影響[J]. 生態(tài)學(xué)報(bào), 2018, 38(4): 1411-1417.

    LI Na, ZHAO Yuchao, WANG Renjie, et al. Effects of high salinity on digestive and immunity-related enzy-mes in[J]. Acta Ecologica Sinica, 2018, 38(4): 1411-1417.

    [9] CHEN Y H, LI H Y, DING H B, et al. Heritability estimation and path analysis for growth traits of the razor clam Sinonovacula constricta under high salinity[J]. Aquaculture, 2021, 545: 737175.

    [10] 陳澤輝. 群體與數(shù)量遺傳學(xué)[M]. 貴陽: 貴州科技出版社. 2009: 123-125.

    CHEN Zehui. Population and quantitative genetics[M]. Guiyang: Guizhou Science and Technology Publishing House Co., Ltd. 2009: 123-125.

    [11] 孟德龍, 薛寶寶, 申奔龍, 等. 縊蟶抗高溫新品系的早期選育與篩選[J]. 大連海洋大學(xué)學(xué)報(bào), 2020, 35(1): 71-74.

    MENG Delong, XUE Baobao, SHEN Benlong, et al. Early selective breeding and screening of new high temperature resistant lines of razor clam Sinonovacula constrict[J]. Journal of Dalian Fisheries University, 2020, 35(1): 71-74.

    [12] 王慶志, 李石磊, 付成東, 等. 蝦夷扇貝耐高溫育種家系的建立與早期篩查[J]. 水產(chǎn)學(xué)報(bào), 2014, 38(3): 371-377.

    WANG Qingzhi, LI Shilei, FU Chengdong, et al. Estab-lishment of high temperature resistance families and use of laboratory assays to predict subsequent survival in juvenile stage of the Japanese scallop ()[J]. Journal of Fishries China, 2014, 38(3): 371-377.

    [13] 張興志, 官俊良, 何蘋萍, 等. 香港牡蠣野生群體家系建立及耐高鹽性篩選[J]. 南方農(nóng)業(yè)學(xué)報(bào), 2019, 50(2): 385-390.

    ZHANG Xingzhi, GUAN Junliang, HE Pingping, et al. Establishment of families and subsequent high salinity resistance testing for wild population of[J]. Journal of Southern Agriculture, 2019, 50(2): 385-390.

    [14] 王化敏, 丁鑒鋒, 楊東敏, 等. 4種殼色菲律賓蛤仔在低氧脅迫下的耐受能力比較研究[J]. 大連海洋大學(xué)學(xué)報(bào), 2018, 33(2): 181-189.

    WANG Huamin, DING Jianfeng, YANG Dongmin, et al. Comparison of hypoxia stress tolerance among Manila clamwith four shell-colors[J]. Journal of Dalian Fisheries University, 2018, 33(2): 181-189.

    [15] NIU D H DU Y C, WANG Z, et al. Construction of the first high-density genetic linkage map and analysis of quantitative trait loci for growth-related traits in[J]. Marine Biotechnology, 2017, 19(5): 488-496.

    [16] 李煉星, 杜文俊, 王成東, 等. 縊蟶家系生長和耐熱、耐高鹽性能的對(duì)比研究[J]. 上海海洋大學(xué)學(xué)報(bào), 2016(4): 515-521.

    LI Lianxing, DU Wenjun, WANG Chengdong, et al. Comparative analysis of growth and heat tolerance, salt tolerance traits amongfamilies[J]. Journal of Dalian Fisheries University, 2016(4): 515-521.

    [17] CHEN Y H, YE B, NIU D H, et al. Changes in metabolism and immunity in response to acute salinity stress in Chinese razor clams from different regions[J]. Aqua-culture Reports, 2021, 19(6): 100624.

    [18] MAS-MU?OZ J, BLONK R, SCHRAMA J W, et al. Genotype by environment interaction for growth of sole () reared in an intensive aquaculture system and in a semi-natural environment[J]. Aquaculture, 2013, 410-411: 230-235.

    [19] RAN Z S, LI Z Z, YAN X J, et al. Chromosome-level genome assembly of the razor clam(Lamarck, 1818)[J]. Molecular Ecology Resources, 2019, 19(6): 1647-1658.

    [20] 王金鵬. 牡蠣生長與高溫耐受性狀的遺傳解析[D]. 北京: 中國科學(xué)院大學(xué), 2016.

    WANG Jinpeng. Genetic architecture dissection of growth-related and thermo-tolerant traits in oysters[D]. Beijing: University of Chinese Academy of Sciences, 2016.

    [21] 劉寶鎖, 張?zhí)鞎r(shí), 孔杰, 等. 大菱鲆生長和耐高溫性狀的遺傳參數(shù)估計(jì)[J]. 水產(chǎn)學(xué)報(bào), 2011, 35(11): 1601- 1606.

    LIU Baosuo, ZHAN Tianshi, KONG Jie, et al. Estimation of genetic parameters for growth and upper thermal tolerance traits in turbot[J]. Journal of Fisheries of China, 2011, 35(11): 1601-1606.

    [22] 林思恒, 吳富村, 張國范. 高溫脅迫下皺紋盤鮑不同養(yǎng)殖群體心率變化比較[J]. 海洋科學(xué), 2016, 40(11): 84-90.

    LIN Siheng, WU Fucun, ZHANG Guofan. Variation in cardiac response to thermal stress in two different cultured populations of[J]. Marine Scien-ces, 2016, 40(11): 84-90.

    [23] 曹學(xué)順, 丁君, 常亞青, 等. 加州扁鳥蛤()數(shù)量性狀的相關(guān)性及通徑分析[J]. 中國農(nóng)業(yè)科技導(dǎo)報(bào), 2014, 16(5): 139-146.

    CAO Xueshun, DING Jun, CHANG Yaqing, et al. Correlation and Path Analysis on Quantitative Traits of[J]. Journal of Agricultural Science and Technology, 2014, 16(5): 139-146.

    [24] 呂昊澤, 劉健, 陳錦輝, 等. 鹽度對(duì)長江口3種濾食性貝類濾水率、攝食率、同化率的影響[J]. 海洋科學(xué), 2016, 40(8): 10-17.

    LU Haoze, LIU Jian, CHEN Jinhui, et al. Effects of salinity on filtration, ingestion, and assimilation rates of three filter-feeding bivalves in the Yangtze River estuary[J]. Marine Sciences, 2016, 40(8): 10-17.

    [25] 張曉娟, 周莉, 桂建芳. 遺傳育種生物技術(shù)創(chuàng)新與水產(chǎn)養(yǎng)殖綠色發(fā)展[J]. 中國科學(xué): 生命科學(xué), 2019(11): 1409-1429.

    ZHANG Xiaojuan, ZHOU Li, GUI Jianfang. Biotech-nological innovation in genetic breeding and sustainable green development in Chinese aquaculture[J]. Scientia Sinica(Vitae), 2019(11): 1409-1429.

    [26] 馬愛軍, 黃智慧, 王新安, 等. 大菱鲆()耐高溫品系選育及耐溫性能評(píng)估[J]. 海洋與湖沼, 2012, 43(4): 797-804.

    MA Aijun, HUANG Zhihui, WANG Xinan, et al. The selective breeding of thermal tolerance family and appraisal of performance in turbot[J]. Oceanologia et Limnologia Sinica, 2012, 43(4): 797-804.

    [27] 楊淞, 韋其鋒, 嚴(yán)太明, 等. 孔雀魚、月光魚耐溫限度的初步研究[J]. 四川農(nóng)業(yè)大學(xué)學(xué)報(bào), 2009, 27(1): 106-110.

    YANG Song, WEI Qifeng, YAN Taiming, et al. The preliminary study on temperature resistance ofand[J]. Journal Of Sichuan Agricultural University, 2009, 27(1): 106-110.

    [28] 桑士田, 閆喜武, 霍忠明, 等. 家系內(nèi)大、小兩種規(guī)格菲律賓蛤仔的雙列雜交[J]. 水產(chǎn)學(xué)報(bào), 2012, 36(6): 832-837.

    SANG Shitian, YAN Xiwu, HUO Zhongming, et al. The diallel cross between large and small sizes within family of Manila clam ()[J]. Journal Of Fishries China, 2012, 36(6): 832-837.

    [29] 戶國, 谷偉, 白慶利, 等. 鮭鱒魚類育種中常用策略、方法及其應(yīng)用概述[J]. 水產(chǎn)學(xué)雜志, 2014, 27(6): 66-72.

    HU Guo, GU Wei, BAI Qingli, et al. An overview of common strategies, methods and applications in salmon and trout breeding[J]. Chinese Journal of Fisheries, 2014, 27(6): 66-72.

    [30] 戴灼華, 王亞馥, 粟翼玟. 遺傳學(xué).第2版[M]. 北京: 高等教育出版社. 2008: 215-220.

    DAI Zhuohua, WANG Yafu, LI Yiwen. Genetics[M]. Beijing: Higher Education Press, 2008: 215-220.

    [31] 欒生, 孔杰, 張?zhí)鞎r(shí), 等. 基于表型值和育種值的中國對(duì)蝦生長、抗逆性狀相關(guān)分析[J]. 漁業(yè)科學(xué)進(jìn)展, 2008, 29(3): 14-20.

    LUAN Sheng, KONG Jie, ZHANG Tianshi, et al. Phenotype and breeding value based correlation analyses of arowth and stress-resistance traits in[J]. Progress in Fishery Sciences, 2008, 29(3): 14-20.

    [32] 李加納. 數(shù)量遺傳學(xué)概論[M]. 重慶: 西南師范大學(xué)出版社. 2007.

    LI Jiana. An outline of quantitative genetics[M]. Chongqing: Southwest Normal University Press. 2007.

    [33] WANG B, WANG R D, WANG Y M. Compatible matrices of Spearman’s rank correlation[J]. Statistics & Probability Letters, 2019, 151: 67-72.

    [34] 陳敏嫣, 水珊珊, 周曉嬌, 等. 沿海紫貽貝體內(nèi)重金屬含量及相關(guān)性分析[J]. 環(huán)境科學(xué)與技術(shù), 2020, 43(6): 24-29.

    CHEN Minyan, SHUI Shanshan, ZHOU Xiaojiao, et al. Analysis on the content and correlation of heavy metals in[J]. Environmental Science & Technology, 2020, 43(6): 24-29.

    [35] 陳平雁, 黃浙明. IBM SPSS 19統(tǒng)計(jì)軟件應(yīng)用教程[M]. 北京: 人民衛(wèi)生出版社. 2012.

    CHEN Pingyan, HUANG Zheming. IBM SPSS 19 statistical software application tutorial[M]. Beijing: People’s Medical Publishing House, 2012.

    Correlation between high salt tolerance and growth perfor-mance of the razor clam

    CHEN Yi-hua1, 3, NIU Dong-hong1, 2, DONG Zhi-guo3, LI Jia-le1, 2

    (1. Shanghai Collaborative Innovation Center for Cultivating Elite Breeds and Green-culture of Aquaculture animals, Shanghai 201306, China; 2. Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Ministry of Education, Shanghai Ocean University, Shanghai 201306, China; 3. Jiangsu Key Laboratory of Marine Biotechnology, Jiangsu Ocean University, Lianyungang 222005, China)

    The razor clam,, is an aquatic species of economic importance in China. Breeding for a new razor clam variety with improved salt tolerance can balance the yield of clams in different regions. In this study, the young razor clams from different families with a 2.67±0.38-cm mean shell length and 1.48±0.26-g mean body weight were used for the high salt stress experiment. Data were tested for correlation between high salt tolerance and growth performance of razor clams. The results showed that there was a wide variation in the phenotypic values of shell length, shell height, shell width, and body weight among different families. The coefficient of variation of body weight was 15.53%. The mortality among different families ranged from 0% to 100%, indicating that tolerance to high salt stress differed among the families. The family with the highest mortality had a large degree of variation in time to death, which indicated that there were great differences in high salt tolerance among individuals within the family. The stress experiment lasted for 9 days, and the razor clams died in large numbers on the 8th day. This showed that there was a threshold for high salt tolerance in razor clams. Pearson and Spearman correlation analysis showed that the survival time of the razor clam was significantly and positively correlated with shell length and body weight (<0.05). The survival rate of razor clams in each family was significantly and positively correlated with body weight and specific growth rate (<0.05). Overall, razor clams with longer shell lengths and heavier body weights were more tolerant to high salt stress and had a higher growth rate in the high salt environment. This study provides a baseline for the selection and breeding of the razor clam for high salt tolerance.

    ; high salt tolerance; growth performance; correlation analysis

    Sep. 28, 2021

    S968

    A

    1000-3096(2022)10-0051-08

    10.11759/hykx20210928004

    2021-09-28;

    2022-03-02

    國家重點(diǎn)研發(fā)計(jì)劃“藍(lán)色糧倉科技創(chuàng)新”專項(xiàng)項(xiàng)目(2019YFD0900400)

    [National Key R&D Plan “Blue Granary Science and Technology Innovation” special project, No. 2019YFD0900400]

    陳義華 (1993—), 男, 山東濟(jì)南人, 博士研究生, 主要從事海洋貝類遺傳育種研究, E-mail: cyihua@163.com; 牛東紅,通信作者, E-mail: dhniu@shou.edu.cn

    (本文編輯: 康亦兼)

    猜你喜歡
    縊蟶耐鹽性家系
    蘇北沿??O蟶養(yǎng)殖死亡的原因分析及防控措施
    沿??O蟶養(yǎng)殖春季肥水技術(shù)
    郁金香耐鹽性鑒定方法研究
    源于大麥小孢子突變體的苗期耐鹽性研究
    三個(gè)大豆品種萌發(fā)期和苗期的耐鹽性比較
    馬氏珠母貝紅色閉殼肌F1代的家系選育及家系評(píng)定
    一個(gè)非綜合征型聾家系的分子病因?qū)W研究
    肝豆?fàn)詈俗冃?個(gè)家系的基因突變分析
    寧波近岸海域縊蟶有機(jī)氯農(nóng)藥風(fēng)險(xiǎn)評(píng)估
    兩個(gè)雄激素不敏感綜合征家系中AR基因突變檢測
    听说在线观看完整版免费高清| 丝瓜视频免费看黄片| 亚洲伊人久久精品综合| 街头女战士在线观看网站| 91午夜精品亚洲一区二区三区| 真实男女啪啪啪动态图| 高清欧美精品videossex| 亚洲高清免费不卡视频| 99热全是精品| 男女无遮挡免费网站观看| 久久久精品免费免费高清| 精品酒店卫生间| 国产爱豆传媒在线观看| 精品国产三级普通话版| 久久久成人免费电影| 哪个播放器可以免费观看大片| 我的女老师完整版在线观看| 久久精品国产亚洲av涩爱| 高清欧美精品videossex| 天天躁日日操中文字幕| 国产黄片视频在线免费观看| 女人十人毛片免费观看3o分钟| 国产成人精品婷婷| 日本猛色少妇xxxxx猛交久久| 日韩伦理黄色片| 成人黄色视频免费在线看| 精品人妻视频免费看| 哪个播放器可以免费观看大片| av在线观看视频网站免费| 国产高潮美女av| 22中文网久久字幕| 熟女电影av网| 欧美高清成人免费视频www| 亚洲欧美日韩另类电影网站 | 国产日韩欧美亚洲二区| 乱码一卡2卡4卡精品| videos熟女内射| 久久久a久久爽久久v久久| 国产成人精品一,二区| 日韩欧美一区视频在线观看 | 国产黄片视频在线免费观看| 久久韩国三级中文字幕| 最后的刺客免费高清国语| 美女cb高潮喷水在线观看| 国产白丝娇喘喷水9色精品| 中文字幕久久专区| 国产综合懂色| 网址你懂的国产日韩在线| 91aial.com中文字幕在线观看| 在线观看三级黄色| 中文乱码字字幕精品一区二区三区| 亚洲久久久久久中文字幕| 天天一区二区日本电影三级| 3wmmmm亚洲av在线观看| av卡一久久| 欧美成人精品欧美一级黄| 国产精品不卡视频一区二区| 国产午夜精品久久久久久一区二区三区| 一本色道久久久久久精品综合| av国产精品久久久久影院| 乱系列少妇在线播放| 亚洲三级黄色毛片| 欧美zozozo另类| 少妇人妻一区二区三区视频| av在线亚洲专区| 国产精品福利在线免费观看| 午夜福利在线在线| 一级毛片黄色毛片免费观看视频| av.在线天堂| 女人久久www免费人成看片| 欧美老熟妇乱子伦牲交| 成人高潮视频无遮挡免费网站| 国产成人午夜福利电影在线观看| 国产伦精品一区二区三区四那| 国产精品一及| 日本wwww免费看| 日韩av不卡免费在线播放| 欧美xxxx性猛交bbbb| 一级黄片播放器| 麻豆成人午夜福利视频| 熟女人妻精品中文字幕| 老司机影院成人| 免费观看的影片在线观看| 黄色视频在线播放观看不卡| 久久久久国产精品人妻一区二区| 七月丁香在线播放| 国产成人91sexporn| www.色视频.com| 日本免费在线观看一区| 精品酒店卫生间| 久久久久久久久久久免费av| 国产成人freesex在线| 亚洲四区av| 日韩免费高清中文字幕av| 男插女下体视频免费在线播放| 午夜亚洲福利在线播放| 亚洲aⅴ乱码一区二区在线播放| 少妇裸体淫交视频免费看高清| 免费av不卡在线播放| 欧美最新免费一区二区三区| 久久热精品热| 直男gayav资源| 国产综合懂色| av在线老鸭窝| 伊人久久国产一区二区| av线在线观看网站| 中文在线观看免费www的网站| 美女高潮的动态| 一区二区三区免费毛片| 国产精品爽爽va在线观看网站| 少妇被粗大猛烈的视频| av福利片在线观看| 精品人妻一区二区三区麻豆| 日本三级黄在线观看| 极品少妇高潮喷水抽搐| 高清毛片免费看| 成年av动漫网址| 亚洲不卡免费看| 欧美3d第一页| 中文乱码字字幕精品一区二区三区| 久久6这里有精品| 97热精品久久久久久| 久久久精品94久久精品| 少妇人妻久久综合中文| 欧美日韩视频精品一区| 国产成人精品久久久久久| 91精品一卡2卡3卡4卡| 成人一区二区视频在线观看| 青春草视频在线免费观看| 日韩一区二区视频免费看| av线在线观看网站| 亚洲成人一二三区av| 久久久久久久亚洲中文字幕| 午夜精品一区二区三区免费看| 亚洲人成网站高清观看| 成人亚洲精品av一区二区| 少妇 在线观看| 亚洲av福利一区| 国产视频首页在线观看| 亚洲色图综合在线观看| 男女国产视频网站| 欧美精品国产亚洲| 国产精品一及| 综合色av麻豆| av在线天堂中文字幕| 久久久久久国产a免费观看| 国产精品精品国产色婷婷| 久久ye,这里只有精品| 亚洲成人一二三区av| 成人亚洲精品一区在线观看 | 最近中文字幕2019免费版| 18禁裸乳无遮挡动漫免费视频 | 日韩av不卡免费在线播放| 深爱激情五月婷婷| 日日啪夜夜爽| 久久久午夜欧美精品| 观看美女的网站| 蜜桃久久精品国产亚洲av| 精品熟女少妇av免费看| 一级毛片久久久久久久久女| 真实男女啪啪啪动态图| 深夜a级毛片| 嫩草影院入口| 国产视频内射| 亚洲国产成人一精品久久久| 亚洲精品第二区| 国产在视频线精品| 亚洲精品日本国产第一区| 日本av手机在线免费观看| 久久久久久久精品精品| 视频中文字幕在线观看| 2018国产大陆天天弄谢| 精品一区二区免费观看| 亚洲精品日韩在线中文字幕| 老师上课跳d突然被开到最大视频| 亚洲av电影在线观看一区二区三区 | 国产高清三级在线| 国产亚洲精品久久久com| 亚洲经典国产精华液单| 成人一区二区视频在线观看| 久久久a久久爽久久v久久| 国产免费又黄又爽又色| 乱码一卡2卡4卡精品| 亚洲精品456在线播放app| 一级av片app| 久久精品国产亚洲网站| 欧美激情在线99| 夫妻午夜视频| 黑人高潮一二区| 波野结衣二区三区在线| 大话2 男鬼变身卡| 色吧在线观看| 亚洲国产最新在线播放| 久久精品久久久久久久性| 国产精品不卡视频一区二区| 亚洲精品日韩在线中文字幕| 色吧在线观看| 亚洲,一卡二卡三卡| 日本一二三区视频观看| 五月玫瑰六月丁香| 女人十人毛片免费观看3o分钟| 国产成人a∨麻豆精品| 久久99热6这里只有精品| 亚洲av中文字字幕乱码综合| 麻豆成人av视频| 在现免费观看毛片| 在线观看免费高清a一片| 久久久久国产精品人妻一区二区| 中文字幕亚洲精品专区| 人妻一区二区av| av在线观看视频网站免费| 精品一区二区免费观看| 校园人妻丝袜中文字幕| 水蜜桃什么品种好| 国产淫语在线视频| 国产成人精品一,二区| 成人欧美大片| 免费大片黄手机在线观看| 日韩电影二区| 成人毛片60女人毛片免费| 亚洲国产欧美在线一区| 成年人午夜在线观看视频| 久久综合国产亚洲精品| 中国国产av一级| 麻豆精品久久久久久蜜桃| 国产日韩欧美在线精品| 黄色一级大片看看| 亚洲真实伦在线观看| 97在线视频观看| 一本久久精品| 如何舔出高潮| 不卡视频在线观看欧美| 免费看a级黄色片| 日韩强制内射视频| 少妇人妻精品综合一区二区| 精品久久久噜噜| 老司机影院毛片| 97超视频在线观看视频| 亚洲精品国产av成人精品| 精品亚洲乱码少妇综合久久| 国产成人freesex在线| 亚洲美女视频黄频| 热99国产精品久久久久久7| 可以在线观看毛片的网站| 欧美精品人与动牲交sv欧美| 搡老乐熟女国产| 91精品国产九色| 岛国毛片在线播放| 亚洲成人中文字幕在线播放| 在线观看国产h片| 黄色怎么调成土黄色| 国产一区亚洲一区在线观看| 2022亚洲国产成人精品| 精华霜和精华液先用哪个| 爱豆传媒免费全集在线观看| 久久久久久久大尺度免费视频| 亚洲av不卡在线观看| 欧美成人a在线观看| 性色avwww在线观看| 国产精品熟女久久久久浪| 日本一本二区三区精品| 亚洲精品亚洲一区二区| 国产av国产精品国产| 欧美日韩亚洲高清精品| 一级毛片黄色毛片免费观看视频| 成人国产av品久久久| 精品亚洲乱码少妇综合久久| 日韩欧美 国产精品| 亚洲国产欧美在线一区| 一级黄片播放器| 国产免费福利视频在线观看| 一级毛片我不卡| 日韩国内少妇激情av| 亚洲国产精品专区欧美| 国产黄片视频在线免费观看| 在线观看三级黄色| 免费观看无遮挡的男女| 久久精品国产鲁丝片午夜精品| 欧美激情国产日韩精品一区| 亚洲国产欧美人成| 日韩视频在线欧美| 精品午夜福利在线看| 国产精品久久久久久久久免| 神马国产精品三级电影在线观看| 2021少妇久久久久久久久久久| 久久亚洲国产成人精品v| 亚洲国产欧美在线一区| 国产大屁股一区二区在线视频| 人妻制服诱惑在线中文字幕| 亚洲欧美成人综合另类久久久| 亚洲av中文字字幕乱码综合| 夜夜看夜夜爽夜夜摸| 精品久久久久久久久av| 亚洲精品第二区| 国产探花在线观看一区二区| 亚洲成人精品中文字幕电影| 一区二区三区乱码不卡18| 99久久精品国产国产毛片| 亚洲丝袜综合中文字幕| 国产精品国产三级专区第一集| 国产男人的电影天堂91| 国产精品99久久99久久久不卡 | 国产 一区 欧美 日韩| 熟女人妻精品中文字幕| 国产一区有黄有色的免费视频| 80岁老熟妇乱子伦牲交| 国产黄频视频在线观看| 网址你懂的国产日韩在线| 欧美区成人在线视频| www.色视频.com| 亚洲高清免费不卡视频| 成人二区视频| 日韩国内少妇激情av| 丰满乱子伦码专区| 欧美极品一区二区三区四区| 久久99热6这里只有精品| 大片免费播放器 马上看| 中文字幕制服av| 波野结衣二区三区在线| 亚洲国产av新网站| 国产精品三级大全| 亚洲av免费在线观看| 美女内射精品一级片tv| 亚洲美女搞黄在线观看| 久久久久久久久久人人人人人人| 国产成人91sexporn| 晚上一个人看的免费电影| 国产极品天堂在线| 免费在线观看成人毛片| 久久久亚洲精品成人影院| 亚洲国产色片| 成人毛片a级毛片在线播放| 日本欧美国产在线视频| 国产免费一级a男人的天堂| 美女cb高潮喷水在线观看| 伊人久久精品亚洲午夜| 97超视频在线观看视频| 男人和女人高潮做爰伦理| 一级毛片aaaaaa免费看小| 我的老师免费观看完整版| 大片电影免费在线观看免费| 美女视频免费永久观看网站| 亚洲精品亚洲一区二区| 校园人妻丝袜中文字幕| 禁无遮挡网站| tube8黄色片| 成年女人看的毛片在线观看| 我要看日韩黄色一级片| 国产色爽女视频免费观看| 边亲边吃奶的免费视频| 久久人人爽av亚洲精品天堂 | 精品久久久精品久久久| 国产亚洲精品久久久com| 国产欧美日韩精品一区二区| 午夜日本视频在线| 久久99热这里只频精品6学生| 色视频www国产| 三级经典国产精品| 91久久精品国产一区二区成人| 精品视频人人做人人爽| 欧美日韩在线观看h| 综合色丁香网| 亚洲三级黄色毛片| 美女主播在线视频| av国产精品久久久久影院| 久久精品国产a三级三级三级| 亚洲av日韩在线播放| 久久久国产一区二区| 2022亚洲国产成人精品| 视频中文字幕在线观看| 少妇熟女欧美另类| av卡一久久| 亚洲精品国产av蜜桃| 女人久久www免费人成看片| 国产精品人妻久久久影院| 国产精品久久久久久精品电影小说 | 国产精品99久久99久久久不卡 | 国产69精品久久久久777片| freevideosex欧美| 国产精品嫩草影院av在线观看| 91精品一卡2卡3卡4卡| 久久精品国产自在天天线| 国产精品久久久久久av不卡| 国产一区有黄有色的免费视频| 日韩av在线免费看完整版不卡| 久久久精品欧美日韩精品| 在线观看一区二区三区激情| 久久人人爽人人爽人人片va| 精品久久久久久久人妻蜜臀av| 精品人妻视频免费看| 精品久久久精品久久久| 麻豆久久精品国产亚洲av| 免费观看的影片在线观看| 国产免费福利视频在线观看| 国产 精品1| 久久久久国产网址| 国产精品秋霞免费鲁丝片| av女优亚洲男人天堂| videos熟女内射| 永久免费av网站大全| 免费观看无遮挡的男女| 久久久欧美国产精品| 97在线视频观看| 成人一区二区视频在线观看| 色吧在线观看| 免费观看的影片在线观看| 日日摸夜夜添夜夜爱| 国产精品久久久久久久久免| 成人免费观看视频高清| 国产精品爽爽va在线观看网站| 又爽又黄无遮挡网站| 内射极品少妇av片p| 尤物成人国产欧美一区二区三区| 日本熟妇午夜| 国产精品久久久久久精品电影| 国产精品嫩草影院av在线观看| 一级毛片aaaaaa免费看小| 精品国产乱码久久久久久小说| 成人鲁丝片一二三区免费| 日本色播在线视频| 亚洲精品一区蜜桃| 午夜激情久久久久久久| 欧美高清成人免费视频www| 青青草视频在线视频观看| 人妻制服诱惑在线中文字幕| 能在线免费看毛片的网站| 日本-黄色视频高清免费观看| videossex国产| 特大巨黑吊av在线直播| 欧美高清成人免费视频www| 精品一区二区免费观看| 美女高潮的动态| 亚洲精品成人久久久久久| 日本三级黄在线观看| 亚洲精品乱久久久久久| 欧美变态另类bdsm刘玥| av在线播放精品| 三级经典国产精品| 久久久色成人| 只有这里有精品99| 亚洲美女搞黄在线观看| 亚洲国产欧美在线一区| 亚洲精华国产精华液的使用体验| 18禁在线无遮挡免费观看视频| 亚洲在线观看片| 在线看a的网站| 国精品久久久久久国模美| 搞女人的毛片| 六月丁香七月| 在线天堂最新版资源| 精品一区二区三区视频在线| 国产精品福利在线免费观看| 内地一区二区视频在线| 成人国产av品久久久| 久久久色成人| 十八禁网站网址无遮挡 | 日本一二三区视频观看| videos熟女内射| 精品人妻视频免费看| 丝袜脚勾引网站| 国产白丝娇喘喷水9色精品| 国产精品久久久久久久电影| 97在线视频观看| 成年女人在线观看亚洲视频 | 久久97久久精品| 91狼人影院| 欧美性猛交╳xxx乱大交人| 国产一区二区三区综合在线观看 | 99热这里只有精品一区| 日韩三级伦理在线观看| 欧美日本视频| 色播亚洲综合网| 亚洲第一区二区三区不卡| 午夜激情福利司机影院| 精品人妻视频免费看| 日日啪夜夜爽| 一级毛片电影观看| www.av在线官网国产| 午夜福利视频1000在线观看| 国产极品天堂在线| 直男gayav资源| 国产av不卡久久| eeuss影院久久| 一级二级三级毛片免费看| 午夜精品一区二区三区免费看| 久久精品国产a三级三级三级| 在线a可以看的网站| 超碰av人人做人人爽久久| 国产午夜精品一二区理论片| 噜噜噜噜噜久久久久久91| 亚洲,欧美,日韩| 亚洲婷婷狠狠爱综合网| 日韩av不卡免费在线播放| 女的被弄到高潮叫床怎么办| 自拍欧美九色日韩亚洲蝌蚪91 | 天堂俺去俺来也www色官网| 搡女人真爽免费视频火全软件| 免费黄频网站在线观看国产| 中国国产av一级| 人人妻人人澡人人爽人人夜夜| 亚洲综合精品二区| freevideosex欧美| 亚洲精品一区蜜桃| 91久久精品国产一区二区三区| 尾随美女入室| 能在线免费看毛片的网站| 91午夜精品亚洲一区二区三区| 亚洲成人一二三区av| 高清日韩中文字幕在线| 午夜免费男女啪啪视频观看| 18禁在线无遮挡免费观看视频| 免费观看在线日韩| 亚洲第一区二区三区不卡| 国产乱来视频区| 亚洲伊人久久精品综合| 国产成人午夜福利电影在线观看| 国产高清有码在线观看视频| 国产成人精品久久久久久| 热re99久久精品国产66热6| 日日撸夜夜添| 亚洲国产精品999| 国语对白做爰xxxⅹ性视频网站| 国产爱豆传媒在线观看| 国产91av在线免费观看| 国产亚洲午夜精品一区二区久久 | 一级毛片我不卡| 2022亚洲国产成人精品| 精品人妻一区二区三区麻豆| 激情五月婷婷亚洲| 黄色日韩在线| 最近的中文字幕免费完整| 国产精品秋霞免费鲁丝片| 精品久久久久久久人妻蜜臀av| 色视频在线一区二区三区| 成年女人看的毛片在线观看| 啦啦啦啦在线视频资源| 免费不卡的大黄色大毛片视频在线观看| 国产伦精品一区二区三区四那| 校园人妻丝袜中文字幕| 三级国产精品欧美在线观看| 女的被弄到高潮叫床怎么办| 国产男女内射视频| 99久久九九国产精品国产免费| 国内精品宾馆在线| 欧美日韩精品成人综合77777| 亚洲天堂国产精品一区在线| 我要看日韩黄色一级片| 亚洲国产高清在线一区二区三| 久久人人爽人人片av| 少妇人妻久久综合中文| 久久亚洲国产成人精品v| 亚洲精品国产av成人精品| 国产午夜精品久久久久久一区二区三区| 久久精品国产鲁丝片午夜精品| 麻豆国产97在线/欧美| 久久鲁丝午夜福利片| 国产人妻一区二区三区在| 成人国产麻豆网| 人人妻人人看人人澡| 性色avwww在线观看| 中文欧美无线码| 国产日韩欧美在线精品| 啦啦啦中文免费视频观看日本| 国产午夜福利久久久久久| 女人被狂操c到高潮| 大码成人一级视频| 国产一区二区三区综合在线观看 | 欧美97在线视频| 一二三四中文在线观看免费高清| 欧美xxⅹ黑人| 中文字幕免费在线视频6| 激情五月婷婷亚洲| 国产有黄有色有爽视频| 嘟嘟电影网在线观看| 一级a做视频免费观看| 波野结衣二区三区在线| 久久99热6这里只有精品| 国产成年人精品一区二区| 亚洲欧美日韩东京热| 国内精品宾馆在线| 久久精品久久久久久噜噜老黄| 成人一区二区视频在线观看| 天天躁日日操中文字幕| av在线观看视频网站免费| 成人美女网站在线观看视频| 国产免费视频播放在线视频| 亚洲精品久久久久久婷婷小说| 又爽又黄a免费视频| 午夜免费鲁丝| 在线观看av片永久免费下载| 亚洲天堂国产精品一区在线| 天美传媒精品一区二区| 国产永久视频网站| 老女人水多毛片| 国产高清有码在线观看视频| 日韩av不卡免费在线播放| 国产综合懂色| 欧美激情在线99| 国产极品天堂在线| 夜夜爽夜夜爽视频| 亚洲婷婷狠狠爱综合网| 欧美日韩在线观看h| 日韩一区二区视频免费看| 一区二区三区乱码不卡18| 亚洲成人中文字幕在线播放| 少妇的逼水好多| 毛片女人毛片| 日本色播在线视频| 午夜免费男女啪啪视频观看| 国内精品宾馆在线| 国产精品秋霞免费鲁丝片| 在线精品无人区一区二区三 | 成人欧美大片| 日韩欧美精品免费久久| 精品久久久久久久久亚洲|