《山羊褪黑激素受体1A基因HinfⅠ酶切多态性及其与产羔数性状的关联分析本科毕业论文.doc》由会员分享，可在线阅读，更多相关《山羊褪黑激素受体1A基因HinfⅠ酶切多态性及其与产羔数性状的关联分析本科毕业论文.doc（18页珍藏版）》请在三一办公上搜索。
1、本科毕业论文题 目：山羊褪黑激素受体1A基因Hinf酶切多态性及其与产羔数性状的关联分析姓 名：学 号专 业：动 物 科 学指导教师：职 称分类号 密级华中农业大学本科毕业论文 山羊褪黑激素受体1A基因Hinf酶切多态性及其与产羔数性状的关联分析Hinf-RFLP of Melatonin receptor 1 A(MTNR1A) gene and the association with the litter size traits in goat目 录摘 要4ABSTRACT5前 言61.材料与方法71.1 材料71.1.1 试验动物71.1.2 实验材料71.1.3主要试剂的来源71.1
2、.4主要仪器设备81.1.5主要试剂的配置81.1.6主要分子生物学软件101.2 方法101.2.1.山羊基因组DNA的提取101.2.2.引物设计及其序列101.2.3.PCR扩增及其电泳检测111.2.4. 限制性内切酶消化扩增产物以及电泳检测131.2.5.统计分析142 结果与分析152.1 PCR扩增结果152.2酶切结果152.3 褪黑激素基因频率和基因型频率分析结果162.4 MTNR1A基因型与波尔山羊头胎产羔数性状的关联分析172.5 MTNR1A基因型与经产产羔数性状的关联分析173 讨论1831 不同山羊品种基因型频率及基因频率差异1832 MTNR1A不同基因型对波尔
3、山羊头胎和经产羔数的影响18参考文献19致 谢20山羊褪黑激素受体1A基因HINF酶切多态性及其与产羔数性状的关联分析摘 要本研究以褪黑激素受体1A 基因为候选基因，分析该基因多态性对山羊产羔数性状的影响。采用限制性片段长度多态性（Restriction Fragment Length polymorphism， PCR-RFLP）方法，检测了65只波尔山羊和51马头山羊褪黑激素受体1A基因Hinf酶切多态性及其与产羔数性状的关联分析。结果表明，褪黑激素受体1A基因在波尔山羊中存在Hinf酶切多态性，而在马头山羊样本中则无多态性。在波尔山羊中，检测到AA、GG两种基因型，未检测到AG型；其中A
4、A基因型个体在出现频率最高，为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等位基因可能与波尔山羊产羔数呈正相关，该位点主要以加性方式起作用。关键词：山羊；繁殖性状；褪黑激素受体；PCRRFLPHinf-RFLP of Melatonin receptor 1 A（MTNR1A） gene and the association with
5、 the litter size traits in goatAbstractIn 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
6、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 th
7、at 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 w
8、ere 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前 言褪黑激素主要是由松果体合成的，其它合成部位还有:视网膜、眼眶腔的副泪腺、唾液腺、肠的嗜铬细胞及红细胞，松果体把对外界光照形成的神经信息经下丘脑视交叉上核转变成体液信息ME
9、L。1.褪黑激素生理功能、作用机制；褪黑激素对生殖系统作用的机理主要是通过光作用于视网膜，然后通过视网膜- 下视丘至视交叉上核，再经过一系列复杂的神经交换至颈上神经节，其节后交感神经作用于松果体。松果体分泌褪黑激素，通过下丘脑垂体性腺轴影响生殖系统。褪黑激素对生殖系统的作用因动物种类、生理状态不同而表现出抑制、促进或无作用的多重性（焦传珍，2005）。多项研究结果表明，褪黑激素虽然对牛、鼠类、禽等动物和人的生殖系统的发育有抑制作用，但是对绵羊、山羊、鹿等动物则明显表现出促进作用。2褪黑激素受体类型、作用机制；褪黑激素受体属于G蛋白耦联受体家族。根据其分子生物学方法又可分为MTNR1A、MTNR
10、1B、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多态
11、片段， R sa都存在295 bp 和290 bp 多态片段。Pelletier等（2000）用M n l限制性内切酶消化2组发情行为有明显差异（常年发情和季节性发情）的Merinos dArles母绵羊以及季节性发情明显的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个绵羊
12、品种中都存在303 bp和236 bp多态片段， R sa都以上迹象表明，这是直接从网上下载的。存在290 bp和267 bp 多态片段。在绵羊中的研究结果表明，该基因对绵羊繁殖性能存在影响。但是，对于山羊，只有对产绒或营养调控的研究，目前尚未发现有对山羊繁殖性能影响的相关报道。4本实验目的意义本实验室滑国华等（2006）在山羊MTNR1A中共发现7个SNPs，对其中4个SNPs分析结果显示：MTNR1A不同基因型对山羊产羔数性状有一定的影响。为进一步研究该基因对山羊产羔数性状的影响，本试验选取G493A突变位点，采用引入酶切位点方法，即限制性片段长度多态性Forced PCR-RFLP（Po
13、lymerase Chain Reaction-Restriction Fragment Length Polymorphsim），分析该位点在马头山羊和波尔山羊中的分布情况，并分析该位点对波尔山羊产羔数影响，为山羊标记辅助选择提供理论依据与滑国华的区别？更确切地说，创新点？或特色？。1.材料与方法1.1 材料1.1.1 试验动物 马头山羊51头，采自湖北省十堰市；波尔山羊65头，采自湖北省宜都市波尔山羊种羊场。颈静脉采血10ml/只，EDTA抗凝，-20保存。1.1.2 实验材料 1.5ml Eppendorf管、双面板、微量移液枪（1000l、200l、40l、10l、2.5l）及枪头、
14、一次性PE塑料手套。1.1.3主要试剂的来源TagDNA聚合酶购自上海生工生物工程技术有限公司；DNA限制性内切酶Hinf购自晶美生物工程有限公司；分子量标准100bp DNA Ladder和PBR322DNA/MspI购自上海生工生物工程技术有限公司；三羟甲基氨基甲烷（Tris）、苯酚、蛋白酶K、十二烷基硫酸钠（SDS）、甲酰胺、丙烯酰胺（Acrylamide，Acr）、NN亚甲基丙烯酰胺（Bisacrylamide）、琼脂糖、溴化乙锭（EB）、琼脂糖、乙二胺乙二酸二钠盐（EDTANa2）等均购自武汉凌飞科技有限公司。1.1.4主要仪器设备表1 主要仪器设备Table 1 Laborator
15、y 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美国LABC
16、ONCO公司移液器1ml；200l；40l； 10l；2.5l芬兰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 Na2EDTA2H2O溶于800ml双蒸水中，加入NaOH（约20g）调至pH=8.0，定容到1000ml。高压灭菌，室温保存；3) 10%SDS：称取50g SDS加热溶于400ml双蒸水中，定容至500ml。高压灭菌，室温保存；
17、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：称取NaAc3H2O溶于400ml双蒸水中，用HAc调至pH=5.2，加水定容至500ml。高压灭菌，室温保存；4) TE缓冲液：100mM Triscl（pH8.0），1mM EDTA（pH8.0）。1.1.5.3 用于琼脂糖凝胶电泳及检测溶液的配制1) 50TAE
18、储备液： 称取242g Tris碱溶于750ml双蒸水中，加入57.1ml HAc、100ml 0.5M EDTA（pH8.0），加水定容至1000ml；2) 5TBE储备液：称取54g Tris碱和27.5g硼酸溶于750ml双蒸水中，加入20ml 0.5M EDTA（pH8.0），加水定容至1000ml；3) 6上样缓冲液：0.25%溴酚兰，0.25%二甲苯青，15%Ficoll聚蔗糖；4) 溴化乙锭（EB，10mg/ml）:100ml双蒸水中加入0.1g EB，搅拌使完全溶解，转入棕色试剂瓶，锡箔包裹。1.1.5.4 用于聚丙烯酰胺凝胶电泳和银染的溶液的配制1) 30% 丙烯酰胺（丙烯酰
19、胺:甲叉双丙烯酰胺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：量取
20、30ml冰HAc溶于970ml双蒸水，混匀。1.2 方法1.2.1.山羊基因组DNA的提取从山羊颈静脉采血10ml，放入装有抗凝剂（0.5mol/L EDTANa2）的离心管中，6000r/min离心15min，弃上清，吸取白细胞沉淀（参杂少量红细胞）于另一离心管中；加入两倍体积双蒸馏水（灭菌），摇匀沉淀，静置10min15min，破碎红细胞；6000r/min离心15min，轻轻弃去上清，留白细胞沉淀；加入1SET 5ml，蛋白酶K（10mg/ml）0.1ml至终浓度200ng/ml，10% SDS 250l至终浓度0.5%，55水浴消化过夜；加入等体积饱和平衡酚轻轻摇匀，12000r/mi
21、n离心15min，轻轻吸取上层水相于另一离心管中，弃下相苯酚；重复上步；加入等体积苯酚/氯仿/异戊醇（25:24:1），轻轻摇动10min，12000r/min离心15min，轻轻吸取上层水相于另一离心管中，弃下层有机相；加入等体积氯仿/异戊醇（24:1），轻轻摇动10min，12000r/min离心15min，轻轻吸取上层水相于另一离心管中，弃下层有机相；加入两倍体积预冷无水乙醇，轻轻转动管子，DNA呈絮状析出，用1ml吸头（灭菌）吸出DNA沉淀，放入1.5ml离心管中，加入预冷无水乙醇洗涤一次，离心弃去无水乙醇，再用75%乙醇漂洗2次，小心吸去乙醇，待DNA沉淀自然干燥后，加入适量TE（1
22、00l300l）溶解。1.2.2.引物设计及其序列引物参考滑国华等（2008）进行设计，序列如下：上游引物: 5ACGACCCGAGGATCTATTCC 3 下游引物: 5ATATAAGTCCTGGGGCTTCAGATT 3，在下游引物中引入错配碱基A，当突变位点为G（即反义链为C）时，可形成Hinf（GATTC）酶切位点；当突变位点为A（即反义链为T）时，该酶切位点缺失，由此可进行基因型判定。1.2.3.PCR扩增及其电泳检测1.2.3.1 PCR扩增反应体系（表2）表2 PCR扩增反应体系（20l）Table 2 PCR amplification of the response syst
23、em （20l）反应体系response system反应体积V（l） response volume（l）灭菌蒸馏水14.610Buffer2Mg2+1.2dNTP0.2上游引物0.4下游引物0.4Taq DNA聚合酶（5 U/ml）0.2模板DNA1Total201.2.3.2 退火温度的优化利用PCR机器上都有一个Gradient功能，在一次PCR过程中，对机器中不同位点的反应采用不同的退火温度，一般在第一次实验中，采用引物溶点上下10度设置退火温度梯度。每个温度梯度一个反应，最后电泳检测并比较各个退火温度的条带效果，得到最好结果。本试验中退火温度设定为5565，结果表明，退火温度60时
24、，无非特异性带，目的条带最为清晰。1.2.3.3 PCR反应程序（表3）表3 PCR反应程序Table 3 PCR amplification of the response term步骤stepsPCR反应条件PCR reaction conditions温度temperature时间time1预变性955min2C变性9530sec3C退火6030sec4C延伸7230sec5延伸727min6保温1610min第二步至第四步循环35次。1.2.3.4 琼脂糖凝胶电泳检测 PCR扩增产物用1.5%琼脂糖凝胶电泳检测， EB染色，利用BIORAD凝胶成相系统凝胶成像系统进行成像分析，记录结果
25、并拍照。1.2.3.4.1 琼脂糖凝胶电泳方法（1） 灌胶将透明胶带粘于灌胶模槽的两端，制成灌胶模。放入梳子。将模放于水平台上。取一烧杯，放入：琼脂糖0.75 g； 1TAE 50 ml。混匀，加热，至琼脂糖彻底溶解。待稍冷后，加入5ul饱和的溴乙锭，轻晃混匀，然后轻轻倒入模子。用吸头吸去表面小泡。待冷却成胶。（2）上样电泳 轻轻拔出梳子，将胶两头的透明胶带揭去，移入电泳槽，点样孔一侧靠近操作者，注入1TAE缓冲液，浸没凝胶。使用微量可调采样器取3l DNA充分混匀于适量溴酚兰，点样于1.5琼脂糖凝胶孔内，并点3l 600bp Marker作为参照。点样时，加样头尖端插入上样孔内1/3高度，轻
26、轻将样品推入，避免刺入前后凝胶和刺穿凝胶底部。使比重较大的样品自然沉入上样孔底。接上电源：阳极远离上样孔，阴极靠近上样孔（DNA向阳极移动）。开启电源，电压调至100V。2030min后检查。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.710Buf
27、fer1Hinf限制性内切酶0.3PCR产物4Total10将上述酶切体系置37水浴锅或恒温箱中酶切过夜。1.2.4.2非变性聚丙烯酰胺凝胶电泳检测1.2.4.2.1聚丙烯酰胺凝胶电泳操作步骤不同的实验需要优化聚丙烯酞胺凝胶的交联度和浓度，才能达到最好的实验效果，本实验采取一下方法：1. 在小烧杯中加入18.43m1双蒸水，9.33ml 30%丙烯酞胺（29:1），7ml 的5TBE混匀；2. 清洗制胶用玻璃板，烘干后用透明胶封好底边和侧边，固定在制胶用的有机玻璃板上；3. 向小烧杯中加入245ul的10%AP，13ul TEMED混匀后迅速灌胶；4. 当胶灌至离玻璃板上边缘0.5cm时停止灌
28、胶，插入梳子，注意梳齿下不要有气泡，室温聚合1小时；5. 凝胶聚合好后，拔出梳子，注射器吸出梳齿孔中未聚合的液体；6. 取下玻璃板，撕下胶线，固定在电泳槽上，倒入电泳缓冲液1TBE；7. 在PCR产物中加入1/6体积上样缓冲液，用微量进样器取5u1酶切样点样，同时点分子量标记物PBR322/Mspl； 8. 上样完毕后，打开电泳议以80V电泳20min，再调到120V电泳45h。1.2.4.2.2聚丙烯酰胺凝胶染色步骤1. 拆下电泳装置，剥下聚丙烯酰胺凝胶，放到一个塑料盘内，用双蒸水冲漂洗2次； 2. 加入500ml 10乙醇固定液，在摇床上固定810 min，弃去反应液； 3. 加入500m
29、l 1HNO3氧化，在摇床上轻摇810 min，弃去反应液，用双蒸水漂洗2次；4. 倒入500ml 0.2AgNO3银染液，染色1530 min，弃去反应液，用双蒸水漂洗2次； 5. 加入500ml 3Na2CO3显色液，手摇显色至出现清晰的银染带，迅速弃去反应液；6. 用500ml 3HAC停止显色；7. 将已经染好的聚丙烯酰胺凝胶置于透射白光板上，观察酶切反应结果、记录并拍照。1.2.5.统计分析基因频率与基因型频率使用软件SPSS10.0进行统计分析，不同基因型对产羔数的最小二乘均值分析使用SAS（8.0版）GLM（General Linear Model）软件包进行，模型如下：Yijk
30、l=+Gi+Pi+eijkYijkl为产羔数，为群体均值，Gi为基因型效应，Pi为胎次效应，eijk为随机效应。基因效应分析：加性效应=（GG-AA）/2显性效应=AG-（AA+GG）/22 结果与分析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
31、 PCR product.M，1-3，4-6: 100bp DNA Ladder， PCR product of ?goat， PCR product of ?goat2.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 ；而马头山羊
32、在118处的酶切位点则无多态性，基因型均为AA型。酶切图谱见图3 。图2 波尔山羊酶切图谱注：M， pBR322/Msp 1 Marher；113为酶切产物条带Fig.2 genotypes of the M TNR1A fragment defined by PCR-RFLP图3 马头山羊酶切图谱注：M， pBR322/Msp 1 Marher18为酶切产物条带Fig.3 genotypes of the M TNR1A fragment defined by PCR-RFLP2.3 褪黑激素基因频率和基因型频率分析结果分析两种山羊品种的褪黑激素基因型频率发现，在该酶切位点，马头山羊只出现A
33、A基因型，波尔山羊虽然出现两种基因型，但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波尔山羊65AA0.554A 0.554GG0.446G 0.446马头山羊51AA1A 12.4 MTNR1A基因型与波尔山羊头胎产羔数性
34、状的关联分析由表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
35、 meansstandard error加性效应additive effec显性效应dominant effects波尔山羊65AA1.45450.52220.035-1.42GG1.38460.50642.5 MTNR1A基因型与经产产羔数性状的关联分析由表8可知，波尔山羊的AA型经产羔数比GG型多0.09只；加性效应值为0.046，显性效应值为-1.84。结果表明：A等位基因可能与波尔山羊产羔数呈正相关，该位点主要以加性方式起作用。表8 波尔山羊MTNR1A基因型对经产产羔数性状的最小二乘均值及基因效应分析Table 8 Least squares means and allele sub
36、stitution effects for litter size in multiparous boar goat基因型Genetype最小二乘均值标准误Least squares meansstandard error加性效应additive effects显性效应dominant effectsAA1.88270.58710.046-1.84GG1.79050.59883 讨论31 不同山羊品种基因型频率及基因频率差异利用荧光原位杂交方法和遗传连锁分析方法，褪黑激素受体1A已经在绵羊等多种动物中得到精细定位。并有多个研究结果表明，褪黑激素受体1A基因与绵羊的繁殖性状有着紧密的关联。在本研
37、究所检测的两个山羊群体中，基因型频率及基因频率存在较大差异。本研究结果显示：马头山羊是AA基因型纯合群体，而波尔山羊中存在AA和GG两种基因型。这可能与两个品种的育成史有关。马头山羊是湖北省本地品种，长期闭锁选育，湖北马头山羊具有个体大、生长快、繁殖率高等特点，其全年均可发情，母羊一年产2胎，年产羔率428%，平均窝产羔数2.14头，马头山羊繁殖力远远高于全国其它地方良种山羊品种（杨利国，2004）。波尔山羊是引入品种，平均窝产羔数为1.93头，低于马头山羊。这也与本研究结果符合，本研究结果显示AA基因型头胎和经产羔数均大于GG型，而马头山羊群体均为AA基因型。32 MTNR1A不同基因型对波
38、尔山羊产羔数的影响本研究结果显示：在山羊样本中，无论是头胎还是经产山羊，褪黑激素受体1A不同的基因型对性状的影响表现都为AAGG，由于受到样本含量限制，该差异并未达到显著水平，但是该结果初步揭示褪黑激素受体1A基因可能与山羊繁殖性状相关，而该基因是否能够作为分子标记应用于山羊繁殖性状标记辅助选择，则需进一步扩大样本含量，增加品种个数进行验证。参考文献1 滑国华山羊繁殖性状五个候选基因多态性及其与产羔数性状关联分析硕士学位论文武汉：华中农业大学图书馆，20062 吴伟生抑制素基因作为山羊多胎性侯选基因的研究硕士学位论文武汉：华中农业大学图书馆，20063 周文然，储明星，孙少华等小尾寒羊高繁殖力
39、候选基因INHA的研究农业生物技术学报,2007,15(1)：32-364 陈桂芳，李齐发，强巴央宗等西藏牦牛和黑白花奶牛生长激素基因Alu多态性的比较研究畜牧与兽医，2002，35：115 储明星，程笃学，刘文忠褪黑激素受体1A基因外显子2与小尾寒羊高繁殖力的关联安徽农业大学学报， 2008，35 (1) ：1246 王林枫光照和埋植褪黑激素对内蒙古白绒山羊含氮物质分配和产绒性能的影响及调控的研究硕士学位论文内蒙古：内蒙古农业大学图书馆，20047 闫守庆 ,王慧芳 ,祝万菊蓝狐与银狐褪黑激素受体 1A基因的克隆及序列分析中国畜牧兽医， 2008，35：3 8 季从亮，储明星，陈国宏4 个绵
40、羊品种褪黑激素受体 1a 基因第二外显子PCRRFLP 分析华中农业大学学报，2003，4（22）：29 郭玉强 ,蔡国宝 ,刘艳婷等褪黑激素对绵羊季节性繁殖的调控机理甘肃畜牧兽医，2008，1：19810 张凯，文勇立，李学伟等褪黑激素受体基因表达的研究进展四川生理科学杂志，2006;28（2）77-7911 焦传珍褪黑激素及其生理功能研究进展生物学教学，2005，30（2）2-312 杨利国湖北马头山羊20032004 年全国养羊生产与学术研讨会议论文集，111-11313 6 季从亮，储明星，陈国宏等. 4个绵羊品种褪黑激素受体 1a基因第二外显子 PCR-RFLP分析 J . 华中农业
41、大学学报 ,2003,22（2）: 105-10914 Barrett P, Conway S, Jockers R, et al. Cloning and functional analysis of a polymorphic variant of the ovine Mel 1amelatonin receptorJ. Biochimica etBiophysica Acta,1997,1356（3）: 299-30715 ChuM X, Ji C L, Chen G H. Association between PCR-RFLP of melatonin receptor 1a gen
42、e 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 cattleJ . Mammalian Genome, 1997, 8（5）:368-37016 NotterD R,CockettN E,Hadfield T S. Evalua
43、tion of melatonin receptor 1a as a candidate gene influencing reproduction in an autumn-lambing sheep flock J . J Anim Sci,2003,81 4 :912-91717 PelletierJ,BodinL,Hanocq E, et al. Association between expression of reproductive seasonality and alleles of the gene for Mel1a receptor in the eweJ. Biology of Repro-duction,2000,62 4 : 1096-1101致 谢本实验是在杨利国导师的悉心指导下顺利完成的。杨老师渊博的学识和严谨的治学态度潜移默化的影响着我，使我受益匪浅。在实验的过程中，滑国华师姐、于俊娜师姐和杨菲菲师姐在实验过程中教给我很多知识和技能，给了我很好的指导和很大的帮助，在此表示感谢。同时还感谢实验室的各位师兄师姐给我耐心的指导和讲解，感谢与我一起做实验的毕业生同学的关心和鼓励！
