(57)【要約】【課題】多くのリード数を得られる次世代シーケンサー用プライマーを提供する。【解決手段】5’-CAAGCAGAAGACGGCATACGAGAT-N5〜15-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3’（N5〜15は５〜１５塩基のインデックス配列）において、リード数を目的変数としインデックス配列における塩基の種類を説明変数とする推定式に基づいて算出したリード数の推定値が所定の値を超える塩基配列としてインデックス配列を設計する。【選択図】図１
2019110840.tif 000002
【特許請求の範囲】
【請求項１】
5’-CAAGCAGAAGACGGCATACGAGAT(配列番号１)-N5〜15-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG(配列番号２)-3’（N5〜15は５〜１５塩基のインデックス配列）の塩基配列を含み、リード数を目的変数としインデックス配列における塩基の種類を説明変数とする推定式に基づいて算出したリード数の推定値が所定の値を超える塩基配列として設計されたインデックス配列を含む次世代シーケンサー用プライマー。
【請求項２】
上記インデックス配列は８塩基（N=8、配列番号６７）であることを特徴とする請求項１記載の次世代シーケンサー用プライマー。
【請求項３】
上記推定式は、インデックス配列を構成するN個の塩基について、塩基の種類とそれに応じた係数とからなる項を含むことを特徴とする請求項１記載の次世代シーケンサー用プライマー。
【請求項４】
上記所定の値は、15000〜25000の値とすることを特徴とする請求項１記載の次世代シーケンサー用プライマー。
【請求項５】
配列番号２６２〜９６３からなる群から選ばれる１つの塩基配列からなることを特徴とする請求項１記載の次世代シーケンサー用プライマー。
【請求項６】
5’-CAAGCAGAAGACGGCATACGAGAT(配列番号１)-N5〜15-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG(配列番号２)-3’（ N5〜15は５〜１５塩基のインデックス配列）の塩基配列を含む次世代シーケンサー用プライマーについて、リード数を目的変数としインデックス配列における塩基の種類を説明変数とする推定式に基づいて、インデックス配列の塩基配列からリード数の推定値を算出し、算出したリード数の推定値が所定の値を超える塩基配列をインデックス配列の塩基配列として設計する工程と
上記工程で設計したインデックス配列を含むヌクレオチドを合成する工程とを含む次世代シーケンサー用プライマーの製造方法。
【請求項７】
上記インデックス配列は８塩基（N=8、配列番号６７）であることを特徴とする請求項６記載の次世代シーケンサー用プライマーの製造方法。
【請求項８】
上記推定式は、インデックス配列を構成するN個の塩基について、塩基の種類とそれに応じた係数とからなる項を含むことを特徴とする請求項６記載の次世代シーケンサー用プライマーの製造方法。
【請求項９】
上記所定の値は、15000〜25000の値とすることを特徴とする請求項６記載の次世代シーケンサー用プライマーの製造方法。
【請求項１０】
配列番号２６２〜９６３からなる群から選ばれる１つの塩基配列における25〜32番目をインデックス配列の塩基配列として設計することを特徴とする請求項６記載の次世代シーケンサー用プライマーの製造方法。
【請求項１１】
請求項１〜５いずれか一項記載の次世代シーケンサー用プライマーの塩基配列を、解析対象DNAの一方端部に有するDNA断片を含むDNAライブラリー。
【請求項１２】
上記解析対象DNAは、核酸増幅反応により得られた断片或いはゲノムDNAを断片化して得られた断片であることを特徴とする請求項１１記載のDNAライブラリー。
【請求項１３】
上記解析対象DNAは、ゲノムDNA及び高濃度のランダムプライマーを含む反応液にて核酸増幅反応することで得られる断片であって、当該ランダムプライマーは、上記次世代シーケンサー用プライマー3’末端側の塩基配列に対して相補的な塩基配列を5’末端側に有することを特徴とする請求項１１記載のDNAライブラリー。
【請求項１４】
請求項１〜５いずれか一項記載の次世代シーケンサー用プライマーと、上記次世代シーケンサー用プライマー3’末端側の塩基配列に対して相補的な塩基配列を5’末端側に有する解析対象DNAとを用いて核酸増幅反応を行う工程を含むDNAライブラリーの製造方法。
【請求項１５】
上記解析対象DNAは、核酸増幅反応により得られた断片或いはゲノムDNAを断片化して得られた断片であることを特徴とする請求項１４記載のDNAライブラリーの製造方法。
【請求項１６】
上記解析対象DNAは、ゲノムDNA及び高濃度のランダムプライマーを含む反応液にて核酸増幅反応することで得られる断片であって、当該ランダムプライマーは、上記次世代シーケンサー用プライマー3’末端側の塩基配列に対して相補的な塩基配列を5’末端側に有することを特徴とする請求項１４記載のDNAライブラリーの製造方法。
【請求項１７】
請求項１１〜１３いずれか一項記載のDNAライブラリーを次世代シーケンサーにより解析し、当該DNAライブラリーに含まれるDNA断片の塩基配列を決定することを特徴とするDNA解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数検体の同時解析を可能とするインデックスを有する次世代シーケンサー用プライマー並びにその製造方法、当該次世代シーケンサー用プライマーを用いたDNAライブラリー並びにその製造方法、及びDNAライブラリーを用いたゲノムDNA解析方法に関する。
【背景技術】
【０００２】
次世代シーケンサー（Next Generation Sequencer：NGS）は、多数のDNA断片に関する塩基配列を並列に読み出すことができる装置である。一例として、イルミナ社の次世代シーケンサーでは、ランダムに切断された数千万〜数億のDNA断片の両端部にそれぞれアダプターを連結し、アダプターを介して5’末端側をフローセル上に固定させる。次に、フローセル上に予め固定された5’末端側のアダプターと、DNA断片の3’末端側のアダプター配列とをアニールさせて、ブリッジ状のDNA断片を形成させる。この状態でDNAポリメラーゼによる核酸増幅反応を行うことで、多数の1本鎖DNA断片を局所的に増幅して固定することができる。そして、次世代シーケンサーでは、得られた1本鎖DNAを鋳型として、シーケンシングを行うことで、1回の解析において40〜200Gbという膨大な配列情報を得ることができる。
【０００３】
次世代シーケンサーにおけるシーケンシングは、蛍光標識したdNTPの取り込みを蛍光顕微鏡により解析するという方法を採用している。具体的には、3’末端側が保護基でブロックされるとともに蛍光標識されたdNTPを用いる。DNAポリメラーゼにより、1本鎖DNA断片に対して相補的なdNTPが取り込まれ、これをレーザー光により励起し、蛍光顕微鏡により蛍光を読み取る。そして、dNTPから保護基を除去し、同様にして次の塩基を解析する。このように、次世代シーケンサーでは、フローセルに固定された1本差DNAについて1塩基ずつ連続して解析する。
【０００４】
特に、次世代シーケンサーでは、解析対象のDNA断片に連結するアダプターにインデックス（インデックス配列やバーコード配列とも称する）を付与することで、複数サンプル由来のDNA断片を識別することができる。すなわち、上述のように、一回の解析において膨大な配列情報を得ることができるが、配列情報に含まれるインデックス配列を指標して如何なるサンプルに由来するDNA断片に関する配列情報であるか判断することができる。
【０００５】
ところが、非特許文献１に記載されるように、インデックス配列を利用した次世代シーケンサーによる解析では、インデックス配列によっては、リード数が大きくばらつくことがあるといった問題があった。しかし、非特許文献１では、インデックス配列の性能に関する相違について体系的に解析されておらず、インデックス配列を利用した次世代シーケンサーによる解析精度が不十分であるといった問題があった。
【０００６】
これまで、一対のプライマーにそれぞれ異なる塩基配列からなるユニバーサル・テール配列を付与し、これら一対のプライマーを用いてマルチプレックスPCRを行い、次世代シーケンサーに供するアンプリコン群（同じインデックス配列を有する）を得る方法が知られている（特許文献１）。また、大量の試料を分析する際の効率化を図る目的で、アダプターとインデックスと標的DNA特異的配列を全て含む一対のプライマーを用いて、次世代シーケンサー用のDNAライブラリーを作製する方法が知られている（特許文献２）。特許文献２に開示されるプライマーは、ヒトミトコンドリアＤＮＡの超可変領域等の標的DNAに特異的に結合するプライマーとNGSライブラリーの作製に必須なアダプタープライマー、インデックスプライマー及びシーケンシングプライマーを全て含む統合プライマーである。
【先行技術文献】
【非特許文献】
【０００７】
【非特許文献１】David W. Craig et al., Nat Methods. 2008 October; 5(10): 887-893
【特許文献】
【０００８】
【特許文献１】米国特許公開公報  US 2016/0326572 A1
【特許文献２】特開2017-97935号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかしながら、インデックスを含む次世代シーケンサー用プライマーを使用したときに、インデックスの配列によってリード数が著しく減少するといった問題点は未解明のままであり、当該問題点を解決する技術は知られていないのが現状である。そこで、本発明は、このような実情に鑑みて、インデックスの塩基配列とリード数との関係を解明することで、多くのリード数を得られる次世代シーケンサー用プライマー及びその製造方法を提供することを目的とする。また、本発明は、多くのリード数を得られる次世代シーケンサー用プライマーを用いたDNAライブラリー並びにその製造方法、及びDNAライブラリーを用いたゲノムDNA解析方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明者らは、上述した目的を達成するため鋭意検討した結果、インデックスを含む次世代シーケンサー用プライマーにおけるインデックス配列とリード数との関係を明らかにし、インデックスの塩基配列からリード数を推定できることを明らかにした。そして、推定されたリード数が所定の値を超えるインデックス配列を設計することで、多くのリード数を得られる次世代シーケンサー用プライマーを製造できることを見いだし、本発明を完成するに至った。
【００１１】
本発明は以下を包含する。
（１）5’-CAAGCAGAAGACGGCATACGAGAT(配列番号１)-N5〜15-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG(配列番号２)-3’（N5〜15は５〜１５塩基のインデックス配列）の塩基配列を含み、リード数を目的変数としインデックス配列における塩基の種類を説明変数とする推定式に基づいて算出したリード数の推定値が所定の値を超える塩基配列として設計されたインデックス配列を含む次世代シーケンサー用プライマー。
（２）上記インデックス配列は８塩基（N=8、配列番号６７）であることを特徴とする（１）記載の次世代シーケンサー用プライマー。
（３）上記推定式は、インデックス配列を構成するN個の塩基について、塩基の種類とそれに応じた係数とからなる項を含むことを特徴とする（１）記載の次世代シーケンサー用プライマー。
（４）上記所定の値は、15000〜25000の値とすることを特徴とする（１）記載の次世代シーケンサー用プライマー。
（５）配列番号２６２〜９６３からなる群から選ばれる１つの塩基配列からなることを特徴とする（１）記載の次世代シーケンサー用プライマー。
（６）5’-CAAGCAGAAGACGGCATACGAGAT(配列番号１)-N5〜15-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG(配列番号２)-3’（N5〜15は５〜１５塩基のインデックス配列）の塩基配列を含む次世代シーケンサー用プライマーについて、リード数を目的変数としインデックス配列における塩基の種類を説明変数とする推定式に基づいて、インデックス配列の塩基配列からリード数の推定値を算出し、算出したリード数の推定値が所定の値を超える塩基配列をインデックス配列の塩基配列として設計する工程と
上記工程で設計したインデックス配列を含むヌクレオチドを合成する工程とを含む次世代シーケンサー用プライマーの製造方法。
（７）上記インデックス配列は８塩基（N=8、配列番号６７）であることを特徴とする（６）記載の次世代シーケンサー用プライマーの製造方法。
（８）上記推定式は、インデックス配列を構成するN個の塩基について、塩基の種類とそれに応じた係数とからなる項を含むことを特徴とする（６）記載の次世代シーケンサー用プライマーの製造方法。
（９）上記所定の値は、15000〜25000の値とすることを特徴とする（６）記載の次世代シーケンサー用プライマーの製造方法。
（１０）配列番号２６２〜９６３からなる群から選ばれる１つの塩基配列における25〜32番目をインデックス配列の塩基配列として設計することを特徴とする（６）記載の次世代シーケンサー用プライマーの製造方法。
（１１）上記（１）〜（５）いずれか記載の次世代シーケンサー用プライマーの塩基配列を、解析対象DNAの一方端部に有するDNA断片を含むDNAライブラリー。
（１２）上記解析対象DNAは、核酸増幅反応により得られた断片或いはゲノムDNAを断片化して得られた断片であることを特徴とする（１１）記載のDNAライブラリー。
（１３）上記解析対象DNAは、ゲノムDNA及び高濃度のランダムプライマーを含む反応液にて核酸増幅反応することで得られる断片であって、当該ランダムプライマーは、上記次世代シーケンサー用プライマー3’末端側の塩基配列に対して相補的な塩基配列を5’末端側に有することを特徴とする（１１）記載のDNAライブラリー。
（１４）上記（１）〜（５）いずれか記載の次世代シーケンサー用プライマーと、上記次世代シーケンサー用プライマー3’末端側の塩基配列に対して相補的な塩基配列を5’末端側に有する解析対象DNAとを用いて核酸増幅反応を行う工程を含むDNAライブラリーの製造方法。
（１５）上記解析対象DNAは、核酸増幅反応により得られた断片或いはゲノムDNAを断片化して得られた断片であることを特徴とする（１４）記載のDNAライブラリーの製造方法。
（１６）上記解析対象DNAは、ゲノムDNA及び高濃度のランダムプライマーを含む反応液にて核酸増幅反応することで得られる断片であって、当該ランダムプライマーは、上記次世代シーケンサー用プライマー3’末端側の塩基配列に対して相補的な塩基配列を5’末端側に有することを特徴とする（１４）記載のDNAライブラリーの製造方法。
（１７）上記（１１）〜（１３）いずれか記載のDNAライブラリーを次世代シーケンサーにより解析し、当該DNAライブラリーに含まれるDNA断片の塩基配列を決定することを特徴とするDNA解析方法。
【発明の効果】
【００１２】
本発明によれば、インデックス配列に起因してリード数が減少するといった不都合を回避でき、多くのリード数を得ることができ次世代シーケンサー用プライマー及びその製造方法、並びに当該次世代シーケンサー用プライマーを用いて作製したDNAライブラリー及びその製造方法を提供することができる。
【００１３】
また、本発明に係る次世代シーケンサー用プライマーを利用して作製したDNAライブラリーを利用することで、サンプル毎のデータ量（リード数）がばらつくことを防止し、高精度なDNA解析を行うことができる。
【図面の簡単な説明】
【００１４】
【図１】次世代シーケンサー用プライマーを用いた次世代シーケンサーによる解析スキームを模式的に示す特性図である。
【図２】次世代シーケンサー用プライマーP7におけるインデックス配列の種類と、リード数との関係を示す特性図である。
【図３】次世代シーケンサー用プライマーP5におけるインデックス配列の種類と、リード数との関係を示す特性図である。
【図４】次世代シーケンサー用プライマーP7におけるインデックス配列を構成する塩基の種類に基づいて推定式で算出した推定リード数と、実際に測定したリード数との関係を示す特性図である。
【図５】次世代シーケンサー用プライマーP5におけるインデックス配列を構成する塩基の種類に基づいて推定式で算出した推定リード数と、実際に測定したリード数との関係を示す特性図である。
【発明を実施するための形態】
【００１５】
以下、本発明を詳細に説明する。
【００１６】
本発明に係る次世代シーケンサー用プライマーは、5’-CAAGCAGAAGACGGCATACGAGAT(配列番号１)-N5〜15-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG(配列番号２)-3’（N5〜15は５〜１５塩基のインデックス配列）の塩基配列を含んでいる。上記塩基配列の表記において、N5〜15で示されるインデックス配列とは、詳細を後述するように設計された配列であり、例えば、サンプルを識別するための索引として利用できる配列である。すなわち、インデックス配列を後述のように設計すれば、インデックス配列は如何なる配列であっても良い。例えば、複数のサンプルに対してそれぞれ異なるインデックス配列を設計することによって、次世代シーケンサーで解析した結果の塩基配列データについて、インデックス配列に基づいてサンプルへの帰属を明確に理解することができる（マルチプレックス解析）。
【００１７】
ここで、次世代シーケンサー用プライマーを用いた次世代シーケンサーによる解析スキームを模式的に図１に示す。本発明に係る次世代シーケンサー用プライマーは、所謂、イルミナ社が提供する次世代シーケンサー（NGS）においてP7と呼称されるプライマーに基づいている。次世代シーケンサーに供するDNAライブラリーを調製する際、次世代シーケンサー用プライマーP7と次世代シーケンサー用プライマーP5とを用いてPCRを行う。なお、次世代シーケンサー用プライマーP5は、図１においてインデックス配列を有する構成を示したが、インデックス配列を有しない構成であっても良い。
【００１８】
図１に示すように、次世代シーケンサー用プライマーP7と次世代シーケンサー用プライマーP5とを用いたPCRにより、次世代シーケンサー用プライマーP7と次世代シーケンサー用プライマーP5との間に解析対象DNAが配置されたDNA断片を得ることができる。当該PCRにより得られたDNA断片の群をDNAライブラリー又は次世代シーケンサー用DNAライブラリーと称する。
【００１９】
図１に示すように、上述のように得られたDNAライブラリーは、イルミナ社の次世代シーケンサーにて解析される。次世代シーケンサーは、DNAライブラリーを解析した結果として多数のリードに関する塩基配列データ（すなわち、解析対象DNAの塩基配列を含む塩基配列データ）を出力する。
【００２０】
ここで解析対象DNAとは、特に限定されず、ゲノム解析の対象である生物由来のゲノムDNA、エピゲノム解析のために調製したDNA又はトランスクリプト解析のたに調製されたDNA等に由来するDNA断片、若しくはゲノムDNA等を鋳型として得られる増幅産物（アンプリコン）とすることができる。
【００２１】
〔インデックス配列の設計〕
本発明に係る次世代シーケンサー用プライマーにおけるインデックス配列は、次世代シーケンサーを用いてDNAライブラリーの塩基配列を解析する際に十分に多くのリード数を得ることができるように設計される。具体的に、インデックス配列を設計する際には、リード数を目的変数とし、インデックス配列における塩基の種類を説明変数とする推定式を準備する。当該推定式を準備するには、先ず、特定の塩基配列からなるインデックス配列を有する次世代シーケンサー用プライマーを実際に合成し、この次世代シーケンサー用プライマーを使用したときのリード数を解析する。そして、インデックス配列の塩基配列と得られたリード数とから上記推定式を算出することができる。
【００２２】
具体的に、推定式を算出する際には、特に限定されないが、一般的に相関解析に利用される各種のアルゴリズムを適用することができる。具体的には、Lasso（Least Absolute Shrinkage and Selection Operator）等のスパースモデリング方法を適用することができる。スパースモデリング方法のうちL1型正則化推定法によれば、重要でないパラメータの係数をOに推定することができ、適切なパラメータのみからなる推定式を得ることができる。
【００２３】
例えば、推定式は、インデックス配列における所定の位置について塩基毎に算出された係数を含む項と切片とを含み、インデックス配列の全ての位置について特定の塩基を指定するとリード数の推定値を算出できる。このとき、LASSO等のスパースモデリング方法を適用することで、重要でないパラメータ、すなわちリード数の増減に対する関与の低い塩基については係数を0とした推定式を算出することができる。
【００２４】
更に具体的に、例えば、８塩基長のインデックス配列（N=8）、すなわち5’末端側から3’末端側に向かって1番目から8番目の塩基で構成されるインデックス配列を設計する場合、推定式は、1番目の塩基がアデニン或いはグアニンの場合に推定リード数が低減し、シトシンの場合に推定リード数が増加するような式とすることができる。同様に、推定式は、2番目の塩基がアデニン又はチミンの場合に推定リード数が増加し、グアニンの場合に推定リード数が低減するような式とすることができる。同様に、推定式は、3番目の塩基がアデニン又はチミンの場合に推定リード数が増加し、グアニンの場合に推定リード数が低減するような式とすることができる。同様に、推定式は、4番目の塩基がアデニンの場合に推定リード数が増加し、シトシン又はグアニンの場合に推定リード数が低減するような式とすることができる。同様に、推定式は、5番目の塩基がアデニンの場合に推定リード数が増加し、グアニンの場合に推定リード数が低減するような式とすることができる。同様に、推定式は、6番目の塩基がシトシンの場合に推定リード数が減少し、チミンの場合に推定リード数が増加するような式とすることができる。同様に、推定式は、7番目の塩基がアデニンの場合に推定リード数が減少し、グアニンの場合に推定リード数が増加するような式とすることができる。同様に、推定式は、8番目の塩基がグアニンの場合に推定リード数が減少し、チミンの場合に推定リード数が増加するような式とすることができる。
【００２５】
なかでも、推定式は、1番目の塩基がアデニン又はグアニンである場合、2番目の塩基がグアニンである場合及び8番目の塩基がグアニンである場合に、特に推定リード数が大きく減少するような式とすることができる。また、推定式は、2番目の塩基がアデニン又はチミンである場合、3番目の塩基がアデニンである場合、5番目の塩基がアデニンである場合に、特に推定リード数が大きく増加するような式とすることができる。
【００２６】
より具体的には、インデックス配列における塩基の種類毎の係数及び切片を下記表のように設定することができる。
【００２７】
【表１】
2019110840.tif 000003
【００２８】
表１において、A1は、インデックス配列の5’末端側から3’末端側の方向において1番目の塩基がアデニンである場合を意味している。他の表記も同様に、インデックス配列の5’末端側から3’末端側の方向における位置と塩基の種類とを意味している。
【００２９】
一例として上記表に示した推定式を用いることで、インデックス配列を所定の塩基配列としたときに得られるリード数の推定値を算出することができる。上述したように得られた推定式は、表１に具体的に示した推定式に限定されず、実測値との相関が非常に高いといった特徴を有する（相関係数が0.9以上）。したがって、当該推定式を使用してインデックス配列の塩基配列毎に推定リード数を計算し、推定リード数が所定の値以上になるような塩基配列を選択することで、実際に多くのリード数を得ることができる。
【００３０】
ここで、インデックス配列を設計するにあたり、閾値とする推定リード数は、特に限定されず、目的とするデータ解析に応じて適宜設定することができる。例えば、推定リード数の閾値として15000〜25000の値を設定することができ、17000〜23000の値を設定することが好ましく、19000〜21000の値を設定することがより好ましく、19500〜20500の値を設定することが更に好ましい。詳細を後述する実施例に示すように、インデックス配列の塩基配列に関して、リード数が著しく低い一群の塩基配列が同定されている。この一群の塩基配列について推定リード数を計算したところ、最大でリード数が20000程度、具体的には20051.8であることが判明している。したがって、上述した閾値を20000程度、例えば20052とすることで、実際に多くのリード数を得ることができるインデックス配列を設計することができる。
【００３１】
上記表１に示した推定式を用いて、推定リード数が20052を超えるインデックス配列を有する次世代シーケンサー用プライマーの塩基配列を特定している（詳細は後述の実施例参照）。すなわち、本発明に係る次世代シーケンサー用プライマーとしては、例えば、配列番号２６２〜９６３からなる群から選ばれる１つの塩基配列からなるものを挙げることができる。
【００３２】
〔次世代シーケンス技術への応用〕
以上で説明した本発明に係る次世代シーケンサー用プライマーを用いた核酸増幅反応により、次世代シーケンサーで使用するDNAライブラリー（次世代シーケンサー用DNAライブラリー）を作製することができる。具体的には、図１に示すように、本発明に係る次世代シーケンサー用プライマーP7、解析対象DNA及び次世代シーケンサー用プライマーP5を含む反応液でPCRを行う。解析対象DNAの両末端は、それぞれ次世代シーケンサー用プライマーP7と共通する配列及び次世代シーケンサー用プライマーP5と共通する配列を有する。このため、核酸増幅反応においては、次世代シーケンサー用プライマーP7及びP5が解析対象DNAの両末端にアニールして解析対象DNAを鋳型とした核酸増幅反応が進行する。その結果、図１に示すように、次世代シーケンサー用プライマーP5、解析対象DNA及び次世代シーケンサー用プライマーP7をこの順で有する核酸断片の群（DNAライブラリー）を得ることができる。
【００３３】
得られたDNAライブラリーは、その両末端に次世代シーケンサー用プライマーP5及びP7を有するため、イルミナ社の次世代シーケンサーに供することができる。そして、得られたDNAライブラリーは、上述した本発明に係る次世代シーケンサー用プライマーP7を有するため、次世代シーケンサーにより多くのリード数を解析することができる。具体的には、上述した推定式で算出される推定リード数に近い値のリード数の解析が可能となる。
【００３４】
なお、使用される次世代シーケンサーは、ブリッジPCR法とSequencing-by-synthesis法により、フローセル上で目的DNAを増幅させ、合成しながらシーケンシングを行うといった原理である。
【００３５】
本発明に係る次世代シーケンサー用プライマーP7、解析対象DNA及び次世代シーケンサー用プライマーP5を含む反応液を用いた核酸増幅反応は、特に限定されず、通常の核酸増幅反応の条件を適用することができる。例えば、この反応液は、鋳型となる解析対象DNA、上述した次世代シーケンサー用プライマーP5及びP7、DNAポリメラーゼ、基質となるデオキシヌクレオチド三リン酸（dNTP、すなわちdATP、dCTP、dTTP及びdGTPの混合物）及びバッファーを含む。
【００３６】
特に次世代シーケンサー用プライマーP5及びP7の濃度としては、０．０１〜５．０μＭとすることができ、０．１〜２．５μＭとすることが好ましく、０．３〜０．７μＭとすることが最も好ましい。
【００３７】
また、核酸増幅反応において鋳型となる解析対象DNAは、特に限定されないが、反応液の量を５０μｌとしたときに、０．１〜１０００ｎｇとすることが好ましく、１〜５００ｎｇとすることがより好ましく、５〜２００ｎｇとすることが更に好ましく、１０〜１００ｎｇとすることが最も好ましい。
【００３８】
鋳型となる解析DNA断片の調整方法は、特に限定されず、詳細を後述するランダムプライマーを用いた核酸増幅反応が終わった後の反応液をそのまま使用しても良いし、当該反応液から解析DNAを精製したものを使用しても良い。
【００３９】
また、核酸増幅反応に使用するDNAポリメラーゼの種類、基質となるデオキシヌクレオチド三リン酸（dNTP、すなわちdATP、dCTP、dTTP及びdGTPの混合物）の濃度、バッファー組成、温度サイクル条件については、通常の核酸増幅反応に採用される条件とすることができる。また、次世代シーケンサー用プライマーを用いた核酸増幅反応において、ホットスタート法を採用しても良いし、核酸増幅反応によって増幅断片を取得するものでも良い。
【００４０】
以上のように、ランダムプライマーを用いて取得した第1のDNA断片を鋳型とし、次世代シーケンサー用プライマーを用いて増幅した第2のDNA断片を使用することで、次世代シーケンス装置に適用可能なDNAライブラリーを簡便に作製することができる。
【００４１】
〔解析対象DNA〕
解析対象DNAの一例として、ゲノムDNA等を鋳型として得られる増幅産物（アンプリコン）を挙げることができる。ここで、増幅産物を得るには、例えば、任意の塩基配列を有するプライマー（以下、ランダムプライマー）を高濃度となるように調整した反応液で核酸増幅反応を行う方法を挙げることができる。ここで、高濃度とは、通常の核酸増幅反応におけるプライマー濃度と比較して高濃度であることを意味する。すなわち、本方法は、通常の核酸増幅反応におけるプライマー濃度と比較して高濃度のランダムプライマーを使用することに特徴を有している。ここで、反応液に含まれる鋳型としては、解析対象の生物から調整したゲノムDNAを使用することができる。
【００４２】
なお、本方法において、対象の生物種には何ら限定されず、ヒトを含む動物、植物、微生物、ウイルス等いかなる生物種も対象とすることができる。すなわち、本方法によれば、如何なる生物種からも解析対象DNAとして多数の増幅産物を得ることができる。
【００４３】
本方法では、ランダムプライマーの濃度を上記のように規定することによって、高い再現性で核酸断片（核酸断片群）を増幅することができる。ここで、再現性とは、同一の鋳型及び同一のランダムプライマーを用いて複数回の核酸増幅反応を行った場合に、複数回の核酸増幅反応の間で増幅される核酸断片が一致する程度を意味する。つまり、高い再現性（再現性が高い）とは、同一の鋳型及び同一のランダムプライマーを用いて複数回の核酸増幅反応を行った場合に、複数回の核酸増幅反応の間で増幅される核酸断片の一致度が高いことを意味する。
【００４４】
再現性の高低については、例えば、同一の鋳型及び同一のランダムプライマーを用いて複数回の核酸増幅反応を行い、各回で得られた増幅断片を電気泳動し得られた蛍光ユニット(Fluorescence Unit: FU)についてスピアマンの順位相関係数を算出し、当該係数に基づいて評価することができる。スピアマンの順位相関係数とは、一般的にρで表され、一例としてはρ＞０．９をもって再現性有りと評価することができる。
【００４５】
本方法に使用できるランダムプライマーとしては、その配列には何ら限定されず、例えば９〜３０塩基長のヌクレオチドを使用することができる。特に、ランダムプライマーとは、任意の配列を有する、９〜３０塩基長のヌクレオチドであって、ヌクレオチドの種類（配列の種類）は特に限定されず、１種類以上のヌクレオチド、好ましくは１〜１００００種類のヌクレオチド、より好ましくは１〜１０００種類のヌクレオチド、より好ましくは１〜１００種類のヌクレオチド、最も好ましくは１〜９６種類のヌクレオチドを意味する。ランダムプライマーとして上述の範囲のヌクレオチド（ヌクレオチド群）を使用することによって、より高い再現性で増幅核酸断片を得ることができる。なお、ランダムプライマーとして、複数のヌクレオチドを含む場合、全てのヌクレオチドが同じ塩基長（９〜３０塩基長）である必要はなく、異なる塩基長の複数のヌクレオチドを含んでいても良い。
【００４６】
通常、核酸増幅反応を用いて特定のアンプリコンを得るためには、当該アンプリコンに応じてプライマーの塩基配列を設計する。例えば、ゲノムDNA等の鋳型DNAにおけるアンプリコンに対応する位置を挟み込むように一対のプライマーを設計する。この場合、プライマーは、鋳型に含まれる特定の領域にハイブリダイズするように設計されるため「特異的プライマー」と呼称することができる。
【００４７】
これに対して、ランダムプライマーは、特定のアンプリコンを得る目的で設計されるプライマーとは異なり、鋳型DNAにおける特定の領域にハイブリダイズするように設計されるのではなく、ランダムなアンプリコンを得るために設計される。ランダムプライマーは、その塩基配列が如何なる配列であってもよく、鋳型DNAに含まれる相補的な領域に偶発的にハイブリダイズすることでランダムなアンプリコン増幅に関与できる。
【００４８】
すなわち、ランダムプライマーとは、上述のように、ランダムなアンプリコン増幅に関与する任意配列を有するヌクレオチドということができる。ここで任意配列とは、何ら限定されないが、例えば、アデニン、グアニン、シトシン及びチミンの群から無作為に選択された塩基配列として設計しても良いし、特定の塩基配列として設計しても良い。特定の塩基配列としては、例えば、制限酵素認識配列を含む塩基配列や、次世代シーケンサーに使用するアダプター配列を有する塩基配列を挙げることができる。
【００４９】
ランダムプライマーとして複数種類のヌクレオチドを設計する場合、アデニン、グアニン、シトシン及びチミンの群から無作為に選択して所定の長さの塩基配列を複数設計する方法が適用できる。また、ランダムプライマーとして複数種類のヌクレオチドを設計する場合、特定の塩基配列からなる共通部分と、任意の塩基配列からなる非共通部分とからなる塩基配列を複数設計する方法も適用できる。ここで、非共通部分は、アデニン、グアニン、シトシン及びチミンの群から無作為に選択された塩基配列としても良いし、アデニン、グアニン、シトシン及びチミンからなる４種類の塩基の全ての組み合わせ、又はこれら全ての組み合わせから選ばれる一部の組み合わせとすることができる。共通部分は、特に限定されず如何なる塩基配列でもよいが、例えば、制限酵素認識配列を含む塩基配列や、次世代シーケンサーに使用するアダプター配列を有する塩基配列、特定の遺伝子ファミリーに共通する塩基配列とすることができる。
【００５０】
複数のランダムプライマーとして、４種類の塩基から無作為に選択して所定の長さの塩基配列を複数設計する場合、全体の３０％以上、好ましくは５０％以上、より好ましくは７０％以上、更に好ましくは９０％以上が、７０％以下の同一性、好ましくは６０％以下の同一性、より好ましくは５０％以下の同一性、最も好ましくは４０％以下の同一性となるように設計することが好ましい。複数のランダムプライマーとして、４種類の塩基から無作為に選択して所定の長さの塩基配列を複数設計する場合であって、上述の範囲のヌクレオチドについて上記範囲の同一性となるように設計することで、対象生物種のゲノムDNA全体に亘って増幅断片を得ることができる。すなわち、増幅断片の均一性を高めることができる。
【００５１】
複数のランダムプライマーとして、特定の塩基配列からなる共通部分と、任意の塩基配列からなる非共通部分とからなる塩基配列を複数設計する場合、例えば、３’末端側の数塩基を非共通部分とし、残りの５’末端側を共通部分とするように設計することができる。３’末端側のｎ個の塩基を非共通部分とすれば、４ｎ種類のランダムプライマーを設計することができる。ここで、ｎ個としては、１〜５個とすることができ、好ましくは２〜４個、より好ましくは２〜３個である。
【００５２】
例えば、共通部分と非共通部分とからなるランダムプライマーとしては、５’末端側を次世代シーケンサーに使用するアダプター配列（共通部分）とし、３’末端側を２塩基（非共通部分）とした、合計１６種類のランダムプライマーを設計することができる。なお、３’末端側を３塩基（非共通部分）とすれば、合計６４種類のランダムプライマーを設計することができる。ランダムプライマーの種類が多くなるほど、対象生物種のゲノムDNA全体に亘って増幅断片をより網羅的に得ることができる。したがって、共通部分と非共通部分とからなるランダムプライマーを設計する場合、３’末端側の塩基は３塩基とすることが好ましい。
【００５３】
ただし、例えば、共通部分と３塩基の非共通部分とからなる６４種類の塩基配列を設計した後、これら６４種類の塩基配列から選ばれる６３種類以下のランダムプライマーを使用してもよい。言い換えると、６４種類のランダムプライマーの全てを使用した場合と比較して、６３種類以下のランダムプライマーを使用した場合のほうが核酸増幅反応や、次世代シーケンサーを用いた解析において優れた結果を示す場合がある。具体的には、６４種類のランダムプライマーを使用した場合、特定の核酸増幅断片のリード数が著しく多くなるような場合がある。この場合には、６４種類のランダムプライマーのなかから当該特定の核酸増幅断片の増幅に関与する１又は複数のランダムプライマーを除いた残りの６３種類以下のランダムプライマーを使用したほうが良好な解析結果が得られることになる。
【００５４】
なお、共通部分と２塩基の非共通部分とからなる１６種類のランダムプライマーを設計した場合も同様に、１６種類のランダムプライマーから選ばれる１５種類以下のランダムプライマーを使用した場合に、核酸増幅反応や次世代シーケンサーを用いた解析において優れた結果を示す場合がある。
【００５５】
一方、ランダムプライマーとして使用するヌクレオチドは、特に、GC含量が５〜９５％の範囲となるように設計することが好ましく、１０〜９０％の範囲となるように設計することがより好ましく、１５〜８０％の範囲となるように設計することが更に好ましく、２０〜７０％の範囲となるように設計することが最も好ましい。ランダムプライマーとしてGC含量が上述の範囲のヌクレオチドの集合を使用することによって、より高い再現性で増幅核酸断片を得ることができる。なお、GC含量とは、ヌクレオチド鎖全体に含まれるグアニン及びシトシンの割合である。
【００５６】
さらに、ランダムプライマーとして使用するヌクレオチドは、特に、全体の長さに対して連続塩基が８０％以下となるように設計することが好ましく、７０％以下となるように設計することがより好ましく、６０％以下となるように設計することが更に好ましく、５０％以下となるように設計することが最も好ましい。或いは、ランダムプライマーとして使用するヌクレオチドは、特に、連続塩基の数が８個以下となるように設計することが好ましく、７個以下となるように設計することがより好ましく、６個以下となるように設計することが更に好ましく、５個以下となるように設計することが最も好ましい。ランダムプライマーとして連続塩基数が上述の範囲のヌクレオチドの集合を使用することによって、より高い再現性で増幅核酸断片を得ることができる。
【００５７】
さらにまた、ランダムプライマーとして使用するヌクレオチドは、特に、分子内に６塩基長以上、好ましくは５塩基以上、より好ましくは４塩基以上の相補領域を有しないように設計することが好ましい。ヌクレオチドに上記範囲の相補領域を有しないように設計することで分子内における二本鎖形成を防止でき、より高い再現性で増幅核酸断片を得ることができる。
【００５８】
さらにまた、ランダムプライマーとして複数種類のヌクレオチドを設計する場合、特に、複数のヌクレオチド間において６塩基長以上、好ましくは５塩基以上、より好ましくは４塩基以上の相補領域を有しないように設計することが好ましい。複数のヌクレオチド間に上記範囲の相補領域を有しないように設計することでヌクレオチド間の二本鎖形成を防止でき、より高い再現性で増幅核酸断片を得ることができる。
【００５９】
さらにまた、ランダムプライマーとして複数種類のヌクレオチドを設計する場合、特に、３’末端側の６塩基以上、好ましくは５塩基以上、より好ましくは４塩基以上が相補的な配列にならないように設計することが好ましい。複数のヌクレオチドの３’末端側の上記範囲において相補的な配列を有しないように設計することでヌクレオチド間の二本鎖形成を防止でき、より高い再現性で増幅核酸断片を得ることができる。
【００６０】
なお、相補領域及び相補的な配列とは、例えば８０〜１００％の同一性を有する領域及び配列（例えば５塩基長の領域であれば、４塩基又は５塩基が相補的である領域及び配列）、或いは９０〜１００％の同一性を有する領域及び配列（例えば５塩基長の領域であれば５塩基が相補的である領域及び配列）を意味する。
【００６１】
さらにまた、ランダムプライマーとして使用するヌクレオチドは、核酸増幅反応における温度サイクル条件（特に、アニーリング温度）に適したTm値となるように設計することが好ましい。特に限定されないが、Tm値は、最近接塩基対法、Wallace法及びGC％法等の公知の計算方法により算出することができる。具体的に、ランダムプライマーとして使用するヌクレオチドは、特にTm値が１０〜８５℃、好ましくは１２〜７５℃、より好ましくは１４〜７０℃、最も好ましくは１６〜６５℃となるように設計することが好ましい。ヌクレオチドのTm値を上記範囲となるように設計することで、核酸増幅反応における所定の温度サイクル条件（特に、所定のアニーリング温度）下において、より高い再現性で増幅核酸断片を得ることができる。
【００６２】
さらにまた、ランダムプライマーとして複数種類のヌクレオチドを設計する場合、特に、複数のヌクレオチド間において各ヌクレオチドのTm値のばらつきが５０℃以下、好ましくは４５℃以下、より好ましくは４０℃以下、最も好ましくは３５℃以下となるように設計することが好ましい。複数のヌクレオチド間でTm値のばらつきが上記範囲となるように設計することで、核酸増幅反応における所定の温度サイクル条件（特に、所定のアニーリング温度）下において、より高い再現性で増幅核酸断片を得ることができる。
【００６３】
本方法では、上述したランダムプライマー及び鋳型としてのゲノムDNAを用いた核酸増幅反応によって、多数の増幅断片を取得する。特に、核酸増幅反応において、反応液中のランダムプライマーの濃度を、通常の核酸増幅反応におけるプライマー濃度と比較して高濃度とする。これにより、高い再現性を達成しながら多数の増幅断片をゲノムDNAを鋳型として得ることができる。
【００６４】
ここで、核酸増幅反応とは、鋳型となるゲノムDNA、上述したランダムプライマー、DNAポリメラーゼ、基質となるデオキシヌクレオチド三リン酸（dNTP、すなわちdATP、dCTP、dTTP及びdGTPの混合物）及びバッファーからなる反応液で、所定の温度サイクル条件を加えることで増幅断片を合成する反応である。なお、核酸増幅反応には、反応液中に所定濃度のMg2+が必要であり、上述した組成においてバッファーにMgCl2が含まれている。バッファーにMgCl2が含まれていない場合、上記組成に加えてMgCl2が含まれることとなる。
【００６５】
特に核酸増幅反応において、ランダムプライマーの濃度は、ランダムプライマーの塩基長に応じて適宜設定することが好ましい。ここで、ランダムプライマーの塩基長は、異なる塩基長の複数種類のヌクレオチドをランダムプライマーとして使用する場合には、その平均値（単純平均でもよいし、ヌクレオチド量を加味した加重平均でもよい）とすることができる。
【００６６】
具体的には、９〜３０塩基長のランダムプライマーを用い、当該ランダムプライマー濃度を４〜２００μＭとする条件、好ましくは４〜１００μＭとする条件で核酸増幅反応を行う。この条件であれば、核酸増幅反応により、高い再現性を達成しながら多数の増幅断片、特に１００〜５００塩基長の多数の増幅断片を得ることができる。
【００６７】
また、核酸増幅反応において鋳型となるゲノムDNAは、特に限定されないが、反応液の量を５０μｌとしたときに、０．１〜１０００ｎｇとすることが好ましく、１〜５００ｎｇとすることがより好ましく、５〜２００ｎｇとすることが更に好ましく、１０〜１００ｎｇとすることが最も好ましい。鋳型となるゲノムDNAの量をこの範囲とすることで、ランダムプライマーからの増幅反応が阻害されることなく、高い再現性を達成しながら多数の増幅断片を得ることができる。
【００６８】
ゲノムDNAの調整方法は、特に限定されず、従来公知の方法を適用することができる。また、市販されているキットを利用することで、対象の生物種から簡便にゲノムDNAを調整することができる。なお、ゲノムDNAとしては、従来公知の方法や市販のキットにより生物から抽出されたものをそのまま使用しても良いし、生物から抽出したものを精製したものでも良いし、制限酵素処理や超音波処理した後のものを使用しても良い。
【００６９】
また、核酸増幅反応においてDNAポリメラーゼとしては、特に限定されず、核酸増幅反応のための温度サイクル条件下でDNAポリメラーゼ活性を有する酵素を使用することができる。具体的には、通常の核酸増幅反応に使用される耐熱性DNAポリメラーゼを使用することができる。例えば、DNAポリメラーゼとしては、Taq DNAポリメラーゼ等の好熱細菌由来DNAポリメラーゼ、KOD DNAポリメラーゼやPfu DNAポリメラーゼ等の超好熱Archaea由来DNAポリメラーゼを挙げることができる。特に、核酸増幅反応においては、上述したランダムプライマーとともに、DNAポリメラーゼとしてはPfu DNAポリメラーゼを使用することが好ましい。これらDNAポリメラーゼを使用することで、高い再現性を達成しながら多数の増幅断片をより確実に得ることができる。
【００７０】
さらに、核酸増幅反応において、基質となるデオキシヌクレオチド三リン酸（dNTP、すなわちdATP、dCTP、dTTP及びdGTPの混合物）の濃度は、特に限定されないが、５μＭ〜０．６ｍＭとすることができ、１０μＭ〜０．４ｍＭとすることが好ましく、２０μＭ〜０．２ｍＭとすることがより好ましい。基質となるdNTPの濃度をこの範囲とすることで、DNAポリメラーゼによる誤った取り込みによるエラーの発生を防止し、高い再現性を達成しながら多数の増幅断片を得ることができる。
【００７１】
さらに、核酸増幅反応において、バッファーとしては、特に限定されないが、上述のようにMgCl2を含み、例えばTris-HCl（pH8.3）及びKClを含む溶液を挙げることができる。ここで、Mg2+の濃度としては、特に限定されないが、例えば０．１〜４．０ｍＭとすることができ、０．２〜３．０ｍＭとすることが好ましく、０．３〜２．０ｍＭとすることがより好ましく、０．５〜１．５ｍＭとすることが更に好ましい。反応液中のMg2+濃度をこの範囲に設定することで、高い再現性を達成しながら多数の増幅断片を得ることができる。
【００７２】
さらにまた、核酸増幅反応における温度サイクル条件としては、特に限定されず、通常の温度サイクルを採用することができる。具体的に温度サイクルとは、先ず、鋳型のゲノムDNAを一本鎖に乖離するための最初の熱変性温度、その後、「熱変性温度→アニーリング温度→伸長反応温度」を複数回（例えば、２０〜４０回）行い、最後に必要であれば、所定時間伸長反応温度とし、最後に保存のための温度とするサイクルを例示することができる。
【００７３】
熱変性温度としては、例えば９３〜９９℃、好ましくは９５〜９８℃、より好ましくは９７〜９８℃とすることができる。アニーリング温度としては、上述したランダムプライマーのTm値にもよるが、例えば３０〜７０℃、好ましくは３５〜６８℃、より好ましくは３７〜６５℃とすることができる。伸長反応温度としては、例えば７０〜７６℃、好ましくは７１〜７５℃、より好ましくは７２〜７４℃とすることができる。また、保存のための温度としては例えば４℃とすることができる。
【００７４】
また、最初の熱変性は、上述した温度範囲において、例えば５秒〜１０分、好ましくは１０秒〜５分、より好ましくは３０秒〜２分とすることができる。「熱変性温度→アニーリング温度→伸長反応温度」のサイクルにおける熱変性は、上述した温度範囲において、例えば２秒〜５分、好ましくは５秒〜２分、より好ましくは１０秒〜１分とすることができる。「熱変性温度→アニーリング温度→伸長反応温度」のサイクルにおけるアニーリングは、上述した温度範囲において、例えば１秒〜３分、好ましくは３秒〜２分、より好ましくは５秒〜１分とすることができる。「熱変性温度→アニーリング温度→伸長反応温度」のサイクルにおける伸長反応は、上述した温度範囲において、例えば１秒〜３分、好ましくは３秒〜２分、より好ましくは５秒〜１分とすることができる。
【００７５】
また、本方法としては、ホットスタート法を採用した核酸増幅反応によって増幅断片を取得するものでも良い。ホットスタート法とは、「熱変性温度→アニーリング温度→伸長反応温度」のサイクル前のミスプライミングやプライマーダイマーに由来する非特異増幅を防ぐ方法である。ホットスタート法では、抗DNAポリメラーゼ抗体を結合させたり、化学修飾を行うことで、DNAポリメラーゼ活性を抑制させた状態の酵素を使用する。この状態では、DNAポリメラーゼ活性が抑制され、温度サイクル前の非特異的な反応を防止することができる。ホットスタート法では、最初の温度サイクル時に温度を高く設定することでDNAポリメラーゼ活性が回復し、その後の核酸増幅反応が進行することとなる。
【００７６】
以上のように、上述した９〜３０塩基長のランダムプライマーを使用し、反応液中の当該ランダムプライマー濃度を４〜２００μＭとして核酸増幅反応を行うことで、ゲノムDNAを鋳型としてランダムプライマーが多数の増幅断片を得ることができる。上述した９〜３０塩基長のランダムプライマーを使用し、反応液中の当該ランダムプライマー濃度を４〜２００μＭとした場合、非常に再現性の高い核酸増幅反応となる。すなわち、上述した核酸増幅反応によれば、非常に高い再現性を達成しながら多数の増幅断片を得ることができる。
【００７７】
また、上述した９〜３０塩基長のランダムプライマーを使用し、反応液中の当該ランダムプライマー濃度を４〜２００μＭとして核酸増幅反応を行うことで、特に、ゲノムDNAを鋳型として約１００〜５００塩基長の多数の増幅断片を得ることができる。この約１００〜５００塩基長の多数の増幅断片は、例えば次世代シーケンサーによる塩基配列の大量解析に適したサイズであり、高精度な配列情報を得ることができる。すなわち、本発明によれば、約１００〜５００塩基長といったDNA断片を作製することができる。
【００７８】
さらに、上述した９〜３０塩基長のランダムプライマーを使用し、反応液中の当該ランダムプライマー濃度を４〜２００μＭとして核酸増幅反応を行うことで、特に、ゲノムDNAの全体に亘って均一に増幅断片を得ることができる。言い換えると、当該ランダムプライマーを用いた核酸増幅反応では、ゲノムDNAの所定の領域に偏ってDNA断片が増幅されるのではなく、ゲノム全体に分散してDNA断片が増幅される。
【００７９】
なお、上述したランダムプライマーを使用して核酸増幅反応を行った後、得られた増幅断片に対して制限酵素処理、サイズセレクション処理及びシーケンスキャプチャー処理等を行うことができる。増幅断片に対してこれら制限酵素処理、サイズセレクション処理及びシーケンスキャプチャー処理を行うことで、得られた増幅断片のなかから特定の増幅断片（特定の制限酵素部位を有する断片、特定のサイズ範囲の増幅断片、特定の配列を有する増幅断片）を得ることができる。
【００８０】
〔ゲノムDNA解析方法〕
上述のように作製されたアンプリコンを使用することで、遺伝子型解析等のゲノムDNA解析を行うことができる。上述のように作製したアンプリコンは、非常に高い再現性を有しており、次世代シーケンサーに適したサイズを有しており、ゲノム全体に亘って均一性を有している。したがって、上述のように作製したアンプリコンは、DNAマーカー（遺伝マーカー、遺伝子マーカーとも称される）として使用することができる。ここで、DNAマーカーとは、ゲノムDNA内に存在する特徴的な塩基配列を広く意味する。また、DNAマーカーとしては、特に、遺伝的形質に関連する目印となるゲノム上の塩基配列とすることもできる。DNAマーカー、例えば遺伝子型同定、連鎖地図、遺伝子マッピング、マーカーを利用した選抜工程を含む育種、マーカーを利用した戻し交配、量的形質遺伝子座のマッピング、バルクセグリガント分析、品種識別、又は連鎖不均衡マッピング等に利用することができる。
【００８１】
特に、本発明に係る次世代シーケンサー用プライマーを使用することで、次世代シーケンサー等を使用して、上述のように作製されたアンプリコンの塩基配列を決定し、得られた塩基配列に基づいてDNAマーカーの存否を確認することができる。
【００８２】
一例としては、得られた塩基配列のリード数からDNAマーカーの存否を確認することができる。ここで、次世代シーケンサーとは、特に限定されないが、上述した本発明に係る次世代シーケンサー用プライマーを利用することができ、第2世代シーケンサーとも呼称され、数千万のDNA断片の塩基配列を同時並行的に決定できる塩基配列決定装置を意味する。次世代シーケンサーにおけるシーケンシング原理としては、特に限定されず、例えばブリッジPCR法とSequencing-by-synthesis法により、フローセル上で目的DNAを増幅させ、合成しながらシーケンシングを行うといった原理を挙げることができる。より具体的に、次世代シーケンサーとしては、イルミナ社（Illumina）のMiniSeq、MiSeq、NextSeq、HiSeq及びHiSeq Xシリーズを挙げることができる。
【００８３】
また、他の例としては、上述のように作製されたアンプリコンについて得られた塩基配列を参照用の塩基配列と比較することでDNAマーカーの存否を確認することができる。ここで、参照用の塩基配列とは、基準となる既知の配列を意味し、例えば、データベースに格納された既知配列とすることができる。すなわち、所定の生物について、上述のようにアンプリコンを作製し、その塩基配列を決定し、参照用の塩基配列と比較する。そして、参照用の塩基配列と相違する塩基配列を、当該所定の生物に関するDNAマーカー（ゲノムDNA内に存在する特徴的な塩基配列）とすることができる。また、特定したDNAマーカーについては、定法に従って更に解析することによって、遺伝的形質（表現型）に関連性を決定することができる。すなわち、上述のように特定したDNAマーカーの中から、表現型に関連するDNAマーカー（選抜マーカーと称する場合もある）を特定することができる。
【００８４】
さらに、他の例としては、上述のように作製されたアンプリコンについて得られた塩基配列を、他の生物由来のゲノムDNA又は他の組織由来のゲノムDNAを用いて作製したアンプリコンの塩基配列と比較することでDNAマーカーの存否を確認することができる。すなわち、２以上の生物又は２つの異なる組織について、それぞれ上述のようにアンプリコンを作製し、それらの塩基配列を決定し、塩基配列同士を比較する。そして、相違する塩基配列を、供試した生物又は組織に関するDNAマーカー（ゲノムDNA内に存在する特徴的な塩基配列）とすることができる。また、特定したDNAマーカーについては、定法に従って更に解析することによって、遺伝的形質（表現型）に関連性を決定することができる。すなわち、上述のように特定したDNAマーカーの中から、表現型に関連するDNAマーカー（選抜マーカーと称する場合もある）を特定することができる。
【００８５】
或いは、上述のように本発明に係る次世代シーケンサー用プライマーを用いて解析された塩基配列データは、微生物等の多種多様性を調べるメタゲノム解析や腫瘍組織などの体細胞ゲノム変異解析、マイクロアレイを利用した遺伝子型解析、倍数性の判定解析、染色体数の算出解析、染色体の増減解析、染色体の部分的挿入・欠失・複製・転座解析、外来ゲノムの混入解析、親子判別解析、交配種子純度検定解析といった解析に利用することができる。
【実施例】
【００８６】
以下、実施例を用いて本発明を更に詳細に説明するが、本発明の技術的範囲は以下の実施例に限定されるものではない。
【００８７】
〔実施例１〕
本実施例では、イルミナ社の次世代シーケンサーに使用される次世代シーケンサー用プライマーP7及びP5について、インデックス配列と次世代シーケンサーから得られるデータ量との関係を明らかにし、次世代シーケンサーを用いた解析に適した次世代シーケンサー用プライマーを開発した。
【００８８】
１．材料
本実施例では、イネ品種日本晴からDNeasy Plant Mini kit(QIAGEN)を用いてゲノムDNAを抽出及び精製したものを使用した。
【００８９】
２．方法
２．１.ランダムプライマーの設計
ランダムプライマーは、イルミナ社の次世代シーケンサー用のアゲプターNextera adapterにおける3'末端の10塩基（GTTACACACG）と、当該10塩基の3'末端にTGCを除く3塩基の任意の配列を付加した全長13塩基からなる63種類の塩基配列を設計した（表２）。
【００９０】
【表２】
2019110840.tif 000004
【００９１】
３．２．解析対象DNAの作製
〔１．材料〕で説明したイネのゲノムDNA（30ng）に最終濃度0.2mM dNTP mixture、10mM MgC12及び1.25unitのDNA Polymerase（TAKARA、PrimeSTAR）に、40μMのランゲムプライマーをそれぞれ加え、最終反応量50μlでPCRを行った。PCRの温度条件は、98℃を2分の後、98℃を10秒間、50℃を15秒間、72℃を20秒間を1サイクルとして30サイクル行い、その後、4℃で保存する条件とした。このPCRにより、ランダムプライマーによってイネゲノムDNAを鋳型として多数の解析対象DNAを作製した。
【００９２】
３．３．次世代シーケンサー用DNAライブラリーの作製
上記３．２で作製した解析対象DNA（１μl）に最終濃度0.2mM dNTP mixture、1.0mM MgC12及び1.25unitｎDNAPolymerase(TAKARA、PrimeSTAR)を加え、更に、0.25μMの次世代シーケンサー用プライマーP7及びP5をそれぞれ加えた反応液（50μl）を準備した。
【００９３】
次世代シーケンサー用プライマーP7及びP5は、表３に示したイルミナ社のNextera adaptorの配列情報をもとに作製した。なお、下記表においてアスタリスクで示した８塩基がインデックス配列である。ただし、P5のインデックス配列に対して、P7のインデックス配列は逆相補鎖である。
【００９４】
【表３】
2019110840.tif 000005
【００９５】
詳細には、次世代シーケンサー用プライマーP7において、96種のインデックス配列を有する次世代シーケンサー用プライマーP7を設計した（表４）。表４に示した次世代シーケンサー用プライマーP7を使用する場合、次世代シーケンサー用プライマーP5としてはAATGATACGGCGACCACCGAGATCTACACCGCGCAGATCGTCGGCAGCGTCAGATGTGTATAAGAGACAGを用いた。表４に示した次世代シーケンサー用プライマーP7を用いて作製したDNAライブラリーをDNAライブラリー1と称する。
【００９６】
【表４】
2019110840.tif 000006
【００９７】
2019110840.tif 000007
【００９８】
2019110840.tif 000008
【００９９】
一方、次世代シーケンサー用プライマーP5において、98種のインデックス配列を有する次世代シーケンサー用プライマーP5を設計した（表５）。表５に示した次世代シーケンサー用プライマーP5を使用する場合、次世代シーケンサー用プライマーP7としてはCAAGCAGAAGACGGCATACGAGATTCGTCAGAGTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGを用いた。表５に示した次世代シーケンサー用プライマーP5を用いて作製したDNAライブラリーをDNAライブラリー2と称する。
【０１００】
【表５】
2019110840.tif 000009
【０１０１】
2019110840.tif 000010
【０１０２】
2019110840.tif 000011
【０１０３】
DNAライブラリー1及び2を作製するための核酸増幅反応における温度条件は、95℃を2分間の後、98℃を15秒間、55℃を15秒間及び72℃を20秒間を1サイクルとして25サイクル行い、その後、72℃を1分後とし、4℃で保存する条件とした。そして、核酸増幅反応により作製したDNAライブラリーはMinElute PCR Purification Kit（QIAGEN）で精製した。
【０１０４】
３．４．次世代シーケンサーによる解析
上記３．３で作製したDNAライブラリー1及び2については、MiSeq Reagent Kit V2 500 Cycte（イルミナ）を用いてリード長100塩基のペアエンド条件のもと、それぞれMiSeqにより解析した。DNAライブラリー1及び2について解析の結果として得られたリード数をそれぞれ表４及び５に示した。
【０１０５】
３．５．インデックスとデータ量の解析
上記３．４で解析した結果として得られた次世代シーケンサー用プライマー毎のリード数をもとに、次世代シーケンサー用プライマーに含まれるインデックス配列の８塩基とデータ量の関係をGLMNET LASSO法により解析し、インデックス配列に含まれる塩基の種類とによるリード数の推定式を作成した。作成した推定式より算出した推定リード数と、実際の測定されたリード数との相関係数を算出した。
【０１０６】
４．結果
４．１．次世代シーケンサー用プライマー毎のリード数の分布
次世代シーケンサー用プライマーのインデックスの違いによるデータ量への影響を評価するため、インデックス配列が異なる次世代シーケンサー用プライマーP7及びP5について、プライマー毎のリード数の分布を調査した。次世代シーケンサー用プライマーP7のインデックス配列が異なるDNAライブラリー1では、全体のリード数の平均値66,9617に対し、リード数が15,000以下のインデックスが9種(94%)存在し、その平均値は2,492.2と全体平均値の3.7%だった(図２)。また、リード数が40,000以下のインデックスが22種(220%)存在し、その平均値は16,2376と全体平均値の24.3%だった。一方、次世代シーケンサー用プライマーP5のインデックス配列が異なるDNAライブラリー2では、全体のリード数の平均値191,5231で、すべて100,000以下のリード数を示した（図３）。
【０１０７】
４．２．インデックス配列とリード数の関係
上述したDNAライブラリー1及び2それぞれを次世代シーケンサーで解析して得られたデータをもとに、インデックス配列とリード数との関係をGLMNET LASSO法で解析し、インデックス配列における塩基の種類とリード数との推定式を作成した。すなわち、リード数を目的変数とし、インデックス配列における塩基の種類を説明変数とする推定式を作成した。DNAライブラリー1について作成した推定式により算出した推定リード数と実測したリード数との関係を図４に示した。図４に示すように、インデックス配列を構成する塩基の種類に基づいて推定式で算出した推定リード数と、実際に測定したリード数との相関係数は
r=0.94069だった。また、図４に示しように、推定リード数と実測したリード数のプロットは、リード数が極めて少ないグループ1、リード数が多いグループ3、その中間のグループ2の3グループに分類できた。グループ1に属する次世代シーケンサー用プライマーP7のインデックス配列に基づいて算出した推定リード数の最大値は20,051.8だった。また、グループ3に属する次世代シーケンサー用プライマーP7のインデックス配列に基づいて算出した推定リード数は、いずれも50,000以上であった。
【０１０８】
一方、同様に、DNAライブラリー2についても推定式を算出し、インデックス配列を構成する塩基の種類に基づいて算出した推定リード数と、実測したリード数との相関係数を算出した結果、相関係数はr=0.57295であった（図５）。すなわち、次世代シーケンサー用プライマーP5については、インデックス配列を構成する塩基の種類と、リード数との間に相関が見られなかった。
【０１０９】
４．３次世代シーケンサー用プライマーP7におけるインデックス配列の設計
次世代シーケンサー用プライマーP7のインデックスが異なるDNAライブラリー1では、リード数とGLMNET LASSO法で得られた推定値との間での相関係数が0.9以上を示した。この結果から、ライブラリー1についてGLMNET LASSO法により作成した推定式は、インデックス配列を構成する塩基の種類からリード数の推定値を高精度に算出できるものと評価できた。具体的に、GLMNET LASSO法により作成した推定式は、以下のように、インデックス配列における所定の位置について塩基毎に算出された係数を含む項と切片とを含んだ式である。
推定値=-65033.1×(A1)+1326.4×(C1)-16997×(G1)+10936.3×(A2)-12399.2×(G2)+11712.9×(T2)+12112.2×(A3)-623.5×(G3)+5964.4×(T3)+6884.5×(A4)-5664.4×(C4)-6049.9×(G4)+9257×(A5)-6210.8×(G5)-644×(C6)+3.2×(T6)-3575.9×(A7)+1013.l×(G7)-8607.7×(G8)+6490.3×(T8)+81720.7
上記推定式において（A1）は、インデックス配列の5’末端側から3’末端側の方向において1番目の塩基がアデニンである場合に"1"を代入し、それ以外は"0"を代入するパラメータである。他の表記も同様に、インデックス配列の5’末端側から3’末端側の方向における位置と塩基の種類とを意味し、"1"又は"0"が代入されるパラメータである。
【０１１０】
以上のように作成した推定式を用い、推定値が20,052以上（グループ2又は3と推定される）を示す次世代シーケンサー用プライマーP7と、推定値が50,000以上（グループ3と推定される）を示す次世代シーケンサー用プライマーP7とを選抜した(それぞれ表６及び７)。なお、表６及び７中の塩基配列において、Nは、アデニン、シトシン、グアニン及びチミンのいずれか任意の塩基を意味する。
【０１１１】
【表６】
2019110840.tif 000012
【０１１２】
2019110840.tif 000013
【０１１３】
2019110840.tif 000014
【０１１４】
2019110840.tif 000015
【０１１５】
2019110840.tif 000016
【０１１６】
2019110840.tif 000017
【０１１７】
2019110840.tif 000018
【０１１８】
2019110840.tif 000019
【０１１９】
2019110840.tif 000020
【０１２０】
2019110840.tif 000021
【０１２１】
2019110840.tif 000022
【０１２２】
2019110840.tif 000023
【０１２３】
2019110840.tif 000024
【０１２４】
2019110840.tif 000025
【０１２５】
2019110840.tif 000026
【０１２６】
2019110840.tif 000027
【０１２７】
2019110840.tif 000028
【０１２８】
2019110840.tif 000029
【０１２９】
2019110840.tif 000030
【０１３０】
【表７】
2019110840.tif 000031
【０１３１】
2019110840.tif 000032
【０１３２】
2019110840.tif 000033
【０１３３】
2019110840.tif 000034
【０１３４】
2019110840.tif 000035
【０１３５】
2019110840.tif 000036
【０１３６】
2019110840.tif 000037
【０１３７】
2019110840.tif 000038
【０１３８】
2019110840.tif 000039
【０１３９】
2019110840.tif 000040
【０１４０】
2019110840.tif 000041
【０１４１】
2019110840.tif 000042
【０１４２】
2019110840.tif 000043
【０１４３】
2019110840.tif 000044
【０１４４】
2019110840.tif 000045
【０１４５】
５．考察
本実施例では、イルミナ社の次世代シーケンサー用プライマーにおけるインデックス配列とリード数について解析した結果、次世代シーケンサー用プライマーP5では明確な関係は見られなかったものの、次世代シーケンサー用プライマーP7では相関係数0.9以上と明確な相関関係が見られた。特に、GLMNET LASSO法により作製した推定式を用いてインデックス配列を構成する塩基の種類から推定リード数を算出し、推定リード数と実際のリード数とから、インデックス配列の異なる次世代シーケンサー用プライマーP7は、3つのグループに分類できることが明らかとなった。すなわち、リード数が15,000以下と極めて少なく、次世代シーケンサーを用いた解析への利用が困難と考えられるグループ１を特定することができた。このグループ１に属する次世代シーケンサー用プライマーP7の推定リード数の最大値が20,051.8であった。したがって、本実施例では、インデックス配列を構成する塩基の種類をパラメータとする推定式にて算出した推定リード数が20,052以上となる全てのインデックス配列を選抜した（表６）。選抜したインデックス配列を有する次世代シーケンサー用プライマーP7（表６）は、次世代シーケンサーに利用したときに多くのリード数を得ることができると考えられる。
【０１４６】
また、グループ３に属する次世代シーケンサー用プライマーP7の推定リード数は50,000以上であったことから、インデックス配列を構成する塩基の種類をパラメータとする推定式にて算出した推定リード数が50,000以上となる全てのインデックス配列を選抜した（表７）。選抜したインデックス配列を有する次世代シーケンサー用プライマーP7（表７）は、次世代シーケンサーに利用したときに更に多くのリード数を得ることができると考えられる。
【０１４７】
本実施例で設計された次世代シーケンサー用プライマーP7を利用することで、次世代シーケンサーからより安定したデータが得られることが期待できる。
【図１】
2019110840.tif 000046
【図２】
2019110840.tif 000047
【図３】
2019110840.tif 000048
【図４】
2019110840.tif 000049
【図５】
2019110840.tif 000050
【配列表】
2019110840000001.app