论文部分内容阅读
我国短脂尾绵羊品种资源丰富,分布广泛于中北部区域。据考古学和遗传学研究,蒙古羊是短脂尾绵羊品种的共同祖先,在经过迁移、驯化和适应之后,许多演化品种在表型性状上发生了明显变化,特别是在繁殖、生长等性状上差异较大,但由进化所导致表型差异的遗传机制尚不清楚,缺乏全基因组水平的系统研究。本研究以蒙古羊、小尾寒羊和多浪羊为研究对象,采用全基因组重测序、生物信息学分析等技术方法开展了相关研究。取得结果如下：(1)本研究通过全基因组重测序,总共获得三个绵羊品种基因组测序数据78 G,深度达到已公布的绵羊参考基因组的43.2×,平均每个品种14.4×；与参考基因组比对,共获得13,209,060个高质量SNPs(单核苷酸多态性)和912,131个InDels(插入缺失位点),通过Sanger测序验证SNP准确率达到95.8%(115/120)。(2)位于小尾寒羊蛋白编码区(CDS)的SNP共90,853个,其中同义SNP有36,820个,错义SNP有54,033个(分别位于11,063个基因内),无义SNP有2,474个(分别位于1,789个基因内)。位于蒙古羊CDS的SNP共90,720个,其中同义SNP有37,449个,错义SNP有53,271个(分别位于11,421个基因内),无义SNP有2,421个(分别位于1,810个基因内)。位于多浪羊CDS的SNP共56,804个,其中同义SNP有21,029个,错义SNP有35,775个(分别位于7,724个基因内),无义SNP有1,786个(分别位于1,192个基因内)。(3)由于小尾寒羊和多浪羊都具有多胎、常年发情的表型,而蒙古羊为单胎、季节性发情,因此我们选取了同时在小尾寒羊和多浪羊中发生错义突变而不在蒙古羊中发生错义突变的基因共797个,进行GO功能注释,主要富集于6大类生物学功能,其中15个基因富集于繁殖条目(其中12个基因富集于GO：0007283精子发生),KEGG通路分析共注释到5条生化通路,其中Progesterone-mediated oocyte maturation(孕激素介导的卵母细胞成熟)通路包含7个基因(RPS6KA3、MAD2L1、CCNB2、GNAI2、ADCY5、 PIK3R5、CDC25B),我们推测这些基因与短脂尾绵羊的繁殖性状有关。(4)我们在小尾寒羊基因组上共鉴定出选择进化区域145个,总长度45.3 Mb,占整个基因组长度的1.74%。其中常染色体上,以Z(HP)<-4和Z(FST)>4为阂值,有108个高度纯合的区域和65个与蒙古羊的遗传距离极远的区域。合并后共136个区域,包含774个蛋白编码基因。X染色体上,以Z(Hp)<-3和Z(FST)>3为阈值,共有9个区域包含24个蛋白编码基因被鉴定。(5)对小尾寒羊选择进化区域内的基因进行GO功能注释,主要富集到8大类生物学功能(发育过程、细胞过程、多细胞生物过程、生物调节、代谢过程、繁殖、定位和生长)的32个功能条目。KEGG通路分析注释到11个生化通路,其中包括Endometrial cancer(子宫内膜癌)、Prostate cancer(前列腺癌)等,可能与短脂尾绵羊的繁殖性状具有潜在的关系。(6)小尾寒羊基因组选择进化区域内繁殖功能相关的基因共25个,其中只在小尾寒羊发生进化的有8个：MMP14、PMCH、HMGCR、SLC6A3、PRL、OXT、TSSK3、FLT1；同时在小尾寒羊和多浪羊中发生进化的有5个:ADAM2、GAMT、NR3C1、DAZAP1、 ADAM29;只在蒙古羊发生进化的有1个：HOXA10;三个品种都发生进化的有5个：DMRTB1、RQCD1、FGF9、ABCB9、YBX2。