JP6417465B2 - 物質の遺伝毒性の評価方法 - Google Patents

Description

本発明は、変異の解析、又は物質の遺伝毒性を評価する方法に関する。

遺伝毒性とは、ＤＮＡを中心とする細胞内の遺伝物質に対する毒性の総称であり、より狭義には、ＤＮＡの損傷や突然変異等を引き起こし、その遺伝情報を変化させる性質（変異原性）をいう。ＤＮＡの遺伝情報は、Ａ（アデニン）、Ｔ（チミン）、Ｇ（グアニン）及びＣ（シトシン）の４種の塩基で構成された塩基配列に保持されている。遺伝毒性物質は、直接的あるいは間接的にＤＮＡに作用し、その塩基配列に質的又は量的に変化を及ぼして、遺伝情報を変化させる。遺伝毒性物質による遺伝情報の変化は、発がんや生殖発生毒性の原因となることが知られており、医薬品、化粧料、各種化学物質等の遺伝毒性を評価することは、公衆の安全のために重要である。

遺伝毒性のメカニズムは多様であり、大きく分けて、ＤＮＡの塩基対情報を別の塩基対に変化させる塩基対置換型変異、ＤＮＡの配列中に短い塩基配列の挿入や欠失を引き起こす短い挿入・欠失変異、及びゲノム配列全体に比較的長い塩基配列の挿入、欠失、転座、逆位などを引き起こし、ゲノム構造を変化させるゲノム構造変化が存在する。特に、短い挿入・欠失変異は、遺伝子のコードするタンパク質の読み枠を変化させる場合、フレームシフト型変異とも呼ばれる。化学物質による遺伝毒性においては、塩基対置換型変異又は短い挿入・欠失変異を引き起こすものが多いと言われる。

これまで、物質の遺伝毒性を評価する方法として、インビトロやインビボのモデルを用いた様々な遺伝毒性試験が開発されてきた。例えば、前述の塩基対置換型変異、あるいはフレームシフト型変異を検出する試験として、Ｂｒｕｃｅ Ｎ．Ａｍｅｓ博士が開発したＡｍｅｓ試験がある（非特許文献１）。Ａｍｅｓ試験では、ヒスチジン生合成遺伝子に変異があり、ヒスチジンを含まない培地では生育できないサルモネラ菌株を使用する。物質曝露によって当該遺伝子に変異が起き、ヒスチジンを合成可能になると、ヒスチジンを含まない培地上でコロニーを形成できるようになる。生じたコロニーを計数することにより、物質の変異原性を確認する。この他に、ゲノム構造変化の有無を検出する試験として、哺乳動物細胞を用いた小核試験などが用いられている（非特許文献２）。さらに、複数の遺伝毒性試験を組み合わせることで、高感度に物質の遺伝毒性の有無を確認することができる。

しかし、前述したような従来の遺伝毒性試験では、遺伝毒性の検出は、変異の量や質を直接反映しない間接的な指標に依存している。そのため、従来の試験では、どのような変異が起きたか、それが何塩基に一つの割合で起きたか等、変異の質的及び量的な情報が詳細には得られない。また、各種試験の間に統一的な指標が存在しないため、異なる試験間の結果の比較は困難である。したがって、複数の変異原間での強度の比較や、メカニズムによる遺伝毒性の分類といった遺伝毒性の体系的理解を進める上で、従来の遺伝毒性試験は十分な情報を提供しない。

次世代シーケンサー等を用いた高スループットシーケンシング技術の遺伝毒性評価への応用が提案されている。Ｍａｓｌｏｖら（非特許文献３）は、高スループットシーケンシングの遺伝毒性評価への応用の方法論を開示している。その方法の一つとして、細胞を物質に曝露した後、単一細胞に由来するゲノム情報が均一な集団を作製し、そのゲノム情報を次世代シーケンサーで取得することで変異部位を同定するものがある。この方法論の実用例として、松田ら（非特許文献４）は、変異原に曝露したＳａｌｍｏｎｅｌｌａ Ｔｙｐｈｉｍｕｒｉｕｍ由来の菌株ＴＡ１００株のシングルコロニーを単離し、その全ゲノム配列を次世代シーケンサーで取得し、リファレンスとなる配列とリード配列間で配列を比較して、一定の頻度で複数のリード配列に同じ塩基の変化が生じている部位を変異部位として検出することで、個々の変異及びその位置を同定したことを報告している。さらに松田ら（非特許文献５）は、シングルコロニーを単離する代わりに希釈した菌株培養液を微量採取し、追培養したものを用いて、非特許文献４と同様の手法で変異を検出する方法を報告している。また別の方法として、次世代シーケンサーを用いて放射線等によるＤＮＡの変異の蓄積を評価する方法が報告されている。具体的には制限酵素サイト等に特異的な配列（タグ配列）に着目し、その出現頻度に基づき、ゲノム中の変異頻度を推定する評価を行っている（特許文献１）。また、セルフリー（ｃｆ）ＤＮＡの各分子に固有のタグ配列を付加し、同一分子から得られる複数のリード配列のコンセンサス配列を得た上で、ゲノム上の同一箇所に複数のリード配列を整列させて比較する変異検出方法が開示されている（非特許文献６、７）。

国際公開公報第２０１４／１７５４２７号

Mortelmans et al., Mutation Research, 2000, 455:29-60 Matsushima et al., Mutagenesis, 1999, 14:569-580 Maslov et al., Mutation Research 2015, 776:136-143 Matsuda, Genes and Environment, 2013, 35:53-56 Matsuda et al., Genes and Environment, 2015、37：15-24 Nucleic Acids Research, 2016, 44(11):e105 Clinical Oncology, 2016, 28:735-738

本発明者は、参照配列の特定の部位を含む複数のリード配列間での比較により、特定部位における塩基配列の変化が、複数のリード配列中で一定の頻度で生じている部位を変異部位として検出する従来法の代わりに、リード配列の各々を参照配列と比較して１つ１つのリード配列から塩基の変異を検出し、その結果を解析することで、変異のパターン及びその頻度を算出する解析方法を見出した。この解析方法は、従来法よりも効率的に、多量の塩基配列情報に基づく変異検出、または高速かつ高感度な変異検出を可能にし、ひいては、細胞集団全体としての量的及び質的な変異の傾向を反映するデータをもたらすことができる。

一実施形態において、本発明は、試験物質の遺伝毒性の評価方法であって、
（１）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
（２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５）取得した各変異を、塩基対の変異パターンに従って分類すること；
（６）該（５）で得られた変異パターンの各々の変異頻度を決定すること、
を含む、方法を提供する。

別の一実施形態において、本発明は、試験物質の遺伝毒性の評価方法であって、
（１’）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
（２’）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３’）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４’）該（３’）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５’）取得した変異の各々について、該参照配列に基づいて、変異前の塩基と、該変異前の塩基の上流及び下流に隣接する塩基とを含むコンテクスト配列を決定すること；
（６’）該（４’）で取得した各変異を、該（５’）で決定したコンテクスト配列及び変異後の塩基の種類に従ってタイプ分けすること；
（７’）該（６’）で得られた変異タイプの各々の変異頻度を決定すること、
を含む、方法を提供する。

さらに別の一実施形態において、本発明は、試験物質の遺伝毒性の評価方法であって、
（１”）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
（２”）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３”）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列における該参照配列に対して塩基が挿入もしくは欠失した部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４”）該（３”）で検出した部位を、挿入もしくは欠失変異を有する変異部位として取得すること；
（５”）取得した変異の各々について、挿入もしくは欠失の塩基長、及び／又は挿入された塩基の種類を決定すること；
（６”）該（５”）で決定された挿入もしくは欠失部位の塩基長及び／又は挿入された塩基の種類ごとの変異頻度を決定すること、
を含む、方法を提供する。

なお別の一実施形態において、本発明は、がん細胞における変異の評価方法であって、
（１）がん細胞集団を試験群とし、そのＤＮＡを取得すること；
（２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５）取得した各変異を、塩基対の変異パターンに従って分類すること；
（６）該（５）で得られた変異パターンの各々の変異頻度を決定すること、
を含む、方法を提供する。

なお別の一実施形態において、本発明は、培養細胞の遺伝情報の評価方法であって、
（１）培養細胞集団を試験群とし、そのＤＮＡを取得すること；
（２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５）取得した各変異を、塩基対の変異パターンに従って分類すること；
（６）該（５）で得られた変異パターンの各々の変異頻度を決定すること、
を含む、方法を提供する。

本発明の方法によれば、一回の解析で、細胞集団内における変異についての量的及び質的な情報を得ることができる。したがって、本発明の方法によれば、従来法と比べてはるかに簡便かつ低コストに、細胞集団レベルでの変異解析が可能になる。本発明の方法は、物質の遺伝毒性評価や癌の評価など、遺伝情報がヘテロな細胞集団における変異の傾向を把握したい場合に特に有効である。

合成ＤＮＡサンプルにおける塩基対の変異パターンの変異コール割合。Ａ：ＧＣ塩基対の変異パターン、Ｂ：ＡＴ塩基対の変異パターン。 合成ＤＮＡサンプルにおける塩基対の変異パターンの変異頻度増加量。Ａ：ＧＣ塩基対、Ｂ：ＡＴ塩基対。 合成ＤＮＡサンプルにおける挿入変異の変異頻度増加量。Ａ：挿入塩基数、Ｂ：挿入塩基の種類。 変異原曝露サンプルにおける塩基対の変異パターンの変異頻度増加量。Ａ：ＧＣ塩基対、Ｂ：ＡＴ塩基対。＊：ｐ＜０．０５、＊＊：ｐ＜０．０１、＊＊＊：ｐ＜０．００１（Ｄｕｎｎｅｔｔ’ｓ ｔｅｓｔ）。 オリジナル塩基の種類に依存した、変異原曝露サンプルにおける塩基対の変異パターンの変異頻度増加量の違い。 変異原曝露サンプルにおける塩基対変異パターンのスペクトル解析結果。 変異原曝露サンプルにおける塩基対変異のシーケンスコンテクスト解析結果。Ｓｉｇｎｉｔｕｒｅ １１：公知のアルキル化剤による変異シグニチャーのパターン、ＥＮＵ：実施例２におけるＥｔｈｙｌｎｉｔｒｏｓｏｕｒｅａ処理による変異パターン。 変異原曝露サンプルにおける塩基対変異のシーケンスコンテクスト解析結果（図７の続き）。Ｓｉｇｎｉｔｕｒｅ １１：公知のアルキル化剤による変異シグニチャーのパターン、ＥＮＵ：実施例２におけるＥｔｈｙｌｎｉｔｒｏｓｏｕｒｅａ処理による変異パターン。

（１．定義）
本明細書において、「変異（又は突然変異）」（mutation）とは、ＤＮＡに生じる突然変異をいい、例えば、ＤＮＡにおける塩基又は配列の欠失、挿入、置換、付加、逆位、及び転座が挙げられる。本明細書における変異は、１塩基の欠失、挿入、置換、付加、ならびに２以上の塩基からなる配列の欠失、挿入、置換、付加、逆位、及び転座を包含する。また本明細書における変異には、コード領域及び非コード領域における変異が含まれ、また発現するアミノ酸の変化を伴う変異及び伴わない変異（サイレント変異）が含まれる。

本発明において評価される物質の「遺伝毒性」とは、該物質が変異を引き起こす性質（いわゆる変異原性）をいう。

本明細書において、「オリジナルフラグメント」とは、解析対象ＤＮＡのフラグメントであって、シーケンシング反応によって配列を読み取られる側の一本鎖ＤＮＡ断片をいう。また本明細書において、変異部位の塩基に関する「オリジナル塩基」とは、それぞれ、オリジナルフラグメントの変異箇所における変異前の塩基をいう。

本明細書中で引用された全ての特許文献、非特許文献、およびその他の刊行物は、その全体が本明細書中において参考として援用される。

（２．細胞集団における変異の解析方法）
高スループットシーケンシングを用いた遺伝毒性評価方法により、変異原により引き起こされた変異の量や質を直接的に評価できると期待される。また、このような方法は、ゲノム配列が利用可能であれば、基本的には全ての生物種に適用可能と考えられる。一方で、非特許文献４のような従来の高スループットシーケンシングを用いた遺伝毒性評価方法では、一部位での塩基の変異の検出のために、その部位を含む多数のリード配列を全て整列させて比較するという方法をとるため、大量の配列情報が必要となり、配列情報の取得及び変異部位の検出に多大な時間、労力及びコストを要する。また、物質曝露による変異は、一般的には非常に低頻度であり集団内の個々の細胞に均等に発生するわけではないため、非特許文献４のような単一細胞の解析では、細胞集団内での変異頻度や、変異原が細胞集団に及ぼす影響を正確に評価することは難しいと考えられる。非特許文献５記載の方法でも、解析に供したサンプルの遺伝情報は依然として均質であったことから、その結果は細胞集団の情報を反映していたとは言い難い。遺伝毒性の評価では、個々の細胞に入った低頻度な変異を如何に検出するか、及びそれに基づいて、集団に対する遺伝毒性をどのように評価するかが重要なポイントである。

非特許文献４又は５のような解析を複数の異なる単一細胞由来のサンプルに対して繰り返せば、細胞集団での事象を反映した情報を得ることができるが、そのための時間、労力やコストは膨大である。特に複数の物質について評価する場合、また用量反応関係を評価する場合においては、この方法は現実的ではない。一方、特許文献１の方法は、比較的コストの低い方法であるが、制限酵素サイト等の特定の配列を利用した変異頻度の推定にとどまっており、ゲノム全体における変異に対しての定性及び定量性が低く、正確な変異情報に基づく遺伝毒性評価を実現可能な方法ではない。

また、従来の高スループットシーケンシングを用いた遺伝毒性評価方法は、哺乳動物細胞の単離培養における技術的ハードルや、ゲノムサイズに起因する変異解析コストの増大を考えた場合、複数の物質の評価への適用は極めて困難である。様々な変異情報を保持した細胞集団における変異原に起因する変異の量的及び質的情報を、簡便かつ低コストに検出する方法の開発が望まれている。

本発明者は、参照配列の特定の部位を含む複数のリード配列間での比較により、特定部位における塩基配列の変化が、複数のリード配列中で一定の頻度で生じている部位を変異部位として検出する従来法の代わりに、リード配列の各々を参照配列と比較して１つ１つのリード配列から塩基の変異を検出し、その結果を解析することで、変異のパターン及びその頻度を算出する解析方法を見出した。この解析方法は、従来法よりも効率的に、多量の塩基配列情報に基づく変異検出、または高速かつ高感度な変異検出を可能にし、ひいては、細胞集団全体としての量的及び質的な変異の傾向を反映するデータをもたらすことができる。

本発明の方法によれば、一回の解析で、細胞集団内における変異についての量的及び質的な情報を得ることができる。したがって、本発明の方法によれば、従来法と比べてはるかに簡便かつ低コストに、細胞集団レベルでの変異解析が可能になる。本発明の方法は、物質の遺伝毒性評価や癌の評価など、遺伝情報がヘテロな細胞集団における変異の傾向を把握したい場合に特に有効である。

したがって、一実施形態において、本発明は、細胞集団における変異を解析する方法を提供する。本発明の方法の基本的な手順は、以下のとおりである：
（Ａ）細胞集団由来のＤＮＡを取得する；
（Ｂ）該ＤＮＡのフラグメント（すなわちオリジナルフラグメント）をシーケンシングし、各フラグメントにつき１つ以上のリード配列を得る；
（Ｃ）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出する；
（Ｄ）該１つ以上のリード配列について検出した部位を、変異部位として取得する；
（Ｅ）該変異部位における変異の情報を取得し、該情報に基づいて変異の傾向を解析する。

本発明の方法で用いる細胞集団は、ホモ集団（例えば単一コロニー由来の細胞集団など、遺伝情報が均一なもの）であっても、ヘテロ集団（遺伝情報が不均一なもの）であっても、遺伝情報の均一性が不明な集団であってもよいが、好ましくは、ヘテロ集団又はそれと推定される細胞集団である。本発明の方法で用いる細胞集団の例としては、動物もしくは植物から採取した検体、及び動物、植物もしくは微生物由来の培養細胞の集団などが挙げられ、好ましくは、動物、植物もしくは微生物の株由来の培養細胞の集団が挙げられる。動物の例としては、好ましくはヒト等の哺乳動物や、カイコ、線虫などが挙げられ、微生物の例としては、好ましくは大腸菌、サルモネラ菌、酵母などが挙げられる。本発明の方法で用いる細胞集団の別の例としては、生体から採取した検体もしくはその培養物、変異原もしくはその候補物質に曝した培養細胞、又は薬物もしくはその候補物質を投与した培養細胞などが挙げられる。

本発明で用いる細胞集団由来のＤＮＡは、該細胞集団から当該分野における通常の方法を用いて抽出又は単離することによって、取得することができる。該抽出又は単離には、例えば、市販のＤＮＡ抽出キットなどを用いることができる。あるいは、抽出又は単離後保存されている該細胞集団由来のＤＮＡを取得し、本発明の方法で使用してもよい。本発明で用いる細胞集団由来のＤＮＡとしては、ゲノムＤＮＡ、ミトコンドリアゲノムＤＮＡ、葉緑体ゲノムＤＮＡ、プラスミドＤＮＡ、ウイルスゲノムＤＮＡなどが挙げられる。このうち、ゲノムＤＮＡが好ましい。

あるいは、細胞集団中のＲＮＡウイルスについて解析する場合は、ＤＮＡの代わりにＲＮＡを取得および解析してもよい。細胞中のＲＮＡは、市販のＲＮＡ抽出キットなど、当該分野における通常の方法で抽出又は単離することができる。あるいは、抽出又は単離後保存されている該細胞集団由来のＲＮＡを取得し、本発明の方法で使用してもよい。本発明の方法においてＲＮＡを取得および解析する場合、本明細書中における「ＤＮＡ」は「ＲＮＡ」と読み替えられ、また塩基Ｔは塩基Ｕと読み替えられる。

本発明の方法で用いる参照配列とは、変異解析の対象であるＤＮＡ中に含まれる既知の配列である。ここで、既知の配列としては、公共のデータベース等に登録されている配列を使用することが好ましいが、上記本発明の方法の（Ｃ）に先立って予めシーケンサー等で配列決定した解析対象ＤＮＡ中の配列であってもよい。該参照配列の領域や長さ、その数は特に限定されず、変異解析の目的に応じてＤＮＡ中から適宜選択され得る。例えば、本発明による変異解析の目的が物質の遺伝毒性評価である場合、参照配列の長さは、特に限定されないが、合計で１，０００ｂｐ以上が好ましく、１０，０００ｂｐ以上がより好ましく、１００，０００ｂｐ以上がさらに好ましく、１，０００，０００ｂｐ以上がなお好ましい。

ＤＮＡのフラグメント化は、超音波処理、酵素処理など、当該分野における通常の方法を用いて行うことができる。調製するフラグメントの長さは、シーケンサーが精度よく読み取れる長さに応じて適宜選択され得る。一般的には、１００〜１０，０００ｂｐが選択され得るが、シーケンサーが精度よく読み取れる限りは１０，０００ｂｐ以上の長さのフラグメントが調製されてもよく、シーケンサーの種類に依存してより適切な範囲が選択され得る。例えば、フラグメントの増幅を行うシーケンシング反応用のシーケンサーにかける場合は、フラグメントの長さは平均長１００〜５００ｂｐが好ましく、平均長１５０〜２００ｂｐがより好ましい。また例えば、一分子リアルタイムシーケンシングを行うシーケンサーにかけるためには、フラグメントの長さは平均長１５０〜１０，０００ｂｐが好ましく、平均長２００〜１，０００ｂｐがより好ましく、平均長２００〜５００ｂｐがより好ましい。

次いで、得られたフラグメントをシーケンシングする。フラグメントのシーケンシングは、後述する参照配列との配列比較に使用すべき部分について行えば足りる。例えば、その配列の少なくとも一部、好ましくは全体が、参照配列のＤＮＡ領域に対応するフラグメントをシーケンシングすればよい。哺乳動物細胞等の場合には、エクソン領域等を選択的にシーケンシングしてもよい。領域の選択には、ＳｕｒｅＳｅｌｅｃｔ（アジレント・テクノロジー社製）等のキットが上市されている。

フラグメントの増幅を行うシーケンシング反応では、フラグメントをＰＣＲ等で増幅するとともに、増幅された各フラグメントの配列を決定する。一分子リアルタイムシーケンシング反応では、フラグメントを増幅せずにその配列を決定する。フラグメントのシーケンシングには、公知のシーケンサーを使用すればよいが、好ましくは高スループットシーケンサー（いわゆる次世代シーケンサー）が用いられる。フラグメントの増幅を行う高スループットシーケンサーとしては、ＨｉＳｅｑ（イルミナ社製）、ＭｉＳｅｑ（イルミナ社製）などが上市されている。また、一分子リアルタイムシーケンスを行う高スループットシーケンサーとしては、Ｐａｃ Ｂｉｏ ＲＳＩＩ（Ｐａｃｆｉｃ Ｂｉｏｓｃｉｅｎｃｅｓ社製）、Ｐａｃ Ｂｉｏ Ｓｅｑｕｅｌ Ｓｙｓｔｅｍ（Ｐａｃｉｆｉｃ Ｂｉｏｓｃｉｅｎｃｅｓ社製）などが上市されている。

フラグメントのシーケンシング方法の詳細な手順は特に限定されないが、より精度の高いシーケンシング方法であることが好ましい。そのような方法の一つとして、フラグメントの増幅を行うシーケンシングにおいては、後述するようにフラグメントの両側から配列を取得し、共通する部分を活用する方法が挙げられる。また、互いに相補的な２つのオリジナルフラグメントから取得した配列を相補鎖間で比較することにより、さらに比較精度を向上させることもできる。そのような技術としては、ＨｉＳｅｑやＭｉＳｅｑのライブラリ調製時に、各フラグメント分子固有のアダプターを付加し、シーケンシング後に該アダプター配列情報をもとに相補鎖間の塩基配列を参照する方法（Proc Natl Acad Sci U S A, 2012, 109(36):14508-14513）などが挙げられる。一分子リアルタイムシーケンシングでは、２本鎖フラグメントの両末端にヘアピン配列のアダプターを付加することで一本の環状にし、該１分子の環状フラグメントを連続で数回シーケンシングし、その配列情報を統合する方法（Nucleic Acids Res, 2010, 38(15):e159）などが挙げられる。

以上のシーケンシングの結果、上述したＤＮＡのフラグメント化により得られたオリジナルフラグメントの各々について、読み取り結果（リード配列）が取得される。各オリジナルフラグメントについて取得するリード配列の数は１つ以上であればよいが、後述する変異傾向の解析における精度向上の観点からは、好ましくは、各フラグメントについてそれぞれ複数のリード配列が取得される。１フラグメントから取得されるリード配列の数は、好ましくは２本以上、より好ましくは１０本以上である。他方、解析の効率の観点からは、リード配列の数は５本以下であることが好ましく、２本以下であることがより好ましい。

取得した該１つ以上のリード配列は、そのまま以降の参照配列との比較に用いることができる。変異解析の精度向上の観点からは、該取得したリード配列の塩基の中から、シーケンシングによる読み取りの信頼度の高いものを抽出することが好ましい。好ましくは、同じ１つのフラグメントから得られた２つ以上のリード配列の対応する塩基の中から、信頼度の高い塩基を抽出する。抽出された信頼度の高い塩基は、本明細書において“コンセンサス塩基”と称される。“コンセンサス塩基”の抽出は、次世代シーケンサーに付属のプログラム等を用いて行うことができる。より詳細な手順としては、例えば、同じ１つのフラグメントから得られた２つのリード配列の間で、対応する位置の塩基の種類が、同じであるか、又は２つのリード配列が相補鎖の場合には相補的である塩基を“コンセンサス塩基”として抽出する方法；同じ１つのフラグメントから得られた２つ以上のリード配列間で塩基を比較し、配列上の各位置について、相補的な鎖を含む場合は相補的な塩基も含めて最大の頻度で出現する塩基を決定し、“コンセンサス塩基”として抽出する方法；リード配列間の対応する位置にある塩基の中でシーケンサーでの読み取り精度（クオリティ値）の最も高い塩基を“コンセンサス塩基”として採用する方法；クオリティ値や塩基の出現頻度等を基に、確率論的に“コンセンサス塩基”を決定する方法；同じ１つのフラグメントから得られた全てのリード配列の対応する塩基が一致した場合にその塩基を“コンセンサス塩基”とする方法；またはこれらを組み合わせた方法、などが挙げられる。本発明の方法において、該“コンセンサス塩基”の抽出は、必要に応じて、後述するリード配列の参照配列へのマッピングの前に実施してもよく、又はマッピング後に実施してもよい。

次いで、シーケンシングにより得られた各オリジナルフラグメントについての該１つ以上のリード配列を、それぞれ参照配列に対してマッピングし、配列を比較する。比較の結果、該各リード配列と参照配列とで塩基がマッチしない部位を検出する。好ましくは、当該比較は、上述したリード配列上の“コンセンサス塩基”についてのみ行われ得る。該「マッチしない部位」のタイプとしては、例えば、参照配列に対してリード配列上の塩基の種類が異なる部位（置換部位）、参照配列に対してリード配列上の塩基が欠失している部位（欠失部位）、参照配列に対してリード配列上に塩基が挿入されている部位（挿入部位）などが挙げられる。

検出した「マッチしない部位」は、解析対象ＤＮＡにおける変異部位として取得される。より詳細には、変異部位における変異の情報を取得する。取得される情報としては、例えば、該マッチしない部位（変異部位）のタイプ（例えば、置換部位か、欠失部位か、挿入部位か）、該部位における塩基の種類及び変異前の塩基（又は塩基対）の種類（例えば、参照配列上で該部位に対応する位置にある塩基の種類）、該部位の両隣の塩基の種類（例えば、参照配列上で該部位に対応する位置の両隣にある塩基の種類）などが挙げられるが、これらに限定されない。好ましい実施形態においては、該マッチしない部位が置換部位であれば、該部位における塩基の種類、及び変異前の塩基の種類の情報が取得され；マッチしない部位が欠失部位であれば、変異前の塩基の種類、及び該欠失部位の両隣の塩基の種類の情報が取得され；マッチしない部位が挿入部位であれば、該部位における塩基の種類、及び該挿入部位の両隣の塩基の種類の情報が取得される。

リード配列と参照配列との比較からは、リード配列と参照配列とで塩基がマッチした部位も検出され得る。これらの「マッチした部位」は、解析対象ＤＮＡにおいて変異がない部位として取得され得る。これらの変異がない部位に関する情報を取得することができる。取得される情報としては、例えば、該部位における塩基の種類や、該部位の両隣の塩基の種類などが挙げられる。

該１つ以上のリード配列の各々について、上述した変異部位における変異の情報、又変異がない部位に関する情報の取得を行う。取得した情報を集めて、変異解析用データベースを作成することができる。例えば、各リード配列から取得した変異部位に関する情報を全て集めたデータベースを作成してもよく；各リード配列から取得した変異情報を、変異部位のタイプごとに分類したデータベースを作成してもよく；各リード配列から取得した変異情報を、変異部位における変異前（例えば参照配列）もしくは変異後の塩基の種類ごとに分類したデータベースを作成してもよく；各リード配列から取得した変異情報を、変異部位の塩基長（例えば挿入、欠失又は置換部位の長さ）ごとに分類したデータベースを作成してもよく；さらに、これらの分類を組み合わせたデータベースを作成してもよい。あるいは、変異部位に関する情報と変異がない部位に関する情報とを統合したデータベースを作成してもよい。例えば、変異部位と変異がない部位の情報とを、変異前（例えば参照配列）の塩基の種類（Ａ、Ｔ、Ｇ及びＣ）ごとにまとめたデータベースを作成してもよい。あるいは、変異部位のゲノム上での位置を特定し、遺伝子のコード領域または非コード領域に該当するか、さらに遺伝子のコード領域の場合には、イントロンかエクソンか、ＲＮＡに転写される側の鎖か、そうでないかなどの情報とともにデータベースを作成してもよい。

上述した「マッチしない部位」の検出、ならびに変異部位又は変異がない部位に関する情報の取得は、シーケンシングで得られた全リード配列に対して行ってもよいが、一部のリード配列に対して行ってもよい。当該検出に用いるリード配列の総量（用いるリード配列の合計長）は、その後の変異の傾向の解析を可能にする量であれば特に限定されないが、変異頻度の逆数以上の塩基長であることが好ましく、変異頻度の逆数の１００倍以上の塩基長であることがより好ましい。例えば、後述する実施例２における変異の頻度が約１／１０⁵ｂｐのオーダーであることから、用いるリード配列の長さの合計は、１×１０⁵ｂｐ以上であることが好ましく、１×１０⁷ｂｐ以上あることがより好ましく、１×１０⁹ｂｐ以上であることがさらに好ましい。また、変異の頻度が１／１０⁶ｂｐのオーダーである変異の傾向を解析するためには、検出に用いるリード配列の総量は、前記量の１０倍であることが好ましく、変異の頻度が１／１０⁴ｂｐのオーダーである変異の傾向を解析するためには、検出に用いるリード配列の総量は、前記量の１／１０倍であることが好ましい。他方、解析の効率の点から、検出に用いるリード配列の総量は、変異頻度の逆数の１万倍以下の塩基長であることが好ましく、１０００倍以下の塩基長であることがより好ましく、１００倍以下の塩基長であることがより好ましい。例えば、検出に用いるリード配列の総量は、１×１０¹⁰ｂｐ以下であることが好ましく、１×１０⁹ｂｐ以下であることがより好ましく、１×１０⁸ｂｐ以下であることがさらに好ましく、１×１０⁷ｂｐ以下であることがよりさらに好ましく、１×１０⁶ｂｐ以下であることがなお好ましい。また、該変異解析用データベースは、取得した変異部位又は変異がない部位の全てに関する情報に基づいて作成されてもよいが、その後の変異の傾向の解析を可能にする限り、一部の部位の情報のみに基づいて作成されてもよい。

従来の方法（例えば、非特許文献４、５）では、参照配列の特定の部位に対応する複数のリード配列を取得した後、該複数のリード配列間で同じ部位に同じタイプのマッチしない塩基が一定頻度で見られたときに、当該部位を解析対象ＤＮＡにおける変異部位として決定していた。この方法では、低頻度の変異を見逃す可能性がある。またこの方法では、変異検出するＤＮＡ領域が該複数のリード配列がオーバーラップする限定された長さの領域に制限され、かつリードをオーバーラップさせるために多くのデータを必要とするため、ＤＮＡの広い領域にわたって解析し、全体での変異の傾向を把握するには膨大な時間と労力を要する。

他方、本発明の方法では、そのようなリード配列間でのマッチしない塩基の出現頻度に基づく変異部位の決定は行われない。本発明の方法では、基本的には参照配列の特定の部位に対応する１つ以上のリード配列についてそれぞれ、参照配列との比較に基づいた変異情報を取得し、必要に応じてそれらの情報を分類して、データベースを作成する。該データベースに基づいて、解析対象ＤＮＡにおける変異の傾向を解析する。例えば、該データベースに含まれる任意の要素を標本集団とした統計解析（例えば変異頻度解析、変異のパターン解析等）を行うことができる。

本発明の方法では、各リード配列についてそれぞれ変異情報を取得するため、低頻度の変異を見逃すことなく検出することができる。また本発明の方法では、変異検出に用いたリード配列のいずれかに対応するＤＮＡ上の広範な領域についての変異の検出と解析が可能になる。したがって、本発明の方法によれば、従来の方法よりも高速かつ高感度に変異を検出することができるので、より効率よくかつより正確な変異解析を行うことができる。

（２−１．次世代シーケンサーを用いた各フラグメントにおける変異の検出）
本発明の方法の好ましい実施形態として、ＰＣＲを行う次世代シーケンサーを用いたＤＮＡフラグメントのシーケンシング、及び参照配列との比較による解析対象ＤＮＡにおける変異の解析の手順を以下に詳述する。

シーケンシングされる解析対象ＤＮＡに由来するフラグメントは、シーケンシングのためのアダプター配列を両端に付加される。アダプターの付加されたフラグメントは、次世代シーケンサーでのシーケンシングに際し、ＰＣＲ法で検出可能な量まで増幅される。増幅されたフラグメントは配列を読み取られ、読み取られた配列はリード配列として出力される。本発明の好ましい実施形態においては、各増幅されたフラグメントに対して、２本のリード配列（リード１、リード２）を取得する。このとき、シーケンシングで読み取られたオリジナルフラグメントの配列側に対応するのがリード１であり、その相補鎖側に対応するのがリード２である。シーケンサーのリード長の２倍未満のサイズでフラグメントを調製した場合、各増幅されたフラグメントのリード１及びリード２は、共通領域として該フラグメントの少なくとも一部を含み、さらに、それぞれその上流領域又は下流領域を含む。読み取ったリード１とリード２を共通領域で重ね合わせることにより、各増幅されたフラグメントについて１本の合成リード配列を構築する。２本のリード配列からの合成リード配列の構築は、ＰＥＡＲ（Bioinformatics, 2014, 30(5):614-620）、ＦＬＡＳＨ（Bioinformatics, 2011, 27(21):2957-2963）、ＰＡＮＤＡｓｅｑ（BMC Bioinformatics, 2012, 13:31）等のソフトウェアにより実行することができる。

次いで、各々のリード配列を参照配列上にマッピングし、比較する。好ましい実施形態においては、比較精度の向上のため、当該比較は、リード配列上の上述した“コンセンサス塩基”についてのみ行われ得る。別の好ましい実施形態においては、比較精度の向上のため、比較に供するリード配列として合成リード配列を用いる。より好ましくは、合成リード配列における参照配列と比較する領域を、リード１とリード２の重複部分に限定し、かつ比較に用いる合成リード配列上の塩基を、リード１とリード２の間で塩基の種類が相補的であるもの（すなわち“コンセンサス塩基”）に限定する。これにより、シーケンシング時のエラーが配列比較に及ぼす悪影響を低減することができる。該“コンセンサス塩基”への限定は、参照配列へのマッピングを行う前に実施してもよいが、マッピング後に実施することもできる。

該参照配列へのマッピングと比較によって、各々のリード配列における、参照配列と塩基がマッチしない部位（変異部位）を検出することができる。さらに、該変異部位のタイプ（置換部位か、欠失部位か、挿入部位か）、該部位における塩基の種類、変異前の塩基の種類、該部位の両隣の塩基の種類、などの変異情報を取得することができる。これらの手順を、シーケンシングで得られた各オリジナルフラグメント由来のリード配列について行うことにより、変異情報を集めて、各リード配列からの変異情報のデータベースを作成することができる。上記の手順は、１つ１つのリード配列について順次行うこともできるが、複数のリード配列について並行して行ってもよい。

上述したリード配列の参照配列へのマッピング、比較領域の絞り込み、参照配列と塩基がマッチしない部位の表示、及び該部位における変異情報の取得の手順は、例えば、マッピングはＢｏｗｔｉｅ２ソフトウェア（Nature Methods, 2012, 9(4):357-359）や、ＢＷＡソフトウェア（Bioinformatics, 2009, 25(14):1754-1760）、比較領域の絞り込み、参照配列とマッチしない塩基の表示はＳａｍｔｏｏｌｓソフトウェア（Bioinformatics, 2009, 25(16):2078-2079）、及び該部位における変異情報の取得はＰｙｔｈｏｎ等のプログラミング言語を用いて作成した参照配列と異なる塩基を検出するプログラム等により実行することができる。但し、上記本発明の方法の手順を実行するためのソフトウェア又はプログラミング言語はこれらに限定されない。

（２−２．細胞集団における変異の解析）
さらに、取得された変異情報のデータベースに基づいて、細胞集団における変異の傾向を調べることができる。好ましくは、本発明により解析される変異としては、ＤＮＡの塩基対を別の塩基対に変化させる塩基対置換型変異、及びＤＮＡの配列中に短い塩基配列の挿入や欠失を引き起こす短い挿入・欠失変異が挙げられる。塩基対置換型変異としては、１塩基対置換型変異、及び２塩基対又は３塩基対以上が置換した多塩基対置換型変異が挙げられる。このうち本発明では、好ましくは１塩基対置換型変異が解析される。本発明により、これらの変異の変異パターン及び変異頻度を決定することができる。以下、その解析手順を詳述する。

（２−２−１．塩基対置換型変異の解析）
一実施形態においては、１塩基対置換型変異が解析される。本実施形態においては、上述したように各オリジナルフラグメントからの１以上のリード配列をそれぞれ参照配列と比較し、各リード配列における該参照配列と塩基がマッチしない部位を検出する。これらの部位は、参照配列に対して塩基対置換型変異を有する変異部位として取得される。次いで、検出した変異部位の塩基の種類と変異前の塩基に基づいて、各変異を塩基の変異パターンに従って分類する。次いで、該塩基の変異パターンの各々について、出現頻度を決定する。これらの手順は、上述したＰｙｔｈｏｎ等のプログラミング言語を用いて作成したプログラム等を用いて行うことができる。

より詳細な例においては、リード配列に含まれる各塩基を下記(i)〜(iv)に分ける。
(i) 参照配列上の塩基がＡである位置に存在する塩基
(ii) 参照配列上の塩基がＴである位置に存在する塩基
(iii)参照配列上の塩基がＧである位置に存在する塩基
(iv) 参照配列上の塩基がＣである位置に存在する塩基
上記(i)及び(ii)は、参照配列の塩基対がＡＴであった部位に存在する塩基であり、上記(iii)及び(iv)は、参照配列の塩基対がＧＣであった部位に存在する塩基である。これらの塩基の中から、参照配列と塩基がマッチしない（すなわち塩基対置換変異している）ものを検出する。次いで、検出された変異した塩基の各々について、上記(i)〜(iv)の分類に基づいて、参照配列の塩基情報から該変異部位に存在していた変異前の塩基対を求め、また各リード配列の塩基情報に基づいて変異後の塩基対を求める。これらのデータから、各変異を、変異前の塩基対がＡＴであった場合について[ＡＴ→ＴＡ、ＡＴ→ＣＧ、及びＡＴ→ＧＣ]の３パターン、変異前の塩基対がＧＣであった場合について[ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴ]の３パターンの、全部で６つの塩基対の変異パターンに分類することができる。さらに、各変異パターンに属する変異の総数、及び解析した塩基の総数に基づいて、各変異パターンの出現頻度を決定することができる。例えば、ＡＴ、ＧＣ塩基対それぞれについての解析した塩基の総数に基づいて、各々の塩基対ごとに３種類の変異パターンの出現頻度を算出することができる。

さらに、上記で得られた各変異の変異パターンを、オリジナル塩基に従ってさらに分類することができる。変異前の塩基対がＡＴであった部位の変異のオリジナル塩基はＡ又はＴであり、変異前の塩基対がＧＣであった部位の変異のオリジナル塩基はＧ又はＣである。したがって、上記６つの塩基対の変異パターンのそれぞれは、オリジナル塩基に従ってさらに２つに分けることができる。このような分類は、細胞からのＤＮＡの抽出又は単離の過程において生じた塩基の修飾によって生じたシーケンシングの読み取りエラーの除去に有用である。特に、Ｇ塩基は、ＤＮＡ調製過程で酸化による化学修飾を受け、ＧをＴとして読み間違えるエラーを起こしやすいことが知られている。通常であれば、塩基対の変異により該塩基対を構成する双方の塩基が変化するため、オリジナル塩基に従って２つに分けた塩基対の変異パターンは、同等の変異頻度を示すはずである。変異頻度がどちらかのオリジナル塩基に偏っていれば、塩基の修飾に起因するシーケンシングエラーが示唆される。

本発明の別の一実施形態においては、多塩基対置換型変異が解析される。多塩基対置換型変異としては、例えば、２塩基対置換型変異及び３塩基対置換型変異が挙げられる。多塩基対置換型変異の解析の場合には、例えば、変異前の塩基配列に応じて変異パターンを分類し（例えば２塩基対置換型においては４×４＝１６通り）、次いで、各変異パターンに属する変異の総数、及び解析した変異の総数に基づいて、各変異パターンの出現頻度を決定することができる。

（２−２−２．シーケンスコンテクスト解析）
近年、がん細胞のゲノム上に蓄積された変異情報から、変異原が引き起こした変異の要素（変異シグニチャー）を数学的手法で抽出するアプローチが提案されている。ヒトの様々ながんゲノムに蓄積された変異情報から、様々な変異シグニチャーが同定されている（Cell Rep, 2013, 3:246-259）。変異シグニチャーの理論においては、塩基対置換型変異の変異パターンを該塩基対が位置する前後のシーケンスコンテクストに基づいて分類する。シーケンスコンテクスト解析を行うことにより、塩基対置換型変異のより詳細な解析が可能になる。

したがって、別の一実施形態として、１塩基対置換型変異のシーケンスコンテクスト解析の手順を示す。本実施形態においては、まず、上述したようにリード配列をそれぞれ参照配列と比較することによって、各リード配列における１塩基対置換型変異を検出する。次いで、検出した各変異について、該参照配列に基づいて、変異前の塩基と、該変異前の塩基の上流及び下流に隣接する塩基とを含む配列（いわゆるコンテクスト）を決定する。続いて、各変異を、塩基対の変異パターン及び該コンテクストに従ってタイプ分けする。すなわち、検出した変異を、上述した（２−２−１．）と同様の手順で、６つの塩基対の変異パターン［ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴ］に分ける。一方で、検出した各変異を、コンテクストに従って分類する。例えば、変異部位の両隣の１塩基ずつを含めた３塩基長のコンテクストは、４×４の１６群［例えば、Ｃからの変異の場合、ＡＣＡ、ＡＣＣ、ＡＣＧ、ＡＣＴ、ＣＣＡ、ＣＣＣ、ＣＣＧ、ＣＣＴ、ＧＣＡ、ＧＣＣ、ＧＣＧ、ＧＣＴ、ＴＣＡ、ＴＣＣ、ＴＣＧ、及びＴＣＴ］に分類される。結果、各変異は、塩基対の変異パターンとコンテクストに従って、全部で９６（４×６×４）のタイプに分類される。この解析に使用するコンテクストの配列は、変異前の塩基と、該変異前の塩基の上流に隣接する１以上の塩基と、該変異前の塩基の下流に隣接する１以上の塩基からなるものであればよい。またコンテクストの長さは、３塩基以上であればよいが、それに限定されるものではなく、必要に応じてさらに長いコンテクストを解析することも可能である。例えば、変異部位の両隣の２塩基ずつを含めた５塩基長のコンテクストに従うと、各変異は２５６群（４×４×４×４）に分類され、この分類と６つの塩基対パターンにより、各変異は最終的に全部で１５３６（４×４×６×４×４）のタイプに分類される。さらに変異部位の両隣のｎ塩基ずつを含めた２ｎ＋１塩基長のコンテクストに従うと、各変異は４²ⁿ群に分類され、この分類と６つの塩基対パターンにより、各変異は最終的に全部で４²ⁿ×６個のタイプに分類される。次いで、各変異タイプに属する変異の総数、及び解析した塩基の総数に基づいて、上記変異タイプの各々の変異頻度を決定することができる。

（２−２−３．短い挿入・欠失変異の解析）
さらに別の一実施形態においては、短い挿入・欠失変異が解析される。本実施形態においては、上述したようにリード配列をそれぞれ参照配列と比較することによって、各リード配列における該参照配列に対して塩基が挿入又は欠失されている部位を検出する。これらの部位は、参照配列に対して挿入又は欠失変異を有する変異部位として取得される。さらに、取得された各変異について、変異のタイプ（挿入変異か又は欠失変異か）、該挿入もしくは欠失部位の塩基長、及び／又は挿入もしくは欠失した塩基の種類を決定する。本実施形態で検出される挿入もしくは欠失部位としては、好ましくは挿入もしくは欠失した塩基の長さが１０ｂｐ以下、より好ましくは１〜５ｂｐである部位がよいが、これに限定されない。特定の塩基長の挿入もしくは欠失部位を検出する手順は、上述したＰｙｔｈｏｎ等のプログラミング言語を用いて作成したプログラムを用いて行うことができる。さらに、各リード配列と参照配列との比較によって、挿入もしくは欠失した塩基の種類を同定することができる。これらにより、各リード配列における挿入もしくは欠失部位の塩基長、及び／又は挿入もしくは欠失部位の塩基の種類を決定することができる。次いで、挿入もしくは欠失の頻度を、塩基長及び／又は塩基の種類ごとに決定する。例えば、各リード配列について取得した挿入もしくは欠失変異を塩基長ごとに分類し、それぞれの頻度を決定することができる。また例えば、挿入もしくは欠失した塩基をその種類（Ａ、Ｔ、Ｇ、及びＣ）ごとに分類し、それぞれの頻度を決定することができる。さらに、該塩基長及び塩基の種類による分類を組み合わせたより細かい変異の分類を行い、それぞれの頻度を決定することができる。

（２−３．変異頻度増加量の解析）
上述の手順に従ってリード配列と参照配列との比較により検出された変異は、シーケンシングにおける塩基読み取りエラーを含み得る。より高精度な変異解析のためには、このエラー成分を除去することが好ましい。エラー成分の除去は、解析対象の細胞集団の変異頻度から、対照の細胞集団の変異頻度を差し引くことによって行うことができる。さらに、解析対象の細胞集団が特定の条件に曝された細胞集団である場合、解析対象と対照との変異頻度の差を求めることによって、当該特定の条件が変異頻度に及ぼす影響の解析が可能になる。例えば、解析対象とする細胞集団が、特定の条件に曝された細胞集団、例えば、変異原に曝された細胞集団、薬物を投与された細胞集団などである場合、これらの条件に曝されていない同じ細胞集団を対照の細胞集団とする。この対照の細胞集団について、上記と同様の手順で塩基対置換型変異の解析、シーケンスコンテクスト解析、又は挿入・欠失変異の解析を行い、変異頻度を決定する。得られた対照の細胞集団の変異頻度を、解析対象の細胞集団の変異頻度から差し引く。これにより、解析対象の細胞集団の変異頻度からシーケンシングエラーの成分を除去するとともに、当該特定の条件による変異頻度の増加の有無、又は当該特定の条件下での対照に対する変異頻度増加量を調べることができる。さらに、変異頻度の解析を、上述したオリジナル塩基に従った分類に基づいて行うと、塩基の修飾に起因するシーケンシングエラーを検出することができるため好ましい。

（３．応用）
上述した本発明による変異解析方法により、細胞集団に生じた変異を定量的かつ定性的に解析することができる。本発明の解析方法は、変異に関連する様々な解析又は評価に応用することができる。代表的な応用例としては、物質の遺伝毒性の評価、腫瘍発生に伴う変異の評価方法（例えば、がん細胞における変異の評価、及びｃｆＤＮＡにおける変異の評価方法）、ならびに培養細胞の品質管理（例えば変異の有無又は変異タイプの評価等の遺伝情報の評価）である。

したがって、本発明の好ましい実施形態として、試験物質の遺伝毒性の評価方法が提供される。該方法の具体的手順を以下に説明する。

（３−１．物質の遺伝毒性の評価方法）
（３−１−１．塩基対置換型変異の解析に基づく評価）
本発明による試験物質の遺伝毒性の評価方法の一実施形態は、上述した塩基対置換型変異の解析に基づく。当該方法は、
（１）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
（２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５）取得した各変異を、塩基対の変異パターンに従って分類すること；
（６）該（５）で得られた変異パターンの各々の変異頻度を決定すること、
を含む。

上記工程（１）〜（６）については、上記（２．）、特に（２−２−１．）で述べたとおりである。より詳細には、１塩基対置換型変異を解析する場合、工程（５）では、各変異を、参照配列の塩基対がＡＴである部位については[ＡＴ→ＴＡ、ＡＴ→ＣＧ、及びＡＴ→ＧＣ]の３パターンの変異パターンに分類し、また参照配列の塩基対がＧＣである部位については[ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴ]の３パターンに分類する。好ましくは、これらを組み合わせて全部で６つの塩基対の変異パターンに分類する。さらに、これら６つの塩基対の変異パターンを、オリジナル塩基の種類（ＡＴであればＡであるかＴであるか、ＧＣであればＧであるかＣであるか）に基づいて、それぞれさらに２群に分け、全部で１２の変異パターンに分類してもよい。次いで工程（６）において、工程（５）で決定した変異パターンの各々について、変異頻度を決定する。これにより、各塩基対の変異パターンの変異頻度を決定することができる。

好ましい実施形態においては、上記方法はさらに以下を包含する：
（７）該試験物質に曝露していない該細胞集団を対照群とし、該（１）〜（６）と同様の手順で、該対照群における各塩基対の変異パターンの変異頻度を決定すること；
（８）該（６）で得られた試験群における各変異パターンの変異頻度から、該（７）で得られた対照群における各変異パターンの変異頻度を引き算すること。
これにより、シーケンシングエラーの影響を除去した試験群における変異頻度増加量を求めることができる。

（３−１−２．シーケンスコンテクスト解析に基づく評価）
本発明による試験物質の遺伝毒性の評価方法のさらなる一実施形態は、上述したシーケンスコンテクスト解析に基づく。当該方法は、
（１’）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
（２’）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３’）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４’）該（３’）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５’）取得した変異の各々について、該参照配列に基づいて、変異前の塩基と、該変異前の塩基の上流及び下流に隣接する塩基とを含むコンテクスト配列を決定すること；
（６’）該（４’）で取得した各変異を、該（５’）で決定したコンテクスト配列及び変異後の塩基の種類に従ってタイプ分けすること；
（７’）該（６’）で得られた変異タイプの各々の変異頻度を決定すること、
を含む。

上記工程（１’）〜（７’）については、上記（２．）、特に（２−２−２．）で述べたとおりである。より詳細には、工程（６’）では、各変異を、上述した塩基対の６つの変異パターン[ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴ]と、変異部位の両隣の塩基の種類に従う１６群[たとえば、Ｃからの変異の場合、ＡＣＡ、ＡＣＣ、ＡＣＧ、ＡＣＴ、ＣＣＡ、ＣＣＣ、ＣＣＧ、ＣＣＴ、ＧＣＡ、ＧＣＣ、ＧＣＧ、ＧＣＴ、ＴＣＡ、ＴＣＣ、ＴＣＧ、及びＴＣＴ]に基づいて、全部で９６のタイプに分類する。次いで工程（７’）において、工程（６’）で決定した変異タイプの各々について、変異頻度を決定する。これにより、変異のタイプ及び変異頻度を決定することができる。

好ましい実施形態においては、上記方法はさらに以下を包含する：
（８’）該試験物質に曝露していない該細胞集団を対照群とし、該（１’）〜（７’）と同じ手順で、対照群における各変異タイプの変異頻度を決定すること；
（９’）該（７’）で得られた試験群における各変異タイプの変異頻度から、該（８’）で得られた対照群における各変異タイプの変異頻度を引き算すること。
これにより、シーケンシングエラーの影響を除去した試験群における変異頻度増加量を求めることができる。

（３−１−３．短い挿入もしくは欠失変異の解析に基づく評価）
本発明による試験物質の遺伝毒性の評価方法のさらなる一実施形態は、上述した短い挿入もしくは欠失変異の解析に基づく。当該方法は、
（１”）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
（２”）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３”）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列における該参照配列に対して塩基が挿入もしくは欠失した部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４”）該（３”）で検出した部位を、挿入もしくは欠失変異を有する変異部位として取得すること；
（５”）取得した変異の各々について、挿入もしくは欠失の塩基長、及び／又は挿入された塩基の種類を決定すること；
（６”）該（５”）で決定された挿入もしくは欠失部位の塩基長及び／又は挿入された塩基の種類ごとの変異頻度を決定すること、
を含む。

上記工程（１”）〜（６”）については、上記（２．）、特に（２−２−３．）で述べたとおりである。より詳細には、工程（３”）では、参照配列との比較により、各リード配列における挿入もしくは欠失部位であって、その挿入もしくは欠失した塩基の長さが好ましくは１０ｂｐ以下、より好ましくは１〜５ｂｐである部位を検出する。さらに挿入又は欠失が検出された場合、各リード配列の配列と参照配列との比較によって、挿入又は欠失した塩基の種類を同定することができる。これにより、各リード配列における挿入もしくは欠失部位の塩基長、及び／又は挿入もしくは欠失した塩基の種類を決定することができる（工程（５”））。次いで工程（６”）において、工程（５”）で決定された挿入もしくは欠失部位の塩基長、及び／又は挿入もしくは欠失した塩基の種類の各々について、変異頻度を決定する。これにより、変異パターン及び変異頻度を決定することができる。

好ましい実施形態においては、上記方法はさらに以下を包含する：
（７”）該試験物質に曝露していない該細胞集団を対照群とし、該（１”）〜（６”）と同様の手順で、対照群における挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度を決定すること；
（８”）前記（６”）で得られた試験群における挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度から、該（７”）で得られた対照群における変異頻度を引き算すること。
これにより、シーケンシングエラーの影響を除去した試験群における変異頻度増加量を求めることができる。

本発明による試験物質の遺伝毒性の評価方法によれば、試験物質への曝露によってそれぞれ異なる変異を生じている可能性のある細胞の集団について、該細胞集団レベルでの変異の傾向を解析することができる。したがって本発明によれば、従来の単一細胞ベースでの解析方法（例えば非特許文献４、５）と比べて、生体内における該試験物質による実際の影響により近い情報を得ることができる。また本発明によれば、細胞集団に対して試験物質が及ぼす影響を、ＤＮＡ上の個々の塩基のレベルで定量的かつ定性的に解析することが可能であるため、試験物質の遺伝毒性に関するより詳細な情報を得ることができる。

（３−２．腫瘍発生に伴う変異の評価方法）
（３−２−１．がん細胞における変異の評価方法）
本発明の別の好ましい実施形態として、がん細胞における変異の評価方法が提供される。該方法の具体的手順は、基本的には上述した試験物質の遺伝毒性の評価方法と同じである。ただし試験群としては、がん細胞の集団を使用する。あるいは、がん細胞集団の代わりに、がんが疑われる細胞や、がん化リスクを評価したい細胞の集団を試験群として用いてもよい。対照群としては、非がん細胞（例えば正常細胞）の集団やがん化リスクの低い細胞の集団を使用する。当該方法は、がん種に特有の変異の傾向を特定したり、がん化リスクを評価したり、がんの進行度や悪性度を確認するために有用である。がん細胞における変異の評価に関しては、（２−２−２．）に示したようなシーケンスコンテクスト解析により、より詳細な解析が可能になる。これまでにも、ヒトのがんゲノム解析において、変異を、３塩基のコンテクストに基づいて９６通り（４×６×４）、あるいは５塩基のコンテクストに基づいて１５３６通り（４×４×６×４×４）に分類したことが報告されている（Cell Rep, 2013, 3:246-259）。

本発明によるがん細胞における変異の評価方法によれば、従来のゲノム上の特定部位での変異の解析方法、すなわち、ゲノム上の同一箇所に対する複数のリード配列を整列させて比較する方法（例えば非特許文献４、５）のようながん組織中で選択を受けた結果としてがん組織中の細胞に一定以上の割合で存在する変異のみを抽出する方法と異なり、がん細胞集団全体に生じている当該がん種特有の変異の傾向を特定することができる。がん細胞中の変異の量や質は、がんの種類によって様々であることが報告されている（Nature, 2013, 500(7463):415-421）。本発明の方法によれば、がん細胞集団における変異の傾向を、定量的かつ定性的に解析することが可能であるため、がんの進行度や、種類の診断に有用である可能性が考えられる。

さらに、これまでの変異シグニチャーの解析では、特定の原因によって腫瘍を誘発したと考えられる集団に対して、その集団内の各人から得られたゲノムＤＮＡを解読し、従来法（例えば非特許文献４、５）に従ってがん組織中に一定以上の割合存在する変異を抽出した後、抽出した変異に対してシーケンスコンテクスト解析を行い、次いで各人から得られた変異情報を集計して、集団内の変異シグニチャーを同定する。しかし、個人の診断において変異シグニチャー解析を行う場合、同定される変異数が少ないため、得られたシーケンスコンテクストと既知の変異シグニチャーとの類似性の判断が難しい可能性がある。一方で、本発明によれば、図７及び８に示すように、１回の次世代シーケンシング解析によって、変異シグニチャーとの類似性を確認可能な品質のシーケンスコンテクストデータを得ることが可能である。これは、従来の方法と異なり、一定以上の割合を占めていない変異を解析対象に加えることによって、シーケンスコンテクストの解析に十分な量の変異数を確保できるためである。従って本発明は、個人の変異におけるシーケンスコンテクストのハイスループット解析に活用できると考えられる。

（３−２−２．ｃｆＤＮＡにおける変異の評価方法）
本発明の別の好ましい実施形態として、血中セルフリーＤＮＡ（ｃｆＤＮＡ）の解析が提供される。ｃｆＤＮＡは、ヒトの腫瘍形成の低侵襲な診断方法として着目されている。ｃｆＤＮＡは血漿や血清、尿などのリキッドバイオプシーから得られる。該方法のｃｆＤＮＡ解析への適用に関する具体的手順は、基本的には上述した試験物質の遺伝毒性の評価方法と同じである。ただし試験群のＤＮＡとしては、がん細胞集団から取得したＤＮＡの代わりに、がん患者のヒト等から採取したリキッドバイオプシー中のｃｆＤＮＡを使用する。対照群としては、例えば健常なヒト由来のｃｆＤＮＡや、予め採取しておいたがんを患う前の時点における同一のヒトのｃｆＤＮＡを用いる。当該方法は、がん患者に特有のｃｆＤＮＡ中の変異の傾向を特定したり、がんの進行度や悪性度を確認したりするのに有用である。

本発明によるｃｆＤＮＡにおける変異の評価方法によれば、従来の解析方法、すなわち、ｃｆＤＮＡの各分子に固有のタグ配列を付加し、同一分子から得られる複数のリード配列のコンセンサス配列を得た上で、ゲノム上の同一箇所に複数のリード配列を整列させて比較する方法（非特許文献６、７を参照）のような、ｃｆＤＮＡ中で一定の割合を占め、がん組織中の細胞に一定の割合で存在すると推定される変異のみを抽出する方法と異なり、がん細胞集団全体に生じている当該がん種特有の変異の傾向を特定することができる。ｃｆＤＮＡにおける変異の評価に関しては、（２−２−２．）に示したようなシーケンスコンテクスト解析により、より詳細な解析が可能になる。本発明の方法は、低侵襲性であり、かつがん細胞特有の変異の傾向の特定や、がんの進行度や悪性度の確認、低進行度の微小腫瘍の検出等に有用である。したがって、本発明の方法のがんの定期検診や健康診断等への適用が考えられる。

（３−３．培養細胞の遺伝情報の評価方法）
本発明の別の好ましい実施形態として、培養細胞の遺伝情報の評価方法が提供される。該方法の具体的手順は、基本的には上述した試験物質の遺伝毒性の評価方法と同じである。ただし試験群としては、変異の有無を調べたい培養細胞の集団を使用する。例えば、試験群としては、ある一定期間継代した細胞であって、その変異の傾向を確認したいものなどが挙げられる。対照群としては、同じ種類の培養細胞であって、遺伝情報既知の（例えば変異の有無及びその変異タイプが確認されている）細胞の集団を使用する。例えば、対照群としては、継代に供する前の細胞などが挙げられる。当該方法によって、培養細胞の変異の有無又はその変異タイプを評価することができる。当該方法は、培養細胞の品質管理のために有用である。

本発明による培養細胞の遺伝情報の評価方法によれば、個別の細胞に生じた変異ではなく、当該培養細胞集団に生じている当該培養細胞特有の変異の傾向を特定することができる。当該方法によれば、ｉＰＳ細胞等の培養細胞において遺伝的な品質が保持されているか否か（変異が起こっていないか否か）を評価することができる。例えば、ヒト由来のｉＰＳ細胞を作成した場合に、その臨床適用をするにあたり遺伝的な品質管理を行うことは極めて重要である。ｉＰＳ細胞のゲノムにおいては、その樹立の過程で、様々な変異を生じることが報告されている。これらは、患者への移植後に発がん等につながる可能性があり、その遺伝的な品質の管理は必須である（Nature, 2011, 471(7336):63-67）。本発明の方法を用いることによって、ｉＰＳ細胞の集団において生じている変異の傾向を簡便にとらえることが可能になる。さらに本発明の方法は、ＰＣＲを用いた従来一般的なｉＰＳ細胞の品質評価方法に比べて網羅性の高い方法であり、かつ別の従来法であるＳＣＩＤマウスを用いた腫瘍形成法（PLoS One, 2012, 7(5):e37342）に比べて非常に安価であり得る。従って、本発明の方法は、ｉＰＳ細胞の遺伝的な品質管理のための、簡便かつ安価なスクリーニング手法としても有用であろう。

（３．４．各種条件）
上記実施形態のいずれにおいても、参照配列としては、試験群の細胞集団のＤＮＡ中の既知配列を使用することができる。該参照配列は、公共のデータベース等に登録されている配列を使用することが好ましいが、上記本発明の方法に先立って予めシーケンサー等で配列決定した該細胞集団のゲノムＤＮＡ中の配列であってもよい。

本発明による試験物質の遺伝毒性の評価方法に使用される試験物質の例としては、その遺伝毒性を評価したい物質であれば特に制限されない。例えば、遺伝毒性を有すると疑われる物質、又は遺伝毒性の有無を確認したい物質、どのような変異を誘発するかを調べたい物質などが挙げられる。試験物質は、天然に存在する物質であっても、化学的もしくは生物学的方法等で人工的に合成した物質であってもよく、又は化合物であっても、組成物もしくは混合物であってもよい。あるいは、該試験物質は、紫外線や放射線などであってもよい。

細胞集団を試験物質に曝露する手段は、試験物質の種類に応じて適宜選択すればよく、特に限定されない。例えば、細胞集団を含む培地に試験物質を添加する方法、細胞集団を試験物質の存在する雰囲気下に置く方法などが挙げられる。

本発明の方法に使用される細胞集団の例としては、動物もしくは植物から採取した検体、及び動物、植物もしくは微生物由来の培養細胞の集団などが挙げられ、好ましくは、動物、植物もしくは微生物の株由来の培養細胞の集団が挙げられる。動物の例としては、好ましくはヒト等の哺乳動物や、カイコ、線虫などが挙げられ、微生物の例としては、好ましくは大腸菌、サルモネラ菌、酵母などが挙げられる。

上記に列挙した細胞集団のうち、試験物質の遺伝毒性の評価方法においては、好ましくは微生物株由来の培養細胞の集団が使用され、より好ましくは大腸菌細胞の集団及びサルモネラ菌細胞の集団からなる群より選択される少なくとも１種が使用される。サルモネラ菌の好ましい例としては、Ｓ．Ｔｙｐｈｉｍｕｒｉｕｍ ＬＴ−２株や、Ａｍｅｓ試験に使用されるＳ．Ｔｙｐｈｉｍｕｒｉｕｍ ＴＡ１００株、ＴＡ９８株、ＴＡ１５３５株、ＴＡ１５３８株、ＴＡ１５３７株等が挙げられる。大腸菌の好ましい例としては、同じくＡｍｅｓ試験に使用されるＷＰ２株、ＷＰ２ ｕｖｒＡ株等が挙げられる。

本発明に使用することができるがん細胞の種類は、特に限定されないが、例えば、肺がん、乳がん、前立腺がん、舌がん、喉頭もしくは咽頭がん、消化器がん（例えば食道がん、胃がん、十二指腸がん、大腸がん、結腸もしくは直腸がん等）、肝臓がん、膵臓がん、子宮頸がん、子宮体がん、腎細胞がん、腎盂がん、膀胱がん、脳腫瘍、骨腫瘍、白血病、リンパ腫、骨髄腫、皮膚がん、悪性黒色腫などが挙げられる。これらのがん細胞は、動物から採取した検体由来のものであってもよく、又は培養がん細胞株であってもよい。

本発明の例示的実施形態として、さらに以下の物質、製造方法、用途、方法等を本明細書に開示する。ただし、本発明はこれらの実施形態に限定されない。

〔１〕試験物質の遺伝毒性の評価方法であって、
（１）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
（２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５）取得した各変異を、塩基対の変異パターンに従って分類すること；
（６）該（５）で得られた変異パターンの各々の変異頻度を決定すること、
を含む、方法。

〔２〕好ましくは、前記（２）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
前記（３）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
〔１〕記載の方法。

〔３〕好ましくは、前記（３）におけるリード配列と参照配列との比較により検出された部位の塩基のうち、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含む、〔１〕記載の方法。

〔４〕好ましくは、前記（３）〜（５）が、
前記リード配列に含まれる塩基を下記(i)〜(iv)に分けること、
(i) 参照配列上の塩基がＡである位置に存在する塩基
(ii) 参照配列上の塩基がＴである位置に存在する塩基
(iii)参照配列上の塩基がＧである位置に存在する塩基
(iv) 参照配列上の塩基がＣである位置に存在する塩基
該リード配列に含まれる塩基の中から、参照配列と塩基がマッチしないものを検出し、当該塩基の存在する部位を、塩基対置換型変異を有する変異部位として取得すること、
検出されたマッチしない塩基の各々について、該変異部位における変異前及び変異後の塩基対を取得すること、及び
該変異部位の塩基対置換型変異を、変異前の塩基対と変異後の塩基対の種類に従って、ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴの６つの塩基対の変異パターンに分類すること、
を含む、〔１〕〜〔３〕のいずれか１項記載の方法。

〔５〕好ましくは、前記６つの塩基対の変異パターンのそれぞれを、オリジナル塩基に従って２群に分類することをさらに含む、〔４〕記載の方法。

〔６〕好ましくは、さらに、
（７）前記試験物質に曝露していない前記細胞集団を対照群とし、前記（１）〜（６）と同様の手順で、該対照群における各塩基対の変異パターンの変異頻度を決定すること；
（８）前記（６）で得られた試験群における各変異パターンの変異頻度から、該（７）で得られた対照群における各変異パターンの変異頻度を引き算すること、
を含む、〔１〕〜〔５〕のいずれか１項記載の方法。

〔７〕試験物質の遺伝毒性の評価方法であって、
（１’）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
（２’）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３’）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４’）該（３’）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５’）取得した変異の各々について、該参照配列に基づいて、変異前の塩基と、該変異前の塩基の上流及び下流に隣接する塩基とを含むコンテクスト配列を決定すること；
（６’）該（４’）で取得した各変異を、該（５’）で決定したコンテクスト配列及び変異後の塩基の種類に従ってタイプ分けすること；
（７’）該（６’）で得られた変異タイプの各々の変異頻度を決定すること、
を含む、方法。

〔８〕好ましくは、前記（２’）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
前記（３’）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
〔７〕記載の方法。

〔９〕好ましくは、前記（３’）におけるリード配列と参照配列との比較により検出された部位の塩基のうち、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含む、〔７〕記載の方法。

〔１０〕好ましくは、前記（３’）〜（６’）が、
前記リード配列に含まれる塩基を下記(i)〜(iv)に分けること、
(i) 参照配列上の塩基がＡである位置に存在する塩基
(ii) 参照配列上の塩基がＴである位置に存在する塩基
(iii)参照配列上の塩基がＧである位置に存在する塩基
(iv) 参照配列上の塩基がＣである位置に存在する塩基
該リード配列に含まれる塩基の中から、参照配列と塩基がマッチしないものを検出し、当該塩基の存在する部位を、塩基対置換型変異を有する変異部位として取得すること、
検出されたマッチしない塩基の各々について、該変異部位における変異前及び変異後の塩基対を取得すること、
該変異部位の塩基対置換型変異を、変異前の塩基対と変異後の塩基対の種類に従って、ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴの６つの塩基対の変異パターンに分類すること、
該変異部位における変異前の塩基と、該変異前の塩基の上流に隣接する１以上の塩基と、該変異前の塩基の下流に隣接する１以上の塩基とからなるコンテクスト配列を決定すること、及び
該６つの塩基対の変異パターンと該コンテクスト配列に従って、該塩基対置換型変異をタイプ分けすること、
を含む、〔７〕〜〔９〕のいずれか１項記載の方法。

〔１１〕好ましくは、前記コンテクスト配列が、変異部位における変異前の塩基と、その両隣の１塩基ずつを含めた３塩基長の配列であり、かつ前記６つの塩基対の変異パターンと該３塩基長のコンテクスト配列に従って、前記塩基対置換型変異が９６にタイプ分けされる、〔１０〕記載の方法。

〔１２〕好ましくは、さらに、
（８’）前記試験物質に曝露していない前記細胞集団を対照群とし、前記（１’）〜（７’）と同じ手順で、対照群における各変異タイプの変異頻度を決定すること；
（９’）前記（７’）で得られた試験群における各変異タイプの変異頻度から、該（８’）で得られた対照群における各変異タイプの変異頻度を引き算すること、
を含む、〔７〕〜〔１１〕のいずれか１項記載の方法。

〔１３〕試験物質の遺伝毒性の評価方法であって、
（１”）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
（２”）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３”）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列における該参照配列に対して塩基が挿入もしくは欠失した部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４”）該（３”）で検出した部位を、挿入もしくは欠失変異を有する変異部位として取得すること；
（５”）取得した変異の各々について、挿入もしくは欠失の塩基長、及び／又は挿入された塩基の種類を決定すること；
（６”）該（５”）で決定された挿入もしくは欠失部位の塩基長及び／又は挿入された塩基の種類ごとの変異頻度を決定すること、
を含む、方法。

〔１４〕好ましくは、前記（２”）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
前記（３”）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
〔１３〕記載の方法。

〔１５〕好ましくは、前記（３”）におけるリード配列と参照配列との比較により検出された部位の塩基のうち、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含む、〔１３〕記載の方法。

〔１６〕検出される挿入もしくは欠失した部位の塩基長が、好ましくは１０ｂｐ以下、より好ましくは１〜５ｂｐである、〔１３〕〜〔１５〕のいずれか１項記載の方法。

〔１７〕好ましくは、さらに、
（７”）前記試験物質に曝露していない前記細胞集団を対照群とし、前記（１”）〜（６”）と同様の手順で、対照群における挿入もしくは欠失部位の塩基長及び／又は挿入された塩基の種類ごとの変異頻度を決定すること；
（８”）前記（６”）で得られた試験群における挿入もしくは欠失部位の塩基長及び／又は挿入された塩基の種類ごとの変異頻度から、該（７”）で得られた対照群における該変異頻度を引き算すること、
を含む、〔１３〕〜〔１６〕のいずれか１項記載の方法。

〔１８〕好ましくは、前記塩基対置換型変異が、１塩基対置換型変異、２塩基対置換型変異、又は３塩基対置換型変異である、〔１〕〜〔１２〕のいずれか１項記載の方法。

〔１９〕好ましくは、前記細胞集団が、サルモネラ菌細胞集団および大腸菌細胞集団からなる群より選択される少なくとも１種である、〔１〕〜〔１８〕のいずれか１項記載の方法。

〔２０〕好ましくは、前記サルモネラ菌がＳ．Ｔｙｐｈｉｍｕｒｉｕｍ ＬＴ−２株、ＴＡ１００株、ＴＡ９８株、ＴＡ１５３５株、ＴＡ１５３８株又はＴＡ１５３７株である、〔１９〕記載の方法。

〔２１〕がん細胞における変異の評価方法であって、
（１）がん細胞集団を試験群とし、そのＤＮＡを取得すること；
（２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５）取得した各変異を、塩基対の変異パターンに従って分類すること；
（６）該（５）で得られた変異パターンの各々の変異頻度を決定すること、
を含む、方法。

〔２２〕培養細胞の遺伝情報の評価方法であって、
（１）培養細胞集団を試験群とし、そのＤＮＡを取得すること；
（２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５）取得した各変異を、塩基対の変異パターンに従って分類すること；
（６）該（５）で得られた変異パターンの各々の変異頻度を決定すること、
を含む、方法。

〔２３〕好ましくは、前記（２）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
前記（４）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
〔２１〕又は〔２２〕記載の方法。

〔２４〕好ましくは、前記（３）におけるリード配列と参照配列との比較により検出された部位の塩基のうち、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含む、〔２１〕又は〔２２〕記載の方法。

〔２５〕好ましくは、前記（３）〜（５）が、
前記リード配列に含まれる塩基を下記(i)〜(iv)に分けること、
(i) 参照配列上の塩基がＡである位置に存在する塩基
(ii) 参照配列上の塩基がＴである位置に存在する塩基
(iii)参照配列上の塩基がＧである位置に存在する塩基
(iv) 参照配列上の塩基がＣである位置に存在する塩基
該リード配列に含まれる塩基の中から、参照配列と塩基がマッチしないものを検出し、当該塩基の存在する部位を、塩基対置換型変異を有する変異部位として取得すること、
検出されたマッチしない塩基の各々について、該変異部位における変異前及び変異後の塩基対を取得すること、及び
該変異部位の塩基対置換型変異を、変異前の塩基対と変異後の塩基対の種類に従って、ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴの６つの塩基対の変異パターンに分類すること、
を含む、〔２１〕〜〔２４〕のいずれか１項記載の方法。

〔２６〕好ましくは、前記６つの塩基対の変異パターンのそれぞれを、オリジナル塩基に従って２群に分類することをさらに含む、〔２５〕記載の方法。

〔２７〕好ましくは、さらに、
（７）前記（１）〜（６）と同様の手順で、対照群における各塩基対の変異パターンの変異頻度を決定すること；
（８）前記（６）で得られた試験群における各変異パターンの変異頻度から、該（７）で得られた対照群における各変異パターンの変異頻度を引き算すること、
を含む、〔２１〕〜〔２６〕のいずれか１項記載の方法。

〔２８〕好ましくは、
前記方法が、がん細胞における変異の評価方法であり、前記対照群が、非がん細胞集団であるか、又は、
前記試験群が、がんが疑われる細胞又はがん化リスクを評価したい細胞の集団であり、前記対照群が、非がん細胞集団又はがん化リスクの低い細胞の集団であり、かつ前記方法によって細胞のがん化リスクが評価される、
〔２７〕記載の方法。

〔２９〕好ましくは、前記方法が、培養細胞の遺伝情報の評価方法であり、前記対照群が、前記試験群と同じ種類の培養細胞であって、遺伝情報既知の細胞の集団である、〔２７〕記載の方法。

〔３０〕がん細胞における変異の評価方法であって、
（１’）がん細胞集団を試験群とし、そのＤＮＡを取得すること；
（２’）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３’）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４’）該（３’）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５’）取得した変異の各々について、該参照配列に基づいて、変異前の塩基と、該変異前の塩基の上流及び下流に隣接する塩基とを含むコンテクスト配列を決定すること；
（６’）該（４’）で取得した各変異を、該（５’）で決定したコンテクスト配列及び変異後の塩基の種類に従ってタイプ分けすること；
（７’）該（６’）で得られた変異タイプの各々の変異頻度を決定すること、
を含む、方法。

〔３１〕培養細胞の遺伝情報の評価方法であって、
（１’）培養細胞集団を試験群とし、そのＤＮＡを取得すること；
（２’）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３’）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４’）該（３’）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
（５’）取得した変異の各々について、該参照配列に基づいて、変異前の塩基と、該変異前の塩基の上流及び下流に隣接する塩基とを含むコンテクスト配列を決定すること；
（６’）該（４’）で取得した各変異を、該（５’）で決定したコンテクスト配列及び変異後の塩基の種類に従ってタイプ分けすること；
（７’）該（６’）で得られた変異タイプの各々の変異頻度を決定すること、
を含む、方法。

〔３２〕好ましくは、前記（２’）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
前記（３’）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
〔３０〕又は〔３１〕記載の方法。

〔３３〕好ましくは、前記（３’）におけるリード配列と参照配列との比較により検出された部位の塩基のうち、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含む、〔３０〕又は〔３１〕記載の方法。

〔３４〕好ましくは、前記（３’）〜（６’）が、
前記リード配列に含まれる塩基を下記(i)〜(iv)に分けること、
(i) 参照配列上の塩基がＡである位置に存在する塩基
(ii) 参照配列上の塩基がＴである位置に存在する塩基
(iii)参照配列上の塩基がＧである位置に存在する塩基
(iv) 参照配列上の塩基がＣである位置に存在する塩基
該リード配列に含まれる塩基の中から、参照配列と塩基がマッチしないものを検出し、当該塩基の存在する部位を、塩基対置換型変異を有する変異部位として取得すること、
検出されたマッチしない塩基の各々について、該変異部位における変異前及び変異後の塩基対を取得すること、
該変異部位の塩基対置換型変異を、変異前の塩基対と変異後の塩基対の種類に従って、ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴの６つの塩基対の変異パターンに分類すること、
該変異部位における変異前の塩基と、該変異前の塩基の上流に隣接する１以上の塩基と、該変異前の塩基の下流に隣接する１以上の塩基とからなるコンテクスト配列を決定すること、及び
該６つの塩基対の変異パターンと該コンテクスト配列に従って、該塩基対置換型変異をタイプ分けすること、
を含む、〔３０〕〜〔３３〕のいずれか１項記載の方法。

〔３５〕好ましくは、前記コンテクスト配列が、変異部位における変異前の塩基と、その両隣の１塩基ずつを含めた３塩基長の配列であり、かつ前記６つの塩基対の変異パターンと該３塩基長のコンテクスト配列に従って、前記塩基対置換型変異が９６にタイプ分けされる、〔３４〕記載の方法。

〔３６〕好ましくは、さらに、
（８’）前記（１’）〜（７’）と同じ手順で、対照群における各変異タイプの変異頻度を決定すること；
（９’）前記（７’）で得られた試験群における各変異タイプの変異頻度から、該（８’）で得られた対照群における各変異タイプの変異頻度を引き算すること、
を含む、〔３０〕〜〔３５〕のいずれか１項記載の方法。

〔３７〕好ましくは、
前記方法が、がん細胞における変異の評価方法であり、前記対照群が、非がん細胞集団であるか、又は、
前記試験群が、がんが疑われる細胞又はがん化リスクを評価したい細胞の集団であり、前記対照群が、非がん細胞集団又はがん化リスクの低い細胞の集団であり、かつ前記方法によって細胞のがん化リスクが評価される、
〔３６〕記載の方法。

〔３８〕好ましくは、前記方法が、培養細胞の遺伝情報の評価方法であり、前記対照群が、前記試験群と同じ種類の培養細胞であって、遺伝情報既知の細胞の集団である、〔３６〕記載の方法。

〔３９〕がん細胞における変異の評価方法であって、
（１”）がん細胞集団を試験群とし、そのＤＮＡを取得すること；
（２”）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３”）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列における該参照配列に対して塩基が挿入もしくは欠失した部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４”）該（３”）で検出した部位を、挿入もしくは欠失変異を有する変異部位として取得すること；
（５”）取得した変異の各々について、挿入もしくは欠失の塩基長、及び／又は挿入された塩基の種類を決定すること；
（６”）該（５”）で決定された挿入もしくは欠失部位の塩基長及び／又は挿入された塩基の種類ごとの変異頻度を決定すること、
を含む、方法。

〔４０〕培養細胞の遺伝情報の評価方法であって、
（１”）培養細胞集団を試験群とし、そのＤＮＡを取得すること；
（２”）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
（３”）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列における該参照配列に対して塩基が挿入もしくは欠失した部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
（４”）該（３”）で検出した部位を、挿入もしくは欠失変異を有する変異部位として取得すること；
（５”）取得した変異の各々について、挿入もしくは欠失の塩基長、及び／又は挿入された塩基の種類を決定すること；
（６”）該（５”）で決定された挿入もしくは欠失部位の塩基長及び／又は挿入された塩基の種類ごとの変異頻度を決定すること、
を含む、方法。

〔４１〕好ましくは、前記（２”）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
前記（３”）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
〔３９〕又は〔４０〕記載の方法。

〔４２〕好ましくは、前記（３”）におけるリード配列と参照配列との比較により検出された部位の塩基のうち、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含む、〔３９〕又は〔４０〕記載の方法。

〔４３〕検出される挿入もしくは欠失した部位の塩基長が、好ましくは１０ｂｐ以下、より好ましくは１〜５ｂｐである、〔３９〕〜〔４２〕のいずれか１項記載の方法。

〔４４〕好ましくは、さらに、
（７”）前記（１”）〜（６”）と同様の手順で、対照群における挿入もしくは欠失部位の塩基長及び／又は挿入された塩基の種類ごとの変異頻度を決定すること；
（８”）前記（６”）で得られた試験群における挿入もしくは欠失部位の塩基長及び／又は挿入された塩基の種類ごとの変異頻度から、該（７”）で得られた対照群における該変異頻度を引き算すること、
を含む、〔３９〕〜〔４３〕のいずれか１項記載の方法。

〔４５〕好ましくは、
前記方法が、がん細胞における変異の評価方法であり、前記対照群が、非がん細胞集団であるか、又は、
前記試験群が、がんが疑われる細胞又はがん化リスクを評価したい細胞の集団であり、前記対照群が、非がん細胞集団又はがん化リスクの低い細胞の集団であり、かつ前記方法によって細胞のがん化リスクが評価される、
〔４４〕記載の方法。

〔４６〕好ましくは、前記方法が、培養細胞の遺伝情報の評価方法であり、前記対照群が、前記試験群と同じ種類の培養細胞であって、遺伝情報既知の細胞の集団である、〔４４〕記載の方法。

〔４７〕好ましくは、前記塩基対置換変異が、１塩基対置換型変異、２塩基対置換型変異、又は３塩基対置換型変異である、〔２１〕〜〔３８〕のいずれか１項記載の方法。

〔４８〕前記検出に用いるリード配列の総量が、好ましくは１×１０¹⁰ｂｐ以下、より好ましくは１×１０⁹ｂｐ以下、さらに好ましくは１×１０⁸ｂｐ以下、よりさらに好ましくは１×１０⁷ｂｐ以下、なお好ましくは１×１０⁶ｂｐ以下である、〔１〕〜〔４７〕のいずれか１項記載の方法。

以下、実施例を示し、本発明をより具体的に説明する。

実施例１ 解析方法のバリデーション
本実施例では、変異頻度既知の合成ＤＮＡサンプルを本発明の試験方法で解析して変異を定性及び定量評価することによって、本発明の解析方法の有効性を検証した。

１．ＤＮＡサンプルの調製
様々な変異パターンを既知の量で含む合成ＤＮＡサンプルを調製した。模式図１に、ＤＮＡサンプル調製手順の概念図を示す。１０００ｂｐのランダムな配列を有する合成ＤＮＡ配列（配列番号１；以下、ランダムＤＮＡ配列と称する）を合成した。このランダムＤＮＡ配列には、ＧＣ、ＡＴ塩基対が約５０％ずつ存在した。このランダムＤＮＡ配列をベースに、変異（塩基対置換型変異又は短い挿入・欠失変異）を導入したＤＮＡ配列（以下、変異ＤＮＡ配列とする）を調製した。以下に詳細を説明する。

塩基対置換型変異を含む変異ＤＮＡ配列については、ランダムＤＮＡ配列の中心に位置するＧＣ塩基対（５０１番目）を他の塩基対（ＧＣ→ＴＡ、ＣＧ又はＡＴ、表１参照）で置換した３種類の変異配列を作製し、それぞれｐＴＡＫＮ−２ベクターに組み込んだ。得られたベクターをＴＥバッファー（ｐＨ８．０、和光純薬工業社製）に溶解し、１００ｎｇ／μＬの濃度に調整し、各変異ＤＮＡ配列を含む溶液を等量混合した。同様に、ＡＴ塩基対（５０２番目）を他の塩基対（ＡＴ→ＴＡ、ＣＧ又はＧＣ、表１参照）で置換した３種類の変異配列を作製し、それぞれｐＴＡＫＮ−２ベクターに組み込み、ベクターのＴＥバッファー溶液（１００ｎｇ／μＬ）を調製し、得られた３種類の溶液を等量混合した。各々の等量混合溶液を変異ＤＮＡ溶液とした。また、該変異ＤＮＡ溶液１０μＬを９０μＬのＴＥバッファーと混合して、１０倍希釈変異ＤＮＡ溶液を調製した。さらに該１０倍希釈変異ＤＮＡ溶液１０μＬを９０μＬのＴＥバッファーを混合して、１００倍希釈変異ＤＮＡ溶液を調製した。別途、ランダムＤＮＡ配列をｐＴＡＫＮ−２ベクターに組み込み、ベクターのＴＥバッファー溶液（１００ｎｇ／μＬ）を調製した（ランダムＤＮＡ溶液）。ランダムＤＮＡ溶液に、変異ＤＮＡ溶液、１０倍希釈変異ＤＮＡ溶液、又は１００倍希釈変異ＤＮＡ溶液を混合し、各塩基対置換が均等の頻度で認められ、かつ総変異頻度が１／１０³、１／１０⁴、１／１０⁵、１／１０⁶ｂｐのＤＮＡサンプルを調製した（表２参照）。

短い挿入・欠失変異を含む変異ＤＮＡ配列については、５０１番塩基対の前に一塩基（Ａ、すなわちＡＴ塩基対）を挿入した変異配列を作製した（表１参照）。これを上記と同様にｐＴＡＫＮ−２ベクターに組み込み、ベクターのＴＥバッファー溶液（１００ｎｇ／μＬ）を調製して変異ＤＮＡ溶液とした。また、該変異サンプル溶液の１０倍及び１００倍希釈変異ＤＮＡ溶液を調製した。ランダムＤＮＡ溶液に、変異ＤＮＡ溶液、１０倍希釈変異ＤＮＡ溶液、又は１００倍希釈変異ＤＮＡ溶液を混合し、総変異頻度が１／１０³、１／１０⁴、及び１／１０⁵、１／１０⁶ｂｐのＤＮＡサンプルを調製した（表２参照）。各変異ＤＮＡ溶液を含むＤＮＡサンプルを変異サンプル、変異ＤＮＡ溶液を含まない（ランダムＤＮＡ溶液のみの）ＤＮＡサンプルを対照サンプルとして、以下のシーケンシングを行った。

２．高スループットシーケンシング
１．で調製した各変異サンプル及び対照サンプルについて、次世代シーケンサーＨｉＳｅｑ２５００（イルミナ社製、以下、ＨｉＳｅｑとも称する）を用いて標準プロトコールに従って、塩基配列を解読した。その際、ＤＮＡは超音波処理により平均約１５０ｂｐの長さにフラグメント化し、各フラグメント両端にアダプターを付加して、２×１２５ｂｐのリード長でシーケンシングした。サンプルあたり、平均で１．９Ｇｂｐの塩基配列情報を得た。

３．リード配列の編集及び変異解析用フォーマットの作製
シーケンシングによって得られたリード配列の編集及び解析フローの概念図を模式図２に示す。まず、Ｃｕｔａｄａｐｔソフトウェア（Ｍａｒｔｉｎ，２０１１）を用い、各リード末端のアダプター配列の除去及び、クオリティの低い塩基の除去を行った。
i）ＨｉＳｅｑでは、シーケンシングに供されたＤＮＡフラグメント（オリジナルフラグメント）一つにつき、片側からリード１配列を読み取った後、反対側からそのペアとなるリード２配列を取得する。そこで、Ｃｕｔａｄａｐｔの選定をパスした各ペアのリード配列について、ＰＥＡＲソフトウェアを用い、両リードの配列が対を形成する部分を重ね合わせ、一本の合成リード（Ｃｏｎｊｕｇａｔｅｄ Ｒｅａｄ）を構築した。合成リードの配列の塩基は、全てリード１にあわせた。合成リードにおける各塩基のクオリティ値は、リード１とリード２の塩基が相補的である場合には両塩基のクオリティ値の和、両リードの塩基が相補的でない場合にはクオリティ値の大きい方の値から小さい方のクオリティ値を引いた値を採用した。これにより、クオリティ値の違いにしたがって、合成リードにおける各塩基のうち、リード１と２の塩基が対形成していたものを選別することができる。
ii）作製した各オリジナルフラグメントについての合成リードを、Ｂｏｗｔｉｅ２ソフトウェアを用いて、参照配列（ランダムＤＮＡ配列を挿入したｐＴＡＫＮ−２ベクターの配列）にマッピングし、Ｓａｍフォーマットのファイルを作成した。
iii）得られたＳａｍファイルを、Ｓａｍｔｏｏｌｓソフトウェアを用いてｐｉｌｅｕｐフォーマットに変換した。この際、塩基のクオリティ値をもとに、解析対象の塩基情報を、ペアのリードの重なり領域において両リードの塩基が相補的に対形成していた範囲に限定した。
iv）得られたｐｉｌｅｕｐフォーマットを、プログラミング言語Ｐｙｔｈｏｎを用いて作成したプログラムを用いた変異解析に供した。

４．変異解析
１）塩基対置換型変異の検出
１塩基対置換型変異についての変異解析アルゴリズムの概念図を模式図３に示す。解析に供したｐｉｌｅｕｐフォーマットから、プログラミング言語Ｐｙｔｈｏｎを用いて作成したプログラムを用いて、リード配列中の全解析対象塩基を、対応する参照配列の塩基がＡである群、Ｔである群、Ｇである群、及びＣである群の４群に分類した。次いで、各群に振り分けられた塩基の総数、及び変異した塩基を検出した。得られたデータから、変異前のＡＴ塩基対１０⁶ｂｐあたりに占めるＡＴ塩基対の各変異パターン（ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ）の変異コール割合、及び変異前のＧＣ塩基対１０⁶ｂｐあたりに占めるＧＣ塩基対の各変異パターン（ＧＣ→ＴＡ、ＧＣ→ＣＧ、ＧＣ→ＡＴ）の変異コール割合を算出した。

塩基対置換型変異を含む変異サンプルにおける、ＧＣ及びＡＴ塩基対それぞれの各変異パターンの変異コール割合を図１に示す。いずれの変異パターンについても、サンプル中の変異頻度に依存して変異コール割合が上昇した。対照サンプルでも変異が検出されたが、これは背景エラー（シーケンシングエラーを含むサンプル調製からシーケンシングの過程で生じたエラー）を表す。さらに、対照サンプルでも変異コール割合は変異パターンごとに異なっており、ＧＣ塩基対は、ＡＴ塩基対に比べて変異コール割合が高い傾向にあった。これは、ＧＣ塩基対がＤＮＡ抽出等のライブラリ作製過程で、酸化等の化学修飾の影響を受けやすいことが原因と考えられる。

２）変異頻度増加量の算出
次に、変異サンプルの変異コール割合から、対照サンプルの変異コール割合を差し引くことによって、シーケンスエラーを含む背景エラーを除外し、対照サンプルに対する変異サンプルでの変異頻度の増加量（以下、変異頻度増加量という）を算出した。算出した変異頻度増加量を図２に示す。いずれの変異パターンにおいても、変異頻度増加量と、導入した変異頻度は概ね一致しており、本方法により約１０⁵ｂｐに一つの頻度の塩基対置換型変異を検出できたことが示された。

３）短い挿入・欠失変異の検出
短い挿入・欠失変異については、プログラミング言語Ｐｙｔｈｏｎを用いて作成したプログラムを用いて、ランダムＤＮＡ配列に対して１０ｂｐ以下の長さで挿入又は欠失した塩基を全て検出し、その挿入又は欠失長さ（ｂｐ）、及び挿入又は欠失塩基の種類ごとに出現頻度を計数した。さらに、前述した塩基対置換型変異の解析と同様に、対照サンプルの変異頻度を差し引いて背景エラーを除外し、変異頻度増加量を算出した。挿入変異の変異頻度増加量の結果を図３に示す。図３Ａは挿入された塩基の長さ（ｂｐ）、図３Ｂは挿入された塩基の種類について、変異頻度増加量を示している。本方法により、約１０⁵ｂｐに一つの頻度の短い挿入変異を検出することができた。また、挿入変異と同様の方法により一塩基欠失についても検討し、約１０⁵ｂｐに一つの頻度の短い欠失変異を検出することができた（データは示さない）。

５．結論
本実施例では、様々な変異情報を含むＤＮＡにおける塩基対置換型変異、及び短い挿入・欠失変異について、量（頻度）及び質（変異パターン）を含む総合的な変異情報を高感度に取得することができた。これらの結果から、本発明の解析方法によりＤＮＡ中に存在する低頻度の変異を定性的かつ定量的に解析できることが示された。

実施例２ 変異原による遺伝毒性の解析
本実施例では、本発明の解析方法による、変異原への曝露によって生物のゲノムに生じた変異パターンの定性及び定量解析により、該変異原の遺伝毒性を解析した。変異原としてＥｔｈｙｌｎｉｔｒｏｓｏｕｒｅａ（ＥＮＵ、ＣＡＳ Ｎｏ．７５９−７３−９）を用いた。変異原に曝露する生物として、塩基対置換型の変異を検出可能な、Ａｍｅｓ試験に汎用されるサルモネラ菌ＴＡ１００株を使用した。実験は、独立した３回の操作を行い、サンプルを調製した（ｎ＝３）。

１．ＴＡ１００株の変異原への曝露
変異原への曝露は、Ａｍｅｓ試験のプレインキュベーション法（K. Mortelmans et al., Mutat. Res. - Fundam. Mol. Mech. Mutagen., 2000, 455, 29-60）に準拠して実施した。ＴＡ１００株を２ｍＬのニュートリエントブイヨン Ｎｏ．２（Ｏｘｏｉｄ社製）に植菌し、３７℃、１８０ｒｐｍで４時間振とう培養し、Ｏ．Ｄ．６６０値が１．０以上の前培養液を得た。ＥＮＵ（５４％；シグマアルドリッチ社製）は、ジメチルスルホキシド（ＤＭＳＯ；和光純薬工業製）で希釈した。試験管内に、適切な濃度に希釈したＥＮＵ溶液１００μＬ、０．１Ｍリン酸バッファー５００μＬ、及び前培養液１００μＬを添加し（ＥＮＵ濃度：６７．５、１３５、２７０、４０５、５４０、８１０、及び１０８０μｇ／ｔｕｂｅ、３７℃のウォーターバス中で２０分間、１００ｒｐｍで振とう培養した。対照群としては、ＥＮＵ溶液の代わりに溶媒（ＤＭＳＯ）１００μＬを添加した。２０分間振とう培養後、培養液を含む試験管をウォーターバスから取り出し、予め分注しておいた２ｍＬのＮｕｔｒｉｅｎｔ Ｂｒｏｔｈ溶液に培養液５０μＬを添加し、３７℃、１８０ｒｐｍで１４時間追培養した後、菌懸濁液を１ｍＬ回収し、７５００ｒｐｍで５分間遠心し、上清を除去し、菌体を回収した。

また、Ａｍｅｓ試験用に、上記と同様の条件でＥＮＵを曝露した菌懸濁液を作製し、４５℃に加温した２ｍＬのｔｏｐ ａｇａｒ（１％ＮａＣｌ、１％ａｇａｒ、０．０５ｍＭ Ｈｉｓｔｉｄｉｎｅ及び０．０５ｍＭ Ｂｉｏｔｉｎｅを含む）を添加し、ボルテックスで懸濁した後、最小グルコース寒天培地（テスメディア（登録商標）ＡＮ；オリエンタル酵母工業製）上に重層した。得られたプレートを３７℃で４８時間培養後、観察されたコロニーを計数した。

２．Ｔｏｔａｌ ＤＮＡの回収
１．で得られた菌体から、ＤＮｅａｓｙ Ｂｌｏｏｄ ＆ Ｔｉｓｓｕｅ Ｋｉｔ（キアゲン社製）を用い、推奨プロトコールに従って、Ｔｏｔａｌ ＤＮＡを回収した。

３．高スループットシーケンス
２．で回収した対照群及びＥＮＵ処理群からのＴｏｔａｌ ＤＮＡ液を用いて、次世代シーケンサーＨｉＳｅｑ２５００（イルミナ社製）により、標準プロトコールに従って塩基配列を解読した。その際、ＤＮＡは超音波処理により平均約１５０ｂｐの長さにフラグメント化し、各フラグメント両端にアダプターを付加して、２×１２５ｂｐのリード長でシーケンシングした。サンプルあたり、平均で５．０Ｇｂｐの塩基配列情報を得た。

４．リード配列の編集及び変異解析
シーケンシングによって得られたリード配列の編集及び変異解析は、実施例１と同様に、模式図２に示す解析フローの概念図に従って、Ｃｕｔａｄａｐｔソフトウェア、ＰＥＡＲソフトウェア、Ｂｏｗｔｉｅ２ソフトウェア、Ｓａｍｔｏｏｌｓソフトウェア、及びプログラミング言語Ｐｙｔｈｏｎを用いて作成したプログラムを用いて行った。本実施例においては、Ｂｏｗｔｉｅ２ソフトウェアでのマッピング後、Ｐｉｃａｒｄ ｔｏｏｌｓ（broadinstitute.github.io/picard/）を用いて、ＰＣＲ ｄｕｐｌｉｃａｔｅｓの除去を行った。

Ｂｏｗｔｉｅ２ソフトウェアでマッピングする参照配列は、Ｓ．Ｔｙｐｈｉｍｕｒｉｕｍ ＴＡ１００株のゲノム配列をもとに構築した。まず、ＴＡ１００株からＤＮＡを抽出し、次世代シーケンサーＨｉＳｅｑ２５００（イルミナ社製）により、標準プロトコールに従って塩基配列を解読した。その際、ＤＮＡは超音波処理により平均約３００ｂｐの長さにフラグメント化し、各フラグメント両端にアダプターを付加して、２×１２５ｂｐのリード長でシーケンシングした。得られたリード配列をＳ．Ｔｙｐｈｉｍｕｒｉｕｍ ＬＴ−２株のゲノム配列とｐＳＬＴプラスミド配列（ＧＣＡ０００００６９４５．２）、およびＲ４６プラスミド配列（ＮＣ＿００３２９２．１）にマッピング後、Ｓａｍｔｏｏｌｓソフトウェアを用いて変異検出を行った。得られた変異情報を反映したＴＡ１００株のゲノム配列に基づく参照配列を作製し、変異解析に使用した。ＴＡ１００株のゲノム配列を配列番号２に示す。

５．変異頻度増加量の算出
実施例１と同様の手順で、プログラミング言語Ｐｙｔｈｏｎを用いたプログラムを用いて、参照配列に対してマッピングされた全リード配列中の全解析対象塩基を、対応する参照配列の塩基によって４群に振り分け、次いで各群の塩基の総数、及び参照配列に対して変異した塩基を検出した。変異した塩基を参照配列の塩基と比較することで、ＥＮＵ処理群及び対照群それぞれについて、解析対象塩基中におけるＡＴ塩基対、ＧＣ塩基対の各１０⁶ｂｐにおける各変異パターン（ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、及びＧＣ→ＴＡ、ＧＣ→ＣＧ、ＧＣ→ＡＴ）、及び各変異パターンの変異コール割合を算出した。対照群とＥＮＵ処理群間における各変異パターンの頻度についての統計学的検定は、変異パターンごとにＤｕｎｎｅｔｔの多重比較検定にて行った。

次いで、実施例１の２）と同様の手順で各塩基対の変異パターンごとに、ＥＮＵ処理群の変異コール割合から対照群の変異コール割合を差し引くことにより、ＥＮＵ曝露による変異頻度増加量を算出した。

６．シーケンスコンテクスト解析
本検討では、シーケンスコンテクストに基づいて変異コール割合を解析した。すなわち、５．で得られた変異コール情報及び参照配列の情報をもとに、各々の変異コール箇所の塩基対の変異タイプを検出した。さらに参照配列の情報をもとに、各変異コール箇所及びその両隣の塩基を含む３塩基の情報を収集した。各変異コール箇所の変異を、６通りの塩基対変異タイプに、その両隣の塩基情報（４×４＝１６通り）をかけ合わせた９６通りのタイプに分類した。このシーケンスコンテクストに基づく９６通りの変異タイプについて、各々の変異コール割合（／１０⁶ｂｐ）を算出した。各変異タイプについて、ＥＮＵ処理群における変異コール割合から対照群の変異コール割合を差し引くことにより、ＥＮＵ曝露による変異頻度増加量を算出した。

７．結果
１）Ａｍｅｓ試験の復帰突然変異体数
表３にＥＮＵ曝露後の復帰突然変異コロニー数を示す。データは３枚のプレートの測定値の平均及び標準偏差を示す。ＥＮＵ曝露により、復帰変異突然変異体数の増加が認められたことから、ＥＮＵ曝露によりＴＡ１００株のゲノム中に変異が導入されたことが示された。

２）ＥＮＵ曝露による変異コール割合の変化の算出
ＥＮＵ処理群（ＥＮＵ濃度１３５、２７０、４０５及び５４０μｇ／ｔｕｂｅ）について、５．で算出した変異頻度増加量を図４に示す。ＥＮＵ曝露によって、複数の塩基対の変異パターンの頻度に増加が認められた。ＧＣ塩基対においては、ＧＣ＞ＡＴ変異の頻度の増加が認められた（図４Ａ）。一方、ＡＴ塩基対においては、主にＡＴ＞ＴＡ、及びＡＴ＞ＧＣ変異の頻度の増加が認められた（図４Ｂ）。

３）各変異パターンのオリジナル塩基による分類
ＨｉＳｅｑにおいては、リード１の配列がシーケンシング反応に供された元のＤＮＡ断片（オリジナルフラグメント）に対応している。したがって、リード１配列の塩基が参照配列のＡ、Ｔ、Ｇ及びＣのいずれの塩基にマップされたかを調べることにより、オリジナルフラグメントの変異箇所における変異前の塩基（すなわち、オリジナル塩基）を確認できる。
背景エラー頻度はオリジナル塩基によって異なり得るため、変異パターンをオリジナル塩基に従ってさらに分類し、各分類の変異頻度増加量を求めた。すなわち、上記５．で求めた対照サンプル及び変異サンプルにおける６つの塩基対変異パターンのそれぞれを、オリジナル塩基の種類によってさらに２つずつに分類し、各分類の変異頻度を求めた。次いで、変異サンプルの変異頻度から対応する対照サンプルの変異頻度を差し引いて、変異頻度増加量を算出した。

図５はＥＮＵ暴露サンプルにおけるオリジナル塩基の種類ごとに分けた塩基対の変異パターンの変異頻度増加量を示している。通常、ＥＮＵの曝露によってゲノム中に固定された変異の場合、塩基対の双方の塩基が変化するため、塩基対を形成するどちらの塩基においても同等の頻度で変異頻度の増加が認められる。しかし、図５中のＧＣ＞ＴＡ変異においては、明らかにオリジナル塩基がＧの場合にオリジナル塩基がＣの場合と比べて変異頻度が高い傾向が認められた。

上記のバイアスは、ＧＣ＞ＴＡ変異の大部分が、化学修飾されたＧによって引き起こされたことを強く示唆している。すなわち、本検討で認められたＧＣ＞ＴＡ変異の増加は、ゲノム内に固定された変異ではなく、シーケンシングにおける化学修飾されたＧ塩基の読み取りエラーを反映しているものと考えられた。Ｇ塩基は、ＤＮＡサンプルの調製過程で酸化による化学修飾を受けやすく、シーケンシングにおけるＧをＴとして読み間違えるエラーを起こしやすいことが知られている。したがって、本検討においていくつかのサンプルで認められたＧＣ＞ＴＡ変異頻度の増加は、ＤＮＡ調製過程にＧが酸化されたことにより生じた人為的影響であると考えられた。よって本解析方法により、細胞からのＤＮＡサンプル調製過程で生じるエラーを原因としたシーケンシングエラーを取り除くことができることが示唆された。

４）変異スペクトル解析
次に、ＥＮＵ５４０μｇ／ｔｕｂｅ群で増加の認められた各変異パターンの変異頻度増加量をもとに、変異スペクトルの解析を行った。具体的には、各変異パターンの１０⁶ｂｐにおける変異頻度増加量を合計し、各変異パターンの全体に占める割合を算出した。結果を図６及び表４に示す。最も多く頻度増加の認められた変異パターンはＧＣ＞ＡＴ変異であり、次に割合の高かった変異パターンはＡＴ＞ＧＣ変異であった。これは非特許文献５において、Ａｍｅｓ試験菌株の一つであるＹＧ７１０８株において認められたＥＮＵの変異パターンと同様の結果であったことから、本方法により、ＥＮＵによる変異パターンを正確に検出できたことが示唆された。

５）シーケンスコンテクスト解析
シーケンスコンテクスト解析で算出された変異頻度増加量を図７及び８に示す。図中の変異パターンの表記は、変異した塩基対のうちのピリミジン塩基の変異パターン（Ｃ＞Ａ、Ｃ＞Ｇ、Ｃ＞Ｔ、Ｔ＞Ａ、Ｔ＞Ｃ及びＴ＞Ｇ）と、該ピリミジン塩基とその両隣の塩基を含む３塩基配列（例えば、Ｃ＞Ｔにおいて該Ｃ塩基がＡとＴに挟まれていた場合は、ＡＣＴと表記される）で表されている。本解析において最も変異頻度増加量の大きかったＣ＞Ｔ変異においては、Ｃの３’塩基側にピリミジン塩基（Ｃ又はＴ）の位置するコンテクストにおいて、変異頻度増加量が高くなる傾向が認められた。これは、アルキル化剤による変異シグニチャーの示したパターン（図７〜８中のＳｉｇｎｉｔｕｒｅ １１、Nature, 2013, 500(7463):415-421参照）と類似しており、ＥＮＵがアルキル化剤であることと矛盾しない結果となった。この結果から、シーケンスコンテクスト解析は、変異原の遺伝毒性メカニズムの推定や、ヒトの発がんにおける役割を簡便に推定することができる有用な方法であると考えられた。

Claims (28)

1. 試験物質の遺伝毒性の評価方法であって、
   （１）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
   （２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
   （３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
   （４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
   （５）取得した各変異を、塩基対の変異パターンに従って分類すること；
   （６）該（５）で得られた変異パターンの各々の変異頻度を決定すること；
   （７）該試験物質に曝露していない該細胞集団を対照群とし、該（１）〜（６）と同様の手順で、該対照群における各塩基対の変異パターンの変異頻度を決定すること；
   （８）該（６）で得られた試験群における各変異パターンの変異頻度から、該（７）で得られた対照群における各変異パターンの変異頻度を引き算すること、
   を含む、方法。
2. 前記（２）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
   前記（３）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
   請求項１記載の方法。
3. 前記（３）〜（５）が、
   前記リード配列に含まれる塩基を下記(i)〜(iv)に分けること、
   (i) 参照配列上の塩基がＡである位置に存在する塩基
   (ii) 参照配列上の塩基がＴである位置に存在する塩基
   (iii)参照配列上の塩基がＧである位置に存在する塩基
   (iv) 参照配列上の塩基がＣである位置に存在する塩基
   該リード配列に含まれる塩基の中から、参照配列と塩基がマッチしないものを検出し、当該塩基の存在する部位を、塩基対置換型変異を有する変異部位として取得すること、
   検出されたマッチしない塩基の各々について、該変異部位における変異前及び変異後の塩基対を取得すること、及び
   該変異部位の塩基対置換型変異を、変異前の塩基対と変異後の塩基対の種類に従って、ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴの６つの塩基対の変異パターンに分類すること、
   を含む、請求項１又は２記載の方法。
4. 試験物質の遺伝毒性の評価方法であって、
   （１'）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
   （２'）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
   （３'）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
   （４'）該（３'）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
   （５'）取得した変異の各々について、該参照配列に基づいて、変異前の塩基と、該変異前の塩基の上流及び下流に隣接する塩基とを含むコンテクスト配列を決定すること；
   （６'）該（４'）で取得した各変異を、該（５'）で決定したコンテクスト配列及び変異後の塩基の種類に従ってタイプ分けすること；
   （７'）該（６'）で得られた変異タイプの各々の変異頻度を決定すること；
   （８'）該試験物質に曝露していない該細胞集団を対照群とし、該（１'）〜（７'）と同じ手順で、対照群における各変異タイプの変異頻度を決定すること；
   （９'）該（７'）で得られた試験群における各変異タイプの変異頻度から、該（８'）で得られた対照群における各変異タイプの変異頻度を引き算すること、
   を含む、方法。
5. 前記（２'）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
   前記（３'）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
   請求項４記載の方法。
6. 前記（３'）〜（６'）が、
   前記リード配列に含まれる塩基を下記(i)〜(iv)に分けること、
   (i) 参照配列上の塩基がＡである位置に存在する塩基
   (ii) 参照配列上の塩基がＴである位置に存在する塩基
   (iii)参照配列上の塩基がＧである位置に存在する塩基
   (iv) 参照配列上の塩基がＣである位置に存在する塩基
   該リード配列に含まれる塩基の中から、参照配列と塩基がマッチしないものを検出し、当該塩基の存在する部位を、塩基対置換型変異を有する変異部位として取得すること、
   検出されたマッチしない塩基の各々について、該変異部位における変異前及び変異後の塩基対を取得すること、
   該変異部位の塩基対置換型変異を、変異前の塩基対と変異後の塩基対の種類に従って、ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴの６つの塩基対の変異パターンに分類すること、
   該変異部位における変異前の塩基と、該変異前の塩基の上流に隣接する１以上の塩基と、該変異前の塩基の下流に隣接する１以上の塩基とからなるコンテクスト配列を決定すること、及び
   該６つの塩基対の変異パターンと該コンテクスト配列に従って、該塩基対置換型変異をタイプ分けすること、
   を含む、請求項４又は５記載の方法。
7. 試験物質の遺伝毒性の評価方法であって、
   （１"）試験物質に曝露した細胞集団を試験群とし、そのＤＮＡを取得すること；
   （２"）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
   （３"）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列における該参照配列に対して塩基が挿入もしくは欠失した部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
   （４"）該（３"）で検出した部位を、挿入もしくは欠失変異を有する変異部位として取得すること；
   （５"）取得した変異の各々について、挿入もしくは欠失の塩基長、及び／又は挿入もしくは欠失した塩基の種類を決定すること；
   （６"）該（５"）で決定された挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度を決定すること；
   （７"）該試験物質に曝露していない該細胞集団を対照群とし、該（１"）〜（６"）と同様の手順で、対照群における挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度を決定すること；
   （８"）該（６"）で得られた試験群における挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度から、該（７"）で得られた対照群における該変異頻度を引き算すること、
   を含む、方法。
8. 前記（２"）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
   前記（３"）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
   請求項７記載の方法。
9. 前記参照配列とマッチしない塩基を検出したリード配列上の部位を、該部位における該塩基のリード配列間での出現頻度に基づいて該部位を変異部位として決定する工程を行うことなく、塩基対置換型変異を有する変異部位として取得する、請求項１〜６のいずれか１項記載の方法。
10. 前記参照配列に対して塩基が挿入もしくは欠失したリード配列上の部位を、該挿入もしくは欠失部位のリード配列間での出現頻度に基づいて該部位を変異部位として決定する工程を行うことなく、挿入もしくは欠失変異を有する変異部位として取得する、請求項７又は８記載の方法。
11. 前記塩基対置換型変異が、１塩基対置換型変異、２塩基対置換型変異、又は３塩基対置換型変異である、請求項１〜６及び９のいずれか１項記載の方法。
12. 前記細胞集団が、サルモネラ菌細胞集団および大腸菌細胞集団からなる群より選択される少なくとも１種である、請求項１〜１１のいずれか１項記載の方法。
13. 前記サルモネラ菌がＳ．Ｔｙｐｈｉｍｕｒｉｕｍ ＬＴ−２株、ＴＡ１００株、ＴＡ９８株、ＴＡ１５３５株、ＴＡ１５３８株又はＴＡ１５３７株である、請求項１２記載の方法。
14. がん細胞における変異の評価方法であって、
    （１）がん細胞集団を試験群とし、そのＤＮＡを取得すること；
    （２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
    （３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
    （４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
    （５）取得した各変異を、塩基対の変異パターンに従って分類すること；
    （６）該（５）で得られた変異パターンの各々の変異頻度を決定すること；
    （７）該（１）〜（６）と同様の手順で、対照群における各塩基対の変異パターンの変異頻度を決定すること；
    （８）該（６）で得られた試験群における各変異パターンの変異頻度から、該（７）で得られた対照群における各変異パターンの変異頻度を引き算すること、
    を含む、方法。
15. 培養細胞の遺伝情報の評価方法であって、
    （１）培養細胞集団を試験群とし、そのＤＮＡを取得すること；
    （２）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
    （３）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
    （４）該（３）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
    （５）取得した各変異を、塩基対の変異パターンに従って分類すること；
    （６）該（５）で得られた変異パターンの各々の変異頻度を決定すること；
    （７）該（１）〜（６）と同様の手順で、対照群における各塩基対の変異パターンの変異頻度を決定すること；
    （８）該（６）で得られた試験群における各変異パターンの変異頻度から、該（７）で得られた対照群における各変異パターンの変異頻度を引き算すること、
    を含む、方法。
16. 前記（２）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
    前記（３）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、請求項１４又は１５記載の方法。
17. 前記（３）〜（５）が、
    前記リード配列に含まれる塩基を下記(i)〜(iv)に分けること、
    (i) 参照配列上の塩基がＡである位置に存在する塩基
    (ii) 参照配列上の塩基がＴである位置に存在する塩基
    (iii)参照配列上の塩基がＧである位置に存在する塩基
    (iv) 参照配列上の塩基がＣである位置に存在する塩基
    該リード配列に含まれる塩基の中から、参照配列と塩基がマッチしないものを検出し、当該塩基の存在する部位を、塩基対置換型変異を有する変異部位として取得すること、
    検出されたマッチしない塩基の各々について、該変異部位における変異前及び変異後の塩基対を取得すること、及び
    該変異部位の塩基対置換型変異を、変異前の塩基対と変異後の塩基対の種類に従って、ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴの６つの塩基対の変異パターンに分類すること、
    を含む、請求項１４〜１６のいずれか１項記載の方法。
18. がん細胞における変異の評価方法であって、
    （１'）がん細胞集団を試験群とし、そのＤＮＡを取得すること；
    （２'）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
    （３'）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
    （４'）該（３'）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
    （５'）取得した変異の各々について、該参照配列に基づいて、変異前の塩基と、該変異前の塩基の上流及び下流に隣接する塩基とを含むコンテクスト配列を決定すること；
    （６'）該（４'）で取得した各変異を、該（５'）で決定したコンテクスト配列及び変異後の塩基の種類に従ってタイプ分けすること；
    （７'）該（６'）で得られた変異タイプの各々の変異頻度を決定すること；
    （８'）該（１'）〜（７'）と同じ手順で、対照群における各変異タイプの変異頻度を決定すること；
    （９'）該（７'）で得られた試験群における各変異タイプの変異頻度から、該（８'）で得られた対照群における各変異タイプの変異頻度を引き算すること、
    を含む、方法。
19. 培養細胞の遺伝情報の評価方法であって、
    （１'）培養細胞集団を試験群とし、そのＤＮＡを取得すること；
    （２'）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
    （３'）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列と該参照配列とで塩基がマッチしない部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
    （４'）該（３'）で検出した部位を、塩基対置換型変異を有する変異部位として取得すること；
    （５'）取得した変異の各々について、該参照配列に基づいて、変異前の塩基と、該変異前の塩基の上流及び下流に隣接する塩基とを含むコンテクスト配列を決定すること；
    （６'）該（４'）で取得した各変異を、該（５'）で決定したコンテクスト配列及び変異後の塩基の種類に従ってタイプ分けすること；
    （７'）該（６'）で得られた変異タイプの各々の変異頻度を決定すること；
    （８'）該（１'）〜（７'）と同じ手順で、対照群における各変異タイプの変異頻度を決定すること；
    （９'）該（７'）で得られた試験群における各変異タイプの変異頻度から、該（８'）で得られた対照群における各変異タイプの変異頻度を引き算すること、
    を含む、方法。
20. 前記（２'）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
    前記（３'）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
    請求項１８又は１９記載の方法。
21. 前記（３'）〜（６'）が、
    前記リード配列に含まれる塩基を下記(i)〜(iv)に分けること、
    (i) 参照配列上の塩基がＡである位置に存在する塩基
    (ii) 参照配列上の塩基がＴである位置に存在する塩基
    (iii)参照配列上の塩基がＧである位置に存在する塩基
    (iv) 参照配列上の塩基がＣである位置に存在する塩基
    該リード配列に含まれる塩基の中から、参照配列と塩基がマッチしないものを検出し、当該塩基の存在する部位を、塩基対置換型変異を有する変異部位として取得すること、
    検出されたマッチしない塩基の各々について、該変異部位における変異前及び変異後の塩基対を取得すること、
    該変異部位の塩基対置換型変異を、変異前の塩基対と変異後の塩基対の種類に従って、ＡＴ→ＴＡ、ＡＴ→ＣＧ、ＡＴ→ＧＣ、ＧＣ→ＴＡ、ＧＣ→ＣＧ、及びＧＣ→ＡＴの６つの塩基対の変異パターンに分類すること、
    該変異部位における変異前の塩基と、該変異前の塩基の上流に隣接する１以上の塩基と、該変異前の塩基の下流に隣接する１以上の塩基とからなるコンテクスト配列を決定すること、及び
    該６つの塩基対の変異パターンと該コンテクスト配列に従って、該塩基対置換型変異をタイプ分けすること、
    を含む、請求項１８〜２０のいずれか１項記載の方法。
22. がん細胞における変異の評価方法であって、
    （１"）がん細胞集団を試験群とし、そのＤＮＡを取得すること；
    （２"）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
    （３"）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列における該参照配列に対して塩基が挿入もしくは欠失した部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
    （４"）該（３"）で検出した部位を、挿入もしくは欠失変異を有する変異部位として取得すること；
    （５"）取得した変異の各々について、挿入もしくは欠失の塩基長、及び／又は挿入もしくは欠失した塩基の種類を決定すること；
    （６"）該（５"）で決定された挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度を決定すること；
    （７"）該（１"）〜（６"）と同様の手順で、対照群における挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度を決定すること；
    （８"）該（６"）で得られた試験群における挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度から、該（７"）で得られた対照群における該変異頻度を引き算すること、
    を含む、方法。
23. 培養細胞の遺伝情報の評価方法であって、
    （１"）培養細胞集団を試験群とし、そのＤＮＡを取得すること；
    （２"）該ＤＮＡのフラグメントをシーケンシングし、各フラグメントにつき１つ以上のリード配列を得ること；
    （３"）該１つ以上のリード配列をそれぞれ参照配列と比較して、該リード配列における該参照配列に対して塩基が挿入もしくは欠失した部位を検出すること、ここで該参照配列は、該ＤＮＡ中の既知配列である；
    （４"）該（３"）で検出した部位を、挿入もしくは欠失変異を有する変異部位として取得すること；
    （５"）取得した変異の各々について、挿入もしくは欠失の塩基長、及び／又は挿入もしくは欠失した塩基の種類を決定すること；
    （６"）該（５"）で決定された挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度を決定すること；
    （７"）該（１"）〜（６"）と同様の手順で、対照群における挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度を決定すること；
    （８"）該（６"）で得られた試験群における挿入もしくは欠失部位の塩基長及び／又は挿入もしくは欠失した塩基の種類ごとの変異頻度から、該（７"）で得られた対照群における該変異頻度を引き算すること、
    を含む、方法。
24. 前記（２"）で得られたリード配列の塩基から、シーケンシングによる読み取りの信頼度の高い塩基を抽出することをさらに含み、かつ
    前記（３"）において、該抽出されたリード配列上の塩基を参照配列の塩基と比較する、
    請求項２２又は２３記載の方法。
25. 前記参照配列とマッチしない塩基を検出したリード配列上の部位を、該部位における該塩基のリード配列間での出現頻度に基づいて該部位を変異部位として決定する工程を行うことなく、塩基対置換型変異を有する変異部位として取得する、請求項１４〜２１のいずれか１項記載の方法。
26. 前記参照配列に対して塩基が挿入もしくは欠失したリード配列上の部位を、該挿入もしくは欠失部位のリード配列間での出現頻度に基づいて該部位を変異部位として決定する工程を行うことなく、挿入もしくは欠失変異を有する変異部位として取得する、請求項２２〜２４のいずれか１項記載の方法。
27. 前記塩基対置換変異が、１塩基対置換型変異、２塩基対置換型変異、又は３塩基対置換型変異である、請求項１４〜２１及び２５のいずれか１項記載の方法。
28. 前記検出に用いるリード配列の総量が、１×１０¹⁰ｂｐ以下である、請求項１〜２７のいずれか１項記載の方法。