(57)【要約】  ｃＤＮＡ内でヌクレオチド配列、ｃＤＮＡ試料中のヌクレオチド配列の頻度を決定するための方法、並びに複雑性の低下及び（高処理量）配列決定を使用する、遺伝子の配列情報を必要とせずに、前記遺伝子の相対的な転写産物レベルを（偏向せずに）決定するための方法が記載されている。
2009520500.tif 000002
【特許請求の範囲】
【請求項１】
（ａ）ｃＤＮＡを準備する段階；
（ｂ）前記ｃＤＮＡの少なくとも部分について複雑性の低下を実施してｃＤＮＡ断片を含むｃＤＮＡの第一のライブラリを得る段階；
（ｃ）高処理量配列決定によって、第一のライブラリの前記ｃＤＮＡ断片のヌクレオチド配列の少なくとも一部を決定する段階；
（ｄ）段階ｄ）の第一のライブラリの前記ｃＤＮＡ断片のヌクレオチド配列を整列させて第一のライブラリの連結断片を生成する段階；及び
（ｅ）前記ｃＤＮＡのヌクレオチド配列を決定する段階
を含む、ｃＤＮＡのヌクレオチド配列を決定するための方法。
【請求項２】
（ａ）ｃＤＮＡを準備する段階；
（ｂ）前記ｃＤＮＡの少なくとも部分について複雑性の低下を実施してｃＤＮＡ断片を含むｃＤＮＡの第一のライブラリを得る段階；
（ｃ）配列決定によって、第一のライブラリの前記ｃＤＮＡ断片のヌクレオチド配列の少なくとも一部を決定する段階；及び
（ｄ）ヌクレオチド配列の頻度を決定する段階
を含む、ヌクレオチド配列の頻度を決定するための方法。
【請求項３】
（ａ）第一のｃＤＮＡ試料に対して請求項２に記載の方法を実施することによって、第一のｃＤＮＡ試料中のヌクレオチド配列の頻度を決定する段階；
（ｂ）第二の及び／又はさらなるｃＤＮＡ試料に対して請求項２に記載の方法を実施することによって、第二の及び／又はさらなるｃＤＮＡ試料において同一のヌクレオチド配列の頻度を決定する段階；並びに
（ｃ）第一のｃＤＮＡ試料中のヌクレオチド配列の頻度を、第二の及び／又はさらなるｃＤＮＡ試料中の同一のヌクレオチド配列の頻度と比較して前記ヌクレオチド配列の相対的な転写レベルを得る段階
を含む、ｃＤＮＡ試料中のヌクレオチド配列の相対的な転写レベルを決定するための方法。
【請求項４】
（ａ）第一のｃＤＮＡ試料を準備する段階；
（ｂ）第一のｃＤＮＡ試料に対して複雑性の低下を実施して第一のライブラリを得る段階；
（ｃ）第一のライブラリにタグ付けして第一のタグ付きライブラリを得る段階；
（ｄ）好ましくは各ｃＤＮＡ試料に対して異なるタグを使用して、第二の及び／又はさらなるｃＤＮＡ試料を用いて、連続して又は同時に段階（ａ）及び（ｂ）を実施して第二の及び／又はさらなるタグ付きライブラリを得る段階；
（ｅ）第一のタグ付きライブラリ並びに第二の及び／又はさらなるタグ付きライブラリを組み合わせて組み合わされたライブラリを得る段階；
（ｆ）前記組み合わされたライブラリのヌクレオチド配列の少なくとも一部を配列決定によって決定する段階；
（ｇ）第一のｃＤＮＡ試料並びに第二の及び／又はさらなるＤＮＡ試料中のヌクレオチド配列の頻度を決定する段階；並びに
（ｈ）第一のｃＤＮＡ試料中のヌクレオチド配列の頻度を、第二の及び／又はさらなるｃＤＮＡ試料中のヌクレオチド配列の頻度と比較して前記ｃＤＮＡ試料中のヌクレオチド配列の相対的な転写レベルを得る段階
を含む、ｃＤＮＡ試料中のヌクレオチド配列の相対的な転写レベルを決定するための方法。
【請求項５】
複雑性の低下が、増幅化断片長多形性技術、指標付き連結、ゲノム分割、遺伝子発現及びその修飾の連続分析、大量並行サイン配列決定、リアルタイム多重連結依存的プローブ増幅、高適用範囲発現特徴づけ、普遍的マイクロアレイシステム、トランスクリプトーム減算法、断片ディスプレイ、ディファレンシャルディスプレイ及び順序化ディファレンシャルディスプレイからなる群から選択される方法によって実施される、請求項１から４の何れか一項に記載の方法。
【請求項６】
配列決定が、高処理量配列決定によって実施される、請求項２から５の何れか一項に記載の方法。
【請求項７】
高処理量配列決定が、ビーズ等の固体支持体上で実施される、請求項１から６の何れか一項に記載の方法。
【請求項８】
高処理量配列決定が、合成による配列決定、好ましくはピロシーケンスに基づいている、請求項６又は７の何れか一項に記載の方法。
【請求項９】
高処理量配列決定が、
（ｃ１）配列決定アダプターを断片に連結する段階；
（ｃ２）各々、単一の断片と徐冷するビーズに、配列決定アダプター連結断片を徐冷する段階；
（ｃ３）各々、単一ビーズを含む油中水微小反応器中で前記ビーズを乳化する段階；
（ｃ４）乳濁液ＰＣＲを実施してビーズの表面上で配列決定アダプター連結断片を増幅する段階；
（ｃ５）増幅された配列決定アダプター連結断片を含有するビーズを選択し／濃縮する段階；
（ｃ６）各々単一ビーズを含むウェル中にビーズを負荷する段階；及び
（ｃ７）ピロリン酸塩シグナルを発生させる段階
を含む、請求項６から７の何れか一項に記載の方法。
【請求項１０】
（ａ）少なくとも１つの制限エンドヌクレアーゼでｃＤＮＡを消化して制限断片へと断片化する段階；
（ｂ）制限断片の一方又は両方の末端と適合性のある１つの末端を有する少なくとも１つの二本鎖合成オリゴヌクレオチドアダプターと前記制限断片を連結してアダプターが連結された制限断片を生成する段階；
（ｃ）前記アダプターが連結された制限断片を、少なくとも１つのアダプターの一部と、及び制限エンドヌクレアーゼの認識配列の残りの部分の一部と相補的なヌクレオチド配列部分を含むプライマー配列を有する１つ又はそれ以上のオリゴヌクレオチドプライマーと、ハイブリッド形成条件下で接触させる段階；及び
（ｄ）ハイブリッド形成された１つ又はそれ以上のオリゴヌクレオチドプライマーの伸長によって前記アダプターが連結された制限断片を増幅する段階
を含む方法によって複雑性の低下が実施される、請求項１から９の何れかにに記載の方法。
【請求項１１】
プライマーが、プライマー配列の３’末端に、選択された配列をさらに含み、１から１０個の選択的ヌクレオチドを含む前記選択された配列が、制限エンドヌクレアーゼの認識配列の残りの部分に直接隣接して配置される部分と相補的である、請求項１０に記載の方法。
【請求項１２】
プライマー配列の３’末端の選択された配列が、１から８個、好ましくは１から５個、より好ましくは１から３個の選択的ヌクレオチドを含む、請求項１０又は１１の何れか一項に記載の方法。
【請求項１３】
前記アダプターが、同定子配列をさらに含む、請求項１０から１２の何れか一項に記載の方法。
【請求項１４】
タグが、同定子配列である、請求項４に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、分子生物学及び遺伝学の分野に関する。本発明は、高処理量配列決定技術の使用に基づいて転写産物の配列を決定するための改良された戦略に関する。本発明はさらに、偏向されていない転写産物の特徴づけのための改良された戦略に関する。
【背景技術】
【０００２】
転写産物の特徴づけは、現代のバイオテクノロジー研究において使用される基礎技術の１つである。転写産物の特徴づけの主要な適用領域は、複雑な形質にかかわる遺伝子の発見である。これには、医薬の開発のための標的の同定（標的発見）、生体分子の合成を調節する生化学的経路（発酵産業）の解明、植物及び動物の育種のための複雑な形質の解析（遺伝子発見）及びその他の多くの目的のために、（ヒトの）疾病に関与する遺伝子を発見すること等、生物学的現象の広範な範囲が含まれる。
【０００３】
第二の適用領域は、逆の経路に従う。すなわち、複雑な表現型を予測するために、遺伝子（の選択されたサブセット）の転写産物特性の定型的な診断決定のために転写産物の特徴づけを使用する。本カテゴリーにおける例は、ヒト乳癌（Ｖａｎ  ｄｅ  Ｖｉｊｖｅｒ  ｅｔ  ａｌ．，２００２，Ｎ．Ｅｎｇｌ．Ｊ．Ｍｅｄ．，ｖｏｌ．３４７，２５：１９９９−２００９；ｖａｎ’ｔ  Ｖｅｅｒ  ｅｔ  ａｌ．，２００２，Ｂｒｅａｓｔ  Ｃａｎｃｅｒ  Ｒｅｓ．，ｖｏｌ．５（１）：５７−８；ｗｗｗ．ａｇｅｎｄｉａ．ｃｏｍ）及び乳頭状腎細胞癌（Ｙａｎｇ  ｅｔ  ａｌ．，２００５）の臨床的な予後の分子上の分類、診断及び予測である。分離された集団において回収された転写産物特徴づけデータに基づいて、関連遺伝子を同定するためのアプローチが、Ｓｃｈａｄｔ及び共同研究者（２００５，Ｓｃｉ．ＳＴＫＥ，ｖｏｌ．２９６：ｐe４０）によって記載されている。簡潔には、特徴づけは、ライフサイエンス研究において極めて重要である。
【０００４】
転写産物の特徴づけのための技術は、過去１０年にわたって急速に進化してきた。１９９０年代の初期（ＰＣＲが広範に利用可能になった直後）まで、転写産物の特徴づけは、ノザンブロット分析又はＲＮＡｓｅ保護アッセイによって実施されていた。これらの技術は、かなり特異的であり、且つ感受性が高い（特にＲＮＡｓｅ保護アッセイ）が、これらの技術の限界は、唯一の又は数個の遺伝子を同時に分析できる（低処理量）が、手法が退屈で、時間がかかることである。さらに、両方法は、健康被害をもたらす放射性標識技術の使用を必要とする。
【０００５】
１９９２年にディファレンシャルディスプレイ（ＤＤ）技術が登場すると（Ｌｉａｎｇ  ＆  Ｐａｒｄｅｅ，１９９２，Ｓｃｉｅｎｃｅ，ｖｏｌ．２５７（５０７２）：９６７−７１）、ＤＤに関して多くの改変及び改良がなされ（例えば、Ｏｒｄｅｒｅｄ  Ｄｉｆｆｅｒｅｎｔｉａｌ  Ｄｉｓｐｌａｙ，Ｍａｔｚ  ｅｔ  ａｌ．，１９９７，Ｎｕｃｌ．Ａｃｉｄｓ．Ｒｅｓ．，ｖｏｌ．２５（１２）：２５４１−２）、多重転写産物の特徴づけを目標として、第一歩が踏み出された。ＤＤの特徴は、無作為に設計されたＰＣＲプライマーを、分析されるべきｃＤＮＡ試料に低い厳密度で徐冷することによって、遺伝子の無作為なサブセットが標的とされ、使用されるＰＣＲプライマーに対して高い相同性を有する配列を含有する発現された転写産物を優先的に増幅させることである。次に、増幅産物が配列ゲル上で分離され、転写された遺伝子のサブセットを表す指紋パターンを生じる。ＤＤ法は、ノザンブロット及びＲＮＡｓｅ保護アッセイと比較して処理量は高いが、ＤＤ法の制約は、これらの技術の再現性／頑強性がかなり低いことである。これは、使用されるランダムＰＣＲプライマーの非特異的な徐冷が１つの原因である。その結果、異なるランダムプライマーを用いて生成された指紋パターンは、転写産物の異なる（相補的な）サブセットを系統的に標的化しない。さらなる欠点は、ＤＤ法が、スラブゲルの調製又は毛細管ゲル電気泳動による検出を必要とすることである。さらに別の制約は、指紋中の観察されたバンドの遺伝子源が不明であり、遺伝資源を明らかにするために、バンドの切り出し、溶出、再増幅及びＤＮＡ配列決定を必要とすることである。後者の制約は、指紋をベースとした他の転写産物特徴づけ法と共通である。最後に、ゲル上のレーン／毛細管トレースあたり５０から１００個の断片を検出する本技術は、中程度に多重である。
【０００６】
ｃＤＮＡ−ＡＦＬＰ法（Ｂａｃｈｅｍ  ｅｔ  ａｌ．，１９９６，Ｐｌａｎｔ  Ｊ．，ｖｏｌ．９（５）：７４５−５３）は、ＤＤ技術の主要な制約の２つ、すなわち再現性／頑強性及び異なるＰＣＲプライマーを用いて生成された指紋中に得られる情報の相補性に対処する。選択的ＡＦＬＰ（登録商標）（Ｋｅｙｇｅｎｅ  Ｎ．Ｖ．，ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ；例えば欧州特許第０，５３４，８５８号及びＶｏｓ  Ｐ．，ｅｔ  ａｌ．（１９９５）．ＡＦＬＰ：ａ  ｎｅｗ  ｔｅｃｈｎｉｑｕｅ  ｆｏｒ  ＤＮＡ  ｆｉｎｇｅｒｐｒｉｎｔｉｎｇ．Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．２３，Ｎｏ．２１，ｐ．４４０７−４４１４参照）プライマーを用いた、アダプターが連結された制限断片の増幅は、厳密度の高い条件下で起こり再現性が高い指紋パターンが得られるので、ｃＤＮＡ−ＡＦＬＰ法の頑強性及び再現性は非常に高い。さらに、異なる選択的ヌクレオチドを有する選択的ＡＦＬＰプライマーを使用することによって、相補的情報を含有する指紋が確実に得られる。それゆえ、ｃＤＮＡ−ＡＦＬＰ技術によって、トランスクリプトームのサブセットのサンプリングが再現可能となる。（ｃＤＮＡ−）ＡＦＬＰ（及びＤＤ）の別の利点は、配列情報が事前に必要ではなく、それゆえ生物の広範囲に本技術を適用できることである。ｃＤＮＡ−ＡＦＬＰの制約は、レーン／トレースあたりの多重レベルが中程度であること、及びバンドの遺伝子源が直接には不明であるという事実である（ＤＤも参照）。
【０００７】
上述の転写産物特徴付け方法の多重レベルにおける制約は、ＳＡＧＥ（遺伝子発現の連続分析；Ｖｅｌｃｕｌｅｓｃｕ  ｅｔ  ａｌ．，１９９５，Ｓｃｉｅｎｃｅ，ｖｏｌ．２７０（５２３５）：４８４−７）及び大量並行サイン配列決定（ＭＰＳＳ（Ｍａｓｓｉｖｅｌｙ  Ｐａｒａｌｌｅｌ  Ｓｉｇｎａｔｕｒｅ  Ｓｅｑｕｅｎｃｉｎｇ）：Ｂｒｅｎｎｅｒ  ｅｔ  ａｌ．，２０００，Ｎａｔｕｒｅ  Ｂｉｏｔｅｃｈｎｏｌｏｇｙ，ｖｏｌ．１８（６）：６３０−４；Ｍｅｙｅｒｓ  ｅｔ  ａｌ．，２００４，Ｎａｔｕｒｅ  Ｂｉｏｔｅｃｈｎｏｌｏｇｙ，ｖｏｌ．２２（８）：１００６−１１）の両者によって対処されてきた。ｃＤＮＡ−ＡＦＬＰと同様、両方法は、試料ｃＤＮＡを切り出した後にアダプターを連結するために、ＩＩＳ型制限酵素を使用する。
【０００８】
ＳＡＧＥでは、アダプターが連結された断片は、その後、連鎖状になり、サンガー配列決定法によって配列決定される。サンガー配列トレースから、１４から２０ｂｐの短い配列タグが抽出され、転写された遺伝子に関する定量的情報を提供する（「デジタルノザン」）。試料間でタグの頻度を比較することによって、予め配列情報を必要とせずに、研究される試料間での相対的な発現レベルに関する情報が得られる。これにより異なる試料中での相対的な転写産物量が（正確に）決定されるが、得られた配列タグが短いことを考えると、研究される生物の巨大なＥＳＴ集合物又は全ゲノム配列が利用可能でなく、及び、タグ配列がＢＬＡＳＴ（Ｂａｓｉｃ  Ｌｏｃａｌ  Ａｌｉｇｎｍｅｎｔ  Ｓｅａｒｃｈ  Ｔｏｏｌ）分析等の相同性研究へ供することができなければ、タグがどの遺伝子に由来するかを評価することは困難である。それゆえ、ＳＡＧＥは、多重性が高く、再現性があり、頑強であるが、その価値は、配列決定されたゲノムを有する生物に限定される。別の制約は、前記方法が、大規模なサンガー配列決定に要する経費のため、多量の試料を処理することにあまり適していない（低処理量）ということである。
【０００９】
ＳＡＧＥとは対照的に、ＭＰＳＳは、固相配列決定反応に基づいている。しかしながら、ＭＰＳＳには、本質的に、ＳＡＧＥと同一の制約がある。すなわち、非常に短い配列タグ（約２０ｂｐ）が得られ、そのため、限定された（ゲノム）配列が利用することができる生物中の興味深い配列タグをさらに追跡する（遺伝子同定／アッセイ転換）ことが大きく妨げられる。要約すると、ＳＡＧＥ及びＭＰＳＳは、使用するために、配列情報を予め必要としない、頑強で、高度に多重化された転写産物特徴づけ技術であり、これらの価値は、全ゲノム配列が決定されている生物又は配列タグを遺伝子に接続するために、巨大なＥＳＴ集合物を利用できる生物に事実上限定される。両方法は、低処理量であり、技術的に複雑である。
【００１０】
概念的な長所は、両方法が、多くの遺伝子の（相対的な）転写レベルの非偏向概算値を同時に提供する正確な配列決定と組み合わせた（「デジタルノザン」を生じる、）転写産物ライブラリの統計的なサンプリングに依存すること、及び転写産物の特徴づけが、固体支持体上でのプローブとの交差ハイブリッド形成に弱点がないという事実である。
【００１１】
１９９５年に、転写産物特徴づけの分野にパラダイム変化を与えた遺伝子発現マイクロアレイが導入された（Ｓｃｈｅｎａ  ｅｔ  ａｌ．，１９９５，Ｓｃｉｅｎｃｅ，ｖｏｌ．２７０（５２３５）：４６７−７０）。初期には、プローブとしてＥＳＴ由来のＰＣＲ産物を含有する、いわゆる「点状の」マイクロアレイが使用されたが、その後、頑強性が高く、拡大縮小が柔軟であるために、焦点は、オリゴヌクレオチドＤＮＡチップに移った（Ｐｅａｓｅ  ｅｔ  ａｌ．，１９９４，Ｐｒｏｃ．Ｎａｔ．Ａｃ．Ｓｃｉ．ＵＳＡ，ｖｏｌ．９１（１１）：５０２２−６）。現在、転写産物の特徴づけの市場は、多様な供給元（例えば、Ａｆｆｙｍｅｔｒｉｘ、Ｎｉｍｂｌｅｇｅｎ、Ａｇｉｌｅｎｔ等）から得られるオリゴヌクレオチドＤＮＡチップが支配的である。ＤＮＡチップの威力は、その表面に付着／合成できるＤＮＡ配列数が多いことにあり、これによって、転写産物の特徴づけを大量に並行処理することが可能となり、例えば、既知の全てのヒト遺伝子（＝高い遺伝子多重化レベル）に対して転写産物の特徴づけを実施することができる。さらに、チップの製造過程及びハイブリッド形成は、自動化及び調節が可能であり、これにより、それぞれ高処理量及び頑強性が可能となる。その結果、ＤＮＡチップは、２００５年において、転写産物の特徴づけに関する最先端技術である。しかしながら、多重化能、処理量及び頑強性は、ＤＮＡチップの非常に重要な利点であるが、チップを基礎とする転写産物特徴づけに関する２つの重要な制約は、チップを構築できるようにするために配列情報が必要とされること、及び、複製された遺伝子ファミリーのメンバー由来のものなど、相同性の高い配列間の交差ハイブリッド形成が、結果の正確性に影響を及ぼし得ることである。ハイブリッド形成をベースとした検出の本質的な特徴であるあるので、後者は、モニターし／排除することが非常に困難である。これらの事実により、（異なる製造技術及び応用プロトコールを反映する可能性を秘めた）異なる供給元から得られたＤＮＡチップを使用して得られた結果の比較は、実施が困難である（Ｙａｕｋ  ｅｔ  ａｌ．，２００５，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．３２（１５）：ｅ１２４）。１つのプラットフォーム内で、リアルタイムＰＣＲアッセイ（例えばＴａｑＭａｎ、Ｉｎｖａｄｅｒ）等の独立した方法による結果の検証が必要とされる。従って、同一プラットフォームが全ての試料に関して使用される場合、ＤＮＡチップは、デジタルノザンの概念に適合するデータを提供しないが、相対的な発現レベルを決定するのに有用である。
【発明の開示】
【発明が解決しようとする課題】
【００１２】
理想的には、転写産物の特徴づけ技術は、高度に多重化されており、すなわち、多くの遺伝子が同時に研究され、高処理量であり、非常に頑強で、再現性があり、非常に正確で（交差ハイブリッド形成という弱点がない。）、あらかじめ配列情報を必要とせずに適用可能であり得る。以下に記載される本発明は、このような基準に適合する方法を提供する。
【課題を解決するための手段】
【００１３】
発明の要旨
本発明者らは、異なる戦略を使用して、この問題を解決できること、及び転写産物の特徴づけにおいて、高処理量配列決定技術を効率的に使用できることをここに発見した。
【００１４】
本発明は、好ましくは、トランスクリプトーム（ｔｒａｎｓｃｒｉｐｔｏｍｅ）を再現性のあるサブセットに分割する技術を使用することを含む。これらのサブセットは、配列決定され、個々の転写産物に相当する連結断片へと構築される。異なる再現可能なサブセットが提供されるようにこの段階を反復することによって、連結断片の異なるセットが得られる。これらの異なる連結断片は、転写産物のドラフト配列を構築するために使用される。本発明は、配列に関する知識を一切必要とせず、いかなる複雑性の転写産物に対しても適用できる。本発明は、例えば、同一の生物又は別異の生物の異なる組織に由来する転写産物の組み合わせにも適用することができる。本発明は、目的のあらゆる転写産物へのアクセスを、より迅速で、信頼性が高く、且つより高速なものとし、それにより、転写産物の分析を促進する。
【００１５】
本発明は、これらの遺伝子の配列情報を必要とせずに、遺伝子の相対的な転写産物レベルを（偏向せずに）決定することにも関する。この目的のために、ｃＤＮＡ試料内での配列の頻度は、前記ｃＤＮＡ試料の複雑性が低下したライブラリを配列決定し、及び配列がライブラリ中で同定される回数を決定するために配列を整列することによって決定される。これは、第二のｃＤＮＡ試料に対して反復され得、２つのｃＤＮＡ試料の頻度は、必要であれば相対的な転写レベルを決定するために標準化され、比較され得る。
【００１６】
定義
以下の記載及び例において、多くの用語が使用される。このような用語に付与されるべき範囲を含め、明細書及び特許請求の範囲を明確且つ一貫して理解するために、以下の定義が提供される。本明細書に別段の定義がなければ、使用される全ての技術用語及び科学用語は、本発明が属する分野の当業者によって一般的に理解されるのと同一の意味を有する。全ての公報、特許出願、特許及び他の参考文献の開示は、参照により、その全ての内容が全体として本明細書に組み込まれる。
【００１７】
核酸：本発明に係る核酸には、ピリミジン及びプリン塩基、好ましくは、シトシン、チミン及びウラシル、並びに、それぞれアデニン及びグアニンの全てのポリマー又はオリゴマーが含まれ得る（Ａｌｂｅｒｔ  Ｌ．Ｌｅｈｎｉｎｇｅｒ，Ｐｒｉｎｃｉｐｌｅｓ  ｏｆ  Ｂｉｏｃｈｅｍｉｓｔｒｙ，ａｔ  ７９３−８００（Ｗｏｒｔｈ  Ｐｕｂ．１９８２）（全ての目的のために、その全ての内容が、参照により、本明細書に組み込まれる））。本発明は、全てのデオキシリボヌクレオチド、リボヌクレオチド又はペプチド核酸成分、及びこれらの塩基のメチル化、ヒドロキシメチル化又はグリコシル化された形態等のあらゆる化学的変種等を想定する。ポリマー又はオリゴマーの組成は、不均一又は均一であり得、天然源から単離され得、又は人工的に若しくは合成的に作製され得る。さらに、核酸は、ＤＮＡ若しくはＲＮＡ又はこれらの混合物であり得、ホモ二本鎖、ヘテロ二本鎖及びハイブリッド状態を含む一本鎖又は二本鎖の形態で、永続的に又は一過的に存在し得る。
【００１８】
複雑性の低下：複雑性の低下という用語は、試料のサブセットの生成によって、ゲノムＤＮＡ等の核酸試料の複雑性を低下させる方法を表すために使用される。このサブセットは、完全な（すなわち、複雑な）試料の代表であり得、好ましくは、再現性のあるサブセットである。本文脈において、再現性があるとは、同一方法を使用して、同一試料の複雑性を低下させる場合に、同一のサブセット又は少なくとも同程度のサブセットが得られることを意味する。複雑性の低下のために使用される方法は、本分野で公知の複雑性を低下させるための全ての方法であり得る。複雑性を低下させるための方法の非限定的な例には、ＡＦＬＰ（登録商標）（Ｋｅｙｇｅｎｅ  Ｎ．Ｖ．，ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ；例えば欧州特許第０，５３４，８５８号参照）、Ｄｏｎｇによって記載されている方法（例えば、国際公開第０３／０１２１１８号、国際公開第００／２４９３９号参照）、指標付き連結（Ｕｎｒａｕ，ｅｔ  ａｌ．，１９９４，Ｇｅｎｅ，１４５：１６３−１６９）、米国特許出願第２００５／２６０６２８号、国際公開第０３／０１０３２８号、米国特許出願第２００４／１０１５３号に開示されているもの、ゲノム分割（例えば、国際公開第２００４／０２２７５８号参照）、遺伝子発現の連続分析（ＳＡＧＥ；例えば、Ｖｅｌｃｕｌｅｓｃｕ  ｅｔ  ａｌ．，１９９５、上述参照、及びＭａｔｓｕｍｕｒａ  ｅｔ  ａｌ．，１９９９，Ｔｈｅ  Ｐｌａｎｔ  Ｊｏｕｒｎａｌ，ｖｏｌ．２０（６）：７１９−７２６）及びＳＡＧＥの変法（例えば、Ｐｏｗｅｌｌ，１９９８，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．２６（１４）：３４４５−３４４６；及びＫｅｎｚｅｌｍａｎｎ  ａｎｄ  Ｍｕｈｌｅｍａｎｎ，１９９９，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．２７（３）：９１７−９１８参照）、微量ＳＡＧＥ（例えば、Ｄａｔｓｏｎ  ｅｔ  ａｌ．，１９９９，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．２７（５）：１３００−１３０７参照）、大量並行サイン配列決定（ＭＰＳＳ；例えば、Ｂｒｅｎｎｅｒ  ｅｔ  ａｌ．，２０００，Ｎａｔｕｒｅ  Ｂｉｏｔｅｃｈｎｏｌｏｇｙ，ｖｏｌ．１８：６３０−６３４及びＢｒｅｎｎｅｒ  ｅｔ  ａｌ．，２０００，ＰＮＡＳ，ｖｏｌ．９７（４）：１６６５−１６７０参照）、自己減算型ｃＤＮＡライブラリ（Ｌａｖｅｄｅｒ  ｅｔ  ａｌ．，２００２，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．３０（９）：ｅ３８）、リアルタイム多重連結依存的プローブ増幅（ＲＴ−ＭＬＰＡ；例えば、Ｅｌｄｅｒｉｎｇ  ｅｔ  ａｌ．，２００３，ｖｏｌ．３１（２３）：ｅ１５３参照）、高適用範囲発現特徴づけ（ＨｉＣＥＰ；例えば、Ｆｕｋｕｍｕｒａ  ｅｔ  ａｌ．，２００３，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．３１（１６）：ｅ９４）、Ｒｏｔｈ  ｅｔ  ａｌ．，Ｎａｔｕｒｅ  Ｂｉｏｔｅｃｈｎｏｌｏｇｙ，ｖｏｌ．２２（４）：４１８−４２６に開示されている普遍的マイクロアレイシステム、トランスクリプトーム減算法（例えば、Ｌｉ  ｅｔ  ａｌ．，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．３３（１６）：ｅ１３６参照）、及び断片ディスプレイ（例えば、��ｅｔｓｉｓ  ｅｔ  ａｌ．，２００４，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．３２（１６）：ｅ１２７参照）が含まれる。本発明において使用される複雑性低下法は、再現性があることで共通している。同一の試料の複雑性を同じ様式で低下させる場合に、より無作為な複雑性の低下（顕微解剖又は選択された組織中で転写されたゲノムの一部に相当するｍＲＮＡ（ｃＤＮＡ）の使用など）とは異なり、試料の同一のサブセットが得られるという意味において再現性があり、その再現性に関しては、組織の選択、単離時間などに依存する。
【００１９】
タグ付け：タグ付けという用語は、核酸試料を第二の又はさらなる核酸試料と識別できるようにするために、核酸試料へタグを付加することを指す。タグ付けは、例えば、複雑性を低下させる間に配列同定子を添加することによって、又は本分野で公知のその他の手段によって、実施され得る。このような配列同定子は、様々な長さの、但し、特異的な核酸試料を同定するために一意的に使用される確定された長さの特有の塩基配列とすることができる。その典型的な例は、例えば、ＺＩＰ配列である。このようなタグを使用して、さらなる加工を行って、試料の起源を決定することができる。異なる核酸試料に由来する加工された生成物を組み合わせるには、異なる核酸試料は、異なるタグを使用して同定すべきである。
【００２０】
タグ付きライブラリ：タグ付きライブラリという用語は、タグ付き核酸のライブラリを指す。
【００２１】
配列決定：  配列決定という用語は、核酸試料、例えば、ＤＮＡ又はＲＮＡ中でヌクレオチド（塩基配列）の順序を決定することを指す。
【００２２】
整列させること及び整列化：「整列させること」及び「整列化」という用語は、同一又は類似のヌクレオチドの短い又は長い伸長の存在に基づいて、２つ又はそれ以上のヌクレオチド配列を比較することを意味する。以下にさらに説明されているように、ヌクレオチド配列を整列化するための幾つかの方法が、本分野で公知である。「構築すること」又は「クラスター形成すること」という用語が、同義語として使用される場合があるが、これらの用語は、技術的には同一ではない。整列化は、最大の相同性を比較することに基づいて行われるのに対して、構築とは、重複に基づいて連結断片を調製することを意味する。
【００２３】
高処理量スクリーニング：しばしばＨＴＳと略記される高処理量スクリーニングは、生物学及び化学の分野と特に関連する科学実験のための方法である。現代のロボット工学と他の特殊化された研究室用ハードウェアとの組み合わせを通じて、研究者は、高処理量スクリーニングによって、多量の試料を同時に効果的にスクリーニングする仔おtが可能となる。
【００２４】
高処理量配列決定：高処理量技術を使用してヌクレオチド配列の配列を決定すること。
【００２５】
制限エンドヌクレアーゼ：制限エンドヌクレアーゼ又は制限酵素は、二本鎖ＤＮＡ分子中の特異的なヌクレオチド配列（標的部位）を認識する酵素であり、あらゆる標的部位において、ＤＮＡ分子の両鎖を切断する。
【００２６】
制限断片：制限エンドヌクレアーゼによる消化によって生成されるＤＮＡ分子は、制限断片と呼ばれる。あらゆる所与のゲノム（又は核酸（その由来に関わらない。））は、特定の制限エンドヌクレアーゼによって制限断片の分離したセットへと消化される。制限エンドヌクレアーゼの切断によって生じるＤＮＡ断片は、様々な技術においてさらに使用することが可能であり、例えば、ゲル電気泳動によって検出することができる。
【００２７】
ゲル電気泳動：制限断片を検出するために、大きさに基づいて二本鎖ＤＮＡ分子を分画するための分析方法が必要とされ得る。このような分画を達成するために最も一般的に使用される技術は、（毛細管）ゲル電気泳動である。ＤＮＡ断片がこのようなゲル中で移動する速度は、ＤＮＡ断片の分子量に依存し、従って、移動する距離は、断片長が増大するにつれて減少する。ゲル電気泳動によって分画されるＤＮＡ断片は、パターン中に含まれる断片の数が十分小さければ、染色手法、例えば銀染色又は臭化エチジウムを用いた染色によって直接可視化することができる。あるいは、ＤＮＡ断片のさらなる処理によって、フルオロフォア又は放射性標識等の、検出可能な標識を断片中に取り込ませ得る。
【００２８】
連結：二本鎖ＤＮＡ分子が互いに共有結合される、リガーゼ酵素によって触媒される酵素反応は、連結と呼ばれる。一般に、両ＤＮＡ鎖は、互いに共有結合しているが、鎖の末端の１つの化学的若しくは酵素的修飾を通じて２つの鎖のうちの１つが連結しないようにすることも可能である。その場合、共有結合は、２つのＤＮＡ鎖のうちの１つのみに生じる。
【００２９】
合成オリゴヌクレオチド：化学的に合成できる好ましくは約１０から約５０までの塩基を有する一本鎖ＤＮＡ分子は、合成オリゴヌクレオチドと呼ばれる。一般に、これらの合成ＤＮＡ分子は、固有のヌクレオチド配列又は望ましいヌクレオチド配列を有するように設計されるが、関連配列を有し、ヌクレオチド配列内の特定の位置で異なるヌクレオチド組成を有する分子のファミリーを合成することが可能である。合成オリゴヌクレオチドという用語は、設計された又は望ましいヌクレオチド配列を有するＤＮＡ分子を指すために使用される。
【００３０】
アダプター：塩基対の限定された数、例えば、約１０から約３０塩基対の長さの短い二本鎖ＤＮＡ分子であり、制限断片の末端に連結できるように設計されている。アダプターは、一般的に、互いに部分的に相補的であるヌクレオチド配列を有する２つの合成オリゴヌクレオチドからなる。２つの合成オリゴヌクレオチドを、溶液中で、適切な条件下で混合すると、互いに徐冷して二本鎖構造を形成する。徐冷後に、アダプター分子の一方の末端は、制限断片の末端に対して適合的であり、一方の末端に連結できるように設計することが可能であり、アダプターの他方の末端は、連結できないように設計することが可能であるが、必ずしも、このようにする必要はない（二重連結アダプター）。
【００３１】
アダプターが連結された制限断片：連結の結果として、アダプターによってキャップ形成された制限断片。
【００３２】
プライマー：一般に、プライマーという用語は、ＤＮＡの合成を開始させることができるＤＮＡ鎖を指す。ＤＮＡポリメラーゼは、プライマーなしに、ＤＮＡを新規に合成することはできない。ＤＮＡポリメラーゼは、構築されるべきヌクレオチドの順序にを指定するために、相補鎖がテンプレートとして使用される反応において、既存のＤＮＡ鎖を伸長することができるにすぎない。本発明者らは、ポリメラーゼ連鎖反応（ＰＣＲ）において使用される合成オリゴヌクレオチド分子をプライマーと称する。
【００３３】
ＤＮＡ増幅：ＤＮＡ増幅という用語は、典型的には、ＰＣＲを使用して、二本鎖ＤＮＡ分子をインビトロで合成することを指すために使用される。他の増幅方法が存在すること、及び、本旨から逸脱せずに、本発明において、他の増幅方法を使用し得ることが注目される。
【００３４】
発明の詳細な説明
本発明は、
（ａ）ｃＤＮＡを準備する段階；
（ｂ）前記ｃＤＮＡの少なくとも１つに対して複雑性の低下を実施してｃＤＮＡ断片を含むｃＤＮＡの第一のライブラリを得る段階；
（ｃ）高処理量配列決定によって、第一のライブラリの前記ｃＤＮＡ断片のヌクレオチド配列の少なくとも一部を決定する段階；
（ｄ）段階ｄ）の第一のライブラリの前記ｃＤＮＡ断片の前記ヌクレオチド配列を整列させて第一のライブラリの連結断片を生成する段階；及び
（ｅ）前記ｃＤＮＡの前記ヌクレオチド配列を決定する段階
を含む、ｃＤＮＡのヌクレオチド配列を決定するための方法を提供する。
【００３５】
これまで、配列決定技術の分野において、転写産物を示すための、ｃＤＮＡの高処理量配列決定と組み合わせた前記複雑性の低下の使用は、開示されておらず、示唆されてもいない。
【００３６】
方法の段階（ａ）において、ｃＤＮＡが準備される。ｃＤＮＡの調製方法は、本分野で周知である。調製のための方法が、以下に示されている。しかしながら、ｃＤＮＡの調製のためのあらゆる方法が使用され得る。
【００３７】
ｃＤＮＡ（相補的ＤＮＡ）は、通常、逆転写酵素を使用してｍＲＮＡから調製される。この場合、逆転写酵素は、ＲＮＡに対して塩基対形成され、遊離の３’−Ｏｈ基を含有するプライマーと共に準備されると、ＲＮＡテンプレートに相補的なＤＮＡ鎖を合成する。このようなプライマーは、例えば、多くの真核生物ｍＲＮＡ分子の３’末端に存在するポリＡ配列と対形成するオリゴｄＴプライマーであり得る。次に、４つのデオキシリボヌクレオシド三リン酸塩の存在下で、ｃＤＮＡ鎖の残りを合成できる。得られたＲＮＡ−ＤＮＡハイブリッドのＲＮＡ鎖は、その後、例えば、ｐＨを上昇させることによって加水分解される。ＲＮＡとは異なり、ＤＮＡは、アルカリ加水分解耐性であり、そのため、ＤＮＡ鎖は未変化のままである。これに代わるプライマーは、ランダムプライマーであり得る。ｃＤＮＡの無作為な合成開始は、逆転写酵素がｍＲＮＡテンプレートを完全に転写することができない場合に、又は二次構造が存在する場合に、利点があり得る。さらに、これに代わるプライマーは、配列特異的プライマーであり得る。
【００３８】
生物体の組織の細胞又は生物自体からＲＮＡを単離する方法は、分子生物学の分野で周知である。さらに、ｃＤＮＡ合成のための多くの市販のキットが、例えば、ＡＢｇｅｎｅ、Ａｍｂｉｏｎ、Ａｐｐｌｉｅｄ  Ｂｉｏｓｙｓｔｅｍｓ、ＢｉｏＣｈａｉｎ、Ｂｉｏ−Ｒａｄ、Ｃｌｏｎｔｅｃｈ、ＧＥ  Ｈｅａｌｔｈｃａｒｅ、ＧｅｎｅＣｈｏｉｃｅ、Ｉｎｖｉｔｒｏｇｅｎ、Ｎｏｖａｇｅｎ、Ｑｉａｇｅｎ、Ｒｏｃｈｅ  Ａｐｐｌｉｅｄ  Ｓｃｉｅｎｃｅ、Ｓｔｒａｔａｇｅｎｅ等から購入できる。このような方法は、例えば、Ｓａｍｂｒｏｏｋ  ｅｔ  ａｌ．（Ｓａｍｂｒｏｏｋ，Ｊ．，Ｆｒｉｔｓｃｈ，Ｅ．Ｆ．，ａｎｄ  Ｍａｎｉａｔｉｓ，Ｔ．，ｉｎ  Ｍｏｌｅｃｕｌａｒ  Ｃｌｏｎｉｎｇ：Ａ  Ｌａｂｏｒａｔｏｒｙ  Ｍａｎｕａｌ．Ｃｏｌｄ  Ｓｐｒｉｎｇ  Ｈａｒｂｏｒ  Ｌａｂｏｒａｔｏｒｙ  Ｐｒｅｓｓ，ＮＹ，Ｖｏｌ．１，２，３（１９８９））に記載されている。ＲＮＡは、細胞培養物、組織等の幾つもの源から単離され得る。
【００３９】
本発明に係る方法の段階（ｂ）において、複雑性の低下は、ｃＤＮＡ断片を含むｃＤＮＡの第一のライブラリを得るために、ｃＤＮＡの少なくとも一部に対して実施される。複雑性の低下のための多くの方法は、定義の節で示されるように、本分野において公知である。
【００４０】
本発明の一実施形態において、核酸試料の複雑性を低下させる段階は、制限断片中で核酸試料を酵素で切断すること、制限断片を分離すること及び制限断片の特定のプールを選択することを含む。場合により、選択された断片は次に、ＰＣＲプライマーテンプレート／結合配列を含有するアダプター配列に連結される。
【００４１】
複雑性の低下に関する一実施形態において、ＩＩｓ型エンドヌクレアーゼは、核酸試料を消化するために使用され、制限断片は、アダプター配列へ選択的に連結される。アダプター配列は、連結されるべき突出中に多様なヌクレオチドを含有することができ、突出中のヌクレオチドの適合するセットを有するアダプターのみが、断片に連結され、その後増幅される。本技術は、「標識付けするリンカー」として本分野で示されている。本原理の例は、とりわけ、Ｕｎｒａｕ  ａｎｄ  Ｄｅｕｇａｕ（１９９４）Ｇｅｎｅ  １４５：１６３−１６９中に見出され得る。
【００４２】
一実施形態において、複雑性を低下させるための方法は、異なる標的部位及び頻度を有する２つの制限エンドヌクレアーゼ並びに２つの異なるアダプター配列を利用して、ＡＦＬＰ等においてアダプターが連結された制限断片を準備する。
【００４３】
本発明の一実施形態において、複雑性を低下させる段階は、任意に開始されるＰＣＲ（Ａｒｂｉｔｒａｒｉｌｙ  Ｐｒｉｍｅｄ  ＰＣＲ）を試料に対して実施することを含む。
【００４４】
本発明の一実施形態において、複雑性を低下させる段階は、ＤＮＡを変性及び再徐冷した後、二本鎖を除去することによって、反復された配列を除去することを含む。
【００４５】
本発明のある種の実施形態において、複雑性の低下の段階は、望ましい配列を含有するオリゴヌクレオチドプローブに結合されている磁気ビーズに、核酸試料をハイブリッド形成することを含む。本実施形態はさらに、ハイブリッド形成された試料を一本鎖ＤＮＡヌクレアーゼへ曝露して、一本鎖ＤＮＡを除去すること、ＩＩｓクラスの制限酵素を含有するアダプター配列を連結して磁気ビーズを放出することを含み得る。本実施形態は、単離されたＤＮＡ配列の増幅を含んでもよく、又は含まなくてもよい。さらに、アダプター配列は、ＰＣＲオリゴヌクレオチドプライマーのためのテンプレートとして使用されてもよく、又は使用しなくてもよい。本実施形態において、アダプター配列は、配列同定子又はタグを含有してもよく、又は含有しなくてもよい。
【００４６】
本発明の特定の実施形態において、複雑性の低下は、ディファレンシャルディスプレイ技術又はＲＥＡＤＳ（Ｇｅｎｅ  Ｌｏｇｉｃ）技術を利用する。
【００４７】
本発明の特定の実施形態において、複雑性を低下させる方法は、ＤＮＡ試料を不適正塩基対の結合タンパク質へ曝露すること、及び試料を３’から５’のエキソヌクレアーゼで消化した後、一本鎖ヌクレアーゼで消化することを含む。本実施形態は、不適正塩基対結合タンパク質に付着した磁気ビーズの使用を含んでもよく、又は含まなくてもよい。
【００４８】
本発明の一実施形態において、複雑性の低下は、本明細書の別の箇所に記載されているＣＨＩＰ法、又はＳＳＲ、ＮＢＳ領域（ヌクレオチド結合領域）、プロモーター／エンハンサー配列、テロメア共通配列、ＭＡＤＳボックス遺伝子、ＡＴＰａｓｅ遺伝子ファミリー及び他の遺伝子ファミリー等の保存モチーフに対して配向するＰＣＲプライマーの設計を含む。
【００４９】
段階（ｃ）において、第一のライブラリのｃＤＮＡ断片のヌクレオチド配列の少なくとも一部が、高処理量配列決定によって決定される。高処理量配列決定法の非限定的な例は、国際公開第０３／００４６９０号、国際公開第０３／０５４１４２号、国際公開２００４／０６９８４９号、国際公開第２００４／０７０００５号、国際公開第２００４／０７０００７号及び国際公開第２００５／００３３７５号において開示される方法、Ｓｅｏ  ｅｔ  ａｌ．（２００４）Ｐｒｏｃ．Ｎａｔｌ．Ａｃａｄ．Ｓｃｉ．ＵＳＡ  １０１：５４８８−９３による方法、及びＨｅｌｉｏｓ，Ｓｏｌｅｘａ，ＵＳ  Ｇｅｎｏｍｉｃｓ等の技術であり、これらは、参照により本明細書に組み込まれる。配列決定は、国際公開第０３／００４６９０号、国際公開第０３／０５４１４２号、国際公開２００４／０６９８４９号、国際公開第２００４／０７０００５号、国際公開第２００４／０７０００７号及び国際公開第２００５／００３３７５号（参照により本明細書に組み込まれる。）（全て、４５４法人の名においてである。）において開示される装置及び／又は方法を使用して実施されるのが最も好ましい。記載されている技術によって、単一試行で４００万個の塩基の配列決定が可能であり、サンガー配列決定に基づいた競合する技術及びＭｅｇａＢＡＣＥ（ＧＥ  Ｈｅａｌｔｈｃａｒｅ）若しくはＡＢＩ３７００（ｘ１）（Ａｐｐｌｉｅｄ  Ｂｉｏｓｙｓｔｅｍｓ）よりも１００倍高速で安価である。配列決定技術は、大まかに、４段階、すなわち：１）ＤＮＡの断片化及び一本鎖ＤＮＡ（ｓｓＤＮＡ）のライブラリへの特異的アダプターが連結；２）ｓｓＤＮＡのビーズへの徐冷及び油中水微小反応器中でのビーズの乳化；３）ＰｉｃｏＴｉｔｅｒＰｌａｔｅ（登録商標）中へのＤＮＡ運搬ビーズの配置；並びに４）ピロリン酸塩光信号の発生による多重ウェル中での同時配列決定からなる。本方法は、後により詳細に説明される。
【００５０】
段階（ｄ）において、第一のライブラリの連結断片を生成させるために、段階（ｄ）の第一のライブラリのｃＤＮＡ断片のヌクレオチド配列が整列される。
【００５１】
配列から連結断片を構築することによって、集合段階は、計算上あまり複雑ではなく、それゆえ、より迅速に実施される。ライブラリ中で配列を整列させることによって、制限断片のセットの各制限断片に対する連結断片は、各プライマーの組み合わせに対して構築できる。これにより、連結断片のセットが生じ、各セットは、特定の制限断片に対応する。その結果、少なくとも１つの制限エンドヌクレアーゼによるｃＤＮＡの制限から得られる各断片は、この段階で、決定された（連結断片の）配列を有する。
【００５２】
比較を目的とする配列の整列化の方法は、本分野において周知である。非限定的な多様なプログラム及び整列化アルゴリズムが、Ｓｍｉｔｈ  ａｎｄ  Ｗａｔｅｒｍａｎ（１９８１）Ａｄｖ．Ａｐｐｌ．Ｍａｔｈ．２：４８２；Ｎｅｅｄｌｅｍａｎ  ａｎｄ  Ｗｕｎｓｃｈ（１９７０）Ｊ．Ｍｏｌ．Ｂｉｏｌ．４８：４４３；Ｐｅａｒｓｏｎ  ａｎｄ  Ｌｉｐｍａｎ（１９８８）Ｐｒｏｃ．Ｎａｔｌ．Ａｃａｄ．Ｓｃｉ．ＵＳＡ  ８５：２４４４；Ｈｉｇｇｉｎｓ  ａｎｄ  Ｓｈａｒｐ（１９８８）Ｇｅｎｅ  ７３：２３７−２４４；Ｈｉｇｇｉｎｇ  ａｎｄ  Ｓｈａｒｐ（１９８９）ＣＡＢＩＯＳ  ５：１５１−１５３；Ｃｏｒｐｅｔ  ｅｔ  ａｌ．（１９８８）Ｎｕｃｌ．Ａｃｉｄｓ  Ｒｅｓ．１６：１０８８１−９０；Ｈｕａｎｇ  ｅｔ  ａｌ．（１９９２）Ｃｏｍｐｕｔｅｒ  Ａｐｐｌ．ｉｎ  ｔｈｅ  Ｂｉｏｓｃｉ．８：１５５−６５；及びＰｅａｒｓｏｎ  ｅｔ  ａｌ．（１９９４）Ｍｅｔｈ．Ｍｏｌ．Ｂｉｏｌ．２４：３０７−３１に記載されており、それらは参照により本明細書に組み込まれる。Ａｌｔｓｃｈｕｌ  ｅｔ  ａｌ．（１９９４）Ｎａｔｕｒｅ  Ｇｅｎｅｔ．６：１１９−２９（参照により本明細書に組み込まれる。）は、配列整列化方法及び相同性算出に関する詳細な考えを示している。
【００５３】
ＮＣＢＩのＢａｓｉｃ  Ｌｏｃａｌ  Ａｌｉｇｎｍｅｎｔ  Ｓｅａｒｃｈ  Ｔｏｏｌ（ＢＬＡＳＴ）（Ａｌｔｓｃｈｕｌ  ｅｔ  ａｌ．，１９９０）は、Ｎａｔｉｏｎａｌ  Ｃｅｎｔｅｒ  ｆｏｒ  Ｂｉｏｌｏｇｉｃａｌ  Ｉｎｆｏｒｍａｔｉｏｎ（ＮＣＢＩ，Ｂｅｔｈｅｓｄａ，Ｍｄ．）及びインターネットを含む幾つかの供給源から利用可能であり、配列分析プログラムｂｌａｓｔｐ、ｂｌａｓｔｎ、ｂｌａｓｔｘ、ｔｂｌａｓｔｎ及びｔｂｌａｓｔｘと共に使用される。ＢＬＡＳＴは、「ｈｔｔｐ／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／ＢＬＡＳＴ／」においてアクセス可能である。本プログラムを使用して配列の同一性を決定する方法に関する記載は、「ｈｔｔｐ：／／ｗｗｗ．ｎｃｂｉ．ｎｌｍ．ｎｉｈ．ｇｏｖ／ＢＬＡＳＴ／ｂｌａｓｔ＿ｈｅｌｐ．ｈｔｍｌ」において利用可能である。マイクロサテライトマイニングにおいて、さらに応用され得る（Ｖａｒｓｈｎｅｙ  ｅｔ  ａｌ．（２００５）Ｔｒｅｎｄｓ  ｉｎ  Ｂｉｏｔｅｃｈｎ．２３（１）：４８−５５参照）。
【００５４】
一実施形態において、整列化は、アダプター／プライマーに関して整えられた配列データに対して、及び／又は再構築された制限酵素認識配列を使用して、すなわちｃＤＮＡ由来の断片から配列データのみを使用することによって実施される。典型的には、得られた配列データは、断片の起源（すなわち、どの試料由来か）を同定するために使用され、アダプター及び／又は同定子配列に由来する配列は、データから除去され、この整えられたセットに対して整列化が実施される。
【００５５】
段階（ｅ）において、ｃＤＮＡのヌクレオチド配列は、例えば配列の集合化によって決定される。
【００５６】
前記方法は、例えば、ｃＤＮＡ中又は前記ｃＤＮＡの複雑性の低下した画分中に存在する異なる配列の数を決定するために、又は特定の遺伝子の発現を発見するために有用である。
【００５７】
一実施形態において、段階（ａ）は、ｉ）生物試料を準備する段階、ｉｉ）前記生物試料から全ＲＮＡ又はｍＲＮＡを単離する段階、ｉｉｉ）全ＲＮＡ又はｍＲＮＡからｃＤＮＡを合成する段階を含む。
【００５８】
一実施形態において、高処理量配列決定は、ビーズ等の固体支持体上で実施され（例えば、国際公開第０３／００４６９０号、国際公開第０３／０５４１４２号、国際公開２００４／０６９８４９号、国際公開第２００４／０７０００５号、国際公開第２００４／０７０００７号及び国際公開第２００５／００３３７５号参照）（全て、４５４  Ｃｏｒｐｏｒａｔｉｏｎの名で公開されている。）、これらは、参照により本明細書に組み込まれる。このような配列決定方法は、同時に多くの試料を安価で効率的に配列決定するのに特に適している。
【００５９】
さらなる実施形態において、高処理量配列決定は、合成による配列決定、好ましくはピロシーケンスに基づいている。ピロシーケンスは、本分野で公知であり、とりわけ、ｗｗｗ．ｂｉｏｔａｇｅｂｉｏ．ｃｏｍ；ｗｗｗ．ｐｙｒｏｓｅｑｕｅｎｃｉｎｇ．ｃｏｍ／技術部分に記載されている。前記技術は、例えば、国際公開第０３／００４６９０号、国際公開第０３／０５４１４２号、国際公開２００４／０６９８４９号、国際公開第２００４／０７０００５号、国際公開第２００４／０７０００７号及び国際公開第２００５／００３３７５号（参照により本明細書に組み込まれる。）において開示される装置及び／又は方法（全て、４５４  Ｌｉｆｅ  Ｓｃｉｅｎｃｅｓの名で公開されている。）においてさらに応用され、これらは参照により本明細書に組み込まれる。前記技術は、高処理量配列決定に特に適した高速で、非常に再現性のある技術である。
【００６０】
好ましい実施形態において、高処理量配列決定は、
（ｃ１）配列決定アダプターを断片に連結する段階；
（ｃ２）、各々、単一の断片と徐冷するビーズに、配列決定アダプター連結断片を徐冷する段階；
（ｃ３）各々、単一ビーズを含む油中水微小反応器中でビーズを乳化する段階；
（ｃ４）ビーズの表面上で配列決定アダプター連結断片を増幅するために、乳濁液ＰＣＲを実施する段階；
（ｃ５）増幅された配列決定アダプター連結断片を含有するビーズを選択／濃縮する段階；
（ｃ６）各々、単一ビーズを含むウェル中にビーズを負荷する段階；及び
（ｃ７）ピロリン酸塩シグナルを発生させる段階
を含む。
【００６１】
段階ｃ１）において、配列決定アダプターは、ライブラリ内の断片に連結される。前記配列決定アダプターには、ビーズに徐冷するための少なくとも１つの「鍵」領域、配列決定プライマー領域及びＰＣＲプライマー領域が含まれる。従って、適宜に改変された断片が得られる。
【００６２】
段階ｃ２）において、配列決定アダプター連結断片をビーズ（各ビーズは、単一の断片と徐冷する。）に徐冷する。配列決定アダプター連結断片のプールにビーズを過剰量で添加し、ビーズの大部分にについて、ビーズ当たり単一の適応断片の徐冷が確実に起こるようにする（ポワソン分布）。
【００６３】
段階ｃ３）において、各々単一のビーズを含む油中水微小反応器中でビーズを乳化する。
【００６４】
段階ｃ４）において、乳濁液ＰＣＲを実施し、ビーズの表面上で配列決定アダプター連結断片を増幅する。ＰＣＲ試薬は、油中水微小反応器中に存在し、ＰＣＲ反応が前記微小反応器内で起きる。
【００６５】
段階ｃ５）において、増幅された配列決定アダプター連結断片を含有するビーズを選択／濃縮する。
【００６６】
段階ｃ６）において、各々、単一のビーズを含むウェル中にビーズを負荷する。ウェルは、好ましくは、断片を多量に同時配列決定することが可能であるＰｉｃｏＴｉｔｅｒ（商標）Ｐｌａｔｅの一部である。酵素搭載ビーズを添加した後、ピロシーケンスを使用して、断片の配列を決定する。
【００６７】
段階ｃ７）において、ピロリン酸塩信号を発生させる。連続的段階において、ＰｉｃｏＴｉｔｅｒ（商標）Ｐｌａｔｅ及びビーズ及びその中の酵素ビーズは、慣用の配列決定試薬の存在下において、異なるデオキシリボヌクレオチドへ供され、デオキシリボヌクレオチドが組み込まれると、光信号が発生し、記録される。正しいヌクレオチドの組み込みによって、本分野で公知の手段によって検出することができるピロシーケンス信号が発生する。
【００６８】
本発明に係る方法の好ましい実施形態において、複雑性の低下は、
ｉ）．少なくとも１つの制限エンドヌクレアーゼでｃＤＮＡを消化して制限断片へと断片化する段階；
ｉｉ）．制限断片の一方又は両方の末端と適合性のある１つの末端を有する少なくとも１つの二本鎖合成オリゴヌクレオチドアダプターと前記制限断片を連結してアダプターが連結された制限断片を生成する段階；
ｉｉｉ）．前記アダプターが連結された制限断片を、１つ又はそれ以上のオリゴヌクレオチドプライマーと、ハイブリッド形成条件下で接触させる段階（前記１つ又はそれ以上のオリゴヌクレオチドプライマーは、少なくとも１つのアダプターの一部と、及び制限エンドヌクレアーゼの認識配列の残りの部分の一部とに相補的なヌクレオチド配列部分を含むプライマー配列を有する。）；並びに
ｉｖ）．ハイブリッド形成された１つ又はそれ以上のオリゴヌクレオチドプライマーの伸長によって、前記アダプターが連結された制限断片を増幅する段階
を含む方法によって実施される。
【００６９】
複雑性を低下させるための上述の方法は、ＡＦＬＰ（登録商標）（Ｋｅｙｇｅｎｅ  Ｎ．Ｖ．，ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓ；例えば、その全ての内容が全体として参照により本明細書に組み込まれる、欧州特許第０，５３４，８５８号及びＶｏｓ  ｅｔ  ａｌ．（１９９５）．ＡＦＬＰ：ａ  ｎｅｗ  ｔｅｃｈｎｉｑｕｅ  ｆｏｒ  ＤＮＡ  ｆｉｎｇｅｒｐｒｉｎｔｉｎｇ，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．２３，ｎｏ．２１，４４０７−４４１４参照）とも呼ばれる。ＡＦＬＰは、複雑性を低下させるための非常に再現性のある方法であり、それゆえ、本発明に係る方法に特に適している。ＡＦＬＰは、選択的制限断片増幅のための方法である。ＡＦＬＰは、配列情報を事前に一切必要とせず、全ての開始ｃＤＮＡに対して実施することが可能である。
【００７０】
従って、ＡＦＬＰは、アダプター連結断片の再現性のあるサブセットを提供する。ＡＦＬＰ技術の１つの有用な変種は、選択的ヌクレオチドを使用せず（すなわち、＋０／＋０プライマー）、時に、リンカー−ＰＣＲと呼ばれる。これは、特に転写産物及びそこから得られるｃＤＮＡに非常に適した複雑性の低下も提供する。
【００７１】
段階ｉ）において、制限断片へと断片化するために、少なくとも１つの制限エンドヌクレアーゼでｃＤＮＡを消化する。特定の実施形態において、少なくとも２つの制限エンドヌクレアーゼが使用される。他の実施形態において、３つ又はそれ以上の制限エンドヌクレアーゼが使用できる。制限エンドヌクレアーゼは、頻度の高い切断剤であり得（すなわち、典型的には４及び５個の切断剤、すなわち、それぞれ４個又は５個のヌクレオチドの認識配列を有する制限エンドヌクレアーゼ）、又は頻度の低い切断剤であり得（すなわち、典型的には、それぞれ６個又はそれ以上のヌクレオチドの認識部位を有する。）、又はそれらの組み合わせであり得る。特定の実施形態において、頻度の低い又は頻度の高い切断剤の組み合わせが使用され得る。制限エンドヌクレアーゼは、認識配列の片側又は両側の何れかで、認識配列の外側のｃＤＮＡを切断するＩＩｓ型及びＩＩｓａ型を含む全ての種類であり得る。
【００７２】
段階ｉｉ）において、アダプターが連結された制限断片を生成するために、制限断片の一方又は両方の末端と適合性のある１つの末端を有する少なくとも１つの二本鎖合成オリゴヌクレオチドアダプターと前記制限断片を連結する。好ましくは、アダプターは、エンドヌクレアーゼ認識部位が、アダプターが連結の際に修復されないようになっている。例えば、段階ｉ）における２つ又はそれ以上の制限エンドヌクレアーゼを使用する場合、２つ又はそれ以上の異なるアダプターを使用することも可能である。本連結段階は、アダプターが連結された制限断片を生じる。アダプターは、平滑末端化でき、又は段階ｉ）において使用される制限エンドヌクレアーゼに依存して、突出を含有し得る。
【００７３】
特定の実施形態において、アダプターは、標識付けリンカーとして知られるアダプターのセットであり得る（Ｕｎｒａｕ，ｅｔ  ａｌ．，１９９４，Ｇｅｎｅ，１４５：１６３−１６９）。
【００７４】
段階ｉｉｉ）において、前記アダプターが連結された制限断片を、１つ又はそれ以上のオリゴヌクレオチドプライマーと、ハイブリッド形成条件下で接触させる。１つ又はそれ以上のオリゴヌクレオチドプライマーは、少なくとも１つのアダプターの一部と、及び制限エンドヌクレアーゼの認識配列の残りの部分の一部と相補的なヌクレオチド配列部分を含むプライマー配列を有する。
【００７５】
標準的なハイブリッド形成条件は、選択的ハイブリッド形成のための条件である。選択的ハイブリッド形成とは、厳密なハイブリッド形成条件下で、非標的核酸配列とのハイブリッド形成より検出可能に大きな程度（例えば、バックグラウンドより少なくとも２倍）まで、核酸配列を、指定された核酸標的配列とハイブリッド形成すること、及び非標的核酸を実質的に除外することに関する。「厳密な条件」又は「厳密なハイブリッド形成条件」という用語には、プローブが、他の配列（例えば、バックグラウンドより少なくとも２倍）より検出可能に大きな程度まで、その標的配列とハイブリッド形成する条件に対する表記が含まれる。厳密な条件は、配列依存的であり、異なる環境において異なる。ハイブリッド形成及び／又は洗浄の条件の厳密性を調節することによって、プローブと１００％相補的である標的配列を同定することができる（相同的探索）。あるいは、類似性のより低い程度が検出できる（不均一性探索）ように、厳密性の条件は、配列中に幾つかの誤塩基対形成が許容されるように調整することができる。一般的に、プローブは、長さ約１００ヌクレオチド未満であり、場合により、長さ５０以下、又は２５ヌクレオチドである。典型的には、厳密な条件は、塩濃度が、約１．５Ｍ未満のＮａイオンである条件であり、典型的には、ｐＨ７．０から８．３において、約０．０１から１．０ＭのＮａイオン濃度（又は他の塩）であり、温度は、短いプローブ（例えば、１０から５０ヌクレオチド）に関しては、少なくとも約３０℃であり、長いプローブ（例えば、５０超のヌクレオチド）に関しては、少なくとも約６０℃である。厳密な条件は、ホルムアミド等の脱安定化剤の添加によっても達成され得る。典型的な低い厳密度の条件には、３７℃で、３０から３５％ホルムアミド、１Ｍ  ＮａＣｌ、１％ＳＤＳ（ドデシル硫酸ナトリウム）の緩衝溶液によるハイブリッド形成、及び５０から５５℃での１×から２×ＳＳＣ（２０×ＳＳＣ＝３．０Ｍ  ＮａＣｌ／０．３Ｍクエン酸三ナトリウム）中での洗浄が含まれる。典型的な中程度の厳密度の条件には、３７℃で、４０から４５％ホルムアミド、１Ｍ  ＮａＣｌ、１％ＳＤＳ中でのハイブリッド形成、及び５５から６０℃で、０．５×から１×ＳＳＣ中での洗浄が含まれる。典型的な高い厳密度の条件には、３７℃で、５０％ホルムアミド、１Ｍ  ＮａＣｌ、１％ＳＤＳ中でのハイブリッド形成、及び６０から６５℃での０．１×ＳＳＣ中での洗浄が含まれる。特異性は典型的に、ハイブリッド形成後の洗浄の作用であり、決定的因子は、最終洗浄溶液のイオン強度及び温度である。ＤＮＡ−ＤＮＡハイブリッドに関し、Ｔｍは、Ｍｅｉｎｋｏｔｈ  ａｎｄ  Ｗａｈｌ，Ａｎａｌ．Ｂｉｏｃｈｅｍ．，１３８：２６７−２８４（１９８４）の式：Ｔｍ＝８１．５℃＋１６．６（ｌｏｇＭ）＋０．４１（％ＧＣ）−０．６１（％ホルムアミド）−５００／Ｌから概算でき、式中Ｍは、一価の陽イオンのモル濃度、％ＧＣは、ＤＮＡ中のグアノシン及びシトシンのヌクレオチドの％であり、％ホルムアミドは、ハイブリッド形成溶液中のホルムアミドの％であり、Ｌは、塩基対中のハイブリッドの長さである。Ｔｍは、相補的標的配列が、完全に適合したプローブとハイブリッド形成する（定義されたイオン強度及びｐＨのもとでの）温度である。Ｔｍは、不適正塩基対の各１％につき約１℃低下し、従って、Ｔｍ、ハイブリッド形成及び／又は洗浄条件は、望ましい同一性の配列とハイブリッド形成するように調整できる。例えば、９０％超の同一性を有する配列が探索される場合、Ｔｍは１０℃低下し得る。一般的に、厳密な条件は、定義されたイオン強度及びｐＨで、具体的な配列及びその相補対に関する融点（Ｔｍ）よりも約５℃低く選択される。しかしながら、非常に厳密な条件は、融点（Ｔｍ）よりも１、２、３又は４℃低いハイブリッド形成及び／又は洗浄を利用し得、中程度に厳密な条件は、融点（Ｔｍ）よりも６、７、８、９又は１０℃低いハイブリッド形成及び／又は洗浄を利用し得、低い厳密度の条件は、融点（Ｔｍ）よりも１１、１２、１３、１４、１５又は２０℃低いハイブリッド形成及び／又は洗浄を利用し得る。等式を使用して、ハイブリッド形成及び洗浄組成物、並びに望ましいＴｍを使用して、当業者は、ハイブリッド形成及び／又は洗浄溶液の厳密性における変動が、内在的に記載されていることを理解する。誤塩基対の望ましい程度が、４５℃未満（水溶液）又は３２℃（ホルムアミド溶液）のＴｍをもたらす場合、ＳＳＣ濃度を上昇させることが好ましく、これによりさらに高い温度が使用できる。核酸のハイブリッド形成に対する広範囲に及ぶ指針は、Ｔｊｊｓｓｅｎ，Ｌａｂｏｒａｔｏｒｙ  Ｔｅｃｈｎｉｑｕｅｓ  ｉｎ  Ｂｉｏｃｈｅｍｉｓｔｒｙ  ａｎｄ  Ｍｏｌｅｃｕｌａｒ  Ｂｉｏｌｏｇｙ−Ｈｙｂｒｉｄｉｚａｔｉｏｎ  ｗｉｔｈ  Ｎｕｃｌｅｉｃ  Ａｃｉｄ  Ｐｒｏｂｅｓ，Ｐａｒｔ  １，Ｃｈａｐｔｅｒ  ２  “Ｏｖｅｒｖｉｅｗ  ｏｆ  ｐｒｉｎｃｉｐｌｅｓ  ｏｆ  ｈｙｂｒｉｄｉｚａｔｉｏｎ  ａｎｄ  ｔｈｅ  ｓｔｒａｔｅｇｙ  ｏｆ  ｎｕｃｌｅｉｃ  ａｃｉｄ  ｐｒｏｂｅ  ａｓｓａｙｓ”，Ｅｌｓｅｖｉｅｒ，Ｎ．Ｙ．（１９９３）；及びＣｕｒｒｅｎｔ  Ｐｒｏｔｏｃｏｌｓ  ｉｎ  Ｍｏｌｅｃｕｌａｒ  Ｂｉｏｌｏｇｙ，Ｃｈａｐｔｅｒ  ２，Ａｕｓｕｂｅｌ，ｅｔ  ａｌ．，Ｅｄｓ．，Ｇｒｅｅｎｅ  Ｐｕｂｌｉｓｈｉｎｇ  ａｎｄ  Ｗｉｌｅｙ−Ｉｎｔｅｒｓｃｉｅｎｃｅ，Ｎｅｗ  Ｙｏｒｋ（１９９５）中に見出される。
【００７６】
２つ又はそれ以上の制限エンドヌクレアーゼが使用される場合、２つ又はそれ以上のオリゴヌクレオチドプライマーが、エンドヌクレアーゼの認識部位に依存して、段階ｉｉｉ）において使用される可能性がある。欧州特許第０，５３４，８５８号及びＶｏｓ  ｅｔ  ａｌ．（（１９９５）．ＡＦＬＰ：ａ  ｎｅｗ  ｔｅｃｈｎｉｑｕｅ  ｆｏｒ  ＤＮＡ  ｆｉｎｇｅｒｐｒｉｎｔｉｎｇ，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，ｖｏｌ．２３，ｎｏ．２１，  ４４０７−４４１４）においてさらに説明されるように、オリゴヌクレオチドプライマーは、少なくとも１つのアダプターに対して、並びに制限エンドヌクレアーゼの認識配列の残りの部分、及び場合により制限エンドヌクレアーゼの認識配列の残りの部分の一部に対して相補的なヌクレオチド配列部分を含むプライマー配列を有する。典型的には、認識配列の一部は、制限エンドヌクレアーゼによる配列の制限の後に残存しているその一部である。従って、要約すると、プライマーは、アダプターが連結された制限断片の既知の一部に対して少なくとも相補的である。
【００７７】
段階ｉｖ）において、ハイブリッド形成された１つ又はそれ以上のオリゴヌクレオチドプライマーの伸長によって、前記アダプターが連結された制限断片を増幅する増幅は、好ましくは、本分野で周知のＰＣＲを使用して実施される。
【００７８】
本発明の好ましい実施形態において、プライマーは、プライマー配列の３末端に選択された配列をさらに含み、１から１０個の選択的ヌクレオチドを含む前記選択された配列が、制限エンドヌクレアーゼの認識配列の残りの部分に直接隣接して配置される部分と相補的である。典型的には、認識配列の一部は、制限エンドヌクレアーゼによる配列の制限の後に残存しているその一部である。その３’末端で、プライマーは、好ましくは、選択された配列を含有する。選択された配列は、１から１０個のヌクレオチドの既に選択されたセット、好ましくは１から８個、好ましくは１から５個、より好ましくは１から３個の選択されたヌクレオチドを含む。典型的なプライマーは、（２つの選択的なヌクレオチド（ＡＣ）に関して）実例となる次の構造、すなわち「５’アダプター特異的領域−制限配列特異的領域−ＡＣ−３’」を有し得る。従って、この典型的なプライマーは、２つの選択的ヌクレオチドＡＣを含有し、アダプターが連結された制限断片の既知の部分の後に（すなわち、制限エンドヌクレアーゼの制限部位の残りの後に）、最初の２つのヌクレオチドとして相補的ＴＧを含有するアダプター連結断片を増幅するにすぎない。
【００７９】
ＡＦＬＰ、その利点、その実施形態、並びに技術、酵素、アダプター、プライマー及びさらなる化合物及びＡＦＬＰに使用されるツールに関するさらなる記載については、その全ての内容が、全体として、本明細書により組み込まれる米国特許第６，０４５，９９４号、欧州特許（Ｂ）第０，５３４，８５８号、欧州特許第９７６８３５号及び欧州特許第９７４６７２号、国際公開第０１／８８１８９号及びＶｏｓ  ｅｔ  ａｌ．Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓｅａｒｃｈ，１９９５，  ２３，  ４４０７−４４１４に表記される。
【００８０】
一実施形態において、前記アダプターは、同定子配列をさらに含む。このような同定子配列は、例えば、複雑性の低下によって得られるライブラリの起源を示すのに使用される変動する長さの特有の塩基配列であり得る。
【００８１】
本発明は、
（ａ）ｃＤＮＡを準備する段階；
（ｂ）前記ｃＤＮＡの少なくとも部分について複雑性の低下を実施してｃＤＮＡ断片を含む前記ｃＤＮＡの第一のライブラリを得る段階；
（ｃ）配列決定によって、第一のライブラリの前記ｃＤＮＡ断片のヌクレオチド配列の少なくとも一部を決定する段階；及び 
（ｄ）ヌクレオチド配列の頻度を決定する段階
を含む、ヌクレオチド配列の頻度を決定するための方法にも関する。
【００８２】
方法の段階（ａ）において、ｃＤＮＡが準備される。ｃＤＮＡを調製する方法は本分野で周知であり、適切な方法が上述に提供されている。ｃＤＮＡは、上述に説明されるように、全ての起源に由来し得る。
【００８３】
方法の段階（ｂ）において、複雑性の低下は、ｃＤＮＡ断片を含むｃＤＮＡの第一のライブラリを得るために、ｃＤＮＡの少なくとも一部に対して実施される。複雑性の低下は、上述に説明されるように、本分野で公知の何れかの方法によって実施され得る。
【００８４】
本発明に従った方法の段階（ｃ）において、第一のライブラリのｃＤＮＡ断片のヌクレオチド配列の少なくとも一部が、配列決定によって決定される。配列決定は、周知のサンガー（ジデオキシ）法を含む、本分野で公知の何れかの方法によって実施できる。好ましい実施形態において、配列決定は、複数の試料の同時配列決定が可能な高処理量配列決定を使用して実施される。高処理量配列決定のための好ましい方法は、上述に説明される。
【００８５】
本発明に係る方法の段階（ｄ）において、ヌクレオチド配列の頻度が決定される。ヌクレオチド配列の頻度は、例えば、次の方法によって決定され得る。ｃＤＮＡ断片のヌクレオチド配列の整列化は、同一の転写された遺伝子に由来するヌクレオチド配列を回収するために、及び前記ヌクレオチド配列を計数するために使用され得る。ヌクレオチド配列が、同一の転写された遺伝子に由来するかどうかは、配列間の相同性によってまだ確立されていない。本発明において、ヌクレオチド配列が、少なくとも１０の、好ましくは少なくとも１５の、より好ましくは少なくとも２０の、さらにより好ましくは少なくとも２５、３０、４０、５０、１００、１５０、２００のヌクレオチドの長さにわたって、少なくとも９５、９６、９７、９８、９９、１００％相同である場合に、前記ヌクレオチド配列が、同一の転写された遺伝子に由来すると仮定される。前記方法は、統計的に異なる頻度を示すために、Ｔ検定等の統計的解釈によって補助され得る。配列の同定された数に基づいて、単純な順位決定を実施することも可能である。試料１において、（不明な）遺伝子「Ｘ」のヌクレオチド配列が１０回測定され（１０は、例えば、９８％の配列相同性を有するヌクレオチド配列の数である。）、試料２において、同一の配列が２０回測定されると仮定する。この場合、試料２における遺伝子Ｘの転写レベルは、試料１の２倍であると思われる。但し、試料１及び２に対して、決定された配列の総数は同一であり、従って、転写産物の正確な特徴づけは、試料間の標準化を必要とし得、及び／又は遺伝子「Ｘ」に由来する配列の頻度を、相対的な転写レベルが複数の試料間で一定であると仮定されるいわゆるハウスキーピング遺伝子の頻度と比較し得る。試料の表現型特徴に関する試料間での相対的な転写特性の順位決定は、異なる表現型の発生にどの遺伝子が影響するかに関する情報を提供する。表現型という用語には、生物体の特徴の全ての種類、例えば、疾病状態等が含まれる。
【００８６】
遺伝子あたりのヌクレオチド配列数（すなわち、デジタルノザン）の統計的評価に関して、ｃＤＮＡ断片の重複的配列決定を確保することが重要である。このように、実験が実施される前に配列ライブラリの複雑性を確立し、十分な配列を得るのに必要な配列読み取り数を調整することは有用であり得る。例えば、典型的なｃＤＮＡ試料は、８，０００から１６，０００の異なる転写産物を含む。＋０／＋１のｃＤＮＡ−ＡＦＬＰが使用される場合、転写産物の総数の約８０％を標的とする４ヌクレオチドの配列を認識する２つの制限エンドヌクレアーゼが使用されると想定すると、複雑性の低下した試料は、約１，６００から３，２００の転写産物を含む。２０倍の重複的配列決定を用いると、これは、試料あたり必要とされる３２，０００から６４，０００の読み取りに相当する。これは、比較的低レベルで発現する遺伝子の転写産物レベルも決定し得るのに十分である。
【００８７】
配列ライブラリの複雑性を決定するための極めて適切な方法が、国際公開第０３／０１０３２８号に記載されており、参照により、本明細書に組み込まれる。
【００８８】
本発明は、
（ａ）第一のｃＤＮＡ試料に対して請求項２に記載の方法を実施することによって、前記第一のｃＤＮＡ試料中のヌクレオチド配列の頻度を決定する段階；
（ｂ）第二の及び／又はさらなるｃＤＮＡ試料に対して請求項２に記載の方法を実施することによって、第二の及び／又はさらなるｃＤＮＡ試料中の同一のヌクレオチド配列の頻度を決定する段階；及び
（ｃ）第一のｃＤＮＡ試料中のヌクレオチド配列の頻度を、第二の及び／又はさらなるｃＤＮＡ試料中の同一のヌクレオチド配列の頻度と比較して前記ヌクレオチド配列の相対的な転写レベルを得る段階
を含む、ｃＤＮＡ試料中のヌクレオチド配列の相対的な転写レベルを決定するための方法にも関する。
【００８９】
方法の段階（ａ）において、第一のｃＤＮＡ試料に対して請求項２に記載の方法を実施することによって、前記第一のｃＤＮＡ試料中でヌクレオチド配列の頻度を決定する。
【００９０】
方法の段階（ｂ）において、第二の及び／又はさらなるｃＤＮＡ試料に対して請求項２に記載の方法を実施することによって、前記第二の第二の及び／又はさらなるｃＤＮＡ試料中で同一のヌクレオチド配列の頻度を決定する。
【００９１】
段階（ｃ）において、前記第一のｃＤＮＡ試料中のヌクレオチド配列の頻度を、前記第二の及び／又はさらなるｃＤＮＡ試料中の同一のヌクレオチド配列の頻度と比較して前記ヌクレオチド配列の相対的な転写レベルを得る。
【００９２】
このような相対的な転写レベルに関する知識は、上述に論議されるように、特定の表現型に関して重要な転写産物を確立するために重要であり得る。
【００９３】
本発明は、
ａ）第一のｃＤＮＡ試料を準備する段階；
ｂ）第一のｃＤＮＡ試料に対して複雑性の低下を実施して第一のライブラリを得る段階；
ｃ）第一のライブラリにタグ付けして第一のタグ付きライブラリを得る段階；
ｄ）好ましくは各ｃＤＮＡ試料に対して異なるタグを使用して、第二の及び／又はさらなるｃＤＮＡ試料を用いて連続して又は同時に段階（ａ）及び（ｂ）を実施して第二の及び／又はさらなるタグ付きライブラリを得る段階；
ｅ）前記第一のタグ付きライブラリ並びに第二の及び／又はさらなるタグ付きライブラリを組み合わせて組み合わされたライブラリを得る段階；
ｆ）前記組み合わされたライブラリのヌクレオチド配列の少なくとも一部を配列決定によって決定する段階；
ｇ）第一のｃＤＮＡ試料並びに第二の及び／又はさらなるＤＮＡ試料中のヌクレオチド配列の頻度を決定する段階；並びに
ｈ）第一のｃＤＮＡ試料中のヌクレオチド配列の頻度を、第二の及び／又はさらなるｃＤＮＡ試料中のヌクレオチド配列の頻度と比較して前記ｃＤＮＡ試料中のヌクレオチド配列の相対的な転写レベルを得る段階
を含む、前記ｃＤＮＡ試料中のヌクレオチド配列の相対的な転写レベルを決定するための方法にも関する。
【００９４】
段階（ａ）において、第一のｃＤＮＡ試料を準備する。ｃＤＮＡ試料は、上述に論議されるように得られ得る。
【００９５】
段階（ｂ）において、第一のｃＤＮＡ試料に対して複雑性の低下を実施して第一のライブラリを得る。複雑性の低下は、何れかの技術によって実施され得るが、好ましくはＫｅｙｇｅｎｅのＡＦＬＰ（登録商標）技術によって実施される。
【００９６】
段階（ｃ）において、第一のタグ付きライブラリを得るため、第一のライブラリをタグ付けする。タグ付けは、複雑性の低下の段階（ｂ）と同時に起こり得る。このような同時タグ付けは、例えば、各試料のための固有の（ヌクレオチド）同定子を含むアダプターを使用して、ＡＦＬＰによって達成できる。
【００９７】
タグ付けは、組み合わされたライブラリを得るために、２つ又はそれ以上のｃＤＮＡ試料の２つ又はそれ以上の複雑性低下ライブラリが組み合わされる場合に、例えば、異なる植物の系から得られた異なる源の試料を識別することを目的とする。従って、好ましくは、第一のｃＤＮＡ試料及び第二又はさらなるｃＤＮＡ試料のタグ付きライブラリを調製するために、異なるタグが使用される。例えば、５個の核酸試料を使用する場合には、５個の異なるタグ付きライブラリを取得することが予定され、５個の異なるタグは、個々のもとの試料を表す。
【００９８】
タグは、核酸試料を識別するための、本分野で公知の何れのタグでもあり得るが、好ましくは短い同定子配列である。このような同定子配列は、例えば、複雑性の低下によって得られるライブラリの起源を示すために使用される変動する長さの固有の塩基配列であり得る。ライブラリをタグ付けするさらなる段階が必要ではないので、アダプター又はプライマー中にオリゴヌクレオチドタグを組み込むことは、非常に簡便である。このような同定子配列は、比較されるべき核酸試料の数に応じて長さが変動し得る。約試料の限定された数の源間を識別するには、４個の塩基の長さ（４４＝起こり得る２５６個の異なるタグ配列）が十分である（２５６個まで）が、タグ配列は、識別されるべき試料間で２つ以上の塩基が異なることが好ましい。必要に応じて、タグ配列の長さは、適宜調整することができる。
【００９９】
段階ｄ）において、好ましくは、各ｃＤＮＡ試料に対して異なるタグを使用し、第二の及び／又はさらなるｃＤＮＡ試料を用いて、連続して又は同時に段階（ａ）及び（ｂ）を実施して第二の及び／又はさらなるタグ付きライブラリを得る。ｃＤＮＡ試料は、例えば、異なる起源であり得、例えば、異なる植物の系であり得、これにより、このような植物の系のこのような転写産物の特性が比較され得る。あるいは、ｃＤＮＡ試料は、例えば、植物の発達中の転写産物の特性を比較するために、発達の異なる段階における単一の植物の系に由来し得る。単に有効性を調べるために、完全に無関係のｃＤＮＡ試料に対して、本発明に係る方法を実施することも可能である。
【０１００】
段階（ｅ）において、第一のタグ付きライブラリ並びに第二の及び／又はさらなるタグ付きライブラリを組み合わせて組み合わされたライブラリを得る。このような組み合わされたライブラリを同時配列決定に供して、非常に効果的なプロセスを提供し得る。
【０１０１】
段階（ｆ）において、組み合わされたライブラリのヌクレオチド配列の少なくとも一部を、配列決定、好ましくは高処理量配列決定によって、好ましくは上述のとおり決定する。
【０１０２】
段階（ｇ）において、第一のｃＤＮＡ試料並びに前記第二の及び／又はさらなるＤＮＡ試料中のヌクレオチド配列の頻度を決定する。タグによって、第一のライブラリのヌクレオチド配列を、第二の及び／又はさらなるライブラリのヌクレオチド配列と識別する。この場合、整列化は、アダプター／プライマーに関して整えられた配列データに対して、及び／又は同定子に関して整えられた配列データに対して行われるが、再構築された制限酵素認識配列を使用することによって、すなわちｃＤＮＡ由来の断片から配列データのみを使用することによって実施される。典型的には、得られた配列データは、断片の起源（すなわち、どの試料由来か）を同定するために使用され、アダプター及び／又は同定子配列に由来する配列は、データから除去され、この整えられたセットに対して整列化が実施される。
【０１０３】
段階（ｈ）において、第一のｃＤＮＡ試料中のヌクレオチド配列の頻度を、第二の及び／又はさらなるｃＤＮＡ試料中のヌクレオチド配列の頻度と比較してｃＤＮＡ試料中のヌクレオチド配列の相対的な転写レベルを得る。
【０１０４】
タグ付け戦略により、異なるｃＤＮＡ試料に関するヌクレオチド配列の転写レベルの決定は、同時に実施でき、同時実施が非常に有利である。本方法は、上で論述されているように、ある種の表現型形質に関与する転写産物を迅速に同定するのに非常に適している。
【０１０５】
好ましい実施形態において、第一のライブラリ及び第二の又はさらなるライブラリのタグ付けは、異なるタグを使用して実施される。上で論述されているように、、ｃＤＮＡ試料の各ライブラリが、それ自体のタグによって同定されることが好ましい。
【実施例】
【０１０６】
実施例
ノザンハイブリッド形成又はＤＮＡマイクロアレイ発現適用等のアプローチを使用して、高等植物中での遺伝子発現の時間的及び空間的制御の多数の例が蓄積されている。後者の技術は、遺伝子千個の発現を同時にモニターすることを可能とする。これらの分析方法とは異なり、遺伝子発現の特徴づけのデジタル分析は、高処理量配列技術を直接使用して、タグ付き転写産物を配列決定することによって実現され得る。試料中の特異的転写産物から得られた配列数は、この特定の配列の転写レベルを反映する。配列決定の深度を考慮しながら、複数の試料間でこれらの数を比較することによって、これらの試料間で転写レベルを正確に測定することが可能である。前記技術は、ある種の発現特性と関連した新たな未知の性質マーカーを発見するための強力なツールと思われる。
【０１０７】
本明細書において、本発明者は、２つのコショウ系からのｍＲＮＡ画分に由来するＡＦＬＰ技術を使用して複雑性を低下させる、ｃＤＮＡの高処理量配列決定を記載する。タグ付きｃＤＮＡ断片を直接配列決定することによって、発現特性が生じ得る。
【０１０８】
方法
全ＲＮＡ／ポリ（Ａ）＋ＲＮＡ単離
コショウ系ＰＳＰ１１及びＰＩ２０１２３４から、ＲＮｅａｓｙミニキット（カタログ番号７４１０４）を使用するＱＩＡＧＥＮ製ＲＮｅａｓｙ  Ｐｌａｎｔ  Ｍｉｎｉ  Ｐｒｏｔｏｃｏｌを使用して、葉の材料から全ＲＮＡを単離した。入力として、試料あたり１００ｍｇの葉の材料を使用した。
【０１０９】
本プロトコールの後、試料あたり２．５から３μｇの全ＲＮＡ収量を得た。その後、ＱＩＡＧＥＮ製Ｏｌｉｇｏｔｅｘ  ｍＲＮＡ  Ｍｉｎｉ  Ｋｉｔ（カタログ番号７００２２）を使用して、全ＲＮＡ試料１μｇからのポリ（Ａ）＋ＲＮＡ画分を単離した。１５０から２００ｎｇのポリ（Ａ）＋ＲＮＡ収量を得た。これらの試料の濃度は、５から１０ｎｇ／μＬであった。全ＲＮＡ及び全ポリ（Ａ）＋ＲＮＡの両者をアガロースゲルで分析し、ＲＮＡの性質をチェックした。
【０１１０】
ｃＤＮＡ合成
以下のプロトコールに従って、ｃＤＮＡを作製した。
【０１１１】
第一の鎖のｃＤＮＡ合成
以下を共に添加する。
１０μＬポリ（Ａ）＋ＲＮＡ（５０から１００ｎｇ）
５μＬオリゴ−ｄＴ２５（７０ｎｇ／μＬ）
その後、以下を添加する。
５μＬの５×第一の鎖の緩衝液（Ｓｕｐｅｒｓｃｒｉｐｔ  ＩＩ  ＲＴとともに供給）
２．５μＬの０．１Ｍ  ＤＴＴ
１μＬの１０ｍＭ  ｄＮＴＰ
０．５μＬのＳｕｐｅｒｓｃｒｉｐｔ  ＩＩ（２００Ｕ／μＬ）
１μＬのＭＱ水で２５μＬの最終容積にする。
４２℃で２時間定温放置する。
【０１１２】
第二の鎖のｃＤＮＡ合成
以下のものを、一緒に添加する。
２５μＬの第一の鎖の反応混合物
８μＬの１０×第二の鎖の緩衝液
１．５μＬの１０ｍＭ  ｄＮＴＰ
７．５単位のイー・コリ（Ｅ．ｃｏｌｉ）ＤＮＡリガーゼ
２５単位のイー・コリポリメラーゼ
０．８単位のＲＮａｓｅ−Ｈ（１Ｕ／μＬ）
ＭＱ水を添加して、８０μＬの最終容積にする。
１２℃で１時間定温放置する。
２２℃で１時間定温放置する。
【０１１３】
その後、ＱＩＡＧＥＮ製Ｑｉａｑｕｉｃｋ  ＰＣＲメンブレン精製キット（カタログ番号２８１０４）を使用して、ｃＤＮＡ試料を精製した。３０μＬの溶出緩衝液を使用して、溶出を実施した（５ｍＭトリス−ＨＣｌ、ｐＨ８．５）。
【０１１４】
タグ付きＡＦＬＰプライマーを使用するｃＤＮＡ−ＡＦＬＰテンプレート調製
Ｚａｂｅａｕ  ＆  Ｖｏｓ，１９９３：Ｓｅｌｅｃｔｉｖｅ  ｒｅｓｔｒｉｃｔｉｏｎ  ｆｒａｇｍｅｎｔ  ａｍｐｌｉｆｉｃａｔｉｏｎ；a  ｇｅｎｅｒａｌ  ｍｅｔｈｏｄ  ｆｏｒ  ＤＮＡ  ｆｉｎｇｅｒｐｒｉｎｔｉｎｇ．欧州特許出願公開（Ａ１）第０５３４８５８号；米国特許第６０４５９９４号及びＶｏｓら（Ｖｏｓ，Ｐ．，Ｈｏｇｅｒｓ，Ｒ．，Ｂｌｅｅｋｅｒ，Ｍ．，Ｒｅｉｊａｎｓ，Ｍ．，ｖａｎ  ｄｅ  Ｌｅｅ，Ｔ．，Ｈｏｒｎｅｓ，Ｍ．，Ｆｒｉｊｔｅｒｓ，Ａ．，Ｐｏｔ，Ｊ．，Ｐｅｌｅｍａｎ，Ｊ．，Ｋｕｉｐｅｒ，Ｍ．ｅｔ  ａｌ．（１９９５）ＡＦＬＰ：ａ  ｎｅｗ  ｔｅｃｈｎｉｑｕｅ  ｆｏｒ  ＤＮＡ  ｆｉｎｇｅｒｐｒｉｎｔｉｎｇ．Ｎｕｃｌ．Ａｃｉｄｓ  Ｒｅｓ．，２１，４４０７−４４１４）によって記載されているとおり、制限エンドヌクレアーゼの組み合わせＴａｑＩ／ＭｓｅＩを使用して、コショウ親系ＰＳＰ１１及びＰＩ−２０１２３４の作製されたｃＤＮＡのＡＦＬＰテンプレートを調製した。
【０１１５】
ｃＤＮＡの制限及び連結手法
２つの段階、すなわち、ＴａｑＩによる第一段階（最高定温放置温度）、ＭｓｅＩによるその後の段階（最低定温放置温度）において消化を実施した。
【０１１６】
ＴａｑＩ及びＭｓｅＩによるｃＤＮＡの制限を以下のとおり実施した。
【０１１７】
ＤＮＡ制限処理
以下のものを、一緒に添加する。
２５０ｎｇ  ｃＤＮＡ
１０単位のＴａｑＩ
８μＬの５×ＲＬ緩衝液  ５×ＲＬ緩衝液は、５０ｍＭトリス−ＨＡｃ、５０ｍＭ  ＭｇＡｃ、２５０ｍＭ  ＫＡｃ、２５ｍＭ  ＤＴＴ、２５０ｎｇ／ＢＳＡμＬ；ｐＨ７．５である。
ＭＱ水を添加して、４０μＬの最終容積にする。
６５℃で２時間定温放置する。
【０１１８】
ＴａｑＩによる制限の後、
以下を添加する。
１０単位のＭｓｅＩ
２μＬの５×ＲＬ緩衝液
ＭＱ水を添加して、５０μＬの最終容積にする。
３７℃で２時間定温放置する。
【０１１９】
アダプターの連結
消化混合物に以下の成分を添加した。
１μＬの１０ｍＭ  ＡＴＰ
１μＬのＴ４  ＤＮＡリガーゼ
１μＬのＴａｑＩアダプター（５０ｐｍｏｌ／μＬ）
【０１２０】
【化１】
2009520500.tif 000003
１μＬのＴａｑＩアダプター（５０ｐｍｏｌ／μＬ）
【０１２１】
【化２】
2009520500.tif 000004
２μＬの５×ＲＬ緩衝液。
ＭＱ水を添加して、６０μＬの最終容積にする。
３７℃で３時間定温放置する。
【０１２２】
ｃＤＮＡ−ＡＦＬＰ増幅
制限−連結の後、非選択的増幅段階において、この制限／連結反応性生物をテンプレートとして使用した。選択的増幅（＋１／＋１）のために、これらの非選択的ＡＦＬＰ生成物をその後テンプレートとして使用した。＋２／＋３の選択的増幅を実施することによって、前記＋１／＋１生成物に関して性質チェックを実施した。後者の増幅の生成物を４．５％配列ゲルでチェックした。
【０１２３】
非選択的ｃＤＮＡ−ＡＦＬＰ増幅を以下のとおり実施した。
５μＬの希釈されていない制限−連結混合物
１．５μＬのＴａｑＩ−プライマー（５０ｎｇ／μＬ）
【化３】
2009520500.tif 000005
１．５μＬのＭｓｅＩ−プライマー（５０ｎｇ／μＬ）
【化４】
2009520500.tif 000006
２μＬの５ｍＭ  ｄＮＴＰ
１単位のＴａｑ．ポリメラーゼ
５μＬの１０×ＰＣＲ緩衝液
ＭＱ水を添加して、５０μＬの最終容積にする。
【０１２４】
次の条件、すなわち３０周期（９４℃で３０秒、５６℃で６０秒及び７２℃で１２０秒）を使用して、金又は銀のブロックとともにＰＥ９７００を使用して、ＰＣＲ増幅を実施した。
【０１２５】
タグ配列を使用する選択的ｃＤＮＡ−ＡＦＬＰ増幅を以下のとおり実施した。
コショウ系ＰＳＰ１１に由来する非選択的ｃＤＮＡ−ＡＦＬＰ生成物用
５μＬの６００×希釈した非選択的生成物
１．５μＬのＴｒ０１ＡＣＡＣプライマー（＋Ａ）＊（５０ｎｇ／μｇ）
【０１２６】
【化５】
2009520500.tif 000007
１．５μＬのＭ０２ＡＣＡＣプライマー（＋Ｃ）＊（５０ｎｇ／μｇ）
【０１２７】
【化６】
2009520500.tif 000008
２μＬの５ｍＭ  ｄＮＴＰ
１．５単位のＡｍｐｌｉＴａｑ−Ｇｏｌｄポリメラーゼ
５μＬの１０×ＰＣＲ緩衝液
ＭＱ水を添加して、５０μＬの最終容積にする。
【０１２８】
コショウ系ＰＩ２０１２３４に由来する非選択的ｃＤＮＡ−ＡＦＬＰ０／０生成物用
５μＬの６００×希釈した非選択的生成物
１．５μＬのＴｒ０１ＡＧＣＴプライマー（＋Ａ）＊（５０ｎｇ／μｇ）
【０１２９】
【化７】
2009520500.tif 000009
１．５μＬのＭ０２ＡＧＣＴプライマー（＋Ｃ）＊（５０ｎｇ／μｇ）
【０１３０】
【化８】
2009520500.tif 000010
２μＬの５ｍＭ  ｄＮＴＰ
１．５単位のＡｍｐｌｉＴａｑ−Ｇｏｌｄポリメラーゼ
５μＬの１０×ＰＣＲ緩衝液
ＭＱ水を添加して、５０μＬの最終容積にする。
【０１３１】
次の条件、すなわち９４℃で１２秒（高温開始）、９４℃で３０秒、６５℃で３０秒、７２℃で６０秒を１周期；２３周期−より低温の徐冷温度であって１２周期中各周期０．７℃−１３周期のタッチダウン相−９４℃で３０秒、５６℃で３０秒、７２℃で６０秒を使用して、金ブロックとともにＰＥ９７００を使用して、ＰＣＲ増幅を実施した。１００塩基対ラダーを使用する１％アガロースゲル上で、生成された＋１/＋１産物の品質をチェックし、断片長の分布をチェックした（図１参照）。
【０１３２】
選択的プライマーは、その５’末端で４ｂｐタグ（上述の下線部）を含有し、配列決定処理の終了時に個々のコショウ系由来の増幅産物を識別する。本発明に従ったタグ付きｃＤＮＡ−ＡＦＬＰ  ＰＣＲ産物を生成する原理を図２に示す。
【０１３３】
配列ライブラリ調製及び高処理量配列決定
Ｍａｒｇｕｌｉｅｓら（Ｍａｒｇｕｌｉｅｓ  ｅｔ  ａｌ．，Ｎａｔｕｒｅ  ４３７，ｐｐ．３７６−３８０及びオンライン補遺）によって記載されている４５４  Ｌｉｆｅ  Ｓｃｉｅｎｃｅｓ／Ｒｏｃｈｅ  ＧＳ２０配列決定技術を使用する高処理量配列決定に、両コショウ系由来のタグ付きｃＤＮＡ  ＡＦＬＰ生成物を供した。
タグ付きｃＤＮＡ−ＡＦＬＰ  ＰＣＲ生成物をまず生成し、修飾されたアダプター（
【０１３４】
【化９】
2009520500.tif 000011
及び
【０１３５】
【化１０】
2009520500.tif 000012
）に連結し、乳濁液−ＰＣＲ増幅並びにＭａｒｇｕｌｉｅｓ及び共同研究者によって記載されているその後の断片配列決定を容易にした。
【０１３６】
乳濁液ＰＣＲプライマー、配列プライマー及び配列実行条件は全て、Ｍａｒｇｕｌｉｅｓ及び共同研究者によって記載されるとおりであった。配列ライブラリ調製手法を図３に示す。Ｋｅｙｇｅｎｅ  ＮＶ，Ｗａｇｅｎｉｎｇｅｎ，Ｔｈｅ  Ｎｅｔｈｅｒｌａｎｄｓの研究室で、高処理量ＧＳ２０配列実行を実施した。
【０１３７】
ＧＳ２０配列実行データ処理
バイオインフォマティクスパイプライン（Ｋｅｙｇｅｎｅ  Ｎ．Ｖ．）を使用して、ＧＳ２０配列実行の半分（すなわち、ＧＳ２０ＰｉｃｏＴｉｔｅｒＰｌａｔｅ上で利用可能な２チャネルのうちの１チャネル）から得られる配列データを処理した。具体的には、塩基が呼び出された（ｂａｓｅｃａｌｌｅｄ）生配列の読み取りをＦＡＳＴＡフォーマット中で変換し、ＢＬＡＳＴアルゴリズムを使用して、タグ付きＡＦＬＰアダプター配列の存在に関して検査した。既知のタグ付きＡＦＬＰプライマー配列に対して信頼度の高い合致が得られたら、配列を整え、制限エンドヌクレアーゼ部位を復元し、適切なタグを割り当てた。その後、全体的な配列相同性に基づいたメガＢＬＡＳＴ手法を使用して、３３塩基よりも大きな整えられた配列を全てクラスター化した。次に、ＣＡＰ３多重整列化アルゴリズムを使用して、クラスターあたり１つ又はそれ以上の連結断片中にクラスターを構築した。
【０１３８】
１３個の配列読み取りの出力例
クラスター３８７
試料２のＩＤタグ（ＡＧＴＣ）を太字で示している。試料１のＩＤタグ（ＡＣＡＣ）を下線で示している。図４を参照されたい。
【０１３９】
配列実行の全体的な統計を表１に示す。
【０１４０】
【表１】
2009520500.tif 000013
【０１４１】
解釈：
段階１）「試料配列決定深度標準化因子」は、２．４５であり、試料２から得られた読み取り全体を試料１に由来する読み取り総数で除したものとして定義される（１２３８２２／５０５９９＝２．４５）。連結断片あたりの試料２に由来する読み取り数を２．４５で除して、転写レベルを試料１の転写レベルと比較した。
【０１４２】
段階２）内部標準物質として供する「ハウスキーピング」遺伝子の「発現」を決定することによって、第二の「ハウスキーピング遺伝子の標準化」段階を実施した。このため、リコペルシコン・エスキュレンタム（Ｌｙｃｏｐｅｒｓｉｃｏｎ  ｅｓｃｕｌｅｎｔｕｍ）アルギニンデカルボキシラーゼ遺伝子を選択した。ＣＡＰ３多重整列化を使用して得られた連結断片の配列に対して、リコペルシコン・エスキュレンタムアルギニンデカルボキシラーゼの配列を「ＢＬＡＳＴ処理」し、コショウアルギニンデカルボキシラーゼ遺伝子の転写産物が、試料１及び２においてどの程度頻繁に観察されたかを決定した。その後、「試料配列決定深度標準化因子」を最初に適用した後、これらの転写産物が試料１及び２において観察された比を算出した（段階１）。本実施例において、前記比（＝ハウスキーピング遺伝子標準化因子）は、試料１／試料２に関して１７／１４＝１．２であった（表１）。
【０１４３】
連結断片プールに対するハウスキーピング遺伝子（リコペルシコン・エスキュレンタムアルギニンデカルボキシラーゼ）のＢＬＡＳＴ検索例
参考文献：Ａｌｔｓｃｈｕｌ，Ｓｔｅｐｈｅｎ  Ｆ．，Ｔｈｏｍａｓ  Ｌ．Ｍａｄｄｅｎ，Ａｌｅｊａｎｄｒｏ  Ａ．Ｓｃｈａｆｆｅｒ，Ｊｉｎｇｈｕｉ  Ｚｈａｎｇ，Ｚｈｅｎｇ  Ｚｈａｎｇ，Ｗｅｂｂ  Ｍｉｌｌｅｒ，ａｎｄ  Ｄａｖｉｄ  Ｊ．Ｌｉｐｍａｎ（１９９７），“Ｇａｐｐｅｄ  ＢＬＡＳＴ  ａｎｄ  ＰＳＩ−ＢＬＡＳＴ：ａ  ｎｅｗ  ｇｅｎｅｒａｔｉｏｎ  ｏｆ  ｐｒｏｔｅｉｎ  ｄａｔａｂａｓｅ  ｓｅａｒｃｈ  ｐｒｏｇｒａｍｓ”，Ｎｕｃｌｅｉｃ  Ａｃｉｄｓ  Ｒｅｓ．２５：３３８９−３４０２。
【０１４４】
クエリ＝ｇｉ｜２９３５４９｜ｇｂ｜Ｋ１６５８２｜ＴＯＭＡＲＧＤＥＣＡ  リコペルシコン・エスキュレンタムアルギニンデカルボラーゼｍＲＮＡ、完全コード（２０６０文字）
データベース：ｔａｇｇｅｄＲｅａｄｓ．ｆｎａ
１７４，４２１配列；１５，４０８，１９２全文字。結果を図５に示す。
【０１４５】
【表２】
2009520500.tif 000014
【０１４６】
段階３）実際の発現特徴づけに関し、１０超の読み取りを含有する連結断片のみを考慮した。連結断片あたりの１０個の読み取りの最小レベルを選択して、不十分な配列決定深度による不正確な転写産物の特徴づけの結果を回避した。表２は、上述に概略された３段階手法後のＰＳＰ１１（試料１）対ＰＩ２０１２３４（試料２）において差動的に発現した２つの転写産物の相対的なｍＲＮＡ発現レベルを示す。具体的には、クラスター２２１５は、試料１において上方制御された転写産物を表し、クラスター８４７は、試料１において下方制御された転写産物を表し；これらの転写産物の相対的な転写レベルの算出を表３に示す。最後に、表４は、上述の原理の基づいた全データセットにおける差動的に転写された遺伝子数の概略を含有する。
【０１４７】
上方制御試料１の例−生データ。クラスター２２１５。試料２のＩＤタグ（ＡＧＴＣ）を太字で示している。試料１のＩＤタグ（ＡＣＡＣ）を図６において下線で示している。
【０１４８】
下方制御試料１の例−生データ。クラスター８４７。試料２のＩＤタグ（ＡＧＴＣ）を太字で示している。試料１のＩＤタグ（ＡＣＡＣ）を図７において下線で示している。
【０１４９】
【表３】
2009520500.tif 000015
【０１５０】
【表４】
2009520500.tif 000016
【図面の簡単な説明】
【０１５１】
【図１】図１：タグ付き（Ａ／Ｃ）ｃＤＮＡ−ＡＦＬＰ生成物は、コショウ系ＰＳＰ１１及びＰＩ２０１２３４を形成する。両系由来の２つの試料を１％アガロースゲル上で二重に負荷する。  Ｍ：１００ｂｐマーカー  １：ｃＤＮＡ−ＡＦＬＰ  ＰＳＰ１１試料１  ２：ｃＤＮＡ−ＡＦＬＰ  ＰＳＰ１１試料１  ３：ｃＤＮＡ−ＡＦＬＰ  ＰＳＰ１１試料２  ４：ｃＤＮＡ−ＡＦＬＰ  ＰＳＰ１１試料２  ５：ｃＤＮＡ−ＡＦＬＰ  ＰＩ２０１２３４−試料１  ６：ｃＤＮＡ−ＡＦＬＰ  ＰＩ２０１２３４−試料１  ７：ｃＤＮＡ−ＡＦＬＰ  ＰＩ２０１２３４−試料２  ８：ｃＤＮＡ−ＡＦＬＰ  ＰＩ２０１２３４−試料２
【図２】図２：４ｂｐの５プライムタグ配列を含有するＡＦＬＰプライマーによる増幅後のコショウＡＦＬＰ＋１／＋１増副産物の模式的呈示
【図３】図３：配列ライブラリの調製の作業の流れ
【図４】図４：１３個の配列読み取りの出力例
【図５】図５：ＢＬＡＳＴの結果
【図６−１】図６−１：上方制御の生データの呈示
【図６−２】図６−２：上方制御の生データの呈示
【図６−３】図６−３：上方制御の生データの呈示
【図６−４】図６−４：上方制御の生データの呈示
【図７−１】図７−１：上方制御の生データの呈示
【図７−２】図７−２：上方制御の生データの呈示
【図７−３】図７−３：上方制御の生データの呈示
【図７−４】図７−４：上方制御の生データの呈示
【図７−５】図７−５：上方制御の生データの呈示
【図７−６】図７−６：上方制御の生データの呈示
【図１】
2009520500.tif 000017
【図２】
2009520500.tif 000018
【図３】
2009520500.tif 000019
【図４】
2009520500.tif 000020
【図５】
2009520500.tif 000021
【図６−１】
2009520500.tif 000022
【図６−２】
2009520500.tif 000023
【図６−３】
2009520500.tif 000024
【図６−４】
2009520500.tif 000025
【図７−１】
2009520500.tif 000026
【図７−２】
2009520500.tif 000027
【図７−３】
2009520500.tif 000028
【図７−４】
2009520500.tif 000029
【図７−５】
2009520500.tif 000030
【図７−６】
2009520500.tif 000031
【配列表】
2009520500000001.app
【国際調査報告】
2009520500.tif 000032
2009520500.tif 000033
2009520500.tif 000034
2009520500.tif 000035