使用geneHapR进行基因单倍型分析(以vcf文件为例)

前记

       在群体基因组学研究中,我们常常需要知道一些位点的变异情况,以便于根据对应的表型信息估算这些位点的效应,同时了解这些位点在不同亚群之间的变化情况。这个时候我们就需要进行单倍型分析(Haplotype Analysis),单倍型分析是研究基因组中特定区域的单倍型状态的方法,通过确定个体的单倍型组合及其频率,了解基因座间的关联性以及对特定性状影响。

       geneHapR是2023年发布的一款可用于单倍型分析的R包,由中国农业科学院刁现民老师课题组开发,发表于BMC旗下Bioinformatics期刊,该包的使用方法较为简单,适用于新手入门分析。

       本文以vcf文件为例,演示如何进行单倍型分析。

一、geneHap包简介

1、geneHap的工作流程

如下图所示,可以看到,用户需要准备的包括:基因型文件、注释文件和样本信息文件,以及后续会用的表型文件。输出的主要内容有:变异位点的可视化、单倍型网络、地理分布、表型比较和连锁不平衡热图。

geneHapR的工作流程

2、geneHapR的安装

打开R或Rstudio,输入以下代码进行安装:

library(BiocManager)
BiocManager::install(c("Biostrings", "GenomicRanges", "muscle", "IRanges", "rtracklayer", "trackViewer"))
install.packages("geneHapR")
library(geneHapR)

二、输入文件的设置

        需事先查询好目的基因的物理位置信息,使用bcftools工具提取目的基因的vcf变异信息文件,代码如下:

#bcftools将vcf生成bgzip和index格式
bcftools view my.vcf -Oz -o my.vcf.b.gz
bcftools index my.vcf.b.gz#提取基因的vcf文件
bcftools filter my.vcf.b.gz --regions chr1:3478748-3480748 > mygene.vcf

在R中读取该文件:

vcf <- import_vcf("mygene.vcf")

单纯地进行基因的单倍型分析(即获取基因的不同单倍型分类情况),只需基因的vcf文件!

三、进行单倍型分析

R或Rstudio,输入以下代码:

library(geneHapR)
setwd("D:/Working-Folder/R-work/geneHapR/")vcf <- import_vcf("mygene.vcf")
geneID <- "mygene"      # 基因ID
Chr <- "num"           # 基因所处的染色体名称
start <- start        # 基因的起始位置(染色体坐标)
end <- end          # 基因的终止位置(染色体坐标)
hapPrefix <- "Hap"        # 单倍型名称的前缀
# 从VCF开始单倍型鉴定
hapResult <- vcf2hap(vcf, hapPrefix = hapPrefix,hetero_remove = TRUE, # 移除包含杂合位点的样本na_drop = TRUE) # 移除包含基因型缺失的样本
# 对单倍型结果进行汇总整理
hapSummary <- hap_summary(hapResult, hapPrefix = hapPrefix)# 将单倍型鉴定结果保存到硬盘
write.hap(hapResult, file = "mygene.hapResult")
write.hap(hapSummary, file = "mygene.hapSummary")# 导入之前的单倍型分析结果
hapResult <- import_hap(file = "mygene.hapResult")
hapSummary <- import_hap(file = "mygene.hapSummary")# 以表格形式展示各单倍型的基因型
plotHapTable(hapSummary,             # 单倍型结果hapPrefix = hapPrefix,  # 单倍型名称前缀angle = 45,             # 物理位置的角度displayIndelSize = 0,   # 图中展示最大的Indel大小title = geneID)         # 图片标题

运行结束后,会产生一个表格型图片,显示该基因的单倍型信息。

可以看到,该基因有6种单倍型,最后面一列是每种单倍型对应的样本数目。

四、优异单倍型挖掘

       获取单倍型及对应的样本信息后,再结合对应的表型,进行多重比较,即可获知哪一种单倍型的表型最好。此过程可在Excel活GraphPad Prism等软件进行,过程较为简单,在此不再演示。

五、参考信息

geneHapR做基因单倍型分析-CSDN博客icon-default.png?t=N7T8https://blog.csdn.net/zhang_rl/article/details/130831155

GitHub - ZhangRenL/geneHapRContribute to ZhangRenL/geneHapR development by creating an account on GitHub.icon-default.png?t=N7T8https://github.com/ZhangRenL/geneHapRgeneHapR: an R package for gene haplotypic statistics and visualization | BMC Bioinformatics | Full Text (biomedcentral.com)icon-default.png?t=N7T8https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-023-05318-9

后记

以上是简单使用geneHapR包进行单倍型分析的操作,后续讲述该包的其它功能(基因结构、LD分析、单倍型网络和地理分布等)。

--------2024.4.12

--------CXGG

千里之行,始于足下。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/811905.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Python学习之-Property详解

前言&#xff1a; 在Python中&#xff0c;property 是一种内置的装饰器&#xff0c;它可以将类的方法转换为属性&#xff0c;让你在不改变类接口的情况下添加额外的逻辑&#xff0c;如输入值的验证、取值的计算等。property 可以作为一种方式让你的类接口保持清晰且易于使用。…

Spring Boot | SpringBoot 对 SpringMVC的 “整合支持“

目录: SpringMVC 的 “整合支持” ( 引入"Web依赖启动器"&#xff0c;几乎可以在无任何额外的配置的情况下进行"Web开发")1.SpringMVC "自动配置" 介绍 ( 引入Web依赖启动器"后&#xff0c;SpringBoot会自动进行一些“自动配置”&#xff0…

Linux下mysql的彻底卸载

Linux下mysql的彻底卸载 1、查看mysql的安装情况2、删除上图安装的软件3、都删除成功之后&#xff0c;查找相关的mysql的文件4、删除全部文件5、再次执行命令 1、查看mysql的安装情况 rpm -qa | grep -i mysql2、删除上图安装的软件 rpm -ev mysql-community-libs-5.7.27-1.e…

部署Nginx+keepalived+Tomcat集群架构

部署Nginx+keepalived+Tomcat集群架构 集群架构背景Nginx的高可用负载均衡架构Nginx和LVS对比的总结:搭建Nginx+keepalived+Tomcat的高可用负载均衡架构准备条件: 测试服务器: 系统:centos7.6服务器1: 172.16.2.203 (Nginx+keepalived+Tomcat) 服务器2: 172.16.2.206 (N…

[iOS]协议中如何添加属性?

在 Objective-C 中&#xff0c;协议&#xff08;Protocol&#xff09;本身不直接支持存储属性&#xff0c;因为协议的目的是定义一个接口&#xff0c;即一组方法声明&#xff0c;它们可以被任何类实现。协议主要用于声明方法&#xff08;包括required 和 optional 方法&#xf…

【技术支持】禁止html中referer

如果页面中包含了如下 meta 标签&#xff0c;所有从当前页面中发起的请求将不会携带 referer&#xff1a; <meta name"referrer" content"never"> 如果页面中包含了如下 meta 标签&#xff0c;则从当前页面中发起的 http请求将只携带 origin 部分&a…

蓝桥杯 前一晚总结 模板 新手版

《准备实足&#xff0c;冲冲冲 省一》https://www.yuque.com/lenyan-svokd/hi7hp2/hfka297matrtsxy2?singleDoc# 《准备实足&#xff0c;冲冲冲 省一》 #include<bits/stdc.h> // 包含标准库头文件using namespace std; using ll long long; // 定义 long long 数据类…

【opencv】示例-opencv_version.cpp 输出OpenCV的版本和构建配置的示例

#include <opencv2/core/utility.hpp> // 引入OpenCV核心工具库 #include <iostream> // 引入标准输入输出流库// 定义一个包含命令行参数的字符串 static const std::string keys "{ b build | | print complete build info }" // 定义参数b&#xff…

【数据结构】06图

图 1. 定义1.1 无向图和有向图1.2 度、入度和出度1.3 图的若干定义1.4 几种特殊的图 2. 图的存储2.1 邻接矩阵-顺序存储&#xff08;数组&#xff09;2.2 邻接表-顺序存储链式存储&#xff08;数组链表&#xff09;2.3 十字链表-适用于有向图2.4 邻接多重表-适用于无向图 3. 图…

设计模式代码实战-建造者模式

1、问题描述 小明家新开了一家自行车工厂&#xff0c;用于使用自行车配件&#xff08;车架 frame 和车轮 tires &#xff09;进行组装定制不同的自行车&#xff0c;包括山地车和公路车。 山地车使用的是Aluminum Frame&#xff08;铝制车架&#xff09;和 Knobby Tires&#x…

【随笔】Git 高级篇 -- 管理多分支 git rebase(二十二)

&#x1f48c; 所属专栏&#xff1a;【Git】 &#x1f600; 作  者&#xff1a;我是夜阑的狗&#x1f436; &#x1f680; 个人简介&#xff1a;一个正在努力学技术的CV工程师&#xff0c;专注基础和实战分享 &#xff0c;欢迎咨询&#xff01; &#x1f496; 欢迎大…

vue elementUI form组件动态添加el-form-item rules且支持添加自定义校验方法

vue elementUI form组件动态添加el-form-item rules且支持添加自定义校验方法 组件动态添加el-form-item并且动态添加rules的方法可以参考博客&#xff1a;添加自定义校验方法validatePassFun 组件动态添加el-form-item并且动态添加rules的方法可以参考博客&#xff1a; vue e…

PyTorch环境配置问题

为什么深度学习都是用英伟达的显卡&#xff1f; 首先我们需要了解什么是CUDA&#xff1f; CUDA&#xff08;Compute Unified Device Architecture&#xff09;&#xff0c;是显卡厂商 NVIDIA 推出的运算平台。 CUDA就类似于编程语言&#xff0c;开发者和显卡通过CUDA进行交流…

Android网络抓包--Charles

一、Android抓包方式 对Https降级进行抓包&#xff0c;降级成Http使用抓包工具对Https进行抓包 二、常用的抓包工具 wireshark&#xff1a;侧重于TCP、UDP传输层&#xff0c;HTTP/HTTPS也能抓包&#xff0c;但不能解密HTTPS报文。比较复杂fiddler&#xff1a;支持HTTP/HTTPS…

【SpringBoot】mybatis-plus实现增删改查

mapper继承BaseMapper service 继承ServiceImpl 使用方法新增 save,updateById新增和修改方法返回boolean值,或者使用saveOrUpdate方法有id执行修改操作,没有id 执行新增操作 案例 Service public class UserService extends ServiceImpl<UserMapper,User> {// Au…

用于大规模数据集(大于1TB)的并行运算的MapReduce是怎么实现的?

MapReduce 是一种编程模型&#xff0c;用于处理和生成大数据集。MapReduce 分为两个阶段&#xff1a;Map 阶段和 Reduce 阶段。 Map 阶段&#xff1a;在这个阶段&#xff0c;输入数据被拆分成不同的数据块&#xff0c;这些数据块被分发到各个 Map 任务上。每个 Map 任务对输入的…

Golang | Leetcode Golang题解之第24题两两交换链表中的节点

题目&#xff1a; 题解&#xff1a; func swapPairs(head *ListNode) *ListNode {dummyHead : &ListNode{0, head}temp : dummyHeadfor temp.Next ! nil && temp.Next.Next ! nil {node1 : temp.Nextnode2 : temp.Next.Nexttemp.Next node2node1.Next node2.Nex…

数据结构(算法)

总结&#xff0c;建议看EXCEL的《算法》页签&#xff0c;不然感觉有点乱 备注原理/步骤时间复杂度空间复杂度串的应用模式匹配简单/暴力O(mn) KMP  O(mn) 树的应用树哈夫曼树1、带权路径长度WPL 2、外部排序-最佳归并树1、哈夫曼树的度&#xff0c;只有0和m&#xff08;m叉…

react17+antd4 Menu 点击菜单收起其他展开的所有菜单、页面刷新时设置菜单的选中状态和展开状态

菜单栏展开回收、页面刷新时确保菜单状态与当前页面匹配 1.菜单栏展开和回收事件2. 刷新时默认当前选中项样式的处理2.1 刷新页面菜单保持用户之前的选中状态2.2 配置展开项openKeys的初始值 export const asyncRouterMap [{path: /page1,title: page1,icon: HomeOutlined,}, …

Linux上下载部署zentao v15.5及具体的使用

1.先查询一下Linux的操作系统的位数&#xff0c;确保下载的文件位数与os的一致 [rootlocalhost xiaoming]# uname -m x86_64 [rootlocalhost xiaoming]# getconf LONG_BIT 64 2.下载zentao的Linux压缩包 wget https://www.zentao.net/dl/zentao/15.5/ZenTaoPMS.15.5.zbox…