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

    基于矢耳石形態(tài)的平鲉屬(Sebastes)魚(yú)類(lèi)鑒別及研究方法比較

    2015-06-24 14:10:05張鳳俠張秀梅覃樂(lè)政李文濤
    海洋學(xué)報(bào) 2015年4期
    關(guān)鍵詞:許氏判別函數(shù)耳石

    張鳳俠, 張秀梅*, 覃樂(lè)政, 李文濤

    (1.中國(guó)海洋大學(xué) 教育部海水養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室,山東 青島 266003)

    基于矢耳石形態(tài)的平鲉屬(Sebastes)魚(yú)類(lèi)鑒別及研究方法比較

    張鳳俠1, 張秀梅1*, 覃樂(lè)政1, 李文濤1

    (1.中國(guó)海洋大學(xué) 教育部海水養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室,山東 青島 266003)

    分析了青島膠州灣鎧平鲉(Sebasteshubbsi)、花斑平鲉(Sebastesnigricans)和許氏平鲉(Sebastesschlegelii) 3種平鲉屬魚(yú)類(lèi)的矢耳石形態(tài)特征,探究基于矢耳石形態(tài)鑒別3種魚(yú)類(lèi)的可行性。以8個(gè)耳石形狀指標(biāo)和37個(gè)橢圓傅里葉系數(shù)進(jìn)行的單因素方差分析結(jié)果基本一致,顯示鎧平鲉與許氏平鲉的耳石形態(tài)差異最大,鎧平鲉與花斑平鲉間差異最小。利用耳石形狀指標(biāo)、橢圓傅里葉系數(shù)、耳石形狀指標(biāo)與橢圓傅里葉系數(shù)相結(jié)合的方法,對(duì)3種魚(yú)類(lèi)的綜合判別率均可達(dá)95%以上,判別分析散點(diǎn)圖可明顯將3種魚(yú)類(lèi)區(qū)分開(kāi)。耳石長(zhǎng)介于4~7 mm的3種魚(yú)類(lèi)小規(guī)格組與2齡及以上成魚(yú)組耳石樣本的判別結(jié)果一致,表明針對(duì)不同發(fā)育階段平鲉屬魚(yú)類(lèi),基于耳石形態(tài)的種類(lèi)鑒別均是可行的;僅使用橢圓傅里葉系數(shù)鑒別3種魚(yú)類(lèi)時(shí)綜合判別率雖高達(dá)100%,但難以快速定位3種魚(yú)類(lèi)耳石形態(tài)間,在統(tǒng)計(jì)上差異顯著而實(shí)際上微小的形狀差異所在;而使用耳石形狀指標(biāo)和橢圓傅里葉系數(shù)聯(lián)合判別,綜合判別率亦高達(dá)99%,且能快速找出3種魚(yú)類(lèi)耳石在矩形趨近率、形態(tài)因子、面密度等的形態(tài)差異。研究表明,耳石形狀指標(biāo)和橢圓傅里葉分析聯(lián)合判別是平鲉屬魚(yú)類(lèi)種類(lèi)鑒別及耳石形態(tài)差異分析的有效手段。

    平鲉屬;耳石形態(tài);耳石形狀指標(biāo);橢圓傅里葉系數(shù);鑒別

    1 引言

    耳石是魚(yú)類(lèi)生長(zhǎng)過(guò)程中沉淀在內(nèi)耳中的鈣質(zhì)結(jié)構(gòu),分矢耳石(sagittae)、微耳石(lapillus)和星耳石(asteriscus)共3對(duì)。耳石隨著魚(yú)類(lèi)的生長(zhǎng)而生長(zhǎng)[1],一般認(rèn)為達(dá)到性成熟后耳石形態(tài)趨于穩(wěn)定[2-3]。由于耳石具有不易消化的特點(diǎn),被廣泛應(yīng)用于食魚(yú)性動(dòng)物的食性分析中[4]。研究表明,耳石形態(tài)主要受遺傳[5]、環(huán)境[6]、攝食[7]等因素影響,其中遺傳占主導(dǎo)因素[8-9]。Morales-Nin[10]認(rèn)為魚(yú)類(lèi)耳石在不同空間方向上的沉積導(dǎo)致耳石形狀不規(guī)則;Maisey[11]、葉振江[12]指出耳石形態(tài)具有高度的物種特異性和顯著的種間差異性;Begg等[13]指出采用先進(jìn)的數(shù)碼圖像分析耳石形狀為種群識(shí)別提供可能。眾多研究表明,耳石二維圖像的數(shù)學(xué)分析可以定量解釋魚(yú)類(lèi)種間、群體間的耳石形態(tài)差異,可用于魚(yú)類(lèi)種間、群體間的識(shí)別研究。

    依據(jù)耳石形態(tài),在鮨屬(Serranus)[14]、鯉科(Cyprinidae)[15]、鱭屬(Coilia)[16]、鱒屬(Salmo)[17]、四大家魚(yú)仔稚魚(yú)[18]以及熱帶鳀科(Engraulididae)[19]等魚(yú)類(lèi)的種類(lèi)鑒別中取得了較好的鑒別效果。但在耳石形態(tài)分析時(shí),參數(shù)選擇和數(shù)據(jù)分析方法不同,會(huì)顯著影響鑒別結(jié)果。目前,常用的耳石分析方法包括文字描述法[2]、可量性狀比較法[16]、地標(biāo)法[20]、耳石形狀指標(biāo)法[9]與橢圓傅里葉分析法[21]。耳石形狀指標(biāo)法運(yùn)用于九帶鮨(Serranuscabrilla)[22]、鱚屬(Sillago)[20]等的鑒別研究中;耳石橢圓傅里葉分析已成功進(jìn)行竹筴魚(yú)(Trachurustrachurus)[23]和條紋婢(Latrislineata)[24]的群體鑒別;耳石形狀指標(biāo)法結(jié)合橢圓傅里葉分析,已在黑線(xiàn)鱈(Melanogrammusaeglefinus)[25]和大西洋鱈(Gadusmorhua)[26]的鑒別中應(yīng)用。

    平鲉屬(Sebastes)隸屬于鲉形目(Scorpaeniformes)鲉亞目(Scorpaenoidei)的鲉科(Scorpaenidae),我國(guó)常見(jiàn)分布種為鎧平鲉(Sebasteshubbsi)、花斑平鲉(Sebastesnigricans)和許氏平鲉(Sebastesschlegelii),主要分布在黃海、渤海海域,為巖礁區(qū)底層經(jīng)濟(jì)魚(yú)類(lèi)。目前3種魚(yú)類(lèi)的鑒別主要依據(jù)背鰭鰭棘數(shù)、鰭條數(shù)和體側(cè)斑點(diǎn)等生物學(xué)特征。依據(jù)外部形態(tài)特征的種類(lèi)鑒別,當(dāng)樣品保存不當(dāng)時(shí),種類(lèi)鑒別則受到限制,此時(shí)形態(tài)穩(wěn)定性較好的耳石則可作為分類(lèi)的輔助手段。截至目前,基于耳石形態(tài)的平鲉屬魚(yú)類(lèi)的種類(lèi)鑒別尚無(wú)報(bào)道,本研究以青島膠州灣的平鲉屬魚(yú)類(lèi)為研究對(duì)象,旨在探究基于耳石形態(tài)鑒別平鲉屬魚(yú)類(lèi)的可行性,為傳統(tǒng)形態(tài)分類(lèi)提供佐證和輔助資料;并比較分析不同方法對(duì)鑒別結(jié)果的影響,嘗試找出適合平鲉屬魚(yú)類(lèi)耳石分類(lèi)的研究方法,為今后開(kāi)展相關(guān)耳石研究提供基礎(chǔ)資料。

    2 材料和方法

    2.1 實(shí)驗(yàn)材料

    實(shí)驗(yàn)所用平鲉屬樣本于2013年5-8月采自青島膠州灣海域地籠網(wǎng)和拖網(wǎng)漁獲物。種類(lèi)鑒定后測(cè)量體長(zhǎng)(精確到1 mm)、體質(zhì)量(精確到0.01 g)(表1)。取出左、右矢耳石(以下皆以耳石簡(jiǎn)稱(chēng)),洗凈后保存于離心管中待進(jìn)一步分析。

    表1 膠州灣3種平鲉屬魚(yú)類(lèi)樣品的體長(zhǎng)與體質(zhì)量

    2.2 研究方法

    2.2.1 年齡鑒定及耳石樣本的選取

    將洗凈的耳石浸入盛水的培養(yǎng)皿中,凹面向上,置于Nikon SMZ8000型體式顯微鏡下,透射光觀(guān)察,可見(jiàn)清晰的輪紋,其中一個(gè)完整的明帶和暗帶計(jì)為一個(gè)年輪。左、右耳石的年齡鑒定結(jié)果一致,實(shí)驗(yàn)首先選取2齡及以上平鲉屬魚(yú)類(lèi)耳石用于后續(xù)分析。

    2.2.2 耳石形態(tài)參數(shù)的選取

    將選取的耳石置于50℃烘箱中干燥12 h以上至恒重。干燥后的耳石用精密電子天平稱(chēng)重,單位精確到0.01 mg。耳石凸面向上、凹面向下,用連接于體式顯微鏡的Nikon DN100數(shù)碼相機(jī)拍攝耳石的二維圖像(見(jiàn)圖1)。

    采用Image-ProPlus6.0圖像分析軟件分析耳石圖像,測(cè)量耳石形態(tài)學(xué)參數(shù),共得出8個(gè)尺寸參數(shù):耳石面積、耳石周長(zhǎng)、耳石長(zhǎng)、耳石寬、最大Feret直徑、最小Feret直徑、最大半徑和最小半徑。其中半徑測(cè)量的是通過(guò)耳石幾何中心的半徑長(zhǎng),F(xiàn)eret直徑指耳石按照某一角度擺放時(shí),經(jīng)過(guò)形心的耳石外緣兩點(diǎn)間的距離,面積為耳石的二維投影面積。根據(jù)潘曉哲和高天翔[20]、徐勝勇等[27]、Tuset等[28]、Canas等[29]相關(guān)研究,結(jié)合平鲉屬魚(yú)類(lèi)的耳石形態(tài)特征,將8個(gè)耳石尺寸參數(shù)和耳石質(zhì)量通過(guò)計(jì)算公式的轉(zhuǎn)換,得出平鲉屬魚(yú)類(lèi)的8個(gè)耳石形狀指標(biāo)(Otolith shape indices)(見(jiàn)表2)。耳石形狀指標(biāo)在一定程度上可消除因耳石放置位置不同、耳石大小不同等因素對(duì)圖像分析結(jié)果的影響。

    圖1 膠州灣3種平鲉屬魚(yú)類(lèi)左耳石內(nèi)側(cè)面照片(BL代表體長(zhǎng),OL代表耳石長(zhǎng))Fig.1 Left sagittal otoliths of three Sebastes collected in the Jiaozhou Bay shown with the proximal surface up (BL represents body length, OL represents otolith length)

    表2 平鲉屬魚(yú)類(lèi)耳石尺寸參數(shù)和形狀指標(biāo)

    Tab.2 Size parameters and size based shape indices with calculation formulas forSebastesspecies

    尺寸參數(shù)耳石形狀指標(biāo)面積(A)圓度=(4A)/(πOL2)周長(zhǎng)(P)矩形趨近率=A/(OL×OW)耳石長(zhǎng)(OL)橢圓率=(OL-OW)/(OL+OW)耳石寬(OW)形態(tài)因子=(4πA)/P2耳石質(zhì)量(OW′)幅形比=OL/OW最大Feret直徑(Fmax)Feret直徑比=Fmax/Fmin最小Feret直徑(Fmin)半徑比=Rmax/Rmin最大半徑(Rmax)面密度=OW′/A最小半徑(Rmin)

    注:文中以O(shè)W′代表耳石質(zhì)量。面積度量單位為mm2,長(zhǎng)度度量單位為mm,質(zhì)量度量單位為mg。

    橢圓傅里葉分析(Elliptic Fourier analysis)能夠有效對(duì)耳石輪廓進(jìn)行描述,定量獲取耳石外輪廓的有效信息[30]。研究表明,傅里葉諧值越多越能精確地描述物體的外輪廓[31],但前10~20個(gè)傅里葉特征系數(shù)能夠描述耳石的主要形態(tài)輪廓[32]。本研究使用SHAPE1.3軟件包前10組傅里葉諧波對(duì)耳石圖像進(jìn)行分析。每組諧波由4個(gè)形態(tài)變量an、bn、cn、dn組成,因此,每個(gè)耳石共由40個(gè)傅里葉特征系數(shù)描述。對(duì)10個(gè)諧值進(jìn)行標(biāo)準(zhǔn)化后,前3個(gè)特征系數(shù)為定值:a1=1,b1=c1=0,故有效傅里葉特征系數(shù)為37個(gè)。

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

    本研究分別對(duì)3種平鲉屬魚(yú)類(lèi)左、右耳石各原始測(cè)量指標(biāo)進(jìn)行獨(dú)立樣本T檢驗(yàn),以檢驗(yàn)3種魚(yú)類(lèi)種內(nèi)左、右耳石間的差異性。同樣使用獨(dú)立樣本T檢驗(yàn),檢測(cè)3種魚(yú)類(lèi)雌、雄個(gè)體間耳石形態(tài)的差異情況。顯著性水平取α=0.05。

    鑒于魚(yú)類(lèi)個(gè)體生長(zhǎng)差異可能對(duì)檢驗(yàn)結(jié)果產(chǎn)生影響,為實(shí)現(xiàn)無(wú)偏分析,在樣本選擇時(shí)除選取2齡及以上的性成熟個(gè)體外,還在數(shù)據(jù)分析時(shí)引入?yún)f(xié)變量(體長(zhǎng))以降低或消除個(gè)體生長(zhǎng)差異對(duì)耳石形態(tài)參數(shù)的影響。本文采用Pearson相關(guān)性分析確定耳石形態(tài)參數(shù)(8個(gè)耳石形狀指標(biāo)和37個(gè)橢圓傅里葉系數(shù))與體長(zhǎng)間關(guān)系。當(dāng)體長(zhǎng)與耳石形態(tài)參數(shù)存在線(xiàn)性相關(guān)時(shí),參照竇碩增[33]的協(xié)方差校正方法,利用回歸斜率按照體長(zhǎng)校正公式校正各耳石形態(tài)參數(shù):

    (1)

    以單因素方差分析檢驗(yàn)各個(gè)耳石形態(tài)參數(shù)在3種魚(yú)類(lèi)間的差異性,用多重比較驗(yàn)證差異來(lái)源。以上統(tǒng)計(jì)分析過(guò)程中,對(duì)具有方差齊性的變量采用LSD法,不具方差齊性的采用Tamhane’s T2法,取顯著性水平α=0.05。

    利用SPSS18.0對(duì)3種平鲉屬魚(yú)類(lèi)作判別分析,本文采用3種方式進(jìn)行對(duì)比分析:(1)僅使用8個(gè)耳石形狀指標(biāo);(2)僅使用37個(gè)橢圓傅里葉特征系數(shù);(3)8個(gè)耳石形狀指標(biāo)并結(jié)合37個(gè)橢圓傅里葉特征系數(shù)。

    2.2.4 對(duì)比分析樣本的補(bǔ)充

    2齡及以上3種平鲉屬魚(yú)類(lèi)間體長(zhǎng)、耳石長(zhǎng)存在較大差異,而近岸捕撈生產(chǎn)中常見(jiàn)低于1齡的許氏平鲉幼魚(yú)。因此,為驗(yàn)證鑒別方法的有效性,本實(shí)驗(yàn)對(duì)比加入不足1齡的許氏平鲉幼魚(yú)樣本,即補(bǔ)充加入耳石長(zhǎng)介于4~7 mm的小個(gè)體許氏平鲉樣本。將耳石長(zhǎng)介于4~7 mm的3種平鲉屬魚(yú)類(lèi)單獨(dú)進(jìn)行判別分析,以充分探討耳石形態(tài)在3種魚(yú)類(lèi)種類(lèi)鑒別中的可行性。在此,采用耳石長(zhǎng)作為協(xié)變量引入?yún)f(xié)方差校正方程,校正與耳石長(zhǎng)顯著相關(guān)的參數(shù),校正后的參數(shù)用于后續(xù)判別分析,進(jìn)而比較交叉檢驗(yàn)結(jié)果。

    為有效區(qū)分,下文中2齡及以上許氏平鲉樣本統(tǒng)一稱(chēng)為許氏平鲉;耳石長(zhǎng)介于4~7 mm的小個(gè)體許氏平鲉統(tǒng)一稱(chēng)為許氏平鲉(小)。下文中,若未詳細(xì)說(shuō)明,則3種平鲉屬魚(yú)類(lèi)均指2齡及以上、耳石形態(tài)趨于穩(wěn)定的樣本;耳石長(zhǎng)介于4~7 mm的3種魚(yú)類(lèi)樣本信息為:鎧平鲉(n=61,體長(zhǎng)范圍73~119 mm,耳石長(zhǎng)范圍4.45~6.58 mm)、花斑平鲉(n=14,體長(zhǎng)范圍102~165 mm,耳石長(zhǎng)范圍4.97~6.64 mm)、許氏平鲉(小)(n=45,體長(zhǎng)范圍101~131 mm,耳石長(zhǎng)范圍4.29~5.93 mm)。

    3 結(jié)果

    3.1 種內(nèi)耳石形態(tài)的差異性

    對(duì)左、右耳石的8個(gè)測(cè)量指標(biāo)和耳石質(zhì)量分別進(jìn)行T檢驗(yàn),結(jié)果顯示3種平鲉屬魚(yú)類(lèi)的左、右耳石在這些測(cè)量指標(biāo)上均無(wú)顯著性差異。

    對(duì)雌、雄個(gè)體耳石的8個(gè)測(cè)量指標(biāo)和耳石質(zhì)量分別進(jìn)行T檢驗(yàn),結(jié)果顯示3種平鲉屬魚(yú)類(lèi)雌、雄個(gè)體的耳石在這些測(cè)量指標(biāo)上均無(wú)顯著性差異。

    因此,下文中的耳石形態(tài)研究不分性別統(tǒng)一使用左耳石。表3為以鎧平鲉為例,左、右耳石以及雌、雄個(gè)體間耳石差異性的檢驗(yàn)結(jié)果。

    表3 鎧平鲉左右耳石、雌雄個(gè)體間耳石差異性比較

    注:n♀表示鎧平鲉雌性個(gè)體的數(shù)量,n♂表示雄性個(gè)體的數(shù)量。

    3.2 單因素方差分析

    Pearson相關(guān)分析得出3個(gè)耳石形狀指標(biāo)(橢圓率、形態(tài)因子、面密度)及14個(gè)橢圓傅里葉系數(shù)(h4、h6、h11、h12、h14、h17、h18、h19、h22、h23、h25、h29、h30、h37)與體長(zhǎng)顯著相關(guān),這些變量使用協(xié)方差校正方法校正后進(jìn)行后續(xù)的數(shù)據(jù)分析。

    8個(gè)耳石形狀指標(biāo)及單因素方差分析結(jié)果如表4所示。鎧平鲉和許氏平鲉,花斑平鲉和許氏平鲉間均有7個(gè)耳石形狀指標(biāo)存在顯著差異;鎧平鲉和花斑平鲉有3個(gè)指標(biāo)存在顯著差異。圓度在3種平鲉屬魚(yú)類(lèi)間差異均不顯著,3種魚(yú)類(lèi)兩兩之間在矩形趨近率、形態(tài)因子、面密度間均差異顯著。

    表4 平鲉屬魚(yú)類(lèi)8個(gè)耳石形狀指標(biāo)均值及單因素方差分析結(jié)果

    注:表中耳石形狀信息為原始數(shù)據(jù),顯著性判定為依據(jù)校正后的數(shù)據(jù)。上述字母不同表示差異顯著,字母相同表示差異不顯著(p<0.05)。

    對(duì)3種魚(yú)類(lèi)的37個(gè)傅里葉特征系數(shù)進(jìn)行單因素方差分析,將兩兩比較具有差異性顯著的特征系數(shù)的個(gè)數(shù)及其所占總系數(shù)(37)的比例列入表5??梢钥闯觯?種魚(yú)類(lèi)在多個(gè)特征系數(shù)上的差異具有明顯的統(tǒng)計(jì)意義。其中,鎧平鲉和許氏平鲉間的差異最大,有16個(gè)特征系數(shù)存在顯著性差異,占總系數(shù)的43.2%;鎧平鲉和花斑平鲉之間的差異較小,有12個(gè)特征系數(shù)存在差異,占總系數(shù)的32.4%。

    表5 平鲉屬魚(yú)類(lèi)耳石傅里葉特征系數(shù)顯著個(gè)數(shù)及比例

    3.3 種類(lèi)識(shí)別

    分別采用3種方式進(jìn)行平鲉屬魚(yú)類(lèi)種類(lèi)識(shí)別的判別分析結(jié)果如表6所示。交叉驗(yàn)證結(jié)果表明:

    僅耳石形狀指標(biāo),典則判別函數(shù)一能夠解釋91.9%的魚(yú)種間差異,可以有效地將許氏平鲉與鎧平鲉和花斑平鲉區(qū)分開(kāi);典則判別函數(shù)二能夠解釋8.1%的魚(yú)種間差異,可以將鎧平鲉與花斑平鲉區(qū)分開(kāi)(見(jiàn)圖2a)。對(duì)判別函數(shù)一貢獻(xiàn)最大的參數(shù)為圓度,對(duì)判別函數(shù)二貢獻(xiàn)最大的參數(shù)為形態(tài)因子。僅通過(guò)耳石形狀指標(biāo)進(jìn)行的種類(lèi)劃分,鎧平鲉與花斑平鲉之間有少許的誤判,而許氏平鲉與鎧平鲉、花斑平鲉之間無(wú)誤判出現(xiàn)。表6可以看出,3種魚(yú)類(lèi)判別正確率為93.3%~100%,綜合判別正確率為96.2%。

    僅橢圓傅里葉系數(shù),典則判別函數(shù)一能夠解釋60.2%的魚(yú)種間差異,可以有效地將鎧平鲉和許氏平鲉區(qū)分開(kāi);典則判別函數(shù)二能夠解釋39.8%的魚(yú)種間差異,可以有效地將花斑平鲉與鎧平鲉和許氏平鲉區(qū)分開(kāi)(見(jiàn)圖2b)。對(duì)判別函數(shù)一貢獻(xiàn)最大的參數(shù)為第12號(hào)傅里葉系數(shù),對(duì)判別函數(shù)二貢獻(xiàn)最大的參數(shù)為第17號(hào)傅里葉系數(shù)。通過(guò)橢圓傅里葉系數(shù)進(jìn)行的判別分析,3種魚(yú)類(lèi)無(wú)任何誤判,各魚(yú)種判別正確率均為100%(見(jiàn)表6)。

    結(jié)合耳石形狀指標(biāo)和橢圓傅里葉系數(shù),典則判別函數(shù)一能夠解釋97.3%的魚(yú)種間差異,可以非常有效地將許氏平鲉與鎧平鲉和花斑平鲉區(qū)分開(kāi);典則判別函數(shù)二能夠解釋2.7%的魚(yú)種間差異,可以有效地將花斑平鲉與鎧平鲉和許氏平鲉區(qū)分開(kāi)(見(jiàn)圖2c)。對(duì)判別函數(shù)一貢獻(xiàn)最大的參數(shù)為圓度,對(duì)判別函數(shù)二貢獻(xiàn)最大的參數(shù)為Feret比。結(jié)合耳石形狀指標(biāo)和橢圓傅里葉系數(shù)的判別分析,3種魚(yú)類(lèi)間無(wú)錯(cuò)判,判別正確率均為100%(見(jiàn)表6)。

    綜上,僅傅里葉分析以及耳石形狀指標(biāo)結(jié)合傅里葉系數(shù)的判別結(jié)果,綜合判別正確率均高達(dá)100%;僅采用耳石形狀指標(biāo)進(jìn)行的判別分析,綜合判別率雖未達(dá)100%,但亦達(dá)較高正確率。

    表6 基于平鲉屬3種魚(yú)類(lèi)耳石形態(tài)的交叉檢驗(yàn)結(jié)果

    圖2 基于平鲉屬耳石形態(tài)的前2個(gè)判別函數(shù)值的分析結(jié)果散點(diǎn)圖Fig.2 Scatter plot of scores based on the first two DFA of Sebastes sagittae

    3.4 耳石長(zhǎng)介于4~7 mm樣本的對(duì)比鑒別結(jié)果

    對(duì)耳石長(zhǎng)介于4~7 mm的3種魚(yú)類(lèi)進(jìn)行Pearson相關(guān)分析,校正后的耳石參數(shù)進(jìn)行后續(xù)的判別分析。亦采用3種判別方式進(jìn)行3種魚(yú)類(lèi)的交叉檢驗(yàn),結(jié)果如下:

    僅耳石形狀指標(biāo),典則判別函數(shù)一能夠解釋65.3%的魚(yú)種間差異,可以有效地將鎧平鲉與許氏平鲉(小)區(qū)分開(kāi);典則判別函數(shù)二能夠解釋34.7%的魚(yú)種間差異,可以將花斑平鲉與鎧平鲉和許氏平鲉(小)區(qū)分開(kāi)(見(jiàn)圖3a)。對(duì)判別函數(shù)一貢獻(xiàn)最大的參數(shù)為圓度,對(duì)判別函數(shù)二貢獻(xiàn)最大的參數(shù)為橢圓率。由表7可以看出,僅通過(guò)耳石形狀指標(biāo)進(jìn)行的種類(lèi)劃分,3種魚(yú)類(lèi)間均出現(xiàn)了誤判。61尾鎧平鲉,有1尾誤判為花斑平鲉,2尾誤判為許氏平鲉(小);14尾花斑平鲉,1尾誤判為許氏平鲉(小);45尾許氏平鲉(小),1尾誤判為花斑平鲉;可見(jiàn)花斑平鲉與許氏平鲉(小)耳石較鎧平鲉更相似。3種魚(yú)類(lèi)的判別正確率為92.9%~97.8%,綜合判別正確率為95.3%。

    僅橢圓傅里葉系數(shù),典則判別函數(shù)一能夠解釋59.5%的魚(yú)種間差異,可以有效地將許氏平鲉(小)與鎧平鲉和花斑平鲉區(qū)分開(kāi);典則判別函數(shù)二能夠解釋40.5%的魚(yú)種間差異,可以有效地將鎧平鲉與花斑平鲉區(qū)分開(kāi)(見(jiàn)圖3b)。對(duì)判別函數(shù)一貢獻(xiàn)最大的參數(shù)為第4號(hào)傅里葉系數(shù),對(duì)判別函數(shù)二貢獻(xiàn)最大的參數(shù)為第15號(hào)傅里葉系數(shù)。通過(guò)橢圓傅里葉系數(shù)進(jìn)行的判別分析,3種魚(yú)類(lèi)無(wú)任何誤判,各魚(yú)種判別正確率均為100%(表7)。

    結(jié)合耳石形狀指標(biāo)和橢圓傅里葉系數(shù),典則判別函數(shù)一能夠解釋60.8%的魚(yú)種間差異,可以非常有效地將鎧平鲉與花斑平鲉和許氏平鲉(小)區(qū)分開(kāi);典則判別函數(shù)二能夠解釋39.2%的魚(yú)種間差異,可以有效地將花斑平鲉與許氏平鲉(小)區(qū)分開(kāi)(見(jiàn)圖3c)。對(duì)判別函數(shù)一貢獻(xiàn)最大的參數(shù)為圓度,對(duì)判別函數(shù)二貢獻(xiàn)最大的參數(shù)為橢圓率。結(jié)合耳石形狀指標(biāo)和橢圓傅里葉系數(shù)的判別分析:60尾鎧平鲉有2尾錯(cuò)判為許氏平鲉(小),鎧平鲉的判別正確率為96.7%;花斑平鲉和許氏平鲉(小)的判別正確率均為100%(表7)。

    綜上,僅使用橢圓傅里葉系數(shù)的判別分析,綜合判別正確率高達(dá)100%;僅采用耳石形狀指標(biāo)以及耳石形狀指標(biāo)和傅里葉系數(shù)聯(lián)合判別分析,綜合判別率雖未達(dá)100%,但亦達(dá)較高正確率。

    表7 耳石長(zhǎng)4~7 mm的3種魚(yú)類(lèi)耳石形態(tài)的交叉檢驗(yàn)結(jié)果

    注:文中以S.schlegelii′表示許氏平鲉(小)的拉丁名。

    圖3 耳石長(zhǎng)4~7 mm的3種魚(yú)類(lèi)判別分析結(jié)果散點(diǎn)圖Fig.3 Scatter plot of scores based on the first two DFA of the three Sebastes whose sagittae length were 4 to 7 mm

    4 討論

    4.1 樣本選擇及耳石形態(tài)參數(shù)的選擇與校正

    研究發(fā)現(xiàn),鎧平鲉在2齡時(shí)初次達(dá)到性成熟[34]、許氏平鲉在2~3齡時(shí)達(dá)性成熟[35]。本文首先選取2齡及以上樣本探究基于耳石形態(tài)鑒別平鲉屬魚(yú)類(lèi)的可行性,原則上排除了幼魚(yú)耳石形狀不穩(wěn)定對(duì)實(shí)驗(yàn)結(jié)果的影響,且3種魚(yú)類(lèi)不同體長(zhǎng)間耳石形狀指標(biāo)的單因素方差分析結(jié)果顯示,2齡及以上3種魚(yú)類(lèi)的耳石形態(tài)已趨于穩(wěn)定;另外,根據(jù)實(shí)驗(yàn)采集的樣本和文獻(xiàn)報(bào)道[27]發(fā)現(xiàn),近岸捕撈的許氏平鲉多為體長(zhǎng)80~180 mm不足1齡的幼魚(yú),此體長(zhǎng)階段的耳石長(zhǎng)恰為4~7 mm,因此文中對(duì)比加入耳石長(zhǎng)介于4~7 mm的小個(gè)體許氏平鲉樣本,將耳石長(zhǎng)介于4~7 mm的3種平鲉屬魚(yú)類(lèi)單獨(dú)進(jìn)行基于耳石形態(tài)的判別分析,充分探討了針對(duì)不同生長(zhǎng)階段3種魚(yú)類(lèi)種類(lèi)鑒別的可行性及分析方法的有效性。

    同時(shí),研究樣本均采自青島膠州灣海域,一定程度上克服了地理環(huán)境不同對(duì)耳石形態(tài)的影響,降低了環(huán)境因素對(duì)實(shí)驗(yàn)結(jié)果的干擾。由于花斑平鲉在我國(guó)僅分布在黃海海域且資源量較少,故采集到的樣本數(shù)量有限,今后應(yīng)進(jìn)一步增加花斑平鲉耳石樣本數(shù)量,以提高實(shí)驗(yàn)結(jié)果的可靠性。

    耳石形態(tài)特征主要受遺傳因素的影響,魚(yú)類(lèi)的生物學(xué)特征(性別、耳石位置、個(gè)體大小等)亦對(duì)耳石的形態(tài)特征產(chǎn)生一定影響,秦巖等[2]報(bào)道了鲬?hù)~(yú)(Platycephalusindicus)雌雄個(gè)體間矢耳石形態(tài)特征存在差異,陸化杰等[36]指出不同個(gè)體、不同性別間阿根廷滑柔魚(yú)(Illexargentinus)耳石形狀參數(shù)存在顯著差異。本文逐一檢驗(yàn)了3種平鲉屬魚(yú)類(lèi)的左、右耳石以及雌、雄個(gè)體間耳石的形態(tài)差異性。相關(guān)研究表明,當(dāng)影響因素如性別、耳石位置等對(duì)耳石形態(tài)影響水平不顯著時(shí),若樣本量充足則可去除該因素,若樣本量不足則應(yīng)匯總樣本進(jìn)行分析[37]。鑒于本研究樣本數(shù)較少,且3種魚(yú)類(lèi)耳石在雌雄個(gè)體間均無(wú)顯著差異,故將雌雄個(gè)體的左耳石一起作為實(shí)驗(yàn)樣本納入后續(xù)分析中。

    對(duì)于耳石參數(shù)的選擇,包含了稱(chēng)量指標(biāo)耳石質(zhì)量、由耳石大小變量進(jìn)行轉(zhuǎn)換的耳石形狀指標(biāo)和耳石輪廓變量傅里葉系數(shù),引入耳石質(zhì)量信息增加了耳石的三維特征,結(jié)果表明耳石質(zhì)量信息的加入可以提高判別率[20]。以耳石形態(tài)鑒別3種魚(yú)類(lèi),為實(shí)現(xiàn)無(wú)偏分析,鑒于2齡及以上的3種平鲉屬魚(yú)類(lèi)體長(zhǎng)差異較大,故對(duì)其進(jìn)行體長(zhǎng)效應(yīng)的校正;而耳石長(zhǎng)介于4~7 mm的3種魚(yú)類(lèi),體長(zhǎng)差異較小,故選擇通過(guò)耳石長(zhǎng)進(jìn)行校正。校正后的相關(guān)指標(biāo)用于后續(xù)的統(tǒng)計(jì)分析中,排除了魚(yú)體大小、耳石大小對(duì)判別結(jié)果的影響。

    4.2 多元統(tǒng)計(jì)結(jié)果的分析

    8個(gè)耳石形狀指標(biāo)的單因子方差分析結(jié)果(見(jiàn)表4)顯示,3種平鲉屬魚(yú)類(lèi)間耳石形狀指標(biāo)存在顯著差異,尤其是鎧平鲉與許氏平鲉、花斑平鲉與許氏平鲉的耳石,有7個(gè)形狀指標(biāo)存在顯著差異。對(duì)3種魚(yú)類(lèi)的37個(gè)橢圓傅里葉特征系數(shù)進(jìn)行兩兩比較,傅里葉系數(shù)差異顯著的個(gè)數(shù)占總特征系數(shù)的比例均達(dá)30%以上,其中鎧平鲉與許氏平鲉之間的差異高達(dá)43.2%(見(jiàn)表5)。綜合可見(jiàn),2種方法的方差分析結(jié)果基本一致,鎧平鲉與許氏平鲉的耳石形態(tài)差異最大,鎧平鲉與花斑平鲉則差異最小。其中,差異不顯著的參數(shù)為平鲉屬魚(yú)類(lèi)耳石的共有特征,如3種魚(yú)類(lèi)間耳石圓度差異不顯著,表明該圓度范圍即為3種魚(yú)類(lèi)耳石的共有特征;而對(duì)比分析中差異顯著的參數(shù)為耳石的種間差異所在,即為基于耳石形態(tài)鑒別平鲉屬魚(yú)類(lèi)提供了可能。

    本研究采用耳石形狀指標(biāo)、橢圓傅里葉系數(shù)以及耳石形狀指標(biāo)與橢圓傅里葉系數(shù)聯(lián)合判別法進(jìn)行3種平鲉屬魚(yú)類(lèi)的鑒別,結(jié)果顯示:(1)3種判別方式對(duì)于3種魚(yú)類(lèi)(包括大小規(guī)格的許氏平鲉)的綜合判別率均高達(dá)95%以上,且判別函數(shù)所作散點(diǎn)圖可明顯將3種魚(yú)類(lèi)區(qū)分開(kāi)來(lái)(見(jiàn)圖2,圖3),可見(jiàn)基于耳石形態(tài)進(jìn)行平鲉屬魚(yú)類(lèi)的鑒別是可行的,3種判別方式均可作為平鲉屬魚(yú)類(lèi)種類(lèi)鑒別的輔助方法。(2)僅使用耳石形狀指標(biāo)進(jìn)行的判別分析結(jié)果顯示,大個(gè)體許氏平鲉與鎧平鲉和花斑平鲉間無(wú)誤判出現(xiàn),表明其與鎧平鲉和花斑平鲉間的耳石形態(tài)差異較大;而小個(gè)體的許氏平鲉與鎧平鲉、花斑平鲉均有一定程度的誤判(最大誤判率為3.3%),說(shuō)明低于1齡的許氏平鲉(小個(gè)體)與鎧平鲉、花斑平鲉的耳石形態(tài)差異較之大個(gè)體許氏平鲉相對(duì)較小。這表明魚(yú)體生長(zhǎng)過(guò)程中許氏平鲉耳石發(fā)生一定的形態(tài)變化,亦驗(yàn)證了隨著魚(yú)體的生長(zhǎng),耳石亦經(jīng)歷生長(zhǎng)、發(fā)育和穩(wěn)定過(guò)程的論點(diǎn)。而圖2、圖3的散點(diǎn)圖顯示,無(wú)論許氏平鲉處于哪個(gè)生長(zhǎng)階段,依據(jù)耳石均能將3種平鲉屬魚(yú)類(lèi)區(qū)分開(kāi)來(lái),表明基于耳石形態(tài)進(jìn)行不同生長(zhǎng)階段的平鲉屬魚(yú)類(lèi)種類(lèi)鑒別是可行性的。(3)比較3種判別方式對(duì)平鲉屬3種魚(yú)類(lèi)的判別結(jié)果,由表6、表7(包含大小規(guī)格的許氏平鲉)可知,僅使用橢圓傅里葉系數(shù)的判別分析均能實(shí)現(xiàn)100%的判別正確率,相較于其他2種判別方法,該研究方法在探討基于耳石形態(tài)的3種平鲉屬魚(yú)類(lèi)的鑒別中更具優(yōu)勢(shì)。如李輝華等[38]曾證實(shí)對(duì)于長(zhǎng)江鳳鱭(Coilianasus)和刀鱭(Coiliamystus)的種間和種內(nèi)鑒別,傅里葉分析法較傳統(tǒng)形態(tài)測(cè)量法更具優(yōu)勢(shì)。然而,傅里葉系數(shù)雖可以準(zhǔn)確的定量表達(dá)耳石外輪廓,但是與魚(yú)類(lèi)種間的實(shí)際耳石形狀差異較難聯(lián)系起來(lái),難以快速找出不同種間的耳石在統(tǒng)計(jì)上差異顯著而實(shí)際上微小的形狀差異所在[33,39]。本研究顯示,結(jié)合耳石形狀指標(biāo)法和橢圓傅里葉分析2種方法,對(duì)3種平鲉屬魚(yú)類(lèi)的綜合判別率亦高達(dá)99%,且由耳石形狀指標(biāo)所得的單因素方差分析結(jié)果可知(見(jiàn)表4),3種魚(yú)類(lèi)耳石在矩形趨近率、形態(tài)因子、面密度等耳石形狀指標(biāo)上存在顯著差異,因此,耳石形狀指標(biāo)法和橢圓傅里葉分析的聯(lián)合使用是基于耳石形態(tài)進(jìn)行平鲉屬魚(yú)類(lèi)種類(lèi)鑒別的有效手段,該方法既能快速找出3種魚(yú)類(lèi)間耳石形態(tài)的差異所在,又能實(shí)現(xiàn)較高的種間判別正確率。研究表明,耳石形狀指標(biāo)結(jié)合橢圓傅里葉系數(shù)對(duì)鮟鱇(Lophiuspiscatorius)[29]、鲯鰍(Coryphaenahippurus)[40]、黃姑魚(yú)(Argyrosomusjaponicus)[41]等種群鑒別時(shí)效果良好,2種方法聯(lián)合使用在魚(yú)類(lèi)耳石的相關(guān)研究中可靠性更高。

    致謝:本文依托國(guó)家自然科學(xué)基金“許氏平鲉資源增殖效果評(píng)估新方法及其應(yīng)用基礎(chǔ)研究”。中國(guó)海洋大學(xué)漁業(yè)生態(tài)學(xué)實(shí)驗(yàn)室為本項(xiàng)目的開(kāi)展給予了大力支持與幫助,作者在此深表感謝。

    [1] 張曉霞,葉振江,王英俊,等. 青島海域小眼綠鰭魚(yú)耳石形態(tài)的初步研究[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,39(4): 622-626.

    Zhang Xiaoxia,Ye Zhenjiang,Wang Yingjun,et al. A preliminary study on otolith morphology ofChelidonichthysspinosusin the Qingdao coastal waters[J]. Periodical of Ocean University of China(Natural Science),2009,39(4): 622-626.

    [2] 秦巖,潘曉哲,高天翔. 鲬?hù)~(yú)耳石形態(tài)多樣性研究[J]. 水生生物學(xué)報(bào),2013,37(4): 759-765.

    Qin Yan,Pan Xiaozhe,Gao Tianxiang. Study on morphological diversity of sagittae ofPlatycephalusindicus[J]. Aquatic Biology,2013,37(4): 759-765.

    [3] Campana S E. Photographic atlas of fish otoliths of the Northwest Atlantic Ocean[M]. Ottawa:NRC Research Press,2004.

    [4] Fitch J E,Brownell Jr R L. Fish otoliths in cetacean stomachs and their importance in interpreting feeding habits[J]. Journal of the Fisheries Board of Canada,1968,25(12): 2561-2574.

    [5] Lychakov D V,Rebane Y T. Otolith regularities[J]. Hearing Research,2000,143(1): 83-102.

    [6] Vignon M. Ontogenetic trajectories of otolith shape during shift in habitat use: Interaction between otolith growth and environment[J]. Journal of Experimental Marine Biology and Ecology,2012,420: 26-32.

    [7] Gagliano M,McCormick M I. Feeding history influences otolith shape in tropical fish[J]. Marine Ecology Progress Series,2004,278: 291-296.

    [8] Vignon M,Morat F. Environmental and genetic determinant of otolith shape revealed by a non-indigenous tropical fish[J]. Marine Ecology Progress Series,2010,411: 231-241.

    [9] Agüera A,Brophy D. Use of saggital otolith shape analysis to discriminate Northeast Atlantic and Western Mediterranean stocks of Atlantic saury,Scomberesoxsaurussaurus(Walbaum)[J]. Fisheries Research,2011,110(3): 465-471.

    [10] Morales-Nin B. Review of the growth regulation processes of otolith daily increment formation[J]. Fisheries Research,2000,46(1): 53-67.

    [11] Maisey J G. Notes on the structure and phylogeny of vertebrate otoliths[J]. Copeia,1987(2): 495-499.

    [12] 葉振江. 中國(guó)海洋魚(yú)類(lèi)耳石形態(tài)學(xué)分析及應(yīng)用研究[D]. 青島:中國(guó)海洋大學(xué),2010.

    Ye Zhenjiang. Study on otolith shape,taxonomy and application for fishes in China sea[D]. Qingdao:Ocean University of China,2010.

    [13] Begg G A,Overholtz W J,Munroe N J. The use of internal otolith morphometrics for identification of haddock (Melanogrammusaeglefinus) stocks on Georges Bank[J]. Fishery Bulletin,2001,99(1):1-14.

    [14] Tuset V M,Rosin P L,Lombarte A. Sagittal otolith shape used in the identification of fishes of the genusSerranus[J]. Fisheries Research,2006,81(2): 316-325.

    [15] 張國(guó)華,但勝?lài)?guó),苗志國(guó),等. 六種鯉科魚(yú)類(lèi)耳石形態(tài)以及在種類(lèi)和群體識(shí)別中的應(yīng)用[J]. 水生生物學(xué)報(bào),1999,23(6): 683-688.

    Zhang Guohua,Dan Shengguo,Miao Zhiguo,et al. Otolith morphology of six Cyprinid species with the use in species and stock discriminations[J]. Acta Hydrobiological Sinica,1999,23(6): 683-688.

    [16] 郭弘藝,魏凱,唐文喬,等. 基于矢耳石形態(tài)特征的中國(guó)鱭屬魚(yú)類(lèi)種類(lèi)識(shí)別[J]. 動(dòng)物分類(lèi)學(xué)報(bào),2010 (1): 127-134.

    Guo Hongyi,Wei Kai,Tang Wenqiao,et al. Sibling species discrimination for Chinese genus ofCoiliafishes based on sagittal otolith morphology[J]. Acta Zootaxonomica Sinica,2010 (1): 127-134.

    [17] L’Abée-Lund J H. Otolith shape discriminates between juvenile Atlantic salmon,SalmosalarL.,and brown trout,SalmotruttaL.[J]. Journal of fish biology,1988,33(6): 899-903.

    [18] 曾祥波,張國(guó)華. 基于傅立葉耳石形態(tài)分析法對(duì)幾種家魚(yú)仔稚魚(yú)的種類(lèi)判別分析[J]. 中國(guó)水產(chǎn)科學(xué),2012,6: 970-977.

    Zeng Xiangbo,Zhang Guohua. Species identification at the larval and juvenile stages for several Chinese domestic fishes by elliptical Fourier analysis of otolith form[J]. Journal of Fishery Sciences of China,2012,6: 970-977.

    [19] Ponton D. Is geometric morphometrics efficient for comparing otolith shape of different fish species?[J]. Journal of Morphology,2006,267(6): 750-757.

    [20] 潘曉哲,高天翔. 基于耳石形態(tài)的鱚屬魚(yú)類(lèi)鑒別[J]. 動(dòng)物分類(lèi)學(xué)報(bào),2010,35(4): 799-805.

    Pan Xiaozhe,Gao Tianxiang. Sagittal otolith shape used in the discrimination of fishes of the GenusSillagoin China[J]. Acta Zootaxonomica Sinica,2010,35(4): 799-805.

    [21] Lord C,Morat F,Lecomte-Finiger R,et al. Otolith shape analysis for threeSicyopterus(Teleostei: Gobioidei: Sicydiinae) species from New Caledonia and Vanuatu[J]. Environmental Biology of Fishes,2012,93(2): 209-222.

    [22] Tuset V M,Lozano I J,Gonzalez J A,et al. Shape indices to identify regional differences in otolith morphology of comber,Serranuscabrilla(L.,1758)[J]. Journal of Applied Ichthyology,2003,19(2): 88-93.

    [23] Stransky C,Murta A G,Schlickeisen J,et al. Otolith shape analysis as a tool for stock separation of horse mackerel (Trachurustrachurus) in the Northeast Atlantic and Mediterranean[J]. Fisheries Research,2008,89(2): 159-166.

    [24] Tracey S R,Lyle J M,Duhamel G. Application of elliptical Fourier analysis of otolith form as a tool for stock identification[J]. Fisheries Research,2006,77(2): 138-147.

    [25] Begg G A,Brown R W. Stock identification of haddockMelanogrammusaeglefinuson Georges Bank based on otolith shape analysis[J]. Transactions of the American Fisheries Society,2000,129(4): 935-945.

    [26] Galley E A,Wright P J,Gibb F M. Combined methods of otolith shape analysis improve identification of spawning areas of Atlantic cod[J]. ICES Journal of Marine Science: Journal du Conseil,2006,63(9): 1710-1717.

    [27] 徐勝勇,張輝,潘曉哲,等. 許氏平鲉群體耳石形態(tài)學(xué)比較[J]. 中國(guó)海洋大學(xué)學(xué)報(bào) (自然科學(xué)版),2012,42(11):54-61.

    Xu Shengyong,Zhang Hui,Pan Xiaozhe,et al. Morphological comparison of the otolith ofSebastesschlegeliipopulations[J]. Periodical of Ocean University of China(Natural Science),2012,42(11):54-61.

    [28] Tuset V M,Lombarte A,Gonzalez J A,et al. Comparative morphology of the sagittal otolith inSerranusspp.[J]. Journal of fish biology,2003,63(6): 1491-1504.

    [30] 王英俊. 傅里葉分析在魚(yú)類(lèi)耳石形態(tài)學(xué)中的應(yīng)用研究[D]. 青島:中國(guó)海洋大學(xué),2010.

    Wang Yingjun. The application of Fourier Analysis in the research of otolith morphology[D]. Qingdao:Ocean University of China,2010.

    [31] Kuhl F P,Giardina C R. Elliptic Fourier features of a closed contour[J]. Computer Graphics and Image Processing,1982,18(3): 236-258.

    [32] 王英俊,葉振江,劉群,等. 細(xì)條天竺魚(yú) (Apogonichthyslineatus) 與黑鰓天竺魚(yú) (Apogonichthysarafurae) 耳石形態(tài)識(shí)別的初步研究[J]. 海洋與湖沼,2010,41(2): 281-284.

    Wang Yingjun,Ye Zhenjiang,Liu Qun,et al. Using otolith to identify two fish species:ApogonichthyslineatusandApogonichthysarafurae[J]. Oceanologica et Limnological Sinica,2010,41(2): 281-284.

    [33] 竇碩增,于鑫,曹亮. 魚(yú)類(lèi)矢耳石形態(tài)分析及其在群體識(shí)別中的應(yīng)用實(shí)例研究[J]. 海洋與湖沼,2012,43(4): 702-712.

    Dou Shuozeng,Yu Xin,Cao Liang. Otolith shape analysis and its application in fish stock discrimination: a case study[J]. Oceanologia et Limnologia Sinica,2012,43(4): 702-712.

    [34] 王英俊,葉振江,張弛,等. 車(chē)牛山島鎧平鲉繁殖群體生物學(xué)研究[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,41(12): 46-52.

    Wang Yingjun,Ye Zhenjiang,Zhang Chi,et al. The fisheries biology ofSebasteshubbsiin Che-Niushan island[J]. Periodical of Ocean University of China(Natural Science),2011,41(12): 46-52.

    [35] 張輝. 西北太平洋兩種卵胎生魚(yú)類(lèi) (許氏平鲉和褐菖鲉) 的分子系統(tǒng)地理學(xué)研究[D]. 青島:中國(guó)海洋大學(xué),2013.

    Zhang Hui. Molecular phylogeography of two marine ovoviviparous fishes in northwestern Pacific[D]. Qingdao:Ocean University of China,2013.

    [36] 陸化杰,陳新軍,劉必林. 個(gè)體差異對(duì)西南大西洋阿根廷滑柔魚(yú)耳石形態(tài)的影響[J]. 水產(chǎn)學(xué)報(bào),2011,35(2): 247-253.

    Lu Huajie,Chen Xinjun,Liu Bilin. Effects of individual difference on statolith morphology ofIllexargentinusin the southwest Atlantic Ocean[J]. Journal of Fisheries of China,2011,35(2): 247-253.

    [37] 于鑫,曹亮,南鷗,等. 基于矢耳石形態(tài)分析的鳳鱭 (Coiliamystus) 群體識(shí)別研究[J]. 海洋與湖沼,2013,44(3):768-774.

    Yu Xin,Cao Liang,Nan Ou,et al. Stock identification ofCoiliamystususing otolith shape analysis[J]. Oceanologia et Limnologia Sinica,2013,44(3):768-774.

    [38] 李輝華,郭弘藝,唐文喬,等. 兩種耳石分析法在鱭屬種間和種群間識(shí)別效果的比較研究[J]. 淡水漁業(yè),2013,43(1): 14-18.

    Li Huihua,Guo Hongyi,Tang Wenqiao,et al. Comparative study of two otolith shape analysis for genusCoillaspecies and stocks identification[J]. Freshwater Fisheries,2013,43(1): 14-18.

    [39] DeVries D A,Grimes C B,Prager M H. Using otolith shape analysis to distinguish eastern Gulf of Mexico and Atlantic Ocean stocks of king mackerel[J]. Fisheries Research,2002,57(1): 51-62.

    [40] Duarte-Neto P,Lessa R,Stosic B,et al. The use of sagittal otoliths in discriminating stocks of common dolphinfish (Coryphaenahippurus) off northeastern Brazil using multishape descriptors[J]. ICES Journal of Marine Science: Journal du Conseil,2008,65(7): 1144-1152.

    [41] Ferguson G J,Ward T M,Gillanders B M. Otolith shape and elemental composition: Complementary tools for stock discrimination of mulloway (Argyrosomusjaponicus) in southern Australia[J]. Fisheries Research,2011,110(1): 75-83.

    Sagittae morphology used in the discrimination of the fish ofSebastesand the comparison of sagittal research method

    Zhang Fengxia1,Zhang Xiumei1,Qin Lezheng1,Li Wentao1

    (1.TheKeyLaboratoryofMariculture,MinistryofEducation,OceanUniversityofChina,Qingdao266003,China)

    The sagittae morphological characteristics of three speciesSebasteshubbsi,S.nigricansandS.schlegeliifrom Jiaozhou Bay,Qingdao,were analyzed to test the hypothesis that the threeSebastesspecies could be identified by their sagittae morphology and to find out which method was more appropriate. Eight otolith shape indices and 37 elliptical Fourier coefficients (eFcs) were analyzed,and results showed that the difference in the sagittae morphology betweenS.hubbsiandS.schlegeliiwere the largest,and those betweenS.hubbsiandS.nigricanswere the least. In the discriminant analyses of the three species,discriminating rates of more than 95% were obtained and they were clearly sorted in scatter plots when using the otolith shape indices,the eFcs,as well as a combination of the two. When applied to the three species of which otolith was 4 to 7 mm and those were more than two years old,the species identification based on the otolith morphology was proved to be feasible: when only using the eFcs,a 100% of discriminating rate was obtained but this method couldn’t quickly locate the shape which was tiny in fact but significant difference among the sagittae of the threeSebastesin statistics. Combing otolith shape indices and the eFcs could not only obtain a high classification rate of about 99%,but also identify the differences in the rectangularity,form-factor and surface density among the threeSebastessagittae. The results showed that combing shape indices analysis and elliptical Fourier analysis was the effective method for species discrimination of genusSebastesand it can find out the differences of the sagittae of the threeSebastes.

    Sebastes; sagittae morphology; otolith shape indices; elliptical Fourier coefficients; species discrimination

    10.3969/j.issn.0253-4193.2015.04.003

    2014-04-20;

    2014-12-28。

    國(guó)家自然科學(xué)基金(31172447,41176117);國(guó)家海洋公益性行業(yè)科研專(zhuān)項(xiàng)(201305043)。

    張鳳俠(1990—),女,安徽省阜陽(yáng)市人,從事資源生物學(xué)研究。 Email: fxzhang1003@163.com

    *通信作者:張秀梅,教授。E-mail: xiumei1227@163.com

    S917.4

    A

    0253-4193(2015)04-0028-12

    張鳳俠, 張秀梅, 覃樂(lè)政,等.基于矢耳石形態(tài)的平鲉屬(Sebastes)魚(yú)類(lèi)鑒別及研究方法比較[J].海洋學(xué)報(bào),2015,37(4):28—39,

    Zhang Fengxia, Zhang Xiumei, Qin Lezheng, et al. Sagittae morphology used in the discrimination of the fish ofSebastesand the comparison of sagittal research method[J]. Haiyang Xuebao,2015,37(4):28—39, doi: 10.3969/j.issn.0253-4193.2015.04.003

    猜你喜歡
    許氏判別函數(shù)耳石
    耳石癥復(fù)位后避免側(cè)臥
    水溫變化對(duì)近海網(wǎng)箱養(yǎng)殖中許氏平鲉攝食及生長(zhǎng)的影響
    耳石癥患者后管復(fù)位治療成功后頭暈癥狀臨床觀(guān)察
    游樂(lè)設(shè)施事故與危險(xiǎn)量化判別函數(shù)的構(gòu)建
    你了解耳石癥嗎
    老友(2018年7期)2018-08-08 15:32:30
    探究上市公司財(cái)務(wù)預(yù)警的數(shù)學(xué)模型
    名醫(yī)妙用單方治好胃下垂
    耳朵生病也會(huì)眩暈
    基于Fisher判別函數(shù)的酒店員工離職預(yù)警研究
    許氏平鲉早期異速生長(zhǎng)模式的研究
    亚洲av男天堂| 国产成人免费无遮挡视频| 天天添夜夜摸| 欧美日韩av久久| 国产精品成人在线| 国产又色又爽无遮挡免| 欧美日韩av久久| 老司机午夜福利在线观看视频 | av一本久久久久| 黄色视频,在线免费观看| 亚洲精品国产av蜜桃| 一区二区三区四区激情视频| 国产精品成人在线| 深夜精品福利| 黄色 视频免费看| 国产日韩一区二区三区精品不卡| 电影成人av| 天天影视国产精品| 精品少妇一区二区三区视频日本电影| 黑人猛操日本美女一级片| 国产免费福利视频在线观看| 亚洲国产欧美日韩在线播放| 国产成人精品久久二区二区免费| 亚洲伊人久久精品综合| 丝袜美腿诱惑在线| 国产在线观看jvid| 国产精品久久久久久精品电影小说| 日韩一卡2卡3卡4卡2021年| 欧美激情极品国产一区二区三区| 亚洲国产成人一精品久久久| 午夜福利影视在线免费观看| 精品欧美一区二区三区在线| 老鸭窝网址在线观看| 99国产精品免费福利视频| 日韩中文字幕欧美一区二区| 色综合欧美亚洲国产小说| 在线看a的网站| 精品亚洲乱码少妇综合久久| 91精品国产国语对白视频| 久久这里只有精品19| 日本a在线网址| 日韩精品免费视频一区二区三区| 国产精品一区二区精品视频观看| 国产精品 欧美亚洲| 在线观看免费午夜福利视频| 中文字幕人妻丝袜一区二区| 天堂中文最新版在线下载| 久久综合国产亚洲精品| 亚洲国产精品成人久久小说| 亚洲国产精品成人久久小说| 两个人看的免费小视频| 亚洲精品久久久久久婷婷小说| 91成人精品电影| 日韩 亚洲 欧美在线| 国产人伦9x9x在线观看| 国产成人av教育| 久久人妻熟女aⅴ| av欧美777| 精品国产超薄肉色丝袜足j| 日本av手机在线免费观看| av福利片在线| 国产精品一区二区在线观看99| 精品人妻一区二区三区麻豆| 不卡av一区二区三区| 国产高清videossex| 精品熟女少妇八av免费久了| 亚洲av电影在线进入| av视频免费观看在线观看| 亚洲国产成人一精品久久久| 国产精品久久久久久精品电影小说| av在线播放精品| 大码成人一级视频| 国产熟女午夜一区二区三区| 精品熟女少妇八av免费久了| 国产亚洲精品久久久久5区| 久久天堂一区二区三区四区| 啦啦啦啦在线视频资源| 99九九在线精品视频| 国产一区二区 视频在线| 国产成人免费观看mmmm| 日韩有码中文字幕| 中文欧美无线码| 国产精品九九99| 日本欧美视频一区| 亚洲欧美激情在线| 99久久综合免费| 色综合欧美亚洲国产小说| 国产精品成人在线| 精品免费久久久久久久清纯 | 操美女的视频在线观看| 欧美精品一区二区大全| 黄色视频不卡| 久久久精品免费免费高清| 女人高潮潮喷娇喘18禁视频| 国产无遮挡羞羞视频在线观看| 777久久人妻少妇嫩草av网站| 黄色怎么调成土黄色| 成年女人毛片免费观看观看9 | 老司机午夜十八禁免费视频| 精品福利永久在线观看| 午夜福利乱码中文字幕| av免费在线观看网站| 国产成人a∨麻豆精品| 亚洲国产中文字幕在线视频| 亚洲人成电影免费在线| 五月开心婷婷网| 亚洲第一av免费看| 中文字幕制服av| 天天躁日日躁夜夜躁夜夜| 黑人猛操日本美女一级片| 亚洲伊人久久精品综合| 亚洲精品久久久久久婷婷小说| 午夜福利免费观看在线| 美女午夜性视频免费| 蜜桃国产av成人99| 亚洲欧美清纯卡通| 国产人伦9x9x在线观看| 精品欧美一区二区三区在线| 亚洲天堂av无毛| 亚洲av日韩在线播放| 99国产精品免费福利视频| 国产主播在线观看一区二区| 999久久久精品免费观看国产| 手机成人av网站| 50天的宝宝边吃奶边哭怎么回事| 国产精品 欧美亚洲| 亚洲七黄色美女视频| 久久亚洲国产成人精品v| 嫩草影视91久久| 亚洲全国av大片| 麻豆乱淫一区二区| 欧美 亚洲 国产 日韩一| 啦啦啦视频在线资源免费观看| 国产激情久久老熟女| 国产精品欧美亚洲77777| 日本wwww免费看| 制服诱惑二区| 成人黄色视频免费在线看| 亚洲精华国产精华精| 天堂中文最新版在线下载| 亚洲国产精品999| 国产不卡av网站在线观看| 日韩 欧美 亚洲 中文字幕| av在线老鸭窝| 老司机靠b影院| 狂野欧美激情性bbbbbb| 午夜两性在线视频| 亚洲欧美日韩另类电影网站| 黑人巨大精品欧美一区二区mp4| av免费在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产中文字幕在线视频| 一区在线观看完整版| 成人手机av| 丝袜人妻中文字幕| 久久久久国内视频| a级毛片在线看网站| 国产成人一区二区三区免费视频网站| 亚洲精品国产精品久久久不卡| 国产精品久久久人人做人人爽| 少妇精品久久久久久久| 国产在线一区二区三区精| 嫩草影视91久久| 三上悠亚av全集在线观看| 女人被躁到高潮嗷嗷叫费观| 国产精品国产三级国产专区5o| 国产精品久久久久久精品电影小说| 国产免费一区二区三区四区乱码| 国产高清国产精品国产三级| 狂野欧美激情性xxxx| 久久香蕉激情| 亚洲美女黄色视频免费看| 亚洲精品日韩在线中文字幕| 亚洲欧美一区二区三区黑人| 亚洲人成77777在线视频| 国产成人精品在线电影| 免费观看人在逋| 最近最新中文字幕大全免费视频| 1024视频免费在线观看| 1024视频免费在线观看| 18在线观看网站| 最近最新免费中文字幕在线| 青草久久国产| 久久久国产成人免费| 国产国语露脸激情在线看| 久9热在线精品视频| 在线观看免费午夜福利视频| 美女主播在线视频| 欧美日韩一级在线毛片| 日本vs欧美在线观看视频| 色婷婷av一区二区三区视频| 18禁黄网站禁片午夜丰满| 午夜老司机福利片| 三级毛片av免费| 99久久精品国产亚洲精品| 久久这里只有精品19| 亚洲国产av影院在线观看| 中文字幕人妻熟女乱码| 美女主播在线视频| 成年人黄色毛片网站| 中国美女看黄片| 十八禁人妻一区二区| 91国产中文字幕| 别揉我奶头~嗯~啊~动态视频 | 女性被躁到高潮视频| 女性被躁到高潮视频| 国产精品99久久99久久久不卡| av不卡在线播放| 国产精品久久久久久精品电影小说| 亚洲一区中文字幕在线| av免费在线观看网站| 国产精品一二三区在线看| 悠悠久久av| av在线app专区| 日韩精品免费视频一区二区三区| 国产又色又爽无遮挡免| 黄片小视频在线播放| 久久久欧美国产精品| 丝袜美腿诱惑在线| 亚洲精华国产精华精| 法律面前人人平等表现在哪些方面 | 人人澡人人妻人| 久久青草综合色| 一级毛片精品| a 毛片基地| 建设人人有责人人尽责人人享有的| 少妇 在线观看| 婷婷色av中文字幕| 自线自在国产av| 男人添女人高潮全过程视频| 国产免费av片在线观看野外av| 精品少妇黑人巨大在线播放| 久久久国产精品麻豆| 色94色欧美一区二区| 99国产精品一区二区蜜桃av | 男男h啪啪无遮挡| 无限看片的www在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 乱人伦中国视频| 波多野结衣av一区二区av| 久热这里只有精品99| 成年人免费黄色播放视频| 人人妻人人澡人人爽人人夜夜| 99国产综合亚洲精品| 亚洲精华国产精华精| 精品一区二区三区四区五区乱码| e午夜精品久久久久久久| 久久久久久亚洲精品国产蜜桃av| 欧美中文综合在线视频| 亚洲情色 制服丝袜| 99久久综合免费| 一二三四在线观看免费中文在| 久久久久久免费高清国产稀缺| 韩国精品一区二区三区| 欧美av亚洲av综合av国产av| 国产91精品成人一区二区三区 | a 毛片基地| 一区二区三区激情视频| 91字幕亚洲| 亚洲欧美色中文字幕在线| 精品人妻在线不人妻| 精品乱码久久久久久99久播| 在线天堂中文资源库| 成在线人永久免费视频| 久久精品亚洲av国产电影网| 欧美黄色淫秽网站| 淫妇啪啪啪对白视频 | 久久久久久免费高清国产稀缺| 老熟女久久久| 少妇被粗大的猛进出69影院| 久久中文字幕一级| 伊人亚洲综合成人网| 国产成人一区二区三区免费视频网站| 国产亚洲午夜精品一区二区久久| 亚洲 欧美一区二区三区| 成年人午夜在线观看视频| 国产深夜福利视频在线观看| 国产精品国产av在线观看| 久久女婷五月综合色啪小说| 俄罗斯特黄特色一大片| 51午夜福利影视在线观看| 国产欧美亚洲国产| 飞空精品影院首页| 精品国产一区二区久久| 乱人伦中国视频| 欧美亚洲 丝袜 人妻 在线| 亚洲av国产av综合av卡| 亚洲天堂av无毛| 色综合欧美亚洲国产小说| 少妇猛男粗大的猛烈进出视频| 色老头精品视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产主播在线观看一区二区| 欧美日本中文国产一区发布| 亚洲av欧美aⅴ国产| 极品人妻少妇av视频| 青青草视频在线视频观看| 一边摸一边抽搐一进一出视频| 久久人人97超碰香蕉20202| 欧美午夜高清在线| 免费女性裸体啪啪无遮挡网站| 操出白浆在线播放| av天堂久久9| 99久久99久久久精品蜜桃| 性色av一级| 精品人妻一区二区三区麻豆| 老汉色av国产亚洲站长工具| 99精国产麻豆久久婷婷| 超色免费av| 各种免费的搞黄视频| 久久久久久久精品精品| 一边摸一边抽搐一进一出视频| 超碰97精品在线观看| 精品亚洲成a人片在线观看| 巨乳人妻的诱惑在线观看| 老汉色av国产亚洲站长工具| 亚洲国产精品一区二区三区在线| tocl精华| 法律面前人人平等表现在哪些方面 | 三级毛片av免费| 免费久久久久久久精品成人欧美视频| 欧美中文综合在线视频| 日日摸夜夜添夜夜添小说| 手机成人av网站| 国产欧美日韩一区二区三区在线| 考比视频在线观看| 成年美女黄网站色视频大全免费| 午夜福利在线观看吧| 婷婷丁香在线五月| 国产高清国产精品国产三级| 亚洲欧美清纯卡通| 999久久久精品免费观看国产| 一区二区日韩欧美中文字幕| 成人亚洲精品一区在线观看| 后天国语完整版免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看免费高清a一片| 激情视频va一区二区三区| 亚洲国产看品久久| 午夜福利免费观看在线| 18在线观看网站| 成人三级做爰电影| 搡老熟女国产l中国老女人| 妹子高潮喷水视频| 乱人伦中国视频| 精品少妇黑人巨大在线播放| 一区二区三区精品91| av网站在线播放免费| 天天躁日日躁夜夜躁夜夜| 亚洲欧美激情在线| 啦啦啦免费观看视频1| 在线天堂中文资源库| 如日韩欧美国产精品一区二区三区| 精品欧美一区二区三区在线| 国产免费一区二区三区四区乱码| 狠狠精品人妻久久久久久综合| 80岁老熟妇乱子伦牲交| 中文字幕人妻丝袜制服| 热99re8久久精品国产| 久久综合国产亚洲精品| 一二三四在线观看免费中文在| 精品人妻1区二区| 女人久久www免费人成看片| 香蕉国产在线看| 法律面前人人平等表现在哪些方面 | 9热在线视频观看99| 国产精品影院久久| 欧美黄色片欧美黄色片| 亚洲av国产av综合av卡| 丰满人妻熟妇乱又伦精品不卡| 老司机亚洲免费影院| 精品亚洲乱码少妇综合久久| 啦啦啦 在线观看视频| 亚洲人成电影免费在线| 国产99久久九九免费精品| 午夜两性在线视频| 国产欧美日韩一区二区三 | 欧美日韩一级在线毛片| 久久热在线av| 性色av一级| 日本精品一区二区三区蜜桃| 亚洲一区二区三区欧美精品| 欧美另类亚洲清纯唯美| 亚洲自偷自拍图片 自拍| 女人高潮潮喷娇喘18禁视频| 午夜福利免费观看在线| 精品一区二区三区四区五区乱码| 99久久99久久久精品蜜桃| 日韩制服丝袜自拍偷拍| 国产男女超爽视频在线观看| 亚洲伊人久久精品综合| 中文字幕精品免费在线观看视频| 日本欧美视频一区| 亚洲激情五月婷婷啪啪| 黄片播放在线免费| 亚洲全国av大片| 国产男女超爽视频在线观看| 日本vs欧美在线观看视频| 亚洲成人手机| 97在线人人人人妻| 免费观看av网站的网址| 欧美日韩亚洲国产一区二区在线观看 | 国产男女内射视频| 国产精品熟女久久久久浪| 国产在线观看jvid| 国产亚洲精品久久久久5区| 久久久久视频综合| 欧美国产精品va在线观看不卡| 日本欧美视频一区| 电影成人av| 国产日韩一区二区三区精品不卡| 中文字幕制服av| 天天添夜夜摸| 成人国语在线视频| 夜夜骑夜夜射夜夜干| 美女中出高潮动态图| 美女大奶头黄色视频| 国产精品欧美亚洲77777| 国产精品国产三级国产专区5o| 亚洲久久久国产精品| 国产黄频视频在线观看| 久久99热这里只频精品6学生| 91成人精品电影| 亚洲欧美精品自产自拍| 国产精品 欧美亚洲| 黄网站色视频无遮挡免费观看| 日韩,欧美,国产一区二区三区| 国产精品秋霞免费鲁丝片| 国产精品秋霞免费鲁丝片| 老司机午夜福利在线观看视频 | 一区二区av电影网| 91九色精品人成在线观看| 亚洲少妇的诱惑av| 久久精品人人爽人人爽视色| 精品久久久久久电影网| 亚洲中文日韩欧美视频| 涩涩av久久男人的天堂| 久久99一区二区三区| 国产亚洲欧美精品永久| 亚洲成人免费电影在线观看| 欧美日本中文国产一区发布| av福利片在线| 精品高清国产在线一区| 自线自在国产av| 亚洲性夜色夜夜综合| 99热国产这里只有精品6| 大码成人一级视频| 欧美性长视频在线观看| av欧美777| 久久性视频一级片| 91老司机精品| 极品少妇高潮喷水抽搐| 亚洲一码二码三码区别大吗| 亚洲精品国产一区二区精华液| 久久久精品区二区三区| 色94色欧美一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 妹子高潮喷水视频| 人人妻人人澡人人爽人人夜夜| 多毛熟女@视频| 亚洲 国产 在线| 男女免费视频国产| 91大片在线观看| 天天躁夜夜躁狠狠躁躁| 中文字幕另类日韩欧美亚洲嫩草| 黑人猛操日本美女一级片| 啦啦啦啦在线视频资源| 国产成人精品在线电影| 夫妻午夜视频| 黄色 视频免费看| 久久人妻熟女aⅴ| 丰满人妻熟妇乱又伦精品不卡| 国产无遮挡羞羞视频在线观看| 日韩熟女老妇一区二区性免费视频| 老司机福利观看| 黄网站色视频无遮挡免费观看| 91精品三级在线观看| 亚洲天堂av无毛| av电影中文网址| 欧美老熟妇乱子伦牲交| 性少妇av在线| 亚洲精品一卡2卡三卡4卡5卡 | 母亲3免费完整高清在线观看| 少妇猛男粗大的猛烈进出视频| 日韩视频在线欧美| 一区在线观看完整版| 在线 av 中文字幕| 国产欧美日韩精品亚洲av| 国产男人的电影天堂91| 妹子高潮喷水视频| 亚洲黑人精品在线| 免费在线观看日本一区| 美女中出高潮动态图| 国产成人一区二区三区免费视频网站| 久久久久网色| 高清视频免费观看一区二区| 国产区一区二久久| 国产精品国产三级国产专区5o| 国产精品自产拍在线观看55亚洲 | 新久久久久国产一级毛片| a级毛片在线看网站| 亚洲av片天天在线观看| 在线永久观看黄色视频| 日本撒尿小便嘘嘘汇集6| 亚洲精品中文字幕在线视频| 欧美xxⅹ黑人| 国产av国产精品国产| 国产精品国产av在线观看| 国产精品亚洲av一区麻豆| 欧美亚洲日本最大视频资源| 狠狠婷婷综合久久久久久88av| 久久国产精品人妻蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 中文欧美无线码| 一边摸一边抽搐一进一出视频| 久久精品人人爽人人爽视色| 成人影院久久| 少妇精品久久久久久久| 午夜精品久久久久久毛片777| 丝袜脚勾引网站| 黄色视频,在线免费观看| 欧美另类一区| 国产在线一区二区三区精| 在线观看免费日韩欧美大片| 在线 av 中文字幕| 午夜福利免费观看在线| 久久人人爽av亚洲精品天堂| 欧美在线黄色| 精品一区二区三区四区五区乱码| 在线观看www视频免费| 黄频高清免费视频| 窝窝影院91人妻| 他把我摸到了高潮在线观看 | bbb黄色大片| 国产精品国产三级国产专区5o| 高清黄色对白视频在线免费看| 国产成人精品久久二区二区免费| 久久中文看片网| 国产男女内射视频| 亚洲成人免费av在线播放| 中文字幕人妻熟女乱码| 亚洲精品成人av观看孕妇| 久久人妻熟女aⅴ| 精品少妇内射三级| 日韩欧美一区二区三区在线观看 | 中文字幕精品免费在线观看视频| 啦啦啦 在线观看视频| 中文字幕最新亚洲高清| 精品一区二区三卡| 首页视频小说图片口味搜索| 伊人久久大香线蕉亚洲五| 黄片大片在线免费观看| www日本在线高清视频| 久久女婷五月综合色啪小说| 天堂8中文在线网| 天天躁狠狠躁夜夜躁狠狠躁| 99热国产这里只有精品6| 丰满饥渴人妻一区二区三| 亚洲精品中文字幕一二三四区 | 午夜激情久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成电影观看| av不卡在线播放| 国产在线免费精品| 亚洲情色 制服丝袜| avwww免费| 最新在线观看一区二区三区| 悠悠久久av| 亚洲熟女精品中文字幕| 国产精品一区二区精品视频观看| 51午夜福利影视在线观看| 超碰成人久久| 俄罗斯特黄特色一大片| 91精品伊人久久大香线蕉| 午夜福利,免费看| 男女国产视频网站| 午夜免费观看性视频| 欧美一级毛片孕妇| 日本91视频免费播放| 亚洲精品美女久久av网站| 成人三级做爰电影| 日本欧美视频一区| 老司机午夜福利在线观看视频 | 亚洲精品国产av成人精品| 成人亚洲精品一区在线观看| 中文字幕高清在线视频| 老司机影院毛片| 午夜福利影视在线免费观看| 老司机影院毛片| 久久国产亚洲av麻豆专区| www.999成人在线观看| 久久99热这里只频精品6学生| 国产不卡av网站在线观看| 成人18禁高潮啪啪吃奶动态图| av视频免费观看在线观看| 2018国产大陆天天弄谢| 亚洲欧美色中文字幕在线| 久久久久视频综合| 一个人免费看片子| 自拍欧美九色日韩亚洲蝌蚪91| 搡老岳熟女国产| 国产亚洲午夜精品一区二区久久| 精品人妻在线不人妻| 日本猛色少妇xxxxx猛交久久| 巨乳人妻的诱惑在线观看| 亚洲免费av在线视频| 欧美亚洲日本最大视频资源| 国产成人精品久久二区二区91| 99国产精品99久久久久| 蜜桃在线观看..| 国产成人欧美在线观看 | 免费在线观看影片大全网站| 黑丝袜美女国产一区| netflix在线观看网站| 99精国产麻豆久久婷婷|