西藏牦牛群体遗传多样性及遗传结构分析

采集西藏地区5个牦牛类群（西藏高山牦牛,类乌齐牦牛,帕里牦牛,斯布牦牛,娘亚牦牛）共49头牦牛的血样,利用全基因组重测序技术对其遗传多样性及遗传结构进行了分析,以期为西藏牦牛遗传资源的保护、开发和利用提供参考。SNP变异检测结果表明,经过滤后得到7 655 923个SNP位点;群体遗传多样性分析发现,西藏地区不同牦牛群体之间的基因杂合度差异并不显著,且均存在不同程度的杂合度下降现象,类乌齐牦牛的观测杂合度最高,为0.336 4;群体结构分析发现,5个牦牛类群之间存在基因交流现象。说明目前西藏牦牛各群体整体存在遗传多样性下降的潜在风险,需要进一步加强群体保种工作,以确保不同品种资源的多样性得到有效保护。     

牦牛MDHⅠ基因多态性及其与生长性状的关联分析

本研究旨在探讨牦牛细胞质苹果酸脱氢酶（cytoplasmic malate dehydrogenase,MDHⅠ）基因的遗传多态性及其与生长性状的关联性,以期为牦牛优良性状选育提供参考依据。选取5个牦牛类群（共178头）,采集耳组织样提取DNA并构建DNA池,对MDHⅠ基因所有外显子设计特异性引物进行扩增测序,用Mega 5.2筛查SNPs位点,直接测序法鉴定其基因型,利用PopGene 32进行χ²独立性检测、基因纯合度（Ho）、基因杂合度（He）、有效等位基因频率（Ne）和多态信息含量（PIC）分析,利用SPSS 24.0分析MDHⅠ基因与牦牛生长性状间的关联性。结果表明,申扎牦牛、类乌齐牦牛、斯布牦牛、帕里牦牛和麦洼牦牛均在外显子8区域发现1个突变位点（G23093C）,属同义突变,有3种基因型：GC、CC和GG,5个牦牛类群中,GC基因型为优势基因型（类乌齐牦牛中GG为优势基因型）,G为优势等位基因（斯布牦牛除外）,在斯布牦牛中,G和C基因频率相等。χ²适应性检验表明,申扎牦牛、类乌齐牦牛、斯布牦牛、帕里牦牛和麦洼牦牛均符合Hardy-Weinberg平衡（P>0.05）。MDHⅠ基因在申扎牦牛、帕里牦牛和斯布牦牛中属于高度多态（PIC>0.5）,在麦洼牦牛和类乌齐牦牛中属于中度多态（0.25 < PIC < 0.5）;纯合度最高的是申扎牦牛和类乌齐牦牛。关联性分析结果表明,MDHⅠ基因不同基因型对类乌齐牦牛、帕里牦牛和麦洼牦牛体重有显著性影响（P<0.05）;且种间分析发现,MDHⅠ基因不同基因型对牦牛体重有显著性影响（P<0.05）,说明该基因可作为人工选育的分子标记。     

CTGF、Gprin3基因多态性及其与牦牛体尺性状的关联分析

【目的】 检测牦牛结缔组织生长因子（CTGF）和G蛋白调节诱导神经突增生家族亚型因子3（Gprin3）基因的单核苷酸多态性（SNP），并分析多态位点与牦牛体尺性状的关联性。【方法】 以西藏的斯布牦牛、帕里牦牛、类乌齐牦牛、申扎牦牛及四川的麦洼牦牛5个类群260头个体为研究对象，利用基因池筛选技术检测CTGF、Gprin3基因SNP位点，分析其多态信息含量、有效等位基因数、纯合度及杂合度等指标，利用最小二乘线性模型分析不同基因型与牦牛体尺性状的关联性。【结果】 牦牛CTGF基因存在3个SNPs位点，分别为T1537A、C2195T和C2421T，均位于内含子区，其中T1537A位点在斯布牦牛和麦洼牦牛中均处于Hardy-Weinberg平衡状态，而在帕里牦牛、类乌齐牦牛和申扎牦牛中均偏离Hardy-Weinberg平衡状态，该位点与牦牛体高呈显著相关（P<0.05）；C2195T位点在5个牦牛类群中均处于Hardy-Weinberg平衡状态，且处于中度多态（0.25<PIC<0.5）；C2421T位点在5个牦牛类群中均处于Hardy-Weinberg平衡状态，该位点与牦牛体重、体斜长、胸围、管围和前肢长均呈显著相关（P<0.05）。牦牛Gprin3基因共发现4个SNPs，分别为G310C、C414T、C1100T和G1213A，其中G310C、C414T及C1100T位点在5个牦牛类群中均处于Hardy-Weinberg平衡状态；G1213A位点在类乌齐牦牛中偏离Hardy-Weinberg平衡状态，在其余4个类群中均处于Hardy-Weinberg平衡状态；G310C位点在5个牦牛类群中均处于低度多态（PIC<0.25），C414T和C1100T位点在类乌齐牦牛、斯布牦牛中均处于中度多态，G1213A位点在类乌齐牦牛中处于中度多态，在其余牦牛类群中为低度多态；C414T位点与牦牛体重、体高、体斜长和胸围均呈显著相关，而在斯布牦牛和申扎牦牛中C1100T位点与其胸围和前肢长均呈显著相关（P<0.05）。【结论】 牦牛CTGF和Gprin3基因具有丰富的多态性，CTGF基因T1537A、C2421T位点及Gprin3基因C414T、C1100T位点均与牦牛体尺性状有关，可作为影响生长发育性状的候选基因及基因座用于牦牛定向选育。     

西藏牦牛的RAPD遗传多样性及其分类研究

为了解西藏地区牦牛品种或类群的遗传多样性和亲缘关系,本研究从33个RAPD多态性引物中筛选出8个条带清晰且多态性丰富的引物对西藏地区的巴青牦牛、类乌齐牦牛、丁青牦牛、桑日牦牛、工布江达牦牛、江达牦牛、康布牦牛、桑桑牦牛、嘉黎牦牛、帕里牦牛、斯布牦牛等11个类群的核基因组DNA进行了RAPD分析,并用Nei氏标准距离和UPGMA聚类法分析了类群间的亲缘关系.结果表明:(1)西藏牦牛类群的遗传多样性指数变异范围在0,185 7～0.405 3之间,其中帕里牦牛最小(0.185 7),说明相对较纯,群体较整齐;而工布江达牦牛最大(0.405 3),显示该群体内部具有较多的遗传变异.(2)在11个类群中,其遗传多样性指数大小分别为:工布江达牦牛(0.405 3)＞江达牦牛(0.353 6)＞斯布牦牛(0.344 8)＞康布牦牛(0.342 8)＞嘉黎牦牛(0.332 3)＞桑日牦牛(0.282 3)＞巴青牦牛(0.279 3)＞桑桑牦牛(0.269 8)＞丁青牦牛(0.259 7)＞类乌齐牦牛(0.224 1)＞帕里牦牛(0.185 7),具有西藏东部牦牛类群遗传多样性相对较高,而西部牦牛类群遗传多样性相对较低的趋势,预示着西藏东部可能是牦牛的起源地之一.(3)遗传距离构建的分子聚类关系图表明:西藏11个牦牛类群可分为2大类,帕里牦牛(PL)为一类,其余10个牦牛类群为另一类.综上所述,西藏牦牛具有较丰富的遗传多样性,品种或种群内的遗传分化显著,这是西藏牦牛业持续发展和牦牛适应外界环境的遗传基础,是将来培养牦牛新品种或品系的重要基因资源;西藏牦牛品种可分为2大类群.     

青海省同德牦牛群体的父系遗传多样性及遗传背景探析

【目的】为从分子水平上揭示青海省同德牦牛的父系遗传多样性、群体遗传结构和遗传背景。【方法】本研究对32头同德公牦牛使用5个Y-SNPs标记（SRY4、USP9Y、UTY19、AMELY3和OFD1Y10）和1个Y-STR标记（INRA189）进行PCR扩增、测序和分型,使用BioEdit、Arlequin和Network等生物信息学软件综合分析同德牦牛的父系遗传多样性、群体遗传结构及父系起源。【结果】32头同德牦牛5个Y-SNPs标记即SRY4（969bp）、USP9Y（470bp）、UTY19（290bp）、AMELY3（971bp）和OFD1Y10（763bp）的PCR扩增产物测序长度与先前研究结果一致,INRA189标记分型分析共检测到155bp、157bp和159bp 3个等位基因;在同德牦牛群体AMELY3标记中检测到一个新的Y-SNP位点（即g.719 C>T）;基于5个Y-SNPs标记11个Y-SNPs位点和INRA189标记3个等位基因的联合分型,共确定了6种Y染色体单倍型（即Y1H1、Y1H2、Y1H3、Y1H4、Y1H5和Y2H6）,Y染色体单倍型多样度（Hd）为0.714±0.060,表明同德牦牛具有丰富的父系遗传多样性;系统发育分析显示同德牦牛由2个父系支系组成（即Y1和Y2）,提示其拥有2个父系起源。【结论】同德牦牛拥有特殊的父系遗传信息,具有丰富的父系遗传多样性,由2个父系支系组成,拥有2个父系起源。     

西藏牦牛mtDNA cytb基因的序列多态性及其系统进化分析

为从分子水平上探究西藏牦牛类群的亲缘关系、分类地位和遗传多样性,本研究测定了嘉黎牦牛、桑桑牦牛、桑日牦牛、工布江达牦牛、斯布牦牛、帕里牦牛、康布牦牛、江达牦牛、类乌齐牦牛、丁青牦牛、巴青牦牛等11个西藏牦牛类群共110头牦牛的细胞色素b基因全序列,分析了其多态性,并构建了11个类群的系统进化树.结果表明:11个西藏牦牛类群的细胞色素b基因全序列长均为1140 bp,共有单倍型53种,其中新发现的有49种,序列间共有14个SNPs多态位点,核苷酸变异类型包括转换和颠换,无插入和缺失,以同义突变为主,说明西藏牦牛具有较丰富的遗传多样性.西藏11个牦牛类群可分为:帕里牦牛系、江达牦牛系、巴青牦牛系、桑日牦牛系、类乌齐牦牛系等5大系.     

湖羊的遗传多样性及保种效果评价研究

本研究采用微卫星标记技术,评估了湖州市5个保种场的湖羊遗传多样性,并对保种效果进行了分析。结果显示,5个湖羊群体间的等位基因数差异不显著,目前的保种方法保留了湖羊等位基因;5个湖羊群体各微卫星位点观测杂合度（Ho）差异极其显著,其中3个群体遗传多样性较高;期望杂合度（He）差异不显著;多态信息含量（PIC）差异不显著,且均值都高于0.5,处于高度多态,表明群体具有丰富的遗传多样性和选择潜力。进一步计算了5个湖羊群体的遗传距离,进行了聚类分析,构建了各群体的进化树,结果表明,5个群体的湖羊家系数均达到保种要求,其中4个群体杂合个体较多,可能受外源基因影响。本研究对湖羊的遗传多样性及保种效果进行了有效评估,并提出了存在的问题,有利于湖羊的种质保护和开发利用。     

玉树牦牛与四个牦牛群体线粒体DNA D-Loop区遗传多样性分析

利用线粒体DNA标记方法从分子水平上研究牦牛的遗传多样性。根据牦牛线粒体DNA(mtDNA)序列设计引物,扩增出长度为1 000 bp左右的片段,序列对比分析。结果表明:mtDNA D-Loop区长度为945 bp,共发现变异位点127个,单倍型71个,单倍型多样度(Hd)为0.909±0.016,核苷酸多样性(pi)为0.082;其中玉树牦牛发现变异位点74个,单倍型27个,单倍型多样度(Hd)为0.885±0.034,核苷酸多样性(pi)为0.0134;申扎牦牛单倍型多样度、核苷酸多样性和平均核苷酸差异数均最高;5个牦牛群体D-Loop区序列具有一定的A/T碱基偏好性。玉树牦牛不遵循中性进化(P0.05),其余4个牦牛群体符合中性进化(P0.05);玉树牦牛与类乌齐牦牛和帕里牦牛的遗传距离较近(0.023、0.018),申扎牦牛与斯布牦牛的遗传距离最远(0.533)。发现玉树牦牛与其余4个群体的牦牛在同一个分支。斯布牦牛和申扎牦牛群体出现两个分支,说明玉树牦牛在进化的过程中可能不存在相互交流的情况。     

牦牛MDHⅠ基因多态性及其与生长性状的关联分析

本研究旨在探讨牦牛细胞质苹果酸脱氢酶(cytoplasmic malate dehydrogenase,MDHⅠ)基因的遗传多态性及其与生长性状的关联性,以期为牦牛优良性状选育提供参考依据。选取5个牦牛类群(共178头),采集耳组织样提取DNA并构建DNA池,对MDHⅠ基因所有外显子设计特异性引物进行扩增测序,用Mega 5.2筛查SNPs位点,直接测序法鉴定其基因型,利用PopGene 32进行χ~2独立性检测、基因纯合度(Ho)、基因杂合度(He)、有效等位基因频率(Ne)和多态信息含量(PIC)分析,利用SPSS 24.0分析MDHⅠ基因与牦牛生长性状间的关联性。结果表明,申扎牦牛、类乌齐牦牛、斯布牦牛、帕里牦牛和麦洼牦牛均在外显子8区域发现1个突变位点(G23093C),属同义突变,有3种基因型:GC、CC和GG,5个牦牛类群中,GC基因型为优势基因型(类乌齐牦牛中GG为优势基因型),G为优势等位基因(斯布牦牛除外),在斯布牦牛中,G和C基因频率相等。χ~2适应性检验表明,申扎牦牛、类乌齐牦牛、斯布牦牛、帕里牦牛和麦洼牦牛均符合Hardy-Weinberg平衡(P0.05)。MDHⅠ基因在申扎牦牛、帕里牦牛和斯布牦牛中属于高度多态(PIC0.5),在麦洼牦牛和类乌齐牦牛中属于中度多态(0.25PIC0.5);纯合度最高的是申扎牦牛和类乌齐牦牛。关联性分析结果表明,MDHⅠ基因不同基因型对类乌齐牦牛、帕里牦牛和麦洼牦牛体重有显著性影响(P0.05);且种间分析发现,MDHⅠ基因不同基因型对牦牛体重有显著性影响(P0.05),说明该基因可作为人工选育的分子标记。     

运用SNP芯片评估马身猪保种群体的遗传结构

旨在研究马身猪保种群体的遗传多样性、亲缘关系和家系结构。本研究利用Illumina CAUPorince 50 K SNP芯片检测39头保种马身猪的单核苷酸多态性（single nucleotide polymorphism，SNP）；采用Plink软件计算最小等位基因频率、多态信息含量、观察杂合度和期望杂合度，分析保种群体的遗传多样性；采用Plink软件构建状态同源（identity by state，IBS）距离矩阵和分析连续性纯合片段（runs of homozygosity，ROH），采用Gmatix软件构建G矩阵，分析保种群体的亲缘关系；采用Mega X软件构建群体进化树，分析保种群体的家系结构。结果显示，39头保种马身猪中共检测到43 832个SNPs位点，平均基因型检出率为0.980 1；通过质控的SNPs位点有28 859个，其中72.4%具有多态性，该款SNP芯片适用于分析马身猪的遗传多样性。有效等位基因数为1.563 4，多态信息含量为0.412，最小等位基因频率为0.258，表明马身猪保种群体的遗传多样性比较丰富；平均观察杂合度为0.354 1，平均期望杂合度为0.349 9，说明马身猪保种群体出现了分化；平均IBS遗传距离为0.284 2，其中公猪为0.285 2，IBS距离矩阵和G矩阵结果均表明部分种猪之间存在亲缘关系；ROH共有8 131个，其中46.15%的长度在400～600 Mb之间，平均近交系数为0.237，说明保种群体的近交程度高；群体进化树结果表明，马身猪保种群体来源于3个家系，各家系的个体数量差异明显。马身猪保种群体的遗传多样性较丰富，但近交程度高，家系少，各家系的个体数量差异大，容易引起遗传多样性的丢失，因此，需从原种场引入新的血统，扩大保种群体数量，降低近交系数。     

肃南牦牛群体遗传结构、选择信号分析和ROH检测

【目的】试验旨在对肃南牦牛进行群体遗传结构、选择信号分析和连续纯合片段(ROH)检测,挖掘肃南牦牛的种质特性相关基因。【方法】选择肃南牦牛、巴州牦牛、斯布牦牛、九龙牦牛和天祝白牦牛5个牦牛品种共计48头牦牛进行全基因组重测序,采用基因组变异检测流程(GATK)获得高质量的单核苷酸多态性(SNP)标记进行下游分析;对5个牦牛品种的基因型数据进行主成分分析、祖先成分分析和系统进化树分析以确定其群体结构;对5个牦牛群体进行ROH检测以及近交系数计算;采用综合单倍型评分(iHS)方法筛选共有高频ROH区域内受选择位点。【结果】5个牦牛品种共鉴定出15 092 883个SNPs,主要分布于内含子区域。亲缘关系计算结果显示,所有个体均不存在三代以内亲缘关系,满足后续分析要求。主成分分析表明,5个牦牛品种可以显著地分为5个类群,其中肃南牦牛群体内遗传变异较小。群体结构分析显示,肃南牦牛包含其他4种牦牛祖先成分,其中天祝白牦牛祖先成分所占的比例最大。ROH分析显示,5个牦牛品种共检测到8 426个ROHs片段,其中高频ROH区域421个。相比于其他4个牦牛品种,肃南牦牛的ROH片段总长度和数目最多,具有较高的连锁程度,其近交系数远大于其他牦牛群体。在肃南牦牛高频ROH区域上注释得到465个候选基因,主要富集到与机体发育、大脑形态形成、脂肪氧化相关的通路,包括胚胎骨骼系统形态发生、前后模式规范、甲状腺发育通路和Hippo通路,其中与胚胎发育及组织分化过程相关的基因有同源框A3(HOXA3)、HOXA5、HOXD3,与肉品质相关的基因有花生四烯酸12B (ALOX12B)、ALOX15B、ALOXE3,与中脑发育相关的基因为成纤维细胞生长因子8(FGF8)。在7号染色体1个高频ROH区域中(Chr7:12 661 870-13 045 935)鉴定到3个共享基因(核内不均一性核糖核蛋白K (HNRNPK)、驱动蛋白27(KIF27)、G激酶锚定蛋白1(GKAP1))与牦牛共有的高原适应性有关。对共有高频ROH区域内的位点进行iHS分析,鉴定到的受选择基因主要与抗病性、内质网分泌蛋白加工及细胞周期调节相关,包括N-α乙酰转移酶25(NAA25)、内质网分子伴侣29(ERP29)、跨膜蛋白16(TMEM116)、TRAF型锌指结构域蛋白1(TRAFD1)、HECT结构域E3泛素连接酶4(HECTD4)基因。【结论】本研究从全基因组水平系统评估肃南牦牛的遗传多样性和遗传背景,鉴定了肃南牦牛ROH区域内与表型相关的基因,并重点挖掘了与其他牦牛品种共享高频ROH区域内的受选择基因,为肃南牦牛种质资源开发、利用提供了重要理论依据。     

牦牛EPO、PPARα基因SNP与高原低氧适应的相关性分析

试验旨在探究促红细胞生成素基因（EPO）、过氧化物酶体增殖剂激活受体α（PPARα）的遗传多样性,并分析其多态性与牦牛高原低氧适应性的相关性。采集不同海拔高度的6个牦牛类群（中甸牦牛、麦洼牦牛、斯布牦牛、类乌齐牦牛、帕里牦牛、申扎牦牛）以及三江黄牛共375头耳样,提取DNA并分别构建DNA池,采用直接测序法结合PCR-RFLP检测分析EPO、PPARα基因的多态性,最后应用SHEsis软件统计分析候选基因SNPs与牦牛高原适应性的相关性。结果表明,EPO基因存在3个SNPs位点：rs527G→A、rs1031A→T、rs1192T→C;PPARα基因存在3个SNPs位点：rs77363C→T、rs77471C→A和rs77534C→T。χ²适合性检验结果显示,EPO基因3个SNPs位点均符合Hardy-Weinberg平衡状态;PPARα基因的rs77363C→T位点上,6个牦牛类群都处于平衡状态（P>0.05）,在PPARα基因的rs77471C→A和rs77534C→T位点,麦洼牦牛偏离Hardy-Weinberg平衡（P<0.05）。单倍型分析得出,EPO基因的ATC单倍型在高海拔地区牦牛中的分布频率随海拔升高而升高;PPARα基因的TAC单倍型在6个牦牛类群中分布频率显著高于其他单倍型。研究表明,EPO、PPARα基因可作为牦牛适应高原环境的分子标记,为进一步探讨牦牛高原低氧适应性机制提供一定理论帮助。     

基于GBS简化基因组测序评估3个麦洼牦牛保种群的遗传结构研究

麦洼牦牛是青藏高原地区优良的乳肉兼用型牦牛品种,本研究旨在探究四川省龙日种畜场麦洼牦牛保种群的遗传多样性和遗传结构,评价3个不同保种群的保种效果并挖掘重要种质特性基因。对麦洼牦牛3个保种群粉嘴群（n=140）、全黑群（n=211）、弗洛群（n=55）进行GBS简化基因组测序,基于检测到的126122个单核苷酸多态性（SNPs）标记计算遗传统计量,结果表明,整个牦牛群的平均观测杂合度（Ho）和平均期望杂合度（He）为0.3038和0.3036,麦洼牦牛的遗传多样性较丰富。全黑群、粉嘴群、弗洛群的观测杂合度Ho分别为0.3029、0.3042、0.3044,近交系数Fis分别为0.0144、0.0152、0.0209,弗洛群和粉嘴群受人工选择的强度大于全黑群,较低的近交水平说明3个群的保种效果较好。Structure分析中全黑、弗洛群部分个体血缘较纯正,而其他个体血缘关系非常混杂;粉嘴群和全黑群的遗传分化系数（Fst）和遗传距离（DR）最大为0.03513、0.0358,结合系统进化树表明两者亲缘关系最远,有遗传分化趋势。利用Fst和π法对3个保种群进行选择信号分析,发现有104个受选择基因广泛参与生殖机能、免疫系统、胚胎发育、脂质代谢等条目以及生殖激素、内/外分泌、信号传递等通路,其中部分基因提示麦洼牦牛的繁殖、肉质、毛色性状以及应激反应得到了人工选择,如PPP3CC、KCNMA1、ROCK2、GNAQ、MEF2C、KIT等。现有的麦洼牦牛保种策略是可行的,研究结果为未来麦洼牦牛的保种选育和遗传改良提供了参考依据。     

近交陆川猪的遗传稳定性研究

旨在探讨高度近交陆川猪的遗传稳定性。本试验采集20头近交和6头非近交陆川猪的耳组织样品进行个体基因组重测序,通过扫描全基因组SNPs分析近交陆川猪群体的遗传多样性,并基于基因组纯合度挖掘近交陆川猪群体中特异的纯合区域（region of homozygosity,ROHs）及相关基因。结果表明,近交陆川猪群体的观察杂合度（0.373）小于期望杂合度（0.491）,推测其发生了选择或者近交;系统进化树和群体遗传结构分析发现,近交陆川猪群体在近交过程中分化成为两个近交家系,且230、236、239、132、130及126号个体亲缘关系最近;在近交群体中,基因组纯合度和基因组长纯合片段比例的均值分别为59.64%和0.12,都极显著高于非近交群体53.18%和0.05（P<0.01）,说明该近交陆川猪群体已达到高度近交;近交群体和非近交群体中分别有1 250和633个ROHs,重叠分析结果显示,特异存在于近交群体中的阳性纯合区域（positive region of homozygosity,pROHs）有224个,其中由于近交产生的纯合突变可能影响参与激素合成过程、胰液分泌、膜脂代谢过程、IgA产生的肠道免疫网络及纹肌细胞分化等相关的1 322个基因。结果提示,近交陆川猪群体的遗传多样性低、基因组纯合度高且已达到高度近交,说明该群体的遗传稳定性高。此外,近交陆川猪群体在近交过程中产生的纯合突变可能与机体的免疫、生长等生产性状相关。     

亚丁牦牛和拉日马牦牛遗传多样性及遗传结构分析

亚丁牦牛和拉日马牦牛均为肉乳兼用型优良地方牦牛资源,本研究旨在探究亚丁牦牛和拉日马牦牛的遗传多样性及各牦牛群体间的亲缘关系,对亚丁牦牛、拉日马牦牛以及其他7个地方牦牛群体（九龙牦牛、麦洼牦牛、金川牦牛、昌台牦牛、中甸牦牛、玉树牦牛、类乌齐牦牛）进行RAD简化基因组测序,基于检测到的SNP信息计算遗传统计量。结果表明,亚丁牦牛和拉日马牦牛的观测杂合度分别为0.2186、0.2233,亚丁牦牛遗传多样性较贫乏。亚丁牦牛的平均遗传分化指数Fst值最高为0.0653,与其他8个牦牛群体存在中度分化,Structure分析显示其血统构成纯正,可列为一个独立遗传资源。拉日马牦牛Fst值最低为0.0443,其群体遗传结构较复杂,与九龙牦牛和麦洼牦牛存在基因交流,属麦洼牦牛的祖先群之一。亚丁牦牛和拉日马牦牛均值得进一步研究和保护,研究结果对牦牛遗传资源的保护与利用具有重要意义。     

Screening the Polymorphic Tri-nucleotide Repeat Microsatellites from the Genome of Maiwa Yak and Their Genetic Diversity Analysis

In order to study the genetic diversity and differentiation characteristics of the existing population of Maiwa yak and provide a theoretical basis for utilization in the future, a total of 9 466 tri-nucleotide repeats microsatellites loci were screened based on the release of genome of yaks by bioinformatics method. The tri-nucleotide repeats microsatellites were amplified by PCR using the genomic DNA of 191 Maiwa yaks as templates,and the polymorphic microsatellite loci were genotyped through polyacrylamide electrophoresis experiment. After the Hardy-Weinberg equilibrium test, 5 tri-nucleotide microsatellite loci were deviated from the Hardy-Weinberg equilibrium (P <0.05). 15 microsatellite loci were highly polymorphic (PIC >0.5), the mean observed heterozygosity (Ho),mean expected heterozygosity(He) and mean polymorphic information content(PIC) were 0.5041,0.8152 and 0.7831,respectively. Maiwa yak population exhibited no characteristics of differentiation by Structure 3.0 software analysis. Therefore, Maiwa yak had plentiful genetic polymorphism while no signs of genetic differentiation were observed within the population.     

麦洼牦牛基因组三碱基重复微卫星的挖掘及其遗传多态性分析

为了解麦洼牦牛现有群体的遗传多样性及种群遗传分化特征,试验通过生物信息学方法从已公布的牦牛基因组中筛选到9 466个三碱基重复微卫星位点,以191个麦洼牦牛基因组DNA样品为模板对各个位点进行PCR扩增和聚丙烯酰胺凝胶电泳基因分型。经过哈迪-温伯格平衡检验,在15个三碱基微卫星位点中,有5个位点偏离哈迪-温伯格平衡（P<0.05）。15个微卫星位点均为高度多态性（PIC>0.5）,群体平均观测杂合度（Ho）、平均期望杂合度（He）和平均多态信息含量（PIC）分别达到0.5041、0.8152和0.7831。进一步通过Structure 3.0软件分析发现,麦洼牦牛群体结构单一,无分化特征,表明麦洼牦牛品种虽然具有较高的遗传多态性但是无分化迹象。