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

    如皋黃雞開產(chǎn)性狀加性遺傳及母體遺傳效應(yīng)分析

    2016-03-18 11:29:48沈曼曼竇套存胡玉萍王克華
    華南農(nóng)業(yè)大學學報 2016年2期
    關(guān)鍵詞:遺傳力動物模型

    郭 軍, 曲 亮, 沈曼曼, 竇套存, 胡玉萍, 王克華

    (中國農(nóng)業(yè)科學院 家禽研究所, 江蘇 揚州 225125)

    ?

    如皋黃雞開產(chǎn)性狀加性遺傳及母體遺傳效應(yīng)分析

    郭軍, 曲亮, 沈曼曼, 竇套存, 胡玉萍, 王克華

    (中國農(nóng)業(yè)科學院 家禽研究所, 江蘇 揚州 225125)

    摘要:【目的】以動物模型分析加性遺傳、母體遺傳對如皋黃雞開產(chǎn)性狀的影響,并計算開產(chǎn)性狀遺傳參數(shù)?!痉椒ā渴占绺撄S雞2009—2012年生產(chǎn)數(shù)據(jù),包括10 454個動物個體、31 362條生產(chǎn)記錄。依據(jù)是否考量母體效應(yīng),設(shè)計4種統(tǒng)計模型。利用WOMBAT軟件以平均信息約束最大似然法(AIREML)剖分表型方差,以似然比統(tǒng)計值檢驗模型適用性。以單性狀模型分析遺傳力,以多性狀模型分析遺傳相關(guān)系數(shù)和表型相關(guān)系數(shù)?!窘Y(jié)果】母體環(huán)境對開產(chǎn)日齡有顯著的影響,模型II為開產(chǎn)日齡最適模型,開產(chǎn)日齡遺傳力為0.28±0.02;母體遺傳對平均開產(chǎn)蛋質(zhì)量影響顯著,模型III適合計算開產(chǎn)蛋質(zhì)量遺傳力,開產(chǎn)蛋質(zhì)量加性遺傳力為0.23±0.02、母體遺傳力為0.03±0.01;母體遺傳和母體環(huán)境對開產(chǎn)體質(zhì)量有顯著影響,開產(chǎn)體質(zhì)量加性遺傳力為0.34±0.03、母體遺傳力為0.04±0.01,模型IV為開產(chǎn)體質(zhì)量最適模型。如皋黃雞開產(chǎn)蛋質(zhì)量、開產(chǎn)體質(zhì)量加性遺傳與母體遺傳存在負相關(guān)。開產(chǎn)體質(zhì)量、開產(chǎn)蛋質(zhì)量與開產(chǎn)日齡之間的遺傳相關(guān)系數(shù)分別為0.36±0.07、0.33±0.07,開產(chǎn)蛋質(zhì)量與開產(chǎn)體質(zhì)量之間的遺傳相關(guān)系數(shù)為0.65±0.05?!窘Y(jié)論】如皋黃雞開產(chǎn)日齡、開產(chǎn)體質(zhì)量、開產(chǎn)蛋質(zhì)量等開產(chǎn)性狀均受母體效應(yīng)影響,且遺傳力較低。

    關(guān)鍵詞:如皋黃雞; 遺傳力; 繁殖性狀; 母體效應(yīng); 動物模型

    遺傳力是預測群體選擇反應(yīng)、估計個體育種值的關(guān)鍵參數(shù),當前遺傳力研究熱點之一是分析母體效應(yīng)對加性遺傳的影響。開產(chǎn)性狀是蛋雞育種研究熱點之一,不僅反映性成熟狀況,還與早期產(chǎn)蛋性能密切相關(guān)。如皋黃雞主要分布于長江下游流域,屬于肉蛋兼用型品種,具有開產(chǎn)早、蛋黃比例高、肉質(zhì)優(yōu)良的特性。因此,系統(tǒng)研究如皋黃雞繁殖性能遺傳力,對于如皋黃雞開發(fā)利用具有重要意義。近年來,生態(tài)學、進化生物學以及基因組學研究表明鳥類生長、生殖性狀受母體效應(yīng)影響。據(jù)Mousseau等[1]研究,雌鳥可依據(jù)環(huán)境狀況改變后代數(shù)量和初生體質(zhì)量。雌鳥可以改變禽蛋中非DNA物質(zhì)含量,從而影響后代的生長發(fā)育。例如,蛋黃中生殖激素、類胡蘿卜素以及維生素影響后代生長、免疫、行為、形態(tài)以及生活力,母體依據(jù)環(huán)境條件調(diào)節(jié)這些物質(zhì)在禽蛋中的含量[2-3]。據(jù)報道雞基因組存在母體遺傳現(xiàn)象[4-5]。雞1號染色體蛋質(zhì)量、蛋品質(zhì)QTL只是母系來源的等位基因表達。雞3號染色體也存在基因組印記,屠宰性狀以及蛋質(zhì)量QTL是單個等位基因表達。研究表明,蛋雞[6]、肉雞[7]、火雞[8]、 鴨子[9]、鵝[10]、 鵪鶉[11]體質(zhì)量性狀遺傳力受母體效應(yīng)影響。綜上所述,多個研究團隊應(yīng)用不同的研究手段證明禽蛋在介導鳥類母體效應(yīng)中發(fā)揮重要作用。然而,目前的家禽開產(chǎn)性狀遺傳力計算中很少考慮母體效應(yīng),因而可能導致遺傳力估值過高。本研究以如皋黃雞選育群體為材料,收集2009—2012年生產(chǎn)數(shù)據(jù),檢驗加性遺傳、母體效應(yīng)對如皋黃雞繁殖性狀的影響并估計開產(chǎn)日齡、開產(chǎn)蛋質(zhì)量及開產(chǎn)體質(zhì)量遺傳力。研究結(jié)果將為如皋黃雞選育提供支持,為其他地方雞種遺傳參數(shù)研究提供參考。

    1材料與方法

    1.1試驗動物

    試驗數(shù)據(jù)集來自揚州翔龍禽業(yè)有限公司。如皋黃雞產(chǎn)蛋期飼養(yǎng)于指定雞舍,育成期光照條件為8 h光照∶16 h黑暗,產(chǎn)蛋期光照條件為16 h光照∶8 h黑暗。全程籠養(yǎng),常規(guī)蛋雞飼料。試驗群體測定開產(chǎn)日齡、開產(chǎn)蛋質(zhì)量和開產(chǎn)體質(zhì)量數(shù)據(jù),其中開產(chǎn)蛋質(zhì)量是最初3枚蛋平均值。選取2009—2012年連續(xù)4代有系譜記錄和生產(chǎn)數(shù)據(jù)的資料,原始數(shù)據(jù)集由來自10 454只雞的31 362條生產(chǎn)記錄組成,經(jīng)分析整理,去除離群值,其中,開產(chǎn)日齡、開產(chǎn)蛋質(zhì)量、開產(chǎn)體質(zhì)量的記錄數(shù)分別為10 352、10 335、10 378條,父系數(shù)為242個,開產(chǎn)日齡、開產(chǎn)蛋質(zhì)量數(shù)據(jù)來自1 494個母系,開產(chǎn)體質(zhì)量數(shù)據(jù)來自1 493個母系。數(shù)據(jù)統(tǒng)計信息詳見表1。

    表1如皋黃雞開產(chǎn)性狀數(shù)據(jù)描述

    Tab.1Date description statistics of reproduction traits in Rugao yellow chickens

    項目開產(chǎn)日齡/d開產(chǎn)蛋質(zhì)量/g開產(chǎn)體質(zhì)量/g最大值18144.001948.00最小值12322.00857.00平均值±標準差150.69±9.3432.99±3.19 1350.93±161.38CV/%6.209.6711.95

    1.2統(tǒng)計模型

    首先,應(yīng)用SAS GLM分析表型方差,以生產(chǎn)記錄包含的年代、批次、籠位信息對表型值進行方差分析,確定進入固定效應(yīng)的因素。其次,應(yīng)用WOMBAT軟件,采用平均信息最大約束似然法,分析各性狀(協(xié))方差組分以及遺傳力[12]。以是否包含母體遺傳、母體環(huán)境效應(yīng)為區(qū)分標準,建立4個動物模型:

    模型I(model I):Y=Xb+Zaa+e,

    模型II(modelII):Y=Xb+Zaa+Zcc+e,

    模型III(modelIII):Y=Xb+Zaa+Zmm+e,

    模型IV(modelIV):Y=Xb+Zaa+Zcc+Zmm+e,

    式中:Y是所有個體性狀表型值向量;a、b、c、m、e分別是加性遺傳效應(yīng)、固定效應(yīng)、母體環(huán)境效應(yīng)、母體遺傳效應(yīng)以及殘差效應(yīng)的向量。X、Za、Zm、Zc分別是固定效應(yīng)、加性遺傳效應(yīng)、母體遺傳效應(yīng)以及母體環(huán)境效應(yīng)的指定矩陣。最后,以似然比比較各性狀適用模型[13-14]。似然比檢驗統(tǒng)計量為:

    -2/log(e)×[logLmax|modelⅠ-logLmax|modelⅡ]=

    -4.605 2×[logLmax|modelⅠ-logLmax|modelⅡ],

    2結(jié)果與分析

    2.1影響開產(chǎn)性狀表型方差的非遺傳因素

    表型值具有一定的彈性空間,其中環(huán)境是影響性狀表現(xiàn)的重要因素。如表2所示,如皋黃雞開產(chǎn)性狀按環(huán)境因素分組后平均值和標準差存在顯著差異。批次指定依照時間順序,距離當下最近的為第1批,最遠的是第10批。批次之間,飼養(yǎng)管理、氣溫、濕度條件存在差異,因此各批次開產(chǎn)性狀的平均值表現(xiàn)顯著差異。開產(chǎn)早的批次,平均開產(chǎn)蛋質(zhì)量較小。1~3批是2012年7月13日—2012年8月13日出雛。4~6批是2011年8月6日—2011年9月8日出雛,批次之間時間間隔較短,育雛、育成管理較好,因而開產(chǎn)體質(zhì)量差異不顯著,且平均值較大。7~9批是2010年6月30日—2010年8月16日出雛,批次之間時間間隔長,所以批次間開產(chǎn)體質(zhì)量差異顯著。 綜合考慮,(協(xié))方差組分分析時將批次因素列入固定效應(yīng)。

    表2不同批次如皋黃雞開產(chǎn)性狀平均值和標準差1)

    Tab.2Means and standard deviations of reproduction traits of Rugao yellow chickens from different batches

    批次開產(chǎn)日齡/d開產(chǎn)蛋質(zhì)量/g開產(chǎn)體質(zhì)量/g1146.79±9.29g31.97±2.91g1288.76±145.18ef2149.00±9.22f32.69±3.03e1301.19±143.31e3153.88±9.30c32.96±3.00ed1283.50±145.38f4151.39±6.25d33.04±3.07d1418.14±173.35a5150.47±8.29e33.71±2.97b1427.96±169.31a6155.12±8.80b34.27±3.18a1421.30±165.71a7156.14±7.28a33.71±3.26b1412.78±166.21a8144.65±8.89h32.38±3.45f1325.16±137.80d9154.73±7.10b33.34±2.91c1376.48±146.43b10144.21±7.35h32.34±3.49f1343.40±156.06c

    1)同列數(shù)據(jù)后凡是有一個相同小寫字母者,表示在0.05水平上差異不顯著(單因素方差分析SNK法)。

    2.2開產(chǎn)性狀遺傳參數(shù)

    表3如皋黃雞開產(chǎn)性狀不同模型的方差組分和遺傳力1)

    Tab.3Variance components and heritabilities of reproduction traits of Rugao yellow chickens estimated with different models

    性狀模型σ2aσ2mσ2cσ2eσ2ph2a±SEh2m±SEc2a±SEram開產(chǎn)日齡Ⅰ21.7349.3171.050.31±0.02Ⅱ19.741.2149.7270.670.28±0.020.02±0.01Ⅲ20.600.6649.7070.960.29±0.020.01±0.00-0.04Ⅳ19.670.191.1849.6370.670.28±0.0300.02±0.01-0.12開產(chǎn)蛋質(zhì)量Ⅰ2.397.6110.000.24±0.02Ⅱ2.270.097.629.980.23±0.020.01±0.01Ⅲ2.310.267.4410.010.23±0.020.03±0.01-0.43Ⅳ2.250.210.107.429.990.23±0.020.02±0.010.01±0.01-0.57開產(chǎn)體質(zhì)量Ⅰ9735.2014970.0024705.000.39±0.02Ⅱ8256.60875.3315289.0024421.000.34±0.030.04±0.01Ⅲ8960.051186.9514533.0024680.000.36±0.030.05±0.01-0.35Ⅳ8298.22861.181092.0014225.0024476.000.34±0.030.04±0.010.04±0.01-0.92

    如皋黃雞開產(chǎn)日齡加性遺傳效應(yīng)與母體遺傳效應(yīng)之間存在較弱負相關(guān)(ram為-0.04~-0.12,見表3),表明選育提高開產(chǎn)日齡加性遺傳效應(yīng)對母體遺傳效應(yīng)幾乎沒有影響。開產(chǎn)蛋質(zhì)量的ram為-0.43~-0.57,表明選育提高開產(chǎn)蛋質(zhì)量加性遺傳將降低母體遺傳效應(yīng)。對于開產(chǎn)體質(zhì)量而言,加性遺傳效應(yīng)與母體遺傳效應(yīng)相關(guān)系數(shù)為-0.35~-0.92。

    以多性狀模型計算的遺傳力與表3列出的相差不大。由表4可知,表型相關(guān)系數(shù)小于遺傳相關(guān)系數(shù),表明受環(huán)境相關(guān)影響。

    表4開產(chǎn)性狀之間遺傳相關(guān)(左下)和表型相關(guān)(右上)1)

    Tab.4Estimated values of genetic (lower left) and phenotypic (upper right) correlations between traits at the first egg

    性狀開產(chǎn)日齡開產(chǎn)體質(zhì)量開產(chǎn)蛋質(zhì)量開產(chǎn)日齡0.28±0.031)0.19±0.010.23±0.01開產(chǎn)體質(zhì)量0.36±0.070.35±0.030.31±0.01開產(chǎn)蛋質(zhì)量0.33±0.070.65±0.050.23±0.02

    1) 對角線上標示下劃線數(shù)值為遺傳力。

    2.3動物模型比較

    表5遺傳評估模型最大似然值對數(shù)

    Tab.5The maximum logarithmic likelihood value for each model

    性狀模型Ⅰ模型Ⅱ模型Ⅲ模型Ⅳ開產(chǎn)日齡-26751.30-26748.28-26749.80-26748.27開產(chǎn)蛋質(zhì)量-16731.16-16730.40-16729.22-16728.67開產(chǎn)體質(zhì)量-56936.22-56924.53-56929.16-56920.24多性狀分析-97958.88-97938.14-97943.16-97934.66

    表6模型之間似然比值及2檢驗結(jié)果1)

    Tab.6Likelihood ratios and2 test results among different models

    模型比較開產(chǎn)日齡開產(chǎn)蛋質(zhì)量開產(chǎn)體質(zhì)量多性狀分析Ⅰ∶Ⅱ13.91**3.4853.85**95.51**Ⅰ∶Ⅲ6.91*8.94*32.52**72.41**Ⅰ∶Ⅳ13.95**11.46**73.60**111.55**Ⅱ∶Ⅲ|-7.00|**5.46*|-21.33|**Ⅱ∶Ⅳ0.057.98*19.75**16.04*Ⅲ∶Ⅳ7.05**2.2541.07**39.14**

    1)*表示0.05水平差異顯著;**表示0.01水平差異顯著。

    3討論與結(jié)論

    3.1開產(chǎn)性狀遺傳力

    對于開產(chǎn)日齡,2000年之前的研究認為其遺傳力較高,范圍從0.31~0.48[15-17]。近年,以亞洲地方雞種為試材遺傳分析表明開產(chǎn)日齡屬于低遺傳力性狀,其值從0.12~0.32[18-20]。本研究中,如皋黃雞開產(chǎn)日齡遺傳力估值除與東亞地方雞種的近似外,低于其他研究估值。其原因除了本研究考慮了母體效應(yīng)對遺傳力影響外,遺傳背景、統(tǒng)計方法及模型等因素也可能會造成遺傳力估值上的不同。最近,全基因組關(guān)聯(lián)研究解析出的原因性變異較少,預示以往開產(chǎn)日齡遺傳力可能被過高估算[21- 22]。

    對于如皋黃雞,開產(chǎn)蛋質(zhì)量遺傳力為0.23左右。本研究結(jié)果略高于Wolc等[23]、Sang等[19]、Hosseini 等[24]以及Niknafs等[25]的研究結(jié)果,略低于吳桂琴等[26]的研究結(jié)果。與產(chǎn)蛋高峰期相比,開產(chǎn)階段蛋質(zhì)量遺傳力較低,隨著周齡增加,蛋質(zhì)量遺傳力逐漸升高,并趨于穩(wěn)定。產(chǎn)蛋初期,卵巢、下丘腦、垂體以及能量供給器官仍以性成熟為主題發(fā)揮自身功能,在30周齡左右,產(chǎn)蛋雞完成性成熟過程,一些與產(chǎn)蛋相關(guān)基因開始表達,加性遺傳力逐步增加[27]。

    生長性狀遺傳力一般高于生殖性狀遺傳力。如皋黃雞開產(chǎn)體質(zhì)量遺傳力為0.34±0.03,高于開產(chǎn)日齡和開產(chǎn)蛋質(zhì)量遺傳力。如皋黃雞開產(chǎn)體質(zhì)量遺傳力低于伊朗地方品種Far雞(h2=0.57)、Mazandaran雞(h2=0.48),但3個亞洲地方品種均檢測到母體遺傳的影響[20, 28]。如皋黃雞開產(chǎn)體質(zhì)量遺傳力低于峪口褐殼蛋雞父系(h2=0.60)[26],高于新楊白殼蛋雞(h2=0.20)[29],與韓國地方雞種的相近(h2=0.39)[19]。

    3.2動物模型比較

    以模型I作為子模時,似然比值與其他3個模型存在極顯著差異,表明模型I不適于如皋黃雞開產(chǎn)性狀遺傳力估計。對于開產(chǎn)日齡,模型III和IV之間存在顯著差異,表明母體環(huán)境效應(yīng)有著顯著影響;模型II與III之間為負值,如果調(diào)換子模位置,則模型之間存在顯著差異,支持母體環(huán)境效應(yīng)影響下顯著。兩者綜合,可選擇模型II為開產(chǎn)日齡適用動物模型。

    對于開產(chǎn)蛋質(zhì)量,除去模型I之外,模型II與模型III、IV之間存在顯著差異,表明模型II并非適用模型;模型III、IV之間似然比值差異不顯著,殘差比值以及母體遺傳效應(yīng)相差不大,參數(shù)少的模型III是開產(chǎn)蛋質(zhì)量適用模型。

    對于開產(chǎn)體質(zhì)量,模型II與III互換子模位置,則2個模型之間存在顯著差異,模型II好于模型III;除模型I之外,3種模型兩兩間存在顯著差異,所以選擇模型IV為最適模型。

    對肉蛋兼用型品種如皋黃雞來說,開產(chǎn)性狀選育具有重要經(jīng)濟意義。研究結(jié)果表明如皋黃雞開產(chǎn)性狀屬于低遺傳力性狀,受環(huán)境因素影響較大,因而應(yīng)用傳統(tǒng)表型選擇很難獲得良好遺傳進展。母體效應(yīng)影響如皋黃雞開產(chǎn)性狀,計算育種值、選擇反應(yīng)時應(yīng)考慮母體效應(yīng)對遺傳力的影響。如皋黃雞開產(chǎn)日齡、開產(chǎn)蛋質(zhì)量、開產(chǎn)體質(zhì)量遺傳力分別為0.28、0.23、0.34。

    參考文獻:

    [1]MOUSSEAU T A, FOX C W. The adaptive significance of maternal effects[J]. Trends Ecol Evol, 1998, 13(10): 403-407.

    [2]BIARD C, GIL D, KARADAF, et al. Maternal effects mediated by antioxidants and the evolution of carotenoid-based signals in birds[J]. Am Nat, 2009, 174(5): 696-708.

    [3]NEWBREY J L, REED W L. Growth of yellow-headed blackbirdXanthocephalusxanthocephalusnestlings in relation to maternal body condition, egg mass, and yolk carotenoids concentrations[J]. J Avian Biol, 2009, 40(4): 419-429.

    [4]TUISKULA-HAAVISTO M, VILKKI J. Parent-of-origin specific QTL--a possibility towards understanding reciprocal effects in chicken and the origin of imprinting[J]. Cytogenet Genome Res, 2007, 117(1/2/3/4): 305-312.

    [5]MCELROY J P, KIM J J, HARRY D E, et al. Identification of trait loci affecting white meat percentage and other growth and carcass traits in commercial broiler chickens[J]. Poult Sci, 2006, 85(4): 593-605.

    [6]HARTMANN C, JOHANSSON K, STRANDBERG E, et al. Genetic correlations between the maternal genetic effect on chick weight and the direct genetic effects on egg composition traits in a White Leghorn line[J]. Poult Sci, 2003, 82(1): 1-8.

    [7]MANIATIS G, DEMIRIS N, KRANIS A, et al. Model comparison and estimation of genetic parameters for body weight in commercial broilers[J]. Can J Anim Sci, 2013, 93(1): 67-77.

    [8]ASLAM M L, BASTIAANSEN J W M, CROOIJMANS R P M A, et al. Genetic variances, heritabilities and maternal effects on body weight, breast meat yield, meat quality traits and the shape of the growth curve in turkey birds[J]. BMC Genet, 2011,[2015-06-06].http://www.biomedcentral.com/1471-2156/12/14. doi:10.1186/1471-2156-12-14.

    [9]SZWACZKOWSKI T, GRZECH M, BOROWSKA A, et al. Maternal genetic effects on body weight and breast morphological traits in duck population under selection[J]. Arch Tierz-Arch Anim Breed, 2010, 53(5): 600-608.

    [10]SZWACZKOWSKI T, WEZYK S, STANISLAWSKA-BARCZAK E, et al. Genetic variability of body weight in two goose strains under long-term selection[J]. J Appl Genet, 2007, 48(3): 253-260.

    [11]LOTFI E, ZEREHDARAN S, AZARI M A. Direct and maternal genetic effects of body weight traits in Japanese quail (Coturnixcoturnixjaponica)[J]. Arch Geflugelkd, 2012, 76(3): 150-154.

    [12]MEYER K. WOMBAT:A tool for mixed model analyses in quantitative genetics by restricted maximum likelihood (REML)[J]. J Zhejiang Univ Sci B, 2007, 8(11): 815-821.

    [13]白俊艷,李金泉,賈小平,等. 母體遺傳效應(yīng)對絨山羊生產(chǎn)性狀遺傳參數(shù)估計的影響[J]. 遺傳, 2006, 28(9): 1083-1086.

    [14]MEYER K. Likelihood calculations to evaluate experimental designs to estimate genetic variances[J]. Heredity, 2008, 101(3): 212-221.

    [15]POGGENPOEL D G, DUCKITT J S. Genetic basis of the increase in egg weight with pullet age in a White Leghorn flock[J]. Brit Poultry Sci, 1988, 29(4): 863-867.

    [16]KOERHUIS A N M, MCKAY J C. Restricted maximum likelihood estimation of genetic parameters for egg production traits in relation to juvenile body weight in broiler chickens[J]. Livest Prod Sci, 1996, 46(2): 117-127.

    [17]JEYARUBAN M, GIBSON J. Estimation of additive genetic variance in commercial layer poultry and simulated populations under selection[J]. Theor Appl Genet, 1996, 92(3/4): 483-491.

    [18]張龍,劉益平,趙小玲,等. 二郎山山地雞兩品系主要繁殖性狀遺傳參數(shù)的估計及分析[J]. 中國家禽, 2012, 34(14): 69-71.

    [19]SANG B D, KONG H S, KIM H K, et al. Estimation of genetic parameters for economic traits in Korean native chickens[J]. Asian Australas J Anim Sci, 2006, 19(3): 319-323.

    [20]GHORBANI S, KAMALI M A, ABBASI M A, et al. Estimation of maternal effects on some economic traits of north iranian native fowls using different models[J]. J Agric Sci Technol, 2012, 14(1): 95-103.

    [21]SHEN X, ZENG H, XIE L, et al. The GTPase activating Rap/RanGAP domain-like 1 gene is associated with chicken reproductive traits[J]. PLoS One, 2012, 7(4): e33851.

    [22]LIU W, LI D, LIU J, et al. A genome-wide SNP scan reveals novel loci for egg production and quality traits in white leghorn and brown-egg dwarf layers[J]. PLoS One, 2011, 6(12): e28600.

    [23]WOLC A, ARANGO J, SETTAR P, et al. Genetic parameters of egg defects and egg quality in layer chickens[J]. Poult Sci, 2012, 91(6): 1292-1298.

    [24]HOSSEINI S, TAHMOORESPUR M. Estimation of genetic parameters of various economic traits in west Iranian native chickens[J]. Indian J Anim Res, 2013, 47(3): 205-211.

    [25]NIKNAFS S, NEJATI-JAVAREMI A, MEHRABANI-YEGANEH H, et al. Estimation of genetic parameters for body weight and egg production traits in Mazandaran native chicken[J]. Trop Anim Health Prod, 2012, 44(7): 1437-1443.

    [26]吳桂琴, 李花妮, 劉會英, 等. 峪口褐殼純系蛋雞父系主要性狀的遺傳參數(shù)分析[J]. 中國家禽, 2015, 37(9): 49-51.

    [27]WOLC A, SZWACZKOWSKI T. Estimation of genetic parameters for monthly egg production in laying hens based on random regression models[J]. J Appl Genet, 2009, 50(1): 41-46.

    [28]GHORBANI S, TAHMOORESPUR M, MAGHSOUDI A, et al. Estimates of (co)variance components for production and reproduction traits with different models in Fars native fowls[J]. Livest Sci, 2013, 151(2/3): 115-123.

    [29]田春雷, 嚴華祥, 楊長鎖, 等. 新楊白殼純系蛋雞主要性狀的遺傳參數(shù)估計[J]. 上海畜牧獸醫(yī)通訊, 2009(5): 9-11.

    【責任編輯柴焰】

    Direct and maternal genetic effects for traits at the first egg in Rugao yellow chickens

    GUO Jun, QU Liang, SHEN Manman, DOU Taocun, HU Yuping, WANG Kehua

    (Institute of Poultry Science, Chinese Academy of Agricultural Science, Yangzhou 225125, China)

    Abstract:【Objective】 The aims of this study were to evaluate direct and maternal genetic effects for traits at the first egg of Rugao yellow chickens with animal models, and to estimate the genetic parameters for these traits. 【Method】Reproduction data, including data for 10 454 individual animals and 31 362 reproduction records, of Rugao yellow chickens from 2009 to 2012 were collected. Four animal models with and without maternal effects were used to estimate the genetic parameters. Phenotypic variance of reproduction traits was partitioned using average information restricted maximum likelihood (AIREML) method with WOMBAT software. The models were evaluated with the likelihood ratio test. Univariate models were used to estimate heritabilities, and multivariate models were used to analyze genetic and phenotypic correlations.【Result】Maternal environment had significant influence on age at the first egg (AFE), for which the optimum model was model II and the heritability was 0.28±0.02. Maternal genetic effects had significant influence on egg mass at the first egg (EM1), for which the optimum model was model III, the direct heritability was 0.23±0.02 and the maternal genetic heritability was 0.03±0.01. Both maternal environment and maternal genetic effects had significant influence on body mass at the first egg (BM1), for which the optimum model was model IV, the heritability was 0.34±0.03 and the maternal genetic heritability was 0.04±0.01. There were negative correlations between the direct and maternal genetic effects for EM1 and BM1. The genetic correlation coefficients between AFE and BM1, AFE and EM1, and EM1 and BM1 were 0.36±0.07, 0.33±0.07, and 0.65±0.05, respectively.【Conclusion】The traits at the first egg of Rugao yellow chickens are influenced by maternal effects, and are lowly heritable.

    Key words:Rugao yellow chicken; heritability; reproduction trait; maternal effect; animal model

    中圖分類號:S831

    文獻標志碼:A

    文章編號:1001- 411X(2016)02- 0001- 06

    基金項目:現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項資金(CARS-41-K02);江蘇省自然科學基金(BK20131237);江蘇省科技成果轉(zhuǎn)化專項資金(BA2014132)

    作者簡介:郭軍(1971—),男,副研究員,博士,E-mail:guojun.yz@gmail.com;通信作者:王克華(1962—),男,研究員,博士,E-mail:sqbreeding@126.com

    收稿日期:2015- 06- 27優(yōu)先出版時間:2016- 01- 18

    優(yōu)先出版網(wǎng)址:http://www.cnki.net/kcms/detail/44.1110.s.20160118.1657.034.html

    郭軍, 曲亮, 沈曼曼,等.如皋黃雞開產(chǎn)性狀加性遺傳及母體遺傳效應(yīng)分析[J].華南農(nóng)業(yè)大學學報,2016,37(2):1- 6.

    猜你喜歡
    遺傳力動物模型
    某國家生豬核心育種場加系大白母豬繁殖性狀遺傳參數(shù)估計
    肥胖中醫(yī)證候動物模型研究進展
    美系大白種豬生長性狀和繁殖性狀遺傳力估計
    養(yǎng)豬(2022年4期)2022-08-17 07:07:02
    胃癌前病變動物模型復制實驗進展
    潰瘍性結(jié)腸炎動物模型研究進展
    不同類型煙草煙葉氯含量性狀的遺傳分析
    通過胎盤效率改良母豬繁殖性能的研究進展
    糖尿病性視網(wǎng)膜病變動物模型研究進展
    抗蟲棉品種系主要農(nóng)藝性狀配合力與遺傳力分析
    血管性癡呆動物模型的建立與中醫(yī)藥干預進展
    超碰av人人做人人爽久久| 少妇的逼水好多| 国产精品熟女久久久久浪| 男插女下体视频免费在线播放| 成人午夜高清在线视频| 97热精品久久久久久| 免费观看性生交大片5| 亚洲精品乱久久久久久| 亚洲精品成人久久久久久| 中文欧美无线码| 人人妻人人澡人人爽人人夜夜 | 欧美高清成人免费视频www| 久久久久久久亚洲中文字幕| av黄色大香蕉| av在线老鸭窝| 久久久久久久亚洲中文字幕| 国产淫语在线视频| 欧美97在线视频| 真实男女啪啪啪动态图| 亚洲av在线观看美女高潮| 一本久久精品| 国产一级毛片七仙女欲春2| 欧美激情久久久久久爽电影| 搞女人的毛片| 精品酒店卫生间| 国内精品美女久久久久久| 成年女人看的毛片在线观看| 精品少妇黑人巨大在线播放| 好男人在线观看高清免费视频| 精品国产露脸久久av麻豆 | 老师上课跳d突然被开到最大视频| 2021天堂中文幕一二区在线观| 18+在线观看网站| 人妻少妇偷人精品九色| 精品人妻视频免费看| 男女边摸边吃奶| 小蜜桃在线观看免费完整版高清| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美日韩在线观看h| 能在线免费看毛片的网站| 综合色丁香网| 国产亚洲最大av| 97人妻精品一区二区三区麻豆| 亚洲内射少妇av| 26uuu在线亚洲综合色| 大香蕉久久网| 欧美日本视频| 久久久久久久大尺度免费视频| 精品人妻熟女av久视频| 亚洲四区av| 男人爽女人下面视频在线观看| 日韩 亚洲 欧美在线| 亚洲精品影视一区二区三区av| 日韩欧美精品免费久久| 简卡轻食公司| 极品教师在线视频| av黄色大香蕉| 看黄色毛片网站| 看免费成人av毛片| 国产成人精品久久久久久| 在线天堂最新版资源| 高清午夜精品一区二区三区| 网址你懂的国产日韩在线| 国产精品一区二区三区四区免费观看| 看免费成人av毛片| 身体一侧抽搐| 色网站视频免费| 特级一级黄色大片| 国产高清不卡午夜福利| 肉色欧美久久久久久久蜜桃 | 国产毛片a区久久久久| 99视频精品全部免费 在线| 男的添女的下面高潮视频| 自拍偷自拍亚洲精品老妇| 精品国产三级普通话版| 2018国产大陆天天弄谢| 少妇裸体淫交视频免费看高清| 大片免费播放器 马上看| 好男人在线观看高清免费视频| 亚洲精品中文字幕在线视频 | 又爽又黄a免费视频| 亚洲人成网站在线观看播放| 国产成人精品福利久久| 日韩欧美 国产精品| 亚洲最大成人中文| av在线蜜桃| 国产精品福利在线免费观看| 搡女人真爽免费视频火全软件| 亚洲精品色激情综合| 床上黄色一级片| 校园人妻丝袜中文字幕| 青春草亚洲视频在线观看| 国产av不卡久久| 搡老妇女老女人老熟妇| 婷婷色综合www| 韩国高清视频一区二区三区| 日日啪夜夜爽| 国产免费又黄又爽又色| 久久精品国产亚洲av涩爱| 99re6热这里在线精品视频| 亚洲电影在线观看av| 最后的刺客免费高清国语| 久久精品夜夜夜夜夜久久蜜豆| xxx大片免费视频| 一区二区三区免费毛片| 久久久久久久大尺度免费视频| 欧美丝袜亚洲另类| 一级毛片电影观看| 菩萨蛮人人尽说江南好唐韦庄| 国产69精品久久久久777片| 观看美女的网站| 一区二区三区四区激情视频| 成人漫画全彩无遮挡| videossex国产| 人妻少妇偷人精品九色| 免费黄网站久久成人精品| 99久国产av精品国产电影| 国产精品日韩av在线免费观看| 久久6这里有精品| 日本欧美国产在线视频| 国产黄色视频一区二区在线观看| 在线 av 中文字幕| 午夜精品在线福利| 成人漫画全彩无遮挡| 亚洲内射少妇av| 亚洲国产最新在线播放| 大片免费播放器 马上看| 亚洲av免费高清在线观看| 精品欧美国产一区二区三| 午夜精品一区二区三区免费看| 亚洲精品国产av成人精品| 简卡轻食公司| 男女下面进入的视频免费午夜| 成人漫画全彩无遮挡| 男女视频在线观看网站免费| 97精品久久久久久久久久精品| 偷拍熟女少妇极品色| 一区二区三区免费毛片| 欧美zozozo另类| 91久久精品国产一区二区三区| 欧美成人午夜免费资源| 22中文网久久字幕| 乱人视频在线观看| 国产亚洲av片在线观看秒播厂 | .国产精品久久| 久久久亚洲精品成人影院| 国产成人91sexporn| 久久久久久国产a免费观看| 国产精品久久久久久av不卡| 亚洲一区高清亚洲精品| 久久这里有精品视频免费| 中文在线观看免费www的网站| 不卡视频在线观看欧美| 亚洲国产日韩欧美精品在线观看| 国产色爽女视频免费观看| 成年女人在线观看亚洲视频 | 直男gayav资源| 天天一区二区日本电影三级| 国内少妇人妻偷人精品xxx网站| 免费播放大片免费观看视频在线观看| 国产精品爽爽va在线观看网站| 日本一二三区视频观看| 久久精品夜夜夜夜夜久久蜜豆| 国产av不卡久久| 一本一本综合久久| 观看美女的网站| 国产69精品久久久久777片| 久久久久网色| 成人二区视频| 亚洲图色成人| 国产黄片美女视频| 日韩中字成人| 国产精品一及| 亚洲精品色激情综合| 最新中文字幕久久久久| 九九在线视频观看精品| 永久网站在线| 99re6热这里在线精品视频| 日韩av在线免费看完整版不卡| 国内精品美女久久久久久| 黄片无遮挡物在线观看| 久久久精品94久久精品| 国产黄频视频在线观看| 综合色av麻豆| 男女视频在线观看网站免费| 亚洲国产成人一精品久久久| 久久国内精品自在自线图片| 伊人久久国产一区二区| 五月天丁香电影| 中文字幕久久专区| 国产精品不卡视频一区二区| 国产黄频视频在线观看| 十八禁网站网址无遮挡 | 夫妻午夜视频| 国产一区有黄有色的免费视频 | 免费观看在线日韩| 亚洲欧美成人综合另类久久久| 国产精品av视频在线免费观看| 欧美zozozo另类| 亚洲熟女精品中文字幕| 亚洲熟妇中文字幕五十中出| 欧美成人一区二区免费高清观看| 欧美成人精品欧美一级黄| 赤兔流量卡办理| 一级av片app| 黄色欧美视频在线观看| 成人亚洲精品av一区二区| 搞女人的毛片| 亚洲国产高清在线一区二区三| 狂野欧美白嫩少妇大欣赏| 久久国内精品自在自线图片| 搡老妇女老女人老熟妇| 韩国av在线不卡| 日日啪夜夜爽| 亚洲av不卡在线观看| 欧美变态另类bdsm刘玥| 夜夜看夜夜爽夜夜摸| 国产熟女欧美一区二区| 黄片无遮挡物在线观看| 国产精品女同一区二区软件| 亚洲精品影视一区二区三区av| 性插视频无遮挡在线免费观看| 特大巨黑吊av在线直播| 午夜激情欧美在线| 久久久a久久爽久久v久久| 精品国产一区二区三区久久久樱花 | 免费少妇av软件| 黄色配什么色好看| av福利片在线观看| 亚洲怡红院男人天堂| 蜜臀久久99精品久久宅男| 国产伦精品一区二区三区视频9| 搡老乐熟女国产| 69av精品久久久久久| 两个人视频免费观看高清| 亚洲国产成人一精品久久久| 国产精品不卡视频一区二区| 久久久亚洲精品成人影院| 欧美人与善性xxx| 久久99精品国语久久久| 国产精品三级大全| 久久久久精品性色| 伊人久久国产一区二区| 美女主播在线视频| 亚洲精品国产av成人精品| 99久国产av精品| 国精品久久久久久国模美| 国产单亲对白刺激| 男女视频在线观看网站免费| 国产亚洲91精品色在线| 高清日韩中文字幕在线| 我要看日韩黄色一级片| 国产在视频线精品| 麻豆成人午夜福利视频| 久久草成人影院| 99re6热这里在线精品视频| 国产亚洲5aaaaa淫片| 成人毛片60女人毛片免费| 高清av免费在线| eeuss影院久久| av在线播放精品| 人妻系列 视频| 亚洲欧美一区二区三区国产| 亚洲四区av| 建设人人有责人人尽责人人享有的 | 国产白丝娇喘喷水9色精品| 边亲边吃奶的免费视频| 久久精品久久久久久噜噜老黄| 九九在线视频观看精品| 成人美女网站在线观看视频| 天堂网av新在线| 黄色日韩在线| 人人妻人人澡欧美一区二区| av福利片在线观看| 久久99热这里只有精品18| 噜噜噜噜噜久久久久久91| 51国产日韩欧美| 久久亚洲国产成人精品v| 亚洲国产高清在线一区二区三| 亚洲国产av新网站| 色哟哟·www| 一本一本综合久久| 亚洲在线观看片| 欧美成人精品欧美一级黄| 伊人久久精品亚洲午夜| 国产高潮美女av| 色综合站精品国产| 国产美女午夜福利| 国产在线男女| 国产淫语在线视频| 在线天堂最新版资源| 久久这里只有精品中国| 国产精品久久久久久精品电影| 婷婷色综合www| 国产一区有黄有色的免费视频 | 一个人观看的视频www高清免费观看| 大陆偷拍与自拍| 国产免费福利视频在线观看| 美女高潮的动态| 久久久午夜欧美精品| .国产精品久久| 99热网站在线观看| videossex国产| 大片免费播放器 马上看| 别揉我奶头 嗯啊视频| 久久久久久久亚洲中文字幕| 丝袜喷水一区| 国内少妇人妻偷人精品xxx网站| 日本午夜av视频| 一区二区三区免费毛片| 在线免费观看不下载黄p国产| 老司机影院成人| 国产伦理片在线播放av一区| av网站免费在线观看视频 | 极品少妇高潮喷水抽搐| 肉色欧美久久久久久久蜜桃 | 天天躁夜夜躁狠狠久久av| 亚洲av男天堂| 男人舔女人下体高潮全视频| 美女内射精品一级片tv| 欧美一级a爱片免费观看看| 日韩av不卡免费在线播放| 99九九线精品视频在线观看视频| 成人亚洲精品av一区二区| .国产精品久久| 日韩不卡一区二区三区视频在线| 嫩草影院新地址| 精品久久国产蜜桃| 成年人午夜在线观看视频 | 亚洲av在线观看美女高潮| 亚洲精品一二三| av在线观看视频网站免费| 亚洲精品自拍成人| 国产免费视频播放在线视频 | 久久久久久九九精品二区国产| 淫秽高清视频在线观看| 狠狠精品人妻久久久久久综合| 国产精品99久久久久久久久| 亚洲精品自拍成人| 久久久久久伊人网av| 亚洲欧美中文字幕日韩二区| 丝袜喷水一区| 久久久午夜欧美精品| 少妇丰满av| 99久国产av精品国产电影| 国产一区有黄有色的免费视频 | 久久久久精品久久久久真实原创| 禁无遮挡网站| 亚洲美女搞黄在线观看| 女人久久www免费人成看片| 国产男人的电影天堂91| 日本免费a在线| 久久97久久精品| 乱人视频在线观看| 白带黄色成豆腐渣| 午夜福利在线观看免费完整高清在| 全区人妻精品视频| 亚洲精品中文字幕在线视频 | 老师上课跳d突然被开到最大视频| 国产男女超爽视频在线观看| 国产av在哪里看| 亚洲欧洲日产国产| 国产av在哪里看| 97热精品久久久久久| 国产精品麻豆人妻色哟哟久久 | av在线天堂中文字幕| 最后的刺客免费高清国语| 最近的中文字幕免费完整| 亚洲在线自拍视频| 成人毛片a级毛片在线播放| 亚洲av.av天堂| 国产精品一区二区三区四区久久| 毛片一级片免费看久久久久| 欧美日韩亚洲高清精品| 精品亚洲乱码少妇综合久久| 在线a可以看的网站| 国产成人免费观看mmmm| 国产高清国产精品国产三级 | 欧美+日韩+精品| 国产高潮美女av| 亚洲18禁久久av| 在线免费十八禁| 国产成人一区二区在线| 午夜福利高清视频| 久久综合国产亚洲精品| 一级毛片我不卡| 成人欧美大片| 2021天堂中文幕一二区在线观| 肉色欧美久久久久久久蜜桃 | 赤兔流量卡办理| 日韩欧美 国产精品| 中文精品一卡2卡3卡4更新| 中文字幕av成人在线电影| 国产精品三级大全| 99热6这里只有精品| 99久久精品国产国产毛片| 丝瓜视频免费看黄片| 国内精品美女久久久久久| 永久网站在线| 国国产精品蜜臀av免费| 色综合亚洲欧美另类图片| 亚洲精品日本国产第一区| 久久精品国产亚洲av天美| 麻豆乱淫一区二区| 久久久久久久久久久免费av| 久久久久国产网址| 亚洲av成人av| 欧美3d第一页| 国产高清不卡午夜福利| 亚洲欧洲国产日韩| 午夜激情欧美在线| 少妇的逼水好多| 人人妻人人看人人澡| 国产精品av视频在线免费观看| 免费看美女性在线毛片视频| 久久精品国产自在天天线| 国产 亚洲一区二区三区 | 欧美激情在线99| 亚洲在线自拍视频| 午夜福利在线在线| 综合色丁香网| 日本三级黄在线观看| 大陆偷拍与自拍| 一级a做视频免费观看| 亚洲国产色片| 免费av观看视频| 中文天堂在线官网| 欧美成人a在线观看| 欧美变态另类bdsm刘玥| 联通29元200g的流量卡| 美女脱内裤让男人舔精品视频| 久久精品综合一区二区三区| 插阴视频在线观看视频| 精品亚洲乱码少妇综合久久| 三级国产精品欧美在线观看| 日韩三级伦理在线观看| 久久久久久久大尺度免费视频| 欧美高清性xxxxhd video| 午夜日本视频在线| 少妇猛男粗大的猛烈进出视频 | 久久国内精品自在自线图片| 亚洲成色77777| 舔av片在线| 国国产精品蜜臀av免费| 老师上课跳d突然被开到最大视频| 国内精品美女久久久久久| 亚洲怡红院男人天堂| 国产单亲对白刺激| 中国美白少妇内射xxxbb| 成人特级av手机在线观看| 亚洲美女搞黄在线观看| 少妇人妻一区二区三区视频| 亚洲欧美清纯卡通| 国产伦在线观看视频一区| 亚洲自偷自拍三级| 天堂俺去俺来也www色官网 | 国产国拍精品亚洲av在线观看| 国产av国产精品国产| 看黄色毛片网站| www.色视频.com| 美女内射精品一级片tv| 久久久久久久久中文| 一级片'在线观看视频| 久久久久久久久久人人人人人人| 2021天堂中文幕一二区在线观| 亚洲国产色片| 夜夜看夜夜爽夜夜摸| 一夜夜www| 国产成人午夜福利电影在线观看| av免费在线看不卡| 亚洲av免费高清在线观看| 欧美激情国产日韩精品一区| 床上黄色一级片| 国产乱人偷精品视频| 亚洲欧美一区二区三区黑人 | 日本av手机在线免费观看| 中文字幕久久专区| 亚洲精品乱久久久久久| 最近手机中文字幕大全| 亚洲经典国产精华液单| 国产亚洲午夜精品一区二区久久 | 中文在线观看免费www的网站| 久久精品夜夜夜夜夜久久蜜豆| 欧美 日韩 精品 国产| freevideosex欧美| 简卡轻食公司| 全区人妻精品视频| 一区二区三区免费毛片| 精品少妇黑人巨大在线播放| 欧美性感艳星| 国产精品麻豆人妻色哟哟久久 | 极品少妇高潮喷水抽搐| 国内揄拍国产精品人妻在线| 少妇熟女欧美另类| 中文天堂在线官网| 国产老妇伦熟女老妇高清| 夜夜看夜夜爽夜夜摸| 一本久久精品| 免费av不卡在线播放| 久久精品久久精品一区二区三区| 日韩一区二区视频免费看| 亚洲精品乱久久久久久| 精品人妻一区二区三区麻豆| 欧美bdsm另类| 中文字幕制服av| 亚洲熟女精品中文字幕| 欧美3d第一页| 亚洲国产精品sss在线观看| 久久99热6这里只有精品| 国内少妇人妻偷人精品xxx网站| 好男人在线观看高清免费视频| 九九在线视频观看精品| 久久午夜福利片| 久久精品国产亚洲av涩爱| 国产一区二区亚洲精品在线观看| 亚洲在线观看片| 99久久精品热视频| 精品久久久久久成人av| 免费播放大片免费观看视频在线观看| 晚上一个人看的免费电影| 久久99蜜桃精品久久| 国产亚洲一区二区精品| 久久久午夜欧美精品| 国产精品嫩草影院av在线观看| 精品酒店卫生间| 亚洲国产精品成人综合色| 亚洲无线观看免费| 国产国拍精品亚洲av在线观看| 亚洲在线自拍视频| 亚洲av.av天堂| 亚洲不卡免费看| 免费av毛片视频| 日韩欧美国产在线观看| 免费观看a级毛片全部| 久久99热6这里只有精品| 国产一区二区三区av在线| 亚洲精品456在线播放app| 国产一区二区在线观看日韩| 亚洲精品国产av蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av在线观看美女高潮| 哪个播放器可以免费观看大片| 日本黄色片子视频| 少妇人妻一区二区三区视频| 欧美成人精品欧美一级黄| 男女下面进入的视频免费午夜| 国产单亲对白刺激| 街头女战士在线观看网站| 亚洲成人久久爱视频| 国产精品一区二区在线观看99 | 韩国高清视频一区二区三区| 深夜a级毛片| 亚洲成人一二三区av| 国产 一区精品| 亚洲国产成人一精品久久久| 26uuu在线亚洲综合色| 极品教师在线视频| 亚洲精品日韩av片在线观看| 18禁动态无遮挡网站| 精华霜和精华液先用哪个| 晚上一个人看的免费电影| 国产伦精品一区二区三区视频9| 直男gayav资源| 一级毛片 在线播放| 国产探花在线观看一区二区| 夫妻午夜视频| 麻豆成人av视频| 校园人妻丝袜中文字幕| 99久久九九国产精品国产免费| 午夜老司机福利剧场| 大陆偷拍与自拍| 成年人午夜在线观看视频 | 九色成人免费人妻av| 看黄色毛片网站| 偷拍熟女少妇极品色| 久久久久国产网址| 精品久久久久久久末码| 九九在线视频观看精品| 老司机影院成人| 午夜日本视频在线| 韩国高清视频一区二区三区| 久久精品国产自在天天线| 亚洲成人一二三区av| 久久精品国产自在天天线| 极品少妇高潮喷水抽搐| 国产亚洲91精品色在线| 欧美高清成人免费视频www| 激情 狠狠 欧美| 久久精品国产亚洲av涩爱| 乱人视频在线观看| 成人一区二区视频在线观看| 精品午夜福利在线看| 亚洲三级黄色毛片| 精品一区在线观看国产| 亚洲欧美日韩东京热| 99九九线精品视频在线观看视频| 99re6热这里在线精品视频| 免费观看性生交大片5| 三级男女做爰猛烈吃奶摸视频| 午夜福利在线观看吧| 身体一侧抽搐| 一区二区三区乱码不卡18| 一个人看的www免费观看视频| 亚洲av免费高清在线观看| 国产极品天堂在线| 亚洲av免费高清在线观看| 国产精品国产三级专区第一集| 国产91av在线免费观看| 国产亚洲5aaaaa淫片| 日本三级黄在线观看| 国产国拍精品亚洲av在线观看| 尾随美女入室| 女人十人毛片免费观看3o分钟| 亚洲丝袜综合中文字幕| 日韩成人av中文字幕在线观看|