JP2019508781A - 染色体異常の検出 - Google Patents

染色体異常の検出 Download PDF

Info

Publication number
JP2019508781A
JP2019508781A JP2018533694A JP2018533694A JP2019508781A JP 2019508781 A JP2019508781 A JP 2019508781A JP 2018533694 A JP2018533694 A JP 2018533694A JP 2018533694 A JP2018533694 A JP 2018533694A JP 2019508781 A JP2019508781 A JP 2019508781A
Authority
JP
Japan
Prior art keywords
parameter
chromosome
determining
fetal
probability
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
JP2018533694A
Other languages
English (en)
Other versions
JP7079433B2 (ja
JP2019508781A5 (ja
Inventor
チャールズ フォアマン、マシュー
チャールズ フォアマン、マシュー
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.)
Premaitha Ltd
Original Assignee
Premaitha 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 Premaitha Ltd filed Critical Premaitha Ltd
Publication of JP2019508781A publication Critical patent/JP2019508781A/ja
Publication of JP2019508781A5 publication Critical patent/JP2019508781A5/ja
Application granted granted Critical
Publication of JP7079433B2 publication Critical patent/JP7079433B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6813Hybridisation assays
    • C12Q1/6827Hybridisation assays for detection of mutation or polymorphism
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6813Hybridisation assays
    • C12Q1/6841In situ hybridisation
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Organic Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Engineering & Computer Science (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Genetics & Genomics (AREA)
  • Analytical Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Physics & Mathematics (AREA)
  • Biotechnology (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • General Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • Microbiology (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Medical Informatics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Apparatus Associated With Microorganisms And Enzymes (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

【課題】胎児の染色体の異常の確率を判定するコンピュータ実装方法を提供すること。【解決手段】この方法は、女性被験者から取得された生体試料から標的染色体に関する第1のパラメータおよび染色体配列密度を表す第2のパラメータを示すデータを決定することと、胎児の染色体の異常を示す尤度比を決定することであって、この尤度比は、第1のパラメータおよび第2のパラメータに基づいてそれぞれの異常モデルおよび正常モデルに従って染色体の異常の確率と染色体の正常の確率の比として決定することと、1つまたは複数の性能パラメータ閾値を決定することと、試料に関連付けられた1つまたは複数の性能パラメータの推定値を1つまたは複数の性能パラメータ閾値と比較することとを含む。【選択図】図2

Description

本発明は、染色体の異常の検出に関する。より詳細には、本発明は、生体試料から胎児の染色体の異常の確率を決定することに関する。
ダウン症候群は、比較的一般的な遺伝的疾患である。この症候群は、過剰な染色体21(トリソミー21すなわちT21)の存在、または、それほど頻繁ではないが、その染色体の過剰な実質的部分の存在によって引き起こされる。T13またはT18などの他のトリソミーも知られている。母体血漿からのDNA分子のDNA塩基配列決定法に基づいた、そのような状態の出生前診断の方法が知られている。生体試料から異常染色体を確実に判定することは、DNAの処理に関連付けられた変数を考えると、問題がある。
本発明の実施形態の目的は、従来技術の問題のうちの1つまたは複数を少なくとも軽減することである。
本発明の実施形態は、胎児の染色体の異常の確率を決定することに関する。本発明の実施形態では、標的染色体に対する胎児の染色体の異常を示す尤度比が決定される。尤度比は、後で説明するように、第1のモデルおよび第2のモデルに基づいて決定される。そのような異常を検出する際、誤った結果が判定されないことを可能な限り保証することが重要である。特に、偽陰性結果が決定される確率を低下させることが特に望ましい。しかしながら、データが効率的に使用されることを保証することも重要である。それは、結果が、許容可能な数のケースまたは検査において生成されるということである。検査が、その検査に関連付けられた1つまたは複数のパラメータにより結果が信頼できないことを示すのではなく、検査結果は、理想的には、可能ならば宣言されるべきである。たとえば、検査パラメータに関連付けられた所定のカットオフ値の使用によって、検査結果が生成されない、または検査結果が信頼できないと示され、それによって、検査を繰り返す、またはさらなる生体試料が必要とされる。
次に、本発明の実施形態について、添付の図を参照しながら、単に例として説明する。
本発明の一実施形態による装置の説明図である。 本発明の一実施形態による方法を示す流れ図である。 本発明の一実施形態によるモデルを形成するデータの図示である。 本発明の一実施形態によるモデルを形成するデータの図示である。 本発明の一実施形態による特定の染色体配列密度におけるモデルに対する密度曲線を示すグラフである。 本発明の一実施形態による異なる特定の染色体配列密度におけるモデルに対する密度曲線を示すグラフである。 標的染色体に関する第1のモデルおよび第2のモデルのための密度等高線(density contour)を示すグラフである。
図1は、本発明の一実施形態による装置100を示す。装置100は、メモリ・ユニット120に通信可能に結合された処理ユニット110を備える。処理ユニット110は、コンピュータ実行可能コードまたはソフトウェアを形成する命令を実行するように動作可能である。処理ユニット110は、1つまたは複数のプロセッサを備えてよい。処理ユニット110は、後で説明するように、コンピュータ実行可能コードに従って、生体試料から決定されたデータに対する動作を実行するように動作可能である。コンピュータ実行可能コードは、図2を参照して説明される本発明の一実施形態による方法を実行するように構成されてよい。メモリ・ユニット120は、ROMデバイスおよびRAMデバイスのうちの1つまたは複数などの1つまたは複数のメモリ・デバイスによって形成されてよい。メモリ・ユニット120は、記憶されたデータを読み出すために、およびデータをメモリ・ユニット120に記憶するために、処理ユニット110によってアクセス可能である。コンピュータ実行可能コードは、メモリ・ユニット120に記憶されてよい。メモリ・ユニット120は、後で説明するように、第1のモデルおよび第2のモデル121、122を表すデータを記憶してよい。
装置100は、データを受け取るためのインタフェース・ユニット130をさらに備える。特に、データは、DNA処理システム150から受け取られる。データ通信経路151が、インタフェース・ユニット130とDNA処理システム150の間に存在してよい。通信経路151は、インターネットを含む1つまたは複数の通信ネットワークを備えてよい。装置100は、後で説明するように、結果を示すデータを出力するための出力ユニット140をさらに備えてよい。
図2は、本発明の一実施形態による方法200を示す。方法200は、胎児の染色体の異常の確率を決定する方法である。
方法200は、女性被験者から取得された生体試料から標的染色体に関する第1のパラメータおよび染色体配列密度を表す第2のパラメータを示すデータを決定するステップ201を含む。このデータは、データをDNA処理システム150から受け取ることによって決定されてよい。データは、通信経路151を介してインタフェース・ユニット130によって受け取られてよい。
例として、染色体の異常を検出する方法を開示し、あらゆる目的のために参照により本明細書に組み込まれるWO2014/033455に対する参照がなされる。WO‘455に対する参照がなされるが、本発明の実施形態はこの点に限定されないことが理解されるであろう。
WO‘455では、配列データが、生体試料内の核酸分子に関して取得され、マッチング解析が、配列データ内の各核酸配列と参照ゲノムの一意の部分に対応する配列との間で実行され、したがって、各マッチングされた核酸は、参照ゲノム内の特定の染色体または前記染色体の一部に割り当てられる。1つまたは複数の参照染色体の各々に割り当てられたマッチングされた核酸の総数に対する標的染色体に割り当てられたマッチングされた核酸の総数の比が測定される。正常妊娠における比に対する、測定された比における統計的に有意な差は、胎児の標的染色体の異常を示す。標的染色体はT21に相当し得るが、標的染色体は他の染色体にも相当してよい。
標的染色体に関する第1のパラメータは、いくつかの実施形態では、染色体21、染色体18、または染色体13のうちの1つまたは複数などの1つまたは複数の標的染色体に割り当てられたマッチングされた核酸の総数を示し得るが、他の染色体が使用されてよいことが諒解されるであろう。
いくつかの実施形態では、第1のパラメータは、染色体カウント比、すなわち、1つまたは複数の参照染色体の各々に割り当てられたマッチングされた核酸の総数に対する各標的染色体に割り当てられたマッチングされた核酸の総数の比に対応してよい。しかしながら、いくつかの実施形態では、この比は、後で説明するように、第2のパラメータに基づいて装置100によって計算されてよい。
生体試料からの染色体配列密度を表す第2のパラメータは、いくつかの実施形態では、1つまたは複数の参照染色体の各々に割り当てられたマッチングされた核酸の総数を示してよい。いくつかの実施形態では、第2のパラメータは、常染色体片カウント数であってよい。常染色体片カウント数は、DNA処理システム150によってカウントされた常染色体片の数である。他の実施形態では、第2のパラメータは、参照染色体のシーケンシング・カバレッジの尺度であってよい。一実施形態では、第2のパラメータは、カバレッジの平均深度を示してよい。
方法200は、胎児の染色体の異常を示す尤度比を決定するステップ220を含む。尤度比は、それぞれのモデル121、122による染色体の異常の確率と染色体の正常の確率との比として決定されてよい。それぞれのモデルは、メモリ・ユニット120に記憶されてよい。
前述のように、装置100のメモリ・ユニット120は、第1のモデルおよび第2のモデル121、122を表すデータを記憶する。第1のモデル121は、染色体異常の確率を示す。第2のモデル122は、染色体正常を示す。第1のモデルおよび第2のモデル121、122は各々、第1のパラメータおよび第2のパラメータに基づいてパラメータ化される。
図3は、染色体異常の確率を示す第1のモデル121の説明図300である。図4は、染色体正常を示す第2のモデル122の説明図400である。説明図300、400は、トリソミー21に関連するモデルであるが、これは制限的でないことを理解されたい。各モデル121、122は、x軸に沿って示される染色体カウント比と、y軸に沿って示される、常染色体片カウント数すなわちアライメントされた片の総数であってよい染色体配列密度によってパラメータ化される。第1のモデルおよび第2のモデル121、122は、確率密度関数を表す。確率密度関数は、複合モデルを使用して決定されてよい。
第1のモデルおよび第2のモデル121、122の各々は、最小パラメータ値から最大パラメータ値の間で、すなわち、最小染色体カウント比値から最大染色体カウント比値の間、および最小常染色体片カウント数から最大常染色体片カウント数の間で、パラメータに対するそれぞれの確率を決定し得る。第1のモデルおよび第2のモデル121、122は、連続確率密度関数または不連続確率密度関数であってよい、すなわち、確率密度は、モデル121、122によって、最小パラメータ値から最大パラメータ値の間の範囲内の複数のパラメータ値に指定されてよい。パラメータ値が、指定された値の間にある場合、密度値は、特定のパラメータ値に対して補間されてよい。補間は、隣接する指定されたパラメータ値間の線形補間であってよい。
したがって、異常の確率染色体を示す第1のモデル121の場合、密度値は、
Non(Rin,A)
と決定され得、ここで、Rinは入力された染色体カウント比、Aは入力された常染色体片カウント数である。同様に、染色体正常を示す第2のモデル122の場合、密度値は、
Tri(Rin,A)
と決定され得る。
いくつかの実施形態では、モデル121、122の一方または両方は、
それぞれの染色体またはトリソミーに関する平均無影響(unaffected)染色体カウント比を示す値を含んでよい。
の値は、
1.トリソミー13:0.037036
2.トリソミー18:0.030858
3.トリソミー21:0.012843
であってよい。
上記の値は単に例に過ぎず、他の値も使用されてよいことが諒解されるであろう。
いくつかの実施形態では、モデル121、122の一方または両方は、第1の閾値
thresh1を含んでよい。第1の閾値は、それぞれの染色体またはトリソミーに関連付けられてよい。第1の閾値Rthresh1は、それより下では尤度比がゼロにきわめて近いと考えられ、非常に小さい値に設定されるべき染色体カウント比に対応してよい。
thresh1の値は、いくつかの実施形態では、
1.トリソミー13:0.037036
2.トリソミー18:0.030858
3.トリソミー21:0.012843
のうちの1つまたは複数であってよい。
上記の値は単に例に過ぎず、他の値も使用されてよいことが諒解されるであろう。
いくつかの実施形態では、モデル121、122の一方または両方は、第2の閾値
thresh2を含んでよい。第2の閾値は、それぞれの染色体またはトリソミーに関連付けられてよい。第2の閾値Rthresh2は、それより上では尤度比が事実上無限と考えられ、非常に大きな値に設定されるべき影響された染色体カウント比を示す。第2の閾値Rthresh2の値は、いくつかの実施形態では、
4.トリソミー13:0.039153
5.トリソミー18:0.032622
6.トリソミー21:0.013577
のうちの1つまたは複数であってよい。
上記の値は単に例に過ぎず、他の値も使用されてよいことが諒解されるであろう。
いくつかの実施形態では、モデル121、122の一方または両方は、任意選択の尤度比較正係数kLRを含んでよい。尤度比較正係数は、1というデフォルト値を有してよい。
図4から特に諒解され得るように、染色体正常に関する確率密度は、常染色体カウント数を減少させるにつれて、著しく減少する。この減少は、約300万の常染色体カウント数よりも下で特に観測可能である。すなわち、常染色体カウントの数が特に3e未満に減少するにつれて、同じ染色体カウント比値の場合でも、第2のモデル122によって決定される確率密度値は減少する。類似の影響は、第1のモデル121に対しても発生する。確率密度が減少するにつれて、各曲線の分散が増加してよく、すなわち、曲線は、より広い裾(tail)を有し、したがって、各曲線下面積は一定のままであってよい。
特定の入力パラメータ、すなわち、標的染色体に関する特定のカウント比および特定の常染色体片カウント数に対して、確率密度値は、第1のモデルおよび第2のモデル121、122の各々を使用して決定されてよい。たとえば、図5および図6は、第1の常染色体片カウント数および第2の常染色体片カウント数における、第1のモデルおよび第2のモデル121、122の各々に対する確率密度曲線を示す。図5は、2372258カウントの常染色体片カウント数の場合であり、図6は、5588457カウントの常染色体片カウント数の場合である。第2のモデルすなわち正常モデル122のための密度関数510、610と第1のモデルすなわち異常モデル121のための密度関数520、620との形状の差が観測可能である。
図5および図6から諒解され得るように、常染色体片カウント数における標的染色体のためのカウント比に基づいて、染色体の正常の確率は、正常確率分布510、610から決定可能であり、染色体の異常の確率は、異常確率分布520、620から決定可能である。
0.012715の染色体カウント比における、すなわち染色体片の約1.27%が染色体21に由来する、ライン530によって示される図5に示すように、染色体の正常の確率密度は4922.1であるが、染色体の異常の確率密度は0.00054354である。
尤度比(LR)は、
と決定され得る。
したがって、図5の場合、尤度比(LR)は1:9.05551eである。
0.013218の染色体カウント比における、すなわち染色体片の約1.32%が染色体21に由来する、ライン630によって示される図6を参照すると、染色体の正常の確率密度は1.2288e−5であるが、染色体の異常の確率密度は6.68444eである。したがって、LRは、6.68444e:1と決定される。
いくつかの実施形態では、Rin>Rthresh2(ここで、上記で説明したように、Rthresh2は第2の閾値である)であり、かつDTri=0である場合、LRは、比較的に高い値などの所定の値を有するように決定され得る。LRの所定の値は、メモリ・ユニット120によって記憶され得る最大値などの、無限を表す値であってよいが、他の所定の値が使用されてもよいことが理解されよう。
いくつかの実施形態では、Rin<Rthresh1(ここで、上記で説明したように、Rthresh1は第1の閾値である)かつDNon=0である場合、LRは、比較的低い値などの所定の値を有するように決定され得る。LRの所定の値は、ゼロの値であってよい。LRの所定の値は、メモリ・ユニット120によって記憶され得る最小値であってよいが、他の所定の値が使用されてもよいことが理解されよう。値は、最小の非ゼロ値であってもよい。
上記の第1の閾値および/または第2の閾値の状況のどちらも発生せず、ゼロの値にされたDNonによりLRが不定の場合、いくつかの実施形態では、LRは、メモリ・ユニット120に記憶され得る、機械が表すことができる最大値であるように決定されることがある。
図7は、第1のモデルおよび第2のモデル121、122に基づいた密度等高線を示す。各モデル121、122に対して、複数の特定の密度を表す密度等高線が示されている。図7に示されている特定の密度等高線は、各モデル121、122に対する総密度の50、80、90、および99.5%を囲む。染色体異常の確率を示す第1のモデル121および染色体正常の確率を示す第2のモデル122に対する密度等高線710が示されているが、わかりやすくするために、図7では、すべてにラベルが付されているとは限らない。図7では、重複領域が、点線の境界線730とともに示されている。この重複領域は、第1のモデルおよび第2のモデル121、122が重複する領域を示す。重複は、第1の所定の尤度比から第2の所定の尤度比の間にあり、図7では、第1の所定の尤度比および第2の所定の尤度比は1:100および100:1であるが、重複領域は、他の所定の尤度比に対しても決定されてよい。同様に、1の尤度比(1:1)の線が、参照740とともに示されている。諒解され得るように、重複領域は、湾曲した境界線を有し、1のLRの線740も湾曲している。特に、重複領域は、常染色体カウント数の減少とともに広がることが諒解され得る。したがって、図7は、決定された密度およびLR上で常染色体カウント数および染色体カウント比を変化させる影響を示す。図7は染色体21に関するが、これは限定的でないことが理解されるであろう。
図2に戻ると、ステップ230では、いくつかの誤った結果を減少させるための性能パラメータ閾値が決定される。一実施形態では、決定される性能パラメータは胎児由来断片閾値(FFT)である。いくつかの実施形態では、FFTはFFN_minと呼ばれる。FFTは、偽陰性結果すなわち検査が染色体の異常に対して陰性であると間違って決定される場合の数を減少させるために使用される。FFTは、本発明の実施形態では、動的に決定される。動的決定は、過剰な数の検査を信頼できないと宣言させ得る、所定の胎児由来断片閾値すなわち固定された胎児由来断片閾値の使用を回避することができる。FFTの動的決定は、常染色体片カウントなどの染色体配列密度を表す第2のパラメータまたはシーケンシング・カバレッジを示す他のパラメータに基づく。
性能パラメータの推定値は、DNA処理システムによって提供される。特に、いくつかの実施形態では、胎児由来断片の推定値
は、DNA処理システム150によって決定される。胎児由来断片推定値は、胎児から発生した常染色体片の断片の推定値である。本発明の実施形態は、胎児の染色体の異常を決定することを目的とするので、十分な常染色体片が母親ではなく胎児から発生することが必要である。胎児由来断片の推定値は、この推定値に関連付けられた不確定値の指示を伴うことがある。
動的胎児由来断片推定をサポートするために必要とされる最小カウント比RFN_minはカウント比の値として決定され、この値より下では、第1のモデル121からの比例する正常確率密度関数下面積DNon(Rin,A)は、所定の最低感度すなわちMinSensitivityであり、上記のように、Aは総常染色体片カウントである。所定の最低感度すなわちMinSensitivityは、各それぞれの染色体またはトリソミーに対して決定されてよい。MinSensitivityは、
・トリソミー13:0.80
・トリソミー18:0.98
・トリソミー21:0.99
であってよいが、他の値が選択されてもよいことが理解されよう。最小カウント比
FN_minは、次のように、最小胎児由来断片値に変換されてよい。
ここで、上記のように、
は、それぞれの染色体またはトリソミーに関する平均無影響染色体カウント比を示す。
offsは、胎児由来断片妥当性検査閾値オフセットであり、
offs=aFT+bFTσ
と定義され得、ここで、σは胎児由来断片推定段階によって生成される胎児由来断片推定分布の標準偏差、aFTおよびbFTは妥当性検査閾値決定パラメータである。胎児由来断片推定分布の標準偏差は、経験的誤差分布に適合するなどのために、あらかじめ決定されている可能性がある。これらの決定パラメータは、aFT=0および
の一方または両方の初期値を有してよいが、他の初期値が使用されてもよいことが理解されよう。
ステップ230では、試料に対する胎児由来断片の推定値
が動的FFTFFN_minよりも小さいかどうかが判定される。いくつかの実施形態では、胎児由来断片推定値は、胎児の数に応じてスケーリングされる。ステップ230は、
であるかどうかを判定することを含んでよく、ここで、τは多胎妊娠補正である。生体試料が多胎妊娠からの試料である場合、τは、妊娠において身ごもる胎児の数に等しくなるべきであり、そうでない場合、τは1であるべきである。
ステップ240において、
または
閾値よりも小さいと決定される場合、検査は信頼できない可能性があると判定され、方法200はステップ250に進む。ステップ250では、生体試料中の胎児由来断片が少ないことを示す適切な警告が出力されてよい。そうでない場合、ステップ260では、検査の出力が出力されてよい。出力は、試料に対して決定されたLRであってよい。
有利には、本発明の実施形態は、片カウントまたは胎児由来断片検査が検査結果を否認する正確な検査結果が生成されることを可能にすることによって、たとえばプレフィルタリング手法よりも効率的なデータの使用をもたらす。たとえば、本発明の実施形態は、片カウントが少ないまたは胎児由来断片が少ないにもかかわらず、片カウント割合は十分に決定的な「陽性」結果を表す結果を生じさせる。
本発明の実施形態は、ハードウェア、ソフトウェア、またはハードウェアとソフトウェアの組み合わせの形で実現可能であることが諒解されるであろう。任意のそのようなソフトウェアは、消去可能または書き換え可能であろうとなかろうと、たとえばROMのような記憶デバイスなどの、揮発性記憶または不揮発性記憶の形で記憶されてもよいし、たとえば、RAM、メモリ・チップ、デバイス、もしくは集積回路などのメモリの形で記憶されてもよいし、たとえば、CD、DVD、磁気ディスク、または磁気テープなどの、光学的または磁気的に可読な媒体上に記憶されてもよい。記憶デバイスおよび記憶媒体は、実行されると本発明の実施形態を実施する1つまたは複数のプログラムを記憶するのに適した機械可読媒体の実施形態であることが諒解されるであろう。したがって、実施形態は、任意の先行請求項において特許請求されるシステムまたは方法を実施するためのコードを備えるプログラムと、そのようなプログラムを記憶する機械可読記憶とを提供する。またさらに、本発明の実施形態は、ワイヤード接続またはワイヤレス接続上で搬送される通信信号などの任意の媒体を介して電子的に搬送されてよく、実施形態は、それらを適切に包含する。
本明細書(あらゆる添付の特許請求の範囲、要約書、および図面を含む)に開示されている特徴のすべて、および/またはそのように開示されるあらゆる方法もしくはプロセスのステップのすべては、そのような特徴および/またはステップのうちの少なくともいくつかが相互に排他的である組み合わせを除いて、任意の組み合わせで組み合わされてよい。
本明細書(あらゆる添付の特許請求の範囲、要約書、および図面を含む)に開示されている各特徴は、別段に明記されていない限り、同じ目的、等価な目的、または類似の目的にかなう代替特徴によって置き換えられてよい。したがって、別段に明記されていない限り、開示されている各特徴は、一般的な一連の等価な特徴または類似の特徴の一例に過ぎない。
本発明は、任意の前述の実施形態の詳細に制限されない。本発明は、本明細書(あらゆる添付の特許請求の範囲、要約書、および図面を含む)に開示されている特徴の任意の新規な特徴もしくは任意の新規な組み合わせ、またはそのように開示される任意の方法もしくはプロセスのステップの任意の新規なステップもしくは任意の新規な組み合わせに拡張される。特許請求の範囲は、単に前述の実施形態だけでなく、特許請求の範囲に含まれる任意の実施形態も包含するように解釈されるべきである。

Claims (22)

  1. 胎児の染色体の異常の確率を判定するコンピュータ実装方法であって、
    女性被験者から取得された生体試料から標的染色体に関する第1のパラメータおよび染色体配列密度を表す第2のパラメータを示すデータを決定することと、
    胎児の染色体の異常を示す尤度比を決定することであって、該尤度比が、該第1のパラメータおよび該第2のパラメータに基づいてそれぞれの異常モデルおよび正常モデルに従って染色体の異常の確率と染色体の正常の確率の比として決定することと、
    1つまたは複数の性能パラメータ閾値を決定することと、
    該試料に関連付けられた1つまたは複数の性能パラメータの推定値を該1つまたは複数の性能パラメータ閾値と比較することと
    を含む方法。
  2. 前記性能パラメータ閾値のうちの少なくとも1つが、前記異常モデルおよび前記正常モデルの一方または両方から決定される、請求項1記載の方法。
  3. 前記1つまたは複数の性能パラメータ閾値が胎児由来断片閾値を含む、請求項1記載の方法。
  4. 前記胎児由来断片閾値が前記正常モデルに基づいて決定される、請求項3記載の方法。
  5. 前記比較することが、前記試料に対する胎児由来断片の推定値を前記胎児由来断片閾値と比較することを含む、請求項3または4記載の方法。
  6. 前記性能パラメータの前記推定値が前記性能パラメータ閾値よりも小さいとき、前記生体試料を、低い性能パラメータを有すると判定することを含む、請求項1から5のいずれか記載の方法。
  7. 前記胎児由来断片の推定値が前記胎児由来断片閾値よりも小さいとき、前記生体試料を、低い胎児由来断片を有すると判定することを含む、請求項6記載の方法。
  8. 前記性能パラメータ閾値が、前記第1のパラメータおよび前記第2のパラメータに基づいて動的に決定される、請求項1から7のいずれか記載の方法。
  9. 前記胎児由来断片閾値が、前記第1のパラメータおよび前記第2のパラメータに基づいて動的に決定される、請求項3またはそれに従属するいずれかの請求項に従属する場合の、請求項8記載の方法。
  10. 前記胎児由来断片閾値が最小カウント比に基づいて決定される、請求項3またはそれに従属するいずれかの請求項記載の方法。
  11. 前記最小カウント比が、前記正常モデルからの確率密度関数下面積が所定の最低感度であることに基づく、請求項10記載の方法。
  12. 前記胎児由来断片閾値RFN_minが、
    として決定され、ここで、RFN_minは最小カウント比であり、
    は、前記標的染色体に関する平均無影響染色体カウント比である、請求項3またはそれに従属するいずれかの請求項記載の方法。
  13. 前記尤度比(LR)が、
    として決定され、ここで、DTri(Rin,A)は前記異常モデルに従って決定される染色体の異常の確率であり、DNon(Rin,A)は前記正常モデルに従って決定される染色体の正常の確率であり、Rinは前記標的染色体に関する前記第1のパラメータであり、Aは染色体配列密度を表す前記第2のパラメータである、請求項1から12のいずれか記載の方法。
  14. 前記染色体カウント比が、1つまたは複数の参照染色体の各々に割り当てられたマッチングされた核酸の総数に対する標的染色体に割り当てられたマッチングされた核酸の総数の比である、請求項1から13のいずれか記載の方法。
  15. 前記第1のパラメータが染色体カウント比を示す、請求項1から14のいずれか記載の方法。
  16. 前記第2のパラメータが常染色体カウント数を示す、請求項1から15のいずれか記載の方法。
  17. 請求項1から16のいずれか記載の方法を実行するように構成された装置。
  18. 処理ユニットと、
    異常モデルおよび正常モデルを表すデータと、該処理ユニットによって実行されると、
    女性被験者から取得された生体試料から標的染色体に関する第1のパラメータおよび染色体配列密度を表す第2のパラメータを示すデータを決定するステップと、
    胎児の染色体の異常を示す尤度比を決定するステップであって、該尤度比が、該第1のパラメータおよび該第2のパラメータに基づいて該それぞれの異常モデルおよび正常モデルに従って染色体の異常の確率と染色体の正常の確率の比として決定するステップと、
    1つまたは複数の性能パラメータ閾値を決定するステップと、
    該試料に関連付けられた1つまたは複数の性能パラメータの推定値を該1つまたは複数の性能パラメータ閾値と比較するステップと
    を実行するコンピュータ実行可能命令とを記憶するメモリ・ユニットと
    を備える装置。
  19. DNA処理システムからデータを受け取るためのインタフェース・ユニットを備える、請求項18記載の装置。
  20. コンピュータによって実行されると、請求項1から16のいずれか記載の方法を実行するように構成されるコンピュータ・ソフトウェア。
  21. コンピュータ可読媒体上に記憶された、請求項20記載のコンピュータ・ソフトウェア。
  22. コンピュータ可読媒体上に有形に記憶された、請求項20記載のコンピュータ・ソフトウェア。
JP2018533694A 2015-12-22 2016-12-21 染色体異常の検出 Active JP7079433B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
GBGB1522665.7A GB201522665D0 (en) 2015-12-22 2015-12-22 Detection of chromosome abnormalities
GB1522665.7 2015-12-22
PCT/GB2016/054019 WO2017109487A1 (en) 2015-12-22 2016-12-21 Detection of chromosome abnormalities

Publications (3)

Publication Number Publication Date
JP2019508781A true JP2019508781A (ja) 2019-03-28
JP2019508781A5 JP2019508781A5 (ja) 2020-02-06
JP7079433B2 JP7079433B2 (ja) 2022-06-02

Family

ID=55311454

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018533694A Active JP7079433B2 (ja) 2015-12-22 2016-12-21 染色体異常の検出

Country Status (6)

Country Link
US (1) US20190002966A1 (ja)
EP (1) EP3394776A1 (ja)
JP (1) JP7079433B2 (ja)
CN (1) CN109074427B (ja)
GB (1) GB201522665D0 (ja)
WO (1) WO2017109487A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115001997B (zh) * 2022-04-11 2024-02-09 北京邮电大学 基于极值理论的智慧城市网络设备性能异常阈值评估方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070148631A1 (en) * 2003-07-25 2007-06-28 University Of Plymouth Methods and apparatus for screening for chromosomal abmormalities
US20100112590A1 (en) * 2007-07-23 2010-05-06 The Chinese University Of Hong Kong Diagnosing Fetal Chromosomal Aneuploidy Using Genomic Sequencing With Enrichment
US20120190557A1 (en) * 2011-01-25 2012-07-26 Aria Diagnostics, Inc. Risk calculation for evaluation of fetal aneuploidy
WO2013050383A2 (en) * 2011-10-04 2013-04-11 Cezanne S.A.S In vitro method to determine the risk of a pregnant woman of carrying a fetus with fetal aneuploidy
US20150267255A1 (en) * 2012-08-30 2015-09-24 Premaitha Health Ltd. Method of detecting chromosomal abnormalities

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EA035451B9 (ru) * 2007-07-23 2020-09-09 Те Чайниз Юниверсити Ов Гонконг Способ диагностики рака с использованием геномного секвенирования
DK2728014T3 (en) * 2012-10-31 2016-01-25 Genesupport Sa A non-invasive method for the detection of fetal chromosomal aneuploidy
WO2015035555A1 (zh) * 2013-09-10 2015-03-19 深圳华大基因科技有限公司 用于确定胎儿是否存在性染色体数目异常的方法、系统和计算机可读介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070148631A1 (en) * 2003-07-25 2007-06-28 University Of Plymouth Methods and apparatus for screening for chromosomal abmormalities
US20100112590A1 (en) * 2007-07-23 2010-05-06 The Chinese University Of Hong Kong Diagnosing Fetal Chromosomal Aneuploidy Using Genomic Sequencing With Enrichment
US20120190557A1 (en) * 2011-01-25 2012-07-26 Aria Diagnostics, Inc. Risk calculation for evaluation of fetal aneuploidy
WO2013050383A2 (en) * 2011-10-04 2013-04-11 Cezanne S.A.S In vitro method to determine the risk of a pregnant woman of carrying a fetus with fetal aneuploidy
US20150267255A1 (en) * 2012-08-30 2015-09-24 Premaitha Health Ltd. Method of detecting chromosomal abnormalities

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"Optimal Detection of Fetal Chromosomal Abnormalities byMassively Parallel DNA Sequencing of Cell-Fre", CLINICAL CHEMISTRY, vol. 57, no. 7, JPN7021000646, 2011, pages 1042 - 1049, ISSN: 0004457646 *

Also Published As

Publication number Publication date
US20190002966A1 (en) 2019-01-03
EP3394776A1 (en) 2018-10-31
JP7079433B2 (ja) 2022-06-02
WO2017109487A1 (en) 2017-06-29
GB201522665D0 (en) 2016-02-03
CN109074427B (zh) 2023-03-31
CN109074427A (zh) 2018-12-21

Similar Documents

Publication Publication Date Title
CN109239360B (zh) 一种反应曲线异常检测方法及装置
EP3576093A1 (en) Method and system for clinical effectiveness evaluation of artificial intelligence based medical device
US9607128B2 (en) Detection and correction of jumps in real-time PCR signals
CN112365927B (zh) Cnv检测装置
US20130173215A1 (en) Adaptive trend-change detection and function fitting system and method
JP2018514234A (ja) 染色体異数性を測定するためのキット、装置及び方法
CN116386718B (zh) 检测拷贝数变异的方法、设备和介质
CN109101661A (zh) 一种数据样本集合中异常点的检测方法和装置
CN111584002B (zh) 用于检测肿瘤突变负荷的方法、计算设备和计算机存储介质
KR101678962B1 (ko) 대규모 병렬형 게놈서열분석 방법을 이용한 비침습적 산전검사 장치 및 방법
JP5438466B2 (ja) 画像処理装置および画像処理方法
JP2019508781A (ja) 染色体異常の検出
EP2977466B1 (en) Detecting chromosomal aneuploidy
CN107075567B (zh) 用于荧光数据校正的方法
CN114708908B (zh) 检测实体瘤微小残留病灶的方法、计算设备和存储介质
CN113990492B (zh) 确定关于实体瘤微小残留病灶的检测参数的方法、设备和存储介质
CN115394359A (zh) 一种鉴定人类胚胎细胞染色体变异的方法及应用
CN107109324A (zh) 确定胎儿核酸含量的方法和装置
US20170323051A1 (en) Devices, systems, and methods for high-resolution melt analysis
KR102349023B1 (ko) 뉴클레오티드 서열 변이의 빈도 분포 결정
JP2016184185A5 (ja)
KR20150082976A (ko) 반도체 생산 공정에서 센서 데이터들을 이용하여 웨이퍼의 수율을 분석하는 방법
CN113851189A (zh) 一种自适应检测pcr荧光基线的方法、装置及其设备
US20160265051A1 (en) Methods for Detection of Fetal Chromosomal Abnormality Using High Throughput Sequencing
CN108733984B (zh) Nipt的z值结果校正方法及装置、计算机可读存储介质

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20191217

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191217

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210309

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20210531

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210804

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20211214

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211216

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220412

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220511

R150 Certificate of patent or registration of utility model

Ref document number: 7079433

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150