Defining window-boundaries for genomic analyses using smoothing spline techniques.
Defining window-boundaries for genomic analyses using smoothing spline techniques.
复制标题
DOI:
10.1186/s12711-015-0105-9
复制
发表时间:
2015-04-17
期刊:
影响因子:
--
通讯作者:
de Leon N
中科院分区:
文献类型:
--
作者:
Beissinger TM;Rosa GJ;Kaeppler SM;Gianola D;de Leon N
High-density genomic data is often analyzed by combining information over windows of adjacent markers. Interpretation of data grouped in windows versus at individual locations may increase statistical power, simplify computation, reduce sampling noise, and reduce the total number of tests performed. However, use of adjacent marker information can result in over- or under-smoothing, undesirable window boundary specifications, or highly correlated test statistics. We introduce a method for defining windows based on statistically guided breakpoints in the data, as a foundation for the analysis of multiple adjacent data points. This method involves first fitting a cubic smoothing spline to the data and then identifying the inflection points of the fitted spline, which serve as the boundaries of adjacent windows. This technique does not require prior knowledge of linkage disequilibrium, and therefore can be applied to data collected from individual or pooled sequencing experiments. Moreover, in contrast to existing methods, an arbitrary choice of window size is not necessary, since these are determined empirically and allowed to vary along the genome. Simulations applying this method were performed to identify selection signatures from pooled sequencing FST data, for which allele frequencies were estimated from a pool of individuals. The relative ratio of true to false positives was twice that generated by existing techniques. A comparison of the approach to a previous study that involved pooled sequencing FST data from maize suggested that outlying windows were more clearly separated from their neighbors than when using a standard sliding window approach. We have developed a novel technique to identify window boundaries for subsequent analysis protocols. When applied to selection studies based on FST data, this method provides a high discovery rate and minimizes false positives. The method is implemented in the R package GenWin, which is publicly available from CRAN.
登录
查看更多内容
影响因子:
3.7
作者:
Qanbari S;Strom TM;Haberer G;Weigend S;Gheyas AA;Turner F;Burt DW;Preisinger R;Gianola D;Simianer H
通讯作者:
Simianer H
影响因子:
2.1
作者:
REINSCH, CH
通讯作者:
REINSCH, CH
影响因子:
64.8
作者:
Rubin, Carl-Johan;Zody, Michael C.;Andersson, Leif
通讯作者:
Andersson, Leif
影响因子:
3.4
作者:
Hider JL;Gittelman RM;Shah T;Edwards M;Rosenbloom A;Akey JM;Parra EJ
通讯作者:
Parra EJ
影响因子:
3.3
作者:
Kelly JK;Koseva B;Mojica JP
通讯作者:
Mojica JP