| 超详细图文版：如何实现m6A测序数据可视化？ | 您所在的位置：网站首页 › dna测序图怎么分析 › 超详细图文版：如何实现m6A测序数据可视化？ |
超详细图文版：如何实现m6A测序数据可视化？
转录本上线性为内含子区域，蓝色矩形为外显子区域，其中两端的较窄的矩形为UTR区。转录本上箭头方向向右->->->-，则该基因在正链上，左侧的UTR区为5’UTR；转录本上箭头方向向左-
有箭头间隔不均匀的，说明该基因有多个转录本，在此重叠了。可右键点击左侧track名（RefSeq Genes）选择“Expanded”模式查看多个转录本的情况。
有时分析所用的注释文件与IGV默认加载的注释信息有出入，想要对照自己所用的注释信息来查看，可将自己使用的注释文件（.gtf）拖入红框区域。
（云序生物提供的测序服务，使用的注释文件见：Report>Sequence_Results>Alignment）
*出现建议建立索引文件的提示，可以点Go等待自动生成索引文件后再次拖动gtf文件进来（推荐），但如果gtf文件是没有sorted过的（文件名中没有.sorted.字样）则无法建立索引文件，如果gtf文件是由云序生物提供的，可联系销售。也可以点cancel忽略继续进行（可能加载会较慢）。
5、导入样品比对文件
前面准备工作做完，要考察样品的可视化情况，需要样品比对到基因组后得到的bam文件，并且其对应的索引文件bai文件需要与其存放在同一文件夹中，勿随意改名，如果一定要改，注意bam与其bai文件前面部分名称要一致。进行可视化时可以直接使用bam文件，也可以将其转化为tdf格式文件后使用。
方式一：使用bam文件进行可视化
将bam文件拖入红框区域。
样品刚拖入右侧时，有可能只会显示灰色字“zoom in to see features”。这只是提示需要放大到一定程度才会显示出峰。
界面上方可以选择要观察的染色体、输入染色体上具**置或者直接输入基因名，也可以点“-”和“+”来放大和缩小区域。Bam文件在缩放区域时可能会较慢需要耐心等待。
Bam文件加载好后，出现两条track，上面一条coverage即通常可视化图中展示的峰形，下面一条track展现的是每一条read在基因组上的比对情况。
作图时可以右键点击track名选择隐藏这条alignment track，以便让出位置展示多个样品中的甲基化峰。
方式二：使用tdf文件进行可视化
由于bam文件较大，读取和调整较慢，可以通过IGV内自带的IGV Tools将bam文件转化为tdf文件来进行可视化。相较bam文件，tdf文件去除了alignment track的信息只保留了coverage信息，但文件变小便于快速操作，并且用tdf文件操作时可以进行多样品叠加展示（见后文6调整图像）。
在菜单选Tools>Run igvtools...
点击Input File后面的browse，选择要转化的bam文件。
选择好以后，Output File会自动命名，不需要再选择。点击Run进行转化。
*有时出现WARNING: BAM index file...bam.bai is older than ...bam，是提示这个索引文件bai文件比bam文件生成的时间要早，而正常情况下理当先有bam再建立它的bai文件，因此软件提醒注意是否文件有误。一般可能只是由于涉及到原始数据二次下载导致文件生成时间改变，可以忽略这个提示。
点击Run之后等待其转换，需要一定时间。此时可以看到Run变灰，并且等待一段时间后在Message信息框可以观察到“..”点逐渐增加，说明正在转换中。完成后会显示“100%”和“Done”。
生成的tdf文件被保存在bam同文件夹中。在igvtools中重新选择其他Input File可逐个进行转换。
*若tdf文件大小为0或特别小，说明转换未完成或不完整，无法使用，尝试重新进行转换。
tdf文件的使用方式同bam一样，直接将文件拖入IGV。
6、调整图像
导入多个样本后，可以右键track名，进行各种调整和操作，将图像调整到满意的形式展示。可以自己多多尝试。下面例举一些常用的操作供参考。
右键样品名可以改样品名（Rename Track）、更改track的高度、颜色、名字字体大小（Change...）等基本操作。
要对比多个样品中甲基化峰时，会需要将各样品纵坐标scale调整至一致才便于直观比较。
可以 1)根据甲基化峰信号值手动调整个别样品的scale：右键点样本名选Set Data Range后输入*大值；
2）单个样品自动调节至合适的scale：右键点样本名选Autoscale，会自动调整到一个可以显露出该视窗内该样本**信号值的scale；
3）也可以同时选中要统一的所有样品名（按住shift或者control点选）然后右键选择Group Autoscale，则全部样品会被统一采用一个相同的scale。
如果是用tdf文件进行的可视化（见5 方式二），可以将多个样品overlay后用不同的颜色标示，常用在将样品的IP与Input重叠展示。同时选中要重叠的track，右键选择Overlay Tracks。Overlay后注意及时重新命名以免混淆。可以分别Change两个Track的Color （Positive和Negative两条track）.
7、导出图像
在菜单File>Save Image中导出图像。
有三种文件格式可以选择。
导出的图片为当前视窗下的全部内容。
或右键track名，选择Save Image则*导出样品track的图像。
**有文献IGV可视化图上下方有其他个性化的注释信息、图形，甲基化区域有颜色高亮或框，为IGV外画图软件（PS等）自行添加。
**IGV中其他详细信息，可参阅IGV官网http://www.broadinstitute.org/igv/
云序相关产品
m6A RNA甲基化测序 m5C RNA甲基化测序 m1A RNA甲基化测序 m7G RNA甲基化测序 ac4C RNA乙酰化测序 O8G RNA氧化修饰测序 2’-O-RNA甲基化测序 m6Am RNA甲基化测序 5mC DNA甲基化测序 5hmC DNA甲基化测序 6mA DNA甲基化测序 ChIP测序 ATAC测序 云序生物m6A修饰研究五大模块 01 m6A RNA修饰测序 m6A RNA修饰测序（m6A-MeRIP-seq） 对m6A RNA甲基化，目前***的检测手段为m6A-Seq技术，适用于m6A RNA甲基化谱研究，快速筛选m6A RNA甲基化靶基因。云序可提供mRNA和多种非编码RNA的m6A测序： m6A 全转录组测序（涵盖mRNA，LncRNA，circRNA） m6A LncRNA测序（涵盖LncRNA和mRNA） m6A Pri-miRNA测序（涵盖Pri-miRNA和mRNA） m6A mRNA测序 m6A miRNA测序 02 检测整体m6A RNA修饰水平 LC-MS/MS检测整体RNA修饰水平 **高效，可以实现一次检测，9类修饰水平检测，一步到位。 比色法检测整体RNA修饰水平 快速检测m6A整体甲基化水平 03 m6A RNA修饰上游酶的筛选 比色法检测整体RNA修饰水平 寻找上游直接调控m6A RNA甲基化的甲基转移酶。 04 m6A RNA修饰靶基因验证 MeRIP-qPCR 云序提供各类不同修饰的meRIP-qPCR服务，可针对mRNA，lncRNA，环状RNA等不同类型的RNA分子进行检测，低通量验证RNA修饰靶基因表达水平。 05机制互作研究 RIP-seq/qPCR 筛选或验证RNA修饰直接靶点，研究RNA修饰靶基因的调控机制。 RNA pull down -MS/WB 筛选或验证目标RNA互作基因或蛋白，研究相应的分子调控机制。 双荧光素酶实验 验证两基因互作，研究相应的分子调控机制。 ChIP-seq 筛选或验证目标蛋白与DNA互作，研究相应的分子调控机制。 云序生物服务优势 优势一：发表10分以上文章*多的m6A RNA甲基化测序服务平台。云序已累计支持客户发表57篇高水平RNA修饰文章，合计影响因子520分+，是国内支持发文*多、累计影响因子**的公司。 优势二：至今完成4000+例 m6A测序样本，**覆盖医口、农口等各类样本。 优势三：**检测mRNA和各类非编码RNA（circRNA，lncRNA，Pri-miRNA等）。 优势四：**提供m6A一站式服务：m6A整体水平检测、m6A测序、MeRIP-qPCR验证、RIP和RNA pull-down等。 优势五：率先研发超微量MeRIP测序技术，RNA量低至500ng起。 优势六：国内*全的RNA修饰测序平台，提供m6A、m5C、m1A、m7G、m3C、m6Am、O8G、ac4C乙酰化和2'-O-甲基化测序。云序客户RNA修饰部分文章列表
往期回顾 云序生物 | 2021年客户m6A项目文章汇总9-12月m6A**高分文章汇总 1区 IF: 21|云序MeRIP-Seq & RIP-Seq助力肺动脉高压中m6A阅读蛋白YTHDF1机制研究 云序生物**“RNA”甲基化研究汇总——m6A修饰之植物篇 云序项目文章纯测序 |
【本文地址】
公司简介
联系我们