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

    長(zhǎng)牡蠣殼橙品系形態(tài)性狀與體質(zhì)量的相關(guān)及通徑分析*

    2017-11-01 15:02:31韓自強(qiáng)
    關(guān)鍵詞:通徑活體牡蠣

    韓自強(qiáng), 李 琪

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

    長(zhǎng)牡蠣殼橙品系形態(tài)性狀與體質(zhì)量的相關(guān)及通徑分析*

    韓自強(qiáng), 李 琪**

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

    為研究長(zhǎng)牡蠣(Crassostreagigas)殼橙品系的殼形態(tài)性狀對(duì)體質(zhì)量性狀的影響。本實(shí)驗(yàn)隨機(jī)選取144只二齡殼橙品系長(zhǎng)牡蠣為研究對(duì)象,測(cè)量了殼高(XH)、殼長(zhǎng)(XL)、殼寬(XW)、活體總重(YL)和軟體部重(YS),以殼形態(tài)性狀(XH,XL,XW)為自變量,體質(zhì)量性狀(YL,YS)為因變量進(jìn)行了相關(guān)分析、通徑分析和多元回歸分析。相關(guān)分析結(jié)果表明,所測(cè)量的長(zhǎng)牡蠣橙色品系各性狀間均呈極顯著的表型相關(guān)(P<0.01)。通徑分析結(jié)果表明,殼形態(tài)性狀對(duì)活體總重直接影響依次為殼高(0.477)>殼寬(0.336)>殼長(zhǎng)(0.326),殼高對(duì)活體總重的直接決定系數(shù)最大(0.228),是影響活體總重的主要因素;然而,軟體部重的相關(guān)系數(shù)R2<0.850,表明可能還存在對(duì)軟體部重影響較大的其他因素。通過(guò)多元回歸分析,建立了殼形態(tài)性狀對(duì)體質(zhì)量性狀的最優(yōu)回歸方程。

    長(zhǎng)牡蠣;殼橙品系;殼形態(tài)性狀;體質(zhì)量性狀;相關(guān)分析;通徑分析;多元回歸分析

    長(zhǎng)牡蠣(Crassostreagigas)又稱(chēng)太平洋牡蠣,具有生長(zhǎng)速度快、抗逆性強(qiáng)、營(yíng)養(yǎng)豐富和肉質(zhì)鮮美等特點(diǎn),是世界上養(yǎng)殖產(chǎn)量最大的經(jīng)濟(jì)貝類(lèi)。2015年中國(guó)牡蠣養(yǎng)殖總產(chǎn)量達(dá)457.3萬(wàn)t,占全國(guó)所有海水養(yǎng)殖品產(chǎn)量的24.38%[1],牡蠣養(yǎng)殖在我國(guó)海水養(yǎng)殖業(yè)中占有重要地位。為保障長(zhǎng)牡蠣養(yǎng)殖業(yè)的健康發(fā)展,國(guó)內(nèi)學(xué)者開(kāi)展了長(zhǎng)牡蠣遺傳育種研究。殼色作為一種影響消費(fèi)者選擇的重要因素,逐漸成為貝類(lèi)育種關(guān)注的熱點(diǎn)。目前,海產(chǎn)貝類(lèi)中結(jié)合殼色和快速生長(zhǎng)已選育出多個(gè)殼色新品種[2-5],長(zhǎng)牡蠣新殼色快速生長(zhǎng)品系的選育是改善長(zhǎng)牡蠣品質(zhì)和提高長(zhǎng)牡蠣商品價(jià)值的一條重要途徑。在長(zhǎng)期的育種實(shí)踐中,我們選育出金黃色、白色、黑色和紫色4種長(zhǎng)牡蠣殼色品系[6],通過(guò)黑殼色和紫殼色長(zhǎng)牡蠣雜交形成紫黑殼色長(zhǎng)牡蠣個(gè)體,并在紫黑殼色長(zhǎng)牡蠣個(gè)體的自交后代中發(fā)現(xiàn)橙殼色突變個(gè)體,結(jié)合家系和群體選育技術(shù),培育出穩(wěn)定遺傳的殼橙長(zhǎng)牡蠣品系,具有重要的經(jīng)濟(jì)和科研價(jià)值。

    貝類(lèi)的體質(zhì)量性狀作為最直接的育種目標(biāo),相比殼高、殼長(zhǎng)、殼寬等殼形態(tài)性狀不夠直觀,因此需要通過(guò)多元回歸分析和通徑分析查明殼形態(tài)性狀(殼高、殼長(zhǎng)和殼寬)與體質(zhì)量性狀(活體總重和軟體部重)之間的相關(guān)性,以及對(duì)體質(zhì)量性狀的影響和決定程度[7]。目前,多元回歸與通徑分析在魚(yú)類(lèi)[8-9]、蝦類(lèi)[10]、蟹類(lèi)[11]、貝類(lèi)[12-17]等中已有較多報(bào)道,并且雙殼貝類(lèi)的同一物種的不同地理群體[18]、不同養(yǎng)殖方式群體[19]、不同性別群體[20]、不同生長(zhǎng)階段群體[21-22]、不同殼色群體[23-25]的殼形態(tài)性狀對(duì)體質(zhì)量性狀影響程度均表現(xiàn)出差異。

    在長(zhǎng)牡蠣生長(zhǎng)性狀方面,王慶志等比較了中國(guó)、日本和韓國(guó)3個(gè)地理群體快速生長(zhǎng)選育系的生長(zhǎng)性狀、殼型指數(shù)和顏色形狀的差異[26],以及長(zhǎng)牡蠣成體殼高、殼長(zhǎng)、殼寬、殼重、肉重和出肉率的遺傳參數(shù)[27];孔寧等研究了長(zhǎng)牡蠣快速生長(zhǎng)選育群體不同時(shí)期生長(zhǎng)性狀的發(fā)育模型[28];林青等比較了長(zhǎng)牡蠣和福建牡蠣(葡萄牙牡蠣)殼形態(tài)性狀對(duì)活體總重的影響[29]。然而,關(guān)于長(zhǎng)牡蠣殼色品系的殼形態(tài)性狀對(duì)體質(zhì)量性狀影響程度的分析還未見(jiàn)報(bào)道。本研究中,以長(zhǎng)牡蠣殼橙新品系為研究對(duì)象,通過(guò)通徑分析和多元回歸分析,研究殼形態(tài)性狀對(duì)體質(zhì)量性狀的直接和間接影響,旨在為長(zhǎng)牡蠣殼橙優(yōu)良品系的選育提供基礎(chǔ)資料。

    1 材料與方法

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

    2016年3月,于山東省榮成市愛(ài)蓮灣海區(qū)采集2齡長(zhǎng)牡蠣殼橙品系第4代選育群體作為實(shí)驗(yàn)材料,進(jìn)行殼形態(tài)性狀和體質(zhì)量性狀的測(cè)量。

    1.2 樣品測(cè)量

    隨機(jī)選取144個(gè)殼橙品系長(zhǎng)牡蠣,清除殼表面附著物,用紗布吸收殼內(nèi)外表面水分,使用數(shù)顯游標(biāo)卡尺(精確到0.01mm)測(cè)量其殼高(XH)、殼長(zhǎng)(XL)、殼寬(XW),用電子天平(精確到0.01g)測(cè)量活體總重(YL)和殼質(zhì)量。

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

    使用Excel對(duì)數(shù)據(jù)進(jìn)行整理,計(jì)算軟體部重(YS)=活體總重(YL)殼質(zhì)量。參考杜家菊等[30]和孫澤偉[16]的方法,使用軟件SPSS19.0統(tǒng)計(jì)殼形態(tài)性狀和體質(zhì)量性狀的均值(Mean)、標(biāo)準(zhǔn)差(SD)和變異系數(shù)(CV%),對(duì)測(cè)量性狀進(jìn)行K-S正態(tài)性檢驗(yàn)和相關(guān)分析;將殼形態(tài)性狀(殼高、殼長(zhǎng)和殼寬)作為自變量,體質(zhì)量性狀(活體總重和軟體部重)作為因變量進(jìn)行通徑分析并計(jì)算決定系數(shù),分析每個(gè)殼形態(tài)性狀分別對(duì)活體總重和軟體部重的直接作用和間接作用;通過(guò)逐步回歸分析法確立殼形態(tài)性狀估計(jì)體質(zhì)量性狀最優(yōu)回歸方程。顯著性水平設(shè)定P<0.05,極顯著水平設(shè)定P<0.01。

    活體總重YL和軟體部重YS的回歸采用的線(xiàn)性模型為:

    Y=a+b1X1+b2X2+b3X3+……biXi。

    式中:Y(YL或YS)是因變量;a是常數(shù)項(xiàng);bi是偏回歸系數(shù);Xi是自變量。根據(jù)下列公式由偏回歸系數(shù)進(jìn)一步計(jì)算通徑系數(shù)[22]。

    PXiXj=rijPj,Y(i≠j)。

    式中:PY.Xi是Xi對(duì)Y的直接通徑系數(shù);SXi是自變量Xi的標(biāo)準(zhǔn)差;SY是因變量Y的標(biāo)準(zhǔn)差;PXiYj是Xi通過(guò)Xj對(duì)Y的間接通徑系數(shù);rij是相關(guān)系數(shù);Pj,Y是Xj對(duì)Y的通徑系數(shù)。根據(jù)下列公式計(jì)算直接決定系統(tǒng)和間接決定系數(shù)[22]。

    di=PY.Xi2;

    dij=2rijPY.XiPY.Xj。

    式中:di是直接決定系數(shù);dij是間接決定系數(shù)。

    2 結(jié)果

    2.1 殼形態(tài)和體質(zhì)量性狀的參數(shù)分析

    所測(cè)性狀的相關(guān)統(tǒng)計(jì)量見(jiàn)表1,其中殼高、殼長(zhǎng)和殼寬3個(gè)殼形態(tài)性狀的變異系數(shù)相對(duì)較小,而2個(gè)體質(zhì)量性狀的變異系數(shù)相對(duì)較大,軟體部重變異系數(shù)最大,為34.034%。K-S正態(tài)性檢驗(yàn)結(jié)果表明,5個(gè)性狀原始數(shù)據(jù)均符合正態(tài)分布,可進(jìn)行通徑分析。

    表1 殼橙品系長(zhǎng)牡蠣各性狀統(tǒng)計(jì)量分析Table 1 Statistics of various traits in the orange-shell strain of the Pacific oyster Crassostrea gigas

    2.2 殼形態(tài)和體質(zhì)量性狀間的相關(guān)性分析

    5個(gè)性狀間的表型相關(guān)均呈極顯著(P<0.01),其中以活體總重與軟體部重間的相關(guān)系數(shù)最大,活體總重和軟體部重與殼形態(tài)性狀間的相關(guān)系數(shù)大小關(guān)系相同,依次為:殼高>殼長(zhǎng)>殼寬(見(jiàn)表2)。

    2.3 殼形態(tài)性狀對(duì)體質(zhì)量性狀的通徑分析和相關(guān)指數(shù)分析

    以殼高、殼長(zhǎng)、殼寬為自變量,活體總重和軟體部重為因變量進(jìn)行通徑分析,結(jié)果見(jiàn)表3。通徑分析反映自變量對(duì)因變量的直接作用,殼高對(duì)活體總重和軟體部重的通徑系數(shù)最大(0.447和0.451),殼寬次之(0.336和0.310),且均為極顯著(P<0.01)。殼形態(tài)性狀對(duì)活體總重和軟體部重的相關(guān)指數(shù)分別為0.860和0.683。

    表2 殼橙品系長(zhǎng)牡蠣各性狀間表型相關(guān)系數(shù)Table 2 Phenotypic correlation coefficients among the traits in the orange-shell strain of the Pacific oyster

    注:**表示極顯著性相關(guān)(P<0.01),下同。

    Note:**means very significantly different (P<0.01), et sequentia.

    表3 殼橙品系長(zhǎng)牡蠣殼形態(tài)性狀對(duì)體質(zhì)量性狀的通徑分析和相關(guān)指數(shù)分析Table 3 Path coefficient and correlation index of shell morphological traits on body weight traits in the orange-shell strain of the Pacific oyster

    Note:①Weight traits;②Path coefficient;③Correlation indexR2

    2.4 殼形態(tài)性狀對(duì)體質(zhì)量性狀的作用

    殼形態(tài)性狀分別對(duì)活體總重和軟體部重的相關(guān)系數(shù),及分解后的直接作用和間接作用見(jiàn)表4。3個(gè)殼形態(tài)性狀對(duì)活體總重的作用中,殼高對(duì)活體總重的直接作用最大(0.477,P<0.01),殼寬次之(0.336,P<0.01);通過(guò)分析各間接作用發(fā)現(xiàn),殼長(zhǎng)對(duì)活體總重的間接作用最大(0.441)。而3個(gè)殼形態(tài)性狀對(duì)軟體部重的作用中,直接作用最大的同樣是殼高(0.451,P<0.01),殼寬次之(0.310,P<0.01),殼長(zhǎng)對(duì)活體總重的間接作用最大(0.413)。

    2.5 殼形態(tài)性狀對(duì)體質(zhì)量性狀的決定系數(shù)

    單一殼形態(tài)性狀(殼高、殼長(zhǎng)或殼寬)對(duì)體質(zhì)量性狀(活體總重和軟體部重)的直接決定系數(shù)和兩兩殼形態(tài)性狀對(duì)體質(zhì)量性狀的共同決定系數(shù)見(jiàn)表5。殼高對(duì)活體總重和軟體部重的直接決定系數(shù)最大,分別為0.228和0.203。殼高和殼長(zhǎng)對(duì)活體總重和軟體部重的共同決定系數(shù)最大,分別為0.176和0.128。殼橙品系長(zhǎng)牡蠣殼形態(tài)性狀對(duì)體質(zhì)量性狀的決定系數(shù)的總和∑di與各自的相關(guān)指數(shù)R2相近,所得結(jié)果與通徑分析和相關(guān)指數(shù)分析結(jié)果一致。

    表4 殼橙品系長(zhǎng)牡蠣殼形態(tài)性狀對(duì)體質(zhì)量性狀的作用

    Note:①Weight traits;②Shell morphological traits;③Correlation coefficient;④Direct effect;⑤Indirect effect

    表5 殼形態(tài)性狀對(duì)體質(zhì)量性狀的決定系數(shù)

    Note:①Weight traits;②Shell morphological traits;③Correlation indexR2

    2.6 殼形態(tài)性狀對(duì)體質(zhì)量性狀的多元回歸分析

    通過(guò)逐步回歸分析法得到殼形態(tài)性狀對(duì)體質(zhì)量性狀的作用結(jié)果見(jiàn)表6。殼橙品系長(zhǎng)牡蠣殼高、殼長(zhǎng)、殼寬與活體總重和軟體部重呈極顯著相關(guān)(P<0.01),多元回歸方程為:

    YL=-88.783+0.747XH+1.399XW+0.942XL,R2=0.860;

    YS=-28.399+0.258XH+0.471XW+0.266XL,R2=0.683。

    表6 殼橙品系長(zhǎng)牡蠣殼形態(tài)性狀的偏回歸系數(shù)檢驗(yàn)

    Note:①Weight traits;②Shell morphological traits;③Partial regression coefficient;④Standard error;⑤Ttest value;⑥Pvalue;⑦Constant

    由表7可知,3個(gè)殼形態(tài)性狀的偏回歸系數(shù)對(duì)2個(gè)體質(zhì)量性狀均達(dá)到極顯著水平(P<0.01),能夠客觀的反映出殼形態(tài)性狀與體質(zhì)量性狀間的相互影響,可應(yīng)用于實(shí)際統(tǒng)計(jì)操作。

    表7 殼橙品系長(zhǎng)牡蠣殼形態(tài)性狀對(duì)體質(zhì)量性狀的多元回歸方差分析

    Note:①Weight traits;②Sources of variation;③Sum of square;④Degree of freedom;⑤Mean square;⑥Ftest value;⑦Pvalue;⑧Regression;⑨Residual;⑩Total

    3 討論

    優(yōu)良性狀親本的選擇是通過(guò)挑選滿(mǎn)足期望表型性狀的個(gè)體實(shí)現(xiàn)的,因此,在開(kāi)始育種前,通過(guò)估算表型的相關(guān)性,提高選擇的準(zhǔn)確性可有效提高性狀改良的速率[31]。本研究測(cè)量了殼橙品系長(zhǎng)牡蠣的5個(gè)主要經(jīng)濟(jì)性狀,其中體質(zhì)量性狀作為最直接的選育目標(biāo),相對(duì)于3個(gè)殼形態(tài)性狀(殼高、殼長(zhǎng)、殼寬)在測(cè)量上不夠直觀,并且3個(gè)殼形態(tài)性狀與2個(gè)體質(zhì)量性狀間存在不同程度的直接和間接的影響,因此,需要利用通徑分析法確定影響殼橙品系長(zhǎng)牡蠣體質(zhì)量性狀的主要?dú)ば螒B(tài)性狀。

    本研究中,殼高、殼長(zhǎng)、殼寬、活體總重和軟體部重5個(gè)數(shù)量性狀的變異范圍與香港牡蠣[20]和近江牡蠣[16]的研究結(jié)果相近,其中殼形態(tài)性狀(殼高、殼長(zhǎng)和殼寬)的變異系數(shù)較小,體質(zhì)量性狀(活體總重和軟體部重)的變異系數(shù)較大,這可能與長(zhǎng)牡蠣個(gè)體間肥滿(mǎn)度差異較大有關(guān),具備較大的改良潛力。長(zhǎng)牡蠣殼橙品系的殼高、殼長(zhǎng)、殼寬、活體總重和軟體部重5個(gè)數(shù)量性狀間均存在極顯著的表型相關(guān)關(guān)系,活體總重和軟體部重與殼形態(tài)性狀間的相關(guān)系數(shù)大小關(guān)系相同,依次為:殼高>殼長(zhǎng)>殼寬。在進(jìn)行通徑分析和計(jì)算決定系數(shù)時(shí),當(dāng)相關(guān)指數(shù)R2或各自變量對(duì)因變量的單獨(dú)決定系數(shù)與兩兩共同決定系數(shù)的總和∑di(在數(shù)值上R2=∑di)大于或等于0.850,方可確定影響因變量的主要自變量[7, 15]。本研究中,活體總重的R2(0.861)與∑d(0.860)近似相等且大于0.850,表明殼高、殼長(zhǎng)、殼寬是影響活體總重的主要性狀[17]。通徑分析結(jié)果中,對(duì)活體總重直接作用最大的是殼高(0.477,P<0.01),大于間接作用(0.317,P<0.01);殼寬(0.336,P<0.01)和殼長(zhǎng)(0.326,P<0.01)對(duì)活體總重的直接作用相近,且殼寬和殼長(zhǎng)的間接作用均大于兩者的直接作用,表明殼寬和殼長(zhǎng)主要是通過(guò)影響殼高從而間接影響活體總重,殼長(zhǎng)間接作用大于殼寬,而殼高是影響活體總重的主要?dú)ば螒B(tài)性狀。結(jié)合決定系數(shù)的分析結(jié)果,得出殼高為影響活體總重的主要因素,殼長(zhǎng)為次要因素。在櫛孔扇貝(Chlamysfarreri)[15]、日本鏡蛤(Dosiniajaponica)[32]、蝦夷扇貝(Patinopectenyessoensis)[19]、黑蝶貝(Pinctadamargaritifera)[33]的相關(guān)研究中,同樣得到類(lèi)似的結(jié)論。同樣,近江牡蠣對(duì)活體總重直接作用最大的是殼高,間接作用最大的是殼長(zhǎng)[16]。在對(duì)長(zhǎng)牡蠣和福建牡蠣(葡萄牙牡蠣)(Crossostreaangulata)的比較研究中,林青[29]也得出殼高對(duì)活體總重影響作用最大的結(jié)果。

    軟體部重的R2<0.850,說(shuō)明影響軟體部重的主要因素不止是本研究所測(cè)的3個(gè)殼形態(tài)性狀,還存在其他影響較大的因素,在其他研究中也存在與軟體部重的殼形態(tài)性狀R2<0.850的現(xiàn)象。在蝦夷扇貝中,軟體部重可能與貝殼的凹凸、兩殼的絞合度、貝殼表面積和水環(huán)境等因素有關(guān);在菲律賓蛤仔(Ruditapesphilippinarum)中,貝齡、性別、發(fā)育時(shí)期、活體總重、殼質(zhì)量、養(yǎng)殖方式等是軟體部重潛在的影響因素[17, 21, 23]。在不同牡蠣種類(lèi)的軟體部重通徑分析中同樣得到R2<0.850[16, 20, 29],推測(cè)殼形態(tài)性狀之外的因素對(duì)軟體部重有較大影響。但是軟體部重的R2與∑d近似相等,決定系數(shù)分析的結(jié)果和通徑分析結(jié)果相同,表明在這3個(gè)殼形態(tài)性狀中,通過(guò)選育殼高可使軟體部重獲得較好的選育效果[27]。

    4 結(jié)語(yǔ)

    本研究確定了殼橙品系長(zhǎng)牡蠣活體總重的主要影響因素為殼高,次要因素為殼長(zhǎng);殼高對(duì)殼橙品系長(zhǎng)牡蠣軟體部重有較大的影響;獲得了殼高、殼長(zhǎng)、殼寬與活體總重和軟體部重的多元回歸方程,為殼橙品系長(zhǎng)牡蠣體質(zhì)量性狀的遺傳改良提供了重要的基礎(chǔ)資料。

    [1] 農(nóng)業(yè)部漁業(yè)漁政管理局. 中國(guó)漁業(yè)統(tǒng)計(jì)年鑒[M]. 北京: 中國(guó)農(nóng)業(yè)出版社, 2016: 28-29.

    Ministry of Agriculture Fisheries. China Fishery Statistical Year Book[M]. Beijing: China Agriculture Press, 2016: 28-29.

    [2] 林志華, 董迎輝. 文蛤“萬(wàn)里紅”[J]. 中國(guó)水產(chǎn), 2015, 479(10): 50-52.

    Li Z H, Dong Y W. “Red Forever” ClamMeretrixmeretrix[J]. China Fisheries, 2015, 479(10): 50-52.

    [3] 許飛, 鄭懷平, 張海濱, 等. 海灣扇貝“中科紅”品種與普通養(yǎng)殖群體不同溫度下早期性狀的比較[J]. 水產(chǎn)學(xué)報(bào), 2008, 32(6): 876-83.

    Xu F, Zheng H P, Zhang H B, et al. Comparison on the larval traits between strain“Zhongkehong” and common cultured population of bay scallopArgopectenirradiansunder different temperature[J]. Journal of Fisheries of China, 2008, 32(6): 876-83.

    [4] 閆喜武, 霍忠明, 楊鳳, 等. 菲律賓蛤仔“斑馬蛤”[J]. 中國(guó)水產(chǎn), 2015, 478(9): 52-53.

    Yan X W, Huo ZM, Yang F, et al. “Zebra Clam” Manila ClamRuditapesphilippinarum[J]. China Fisheries, 2015, 478(9): 52-53.

    [5] 鄭懷平, 劉合露, 陳興強(qiáng), 等. 華貴櫛孔扇貝“南澳金貝”[J]. 中國(guó)水產(chǎn), 2015, 479(10): 57-58.

    Zheng H P, Liu H L, Chen X Q, et al. “Nan'ao golden Scallop”Noble ScallopChlamysnobilis[J]. China Fisheries, 2015, 479(10): 57-58.

    [6] 叢日浩, 李琪, 葛建龍, 等. 長(zhǎng)牡蠣4種殼色家系子代的表型性狀比較[J]. 中國(guó)水產(chǎn)科學(xué), 2014, 21(3): 494-502.

    Cong R H, Li Q, Ge J L, et al. Comparison of phenotypic traits of four shell color families of the Pacific oyster (Crassostreagigas)[J]. Journal of Fishery Sciences of China, 2014, 21(3): 494-502.

    [7] 常亞青, 張存善, 曹學(xué)彬, 等. 1齡蝦夷扇貝形態(tài)性狀對(duì)重量性狀的影響效果分析[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào), 2008, 23(5): 330-334.

    Chang Y Q, Zhang C S, Cao X B, et al. Effect of morphometrical traits on weght traits in one-year old yesso scallopPatinopectenyessoensis[J]. Journal of Dalian Ocean University, 2008, 23(5): 330-334.

    [8] 劉峰, 陳琳, 樓寶, 等. 小黃魚(yú)(Pseudosciaenapolyactis)形態(tài)性狀與體質(zhì)量的相關(guān)性及通徑分析[J]. 海洋與湖沼, 2016, 47(3): 655-662.

    Liu F, Chen L, Lou B, et al. Correlation and path coefficient analysis on body weight and morphometric traits of small yellow croakerPseudosciaenapolyactis[J]. Oceanologia et Limnologia Sinica, 2016, 47(3): 655-662.

    [9] 劉賢德, 蔡明夷, 王志勇, 等. 閩-粵東族大黃魚(yú)生長(zhǎng)性狀的相關(guān)與通徑分析[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2008, 38(6): 916-920.

    Liu X D, Cai M Y, Wang Z Y, et al. The Correlation and Path Analysis for Growth-Related Traits of Large Yellow CroakerPseudosciaenacroceafrom Min-Yuedong Tribe[J]. Periodical of Ocean University of China, 2008, 38(6): 916-920.

    [10] 劉小林, 吳長(zhǎng)功, 張志懷, 等. 凡納對(duì)蝦形態(tài)性狀對(duì)體重的影響效果分析[J]. 生態(tài)學(xué)報(bào), 2004, 24(4): 857-862.

    Liu X L, Wu C G, Zhang Z H, et al. Mathematical analysis of effects of morphometric attributes on body weight forPenaeusvannamei[J]. Acta Ecologica Sinica, 2004, 24(4): 857-862.

    [11] Ma H, MA C, Ma L, et al. Correlation of growth-related traits and their effects on body weight of the mud crab (Scyllaparamamosain)[J]. Genet Mol Res, 2013, 12(4): 4127-4136.

    [12] Deng Y W, Du X D, Wang Q H, et al. Correlation and path analysis for growth traits in F1 population of pearl oyster Pinctada martensii[J]. Marine Science Bulletin, 2008, 10(2): 68-73.

    [13] Luo X, Ke C, You W. Estimates of correlations for shell morphological traits on body weight of interspecific hybrid abalone (HaliotisdiscushannaiandHaliotisgigantea)[J]. Journal of Shellfish Research, 2013, 32(1): 115-118.

    [14] Zhao L, He Y, Yang F, et al. Correlation and path analysis of morphological and weight traits in marine gastropod Glossaulax reiniana[J]. Chinese Journal of Oceanology and Limnology, 2014, 32(4): 821-827.

    [15] 劉小林, 常亞青, 相建海, 等. 櫛孔扇貝殼尺寸性狀對(duì)活體重的影響效果分析[J]. 海洋與湖沼, 2002, 33(6): 673-678.

    Liu X L, Chang Y Q, Xiang J H, et al. Analysis of effects of shell size characters on live weight in Chines ScallopChlamysFarreri[J]. Oceanologia et Limnologia Sinica, 2002, 33(6): 673-678.

    [16] 孫澤偉, 鄭懷平, 楊彥鴻, 等. 近江牡蠣養(yǎng)殖群體數(shù)量性狀間的相關(guān)及通徑分析[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2010, 26(6): 332-336.

    Sun Z W, Zheng H P, Yang Y H, et al. Correlation and path analysis to quantitative traits for a cultured population of Jinjiang oysterCrassostreahongkongensis[J]. Chinese Agricultural Science Bulletin, 2010, 26(6): 332-336.

    [17] 閆喜武, 王琰, 郭文學(xué), 等. 四角蛤蜊形態(tài)性狀對(duì)重量性狀的影響效果分析[J]. 水產(chǎn)學(xué)報(bào), 2011, 35(10): 1513-1518.

    Yan X W, Wang Y, Guo W X, et al. Effects of shell morphological traits on the weight traits of clamMactraveneriformisReeve along northern coast in China[J]. Journal of Fisheries of China, 2011, 35(10): 1513-1518.

    [18] 王沖, 孫同秋, 王玉清, 等. 不同群體毛蚶形態(tài)性狀對(duì)重量性狀的影響效果分析[J]. 海洋漁業(yè), 2015, 37(5): 427-433.

    Wang C, Sun T Q, Wang Y Q, et al. Effects of shell morphological traits on the weight traits ofScapharcasubcrenataLischke in different populations[J]. Marine Fisheries, 2015, 37(5): 427-433.

    [19] 杜美榮, 劉毅, 蔣增杰, 等. 底播蝦夷扇貝數(shù)量性狀的相關(guān)性和通徑分析[J]. 水產(chǎn)科學(xué), 2015, 34(1): 8-13.

    Du M R, Liu Y, Jiang Z J, et al. Correlation and path analysis of quantitative traits in bottom-cultured Yesso scallopPatinopectenyessoensis[J]. Fisheries Science, 2015, 34(1): 8-13.

    [20] 肖述, 符政君, 喻子牛. 香港巨牡蠣雌雄群體的數(shù)量性狀通徑分析[J]. 南方水產(chǎn)科學(xué), 2011, 7(4): 1-9.

    Xiao S, Fu Z J, Yu Z N. Path analysis of quantitative traits of male and female Hong KongoysterCrassostreahongkongensis[J]. South China Fisheries Science, 2011, 7(4): 1-9.

    [21] HUO Zhong-Ming, YAN Xi-Wu, ZHAO Li-Qiang, et al. Effects of shell morphological traits on the weight traits of Manila clam (Ruditapesphilippinarum)[J]. Acta Ecologica Sinica, 2010, 30(5): 251-256.

    [22] 張根芳, 張文府, 方愛(ài)萍, 等. 養(yǎng)殖條件下不同年齡背瘤麗蚌(Lamprotulaleai)數(shù)量性狀的相關(guān)與通徑分析[J]. 海洋與湖沼, 2014, 45(5): 1115-1121.

    Zhang G F, Zhang W F, Fang A P, et al. Correlation and path analysis of quantitaitve traits of different-ageLamprotulaleaiin artificial breeding[J]. Oceanologia et Limnologia Sinica, 2014, 45(5): 1115-1121.

    [23] 劉輝, 張興志, 鹿瑤, 等. 菲律賓蛤仔橙色品系殼形態(tài)性狀對(duì)質(zhì)量性狀的通徑及多元回歸分析[J]. 大連海洋大學(xué)學(xué)報(bào), 2015, 30(5): 514-518.

    Liu H, Zhang X Z, Lu Y, et al. Effects of shell morphological traits on the weight traits of orange strain of Manila clamRuditapesphilippinarum[J]. Journal of Dalian Ocean University, 2015, 30(5): 514-518.

    [24] 孫秀俊, 楊?lèi)?ài)國(guó), 劉志鴻, 等. 2種殼色蝦夷扇貝的形態(tài)學(xué)指標(biāo)比較分析[J]. 安徽農(nóng)業(yè)科學(xué), 2008, 36(23): 10008-10160.

    Sun X J, Yang A G, Liu Z H, et al. Comparative analysis of morphological indices of Papanese scallops with 2 shell colors[J]. Journal of Anhui Agricultural Sciences, 2008, 36(23): 10008-10160.

    [25] 趙鵬, 丁君, 常亞青. 兩種殼色蝦夷扇貝殼體尺性狀對(duì)活體重影響效果的分析[J]. 大連海洋大學(xué)學(xué)報(bào), 2011, 26(1): 1-5.

    Zhao P, Ding J, Chang Y Q. Effect of shell size characters on live body weight in Japanese scallopMizuhopectenyessoensiswith two colors[J]. Journal of Dalian Ocean University, 2011, 26(1): 1-5.

    [26] 王慶志, 李琪, 劉士凱, 等. 長(zhǎng)牡蠣不同地理群體選育系數(shù)量性狀的比較[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 41(Z2): 36-41.

    Wang Q Z, Li Q, Liu S K, et al. Comparison of quantitative traits among the breeding lines of different geographic populations ofCrassostreagigas[J]. Periodical of Ocean University of China, 2011, 41(Z2): 36-41.

    [27] 王慶志, 李琪, 劉世凱, 等. 長(zhǎng)牡蠣成體生長(zhǎng)性狀的遺傳參數(shù)估計(jì)[J]. 中國(guó)水產(chǎn)科學(xué), 2012, 19(4): 700-706.

    Wang Q Z, Li Q, Liu S K, et al. Estimates of genetic parameters for growth-related traits in adultCrassostreagigas[J]. Journal of Fishery Sciences of China, 2012, 19(4): 700-706.

    [28] 孔寧, 李琪, 叢日浩, 等. 長(zhǎng)牡蠣F3代快速生長(zhǎng)選育群體生長(zhǎng)特性的研究[J]. 海洋科學(xué), 2015, 39(3): 7-11.

    Kong N, Li Q, Cong R H, et al. Study on growth characteristics of the selected third generationof Pacific oysterCrassostreagigaswith rapid growth[J]. Marine Sciences, 2015, 39(3): 7-11.

    [29] 林清, 王亞駿, 王迪文, 等. 太平洋牡蠣和葡萄牙牡蠣養(yǎng)殖群體數(shù)量性狀比較分析[J]. 海洋通報(bào), 2014, 33(1): 106-111.

    Lin Q, Wang Y J, Wang D W, et al. Comparison analysis of quantitative traits and path between twocultured populations ofCrassostreagigasandC.angulata[J]. Marine Science Bulletin, 2014, 33(1): 106-111.

    [30] 杜家菊, 陳志偉. 使用SPSS線(xiàn)性回歸實(shí)現(xiàn)通徑分析的方法[J]. 生物學(xué)通報(bào), 2010, 45(2): 4-6.

    Du J J, Chen Z W. Method of path analysis using SPSS regression[J]. Bulletin of Biology, 2010, 45(2): 4-6.

    [31] 張國(guó)范, 李霞, 薛真福. 我國(guó)養(yǎng)殖貝類(lèi)大規(guī)模死亡的原因分析及防治對(duì)策[J]. 中國(guó)水產(chǎn), 1999(9): 34-39.

    Zhang G F, Li X, Xue Z F. Cause analysis and prevention methods of mass cutured shellfish death in China[J]. China Fisheries, 1999(9): 34-39.

    [32] 張偉杰, 常亞青, 丁君, 等. 日本鏡蛤(DosiniajaponicaReeve)殼尺寸與重量性狀的相關(guān)與回歸分析[J]. 海洋與湖沼, 2013, 44(3): 796-800.

    Zhang W J, Chang Y Q, Ding J, et al. Correlation and regression of shell size and weight ofDosinlajaponicaReeve[J]. Oceanologia et Limnologia Sinica, 2013, 44(3): 796-800.

    [33] 嚴(yán)俊賢, 劉寶鎖, 李有寧, 等. 野生黑蝶貝表型性狀對(duì)體質(zhì)量的影響分析[J]. 水產(chǎn)科學(xué), 2015, 34(9): 560-564.

    Yan J X, Liu B S, Li Y N, et al. Effects of phenotypic traits on body weight on wild population of Pearl oysterPinctadamargaritifera[J]. Fisheries Science, 2015, 34(9): 560-564.

    MultipleRegressionandPathAnalysisofMorphologicalandWeightTraitsofPacificOyster(Crassostreagigas)Orange-ShellStrain

    HAN Zi-Qiang, LI Qi

    (The Key Laboratory of Mariculture (Ocean University of China), Ministry of Education, Qingdao 266003, China)

    To evaluate the effect of shell morphological traits on weight traits in the orange-shell strain of Pacific oysterCrassostreagigas, shell height (XH), shell length (XL), shell width (XW), live body weight (YL) and soft tissue weight (YS) were measured for 144 individuals of two-year-old orange-shell strain of Pacific oyster. Morphological traits (XH,XLandXW) were used as independent variables and weight traits (YL,YS) as dependent variables for path coefficient and multiple regression analyses. Correlation indices showed that the direct effects of correlation coefficient between each shell morphological traits (XH,XLandXW) and weight traits (YLandYS) were all significantly different (P<0.01). In path analysis, the direct effects of the shell morphological traits on live body weight were in order of shell height (0.477) > shell length (0.336) > shell width (0.326). Shell height had the maximum direct effect (0.228) on the live body weight, and was the key influencing factor. However, the correlation indices (R2) of morphological traits against the edible tissue weight was less than 0.85, indicating that some other factors might associate with the edible tissue weight. The optimum multiple regression equations were also obtained through stepwise multiple regression analysis.

    Crassostreagigas; orange-shell strain; morphological trait; weight trait; correlation analysis; path analysis; multiple regression analysis

    S968.31

    A

    1672-5174(2017)12-046-07

    責(zé)任編輯 朱寶象

    10.16441/j.cnki.hdxb. 20170001

    韓自強(qiáng), 李琪. 長(zhǎng)牡蠣殼橙品系形態(tài)性狀與體質(zhì)量的相關(guān)及通徑分析[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 47(12): 46-52.

    HAN Zi-Qiang, LI Qi. Multiple regression and path analysis of morphological and weight traits of pacific oyster (Crassostreagigas) orange-shell strain [J]. Periodical of Ocean University of China, 2017, 47(12): 46-52.

    泰山學(xué)者種業(yè)計(jì)劃專(zhuān)家項(xiàng)目;農(nóng)業(yè)科技成果轉(zhuǎn)化資金項(xiàng)目(2014GB2B020029);山東省科技發(fā)展計(jì)劃項(xiàng)目(2014GHY115002);鰲山科技創(chuàng)新計(jì)劃項(xiàng)目(2015ASKJ02)資助

    Supported by Taishan Scholar Seed Industry Experts Project;Agriculture Science Technology Achievement Transformation Fund (2014GB2B020029);Science and Technology Development Project of Shandong Province (2014GHY115002);Aoshan Science and Technology Innovation Program (2015ASKJ02)

    2017-01-01;

    2017-03-29

    韓自強(qiáng)(1993-),男,碩士生。E-mail: hanziqiang1993@gmail.com

    ** 通訊作者:E-mail: qili66@ouc.edu.cn

    猜你喜歡
    通徑活體牡蠣
    牡丹江市氣象因子與PM2.5濃度影響通徑分析
    黑龍江氣象(2021年2期)2021-11-05 07:06:58
    告別自汗用牡蠣,四季都輕松
    新型全自動(dòng)氣動(dòng)通徑機(jī)設(shè)計(jì)要點(diǎn)探究
    張帆:肝穿刺活體組織學(xué)檢查
    肝博士(2021年1期)2021-03-29 02:32:08
    讓活體肝移植研究走上世界前沿
    氣動(dòng)通徑機(jī)設(shè)計(jì)及使用中應(yīng)注意的幾個(gè)問(wèn)題
    活體器官移植,你怎么看?
    “汪星人”要打“飛的”——話(huà)說(shuō)活體空運(yùn)
    曇石山文化的牡蠣器
    大眾考古(2015年6期)2015-06-26 08:27:16
    新型鋼管氣動(dòng)通徑機(jī)的研發(fā)
    一区二区三区四区激情视频 | 午夜福利高清视频| 两个人视频免费观看高清| 不卡视频在线观看欧美| 亚洲无线在线观看| 欧美最黄视频在线播放免费| 国产亚洲91精品色在线| 国产精品永久免费网站| 99久久精品热视频| 在线播放国产精品三级| 日韩欧美一区二区三区在线观看| 久久精品国产亚洲av天美| 少妇丰满av| 久久精品国产亚洲av天美| 国产精品久久久久久精品电影| av黄色大香蕉| 久久中文看片网| 久久久久久久久久久丰满 | 国产黄a三级三级三级人| 九色国产91popny在线| 欧美黑人欧美精品刺激| 免费观看精品视频网站| 九色成人免费人妻av| 成人精品一区二区免费| 日韩强制内射视频| 亚洲,欧美,日韩| 亚洲av熟女| 三级男女做爰猛烈吃奶摸视频| 免费搜索国产男女视频| 天堂网av新在线| www日本黄色视频网| 亚洲avbb在线观看| 草草在线视频免费看| 3wmmmm亚洲av在线观看| 搡女人真爽免费视频火全软件 | 久久久久久久午夜电影| 免费av观看视频| 校园人妻丝袜中文字幕| 老熟妇乱子伦视频在线观看| 久久午夜福利片| 两个人的视频大全免费| 非洲黑人性xxxx精品又粗又长| 久久国产乱子免费精品| 欧美高清性xxxxhd video| 97超视频在线观看视频| 国产色婷婷99| or卡值多少钱| 嫩草影院精品99| 他把我摸到了高潮在线观看| av.在线天堂| 久久久久久久久久成人| 久久国产精品人妻蜜桃| 啦啦啦观看免费观看视频高清| or卡值多少钱| 99视频精品全部免费 在线| 国产成人av教育| 最新中文字幕久久久久| 精品久久久噜噜| 欧美一区二区亚洲| 国产精品亚洲一级av第二区| 国产成人aa在线观看| 精品人妻视频免费看| 色5月婷婷丁香| 人人妻人人看人人澡| 又爽又黄无遮挡网站| 午夜老司机福利剧场| 舔av片在线| 国产精品一及| 美女大奶头视频| 午夜日韩欧美国产| 俺也久久电影网| 精品人妻视频免费看| 人妻丰满熟妇av一区二区三区| 超碰av人人做人人爽久久| 亚洲人成网站在线播放欧美日韩| 91久久精品国产一区二区三区| 国产高清三级在线| 亚洲精品色激情综合| 色5月婷婷丁香| 午夜精品在线福利| 精品人妻视频免费看| 午夜精品一区二区三区免费看| 国产69精品久久久久777片| 免费av毛片视频| 在线播放无遮挡| 性色avwww在线观看| 久久人人精品亚洲av| 精品人妻视频免费看| 国产激情偷乱视频一区二区| 夜夜爽天天搞| 婷婷精品国产亚洲av在线| 乱码一卡2卡4卡精品| 欧美日韩乱码在线| 麻豆久久精品国产亚洲av| 如何舔出高潮| 色哟哟·www| 国产老妇女一区| 欧美中文日本在线观看视频| 亚洲自偷自拍三级| 国产精品一区二区三区四区免费观看 | 久久精品国产鲁丝片午夜精品 | 免费黄网站久久成人精品| 久久久色成人| 午夜福利在线在线| 久久久久久久亚洲中文字幕| 亚洲欧美日韩卡通动漫| 国产精品综合久久久久久久免费| 亚洲久久久久久中文字幕| or卡值多少钱| 人妻夜夜爽99麻豆av| 别揉我奶头 嗯啊视频| 欧美性感艳星| 日韩欧美 国产精品| 国产aⅴ精品一区二区三区波| 午夜亚洲福利在线播放| 亚洲成人久久性| 欧洲精品卡2卡3卡4卡5卡区| 偷拍熟女少妇极品色| 午夜福利高清视频| 亚洲av成人精品一区久久| 嫩草影院入口| 亚洲国产精品sss在线观看| a级一级毛片免费在线观看| 午夜福利视频1000在线观看| 久久精品国产自在天天线| 欧美高清成人免费视频www| 一本一本综合久久| 丰满的人妻完整版| 国产探花极品一区二区| 精品无人区乱码1区二区| 日本精品一区二区三区蜜桃| 婷婷精品国产亚洲av| 老司机午夜福利在线观看视频| 久久久久久久午夜电影| 色哟哟哟哟哟哟| 九色国产91popny在线| 人妻久久中文字幕网| 欧美bdsm另类| 亚洲国产欧洲综合997久久,| 亚洲欧美日韩无卡精品| 亚洲av五月六月丁香网| 在线播放国产精品三级| 中文字幕人妻熟人妻熟丝袜美| 国内精品久久久久精免费| 国产欧美日韩精品一区二区| 日本五十路高清| 一个人看的www免费观看视频| 99久久九九国产精品国产免费| av福利片在线观看| 日韩欧美 国产精品| 神马国产精品三级电影在线观看| 婷婷色综合大香蕉| 日日啪夜夜撸| 国产午夜精品久久久久久一区二区三区 | 有码 亚洲区| 搡老妇女老女人老熟妇| 午夜福利在线观看吧| 欧美色欧美亚洲另类二区| 免费看日本二区| 中文字幕久久专区| 国产亚洲精品综合一区在线观看| 18+在线观看网站| 麻豆成人av在线观看| 99国产极品粉嫩在线观看| 在线免费观看的www视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av第一区精品v没综合| 禁无遮挡网站| 国产单亲对白刺激| 天堂网av新在线| 国产黄片美女视频| 永久网站在线| 婷婷色综合大香蕉| 亚洲av成人av| 国模一区二区三区四区视频| 亚洲成人久久爱视频| 亚洲精品乱码久久久v下载方式| 天堂网av新在线| 婷婷精品国产亚洲av| 国产爱豆传媒在线观看| 亚洲精品成人久久久久久| 亚洲av免费在线观看| 一区二区三区四区激情视频 | 国产日本99.免费观看| 在现免费观看毛片| 少妇人妻精品综合一区二区 | 亚洲午夜理论影院| 国产精品亚洲一级av第二区| 身体一侧抽搐| 人人妻人人澡欧美一区二区| 97超视频在线观看视频| 狠狠狠狠99中文字幕| 久9热在线精品视频| 黄色丝袜av网址大全| 亚洲国产精品合色在线| ponron亚洲| 丰满人妻一区二区三区视频av| 淫秽高清视频在线观看| 99久久精品国产国产毛片| 夜夜夜夜夜久久久久| 成人欧美大片| 国语自产精品视频在线第100页| 男人舔奶头视频| 3wmmmm亚洲av在线观看| a级毛片a级免费在线| 免费搜索国产男女视频| 麻豆精品久久久久久蜜桃| 中文在线观看免费www的网站| 久久久久久久久久成人| 亚洲中文字幕一区二区三区有码在线看| 亚洲国产精品sss在线观看| 亚洲精品国产成人久久av| 97热精品久久久久久| 婷婷丁香在线五月| 最近在线观看免费完整版| 日本撒尿小便嘘嘘汇集6| 蜜桃久久精品国产亚洲av| 午夜日韩欧美国产| 精品不卡国产一区二区三区| 久久欧美精品欧美久久欧美| 国产精品国产三级国产av玫瑰| av视频在线观看入口| 性欧美人与动物交配| av专区在线播放| 亚洲一区二区三区色噜噜| 日韩精品有码人妻一区| 国内毛片毛片毛片毛片毛片| 一边摸一边抽搐一进一小说| 欧美3d第一页| 人妻丰满熟妇av一区二区三区| 亚洲专区国产一区二区| 日本 av在线| 男女之事视频高清在线观看| 久久久久九九精品影院| 亚洲av电影不卡..在线观看| 国产 一区 欧美 日韩| 99视频精品全部免费 在线| 日韩欧美精品v在线| 亚洲最大成人av| 久久久国产成人免费| 日本在线视频免费播放| av在线亚洲专区| 国产精品久久久久久亚洲av鲁大| 99国产极品粉嫩在线观看| 日韩中字成人| 欧美一区二区亚洲| 两人在一起打扑克的视频| 成人欧美大片| 人人妻人人澡欧美一区二区| 亚洲欧美清纯卡通| av福利片在线观看| 91久久精品国产一区二区成人| 亚洲av中文字字幕乱码综合| 男人舔女人下体高潮全视频| 少妇猛男粗大的猛烈进出视频 | 精品人妻熟女av久视频| 成熟少妇高潮喷水视频| 波多野结衣高清作品| 久久6这里有精品| 老司机午夜福利在线观看视频| 精品久久久噜噜| 国产精品国产高清国产av| 亚洲精品一区av在线观看| 国产精品自产拍在线观看55亚洲| 1024手机看黄色片| 精品国内亚洲2022精品成人| 国产高清不卡午夜福利| avwww免费| 美女xxoo啪啪120秒动态图| 能在线免费观看的黄片| 国产aⅴ精品一区二区三区波| 国产女主播在线喷水免费视频网站 | 国产伦精品一区二区三区四那| 亚洲电影在线观看av| 国产一区二区三区av在线 | 人妻制服诱惑在线中文字幕| 午夜影院日韩av| 国模一区二区三区四区视频| 免费在线观看日本一区| 久久热精品热| 99九九线精品视频在线观看视频| 日本黄色片子视频| 窝窝影院91人妻| 99久久久亚洲精品蜜臀av| 国产欧美日韩精品亚洲av| 真人一进一出gif抽搐免费| 日韩欧美在线二视频| .国产精品久久| 久久精品久久久久久噜噜老黄 | 美女xxoo啪啪120秒动态图| 三级毛片av免费| 真人做人爱边吃奶动态| 床上黄色一级片| 久久国产精品人妻蜜桃| 热99在线观看视频| 色视频www国产| 国产高潮美女av| 18禁裸乳无遮挡免费网站照片| 免费看光身美女| 欧美一区二区国产精品久久精品| 久久人妻av系列| 国产一区二区在线观看日韩| 极品教师在线视频| 日韩欧美精品v在线| 国产成人av教育| 亚洲第一区二区三区不卡| 欧美zozozo另类| 国产精品不卡视频一区二区| 亚洲熟妇熟女久久| 亚洲第一电影网av| 国产精品av视频在线免费观看| 日本黄色片子视频| 国产午夜精品久久久久久一区二区三区 | 女人十人毛片免费观看3o分钟| 韩国av一区二区三区四区| 99在线人妻在线中文字幕| 天堂影院成人在线观看| 高清日韩中文字幕在线| 尤物成人国产欧美一区二区三区| 搡女人真爽免费视频火全软件 | 日本一本二区三区精品| 深夜精品福利| 国产色婷婷99| 亚洲国产欧洲综合997久久,| 免费在线观看日本一区| 亚洲精品粉嫩美女一区| 波多野结衣巨乳人妻| 久久99热这里只有精品18| 国产欧美日韩一区二区精品| 在线播放无遮挡| 最新中文字幕久久久久| 午夜精品久久久久久毛片777| 国产午夜精品论理片| 亚洲精品日韩av片在线观看| 男人的好看免费观看在线视频| 亚洲国产精品久久男人天堂| 国产午夜精品论理片| 神马国产精品三级电影在线观看| 在线观看美女被高潮喷水网站| 两人在一起打扑克的视频| 亚洲av成人精品一区久久| 国产免费一级a男人的天堂| 亚洲av日韩精品久久久久久密| 97超级碰碰碰精品色视频在线观看| .国产精品久久| 色播亚洲综合网| 人妻制服诱惑在线中文字幕| 免费人成视频x8x8入口观看| 国语自产精品视频在线第100页| 亚洲aⅴ乱码一区二区在线播放| 亚洲,欧美,日韩| 超碰av人人做人人爽久久| av在线亚洲专区| 99九九线精品视频在线观看视频| 欧美激情国产日韩精品一区| 亚洲性夜色夜夜综合| 看黄色毛片网站| 国产精品日韩av在线免费观看| 可以在线观看毛片的网站| 亚洲内射少妇av| 亚洲三级黄色毛片| 国产成人福利小说| 日韩欧美在线乱码| 国产精品女同一区二区软件 | 男女视频在线观看网站免费| av福利片在线观看| 美女cb高潮喷水在线观看| 国产熟女欧美一区二区| 亚洲欧美日韩高清专用| 国产乱人视频| 日本黄色片子视频| 精品久久久久久久久亚洲 | 国产精品国产三级国产av玫瑰| 国产又黄又爽又无遮挡在线| 五月玫瑰六月丁香| 精品久久久久久久人妻蜜臀av| 欧美日韩中文字幕国产精品一区二区三区| 欧美日本视频| 少妇裸体淫交视频免费看高清| 亚洲,欧美,日韩| 国产黄色小视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产高清有码在线观看视频| 波野结衣二区三区在线| 99久久精品国产国产毛片| 久久精品91蜜桃| 国产亚洲av嫩草精品影院| 成年女人看的毛片在线观看| 1024手机看黄色片| 欧美又色又爽又黄视频| 色视频www国产| 亚洲欧美日韩无卡精品| 蜜桃久久精品国产亚洲av| 日韩欧美一区二区三区在线观看| 亚洲av成人精品一区久久| 国产欧美日韩精品亚洲av| 成人欧美大片| 欧美一区二区亚洲| 成年版毛片免费区| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在 | 99热网站在线观看| 日韩大尺度精品在线看网址| 男女边吃奶边做爰视频| 91在线观看av| 精品免费久久久久久久清纯| 欧美不卡视频在线免费观看| 桃红色精品国产亚洲av| 日本免费一区二区三区高清不卡| aaaaa片日本免费| 亚洲乱码一区二区免费版| 男人狂女人下面高潮的视频| 中文在线观看免费www的网站| 久久精品国产亚洲av香蕉五月| 精品人妻一区二区三区麻豆 | 久久精品夜夜夜夜夜久久蜜豆| 99久久精品国产国产毛片| 91狼人影院| 亚洲av成人av| 成人二区视频| 久久精品91蜜桃| 丰满乱子伦码专区| 麻豆av噜噜一区二区三区| 国产精品一区二区三区四区免费观看 | 啦啦啦观看免费观看视频高清| 91午夜精品亚洲一区二区三区 | 伦精品一区二区三区| 99在线人妻在线中文字幕| 欧美国产日韩亚洲一区| 日韩精品青青久久久久久| 中亚洲国语对白在线视频| 中文资源天堂在线| 国产精品爽爽va在线观看网站| 国产伦一二天堂av在线观看| 一个人观看的视频www高清免费观看| 麻豆成人av在线观看| 精品欧美国产一区二区三| 亚洲精品亚洲一区二区| 精品午夜福利视频在线观看一区| 免费av毛片视频| 淫妇啪啪啪对白视频| 国产成年人精品一区二区| 成人一区二区视频在线观看| 我要搜黄色片| 小说图片视频综合网站| 亚洲成人久久爱视频| 午夜福利18| 免费黄网站久久成人精品| 日日撸夜夜添| 国产精品野战在线观看| 欧美一区二区亚洲| 在线观看66精品国产| 免费观看的影片在线观看| 日韩人妻高清精品专区| 亚洲人成伊人成综合网2020| 91av网一区二区| 日本三级黄在线观看| 国产高清视频在线观看网站| 中国美白少妇内射xxxbb| 亚洲国产精品成人综合色| 亚洲国产高清在线一区二区三| 精品一区二区三区视频在线| 国产精品福利在线免费观看| 国内久久婷婷六月综合欲色啪| 国产精品一区二区三区四区久久| 中文字幕高清在线视频| 日韩av在线大香蕉| 国产一级毛片七仙女欲春2| 精品无人区乱码1区二区| 欧美日韩精品成人综合77777| 中文在线观看免费www的网站| 床上黄色一级片| 国产精品99久久久久久久久| 久久天躁狠狠躁夜夜2o2o| 免费在线观看影片大全网站| 成人鲁丝片一二三区免费| 99久久精品热视频| 99热精品在线国产| 国内毛片毛片毛片毛片毛片| 在线观看免费视频日本深夜| 三级毛片av免费| 午夜福利欧美成人| 人妻丰满熟妇av一区二区三区| 黄色视频,在线免费观看| 国产白丝娇喘喷水9色精品| 真人一进一出gif抽搐免费| 最好的美女福利视频网| av黄色大香蕉| 丰满乱子伦码专区| 精品福利观看| 能在线免费观看的黄片| 少妇被粗大猛烈的视频| 国产成人a区在线观看| 乱人视频在线观看| 欧美日韩综合久久久久久 | 国产一级毛片七仙女欲春2| 国产日本99.免费观看| 草草在线视频免费看| 啦啦啦观看免费观看视频高清| 深爱激情五月婷婷| 中文亚洲av片在线观看爽| 午夜激情福利司机影院| 国产伦精品一区二区三区视频9| 啦啦啦韩国在线观看视频| 国国产精品蜜臀av免费| 窝窝影院91人妻| 床上黄色一级片| 精品一区二区三区av网在线观看| 久久久久久伊人网av| 国产精品久久电影中文字幕| 精品久久久久久,| 两人在一起打扑克的视频| ponron亚洲| 久9热在线精品视频| 国产午夜福利久久久久久| 少妇裸体淫交视频免费看高清| 国产精品野战在线观看| 我的女老师完整版在线观看| 色av中文字幕| 日本五十路高清| 午夜激情欧美在线| 日韩精品青青久久久久久| 成人鲁丝片一二三区免费| 丰满人妻一区二区三区视频av| 日韩在线高清观看一区二区三区 | 18+在线观看网站| 久久天躁狠狠躁夜夜2o2o| 成人国产一区最新在线观看| 日韩欧美在线二视频| 欧美一区二区国产精品久久精品| 黄色日韩在线| 免费在线观看成人毛片| 99热网站在线观看| 久久久久免费精品人妻一区二区| 国模一区二区三区四区视频| av天堂在线播放| 国产综合懂色| 少妇裸体淫交视频免费看高清| 成年女人看的毛片在线观看| 久久草成人影院| 国产欧美日韩精品亚洲av| 欧美日韩瑟瑟在线播放| 精品人妻熟女av久视频| 国产精品av视频在线免费观看| 精品久久久久久久人妻蜜臀av| 日本 欧美在线| 欧美又色又爽又黄视频| 国产综合懂色| 国产伦人伦偷精品视频| 两个人的视频大全免费| 亚洲人成网站在线播| 琪琪午夜伦伦电影理论片6080| 一个人看的www免费观看视频| 少妇人妻精品综合一区二区 | 永久网站在线| 老司机福利观看| 男女边吃奶边做爰视频| 变态另类成人亚洲欧美熟女| xxxwww97欧美| 国产高清视频在线观看网站| 黄色配什么色好看| 可以在线观看的亚洲视频| 黄色女人牲交| 91狼人影院| 真人一进一出gif抽搐免费| 少妇丰满av| 99久久精品一区二区三区| 中亚洲国语对白在线视频| 亚洲精品一区av在线观看| 国产成人a区在线观看| 亚洲av熟女| 亚洲精品亚洲一区二区| 老女人水多毛片| 69人妻影院| 一区二区三区高清视频在线| 久久精品国产亚洲av天美| 舔av片在线| 中文亚洲av片在线观看爽| 蜜桃亚洲精品一区二区三区| 国产一区二区在线av高清观看| 欧美最新免费一区二区三区| 麻豆av噜噜一区二区三区| 免费高清视频大片| 成年版毛片免费区| 国产视频内射| 久久人人爽人人爽人人片va| 国产精品98久久久久久宅男小说| 夜夜夜夜夜久久久久| 长腿黑丝高跟| 看十八女毛片水多多多| 亚洲成a人片在线一区二区| 十八禁网站免费在线| 亚洲经典国产精华液单| 日本黄色片子视频| 亚洲熟妇熟女久久| 亚洲专区国产一区二区| 国产单亲对白刺激| 色综合站精品国产| 国产视频一区二区在线看| 免费看av在线观看网站| 精品久久国产蜜桃| 国产亚洲精品久久久久久毛片| 亚洲欧美精品综合久久99| a级毛片a级免费在线| 亚洲中文字幕一区二区三区有码在线看| 午夜福利成人在线免费观看| 我要看日韩黄色一级片| 欧美日韩黄片免| 国产精品永久免费网站| 欧美一区二区亚洲| 91久久精品电影网| av黄色大香蕉| av视频在线观看入口| 精品人妻1区二区| 老女人水多毛片|