JP4034740B2 - Ｄｎａ合成用のプライマーの選定方法 - Google Patents

Description

技術分野
本発明は，選択的にＤＮＡを合成するためのプライマーを選定する方法に関する。
背景技術
病原体およびその他の疾病の診断を遺伝子を用いて行う場合（遺伝子診断），検査試料中に含まれるＲＮＡからｃＤＮＡを合成し，さらにＰＣＲ法等を用いてＤＮＡを増幅した後にこれを検出する方法が用いられる。ＤＮＡの増幅はＲＮＡから逆転写反応によりｃＤＮＡを生成させた後に行うため，逆転写反応を開始させるプライマーの塩基配列が検出感度および効率に大きく影響を与える。
従来は，このプライマーとして，検査対象の遺伝子・病原体に特異的なプライマー（遺伝子特異的プライマー），哺乳動物のメッセンジャーＲＮＡ（ｍＲＮＡ）の多くが持つｐｏｌｙ（Ａ）配列に結合するｍＲＮＡ用プライマー（オリゴｄＴプライマー）および全てのＲＮＡと非特異的に結合するプライマー（ランダムヘキサマー）の三種が用いられてきた。これらのうち，遺伝子特異的プライマーは特定の遺伝子を逆転写させる効率が高いが，対象以外の遺伝子は逆転写されないため，多くの遺伝子を対象とする検査を行う場合には，多種の遺伝子特異的プライマーを調製する必要があるため多額のコストを要するという問題点がある。例えば，大腸菌の特定の株に特異的なプライマーセットとして４０００種類のプライマーの混合物が用いられている。このような特異的プライマーを用いる検出方法は，特に食中毒患者からの試料の検査など対象微生物が不明である場合には実用的ではない。また，遺伝子特異的プライマーを用いた場合には，検出可能な塩基配列が限定されるため，塩基配列が若干異なる可能性のある多数の遺伝子を検出することが困難である。特に，細菌類やウイルスは種間で遺伝子の配列が若干異なる上に，頻繁に遺伝子の突然変異を起こすため，これらの変種を検出しうる遺伝子特異的プライマーの選定は非常に難しい。一方，オリゴｄＴプライマーは，ｐｏｌｙ（Ａ）配列を持つｍＲＮＡに対しては比較的高い効率で検出が可能であるが，ｐｏｌｙ（Ａ）配列を持たないｍＲＮＡとして転写される遺伝子を検出することができないという問題点がある。また，非特異的にＲＮＡからｃＤＮＡを合成するランダムヘキサマーを用いた場合には，細胞内のＲＮＡ中で大きな割合を占めるリボゾームＲＮＡや高度リピート配列に対応するｃＤＮＡが大量に合成されてしまうため，目的とする遺伝子・病原体の検出効率が低いという問題点がある。
したがって，本発明の目的は，ランダムヘキサマーと同様に，細菌，ウイルス，寄生虫，真菌などに由来するＲＮＡを非特異的に逆転写するが，リボゾームＲＮＡおよび高度リピート配列を逆転写する頻度が低いプライマーの混合物（制限プライマーミックス）を選定する方法を提供することを目的とする。
発明の開示
本発明は，ＤＮＡ合成用のプライマーを選定するプライマー選定方法を提供する。該方法は，プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，該頻度情報に基づき，該非標的配列に出現する頻度が所定の低頻度閾値以下のモチーフをプライマーとして選定するプライマー選定ステップと，を含むことを特徴とする。
別の観点においては，本発明のプライマー選定方法は，プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度および標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，該頻度情報に基づき，該非標的配列に出現する頻度が低く，該標的配列に出現する頻度が高い所定の選択特性をもつモチーフをプライマーとして選定するプライマー選定ステップと，を含むことを特徴とする。
好ましくは，該プライマー選定ステップは，モチーフが該非標的配列に出現する頻度が所定の低頻度閾値以下であり，該標的配列に出現する頻度が所定の高頻度閾値以上であるとき，そのモチーフが該所定の選択特性をもつと判定することを特徴とする。また好ましくは，該プライマー選定ステップは，モチーフが該非標的配列に出現する頻度と該標的配列に出現する頻度の相違に基づいて，そのモチーフが所定の選択特性をもつか否かを判定することを特徴とする。
本発明のプライマー選定方法の特に好ましい観点においては，非標的配列はリボゾームまたは高度リピート配列の塩基配列である。また好ましい観点においては，標的配列は正常細胞に存在しない異常ＲＮＡ，例えば，細菌，ウイルス，寄生虫または真菌に由来するＲＮＡの塩基配列である。
別の観点においては，本発明は，ＤＮＡ合成用のプライマーを選定するプライマー選定装置であって，プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度を示す頻度情報を求める頻度情報生成部と，該頻度情報に基づき，該非標的配列に出現する頻度が所定の低頻度閾値以下のモチーフをプライマーとして選定するプライマー選定部と，を含むことを特徴とするプライマー選定装置を提供する。好ましくは，本発明のプライマー選定装置は，プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度および標的配列に出現する頻度を示す頻度情報を求める頻度情報生成部と，該頻度情報に基づき，該非標的配列に出現する頻度が低く，該標的配列に出現する頻度が高い所定の選択特性をもつモチーフをプライマーとして選定するプライマー選定部と，を含む。
さらに別の観点においては，本発明は，ＤＮＡ合成用のプライマーを選定するための処理をコンピュータに実行させるためのプログラムであって，プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，該頻度情報に基づき，該非標的配列に出現する頻度が所定の低頻度閾値以下のモチーフをプライマーとして選定するプライマー選定ステップと，を該コンピュータに実行させることを特徴とするプログラムを提供する。好ましくは，本発明のプログラムは，プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度および標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，該頻度情報に基づき，該非標的配列に出現する頻度が低く，該標的配列に出現する頻度が高い所定の選択特性をもつモチーフをプライマーとして選定するプライマー選定ステップと，を該コンピュータに実行させることを特徴とする。
別の観点においては，本発明は，ＤＮＡ合成用のプライマーを選定するための処理をコンピュータに実行させるためのプログラムを記録したコンピュータ読取可能記録媒体を提供する。
さらに別の観点においては，本発明は，上述の本発明のプライマー選定方法により選定されたプライマーの配列を有するオリゴヌクレオチド混合物，ならびに該混合物を含むｃＤＮＡ合成用キットを提供する。本発明はまた，上述の本発明のオリゴヌクレオチド混合物をプライマーとして，逆転写酵素またはＤＮＡ合成酵素を用いてＤＮＡを合成することを特徴とするＤＮＡの製造方法を提供する。このようにして合成されたＤＮＡは，ＰＣＲのテンプレートとして，およびオリゴヌクレオチドのアレイを用いるハイブリダイゼーションアッセイ系のプローブとして用いることができる。本発明はまた，試料において正常細胞に存在しないＲＮＡを検出する方法であって，試料からＲＮＡを抽出し，請求項２３記載のオリゴヌクレオチド混合物をプライマーとして逆転写によりｃＤＮＡを合成し，そして，該ｃＤＮＡ中に正常細胞に存在しないＲＮＡに対応する配列が存在するか否かを判定することを特徴とする方法を提供する。正常細胞に存在しないＲＮＡは，好ましくは，細菌，ウイルス，寄生虫または真菌に由来するＲＮＡである。
本発明にしたがえば，ランダムヘキサマーと同様に，細菌，ウイルス，寄生虫，真菌等に由来するＲＮＡを非特異的に逆転写するが，リボゾームＲＮＡおよび高度リピート配列を逆転写する頻度が低いプライマーを容易に選定することができる。本発明にしたがって選定されたプライマーの混合物（制限プライマーミックス）を用いて逆転写反応を行うことにより，逆転写産物中に占める検出対象のｃＤＮＡの比率が増加するため，ヒトや家畜に感染した細菌，ウイルス，寄生虫，真菌等を高い感度で検出することができる。また，本発明の制限プライマーミックスを用いることにより，このような細菌やウイルスの遺伝子配列が不明の場合であっても，正常細胞に存在しない異常ＲＮＡを高感度で検出することができる。本発明の制限プライマーミックスは，特に，食中毒患者からの糞便試料を用いて未知の病原体を検出するアッセイにおいて有用である。
発明を実施するための最良の形態
以下，本発明の好適な実施の形態を図面を参照して説明する。
プライマーの選定方法
図１は，本実施の形態におけるプライマー選定装置のハードウエア構成例を示している。プライマー選定装置１は，本発明のプライマー選定方法を実行する。
図１に示されるように，プライマー選定装置１は，ＣＰＵ３，ＲＯＭ５，ＲＡＭ７，ハードディスク９，ＣＤ−ＲＯＭ等の媒体装着部１１，キーボード等の入力装置１３，ディスプレイ，プリンタ等の出力装置１５および通信装置１７を含む。
図１から明らかなように，プライマー選定装置１は，汎用のコンピュータで構成されてよい。本発明のプライマー選定方法を実現するためのプログラムがコンピュータにインストールされる。このプログラムがＣＰＵ３により実行されて，プライマー選定装置１が実現される。プライマー選定装置１は，インターネット等の通信により他のコンピュータからアクセスされてもよい。そして，他のコンピュータからの要求に応えてプライマーを選定してもよい。
図２は，プライマー選定装置１の機能ブロック図である。図２において，プライマー長さ取得部２１は，プライマーの長さを取得する。この取得された長さをもつ適切なプライマーが，以降の解析処理にて求められる。プライマー長さは，例えば，オペレータによりキーボードを用いて入力される。また，プライマー長さは，通信により取得されてもよい。プライマー長さは，塩基長で表される。本実施形態では，プライマー長さを「分析塩基単位」とも称する。プライマー長さは，標的配列および非標的配列の性質に応じて，当業者が適宜選択することができる。好ましくはプライマー長さは４−３０塩基である。４塩基より短いとプライミングの効率が低く，３０塩基を越えるとプライマー合成のコストが高くなる。より好ましくはプライマー長さは６−２０塩基であり，最も好ましくは６−９塩基である。
モチーフ生成部２３は，プライマー長さに基づき，プライマー候補モチーフ群を生成する。モチーフは，プライマー長さ（分析塩基単位）をもつ塩基配列のパターンである。そして，プライマー候補モチーフ群は，与えられたプライマー長さにて実現可能なモチーフの集合である。例えば，塩基がＡ，Ｃ，Ｇ，Ｔの４種類であり，プライマー長さが６である場合，プライマー候補モチーフ群は，すべての可能な塩基配列を有する６塩基のモチーフ４^６（＝４０９６）個の集合である。プライマー候補モチーフ群は，以降の処理にて，プライマー候補のモチーフの集合として使われる。
配列情報取得部２５は，プライマー選定処理に使う塩基配列を取得する。ここでは，非標的配列および標的配列が取得される。非標的配列は，逆転写されるべきでない配列である。一方，標的配列は，優先的に逆転写されるべき配列である。例えば，細胞に感染したウイルスを検出するためのプライマーが選定されるき，非標的配列の例が細胞のリボゾーム配列であり，標的配列の例が検出すべきウイルスの配列である。配列の情報は，例えば，媒体装着部に装着された記録媒体（ＣＤ−ＲＯＭ等）から取得される。配列は，通信等の他の手法により，公共のデータベース等から取得されてもよく，ＤＮＡ配列決定装置に装備されるコンピュータから取得されてもよい。
頻度情報生成部２７は，プライマー候補モチーフ群の各モチーフが配列に出現する頻度を示す頻度情報を求める。頻度情報生成部２７は，非標的頻度情報生成部２９および標的頻度情報生成部３１を含む。非標的頻度情報生成部２９は，非標的配列（例えばリボゾーム）におけるモチーフ出現頻度を求める。標的頻度情報生成部３１は，標的配列（例えば検出すべきウイルス）におけるモチーフ出現頻度を求める。
プライマー選定部３３は，頻度情報生成部２７により求められた頻度情報に基づき，プライマー候補モチーフ群からプライマーを選定する。標的頻度情報生成部３１は，所定の選択特性をもつモチーフをプライマーとして選定する。選択特性は，後述するように，非標的配列に出現する頻度が低く，前記標的配列に出現する頻度が高いことを示すプライマーの特性である。典型例では，あるモチーフの非標的配列での出現頻度が所定の低頻度閾値以下であり，同モチーフの標的配列での出現頻度が所定の高頻度閾値以上のとき，そのモチーフが本実施形態の選択特性をもつと判断される。低頻度閾値および高頻度閾値の値は，検出の目的および対象の特性を考慮して，当業者が容易に設定するか，または実験的に求めることができる。
本実施形態では，プライマー選定部３３により選定されたプライマーの集合を，制限プライマーミックスという。プライマーセット出力処理部３５は，制限プライマーミックスの情報を出力するための処理を行う。典型的には，制限プライマーミックスの情報は，プリンターおよびディスプレイ等の出力装置から出力される。制限プライマーミックスの情報は，通信により他のコンピュータ等へ出力されてもよい。好ましくは，制限プライマーミックスの情報は，ＤＮＡ合成機の運転を制御するコンピュータに送られる。
図３は，図２のプライマー選定装置１によるプライマー選定処理の概要を示している。この例では，プライマー長さ（分析塩基単位）が６であり，非標的配列がリボゾームであり，標的配列がウイルスである。
図３に示されるように，モチーフ生成部２３は，６塩基の長さを有するモチーフの群を生成する。モチーフ群は，典型的には４種類の塩基のすべての組み合わせを含むが，プライマーの推定二次構造を考慮して，または非標的配列または標的配列に関する従来の知見に基づいて，モチーフ群から特定の配列を排除してもよい。頻度情報生成部２７の非標的頻度情報生成部２９は，モチーフ群の各モチーフがリボゾーム配列中に出現する頻度を求める。ここでは，リボゾーム配列中の任意の連続する６塩基からなる部分配列がモチーフ群の各モチーフと比較される。リボゾーム配列の始端から６文字の部分配列が選択され，モチーフ群を参照して，選択された６文字に該当するモチーフを選択し，該当するモチーフの頻度に「１」が加えられる。すなわち，頻度の値が１ずつ増加する。次に，１塩基ずらして部分配列が選択され，同様の処理が行われる。このようにして，非標的頻度情報生成部２９は，リボゾーム配列全体における各モチーフの出現頻度を求める。一方，頻度情報生成部２７の標的頻度情報生成部３１は，モチーフ群の各モチーフがウイルス配列中に出現する頻度を求める。この処理は，上述のリボゾームに関する処理と同様でよい。このようにして，プライマー選定装置１は，モチーフ群の生成および塩基配列中における各モチーフの出現頻度の集計を行う。
図４は，頻度の集計結果の一例を表す表である。プライマー選定部３３が，図４の頻度情報に基づき，モチーフ群からモチーフを選定する。プライマー選定部３３は，リボゾームに出現する頻度が低く，かつ，ウイルスに出現する頻度が高いモチーフを求める。
例えば，図４においては，モチーフ１は，リボゾームでの出現頻度が高く，ウイルスでの出現頻度が低い。また，モチーフ２は，リボゾームおよびウイルスでの出現頻度が高い。また，モチーフ３は，リボゾームおよびウイルスでの出現頻度が低い。これらのモチーフ１〜３は，プライマーとして選ばれない。一方，モチーフ４は，リボゾームでの出現頻度が低く，ウイルスでの出現頻度が高い。そこで，モチーフ４は，プライマーとして選択される。同種の特性をもつモチーフがさらに選択され，プライマーのセットが得られる。
上述のように，プライマー選定部３３は，リボゾームに出現する頻度が低く，ウイルスに出現する頻度が高いモチーフをプライマーの配列として選定する。出現頻度の大きさは，このプライマーを用いて逆転写反応を実施したときの逆転写の頻度の大きさと相関する。すなわち，あるモチーフのリボゾーム中での出現頻度が低いとき，そのモチーフを有するプライマーからリボゾーム配列が逆転写される頻度も低い。あるモチーフのウイルス中での出現頻度が高いとき，そのモチーフを有するプライマーからウイルス配列が逆転写される頻度も高い。したがって，上述の処理により，プライマー選定部３３は，ウイルスに由来する配列を優先的に逆転写するプライマーで構成される適切なプライマーセットを選定することができる。
図５は，頻度情報の生成および処理の例を示すフローチャートである。図５の処理は，非標的配列および標的配列のどちらの頻度を求めるときにも適用される。
まず，プライマー長取得部２１が分析塩基単位（プライマー長）を取得し，これにより分析塩基単位が決定される（Ｓ１０）。そして，モチーフ生成部２３が，分析塩基単位をもつ可能な全てのモチーフを自動生成する（Ｓ１２）。これにより，プライマー候補モチーフ群（モチーフ一覧）が得られる。
Ｓ１４では，配列情報取得部２５が，分析対象の配列の情報を取得する。分析対象の配列は，非標的配列または標的配列である。配列情報取得部２５は，二重鎖の片側の情報を取得する。Ｓ１６では，配列情報がテキストデータに変換される。そして，Ｓ１８では，テキストの始端から６文字の部分配列が選択される（分析塩基単位が６塩基である場合）。ここでは，この選択された部分配列を選択配列と呼ぶ。Ｓ２０では，モチーフ一覧を参照して，選択された６文字に該当するモチーフを選択する。Ｓ２２では，該当するモチーフの頻度に「１」が加えられる。すなわち，頻度の値が１ずつ増加する。
Ｓ２４では，選択配列が１文字シフトされる。Ｓ２６にて，選択配列の文字数が分析塩基単位（６）以上であるか否かが判定される。Ｓ２６がＹＥＳの場合，Ｓ２０に戻る。Ｓ２６がＮＯの場合，塩基配列の全範囲について頻度の取得が終了したことを意味する。そこで，Ｓ２６がＮＯの場合，Ｓ２８に進む。
次に，Ｓ２８では，反対鎖テキストが生成される。より詳細には，Ｓ１６では，頻度情報生成部２７は，二重鎖の片側の配列情報をテキストデータに変換した。Ｓ２８では，Ｓ１６のテキストデータを，反対鎖のテキストデータに変換する。ＡがＴに，ＴがＡに，ＣがＧに，ＧがＣに変換される。さらに，テキストの方向が反転される。
Ｓ３０〜Ｓ３８は，Ｓ１８〜Ｓ２６と同様の処理である。Ｓ３０〜Ｓ３８により，頻度情報生成部２７は，反対鎖テキストにおけるモチーフ出現頻度の集計結果を得る。Ｓ３８がＮＯであれば，反対鎖の全範囲を対象として頻度が求められている。すなわち，二重鎖の両側の配列について頻度の集計が得られている。そこで，Ｓ３８がＮＯであれば，頻度情報生成処理が終了する。なお，Ｓ２８〜Ｓ３８は，二重鎖を想定する必要がある場合に行われればよい。
以上，図５の頻度情報生成処理を説明した。上記の処理によれば，頻度情報生成部２７において，モチーフ一覧中の各モチーフの出現回数の集計結果が，頻度情報として得られる。既に述べたように，図５の処理は，非標的配列における頻度を求めるときにも，標的配列における頻度を求めるときにも適用される。また，図示されないが，頻度情報は，プライマー選定装置１の記憶装置に記憶される。
図６は，プライマー選定部３３がプライマーを選定する処理の一例を示すフローチャートである。プライマー選定部３３は，非標的配列における頻度の情報を参照し（Ｓ５０），モチーフ一覧に示されるモチーフ群から，所定の低頻度閾値以下の出現頻度をもつモチーフを選ぶ（Ｓ５２）。さらに，プライマー選定部３３は，標的配列における頻度の情報を参照し（Ｓ５４），所定の高頻度閾値以上の出現頻度をもつモチーフを選ぶ（Ｓ５６）。プライマー選定部３３は，Ｓ５２およびＳ５６の両方で選択されたモチーフのセットを，本実施形態の制限プライマーセットとして特定する（Ｓ５８）。このプライマーセットに含まれるプライマーは，非標的配列を逆転写しにくく，標的配列を逆転写しやすい傾向をもつ。このようにして，本実施形態の好適な制限プライマーセットが得られる。
図７は，プライマー選定部３３がプライマーを選定する処理のもう一つの例を示すフローチャートである。プライマー選定部３３は，各モチーフについて非標的配列および標的配列における頻度の情報を参照し（Ｓ７０），相違パラメータを計算する（Ｓ７２）。相違パラメータは，標的配列での出現頻度と非標的配列での出現頻度の相違を表すパラメータである。相違パラメータは，非標的配列での出現頻度に対する標的配列での出現頻度の比であってもよく，標的配列での出現頻度と非標的配列での出現頻度の差であってもよい。後者の場合，差を求める前に標的配列と非標的配列の長さの相違にしたがって出現頻度の値を出現確率の値に補正する必要がある。あるいは，相違パラメータは，標的配列および非標的配列における出現頻度の順位の差または比であってもよく，その他の関数により求めてもよい。
図７の例では，相違パラメータは，非標的配列での出現頻度に対する標的配列での出現頻度の比である。Ｓ７２では，プライマー選定部３３は，相違パラメータである出現頻度の比に基づき，モチーフ一覧に示されるモチーフ群からプライマーを選ぶ（Ｓ７４）。出現頻度の比が大きいとき，非標的配列での頻度が低く，かつ，標的配列での頻度が高い。すなわちこのモチーフは本実施形態の選択特性をもつ。そこで，Ｓ７４では，プライマー選定部３３は，所定の頻度比閾値以上の出現頻度比をもつモチーフを選ぶ。選ばれたモチーフのセットが，制限プライマーセットとして特定される（Ｓ７６）。このようにして，プライマー選定部３３は，非標的配列を逆転写しにくく，標的配列を転写しやすい傾向をもつプライマーを選択できる。したがって，図７の処理によっても，本実施形態の好適な制限プライマーセットが得られる。
ここで，図６および図７において，プライマー選定部３３は，各モチーフが選択特性を有するか否かを判定している。ここでいう選択特性とは，非標的配列に出現する頻度が低く，標的配列に出現する頻度が高い特性である。この選択特性が本実施形態ではプライマー選定基準として用いられている。そして図６では，低頻度閾値および高頻度閾値を用いて選択特性の判定が行われた。また，図７では，相違パラメータを用いて選択特性の判定が行われた。
しかし，本発明はこれらの処理に限定されない。プライマー選定装置１は，以下に例示される処理により，適当な選択特性をもつモチーフを求めてもよい。
例えば，出現頻度は，配列中での出現回数以外の数値で表されてもよい。例えば，出現頻度は出現確率で表されてもよい。また，出現頻度は，正規分布の偏差値のような別の統計的数値で表されてもよい。その他にも，出現頻度は，頻度の高低を示す別のパラメータで表されてもよい。
また，本実施形態では，モチーフの選定のために，頻度の順位を用いることも可能である。例えば，非標的配列での頻度順位が所定の上限値以下であり，標的配列での頻度順位が所定の下限値以上のとき，プライマー選定装置１は，モチーフが適当な選択特性をもつと判断する。
また，図７の処理は，標的頻度と非標的頻度での出現頻度の相違を表す相違パラメータを用いた。図７の例では，相違パラメータは，出現頻度の比であるが，例えば，相違パラメータは出現頻度の差でもよい。標的配列での出現頻度と非標的配列での出現頻度の差が所定の値以上のとき，プライマー選定装置１が，モチーフが適当な選択特性をもつと判断する。ここでは，出現確率の差が好適に用いられる。あるいは，頻度順位に基づいて相違パラメータを決定してもよい。
図８は，本発明の別の実施形態におけるプライマー選定装置１の動作を示すフローチャートである。上述の実施形態では，予めプライマーの長さが決められていた。一方，本実施形態では，プライマー選定装置１は，適切な塩基長をもつプライマーを求める。プライマー選定装置１の構成は，既に述べた実施の形態（図１および図２と同様でよい）。
図８を参照すると，プライマー選定装置１は，まず，プライマー長さを初期値ｎ０に設定する（Ｓ８０）。初期値ｎ０は，例えば５である。そして，プライマー選定装置１は，頻度情報を生成する（Ｓ８２）。Ｓ８２は，図５に示された処理でよい。さらに，プライマー選定装置１は，頻度情報からプライマーを選定する（Ｓ８４）。Ｓ８４は，図６または図７に示された処理でよい。
次に，プライマー選定装置１は，Ｓ８４で得られたプライマーの数を，所定の要求プライマー数と比較する（Ｓ８６）。そして，プライマー選定装置１は，プライマーの数が，要求プライマー数以上か否かを判定する。Ｓ８６がＮＯであれば，プライマー選定装置１は，プライマー長さを一つ延ばし（Ｓ８８），Ｓ８２に戻る。Ｓ８６がＹＥＳであれば，処理を終了する。
このように，本実施形態によれば，プライマー選定装置１は，複数のプライマー長さの各々について，頻度情報に基づきプライマーセットを選定する処理を行う。そして，プライマー選定装置１は，要求数を満たす数のプライマーが得られる範囲で，最小の長さのプライマーで構成されるプライマーセットを求める。
したがって，本実施形態によれば，短いプライマーで構成され，かつ，適当な数のプライマーをもつプライマーセットが得られる。プライマーが短いので，オリゴヌクレオチド合成のコストを削減でき，かつ，適当な数のプライマーが得られるので，逆転写の効率もよい。
以上，非標的配列としてリボゾームを例として本発明のプライマー選定方法について説明してきたが，当業者は，非標的配列として例えば高度リピート配列を用いた場合にも同様にして適切なプライマーを選定しうることを理解するであろう。
さらに，本発明のプライマー選定方法は，標的配列，例えばウイルス配列が未知である場合にも適用することができる。この場合，頻度情報生成部２７は非標的頻度情報生成部２９からなり，ここでは非標的配列の頻度情報のみが求められる。次に，プライマー選定部３３は，非標的頻度情報のみに基づいてプライマーを選定する。選定は，例えば，モチーフが非標的配列に出現する頻度が所定の低頻度閾値以下である場合に，このモチーフをプライマーとして選定することにより行うことができる。
また，本発明のプライマー選定方法は，標的配列，例えばウイルス配列が複数である場合にも適用することができる。この場合，プライマーは，上述の配列未知の場合と同様に，非標的頻度情報のみに基づいて選定されてもよく，あるいは，それぞれの標的に対応する複数の標的頻度情報を参照することにより選定されてもよい。例えば，各標的について図５−８に記載されるステップにより適切なプライマーを選定した後にこれらのプライマーを混合してもよく，あるいは，各標的についての標的頻度情報を統合して数学的または統計学的に処理した後に，非標的頻度情報と比較することによりプライマーを選定してもよい。
プライマーの製造および使用
上述のようにして選定されたプライマーの配列にしたがって，それぞれの配列を有するオリゴヌクレオチドの混合物（制限プライマーミックス）を製造する。オリゴヌクレオチドの合成方法は当該技術分野においてよく知られており，例えば汎用のＤＮＡ合成装置（例えば３９４合成器，Ａｐｐｌｉｅｄ Ｂｉｏｓｙｓｔｅｍｓ社）を用いて製造することができる。
次に，この制限プライマーミックスを用いて生化学的手段によりｃＤＮＡを合成することができる。被検材料とその対照材料からＲＮＡを抽出後，上述のように作成した制限プライマーミックスをプライマーとして逆転写酵素によりｃＤＮＡを合成する。合成されたｃＤＮＡは，効率の高いＰＣＲ検査において用いることができる。以下に，ｃＤＮＡ合成の詳細を説明する。
ＲＮＡを抽出する材料としては，細菌が感染した組織，ウイルスが感染した細胞およびｐｏｌｙ（Ａ）配列を持たないＲＮＡを検出する必要がある細胞などがあげられる。
ＲＮＡを抽出するためには，グアニジン−イソチオシアネート混合液または適当な界面活性剤，例えばＳＤＳ，ＮＰ−４０，トリトンＸ−１００，デオキシコール酸等を用いて，ホモジナイサーを用いた方法や凍結融解等の物理学的方法によって，細胞を部分的または完全に破壊，可溶化する。次に，染色体ＤＮＡを，ポリトロン等のミキサーもしくは注射筒を用いて，ある程度せん断し，その後，蛋白質と核酸分画とを分別する。この分別操作には，特にフェノール・クロロホルム抽出もしくは超遠心を用いるＣｓＣｌ重層法等が一般的に用いられる。
また，上記方法においては，ＲＮａｓｅによるＲＮＡの分解を防ぐために，ＲＮａｓｅインヒビター，例えばＬｉｆｅ Ｔｅｃｈｎｏｌｏｇｉｅｓ社のＲＮａｓｅ Ｉｎｈｉｂｉｔｏｒ，Ａｍｂｉｏｎ社のＳｕｐｅｒ ＲＮａｓｅＩｎｎ等を添加しておくのが良い。
抽出されたＲＮＡを鋳型としたｃＤＮＡの合成は，鋳型としてＲＮＡ，プライマーとして上述のように合成した制限プライマーミックス，ｄＮＴＰ（ｄＡＴＰ，ｄＧＴＰ，ｄＣＴＰおよびｄＴＴＰ）混合物の存在下で，逆転写酵素を用いてＲＮＡに相補的な一本鎖ｃＤＮＡを合成する。逆転写酵素としては，例えば，ＳｕｐｅｒＳｃｒｉｐｔＩＩ（ＲＮａｓｅＨ⁻）（Ｉｎ Ｖｉｔｒｏｇｅｎ社）を用い，添付の指針に記載の条件下で逆転写反応を行うことができる。
このようにして選択的に合成されたｃＤＮＡを用いて，ＰＣＲ法等により，細菌・ＲＮＡウイルス等の病原体を検出することができる。また，細菌・ＲＮＡウイルス等と被検材料の遺伝子を含む複数の対象遺伝子についてＰＣＲによる検査を実施することができる。たとえば，このように合成したｃＤＮＡの一部を用いて，口腔内の炎症を起こした細胞の標本について癌関連の遺伝子の異常の有無を調べた後，その結果として細菌・ＲＮＡウイルス等を検査する必要が生じた場合に残りのｃＤＮＡを用いてこれらを検査することができる。
また，本発明のプライマーミックスを用いて作成されたｃＤＮＡは，非特定の塩基配列を検出するためにも用いることができる。細菌が感染した組織やＲＮＡウイルスの感染細胞と非感染の対照細胞について，それぞれ上述のようにｃＤＮＡを合成した後，二重鎖のｃＤＮＡを合成し，サブトラクション法の１つとして知られるＲＤＡ法（Ｒｅｐｒｅｓｅｎｔａｔｉｏｎａｌ Ｄｉｆｆｅｒｅｎｃｅ Ａｎａｌｙｓｉｓ）に基づく方法を実施することにより，細菌やウイルスの特定の塩基配列に依存せずにこれらの塩基配列を増幅することができ，このことにより未同定の病原体を検出することが可能となる。
上述の逆転写反応終了時のｃＤＮＡは，ＲＮＡと結合した状態であるが，これにＲＮａｓｅＨおよびＤＮＡポリメラーゼを反応させることによりＲＮＡをＤＮＡに置換して二重鎖のｃＤＮＡとする。この二重鎖のｃＤＮＡについてＲＤＡ法を実施する。ＲＤＡ法の概要を図９に示す。ＲＤＡ法とは，二重鎖のｃＤＮＡを，パリンドローム配列を認識し５’突出の断端を生成する制限酵素（例えばＡＧＣＴの配列を認識するＤｐｎＩＩ）で切断し，二種類のリンカーをＴ４ＤＮＡリガーゼおよび制限酵素で順次結合・切断を繰り返した後、二種の配列のヘテロ二重鎖をＰＣＲにより増幅させることにより、一方のｃＤＮＡにのみ存在する配列を選択的に増幅させる方法である。すなわち、図９に示すように，ＤｐｎＩＩ切断断片と結合する粘着末端を生じる２４塩基および１２塩基の二種のオリゴマーから構成されるリンカー１（５’−ＡＧＣＡＣＴＣＴＣＣＣＡＧＣＣＴＣＴＣＡＣＣＧＡＧ−３’（リンカープライマー１）および５’−ＧＡＴＣＣＴＣＧＧＴＧＡ−３’（リンカープライマー２））（Ｌｉｓｉｔｓｙｎ，ｅｔ ａｌ．（１９９３）Ｓｃｉｅｎｃｅ ２５９，９４６）を、ＤｐｎＩＩで切断されたｃＤＮＡとリガーゼ反応液中で混合し，リガーゼを加えて１４℃で保温することにより，ｃＤＮＡとリンカー１を共有結合（リンカーの連結）させる。このＤｐｎＩＩによる切断とリンカーの連結を、感染細胞由来のｃＤＮＡと非感染細胞由来のｃＤＮＡの双方について実施する（感染細胞および非感染細胞由来リンカー付ｃＤＮＡ）。次に、感染細胞および非感染細胞由来リンカー付ｃＤＮＡをＴａｑポリメラーゼを含むＴａｑポリメラーゼ合成混合液と混合し、７２℃で１０分間反応させることにより，２４塩基のリンカープライマー１の相補鎖をｃＤＮＡの３’末端から合成させる（感染細胞および非感染細胞由来リンカー配列合成ｃＤＮＡ）。感染細胞および非感染細胞由来リンカー配列合成ｃＤＮＡについて，２４塩基のリンカープライマー１をプライマーとしたＰＣＲ反応を行うことにより、感染細胞および非感染細胞由来のｃＤＮＡのＤｐｎＩＩ切断断片を増幅する（感染細胞および非感染細胞由来アンプリコン）。
ここで増幅された感染細胞および非感染細胞由来アンプリコンの両端にはＤｐｎＩＩで切断可能なリンカーが含まれているため，これらを切断する。切断後ｃＤＮＡと切り離されたリンカー１に相当する配列は，ミニカラムクロマトグラフィー（Ａｍｅｒｓｈａｍ社）またはシリカゲルを利用した低分子ＤＮＡ除去キット（ＢＩＯ１０１ ＧｅｎｅＣｌｅａｎＩＩ，Ｂｉｏｇｅｎ社）等によって除去する（ＤｐｎＩＩ切断感染細胞由来アンプリコンおよびＤｐｎＩＩ切断非感染細胞由来アンプリコン）。ＤｐｎＩＩ切断感染細胞由来アンプリコンについては，末端にリンカー２（５’−ＡＣＣＧＡＣＧＴＣＧＡＣＴＡＴＣＣＡＴＧＡＡＣＧ−３’（リンカープライマー３）および５’−ＧＡＴＣＣＧＴＴＣＡＴＧ−３’（リンカープライマー４））を、リンカー１を連結したときと同様の方法により、改めて連結する（リンカー２接続感染細胞由来アンプリコン）。リンカー２接続感染細胞由来アンプリコンを過剰量のＤｐｎＩＩ切断非感染細胞由来アンプリコンと混合した後，９９℃で加熱後緩やかに冷却することにより，感染細胞と非感染細胞由来のｃＤＮＡ（アンプリコン）についてヘテロ二重鎖を形成させる。このヘテロ二重鎖の形成においては，非感染細胞由来のｃＤＮＡ（アンプリコン）のモル比率が大きいため，感染細胞由来のｃＤＮＡ（アンプリコン）のうち非感染細胞に存在するｃＤＮＡ（アンプリコン）と塩基配列が相同の断片のほとんどは非感染細胞のｃＤＮＡ（アンプリコン）とヘテロ二重鎖を形成する。感染細胞由来アンプリコンのうち、ウイルスのＲＮＡから合成されたアンプリコンと相同の配列は、非感染細胞由来のアンプリコンには存在しないため、リンカー２接続感染細胞由来アンプリコン同士が二重鎖を形成する。結果としてウイルスに由来するアンプリコンのみがその両端にリンカープライマー３が連結している二重鎖を形成する。
つづいて、感染細胞および非感染細胞由来ｃＤＮＡのリンカー配列を合成した際と同様にＴａｑポリメラーゼによりリンカープライマー３の反対鎖を合成した後に，２４塩基のリンカープライマー３をプライマーとしてＰＣＲ反応を実施する。ＰＣＲ反応で効率良く増幅反応が進行するためには二重鎖ＤＮＡの両端にプライマーと相同の配列が存在する必要があるが、上記反応の場合には二重鎖の両端にリンカープライマー３と相同な配列を持つｃＤＮＡ（アンプリコン）のみが増幅する。非感染細胞には存在せず感染細胞にのみ存在するＲＮＡのうち主要なものは，細菌・ＲＮＡウイルス等に由来するＲＮＡであるため，このＲＤＡ法によって細菌・ＲＮＡウイルス等の断片が増幅される。これをアガロースゲル中で電気泳動した後，エチジウムブロミドにより染色することによりＤＮＡ断片を可視化して検出する。あるいは，特異的プライマーを用いてハイブリダイゼーションにより検出してもよく，適切なプラスミドに挿入した後シークエンシングしてもよい。このようにして，塩基配列特異的プライマーを用いることなく感染細胞中の細菌・ＲＮＡウイルス等を検出することができる。
すなわち，本発明の制限プライマーミックスを用いれば，ウイルス種を特定しないまま細菌・ＲＮＡウイルス等の病原体を検出・解析することが可能になる。たとえば，動物の正常細胞ＲＮＡ由来のｃＤＮＡ断片を準備しておくことにより，原因不明の疾病によって死亡した動物の材料について感染していた細菌やＲＮＡウイルスを迅速に同定することが可能となる。また，ＨＩＶのように高頻度で塩基配列の置換を起こすため，ウイルス特異的なプライマーオリゴマーを選定することが困難である病原体についても，検出・解析が容易となる。このように，本発明によるＲＮＡを検出することにより，非特定のウイルスの検査を行う際の検出および作業効率を飛躍的に向上させることができる。
さらに，本発明のプライマーミックスは，ＲＮＡをテンプレートとする場合のみならず，ＤＮＡをテンプレートとしてＤＮＡを合成する場合にも用いることができる。例えば，本発明のプライマーミックスを用いてＰＣＲ反応を実施することにより，目的とするＤＮＡを優先的に増幅させることができる。
また，塩基配列決定において，ランダムプライマーの代わりに本発明のプライマーミックスを用いることにより，目的とする遺伝子の塩基配列決定を効率的に行うことができる。
実施例
以下の実施例により本発明をさらに具体的に説明するが，本発明は実施例によって限定されることはない。
実施例１
ヒトリボゾームＲＮＡの塩基配列（ＧｅｎＢａｎｋ受託番号Ｕ１３３６９）を用いて，プライマー長さを６塩基として生成したモチーフ群の各モチーフについて，リボゾーム配列における出現頻度を求めた。頻度の低い９６種のモチーフを選択し，ＤＮＡ合成機（Ａｂａｃｕｓ、シグマジェノシスＩｎｃ）を用いて選択されたモチーフの配列を有するオリゴデオキシリボヌクレオチドを合成した。
実施例２
ヒトの細胞としてＨＥｐ−２細胞（ヒト咽喉頭がん由来細胞株，アカゲサルの細胞としてＬＬＣ細胞（アカゲザル，腎臓由来細胞株），ウシの細胞としてＭＤＢＫ細胞（ウシ，腎臓由来細胞株）を使用した。また，ＲＮＡウイルスの感染例として，ヒトの細胞についてはヒトＲＳウイルス，アカゲザルの細胞についてはヒトパラインフルエンザウイルス３型を，ウシの細胞についてはウシパラインフルエンザウイルス３型を用いた。感染した細胞各１０^６個を培養プレートより剥離し，遠心管により集めた後，５ｍｌのリン酸緩衝液に再懸濁後，遠心によって試験管底部に細胞を集め，上清を除くことにより，培養液などの混入の少ない細胞を得た。
これらの細胞からＲＮＡを抽出するため，遠心後の細胞を１ｍｌのＴｒｉｚｏｌ液（Ｌｉｆｅ Ｔｅｃｈｎｏｌｏｇｉｅｓ社）に懸濁した後，ポリトロンホモジナイザー（Ｋｉｎｅｍａｔｉｃａ社）により細胞をせん断・溶解した。続けて，各細胞のせん断・溶解液に対して０．６ｍｌのクロロホルム溶液を加えて撹拌した後に遠心分離することにより，蛋白質を有機溶媒層に，核酸を水層に分離した。各試料の水層に，０．６ｍｌのイソプロパノールを加えて撹拌後遠心（１２，０００×ｇ，１０分）をすることによりＲＮＡを沈殿させ，上清をマイクロピペットにより除去することによりＲＮＡを分離した。このＲＮＡには抽出途中に混入した塩類が含まれているため，１ｍｌの７０％エタノールを加え，ＲＮＡを再懸濁後遠心することによりＲＮＡを沈殿させ，上清を除去して逆転写の鋳型に用いるＲＮＡを調製した。
ｃＤＮＡの合成は，鋳型として精製したＲＮＡ１μｇ，プライマーとして実施例１の制限プライマーミックス各０．５ｐｍｏｌ，１００ｎｍｏｌのｄＮＴＰ（ｄＡＴＰ，ｄＧＴＰ，ｄＣＴＰまたはｄＴＴＰ，Ｓｉｇｍａ−Ａｌｄｒｉｃｈ社），０．１ＭＢｑの［α−^３２Ｐ］ｄＣＴＰ（ＨＡＳ社製造，日本アイソトープ協会販売），２００ユニットの逆転写酵素（Ｓｕｐｅｒ ＳｃｒｉｐｔＩＩ，ＩｎＶｉｔｒｏｇｅｎ社），５倍濃度の逆転写用緩衝液（Ｓｕｐｅｒ ＳｃｒｉｐｔＩＩに添付），２μｌのＲＮａｓｅインヒビター（Ｓｕｐｅｒ ＲＮａｓｅＩｎ，Ａｍｂｉｏｎ社）を混合後，超純水（ミリポア社逆浸透純水製造装置により作成）により２０μｌに容量をあわせた反応液を試料ごとに氷上で作成した。従来の方法と比較するための対照として，制限プライマーミックスの代わりに５０ｎｇのオリゴｄＴプライマーまたは５０ｐｍｏｌのランダムヘキサマーを添加した反応液も氷上で作成した。これらの反応液を４２℃に設定したヒートブロック（ＴＡＨ−１Ｂ，ＴＡＩＴＥＣ）に入れて３０分間保温することにより，ＲＮＡに相補的な一本鎖ｃＤＮＡを合成した。
これらのｃＤＮＡの分子量を可視化するため，^３２Ｐで標識されているｃＤＮＡ２μｌを０．１％ブロモフェノールブルー，５０％ホルムアミド，１０％ ホルムアルデヒドで８０℃で５分間保温後，７％ホルムアルデヒドと１×ＭＯＰＳ緩衝液（２０ｍＭ ＭＯＰＳ，５ｍＭ 酢酸ナトリウム，１ｍＭ ＥＤＴＡ）を含む１．５％アガロースゲル（ＡｇａｒｏｓｅＳ，ニッポンジーン株式会社）中を，１００Ｖ電位，１×ＭＯＰＳ緩衝液中で電気泳動した。
泳動後のゲル中のｃＤＮＡは，アガロースゲルは、ナイロン膜（ＢｙｏｄｉｎｅＢ，Ｐａｌｌ社）の上に密着後，ゲル乾燥機（アトー ＡＥ−３７０１）中で吸引しつつ乾燥された。アガロースゲル内のｃＤＮＡは、乾燥したアガロース膜上に残るかナイロン膜に移動するとともに固定された。この乾燥アガロース膜とナイロン膜を密着させたまま、放射線検出機（Ｐｈｏｐｈｏｒ Ｉｍａｇｅｒ Ｃｙｃｌｏｎｅ，Ｐａｃｋａｒｄ社）のＩＰプレートに感光することにより泳動像を可視化した。結果を図１０に示す。図中，左端より６レーン，次の６レーンおよび右端の５レーンは，それそれランダムヘキサマー，オリゴｄＴプライマーおよび制限プライマーミックスを使用して合成されたｃＤＮＡについての泳動像を示す。各６レーンまたは５レーンは，使用した細胞ごとにウシ，ヒトおよびサルの標記に分かれるが，これらは，それぞれ細胞株としてＭＤＢＫ細胞（ウシ，腎臓由来細胞株），ＨＥｐ−２細胞（ヒト咽喉頭がん由来細胞株）またはＬＬＣ細胞（アカゲザル，腎臓由来細胞株）を用いたことを示す。”ウイルス感染の有無”の欄に”＋”と記載されている場合には，ＲＮＡウイルスの感染例として，ウシ細胞の場合ウシパラインフルエンザウイルス３型，ヒト細胞の場合ヒトＲＳ−ウイルス，サル細胞の場合ヒトパラインフルエンザウイルス３型で細胞を感染したことを示す。ただし，サルのＲＮＡウイルス感染細胞を制限プライマーミックスを用いて逆転写した事例はここには示されていない。
図１０から明らかなように，オリゴｄＴプライマーと同様に，ランダムヘキサマーと比較して高分子量のｃＤＮＡが得られたことが示された。
実施例３
感染細胞としては，牛パラインフルエンザウイルスに感染したＭＤＢＫ細胞を使用し，非感染細胞にはＭＤＢＫ細胞を使用した。実施例１と同じ方法および試薬によりＲＮＡを抽出し，ｃＤＮＡを合成した。ｃＤＮＡに相補的なＤＮＡは，２６μｌのｃＤＮＡ合成反応溶液中に１０μｌの５×第２鎖合成緩衝液（２５０ｍＭ Ｔｒｉｓ（ｐＨ７．５），２５０ｍＭ ＮａＣｌ，２００ｍＭ（ＮＨ_４）_２ＳＯ_４，３５ｍＭ ＭｇＣｌ_２，２５ユニットのＥ．ｃｏｌｉ ＤＮＡポリメラーゼ（東洋紡績）および２．２５ユニットのＲＮａｓｅＨ（東洋紡績），５μｌの０．１Ｍ ＤＴＴを含む５０μｌの溶液を調整し１６℃で３時間反応させる反応によりｃＤＮＡと相補的なＲＮＡをＤＮＡに置換した。結果としてＲＮＡの塩基配列情報を持つ二重鎖のＤＮＡが感染細胞および非感染細胞のそれぞれについて合成された。
感染および非感染の細胞について得られたｃＤＮＡ全量を１００μｌの溶液中で２０ユニットの制限酵素ＤｐｎＩＩ（Ｎｅｗ Ｅｎｇｌａｎｄ ＢｉｏＬａｂ，Ｉｎｃ．）で切断したのち，シリカマトリックスを用いたリンカープライマー除去キット（ＤＮＡ Ｃｌｅａｎ ａｎｄ Ｃｏｎｃｅｎｔｒａｔｏｒ５，Ｚｙｍｏ Ｒｅｓｅａｒｃｈ，Ｉｎｃ．）によりＤｐｎＩＩ切断断片のうち４０塩基対以下の断片を除去した（ＤｐｎＩＩ切断ｃＤＮＡおよびＤｐｎＩＩ切断非感染細胞由来ｃＤＮＡ）。２０μｌの溶出液（ＤＮＡ Ｃｌｅａｎ ａｎｄ Ｃｏｎｃｅｎｔｒａｔｏｒ５に添付，ＺｙｍｏＲｅｓｅａｒｃｈＩｎｃ．）中のＤｐｎＩＩ切断ｃＤＮＡおよびＤｐｎＩＩ切断非感染細胞由来ｃＤＮＡは，それぞれ１０μｌずつを２μｌのリンカー１（リンカープライマー１（５’−ＡＧＣＡＣＴＣＴＣＣＣＡＧＣＣＴＣＴＣＡＣＣＧＡＧ−３’）およびリンカープライマー２（５’−ＧＡＴＣＣＴＣＧＧＴＧＡ−３’））１ｎｍｏｌと混合後，さらに１２μｌのライゲーション反応液（Ｐｒｏｍｅｇａ社 Ｒａｐｉｄ Ｌｉｇａｔｉｏｎ Ｋｉｔ）および１μｌのリガーゼ（Ｐｒｏｍｅｇａ社 Ｒａｐｉｄ Ｌｉｇａｔｉｏｎ Ｋｉｔ）と混合し、１４℃で３０分間保温することによりリンカー１をＤｐｎＩＩ切断端に付加した（感染細胞および非感染細胞由来リンカー付ｃＤＮＡ）。次に、５ユニットのＴａｑポリメラーゼ（Ｓｉｇｍａ−Ａｌｄｒｉｃｈ），１０×Ｔａｑポリメラーゼ緩衝液（Ｓｉｇｍａ−Ａｌｄｒｉｃｈ），０．３ｍＭのｄＮＴＰを含む１００μｌの反応溶液を７２℃で１０分保温することにより、２４塩基のリンカープライマー１の相補鎖をｃＤＮＡの３‘末端から合成させた（感染細胞および非感染細胞由来リンカー配列合成ｃＤＮＡ）。この反応液に５０ｐｍｏｌの２４塩基のリンカープライマー１を加えてＰＣＲ反応（９５℃１分・７０℃３分を２０サイクル）を実施することにより，感染細胞および非感染細胞由来のｃＤＮＡのＤｐｎＩＩ切断断片を増幅した（感染細胞および非感染細胞由来アンプリコン）。
ここで増幅された感染細胞および非感染細胞由来アンプリコンの両端にはＤｐｎＩＩで切断可能なリンカーが含まれているため，これらを１００μｌの溶液中で２０ユニットの制限酵素ＤｐｎＩＩにより切断した。切断後ｃＤＮＡと切り離されたリンカー１に相当する配列は，シリカゲルを利用した低分子ＤＮＡ除去キット（ＢＩＯ１０１ ＧｅｎｅＣｌｅａｎＩＩ，Ｂｉｏｇｅｎ社）によって除去された（ＤｐｎＩＩ切断感染細胞由来アンプリコンおよびＤｐｎＩＩ切断非感染細胞由来アンプリコン）。ＤｐｎＩＩ切断感染細胞由来アンプリコンについては，前述の方法に従い末端にリンカー２（５’−ＡＣＣＧＡＣＧＴＣＧＡＣＴＡＴＣＣＡＴＧＡＡＣＧ−３’（リンカープライマー３）および５’−ＧＡＴＣＣＧＴＴＣＡＴＧ−３’（リンカープライマー４））が、リンカー１を連結したときと同様の方法により、改めて連結された（リンカー２接続感染細胞由来アンプリコン）。
２０ｎｇのリンカー２接続感染細胞由来アンプリコンを２μｇのＤｐｎＩＩ切断非感染細胞由来アンプリコンと混合した後，１０分の１容量の３Ｍ酢酸ナトリウムおよび２．５倍量のエタノールを加えることにより、混合液中のｃＤＮＡ（アンプリコン）をすべて沈殿させた。この沈殿を乾燥後、４μｌの３×ＥＥ緩衝液（１０ｍＭ Ｎ−（２−ヒドロキシエチル）ピペラジン−Ｎ’−（３−プロパンスルホン酸），１ｍＭ ＥＤＴＡ，ｐＨ８．０）中に溶解し，９９℃で３分加熱後急冷し，ここに５Ｍ ＮａＣｌ溶液を１μｌ加えた後，６７℃で２０時間保温することにより感染細胞由来および非感染細胞由来の二種のｃＤＮＡ（アンプリコン）をハイブリダイズさせた。
ハイブリダイズ終了後、その溶液をＴＥ緩衝液（１０ｍＭ Ｔｒｉｓ，１ｍＭＥＤＴＡ，ｐＨ８，０）で２００倍希釈した後、その溶液の１μｌを５ユニットのＴａｑポリメラーゼ（Ｓｉｇｍａ−Ａｌｄｒｉｃｈ），１０×Ｔａｑポリメラーゼ緩衝液（Ｓｉｇｍａ−Ａｌｄｒｉｃｈ），０．３ｍＭのｄＮＴＰを含む１００μｌの反応溶液に加え７２℃で１０分保温することによりリンカープライマー３の反対鎖を合成した。その後に，５０ｐｍｏｌの２４塩基のリンカープライマー３（５’−ＡＣＣＧＡＣＧＴＣＧＡＣＴＡＴＣＣＡＴＧＡＡＣＧ−３’）を加えてＰＣＲ反応（９５℃１分・７０℃３分を２０サイクル）を実施することにより，ＤｐｎＩＩ切断断片を増幅した。
この反応では両端にリンカー配列を持つウイルス由来ｃＤＮＡの制限酵素切断断片のみが増幅されることが期待された。２％アガロースゲル中（１×ＴＡＥ緩衝液中）で非変性条件で泳動した結果を図１１に示す。レーン１および２：ウシパラインフルエンザウイルス３型で感染されたＭＤＢＫ細胞（図中ウイルス感染の有無”＋”で示す）および非感染ＭＤＢＫ細胞（図中ウイルス感染の有無”−”で示す）から抽出したＲＮＡ（図中”ＲＮＡ”で示す）；レーン４および５：逆転写することにより合成されたｃＤＮＡを二重鎖ｃＤＮＡとした反応産物（図中”二本鎖ｃＤＮＡｓ”で示す）；レーン６および７：二重鎖ｃＤＮＡを制限酵素ＤｐｎＩＩで切断後リンカーを接続した後，リンカー配列をプライマーとして非特異的に増幅した反応産物（図中”アンプリコン”で示す）；レーン８：ランダムプライマーを用いて作成した二重鎖ｃＤＮＡを非特異的に増幅した反応産物をＲＤＡの原法に従って処理してウイルス感染細胞に特異的なＤｐｎＩＩ切断断片を増幅した反応産物（図中”ＲＤＡ（アンプリコン（Ｉｎｆ）−アンプリコン（Ｈｏｓｔ））”で示す）；レーン９：ＲＤＡ法の一部を利用して増幅した感染細胞特異的なＲＮＡ由来のＤＮＡ断片（図中，”ＲＤＡｃＤＮＡ（Ｉｎｆ）−アンプリコン（Ｈｏｓｔ））”で示す）。
ＲＤＡ法の際の感染細胞由来リンカーつきＤＮＡ断片と非感染細胞由来ＤＮＡ断片の混合比率は，１：１００（”ＲＤＡ（アンプリコン（Ｉｎｆ）−アンプリコン（Ｈｏｓｔ））”および”ＲＤＡ（ｃＤＮＡ（Ｉｎｆ））−アンプリコン（Ｈｏｓｔ））”のうちの右側のレーン）または１：１０（”ＲＤＡ（ｃＤＮＡ（Ｉｎｆ））−Ａｍｐｌｉｃｏｎ（Ｈｏｓｔ））”のうちの左側のレーン）である。図中”Ｍ”で示す分子量マーカーは，低分子量より１００，２００，３００，４００，５００（高濃度），６００，７００，８００，９００，１０００，２０００（高濃度）および３０００（高濃度）塩基対のバンドから構成される。二重鎖ｃＤＮＡ中には，ｃＤＮＡ合成中分解されなかったリボゾームＲＮＡが含まれている。
図１１から明らかなように，この制限プライマーミックスを用いて合成されたｃＤＮＡについてのＲＤＡ法によって，感染細胞にのみ存在するウイルスＲＮＡが増幅された。
産業上の利用性
本発明にしたがって選定されたプライマーの混合物を用いることにより，ヒトや家畜の細胞中において細菌やウイルスの遺伝子を容易に検出することができる。
【図面の簡単な説明】
図１は，本実施の形態におけるプライマー選定装置のハードウエア構成例を示す。
図２は，プライマー選定装置の機能ブロック図である。
図３は，図２のプライマー選定装置によるプライマー選定処理の概要を示す。
図４は，頻度の集計結果の一例を表す表である。
図５は，頻度情報の生成および処理の例を示すフローチャートである。
図６は，プライマー選定部がプライマーを選定する処理の一例を示すフローチャートである。
図７は，プライマー選定部がプライマーを選定する処理の別の例を示すフローチャートである。
図８は，本発明の別の実施形態におけるプライマー選定装置１の動作を示すフローチャートである。
図９は，ＲＤＡ法（Ｒｅｐｒｅｓｅｎｔａｔｉｏｎａｌ Ｄｉｆｆｅｒｅｎｃｅ Ａｎａｌｙｓｉｓ）の概要を示す。
図１０は，本発明のプライマーミックスを用いてＲＮＡから逆転写した一本鎖ｃＤＮＡの電気泳動写真を示す。
図１１は，本発明のプライマーミックスを用いて製造した二本鎖ＤＮＡおよびこれを増幅した反応産物の電気泳動写真を示す。

Claims (25)

1. ＤＮＡ合成用のプライマーを選定するプライマー選定方法であって，
   プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，
   前記頻度情報に基づき，前記非標的配列に出現する頻度が所定の低頻度閾値以下のモチーフをプライマーとして選定するプライマー選定ステップと，
   を含むことを特徴とし、ここで、前記モチーフの長さは６塩基以上である、プライマー選定方法。
2. ＤＮＡ合成用のプライマーを選定するプライマー選定方法であって，
   プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度および標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，
   前記非標的配列に出現する頻度が所定の低頻度閾値以下であり，かつ前記標的配列に出現する頻度が所定の高頻度閾値以上であるモチーフをプライマーとして選定するプライマー選定ステップと，
   を含むことを特徴とし、ここで、前記モチーフの長さは６塩基以上である、プライマー選定方法。
3. ＤＮＡ合成用のプライマーを選定するプライマー選定方法であって，
   プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度および標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，
   前記頻度情報に基づき，モチーフが前記非標的配列に出現する頻度と前記標的配列に出現する頻度の相違に基づいてプライマーを選定するプライマー選定ステップと，
   を含むことを特徴とし、ここで、前記モチーフの長さは６塩基以上である、プライマー選定方法。
4. 非標的配列がリボゾームの塩基配列である，請求項１−３のいずれかに記載の方法。
5. 非標的配列がリピート配列の塩基配列である，請求項１−３のいずれかに記載の方法。
6. 標的配列が細菌，ウイルス，寄生虫または真菌の塩基配列である，請求項２または３に記載の方法。
7. ＤＮＡ合成用のプライマーを選定するプライマー選定装置であって，
   プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度を示す頻度情報を求める頻度情報生成部と，
   前記頻度情報に基づき，前記非標的配列に出現する頻度が所定の低頻度閾値以下のモチーフをプライマーとして選定するプライマー選定部と，
   を含むことを特徴とし、ここで、前記モチーフの長さは６塩基以上である、プライマー選定装置。
8. ＤＮＡ合成用のプライマーを選定するプライマー選定装置であって，
   プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度および標的配列に出現する頻度を示す頻度情報を求める頻度情報生成部と，
   前記頻度情報に基づき，前記非標的配列に出現する頻度が所定の低頻度閾値以下であり，かつ前記標的配列に出現する頻度が所定の高頻度閾値以上であるモチーフをプライマーとして選定するプライマー選定部と，
   を含むことを特徴とし、ここで、前記モチーフの長さは６塩基以上である、プライマー選定装置。
9. 非標的配列がリボゾームの塩基配列である，請求項７または８に記載の装置。
10. 非標的配列がリピート配列の塩基配列である，請求項７または８に記載の装置。
11. 標的配列が細菌，ウイルス，寄生虫または真菌の塩基配列である，請求項８記載の装置。
12. ＤＮＡ合成用のプライマーを選定するための処理をコンピュータに実行させるためのプログラムであって，
    前記プログラムは，
    プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，
    前記頻度情報に基づき，前記非標的配列に出現する頻度が所定の低頻度閾値以下のモチーフをプライマーとして選定するプライマー選定ステップと，
    を前記コンピュータに実行させることを特徴とし、ここで、前記モチーフの長さは６塩基以上である、プログラム。
13. ＤＮＡ合成用のプライマーを選定するための処理をコンピュータに実行させるためのプログラムであって，
    前記プログラムは，
    プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度および標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，
    前記頻度情報に基づき，前記非標的配列に出現する頻度が所定の低頻度閾値以下であり，かつ前記標的配列に出現する頻度が所定の高頻度閾値以上であるモチーフをプライマーとして選定するプライマー選定ステップと，
    を前記コンピュータに実行させることを特徴とし、ここで、前記モチーフの長さは６塩基以上である、プログラム。
14. 非標的配列がリボゾームの塩基配列である，請求項１２または１３に記載のプログラム。
15. 非標的配列がリピート配列の塩基配列である，請求項１２または１３に記載のプログラム。
16. 標的配列が細菌，ウイルス，寄生虫または真菌の塩基配列である，請求項１３記載のプログラム。
17. ＤＮＡ合成用のプライマーを選定するための処理をコンピュータに実行させるためのプログラムを記録したコンピュータ読取可能記録媒体であって，
    前記プログラムは，
    プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，
    前記頻度情報に基づき，前記非標的配列に出現する頻度が所定の低頻度閾値以下のモチーフをプライマーとして選定するプライマー選定ステップと，
    を前記コンピュータに実行させることを特徴とし、ここで、前記モチーフの長さは６塩基以上である、コンピュータ読取可能記録媒体。
18. ＤＮＡ合成用のプライマーを選定するための処理をコンピュータに実行させるためのプログラムを記録したコンピュータ読取可能記録媒体であって，
    プライマー候補モチーフ群の各モチーフが非標的配列に出現する頻度および標的配列に出現する頻度を示す頻度情報を求める頻度情報生成ステップと，
    前記頻度情報に基づき，前記非標的配列に出現する頻度が所定の低頻度閾値以下であり，かつ前記標的配列に出現する頻度が所定の高頻度閾値以上であるモチーフをプライマーとして選定するプライマー選定ステップと，
    を前記コンピュータに実行させることを特徴とし、ここで、前記モチーフの長さは６塩基以上である、コンピュータ読取可能記録媒体。
19. 非標的配列がリボゾームの塩基配列である，請求項１７または１８に記載の記録媒体。
20. 非標的配列がリピート配列の塩基配列である，請求項１７または１８に記載の記録媒体。
21. 標的配列が細菌，ウイルス，寄生虫または真菌の塩基配列である，請求項１８記載の記録媒体。
22. 請求項４記載の方法により、プライマー長さを６塩基としてプライマーを選定し、選定されたプライマーの配列からなるオリゴヌクレオチドの混合物を製造し、前記オリゴヌクレオチドの混合物をプライマーとして，前記非標的配列および前記標的配列の存在下で、逆転写酵素またはＤＮＡ合成酵素を用いてＤＮＡを合成することを特徴とする，ＤＮＡの製造方法。
23. 試料において正常細胞に存在しないＲＮＡを検出する方法であって，
    請求項１−６のいずれかに記載の方法によりプライマーを選定し、選定されたプライマーの配列からなるオリゴヌクレオチドの混合物を製造し、ここで非標的配列はリボゾームの塩基配列であり、
    試料からＲＮＡを抽出し，
    前記オリゴヌクレオチドの混合物をプライマーとして逆転写により前記ＲＮＡからｃＤＮＡを合成し，そして，
    該ｃＤＮＡ中に正常細胞に存在しないＲＮＡに対応する配列が存在するか否かを判定することを特徴とする方法。
24. 正常細胞に存在しないＲＮＡが，細菌，ウイルス，寄生虫または真菌に由来するＲＮＡである，請求項２３記載の方法。
25. 試料が食中毒患者の糞便である，請求項２３記載の方法。

Images