(57)【要約】【課題】配列間の変位部位や共通部位などに関する情報へのアクセスが容易で、検索された複数の配列をそのまま結果表示できる配列表示方法および相同性検索方法を提供する。【解決手段】本発明の配列表示方法は、類似した複数の塩基配列またはアミノ酸配列における変異性および／または共通性に応じた視覚的特徴を付加する。視覚的特徴は、変異部位および／または共通部位に付加され、変異性および／または共通性の程度、変異部位における変異の頻度に応じて付加される。さらに、視覚的特徴は、塩基配列をアミノ酸配列に変換し、変異部位のコドンに対応するアミノ酸名および性質を含むアミノ酸情報に基づいて付加される。類似した複数の塩基配列またはアミノ酸配列に関する情報へのリンクも設けられる。本発明の相同性検索方法は、この配列表示方法を用いて検索結果を表示する。【選択図】  図７
2004192257.tif 000002
【特許請求の範囲】
【請求項１】
類似した複数の塩基配列またはアミノ酸配列の表示方法において、
前記類似した複数の塩基配列またはアミノ酸配列における変異性および／または共通性に応じた視覚的特徴を付加することを特徴とする配列表示方法。
【請求項２】
前記視覚的特徴を前記類似した複数の塩基配列またはアミノ酸配列における変異部位および／または共通部位に付加する、請求項１に記載の配列表示方法。
【請求項３】
前記視覚的特徴を前記類似した複数の塩基配列またはアミノ酸配列における変異性および／または共通性の程度に応じて付加する、請求項１または２に記載の配列表示方法。
【請求項４】
前記視覚的特徴を前記類似した複数の塩基配列またはアミノ酸配列の前記変異部位における変異の頻度に応じて付加する、請求項２または３に記載の配列表示方法。
【請求項５】
前記視覚的特徴を、前記類似した複数の塩基配列におけるコドンに対応する、アミノ酸名および性質を含むアミノ酸情報に基づいて付加する、請求項２から４のいずれか１項に記載の配列表示方法。
【請求項６】
前記視覚的特徴を、表示色を用いて表現する、請求項１から５のいずれか１項に記載の配列表示方法。
【請求項７】
前記視覚的特徴を、表示色の変化を用いて表現する、請求項１から６のいずれか１項に記載の配列表示方法。
【請求項８】
前記視覚的特徴を、文字の種類を用いて表現する、請求項１から５のいずれか１項に記載の配列表示方法。
【請求項９】
前記類似した複数の塩基配列またはアミノ酸配列を、各配列中の対応する塩基またはアミノ酸を揃えて表示する、請求項１から８のいずれか１項に記載の配列表示方法。
【請求項１０】
前記類似した複数の塩基配列またはアミノ酸配列どうしのリンクおよび／または関連する情報へのリンクを付加する、請求項１から９のいずれか１項に記載の配列表示方法。
【請求項１１】
類似した複数の塩基配列またはアミノ酸配列の表示装置において、
前記類似した複数の塩基配列またはアミノ酸配列における変異性および／または共通性に応じた視覚的特徴を付加する手段を有することを特徴とする配列表示装置。
【請求項１２】
前記手段は、前記視覚的特徴を前記類似した複数の塩基配列またはアミノ酸配列における変異部位および／または共通部位に付加する、請求項１１に記載の配列表示装置。
【請求項１３】
前記手段は、前記視覚的特徴を前記類似した複数の塩基配列またはアミノ酸配列における変異性および／または共通性の程度に応じて付加する、請求項１１または１２に記載の配列表示装置。
【請求項１４】
前記手段は、前記視覚的特徴を前記類似した複数の塩基配列またはアミノ酸配列の前記変異部位における変異の頻度に応じて付加する、請求項１２または１３に記載の配列表示装置。
【請求項１５】
前記手段は、前記視覚的特徴を、前記類似した複数の塩基配列におけるコドンに対応する、アミノ酸名および性質を含むアミノ酸情報に基づいて付加する、請求項１２から１４のいずれか１項に記載の配列表示装置。
【請求項１６】
前記手段は、前記視覚的特徴を、表示色を用いて表現する、請求項１１から１５のいずれか１項に記載の配列表示装置。
【請求項１７】
前記手段は、前記視覚的特徴を、表示色の変化を用いて表現する、請求項１１から１６のいずれか１項に記載の配列表示装置。
【請求項１８】
前記手段は、前記視覚的特徴を、文字の種類を用いて表現する、請求項１１から１５のいずれか１項に記載の配列表示装置。
【請求項１９】
前記類似した複数の塩基配列またはアミノ酸配列を、各配列中の対応する塩基またはアミノ酸を揃えて表示する手段をさらに有する、請求項１１から１８のいずれか１項に記載の配列表示装置。
【請求項２０】
前記類似した複数の塩基配列またはアミノ酸配列どうしのリンクおよび／または関連する情報へのリンクを付加する手段をさらに有する、請求項１１から１９のいずれか１項に記載の配列表示装置。
【請求項２１】
請求項１から１０のいずれか１項に記載の配列表示方法をコンピュータに実行させる配列表示プログラム。
【請求項２２】
請求項１から１０のいずれか１項に記載の配列表示方法をコンピュータに実行させる配列表示プログラムを記録した、コンピュータで読み取りが可能な記録媒体。
【請求項２３】
ユーザによって入力された生物情報データベースへの問い合わせを解析するステップと、
該問い合わせの解析をもとに、前記生物情報データベースに適合する検索条件を生成して該生物情報データベースを検索するステップと、
該検索結果を解析するステップと、
該解析結果と請求項１から１０のいずれか１項に記載の配列表示方法を用いて、前記検索結果を表示するステップを有する相同性検索方法。
【請求項２４】
ユーザによって入力された生物情報データベースへの問い合わせを解析するクエリ解析手段と、
該問い合わせの解析をもとに、前記生物情報データベースに適合する検索条件を生成して該生物情報データベースを検索する通信手段と、
該検索結果を解析する結果解析手段と、
該解析結果と請求項１１から２０のいずれか１項に記載の配列表示装置を用いて、前記検索結果を表示する表示手段を有する相同性検索装置。
【請求項２５】
前記クエリ解析手段と前記通信手段と前記結果解析手段と前記表示手段の動作を制御する制御手段をさらに有する、請求項２４に記載の相同性検索装置。
【請求項２６】
請求項２３に記載の相同性検索方法をコンピュータに実行させる相同性検索プログラム。
【請求項２７】
請求項２３に記載の相同性検索方法をコンピュータに実行させる相同性検索プログラムを記録した、コンピュータで読み取りが可能な記録媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、塩基配列またはアミノ酸配列の相同性検索に関する。
【０００２】
【従来の技術】
現在、ヒトゲノム配列（ヒトのすべてのＤＮＡ配列）がほぼ解読され、アメリカ国立健康研究所ＮＩＨ（Ｎａｔｉｏｎａｌ  Ｉｎｓｔｉｔｕｔｅｓ  ｏｆ  Ｈｅａｌｔｈ）に所属する一機関であるアメリカ国立医学図書館ＮＬＭ（Ｎａｔｉｏｎａｌ  Ｌｉｂｒａｒｙ  ｏｆ  Ｍｅｄｉｃｉｎｅ）の生物工学情報センターＮＣＢＩ（Ｎａｔｉｏｎａｌ  Ｃｅｎｔｅｒ  ｆｏｒ  Ｂｉｏｔｅｃｈｎｏｌｏｇｙ  Ｉｎｆｏｍａｔｉｏｎ）、日本の国立遺伝学研究所、欧州の生物情報研究機関であるＥＢＩ（Ｅｕｒｏｐｅａｎ  Ｂｉｏｉｎｆｏｒｍａｔｉｃｓ  Ｉｎｓｔｉｔｕｔｅ）およびヨーロッパ分子生物学研究所ＥＭＢＬ（Ｅｕｒｏｐｅａｎ  Ｍｏｌｅｃｕｌａｒ  Ｂｉｏｌｏｇｙ  Ｌａｂｏｒａｔｏｒｙ）、日本の東京大学医科学研究所等によって、ＤＮＡ配列およびＤＮＡ配列に由来するアミノ酸配列などがデータベース（ＤＢ）として構築されている。
【０００３】
他の生物種についても、マウス、ラット、フグ、ゼブラフィッシュ、ショウジョウバエ、線虫などのゲノムが精力的に解読されつつあるか、すでに解読が完了している。
【０００４】
さらに次の段階として、ゲノム配列における個人差を調べる試みも進みつつある。
【０００５】
このような状況において、生物学および医学の分野では、研究者がその研究の過程で得た塩基配列またはアミノ酸配列と既知の配列との相同性比較を行うことが非常に重要である。なぜなら、他の研究者によって相同性の非常に高い配列がすでに登録されている場合、その配列に由来する生体機能がある程度解析済みであることを意味し、自身の研究の方向性を変更しなければならないからである。逆に、まだ登録されていない場合には、既知の配列から機能を予測することも必要となる。
【０００６】
ＢＬＡＳＴ（Ｂａｓｉｃ  Ｌｏｃａｌ  Ａｌｉｇｎｍｅｎｔ  Ｓｅａｒｃｈ  Ｔｏｏｌ、ｈｔｔｐ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／ｂｌａｓｔ／）は、ＮＣＢＩが開発した塩基配列およびアミノ酸配列の相同性検索システムである（例えば、非特許文献１参照）。これを用いると、ＮＣＢＩが提供する塩基配列およびアミノ酸配列をはじめとするＤＢであるＥｎｔｒｅｚ（ｈｔｔｐ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／ｅｎｔｒｅｚ／）などから、相同性の高い配列情報を取得することができる。このため、全世界の生物医学系の研究者により日常的に幅広く利用されている。
【０００７】
ＢＬＡＳＴ以外の塩基配列およびアミノ酸配列の相同性検索システムとしては、ＦＡＳＴＡ（Ｆａｓｔ  Ａｌｉｇｎｍｅｎｔ）がよく利用されている（例えば、非特許文献２参照）。配列ＤＢに関しても、ＮＣＢＩのＤＢ以外に、日本の国立遺伝学研究所が運営する日本ＤＮＡデータバンクＤＤＢＪ（ＤＮＡ  Ｄａｔａ  Ｂａｎｋ  ｏｆ  Ｊａｐａｎ）、ＥＢＩ／ＥＭＢＬが運営する塩基配列およびアミノ酸配列ＤＢなどがある。さらに、日本の東京大学医科学研究所のヒトの一塩基置換ＤＢであるＪＳＮＰｓなど、その目的に応じて様々な塩基配列およびアミノ酸配列ＤＢが構築され、各研究者の目的に応じて利用されている。
【０００８】
上記各種の配列情報ＤＢに対して、塩基配列またはアミノ酸配列の相同性検索を行う場合、研究者は主に上記のＢＬＡＳＴまたはＦＡＳＴＡシステムを用いることが多い。
【０００９】
さらに、相同性検索を用いた配列解析では、３つ以上の塩基配列またはアミノ酸配列について、相同性に基づいて対応する部分をそろえて並べる機能（以下、多重配列比較またはマルチプルアライメントと記す）が必要なことも多い。マルチプルアライメントを実行するソフトウェアとしては例えば、ＥＢＩ／ＥＭＢＬで開発されたＣｌｕｓｔａｌ  Ｗ、Ｃｌｕｓｔａｌ  Ｘなどがあり、その他にも日本国内において開発されたソフトウェアを含めて数多くの多重配列比較ソフトウェアが存在し、研究者に利用されている。
【００１０】
【非特許文献１】
ジャーナル・オブ・モリキュラー・バイオロジー、２１５（１９９０）４０３〜４１０ページ（Ｊｏｕｎａｌ  ｏｆ  Ｍｏｌｅｃｕｌａｒ  Ｂｉｏｌｏｇｙ、２１５（１９９０）ｐｐ４０３〜４１０）
【非特許文献２】
ドゥーリトル（Ｄｏｏｌｉｔｔｌｅ）編、メソッド・イン・エンザイモロジー、１８３（１９９０）６３〜９８ページ（Ｍｅｔｈｏｄｓ  ｉｎ  Ｅｎｚｙｍｏｌｏｇｙ、１８３（１９９０）ｐｐ６３〜９８）、アカデミックプレス（Ａｃａｄｅｍｉｃ  Ｐｒｅｓｓ）
【００１１】
【発明が解決しようとする課題】
塩基配列またはアミノ酸配列の相同性検索において、ユーザが入力した塩基配列またはアミノ酸配列、あるいは、ユーザが入力した情報により導かれた塩基配列またはアミノ酸配列が、ＤＢへの問い合わせ（以下、クエリと呼ぶ）となり、類似の配列が検索される。検索結果においては、相同性により対応付けられる部位で、クエリとした塩基配列またはアミノ酸配列と、検索結果の塩基配列またはアミノ酸配列との間で、塩基またはアミノ酸が異なる（以下、変異と呼ぶ）箇所および同じである（以下、共通と呼ぶ）箇所が重要な意味を持っている。ＢＬＡＳＴまたはＦＡＳＴＡシステムを用いる場合、以下のような問題点がある。
【００１２】
１．クエリとした配列と検索された配列が対で表示されるため、多数の配列間における変異の部位および頻度などの重要な意味を持つ情報を得るには、ユーザがすべての対の比較結果を相互参照しなければならない。
【００１３】
２．検索の結果は、相同性の高い部位だけが切り取られて部分的に表示されるので、全体の配列からみてどの部分が表示されているのかが非常にわかりにくい。
【００１４】
３．比較した部位の塩基またはアミノ酸が同じ場合は「｜」や「．」で表され、結果を一目で把握することが困難である。
【００１５】
４．検索結果が単純なテキスト形式でのみ表示されるため、変異部位を見落としやすい。
【００１６】
以上まとめると、ＢＬＡＳＴまたはＦＡＳＴＡシステムの検索結果表示では、ユーザが最も必要とする、配列間の変位部位または共通部位などの重要な情報へのアクセスが容易ではない点が問題となる。
【００１７】
上記問題点に対処するために、研究者はしばしば多重配列比較（マルチプルアライメント）を行う場合がある。しかしながら、多重配列比較を行う上記のソフトウェアを用いても、上記問題点は十分に解決されない。
【００１８】
すなわち、例えば上記のＣｌｕｓｔａｌ  Ｘは、アミノ酸の変異部位を色付き表示する機能もあるが、色付けの基準が一定ではなく、変異があることを強調する以上の意味を有していない。
【００１９】
さらに、用いるソフトウェアによっては相同性の求め方がＢＬＡＳＴやＦＡＳＴＡのそれとは異なっており、ＢＬＡＳＴやＦＡＳＴＡの検索結果に対して多重配列比較を行うと、しばしば塩基位置の対応がずれる。このような場合、例えばＢＬＡＳＴで得た検索結果を基準として塩基位置の対応を合わせる手作業が必要となる。
【００２０】
加えて、これらのソフトウェアを用いて相同性を求めるには、比較を行う全ての配列情報を予め収集しておく必要があり、収集すべき配列情報が多量にある場合、収集漏れが生じたり、漏れがないか検証する作業が研究者への大きな負担となっていたりした。
【００２１】
そこで本発明の目的は、配列間の変位部位や共通部位などに関する情報へのアクセスが容易で、なおかつ、検索された複数の配列がそのまま結果表示される配列表示方法および相同性検索方法を提供する。
【００２２】
【課題を解決するための手段】
上記目的を達成するために、本発明の配列表示方法は、類似した複数の塩基配列またはアミノ酸配列における変異性および／または共通性に応じた視覚的特徴を付加することを特徴としている。
【００２３】
視覚的特徴は、類似した複数の塩基配列またはアミノ酸配列における変異部位および／または共通部位に付加される。
【００２４】
さらに、視覚的特徴は、類似した複数の塩基配列またはアミノ酸配列における変異性および／または共通性の程度に応じて付加される。
【００２５】
さらに、視覚的特徴は、類似した複数の塩基配列またはアミノ酸配列の変異部位における変異の頻度に応じて付加される。
【００２６】
以上により、ユーザは、ユーザ自身による相互参照なしに、変異性および共通性などの重要な情報を即座に得ることができる。
【００２７】
さらに、視覚的特徴は、塩基配列をアミノ酸配列に変換して、塩基配列における変異部位のコドンに対応する、アミノ酸名および性質を含むアミノ酸情報に基づいて付加される。これにより、塩基配列の変異がアミノ酸生成に及ぼす影響を評価することが可能になる。
【００２８】
これらの視覚的特徴は、表示色、さらに表示色の変化、または文字の種類を用いて表現される。
【００２９】
また、類似した複数の塩基配列またはアミノ酸配列は、各配列中の対応する塩基またはアミノ酸を揃えて表示される。これにより、マルチプルアライメントの結果においても、塩基またはアミノ酸の変異および共通に関する部位、程度、頻度等をユーザが容易に認識することができる。
【００３０】
また、類似した複数の塩基配列またはアミノ酸配列に関する情報へのリンクを設ける。これにより、相同性検索結果が含む情報どうし、および、相同性検索結果が含む情報とインターネット上の情報のような外部の情報との関連付けも、表示情報が含むことができるようになり、相互参照や参照情報の閲覧も容易に行えるようになる。
【００３１】
本発明の相同性検索方法は、ユーザによって入力された生物情報データベースへの問い合わせを解析するステップと、この問い合わせの解析をもとに、生物情報データベースに適合する検索条件を生成して生物情報データベースを検索するステップと、検索結果を解析するステップと、解析結果と上述の配列表示方法を用いて、検索結果を表示するステップからなる。以上により、検索から表示が一体に行え、多重配列比較を行う場合にも、収集漏れが生じたり、漏れがないか検証したりする必要はない。
【００３２】
また、相同性検索装置に制御部を設け、各上記ステップの処理を制御させるようにしてもよい。
【００３３】
【発明の実施の形態】
本発明の実施の形態について、図面を参照して詳細に説明する。
【００３４】
図１を参照すると、本発明の第１の実施形態の相同性検索システム１０は、クエリ解析部１１と、通信部１２と、結果解析部１３と、表示部１４を有する。
【００３５】
クエリ解析部１１は、ユーザが入力した生物情報ＤＢのクエリ情報（クエリ）を受け取って検索条件の解析を行い、生物情報ＤＢの検索条件情報を通信部１２に出力する。入力されるクエリは、塩基配列またはアミノ酸配列が主な情報として含まれる。他に、生物学的情報としてアクセッション番号、実験データ、ファイルなどを情報として含めることができる。加えて、ＤＢの指定、結果を絞り込むための生物種、塩基配列ならばアミノ酸に翻訳するか、アミノ酸配列ならばＰＳＩ−ＢＬＡＳＴを行うかの指定の情報を含めることができる。
【００３６】
通信部１２は、クエリ解析部１１からの検索条件情報を受け取り、サーバ２１を介して配列情報ＤＢ２２に送り、配列情報ＤＢ２２からサーバ２１を介して検索結果を受け取って検索結果および検索条件情報を結果解析部１３に出力する。
【００３７】
結果解析部１３は、通信部１２から検索条件情報および検索結果を受け取り、検索結果の解析を行って、表示部１４に解析結果を出力する。検索処理によって生成された結果を解析し、例えば変異の部位、頻度、程度といった情報を得る。
【００３８】
表示部１４は、結果解析部１３から解析結果を受け、視覚的情報（例えば文字の色、文字の種類など）を付けた上で表示情報を生成し、出力する。結果解析部１３によって機能づけされた情報を、視覚的特性（例えば文字の色）を付ける、またはリンクを生成するなど行い、適切な表示形式に変換し、表示を行う。
【００３９】
図２を参照すると、クエリ解析部１１における処理のフローチャートが示されている。
【００４０】
ステップ２０１で、入力されたクエリに基づいて相同性検索および検索結果の解析をどのように行うかを決定する。クエリが配列データまたはアクセッション番号等であった場合、相同性検索および検索結果の解析を行うようにフラグを設定する。クエリが検索結果ファイルを指定するものであった場合、相同性検索を行わないようにフラグを設定する。
【００４１】
ステップ２０２で、入力されたクエリに基づき、クエリが含む情報から配列情報への変換を行うかどうかのフラグを設定する。例えば、クエリがアクセッション番号を含み、ユーザが配列情報への変換を行うことを指定したような場合、ステップ２０３で、配列情報変換フラグの値を真にする。
【００４２】
ステップ２０４で、入力されたクエリに基づき、入力された、または入力された生物情報により導かれた配列が、塩基配列かアミノ酸配列かを判定する。
【００４３】
ステップ２０４で塩基配列と判定された場合、ステップ２０５で、ユーザが比較すべき塩基配列ＤＢとして何を指定したかの情報を取得する。塩基配列ＤＢとしては例えば、ｎｒ（Ｎｏｎ−Ｒｅｄｕｎｄａｎｔ、標準となる非冗長配列ＤＢ）、ＥＳＴ（Ｅｘｐｒｅｓｓｅｄ  Ｓｅｑｕｅｎｃｅ  Ｔａｇ、発現している遺伝子の配列ＤＢ）、ＳＮＰｓ（Ｓｉｎｇｌｅ  Ｎｕｃｌｅｏｔｉｄｅ  Ｐｏｌｙｍｏｒｐｈｉｓｍｓ、一塩基多型配列ＤＢ）などである。
【００４４】
ステップ２０６で、入力されたクエリに基づき、入力されたクエリの、または入力されたクエリによって導かれる塩基配列をアミノ酸配列に変換して、さらなる解析を行うかどうかを判断する。変換する場合、ステップ２０７でアミノ酸変換フラグの値を真にする。
【００４５】
ステップ２０４で配列がアミノ酸配列だったと判定された場合、ステップ２０８で、ユーザが比較すべきアミノ酸配列ＤＢとして何を指定したのかの情報を取得する。アミノ酸配列ＤＢとしては例えば、ｎｒ（ｎｏｎ−Ｒｅｄｕｎｄａｎｔ、標準となる非冗長配列ＤＢ）、ＰＤＢ（Ｐｒｏｔｅｉｎ  Ｄａｔａ  Ｂａｎｋ、タンパク質立体構造ＤＢ）、ｓｗｉｓｓｐｒｏｔ（アミノ酸配列ＤＢ）、ｐａｔｅｎｔ（特許における配列ＤＢ）、ｙｅａｓｔ（酵母の配列ＤＢ）、ｍｏｎｔｈ（直近に追加された配列ＤＢ）などである。
【００４６】
ステップ２０９で、入力されたクエリに基づき、ＰＳＩ−ＢＬＡＳＴ（Ｐｏｓｉｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｔｅｒａｃｔｉｖｅ  Ｂａｓｉｃ  Ｌｏｃａｌ  Ａｌｉｇｎｍｅｎｔ  Ｓｅａｒｃｈ  Ｔｏｏｌ）を行うかどうかを判断する。行う場合、ステップ２１０でＰＳＩ−ＢＬＡＳＴフラグの値を真にする。
【００４７】
ステップ２１１で、入力されたクエリに基づき、検索結果を生物種によって絞り込むかどうかを判断する。ここで生物種とは例えば、ウイルス、バクテリア、菌類、真核生物、哺乳類、げっ歯類、シロイヌナズナ、枯草菌、線虫、ゼブラフィッシュ、キイロショウジョウバエ、大腸菌、ヒト、マウス、ラット、アフリカツメガエルなどである。ユーザが絞り込むと指定している場合、ステップ２１２で生物種情報フラグの値を真にする。
【００４８】
図３を参照すると、通信部１２における処理のフローチャートが示されている。
【００４９】
ステップ３０１で、クエリ解析部１１で解析されたクエリに関し、あらかじめ生物情報ＤＢの情報を保持しておき、どの生物情報ＤＢに問い合わせるかを決定する。もちろん、ユーザによってクエリを行う生物情報ＤＢの情報を設定するようにすることも容易に実現できる。
【００５０】
ＤＢの場所は、大きく分ければ、ローカルのみ、ローカルとインターネットの両方、インターネットのみである。ここでローカルとは、本実施形態を実装するコンピュータ自身およびＬＡＮ（ローカルエリアネットワーク）で接続された範囲を指す。インターネットとは、ＷＡＮ（ワイドエリアネットワーク、電話回線などを利用した接続）の範囲をさす。各々、複数のＤＢにクエリを振り分けることができる。
【００５１】
ローカルのみの場合、ステップ３０２で、クエリ解析部１１で解析された塩基配列またはアミノ酸配列の情報をＤＢに見合った形式で作成する。
【００５２】
ステップ３０３で、ステップ３０２で作成されたクエリをローカルのサーバに渡し、相同性検索を実行させる。
【００５３】
ステップ３０４で、ステップ３０３で実行させた相同性検索の結果を統合する。これは、ローカルのみ、ローカルとインターネット、インターネットのみのそれぞれで、複数のＤＢにクエリを渡し、実行させるからである。
【００５４】
インターネットとローカルの両方の場合、ステップ３０５で、クエリ解析部１１で解析された塩基配列またはアミノ酸配列の情報をＤＢに見合った形式で作成する。
【００５５】
ステップ３０６で、ステップ３０５で作成されたクエリをローカルのサーバに渡し、相同性検索を実行させる。ステップ３０７で、ステップ３０５で作成されたクエリをインターネット上のサーバに渡し、相同性検索を実行させる。
【００５６】
ステップ３０８で、ステップ３０６およびステップ３０７で実行させた相同性検索の結果を統合する。
【００５７】
インターネットのみの場合、ステップ３０９で、クエリ解析部１１で解析された塩基配列またはアミノ酸配列の情報をＤＢに見合った形式で作成する。
【００５８】
ステップ３１０で、ステップ３０９で作成されたクエリをインターネット上のサーバに渡し、相同性検索を実行させる。
【００５９】
ステップ３１１で、ステップ３１０で実行させた相同性検索の結果を統合する。
【００６０】
ステップ３１２で、ステップ３０４またはステップ３０８またはステップ３１１で統合された検索結果を保存する。
【００６１】
図４を参照すると、結果解析部１３における処理のフローチャートが示されている。
【００６２】
ステップ４０１で、通信部１２から受信した検索結果をとりこむ。もしも、相同性検索を行わない場合、ユーザによって指定されたファイルをとりこむ。
【００６３】
ステップ４０２で、ステップ４０１でとりこんだ検索結果を解析し、クエリの配列での変異部位を同定し、その情報を取得する。
【００６４】
ステップ４０３で、ステップ４０２で得られた変異部位の情報を統合し、その変異が比較した配列の中でどのくらいの頻度で起きているかを計算する。
【００６５】
ステップ４０４で、ステップ４０３で得られた変異部位において、変異の程度を計算する。クエリが塩基配列で、図２のステップ２０７で塩基配列をアミノ酸に翻訳して解析するフラグが設定された場合、入力された塩基配列と変異による塩基配列をそれぞれアミノ酸配列に変換し、これらの比較を行う。
【００６６】
ここでいう変異の程度とは例えば、塩基配列ならばアミノ酸置換を伴うか、または、アミノ酸配列ならば疎水性度や酸性／塩基性度などの性質がどの程度変わるのかをいう。アミノ酸配列またはアミノ酸に翻訳した塩基配列の場合、例えば、ＢＬＯＳＵＭ６２という変換行列により、その性質の近さを計算できる。ただし、これらはあくまでも例であって、酸性度、疎水性度、等電点によって計算するなど、これらに限定されるものではない。
【００６７】
ステップ４０５で、表示部１４に解析結果を出力する。
【００６８】
図４を参照すると、表示部１４における処理のフローチャートが示されている。
【００６９】
ステップ５０１では、図４のステップ４０２によって解析された変異部位の情報をもとに、対応する部位をそろえるよう並べなおす。
【００７０】
ステップ５０２では、図４のステップ４０２によって解析された変異部位の情報をもとに、変異部位の視覚的特性（例えば文字の色や文字の種類など）を決定する。
【００７１】
ここで視覚的特性とは例えば、文字についてフォント、大文字／小文字、色など、地について色、パターン、テクスチャ、アニメーションなどがある。ただし、これらはあくまでも例であって、文字のサイズ、太字、下線、斜体、点滅、影つき、白抜き、箱囲み、地のテクスチャ、アニメーションなども含まれる。視覚的特性に関しても同様である。
【００７２】
ステップ５０３で、図４のステップ４０３によって解析された変異の頻度の情報をもとに、変異の頻度の視覚的特性（例えば地の色）を決定する。ここでは、表示色の変化を用いて頻度を表現することが望ましい。ここでいう表示色の変化とは、例えば色の濃淡や色味の変化（一例として、赤から青へなど）などをさす。
【００７３】
ステップ５０４で、図４のステップ４０４によって解析された変異の程度の情報をもとに、変異の程度の視覚的特性（例えば地の色）を決定する。ここでは、表示色の変化を用いて頻度を表現することが望ましい。ここでいう表示色の変化とは、例えば色の濃淡や色味の変化（一例として、赤から青へなど）などをさす。
【００７４】
ステップ５０５で、相同性検索の結果が含む情報をリンクにより結ぶ。
【００７５】
ステップ５０６で、図５のステップ５０１からステップ５０５によってつけられた視覚的特性の情報をもとに、出力に合わせて適切な表示形式に変換する。ここでいう出力とは例えば、画面やプリンタなどがあり、表示形式とは例えば、ＨＴＭＬ（Ｈｙｐｅｒ  Ｔｅｘｔ  Ｍａｒｋｕｐ  Ｌａｎｇｕａｇｅ）、ＸＭＬ（ｅＸｔｅｎｓｉｂｌｅ  Ｍａｒｋｕｐ  Ｌａｎｇｕａｇｅ）などがある。ただし、これらはあくまでも例であって、ＴｅＸ、ｂｍｐ、ｇｉｆ、ｊｐｅｇ、ＰＮＧ、ＴＩＦＦ、ＰＩＣＴ、ＰＤＦ（Ｐｏｒｔａｂｌｅ  Ｄｏｃｕｍｅｎｔ  Ｆｏｒｍａｔ）、ＰｏｓｔＳｃｒｉｐｔを含む。
【００７６】
次に、第１の実施形態の具体的な実施例について、図面を参照して詳細に説明する。以下は、図２〜図５のフローチャートをコンピュータプログラムとして作成して記録媒体に記録し、例えばパーソナルコンピュータに読み取らせて動作させる場合を示す。
【００７７】
図６は、塩基配列について相同性検索を行う際の入力画面を示している。ここでは、ユーザによって、塩基配列の入力、検索を行うＤＢの指定、アミノ酸に翻訳して解析するかどうかの指定、検索結果を限定するための生物種の指定、まだ解析が行われていない検索結果ファイルの指定が行える。
【００７８】
ユーザが塩基配列入力フォームに塩基配列を入力した後、Ｓｕｂｍｉｔボタンをクリックすると、クエリ解析部１１が塩基配列を受け取り、続いて上述のような各部の処理が行われる。
【００７９】
本実施例では、ユーザが、アミノ酸に翻訳して解析するかどうかを指定できるようになっており、ここでの指定をアミノ酸変換フラグとして保持する。フラグが真ならば、アミノ酸変換手段により入力された塩基配列をアミノ酸配列に変換したうえで、変異の程度をアミノ酸配列も考慮に入れて計算する。これにより、ユーザは、その塩基の変異がアミノ酸置換を伴うかどうかを容易に把握し、実験操作によりその変異を元に戻さねばならないか否かの示唆が得られる。
【００８０】
また、本実施例では、比較するＤＢをユーザが指定できるようになっており、これにより、あらかじめ保持してある生物情報ＤＢの情報をもとに、クエリを振り分ける。例えば、ユーザがクエリとしてある塩基配列、および生物情報ＤＢとしてｎｒおよびＥＳＴを指定すると、この塩基配列をｎｒとＥＳＴ、それぞれで相同性検索を行うように検索条件情報を通信部１２に送る。これにより、ユーザは、いくつかのＤＢを指定することで必要な情報を網羅的に得られる一方、用途に合わせたＤＢを選ぶことで冗長な検索結果を排除することができ、検索結果を容易に把握することができる。
【００８１】
また、本実施例では、検索結果を生物種によって絞り込むかどうかをユーザが指定することができるようになっている。ここでの指定を生物種情報フラグに保持する。フラグが真ならば、特定の生物種の配列との比較結果のみを表示する。例えば、ユーザは生物種として、ヒト、マウス、ラットなどが指定できる。これにより、ユーザは、冗長な検索結果を排除することができ、検索結果を容易に把握することができる。
【００８２】
また、本実施例では、ユーザが従来の相同性検索によって得られた結果を解析対象として選択することができるようになっている。例えば、結果がまだ未解析の解析結果ファイルをユーザが指定できる。これにより、相同性検索を行わず、結果解析のみを行うようにフラグを設定する。そのためユーザは、本発明によらない検索結果に対しても、本発明を用いた解析によりこの検索結果を把握することができる。
【００８３】
図７は、塩基配列について相同性検索を行う際の結果画面の一部を示している。入力された塩基配列のうち、変異のある部位は赤字で表している。塩基配列の下に、対応するアミノ酸配列が記載されており、特にアミノ酸置換を伴う変異部位は小文字のアルファベットで示されている。これにより、その塩基の変異が、アミノ酸置換を伴うかどうかをユーザは容易に把握することができ、実験操作により、その変異を元に戻さねばならないかどうかをユーザは判定できる。
【００８４】
また、各々の変異部位に、図８にて後述する詳細な配列比較結果へのリンクの情報がつけられている。これにより、ユーザは入力した配列の全体的な変異の状況と各変異の詳細な状況を相互参照することが可能である。
【００８５】
図８は、塩基配列について相同性検索を行う際の結果画面の別の一部を示している。各々の配列は６０文字ごとに区切られている。文字数は、例えば９０文字、１２０文字、または表示できるだけの文字数であってもよい。本実施例では、塩基３つでアミノ酸１つに対応して翻訳されることを考慮して、３で割り切れ、なおかつ、ユーザへの分かりやすさから１０で割り切れる、すなわち３０の倍数の文字数ごとに区切っている。
【００８６】
また、検索結果の中に、変異の中でも、特に塩基の挿入が見られた部位では、地の色を黄色にして特に強調している。これは、塩基の挿入または欠落は、その部位以降のアミノ酸配列がまったく変わってしまうことを意味し、研究者にとって、特に重要な情報だからである。
【００８７】
また、各々の変異において、それがどのくらいの頻度でおきているかを地の色の濃淡で表現している。これにより、ユーザは、その変異が実験上、避けられない塩基置換であるか、または、配列特有の特徴（例えば、その変異が、ある疾患の原因となる変異である場合など）であるか、実験によりその変異を元に戻さなければならないかの示唆が得ることができる。
【００８８】
また、相同性に基づいて対応する部分をそろえて配列を表示している。既存の配列比較では、図１２に示すように、入力された配列と比較する配列を個別に１対１で比較していた。これにより、ユーザは入力された配列と検索された配列で網羅的に変異の情報について把握することができ、実験操作によりその変異を元に戻すかの示唆を得ることができる。
【００８９】
また、既存の配列比較結果である図１２にリンクの情報がつけられている。また、ユーザが配列の何塩基目を見たいかを指定すると、リンクによりその部位における配列比較を表示することができる。これにより、ユーザは変異情報を網羅的に相互参照することが可能である。
【００９０】
図９は、アミノ酸配列について相同性検索を行う際の入力画面を示している。ここでは、ユーザによって、アミノ酸配列の入力、検索を行うためのＤＢの指定、ＰＳＩ−ＢＬＡＳＴを行うかの指定、検索結果を限定するための生物種の指定、まだ解析の行われていない検索結果ファイルの指定が行える。
【００９１】
ユーザがアミノ酸配列入力フォームにアミノ酸配列を入力した後、Ｓｕｂｍｉｔボタンをクリックすると、クエリ解析部１１がアミノ酸配列を受け取り、上述の各部の処理が行われる。
【００９２】
本実施例では、ＰＳＩ−ＢＬＡＳＴを用いて、検索結果を得るかどうかをユーザが指定できるようになっており、ここでの指定を、ＰＳＩ−ＢＬＡＳＴフラグとして保持する。これにより、ユーザは検索結果を得る手段として、ＰＳＩ−ＢＬＡＳＴを使うことができ、ＰＳＩ−ＢＬＡＳＴでないと得られないような相同性をもつ配列について、比較することが可能である。
【００９３】
また、結果を生物種によって絞り込むかどうかをユーザが指定できるようになっており、ここでの指定を生物種情報フラグに保持する。
【００９４】
また、ユーザが従来の相同性検索によって得られた結果を解析対象として選択することが可能である。
【００９５】
図１０は、アミノ酸配列について相同性検索を行う際の結果画面を示している。入力されたアミノ酸配列のうち、変異のある部位は赤い文字で表している。これにより、ユーザは入力した配列について変異部位についての全体的な状況を容易に把握することができる。
【００９６】
また、比較した各々の配列は６０字ごとに区切られている。
【００９７】
また、相同性に基づいて対応する部分をそろえて配列を表示している。
【００９８】
また、各々の変異において、それがどの程度であるかを地の色の濃淡で表現している。
【００９９】
また、アミノ酸の挿入が見られた場合には、地の色を黄色にして表示している。
【０１００】
アミノ酸は、生物について２０種類が用いられているが、各々、疎水性度や酸性／塩基性度などの性質が決まっている。変異に伴って、これらの性質が変わると、配列全体としての性質が変わることが予想されるだけでなく、その変異箇所が酵素の反応中心の部位だった場合、酵素の機能が失われることとなる。それをふまえて、アミノ酸どうしの性質の違いを記述した変換行列を用いて、アミノ酸どうしの性質の近さを計算し、その値によって地の色の濃淡を変えて表現している。本実施例では、この変換行列として、図１１に示すＢＬＯＳＵＭ６２なる変換行列を用いている。
【０１０１】
また、図１０において既存の配列比較結果にリンクの情報がつけられている。また、ユーザが配列の何塩基目を見たいかを指定すると、リンクによりその部位を表示することができる。このことにより、ユーザは変異情報を網羅的に相互参照することが可能である。
【０１０２】
図１３を参照すると、本発明の第２の実施形態が示されている。本発明の第２の実施形態の相同性検索システムは、制御部１５を有する点が第１の実施の形態と異なっている。よって、以下では制御部１５に関する説明のみ行う。
【０１０３】
制御部１５は、クエリ解析部１１と、通信部１２と、結果解析部１３と、表示部１４の動作を制御するものである。制御部１５を有することにより、クエリを連続的に処理する、または、複数のＤＢを検索するなどの際、通信回線での混雑状況のモニター、または、解析結果の生成のモニターを行うなどによって、各部における処理のタイミングを制御することが可能となる。
【０１０４】
制御部１５は、クエリ解析部１１と、通信部１２と、結果解析部１３と、表示部１４から処理状況の通知を受け、その通知に従って各部の処理が順次行われるように制御する。よって、第２の実施形態におけるクエリ解析部１１と、通信部１２と、結果解析部１３と、表示部１４は、各々の処理の終了時に制御部１５に制御終了を通知する機能をさらに有し、制御部１５からの動作制御に従って動作する。クエリ解析部１１と、通信部１２と、結果解析部１３と、表示部１４における処理状況の通知出力および動作制御入力、制御部の処理をコンピュータプログラムとして作成し、第１の実施形態に付加することによって、第２の実施形態を実現することができる。これにより、例えばクエリが多数あって自動的に検索および解析を行いたい場合でも、ユーザがすべてのクエリを入力する必要がなくなり、ユーザの負担を大幅に削減することができる。
【０１０５】
【発明の効果】
以上説明したように、本発明によれば、次のような効果を有する。
【０１０６】
第１の効果は、塩基配列およびアミノ酸配列の相同性検索結果における変異の部位、頻度など重要な意味を持つ情報が、ユーザによる相互参照や再検索なしに得られることである。その理由は、複数の類似した塩基配列またはアミノ酸配列を含む相同性検索結果を入力とし、類似した複数の塩基配列またはアミノ酸配列における変異性と共通性に応じた視覚的な特徴を付加して表示情報として出力することができるからである。
【０１０７】
第２の効果は、相同性検索結果である複数の配列における変異部分が一目でわかることである。その理由は、類似する複数の塩基配列またはアミノ酸配列における変異部位に関する視覚的特徴を有する表示情報を出力できるからである。
【０１０８】
第３の効果は、相同性検索結果である複数の配列における共通部分が一目でわかることである。その理由は、類似する複数の塩基配列またはアミノ酸配列における共通部位に関する視覚的特徴を有する表示情報を出力できるからである。
【０１０９】
第４の効果は、相同性検索結果である複数の配列全体とともに相同性の高い部位を表示できることである。その理由は、本発明の表示情報が、配列全体を含むことができるからである。
【０１１０】
第５の効果は、マルチプルアライメントの結果における塩基位置のずれをなくして結果を表示できることである。その理由は、塩基およびアミノ酸の位置対応関係に従って整列された形式の表示情報を出力できるからである。
【０１１１】
第６の効果は、比較を行いたい配列情報を予め収集しておくことや、所望とする配列情報が手元に無い場合に、配列と関連する情報、例えば、遺伝子名、アミノ酸名、疾患名、タンパク質名、等を用いて配列情報を取得しておくことができることである。その理由は、複数の生物情報ＤＢに検索条件情報を送る手段を有するからである。
【０１１２】
第７の効果は、相同性検索結果およびマルチプルアライメントの結果である配列において、変異部分および共通部分が一目でわかるような表示を行えることである。その理由は、本発明が上記第１の効果から第６の効果を併せ持つことができるからである。
【０１１３】
第８の効果は、相同性検索結果を得た前後にユーザが行わなければならなかった作業を低減し、ユーザの負担を軽減することである。その理由は、マルチプルアライメントにおける位置合わせが不要であり、リンクを用いて相互参照も表示情報に包含可能であり、配列と関連する情報を用いて配列情報を取得可能であるからである。
【０１１４】
第９の効果は、検索条件を変えて何度も検索する場合でも、研究者が条件を変えては検索しなおすという作業を不要とし、効率的な相同性検索を可能にすることである。その理由は、検索の一連の手順を制御する制御部を有するからである。このような制御部を有することによって、公共の生物情報ＤＢへのクエリのタイミングを制御することもでき、公共の生物情報ＤＢサーバの負荷を軽減することにも有効である。
【図面の簡単な説明】
【図１】本発明の第１の実施形態の相同性検索システムの構成を示した図である。
【図２】クエリ解析部１１において実施される処理の流れを説明した図である。
【図３】通信部１２において実施される処理の流れを説明した図である。
【図４】結果解析部１３において実施される処理の流れを説明した図である。
【図５】表示部１４において実施される処理の流れを説明した図である。
【図６】塩基配列における相同性検索の入力画面の一例である。
【図７】塩基配列における相同性検索の結果画面の一例の一部である。
【図８】図７の塩基配列における相同性検索の結果画面の別の一部である。
【図９】アミノ酸配列における相同性検索の入力画面の一例である。
【図１０】アミノ酸配列における相同性検索の結果画面の一例である。
【図１１】アミノ酸置換行列ＢＬＯＳＵＭ６２である。
【図１２】遺伝子のクエリ配列に対する従来のＢＬＡＳＴ検索結果例である。
【図１３】本発明の第２の実施形態の相同性検索システムの構成を示した図である。
【符号の説明】
１０  相同性検索システム
１１  クエリ解析部
１２  通信部
１３  結果解析部
１４  表示部
１５  制御部
２１  サーバ
２２  データベース
２０１〜２１２  ステップ
３０１〜３１２  ステップ
４０１〜４０５  ステップ
５０１〜５０６  ステップ
【図１】
2004192257.tif 000003
【図２】
2004192257.tif 000004
【図３】
2004192257.tif 000005
【図４】
2004192257.tif 000006
【図５】
2004192257.tif 000007
【図６】
2004192257.tif 000008
【図７】
2004192257.tif 000009
【図８】
2004192257.tif 000010
【図９】
2004192257.tif 000011
【図１０】
2004192257.tif 000012
【図１１】
2004192257.tif 000013
【図１２】
2004192257.tif 000014
【図１３】
2004192257.tif 000015