(57)【要約】【課題】新規な機能性配列選択方法及び機能性配列選択システムを提供すること。【解決手段】  遺伝子を細胞で発現させるための組換え遺伝子を作製するための機能性配列選択方法であって、データベースにおいて、目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、相同性を有する塩基配列を検索する工程と、得られた塩基配列のうち、全体がゲノム由来である塩基配列以外の塩基配列を抽出し、（１）抽出した塩基配列に対し、上流または下流に塩基配列が存在する場合、機能性配列かどうか解析して第１の機能性配列を抽出し、（２）存在しない場合、遺伝子情報の記載に、機能性配列を示す記載があるかどうかを検索して第２の機能性配列を抽出する工程と、得られた機能性配列に対してスコアリングを行う工程と、スコアリングの結果の上位の機能性配列を選択する工程と、を含む機能性配列選択方法とする。【選択図】図１
2020150919.tif 000002
【特許請求の範囲】
【請求項１】
遺伝子もしくはその一部の配列、遺伝子がコードするアミノ酸の配列もしくはその一部の配列、および／または前記遺伝子の遺伝子情報を含むデータを含むデータベースを用いた、目的遺伝子を細胞で発現させるための組換え遺伝子を作製するための機能性配列選択方法であって、
前記目的遺伝子のコーディング領域の塩基配列もしくはアミノ酸配列またはそれらの一部を含む配列をクェリとして入力するための入力装置と、
機能性配列を選択するための選択装置と、
選択した機能性配列を出力するための出力装置と、
を備えた機能性配列選択システムにおいて、
前記入力装置が前記クェリを入力し、
前記選択装置において、
相同性配列検索手段が、前記データベースにおいて、前記目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、当該塩基配列または当該アミノ酸配列をコードする塩基配列と相同性を有する塩基配列を検索し、
第１の配列リスト記憶手段が、検索して得られた塩基配列を第１の配列リストとして記憶し、
機能性配列抽出手段が、第１の配列リストから、塩基配列全体がゲノム由来である塩基配列を除去して、塩基配列全体がゲノム由来である塩基配列以外の塩基配列を抽出し、（１）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流または下流に塩基配列が存在する場合、上流または下流に存在する該塩基配列が機能性配列かどうか解析して１つまたはそれ以上の第１の機能性配列を抽出し、
（２）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流かつ下流に塩基配列が存在しない場合、前記データ内の前記遺伝子情報に、機能性配列を示す記載があるかどうかを検索して１つまたはそれ以上の第２の機能性配列を抽出し、
第２の配列リスト記憶手段が、第１の機能性配列及び第２の機能性配列を含む第２の配列リストを記憶し、
スコアリング手段が、第２の配列リストに対し、スコアリングを行い、
機能性配列選択手段が、前記スコアリングの結果の上位の機能性配列を選択し、
機能性配列記憶手段が、選択された前記上位の機能性配列を記憶し、
前記出力装置が、選択された前記上位の機能性配列を出力する、
機能性配列選択方法。
【請求項２】
前記遺伝子情報に、ゲノム由来である塩基配列を示すキーワードがあるかどうか検索することによって、全体がゲノム由来である塩基配列以外の塩基配列が抽出される、請求項１に記載の機能性配列選択方法。
【請求項３】
前記機能性配列は、プロモーター配列、ターミネーター配列、及び終止コドンを含む、請求項１または２に記載の機能性配列選択方法。
【請求項４】
（１）の場合において、機能性配列ライブラリーまたは前記遺伝子情報を用いて、上流または下流に存在する該塩基配列が機能性配列かどうか解析される、請求項１〜３に記載の機能性配列選択方法。
【請求項５】
第１の機能性配列及び第２の機能性配列に対し、出現頻度を基準として、または各機能性配列が由来する前記データ内の前記遺伝子情報に出現する特定のキーワードの頻度を基
準としてスコアリングが行われる、請求項１〜４に記載の機能性配列選択方法。
【請求項６】
第１の機能性配列及び第２の機能性配列に対し、スコアリングを行うまえに、
第１の機能性配列及び第２の機能性配列から、前記目的遺伝子が天然に有する機能性配列を除外する工程を有する、請求項１〜５に記載の機能性配列選択方法。
【請求項７】
遺伝子もしくはその一部の配列、遺伝子がコードするアミノ酸の配列もしくはその一部の配列、および／または前記遺伝子の遺伝子情報を含むデータを含むデータベースを用いた、目的遺伝子を細胞で発現させるための組換え遺伝子を作製するための機能性配列選択システムであって、
前記目的遺伝子のコーディング領域の塩基配列もしくはアミノ酸配列またはそれらの一部を含む配列をクェリとして入力するための入力装置と、
機能性配列を選択するための選択装置と、
選択した機能性配列を出力するための出力装置と、
を備え、前記選択装置は、
前記データベースにおいて、前記目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、当該塩基配列または当該アミノ酸配列をコードする塩基配列と相同性を有する塩基配列を検索するための相同性配列検索手段と、
検索して得られた塩基配列を第１の配列リストとして格納するための第１の配列リスト記憶手段と、
第１の配列リストから、塩基配列全体がゲノム由来である塩基配列を除去して、塩基配列全体がゲノム由来である塩基配列以外の塩基配列を抽出し、
（１）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流または下流に塩基配列が存在する場合、上流または下流に存在する該塩基配列が機能性配列かどうか解析して１つまたはそれ以上の第１の機能性配列を抽出し、
（２）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流かつ下流に塩基配列が存在しない場合、前記データ内の前記遺伝子情報に、機能性配列を示す記載があるかどうかを検索して１つまたはそれ以上の第２の機能性配列を抽出するための機能性配列抽出手段と、
第１の機能性配列及び第２の機能性配列を含む第２の配列リストを格納するための第２の配列リスト記憶手段と、
第２の配列リストに対し、スコアリングを行うためのスコアリング手段と、
前記スコアリングの結果の上位の機能性配列を選択するための機能性配列選択手段と、
選択された前記上位の機能性配列を格納するための機能性配列記憶手段と、
を備える機能性配列選択システム。
【請求項８】
前記遺伝子情報に、ゲノム由来である塩基配列を示すキーワードがあるかどうか検索することによって、全体がゲノム由来である塩基配列以外の塩基配列が抽出される、請求項９に記載の機能性配列選択システム。
【請求項９】
前記機能性配列は、プロモーター配列、ターミネーター配列、及び終止コドンを含む、請求項７または８に記載の機能性配列選択システム。
【請求項１０】
（１）の場合において、機能性配列ライブラリーまたは前記遺伝子情報を用いて、上流または下流に存在する該塩基配列が機能性配列かどうか解析される、請求項７〜９に記載の機能性配列選択システム。
【請求項１１】
第１の機能性配列及び第２の機能性配列に対し、出現頻度を基準として、または各機能性配列が由来する前記遺伝子情報に出現する特定のキーワードの頻度を基準としてスコアリングが行われる、請求項７〜１０に記載の機能性配列選択システム。
【請求項１２】
第１の機能性配列及び第２の機能性配列に対し、スコアリングを行う前に、
第１の機能性配列及び第２の機能性配列から、前記目的遺伝子が天然に有する機能性配列が除外される、請求項７〜１１に記載の機能性配列選択システム。
【請求項１３】
遺伝子もしくはその一部の配列、遺伝子がコードするアミノ酸の配列もしくはその一部の配列、および／または前記遺伝子の遺伝子情報を含むデータを含むデータベースを用いた、目的遺伝子を細胞で発現させるための組換え遺伝子を作製するための機能性配列選択方法であって、
前記データベースにおいて、前記目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、当該塩基配列または当該アミノ酸配列をコードする塩基配列と相同性を有する塩基配列を検索する工程と、
検索して得られた塩基配列またはアミノ酸配列をコードする塩基配列のうち、全体がゲノム由来である塩基配列を除去して、全体がゲノム由来である塩基配列以外の塩基配列を抽出し、
（１）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流または下流に塩基配列が存在する場合、上流または下流に存在する該塩基配列が機能性配列かどうか解析して１つまたはそれ以上の第１の機能性配列を抽出し、
（２）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流かつ下流に塩基配列が存在しない場合、前記遺伝子情報の記載に、機能性配列を示す記載があるかどうかを検索して１つまたはそれ以上の第２の機能性配列を抽出する工程と、
第１の機能性配列及び第２の機能性配列に対し、スコアリングを行う工程と、
前記スコアリングの結果の上位の機能性配列を選択する工程と、
を含む機能性配列選択方法。
【請求項１４】
請求項１〜６または１３に記載の機能性配列選択方法を前記機能性配列選択システムに実行させるためのプログラム。
【請求項１５】
請求項１４のプログラムを格納したコンピュータ可読記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、機能性配列選択方法及び機能性配列選択システムに関する。
【背景技術】
【０００２】
細胞に遺伝子を導入し、細胞内で発現させるために必要な機能性配列であるプロモーターやターミネーター配列は、導入する遺伝子配列や細胞の由来である生物種により異なっているため、最適な配列を選択するためには試行錯誤が必要となる。このような課題を解決するために、遺伝子設計、特に至適機能性配列選択を行なうことが容易となる遺伝子配列設計ツールが必要とされており、現在複数開発されている。
【０００３】
例えば、目的の遺伝子名を検索すると、米国のNCBI（ National Center For Biotechnology Information）とiGEM（The International Genetically Engineered Machine competition）のデータベースの情報のうち、登録されている文献及びプロモーター、ターミネーター配列を抽出することができるツールが開発されている（特許文献１参照）。このツールでは、ユーザーは、抽出した情報の中から、目的に近い文献を自身で選択することで、その文献で取り扱われたプロモーターとターミネーター配列を抽出することができる。このツールにより、登録文献から、導入遺伝子配列に適する機能性配列情報を得ることができるが、該当配列を選択するためには研究者が文献を解読し判断することが必要となる。また、目的に則した機能性配列を自動で選択できるツールも開発されている（特許文献２参照）。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Bates M et al., ACS Synth Biol., 6,12. (2017)
【非特許文献２】Nealsen A A K et al., Science,352, 6281.(2016)
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明の目的は、新規な機能性配列選択方法及び機能性配列選択システムを提供することである。
【課題を解決するための手段】
【０００６】
本発明は以下の態様を包含する。
【０００７】
本発明の一実施態様は、遺伝子もしくはその一部の配列、遺伝子がコードするアミノ酸の配列もしくはその一部の配列、および／または前記遺伝子の遺伝子情報を含むデータを含むデータベースを用いた、目的遺伝子を細胞で発現させるための組換え遺伝子を作製するための機能性配列選択方法であって、前記目的遺伝子のコーディング領域の塩基配列もしくはアミノ酸配列またはそれらの一部を含む配列をクェリとして入力するための入力装置と、機能性配列を選択するための選択装置と、選択した機能性配列を出力するための出力装置と、を備えた機能性配列選択システムにおいて、前記入力装置が前記クェリを入力し、前記選択装置において、相同性配列検索手段が、前記データベースにおいて、前記目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、当該塩基配列または当該アミノ酸配列をコードする塩基配列と相同性を有する塩基配列を検索し、第１の配列リスト記憶手段が、検索して得られた塩基配列を第１の配列リストとして記憶し、機能性配列抽出手段が、第１の配列リストから、塩基配列全体がゲノム由来である塩基配列を除去して、塩基配列全体がゲノム由来である塩基
配列以外の塩基配列を抽出し、
（１）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流または下流に塩基配列が存在する場合、上流または下流に存在する該塩基配列が機能性配列かどうか解析して１つまたはそれ以上の第１の機能性配列を抽出し、
（２）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流かつ下流に塩基配列が存在しない場合、前記データ内の前記遺伝子情報に、機能性配列を示す記載があるかどうかを検索して１つまたはそれ以上の第２の機能性配列を抽出し、
第２の配列リスト記憶手段が、第１の機能性配列及び第２の機能性配列を含む第２の配列リストを記憶し、スコアリング手段が、第２の配列リストに対し、スコアリングを行い、機能性配列選択手段が、前記スコアリングの結果の上位の機能性配列を選択し、機能性配列記憶手段が、選択された前記上位の機能性配列を記憶し、前記出力装置が、選択された前記上位の機能性配列を出力する、機能性配列選択方法である。前記遺伝子情報に、ゲノム由来である塩基配列を示すキーワードがあるかどうか検索することによって、全体がゲノム由来である塩基配列以外の塩基配列が抽出されてもよい。前記機能性配列は、プロモーター配列、ターミネーター配列、及び終止コドンを含んでもよい。（１）の場合において、機能性配列ライブラリーまたは前記遺伝子情報を用いて、上流または下流に存在する該塩基配列が機能性配列かどうか解析されてもよい。第１の機能性配列及び第２の機能性配列に対し、出現頻度を基準として、または各機能性配列が由来する前記データ内の前記遺伝子情報に出現する特定のキーワードの頻度を基準としてスコアリングが行われてもよい。第１の機能性配列及び第２の機能性配列に対し、スコアリングを行うまえに、第１の機能性配列及び第２の機能性配列から、前記目的遺伝子が天然に有する機能性配列を除外する工程を有してもよい。
【０００８】
本発明の他の実施態様は、遺伝子もしくはその一部の配列、遺伝子がコードするアミノ酸の配列もしくはその一部の配列、および／または前記遺伝子の遺伝子情報を含むデータを含むデータベースを用いた、目的遺伝子を細胞で発現させるための組換え遺伝子を作製するための機能性配列選択システムであって、前記目的遺伝子のコーディング領域の塩基配列もしくはアミノ酸配列またはそれらの一部を含む配列をクェリとして入力するための入力装置と、機能性配列を選択するための選択装置と、選択した機能性配列を出力するための出力装置と、を備え、前記選択装置は、前記データベースにおいて、前記目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、当該塩基配列または当該アミノ酸配列をコードする塩基配列と相同性を有する塩基配列を検索するための相同性配列検索手段と、検索して得られた塩基配列を第１の配列リストとして格納するための第１の配列リスト記憶手段と、第１の配列リストから、塩基配列全体がゲノム由来である塩基配列を除去して、塩基配列全体がゲノム由来である塩基配列以外の塩基配列を抽出し、
（１）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流または下流に塩基配列が存在する場合、上流または下流に存在する該塩基配列が機能性配列かどうか解析して１つまたはそれ以上の第１の機能性配列を抽出し、
（２）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流かつ下流に塩基配列が存在しない場合、前記データ内の前記遺伝子情報に、機能性配列を示す記載があるかどうかを検索して１つまたはそれ以上の第２の機能性配列を抽出するための機能性配列抽出手段と、
第１の機能性配列及び第２の機能性配列を含む第２の配列リストを格納するための第２の配列リスト記憶手段と、第２の配列リストに対し、スコアリングを行うためのスコアリング手段と、前記スコアリングの結果の上位の機能性配列を選択するための機能性配列選択手段と、選択された前記上位の機能性配列を格納するための機能性配列記憶手段と、を備
える機能性配列選択システムである。前記遺伝子情報に、ゲノム由来である塩基配列を示すキーワードがあるかどうか検索することによって、全体がゲノム由来である塩基配列以外の塩基配列が抽出されてもよい。前記機能性配列は、プロモーター配列、ターミネーター配列、及び終止コドンを含んでもよい。（１）の場合において、機能性配列ライブラリーまたは前記遺伝子情報を用いて、上流または下流に存在する該塩基配列が機能性配列かどうか解析されてもよい。第１の機能性配列及び第２の機能性配列に対し、出現頻度を基準として、または各機能性配列が由来する前記遺伝子情報に出現する特定のキーワードの頻度を基準としてスコアリングが行われてもよい。第１の機能性配列及び第２の機能性配列に対し、スコアリングを行う前に、第１の機能性配列及び第２の機能性配列から、前記目的遺伝子が天然に有する機能性配列が除外されてもよい。
【０００９】
本発明のさらなる実施態様は、遺伝子もしくはその一部の配列、遺伝子がコードするアミノ酸の配列もしくはその一部の配列、および／または前記遺伝子の遺伝子情報を含むデータを含むデータベースを用いた、目的遺伝子を細胞で発現させるための組換え遺伝子を作製するための機能性配列選択方法であって、前記データベースにおいて、前記目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、当該塩基配列または当該アミノ酸配列をコードする塩基配列と相同性を有する塩基配列を検索する工程と、検索して得られた塩基配列またはアミノ酸配列をコードする塩基配列のうち、全体がゲノム由来である塩基配列を除去して、全体がゲノム由来である塩基配列以外の塩基配列を抽出し、
（１）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流または下流に塩基配列が存在する場合、上流または下流に存在する該塩基配列が機能性配列かどうか解析して１つまたはそれ以上の第１の機能性配列を抽出し、
（２）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流かつ下流に塩基配列が存在しない場合、前記遺伝子情報の記載に、機能性配列を示す記載があるかどうかを検索して１つまたはそれ以上の第２の機能性配列を抽出する工程と、第１の機能性配列及び第２の機能性配列に対し、スコアリングを行う工程と、前記スコアリングの結果の上位の機能性配列を選択する工程と、を含む機能性配列選択方法である。
【００１０】
本発明のさらなる実施態様は、上記いずれかの機能性配列選択方法を前記機能性配列選択システムに実行させるためのプログラムである。また、本発明のさらなる実施態様は、このプログラムを格納したコンピュータ可読記憶媒体である。
【発明の効果】
【００１１】
本発明によれば、新規な機能性配列選択方法及び機能性配列選択システムを提供することができるようになった。
【図面の簡単な説明】
【００１２】
【図１】本発明の一実施形態にかかる機能性配列選択方法における各工程の概略図である。
【図２】本発明の一実施形態にかかる機能性配列選択方法において、第３工程から第５工程を詳細に示すフローチャートである。
【図３】本発明の一実施形態にかかる遺伝子配列設計システムの模式図である。
【図４】本発明の一実施形態における検索インタフェースの模式図である。
【図５】本発明の一実施形態における第１工程から第４−２工程までのフローチャーである。
【図６】本発明の一実施形態における第４−２工程から第５工程までのフローチャートである。
【図７】本発明の一実施形態における第４−２工程において、検索配列の上流及び下流に塩基配列が存在しないデータの相同性検索結果の一例である。
【図８Ａ】本発明の一実施形態において、図７に示す特許文献からの機能性配列領域の抽出方法を示す図である。
【図８Ｂ】本発明の一実施形態で用いることのできる機能性配列ライブラリーの一例である。
【図９】本発明の一実施形態における第５工程から第７工程までのフローチャートである。
【図１０】本発明の一実施形態における出力結果表示インタフェースを示す模式図である。
【図１１Ａ】本発明の一実施形態における出力結果インタフェースのプロモーター配列のポップアップ画面を示す。
【図１１Ｂ】本発明の一実施形態における出力結果インタフェースのプロモーター配列のポップアップ全件表示画面を示す。
【図１２Ａ】本発明の一実施形態における出力結果インタフェースのターミネーター配列のポップアップ画面を示す。
【図１２Ｂ】本発明の一実施形態における出力結果インタフェースのターミネーター配列のポップアップ全件表示画面を示す。
【図１３】本発明の一実施形態における出力結果インタフェースのプロモーター・ターミネーターペア配列のポップアップ画面を示す。
【図１４】本発明の一実施形態における出力結果インタフェースのクェリの遺伝子塩基配列のポップアップ画面を示す。
【図１５】本発明の一実施形態において、図１０で示す出力結果インタフェースのダウンロード表示結果を示す。
【発明を実施するための形態】
【００１３】
本発明の目的、特徴、利点、およびそのアイデアは、本明細書の記載により、当業者には明らかであり、本明細書の記載から、当業者であれば、容易に本発明を再現できる。以下に記載された発明の実施の形態および具体的な実施例などは、本発明の好ましい実施態様を示すものであり、例示または説明のために示されているのであって、本発明をそれらに限定するものではない。本明細書で開示されている本発明の意図並びに範囲内で、本明細書の記載に基づき、様々な改変並びに修飾ができることは、当業者にとって明らかである。
【００１４】
＝＝機能性配列選択方法＝＝
本発明の一実施形態にかかる機能性配列選択方法を用いる。すなわち、遺伝子もしくはその一部の配列、遺伝子がコードするアミノ酸の配列もしくはその一部の配列、および／または遺伝子の遺伝子情報を含むデータを含むデータベースを用いた、目的遺伝子を細胞で発現させるための組換え遺伝子を作製するための機能性配列選択方法であって、データベースにおいて、目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、当該塩基配列または当該アミノ酸配列をコードする塩基配列と相同性を有する塩基配列を検索する工程と、検索して得られた塩基配列またはアミノ酸配列をコードする塩基配列のうち、全体がゲノム由来である塩基配列を除去して、全体がゲノム由来である塩基配列以外の塩基配列を抽出し、
（１）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流または下流に塩基配列が存在する場合、上流または下流に存在する該塩基配列が機能性配列かどうか解析して１つまたはそれ以上の第１の機能性配列を抽出し、
（２）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流かつ下流に塩基配列が存在しない場合、前記遺伝子情報の記載に、
機能性配列を示す記載があるかどうかを検索して１つまたはそれ以上の第２の機能性配列を抽出する工程と、第１の機能性配列及び第２の機能性配列に対し、スコアリングを行う工程と、スコアリングの結果の上位の機能性配列を選択する工程と、を含む機能性配列選択方法である。
【００１５】
本発明の機能性配列選択方法は、後述する機能性配列選択システムにおいて実行される。すなわち、遺伝子もしくはその一部の配列、遺伝子がコードするアミノ酸の配列もしくはその一部の配列、および／または前記遺伝子の遺伝子情報を含むデータを含むデータベースを用いた、目的遺伝子を細胞で発現させるための組換え遺伝子を作製するための機能性配列選択方法であって、前記目的遺伝子のコーディング領域の塩基配列もしくはアミノ酸配列またはそれらの一部を含む配列をクェリとして入力するための入力装置と、機能性配列を選択するための選択装置と、選択した機能性配列を出力するための出力装置と、を備えた機能性配列選択システムにおいて、入力装置がクェリを入力し、選択装置において、相同性配列検索手段が、データベースにおいて、目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、当該塩基配列または当該アミノ酸配列をコードする塩基配列と相同性を有する塩基配列を検索し、第１の配列リスト記憶手段が、検索して得られた塩基配列を第１の配列リストとして記憶し、機能性配列抽出手段が、第１の配列リストから、塩基配列全体がゲノム由来である塩基配列を除去して、塩基配列全体がゲノム由来である塩基配列以外の塩基配列を抽出し、
（１）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流または下流に塩基配列が存在する場合、上流または下流に存在する該塩基配列が機能性配列かどうか解析して１つまたはそれ以上の第１の機能性配列を抽出し、
（２）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の上流かつ下流に塩基配列が存在しない場合、前記データ内の前記遺伝子情報に、機能性配列を示す記載があるかどうかを検索して１つまたはそれ以上の第２の機能性配列を抽出し、
第２の配列リスト記憶手段が、第１の機能性配列及び第２の機能性配列を含む第２の配列リストを記憶し、スコアリング手段が、第２の配列リストに対し、スコアリングを行い、機能性配列選択手段が、スコアリングの結果の上位の機能性配列を選択し、機能性配列記憶手段が、選択された上位の機能性配列を記憶し、出力装置が、選択された前記上位の機能性配列を出力する、機能性配列選択方法である。
【００１６】
以下、本方法について、詳細に説明する。
【００１７】
本実施形態にかかる方法は図1に示す第1〜7で構成される工程からなる。以下、各工程
について詳細に説明する。
【００１８】
（第１工程）
第１工程であるクェリ入力工程１は、細胞に導入したい目的遺伝子のコーディング領域の塩基配列もしくはアミノ酸配列またはそれらの一部を含む配列をクェリとして入力し、同方法を開始する工程である。すなわち、クェリは、以下の（１）もしくは（２）からなる配列、または（１）もしくは（２）を含む配列である。
（１）目的遺伝子のコーディング領域の塩基配列またはその一部
（２）目的遺伝子がコードするアミノ酸配列またはその一部
配列が登録されている名前を用いる場合は、例えば遺伝子名から遺伝子配列を出力する前処理を行なった後に、その配列をクェリとして入力すればよい。クェリとして入力する配列数は複数でもよいが、その場合は複数の配列をそれぞれ識別した状態で入力することが好ましい。
【００１９】
また、クェリに付随して、図1に示す第２工程（相同性検索工程２）、第４工程（検索
結果解析工程４）、第６工程（スコアリング実行工程６）、第７工程（機能性配列領域出力工程７）における各種条件を入力してもよい。
【００２０】
（第２工程）
第２工程である相同性検索工程２は、第１工程でクェリとして入力した目的遺伝子のコーディング領域の塩基配列もしくはアミノ酸配列またはそれらの一部を含む配列を検索配列として、当該検索配列と相同性の高い塩基配列またはアミノ酸配列を検索する工程である。このような遺伝子の塩基配列またはアミノ酸配列の相同性検索は、遺伝子もしくはその一部の配列、遺伝子がコードするアミノ酸の配列もしくはその一部の配列、および／または遺伝子情報を含む生物情報データベースに対して実行する。このようなデータベースとして、研究者により頻繁に用いられるものとしては、NCBIのEntrez （http://www．ncbi．nlm．nih．gov/entrez/）や日本の国立遺伝学研究所（DDBJ:https://www.ddbj.nig.ac.jp/index.html）、ヨーロッパ分子生物学研究所（EMBL: European Molecular Biology Laboratory）の中の EBI（European Bioinformatics Institute）、NCBIが提供するLocusLink （http：//www．Ncbi．nlm．nih．gov/LocusLink/）、タンパク質情報中心のSWISS-PROT （http：//kr．Expasy．org/sprot/） 等が良く知られている。どのデータベースを
用いるかは、あらかじめ決められていても良いし、第１工程のクェリに付随して入力してもよい。
【００２１】
データベースを用いて相同性検索を行なう方法としては、例えば塩基配列あるいはアミノ酸配列の相同性を調べるツールであるBLAST（Basic Local Alignment Search Tool）やSSEARCH（http://www.biology.wustl.edu/gcg/ssearch.html）などが知られており、同ツールを用いてクェリと相同性の高い遺伝子配列情報を抽出してもよい。このとき、相同性判定基準を、第１工程のクェリに付随して入力してもよい。相同性判定基準に関しては、例えばBLASTなどではe-value値を設定することができる。
【００２２】
ここで、遺伝子情報とは、遺伝子の特徴の記載であって、例えば、データベースに登録されているデータ名やデータＩＤ、データが記載されている文献情報（例えば、文献発行年月日、文献要旨、著者、文献内に含まれる塩基配列やアミノ酸配列とそれらの由来、などの一部または全部）、などを含む。
【００２３】
（第３工程）
第３工程である第一配列リスト記憶工程３は、第２工程で検索した相同性の高い配列情報からリストを作成し、記憶する工程である。ここでリスト化する情報としては、検索より得られた高い相同性を持つ塩基配列もしくはアミノ酸配列またはアミノ酸配列をコードする塩基配列だけでなく、塩基配列上やアミノ酸配列上で検索配列と相同性の高い領域の位置情報の他、データベースに登録されているデータ名やデータＩＤ、データが記載されている文献情報などが挙げられる。検索配列と相同性の高い領域の位置情報とは、例えば、登録されている塩基配列の５’末端からの塩基数などが挙げられるが、その限りではない。
【００２４】
（第４工程）
第４工程である検索結果解析工程４は、第３工程でリスト化された塩基配列及び情報を、天然配列とその他の配列へ分類し（第４-１工程）、機能性配列領域を抽出・解析した
後に（第４-２工程）、天然配列と人工配列に再分類する（第４−３工程）工程である。
本工程のフローチャートを図２に示し、以下説明する。
【００２５】
第４-１工程：天然配列と人工配列の分類
本工程では、第３工程にて得られた第一配列リスト１０のデータを、天然配列分類工程
１１により、天然配列と人工配列の候補であるその他の配列に分類する。ここで言う天然配列とは、全体がゲノム由来である、すなわち生物のゲノム上に存在する配列のみが含まれる塩基配列を意味し、人工配列とは、全体がゲノム由来である配列以外の配列、すなわち天然には存在しない配列が含まれる塩基配列のことを示す。このような分類には、天然配列の場合、例えば、相同性検索を行なう際にヒットしたデータに含まれる遺伝子情報に、天然に存在する配列であることを示すキーワードを含むものを使用することができる。このキーワードとして、例えば「complete genome」や「genome」または「chromosome」
などが該当するが、その限りではない。また、その他の配列は上記キーワードを含まないデータが挙げられる。そのため、分類としては、天然配列以外のデータをその他の配列として分けることが可能である。また、人工配列は、天然配列以外の配列としたり、データ情報に人工的に設計した配列であることを示すキーワードを含むものとしたりして、分類することができる。このキーワードとして、例えば「synthesis」や「mutant」やその類
義語などが上げられるが、その限りではない。第一配列リスト１０から天然配列分類工程１１により天然配列とその他の配列として分類した後に、第４−２工程で機能性配列領域を抽出し、解析することが望ましいが、工程の順序はこの限りではない。
【００２６】
第４−２工程：機能性配列領域の抽出・解析
機能性配列領域の抽出・解析は、第４−１工程で分類されたその他の配列リスト１３において、遺伝子が生物内で機能するために必要な機能性配列を含む領域を抽出し、機能性配列を解析する工程である。このような機能性配列としては、プロモーター配列、ターミネーター配列、終止コドンなどが挙げられるが、それらに限定されない。ここで、プロモーター配列とは、存在する転写開始点の上流、たとえば転写開始点から上流約３００ｂｐまで、１００ｂｐまで、あるいは６０ｂｐまでに存在し、転写の開始に必要な塩基配列のことであって、RNAポリメーラーゼ等の基本転写因子の結合に寄与する領域の配列を意味
する。
【００２７】
具体的な機能性配列領域の解析方法は、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同領域の下流に塩基配列が存在するかどうかを確認するが（機能性配列領域抽出工程１４）、その結果によって解析方法が異なる。
【００２８】
検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同性領域の上流または下流に塩基配列が存在する場合は、図２の機能性配列領域抽出・解析工程１５に示すように、コーディング配列の端から一定塩基長を抽出した後に、同抽出領域内に機能性配列が存在するかを解析する。ここで、上流及び下流とは、それぞれクェリの配列の相同性領域と隣接した、または離れた上流側及び下流側を意味し、上流は、クェリの配列の相同性領域より上流でもよく、開始コドンの上流でもよく、転写開始点から上流でもよい。また、下流は、クェリの配列との相同性領域より下流でもよく、終始コドンを含んで、そこから下流でもよく、終始コドンの下流でもよい。また一定塩基長とは３塩基以上１００００塩基以下であればよく、好ましくは１０塩基以上５００塩基以下であり、より好ましくは、２０塩基以上２００塩基以下である。この抽出塩基長は、上流と下流で、独立して選択することができる。このように抽出した配列において、通常上流配列にプロモーター配列が、下流配列にターミネーター配列及び／又は終止コドンが含まれている。これらの配列を解析する方法としては、プロモーター配列やターミネーター配列の機能性配列データベース、コマーシャルに用いられている配列を検索する方法（図２において機能性配列領域抽出・解析工程１５中ＤＢと記載）、機能性配列推測プログラムを用いてプロモーター、ターミネーター配列を推測する方法（図２において機能性配列領域抽出・解析工程１５中推定ツールと記載）、抽出配列中に他の遺伝子コード配列が存在するかを確認し、該当領域以外の配列を機能性配列と設定する方法（機能性配列領域抽出・解析工程１５中その他と記載）がある。ＤＢにより同定する方法において用いるデータベースとしては、例えばRegulon DBやiGEMなどが知られている他、コマーシャルに用いられている配列としては
AddgeneやGene Scan、TaKaRaバイオ、Thermo Scientific、Invitrogen等で取り扱ってい
る製品の配列情報を確認することができる。この場合は、各データベースやコマーシャル情報のライブラリーである機能性配列ライブラリーを作製し、一致する、または類似する配列を同定することができる。また、推測ツールにより同定する方法としては、機械学習手法を用いた推測方法等を用いることができる。このような方法としては、CNNpromoter(http://www.softberry.com/berry.phtml？topic=index&group=programs&subgroup=deeplearn)や、Findtermtopic=index&group=programs&subgroup=deeplearn)や、Findterm（http://www.softberry.com/berry.phtml?topic=findterm&group=programs&subgroup=gfindb）などが知られており、同技術、または同技術を修飾したプログラムにて推測することが可能となる。さらにその他の方法として、抽出配列中に他の遺伝子配列が含まれているかを調べる方法としては、抽出配列をクェリとしてBLASTやSSEARCHにて検索する方法がある。同検索により遺伝子配列をコードしている領域を同定し、同領域以外の配列を機能性配列領域として抽出することができる。また、遺伝子情報に記載があれば、その情報を用いてもよく、文献から、テキストマイニング技術を用いて機能性配列を抽出してもよい。テキストマイニング技術については、後述する。
【００２９】
検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同性領域の上流および下流の配列が存在しない場合、図２の機能性配列領域解析工程１６に示すように、文献からテキストマイニング技術を用いて機能性配列を抽出する方法を用いる。例えば、特許文献内に相同性の高い遺伝子配列が登録されていた場合は、遺伝子情報には、該当する特許文献情報と、同文献内で登録されている遺伝子配列のＩＤ情報が記載されている。そこで、特許本文中の同配列ＩＤが記載されている文章、及び段落内に、関連するキーワードが存在するかを確認し、同キーワードを抽出することでプロモーターまたはターミネーター配列情報を得ることができる。このようなキーワードとしては、「promoter」、「terminator」などがあり、同単語と同一文章中の塩基配列を抽出することができる。また、上記以外にも、プロモーター名、ターミネーター名、ベクター名やプラスミド名が記載されている場合もある。このような場合は、機能性配列領域抽出・解析工程１５のＤＢでの検索と同様に、データベースやコマーシャルの製品情報を基に、上記の名前と塩基配列情報が格納された機能性配列ライブラリーを作製することで、配列情報を抽出することが可能となる。  
本工程において、機能性配列として終止コドンも検索することが好ましい。終止コドンは、機能性配列領域抽出・解析工程１５及び機能性配列領域解析工程１６のいずれにおいても、抽出した下流の塩基配列の５’末端の３塩基において、終止コドンとして登録されているＤＮＡ配列（ＴＡＡ、ＴＡＧ、ＴＧＡ）を検索し、これらのいずれかが存在する場合に第６工程にて出力することができる。この際、もし上記３種類の配列が存在しない場合は、終止コドンを出力として提示しないか、もしくはＴＡＡを付与して提示してもよい。
【００３０】
機能性配列のうち、開始コドンは概ねクェリとして入力した塩基配列の５’末端に３塩基存在している。そのため、クェリの５’末端の３塩基が後述する開始コドンから始まっているかを確認し、開始コドンがあれば、それを選択する。開始コドンをコードする塩基配列としては、真核生物では原則的にＡＴＧであるが、原核生物ではＡＴＧ以外にＧＴＧ、ＴＴＧが知られている。
【００３１】
第４−３工程：天然配列と人工配列の再分類
検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同性領域の上流または下流の配列が存在する場合と上流および下流の配列が存在しない場合とで、それぞれ抽出された機能性配列を合わせて、データベースを作成する。
【００３２】
このデータベースに対し、オプションとして、第４−２工程で抽出・解析された機能性
配列領域及びそれらの情報は、その他の配列に分類された機能性配列領域が天然配列に分類されたものと同じ配列を持つ場合は天然配列として天然配列リスト１８へ再分類し（再分類工程１７）、該当しないその他の配列を人工配列リスト１９へ分類してもよいが、再分類することなく、第二配列リスト２０としてもよい。
【００３３】
（第５工程）
第５工程である第二配列リスト記憶工程５は、図２に示すように、第４工程で得られた人工配列リスト１９に分類された機能性配列領域情報から第二配列リスト２０を作成し、記憶する工程である。このとき、配列に付随して、第一配列リストから紐付く各データの情報や、機能性配列の抽出及び解析方法をリストに入れてもよい。特に機能性配列解析において、ＤＢによる同定を行なった際には、由来となったDB名や配列名も情報としてリスト化されることが好ましい。
【００３４】
（第６工程）
第６工程であるスコアリング実行工程６は、第５工程で得られる第二配列リストにおいて、所望の観点よりスコアリングを実行する工程である。所望の観点とは、例えば頻度が挙げられる。この場合は、第二配列リスト２０において登録されている数をスコアとして用い、最も数が多いものをスコア上位としてプロモーター配列、ターミネーター配列、または配列ペアを提示する。別の所望の観点としては、機能性配列へ求める機能が挙げられる。例えば、ある物質合成を生物へ付与するための遺伝子に必要な機能性配列情報を調べる場合、１、遺伝子を導入する生物種、２、目的とする物質合成の実績の有無、３、導入遺伝子の生物内の量（発現量）、などの機能が関連する。このような場合は、１、２、３、に関連するキーワードを第１工程においてクェリに付随させ、第４工程でこれらのキーワードが該当情報中にどの程度存在するかをカウントし、得られたカウント数をスコアとして上位から提示してもよい。ここでクェリとするキーワードとしては、上記１、２、３に該当するもの以外でもよい。またキーワードを検索する該当情報としては、例えばデータ名や文献などが挙げられるが、それらに限定されない。
【００３５】
（第７工程）
第７工程である機能性配列領域出力工程７は、第６工程で実施するスコアリングを基に、上位となった機能性配列領域を選択し、出力する工程を示す。機能性配列領域は、リストとして出力されてもよい。ここで選択される塩基配列は、プロモーター配列または領域、クェリとした遺伝子配列（開始コドンを含む）、終止コドン、ターミネーター配列または領域から選択される１つ以上の機能性配列が選択されることが望ましい。ここで言うプロモーター配列領域は、プロモーター配列と開始コドンの間、もしくはプロモーター配列より上流に、更なる配列が付与されている領域を示す。更なる配列があらかじめ決められている場合は、出力の際にプロモーター配列の前後に付与する設定とすることもできる。このような条件は、第１工程のクェリに付随して入力しておくことができる。例えば、真核生物では、プロモーター配列と開始コドンの間にKozak配列やTATAbox配列を付与することができ、原核生物にもコンセンサス配列が存在することが知られており、生物種に応じて適宜必要配列を付与してもよい。ターミネーター配列領域も、プロモーター配列領域と同様に、ターミネーター配列の前後に更なる配列が付与されてもよく、更なる配列があらかじめ決められている場合は、出力の際にターミネーター配列の前後に付与してもよい。
【００３６】
＝＝発現ベクター設計方法＝＝
機能性配列選択方法で選択された機能性配列を、目的遺伝子の前後に結合し、発現ベクターのエンハンサーの下流の適正の位置に挿入する。それによって、目的遺伝子の発現レベルの高い発現ベクターを作製することができる。
【００３７】
このように、発現ベクターを作製するにあたって、本開示の機能性配列選択方法によっ
て、データベースから配列情報抽出を行なうことで、自動的に、そして遺伝子配列や生物種を限定することなく、網羅的に機能性配列領域情報を抽出することができる。さらに、遺伝子導入目的が定まっている場合は、目的に適した機能性配列情報をあらかじめクェリに付随して入力してスコアリングすることで、目的に最も適した機能性配列を出力として得ることができ、至適機能性配列情報を短時間で自動的に得ることが可能になる。得られた機能性配列を用いて発現ベクターを作製することにより、遺伝子導入に最適の発現ベクターを構築することができる。
【００３８】
＝＝機能性配列選択システム＝＝
本実施形態の機能性配列選択システムは、機能性配列選択方法の第１から７の工程を行うシステムである。目的遺伝子のコーディング領域の塩基配列もしくはアミノ酸配列またはそれらの一部を含む配列をクェリとして入力するための入力装置と、機能性配列を選択するための選択装置と、選択した機能性配列を出力するための出力装置と、を備え、選択装置は、データベースにおいて、目的遺伝子のコーディング領域の塩基配列、アミノ酸配列をコードする塩基配列、またはアミノ酸配列を用いて、塩基配列またはアミノ酸配列をコードする塩基配列と相同性を有する塩基配列を検索するための相同性配列検索手段と、検索して得られた塩基配列を第１の配列リストとして格納するための第１の配列リスト記憶手段と、第１の配列リストから、塩基配列全体がゲノム由来である塩基配列を除去して、塩基配列全体がゲノム由来である塩基配列以外の塩基配列を抽出し、
（１）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同性領域の上流または下流に塩基配列が存在する場合、上流または下流に存在する該塩基配列が機能性配列かどうか解析して１つまたはそれ以上の第１の機能性配列を抽出し、
（２）抽出した塩基配列に、検索に用いた塩基配列またはアミノ酸配列をコードする塩基配列の相同性領域の上流かつ下流に塩基配列が存在しない場合、データ内の前記遺伝子情報に、機能性配列を示す記載があるかどうかを検索して１つまたはそれ以上の第２の機能性配列を抽出するための機能性配列抽出手段と、
第１の機能性配列及び第２の機能性配列を含む第２の配列リストを格納するための第２の配列リスト記憶手段と、第２の配列リストに対し、スコアリングを行うためのスコアリング手段と、スコアリングの結果の上位の機能性配列を選択するための機能性配列選択手段と、選択された前記上位の機能性配列を格納するための機能性配列記憶手段と、を備える。以下、具体的なシステムについて、詳細に述べる。
【００３９】
図３に一実施形態の遺伝子設計システム３０のブロック構成図を示す。本実施形態の機能性配列選択システムは、入力装置としてユーザーがクェリを入力する入出力端末３１と、本システムへアクセスするために利用するインターネット３２と、出力装置としてクェリ及び出力情報を提示するインタフェースとなる情報制御部３３と、クェリをシステムへ導入するクェリ導入部３４と、相同性配列検索手段として生物情報データベース３５を用いてクェリの相同性検索を行なうデータベース検索部３６と、第１の配列リスト記憶手段としてクェリの相同性解析結果を格納する検索結果保持部３７と、機能性配列抽出手段として、解析に用いる機能性配列ライブラリー３８及び検索結果解析部３９と、第２の配列リスト記憶手段として同解析結果を保持する解析結果保持部４０（第ニ配列リスト）と、スコアリング手段として解析結果保持部４０のデータを基にスコアリングを行うスコア解析部４１と、機能性配列選択手段及び機能性配列記憶手段としてスコア解析結果を基に出力情報を抽出する出力制御部４２からなる。出力制御部４２により、情報制御部３３を介して解析結果が表示される。
【００４０】
なお、データベース検索部３６、検索結果保持部３７、検索結果解析部３９、解析結果保持部４０、スコア解析部４１において、順に図1中の第２工程（相同性解析工程２）、
第３工程（第一配列リスト作製工程３）、第４工程（検索結果解析工程４）、第５工程（
第二配列リスト作製工程５）、第６工程（スコアリング実行工程６）が実行され、詳細は上述した各工程の説明に記載されている。機能性配列ライブラリーは、プロモーター名、ターミネーター名、ベクター名やプラスミド名とそれらの配列がリスト化されているライブラリーを意味し、詳細は第４−２工程の説明に記載されている。
【００４１】
以下、本システムを用いた機能性配列選択方法を概略する。まず、本実施形態の機能性配列選択システムにおいて、ユーザーが入出力端末３１より入力した遺伝子配列情報と検索、解析、スコアリング、出力条件をクェリ導入部３４によって受け付ける。クェリ導入部３４はデータベース検索部３６、検索結果解析部３９、スコア解析部４１、出力制御部４２へクェリ及び付随情報を導入する。クェリ導入部３４を介して、遺伝子配列情報及び検索条件をキーとして生物情報データベース３５に対してデータベース検索部３６にて相同性検索を実行し、その結果を第一配列リストとして検索結果保持部３７へ格納する。この時、検索条件として、生物情報データベース３５の種類や相同性判定基準を入力することができる。次に、クェリ導入部３４を介して導入された解析条件をキーとして、検索結果保持部３７のデータを天然配列及び人工配列に分類し、機能性配列領域を抽出・解析し、必要があれば天然配列の再分類を行なうことを検索結果解析部３９にて実行し、その結果を第二配列リストとして解析結果保持部４０に格納する。この時、解析条件としては、機能性配列領域の抽出塩基長を入力することができる。さらに、機能性配列の解析において参考とするDBは機能性配列ライブラリー３８の情報を参考に行なうこともできる。さらに、クェリ導入部３４を介して導入されたスコア条件をキーとして、解析結果保持部４０のデータをスコア解析部４１にてスコアリングを実行する。この時、スコア条件としては、頻度、または所望の機能に関するキーワードを入力することができる。ここで得られたスコア上位の機能性配列または領域及びクェリの遺伝子配列を、出力制御部４２から情報制御部３３を介して表示することができる。ここで、出力条件として、機能性配列領域に付与したい配列がある場合は、クェリとして設定することができ、その場合は入出力端末３１へ当該条件を入力することで設定することができる。
【００４２】
＝＝具体的な機能性配列選択方法＝＝
［１］本実施形態では、図１に記載する遺伝子設計方法を基本構成とした図３に記載の遺伝子設計システムを用いて、第１工程から第５工程である第二配列リスト作製工程までを実施した一例を示す。
【００４３】
以下、図４〜６を参照しながら詳細に説明する。図４に検索インタフェースと、図５、６に第１工程から第４−１及び４−２工程の機能性配列領域抽出までと、第４−２工程の機能性配列領域解析から第５工程までのフローチャートを示す。
【００４４】
まず、検索インタフェース５０へ入力を行なうことによりクェリ入力工程６０を実行する。本実施例のクェリである遺伝子のアミノ酸配列を、遺伝子配列入力部５１へ直接入力するか、または５２へテキストファイルとして入力できる。また、相同性検索工程のためにBLAST検索６１を実行するが、クェリの付随情報の検索条件として、本工程で用いる生
物情報データベース６２の選択および相同性判定基準をインタフェースに入力する。本実施例では、BLAST検索６１としてtBLASTnを用い、生物情報データベース６２にはNCBIに登録されているNucleotide collectionと、Patent sequenceのデータベースの選択をデータベース選択５３にて実行し、且つ相同性判定基準５４もe-value値を入力する。さらに、
第４−２工程で行なう機能性配列領域の塩基配列長も同画面の塩基長選択５５にて選択することができる。これらの条件を入力・選択した後、インタフェースの検索開始ボタン５６をクリックすることで、相同性検索が開始され、以降の工程が自動で実施されるようになる。
【００４５】
BLAST検索６１を実施し、検索結果を第一配列リスト６３へ格納する。格納されたデー
タ情報のうち、データ名に「complete genome」、または「chromosome」が含まれている
かをキーワード検索６４にて実行し、含まれる場合は天然配列リスト６５へ、含まれない場合はその他の配列リスト６６へデータを格納することで、第４−１工程を実施する。
【００４６】
次に、抽出工程６７において、その他の配列リスト６６を用いて相同性解析を行い、及びリストの配列と検索配列のアライメントを行い、クェリのアミノ酸配列がコードする塩基配列の相同性領域の上流、下流の配列を抽出する。このとき、フレーム確認工程６９として、抽出した配列情報のうち、フレームがクェリのアミノ酸配列をコードする側と同じ（プラス）か、反対（マイナス）かを生物情報データベース６８より確認する。そして、第一配列リストを参照して、フレーム情報を第一配列リストに格納してもよい。フレームがマイナスであった場合、クェリのアミノ酸配列をコードする塩基配列と逆の塩基配列が表示され、且つ上流と下流の配列が逆の位置となるため、抽出した配列において、相補配列対応工程７０を実施する。フレームの統一を行なった後、下流の抽出配列の5’末端に
終止コドンが存在しているかどうかを確認し（７１）、これらの変更及び追加情報をその他の配列リストに追加することで（７２）、第４−２工程の機能性配列領域の抽出を実行する。
【００４７】
続いて、その他の配列リスト８０へ機能性配列領域が抽出・格納された後に、それぞれのデータにおいて上流と下流の抽出配列があるか（８１）を実行する。以下、抽出配列がある場合と無い場合の工程に関して記載する。
【００４８】
抽出配列が存在する場合は、抽出配列をデータベース等に登録されている情報を機能性配列ライブラリー８２に格納し、これらのデータを基に検索する（８３）。本実施例では、機能性配列ライブラリー８２へ、RegulonDB、iGEM、Addgeneに登録されているプロモーター、ターミネーター、ベクター、プラスミド名及び配列を格納するが、本ライブラリーには追加でデータを入力でき、たとえば、コマーシャルにて使用されるSnap Gene, Invitrogen、Takaraバイオ等で取り扱っている上記情報を格納できる。DB検索８３の検索結果
チェック工程８４においてヒットした配列は、その他の配列リストへ変更データを追加する（８５）が、ヒットしなかった配列は、次に推定ツール検索８６により機能性配列の同定を行なう。本ツールにはCNNpromoterやFindtermを用いる。推定ツール検索８６の検索
結果チェック工程８７においてヒットしなかった配列は、次にその他の検索８８により機能性配列の同定を行なう。その他の検索にはBLASTを用いて抽出配列中に遺伝子配列とし
てコードされている領域があるかを確認する。その他の検索８８の検索結果チェック工程８９において、遺伝子配列領域と相同性を持つ領域がある場合は、遺伝子配列領域以外の配列領域を機能性配列領域としてその他の配列リストへ変更データを追加するが、存在しなかった配列は、その旨の変更データをその他の配列リストへ追加する。
【００４９】
抽出配列が存在しない場合は、第一配列リストに登録されている遺伝子情報を基に、機械学習技術により機能性配列領域情報を抽出する（９０）。具体的には、データが記載されている文献情報を用いて、特許の場合は配列ＩＤを、論文の場合はデータ名を検索し、同情報が記載されている文章またはセンテンスより、プロモーター、ターミネーター、ベクター、プラスミド名の記述を抽出する。その一例を図７に示す。
【００５０】
図７では、大腸菌のＭＥＰ代謝経路に関与する酵素であるispE遺伝子のアミノ酸配列をクェリとして本ツールを用い特許文献より配列情報を抽出する。まず、図５の相同性検索６１により、たとえば図７に示す特許文献情報１００が得られる。ここでは、相同性を持つデータ名１０１や相同性配列１０２、データベースsequence ID１０３、フレーム情報
１０４等が得られ、これらの情報は第一配列リストへ格納される。ここで得られたデータ名１０１やデータベース上のsequence ID １０３より、特許文献情報と、その特許内にて登録されているSequence ID情報が得られる。図７に示すデータの場合は、データベース
名１０１から、特許文献情報US6234514と文献中のSequence ID(SEQ ID)No:5にクェリの配列が登録されていることが示されている。ここでの検索式としては、SEQ ID No:5として
いる。これらの情報を基に、実際の特許文献中にSEQ ID No:5が記載されている文章を検
索する。図８Ａに検索イメージを示すが、特許本文１０９中、SEQ ID NO:5が記載されて
いる文章の前後を抽出し、抽出領域１１０中に、機能性配列ライブラリーに登録されている文言の有無を検索する。図８Ａの場合、抽出領域１１０中に含まれるpBAD TOPO TAという単語が機能性ライブラリーにベクター名として登録されているため、この情報よりプロモーター名とターミネーター名の抽出結果１１２を得ることができ、特許文献から機能性配列領域を抽出することができる。このような機能性配列ライブラリー１１３の例を図８Ｂに示す。機能性配列ライブラリー１１３はベクターまたはプラスミド名と同領域に含まれるプロモーター名及び配列、ターミネーター名及び配列とこれらの配列の天然、人工配列の分類、且つ出典が記載されている。このようなライブラリーは既存のプロモーター、ターミネーター、ベクター、プラスミド配列データベースや、コマーシャルで使用されているこれらの製品情報を記載する媒体より抽出して作製する。このように得られた機能性配列領域情報を、抽出配列が存在する場合と同様に、その他の配列リストへ変更データとして追加する（８５）。本工程により第４−２工程の機能性配列領域の解析を実行する。
【００５１】
次に、ゲノムデータベースなどを用いて、その他の配列リストの中に天然配列リストと同じ機能性配列及び領域が存在するかを検索結果チェック工程９１により確認する。同じ配列が存在する場合は、その他の配列リストから除去し（９２）、同じ配列が存在しない場合は、人工配列リストとして分類する（９３）。得られた人工配列リストを第二配列リスト９４として格納する。このように、導入したい遺伝子のアミノ酸配列をクェリとして入力することにより、遺伝子導入の際に使用される機能性配列領域を第二配列リストとして得ることができる。
【００５２】
［２］本実施形態では、図１に記載する遺伝子設計方法を基本構成とした図３に記載の遺伝子設計システムを用いて、第５工程から第７工程の機能性配列領域出力工程までを実施する例を示す。
【００５３】
図９に第５工程から第６工程のスコアリング工程と第７工程の機能性配列領域出力工程までのフローチャートと、図１０に出力結果表示インタフェースを示す。スコアリング条件は第１工程のクェリ入力の際に、頻度で行なうか、所望の条件におけるキーワードにおいて行なうかを選択する。この選択は図４の検索インタフェースのスコアリング条件選択５７にて行うことができる。
【００５４】
図９のスコアリング条件選択工程１２１において、頻度でのスコアリングを選択した場合、第二配列リスト１２０の情報より登録数をスコアリング対象として頻度別スコアリングを実施し、結果を頻度別スコア結果リストへ格納する（１２２）。また、頻度でのスコアリングを選択しなかった場合、所望のスコアリングキーワードの存在率をスコアとする。本スコアのキーワードは図４のスコアリング条件選択５７で入力できる。ここで入力されたキーワードは、キーワードスコアリング実行工程１２３において、遺伝子情報に記載の文献またはデータ名に、同キーワードが存在するかを検索しその数をスコアとしてカウントする。同スコアに基づいてキーワードスコア結果リストを作製する（１２４）。いずれのスコアリング方法で得られた結果でも、スコア上位を出力対象とし、機能性配列、つまりプロモーター、ターミネーター、終止コドン、開始コドン等を含む、クェリで入力した遺伝子配列の塩基配列を機能性配列領域出力工程１２５時に出力データとして選択し、その結果を図１０に示すインタフェースにより表示する（１２６）。図１０は検索インタフェースの下部にスコアリング上位の出力結果１３０が表示される例を示す。それぞれプロモーター配列、遺伝子の塩基配列、ターミネーター配列、プロモーター・ターミネーターペア配列が表示されるが、それぞれの表示領域のポップアップボタン１３１、１３２、
１３３、１３４をクリックすると、図１１、１２、１３、１４に示すようにそれぞれのポップアップ画面が表示される。
【００５５】
まずプロモーター配列において、ポップアップボタン１３１をクリックすると、図１１Aに示すように、各条件にて同定された際に得られたプロモーター名１４０と、プロモー
ター領域配列１４１と、天然配列（図中Ｎと表示）か人工配列（図中Ａと表示）か、両方に存在するか（図中ＮＡと表示）を示す分類欄１４２と、配列同定方法及び同定に用いたデータベース情報１４３と、頻度によるスコア１４４からなる。配列同定方法としては機能性配列ライブラリーを参照するＤＢ法と推測ツールを用いる方法、その他の方法がある。ＤＢ法で同定した場合は機能性配列ライブラリーに登録したデータベース名が、推定ツール及びその他の方法で推定された領域に関しては、推定ツールまたはその他（図中ハイフンで記載）と表示される。そして、人工配列のうち最もスコアが高いものがスコアリングの上位として表示されている。図１１のチェック箇所１４５により選択し、セレクトボタン１４６をクリックすることで、図１０に表示される結果の配列を変更することも可能である。さらに全件表示ボタン１４７をクリックすることにより、図１１Ｂに示すように、検索により得られた全結果を表示する（全結果表示１４８）ことができる。この操作は、図１２に示すターミネーター配列及び図１３に示すプロモータ・ターミネータ配列ペアのポップアップ画面においても可能となる。さらに、図１４には、本解析のクェリとして入力した遺伝子の塩基配列１７０のポップアップ画面を示すが、５’末端に開始コドン１７１、３’末端に終止コドン１７２が付与されている。クェリカバー率１７３も表示されており、クェリカバー率１００％となる塩基配列が図１０及び本ポップアップ画面に表示される。このように、各配列において、スコア上位もしくは任意で選択した配列を図１０中に表示し、ダウンロードボタン１３５をクリックすることにより、図１５の出力配列１８０に示すように、プロモーター配列１８１、開始コドン１８２、終止コドン１８３を含む遺伝子配列、ターミネーター配列１８４を一括で表示することができる。そして、細胞に導入したい遺伝子のアミノ酸配列の機能性配列領域が格納された第二配列リストより、機能性配列であるプロモーター、ターミネーター配列を同定または推定し、所望の条件でスコアリングを行なうことにより、自動で至適機能性配列情報を得ることが可能となる。以上のように、本ツールを用いることで遺伝子配列設計工程の短時間化を可能とすることができる。
＝＝プログラム及びコンピュータ可読記憶媒体＝＝
本発明の一実施態様は、上記の機能性配列選択システムに、上記の機能性配列選択方法を行わせるためのプログラムである。また、このプログラムをコンピュータによって読み取り可能に格納する記録媒体も、本発明の実施形態の一つである。これらによって、上記の機能性配列選択方法が汎用できるようになる。
【実施例】
【００５６】
以下、本発明を実施例に基づいてさらに具体的に説明するが、本発明はこれに限定されるものではない。当業者は、本発明の意義を逸脱することなく様々な態様に本発明を変更することができ、これらの変更も本発明の範囲に含まれる。
【００５７】
本実施例では、クェリを特定し、図６の工程を行ったところ、機能性配列が得られたことを示す。
【００５８】
［実施例１］
以下の条件で図６のデータベース検索（８３）を行ったところ、プロモーターとして、RegulonDB: ispFp2、ターミネーターとして、（iGEM:BBa_B0010)が得られた。
【００５９】
クェリの一つである遺伝子のアミノ酸配列：非メバロン酸経路関連遺伝子ispF
生物情報データベース：nucleotide collection , patent sequenceを選択
上流・下流の各抽出配列長：1000bp
［実施例２］
以下の条件で図６のデータベース検索（８３）を行ったところ、機能性配列が得られなかったので、blast search（８８）を行ったところ、遺伝子コード領域として「hypothetical protein AC239_22870」が、上流抽出塩基（1000nt)の5’側から196〜804ntに存在することが検出され、残りの196ntにlac promoter配列が含まれていることが検出された。
【００６０】
クェリの一つである遺伝子のアミノ酸配列：非メバロン酸経路関連遺伝子ispG
生物情報データベース６２：nucleotide collection , patent sequenceを選択
上流・下流の各抽出配列長：1000bp
【符号の説明】
【００６１】
１、６０…クェリ入力工程
２…相同性検索工程
３…第一配列リスト作製工程
４…検索結果解析工程
５…第二配列リスト作製工程
６…スコアリング実行工程
７、１２５…機能性配列領域出力工程
１０、６３…第一配列リスト
１１…天然配列分類工程
１２、１８…天然配列リスト
１３、６６…その他の配列リスト
１４…機能性配列領域抽出工程
１５…機能性配列領域抽出・解析工程
１６…機能性配列領域解析工程
１７…再分類工程
１９、６５…人工配列リスト
２０、９４、１２０…第二配列リスト
３０…遺伝子設計システム
３１…入出力端末
３２…インターネット
３３…情報制御部
３４…データベース保持部
３５、６２、６８…生物情報データベース
３６…データベース検索部
３７…検索結果保持部
３８、８２、１１１、１１３…機能性配列ライブラリー
３９…検索結果解析部
４０…検索結果保持部
４１…スコア化遺跡部
４２…出力制御部
５０…検索インタフェース
５１、５２…遺伝子配列入力部
５３…データベース選択
５４…相同性判定基準
５５…塩基長選択
５６…検索開始ボタン
５７…スコアリング条件選択
６１…BLAST検索
６４…キーワード検索
６７…抽出工程
６９…フレーム確認工程
７０…相補配列対応工程
７１…終止コドン確認工程
７２、８５…変更情報追加工程
８０…天然配列リストとその他の配列リスト
８１…抽出配列確認工程
８３…DB検索
８４、８７、８９、９１…検索結果チェック工程
８６…推定ツール検索
８８…その他の検索
９０…機械学習技術による抽出工程
９２…天然配列リスト再分類工程
９３…人工配列リスト分類工程
１００…特許文献情報
１０１…データ名
１０２…相同性配列
１０３…データベース中のsequence ID
１０４…フレーム情報
１０９…特許本文
１１０…抽出領域
１１２…検索結果
１２１…スコアリング条件選択工程
１２２…頻度別スコア結果リスト作製工程
１２３…キーワードスコアリング実行工程
１２４…キーワードスコア結果作製工程
１２６…結果表示工程
１３０…出力結果
１３１、１３２、１３３、１３４…ポップアップボタン
１３５…ダウンロードボタン
１４０、１６０…プロモーター名
１４１、１６１、１９１…プロモーター領域配列
１４２、１５２、１６５…分類欄
１４３、１５３、１６２、１６６…配列同定方法及び同定に用いたデータベース情報
１４４、１５４、１６７…スコア
１４５、１５５、１６８…チェック箇所
１４６、１５６、１７０…セレクトボタン
１４７、１５７、１６９…全件表示ボタン
１４８、１５８…全結果表示
１５０、１６３…ターミネーター名
１５１、１６４、１９４…ターミネーター領域配列
１８０…クェリとして入力した遺伝子の塩基配列
１８１、１９２…開始コドン
１８２、１９３…終止コドン
１８３…クェリカバー率
１９０…出力配列
【図１】
2020150919.tif 000003
【図２】
2020150919.tif 000004
【図３】
2020150919.tif 000005
【図４】
2020150919.tif 000006
【図５】
2020150919.tif 000007
【図６】
2020150919.tif 000008
【図７】
2020150919.tif 000009
【図８Ａ】
2020150919.tif 000010
【図８Ｂ】
2020150919.tif 000011
【図９】
2020150919.tif 000012
【図１０】
2020150919.tif 000013
【図１１Ａ】
2020150919.tif 000014
【図１１Ｂ】
2020150919.tif 000015
【図１２Ａ】
2020150919.tif 000016
【図１２Ｂ】
2020150919.tif 000017
【図１３】
2020150919.tif 000018
【図１４】
2020150919.tif 000019
【図１５】
2020150919.tif 000020
【配列表】
2020150919000001.app