新疆褐牛乳中体细胞数与产奶性状的影响因素分析

【目的】分析影响新疆褐牛乳中体细胞数（SCC）与产奶性状的非遗传因素,为有效指导并提高新疆褐牛的生产性能提供科学依据。【方法】构建协方差分析模型,利用SAS 9.2软件分析不同场、胎次、泌乳月及两两间互作效应等非遗传因素对262头新疆褐牛SCC与产奶性状的影响,并阐述各性状间的相互影响关系;采用Excel 2 013对新疆褐牛1 510条DHI数据进行整理,利用相邻两月SCC差值等级划分方法,绘制了不同牛场SCC变化趋势及其与产奶性状间的变化规律。【结果】场和胎次对新疆褐牛SCC及产奶性状均有极显著影响（P＜0.01）;泌乳月效应对产奶量和乳中SCC影响极显著（P＜0.01）,对乳糖率有显著影响（P＜0.05）,对其它乳成分无显著影响（P＞0.05）;不同场、胎次、泌乳月两两间互作效应对新疆褐牛SCC均无显著影响（P＞0.05）,而场和胎次互作效应对产奶量和乳成分有极显著影响（P＜0.01）,场和泌乳月互作对乳脂率和乳中干物质有极显著影响（P＜0.01）,对乳糖率有显著影响（P＜0.05）,对其它性状无显著影响（P＞0.05）;胎次和泌乳月互作效应对新疆褐牛各产奶性状均无显著影响（P＞0.05）。随着新疆褐牛乳中SCC增加,乳中干物质含量增加,而产奶量、乳脂率、乳蛋白率和乳糖率下降;绘制泌乳期SCC差值与产奶量和乳糖率变化曲线发现,随着SCC差值升高,产奶量及乳糖率均降低,表明使用SCC差值等级法能够有效评价奶牛产奶量变化并及时判定奶牛乳房炎的发病趋势。通过比较新疆褐牛两个牛场的产奶量、乳成分以及SCC控制情况,发现二号牛场各性状均优于一号牛场,其产奶量较高、SCC较低,且不同胎次水平的乳脂率、乳糖率及乳中干物质均较高,体现较好的生产管理水平。【结论】针对新疆褐牛非遗传因素的研究表明,为有效提高和改善该乳肉兼用牛品种的产奶性能和奶品质,应着力改善牛场管理水平、优化环境卫生及规范挤奶设备等;一号牛场应加强生产管理水平和DHI监测记录,通过控制新疆褐牛SCC来降低和避免泌乳牛的乳房炎隐患。     

基于全基因组重测序的山羊产羔数性状关键调控基因的筛选

目的 对产羔数不同的山羊进行全基因组重测序分析，挖掘参与调控川中黑山羊产羔数性状关键调控基因，为解析山羊产羔数性状遗传机制及分子遗传改良提供理论依据。方法 选择6只产4—6羔的川中黑山羊为高繁组（high fecundity, HF）和6只产单羔的川中黑山羊为低繁组（low fecundity, LF），采集颈静脉血液样本提取基因组DNA，构建350 bp双末端测序文库，利用Illumina HiSeq PE150平台对12个文库进行全基因组重测序。测序产出的净数据经BWA软件比对至山羊参考基因组ARS1，所获得的高质量SNPs通过两种全基因组扫描分析方法（Fst、Hp）的综合分析确定候选区域，候选区域的注释基因分别利用g:Profiler和KOBAS在线数据库进行GO分析与KEGG通路分析，以筛选调节川中黑山羊产羔数性状候选基因。为了进一步鉴定调节山羊产羔数目的关键遗传标记，根据基因组重测序变异分析，对繁殖候选基因的同义与非同义SNPs进行定位筛选，后续将12个山羊样本的扩增产物进行Sanger测序以验证重测序结果。结果 12只山羊共获得431.50 Gb 净数据，经变异检测与注释发现，HF组山羊共发现7 771 417个单核苷酸多态性(single nucleotide polymorphism, SNPs)，LF组山羊检测到8 935 907个SNPs，且LF组各类SNPs 均多于HF组。设置同时达到top 5%最大ZFst值和top 5%最小ZHp值的窗口为候选区域，在低杂合性、高遗传分化的区域共注释130个强选择信号，其中HF组、LF组以及共享窗口的注释基因分别为84、59和13个，经GO富集与KEGG通路分析发现，19个候选基因参与川中黑山羊的繁殖、繁殖过程和胚胎发育等调控，包括11个HF组特异性候选基因（ADCY10、DRD1、HS6ST1、IGFBP7、MSX2、NOG、NPHP4、PAPPA、PRLHR、TDRP和XYLT1），5个LF组特异性候选基因（ANXA5、IGF1、EDNRA、FANCL 和TAC1）和3个HF组与LF组共享窗口基因（AKR1B3、HDAC4和OPRM1）。同时，大多数GO分析，如G蛋白偶联受体活性、激素反应和神经肽信号通路等，都包含这19个候选基因。此外，14个HF候选基因有9个显著富集在代谢途径、神经活性配体-受体相互作用、糖胺聚糖-硫酸乙酰肝素/肝素的生物合成、钙离子信号通路、cAMP信号通路和叶酸生物合成等KEGG通路中（P<0.05）。19个繁殖候选基因中共有2个同义突变（MSX2 G771T，ADCY10 A4662G）与2个非同义突变（PRLHR G529A，DRD1 A281T），且仅定位于HF候选基因中。Sanger测序发现，MSX2、PRLHR和DRD1突变位点均可检测到多态性，与基因组重测序结果一致，其中PRLHR G529A多态性导致第177位丙氨酸突变为苏氨酸，DRD1 A281T多态性导致翻译提前终止。结论 本研究共发现11个HF组特异性候选基因，推测是川中黑山羊多羔性状的关键调控基因，PRLHR外显子G529A与DRD1外显子A281T突变可能是调控山羊多羔性状的关键遗传标记，在改良山羊繁殖性能方面具有较大的应用价值。     

新疆褐牛基因外显子多态性及其与SCS和泌乳性状关联性分析

【目的】探索ERAP2,Cwc27,PPFIBP2,THADA,ZNF804B和AMHR2基因多态性与新疆褐牛体细胞评分及泌乳性状的相关性,旨在寻找与新疆褐牛体细胞评分及泌乳性状相关的分子标记。【方法】以新疆乌鲁木齐种牛场、新疆天山畜牧生物工程股份有限公司共169头新疆褐牛母牛为试验对象,以2008—2015年间的生产性能测定记录(包括乳脂率、乳蛋白率、总固体含量、乳糖率和体细胞计数等)和305 d产奶量为试验数据,基于课题组前期15头中国荷斯坦牛DNA混池重测序的结果筛选外显子区域10个SNPs,采用基质辅助激光电离飞行时间质谱法(sequenom mass Array genotype)分型技术验证10个SNP位点的遗传多态性,运用SAS8.1软件GLM过程对这10个SNP突变与新疆褐牛体细胞评分(SCS)、305 d产奶量、乳脂率、乳蛋白率、乳糖率和总固体的关联程度进行统计分析。【结果】10个SNP位点均具有多态性,χ~2检验表明群体中9个SNP位点符合Hardy-Weinberg平衡状态。关联分析结果表明ERAP2基因两个位点(T98741711C和G98736141A)均与新疆褐牛305 d产奶量达到极显著关联(P0.001),其中TT型和GG型个体的产奶量最高;PPFIBP2(C45667492G)位点与新疆褐牛乳糖率达到显著关联(P0.005),其中GG型个体的乳糖率最高;Cwc27(T14533269A)位点与新疆褐牛乳脂率和总固体含量达到显著关联(P0.005),与SCS达到极显著关联(P0.001),其中TT型个体的乳脂率和总固体含量均最高,AA型个体的SCS极显著高于AT和TT型;AMHR2(C26758055G)位点与新疆褐牛乳蛋白率和总固体含量达到极显著关联(P0.001),其中GG型个体显著高于CC和GC型个体。连锁不平衡分析与单倍型构建的结果发现10个SNPs构建了两个单倍型模块,其中SNP1和SNP2位点之间处于连锁不平衡状态(r20.3),SNP4和SNP10位点之间处于强连锁不平衡状态(r20.6),随后分析了单倍型与SCS和泌乳性状之间的相关性,发现单倍型与新疆褐牛SCS和泌乳性状无显著关联(P0.001)。【结论】初步发现了ERAP2,Cwc27,PPFIBP2和AMHR2与新疆褐牛SCS和泌乳性能相关,结果可为新疆褐牛产奶性状的分子标记辅助选育提供参考和依据。     

利用加权基因共表达网络分析牛口蹄疫病毒感染通路变化

从基因表达综合数据库下载GSE83514数据集,利用加权基因共表达网络分析(WGCNA)筛选与牛口蹄疫病毒感染的相关模块及枢纽基因,并对整个模块内的基因进行GO功能注释和KEGG通路分析。结果显示,通过构建WGCNA共表达网络,确定与牛口蹄疫病毒感染显著负相关的Darkred模块和显著正相关的Green模块为枢纽模块,它们的枢纽基因分别为MFSD4和RHOH。这两个模块基因共富集到20个生物学过程及14条KEGG信号通路,并且两个模块同时富集到DNA复制通路。本研究为探究牛口蹄疫病毒感染通路变化提供了生物信息学依据,筛选的枢纽基因MFSD4和RHOH有望成为抑制牛口蹄疫病毒感染的治疗靶标。     

基于比较转录组学分析揭示中华蜜蜂及意大利蜜蜂幼虫的球囊菌抗性差异机制

为探究中华蜜蜂(Apis cerana cerana,以下简称中蜂)与意大利蜜蜂(Apis mellifera ligustica,以下简称意蜂)幼虫的球囊菌(Ascospahaera apis)抗性差异产生的原因,利用Venn分析、GO分类分析、KEGG代谢通路分析以及Ka/Ks分析对二者的转录组进行比较研究。Venn分析结果显示,中蜂与意蜂的同源基因有17 656个,归属于8 111个基因家族,二者的特有基因分别有1 158和241个,分别归属于468和86个基因家族。GO分类结果显示,中蜂和意蜂的特有基因分别富集在38和28个GO条目。KEGG代谢通路富集分析结果显示,中蜂和意蜂的特有基因分别富集于96和21个代谢通路;进一步分析发现中蜂的特有基因富集在内吞作用、溶酶体、泛素介导的蛋白水解和黑化作用等4个细胞免疫通路,以及MAPK信号通路和Toll-like受体信号通路等2个体液免疫通路,而意蜂仅有1个特有基因富集在MAPK信号通路。Ka/Ks分析结果显示受到强烈正向选择、弱正向选择、中性选择和负向选择的单拷贝同源基因分别有37、281、2 630和3 269个。对受到强烈正向选择的基因进行GO分类和KEGG代谢通路富集分析,GO结果显示中蜂和意蜂的单拷贝同源基因富集的GO条目类型相同;KEGG结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一,中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。     

基于RNA-seq技术对不同品种猪背最长肌差异表达基因的筛选与注释

【目的】筛选皖南花猪和大约克猪背最长肌组织差异表达基因。【方法】采集皖南花猪和大约克猪背最长肌(各3头),测定其肉质性状指标,并提取其RNA,采用高通量转录组测序(RNA-seq)技术进行测序,对所获序列进行GO和Pathway显著性富集分析。【结果】测序后皖南花猪3个样本得到的总序列数分别为65 584 126,64 327 738和59 244 502,其中有效序列达到94.4%,94.3%和94.2%;大约克猪3个样本得到的总序列数分别为57 969 134,68 254 168和72 298 142,其中有效序列达到93.8%,93.7%和93.8%。测序饱和度良好,证明测序数据真实可靠。差异表达分析结果显示,共筛选到347个差异表达基因,其中包含94个上调基因,253个下调基因。GO功能显著性富集和Pathway显著性富集分析发现,差异表达基因富集在与糖酵解代谢和肌肉发育相关的生物学过程中以及与脂肪酸代谢和生长性状、胴体性状相关的信号通路上。【结论】获得了猪背最长肌组织基因表达谱,筛选出了一些与肉质性状和生长性状有关的候选基因。     

育成期高能饮食对蛋鸡肝脏转录组的影响

[目的]明确相关基因在育成期蛋鸡肝脏脂质代谢调控过程中的作用机理,为改善蛋鸡的生产性能打下基础.[方法]分别选取高能量(12.14 MJ/kg)和低能量(10.88 MJ/kg)饲料饲喂育成期蛋鸡,利用Illumina NextSeq 500平台对蛋鸡肝脏进行转录组测序,通过HTSeq和DESeq等生物信息学方法筛选出差异表达基因,并以超几何分布对差异表达基因进行GO功能富集分析和KEGG信号通路富集分析.[结果]从低能组和高能组蛋鸡肝脏样品RNA-seq测序获得的原始序列均超过8000万条,且能比对上基因组和基因的序列均在80.00%以上.其中,基因序列占87.32%~89.78%,基因间序列占10.22%~12.68%;外显子序列在基因序列中的所占比例很高,为78.10%~82.88%,内含子序列所占比例在17.12%~21.90%.相对于低能组,高能组蛋鸡肝脏中有82个基因上调表达、106个基因下调表达.差异表达基因GO功能富集分析发现有314个显著富集的GO功能(P<0.05,下同),其中,241个富集在生物过程(BP),18个富集在细胞组分(CC),55个富集在分子功能(MF).脂质代谢过程和细胞脂质代谢过程富集的基因分别有16和14个,占差异表达基因总数的8.51%和7.45%;脂质代谢过程较细胞脂质代谢过程多出的基因是PLIN1和FGF19.差异表达基因KEGG信号通路富集分析发现有4条显著富集的KEGG信号通路,且均与蛋鸡肝脏的脂质代谢相关.[结论]育成期对蛋鸡进行高能饮食饲喂,主要是通过上调或下调脂质代谢相关基因表达而影响肝脏脂质代谢,最终实现蛋鸡生产性能调控.     

基于BSA-seq技术挖掘糙皮侧耳抗螨候选基因

根据前期不同糙皮侧耳菌株抗螨性差异调查结果,构建了糙皮侧耳孢子单核体群体,群体内479个菌株抗螨性呈正态分布。利用BSA-seq技术对高抗和高感混池进行SNP和InDel位点的欧氏距离关联分析,将抗螨性候选基因定位到一条染色体上,2个相邻候选区域总长度为1.75 Mb,内有基因605个,其中非同义突变基因353个,移码突变基因89个。经功能注释以及GO和KEGG通路富集等分析,发现26个候选基因参与了信号传导、防御过程和次级代谢相关通路,推测这些候选基因可能与糙皮侧耳抗螨性相关。     

卵形鲳鲹卵巢不同发育时期的转录组测序分析

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的cDNA文库,采用Illumina HiSeqTM 2500进行转录组测序,经过滤、质量控制及拼接组装后获得的Unigenes在七大数据库（Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO）中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行功能注释及信号通路富集分析,并采用MISA和GATK3进行SSR鉴定及SNP分析。【结果】卵形鲳鲹卵巢组织转录组测序获得的325156432条Raw reads,经过滤筛选得到317206752条Clean reads,拼接组装后得到59554条Unigenes;69.65%的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等七大数据库中注释成功,其中有24599条Unigenes被注释到GO数据库,15997条Unigenes被注释到KEGG数据库。在卵形鲳鲹卵巢组织的2个发育时期共鉴定获得56115个基因,经差异表达分析后获得17737个差异基因,其中8169个基因在卵巢Ш期上调表达、9568个基因在卵巢Ш期下调表达。GO功能注释分析发现,卵形鲳鲹卵巢差异表达基因主要注释在细胞过程、氮化合物代谢过程、初级代谢过程、核、核部分、离子结合及水解酶活性等条目上;而KEGG信号通路富集分析结果显示,17737个差异表达基因显著富集在318条代谢途径上,其中前20条KEGG信号通路包括2-氧代羧酸代谢、PI3K-Akt信号通路、甲状腺激素信号通路、磷脂酶D信号通路、Fc εRI信号通路和细胞周期等。卵形鲳鲹卵巢转录组（59554条Unigenes）中共存在30133个SSRs和82490个SNPs。【结论】GnRHR、FSHR、FSHβ、CYP11A、SIRT3和PEG3等差异表达基因及PI3K-Akt信号通路和VEGF信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

意大利蜜蜂4、5和6日龄幼虫肠道发育过程中差异表达基因的趋势分析

为探讨意大利蜜蜂(简称意蜂)幼虫肠道发育过程中的基因表达谱和基因差异表达信息,研究基于前期已获得的高质量的意蜂幼虫肠道转录组数据,利用STEM软件对意蜂4、5和6日龄幼虫肠道的差异表达基因(DEGs)进行趋势分析,进一步通过相关生物学软件对显著上调趋势包含的DEGs进行GO与KEGG代谢通路富集分析及表达量聚类分析。趋势分析结果显示,5 920个DEGs共聚类在8个基因表达模式,其中935个DEGs聚类在1个显著上调趋势。GO富集分析结果显示,显著上调趋势包含DEGs富集于39个GO条目,富集基因数最多的为结合,其次是细胞进程和单组织进程。KEGG代谢通路富集分析结果显示,显著上调趋势包含DEGs富集于94个代谢通路,富集基因数最多的是嘌呤代谢、内质网蛋白质加工和内吞作用。表达量聚类分析结果显示,随着意蜂幼虫肠道发育时间的延长,代谢和免疫相关通路富集基因的表达水平逐渐增强。研究提供了意蜂幼虫肠道发育过程的基因表达谱和基因差异表达信息,揭示了基因表达规律,为深入研究意蜂幼虫肠道发育的分子机制打下了一定基础。     

松乳菇子实体两个发育时期的转录组分析

以松乳菇子实体为研究材料,采用de novo测序对松乳菇子实体进行测序与生物信息学分析,探究松乳菇子实体生长发育过程中幼菇期(LDG-1)和成熟期(LDG-2)基因表达的差异,并探寻影响松乳菇生长发育相关的主要基因和代谢通路。测序结果显示,共获得7.44G(LDG-1)和7.21G(LDG-2)的分析数据(clean reads)。KOG功能注释结果显示,松乳菇在幼菇期(LDG-1)和成熟期(LDG-2)具有全面而复杂的基因功能类别。GO富集结果提示,松乳菇在不同生长阶段(LDG-1和LDG-2)的生物代谢功能存在一定的差异。基因表达水平与聚类分析显示,ATP酶、多功能伴侣蛋白等家族基因为LDG-1和LDG-2主要高表达基因,在能量代谢与参与调控细胞周期中起重要作用。分析差异表达基因发现,相较于LDG-1期,LDG-2期有16 789个差异表达基因,其中7 635个基因上调,9 154个基因下调。差异表达基因KEGG富集通路分析结果显示,氧化磷酸化途径是松乳菇幼菇期生长发育过程中能量代谢的关键途径。进一步分析显示,MAPK信号通路是影响松乳菇子实体发育的关键信号通路,在促进松乳菇子实体进行有性生殖和调控菌丝体的极性生长中起重要作用。     

藏黄牛与宣汉黄牛心脏miRNA表达谱比较

【目的】miRNA作为一类非编码RNA,广泛参与机体多种生命活动,研究旨在挖掘miRNA在藏黄牛和宣汉黄牛心脏组织中的差异miRNA,为进一步研究藏黄牛低氧适应的分子机制提供基础数据。【方法】随机选取健康的藏黄牛和宣汉黄牛各3头,迅速采集其心脏组织,使用Trizol法提取RNA,琼脂糖凝胶电泳切胶选取18—30nt的片段,连接3'端和5'端,反转录后进行扩增,凝胶电泳切胶纯化后建立藏黄牛和黄牛各3个文库,利用Illumina HiSeq4000测序平台进行高通量测序;将测序得到的序列进行过滤,比对GenBank和Rfam数据库,筛选出藏黄牛和宣汉黄牛的差异miRNA,进行功能注释及信号通路富集分析;随机选择8个miRNAs,利用实时荧光定量PCR检测其在心脏组织的表达量,以验证测序数据的准确性。【结果】藏黄牛和宣汉黄牛心脏组织的高质量读值的序列分别为17 463 446条和13 662 812条,干净读值为16 552 296条和12 055 304条,且在藏黄牛和宣汉黄牛中高质量核酸序列长度富集最多的均是21nt,分别为37.5%和32.1%;且共筛选出219个差异miRNAs,其中48个显著上调,171个显著下调;GO功能注释得到差异miRNA靶基因分子功能中显著富集的条目有22条,如,GO:0005488（结合）、GO:0005515（蛋白质结合）和GO:0043167（离子结合）;细胞组分中显著富集的条目有20条,如,GO:0005623（细胞）、GO:0044464（细胞组分）和GO:0005622（细胞内）;生物过程中显著富集的条目有13条,如,GO:0035556（细胞内信号转导）、GO:0032774（RNA生物合成过程）和GO:0006351（转录,DNA模板化）,KEGG信号通路分析得到差异表达miRNAs靶基因显著富集到胰岛素信号通路（139个靶基因）、mTOR信号通路（38个靶基因）和HIF-1 信号通路（92个靶基因）等232个信号通路中,其中有12个靶基因共同作用于这3个信号通路,说明miRNAs可能通过这3个信号通路,共同参与藏黄牛低氧适应性的调控;随机选取8个miRNAs进行荧光定量验证,其表达趋势与测序结果表达趋势基本一致,表明高通量测序数据可信度较高。【结论】得到了miRNA在藏黄牛与宣汉黄牛心脏组织中的表达谱,为揭示藏黄牛低氧适应性的调控机制研究奠定了基础。     

中国荷斯坦奶牛PPARGC1A基因内含子9的SNP与产奶性状的相关分析

采用PCR产物直接测序的方法,对213头中国荷斯坦奶牛的PPARGC1A基因进行单核苷酸多态性分析,并研究不同基因型与产奶性状的相关关系。结果表明,PPARGC1A基因内含子9存在C>T SNP位点,可分为有CC、CT和TT这3种基因型,基因型频率分别为0.3568、0.6291和0.0141;该位点突变偏离Hardy-Weinberg平衡状态且中度多态;CC基因型与中国荷斯坦奶牛的乳脂肪率显著相关。PPARGC1A基因可以作为研究中国荷斯坦奶牛泌乳性状的候选基因。     

新疆褐牛瘦素基因在不同部位中表达水平与肉质指标的关联性

[目的]肉品质与物种的遗传因素具有一定的关联性.新疆褐牛瘦素基因在不同部位中表达水平与肉质指标的关联性,为牛肉属性差异的品质标识提供参考.[方法]选取新疆褐牛为研究对象,采用荧光定量聚合酶链式反应(qRT-PCR)技术测定新疆褐牛皮下、肌间、肾周、心周4个部位脂肪组织中的瘦素基因(lep-tin)的表达水平,检测活重、...     

A post-GWAS replication study confirming the association of C1<?Pub Caret?>4H8orf33 gene with milk production traits in dairy cattle

Genome-wide association studies with an Illumina Bovine50K chip have detected 105 SNPs associated with one or multiple milk production traits in the Chinese Holstein population. Of these, 38 significant SNPs detected with high confidence by both L1-TDT and MMRA methods were selected to further mine potential key genes affecting milk yield and milk composition. By blasting the flanking sequences of these 38 SNPs with the bovine genome sequence combined with comparative genomics analysis, 26 genes were found to contain or be near to such SNPs. Among them, the C14H8orf33 gene is merely 87 bp away from the significant SNP, Hapmap30383-BTC-005848. Hence, we report herein genotype-phenotype associations to further validate the genetic effects of the C14H8orf33 gene. By pooled DNA sequencing of 14 unrelated Holstein sires, a total of 18 with seven novel SNPs were identified. Among them, nine SNPs were in the 5′ regulatory region, one in exon 6 and the other in the 3′ UTR and 3′ regulatory region. A total of nine of these identified SNPs were successfully genotyped and analyzed by mass spectrometry for association with five milk production traits in an independent resource population. The results showed that these SNPs were statistically significant for more than two traits [P<(0.0001-0.0267)]. In addition, mRNA expression analyses revealed that C14H8orf33 was ubiquitous in eight different tissues, with a relatively higher expression level in the mammary gland than in other tissues. These findings, therefore, provide strong evidence for association of C14H8orf33 variants with milk yield and milk composition traits and may be applied in Chinese Holstein breeding programs.     

三河牛目标性状边际效益研究

【目的】三河牛是1983年验收的中国自主培育的乳肉兼用牛品种,因起源于内蒙古呼伦贝尔市三河地区而得名,具有耐粗饲、宜牧、抗寒、适应性强的特点。研究旨在估计三河牛目标性状的边际效益,为制定三河牛育种规划提供基础。【方法】依据内蒙古谢尔塔拉种牛场三河牛中心产区2011年的生产和育种数据,参照奶牛和肉牛营养需要获得三河牛的生产性能、市场经济、营养学方面的参数,根据三河牛生产和育种体系,确定三河牛育种目标性状,以效益方程法建立三河牛的生物经济利润模型,计算三河牛乳用性状、肉用性状、功能性状的边际效益,并对影响经济效益的各个因素进行了敏感性分析。【结果】（1）确定了三河牛10个育种目标性状,包括3个乳用性状（产奶量,乳蛋白率和乳脂率）,3个肉用性状（初生重,育肥期日增重和屠宰率）和4个功能性状（首次产犊日龄,产犊间隔,体细胞数和母牛使用年限）：（2）牛群规模固定时,三河牛的总效益为4 808.26元。对收益进行分类,产奶和淘汰母牛的收益分别占总收益的67.08%和7.8%,育肥牛的收益占总收益的20.58%。成本方面,可变成本占总成本的96.89%,其中饲养费用是影响可变成本的主要因素,青年母牛、泌乳母牛和育肥公牛的饲养成本分别占总成本的45.69%、20.46%和13.15%;（3）以乳成分和质量（乳脂率、乳蛋白率、体细胞数）定价时,三河牛10个目标性状的边际效益为：产奶量（kg）2.26元、乳脂率（%）40.48元、乳蛋白率（%）61.39元、体细胞数-231.52元、初生重（kg）6.19元、育肥期日增重（g•d-1）1.39元、屠宰率（%）76.70元、使用寿命（天）0.25元、首次产犊日龄（天）-1.62元、产犊间隔（天）-2.71元;（4）乳用性状、肉用性状、功能性状的相对经济权重比例为62﹕14﹕22,近似于3﹕1﹕1的比例。（5）牛奶价格提高10%,产奶量和产犊间隔的边际效益分别增加0.29元和0.17元;牛肉价格提高10%,初生重、育肥期日增重、屠宰率、产犊间隔的边际效益分别增加0.81元、0.15元、7.86元,使用年限的边际效益降低0.03元;饲料价格提高10%,产奶量、乳脂率、乳蛋白率、初生重、首次产犊日龄、育肥期日增重的边际效益分别降低0.05元、0.58元、0.34元、0.21元、0.16元、0.01元,产犊间隔和生产寿命的边际效益分别增加0.88元、0.06元。【结论】市场价格的变动对三河牛乳用性状和肉用性状的边际效益影响较大。研究构建的生物经济利润模型可为在中国实施兼用牛选育提供参考,为制定三河牛育种体系的育种目标提供依据。