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

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

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

    ?!〖t,寶金山,吳 洋,朱 波,張文剛,夏江威,宋禹昕,郭 鵬,徐凌洋,陳 燕,高 雪,張路培,高會江,李俊雅*

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

    ?

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

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

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

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

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

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

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

    1 材料與方法

    1.1材料

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

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

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

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

    1.3影響因素顯著性檢驗

    影響因素的處理和顯著性檢驗是通過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ìn)場重和育肥天數(shù)協(xié)變量,而出欄重、胴體重、屠宰率和凈肉率影響因素的顯著性檢驗分析模型中則包括進(jìn)場重和育肥天數(shù)協(xié)變量??傊?,對于表型性狀的校正,要把顯著的影響因素加入到校正模型中,而不顯著的影響因素要去掉。

    1.4遺傳參數(shù)估計

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

    y=1nu+Zg+e

    (2)

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

    (3)

    (4)

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

    表1出欄重、胴體重、屠宰率和凈肉率性狀影響因素不同水平的觀察個數(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觀察值個數(shù)N水平Level觀察值個數(shù)N水平Level觀察值個數(shù)N水平Level觀察值個數(shù)N水平Level觀察值個數(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)計量

    表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)條件下有較大的波動。而屠宰率和凈肉率的變異系數(shù)均在10%以下,表現(xiàn)出相對的穩(wěn)定性。圖1為各個性狀表型值每年均值的變化圖。

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

    如表3中所示,通過顯著性檢驗可以得知,母親年齡效應(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)的時候,把顯著的因素加入到模型進(jìn)行校正,而不顯著的因素去掉。

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

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

    Table 2Summary of important economic traits

    性狀Trait記錄個數(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)中的兩個影響因素(場效應(yīng)和年效應(yīng))對初生重、斷奶重、出欄重、胴體重、屠宰率和凈肉率的影響,見表4,可以看出,除初生重和出欄重不隨場的變化而變化外(P>0.05),斷奶重、胴體重、屠宰率和凈肉率4個性狀均呈現(xiàn)隨場效應(yīng)顯著變化的趨勢(P<0.05),而6個性狀均呈現(xiàn)隨年效應(yīng)顯著變化的趨勢(P<0.05)。

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

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

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

    表3不同性狀影響因素的顯著性檢驗結(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)濟性狀的遺傳參數(shù)估計

    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)濟性狀兩兩間的遺傳相關(guān)估計

    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]。由于西門塔爾牛具有生長速度快、耐粗飼、繁殖力強等特性,這使得中國西門塔爾牛肉用群體的建立顯得十分必要。人們可以立足已建立的中國西門塔爾?;蚪M選擇肉用參考群體的基礎(chǔ)上,完善肉牛的育種技術(shù)和建立全基因組選擇技術(shù)平臺,開展肉牛全基因組選擇、新基因發(fā)掘和功能解析等方面的研究[10-11],提供使肉牛種業(yè)可持續(xù)發(fā)展的技術(shù)支撐平臺,并培育我國自主的肉用西門塔爾牛新品系,擺脫我國肉牛種業(yè)依賴進(jìn)口的局面,為全國肉牛進(jìn)一步遺傳改良提供借鑒和參考。

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

    本研究估計的西門塔爾牛的初生重和斷奶重遺傳相關(guān)為0.57,小于M.B.Long等[23]報道的海福特牛初生重和斷奶重的遺傳相關(guān)值0.64。同時本研究估計的出欄重和胴體重的遺傳相關(guān)為0.94,與D.M.Marshall[24]報道的海福特牛出欄重和胴體重遺傳相關(guān)0.93相一致。而屠宰率和凈肉率的遺傳相關(guān),國外未見相關(guān)報道。本研究中所估計的遺傳參數(shù)與國外相比,其差異可能與數(shù)據(jù)量和估計方法有關(guān)。在育種工作中,若兩個性狀有較高的遺傳相關(guān),選擇其中一個性狀的同時會間接地改進(jìn)和提高另一性狀的選擇效果。因此,遺傳參數(shù)的估計在肉牛生產(chǎn)中具有較高實用價值,對進(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)濟性狀的遺傳參數(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].安徽農(nóng)業(yè)大學(xué)學(xué)報,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]范大有,許尚忠,李俊雅,等.中國西門塔爾牛次級性狀與生產(chǎn)性狀的遺傳統(tǒng)計分析[J].畜牧獸醫(yī)學(xué)報,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]霍靈光,田露,張越杰.中國牛肉需求量中長期預(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

    “十二五”國家科技支撐計劃(2011BAD28B04);國家自然科學(xué)基金(31372294);國家高技術(shù)研究發(fā)展計劃(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

    午夜福利在线观看吧| 我要搜黄色片| 色噜噜av男人的天堂激情| 国产一区二区三区视频了| 国产精品1区2区在线观看.| 少妇的逼水好多| 夜夜夜夜夜久久久久| 禁无遮挡网站| 精品国产乱码久久久久久男人| 美女高潮喷水抽搐中文字幕| 亚洲真实伦在线观看| 小说图片视频综合网站| 91在线观看av| 黑人欧美特级aaaaaa片| 真人做人爱边吃奶动态| 国产成人av激情在线播放| 久久人妻av系列| 国产精品永久免费网站| 99热这里只有精品一区 | 成人亚洲精品av一区二区| 国产日本99.免费观看| 最近在线观看免费完整版| 久久性视频一级片| 亚洲精品中文字幕一二三四区| 叶爱在线成人免费视频播放| 国内毛片毛片毛片毛片毛片| 久久精品国产亚洲av香蕉五月| 精品国产亚洲在线| 日韩中文字幕欧美一区二区| 欧美不卡视频在线免费观看| 日本黄色视频三级网站网址| 日本黄大片高清| 日韩人妻高清精品专区| 亚洲自拍偷在线| av国产免费在线观看| 亚洲成人中文字幕在线播放| 欧美色视频一区免费| www.www免费av| 给我免费播放毛片高清在线观看| 亚洲成a人片在线一区二区| 可以在线观看毛片的网站| 亚洲精品在线美女| 黄片小视频在线播放| 久久天躁狠狠躁夜夜2o2o| 亚洲avbb在线观看| 在线a可以看的网站| 国产精品香港三级国产av潘金莲| 丝袜人妻中文字幕| 怎么达到女性高潮| 好男人电影高清在线观看| 成年女人永久免费观看视频| 精品无人区乱码1区二区| 麻豆一二三区av精品| 午夜福利高清视频| 极品教师在线免费播放| 人妻久久中文字幕网| 日本熟妇午夜| 在线永久观看黄色视频| 一边摸一边抽搐一进一小说| 精品久久久久久久久久免费视频| 色哟哟哟哟哟哟| 国产精品一区二区三区四区久久| 欧美乱妇无乱码| 黑人欧美特级aaaaaa片| 又粗又爽又猛毛片免费看| 久久精品综合一区二区三区| 美女高潮的动态| 一本综合久久免费| 久久精品国产亚洲av香蕉五月| a级毛片在线看网站| 黄色丝袜av网址大全| 亚洲人成网站在线播放欧美日韩| 床上黄色一级片| 国产麻豆成人av免费视频| 中文字幕高清在线视频| 日韩欧美国产一区二区入口| 欧美激情在线99| 99热这里只有精品一区 | 噜噜噜噜噜久久久久久91| 国产成人av教育| 亚洲中文日韩欧美视频| 国产一区二区在线观看日韩 | 日韩高清综合在线| 国产精品野战在线观看| 久久国产乱子伦精品免费另类| 欧美日韩乱码在线| 国产aⅴ精品一区二区三区波| 日韩中文字幕欧美一区二区| 综合色av麻豆| 成人亚洲精品av一区二区| 欧美日韩精品网址| 精品国产美女av久久久久小说| 91在线精品国自产拍蜜月 | 国产三级黄色录像| a级毛片a级免费在线| 亚洲成av人片免费观看| 色av中文字幕| 亚洲精品乱码久久久v下载方式 | 亚洲成av人片免费观看| 很黄的视频免费| 久久久精品大字幕| 精品人妻1区二区| 国内毛片毛片毛片毛片毛片| 在线播放国产精品三级| 日本免费一区二区三区高清不卡| 国产伦人伦偷精品视频| 一卡2卡三卡四卡精品乱码亚洲| 国产精品98久久久久久宅男小说| 床上黄色一级片| 老司机午夜十八禁免费视频| 久久国产精品影院| 久久久久久久午夜电影| 成人三级做爰电影| 999久久久精品免费观看国产| 九色成人免费人妻av| 午夜激情福利司机影院| 国产精品1区2区在线观看.| 日韩欧美三级三区| xxxwww97欧美| 在线十欧美十亚洲十日本专区| 搡老岳熟女国产| 啪啪无遮挡十八禁网站| 观看免费一级毛片| 日本与韩国留学比较| 久久久久久久精品吃奶| 怎么达到女性高潮| av国产免费在线观看| 国产精品久久视频播放| 在线国产一区二区在线| 婷婷六月久久综合丁香| 在线观看免费午夜福利视频| 色老头精品视频在线观看| 国产免费av片在线观看野外av| 男插女下体视频免费在线播放| 精品国内亚洲2022精品成人| 国产不卡一卡二| 99riav亚洲国产免费| 99精品在免费线老司机午夜| 亚洲国产看品久久| 熟女人妻精品中文字幕| 国产乱人伦免费视频| 岛国视频午夜一区免费看| 亚洲人成电影免费在线| 亚洲色图av天堂| 国产精品免费一区二区三区在线| 欧美色欧美亚洲另类二区| 欧美精品啪啪一区二区三区| 久久99热这里只有精品18| 最近在线观看免费完整版| 亚洲成av人片在线播放无| 一a级毛片在线观看| 午夜福利视频1000在线观看| 一本综合久久免费| 十八禁网站免费在线| 欧美日韩福利视频一区二区| 国语自产精品视频在线第100页| 一级a爱片免费观看的视频| 长腿黑丝高跟| 女人被狂操c到高潮| 国产探花在线观看一区二区| 婷婷精品国产亚洲av在线| 欧美激情在线99| 亚洲人成网站高清观看| 18禁黄网站禁片免费观看直播| 亚洲黑人精品在线| 小蜜桃在线观看免费完整版高清| 久久国产乱子伦精品免费另类| 国产精品久久视频播放| 搡老熟女国产l中国老女人| 久久婷婷人人爽人人干人人爱| 99在线视频只有这里精品首页| 国产伦在线观看视频一区| 久久精品aⅴ一区二区三区四区| 成人特级av手机在线观看| 欧美成人免费av一区二区三区| 夜夜看夜夜爽夜夜摸| 一个人免费在线观看的高清视频| 最新美女视频免费是黄的| 欧美一区二区精品小视频在线| 亚洲av电影在线进入| 在线观看午夜福利视频| 18禁国产床啪视频网站| 久久天堂一区二区三区四区| 嫩草影院精品99| 亚洲五月婷婷丁香| 欧洲精品卡2卡3卡4卡5卡区| 中文资源天堂在线| 日韩欧美国产一区二区入口| 欧美黑人欧美精品刺激| 2021天堂中文幕一二区在线观| 精品乱码久久久久久99久播| 岛国在线免费视频观看| 免费在线观看亚洲国产| 黑人欧美特级aaaaaa片| 在线观看66精品国产| 人人妻,人人澡人人爽秒播| 天堂√8在线中文| 成人无遮挡网站| 校园春色视频在线观看| 亚洲无线观看免费| 午夜视频精品福利| 亚洲,欧美精品.| 国产精品 国内视频| 中文字幕熟女人妻在线| 最新美女视频免费是黄的| 女人高潮潮喷娇喘18禁视频| 久久久色成人| 国产黄片美女视频| 99国产综合亚洲精品| 精品久久久久久久久久久久久| 97碰自拍视频| 中亚洲国语对白在线视频| 亚洲av中文字字幕乱码综合| 色综合亚洲欧美另类图片| 18禁黄网站禁片午夜丰满| h日本视频在线播放| 九九热线精品视视频播放| 色尼玛亚洲综合影院| 国产爱豆传媒在线观看| 十八禁人妻一区二区| 在线观看舔阴道视频| 色吧在线观看| 美女扒开内裤让男人捅视频| 亚洲自拍偷在线| 欧美成狂野欧美在线观看| 天天一区二区日本电影三级| 国产精品爽爽va在线观看网站| 丰满人妻熟妇乱又伦精品不卡| 日韩精品青青久久久久久| 中国美女看黄片| 99久久99久久久精品蜜桃| 精品国产美女av久久久久小说| 国产日本99.免费观看| 全区人妻精品视频| 免费搜索国产男女视频| 亚洲一区二区三区不卡视频| 19禁男女啪啪无遮挡网站| 午夜成年电影在线免费观看| av中文乱码字幕在线| 午夜免费激情av| 国产伦在线观看视频一区| 亚洲精品色激情综合| 欧美成人免费av一区二区三区| 欧美一区二区精品小视频在线| 嫩草影院入口| 日本熟妇午夜| 叶爱在线成人免费视频播放| 可以在线观看毛片的网站| 99久久精品热视频| 色综合欧美亚洲国产小说| 精品国产超薄肉色丝袜足j| 午夜亚洲福利在线播放| 人人妻人人看人人澡| 久久久久久久久久黄片| svipshipincom国产片| 国产三级中文精品| 欧美日韩福利视频一区二区| 国产一级毛片七仙女欲春2| 国产精品av久久久久免费| 日本一本二区三区精品| 草草在线视频免费看| 精品熟女少妇八av免费久了| 不卡一级毛片| 高清毛片免费观看视频网站| 婷婷精品国产亚洲av| 国产精品精品国产色婷婷| 欧美一级a爱片免费观看看| 亚洲国产看品久久| 狂野欧美激情性xxxx| ponron亚洲| 在线观看午夜福利视频| 免费无遮挡裸体视频| 露出奶头的视频| 精品一区二区三区av网在线观看| 国产私拍福利视频在线观看| 久久伊人香网站| 人人妻人人澡欧美一区二区| x7x7x7水蜜桃| 亚洲 国产 在线| 亚洲五月婷婷丁香| 国产精品亚洲一级av第二区| 日本 欧美在线| 亚洲专区字幕在线| 久久久国产精品麻豆| 国产精品一区二区免费欧美| 色综合亚洲欧美另类图片| 亚洲精品456在线播放app | 国产亚洲av嫩草精品影院| 亚洲色图av天堂| 精品国产乱码久久久久久男人| 免费观看精品视频网站| 波多野结衣高清作品| 亚洲国产精品久久男人天堂| 国产主播在线观看一区二区| 91麻豆av在线| 天天躁狠狠躁夜夜躁狠狠躁| 99久久精品国产亚洲精品| 欧美一区二区精品小视频在线| 嫁个100分男人电影在线观看| 日韩精品中文字幕看吧| 18禁美女被吸乳视频| 亚洲国产欧美网| 国产成人福利小说| 99热精品在线国产| 99在线人妻在线中文字幕| 国产又黄又爽又无遮挡在线| 在线永久观看黄色视频| 99久久成人亚洲精品观看| 色噜噜av男人的天堂激情| 日本黄色片子视频| 欧美日本亚洲视频在线播放| 老鸭窝网址在线观看| 国产免费男女视频| 性色av乱码一区二区三区2| 色综合婷婷激情| 深夜精品福利| 欧美丝袜亚洲另类 | 夜夜爽天天搞| 欧美另类亚洲清纯唯美| 亚洲天堂国产精品一区在线| 精品电影一区二区在线| 欧美性猛交黑人性爽| 色综合欧美亚洲国产小说| 免费在线观看视频国产中文字幕亚洲| 日韩 欧美 亚洲 中文字幕| 久久精品91无色码中文字幕| 日韩大尺度精品在线看网址| ponron亚洲| 岛国视频午夜一区免费看| 亚洲欧美激情综合另类| 亚洲专区字幕在线| 九九久久精品国产亚洲av麻豆 | 淫妇啪啪啪对白视频| 1000部很黄的大片| a级毛片a级免费在线| 中出人妻视频一区二区| 又紧又爽又黄一区二区| 久久久久久国产a免费观看| 欧美丝袜亚洲另类 | 精品久久蜜臀av无| 99re在线观看精品视频| 国产成人欧美在线观看| 欧美在线黄色| 99热这里只有精品一区 | 桃色一区二区三区在线观看| 色在线成人网| 日本撒尿小便嘘嘘汇集6| 嫩草影视91久久| 99精品久久久久人妻精品| 三级国产精品欧美在线观看 | 亚洲真实伦在线观看| 99热6这里只有精品| 国产爱豆传媒在线观看| 精品一区二区三区四区五区乱码| 人人妻人人澡欧美一区二区| 免费看日本二区| 亚洲欧美日韩无卡精品| 老司机午夜福利在线观看视频| 亚洲欧美日韩无卡精品| 欧美极品一区二区三区四区| 日本 av在线| 校园春色视频在线观看| 国产精品久久久人人做人人爽| 色尼玛亚洲综合影院| 久久婷婷人人爽人人干人人爱| 法律面前人人平等表现在哪些方面| 巨乳人妻的诱惑在线观看| 精品久久蜜臀av无| 久久精品夜夜夜夜夜久久蜜豆| 熟女少妇亚洲综合色aaa.| 啦啦啦免费观看视频1| 欧美中文日本在线观看视频| 国产精品香港三级国产av潘金莲| 成人18禁在线播放| 国产精品av视频在线免费观看| 午夜成年电影在线免费观看| 免费电影在线观看免费观看| 久久久精品大字幕| 久久久久久国产a免费观看| 国产1区2区3区精品| 国产精品日韩av在线免费观看| 一区二区三区国产精品乱码| 欧美性猛交╳xxx乱大交人| 香蕉久久夜色| 99国产精品一区二区蜜桃av| 国产精品99久久99久久久不卡| a级毛片在线看网站| 麻豆成人av在线观看| 又大又爽又粗| 国产一级毛片七仙女欲春2| 欧美色视频一区免费| 国内少妇人妻偷人精品xxx网站 | 日本 av在线| 国产精品永久免费网站| 男女床上黄色一级片免费看| 精品久久久久久久人妻蜜臀av| 亚洲av电影在线进入| 国产亚洲精品一区二区www| 国产黄色小视频在线观看| 国产精品乱码一区二三区的特点| 母亲3免费完整高清在线观看| 又紧又爽又黄一区二区| 欧美色视频一区免费| 18禁观看日本| 女人被狂操c到高潮| 五月伊人婷婷丁香| 级片在线观看| 九色成人免费人妻av| 国产成人精品久久二区二区免费| 男女之事视频高清在线观看| 亚洲专区国产一区二区| 亚洲欧美日韩高清专用| 午夜福利在线在线| 成熟少妇高潮喷水视频| 丰满人妻一区二区三区视频av | 久久精品亚洲精品国产色婷小说| 亚洲aⅴ乱码一区二区在线播放| 黄色日韩在线| 真人一进一出gif抽搐免费| 亚洲专区中文字幕在线| 又爽又黄无遮挡网站| 午夜a级毛片| 国产又黄又爽又无遮挡在线| 女生性感内裤真人,穿戴方法视频| 亚洲人成网站高清观看| 日本 av在线| 国产又黄又爽又无遮挡在线| 国产三级黄色录像| 老司机深夜福利视频在线观看| 国产高清有码在线观看视频| 99视频精品全部免费 在线 | 国产又色又爽无遮挡免费看| 日韩欧美国产在线观看| 床上黄色一级片| 91在线观看av| 夜夜爽天天搞| 我的老师免费观看完整版| 在线免费观看不下载黄p国产 | 精品一区二区三区av网在线观看| 90打野战视频偷拍视频| 在线十欧美十亚洲十日本专区| 天天躁日日操中文字幕| 免费在线观看影片大全网站| 国产亚洲欧美98| 18禁黄网站禁片免费观看直播| 亚洲第一电影网av| 国产午夜精品久久久久久| 韩国av一区二区三区四区| 亚洲欧美日韩卡通动漫| 动漫黄色视频在线观看| 精品一区二区三区四区五区乱码| 男人舔女人下体高潮全视频| 在线播放国产精品三级| 日韩三级视频一区二区三区| 欧美绝顶高潮抽搐喷水| 18禁裸乳无遮挡免费网站照片| 99久久久亚洲精品蜜臀av| 三级毛片av免费| 欧美av亚洲av综合av国产av| 天堂影院成人在线观看| 欧美一区二区国产精品久久精品| 757午夜福利合集在线观看| 性色avwww在线观看| 午夜福利高清视频| 国产成人欧美在线观看| 国产精品久久久久久精品电影| 99热精品在线国产| 国内毛片毛片毛片毛片毛片| 国产精品99久久99久久久不卡| 亚洲色图 男人天堂 中文字幕| 人人妻,人人澡人人爽秒播| 人妻久久中文字幕网| 国产精品免费一区二区三区在线| 真人做人爱边吃奶动态| 51午夜福利影视在线观看| 九九在线视频观看精品| 欧美一区二区国产精品久久精品| 国产亚洲精品av在线| 91久久精品国产一区二区成人 | 午夜福利在线在线| 熟女少妇亚洲综合色aaa.| 五月玫瑰六月丁香| 黑人操中国人逼视频| 啪啪无遮挡十八禁网站| 色哟哟哟哟哟哟| 成年免费大片在线观看| 久久亚洲精品不卡| 国产探花在线观看一区二区| 免费观看人在逋| 国产亚洲av高清不卡| 一进一出抽搐动态| 亚洲国产看品久久| 亚洲在线观看片| 免费看美女性在线毛片视频| 久久国产精品影院| 黄色片一级片一级黄色片| 国产真人三级小视频在线观看| 最近最新中文字幕大全电影3| 国产一区二区激情短视频| 在线十欧美十亚洲十日本专区| 国产精华一区二区三区| 国产三级中文精品| 欧美中文综合在线视频| 国产黄a三级三级三级人| 日本熟妇午夜| 欧美成狂野欧美在线观看| 欧美日本亚洲视频在线播放| 久久久久精品国产欧美久久久| 久久久久久久午夜电影| 亚洲午夜精品一区,二区,三区| 国产欧美日韩一区二区精品| 亚洲欧美日韩无卡精品| 最近最新中文字幕大全电影3| 久久久精品欧美日韩精品| 黄色成人免费大全| 18禁黄网站禁片午夜丰满| 国产精品香港三级国产av潘金莲| 一级毛片女人18水好多| 日本a在线网址| 不卡一级毛片| 免费看日本二区| 国产成年人精品一区二区| 国产亚洲精品一区二区www| 十八禁人妻一区二区| 亚洲真实伦在线观看| 国产人伦9x9x在线观看| 嫩草影院精品99| 小说图片视频综合网站| 免费在线观看亚洲国产| 久久天堂一区二区三区四区| 亚洲午夜理论影院| 欧美xxxx黑人xx丫x性爽| 久久欧美精品欧美久久欧美| 蜜桃久久精品国产亚洲av| 国产一区二区激情短视频| 亚洲专区字幕在线| 两性夫妻黄色片| 午夜福利18| 99热这里只有是精品50| 一夜夜www| 黄色 视频免费看| 精品午夜福利视频在线观看一区| 久久精品影院6| 国产精品国产高清国产av| 欧美日韩一级在线毛片| 久久精品夜夜夜夜夜久久蜜豆| 国产激情偷乱视频一区二区| 99国产综合亚洲精品| 欧美不卡视频在线免费观看| 性色avwww在线观看| 免费在线观看影片大全网站| 99在线视频只有这里精品首页| 亚洲精品在线观看二区| 老司机福利观看| 老汉色av国产亚洲站长工具| 欧美日韩瑟瑟在线播放| 久99久视频精品免费| 国产伦精品一区二区三区视频9 | 亚洲欧美精品综合一区二区三区| 国产成人一区二区三区免费视频网站| 丰满的人妻完整版| netflix在线观看网站| 性色avwww在线观看| 亚洲乱码一区二区免费版| 淫妇啪啪啪对白视频| 成人午夜高清在线视频| 国产爱豆传媒在线观看| 在线播放国产精品三级| 香蕉丝袜av| 人妻夜夜爽99麻豆av| 99热这里只有精品一区 | 国产探花在线观看一区二区| 色播亚洲综合网| 欧美又色又爽又黄视频| 久久热在线av| а√天堂www在线а√下载| 午夜福利高清视频| 高清在线国产一区| 我要搜黄色片| 人妻夜夜爽99麻豆av| 日韩欧美在线二视频| 亚洲性夜色夜夜综合| 亚洲中文字幕一区二区三区有码在线看 | netflix在线观看网站| 国产精华一区二区三区| 午夜福利欧美成人| 午夜日韩欧美国产| 久久性视频一级片| www日本黄色视频网| 亚洲av中文字字幕乱码综合| 十八禁人妻一区二区| 啦啦啦免费观看视频1| 亚洲中文日韩欧美视频| 国产亚洲av嫩草精品影院| 国产精品香港三级国产av潘金莲| 国产亚洲精品综合一区在线观看| 青草久久国产| 99re在线观看精品视频| 99久久99久久久精品蜜桃| 成年人黄色毛片网站| av欧美777| 国产亚洲精品一区二区www| 麻豆国产av国片精品| 色综合亚洲欧美另类图片| 免费在线观看视频国产中文字幕亚洲| 操出白浆在线播放| 亚洲,欧美精品.| 欧美不卡视频在线免费观看| av片东京热男人的天堂| 变态另类丝袜制服| 亚洲人成网站在线播放欧美日韩| 色综合欧美亚洲国产小说| 中文字幕熟女人妻在线| 黄色 视频免费看| 国产高清videossex|