标签：判断 chry chrx WES reads test bam 性别
在分析全外数据的时候，无论是在样本鉴定还是CNV分析的过程中，都需要对样本进行性别判断，那么我们如何对WES数据进行性别判断那？
一、性别判断思路
对于WES数据，个人认为性别判断主要有三个思路：
1、根据Y/X染色体的reads比，这个最好理解，男性的Y/X染色体的reads比会显著大于女性的；
2、根据X染色体上SNP的杂合比例，由于女性有两条X染色体，所以女性X染色体SNP的杂合位点的比例会显著大于男性；
3、根据SRY基因，男性的SRY基因上有reads覆盖，而女性没有。
当然，以上的三种方法使用的前提是样本没有性染色体大片段拷贝数异常，比如患有特纳综合征或者克氏综合征的样本用到以上的一些方法可能会出现问题。
基于以上思路我们可以自己写脚本进行性别判定，如果您不会写脚本也没有关系，有很多开源的软件可以实现这个功能。
二、性别判定软件
今天给大家分享的是ngs-bits下的一个工具SampleGender，使用方法如下：
SampleGender -in test.bam -method xy
SampleGender -in test.bam -method hetx
SampleGender -in test.bam -method sry
-in：输入bam文件；
-method：性别判定的方法，xy、hetx、sry分别是基于我上面三个思路的方法，Y/Xreads比例(xy), X染色体杂合比(hetx), SRY基因 (sry)；
软件的分析结果如下，该样本三种方法的判定结果都是女性：
#file gender reads_chry reads_chrx ratio_chry_chrx
test.bam female 9010 3274631 0.0028
#file gender snps_usable hom_count het_count het_fraction
test.bam female 491 of 525 305 186 0.3788
#file gender coverage_sry
test.bam female 0.00
有的时候我们的分析结果会处在灰区，如下：
#file gender reads_chry reads_chrx ratio_chry_chrx
test.bam unknown (ratio in gray area) 145661 1714764 0.0849
这个时候我们可以修改该方法的阈值，比如将xy方法的男性cutoff调整为0.07，再重新判定。
SampleGender -in test.bam -method xy -min_male 0.07
结果如下，性别被判定为男性：
#file gender reads_chry reads_chrx ratio_chry_chrx
test.bam male 145661 1714764 0.0849
对于不同方法cutoff值的选择，建议可以选择100例以上的训练样本进行训练，得到合理的cutoff值。
欢迎大家关注我的公众号
我爱学生信 公众号 主要分享：CNV、转录组、全外显子和全基因组生信分析流程，欢迎各位交流与讨论！
标签：判断,chry,chrx,WES,reads,test,bam,性别 来源： https://blog.csdn.net/qq_42189072/article/details/116796921
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。