你的基因名被EXCEL自动转换成数字了吗?

1.充满恶意的Excel自动转换

做基因组学分析时发现,在默认设置下将数据导入Microsoft Excel后,部分基因名会自动转换为日期或浮点1。2016年发表于Genome Biology上的一项研究显示,补充材料里包含Excel格式基因名的文献中,约有20%出现了这种错误的转换2

The spreadsheet software Microsoft Excel, when used with default settings, is known to convert gene names to dates and floating-point numbers. A programmatic scan of leading genomics journals reveals that approximately one-fifth of papers with supplementary Excel gene lists contain erroneous gene name conversions.

举例来说,将下表中第一列的数据导入Excel后,将自动转换为第二列的内容:

源格式 导入Excel后的格式
MARC1 1-Mar
MARCH1 1-Mar
DEC1 1-Dec
231012E4 2.31E+09
123E1 1.23E+03

*请注意,基因MARC1和基因MARCH1均被转换成了1-Mar。

一旦导入Excel后,这种转换是不可逆的,也就是说无法通过将单元格格式设置为文本格式修复

如果我们尝试修改单元格格式为文本,将出现下图第三列所示的内容:

源格式 导入Excel后的格式 导入后再将单元格设置为文本格式
MARC1 1-Mar 43525
MARCH1 1-Mar 43525
DEC1 1-Dec 43800
231012E4 2.31E+09 2.31E+09
123E1 1.23E+03 1.23E+03

看到43525,43800这种数字,是不是感觉很眼熟?我相信很多人的基因Excel表里都有他们的存在。那么这个数字是怎么来的呢?
Excel的日期默认是从1990/1/1开始的,也就是说1900/1/1转化为数字格式将是1。而笔者写本文的时间是2019年,所以MARC1基因转换为1-Mar,也就是2019/3/1,它的数字格式为2019/3/1减去1900/1/1,即43525。

这一问题困扰了非常多的的科研人员,在Biostar上相关问题的点击已经超过3.5K3

2.如何避免这种转换?

最简单的方法,新建Excel表,单元格全选,更改为文本格式,然后再导入数据。但这种方式不利于数字计算。

此外,Eric A. Welsh开发了一个小工具escape_excel以阻止这种转换,其原理十分简单,就是在原始数据前后加上引号。这项工作还在2017年发表到了Plos One上4,5。(不得不说,想发表文章,idea非常重要)

该工具提供在线版本,Windows系统下的Excel插件,Perl脚本和Galaxy访问接口,教程请参见作者Github6

3.已经转换的数据如何处理?

建议从测序原始数据中提取正确的基因名,毕竟测序公司输出的文件一般都不是Excel格式。
如果找不到原始数据,也没有NCBI_GENE_ID或其他数据库基因ID,我们还是能把转换成日期或数字的基因名尽量还原。

日期格式的基因名:

以1-Mar为例,由于Excel日期转换只针对3个或以上字母的匹配,所以它的原始基因名可能为MAR1或MARC1或MARCH1,检索得知不存在MAR1基因,所以范围缩小到MARC1基因或MARCH1基因。

数字格式的基因名:

以43525为例,先将单元格转化为日期格式2019/3/1,忽略年份得到1-Mar,后续步骤同上。
因此,43160转化后的结果也为1-Mar,这与转换发生时的年份有关,但不影响还原方式。

以上就是基因名被EXCEL自动转换的解决方案了。最后留一个小建议,在操作每一个文件前,请先留好备份,以免发生意外情况。

参考资料:

1. Zeeberg B R , Riss J , Kane D W , et al. Mistaken Identifiers: Gene name errors can be introduced inadvertently when using Excel in bioinformatics[J]. BMC Bioinformatics, 2004, 5(1):80-0.
2. Ziemann M , Eren Y , El-Osta A . Gene name errors are widespread in the scientific literature[J]. Genome Biology, 2016, 17(1):177.
3. Question: How to avoid conversion of gene symbols to date format in Excel
4. Welsh E A , Stewart P A , Kuenzi B M , et al. Escape Excel: A tool for preventing gene symbol and accession conversion errors[J]. PLOS ONE, 2017, 12(9):e0185207-.
5. Escape gene name-mangling with ‘Escape Excel’
6. escape_excel插件的github地址

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 212,185评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,445评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 157,684评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,564评论 1 284
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,681评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 49,874评论 1 290
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,025评论 3 408
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,761评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,217评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,545评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,694评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,351评论 4 332
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,988评论 3 315
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,778评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,007评论 1 266
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,427评论 2 360
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,580评论 2 349

推荐阅读更多精彩内容