论文部分内容阅读
本研究以秦川牛、南阳牛、鲁西牛和安格斯牛4个肉用品种共计630头个体为研究对象,利用PCR-SSCP、PCR-RFLP、DNA随机个体直接测序、创造酶切位点酶切(ACRS-PCR-RPLP)、飞行质谱MALDI-TOF技术研究分析了牛AMPK家族中7个基因的遗传变异,分析了其遗传结构和遗传多样性。对7个基因的多态位点与200头秦川牛的生长及肉质性状进行了关联分析,同时应用生物信息学方法,分析了这7个基因的DNA和蛋白质序列特征,从分子水平上深入了解了它们的作用规律和机理,以期筛选对肉牛重要经济性状具有显著效应的功能基因或分子标记,为肉牛的高效选育和分子标记数据库的建立、种质资源保护和利用提供遗传学依据。本研究主要获得了以下结果：1、系统分析了牛PRKAA1、PRKAA2、PRKAB1、PRKAB2、PRKAG1、PRKAG2、 PRKAG3共7个基因序列和蛋白序列,借助生物信息学理论和方法,分析归纳了7个基因的组织表达差异,预测了各基因蛋白二级和三级结构、信号肽、跨膜螺旋、结构域特征和亚细胞定位,并将牛各基因蛋白与其它物种进行同源比较,构建了系统发生树,阐释了这一基因家族的高度保守性和在各物种起源进化过程中的重要作用。2、采用PCR-SSCP与基因测序相结合的方法,研究了秦川牛、南阳牛、鲁西牛、安格斯牛四个品种PRKAA1基因5、7、8、9、10外显子的多态性,发现了PRKAA1基因第7外显子上的213T>A多态位点与秦川牛肌内脂肪含量极显著关联(P<0.01)。3、将随机个体测序与飞行质谱技术相结合,发现存在于牛PRKAA2基因第4外显子27G>C、60T>C、第8内含子1046G>C、1188T>A四个突变。27G>C位点与秦川牛背膘厚和眼肌面积显著关联(P<0.05)；60T>C位点与秦川牛体斜长、体高、尻长、腰角宽、坐骨端宽、眼肌面积显著关联(P<0.05)；1046G>C位点与秦川牛胸深、背膘厚显著关联(P<0.05)；1188T>A位点与秦川牛坐骨端宽、背膘厚及眼肌面积显著关联(P<0.05)。单倍型和连锁不平衡分析,以上4个SNPs存在5个频率大于0.03的单倍型。27G>C与60T>C位点合并效应分析,CCTT、CGCT、CGTT、GGCC、GGCT、GGTT6种组合基因型中GGCC基因型优于其他基因型组合效应；1046G>C、1188T>A位点合并分析,CCTT、CGAT、CGTT、GGAA、GGAT、GGTT6种组合基因型中GGAT基因型优于其他基因型组合效应。4、采用PCR-SSCP方法结合直接测序技术,在秦川牛PRKAB1基因上共检测到18个多态位点。其中第一外显子129T/C多态位点与秦川牛肌内脂肪含量、眼肌面积显著关联(P<0.05)。第二外显子63T/C多态位点与秦川牛体斜长、腰高、眼肌面积显著关联(P<0.05)。第3内含子、第4外显子、第4内含子上15个SNPs中12个位点与秦川牛体高、腰高、尻长、腰角宽、胸深、胸围六个体尺性状显著关联(P<0.05), In444G/A位点多态性与眼肌面积显著关联(P<0.05)。第5内含子57C／T位点与胸围、眼肌面积极显著关联(P<0.01)。单倍型和连锁不平衡分析表明,存在于intron3、exon4、intron4上的15个SNPs中仅有5个频率大于0.03的单倍型。In423A/G与In424C/T位点合并基因型组合效应分析,AACT基因型优于其他组合基因型。5、SSCP结合测序方法在牛PRKAB2基因上共检测到10个突变。其中第三内含子上的In336C/A和In385T/G位点与秦川牛体尺和肉质性状没有关联,第五内含子上的8个SNPs中,In5130C/G、In5410T/C三个位点与秦川牛体尺性状没有关联(P>0.05),其他5个位点In560A/G、In5118A/G、In5396A/G、In5421A/C)与秦川牛体尺性状都有不同程度的关联性(P<0.05),In548T/C、In561A/G、 Ins118A/G、In5130C/G、In5421A/C位点与肌内脂肪含量显著关联(P<0.05)。第五内含子8个SNPs位点可能形成5个频率大于0.03的单倍型,In560A/G与In561A/G位点组合效应值分析,AAAA、AAAG组合基因型为最优组合型。6、测序结合创造酶切法和飞行质谱技术在秦川牛PRKAG1基因上共检测到In21174T>G和In51248T>C两个多态性位点,均为中度多态,秦川牛群体在两位点处于Hardy-Weinberg平衡(P>0.05)。 In211746T>G位点与秦川牛体尺和肉质没有相关(P>0.05),In512482T>C位点与体高、腰角宽、胸围、坐骨端宽、眼肌面积、肌内脂肪含量显著关联(P<0.05)。7、采用同样方法在秦川牛PRKAG2基因上发现了8个SNPs,其中En3223A>G突变为错义突变,引起编码氨基酸由丝氨酸到甘氨酸的改变。生物信息学分析,该位点氨基酸突变引起了PRKAG2三级结构较大变化。关联分析显示该位点与秦川牛体高、腰高、腰角宽、尻长、胸深、胸围、眼肌面积等性状显著关联(P<0.05),In516C>T位点只与秦川牛体高性状显著关联(P<0.05), En,3112T>C位点与秦川牛体高、尻长及腰角宽性状显著关联(P<0.05)。其他5个位点与秦川牛体尺和肉质性状并不相关(P>0.05)。8、采用直接测序、飞行质谱、ACRS-RFLP-PCR多种方法,在牛PRKAG3基因上共检出45个SNPs位点,其中编码区10个,3’UTR区18个,部分内含子上17个。关联分析表明,4个SNPs位点与秦川牛体斜长显著相关,8个位点与体高显著相关,21个位点与腰高和尻长相关,17个位点与腰角宽和胸深相关,6个位点与胸围相关,5个位点与坐骨端宽相关,2个位点与肌内脂肪含量相关,7个位点与眼肌面积相关,3个位点与背膘厚相关。