特開 2023-31319 分析装置、分析方法及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2023031319

(43)【公開日】2023-03-08

(54)【発明の名称】分析装置、分析方法及びプログラム

(51)【国際特許分類】

   C12M 1/00 20060101AFI20230301BHJP

   G16B 30/10 20190101ALI20230301BHJP

   C12Q 1/6869 20180101ALN20230301BHJP

【ＦＩ】

C12M1/00 A

G16B30/10

C12Q1/6869 Z

【審査請求】未請求

【請求項の数】11

【出願形態】ＯＬ

(21)【出願番号】P 2022194453

(22)【出願日】2022-12-05

(62)【分割の表示】P 2021520435の分割

【原出願日】2020-10-02

(31)【優先権主張番号】P 2019185444

(32)【優先日】2019-10-08

(33)【優先権主張国・地域又は機関】JP

(71)【出願人】

【識別番号】504137912

【氏名又は名称】国立大学法人 東京大学

(71)【出願人】

【識別番号】320008650

【氏名又は名称】株式会社Ｌｉｑｕｉｄ Ｍｉｎｅ

(74)【代理人】

【識別番号】100166006

【弁理士】

【氏名又は名称】泉 通博

(74)【代理人】

【識別番号】100154070

【弁理士】

【氏名又は名称】久恒 京範

(74)【代理人】

【識別番号】100153280

【弁理士】

【氏名又は名称】寺川 賢祐

(72)【発明者】

【氏名】横山 和明

(72)【発明者】

【氏名】清水 英悟

(72)【発明者】

【氏名】近藤 幹也

(57)【要約】

【課題】病気の発生や進行に影響する変異である可能性の程度を提示する。
【解決手段】シーケンスアライメントにより分析の対象となる検体の遺伝情報から抽出された、塩基配列の変異箇所及び変異の内容を含む変異状態を表す変異塩基配列情報を受け入れ、当該受け入れた変異塩基配列情報が表す変異状態ごとに、予め定められた複数の条件のそれぞれを満足するか否かにより設定したスコアを合成し、当該スコアに基づいて、変異状態ごとの病的である可能性の程度を表すランクを推定して出力する分析装置である。複数の変異状態それぞれが、複数の前記条件のいずれの条件を満足したかを表す記録情報をさらに生成し、ランク情報に関連付けて記録情報を出力してもよい。
【選択図】図１

【特許請求の範囲】

【請求項1】

シーケンスアライメントにより分析の対象となる検体の遺伝情報から抽出された、塩基配列の変異箇所及び変異の内容を含む変異状態を表す変異塩基配列情報を受け入れる受入部と、
前記受入部が受け入れた前記変異塩基配列情報が表す変異状態が予め定められた条件を満足するか否かにより決定したスコアを出力するフィルタ処理部と、
前記フィルタ処理部が出力したスコアに基づいて、前記変異状態が病的である可能性の程度を表すランク情報を作成し、作成した前記ランク情報を出力する出力部と、
を含む分析装置。

【請求項2】

前記フィルタ処理部は、複数の前記変異状態それぞれが、複数の前記条件のいずれの条件を満足したかを表す記録情報をさらに生成し、
前記出力部は、前記ランク情報に関連付けて前記記録情報を出力する、
請求項１に記載の分析装置。

【請求項3】

前記受入部が受け入れた前記変異塩基配列情報を保持する保持部をさらに備え、
前記フィルタ処理部は、前記受入部が今回受け入れた前記変異塩基配列情報と、共通の個体の遺伝情報から抽出された、過去に前記受入部が受け入れた前記変異塩基配列情報と、に基づいて得られる、前記変異塩基配列情報の時系列変化に基づく前記条件を満足するか否かにより前記スコアを決定する、
請求項１又は２に記載の分析装置。

【請求項4】

前記フィルタ処理部は、他の遺伝子と融合した融合遺伝子においてドライバー変異候補となる候補遺伝子の塩基配列と類似する塩基配列の遺伝子が他の遺伝子と融合した融合遺伝子が前記変異状態に対応する塩基配列に含まれているか否かにより前記スコアを決定する、
請求項１から３のいずれか一項に記載の分析装置。

【請求項5】

前記フィルタ処理部は、特定の組み合わせの前記候補遺伝子が融合した融合遺伝子においてドライバー変異候補となる２つの当該候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が、前記変異状態に対応する塩基配列に含まれているか否かにより前記スコアを決定する、
請求項４に記載の分析装置。

【請求項6】

前記フィルタ処理部は、異なる生物種のゲノムの間において保存された塩基配列である保存配列の位置が前記変異状態の前記変異箇所に含まれるか否かにより前記スコアを決定する、
請求項１から５のいずれか一項に記載の分析装置。

【請求項7】

前記フィルタ処理部は、前記変異状態が染色体の転座であるか否かにより前記スコアを決定する、
請求項１から６のいずれか一項に記載の分析装置。

【請求項8】

前記フィルタ処理部は、前記変異状態が複数の遺伝子に及ぶ変異であるか否かにより前記スコアを決定する、
請求項１から７のいずれか一項に記載の分析装置。

【請求項9】

前記フィルタ処理部は、前記変異状態が予め登録された登録遺伝子の欠失であるか否かにより前記スコアを決定する、
請求項１から８のいずれか一項に記載の分析装置。

【請求項10】

コンピュータが、
シーケンスアライメントにより分析の対象となる検体の遺伝情報から抽出された、塩基配列の変異箇所及び変異の内容を含む変異状態を表す変異塩基配列情報を受け入れる工程と、
受け入れた変異塩基配列情報が表す変異状態が予め定められた条件を満足するか否かによりスコアを決定する工程と、
決定したスコアに基づいて、前記変異状態が病的である可能性の程度を表すランク情報を作成する工程と、
を実行する分析方法。

【請求項11】

コンピュータに、
シーケンスアライメントにより分析の対象となる検体の遺伝情報から抽出された、塩基配列の変異箇所及び変異の内容を含む変異状態を表す変異塩基配列情報を受け入れる受入ステップと、
当該受け入れた変異塩基配列情報が表す変異状態が予め定められた条件を満足するか否かにより決定したスコアを出力するフィルタステップと、
前記フィルタステップにおいて出力したスコアに基づいて、前記変異状態が病的である可能性の程度を表すランク情報を作成し、作成した前記ランク情報を出力する出力ステップと、
を実行させるプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、分析装置、分析方法及びプログラムに関する。

【0002】

体細胞の遺伝情報に含まれる塩基配列の突然変異により、疾病が生じる場合があることは広く知られている。近年では、種々の体細胞変異に関して、どのような疾病に関係するかといった情報が収集されてデータベースに記録されており、広く利用されている（非特許文献１）。

【0003】

また、近年、網羅的な塩基配列解析技術（例えば次世代シークエンサー）の進歩により、一回の変異の分析で検出される変異は、検体あたり数百から数百万と、膨大な量となっており、当該変異一つ一つにつき、人為的にその結果の解釈を行うのは効率的でなく、現実的でもない。そこで分析結果の人間による解釈を補助する装置が要望されている。

【先行技術文献】

【非特許文献】

【0004】

【非特許文献1】inet：COSMIC v90、[online]、2019年９月５日、[令和１年９月３０日検索]、インターネット<URL: https://cancer.sanger.ac.uk/cosmic>

【発明の概要】

【発明が解決しようとする課題】

【0005】

しかしながら、上記従来のデータベースを利用した塩基配列の変異を分析するだけでは、その検体において上記データベースに記録された変異が生じているか否かを判断できるものの、そのような変異が生じている、という理由だけで登録された症例におけるがんなどの病気の形成や進行に直接影響する変異（例えばがんに対するドライバー変異）であると断定できない。変異の分析結果の解釈には、この他に考慮すべき事項が多岐に亘るためである。ところが、従来は、そのような多岐に亘る項目を組み合わせて、病的である可能性がどの程度であるかを分析することは行われていなかった。

【0006】

本発明は上記実情に鑑みて為されたもので、病気の発生や進行に影響する変異である可能性の程度を提示できる分析装置、分析方法及びプログラムを提供することを、その目的の一つとする。

【課題を解決するための手段】

【0007】

上記従来例の問題点を解決するための本発明の第１の態様の分析装置は、シーケンスアライメントにより分析の対象となる検体の遺伝情報から抽出された、塩基配列の変異箇所及び変異の内容を含む変異状態を表す変異塩基配列情報を受け入れる受入部と、前記受入部が受け入れた前記変異塩基配列情報が表す変異状態が予め定められた条件を満足するか否かにより決定したスコアを出力するフィルタ処理部と、前記フィルタ処理部が出力したスコアに基づいて、前記変異状態が病的である可能性の程度を表すランク情報を作成し、作成した前記ランク情報を出力する出力部と、を含むものである。

【0008】

前記フィルタ処理部は、複数の前記変異状態それぞれが、複数の前記条件のいずれの条件を満足したかを表す記録情報をさらに生成し、前記出力部は、前記ランク情報に関連付けて前記記録情報を出力してもよい。

【0009】

前記分析装置は、前記受入部が受け入れた前記変異塩基配列情報を保持する保持部をさらに備え、前記フィルタ処理部は、前記受入部が今回受け入れた前記変異塩基配列情報と、共通の個体の遺伝情報から抽出された、過去に前記受入部が受け入れた前記変異塩基配列情報と、に基づいて得られる、前記変異塩基配列情報の時系列変化に基づく前記条件を満足するか否かにより前記スコアを決定してもよい。

【0010】

前記フィルタ処理部は、他の遺伝子と融合した融合遺伝子においてドライバー変異候補となる候補遺伝子の塩基配列と類似する塩基配列の遺伝子が他の遺伝子と融合した融合遺伝子が前記変異状態に対応する塩基配列に含まれているか否かにより前記スコアを決定してもよい。

【0011】

前記フィルタ処理部は、特定の組み合わせの前記候補遺伝子が融合した融合遺伝子においてドライバー変異候補となる２つの当該候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が、前記変異状態に対応する塩基配列に含まれているか否かにより前記スコアを決定してもよい。

【0012】

前記フィルタ処理部は、異なる生物種のゲノムの間において保存された塩基配列である保存配列の位置が前記変異状態の前記変異箇所に含まれるか否かにより前記スコアを決定してもよい。

【0013】

前記フィルタ処理部は、前記変異状態が染色体の転座であるか否かにより前記スコアを決定してもよい。前記フィルタ処理部は、前記変異状態が複数の遺伝子に及ぶ変異であるか否かにより前記スコアを決定してもよい。

【0014】

前記フィルタ処理部は、前記変異状態が予め登録された登録遺伝子の欠失であるか否かにより前記スコアを決定してもよい。

【0015】

本発明の第２の実施態様の分析方法は、コンピュータが、シーケンスアライメントにより分析の対象となる検体の遺伝情報から抽出された、塩基配列の変異箇所及び変異の内容を含む変異状態を表す変異塩基配列情報を受け入れる工程と、受け入れた変異塩基配列情報が表す変異状態が予め定められた条件を満足するか否かによりスコアを決定する工程と、決定したスコアに基づいて、前記変異状態が病的である可能性の程度を表すランク情報を作成する工程と、を実行する。

【0016】

本発明の第３の態様のプログラムは、コンピュータに、シーケンスアライメントにより分析の対象となる検体の遺伝情報から抽出された、塩基配列の変異箇所及び変異の内容を含む変異状態を表す変異塩基配列情報を受け入れる受入ステップと、当該受け入れた変異塩基配列情報が表す変異状態が予め定められた条件を満足するか否かにより決定したスコアを出力するフィルタステップと、前記フィルタステップにおいて出力したスコアに基づいて、前記変異状態が病的である可能性の程度を表すランク情報を作成し、作成した前記ランク情報を出力する出力ステップと、を実行させる。

【発明の効果】

【0017】

本発明によると、病気の発生や進行に影響する変異である可能性の程度を提示できる。

【図面の簡単な説明】

【0018】

【図1】本発明の実施の形態に係る分析装置の構成例を示すブロック図である。

【図2】本発明の実施の形態に係る分析装置の例を表す機能ブロック図である。

【図3】本発明の実施の形態に係る分析装置に入力される変異塩基配列情報の例を表す説明図である。

【図4】本発明の実施の形態に係る分析装置が出力する出力情報の一例を表す説明図である。

【図5】本発明の実施の形態に係る分析装置の動作例を表す流れ図である。

【図6】本発明の実施の形態に係る分析装置が出力する出力情報のもう一つの例を表す説明図である。

【図7】第２の実施形態の制御部の構成を示す図である。

【図8】分析装置が病原性の疑いがあるか否かを判定する処理手順を示すフローチャートである。

【発明を実施するための形態】

【0019】

＜第１の実施形態＞
本発明の第１の実施形態について図面を参照しながら説明する。本発明の実施の形態に係る分析装置１は、図１に例示するように、制御部１１、記憶部１２、操作部１３、表示部１４、入出力部１５、及び通信部１６を含む一般的なコンピュータ装置を用いて実現される。

【0020】

ここで制御部１１は、ＣＰＵ等のプログラム制御デバイスであり、記憶部１２に格納されたプログラムに従って動作する。本実施の形態では、この制御部１１は、分析の対象となる個体と、個体から得た検体とを識別する検体識別情報とともに、当該検体の遺伝情報からシーケンスアライメントにより抽出された、塩基配列の変異箇所及び変異の内容を含む変異状態を表す変異塩基配列情報を受け入れる。変異状態は、一塩基の変異であってもよく、複数の遺伝子に及ぶ染色体の転座等の構造変異であってもよい。具体的に変異箇所及び変異の内容とは、変異のあった位置（遺伝子と、その遺伝子上の一方側から何番目の塩基であるかを示す情報等）と、あるべき塩基がどの塩基に変異しているのかを表す情報を含むものである。

【0021】

そして制御部１１は、当該受け入れた変異塩基配列情報が表す変異状態ごとに、予め定められた複数の条件のそれぞれを満足するか否かにより、受け入れた変異塩基配列情報に含まれる変異状態ごとにスコアを設定し、さらに当該変異状態ごとに、設定したスコアを合成した合成スコア情報を生成する。

【0022】

制御部１１は、この変異状態ごとに得た合成スコア情報に基づいて、変異状態ごとの病的である可能性の程度を表すランクを推定し、当該推定したランクの情報を出力する。この制御部１１の動作については後に詳しく述べる。

【0023】

記憶部１２は、メモリデバイスやディスクデバイス等であり、制御部１１によって実行されるプログラムを保持する。この記憶部１２は、また、制御部１１のワークメモリとしても動作する。

【0024】

操作部１３は、キーボードやマウス等であり、ユーザからの操作を受け入れて、当該操作を表す情報を制御部１１に出力する。表示部１４は、ディスプレイ等であり、制御部１１から入力される指示に従って情報を表示する。

【0025】

入出力部１５は、例えばＵＳＢ（Universal Serial Bus）インタフェース等であり、他の機器で生成された情報を受け入れて制御部１１に出力する等の処理を行う。通信部１６は、ネットワークインタフェースであり、制御部１１から入力される指示に従って、ネットワークを介して接続された他のコンピュータ装置との間で通信してデータの送受を行う。

【0026】

次に、制御部１１の具体的な動作の内容について説明する。本実施の形態では、この制御部１１は、記憶部１２に格納されたプログラムを実行することで、図２に例示する構成を機能的に実現する。すなわち制御部１１は、データ受入部２１と、設定受入部２２と、フィルタ処理部２３と、ランク推定部２４と、出力部２５とを含む構成を実現する。

【0027】

ここでデータ受入部２１は、分析の対象となる検体の遺伝情報からシーケンスアライメントにより抽出された、塩基配列の変異状態を表す変異塩基配列情報を受け入れる。この情報は具体的には図３に例示するように、変異状態ごとに、当該変異状態の塩基配列が見出された染色体の番号（Chr）と、開始位置（Start）と、終了位置（End）と、本来あるべき塩基配列（Ref）と、抽出された変異している塩基配列（Alt）と、変異している塩基配列の割合（アレル頻度：ＡＦ）とを関連付けたものである。ここでの例ではさらに、これらの情報に、深度（depth）や変異状態のカウント数（AltCount）など、品質に関係する指標等を関連付けたものとしている。なお、塩基配列の長さは「１」（この場合、塩基配列の情報はＡ，Ｔ，Ｃ，Ｇのいずれかの塩基を表す情報となる）であってもよい。

【0028】

またこの変異塩基配列情報には、個体の症例等に関する情報（病名や治療歴、腫瘍割合などの情報）が含まれてもよい。

【0029】

本実施の形態の一例では、このデータ受入部２１は、同じ個体から、分析の対象とする変異塩基配列情報を抽出したタイミングとは異なるタイミング（複数あってもよい）で抽出された変異塩基配列情報（時系列情報）を受け入れてもよい。この場合は、データ受入部２１は、分析の対象とする時点の変異塩基配列情報の指定の入力を受けておく。

【0030】

設定受入部２２は、分析の設定を受け入れる。この設定は、例えば、どのフィルタを用いるかという設定、又はパラメータ情報の設定である。この設定の具体例は、次のフィルタ処理部２３の構成とともに説明する。

【0031】

本実施の形態では、このフィルタ処理部２３の動作により、変異の分析結果の解釈に影響する種々の情報に基づく、病的である可能性に関する評価が行われる。ここで、解釈に影響する情報には、（１）分析の際に得られる当該変異の付帯情報、及び（２）文献やデータベースに収載された変異に関連する情報、が含まれる。このうち（１）分析の際に得られる当該変異の付帯情報には、（ａ）検出精度や信頼性の情報（変異が検出エラーでない確率）、（ｂ）変異のアレル頻度（同一変異を持つ細胞集団の全体に占める割合に関連する指標）、（ｃ）時系列情報、すなわち同一症例の他の時点での検体において繰り返して当該変異が検出されているか否か、などがある。

【0032】

また、（２）文献やデータベースに収載された変異に関連する情報には、当該変異が疾患のドライバー変異として記載されているか否か（あるいはどの程度の頻度で記述がされているか）を表す情報が含まれる。ＳＮＰ（single nucleotide polymorphism）データベースにも登録がある場合、変異アレルはどの程度の対立アレル頻度で、当該人種においてＳＮＰとして報告されているかという情報が文献やデータベースに収載されていてもよい。さらには、機能予測として、当該変異がコード化されたたんぱく質の立体構造や機能に影響を与えるか、例えばがんの病態形成に関わるとして実験などにより示されているか予測されているか否かを表す情報などが文献やデータベースに収載されていてもよい。

【0033】

フィルタ処理部２３は、データ受入部２１が受け入れた変異状態（時系列情報を受け入れたときには、そのうち分析の対象として指定された変異塩基配列情報に含まれる変異状態、以下、分析の対象とする変異状態と呼ぶ）ごとに、予め定められた複数の条件を満足するか否かを判断する。具体的に本実施の形態では、図２に例示したように、フィルタ処理部２３は、基本フィルタ２３１と、時系列フィルタ２３２と、データベースフィルタ２３３と、機能予測フィルタ２３４と、クオリティフィルタ２３５とを含んで構成される。

【0034】

ここで基本フィルタ２３１は、分析の対象とする変異状態が、良性であると判断できる場合に、良性変異であることを表すスコア（例えば「４」）を設定して、ランク推定部２４にその結果を出力する。また、基本フィルタ２３１は、分析の対象とする変異状態が、良性であると判断できなければ、良性変異でないことを表すスコア（例えば「３」）を設定し、次のフィルタとして設定されたフィルタに処理を受け渡す（処理をパスする）。

【0035】

ここで良性と判断できる場合とは、がん化等を引き起こす既知の変異の塩基配列と変異状態に対応する変異した塩基配列との重複部分が比較的短い重複部分である場合、変異状態が表す変異の位置する領域がイントロン領域である場合、ＳＮＰデータベースなどのように異常がない変異を蓄積したデータベースに変異状態が登録されている場合、又はＧＤＩ（Gene Damage Index）に基づいて変異状態が良性と判断できる場合等が相当する。

【0036】

ここでＧＤＩは、各遺伝子について健常人にどれだけダメージが蓄積されているかを表す指標であり、人によって大きいダメージを受けていても（多様性があっても）、変異によって病的になると考えられない遺伝子である可能性を示す。

【0037】

基本フィルタ２３１は、設定受入部２２から、がん化等を引き起こす既知の変異の塩基配列と変異状態に対応する変異した塩基配列との重複部分の長さの閾値、ＳＮＰであるか否かを判断するデータベースを特定する情報、データベースごとのパラメータ（良性か否かを判断する基準となる良性判断閾値、又はデータベースにＳＮＰである確率などとして登録された値と比較される）の少なくともいずれかの設定を受け入れる。基本フィルタ２３１は、受け入れた設定に基づいて、分析の対象とする変異状態が良性であるか否かを判断する。

【0038】

具体的に基本フィルタ２３１は、がん化等を引き起こす既知の変異の塩基配列と変異状態に対応する変異した塩基配列との重複部分が予め定めた長さの閾値より短い重複部分である場合は良性変異であることを表すスコアを設定する。また基本フィルタ２３１は、重複部分が予め定めた長さの閾値より短い重複部分である場合は良性変異でなくても、変異状態が表す変異の位置する領域がイントロン領域であれば、良性変異であることを表すスコアを設定する。

【0039】

さらに基本フィルタ２３１は、上記２つの条件を満足しなくても、指定されたＳＮＰデータベースを検索した結果に基づいて、良性変異であることを表すスコアを設定してもよい。基本フィルタ２３１は、例えば、検索によって変異状態が表す変異がＳＮＰデータベースに登録され、かつ、そのＳＮＰである確率として登録された値が、当該ＳＮＰデータベースについて予め定められた良性判断閾値を超えている場合、良性変異であることを表すスコアを設定する。

【0040】

また基本フィルタ２３１は、ここまでの条件を満足しない場合であっても、当該変異状態が存在する遺伝子のＧＤＩを参照して、予め定められたＧＤＩ閾値より大きい場合に、良性変異であることを表すスコアを設定する。

【0041】

これにより本実施の形態では、例えばがんのドライバー変異となり得ない（あるいはその可能性が十分低い）遺伝子を予めふるいわけすることが可能となる。

【0042】

またこの基本フィルタ２３１では、設定受入部２２から予め定めた、良性と判断するための複数の条件のうち、どの条件を利用するか（あるいはすべての条件を利用せず、基本フィルタ２３１としての動作をせずに、すべての変異状態についてスコアを「３」に設定して処理をパスするか否か）の設定を受け入れてもよい。

【0043】

この例では基本フィルタ２３１は、利用すると設定された条件に限って、当該条件を満足するか否かを判断することとなる。

【0044】

時系列フィルタ２３２は、基本フィルタ２３１が処理をパスした（スコア「３」が設定された）場合に、分析の対象とする変異状態に対応する、時系列情報に含まれる変異状態の情報を参照して、異なるタイミングで抽出した時系列情報においても同じ変異があったか否かを判断する。

【0045】

時系列フィルタ２３２は、分析の対象とする変異状態と、時系列情報に含まれる対応する変異状態とを用い、同じ変異が存在する場合に、問題とするべき変異があるものとしてスコア（例えば現在のスコアから第１所定量として「１」を引く）を設定してクオリティフィルタ２３５に処理をパスする。第１所定量は、例えば、変異状態にかかるスコアから１回の演算において減算又は加算される最小値である。ここでの例では基本フィルタ２３１が処理をパスしているので、当初のスコアは「３」であり、ここで時系列フィルタ２３２が問題とするべき変異があるものとしたときには、このスコア「３」から第１所定量として「１」を引いてスコアを「２」と設定することとなる。

【0046】

一方、時系列フィルタ２３２は、分析の対象とする変異状態と、時系列情報に含まれる対応する変異状態とを用い、同じ変異が存在しないときには、スコアをそのままに設定（ここでは当初のスコアが「３」であるので、そのまま「３」に設定）して、データベースフィルタ２３３へ処理をパスする。

【0047】

なお、時系列フィルタ２３２は、設定受入部２２から深度や、その他のシーケンスクオリティ、変異アレル頻度等に関する閾値の設定を受けていてもよい。例えば時系列情報に含まれる対応する変異状態に係る深度がここで設定された閾値（例えば「２０」）を超えない場合には、時系列フィルタ２３２は、同じ変異状態があったか否かを判断することなく、スコアをそのままに設定（ここでは当初のスコアが「３」であるので、そのまま「３」に設定）して、データベースフィルタ２３３へ処理をパスする。

【0048】

さらに本実施の形態の例では、この時系列フィルタ２３２は、データ受入部２１が時系列情報を受け入れていない場合（変異塩基配列情報としては、分析の対象とする変異塩基配列情報だけを受け入れている場合）には、同じ変異状態があるか否かを判断することなく、スコアをそのままに設定（ここでは当初のスコアが「３」であるので、そのまま「３」に設定）して、データベースフィルタ２３３へ処理をパスしてもよい。

【0049】

また設定受入部２２から時系列フィルタ２３２を利用しない設定が入力されている場合、時系列フィルタ２３２は、同じ変異状態があるか否かを判断することなく、スコアをそのままに設定（ここでは当初のスコアが「３」であるので、そのまま「３」に設定）して、データベースフィルタ２３３へ処理をパスする。

【0050】

データベースフィルタ２３３は、分析の対象とする変異状態が、予め定めた、問題とするべき変異に関する情報を蓄積したデータベース（例えばCOSMIC Cancer Databaseなど）に登録されているか否かを、当該データベースのサーバへ変異状態に関する情報を送信して調べ、登録されている場合に、問題とするべき変異があるものとしてスコア（例えば現在のスコアから第１所定量として「１」を引く）を設定して、クオリティフィルタ２３５に処理をパスする。ここでの例では、分析の対象とする変異状態について、基本フィルタ２３１が処理をパスし、さらに時系列フィルタ２３２においてスコアをそのままとした状態で処理がパスされたときに、このデータベースフィルタ２３３が判断をすることとなるので、データベースフィルタ２３３は、その際のスコア「３」から第１所定量として「１」を引いて、スコアを「２」に設定した上で、クオリティフィルタ２３５に処理をパスすることとなる。

【0051】

またこのデータベースフィルタ２３３は、分析の対象とする変異状態が、上記問題とするべき変異に関する情報を蓄積したデータベースに登録されていなかった場合に、スコアをそのままに設定して機能予測フィルタ２３４に処理をパスする。ここでの例では、このときのスコアは「３」のままとなる。

【0052】

なお、このデータベースフィルタ２３３は、上記問題とするべき変異に関する情報を蓄積したデータベースとして、どのようなデータベースを利用するかの設定を、設定受入部２２から受け入れておく。

【0053】

この設定では、複数のデータベースを用いるべき旨の指示がされてもよく、この場合、データベースフィルタ２３３は、分析の対象とする変異状態が、上記問題とするべき変異に関する情報を蓄積したデータベースのいずれかに登録されていた場合に、問題とするべき変異があるものとしてスコアを設定する。

【0054】

機能予測フィルタ２３４は、変異の病原性を評価したデータベースを参照して、分析の対象とする変異状態に係る変異が、病原性のあるものとしてデータベースに登録されている場合に、病原性の変異があるものとしてスコア（例えば現在のスコアから第１所定量として「１」を引く）を設定して、クオリティフィルタ２３５に処理をパスする。

【0055】

ここで変異の病原性を評価したデータベースとしては、ＳＩＦＴや、PolyPhen2など広く知られたものがある。また、これらのデータベースのうちには、病原性の有無について多段階で評価しているものがあるが、このような場合は、この機能予測フィルタ２３４は、例えば病原性の疑いがあるとの判断段階にある場合に、病原性の変異があるものとしてスコア（例えば現在のスコアから第１所定量として「１」を引く）を設定して、クオリティフィルタ２３５に処理をパスするものとする。

【0056】

ここでの例では、分析の対象とする変異状態について、基本フィルタ２３１が処理をパスし、時系列フィルタ２３２においてスコアをそのままとした状態で処理がパスされ、さらにデータベースフィルタ２３３においてもスコアがそのままの状態で処理がパスされたときに、この機能予測フィルタ２３４が判断をすることとなるので、機能予測フィルタ２３４は、その際のスコア「３」から第１所定量として「１」を引いて、スコアを「２」に設定した上で、クオリティフィルタ２３５に処理をパスすることとなる。

【0057】

またこの機能予測フィルタ２３４は、変異の病原性を評価したデータベースを参照して、分析の対象とする変異状態に係る変異が、病原性のあるものとしてデータベースに登録されていなければ（または登録されていても不明であるとか、良性ないし良性と推定される場合として登録されている場合）、スコアをそのままに設定してクオリティフィルタ２３５に処理をパスする。ここでの例では、このときのスコアは「３」のままとなる。

【0058】

なお、この機能予測フィルタ２３４においても、どのようなデータベースを利用するかの設定を、設定受入部２２から受け入れておくものとする。

【0059】

クオリティフィルタ２３５は、分析の対象とする変異状態をシーケンスしたときの深度や、その他の、分析の対象とする変異状態のシーケンス処理の品質を評価する。この品質の指標については、深度のほか、変異状態のカウント数など、広く知られた指標があり、クオリティフィルタ２３５は、これらを組み合わせて（あるいはその組み合わせを、設定受入部２２から受け入れ、当該受け入れた指標の組み合わせに従って）品質の評価を行う。なお、クオリティフィルタ２３５は、複数の指標を組み合わせる場合は、すべての指標により、品質が十分高いとの条件を満足する場合に、品質が十分であると判断することとする。

【0060】

クオリティフィルタ２３５は、この評価により、分析の対象とする変異状態のシーケンス処理の品質が十分である（十分高い）と判断したときに、判断が適性であるとしてスコア（例えば現在のスコアから第１所定量として「１」を引く）を設定して、ランク推定部２４に当該スコアを出力する。またこのクオリティフィルタ２３５は、分析の対象とする変異状態のシーケンス処理の品質が十分である（十分高い）と判断できないときには、スコアをそのままに設定してランク推定部２４に当該スコアを出力する。

【0061】

本実施の形態において特徴的なことの一つは、このように、フィルタ処理部２３が、複数の条件に対応するフィルタを、予め定められた順序で適用しつつ、変異状態ごとに、各条件を満足するか否かの判断結果に基づくスコアを設定することである。

【0062】

ランク推定部２４は、フィルタ処理部２３が出力する、変異状態ごとのスコアに従い、変異状態ごとの病的である可能性の程度を表すランク値を得て、当該ランク値を対応する変異状態を特定する情報に関連付けて出力する。この病的である可能性の程度を表すランク値は、各変異状態に対応してフィルタ処理部２３が出力したスコアの値そのものであってもよい。

【0063】

出力部２５は、ランク推定部２４が出力する、変異塩基配列情報に含まれる変異状態のそれぞれを特定する情報と、それぞれに関連付けられたランク値とを集約して出力する。出力される情報（出力情報）における、変異状態のそれぞれを特定する情報は、例えば図４に例示するように、データ受入部２１が受け入れた、変異塩基配列情報に含まれる変異状態の塩基配列が見出された染色体の番号（Chr）と、開始位置（Start）と、終了位置（End）と、本来あるべき塩基配列（Ref）と、抽出された変異している塩基配列（Alt）とであってもよい。この場合、出力情報は、図４に例示するように、これらの変異状態のそれぞれを特定する情報と、ランク値（Rank）とを関連付けた行が配列されたものとなる。

【0064】

［動作］
本実施の形態は以上の構成を備えており、基本的に次のように動作する。本実施の分析装置１は、例えば入出力部１５に接続された記憶デバイスから、図３に例示したような、分析の対象となる検体（被験者Ａ）の遺伝情報からシーケンスアライメントにより抽出された、塩基配列の変異状態を表す変異塩基配列情報を読み出して受け入れる。

【0065】

本実施の形態のここでの例では、同じ個体である被験者Ａから、分析の対象とする変異塩基配列情報を抽出した日時（例えば病気治療後の日時）より前の日時（例えば病気治療前の日時）で抽出された変異塩基配列情報を時系列情報として受け入れるものとする。

【0066】

この分析装置１のユーザは、病気治療後の日時に抽出された変異塩基配列情報を分析の対象として設定する操作を行い、分析装置１は、この操作を受け入れて、分析の対象とする変異塩基配列情報と、時系列情報とを識別する。

【0067】

またユーザは、分析の設定を行う。本実施の形態のここでの例では、病的である可能性の程度であるランクを決定するための条件を判断する処理として、
（１）良性変異を登録したデータベースに登録があるか等を調べて条件を判断する基本フィルタ、
（２）時系列情報を参照した条件の判断を行う時系列フィルタ、
（３）問題のある変異（病原性のある変異）であることを登録したデータベースを参照して条件を判断するフィルタ（データベースフィルタ）、
（４）変異の病原性を評価したデータベースを参照して条件を判断する機能予測フィルタ、
（５）シーケンスの品質を評価して条件を判断するクオリティフィルタ
を用いるものとする。

【0068】

そこで、ここでは、（１）基本フィルタについては、良性と判断するための公知の条件のうち、どのような条件を用いるかを設定受入部２２が設定する。具体的には、本実施形態においては、がん化等を引き起こす既知の変異の塩基配列と変異状態に対応する変異した塩基配列との重複部分の長さの閾値（これより短い場合良性と判断する）αと、良性変異を登録したデータベースとしてどのデータベースを利用するかを設定受入部２２が設定するものとする。

【0069】

また、（２）時系列フィルタについては、例えば過去または今回の変異状態に係る深度や変異アレル頻度等に関する閾値（ここでは深度及び変異アレル頻度の閾値とする）β，γを設定する。そして、（３）データベースフィルタ及び（４）機能予測フィルタについては、それぞれにおいてどのデータベースを利用するかを設定受入部２２が設定する。この設定はそれぞれのフィルタに対応して予め定められた選択肢から、それぞれのフィルタで使用するデータベースを選択することで行う。

【0070】

（５）クオリティフィルタについては、広く知られた指標のうち、どの指標を使用するかを、予め定められた選択肢のうちから選択するとともに、それぞれの指標における閾値を設定受入部２２が設定する。以下の説明の例では、深度と、変異状態のカウント数とを用いるものとし、深度については「２０」、変異状態のカウントは「１０」と設定されたものとする。

【0071】

分析装置１は、処理開始の指示をユーザから受けると、図５に例示するように、分析の対象とする変異塩基配列情報に含まれる、各変異状態を順次選択する（Ｓ１１）。分析装置１は、選択した変異状態について、基本フィルタの表す条件により、病原性の疑いがあるか否かを判断する（Ｓ１２）。具体的には、ここでの例では、がん化等を引き起こす既知の変異の塩基配列と変異状態に対応する変異した塩基配列との重複部分の長さの閾値αと、良性変異を登録したデータベースとしてどのデータベースを利用するかが設定されているので、分析装置１は、このステップＳ１２では、選択した変異状態のがん化等を引き起こす変異と変異状態に対応する変異した塩基配列との重複の長さがαを超えるか否かを判断するとともに、良性変異を登録したデータベースに登録されているか否かを判断する。

【0072】

そして分析装置１は、選択した変異状態のがん化等を引き起こす変異と変異状態に対応する変異した塩基配列との重複の長さがαを超え、かつ、良性変異を登録したデータベースに登録されていなければ、病原性の疑いがあるとして（Ｓ１２：Ｙｅｓ）、この選択した変異状態に係るスコアを仮に「３」とし、ステップＳ１３に移行する。

【0073】

すなわちステップＳ１２で病原性の疑いがあるとされると、分析装置１は、時系列フィルタの表す条件により、病原性の疑いがあるか否かを判断する（Ｓ１３）。具体的に分析装置１は、ステップＳ１１で選択した変異状態と、当該選択された変異状態に対応して、時系列情報に同じ変異状態があり（当該対応する情報があり）、かつ、選択された変異状態または当該時系列情報に記録された、対応する情報に係る深度の情報が、予め定められた閾値βを超え、かつ、選択された変異状態のアレル頻度の値と、対応する時系列情報に係るアレル頻度の値とのうち少なくともいずれか一つが設定された閾値γを上回る場合に、病原性の疑いがあると判断する（Ｓ１３：Ｙｅｓ）。分析装置１は、この時点で選択された変異状態に係るスコアとして仮に設定されたスコアから「１」を差し引いて、選択された変異状態に係るスコアを仮に「２」に設定し、ステップＳ１４に移行する。

【0074】

分析装置１は、ステップＳ１３において病原性の疑いがあると判断されると、クオリティフィルタの表す条件により、ステップＳ１１で選択した変異状態のシーケンス処理の品質を評価する（Ｓ１４）。

【0075】

具体的にこのステップＳ１４において分析装置１は、選択された変異状態に関連付けて、分析の対象とする変異塩基配列情報に記録されている、品質に関連する情報（ここでの例では深度と変異状態のカウント）を読み出して、それぞれが、設定された閾値と比較して、品質が十分であるか否かを調べる（ここでの例ではいずれも閾値を超える場合に品質が十分高いものと判断される）。

【0076】

分析装置１は、品質が十分であると判断すると（Ｓ１４：Ｙｅｓ）、この時点で選択された変異状態に係るスコアとして仮に設定されたスコアから「１」を差し引いて、選択された変異状態に係るスコアを設定して、選択された変異状態に対応するランク値を、当該設定したスコアとして、選択された変異状態を特定する情報に関連付けて記録する（Ｓ１５）。

【0077】

具体的な例として、分析装置１が、ステップＳ１３で病原性の疑いがあると判断し、選択された変異状態に係るスコアを「２」と設定した状態でステップＳ１４を実行して、十分な品質であると判断すると、スコアは「１」（ランク値は「１」）となる。

【0078】

また分析装置１が後に述べる動作により、変異状態に係るスコアを「３」と設定した状態でステップＳ１４を実行して、十分な品質であると判断すると、スコアは「２」（ランク値は「２」）となる。

【0079】

一方、分析装置１は、品質が十分でないと判断する場合（ここでの例では、品質に関連する深度と変異状態のカウントのいずれかが、設定された閾値を超えない場合）、この時点で選択された変異状態に係るスコアとして仮に設定されたスコアをそのままとして、ステップＳ１５に移行し、選択された変異状態を特定する情報に、当該スコアをランク値として関連付けて記録する。

【0080】

具体的な例として、分析装置１が、ステップＳ１３で病原性の疑いがあると判断し、選択された変異状態に係るスコアを「２」と設定した状態でステップＳ１４を実行して、十分な品質でないと判断すると、スコアは「２」（ランク値は「２」）となる。

【0081】

また分析装置１が後に述べる動作により、変異状態に係るスコアを「３」と設定した状態でステップＳ１４を実行して、十分な品質でないと判断すると、スコアは「３」（ランク値は「３」）となる。

【0082】

また分析装置１は、ステップＳ１３において、選択された変異状態または時系列情報において対応する情報に係る深度の情報が、予め定められた閾値βを下回る場合、あるいは、選択された変異状態のアレル頻度の値と、対応する時系列情報に係るアレル頻度の値とがともに、設定された閾値γを下回る場合に、病原性の疑いがないとして（Ｓ１３：Ｎｏ）、この時点で選択された変異状態に係るスコアとして仮に設定されたスコアをそのまま（ここでの例では「３」）として、ステップＳ１６に移行する。

【0083】

分析装置１は、ステップＳ１３において病原性の疑いがないと判断されると、データベースフィルタの表す条件により、ステップＳ１１で選択した変異状態に病原性の疑いがあるか否かを判断する（Ｓ１６）。

【0084】

すなわち分析装置１は、利用するデータベースとして設定されたデータベースに、通信部１６を介してアクセスし、ステップＳ１１で選択した変異状態が当該データベースに登録されているか否かを問い合わせる。ここで当該問い合わせの結果、選択した変異状態が当該データベースに登録されている旨の応答があった場合（Ｓ１６：Ｙｅｓ）は、この時点での、選択した変異状態に係るスコア「３」から「１」を差し引いて、選択された変異状態に係るスコアを「２」に設定し、ステップＳ１４に移行して処理を続ける。

【0085】

一方、分析装置１は、ステップＳ１６において、問い合わせの結果、選択した変異状態が当該データベースに登録されていない旨の応答があった場合（Ｓ１６：Ｎｏ）、病原性の疑いがないと判断して、この時点での、選択した変異状態に係るスコア「３」をそのままとしてステップＳ１７に移行する。

【0086】

分析装置１は、ステップＳ１６において、病原性の疑いがないと判断したときには、機能予測フィルタの表す条件により、ステップＳ１１で選択した変異状態に病原性の疑いがあるか否かを判断する（Ｓ１７）。

【0087】

すなわち分析装置１は、機能予測フィルタに関して利用するデータベースとして設定されたデータベースに、通信部１６を介してアクセスし、ステップＳ１１で選択した変異状態が当該データベースにおいてどのような病原性の程度として登録されているかを問い合わせる。ここで当該問い合わせの結果、選択した変異状態に係る病原性の評価が予め定めた基準に比べて高いものであるときには（Ｓ１７：Ｙｅｓ）、分析装置１は、この時点での、選択した変異状態に係るスコア「３」から「１」を差し引いて、選択された変異状態に係るスコアを「２」に設定し、ステップＳ１４に移行して処理を続ける。

【0088】

また分析装置１は、ステップＳ１７において、問い合わせの結果、選択した変異状態が当該データベースに登録されていない旨の応答があった場合、あるいは、選択した変異状態に係る病原性の評価が予め定めた基準に比べて高いものでなかったときには（Ｓ１７：Ｎｏ）、病原性の疑いがないと判断して、この時点での、選択した変異状態に係るスコア「３」をそのままとしてステップＳ１４に移行する。

【0089】

また分析装置１は、ステップＳ１２において、選択した変異状態においてがん化等を引き起こす変異と変異状態に対応する変異した塩基配列との重複の長さがαを超えない場合、または良性変異を登録したデータベースに登録されている場合に、病原性の疑いがない（良性変異である）と判断する（Ｓ１２：Ｎｏ）。分析装置１は、選択した変異状態に係るスコアを「４」として、ステップＳ１５に移行する。この場合、分析装置１は、選択された変異状態を特定する情報に当該スコア「４」をランク値として関連付けて記録する。

【0090】

分析装置１は、分析の対象とする変異塩基配列情報に、未だ選択されていない変異状態があれば、ステップＳ１１に戻って当該未選択のものから一つを選択して処理を続ける。一方、未選択のものがなければ（分析の対象とする変異塩基配列情報に含まれるすべての変異状態についてランク値が関連付けられたならば）、各変異状態を特定する情報とそのランク値とを関連付けた、出力情報を出力する。この出力情報は、例えば図４に例示したものとなる。

【0091】

ユーザは、この出力情報を参照し、例えばランク値が「１」または「２」となる変異に着目して病的であるか否かを判断することができる。またここでの例では、ランク値が「３」である場合は、シーケンスエラーの可能性があることから、ユーザは、ランク値が「３」の変異が多い場合は再検査の必要性などを判断できる。

【0092】

なお、ここで例示した分析装置１の動作の流れは一例であり、その順序はユーザが任意に定め得る。例えばステップＳ１７とＳ１６とは逆順であってもよい。また、いずれかのフィルタに係るステップの処理は省略されてもよい。

【0093】

［判断に関する記録］
また本実施の形態の例において、分析装置１の制御部１１は、フィルタ処理部２３の動作において、基本フィルタ２３１と、時系列フィルタ２３２と、データベースフィルタ２３３と、機能予測フィルタ２３４と、クオリティフィルタ２３５とのうち、変異状態ごとに実行したフィルタの動作（複数の条件の判断）において、どの条件がどのように判断されたか（いずれの条件を満足したかなど）を表す記録情報をさらに生成して、変異状態ごとのランク値とともに、併せて出力情報に含めて出力してもよい。

【0094】

例えば制御部１１は、基本フィルタ２３１として機能する際に、例えばＳＮＰデータベース（異常がない変異を蓄積したデータベースの例）に登録されていると判断すると、スコアを「４」と設定するとともに、ＳＮＰデータベースに記録があった旨の記録情報を出力する。そして出力部２５として機能するときに、出力情報において、スコアの値に基づくランク値とともに変異状態を特定する情報に関連付けて当該記録情報を出力する。

【0095】

この場合は、当該変異状態については他のフィルタとしての機能は実行されないため、当該実行されない他のフィルタでの記録情報は含まれないこととなる。

【0096】

また、制御部１１は、基本フィルタ２３１として機能したときに、良性であると判断できず、スコアを「３」に設定して処理をパスするときには、基本フィルタ２３１としては記録情報を生成しなくてもよい。あるいは、スコアを設定するとともに、ＳＮＰデータベースに記録がなかったことなど、良性と判断できなかった理由となる記録情報を出力してもよい。このときには出力部２５として機能するときに、出力情報において、スコアの値に基づくランク値とともに変異状態を特定する情報に関連付けてこの記録情報を出力する。

【0097】

そして上述の例では、このときには、制御部１１は、時系列フィルタ２３２としての機能も実行することとなるので、分析の対象とする変異状態のアレル頻度と設定された閾値とを比較してどのような結果であったかなどの記録情報を生成して出力する。

【0098】

このような例によると、例えば（図４に対応する）図６に例示するように、染色体の番号（Chr）と、開始位置（Start）と、終了位置（End）と、本来あるべき塩基配列（Ref）と、抽出された変異している塩基配列（Alt）とを含む、変異状態を特定する情報に関連付けて、ランク値（Rank）と、判断の理由となる記録情報（図６の例では、判断の理由ごとにカラムを分けて記録した例としている；Ｒ）を出力する。

【0099】

ユーザは、ランク値のみならず、この判断の理由となる記録情報を併せて参照することで、例えば病原性ありと評価されたときに（ランク値が「１」であるときに）、その判断の理由が、主に時系列フィルタの判断に基づくものであるか、データベースフィルタ等の判断に基づくものであるかを区別して扱うことができる。このことは、変異の性質の理解に役立つ。

【0100】

［回帰性の分析］
また、第１の実施形態の分析装置１は、変異塩基配列情報の分析処理を行うごとに（検体，個体によらずに）、当該変異塩基配列情報に含まれる変異状態ごとの出現頻度を累算して記録してもよい。またこの分析装置１は、当該記録した変異状態ごとの出現頻度をユーザに提示してもよい。例えば、分析の対象とした変異塩基配列情報に含まれる変異状態ごとに、上記記録した、当該変異状態の出現頻度の累算値を出力してもよい。

【0101】

この累算値の情報を参照すると、ユーザは、病的である可能性の程度を表すランク値が高くない変異状態において、検体，個体によらずに数多く出現する変異状態であるか否かが判断でき、検体，個体によらずに数多く出現する変異状態であれば、次世代シークエンサーの読み取り誤りが多く発生しやすい部位であるなどといった判断が可能となる。

【0102】

またユーザは、この累算値の情報を参照して、病的である可能性のある変異状態において、検体，個体によらずに数多く出現する変異状態であるか否かが判断でき、検体，個体によらずに数多く出現する変異状態であれば、当該変異状態の登録がデータベースに登録されていない場合でも、いままで知られていなかった病的な変異である可能性を知ることができる。

【0103】

また、分析装置１は、変異塩基配列情報に含まれる症例の情報（例えば病名）ごとに、検体，個体によらない変異状態ごとの出現頻度の累算値をさらに得て記録してもよい。また、変異状態ごとの出現頻度の累算値として、さらに変異状態のランク値ごとの出現頻度の累算値を得て記録してもよい。

【0104】

さらに、変異塩基配列情報に、シーケンスを行った次世代シークエンサーのメーカー及び機種名の情報を含めておき、当該シークエンサーのメーカーや機種名の情報ごとに上記出現頻度の集計を行ってもよい。この場合、シークエンサーのメーカーや機種ごとのシーケンスエラーの癖を知ることができる。

【0105】

このように第１の実施形態によると、病気の発生や進行に影響する変異である可能性、すなわち病的である可能性の程度を提示でき、各変異が、病気の発生や進行に直接影響する変異（例えばがんに対するドライバー変異）である可能性がどの程度であるかといった判断に資することができる。

【0106】

＜第２の実施形態＞
第１の実施形態では、検体の一塩基多型が病的である可能性の程度を提示する場合の例について説明した。これに対し、第２の実施形態では、検体の二塩基以上の変異が病的である可能性の程度を提示する場合の例について説明する。

【0107】

図７は、第２の実施形態の制御部１１の構成を示す図である。制御部１１は、記憶部１２に格納されたプログラムを実行することにより、図７に例示する構成を機能的に実現する。図７に示す制御部は、図２の構成と比較すると、データベースフィルタ２３３及び機能予測フィルタ２３４を備えず、候補配列取得部３０１、保存位置取得部３０２、融合遺伝子フィルタ３０３、保存位置フィルタ３０４及び構造フィルタ３０５を備える点において異なる。図２と同様の機能ブロックについては同じ符号を付して説明を省略する。

【0108】

特定の組み合わせの２つの遺伝子が染色体の転座や逆位等に起因して融合することにより、がん細胞の増殖を引き起こすことが知られている。例えば、ＢＣＲ遺伝子とＡＢＬ遺伝子とが染色体の転座により融合したＢＣＲ－ＡＢＬ融合遺伝子は、白血病細胞を増殖させることが知られている。記憶部１２には、特定の組み合わせの２つの候補遺伝子が融合した融合遺伝子においてドライバー変異を引き起こすことが知られている複数の組み合わせの候補遺伝子がコード化された塩基配列が融合遺伝子ごとにそれぞれ記憶されている。例えば、ＢＣＲ遺伝子とＡＢＬ遺伝子とがコード化された塩基配列が記憶部１２に記憶されている。

【0109】

候補配列取得部３０１は、特定の組み合わせの候補遺伝子が融合した融合遺伝子（以下、第１融合遺伝子）においてドライバー変異候補となる２つの候補遺伝子の塩基配列を第１融合遺伝子ごとに取得する。第２の実施形態の例では、候補配列取得部３０１は、記憶部１２に記憶されている複数の第１融合遺伝子に含まれる２つの候補遺伝子のそれぞれの塩基配列を第１融合遺伝子ごとに記憶部１２から取得する。

【0110】

また、外部サーバ（不図示）が複数の第１融合遺伝子の候補遺伝子がコード化された塩基配列を記憶していてもよい。候補配列取得部３０１は、通信部１６を介して、外部サーバから第１融合遺伝子の２つの候補遺伝子がコード化された塩基配列を第１融合遺伝子ごとに取得してもよい。

【0111】

特定の候補遺伝子と他の遺伝子とが融合した融合遺伝子ががん細胞の増殖を引き起こすことがある。例えば、ＡＬＫ遺伝子が他の遺伝子と融合した融合遺伝子は、がん細胞の増殖を引き起こすことが知られている。記憶部１２には、他の遺伝子と融合した融合遺伝子（以下、第２融合遺伝子ともいう）においてドライバー変異候補となる複数の候補遺伝子の塩基配列が記憶されている。

【0112】

候補配列取得部３０１は、他の遺伝子と融合した第２融合遺伝子においてドライバー変異候補となる候補遺伝子の塩基配列を取得する。例えば、候補配列取得部３０１は、複数の第２融合遺伝子の候補遺伝子の塩基配列を記憶部１２から取得する。候補配列取得部３０１は、通信部１６を介して、複数の第２融合遺伝子の候補遺伝子の塩基配列を外部サーバから取得してもよい。

【0113】

保存位置取得部３０２は、異なる生物種のゲノムの間において保存された塩基配列である保存配列の位置を示す保存配列位置情報を取得する。例えば、保存位置取得部３０２は、保存配列位置情報を記憶部１２から取得する。保存位置取得部３０２は、通信部１６を介して、保存配列位置情報を外部サーバから取得してもよい。

【0114】

基本フィルタ２３１は、一塩基多型に特有の処理を実行しない点を除いて、図２と同様である。基本フィルタ２３１は、分析の対象とする変異状態が、良性であると判断できる場合に、良性変異であることを表すスコア（例えば「４」）を設定して、ランク推定部２４にその結果を出力する。また、基本フィルタ２３１は、分析の対象とする変異状態が、良性であると判断できなければ、良性変異でないことを表すスコア（例えば「３」）を設定し、次のフィルタとして設定されたフィルタに処理をパスする。

【0115】

基本フィルタ２３１は、設定受入部２２からがん化等を引き起こす既知の変異の塩基配列と変異状態に対応する変異した塩基配列との重複部分の長さの閾値を特定する情報と、データベースごとのパラメータ（良性か否かを判断する基準となる良性判断閾値などとして登録された値と比較される）の設定を受け入れて、当該設定に基づいて分析の対象とする変異状態が、良性であるか否かを判断する。

【0116】

具体的に基本フィルタ２３１は、がん化等を引き起こす既知の変異の塩基配列と変異状態に対応する変異した塩基配列との重複部分が予め定めた長さの閾値より短い重複部分である場合は良性変異であることを表すスコアを設定する。また基本フィルタ２３１は、そうでなくても、変異状態が表す、変異の位置する領域がイントロン領域であれば、良性変異であることを表すスコアを設定する。

【0117】

さらに基本フィルタ２３１は、上記２つの条件を満足しなくても、指定されたデータベースを検索し、検索によって変異状態が表す変異がデータベースに登録され、かつ、その変異である確率として登録された値が、当該データベースについて予め定められた良性判断閾値を超えている場合、良性変異であることを表すスコアを設定する。

【0118】

時系列フィルタ２３２は、分析の対象とする変異状態に対応するスコアから減算する値が第１の実施形態と異なることや、時系列フィルタ２３２による演算後のスコアの出力先が第１の実施形態と異なることを除いて、第１の実施形態と同様である。時系列フィルタ２３２は、分析の対象とする変異状態に対応する、時系列情報に含まれる変異状態の情報を参照して、異なるタイミングで抽出した時系列情報においても同じ変異があったか否かを判断する。

【0119】

時系列フィルタ２３２は、分析の対象とする変異状態と、時系列情報に含まれる対応する変異状態とを用い、同じ変異が存在する場合に、病的である可能性があるものとして分析の対象とする変異状態に対応するスコア（例えばスコアから第２所定量として「２」を引く）を決定してクオリティフィルタ２３５に処理をパスする。ここでの例では基本フィルタ２３１が処理をパスしているので、当初のスコアは「３」であり、ここで時系列フィルタ２３２が病的である可能性があるものとしたときには、このスコア「３」から第２所定量として「２」を引いてスコアを「１」と設定することとなる。第２所定量は、第１所定量より大きい値である。

【0120】

一方、時系列フィルタ２３２は、分析の対象とする変異状態と、時系列情報に含まれる対応する変異状態とを用い、同じ変異が存在しないときには、スコアをそのままに設定（ここでは当初のスコアが「３」であるので、そのまま「３」に設定）して、データベースフィルタ２３３へ処理をパスする。

【0121】

なお、時系列フィルタ２３２は、設定受入部２２から深度や、その他のシーケンスクオリティ、変異アレル頻度等に関する閾値の設定を受けていてもよい。例えば時系列情報に含まれる対応する変異状態に係る深度がここで設定された閾値（例えば「２０」）を超えない場合には、時系列フィルタ２３２は、同じ変異状態があったか否かを判断することなく、スコアをそのままに設定（ここでは当初のスコアが「３」であるので、そのまま「３」に設定）して、データベースフィルタ２３３へ処理をパスする。

【0122】

さらに、この時系列フィルタ２３２は、第１の実施形態と同様に、データ受入部２１が時系列情報を受け入れていない場合（変異塩基配列情報としては、分析の対象とする変異塩基配列情報だけを受け入れている場合）には、同じ変異状態があるか否かを判断することなく、スコアをそのままに設定（ここでは当初のスコアが「３」であるので、そのまま「３」に設定）して、データベースフィルタ２３３へ処理をパスしてもよい。

【0123】

また設定受入部２２から時系列フィルタ２３２を利用しない設定が入力されている場合、時系列フィルタ２３２は、同じ変異状態があるか否かを判断することなく、スコアをそのままに設定（ここでは当初のスコアが「３」であるので、そのまま「３」に設定）して、融合遺伝子フィルタ３０３へ処理をパスする。

【0124】

［融合遺伝子に類似する塩基配列を含むか否かの判定］
以下、変異塩基配列情報に含まれるいずれかの変異状態に対応する変異している塩基配列を変異塩基配列ともいう。融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が変異塩基配列に含まれているか否かを判定する。より詳しくは、融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した複数の第１融合遺伝子について、第１融合遺伝子の２つの候補遺伝子がコード化された２つの塩基配列と、変異塩基配列に含まれる少なくとも一部の塩基配列との類似度が２つとも閾値以上であるか否かを第１融合遺伝子ごとに判定する。類似度は、例えば２つの塩基配列のアライメントが一致する割合により表される。２つの塩基配列のアライメントが一致する割合が閾値以上である場合に、２つの塩基配列が類似すると判定される。

【0125】

一例としては、融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得したＢＣＲ遺伝子とＡＢＬ遺伝子とが融合したＢＣＲ－ＡＢＬ第１融合遺伝子においてＢＣＲ遺伝子がコード化された塩基配列と、変異塩基配列における対応する塩基配列との類似度を求める。次に、融合遺伝子フィルタ３０３は、ＢＣＲ－ＡＢＬ第１融合遺伝子においてＡＢＬ遺伝子がコード化された塩基配列と、変異塩基配列における対応する塩基配列との類似度を求める。

【0126】

融合遺伝子フィルタ３０３は、求めた２つの類似度が２つとも閾値以上であるか否かを判定する。閾値は、例えば、第１融合遺伝子がコード化されたタンパク質の活性と、変異塩基配列が示すタンパク質の活性とが同様であることが想定される値である。

【0127】

融合遺伝子フィルタ３０３は、求めた２つの類似度が２つとも閾値以上である場合に、第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が変異塩基配列に含まれていると判定する。

【0128】

一方、融合遺伝子フィルタ３０３は、求めた２つの類似度のうち、少なくとも一方の類似度が閾値未満である場合に、候補配列取得部３０１が取得した別の第１融合遺伝子について同様の判定を繰り返す。融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した全ての第１融合遺伝子について、求めた２つの類似度の少なくとも一方が閾値未満である場合に、どの第１融合遺伝子についても、第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が変異塩基配列に含まれていないと判定する。

【0129】

また、融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した第１融合遺伝子の２つの候補遺伝子の塩基配列と、変異塩基配列に含まれている融合遺伝子の２つの遺伝子の塩基配列との間の類似度がそれぞれ６５％以上１００％以下である場合に、第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が変異塩基配列に含まれていると判定してもよい。好ましくは、融合遺伝子フィルタ３０３は、第１融合遺伝子の２つの候補遺伝子の塩基配列と、変異塩基配列に含まれている融合遺伝子の２つの遺伝子の塩基配列との間の類似度がそれぞれ８０％以上１００％以下である場合に、第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が変異塩基配列に含まれていると判定してもよい。

【0130】

また、融合遺伝子フィルタ３０３は、複数の第１融合遺伝子の候補遺伝子の組み合わせを記憶している外部サーバへ、分析対象の変異状態に対応する変異塩基配列を送信してもよい。融合遺伝子フィルタ３０３は、外部サーバのデータベースに登録されている第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子の融合遺伝子が変異塩基配列に含まれているか否かを調べる。融合遺伝子フィルタ３０３は、外部サーバのデータベースに登録されている複数の第１融合遺伝子のうち、いずれかの第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子の融合遺伝子が変異塩基配列に含まれていることを示す通知を外部サーバから受信した場合に、第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が変異塩基配列に含まれていると判定してもよい。

【0131】

融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した第２融合遺伝子の候補遺伝子の塩基配列に類似する塩基配列の遺伝子と他の遺伝子とが融合した融合遺伝子が変異塩基配列に含まれているか否かを判定する。より詳しくは、融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した複数の第２融合遺伝子について、第２融合遺伝子の候補遺伝子の塩基配列と、変異塩基配列に含まれる融合遺伝子の一方の遺伝子の塩基配列との類似度を第２融合遺伝子ごとに求める。融合遺伝子フィルタ３０３は、求めた類似度が閾値以上であるか否かを判定する。閾値は、第２融合遺伝子がコード化されたタンパク質の活性と、変異塩基配列が示すタンパク質の活性とが同様であることが想定される値である。

【0132】

融合遺伝子フィルタ３０３は、求めた類似度が閾値以上である場合に、候補配列取得部３０１が取得した第２融合遺伝子の候補遺伝子と類似する遺伝子の融合遺伝子を変異塩基配列が含むと判定する。融合遺伝子フィルタ３０３は、求めた類似度が閾値未満である場合に、候補配列取得部３０１が取得した別の第２融合遺伝子の候補遺伝子について同様の判定を繰り返す。融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した全ての第２融合遺伝子について、求めた類似度が閾値未満である場合に、どの第２融合遺伝子の候補遺伝子と類似する遺伝子の融合遺伝子も変異塩基配列が含んでいないと判定する。

【0133】

また、融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した第２融合遺伝子の候補遺伝子の塩基配列と、変異塩基配列に含まれている融合遺伝子における一方の遺伝子の塩基配列との間の類似度が６５％以上１００％以下である場合に、第２融合遺伝子の候補遺伝子の塩基配列に類似する塩基配列の遺伝子と他の遺伝子とが融合した融合遺伝子が変異塩基配列に含まれていると判定してもよい。好ましくは、融合遺伝子フィルタ３０３は、第２融合遺伝子の候補遺伝子の塩基配列と、変異塩基配列に含まれている融合遺伝子における一方の遺伝子の塩基配列との間の類似度が８０％以上１００％以下である場合に、第２融合遺伝子の候補遺伝子の塩基配列に類似する塩基配列の遺伝子と他の遺伝子とが融合した融合遺伝子が変異塩基配列に含まれていると判定してもよい。

【0134】

また、融合遺伝子フィルタ３０３は、複数の第２融合遺伝子を記憶している外部サーバへ、変異塩基配列を送信してもよい。融合遺伝子フィルタ３０３は、外部サーバのデータベースに登録されている複数の第２融合遺伝子の候補遺伝子のいずれかと類似する遺伝子の融合遺伝子を変異塩基配列が含んでいるか否かを調べる。融合遺伝子フィルタ３０３は、登録されている複数の第２融合遺伝子の候補遺伝子のいずれかと類似する遺伝子の融合遺伝子を変異塩基配列が含んでいることを示す通知を外部サーバから受信した場合に、第２融合遺伝子の候補遺伝子と類似する遺伝子を変異塩基配列が含んでいると判定してもよい。

【0135】

融合遺伝子フィルタ３０３は、第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が変異塩基配列に含まれているか否かの判定結果によりスコアを決定する。例えば、融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した複数の第１融合遺伝子のいずれかについて、第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する２つの遺伝子が融合した融合遺伝子が変異塩基配列に含まれていると判定した場合に、病的である可能性があるものと、分析の対象とする変異状態に対応するスコアを決定して（例えばスコアから第２所定量として「２」を引く）クオリティフィルタ２３５に処理をパスする。

【0136】

このようにして、融合遺伝子フィルタ３０３は、ドライバー変異である可能性が比較的高いことが知られている第１融合遺伝子の２つの候補遺伝子の塩基配列を参照して、変異状態が病的である可能性の程度をスコアにより精度よく推定することができる。

【0137】

融合遺伝子フィルタ３０３は、第２融合遺伝子の候補遺伝子の塩基配列と類似する塩基配列の遺伝子が他の遺伝子と融合した融合遺伝子を変異塩基配列が含むか否かの判定結果によりスコアを決定する。例えば、融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した複数の第２融合遺伝子のいずれかの候補遺伝子と類似する遺伝子を変異塩基配列が含むと判定した場合に、病的な可能性があるものと分析の対象とする変異状態に対応するスコアを決定して（例えばスコアから第１所定量として「１」を引く）保存位置フィルタ３０４に処理をパスする。

【0138】

融合遺伝子フィルタ３０３は、候補配列取得部３０１が取得した第１融合遺伝子の２つの候補遺伝子とそれぞれ類似する候補遺伝子の融合遺伝子が変異塩基配列に含まれていないと判定した場合や、第２融合遺伝子の候補遺伝子と類似する遺伝子の融合遺伝子を変異塩基配列が含んでいないと判定した場合に、スコアをそのままに設定（ここでは当初のスコアが「３」であるので、そのまま「３」に設定）して、保存位置フィルタ３０４へ処理をパスする。

【0139】

融合遺伝子の２つの候補遺伝子の組み合わせの一方が記憶部１２に登録されていない場合であっても、特定の候補遺伝子を含む第２融合遺伝子についてはドライバー変異となる可能性があることが知られている。融合遺伝子フィルタ３０３は、第２融合遺伝子の候補遺伝子の塩基配列を参照することにより、変異状態が病的である可能性の程度をスコアにより精度よく提示することができる。

【0140】

［保存配列の位置が変異箇所に含まれるか否かの判定］
異なる生物種のゲノムの間において保存された保存配列は、細胞の生理活性に重要な役割を果たしていることが多い。このため、保存配列の位置に変異が生じている場合、変異状態が病的である可能性が比較的高くなる。保存位置フィルタ３０４は、異なる生物種のゲノムの間において保存された塩基配列である保存配列の位置が、変異状態の変異箇所に含まれるか否かによりスコアを決定する。より詳しくは、保存位置フィルタ３０４は、保存位置取得部３０２が取得した保存配列位置情報が示す保存配列の位置が、変異箇所に含まれるか否かを判定する。

【0141】

保存位置フィルタ３０４は、変異箇所に保存配列の位置が含まれると判定した場合に、病的である可能性があるものと分析の対象とする変異状態に対応するスコアを決定して（例えばスコアから第１所定量として「１」を引く）、構造フィルタ３０５に処理をパスする。一方、保存位置フィルタ３０４は、変異箇所に保存配列の位置が含まれていないと判定した場合に、スコアをそのままに設定して、構造フィルタ３０５へ処理をパスする。このようにして、保存位置フィルタ３０４は、保存配列の位置を示す情報を利用して、この変異箇所に対応する変異状態が病的である可能性の程度をスコアにより精度よく提示することができる。

【0142】

［構造変異の有無の判定］
染色体の転座や重要な遺伝子の欠失、複数の遺伝子に及ぶ変異等の構造変異が生じている場合、これらの構造変異が病的である可能性は比較的高いことが知られている。構造フィルタは、変異塩基配列情報が表す変異状態が染色体の転座等の構造変異であるか否かを判定する。

【0143】

構造フィルタ３０５は、変異塩基配列情報が表す変異状態が染色体の転座であるか否かを判定し、この判定結果によりスコアを決定する。構造フィルタ３０５は、変異塩基配列情報が示す変異状態に含まれる変異の内容や変異箇所を参照して、染色体の転座が生じているか否かを判定する。また、構造フィルタ３０５は、変異状態に対応する変異塩基配列を複数の塩基配列に分割し、分割した塩基配列ごとにゲノム上の位置を特定することにより、変異状態が染色体の転座であるか否かを判定してもよい。

【0144】

構造フィルタ３０５は、変異塩基配列情報が表す変異状態が複数の遺伝子に及ぶ変異であるか否かを判定し、この判定結果によりスコアを決定する。構造フィルタ３０５は、変異塩基配列情報が示すいずれかの変異状態に含まれる変異の内容や変異箇所を参照して、複数の遺伝子に及ぶ変異が生じているか否かを判定する。構造フィルタ３０５は、変異状態に対応する変異塩基配列を複数の塩基配列に分割し、分割した塩基配列ごとにゲノム上の位置を特定することにより、変異状態が複数の遺伝子に及ぶ変異であるか否かを判定してもよい。

【0145】

記憶部１２には、細胞のがん化等に関与する複数の登録遺伝子を示す情報が予め登録されている。登録遺伝子を示す情報は、例えば、登録遺伝子を識別するための識別情報や登録遺伝子の染色体上の位置を示す情報である。構造フィルタ３０５は、変異塩基配列情報が表す変異状態が登録遺伝子の欠失であるか否かを判定し、この判定結果によりスコアを決定してもよい。構造フィルタ３０５は、変異塩基配列情報が示すいずれかの変異状態に含まれる変異の内容や変異箇所を参照して、記憶部１２に登録されている複数の登録遺伝子のいずれかが欠失したか否かを判定する。

【0146】

記憶部１２には、細胞のがん化等に関与する遺伝子の発現を制御するエンハンサーの染色体上の位置情報が予め登録されている。構造フィルタ３０５は、転座、逆位、欠失等が生じていると判定した場合において、変異塩基配列情報が表す変異状態が記憶部１２に登録されているがん遺伝子が、記憶部１２に登録されているエンハンサーの近傍に位置する脱制御異常であるか否かを判定し、この判定結果によりスコアを決定してもよい。

【0147】

記憶部１２には、遺伝子領域のゲノムにおける向き（５’→３’, ３’→５’）の情報が予め登録されている。構造フィルタ３０５は、転座や欠失等により、変異塩基配列情報が表す変異状態が第１融合遺伝子や第２融合遺伝子等の融合遺伝子を形成すると判定した場合において融合遺伝子を形成する２つの遺伝子を第一候補遺伝子及び第二候補遺伝子とすると、第一候補遺伝子と第二候補遺伝子の向きがそれぞれ同一の方向であるか（例: 第一候補遺伝子５’→３’で第二候補遺伝子も ５’→３’方向、もしくは、第一候補遺伝子３’→５’で第二候補遺伝子 ３’→５’の組み合わせであるか）を判定し、機能的な融合遺伝子が形成されるかの有無を判定し、この判定結果によりスコアを決定してもよい。

【0148】

記憶部１２には、遺伝子領域のアミノ酸翻訳（コドン）やＲＮＡのスプライシングに関わる配列情報が予め登録されている。構造フィルタ３０５は、転座や欠失等により、変異塩基配列情報が表す変異状態が融合遺伝子を形成すると判定した場合において、上記項目の情報を元に機能的な融合遺伝子が形成されるかの有無を判定し、この判定結果によりスコアを決定してもよい。

【0149】

また、構造フィルタ３０５は、変異塩基配列を複数の塩基配列に分割し、分割した塩基配列ごとにゲノム上の位置を特定する。構造フィルタ３０５は、特定した塩基配列のゲノム上の位置と、記憶部１２に登録されている複数の登録遺伝子の位置とを比較することにより、いずれかの登録遺伝子の欠失が生じたか否かを判定してもよい。

【0150】

構造フィルタ３０５は、転座が生じていると判定した場合に、病的である可能性があるものとして分析の対象とする変異状態に対応するスコアを決定する。例えば、構造フィルタ３０５は、変異状態に対応するスコアから第１所定量として「１」を引く。一方、転座が生じていないと判定した場合に、分析の対象とする変異状態に対応するスコアをそのままとする。

【0151】

構造フィルタ３０５は、複数の遺伝子に及ぶ変異が生じていると判定した場合に、病的である可能性があるものとして分析の対象とする変異状態に対応するスコア（例えば変異状態に対応するスコアから第１所定量として「１」を引く）を決定する。一方、構造フィルタ３０５は、複数の遺伝子に及ぶ構造変異が生じていないと判定した場合に、変異状態に対応するスコアをそのままとする。

【0152】

構造フィルタ３０５は、記憶部１２に登録されている複数の登録遺伝子のいずれかが欠失していると判定した場合に、分析の対象とする変異状態に対応するスコアから第１所定量をさらに引いてクオリティフィルタ２３５に処理をパスする。一方、構造フィルタ３０５は、記憶部１２に登録されている複数の遺伝子がいずれも欠失していないと判定した場合に、分析の対象とする変異状態に対応するスコアをそのままとし、クオリティフィルタ２３５に処理をパスする。このようにして、構造フィルタ３０５は、染色体の転座や複数の遺伝子に及ぶ変異、細胞のがん化等に関与する遺伝子の欠失等の構造変異が生じているか否かを判定することにより、変異状態が病的である可能性の程度をスコアにより精度よく提示することができる。

【0153】

［変異の状態のスコア評価の処理手順］
図８は、変異塩基配列情報に含まれる複数の変異状態が病的である可能性の程度を分析装置１が提示する処理手順を示すフローチャートである。この処理手順は、例えば、変異状態が病的である可能性の程度を提示する処理を開始させるユーザの指示を操作受付部（不図示）が受け付けたときに開始する。

【0154】

まず、基本フィルタ２３１は、分析の対象とする変異塩基配列情報に含まれる、各変異状態を順次選択し（Ｓ２１）、選択した変異状態について、図５のＳ１２と同様にして、基本フィルタの表す条件により、病原性の可能性があるか否かを判定する（Ｓ２２）。

【0155】

そして基本フィルタ２３１は、選択した変異状態が病原性の可能性があると判定した場合（Ｓ２２のＹｅｓ）、この選択した変異状態に対応するスコアを良性変異でないことを示すスコアとして仮設定し、ステップＳ２３に移行する。

【0156】

時系列フィルタ２３２は、図５のＳ１３と同様にして、時系列フィルタの表す条件により病原性の可能性があるか否かを判定し（Ｓ２３）、判定結果によりスコアを設定する。時系列フィルタ２３２は、変異状態に病原性の可能性があると判定した場合に（Ｓ２３のＹｅｓ）、クオリティフィルタ２３５の表す条件により、クオリティが十分であるか否かを判定し（Ｓ２４）、判定結果によりスコアを決定する。クオリティフィルタ２３５は、決定されたスコアを選択された変異状態に対応するランク値とし、選択された変異状態を特定する情報にこのランク値を関連付けて記録する（Ｓ２５）。

【0157】

ランク推定部２４は、変異塩基配列情報に含まれる全ての変異状態を特定する情報にランク値を関連付けて記録するまで、Ｓ２２以降の処理を繰り返す。出力部２５は、変異塩基配列情報に含まれる変異状態のそれぞれを特定する情報と、それぞれに関連付けられたランク値とを集積して出力し（Ｓ２６）、処理を終了する。

【0158】

融合遺伝子フィルタ３０３は、ステップＳ２３において、病原性の可能性がないと時系列フィルタ２３２が判定した場合（Ｓ２３のＮｏ）、第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する遺伝子を融合させた融合遺伝子を、選択された変異状態に対応する変異塩基配列が含むか否かを判定する（Ｓ２７）。融合遺伝子フィルタ３０３は、この判定結果によりスコアを決定する。

【0159】

融合遺伝子フィルタ３０３は、第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する遺伝子を融合させた融合遺伝子を変異塩基配列が含まないと判定した場合に（Ｓ２７のＮｏ）、第２融合遺伝子の候補遺伝子に類似する遺伝子の融合遺伝子を変異塩基配列が含むか否かを判定する（Ｓ２８）。融合遺伝子フィルタ３０３は、この判定結果によりスコアを決定する。保存位置フィルタ３０４は、選択された変異状態に対応する変異箇所に保存配列の位置が含まれるか否かを判定し（Ｓ２９）、この判定結果によりスコアを決定する。

【0160】

構造フィルタ３０５は、選択された変異状態が各種の構造変異を含むか否かを判定する（Ｓ３０）。例えば、構造フィルタ３０５は、選択された変異状態が染色体の転座であるか否かを判定し、判定結果によりスコアを決定する。構造フィルタ３０５は、選択された変異状態が複数の遺伝子に及ぶ変異であるか否かを判定し、判定結果によりスコアを決定する。構造フィルタ３０５は、選択された変異状態が記憶部１２に登録された複数の登録遺伝子のいずれかを欠失する変異であるか否かを判定し、判定結果によりスコアを決定する。構造フィルタ３０５は、Ｓ３０においてそれぞれの判定を実行した後、Ｓ２４の処理に移る。

【0161】

融合遺伝子フィルタ３０３は、Ｓ２７の判定において第１融合遺伝子の２つの候補遺伝子にそれぞれ類似する遺伝子を融合させた融合遺伝子を変異塩基配列が含むと判定した場合に（Ｓ２７のＹｅｓ）、Ｓ２４の判定に移る。

【0162】

以上、本発明を実施の形態を用いて説明したが、本発明の技術的範囲は上記実施の形態に記載の範囲には限定されず、その要旨の範囲内で種々の変形及び変更が可能である。例えば、装置の全部又は一部は、任意の単位で機能的又は物理的に分散・統合して構成することができる。また、複数の実施の形態の任意の組み合わせによって生じる新たな実施の形態も、本発明の実施の形態に含まれる。組み合わせによって生じる新たな実施の形態の効果は、もとの実施の形態の効果を併せ持つ。

【符号の説明】

【0163】

１ 分析装置
１１ 制御部
１２ 記憶部
１３ 操作部
１４ 表示部
１５ 入出力部
１６ 通信部
２１ データ受入部
２２ 設定受入部
２３ フィルタ処理部
２４ ランク推定部
２５ 出力部
２３１ 基本フィルタ
２３２ 時系列フィルタ
２３３ データベースフィルタ
２３４ 機能予測フィルタ
２３５ クオリティフィルタ
３０１ 候補配列取得部
３０２ 保存位置取得部
３０３ 融合遺伝子フィルタ
３０４ 保存位置フィルタ
３０５ 構造フィルタ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】