山羊褪黑激素受体1a基因hinfⅰ酶切多态性及其与产羔数性状的关联分析-学位论文.doc

华中农业大学本科毕业论文（或设计）

本科毕业论文

题 目：
山羊褪黑激素受体1A基因HinfⅠ酶切

多态性及其与产羔数性状的关联分析
姓 名：
贾 征
学 号
040801030
专 业：
动 物 科 学
指导教师：
杨 利 国
职 称
教 授

中国·武汉
二○○八年六月

分类号 密级

华中农业大学本科毕业论文

山羊褪黑激素受体1A基因HinfⅠ酶切多态性及其与产羔数性状的关联分析
HinfⅠ-RFLP of Melatonin receptor 1 A(MTNR1A) gene and the association with the litter size traits in goat

学生姓名：贾 征
学生学号：040801030
学生专业：动物科学
指导教师：杨利国 教 授

华中农业大学动物科学学院
二○○八年六月
目 录
摘 要 4
ABSTRACT 5
前 言 6
1.材料与方法 7
1.1 材料 7
1.1.1 试验动物 7
1.1.2 实验材料 7
1.1.3主要试剂的来源 7
1.1.4主要仪器设备 8
1.1.5主要试剂的配置 8
1.1.6主要分子生物学软件 10
1.2 方法 10
1.2.1.山羊基因组DNA的提取 10
1.2.2.引物设计及其序列 10
1.2.3.PCR扩增及其电泳检测 11
1.2.4. 限制性内切酶消化扩增产物以及电泳检测 13
1.2.5.统计分析 14
2 结果与分析 15
2.1 PCR扩增结果 15
2.2酶切结果 15
2.3 褪黑激素基因频率和基因型频率分析结果 16
2.4 MTNR1A基因型与波尔山羊头胎产羔数性状的关联分析 17
2.5 MTNR1A基因型与经产产羔数性状的关联分析 17
3 讨论 18
3．1 不同山羊品种基因型频率及基因频率差异 18
3．2 MTNR1A不同基因型对波尔山羊头胎和经产羔数的影响 18
参考文献 19
致 谢 20

山羊褪黑激素受体1A基因HINFⅠ酶切多态性及其与产羔数性状的关联分析

摘 要
本研究以褪黑激素受体1A 基因为候选基因，分析该基因多态性对山羊产羔数性状的影响。采用限制性片段长度多态性（Restriction Fragment Length polymorphism， PCR-RFLP）方法，检测了65只波尔山羊和51马头山羊褪黑激素受体1A基因HinfⅠ酶切多态性及其与产羔数性状的关联分析。结果表明，褪黑激素受体1A基因在波尔山羊中存在HinfⅠ酶切多态性，而在马头山羊样本中则无多态性。在波尔山羊中，检测到AA、GG两种基因型，未检测到AG型；其中AA基因型个体在出现频率最高，为0.554，GG型为0.446；A、G等位基因频率分别为0.554和0.446。性状关联分析结果表明：波尔山羊的AA基因纯合型头胎平均产羔数比GG基因杂合型多0.07只，加性效应值为0.035，显性效应值为-1.42；经产羔数的AA型比GG型多0.09只；加性效应值为0.046，显性效应值为-1.84。A等位基因可能与波尔山羊产羔数呈正相关，该位点主要以加性方式起作用。

关键词山羊；繁殖性状；褪黑激素受体；PCR—RFLP

HinfⅠ-RFLP of Melatonin receptor 1 A（MTNR1A） gene and the association with the litter size traits in goat

Abstract
In the present study， Melatonin receptor 1 A gene （MTNR1A） was analysed to find the assocation with litter size in goat as a candidate gene. HinfⅠ-RFLP （Restriction Fragment Length polymorphism， RFLP） of MTNR1A was found in 65 Boer goats and 51 goats Ma Tou. No polymorphism was detected in Ma Tau goats . Two genotypes AA， GG were found in Boer goats. The frequency of AA genotype was higher than GG genotype， which were 0.554 and 0.446， respectively. The allele frequencies for A and G were 0.554 and 0.446. The results of traits association analysis showed that individuals with AA genotype had additional 0.07 average litter size at first kidding in Boer Goat. The additive effect value was 0.035， while the dominant effects value -1.42. Moreover， AA genotype had additional 0.09 average litter size for multiparity， the additive and dominant effects value were 0.046 and -1.84， respectively. Allele A was possiblely positively correlated with litter size in Boer goats. The additive effect was dominant at the site.

Key words: Goat； Reproductive traits； MTNR1A； PCR-RFLP

前 言
褪黑激素主要是由松果体合成的，其它合成部位还有:视网膜、眼眶腔的副泪腺、唾液腺、肠的嗜铬细胞及红细胞，松果体把对外界光照形成的神经信息经下丘脑视交叉上核转变成体液信息—MEL。
1.褪黑激素生理功能、作用机制；
褪黑激素对生殖系统作用的机理主要是通过光作用于视网膜，然后通过视网膜- 下视丘至视交叉上核，再经过一系列复杂的神经交换至颈上神经节，其节后交感神经作用于松果体。松果体分泌褪黑激素，通过下丘脑—垂体性腺轴影响生殖系统褪黑激素对生殖系统的作用因动物种类、生理状态不同而表现出抑制、促进或无作用的多重性（焦传珍，2005）。多项研究结果表明，褪黑激素虽然对牛、鼠类、禽等动物和人的生殖系统的发育有抑制作用，但是对绵羊、山羊、鹿等动物则明显表现出促进作用。
2．褪黑激素受体类型、作用机制；
褪黑激素受体属于G蛋白耦联受体家族。根据其分子生物学方法又可分为MTNR1A、MTNR1B、MTNR1C三种亚型（张凯等，2006）。目前在哺乳动物中还没有发现MTNR1C受体，而具有明显季节性繁殖特征的动物仅在垂体结节部发现MTNR1A受体，提示垂体结节部可能是褪黑激素繁殖效应的作用位点。
3．褪黑激素受体分子标记在其它动物中的研究进展
目前国内外对绵羊M TNR1A 基因已经有所研究已经很多了。Messer等（1997）参考文献中未列出
和Notter等（2003）用M n lⅠ和R saⅠ酶分别对白面杂种羊、萨福克羊、特克塞尔羊、柯泊华斯羊和合成系绵羊群体M TNR1A 基因外显子2 的824 bp扩增产物进行消化，发现M n lⅠ在5 个绵羊群体中都存在286 bp和236 bp多态片段， R saⅠ都存在295 bp 和290 bp 多态片段。Pelletier等（2000）用M n lⅠ限制性内切酶消化2组发情行为有明显差异（常年发情和季节性发情）的Merinos d′Arles母绵羊以及季节性发情明显的Ile2de2France 母绵羊M TNR1A基因外显子2的824 bp扩增产物，发现存在303 bp和236 bp 多态片段。季从亮等（2003）和Chu等（2003）用M n lⅠ和R saⅠ对常年发情的小尾寒羊和湖羊以及季节性发情的多赛特羊和萨福克羊共146只绵羊M TNR1A 基因外显子2的824 bp 扩增产物进行PCR-RFLP多态性分析，发现M n lⅠ在4个绵羊品种中都存在303 bp和236 bp多态片段， R saⅠ都以上迹象表明，这是直接从网上下载的。
存在290 bp和267 bp 多态片段。在绵羊中的研究结果表明，该基因对绵羊繁殖性能存在影响。但是，对于山羊，只有对产绒或营养调控的研究，目前尚未发现有对山羊繁殖性能影响的相关报道。
4．本实验目的意义
本实验室滑国华等（2006）在山羊MTNR1A中共发现7个SNPs，对其中4个SNPs分析结果显示：MTNR1A不同基因型对山羊产羔数性状有一定的影响。为进一步研究该基因对山羊产羔数性状的影响，本试验选取G493A突变位点，采用引入酶切位点方法，即限制性片段长度多态性Forced PCR-RFLP（Polymerase Chain Reaction-Restriction Fragment Length Polymorphsim），分析该位点在马头山羊和波尔山羊中的分布情况，并分析该位点对波尔山羊产羔数影响，为山羊标记辅助选择提供理论依据与滑国华的区别？更确切地说，创新点？或特色？
。

1.材料与方法
1.1 材料
1.1.1 试验动物
马头山羊51头，采自湖北省十堰市；波尔山羊65头，采自湖北省宜都市波尔山羊种羊场。颈静脉采血10ml/只，EDTA抗凝，-20℃保存。
1.1.2 实验材料
1.5ml Eppendorf管、双面板、微量移液枪（1000μl、200μl、40μl、10μl、2.5μl）及枪头、 一次性PE塑料手套。
1.1.3主要试剂的来源
TagDNA聚合酶购自上海生工生物工程技术有限公司；
DNA限制性内切酶HinfⅠ购自晶美生物工程有限公司；
分子量标准100bp DNA Ladder和PBR322DNA/MspI购自上海生工生物工程技术有限公司；
三羟甲基氨基甲烷（Tris）、苯酚、蛋白酶K、十二烷基硫酸钠（SDS）、甲酰胺、丙烯酰胺（Acrylamide，Acr）、N－N’－亚甲基丙烯酰胺（Bisacrylamide）、琼脂糖、溴化乙锭（EB）、琼脂糖、乙二胺乙二酸二钠盐（EDTANa2）等均购自武汉凌飞科技有限公司。

1.1.4主要仪器设备
表1 主要仪器设备
Table 1 Laboratory equipment and instruments
仪器设备名称
型号
生产厂家
电热恒温水浴锅
GD120
英国Grant公司
台式离心机
TGL-16G
上海安亭科学仪器厂
梯度PCR仪
T1 Thermo cyler
德国Biometra公司
冷冻离心机
Centrfuge 5810
德国eppendorf公司
恒温摇床
HS80
中国科学院武汉科学仪器厂
自动双重纯水蒸馏器
SZ-93
上海亚荣生化仪器厂
紫外分析仪
UV-IV
北京市新科技应用研究所
电泳槽
DYY-III
北京六一仪器厂
稳压稳流电泳议
DYY-2
北京六一仪器厂
手提式高压灭菌锅
YXQ-SG46-280A
上海博讯实业有限公司医疗设备厂
超纯水仪
WaterPro
美国LABCONCO公司
移液器
1ml；200μl；40μl； 10μl；2.5μl
芬兰Labsystems公司

1.1.5主要试剂的配置
1.1.5.1用于血液样品采集的溶液的配制
1) 1M Tris.cl储备液：称取121.4g Tris碱溶于800ml双蒸水中，加浓盐酸调至pH=8.0，定容到1000ml。高压灭菌，室温保存；
2) 0.5M EDTA（pH=8.0）：称取186.1g Na2EDTA·2H2O溶于800ml双蒸水中，加入NaOH（约20g）调至pH=8.0，定容到1000ml。高压灭菌，室温保存；
3) 10%SDS：称取50g SDS加热溶于400ml双蒸水中，定容至500ml。高压灭菌，室温保存；
4) 抗凝剂（0.2MEDTA）：量取200ml 0.5M EDTA加入双蒸水定容至500ml。
1.1.5.2 用于基因组DNA提取的溶液的配制
1) 蛋白酶K：称取200mg蛋白酶K溶于10ml双蒸水，以1ml分装，-20℃冷冻保存；
2) 氯仿/异戊醇（24:1）：量取480ml氯仿，加入20ml异戊醇，混匀；
3) 3M NaAc：称取NaAc·3H2O溶于400ml双蒸水中，用HAc调至pH=5.2，加水定容至500ml。高压灭菌，室温保存；
4) TE缓冲液100mM Tris·cl（pH＝8.0），1mM EDTA（pH＝8.0）。
1.1.5.3 用于琼脂糖凝胶电泳及检测溶液的配制
1) 50×TAE储备液： 称取242g Tris碱溶于750ml双蒸水中，加入57.1ml HAc、100ml 0.5M EDTA（pH＝8.0），加水定容至1000ml；
2) 5×TBE储备液：称取54g Tris碱和27.5g硼酸溶于750ml双蒸水中，加入20ml 0.5M EDTA（pH＝8.0），加水定容至1000ml；
3) 6×上样缓冲液：0.25%溴酚兰，0.25%二甲苯青，15%Ficoll聚蔗糖；
4) 溴化乙锭（EB，10mg/ml）:100ml双蒸水中加入0.1g EB，搅拌使完全溶解，转入棕色试剂瓶，锡箔包裹。
1.1.5.4 用于聚丙烯酰胺凝胶电泳和银染的溶液的配制
1) 30% 丙烯酰胺（丙烯酰胺:甲叉双丙烯酰胺＝29:1）:700ml双蒸水中加入290g丙烯酰胺，10g甲叉双丙烯酰胺，搅拌使完全溶解，加水定容至1000ml，4℃保存；
2) 10% 过硫酸铵：称取1g过硫酸铵溶于8ml双蒸水，定容至10ml；
3) 10% 乙醇：量取50ml无水乙醇溶于450ml水中，混匀；
4) 1% HNO3：量取15ml浓HNO3溶于985ml双蒸水中，混匀；
5) 0.2% AgNO3：称取1g AgNO3溶于500ml双蒸水中，置于棕色试剂瓶；
6) 3% Na2CO3（0.05%甲醛）：称取30g无水Na2CO3溶于900ml双蒸水中，加入1.5ml甲醛，定容至1000ml；
7) 3% HAc：量取30ml冰HAc溶于970ml双蒸水，混匀。
1.2 方法
1.2.1.山羊基因组DNA的提取
从山羊颈静脉采血10ml，放入装有抗凝剂（0.5mol/L EDTANa2）的离心管中，6000r/min离心15min，弃上清，吸取白细胞沉淀（参杂少量红细胞）于另一离心管中；加入两倍体积双蒸馏水（灭菌），摇匀沉淀，静置10min～15min，破碎红细胞；6000r/min离心15min，轻轻弃去上清，留白细胞沉淀；加入1×SET 5ml，蛋白酶K（10mg/ml）0.1ml至终浓度200ng/ml，10% SDS 250μl至终浓度0.5%，55℃水浴消化过夜；加入等体积饱和平衡酚轻轻摇匀，12000r/min离心15min，轻轻吸取上层水相于另一离心管中，弃下相苯酚；重复上步；加入等体积苯酚/氯仿/异戊醇（25:24:1），轻轻摇动10min，12000r/min离心15min，轻轻吸取上层水相于另一离心管中，弃下层有机相；加入等体积氯仿/异戊醇（24:1），轻轻摇动10min，12000r/min离心15min，轻轻吸取上层水相于另一离心管中，弃下层有机相；加入两倍体积预冷无水乙醇，轻轻转动管子，DNA呈絮状析出，用1ml吸头（灭菌）吸出DNA沉淀，放入1.5ml离心管中，加入预冷无水乙醇洗涤一次，离心弃去无水乙醇，再用75%乙醇漂洗2次，小心吸去乙醇，待DNA沉淀自然干燥后，加入适量TE（100μl～300μl）溶解。
1.2.2.引物设计及其序列
引物参考滑国华等（2008）进行设计，序列如下：
上游引物: 5’ACGACCCGAGGATCTATTCC 3’
下游引物: 5’ATATAAGTCCTGGGGCTTCAGATT 3’，
在下游引物中引入错配碱基A，当突变位点为G（即反义链为C）时，可形成HinfⅠ（GATTC）酶切位点；当突变位点为A（即反义链为T）时，该酶切位点缺失，由此可进行基因型判定。
1.2.3.PCR扩增及其电泳检测
1.2.3.1 PCR扩增反应体系（表2）
表2 PCR扩增反应体系（20μl）
Table 2 PCR amplification of the response system （20μl）
反应体系
response system
反应体积V（μl）
response volume（μl）
灭菌蒸馏水
14.6
10×Buffer
2
Mg2+
1.2
dNTP
0.2
上游引物
0.4
下游引物
0.4
Taq DNA聚合酶（5 U/ml）
0.2
模板DNA
1
Total
20

1.2.3.2 退火温度的优化
利用PCR机器上都有一个Gradient功能，在一次PCR过程中，对机器中不同位点的反应采用不同的退火温度，一般在第一次实验中，采用引物溶点上下10度设置退火温度梯度。每个温度梯度一个反应，最后电泳检测并比较各个退火温度的条带效果，得到最好结果。本试验中退火温度设定为55－65℃，结果表明，退火温度60℃时，无非特异性带，目的条带最为清晰。
1.2.3.3 PCR反应程序（表3）
表3 PCR反应程序
Table 3 PCR amplification of the response term
步骤
steps
PCR反应条件
PCR reaction conditions
温度
temperature
时间
time
1
预变性
95℃
5min
2C
变性
95℃
30sec
3C
退火
60℃
30sec
4C
延伸
72℃
30sec
5
延伸
72℃
7min
6
保温
16℃
10min
第二步至第四步循环35次。

1.2.3.4 琼脂糖凝胶电泳检测
PCR扩增产物用1.5%琼脂糖凝胶电泳检测， EB染色，利用BIORAD凝胶成相系统凝胶成像系统进行成像分析，记录结果并拍照。
1.2.3.4.1 琼脂糖凝胶电泳方法
（1） 灌胶
将透明胶带粘于灌胶模槽的两端，制成灌胶模。放入梳子。将模放于水平台上。取一烧杯，放入：琼脂糖0.75 g； 1×TAE 50 ml。混匀，加热，至琼脂糖彻底溶解。待稍冷后，加入5ul饱和的溴乙锭，轻晃混匀，然后轻轻倒入模子。用吸头吸去表面小泡。待冷却成胶。
（2）上样电泳
轻轻拔出梳子，将胶两头的透明胶带揭去，移入电泳槽，点样孔一侧靠近操作者，注入1×TAE缓冲液，浸没凝胶。使用微量可调采样器取3μl DNA充分混匀于适量溴酚兰，点样于1.5％琼脂糖凝胶孔内，并点3μl 600bp MarkerⅠ作为参照。点样时，加样头尖端插入上样孔内1/3高度，轻轻将样品推入，避免刺入前后凝胶和刺穿凝胶底部。使比重较大的样品自然沉入上样孔底。接上电源：阳极远离上样孔，阴极靠近上样孔（DNA向阳极移动）。开启电源，电压调至100V。20～30min后检查。
1.2.3.4.2 成像分析
利用BIORAD凝胶成相系统凝胶成像系统进行成像分析，记录结果并拍照。
1.2.4. 限制性内切酶消化扩增产物以及电泳检测
1.2.4.1 酶切反应体系（表4）
表4 酶切反应体系
Table 4 PCR and the amount of the components
反应体系
response system
反应体积V（μl）
response volume（μl）
灭菌蒸馏水
4.7
10×Buffer
1
HinfⅠ限制性内切酶
0.3
PCR产物
4
Total
10
将上述酶切体系置37℃水浴锅或恒温箱中酶切过夜。
1.2.4.2非变性聚丙烯酰胺凝胶电泳检测
1.2.4.2.1聚丙烯酰胺凝胶电泳操作步骤
不同的实验需要优化聚丙烯酞胺凝胶的交联度和浓度，才能达到最好的实验效果，本实验采取一下方法：
1. 在小烧杯中加入18.43m1双蒸水，9.33ml 30%丙烯酞胺（29:1），7ml 的5×TBE混匀；
2. 清洗制胶用玻璃板，烘干后用透明胶封好底边和侧边，固定在制胶用的有机玻璃板上；
3. 向小烧杯中加入245ul的10%AP，13ul TEMED混匀后迅速灌胶；
4. 当胶灌至离玻璃板上边缘0.5cm时停止灌胶，插入梳子，注意梳齿下不要有气泡，室温聚合1小时；
5. 凝胶聚合好后，拔出梳子，注射器吸出梳齿孔中未聚合的液体；
6. 取下玻璃板，撕下胶线，固定在电泳槽上，倒入电泳缓冲液1×TBE；
7. 在PCR产物中加入1/6体积上样缓冲液，用微量进样器取5u1酶切样点样，同时点分子量标记物PBR322/Mspl；
8. 上样完毕后，打开电泳议以80V电泳20min，再调到120V电泳4～5h。
1.2.4.2.2聚丙烯酰胺凝胶染色步骤
1. 拆下电泳装置，剥下聚丙烯酰胺凝胶，放到一个塑料盘内，用双蒸水冲漂洗2次；
2. 加入500ml 10％乙醇固定液，在摇床上固定8～10 min，弃去反应液；
3. 加入500ml 1％HNO3氧化，在摇床上轻摇8～10 min，弃去反应液，用双蒸水漂洗2次；
4. 倒入500ml 0.2％AgNO3银染液，染色15～30 min，弃去反应液，用双蒸水漂洗2次；
5. 加入500ml 3％Na2CO3显色液，手摇显色至出现清晰的银染带，迅速弃去反应液；
6. 用500ml 3％HAC停止显色；
7. 将已经染好的聚丙烯酰胺凝胶置于透射白光板上，观察酶切反应结果、记录并拍照。
1.2.5.统计分析
基因频率与基因型频率使用软件SPSS10.0进行统计分析，不同基因型对产羔数的最小二乘均值分析使用SAS（8.0版）GLM（General Linear Model）软件包进行，模型如下：
Yijkl=μ+Gi+Pi+eijk
Yijkl为产羔数，μ为群体均值，Gi为基因型效应，Pi为胎次效应，eijk为随机效应。
基因效应分析：
加性效应=（GG-AA）/2
显性效应=AG-（AA+GG）/2

2 结果与分析
2.1 PCR扩增结果
扩增波尔山羊和马头山羊基因组，产物用2.0%琼脂糖凝胶电泳检测，得到的片段与预期片断大小一致（197bp），条带清晰（图1），特异性高，满足后续分析要求。

图1 PCR扩增结果
注：M， MarherⅠ 100bp DNA Ladder；1、2、3为波尔山羊DNA扩增产物；4、5、6为马头山羊DNA扩增产物
Fig. 1 Agarose gel electrophoresis of PCR product.
M，1-3，4-6: 100bp DNA Ladder， PCR product of ?goat， PCR product of ?goat

2.2酶切结果
利用限制性内切酶HinfⅠ对PCR扩增产物进行酶切，使用8%的聚丙烯酰胺凝胶（29：1）电泳检测酶切产物。
在197 bp的扩增片段在118bp处存在另一HinfⅠ酶切位点，酶切后产生118 bp和79bp片段，其中波尔山羊在118处的酶切位点存在多态性，酶切位点的存在与缺失产生2 种基因型，分别命名为AA和GG:AA （ 118bp /118 bp ） 、GG（118bp /79bp ） ，酶切图谱见图2 ；而马头山羊在118处的酶切位点则无多态性，基因型均为AA型。酶切图谱见图3 。

图2 波尔山羊酶切图谱
注：M， pBR322/Msp 1 Marher；1～13为酶切产物条带
Fig.2 genotypes of the M TNR1A fragment defined by PCR-RFLP

图3 马头山羊酶切图谱
注：M， pBR322/Msp 1 Marher1～8为酶切产物条带
Fig.3 genotypes of the M TNR1A fragment defined by PCR-RFLP

2.3 褪黑激素基因频率和基因型频率分析结果
分析两种山羊品种的褪黑激素基因型频率发现，在该酶切位点，马头山羊只出现AA基因型，波尔山羊虽然出现两种基因型，但AA型频率高于GG型。在所检测的两个群体中，A等位基因是优势等位基因。
表6 山羊MTNR1A基因HinfⅠ酶切基因型和基因频率分布表
Table 6 Genotype and allele frequencies of HinfⅠ-RFLP for the Goat M TNR1A gene

品种
Breed
数量
N
基因型
Genetype
基因型频率
Genetype frequency
等位基因频率
Allele frequency
波尔山羊
65
AA
0.554
A 0.554
GG
0.446
G 0.446
马头山羊
51
AA
1
A 1
2.4 MTNR1A基因型与波尔山羊头胎产羔数性状的关联分析
由表7可知，波尔山羊的AA基因纯合型头胎平均产羔数比GG基因杂合型多0.07只，加性效应值为0.035，显性效应值为-1.42。研究表明：波尔山羊的A等位基因可能与山羊头胎产羔数呈正相关，该位点主要以加性方式起作用。
表7 波尔山羊MTNR1A基因型对头胎产羔数性状的最小二乘均值及基因效应分析
Table 7 Least squares means and allele substitution effects for litter size at first kidding in boar goat
品种
Breed
数量
N
基因型
Genetype
最小二乘均值±标准误
Least squares means±standard error

加性效应
additive effec
显性效应
dominant effects
波尔山羊
65
AA
1.4545±0.5222
0.035
-1.42
GG
1.3846±0.5064

2.5 MTNR1A基因型与经产产羔数性状的关联分析
由表8可知，波尔山羊的AA型经产羔数比GG型多0.09只；加性效应值为0.046，显性效应值为-1.84。结果表明：A等位基因可能与波尔山羊产羔数呈正相关，该位点主要以加性方式起作用。
表8 波尔山羊MTNR1A基因型对经产产羔数性状的最小二乘均值及基因效应分析
Table 8 Least squares means and allele substitution effects for litter size in multiparous boar goat
基因型
Genetype
最小二乘均值±标准误
Least squares means±standard error
加性效应
additive effects
显性效应
dominant effects
AA
1.8827±0.5871
0.046
-1.84
GG
1.7905±0.5988

3 讨论
3．1 不同山羊品种基因型频率及基因频率差异
利用荧光原位杂交方法和遗传连锁分析方法，褪黑激素受体1A已经在绵羊等多种动物中得到精细定位。并有多个研究结果表明，褪黑激素受体1A基因与绵羊的繁殖性状有着紧密的关联。在本研究所检测的两个山羊群体中，基因型频率及基因频率存在较大差异。本研究结果显示：马头山羊是AA基因型纯合群体，而波尔山羊中存在AA和GG两种基因型。这可能与两个品种的育成史有关。马头山羊是湖北省本地品种，长期闭锁选育，湖北马头山羊具有个体大、生长快、繁殖率高等特点，其全年均可发情，母羊一年产2胎，年产羔率428%，平均窝产羔数2.14头，马头山羊繁殖力远远高于全国其它地方良种山羊品种（杨利国，2004）。波尔山羊是引入品种，平均窝产羔数为1.93头，低于马头山羊。这也与本研究结果符合，本研究结果显示AA基因型头胎和经产羔数均大于GG型，而马头山羊群体均为AA基因型。
3．2 MTNR1A不同基因型对波尔山羊产羔数的影响
本研究结果显示：在山羊样本中，无论是头胎还是经产山羊，褪黑激素受体1A不同的基因型对性状的影响表现都为AA＞GG，由于受到样本含量限制，该差异并未达到显著水平，但是该结果初步揭示褪黑激素受体1A基因可能与山羊繁殖性状相关，而该基因是否能够作为分子标记应用于山羊繁殖性状标记辅助选择，则需进一步扩大样本含量，增加品种个数进行验证。

参考文献
[1] 滑国华．山羊繁殖性状五个候选基因多态性及其与产羔数性状关联分析．［硕士学位论文］．武汉：华中农业大学图书馆，2006
[2] 吴伟生．抑制素基因作为山羊多胎性侯选基因的研究．［硕士学位论文］．武汉：华中农业大学图书馆，2006
[3] 周文然，储明星，孙少华等．小尾寒羊高繁殖力候选基因INHA的研究．农业生物技术学报,2007,15(1)：32-36
[4] 陈桂芳，李齐发，强巴央宗等．西藏牦牛和黑白花奶牛生长激素基因AluⅠ多态性的比较研究．畜牧与兽医，2002，35：11
[5] 储明星，程笃学，刘文忠．褪黑激素受体1A基因外显子2与小尾寒羊高繁殖力的关联．安徽农业大学学报， 2008，35 (1) ：124
[6] 王林枫．光照和埋植褪黑激素对内蒙古白绒山羊含氮物质分配和产绒性能的影响及调控的研究．［硕士学位论文］．内蒙古：内蒙古农业大学图书馆，2004
[7] 闫守庆 ,王慧芳 ,祝万菊．蓝狐与银狐褪黑激素受体 1A基因的克隆及序列分析．中国畜牧兽医， 2008，35：3
[8] 季从亮，储明星，陈国宏．4 个绵羊品种褪黑激素受体 1a 基因第二外显子PCR—RFLP 分析．华中农业大学学报，2003，4（22）：2
[9] 郭玉强 ,蔡国宝 ,刘艳婷等．褪黑激素对绵羊季节性繁殖的调控机理．甘肃畜牧兽医，2008，1：198
[10] 张凯，文勇立，李学伟等．褪黑激素受体基因表达的研究进展．四川生理科学杂志，2006;28（2）77-79
[11] 焦传珍．褪黑激素及其生理功能研究进展．生物学教学，2005，30（2）2-3
[12] 杨利国．湖北马头山羊．2003～2004 年全国养羊生产与学术研讨会议论文集，111-113
[13] [ 6 ] 季从亮，储明星，陈国宏等. 4个绵羊品种褪黑激素受体 1a基因第二外显子 PCR-RFLP分析 [J ]. 华中农业大学学报 ,2003,22（2）: 105-109
[14] Barrett P, Conway S, Jockers R, et al. Cloning and functional analysis of a polymorphic variant of the ovine Mel 1amelatonin receptor[J]. Biochimica etBiophysica Acta,1997,1356（3）: 299-307
[15] ChuM X, Ji C L, Chen G H. Association between PCR-RFLP of melatonin receptor 1a gene and high prolificacy in Small Tail Han sheep [J ]. Asian2Aust J Anim Sci,2003,16(12 ): 1701-1704MesserL A, LizhengWang, Tuggle C K, et al Mapping of the melatonin receptor 1a （M TNR1A） gene in pigs, sheep, and cattle[J ]. Mammalian Genome, 1997, 8（5）:368-370
[16] NotterD R,CockettN E,Hadfield T S. Evaluation of melatonin receptor 1a as a candidate gene influencing reproduction in an autumn-lambing sheep flock [J ]. J Anim Sci,2003,81 4 :912-917
[17] PelletierJ,BodinL,Hanocq E, et al. Association between expression of reproductive seasonality and alleles of the gene for Mel1a receptor in the ewe[J]. Biology of Repro-duction,2000,62 4 : 1096-1101
致 谢
本实验是在杨利国导师的悉心指导下顺利完成的。杨老师渊博的学识和严谨的治学态度潜移默化的影响着我，使我受益匪浅。在实验的过程中，滑国华师姐、于俊娜师姐和杨菲菲师姐在实验过程中教给我很多知识和技能，给了我很好的指导和很大的帮助，在此表示感谢。
同时还感谢实验室的各位师兄师姐给我耐心的指导和讲解，感谢与我一起做实验的毕业生同学的关心和鼓励！
18