bcftools csq分析基因突变对蛋白水平的影响

欢迎关注"生信修炼手册"！

csq命令可以分析SNP位点在基因组上的位置，同时还会预测基因突变对编码蛋白的影响。

和其他预测基因突变对蛋白质影响的软件不同，bcftools 将基因组划分为不同的独立区域（和单倍型区域概念类似)，在分析蛋白质变化时，会综合考虑该区域内的所有突变位点，示意图如下

在A图中，该区域包含两个SNP位点，如果单独考虑每个位点，只能预测到氨基酸替换，由精氨酸替换为色氨酸或者谷氨酰胺，综合考虑两个SNP位点时，对应的DNA序列变成了一个终止密码子，蛋白质长度都发生了变化。

在B图中，该区域包含了2个indel 位点，单独考虑每个indel位点时，都是发生了移码突变，氨基酸长度发生了变化，综合考虑两个SNP位点时，氨基酸变化和单独分析一个位点时，又大不一样。

在C图中，两个SNP位点发生在剪切位点两侧，单独考虑每个SNP位点，氨基酸由天冬氨酸替换为天冬酰胺或者谷氨酸，综合考虑两个突变位点时，氨基酸由天冬酰胺替换成赖氨酸。

从示意图可以发现，单独考虑每个SNP位点对于蛋白质的影响，其结果是有偏差的，只有综合考虑邻近范围内所有的突变位点，预测到的蛋白质变化结果才更加可靠。

csq 运行命令如下

bcftools csq -f csq.fa -g csq.gff3 csq.vcf > csq.out

-f参数指定参考基因组的fasta文件，-g参数指定参考基因组的gff3文件，csq.vcf为输入的VCF文件，csq.out为输出文件。

输出文件的格式也是VCF格式，会在INFO列中新增一个BCSQ字段，用来描述突变位点在基因组上的位置和蛋白质序列的变化,示例如下

BCSQ=synonymous|XYZ|ENST00000000001|protein_coding|+|1Y|102C>T

BCSQ的信息由多个字段构成，中间用|连接，包含以下字段

consequence type
基因突变对蛋白影响的类型，包括synonymous, missense, inframe_deletion等类型
gene
基因名称
transcript
转录本名称
biotype
基因类型
strand
正负链信息
amino acid positon
氨基酸的位置
variants list
预测氨基酸变化时，考虑的突变位点的集合

由于bcftools是综合考虑多个突变位点对蛋白质的共同作用，在实际分析时，应该尽可能的过滤掉假阳性的突变位点，然后再分析蛋白水平的影响，这样的分析结果，可信度会更高。

扫描关注微信号，更多精彩内容等着你！

bcftools csq分析基因突变对蛋白水平的影响相关推荐

SIFT score 预测基因突变对蛋白质功能的影响
欢迎关注"生信修炼手册"! SIFT是Sortig Intolerant From Tolerant的缩写, 是一款软件,基于氨基酸序列的同源性和物理性质来预测氨基酸的替换对蛋白质 ...
《基因突变和基因重组》在线教学案例
<基因突变和基因重组>在线教学案例石家庄市第四十九中学赵志国一.教材内容分析 <基因突变和基因重组>是必修二第五章第一节的内容.这节课的内容时主要是通过常态下的教学教材, ...
伴髓系恶性转化的Shwachman-Diamond综合征患儿的临床特征及基因突变分析
摘要:目的探讨伴髓系恶性转化的Shwachman-Diamond综合征(SDS)患儿的临床特征及基因突变情况.方法采用二代测序方法检测11例伴髓系恶性转化的SDS患儿的基因突变情况,分析患儿的临床 ...
研究表明：喝酒“上脸”是基因突变，不仅容易老年痴呆，还容易得胃癌
全世界只有3.14 % 的人关注了爆炸吧知识本文转自科研大匠 "喝酒上脸的人能喝!" 这句话,不管来自天南还是海北的,在酒桌上,肯定都耳熟能详有没有? 其实,喝酒"上 ...
Science子刊：喝酒脸红的人，患胃癌风险大增，他们都有同一个基因突变
撰文:nagashi 导读:喝酒脸红的现象往往只出现在亚洲人中,因此这种现象也被称为"亚洲红".东京大学的研究人员在Science Advances杂志发表研究论文,表明亚洲人常见 ...
双基因突变患者_双任务干预对携带LRRK2基因突变的帕金森病患者手灵活性的影响...
双任务干预对携带LRRK2基因突变的帕金森病患者手灵活性的影响目的:观察执行单.双任务时LRRK2(+)帕金森病(PD)患者手灵活性变化及双任务干预情况.方法:受试者总数为122例,PD患者46例纳 ...
眼液蛋白水平或许可预测阿尔茨海默症
2019年3月8日发表在阿尔茨海默病杂志上的一项新研究发现,阿尔茨海默病(AD)的生物标记物淀粉样蛋白-β和tau蛋白在眼液中的含量降低与低认知评分显著相关.在波士顿医学中心的研究人员的带领下,该研究 ...
基因突变PHP6,基因突变中那些“披着狼皮的羊” 很多“致命性”基因突变正在被证实无害...
在若干年前,很少有人会提出这样一个问题.不过,近年来,医学遗传学经历了一些"自我反省".自本世纪开始以来,基因组研究的快速发展使相关文献充满了上千个同疾病和残疾相关的基因突变.虽然 ...
风湿与类风湿的关系以及风湿病涉及基因突变与基因甲基化的文献收集（仍在更新。。。）
风湿风湿是什么可分为哪几类发病原因风湿与类风湿的关系类风湿关节炎RA血清生物学标志物诊断生物标志物Diagnostic biomarkers for RA 常用的血清生物学标志物新的血清 ...

bcftools csq分析基因突变对蛋白水平的影响

bcftools csq分析基因突变对蛋白水平的影响相关推荐

最新文章

热门文章