JPS63191959A - 核酸の塩基配列決定のための信号処理方法 - Google Patents

核酸の塩基配列決定のための信号処理方法

Info

Publication number
JPS63191959A
JPS63191959A JP62025010A JP2501087A JPS63191959A JP S63191959 A JPS63191959 A JP S63191959A JP 62025010 A JP62025010 A JP 62025010A JP 2501087 A JP2501087 A JP 2501087A JP S63191959 A JPS63191959 A JP S63191959A
Authority
JP
Japan
Prior art keywords
peak
image data
autoradiograph
signal processing
base
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.)
Granted
Application number
JP62025010A
Other languages
English (en)
Other versions
JPH0664058B2 (ja
Inventor
Makoto Hara
誠 原
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP62025010A priority Critical patent/JPH0664058B2/ja
Priority to US07/030,062 priority patent/US4885696A/en
Priority to DE8787104502T priority patent/DE3781307T2/de
Priority to DE3751366T priority patent/DE3751366T2/de
Priority to EP91118244A priority patent/EP0476712B1/en
Priority to EP87104502A priority patent/EP0242629B1/en
Publication of JPS63191959A publication Critical patent/JPS63191959A/ja
Publication of JPH0664058B2 publication Critical patent/JPH0664058B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Investigating Or Analysing Biological Materials (AREA)
  • Measurement Of Radiation (AREA)
  • Saccharide Compounds (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 [発明の分野] 本発明は、核酸の塩基配列決定のための信号処理方法に
関するものである。
[発明の背景] 近年、急速に発達して来た分子生物学の分野においては
、生物体の機能や複製のメカニズムを解明するために、
生物体のもつ遺伝情報を明らかにすることが必須のこと
となっている。とりわけ、特定の遺伝情報を担うDNA
 (もしくはDNA断片物、以下同様)などの核酸の塩
基配列を決定することが必要不可欠なこととなっている
DNA、RNAなとの核酸の塩基配列を決定するための
代表的な方法として、オートラジオグラフィーを利用す
るマキサム・ギルバー) (14aza■−Gilbe
rt )法およびサンガー・クールンン(Sanger
−C:oulson)法が知られている。前者のマキサ
ム・ギルバート法は、まず、塩基配列を決定しようとし
ているDNAあるいはDNA断片物の鎖状分子の一方の
端部に32 p等の放射性同位元素を含む基を結合させ
ることにより、その対象物を放射性標識物質としたのち
、化学的な手段を利用して鎖状分子の各構成単位間の結
合を塩基特異的に切断する0次に、この操作により得ら
れた塩基特異的DNA切断分解物の混合物をゲル電気泳
動法により分離展開し、多数の切断分解物がそれぞれ分
離展開されて形成された分離展開パターン(ただし、視
覚的には見ることができない)を得る。この分離展開パ
ターンをたとえばX線フィルム上に可視化してそのオー
トラジオグラフを得、得られたオートラジオグラフと各
々の塩基特異的切断手段とから、放射性同位元素が結合
された鎖状分子の端部から一定の位置関係にある塩基を
順次決定し、これにより対象物全ての塩基配列を決定す
ることができる。
また、後者のサンガーやクールラン法は、DNAあるい
はDNA断片物の鎖状分子と相補的であって、かつ放射
性標識が付与されたDNA合成物を化学的な手段を利用
して塩基特異的に合成し、この塩基特異的DNA合成物
の混合物を用いて上記と同様にしてそのオートラジオグ
ラフから塩基配列を決定する方法である。
本出願人は、上記核酸の塩基配列決定を簡易かつ高精度
で行なうことを目的として、それに利用されるオートラ
ジオグラフ測定操作において、上記X線フィルム等の写
真感光材料を用いる従来の放射線写真法の代りに、蓄積
性蛍光体シートを用いる放射線像変換方法を利用する方
法について既に特許出願している(特開昭59−830
57号、特願昭58−201231号)、ここで、蓄積
性蛍光体シートは輝尽性蛍光体からなるものであり、放
射線エネルギーを該蛍光体シートの輝尽性蛍光体に吸収
させたのち、可視乃至赤外領域の電磁波(励起光)で励
起することにより、放射線エネルギーを蛍光として放出
させることができるものである。この方法によれば、露
光時間を大幅に短縮化することができ、また従来より問
題となっていた化学カブリ等が発生することがない、さ
らに、放射性標識物質のオートラジオグラフは、一旦放
射線エネルギーとして蛍光体シートに蓄積されたのち輝
尽光として光電的に読み出されるから、直接にデジタル
信号として得たのち適当な記録媒体に保存することがで
きる。
従来より、核酸の塩基配列は、可視化されたオートラジ
オグラフについて、放射性標識が付与された核酸の塩基
特異的切断分解物もしくは塩基特異的合成物(以下、単
に核酸の塩基特異的断片物と称する)の分離展開位δ(
バンド)を視覚的に判断し、それらバンドの位置を相互
に比較することにより決定されている。よって、オート
ラジオグラフの解析は通常人間の視覚を通して行なわれ
ており、そのために多大な時間と労力が費されている。
また、人間の目に依存しているため、オートラジオグラ
フを解析して得られる核酸の塩基配列が解析者によって
異なるなど塩基配列情報の精度には限界がある。
そこで1本出願人は、上記オートラジオグラフをデジタ
ル画像データとして得た模このデジタル画像データに連
出なデジタル信号処理を施すことにより、DNAの塩基
配列を自動的に決定する方法についても既に特許出願し
ている(特開昭59−126527号、特願昭59−8
9615号、特願昭60−226091号、特願昭60
−226092号等)、オートラジオグラフに対応する
デジタル画像データは、従来の放射線フィルムを利用す
る場合には一旦オートラジオグラフを該フィルム上に可
視画像化したのち、反射光または透過光を利用して光電
的に読み取ることにより得られる。また、蓄積性蛍光体
シートを用いる場合には、□オートラジオグラフが蓄積
記録された蛍光体シートを直接に読み出すことにより得
られる。
しかしながら、実際に放射性標識物質を電気泳動法など
により支持媒体上に分離展開させて得られた分離展開パ
ターンには種々の歪みおよびノイズが生じがちである。
たとえば、試料の作成時において核酸の塩基特異的断片
物の調製、分離が不十分であったり、試料に放射性不純
物が混入したり、あるいは露光時に自然放射線等の照射
を受けたりすることにより、そのオートラジオグラフ上
にノイズが発生することがある。この結果、ノイズも含
めてバンドの比較同定が行なわれるために、塩基配列決
定に誤差が生じて得られる情報の精度が低下してしまう
このようなノイズが発生した場合であっても、そのオー
トラジオグラフを包含するデジタル画像データに効率良
く信号処理を施して核酸の塩基配列を高精度で自動決定
することが望まれる。
[発明の要旨] 本発明者は、放射性標識物質からなる分離展開パターン
のオートラジオグラフにノイズが生じている場合であっ
てもそのデジタル画像データに好適なデジタル信号処理
を施すことにより、核酸の塩基配列を簡易かつ高精度で
自動決定することを実現した。
すなわち、本発明は、放射性標識が付与された塩基特異
的DNA断片物もしくは塩基特異的RNA断片物の混合
物が支持媒体上に一次元的方向に分S展開されて形成さ
れた分離展開パターンのオートラジオグラフを包含する
デジタル画像データについてデジタル信号処理を行なう
ことにより、核酸の塩基配列を決定する方法において、
1)一つの分離展開列について、分離展開方向に沿った
位置と信号のレベルとからなる一次元波形を少なくとも
二つ得る工程、 2)全ての一次元波形上でピークを検出する工程、 3)一つの一次元波形上の一つのピークについて、隣接
する一次元波形上の該ピークの位置を中心とする一定範
囲内に、ピークが存在するか否かを探索する工程、 4)第三工程においてピークが存在する場合には、順次
隣接する一次元波形に対して第三工程を繰り返す工程、 5)第三および第四工程で、互いに隣接する複数個の一
次元波形上でピークが検出された場合には、一連のピー
クの位置にバンドが存在すると決定し、それ以外の場合
には該ピークの位置にバンドは存在しないと決定する工
程、および6)上記第三乃至第五工程を順次繰り返すこ
とにより、分離展開列上の全てのバンドを検出する工程
、 を含むことを特徴とする核酸の塩基配列決定のための信
号処理方法を提供するものである。
本発明によれば、核酸の塩基特異的断片物の混合物を支
持媒体上で分離展開させて得られた分離U13t!パタ
ーンのオートラジオグラフを包含するデジタル画像デー
タにおいて、分離展開パターンにノイズが生じている場
合でも、ノイズを排除して真正バンドのみを検出しうる
信号処理機能を有する処理回路を通すことにより、核酸
の塩基配列を簡易かつ高精度で得ることができる。
分離展開列についてその分離展開方向に沿った位置と信
号のレベルとからなる一次元波形(ラスタ)上で単純に
ピークを検出しただけでは、真のバンドのほかにノイズ
をもバンドと誤認しがちであり、正確な塩基配列情報を
得ることができない。通常、バンドはスロットの大きざ
に応じて分離展開方向に垂直な方向に幅を持っている。
一方、試料の作成過程における放射性同位元素の不完全
な導入あるいは露光過程における自然放射線の混入等が
原因となってオートラジオグラフ上に現われるノイズは
スポット状であるか、あるいはバンド状であるとしても
真のバンドよりもその幅が狭い。
本発明によれば、一つの分離展開列について多数のラス
タを作成し、分離展開方向に垂直な方向にこれらのラス
ク上でピークが連続的に存在するかどうかを判断基準と
して、検出されたピークが真のバンドを示すピークであ
るのかあるいはノイズであるのかを決定することにより
、ノイズを排除して真正バンドのみを検出することがで
きる。
そして、検出されたバンドの位置に基づいて分離展開列
間でバンドを比較することにより、核酸の塩基配列を簡
便かつ高精度に決定することができる。
[発明の構成] 本発明において用いられる試料の例としては、放射性標
識が付与されたDNA、RNA等の核酸の塩基特異的断
片物の混合物を挙げることができる。ここで、核酸の断
片物とは長鎖状の分子の一部分を意味する。たとえば、
塩基特異的DNA断片物混合物の一種である塩基特異的
DNA切断分解物混合物は、前述のマキサム・ギルバー
ト法に従って、放射性標識が付与されたDNAを塩基特
異的に切断分解することにより得られる。
また、塩基特異的DNA合成物混合物は前述のサンガー
・クールラン法に従って、DNAをテンプレート(鋳型
)として、放射性標識が付与されたデオキシヌクレオシ
ドトリフオスフェートとDNA合成酵素とを用いて合成
することにより得られる。
さらに、塩基特異的RNA断片物の混合物も上記と同様
の方法により、切断分解物混合物としてまたは合成物混
合物として得ることができる。なお、DNAはその構成
単位としてアデニン、グアニン、チミン、シトシンの四
種類の塩基からなるが、一方RNAはアデニン、グアニ
ン、ウラシル、シトシンの四種類の塩基からなる。
放射性標識は、これらの物質に適当な方法で32p、+
4q、418S、3 H,12sIなどの放射性同位元
素を保持させることによって付与される。
試料である放射性標識が付与された核酸の塩基特異的断
片物の混合物はゲル状支持媒体など公知の各種の支持媒
体を用いて、電気泳動法、薄層クロマトグラフィー、カ
ラムクロマトグラフィー、ペーパーグロマトグラフィー
など種々の分離展開方法により支持媒体上に分#展開さ
れる。
次に、放射性標識物質が分離展開された支持媒体につい
て、従来の写真感光材料を用いる放射線写真法により、
あるいは蓄積性蛍光体シートを用いる放射線像変換方法
により是のオートラジオグラフが得られ、次いで適当な
読取り(読出し)系を介して放射性標識物質のオートラ
ジオグラフを包含するデジタル画像データが得られる。
前者の放射線写真法を利用する場合には、まず支持媒体
とX&litフィルム等の写真感光材料とを低温もしく
は常温で髪時間(数時間〜数十時間)重ね合わせて放射
線フィルムを感光させたのち、現像して放射性標識物質
のオートラジオグラフを放射線フィルム上に可視画像化
する。次いで、画像読取装置を用いて放射線フィルム上
に可視化されたオートラジオグラフを読み取る。たとえ
ば、放射線フィルムに光ビームを照射してその透過光ま
たは反射光を光電的に検出することにより、オートラジ
オグラフは電気信号として得られる。さらに、この電気
信号なA/D変換することにより、オートラジオグラフ
を包含するデジタル画像データを得ることができる。
後者の放射線像変換方法を利用する場合には、まず、支
持媒体と蓄積性蛍光体シートとを常温で短時間(数秒〜
数十分間)重ね合わせて蛍光体シートに放射性標識物質
から放出される放射線エネルギーを蓄積させることによ
り、そのオートラジオグラフを蛍光体シートに一種の潜
像として記録する。ここで、蓄積性蛍光体シートは、た
とえばプラスチックフィルムからなる支持体、二価ユー
ロピウム賦活弗化臭化バリウム(BaFBr:Eu2ゝ
)等の輝尽性蛍光体からなる蛍光体層、および透明な保
護膜がこの順に積層されたものである。蓄積性蛍光体シ
ートに含有されている輝尽性蛍光体は、X線等の放射線
が照射されるとその放射線エネルギーを吸収して蓄積し
、そののち可視乃至赤外領域の光で励起すると蓄積して
いた放射線エネルギーを輝尽光として放出するという特
性を有する。
次いで、読出装置を用いて蓄積性蛍光体シートに蓄積記
録されたオートラジオグラフを読み出す。具体的には、
たとえば蛍光体シートをレーザー光で走査して放射線エ
ネルギーを輝尽光として放出させ、この輝尽光を光電的
に検出することにより、放射性標識物質の十−トラジオ
グラフは可視画像化することなく直接に電気信号として
得られる。さらに、この電気信号をA/D変換すること
により、オートラジオグラフを包含するデジタル画像デ
ータを得ることができる。
上述のオートラジオグラフ測定操作およびオートラジオ
グラフを包含するデジタル画像データを得る方法の詳細
については、前記特開昭59−83057号、特開昭5
9−126527号、特開昭59−126278号等の
各公報に記載されている。
なお、上記においては、支持媒体上に分離展開された放
射性標識物質のオートラジオグラフを包含するデジタル
画像データを得る方法として、従来の放射線写真法およ
び放射線像変換方法を利用する方法について述べたが、
これらの方法に限定されるものではなく、それ以外の如
何なる方法により得られたデジタル画像データであって
も放射性標識物質のオートラジオグラフを包含するもの
である限り、本発明の信号処理方法を適用することが可
1走である。
また、上記いずれの方法においてもオートラジオグラフ
の読取り(または読出し)は、放射線フィルム(または
蓄積性蛍光体シー]・)の全面に亘って行なう必要はな
く、画像領域のみについて行なうことも勿論可能である
さらに、予め各分離展開列の位置およびバンドの幅等に
ついての情報を入力して読取り(読出し)条件を設定し
ておき、読取り−(読出し)操作においては各バンド上
を二本以上の走査線が通過するような走査線密度で光ビ
ームによる走査を行なうことにより、読取り(読出し)
時間を短縮化して必要な情報を効率良く得ることができ
る。なお1本発明においてオートラジオグラフを包含す
るデジタル画像データには、このようにして得られたデ
ジタル信号をも含まれる。
得られたデジタル画像データは、放射線フィルム(また
は蛍光体シート)に固定された座標系で表わされた座標
(x、y)とその座標における信号のレベル(2)とか
らなるデジタル信号oxyの集合であり、一つの信号は
一つの画素に対応している。信号のレベルはその座標に
おける画像濃度、すなわち放射性標識物質の量を表わし
ている。従って、デジタル画像データは放射性標識物質
の二次元的な位置情報を有している。
このようにして得られた支持媒体上の放射性標識物質の
オートラジオグラフを包含するデジタル画像データには
、以下に述べるような本発明の方法により信号処理が施
されて、目的の核酸の塩基配列の決定が行なわれる。
本発明の信号処理方法の実施の態様を、次の四種類の放
射性標識が付与された塩基特異的DNA断片物の組合せ
により形成された泳動列(分離展開列)からなる場合に
ついて説明する。
1)グアニン(G)特異的DNA合成物2)アデニン(
A)特異的DNA合成物3)チミン(T)特異的DNA
合成物 4)シトシン(C)特異的DNA合成物ここで、各塩基
特異的DNA合成物は、塩基特異的に切断分解もしくは
合成された、すなわち末端の塩基を同じくする種々の長
さのDNA断片物からなる。
第1図は、上記四種類の塩基特異的DNA断片物がそれ
ぞれ四つのスロットに電気泳動されてなル泳動パターン
のオートラジオグラフの例を示す、泳動方向はX方向で
ある。
このオートラジオグラフを包含するデジタル画像データ
は、信号処理回路において一旦メモリ(バッファーメモ
リ、または磁気ディスク等の不揮発性メモリ)に記憶さ
れる。
通常、放射線フィルムもしくは蓄積性蛍光体シートを読
み取ることにより得られたデジタル画像データは、泳動
パターンだけではなくフィルムもしくは蛍光体シート全
面についての画像情報を含んでおり、画像領域(パター
ン領域)以外の余分な部分のデジタルデータを含んでい
る。従って、本発明の塩基配列決定のための信号処理に
先立って、該処理の質的向上および処理効率の向上のた
めには、得られたデジタル画像データから画像領域に対
応する画像データのみを抽出すること(いわゆる「画像
領域の切出しJ)が望ましい。
たとえば、泳動パターンに対応する画像データは、次の
ような信号処理により抽出することができる。まず、得
られた画像データを、泳動方向に垂直な方向(X方向)
に沿って二つ以上のブロックに分割したのち、ブロック
ごとに分離展開方向(X方向)に沿って画像データを加
算することによりX方向の位置と信号のレベルとからな
る一次元波形(射影)を得る0次いで、ブロックごとに
その射影から泳動パターン領域の境界、更には各泳動列
(レーン)の境界を検出し、ブロック間でこれらの境界
を補間することにより、@方向(X方向)における泳動
パターン領域全体およびレーン全体の境界を決定するこ
とができる。同様にして、長さ方向(X方向)における
境界も決定することができる。
なお、上記画像領域の切出しのための信号処理方法の詳
細については、本出願人による昭和61年12月27日
出願の特願昭61−    号明細書に記載されている
あるいは、得られたデジタル画像データに基づいてオー
トラジオグラフをCRT等の画面に電気的に表示したの
ち、画面上でカーソル、マウス等を用いて操作すること
によりあるいはキーボード操作によって、泳動パターン
領域は情報として入力することができる。
なお、長さ方向における境界は試料の注入位置(スロッ
ト)および泳動条件等によっておおよその位置が判明し
ているから、予め固定座標として設定されてもよい。
まず第一に、泳動パターン領域に対応する画像データに
基づいて、各泳動夕1(レーン)について泳動方向に沿
った位置と信号のレベルとからなる一次元波形(ラスタ
)を二つ以上作成する。
デジタル画像デー”夕が、前記のように各バンドについ
て少なくとも二本の走査線がかかるような走査線密度で
泳動方向に沿って走査することによって読み取られたも
のである場合には、直接に各走査線について位!(y)
と信号のレベル(2)とからなるラスタを作成すること
ができる。また、オートラジオグラフ全面に渡って読み
取られたものである場合には、画像データ上で上記と同
様の走査を行なうことにより、各レーンに沿って信号を
抽出した後ラスタを作成する。
一つのレーン当りのラスタの数はスロットの形状(幅)
、試料の量および画素数などによっても異なるが、一般
にバンド検出の精度の点からは5〜40個の範囲にある
。−例として、スロットの幅が5mm、画素の間隔が0
.2mmであるとき、ラスタ間隔も0.2mm(画素単
位でラスタ作成)とするとレーン当りのラスタ数は25
個である。
第2図は、第一レーンについて作成された多数のラスタ
を部分的に示す、ラスタは換言すれば、第一レーンをX
方向の各位置で切断したときの濃度断面図に相当する。
第二に、全てのラスク上でピークを検出する。
ピークはたとえば、信号のレベルの差分の符号が反転す
る点を求めることにより検出することができる。
第2図において、矢印(↓)はそれぞれピークを示す。
第三に、一つのラスク上の一つのピークに着目して、隣
接するラスク上の該ピークの位置を中心とする一定範囲
内に、ピークが存在するか否かを探索する。次いで第四
に、ピークが探し出される限り、順次隣接するラスク上
でピーク探索を繰り返し行なう。
ここで、ピークの位置を中心とする一定範囲とは、通常
はピーク位置Vaを中心としてYa±σmmの範囲であ
り、σはラスタ間隔に依存して決定され゛る。−例とし
て、ラスタ間隔が0.2mmであるときσ=0.4であ
る。ピーク探索の範囲を限定することにより、複数のラ
スタ」二で泳動方向に大きなずれのないピークが存在す
る場合にのみバンドと判断され、バンド検出の精度を高
めることができる。
また、ピーク探索は、隣接する一つのラスタについて行
なってもよいし、あるいは同時に二つ以りのラスタにつ
いて行なってもよい、たとえば、2〜5個のラスタにつ
いて同時に探索することができる。この複数のラスタに
ついての同時探索は後記するピークの途切れがある場合
に好適な方法である。
第3図は、第2図の各ラスタについて検索されたピーク
の存在を示す。X印はピーク位置を表わす。
第五に、検索されたピークの連続性を判定し、この判定
に基づいて一定範囲(ya±σ)内にある一連のピーク
がバンドであるか否かを決定する。
ピークの連続性の判定は、バンドの幅(スロットの幅)
、パターンの形状、ラスタ間隔などによっても異なるが
、一般にN個以上のラスタでビーりが連続して検出され
た場合に、これら一連のピークの位置にバンドが存在す
ると決定することができる。ここで、Nは、スロットの
幅が5mmおよびラスタ間隔が0.2mmである場合(
レーン当りのラスク数:25個)には、バンドの幅が電
気泳動の際に伸縮することを考慮して通常は13≦N≦
37の範囲の整数である。
好ましくは、N個以上M個以下のラスタでピークが連続
して検出された場合にバンドが存在すると決定する。」
−記条件の場合には、たとえばN=15およびM= 3
0とすることができる。
さらに、ピークがなまっていたり、信号レベルが局部的
に低くなっている場合には、一連のピークが必ずしもラ
スク間で連続して一検出されないことがある。従って、
ピーク探索がなされた一連のラスタにおいて、ピークの
途切れ(ピークの検出されないラスタの存在)が一定範
囲内であれば、これもバンドと判断決定するのが望まし
い。たとえば、ピークの途切れがL箇所以下であって、
かつ各途切れた箇所においてピークの検出されないラス
タの数かに個以下であると3にも、該一連のビーク位置
にバンドが存在すると決定する。ここで、KおよびLは
それぞれラスタ間隔に依存して決定され、たとえばラス
タ間隔が0.2mmである場合には1≦に≦3および1
≦L≦3の範囲の整数である。
そして、上記以外の場合にはピークの位置にバンドは存
在しないと決定する。
第4図は、複数のラスタについて検索されたピークの存
在を示す別の例である(×印:ピーク位置)。
第4図において、一連のピークAおよびEはそれぞれバ
ンドであり、残りのピークB、CおよびDはバンドでは
なくノイズであると決定される。
以上の操作を全てのピークについて行ない、泳動パター
ン上の全てのバンドを検出する。このようにして、泳動
パターンにノイズが生じている場合であっても、高精度
でバンドを検出することができる。
なお、上記においては画像領域の切出しを泳動パターン
単位で行なった場合について説明したが、始めにレーン
単位で画像領域を切り出し、切り出されたレーンごとに
上述の操作を繰り返し行なってもよい。
なお、泳動パターンにノイズの発生のほかに、スマイリ
ング現象、オフセット歪みあるいはバンドの融合などの
種々の歪みが発生している場合には、上記バンドの検出
後にこれらの補正のための信号処理を行なってもよい。
ここで、スマイリング現象は、支持媒体の中央部のスロ
ットの泳動距離に比べて両端部のスロットの泳動距離が
短くなる現象であり、泳動過程における放熱効果(いわ
ゆるエツジ効果)などが原因となって生じる。オフセッ
ト歪みとは、レーン間相互の全体的な位置ズレをいい、
スロットの形状の相違等により試料の電気泳動の開始位
置、開始時間が各スロットで異なることなどが原因とな
って生じる。また、バンドの融合は、泳動が十分でない
ために、二乃至三個のバンドが連結して一個の幅広なバ
ンドを形成していることをいう、一般にパターン上部の
泳動開始位置に近い領域で発生しやすい。
信号処理によるこれらの補正の詳細については、本出願
人による特願昭60−74899号、特願昭60−74
900号、特願昭60−85275号、特願昭60−8
5276号、特願昭60−111186号、特願昭60
−111187号の各明細書に記載されている。
DNAの塩基配列は検出されたバンドの位置を相互に比
較し、バンドに序列を付けることにより直ちに決定する
ことができる。このとき、上記四種類の塩基特異的DN
A断片物の組合せが排他的な組合せであることから、異
なるレーン上の同じ位置に二つ以」二のバンドは存在し
えないことを利用して、容易に序列を決定することがで
きる。上記(1)〜(4)のスロットはそれぞれ(G)
、(A)、(T)、(C)からなる末端塩基についての
情報を有するから、各バンドの属するスロットに対応す
る塩基で置換することにより、DNAの塩基配列(例え
ばA−G−C−T−A−A−G−・・・)を得ることが
できる。
このようにして、DNAの片方の鎖状分子についての塩
基配列を決定することができる。なお、DNAの塩基配
列についての情報は、上記の表示形態に限られるもので
はなく、たとえば所望により同時に各バンドの強度(2
′)を放射性標識物質の相対量として表示することも可
能である。さらに、I)NAの二本の鎖状分子両方につ
いての塩基配列を表示することもできる。
あるいはまた、DNAの塩基配列情報は、上記の信号処
理がなされたデジタル画像データに基づいて画像として
表示してもよい。すなわち、検出された各バンドの位置
をオリジナルのオートラジオグラフとともに可視画像化
して表示することができる。この場合には、最終的な塩
基配列決定を解析者自身がこの表示画像に基づいて行な
うことが回部である。
なお、上記においては、試料である塩基特異的DNA断
片物の8合物とl、て(G、A、T、C)の排他的組合
せを利用した場合について説明したが、本発明の信号処
理方法はこの組合せに限定されるものではなく、例えば
(G、G+A、T+C,C)などの種々の組合せに適用
することができる。また同様に、塩基特異的RNA断片
物の混合物(例えば、G、A、U、Cの組合せ)につい
ても本発明の方法を適用することができる。さらに、バ
ンドの検出は一組の核酸の塩基特異的断片物の分離展開
列に限定されるものではなく、支持媒体上に同時に分離
展開された全ての分離展開列について行なうことが可能
である。
このようにして得られた塩基配列情報についてはこのほ
かにも、たとえば、既に記録保存されている他の核酸の
塩基配列と照合するなどの遺伝言語学的情報処理を行な
うことも可能である。
上述の信号処理により決定された核酸の塩基配列につい
ての情報は、信号処理回路から出力されたのち、次いで
直接的に、もしくは必要により磁気ディスクや磁気テー
プなどの記憶保存手段を介して記録装置に伝送される。
記録装置としては、たとえば、感光材料上をし一ザー光
等で走査して光学的に記録するもの、CRT等に表示さ
れた記号−数値をビデオ拳プリンター等に記録するもの
、熱線を用いて感熱記録材料上に記録するものなど種々
の原理に基づいた記録装置を用いることができる。
【図面の簡単な説明】
第1図は、ノイズが生じている泳動パターンのオートラ
ジオグラフの例を示す図である。 第2図は、第一レーンについての複数のラスタを示す部
分図である。 第3図は、第2図に示されたラスタについて検索された
ピークの存在を示す部分図である。 第4図は、別の複数のラスタについて検索されたピーク
の存在を示す部分図である。 矢印(↓):ビーク、×印:ビーク位置、特許出願人 
 富士写真フィルム株式会社代 理 人  弁理士  
柳 川 泰 男第1図 (+)  (2)  (3)  (4)略、=、   
・ e 第3図 ABC:DE

Claims (1)

  1. 【特許請求の範囲】 1、放射性標識が付与された塩基特異的DNA断片物も
    しくは塩基特異的RNA断片物の混合物が支持媒体上に
    一次元的方向に分離展開されて形成された分離展開パタ
    ーンのオートラジオグラフを包含するデジタル画像デー
    タについてデジタル信号処理を行なうことにより、核酸
    の塩基配列を決定する方法において、 1)一つの分離展開列について、分離展開方向に沿った
    位置と信号のレベルとからなる一次元波形を少なくとも
    二つ得る工程、 2)全ての一次元波形上でピークを検出する工程、 3)一つの一次元波形上の一つのピークについて、隣接
    する一次元波形上の該ピークの位置を中心とする一定範
    囲内に、ピークが存在するか否かを探索する工程、 4)第三工程においてピークが存在する場合には、順次
    隣接する一次元波形に対して第三工程を繰り返す工程、 5)第三および第四工程で、互いに隣接する複数個の一
    次元波形上でピークが検出された場合には、一連のピー
    クの位置にバンドが存在すると決定し、それ以外の場合
    には該ピークの位置にバンドは存在しないと決定する工
    程、および 6)上記第三乃至第五工程を順次繰り返すことにより、
    分離展開列上の全てのバンドを検出する工程、 を含むことを特徴とする核酸の塩基配列決定のための信
    号処理方法。 2、上記第一工程の前に、デジタル画像データから分離
    展開パターン領域に対応する画像データを抽出すること
    により、該パターン領域に対応する画像データについて
    第一工程乃至第六工程を行なうことを特徴とする特許請
    求の範囲第1項記載の核酸の塩基配列決定のための信号
    処理方法。 3、上記第一工程の前に、デジタル画像データから各分
    離展開列に対応する画像データを抽出することにより、
    分離展開列ごとに該分離展開列に対応する画像データに
    ついて第一工程乃至第六工程を繰り返し行なうことを特
    徴とする特許請求の範囲第1項記載の核酸の塩基配列決
    定のための信号処理方法。 4、上記塩基特異的DNA断片物の混合物が、 (1)グアニン特異的DNA断片物、 (2)アデニン特異的DNA断片物、 (3)チミン特異的DNA断片物、 (4)シトシン特異的DNA断片物、 の四種類からなり、分離展開パターンが、これら四種類
    の塩基特異的DNA断片物がそれぞれ支持媒体上に分離
    展開されて形成された四列の分離展開列からなることを
    特徴とする特許請求の範囲第1項記載の信号処理方法。 5、上記オートラジオグラフを包含するデジタル画像デ
    ータが、支持媒体と輝尽性蛍光体を含有する蓄積性蛍光
    体シートとを重ね合わせて、支持媒体上の放射性標識物
    質のオートラジオグラフを該蛍光体シートに蓄積記録し
    たのち、該蛍光体シートに励起光を照射して該オートラ
    ジオグラフを輝尽光として光電的に読み出すことにより
    得られたものであることを特徴とする特許請求の範囲第
    1項記載の信号処理方法。 6、上記オートラジオグラフを包含するデジタル画像デ
    ータが、支持媒体と写真感光材料とを重ね合わせて、支
    持媒体上の放射性標識物質のオートラジオグラフを該感
    光材料に感光記録したのち、該感光材料上に可視化され
    たオートラジオグラフを光電的に読み取ることにより得
    られたものであることを特徴とする特許請求の範囲第1
    項記載の信号処理方法。
JP62025010A 1986-03-26 1987-02-04 核酸の塩基配列決定のための信号処理方法 Expired - Fee Related JPH0664058B2 (ja)

Priority Applications (6)

Application Number Priority Date Filing Date Title
JP62025010A JPH0664058B2 (ja) 1987-02-04 1987-02-04 核酸の塩基配列決定のための信号処理方法
US07/030,062 US4885696A (en) 1986-03-26 1987-03-25 Signal processing method for determining base sequence of nucleic acid
DE8787104502T DE3781307T2 (de) 1986-03-26 1987-03-26 Signalverarbeitungsmethode zur bestimmung des grundaufbaus von nukleinsaeuren.
DE3751366T DE3751366T2 (de) 1986-03-26 1987-03-26 Signalverarbeitungsmethode zur Bestimmung des Grundaufbaus von Nukleinsäuren.
EP91118244A EP0476712B1 (en) 1986-03-26 1987-03-26 Signal processing method for determining base sequence of nucleic acid
EP87104502A EP0242629B1 (en) 1986-03-26 1987-03-26 Signal processing method for determining base sequence of nucleic acid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP62025010A JPH0664058B2 (ja) 1987-02-04 1987-02-04 核酸の塩基配列決定のための信号処理方法

Publications (2)

Publication Number Publication Date
JPS63191959A true JPS63191959A (ja) 1988-08-09
JPH0664058B2 JPH0664058B2 (ja) 1994-08-22

Family

ID=12153953

Family Applications (1)

Application Number Title Priority Date Filing Date
JP62025010A Expired - Fee Related JPH0664058B2 (ja) 1986-03-26 1987-02-04 核酸の塩基配列決定のための信号処理方法

Country Status (1)

Country Link
JP (1) JPH0664058B2 (ja)

Also Published As

Publication number Publication date
JPH0664058B2 (ja) 1994-08-22

Similar Documents

Publication Publication Date Title
EP0476712B1 (en) Signal processing method for determining base sequence of nucleic acid
US4802101A (en) Signal processing method for determining base sequence of nucleic acid
US4972325A (en) Signal processing method for determining base sequence of nucleic acid
JPS63169588A (ja) オ−トラジオグラフ解析のための信号処理方法
JPH0462345B2 (ja)
JPH0570792B2 (ja)
JPH0467917B2 (ja)
JPS63191959A (ja) 核酸の塩基配列決定のための信号処理方法
JPS62225956A (ja) 核酸の塩基配列決定のための信号処理方法
JPH0467630B2 (ja)
JPH0529068B2 (ja)
JPS61243361A (ja) 核酸の塩基配列決定のための信号処理方法
JPH0529067B2 (ja)
JPS61269066A (ja) 核酸の塩基配列決定のための信号処理方法
JPH0467629B2 (ja)
JPS61243362A (ja) 核酸の塩基配列決定のための信号処理方法
JPS61269065A (ja) 核酸の塩基配列決定のための信号処理方法
JPH0533747B2 (ja)
JPH0528785B2 (ja)
JPH0529069B2 (ja)
JPH0529072B2 (ja)
JPH0570791B2 (ja)
JPH0513474B2 (ja)
JPS62126348A (ja) 核酸の塩基配列決定のための信号処理方法

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees