长江刀鲚 GHRH基因遗传多态性及其与生长性状的关联分析

生长激素释放激素（GHRH）是鱼类生长激素合成和分泌的重要调控因子,对细胞的生长、增殖和分化极为重要。基于长江刀鲚转录组测序数据,利用基因克隆和PCR扩增技术,获得其 GHRH基因的cDNA全长和基因组DNA序列,并开展基因遗传多态性与生长性状的关联性研究,为其分子标记辅助育种提供技术支持。结果显示：长江刀鲚 GHRH基因cDNA全长1 329 bp,由5′-UTR（300 bp）、ORF（465 bp）和3′-UTR（564 bp）组成,推导的蛋白序列由154个氨基酸组成,具有一个信号肽和GLUCA结构域;该基因的DNA序列长2 983 bp,由4个内含子和5个外显子组成,经筛选,在第1、2内含子和第3外显子上分别筛选到9、2和1个SNP位点;12个突变位点共36种基因型,分别与100尾长江刀鲚的体质量、体长和体高性状进行关联分析,获得2个与其生长性状显著相关的位点（g.1636C>T和g.1882A>C）（P<0.05）,CC基因型均是优势基因型;利用12个SNP位点分析刀鲚养殖群体的遗传多样性,获得的观测杂合度（HO）、期望杂合度（HE）和多态信息含量（PIC）分别为0.270~0.620、0.338~0.490和0.281~0.370;哈迪-温伯格平衡（HWE）检测结果显示,12个SNP位点中有4个位点显著偏离哈迪-温伯格平衡（P<0.05）。结果表明,长江刀鲚 GHRH基因上位点g.1636C>T和g.1882A>C可以作为优选长江刀鲚繁育亲本的重要候选分子标记。     

牦牛DKK1基因多态性及其与生长性状的相关性

选取麦洼牦牛、类乌齐牦牛、申扎牦牛、帕里牦牛和斯布牦牛5个群体共287头个体,采用直接测序法检测牦牛DKK1基因的SNP位点,分析其与生长性状的相关性。结果表明,牦牛DKK1基因存在2个突变位点,g.983A→G位点和g.1075C→G位点;位点g.983A→G在斯布牦牛中处于Hardy-Weinberg不平衡状态,其余位点处于平衡状态。在遗传多态性研究中,2个SNP位点均为中度多态。2个位点不同基因型与体高、体斜长、胸围、管围和体质量的相关分析结果显示,位点g.983A→G与体高、体斜长、胸围和体质量均具有相关性;位点g.1075C→G与生长性状不相关。     

大别山牛PAX3基因多态性及其与生长性状的关联分析

【目的】探究PAX3基因多态性与大别山牛生长性状的相关性。【方法】采集292头安徽大别山成年（24～48月龄）母牛的耳组织,并测定其体尺数据,提取DNA后,通过PCR和测序技术鉴定出大别山牛群体中PAX3基因的多态性,利用POPGENE、Haploview和SHEsis软件,分别对多态位点进行遗传多样性、哈代 温伯格平衡和单倍型分析,利用SPSS 23.0软件分析其与生长性状的相关性。【结果】在大别山牛PAX3基因第4内含子中检测到3个SNPs位点（g.4718G>C、g.4744T>C和g.4757C>A）,第6内含子中检测到1个SNPs位点（g.78920C>T）,4个SNPs位点均存在3种基因型。遗传多样性分析发现,g.4718G>C位点属于低度多态（PIC≤0.25）,g.4744T>C、g.4757C>A和g.78920C>T位点均属于中度多态（0.25A外,其余3个位点均处于哈代-温伯格平衡状态（P>0.05）。单倍型分析发现,4个SNPs位点间无强连锁相关（r²<0.33）,形成13种单倍型,其中频率大于0.03的单倍型有6个。相关性分析发现,g.4718G>C位点与大别山牛的十字部高显著相关(P<0.05);g.4744T>C位点与十字部高、腹围和腰角宽极显著相关(P<0.01);g.4757C>A位点与腰角宽显著相关(P<0.05);g.78920C>T位点与腹围极显著相关(P<0.01),与管围显著相关(P<0.05)。【结论】PAX3基因第4内含子的g.4718G>C、g.4744T>C 和g.4757C>A位点以及第6内含子的g.78920C>T位点多态性与大别山牛群体部分生长性状具有显著或极显著的相关性,可以作为大别山牛遗传选育的候选分子标记。     

欧拉型藏羊UCP2和UCP3基因多态性及其与生长性状的关联分析

【目的】探索欧拉型藏羊解偶联蛋白2（UCP2）和解偶联蛋白3（UCP3）基因的多态性及其与生长性状的关联性,为运用分子标记辅助育种方法提高欧拉型藏羊的生长性状提供理论依据。【方法】以239头成年欧拉型藏羊为试验对象,采用PCR产物直接测序技术检测UCP2和UCP3基因的遗传多态性,并将其与体质量、体高、体斜长和胸围生长性状进行关联性分析。【结果】在欧拉型藏羊UCP2基因上筛选出2个SNPs位点,分别为g.52130414 C>T和g.52130584 G>A,其中g.52130414 C>T位点产生TT、CT和CC 3种基因型,g.52130584 G>A位点产生AA、AG和GG 3种基因型;在UCP3基因上筛选出1个SNP位点g.52157335 C>T,产生CC、CT和TT 3种基因型。基因多态性分析表明,3个SNPs位点均为错义突变位点,属于中度多态,且均符合Hard-Weinberg平衡适应性检验（P>0.05）。关联分析结果显示,UCP2基因g.52130414 C>T位点TT基因型欧拉型藏羊体质量显著高于CT基因型(P<0.05);g.52130584 G>A位点AA基因型欧拉型藏羊体质量和体斜长均极显著高于AG和GG基因型(P<0.01),AA和GG基因型欧拉型藏羊胸围均显著高于AG基因型（P<0.05）;UCP3基因g.52157335 C>T位点CC基因型欧拉型藏羊体质量显著高于CT基因型(P<0.05),TT基因型欧拉型藏羊胸围显著高于CT基因型(P<0.05)。3个多态位点基因型组合分析发现,5个发生频率较高的双倍型中,D5型欧拉型藏羊体质量极显著高于D1、D2、D3型(P<0.01),体高显著高于D3型(P<0.05),体斜长显著高于D1、D3、D4型(P<0.05)且极显著高于D2型(P<0.01)。【结论】UCP2和UCP3基因多态性与欧拉型藏羊的部分生长性状显著相关,UCP2和UCP3基因可作为欧拉型藏羊生长性状选育的候选基因。D5双倍型为优势基因型,可在选育工作中优先考虑。     

牦牛SMAD4基因SNPs检测及其与生长性状的关联分析

【目的】研究牦牛SMAD4基因SNPs与生长性状的关系,探寻与牦牛生长性状相关的分子标记。【方法】采用DNA测序和单倍型分型技术,对青海牦牛SMAD4基因进行SNPs检测及其基因分型、连锁不平衡和单倍型分析,并对多态位点的基因型及组合单倍型与牦牛生长性状的关联性进行分析。【结果】牦牛SMAD4基因在内含子区域存在6个突变位点,其中g.393A>G、g.555A>G、g.6525A>G和g.18360C>T均存在3种基因型,g.582A>T和g.6492C>T均存在2种基因型。连锁不平衡分析发现,6个位点间不存在强连锁不平衡效应。单倍型分析发现,在14种不同的单倍型中,单倍型H₁的发生频率最高。关联性分析表明,g.393A>G位点与体斜长、胸围和管围均显著相关（P<0.05）,g.555A>G位点与体斜长显著相关（P<0.05）,g.582A>T位点与体质量和体高均显著相关（P<0.05）,g.6492C>T位点与体高显著相关（P<0.05）,g.6525A>G位点与体质量显著相关（P<0.05）,g.18360C>T位点与胸围显著相关（P<0.05）。基因型组合发现,单倍型H₁H₁₄可能是影响牦牛生长性状的最优组合。【结论】牦牛SMAD4基因内含子区域上的6个位点与生长性状显著关联,可作为牦牛分子标记辅助选择的候选基因。     

无角牦牛FTO基因多态性与生长性状的关联性分析

【目的】探究无角牦牛脂肪和肥胖相关基因(FTO)多态性与生长性状之间的关联性,以期找到与无角牦牛生长相关的分子标记。【方法】参考GenBank中野牦牛的FTO基因序列,针对其第4内含子和第8内含子设计引物,采用PCR-SSCP和DNA测序技术对354头无角牦牛的FTO基因进行多态性分析,并运用SPSS 20.0软件分析各基因与体质量、体高、体斜长和胸围等生长性状的相关性。【结果】在无角牦牛FTO基因的第4内含子区域检测出A164506C突变位点,有2个等位基因A和C,有3种基因型(AA、AC和CC),优势等位基因和基因型分别是A和AA,其频率分别为0.694 9和0.485 9;在第8内含子区域检测出G309000A突变位点,共有2个等位基因A和G,有3种基因型(GG、GA和AA),优势等位基因和基因型分别是G和GG,其频率分别为0.926 3和0.858 4。A164506C突变位点为中度多态位点,纯合度较低,处于Hardy-Weinberg平衡状态;G309000A突变位点为低度多态位点,纯合度较高,处于Hardy-Weinberg不平衡状态。A164506C位点与无角牦牛的体高和体斜长具有显著或极显著相关性,AC基因型个体优于AA基因型个体。G309000A位点与无角牦牛的胸围和体斜长具有显著或极显著相关性,GA基因型个体优于GG基因型个体。【结论】FTO基因第4内含子区域A164506C突变和第8内含子区域G309000A突变均与无角牦牛生长性状具有相关性。     

牛MyoG基因单核苷酸多态性分析

以皖东牛和广丰牛共128头牛为研究对象,采用PCR-RFLP方法检测牛MyoG基因的多态性。结果表明,在牛MyoG基因的外显子466碱基处发生G>C突变,检测到AA和AB 2种基因型,该位点的突变未引起氨基酸序列的变化。在皖东牛群体中AA基因型频率为23.66%,AB基因型频率为76.34%;在广丰牛群体中AA基因型频率为17.24%,AB基因型频率为82.76%。皖东牛的多态信息含量（PIC）为0.3606,广丰牛为0.3702,在该位点均处于中度多态。卡方检验表明皖东牛群体处于哈代-温伯格不平衡状态(P<0. 05);广丰牛处于哈代-温伯格平衡状态(P>0. 05)。     

3个安徽地方猪品种TAP1基因编码区遗传多态性分析

以3个安徽地方猪品种皖南黑猪、圩猪和霍寿黑猪共324头猪为研究对象,采用PCR-SSCP方法检测猪TAP1基因编码区的遗传多态性。结果表明,在皖南黑猪、圩猪和霍寿黑猪TAP1基因的第10外显子41 bp碱基处发生G>A突变,检测到AA、GA、GG 3种基因型,其中AA基因型频率在这3个群体中分别为59.53%、41.27%和54.91%,为优势基因型;皖南黑猪、圩猪和霍寿黑猪的多态信息含量（PIC）分别为0.2991、0.3611和0.3242,在该位点均处于中度多态。卡方检验表明,皖南黑猪、圩猪和霍寿黑猪群体均处于哈迪-温伯格平衡状态(P>0.05)。     

朝鲜鹌鹑OCA2基因克隆及其多态性与羽色的相关性

【目的】克隆朝鲜鹌鹑(Korean quail)OCA2基因的cDNA序列,检测其在不同羽色(白羽和栗羽)群体胚胎发育不同阶段翅组织中mRNA的表达情况,探讨OCA2基因多态性对朝鲜鹌鹑羽色的影响,筛选可能与羽色相关的多态性位点。【方法】取白羽、栗羽朝鲜鹌鹑孵化6,8,10,12和14 d的胚胎翅组织,提取总RNA,RT-PCR扩增OCA2基因cDNA序列,利用实时荧光定量PCR检测OCA2基因mRNA在朝鲜鹌鹑胚胎发育不同阶段和不同羽色鹌鹑胚胎翅组织中的表达水平;用胚胎发育10 d时的cDNA样品扩增OCA2基因cDNA序列并进行生物信息学分析。以300羽白羽朝鲜鹌鹑和栗羽朝鲜鹌鹑孵化10 d时的胚胎为研究材料,采集翅组织样品,提取基因组DNA, PCR扩增OCA2基因,通过不对称PCR-SSCP方法结合测序确定基因型,分析不同基因型对朝鲜鹌鹑羽色形成的影响。【结果】朝鲜鹌鹑OCA2基因胚胎发育6～8 d时的表达量低于胚胎发育10～14 d时的表达量,且在10 d时达到最高。6～8 d时,OCA2基因的表达水平在2种羽色朝鲜鹌鹑胚胎中差异不显著; 10～14 d时,栗羽朝鲜鹌鹑OCA2基因mRNA的表达水平显著高于白羽朝鲜鹌鹑(P0.05)。OCA2基因CDS区片段大小为2 660 bp,含有24个外显子,开放阅读框长2 556 bp,编码含851个氨基酸的稳定疏水性蛋白。在朝鲜鹌鹑OCA2基因中检测到3个SNP位点,分别为外显子2上的c.93AG、外显子9上的c.999AG、3′UTR上的c.2569CT;群体遗传学分析显示,3个突变位点均表现为中度多态(0.25PIC0.5),χ~2检验表明朝鲜鹌鹑群体在3个SNP位点上均处于Hardy-Weinberg平衡(P0.05)。在3′UTR上的c.2569CT位点,栗羽朝鲜鹌鹑3种基因型(CC、 CT和TT)的频率分布与白羽朝鲜鹌鹑存在极显著差异(P0.01),说明该突变位点与鹌鹑羽色之间存在一定的关联性。【结论】白羽朝鲜鹌鹑翅组织中OCA2基因低表达可能与其羽毛白化有一定的关系,OCA2基因3′UTR上的c.2569CT突变与朝鲜鹌鹑羽色性状之间具有一定的关联性,可以作为研究朝鲜鹌鹑羽色的候选基因。     

秦川牛 MYH8基因多态性及其与生长性状的关联分析

为研究秦川牛MYH8基因的多态性及其与秦川牛生长性状的关联性,选取163头健康的30月±2月龄秦川牛,采用DNA测序技术发现MYH8基因外显子上的SNP位点,并以PCR-RFLP技术分析秦川牛群MYH8基因型,研究其遗传多样性,再通过SPSS 19.0软件分析SNP位点与秦川牛体尺性状和胴体性状的关联性。研究结果表明,在MYH8基因第33外显子上存在3个SNP位点,分别在基因序列26 064、26 094和26 100bp 3处,且26 064bp处(TC)为同义突变,26 094bp处(TC)为错义突变,26 100bp处(CG)为错义突变。关联分析显示,C26100G突变发生在体斜长指标,表现为GG基因型个体的群体均值显著高于GC、CC基因型个体(P0.05);而其他突变中不存在显著性差异。某些基因型(T26064C突变处TC基因型和T26094C突变处CT基因型)可能为优势基因型,但还需通过更多个体验证。说明:MYH8基因对秦川肉牛生长性状存在一定影响。     

草鱼柠檬酸合酶基因SNP筛选及与生长性状的关联分析

根据草鱼(Ctenopharyngodon idella)EST-SNP库中预测的柠檬酸合酶基因的1个Contig序列,设计引物扩增该基因的部分序列,然后采用PCR产物直接测序法筛选SNP突变点。测序的序列经比对后,共筛选到2个在内含子10上的SNP突变点:A-386G和C-499G。随机选择同批繁殖和同塘混养的144尾草鱼对2个SNP位点用Snapshot方法进行检测和分型,并统计基因型频率。A-386G位点的AA基因型占47.10%,AG基因型占38.41%、GG基因型占14.49%;C-499G位点的CC基因型占31.85%,CG基因型占46.67%,GG基因型占21.48%,均属于中度多态位点。利用一般线性模型(GLM)分析2个SNP位点与草鱼6个生长性状(体质量、体长、体高、头长、尾柄长和尾柄高)的相关性,结果显示,A-386G和C-499G位点的不同基因型在这6个生长性状均值有差异,但均不存在显著差异。将这2个SNP不同基因型两两组合,一共组成7种双倍型(去掉频率小于3%的组合D2和D5),GLM相关分析表明5种双倍型在5个生长性状(体质量、体长、体高、尾柄长和尾柄高)上均存在差异,其中在头长上存在显著差异,D6双倍型(A-386G/AG和C-499G/GG)在6个生长性状上均值均最大,属于生长优势双倍型。本研究显示,筛选到的柠檬酸合酶基因上的2个SNP位点具有用于草鱼生长性状的分子辅助育种的潜力。     

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

【目的】对秦川牛、鲁西牛、鲁西与西门塔尔杂种牛(简称鲁杂牛)、南阳牛、夏南牛和郏县红牛6个品种中国黄牛Somatostatin基因(SST)的第1外显子和第2外显子进行SNPs检测,研究其与6个品种中国黄牛部分生长性状的相关性。【方法】采用DNA直接测序技术和PCR-SSCP方法,对669头中国黄牛SST基因第1和第2外显子进行多态性分析,运用最小二乘线性模型对SST基因的不同基因型与6个品种中国黄牛生长性状进行关联性分析。【结果】在SST基因第126个碱基处找到1处多态位点,除秦川牛和鲁杂牛有GG和AG 2种基因型外,其他4个品种牛中均存在AA、AG、GG 3种基因型。运用卡方检验模型对6个品种黄牛中等位基因A/G的频率、基因型频率、多态信息含量等进行了统计分析,结果表明,各供试品种在这一位点上均处于Hardy-Weinberg平衡状态(P>0.05);该位点的A、G等位基因频率和多态信息含量在6个品种中国黄牛中分别为0.053 8～0.218 8,0.781 2～0.946 2和0.096 6～0.283 4。SST基因的多态性与6个品种黄牛的8个生长指标的相关性分析结果显示,该位点对669个个体中的某些生长性状指标具有显著影响;其中24月龄秦川牛品种中,GG型个体的体斜长、腰角宽和坐骨端宽均高于AG型个体(P<0.05);12月龄鲁西牛品种中,GG型个体的体斜长均高于AA型个体(P<0.05);12月龄鲁杂牛牛群体中,GG型个体的体斜长和体高均高于AG型个体(P<0.05);其他品种牛各基因型个体之间的生长性状指标差异均不显著。【结论】SST基因可能与秦川牛、鲁西牛和鲁杂牛的生长性状有密切的相关性,可作为秦川牛、鲁西牛和鲁杂牛生长性状的候选基因。     

青海高原型牦牛GHR基因多态性与生长发育的关联性

【目的】探索青海高原型牦牛（Bos grunniens）生长激素受体（growth hormone receptor）基因GHR多态性及其与生长性状的关联性。【方法】以440头同等放牧条件下的30~36月龄健康牦牛为试验群体,PCR扩增其GHR基因,测序后用DNASTAR 7.1软件分析单核苷酸多态性（SNP）位点,研究不同基因型及其合并基因型与生长发育指标（体质量、体高、体斜长、胸围和管围）的关联性。【结果】青海高原型牦牛GHR基因上存在g.732091C>G、g.732195A>G和g.732373G>A 3个SNP位点,其中g.732195A>G上存在2种基因型（AA,AG）,g.732091C>G和g.732373G>A上各存在3种基因型（分别为CC、CG、GG和GA、AA、GG）;卡方检验表明,g.732091C>G、g.732195A>G和g.732373G>A均处于Hardy-Weinberg极度不平衡状态,且g.732195A>G为低度多态（PIC＜0.25）,g.732091C>G和g.732373G>A为中度多态（0.25G、g.732195A>G和g.732373G>A能够显著或极显著影响高原型牦牛的体质量和胸围,优势基因型分别为g.732091C>G和g.732373G>A的GG以及g.732195A>G的AA。对合并基因型分析发现,合并基因型 GG-AA-GG个体的生长发育尤其是体斜长表现最佳。【结论】青海高原型牦牛GHR基因有3个SNP位点,因其与青海高原型牦牛的部分生长性状相关,可作为牦牛分子标记辅助选择的候选基因。     

同羊ETAA1基因InDel检测及其与生长性状的关联分析

为寻找同羊中与生长性状相关的分子标记,试验对ETAA1基因InDel的多态性与同羊生长性状进行了关联分析.以166只无亲缘关系的同羊为研究对象,通过Ensembl数据库筛选ETAA1基因潜在的InDel位点,利用PCR和琼脂糖凝胶电泳对筛选到的InDel位点进行检测和分型后,进一步对InDel位点不同基因型的同羊生长性状进行关联分析.结果 显示,位于绵羊3号染色体的ETAA1基因第1内含子存在一个插入12 bp的InDel突变;根据插入的类型,可以分为II,DD和ID 3种基因型;检测出的InDel位点多态性对同羊公羊臀端高有显著影响,该位点Ⅱ基因型个体臀端高显著高于ID型与DD型(P＜0.05).因此,可以把ETAA1基因此位点InDel多态作为同羊臀端高性状选择的候选分子标记.     

山羊SKIP基因的分离、鉴定及与生长性状的关联分析

为研究SKIP基因对山羊生长发育的影响,并寻找与山羊生长性状相关的分子标记,利用绵羊与牛SKIP基因的同源序列设计引物,RT-PCR扩增得到一段长1 350bp的山羊SKIP基因表达序列,包括其完整编码区序列,编码449个氨基酸的蛋白。通过序列Blast比较和系统发育树的构建,发现山羊SKIP与绵羊SKIP序列相似性最高(99.52%)。利用生物信息学分析SKIP蛋白的理化性质,预测SKIP蛋白分子质量为52ku,等电点为5.18,包含典型的疏水区域和磷脂酰肌醇5-磷酸磷酸酶结构域,定位在细胞质(60.9%)和细胞核(26.1%)中。间接免疫荧光试验验证了SKIP的细胞定位,并发现位于细胞核中的SKIP在成肌细胞分化过程中转移到细胞质。通过波尔山羊和麻城黑山羊的DNA混池测序,发现SKIP内含子2有一处G/A突变,利用PCRRFLP方法对波尔×麻城黑山羊F2代杂交后代进行基因分型,性状关联分析发现该单核苷酸多态性位点(SNP)与出生体质量、胸围、管围、身高等性状显著相关,等位基因G与生长呈正相关。     

威宁黄牛SOCS4基因多态性及其与生长性状的关联分析

[目的]了解威宁黄牛SOCS4基因多态性及其与生长性状的关联性,为深入研究SOCS4基因遗传效应及其功能和调控机制奠定基础.[方法]随机抽取106头威宁黄牛,采集血样后提取总DNA,利用PCR-SSCP检测SNP位点及基因型,并分析SOCS4基因SNP位点与不同个体生长性状的相关性.[结果]在SOCS4基因中共发现两个SNP位点,分别是第1内含子的T+435C和3+UTR的G+13700C,均有3种基因型.在威宁黄牛群体中,T+435C位点是以T为优势等位基因,G+13700C位点则以C为优势等位基因,但只有T+435C位点与威宁黄牛的生长性状存在显著关联(P＜0.05,下同),且在公牛和母牛中呈不同关联性特征.在公牛中,T+435C位点与胸围显著相关,TT型个体胸围均值显著高于CC型个体;而在母牛中,T+435C位点与管围、胸深、腰高和坐骨端宽等生长性状显著关联,且均以CC型个体值高于其他基因型个体.[结论]威宁黄牛SOCS4基因多态性不丰富,且只有内含子区的T+435C位点与其生长性状存在明显关联,即T+435C位点可作为威宁黄牛生长性能的辅助选择分子标记.     

牛PRKAA2基因SNP检测及与其生长和肉质性状的关联分析

【目的】研究秦川牛、南阳牛、鲁西牛和安格斯牛共630头个体的单磷酸腺苷激活的蛋白激酶α2亚基(AMPKα2)基因PRKAA2的多态性,并分析其多态性与200头秦川牛生长及肉质性状的关联性。【方法】利用PCRSSCP、DNA测序、飞行质谱(MALDI-TOF)等技术,分析检测4个品种牛PRKAA2基因6个外显子的单核酸多态性(SNPs);采用生物信息学方法结合统计分析软件,分析秦川牛各多态位点的遗传变异特性与体尺及肉质性状的关联性。【结果】PRKAA2基因第2、4、6、7、8、9外显子中,只检测到第4外显子27GC和60TC 2个突变位点;关联分析表明,27GC位点与秦川牛背膘厚和眼肌面积显著关联(P0.05);60TC位点与秦川牛体斜长极显著相关(P0.01),与体高、尻长、腰角宽、坐骨端宽、眼肌面积显著相关(P0.05)。【结论】PRKAA2基因对秦川牛部分体尺及肉质性状有极显著或显著影响,可作为秦川肉牛新品种早期选择的候选基因和分子辅助标记。     

绵羊Izumo1基因多态性及其与产羔数关联分析

为研究Izumo1基因多态性及与产羔数关系,本研究利用PCR和直接测序法对小尾寒羊和苏尼特羊的Izumo1基因DNA序列进行扩增、测序和BLAST分析,并和本课题组前期绵羊重测序数据进行比对,筛选出Izumo1基因SNP位点。同时,采用Sequenom MassARRAY~?技术进行基因分型,并对其SNP位点的基因型和等位基因频率在各群体中的分布进行研究。结果表明:小尾寒羊Izumo1基因DNA全长序列3 385 bp,苏尼特羊Izumo1基因DNA全长序列3 382 bp;筛选出8个Izumo1基因SNP位点,经过初步筛选,将在小尾寒羊、苏尼特羊、滩羊、萨福克羊、杜泊羊、草原藏羊这6个绵羊品种中位点分布没有差异的SNP位点排除,最终筛选获得的g.54412135AG和g.54412107CA 2个位点。Izumo1基因g.54412135AG位点在多羔和单羔绵羊品种中存在GG、GA和AA 3种基因型,G基因为优势等位基因;g.54412107CA在多羔品种中存在CC和CA 2种基因型,而在单羔品种中存在CC、CA和AA 3种基因型,C基因为优势等位基因,其基因型频率和基因频率在单、多羔绵羊群体间的分布差异均极显著(P0.01);群体遗传学分析得出g.54412135AG多态位点在6个品种中的多态信息含量(PIC)都属于低度多态(PIC0.25);g.54412107CA多态位点在杜泊羊中属于中度多态(0.25PIC0.5),在其他5个品种中都属于低度多态(PIC0.25),然而,关联分析发现Izumo1基因g.54412135AG和g.54412107CA位点的不同基因型与小尾寒羊不同胎次产羔数之间不存在显著关联(P0.05)。