男性福音:吾射不亦精乎?

您还在为妻子的埋怨眼神疲倦吗?您还在为家人的催促苦恼吗?你还在糊里糊涂晕了吧唧不知所以然嘛?咳咳,今天就为大家普及一个男性福音的羞羞哒技能! 咳咳,男科疾病怎能困住我的生信男......

最近!一个偶然的机会,发现了一个男性福音,咳咳,不(jiu)是故意的!!!

故事还要从那个风雨交加电闪雷明而又孤独的夜说起......

那天

夜里

读了一篇文献

这里插播一条广告:早早先,我把寄几的vcf文件全部注释到了dbSNP的rs ID号,检索非常方便。一般基因有rs ID号的位点代表这个基因是研究比较多的,这个大家都知道。 如果不清楚的可以看我们前文的直播基因组系列:

【直播】我的基因组 31:vcf文件标记dbSNP的rsID号

好了,故事又回到那个凄苦的夜晚:作为苦逼生物科研狗的我在paper里面看到了这个类DNA甲基转移酶蛋白3(DNMT3L)基因,想找一些中文资料看看这个基因被研究的如何。

纳尼!排在首位的文献里居然提到了DNMT3L基因上面的rs2070565位点的特殊功能,我马不停蹄滴特意去SNPmedia里面搜索了:

[PMID 22116073] Association between single-nucleotide polymorphisms of DNMT3L and infertility with azoospermia in Chinese men

故事朝着无法逆转的结局奔去!

人生一大憾事是其始也至善,一日不如一日,至其终则坏无可坏。by 马克·吐温

控制不住寄几,我拿这个位点去我的VCF文件里面搜索了一下。结果:

21 45681153 rs2070565 T C 228.0 . DP=78;VDB=0.451144;SGB=-0.693147;MQSB=1;MQ0F=0;AC=2;AN=2;DP4=0,0,37,34;MQ=60;ASP;CAF=0.2139,0.7861;COMMON=1;G5;GENEINFO=DNMT3L:29947;GNO;HD;INT;KGPhase1;KGPhase3;PM;PMC;RS=2070565;RSPOS=45681153;SAO=0;SLO;SSR=0;VC=SNV;VLD;VP=0x05012808000515053f000101;WGT=1;dbSNPBuildID=96 GT:PL 1/1:255,214,0

吓死宝宝了!

DP4=0,0,37,34 ——说明我这个变异位点很可靠,没有正负链的偏差。

上面标记着G5,说明5% minor allele frequency in each and all populations 这是个常见的variation。

标记KGPhase1;KGPhase3;说明在千人基因组计划里面也出现了。

其它的TAG不那么重要,我就不一一解读了!

VCF格式的具体含义大家可以看:

【直播】我的基因组28-必须要理解vcf格式记录的变异位点信息

我抱着求学上进,求知若渴的心态仔细阅读了文章:

方法:应用PCR-RFLP分析方法,在198名少精症患者(少精症组)和249 名正常的男性个体(正常对照组)中,对DNMT3L基因rs2070565位点的基因频率和基因型频率分布进行调查。结果:在少精症和正常男性中 rs2070565位点的等位基因频率分布存在差异,少精症组中等位基因A的频率显著高于正常对照组(20.7%vs 14.9%,P=0.03)。

很明显,参考基因组在这里是T,而我变成了C。难道我.......????

幸亏!!我的求知精神再次感动我寄几:文章中的研究者得出的结论是——如果该位点是A(T)才是少精症概率高,那么我是C!我是正常的!!!哈哈哈哈.....喜悦之情溢于言表!而且,即使该位点是A(T),那也只是风险稍微高那么一丁点而已啦!

哦,虚惊一场~

然后,我又马不停蹄地去了千人基因组计划里面看人群频率。变成C的是大多数人,反而跟参考基因组一样的才是少数人。所以,怕什么!?不过,你不担心你就是那少数人吗?(想知道自己在该位点是什么吗?送你3个字,联系我!)

编者导语

怎地?要算一卦么?

解释一下:题目那句文言取自文言文《卖油翁》不是你们想得那个样子啦~~~

原文发布于微信公众号 - 生信技能树(biotrainee)

原文发表时间:2017-01-14

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏生信宝典

肿瘤化疗无效是对预先存在的突变的选择还是诱发新突变,Cell给你答案

1094
来自专栏思影科技

揭示受伤大脑隐藏认知的EEG信号特征

来自美国哈佛医学院的William H. Curley试图探索能够作为意识障碍患者认知能力指标的EEG信号特征,并检验将来用于重建沟通的可靠性。该研究发表在《B...

2648
来自专栏生信技能树

【直播】我的基因组56:探索遗传起源

首先,节日快乐!在这个众人狂欢的节日里,我冷静了冷静,听说知识量储备差太多的人做不了朋友,于是默默的搬起了板凳专心学习。 ? 昨天我们看了看千人基因组计划的公共...

2807
来自专栏大数据文摘

[可视化]数据原来这样美

2214
来自专栏机器人网

史上最全汽车图解,老司机也未必认得完 (收藏版)

这份史上最全汽车各部件图解,非常值得收藏!就算是老师傅,有很多部件的名字你肯定听说过但不一定都知道在哪个位置吧。 ? 打开发动机盖,就是这个样子了,4A13...

2758
来自专栏镁客网

脑细胞不够用?Gladstone研究所推两步法研制出脑细胞 | 黑科技

960
来自专栏云计算D1net

NSA在两个RSA加密产品中植入了后门

根据路透社的独家报道,NSA 在安全行业领导企业RSA的两个加密产品都植入了随机数生成器后门,而不只是此前斯诺登爆料的一个。 2013年12月路透社曾...

2523
来自专栏数说工作室

这个细胞是单身狗最后的福音,基因编辑让它“栩栩如生”,机器学习让它造福人类

干细胞——单身狗最后的福音 干细胞,被称为“万能细胞”,之前在 Science发布2016年度科学突破 里说过这个神奇的东西: 日本的科学家将小鼠的干细胞制作...

4048
来自专栏顶级程序员

全球顶尖公司的烧脑面试题,普通人一道都答不出来!

一些世界知名的企业在招聘时,可能会提供面试智力题,来筛选应聘者。那些越是大牌,越是有名的国内500强,乃至世界500强,给出的面试真是一般人都答不出来。 谷歌...

2577
来自专栏生信技能树

【直播我的基因组66:大多数性状往往是多个基因控制的

前面我们说到了那些简单的由单个基因决定的性状,这东西不需要预测,其中的生物学机制已经研究的非常透彻,只要拿到你的基因信息,很容易推断你的性状,比如人的乙醇脱氢酶...

3348

扫码关注云+社区