BioJava – 开源的生物信息学基础库

类别:Java 点击:0 评论:0 推荐:

转自醒客视点 autoasm.blog-city.com

 

BioJava – 开源的生物信息学基础库

 

 

Alex Dou

[email protected]

 

 

什么是生物信息学?

 

生物信息学,这确实是个很酷的名字。

从字面上理解,这是一门与生命科学和信息科学这两个当今热点领域都相关的学科。

那到底什么是生物信息学?在这里,我只能很遗憾地告诉大家,对于这样一门新兴的、而且正在不断发展和变化中的学科,我很难给出一个精确的定义。

 

广义地说,生物信息学从事对基因组研究相关生物信息的获取、加工、储存、分配、分析和解释。这一定义包括了两层含义,一是对海量数据的收集、整理与服务,也就是管好这些数据;另一个是从中发现新的规律,也就是用好这些数据。

 

具体地说,生物信息学是把基因组 DNA序列信息分析作为源头,找到基因组序列中代表蛋白质和 RNA基因的编码区;同时,阐明基因组中大量存在的非编码区的信息实质,破译隐藏在 DNA序列中的遗传语言规律;在此基础上,归纳、整理与基因组遗传信息释放及其调控相关的转录谱和蛋白质谱的数据,从而认识代谢、发育、分化、进化的规律。

 

信息技术的发展和应用使得几乎所有的人受益,研究DNA、RNA和蛋白质的分子生物学家也不例外。很难相信,离开信息技术,分子生物学家可以完成人类基因组的测序(事实上,大规模测序中所使用的鸟枪法确实依赖于计算技术),如果没有信息技术和计算分子生物学的理论支持,研究SARS的病毒学家也不可能在短时间内判断SARS病毒的种类归属(除非他们火眼精精能够从数以百万计的AGCT这几个符号的杂乱无章的排列中找出规律)。

 

生物信息学的本质即使用各种信息技术来支持生命科学的研究与开发。

 

 

BioJava简介

生物信息学在理论和工程上面临许多挑战。而要开发出复杂的生物序列分析系统,需要一些基础库的支持,BioJava就是这样的一套基础库。

BioJava是使用Java语言开发的、用于分析和表示生物序列(如DNA、RNA和蛋白质)的基础库。BioJava提供了生物序列处理功能(如转录与翻译)、文件格式转换功能和一些简单的科学计算(如隐马尔科夫模型)。

 

读者可以从http://www.biojava.org 获得更多的BioJava的介绍。另外,感谢北京大学生物信息学中心的吴欣(音译),他将BioJava的入门文档翻译成了中文。

 

 

转录—一个简单的例子

 

绝大多数生物都使用DNA来表达遗传信息(诸如SARS病毒和艾兹病毒使用RNA表示遗传信息)。但直接指导合成蛋白质的确是信使RNA。在分子生物学上,复制过程中从DNA到RNA的过程被称作转录(类似的,从RNA到DNA的过程叫做逆转录,逆转录过程通常出现在以SARS病毒和艾兹病毒为代表的RNA病毒复制过程中),如图1。

 

图1.转录与翻译

 

下面这段代码使用BioJava库来获得一段DNA序列所对应的RNA序列。

 

 

import org.biojava.bio.symbol.*;

import org.biojava.bio.seq.*;

 

public class TranscribeDNAtoRNA {

   

    public static void main(String[] args) {

    try {

            //make a DNA SymbolList

            SymbolList symL = DNATools.createDNA("atgccgaatcgtaa");

 

            //transcribe it to RNA

            symL = RNATools.transcribe(symL);

 

            //just to prove it worked

            System.out.println(symL.seqString());

    }

    catch (IllegalSymbolException ex) {

            //this will happen if you try and make the DNA seq using non          IUB symbols

            ex.printStackTrace();

    }

    catch (IllegalAlphabetException ex) {

            //this will happen if you try and transcribe a non DNA SymbolList

            ex.printStackTrace();

    }

    }

}

 

上面代码中,我们首先创建了一个DNA序列对象,该序列的内容为atgccgaatcgtaa。然后调用RNATools.transcribe()方法获得转录后的RNA序列。由于转录过程中DNA序列与RNA序列是一一对应的,且两个字符集中的符号映射关系也是确定的,所以这个方法的实现其实十分简单。当然,BioJava中也提供了诸如HMMER的更复杂的功能,有兴趣的读者可以参考BioJava网站。

 

 

相关阅读材料

如果你是一名软件工程师并且对生物信息学有兴趣的话,你可能需要了解更多的分子遗传学方面的知识。科学出版社的《基因组》是一个不错的选择,影印版的《分子生物学精要》也相当不错,前提是你要做好背专业词汇的准备。全美经典的《分子生物学和细胞生物学》则是不错的参考书。

 

如果你有生命科学领域的背景,你或许需要掌握更多的计算机科学方面的技能。《生物信息学中的计算机技术》是很好的入门材料。

 

另外,《计算分子生物学导论》是本相当经典的理论书籍,不过该书的读者需要有一定的计算机科学理论知识。

 

本文地址:http://com.8s8s.com/it/it16482.htm