アセンブリ誤り検出のための方法およびシステム（アセンブリ誤り検出）
【課題】読取りデータの再アセンブリはアセンブリにおける配列誤りを含み得る。なぜならセグメントを正確な元の順序に戻すことは難しいことがあるからである。
【解決手段】遺伝子配列アセンブリの誤りを検出するための方法は、遺伝子データの配列のアセンブリ（Ａ）を規定するステップと、読取りデータを読取りのライブラリ（Ｌ）に集めるステップと、読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットするステップと、分布（Ｄ）を被覆度Ｃで正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留するステップと、ＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集めるステップと、Ｓｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算するステップと、ユーザに対してディスプレイ上に結果を出力するステップとを含む。
【解決手段】遺伝子配列アセンブリの誤りを検出するための方法は、遺伝子データの配列のアセンブリ（Ａ）を規定するステップと、読取りデータを読取りのライブラリ（Ｌ）に集めるステップと、読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットするステップと、分布（Ｄ）を被覆度Ｃで正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留するステップと、ＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集めるステップと、Ｓｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算するステップと、ユーザに対してディスプレイ上に結果を出力するステップとを含む。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、デオキシリボ核酸（ｄｅｏｘｙｒｉｂｏｎｕｃｌｅｉｃ ａｃｉｄ：ＤＮＡ）におけるアセンブリ誤り検出、ならびにリボ核酸（Ｒｉｂｏｎｕｃｌｅｉｃ ａｃｉｄ：ＲＮＡ）における過剰発現および低発現検出に関する。
【背景技術】
【０００２】
デオキシリボ核酸（ＤＮＡ）ゲノム配列は、ＤＮＡをいくつかの塩基の配列を有するいくつかのセグメントまたは小片に分割する方法を用いて決定されることがある。各セグメント内の塩基の配列決定と、セグメントの順序の決定とをともに用いて、ＤＮＡ全体の配列を決定してもよい。セグメントの順序の決定は、バイオインフォマティクス・アセンブリ法を用いてイン・シリコ（ｉｎ−ｓｉｌｉｃｏ）で行なわれてもよい。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
読取りデータの再アセンブリはアセンブリにおける配列誤りを含み得る。なぜならセグメントを正確な元の順序に戻すことは難しいことがあるからである。
【課題を解決するための手段】
【０００４】
本発明の一局面において、遺伝子配列アセンブリの誤りを検出するための方法は、遺伝子データの配列のアセンブリ（Ａ）を規定するステップと、読取りデータを読取りのライブラリ（Ｌ）に集めるステップと、読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットするステップと、分布（Ｄ）を被覆度（ｃｏｖｅｒａｇｅ）（Ｃ）で正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留するステップと、ＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集めるステップと、Ｓｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算するステップと、ユーザに対してディスプレイ上に結果を出力するステップとを含む。
【０００５】
本発明の別の局面において、遺伝子配列の誤りを検出するためのシステムは、メモリと、ディスプレイと、遺伝子データの配列のアセンブリ（Ａ）を規定し、読取りデータを読取りのライブラリ（Ｌ）に集め、読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットし、分布（Ｄ）を被覆度（Ｃ）で正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留し、ＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集め、Ｓｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算し、ユーザに対してディスプレイ上に結果を出力するために動作するプロセッサとを含む。
【０００６】
本発明の技術によって、付加的な特徴および利点が実現される。本明細書には本発明の他の実施形態および局面が詳細に記載されており、請求される発明の一部とみなされる。利点および特徴を伴う本発明をより良く理解するために、説明および図面を参照されたい。
【０００７】
本発明であるとみなされる主題は、請求項において特定的に示され明確に主張されている。本発明の前述およびその他の特徴および利点は、以下の詳細な説明および添付の図面から明らかである。
【図面の簡単な説明】
【０００８】
【図１】複数のＤＮＡ配列およびそれらの配列のセグメントへの分割を示す図である。
【図２】配列中の誤りを定めるためのシステム２００の例示的な実施形態を示す図である。
【図３】図２のシステムによって行なわれ得る例示的な処理方法を示すブロック図である。
【図４】図２のシステムによって行なわれ得る例示的な処理方法を示すブロック図である。
【図５】読取りの頻度のヒストグラムを示す図である。
【発明を実施するための形態】
【０００９】
デオキシリボ核酸（ＤＮＡ）ゲノム配列は、たとえば圧縮空気装置（噴霧器）または制限酵素などを用いることによって、ＤＮＡをいくつかの塩基の配列を有するいくつかのセグメントまたは小片に分割することによって決定されることがある。図１は、複数の類似のＤＮＡ配列、およびその配列のセグメントへの分割を示している。これに関して、いくつかの類似のＤＮＡ鎖１０２（例、５０個またはそれ以上の鎖）が、たとえば５０塩基から５００塩基などの範囲のいくつかの塩基１０６を有する複数のセグメント１０４に分裂または切断されてもよい。セグメント１０４は必ずしも等しい長さに切断されるわけではない。セグメント１０４が切断されると、セグメント１０４を読取って塩基１０６が識別され、識別された塩基１０６の各セグメントにおける位置が定められることによって、各セグメント１０４に対する読取りデータがもたらされる。代替的には、セグメントの端部（例、各端部から１００塩基）を読取って塩基を識別してもよい。セグメントの読取りは、たとえばヌクレオチドの蛍光ラベリングおよび高分解能レーザ・イメージングを含む、合成による配列決定（ｓｅｑｕｅｎｃｉｎｇ−ｂｙ−ｓｙｎｔｈｅｓｉｓ）のプロセスなどによって行なわれてもよい。結果として生じるデータは複数の読取りを含み、各読取りは塩基１０６と、各セグメント１０４におけるその塩基１０６の位置とを識別する。読取りデータはグループ化されて、特定の長さにおける読取りの頻度（すなわち特定の長さの塩基を有する読取りの数）を含む読取りのライブラリ（Ｌ）に入れられる。被覆度（Ｃ）とは、配列決定されたＤＮＡにおける位置が重複しているセグメント１０４の平均コピー数である。配列決定されたセグメント１０４の長さに加えてＤＮＡ配列の長さが既知であれば、被覆度Ｃが分かる。ＤＮＡゲノム配列の長さが未知であるときには、ユーザが推定の長さを与えてもよい。読取りデータの「再アセンブリ」を行なうことによって、ＤＮＡゲノム配列の一部または全体を表わすアセンブリ（Ａ）データをもたらしてもよい。アセンブリは、たとえばアセンブラ（イン・シリコのバイオインフォマティクス・ツール）などを用いて、読取りにおける塩基間の重複を考慮し、可能なところで重複した読取りを連結することによって行なわれてもよい。アセンブリ・データは、所与の位置ｉにおける読取りカウントｃｉおよび読取り長さｌを含むベクトルＶ＝＜ｉ，ｃｉ，ｌ１，ｌ２，…，ｌｃｉ＞を含む。ベクトルの一例はＶ＝＜３４，３，１０，１２，１０２＞を含み、これは位置３４がそれぞれ長さ１０、１２および１０２の３つの読取りによって重複することを示す。読取りデータの再アセンブリはアセンブリにおける配列誤りを含み得る。なぜならセグメントを正確な元の順序に戻すことは難しいことがあるからである。以下に説明される例示的な方法およびシステムは、アセンブリにおける誤りの検出を改善するものである。
【００１０】
これに関して、図２は配列中の誤りを定めるためのシステム２００の例示的な実施形態を示す。示される実施形態は、ディスプレイ装置２０４と、入力装置２０６と、読取りデータ２０１およびアセンブリ２０３を保存するメモリ２０８とに通信的に接続されたプロセッサ２０２を含む。
【００１１】
図３および図４は、システム２００によって行なわれ得る例示的な処理方法のブロック図を示す。図３を参照すると、ブロック３０２において読取りデータを含むアセンブリ（Ａ）が規定される。ブロック３０４において、読取りデータが読取りのライブラリ（Ｌ）に集められる。ブロック３０６において、Ｌから読取りのサイズ対サイズ当りの読取り数のヒストグラムがプロットされる。ヒストグラムの一例を図５に示す。ブロック３０８において、被覆度Ｃを用いて分布Ｄを正規化して（Ｄ’）を得、ここでＤ’はＬの期待標準分布であり、平均μおよび標準偏差σを有する。正規化は、Ａに対する被覆度Ｃを用い、（ユーザによって与えられる上側および下側のカットオフを用いて）被覆度Ｃを表わすとは考えられないベクトルＶをフィルタリングして除去することによって行なわれる。前のステップの出力を用いてライブラリが再び計算される。Ｄ’を得るために用いられない位置（ｉ）が保留される。ブロック３１０において、アセンブリＡ内の各位置（ｉ）に対して、位置ｉが重複している読取りの部分集合Ｓｉ⊂ＬがベクトルＶｉに集められる。ブロック３１２において、Ｓｉから平均（μｉ）および標準偏差（√ｃｉ・σｉ）が算出される。（図４の）ブロック３１４において、ライブラリのμからのμｉの偏差が計算される。ブロック３１６において、ライブラリのσからの（√ｃｉ・σｉ）の偏差が定められる。ブロック３１８において、閾値を用いて、μｉおよび（√ｃｉ・σｉ）における異常な偏差（すなわち閾値の外側の偏差）が定められる。
【００１２】
ブロック３２０において、ユーザ分析のために結果がディスプレイ装置に出力されてもよい。アセンブリにおける各位置ｉに対して、平均（μｉ）が期待値から所与の閾値よりも多く逸脱しているか、または標準偏差（√ｃｉ・σｉ）が所与の閾値よりも大きいとき、その位置ｉはアセンブリの誤りの可能性があるとしてフラグを立てられる。次いでユーザは、別の方法によってデータを再アセンブリするか、追加の読取りを生成して再アセンブリするか、または配列情報の代替的供給源を用いることによって、これらのフラグを立てられた領域におけるアセンブリの間違いの可能性を修正することに集中できる。
【００１３】
ＲＮＡデータに対しても同様のプロセスを用いることができるが、フラグを立てられる位置は過剰発現または低発現に関連付けられる。
【００１４】
本明細書において用いられる用語は特定の実施形態を説明する目的のためのみのものであって、本発明を限定することは意図されない。本明細書において用いられる単数形「ａ」、「ａｎ」および「ｔｈｅ」は、状況が明らかに別様を示していない限り、複数形をも含むことが意図される。さらに、「含む（ｃｏｍｐｒｉｓｅｓ）」もしくは「含む（ｃｏｍｐｒｉｓｉｎｇ）」という用語またはその両方が本明細書において用いられるとき、それは述べられる特徴、完全体、ステップ、動作、構成要素もしくは成分またはその組合わせの存在を特定するが、１つまたはそれ以上の他の特徴、完全体、ステップ、動作、構成要素成分もしくはそのグループまたはその組合わせの存在または追加を排除するものではないことが理解されるだろう。
【００１５】
以下の請求項におけるすべての手段またはステップ・プラス機能（ｍｅａｎｓ ｏｒ ｓｔｅｐ ｐｌｕｓ ｆｕｎｃｔｉｏｎ）要素に対応する構造、材料、動作、および同等物は、特定的に請求される他の請求要素と組合わせてその機能を行なうためのあらゆる構造、材料または動作を含むことが意図される。本発明の説明は例示および説明の目的のために提供されたものであるが、網羅的になったり、開示される形に本発明を制限したりすることは意図されない。本発明の範囲および趣旨から逸脱することなく、通常の当業者には多くの修正および変更が明らかになるだろう。実施形態は、本発明の原理および実際の適用を最も良く説明し、他の通常の当業者が予期される特定の使用に好適であるようなさまざまな修正を伴うさまざまな実施形態に対して本発明を理解できるようにするために選択されて記載されたものである。
【００１６】
本明細書に示される図面は単なる一例である。本発明の趣旨から逸脱することなく、この図面または本明細書に記載されるステップ（もしくは動作）には多くの変更形が存在するだろう。たとえば、これらのステップが異なる順序で行なわれてもよいし、ステップが追加、削除または変更されてもよい。これらの変更形はすべて、請求される本発明の一部であるとみなされる。
【００１７】
本発明に対する好ましい実施形態を説明したが、現在および将来にわたり、当業者が以下の請求項の範囲内のさまざまな改善および強化を行ない得ることが理解されるだろう。これらの請求項は、最初に記載された本発明に対する適切な保護を維持するものと解釈されるべきである。
【符号の説明】
【００１８】
３０２ 読取りデータからアセンブリ（Ａ）を規定する
３０４ 読取りデータを読取りのライブラリ（Ｌ）に集める
３０６ Ｌから読取りサイズ対サイズ当りの読取り数のヒストグラムをプロットする
３０８ 被覆度Ｃを用いてＤを正規化してＤ’を得、Ｄ’を得るために用いられない位置ｉを保留する
３１０ ＡおよびＤ’を用いて読取りの部分集合Ｓｉ⊂Ｌを集める
３１２ 平均μｉおよび標準偏差（√ｃｉ・σｉ）を計算する
【技術分野】
【０００１】
本発明は、デオキシリボ核酸（ｄｅｏｘｙｒｉｂｏｎｕｃｌｅｉｃ ａｃｉｄ：ＤＮＡ）におけるアセンブリ誤り検出、ならびにリボ核酸（Ｒｉｂｏｎｕｃｌｅｉｃ ａｃｉｄ：ＲＮＡ）における過剰発現および低発現検出に関する。
【背景技術】
【０００２】
デオキシリボ核酸（ＤＮＡ）ゲノム配列は、ＤＮＡをいくつかの塩基の配列を有するいくつかのセグメントまたは小片に分割する方法を用いて決定されることがある。各セグメント内の塩基の配列決定と、セグメントの順序の決定とをともに用いて、ＤＮＡ全体の配列を決定してもよい。セグメントの順序の決定は、バイオインフォマティクス・アセンブリ法を用いてイン・シリコ（ｉｎ−ｓｉｌｉｃｏ）で行なわれてもよい。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
読取りデータの再アセンブリはアセンブリにおける配列誤りを含み得る。なぜならセグメントを正確な元の順序に戻すことは難しいことがあるからである。
【課題を解決するための手段】
【０００４】
本発明の一局面において、遺伝子配列アセンブリの誤りを検出するための方法は、遺伝子データの配列のアセンブリ（Ａ）を規定するステップと、読取りデータを読取りのライブラリ（Ｌ）に集めるステップと、読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットするステップと、分布（Ｄ）を被覆度（ｃｏｖｅｒａｇｅ）（Ｃ）で正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留するステップと、ＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集めるステップと、Ｓｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算するステップと、ユーザに対してディスプレイ上に結果を出力するステップとを含む。
【０００５】
本発明の別の局面において、遺伝子配列の誤りを検出するためのシステムは、メモリと、ディスプレイと、遺伝子データの配列のアセンブリ（Ａ）を規定し、読取りデータを読取りのライブラリ（Ｌ）に集め、読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットし、分布（Ｄ）を被覆度（Ｃ）で正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留し、ＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集め、Ｓｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算し、ユーザに対してディスプレイ上に結果を出力するために動作するプロセッサとを含む。
【０００６】
本発明の技術によって、付加的な特徴および利点が実現される。本明細書には本発明の他の実施形態および局面が詳細に記載されており、請求される発明の一部とみなされる。利点および特徴を伴う本発明をより良く理解するために、説明および図面を参照されたい。
【０００７】
本発明であるとみなされる主題は、請求項において特定的に示され明確に主張されている。本発明の前述およびその他の特徴および利点は、以下の詳細な説明および添付の図面から明らかである。
【図面の簡単な説明】
【０００８】
【図１】複数のＤＮＡ配列およびそれらの配列のセグメントへの分割を示す図である。
【図２】配列中の誤りを定めるためのシステム２００の例示的な実施形態を示す図である。
【図３】図２のシステムによって行なわれ得る例示的な処理方法を示すブロック図である。
【図４】図２のシステムによって行なわれ得る例示的な処理方法を示すブロック図である。
【図５】読取りの頻度のヒストグラムを示す図である。
【発明を実施するための形態】
【０００９】
デオキシリボ核酸（ＤＮＡ）ゲノム配列は、たとえば圧縮空気装置（噴霧器）または制限酵素などを用いることによって、ＤＮＡをいくつかの塩基の配列を有するいくつかのセグメントまたは小片に分割することによって決定されることがある。図１は、複数の類似のＤＮＡ配列、およびその配列のセグメントへの分割を示している。これに関して、いくつかの類似のＤＮＡ鎖１０２（例、５０個またはそれ以上の鎖）が、たとえば５０塩基から５００塩基などの範囲のいくつかの塩基１０６を有する複数のセグメント１０４に分裂または切断されてもよい。セグメント１０４は必ずしも等しい長さに切断されるわけではない。セグメント１０４が切断されると、セグメント１０４を読取って塩基１０６が識別され、識別された塩基１０６の各セグメントにおける位置が定められることによって、各セグメント１０４に対する読取りデータがもたらされる。代替的には、セグメントの端部（例、各端部から１００塩基）を読取って塩基を識別してもよい。セグメントの読取りは、たとえばヌクレオチドの蛍光ラベリングおよび高分解能レーザ・イメージングを含む、合成による配列決定（ｓｅｑｕｅｎｃｉｎｇ−ｂｙ−ｓｙｎｔｈｅｓｉｓ）のプロセスなどによって行なわれてもよい。結果として生じるデータは複数の読取りを含み、各読取りは塩基１０６と、各セグメント１０４におけるその塩基１０６の位置とを識別する。読取りデータはグループ化されて、特定の長さにおける読取りの頻度（すなわち特定の長さの塩基を有する読取りの数）を含む読取りのライブラリ（Ｌ）に入れられる。被覆度（Ｃ）とは、配列決定されたＤＮＡにおける位置が重複しているセグメント１０４の平均コピー数である。配列決定されたセグメント１０４の長さに加えてＤＮＡ配列の長さが既知であれば、被覆度Ｃが分かる。ＤＮＡゲノム配列の長さが未知であるときには、ユーザが推定の長さを与えてもよい。読取りデータの「再アセンブリ」を行なうことによって、ＤＮＡゲノム配列の一部または全体を表わすアセンブリ（Ａ）データをもたらしてもよい。アセンブリは、たとえばアセンブラ（イン・シリコのバイオインフォマティクス・ツール）などを用いて、読取りにおける塩基間の重複を考慮し、可能なところで重複した読取りを連結することによって行なわれてもよい。アセンブリ・データは、所与の位置ｉにおける読取りカウントｃｉおよび読取り長さｌを含むベクトルＶ＝＜ｉ，ｃｉ，ｌ１，ｌ２，…，ｌｃｉ＞を含む。ベクトルの一例はＶ＝＜３４，３，１０，１２，１０２＞を含み、これは位置３４がそれぞれ長さ１０、１２および１０２の３つの読取りによって重複することを示す。読取りデータの再アセンブリはアセンブリにおける配列誤りを含み得る。なぜならセグメントを正確な元の順序に戻すことは難しいことがあるからである。以下に説明される例示的な方法およびシステムは、アセンブリにおける誤りの検出を改善するものである。
【００１０】
これに関して、図２は配列中の誤りを定めるためのシステム２００の例示的な実施形態を示す。示される実施形態は、ディスプレイ装置２０４と、入力装置２０６と、読取りデータ２０１およびアセンブリ２０３を保存するメモリ２０８とに通信的に接続されたプロセッサ２０２を含む。
【００１１】
図３および図４は、システム２００によって行なわれ得る例示的な処理方法のブロック図を示す。図３を参照すると、ブロック３０２において読取りデータを含むアセンブリ（Ａ）が規定される。ブロック３０４において、読取りデータが読取りのライブラリ（Ｌ）に集められる。ブロック３０６において、Ｌから読取りのサイズ対サイズ当りの読取り数のヒストグラムがプロットされる。ヒストグラムの一例を図５に示す。ブロック３０８において、被覆度Ｃを用いて分布Ｄを正規化して（Ｄ’）を得、ここでＤ’はＬの期待標準分布であり、平均μおよび標準偏差σを有する。正規化は、Ａに対する被覆度Ｃを用い、（ユーザによって与えられる上側および下側のカットオフを用いて）被覆度Ｃを表わすとは考えられないベクトルＶをフィルタリングして除去することによって行なわれる。前のステップの出力を用いてライブラリが再び計算される。Ｄ’を得るために用いられない位置（ｉ）が保留される。ブロック３１０において、アセンブリＡ内の各位置（ｉ）に対して、位置ｉが重複している読取りの部分集合Ｓｉ⊂ＬがベクトルＶｉに集められる。ブロック３１２において、Ｓｉから平均（μｉ）および標準偏差（√ｃｉ・σｉ）が算出される。（図４の）ブロック３１４において、ライブラリのμからのμｉの偏差が計算される。ブロック３１６において、ライブラリのσからの（√ｃｉ・σｉ）の偏差が定められる。ブロック３１８において、閾値を用いて、μｉおよび（√ｃｉ・σｉ）における異常な偏差（すなわち閾値の外側の偏差）が定められる。
【００１２】
ブロック３２０において、ユーザ分析のために結果がディスプレイ装置に出力されてもよい。アセンブリにおける各位置ｉに対して、平均（μｉ）が期待値から所与の閾値よりも多く逸脱しているか、または標準偏差（√ｃｉ・σｉ）が所与の閾値よりも大きいとき、その位置ｉはアセンブリの誤りの可能性があるとしてフラグを立てられる。次いでユーザは、別の方法によってデータを再アセンブリするか、追加の読取りを生成して再アセンブリするか、または配列情報の代替的供給源を用いることによって、これらのフラグを立てられた領域におけるアセンブリの間違いの可能性を修正することに集中できる。
【００１３】
ＲＮＡデータに対しても同様のプロセスを用いることができるが、フラグを立てられる位置は過剰発現または低発現に関連付けられる。
【００１４】
本明細書において用いられる用語は特定の実施形態を説明する目的のためのみのものであって、本発明を限定することは意図されない。本明細書において用いられる単数形「ａ」、「ａｎ」および「ｔｈｅ」は、状況が明らかに別様を示していない限り、複数形をも含むことが意図される。さらに、「含む（ｃｏｍｐｒｉｓｅｓ）」もしくは「含む（ｃｏｍｐｒｉｓｉｎｇ）」という用語またはその両方が本明細書において用いられるとき、それは述べられる特徴、完全体、ステップ、動作、構成要素もしくは成分またはその組合わせの存在を特定するが、１つまたはそれ以上の他の特徴、完全体、ステップ、動作、構成要素成分もしくはそのグループまたはその組合わせの存在または追加を排除するものではないことが理解されるだろう。
【００１５】
以下の請求項におけるすべての手段またはステップ・プラス機能（ｍｅａｎｓ ｏｒ ｓｔｅｐ ｐｌｕｓ ｆｕｎｃｔｉｏｎ）要素に対応する構造、材料、動作、および同等物は、特定的に請求される他の請求要素と組合わせてその機能を行なうためのあらゆる構造、材料または動作を含むことが意図される。本発明の説明は例示および説明の目的のために提供されたものであるが、網羅的になったり、開示される形に本発明を制限したりすることは意図されない。本発明の範囲および趣旨から逸脱することなく、通常の当業者には多くの修正および変更が明らかになるだろう。実施形態は、本発明の原理および実際の適用を最も良く説明し、他の通常の当業者が予期される特定の使用に好適であるようなさまざまな修正を伴うさまざまな実施形態に対して本発明を理解できるようにするために選択されて記載されたものである。
【００１６】
本明細書に示される図面は単なる一例である。本発明の趣旨から逸脱することなく、この図面または本明細書に記載されるステップ（もしくは動作）には多くの変更形が存在するだろう。たとえば、これらのステップが異なる順序で行なわれてもよいし、ステップが追加、削除または変更されてもよい。これらの変更形はすべて、請求される本発明の一部であるとみなされる。
【００１７】
本発明に対する好ましい実施形態を説明したが、現在および将来にわたり、当業者が以下の請求項の範囲内のさまざまな改善および強化を行ない得ることが理解されるだろう。これらの請求項は、最初に記載された本発明に対する適切な保護を維持するものと解釈されるべきである。
【符号の説明】
【００１８】
３０２ 読取りデータからアセンブリ（Ａ）を規定する
３０４ 読取りデータを読取りのライブラリ（Ｌ）に集める
３０６ Ｌから読取りサイズ対サイズ当りの読取り数のヒストグラムをプロットする
３０８ 被覆度Ｃを用いてＤを正規化してＤ’を得、Ｄ’を得るために用いられない位置ｉを保留する
３１０ ＡおよびＤ’を用いて読取りの部分集合Ｓｉ⊂Ｌを集める
３１２ 平均μｉおよび標準偏差（√ｃｉ・σｉ）を計算する
【特許請求の範囲】
【請求項１】
コンピュータ・システムにより遺伝子配列アセンブリの誤りを検出するための方法であって、前記コンピュータ・システムはメモリとディスプレイとプロセッサとを含み、
前記方法は、
前記プロセッサが前記メモリ上で遺伝子データの配列のアセンブリ（Ａ）を規定するステップと、
前記プロセッサが前記メモリ上で読取りデータを読取りのライブラリ（Ｌ）に集めるステップと、
前記プロセッサが読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットするステップと、
前記プロセッサが分布（Ｄ）を被覆度Ｃで正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留するステップと、
前記プロセッサがＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集めるステップと、
前記プロセッサがＳｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算するステップと、
前記プロセッサがユーザに対して前記ディスプレイ上に結果を出力するステップと
を含む、方法。
【請求項２】
前記方法は、前記プロセッサが読取りの前記ライブラリから各位置（ｉ）に対するμからのμｉの偏差を計算するステップをさらに含む、請求項１に記載の方法。
【請求項３】
前記方法は、前記プロセッサが読取りの前記ライブラリから各位置（ｉ）に対するσからの√ｃｉ・σｉの偏差を定めるステップをさらに含む、請求項１に記載の方法。
【請求項４】
前記方法は、前記偏差を閾値と比較して、前記閾値よりも大きいかまたは小さい偏差を識別するステップをさらに含む、請求項２に記載の方法。
【請求項５】
前記方法は、前記プロセッサが前記偏差を閾値と比較して、前記閾値よりも大きいかまたは小さい偏差を識別するステップをさらに含む、請求項３に記載の方法。
【請求項６】
前記方法は、前記プロセッサが前記識別された偏差の位置ｉをユーザに対して前記ディスプレイ上に出力させるステップを含む、請求項４に記載の方法。
【請求項７】
前記方法は、前記プロセッサが前記識別された偏差の位置ｉをユーザに対して前記ディスプレイ上に出力させるステップを含む、請求項５に記載の方法。
【請求項８】
前記アセンブリは、配列アセンブリのためのイン・シリコのバイオインフォマティクス法によって規定される、請求項１に記載の方法。
【請求項９】
前記読取りデータは、デオキシリボ核酸（ＤＮＡ）のセグメント内の複数の塩基の位置および識別子を含む、請求項１に記載の方法。
【請求項１０】
読取りの前記ライブラリは複数の読取りデータを含む、請求項１に記載の方法。
【請求項１１】
遺伝子配列の誤りを検出するためのシステムであって、前記システムは、
メモリと、
ディスプレイと、
遺伝子データの配列のアセンブリ（Ａ）を規定し、読取りデータを読取りのライブラリ（Ｌ）に集め、読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットし、分布（Ｄ）を被覆度Ｃで正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留し、ＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集め、Ｓｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算し、ユーザに対して前記ディスプレイ上に結果を出力するために動作するプロセッサと
を含む、システム。
【請求項１２】
前記プロセッサは、読取りの前記ライブラリから各位置（ｉ）に対するμからのμｉの偏差を計算するためにさらに動作する、請求項１１に記載のシステム。
【請求項１３】
前記プロセッサは、読取りの前記ライブラリから各位置（ｉ）に対するσからの√ｃｉ・σｉの偏差を定めるためにさらに動作する、請求項１１に記載のシステム。
【請求項１４】
前記プロセッサは、前記偏差を閾値と比較して、前記閾値よりも大きいかまたは小さい偏差を識別するためにさらに動作する、請求項１２に記載のシステム。
【請求項１５】
前記プロセッサは、前記偏差を閾値と比較して、前記閾値よりも大きいかまたは小さい偏差を識別するためにさらに動作する、請求項１３に記載のシステム。
【請求項１６】
前記プロセッサは、前記識別された偏差の位置ｉをユーザに対して前記ディスプレイ上に出力するためにさらに動作する、請求項１４に記載のシステム。
【請求項１７】
前記プロセッサは、前記識別された偏差の位置ｉをユーザに対して前記ディスプレイ上に出力するためにさらに動作する、請求項１５に記載のシステム。
【請求項１８】
前記アセンブリは、配列アセンブリのためのイン・シリコのバイオインフォマティクス法によって規定される、請求項１１に記載のシステム。
【請求項１９】
前記読取りデータは、デオキシリボ核酸（ＤＮＡ）のセグメント内の複数の塩基の位置および識別子を含む、請求項１１に記載のシステム。
【請求項２０】
読取りの前記ライブラリは複数の読取りデータを含む、請求項１１に記載のシステム。
【請求項１】
コンピュータ・システムにより遺伝子配列アセンブリの誤りを検出するための方法であって、前記コンピュータ・システムはメモリとディスプレイとプロセッサとを含み、
前記方法は、
前記プロセッサが前記メモリ上で遺伝子データの配列のアセンブリ（Ａ）を規定するステップと、
前記プロセッサが前記メモリ上で読取りデータを読取りのライブラリ（Ｌ）に集めるステップと、
前記プロセッサが読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットするステップと、
前記プロセッサが分布（Ｄ）を被覆度Ｃで正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留するステップと、
前記プロセッサがＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集めるステップと、
前記プロセッサがＳｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算するステップと、
前記プロセッサがユーザに対して前記ディスプレイ上に結果を出力するステップと
を含む、方法。
【請求項２】
前記方法は、前記プロセッサが読取りの前記ライブラリから各位置（ｉ）に対するμからのμｉの偏差を計算するステップをさらに含む、請求項１に記載の方法。
【請求項３】
前記方法は、前記プロセッサが読取りの前記ライブラリから各位置（ｉ）に対するσからの√ｃｉ・σｉの偏差を定めるステップをさらに含む、請求項１に記載の方法。
【請求項４】
前記方法は、前記偏差を閾値と比較して、前記閾値よりも大きいかまたは小さい偏差を識別するステップをさらに含む、請求項２に記載の方法。
【請求項５】
前記方法は、前記プロセッサが前記偏差を閾値と比較して、前記閾値よりも大きいかまたは小さい偏差を識別するステップをさらに含む、請求項３に記載の方法。
【請求項６】
前記方法は、前記プロセッサが前記識別された偏差の位置ｉをユーザに対して前記ディスプレイ上に出力させるステップを含む、請求項４に記載の方法。
【請求項７】
前記方法は、前記プロセッサが前記識別された偏差の位置ｉをユーザに対して前記ディスプレイ上に出力させるステップを含む、請求項５に記載の方法。
【請求項８】
前記アセンブリは、配列アセンブリのためのイン・シリコのバイオインフォマティクス法によって規定される、請求項１に記載の方法。
【請求項９】
前記読取りデータは、デオキシリボ核酸（ＤＮＡ）のセグメント内の複数の塩基の位置および識別子を含む、請求項１に記載の方法。
【請求項１０】
読取りの前記ライブラリは複数の読取りデータを含む、請求項１に記載の方法。
【請求項１１】
遺伝子配列の誤りを検出するためのシステムであって、前記システムは、
メモリと、
ディスプレイと、
遺伝子データの配列のアセンブリ（Ａ）を規定し、読取りデータを読取りのライブラリ（Ｌ）に集め、読取りのサイズ対サイズ当りの読取り数のヒストグラムをプロットし、分布（Ｄ）を被覆度Ｃで正規化することによって、平均（μ）および標準偏差（σ）を有するＤ’を得て、Ｄ’を得るために用いられない位置（ｉ）を保留し、ＡおよびＤ’を用いて読取りの部分集合（Ｓｉ⊂Ｌ）を集め、Ｓｉを用いて平均（μｉ）および標準偏差（√ｃｉ・σｉ）を計算し、ユーザに対して前記ディスプレイ上に結果を出力するために動作するプロセッサと
を含む、システム。
【請求項１２】
前記プロセッサは、読取りの前記ライブラリから各位置（ｉ）に対するμからのμｉの偏差を計算するためにさらに動作する、請求項１１に記載のシステム。
【請求項１３】
前記プロセッサは、読取りの前記ライブラリから各位置（ｉ）に対するσからの√ｃｉ・σｉの偏差を定めるためにさらに動作する、請求項１１に記載のシステム。
【請求項１４】
前記プロセッサは、前記偏差を閾値と比較して、前記閾値よりも大きいかまたは小さい偏差を識別するためにさらに動作する、請求項１２に記載のシステム。
【請求項１５】
前記プロセッサは、前記偏差を閾値と比較して、前記閾値よりも大きいかまたは小さい偏差を識別するためにさらに動作する、請求項１３に記載のシステム。
【請求項１６】
前記プロセッサは、前記識別された偏差の位置ｉをユーザに対して前記ディスプレイ上に出力するためにさらに動作する、請求項１４に記載のシステム。
【請求項１７】
前記プロセッサは、前記識別された偏差の位置ｉをユーザに対して前記ディスプレイ上に出力するためにさらに動作する、請求項１５に記載のシステム。
【請求項１８】
前記アセンブリは、配列アセンブリのためのイン・シリコのバイオインフォマティクス法によって規定される、請求項１１に記載のシステム。
【請求項１９】
前記読取りデータは、デオキシリボ核酸（ＤＮＡ）のセグメント内の複数の塩基の位置および識別子を含む、請求項１１に記載のシステム。
【請求項２０】
読取りの前記ライブラリは複数の読取りデータを含む、請求項１１に記載のシステム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図２】
【図３】
【図４】
【図５】
【公開番号】特開２０１２−１５５７１５（Ｐ２０１２−１５５７１５Ａ）
【公開日】平成２４年８月１６日（２０１２．８．１６）
【国際特許分類】
【出願番号】特願２０１２−７７６４（Ｐ２０１２−７７６４）
【出願日】平成２４年１月１８日（２０１２．１．１８）
【出願人】（３９０００９５３１）インターナショナル・ビジネス・マシーンズ・コーポレーション (4,084)
【氏名又は名称原語表記】ＩＮＴＥＲＮＡＴＩＯＮＡＬ ＢＵＳＩＮＥＳＳ ＭＡＳＣＨＩＮＥＳ ＣＯＲＰＯＲＡＴＩＯＮ
【Ｆターム（参考）】
【公開日】平成２４年８月１６日（２０１２．８．１６）
【国際特許分類】
【出願日】平成２４年１月１８日（２０１２．１．１８）
【出願人】（３９０００９５３１）インターナショナル・ビジネス・マシーンズ・コーポレーション (4,084)
【氏名又は名称原語表記】ＩＮＴＥＲＮＡＴＩＯＮＡＬ ＢＵＳＩＮＥＳＳ ＭＡＳＣＨＩＮＥＳ ＣＯＲＰＯＲＡＴＩＯＮ
【Ｆターム（参考）】
[ Back to top ]