説明

塩基配列判定方法、塩基配列判定システム及び塩基配列判定プログラム

【課題】 判定精度が高く実態に即した判定ができる塩基配列判定方法、塩基配列判定システム及びプログラムを提供する。
【解決手段】 第1検出用電極、第2検出用電極、コントロール電極を備えるチップカートリッジ11、第1検出用電極から第1検出信号、第2検出用電極から第2検出信号、コントロール電極からコントロール信号を取得する測定系12、チップカートリッジ11に薬液を送液する送液系13、測定系12及び送液系13を制御する制御機構15、第1検出信号の平均値からコントロール信号の平均値を減算し、コントロール信号の平均値で除算した値をX座標、第2検出信号の平均値からコントロール信号の平均値を減算し、コントロール信号の平均値で除算した値をY座標とし、原点からX座標及びY座標が定める点までのベクトルが正のX軸となす角度を角度判定基準と比較しDNAの型を判定する型判定モジュールを含むコンピュータ16を備える。

【発明の詳細な説明】
【技術分野】
【0001】
本発明は、核酸の塩基配列を判定する塩基配列判定方法、この塩基配列判定方法を実現する塩基配列判定システム、及びこの塩基配列判定システムを自動で制御し、測定信号を自動解析される塩基配列判定プログラムに関する。
【背景技術】
【0002】
ヒトゲノムは、約30億の遺伝子暗号(塩基)により構成されている。現在進行している「ヒトゲノムプロジェクト」は、この遺伝子暗号(塩基配列)をすべて解き明かすことになる。これらの経緯の中で、個人間において多くの遺伝子暗号(塩基配列)の違いが存在することが明らかになってきている。ヒトゲノムの塩基配列の違い(多型)は、1個の塩基が他の塩基に置き換わっている一塩基多型(SNP)、1塩基から数千塩基に亘る欠損や挿入によるもの、繰り返し配列の繰り返し回数の違い(VNTR,マイクロサテライト多型)などに分類されているが、現在は、これらの多型の中で、特に一塩基多型(SNP)が注目を集めている。一塩基多型(SNP)とは、DNAの塩基配列の中でたった1つの塩基が違うことであり、酒に強い弱いや、薬が効きやすいかどうかなど、人間の個性の最小単位にあたるものである。ヒトゲノム30億塩基対の内、一塩基多型(SNP)は約300万個(500〜1000塩基対に1個の割合)〜1,000万個存在すると示唆され、特定のタンパク質が作られなかったり、他人と違うものを作ったりして個人差(体質)、人種差などの違いをもたらしている。ヒトにおける遺伝子の個人差の研究では、一塩基多型を解析して、病気に対する感受性や薬物への応答を調べ、そのヒトにあった副作用の少ない薬を投薬するようなテーラーメイド医療が可能になると言われ、一塩基多型(SNP)の解析の研究が進められている。植物であれば、植物が本来持っている病気や害虫に対する抵抗性の仕組みを解明し、その機能を高めたりすることができる。
【0003】
一塩基多型(SNP)が注目を集める理由としては、その解析技術向上により多数のSNPの解析が可能になったことによる疾患とSNPとの関連についての関心の高まりが上げられる。その研究目的は、疾患関連遺伝子や、薬剤代謝の個人差の解析、慢性疾患など広範囲に亘っている。一部の薬物代謝、脂質代謝については、SNPとの関連も明らかになってきている。今後も、SNPとこれらの関連が徐々に明らかになってくるものと期待されている。
【0004】
SNPの解析のような分子生物学的技術には、非常に多くの試料につき膨大な数の操作を行うことが含まれている。それらは、しばしば、複雑で時間が掛かり、一般に高い精度が要求される。多くの技術は、感度、特異性又は再現性の欠如によってその適用が制限されている。
【0005】
例えば、感度及び特異性に伴う問題は、これまで、核酸ハイブリダイゼーションの実際的な適用を制限してきた。「ハイブリダイゼーション」は核酸のハイブリッド形成又は核酸雑種分子形成のことを言い、核酸の1次構造、即ち塩基配列の相同性を調べたり、相同の塩基配列を持つ核酸を検出したりする方法として用いられている。1本鎖にした核酸同士が、相補性を持つ塩基対間、即ち、アデニン(A)−チミン(T)間、或いはグアニン(G)−シトシン(C)間で水素結合ができ、二重らせんの2本鎖核酸を形成する性質を利用するものである。核酸ハイブリダイゼーshン分析には、一般的に、大量の被標的核酸の中から非常に少数の特異的標的核酸(DNA又はRNA)をプローブで検出することが含まれる。高い特性を維持するためには、通常は、温度、塩、洗浄剤、溶媒、カオトロピック剤及び変性剤の種々の組合せにより達成される最も厳格な条件下でハイブリダイゼーションを行う。大部分の試料、特にヒトゲノムDNA試料中のDNAが非常に複雑なことに関連する。試料が、特異的標的配列に酷似した膨大な数の配列よりなる場合には、最もユニークなプローブ配列でさえ、被標的配列と多数の部分的ハイブリダイゼーションを起こす。又、プローブDNAとの特異的標的(サンプルDNA)との間の好ましくないハイブリダイゼーション動力学が関連する場合もある。最良の条件下でさえ、大部分のハイブリダイゼーション反応は、相対的に低濃度のプローブDNA及び標的分子(サンプルDNA)で行われる。加えて、プローブDNAはしばしばサンプルDNAに対する相補鎖と競合する。又、高レベルの非特異的バックグラウンド・シグナルが、ほとんどいずれの物質にもプローブDNAが親和性を有することによって引き起こされるという問題もある。したがって、これらの問題が個々に、又は組み合わさって、核酸ハイブリダイゼーションの感度や特異性を喪失させる。
【0006】
この様な事情のもと、本発明者は、既に、例えば、SNP(一塩基多型)の判定(2種の塩基のホモ/へテロ判定)をするのに、それぞれの信号の大きさに対するt検定により、有意差判定を行う方法や(特許文献1参照)や、それぞれの信号の大きさの比を解析する方法を提案している(特許文献2参照)。特許文献1に記載された方法では、ほぼ完全に2種の信号値が一致していない限り、ヘテロ判定をすることはないが、実際の測定では、この様なことは不可能である。又、特許文献2に記載された方法では、その課題を克服すべく、信号の比によって解析する方法を提案しているが、一方の信号増加量が「負」若しくは非常に小さい場合の判定精度が悪いという問題がある。
【0007】
この様に、従来も、核酸の塩基配列を判定するアルゴリズムは存在していたが、判定の精度に課題があった。
【特許文献1】特開2004−125777号公報
【特許文献2】特開2006−258702号公報
【発明の開示】
【発明が解決しようとする課題】
【0008】
本発明は、判定の精度が高く、実態に即した判定ができる、塩基配列判定方法、塩基配列判定システム及び塩基配列判定プログラムを提供することを目的とする。
【課題を解決するための手段】
【0009】
本発明の第1の態様は、
第1プローブがそれぞれ固定された複数の第1検出用電極、前記第1プローブとは塩基配列が異なる第2プローブがそれぞれ固定された複数の第2検出用電極、前記複数の第1検出用電極及び第2検出用電極に対応して配置され、前記第1プローブ及び前記第2プローブとは塩基配列が異なるコントロールDNAがそれぞれ固定された複数のコントロール電極とを備えるチップに、サンプル核酸を含む薬液を供給する段階と、
前記複数の第1検出用電極を介してそれぞれ第1検出信号を、前記複数の第2検出用電極を介してそれぞれ第2検出信号を、前記複数のコントロール電極を介してそれぞれコントロール信号を、それぞれ取得する段階と、
前記第1検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をX座標とする段階と、
前記第2検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をY座標とする段階と、
原点から前記X座標及び前記Y座標が定める点までのベクトルが正のX軸となす角度を算出する段階と、
前記角度を、予め定めた角度判定基準と比較する段階
とを含み、前記角度と前記角度判定基準との大小関係により、前記サンプル核酸の型を判定することを特徴とする塩基配列判定方法である。
【0010】
本発明の第2の態様は、
第1プローブがそれぞれ固定された複数の第1検出用電極、前記第1プローブとは塩基配列が異なる第2プローブがそれぞれ固定された複数の第2検出用電極、前記第1プローブ及び前記第2プローブとは塩基配列が異なるコントロールDNAがそれぞれ固定された複数のコントロール電極とを備えるチップを装着して使用され、
前記複数の第1検出用電極を介してそれぞれ第1検出信号を、前記複数の第2検出用電極を介してそれぞれ第2検出信号を、前記複数のコントロール電極を介してそれぞれコントロール信号を、それぞれ取得する測定系と、
前記チップに、薬液を送液する送液系と、
前記測定系及び送液系を制御する制御機構と、
前記第1検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をX座標とし、前記第2検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をY座標とし、原点から前記X座標及び前記Y座標が定める点までのベクトルが正のX軸となす角度を算出し、前記角度を、予め定めた角度判定基準と比較し、前記角度と前記角度判定基準との大小関係により、前記サンプル核酸の型を判定する型判定モジュールを含むコンピュータ
とを備えることを特徴とする塩基配列判定システムである。
【0011】
本発明の第3の態様は、
第1プローブがそれぞれ固定された複数の第1検出用電極、前記第1プローブとは塩基配列が異なる第2プローブがそれぞれ固定された複数の第2検出用電極、前記第1プローブ及び前記第2プローブとは塩基配列が異なるコントロールDNAがそれぞれ固定された複数のコントロール電極とを備えるチップにサンプル核酸を含む薬液を供給された後、
前記複数の第1検出用電極を介してそれぞれ第1検出信号を、前記複数の第2検出用電極を介してそれぞれ第2検出信号を、前記複数のコントロール電極を介してそれぞれコントロール信号を、それぞれ取得する手順と、
前記第1検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をX座標とし、前記第2検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をY座標とし、原点から前記X座標及び前記Y座標が定める点までのベクトルが正のX軸となす角度を算出する手順と、
前記角度を、予め定めた角度判定基準と比較する手順
とを含む処理をコンピュータに実行させ、これにより前記コンピュータに、前記サンプル核酸の型が、野生型、ヘテロ型、変異型、若しくは判定不能と判断させるための塩基配列判定プログラムである。
【発明の効果】
【0012】
本発明によれば、SNPのホモ・ヘテロの型判定を、高い精度で実態に即してできる塩基配列判定方法、塩基配列判定システム及び塩基配列判定プログラムを提供することができる。
【発明を実施するための最良の形態】
【0013】
次に、図面を参照して、本発明の実施の形態を説明する。以下の図面の記載において、同一又は類似の部分には同一又は類似の符号を付している。ただし、図面は模式的なものであり、厚みと平面寸法との関係、各層の厚みの比率等は現実のものとは異なることに留意すべきである。したがって、具体的な厚みや寸法は以下の説明を参酌して判断すべきものである。又、図面相互間においても互いの寸法の関係や比率が異なる部分が含まれていることは勿論である。又、以下に示す実施の形態は、本発明の技術的思想を具体化するための装置や方法を例示するものであって、本発明の技術的思想は、構成部品の材質、形状、構造、配置等を下記のものに特定するものでない。本発明の技術的思想は、特許請求の範囲に記載された技術的範囲内において、種々の変更を加えることができる。
(塩基配列判定システム)
本発明の実施の形態に係る塩基配列判定システムは、図7に示すように、チップカートリッジ11と、このチップカートリッジ11と電気的に接続される測定系12、チップカートリッジ11に設けられた流路とインタフェース部を介して物理的に接続される送液系13及びチップカートリッジ11の温度制御を行う温度制御機構14等を備える。図7の測定系12は、図1に示すように、対極502の入力に対して参照極561,562の電圧をフィードバック(負帰還)させることにより、セル内の電極や溶液などの各種条件の変動によらずに溶液中に所望の電圧を印加する3電極方式のポテンシオ・スタットであり、検出チップ21の端子C,R,Wに接続されている。図1の検出チップ21は、図7のチップカートリッジ11に備えられている。
検出チップ21は、図2に示すように、SNP1検出用電極である作用極551,SNP2検出用電極である作用極552,コントロール電極である作用極553と、これらの作用極551,552,553に対する参照極561,562及び対極502が検出チップ上に配列された電極ユニットを用いている。即ち、測定系12は、作用極551,552,553に対する参照極561,562の電圧を、ある所定の特性に設定されるように対極502の電圧を変化させ、挿入剤の電気化学的反応による電流(以下において「電気化学的電流」と言う。)を電気化学的に測定する。
本発明の実施の形態に係る塩基配列判定方法は、先ず、塩基配列判定の目的とされる標的塩基配列(サンプルDNA)581,582,583(図9〜11参照。)と相補的な塩基配列を有するプローブDNA571,572について、それぞれ作用極551,552に固定化させる。即ち、図3(a)に示すように、作用極551は、図9に示した塩基配列CTGAG……を有する標的塩基配列(サンプルDNA)581とは相補的な塩基配列GACTC……を有するプローブDNA571が固定化される電極である。なお、本実施形態では、サンプルDNA,プローブDNAと記載しているが、DNA以外にPNA,RNA等、他の核酸をサンプルもしくはプローブとして用いても良い。図3(a)では、SNP位置の塩基を、下から3番目のCとしたSNP="G"検出用電極である。
又、図3(b)に示すように、作用極552は、図10に示した塩基配列CTTAG……を有する標的塩基配列(サンプルDNA)582とは相補的な塩基配列GAATC……を有するプローブDNA572が固定化される電極である。図3(b)も、SNP位置の塩基を、下から3番目のAとしたSNP="T"検出用電極である。
一方、図3(c)に示すように、作用極553は、標的塩基配列(サンプルDNA)581,582とは相補的関係のない塩基配列CAGTG……を有するプローブDNA(コントロールDNA)573が固定化されるコントロール電極である。これらの作用極551,552,553は、セル内の反応電流を検出する電極として機能する。なお、作用極551,552,553に固定化されるプローブDNA571,572,573の種類は 例示に過ぎないが、作用極551,552,553に固定化されるプローブDNAの種類は、それぞれ異なる塩基配列を有するサンプルDNAを検出するように設計されていれば良く、それぞれの作用極に固定化されるプローブの種類は単一である必要はない。
【0014】
対極502は、作用極551,552,553との間に所定の電圧を印加してセル内に電流を供給する電極である。参照極561,562は、参照極561,562と作用極551,552,553との間の電圧を所定の電圧特性に制御すべく、その電極電圧を対極502にフィードバックさせる電極であり、これにより対極502による電圧が制御され、セル内の各種検出条件に左右されない精度の高い電気化学的電流の検出が行える。
【0015】
図1に示すように、本発明の実施の形態に係る塩基配列判定システムの測定系12は、電極間を流れる電流を検出するための電圧パターンを発生させる電圧パターン発生回路510を備え、この電圧パターン発生回路510は、配線512bを介して参照極561,562の参照電圧制御用の反転増幅器512(OP)の反転入力端子に接続されている。電圧パターン発生回路510は、図7に示す制御機構15から入力されるデジタル信号をアナログ信号に変換して電圧パターンを発生させる回路であり、DA変換器を備える。配線512bには抵抗Rが接続されている。反転増幅器512の非反転入力端子は接地され、出力端子には配線502aが接続されている。反転増幅器512の反転入力端子側の配線512bと出力端子側の配線502aは配線512aで接続されている。この配線512aには、フィードバック抵抗Rff及びスイッチSWfからなる保護回路500が設けられている。線502aは検出チップ21の端子Cに接続されている。端子Cは、検出チップ21上の対極502に接続されている。対極502が複数設けられている場合には、各々に対して並列に端子Cが接続される。これにより、1つの電圧パターンにより複数の対極502に同時に電圧を印加することができる。線502aには、端子Cへの電圧印加のオンオフ制御を行うスイッチSWが設けられている。
【0016】
反転増幅器512に設けられた保護回路500により、対極502に過剰な電圧がかからないような構成となっている。したがって、測定時に過剰な電圧が印加され、溶液が電気分解されてしまうことにより、所望の挿入剤の電気化学的電流の検出に影響を及ぼすことがなく、安定した測定が可能となる。
【0017】
検出チップ21の端子Rは、配線503aにより電圧フォロア増幅器(OPr)513の非反転入力端子に接続されている。電圧フォロア増幅器513の反転入力端子は、その出力端子に接続された配線513bと配線513aにより短絡している。配線513bには抵抗Rfが設けられており、配線512bの抵抗と、配線512aと配線512bの交点との間に接続されている。これにより、電圧パターン発生回路510により生成される電圧パターンに、参照極561,562の電圧をフィードバックさせた電圧を反転増幅器512に入力させ、そのような電圧を反転増幅した出力に基づき対極502の電圧を制御する。
【0018】
検出チップ21の端子Wは、配線501aによりトランス・インピーダンス増幅器(OPw)511の反転入力端子に接続されている。トランス・インピーダンス増幅器511の非反転入力端子は接地され、その出力端子に接続された配線511bと配線501aとは配線511aにより接続されている。配線511aには抵抗Rwが設けられている。このトランス・インピーダンス増幅器511の出力側の端子Oの電圧をVw、電流をIwとすると:
w=Iw・Rw ……(1)
となる。この端子Oから得られる電気化学信号は、図7に示す制御機構15に出力される。作用極551,552,553は複数あるため、端子W及び端子Oは、作用極551,552,553のそれぞれに対応して複数設けられる。複数の端子Oからの出力は後述する信号切替部により切り替えられ、AD変換されることにより各作用極551,552,553からの電気化学信号をデジタル値としてほぼ同時に取得することができる。なお、端子W及び端子Oの間のトランス・インピーダンス増幅器511などの回路は、複数の作用極551,552,553で共有しても良い。この場合、配線501aに複数の端子Wからの配線を切り替えるための信号切替部を備えれば良い。
図7の塩基配列判定システムを構成するチップカートリッジ11は、図4及び5に示すように、カセット上蓋711、カセット下蓋712、パッキン713(シール部材)及び基板714からなるカセットを構成している。カセット上蓋711及びカセット下蓋712の内表面同士を対向させ、かつこれらカセット上蓋711及びカセット下蓋712の間にパッキン713及び基板714を狭んだ状態で固定している。カセット上蓋711の外表面から内表面にかけて、断面が略円形のノズル差込孔722及び723が貫通して形成されている。このノズル差込孔722及び723の内径は、送入及び送出ポート752,753、及び図6のバルブユニット705のノズル707,708、の外径よりも若干大きく設定されており、例えば3.2mm程度である。図4に示すように、カセット上蓋711の外表面から内表面にかけて断面が略長方形の電気コネクタ用ポート724及び725が貫通して形成されている。この電気コネクタ用ポート724及び725の各々には、後述する電気コネクタが挿着され使用される。更に、外表面には内表面にかけてシール検出孔726が貫通して形成されている。このシール検出孔726は、シールの有無の検出に用いられる。即ち、カセット(検出チップ)21の外表面のシール検出孔726の表面から電気コネクタ用ポート724及び725の表面にかけてシールが貼付された状態で、カセット(検出チップ)21へ薬液(試料)を注入し、カセット(検出チップ)21に薬液(試料)を注入した後にシールが剥がすが、そのシールの有無の検出する。カセット(検出チップ)21に、シールが貼付された状態で、薬液(試料)を注入することにより、たとえ薬液(試料)溶液が誤って電気コネクタポート724若しくは725上に滴下しても、シールで被覆されているため、実際のポート724若しくは725内には液が侵入することがなく、電気的に短絡してしまう等の不具合を発生させる心配を生じないようにしている。
図5に示すように、カセット上蓋711の内表面側には、所定の深さでかつ基板714の断面形状とほぼ同じ断面形状を有する基板位置決め溝が設けられており、この基板位置決め溝の周囲は内表面で囲まれている。基板位置決め溝はノズル差込孔722及び723、電気コネクタ用ポート724及び725にオーバーラップして形成されている。基板位置決め溝にあわせて基板714をはめ込むことにより、基板714をカセット上蓋711に位置決め配置することができる。基板位置決め溝の深さは基板714の厚さとほぼ同じになるように形成されている。
【0019】
図5に示すように、カセット上蓋711の内表面側で基板位置決め溝にオーバーラップして、基板位置決め溝よりも更に深いパッキン位置決め溝が設けられており、このパッキン位置決め溝の周囲は基板位置決め溝で囲まれている。パッキン位置決め溝はノズル差込孔722及び723にオーバーラップして形成されている。パッキン位置決め溝にあわせてパッキン713をはめ込むことにより、パッキン713をカセット上蓋711に位置決め配置することができる。パッキン位置決め溝の基板位置決め溝に対する深さは後述するパッキン713の厚さとほぼ同じになるように形成されている。したがって、パッキン位置決め溝の内表面に対する深さは、パッキン713の厚さに基板714の厚さを加算した厚さとほぼ同じになるように形成されている。
【0020】
カセット上蓋711の内表面の周縁部には4つのねじ孔727a、727b、727c及び727dが設けられている。これらねじ孔727a〜727dにより、カセット上蓋711とカセット下蓋712をねじ止めできる。カセット上蓋711の内表面の周縁部には2つのカセット位置決め孔728a及び728bが設けられている。このカセット位置決め孔728a及び728bを、スライドステージの上に設けられた2つの位置決め用ピンにあわせてカセット(検出チップ)21を配置することにより、スライドステージに対してカセット(検出チップ)21を位置決め配置することができる。
【0021】
又、外表面にはシール検出孔746が貫通して形成されている。又、カセット下蓋712のシール検出孔746は、カセット上蓋711とカセット下蓋712が止着された状態で、カセット上蓋711のシール検出孔726と通じる位置に形成されている。これにより、カセット上蓋711とカセット下蓋712が止着された際に、カセット上蓋711からカセット下蓋712にかけて貫通したシール検出孔726が設けられ、このシール検出孔726に対して検出光を照射することにより、シールの有無を判別できる。
【0022】
カセット下蓋712の外表面の周縁部には4つのねじ孔747a、747b、747c及び747dが設けられている。これらねじ孔747a〜747dとカセット上蓋711に設けられた対応するねじ孔727a〜727dの各々をねじ止めすることにより、カセット下蓋712をカセット上蓋711に止着することができる。カセット下蓋712の外表面の周縁部には、2つのカセット位置決め孔748a及び748bが設けられている。このカセット位置決め孔748a及び748bには、カセット上蓋711のカセット位置決め孔728a及び728bがそれぞれ貫通する。スライドステージに対するカセット(検出チップ)21の位置決めは、スライドステージの上に設けられた2つの位置決め用ピンと、カセット下蓋712のカセット位置決め孔748a及び748bを貫通しているカセット上蓋711の2つの位置決め孔728a及び728bによって規定される構成になっている。カセット下蓋712には、更にカセット種類判別用孔749が設けられ、このカセット種類判別用孔749の有無により、カセット(検出チップ)21の種類を判別できる。なお、種類判別は、カセット種類判別用ピン(図示省略)の押し下げの有無により自動で行うことができる。カセット種類判別用ピン(図示省略)のピンの押し下げ状態は、制御機構15で検知される。
【0023】
なお、カセット種類判別用孔749が設けられていないカセット(検出チップ)21を用いても、判別されるカセット(検出チップ)21の種類が異なるのみで同様の測定が可能である。若しくは、カセット(検出チップ)21の種類が異なることを表示部(図示省略)に制御機構15が表示させて警告して、測定段階に進まないように設計しておくこともできる。或いは、カセット種類判別用ピン(図示省略)は固定された固定ピンを用いて、カセット種類判別用孔749が設けられていないカセット(検出チップ)21は装着できないようにすることにより、誤ったカセット(検出チップ)21をセットしないようにすることもできる。 パッキン713は、ほぼ長方形で所定の厚さを有し、その4隅が切欠いて形成された板状部分と、この板状部分の主面上であって長辺の両端近傍に位置し、かつ短辺の中央近傍に設けられた円筒形状の送入ポート752及び送出ポート753からなる。送入ポート752及び送出ポート753の先端には開口部が設けられている。この開口部から板状部分にかけて、送入ポート752及び送出ポート753の軸心には、板状部分の主面に対して垂直な方向に流路が設けられている。板状部分の裏面には、図5に示すように、送入ポート752の形成位置から送出ポート753の形成位置にかけて折れ曲がり形状の溝が形成されている。溝は、流路を構成している。溝は、複数回交互に進行するように形成され、溝の折り返し地点を所定の曲率のカーブとすることにより、折り返し地点のコーナーなどが設けられた場合に発生する薬液やエアの滞留を抑制できる。
【0024】
図4に示すように、基板714の主面には、電極ユニット761、パッド762及び763が形成されている。電極ユニット761は、図2に示すように、対極502、作用極551,552,553及び参照極561,562の組合せからなる3電極系のユニットである。電極ユニット761の作用極551,552,553にはプローブDNA571,572,573が固定化される。電極ユニット761とパッド762、電極ユニット761とパッド763は、図示しない配線により接続されている。図4では、基板714に対して、プローブDNAが固定されている電極ユニット761と、パッド762及び763が同一面上に形成させている場合を例示しているが、パッド762及び763が、基板714の電極ユニット761が形成された面に対して反対側の面に形成されている場合には、バルブユニット705はカセット(検出チップ)21に対して上側に、そしてプローブユニット710は下側に設置されることになる。その場合には、必然的に、バルブユニット705とプローブユニット710は一体的に形成される訳ではなくなる。
【0025】
電極ユニット761の配列は、パッキン配置位置でかつ折れ曲がり形状の流路形成位置にあわせて形成されている。これにより、パッキン713と基板714がカセット上蓋711及びカセット下蓋712に位置決めされた状態で狭着固定された際には、折れ曲がり形状の溝と基板714表面により折れ曲がり形状の流路が形成され、かつその流路表面には電極ユニット761が露出した構成となる。即ち、電極ユニット761上は、折れ曲がり形状の溝により間隙が設けられ、この間隙により折れ曲がり形状の流路が形成される。又、この状態では、パッキン713と基板714とはシールが保持される。
先ず、カセット上蓋711の内表面のパッキン位置決め溝にあわせて、かつノズル差込孔722及び723に送入ポート752及び送出ポート753が挿通するように、パッキン713をパッキン位置決め溝に嵌挿する。次に、基板714を、その主面、即ち電極ユニット761、パッド762及び763が形成された面がカセット上蓋711側に向くように、基板位置決め溝に位置決め配置する。次に、カセット下蓋712を、その内表面742がカセット上蓋711側に向くように、又ねじ孔747a〜747d及びねじ孔727a〜727dの位置が対応するようにカセット上蓋711上に載置する。そして、ねじ孔747a〜747d及びねじ孔727a〜727dにねじ770a〜770dを螺挿する。これにより、カセット上蓋711及びカセット下蓋712が螺着され、かつカセット上蓋711及びカセット下蓋712間にはパッキン713及び基板714が狭着固定され、カセット(検出チップ)21が完成する。この完成した状態では、ノズル差込孔722からノズル差込孔723にかけて、折れ曲がり形状の流路が形成される。
【0026】
図4及び5では、ねじ止めによりカセット上蓋711及びカセット下蓋712を固定する例を示したが、これに限定されない。例えば凹凸の部材を用いた係着手法を用いても良い。 図6は、本発明の実施の形態に係る塩基配列判定システムを構成する送液系13に備えられバルブユニット705の全体構成を示す図である。なお、図6では、プローブユニットの構成は省略して記載してあるが、バルブユニット705にはプローブユニットが一体的に形成され、バルブユニット・プローブユニット駆動機構によりこれらバルブユニット及びプローブユニットが同時に駆動される。例えばガラスエポキシ基板などからなるプローブユニットに、2つの電気コネクタが所定の間隔をおいて配置される。電気コネクタの先端には、複数の凸状電極が基板714上のパッドと同じ配列によりマトリクス状に配置されており、これら各凸状電極が図4に示した基板714上のパッド762,763と接触することにより、基板714とプローブユニットとの電気的接続が確保される。又、電気コネクタ内には配線が設けられており、各凸状電極と制御機構15がこの配線により電気的に接続される。
【0027】
バルブユニット・プローブユニット駆動機構は制御機構15からの指示により駆動する。バルブユニット・プローブユニット駆動機構は、昇降方向の駆動方向を有する。これにより、スライドステージ側のカセット(検出チップ)21の上部に対してノズル707、708及び電気コネクタ703が降下することにより、図6に示すように、ノズル707及び708がノズル差込孔722及び723に位置決めされ、かつ電気コネクタ703が電気コネクタ用ポート724及び725に位置決めされる。これにより、送液系とカセット(検出チップ)21内の流路が連通する。又、電気コネクタがカセット(検出チップ)21のパッドに位置決めされ、パッドと電気コネクタが電気的に接続される。
【0028】
バルブユニット705は、バルブボディ781が連結固定して用いられる。バルブボディ781には2方電磁弁403、3方電磁弁413、423及び433が設けられており、又バルブボディ782には3方電磁弁441及び445が設けられている。バルブボディ781は、例えばPEEK樹脂により作製される。なお、バルブボディ781が別個に作製され、その両者をつなぎ合わせる場合には、その継ぎ目部分には例えばPTFEをパッキンとして使用する。したがって、バルブボディ781の両者では、薬液に接する部分の材質はPEEK及びPTFEからなる。又、バルブボディ781にはほぼ一定の断面形状の空洞が設けられている。この空洞は、後述する各電磁弁の間やパッキン713等との間を接続する配管として機能する。又、バルブボディ782に設けられた空洞には、ノズル707及びノズル708が連通している。これらノズル707及びノズル708はPEEK樹脂からなる。
【0029】
3方電磁弁413は、エアと純水を切り換えて下流側の3方電磁弁423に供給する。3方電磁弁423は、バッファと3方電磁弁413からのエア又は純水を切り換えて下流側の3方電磁弁433に供給する。3方電磁弁433は挿入剤と3方電磁弁423からのエア、純水又はバッファを切り換えて下流側のバルブボディ782に供給する。3方電磁弁441は、バルブボディ781からのエア又は薬液のノズル707への供給又はバイパス配管を介した3方電磁弁445への供給を切り換える。3方電磁弁445は、3方電磁弁441からのエア又は薬液の供給又はカセット(検出チップ)21からのノズル708を介した薬液又はエアの送出を切り換える。
【0030】
図6に示すバルブユニット705において、カセット(検出チップ)21内にバッファを送液するためには、3方電磁弁423、441、445、及び送液ポンプ454をONにすれば良い。これにより、バッファが吸い上げられ、薬液はノズル707側に切り替わり、ノズル707からカセット(検出チップ)21へ、更にカセット(検出チップ)21からノズル708に吸い出され、3方電磁弁445を経由して廃液することができる。カセット(検出チップ)21内に純水を送液するためには、3方電磁弁423にかえて3方電磁弁413をONにすれば良い。カセット(検出チップ)21内に挿入剤を送液するためには、3方電磁弁423にかえて3方電磁弁433をONにすれば良い。カセット(検出チップ)21内にエアを供給するためには、3方電磁弁403をONにし、3方電磁弁413,423及び433のいずれもOFFにすれば良い。バルブユニット705のバルブボディ781内に設けられた空洞部分の配管の内部容量は、バルブ内容量を含めても100μL程度である。本実施の形態とは異なり各3方弁をチューブで接続して同じフローを構成する場合には、500μL程度の内部容量が必要であるが、これと比較すると大幅に試薬量を節減することができる。更に、バルブユニット705とカセット(検出チップ)21間の内部容量も、本実施の形態とは異なる例では100μL以上あるが、本実施の形態では10μLと大幅に低減できる。この様な構造により、試薬切換え後に、更に不本意にカセット(検出チップ)21内を流れる溶液若しくは空気の量を大幅に減らすことができる。その結果、反応及び測定のバラつきを低減し、結果の再現性も大幅に向上する。
【0031】
又、図示していない、薬液揺動装置を備えており、チップカセット内の薬液揺動を行うことができる。薬液の揺動は、(1)サンプルDNAとプローブDNAとのハイブリダイゼーション工程、(2)洗浄工程、(3)挿入剤供給工程などで行うのが有効である。(1)のハイブリダイゼーション工程でサンプルDNAを揺動させることにより、ハイブリダイゼーション効率を向上させ、ハイブリダイゼーション時間を短縮させることができる。又、(2)洗浄工程でバッファ液を揺動させることにより、非特異吸着DNAを引き剥がす効率を向上させることにより、洗浄時間を短縮させることができる。又、(3)挿入剤供給工程で挿入剤を揺動させることにより、挿入剤の濃度の均一性を向上させ、又挿入剤の吸着均一性も向上させることにより、信号のバラつき、S/N比を改善することができる。これら(1)〜(3)のすべての工程に薬液揺動工程を適用しても、一部の工程に適用しても、薬液揺動の効果が得られる。
冒頭で、一部言及したが、本発明の実施の形態に係る塩基配列判定システムは、図7に示すように、測定ユニット10と、測定ユニット10に接続された制御機構15と、制御機構15に接続されたコンピュータ16とを備える。測定ユニット10は、チップカートリッジ11と、このチップカートリッジ11と電気的に接続される測定系12、チップカートリッジ11に設けられた流路とインタフェース部を介して物理的に接続される送液系13及びチップカートリッジ11の温度制御を行う温度制御機構14等を備える。
図7に示したコンピュータ16は、図8に示すように、操作者からのデータや命令などの入力を受け付ける入力装置304と、SNPの型が2種の内のどちらであるか、又、それがホモ型であるかヘテロ型であるかを判定するか、もしくは核酸の存在の有無を検出する中央処理装置(CPU)300と、判定結果を出力する出力装置305及び表示装置306と、塩基配列の判定に必要な所定のデータなどを格納したデータ記憶部(図示省略)と、塩基配列判定プログラムなどを格納したプログラム記憶部(図示省略)とを備える。
CPU300は、ノイズ除去モジュール301,電流波形判定モジュール302,ピーク電流検出モジュール310,グループ正常判定モジュール320,検査有効性判定モジュール325、有無判定モジュール330及び型判定モジュール340を備える。
【0032】
ノイズ除去モジュール301は、図2に示したSNP1検出用電極551,SNP2検出用電極552,コントロール電極553を介して測定される電流から、「単純移動平均法」を用いて平準化(スムージング)により、ノイズを除去する。スムージングは、例えば、図15に示すような、単純移動平均法を用いれば良い。「単純移動平均法」は文字通り、いくつかの実績値を単純平均するものであり、図15(a)に示すような時系列データを、その規則性に着目して平均するものである。図15では、移動平均値の区間をm=5として、それぞれのデータの後5点のデータを5で割った数値:
y[n]=(x[n]+x[n+1]+x[n+2]+x[n+3]+x[n+4])/5 ……(2)
が、図15(b)に示す移動平均値となる。移動平均により、図15(a)に示すような変動が図15(b)に示すように平準化(スムージング)され、全体的な傾向を分析するが容易になる。
電流波形判定モジュール302は、図16のフローチャートに示した手順に従って、図2に示したSNP1検出用電極551,SNP2検出用電極552,コントロール電極553を介してそれぞれ測定された電流波形(電流−電圧特性)のベースラインの傾きをそれぞれ求め、それぞれのベースラインの傾きにより、それぞれの検出信号(電流波形)の正常・異常を判定し、異常な検出信号を演算対象外とする。
ピーク電流検出モジュール310は、電圧値算出部311,ベースライン近似部312,ピーク電流値演算部313を備え、図18及び図19のフローチャートに示した手順に従い、SNP1検出用電極551,SNP2検出用電極552,コントロール電極553を介して測定される電流(検出信号)から、バックグラウンド電流を差し引いて、挿入剤591に起因する真の電気化学的電流(真の検出信号)のピーク値(ピーク電流値)を求める。
【0033】
電圧値算出部311は、図18のステップS221〜S223に示す手順に従い、チップカートリッジ11で測定された電気化学的電流の電流(i)−電圧(v)特性の波形を電圧値(v)に対して微分演算し、予め定めた下限値V1〜上限値V2の範囲において、微分値(di/dv)が「ゼロクロス」する点の電圧値Vpk1と電流値Ipk1を、各電極ユニット761のそれぞれの電流−電圧特性に対して算出する(図20参照。)。「ゼロクロス」する点とは、電気化学的電流の微分値(di/dv)が正値から負値に変化する点、又は負値から正値に変化する点であり、電流ピークを与える電圧値Vpk1及び電流値Ipk1に対応する。図20には、電圧値の増大に伴い、微分値(di/dv)が負値から正値に変化する点の電圧値Vpk1と電流値Ipk1が示されている。「ゼロクロス点」は、ノイズによる影響やピーク検出範囲(下限値V1〜上限値V2)設定による悪影響を受けないために、負値から正値に転じて連続3点が正値を保持している場合の最初に正値となった点の電圧値を、電圧値Vpk1として採用するのが好ましい。
【0034】
ベースライン近似部312は、図18のステップS224〜図19のS228に示す手順に従い、電気化学的電流の電流−電圧特性のベースライン(バックグラウンド)を、各電極ユニット761のそれぞれの電流−電圧特性に対して近似する(図21及び図22参照。)。ピーク電流値演算部313は、図19に示すステップS229〜S230の手順に従い、電気化学的電流の電流−電圧特性曲線における真のピーク電流値Ipk2を、電圧値算出部311が算出したゼロクロス電流値Ipk1からベースライン近似部312が算出したバックグラウンド(ベースライン)の電流値Ibgを差し引いて求める。
【0035】
グループ正常判定モジュール320は、図25のフローチャートに示す一連の手順を用いて、ピーク電流検出モジュール310が算出した真の電気化学的電流(真の検出信号)のピーク値(ピーク電流値)のデータ群中から、異常データを除外し、データ群がその後の型判定を行うに相応しいデータ群であるかどうかを判定する。即ち、グループ正常判定モジュール320は、図4に示す基板714上に複数配列される電極ユニット761にそれぞれ含まれる複数のSNP1検出用電極551,SNP2検出用電極552,コントロール電極553からそれぞれ得られる一群の電流値Ipk2のデータから、所定の基準を満たさないデータを異常値として排除し(歯抜け異常判定S31)、又グループ単位で所定の基準を満たさない場合にはそのグループを異常として、その後の判定モジュールでの演算対象外とする(バラツキ異常判定S32)。
検査有効性判定モジュール325は、図32に示すフローチャートに示す一連の手順に従い、検査が有効であるかどうかを判定する。即ち、図4に示す基板714上に複数配列される電極ユニット761にそれぞれ含まれる複数のポジティブ・コントロール用電極(作用極)554から検出されるポジティブ・コントロール電流と、ネガティブ・コントロール用電極(作用極)555から検出されるネガティブ・コントロール電流とを所定の基準値と大小関係比較を行い、ネガティブ・コントロール電流が所定の基準よりも小さく、ポジティブ・コントロール電流が所定の基準よりも大きい場合に、当該検査が有効である、との判断を行い、その後の判定モジュールでの演算に移行する。
【0036】
有無判定モジュール330は、図27に示すフローチャートに示す一連の手順に従い、対象とする核酸が存在するかどうかを判定する。詳細は、後述する。
【0037】
型判定モジュール340は、SNPの型が、SNP="G",SNP="T"の2種の内のどちらであるか、又、それがG/Gホモ型であるか、G/Tヘテロ型であるか、又はT/Tホモ型であるか等を、図28に示すフローチャートに示す一連の手順で判断する。詳細は、後述する。
【0038】
図8に示したように、CPU300には、バス303を介して、波形判定用電圧範囲記憶部351,傾き許容範囲記憶部352,ピーク探索電圧値範囲記憶部353,ゼロクロス値記憶部354,変曲点記憶部355,交点電圧記憶部356,オフセット電圧記憶部357,ベースライン電流値記憶部358,平均値・標準偏差記憶部360,グループ正常判定記憶部361,SLL記憶部362,ESLL記憶部363,MIR記憶部364,規格化座標記憶部365、ベクトル長記憶部366、角度記憶部367、角度パラメータ記憶部368,及び分類結果記憶部369が接続されている。
【0039】
波形判定用電圧範囲記憶部351は、電流波形判定モジュール302が、SNP1検出用電極551,SNP2検出用電極552,コントロール電極553を介して測定される電流(検出信号)のベースラインの傾きを計算する範囲としての「範囲下限電圧VLo」と「範囲上限電圧VHi(VLo < VHi)」を記憶する。傾き許容範囲記憶部352は、電流波形判定モジュール302が、検出信号のベースラインの傾きの許容範囲を規定するパラメータとしての「傾き下限値(Coef Lo)」と「傾き上限値(Coef Hi)」を記憶する。
【0040】
ピーク探索電圧値範囲記憶部353は、予め決定されたピーク探索電圧値範囲[V1,V2]を設定パラメータとして格納しておき、電圧値算出部311がピーク探索電圧値範囲[V1,V2]を読み出し可能にしている。電気化学的電流の電流−電圧特性が示す電流ピークは、測定条件を固定すればほぼ一定の電圧範囲に現れるので、ピーク探索電圧値範囲[V1,V2]を設定パラメータとして決めておくのである。ゼロクロス値記憶部354は、電圧値算出部311が算出した「ゼロクロス値(ゼロクロス電圧値Vpk1,ゼロクロス電流値Ipk1)」を、図4に示した基板714上の、すべての電極ユニット761毎に整理して格納する。
【0041】
変曲点記憶部355は、ベースライン近似部312の演算に必要な変曲点電圧Vifpを記憶する。「変曲点電圧Vifp」は、図21に示すように、電流ピークを与えるゼロクロス電圧値Vpk1から電圧を負の方向に遡って(電圧を減少させて)微分値が最小となる電圧である。交点電圧記憶部356は、交点電圧Vcrsを記憶する。「交点電圧Vcrs」とは、図22に示すように、電気化学的電流の電流−電圧特性曲線とこの電流−電圧特性曲線の接線との交点が与える電圧である。オフセット電圧記憶部357は、図22のように、交点電圧Vcrsを起点として、設定パラメータとして規定されるオフセット値だけ電圧を負の方向に遡った(電圧を減少させた)オフセット電圧Vofsを記憶する。
【0042】
ベースライン電流値記憶部358は、ピーク電流値演算部313の演算に必要なバックグラウンド(ベースライン)の電流値Ibgを記憶する。「バックグラウンド(ベースライン)の電流値Ibg」とは、図22に示すように、ベースライン近似部312が算出したベースラインの近似一次式に、電圧値算出部311が算出したゼロクロス電圧値Vpk1を代入して得られる電流値である。
【0043】
平均値・標準偏差記憶部360は、グループ正常判定モジュール320、検査有効性判定モジュール325、有無判定モジュール330、及び型判定モジュール340が、それぞれ算出した算出したポジティブ・コントロール電極554の電流値の平均値X、ネガティブ・コントロール電極555の電流値の平均値X、有無検出用電極550の電流値の平均値X、有無検出用電極550に対応するコントロール電極(C)553の電流の平均値X、SNP1検出用電極551の電流の平均値X1,SNP1検出用電極551に対応するコントロール電極(C1)553の電流の平均値Xc1,SNP2検出用電極552の電流の平均値X2,SNP2検出用電極552に対応するコントロール(C2)電極553の電流の平均値Xc2、ポジティブ・コントロール電極554の電流値の標準偏差σ、ネガティブ・コントロール電極555の電流値の標準偏差σ、有無検出用電極550からのピーク電流値の標準偏差σ、対応するコントロール電極53からのピーク電流値の標準偏差σc、SNP1検出用電極551からのピーク電流値の標準偏差σ1、対応するコントロール電極553からのピーク電流値の標準偏差σc1、SNP2検出用電極552からのピーク電流値の標準偏差σ2、対応するコントロール電極553からのピーク電流値の標準偏差σc2、等を記憶する。これらの平均値X、X,X、X1,Xc1,X2,Xc2や標準偏差σ、σp、σn、σ1,σc1,σ2,σc2等は、グループ正常判定モジュール320、検査有効性判定モジュール325、有無判定モジュール330、及び型判定モジュール340の演算の要請に応じて随時読み出される。
【0044】
グループ正常判定記憶部361は、グループ正常判定モジュール320において演算された歯抜けデータ数Nr、CV値CV及び、グループ正常判定モジュール320の演算に必要な歯抜け電流基準値MS(Minimum Signal)、歯抜け許容割合P、基準CV値CV0(%)等を記憶する。「CV(変動係数)値」とは、対象とする一群のデータの標準偏差を対応する平均値で除した値を100倍して百分率表示したものである。分散や変動は標本の単位を持っているので、複数の標本群間のバラツキの程度を比較できないので、それぞれの平均値で割って無名数としたものである。SLL記憶部362は、検査有効性判定モジュール325の演算に必要なネガティブ・コントロール基準値NCLL(Negative Control Lower Limit),NCUL(Negative Control Upper Limit),ポジティブ・コントロール基準値PCLL(Positive Control Lower Limit)、有無判定モジュール330の演算に必要な信号増加量判定基準SL(+)、SL(−)を記憶している。信号判定基準SLL(Signal Limit Level)記憶部362に記憶されている各パラメータは、信号の絶対値若しくは、コントロール電極553に対する信号増加量の判定アルゴリズムの判定基準を与える設定パラメータである。ESLL記憶部363は、検査有効性判定モジュール325の演算に必要なポジティブ・コントロール実効変動係数PESL、有無判定モジュール330の演算に必要な実効変動係数ESLLを記憶している。ポジティブ・コントロール実効変動係数PESL(Positive Control Effective Signal Lower Limit)及び実効変動係数ESLL(Effective Signal Lower Limit)は、信号増加量が標準偏差σに対して何倍あるかの判定の下限を与える設定パラメータである。即ち、信号増加の信頼性を表す指標となる。
【0045】
MIR記憶部364は、型判定モジュール340の演算に必要な電流増加割合基準値MIRを記憶している。電流増加割合基準値MIR(Minimum Increase Ratio)は型判定で必要な電流増加量のコントロール電流値に対する比の下限を与える設定パラメータである。図31に示した判定基準となる種々の設定パラメータを示す概念イメージにおいては、原点からの距離としてMIRが記載されている。これは、SNP1検出用電極551の電流の平均値X1のSNP1検出用電極551に対応するコントロール電極(C1)553の電流の平均値Xc1からの増加量、及びSNP2検出用電極552の電流の平均値X2のSNP2検出用電極552に対応するコントロール(C2)電極553の電流の平均値Xc2からの増加量が、小さい場合(半径MIRの円よりも原点に近い位置にデータ点がある場合)には、システムとして期待される電流値が得られておらず、型判定を行うに値しない(判定不能)、ということを意味している。
【0046】
角度パラメータ記憶部368は、型判定モジュール340の演算に必要な野生型角度下限値Wmin、野生型角度上限値Wmax、ヘテロ型角度下限値Hmin、ヘテロ型角度上限値Hmax、変異型角度下限値Mmin、変異型角度上限値Mmaxを記憶している。図31に示した判定基準となる種々の設定パラメータを示す概念イメージにおいて、原点からデータ点までのベクトルが正のX軸となす角度が、野生型角度下限値Wminと野生型角度上限値Wmaxの間の角度の場合には、野生型と判定し、ヘテロ型角度下限値Hminとヘテロ型角度上限値Hmaxの間の角度の場合には、ヘテロ型と判定し、変異型角度下限値Mminと変異型角度上限値Mmaxの間の角度の場合には、変異型と判定する。図31に示している点では、ヘテロ型と判定されることになる。
【0047】
分類結果記憶部369は、有無判定モジュール330及び型判定モジュール340が分類した種々の分類結果を記憶している。
【0048】
図示を省略しているが、CPU300には、バス303を介してとインタフェースが接続され、このインタフェースを介してローカルバス(図示省略)を通じて図7に示す制御機構15との間でデータの送受信を行うことができる。
【0049】
図8において、入力装置304はキーボード、マウス、ライトペン又はフレキシブルディスク装置などで構成される。入力装置304より塩基配列判定実行者は、入出力データを指定したり、設定パラメータや許容誤差の値及び誤差の程度を設定できる。更に、入力装置304より出力データの形態等を設定することも可能で、又、演算の実行や中止等の指示の入力も可能である。又、出力装置305及び表示装置306は、例えば、それぞれプリンタ装置及びディスプレイ装置等により構成されることが可能である。表示装置306は入出力データや判定結果や判定パラメータ等を表示する。データ記憶部(図示省略)は入出力データや判定パラメータ及びその履歴や演算途中のデータ等を記憶する。
【0050】
以上説明したように、本発明の実施の形態に係る塩基配列システムによれば、核酸の存在の有無判定、SNPのホモ・ヘテロの型判定を、高い精度で、実態に即してできる。
【0051】
(塩基配列判定方法)
図14に示すフローチャートを用いて、本発明の実施の形態に係る塩基配列判定方法を説明する。なお、以下に述べる塩基配列判定方法は、一例であり、この変形例を含めて、これ以外の種々の塩基配列判定方法により、実現可能であることは勿論であるが、いずれにせよ、図3に示すようなプローブDNA571,572,573が固定化されたチップカートリッジ11にサンプルDNA581,582,583を含む薬液(試料)を注入してハイブリダイゼーション反応などを生じさせ、バッファによる洗浄、挿入剤の導入による電気化学反応を経て測定系12を用いて図13に示すような電気化学的電流の電流−電圧特性を求めることが基礎となる。電気化学的電流の電流−電圧特性から、各プローブDNA571,572,573のハイブリダイゼーション反応に定量的に対応するピーク電流値を算出する。そして算出されたピーク電流値を統計的にデータ処理し、これにより核酸の存在の有無を判定し、或いは核酸の一塩基多型の型判定をする。
【0052】
図14に示すフローチャートの説明の前に、先ず、図9〜図11を用いて、プローブDNAとサンプルDNAとのハイブリダイゼーションについて説明する。ハイブリダイゼーション工程は、図4及び図5に示したチップカートリッジ11を用いれば良い。
【0053】
図9では、サンプルDNA581が塩基配列CTGAG……を有する標的塩基配列(SNP="G"のサンプルDNA)である場合を示す。図9(a)に示すように、塩基配列GACTC……を有するプローブDNA571が固定化された作用極(SNP1検出用電極)551は、塩基配列CTGAG……を有する標的塩基配列(SNP="G"のサンプルDNA)581と配列が完全にマッチングしているので2本鎖を形成する。しかし、1塩基でも配列が異なると2本鎖を形成しないので、図9(b)に示すように、塩基配列GAATC……を有するプローブDNA572が固定化された作用極(SNP2検出用電極)552は、標的塩基配列(SNP="G"のサンプルDNA)581とは、2本鎖を形成できない。又、完全に配列が異なれば、当然に2本鎖を形成できないので、図9(c)に示すように、塩基配列CAGTG……を有するプローブDNA(コントロールDNA)573が固定化された作用極(コントロール電極)553は、標的塩基配列(SNP="G"のサンプルDNA)581とは、2本鎖を形成できない。
【0054】
図10では、サンプルDNA582が塩基配列CTTAG……を有する標的塩基配列(SNP="T"のサンプルDNA)である場合を示す。同様に、チップカートリッジ11によるハイブリダイゼーション工程により、図10(b)に示すように、塩基配列GAATC……を有するプローブDNA572が固定化された作用極(SNP2検出用電極)552は、塩基配列CTTAG……を有する標的塩基配列(SNP="T"のサンプルDNA)582と配列が完全にマッチングしているので2本鎖を形成する。しかし、1塩基でも配列が異なると2本鎖を形成しないので、図10(a)に示すように、塩基配列GACTC……を有するプローブDNA571が固定化された作用極(SNP1検出用電極)551は、標的塩基配列(SNP="T"のサンプルDNA)582とは、2本鎖を形成できない。又、完全に配列が異なれば、当然に2本鎖を形成できないので、図10(c)に示すように、塩基配列CAGTG……を有するプローブDNA(コントロールDNA)573が固定化された作用極(コントロール電極)553は、標的塩基配列(SNP="T"のサンプルDNA)582とは、2本鎖を形成できない。
【0055】
図11では、サンプルDNA583が塩基配列CTGAG……を有する標的塩基配列(SNP="G"のサンプルDNA)と塩基配列CTTAG……を有する標的塩基配列(SNP="T"のサンプルDNA)とのヘテロタイプである場合を示す。同様に、チップカートリッジ11によるハイブリダイゼーション工程により、図11(a)に示すように、塩基配列GACTC……を有するプローブDNA571が固定化された作用極(SNP1検出用電極)551は、塩基配列CTGAG……と塩基配列CTTAG……との両方を有するG/Tヘテロの標的塩基配列(SNP="G/T"ヘテロのサンプルDNA)583と配列が完全にマッチングしているので2本鎖を形成する。又、図11(b)に示すように、塩基配列GAATC……を有するプローブDNA572が固定化された作用極(SNP2検出用電極)552は、塩基配列CTTAG……を有するG/Tヘテロの標的塩基配列(SNP="G/T"ヘテロのサンプルDNA)583と配列が完全にマッチングしているので2本鎖を形成する。しかし、完全に配列が異なれば、当然に2本鎖を形成できないので、図11(c)に示すように、塩基配列CAGTG……を有するプローブDNA(コントロールDNA)573が固定化された作用極(コントロール電極)553は、G/Tヘテロの標的塩基配列(SNP="G/T"ヘテロのサンプルDNA)583とは、2本鎖を形成できない。
【0056】
図4及び図5に示した構成とは異なる場合、即ち基板(チップ)上に平面状のパッキンを搭載し、カセット(チップカートリッジ)内に流路を形成する構成をとる場合には、カセット(検出チップ)21内の流路が長くなり、不必要な試薬量が多くなる。又、カセット(検出チップ)21に薬液(試料)を注入する場合には、流路が基板上のみならずカセット(検出チップ)21内にも長く存在するため、基板以外の不要な部分に薬液(試料)が流れ込み、無駄となってしまう。又、パッキンに対するカセット(検出チップ)21の密着性が難しく、パッキンとカセット(検出チップ)21との間でリークが生じてしまい、送液の不具合が発生することが多い。本実施の形態のような構成により、不必要な試薬量が少なくなり、又パッキン、基板及びカセット(検出チップ)21の密着性が高くなり、送液の安定性が増す。
【0057】
図12は、前述の図9と同様に、塩基配列CTGAG……を有する標的塩基配列(サンプルDNA)581を用いて、ハイブリダイゼーションした各作用極551,552,553に、挿入剤591を導入した状態を示す。図12(a)に示すように、塩基配列GACTC……を有するプローブDNA571が固定化された作用極(SNP1検出用電極)551は、塩基配列CTGAG……を有する標的塩基配列(サンプルDNA)581と配列が完全にマッチングし、2本鎖を形成しているので、2本鎖DNAに挿入剤591が結合する。しかし、図12(b)に示すように、塩基配列GAATC……を有するプローブDNA572が固定化された作用極(SNP2検出用電極)552は、標的塩基配列(サンプルDNA)581とは、2本鎖を形成できないので、挿入剤591が結合できない。又、図12(c)に示すように、塩基配列CAGTG……を有するプローブDNA(コントロールDNA)573が固定化された作用極(コントロール電極)553には、標的塩基配列(サンプルDNA)581とは、2本鎖を形成できないので、挿入剤591が結合できない。
【0058】
図13は、各作用極551,552,553上に固定化されたプローブDNA571,572,573にハイブリダイゼーションした2本鎖DNAに結合した挿入剤591からの電気化学的電流、又は2本鎖DNAに挿入剤591が結合できなかった場合の電流と電圧との関係を示す図である。図13の曲線(a)は、図12(a)に対応する。即ち、曲線(a)は、プローブDNA571と標的塩基配列(サンプルDNA)581と配列が完全にマッチングし、2本鎖を形成し、この2本鎖DNAに挿入剤591が結合した場合の電気化学的電流の電流−電圧特性であり、大きな電流値のピークを示す。しかし、図13の曲線(b)は、図12(b)に対応し、プローブDNA572と標的塩基配列(サンプルDNA)581とが2本鎖を形成できず、挿入剤591が結合できない場合の電流−電圧特性であり、曲線(a)に比して小さな電流値のピークを示す。又、図13の曲線(c)は、図12(c)に対応し、コントロールDNA573と標的塩基配列(サンプルDNA)581とが2本鎖を形成できず、挿入剤591が結合できない場合の電流−電圧特性であり、曲線(b)に比して、更に小さな電流値のピークを示す。図13(b)、(c)で見られる小さな電流値ピークは、図12(b)、(c)で示したように、電極552、553表面にわずかに吸着した挿入剤591起因の電流である。
【0059】
前置きが長くなったが、図14に示すフローチャートを用いて、本発明の実施の形態に係る塩基配列判定方法を説明する:
(イ)先ず、チップカートリッジ11にサンプルDNAを含む薬液(試料)を注入してハイブリダイゼーション反応などを生じさせ、挿入剤の導入による電気化学反応による電流−電圧特性を測定系12を用いて、各電極毎に測定する。図4に基板714上に多数の電極ユニット761を模式的に示したが、この多数の電極ユニット761に対応して、プローブDNA571,572,573がそれぞれ固定化されるSNP1検出用電極551,SNP2検出用電極552,コントロール電極553は、等価な電極として多数存在し、それに対応して多数のデータが取得される。ステップS101において、ノイズ除去モジュール301が、プローブDNA571,572,573が固定化された各電極毎に、測定された一群のデータのそれぞれを、平準化(スムージング)により、ノイズを除去する。スムージングは、既に説明したように、図15に示すような、単純移動平均法を用いれば良い。 以降の処理は、すべて平準化処理後のデータに対して演算を行う。
【0060】
(ロ)次に、ステップS102において、電流波形判定モジュール302が、各電極毎に測定された電流波形(電流−電圧特性)のベースラインの傾きをそれぞれ求め、それぞれのベースラインの傾きにより、それぞれの検出信号(電流波形)の正常・異常を判定し、異常な検出信号を演算対象外とする。ステップS102における電流波形判定モジュール302の処理の詳細は、図16のフローチャートを用いて後述する。
【0061】
(ハ)次に、ステップS103において、ピーク電流検出モジュール310が、各電極毎に測定された検出信号のピーク値(ピーク電流値)をそれぞれ検出する。ステップS103におけるピーク電流検出モジュール310の処理の詳細は、図18及び図19のフローチャートを用いて後述するが、ステップS103における処理で、図12に示したような挿入剤591からの電気化学的電流から、これ以外のバックグラウンド電流を差し引いて、真の検出信号のピーク値(ピーク電流値)が、等価な各電極毎に、それぞれ一群のデータとして求められる。
【0062】
(ニ)ステップS103でのバックグラウンド電流成分の除去後、一群のデータのそれぞれに、ステップS104における信号処理がなされる。即ち、ステップS104において、グループ正常判定モジュール320が、一群のデータが正常であるか否かを判定する。ステップS104におけるグループ正常判定モジュール320の処理の詳細は、図25のフローチャートを用いて後述する。
【0063】
(ホ)その後、ステップS105において、検査有効性判定モジュール325が、検査が有効であるか無効であるかを判定する。検査有効性判定モジュール325の処理の詳細は、図32のフローチャートを用いて後述する。
【0064】
(ヘ)次に、ステップS106において、有無判定モジュール330が、核酸の存在の有無を判定、又は型判定モジュール340が、核酸の一塩基多型(SNP)の型判定を行う。ステップS106における有無判定モジュール330及び型判定モジュール340のそれぞれ処理の詳細は、図26〜図28のフローチャートを用いて後述する。
【0065】
図14に示すフローチャートに示す本発明の実施の形態に係る塩基配列判定方法によれば、チップカートリッジ11や測定系12に異常があった場合や、データにバラつきがあった場合等でも、それらの異常なデータを除外して、正確にある核酸が存在するかどうか、SNPの型が、どの型であるか、又、それがホモ型であるか、ヘテロ型であるか、を判断できる。以下具体的に、図14に示すフローチャートの各ステップを説明する。
【0066】
[ステップS102:電流波形の正常・異常の判定]
図1の測定系12が測定したチップカートリッジ11における電気化学的電流の信号(電流波形)は、特徴的に、図17に示したような電流−電圧特性の波形をとる。電気信
を発生する物質(挿入剤)に特有の電圧において、図17に示したようなピーク形状を有する電流−電圧特性の波形を持っている。図17では、データ1とデータ2の2種類の電流−電圧特性を示しているが、データ1の電流−電圧特性が示すベースラインの傾きに比して、データ2の電流−電圧特性が示すベースラインの傾きの方が大きく、データ2の電流−電圧特性では、ピーク形状が不明確で肩(ショルダー)的な変化を示している。
【0067】
理想的には、電気化学的電流の電流−電圧特性は、ピーク電流を発生する電圧以下ではほぼ「零」の電流値となる。ところが、基板714に何らかの不具合が生じた場合には、しばしば、データ2のように、ベースラインが傾いた電流−電圧特性となる。データ2のような場合には、ピーク電流値を正しく検出することができないため、図14のステップS102において、「異常」として排除する必要がある。
【0068】
ステップS102における電流波形判定モジュール302の処理の詳細は、図16のフローチャートに示す通りであり:
(a)先ず、ステップS201において、各電極毎に測定された電流波形(電流−電圧特性)のベースラインの傾きを計算する範囲を抽出し、設定する。この範囲の抽出では、設定パラメータとして、範囲下限電圧VLoと範囲上限電圧VHi(VLo < VHi)を入力装置304を用いて設定し、波形判定用電圧範囲記憶部351に格納する。更に、ベースラインの傾きの許容範囲を規定するパラメータとして「傾き下限値(Coef Lo)」と「傾き上限値(Coef Hi)」を設定し、傾き許容範囲記憶部352に格納する。
【0069】
(b)次に、ステップS202において、電流波形判定モジュール302は、波形判定用電圧範囲記憶部351に格納された範囲下限電圧VLo及び範囲上限電圧VHiとを読み出し、この設定範囲で、ベースラインの傾きの近似式を導出する。範囲下限電圧VLo及び範囲上限電圧VHiは、ベースラインの傾きを算出するための領域を規定するパラメータであり、範囲下限電圧VLoと範囲上限電圧VHiとの間の電圧値で、各電極毎に測定された電流−電圧特性の波形に対して最小二乗近似によって、直線(ベースライン)を近似する。
【0070】
(c)次に、ステップS203において、電流波形判定モジュール302は、各電極毎に測定された電流波形(電流−電圧特性)のそれぞれについて、範囲下限電圧VLo及び範囲上限電圧VHiを読み出し、それぞれ開始点及び終了点とし、これらの開始点及び終了点を、ステップS202において導出した近似式に代入し、各電極毎に測定された電流波形(電流−電圧特性)のベースラインの傾き(b)を算出する。
【0071】
(d)次に、ステップS204において、電流波形判定モジュール302は、傾き許容範囲記憶部352から、「傾き下限値(Coef Lo)」と「傾き上限値(Coef Hi)」を読み出し、ステップS203において、各電極毎に算出したベースラインの傾き(b)が、「傾き下限値(Coef Lo)」と「傾き上限値(Coef Hi)」の間に存在するか、それぞれ判定する。
【0072】
(e)ステップS204において、特定の電極で測定された電流波形(電流−電圧特性)のベースラインの傾き(b)が、「傾き下限値(Coef Lo)」と「傾き上限値(Coef Hi)」の間に存在すれば、「正常波形」と判断して、図14のステップS103に進む。ステップS204において、特定の電極で測定された電流波形(電流−電圧特性)が、「傾き下限値(Coef Lo)」と「傾き上限値(Coef Hi)」の範囲外であると判断された場合には、その電極で測定された電流波形(電流−電圧特性)は、「異常波形」であると判断して、ステップS205において、電流波形判定モジュール302は、その電極で測定された電流波形(電流−電圧特性)に対し「エラー判定」を下し、適宜、表示装置306及び出力装置305に出力させる。
【0073】
図14に示したステップS103の手順は、図4に示した基板714上の、すべての電極ユニット761について実施される。
【0074】
[ステップS103:ピーク電流値の検出]
ステップS103におけるピーク電流検出モジュール310の処理の詳細を、図18及び図19のフローチャートを用いて説明する。ステップS103において、測定系12が測定したそれぞれの電極ユニット761による電流−電圧特性の波形から、それぞれの正味のピーク電流値を検出する段階(ステップ)は、ステップS221〜S223における電流ピークを与える電圧値を算出する段階(図20参照。);ステップS224〜S228におけるベースラインを近似する段階(図21及び図22参照。);及びステップS229〜S230におけるピーク電流値を算出する段階(図23参照。)を、各電極ユニット761のそれぞれの電流−電圧特性に対して実行する手順からなる:
(a)チップカートリッジ11で測定された電気化学的電流の電流−電圧特性が示す電流ピークは、ほぼ一定の電圧範囲に現れる。このため、ステップS221においては、図8に示した入力装置304を用い、予め、ピーク探索電圧値範囲[V1,V2]を設定パラメータとしてピーク探索電圧値範囲記憶部353に格納しておく。即ち、図20に示すように、電気化学的電流のピーク電流探索は、下限値V1〜上限値V2の範囲で行う。先ず、ステップS222において、ピーク電流検出モジュール310の電圧値算出部311は、電気化学的電流の電流(i)−電圧(v)特性の波形を電圧値(v)に対して微分演算する。そして、ステップS223において、電圧値算出部311は、下限値V1〜上限値V2の範囲において、電気化学的電流の微分値(di/dv)が「ゼロクロス」する点の電圧値Vpk1と電流値Ipk1を算出する(図20参照。)。「ゼロクロス」する点とは、電気化学的電流の微分値(di/dv)が正値から負値に変化する点、又は負値から正値に変化する点であり、電流ピークを与える電圧値Vpk1及び電流値Ipk1に対応する。図20には、電圧値の増大に伴い、微分値(di/dv)が負値から正値に変化する点の電圧値Vpk1と電流値Ipk1が示されている。「ゼロクロス値」は、ノイズによる影響やピーク検出範囲(下限値V1〜上限値V2)設定による悪影響を受けないために、負値から正値に転じて連続3点が正値を保持している場合の最初に正値となった点の電圧値を、電圧値Vpk1として採用する。ステップS223において算出された「ゼロクロス値(ゼロクロス電圧値Vpk1,ゼロクロス電流値Ipk1)」は、ゼロクロス値記憶部354に、図4に示した基板714上の、すべての電極ユニット761毎に整理して、格納する。
【0075】
(b)ステップS224において、ピーク電流検出モジュール310のベースライン近似部312は、電流ピークを与えるゼロクロス電圧値Vpk1から電圧を負の方向に遡って(電圧を減少させて)、図21に示すように、微分値が最小となる電圧を変曲点電圧Vifpとする。変曲点電圧Vifpは、変曲点記憶部355に格納する。そして、ベースライン近似部312は、ゼロクロス電圧値Vpk1と変曲点電圧Vifpをゼロクロス値記憶部354及び変曲点記憶部355からそれぞれ読み出し、ステップS225において、電流−電圧特性曲線への接線の起点を求めるための一次式:
y = ax+b ……(3)
を、ゼロクロス電圧値Vpk1と変曲点電圧Vifpの間で求める。例えば、ゼロクロス電圧値Vpk1と変曲点電圧Vifpの間における電流−電圧特性の波形データを、最小二乗法により式(3)のように直線近似する。図21に示す例では、一次式の係数a=-1.397×10-10、定数b=3.396×10-8と求められ、図22に示す例では、一次式の係数a=-2.899×10-11、定数b=4.504×10-9と求められる。更に、ベースライン近似部312はステップS226において、電流−電圧特性の波形と式(3)の近似直線との交点を与える交点電圧Vcrsを図22に示すように求め、この交点電圧Vcrsを交点電圧記憶部356に格納する。そして、ベースライン近似部312はステップS227において、交点電圧Vcrsを起点として、設定パラメータとして規定しているオフセット値だけ電圧を負の方向に遡った(電圧を減少させた)オフセット電圧Vofsを、図22のように求める。求められたオフセット電圧Vofsは、オフセット電圧記憶部357に格納する。更に、ベースライン近似部312は、オフセット電圧Vofsをオフセット電圧記憶部357から、交点電圧Vcrsを交点電圧記憶部356からそれぞれ読み出し、ステップS228において、オフセット電圧Vofsと交点電圧Vcrsの間の電流−電圧特性の波形データのバックグラウンドの接線(ベースライン)となる近似一次式を最小二乗法により、図23に示すように求める。この近似一次式は、式(3)と同様な形式で表現可能で、図23に示す直線近似では、一次式の係数a=-6.072×10-12、定数b=-5.902×10-9と求められる。
【0076】
(c)この後、ピーク電流検出モジュール310のピーク電流値演算部313は、ゼロクロス電圧値Vpk1をゼロクロス値記憶部354から読み出す。そして、ステップS229において、ピーク電流値演算部313は、ステップS228で求めたベースラインの近似一次式に、ゼロクロス電圧値Vpk1を代入し、バックグラウンド(ベースライン)の電流値Ibgを求める。バックグラウンド(ベースライン)の電流値Ibgは、ベースライン電流値記憶部358に格納する。更に、ピーク電流値演算部313は、電流−電圧特性の波形のピークを示すゼロクロス電流値Ipk1を、ゼロクロス値記憶部354から読み出し、ステップS230において、ゼロクロス電流値Ipk1から、バックグラウンド(ベースライン)の電流値Ibgを、式(4)に示すように差し引く:
pk2 =abs(Ipk1−Ibg) ……(4)
式(4)を各電極ユニット761のそれぞれのSNP1検出用電極(SNP="G"検出用電極)551,SNP2検出用電極(SNP="T"検出用電極)552,コントロール電極553のそれぞれの電流−電圧特性に対して適用することにより、各電極ユニット761のそれぞれの電極551,552,553から、それぞれ真の電流値Ipk2が算出される。
【0077】
[ステップS104:グループ正常判定]
以上説明したように、図14のステップS103では、ピーク電流検出モジュール310が、測定系12が測定した各電極ユニット761によるそれぞれの電流−電圧特性のピークを示すゼロクロス電流値Ipk1から、バックグラウンド(ベースライン)の電流値Ibgを差し引いて、各電極ユニット761のそれぞれのSNP1検出用電極(SNP="G"検出用電極)551,SNP2検出用電極(SNP="T"検出用電極)552,コントロール電極553毎に、それぞれの正味の電流値Ipk2が、種々の判定モードで算出され、分類される。
【0078】
しかし、例えば、複数の同等のSNP1検出用電極(SNP="G"検出用電極)551の内1電極だけ、異常に大きい値、若しくは小さい値の場合、複数の同等のSNP2検出用電極(SNP="T"検出用電極)552の内1電極だけ、異常に大きい値、若しくは小さい値を示す場合は、それら異常値を除外しないと、判定アルゴリズムに混乱を生じさせる。即ち、図4に示すように、基板714上に複数配列した同等の電極ユニット761から得られるピーク電流値Ipk2をそのまま、すべて用いて、図14のステップS106の処理、つまり、ある核酸が存在するかどうかの判定や、SNPの型が、2種の内のどちらであるか、又、それがホモ型であるか、ヘテロ型であるかを判定する場合に、問題が生じてしまう場合がある。
【0079】
そこで、図14のステップS104では、ステップS105に進む前に、グループ正常判定モジュール320が、図4に示す基板714上に複数配列される電極ユニット761によって定義されるすべての電極から得られる電流値Ipk2に対し、図25に示すフローチャートの手順をグループ単位で実施する。図25に示すフローチャートでは、ステップS31において、図24(a)に示すような歯抜け異常を判定し、その後ステップS32において、図24(b)に示すようなバラツキ異常を判定する。ステップS31における歯抜け異常の判定、及びステップS32におけるバラツキ異常の判定は、以下に示すような所定の基準で、グループ単位でその後の判定処理に値するかどうかの有効性を判断する:
(a)ステップS301では、グループ正常判定モジュール320において、ステップS31における歯抜け判定を行うかどうかの選択を行う。歯抜け判定を行う場合には、ステップS31のステップS302に進む。歯抜け判定を行わない場合には、バラツキ異常を判定するステップS32のステップS305に進む。
【0080】
(b)ステップS302において、歯抜けデータ数をカウントするため、歯抜けデータ数Nrの初期値として「0」をセットして、次のステップS303に進む。
【0081】
(c)ステップS303において、グループ正常判定記憶部361に格納されている歯抜け電流基準値MSを読み出し、各電流値Ipk2に対して、歯抜け電流基準値MSとの大小関係を比較する。歯抜け電流基準値MS以上であれば、歯抜けではないと判断して、次のステップS304に進む。歯抜け電流基準値MSよりも小さい場合には、当該データは歯抜けデータと判断し、ステップS311において、「判定対象外」とするエラーセットとする。更に、ステップS312において、歯抜けデータ数Nrをカウントアップし、ステップS303に戻る。グループ内のすべてのデータに対して、このルーチンを繰り返し、最終的なグループ内の歯抜けデータ数をカウントする。ステップS303において、歯抜けでないと判断されたデータは判定対象データとする。
【0082】
(d)ステップS304では、グループ正常判定記憶部361に格納されている歯抜け許容割合Pを読み出し、積算された歯抜けデータ数Nrと、グループに割り当てられている電極数N0をPで除算して値との大小関係を比較する。歯抜けデータ数Nrが、N0/Pよりも小さければ、バラツキ異常を判定するステップS32のステップS305に進む。歯抜けデータ数Nrが、N0/P以上である場合には、当該グループは、「歯抜け多数」につき、型判定には値しないと判断し、「グループ異常」と判定する。
【0083】
(e)ステップS305において、判定対象データグループの電流値Ipk2の平均値が「零」でないか否かの判定を行い、「零」の場合には、ステップS314において、「グループ異常」と判定する。ここで、歯抜け判定を行っている場合には、平均値が「零」になることはないはずであるが、ステップS301において、歯抜け判定を行わない場合には、ステップS305において、「グループ異常」と判定される場合が生じる。「零」でない場合には、次のステップS306に進む。
【0084】
(f)ステップS306では、判定対象データの電流値Ipk2の平均値及び標準偏差を算出し、標準偏差を平均値で除算して、CV値を求め、次のステップS307に進む。
【0085】
(g)ステップS307で、グループ正常判定記憶部361から設定パラメータとして規定した基準CV値CV0と比較し、基準CV値CV0よりも小さい場合には、次のステップS308で当該グループは正常であると判定し、更に、図14のフローチャートのステップS105の検査有効性判定へと進む。基準CV値CV0以上の場合には、図25のフローチャートのステップS315において、当該グループは、「グループ異常」と判定する。
【0086】
以上のルーチンをすべてのグループについて繰り返すことにより、グループが正常かどうかの判定を行う。
【0087】
グループ正常判定モジュール320は、電極から得られる電流値Ipk2に対し演算対象外の際に「グループ異常判定」を下し、適宜、表示装置306及び出力装置305に出力させる。これ以降のステップにおいては、すべて、真の電流値Ipk2に対する演算であるため、特に断らない限り、「真の電流値Ipk2」を単に「電流値」として記述することにする。又、基本的には、判定対象データのみを、今後の演算対象とする(歯抜けデータは除いて処理する)。
【0088】
[ステップS105:検査有効性判定]
ステップS105における検査有効性判定モジュール325の処理の詳細を、図32のフローチャートを用いて説明する。ここでは、図4に示す基板714上に複数配列される電極ユニット761にそれぞれ含まれる複数のポジティブ・コントロール用電極554,ネガティブ・コントロール用電極555からそれぞれ得られる一群の電流値Ipk2のデータについて判定を行う。
【0089】
(a)ステップS401では、先ず、検査有効性判定において、ネガティブ・コントロールでの判定を行うかどうかの選択を行う。ネガティブ・コントロールでの判定を行わない場合には、ステップS435において、検査有効性判定は「有効」として、本処理を終了する。ネガティブ・コントロールでの判定を行う場合には、次のステップS402に進む。
【0090】
(b)ステップS402では、ネガティブ・コントロール用電極555からそれぞれ得られる一群の電流値Ipk2のデータが、図14のフローチャートのステップS104におけるグループ正常判定において、「正常」と判定されていた場合には、次のステップS443に進む。図14のフローチャートのステップS104におけるグループ正常判定において、「グループ異常」と判定された場合には、ステップS441において、「判定不能(N.D.52)」を、表示装置306及び出力装置305に出力させる。
【0091】
(c)ステップS443では、ネガティブ・コントロール用電極555からそれぞれ得られる一群の電流値Ipk2の平均値X、及び標準偏差σをそれぞれ計算し、平均値・標準偏差記憶部360に格納する。
【0092】
(d)ステップS403では、SLL記憶部362から格納されている設定パラメータネガティブ・コントロール判定上限値NCUL、ネガティブ・コントロール判定下限値NCLLを、平均値・標準偏差記憶部360からは、ステップS443で算出したネガティブ・コントロール電流の平均値Xを読み出し、設定パラメータネガティブ・コントロール判定上限値NCUL、ネガティブ・コントロール判定下限値NCLLとネガティブ・コントロール電流値の平均値Xとの大小関係を比較する。ネガティブ・コントロール電流値の平均値Xが、判定下限値NCLL以上、判定上限値NCUL以下である場合には、ステップS404において、ネガティブ・コントロールは妥当とみなして、次のステップS405に進む。ネガティブ・コントロール電流値の平均値Xが、判定下限値NCLLと判定上限値NCULの範囲から逸脱している場合には、ステップS421において、ネガティブ・コントロール異常として、ステップS422に進む(図24(c)は、ネガティブ・コントロール電流値の平均値が、上限値NCULを超えてしまった場合を概念的に示している。)。
【0093】
(e)ステップS405では、ポジティブ・コントロール判定を行うかどうか選択を行う。ポジティブ・コントロールでの判定を行わない場合には、ステップS412において、検査有効性判定は「検査有効」として、本処理を終了する。ポジティブ・コントロールでの判定を行う場合には、次のステップS406に進む。
【0094】
(f)ステップS406では、ネガティブ・コントロールに対して、ステップS402において行ったと同様に、ポジティブ・コントロールに対して行う。図14のフローチャートのステップS104におけるグループ正常判定において、「正常」と判定されていた場合には、次のステップS444に進む。図14のフローチャートのステップS104におけるグループ正常判定において、「グループ異常」と判定された場合には、ステップS442において「判定不能(N.D.51)」とし、分類結果記憶部369に判定結果を分類して格納し、表示装置306及び出力装置305に出力させる。
【0095】
(g)ステップS444では、ポジティブ・コントロール用電極554からそれぞれ得られる一群の電流値Ipk2の平均値X、及び標準偏差σをそれぞれ計算し、平均値・標準偏差記憶部360に格納する。
【0096】
(h)ステップS407では、SLL記憶部362から格納されている設定パラメータポジティブ・コントロール判定下限値PCLL及びポジティブ・コントロール実効変動係数PESLを、平均値・標準偏差記憶部360からは、ステップS443で算出したネガティブ・コントロール電流の平均値X及び標準偏差σ、ステップS444で算出したポジティブ・コントロール電流の平均値X及び標準偏差σを読み出す。ポジティブ・コントロール電流値の平均値Xからネガティブ・コントロール電流値の平均値Xを減算した値とポジティブ・コントロール判定下限値PCLLと大小関係を比較する。更に、ポジティブ・コントロール電流値平均値Xから、ポジティブ・コントロール電流の標準偏差σにポジティブ・コントロール実効変動係数PESLを乗算した値を減算した値と、ネガティブ・コントロール電流値平均値Xに、ネガティブ・コントロール電流の標準偏差σにポジティブ・コントロール実効変動係数PESLを乗算した値を加算した値を比較する。ポジティブ・コントロール電流値の平均値Xからネガティブ・コントロール電流値の平均値Xを減算した値が、ポジティブ・コントロール判定下限値PCLL以上でかつ、ポジティブ・コントロール電流値平均値Xから、ポジティブ・コントロール電流の標準偏差σにポジティブ・コントロール実効変動係数PESLを乗算した値を減算した値がネガティブ・コントロール電流値平均値Xに、ネガティブ・コントロール電流の標準偏差σにポジティブ・コントロール実効変動係数PESLを乗算した値を加算した値以上である場合(ステップS408)には、ステップS409において、「検査有効」と判定し、処理を終了する。いずれか一方でも満たさない場合(ステップS410)には、ステップS411において、「検査無効」と判断する。
【0097】
(i)ステップS403で、ネガティブ・コントロール電流値の平均値が、図24(c)に示すように、判定下限値NCLLと判定上限値NCULの範囲から逸脱している場合には、ステップS421でネガティブ・コントロール異常として、ステップS422に進む。ステップS422では、ステップS406と同様にポジティブ・コントロール判定を行うかどうか選択を行う。ポジティブ・コントロールでの判定を行わない場合には、ステップS434において検査有効性判定は「検査無効」として、本処理を終了する。ポジティブ・コントロールでの判定を行う場合には、次のステップS423に進む。
【0098】
(j)ステップS423では、ステップS406において行ったと同様に、ポジティブ・コントロールに対して、図14のフローチャートのステップS104におけるグループ正常判定において、「正常」と判定されていた場合には、次のステップS445に進む。図14のフローチャートのステップS104におけるグループ正常判定において、「グループ異常」と判定された場合には、ステップS431において「検査有効性判定」は「検査無効」として、本処理を終了する。
【0099】
(k)ステップS445では、ステップS444と同様の処理を行う。即ち、ポジティブ・コントロール用電極554からそれぞれ得られる一群の電流値Ipk2の平均値X、及び標準偏差σをそれぞれ計算し、平均値・標準偏差記憶部360に格納する。
【0100】
(l)ステップS424では、SLL記憶部362から格納されている設定パラメータポジティブ・コントロール判定下限値PCLL及びポジティブ・コントロール実効変動係数PESLを、平均値・標準偏差記憶部360からは、ステップS443で算出したネガティブ・コントロール電流の平均値X及び標準偏差σ、ステップS445で算出したポジティブ・コントロール電流の平均値X及び標準偏差σを読み出す。ポジティブ・コントロール電流値の平均値Xからネガティブ・コントロール電流値の平均値Xを減算した値とポジティブ・コントロール判定下限値PCLLと大小関係を比較する。更に、ポジティブ・コントロール電流値平均値Xから、ポジティブ・コントロール電流の標準偏差σにポジティブ・コントロール実効変動係数PESLを乗算した値を減算した値と、ネガティブ・コントロール電流値平均値Xに、ネガティブ・コントロール電流の標準偏差σnにポジティブ・コントロール実効変動係数PESLを乗算した値を加算した値を比較する。ポジティブ・コントロール電流値の平均値Xからネガティブ・コントロール電流値の平均値Xを減算した値が、ポジティブ・コントロール判定下限値PCLL以上でかつ、ポジティブ・コントロール電流値平均値Xから、ポジティブ・コントロール電流の標準偏差σにポジティブ・コントロール実効変動係数PESLを乗算した値を減算した値がネガティブ・コントロール電流値平均値Xに、ネガティブ・コントロール電流の標準偏差σにポジティブ・コントロール実効変動係数PESLを乗算した値を加算した値以上である場合(ステップS425)は、ステップS426において「検査無効」と判断し、本処理を終了する。ステップS424で、上記条件のいずれか一方でも満たさない場合(ステップS432)も、ステップS433において「検査無効」と判断し、本処理を終了する。
【0101】
[2種類のアルゴリズムから選択]
図14のステップS106は、図26に示すようにステップS332における検査有効性の判定において検査が有効であると判定された場合は、ステップS333に進み、2種類のアルゴリズムの選択を行う。ステップS332における検査有効性の判定において検査が有効でないと判定された場合は、ステップS334において「検査無効」として、分類結果記憶部369に判定結果を分類して格納し、表示装置306及び出力装置305に出力させ、信号処理を終了する。ステップS333においては、ある核酸が存在するかどうかを判定する有無判定アルゴリズムのフローに進むのか、SNPの型が、野生型ホモ型か、ヘテロ型か、変異型ホモ型か、即ち、例えば、G/Gホモ型であるか、G/Tヘテロ型であるか、又はT/Tホモ型であるか等を判定するSNP型判定アルゴリズムのフローに進むのかを決定する。
【0102】
[ステップS106その1:核酸の存在の有無判定]
図26のステップS333において、ある核酸が存在するかどうかを判定するアルゴリズムのフローに進むと決定された場合は、図27に示すフローチャートの手順で有無判定モジュール330が処理を行う。
【0103】
図2において、SNP1検出用電極551,SNP2検出用電極552,コントロール電極553、参照極561,562及び対極502が検出チップ上に配列された電極ユニットを示したが、ある核酸が存在するかどうかを判定するアルゴリズムの場合は、SNP1検出用電極551,SNP2検出用電極552のいずれか一方で良い。即ち、図4に示す基板714上には、SNP1検出用電極551,SNP2検出用電極552のいずれか一方を検出用電極(作用極)として備える電極ユニット761が複数配列されている。ここでは、図2に示したSNP1検出用電極551を、対象核酸を検出する有無検出用電極(作用極)550であると仮定して説明する。
【0104】
(a)図27に示すフローチャートのステップS341では、有無判定モジュール330は、コントロール電極553からの電流値がグループ正常判定モジュール320での処理の結果が「グループ正常」であるか否かを判定する。ステップS341で、「グループ正常」であると判定された場合には、次のステップS342に進む。グループ正常判定モジュール320での処理の結果、「グループ異常」であると判定された場合には、「判定不能(N.D.21)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS348において「判定不能(N.D.21)」を表示装置306及び出力装置305に出力させる。
【0105】
(b)更に、ステップS342において、有無判定モジュール330は、判定対象となる複数の有無検出用電極(作用極)550からの電流値がグループ正常判定モジュール320での処理の結果、「グループ正常」であるか否かを判定する。ステップS342において、「グループ正常」であると判定された場合には、次のステップS343に進む。グループ正常判定モジュール320での処理の結果、「グループ異常」であると判定された場合には、「判定不能(N.D.22)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS349において「判定不能(N.D.22)を表示装置306及び出力装置305に出力させる。
【0106】
(c)次に、ステップS343において判定対象となる複数の有無検出用電極(作用極)550からの電流値の平均値(X)、標準偏差(σ)、対応するコントロール電極553からの電流値の平均値(Xc)及び標準偏差(σc)を算出し、平均値・標準偏差記憶部360に格納する。
【0107】
(d)ステップS344において、有無判定モジュール330は、有無検出用電極(作用極)550からの電流値の平均値(X)と対応するコントロール電極553からの電流値の平均値(Xc)を平均値・標準偏差記憶部360から読み出し、平均値の差(X−X)を演算する。更に、「−」判定信号増加量上限基準SL(−)をSLL記憶部362から読み出す。そして、ステップS344において、平均値の差(X−Xc)と、「−」判定信号増加量上限基準SL(−)との大小関係を比較する。ステップS344で「−」判定信号増加量上限基準SL(−)以下と判断された場合には、有無検出用電極(作用極)550からは対象とする核酸の電気化学的電流の電流値が観測されなかったと判断し、分類結果記憶部369に判定結果を分類して格納し、ステップS350において、表示装置306に「−(検出感度以下)」の判定を表示させる。一方、ステップS344において、「−」判定信号増加量上限基準SL(−)よりも大きい場合には、ステップS345に進む。
【0108】
(e)ステップS345では、有無判定モジュール330は、有無検出用電極(作用極)550からの電流値の平均値(X)と対応するコントロール電極553からの電流値の平均値(Xc)を平均値・標準偏差記憶部360から読み出し、平均値の差(X−X)を演算する。更に、「+」判定信号増加量下限基準SL(+)をSLL記憶部362から読み出す。そして、平均値の差(X−Xc)と、「+」判定信号増加量下限基準SL(+)との大小関係を比較する。ステップS345で「+」判定信号増加量下限基準SL(+)より小さいと判断された場合には、「判定不能(N.D.01)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS351において、「判定不能(N.D.01)」を表示装置306及び出力装置305に出力させる。これは、有無検出用電極(作用極)550からの電気化学的信号の増加量(コントロール電極からの電気化学的信号との差分)が、「+」判定とするには充分でなく、又、「−」判定とするには小さくない、という中途半端な値を示しており、+/−判定を決することができない状態を表している。一方、ステップS345において、「+」判定信号増加量下限基準SL(+)以上の場合には、ステップS346に進む。
【0109】
(f)ステップS346では、有無判定モジュール330は、有無検出用電極(作用極)550からの電流値の平均値(X)、標準偏差(σ)、対応するコントロール電極553からの電流値の平均値(Xc)及び標準偏差(σ)を平均値・標準偏差記憶部360から、更には実効変動係数ESLLをSLL記憶部362から読み出す。有無検出用電極550からの電流の平均値Xに対して、標準偏差σの実効変動係数倍小さい値(X−ESLL*σ)と、コントロール電極からの電流平均値Xに対して、標準偏差σの実効変動係数倍大きい値(X+ESLL*σ)との大小関係を比較する。有無検出用電極550からの電流平均値Xに対して、標準偏差σの実効変動係数倍小さい値(X−ESLL*σ)が、コントロール電極からの電流平均値Xに対して、標準偏差σの実効変動係数倍大きい値(X+ESLL*σ)以上の場合、即ち、X−ESLL*σ≧X+ESLL*σの不等式が成立する場合には、充分な信号増加量が得られているとして、「+」の判定を行い、分類結果記憶部369に判定結果を分類して格納し、ステップS347において、「+」を表示装置306及び出力装置305に出力させる。一方、X−ESLL*σ≧X+ESLL*σの不等式を満足しない場合には、信号バラつきに対して、信号増加量が充分でないと判断し、「判定不能(N.D.02)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS352において、「判定不能(N.D.02)」を表示装置306及び出力装置305に出力させる。
【0110】
即ち、ステップS345において、平均値の差(X−Xc)が「+」判定信号増加量下限基準SL(+)以上の場合でも、直ちに「+」と表示せず、ステップS346において、電流増加量がバラツキに比べて充分であるか否かを判断する訳である。ここで、「+」判定信号増加量下限基準SL(+)及び「−」判定信号増加量上限基準SL(−)については、複数のサンプルに関して測定を行い、上述の手順で、信号増加量を算出し、「+」のサンプル及び「−」のサンプルに対して、それぞれ統計的な値から算出して決定するのが好ましい。より好ましくは、「+」サンプル、「−」サンプルそれぞれの信号増加量の平均値と標準偏差から、定めるのが好ましい。更により好ましくは、「+」判定信号増加量下限基準SL(+)は、「+」サンプルの「平均値−3×標準偏差」とし、「−」判定信号増加量上限基準SL(−)は「−」サンプルの「平均値+3×標準偏差」と定めるのが良い。なお、基準値を決定するのに用いるサンプル数は、それぞれ30サンプル以上であることが望ましい。
【0111】
[ステップS106その2:SNPの型判定]
図26のステップS333において、SNPの型が、野生型ホモ型か、ヘテロ型か、若しくは変異型ホモ型かを判定するSNP型判定アルゴリズムのフローに進むと決定された場合は、図28に示すフローチャートの手順で型判定モジュール340が処理を行う。なお、ここでは仮に、SNP1検出用電極551では、SNP="G"を検出,SNP2検出用電極552では、SNP="T"を検出する場合として説明を行うことにする。又、ここでは、SNP1検出用電極551は、野生型SNPを検出するためのプローブDNAが固定されており、SNP2検出用電極552には、変異型SNPを検出するためのプローブDNAが固定されているものとして、説明を行う。勿論、SNP1検出用電極を、変異型SNP検出用として、SNP2検出用電極を野生型SNP検出用として設定してもなんら問題はない。但し、判定基準パラメータの大小関係には注意が必要である。
【0112】
以下において、型判定モジュール340が、サンプルDNAについて、例えば、あるSNP位置の塩基が、G/Gホモ型であるか、G/Tヘテロ型であるか、又はT/Tホモ型であるかを判定する手順を、図28に示すフローチャートに従って説明する:
(a)先ず、型判定モジュール340は、ステップS361では、標的数が2個であるか否かを判定する。SNP="G",SNP="T"の2種を判定しようとしているのであるから、標的数が0個,1個、3個等では、そもそも初期設定が間違っているので、分類結果記憶部369に判定結果を分類して格納し、ステップS381において、表示装置306及び出力装置305に「設定エラー」の出力させる。
【0113】
(b)ステップS361で、標的数が2個であると判定された場合、型判定モジュール340は、ステップS362に進む。ステップS362では、型判定モジュール340は、SNP1検出用電極551に対応するコントロール電極(C1)553及びSNP2検出用電極552に対応するコントロール電極(C2)556からの電流値がグループ正常判定モジュール320での処理の結果、いずれも「グループ正常」であると判定された場合には、次のステップS363に進む。グループ正常判定モジュール320での処理の結果、いずれかが「グループ異常」であると判定された場合には、「判定不能(N.D.31)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS382において、「判定不能(N.D.31)」を表示装置306及び出力装置305に出力させる。
【0114】
(c)更に、ステップS363において、型判定モジュール340は、判定対象となるSNP1検出用電極551及びSNP2検出用電極552からの電流値がグループ正常判定モジュール320での処理の結果、いずれも「グループ正常」であると判定された場合には、次のステップS364に進む。グループ正常判定モジュール320での処理の結果、いずれかが「グループ異常」であると判定された場合には、「判定不能(N.D.32)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS383において、「判定不能(N.D.32)」を表示装置306及び出力装置305に出力させる。
【0115】
(d)次に、ステップS364では、型判定モジュール340は、SNP1検出用電極(SNP="G"検出用電極)551の電流値の平均値X1,SNP1検出用電極551に対応するコントロール電極(C1)553の電流値の平均値Xc1,SNP2検出用電極(SNP="T"検出用電極)552の電流値の平均値X2,SNP2検出用電極552に対応するコントロール(C2)電極556の電流値の平均値Xc2を算出する。図2では、SNP1検出用電極551とSNP2検出用電極に対応する共通のコントロール電極553を示したが、SNP1検出用電極551に対応するコントロール電極(C1)とSNP2検出用電極552に対応するコントロール(C2)とを別個に設けても構わない。算出したSNP1検出用電極551の電流値の平均値X1,SNP1検出用電極551に対応するコントロール電極(C1)553の電流値の平均値Xc1,SNP2検出用電極552の電流値の平均値X2,SNP2検出用電極552に対応するコントロール(C2)電極556の電流値の平均値Xc2は、平均値・標準偏差記憶部360に格納する。
【0116】
(e)次に、ステップS365では、型判定モジュール340は、SNP1検出用電極551の電流値の平均値X1と対応するコントロール電極(C1)553の電流値の平均値Xc1を平均値・標準偏差記憶部360から読み出し、SNP1検出用電極551の電流値の平均値X1とコントロール電極(C1)553の電流値の平均値Xc1との差(X1−Xc1)をコントロール電極(C1)553の電流値の平均値Xc1で除算した値Z(規格化されたSNP1電流増加量)を算出する:
=(X−Xc1)/Xc1 ……(5)
同様に、SNP2検出用電極552の電流値の平均値Xと対応するコントロール電極(C2)554の電流値の平均値Xc2を平均値・標準偏差記憶部360から読み出し、SNP2検出用電極552の電流値の平均値Xとコントロール電極(C2)556の電流値の平均値Xc2との差(X−Xc2)をコントロール電極(C2)556の電流値の平均値Xc2で除算した値Z(規格化されたSNP2電流増加量)を算出する:
=(X−Xc2)/Xc2 ……(6)
そして、算出したZ及びZを、規格化座標記憶部365に格納し、ステップS366に進む。
【0117】
(f)ステップS366において、型判定モジュール340は、規格化座標記憶部365から、規格化されたSNP1電流増加量Z及び規格化されたSNP2電流増加量Zを読み出す。ステップS366では、更に、規格化されたSNP1電流増加量Z及び規格化されたSNP2電流増加量Zがともに「負」の場合には、SNP1検出用電極551及びSNP2検出用電極552がともに電流増加が得られなかったことにより、「判定不能(N.D.16)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS384において、「判定不能(N.D.16)」を表示装置306及び出力装置305に出力させる。規格化されたSNP1電流増加量Z及び規格化されたSNP2電流増加量Zのいずれかでも「零以上」の場合には、次のステップS367に進む。
【0118】
(g)ステップS367において、型判定モジュール340は、規格化座標記憶部365から、規格化されたSNP1電流増加量Z及び規格化されたSNP2電流増加量Zを読み出す。ステップS367では、更に、規格化されたSNP1電流増加量ZをX座標とし、規格化されたSNP2電流増加量ZをY座標として定まる点の座標(Z,Z)を、極座標(R,A)に変換する。即ち、図29に示すように、X座標ZがY座標Zで定まる点Pに対して、Rは、原点から、点P(Z,Z)までのベクトルの距離(ベクトル長)を示し、Aは、正のX軸と、原点から点P(Z,Z)へのベクトルのなす角度を示す。ここで、ステップS366において、X座標(Z),Y座標(Z)ともに、負の場合は、判定不能として除外されているため、角度Aは、−π/2からπの間の値をとることになる。そして、ステップS367において算出されたベクトル長Rは、ベクトル長記憶部366に、又角度Aは、角度記憶部367にそれぞれ格納され、ステップS368に進む。
【0119】
(h)ステップS368において、型判定モジュール340は、ベクトル長記憶部366から、ステップS367で算出したベクトル長Rを読み出す。更に、MIR記憶部364から、電流増加割合下限基準値MIRを読み出し、ベクトル長Rと電流増加割合下限基準値MIRとの大小関係を比較する。ベクトル長Rが電流増加割合下限基準値MIRよりも小さい場合には、電流増加量が小さすぎると判断し、「判定不能(N.D.15)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS385において、「判定不能(N.D.15)」を表示装置306及び出力装置305に出力させる。ベクトル長Rが電流増加割合下限基準値MIR以上の場合には、次のステップS369に進む。ステップS369以降で、判定対象のSNP型が、野生型ホモ型か、ヘテロ型か、変異型ホモ型かを判定することになる。ここでは、SNP1検出用電極として野生型SNPを検出するためのプローブDNAが固定され、SNP2検出用電極として変異型SNPを検出するためのプローブDNAが固定されている場合で説明していることを想起して欲しい。
【0120】
(i)型判定モジュール340は、角度記憶部367からステップS367で算出した角度Aを読み出し、角度パラメータ記憶部368から、野生型角度下限値Wminを読み出し、ステップS369において、大小関係を比較する。角度Aが、野生型角度下限値Wminよりも小さい場合には、「判定不能(N.D.11)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS391において、「判定不能(N.D.11)」を表示装置306及び出力装置305に出力させる。角度Aが、野生型角度下限値Wmin以上の場合には、次のステップS370に進む。
【0121】
(j)型判定モジュール340は、角度記憶部367からステップS367で算出した角度Aを読み出し、角度パラメータ記憶部368から、野生型角度下限値Wmin及び野生型角度上限値Wmaxを読み出し、ステップS370において、大小関係を比較する。角度Aが、野生型角度下限値Wmin以上、野生型角度上限値Wmax以下の場合には、「SNP1型」(ここでは、G/Gホモ型、即ち、野生型ホモ型)として、分類結果記憶部369に判定結果を分類して格納し、ステップS392において、「G/G型」を表示装置306及び出力装置305に出力させる。角度Aが、野生型角度上限値Wmaxよりも大きい場合には、次のステップS371に進む。
【0122】
(k)型判定モジュール340は、角度記憶部367からステップS367で算出した角度Aを読み出し、角度パラメータ記憶部368から、野生型角度上限値Wmax及びヘテロ型角度下限値Hminを読み出し、ステップS371において、大小関係を比較する。角度Aが、野生型角度上限値Wmaxよりも大きく、ヘテロ型角度下限値Hminよりも小さい場合には、「判定不能(N.D.12)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS393において、「判定不能(N.D.12)」を表示装置306及び出力装置305に出力させる。角度Aが、ヘテロ型角下限値Hmin以上の場合には、次のステップS372に進む。
【0123】
(l)型判定モジュール340は、角度記憶部367からステップS367で算出した角度Aを読み出し、角度パラメータ記憶部368から、ヘテロ型角度下限値Hmin及びヘテロ型角度上限値Hmaxを読み出し、ステップS372において、大小関係を比較する。角度Aが、ヘテロ型角度下限値Hmin以上、へテロ型角度上限値Hmax以下の場合には、「SNP1/SNP2型」(ここでは、G/Tヘテロ型即ち、ヘテロ型)として、分類結果記憶部369に判定結果を分類して格納し、ステップS394において、「G/T型」を表示装置306及び出力装置305に出力させる。角度Aが、ヘテロ型角度上限値Hmaxよりも大きい場合には、次のステップS373に進む。
【0124】
(m)型判定モジュール340は、角度記憶部367からステップS367で算出した角度Aを読み出し、角度パラメータ記憶部368から、ヘテロ型角度上限値Hmax及び変異型角度下限値Mminを読み出し、ステップS373において、大小関係を比較する。角度Aが、ヘテロ型角度上限値Hmaxよりも大きく、変異型角度下限値Mminよりも小さい場合には、「判定不能(N.D.13)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS395において、「判定不能(N.D.13)」を表示装置306及び出力装置305に出力させる。角度Aが、変異型角下限値Mmin以上の場合には、次のステップS374に進む。
【0125】
(n)型判定モジュール340は、角度記憶部367からステップS367で算出した角度Aを読み出し、角度パラメータ記憶部368から、変異型角度下限値Mmin及び変異型角度上限値Mmaxを読み出し、ステップS374において、大小関係を比較する。角度Aが、変異型角度下限値Mmin以上、変異型角度上限値Mmax以下の場合には、「SNP2型」(ここでは、T/Tホモ型即ち、変異型ホモ型)として、分類結果記憶部369に判定結果を分類して格納し、ステップS396において、「T/T型」を表示装置306及び出力装置305に出力させる。角度Aが、変異型角度上限値Mmaxよりも大きい場合には、次のステップS375に進む。
【0126】
(o)型判定モジュール340は、角度記憶部367からステップS367で算出した角度Aを読み出し、角度パラメータ記憶部368から、変異型角度上限値Mmaxを読み出し、ステップS375において、大小関係を比較する。角度Aが、変異型角度上限値Mmaxよりも大きい場合には、「判定不能(N.D.14)」として、分類結果記憶部369に判定結果を分類して格納し、ステップS397において、「判定不能(N.D.14)」を表示装置306及び出力装置305に出力させる。なお、ステップS375での角度Aと変異型角度上限値Mmaxとの大小関係比較は、必ずしも必要なく、ステップS375に進んできた場合には、無条件に「判定不能(N.D.14)」としても良い。
【0127】
ここで、角度判定基準値、即ち、野生型角度下限値Wmin,野生型角度上限値Wmax,ヘテロ型角度下限値Hmin,ヘテロ型角度上限値Hmax,変異型角度下限値Mmin,変異型角度上限値Mmaxは、複数のサンプルに関して測定を行い、上述の手順で、算出した規格化された電流増加量Z、ZをそれぞれX座標、Y座標とする点を原点から見たベクトルの、正のX軸からの角度を算出し、各型(野生型、ヘテロ型、変異型)についての統計的な値から算出して決定するのが好ましい。より好ましくは、各型毎に算出された角度の平均値と標準偏差から、定めるのが好ましい。更により好ましくは、「平均値−3×標準偏差」を角度下限値とし、「平均値+3×標準偏差」を角度上限値と定めるのが良い。なお、基準値を決定するのに用いるサンプル数は、各型30サンプル以上であることが望ましい。
【0128】
電流増加割合下限基準値MIRについても、同様に複数のサンプルに関するデータから、統計的に算出して決めるのが好ましい。MIRについても、角度判定基準値と同様に、複数のサンプルに関する測定から求めたベクトル長の平均値と標準偏差を用いて決定するのが好ましい。MIRは、「平均値−3×標準偏差」と定めるのがより好ましい。なお、基準値を決定するのに用いるサンプル数は、30サンプル以上であることが望ましい。
【0129】
型判定モジュール340での型判定手順は、図31を見ると、容易に理解できる。詳細には、上記に詳述した通りであるが、概要は、角度Aが、野生型角度基準範囲(Wmin〜Wmax)の場合には「野生型」、ヘテロ型角度基準範囲(Hmin〜Hmax)の場合には「ヘテロ型」、変異型角度基準範囲(Mmin〜Mmax)の場合には「変異型」となることを意味している。
【0130】
又、本アルゴリズムをSAA1遺伝子多型に適用し、横軸を角度Aに取った際の頻度分布を図30に示している。野生型、ヘテロ型、変異型で、分布が分離しており、判定基準値として、野生型角度下限値Wminを-0.28,野生型角度上限値Wmaxを0.28,ヘテロ型角度下限値Hminを0.74,ヘテロ型角度上限値Hmaxを1.09,変異型角度下限値Mminを1.39,変異型角度上限値Mmaxを1.85と設定できることが実証されている。
【0131】
以上の説明で理解できるように、本発明の実施の形態に係る塩基配列方法によれば、チップカートリッジ11や測定系12に異常があった場合や、データにバラつきがあった場合等でも、それらの異常なデータを除外して、正確にある核酸が存在するかどうかを高い精度で判定できる。又、本発明の実施の形態に係る塩基配列方法によれば、初期設定にエラーがあれば「設定エラー」と判定し,チップ若しくは試料(バイオサンプル)に異常があれば「判定不能(サンプル・エラー)」と判定し、装置(ハードウェア)に不具合があれば「判定不能(ハードウェア・エラー)」と判定する等の処理をして、それらの異常なデータを除外できる。更に、信号のバラつきにより判定できない場合や信号強度が弱い場合等についてもその状況に応じた判定が可能である。このため、本発明の実施の形態に係る塩基配列方法によれば、種々の測定環境や状況(実状)に適格に応じながら、「野生型ホモ型」,「変異型ホモ型」,「ヘテロ型」等のように、SNPの型が、どの型であるか、又、それがホモ型であるか、ヘテロ型であるか、を高い精度で判断できる。ここで、「野生型」「ヘテロ型」「変異型」と表現したものは、具体的に塩基名を用いて”A”、”T”、”G”、”C”等の表示を行っても良い。
【0132】
(塩基配列判定プログラム)
図14,図16,図18−19,図25−28、図32に示した一連の判定操作は、図14,図16,図18−19,図25−28、図32と等価なアルゴリズムのプログラムにより、図8に示した塩基配列判定システムを制御して実行できる。このプログラムは、本発明の塩基配列判定システムを構成するコンピュータシステムのプログラム記憶装置(図示省略)に記憶させれば良い。又、このプログラムは、コンピュータ読取り可能な記録媒体に保存し、この記録媒体を塩基配列判定システムのプログラム記憶装置に読み込ませることにより、本発明の一連の判定操作を実行することができる。ここで、「コンピュータ読取り可能な記録媒体」とは、例えばコンピュータの外部メモリ装置、半導体メモリ、磁気ディスク、光ディスク、光磁気ディスク、磁気テープなどのプログラムを記録することができるような媒体などを意味する。具体的には、フレキシブルディスク、CD−ROM,MOディスク、カセットテープ、オープンリールテープ、メモリカード、ハードディスク、リムーバブルディスクなどが「コンピュータ読取り可能な記録媒体」に含まれる。
【0133】
例えば、塩基配列判定システムの本体は、フレキシブルディスク装置(フレキシブルディスクドライブ)及び光ディスク装置(光ディスクドライブ)を内蔵若しくは外部接続するように構成できる。フレキシブルディスクドライブに対してはフレキシブルディスクを、又光ディスクドライブに対してはCD−ROMをその挿入口から挿入し、所定の読み出し操作を行うことにより、これらの記録媒体に格納されたプログラムを、塩基配列判定システムを構成するプログラム記憶装置にインストールすることができる。又、所定のドライブ装置を接続することにより、例えばゲームパック等に利用されているメモリ装置としてのROMや、磁気テープ装置としてのカセットテープを用いることもできる。更に、インターネット等の情報処理ネットワークを介して、このプログラムをプログラム記憶装置に格納することが可能である。
【0134】
(その他の実施の形態)
上記のように、本発明は上記の実施の形態によって記載したが、この開示の一部をなす論述及び図面は本発明を限定するものであると理解すべきではない。この開示から当業者には様々な代替実施の形態、実施例及び運用技術が明らかとなろう。
【0135】
既に述べた上記の実施の形態の説明においては、G型、T型或いはG/T型のいずれに該当するかを判定する手法を示したが、A型,G型,C型,T型の内のいずれか2つの型、或いはそれらのヘテロの判定に適用できることは勿論である。又、上記の実施の形態の説明で理解できるように、必ずしもA型,G型,C型,T型のグループの4種類について測定データを取得する必要はなく、SNPの考えられ得る2つの塩基に関する2グループのみについて取得するのみでも良い。
【0136】
この様に、本発明はここでは記載していない様々な実施の形態等を含むことは勿論である。したがって、本発明の技術的範囲は上記の説明から妥当な特許請求の範囲に係る発明特定事項によってのみ定められるものである。
【図面の簡単な説明】
【0137】
【図1】本発明の実施の形態に係る塩基配列判定システムを構成する測定系の一例を説明する回路図である。
【図2】図1の測定系の一部をなす検出チップの構造を説明する模式的な平面図である。
【図3】図3(a)は、SNP="G"検出用電極に塩基配列GACTC……を有するプローブDNAが固定化された状態を説明する模式図で、図3(b)は、SNP="T"検出用電極に塩基配列GAATC……を有するプローブDNAが固定化された状態を説明する模式図で、図3(c)は、コントロール電極に塩基配列CAGTG……を有するプローブDNA(コントロールDNA)が固定化された状態を説明する模式図である。
【図4】本発明の実施の形態に係る塩基配列判定システムに用いられるチップカートリッジの構造の一例を説明する模式的な組み立て鳥瞰図である。
【図5】図4を上下逆にして示す、本発明の実施の形態に係るチップカートリッジの組み立て鳥瞰図である。
【図6】本発明の実施の形態に係る塩基配列判定システムを構成する送液系に備えられるバルブユニットの全体構成を示す模式図である。
【図7】本発明の実施の形態に係る塩基配列判定システムの一例を説明する概念的なブロック図である。
【図8】本発明の実施の形態に係る塩基配列判定システムを構成するコンピュータの構成の一例を説明する概念的なブロック図である。
【図9】図9(a)は、SNP1検出用電極のプローブDNAとSNP="G"のサンプルDNAとが、塩基配列が完全にマッチングしているので2本鎖を形成した状態を、図9(b)は、SNP2検出用電極のプローブDNAとSNP="G"のサンプルDNAとが2本鎖を形成できない状態を、図9(c)は、コントロール電極のコントロールDNAとSNP="G"のサンプルDNAとが2本鎖を形成できない状態を示す模式図である。
【図10】図10(a)は、SNP1検出用電極のプローブDNAとSNP="T"のサンプルDNAとが2本鎖を形成できない状態を、図10(b)は、SNP2検出用電極のプローブDNAとSNP="T"のサンプルDNAとが2本鎖を形成した状態を、図10(c)は、コントロール電極のコントロールDNAとSNP="T"のサンプルDNAとが2本鎖を形成できない状態を示す模式図である。
【図11】図11(a)は、SNP1検出用電極のプローブDNAとSNP="G/T"ヘテロのサンプルDNAとが2本鎖を形成した状態を、図11(b)は、SNP2検出用電極のプローブDNAとSNP="G/T"ヘテロのサンプルDNAとが2本鎖を形成した状態を、図11(c)は、コントロール電極のコントロールDNAとSNP="/T"ヘテロのサンプルDNAとが2本鎖を形成できない状態を示す模式図である。
【図12】図12(a)は、SNP1検出用電極のプローブDNAとSNP="G"のサンプルDNAとが、2本鎖を形成しているので、2本鎖に挿入剤が結合した状態を、図12(b)は、SNP2検出用電極のプローブDNAとSNP="G"のサンプルDNAとが2本鎖を形成できないので、2本鎖に挿入剤が結合できないが電極表面等に吸着している状態を、図12(c)は、コントロール電極のコントロールDNAとSNP="G"のサンプルDNAとが2本鎖を形成できないので、2本鎖に挿入剤が結合できないが電極表面等に吸着している状態を示す模式図である。
【図13】曲線(a)は、図12(a)に対応し、SNP1検出用電極のプローブDNAとSNP="G"のサンプルDNAとが2本鎖を形成し、2本鎖に挿入剤が結合した場合の電気化学的電流の電流−電圧特性であり、大きな電流値のピークを示す。曲線(b)は、図12(b)に対応し、SNP2検出用電極のプローブDNAとSNP="G"のサンプルDNAとが2本鎖を形成できず、2本鎖に挿入剤が結合できない場合の電流−電圧特性であり、曲線(a)に比して小さな電流値のピークを示す。又、図13の曲線(c)は、図12(c)に対応し、コントロール電極のコントロールDNAとSNP="G"のサンプルDNとが2本鎖を形成できず、2本鎖に挿入剤が結合できない場合の電流−電圧特性であり、線(b)に比して、更に小さな電流値のピークを示す。
【図14】本発明の実施の形態に係る塩基配列判定方法の全体としての流れを説明するフローチャートである。
【図15】単純移動平均法による平準化(スムージング)を説明する模式図である。
【図16】本発明の実施の形態に係る塩基配列判定方法において、各電極毎に測定された電気化学的電流の電流波形(電流−電圧特性)のそれぞれのベースラインの傾きにより、それぞれの電流波形(電流−電圧特性)の正常・異常を判定する方法の一例を説明するフローチャートである。
【図17】チップカートリッジにおける電気化学的電流の信号の例として、データ1とデータ2の2種類の電流−電圧特性を示す。データ1の電流−電圧特性が示すベースラインの傾きに比して、データ2の電流−電圧特性が示すベースラインの傾きの方が大きい。
【図18】本発明の実施の形態に係る塩基配列判定方法において、各電極毎に測定された電気化学的電流の電流波形(電流−電圧特性)から、それぞれのバックグラウンド電流を差し引いて、真の検出信号のピーク値(ピーク電流値)を求める方法の一例を説明するフローチャートである。
【図19】図18に示したフローチャートの手順に引き続き、各電極毎に測定された電流波形(電流−電圧特性)から、真の検出信号のピーク値(ピーク電流値)を求める方法の手順の流れを説明するためのフローチャートである。
【図20】本発明の実施の形態に係る塩基配列判定方法において、各電極毎に測定された電気化学的電流の電流波形(電流−電圧特性)から、ピーク位置を決定するために、電気化学的電流の微分値(di/dv)が「ゼロクロス」する点を求める方法を説明する模式図である。
【図21】本発明の実施の形態に係る塩基配列判定方法において、各電極毎に測定された電気化学的電流の電流波形(電流−電圧特性)から、それぞれのバックグラウンド電流を近似するための基準点を算出する手順を説明する模式図である。
【図22】図21の拡大図である。
【図23】本発明の実施の形態に係る塩基配列判定方法において、ゼロクロス電圧におけるピーク電流から、対応するバックグラウンド電流を差し引いて、真の検出信号(ピーク電流値)を求める方法を説明する模式図である。
【図24】図23に示す方法で求めた真の検出信号(ピーク電流値)を、複数のコントロール電極毎、複数のSNP2検出用電極(SNP="T"検出用電極)毎、複数のSNP1検出用電極(SNP="G"検出用電極)毎に整理した棒グラフである。図24(a)に示したモードAは、SNP2検出用電極で歯抜けが発生しているモードであり、図24(b)に示したモードBは、SNP2検出用電極からの信号にバラつきが大きい場合を示したモードであり、図24(c)は、ネガティブ・コントロール電流値の平均値が、ネガティブ・コントロールの上限値NCULを超えてしまった場合を概念的に示す図である。
【図25】本発明の実施の形態に係る塩基配列判定方法において、各電極毎にそれぞれ算出された一群のピーク電流値のデータに対して、グループ正常判定を行う方法の一例を説明するフローチャートである。
【図26】本発明の実施の形態に係る塩基配列判定方法において、ある核酸が存在するかどうかを判定する有無判定アルゴリズムのフローに進むのか、SNPの型が、2種の内のどちらであるか、それがホモ型であるか、ヘテロ型であるかを判定するSNP型判定アルゴリズムのフローに進むのかを決定する方法の一例を説明するフローチャートである。
【図27】本発明の実施の形態に係る塩基配列判定方法において、ある核酸が存在するかどうかを判定する有無判定アルゴリズムの一例を説明するフローチャートである。
【図28】本発明の実施の形態に係る塩基配列判定方法において、SNPの型が、野生型であるか、ヘテロ型であるか、或いは変異型であるかを判定するSNP型判定アルゴリズムの一例を説明するフローチャートである。
【図29】極座標(R,A)系における、ベクトル長R及び角度Aの定義を説明する図である。
【図30】本発明の実施の形態に係る塩基配列判定方法において、SAA1遺伝子に本方法を適用した例を示す図である。
【図31】本発明の実施の形態に係る塩基配列判定方法において、判定に関する基準となる種々の設定パラメータを示す概念イメージ図である。
【図32】本発明の実施の形態に係る塩基配列判定方法において、検査有効性を判定するアルゴリズムの一例を説明するフローチャートである。
【図33】本発明の実施の形態に係る塩基配列判定方法において、核酸の有無を判定する例について、種々の設定パラメータを示す概念イメージ図である。
【符号の説明】
【0138】
10…測定ユニット
11…チップカートリッジ
12…測定系
13…送液系
14…温度制御機構
15…制御機構
16…コンピュータ
21…検出チップ
300…CPU
301…ノイズ除去モジュール
301,…ノイズ除去モジュール
302…電流波形判定モジュール
303…バス
304…入力装置
305…出力装置
306…表示装置
310…ピーク電流検出モジュール
311…電圧値算出部
312…ベースライン近似部
313…ピーク電流値演算部
320…グループ正常判定モジュール
325…有効性判定モジュール
330…有無判定モジュール
340…型判定モジュール
351…波形判定用電圧範囲記憶部
352…許容範囲記憶部
353…ピーク探索電圧値範囲記憶部
354…ゼロクロス値記憶部
355…変曲点記憶部
356…交点電圧記憶部
357…オフセット電圧記憶部
358…ベースライン電流値記憶部
360…平均値・標準偏差記憶部
361…グループ正常判定記憶部
362…SLL記憶部
363…ESLL記憶部
364…MIR記憶部
365…規格化座標記憶部
366…ベクトル長記憶部
367…角度記憶部
368…角度パラメータ記憶部
369…分類結果記憶部
403,413,423,433、441,445,454…電磁弁
454…送液ポンプ
500…保護回路
501a,502a,503a,511a,511b,512a,512b,513a,513b…配線
502…対極
510…電圧パターン発生回路
511…トランス・インピーダンス増幅器
512…反転増幅器
513…電圧フォロア増幅器
550…有無検出用電極(作用極)
551…SNP1検出用電極(作用極)
552…SNP2検出用電極(作用極)
553…コントロール電極(作用極)
554…ポジティブ・コントロール電極(作用極)
555…ネガティブ・コントロール電極(作用極)
556…SNP2検出用電極に対応するコントロール電極(作用極)
561,562…参照極
571,572,573…プローブDNA
581,582,583…サンプルDNA
591…挿入剤
703…電気コネクタ
705…バルブユニット
707,708…ノズル
710…プローブユニット
711…カセット上蓋
712…カセット下蓋
713…パッキン
714…基板
722,723…ノズル差込孔
724,725…電気コネクタ用ポート
726…シール検出孔
727a〜727d,728a,728b,747a〜747d、748a,748b…孔
742…内表面
746…シール検出孔
749…カセット種類判別用孔
752,753…送出ポート
756,757…流路
761…電極ユニット
762,763…パッド
781,782…バルブボディ
789…カセット種類判別用ピン
Rff…フィードバック抵抗
Rf,Rw,Rs…抵抗

【特許請求の範囲】
【請求項1】
第1プローブがそれぞれ固定された複数の第1検出用電極、前記第1プローブとは塩基配列が異なる第2プローブがそれぞれ固定された複数の第2検出用電極、前記複数の第1検出用電極及び第2検出用電極に対応して配置され、前記第1プローブ及び前記第2プローブとは塩基配列が異なるコントロールDNAがそれぞれ固定された複数のコントロール電極とを備えるチップに、サンプル核酸を含む薬液を供給する段階と、
前記複数の第1検出用電極を介してそれぞれ第1検出信号を、前記複数の第2検出用電極を介してそれぞれ第2検出信号を、前記複数のコントロール電極を介してそれぞれコントロール信号を、それぞれ取得する段階と、
前記第1検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をX座標とする段階と、
前記第2検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をY座標とする段階と、
原点から前記X座標及び前記Y座標が定める点までのベクトルが正のX軸となす角度を算出する段階と、
前記角度を、予め定めた角度判定基準と比較する段階
とを含み、前記角度と前記角度判定基準との大小関係により、前記サンプル核酸の型を判定することを特徴とする塩基配列判定方法。
【請求項2】
前記第1プローブは、野生型の遺伝子型を検出するためのプローブとし、前記第2プローブは、変異型の遺伝子型を検出するためのプローブとすることを特徴とする請求項1記載の塩基配列判定方法。
【請求項3】
前記角度判定基準として、野生型下限値、野生型上限値、ヘテロ型下限値、ヘテロ型上限値、変異型下限値、変異型上限値が設定され、
前記角度が、野生型下限値と野生型上限値の間である場合に、前記サンプル核酸の型を野生型、
前記角度が、ヘテロ型下限値とヘテロ型上限値の間である場合に、前記サンプル核酸の型をヘテロ型、
前記角度が、変異型下限値と変異型上限値の間である場合に、前記サンプル核酸の型を変異型と、
それぞれ判定することを特徴とする請求項2記載の塩基配列判定方法。
【請求項4】
前記極座標変換により得られる前記ベクトルの長さを、予め定めたベクトル長判定基準と比較する段階と、
前記ベクトルの長さが前記ベクトル長判定基準よりも小さい場合には、判定不能と判断することを特徴とする請求項1記載の塩基配列判定方法。
【請求項5】
前記角度判定基準が、複数の既知の塩基配列のサンプル核酸の測定データの平均値と標準偏差を用いて定められることを特徴とする請求項3記載の塩基配列判定方法。
【請求項6】
第1プローブがそれぞれ固定された複数の第1検出用電極、前記第1プローブとは塩基配列が異なる第2プローブがそれぞれ固定された複数の第2検出用電極、前記第1プローブ及び前記第2プローブとは塩基配列が異なるコントロールDNAがそれぞれ固定された複数のコントロール電極とを備えるチップを装着して使用され、
前記複数の第1検出用電極を介してそれぞれ第1検出信号を、前記複数の第2検出用電極を介してそれぞれ第2検出信号を、前記複数のコントロール電極を介してそれぞれコントロール信号を、それぞれ取得する測定系と、
前記チップに、薬液を送液する送液系と、
前記測定系及び送液系を制御する制御機構と、
前記第1検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をX座標とし、前記第2検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をY座標とし、原点から前記X座標及び前記Y座標が定める点までのベクトルが正のX軸となす角度を算出し、前記角度を、予め定めた角度判定基準と比較し、前記角度と前記角度判定基準との大小関係により、前記サンプル核酸の型を判定する型判定モジュールを含むコンピュータ
とを備えることを特徴とする塩基配列判定システム。
【請求項7】
前記コンピュータが、
前記第1検出信号、前記第2検出信号、前記コントロール信号をそれぞれ電流−電圧特性として取得し、該電流−電圧特性のベースラインの傾きをそれぞれ求め、それぞれのベースラインの傾きにより、それぞれの電流−電圧特性の正常・異常を判定し、異常な検出信号を演算対象外とする電流波形判定モジュールと、
前記第1検出信号、前記第2検出信号、前記コントロール信号に対応するそれぞれの電流−電圧特性から、それぞれのバックグラウンド電流を差し引いて、それぞれの真の検出信号をピーク電流値として取得するピーク電流検出モジュールと、
前記ピーク電流検出モジュールが算出した真のピーク電流値のデータ群中から、正常グループを抽出するグループ正常判定モジュールと、
抽出された正常データから真のピーク電流値のデータ群を用いて、対象とする核酸の塩基が存在するかどうかを判定する有無判定モジュールもしくはサンプル核酸の一塩基多型を同定する型判定モジュール
とを更に備えることを特徴とする請求項6に記載の塩基配列判定システム。
【請求項8】
第1プローブがそれぞれ固定された複数の第1検出用電極、前記第1プローブとは塩基配列が異なる第2プローブがそれぞれ固定された複数の第2検出用電極、前記第1プローブ及び前記第2プローブとは塩基配列が異なるコントロールDNAがそれぞれ固定された複数のコントロール電極とを備えるチップにサンプル核酸を含む薬液を供給された後、
前記複数の第1検出用電極を介してそれぞれ第1検出信号を、前記複数の第2検出用電極を介してそれぞれ第2検出信号を、前記複数のコントロール電極を介してそれぞれコントロール信号を、それぞれ取得する手順と、
前記第1検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をX座標とし、前記第2検出信号の平均値から前記コントロール信号の平均値を減算し、前記コントロール信号の平均値で除算した値をY座標とし、原点から前記X座標及び前記Y座標が定める点までのベクトルが正のX軸となす角度を算出する手順と、
前記角度を、予め定めた角度判定基準と比較する手順
とを含む処理をコンピュータに実行させ、これにより前記コンピュータに、前記サンプル核酸の型が、野生型、ヘテロ型、変異型、若しくは判定不能と判断させるための塩基配列判定プログラム。

【図1】
image rotate

【図2】
image rotate

【図3】
image rotate

【図4】
image rotate

【図5】
image rotate

【図6】
image rotate

【図7】
image rotate

【図8】
image rotate

【図9】
image rotate

【図10】
image rotate

【図11】
image rotate

【図12】
image rotate

【図13】
image rotate

【図14】
image rotate

【図15】
image rotate

【図16】
image rotate

【図17】
image rotate

【図18】
image rotate

【図19】
image rotate

【図20】
image rotate

【図21】
image rotate

【図22】
image rotate

【図23】
image rotate

【図24】
image rotate

【図25】
image rotate

【図26】
image rotate

【図27】
image rotate

【図28】
image rotate

【図29】
image rotate

【図30】
image rotate

【図31】
image rotate

【図32】
image rotate

【図33】
image rotate


【公開番号】特開2009−148187(P2009−148187A)
【公開日】平成21年7月9日(2009.7.9)
【国際特許分類】
【出願番号】特願2007−327737(P2007−327737)
【出願日】平成19年12月19日(2007.12.19)
【出願人】(000003078)株式会社東芝 (54,554)
【Fターム(参考)】