候補タンパク質をコードするcDNAライブラリーを逆転写してmRNAライブラリーを調製する転写工程、前記転写工程で調製されたmRNAライブラリーからmRNAタンパク質対応付け分子ライブラリーを調製する対応付け工程、前記対応付け工程で調製されたmRNAタンパク質対応付け分子ライブラリーから、標的物質と相互作用するmRNAタンパク質対応付け分子を選択する選択工程、及び、前記選択工程で選択されたmRNAタンパク質対応付け分子のmRNA部分に基づく核酸増幅により、候補タンパク質をコードするcDNAライブラリーを調製する増幅工程、増幅工程で調製されたcDNAライブラリーを用いて工程~を繰り返すことを含む、標的物質と相互作用するタンパク質の検出方法であって、最初のcDNAライブラリーの調製及び繰返しの増幅工程におけるcDNAライブラリーの調製のうち複数回のcDNAライブリーの調製において、その調製回に特異的な配列を有するプライマーを用いて前記cDNAライブラリーを調製し、前記調製回に特異的な配列を有するプライマーを用いて調製されたcDNAライブラリーを混合し、混合されたcDNAライブラリーについて配列決定を行い、決定された配列を、前記調製回に特異的な配列に基づいて、調製回毎における、同一候補タンパク質をコードする分子の数を測定し、前記調製回に応じて有意に増加している分子にコードされる候補タンパク質を、前記標的物質と相互作用するタンパク質として検出することを含む前記方法。工程におけるcDNAライブラリーの混合比において、早い調製回の量を次の調製回の量以上とする請求項1に記載の方法。前記調製回に特異的な配列の長さが4~10塩基長である請求項1又は2に記載の方法。前記プライマーが前記調製回に特異的な配列を5末端に有する請求項1~3のいずれか一項に記載の方法。工程で混合するcDNAライブラリーの数が3以上である請求項1~4のいずれか一項に記載の方法。工程で混合するcDNAライブラリーが最初のcDNAライブラリーを含む請求項1~5のいずれか一項に記載の方法。工程において分子の増加を統計学的方法により評価することを含む請求項1~6のいずれか一項に記載の方法前記標的物質がタンパク質である請求項1~7のいずれか一項に記載の方法。IVV法と次世代シーケンサーと組み合わせることで、10万から1億リードの配列解析が可能となり、これまでの偽陰性問題ライブラリーの中の配列解析出来ていなかったデータの問題が著しく改善することが期待できる。しかしながら、大量の配列解析は、大量の偽陽性データも生産することになり、偽陰性問題のみだけでなく、偽陽性問題も同時に解決することが大きな課題となる。本発明は、偽陰性問題のみだけでなく、偽陽性問題も同時に解決することが可能な相互作用の検出方法を提供することを目的とする。本発明者らは、ライブラリー中における目的タンパク質の絶対量による影響を受けずに、ラウンド毎の濃縮を検出できれば、偽陽性問題を解決し得るという仮説をたて、このような検出が、IVV法の特性を利用して特定の様式で調製した試料の配列解析を行うことにより可能になることを、次世代シーケンサーにより解析することにより見出し、本発明を完成した。本発明は、以下の方法を提供する。候補タンパク質をコードするcDNAライブラリーを逆転写してmRNAライブラリーを調製する転写工程、前記転写工程で調製されたmRNAライブラリーからmRNAタンパク質対応付け分子ライブラリーを調製する対応付け工程、前記対応付け工程で調製されたmRNAタンパク質対応付け分子ライブラリーから、標的物質と相互作用するmRNAタンパク質対応付け分子を選択する選択工程、及び、前記選択工程で選択されたmRNAタンパク質対応付け分子のmRNA部分に基づく核酸増幅により、候補タンパク質をコードするcDNAライブラリーを調製する増幅工程、増幅工程で調製されたcDNAライブラリーを用いて工程~を繰り返すことを含む、標的物質と相互作用するタンパク質の検出方法であって、最初のcDNAライブラリーの調製及び繰返しの増幅工程におけるcDNAライブラリーの調製のうち複数回のcDNAライブリーの調製において、その調製回に特異的な配列を有するプライマーを用いて前記cDNAライブラリーを調製し、前記調製回に特異的な配列を有するプライマーを用いて調製されたcDNAライブラリーを混合し、混合されたcDNAライブラリーについて配列決定を行い、決定された配列を、前記調製回に特異的な配列に基づいて、調製回毎における、同一候補タンパク質をコードする分子の数を測定し、前記調製回に応じて有意に増加している分子にコードされる候補タンパク質を、前記標的物質と相互作用するタンパク質として検出することを含む前記方法。本発明方法においては、cDNAのライブラリーを混合する量比において、早い調製回の量を次の調製回の量以上とすることが好ましい。調製回に特異的な配列の長さは通常には4~10塩基長である。前記プライマーは、通常には、調製回に特異的な配列を5末端に有する。工程で混合するcDNAライブラリーの数は通常には3以上である。工程で混合するcDNAライブラリーが最初のcDNAライブラリーを含んでもよい。工程において分子の増加を統計学的方法により評価することが好ましい。標的物質は、例えば、タンパク質である。タンパク質の場合は、タンパク質間相互作用を検出できる。本発明により、偽陰性及び偽陽性の問題の軽減された検出方法が提供される。本発明の方法によれば、ライブラリーに含まれる絶対量の小さいタンパク質をコードする分子も検出できる。また、従来の方法では、偽陽性を排除するために、大量のリアルタイムPCRなどでの検証実験を行っていたが、それらが必要なくなる。データの信頼性が向上するだけでなく、実験時間やコストの大幅な削減になる。本発明の方法は、IVV法によるスクリーニングと、配列の解析とを組み合わせたものである図1。第1の部分は、特定のプライマーを用いる増幅によりcDNAライブラリーを調製する他は、従来のIVV法による方法と同様でよい。第2の部分は、特定のプライマーに含まれる配列を利用して配列を解析するものである。先ず、第1のスクリーニング部分について説明する。この部分は、基本的にIVVを用いたPPI検出のための工程と同様でよい。この工程は更にいくつかの部分に分けることができる。すなわち、cDNAライブラリーからのmRNAの転写、通常にはPEGPuroスペーサーのmRNAへのライゲーション及び翻訳からなるIVV形成対応付け、選択、そして、次のラウンドで用いるcDNAライブラリーを調製する増幅である。これらのinvitro工程を繰り返す通常には3~5ラウンド程度ことによって、ベイトに結合するプレイ分子は徐々にライブラリー中で濃縮されていき、RTPCRプロダクトのシーケンシングによって同定することができるようになる。以下、工程毎に説明する。転写工程本工程は、候補タンパク質をコードするcDNAライブラリーを逆転写してmRNAライブラリーを調製する転写工程である。候補タンパク質をコードするcDNAライブラリーは、特許文献3に記載のcDNAライブラリーの製造方法により調製できる。cDNAライブラリーの調製に使用されるRNA又はmRNAライブラリーは、原核、真核生物、ウイルスなどあらゆる種のいかなる組織から抽出したRNA又はmRNAライブラリーでも構わない。また、解読したゲノムやcDNAライブラリーを転写したRNAライブラリーやそれを再現した人工のRNAライブラリー、あるいは自然には存在しない配列を含む人工のcDNAライブラリーを転写したRNAライブラリーでも構わない。cDNAライブラリーを逆転写してmRNAライブラリーを調製する転写工程は、特許文献3に記載のRNAライブラリーの製造方法により行うことができる。対応付け工程本工程は、転写工程で調製されたmRNAライブラリーからmRNAタンパク質対応付け分子ライブラリーを調製する対応付け工程である。この工程は、特許文献2又は特許文献3に記載の方法により行なうことができる。以下、具体例を説明する。対応付け分子本明細書において、対応付け分子とは、表現型と遺伝子型と対応付ける分子を意味する。対応付け分子は、遺伝子型を反映する塩基配列を有する核酸を含む遺伝子型分子と、表現型の発現に関与するタンパク質を含む表現型分子とが結合してなる。遺伝子型分子は、遺伝子型を反映する塩基配列を、その塩基配列が翻訳され得るような形態で有するコード分子と、スペーサー部とが結合してなる。このような対応付け分子は、例えば、特許文献2に記載されているので特許文献2を参照して説明する。対応付け分子における、表現型分子に由来する部分、スペーサー分子に由来する部分、及び、コード分子に由来する部分をそれぞれ、デコード部、スペーサー部及びコード部と呼ぶ。また、遺伝子型分子における、スペーサー分子に由来する部分、及び、コード分子に由来する部分をそれぞれ、スペーサー部及びコード部と呼ぶ。特許文献2図8に、対応付け分子、スペーサー分子及びコード分子の一例の大まかな構成が示されている。この対応付け分子は、ピューロマイシンを含むスペーサーと表現型のコードを反映する塩基配列からなる。この対応付け分子は、コード分子に何らかの方法によってピューロマイシンを含むスペーサー部を結合して遺伝子型分子とし、無細胞翻訳系において、リボソーム上で表現型分子と連結した構成をもつ。スペーサー分子は、ポリエチレングリコールを主成分としたPEG領域、少なくともピューロマイシンあるいはピューロマイシンと1残基以上のDNAあるいはまたはRNAからなるCCA領域、少なくとも1残基以上のDNAあるいはまたはRNAを含むドナー領域、さらに、少なくとも1残基のDNAあるいはまたはRNAの塩基に機能修飾を施した機能付与ユニットからなる。コード分子は、デコード部の一部の配列からなるDNAあるいはまたはRNAのポリA配列を含む3末端領域、および、DNAあるいはまたはRNAからなる転写プロモーターおよび翻訳エンハンサーを含んだ5UTR、さらに、主として表現型分子の配列からなるORF領域から構成される。以下、この例を参照して説明するが、本発明はこれに限定されるものではない。スペーサー分子スペーサー分子は、核酸の3末端に結合できるドナー領域と、ドナー領域に結合した、ポリエチレングリコールを主成分としたPEG領域と、PEG領域に結合した、ペプチド転移反応によってペプチドと結合し得る基を含むペプチドアクセプター領域とを含む。核酸の3末端に結合できるドナー領域は、通常、1以上のヌクレオチドからなる。ヌクレオチドの数は、通常には1~15、好ましくは1~2である。ヌクレオチドはリボヌクレオチドでもデオキシリボヌクレオチドでもよい。ドナー領域の5末端の配列は、ライゲーション効率を左右する。コード部とスペーサー部をライゲーションさせるためには、少なくとも1残基以上を含むことが必要であり、ポリA配列をもつアクセプターに対しては、少なくとも1残基のdCあるいは2残基のdCdCが好ましい。塩基の種類としては、CU又はTGAの順で好ましい。PEG領域はポリエチレングリコールを主成分とするものである。ここで、主成分とするとは、PEG領域に含まれるヌクレオチドの数の合計が20bp以下、又は、ポリエチレングリコールの平均分子量が400以上であることを意味する。好ましくは、ヌクレオチドの合計の数が10bp以下、又は、ポリエチレングリコールの平均分子量が1000以上であることを意味する。PEG領域のポリエチレングリコールの平均分子量は、通常には、400~30000、好ましくは1000~10000、より好ましくは2000~8000である。ここで、ポリエチレングリコールの分子量が約400より低いと、このスペーサー分子に由来するスペーサー部を含む遺伝子型分子を対応付け翻訳したときに、対応付け翻訳の後処理が必要となることがあるがLiuRBarrickESzostakJWRobertsRWMethodsinEnzymologyvol318268293、分子量1000以上、より好ましくは2000以上のPEGを用いると、対応付け翻訳のみで高効率の対応付けができるため、翻訳の後処理が必要なくなる。また、ポリエチレングリコールの分子量が増えると、遺伝子型分子の安定性が増す傾向があり、特に分子量1000以上で良好であり、分子量400以下ではDNAスペーサーと性質がそれほどかわらず不安定となることがある。ペプチドアクセプター領域は、ペプチドのC末端に結合できるものであれば特に限定されないが、例えば、ピューロマイシン、3Nアミノアシルピューロマイシンアミノヌクレオシド3NAminoacylpuromycinaminonucleosidePANSアミノ酸、例えばアミノ酸部がグリシンのPANSGly、バリンのPANSVal、アラニンのPANSAla、その他、全アミノ酸に対応するPANS全アミノ酸が利用できる。また、化学結合として3アミノアデノシンのアミノ基とアミノ酸のカルボキシル基が脱水縮合した結果形成されたアミド結合でつながった3Nアミノアシルアデノシンアミノヌクレオシド3AminoacyladenosineaminonucleosideAANSアミノ酸、例えばアミノ酸部がグリシンのAANSGly、バリンのAANSVal、アラニンのAANSAla、その他、全アミノ酸に対応するAANS全アミノ酸が利用できる。また、ヌクレオシドあるいはヌクレオシドとアミノ酸のエステル結合したものなども利用できる。その他、ヌクレオシドあるいはヌクレオシドに類似した化学構造骨格を有する物質と、アミノ酸あるいはアミノ酸に類似した化学構造骨格を有する物質を化学的に結合可能な結合様式のものなら全て利用することができる。ペプチドアクセプター領域は、好ましくは、ピューロマイシンもしくはその誘導体、又は、ピューロマイシンもしくはその誘導体と1残基もしくは2残基のデオキシリボヌクレオチドもしくはリボヌクレオチドからなることが好ましい。ここで、誘導体とはタンパク質翻訳系においてペプチドのC末端に結合できる誘導体を意味する。ピューロマイシン誘導体は、ピューロマイシン構造を完全に有しているものに限られず、ピューロマイシン構造の一部が欠落しているものも包含する。ピューロマイシン誘導体の具体例としては、PANSアミノ酸、AANSアミノ酸などが挙げられる。ペプチドアクセプター領域は、ピューロマイシンのみの構成でもかまわないが、5側に1残基以上のDNAあるいはまたはRNAからなる塩基配列を持つことが好ましい。配列としては、dCピューロマイシンrCピューロマイシンなど、より好ましくはdCdCピューロマイシンrCrCピューロマイシンrCdCピューロマイシンdCrCピューロマイシンなどの配列で、アミノアシルtRNAの3末端を模倣したCCA配列PhilippsGRNature223374377が適当である。塩基の種類としては、CU又はTGAの順で好ましい。コード分子コード分子は、転写プロモーターおよび翻訳エンハンサーを含む5非翻訳領域と、5非翻訳領域の3側に結合した、タンパク質をコードするORF領域と、ORF領域の3側に結合した、ポリA配列及び、必要によりその5側に翻訳増強配列例えば制限酵素XhoIが認識する配列を含む3末端領域を含む核酸である。コード分子は、前記転写工程で調製されたmRNAライブラリーを構成する。3末端領域は、好ましくはSNNS配列例えばXhoI配列とその下流にポリA配列を含む。スペーサー分子とコード分子とのライゲーション効率に影響を与える要素としては3末端領域のポリA配列が重要であり、ポリA配列は、少なくとも2残基以上のdAあるいはまたはrAの混合あるいは単一のポリA連続鎖であり、好ましくは、3残基以上、より好ましくは6以上、さらに好ましくは8残基以上のポリA連続鎖である。コード分子の翻訳効率に影響する要素としては、転写プロモーターと翻訳エンハンサーからなる5UTR、および、ポリA配列を含む3末端領域の組み合わせがある。3末端領域のポリA配列の効果は通常には10残基以下で発揮される。5UTRの転写プロモーターはT7T3あるいはSP6などが利用でき、特に制限はない。好ましくはSP6であり、特に、翻訳のエンハンサー配列としてオメガ配列やオメガ配列の一部を含む配列Ω様配列を利用する場合はSP6を用いることが特に好ましい。翻訳エンハンサーは好ましくはオメガ配列の一部であり、オメガ配列の一部としては、TMVのオメガ配列GallieDRWalbotVNucleicAcidsResvol2046314638の一部を含んだものが好ましい。また、翻訳効率に関し、3末端領域においては、XhoI配列とポリA配列の組み合わせが重要となる。また、ORF領域の下流部分、すなわちXhoI配列の上流に親和性タグがついたものとポリA配列の組み合わせも重要となる。親和性タグ配列としては、抗原抗体反応など、タンパク質を検出できるいかなる手段を用いるための配列であればよく、制限はない。好ましくは、抗原抗体反応によるアフィニティー分離分析用タグであるFlagtag配列である。ポリA配列効果としては、Flagtag等の親和性タグにXhoI配列がついたものとそこへさらにポリA配列がついたものの翻訳効率が上昇する。上記の翻訳効率に関し効果のある構成は、対応付け効率にも有効である。ORF領域については、遺伝子配列、エキソン配列、イントロン配列、ランダム配列、あるいは、いかなる自然界の配列、人為的配列が可能であり、配列の制限はない。また、コード分子の5UTRをSP6O29とし、3末端領域を、たとえば、FlagXhoIAnとすることで、各長さは、5UTRで約60bp、3末端領域で約40bpであり、PCRのプライマーにアダプター領域として組み込める長さである。このため、あらゆるベクターやプラスミドやcDNAライブラリーからPCRによって、5UTRと3末端領域をもったコード分子を簡単に作成できる。コード分子において、翻訳はORF領域を超えてされてもよい。すなわち、ORF領域の末端に終止コドンがなくてもよい。特許文献2図10に、コード分子の一例の詳細な構成が示されている。コード分子は、3末端領域と、転写プロモーターおよび翻訳エンハンサーを含む5UTRと、デコード部の配列情報からなる、すなわち表現型タンパク質をコードするORF領域とからなる。ここでは、3末端領域として、親和性タグ配列、XhoI配列、ポリA配列を含み、Flagtag配列を用いている。5UTRとして、転写プロモーターのSP6、翻訳エンハンサーのオメガ配列の一部であるO29を含む配列を用いている。遺伝子型分子およびその製造方法遺伝子型分子は、転写プロモーターおよび翻訳エンハンサーを含む5非翻訳領域と、5非翻訳領域の3側に結合した、タンパク質をコードするORF領域と、ORF領域の3側に結合した、ポリA配列を含む3末端領域を含む核酸であるコード分子の3末端と、スペーサー分子のドナー領域とが結合してなる。遺伝子型分子を構成するコード分子は、上記のコード分子においてXhoI配列が必須ではない他は、コード分子について説明したとおりである。しかしながら、XhoI配列を有することが好ましい。遺伝子型分子は、上記コード分子の3末端と、スペーサー分子のドナー領域を、通常のリガーゼ反応により結合させることにより製造できる。反応条件としては、通常、4~25°Cで4~48時間の条件が挙げられ、PEG領域を含むスペーサー分子のPEG領域内のポリエチレングリコールと同じ分子量のポリエチレングリコールを反応系に添加する場合には、15°Cで05~4時間に短縮することも可能である。スペーサー分子とコード分子の組み合わせはライゲーション効率に重要な効果をもたらす。アクセプターにあたるコード部の3末端領域において、少なくとも2残基以上、好ましくは3残基以上、さらに好ましくは6~8残基以上のDNAあるいはまたはRNAのポリA配列があること、さらに、5UTRの翻訳エンハンサーとしては、オメガ配列の部分配列が好ましく、スペーサー部のドナー領域としては、少なくとも1残基のdCあるいは2残基のdCdCが好ましい。このことによって、RNAリガーゼを用いることでDNAリガーゼのもつ問題点を回避し、かつ効率を60~80に保つことができる。転写プロモーターおよび翻訳エンハンサーを含む5非翻訳領域と、5非翻訳領域の3側に結合した、タンパク質をコードするORF領域と、ORF領域の3側に結合した、ポリA配列を含む3末端領域を含むコード分子の3末端と、1~4のいずれか1項に記載のスペーサー分子のドナー領域であってRNAからなるものとを、スペーサー分子内のPEG領域を構成するポリエチレングリコールと同じ分子量をもつ遊離のポリエチレングリコールの存在下で、RNAリガーゼにより結合させることが好ましい。ライゲーション反応時に、PEG領域を含むスペーサー部のPEG領域と同じ分子量のポリエチレングリコールを添加することによって、スペーサー部のポリエチレングリコールの分子量によらずライゲーション効率が80~90以上に向上し、反応後の分離工程も省略することができる。対応付け分子及びその製造方法mRNAタンパク質対応付け分子は、上記の遺伝子型分子を、ペプチド転移反応で、遺伝子型分子内のORF領域によりコードされたタンパク質である表現型分子と連結してなるものである。対応付け分子は、遺伝子型分子を無細胞翻訳系で翻訳することにより、ペプチド転移反応で、遺伝子型分子内のORF領域によりコードされたタンパク質である表現型分子と連結することを含む。無細胞翻訳系は、好ましくは、小麦胚芽又はウサギ網状赤血球のものである。翻訳の条件は通常に採用される条件でよい。例えば、25~37°Cで15~240分の条件が挙げられる。無細胞翻訳系については、これまで大腸菌、ウサギ網状赤血球、小麦胚芽の系で対応付け分子の形成が検討され、ウサギ網状赤血球の系でのみ対応付け分子が確認されていたがNemotoNMiyamotoSatoEYanagawaHFEBSLett414405RobertsRWSzostakJWProcNatlAcadSciUSA9412297、この態様によれば、PEG領域を含むスペーサー部をもつ対応付け分子として、小麦胚芽の系でも対応付け分子の形成を行うことができる。また、これまでウサギ網状赤血球の系では遺伝子型分子の安定性を欠くために実用性に乏しく、短い鎖長の遺伝子型分子にのみ適用されてきたがRobertsRWSzostakJWProcNatlAcadSciUSA9412297NemotoNMiyamotoSatoEYanagawaHFEBSLett414405、PEG領域を含むスペーサー部をもつ対応付け分子は、小麦胚芽の系ではより安定であり長い鎖長を取り扱える実用的な系である。無細胞共翻訳で相互作用が実現される場合、その無細胞翻訳系については、大腸菌、ウサギ網状赤血球、小麦胚芽の系などいずれでも構わない。invitrovirus法IVV法では、対応付け分子の形成は、大腸菌はかなり不安定であるが、ウサギ網状赤血球の系FEBSLett414405RobertsRWSzostakJWProcNatlAcadSciUSA9412297では安定に確認されており、さらに小麦胚芽の系ではより安定に確認されている。選択工程本工程は、対応付け工程で調製されたmRNAタンパク質対応付け分子ライブラリーから、標的物質と相互作用するmRNAタンパク質対応付け分子を選択する選択工程である。選択の方法は、対応付け分子を用いることができる限り特に制限されず、公知の方法に従って行うことができる。例えば、特許文献1に記載の淘汰プロセスが挙げられる。具体的には次のような方法が挙げられる。標的物質のタンパク質や核酸DNAまたはRNAや他の物質、例えば糖質や脂質などをマイクロプレートやビーズに予め共有結合や非共有結合を介して結合させておき、これに対応付け工程で調製した対応付け分子を加え、ある温度条件で、一定時間反応させた後、洗浄し、標的に結合しない対応付け分子を除去する。その後、標的に結合した対応付け分子を遊離させ、回収する。増幅工程本工程は、選択工程で選択されたmRNAタンパク質対応付け分子のmRNA部分に基づく核酸増幅により、候補タンパク質をコードするcDNAライブラリーを調製する増幅工程である。増幅の方法は、対応付け分子のmRNA部分に基づく核酸増幅が可能な限り特に限定されず、公知の方法に従って行うことができる。例えば、RTPCRにより逆転写及び増幅をする方法が挙げられる。本発明において、最初のcDNAライブラリーの調製及び繰返しの増幅工程におけるcDNAライブラリーの調製のうち複数回のcDNAライブリーの調製において、その調製回に特異的な配列を有するプライマーを用いて前記cDNAライブラリーを調製する。これにより、ライブラリーを構成する核酸部分の配列を決定したときに、いずれのライブラリーに由来するかを判別できるようにする。特異的な配列の長さや位置は、後工程の配列決定後に識別できれば特に制限されない。通常のスクリーニングでは上記~の工程が3~5回繰り返されるので、最初の調製を含めて6回分が識別できる長さであればよく、この場合には3塩基長あれば十分となる。一方、識別を容易にするために、論理的に必要な塩基数より長い配列を用いてもよい。このような観点から、通常には、4~10塩基長の配列を用いる。プライマーにおけるこの配列の位置は、増幅を妨げない限り特に制限されず、通常には、プライマーの5側に付加される。調製回に特異的な配列を有するプライマーを用いるcDNAライブラリーの調製は、最初及び全てのラウンド毎に行う必要はなく、後述する検出において有意な結果が得られる数だけ行えばよい。通常には、調製回に特異的な配列を含む3~6のcDNAライブラリーを準備すればよい。次に、特定のプライマーに含まれる配列を利用して配列を解析する部分について説明する。本発明では、調製回に特異的な配列を有するプライマーを用いて調製されたcDNAライブラリーを混合し、混合されたcDNAライブラリーの配列決定を行う。cDNAライブラリーの混合比は、配列決定に用いるシークエンサーの性能を考慮して決定すればよいが、目的配列が濃縮されていない早い調製回のライブラリーの量を多くすることが好ましい。すなわち、cDNAのライブラリーを混合する量比において、早い調製回の量を次の調製回の量以上とすることが好ましい。調整回ごとのプレイの濃縮率は通常10~100倍程度の範囲に収まるため、通常は配列決定を行う規模の大きさを考慮しながら、先の調製回の量を次の調製回の量の1~100倍とすることが好ましい。このようにすることにより解析を効率的に行うことができる。配列決定の方法は、混合物に含まれる配列調製回に特異的な配列を含めて、タンパク質をコードする配列を決定できる限り、特に制限されないが、混合により対象となる分子の数が増えるため、いわゆる次世代シークエンサー又はハイスループットシーケンシングに使用されるシーケンサーを用いる方法が好ましい。このように決定された配列を、調製回に特異的な配列に基づいて、調製回毎における、同一候補タンパク質をコードする配列の数分子の数に相当するを測定する。これは、いわゆるインシリコ解析と知られているコンピューターによる方法により行うことができる。このように調製回毎に測定された数に基づき、調製回に応じて有意に増加している配列分子を特定し、その配列にコードされる候補タンパク質を、前記標的物質と相互作用するタンパク質として検出する。増加の判定においては、調製回の一部のデータを用いてもよい。例えば、最初及び1~4ラウンドのデータがある場合に、後半の2~4ラウンドのデータのみを用いてもよいし、最初、2及び4のように中間を省略したデータを用いてもよい。配列分子の増加は、統計学的方法により評価することが好ましい。このような統計学的手法としては、フィッシャーの正確確率検定が挙げられる。標的物質は、例えば、タンパク質である。タンパク質の場合は、タンパク質間相互作用を検出できる。従来のIVV法などのmRNAディスプレイ技術を用いたタンパク質間相互作用解析技術では、偽陽性の判別のために、定量リアルタイムPCRアッセイなどの追加的な検証実験に大きなコストを払う必要があったが、本発明はそのコストを軽減できる。従来の方法では、プレイタンパク質の種別決定にラウンド毎の配列決定の結果を用いているため、その際のクローニング工程が解析効率を下げることになり、多くの偽陰性が生じていたが、本発明ではそれらの偽陰性を軽減できる。本発明によれば、ラウンド毎に配列決定をするのではなく、複数のラウンドのライブラリーを混合して配列決定をするとともに、ラウンド毎に配列を評価できるようにすることにより、ラウンド毎の変化の検出の誤差を最小にすることができ、それにより、IVV法によるスクリーニングにおける濃縮を正確に評価できるようになると考えられる。以下、具体的に本発明の実施例を記述するが、下記の実施例は本発明についての具体的認識を得る一助とみなすべきものであり、本発明の範囲は下記の実施例により何ら限定されるものでない。本実施例では、本発明を具体的に適用した事例に加え、従来のラウンド毎に配列決定を行う方法との比較結果について記載する。1概要IVV法によるセレクション過程を4ラウンド繰り返し、Irf7Interferonregulatoryfactor7と相互作用するタンパク質をマウスの脾臓由来のランダムフラグメントライブラリから抽出する試みを行った。2IVV法によるセレクションプレイライブラリーの準備mRNAディスプレイに用いるRNAライブラリーは、以前に報告された方法GenomeRes15710717MiyamotoSatoEetalNucleicAcidsRes31e78に従って、C57BL6マウスから抽出したpolymRNAから作成した。まず、polymRNAのランダムプライミングを行い、ライゲーションを介した増幅、そしてinvitro発現を行うためのマルチステップPCRを行った。結果のPCR産物SP6ΩT7FlagmentKpn1FLAGA8はQIAquickPCRPurificationKitで精製し、m7GpppGRNACappingAnalogと共にRiboMAXLargeScaleRNAProductionSystemSP6を用いて転写した。RNeasy96BioRobot8000Kitを用いてmRNAの精製を行った後、PEGPuroスペーサーを、T4RNAリガーゼPromegaを用いてmRNAの3末端にライゲーションし、再び精製した。ベイト調製ベイト用のDNAは、以前の報告に従い準備したGenomeRes15710717MiyamotoSatoEetalNucleicAcidsRes31e78。マルチステップPCRによって作成されたcDNAコンストラクトの構造T7Irf7CBPzzHisには、マウスのIrf7のコード領域全体が含まれている。QIAquickPCRPurificationKitを用いて生成されたPCR産物は、RiboMaxLargeScaleRNAProductionSystemSP6を用いてm7GpppGRNACappingAnalogと共にmRNAに転写され、RNeasy96BioRobot8000Kitを用いて精製された。Invitro翻訳とセレクションプレイmRNAライブラリーとベイトのIrf7のmRNAは、小麦胚芽抽出を用いて無細胞翻訳系にて共翻訳した。このとき、mRNAの3’末端に付加してあるPuromycinは、mRNAから翻訳されたタンパク質のC末端と共有結合で連結され、IVV分子が完成する。その後、タグ付けされたベイトタンパク質を用いた1段階4ラウンド目のみ2段階の精製RigautGetalNatBiotechnol1710301032を行った。各セレクションラウンドの後、プレイのmRNAはRTPCRによって増幅され、その後転写翻訳反応が行われ、次のラウンドのセレクションが再び行われる。ネガティブコントロール実験は、ベイトタンパク質が不在であること以外は同様の条件で行われた。454GSFLXによる配列決定のためのサンプル準備454GSFLXによるシーケンシングの後に配列の由来ラウンドを見分けるために、初期~4thラウンドのmRNAライブラリーのRTPCRによる増幅においては、次の表1のように4塩基のユニーク配列を持つバーコードプライマーを使用した。バーコードプライマーで増幅されたDNAサンプルは、QIAquickPCRPurificationKitを用いて精製した。DNAの精製後、それぞれのサンプルの濃度を、NanoDropを用いて測定した。最後に、初期及び1~4ラウンド目までのサンプルを2010721濃度比の割合で混合した。ネガティブコントロール実験におけるサンプルは、同様の方法で作成した。最後に、初期ライブラリーのサンプルと2つのミクスチャーサンプルをGSFLX454Rocheシーケンサーにかけ、配列決定を行った。シーケンシングの結果、ベイトライブラリについて206322リード、とベイトライブラリについて304504リード、イニシャルライブラリについて277833リードを得た。これらのうち、177935、278816、238683はバーコード領域によってセレクションラウンドを決定することができた。これらのリードはBLATを用いてマウスゲノムにマッピングされ、結果として、ベイトライブラリについて47849、ベイトライブラリについて63306、イニシャルライブラリにいて102092リードがゲノムに適正にヒットすることが確認され、それらを続くインシリコ解析に用いることとした。22インシリコ解析ゲノムへのマッピングによる領域のカウントと統計解析最初に、cDNA配列の端に存在するラウンド特異的な4塩基の部位のデコードを行い、各リードがどのラウンド由来のものであるのかを判定した。その後、処理されたリードはマスクされたゲノム配列に対してBLATを用いて次の条件の下でマッピングを行った:matchlength≧30bpandidentity≧95。マップされたリードのそれぞれのセレクションラウンドにおける頻度は、ゲノム上の各塩基ポジションについて計算され、初期ライブラリーとネガティブコントロールの結果に対して比較が行われた。さらに、初期ライブラリーおよびネガティブコントロールよりも高い頻度を示したポジションについては、統計的検定が行われたFisherの正確確率検定。統計的に有意性p0001のある差が、初期ライブラリーとネガティブコントロールと比べて認められた場合、その領域をベイト特異的に濃縮された領域とみなした。有意な濃縮を示した領域のうち、連続しているものは一つに統合され、一つの相互作用領域InteractionregionIRとみなした。リアルタイムPCRとの比較による定量性およびその有効性の確認本手法の定量性能を確認するために、複数のラウンドでリードが得られており、RefSeq領域との重なりが認められた21個の領域をランダムに抽出し、下記のようにリアルタイムPCRアッセイを行った。まず、ゲノム上でのリードの頻度がセレクションラウンドごとに計算された。これらの頻度はゲノム配列にマッピングされたリードの数に基づく。図2は、シーケンサーおよびリアルタイムPCRで定量された分子数の間の比較の例を示している。目視においては、本手法によるデータはリアルタイムPCRによる定量結果とよく符合した結果を見せたが、さらなる確証を得るために、それらの二つのデータの間の相関係数を計算し、それが統計的に有意に高いことを確認した。また、本手法で陽性と判定されたもののうち、87578がリアルタイムPCRでも陽性と判定され、リアルタイムPCRで偽陽性と判定されたもののうち、本手法で陽性であると判定されたものは、8112のみであることを確認した。リアルタイムPCRアッセイリアルタイムPCRアッセイは7300RealTimePCRSystemを用い、標準プロトコルに従って行った。その際、プレイライブラリー由来の5ngのDNAテンプレート、10μM05μlのプライマーフォワード側、リバース側それぞれ、125μlのPowerSYBRGreenPCRMasterMixを含む25μlのリアクションmixが使用された。サンプルの測定は2回ずつ行われ、その数値の平均を採用した。ハイスループットな相互作用領域の決定上に示したように、本手法で得られたリード数とリアルタイムPCRで計られた分子数の間に有意な相関が確認された。したがって、本手法の定量性能を生かして、システマティックに統計的に有意な濃縮領域を決定した。濃縮の統計的な有意性は2×2分割表に対するFisherの正確確率検定を用いて行った。それぞれの分割表は、ある領域についてのあるセレクションラウンドにおけるベイト、ベイト実験で得られたリード数、および各実験でよって得られた総リード数から成り立っている。加えて、イニシャルライブラリとあるラウンドとの比較も同様の方法で行われた。結果として、pvalue0001の条件でゲノム上の466個のプレイ領域が決定された。23従来法との比較偽陽性の軽減効果上記のランダムに選択した領域ただし、複数のラウンドでリードが得られていて、RefSeq領域にプライマーの設計が可能という条件付きのリアルタイムPCRの判定結果に基づけば、429921がライブラリーにおける陽性となる領域の割合と考えることができ、ラウンド毎のシーケンスの結果をリアルタイムPCRで判定する場合従来法の結果として考えることが出来る。本実施例では、本発明方法IVVHiTSeq法を介して陽性と判断されたもののうち、定量リアルタイムPCRでも陽性と判断されたものの割合は875であった。したがって、IVVHiTSeq法を用いることで従来法と比較して、高い信頼性を持って相互作用分子を特定することができると考えられる。偽陰性の軽減効果Irf7をベイトとして同一ライブラリーに対して、ラウンド毎の増幅物を用いる配列決定法で決定した740配列のうち、640配列865がIVVHiTSeq法によっても得られていたことから、従来法による結果の殆どを包含していることが示された表2、図4。また、IVVHiTSeq法によって陽性判定を得た479領域のうち、従来法で得られていたのは72領域のみであった。これは、残りの423領域がIVVHiTSeq法でないと得られないことを意味しており、IVVHiTSeq法の検出能力の高さを示している。更に、重なりがあった72領域も、ラウンド毎の増幅物を用いる配列決定によって得られたクローンだけからは、その真偽を確認することは出来ないため、定量リアルタイムPCRアッセイなどの検証を追加的に行う必要がある。本実施例は4790479×ラウンド数×2回以上のリアルタイム定量PCRアッセイを行ったことに相当し、ラウンド毎の配列決定とその後の定量リアルタイムPCRによる解析に比べてポストセレクションの効率が飛躍的に高くなっていることを示している。