java fit 16s,科学网—16s rRNA分析流程和工具的介绍 - 肖斌的博文

2023-12-28 04:20

本文主要是介绍java fit 16s,科学网—16s rRNA分析流程和工具的介绍 - 肖斌的博文,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

“split_libraries.py” 和“split_libraries_fastq.py”实现数据拆分和数据过滤的双重目的。Mothur利用“Trim.seqs”。不过QIIME和Mothur都不能直接处理sff文件(454下机产生的数据格式),不过可各自利用“process_sff.py”和Sffinfo将sff格式转换为FASTA和QUAL文件。

数据过滤考虑的参数有:minimum average quality score allowed in a read、maximum number of ambiguous bases allowed、minimum and maximum sequence length、maximum length of homopolymer allowed、maximum mismatches inprimer or barcode allowed、whether to truncate reverse primer, and so on.

(2)Denoise and chimera checking

16s建库的pcr过程、测序过程均会导致序列出现错误,在分析过程过程中需要有效排除这种错误。测序误差矫正常用的工具有Denoiser(implemented in QIIME)、AmpliconNoise、Acacia、Pre.cluster(implemented in Mothur)。嵌合体查找的工具有ChimeraSlayer、UCHIME、Persus、DECIPHER,ChimeraSlayer、UCHIME、Persus在mothur中均可调用。在这些工具中,存在有待于优化的问题(these different methods often disagree with one another on the list of identified chimeras,probably because of their different mechanisms or algorithms. More efforts are required to evaluate these methods and coordinate their inconsistencies in chimera identification.)

328d5170a2a19271a5e955cef593521f.png

在分析中有个关于古细菌序列的情况需要注意:a very small proportion of archaeal sequences may be generated for 16S rRNA gene amplicon datasets amplified with bacteria-specific primers. These unexpected sequences should be identified after denoising and chimeraremoval, and are advised to be discarded before subsequent data normalization.

(3)Data normalization

测序深度不理想和不均匀的话会对alpha多样性及beta多样性均有影响。Uneven sequencing depth can affect diversity estimates in a single sample (i.e.,alpha diversity), as well as comparisons across different samples (i.e., beta diversity), thus data normalization is required. 对于此问题有两种处理策略,分别是relative abundance and random sampling (i.e., rarefaction),in addition, z-score亦用于normalization的过程中。但不同的方法均会有缺点。

(4)Picking OTUs and representative sequences

OTU的界定主要根据序列的一致性进行,(The OTUs are picked based on sequence identity, and various identity cutoffs of 16S rRNA gene have been used for different taxonomic ranks. For example, identity cutoffs recommended by MEGAN are 99 % for species, 97 % for genus,95 % for family, and 90 % for order level, respectively)。OTU界定时选择的工具与算法对后期分析有很大影响(The OTU picking strategy and algorithms have significant effects in the downstream data interpretation. )根据分析过程中是否使用数据库,OTU界定的策略可分为de novo、closed reference和open reference。在  OTU界定中有很多聚类的方法,There are many clustering or alignment tools available for OTU picking, such as Uclust, cd-hit, BLAST, mothur, usearch, and prefix/suffix. These tools are implemented in QIIME. Among them, the mothur method contains three clustering algorithms to pick de novo OTUs, namely, nearest neighbor, furthest neighbor, or average neighbor. 当序列聚类好后,代表了一个OTU,接下来就是从这个OTU找到代表序列,一种做法是a representative sequence can be a random, the longest, the most abundant(as default in QIIME), 另一种操作方法是the first sequence in an OTU cluster.  还有一种策略是the distance method in mothur identifies the   sequence with the smallest maximum distance to the other sequences as the representative sequence.

(5)Taxonomic assignment

taxonomic assignment的策略有:(1)word match,如RDP classfier,(2)best hit,(3)Lowest Common Ancestor,如MEGAN、SINA Alignment Service.

(6)Phylogenetic analysis

Phylogenetic relationships一般可以用树来表示,phylogenetic relationships主要是通过序列比对来实现的,序列比对的工具有ClustalW, MUSCLE, Clustal Omega, Kalign, T-COFFEE, COBLAT和FastTree. 目前针对16s rRNA NGS数据的分析工具都可以实现,如MEGA,RAxML,MRBAYES,PhyML,TreeView,Clearcut,FitTree. 其中RAxMLand PhyML are the most widely used programs for maximum-likelihood phylogenetic analysis, probably because they are specifically designed and optimized for such purpose.

(7)Alpha- and beta-diversity analyses

alpha多样性有众多指标可以表示,在mothur中有Shannon, Berger-Parker,Simpson, Q statistic; observed richness, Chao1, ACE, and jackknife。而在QIIME中,有phylogenetic diversity (PD)-whole tree, chao1, and observed species.

还有一种物种丰度的比较方法:rarefaction curve. QIIME中主要用“single_rarefaction.py”、 “multiple_rarefaction.py”,在mothur中主要用“Rarefaction.single”和“Rarefaction.shared”.

beta多样性计算主要反映不同样本之间的差异度,several distance metrics, such as Unifrac, Bray-Curtis, Euclidean,Jaccard index, Yue & Clayton, and Morisita-Horn, have been often employed. beta多样性计算根据是否考虑OTU的相对丰度,可分为定量指数和定性指数。

(8)Statistical and network analysis

在Two-sample/group中,多考虑t-test。在其中需要注意,Particularly for independent two-samplet-test, independence and equal variances (which canbe tested by F-test, Levene’s test, etc.) of two populations arerequired. In the case of non-normal distribution of data sets,nonparametric two-sample tests robust to data non-normality,such as Wilcoxon signed-rank test, and Mann-Whitney U testare applicable for significance testing of difference betweengroup medians.

在Multiple-sample/group tests中,ANOVA。

(9)Clustering and classification

clustering可以分析样品之间的亲疏关系。classfication的策略用来对样品进行类别判定。

(10)Ordination analysis

在样本的相似度和距离计算完后,可以利用principal component analysis (PCA), principal coordinates analysis(PCoA, also known as metric multidimensional scaling), Nonmetric multidimensional scaling (NMDS), canonical correspondence analysis (CCA), linear discriminantanalysis (LDA), and redundancy analysis (RDA)等构建样本间的关系。

(11)Network-based modeling

与基因表达、代谢分子、蛋白等数据一起分析共表达网路或者共表达模式(co-occurrence and co-exclusion patterns)

参考文章:JuF,ZhangT.16srRNAgenehigh throughputsequencingdataminingofmicrobiotadiversityandinteractions,ApplMicrobiolBiotechnol.2015,99(10):4119-4129

这篇关于java fit 16s,科学网—16s rRNA分析流程和工具的介绍 - 肖斌的博文的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/545055

相关文章

Java使用Curator进行ZooKeeper操作的详细教程

《Java使用Curator进行ZooKeeper操作的详细教程》ApacheCurator是一个基于ZooKeeper的Java客户端库,它极大地简化了使用ZooKeeper的开发工作,在分布式系统... 目录1、简述2、核心功能2.1 CuratorFramework2.2 Recipes3、示例实践3

Springboot处理跨域的实现方式(附Demo)

《Springboot处理跨域的实现方式(附Demo)》:本文主要介绍Springboot处理跨域的实现方式(附Demo),具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不... 目录Springboot处理跨域的方式1. 基本知识2. @CrossOrigin3. 全局跨域设置4.

springboot security使用jwt认证方式

《springbootsecurity使用jwt认证方式》:本文主要介绍springbootsecurity使用jwt认证方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地... 目录前言代码示例依赖定义mapper定义用户信息的实体beansecurity相关的类提供登录接口测试提供一

Spring Boot 3.4.3 基于 Spring WebFlux 实现 SSE 功能(代码示例)

《SpringBoot3.4.3基于SpringWebFlux实现SSE功能(代码示例)》SpringBoot3.4.3结合SpringWebFlux实现SSE功能,为实时数据推送提供... 目录1. SSE 简介1.1 什么是 SSE?1.2 SSE 的优点1.3 适用场景2. Spring WebFlu

基于SpringBoot实现文件秒传功能

《基于SpringBoot实现文件秒传功能》在开发Web应用时,文件上传是一个常见需求,然而,当用户需要上传大文件或相同文件多次时,会造成带宽浪费和服务器存储冗余,此时可以使用文件秒传技术通过识别重复... 目录前言文件秒传原理代码实现1. 创建项目基础结构2. 创建上传存储代码3. 创建Result类4.

Java利用JSONPath操作JSON数据的技术指南

《Java利用JSONPath操作JSON数据的技术指南》JSONPath是一种强大的工具,用于查询和操作JSON数据,类似于SQL的语法,它为处理复杂的JSON数据结构提供了简单且高效... 目录1、简述2、什么是 jsONPath?3、Java 示例3.1 基本查询3.2 过滤查询3.3 递归搜索3.4

Tomcat版本与Java版本的关系及说明

《Tomcat版本与Java版本的关系及说明》:本文主要介绍Tomcat版本与Java版本的关系及说明,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录Tomcat版本与Java版本的关系Tomcat历史版本对应的Java版本Tomcat支持哪些版本的pythonJ

springboot security验证码的登录实例

《springbootsecurity验证码的登录实例》:本文主要介绍springbootsecurity验证码的登录实例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,... 目录前言代码示例引入依赖定义验证码生成器定义获取验证码及认证接口测试获取验证码登录总结前言在spring

SpringBoot日志配置SLF4J和Logback的方法实现

《SpringBoot日志配置SLF4J和Logback的方法实现》日志记录是不可或缺的一部分,本文主要介绍了SpringBoot日志配置SLF4J和Logback的方法实现,文中通过示例代码介绍的非... 目录一、前言二、案例一:初识日志三、案例二:使用Lombok输出日志四、案例三:配置Logback一

springboot security快速使用示例详解

《springbootsecurity快速使用示例详解》:本文主要介绍springbootsecurity快速使用示例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录创www.chinasem.cn建spring boot项目生成脚手架配置依赖接口示例代码项目结构启用s