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

    中國西門塔爾牛肉用群體重要經(jīng)濟(jì)性狀遺傳參數(shù)估計(jì)

    2016-11-01 11:30:38寶金山張文剛夏江威宋禹昕徐凌洋張路培高會(huì)江李俊雅
    畜牧獸醫(yī)學(xué)報(bào) 2016年9期

    牛 紅,寶金山,吳 洋,朱 波,張文剛,夏江威,宋禹昕,郭 鵬,徐凌洋,陳 燕,高 雪,張路培,高會(huì)江,李俊雅*

    (1.中國農(nóng)業(yè)科學(xué)院北京畜牧獸醫(yī)研究所,北京 100193; 2.錫林郭勒盟烏拉蓋管理區(qū)獸醫(yī)局,烏拉蓋 026321)

    ?

    中國西門塔爾牛肉用群體重要經(jīng)濟(jì)性狀遺傳參數(shù)估計(jì)

    牛紅1,寶金山2,吳洋1,朱波1,張文剛1,夏江威1,宋禹昕1,郭鵬1,徐凌洋1,陳燕1,高雪1,張路培1,高會(huì)江1,李俊雅1*

    (1.中國農(nóng)業(yè)科學(xué)院北京畜牧獸醫(yī)研究所,北京 100193; 2.錫林郭勒盟烏拉蓋管理區(qū)獸醫(yī)局,烏拉蓋 026321)

    本研究旨在對出生于2008-2013年的中國西門塔爾牛肉用群體的重要經(jīng)濟(jì)性狀進(jìn)行遺傳參數(shù)和方差組分估計(jì),該群體大小為2 939頭。采用非求導(dǎo)約束最大似然法估計(jì)遺傳力、遺傳相關(guān)和方差組分。結(jié)果顯示,初生重、斷奶重、出欄重、胴體重、屠宰率和凈肉率的遺傳力估計(jì)值分別為0.48、0.44、0.43、0.38、0.31、0.39。其中,初生重和斷奶重的遺傳相關(guān)估計(jì)結(jié)果為0.57,出欄重和胴體重的遺傳相關(guān)為0.94,屠宰率和凈肉率的遺傳相關(guān)為0.89。該群體的生長和屠宰相關(guān)性狀均屬于中高遺傳力性狀,且性狀之間具有較高的遺傳相關(guān)。本研究對中國西門塔爾牛肉用群體重要經(jīng)濟(jì)性狀的遺傳參數(shù)做了系統(tǒng)評估分析,為將來制定育種方案和遺傳評估奠定基礎(chǔ)。

    中國西門塔爾牛肉用群體;非求導(dǎo)約束最大似然法;方差組分;遺傳參數(shù)

    家畜育種的主要目的是不斷改進(jìn)畜群遺傳品質(zhì),提高個(gè)體生產(chǎn)性能,培育新的品種或品系。然而,一個(gè)育種方案的制定必須以準(zhǔn)確、可靠的遺傳參數(shù)為前提[1-2]。同時(shí),家畜重要經(jīng)濟(jì)性狀的遺傳評估也為確定選種方法,預(yù)測選種效果,制定選擇指數(shù)以及進(jìn)行間接選種提供重要依據(jù)[2]。中國西門塔爾牛于2001通過國家畜禽品種資源委員會(huì)牛品種審定委員會(huì)的審定,成為中國牛培育品種中的一個(gè)新品種。它適應(yīng)性好,抗病力強(qiáng),耐粗飼,分布范圍廣,在我國多種生態(tài)條件下均能表現(xiàn)出良好的生產(chǎn)性能[3]。由于其良好的生產(chǎn)性能,育種工作者對中國西門塔爾牛進(jìn)行大量的基礎(chǔ)性研究工作。汪春乾等[4]估計(jì)中國西門塔爾牛生長性狀(出生重、斷奶重、周歲重、18月齡重和成年體重)的遺傳參數(shù)為0.11~0.47,周歲重與18月齡重、成年母牛體重間的遺傳相關(guān)分別為0.84和0.64。范大有等[5]估計(jì)了中國西門塔爾牛次級(jí)性狀和生產(chǎn)性狀共19個(gè)性狀的遺傳參數(shù)為0.07~0.41,其中泌乳速度與使用年限和整體評分有較高的遺傳相關(guān),分別為0.87和0.61。近年來,隨著冷配技術(shù)的不斷實(shí)施,基礎(chǔ)母牛數(shù)量的不斷增加,中國西門塔爾牛肉用群體已具備一定規(guī)模。

    為提高育種工作效率,加快遺傳進(jìn)展,以及為今后制定具體育種方案提供依據(jù),本研究應(yīng)用非求導(dǎo)約束最大似然法[6-7]對內(nèi)蒙古烏拉蓋管理區(qū)的中國西門塔爾牛肉用群體的重要經(jīng)濟(jì)性狀進(jìn)行了遺傳參數(shù)估計(jì)。

    1 材料與方法

    1.1材料

    試驗(yàn)動(dòng)物為出生于內(nèi)蒙古錫林郭勒盟烏拉蓋管理區(qū)的中國西門塔爾牛肉用群體,該群體自2008年開始建立并記錄數(shù)據(jù)。初生重和斷奶重為實(shí)時(shí)測定,其中初生重?cái)?shù)據(jù)記錄2 939頭,斷奶重?cái)?shù)據(jù)記錄2 910頭。為了便于其他數(shù)據(jù)的測量和采集,將5~9月齡的西門塔爾牛運(yùn)送至北京金維福仁清真食品有限公司肉牛養(yǎng)殖場,在該場期間按照統(tǒng)一的飼養(yǎng)管理方法集中育肥,每3個(gè)月測量一次體尺和體重?cái)?shù)據(jù)。當(dāng)飼養(yǎng)至18~24月齡時(shí)進(jìn)行分批屠宰分割,根據(jù)中華人民共和國國家標(biāo)準(zhǔn)(GB/T 27643-2011)測量收集出欄重、胴體重、屠宰率和凈肉率性狀數(shù)據(jù)。從2008年到2013年共收集1 302頭個(gè)體表型數(shù)據(jù),其中有個(gè)體系譜記錄的為1 086頭。

    1.2影響因素統(tǒng)計(jì)

    影響因素的校正可以降低數(shù)據(jù)資料的不平衡性,提高估計(jì)的準(zhǔn)確性。本研究中考慮的影響因素為母親年齡、性別、出生年和季節(jié)效應(yīng)、出生場效應(yīng),除此之外,出欄重、胴體重、屠宰率和凈肉率還把進(jìn)場重和育肥天數(shù)作為協(xié)變量。試驗(yàn)動(dòng)物群體中母親年齡在9歲和9歲以上數(shù)量較少,因此,將9歲和9歲以上的個(gè)體當(dāng)作同一水平。根據(jù)烏拉蓋地區(qū)的牧戶所在區(qū)域的分布,將出生場效應(yīng)劃分為5個(gè)場,分別命名為1、2、3、4、5場。

    表1為1 302頭西門塔爾牛的出欄重、胴體重、屠宰率和凈肉率性狀的各影響因素不同水平的觀察個(gè)數(shù)。而初生重和斷奶重與其他4個(gè)性狀的影響因素一樣,只是不同水平的觀察個(gè)數(shù)不同。

    1.3影響因素顯著性檢驗(yàn)

    影響因素的處理和顯著性檢驗(yàn)是通過R程序中的GLM進(jìn)行分析,采用的模型是一般線性模型(General linear model)。如下所示:yijklmn=u+dai+sexj+yeark+seasonl+farmm+enterweight+fattenday+eijklmn

    (1)

    其中,yijklmn是性狀表型值,u為群體均值,dai是母親年齡效應(yīng),sexj是性別效應(yīng),yeark是犢牛出生年效應(yīng),seasonl是犢牛出生季節(jié)效應(yīng),farmm是犢牛出生場效應(yīng),enterweight是進(jìn)場重協(xié)變量,fattenday是育肥天數(shù)協(xié)變量,eijklmn是隨機(jī)殘差。在初生重和斷奶重影響因素顯著性檢驗(yàn)分析模型中,不包括進(jìn)場重和育肥天數(shù)協(xié)變量,而出欄重、胴體重、屠宰率和凈肉率影響因素的顯著性檢驗(yàn)分析模型中則包括進(jìn)場重和育肥天數(shù)協(xié)變量??傊瑢τ诒硇托誀畹男U?,要把顯著的影響因素加入到校正模型中,而不顯著的影響因素要去掉。

    1.4遺傳參數(shù)估計(jì)

    遺傳參數(shù)估計(jì)模型:

    y=1nu+Zg+e

    (2)

    其中,y為校正的表型向量;1n是n個(gè)值為1的列向量;u為均值;Z為隨機(jī)效應(yīng)育種值的結(jié)構(gòu)矩陣;g為個(gè)體育種值向量;e為隨機(jī)誤差向量。該模型的期望和方差:

    (3)

    (4)

    其中,E(y),E(g)和E(e)分別為表型向量、個(gè)體育種值向量和隨機(jī)誤差向量的期望,V(g),V(e)和V(y)分別為個(gè)體育種值向量、隨機(jī)誤差向量和表型向量的方差,G為遺傳方差結(jié)構(gòu)矩陣,R為隨機(jī)誤差方差結(jié)構(gòu)矩陣,Z′為矩陣Z的轉(zhuǎn)置矩陣。

    表1出欄重、胴體重、屠宰率和凈肉率性狀影響因素不同水平的觀察個(gè)數(shù)

    Table 1The number of observations at different levels of influencing factors in marketing weight,carcass weight,dressing percentage and meat percentage

    出生年Yearofbirth出生季Seasonofbirth出生場Farmofbirth母親年齡Ageofdam性別Sex水平Level觀察值個(gè)數(shù)N水平Level觀察值個(gè)數(shù)N水平Level觀察值個(gè)數(shù)N水平Level觀察值個(gè)數(shù)N水平Level觀察值個(gè)數(shù)N200813513051135255公11962009286229626553165母10620102433346318841642011296435541035356201225552216236201387710081589~1268

    其中出生年單位為年;出生季1、2、3、4指代春季、夏季、秋季、冬季;母親年齡單位為歲

    The unit of born year is the number of year;The unit of born season 1,2,3,4 means spring,summer,autumn and winter,respectively;The unit of age of dam is the age

    方程組:

    (5)

    (6)

    (7)

    (8)

    其中,h2為遺傳力,rA為遺傳相關(guān),rP為表型相關(guān),a1,a2為兩性狀的育種值,p1,p2為兩性狀的表型值。

    2 結(jié) 果

    2.1主要性狀的基本統(tǒng)計(jì)量

    表2列出了各性狀的均值、標(biāo)準(zhǔn)差和變異系數(shù),從表2中可以看出,初生重、斷奶重、出欄重和胴體重的均值分別為39.98、183.22、508.73和274.02,該群體的屠宰率和凈肉率平均值為53.74%和45.72%。這些與體重和屠宰性狀相關(guān)的群體平均值相對于國內(nèi)外的其他品種均處于較高的水平。但體重相關(guān)性狀的變異系數(shù)為13.74%~24.59%,顯示出這些性狀在不同的場、年、季、性別等不同環(huán)境效應(yīng)條件下有較大的波動(dòng)。而屠宰率和凈肉率的變異系數(shù)均在10%以下,表現(xiàn)出相對的穩(wěn)定性。圖1為各個(gè)性狀表型值每年均值的變化圖。

    2.2影響因素的顯著性結(jié)果

    如表3中所示,通過顯著性檢驗(yàn)可以得知,母親年齡效應(yīng)對初生重和斷奶重性狀影響極顯著(P<0.01),對出欄重影響顯著(P<0.05),而對胴體重、屠宰率和凈肉率影響不顯著。性別效應(yīng)對出欄重性狀影響極顯著(P<0.01),對其他性狀影響均為顯著(P<0.05)。出生場效應(yīng)、年效應(yīng)和季節(jié)效應(yīng)對于不同性狀有不同的影響。進(jìn)場重對出欄重性狀影響為極顯著,對胴體重、屠宰率和凈肉率性狀影響均為顯著,而育肥天數(shù)對出欄重、胴體重、屠宰率和凈肉率性狀影響均為極顯著。因此,在校正固定效應(yīng)的時(shí)候,把顯著的因素加入到模型進(jìn)行校正,而不顯著的因素去掉。

    圖1 重要經(jīng)濟(jì)性狀不同年份的平均表型值變化圖Fig.1 The mean value of important economic traits for each year

    表2重要經(jīng)濟(jì)性狀的表型基本統(tǒng)計(jì)量

    Table 2Summary of important economic traits

    性狀Trait記錄個(gè)數(shù)Number均值Mean標(biāo)準(zhǔn)差Standarderror變異系數(shù)Variablecoefficient初生重/kgBirthweight293939.989.8324.59斷奶重/kgWeaningweight2910183.2238.8824.49出欄重/kgMarketingweight1302508.7369.9113.74胴體重/kgCarcassweight1302274.0244.8916.38屠宰率/%Dressingpercentage130253.742.845.29凈肉率/%Meatpercentage130245.722.946.44

    利用鄧肯多重比較研究固定效應(yīng)中的兩個(gè)影響因素(場效應(yīng)和年效應(yīng))對初生重、斷奶重、出欄重、胴體重、屠宰率和凈肉率的影響,見表4,可以看出,除初生重和出欄重不隨場的變化而變化外(P>0.05),斷奶重、胴體重、屠宰率和凈肉率4個(gè)性狀均呈現(xiàn)隨場效應(yīng)顯著變化的趨勢(P<0.05),而6個(gè)性狀均呈現(xiàn)隨年效應(yīng)顯著變化的趨勢(P<0.05)。

    2.3主要性狀的遺傳參數(shù)估計(jì)

    各性狀的方差組分估計(jì)結(jié)果見表5。將表中的方差組分估計(jì)值帶入相應(yīng)的遺傳力計(jì)算公式,便可以得出遺傳力h2的估計(jì)結(jié)果。從表中可以看出,中國西門塔爾牛肉用群體的生長和屠宰相關(guān)性狀的遺傳力均屬于中高遺傳力(0.31~0.48)。

    初生重、斷奶重、出欄重、胴體重、屠宰率和凈肉率的遺傳相關(guān)和表型相關(guān)的估計(jì)結(jié)果見表6。從表中我們可以看出初生重和斷奶重、出欄重和胴體重、屠宰率和凈肉率之間均呈現(xiàn)出較高的遺傳相關(guān)和表型相關(guān)。

    表3不同性狀影響因素的顯著性檢驗(yàn)結(jié)果(P值)

    Table 3Significance test results of factors affecting different traits

    性狀Trait母親年齡Dam性別Sex場Farm年Year季節(jié)Season進(jìn)場重Enterweight育肥天數(shù)Fattenday初生重Birthweight0.003**0.018*0.0520.018*0.026*--斷奶重Weaningweight0.001**0.024*0.003**0.043*0.014*--出欄重Marketingweight0.044*0.007**0.0580.006**0.027*0.002**0.001**胴體重Carcassweight0.0580.012*0.023*0.005**0.044*0.017*0.004**屠宰率Dressingpercentage0.0680.041*0.016*0.014*0.035*0.026*0.002**凈肉率Meatpercentage0.0720.038*0.028*0.026*0.039*0.014*0.005**

    *.P<0.05;**.P<0.01

    表4不同性狀主要影響因素(場效應(yīng)和年效應(yīng))的多重比較

    Table 4Multiple comparison of the major influence factors (farm and year) for different traits

    因素Factor初生重/kgBirthweight斷奶重/kgWeaningweight出欄重/kgMarketingweight胴體重/kgCarcassweight屠宰率/%Dressingpercentage凈肉率/%Meatpercentage場效應(yīng)Farm137.28a174.84a505.33a316.79a56.26a47.49a237.88a188.35ac501.54a268.71b53.47b45.51b337.54a212.84b509.86a257.06b52.55c44.60c437.65a212.50b508.09a296.16c54.91d46.94ad537.72a197.35bc504.63a288.33c54.61d46.51d年效應(yīng)Year200834.01a153.32a460.38a231.71a50.31a41.23a200937.64b207.57b508.25b269.65b52.26b44.05b201037.02b155.76a465.94a243.76a53.02b45.32c201139.57c209.32b496.76b270.83b54.43c46.17c201239.22c201.39b547.46c302.62c54.61c46.33c201340.51d205.75b578.49d315.93d55.15c47.43d

    各影響因素中同一列數(shù)值右肩的不同大寫字母表示差異顯著(P<0.05)

    The letters on numerical right shoulder in the same line for different influence factors indicate significant difference (P<0.05)

    表5重要經(jīng)濟(jì)性狀的遺傳參數(shù)估計(jì)

    Table 5Genetic parameters of important economic traits

    性狀Trait表型方差σ2p加性方差σ2a殘差方差σ2e遺傳h2初生重Birthweight73.4335.4537.980.48斷奶重Weaningweight1485.97646.47839.490.44出欄重Marketingweight1895.36806.661088.700.43胴體重Carcassweight709.38269.37440.010.38屠宰率Dressingpercentage5.921.814.210.31凈肉率Meatpercentage6.232.423.810.39

    表6重要經(jīng)濟(jì)性狀兩兩間的遺傳相關(guān)估計(jì)

    Table 6Genetic correlation among important economic traits

    性狀Trait初生重Birthweight斷奶重Weaningweight出欄重Marketingweight胴體重Carcassweight屠宰率Dressingpercentage凈肉率Meatpercentage初生重Birthweight10.570.040.040.010.03斷奶重Weaningweight0.5510.130.11-0.04-0.06出欄重Marketingweight0.260.5210.940.230.20胴體重Carcassweight0.240.500.9210.520.45屠宰率Dressingpercentage0.090.230.380.6210.89凈肉率Meatpercentage0.080.280.440.640.751

    其中上三角為遺傳相關(guān),下三角為表型相關(guān)

    Genetic correlations are given above the diagonal and phenotypic correlations below the diagonal

    3 討 論

    西門塔爾牛以其優(yōu)良的乳肉特性和良好的適應(yīng)性逐漸成為中國地方黃牛雜交改良的首選品種。近年來,隨著人們生活水平的提高,對牛肉的需求量正在逐年增加[3,9]。由于西門塔爾牛具有生長速度快、耐粗飼、繁殖力強(qiáng)等特性,這使得中國西門塔爾牛肉用群體的建立顯得十分必要。人們可以立足已建立的中國西門塔爾?;蚪M選擇肉用參考群體的基礎(chǔ)上,完善肉牛的育種技術(shù)和建立全基因組選擇技術(shù)平臺(tái),開展肉牛全基因組選擇、新基因發(fā)掘和功能解析等方面的研究[10-11],提供使肉牛種業(yè)可持續(xù)發(fā)展的技術(shù)支撐平臺(tái),并培育我國自主的肉用西門塔爾牛新品系,擺脫我國肉牛種業(yè)依賴進(jìn)口的局面,為全國肉牛進(jìn)一步遺傳改良提供借鑒和參考。

    遺傳參數(shù)的估計(jì)因世代、畜群結(jié)構(gòu)、畜群飼養(yǎng)管理及參數(shù)估計(jì)時(shí)所采用方法的不同,結(jié)果也有所差異。本研究采用非求導(dǎo)約束最大似然法估計(jì)了中國西門塔爾牛肉用群體的生長和屠宰性狀的遺傳力,均屬于中高遺傳力性狀,所估計(jì)的遺傳力均處于文獻(xiàn)報(bào)道的范圍之內(nèi)[12-13]。對于初生重遺傳力估計(jì)值,K.D.Bullock等[14]估計(jì)的海福特牛的遺傳力為0.49,R.Alenda[15]估計(jì)安格斯牛的遺傳力為0.46,M.G.Jeyaruban等[16]估計(jì)澳大利亞西門塔爾牛的遺傳力為0.36。國內(nèi)對于西門塔爾牛群體初生重的遺傳力估計(jì)值為0.32[4]。斷奶重則反映肉牛早期生長狀況,Winder[17]估計(jì)安格斯牛斷奶重的遺傳力為0.39,L.L.Benyshek[18]估計(jì)的西門塔爾牛群體的斷奶重遺傳力為0.34,D.J.Garrick等[19]運(yùn)用多性狀模型得出西門塔爾牛公犢、母犢遺傳力分別為0.32和0.39,均與本研究相近。在出欄重和胴體重遺傳參數(shù)的報(bào)道中,R.L.Baker[20]估計(jì)海福特群體的成年牛體重遺傳力為0.38,F(xiàn).Mukai等[21]在日本和牛群體中估計(jì)的胴體重遺傳力為0.39。屠宰率和凈肉率的相關(guān)報(bào)道中,L.L.Benyshek[13]和C.A.Morris等[22]估計(jì)群體的遺傳力均為0.31。以上報(bào)道均說明本研究對于該群體的遺傳參數(shù)估計(jì)結(jié)果與國內(nèi)外報(bào)道較為一致。

    本研究估計(jì)的西門塔爾牛的初生重和斷奶重遺傳相關(guān)為0.57,小于M.B.Long等[23]報(bào)道的海福特牛初生重和斷奶重的遺傳相關(guān)值0.64。同時(shí)本研究估計(jì)的出欄重和胴體重的遺傳相關(guān)為0.94,與D.M.Marshall[24]報(bào)道的海福特牛出欄重和胴體重遺傳相關(guān)0.93相一致。而屠宰率和凈肉率的遺傳相關(guān),國外未見相關(guān)報(bào)道。本研究中所估計(jì)的遺傳參數(shù)與國外相比,其差異可能與數(shù)據(jù)量和估計(jì)方法有關(guān)。在育種工作中,若兩個(gè)性狀有較高的遺傳相關(guān),選擇其中一個(gè)性狀的同時(shí)會(huì)間接地改進(jìn)和提高另一性狀的選擇效果。因此,遺傳參數(shù)的估計(jì)在肉牛生產(chǎn)中具有較高實(shí)用價(jià)值,對進(jìn)行下一步的基因組選擇研究和肉牛新品種或新品系的培育有重要意義。

    4 結(jié) 論

    中國西門塔爾牛肉用群體初生重和斷奶重的遺傳力分別為0.48和0.44,出欄重和胴體重的遺傳力分別為0.43和0.38,屠宰率和凈肉率的遺傳力分別為0.31和0.39,這些性狀的遺傳力均為中等偏上水平。初生重和斷奶重、出欄重和胴體重、屠宰率和凈肉率之間遺傳相關(guān)分別為0.57、0.94和0.89。本研究首次對中國西門塔爾牛肉用群體重要經(jīng)濟(jì)性狀的遺傳參數(shù)做了系統(tǒng)評估分析,為將來制定育種方案和遺傳評估奠定基礎(chǔ)。

    [1]AKESSON M,BENSCH S,HASSELQUIST D,et al.Estimating heritabilities and genetic correlations:comparing the ‘a(chǎn)nimal model’ with parent-offspring regression using data from a natural population[J].PLoSOne,2008,3(3):e1739.

    [2]KRUUK L E.Estimating genetic parameters in natural populations using the "animal model"[J].PhilosTransRSocLondBBiolSci,2004,359(1446):873-900.

    [3]方曉敏,許尚忠,張英漢.我國新的牛種資源——中國西門塔爾牛[J].黃牛雜志,2002,28(5):67-69.

    FANG X M,XU S Z,ZHANG Y H.A new cattle resources in our country-Chinese Simmental cattle[J].JournalofYellowCattleScience,2002,28(5):67-69.(in Chinese)

    [4]汪春乾,許尚忠,陳宏權(quán),等.中國西門塔爾牛生長性狀遺傳參數(shù)估計(jì)[J].安徽農(nóng)業(yè)大學(xué)學(xué)報(bào),2004,31(4):417-420.

    WANG C Q,XU S Z,CHEN H Q,et al.Estimating parameters growth traits of Chinese Simmental[J].JournalofAnhuiAgriculturalUniversity,2004,31(4):417-420.(in Chinese)

    [5]范大有,許尚忠,李俊雅,等.中國西門塔爾牛次級(jí)性狀與生產(chǎn)性狀的遺傳統(tǒng)計(jì)分析[J].畜牧獸醫(yī)學(xué)報(bào),2008,39(8):1025-1032.FANG D Y,XU S Z,LI J Y,et al.Genetic and statistical analysis between production traits and secondary traits in Chinese Simmental[J].ActaVeterinariaetZootechnicaSinica,2008,39(8):1025-1032.(in Chinese)

    [6]HENDERSON C R,QUAAS R L.Multiple trait evaluation using relatives records[J].JAnimSci,1976,43(6):1188-1197.

    [7]ZHANG L C,NING Z H,XU G Y,et al.Heritabilities and genetic and phenotypic correlations of egg quality traits in brown-egg dwarf layers[J].PoultSci,2005,84(8):1209-1213.

    [8]VAN VLECK L D,BOLDMAN K.A manual for use of MTDFREML[M].USDA:ARS,1993.

    [9]霍靈光,田露,張?jiān)浇?中國牛肉需求量中長期預(yù)測分析[J].中國畜牧雜志,2010,46(2):43-47.HUO L G,TIAN L,ZHANG Y J.Chinese beef demand in the long-term prediction analysis[J].ChineseJournalofAnimalScience,2010,46(2):43-47.(in Chinese)

    [10]WU Y,F(xiàn)AN H Z,WANG Y H,et al.Genome-wide association studies using haplotypes and individual SNPs in Simmental cattle[J].PLoSOne,2014,9(10):e109330.

    [11]ZHU M,ZHU B,WANG Y H,et al.Linkage disequilibrium estimation of chinese beef simmental cattle using high-density SNP panels[J].Asian-AustralasJAnimSci,2013,26(6):772-779.

    [12]KAUSE A,MIKKOLA L,STRANDEN I,et al.Genetic parameters for carcass weight,conformation and fat in five beef cattle breeds[J].Animal,2015,9(1):35-42.

    [13]BENYSHEK L L.Heritabilities for growth and carcass traits estimated from data on Hereford under commercial conditions[J].JAnimSci,1981,53(1):49-56.

    [14]BULLOCK K D,BERTRAND J K,BENYSHEK L L.Genetic and environmental parameters for mature weight and other growth measures in Polled Hereford cattle[J].JAnimSci,1993,71(7):1737-1741.

    [15]ALENDA R.Genetic parameters and consequences of selection for growth traits in a beef herd selected for yearling weight[J].JAnimSci,1987,64(2):366-372.

    [16]JEYARUBAN M G,JOHNSTON D J,TIER B,et al.Genetic parameters for calving difficulty using complex genetic models in five beef breeds in Australia[J].AnimProdSci,2016,56(5):927-933.

    [17]WINDER.Relationships among factors associated with mothering ability in beef cattle[J].JAnimSci,1990,18(3):938-946.

    [18]BENYSHEK L L.Estimating heritability of Hereford cow weight[J].JAnimSci,1973,36(1):854-861.

    [19]GARRICK D J,POLLAK E J,QUAAS R L,et al.Variance heterogeneity in direct and maternal weight traits by sex and percent purebred for Simmental-sired calves[J].JAnimSci,1989,67(10):2515-2528.

    [20]BAKER R L.Results of selection for yearling or 18-month weight in Angus and Hereford cattle[J].LivestProdSci,1991,29(4):277-296.

    [21]MUKAI F,OYAMA K,KOHNO S.Genetic relationships between performance test traits and field carcass traits in Japanese Black cattle[J].LivestProdSci,1995,44(3):199-205.

    [22]MORRIS C A,CULLEN N G,MCCALL D G.Genetic and phenotypic relationships among carcass measurements in beef cattle[J].NZJAgrRes,1999,42(4):415-421.

    [23]LONG M B,MARSHAL D M.Relationship of Beef sire birth weight and weaning weight expected progeny differences to actual performance of crossbred offspring[J].Cattle,1994,10:31-35.

    [24]MARSHALL D M.Genetic parameters for carcass traits in beef cattle[J].Cattle,1993,8:28-34.

    (編輯郭云雁)

    Estimation of Genetic Parameters for Economic Important Traits in Chinese Simmental Beef Cattle

    NIU Hong1,BAO Jin-shan2,WU Yang1,ZHU Bo1,ZHANG Wen-gang1,XIA Jiang-wei1,SONG Yu-xin1,GUO Peng1,XU Ling-yang1,CHEN Yan1,GAO Xue1,ZHANG Lu-pei1,GAO Hui-jiang1,LI Jun-ya1*

    (1.InstituteofAnimalScience,ChineseAcademyofAgriculturalSciences,Beijing100193,China;2.VeterinaryBureauofUlgaiPrecinctinXilinGolLeague,Ulgai026321,China)

    This experiment was conducted to estimate the variance component and genetic parameters including heritability and genetic correlation for economic important traits of 2 939 Simmental beef cattle that were born between 2008 and 2013.These parameters were estimated using derivative-free restricted maximum likelihood method.We found that the estimation of heritability for birth weight,weaning weight,marketing weight,carcass weight,dressing percentage and meat percentage were 0.48,0.44,0.43,0.38,0.31 and 0.39,respectively.Furthermore,the genetic correlation between birth weight and weaning weight,marketing weight and carcass weight,dressing percentage and meat percentage were 0.57,0.94 and 0.89,respectively.The heritabilities of these development traits and carcass traits in Chinese Simmental beef cattle population belong to medium and high heritability,and these traits had high genetic correlation values.Our study conducted the system analysis on genetic parameters of important economic traits in Chinese Simmental beef cattle.And the results lay the foundation for further breeding scheme and genetic evaluation.

    Chinese Simmental beef cattle population;derivative-free restricted maximum likelihood;variance component;genetic parameters

    10.11843/j.issn.0366-6964.2016.09.009

    2016-01-22

    “十二五”國家科技支撐計(jì)劃(2011BAD28B04);國家自然科學(xué)基金(31372294);國家高技術(shù)研究發(fā)展計(jì)劃(2013AA102505-4)

    牛紅(1990-),女,吉林公主嶺人,碩士生,主要從事分子數(shù)量遺傳學(xué)研究,Tel: 010-62818176,E-mail: nhnh0304@126.com

    李俊雅,研究員,E-mail: JL1@iascaas.net.cn

    S823;S813.3

    A

    0366-6964(2016)09-1817-07

    99久久综合免费| 欧美日韩综合久久久久久| 黄色怎么调成土黄色| 大码成人一级视频| 日本与韩国留学比较| 亚洲欧美精品自产自拍| 人人妻人人爽人人添夜夜欢视频| 午夜91福利影院| 日本爱情动作片www.在线观看| 亚洲综合色网址| 日日爽夜夜爽网站| 国产片特级美女逼逼视频| 免费黄色在线免费观看| 午夜精品国产一区二区电影| 另类精品久久| 久久国产精品男人的天堂亚洲 | 美女cb高潮喷水在线观看| 亚洲av免费高清在线观看| 九色亚洲精品在线播放| 国产精品蜜桃在线观看| 又黄又爽又刺激的免费视频.| 另类亚洲欧美激情| 精品国产露脸久久av麻豆| 亚洲一区二区三区欧美精品| 一本—道久久a久久精品蜜桃钙片| 国产精品三级大全| 美女福利国产在线| 欧美丝袜亚洲另类| 80岁老熟妇乱子伦牲交| 日产精品乱码卡一卡2卡三| 日本wwww免费看| 亚洲情色 制服丝袜| 成年人午夜在线观看视频| 亚洲五月色婷婷综合| 性色avwww在线观看| av卡一久久| 日日爽夜夜爽网站| 18禁在线无遮挡免费观看视频| 久久韩国三级中文字幕| 卡戴珊不雅视频在线播放| 亚洲欧美清纯卡通| 男人爽女人下面视频在线观看| 欧美xxxx性猛交bbbb| 一边摸一边做爽爽视频免费| 能在线免费看毛片的网站| a级毛片黄视频| 日韩一区二区三区影片| 在线亚洲精品国产二区图片欧美 | 人人妻人人澡人人看| 欧美少妇被猛烈插入视频| 王馨瑶露胸无遮挡在线观看| 欧美成人午夜免费资源| 欧美日韩一区二区视频在线观看视频在线| 日日撸夜夜添| 丰满迷人的少妇在线观看| 国产综合精华液| 亚洲情色 制服丝袜| 欧美少妇被猛烈插入视频| 久久国产亚洲av麻豆专区| 五月天丁香电影| 国产精品久久久久久精品古装| 日本91视频免费播放| 91精品一卡2卡3卡4卡| 成人亚洲欧美一区二区av| 国产精品 国内视频| 国产欧美日韩综合在线一区二区| 香蕉精品网在线| 超碰97精品在线观看| 麻豆精品久久久久久蜜桃| 国产精品国产三级国产专区5o| 三级国产精品片| 建设人人有责人人尽责人人享有的| 狂野欧美白嫩少妇大欣赏| 少妇熟女欧美另类| 一级黄片播放器| 亚洲人与动物交配视频| 久久精品国产亚洲网站| 一级黄片播放器| 久久人人爽人人片av| 国产精品一二三区在线看| 秋霞伦理黄片| 最近的中文字幕免费完整| 国产精品一区www在线观看| 久久韩国三级中文字幕| 久久综合国产亚洲精品| 日韩av免费高清视频| 91久久精品国产一区二区三区| 精品久久国产蜜桃| 女人精品久久久久毛片| 免费播放大片免费观看视频在线观看| 2022亚洲国产成人精品| 久久人妻熟女aⅴ| 亚洲国产精品一区二区三区在线| 男女免费视频国产| .国产精品久久| 欧美精品一区二区免费开放| 国产片内射在线| 夜夜骑夜夜射夜夜干| 久久久久久久国产电影| 全区人妻精品视频| 菩萨蛮人人尽说江南好唐韦庄| 日韩成人av中文字幕在线观看| a 毛片基地| 午夜免费男女啪啪视频观看| 欧美三级亚洲精品| 各种免费的搞黄视频| 久久ye,这里只有精品| 插阴视频在线观看视频| 亚洲精品一区蜜桃| 欧美最新免费一区二区三区| 一级片'在线观看视频| 国产成人a∨麻豆精品| 久久久久国产精品人妻一区二区| 美女视频免费永久观看网站| 内地一区二区视频在线| 在线观看www视频免费| 国产成人aa在线观看| 青春草亚洲视频在线观看| 多毛熟女@视频| 午夜激情久久久久久久| 国国产精品蜜臀av免费| 多毛熟女@视频| 精品一区在线观看国产| 青春草亚洲视频在线观看| 99九九线精品视频在线观看视频| 青春草亚洲视频在线观看| 丰满迷人的少妇在线观看| 欧美日韩视频高清一区二区三区二| 啦啦啦中文免费视频观看日本| 国产成人免费无遮挡视频| 国产精品久久久久成人av| av线在线观看网站| 国产白丝娇喘喷水9色精品| 日产精品乱码卡一卡2卡三| 日本猛色少妇xxxxx猛交久久| 精品国产国语对白av| 久久99热6这里只有精品| 精品一区二区三区视频在线| 少妇的逼好多水| 成人手机av| 国产老妇伦熟女老妇高清| 久久 成人 亚洲| 国产精品人妻久久久影院| 少妇丰满av| 亚洲天堂av无毛| 午夜福利,免费看| 中国美白少妇内射xxxbb| 久久毛片免费看一区二区三区| 中文字幕免费在线视频6| 另类精品久久| 亚洲精品国产av蜜桃| 欧美成人午夜免费资源| 亚洲丝袜综合中文字幕| 少妇高潮的动态图| 在线天堂最新版资源| 91aial.com中文字幕在线观看| 免费av中文字幕在线| 有码 亚洲区| 夜夜骑夜夜射夜夜干| 国产在线一区二区三区精| 在线精品无人区一区二区三| 国产精品偷伦视频观看了| 一级a做视频免费观看| 热99国产精品久久久久久7| 老熟女久久久| 久久精品人人爽人人爽视色| 少妇丰满av| 国产精品欧美亚洲77777| 中文字幕精品免费在线观看视频 | 2022亚洲国产成人精品| 在线观看国产h片| 国产精品嫩草影院av在线观看| 内地一区二区视频在线| 精品酒店卫生间| 男人操女人黄网站| 精品午夜福利在线看| 国模一区二区三区四区视频| 涩涩av久久男人的天堂| 国产高清不卡午夜福利| 精品一区二区免费观看| 日韩成人av中文字幕在线观看| 国产乱来视频区| 久热这里只有精品99| 麻豆成人av视频| 免费看不卡的av| 国产精品一区二区三区四区免费观看| 亚洲经典国产精华液单| 老司机影院毛片| 国产白丝娇喘喷水9色精品| 大陆偷拍与自拍| 国产免费福利视频在线观看| a级片在线免费高清观看视频| 午夜影院在线不卡| 老熟女久久久| 自线自在国产av| 亚洲av成人精品一区久久| h视频一区二区三区| 国产探花极品一区二区| 中文字幕人妻丝袜制服| 最新的欧美精品一区二区| 欧美bdsm另类| 一本—道久久a久久精品蜜桃钙片| 久久久国产一区二区| 激情五月婷婷亚洲| 午夜免费鲁丝| 午夜福利视频精品| 麻豆乱淫一区二区| 亚洲人成网站在线播| 国产亚洲av片在线观看秒播厂| 人妻一区二区av| 日本免费在线观看一区| 这个男人来自地球电影免费观看 | 97精品久久久久久久久久精品| 天堂俺去俺来也www色官网| 欧美丝袜亚洲另类| 国产淫语在线视频| 国产成人一区二区在线| 国产在线视频一区二区| 狠狠婷婷综合久久久久久88av| 制服丝袜香蕉在线| 只有这里有精品99| 91午夜精品亚洲一区二区三区| 精品久久久久久电影网| 少妇人妻久久综合中文| 伦理电影免费视频| 王馨瑶露胸无遮挡在线观看| 51国产日韩欧美| 亚洲精品一区蜜桃| 久久这里有精品视频免费| 99视频精品全部免费 在线| 成人手机av| 精品午夜福利在线看| 成人影院久久| 国产成人免费无遮挡视频| av女优亚洲男人天堂| 飞空精品影院首页| 亚洲精品一区蜜桃| 午夜免费男女啪啪视频观看| 80岁老熟妇乱子伦牲交| 99国产综合亚洲精品| 久久国产精品大桥未久av| 亚洲欧美一区二区三区黑人 | 欧美xxⅹ黑人| 美女脱内裤让男人舔精品视频| 大陆偷拍与自拍| 亚洲欧洲国产日韩| 国产在视频线精品| 极品少妇高潮喷水抽搐| 成人亚洲欧美一区二区av| 极品人妻少妇av视频| 国产亚洲精品久久久com| 亚洲三级黄色毛片| 天天躁夜夜躁狠狠久久av| 日韩强制内射视频| 免费看av在线观看网站| 日韩熟女老妇一区二区性免费视频| 国产国语露脸激情在线看| 国模一区二区三区四区视频| 人妻少妇偷人精品九色| 老熟女久久久| 国产成人freesex在线| 男女免费视频国产| 欧美bdsm另类| 欧美另类一区| 国内精品宾馆在线| 亚洲国产成人一精品久久久| 久久亚洲国产成人精品v| 免费观看av网站的网址| 人成视频在线观看免费观看| 天堂中文最新版在线下载| 免费高清在线观看视频在线观看| 伊人亚洲综合成人网| 日本av免费视频播放| 亚洲精品视频女| 日韩中文字幕视频在线看片| 亚洲精品av麻豆狂野| 欧美精品亚洲一区二区| 国产av精品麻豆| 亚洲欧美精品自产自拍| 亚洲欧洲国产日韩| 日韩成人av中文字幕在线观看| 最近2019中文字幕mv第一页| 涩涩av久久男人的天堂| 国产成人freesex在线| 日本与韩国留学比较| 亚洲欧洲日产国产| av卡一久久| 超碰97精品在线观看| av视频免费观看在线观看| 久久人妻熟女aⅴ| 又大又黄又爽视频免费| 国产综合精华液| 高清不卡的av网站| 国产精品国产三级国产专区5o| 久久久久久久久久久丰满| 精品少妇黑人巨大在线播放| 免费观看av网站的网址| 欧美亚洲 丝袜 人妻 在线| 搡女人真爽免费视频火全软件| 色哟哟·www| 亚洲婷婷狠狠爱综合网| 在线观看www视频免费| 99久久精品一区二区三区| 一级毛片黄色毛片免费观看视频| 久久久久精品久久久久真实原创| 激情五月婷婷亚洲| av女优亚洲男人天堂| 国产爽快片一区二区三区| 在线播放无遮挡| 国产精品久久久久久久电影| 哪个播放器可以免费观看大片| 国国产精品蜜臀av免费| 老司机亚洲免费影院| 晚上一个人看的免费电影| 亚洲av在线观看美女高潮| 免费看av在线观看网站| 中文字幕久久专区| 免费大片黄手机在线观看| 亚洲精品久久久久久婷婷小说| 亚洲精品日韩在线中文字幕| 一区二区三区精品91| 国产亚洲最大av| 午夜激情福利司机影院| 日韩视频在线欧美| .国产精品久久| 国产极品粉嫩免费观看在线 | 日本午夜av视频| 成人18禁高潮啪啪吃奶动态图 | 狂野欧美白嫩少妇大欣赏| 欧美一级a爱片免费观看看| 中文字幕人妻丝袜制服| 国产一区二区在线观看av| 午夜福利视频精品| 曰老女人黄片| 欧美97在线视频| 亚洲第一区二区三区不卡| 日韩大片免费观看网站| 国产精品熟女久久久久浪| 亚洲图色成人| 九草在线视频观看| 国产女主播在线喷水免费视频网站| 国产精品国产三级专区第一集| av一本久久久久| 一区二区三区乱码不卡18| 天堂俺去俺来也www色官网| 亚洲精品一区蜜桃| videosex国产| 哪个播放器可以免费观看大片| 免费观看无遮挡的男女| 青春草国产在线视频| 熟女av电影| www.av在线官网国产| 精品久久久久久久久亚洲| 天堂8中文在线网| 国产精品久久久久久久久免| 熟女av电影| 欧美bdsm另类| 国产高清有码在线观看视频| 午夜精品国产一区二区电影| 在线精品无人区一区二区三| 免费高清在线观看日韩| 超色免费av| 国产日韩一区二区三区精品不卡 | 日韩 亚洲 欧美在线| 麻豆精品久久久久久蜜桃| 久久精品久久久久久久性| 高清视频免费观看一区二区| 欧美精品一区二区大全| 日韩三级伦理在线观看| 久久精品国产亚洲网站| 女人久久www免费人成看片| 多毛熟女@视频| 亚洲精品乱码久久久v下载方式| 亚洲av日韩在线播放| 在线观看美女被高潮喷水网站| 中文字幕久久专区| 国产精品欧美亚洲77777| 一区二区日韩欧美中文字幕 | 亚洲av福利一区| 午夜影院在线不卡| 有码 亚洲区| 日韩一区二区视频免费看| 2018国产大陆天天弄谢| 热99久久久久精品小说推荐| 国产成人av激情在线播放 | 精品久久国产蜜桃| 久久毛片免费看一区二区三区| 伦理电影免费视频| 国产精品99久久久久久久久| 久久 成人 亚洲| 国产成人精品在线电影| 亚洲婷婷狠狠爱综合网| 91精品三级在线观看| 一边亲一边摸免费视频| 自线自在国产av| 国产精品 国内视频| 国产精品久久久久久精品电影小说| 国产伦理片在线播放av一区| 免费不卡的大黄色大毛片视频在线观看| 最近手机中文字幕大全| 十八禁网站网址无遮挡| 国产男人的电影天堂91| 中文字幕亚洲精品专区| 国产成人精品无人区| 国产免费视频播放在线视频| 国产精品不卡视频一区二区| 日韩大片免费观看网站| 国产亚洲一区二区精品| 亚洲第一区二区三区不卡| 久久99热这里只频精品6学生| 男女无遮挡免费网站观看| 3wmmmm亚洲av在线观看| 日韩免费高清中文字幕av| 一级毛片aaaaaa免费看小| 亚洲国产精品专区欧美| 久久精品久久久久久久性| 欧美三级亚洲精品| 制服人妻中文乱码| 国产在线一区二区三区精| 亚洲av.av天堂| 日本欧美视频一区| 色吧在线观看| 十八禁高潮呻吟视频| 少妇人妻 视频| 嘟嘟电影网在线观看| 日韩一本色道免费dvd| 亚洲av免费高清在线观看| 人人妻人人澡人人爽人人夜夜| 制服人妻中文乱码| 一级毛片aaaaaa免费看小| 在线观看三级黄色| 丁香六月天网| 成人国产麻豆网| 欧美亚洲日本最大视频资源| 免费高清在线观看视频在线观看| 一级黄片播放器| 成人毛片a级毛片在线播放| 婷婷色综合www| 久久久久久久久久人人人人人人| av在线播放精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产老妇伦熟女老妇高清| 老司机影院毛片| 好男人视频免费观看在线| 久久女婷五月综合色啪小说| 校园人妻丝袜中文字幕| 国产伦理片在线播放av一区| 亚洲内射少妇av| 精品久久久噜噜| 日韩大片免费观看网站| 九九爱精品视频在线观看| 国产精品成人在线| 亚洲欧美日韩另类电影网站| 人妻夜夜爽99麻豆av| 最近手机中文字幕大全| 看免费成人av毛片| 色视频在线一区二区三区| 国内精品宾馆在线| 岛国毛片在线播放| 久久99热这里只频精品6学生| 人妻一区二区av| 激情五月婷婷亚洲| 精品一区二区三卡| 女人精品久久久久毛片| 亚洲不卡免费看| 国产熟女午夜一区二区三区 | 人妻夜夜爽99麻豆av| av电影中文网址| 国产白丝娇喘喷水9色精品| 天堂8中文在线网| 日韩一区二区视频免费看| 十分钟在线观看高清视频www| 亚洲成人一二三区av| 亚洲av中文av极速乱| 黄色视频在线播放观看不卡| 天堂中文最新版在线下载| 好男人视频免费观看在线| 看十八女毛片水多多多| 国产高清有码在线观看视频| 99re6热这里在线精品视频| 国产成人精品婷婷| 一个人免费看片子| 亚洲怡红院男人天堂| 亚洲人成网站在线播| 精品国产乱码久久久久久小说| 999精品在线视频| 这个男人来自地球电影免费观看 | 在线观看人妻少妇| 国产伦精品一区二区三区视频9| 国产精品久久久久久久久免| 在线观看免费日韩欧美大片 | 国产精品国产三级专区第一集| 亚洲欧美一区二区三区国产| 天天躁夜夜躁狠狠久久av| 国产精品99久久99久久久不卡 | 99久久中文字幕三级久久日本| 中文字幕人妻熟人妻熟丝袜美| 九九久久精品国产亚洲av麻豆| 欧美成人精品欧美一级黄| 3wmmmm亚洲av在线观看| 欧美成人精品欧美一级黄| 亚洲天堂av无毛| 亚洲精品456在线播放app| 免费观看a级毛片全部| 国产一级毛片在线| 亚洲熟女精品中文字幕| 欧美3d第一页| 中文欧美无线码| 成人二区视频| 国产 一区精品| 国产精品一二三区在线看| 国产精品成人在线| 亚洲成人一二三区av| 制服诱惑二区| 性色av一级| 如日韩欧美国产精品一区二区三区 | 大片电影免费在线观看免费| 成人综合一区亚洲| 国产精品人妻久久久影院| 日本欧美视频一区| 高清av免费在线| 午夜免费观看性视频| 国产视频首页在线观看| kizo精华| 男女高潮啪啪啪动态图| 日韩在线高清观看一区二区三区| av天堂久久9| 亚洲国产精品999| 国产日韩一区二区三区精品不卡 | 欧美日韩成人在线一区二区| 日韩 亚洲 欧美在线| 亚洲精品自拍成人| 99热6这里只有精品| 亚洲精品,欧美精品| 亚洲av综合色区一区| 亚洲天堂av无毛| 免费黄网站久久成人精品| 搡女人真爽免费视频火全软件| 日日啪夜夜爽| 一级毛片 在线播放| 丝袜喷水一区| 亚洲精品色激情综合| 亚洲三级黄色毛片| 亚洲欧洲日产国产| 国产午夜精品一二区理论片| 欧美性感艳星| 精品人妻熟女毛片av久久网站| 亚洲,一卡二卡三卡| 伦理电影大哥的女人| 亚洲,一卡二卡三卡| videosex国产| 97在线人人人人妻| 最黄视频免费看| 成人二区视频| 3wmmmm亚洲av在线观看| 中文字幕免费在线视频6| 日韩不卡一区二区三区视频在线| 国产日韩欧美亚洲二区| 国产成人精品婷婷| 国产精品国产av在线观看| 久久影院123| 久久久久国产网址| av国产久精品久网站免费入址| 精品熟女少妇av免费看| 下体分泌物呈黄色| av播播在线观看一区| 香蕉精品网在线| 日本vs欧美在线观看视频| 久久久久久久亚洲中文字幕| 有码 亚洲区| 伦精品一区二区三区| 啦啦啦视频在线资源免费观看| 国产高清不卡午夜福利| 老司机亚洲免费影院| 亚洲av中文av极速乱| 美女视频免费永久观看网站| 国产一区二区在线观看av| 亚洲av成人精品一区久久| 国产精品三级大全| 欧美+日韩+精品| 亚洲高清免费不卡视频| 寂寞人妻少妇视频99o| 久久人人爽人人片av| 美女主播在线视频| 丁香六月天网| 有码 亚洲区| 欧美+日韩+精品| 美女福利国产在线| 中文字幕人妻丝袜制服| 黄色视频在线播放观看不卡| 成人手机av| 丝瓜视频免费看黄片| 五月玫瑰六月丁香| 热99久久久久精品小说推荐| xxxhd国产人妻xxx| 大陆偷拍与自拍| 曰老女人黄片| 国产一区二区三区av在线| 纯流量卡能插随身wifi吗| 桃花免费在线播放| 亚洲欧美色中文字幕在线| 欧美另类一区| 一本—道久久a久久精品蜜桃钙片| 极品少妇高潮喷水抽搐| 午夜91福利影院| kizo精华| 中文字幕最新亚洲高清| 又黄又爽又刺激的免费视频.| 免费播放大片免费观看视频在线观看| 建设人人有责人人尽责人人享有的| 天天操日日干夜夜撸| 校园人妻丝袜中文字幕| 欧美日本中文国产一区发布| 久久人人爽av亚洲精品天堂| 欧美激情国产日韩精品一区|