论文部分内容阅读
猪的脂肪代谢是一个十分复杂的过程,涉及多个组织器官的调控,其中甲状腺分泌的甲状腺激素在脂质代谢过程中发挥着重要的作用。目前对甲状腺激素调控的研究主要集中在DNA和mRNA层面,对非编码RNA介导的基因表达调控的研究还比较匮乏。长链非编码RNA (long non-coding RNA, lncRNA)和microRNA都是重要的调控性非编码RNA。约克夏猪是瘦肉型猪,生长快且脂肪率低,而金华猪是脂肪型猪,生长慢但脂肪率高,它们在脂肪代谢方面有很大的差异。本研究以阉割的120日龄的约克夏猪和金华猪为研究对象,比较它们的血清甲状腺激素水平,并利用RNA-seq和small RNA-seq技术筛选出了大量在约克夏猪和金华猪甲状腺组织中差异表达的mRNA、lncRNA和miRNA,并对差异表达的mRNA及差异表达的lncRNA和miRNA的靶基因进行了功能注释和通路分析,利用荧光定量PCR技术对筛选出的部分差异基因进行了验证,此外,还将获得的RNA-seq结果和small RNA-seq结果进行了联合分析。主要研究内容和结果如下:(1)约克夏猪的体重和甲状腺重均显著高于金华猪(P<0.05),但两者比值不变。约克夏猪的血清TT4和TT3水平显著高于金华猪(P<0.05),约克夏猪的血清FT4和FT3水平高于金华猪,但差异不显著(P>0.05)。(2) RNA-seq共获得约351 M (million)高质量的读长为125 nt的双末端原始读段,其中约83%可以比对到猪参考基因组。共鉴定到22,435个已知mRNA在猪甲状腺组织中表达,其中分别有1,492、2,479个在约克夏猪和金华猪特异性表达。并筛选得到1,018个lncRNA,包括760个lincRNA,140个intronic lncRNA和118个anti-sense lncRNA,此外还筛选到1189个新mRNA。本研究鉴定到的lncRNA与猪已知的蛋白编码基因相比,lncRNA具有更短的转录本和ORF长度、更少的外显子和更低的序列保守性。差异表达分析共筛选出492个mRNA的表达量在两个品种间存在显著差异,其中275个mRNA在约克夏猪甲状腺中高表达,217个在金华猪甲状腺中高表达,差异基因被注释到酪氨酸代谢、钙离子信号通路、甲状腺激素合成等生物学通路,主要富集在细胞周期和基于微管的过程。差异表达分析还筛选出48个lncRNA在约克夏猪和金华猪的甲状腺中差异表达,其中10个经cis预测找到13个靶基因。在猪甲状腺组织中共发现了12类可变剪切,其中TSS(转录起始位点)与TTS(转录终止位点)两种剪切类型所占的比例最大,分别约占38%和30%的百分比。(3) small RNA-seq共获得约44 M (million)高质量的单末端原始读段,其中约74%可以比对到猪参考基因组。共鉴定到293条已知的miRNA在猪甲状腺组织中表达,对应266个miRNA的前体,属于146个miRNA家族,此外还鉴定到60个新miRNA.差异表达分析共筛选出18个miRNA在两个品种间差异表达,其中12个miRNA在约克夏猪甲状腺中高表达,6个在金华猪甲状腺中高表达,生物信息学分析显示差异miRNA参与甲状腺激素调控相关通路,如内吞、甲状腺激素合成、钙离子信号通路、过氧化物酶体、溶酶体等。(4)荧光定量PCR检测了18个差异基因在约克夏猪和金华猪甲状腺组织的表达情况,结果表明与测序结果基本一致,说明本实验的分析结果基本可以反映猪甲状腺组织内基因的真实表达情况。(5)对差异表达的mRNA、lncRNA和miRNA进行联合分析,结果得到的网络调控图中包括237个mRNA、18个miRNA和1个lncRNA。