JP3924976B2 - 遺伝子の発現頻度の解析方法 - Google Patents
遺伝子の発現頻度の解析方法 Download PDFInfo
- Publication number
- JP3924976B2 JP3924976B2 JP03853899A JP3853899A JP3924976B2 JP 3924976 B2 JP3924976 B2 JP 3924976B2 JP 03853899 A JP03853899 A JP 03853899A JP 3853899 A JP3853899 A JP 3853899A JP 3924976 B2 JP3924976 B2 JP 3924976B2
- Authority
- JP
- Japan
- Prior art keywords
- restriction enzyme
- sequence
- primer
- vector
- tag
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6809—Methods for determination or identification of nucleic acids involving differential detection
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Analytical Chemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Molecular Biology (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、遺伝子の発現頻度の解析方法に関し、詳しくは、遺伝子発現の動的変化を捉えるために、細胞において蛋白質をコードする全遺伝子より発現しているmRNAの種類と量を、微量生体サンプルを用いて解析することが可能な方法に関する。
【０００２】
【従来の技術】
ゲノム上に存在する蛋白質をコードする遺伝子の総数は、ヒトで約10万種と予想されている。既に全ゲノム構造が明らかとなった酵母では、蛋白質をコードする遺伝子の総数は約5000種と推定されている。
【０００３】
近年、公共の遺伝子データバンクが欧米、日本を中心に設立され、世界中の膨大な量の遺伝子情報がデータバンクに登録されており、さらに日々新たな情報がデータバンクに集まって来る。世界規模で行なわれているヒトゲノムプロジェクト（human genome project）は、2005年を目標にヒトゲノムの全遺伝子の配列を明らかにしつつあり、その遺伝子情報もデータバンクに登録されつつある。ある遺伝子配列についてデータバンクに問い合わせると、その配列と同一又は類似の配列を持つ遺伝子が既に登録されているか否か、さらには登録されていれば、その配列に関する情報、例えば、遺伝子名、機能、関連する文献等を知ることができる。このような検索をホモロジー検索といい、ホモロジー検索を行うソフトウエアは何種類かあるが、多数の検索を行なう場合は検索時間の短いBLASTが通常用いられている。
【０００４】
通常、細胞の遺伝子は全ての遺伝子がmRNAに読みとられ、mRNAから蛋白質が生産されている訳ではなく、ヒトでは１細胞の発現している遺伝子は約15000種と推定されている。このように、細胞中では何種類ものゲノム遺伝子が発現しており、それに応じた種類のメッセンジャーRNA(以下、「mRNA」という)が生成しているが、発現している遺伝子の種類及び量(以下、「遺伝子発現頻度情報」ともいう)は、細胞の種類、状態により変動する。例えば、血液幹細胞からリンパ球前駆細胞、プレＢ細胞、Ｂ細胞、活性化Ｂ細胞と分化するに従い、各細胞で共通に発現する遺伝子もあるが、全く違う遺伝子発現を示す。
【０００５】
上記のような遺伝子発現頻度情報を測定することを、遺伝子発現プロファイル解析という。細胞の生命活動を担っているのは主として蛋白質であり、mRNAから翻訳される蛋白質の種類と量を解析することは、遺伝子発現解析として重要であるが、現状では全蛋白質のプロファイルを得ることは技術的に困難である。一方、mRNAについては、全種類のmRNAを測定することが可能となっている。
【０００６】
遺伝子発現プロファイル解析法として最初に報告されたのは、Body Map法（Gene, 174, 151-158 (1996)）である。Body Map法の概略は以下の通りである。mRNAをベクター上のポリＴ配列とmRNAの3'末端にあるポリAテールを結合させ、ベクターポリＴ配列をプライマーとしてcDNAを合成する。更に、制限酵素MboIでcDNAを切断する。MboIサイトはcDNA上で平均300塩基対に一つあるため、ベクター上のcDNAは平均300塩基対に分断される。このとき、最もポリAテールよりのcDNAは、ベクターに結合したまま残る。このcDNA断片を持つベクターを閉環させ、それを大腸菌に導入してcDNAライブラリーを作製する。ライブラリーから約1000クローンを任意に選択して、それぞれについて平均300塩基対の塩基配列を決定する。それらの配列の中から同じ配列を含むクローン毎にまとめて、それぞれの配列の種類と出現頻度を算出して遺伝子発現プロファイルを得る。各cDNA配列はデータバンクとのホモロジー検索(BLAST検索)を行い、既知の遺伝子と同一の配列を有するクローンにはその遺伝子の名称を与える。配列がデータバンクに登録されていない場合は、その配列に該当する遺伝子は存在しないものとする。
【０００７】
BLAST検索によりホモロジー検索を行なうためには、最低11塩基対の情報が必要である。10塩基からなる配列の種類は約百万種であり、人で存在すると予想される遺伝子の種類10万種を遙かに越える。すなわち、11塩基対の情報があれば、その配列を持つ遺伝子は特定することができ、遺伝子発現プロファイル解析が可能である。したがって、多量のシークエンスを必要とするBody Mapによる遺伝子発現プロファイル解析を効率化するために、Body Mapにおける約300塩基対のcDNA断片を、更に11塩基対以上の短い断片(「タグ（tag）」と呼ばれる)とし、更にこの断片同士を多数連結してベクターに挿入することにより、連結タグのライブラリーを作成し、Body Mapと同様に約1000クローンを任意に選択して、連結タグのDNA配列を決定すれば、Body Mapと同じ手数でより多くの遺伝子発現情報を得ることが期待できる。タグは遺伝子配列を代表し、タグの出現頻度はその遺伝子の発現頻度を現す。通常、１回のシークエンスで判読出来るDNA配列の長さは約600塩基対であるから、１回のシークエンスで最大50個程度のタグのDNA配列を判読出来る。すなわち、Body Map法に比べて最大約50倍の効率で、遺伝子発現プロファイル解析を行うことが可能になる。
【０００８】
上記の考えに基づいた遺伝子発現プロファイル解析法として、Serial analysis of gene expression（SAGE）法がある（米国特許第527,154号、第544,861号、欧州特許公開第0761 822号）。SAGEは、ビオチンが3'末端に結合したポリＴをプライマーとしてcDNAを作製し、Body Mapと同様にMboI等の制限酵素（「アンカリング酵素(anchoring enzyme)」と呼ばれる）でcDNAを切断した後、ビオチンが結合した3'末端を含むcDNA断片をアビジンビーズに吸着させ、ビーズを二分して、それぞれのビーズに吸着したcDNA断片（約13bp）に２種のリンカー（Ａ又はＢ）の一方づつを結合させる。各リンカーにはBsmFIのようなClass II制限酵素（「タグ化酵素(tagging enzyme)」と呼ばれる）のサイトを含ませておく。タグ化酵素でcDNA断片をビーズから切り出し、切断部位を平滑化し、Ａリンカーに接続するタグとＢリンカーに接続するタグを連結させる。これは、ダイタグ（ditag）と称される。ＡリンカーとＢリンカーを認識するプライマーを用いてPCRによりダイ夕グを増幅する。増幅されたダイタグ同士を多数連結させてベクターに組み込み、シークエンスする。１シークエンスで最大50程度のタグシークエンスを得ることができる。このタグシークエンス情報を集計して、遺伝子発現頻度を導き出す。
【０００９】
また、遺伝子の発現頻度を解析する他の方法として、遺伝子チップ法及び遺伝子マイクロアレイ法がある。いずれも、遺伝子断片を何等かの板の上(通常はスライドガラス)に、極めて高密度に（約10個/mm2以上）並べて張り付けたものが用いられる。このチップと蛍光ラベルしたmRNAをハイブリダイズさせて、mRNAの種類と量を測定する。
【００１０】
【発明が解決しようとする課題】
上記のように、遺伝子の発現頻度を解析する方法がいくつか開発され、それなりの成果は得られている。現状では、全ての真核生物の全ての遺伝子の発現頻度を測定するにはSAGE法が最も有効な手段であるが、実際にこの手法を実施しようとすると多くの問題に遭遇して、殆どの研究機関でSAGE法を再現することができなかった。すなわち、SAGE法は、手技が難しく、特別に訓練した者しか実施できない。また、測定にはmRNAが1μg程度必要であり、小量の検体しか入手できない場合、例えば臨床生検材料では測定することや、組織の微細な部分の遺伝子発現の違いを測定するこは、事実上不可能である。さらに、原理的に測定ミスが多い。
【００１１】
SAGE法では、正確にタグの配列を読み取ることは極めて重要である。なぜなら、タグは短く（約13bp）、一カ所でも読み間違いがあると同一物であるにも関わらず別物と判断したり、また、異なるタグを同一物とみなしてしまうことが起こる。ところが、SAGE法にはこの間違いが生じる可能性が高い。なぜなら、SAGE法では２つのタグをつなげてダイタグを形成するが、このときにタグとタグの境界が不明瞭になるからである。タグはBsmFIやFokIなどの制限酵素で切り取られる短い遺伝子断片である。しかし、これらの酵素は切断部位が必ずしも安定しておらず、切り出されたタグの長さはまちまちである。このように長さの異なるタグが混在している状態でタグとタグをつなげてダイタグにすると、タグとタグの連結部位の塩基は、もともとどちらのタグ由来のものなのか判らなくなる。その結果、正確なタグの配列が得られなくなる。このように、SAGE法には原理的に不可避な欠点がある。また、SAGE法にはアビジン、ビオチンビーズを用いてDNAを回収する操作があるが、アビジン、ビオチンビーズを用いてコンタミネーションを生じることなくDNAを回収することは実際には非常に難しく、プロトコール通りの操作では正確なデータを得ることは非常に困難である。また、SAGE法ではデータを得るためには多量のmRNAを必要とし、臨床サンプルなどのサンプル量に限りがあるものの場合、十分なmRNAは得られず、SAGE法は実施困難であった。
【００１２】
また、遺伝子チップ法及び遺伝子マイクロアレイ法では、Body Map法やSAGE法と違って、構造の解っている遺伝子しか測定することができない。したがって、現状では、全ての生物の全ての遺伝子の発現頻度を測定することはできない。
【００１３】
現状では、全ての真核生物の全ての遺伝子の発現頻度を測定するにはSAGE法が最も有効な手段であるが、実際にこの手法を実施しようとすると多くの問題に遭遇して、殆どの研究機関でSAGE法を再現することができなかった。SAGE法の問題点を上げると、（１）手技が難しく、特別に訓練した者しか実施出来ない。（２）測定にmRNA量が1μg程度必要であり、小量の検体しか入手できない場合、例えば臨床生検材料は測定できない。同様に、組織の微細な部分の遺伝子発現の違いを測定することができない。（３）ダイタクを測定するために、原理的に測定ミスが多い。
【００１４】
本発明は、上記現状に鑑みなされたものであり、通常の研究者が容易に実施することができ、また、ごく微量の検体から正確な遺伝子発現頻度解析が行なうことができる方法を提供することを課題とする。
【００１５】
【課題を解決するための手段】
本発明者らは、上記課題を解決するために鋭意検討を行った結果、ポリＴ配列を有するベクタープライマーを用いてmRNAからcDNAを合成し、同ベクター上でcDNA配列をタグ化し、得られたタグを、タグの末端が識別できるような配列を介在させて連結することによりコンカテマーを形成させ、このコンカテマーの塩基配列を解析することによって、効率よく、しかも高精度で遺伝子の発現頻度を解析することができることを見出し、「MAGE（Micro-analysis of Gene Expression）」と名付けた本発明を完成するに至った。
【００１６】
すなわち本発明は、
（１）以下のステップを含む遺伝子の発現頻度を解析する方法；
（ａ）直鎖状のプラスミドベクターの一方の３’末端に一本鎖のポリＴ配列と、その内側に第一の制限酵素認識配列とを有し、他方の末端近傍に第二の制限酵素認識配列と、さらにその内側にタイプIIＳ制限酵素認識配列とを有するベクタープライマーと、遺伝子の発現頻度を解析しようとする細胞に由来するｍＲＮＡとをアニールさせてｃＤＮＡ合成を行い、ｃＤＮＡが結合したベクタープライマーを生成するステップと、
（ｂ）前記ｃＤＮＡが結合したベクタープライマーを、ベクタープライマーを切断せず、かつ、前記第二の制限酵素の切断末端と同じ形の末端を生じさせる第三の制限酵素と、前記第二の制限酵素とを用いて消化してｃＤＮＡの上流側を切除し、ベクタープライマーを閉環するステップと、
（ｃ）閉環したベクタープライマーを、前記第一の制限酵素及び前記タイプIIＳ制限酵素で消化して、ｃＤＮＡの一部からなるタグを残してｃＤＮＡの下流側を切除し、ベクタープライマーを再び閉環するステップと、
（ｄ）前記ベクタープライマーを鋳型とし、ベクタープライマー中のタグの両側のそれぞれの領域に相当する塩基配列を有するオリゴヌクレオチドをプライマーとしてポリメラーゼ・チェイン・リアクション（ＰＣＲ）を行い、タグを増幅するステップと、
（ｅ）前記増幅産物を連結してタグのコンカテマーを形成するステップと、
（ｆ）前記コンカテマーの塩基配列を決定し、その塩基配列中に出現するタグの種類及びその頻度を調べるステップ。
【００１７】
（２）前記ステップ（ｅ）において、連結反応を、一方の末端がタグの末端と同じ形を有するアダプター存在下で行い、それによりコンカテマーの両端にアダプターを配置させ、アダプターの配列に相当する配列を有するオリゴヌクレオチドをプライマーに用いてＰＣＲを行うことによりコンカテマーを増幅することを特徴とする（１）の方法。
（３）前記ステップ（ｅ）の後に、前記コンカテマーの塩基配列の決定を、塩基配列決定用のクローニングベクターにクローニングして行うことを特徴とする（１）又は（２）の方法。
【００１８】
（４）前記第三の制限酵素の認識配列が４塩基である（１）〜（３）のいずれかの方法。
（５）ベクタープライマーは、第二の制限酵素認識配列の切断点の内側と同じ位置又はその内側に切断点を有し、かつ、タイプIIＳ制限酵素による切断によってベクタープライマーから切除されない第四の制限酵素認識配列を有し、
前記（ｄ）のステップで用いるプライマーのうち、タグの下流側のプライマーは、第四の制限酵素切断末端と同じ形の末端を生じさせる第五の制限酵素の認識配列を有し、増幅されたプライマーを第四の制限酵素及び第五の制限酵素で消化した後にコンカテマーを形成させる（１）〜（４）のいずれか一項に記載の方法。
【００１９】
（６）前記ベクタープライマーは、第一の制限酵素認識配列のさらに内側に第五の制限酵素認識配列と一塩基異なる塩基配列を有し、前記タグの下流側のプライマーを用いたＰＣＲにより前記一塩基異なる塩基配列が第五の制限酵素認識配列に変換される（５）の方法。
（７）前記第三の制限酵素、第四の制限酵素及び第五の制限酵素が同一である前記（６）の方法。
【００２０】
（８）ベクタープライマーが、マルチクローニング部位を有するプラスミドを同マルチクローニング部位の二箇所で切断して得られる直鎖状プラスミドと、この直鎖状プラスミドの一方の末端と同じ形の末端を有し、一本鎖のポリＴ配列を有する部分二本鎖ＤＮＡとが連結されたものである（１）〜（７）のいずれかの方法。
【００２１】
【発明の実施の形態】
以下に、本発明をさらに詳細に説明する。
本発明の方法は、上記（ａ）〜（ｆ）のステップを含む。以下、各ステップ毎に説明する。
【００２２】
＜１＞ステップ（ａ）
最初のステップでは、直鎖状のプラスミドベクターの一方の３’末端に一本鎖のポリＴ配列と、その内側に第一の制限酵素認識配列とを有し、他方の末端近傍に第二の制限酵素認識配列と、さらにその内側にタイプIIＳ制限酵素認識配列とを有するベクタープライマー（以下、「逆転写用ベクタープライマー」ともいう）を用いる。ここで、タイプIIＳ制限酵素とは、制限酵素により認識される配列から離れた特定の部位を開裂する制限酵素をいう。
【００２３】
ベクタープライマーの構造の一例を、図２Ｄに示す。このベクタープライマーは、例えば、マルチクローニング部位を有するプラスミドを同マルチクローニング部位の二箇所で切断して得られる直鎖状プラスミドと、この直鎖状プラスミドの一方の末端と同じ形の末端を有し、一本鎖のポリＴ配列を有する部分二本鎖ＤＮＡとを連結することにより、調製することができる。図１に、マルチクローニング部位を有するプラスミドの例として、MAGE/pUC19の構造を示す。このプライマーは、基本構造として公知のクローニングベクターであるpUC19を含み、そのマルチクローニングサイトのEcoRlとHindlll制限酵素部位の間に、種々の制限酵素認識部位を含む配列（Z fragmentと呼ぶ。配列を配列番号１に示す。）が挿入されたものである。
尚、以下の具体例では、MAGE/pUC19は、Damメチレーション系を有するエシェリヒア・コリ宿主、例えばJM109株等を用いて調製したものを使用する例について説明するが、本発明に必須ではない。また、基本構造となるベクターは、pUC19に限られるものでなく、それ以外のベクター、例えばpUC18、pBR322、pHSG299、pHSG298、pHSG399、pHSG398、RSF1010、pMW119、pMW118、pMW219、pMW218等、種々のベクターを用いることができる。尚、用いるベクターの宿主は、形質転換、宿主からのベクターの回収等、通常の遺伝子組換え技術を適用することができるものであれば特に制限されるものではないが、通常、エシェリヒア・コリが用いられる。
【００２４】
以下、マルチクローニング部位を有するプラスミドとしてMAGE/pUC19を用いた例について説明するが、本発明はこれに限られるものではない。
まず、MAGE/pUC19ベクター（図２Ａ）を、制限酵素BstXI及びPstIで処理し、切り取られた小断片を除去する（図２Ｂ）。次にBstXIによる切断末端に、一本鎖のポリＴ配列を有する部分二本鎖ＤＮＡ（図２Ｃに示すポリＴアダプター）をライゲーションする。これにより、MAGE/pUC19ベクターの一端から一本鎖のポリＴが飛び出した形のベクタープライマー（又は、逆転写用ベクタープライマー）が形成される（図２Ｄ）。このベクタープライマーでは、第一の制限酵素はPmeIであり、第二の制限酵素はBglIIであり、タイプIIＳ制限酵素はBsgIである。
【００２５】
第二の制限酵素BglIIの認識配列にはMboI認識配列が重複して含まれており、BglII末端をMboI末端と連結すると、MboIで切断することができる。また、この例では、ベクタープライマーの第一の制限酵素（PmeI）認識配列のさらに内側に、第四の制限酵素（MboI）認識配列と一塩基異なる塩基配列を含むように設計されている（図３ＥのΔMboI）。このΔMboI及びベクタープライマーの第一の制限酵素認識配列は、ポリＴアダプターに由来する。
【００２６】
第一の制限酵素及び第二の制限酵素は、ベクタープライマーを一箇所で切断するものであれば特に限定されない。また、タイプIIＳ制限酵素及びその存在位置は、ベクタープライマーから第四の制限酵素部位を切除せず、かつ、cDNAの上流側の一部をベクタープライマーに残したまま切断するものであれば特に制限されない。具体的には、例えば上記BsgIの他に、BsmFI等が挙げられる。
また、ポリＴ配列の長さは、mRNAのポリＡ配列とアニールできる程度の長さであればよく、通常10〜50塩基程度である。
【００２７】
上記のようなベクタープライマーと、遺伝子の発現頻度を解析しようとする細胞に由来するｍＲＮＡとをアニールさせる。その状態で逆転写反応を行うと、ポリＴがプライマーとなってcDNA合成が始まる（図３Ｅ）。そして、合成された一鎖目のcDNAを鋳型として二鎖目を合成し、二本鎖cDNAが合成できる（図３Ｆ）。多数のmRNA分子から多数のｃＤＮＡが結合したベクタープライマー分子（cDNA-MAGE/pUC19）ができるが、図３Fはその典型例のひとつを示したものである。
【００２８】
mRNAは、遺伝子の発現頻度を解析しようとする細胞から抽出する。発現頻度を解析しようとする細胞は、動物や植物の組織の細胞、酵母等の微生物の細胞等、細胞中のmRNAの３’末端にポリＡ構造を有するものであれば特に制限なく使用することができる。また、原核生物のmRNAは３’末端にポリＡ（poly(A)）構造を持たないために、そのままではベクタープライマーのポリＴにアニールさせることはできないが、mRNAに酵素的にポリＡ構造を付加することにより、真核生物のmRNAと同様にして本発明の方法を実施することができる。
mRNAの調製、及びｃＤＮＡ合成、オリゴヌクレオチドの合成、制限酵素反応、ライゲーション反応、ポリメラーゼ・チェイン・リアクション（ＰＣＲ）、形質転換等の操作は、通常のｃＤＮＡクローニングに用いられるmRNAの調製法（例えばSambrook, J., Fritsch, E. F., and Maniatis, T., "Molecular Cloning A Laboratory Manual, Second Edition", Cold Spring Harbor Laboratory Press,(1989)等参照）と同様にして行うことができる。
【００２９】
＜２＞ステップ（ｂ）
次に、前記のようにして得られるｃＤＮＡが結合したベクタープライマーを、ベクタープライマーを切断せず、かつ、前記第二の制限酵素の切断末端と同じ形の末端を生じさせる第三の制限酵素と、前記第二の制限酵素とを用いて消化してｃＤＮＡの上流側を切除し、ベクタープライマーを閉環させる。閉環されたベクタープライマーは、それを適当な宿主に導入し、得られる形質転換体を培養してプラスミドを回収することにより、増幅させることができる。尚、以降の工程で得られる閉環状のベクタープライマーも同様であるが、回収されたベクタープライマーを切断しようとする制限酵素部位が、Damメチレーションにより切断されない場合には、修飾系を有しない宿主を用いる。
【００３０】
前記第三の制限酵素としては４塩基認識の酵素が好ましい。６塩基認識の酵素であると、cDNA配列中に制限酵素部位が存在しない場合がある。また、ベクタープライマーに残されるcDNAの配列が長いと、後の操作で得られるタグがmRNAのポリＡ配列から遠くなる。そのような場合、データベース中の遺伝子の発現情報（EST:expressed sequence tag）は、通常mRNAの３’末端側の一部であるため、タグ配列を検索してもヒットしない場合がある。本発明に好適な４塩基認識の制限酵素としては、MboI、TaiI等が挙げられる。
以下、第三の制限酵素として、MboIを用いた例を説明する。まず、ｃＤＮＡが結合したベクタープライマー（cDNA-MAGE/pUC19）を制限酵素BglII及びMboIで消化する（図３Ｇ）。このとき、MAGE/pUC19はDamメチレーションを受けているので、MboIでは切断されないが、逆転写により新たに合成されたcDNA中のMbolサイトだけが切断される。図３Ｇは、仮にcDNAが３箇所のMbolサイトを含む例を示している。cDNAに注目すると、ポリＡテールから上流に向かって最初のMbolサイトまでの下流側部分がMAGE/pUC19ベクタープライマーにつながったまま残り、他の部分、すなわちｃＤＮＡの上流側はMbol切断により分離除去される。また、MAGE/pUC19は、BglIIにより一箇所のみが切断される。MboI及びBglIIで切断された末端は同じ形をしており、ライゲーション反応で互いをつなぐことができる。そこで、これらの末端を介して、cDNA-MAGE/pUC19をセルフライゲーション反応により閉環する（図３Ｈ）。
【００３１】
＜３＞ステップ（ｃ）
ステップ（ｂ）で閉環したベクタープライマーを、前記第一の制限酵素及び前記タイプIIＳ制限酵素で消化して、ｃＤＮＡの一部からなるタグを残してｃＤＮＡの下流側を切除し、ベクタープライマーを再び閉環する。
【００３２】
具体的には、閉環されたcDNA-MAGE/pUC19を制限酵素BsgI及びPmeIで消化する（図４Ｉ）。これらの制限酵素の切断により、 cDNAの前記３’末端側部分のうちポリＡテールから最も遠い約１３塩基のみがベクタープライマーに残る。このｃＤＮＡの約１３塩基からなる配列を、タグ（図４Ｊ中の「tag」で示される）と呼ぶ。
【００３３】
上記のBsgIによる切断点は、５’側が２ｂｐ突出しているので、例えばT4 DNAポリメラーゼ処理により平滑化する。また、PmeIによる他方の切断末端は、平滑末端であるので、これらの末端をライゲーションすることができる（図４Ｊ）。そこで、ベクタープライマーをセルフライゲーション反応で閉環すると、タグ以外の部分のcDNAが切り取られ、MAGE/pUC19に短いタグがつながった構造のベクタープライマーが得られる（図４Ｊ、Ｋ）。
【００３４】
＜４＞ステップ（ｄ）
次に、ステップ（ｃ）で閉環されたベクタープライマーを鋳型とし、ベクタープライマー中のタグの両側のそれぞれの領域に相当する塩基配列を有するオリゴヌクレオチドをプライマーとしてＰＣＲを行い、タグを増幅する。
【００３５】
具体的には、タグを含むベクタープライマーを鋳型とし、タグの両側のベクター部分の配列に相当するプライマー（例えば配列番号２及び３）を用いてPCR反応を行う（図４Ｋ）。その結果、中央にタグを含むDNAフラグメント（この場合はプライマーの配列及びタグを含めて110bp）が増幅される（図４Ｌ）。
【００３６】
＜５＞ステップ（ｅ）
前記のＰＣＲ増幅産物を連結してタグのコンカテマーを形成させる。その際、前記（ｄ）のステップで用いるプライマーのうち、タグの下流側のプライマーに、第四の制限酵素切断末端と同じ形の末端を生じさせる第五の制限酵素の認識配列を含ませておき、増幅されたプライマーを第四の制限酵素及び第五の制限酵素で消化すると、効率よくコンカテマーを形成させることができる。第五の制限酵素の種類は特に制限されないが、第四の制限酵素と同じ制限酵素を用いると、操作が簡便となる。また、第四の制限酵素として、第三の制限酵素と同じ制限酵素を用いてもよい。さらに、第三、第四及び第五の制限酵素が同じであってもよい。ただし、第五の制限酵素と第四の制限酵素に同じものを用いる場合は、ステップ（ｂ）でcDNAの全配列がベクタープライマーから切除されてしまうので、それを防ぐために、ベクタープライマーの第一の制限酵素認識配列のさらに内側に、第五の制限酵素認識配列と一塩基異なる塩基配列を含ませておき、タグの下流側に第五の制限酵素認識配列を有するプライマーを用いてＰＣＲを行うことにより、一塩基異なる塩基配列を第五の制限酵素認識配列に変換するようにする必要がある。また、第三の制限酵素と第四の制限酵素が同じである場合には、ベクタープライマーの第四の制限酵素認識部位は、メチレーションにより切断されないようにしておく必要がある。その場合には、第四の制限酵素で消化する前に、ステップ（ｃ）で閉環されたベクタープライマーを修飾系を持たない宿主に導入して脱修飾する。
【００３７】
タグの連結反応を、一方の末端がタグの末端と同じ形を有するアダプター存在下で行い、それによりコンカテマーの両端にアダプターを配置させると、アダプターの配列に相当する配列を有するオリゴヌクレオチドをプライマーに用いてＰＣＲを行うことにより、コンカテマーを増幅することができる。アダプターは、タグの量に比べて少量用いるか、一方の末端がタグの末端と異なる形を有するアダプターを用いると、多数のタグからなるコンカテマーの末端にアダプターを配置させることができる。このようにすることで、ごく微量のコンカテマーでもクローニングすることができ、結果としてごく微量のサンプルmRNAを用いた解析が可能となる。
タグとアダプターの比率は、通常モル比でタグ：アダプター＝１：１〜１：０．０１、好ましくは１：０．２〜１：０．０５である。この範囲でアダプターを用いると、タグが２〜５０程度連結したコンカテマーが得られる。
【００３８】
得られるコンカテマーは、塩基配列決定用のクローニングベクターに挿入すると、シークエンス操作が容易になる。
MAGE/pUC19を用いる場合には、ステップ（ｄ）の増幅産物を制限酵素MboIで消化することにより、タグを含む小断片が得られる（図５Ｍ）。ここで得られたタグ断片は、既知のDNA配列、すなわちGATCとAAACGに挟まれており、どの部分がタグかを明確に識別することができる。このタグ配列は両端にMbolサイトが露出しているので、ライゲーション反応によりタグ同士を連結し、コンカテマーを形成させることができる（図５Ｍ）。
【００３９】
＜６＞ステップ（ｆ）
上記のようにして得られるコンカテマーの塩基配列を決定すると、その中には多数のcDNA分子に由来するタグが含まれているので、その塩基配列中に出現するタグの種類及びその頻度を調べることにより、cDNAが由来する遺伝子の発現頻度を解析することができる。タグの種類は、既知のmRNAの部分配列（EST）の情報に関するデータベースを検索することにより調べることができる。
【００４０】
尚、タグのコンカテマーを増幅することにより、同一配列を持つフラグメントが複数個生じることになる。しかし、クローニング後のシークエンス解析において、先にシークエンスしたものと同一物とみなされるシークエンスは出現した場合は、それを解析から除くことにより、遺伝子発現頻度解析にPCRが及ぼす悪影響を排除することができる。
また、第四と第五の制限酵素として同じ酵素を用いた場合は、コンカテマー中のタグの向きは不定であるので、配列解析においては一方の鎖だけでなく逆鎖の配列も考慮する。
【００４１】
以上説明したように、本発明の方法では、SAGE法のようにダイタグを作らず、個々のタグを既知のDNA断片で挟むことにより、タグが他のタグと連続しないようにすることができる。その結果、タグの境界が曖昧になるという問題点が解決される。また、ＰＣＲを複数回行い、試料を増幅することができるので、微量のmRNAからでも解析が可能となる。さらに、ベクタープライマーを用いることにより、ベクターに融合した形のcDNA合成を行うことで、アビジンやビオチンビーズを用いずに解析することができる。
【００４２】
本発明の遺伝子発現解析法は、生命科学の研究にとって有用であるが、特に、健康人と病気の人における特定臓器、又は細胞における遺伝子発現の差異を解析する事により、病気の治療法開発や診断法開発に有用である。例えば、健康人の肝臓と肝炎の肝臓の発現遺伝子の違いを本法により解析することにより、肝炎で特異的に発現が増加したり減少したりする遺伝子を見つけることができる。これら遺伝子の肝臓での役割を調べることにより、肝炎治療を目的として、この遺伝子の働きを阻害したり、促進させたりする薬剤を開発することができる。又、遺伝子そのもの、遺伝子構造から設計したアンチセンスオリゴヌクレオチド、あるいは遺伝子を発現させて得た蛋白質等を、肝炎治療に用いることができる。
【００４３】
また、発症のメカニズムが解っていない病気に対しても、本発明の方法を用いれば、治療法を開発できる可能性がある。さらに、本発明の方法により疾患特異的に発現が変動する遺伝子を発見することができれば、治療ばかりでなく、疾患の診断法を開発することも可能となる。
【００４４】
さらに、医療における利用ばかりでなく、全ての真核生物を対象に有用遺伝子を発見する手段ともなる。例えば、ビール生産に適した酵母を変異により育種した場合、親株と変異株の遺伝子発現の変化を本発明の方法で解析し、変異により発現が変化した遺伝子を同定することができる。このようにして得たビール生産にとって有利な遺伝子を総合的に操作することにより、よりよいビール生産酵母を創製することができる。
【００４５】
また、例えば、大腸菌やコリネバクテリウムのようなアミノ酸生産菌の解析を行うことにより、より優れたアミノ酸生産菌を創製することができる。
【００４６】
【実施例】
マウスC57BL/6肝臓 1gから、Invitrogen社製FastTrack 2.0キットを用いてmRNA 35μgを抽出した。
得られたmRNA 0.97μgを用いて、上記と同様にしてMAGE/pUC19を用いて逆転写cDNA合成後、遺伝子発現頻度解析を行った。
【００４７】
MAGE/pUC19ベクター（図２Ａ）を、制限酵素BstXI及びPstIで処理し、切り取られた小断片を除去した（図２Ｂ）。次にBstXIによる切断末端に、一本鎖のポリＴ配列を有する部分二本鎖ＤＮＡ（図２Ｃのに示すポリＴアダプター）をライゲーションした（図２Ｄ）。このベクタープライマーでは、第一の制限酵素はPmeIであり、第二の制限酵素はBglIIであり、タイプIIＳ制限酵素はBsgIである。
上記ベクタープライマー0.2μgと、上記のマウス肝臓由来のmRNA 0.97μgとをアニールさせ、cDNA合成を行い、続いて合成された一鎖目のcDNAを鋳型として二鎖目を合成した（図３Ｆ）。そして、得られたcDNAの40分の１に当たる量（mRNA 0.025μg相当）を材料に遺伝子発現頻度解析を行った。
尚、別の実験では、4.5mgの微量臨床サンプルから抽出した、0.05μgの微量mRNAを実験材料にして、逆転写cDNA合成後、以下と同様にして遺伝子発現頻度解析を実施することに成功している。
【００４８】
次に、前記のようにして得られるcDNAが結合したベクタープライマーを、第二の制限酵素BglII及び第三の制限酵素であるMboIで消化した（図３Ｇ）。次に、cDNA-MAGE/pUC19のMboI末端及びBglII末端をセルフライゲーション反応により連結して閉環させた（図３Ｈ）。
【００４９】
前記の閉環したベクタープライマーをBsgI及びPmeIで消化し、T4 DNAポリメラーゼ処理により平滑化した後、ベクタープライマーを再び閉環させた（図４Ｉ〜Ｋ）。この閉環されたベクタープライマーを鋳型とし、ベクタープライマー中のタグの両側のそれぞれの領域に相当する塩基配列を有するオリゴヌクレオチド（配列番号２及び３）をプライマーとして、酵素AmpliTaq Gold（PE Biosystems社）を用いてＰＣＲを行い、タグを増幅した（図４Ｌ）。ＰＣＲ反応は、変性（95℃、0.3分）、アニーリング、およびポリメラーゼによる伸長反応（72℃、1.5分）からなる反応を65サイクル行うことにより行った。その結果、中央にタグを含む110bpのDNAフラグメントが増幅された。
【００５０】
前記のＰＣＲ増幅産物をMboIで消化してタグを切り出し、タグ：アダプターが８：１となるように、配列番号４、５に示す塩基配列を有するオリゴヌクレオチドをアニールさせたアダプターを加えて連結反応を行い、タグのコンカテマーを形成させた（図５Ｍ）。このコンカテマーを、配列番号６に示す塩基配列を有するプライマーを用いたＰＣＲにより増幅した。ＰＣＲ反応は、変性（95℃、0.3分）、アニーリング（40℃、3分）、ポリメラーゼによる伸長反応（72℃、1分）からなる反応を５サイクル行い、引き続き変性（95℃、0.3分）、アニーリング及び伸長反応（72℃、1分）からなる反応を60サイクル行った。増幅産物を制限酵素NotIで消化し、シークエンス用クローニングベクターpKF3のNotI部位に挿入して、塩基配列を決定した。
【００５１】
得られたシークエンスデータ（11000個）から、同一の配列を有するタグをグループ化し、出現頻度順に並べた。頻度上位１０番までを表１に示す。この結果から、マウス肝臓組織中で最も発現頻度が高い遺伝子はウリナリ・プロテインＩ／II遺伝子であることがわかる。以下、発現頻度の高い順に、アルブミン遺伝子、ウリナリ・プロテインIII遺伝子、アルギニノコハク酸合成酵素遺伝子と並んでいる。このように、本発明の方法は、mRNAサンプルが微量であっても、ベクタープライマーを用いることによりcDNAの回収効率が高く、またタグ配列の不明瞭さを除くことが可能であるので、遺伝子発現頻度を高精度で解析することができる。
【００５２】
【表１】
【００５３】
【発明の効果】
本発明により、遺伝子の発現頻度を、簡便に、確実に、かつ精度よく解析することができる。
【００５４】
【配列表】
【００５５】
【００５６】
【００５７】
【００５８】
【００５９】
【００６０】
【００６１】
【００６２】
【００６３】
【００６４】
【００６５】
【００６６】
【００６７】
【００６８】
【００６９】
【００７０】
【図面の簡単な説明】
【図１】 ベク夕ープライマーを作るためのプラスミドDNAの構造の一例（MAGE/pUC19）を示す図。
【図２】 ベクタープライマーの構造の一例とその構築過程を示す図。
【図３】 図３〜図５は、本発明の方法の各ステップを示す。図３は、ステップ（ａ）、（ｂ）を模式的に表す図。
【図４】 本発明の方法のステップ（ｃ）、（ｄ）を模式的に表す図。
【図５】 本発明の方法のステップ（ｅ）と、ステップ（ｅ）により得られる増幅産物をシークエンス用クローニングベクターに挿入するステップを模式的に表す図。
【発明の属する技術分野】
本発明は、遺伝子の発現頻度の解析方法に関し、詳しくは、遺伝子発現の動的変化を捉えるために、細胞において蛋白質をコードする全遺伝子より発現しているmRNAの種類と量を、微量生体サンプルを用いて解析することが可能な方法に関する。
【０００２】
【従来の技術】
ゲノム上に存在する蛋白質をコードする遺伝子の総数は、ヒトで約10万種と予想されている。既に全ゲノム構造が明らかとなった酵母では、蛋白質をコードする遺伝子の総数は約5000種と推定されている。
【０００３】
近年、公共の遺伝子データバンクが欧米、日本を中心に設立され、世界中の膨大な量の遺伝子情報がデータバンクに登録されており、さらに日々新たな情報がデータバンクに集まって来る。世界規模で行なわれているヒトゲノムプロジェクト（human genome project）は、2005年を目標にヒトゲノムの全遺伝子の配列を明らかにしつつあり、その遺伝子情報もデータバンクに登録されつつある。ある遺伝子配列についてデータバンクに問い合わせると、その配列と同一又は類似の配列を持つ遺伝子が既に登録されているか否か、さらには登録されていれば、その配列に関する情報、例えば、遺伝子名、機能、関連する文献等を知ることができる。このような検索をホモロジー検索といい、ホモロジー検索を行うソフトウエアは何種類かあるが、多数の検索を行なう場合は検索時間の短いBLASTが通常用いられている。
【０００４】
通常、細胞の遺伝子は全ての遺伝子がmRNAに読みとられ、mRNAから蛋白質が生産されている訳ではなく、ヒトでは１細胞の発現している遺伝子は約15000種と推定されている。このように、細胞中では何種類ものゲノム遺伝子が発現しており、それに応じた種類のメッセンジャーRNA(以下、「mRNA」という)が生成しているが、発現している遺伝子の種類及び量(以下、「遺伝子発現頻度情報」ともいう)は、細胞の種類、状態により変動する。例えば、血液幹細胞からリンパ球前駆細胞、プレＢ細胞、Ｂ細胞、活性化Ｂ細胞と分化するに従い、各細胞で共通に発現する遺伝子もあるが、全く違う遺伝子発現を示す。
【０００５】
上記のような遺伝子発現頻度情報を測定することを、遺伝子発現プロファイル解析という。細胞の生命活動を担っているのは主として蛋白質であり、mRNAから翻訳される蛋白質の種類と量を解析することは、遺伝子発現解析として重要であるが、現状では全蛋白質のプロファイルを得ることは技術的に困難である。一方、mRNAについては、全種類のmRNAを測定することが可能となっている。
【０００６】
遺伝子発現プロファイル解析法として最初に報告されたのは、Body Map法（Gene, 174, 151-158 (1996)）である。Body Map法の概略は以下の通りである。mRNAをベクター上のポリＴ配列とmRNAの3'末端にあるポリAテールを結合させ、ベクターポリＴ配列をプライマーとしてcDNAを合成する。更に、制限酵素MboIでcDNAを切断する。MboIサイトはcDNA上で平均300塩基対に一つあるため、ベクター上のcDNAは平均300塩基対に分断される。このとき、最もポリAテールよりのcDNAは、ベクターに結合したまま残る。このcDNA断片を持つベクターを閉環させ、それを大腸菌に導入してcDNAライブラリーを作製する。ライブラリーから約1000クローンを任意に選択して、それぞれについて平均300塩基対の塩基配列を決定する。それらの配列の中から同じ配列を含むクローン毎にまとめて、それぞれの配列の種類と出現頻度を算出して遺伝子発現プロファイルを得る。各cDNA配列はデータバンクとのホモロジー検索(BLAST検索)を行い、既知の遺伝子と同一の配列を有するクローンにはその遺伝子の名称を与える。配列がデータバンクに登録されていない場合は、その配列に該当する遺伝子は存在しないものとする。
【０００７】
BLAST検索によりホモロジー検索を行なうためには、最低11塩基対の情報が必要である。10塩基からなる配列の種類は約百万種であり、人で存在すると予想される遺伝子の種類10万種を遙かに越える。すなわち、11塩基対の情報があれば、その配列を持つ遺伝子は特定することができ、遺伝子発現プロファイル解析が可能である。したがって、多量のシークエンスを必要とするBody Mapによる遺伝子発現プロファイル解析を効率化するために、Body Mapにおける約300塩基対のcDNA断片を、更に11塩基対以上の短い断片(「タグ（tag）」と呼ばれる)とし、更にこの断片同士を多数連結してベクターに挿入することにより、連結タグのライブラリーを作成し、Body Mapと同様に約1000クローンを任意に選択して、連結タグのDNA配列を決定すれば、Body Mapと同じ手数でより多くの遺伝子発現情報を得ることが期待できる。タグは遺伝子配列を代表し、タグの出現頻度はその遺伝子の発現頻度を現す。通常、１回のシークエンスで判読出来るDNA配列の長さは約600塩基対であるから、１回のシークエンスで最大50個程度のタグのDNA配列を判読出来る。すなわち、Body Map法に比べて最大約50倍の効率で、遺伝子発現プロファイル解析を行うことが可能になる。
【０００８】
上記の考えに基づいた遺伝子発現プロファイル解析法として、Serial analysis of gene expression（SAGE）法がある（米国特許第527,154号、第544,861号、欧州特許公開第0761 822号）。SAGEは、ビオチンが3'末端に結合したポリＴをプライマーとしてcDNAを作製し、Body Mapと同様にMboI等の制限酵素（「アンカリング酵素(anchoring enzyme)」と呼ばれる）でcDNAを切断した後、ビオチンが結合した3'末端を含むcDNA断片をアビジンビーズに吸着させ、ビーズを二分して、それぞれのビーズに吸着したcDNA断片（約13bp）に２種のリンカー（Ａ又はＢ）の一方づつを結合させる。各リンカーにはBsmFIのようなClass II制限酵素（「タグ化酵素(tagging enzyme)」と呼ばれる）のサイトを含ませておく。タグ化酵素でcDNA断片をビーズから切り出し、切断部位を平滑化し、Ａリンカーに接続するタグとＢリンカーに接続するタグを連結させる。これは、ダイタグ（ditag）と称される。ＡリンカーとＢリンカーを認識するプライマーを用いてPCRによりダイ夕グを増幅する。増幅されたダイタグ同士を多数連結させてベクターに組み込み、シークエンスする。１シークエンスで最大50程度のタグシークエンスを得ることができる。このタグシークエンス情報を集計して、遺伝子発現頻度を導き出す。
【０００９】
また、遺伝子の発現頻度を解析する他の方法として、遺伝子チップ法及び遺伝子マイクロアレイ法がある。いずれも、遺伝子断片を何等かの板の上(通常はスライドガラス)に、極めて高密度に（約10個/mm2以上）並べて張り付けたものが用いられる。このチップと蛍光ラベルしたmRNAをハイブリダイズさせて、mRNAの種類と量を測定する。
【００１０】
【発明が解決しようとする課題】
上記のように、遺伝子の発現頻度を解析する方法がいくつか開発され、それなりの成果は得られている。現状では、全ての真核生物の全ての遺伝子の発現頻度を測定するにはSAGE法が最も有効な手段であるが、実際にこの手法を実施しようとすると多くの問題に遭遇して、殆どの研究機関でSAGE法を再現することができなかった。すなわち、SAGE法は、手技が難しく、特別に訓練した者しか実施できない。また、測定にはmRNAが1μg程度必要であり、小量の検体しか入手できない場合、例えば臨床生検材料では測定することや、組織の微細な部分の遺伝子発現の違いを測定するこは、事実上不可能である。さらに、原理的に測定ミスが多い。
【００１１】
SAGE法では、正確にタグの配列を読み取ることは極めて重要である。なぜなら、タグは短く（約13bp）、一カ所でも読み間違いがあると同一物であるにも関わらず別物と判断したり、また、異なるタグを同一物とみなしてしまうことが起こる。ところが、SAGE法にはこの間違いが生じる可能性が高い。なぜなら、SAGE法では２つのタグをつなげてダイタグを形成するが、このときにタグとタグの境界が不明瞭になるからである。タグはBsmFIやFokIなどの制限酵素で切り取られる短い遺伝子断片である。しかし、これらの酵素は切断部位が必ずしも安定しておらず、切り出されたタグの長さはまちまちである。このように長さの異なるタグが混在している状態でタグとタグをつなげてダイタグにすると、タグとタグの連結部位の塩基は、もともとどちらのタグ由来のものなのか判らなくなる。その結果、正確なタグの配列が得られなくなる。このように、SAGE法には原理的に不可避な欠点がある。また、SAGE法にはアビジン、ビオチンビーズを用いてDNAを回収する操作があるが、アビジン、ビオチンビーズを用いてコンタミネーションを生じることなくDNAを回収することは実際には非常に難しく、プロトコール通りの操作では正確なデータを得ることは非常に困難である。また、SAGE法ではデータを得るためには多量のmRNAを必要とし、臨床サンプルなどのサンプル量に限りがあるものの場合、十分なmRNAは得られず、SAGE法は実施困難であった。
【００１２】
また、遺伝子チップ法及び遺伝子マイクロアレイ法では、Body Map法やSAGE法と違って、構造の解っている遺伝子しか測定することができない。したがって、現状では、全ての生物の全ての遺伝子の発現頻度を測定することはできない。
【００１３】
現状では、全ての真核生物の全ての遺伝子の発現頻度を測定するにはSAGE法が最も有効な手段であるが、実際にこの手法を実施しようとすると多くの問題に遭遇して、殆どの研究機関でSAGE法を再現することができなかった。SAGE法の問題点を上げると、（１）手技が難しく、特別に訓練した者しか実施出来ない。（２）測定にmRNA量が1μg程度必要であり、小量の検体しか入手できない場合、例えば臨床生検材料は測定できない。同様に、組織の微細な部分の遺伝子発現の違いを測定することができない。（３）ダイタクを測定するために、原理的に測定ミスが多い。
【００１４】
本発明は、上記現状に鑑みなされたものであり、通常の研究者が容易に実施することができ、また、ごく微量の検体から正確な遺伝子発現頻度解析が行なうことができる方法を提供することを課題とする。
【００１５】
【課題を解決するための手段】
本発明者らは、上記課題を解決するために鋭意検討を行った結果、ポリＴ配列を有するベクタープライマーを用いてmRNAからcDNAを合成し、同ベクター上でcDNA配列をタグ化し、得られたタグを、タグの末端が識別できるような配列を介在させて連結することによりコンカテマーを形成させ、このコンカテマーの塩基配列を解析することによって、効率よく、しかも高精度で遺伝子の発現頻度を解析することができることを見出し、「MAGE（Micro-analysis of Gene Expression）」と名付けた本発明を完成するに至った。
【００１６】
すなわち本発明は、
（１）以下のステップを含む遺伝子の発現頻度を解析する方法；
（ａ）直鎖状のプラスミドベクターの一方の３’末端に一本鎖のポリＴ配列と、その内側に第一の制限酵素認識配列とを有し、他方の末端近傍に第二の制限酵素認識配列と、さらにその内側にタイプIIＳ制限酵素認識配列とを有するベクタープライマーと、遺伝子の発現頻度を解析しようとする細胞に由来するｍＲＮＡとをアニールさせてｃＤＮＡ合成を行い、ｃＤＮＡが結合したベクタープライマーを生成するステップと、
（ｂ）前記ｃＤＮＡが結合したベクタープライマーを、ベクタープライマーを切断せず、かつ、前記第二の制限酵素の切断末端と同じ形の末端を生じさせる第三の制限酵素と、前記第二の制限酵素とを用いて消化してｃＤＮＡの上流側を切除し、ベクタープライマーを閉環するステップと、
（ｃ）閉環したベクタープライマーを、前記第一の制限酵素及び前記タイプIIＳ制限酵素で消化して、ｃＤＮＡの一部からなるタグを残してｃＤＮＡの下流側を切除し、ベクタープライマーを再び閉環するステップと、
（ｄ）前記ベクタープライマーを鋳型とし、ベクタープライマー中のタグの両側のそれぞれの領域に相当する塩基配列を有するオリゴヌクレオチドをプライマーとしてポリメラーゼ・チェイン・リアクション（ＰＣＲ）を行い、タグを増幅するステップと、
（ｅ）前記増幅産物を連結してタグのコンカテマーを形成するステップと、
（ｆ）前記コンカテマーの塩基配列を決定し、その塩基配列中に出現するタグの種類及びその頻度を調べるステップ。
【００１７】
（２）前記ステップ（ｅ）において、連結反応を、一方の末端がタグの末端と同じ形を有するアダプター存在下で行い、それによりコンカテマーの両端にアダプターを配置させ、アダプターの配列に相当する配列を有するオリゴヌクレオチドをプライマーに用いてＰＣＲを行うことによりコンカテマーを増幅することを特徴とする（１）の方法。
（３）前記ステップ（ｅ）の後に、前記コンカテマーの塩基配列の決定を、塩基配列決定用のクローニングベクターにクローニングして行うことを特徴とする（１）又は（２）の方法。
【００１８】
（４）前記第三の制限酵素の認識配列が４塩基である（１）〜（３）のいずれかの方法。
（５）ベクタープライマーは、第二の制限酵素認識配列の切断点の内側と同じ位置又はその内側に切断点を有し、かつ、タイプIIＳ制限酵素による切断によってベクタープライマーから切除されない第四の制限酵素認識配列を有し、
前記（ｄ）のステップで用いるプライマーのうち、タグの下流側のプライマーは、第四の制限酵素切断末端と同じ形の末端を生じさせる第五の制限酵素の認識配列を有し、増幅されたプライマーを第四の制限酵素及び第五の制限酵素で消化した後にコンカテマーを形成させる（１）〜（４）のいずれか一項に記載の方法。
【００１９】
（６）前記ベクタープライマーは、第一の制限酵素認識配列のさらに内側に第五の制限酵素認識配列と一塩基異なる塩基配列を有し、前記タグの下流側のプライマーを用いたＰＣＲにより前記一塩基異なる塩基配列が第五の制限酵素認識配列に変換される（５）の方法。
（７）前記第三の制限酵素、第四の制限酵素及び第五の制限酵素が同一である前記（６）の方法。
【００２０】
（８）ベクタープライマーが、マルチクローニング部位を有するプラスミドを同マルチクローニング部位の二箇所で切断して得られる直鎖状プラスミドと、この直鎖状プラスミドの一方の末端と同じ形の末端を有し、一本鎖のポリＴ配列を有する部分二本鎖ＤＮＡとが連結されたものである（１）〜（７）のいずれかの方法。
【００２１】
【発明の実施の形態】
以下に、本発明をさらに詳細に説明する。
本発明の方法は、上記（ａ）〜（ｆ）のステップを含む。以下、各ステップ毎に説明する。
【００２２】
＜１＞ステップ（ａ）
最初のステップでは、直鎖状のプラスミドベクターの一方の３’末端に一本鎖のポリＴ配列と、その内側に第一の制限酵素認識配列とを有し、他方の末端近傍に第二の制限酵素認識配列と、さらにその内側にタイプIIＳ制限酵素認識配列とを有するベクタープライマー（以下、「逆転写用ベクタープライマー」ともいう）を用いる。ここで、タイプIIＳ制限酵素とは、制限酵素により認識される配列から離れた特定の部位を開裂する制限酵素をいう。
【００２３】
ベクタープライマーの構造の一例を、図２Ｄに示す。このベクタープライマーは、例えば、マルチクローニング部位を有するプラスミドを同マルチクローニング部位の二箇所で切断して得られる直鎖状プラスミドと、この直鎖状プラスミドの一方の末端と同じ形の末端を有し、一本鎖のポリＴ配列を有する部分二本鎖ＤＮＡとを連結することにより、調製することができる。図１に、マルチクローニング部位を有するプラスミドの例として、MAGE/pUC19の構造を示す。このプライマーは、基本構造として公知のクローニングベクターであるpUC19を含み、そのマルチクローニングサイトのEcoRlとHindlll制限酵素部位の間に、種々の制限酵素認識部位を含む配列（Z fragmentと呼ぶ。配列を配列番号１に示す。）が挿入されたものである。
尚、以下の具体例では、MAGE/pUC19は、Damメチレーション系を有するエシェリヒア・コリ宿主、例えばJM109株等を用いて調製したものを使用する例について説明するが、本発明に必須ではない。また、基本構造となるベクターは、pUC19に限られるものでなく、それ以外のベクター、例えばpUC18、pBR322、pHSG299、pHSG298、pHSG399、pHSG398、RSF1010、pMW119、pMW118、pMW219、pMW218等、種々のベクターを用いることができる。尚、用いるベクターの宿主は、形質転換、宿主からのベクターの回収等、通常の遺伝子組換え技術を適用することができるものであれば特に制限されるものではないが、通常、エシェリヒア・コリが用いられる。
【００２４】
以下、マルチクローニング部位を有するプラスミドとしてMAGE/pUC19を用いた例について説明するが、本発明はこれに限られるものではない。
まず、MAGE/pUC19ベクター（図２Ａ）を、制限酵素BstXI及びPstIで処理し、切り取られた小断片を除去する（図２Ｂ）。次にBstXIによる切断末端に、一本鎖のポリＴ配列を有する部分二本鎖ＤＮＡ（図２Ｃに示すポリＴアダプター）をライゲーションする。これにより、MAGE/pUC19ベクターの一端から一本鎖のポリＴが飛び出した形のベクタープライマー（又は、逆転写用ベクタープライマー）が形成される（図２Ｄ）。このベクタープライマーでは、第一の制限酵素はPmeIであり、第二の制限酵素はBglIIであり、タイプIIＳ制限酵素はBsgIである。
【００２５】
第二の制限酵素BglIIの認識配列にはMboI認識配列が重複して含まれており、BglII末端をMboI末端と連結すると、MboIで切断することができる。また、この例では、ベクタープライマーの第一の制限酵素（PmeI）認識配列のさらに内側に、第四の制限酵素（MboI）認識配列と一塩基異なる塩基配列を含むように設計されている（図３ＥのΔMboI）。このΔMboI及びベクタープライマーの第一の制限酵素認識配列は、ポリＴアダプターに由来する。
【００２６】
第一の制限酵素及び第二の制限酵素は、ベクタープライマーを一箇所で切断するものであれば特に限定されない。また、タイプIIＳ制限酵素及びその存在位置は、ベクタープライマーから第四の制限酵素部位を切除せず、かつ、cDNAの上流側の一部をベクタープライマーに残したまま切断するものであれば特に制限されない。具体的には、例えば上記BsgIの他に、BsmFI等が挙げられる。
また、ポリＴ配列の長さは、mRNAのポリＡ配列とアニールできる程度の長さであればよく、通常10〜50塩基程度である。
【００２７】
上記のようなベクタープライマーと、遺伝子の発現頻度を解析しようとする細胞に由来するｍＲＮＡとをアニールさせる。その状態で逆転写反応を行うと、ポリＴがプライマーとなってcDNA合成が始まる（図３Ｅ）。そして、合成された一鎖目のcDNAを鋳型として二鎖目を合成し、二本鎖cDNAが合成できる（図３Ｆ）。多数のmRNA分子から多数のｃＤＮＡが結合したベクタープライマー分子（cDNA-MAGE/pUC19）ができるが、図３Fはその典型例のひとつを示したものである。
【００２８】
mRNAは、遺伝子の発現頻度を解析しようとする細胞から抽出する。発現頻度を解析しようとする細胞は、動物や植物の組織の細胞、酵母等の微生物の細胞等、細胞中のmRNAの３’末端にポリＡ構造を有するものであれば特に制限なく使用することができる。また、原核生物のmRNAは３’末端にポリＡ（poly(A)）構造を持たないために、そのままではベクタープライマーのポリＴにアニールさせることはできないが、mRNAに酵素的にポリＡ構造を付加することにより、真核生物のmRNAと同様にして本発明の方法を実施することができる。
mRNAの調製、及びｃＤＮＡ合成、オリゴヌクレオチドの合成、制限酵素反応、ライゲーション反応、ポリメラーゼ・チェイン・リアクション（ＰＣＲ）、形質転換等の操作は、通常のｃＤＮＡクローニングに用いられるmRNAの調製法（例えばSambrook, J., Fritsch, E. F., and Maniatis, T., "Molecular Cloning A Laboratory Manual, Second Edition", Cold Spring Harbor Laboratory Press,(1989)等参照）と同様にして行うことができる。
【００２９】
＜２＞ステップ（ｂ）
次に、前記のようにして得られるｃＤＮＡが結合したベクタープライマーを、ベクタープライマーを切断せず、かつ、前記第二の制限酵素の切断末端と同じ形の末端を生じさせる第三の制限酵素と、前記第二の制限酵素とを用いて消化してｃＤＮＡの上流側を切除し、ベクタープライマーを閉環させる。閉環されたベクタープライマーは、それを適当な宿主に導入し、得られる形質転換体を培養してプラスミドを回収することにより、増幅させることができる。尚、以降の工程で得られる閉環状のベクタープライマーも同様であるが、回収されたベクタープライマーを切断しようとする制限酵素部位が、Damメチレーションにより切断されない場合には、修飾系を有しない宿主を用いる。
【００３０】
前記第三の制限酵素としては４塩基認識の酵素が好ましい。６塩基認識の酵素であると、cDNA配列中に制限酵素部位が存在しない場合がある。また、ベクタープライマーに残されるcDNAの配列が長いと、後の操作で得られるタグがmRNAのポリＡ配列から遠くなる。そのような場合、データベース中の遺伝子の発現情報（EST:expressed sequence tag）は、通常mRNAの３’末端側の一部であるため、タグ配列を検索してもヒットしない場合がある。本発明に好適な４塩基認識の制限酵素としては、MboI、TaiI等が挙げられる。
以下、第三の制限酵素として、MboIを用いた例を説明する。まず、ｃＤＮＡが結合したベクタープライマー（cDNA-MAGE/pUC19）を制限酵素BglII及びMboIで消化する（図３Ｇ）。このとき、MAGE/pUC19はDamメチレーションを受けているので、MboIでは切断されないが、逆転写により新たに合成されたcDNA中のMbolサイトだけが切断される。図３Ｇは、仮にcDNAが３箇所のMbolサイトを含む例を示している。cDNAに注目すると、ポリＡテールから上流に向かって最初のMbolサイトまでの下流側部分がMAGE/pUC19ベクタープライマーにつながったまま残り、他の部分、すなわちｃＤＮＡの上流側はMbol切断により分離除去される。また、MAGE/pUC19は、BglIIにより一箇所のみが切断される。MboI及びBglIIで切断された末端は同じ形をしており、ライゲーション反応で互いをつなぐことができる。そこで、これらの末端を介して、cDNA-MAGE/pUC19をセルフライゲーション反応により閉環する（図３Ｈ）。
【００３１】
＜３＞ステップ（ｃ）
ステップ（ｂ）で閉環したベクタープライマーを、前記第一の制限酵素及び前記タイプIIＳ制限酵素で消化して、ｃＤＮＡの一部からなるタグを残してｃＤＮＡの下流側を切除し、ベクタープライマーを再び閉環する。
【００３２】
具体的には、閉環されたcDNA-MAGE/pUC19を制限酵素BsgI及びPmeIで消化する（図４Ｉ）。これらの制限酵素の切断により、 cDNAの前記３’末端側部分のうちポリＡテールから最も遠い約１３塩基のみがベクタープライマーに残る。このｃＤＮＡの約１３塩基からなる配列を、タグ（図４Ｊ中の「tag」で示される）と呼ぶ。
【００３３】
上記のBsgIによる切断点は、５’側が２ｂｐ突出しているので、例えばT4 DNAポリメラーゼ処理により平滑化する。また、PmeIによる他方の切断末端は、平滑末端であるので、これらの末端をライゲーションすることができる（図４Ｊ）。そこで、ベクタープライマーをセルフライゲーション反応で閉環すると、タグ以外の部分のcDNAが切り取られ、MAGE/pUC19に短いタグがつながった構造のベクタープライマーが得られる（図４Ｊ、Ｋ）。
【００３４】
＜４＞ステップ（ｄ）
次に、ステップ（ｃ）で閉環されたベクタープライマーを鋳型とし、ベクタープライマー中のタグの両側のそれぞれの領域に相当する塩基配列を有するオリゴヌクレオチドをプライマーとしてＰＣＲを行い、タグを増幅する。
【００３５】
具体的には、タグを含むベクタープライマーを鋳型とし、タグの両側のベクター部分の配列に相当するプライマー（例えば配列番号２及び３）を用いてPCR反応を行う（図４Ｋ）。その結果、中央にタグを含むDNAフラグメント（この場合はプライマーの配列及びタグを含めて110bp）が増幅される（図４Ｌ）。
【００３６】
＜５＞ステップ（ｅ）
前記のＰＣＲ増幅産物を連結してタグのコンカテマーを形成させる。その際、前記（ｄ）のステップで用いるプライマーのうち、タグの下流側のプライマーに、第四の制限酵素切断末端と同じ形の末端を生じさせる第五の制限酵素の認識配列を含ませておき、増幅されたプライマーを第四の制限酵素及び第五の制限酵素で消化すると、効率よくコンカテマーを形成させることができる。第五の制限酵素の種類は特に制限されないが、第四の制限酵素と同じ制限酵素を用いると、操作が簡便となる。また、第四の制限酵素として、第三の制限酵素と同じ制限酵素を用いてもよい。さらに、第三、第四及び第五の制限酵素が同じであってもよい。ただし、第五の制限酵素と第四の制限酵素に同じものを用いる場合は、ステップ（ｂ）でcDNAの全配列がベクタープライマーから切除されてしまうので、それを防ぐために、ベクタープライマーの第一の制限酵素認識配列のさらに内側に、第五の制限酵素認識配列と一塩基異なる塩基配列を含ませておき、タグの下流側に第五の制限酵素認識配列を有するプライマーを用いてＰＣＲを行うことにより、一塩基異なる塩基配列を第五の制限酵素認識配列に変換するようにする必要がある。また、第三の制限酵素と第四の制限酵素が同じである場合には、ベクタープライマーの第四の制限酵素認識部位は、メチレーションにより切断されないようにしておく必要がある。その場合には、第四の制限酵素で消化する前に、ステップ（ｃ）で閉環されたベクタープライマーを修飾系を持たない宿主に導入して脱修飾する。
【００３７】
タグの連結反応を、一方の末端がタグの末端と同じ形を有するアダプター存在下で行い、それによりコンカテマーの両端にアダプターを配置させると、アダプターの配列に相当する配列を有するオリゴヌクレオチドをプライマーに用いてＰＣＲを行うことにより、コンカテマーを増幅することができる。アダプターは、タグの量に比べて少量用いるか、一方の末端がタグの末端と異なる形を有するアダプターを用いると、多数のタグからなるコンカテマーの末端にアダプターを配置させることができる。このようにすることで、ごく微量のコンカテマーでもクローニングすることができ、結果としてごく微量のサンプルmRNAを用いた解析が可能となる。
タグとアダプターの比率は、通常モル比でタグ：アダプター＝１：１〜１：０．０１、好ましくは１：０．２〜１：０．０５である。この範囲でアダプターを用いると、タグが２〜５０程度連結したコンカテマーが得られる。
【００３８】
得られるコンカテマーは、塩基配列決定用のクローニングベクターに挿入すると、シークエンス操作が容易になる。
MAGE/pUC19を用いる場合には、ステップ（ｄ）の増幅産物を制限酵素MboIで消化することにより、タグを含む小断片が得られる（図５Ｍ）。ここで得られたタグ断片は、既知のDNA配列、すなわちGATCとAAACGに挟まれており、どの部分がタグかを明確に識別することができる。このタグ配列は両端にMbolサイトが露出しているので、ライゲーション反応によりタグ同士を連結し、コンカテマーを形成させることができる（図５Ｍ）。
【００３９】
＜６＞ステップ（ｆ）
上記のようにして得られるコンカテマーの塩基配列を決定すると、その中には多数のcDNA分子に由来するタグが含まれているので、その塩基配列中に出現するタグの種類及びその頻度を調べることにより、cDNAが由来する遺伝子の発現頻度を解析することができる。タグの種類は、既知のmRNAの部分配列（EST）の情報に関するデータベースを検索することにより調べることができる。
【００４０】
尚、タグのコンカテマーを増幅することにより、同一配列を持つフラグメントが複数個生じることになる。しかし、クローニング後のシークエンス解析において、先にシークエンスしたものと同一物とみなされるシークエンスは出現した場合は、それを解析から除くことにより、遺伝子発現頻度解析にPCRが及ぼす悪影響を排除することができる。
また、第四と第五の制限酵素として同じ酵素を用いた場合は、コンカテマー中のタグの向きは不定であるので、配列解析においては一方の鎖だけでなく逆鎖の配列も考慮する。
【００４１】
以上説明したように、本発明の方法では、SAGE法のようにダイタグを作らず、個々のタグを既知のDNA断片で挟むことにより、タグが他のタグと連続しないようにすることができる。その結果、タグの境界が曖昧になるという問題点が解決される。また、ＰＣＲを複数回行い、試料を増幅することができるので、微量のmRNAからでも解析が可能となる。さらに、ベクタープライマーを用いることにより、ベクターに融合した形のcDNA合成を行うことで、アビジンやビオチンビーズを用いずに解析することができる。
【００４２】
本発明の遺伝子発現解析法は、生命科学の研究にとって有用であるが、特に、健康人と病気の人における特定臓器、又は細胞における遺伝子発現の差異を解析する事により、病気の治療法開発や診断法開発に有用である。例えば、健康人の肝臓と肝炎の肝臓の発現遺伝子の違いを本法により解析することにより、肝炎で特異的に発現が増加したり減少したりする遺伝子を見つけることができる。これら遺伝子の肝臓での役割を調べることにより、肝炎治療を目的として、この遺伝子の働きを阻害したり、促進させたりする薬剤を開発することができる。又、遺伝子そのもの、遺伝子構造から設計したアンチセンスオリゴヌクレオチド、あるいは遺伝子を発現させて得た蛋白質等を、肝炎治療に用いることができる。
【００４３】
また、発症のメカニズムが解っていない病気に対しても、本発明の方法を用いれば、治療法を開発できる可能性がある。さらに、本発明の方法により疾患特異的に発現が変動する遺伝子を発見することができれば、治療ばかりでなく、疾患の診断法を開発することも可能となる。
【００４４】
さらに、医療における利用ばかりでなく、全ての真核生物を対象に有用遺伝子を発見する手段ともなる。例えば、ビール生産に適した酵母を変異により育種した場合、親株と変異株の遺伝子発現の変化を本発明の方法で解析し、変異により発現が変化した遺伝子を同定することができる。このようにして得たビール生産にとって有利な遺伝子を総合的に操作することにより、よりよいビール生産酵母を創製することができる。
【００４５】
また、例えば、大腸菌やコリネバクテリウムのようなアミノ酸生産菌の解析を行うことにより、より優れたアミノ酸生産菌を創製することができる。
【００４６】
【実施例】
マウスC57BL/6肝臓 1gから、Invitrogen社製FastTrack 2.0キットを用いてmRNA 35μgを抽出した。
得られたmRNA 0.97μgを用いて、上記と同様にしてMAGE/pUC19を用いて逆転写cDNA合成後、遺伝子発現頻度解析を行った。
【００４７】
MAGE/pUC19ベクター（図２Ａ）を、制限酵素BstXI及びPstIで処理し、切り取られた小断片を除去した（図２Ｂ）。次にBstXIによる切断末端に、一本鎖のポリＴ配列を有する部分二本鎖ＤＮＡ（図２Ｃのに示すポリＴアダプター）をライゲーションした（図２Ｄ）。このベクタープライマーでは、第一の制限酵素はPmeIであり、第二の制限酵素はBglIIであり、タイプIIＳ制限酵素はBsgIである。
上記ベクタープライマー0.2μgと、上記のマウス肝臓由来のmRNA 0.97μgとをアニールさせ、cDNA合成を行い、続いて合成された一鎖目のcDNAを鋳型として二鎖目を合成した（図３Ｆ）。そして、得られたcDNAの40分の１に当たる量（mRNA 0.025μg相当）を材料に遺伝子発現頻度解析を行った。
尚、別の実験では、4.5mgの微量臨床サンプルから抽出した、0.05μgの微量mRNAを実験材料にして、逆転写cDNA合成後、以下と同様にして遺伝子発現頻度解析を実施することに成功している。
【００４８】
次に、前記のようにして得られるcDNAが結合したベクタープライマーを、第二の制限酵素BglII及び第三の制限酵素であるMboIで消化した（図３Ｇ）。次に、cDNA-MAGE/pUC19のMboI末端及びBglII末端をセルフライゲーション反応により連結して閉環させた（図３Ｈ）。
【００４９】
前記の閉環したベクタープライマーをBsgI及びPmeIで消化し、T4 DNAポリメラーゼ処理により平滑化した後、ベクタープライマーを再び閉環させた（図４Ｉ〜Ｋ）。この閉環されたベクタープライマーを鋳型とし、ベクタープライマー中のタグの両側のそれぞれの領域に相当する塩基配列を有するオリゴヌクレオチド（配列番号２及び３）をプライマーとして、酵素AmpliTaq Gold（PE Biosystems社）を用いてＰＣＲを行い、タグを増幅した（図４Ｌ）。ＰＣＲ反応は、変性（95℃、0.3分）、アニーリング、およびポリメラーゼによる伸長反応（72℃、1.5分）からなる反応を65サイクル行うことにより行った。その結果、中央にタグを含む110bpのDNAフラグメントが増幅された。
【００５０】
前記のＰＣＲ増幅産物をMboIで消化してタグを切り出し、タグ：アダプターが８：１となるように、配列番号４、５に示す塩基配列を有するオリゴヌクレオチドをアニールさせたアダプターを加えて連結反応を行い、タグのコンカテマーを形成させた（図５Ｍ）。このコンカテマーを、配列番号６に示す塩基配列を有するプライマーを用いたＰＣＲにより増幅した。ＰＣＲ反応は、変性（95℃、0.3分）、アニーリング（40℃、3分）、ポリメラーゼによる伸長反応（72℃、1分）からなる反応を５サイクル行い、引き続き変性（95℃、0.3分）、アニーリング及び伸長反応（72℃、1分）からなる反応を60サイクル行った。増幅産物を制限酵素NotIで消化し、シークエンス用クローニングベクターpKF3のNotI部位に挿入して、塩基配列を決定した。
【００５１】
得られたシークエンスデータ（11000個）から、同一の配列を有するタグをグループ化し、出現頻度順に並べた。頻度上位１０番までを表１に示す。この結果から、マウス肝臓組織中で最も発現頻度が高い遺伝子はウリナリ・プロテインＩ／II遺伝子であることがわかる。以下、発現頻度の高い順に、アルブミン遺伝子、ウリナリ・プロテインIII遺伝子、アルギニノコハク酸合成酵素遺伝子と並んでいる。このように、本発明の方法は、mRNAサンプルが微量であっても、ベクタープライマーを用いることによりcDNAの回収効率が高く、またタグ配列の不明瞭さを除くことが可能であるので、遺伝子発現頻度を高精度で解析することができる。
【００５２】
【表１】
【発明の効果】
本発明により、遺伝子の発現頻度を、簡便に、確実に、かつ精度よく解析することができる。
【００５４】
【配列表】
【図面の簡単な説明】
【図１】 ベク夕ープライマーを作るためのプラスミドDNAの構造の一例（MAGE/pUC19）を示す図。
【図２】 ベクタープライマーの構造の一例とその構築過程を示す図。
【図３】 図３〜図５は、本発明の方法の各ステップを示す。図３は、ステップ（ａ）、（ｂ）を模式的に表す図。
【図４】 本発明の方法のステップ（ｃ）、（ｄ）を模式的に表す図。
【図５】 本発明の方法のステップ（ｅ）と、ステップ（ｅ）により得られる増幅産物をシークエンス用クローニングベクターに挿入するステップを模式的に表す図。
Claims (8)
- 以下のステップを含む遺伝子の発現頻度を解析する方法；
（ａ）直鎖状のプラスミドベクターの一方の３’末端に一本鎖のポリＴ配列と、その内側に第一の制限酵素認識配列とを有し、他方の末端近傍に第二の制限酵素認識配列と、さらにその内側にタイプIIＳ制限酵素認識配列とを有するベクタープライマーと、遺伝子の発現頻度を解析しようとする細胞に由来するｍＲＮＡとをアニールさせてｃＤＮＡ合成を行い、ｃＤＮＡが結合したベクタープライマーを生成するステップと、
（ｂ）前記ｃＤＮＡが結合したベクタープライマーを、ベクタープライマーを切断せず、かつ、前記第二の制限酵素の切断末端と同じ形の末端を生じさせる第三の制限酵素と、前記第二の制限酵素とを用いて消化してｃＤＮＡの上流側を切除し、ベクタープライマーを閉環するステップと、
（ｃ）閉環したベクタープライマーを、前記第一の制限酵素及び前記タイプIIＳ制限酵素で消化して、ｃＤＮＡの一部からなるタグを残してｃＤＮＡの下流側を切除し、ベクタープライマーを再び閉環するステップと、
（ｄ）前記ベクタープライマーを鋳型とし、ベクタープライマー中のタグの両側のそれぞれの領域に相当する塩基配列を有するオリゴヌクレオチドをプライマーとしてポリメラーゼ・チェイン・リアクション（ＰＣＲ）を行い、タグを増幅するステップと、
（ｅ）前記増幅産物を連結してタグのコンカテマーを形成するステップと、
（ｆ）前記コンカテマーの塩基配列を決定し、その塩基配列中に出現するタグの種類及びその頻度を調べるステップ。 - 前記ステップ（ｅ）において、連結反応を、一方の末端がタグの末端と同じ形を有するアダプター存在下で行い、それによりコンカテマーの両端にアダプターを配置させ、アダプターの配列に相当する配列を有するオリゴヌクレオチドをプライマーに用い、コンカテマーを鋳型としてＰＣＲを行うことによりコンカテマーを増幅することを特徴とする請求項１記載の方法。
- 前記ステップ（ｅ）の後に、前記コンカテマーの塩基配列の決定を、塩基配列決定用のクローニングベクターにクローニングして行うことを特徴とする請求項１又は２記載の方法。
- 前記第三の制限酵素の認識配列が４塩基である請求項１〜３のいずれか一項に記載の方法。
- ベクタープライマーは、第二の制限酵素認識配列の切断点と同じ位置又はその内側に切断点を有し、かつ、タイプIIＳ制限酵素による切断によってベクタープライマーから切除されない第四の制限酵素認識配列を有し、
前記（ｄ）のステップで用いるプライマーのうち、タグの下流側のプライマーは、第四の制限酵素切断末端と同じ形の末端を生じさせる第五の制限酵素の認識配列を有し、増幅されたプライマーを第四の制限酵素及び第五の制限酵素で消化した後にコンカテマーを形成させる請求項１〜４のいずれか一項に記載の方法。 - 前記ベクタープライマーは、第一の制限酵素認識配列のさらに内側に第五の制限酵素認識配列と一塩基異なる塩基配列を有し、前記タグの下流側のプライマーを用いたＰＣＲにより前記一塩基異なる塩基配列が第五の制限酵素認識配列に変換される請求項５記載の方法。
- 前記第三の制限酵素、第四の制限酵素及び第五の制限酵素が同一である請求項６記載の方法。
- ベクタープライマーが、マルチクローニング部位を有するプラスミドを同マルチクローニング部位の二箇所で切断して得られる直鎖状プラスミドと、この直鎖状プラスミドの一方の末端と同じ形の末端を有し、一本鎖のポリＴ配列を有する部分二本鎖ＤＮＡとが連結されたものである請求項１〜７のいずれか一項に記載の方法。
Priority Applications (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03853899A JP3924976B2 (ja) | 1999-02-17 | 1999-02-17 | 遺伝子の発現頻度の解析方法 |
| DK00904009T DK1156107T3 (da) | 1999-02-17 | 2000-02-17 | Fremgangsmåde til at analysere genekspressionsfrekvens |
| EP00904009A EP1156107B1 (en) | 1999-02-17 | 2000-02-17 | Method for analyzing gene expression frequency |
| AT00904009T ATE275196T1 (de) | 1999-02-17 | 2000-02-17 | Verfahren zur analyse der genexpressions-frequenz |
| PCT/JP2000/000902 WO2000049143A1 (fr) | 1999-02-17 | 2000-02-17 | Procede d'analyse d'une frequence d'expression genique |
| US09/926,028 US6806049B1 (en) | 1999-02-17 | 2000-02-17 | Method for analyzing gene expression frequency |
| DE60013410T DE60013410T2 (de) | 1999-02-17 | 2000-02-17 | Verfahren zur analyse der genexpressions-frequenz |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03853899A JP3924976B2 (ja) | 1999-02-17 | 1999-02-17 | 遺伝子の発現頻度の解析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2000232888A JP2000232888A (ja) | 2000-08-29 |
| JP3924976B2 true JP3924976B2 (ja) | 2007-06-06 |
Family
ID=12528071
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP03853899A Expired - Fee Related JP3924976B2 (ja) | 1999-02-17 | 1999-02-17 | 遺伝子の発現頻度の解析方法 |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US6806049B1 (ja) |
| EP (1) | EP1156107B1 (ja) |
| JP (1) | JP3924976B2 (ja) |
| AT (1) | ATE275196T1 (ja) |
| DE (1) | DE60013410T2 (ja) |
| DK (1) | DK1156107T3 (ja) |
| WO (1) | WO2000049143A1 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004097158A (ja) * | 2002-09-12 | 2004-04-02 | Kureha Chem Ind Co Ltd | 発現遺伝子同定用ｃＤＮＡタグの作成方法、及び該ｃＤＮＡタグを用いる遺伝子発現解析方法 |
| JP4034317B2 (ja) | 2005-03-29 | 2008-01-16 | 富士通株式会社 | 画像符号化装置 |
| AU2007341981A1 (en) * | 2006-12-29 | 2008-07-10 | The Salk Institute For Biological Studies | Methods for enhancing exercise performance |
| WO2009091719A1 (en) | 2008-01-14 | 2009-07-23 | Applera Corporation | Compositions, methods, and kits for detecting ribonucleic acid |
| EP2910649A1 (en) | 2009-08-24 | 2015-08-26 | National University Corporation Kanazawa University | Detection of pancreatic cancer by gene expression profiling |
| JP7445334B1 (ja) | 2022-09-05 | 2024-03-07 | 株式会社キュービクス | 膵臓癌に特異的な遺伝子発現パターンの検出及びｃａ１９－９の測定の併用による膵臓癌の検出 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5866330A (en) * | 1995-09-12 | 1999-02-02 | The Johns Hopkins University School Of Medicine | Method for serial analysis of gene expression |
| US5695937A (en) | 1995-09-12 | 1997-12-09 | The Johns Hopkins University School Of Medicine | Method for serial analysis of gene expression |
| US6136537A (en) * | 1998-02-23 | 2000-10-24 | Macevicz; Stephen C. | Gene expression analysis |
-
1999
- 1999-02-17 JP JP03853899A patent/JP3924976B2/ja not_active Expired - Fee Related
-
2000
- 2000-02-17 WO PCT/JP2000/000902 patent/WO2000049143A1/ja active IP Right Grant
- 2000-02-17 DE DE60013410T patent/DE60013410T2/de not_active Expired - Lifetime
- 2000-02-17 US US09/926,028 patent/US6806049B1/en not_active Expired - Fee Related
- 2000-02-17 DK DK00904009T patent/DK1156107T3/da active
- 2000-02-17 EP EP00904009A patent/EP1156107B1/en not_active Expired - Lifetime
- 2000-02-17 AT AT00904009T patent/ATE275196T1/de not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| EP1156107A1 (en) | 2001-11-21 |
| US6806049B1 (en) | 2004-10-19 |
| ATE275196T1 (de) | 2004-09-15 |
| DE60013410T2 (de) | 2005-09-15 |
| EP1156107A4 (en) | 2002-09-25 |
| DE60013410D1 (de) | 2004-10-07 |
| WO2000049143A1 (fr) | 2000-08-24 |
| EP1156107B1 (en) | 2004-09-01 |
| DK1156107T3 (da) | 2004-12-06 |
| JP2000232888A (ja) | 2000-08-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5345857B2 (ja) | 正常および異常遺伝子発現におけるｄｎａの立体構造（ループ構造） | |
| EP0656954B1 (en) | Method for producing tagged genes, transcripts, and proteins | |
| CN110719958B (zh) | 构建核酸文库的方法和试剂盒 | |
| EP0791074B1 (en) | Method for the detection of ras oncogenes, in particular the k-ras ancogene | |
| EP3854873A1 (en) | Compositions and methods for accurately identifying mutations | |
| EP1394272A1 (en) | Method for detection of Ki-ras mutations | |
| JP4644685B2 (ja) | 塩基配列タグの調製方法 | |
| AU2014305994A1 (en) | Urine biomarker cohorts, gene expression signatures, and methods of use thereof | |
| JP7006873B2 (ja) | 突然変異細胞遊離遺伝子分離キット及びそれを用いた突然変異細胞遊離遺伝子分離方法 | |
| WO2019080595A1 (zh) | 一种原位杂交探针的制备方法 | |
| JP3924976B2 (ja) | 遺伝子の発現頻度の解析方法 | |
| KR101913735B1 (ko) | 차세대 염기서열 분석을 위한 시료 간 교차 오염 탐색용 내부 검정 물질 | |
| JP2004154088A (ja) | Ｃｅａ検出のための核酸増幅用プライマーおよび該プライマーを用いた検査方法 | |
| CN108796077B (zh) | 在cfDNA中检测胞嘧啶脱氨酶及相关分子基因修饰差异的引物对组及试剂盒 | |
| US20230104375A1 (en) | Method for multiplexable strand-specific 3' end sequencing of mrna transcriptome primer set, kit and application thereof | |
| CN108841958A (zh) | 一种检测EGFR基因外显子21L858突变的crRNA及其检测方法 | |
| JP2004105013A (ja) | 肝細胞癌検出法 | |
| JP2007228803A (ja) | 転写機構を利用した知能型核酸配列センシングシステム | |
| CN112063690A (zh) | 单分子探针多重靶向捕获文库的构建方法及应用 | |
| JP2004187606A (ja) | 核酸アイソフォームの同定、分析および／またはクローニング方法 | |
| JP4403069B2 (ja) | クローニングおよび分析のためのｍＲＮＡの５’末端の使用方法 | |
| CN113699221B (zh) | HER2 mRNA及环状RNA多重荧光定量PCR检测引物探针及其应用 | |
| WO2008015975A1 (fr) | Procédé d'amplification d'un fragment d'adn | |
| JP5322141B2 (ja) | ５’領域ディファレンシャルディスプレー法 | |
| JP6417603B2 (ja) | Ｒｎａの末端領域に対応する核酸の塩基配列を解読する方法およびｄｎａエレメントの分析方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040917 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070206 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070219 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100309 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100309 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100309 Year of fee payment: 3 |
|
| R255 | Notification that request for automated payment was rejected |
Free format text: JAPANESE INTERMEDIATE CODE: R2525 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110309 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110309 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120309 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120309 Year of fee payment: 5 |
|
| LAPS | Cancellation because of no payment of annual fees |