论文部分内容阅读
环状RNA在三十年前就已被报道。最近随着高通量测序技术的发展,研究人员在越来越多的物种中发现环状RNA,但在硬骨鱼中却仍鲜有报道。大黄鱼是我国一种重要的海水养殖鱼类,本研究通过利用大黄鱼基因组模拟得到的数据,对现有的两种挖掘环状RNA的生物信息学流程进行评估,从中选取合适的流程对实验室原有大黄鱼的真实转录组数据进行环状RNA的挖掘,并进行了大黄鱼雌雄鱼性腺环状RNA的测序、挖掘,比较了两者之间环状RNA的表达差异。主要的研究结果如下:1.利用模拟数据对目前常用的两种环状RNA挖掘流程,从不同的环性与线性转录本测序量方面进行评估,发现在环性转录本较低的条件下,其中的一款流程CIRCexplorer检测到的环状RNA的准确率较高。但在环性转录本较高的情况下,另一款流程CIRI检测到的环状RNA的准确率有所升高,挖掘到的环状RNA的数量也较多,但CIRCexplorer的准确率仍高于CIRI,挖掘到的环状RNA的数量仍低于CIRI。2.根据模拟数据的评估结果,利用CIRCexplorer对大黄鱼多组织混合样本的RNA seq数据进行环状RNA挖掘,检测到了975个环状RNA分子,其中65.95%来自编码基因的外显子;随机挑选3条进行验证,其中2条得到了确认。对挖掘到的外显子环状RNA来源基因进行GO以及KEGG富集分析,结果显示,这些来源基因主要富集在结合、翻译起始因子、细胞质内大分子代谢、翻译等通路中。对其进行microRNA靶点预测的结果显示,大黄鱼环状RNA中有丰富的miRNA的结合位点,其中有363个环状RNA具有超过5个潜在的microRNA结合位点,22个具有10个以上的miRNA-430与let-7家族的结合位点。3.进行了大黄鱼精巢和卵巢环状RNA测序,从测序数据中,共检测到环状RNA 6115个,其中2693个在精巢中特异表达,1053个在卵巢中特异表达;余下2369个在精巢和卵巢都有表达的环状RNA中。表达差异分析结果显示,有1065个环状RNA的表达量在精巢相对上调,而有621个表达量在精巢中相对下调。对这些差异表达的环状RNA的编码基因进行GO富集分析,结果表明这些基因主要富集在神经系统发育(GO:0007399 nervous system development)、蛋白质结合(GO:0005515 protein binding)以及谷氨酸受体活动(GO:0004972 NMDA glutamate receptor activity)等GO术语上。KEGG代谢通过分析结果显示,雌雄性腺环状RNA的编码基因在脂类代谢途径(Ubiquitin meditated proteolysis)、孕激素介导的卵母细胞成熟(Progesterone-mediated oocyte maturation)以及泛素介导的蛋白酶解(TNF signaling pathway)等代谢通路中具有明显差异。性腺中检出的环状RNA具有丰富的mi RNA结合位点,网络分析结果显示大黄鱼性腺中的circRNA-miRNA-mRNA可能具有复杂的调控网络。