2019-11-03 00:10:14 热度:

如何简化美化LEfSe分析结果中的Cladogram图

写在前面

关于LEfSe分析,相信大家早已耳熟能详。网上也有很多指导如何做LEfSe分析流程的文章。可是在实际应用中,仍然会遇到一些问题。LEfSe以出图美观的优势吸引大家用它绘图,然而为什么同样的流程,我们做出来的图总是不如别人发在文章里的漂亮?比如,别人发表的图是这样的:

图1Leastdiscriminantanalysis(LDA)effectsizetaxonomiccladogramcomparingallsamplescategorizedbyfourbacterialprovinces.引自WangKai,EnvironmentalMicrobiology,2015

而我做的图是这样的:

图2我做的cladogram图

美颜攻略

下面就来告诉大家如何将图二美化成图一的样子:

首先,观察第一张图,仔细观察后发现该图漂亮的原因是作者只保留了具有显著差异的分类单元分支,而将无差异点(黄色)进行了过滤去除。

在一般流程下基于LEfSe分析获得的树图分支过多,主要由无差异点(黄色)造成。

于是,提示我们可以从LEfSe流程分析的中间文件.lefse_internal_res入手进行编辑:

将LEfSe分析第二步(LDAEffectSize)的结果文件Galaxy12-[B)LDA_Effect_Size(LEfSe)_on_data_11].lefse_internal_res下载下来,使用notepad编辑器打开,该文件记录了每一个分类单元在各组的统计差异显著性结果,打开如下:

我们要做的就是将其中具有显著差异的微生物挑选出来,以每一个差异分类单元为一行,单独制作成一个文件,这需要用到notepad的编辑功能:首先,使用notepad的查找-标记,将差异显著的单元行进行标记;

然后,删除未标记行;

最后,将编辑好的文件保存,并再次上传至LefSe网站,注意上传文件格式type选择lefse_internal_res;

上传后,在分析界面直接选择Plot_Cladogram绘图即可。

按照上述步骤,我们一开始的(图2)分析结果,经优化后如下:

优化后的cladogram图减少了无差异的分类单元的出现,增大了差异微生物的扇面区,结果更加清晰美观。大家可以通过在线设置进一步调整图片的字体、大小等。

如何简化美化LEfSe分析结果Cladogram图

推荐阅读

nanopore测序技术专题(十三):病原微生物快速鉴定
利用nanopore其便携性,建库方便,实时测序等特点,非常适合现场操作。例如当前利用nanopore测序仪已经在野外,高山,冰川,外太空等复杂环境中进行了实时进行鉴定。通过基因组来快速鉴定微生物。主要是通过测序数据与已知数据之间进行快速的序列比对,寻找二者之间基因组上的差异,以及亲缘关系,是否具有共有以及特有基因等。可以利用全基因组测序也可以测全长16S等,根据不同的实验目的选择合适的方法。 微生物快速鉴定原理 一般来说,序列越长则具有越高的唯一性,利用nanopore的长读长特性,可以不用拼接,直接用...[详细]
2019-11-14
nanopore测序技术专题(十二):fastq文件探索
得到fastq格式的nanopore测序数据就可以开始分析了,但是先别急,我们需要先对fastq格式文件进行一下处理,先不要着急拿过来就开始分析,我见过很多人,拿过来数据之后就开始做拼接,然后就等着错误的结果,然后在使用更多时间来找原因。心有猛虎,也要细嗅蔷薇。 fastq文件格式 fastq格式是一种包含质量值的序列文件,其中的q为quality,一般用来存储原始测序数据,扩展名一般为fastq或者fq。目前illumina测序,BGISEQ,IonTorrent,pacbio,nanopore都以fa...[详细]
2019-11-13
nanopore测序技术专题(十一):下载练习数据
如果想要分析nanopore的数据,又没有测序仪测序,该如何进行练习呢?其实这个不难,目前网上有很多已经公开发布的nanopore测序数据,涵盖多个物种,多种测序平台。nanopore官方社区有数据发布的连接,可以从里面下载到原始的fast5格式,也可以从NCBISRA数据库进行下载。这次内容我们就来介绍一下如何下载练习数据。 SRA数据库 网址:https://www.ncbi.nlm.nih.gov/sra/SRA(SequenceReadArchive)数据库是NCBI用于存储二代测序的原始数据,包...[详细]
2019-11-13
Ahjeong Son:截断适配体纳米适配体法高灵敏度检测双酚A
双酚A(BPA)(2,2-双(4-羟基苯基)丙烷)在化学工业中广泛用于聚碳酸酯塑料和环氧树脂的合成。BPA的广泛使用也引起了人们对其对食品安全和环境健康影响的严重关注。BPA被发现存在人体血液、组织、血清和尿液中。作为一种具有雌激素活性的内分泌干扰化合物,双酚a通过干扰雌激素受体结合过程来干扰内分泌系统,所以设计灵敏的生物传感器去检测双酚A在实际应用中具有很强的意义。 先前基于适配体的检测方法的性能改进主要局限于材料、传感器平台和分析程序的变化,然而,适配体本身的设计缺乏改进。基于适配体的检测方法的性能在...[详细]
2019-11-12
nanopore测序技术专题(十):安装所需软件
纸上得来终觉浅,绝知此事要躬行,做生物信息,必须亲自上手才行,否则就会陷入一看啥都会,一做都不对的怪圈。而且,很多科学研究中的重要内容,都隐藏在细节当中,只有在不断的实践中才能发现科学的奥秘,科学研究就是一个不断探索的过程。从这次内容开始,我们将进入nanopore数据分析实战阶段。Areyouready? 软件列表 目前适合nanopore数据分析的软件已经有不少了,虽然没有二代测序多,但是也足够用,而且很多软件并部分测序平台,只要是固定的文件格式即可。下面给出了一些常用的nanopore数据分析的软件...[详细]
2019-11-11
nanopore测序技术专题(九):利用minion_qc绘制质控图
前面我们提到过利用guppy进行碱基识别后生成一个sequencing_summary.txt文件,这个文件是列表格式的统计结果,R语言最喜欢列表格式了,直接使用R可以对这个列表进行统计绘图。其实如果使用官方的MinKNOW软件自带这些绘图功能,如果是自己进行的碱基识别,可以使用minion_qc来进行统计绘图,非常的容易,其实后面我们介绍的NanoPlot也可以使用这个文件绘图。 软件安装 软件官网:https://github.com/roblanf/minion_qc这个工具其实就是一个R脚本,在R...[详细]
2019-11-11
nanopore测序技术专题(九):利用minion_qc绘制质控图
前面我们提到过利用guppy进行碱基识别后生成一个sequencing_summary.txt文件,这个文件是列表格式的统计结果,R语言最喜欢列表格式了,直接使用R可以对这个列表进行统计绘图。其实如果使用官方的MinKNOW软件自带这些绘图功能,如果是自己进行的碱基识别,可以使用minion_qc来进行统计绘图,非常的容易,其实后面我们介绍的NanoPlot也可以使用这个文件绘图。 软件安装 软件官网:https://github.com/roblanf/minion_qc这个工具其实就是一个R脚本,在R...[详细]
2019-11-10
nanopore测序技术专题(八):利用guppy进行碱基识别
目前,绝大部分的生物信息分析是从原始测序的fastq格式文件开始的。而nanopore可以保存为fast5格式,对于这种格式,上一次推文我们已经介绍过了。MinKNOW软件是可以直接进行basecalling输出fastq格式文件的。但是这个过程比较耗时,如果使用MiniION,是很难在笔记本电脑上进行实时的basecalling的。可以将测序完成的fast5文件传输到更大的计算设备上进行处理。如果你只是要做基因组拼接,直接得到fastq文件就行。但有些人缺乏安全感,管他有用没用,觉得信息越多越好,还是想...[详细]
2019-11-09
双11,选择一款适合做生物信息的笔记本
所谓工欲善其事,必先利其器,选择一款好的笔记本电脑来做生物信息是非常有必要的,何况通哥有句名言用一台比较快的电脑,等于延长了生命。作为混迹数码行业10多年的人来说,对各种科技产品当然是了如指掌,如数家珍。前段时间还有人咨询我如何选择一台好的笔记本,正好最近又看到圈内有人求助,这次就让被生物信息耽误的数码编辑带大家选择称心如意的笔电产品吧。 尽量不要用笔记本来跑数据 题目说是推荐适合做生物信息的笔记本电脑,又说尽量不要用笔记本来跑数据,这不是矛盾吗?是的。虽然不能用来跑数据,但是用来练习总是可以的。跑数据尽...[详细]
2019-11-09
nanopore测序技术专题(七):测序结果文件介绍
前面我们已经介绍了很多关于nanopore测序的一些内容,希望你能够对于nanopore测序有所了解。从这次内容开始对数据分析部分进行介绍,至于如何提取你样品的DNA,这部分涉及内容很多,不同的样品有不同的处理方法,我不了解你研究的项目,也不感兴趣,我要是感兴趣就没你什么事了。所以,这部分自己摸索,科学就是不断探索的过程。至于拿到DNA如何建库测序,可以参加nanopore官方提供的实验培训,好像9000多一天,听起来挺贵,但其实可以自带样品测序。如果有需要,自己去关注吧。 获取nanopore测序数据 ...[详细]
2019-11-08
Xshell6的正确打开方式
远程连接服务器的工具有很多,一般都是支持ssh协议,例如putty,mobaxterm,SSHSecureShellClient,secureCRT,这其中我认为最好用的就是Xshell了,不过xshell是收费的,人民币大概899,加上xftp一起购买应该是1299。收费的东西一般都比较好用。不过xshell对于学生可以免费申请使用。 安装 软件只有windows版本,安装比较容易,直接下一步下一步就可以了。 登录设置 登录服务器需要IP地址,账户,密码 1、选择新建会话,在名称中设置一个别名,便于区分...[详细]
2019-11-08
nanopore测序技术专题(六):测序错误率太高无法使用?
十二年前,你嫌illumina测序错误率高,不愿第一个吃螃蟹;七年前,你嫌pacbio的错误率高,觉得应该在等等;四年前,你担心国产的BGIseq错误率高,不愿尝试;现在你还是害怕nanopore测序错误率高,怕投入失败,决定还是坚守二代,等后面测序质量提高的吧。对不起,君生我未生,我生君已老。而且不仅为你写诗,还为你点歌:爱一个测序仪好难 测序错误太高怎么办 按我对测序仪的评价标准,准确性第一位,测序读长第二位,然后是通量,价格,建库测序简易型。为什么准确性第一位,因为如果碱基准确,会给数据分析带来极大...[详细]
2019-11-07
nanopore测序技术专题(五):建库测序
所谓建库测序就是对测序的DNA进行一些处理,是一个格式化的过程,需要将DNA处理成固定的模式才可以。例如需要给原始的DNA加上A碱基,接头,测序引物,index或者barcode标签等。建库测序是DNA测序过程中非常重要的一个环节,可以说直接影响到测序质量,建库效果不好,测序质量不可能好。 DNA提取 纳米孔可以测序DNA本身长度的读长,这就需要原始DNA具有的长度。因此,在测序之前能够提取到相对完整的DNA比较重要。由于不同物种DNA提取方法不同,这就需要依据一些不同样本提取的经验方法。例如植物基因组具...[详细]
2019-11-07
用于食源性致病菌检测的DNA适体传感器研究进展
食源性致病菌,例如沙门氏菌,链球菌,大肠杆菌和弧菌属,是通过食物污染和扩散而引起传染病的有害物种。常见的感染并发症,包括急性胃肠炎腹泻,头痛,呕吐甚至死亡。世界卫生组织宣布,全球每年有数十亿人与食物中毒爆发有关。 近年来,基于核酸的检测分析方法在食源性致病菌的准确检测中得到了广泛的应用。与基于培养的方法相比,聚合酶链式反应(PCR)具有快速性,敏感性和特异性,它是检测各种病原体的广泛使用的标准方法。基于PCR的方法的潜在缺陷是无法区分存活的病原体和死的病原体,因为DNA可以从活的病原体细胞和死的病原体细胞...[详细]
2019-11-06
nanopore测序技术专题(四):纳米孔测序原理
对DNA测序本质上就是识别ATCG四种碱基,但是一方面四种碱基实在太小了,属于纳米级别,另一方面嘌呤和嘌呤,嘧啶和嘧啶之间化学结构非常相似,不容易区分。 从53年提出DNA双螺旋结构之后,生物学家一直努力通过各种办法识别四种碱基。 目前主流的方法包括将四种碱基转换为光信号,溶液PH值,以及转换为电信号,通过放大后的信号来区分四种碱基。 这也是目前主流测序仪的几种方案。 sanger,illumina,BGIseq,Pacbio等选择光信号,Iontorrent选择溶液PH值,而nanopore选择电信号。...[详细]
2019-11-05
制造一台dna测序仪到底有多难?
造一台测序仪到底有多难,具体我也说不上来,但肯定很难很难,和造原子弹相比呢?和造航母比呢?和登上太空相比呢?这样一比就容易了,首先,从时间上来说,原子弹和航母在二战期间就造出来了,61年加加林就登上太空了;再从空间上来讲,能够造原子弹的国家有很多,包括我们的邻居,都能制造个二踢脚;而一般发达国家也都掌握了航天技术。 制造一台测序仪有多难 目前国际上200多个国家,真正能够制造DNA测序仪的只有三个国家,老牌资本主义国家,曾经的日不落帝国英国;目前国际上最强大的国家美利坚合众国,还有一个就是千年文明古国,并...[详细]
2019-11-05
nanopore测序技术专题(三):选择适合你的测序设备
nanopore是非常适合实验室单独运行的测序设备,无需投入大量资本。而不像二代测序,需要昂贵的购买测序仪成本以及运行成本。也许在不久的将来,测序仪将如PCR仪一样普通飞入寻常百姓家。这次内容我们将给大家介绍一下最新的nanopore测序相关产品。 注:本次内容全部来自nanopore官网产品近期介绍,经常会更新,更多内容大家可以自行去nanopore官方查看,或者联系他们公司人员,本次内容无广告赞助(哈哈哈)。产品介绍:https://nanoporetech.com/cn/products查看价格:h...[详细]
2019-11-04
nanopore测序技术专题(二):一些典型应用
前面一篇推文我们介绍了nanopore测序技术的一些显著优势,简单来说就是长读长、高产出、便携、实时、易用、直接。基于这些特点,在基因组或者转录组分析中可以有很多的应用,这次我们列举一些nanopore测序技术的一些典型应用。 大基因组拼接 nanopore最显著特点就是读长长。长读长对于大基因组的拼接将会产生立竿见影的效果。在以往基于短片段的基因组拼接中,由于一些动植物基因组本身具有多倍体,高度重复,高度杂合的特性,导致基因组拼接是一项异常艰难的工作,有些植物甚至复杂到利用短片段根本无法完成拼接工作,例...[详细]
2019-11-03
nanopore测序技术专题(一):为什么要选择nanopore测序?
为什么要选择nanopore测序技术,这是因为nanopore测序技术具有一些无与伦比的优势,可以解决很多技术难题。主要包括,超长读长,DNA/RNA直接测序,真正的实时性,无需对测序设备的资金投入,可扩展性:便携式或台式测序仪,10分钟文库制备,高保真度,对大基因组的高数据量测序。下面我们具体来介绍一下。 1超长读长 在纳米孔测序中,读长长度可以等于输入片段长度。读长长度不受限于测序设备,用户可以通过所使用的文库制备实验方案来控制片段长度。目前报到处DNA片段长度最高记录为2Mb,直接RNA测序读长最长...[详细]
2019-11-03
microbiomeViz:绘制lefse结果中Cladogram
为啥写这个 平日经常会分析shotgun宏基因组的数据,我们的pipeline使用MetaPhlAn,Kraken等profiler。这种数据经常会产生一个表格,如下 download.file(https://bitbucket.org/biobakery/biobakery/raw/tip/demos/biobakery_demos/data/metaphlan2/output/SRS014459-Stool_profile.txt,SRS014459-Stool_profile.txt)knitr:...[详细]
2019-11-03
如何简化美化LEfSe分析结果中的Cladogram图
写在前面 关于LEfSe分析,相信大家早已耳熟能详。网上也有很多指导如何做LEfSe分析流程的文章。可是在实际应用中,仍然会遇到一些问题。LEfSe以出图美观的优势吸引大家用它绘图,然而为什么同样的流程,我们做出来的图总是不如别人发在文章里的漂亮?比如,别人发表的图是这样的: 图1Leastdiscriminantanalysis(LDA)effectsizetaxonomiccladogramcomparingallsamplescategorizedbyfourbacterialprovinces.引...[详细]
2019-11-03
nanopore测序技术专题开篇词
熟悉我们公众号的朋友应该知道我们喜欢写专题,没错,我们喜欢出系统性大作,包括制作的视频也都是四五十集的大作。虽然这样在营销上不合适,但是看起来整齐呀。前面我们推出过《生物数据分析平台搭建-CentOS版》,《生物数据分析平台搭建-Ubuntu版》,《生物信息神奇网站》以及《生物信息百佳软件》(未完待续)。从这开始我们将开启全新的专题《nanopore测序技术专题》,准备好上车吧,跟着通哥不迷路,通哥带你上高速。 基因学苑文章列表(201909) 纳米孔测序技术一项革命性的技术,三年前,我说这是一项来自未来...[详细]
2019-11-01
鸟枪法宏基因组测序之外我们还能做什么?
摘要 环境微生物组的探索揭示了在自然生态系统中起作用的生态和进化原理,重建群体基因组带来的以基因组为中心的研究进一步加速了其发展。 然而, 计算繁重的短读长组装 、 群落内的菌株异质性 以及 低丰度微生物所需的覆盖深度 仍然是传统鸟枪法宏基因组学要应对的技术挑战。 从这个角度出发,我们提出了未来有希望发展的三个主要方向,包括耦合 稳定性同位素示踪技术 与宏基因组、应用 荧光激活细胞分选技术 在较大的群落中靶向寻找微型宏基因组,以及利用 单分子长读长 和 合成长读长技术 将可移动元件与宿主微生物细胞相链接。...[详细]
2019-10-30
宏基因组测序(mNGS)鉴定脑脊液病原体的实验开发及验证
宏基因组测序(mNGS)检测总病原体(pan-pathogen)已经成功地应用在了对病因不明的急性疾病患者的检测中,提供了在一次试验中准确地鉴定几乎囊括了所有潜在的病原体(病毒、细菌、真菌和寄生虫)的方法。随着测序技术和生物信息学技术的发展,利用mNGS进行病原体检测展现出广泛的应用前景,但是仍面临着很多挑战: (1) 缺乏mNGS临床验证的标准; (2) 如何鉴别定植和感染微生物; (3) 缺少用于临床诊断的定制生物信息学软件; (4) 完善现有数据库的质量和全面性; (5) 临床实验室环境的改善。本文...[详细]
2019-10-29
首个亚洲家猪“超清”基因组图谱发布
近日,中国农业科学院农业基因组研究所猪基因组设计育种创新团队唐中林课题组在BioRxiv上预发表了论文Chromosome-scaledenovoassemblyandphasingofaChineseindigenouspiggenome。 该研究使用基因组denovo技术得到了首个超高质量的雄性陆川猪基因组。 由于中国本土猪与西方商品猪在表型和基因组特征上存在显著差异,构建一个有代表性的高质量中国本土猪参考基因组对探索基因功能、基因组进化和促进猪的遗传改良具有重要意义。 贝瑞基因为该研究提供基因组测序...[详细]
2019-10-25
高通量测序与分析软件、数据格式转换三代测序技术的硬核讨论
关于高通量测序与分析软件 Q:各位大神,我们用Nanopore 测序的数据,和参考基因组比较找片段差异,有什么好的软件推荐 A:片段差异是指SV吗? 做这个的话Sniffles还可以,mapping用他们推荐的NGMLR或者minimap2 A:sniffles是nanopore自己的发文章用的;好像还有个nanosv A:我们用pacbio的数据做sniffles,出来的结果很难解释,各种噪声 Q:感谢各位大神,我就是想用nanopore的长reads map一下基因组,看看特定区段上的序列插入缺失,确...[详细]
2019-10-25
全基因组测序探究与脑膜炎相关的侵入性肺炎球菌突变
革兰氏阳性菌肺炎链球菌(Streptococcuspneumoniae)经常定植在人的咽部,通常会侵入无菌的身体部位引起侵入性肺炎球菌病(IPD),包括菌血症、菌血症性肺炎、脑膜炎。IPD是世界范围内发病和死亡的主要原因,而肺炎球菌性脑膜炎是IPD中更为严重的一种,具有高的死亡和永久性神经后遗症风险。肺炎链球菌的突变有时会引起脑膜炎,增加对脑膜炎相关的肺炎链球菌突变的了解,可以提高人们对发病机制的理解,并为预防策略提供信息。本文报告了IPD病人测试队列(n=2054)、独立的验证队列(n=2518)中,肺...[详细]
2019-10-14
  • Xshell6的正确打开方式
    远程连接服务器的工具有很多,一般都是支持ssh协议,例如putty,mobaxterm,SSHSecureShellClient,secureCRT,这其中我认为最好用的就是Xshell了,不过xshell是收费的,人民币大概899,加上xftp一起购买应该是1299。收费的东西一般都比较好用。不过xshell对于学生可以免费申请使用。 安装 软件只有windows版本,安装比较容易,直接下一步下一步就可以了。 登录设置 登录服务器需要IP地址,账户,密码 1、选择新建会话,在名称中设置一个别名,便于区分...
  • nanopore测序技术专题(六):测序错误率太高无法使用?
    十二年前,你嫌illumina测序错误率高,不愿第一个吃螃蟹;七年前,你嫌pacbio的错误率高,觉得应该在等等;四年前,你担心国产的BGIseq错误率高,不愿尝试;现在你还是害怕nanopore测序错误率高,怕投入失败,决定还是坚守二代,等后面测序质量提高的吧。对不起,君生我未生,我生君已老。而且不仅为你写诗,还为你点歌:爱一个测序仪好难 测序错误太高怎么办 按我对测序仪的评价标准,准确性第一位,测序读长第二位,然后是通量,价格,建库测序简易型。为什么准确性第一位,因为如果碱基准确,会给数据分析带来极大...
  • 鸟枪法宏基因组测序之外我们还能做什么?
    摘要 环境微生物组的探索揭示了在自然生态系统中起作用的生态和进化原理,重建群体基因组带来的以基因组为中心的研究进一步加速了其发展。 然而, 计算繁重的短读长组装 、 群落内的菌株异质性 以及 低丰度微生物所需的覆盖深度 仍然是传统鸟枪法宏基因组学要应对的技术挑战。 从这个角度出发,我们提出了未来有希望发展的三个主要方向,包括耦合 稳定性同位素示踪技术 与宏基因组、应用 荧光激活细胞分选技术 在较大的群落中靶向寻找微型宏基因组,以及利用 单分子长读长 和 合成长读长技术 将可移动元件与宿主微生物细胞相链接。...
  • nanopore测序技术专题(三):选择适合你的测序设备
    nanopore是非常适合实验室单独运行的测序设备,无需投入大量资本。而不像二代测序,需要昂贵的购买测序仪成本以及运行成本。也许在不久的将来,测序仪将如PCR仪一样普通飞入寻常百姓家。这次内容我们将给大家介绍一下最新的nanopore测序相关产品。 注:本次内容全部来自nanopore官网产品近期介绍,经常会更新,更多内容大家可以自行去nanopore官方查看,或者联系他们公司人员,本次内容无广告赞助(哈哈哈)。产品介绍:https://nanoporetech.com/cn/products查看价格:h...
  • nanopore测序技术专题(五):建库测序
    所谓建库测序就是对测序的DNA进行一些处理,是一个格式化的过程,需要将DNA处理成固定的模式才可以。例如需要给原始的DNA加上A碱基,接头,测序引物,index或者barcode标签等。建库测序是DNA测序过程中非常重要的一个环节,可以说直接影响到测序质量,建库效果不好,测序质量不可能好。 DNA提取 纳米孔可以测序DNA本身长度的读长,这就需要原始DNA具有的长度。因此,在测序之前能够提取到相对完整的DNA比较重要。由于不同物种DNA提取方法不同,这就需要依据一些不同样本提取的经验方法。例如植物基因组具...
  • nanopore测序技术专题(二):一些典型应用
    前面一篇推文我们介绍了nanopore测序技术的一些显著优势,简单来说就是长读长、高产出、便携、实时、易用、直接。基于这些特点,在基因组或者转录组分析中可以有很多的应用,这次我们列举一些nanopore测序技术的一些典型应用。 大基因组拼接 nanopore最显著特点就是读长长。长读长对于大基因组的拼接将会产生立竿见影的效果。在以往基于短片段的基因组拼接中,由于一些动植物基因组本身具有多倍体,高度重复,高度杂合的特性,导致基因组拼接是一项异常艰难的工作,有些植物甚至复杂到利用短片段根本无法完成拼接工作,例...
  • nanopore测序技术专题(一):为什么要选择nanopore测序?
    为什么要选择nanopore测序技术,这是因为nanopore测序技术具有一些无与伦比的优势,可以解决很多技术难题。主要包括,超长读长,DNA/RNA直接测序,真正的实时性,无需对测序设备的资金投入,可扩展性:便携式或台式测序仪,10分钟文库制备,高保真度,对大基因组的高数据量测序。下面我们具体来介绍一下。 1超长读长 在纳米孔测序中,读长长度可以等于输入片段长度。读长长度不受限于测序设备,用户可以通过所使用的文库制备实验方案来控制片段长度。目前报到处DNA片段长度最高记录为2Mb,直接RNA测序读长最长...
  • nanopore测序技术专题(十三):病原微生物快速鉴定
    利用nanopore其便携性,建库方便,实时测序等特点,非常适合现场操作。例如当前利用nanopore测序仪已经在野外,高山,冰川,外太空等复杂环境中进行了实时进行鉴定。通过基因组来快速鉴定微生物。主要是通过测序数据与已知数据之间进行快速的序列比对,寻找二者之间基因组上的差异,以及亲缘关系,是否具有共有以及特有基因等。可以利用全基因组测序也可以测全长16S等,根据不同的实验目的选择合适的方法。 微生物快速鉴定原理 一般来说,序列越长则具有越高的唯一性,利用nanopore的长读长特性,可以不用拼接,直接用...
  • nanopore测序技术专题(四):纳米孔测序原理
    对DNA测序本质上就是识别ATCG四种碱基,但是一方面四种碱基实在太小了,属于纳米级别,另一方面嘌呤和嘌呤,嘧啶和嘧啶之间化学结构非常相似,不容易区分。 从53年提出DNA双螺旋结构之后,生物学家一直努力通过各种办法识别四种碱基。 目前主流的方法包括将四种碱基转换为光信号,溶液PH值,以及转换为电信号,通过放大后的信号来区分四种碱基。 这也是目前主流测序仪的几种方案。 sanger,illumina,BGIseq,Pacbio等选择光信号,Iontorrent选择溶液PH值,而nanopore选择电信号。...
  • 高通量测序与分析软件、数据格式转换三代测序技术的硬核讨论
    关于高通量测序与分析软件 Q:各位大神,我们用Nanopore 测序的数据,和参考基因组比较找片段差异,有什么好的软件推荐 A:片段差异是指SV吗? 做这个的话Sniffles还可以,mapping用他们推荐的NGMLR或者minimap2 A:sniffles是nanopore自己的发文章用的;好像还有个nanosv A:我们用pacbio的数据做sniffles,出来的结果很难解释,各种噪声 Q:感谢各位大神,我就是想用nanopore的长reads map一下基因组,看看特定区段上的序列插入缺失,确...