类乌齐牦牛全基因组RNA编辑位点预测及剖析

为揭示牦牛不同组织间受RNA编辑导致转录产物发生的差异变化,进而为牦牛的组织分化、系统遗传调控等研究提供候选基因。本研究以3头4.5岁类乌齐母牦牛的大脑、小脑、臀部脂肪以及肌肉组织作为试验材料,通过Illumina 4000测序平台对各组织表达产物进行扫描,利用SPRINT和JACUSA筛选差异的RNA编辑位点,分析RNA编辑位点的分类、注释以及变异风险评估等。结果发现了总计24 784个RNA编辑位点,其中在4个组织中均发现存在编辑事件的有4 015个位点;4个组织中的RNA编辑位点主要以A-G和T-C编辑类型为主;发现其中一个高风险编辑位点位于SON基因中,使该基因的翻译提前终止,这将导致该组织一些特定位点的RNA与DNA结合出现问题。本研究预测并分析了类乌齐牦牛大、小脑和臀部肌肉、脂肪组织中RNA编辑差异位点的类型和分布,为进一步研究牦牛组织分化和生长发育中重要的调控基因提供了新的参考。     

牦牛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),说明该基因可作为人工选育的分子标记。     

牦牛GPR41基因遗传多态性与生长性状关联性研究

本研究旨在探讨G蛋白偶联受体41(GPR41)基因在牦牛中的遗传多态性,对麦洼牦牛(50头)、申扎牦牛(47头)和类乌齐牦牛(28头)GPR41基因外显子区域片段进行扩增,筛查SNPs位点,利用RFLP分子标记技术鉴定其基因型,开展其与生长性状间的关联分析,为牦牛品种改良和选育提供遗传学依据。结果表明:在麦洼牦牛中未发现突变位点,申扎和类乌齐牦牛在外显子2区域发现1个突变位点,第4570位碱基发生G→A突变,属同义突变;经RFLP鉴定出3种基因型,命名为AA、AG和GG;χ~2适应性检验显示,申扎牦牛基因型分布不符合Hardy-Weinberg平衡(P<0.01),类乌齐牦牛处于Hardy-Weinberg平衡(P>0.05);申扎牦牛和类乌齐牦牛均表现为中度多态;GPR41基因不同基因型与申扎牦牛体重和体高显著相关,且GG型显著高于AA型(P<0.05),而GPR41基因不同基因型对类乌齐牦牛各体尺性状无显著性影响(P>0.05)。综上,申扎牦牛GPR41基因G4570A基因座可为牦牛选育过程中分子标记选择提供参考。     

牦牛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位点均与牦牛体尺性状有关，可作为影响生长发育性状的候选基因及基因座用于牦牛定向选育。     

牦牛酪蛋白基因家族的克隆及组织表达分析

本研究旨在克隆牦牛酪蛋白基因家族(CSN1S1、CSN1S2、CSN2和CSN3)的CDS区序列,鉴定其在牦牛不同组织中的表达水平。选取4岁龄左右处于泌乳期的健康类乌齐母牦牛3头,屠宰后分别采集乳腺、心脏、肝脏、骨骼肌组织,分别提取组织总RNA并反转录为cDNA,设计酪蛋白基因家族特异性引物扩增酪蛋白基因家族序列,进行生物信息学分析,并利用实时荧光定量PCR法分别检测酪蛋白家族基因mRNA水平。结果显示,克隆得到CSN1S1、CSN1S2、CSN2和CSN3基因cDNA序列分别为919、832、805和715bp,其CDS区全长分别为645、669、690和585bp,分别编码214、222、259和194个氨基酸残基。类乌齐牦牛酪蛋白基因家族与黄牛亲缘关系最近,其次是印度水牛,而与单胃动物猪的亲缘关系最远。组织表达结果显示,酪蛋白基因家族在组织中广泛表达,其中在乳腺组织中的表达量最高,其次是骨骼肌组织。在乳腺组织中CSN1S1、CSN1S2、CSN2基因之间表达量差异不显著(P0.05),但CSN2基因表达量显著高于CSN3基因(P0.05)。以上结果为酪蛋白基因家族在牦牛乳腺蛋白质代谢调控机制的研究提供了参考依据。     

肉鸡腹脂双向选择系中RNA编辑位点的鉴定研究

旨在准确获得肉鸡腹脂双向选择系脂肪组织中的RNA编辑位点信息,进而为脂肪相关研究提供更多的信息来源。本研究以东北农业大学肉鸡高、低腹脂双向选择品系第十九世代肉鸡为试验材料,使用转录组测序(RNA-seq)和基因组测序数据识别肉鸡腹部脂肪组织RNA编辑位点。通过生物信息学手段,设定严格的数据过滤与筛选策略,构建RNA编辑位点候选集。结果表明,在候选集中,通过对RNA编辑位点的分类和注释,在低脂系肉鸡脂肪组织中发掘转换型RNA编辑位点28个,颠换型RNA编辑位点101个,插入缺失型RNA编辑位点71个;在高脂系肉鸡脂肪组织中,发掘转换型RNA编辑位点30个,颠换型RNA编辑位点84个,插入缺失型RNA编辑位点77个。绝大多数RNA编辑位点落在内含子、基因间区域及基因上游区域(内含子区41.91%,基因间区域33.08%,基因下游20.59%,基因上游2.94%和外显子区1.47%)。试验验证了4个候选RNA编辑位点,其中3个位点落在与脂肪形成有重要关系的ATP1B3、SLC6A6和FLNB基因区域。本研究鉴定和验证了存在于肉鸡脂肪组织中的RNA编辑位点,为进一步研究RNA编辑影响脂肪组织生长发育的分子机制奠定了基础。     

牦牛酪蛋白基因家族的克隆及组织表达分析

本研究旨在克隆牦牛酪蛋白基因家族(CSN1S1、CSN1S2、CSN2和CSN3)的CDS区序列,鉴定其在牦牛不同组织中的表达水平。选取4岁龄左右处于泌乳期的健康类乌齐母牦牛3头,屠宰后分别采集乳腺、心脏、肝脏、骨骼肌组织,分别提取组织总RNA并反转录为cDNA,设计酪蛋白基因家族特异性引物扩增酪蛋白基因家族序列,进行生物信息学分析,并利用实时荧光定量PCR法分别检测酪蛋白家族基因mRNA水平。结果显示,克隆得到CSN1S1、CSN1S2、CSN2和CSN3基因cDNA序列分别为919、832、805和715bp,其CDS区全长分别为645、669、690和585bp,分别编码214、222、259和194个氨基酸残基。类乌齐牦牛酪蛋白基因家族与黄牛亲缘关系最近,其次是印度水牛,而与单胃动物猪的亲缘关系最远。组织表达结果显示,酪蛋白基因家族在组织中广泛表达,其中在乳腺组织中的表达量最高,其次是骨骼肌组织。在乳腺组织中CSN1S1、CSN1S2、CSN2基因之间表达量差异不显著(P>0.05),但CSN2基因表达量显著高于CSN3基因(P<0.05)。以上结果为酪蛋白基因家族在牦牛乳腺蛋白质代谢调控机制的研究提供了参考依据。     

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

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

类乌齐牦牛SIRT3基因克隆与生物信息学及差异表达分析

本实验旨在丰富牦牛 SIRT3基因的基础研究,通过对西藏类乌齐牦牛 SIRT3基因克隆、生物信息学分析及差异表达分析,进一步探讨 SIRT3基因的生理功能及分子机制。结果表明:类乌齐牦牛 SIRT3基因CDS区长 1 002 bp,编码 333个氨基酸,与普通牛 SIRT3基因 CDS区比对存在碱基突变,编码氨基酸发生错义和同义突变,同义突变发生在碱基第 49、279、342、384 位,错义突变发生在碱基第 149 和 620 位,其编码蛋白属亲水性蛋白,有多个跨膜区域和磷酸化位点,无信号肽片段,主要位于线粒体,二级结构以无规卷曲和α-螺旋为主;SIRT3基因编码蛋白在催化活性、核苷酸结合及辅因子结合等过程发挥主要功能,主要参与生物体代谢过程;类乌齐牦牛 SIRT3基因序列与普通牛一致性较高,亲缘关系最近;QPCR分析显示,SIRT3基因在牦牛肝脏中的表达水平极显著高于其他组织(P<0.01)。     

m6A甲基化酶METTL3和WTAP在牦牛组织及脂肪细胞中的表达

本研究旨在探究m⁶A RNA甲基化酶METTL3和WTAP在牦牛(Bos grunniens)不同组织、前体脂肪细胞增殖与分化过程中的表达模式和前体脂肪细胞分化过程mRNA m⁶A的变化水平。采用qRT-PCR检测牦牛皮下脂肪、肌肉、心、肝、脾、肺、肾和皮下脂肪不同时期(18和30月龄)METTL3及WTAP的mRNA表达水平。应用I型胶原酶消化法获取牦牛前体脂肪细胞,油红O染色和脂肪分化标志基因的检测建立牦牛前体脂肪细胞分化模型,以及qRTPCR检测前体脂肪细胞增殖分化阶段METTL3和WTAP的mRNA表达水平。结果表明,肝脏组织中METTL3表达最高(P<0.05),皮下脂肪组织表达量最低(P<0.05);WTAP在皮下脂肪组织中的表达最为丰富,脾脏中的表达量最低(P<0.05)。30月龄皮下脂肪组织中METTL3和WTAP mRNA表达量高于18月龄。前体脂肪细胞诱导分化12 d时,细胞中出现多而密的脂环,脂肪细胞分化特异性标志基因FABP4、C/EBPα和PPARγ第12天的表达量显著高于第0天(P<0.05)。METTL3和WTAP的表达量在细胞增殖阶段(24、48和72 h)呈现“下降-上升”的表达趋势(P<0.05)。在牦牛前体脂肪细胞分化阶段(0、4、8和12 d),METTL3表达量呈现“上升-下降-上升”的趋势,WTAP呈现“上升-下降”的趋势。细胞分化阶段mRNA m⁶A水平呈现逐渐上升的趋势,在分化12 d时细胞内RNA的m⁶A丰度最高(P<0.05)。本研究获得METTL3和WTAP在牦牛不同组织和前体脂肪细胞增殖分化阶段的变化规律及细胞分化过程中mRNA m⁶A水平变化,初步揭示WTAP和METTL3对牦牛脂代谢具有重要的调控作用。     

九龙牦牛PPARγ基因的克隆及其表达谱分析

根据普通牛(Bos tarus)过氧化物酶体增殖物激活受体γ(PPARγ)基因序列设计引物,用成年九龙牦牛(Bos grunniens)脂肪组织总RNA,经RT-PCR扩增获得了PPARγ基因序列(GenBank登陆号:GU061328),其中cDNA的ORF为1 428 bp,编码475个氨基酸,与普通牛PPARγ氨基酸的同源性达99%,有2个氨基酸发生突变。利用半定量RT-PCR分析九龙牦牛PPARγ基因的mRNA表达特性。结果表明:在脂肪、背最长肌、心、肝、肾、脾和肺脏中均检测到PPARγ基因的表达,并且在脂肪组织中表达量极显著高于其他组织(P0.01),在肝脏和脾脏中亦有较高表达。PPARγ基因在背最长肌中的表达5.5岁九龙牦牛显著高于0.5、3.5岁和9岁以上,其表达与背最长肌的肌内脂肪含量未见显著相关性。     

香猪3种组织环状RNA的表达及其功能分析

本研究旨在探索环状RNA（circular RNA，circRNA）在香猪皮肤、皮下脂肪及背最长肌的表达及其潜在功能。本研究采用RNA-Seq技术和生物信息学方法，分析香猪皮肤、皮下脂肪及背最长肌circRNA的表达，对3种组织特有circRNA亲本基因进行GO和KEGG富集分析，对样品组织结构功能相关亲本基因中circRNA结合的靶miRNA进行预测。从香猪3种组织表达谱中共鉴定出6 483个circRNAs，皮肤、皮下脂肪及背最长肌分别鉴定出4 575、5 180、5 219个circRNAs，其中特异性表达的分别有83、234、586个；circRNAs在染色体上分布广泛，主要来源于外显子，少数来源于内含子和基因间隔区；多数circRNAs表达量较低；皮肤特有circRNAs的亲本基因主要参与蛋白水解、RNA转运、缝隙连接等过程，其中NF1、MAPK1等亲本基因与皮肤性疾病相关，可能以circRNAs的形式发挥作用，其中亲本基因MAPK1中的circRNA（circ-MAPK1）可结合miR-18a、miR-132、miR-411等，可能参与调控胶原蛋白的形成；皮下脂肪特有circRNAs的亲本基因富集于脂肪细胞生长、发育及癌症等信号通路，ABHD5、PTPN11等亲本基因与脂肪代谢有关，circ-PTPN11结合的miR-103、miR-107、miR-199a-5p等参与调控脂肪细胞增殖、分化及脂质沉积；背最长肌特有的circRNAs亲本基因富集于癌症、感染、肌肉发育等过程，ROCK2、PPP1CC等基因与肌细胞分化及肌肉收缩相关，circ-PPP1CC结合的miR-339、miR-181a、miR-181b等参与调控肌细胞分化、增殖及生长。综上，本研究筛选出的circRNAs及其结合的miRNAs可能参与皮肤胶原蛋白形成、脂肪沉积、肌细胞分化成熟等生物学过程。     

Expression and Function Analysis of Circular RNA in Three Tissues of Xiang Pig

The purpose of this study was to explore the expression and potential functions of circular RNA (circRNA) in the skin,subcutaneous fat,and longissimus dorsi muscle of Xiang pig.In this study,RNA-Seq technology and bioinformatics methods were used to analyze the expression of circRNAs in skin,subcutaneous fat,and longissimus dorsi muscle of Xiang pig,GO and KEGG enrichment analysis were performed for the host genes of three tissue-specific circRNAs,then target miRNA prediction was performed for the circRNAs that host genes associated with tissue structure and function.A total of 6 483 circRNAs were identified from the three tissues of Xiang pig,4 575,5 180 and 5 219 circRNAs were identified respectively from skin,subcutaneous fat and longissimus dorsi muscle respectively,among which 83,234 and 586 circRNAs expressed specifically.circRNAs distributed on chromosomes widely and mainly from exons,a few from introns and intergenic regions.The expression levels of most circRNAs were low.The host genes of skin-specific circRNAs were mainly involved inproteolysis,RNA transport,and gap junctions.Among them,host genes such as NF1 and MAPK1 were related to cutaneous diseases and might play a role in the form of circRNAs.The circRNA (circ-MAPK1) that the host gene MAPK1 could bind to miR-18a,miR-132,miR-411,and which might be involved in regulating the formation of collagen.The host genes of circRNAs expressed specifically in subcutaneous fat were enriched in adipocyte growth,development and cancer signaling pathways.ABHD5 and PTPN11 genes were related to fat metabolism,miR-103,miR-107 and miR-199a-5p bound by circ-PTPN11 could be involved in regulating adipocyte proliferation,differentiation and lipidosis.The host genes of circRNAs only expressed in longissimus dorsi muscle associated with cancer,infection,and muscle development.ROCK2 and PPP1CC were involved in myocyte differentiation and muscle contraction,while miR-339,miR-181a and miR-181b bound by circ-PPP1CC were involved in the regulation of myocyte differentiation,proliferation and growth.In conclusion,the circRNAs selected in this study and their combined miRNAs might be involved in skin collagen formation,fat deposition,and muscle cell differentiation and maturation.     

类乌齐牦牛UCP3基因克隆与序列分析

利用RT-PCR方法成功克隆UCP3基因,对其核苷酸序列及氨基酸序列进行生物学在线分析。结果表明,类乌齐牦牛UCP3基因CDS区序列长为936bp,A、G、C和T平均含量分别为20.83%、28.63%、31.09%和19.44%;与普通牛UCP3基因CDS区序列一致性高达99%;共发现碱基突变位点6个,其中同义突变4个,反义突变2个。UCP3基因共编码311个氨基酸,分子质量为34.18ku,理论等电位点9.53;其主要分布在线粒体内膜上,为非分泌蛋白;具有α-螺旋、无规则卷曲、延伸链、β-转角等结构。由系统进化树可知,类乌齐牦牛与野牦牛最接近,与普通牛次之。     

猪CA5B基因的序列特征和表达分析

旨在通过分析猪CA5B(线粒体碳酸酐酶 VB,carbonic anhydrase 5b,mitochondria;CA5B/CA VB/Car5b)基因的序列特征和时空表达,解析CA5B对猪组织器官发育和代谢的影响.本研究采集6头240日龄贵州猪(公、母各3头)的心、肝、脾、肺、肾、皮下脂肪、背最长肌、睾丸、卵巢,采...     

利用微卫星遗传标记评估5个西藏牦牛群体遗传多样性及群体结构

【目的】 评估当前西藏主要牦牛群体遗传多样性，梳理不同地区群体间遗传结构，明确西藏5个牦牛群体（阿里牦牛、斯布牦牛、娘亚牦牛、类乌齐牦牛和帕里牦牛）的保种情况和种群间系统发育关系。【方法】 利用13个微卫星标记（SSR）对5个牦牛群体共计195个个体进行基因分型，并对各群体的等位基因数量、基因杂合度、多态信息含量及群体间遗传距离等遗传参数进行评估。【结果】 阿里牦牛群体等位基因数最多（6.43），类乌齐牦牛等位基因数最少（5.00）；观测杂合度范围为0.5311（娘亚牦牛）~0.5995（类乌齐牦牛）。各群体内偏离Hardy-Weinberg平衡的位点数为5（类乌齐牦牛）~9（阿里牦牛）个；群体内近交系数最高为0.172（阿里牦牛），且4个群体（阿里牦牛、娘亚牦牛、斯布牦牛和帕里牦牛）存在显著近交风险（P<0.05）。从遗传结构来看，所有群体间均为显著遗传分歧（P<0.05），STRUCTURE分析结果显示，5个牦牛群体划分为3个簇，其中阿里牦牛较其他牦牛群体具有更为丰富的遗传背景。系统发育树结果表明，不同群体间系统发育关系相对独立，且与种群栖息地分布不一致。主坐标分析（PCoA）结果显示，不同群体间部分个体存在较近亲缘关系，表明不同群体间存在遗传物质交流。【结论】 5个西藏牦牛群体具有丰富的遗传多样性水平，且种群系统发育关系相对独立，但多数群体存在群体事件风险。本研究不仅有助于明确西藏牦牛地方种群遗传多样性水平，同时可为今后的保种策略提升提供理论参考。