Molecular maker-assisted selection of VRTN gene in two Landrace swine nucleus herds
-
摘要:
VRTN基因是影响猪胸椎数的主效基因，其优势等位基因（Q）比劣势等位基因（q）具有多1对肋骨的效应。本文旨在探索VRTN基因在2个长白种猪群中的分子标记辅助选择方法，以期提高配套肉猪的肋骨数。采用PCR方法检测影响胸椎数（肋骨对数）的VRTN基因在温氏集团420头商品肉猪以及2个专门化长白种猪（913头W51系长白，240头W52系长白）群体中的分布情况。通过商品肉猪的效应分析，验证了在终端商品肉猪中QQ型较qq型猪只多近1对肋骨，说明分子标记辅助选择进行基因改良是可行的。试验结果显示：VRTN基因在2个长白群体中的等位基因频率相差不大，有利等位基因频率分别为0.844和0.905。通过基因型与选育性状表型的相关分析，发现W51系长白猪中，优势基因型猪只体长优势明显，对其他选育性状无影响。W52系中，优势等位基因与达100 kg/115 kg日龄估计育种值呈显著负相关，与其他性状无显著相关性。VRTN优势等位基因型的效应与选育目标一致，2个品系均可开展分子标记辅助选择。结合现场实践育种，W51通过种群更新优化，经过半年就实现了优势等位基因频率100%的完全纯化。W52在扩群的同时，经过2年达到基本纯化。
Abstract:VRTN gene is the major gene affecting the number of thoracic vertebra in pigs, which QQ genotype (ins/ins) individuals have one double ribs more than qq (-/-) genotype ones. The molecular marker-assisted selection method of VRTN gene in two landrace pig populations was explored to improve the number of ribs of mating pigs. In the study, PCR were used to detect the distribution of VRTN genes in 420 commercial pigs of Wens Group and 2 specialized landrace pigs (913 W51 Landrace and 240 W52 Landrace), which affected thoracic vertebra number (rib number). It is verified that QQ type had nearly one pair of ribs more than qq type in Wens terminal pigs by effect analysis of commercial pork pigs, which indicated that gene improvement could be carried out through molecular marker-assisted selection. The allele frequencies of VRTN gene in the W51, W52 Landrace populations were not significantly different, and the favorable allele frequencies were 0.844 and 0.905 respectively. Through the correlation analysis between genotype and phenotypes of breeding traits, it was found that among W51 Landrace pigs, the dominant genotype had obvious body length advantage and had no influence on other breeding traits. In W52 line, the dominant allele had a significant negative correlation with the estimated breeding value of age at 100 kg/115 kg body weight, but had no significant correlation with other traits. The effect of VRTN dominant allele type was consistent with the breeding goal. Molecular marker-assisted selection could be carried out for W51 and W52 lines both. Combining with field practice breeding, W51 achieved 100% complete purification of dominant allele frequency after half a year through population renewal and optimization. W52 was basically purified after 2 years while expanding its population.
-
Keywords:
- VRTN gene /
- pig nucleus herd /
- molecular marker-assisted selection
-
目前，我国食糖生产成本过高、蔗糖产业国际竞争力十分低下、糖价大大高于国际标准[1]，巴西、美国等进口蔗糖的到岸价比我国蔗糖的完税价要低8到10个百分点[2]。造成这一情况的主要原因是我国大部分蔗区的生产收获基本靠人工作业，甘蔗收获机械化程度低。断尾除叶作为甘蔗收获的关键环节，是制约甘蔗收获机械发展的瓶颈[3]。甘蔗断尾是刚性体(机械)和柔性体(甘蔗尾茎)相互作用的过程[4]，现阶段研究[5-9]通常是根据经验估算从而进行断尾机构的设计研究，对甘蔗尾茎力学特性的分析很少，难以从本质上揭示甘蔗断尾的机理。因此，研究甘蔗尾茎的力学特性对于甘蔗断尾机械的设计开发、建立虚拟蔗尾模型进行断尾动力学仿真及有限元分析是非常必要的。在甘蔗尾茎的力学特性研究中，泊松比参数是不可或缺的，它是建立甘蔗茎秆力学模型和虚拟仿真模型的必要参数之一[10]。鲜见明确的测定甘蔗尾茎泊松比的方法[11-12]，构建甘蔗力学模型时所需的泊松比参数是根据横观各向异性材料力学公式及甘蔗尾茎各组分弹性模量推导得出[10]，或者部分泊松比参数根据农业生物力学由其他物料泊松比比拟得出[13]。由于误差的累加以及各种农业物料的微观结构的差异，以上方法间接得到的甘蔗泊松比参数会产生较大误差[14]，而用电测法直接试验得到泊松比参数，过程大为简化，可缩短科研周期[15]。
本研究借鉴竹子、苜蓿秆、苎麻、玉米等茎秆类作物的力学分析方法[16-21]，结合甘蔗茎秆力学特性与复合材料力学理论，建立甘蔗尾茎力学模型；借鉴对农林作物泊松比的相关研究[22-25]，将电测法试验与理论分析相结合，确立甘蔗尾茎的泊松比参数，以期为深入研究甘蔗断尾机理提供理论基础。
1. 甘蔗尾茎结构及力学模型
1.1 甘蔗尾茎结构特点
甘蔗属于实芯秸秆类植物，内部有明显的芯结构，且不同部位具有不同的结构特点。甘蔗尾茎的横截面近似圆形，横截面由内到外分为甘蔗芯和甘蔗皮2部分，蔗茎由若干节和节间组成。甘蔗皮特点是细胞小，排列紧密；甘蔗芯有许多纤维状结构，这些纤维沿轴向整齐排列；甘蔗节的纤维结构较甘蔗芯紧密。因此甘蔗茎秆在微观结构上不连续，并且在空间排列上具有明显的方向性。成熟期蔗茎尾梢被叶鞘包裹，以生长点以下约1~5节包裹最为紧密，这部分蔗尾含糖极少，脆弱易断，前期试验证明蔗尾茎生长点以下1~5节的抗弯强度明显低于6~9节的抗弯强度，差异极显著，蔗尾在第5节位置断裂，可以将甘蔗尾部5~6片包裹紧密的青叶连同茎秆尾梢部分一起断除，大大降低剥叶难度，有效提高原料蔗的质量[26-28]。所以建立甘蔗材料模型进行断尾动力学研究的时候，可将甘蔗整秆简化为蔗尾部分和蔗身部分，中间是蔗节连接，蔗尾部分和蔗身部分及蔗节部分力学性能差异很大，本研究将蔗尾生长点以下第5节的力学参数作为蔗尾部分力学参数，第6节的力学参数作为蔗身部分的力学参数。
1.2 甘蔗尾茎力学模型
甘蔗尾茎皮、芯、节3部分可以抽象为具有正交各向异性，且横观各向同性的材料[29-30]，该材料的结构沿轴向对称，忽略甘蔗节径向尺寸的细微差异，可将甘蔗尾茎试样视为类似圆柱体的几何形状；甘蔗尾茎分别由甘蔗皮、甘蔗芯和甘蔗节3种不同材质组成，各材质的纤维均沿轴向整齐排列，建立如图1所示坐标系，沿甘蔗尾茎轴向建立 Z轴，沿甘蔗尾茎径向分别建立X轴和Y轴。
1.3 甘蔗尾茎泊松比参数的确立
根据复合材料力学理论，甘蔗尾茎各组分泊松比可由9个工程弹性参数来表征，分别为甘蔗皮同性平面泊松比
$ \mu _{{}_{{\rm 1}XY}}$ ；甘蔗皮异性平面泊松比$\mu _{{}_{{\rm 1}XZ}} $ 、$ \mu _{{}_{{\rm 1}YZ}}$ ；甘蔗芯同性平面泊松比$\mu _{{}_{{\rm 2}XY}} $ ；甘蔗芯异性平面泊松比$ \mu _{{}_{{\rm 2}XZ}}$ 、$\mu _{{}_{{\rm 2}YZ}} $ ；甘蔗节轴向泊松比$ \mu _{{}_{{\rm 3}XY}}$ ；甘蔗节径向泊松比$ \mu _{{}_{{\rm 3}XZ}}$ 、$ \mu _{{}_{{\rm 3}YZ}}$ 。根据甘蔗几何形状的假定以及上述推论，甘蔗尾茎属于特殊的复合材料，正交各向异性，且甘蔗尾茎横截面具有横观各向同性的特点，其各组分泊松比参数满足式(1)[31]。
$$ \left\{ \begin{aligned} &{\mu _{{}_{iXZ}}} = {\mu _{{}_{iYZ}}}\\ &{\mu _{{}_{iXY}}} = \displaystyle\frac{{{E_{{}_{iX}}}}}{{2{G_{{}_{iXY}}}}} - 1 \end{aligned} \right. \text{，}$$ (1) 式中，E为弹性模量，MPa；G为剪切模量，MPa；X、Y和Z分别为3个坐标轴方向；i=1, 2, 3分别代表甘蔗皮、甘蔗芯和甘蔗节。
采用电测法，设计专门夹具来测量甘蔗尾茎各组分轴向拉伸泊松比，由试验得出甘蔗尾茎各组分横向应变(
$\varepsilon' $ )和纵向应变($\varepsilon $ )，根据式(2)计算出各组分轴向泊松比[32]。$$ \mu {\text{ = |}}\frac{\varepsilon' }{\varepsilon }{\text{|}}\text{。} $$ (2) 最后根据复合材料力学中关于正交各向异性材料的泊松比关系式(3)，结合甘蔗各组分径向拉伸弹性模量、各组分轴向拉伸弹性模量和各组分异性平面泊松比所得数据，确定甘蔗各组分径向泊松比参数。
$$ {\mu _{{}_{iXY}}} < \frac{{E_{iX}^2}}{{2\mu _{iYZ}^2E_{iZ}^2}}\;\;\;\;\text{。} $$ (3) 2. 力学试验
2.1 试验材料
WD-E精密型微控电子式万能试验机；YJ-4501A/SZ静态数字电阻应变仪；辅助测试工具包括砝码、游标卡尺、直尺、兆欧表、电烙铁、砂轮机等。
本试验根据甘蔗尾茎的结构特点设计专用夹具(图2)，该夹具主要由支座、压块、螺栓等组成。为使甘蔗各组分样本受力均匀，得到较为准确的试验数据，将甘蔗各组分样本制成规则矩形，并在压块和支座上加工有凹槽，上面有斜纹滚花，以增大摩擦。
试验材料为成熟期甘蔗粤糖159，采自广东省广前糖业发展有限公司前进农场，采样时间为2016年1月。将甘蔗尾梢相互重叠的青叶剥除干净，取蔗尾茎顶端生长点以下第5节和第6节制作试验样本。
2.2 方法
2.2.1 拉伸试验
将甘蔗第5节和第6节的甘蔗皮和甘蔗芯制成长60 mm、宽10 mm、厚2 mm的样本，两端夹持尺寸分别为20和15 mm，中间留10 mm标距粘贴应变片，应变片粘贴方法如图3a所示；将甘蔗第5节和第6节中间的蔗节部位制成如图3b所示试样，两端蔗身用于夹具夹持，蔗节部位开槽后长宽高均为10 mm，用于粘贴应变片，沿着试样的轴向进行拉伸试验。试验时为避免拉伸速率太慢，出现滑移现象并造成数据错误，经过测试对比，本试验拉伸速率设定为5 mm·min–1。
为避免甘蔗试样在夹紧的过程中发生滑移和破坏，试验样本在装夹部位用纱布、胶水和橡胶手套固紧。为保证拉伸试验的准确性，防止试样表面不平而引起的应力集中，在拉伸试验前用砂轮机将试验样本装夹表面磨平，并用石蜡抹平试样装夹表面的凹凸处。
2.2.2 泊松比的提取
为了获得较为精确的泊松比数值，每个试样在进行拉伸试验前都要用兆欧表进行绝缘检查，防止应变片与样本短路。YJ-4501A/SZ静态数字电阻应变仪采用半桥测量电路[33]，外接2个固定电阻，电阻应变片的贴片及接桥方式如图4所示。试验时用2台电阻应变仪分别测量试样的横向应变和轴向应变，运用公式算出试样泊松比，防止切换测量通道时引起的测量误差。同一部位的样本做20个试样，去除试验过程中不可抗力因素引起的失败样品，最后取15个成功试样进行统计分析，得出不同组分的泊松比参数值。统计公式如下：
$$ {\mu _{{}_0}} = \frac{1}{n}\sum\limits_{i = 1}^n {(\frac{{\varepsilon {'_i}}}{{{\varepsilon _i}}}} )\;\;\;\;\;\;, $$ (4) 式中，μ0为泊松比；n为取样数；εi'为横向应变；εi为轴向应变。
2.3 试验结果与分析
2.3.1 轴向拉伸泊松比
对甘蔗皮、芯、节的15组试样进行轴向拉伸试验，为保证结果的准确性，分别在10、20、30、40、50 g载荷下测量样本泊松比，最后取其平均值，第5、6节甘蔗皮的泊松比为0.233±0.073和0.238±0.051；第5、6节甘蔗芯的泊松比为0.271±0.045和0.289±0.049；第5、6节间甘蔗节的泊松比为0.344±0.086。
2.3.2 甘蔗尾茎各组分泊松比
采用WD-E精密型微控电子式万能试验机分别对甘蔗第5、6节各组分的15个试样进行轴向和径向拉伸试验，并在万能试验机上得到试样的载荷–位移曲线。为消除尺寸的影响，将试验机上的载荷–位移曲线转化为应力–应变曲线，如图5所示。
![]()
图 5 轴向和径向拉伸应力–应变曲线a~e为轴向拉伸；f~j为径向拉伸。Figure 5. The stress-strain curves of axial and radial tension载荷、位移、应力和应变相互转化关系见以下公式：
$$ \sigma = \frac{P}{A}, $$ (5) $$ E = \frac{{{\rm d}\sigma }}{{{\rm d}\varepsilon }}, $$ (6) $$ \varepsilon = \frac{{\Delta L}}{L}, $$ (7) $$ {\sigma _{\max }} = \frac{{{P_{\max }}}}{A}, $$ (8) 式(5)~(8)中，σ为应力，MPa；P为拉伸载荷，N；A为断面面积，mm2；ε为应变；∆L为伸长量，mm；L为标距原始长度，mm；σmax为抗拉强度，MPa；Pmax为最大拉伸载荷，N。
运用SPSS软件对上述应力–应变曲线进行最小二乘法拟合分析，求取甘蔗尾茎各组分弹性模量平均值得到如下结果：第5节甘蔗皮和甘蔗芯的径向弹性模量分别为6.90和14.30 MPa，相应的轴向弹性模量分别为124.97和89.55 MPa；第6节甘蔗皮和甘蔗芯的径向弹性模量分别为8.14和15.95 MPa，相应的轴向弹性模量分别为705.93和82.66 MPa；第5、6节间甘蔗节的径向弹性模量为63.87 MPa，轴向弹性模量为237.58 MPa。
从图5可以看出，部分应力-应变曲线差异较大，分析其原因主要是：甘蔗各样本结构尺寸差异大，要制作统一的样本很困难；各样本含水率对试验的结果影响很大，样品含水率的差异使应力-应变曲线差异较大；甘蔗自身组织结构的差异，导致试验结果存在较大差异，已有研究证明蔗尾茎秆的力学性能与其微观组织结构有关[34]；试验加载速率也会对试验结果产生影响。
将甘蔗皮、甘蔗芯及蔗尾茎第5、6节各组分轴向弹性模量和径向弹性模量代入式(3)得各组分同性面的泊松比，第5节甘蔗同性面芯和皮泊松比分别小于0.174和0.028；第6节甘蔗同性面芯和皮泊松比分别小于0.223和0.001。甘蔗节同性面泊松比小于0.305。综合试验和计算所到的泊松比参数，结合式(1)，可得甘蔗尾茎各组分泊松比，如表1所示。
表 1 甘蔗尾茎各组分泊松比参数Table 1. Poisson’s ratio parameters of each part of sugarcane tail stalk组分 同性面泊松比(μXY) 异性面泊松比(μYZ) 异性面泊松比(μXZ) 第5节 第6节 第5节 第6节 第5节 第6节 甘蔗皮 <0.028 <0.001 0.233 0.238 0.233 0.238 甘蔗芯 <0.174 <0.223 0.271 0.289 0.271 0.289 甘蔗节 <0.305 0.344 0.344 运用SPSS软件对第5、6节的甘蔗皮、芯试样的同性面泊松比和异性面泊松比分别进行独立样本t检验，分析蔗尾顶端生长点以下5、6节各组分泊松比的差异显著性。分析结果表明：第5、6节的甘蔗皮和芯试样同性面泊松比，P<0.01，由此可以认为第5、6节的甘蔗皮和甘蔗芯的同性面泊松比差异均显著；第5、6节甘蔗皮和芯试样异性面泊松比，P分别为0.096和0.012，由此可以认为第5、6节的甘蔗皮和甘蔗芯的异性面泊松比差异均不显著。
分别对甘蔗第5、6节不同节位各组分对应参数进行单因素方差分析(表2)。从表2可以看出，甘蔗皮同性面泊松比的P=0<0.05，甘蔗芯同性面泊松比的P=0.041<0.05，表明在95%的置信区间内，第5、6节不同节位各组分的同性面泊松比参数有显著差异；甘蔗皮异性面泊松比的P=0.791>0.05，甘蔗芯异性面泊松比的P=0.324>0.05，表明在95%的置信区间内，第5、6节不同节位各组分的异性面泊松比参数差异不显著。分析结果表明，甘蔗尾茎各组分的泊松比参数不一致，差异显著；甘蔗芯部同性面泊松比参数比甘蔗皮部大0.145，在一定范围内，甘蔗从表皮到内部泊松比逐渐变大，说明在建立甘蔗尾茎材料模型的时候，甘蔗皮、芯应视为完全不同的材料，可简化为由皮、芯2种单向复合材料组合而成的复杂复合材料，蔗皮与芯部存在结合力[35]。第5、6节不同节位的各组分泊松比力学参数存在明显差异，证明将蔗尾模型简化为蔗尾部分和蔗身部分，蔗尾生长点以下第5节的力学参数作为蔗尾部分力学参数，第6节的力学参数作为蔗身部分的力学参数这一假设的正确性。
表 2 甘蔗节位对泊松比参数影响方差分析Table 2. Variance analysis of Poisson’s ratio parameters affected by different sugarcane nodes性状 差异来源 平方和 自由度 均方 F P1) 甘蔗皮同性面泊松比 组间 0.009 1 0.009 34.188 0.000** 组内 0.008 28 0.000 总和 0.017 29 甘蔗皮异性面泊松比 组间 0.000 1 0.000 0.072 0.791 组内 0.110 28 0.004 总和 0.110 29 甘蔗芯同性面泊松比 组间 0.024 1 0.024 4.142 0.041* 组内 0.164 28 0.006 总和 0.188 29 甘蔗芯异性面泊松比 组间 0.002 1 0.002 1.008 0.324 组内 0.063 28 0.002 总和 0.065 29 1) *和**分别表示在 0.05、0.01 水平差异显著。 3. 结论
运用复合材料力学理论假定甘蔗力学模型，采用电测法代替传统的比拟法测量分析获得甘蔗各组分试样的横纵向变形量，得到较为准确的甘蔗第5、6节各组分泊松比参数值。
试验得到甘蔗尾茎生长点以下第5节甘蔗皮同性面泊松比小于0.028，异性面泊松比为0.233；第5节甘蔗芯同性面泊松比小于0.174，异性面泊松比为0.271；甘蔗生长点以下第6节甘蔗皮同性面泊松比小于0.001，异性面泊松比为0.238；第6节甘蔗芯同性面泊松比小于0.223，异性面泊松比为0.289；甘蔗节同性面泊松比小于0.305，异性面泊松比为0.344。甘蔗尾茎内部不同部位的同性面泊松比有显著差异，异性面泊松比差异不大；甘蔗芯部同性面泊松比参数比甘蔗皮部大0.145。试验结果表明，在一定范围内，甘蔗从表皮到内部泊松比逐渐变大。甘蔗尾茎第5、6节不同节位的皮、芯同性面泊松比参数有显著差异；第5、6节不同节位的皮、芯异性面泊松比参数差异不显著。甘蔗第5节皮部同性面泊松比参数比第6节大0.027，第5节芯部同性面泊松比参数比第6节小0.050。本研究可为甘蔗断尾机构的设计及数学模型的建立进行动力学仿真提供理论基础。
致谢: 感谢江西农业大学省部共建猪遗传改良与养殖技术国家重点实验室给予的支持。 -
表 1 商品肉猪中VRTN基因频率分布情况1)
基因型 个数 基因型频率 Q基因频率 q基因频率 QQ 141 0.336 0.585 0.415 Qq 209 0.498 qq 70 0.166 1)Q等位基因为多肋有利基因，q等位基因为不利基因 
下载: 导出CSV
表 2 VRTN在商品肉猪中效应分析与肋骨数表型变异系数
基因型 样本数 肋骨平均数/对 标准差 P1) 变异系数/% qq 31 15.16 0.82 < 0.001*** 5 Qq 96 15.57 0.68 QQ 90 15.90 0.58 1)* * *表示在0.001水平显著相关
下载: 导出CSV
表 3 VRTN基因型及基因频率的分布1)
品种 样品数 基因型频率 基因频率 QQ Qq qq Q q W51 913 0.814(743) 0.182(166) 0.004(4) 0.905 0.095 W52 240 0.733(176) 0.221(53) 0.046(11) 0.844 0.156 1)Q等位基因为多肋有利基因，q等位基因为多肋不利基因；括号中的数字为样本数
下载: 导出CSV
表 4 VRTN基因型与W51长白选育性能的关系1)
表型 qq型 Qq型 QQ型 P2) 体长/cm 113.750±2.217 116.910±2.974 117.572±2.996 0.002** 校正30~100 kg日增重/g 939.765±18.105 897.273±75.746 879.760±71.459 0.006** 综合评分/分 8.450±0.100 8.557±0.393 8.647±0.382 0.016* 校正30~115 kg日增重/g 908.072±21.285 886.075±68.357 872.680±66.658 0.041* 终测有效乳头数/个 14.500±1.000 14.266±1.022 14.46±0.974 0.068 健仔数估计育种值 -1.070±1.365 0.020±1.005 0.050±0.965 0.069 繁殖指数 65.528±43.570 100.21±32.206 101.170±30.943 0.072 校正115 kg背膘估计育种值 -1.890±0.807 -0.876±0.835 -0.912±0.902 0.079 校正100 kg背膘估计育种值 -1.677±0.706 -0.794±0.728 -0.825±0.787 0.079 总仔数估计育种值 -1.251±1.549 -0.029±1.160 0.005±1.115 0.081 100 kg时估计瘦肉率/% 53.722±0.536 54.101±1.635 54.387±1.588 0.084 115 kg时估计瘦肉率/% 54.350±0.584 54.777±1.800 55.090±1.750 0.086 校正115 kg眼肌面积估计育种值 1.952±1.185 -0.212±2.070 -0.202±1.990 0.100 校正100 kg体重日龄/d 153.838±8.167 153.325±8.294 154.882±8.591 0.104 校正115 kg体重日龄/d 170.912±9.393 170.324±9.538 172.114±9.880 0.104 校正100 kg背膘厚/mm 12.383±0.619 11.529±1.766 11.272±1.836 0.132 校正115 kg背膘厚/mm 14.197±0.709 13.229±2.017 12.944±2.096 0.144 父系指数 132.135±9.702 123.909±11.223 125.276±11.279 0.167 校正115 kg体重日龄估计育种值 -3.998±2.401 -6.028±3.316 -6.353±3.156 0.176 校正100 kg眼肌面积/cm2 36.500±1.422 36.703±3.393 37.205±3.246 0.190 校正115 kg眼肌面积/cm2 39.718±1.547 39.937±3.692 40.482±3.532 0.190 校正100 kg体重日龄估计育种值 -3.230±2.089 -4.977±2.877 -5.248±2.747 0.191 校正100 kg眼肌面积估计育种值 -0.137±0.587 -1.487±1.676 -1.352±1.612 0.196 母系指数 102.168±23.419 113.246±14.242 114.429±15.596 0.198 外貌指数估计育种值 0.337±0.244 0.306±0.133 0.327±0.140 0.208 同窝仔猪数/头 14.000±0.000 12.108±2.617 12.377±2.717 0.238 活仔数估计育种值 -0.555±0.845 -0.081±0.715 -0.063±0.711 0.375 出生重/kg 1.525±0.472 1.604±0.296 1.633±0.282 0.385 选留整体评分估计育种值 0.031±0.017 0.065±0.058 0.067±0.057 0.420 校正30~100 kg日增重估计育种值 55.914±20.258 46.886±20.396 47.936±18.104 0.547 校正30~115 kg日增重估计育种值 27.232±20.690 35.399±20.334 36.203±19.395 0.595 1)表中数据为平均值±标准差，qq型样本数量为4头，Qq型样本数量为166头，QQ型样本数量为743头；2)*、* *分别表示0.05、0.01水平差异显著
下载: 导出CSV
表 5 VRTN基因型与W52生长性能的关系1)
表型 qq型 Qq型 QQ型 P 2) 校正115kg体重日龄估计育种值 -4.901±1.773 -4.655±2.810 -5.804±2.612 0.016* 校正100 kg体重日龄估计育种值 -4.037±1.557 -3.855±2.437 -4.832±2.275 0.019* 校正100 kg眼肌面积估计育种值 -0.056±1.604 0.014±1.251 0.651±1.460 0.053 校正115 kg体重日龄/d 168.300±4.856 172.033±9.832 168.518±10.475 0.085 校正100 kg体重日龄/d 151.566±4.221 154.811±8.551 151.755±9.109 0.085 校正30~100 kg日增重估计育种值 25.779±15.665 28.86±20.541 34.018±17.065 0.187 校正100 kg背膘估计育种值 -0.685±1.231 -1.074±1.086 -0.816±0.967 0.221 校正115 kg背膘估计育种值 -0.973±1.437 -1.461±1.251 -1.204±1.077 0.251 校正30~115 kg日增重估计育种值 46.243±16.538 48.364±23.411 53.189±21.934 0.264 体长/cm 124.273±5.101 126.208±4.638 126.540±4.741 0.296 父系指数 127.617±12.604 131.98±12.128 132.848±11.527 0.341 终测有效乳头数/个 14.909±0.701 14.925±1.035 15.125±0.978 0.366 同窝仔猪数/头 13.455±2.544 13.151±2.042 12.727±2.582 0.393 出生重/kg 1.600±0.195 1.653±0.269 1.690±0.282 0.442 100 kg时估计瘦肉率/% 52.666±1.963 53.389±1.985 53.324±1.729 0.500 115 kg时估计瘦肉率/% 53.176±2.200 53.975±2.201 53.902±1.911 0.502 校正30~115 kg日增重/g 922.154±80.807 896.160±83.767 910.771±91.245 0.503 母系指数 119.022±14.065 122.837±13.458 123.591±13.271 0.531 综合评分/分 8.000±0.742 7.972±0.710 8.088±0.696 0.549 校正100 kg眼肌面积/cm2 33.783±3.245 35.13±3.723 34.980±3.550 0.550 校正115 kg眼肌面积/cm2 36.761±3.532 38.225±4.051 38.061±3.863 0.550 校正30~100 kg日增重/g 928.909±98.072 900.863±93.354 914.714±101.47 0.575 选留整体评分估计育种值 0.111±0.107 0.159±0.161 0.162±0.163 0.593 校正115 kg眼肌面积估计育种值 0.557±1.608 1.214±1.949 1.095±2.018 0.609 活仔数估计育种值 0.350±0.492 0.268±0.381 0.256±0.418 0.763 校正115 kg背膘厚/mm 14.081±3.807 13.574±2.830 13.760±2.414 0.811 校正100 kg背膘厚/mm 12.206±3.301 11.773±2.458 11.933±2.094 0.816 外貌指数估计育种值 0.454±0.207 0.465±0.220 0.483±0.238 0.830 健仔数估计育种值 0.274±0.917 0.380±0.860 0.395±0.786 0.890 总仔数估计育种值 0.283±1.043 0.377±0.972 0.410±0.889 0.890 繁殖指数 108.512±29.285 111.662±27.417 112.276±25.063 0.891 1)表中数据为平均值±标准差，qq型样本数量为11头，Qq型样本数量为53头，QQ型样本数量为166头；2)*表示0.05水平差异显著
下载: 导出CSV
-
[1] FREEMAN V A. Ariations in the number of vertebrae of swine[J].Hered, 1939, 30:61-64. https://www.researchgate.net/publication/31305538_Variations_in_the_number_of_vertebrae_of_swine
[2] SHAW A W. A method of determining the variations in the vertebral column of the live pig[J]. Science Agriculture, 1930, 10:90-695. http://cn.bing.com/academic/profile?id=2fd4293a46e4fbff5638923059857425&encoded=0&v=paper_preview&mkt=zh-cn
[3] MIKAWA S, MOROZUMI T, SHIMANUKI S.et al. Fine mapping of a swine quantitative trait locus for number of vertebrae and analysis of an orphan nuclear receptor, germ cell nuclear factor(NR6A1)[J].Genome Research, 2007, 17:586-593. http://d.old.wanfangdata.com.cn/OAPaper/oai_pubmedcentral.nih.gov_1855175
[4] FRIETSON G. Why do almost all mammals have seven cervical vertebrae? Developmental constraints, hox genes, and cancer[J].Journal of experimental zoology, 1999, 285:19-26. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ021393319/
[5] KING J W B, ROBERTS R C. Carcass length in the bacon pig:Its association with vertebrae numbers and prediction from radiographs of the young pig[J]. Animal Production, 1960, 2:59-65. http://cn.bing.com/academic/profile?id=7b09c7ac5be4284623e2493a59ceb9f0&encoded=0&v=paper_preview&mkt=zh-cn
[6] BERGE S. Genetical researches on the number of vertebrae in the pig[J]. Journal of Animal Science, 1948, 7:233-238. http://cn.bing.com/academic/profile?id=7c4330533e8d21b0e5ea8b852bb9f8c8&encoded=0&v=paper_preview&mkt=zh-cn
[7] 阮桂凤.猪7号染色体脊椎数QTL的精细定位及候选基因HOX10、PROX2和FOS的SNPs搜寻[D].南昌: 江西农业大学, 2010. [8] 杨广成, 猪复杂性状的遗传解析——以脊椎数和初情期为例[D].南昌: 江西农业大学, 2008. [9] FAN Y, XING Y Y, HUANG L S, et al. A further look at porcine chromosome 7 reveals VRTN variants associated with vertebral number in Chinese and Western pigs[J]. Plos One, 2013, 8(4):e62534.
[10] MIKAWA S, HAYASHI T, NII M, et al. Two quantitative trait loci on Sus scrofa chromosomes 1 and 7 affecting the number of vertebrae[J]. Journal of Animal Science, 2005, 83:2247-2254. http://cn.bing.com/academic/profile?id=6f16fd3e6866ae8ff984c2a3f908440b&encoded=0&v=paper_preview&mkt=zh-cn
[11] YANG J, HUANG L S, YANG M, et al. Possible introgression of the VRTN mutation increasing vertebral number, carcass length and teat number form Chinese pigs into European pigs[J]. Scientific Reports, 2016, 6:19240. https://www.ncbi.nlm.nih.gov/pubmed/26781738
[12] 杨明. MUC13、FUT1基因在2个种猪核心群中的分子标记辅助选择研究[D].广州: 华南农业大学, 2015, 36(6): 1-8. [13] 杨明.氟烷基因在皮特兰猪群体中的分子标记辅助选择研究[D].广州: 华南农业大学, 2018, 39(6): 1-4. -
期刊类型引用(1)
1. 许凯，马英，陈三雄，骆金初，黄伟城，郭微，王龙远. 广东德庆县香山森林公园黄牛木种群特征研究. 安徽农业科学. 2022(13): 125-129 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 744
- HTML全文浏览量: 10
- PDF下载量: 1354
- 被引次数: 1
