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

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

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

相关文章

来自专栏C# 编程

C#使用DataSet类、DataTable类、DataRow类、OleDbConnection类、OleDbDataAdapter类编写简单数据库应用

//注意:请使用VS2010打开以下的源代码。 //源代码地址:http://pan.baidu.com/s/1j9WVR using System; usi...

2460
来自专栏码匠的流水账

聊聊spring cloud的AbstractLoadBalancingClient

本文主要研究一下spring cloud的AbstractLoadBalancingClient

832
来自专栏葡萄城控件技术团队

Table-values parameter(TVP)系列之二: 利用DataTable将其作为参数传给SP

一,回顾         上一部分讲述了“在T-SQL中创建和使用TVP”,通过T-SQL建立如下的对象:         1)Tables ...

2069
来自专栏闻道于事

商城项目整理(三)JDBC增删改查

商品表的增加,修改,删除,订单表的增加,确认,用户表的查看,日志表的增加,查看 商品表建表语句: 1 create table TEST.GOODS_TABL...

5455
来自专栏听雨堂

想修改CSS

      下载了一个“通用”的CSS文件,本来想偷懒的,结果发现有问题,就是它用的颜色是变量定义的,无法识别。我又找不到在哪里可以定义。 BODY{     ...

20410
来自专栏码匠的流水账

聊聊EurekaRibbonClientConfiguration

spring-cloud-netflix-eureka-client-2.0.0.RELEASE-sources.jar!/org/springframewor...

1171
来自专栏DT乱“码”

简单的考勤系统

连接数据库类 package com.lianrui.it; import java.sql.Connection; import java.sql.Driv...

3489
来自专栏跟着阿笨一起玩NET

[C#]工具类—FTP上传下载

  不错的文章:http://www.cnblogs.com/greatverve/archive/2012/03/03/csharp-ftp.html

1181
来自专栏c#开发者

jquery easyui datagrid mvc server端分页排序筛选的实现

1自定义一个ModelBinder public class filterRule { public string field { g...

4189
来自专栏积累沉淀

Hive2.0.0操作HBase 1.2.1报错解决

首先看错  org.apache.hive.service.cli.HiveSQLException: Failed to open new session: ...

2339

扫码关注云+社区