可以使用R将与通过应用C + T PRS方法(例如,使用PLINK或PRSice-2)获得的P值阈值范围相对应的PRS结果可视化
使用流行的遗传分析工具plink计算PRS-尽管plink不是专用的PRS软件,但是可以使用plink执行使用C + T标准方法计算PRS所需的每个步骤,并且可以执行此多步骤过程,是学习计算PRS的过程的好方法(通常由PRS软件自动执行)。
之前体验了一下[impute.me](http://www.impute.me/ “impute.me”)的基因检测分析服务,其内容之丰富和详尽让我感到吃惊。基本上涵盖了几乎所有有影响的GWAS研究结果(GWAS-Catalog等来源),于是我感到很好奇,他的评估计算方法,于是再次打开其开源的github地址,找到了一些内容,学习一下。再次为这个项目点赞!
人类白细胞抗原(英语:human leukocyte antigen,缩写为HLA),是编码人类的主要组织兼容性复合体(MHC)的基因。其位于6号染色体的短臂上(6p21.31),包括一系列紧密连锁的基因座,与人类的免疫系统功能密切相关。其中部分基因编码细胞表面抗原,成为每个人的细胞不可混淆的“特征”,是免疫系统区分自身和异体物质的基础。HLA基因高度多态化,存在许多不同的等位基因, 从而细致调控后天免疫系统。
根据SNPedia(一个snp维基百科式的网站),以及果壳网上的一篇文章,来看一下我的血型,顺便还可以分辨你的血型抗原基因型,比如A和B型血,到底是AA(BB)还是AI(BI)。下面我就来看一下我的,当然准确率有限。
大概对许多人来说,基因检测很大一部分是检测一下自己的祖源,这样想更多的知道自己生物学意义上来自哪里,基因是不会说慌的。当然,从家谱上是可以看到姓氏来自哪里,如果有数千年不间断的族谱的话,你可以知道你的姓氏来自哪里,尽管可能有过继等的存在。
再来看看第二个项目,同样有惊人地数不清的数目,竟然有339个项目之多,同样覆盖了好多不同的领域,很多还是各种生活特征,充满了趣味性的感觉。同样还是和GWAS计算器一样的,是根据每一篇文章来的,后面是附上了各个位点和等位基因频率等相关系数。
除了祖上有不间断的详细族谱的人(这毕竟是少数大家望族才能有的待遇),想必每个人都会对自己的祖上来自哪里充满了好奇,这也是许多人做基因检测的主要目的,探索未知是人类的天性,于是我来探索以下祖源是怎么得出来的。
一直一来,很好奇祖源分析是怎么做的,还尝试过使用一款SPA(Spatial Ancestry analysis)这个软件尝试进行分析,结果差强人意,祖源定位地点总是在海上。。。(在搜索基因型填充时无意中搜到了这个软件,由几个作者写成,第一作者还是华人,有点骄傲,可惜没有针对中国人或者东亚人的模型,只有世界和欧洲人的模型,瞬间又觉得有点失落,好在许多中国公司在积累了几万以上的数据后也开始进行中国人自己的模型构建,而且做得还不错)。
是个可以让你DIY分析你的基因组的网站,我的基因检测结果是没有提供impute(基因型填充)的,这个网站方便地进行了基因型填充,还有各个基因特征的预测,赞一个!我前面也介绍过,这个网站是开源的,使用R的shiny搭建。当然,这个基因型填 充是基于千人基因组计划进行的,数据结果估计不会像国内测了几十万人的准确。这里还有一个小插曲,前面我的原始结果并不是标准的23andme的,上传时并没有识别,作者还热心地联系我帮助我完成了分析,一并表示感谢。
自从高中毕业选择生物工程专业入行,到研究生阶段微生物学的学习,再到工作开始做基因检测方面的工作,很庆幸基本上在生命科学领域学习和工作。每月去一次的图书馆,让我对科技的好奇心得以保持,互联网的畅达,让我得以接触到全世界的资讯。虽然,在工作之余,我的兴趣大部分在电子、电脑和软件方面。这本《生命的语言》就是借自南图,以一个基因检测从业者身份学习科学大家如何对大众进行基因和生命科学的科普。记录一下读书笔记和心得,并摘抄部分我们可以采取的行动,备忘和便及大家。