(57)【要約】【課題】遺伝子の変異の進行状況を容易に把握する。【解決手段】受付部４１が、対象被験者及び対象遺伝子の指定を受け付け、判定部４２が、対象被験者の検体試料について、対象遺伝子に関する融解温度測定により得られ、かつ検体試料の採取時が各々異なる複数の融解曲線データが融解曲線ＤＢ５０に存在するか否かを判定する。複数の融解曲線データが存在すると判定された場合には、表示制御部１４が、複数の融解曲線データの各々が示す融解曲線の各々を、１つの座標系内に表示する。【選択図】図２
2016208949.tif 000002
【特許請求の範囲】
【請求項１】
対象被験者の検体試料について、対象遺伝子に関する融解温度測定により得られ、かつ前記検体試料の採取時が各々異なる複数の融解曲線データが存在するか否かを判定する判定ステップと、
前記判定ステップで、前記複数の融解曲線データが存在すると判定された場合に、前記複数の融解曲線データの各々が示す融解曲線の各々を、１つの座標系内に表示する表示ステップと、
を含む処理をコンピュータに実行させる遺伝子変異分析方法。
【請求項２】
前記表示ステップで、前記融解曲線の各々を１つの座標系内に重ねて表示するか、又は、所定時間毎に切り替えながら表示する請求項１記載の遺伝子変異分析方法。
【請求項３】
前記表示ステップで、前記融解曲線の各々を前記１つの座標系内に表示する際に、前記融解曲線の各々における波形のピークを一致させて表示する請求項１又は請求項２記載の遺伝子変異分析方法。
【請求項４】
前記ピークを、前記対象遺伝子の野生型の融解温度を示すピークとする請求項３記載の遺伝子変異分析方法。
【請求項５】
前記融解曲線の相互の類似度を計算する計算ステップをさらに含む処理を前記コンピュータに実行させる請求項１〜請求項４のいずれか１項記載の遺伝子変異分析方法。
【請求項６】
前記計算ステップで、前記対象遺伝子の野生型及び変異型の少なくとも一方の融解温度を含む所定温度範囲に含まれる融解曲線の部分の類似度を計算する請求項５記載の遺伝子変異分析方法。
【請求項７】
前記計算ステップで、前記類似度として、
融解曲線における前記対象遺伝子の野生型の融解温度を含む第１所定温度範囲の上限温度及び下限温度の各々に対応する融解曲線上の２点を結ぶ線と融解曲線とで囲まれた第１領域の重心の位置、
融解曲線における前記対象遺伝子の変異型の融解温度を含む第２所定温度範囲の上限温度及び下限温度の各々に対応する融解曲線上の２点を結ぶ線と融解曲線とで囲まれた第２領域の重心の位置、
前記第１領域の面積、
前記第２領域の面積、
前記第１領域と前記第２領域との面積比、
前記第１所定温度範囲に含まれる融解曲線の長さ、及び
前記第２所定温度範囲に含まれる融解曲線の長さの少なくとも１つの分析データの類似度を計算する
請求項５又は請求項６記載の遺伝子変異分析方法。
【請求項８】
前記計算ステップで計算された前記対象遺伝子の野生型の融解温度を含む所定温度範囲に含まれる融解曲線の部分の類似度が予め定めた閾値に満たない場合に、表示されている融解曲線を示す融解曲線データの各々の測定時の条件の相違が、融解曲線の相互の類似度に影響していることを示すメッセージを報知する報知ステップをさらに含む処理をコンピュータに実行させる請求項６又は請求項７記載の遺伝子変異分析方法。
【請求項９】
前記表示ステップで表示された融解曲線の各々に対して、所定温度範囲が新たに設定された場合に、新たに設定された所定温度範囲に基づいて、前記類似度を再計算する再計算ステップをさらに含む処理をコンピュータに実行させる請求項６〜請求項８のいずれか１項記載の遺伝子変異分析方法。
【請求項１０】
測定結果が既知のコントロール液に対する融解温度測定により得られた融解曲線データ、及び前記対象遺伝子に関する融解温度測定の測定時に使用された前記対象遺伝子の量の少なくとも１つの情報が取得可能な場合、前記表示ステップで、前記少なくとも１つの情報に基づいて、融解曲線の各々をキャリブレーションした上で表示する請求項１〜請求項９のいずれか１項記載の遺伝子変異分析方法。
【請求項１１】
前記対象遺伝子を、予め定めた種類の遺伝子とする請求項１〜請求項１０のいずれか１項記載の遺伝子変異分析方法。
【請求項１２】
前記予め定めた種類の遺伝子を、該遺伝子の野生型に対する変異型の比が経時的に変化する遺伝子とする請求項１１記載の遺伝子変異分析方法。
【請求項１３】
対象被験者の検体試料について、対象遺伝子に関する融解温度測定により得られ、かつ前記検体試料の採取時が各々異なる複数の融解曲線データが存在するか否かを判定する判定部と、
前記判定部により、前記複数の融解曲線データが存在すると判定された場合に、前記複数の融解曲線データの各々が示す融解曲線の各々を、１つの座標系内に表示する表示制御部と、
を含む遺伝子変異分析装置。
【請求項１４】
請求項１〜請求項１２のいずれか１項記載の遺伝子変異分析方法に含まれる各ステップをコンピュータに実行させるための遺伝子変異分析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、遺伝子変異分析方法、遺伝子変異分析装置、及び遺伝子変異分析プログラムに関する。
【背景技術】
【０００２】
従来、融解温度（Melting Temperature、Ｔｍ）測定により得られた融解曲線のグラフにおいて、シグナル値のピークが現れている温度を含む所定の温度範囲の融解曲線により囲まれた領域の面積に基づいて、遺伝子の野生型と変異型との存在比を判断する技術が存在する（例えば、特許文献１）。また、このようなＴｍ測定の結果を画面表示する技術が存在する。
【０００３】
また、疾患の進行により遺伝子変異量が増加する遺伝子が知られている（例えば、特許文献２）。このような遺伝子に対して、上記の従来技術により、遺伝子の変異の割合を求め、現在の変異の状態を判断することができる。
【０００４】
さらに、遺伝子変異以外の臨床技術の分野において、例えば、血糖値変化、血圧変化などの日々の変化をグラフ化して比較する技術が汎用的に活用されている（例えば、特許文献３及び４）。また、遺伝子系の分野でも、痴呆症に関する遺伝子のペプチドの濃度変化を検出する技術が提案されている（例えば、特許文献５）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１２−２９６８９号公報
【特許文献２】国際公開第２０１１／０５２７５４号パンフレット
【特許文献３】特公平７−３７９４６号公報
【特許文献４】特許第５３６２３０９号明細書
【特許文献５】特表２００４−５３１２５０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
例えば、特許文献１に記載の技術では、１回のＴｍ測定毎に得られた融解曲線を分析して、遺伝子の変異の割合などの結果を得ることが想定されており、複数回のＴｍ測定の結果を比較することは考慮されていない。この背景として、例えば、肥満遺伝子やアルコール関連遺伝子に代表されるように、野生型か変異型かを判断する測定が生涯で一度きりで満たされる遺伝子が多いことが挙げられる。また、遺伝子の測定は、専門部署や外注検査会社などで行われることが主流であり、臨床現場では、専門部署や外注検査会社から、変異あり、又は変異なしといった報告を受ける、というケースが多い。そのため、臨床現場において、Ｔｍ測定の結果を分析するという機会自体が少なく、複数回のＴｍ測定の結果を比較するというニーズが少なかったことも挙げられる。
【０００７】
しかしながら、癌遺伝子に代表されるような疾患の進行により遺伝子変異量が増加する遺伝子については、臨床現場におけるＴｍ測定の結果分析に対する重要度が増してきている。また、リアルタイムなオーダーメイド治療を行うためには、遺伝子の変異の進行状況を容易に把握することが望まれる。
【０００８】
本発明は、遺伝子の変異の進行状況を容易に把握することができる遺伝子変異分析方法、装置、及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
一つの態様では、遺伝子変異分析方法は、対象被験者の検体試料について、対象遺伝子に関する融解温度測定により得られ、かつ前記検体試料の採取時が各々異なる複数の融解曲線データが存在するか否かを判定する判定ステップと、前記判定ステップで、前記複数の融解曲線データが存在すると判定された場合に、前記複数の融解曲線データの各々が示す融解曲線の各々を、１つの座標系内に表示する表示ステップと、を含む処理をコンピュータに実行させる方法である。これにより、遺伝子の変異の進行状況を容易に把握することができる。
【００１０】
また、前記表示ステップで、前記融解曲線の各々を前記１つの座標系内に重ねて表示するか、又は、所定時間毎に切り替えながら表示することができる。
【００１１】
また、前記表示ステップで、前記融解曲線の各々を１つの座標系内に表示する際に、前記融解曲線の各々における波形のピークを一致させて表示することができる。また、前記ピークを、前記対象遺伝子の野生型の融解温度を示すピークとすることができる。これにより、１つの座標系内に表示された複数の融解曲線の比較がより容易になる。
【００１２】
また、前記融解曲線の相互の類似度を計算する計算ステップをさらに含む処理を前記コンピュータに実行させることができる。計算された類似度は、１つの座標系内に表示された複数の融解曲線を比較する際の補助情報として利用することができる。
【００１３】
また、前記計算ステップで、前記対象遺伝子の野生型及び変異型の少なくとも一方の融解温度を含む所定温度範囲に含まれる融解曲線の部分の類似度を計算することができる。融解曲線の中でも、融解温度を含む所定温度範囲の部分の形状に、より特徴が出易いため、所定の温度範囲の部分の類似度を計算することで、より有用な類似度を計算することができる。
【００１４】
また、前記計算ステップで、前記類似度として、融解曲線における前記対象遺伝子の野生型の融解温度を含む第１所定温度範囲の上限温度及び下限温度の各々に対応する融解曲線上の２点を結ぶ線と融解曲線とで囲まれた第１領域の重心の位置、融解曲線における前記対象遺伝子の変異型の融解温度を含む第２所定温度範囲の上限温度及び下限温度の各々に対応する融解曲線上の２点を結ぶ線と融解曲線とで囲まれた第２領域の重心の位置、前記第１領域の面積、前記第２領域の面積、前記第１領域と前記第２領域との面積比、前記第１所定温度範囲に含まれる融解曲線の長さ、及び前記第２所定温度範囲に含まれる融解曲線の長さの少なくとも１つの分析データの類似度を計算することができる。これらの分析データは、融解曲線の形状の特徴を表現するのに適したパラメータであり、かつ簡易な計算で求めることができるため、類似度の計算に用いることに適している。
【００１５】
また、前記計算ステップで計算された前記対象遺伝子の野生型の融解温度を含む所定温度範囲に含まれる融解曲線の部分の類似度が予め定めた閾値に満たない場合に、表示されている融解曲線を示す融解曲線データの各々の測定時の条件の相違が、融解曲線の相互の類似度に影響していることを示すメッセージを報知する報知ステップをさらに含む処理をコンピュータに実行させることができる。これにより、融解曲線の類似度が適切に計算されていない可能性があることを把握することができる。
【００１６】
また、前記表示ステップで表示された融解曲線の各々に対して、所定温度範囲が新たに設定された場合に、新たに設定された所定温度範囲に基づいて、前記類似度を再計算する再計算ステップをさらに含む処理をコンピュータに実行させることができる。これにより、測定時の条件の相違を緩和して、融解曲線の類似度の精度を向上させることができる。
【００１７】
また、測定結果が既知のコントロール液に対する融解温度測定により得られた融解曲線データ、及び前記対象遺伝子に関する融解温度測定の測定時に使用された前記対象遺伝子の量の少なくとも１つの情報が取得可能な場合、前記表示ステップで、前記少なくとも１つの情報に基づいて、融解曲線の各々をキャリブレーションした上で表示することができる。これにより、融解曲線の形状の比較をより容易にすることができる。
【００１８】
また、前記対象遺伝子を、予め定めた種類の遺伝子とすることができる。また、前記予め定めた種類の遺伝子を、該遺伝子の野生型に対する変異型の比が経時的に変化する遺伝子とすることができる。
【００１９】
また、一つの態様では、遺伝子変異分析装置は、対象被験者の検体試料について、対象遺伝子に関する融解温度測定により得られ、かつ前記検体試料の採取時が各々異なる複数の融解曲線データが存在するか否かを判定する判定部と、前記判定部により、前記複数の融解曲線データが存在すると判定された場合に、前記複数の融解曲線データの各々が示す融解曲線の各々を、１つの座標系内に表示する表示制御部と、を含んで構成することができる。
【００２０】
また、一つの態様では、遺伝子変異分析プログラムは、上記の遺伝子変異分析方法に含まれる各ステップをコンピュータに実行させるためのプログラムである。
【発明の効果】
【００２１】
以上説明したように、本発明の遺伝子変異分析方法、装置、及びプログラムによれば、遺伝子の変異の進行状況を容易に把握することができる、という効果が得られる。
【図面の簡単な説明】
【００２２】
【図１】遺伝子変異分析装置の概略構成を示すブロック図である。
【図２】第１実施形態に係る遺伝子変異分析装置の機能ブロック図である。
【図３】第１実施形態における融解曲線データベース（ＤＢ）の一例を示す図である。
【図４】融解曲線の温度範囲を説明するための図である。
【図５】第１実施形態における遺伝子変異分析処理の一例を示すフローチャートである。
【図６】選択画面の一例を示す図である。
【図７】第１実施形態におけるグラフ表示画面の一例を示す図である。
【図８】融解曲線のグラフの重ね合わせについて説明するための図である。
【図９】融解曲線のグラフの重ね合わせについて説明するための図である。
【図１０】第２実施形態に係る遺伝子変異分析装置の機能ブロック図である。
【図１１】第２実施形態における融解曲線ＤＢの一例を示す図である。
【図１２】第２実施形態における遺伝子変異分析処理の一例を示すフローチャートである。
【図１３】第２実施形態におけるグラフ表示画面の一例を示す図である。
【発明を実施するための形態】
【００２３】
以下、図面を参照して、本発明に係る遺伝子変異分析装置の実施形態を詳細に説明する。以下の各実施形態では、融解温度（Melting Temperature、Ｔｍ）測定により得られた融解曲線データに基づいて、遺伝子の変異を分析する遺伝子変異分析装置について説明する。
【００２４】
＜第１実施形態＞
図１に示すように、第１実施形態に係る遺伝子変異分析装置１０は、入力部１２と、表示部１４と、コンピュータ１６とを備える。入力部１２は、キーボード、マウス、タッチパネル、バーコードリーダ等の読取装置等で構成される。また、入力部１２は、ユーザの操作に応じて各種情報をコンピュータ１６へ入力する。表示部１４は、各種情報を表示するためのディスプレイである。
【００２５】
コンピュータ１６は、遺伝子変異分析装置１０全体の制御を司るＣＰＵ２０、後述する遺伝子変異分析処理等の各種プログラムを記憶した記憶媒体としてのＲＯＭ２２、ワークエリアとしてデータを一時的に格納するＲＡＭ２４、各種情報が記憶された記憶手段としてのメモリ２６、入出力ポート（Ｉ／Ｏポート）２８、ネットワークインターフェース（ネットワークＩ／Ｆ）３０、及びこれらを接続するバスを含んで構成されている。Ｉ／Ｏポート２８には、入力部１２及び表示部１４が接続されている。なお、コンピュータ１６は、ＨＤＤ等をさらに含んで構成されてもよい。
【００２６】
図２に、遺伝子変異分析装置１０に含まれるコンピュータ１６により実現される各機能を機能ブロックで表した機能ブロック図を示す。図２に示すように、コンピュータ１６は、機能的には、受付部４１と、判定部４２と、表示制御部４３とを含んだ構成で表すことができる。また、コンピュータ１６の所定の記憶領域（例えば、メモリ２６）には、融解曲線データベース（ＤＢ）５０が記憶される。
【００２７】
受付部４１は、入力部１２を介して入力された各種情報を受け付ける。入力部１２から、どの被験者のどの遺伝子についての融解曲線のグラフを表示するかを指定するための対象被験者及び対象遺伝子の情報が入力された場合には、受付部４１は、受け付けた情報を判定部４２へ受け渡す。また、入力部１２から、Ｔｍ測定に関するデータ（詳細は後述）が入力された場合には、受付部４１は、受け付けたＴｍ測定に関するデータを融解曲線ＤＢ５０に記憶する。
【００２８】
図３に、融解曲線ＤＢ５０の一例を示す。図３の例では、１行（１レコード）が１回分のＴｍ測定に関するデータ（以下、「Ｔｍ測定データ」という）に相当する。Ｔｍ測定データの各々は、被験者の識別情報である被験者ＩＤ、被験者から採取された検査対象の検体試料の識別情報である検体ＩＤ、Ｔｍ測定の測定日、測定対象の遺伝子を示す項目名、測定された融解曲線データ、及び融解曲線を分析した各種分析データの項目が含まれる。なお、本実施形態では、測定日が異なるＴｍ測定データは、それぞれ検体試料の採取日が異なるものとする。すなわち、被験者ＩＤが同一で、かつ測定日が異なるＴｍ測定データは、その被験者ＩＤが示す被験者から、異なる採取日に採取された検体試料についてのＴｍ測定データであることを表す。
【００２９】
融解曲線データは、Ｔｍ測定において検出された各温度における吸光度又は蛍光強度等の検出信号の微分値である。横軸を温度、縦軸を検出信号の微分値とする座標系のグラフに融解曲線データの各値をプロットすることにより、融解曲線のグラフが得られる。
【００３０】
分析データは、融解曲線から得られる遺伝子の野生型及び変異型の各々に関するデータである。図４に示すように、融解曲線に現れたピークを対象遺伝子の野生型の融解温度ＴｍＷ及び変異型の融解温度ＴｍＭとして検出する。そして、野生型の融解温度ＴｍＷを含む所定の温度範囲を野生型温度範囲、変異型の融解温度ＴｍＭを含む所定の温度範囲を変異型温度範囲とする。なお、所定の温度範囲は、例えば、融解温度ＴｍＷ又はＴｍＭの前後数度の範囲（例えば、±３℃）とすることができる。図３に示す融解曲線ＤＢ５０では、野生型温度範囲及び変異型温度範囲の各々について、下限温度、融解温度、及び上限温度の組を、「野生型温度範囲」及び「変異型温度範囲」の項目に記憶している。
【００３１】
また、図３に示す融解曲線ＤＢ５０では、分析データとして、「野生領域面積」、「変異領域面積」、及び「変異野生面積比」の項目を含む。「野生領域面積」は、図４の斜線部Ａに示すように、野生型温度範囲の上限温度及び下限温度の各々に対応する融解曲線上の２点を結ぶ線と融解曲線とで囲まれた領域（以下、「野生領域」）の面積である。同様に、「変異領域面積」は、図４の斜線部Ｂに示すように、変異型温度範囲の上限温度及び下限温度の各々に対応する融解曲線上の２点を結ぶ線と融解曲線とで囲まれた領域（以下、「変異領域」）の面積である。「変異野生面積比」は、野生領域面積に対する変異領域面積の割合である。
【００３２】
なお、分析データは上記の例に限定されず、野生領域の重心、変異領域の重心、野生型温度範囲に含まれる融解曲線の長さ、及び変異型温度範囲に含まれる融解曲線の長さなど、他の分析データを含んでもよい。
【００３３】
判定部４２は、受付部４１から受け渡された対象被験者及び対象遺伝子の情報に基づいて、融解曲線ＤＢ５０から、融解曲線のグラフを表示する対象のＴｍ測定データを検索し、検索したＴｍ測定データを表示制御部４３へ受け渡す。また、判定部４２は、検索されたＴｍ測定データから、複数のＴｍ測定データがユーザにより選択されたか否か、及び対象遺伝子が、遺伝子の野生型に対する変異型の比が経時的に変化（以下、「変異が進行」ともいう）する遺伝子として予め定めた種類の遺伝子か否かを判定し、判定結果を表示制御部４３へ受け渡す。
【００３４】
表示制御部４３は、判定部４２で検索されたＴｍ測定データの一覧を作成し、それぞれのＴｍ測定データを選択可能な状態で表示部１４に表示するように制御する。また、表示制御部４３は、複数のＴｍ測定データが選択され、かつ対象遺伝子が予め定めた種類の遺伝子である旨の判定結果を判定部４２から受け渡された場合には、選択された複数のＴｍ測定データの各々に含まれる融解曲線データが示す融解曲線の各々を、１つの座標系内に重ね合わせて表示するように、表示部１４を制御する。また、表示制御部４３は、１つのＴｍ測定データのみが選択された旨の判定結果を判定部４２から受け渡された場合には、選択されたＴｍ測定データに含まれる融解曲線データが示す融解曲線が表示されるように表示部１４を制御する。また、表示制御部４３は、複数のＴｍ測定データが選択されているが、対象遺伝子が予め定めた種類の遺伝子ではない旨の判定結果を判定部４２から受け渡された場合には、選択された複数のＴｍ測定データの各々に含まれる融解曲線データが示す融解曲線の各々を並べて表示するように、表示部１４を制御する。
【００３５】
次に、第１実施形態に係る遺伝子変異分析装置１０の作用について説明する。Ｔｍ測定データが、入力部１２を介して遺伝子変異分析装置１０に入力されると、受付部４１が、入力されたＴｍ測定データを融解曲線ＤＢ５０に記憶する。なお、通常、病気の進行や投薬による効果を経過観察するために、同一の被験者に対する同一の測定項目（対象遺伝子）についてのＴｍ測定は、半年毎、１年毎のような長期の間隔を経て行われる。本実施形態では、融解曲線ＤＢ５０に記憶されたＴｍ測定データには、同一の被験者に対する同一の測定項目について、上記のように長期の間隔を経て異なる採取日に採取され、かつ異なる測定日の各々に測定されたＴｍ測定データが含まれるものとする。融解曲線ＤＢ５０に、複数のＴｍ測定データが記憶された状態で、融解曲線のグラフを表示する対象となる対象被験者及び対象遺伝子の情報が、入力部１２を介して遺伝子変異分析装置１０に入力されると、遺伝子変異分析装置１０において、図５に示す遺伝子変異分析処理が実行される。なお、遺伝子変異分析装置１０で実行される遺伝子変異分析処理は、本発明の遺伝子変異分析方法の一例である。
【００３６】
図５に示す遺伝子変異分析処理のステップＳ１１で、受付部４１が、入力された対象被験者及び対象遺伝子の情報を受け付け、判定部４２へ受け渡す。
【００３７】
次に、ステップＳ１２で、判定部４２が、受付部４１から受け渡された対象被験者及び対象遺伝子の情報に基づいて、融解曲線ＤＢ５０から、指定された対象被験者及び対象遺伝子についてのＴｍ測定データを検索する。具体的には、判定部４２は、融解曲線ＤＢ５０から、指定された対象被験者を示す被験者ＩＤを有し、かつ項目名が指定された対象遺伝子であるＴｍ測定データを検索する。判定部４２は、検索したＴｍ測定データを表示制御部４３へ受け渡す。そして、表示制御部４３が、例えば、図６に示すような選択画面６０を表示部１４に表示する。
【００３８】
図６に示す選択画面６０には、指定された対象被験者及び対象遺伝子の情報が表示される指定情報領域６２、検索されたＴｍ測定データの一覧が表示される一覧領域６４、一覧に含まれる各Ｔｍ測定データをユーザが選択するための選択ボックス６６、及びグラフの表示を指示する際にユーザに選択されるグラフ表示ボタン６８が含まれる。指定情報領域６２には、被験者に関する各種情報が登録された被験者情報ＤＢ（図示省略）等に記憶された情報を併せて表示することができる。図６の例では、被験者ＩＤに対応する被験者名を表示している。さらに、被験者の性別、生年月日等を表示してもよい。
【００３９】
一覧領域６４には、判定部４２で検索されたＴｍ測定データの各々に含まれる各項目の情報を表示する。なお、図６の例では、Ｔｍ測定データに含まれる測定日、野生領域面積、変異領域面積、及び変異野生面積比の各項目の情報と、融解曲線データが示す融解曲線のグラフのサムネイル画像とを表示している。また、各Ｔｍ測定データに対応して、選択ボックス６６が配置される。ユーザは、表示したい融解曲線を示す融解曲線データを含むＴｍ測定データに対応する選択ボックス６６を選択する。
【００４０】
次に、ステップＳ１３で、判定部４２が、選択画面６０において、グラフ表示ボタン６８が選択されたか否かを判定する。選択されていない場合には、本ステップの判定を繰り返し、選択された場合には、処理はステップＳ１４へ移行する。
【００４１】
ステップＳ１４では、判定部４２が、選択画面６０における選択ボックス６６の各々の選択状態に基づいて、複数のＴｍ測定データが選択されているか否かを判定する。複数のＴｍ測定データが選択されている場合には、処理はステップＳ１５へ移行し、１つのＴｍ測定データのみが選択されている場合には、処理はステップＳ１７へ移行する。なお、判定部４２で検索されたＴｍ測定データが１つしかない場合、すなわち、選択画面６０の一覧領域６４に１つのＴｍ測定データのみが表示されている場合には、そのＴｍ測定データに対応する選択ボックス６６を自動的に選択すると共に、本ステップで否定判定されるようにしてもよい。
【００４２】
ステップＳ１５では、判定部４２が、上記ステップＳ１１で指定された対象遺伝子が、変異が進行する遺伝子として予め定めた種類の遺伝子か否かを判定する。予め定めた種類の遺伝子の場合には、処理はステップＳ１６へ移行し、予め定めた種類の遺伝子ではない場合には、処理はステップＳ１８へ移行する。
【００４３】
ステップＳ１６では、判定部４２が、複数のＴｍ測定データが選択され、かつ対象遺伝子が予め定められた種類の遺伝子である旨の判定結果を、表示制御部４３へ受け渡す。そして、表示制御部４３が、例えば、図７に示すようなグラフ表示画面７０を表示部１４に表示する。図７に示すグラフ表示画面７０には、指定情報領域６２、複数の融解曲線が重ね合わせて表示される第１グラフ領域７２、及び個々の融解曲線のグラフが表示される第２グラフ領域７４が含まれる。
【００４４】
表示制御部４３は、第１グラフ領域７２に、選択された複数のＴｍ測定データの各々に含まれる融解曲線データが示す融解曲線の各々を、１つの座標系内に重ね合わせて表示する。また、表示制御部４３は、第２グラフ領域７４に、選択された複数のＴｍ測定データの各々に含まれる融解曲線データが示す融解曲線のグラフの各々を、それぞれを識別するための情報（図７の例では、測定日）及び各融解曲線の分析データと共に、並べて表示する。
【００４５】
ここで、融解曲線のグラフを重ね合わせて表示することの理由について説明する。
【００４６】
例えば、図８に示すように、ある被験者のある測定項目について、２０１２／１２／１４（前回）にＴｍ測定された融解曲線９０Ａが得られているとする。融解曲線９０Ａからは、微小な変異がみられる。そして、同一の被験者、かつ同一の測定項目について、１年半後の２０１４／０８／２５（今回）に行われたＴｍ測定により、融解曲線９０Ｂが得られたとする。
【００４７】
今回の融解曲線９０Ｂからも、前回の融解曲線９０Ａと同様に、変異の存在は確認することができる。しかし、融解曲線９０Ｂを単体で見た場合には、前回測定時からの変異の進行状況を把握することができない。変異の進行状況を把握するためには、過去の測定結果と比較する必要がある。そこで、前回の融解曲線９０Ａと今回の融解曲線９０Ｂとを比較することを考える。ここで、図８の上段に示すように、前回の融解曲線９０Ａと今回の融解曲線９０Ｂとを並べて表示して比較した場合には、波形が変化した部分を把握し難く、適切に変異の進行状況を読み取ることが困難な場合がある。そこで、各融解曲線９０Ａ及び９０Ｂについての分析データを用いて、数値上で比較してみると、変異領域の面積が１３．３から３５．０に増加しており、一見、変異が進行しているとも考えられる。
【００４８】
しかし、図８の下段に示すように、融解曲線９０Ａと融解曲線９０Ｂとを重ね合わせてみると、変異型温度範囲のピーク付近の波形は、融解曲線９０Ａと融解曲線９０Ｂとで類似している。従って、変異領域の面積の差は、変異領域の面積を求める際に設定した面積の下底のずれが原因であると推定できる。すなわち、変異の進行により変異領域の面積が増加したのではなく、面積を求める際の温度範囲の設定の相違により、面積に差が出たものであり、変異の進行は見られないと判断することができる。
【００４９】
また、別の例として、図９に示すように、ある被験者のある測定項目について、２０１３／０９／１２（前回）にＴｍ測定された融解曲線９０Ｃが得られているとする。融解曲線９０Ｃからは、微小な変異がみられる。そして、同一の被験者、かつ同一の測定項目について、１年後の２０１４／１０／２５（今回）に行われたＴｍ測定により、融解曲線９０Ｄが得られたとする。各融解曲線９０Ｃ及び９０Ｄについての分析データについて比較してみると、変異領域の面積が６７．０から７９．０にわずかに増加している。そこで、図９の下段に示すように、融解曲線９０Ｃと融解曲線９０Ｄとを重ね合わせてみると、融解曲線９０Ｄの変異型温度範囲のピーク付近の波形が、融解曲線９０Ｃに比べ膨らんでいることが確認できる。従って、変異が進行していると判断することができる。
【００５０】
図８の例で説明したように、領域の面積を求める際の温度範囲の設定に、融解曲線毎に差が出る要因として、各測定時の温度計などの装置構成、検体試料の差、試薬のロットの違いなどの測定系の相違により、融解曲線から検出されるピーク位置にずれが生じることが考えられる。また、Ｔｍ測定で温度を上昇させる際の温度の刻み方（例えば、１度ずつ）によっても、融解曲線から検出されるピーク位置にずれが生じる場合もある。ピーク位置にずれが生じることにより、そのピーク位置に基づいて設定される温度範囲も、測定毎にずれが生じる。従って、変異領域の面積などの分析データを比較しただけでは、適切な変異の進行状況を把握することができないため、融解曲線を重ね合わせて表示し、融解曲線の形状の変化を確認することで、変異の進行状況を把握するものである。
【００５１】
図５に示す遺伝子変異分析処理の説明に戻る。以下では、２つのＴｍ測定データが選択されているものとして説明する。
【００５２】
ステップＳ１６で、表示制御部４３は、具体的には、選択された２つのＴｍ測定データに含まれる融解曲線データが示す融解曲線の各々を、１つの座標系内のグラフに展開する。この際、表示制御部４３は、融解曲線の各々を、異なる色で表示する。また、一方の融解曲線における野生領域及び変異領域と、他方の融解曲線における野生領域及び変異領域とを異なる色で表示し、２つの融解曲線の野生領域同士、及び変異領域同士が重なる部分は、さらに異なる色（例えば、一方の色と他方の色とを重ね合わせた色）で表示する。これにより、２つの融解曲線の形状の差をより把握し易くする。
【００５３】
なお、上述したように、測定系の相違に基づく融解曲線毎のピーク位置のずれを考慮して、２つの融解曲線を重ね合わせる際に、横軸のスケールを完全に一致させるのではなく、２つの融解曲線のピークの温度を一致させるように重ね合わせてもよい。なお、変異が複合的に表れる場合などでは、変異型のピークの温度は測定毎に変わる可能性もあるため、野生型のピークの温度を一致させるようにするとよい。これにより、２つの融解曲線の形状をより比較し易くなる。また、測定系の相違を考慮して、２つの融解曲線の最大値を一致させるように、縦軸のスケールを調整してもよい。これにより、より形状に着目して、２つの融解曲線の比較を行うことができる。
【００５４】
一方、上記ステップＳ１４で、１つのＴｍ測定データのみが選択されていると判定されて、ステップＳ１７へ移行した場合には、判定部４２が、１つのＴｍ測定データのみが選択された旨の判定結果を、表示制御部４３へ受け渡す。そして、表示制御部４３が、選択されたＴｍ測定データに含まれる融解曲線データが示す融解曲線のグラフを、表示部１４に表示する。この際の表示画面は、例えば図７に示すグラフ表示画面７０の第１グラフ領域７２を設けず、第１グラフ領域７２の範囲まで第２グラフ領域７４を拡大させたグラフ表示画面とし、この第２グラフ領域７４に、選択された１つの融解曲線のグラフを表示すればよい。
【００５５】
また、上記ステップＳ１５で、対象遺伝子が予め定めた種類の遺伝子ではないと判定されて、ステップＳ１８へ移行した場合には、判定部４２が、複数の融解曲線データが選択されているが、対象遺伝子が予め定めた種類の遺伝子ではない旨の判定結果を、表示制御部４３へ受け渡す。この場合、遺伝子の変異の状況を把握することは重要ではないため、上記ステップＳ１６のような重ね合わせ表示は行わず、表示制御部４３は、選択された複数の融解曲線データの各々が示す融解曲線のグラフ各々を並べて表示するように、表示部１４を制御する。この際の表示画面は、例えば図７に示すようなグラフ表示画面７０の第１グラフ領域７２を設けず、第１グラフ領域７２の範囲まで第２グラフ領域７４を拡大させたグラフ表示画面とし、この第２グラフ領域７４に、選択された複数の融解曲線のグラフを並べて表示すればよい。
【００５６】
以上説明したように、第１実施形態に係る遺伝子変異分析装置によれば、同一の被験者かつ同一の対象遺伝子について、検体試料の採取時が異なる複数の融解曲線が存在する場合には、それら複数の融解曲線を１つの座標系内に重ね合わせて表示する。これにより、遺伝子の変異の進行状況を容易に把握することができる。また、融解曲線を重ね合わせて波形の形状を比較することにより、変異の進行状況等の情報を読み取ることは、融解曲線単体を見て変異の状態等の情報を読み取る場合に比べて容易であるため、融解曲線の波形を用いた判断に慣れていない人（例えば、臨床現場の医師や検査技師など）にとっても、遺伝子の変異の進行状況の把握が容易になる。これにより、臨床現場におけるリアルタイムなオーダーメイド治療にも寄与し得る。
【００５７】
＜第２実施形態＞
次に、第２実施形態について説明する。なお、第２実施形態に係る遺伝子変異分析装置について、第１実施形態に係る遺伝子変異分析装置１０と同様の構成については、同一符号を付して、詳細な説明を省略する。
【００５８】
図１に示すように、第２実施形態に係る遺伝子変異分析装置２１０は、入力部１２と、表示部１４と、コンピュータ２１６とを備える。図１０に、遺伝子変異分析装置２１０に含まれるコンピュータ２１６により実現される各機能を機能ブロックで表した機能ブロック図を示す。図１０に示すように、コンピュータ２１６は、機能的には、受付部４１と、判定部４２と、表示制御部２４３と、計算部４４と、報知部４５と、再計算部４６とを含んだ構成で表すことができる。また、コンピュータ２１６の所定の記憶領域（例えば、メモリ２６）には、融解曲線ＤＢ２５０が記憶される。
【００５９】
図１１に、融解曲線ＤＢ２５０の一例を示す。図１１の例では、第１実施形態における融解曲線ＤＢ５０（図３）の項目に加え、各Ｔｍ測定時に使用された遺伝子の量を示す「遺伝子量」の項目が含まれる。
【００６０】
表示制御部２４３は、第１実施形態における表示制御部４３と同様、Ｔｍ測定データが複数選択され、かつ対象遺伝子が予め定めた種類の遺伝子である場合に、選択されたＴｍ測定データの各々に含まれる融解曲線データが示す融解曲線の各々を重ね合わせて表示する。この際、表示制御部２４３は、選択されたＴｍ測定データの各々に含まれる「遺伝子量」の情報に基づいて、融解曲線データの各々をキャリブレーションする。例えば、一方の融解曲線データの各値（温度毎の検出値の微分値）に、一方の遺伝子量に対する他方の遺伝子量の割合を乗算することにより、キャリブレーションすることができる。これにより、測定時の遺伝子量の差を緩和して、融解曲線の形状の比較をより容易にすることができる。
【００６１】
また、表示制御部２４３は、後述する計算部４４又は再計算部４６で計算される融解曲線の類似度も、重ね合わせた融解曲線のグラフと共に表示する。
【００６２】
計算部４４は、重ね合わせて表示する融解曲線の相互の類似度を計算する。重ね合わせて表示する融解曲線の各々を示す融解曲線データの各々は、異なる採取時に採取され、異なる測定時に測定されたものであるので、計算部４４により計算される類似度の大小が、融解曲線の経時的な変化を表す指標となる。
【００６３】
具体的には、計算部４４は、Ｔｍ測定データに含まれる分析データの同一項目同士の差又は比を、対応する融解曲線同士の類似度として計算することができる。この場合、分析データの同一項目同士の差又は比が大きい程、融解曲線同士の類似度は低いことを表す。例えば、測定日が古い方の融解曲線における野生領域の面積に対する、測定日が新しい方の融解曲線における野生領域の面積の増加量で表される野生領域面積変化量を類似度とすることができる。また、測定日が古い方の融解曲線における変異領域の面積に対する、測定日が新しい方の融解曲線における変異領域の面積の増加量で表される変異領域面積変化量を類似度とすることができる。また、測定日が古い方の融解曲線における変異野生面積比に対する、測定日が新しい方の融解曲線における変異野生面積比の増加量で表される変異野生面積比変化量を類似度とすることができる。
【００６４】
さらに、計算部４４は、Ｔｍ測定データに含まれる分析データ以外の分析データを融解曲線から計算し、計算した分析データを用いて、融解曲線の類似度を計算してもよい。例えば、計算部４４は、融解曲線データに基づいて、野生領域の重心、変異領域の重心、野生型温度範囲に含まれる融解曲線の長さ、変異型温度範囲に含まれる融解曲線の長さ等の分析データを計算する。重心や長さの計算は従来既知の手法を用いることができる。計算部４４は、Ｔｍ測定データに含まれる分析データを用いた類似度の計算と同様に、計算した分析データを用いて類似度を計算することができる。特に、野生領域の重心、及び変異領域の重心は、融解曲線の形状の変化に敏感に影響するパラメータであるため、類似度の計算に用いるパラメータとして有効である。
【００６５】
このように計算された融解曲線の類似度を、重ね合わせた融解曲線のグラフと共に表示することにより、重ね合わせた融解曲線から変異の進行状況を読み取る際の補助的な情報として、融解曲線の類似度を参照することができる。
【００６６】
報知部４５は、野生型に関する分析データを用いて計算部４４で計算された類似度が低い場合に、表示されている融解曲線を示す融解曲線データの各々の測定時の条件の相違が、融解曲線の相互の類似度に影響していることを示すメッセージを報知する。変異が進行する遺伝子であっても、野生型温度範囲に含まれる融解曲線の部分の形状は、経時的な変化が小さいため、野生型に関する分析データを用いて計算された類似度（以下、「野生型の類似度」という）が極端に低くなることは想定し難い。そこで、融解曲線の野生型の類似度が低い場合には、融解曲線データの各測定時における温度計などの装置構成、検体試料、試薬のロットなどの測定系に、融解曲線の類似度を適切に計算できないほどの相違があるものとみなし、その旨を報知するものである。
【００６７】
報知部４５は、具体的には、測定時の条件の相違により、適切な類似度が計算されていない可能性がある旨のメッセージを表示部１４に表示する。また、遺伝子変異分析装置２１０にスピーカーやＬＥＤなどの発光部を設けて、スピーカーからビープ音や音声メッセージを出力したり、ＬＥＤを点滅させたりすることにより、上記のメッセージを報知するようにしてもよい。
【００６８】
再計算部４６は、重ね合わせて表示されている融解曲線に対して、野生型温度範囲及び変異型温度範囲の少なくとも一方が、ユーザにより設定し直された場合に、新たに設定された温度範囲に基づいて、分析データを再計算して、融解曲線の類似度を再計算する。上述したように、融解曲線データの各測定時の条件には相違がある場合がある。このような場合には、各融解曲線に設定される温度範囲にも測定毎にずれが生じる場合がある。これが原因となって、重ね合わせた融解曲線の波形は類似して見えるにも関わらず、計算された類似度は低くなったり、逆に、重ね合わせた融解曲線の波形には変化が見られるにも関わらず、計算された類似度は高くなったりする場合がある。そこで、重ね合わせて表示された融解曲線に基づいて、温度範囲を再設定して類似度を再計算することで、測定時の条件の相違を緩和して、補助情報として表示する融解曲線の類似度の精度を高めることができる。
【００６９】
次に、第２実施形態に係る遺伝子変異分析装置２１０の作用について説明する。第２実施形態では、遺伝子変異分析装置２１０において、図１２に示す遺伝子変異分析処理が実行される。なお、第２実施形態における遺伝子変異分析処理において、第１実施形態における遺伝子変異分析処理（図５）と同様の処理となるステップについては、同一符号を付して、詳細な説明を省略する。
【００７０】
図１２に示す遺伝子変異分析処理のステップＳ１１〜Ｓ１５を経て、複数のＴｍ測定データが選択され、かつ対象遺伝子が予め定めた種類の遺伝子であると判定されると、処理はステップＳ２１へ移行する。
【００７１】
ステップＳ２１では、判定部４２が、複数のＴｍ測定データが選択され、かつ対象遺伝子が予め定められた種類の遺伝子である旨の判定結果を、表示制御部４３へ受け渡す。そして、表示制御部２４３が、選択されたＴｍ測定データの各々に含まれる「遺伝子量」の情報に基づいて、融解曲線データの各々をキャリブレーションする。次に、表示制御部２４３が、例えば、図１３に示すようなグラフ表示画面２７０を表示部１４に表示する。
【００７２】
図１３に示すグラフ表示画面２７０には、指定情報領域６２、複数の融解曲線が重ね合わせて表示される第１グラフ領域７２、及び個々の融解曲線が表示される第２グラフ領域７４が含まれる。さらに、グラフ表示画面２７０には、計算部４４又は再計算部４６で計算された類似度が表示される類似度領域７６、報知部４５によるメッセージが表示される報知領域７８、温度範囲を再設定するための設定ボックス８０、再計算を指示する際にユーザにより選択される再計算ボタン８２が含まれる。
【００７３】
表示制御部２４３は、第１実施形態における表示制御部４３と同様に、第１グラフ領域７２に、選択された融解曲線の各々を重ね合わせて表示し、第２グラフ領域７４に、個々の融解曲線のグラフの各々を並べて表示する。
【００７４】
次に、ステップＳ２２で、計算部４４が、グラフ表示画面２７０の第１グラフ領域７２に重ね合わせて表示された融解曲線同士の類似度を計算し、計算した類似度を表示制御部２４３へ受け渡す。表示制御部２４３は、受け渡された類似度を類似度領域７６に表示する。
【００７５】
次に、ステップＳ２３で、報知部４５が、融解曲線の野生型の類似度が低い否かを判定する。例えば、融解曲線の野生型の類似度として、野生領域面積変化量を用いることができる。この場合、野生領域面積変化量が大きい程、類似度が低いことを表すため、野生領域面積変化量が予め定めた閾値を超えているか否かを判定すればよい。融解曲線の野生型の類似度が低い場合には、処理はステップＳ２４へ移行する。
【００７６】
ステップＳ２４では、報知部４５が、表示されている融解曲線を示す融解曲線データの各々の測定時の条件の相違が、融解曲線の類似度に影響している旨のメッセージ、例えば「測定時の条件の相違により、類似度が適切に計算されていない可能性があります。」等のメッセージを、報知領域７８に表示する。一方、融解曲線の野生型の類似度が低くない場合には、ステップＳ２４の処理はスキップされる。
【００７７】
次に、ステップＳ２５で、再計算部４６が、グラフ表示画面２７０の再計算ボタン８２が選択されたか否かを判定する。ユーザは、重ね合わせて表示された融解曲線、融解曲線の類似度、報知領域７８に表示されたメッセージ等を参考にして、類似度を再計算したい場合には、設定ボックス８０に新たな温度範囲を設定し、再計算ボタン８２を選択する。なお、第１グラフ領域７２に、温度範囲設定のためのガイド線８４を配置し、ユーザがこのガイド線８４を融解曲線のグラフ上で所望の温度の位置に移動させることにより、新たな温度範囲を設定可能にしてもよい。再計算ボタン８２が選択された場合には、処理はステップＳ２６へ移行し、選択されていない場合には、遺伝子変異分析処理を終了する。
【００７８】
ステップＳ２６では、再計算部４６が、新たに設定された温度範囲に基づいて、野生領域の面積、変異領域の面積等の分析データを再計算し、処理はステップＳ２２に戻る。これにより、ステップＳ２２では、再計算部４６が、再計算した分析データに基づいて、融解曲線の類似度を再計算し、表示制御部２４３により、類似度領域７６に表示されている類似度が更新される。
【００７９】
以上説明したように、第２実施形態に係る遺伝子変異分析装置によれば、第１実施形態の効果に加え、融解曲線の野生領域の面積、変異領域の面積等の分析データを用いて計算した融解曲線の類似度を、重ね合わせた融解曲線のグラフから変異の進行状況を読み取る際の補助情報として用いることができる。
【００８０】
また、融解曲線データの各々の測定時の条件の相違を考慮して、温度範囲を再設定して融解曲線の類似度を再計算することで、測定時の条件の相違を緩和した精度の高い類似度を補助情報として用いることができる。
【００８１】
また、複数の融解曲線を重ね合わせて表示する際に、融解曲線データを測定した際に使用された遺伝子の量に基づいてキャリブレーションすることで、融解曲線の形状をより比較し易くすることができる。なお、キャリブレーションに利用する情報は、測定時の遺伝子の量に限定されない。例えば、測定結果が既知のコントロール液に対するＴｍ測定により得られた融解曲線データに基づいて、表示する融解曲線を示す融解曲線データをキャリブレーションしてもよい。
【００８２】
また、第１及び第２実施形態では、１つの座標系内に、検体試料の採取時が異なる複数の融解曲線を重ね合わせて表示する場合について説明したが、これに限定されない。例えば、１つの座標系内に、異なる融解曲線を所定時間毎に交互に表示してもよい。また、１つの座標系内に表示する融解曲線は２つに限定されず、３つ以上であってもよい。３つ以上の融解曲線を所定時間毎に表示する場合には、例えば、測定日が古い順に所定時間ずつ順次切り替えながら表示することができる。この場合、変異の進行状況をより直感的に把握可能なアニメーション表示とすることができる。また、３つ以上の融解曲線を表示する場合における融解曲線の類似度は、３つ以上の融解曲線から任意の２つの融解曲線の組を少なくとも１つ以上選択し、選択した組の融解曲線同士の類似度を計算すればよい。例えば、融解曲線データの測定日順で前後となる融解曲線同士の類似度をそれぞれ計算したり、測定日が最も古い融解曲線と測定日が最も新しい融解曲線との類似度を計算したりすることができる。
【００８３】
また、第１及び第２実施形態では、Ｔｍ測定データに含まれる融解曲線データを展開してグラフを表示する場合について説明したが、Ｔｍ測定データに、展開された融解曲線のグラフの画像データを予め含めておいてもよい。この場合、融解曲線を１つの座標系内に表示する際には、複数の画像データが示すグラフの各々の座標軸を一致させるような画像合成により、重ね合わせたり、交互に表示したりするようにすればよい。また、この場合も、第１実施形態で述べたように、融解曲線の波形のピークを一致させるようにしてもよい。
【００８４】
なお、上記の遺伝子変異分析処理ルーチンを規定したプログラムを記録媒体に記録して提供するようにしてもよい。
【符号の説明】
【００８５】
１０、２１０  遺伝子変異分析装置
１２   入力部
１４   表示部
１６、２１６  コンピュータ
４１   受付部
４２   判定部
４３、２４３  表示制御部
４４   計算部
４５   報知部
４６   再計算部
５０、２５０  融解曲線データベース
７０、２７０  グラフ表示画面
７２   第１グラフ領域
７４   第２グラフ領域
７６   類似度領域
７８   報知領域
８０   設定ボックス
８２   再計算ボタン
８４   ガイド線
２７０ グラフ表示画面
【図１】
2016208949.tif 000003
【図２】
2016208949.tif 000004
【図３】
2016208949.tif 000005
【図５】
2016208949.tif 000006
【図１０】
2016208949.tif 000007
【図１１】
2016208949.tif 000008
【図１２】
2016208949.tif 000009
【図４】
2016208949.tif 000010
【図６】
2016208949.tif 000011
【図７】
2016208949.tif 000012
【図８】
2016208949.tif 000013
【図９】
2016208949.tif 000014
【図１３】
2016208949.tif 000015