奶山羊PITX1基因单核苷酸多态性位点及其检测方法和应用技术领域
本发明属于动物分子遗传育种技术领域,涉及基因单核苷酸多态性(SNP)的检测及其
应用,具体涉及奶山羊PITX1(成对同源异型结构域蛋白转录因子1)基因(以序列
Accession No.NW_003104033为参考)单核苷酸多态性位点及其检测方法和应用。
背景技术
与奶牛相比,奶山羊更适用于在恶劣环境条件下饲养,因为它具有较强抗逆性;与
牛奶相比,羊奶具有更高的营养价值,因为它的脂肪球小,不含过敏原,酸值低,比牛
奶更易为人体吸收,在国际营养界被誉为“奶中之王”。然而目前,我国奶山羊产业相对
滞后,利用优良奶山羊群体获得羊奶的数量与质量,与发达国家相比差距巨大,与目前
国内主流的奶牛业发展“不可同日而语”。为此,作为我国“十二五”规划中奶业发展计划
的重要组成部分,奶山羊业的发展迫在眉睫,但是目前,我国该产业发展滞后,其根本
原因在于我国奶山羊品种种源不足和种群遗传品质较差。因此,如何提高奶山羊种群遗
传品质,从而增加奶山羊的产奶量并改善乳品质,是目前我国奶山羊遗传育种专家尚待
解决的关键问题。
奶山羊泌乳性能主要包括产奶量、乳脂率、乳脂量、乳蛋白率、乳蛋白量、乳固体
物含量等,属于具有重要经济价值的数量性状位点(QTL),由微效多基因控制,且遗传力
较低。在奶山羊育种方面,仅靠传统育种手段对这些泌乳性状进行改良,无疑是非常缓
慢的,而且效果也不会很理想,因此,应用DNA标记的现代育种技术是加快良种选育和
提高种群遗传品质,从而增加奶山羊的产奶量和改善乳品质的先进的有效的方法。因为
DNA标记具有多态性丰富、遗传稳定、不受组织、生理发育阶段及环境影响等诸多优点。
此外,随着DNA标记的种类与数目的增多,基因图谱密度与精度在不断提高,为遗传改
良提供了新的策略,使DNA标记辅助选择(MAS)成为一种具有广阔前景的分子育种技术。
由此可见,利用DNA标记介导的MAS方法和传统育种方法对奶山羊泌乳性状进行改良、
提高,必将极大地加速选择进展,提高选择准确性。
目前,利用DNA标记进行MAS育种的主要步骤如下:首先,在DNA水平上快速
筛查与检测影响奶山羊泌乳性状的密切相关的基因;其次,建立基因多态与泌乳性状的
关联分析,确定与泌乳性状显著关联的DNA标记或SNP位点;最后,实现遗传标记辅
助选择和实现早期诊断选择。具体而言,关于DNA标记与奶山羊泌乳性状之间关系的研
究,主要有两种基本方法:候选基因分析和连锁标记分析。在奶山羊研究中,候选基因
分析比连锁标记分析更易操作,且不需要进行特意实验设计、费用低、统计功效强,被
广泛采用,适用于任何可取得基因型和表型数据的群体。该方法通过统计方法建立相应
的数学模型,从而确定该候选基因与未知基因的关系或效应。通常来说,候选基因是一
些与泌乳性状有关的重要基因,可能是结构基因、调节基因或是在生化代谢途径中影响
该性状表达的基因。
众所周知,生长激素(GH)能够促进哺乳动物生长、泌乳、细胞增殖。泌乳素(PRL)
也能够显著地影响哺乳动物的生长、泌乳以及免疫系统;而且,已有研究表明,GH和
PRL基因突变能够显著影响哺乳动物的泌乳期和泌乳性能。因此,GH基因、PRL基因以
及它们的正调控因子基因均可作为影响泌乳性能的关键候选基因。同源异型结构域蛋白
转录因子(PITX)类转录因子成员在基因活化、垂体发育中的的关键性作用,启发此类
基因与泌乳性能的关联性研究。作为PITX家族的重要成员,同源异型结构域蛋白转录因
子1(PITX1)基因被分别定位在人类第5号染色体、小鼠和原鸡第13号染色体、牛第7
号染色体,但山羊的相应信息尚未公布。序列比对结果表明,PITX1在人类、黑猩猩、
狗、鼠、鸡、斑马鱼和牛中的序列高度保守,在许多发育通路中发挥重要作用,包括影
响脑垂体的器官发生,影响GH、PRL、促甲状腺激素(TSH)的表达从而影响动物的泌乳
性能。值得一提的是,PITX1是在发育中的小鼠、鸡的胚胎的后肢、口腔上皮组织、臂
弓等处被发现的,其决定着肌肉、骨骼等的形态。此外,敲除PITX1基因的小鼠表现出
在垂体、后肢和味觉发育上的缺陷和障碍。因此,PITX1基因是垂体-下丘脑轴相关激素
基因(如:PRL和GH)的重要的调节器,是已知的代表垂体发育的重要标记基因之一。
所谓单核苷酸多态性(SNP)就是指基因组DNA序列中由于单个核苷酸(A/T/C/G)的变
化而引起的多态性,其变异形式主要有颠换、转换、插入和缺失等。近年来,已发展出
了许多用于探寻基因SNP的方法,例如:DNA序列测定法、PCR-SSCP与DNA测序结
合法、PCR-RFLP法、AS-PCR方法、引物延伸法和寡核苷酸连接反应等。在这些检测技
术中,DNA序列测定法是最为准确检测SNP的方法,但是,其检测费用极其昂贵,且需
要有DNA测序仪等大型仪器,同时,在测序过程中需要非常熟练的有经验的技术人员,
所以,DNA序列测定法不是一种应用于生产实践的理想的SNP检测方法;利用PCR-SSCP
与DNA测序结合法检测SNP,可以适当降低检测费用,但是,其实验过程比较长,操作
比较繁琐,且实验过程中存在假阴性问题,所以,也并非理想的SNP检测手段;AS-PCR
方法作为一种新型的SNP检测方法,在未来的应用领域中具有非常广阔的前景,但是该
方法需要设计特别的引物,且只能针对特定的基因位点,同时检测过程中还存在误检的
概率,因此,目前不具有普遍应用的特点;而用引物延伸法和寡核苷酸连接反应技术检
测SNP位点,需要平板读数仪、基因芯片、微球阵列技术和质谱仪等检测平台,对于一
般的分子实验室来说可实施性不强。综上所述,PCR-RFLP法是当前检测SNP最为理想
的遗传标记方法。
最新研究表明,关于PIX1基因SNP研究多见于人、小鼠等生物的体外表达分析,
且常被作为器官机能缺陷和重要的疾病预测相关的候选基因,而鲜见于家畜和家禽SNP
分析。截至目前,关于奶山羊PITX1基因SNP研究尚未见相关报道。中国奶山羊PITX1
基因SNP领域的研究匮乏,基因位点的功能研究及其遗传变异与泌乳性能(如:产奶量、
乳脂率、乳脂量、乳蛋白率、乳蛋白量、乳固体物含量等)的关联研究仍是空白。由于PITX1
基因SNP位点在动物生产实践中对形成胚胎、调控生长激素和促进垂体分泌激素等具有
重要的作用,本发明提供的检测方法为PITX1基因SNP与奶山羊泌乳性状之间关系的建
立奠定了基础,以便用于中国奶山羊泌乳性状的标记辅助选择(MAS),快速建立遗传资
源优良的奶山羊种群,从而提高种群遗传品质,增加奶山羊的产奶量和改善乳品质,为
我国奶山羊业的发展提供理论与实践资料。
发明内容
本发明解决的问题是,利用PCR-RFLP方法快速检测奶山羊PITX1基因SNP,并将
其与奶山羊泌乳性状进行关联分析,验证其作为奶山羊分子育种中辅助选择的DNA标
记,从而加快良种选育速度。
为了实现上述发明的目的,本发明采取如下的技术方案:
奶山羊PITX1基因单核苷酸多态性位点,其中,该基因单核苷酸多态性位点包括奶
山羊PITX1基因第201位为G或A的单核苷酸多态性位点,基因型为GG、GA、AA。
一种检测奶山羊PITX1基因单核苷酸多态性位点的方法,包括下述步骤:
(1)、以包含PITX1基因的DNA序列为模板,以引物对P为引物,PCR扩增奶山羊
PITX1基因;
(2)、然后利用限制性内切酶MspI对步骤(1)获得的PCR产物进行酶切;
(3)、再通过琼脂糖凝胶电泳检测步骤(2)获得的酶切产物即可鉴定奶山羊PITX1基因
的单核苷酸多态性位点包括第201位为G或A的单核苷酸多态性位点,基因型为GG、
GA、AA;
所述的引物对P为:
上游引物:F:5′-CCGCCTTCCACCTAGCCCGCAC-3′;
下游引物:R:5′-TCGTCCATGTCCACGTTCATCG-3′。
上述技术方案中所述的一种检测奶山羊PITX1基因单核苷酸多态性位点的方法,其
中,所述的PCR扩增反应程序为:
94℃预变性5min;35个循环(94℃变性30s,65.2℃退火30s,72℃延伸30s);
72℃延伸10min。
上述技术方案中所述的一种检测奶山羊PITX1基因单核苷酸多态性位点的方法,其
中,所述的琼脂糖凝胶电泳为2.0%琼脂糖凝胶电泳。
上述技术方案中所述的奶山羊PITX1基因单核苷酸多态性位点在作为奶山羊泌乳性
状的标记辅助选择育种中提高乳脂量和乳固体物含量的DNA遗传标记中的应用。
上述技术方案中所述的应用,其中,奶山羊PITX1基因单核苷酸多态性位点的GA
基因型作为提高奶山羊泌乳性能的DNA遗传标记的应用。
在本发明中检测奶山羊PITX1基因单核苷酸多态性位点的方法中,以奶山羊基因组
DNA为模板,在Taq DNA聚合酶、Buffer(缓冲环境)、Mg2+、dNTPs存在的情况下,利
用引物对P进行PCR扩增,然后利用限制性内切酶对扩增得到的DNA片段进行酶切,
再通过琼脂糖凝胶电泳检测即可鉴定奶山羊PITX1基因第201位即内含子1第41位
(NW_003104033:g.201G>A(IVS1+41G>A)的单核苷酸多态性;本发明的buffer中不含有
Mg2+。
本发明中奶山羊PITX1基因第201位为G或A的单核苷酸多态性位点中,A等位基
因表现为230bp,G等位基因表现为:121bp和109bp,相应地,GG基因型表现:121
bp,109bp;GA基因型表现:230bp,121bp,109bp;AA基因型表现:230bp。
本发明根据牛PITX1基因的序列(GenBank Accession number:NW_003104033)设计引
物对P(由于NCBI未公布山羊基因组序列),以关中奶山羊基因组DNA池为模板,进行
PCR扩增,并对PCR产物进行测序,测序后得到奶山羊PITX1基因的部分序列。与NCBI
公布的序列进行比较后,发现在第201位(内含子1的第41位)存在SNP。
针对上述第201位的SNP,本发明还公开了其筛查和检测方法,通过设计特定的引
物,PCR扩增目的片段,然后用特定的限制性内切酶酶切鉴定,能够简单、快速、低成
本、精确地检测其单核苷酸多态性。
本发明具有以下有益效果:
1、本发明公开了奶山羊PITX1基因单核苷酸多态性位点的筛查和检测方法,通过设
计特定的PCR引物扩增片段,能够用RFLP方法简单、快速、成本低、精确的检测其单
核苷酸的多态性;
2、本发明对关中奶山羊该SNP基因型进行了检测和基因频率分析,对上述SNP位
点与奶山羊的泌乳性能(包括羊奶的乳脂含量和乳固体物含量)进行关联分析,结果表明该
位点可以作为提高奶山羊泌乳性能的分子标记;
3、本发明公开的检测方法为奶山羊PITX1基因的SNP与奶山羊的泌乳性状关系的
建立奠定了基础,以便用于奶山羊泌乳性状的标记辅助选择(MAS)育种中用于快速建
立遗传资源优良的奶山羊种群。
附图说明
1、图1为奶山羊PITX1基因第201位GA基因型个体的测序图;
2、图2为MspI PCR-RFLP方法检测山羊PITX1基因第1内含子
SNP(NW_003104033:g.201G>A or IVS1+41G>A)序列分析图;
3、图3为MspI PCR-RFLP方法检测奶山羊PITX1基因第201位(内含子1第41
位)SNP酶切电泳分型图。
具体实施方式:
为使本发明的技术方案便于理解,以下结合具体实施例对本发明作进一步的说明。
本发明下述实施例中所用的主要试剂及其来源如下:
1、生化试剂与生物学试剂:
①Taq DNA聚合酶(购自MBI公司);②限制性内切酶MspI(购自MBI公司);
③dNTPs(购自MBI公司);④Marker I:(购自北京天根生化科技有限公司);⑤蛋白酶
K(购自华美生物工程公司)。
2、普通试剂:
普通试剂从华美生物工程公司购买,为进口分装产品:柠檬酸、柠檬酸钠、葡萄糖、
Tris、EDTA、NaCl、NaOH、KCl、Na2HPO4、KH2PO4、Tris饱和酚、氯仿、异戊醇、无
水乙醇、醋酸钠、十二烷基磺酸钠(SDS)、溴化乙锭(EB)、溴酚蓝、二甲基苯氰FF、乙
酸、蔗糖、硼酸、琼脂糖等。
3、溶液与缓冲液
所有溶液与缓冲液均采用去离子超纯水配制。高压灭菌条件为15bf/in(1.034×105Pa),
25min。试剂配制方法均参考Sambrook等编著的《分子克隆实验指南》。
(1)、样品采样所用溶液:
抗凝剂ACD:柠檬酸0.48g;柠檬酸钠1.32g;葡萄糖1.47g。把它们定容于100mL
水中,高压灭菌。每6mL新鲜血液中加入1mL的ACD液。这一抗凝剂优于EDTA,在
血液贮存过程中能更好地保存高分子的DNA,经其抗凝的血液可以在0℃保存数天或-70
℃长期保存。
(2)、血样基因组DNA分离所用溶液:
①、PBS缓冲液:NaCl 8g,KCl 0.2g,Na2HPO4 1.44g,KH2PO4 0.24g,加超纯水
至1000mL,调pH至7.4,高压灭菌;
②、10%SDS:10g SDS溶解于90mL的超纯水中,68℃水浴溶解,用HCl调pH
至7.2,定容至100mL;
③、0.5mol/L EDTA:EDTA 186.1g,溶于800mL的超纯水中,用NaOH调pH至
8.0,定容至1000mL,高压灭菌,4℃保存;
④、1mol/LTris·Cl:121.14g Tris,溶于800mL超纯水中,HCl调节pH至8.0,定
容至1000mL,高压灭菌,4℃保存;
⑤、5mol/LNaCl:NaCl 292.2g溶于1000mL超纯水中;
⑥、DNA抽提缓冲液:取0.5mmol/L EDTA 40mL,1mmol/L Tris·Cl 10mL;
⑦、5mmol/L NaCl 4mL,10%SDS 10mL定容至100mL。实际浓度为200mmol/L
EDTA,pH 8.0;100mmol/L Tris·HCl,pH 8.0;200mmol/L NaCl,2%SDS;RNase 20μg/mL;
⑧、NaAc缓冲液:NaAc·3H2O 20.4g;超纯水40mL;稀HAc调pH至7.4;定容至
50mL;
⑨、TE缓冲液:Tris·Cl缓冲液(pH 8.0)10mmol/L,EDTA缓冲液(pH 8.0)0.1mmol/L,
高压灭菌,4℃保存;
⑩、蛋白酶K:用超纯水配成20mg/mL,-20℃保存;
(3)、组织样DNA提取所用溶液:
除了基因组DNA提取时的公用溶液外,还配置以下试剂:
①、2mol/L NaCl:11.688g溶于水,定容至100mL,高压灭菌;
②、组织DNA提取液(100mL):1mol/L Tris-Cl(pH 8.0)l mL,0.5mol/L EDTA(pH
8.0)20mL,2mol/L NaCl 5mL,定容至100mL;
(4)、琼脂糖电泳分析所用溶液
①、0.5×TBE缓冲液:取10×TBE 50mL定容至1000mL;
②、上样缓冲液:0.25%溴酚蓝,0.25%二甲苯青FF,40.0%(w/v)蔗糖水溶液。
实施例一:奶山羊PITX1基因部分DNA序列的克降与DNA测序
一、奶山羊样品的采集
提取DNA所用血样(静脉采血)和耳组织样均采自658只健康且无血缘关系的2-4
岁关中奶山羊,采自陕西省三原县奶山羊育种场。
二、血样基因组DNA的提取与分离
(1)、冷冻血样(主要为血细胞)室温解冻,转移500μL至1.5mL Eppendorf管,加入等
体积PBS液,充分混匀,12000r/min离心10min(4℃),弃去上清液,重复上述步骤至上
清液透明、沉淀呈淡黄色。
(2)、在离心管中加入DNA抽提缓冲液500μL,摇动,使血细胞沉淀脱离离心管管
壁,37℃水浴1h。
(3)、加蛋白酶K至30μL(20mg/mL)并混匀,55℃过夜至澄清,尚未澄清者,可补
加10μL蛋白酶K混匀继续消化直至澄清。
(4)、将反应液冷却至室温,加Tris-饱和酚500μL,温和摇动离心管20min,使其充
分混匀;4℃,12000r/min离心10min,将上清液转入另一1.5mL离心管中,重复一次。
(5)、加氯仿500mL,充分混匀20min,4℃,12000r/min离心10min,
将上清液转入另一1.5mL离心管中。
(6)、加氯仿、异戊醇混合液(24∶1)500mL,充分混匀20min,4℃,12000r/min离
心10min,将上清液转入另一1.5mL离心管中。
(7)、加0.1倍体积的NaAc缓冲液及2倍体积的冰冷的无水乙醇,混合转动离心管直
至白色的絮状沉淀析出,-20℃保存30~60min。
(8)、4℃,12000r/min离心10min,弃去上清液,用70%冰冷乙醇漂洗DNA沉淀2
次。
(9)、4℃,12000r/min离心10min,弃去上清液,室温下使乙醇挥发干净。
(10)、干燥后的DNA溶于80~100μL的TE液,4℃保存直至DNA完全溶解,0.8%
琼脂糖凝胶电泳检测其质量,-80℃保存。
三、组织样品DNA的提取与分离
(1)、取奶山羊耳组织样约10mg的组织样品,放于1.5mL的离心管中,用小剪刀尽
量剪碎。
(2)、加入600μL组织提取液,10%SDS至终浓度为1%,蛋白酶K至终浓度为100
μg/mL,55.0℃消化过夜,最好保证组织样较均匀地分布在组织提取液中。
(3)、将溶液冷却至室温,加入等体积的Tris饱和酚,盖紧管盖,缓慢地来回颠倒离
心管,至少持续10min以上,12000r/min离心15min。
(4)、取上清液,加入等体积的酚∶氯仿(1∶1),盖紧管盖,缓慢地来回颠倒离心管,
至少持续10min以上,12000r/min离心15min。
(5)、取上清液,加入等体积的氯仿∶异戊醇(24∶1),盖紧管盖,缓慢地来回颠倒离
心管,至少持续10min以上,12000r/min离心15min。
(6)、取上清液,加入2倍体积的冰冷无水乙醇和1/10体积的3mol/L乙酸钠,盖紧
管盖,缓慢地来回颠倒离心管,直至出现白色絮DNA。
(7)、挑出DNA,放进一个1.5mL的离心管中,加入500μL 70%乙醇,
盖紧管盖,缓慢地来回颠倒离心管,然后12000r/min离心3~5min,小心倒掉乙醇,将
管倒置于吸水纸上。
(8)、再一次向离心管中加入500μL 70%乙醇,盖紧管盖,缓慢颠倒,然后12000r/min
离心3~5min,小心倒掉乙醇,将管倒置于吸水纸上。
(9)、待干燥后,加入60μL灭菌超纯水,4℃保存过夜,待检测。
四、DNA的纯化
(1)、500μL的DNA溶液中加入10%SDS使其终浓度为0.1%,加入蛋白酶K至终
浓度达到100μg/mL。
(2)、5℃保温10h左右。
(3)、等体积苯酚∶氯仿∶异戊醇(25∶24∶1)和氯仿分别抽提一次。
(4)、12000r/min离心5min分相,吸取上层水相至另一离心管中。
(5)、加入1/10体积3mol/L醋酸钠和2倍体积冰冷无水乙醇沉淀DNA。
(6)、倒掉液体,70%乙醇洗涤后凉干,加入60μL灭菌超纯水溶解,4℃待检测。
五、分光光度法检测DNA
用紫外光光度计测定DNA样品在260nm、280nm处的OD值。计算DNA含量和
OD260/OD280的比值。如OD260/OD280比值小于1.6,说明样品中含有较多的蛋白质或酚,
则应进行纯化;若比值大于1.8,则应该考虑去除RNA纯化。
DNA浓度(ng)=50×OD260值×稀释倍数
DNA检测完毕后,取适量稀释至50ng/μL,存于-20℃备用,其余的存放于-80℃。
六、DNA池的构建
DNA检测完毕后,从100个浓度为50ng/μL的关中奶山羊DNA样品中各取10μL
混合,构建关中奶山羊群体DNA池。
七、引物的设计
由于目前美国国立医学图书馆(NCBI)网站上并未公布奶山羊PITX1基因的序列,为
此,以NCBI上公布的牛PITX1基因(GenBank accession No.NW_003104033)序列为参考,
利用软件Primer 5.0设计了扩增奶山羊PITX1基因第1内含子的230bp片段的PCR引物
对P,其引物对P的序列信息如下:
上游引物:F:5′-CCGCCTTCCACCTAGCCCGCAC-3′(nt50-71);
下游引物:R:5′-TCGTCCATGTCCACGTTCATCG-3′(nt 259-279)
八、PCR克隆奶山羊PITX1基因
以构建好的关中奶山羊群体DNA池为模板,用上述设计的引物进行PCR扩增,PCR
反应体系见表1。
表1 PCR反应体系
灭菌超纯水(H2O)
15.75μL
10×PCR缓冲液(MBI,Fermentas)
2.5μL
MgCl2(25mmol/L)
1.5μL
dNTPs(2.5mmol/L)
2.5μL
上游引物(10pmol/L)
0.5μL
下游引物(10pmol/L)
0.5μL
Taq DNA聚合酶(0.5U/μL)
1.25μL
DNA模板(100ng/μL)
0.5μL
总体积
25μL
PCR反应的程序如下:95℃预变性4min,35个循环(94℃ 30s,65.2℃ 30s,72℃30
s),最后72℃延伸10min,4℃保存。
九、PCR产物纯化和DNA池测序
PCR扩增完成之后进行琼脂糖凝胶电泳检测,然后进行PCR产物的切胶回收及纯化:
在紫外灯下从琼脂糖凝胶上切下含目的片段的凝胶,放入1.5mL离心管中,然后用PCR
产物回收纯化试剂盒(北京天根生物公司)纯化PCR产物,按照试剂盒说明书操作。
把以关中奶山羊DNA池为模板的PCR纯化产物送南京金斯瑞有限公司进行双向测
序。对奶山羊PITX1基因第1内含子230bp片段的测序峰图进行分析,其中在同一位点
有两个不同峰表示发生了单核苷酸突变。本发明筛查到了奶山羊PITX1基因的1个SNP,
如图1所示,在其第201位为G或A的单核苷酸多态性,可表示为G>A,其中图1中带
框的字母表示该位点发生G>A的突变(NW_00314033:g.201G>A);绿线代表A碱基
的测序峰,红线代表T碱基的测序峰,蓝线表示C碱基的测序峰,黑线代表G碱基的测
序峰)。然后利用Blastn软件和BioXM(2.6)软件对G等位基因和A等位基因的相应序
列进行序列比对和酶切位点分析,发现G>A的突变(NW_00314033:g.201G>A)存在MspI
多态性,其分析结果如图2所示,图2为MspI PCR-RFLP方法检测山羊PITX1基因第1
内含子SNP(NW_003104033:g.201G>A or IVS1+41G>A)序列分析图,其中带框序列分别表
示上下游引物序列,粗体且下划线的字母或小点且下划则表示SNP位点。
实施例二:奶山羊PITX1基因201位G>A多态性的MspI PCR-RFLP分析
一、多态位点分析
当奶山羊PITX1基因第201位发生G>A突变时,即G突变为A,原来CCG序
列相应地变成CCG从而失去了MspI限制性内切酶识别序列,不能被限制性内切酶
MspI所切割。
二、PCR反应条件
PCR扩增体系和反应条件分别如前面所述,其中模板使用奶山羊个体基因组DNA,
PCR扩增产物用1.0%琼脂糖凝胶电泳检测,可以清晰的看到230bp的条带。
三、PCR产物酶切(MspI)及RFLP检测
对PCR扩增产物首先分别进行MspI酶切,然后根据琼脂糖凝胶电泳结果判定该位
点的单核苷酸多态性。
1、酶切:
(1)、MspI酶切体系(25μL):20μL奶山羊个体DNA PCR产物,2.5μL 10×buffer缓
冲液(含BSA),1.0μL(浓度10U/μL)的MspI限制性内切酶,1.5μL灭菌超纯水(H2O)。
(2)、酶切条件:37℃恒温水浴或空气浴中消化5h。
2、琼脂糖凝胶电泳分型:
制作2.0%的琼脂糖凝胶,100V电压下电泳,当分子量不同的DNA片段分离清晰时,
EB染色,用BIO-RAD Gel Doc 2000凝胶成像分析系统照相分析,人工进行判型。
电泳结果如图3所示,其中泳道1、2、3、5、9、14为GG基因型;泳道6、8、10、
11、12、13为GA基因型;泳道4、7为AA基因型;Marker为Marker I,分别是600bp、
500bp、400bp、300bp、200bp和100bp;由于奶山羊PITX1基因PCR扩增产物的SNP
位点所在序列为CCGG时,能被限制性内切酶MspI识别,230bp的基因片段被切为121bp、
109bp两条带;当上述该位点发生G>A的突变时,不能被限制性内切酶MspI所识别,
即片段不能被限制性内切酶所切割,仅有一条230bp大小的条带;由于山羊为二倍体动
物,当发生G>A突变时,可形成不同的基因型,分别为GG、GA和AA,可以通过其酶
切后的电泳图来进行判别,图3中AA基因型表现为230bp一条带,GA基因型表现为
121bp、109bp、230bp三条带,GG基因型表现为121bp和109bp两条带。
实施例三:本发明制备的DNA标记在关中奶山羊群体多态性中的诊断应用
一、群体多态性中的诊断:
利用上述的SNP检测方法对所有关中奶山羊的DNA样品分别进行PCR-RFLP的基
因型鉴定。
二、SNP位点的频率统计分析:
基因型频率是指一个群体中某一性状的各种基因型之间的比率。PAA=NAA/N,其中
PAA代表某一位点的AA基因型频率;NAA表示群体中具有AA基因型的个体数;N为检
测群体的总数量。
基因频率是指一个群体中某一基因对其等位基因的相对比率。计算的公式可以写成:
PA=(2NAA+NAa1+NAa2+NAa3+NAa4+......+NAan)/2N
式中,PA表示等位基因A频率,NAA表示群体中具有AA基因型的个体数量,NAai表示
群体中具有Aai基因型个体数量,a1-an为等位基因A的n个互不相同的复等位基因。
本研究所涉及的等位基因为A和G,所以具体的基因频率计算公式为:
PA=(2NAA+NAG)/2N
PG=(2NGG+NAG)/2N
式中,PA,PG分别表示等位基因A和G的频率,NAA、NGG和NAG分别表示AA、GG
和AG基因型的个体数量,N表示总群体数目。
如表2所示的基因频率分布,关中奶山羊PITX1基因MspI位点SNP中的A等位基
因频率为0.230,G等位基因频率为0.770,符合等位基因频率大于0.01(1.0%)的SNP的
定义,故本位点属于单核苷酸多态位点。
表2关中奶山羊PITX1-MspI位点等位基因频率及基因型频率
实施例四:奶山羊PITX1基因SNP位点的基因效应关联分析
一、基因型数据:限制性内切酶MspI识别的基因型(AA、AG和GG)
二、生产数据:关中奶山羊奶产量、乳脂含量(包括平均乳脂含量、晨乳脂含量和晚
乳脂含量)和乳固体物含量。
三、关联分析模型:
先对数据进行描述分析,确定是否存在离群值,再利用最小二乘分析对数据校正;
根据数据特征,应用SPSS(18.0)软件分析各基因型间的生产性状效应。在对基因型效应
进行分析时采用了固定模型:
Yijk=μ+Ai+GJ+(AG)ij+eijk,
其中:Yijk为泌乳性状的观察值,μ为总体均值,Ai为年龄效应,Gj为基因型效应,
(AG)ij为年龄和基因型的交互效应,eijk为随机误差。由于所有实验动物(关中奶山羊)都
由同一养殖场提供,因此农场效应不予考虑。
结果表明(见表3):GA基因型的关中奶山羊在晨乳脂含量、晚乳脂含量、平均乳脂
含量方面均显著高于GG和AA基因型个体(P<0.05);GA基因型奶山羊在总固体物含量
方面显著高于AA基因型个体(P<0.05),但与GG基因型个体差异不显著(P>0.05)。这说
明GA基因型可以作为一个提高奶山羊泌乳性能的候选分子遗传标记。
因此,本发明首次成功利用MspI PCR-RFLP方法检测了中国本土奶山羊PITX1基因
内含子1的单核苷酸多态性(SNP)。并通过关联分析验证,这个位点(IVS1+41G>A)的杂
合基因型(GA型)可以作为提高奶山羊泌乳性能的DNA标记。
表3关中奶山羊P1位点(PITX1-MspI)SNP与泌乳性能的关联分析
以上所述,仅为本发明的较佳实施例,并非对本发明作任何形式上和实质上的限制,
凡熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用以上所揭示的技
术内容,而作出的些许更动、修饰与演变的等同变化,均为本发明的等效实施例;同时,
凡依据本发明的实质技术对以上实施例所作的任何等同变化的更动、修饰与演变,均仍
属于本发明的技术方案的范围内。