(57)【特許請求の範囲】
【請求項１】
生物の遺伝子を増幅し、増幅して得たＤＮＡをＴＧＧＥ（Temperature Gradient Gel Electrophoresis）またはＤＧＧＥ（Denaturant Gradient Gel Electrophoresis）により電気泳動させて得られるゲル電気泳動物をゲノムプロフィリング画像として取り込む取込工程と、前記取込工程を介して取り込んだゲノムプロフィリング画像を温度方向および移動度方向に沿ってそれぞれ規格化する規格化工程と、前記規格化工程を介して規格化されたゲノムプロフィリング画像の各バンドに対応して複数の肉付き線分を抽出する肉付き線分抽出工程と、前記肉付き線分抽出工程を介して抽出された複数の肉付き線分の各々について複数の特徴点を抽出する特徴点抽出工程とを含むこと、但し、肉付き線分とは幅のある線分を意味する、並びに前記特徴点抽出工程では、変性前における融解の開始に対応する融解開始点、変性前において移動度が最小値に達する最小移動度点、および変性後の一本鎖状態の移動度に終着する移動度終着点から選択された複数の点を特徴点として抽出すること、および前記ＤＧＧＥによる電気泳動で使用される変性剤が尿素であり、かつ前記ＤＧＧＥにおける温度方向に沿った規格化は尿素濃度を温度換算した還元温度で行うことを特徴とする特徴点抽出方法。
【請求項２】
前記肉付き線分抽出工程は、ゲノムプロフィリング画像の各バンド上の点を局所参照点として指示する局所参照点指示工程と、前記局所参照点指示工程で指示された局所参照点の周囲の点の色の濃淡を判定し、所定の域値内の色の濃さを有する点を新たに局所参照点として設定する二次局所参照点設定工程とを含み、前記二次局所参照点設定工程を介して新たに設定された二次局所参照点にシフトし、該二次局所参照点に対して色の濃さが一定の変動範囲に収まるような新たな局所参照点を設定する工程を繰り返すことにより、各バンドに対応する肉付き線分を抽出することを特徴とする請求項１に記載の方法。
【請求項３】
前記肉付き線分抽出工程を介して抽出された各肉付き線分について、１つの温度点に対して複数の移動度点が対応する状態から、１つの温度点に対して１つの移動度点が対応する状態へスリム化するスリム化工程と、前記スリム化工程を介して得られた各スリム化線分を所定の関数で近似する関数化工程とをさらに含み、前記特徴点抽出工程では、前記関数化工程を介して関数化された各バンド曲線について複数の特徴点を抽出することを特徴とする請求項１乃至２のいずれか１項に記載の方法。
【請求項４】
前記スリム化工程は、１つの温度点に対応する複数の移動度点の重心点を検出する重心点検出工程と、連続する多数の温度点について前記重心点検出工程を繰り返す検出繰返工程と、前記検出繰返工程を介して得られた多数の重心点の連鎖によりスリム化線分を形成するスリム化線分形成工程とを含むことを特徴とする請求項３に記載の方法。
【請求項５】
前記スリム化線分形成工程は、滑らかな曲線からなるスリム化線分を得るために、ある重心点の移動度方向の座標と両側に隣接する重心点の移動度方向の座標との平均座標を、その重心点における移動度方向の座標と設定するスムージング処理工程を含むことを特徴とする請求項４に記載の方法。
【請求項６】
前記関数化工程は、ニュートン・ラプソン（Newton-Raphson）法を用いて各スリム化線分を適当な次数のｎ次斉関数に近似する近似工程を含むことを特徴とする請求項３乃至５のいずれか１項に記載の方法。
【請求項７】
前記規格化工程では、塩基配列のわかっている参照試料であって本来の試料とともに電気泳動させた内部参照試料の特徴点および実験条件パラメータなどに基づいて、移動度方向に沿って０と１との間で無次元化することを特徴とする請求項１乃至６のいずれか１項に記載の方法。
【請求項８】
前記規格化工程に先立って、前記取込工程を介して取り込んだゲノムプロフィリング画像のゆがみを補正するゆがみ補正工程をさらに含むことを特徴とする請求項１乃至７のいずれか１項に記載の方法。
【請求項９】
前記ゆがみ補正工程では、ゲルの４隅を定義し、定義された４隅に基づいて補正処理を行うことを特徴とする請求項８に記載の方法。
【請求項１０】
前記ゆがみ補正工程では、塩基配列のわかっている参照試料であって本来の試料とともに電気泳動させた内部参照試料の特徴点を認識し、認識した特徴点に基づいて補正処理を行うことを特徴とする請求項８に記載の方法。
【請求項１１】
特徴点群に関する登録情報を参照するために比較対象ゲノムの選択を行う選択工程と、前記選択工程を介して選択された比較対象ゲノムの特徴点群に関する登録情報を参照し、請求項１乃至１０のいずれか１項に記載の方法により試料生物について抽出された特徴点群から対応特徴点群を決定する特徴点群対応付け工程と、前記比較対象ゲノムの特徴点群と前記試料生物の対応特徴点群との比較に基づいて、前記試料生物の種の同定を行う種同定工程とを含むことを特徴とする種同定方法。
【請求項１２】
前記種同定工程において、比較対象ゲノムのｎ個の特徴点Ｐ0iの位置ベクトルをＶ0iとし、試料生物のｎ個の特徴点Ｐ1iの位置ベクトルをＶ1iとしたとき、以下の式（１）で規定されるPaSSの値に基づいて、前記試料生物の種の同定を行うことを特徴とする請求項１１に記載の方法。
PaSS＝１−｛Σγ（ｉ）｝／ｎ  （１）
ここで、Σはｉ＝１〜ｎまでの総和記号であり、γ（ｉ）は以下の式（２）で表される。
γ（ｉ）＝｜Ｖ1i−Ｖ0i｜／（｜Ｖ1i｜＋｜Ｖ0i｜）  （２）
【請求項１３】
前記種同定工程において、比較対象ゲノムのｎ個の特徴点Ｐ0iの温度をＴ0iとし、試料生物のｎ個の特徴点Ｐ1iの温度をＴ1iとし、比較対象ゲノムのｎ個の特徴点Ｐ0iの移動度をＭ0iとし、試料生物のｎ個の特徴点Ｐ1iの移動度をＭ1iとし、特徴点の現れる温度範囲幅に相当する温度規格化因子をＴｗとし、移動度変化と温度変化の両方の「重み」の相対値をαとしたとき、以下の式（４）で規定されるPaSSの値に基づいて、前記試料生物の種の同定を行うことを特徴とする請求項１１に記載の方法。
PaSS＝１−｛Σｒ（ｉ）｝／ｎ  （４）
ここで、Σはｉ＝１〜ｎまでの総和記号であり、ｒ（ｉ）は以下の式（５）で表される。
【数１】
0004527841.tif 000004
【請求項１４】
前記PaSSを用いて、以下の式（３）で規定されるゲノム準距離ｄsの値に基づいて、前記試料生物の種の同定を行うことを特徴とする請求項１２または１３に記載の方法。
ｄs＝（１−PaSS）／PaSS  （３）
【請求項１５】
特徴点群に関する登録情報を参照するために比較対象ゲノムの選択を行う選択工程と、前記選択工程を介して選択された比較対象ゲノムの特徴点群に関する登録情報を参照し、請求項１乃至１０のいずれか１項に記載の方法により試料生物について抽出された特徴点群から対応特徴点群を決定する特徴点群対応付け工程と、前記比較対象ゲノムの特徴点群と前記試料生物の対応特徴点群との比較に基づいて、前記試料生物について同一種の個体間における類縁性の同定を行う類縁性同定工程とを含むことを特徴とする種同定方法。
【請求項１６】
前記類縁性同定工程において、比較対象ゲノムのｎ個の特徴点Ｐ0iの位置ベクトルをＶ0iとし、試料生物のｎ個の特徴点Ｐ1iの位置ベクトルをＶ1iとしたとき、以下の式（１）で規定されるPaSSの値に基づいて、前記試料生物について同一種の個体間における類縁性の同定を行うことを特徴とする請求項１５に記載の方法。
PaSS＝１−｛Σγ（ｉ）｝／ｎ  （１）
ここで、Σはｉ＝１〜ｎまでの総和記号であり、γ（ｉ）は以下の式（２）で表される。
γ（ｉ）＝｜Ｖ1i−Ｖ0i｜／（｜Ｖ1i｜＋｜Ｖ0i｜）  （２）
【請求項１７】
前記類縁性同定工程において、比較対象ゲノムのｎ個の特徴点Ｐ0iの温度をＴ0iとし、試料生物のｎ個の特徴点Ｐ1iの温度をＴ1iとし、比較対象ゲノムのｎ個の特徴点Ｐ0iの移動度をＭ0iとし、試料生物のｎ個の特徴点Ｐ1iの移動度をＭ1iとし、特徴点の現れる温度範囲幅に相当する温度規格化因子をＴｗとし、移動度変化と温度変化の両方の「重み」の相対値をαとしたとき、以下の式（４）で規定されるPaSSの値に基づいて、前記試料生物について同一種の個体間における類縁性の同定を行うことを特徴とする請求項１５に記載の方法。
PaSS＝１−｛Σｒ（ｉ）｝／ｎ  （４）
ここで、Σはｉ＝１〜ｎまでの総和記号であり、ｒ（ｉ）は以下の式（５）で表される。
【数２】
0004527841.tif 000005
【請求項１８】
前記PaSSを用いて、以下の式（３）で規定されるゲノム準距離ｄsの値に基づいて、前記試料生物について同一種の個体間における類縁性の同定を行うことを特徴とする請求項１６または１７に記載の方法。
ｄs＝（１−PaSS）／PaSS  （３）
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、ゲノムプロフィリング画像から特徴点を抽出する方法、並びに該方法で得られた特徴点群を用いた遺伝子型による種同定方法および類縁性同定方法に関するものである。
【０００２】
【従来の技術】
たとえば、所定の条件においてランダムＰＣＲで微生物の遺伝子を増幅し、増幅したＤＮＡをＴＧＧＥまたはＤＧＧＥにより電気泳動させることにより、電気泳動物が得られる。この電気泳動物を撮像した画像は、ゲル電気泳動画像あるいはゲノムプロフィリング画像と呼ばれている。
【０００３】
【発明が解決しようとする課題】
２つの微生物が同一種に属する場合、２つの微生物についてそれぞれ得られる２つのゲノムプロフィリング画像の全体的な形態が非常に良く似ていることが知られている。しかしながら、ゲノムプロフィリング画像の全体的な形態に基づく定性的な手法では、遺伝子型による微生物の種の同定を正確に行うことができなかった。
【０００４】
本発明は、前述の課題に鑑みてなされたものであり、遺伝子型による微生物の種の同定を定量的な手法で正確に且つ簡便に行うことを可能にするための中間的な処理方法として、ゲノムプロフィリング画像から特徴点を抽出する方法を提供することを目的とする。
【０００５】
また、本発明は、上述の特徴点抽出方法で得られた特徴点群を用いた定量的な手法で、微生物の種の同定および同一種の微生物における個体間の類縁性の同定を正確に且つ簡便に行うことのできる、遺伝子型による種同定方法および類縁性同定方法を提供することを目的とする。
【０００６】
【課題を解決するための手段】
前記課題を解決するために、本発明の第1発明では、生物の遺伝子を増幅し、増幅して得たＤＮＡをＴＧＧＥ（Temperature Gradient Gel Electrophoresis）またはＤＧＧＥ（Denaturant Gradient Gel Electrophoresis）により電気泳動させて得られるゲル電気泳動物をゲノムプロフィリング画像として取り込む取込工程と、前記取込工程を介して取り込んだゲノムプロフィリング画像を温度方向および移動度方向に沿ってそれぞれ規格化する規格化工程と、前記規格化工程を介して規格化されたゲノムプロフィリング画像の各バンドに対応して複数の肉付き線分を抽出する肉付き線分抽出工程と、前記肉付き線分抽出工程を介して抽出された複数の肉付き線分の各々について複数の特徴点を抽出する特徴点抽出工程とを含むこと、並びに前記特徴点抽出工程では、変性前における融解の開始に対応する融解開始点、変性前において移動度が最小値に達する最小移動度点、および変性後の一本鎖状態の移動度に終着する移動度終着点から選択された複数の点を特徴点として抽出すること、および前記ＤＧＧＥによる電気泳動で使用される変性剤が尿素であり、かつ前記ＤＧＧＥにおける温度方向に沿った規格化は尿素濃度を温度換算した還元温度で行うことを特徴とする特徴点抽出方法を提供する。
【０００８】
また、第１発明の好ましい態様によれば、前記肉付き線分抽出工程は、
ゲノムプロフィリング画像の各バンド上の点を局所参照点として指示する局所参照点指示工程と、
前記局所参照点指示工程で指示された局所参照点の周囲の点の色の濃淡を判定し、所定の域値内の色の濃さを有する点を新たに局所参照点として設定する二次局所参照点設定工程とを含み、
前記二次局所参照点設定工程を介して新たに設定された二次局所参照点にシフトし、該二次局所参照点に対して色の濃さが一定の変動範囲に収まるような新たな局所参照点を設定する工程を繰り返すことにより、各バンドに対応する肉付き線分を抽出する。
【０００９】
さらに、第１発明の好ましい態様によれば、前記肉付き線分抽出工程を介して抽出された各肉付き線分について、１つの温度点に対して複数の移動度点が対応する状態から、１つの温度点に対して１つの移動度点が対応する状態へスリム化するスリム化工程と、前記スリム化工程を介して得られた各スリム化線分を所定の関数で近似する関数化工程とをさらに含み、前記特徴点抽出工程では、前記関数化工程を介して関数化された各バンド曲線について複数の特徴点を抽出する。
【００１０】
この場合、前記スリム化工程は、１つの温度点に対応する複数の移動度点の重心点を検出する重心点検出工程と、連続する多数の温度点について前記重心点検出工程を繰り返す検出繰返工程と、前記検出繰返工程を介して得られた多数の重心点の連鎖によりスリム化線分を形成するスリム化線分形成工程とを含むことが好ましい。また、この場合、前記スリム化線分形成工程は、滑らかな曲線からなるスリム化線分を得るために、ある重心点の移動度方向の座標と両側に隣接する重心点の移動度方向の座標との平均座標を、その重心点における移動度方向の座標と設定するスムージング処理工程を含むことが好ましい。
【００１１】
また、第１発明の好ましい態様によれば、前記関数化工程は、ニュートン・ラプソン（Newton-Raphson）法を用いて各スリム化線分を適当な次数のｎ次斉関数に近似する近似工程を含む。
【００１２】
さらに、第１発明の好ましい態様によれば、前記規格化工程では、塩基配列のわかっている参照試料であって本来の試料とともに電気泳動させた内部参照試料の特徴点および実験条件パラメータなどに基づいて、移動度方向に沿って０と１との間で無次元化するとともに、尿素を温度換算した還元温度で温度方向に沿って規格化する。
【００１３】
また、第１発明の好ましい態様によれば、前記規格化工程に先立って、前記取込工程を介して取り込んだゲノムプロフィリング画像のゆがみを補正するゆがみ補正工程をさらに含む。この場合、前記ゆがみ補正工程では、ゲルの４隅を定義し、定義された４隅に基づいて補正処理を行うことが好ましい。あるいは、前記ゆがみ補正工程では、塩基配列のわかっている参照試料であって本来の試料とともに電気泳動させた内部参照試料の特徴点を認識し、認識した特徴点に基づいて補正処理を行うことが好ましい。
【００１４】
本発明の第２発明では、特徴点群に関する登録情報を参照するために比較対象ゲノムの選択を行う選択工程と、前記選択工程を介して選択された比較対象ゲノムの特徴点群に関する登録情報を参照し、第１発明の方法により試料生物について抽出された特徴点群から対応特徴点群を決定する特徴点群対応付け工程と、前記比較対象ゲノムの特徴点群と前記試料生物の対応特徴点群との比較に基づいて、前記試料生物の種の同定を行う種同定工程とを含むことを特徴とする種同定方法を提供する。
【００１５】
本発明の第３発明では、特徴点群に関する登録情報を参照するために比較対象ゲノムの選択を行う選択工程と、前記選択工程を介して選択された比較対象ゲノムの特徴点群に関する登録情報を参照し、第１発明の方法により試料生物について抽出された特徴点群から対応特徴点群を決定する特徴点群対応付け工程と、前記比較対象ゲノムの特徴点群と前記試料生物の対応特徴点群との比較に基づいて、前記試料生物について同一種の個体間における類縁性の同定を行う類縁性同定工程とを含むことを特徴とする種同定方法を提供する。
【００１６】
第２発明および第３発明において、比較対象ゲノムのｎ個の特徴点Ｐ0iの位置ベクトルをＶ0iとし、試料生物のｎ個の特徴点Ｐ1iの位置ベクトルをＶ1iとしたとき、以下の式（１）で規定されるPaSSの値に基づいて、前記試料生物の種の同定または前記試料生物について同一種の個体間における類縁性の同定を行うことが好ましい。
PaSS＝１−｛Σγ（ｉ）｝／ｎ      （１）
ここで、Σはｉ＝１〜ｎまでの総和記号であり、γ（ｉ）は以下の式（２）で表される。
γ（ｉ）＝｜Ｖ1i−Ｖ0i｜／（｜Ｖ1i｜＋｜Ｖ0i｜）    （２）
【００１７】
あるいは、第２発明および第３発明において、比較対象ゲノムのｎ個の特徴点Ｐ0iの温度をＴ0iとし、試料生物のｎ個の特徴点Ｐ1iの温度をＴ1iとし、比較対象ゲノムのｎ個の特徴点Ｐ0iの移動度をＭ0iとし、試料生物のｎ個の特徴点Ｐ1iの移動度をＭ1iとし、特徴点の現れる温度範囲幅に相当する温度規格化因子をＴｗとし、移動度変化と温度変化の両方の「重み」の相対値をαとしたとき、以下の式（４）で規定されるPaSSの値に基づいて、前記試料生物の種の同定または前記試料生物について同一種の個体間における類縁性の同定を行うことが好ましい。
PaSS＝１−｛Σｒ（ｉ）｝／ｎ      （４）
ここで、Σはｉ＝１〜ｎまでの総和記号であり、ｒ（ｉ）は以下の式（５）で表される。
【数３】
0004527841.tif 000002
【００１８】
この場合、前記PaSSを用いて、以下の式（３）で規定されるゲノム準距離ｄsの値に基づいて、前記試料生物の種の同定または前記試料生物について同一種の個体間における類縁性の同定を行うことがさらに好ましい。
ｄs＝（１−PaSS）／PaSS      （３）
【００１９】
【発明の実施の形態】
本発明の実施形態を、添付図面に基づいて説明する。
図１は、本発明の実施形態にかかる特徴点の抽出方法、すなわちゲノムプロフィリング画像から特徴点を抽出する方法について、その全体的な処理フローを示すフローチャートである。
所定の条件においてランダムＰＣＲ（Random Polymerase Chain Reaction）できればスペルを記入してください）で微生物の遺伝子を増幅し、増幅したＤＮＡをＴＧＧＥ（Temperature Gradient Gel Electrophoresis）またはＤＧＧＥ（Denaturant Gradient Gel Electrophoresis）により電気泳動させることにより、電気泳動物が得られる。
【００２０】
本実施形態では、この電気泳動物を、適当な撮像手段を介して、情報処理手段の適当な電子記録媒体にカラー画像やグレースケール画像として取り込む。具体的には、電気泳動物を、たとえばデジタルカメラを介して撮像した後に、コンピュータのハードディスクなどに、ゲル電気泳動画像すなわちゲノムプロフィリング画像として取り込む（Ｓ１０１）。あるいは、電気泳動物を、アナログカメラを介して撮像し、次いでスキャナーを介してイメージスキャニングした後に、ゲノムプロフィリング画像として取り込む。または、蛍光イメージャーを用い、蛍光性バンドとして得られたゲル中の試料を撮像する。
【００２１】
図２は、典型的なゲノムプロフィリング画像を模式的に示す図である。
図１のゲノムプロフィリング画像は、白黒で模式的に示されている。しかしながら、銀染色法によって得られた実際のゲノムプロフィリング画像の背景（図１中白い部分）は、たとえば濃い黄色である。また、図１のゲノムプロフィリング画像において黒い部分はバンドと呼ばれるが、その色は実際には黒に近い色と薄茶色に近い色との間で変化している。ここで、黒に近い色のバンドが明瞭に現れているバンドである。また、蛍光染色法では、色素ごとに励起・発光波長が異なり、それぞれ異なる色素で染め分けられたバンドを蛍光イメージャーで取り込んだ後に、コンピュータにより多色表示する。
【００２２】
本実施形態では、取り込まれたゲノムプロフィリング画像に対して、必要に応じてゆがみの補正処理を行う（Ｓ１０２）。そして、必要に応じてゆがみの補正処理がなされたゲノムプロフィリング画像に対して、規格化処理を行う（Ｓ１０３）。ここで、ゆがみの補正処理とは、取り込まれたゲノムプロフィリング画像が例えば台形状に歪んでいた場合に、矩形状に形状補正する画像処理である。取り込まれたゲノムプロフィリング画像のゆがみが通常の誤差範囲であれば、すなわち許容範囲内であれば、ゆがみの補正処理を省略することができる。
【００２３】
具体的に、ゆがみの補正処理に際して、ゲルの４隅を定義し、定義された４隅に基づいて補正処理を行う隅４点補正方式を用いることができる。この場合、実験者がゲルの４隅を定義する手動補正方式と、ゲルの４隅を情報処理装置で自動判別して処理する自動補正方式とがある。
【００２４】
また、ゆがみの補正処理に際して、内部参照試料（塩基配列のわかっている参照試料であって本来の試料とともに電気泳動させたもの）の特徴点を認識し、認識した特徴点に基づいて補正処理を行う内部参照特徴点方式を用いることができる。この場合、実験者が特徴点を目視で認識する手動補正方式と、特徴点を情報処理装置で自動認識して処理する自動補正方式とがある。なお、たとえば２種類の参照試料を使用し、各参照試料からそれぞれ２つ以上の特徴点（合計で４つ以上の特徴点）を参照特徴点として得ることが好ましい。
【００２５】
一方、規格化処理とは、内部参照試料の特徴点および実験条件パラメータなどに基づいて縦軸（Ｙ軸：移動度）方向および横軸（Ｘ軸：温度）方向を規格化する画像処理である。規格化処理により、移動度Ｍは０から１までの間で無次元化される。また、温度Ｔとして、たとえば尿素を温度換算した温度すなわち還元温度が適用される。
【００２６】
次いで、必要に応じてゆがみ補正され且つ規格化処理されたゲノムプロフィリング画像の各バンドに対して、肉付き線分抽出処理を行う（Ｓ１０４）。具体的には、肉付き線分抽出処理の開始に際して、実験者がバンド上の１つの点を局所参照点として指示する。情報処理装置では、指示された局所参照点の周囲の点（ピクセル：画素）の色の濃淡を判定し、所定のスレッショルド（閾値）の範囲内（例えば変動１０％）の色の濃さを有する点を新たに二次局所参照点として設定する。
【００２７】
そして、新たに設定された二次局所参照点にシフトし、この二次局所参照点に対して色の濃さが一定の変動範囲に収まるように、その周囲の点の色の濃淡に基づいて、さらに新たな局所参照点を設定する。こうして、順次設定される局所参照点についてその周囲の点に対する色の濃淡判定および局所参照点の設定を繰り返す。その結果、いわゆる連形成参照点シフト方式により、ゲノムプロフィリング画像から各バンドに対応する肉付き線分が抽出される。
【００２８】
なお、変性前から変性後にかけて、すなわち低温領域から高温領域の全体に亘ってバンドが連続的である場合、上述の１点指定方式により所望の肉付き線分が自動抽出される。しかしながら、変性前と変性後とでバンドが不連続になる場合、たとえば変性前のバンド上の１点を指定しただけでは、変性前のバンドに対応する肉付き線分は自動抽出されるが、変性後のバンドに対応する肉付き線分は抽出されないことになる。したがって、変性前と変性後とでバンドが不連続になる場合には、変性前のバンド上において１つの点を指定し且つ変性後のバンド上においてもう１つの点を指定する２点指定方式を用いることにより、変性前のバンドに対応する肉付き線分および変性後のバンドに対応する肉付き線分の双方を自動抽出することができる。
【００２９】
ところで、上述の説明では、連形成参照点シフト方式によりゲノムプロフィリング画像の各バンドから肉付き線分を抽出している。しかしながら、たとえばグレースケールで取り込んだゲノムプロフィリング画像を所定のスレッショルドに基づいて白黒画像に二値化する二値化方式により、ゲノムプロフィリング画像から各バンドに対応する肉付き線分を抽出することもできる。また、ゲノムプロフィリング画像の各バンドには、いくつかの典型的な形状を有するパターンが部分的に現れる。そこで、多数のパターンセグメントから選択された複数のパターンセグメントを組み合わせて、いわゆるパターンイメージ合成法により、ゲノムプロフィリング画像から各バンドに対応する肉付き線分を抽出することもできる。
【００３０】
次いで、ゲノムプロフィリング画像から各バンドに対応して抽出された各肉付き線分に対して、スリム化処理を行う（Ｓ１０５）。図３に示すように、肉付き線分３１では、横軸（Ｘ方向）に沿った１つの温度点Ｔｉに対して縦軸（Ｙ方向）に沿って複数の移動度点Ｍｉ（ｉ＝１〜ｎ）が対応し、いわゆる多価関数となっている。スリム化処理とは、横軸に沿った１つの温度点に対して１つの移動度点が対応するように、肉付き線分をいわゆる１価関数のスリム化線分へ変換する処理である。
【００３１】
肉付き線分のスリム化処理に際しては、図３に示すように、横軸に沿った１つの温度点Ｔｉに対応する複数の移動度点Ｍｉの重心点Ｍｇを求める。ここで、１つの温度点に対応する複数の移動度点がＹ方向に沿って連続して配列されている場合には、複数の移動度点の重心点を求めることは、連続する複数の移動度点の幾何学的中心点を求めることに他ならない。一方、１つの温度点に対応する複数の移動度点がＹ方向に沿って断続的に配列されている場合には、複数の移動度点の重心点はその幾何学的中心点とは異なるものとなるのが一般的である。
【００３２】
複数の移動度点Ｍｉの重心点Ｍｇを求める処理は、横軸に沿った多数の温度点Ｔｉ（ｉ＝１〜ｍ）に対して繰り返される。こうして、図４に示すように、横軸に沿った多数の温度点Ｔｉ（ｉ＝１〜ｍ）に対応して多数の重心点Ｍgi（ｉ＝１〜ｍ）が得られる。ここで、得られた多数の重心点Ｍgiを単に連結すると、滑らかなスリム化線分を得ることができない場合がある。その場合、たとえばランニングアベレージ法を用いて、ある重心点のＹ座標と両側に隣接する重心点のＹ座標との平均座標をその重心点におけるＹ座標と設定するスムージング処理により、多数の重心点からなる滑らかな連鎖を、ひいては滑らかな曲線からなるスリム化線分が得られる。
【００３３】
こうして、滑らかな曲線状に得られたスリム化線分に対して、関数化処理を行う（Ｓ１０６）。関数化処理では、たとえばニュートン・ラプソン（Newton-Raphson）法を用いて、各スリム化線分を適当な次数（たとえば１０次）のｎ次斉関数に近似する。あるいは、各スリム化線分を三角関数近似することもできる。
【００３４】
次いで、スリム化処理および関数化処理を介して得られた各バンド曲線について特徴点抽出処理を行う（Ｓ１０７）。特徴点抽出処理では、特徴点として、たとえば融解開始点Ｐini、最小移動度点Ｐmin、移動度終着点Ｐendを抽出する。ここで、融解開始点Ｐiniは、変性前における融解の開始に対応するバンド曲線上の点である。また、最小移動度点Ｐminは、変性前において移動度が最小値に達するバンド曲線上の点である。さらに、移動度終着点Ｐendは、変性後の一本鎖状態の移動度に終着するバンド曲線上の点である。
【００３５】
図５は、第１の典型的なバンド曲線の形態を示す図である。
図５に示すように、典型的な第１バンド曲線５１では、融解が開始する温度Ｔiniよりも低い温度において移動度Ｍは温度Ｔに関してほぼ線形的に増大する。融解開始点Ｐiniに対応する温度Ｔiniに達して融解が開始すると、移動度Ｍは温度Ｔに関してほぼ線形的に減少した後に、その減少率が漸減しながら、やがて最小移動度点Ｐminに達する。
【００３６】
最小移動度点Ｐminに対応する温度Ｔminに達した後、温度の上昇に対して移動度Ｍが急増し、やがて変性後の一本鎖状態の移動度すなわち移動度終着点Ｐendに達する。移動度終着点Ｐendに対応する温度Ｔendよりも高い温度では、移動度Ｍは再び温度Ｔに関してほぼ線形的に増大する。これは温度上昇に伴う粘度の低下のためであって、形態の変化はほとんどないと考えられている。なお、最小移動度点Ｐminが不連続点となり、最小移動度点Ｐminに対応する温度Ｔminと移動度終着点Ｐendに対応する温度Ｔendよりもわずかに低い温度との間でバンド曲線が不連続になることが多い。
【００３７】
図６は、第２の典型的なバンド曲線の形態を示す図である。
図６に示すように、典型的な第２バンド曲線６１では、第１バンド曲線５１と同様に、融解が開始する温度Ｔiniよりも低い温度において移動度Ｍは温度Ｔに関してほぼ線形的に増大し、融解開始点Ｐiniに対応する温度Ｔiniに達して融解が開始すると移動度Ｍは温度Ｔに関してほぼ線形的に減少する。しかしながら、第１バンド曲線５１の場合とは異なり、やがて最小移動度点Ｐminに達した後、最小移動度点Ｐminに対応する温度Ｔminよりも高い温度では、移動度Ｍは再び温度Ｔに関してほぼ線形的に増大する。すなわち、第２バンド曲線６１では、最小移動度点Ｐminと移動度終着点Ｐendとが一致し、バンド曲線が不連続になることはない。
【００３８】
本実施形態では、関数化処理を介して得られたバンド曲線の一次微分、二次微分、その不連続点などの情報に基づいて、たとえば融解開始点Ｐini、最小移動度点Ｐmin、移動度終着点Ｐendを特徴点として抽出する。
具体的には、融解開始点Ｐiniは、２つの直線部分５２（６２）と５３（６３）とによって挟まれた曲線部分の中間点として抽出される。したがって、融解開始点Ｐiniは、バンド曲線５１（６１）の一次微分値および二次微分値に基づいて求められる。
【００３９】
また、最小移動度点Ｐminは、第１バンド曲線５１の場合、たとえば曲線上の不連続点として抽出される。一方、第２バンド曲線６１の場合、最小移動度点Ｐminは、移動度終着点Ｐendと同じ点として後述のように求められる。
さらに、移動度終着点Ｐendは、直線部分５４（６４）を−Ｘ方向に沿って見ると直線部分から離れる点として、バンド曲線５１（６１）の一次微分値および二次微分値に基づいて求められる。
【００４０】
以上の特徴点抽出処理は、明瞭なバンドすなわち明瞭に線分化されたバンド曲線から順に行われ、全体として所定数の特徴点が得られるまで行われる。なお、試料微生物の種が予め分かっており、特徴点を抽出すべき指定バンドが自明である場合には、それらの指定バンドについて特徴点抽出処理を行う。図７では、３つのバンド曲線から抽出された合計８個の特徴点Ｐｊ（ｊ＝１〜８）すなわち特徴点群を模式的に示している。実際には、３つよりも多くのバンド曲線から８個よりも多くの特徴点を抽出する。
【００４１】
なお、上述の説明では、スリム化処理および関数化処理を介して得られたバンド曲線から特徴点を抽出しているが、スリム化処理および関数化処理を行うことなく、各肉付け線分から直接的に特徴点を抽出することもできる。
【００４２】
図８は、本発明の実施形態にかかる種同定方法、すなわち本発明の特徴点抽出方法で得られた特徴点群を用いた遺伝子型による種同定方法について、その全体的な処理フローを示すフローチャートである。
本実施形態の種同定方法では、比較対象ゲノムの選択を行い、その特徴点群（spiddo: species identification dots）に関するデータベース上の登録情報を参照する（Ｓ２０１）。
【００４３】
まず、試料微生物の種が予め分かっている場合、その種を比較対象とし、その特徴点群に関する登録情報を参照する。
また、試料微生物の種が未知である場合、ゲノムプロフィリング画像の全体的な形態に基づいて予め選定された代表種（たとえば数十種）を仮の比較対象とし、その特徴点群に関する登録情報を参照する。
【００４４】
あるいは、試料微生物の種が未知である場合、上述の特徴点抽出処理により得られた各特徴点について類似する特徴点（たとえば図７において原点Ｏから各特徴点までの距離の差異が５％以内）を有する種を順次リストアップする。そして、リストアップされた種を仮の比較対象とし、その特徴点群に関する登録情報を参照する。
【００４５】
次いで、比較対象ゲノムの特徴点群に関する登録情報を参照しながら、上述の特徴点抽出処理で得られた特徴点群から対応特徴点（たとえば１０点）を決定する（Ｓ２０２）。
第１の決定方法として、実験者指定方式がある。実験者指定方式では、比較対象ゲノムの特徴点群に関する登録情報を参照しながら、実験者がコンピュータディスプレー上で対応特徴点を１つづつ手動で指定する。ここで、対応特徴点とは、比較対象ゲノムの特徴点群の各特徴点に対して位置的に対応する特徴点である。
【００４６】
第２の決定方法として、自動割当方式がある。自動割当方式では、比較対象ゲノムの特徴点群の各特徴点が属する座標ゾーン（ＸＹ座標において二次元的な所定の面積を占める領域）を決定し、その座標ゾーンに属する試料特徴点を対応特徴点として順次自動的に割り当てる。
【００４７】
第３の決定方法として、最適対応付け自動計算法がある。最適対応付け自動計算法では、対応特徴点数を適当に設定し、比較対象ゲノムの特徴点群および抽出された試料特徴点からそれぞれ対応特徴点数だけ任意に選択する。そして、対応付けを組み合わせ論的に行い、後に定義するPaSSが最も大きくなる組み合わせに基づいて対応特徴点を自動的に求める。
【００４８】
こうして、図９に示すように、比較対象ゲノムのｎ個の特徴点Ｐ0i（ｉ＝１〜ｎ）に対応して、試料微生物についてｎ個の特徴点Ｐ1i（ｉ＝１〜ｎ）が得られる。ここで、各特徴点Ｐｉ（比較対象ゲノムのｎ個の特徴点Ｐ0i および試料微生物についてｎ個の特徴点Ｐ1i）は、位置ベクトルＶｉ（比較対象ゲノムのｎ個の特徴点の位置ベクトルＶ0iおよび試料微生物についてｎ個の特徴点の位置ベクトルＶ1i）を有する。なお、位置ベクトルＶｉの原点として、たとえばＭ＝０でＴ＝最小温度の点Ｏを採用することができる。
【００４９】
本実施形態では、各特徴点Ｐｉの位置ベクトルＶｉを用いて、以下の式（１）で規定されるPaSSという概念を導入する。
PaSS＝１−｛Σγ（ｉ）｝／ｎ      （１）
ここで、Σはｉ＝１〜ｎまでの総和記号であり、γ（ｉ）は以下の式（２）で表される。
γ（ｉ）＝｜Ｖ1i−Ｖ0i｜／（｜Ｖ1i｜＋｜Ｖ0i｜）    （２）
【００５０】
なお、上述の説明では、ベクトル（依存）型の計算法にしたがってPaSSの値を求めているが、以下に示すようにスカラー型の計算法に基づいてPaSSの値を求めることもできる。
PaSS＝１−｛Σｒ（ｉ）｝／ｎ      （４）
ここで、Σはｉ＝１〜ｎまでの総和記号であり、ｒ（ｉ）は以下の式（５）で表される。
【数４】
0004527841.tif 000003
【００５１】
ここで、Ｔ0iは比較対象ゲノムのｎ個の特徴点Ｐ0iの温度であり、Ｔ1iは試料生物のｎ個の特徴点Ｐ1iの温度である。また、Ｍ0iは比較対象ゲノムのｎ個の特徴点Ｐ0iの移動度であり、Ｍ1iは試料生物のｎ個の特徴点Ｐ1iの移動度である。さらに、Ｔｗは、特徴点の現れる温度範囲幅（たとえば約６０°Ｃ）に相当する温度規格化因子である。また、αは、移動度変化と温度変化の両方の「重み」の相対値であって、通常は１である。ところで、式（５）の分母部分（α2＋１2）1/2は、第２の規格化因子（ｒ（ｉ）の値を１以下にするための）を構成している。
【００５２】
また、式（１）または（４）に基づいて得られるPaSSを用いて、以下の式（３）で規定されるゲノム準距離ｄsという概念を導入する。
ｄs＝（１−PaSS）／PaSS      （３）
一般に、試料微生物の種と比較対象ゲノムの種とが同じ場合、PaSSの値（スコア）は１に近くなり、ゲノム準距離ｄsの値は０に近くなる。
【００５３】
試料微生物の種が未知である場合、各比較対象ゲノムに対してPaSSおよび必要に応じてゲノム準距離ｄsを計算する（Ｓ２０３）。そして、PaSSの値が標準値（たとえば０．９６）を上回るまで、あるいはゲノム準距離ｄsの値が標準値（たとえば０．０４）を下回るまで、他の対象ゲノムとの比較を繰り返す（Ｓ２０４）。こうして、PaSSの値またはゲノム準距離ｄsの値に基づいて試料微生物の種の同定を行う（Ｓ２０５）。
【００５４】
一方、試料微生物の種が既知である場合、PaSSの値またはゲノム準距離ｄsの値に基づいて同一種における個体間の類縁性を判定することができる。すなわち、同一種に属する２つの試料微生物についてPaSSまたはゲノム準距離ｄsを求め、２つのPaSSの値が十分に近い場合または２つのゲノム準距離ｄsの値が十分に近い場合には２つの試料微生物が同一種であってその個体間の類縁性が高いものと判定される。
【００５５】
さらに、ある一定のゲノム準距離以内にある種や個体の集合を近隣情報としてゲノム配列空間にプロットすることにより保存することが好ましい。
こうして、ゲノムプロフィリング画像、特徴点群、判断情報、近隣情報などは、適当なデータベースに登録され、随時利用可能に構成されることが好ましい。
【００５６】
以上のように、本実施形態では、変性前における融解の開始に対応する融解開始点、変性前において移動度が最小値に達する最小移動度点、および変性後の一本鎖状態の移動度に終着する移動度終着点を特徴点としてゲノムプロフィリング画像から抽出している。その結果、得られた特徴点群を用いた定量的な手法で、微生物の種の同定および同一種の微生物における個体間の類縁性の同定を正確に且つ簡便に行うことができる。
【００５７】
なお、上述の実施形態では、遺伝子型による微生物の種同定方法および類縁性同定方法について説明している。しかしながら、微生物に限定されることなく、一般の生物の種同定方法および類縁性同定方法に本発明を適用することができることはいうまでもない。
【００５８】
【発明の効果】
以上説明したように、本発明のゲノムプロフィリング画像から特徴点を抽出する方法では、遺伝子型による微生物の種の同定を定量的な手法で正確に且つ簡便に行うことを可能にするための中間的な処理を実現することができる。
【００５９】
また、本発明の遺伝子型による種同定方法および類縁性同定方法では、上述の特徴点抽出方法で得られた特徴点群を用いた定量的な手法で、微生物の種の同定および同一種の微生物における個体間の類縁性の同定を正確に且つ簡便に行うことができる。
【図面の簡単な説明】
【図１】本発明の実施形態にかかる特徴点の抽出方法、すなわちゲノムプロフィリング画像から特徴点を抽出する方法について、その全体的な処理フローを示すフローチャートである。
【図２】典型的なゲノムプロフィリング画像を模式的に示す図である。
【図３】肉付き線分では、横軸に沿った１つの温度点に対して縦軸に沿って複数の移動度点が対応し、いわゆる多価関数となっていることを示す図である。
【図４】横軸に沿った多数の温度点に対応して多数の重心点が得られる様子を示す図である。
【図５】第１の典型的なバンド曲線の形態を示す図である。
【図６】第２の典型的なバンド曲線の形態を示す図である。
【図７】複数のバンド曲線からそれぞれ抽出された特徴点群を模式的に示す図である。
【図８】本発明の実施形態にかかる種同定方法、すなわち本発明の特徴点抽出方法で得られた特徴点群を用いた遺伝子型による種同定方法について、その全体的な処理フローを示すフローチャートである。
【図９】比較対象ゲノムにかかる特徴点群と試料微生物にかかる特徴点群とを模式的に示す図である。
【符号の説明】
Ｔ  温度
Ｍ  移動度
Ｏ  原点
Ｐini  融解開始点
Ｐmin  最小移動度点
Ｐend  移動度終着点
３１  肉付き線分
５１  第１バンド曲線
６１  第２バンド曲線
【図１】
0004527841.tif 000006
【図２】
0004527841.tif 000007
【図３】
0004527841.tif 000008
【図４】
0004527841.tif 000009
【図５】
0004527841.tif 000010
【図６】
0004527841.tif 000011
【図７】
0004527841.tif 000012
【図８】
0004527841.tif 000013
【図９】
0004527841.tif 000014