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

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

最近!一个偶然的机会,发现了一个男性福音,咳咳,不(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 条评论
登录 后参与评论

相关文章

来自专栏Data Analysis & Viz

手把手教你完成一个数据科学小项目(8):Emoji提取与分布图谱

请先阅读“中国年轻人正带领国家走向危机”,这锅背是不背? 一文,以对“手把手教你完成一个数据科学小项目”系列有个全局性的了解。

2012
来自专栏FreeBuf

不懂DDoS攻击,不敢去卖煎饼

事件回放 据网友@大咕咕咕鸡爆料: “楼下的煎饼摊子o2o了!只接受微信在线交流,现场拒绝讲话,用一个纸箱子把摊子包起来了!就露一个小洞。买煎饼的在微信订好,提...

2257
来自专栏机器人网

伺服电机,步进电机,同步电动机和异步电动机的区别

同步电动机与异步电动机 原理: 同步电机就是靠励磁电流运行的,如果没有励磁,电机就是异步的。励磁是加在转子上的直流系统,它的旋转速度和极性与定子是一致的...

5847
来自专栏机器人网

什么是伺服电机,伺服电机知识汇总

“伺服”一词源于希腊语“奴隶”的意思。“伺服电机”可以理解为绝对服从控制信号指挥的电机:在控制信号发出之前,转子静止不动;当控制信号发出时,转子立即转动;当控制...

36610
来自专栏java一日一条

每个程序员和设计师必做的10项运动

防止肌肉和颈部麻痹的基本伸展动作。这套舒展动作不但有助于保持体型,还能修复突出的颈椎。

552
来自专栏程序员互动联盟

【程序员故事】三步把U盘变成MP3

程序员的故事,我们的故事。 1、下一关的出口在哪儿 一个电游玩家死后进入了地狱。一个星期之后,撒旦气急败坏地跑来问上帝:“你上周给我送过去的到底是个什么人?”...

2614
来自专栏Java社区

编程江湖的那些事

971
来自专栏机器人网

动图揭秘!身边最平常的机械,它们的工作原理你知道吗?

你一定好奇过自动贩卖机是如何识别硬币的、卷尺是如何把几米长的金属卷起来的?是时候进入他们的内部了,这些原理gif将完美解释你心中的疑问。同时,你会发现这些普通的...

3705
来自专栏小狼的世界

妈呀,狼来了?

听听这充满诱惑的广告语、感受这不断被煽动的情绪,您有何感想呢?目前我们确实已经有了facebook、kaixin001、校内这样的社交网络,我们也有了twitt...

792
来自专栏量子位

视频 | 机器人YuMi双臂齐飞、毛笔手书:祝大家“新年快乐”

下面这段视频,是ABB旗下的双臂机器人YuMi,双臂齐飞、使用毛笔手书“新年快乐”,为大家送上新春的祝福! 视频内容 简单对YuMi介绍两句: - 单臂7轴冗余...

3565

扫码关注云+社区