脂肪组织不仅在个体发育过程中发挥必要作用,并且在牛肉品质的研究中也扮演着重要角色。培育生产优质牛肉的肉牛品种已成为目前分子育种的重要目标。牛肉中脂肪含量决定了牛肉的品质。但由于脂肪组织形成和沉积的分子作用机制复杂,在牛上对脂肪沉积相关基因的研究相对较少,有待挖掘;并且基因的可变剪接丰富了其表达产物的种类,加大了基因作用机制的复杂程度,同时也为分子机制的研究提供了新思路,成为当今研究的热点。为此,本研究对不同发育阶段和不同性别秦川牛的皮下脂肪组织进行转录组测序和基因差异表达分析,并在组织、细胞和分子水平上对差异表达基因逐步筛选得到脂肪沉积相关基因,最后利用电子克隆、基因组织表达谱构建、脂肪细胞分化基因时序表达分析、细胞共定位及亚细胞定位分析等手段进一步研究了可变剪接的组织特异性以及可变剪接对脂肪沉积相关基因功能的影响。主要研究结果如下:1.牛皮下脂肪组织转录组测序及基因表达分析本研究利用二代高通量测序分别对胎牛、成年公牛、成年母牛和成年阉牛皮下脂肪组织进行转录组测序。每个样品获得至少2,600万clean reads,其中82.58%~84.41%的clean reads成功的比对到参考基因组上。我们共检测了12,233个基因在四类皮下脂肪组织内的表达,其中109个基因在四类皮下脂肪组织中共同高表达(≥500 RPKM)。共同高表达的基因中除了含有63个编码核糖体蛋白的基因外,还包含多个与脂肪代谢相关的基因(APOE、GAPDH、FABP4、SPARC和ADIRF等),说明共同高表达的基因可以为筛选脂肪沉积相关基因提供材料。不同性别成年牛皮下脂肪组织基因表达高度相关(r≥0.91),胎牛和各成年牛皮下脂肪组织基因表达相关程度较低(r≤0.64);并且通过对四类皮下脂肪组织两两比较差异基因数目统计分析,发现本研究中年龄相对性别对皮下脂肪组织基因的表达影响更大。此外,我们还在四类牛脂肪组织中检测到大量新转录本,这为将来完善牛基因组注释及非编码RNA的研究提供资源。2.年龄和性别对脂肪组织基因表达差异的影响通过将成年公牛、成年母牛和成年阉牛皮下脂肪组织作为三个重复与胎牛皮下脂肪组织的基因表达进行比较,共检测到2,703个差异表达基因(差异表达倍数>2,P<0.05)。成年牛和胎牛皮下脂肪组织差异基因GO(Gene ontology,基因功能注释)富集分析显示,差异基因显著富集(P<0.05)在与组织发育相关的GO条目;脂肪代谢相关的GO条目没有显著富集,但具有较低的P值。差异基因显著富集在脂肪代谢关键信号通路PPAR pathway,其中促进脂肪沉积的基因在成年牛皮下脂肪组织中显著上调;另外,14个脂肪细胞成脂分化标志基因中有12个基因在成年牛皮下脂肪组织内相对胎牛皮下脂肪组织高表达,说明成年牛皮下脂肪组织相对胎牛皮下脂肪组织进行着更为活跃的脂肪沉积。对成年公牛、成年母牛和成年阉牛皮下脂肪组织两两比较差异基因进行功能和参与的信号通路富集分析发现,差异基因在与甾类激素刺激、脂肪代谢和细胞因子相关的go条目显著富集,在和细胞因子相关的信号通路共同显著富集,这初步证明性别对牛皮下脂肪组织的脂肪代谢有影响,并且可能和细胞因子的分泌相关。3.脂肪沉积相关基因的筛选基于以上对转录组测序数据的分析结果,选取11个在成年牛皮下脂肪组织内均为胎牛皮下脂肪组织内表达量10倍以上的基因和20个在四类牛皮下脂肪组织内共同高表达的基因,作为后续逐步筛选与脂肪沉积相关基因的初始材料。通过构建31个候选基因的组织表达谱,发现其中12个基因(adipor2、cidec、ghr、lipe、s100b、scd、thrsp、tusc5、anxa2、cst3、sparc和vim)在脂肪组织内相对其它组织高表达。通过罗格列酮对牛原代脂肪细胞进行脂肪沉积诱导,检测到其中7个基因(ghr、thrsp、scd、lipe、tusc5、cidec和cst3)在脂肪沉积过程中表达呈上调趋势。利用腺病毒在牛原代脂肪细胞中超表达pparg2,检测到其中3个基因(cidec、tusc5和cst3)表达量显著上调。最终,本实验筛选出3个受pparg2调控且与脂肪沉积相关的基因,同时补充了牛pparg2基因的调控网络。4.脂肪组织内表达基因可变剪接特征分析在以上转录组测序数据中共发现4,753个基因存在可变剪接现象,占总检测基因数目的38.85%。但只有1,319(28%)个基因在四类牛脂肪组织中均出现可变剪接现象,并且共同出现可变剪接的基因只有17%的可变剪接事件在四类牛脂肪组织内同时存在,说明可变剪接在脂肪组织的不同发育时期或不同性别的个体内存在较大差异。本研究选择牛nfix基因对可变剪接的组织特异性进行了进一步验证。通过电子克隆及pcr方法成功预测并验证牛nfix基因存在5种不同转录本。组织表达谱分析显示不同转录本在牛不同发育时期和不同组织内表达出现多样化的差异。以上分析说明同一基因的不同亚型可能会具有不同功能,因此有必要进一步分析可变剪接对脂肪沉积相关基因表达和分布的影响。5.脂肪沉积相关基因可变剪接模式及细胞定位分析本实验分析了三个脂肪沉积相关基因(cebpa、cidec和tusc5)的可变剪接模式,并通过细胞定位分析或细胞共定位分析在239t或3t3-l1细胞中检测可变剪接对以上三个基因在细胞内的表达分布影响。结果显示:cebpa的两种蛋白亚型均分布在细胞核内;cidec的两种蛋白亚型均分布在细胞质内,并且在细胞质内的分布未发现差异;tusc5的两种蛋白亚型均分布在细胞质内,但在细胞质内的分布出现明显差异。对tusc5产生的两种蛋白亚型氨基酸序列进一步分析显示,可变剪接可能通过改变TUSC5的内质网滞留信号影响其在细胞质内的分布。6.可变剪接对脂肪沉积相关基因TUSC5作用的影响对TUSC5不同亚型(TUSC5a和TUSC5b)的组织表达谱和在牛原代脂肪细胞诱导分化不同阶段表达变化趋势分析,均发现TUSC5a和TUSC5b存在一定的差异;并且利用细胞定位载体检测TUSC5a和TUSC5b在293T细胞内的时序表达,发现TUSC5a相对TUSC5b出现明显的迟缓表达。但亚细胞定位分析显示TUSC5a和TUSC5b都在内质网上有分布,CIDEC不在内质网上分布;并且TUSC5两种亚型与CIDEC细胞共定位分析显示都不和CIDEC在细胞质内发生互作。综上,可变剪接会影响TUSC5在细胞质内的表达,但对TUSC5与CIDEC之间的互作关系没有明显影响。本研究对秦川牛皮下脂肪组织进行转录组测序和脂肪沉积相关基因筛选,并将脂肪沉积相关基因的分析深入到可变剪接水平,为牛脂肪组织发育分子机制的研究提供了丰富的资源和新思路。