技术服务描述

circRNA测序及分析报告

生信部

2020年12月29日

1. CircRNA背景及分析流程简介

1.1. 背景简介

  环形RNA是一类在真核生物中广泛存在的具有特殊环状结构的非编码RNA分子。已有文献表明，在生物体内，环形RNA有着miRNA海绵、RBP海绵以及翻译短肽等多项功能，在许多生物学过程中发挥着重要作用。 目前研究表明，大部分环形RNA来源于蛋白编码基因的外显子区域。在pre-mRNA剪接的过程中，除典型的内含子剪接事件外，还可能会发生5’端到3’端的反向剪接事件，从而形成环形RNA。因此，剪接产物中环形RNA所占比例是环形RNA分析的重要指标之一，具有高成环比例的环形RNA分子，可能具有更加重要的生物学功能。 同时，同一基因内部也可能产生多种不同的环形RNA，基因内对环形RNA产生位点的使用偏好，也在一定程度上反映了转录过程对环形RNA产生的调控。因此，环形RNA转录本水平的准确定量，是目前环形RNA分析的重要基础。
  为了解决该问题，赵方庆团队开发了一个新的环形RNA分析算法。根据已有工具鉴定出的环形RNA成环位点信息，研究人员重构了具有反向剪接特征的环形RNA参考序列，简化了复杂的反向剪接位点比对问题，并结合测序读段比对到参考基因组和环形序列的结果，筛选出了高置信度的来自环形RNA的读段，解决了目前环形RNA识别和定量方法中准确度低和假阳性率高的问题。作者在模拟数据和真实转录组数据中，对多种常用环形RNA识别软件的表现进行了综合评估，发现该研究中开发的方法在环形RNA表达量和成环比例的计算中，均取得了最佳的结果。

1.2. 分析流程

  将下机测序数据进行质控，去除接头及各类低质量序列。随后借助于CIRIquant，使用Hisat2与参考基因组比对，Stringtie进行基因水平定量；同时使用bwa-men与参考基因组比对，进行circRNA的鉴定，构建circRNA参考序列；将构建的circRNA序列作为参考基因组使用Hisat2再次进行比对，筛选出高置信度的来自环形RNA的reads；统计circRNA的表达情况，并注释circRNA信息。通过对circRNA差异分析，筛选出具有显著差异的circRNA所对相应的基因进行后续富集分析。 circRNA信息分析简易流程如下所示。

2. 分析结果

2.1. 测序数据质量控制

  对原始测序数据及去除接头后的可用数据进行质量评估。测序数据一般为双端测序，因此，每个测序样本会有两个测序结果。

评估的具体内容：

2.2. CIRIquant分析

2.2.1. CIRIquant分析结果文件

1 . 比对结果文件：

以上结果位于文件夹：result/02.CIRIquant/1.mapping/

2 . CIRIquant鉴定结果文件：

以上结果位于文件夹：result/02.CIRIquant/2.circRNA_detection/

3 . CIRIquant鉴定结果的统计结果文件：

以上结果位于文件夹：result/03.circRNA_info

以上统计图的可视化文件：result/03.circRNA_info/view.html

表头说明： （result/03.circRNA_info/1.circRNA_annotation/*csv 鉴定的circRNA的注释信息表）

2.2.2. 参考基因组比对

  测序片段（fragments）是随机打断的，为了确定这些一段由哪些基因转录来，需要将质控后的clean reads比对到参考基因组上。使用HISAT2软件将Clean Reads与参考基因组进行快速精确的比对，获取Reads在参考基因组上的定位信息[4]。HISAT2软件官方手册。

  如果参考基因组组装的较为完善，而且所测物种与参考基因组一致，且相关实验不存在污染，那么实验所产生的测序reads成功比对到基因组的比例会高于70% (Total Mapped Reads or Fragments)。本项目所用参考基因组为 hg38。

2.2.3. circRNA 预测及鉴定

   使用CIRIquant鉴定 circRNA ，并预测 circRNA 的表达。目前发现的circRNAs主要来源于基因外显子exon，但还有其他类型，比如来源于内含子intron、基因间intergenic、反义链antisense。为了更进一步了解鉴定得到的circRNA详细信息，随后进行circRNA类型，circRNA 的长度分布，circRNA 的染色体分布分别进行分析，统计分析图如下。

result/03.circRNA_info/2.circRNA_length/DYQ-HCT116-input.png	result/03.circRNA_info/2.circRNA_length/DYQ-HCT116-target.png
result/03.circRNA_info/3.circRNA_karyotype/DYQ-HCT116-input.png	result/03.circRNA_info/3.circRNA_karyotype/DYQ-HCT116-target.png
result/03.circRNA_info/4.circRNA_type/DYQ-HCT116-input_barplot.png	result/03.circRNA_info/4.circRNA_type/DYQ-HCT116-target_barplot.png
result/03.circRNA_info/4.circRNA_type/DYQ-HCT116-input_pie.png	result/03.circRNA_info/4.circRNA_type/DYQ-HCT116-target_pie.png

2.3. circRNA 差异分析

  对于无重复样本，使用CIRIquant的CIRI_DE工具鉴定差异表达的circRNA。输出的 DE_score 综合了倍数变化和p值，从而提供了一种有效的方法来对差异表达的circRNA排名。此处我们筛选 |DE_score| > 1 作为显著差异表达结果。

2.3.1. 差异分析结果文件

以上结果位于文件夹：result/04.DE/

表头说明： （result/04.DE/targetVSinput_deg*.xls差异分析结果文件）

2.3.2. 差异circRNA的基因组可视化

  可将比对结果bw文件、CIRIquant鉴定得到的circRNA的bed文件、以及差异circRNA分析结果同时放入IGV查看，如：

2.4. 差异circRNA宿主基因富集分析

  我们将差异circRNA的宿主基因，挑选出仅为 protein coding 的基因，用这些基因进行后续富集分析。

  我们根据基因表达量分析得到差异基因之后，必须进一步落到基因的功能上来。对于差异分析而言，往往涉及到成千上万个基因，这会使分析变得很复杂。解决思路是将一个基因列表分成多个部分，从而减少分析的复杂度。为了解决怎么分成不同类，通常会对基因功能进行富集分析, 期望发现在生物学过程中起关键作用的生物通路, 从而揭示和理解生物学过程的基本分子机制。功能富集分析可以将成百上千个基因、蛋白或者其他分子分到不同的通路中，以减少分析的复杂度。另外，在两种不同实验条件下，激活的通路显然比简单的基因或蛋白列表更有说服力。基因功能富集分析首先要构建基因集（ gene set，如 GO 和 KEGG 数据库等），也就是基因组注释信息进行分类。然后再把我们的目标基因集（差异基因集或者其他基因集）映射到背景基因集上，注意区分注释与富集。

  我们采用 clusterProfiler 软件对差异基因集进行 GO 功能富集分析， KEGG 通路富集分析等。富集分析基于超几何分布原理，其中差异基因集为差异显著分析所得差异基因并注释到 GO 或 KEGG 数据库的基因集，背景基因集为所有进行差异显著分析的基因并注释到 GO 或 KEGG 数据库的基因集。富集分析结果是对每个差异比较组合的所有差异基因集、上调差异基因集、下调差异基因集进行富集。本报告中展示的表格是选取某一个比较组合的富集分析结果，图片是部分富集分析结果。

图 5 基因富集分析原理图

2.4.1. 富集分析结果文件

以上结果位于文件夹：result/05.Enrichment/

网页预览图：

• result/05.Enrichment/ALL-pdf.html

• result/05.Enrichment/DOWN-pdf.html

• result/05.Enrichment/UP-pdf.html

表头说明： （result/05.Enrichment/*/gene.ego_*.csv GO富集结果列表）