知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-12-09
(45)【発行日】2024-12-17
(54)【発明の名称】遺伝子情報処理システム、遺伝子情報処理方法、及び記録媒体
(51)【国際特許分類】
G16B 25/00 20190101AFI20241210BHJP
【ＦＩ】
G16B25/00
【請求項の数】
9
(21)【出願番号】P 2023506774
(86)(22)【出願日】2022-01-06
(86)【国際出願番号】 JP2022000221
(87)【国際公開番号】W WO2022196041
(87)【国際公開日】2022-09-22
【審査請求日】2023-08-15
(31)【優先権主張番号】P 2021041423
(32)【優先日】2021-03-15
(33)【優先権主張国・地域又は機関】JP
(73)【特許権者】
【識別番号】000004237
【氏名又は名称】日本電気株式会社
(74)【代理人】
【識別番号】100104765
【弁理士】
【氏名又は名称】江上 達夫
(74)【代理人】
【識別番号】100107331
【弁理士】
【氏名又は名称】中村 聡延
(74)【代理人】
【識別番号】100131015
【弁理士】
【氏名又は名称】三輪 浩誉
(72)【発明者】
【氏名】麻生川 稔
【審査官】松田 岳士
(56)【参考文献】
【文献】特開２００７－２５９８４７（ＪＰ，Ａ）
【文献】特開２００６－１６３７２０（ＪＰ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ１６Ｂ ５／００－９９／００
Ｃ１２Ｍ １／００－ ３／１０
Ｃ１２Ｑ １／００－ ３／００
(57)【特許請求の範囲】
【請求項1】
ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報を取得する第１取得手段と、前記ＰＣＲの１サイクルごとの前記マイクロサテライトの増幅率を示す第１の割合と、前記ＰＣＲの１サイクルごとに前記マイクロサテライトから発生するスタッターの割合を示す第２の割合と、前記ＰＣＲのサイクル数と、を取得する第２取得手段と、
前記ＰＣＲ後の遺伝子における前記マイクロサテライトの分布情報と、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、に基づいて、前記ＰＣＲ前の前記遺伝子における前記マイクロサテライトの分布を算出する算出手段と
を備え、
前記算出手段は、前記ＰＣＲによる増幅をモデル化した行列であって、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、を成分として含む増幅行列を用いて、前記ＰＣＲ前の前記遺伝子における前記マイクロサテライトの分布を算出する
ことを特徴とする遺伝子情報処理システム。
【請求項2】
前記算出手段は、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、を用いて、繰り返しユニットが１つ少ない又は多いスタッターのアンプリコン量を推定することを特徴とする請求項１に記載の遺伝子情報処理システム。
【請求項3】
前記算出手段は、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、を用いて、繰り返しユニットが２つ少ない又は多いスタッターのアンプリコン量を推定することを特徴とする請求項１又は２に記載の遺伝子情報処理システム。
【請求項4】
前記算出手段は、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、を用いて、前記ＰＣＲで増幅されたマイクロサテライトのアンプリコン量に対する前記スタッターのアンプリコン量の比率であるスタッター比を推定することを特徴とする請求項１から３のいずれか一項に記載の遺伝子情報処理システム。
【請求項5】
前記第２の割合は、繰り返しユニットが少ないマイナススタッターが発生する割合と、繰り返しユニットが多いプラススタッターの割合と、の両方を含むことを特徴とする請求項１から４のいずれか一項に記載の遺伝子情報処理システム。
【請求項6】
前記第２取得手段は、前記ＰＣＲ後の前記遺伝子における前記マイクロサテライトの分布情報から、前記第１の割合、前記第２の割合、及び前記ＰＣＲのサイクル数の少なくとも１つを推定して取得することを特徴とする請求項１から５のいずれか一項に記載の遺伝子情報処理システム。
【請求項7】
前記第２取得手段は、前記第１の割合、前記第２の割合、及び前記ＰＣＲのサイクル数の少なくとも１つを、山登り法又はデータフィッティングにより推定することを特徴とする請求項６に記載の遺伝子情報処理システム。
【請求項8】
少なくとも１つのコンピュータが、
ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報を取得し、
前記ＰＣＲの１サイクルごとの前記マイクロサテライトの増幅率を示す第１の割合と、前記ＰＣＲの１サイクルごとに前記マイクロサテライトから発生するスタッターの割合を示す第２の割合と、前記ＰＣＲのサイクル数と、を取得し、
前記ＰＣＲ後の遺伝子における前記マイクロサテライトの分布情報と、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、に基づいて、前記ＰＣＲ前の前記遺伝子における前記マイクロサテライトの分布を算出し、
前記マイクロサテライトの分を算出する際には、前記ＰＣＲによる増幅をモデル化した行列であって、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、を成分として含む増幅行列を用いて、前記ＰＣＲ前の前記遺伝子における前記マイクロサテライトの分布を算出する
ことを特徴とする遺伝子情報処理方法。
【請求項9】
ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報を取得し、前記ＰＣＲの１サイクルごとの前記マイクロサテライトの増幅率を示す第１の割合と、前記ＰＣＲの１サイクルごとに前記マイクロサテライトから発生するスタッターの割合を示す第２の割合と、前記ＰＣＲのサイクル数と、を取得し、
前記ＰＣＲ後の遺伝子における前記マイクロサテライトの分布情報と、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、に基づいて、前記ＰＣＲ前の前記遺伝子における前記マイクロサテライトの分布を算出し、
前記マイクロサテライトの分を算出する際には、前記ＰＣＲによる増幅をモデル化した行列であって、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、を成分として含む増幅行列を用いて、前記ＰＣＲ前の前記遺伝子における前記マイクロサテライトの分布を算出する
ようにコンピュータを動作させることを特徴とするコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
この開示は、遺伝子情報に関する処理を実行する遺伝子情報処理システム、遺伝子情報処理方法、及び記録媒体の技術分野に関する。
【背景技術】
【0002】
この種のシステムとして、ＰＣＲ（Ｐｏｌｙｍｅｒａｓｅ Ｃｈａｉｎ Ｒｅａｃｔｉｏｎ）において発生するスタッターに関する情報を扱うものが知られている。例えば特許文献１では、ＰＣＲの増幅後に現れるスタッターの相対的な高さを見積もることが記載されている。特許文献２では、ＰＣＲによって得られるピークが真のピークであるかノイズピークであるかを問わずに、ピーク判定を行う技術が開示されている。特許文献３では、パーセントスタッターを算出し、任意のサンプルの値をまとめた平均パーセントスタッターを算出することが開示されている。特許文献４では、ピーク位置、ピーク高、スタッター率を含むモデルパラメータを最適化することが開示されている。
【先行技術文献】
【特許文献】
【0003】
【文献】特開２００６－１６３７２０号公報
【文献】特開２００６－０７９３３４号公報
【文献】特表２００４－５３３２４１号公報
【文献】特表２００４－５１６４５５号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
この開示は、例えば上記各引用文献に鑑みてなされたものであり、遺伝子に関する処理を適切に実行することが可能な遺伝子情報処理システム、遺伝子情報処理方法、及び記録媒体を提供することを目的とする。
【課題を解決するための手段】
【0005】
この開示の遺伝子情報処理システムの一の態様は、ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報を取得する第１取得手段と、前記ＰＣＲの１サイクルごとの前記マイクロサテライトの増幅率を示す第１の割合と、前記ＰＣＲの１サイクルごとに前記マイクロサテライトから発生するスタッターの割合を示す第２の割合と、前記ＰＣＲのサイクル数と、を取得する第２取得手段と、前記ＰＣＲ後の遺伝子における前記マイクロサテライトの分布情報と、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、に基づいて、前記ＰＣＲ前の前記遺伝子における前記マイクロサテライトの分布を算出する算出手段とを備える。
【0006】
この開示の遺伝子情報処理方法の一の態様は、ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報を取得し、前記ＰＣＲの１サイクルごとの前記マイクロサテライトの増幅率を示す第１の割合と、前記ＰＣＲの１サイクルごとに前記マイクロサテライトから発生するスタッターの割合を示す第２の割合と、前記ＰＣＲのサイクル数と、を取得し、前記ＰＣＲ後の遺伝子における前記マイクロサテライトの分布情報と、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、に基づいて、前記ＰＣＲ前の前記遺伝子における前記マイクロサテライトの分布を算出する。
【0007】
この開示の記録媒体の一の態様は、ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報を取得し、前記ＰＣＲの１サイクルごとの前記マイクロサテライトの増幅率を示す第１の割合と、前記ＰＣＲの１サイクルごとに前記マイクロサテライトから発生するスタッターの割合を示す第２の割合と、前記ＰＣＲのサイクル数と、を取得し、前記ＰＣＲ後の遺伝子における前記マイクロサテライトの分布情報と、前記第１の割合と、前記第２の割合と、前記ＰＣＲのサイクル数と、に基づいて、前記ＰＣＲ前の前記遺伝子における前記マイクロサテライトの分布を算出するようにコンピュータを動作させるコンピュータプログラムが記録されている。
【図面の簡単な説明】
【0008】
【図1】第１実施形態に係る遺伝子情報処理システムのハードウェア構成を示すブロック図である。
【図2】第１実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。
【図3】ＰＣＲにおいて発生するスタッターの一例を示すグラフである。
【図4】第１実施形態に係る遺伝子情報処理システムの動作の流れを示すフローチャートである。
【図5】第２実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。
【図6】－１スタッターの一例を示すグラフである。
【図7】スタッターのアンプリコン量を推定する方法を示す概念図である。
【図8】第２実施形態に係る遺伝子情報処理システムの動作の流れを示すフローチャートである。
【図9】－２スタッター及び＋２スタッターの一例を示すグラフである。
【図10】第４実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。
【図11】第４実施形態に係る遺伝子情報処理システムの動作の流れを示すフローチャートである。
【図12】マイナススタッターが発生する割合及びプラススタッターが発生する割合の両方を考慮する例を示す概念図である。
【図13】第６実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。
【図14】ＰＣＲによる増幅をモデル化した増幅行列の一例を示す行列式である。
【図15】第６実施形態に係る遺伝子情報処理システムの動作の流れを示すフローチャートである。
【図16】第７実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。
【図17】第７実施形態に係る遺伝子情報処理システムの動作の流れを示すフローチャートである。
【図18】第８実施形態に係る遺伝子情報処理システムにおいて算出される遺伝子情報の推算量と実際量とを示すグラフである。
【図19】繰り返し回数と－１スタッター比との関係を示すグラフである。
【図20】データフィッティングによりスタッターが発生する割合を推定する方法を示すグラフである。
【発明を実施するための形態】
【0009】
以下、図面を参照しながら、遺伝子情報処理システム、遺伝子情報処理方法、コンピュータプログラム、及び記録媒体の実施形態について説明する。
【0010】
【0011】
（ハードウェア構成）
まず、図１を参照しながら、第１実施形態に係る遺伝子情報処理システム１０のハードウェア構成について説明する。図１は、第１実施形態に係る遺伝子情報処理システムのハードウェア構成を示すブロック図である。
まず、図１を参照しながら、第１実施形態に係る遺伝子情報処理システム１０のハードウェア構成について説明する。図１は、第１実施形態に係る遺伝子情報処理システムのハードウェア構成を示すブロック図である。
【0012】
図１に示すように、第１実施形態に係る遺伝子情報処理システム１０は、プロセッサ１１と、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）１２と、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）１３と、記憶装置１４とを備えている。遺伝子情報処理システム１０は更に、入力装置１５と、出力装置１６とを備えていてもよい。プロセッサ１１と、ＲＡＭ１２と、ＲＯＭ１３と、記憶装置１４と、入力装置１５と、出力装置１６と、カメラ２０とは、データバス１７を介して接続されている。
【0013】
プロセッサ１１は、コンピュータプログラムを読み込む。例えば、プロセッサ１１は、ＲＡＭ１２、ＲＯＭ１３及び記憶装置１４のうちの少なくとも一つが記憶しているコンピュータプログラムを読み込むように構成されている。或いは、プロセッサ１１は、コンピュータで読み取り可能な記録媒体が記憶しているコンピュータプログラムを、図示しない記録媒体読み取り装置を用いて読み込んでもよい。プロセッサ１１は、ネットワークインタフェースを介して、遺伝子情報処理システム１０の外部に配置される不図示の装置からコンピュータプログラムを取得してもよい（つまり、読み込んでもよい）。プロセッサ１１は、読み込んだコンピュータプログラムを実行することで、ＲＡＭ１２、記憶装置１４、入力装置１５及び出力装置１６を制御する。本実施形態では特に、プロセッサ１１が読み込んだコンピュータプログラムを実行すると、プロセッサ１１内には、遺伝子情報に関する処理を実行するための機能ブロックが実現される。また、プロセッサ１１として、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＦＰＧＡ（ｆｉｅｌｄ－ｐｒｏｇｒａｍｍａｂｌｅ ｇａｔｅ ａｒｒａｙ）、ＤＳＰ（Ｄｅｍａｎｄ－Ｓｉｄｅ Ｐｌａｔｆｏｒｍ）、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）のうち一つを用いてもよいし、複数を並列で用いてもよい。
【0014】
ＲＡＭ１２は、プロセッサ１１が実行するコンピュータプログラムを一時的に記憶する。ＲＡＭ１２は、プロセッサ１１がコンピュータプログラムを実行している際にプロセッサ１１が一時的に使用するデータを一時的に記憶する。ＲＡＭ１２は、例えば、Ｄ－ＲＡＭ（Ｄｙｎａｍｉｃ ＲＡＭ）であってもよい。
【0015】
ＲＯＭ１３は、プロセッサ１１が実行するコンピュータプログラムを記憶する。ＲＯＭ１３は、その他に固定的なデータを記憶していてもよい。ＲＯＭ１３は、例えば、Ｐ－ＲＯＭ（Ｐｒｏｇｒａｍｍａｂｌｅ ＲＯＭ）であってもよい。
【0016】
記憶装置１４は、遺伝子情報処理システム１０が長期的に保存するデータを記憶する。記憶装置１４は、プロセッサ１１の一時記憶装置として動作してもよい。記憶装置１４は、例えば、ハードディスク装置、光磁気ディスク装置、ＳＳＤ（Ｓｏｌｉｄ Ｓｔａｔｅ Ｄｒｉｖｅ）及びディスクアレイ装置のうちの少なくとも一つを含んでいてもよい。
【0017】
入力装置１５は、遺伝子情報処理システム１０のユーザからの入力指示を受け取る装置である。入力装置１５は、例えば、キーボード、マウス及びタッチパネルのうちの少なくとも一つを含んでいてもよい。
【0018】
出力装置１６は、遺伝子情報処理システム１０に関する情報を外部に対して出力する装置である。例えば、出力装置１６は、遺伝子情報処理システム１０に関する情報を表示可能な表示装置（例えば、ディスプレイ）であってもよい。
【0019】
（機能的構成）
次に、図２を参照しながら、第１実施形態に係る遺伝子情報処理システム１０の機能的構成について説明する。図２は、第１実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。
次に、図２を参照しながら、第１実施形態に係る遺伝子情報処理システム１０の機能的構成について説明する。図２は、第１実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。
【0020】
図２に示すように、第１実施形態に係る遺伝子情報処理システム１０は、その機能を実現するための処理ブロックとして、第１取得部１１０と、第２取得部１２０と、算出部１３０とを備えて構成されている。第１取得部１１０、第２取得部１２０、及び算出部１３０の各々は、例えば上述したプロセッサ１１（図１参照）において実現されればよい。
【0021】
第１取得部１１０は、ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報を取得可能に構成されている。第１取得部１１０は、例えば事前に行われた実験（即ち、ＰＣＲを用いた増幅処理）の結果として、マイクロサテライトの分布情報を取得すればよい。マイクロサテライトの分布情報は、一の生体（即ち、単数の生体）のマイクロサテライトの分布情報であってもよいし、複数の生体のマイクロサテライトの分布情報（例えば、混合サンプルから得られた分布情報）であってもよい。なお、単数の生体由来でも、１つの遺伝子座のマイクロサテライトは２種であることがある（正確には、１種又は２種である）。マイクロサテライトの分布情報は、このような単数の生体由来の２種のマイクロサテライトの分布情報であってもよい。
【0022】
第２取得部１２０は、ＰＣＲの１サイクルごとのマイクロサテライトの増幅率を示す第１の割合と、ＰＣＲの１サイクルごとにマイクロサテライトから発生するスタッターの割合を示す第２の割合と、ＰＣＲのサイクル数と、を取得可能に構成されている。第２取得部１２０は、第１の割合、第２の割合、及びＰＣＲのサイクル数の少なくとも１つを、例えばユーザ等が入力する既知の値として取得してもよい。或いは、第２取得部１２０は、第２取得部１２０は、第１の割合、第２の割合、及びＰＣＲのサイクル数の少なくとも１つを、他のパラメータから推定して取得してもよい。なお、各パラメータを推定する方法については、後述する他の実施形態で詳しく説明する。
【0023】
算出部１３０は、第１取得部１１０で取得されたＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報と、第２取得部１２０で取得された第１の割合、第２の割合、及びＰＣＲのサイクル数とに基づいて、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を算出可能に構成されている。即ち、算出部１３０は、ＰＣＲの対象となった遺伝子の元々の（即ち、増幅前の）マイクロサテライトの分布を算出可能に構成さている。算出部１３０による具体的な算出の方法については、後述する他の実施形態で詳しく説明する。
【0024】
【0025】
図３に示すように、遺伝子情報をＰＣＲによって増幅すると、本物の遺伝子情報（以下、適宜「真のアリル」と称する）に加えて、スタッターと呼ばれるノイズが検出される。スタッターのピークは、真のアリルのピーク位置に対して、繰り返しユニット数が少ない又は多い位置で検出される。本実施形態に係る遺伝子情報処理システム１０は、このスタッターが発生する割合を第２の割合として取得することで、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を算出する。言い換えれば、ＰＣＲによる増幅でスタッターが生じている場合でも、真のアリルに対応する増幅前の分布情報を得ることができる。
【0026】
（動作の流れ）
次に、図４を参照しながら、第１実施形態に係る遺伝子情報処理システム１０の動作の流れについて説明する。図４は、第１実施形態に係る遺伝子情報処理システムの動作の流れを示すフローチャートである。
次に、図４を参照しながら、第１実施形態に係る遺伝子情報処理システム１０の動作の流れについて説明する。図４は、第１実施形態に係る遺伝子情報処理システムの動作の流れを示すフローチャートである。
【0027】
図４に示すように、第１実施形態に係る遺伝子情報処理システム１０が動作する際には、まず第１取得部１１０が、ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報を取得する（ステップＳ１０１）。
【0028】
続いて、第２取得部１２０が、第１の割合を取得し（ステップＳ１０２）、第２の割合を取得し（ステップＳ１０３）、ＰＣＲのサイクル数を取得する（ステップＳ１０４）。なお、第１の割合、第２の割合、ＰＣＲのサイクル数を取得する順序は、上記の順に限定されるものではなく、互いに前後して取得されてもよいし、同時に取得されてもよい。
【0029】
続いて、算出部１３０が、第１取得部１１０で取得されたＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報と、第２取得部１２０で取得された第１の割合、第２の割合、及びＰＣＲのサイクル数とに基づいて、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を算出する（ステップＳ１０５）。算出部１３０は、算出したＰＣＲ前の遺伝子におけるマイクロサテライトの分布を出力する機能を有していてもよい。算出部１３０は、例えばディスプレイ等に算出したマイクロサテライトの分布を表示するようにしてもよい。
【0030】
（技術的効果）
次に、第１実施形態に係る遺伝子情報処理システム１０によって得られる技術的効果について説明する。
次に、第１実施形態に係る遺伝子情報処理システム１０によって得られる技術的効果について説明する。
【0031】
図１から図４で説明したように、第１実施形態に係る遺伝子情報処理システム１０では、ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報と、第１の割合と、第２の割合と、ＰＣＲのサイクル数とに基づいて、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布が算出される。このようにすれば、ＰＣＲにおいて発生するスタッターの存在を考慮して、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を適切に算出することが可能である。
【0032】
＜第２実施形態＞
第２実施形態に係る遺伝子情報処理システム１０について、図５から図８を参照して説明する。なお、第２実施形態は、上述した第１実施形態と一部の動作が異なるのみであり、例えばハードウェア構成等については、第１実施形態（図１参照）と同一であってよい。このため、以下では、第１実施形態と重複する部分については適宜説明を省略するものとする。
第２実施形態に係る遺伝子情報処理システム１０について、図５から図８を参照して説明する。なお、第２実施形態は、上述した第１実施形態と一部の動作が異なるのみであり、例えばハードウェア構成等については、第１実施形態（図１参照）と同一であってよい。このため、以下では、第１実施形態と重複する部分については適宜説明を省略するものとする。
【0033】
（機能的構成）
まず、図５を参照しながら、第２実施形態に係る遺伝子情報処理システム１０の機能的構成について説明する。図５は、第２実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。なお、図５では、図２で示した構成要素と同様の要素に同一の符号を付している。
まず、図５を参照しながら、第２実施形態に係る遺伝子情報処理システム１０の機能的構成について説明する。図５は、第２実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。なお、図５では、図２で示した構成要素と同様の要素に同一の符号を付している。
【0034】
図５に示すように、第２実施形態に係る遺伝子情報処理システム１０は、その機能を実現するための処理ブロックとして、第１取得部１１０と、第２取得部１２０と、算出部１３０とを備えて構成されている。そして特に、第２実施形態に係る算出部１３０は、アンプリコン量推定部１３１を備えている。
【0035】
アンプリコン量推定部１３１は、第１の割合と、第２の割合と、ＰＣＲ法のサイクル数と、を用いて、繰り返しユニットが１つ少ないスタッター（以下、適宜「－１スタッター」と称する）、又は繰り返しユニットが１つ多いスタッター（以下、適宜「＋１スタッター」と称する）のアンプリコン量を推定可能に構成されている。そして、第２実施形態に係る算出部１３０は、アンプリコン量推定部１３１で推定された－１スタッター又は＋１スタッターのアンプリコン量を用いて、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を算出可能に構成されている。
【0036】
【0037】
図６に示すように、－１スタッターのピークは、真のアリルのピークの位置に対して、繰り返しユニット数が１個分少ない位置で検出される。－１スタッターのピークは、その他のスタッターと比較して発生割合が大きく、ピークの大きさも比較的大きなものとなる。よって、－１スタッターのアンプリコン量を推定することができれば、真のアリルのマイクロサテライトの分布を精度よく算出できるようになる。なお、－１スタッターのアンプリコン量を推定するための具体的な手法については、以下で詳しく説明する。
【0038】
【0039】
図７に示すように、遺伝子情報のうち「５回繰り返し」のものをＰＣＲで増幅する場合を考える。なお、ここでは説明の便宜上、－１スタッターのみが発生する場合を考えるものとする。また、ＰＣＲの１サイクルごとのマイクロサテライトの増幅率（即ち、第１の割合）は“ａ”、ＰＣＲの１サイクルごとにマイクロサテライトから発生するスタッターの割合（即ち、第２の割合）は“ｂ”、ＰＣＲのサイクル数は“ｎ”とする。
【0040】
図７を見ると分かるように、ＰＣＲを３サイクルした場合、「５回繰り返し」のアンプリコン量は、ａ３となる。また、－１スタッター（即ち、「４回繰り返し」）のアンプリコン量は、３×ａ２×ｂとなる。そして、ＰＣＲをｎサイクルした場合、「５回繰り返し」のアンプリコン量は、ａｎとなる。また、－１スタッター（即ち、「４回繰り返し」）のアンプリコン量は、ｎＣ１×ａｎ－１×ｂとなる（Ｃはコンビネーション）。
【0041】
以上のように、第１の割合ａ、第２の割合ｂ、ＰＣＲのサイクル数ｎを用いれば、－１スタッターのアンプリコン量を推定することが可能である。なお、ここでの説明は省略したが、＋１スタッター（即ち、真のアリルのピークの位置に対して、繰り返しユニット数が１個分多い位置で検出されるスタッター）のアンプリコン量を推定してもよい。＋１スタッターのアンプリコン量を推定する場合には、上記の第１の割合ａ、第２の割合ｂに加えて、プラススタッターが発生する第３の割合ｃを考慮すればよい。第３の割合ｃを用いる構成については、後述する他の実施形態（第５実施形態）で詳しく説明する。
【0042】
（動作の流れ）
次に、図８を参照しながら、第２実施形態に係る遺伝子情報処理システム１０の動作の流れについて説明する。図８は、第２実施形態に係る遺伝子情報処理システムの動作の流れを示すフローチャートである。なお、図８では、図４で示した処理と同様の処理に同一の符号を付している。
次に、図８を参照しながら、第２実施形態に係る遺伝子情報処理システム１０の動作の流れについて説明する。図８は、第２実施形態に係る遺伝子情報処理システムの動作の流れを示すフローチャートである。なお、図８では、図４で示した処理と同様の処理に同一の符号を付している。
【0043】
図８に示すように、第２実施形態に係る遺伝子情報処理システム１０が動作する際には、まず第１取得部１１０が、ＰＣＲ後の遺伝子におけるマイクロサテライトの分布情報を取得する（ステップＳ１０１）。
【0044】
続いて、第２取得部１２０が、第１の割合を取得し（ステップＳ１０２）、第２の割合を取得し（ステップＳ１０３）、ＰＣＲのサイクル数を取得する（ステップＳ１０４）。
【0045】
続いて、アンプリコン量推定部１３１が、第１の割合、第２の割合、及びＰＣＲのサイクル数を用いて、－１スタッターのアンプリコン量又は＋１スタッターのアンプリコン量を推定する（ステップＳ２０１）。そして、算出部１３０は、アンプリコン量推定部１３１が推定したアンプリコン量を用いて、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を算出する（ステップＳ２０２）。
【0046】
（技術的効果）
次に、第２実施形態に係る遺伝子情報処理システム１０によって得られる技術的効果について説明する。
次に、第２実施形態に係る遺伝子情報処理システム１０によって得られる技術的効果について説明する。
【0047】
図５から図８で説明したように、第２実施形態に係る遺伝子情報処理システム１０では、－１スタッター又は＋１スタッターのアンプリコン量を推定し、推定したアンプリコン量を用いて、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布が算出される。このようにすれば、ＰＣＲにおいて発生する－１スタッター又は＋１スタッターの存在を考慮して、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を適切に算出することが可能である。
【0048】
＜第３実施形態＞
第３実施形態に係る遺伝子情報処理システム１０について、図９を参照して説明する。なお、第３実施形態は、上述した第２実施形態と比べて、扱うスタッターの種類が異なるのみであり、システムの構成や全体的な動作の流れ等については、第２実施形態と同一であってよい。このため、以下では、すでに説明した部分と重複する部分については適宜説明を省略するものとする。
第３実施形態に係る遺伝子情報処理システム１０について、図９を参照して説明する。なお、第３実施形態は、上述した第２実施形態と比べて、扱うスタッターの種類が異なるのみであり、システムの構成や全体的な動作の流れ等については、第２実施形態と同一であってよい。このため、以下では、すでに説明した部分と重複する部分については適宜説明を省略するものとする。
【0049】
（アンプリコン量推定部の構成）
第３実施形態に係る遺伝子情報処理システム１０では、アンプリコン量推定部１３１が、繰り返しユニットが２個少ないスタッター（以下、適宜「－２スタッター」と称する）、又は繰り返しユニットが２個少ないスタッター（以下、適宜「＋２スタッター」と称する）のアンプリコン量を推定する。アンプリコン量推定部１３１は、典型的には、すでに第２実施形態で説明した－１スタッター又は＋１スタッターのアンプリコン量と共に、－２スタッター又は＋２スタッターのアンプリコン量を推定する。ただし、アンプリコン量推定部１３１は、－１スタッター又は＋１スタッターのアンプリコン量を推定せずに、－２スタッター又は＋２スタッターのアンプリコン量のみを推定してもよい。
第３実施形態に係る遺伝子情報処理システム１０では、アンプリコン量推定部１３１が、繰り返しユニットが２個少ないスタッター（以下、適宜「－２スタッター」と称する）、又は繰り返しユニットが２個少ないスタッター（以下、適宜「＋２スタッター」と称する）のアンプリコン量を推定する。アンプリコン量推定部１３１は、典型的には、すでに第２実施形態で説明した－１スタッター又は＋１スタッターのアンプリコン量と共に、－２スタッター又は＋２スタッターのアンプリコン量を推定する。ただし、アンプリコン量推定部１３１は、－１スタッター又は＋１スタッターのアンプリコン量を推定せずに、－２スタッター又は＋２スタッターのアンプリコン量のみを推定してもよい。
【0050】
（－２スタッター及び＋２スタッター）
次に、図９を参照しながら、第３実施形態に係るアンプリコン量推定部１３１がアンプリコン量を推定する－２スタッター及び＋２スタッターについて具体的に説明する。図９は、－２スタッター及び＋２スタッターの一例を示すグラフである。
次に、図９を参照しながら、第３実施形態に係るアンプリコン量推定部１３１がアンプリコン量を推定する－２スタッター及び＋２スタッターについて具体的に説明する。図９は、－２スタッター及び＋２スタッターの一例を示すグラフである。
【0051】
図９に示すように、－２スタッターのピークは、真のアリルのピークの位置に対して、繰り返しユニット数が２個分少ない位置で検出される。一方、＋２スタッターのピークは、真のアリルのピークの位置に対して、繰り返しユニット数が２個分多い位置で検出される。－２スタッターは、－１スタッターが発生する現象が２回発生することで生ずる。同様に、＋２スタッターは、＋１スタッターが発生する現象が２回発生することで生ずる。このため－２スタッター及び＋２スタッターは、－１スタッター及び＋１スタッターと比べて、その発生割合が小さく、ピークの大きさも比較的小さなものとなる。ただし、－２スタッター及び＋２スタッターのアンプリコン量を推定することができれば、真のアリルのマイクロサテライトの分布を更に精度よく算出できるようになる。なお、－２スタッター及び＋２スタッターのアンプリコン量を推定するための具体的な手法については、以下で詳しく説明する。
【0052】
（アンプリコン量の推定方法）
次に、既出の図７を参照しながら、スタッターのアンプリコン量を推定する方法について具体的に説明する。なお、各種条件は、第２実施形態で説明したものと同様であるものとする。
次に、既出の図７を参照しながら、スタッターのアンプリコン量を推定する方法について具体的に説明する。なお、各種条件は、第２実施形態で説明したものと同様であるものとする。
【0053】
図７において、すでに説明したように、ＰＣＲを３サイクルした場合、「５回繰り返し」のアンプリコン量は、ａ３となり、－１スタッター（即ち、「４回繰り返し」）のアンプリコン量は、３×ａ２×ｂとなる。この場合、－２スタッター（即ち、「３回繰り返し」）のアンプリコン量は、３×ａ×ｂ２となる。
【0054】
そして、ＰＣＲをｎサイクルした場合、「５回繰り返し」のアンプリコン量は、ａｎとなり、－１スタッター（即ち、「４回繰り返し」）のアンプリコン量は、ｎＣ１×ａｎ－１×ｂとなる。この場合、－２スタッター（即ち、「３回繰り返し」）のアンプリコン量は、ｎＣ２×ａｎ－２×ｂ２となる。－１スタッター比ＳＲ－１は、（「４回繰り返し」のアンプリコン量）／（「５回繰り返し」のアンプリコン量）と表現できるため、前述の値を用いると、ＳＲ－１＝（ｎＣ１×ａｎ－１×ｂ）／（ａｎ）＝ｎＣ１×ｂ／ａ＝ｎ×ｂ／ａと表すことができる。同様に、－２スタッター比ＳＲ－２は、（「３回繰り返し」のアンプリコン量）／（「５回繰り返し」のアンプリコン量）と表現できるため、前述の値を用いると、ＳＲ－２＝（ｎＣ２×ａｎ－２×ｂ２）／（ａｎ）＝ｎＣ２×ｂ２／ａ２＝１／２×ｎ×（ｎ－１）×ｂ２／ａ２≒１／２×（ｎ×ｂ／ａ）２＝１／２×（ＳＲ－１）２と表すことができる。
【0055】
以上のように、第１の割合ａ、第２の割合ｂ、ＰＣＲのサイクル数ｎを用いれば、－２スタッターのアンプリコン量を推定することが可能である。なお、ここでの説明は省略したが、＋２スタッターのアンプリコン量についても、同様の手法で推定することが可能である。上述した例では、－２スタッターが、「－１スタッターが発生する現象が２回発生することで生ずる」ことを前提として説明したが、１回のＰＣＲ過程で－２スタッターが生ずる可能性を考慮してもよい。この場合、上記第１の割合ａ、第２の割合ｂに加えて、－２スタッターが発生する割合を用いればよい。この値は、例えば「ｘ回繰り返し」の２次式で近似すればよい。同様に、１回のＰＣＲ過程で－３スタッターが生ずる可能性を考慮してもよい。この場合、－３スタッターが発生する割合を用いればよい。この値は、例えば「ｘ回繰り返し」の３次式（或いは、３次以上の式）で近似すればよい。
【0056】
（技術的効果）
次に、第３実施形態に係る遺伝子情報処理システム１０によって得られる技術的効果について説明する。
次に、第３実施形態に係る遺伝子情報処理システム１０によって得られる技術的効果について説明する。
【0057】
図９で説明したように、第３実施形態に係る遺伝子情報処理システム１０では、－２スタッター又は＋２スタッターのアンプリコン量を推定し、推定したアンプリコン量を用いて、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布が算出される。このようにすれば、ＰＣＲにおいて発生する－２スタッター又は＋２スタッターの存在を考慮して、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を適切に算出することが可能である。
【0058】
なお、ここまでの実施形態では、－１スタッター、＋１スタッター、－２スタッター、及び＋２スタッターのアンプリコン量を推定する例を挙げたが、繰り返しユニット数が３つ以上少ない又は多いスタッター（例えば、－３スタッターや＋３スタッター等）のアンプリコン量を推定して、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を算出するようにしてもよい。
【0059】
＜第４実施形態＞
第４実施形態に係る遺伝子情報処理システム１０について、図１０及び図１１を参照して説明する。なお、第４実施形態は、上述した第１から第３実施形態と比べて、一部の構成及び動作が異なるのみであり、その他の部分は第１から第３実施形態と同一であってよい。このため、以下では、すでに説明した部分と重複する部分については適宜説明を省略するものとする。
第４実施形態に係る遺伝子情報処理システム１０について、図１０及び図１１を参照して説明する。なお、第４実施形態は、上述した第１から第３実施形態と比べて、一部の構成及び動作が異なるのみであり、その他の部分は第１から第３実施形態と同一であってよい。このため、以下では、すでに説明した部分と重複する部分については適宜説明を省略するものとする。
【0060】
（機能的構成）
まず、図１０を参照しながら、第４実施形態に係る遺伝子情報処理システム１０の機能的構成について説明する。図１０は、第４実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。なお、図１０では、図２で示した構成要素と同様の要素に同一の符号を付している。
まず、図１０を参照しながら、第４実施形態に係る遺伝子情報処理システム１０の機能的構成について説明する。図１０は、第４実施形態に係る遺伝子情報処理システムの機能的構成を示すブロック図である。なお、図１０では、図２で示した構成要素と同様の要素に同一の符号を付している。
【0061】
図１０に示すように、第４実施形態に係る遺伝子情報処理システム１０は、その機能を実現するための処理ブロックとして、第１取得部１１０と、第２取得部１２０と、算出部１３０とを備えて構成されている。そして特に、第４実施形態に係る算出部１３０は、スタッター比推定部１３２を備えている。
【0062】
スタッター比推定部１３２は、第１の割合と、第２の割合と、ＰＣＲ法のサイクル数と、を用いて、ＰＣＲで増幅されたマイクロサテライトのアンプリコン量に対するスタッターのアンプリコン量の比率であるスタッター比を推定可能に構成されている。なお、スタッター比の対象となるスタッターは、繰り返しユニット数が少ないマイナススタッター（例えば、－１スタッターや－２スタッター）であってもよいし、繰り返しユニット数が多いプラススタッター（例えば、＋１スタッターや＋２スタッター）であってもよい。そして、第４実施形態に係る算出部１３０は、スタッター比推定部１３２で推定されたスタッター比を用いて、ＰＣＲ前の遺伝子におけるマイクロサテライトの分布を算出可能に構成されている。
【0063】
（スタッター比の算出方法）
スタッター比は、すでに説明した各スタッターのアンプリコン量を用いて算出することができる。例えば、図７に示す例における－１スタッターのスタッター比ｒ１は、「４回繰り返し（即ち、－１スタッター）のアンプリコン量」÷「５回繰り返しのアンプリコン量」として算出できる。具体的には、－１スタッターのスタッター比ｒ１は、（ｎＣ１×ａｎ－１×ｂ）÷ａｎ＝ｎＣ１×（ｂ／ａ）となる。同様に、－２スタッターのスタッター比ｒ２は、「３回繰り返し（即ち、－２スタッター）のアンプリコン量」÷「５回繰り返しのアンプリコン量」として算出できる。具体的には、－２スタッターのスタッター比ｒ２は、（ｎＣ２×ａｎ－２×ｂ２）÷ａｎ＝ｎＣ２×（ｂ／ａ）２となる。
スタッター比は、すでに説明した各スタッターのアンプリコン量を用いて算出することができる。例えば、図７に示す例における－１スタッターのスタッター比ｒ１は、「４回繰り返し（即ち、－１スタッター）のアンプリコン量」÷「５回繰り返しのアンプリコン量」として算出できる。具体的には、－１スタッターのスタッター比ｒ１は、（ｎＣ１×ａｎ－１×ｂ）÷ａｎ＝ｎＣ１×（ｂ／ａ）となる。同様に、－２スタッターのスタッター比ｒ２は、「３回繰り返し（即ち、－２スタッター）のアンプリコン量」÷「５回繰り返しのアンプリコン量」として算出できる。具体的には、－２スタッターのスタッター比ｒ２は、（ｎＣ２×ａｎ－２×ｂ２）÷ａｎ＝ｎＣ２×（ｂ／ａ）２となる。
【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】
（技術的効果）
次に、第７実施形態に係る遺伝子情報処理システム１０によって得られる技術的効果について説明する。
次に、第７実施形態に係る遺伝子情報処理システム１０によって得られる技術的効果について説明する。
【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】
１０ 遺伝子情報処理システム
１１０ 第１取得部
１２０ 第２取得部
１２１ パラメータ推定部
１３０ 算出部
１３１ アンプリコン量推定部
１３２ スタッター比推定部
１３３ 増幅行列算出部
ａ ＰＣＲの増幅率
ｂ マイナススタッターの発生割合
ｃ プラススタッターの発生割合
ｎ ＰＣＲのサイクル数
１１０ 第１取得部
１２０ 第２取得部
１２１ パラメータ推定部
１３０ 算出部
１３１ アンプリコン量推定部
１３２ スタッター比推定部
１３３ 増幅行列算出部
ａ ＰＣＲの増幅率
ｂ マイナススタッターの発生割合
ｃ プラススタッターの発生割合
ｎ ＰＣＲのサイクル数
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】