可乐猪乳头数性状全基因组关联分析

乳头数性状是猪最重要的繁殖性状之一,与养猪业的经济效益密切相关。本试验以Illunima Porcine SNP 60K芯片对22头可乐猪进行基因分型,通过Plink 1.07软件,基于混合线性模型对左乳头数、右乳头数和总乳头数性状进行全基因组关联分析。经过严格的质量控制和多重检验之后,共鉴定出全基因组水平潜在关联的SNPs位点4个（P<2.06E-5）;染色体水平显著关联位点3个,潜在关联位点18个;在关联SNP位点可能连锁的上、下游1 cM区间内检索到304个基因,其中Wnt及Fgf信号通路中的候选基因BTRC、FGF5、FGF8和BMP3、RASGEF1B、HMGB3可能影响猪的乳头数性状或繁殖性状。通过全基因组关联分析策略发掘出的关联SNP位点及潜在候选基因,将为可乐猪的选育保种奠定基础。     

大白母猪乳头数全基因组关联分析

乳头数是评价母猪繁殖性能的重要参数,是大白猪选种的重要指标之一.本研究利用1118头大白猪母猪群体,利用50K SNP芯片对该群体进行基因分型,采用单变量混合线性模型的全基因组关联分析(GWAS)及Fst选择信号分析,筛选该群体的左乳头数、右乳头数和总乳头数性状相关基因.结果显示:GCTA软件计算左、右、总乳头数性状的...     

加系大白猪乳头数性状的全基因组关联分析

乳头数是母猪重要繁殖性状,影响断奶仔猪存活率。研究者对不同品种和品系的猪相继开展了乳头数性状相关数量性状座位(Quantitative Trait Locus,QTL)鉴定,但目前关于加系大白猪乳头数性状相关QTL的研究较少。本研究以944头加系大白猪群体为研究对象,对乳头数性状进行全基因组关联分析(GWAS)。GCTA软件计算左、右、总乳头数性状的遗传力分别为0.1、0.12、0.16。利用rMVP软件在SSC7上分别检测到4、1、4个与左、右和总乳头数基因组水平显著相关(P<1.26×10-6)的SNP,其中1个SNP位于SSC7的14.1 Mb,剩余3个SNP位于SSC7的97.6～102.0 Mb。在SSC3、SSC4、SSC6、SSC7、SSC9、SSC11、SSC13上分别鉴定到4、1、4、5、1、1、5个与乳头数性状潜在显著相关(P<2.52×10^-5)的SNP。进一步扫描25个显著相关SNP上下游200 kb位置处的基因,共获得84个候选基因,其中SSC7上获得39个候选基因,包括FAM71D、MPP5、EIF2S1、PIGH、ARG2...     

儋州鸡体尺性状全基因组关联分析

【目的】 挖掘影响地方鸡体尺性状的有效SNP位点及功能基因, 给儋州鸡育种工作提供有效的数据基础和理论支撑。【方法】 共采集200只儋州鸡血样并提取基因组DNA, 利用10×全基因组重测序技术获得全基因组SNP标记并对试验个体基因型进行分型。使用EMMAX软件基于混合线性模型对70日龄的儋州鸡体尺性状(胫长、胫围、体斜长、胸宽、髋骨宽、胸深、龙骨长)进行全基因组关联分析。【结果】 共发现与胫长性状和胫围性状基因组水平显著相关的SNPs位点有12和8个, 与胫长性状相关SNPs分别定位于1、2、4和8号染色体上; 与胫围性状相关的SNPs定位于2、4、8和13号染色体上。预测与胫长相关的候选基因为KCNA1、TPK1、EZH2、FSTL5和AMY2A基因, 与胫围相关的候选基因为TPK1、FSTL5、AMY2A、TGFBI、LECT2和IL-9。通过KEGG通路分析和GO注释发现, 8个基因参与钾离子跨膜转运、硫胺素新陈代谢、细胞增殖、钙离子结合、骨骼肌卫星细胞维持与骨骼肌再生、细胞受体相互作用、生长因子活性等生物学进程。【结论】 本研究发现了20个与儋州鸡体尺性状关联的SNPs位点, 并筛选到8个目标性状候选基因, 为儋州鸡育种提供候选的分子标记, 为地方鸡标记辅助选择提供新的思路。     

基于重测序数据的湖羊产羔数性状全基因组关联分析

旨在利用全基因组关联分析(genome-wide association study, GWAS)挖掘与湖羊繁殖性状相关的候选基因,解析湖羊高繁殖性状的遗传基础。本研究利用206只湖羊母羊的全基因二代重测序数据(群体平均测序深度为6×),结合第一胎产羔数、第二胎产羔数及2胎平均产羔数的表型性状,基于混合线性模型,并利用质控后的1 604 526个SNP标记进行全基因组关联分析。结果:未发现影响湖羊第一胎产羔数、第二胎产羔数和2胎平均产羔数的显著SNP位点,这可能与羊的繁殖性状是一种低遗传力性状、且受微效多基因控制有关。此外表明,本研究所使用的混合线性模型被验证是可信的。     

畜禽生长发育相关性状的全基因组关联分析研究进展

全基因组关联分析(GWAS)是近年兴起的用于分析复杂性状的重要研究方法。高通量测序技术的成熟发展使得基于全基因组测序技术和基因芯片技术的GWAS解析畜禽复杂性状成为可能,GWAS的运用对畜禽经济性状相关的SNP、QTL和候选功能基因研究起到关键作用。本文主要对GWAS的基本原理和方法、优劣势以及GWAS在畜禽生长发育相关性状中的应用现状进行综述,并对GWAS在今后畜禽育种中的应用前景进行展望,以期为GWAS在畜禽育种中的深入研究提供参考。     

大白公猪精液相关性状全基因组关联分析

本研究旨在通过全基因组关联分析技术挖掘大白公猪原精体积、原精密度、精子畸形率、精子总数和有效精子数等性状的关联SNP位点及候选基因.选取498头大白公猪,提取精液DNA,利用Gene Seek 50K芯片进行基因分型,通过混合线性模型计算个体随机效应值来作为全基因组关联分析表型值,利用Farm-CPU法进行全基因组关联...     

基于全基因组关联分析揭示肉鸡腿病发生的遗传机制

旨在揭示肉鸡腿病发生的遗传机制以降低肉鸡腿病发生率。采用我国自主培育的“广明2号”多品系、多世代共2 331只白羽肉鸡公鸡作为试验素材,于42日龄通过X光鉴定腿部健康状况,利用“京芯一号”55K SNP芯片对全部个体进行基因分型,使用plink运用二分类性状的逻辑回归模型(logistic regression model, LR Model)对腿部健康表型进行全基因组关联分析(genome-wide association study, GWAS)。将不同群体芯片数据合并进行质控后,保留2 330个个体(2 256只健康,74只患病)和30 414个SNPs位点用于后续分析。通过GWAS分析筛选到5个潜在SNPs位点,分布在6号和18号染色体上,分别注释到与腿病相关的3个功能候选基因TBCD、SIRT1和PBLD上。本研究为白羽肉鸡腿部健康表型提供了重要的遗传位点和候选基因,为推进肉鸡抗病育种进程提供遗传基础。     

巴马香猪产活仔数性状全基因组关联分析

为寻找与巴马香猪产活仔数相关的分子标记,试验利用全基因组关联分析（GWAS）定位并筛选了影响产活仔数性状的候选基因,采集297头具有多胎产仔记录的巴马香猪耳组织样品,提取DNA并利用猪50K SNP芯片进行基因分型,分型结果经质控与基因型填充后,使用Tassel软件对巴马香猪产活仔数性状进行全基因组关联分析。结果显示,巴马香猪平均窝产活仔数在1~9胎内随着胎次增加逐渐升高;经质控过滤后共获得32 816个SNPs位点,利用全基因组关联分析共筛选到8个与巴马香猪产活仔数相关的SNPs位点,分别在基因组或染色体水平达到显著;对关联显著SNP位点上下游500 kb内的编码基因进行富集分析,并依据猪繁殖性状相关QTL区域及基因功能,最终筛选到4个基因（CAPZB、MSH3、CITED2和HSD17B7）作为影响巴马香猪产活仔数的候选基因。     

南阳牛生长性状相关基因组区域全基因组关联分析

本研究旨在筛选和鉴定与南阳牛生长性状相关的基因组区域和候选基因,从而更好地了解牛生长性状的遗传机制。试验共采集71头南阳牛母牛血样并提取基因组DNA,利用SLAF-seq（specific-locus amplified fragment sequencing）技术获得全基因组SNP标记并对试验个体基因型进行分型。对每头个体初生重及不同月龄（6、12、18、24、36）的体重、体高、体斜长、胸围和坐骨端宽及每6个月体增重等生长性状进行全基因组关联分析;获得显著相关的基因组区域后,对其中基因进行功能注释以筛选候选基因。结果显示,共获得141 755个筛选后的SNPs,通过全基因组关联分析鉴定出5个分别与12月龄体重（8号染色体：17 320 634~17 347 720 bp）、12月龄胸围（2号染色体：15 063 190~15 155 309 bp）、24月龄体斜长（11号染色体：60 727 342~81 425 987 bp）、36月龄坐骨端宽（14号染色体：15 635 762~15 643 272 bp）和12~18月龄体增重（26号染色体：40 456 192~40 456 477 bp）等生长性状显著相关的基因组区域（LOD≥6.35）。通过对5个基因组区域内的186个基因进行功能注释,共筛选得到11号染色体上的8个基因（BMP10、IFT172、SDC1、TCF23、TRIM54、RAB1A、VPS54和GDF7）与骨生长、肌肉发育和生长调控有关,建议其可优先作为牛生长性状相关候选基因进行进一步验证。     

香猪基因组外显子SNPs检测及繁殖和体型性状候选基因分析

为从基因组水平上探究香猪性成熟早、产仔数少、体型小的分子机理,本研究下载与香猪在繁殖和体型性状方面有明显差异的梅山猪、杜洛克猪的重测序数据作为参照,对26头香猪、24头梅山猪和24头杜洛克猪的重测序数据进行SNPs检测和等位基因频率差值分析,筛选香猪基因组外显子上与梅山猪、杜洛克猪高度分化的SNPs位点,并使用AS-PCR方法和Sanger测序方法统计其中8个SNPs位点的群体等位基因频率。最后,对含有高度分化SNPs位点的基因进行GO功能富集和KEGG通路分析。结果显示,SNPs检测获得香猪基因组外显子上SNPs共236 710个,包括193个终止密码子丢失、1 238个终止密码子获得、90 945个错义SNPs和144 334个同义SNPs。得到1 046个香猪高度分化SNPs,包括429个错义SNPs、2个终止子丢失、6个终止子获得和609个同义SNPs,影响524个蛋白质编码基因;并发现X染色体上44-58 Mb的一段富含高度分化SNPs的热点区域。群体等位基因频率验证结果与重测序结果一致,鉴定了8个位点均为香猪与梅山猪、杜洛克猪高度分化的SNPs。对524个含有高度分化SNPs的基因进行GO和KEGG分析,富集到卵母细胞减数分裂、雌激素信号途径等繁殖相关通路;GO功能注释到细胞内雌激素受体信号、纤毛运动等繁殖相关条目,骨骼系统形态发生、骨骼发育、成骨细胞分化等体型相关条目。得到含有香猪高度分化SNPs的18个繁殖性状候选基因PTGFR、TRO、DNAH17、PKDREJ、KAT8、DNAI2、VDR、TEX14、QSOX1、AR、WNK3、ADCY3、SPEF1、MIGA2、SLC2A8、PSMF1、TBC1D20、IFT172和7个体型相关基因IBSP、NR6A1、HSPG2、TARS、PDZD2、FGF4、AMER1,可能是决定香猪性成熟早、产仔数少、体型小的关键基因。     

雌激素受体基因和骨调素基因对苏钟猪初生重和乳头数量的影响

本文研究了雌激素受体基因和骨调素基因对苏钟猪生长繁殖性能的影响。41头母猪120多窝产仔数记录用来分析雌激素受体和骨调素基因型对仔猪初生重、母猪乳头个数和仔猪乳头个数的影响。在母猪乳头个数方面,骨调素142/152和168/170、142/168基因型之间差异显著(P<0.05);在经产胎次中,仔猪初生重140/152和168/170、164/168、142/168、142/152基因型之间差异显著(P<0.05)。母猪乳头个数雌激素受体AA型和AB型之间差异显著(P<0.05)。     

鲁中肉羊初生体尺体重的遗传参数估计

为估计鲁中肉羊初生重、体高、体长、胸宽、胸深、胸围和管围的遗传参数,利用2018年上半年出生的611只鲁中肉羊的初生体尺体重数据,采用AI-REML (Average information restricted maximum likelihood)算法,借助DMU软件分析以产羔数为固定效应、个体加性遗传效应为随机效应的多性状动物模型。结果表明:鲁中肉羊初生重、胸宽、胸深、体高、体长、胸围和管围的遗传力分别为0.16、0.10、0.22、0.44、0.43、0.46和0.52,各性状之间遗传相关为-0.517～0.773,表型相关为-0.197～0.503。说明鲁中肉羊初生重、胸宽和胸深为低遗传力性状,体高、体长、胸围和管围为中等遗传力性状。     

拷贝数变异全基因组关联分析及数量性状基因座定位联合鉴定猪体高性状候选基因

旨在利用全基因组拷贝数变异区域（copy number variation regions,CNVRs）关联分析以及全基因组数量性状基因座（quantitative trait locus,QTLs）定位联合筛选出影响猪体高性状的候选基因。本研究利用快速检测基因组拷贝数变异软件CNVcaller对本实验室构建的大白×民猪F2代资源群体的重测序数据进行拷贝数变异检测。利用混合线性模型（mixed-linear model,MLM）将性别和胎次作为固定效应对体高性状进行拷贝数变异全基因组关联分析（CNVR-GWAS）。采用软件R/qtl进行QTL分析,并使用置换检验（permutation test,PT）进行检验。将CNVR-GWAS与QTL结果进行联合注释,结合GO富集和KEGG通路分析,对影响猪体高的位点和基因进行挖掘。利用实时荧光定量PCR（qPCR）方法验证候选基因。结果表明,本群体在全基因组范围内共有3 099个CNVRs,其中有两个CNVRs与体高性状在全基因组范围内显著相关,分别位于7号染色体的25 358 001~26 696 400 bp处（CNVR1）和54 087 201~54 090 000 bp处（CNVR2）。在混合线性模型分析的结果中发现,CNVR1拷贝数增加（P<0.01）和CNVR2拷贝数缺失（P<0.01）对猪的体高性状具有显著影响。基因组显著水平可找到2个显著影响猪体高的QTLs,分别为BH-1和BH-2,其中BH-2对体高性状的影响较大。CNVR1和BH-2重叠区存在1个嗅觉受体基因OR12D3和18个未被注释的基因。qPCR验证OR12D3的拷贝数变异与利用混合线性模型统计推断出的结果一致。初步推测,OR12D3基因的拷贝数变异可能与猪体高性状相关。