以水稻为例教你如何使用BSA方法进行遗传定位（下篇）

接着上一篇留下的问题，如何按照以1M为窗口，每次移动10Kb计算均值。我们以KY131是"0/0", 而 "DN422" 是 "1/1"的位点结果为例，我们来增加这条线

原本我的想法是使用split根据每个位点的位置信息进行拆分，但是这里有一个问题是，如果只是每个窗口没有重叠的话，那么每个位点应该都有一个唯一的因子，但这里是以10kb将长度为1Mb的窗口从左往右移动，也就意味着有很多位点时属于多个窗口的，也就是无法直接用split-apply-combine的方法进行处理。

于是乎，我写了一个函数用于根据窗口来计算均值，它的逻辑就是先创建好窗口，然后根据窗口的起始位置和终止位置去筛选符合要求的位点，最后进行求和。同时考虑有些位置可能没有SNP，那么最后还进行了一波过滤。

函数的输入是之前snp-index的位置和对应的值，以及确定窗口的大小和步长，

calcValueByWindow <- function(pos, value,window_size = 1000000,step_size = 100000){# get the max position in the postionmax_pos <- max(pos)# construct the windowwindow_start <- seq(0, max_pos   window_size,step_size)window_end <- window_start   step_sizemean_value <- vector(mode = "numeric", length = length(window_start))# select the value inside the windowfor (j in seq_along(window_start)){pos_in_window <- which(pos > window_start[j] &pos < window_end[j])value_in_window <- value[pos_in_window]mean_value[j] <- mean(value_in_window)}# remove the Not A Number positionnan_pos <-  is.nan(mean_value)mean_value <- mean_value[! nan_pos]window_pos <- ((window_start   window_end)/ 2)[!nan_pos]df <- data.frame(pos   = window_pos,value = mean_value)return(df)
}

得到结果就可以用lines在之前结果上加上均值线

par(mfrow = c(3,4))for (i in paste0("chr", formatC(1:12, width = 2, flag=0)) ){freq_flt <- freq2[grepl(i,row.names(freq2)), ]pos <- as.numeric(substring(row.names(freq_flt), 7))snp_index <- freq_flt[,1] - freq_flt[,2]# bindf <- calcValueByWindow(pos = pos, value = snp_index)plot(x = pos, y =snp_index, ylim = c(-1,1),pch = 20, cex = 0.2,xlab = i,ylab = expression(paste(Delta, " " ,"SNP index")))lines(x = df$pos, y = df$value, col = "red")
}

最后再对文章总结一下。文章并不是只用了BSA的方法进行定位，他们花了几年的时间用SSR分子标记确定了候选基因可能区间，用BSA的方法在原有基础上缩小了定位区间。当然即便如此，候选基因也有上百个，作者通过BLAST的方式，对这些基因进行了注释。尽管中间还加了一些GO富集分析的内容，说这些基因富集在某个词条里，有一个是DNA metabolic processes(GO:0006259)，但我觉得如果作者用clusterProfiler做富集分析，它肯定无法得到任何富集结果。他做富集分析的目的是其实下面这个描述，也就是找到和抗冻相关的基因

LOCOs06g39740 and LOCOs06g39750,were annotated as the function of “response to cold (GO: 0009409)”, suggesting their key roles in regulating cold tolerance in rice. "

当然他还做了qRT-PCR进行了验证，最后推测LOC_Os06g39750应该是目标基因，这个基因里还有8个SNP位点。

以水稻为例教你如何使用BSA方法进行遗传定位（下篇）相关推荐

以水稻为例教你如何使用BSA方法进行遗传定位（上篇）
BSA虽然不是我最早接触的高通量数据分析项目(最早的是RNA-seq),但是却是我最早独立开展的数据分析项目, 我曾经专门写过一篇文章介绍如何使用SHOREMap做拟南芥的EMS诱变群体的BSA分析 ...
如何使用BSA方法进行遗传定位（水稻篇）
BSA虽然不是我最早接触的高通量数据分析项目(最早的是RNA-seq),但是却是我最早独立开展的数据分析项目, 我曾经专门写过一篇文章介绍如何使用SHOREMap做拟南芥的EMS诱变群体的BSA分析 ...
MPB：农科院田健、韩东飞等-水稻根系互作功能微生物的筛选方法
为进一步提高<微生物组实验手册>稿件质量,本项目新增大众评审环节.文章在通过同行评审后,采用公众号推送方式分享全文,任何人均可在线提交修改意见.公众号格式显示略有问题,建议电脑端点击文末阅 ...
两部苹果手机同步照片_怎么恢复苹果手机删除的照片？今天教你三种找回方法...
原标题:怎么恢复苹果手机删除的照片?今天教你三种找回方法怎么恢复苹果手机删除的照片?手机的出现虽然带给我们很大便利,同时却也带来了一些小的麻烦.在手机上很多操作步骤都很简单,因此,难免会遇到手滑误操 ...
机器学习：一步步教你理解反向传播方法
机器学习:一步步教你理解反向传播方法时间 2016-09-13 00:35:59 Yong Yuan's blog 原文 http://yongyuan.name/blog/back-propa ...
iphone屏蔽系统更新_iPhone手机经常提示更新系统，教你一招关闭方法，学到了
用过iPhone手机的都知道,一段时间之后就会不停地有消息更新提醒,让你自动更新手机系统,一段时间的置之不理之后,你发现你的手机已经自动更新到最新系统了. 是不是让人很头大呢?下面就来教教大家怎样关闭 ...
测试电视是不是4k的软件,怎么判断4K电视真假?教你快速检测的方法!
原标题:怎么判断4K电视真假?教你快速检测的方法! 4K电视从进入市场之后一直都受到企业的力捧,随着电视企业对4K电视的大力度宣传和消费环境的逐渐成熟,越来越多的消费者开始认可4K电视,并在购机时表明 ...
java getipaddress_教你java用getAddress方法取得IP地址
本篇教你java用getAddress方法取得IP地址: getAddress方法和getHostAddress类似,它们的唯一区别是getHostAddress方法返回的是字符串形式的IP地址,而g ...
python 暂停程序等待用户输入_遇上Python程序暂停时，不要慌，教你正确的处理方法...
今天为大家带来的内容是:遇上Python程序暂停时,不要慌,教你正确的处理方法! 文章内容主要介绍了Python程序暂停的实现代码,非常不错,具有一定的参考借鉴价值,需要的朋友可以参考下,喜欢的记得点 ...

以水稻为例教你如何使用BSA方法进行遗传定位（下篇）

以水稻为例教你如何使用BSA方法进行遗传定位（下篇）相关推荐

最新文章

热门文章