集全基因组测序、全面代谢组学和先进成像于一体的精准医疗

为了调查全基因组水平的致病基因和变异的价值和临床影响,我们使用了一项前瞻性队列研究设计,招募了同意对其全基因组进行测序并使用临床实验室测试、代谢组学技术和先进的非侵入性成像参与深度表型分析的志愿者。基因组结果与表型结果相结合。

继续阅读 →

ubiome类似数据处理探索8

前面做的许多处理基本上自己拼凑来的,下面再看下完整解决方案。researchgate网站上有人说qiime1版本有这个双向数据配对不拼接的选项?这个没找到。主要发现了有两个方案,一个是有篇文章提出了一个流程*Hybrid-denovo*,还有一篇peer review的文章,几个人评议还有一个人不同意,anyway,都看下。

继续阅读 →

2019新冠状病毒学习笔记

最近这波疫情,重现当年初中非典时期,甚至愈演愈烈,与之前初中时的封校住宿学习不同,已经工作的今天和太多的互联网信息大爆炸让我们有些焦虑,特别是,作为学习生物的人,我们也感到无能为力。官方媒体的科普,已经让大家对这个病毒的具体情况有所了解。我注意到,NJEM也已经把许多文章翻译成了中文版,以正视听。在这个时候,我们不能听信谣言!那么作为有些生物学素养的我们,也应该以自己的知识,学习下这个病毒的信息,以我们自己的理解!

继续阅读 →

ubiome类似数据dada2处理探索7

接着前面的内容,这里再进行下数据库的处理,看看从参考数据库就按测序数据处理是不是能提高物种注释的精度。这里先预报一下,种的分类结果并不能有明显的提升,或许是因为序列长度的缺陷,即使再努力提高技巧,终究不能解决根本的问题,250bp的长度,对比1500bp左右的全长,显然还是太短了,难以实现精确的分类,所以,要想更精确,只有上16S全长,这只能寄望于Pacbio,Oxford Nanopore,和10x linked reads或者类似的技术,比如华大的sLtFR等技术提升读长了。再激进些,等测序成本足够低,上宏基因组,宏转录组了。

继续阅读 →

ubiome类似数据dada2处理探索6

前面我们探索了处理不能拼接的V4 PE150数据,首先双向reads根据质量情况分别切成120bp,然后使用dada2 R包进行了直接+10N拼接,生成ASV表,再分别使用dada2包,decipher包和qiime2进行了物种注释,基本上完成了一个最简单的分析过程。这里填下自己之前挖的坑,比较一下这个含有348条序列的样本,qiime2,dada2和的分类器哪个效果更好。

继续阅读 →

ubiome类似数据dada2处理探索5

前面我们探索了处理不能拼接的V4 PE150数据,首先双向reads根据质量情况分别切成120bp,然后使用dada2 R包进行了直接+10N拼接,生成ASV表,再分别使用dada2包和qiime2进行了物种注释,基本上完成了一个最简单的分析过程,这里,使用比较流行的phyloseq包进行下多样性分析。

继续阅读 →