(57)【要約】  【課題】  一種類の検出用核酸に複数種類の標的核酸がハイブリダイゼーションする場合にも標的遺伝子の発現量を推定できる手段を提供すること。  【解決手段】  基板表面に設けられた反応領域に固定された検出用核酸と、前記検出用核酸とハイブリダイゼーションする標的核酸とのハイブリダイゼーション測定値からサンプル中における標的核酸含有量を推定する方法であって、一つの検出用核酸と複数の標的核酸とのそれぞれの結合強度を用いて、前記ハイブリダイゼーション測定値における、各標的核酸含有量の寄与の度合いを取得することにより、サンプル中における標的核酸含有量を推定する方法を提供する。  【選択図】  図３
2006191844.tif 000002
【特許請求の範囲】
【請求項１】
基板表面に設けられた反応領域に固定された検出用核酸と、前記検出用核酸とハイブリダイゼーションする標的核酸とのハイブリダイゼーション測定値からサンプル中における標的核酸含有量を推定する方法であって、
一つの検出用核酸と複数の標的核酸のそれぞれとの間における各結合強度を用いて、前記ハイブリダイゼーション測定値における、各標的核酸含有量の寄与の度合いを取得することにより、サンプル中における標的核酸含有量を推定する方法。
【請求項２】
複数の反応領域で測定された各ハイブリダイゼーション測定値と、
前記複数の反応領域に固定された各検出用核酸のそれぞれに係わる、前記一つの検出用核酸と複数の標的核酸のそれぞれのとの間における各結合強度と、
を用いて、標的核酸含有量の寄与の度合いを取得することを特徴とする請求項１記載の方法。
【請求項３】
前記結合強度は、一つの検出用核酸と複数の標的核酸との結合性の度合いをそれぞれ測定することにより得られた数値であることを特徴とする請求項１記載の方法。
【請求項４】
前記結合強度は、一つの検出用核酸の塩基配列と複数の標的核酸の塩基配列との間におけるそれぞれの相補性の度合いを探索することにより得られた数値であることを特徴とする請求項１記載の方法。
【請求項５】
前記サンプル中における標的核酸含有量とは、遺伝子発現量であることを特徴とする請求項１記載の方法。
【請求項６】
基板表面に設けられた反応領域に固定された検出用核酸と、前記検出用核酸とハイブリダイゼーションする標的核酸とのハイブリダイゼーション測定値からサンプル中における特定の標的核酸含有量を推定するステップであって、
一つの検出用核酸と複数の標的核酸のそれぞれとの間における各結合強度を用いて、前記ハイブリダイゼーション測定値における、各標的核酸含有量の寄与の度合いを取得するステップに係わる、サンプル中における標的核酸含有量を推定するプログラム。
【請求項７】
複数の反応領域において測定された各ハイブリダイゼーション測定値と、
前記複数の反応領域に固定された各検出用核酸のそれぞれに係わる、前記一つの検出用核酸と複数の標的核酸のそれぞれとの間における結合強度と、
を用いて、標的核酸含有量の寄与の度合いを取得することを特徴とする請求項６記載のプログラム。
【請求項８】
前記結合強度は、一つの検出用核酸と複数の標的核酸との結合性の度合いをそれぞれ測定することにより得られた数値であることを特徴とする請求項６記載のプログラム。
【請求項９】
前記結合強度に係る数値を、一つの検出用核酸の塩基配列と複数の標的核酸の塩基配列との間におけるそれぞれの相補性の度合いを探索することにより得るステップに係わるプログラムを備えることを特徴とする請求項６記載のプログラム。
【請求項１０】
前記サンプル中における標的核酸含有量とは、遺伝子発現量であることを特徴とする請求項６記載のプログラム。
【請求項１１】
基板表面に設けられた反応領域に固定された検出用核酸と、前記検出用核酸とハイブリダイゼーションする標的核酸とのハイブリダイゼーション測定値からサンプル中における標的核酸含有量を推定するシステムであって、
一つの検出用核酸と複数の標的核酸のそれぞれとの間における各結合強度を用いて、前記ハイブリダイゼーション測定値における、各標的核酸含有量の寄与の度合いを取得することにより、サンプル中における標的核酸含有量を推定する標的核酸含有量推定手段を少なくとも備えた標的核酸含有量推定システム。
【請求項１２】
複数の反応領域で測定された各ハイブリダイゼーション測定値と、
前記複数の反応領域に固定された各検出用核酸のそれぞれに係わる、前記一つの検出用核酸と複数の標的核酸のそれぞれとの間における各結合強度と、
からサンプル中における標的核酸含有量を推定することを特徴とする請求項１１記載の標的核酸含有量推定システム。
【請求項１３】
前記結合強度は、一つの検出用核酸と複数の標的核酸との結合性の度合いを測定することにより得られた数値であることを特徴とする請求項１１記載の標的核酸含有量推定システム。
【請求項１４】
前記結合強度に係る数値を、一つの検出用核酸の塩基配列と複数の標的核酸の塩基配列との間におけるそれぞれの相補性の度合いを探索することにより得る結合強度取得手段を少なくとも備えたことを特徴とする請求項１１記載の標的核酸含有量推定システム。
【請求項１５】
前記サンプル中における標的核酸含有量とは、遺伝子発現量であることを特徴とする請求項１１記載の標的核酸含有量推定システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＤＮＡチップなどバイオアッセイ用基板を用いて得られたハイブリダイゼーション測定値などの規格化又は標準化、解析、補正に係わる技術分野に属する。
【背景技術】
【０００２】
近年、ＤＮＡチップ若しくはＤＮＡマイクロアレイ（以下、本願では「ＤＮＡチップ」とする。）の実用化が進んでいる。ＤＮＡチップは、多種・多数のＤＮＡオリゴ鎖を、検出用核酸として基板表面に集積して固定したものである。ＤＮＡチップを用いて、ＤＮＡチップ基板表面の検出用核酸と細胞・組織などより採取した標的核酸とのハイブリダイゼーションを測定することにより、細胞・組織などにおける遺伝子発現などを網羅的に解析することができる。
【０００３】
検出用核酸は、一種類の標的核酸とのみハイブリダイゼーションする塩基配列であることが望ましい。しかしながら、検出用核酸の塩基配列は標的核酸よりも短い場合が多いため、検出用核酸の設計において、目的の標的核酸以外の全ての標的核酸と全く相同性のない配列を探索することは難しい。従って、実際には、ＤＮＡチップを用いたハイブリダイゼーション測定において、一種類の検出用核酸に複数種類の標的核酸がハイブリダイゼーションする場合がある。
【０００４】
それに対し、例えば、目的の標的核酸を認識できる、塩基配列の異なる複数の検出用核酸を基板表面に固定したり、検出用核酸の塩基配列に１〜数塩基分の点変異を加えたりすることにより、目的の標的核酸の認識精度を高める試みが行われている。
【０００５】
なお、本発明に関する先行文献として、例えば、以下のものが挙げられる。
【特許文献１】特表２００１−５１４９０７号公報
【特許文献２】特表２００３−５３２３６７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
上記のように、一種類の検出用核酸に複数種類の標的核酸がハイブリダイゼーションする場合、標的遺伝子の発現量を推定できないという問題があった。例えば、ハイブリダイゼーション測定を蛍光強度に基づいて行う場合、一般的には、蛍光強度と遺伝子発現量が比例していると推定できる。しかしながら、一種類の検出用核酸に複数種類の標的核酸がハイブリダイゼーションする場合、蛍光強度がどの標的核酸に基づくものなのか推定できないため、蛍光強度に基づいた標的遺伝子発現量の推定を行うことが難しかった。
【０００７】
また、一種類の検出用核酸に複数種類の標的核酸がハイブリダイゼーションする場合の中でも、検出用核酸と各標的核酸との結合性にはそれぞれ差があるため、蛍光強度に単純な補正を加えることにより、標的遺伝子の発現量を推定することも難しかった。
【０００８】
そこで、本発明は、一種類の検出用核酸に複数種類の標的核酸がハイブリダイゼーションする場合にも標的遺伝子の発現量を推定できる手段を提供することを主な目的とする。
【課題を解決するための手段】
【０００９】
本発明では、基板表面に設けられた反応領域に固定された検出用核酸と、前記検出用核酸とハイブリダイゼーションする標的核酸とのハイブリダイゼーション測定値からサンプル中における標的核酸含有量を推定する方法であって、一つの検出用核酸と複数の標的核酸のそれぞれとの間における各結合強度を用いて、前記ハイブリダイゼーション測定値における、各標的核酸含有量の寄与の度合いを取得することにより、サンプル中における標的核酸含有量を推定する方法を提供する。
【００１０】
この方法は、例えば、複数の反応領域で測定された各ハイブリダイゼーション測定値と、前記複数の反応領域に固定された各検出用核酸のそれぞれに係わる、前記一つの検出用核酸と複数の標的核酸のそれぞれとの間における各結合強度とを用いて、標的核酸含有量の寄与の度合いを取得することにより、サンプル中における複数の標的核酸含有量を網羅的に推定することができる。
【００１１】
ここで、結合強度は、一つの検出用核酸と複数の標的核酸のそれぞれとの間における結合性（親和性）を示す数値であり、一つの検出用核酸と複数の標的核酸との結合性の度合いをそれぞれ測定することにより得ることができる。また、一つの検出用核酸の塩基配列と複数の標的核酸の塩基配列との間におけるそれぞれの相補性の度合いを探索することにより取得してもよい。
【００１２】
なお、本発明において、サンプル中における標的核酸含有量とは、採取した試料からＲＮＡを取得・調製した場合は遺伝子発現量をいい、採取した試料からＤＮＡゲノムを取得・調製した場合は、標的ゲノムの含有量をいう。
【００１３】
その他、本発明に係る方法は、プログラムで記述することにより、自動化できる。また、本発明は、システム化することができる。
【００１４】
以下、定義づけを行う。
【００１５】
「ハイブリダイゼーション」は、相補的な塩基配列構造を備える核酸間の相補鎖（二本鎖）形成反応を意味する。
【００１６】
「ハイブリダイゼーション測定値」は、ＤＮＡチップの基板表面に固定された検出用核酸と、前記検出用核酸とハイブリダイゼーションする標的核酸とのハイブリダイゼーション量を、蛍光強度などにより測定した実測値（測定データ）、及び、その値に基づいて取得した遺伝子発現量の推定値をいう。
【００１７】
「核酸」とは、プリンまたはピリミジン塩基と糖がグリコシド結合したヌクレオシドのリン酸エステルの重合体（ヌクレオチド鎖）を意味し、プローブＤＮＡを含むオリゴヌクレオチド、ポリヌクレオチド、プリンヌクレオチドとピリミジンヌクレオチオドが重合したＤＮＡ（全長あるいはその断片）、逆転写により得られるｃＤＮＡ（ｃプローブＤＮＡ）、ＲＮＡ、ポリアミドヌクレオチド誘導体（ＰＮＡ）等を広く含む。
【００１８】
「検出用核酸」とは、反応領域に貯留又は保持された媒質中に固定又は遊離の状態で存在し、当該核酸分子と特異的に相互作用する相補的塩基配列を有する核酸分子を検出するための探り針（検出子）として機能する核酸分子である。代表例は、ＤＮＡプローブなどのオリゴヌクレオチド又はポリヌクレオチドである。「標的核酸」は、細胞から取得されたサンプル核酸のうち、前記検出用核酸とハイブリダイゼーションする核酸である。
【発明の効果】
【００１９】
本発明により、一種類の検出用核酸に複数種類の標的核酸がハイブリダイゼーションする場合にもサンプル中における標的核酸含有量（例えば遺伝子発現量など）を推定できる。
【発明を実施するための最良の形態】
【００２０】
本発明を実施するための好適な形態について、添付図面を参照しながら説明する。なお、以下に記載する実施形態は、ＤＮＡチップの基板表面に固定された検出用核酸と採取した試料から取得したサンプル核酸中の標的核酸とのハイブリダイゼーション量を蛍光強度により取得した場合について例示したものであるが、本発明の範囲は、これによって狭く解釈されない。
【００２１】
はじめに、図１から図５を用いて、本発明に係るサンプル中における標的核酸含有量（遺伝子発現量など）の推定方法の概要について説明する。
【００２２】
図１は、ＤＮＡチップを用いたハイブリダイゼーション測定手順の一例を模式的に示した図である。
【００２３】
ＤＮＡチップの基板表面１１には、多種・多数の検出用核酸１２が固定されている。ＤＮＡチップを用いたハイブリダイゼーション測定の際には、まず、個体１３より採取した試料からサンプル核酸１４を抽出・調製などし、次に、サンプル核酸１４をＤＮＡチップの基板表面１１に滴下・注入などする（矢印１５）。そして、検出用核酸１２とサンプル核酸１４中に含まれる標的核酸とのハイブリダイゼーション量を、蛍光強度などにより測定する。
【００２４】
図２は、ハイブリダイゼーション測定値とサンプル中における標的核酸含有量との関係を模式的に示した図である。
【００２５】
検出用核酸２１が標的核酸２２とのみハイブリダイゼーションする場合、検出用核酸２１に係るハイブリダイゼーション測定値ｐ１は「数１」に示す式で表すことができる。即ち、サンプル中における標的核酸２２の含有量（遺伝子発現量など）ｇ１は、ハイブリダイゼーション測定値ｐ１と、検出用核酸２１と標的核酸２２の結合強度ａ１１から取得することができる。
【数１】
2006191844.tif 000003
【００２６】
しかしながら、一般に検出用核酸は標的核酸よりも塩基配列が短いため、例えば、検出用核酸２１の塩基配列と、標的核酸２２以外の他の核酸２３の塩基配列の一部（又は全部）とが相補性を有する場合がある。その場合における、検出用核酸２１に係るハイブリダイゼーション測定値は「数２」に示す式で表される。「数２」の式では、ハイブリダイゼーション測定値ｐ１とサンプル中における標的核酸２２の含有量（遺伝子発現量など）ｇ１とが比例しないため、ハイブリダイゼーション測定値ｐ１と結合強度ａ１１とから、サンプル中における標的核酸２１の含有量（遺伝子発現量など）ｇ１を取得することはできない。
【数２】
2006191844.tif 000004
【００２７】
なお、図２において、核酸２３は検出用核酸２４の標的核酸であり、結合強度ａ２２は検出用核酸２４と標的核酸２３との結合強度を表している。また、ｇ２はサンプル中における標的核酸２３の含有量、ｐ２は検出用核酸２４に係るハイブリダイゼーション測定値である。
【００２８】
図３は、ハイブリダイゼーション測定値を、サンプル中における各標的核酸含有量に結合強度を乗じたものの総和として表すことができることを示した模式図である。
【００２９】
例えば、ＤＮＡチップの基板表面に固定されたｎ個の検出用核酸のうち、ある特定の検出用核酸ｉの塩基配列と、前記ｎ個の検出用核酸のそれぞれに係る標的核酸ｇｎの塩基配列との、それぞれの相補性の度合いを結合強度ａｉｎとした場合、ある特定の検出用核酸ｉに係るハイブリダイゼーション測定値Ｐｉは、「数３」で示す式で表される。
【数３】
2006191844.tif 000005
【００３０】
上記式が示すように、一つの検出用核酸と複数の標的核酸とのそれぞれの結合強度ａｉｎを用いて、ハイブリダイゼーション測定値Ｐｉにおける、各標的核酸含有量の寄与の度合いを取得することにより、サンプル中における標的核酸含有量ｇｎを推定することができる。
【００３１】
図４は、本発明に係るサンプル中における標的核酸含有量（遺伝子発現量など）の推定方法の一例のフローを示す図である。なお、本フローは、プログラムで記述することにより自動化できる。
【００３２】
前記の通り、ある特定の検出用核酸ｉに係るハイブリダイゼーション測定値Ｐｉは、図４中の符号４１に示す式で表すことができる。この式は、全てのハイブリダイゼーション測定値（ｐ１〜ｐｎ）について、一つの行列形式として表すことができる（符号４２）。なお、符号４２に示す式中、ｐはハイブリダイゼーション測定値（ｐ１〜ｐｎ）を、ｓはスケールを、Ａは結合強度行列（ａ１１〜ａｍｎ）を、ｇはサンプル中における各標的核酸含有量（ｇ１〜ｇｎ）を、それぞれ示す。
【００３３】
符号４２に示す式について、両辺の左からそれぞれ結合強度行列Ａの逆行列Ａ−１を乗ずることにより、サンプル中における各標的核酸含有量ｇを取得する行列形式（「数４」で示す式）を得ることができる。即ち、複数のハイブリダイゼーション測定値ｐ（ｐ１〜ｐｎ）と結合強度行列Ａから、サンプル中における各標的核酸含有量ｇ（ｇ１〜ｇｎ）を取得する行列形式を得ることができる。
【数４】
2006191844.tif 000006
【００３４】
但し、結合強度行列Ａは、多くの場合、正方行列（ｍ＝ｎ）ではないため、一般的には逆行列が存在しない。そこで、特異値分解（Ａ＝ＵＤｔＶ）で、結合強度行列Ａの一般逆行列Ａ−１を取得し（符号４３）、サンプル中における各標的核酸含有量ｇを推定する（符号４４）。
【００３５】
以上のように、複数の反応領域で測定された各ハイブリダイゼーション測定値ｐと、前記複数の反応領域に固定された各検出用核酸のそれぞれに係わる、前記一つの検出用核酸と複数の標的核酸とのそれぞれの結合強度Ａとを用いて、標的核酸含有量の寄与の度合いを取得することにより、サンプル中における複数の標的核酸含有量ｇを網羅的に推定することができる。
【００３６】
図５は、図４に示した、サンプル中における標的核酸含有量ｇ（遺伝子発現量など）の推定方法において、ハイブリダイゼーション測定値ｐの誤差によって、標的核酸含有量ｇの推定値が負数になった場合の取り扱いの一例を示すフロー図である。図中、「Ｓ」はフローの始点（スタート）を、「Ｅ」はフローの終点（エンド）を示している（以下同じ）。なお、この方法に係るステップは、プログラムで記述することにより自動化できる。
【００３７】
ハイブリダイゼーション測定値ｐからサンプル中における標的核酸含有量ｇ（遺伝子発現量など）の推定する場合、前記の通り、特異値分解（Ａ＝ＵＤｔＶ）で、結合強度行列Ａの一般逆行列Ａ−１を取得し（符号５１）、サンプル中における各標的核酸含有量ｇを推定する（符号５２）。
【００３８】
符号５２のステップにおいて、標的核酸含有量ｇの最小値ｇｉが負数の場合（符号５３）、符号Ｔ（Ｔｒｕｅ）へ進む。そして、まず、標的核酸含有量ｇから最小値ｇｉを省き（符号５４）、結合強度行列Ａからｉ列の要素を省いた上で（符号５５）、スケールｓの再計算を行う（符号５６）。次に、最初のステップ（符号５１）に戻り、再計算により得られたスケールｓを用いて、再度、サンプル中における各標的核酸含有量ｇを求める（符号５２）。
【００３９】
前記処理は、標的核酸含有量ｇの最小値ｇｉが正数になるまで繰り返す。そして、最小値ｇｉが正数になった時点で、符号Ｆ（Ｆａｌｓｅ）へ進み、上記処理を終了する。標的核酸含有量ｇは、最小値ｇｉが正数になった時点における数値を採用する。
【００４０】
続いて、図６から図９を用いて、結合強度の取得方法について、説明する。結合強度の取得方法として、例えば、（１）ＢＬＡＳＴを用いた方法、（２）標的核酸の塩基配列上に検出用核酸の塩基配列をシフトさせる方法（図６、図７）、（３）測定により結合強度を取得する方法（図８、図９）がある。以下、順に説明する。
【００４１】
まず、（１）ＢＬＡＳＴを用いた方法について説明する。前記の通り、結合強度は、ある特定の検出用核酸と複数の標的核酸との結合性の度合いをそれぞれ数値化したものである。核酸同士の結合は塩基間の相補結合と相関するため、結合強度は、検出用核酸の塩基配列と標的核酸の塩基配列とが相補的になっている割合に比例すると推定される。そこで、例えば、ＢＬＡＳＴなどの公開された塩基配列検索システムを用いて、各検出用核酸の塩基配列と複数の標的核酸の塩基配列との相補性の割合をそれぞれ算出することにより、結合強度を取得することができる。
【００４２】
次に、（２）標的核酸の塩基配列上に検出用核酸の塩基配列をシフトさせる方法について、図６及び図７を用いて説明する。この方法も、前記と同様、結合強度が、検出用核酸の塩基配列と標的核酸の塩基配列とが相補的になっている割合に比例するという推定に基づく方法である。なお、この方法に係るステップは、プログラムで記述することにより自動化できる。
【００４３】
図６は、標的核酸の塩基配列上に検出用核酸の塩基配列をシフトさせる様子を模式的に示した図である。まず、検出用核酸６１の塩基配列と、検出用核酸６１と向かい合う標的核酸６２上の領域６３の塩基配列との相補性を探索し、塩基配列が相補的になっていた数をカウントする。次に、検出用核酸６１を矢印６４方向に一塩基分移動（シフト）させ、また、検出用核酸６１の塩基配列と、シフト後の標的核酸６２上の領域６３の塩基配列が相補的になっていた数をカウントする。順次矢印６４方向にシフトさせながら、両塩基配列の相補性を探索し、塩基配列が相補的になっていた数を全てカウントした後、全カウントを合計して、結合強度の指標とする。
【００４４】
図７は、図６に示した方法を数式化した図である。図７に示すとおり、検出用核酸上の塩基と標的核酸上の塩基の組み合わせが相補的であったとき（ＡとＴ、ＧとＣの組み合わせ）、塩基対の組み合わせに基づく値Ｅが１カウントされ、それ以外の場合はカウントされない。そして、変数ｓにより標的核酸上をシフトさせ、相補的な組み合わせをカウントし、全カウントを合計して、結合強度ａｉｊを得る。なお、Ｗは、長さに依存した重み値である。
【００４５】
次に、（３）測定により結合強度を取得する方法について、図８および図９を用いて説明する。
【００４６】
図８は、特定の標的核酸ｇｊを用いて結合強度を取得する方法の一例を示すフロー図である。
【００４７】
本フローは、特定の標的核酸ｇｊを合成し、ｎ回分に分注する段階（符号８１、８２）、各検出用核酸と特定の標的核酸ｇｊとのハイブリダイゼーションを測定する段階（符号８３）、各検出用核酸に係るハイブリダイゼーション測定値（蛍光強度）を、全ハイブリダイゼーション測定値の総和で規格化する段階（符号８４）、符号８３及び符号８４の段階をｎ回繰り返し、得られた規格化後の各検出用核酸に係るハイブリダイゼーション測定値の平均値を、それぞれの結合強度ａｉｊ・ｓｊとする段階（符号８５）からなる。
【００４８】
図９は、標的核酸ｇ１〜ｇｎを用いて結合強度行列を取得する方法を示した模式図である。
【００４９】
まず、標的核酸ｇ１〜ｇｎ（符号９１）について、それぞれ、各検出用核酸と標的核酸とのハイブリダイゼーションを測定し（符号９２）、標的核酸ｇ１〜ｇｎごとの結合強度ａｉｎを取得する。次に、それぞれの結合強度について、全ハイブリダイゼーション測定値の総和で規格化し（符号９３）、規格化後の結合強度ｓａｉｎを取得する（符号９４）。そして、取得した規格化後の結合強度ｓａｉｎを、結合強度行列（図４中のｓＡ）とする。
【００５０】
続いて、図１０から図１３を用いて、ノイズがある場合における標的核酸含有量の推定方法について説明する。
【００５１】
図１０は、ノイズがある場合のハイブリダイゼーション測定値とサンプル中における核酸の含有量との関係を模式的に示した図である。
【００５２】
例えば、サンプル核酸中に標的核酸１０１、１０２以外の核酸１０３が含有している場合、いずれかの検出用核酸１０４と標的核酸以外の核酸１０３とが一部（又は全長にわたって）ハイブリダイゼーションし、ノイズが発生する場合がある。その場合における、検出用核酸１０４に係るハイブリダイゼーション測定値は「数５」に示す式で表される。
【数５】
2006191844.tif 000007
【００５３】
なお、図１０において、結合強度ａ２２は検出用核酸１０５と標的核酸１０２との結合強度を表している。
【００５４】
図１１は、ノイズがある場合に、ハイブリダイゼーション測定値を、サンプル中における各標的核酸及び標的核酸以外の核酸の含有量に結合強度を乗じたものの総和として表すことができることを示した模式図である。
【００５５】
例えば、ＤＮＡチップの基板表面に固定されたｎ個の検出用核酸のうち、ある特定の検出用核酸ｉの塩基配列と、ノイズとなる核酸ｃｉの塩基配列との、相補性の度合いを結合強度ａｉ０とした場合、ある特定の検出用核酸ｉに係るハイブリダイゼーション測定値Ｐｉは、「数６」で示す式で表される。
【数６】
2006191844.tif 000008
【００５６】
図１２及び図１３は、ノイズがある場合における、本発明に係るサンプル中における標的核酸含有量（遺伝子発現量など）の推定方法の一例のフローを示す図である。なお、本フローに係わるステップは、プログラムで記述することにより自動化できる。
【００５７】
前記の通り、サンプル中にノイズとなる核酸が含有している場合における、ある特定の検出用核酸ｉに係るハイブリダイゼーション測定値Ｐｉは、図１２中の符号１２１に示す式で表すことができる。
【００５８】
符号１２１に示す式は、図４符号４１に示した式と異なり、サンプル中におけるノイズとなる核酸の含有量ｃｉが含まれている。そこで、まず、適当な目的関数ｆ(g)を導入し、ノイズとなる核酸の含有量ｃｉをベクトルとした場合に長さが最小となるｇｊを、サンプル中における各標的核酸含有量と推定する（符号１２２）。次に、目的関数ｆ(g)が最小値になるように、偏微分が０になるｇｊを求め（符号１２３）、これを符号１２１に示す式に代入してノイズとなる核酸の含有量ｃｉを消去し（符号１２４）、符号１２５に示す式を得る。
【００５９】
そして、符号１２５に示す式を行列形式で表し（符号１２６）、両辺から定数ｔＡを消去すると、符号１２７に示す行列形式を得ることができる。符号１２７に示す式は、図４中の符号４２に示した式と同じである。従って、ノイズがある場合も、前記と同様、サンプル中における各標的核酸含有量ｇを推定することができる。
【００６０】
続いて、図１４から図１７を用いて、標的核酸同士がハイブリダイゼーションする場合の、本発明に係るサンプル中における標的核酸含有量の推定方法について説明する。
【００６１】
図１４は、標的核酸同士がハイブリダイゼーションしている場合の、ハイブリダイゼーション測定値とサンプル中における核酸の含有量との関係を模式的に示した図である。
【００６２】
図１４では、図１０に加えて、標的核酸１４１と標的核酸１４２がハイブリダイゼーションしている。これは、例えば、標的核酸１４１の塩基配列の一部又は全長と、標的核酸１４２の塩基配列の一部が相補的になっている場合に起こる可能性がある。なお、図１４中、ｂ１２は、標的核酸１４１における、標的核酸１４１と標的核酸１４２のハイブリダイゼーションに対する寄与率を、ｂ２１は、標的核酸１４２における、標的核酸１４１と標的核酸１４２のハイブリダイゼーションに対する寄与率を、それぞれ示している。
【００６３】
図１５は、標的核酸同士がハイブリダイゼーションする場合における、標的核酸含有量の推定方法のフロー図である。なお、この方法に係るステップは、プログラムで記述することにより自動化できる。
【００６４】
本フローでは、まず、予め、結合強度ａｉｊとハイブリダイゼーションに対する寄与率ｂｉｊを取得しておく（符号１５１）。次に、ハイブリダイゼーション測定値Ｐｉ、結合強度ａｉｊ、ハイブリダイゼーションに対する寄与率ｂｉｊを用いて、標的核酸含有量を推定する（符号１５２、詳細は後述する）。次に、計算打ち切り指標Ｄを０に設定しておく（符号１５３）。
【００６５】
次に、ハイブリダイゼーションに対する寄与率について、操作を行う（符号１５４）。具体的には、まず、標的核酸同士のハイブリダイゼーションの含有量（ｂｉｊ・ｇｊ、ａｊｉ・ｇｉ）の差を求め（符号１５５）、その大小でハイブリダイゼーションに対する寄与率の操作を行う（符号１５６、符号１５７）。そして、前記標的核酸同士のハイブリダイゼーションの含有量の差を計算打ち切り指標Ｄに代入し、１≦ｉ、ｊ≦ｎ、ｉ≠ｊの条件下で、繰り返し演算を行う（符号１５８、符号１５９）。なお、ハイブリダイゼーションに対する寄与率についての変更重みを表すｈは０以上１未満の数値とし、予め設定しておく。
【００６６】
そして、計算打ち切り指標Ｄが閾値よりも小さくなった場合、標的核酸含有量の推定を中止する。一方、計算打ち切り指標Ｄが閾値よりも大きい場合、再び、符号１５２に戻り、標的核酸含有量の推定を行う（符号１６０）。
【００６７】
図１６は、結合強度ａｉｊとハイブリダイゼーションに対する寄与率ｂｉｊを取得する方法（図１５中の符号１５１）の一例を示した模式図である。
【００６８】
まず、ＤＮＡチップの基板表面１６１に、検出用核酸を固定した領域１６２と標的核酸を固定した領域１６３を設ける。次に、基板表面１６１に特定の標的核酸１６４を滴下又は供給し、ハイブリダイゼーションを測定する。そして、検出用核酸を固定した領域１６２におけるハイブリダイゼーション測定値より結合強度ａｉｊを、標的核酸を固定した領域１６２よりハイブリダイゼーションに対する寄与率ｂｉｊを、それぞれ取得する。さらに、全ハイブリダイゼーション測定値の総和を用いて結合強度ａｉｊとハイブリダイゼーションに対する寄与率ｂｉｊを規格化する（符号１６５）。
【００６９】
なお、結合強度ａｉｊとハイブリダイゼーションに対する寄与率ｂｉｊの取得は、図８及び図９に示した方法を用いてもよい。
【００７０】
図１７は、標的核酸含有量を推定する方法（図１５中の符号１５２）の一例のフローを示す図である。なお、この方法に係るステップは、プログラムで記述することにより自動化できる。
【００７１】
本フローでは、符号１７１に示す式の通り、結合強度ａｉｊを取得する際に、ハイブリダイゼーションに対する寄与率ｂｉｊを考慮する。その他のフローについては、図１２と同様のフローにより、標的核酸含有量を推定する。
【００７２】
図１８は、遺伝子発現量規格化の全フローの一例を示した図である。図１８中、フローＡはＲＮＡ処理フローを、フローＢはゲノム処理フローを示している。
【００７３】
ＲＮＡ処理フローＡは、本フローで用いるＤＮＡチップを準備する段階（符号Ａ１、Ａ２）、採取した試料からサンプル核酸を取得し調製する段階（符号Ａ３、Ａ４）、ＤＮＡチップの基板表面に固定された検出用核酸と採取した試料から取得したサンプル核酸中の標的核酸とのハイブリダイゼーション量を蛍光強度によって測定し、遺伝子発現量を取得する段階（符号Ａ５、Ａ６）、測定された遺伝子発現量を規格化するための指標を取得する段階（符号Ａ７）、から構成される。以下、順に説明する。
【００７４】
まず、ＤＮＡチップを準備する段階（符号Ａ１、Ａ２）について説明する。ＲＮＡ処理フローＡで用いるＤＮＡチップには、指標取得に用いる複数の検出用核酸と、遺伝子発現解析を行う検出用核酸とを固定しておく。なお、指標取得に用いる複数の検出用核酸の固定位置は任意であり、例えば、指標取得に用いる複数の検出用核酸を、基板表面の所定位置に集めて固定してもよい。
【００７５】
次に、採取した試料からサンプル核酸を取得し調製する段階（符号Ａ３、Ａ４）について、説明する。ＲＮＡ処理フローＡでは、公知の方法に従い、採取した試料からＲＮＡを抽出後、そのＲＮＡと相補的な配列を持つｃＤＮＡを合成するなどして、サンプル核酸を取得する（符号Ａ３）。サンプル核酸は制限酵素で断片化してもよい（符号Ａ４）。
【００７６】
次に、ＤＮＡチップの基板表面に固定された検出用核酸と採取した試料から取得したサンプル核酸中の標的核酸とのハイブリダイゼーション量を蛍光強度によって測定し、遺伝子発現量を取得する段階（符号Ａ５、Ａ６）について説明する。ＤＮＡチップの基板表面に固定された検出用核酸にサンプル核酸を供給し、検出用核酸とサンプル核酸中の標的核酸とのハイブリダイゼーション量を、蛍光強度などを用いて測定する（符号Ａ５）。そして、その測定データに基づいて、遺伝子発現量（規格化前の推定量）を取得する（符号Ａ６）。なお、本発明に係るサンプル中における標的核酸含有量推定方法は、この段階で適用する。
【００７７】
次に、測定された遺伝子発現量を規格化するための指標を取得する段階（符号Ａ７）について説明する。符号Ａ７の段階では、指標取得のために測定された複数の遺伝子発現量（符号Ａ６）間の相関関係を取得する。そして、取得した相関関係を、遺伝子解析のために測定された遺伝子発現量を規格化するための指標にする。このステップは、プログラムで記述することにより自動化できる。ここで、相関関係とは、前記指標取得のために測定された複数の遺伝子発現量をパラメータとする相関関数から得られた値である。相関関数は、例えば、二以上の実験条件下でそれぞれ採取された細胞から実験条件ごとに取得された複数の遺伝子発現量について、前記実験条件ごとの複数の遺伝子発現量間の相関関係を関数値とし、前記各関数値が一定値に近似する組み合わせを選択することにより得ることができる。
【００７８】
そして、符号Ａ７の段階で得られた指標を用いて（矢印Ａ８）、遺伝子発現解析を行う複数の検出用核酸とサンプル核酸中の標的核酸とのハイブリダイゼーション量（測定データ）に基づく遺伝子発現量（矢印Ａ９）を規格化する（符号Ｃ１）。このステップも、プログラムで記述することにより自動化できる。
【００７９】
ゲノム処理フローＢは、本フローで用いるＤＮＡチップを準備する段階（符号Ｂ１）、採取した試料からサンプル核酸を取得し調製する段階（符号Ｂ３、Ｂ４）、細胞数取得のための検出用核酸と採取した試料から取得したサンプル核酸中の標的核酸とのハイブリダイゼーション量を蛍光強度によって測定し、試料中の細胞数を取得する段階（符号Ｂ５、Ｂ６）、から構成される。以下、順に説明する。
【００８０】
まず、ＤＮＡチップを準備する段階（符号Ｂ１）について説明する。ゲノム処理フローＢで用いるＤＮＡチップの基板表面には、細胞数取得のための検出用核酸を固定しておく。細胞数取得のための検出用核酸には、ゲノム中にほぼ一定の割合で存在する反復配列（例えば、Ａｌｕ配列又はその一部分と同じ配列）をコードする核酸を固定する。
【００８１】
次に、採取した試料からサンプル核酸を取得し調製する段階（符号Ｂ３、Ｂ４）について説明する。ゲノム処理フローＢでは、公知の方法に従い、採取した試料からゲノムＤＮＡを抽出し、サンプル核酸を取得する（符号Ｂ３）。ゲノムＤＮＡから抽出したサンプル核酸は、制限酵素で断片化して用いる（符号Ｂ４）。
【００８２】
次に、細胞数取得のための検出用核酸と採取した試料から取得したサンプル核酸中の標的核酸とのハイブリダイゼーション量を蛍光強度によって測定し、試料中の細胞数を取得する段階（符号Ｂ５、Ｂ６）について説明する。ＤＮＡチップの基板表面に固定された検出用核酸にサンプル核酸を供給し、検出用核酸とサンプル核酸中の標的核酸とのハイブリダイゼーション量を、蛍光強度などを用いて測定する（符号Ｂ５）。そして、標的核酸中に存在する反復配列を、蛍光強度などを用いて定量的に測定することにより、試料中に含まれる細胞数を取得する（符号Ｂ６）。本発明に係るサンプル中における標的核酸含有量推定方法は、この段階で適用する。
【００８３】
そして、遺伝子発現解析を行う複数の検出用核酸と標的核酸とのハイブリダイゼーション量（測定データ）に基づく遺伝子発現量（矢印Ａ９）を、採取した試料中に含まれる細胞数を指標として、単位細胞数当たりの値に換算することにより（矢印Ｂ８）、遺伝子発現量を規格化する（符号Ｃ１）。なお、このステップは、プログラムで記述することにより自動化できる。
【００８４】
その他、符号Ａ７の段階で取得した指標に基づいて規格化した遺伝子発現量と、符号Ｂ６の段階で取得した指標に基づいて規格化した遺伝子発現量と、を比較・検討することにより、測定データ検証を行うことができる（符号Ｃ１）。このステップも、プログラムで記述することにより自動化できる。
【００８５】
続いて、図１９を用いて、本発明に係る標的核酸含有量推定システムの一例について説明する。
【００８６】
本発明に係る標的核酸含有量推定システムは、ハイブリダイゼーション測定値、結合強度などを入力する入力手段１９１、関数を出力する出力手段１９２、一つの検出用核酸と複数の標的核酸とのそれぞれの結合強度を用いて、前記ハイブリダイゼーション測定値における、標的核酸含有量の寄与の度合いを取得することにより、サンプル中における標的核酸含有量を推定する標的核酸含有量推定手段１９３、ＣＰＵ１９４、ＲＡＭ１９５、ＲＯＭ１９６を備える。
【００８７】
また、結合強度に係る数値を、ある一つの検出用核酸の塩基配列と、サンプル中の各標的核酸の塩基配列との相補性の度合いを探索することにより得る結合強度取得手段（図示せず）を備える構成にしてもよい。
【産業上の利用可能性】
【００８８】
本発明により、ＤＮＡチップなどを用いて得られたハイブリダイゼーション測定値から、サンプル中の標的核酸含有量（遺伝子発現量など）に関する高精度な数値を得ることができる。
【００８９】
本発明に係る方法、プログラム、及び、システムは、ＤＮＡチップなどを用いた測定装置に容易に組み込むことができる。
【図面の簡単な説明】
【００９０】
【図１】ＤＮＡチップを用いたハイブリダイゼーション測定手順の一例を模式的に示した図。
【図２】ハイブリダイゼーション測定値とサンプル中における標的核酸含有量との関係を模式的に示した図。
【図３】ハイブリダイゼーション測定値を、サンプル中における各標的核酸含有量に結合強度を乗じたものの総和として表すことができることを示した模式図。
【図４】本発明に係るサンプル中における標的核酸含有量（遺伝子発現量など）の推定方法の一例のフローを示す図。
【図５】図４に示した標的核酸含有量ｇの推定方法において、ハイブリダイゼーション測定値ｐの誤差によって、標的核酸含有量ｇの推定値が負数になった場合の取り扱いの一例を示すフロー図。
【図６】標的核酸の塩基配列上に検出用核酸の塩基配列をシフトさせる様子を模式的に示した図。
【図７】図６に示した方法を数式化した図。
【図８】特定の標的核酸ｇｊを用いて結合強度を取得する方法の一例を示すフロー図。
【図９】標的核酸ｇ１〜ｇｎを用いて結合強度行列を取得する方法を示した模式図。
【図１０】ノイズがある場合のハイブリダイゼーション測定値とサンプル中における核酸の含有量との関係を模式的に示した図。
【図１１】ノイズがある場合に、ハイブリダイゼーション測定値を、サンプル中における各標的核酸及び標的核酸以外の核酸の含有量に結合強度を乗じたものの総和として表すことができることを示した模式図。
【図１２】ノイズがある場合における、本発明に係るサンプル中における標的核酸含有量の推定方法の一例のフローを示す図。
【図１３】図１２の続きのフローを示す図。
【図１４】標的核酸同士がハイブリダイゼーションしている場合の、ハイブリダイゼーション測定値とサンプル中における核酸の含有量との関係を模式的に示した図。
【図１５】標的核酸同士がハイブリダイゼーションする場合における、標的核酸含有量の推定方法のフロー図。
【図１６】結合強度ａｉｊとハイブリダイゼーションに対する寄与率ｂｉｊを取得する方法の一例を示した模式図。
【図１７】標的核酸含有量を推定する方法（図１５中の符号１５２）の一例のフローを示す図。
【図１８】遺伝子発現量規格化の全フローの一例を示した図。
【図１９】標的核酸含有量推定システムの一例を示す図。
【図１】
2006191844.tif 000009
【図２】
2006191844.tif 000010
【図３】
2006191844.tif 000011
【図４】
2006191844.tif 000012
【図５】
2006191844.tif 000013
【図６】
2006191844.tif 000014
【図７】
2006191844.tif 000015
【図８】
2006191844.tif 000016
【図９】
2006191844.tif 000017
【図１０】
2006191844.tif 000018
【図１１】
2006191844.tif 000019
【図１２】
2006191844.tif 000020
【図１３】
2006191844.tif 000021
【図１４】
2006191844.tif 000022
【図１５】
2006191844.tif 000023
【図１６】
2006191844.tif 000024
【図１７】
2006191844.tif 000025
【図１８】
2006191844.tif 000026
【図１９】
2006191844.tif 000027