论文部分内容阅读
水牛分布于我国南方,具有耐粗饲、抗逆性强、奶质优良、役力大以及使用年限长等生物学特性。随着城市化进程以及机械化程度不断提高,我国水牛的用途逐渐以役用为主向乳肉兼用型发展。尤其是营养丰富的水牛奶及其乳制品深受消费者的喜爱,使得水牛奶业具有广阔的市场前景。至今,我国杂交水牛的单产水平已达到15001700 kg,但仍低于河流型水牛品种的单产水平（20002200 kg）。显然,较低的产奶量已严重阻碍了当前水牛奶业的发展。因此,鉴定与水牛产奶性状相关的致因基因对于提高该物种的单产水平至关重要。随着高通量测序技术的迅猛发展,多组学整合分析策略已成畜禽重要经济性状遗传解析的重要手段。本研究拟从基因型数据评估、基因型和转录组数据的整合分析以及候选基因的家族分析与功能验证等角度多层次挖掘和鉴定与水牛产奶性状相关的候选基因,为开展奶水牛分子育种实践奠定基础。本研究共包括四部分内容:试验1:取411头纯种地中海水牛和45头中国杂交水牛基因组数据,探究连锁不平衡模式和评估有效群体大小发现:（1）纯种和杂交水牛群体中常染色体邻近SNP的平均r2值分别为0.13±0.19和0.09±0.13,且标记距离在200 kb时显示出最快的LD衰减速率。（2）LD系数为0.2时,纯种和杂交水牛群体的衰减距离分别为170 kb和50 kb;LD系数为0.3时的衰减距离则分别为60 kb和10 kb。这两个群体的Phase连续性在标记距离小于100 kb时达到最大值(RP,C=0.47),表明这两个群体间的Phase相关性不是很强。（3）纯种和杂交水牛群体中有效群体大小呈下降模式,在前13世代的评估值分别为387和113头。这些研究结果为进一步开展水牛全基因组关联分析和基因组选择研究提供了理论依据。试验2:为探究不同泌乳期的水牛乳腺组织表达谱特征,鉴定与产奶性状相关的候选基因,本试验开展了乳腺组织转录组测序、全基因组关联分析和加权的基因共表达网络分析（WGCNA）等研究发现:（1）基于转录组数据,共鉴定出1420个差异表达的mRNA;（2）基于基因水平的全基因组关联分析,976个基因被纳入到NSGG数据集,筛选出9个与产奶量相关的候选基因(P<10-4);（3）用WGCNA分析,分别从DEGs和NSGG数据集中鉴定出与产奶量相关的天蓝色模块基因544和225个。最后,筛选出12个与产奶量相关的基因（BNIPL,TUBA1C,C2CD4B,DCP1B,MAP3K5,PDCD11,SRGAP1,GDPD5,BARX2,SCARA3,CTU2和RPL27A）作为枢纽基因。这些研究结果为深入了解水牛乳腺组织转录组特征和乳腺发育的生物学功能具有理论指导意义。试验3:基于微管蛋白是微管的主要成分、在许多细胞过程中起着至关重要的作用,本试验从全基因组水平探究了水牛微管蛋白家族及其表达规律发现:（1）共鉴定出33个水牛微管蛋白序列,根据其结构和系统发育特征,可分成五个亚族;（2）同一亚族的微管蛋白序列具有相似的基因结构以及保守的基序组成;（3）水牛微管蛋白家族成员随机分布于10条染色体上,其扩张主要与串联重复事件有关;（4）基于转录组数据和qRT-PCR验证发现,微管蛋白在各种组织中具有不同的表达模式。多数微管蛋白基因在泌乳早期高度表达,可能与乳腺上皮细胞的凋亡有关。这些结果将为进一步验证微管蛋白的功能,鉴定与水牛产奶相关候选基因奠定了理论依据。试验4:基于以上候选基因的筛选结果,本试验用RNA干扰和过表达技术在水牛乳腺上皮细胞（BBMEC）对TUBA1C基因进行功能验证发现,沉默TUBA1C基因可显著影响BBMEC的生长,降低细胞活力,导致凋亡增加,但过表达TUBA1C则对BBMEC的生长影响不显著。此外,干扰或过表达可分别引起11个乳脂代谢和3个乳蛋白代谢通路中相关基因的表达,从而影响BBMEC中甘油三脂的合成和β/κ酪蛋白的分泌。综上,基于系统生物学原理,利用多组学数据和功能验证鉴定水牛产奶性状相关的候选基因,初步阐明了其作用机制,提高了筛选基因的可靠性,为制定水牛分子育种方案、促进其遗传改良提供了有效的分子佐证。