论文部分内容阅读
山羊乳中总脂肪含量,尤其是其中与人类健康密切相关的中短链脂肪酸及不饱和脂肪酸含量均高于牛乳。山羊乳的乳脂特点应与其乳腺乳脂代谢密切相关。然而,在乳脂代谢基因中,至今尚未发现山羊特异性表达或者具有特异功能的基因。据此推测山羊乳脂代谢的特点可能与其在基因不同水平(如转录后、翻译水平)的调控机制有关。MicroRNAs (miRNAs)是一类在转录后水平调控基因表达的小分子RNA,其功能几乎涉及生命活动的各个方面。已证实miRNAs是调控脂肪和肝脏组织中脂肪代谢的重要因子,然而,现有研究缺乏(?)niRNAs调控山羊乳腺乳脂代谢的研究报道。因此,本研究从miRNAs调控乳腺乳脂代谢的角度出发,通过测序技术获得西农萨能奶山羊乳腺组织miRNAs表达谱,分析其序列特征,检测其在乳腺组织不同阶段的表达水平,筛选与山羊乳脂代谢相关的miRNAs,并通过在山羊乳腺上皮细胞(GMEC)中过表达miRNAs,检测细胞中与乳脂代谢相关的指标,证实miRNAs对山羊乳脂合成和组成具有重要影响,同时发现3对可在乳腺组织中协同表达的miRNAs,发现其协同调控山羊乳腺组织乳脂代谢。本研究从miRNAs调控角度阐明山羊乳脂代谢分子调控机制、揭示山羊乳腺泌乳机理,为后期增加乳脂含量及改善乳组成等研究提供重要的理论和试验依据。论文主要结果如下：1.利用Solexa测序技术检测山羊泌乳中期乳腺组织中的小RNAs序列,测序结果经数据库过滤后得到高质量小RNAs序列读数22,084,321个,其中长度为22个核苷酸(nt)的小RNA读数最多,约占所有小RNAs读数的33.4%。将测序所得序列与小RNAs数据库中已知序列进行比对,结果发现山羊乳腺组织中含有5种小RNA,分别为tRNAs、 rRNAs、snoRNAs、snRNAs和miRNAs,其中miRNAs的读数最多(10,534,221),约占乳腺组织小RNAs,总读数的47.7%。对miRNAs序列进行鉴定,共获得1,143个miRNAs,其中931个为保守miRNAs,212个为山羊新的候选miRNAs。将山羊931个保守miRNAs与miRBase (17.0)数据库中已注释的miRNAs序列进行比对,结果发现：(1)山羊和牛的保守niRNAs数目最多,达575个,其中,山羊与牛特有的miRNAs数目为116个；(2)335个miRNAs的391个末端碱基发生变化,其中碱基减少的数目大于增加的数目,且miRNAs3’端的碱基变化数目大于5’端；(3)109个碱基发生替换,包括63个转换和46个颠换。2.利用定量PCR技术检测山羊30个miRNAs在泌乳中期和干奶期乳腺组织中的表达量,获得17个差异表达miRNAs。利用PicTar和TargetScan软件预测miRNAs的靶基因,并通过GO和KEGG软件预测靶基因的功能并进行靶基因富集,结果发现17个差异表达miRNAs的靶基因主要参与乳腺发育、细胞增殖、乳脂合成等代谢过程。从17个差异表达miRNAs及10个功能已知(调控其它动物脂质代谢),miRNAs中,通过催乳素实验获得4个(miR-23a、miR-27a、miR-103和miR-200a)表达量受催乳素调控的miRNAs,推测这4个miRNAs可能与山羊乳腺组织乳脂代谢有关。3.从山羊基因组DNA上扩增出200bp至500bp长的pri-miR-23a、pri-miR-27a、 pri-miR-103和pri-miR-200a序列(皆包括pre-miRNA及侧翼序列),与牛的序列相比,山羊pre-miR-27a3’端第11个碱基由A变为G,山羊其余pre-miRNAs序列与牛相同。以pri-miRNAs为插入片段,成功构建重组腺病毒载体pAd-miRNAs。将腺病毒载体转染HEK293细胞株,病毒载体包装成功后获得Ad-miR-23a、Ad-miR-27a、Ad-miR-103和Ad-miR-200a4种腺病毒,用这4种病毒分别感染GMEC,72h后检测细胞中miRNAs表达水平,结果发现病毒介导miRNAs的表达量依次为对照(感染Ad的GMEC)的3.48、1.94、2.92和2.48倍。4.检测感染Ad-miR-10372h时GMEC中的乳脂合成相关指标,结果显示：miR-103过表达可促进细胞中乳脂积累,与对照(感染Ad的细胞)相比,细胞中甘油三酯的含量增加0.33倍,总脂肪酸含量增加0.16倍,不饱和脂肪酸含量增加0.42倍,其中,不饱和脂肪酸c9-C18：1和c9,t11-C18：2含量大幅上调,较对照相比分别上调0.35倍和1.26倍,同时,脂肪酸和甘油三酯合成及脂滴生成相关基因FASN, DGAT1、ADRP、SLA27A6等表达量升高。在感染病毒0、24、48和72h的GMEC中,miR-103过表达可上调PPARγ、 SREBP-lc和LXRa及它们下游基因mRNAs表达水平,同时下调脂解和氧化相关基因HSL、ATGL、CPT1和ACOX1的表达量。此外,miR-103与PANK3在乳腺组织中的表达量之间具有正相关关系,相关系数达0.891,且在GMEC中进行miR-103过表达和抑制试验时,PANK3的表达量随着miR-103表达量变化而改变。5.检测感染Ad-miR-27a72h时GMEC中的乳脂合成相关指标,结果显示：与对照(感染Ad的细胞)相比,过表达miR-27a能抑制细胞中脂滴形成,下调甘油三酯含量,降低不饱和脂肪酸/饱和脂肪酸的比例,其中,不饱和脂肪酸c9,C18：1和c9,t11-C18：2含量大幅下降,较对照分别下降0.36和0.45倍,而饱和脂肪酸C16：0和C18：0含量升高,较对照分别升高0.19和1.02倍。检测感染病毒0、24、48和72h时细胞中乳脂代谢相关基因的表达量：PPARγ蛋白表达量降低,其下游基因、负责甘油三酯合成的DGAT1的表达量下降,同时,与脂滴形成相关的基因ADRP和TIP47表达量升高,与脂解、氧化相关的基因HSL、ATGL和ACOX1表达量升高。6.检测感染Ad-miR-23a72h时GMEC中的乳脂合成相关指标,结果显示：与对照相比,过表达miR-23a的细胞中脂滴含量增加1.83倍,甘油三酯含量增加2.21倍,同时,乳脂合成相关基因FASN、DGAT1和TIP47的表达量升高,依次为对照的1.82、1.89、2.93倍,而脂解相关基因ATGL的表达量降低,仅为对照的0.6倍。此外,过表达miR-23a还可增加固醇转运相关基因ABCA1和ABCG2的表达水平。检测感染Ad-miR-200a72h时GMEC中乳脂代谢相关基因的表达量,结果发现miR-200a对乳脂合成各阶段相关基因的表达皆有影响,其可降低FASN和ADRP的表达量,提高SCD、 DGTA1和HSL的表达量,同时还可调控PPARγ和SREBP-1c的表达水平。7.利用Pearson相关系数分析miR-23a、miR-27a、miR-103和miR-200a在泌乳中期山羊乳腺组织中表达量之间的相关性,结果发现miR-23a与miR-27a、miR-27a与miR-200a、miR-200a与miR-103三对miRNAs之间的表达量具有正相关关系,它们的相关系数依次为0.79、0.57和0.716。此外,在这三对miRNAs中,miR-23a与miR-27a、 miR-27a与miR-200a之间在过表达时均表现为相互促进作用,而miR-103与miR-200a在过表达时表现为相互抑制作用。