GWAS分析中0-1-2的编码问题

最新推荐文章于 2024-06-25 20:30:56 发布

育种数据分析之放飞自我

最新推荐文章于 2024-06-25 20:30:56 发布

阅读量2.5k

点赞数 1

分类专栏： GWAS

互秀朋友圈呀, 微信: yijiaobani, 备注: 姓名+ 单位

本文链接：https://blog.csdn.net/yijiaobani/article/details/115916562

版权

GWAS 专栏收录该内容

92 篇文章 225 订阅

订阅专栏

昨天一个小伙伴给我写了封信，问我了几个GWAS的问题，我就回信说，答案可以写一下公众号了。

这不，就来了。

邓老师，您好，看了好多博文，学到好多，最近看一些GWAS原理相关的内容，有些疑惑。想请教一下。

1、GWAS分析基本都是基于回归分析，分析时势必会将序列转化为012编码的基因型，这个编码格式什么规格？Major allel
0；杂合1；minor allel ：2?，还是其他规则?仅考虑加性模型。

2、看这个好多软件大多都支持Plink的格式，是否在其他软件（GAPIT、GEMMA等？）进行转化过程中也遵循Plink的编码格式？

3、我是做植物的，GAPIT计算结果输出的Effect
是否就是PlinkGLM模型中的beta值？都表示回归系数？代表单个位点的改变性状的变化？

4、回归系数反应优异等位基因的来源，这和常用的VCF文件中的REF 和ALT 基因型有没有对比，是否Major
allel一定是REF，而Minor allel基因型一定是ALT基因型？

十分感谢！

第一个问题：GWAS中0-1-2的编码？

应该都是采用这样的编码：主等位基因纯合编码为0，杂合为1，次等位基因纯合编码为2。

看一下plink的官网说明：

–recode creates a new text fileset, after applying sample/variant filters and other operations. By default, the fileset includes a .ped
and a .map file, readable with --file.

The ‘12’ modifier causes A1 (usually minor) alleles to be coded as ‘1’
and A2 alleles to be coded as ‘2’, while ‘01’ maps A1→0 and A2→1.
(PLINK forces you to combine ‘01’ with --[output-]missing-genotype
when this is necessary to prevent missing genotypes from becoming
indistinguishable from A1 calls.) The ‘23’ modifier causes a
23andMe-formatted file to be generated. This can only be used on a
single sample’s data (a one-line --keep file may come in handy here).
There is currently no special handling of the XY pseudo-autosomal
region. The ‘AD’ modifier causes an additive (0/1/2) + dominant (het =
1, otherwise 0) component file, suitable for loading from R, to be
generated. ‘A’ is the same, except without the dominance component. By
default, A1 alleles are counted; this can be customized with
–recode-allele. --recode-allele’s input file should have variant IDs in the first column and allele IDs in the second. By default, the
header line for .raw files only names the counted alleles. To include
the alternate allele codes as well, add the ‘include-alt’ modifier.
Haploid additive components are 0/2-valued instead of 0/1-valued, to
maintain a consistent scale on the X chromosome.

测试一下，手动生成一个模拟的plink文件

ped数据：

1 1 0 0 1  0  G G  2 2  C C
1 2 0 0 2  0  A G  0 0  A C
1 3 1 2 1  2  0 0  1 2  A C
2 1 0 0 1  0  G G  2 2  0 0
2 2 0 0 2  2  A A  2 2  0 0
2 3 1 2 1  2  G G  2 2  C C

map数据：

1 snp1 0 1
1 snp2 0 2
1 snp3 0 3

查看一下频率：

plink命令：

plink --file tt --freq --out re1

结果：

$ cat re1.frq
 CHR  SNP   A1   A2          MAF  NCHROBS
   1 snp1    A    G        0.375        8
   1 snp2    1    2            0        6
   1 snp3    A    C         0.25        4

结果中，A1是次等位基因，A2是主等位基因。

转化为0-1-2的形式：

plink命令：

plink --file tt --recodeA --out re2

结果：

FID IID PAT MAT SEX PHENOTYPE snp1_A snp2_1 snp3_A
1 1 0 0 1 -9 0 0 0
1 2 0 0 2 -9 1 NA 1
1 3 1 2 1 2 NA 1 1
2 1 0 0 1 -9 0 0 NA
2 2 0 0 2 2 2 0 NA
2 3 1 2 1 2 0 0 0