JP6522059B2 - 診断のための赤血球動力学 - Google Patents

診断のための赤血球動力学 Download PDF

Info

Publication number
JP6522059B2
JP6522059B2 JP2017144114A JP2017144114A JP6522059B2 JP 6522059 B2 JP6522059 B2 JP 6522059B2 JP 2017144114 A JP2017144114 A JP 2017144114A JP 2017144114 A JP2017144114 A JP 2017144114A JP 6522059 B2 JP6522059 B2 JP 6522059B2
Authority
JP
Japan
Prior art keywords
volume
red blood
hemoglobin
hemoglobin content
ida
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.)
Active
Application number
JP2017144114A
Other languages
English (en)
Other versions
JP2018009992A (ja
Inventor
エム. ヒギンズ,ジョン
エム. ヒギンズ,ジョン
Original Assignee
ザ ジェネラル ホスピタル コーポレイション
ザ ジェネラル ホスピタル コーポレイション
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 ザ ジェネラル ホスピタル コーポレイション, ザ ジェネラル ホスピタル コーポレイション filed Critical ザ ジェネラル ホスピタル コーポレイション
Publication of JP2018009992A publication Critical patent/JP2018009992A/ja
Application granted granted Critical
Publication of JP6522059B2 publication Critical patent/JP6522059B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/02Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61PSPECIFIC THERAPEUTIC ACTIVITY OF CHEMICAL COMPOUNDS OR MEDICINAL PREPARATIONS
    • A61P1/00Drugs for disorders of the alimentary tract or the digestive system
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61PSPECIFIC THERAPEUTIC ACTIVITY OF CHEMICAL COMPOUNDS OR MEDICINAL PREPARATIONS
    • A61P3/00Drugs for disorders of the metabolism
    • A61P3/12Drugs for disorders of the metabolism for electrolyte homeostasis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61PSPECIFIC THERAPEUTIC ACTIVITY OF CHEMICAL COMPOUNDS OR MEDICINAL PREPARATIONS
    • A61P7/00Drugs for disorders of the blood or the extracellular fluid
    • A61P7/06Antianaemics
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/50Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
    • G01N33/68Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids
    • G01N33/6893Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids related to diseases not provided for elsewhere
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/50Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
    • G01N33/72Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving blood pigments, e.g. haemoglobin, bilirubin or other porphyrins; involving occult blood
    • G01N33/721Haemoglobin
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/50Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
    • G01N33/80Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving blood groups or blood types or red blood cells
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/50Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
    • G01N33/94Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving narcotics or drugs or pharmaceuticals, neurotransmitters or associated receptors
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2800/00Detection or diagnosis of diseases
    • G01N2800/06Gastro-intestinal diseases
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2800/00Detection or diagnosis of diseases
    • G01N2800/22Haematology
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2800/00Detection or diagnosis of diseases
    • G01N2800/22Haematology
    • G01N2800/224Haemostasis or coagulation

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Hematology (AREA)
  • Molecular Biology (AREA)
  • Immunology (AREA)
  • Biomedical Technology (AREA)
  • Urology & Nephrology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medicinal Chemistry (AREA)
  • Biochemistry (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Microbiology (AREA)
  • Biotechnology (AREA)
  • Food Science & Technology (AREA)
  • Cell Biology (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Organic Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Pharmacology & Pharmacy (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • General Chemical & Material Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Wood Science & Technology (AREA)
  • Zoology (AREA)
  • Diabetes (AREA)
  • Biophysics (AREA)
  • General Engineering & Computer Science (AREA)
  • Genetics & Genomics (AREA)
  • Obesity (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)

Description

優先権主張
本出願は、参照によりその全体が本明細書に組み込まれている、2010年9月16日
に出願された米国仮出願第61/383,357号の優先権を主張する。
連邦政府の支援による研究または開発
本発明は、国立衛生研究所により付与された認可番号第DK083242号および第H
L091331号の下で、政府の支援により行われた。政府は本発明における確実な権利
を有する。
本発明は、貧血または貧血予備軍の患者を特定し、サラセミア形質貧血を鉄欠乏性貧血
および他の原因の貧血と区別するための方法に関する。この方法を使用して、治療、例え
ば、鉄補給剤および/またはエリスロポエチン(EPO)等の赤血球生成促進剤による療
法を最適化することもできる。また、例えば前記方法に使用するための装置およびシステ
ムについても記載する。
循環するヒト赤血球(RBC)の数、大きさ、ヘモグロビン濃度、その他の特徴を制御
するシステムについては、ほとんど解明されていない。RBCは骨髄から放出された後、
未知の機構(Lew VL他、(1995)、Blood 86、334-341、Waugh RE他、(1992)、Blood 79、13
51-1358)により、体積および総ヘモグロビン含有量が減少し、約120日後に未知の要因
に反応して除去される。
統計物理学による理論および病院の臨床検査室からのデータ(d’Onofrio G他、(1995)
Blood 85、818-823)の一部に基づき、本発明は、RBC成熟およびクリアランスのマス
タ方程式モデルを使用する。モデルは、貧血患者を正確に特定し、サラセミア形質貧血を
鉄欠乏性貧血と区別する。驚いたことに、このモデルは、貧血が臨床的に検出可能になる
数週間前に、多くの貧血予備軍患者を特定する。
一態様では、被験者が鉄欠乏性貧血(IDA)を発症する危険を判定する方法は、被験
者の赤血球を含むサンプル中で、母平均赤血球ヘモグロビン濃度(MCHC)を判定する
ステップと、赤血球体積およびヘモグロビン含有量をMCHC線上に投影することにより
、各赤血球体積およびヘモグロビン含有量を指標に変換するステップとを含む。また方法
は、変換された体積およびヘモグロビン含有量指標がMCHC線上の平均投影位置の閾値
百分率よりも低いサンプル中の赤血球分画を判定して、サンプル分画を提供するステップ
と、サンプル分画を基準分画と比較するステップとを含む。基準分画よりも高いサンプル
分画があれば、被験者がIDAを発症する危険があることを示す。
別の態様では、消化器(GI)疾患をスクリーニングするために被験者をスクリーニン
グまたは選択する方法は、被験者の赤血球を含むサンプル中で、母平均赤血球ヘモグロビ
ン濃度(MCHC)を判定するステップと、赤血球体積およびヘモグロビン含有量をMC
HC線上に投影することにより、各赤血球体積およびヘモグロビン含有量を指標に変換す
るステップとを含む。また方法は、変換された体積およびヘモグロビン含有量指標がMC
HC線上の平均投影位置の閾値百分率よりも低いサンプル中の赤血球分画を判定するステ
ップと、サンプル分画を基準分画と比較するステップとを含む。方法は、GI評価のため
に被験者を選択するか、またはサンプル分画が基準分画よりも低い場合にさらにスクリー
ニングするステップをさらに含む。
別の態様では、本開示は、小球性貧血を有する被験者、例えば、慢性疾患のない被験者の鉄欠乏性貧血(IDA)とサラセミア形質(TT)とを鑑別診断する方法を特徴とする。方法は、被験者の赤血球のサンプル中で、赤血球間のヘモグロビン含有量減少率の変動の大きさの値(D を判定するステップと、Dを基準値と比較するステップとを含む。基準値よりも高いDがあれば、被験者がIDAであるかIDAになる可能性が高いことを示し、基準値よりも低いDがあれば、被験者がTTであるかTTになる可能性が高いことを示す。
別の態様では、健常者の血液ドーピングまたは赤血球生成促進剤の存在または使用を検出するための方法は、ヘモグロビン含有量減少率の変動の大きさの値(D 、および被験者のサンプル中の赤血球体積減少率の変動の大きさの値(D の一方または両方の値を判定するステップを含む。基準値よりも高いDおよび/またはDがあれば、血液ドーピングの存在もしくは使用または赤血球生成促進剤の使用を示す。
別の態様では、被験者のエリスロポエチン(EPO)または他の赤血球生成促進剤(E
SA)の投与量を最適化するための方法は、EPOまたはESA治療を受けている被験者
のサンプル中の正規化臨界体積(v)を判定するステップを含む。また、方法は、v
が下位基準レベルよりも低い場合にEPOの投与量を増加させ、vが上位基準レベルよ
りも高い場合にEPOの投与量を減少させるステップを含む。
別の態様では、鉄剤を使用して被験者を治療する方法が、被験者のサンプル中の正規化
臨界体積(v)を判定するステップと、vが下位基準レベルよりも低い場合に、1回
分の鉄剤を投与するか、または鉄補給剤を処方するステップとを含む。
別の態様では、本開示は、命令を記憶するためのメモリと、記憶された命令を実行可能
な1もしくは複数のプロセッサまたは処理装置とを有するコンピューティング装置を備え
たシステムを特徴とする。記憶された命令は、実行されると、サンプル中の各赤血球の体
積およびヘモグロビン濃度または含有量を表すデータを受信するステップと、平均赤血球
ヘモグロビン濃度(MCHC)を表す線上に赤血球体積およびヘモグロビン含有量を投影
することにより、各赤血球体積およびヘモグロビン含有量を指標に変換するステップとを
含む動作を行う。また、動作は、変換された体積およびヘモグロビン含有量指標が平均投
影の閾値百分率よりも低いサンプル中の赤血球分画を計算するステップと、分画を表す出
力を提供するステップとを含む。
別の態様では、本開示は、1または複数のプロセッサまたは処理装置により実行可能な
命令を記憶するように構成された、1または複数の機械読取可能な記憶装置を特徴とする
。命令は、1または複数のプロセッサにより実行されると、サンプル中の各赤血球の体積
およびヘモグロビン濃度を表すデータを受信するステップと、平均赤血球ヘモグロビン濃
度(MCHC)を表す線上に赤血球体積およびヘモグロビン含有量を投影することにより
、各赤血球体積およびヘモグロビン含有量を指標に変換するステップとを含む動作を行う
。また、動作は、変換された体積およびヘモグロビン含有量指標が平均投影の閾値百分率
よりも低いサンプル中の赤血球分画を計算するステップと、分画を表す出力を提供するス
テップとを含む。
別の態様では、本開示は、(i)母平均赤血球ヘモグロビン濃度(MCHC)、赤血球体積およびヘモグロビン含有量をMCHC線上に投影することによる、各赤血球体積およびヘモグロビン含有量の指標への変換、投影が平均投影の閾値百分率よりも低いサンプル中の赤血球分画、(ii)ヘモグロビン含有量減少率の変動の大きさの値(D 、(iii)赤血球体積減少率の変動の大きさの値(D 、ならびに(iv)正規化臨界体積またはクリアランス閾値(v)の1または複数を計算する旨の機械読取可能な命令を含む、1または複数の機械読取可能な記憶装置を特徴とする。
別の態様では、サラセミア形質(TT)について患者をスクリーニングする方法は、被験者のサンプル中のヘモグロビン含有量減少率の変動の大きさの値(D を判定するステップと、範囲の参照のためにDを比較するステップとを含む。Dが基準範囲外にあれば、患者がTTであるか、TTになる可能性があることを示し、かつ/あるいは、Dが基準範囲外にあれば、遺伝子型決定もしくはさらなる検査のために被験者を選択することを示し、または被験者が妊婦である場合には、その父親をスクリーニングのために選択することを示す。
別の態様では、本開示は、貧血を有する(例えば、貧血と診断された)被験者において、慢性疾患(ACD)に伴う貧血と他の原因の貧血とを鑑別診断するための方法を特徴とする。方法は、(i)ヘモグロビン含有量減少率の変動の大きさの値(D 、(ii)赤血球体積減少率の変動の大きさの値(D 、(iii)正規化臨界体積またはクリアランス閾値(v)、(iv)低速相体積およびヘモグロビン含有量の平均減少率(α)、(v)高速相体積の平均減少率(β)、ならびに(vi)高速相ヘモグロビン含有量の平均減少率(β)の1または複数の値を判定するステップと、これらの値を、ACDの診断確率が基準範囲内の値の数によって判定される基準範囲と比較するステップと、この確率を閾値と比較するステップと、確率が閾値よりも高い場合に、輸血療法をしない選択をし、かつ/または鉄レベルの追加検査を行わないことを選択するステップとを含む。
別の態様では、本開示は、コンピュータ読取可能な命令を記憶するように構成されたコ
ンピュータ読取可能な記憶装置を特徴とする。この命令は、1または複数のプロセッサに
より実行されると、網状赤血球または赤血球の時間依存性体積−ヘモグロビン分布に関連
する第1集合のパラメータを受信するステップと、第1集合のパラメータに基づいて赤血
球の定常状態分布を推定するステップとを含む動作を行わせる。また、動作は、赤血球の
推定定常状態分布と経験分布との相違の程度を表す目的関数を計算するステップを含む。
別の態様では、本開示は、コンピューティング装置を備えたシステムを特徴とする。コ
ンピューティング装置は、命令を記憶するためのメモリと、記憶された命令を実行して種
々の動作を行うことのできる1または複数のプロセッサまたは処理装置を有する検出モジ
ュールとを備える。動作は、網状赤血球または赤血球の時間依存性体積−ヘモグロビン分
布に関連する第1集合のパラメータを受信するステップと、第1集合のパラメータに基づ
いて赤血球の定常状態分布を推定するステップとを含む。また、動作は、赤血球の推定定
常状態分布と経験分布との相違の程度を表す目的関数を計算するステップを含む。
別の態様では、コンピュータにより実施される方法が、網状赤血球または赤血球の時間
依存性体積−ヘモグロビン分布に関連する第1集合のパラメータを受信するステップと、
第1集合のパラメータに基づいて赤血球の定常状態分布を推定するステップとを含む。ま
た、方法は、赤血球の推定定常状態分布と経験分布との相違の程度を表す目的関数を計算
するステップを含む。
前述した方法、システム、コンピュータ読取可能な記憶装置の実施は、以下の特徴の組
合せを含むことができる。
被験者の母平均体積
の75%よりも大きいv、または78〜82%のvを維持するように、EPOの投与
量を調節することができる。被験者の母平均体積
の78%よりも大きいv、または78〜82%のvを維持するように、鉄剤の投与量
を調節することができる。閾値百分率を、MCHC線上の母平均投影の70%、75%、
80%、85%、90%、または95%とすることができる。基準分画は、約0.10、
0.11、0.12、0.13、0.14、または0.15である。サンプルは、被験者
の全血を含むことができる。サンプル中の平均赤血球ヘモグロビン濃度(MCHC)を表
すデータを受信することができる。検出モジュールを、赤血球(RBC)の赤血球体積(
CV)、平均赤血球体積(MCV)、血球ヘモグロビン濃度(CHC)、平均赤血球ヘモ
グロビン濃度(MCHC)、および平均赤血球ヘモグロビン含有量(MCH)の1または
複数、例えば、すべてと、これらの母集団統計とを検出するように構成することができる
。体積およびヘモグロビン含有量座標が、すべての体積およびヘモグロビン含有量座標の
最小二乗線形適合を表す線上に投影されるように、個々の赤血球の体積およびヘモグロビ
ン含有量を変換することができる。
コンピュータまたは機械読取可能な記憶装置は、(i)ヘモグロビン含有量減少率(D
)の変動の大きさ、(ii)赤血球体積減少率(D)の変動の大きさ、(iii)正
規化臨界体積またはクリアランス閾値(v)、(iv)低速相体積およびヘモグロビン
含有量の平均減少率(α)、(v)高速相体積の平均減少率(β)、ならびに(vi)
高速相ヘモグロビン含有量の平均減少率(β)の1または複数を計算する旨の命令を含
むことができる。
目的関数の値が所定の閾値条件を満たす場合に、第1集合のパラメータを調節して第2
集合のパラメータを提供することができる。第2集合のパラメータは、対応する推定定常
状態分布と経験分布との相違の程度が低下するようになっている。第2集合のパラメータ
を、第1集合のパラメータとして1または複数のプロセッサに提供することができる。網
状赤血球または赤血球の時間依存性体積−ヘモグロビン分布は、患者から得た血液サンプ
ルに基づくことができる。定常状態分布を、線形演算子と状態の測定した分布との組合せ
に基づいて推定することができる。線形演算子は、少なくとも1つのヤコビアンおよび少
なくとも1つのラプラシアンを含むことができる。
特に定義のない限り、本明細書で使用されるすべての技術的および科学的用語は、本発
明が属する技術の当業者が通常理解するものと同一の意味を有する。本発明で使用する方
法および材料について本明細書で説明するが、当技術で公知の他の適切な方法および材料
を使用してもよい。材料、方法、および例は例示的なものにすぎず、限定的なものではな
い。本明細書で言及するすべての刊行物、特許出願、特許、シーケンス、データベースエ
ントリ、および他の参考文献は、参照により組み込まれている。矛盾がある場合には、定
義を含めて本明細書に従う。
以下の詳細な説明、図面、および特許請求の範囲から、本発明の他の特徴および利点が
明らかになろう。
図1A及び1Bは、末梢循環における平均RBCの体積およびヘモグロビンの同時制御の経験的測定(1A)および動的モデル(1B)を示す図である。図1Aでは、網状赤血球分布が等確率密度等高線において実線で示され、全RBCの母集団が破線で示される。両パネルの起点へ突出する斜線は、母集団の平均細胞内ヘモグロビン濃度(MCHC)を表す。この線状のいずれかの位置にあるRBCは、MCHCと等しい赤血球ヘモグロビン濃度を有する。まず、高速動力学(β)により、各パネルの右上に示す典型的な大きい未熟網状赤血球についての体積およびヘモグロビンが減少する。次に、低速動力学(α)により、MCHC線に沿った体積およびヘモグロビンが減少する。生物学的プロセスには本質的にノイズがあるため、体積およびヘモグロビンの減少に必要な事象中の小さな不規則変動により、個々の赤血球ヘモグロビン濃度がMCHC線の周りでドリフトし、血球が除去されるときに臨界体積(図1Bのv)に達するまで、図1Bの挿入図に示すようにMCHC線の周りで大きさ(D)が変動し得る。 図1A及び1Bは、末梢循環における平均RBCの体積およびヘモグロビンの同時制御の経験的測定(1A)および動的モデル(1B)を示す図である。図1Aでは、網状赤血球分布が等確率密度等高線において実線で示され、全RBCの母集団が破線で示される。両パネルの起点へ突出する斜線は、母集団の平均細胞内ヘモグロビン濃度(MCHC)を表す。この線状のいずれかの位置にあるRBCは、MCHCと等しい赤血球ヘモグロビン濃度を有する。まず、高速動力学(β)により、各パネルの右上に示す典型的な大きい未熟網状赤血球についての体積およびヘモグロビンが減少する。次に、低速動力学(α)により、MCHC線に沿った体積およびヘモグロビンが減少する。生物学的プロセスには本質的にノイズがあるため、体積およびヘモグロビンの減少に必要な事象中の小さな不規則変動により、個々の赤血球ヘモグロビン濃度がMCHC線の周りでドリフトし、血球が除去されるときに臨界体積(図1Bのv)に達するまで、図1Bの挿入図に示すようにMCHC線の周りで大きさ(D)が変動し得る。 20人の健常者と、3種類の軽度貧血を有する患者、すなわち、慢性疾患(ACD)に伴う貧血患者11人、サラセミア形質(TT)の患者33人、および鉄欠乏性貧血(IDA)の患者27人についてのモデルパラメータの一連の6つの箱型図である。各箱の上縁および下縁は、75パーセンタイルおよび25パーセンタイルに位置する。中央値は箱内部の水平線により示される。垂直線は、箱からの距離が四分位距離の1.5倍未満であるデータ点まで延びる。より極値のデータ点はプラス(+)記号として示される。高速動力学はβ、低速動力学はα、不規則変動はD、およびクリアランス閾値はvにより特徴付けられる。 4か月後にIDAを発症した患者の全血球計算(CBC)の等高線図である。各図は、体積−ヘモグロビン含有量確率密度の35%、60%、75%、および85%を囲む等高線を示す。起点からの破線はMCHCを表す。円は、各赤血球の体積−ヘモグロビン含有量座標のMCHC線上への投影の平均値を示す。MCHC線に垂直な短い実線は、平均投影の85%に対応する線に沿った位置を標示する。図3Aは、患者がIDAを発症する116日前に測定された正常CBCを示す。計算されたP0.85(グレーで塗りつぶされている)は正常である。図3Bは、IDA検出の65日後および51日前に測定された正常CBCを示す。CBCが正常であっても、P0.85は異常である。パネルの図3Cは、IDAが診断されたときのCBCを示す。 4か月後にIDAを発症した患者の全血球計算(CBC)の等高線図である。各図は、体積−ヘモグロビン含有量確率密度の35%、60%、75%、および85%を囲む等高線を示す。起点からの破線はMCHCを表す。円は、各赤血球の体積−ヘモグロビン含有量座標のMCHC線上への投影の平均値を示す。MCHC線に垂直な短い実線は、平均投影の85%に対応する線に沿った位置を標示する。図3Aは、患者がIDAを発症する116日前に測定された正常CBCを示す。計算されたP0.85(グレーで塗りつぶされている)は正常である。図3Bは、IDA検出の65日後および51日前に測定された正常CBCを示す。CBCが正常であっても、P0.85は異常である。パネルの図3Cは、IDAが診断されたときのCBCを示す。 4か月後にIDAを発症した患者の全血球計算(CBC)の等高線図である。各図は、体積−ヘモグロビン含有量確率密度の35%、60%、75%、および85%を囲む等高線を示す。起点からの破線はMCHCを表す。円は、各赤血球の体積−ヘモグロビン含有量座標のMCHC線上への投影の平均値を示す。MCHC線に垂直な短い実線は、平均投影の85%に対応する線に沿った位置を標示する。図3Aは、患者がIDAを発症する116日前に測定された正常CBCを示す。計算されたP0.85(グレーで塗りつぶされている)は正常である。図3Bは、IDA検出の65日後および51日前に測定された正常CBCを示す。CBCが正常であっても、P0.85は異常である。パネルの図3Cは、IDAが診断されたときのCBCを示す。 90日以内の第2の正常CBCを有する患者の20の正常CBCと、90日後までにIDAと診断された患者の20の正常CBCについてのP0.85の箱型図である。P0.85は、実際の診断の90日前まで、75%の感度および100%の特異性でIDAを順調に予測する。さらなる詳細について、以下および実施例3「鉄欠乏性貧血の予測」を参照されたい。 TTを有する5件の訓練ケース、IDAを有する5件の訓練ケース、TTを有する28件の検査ケース、IDAを有する22件の検査ケースについてのDの分布を示す箱型図である。結果は、本方法が、小球性貧血の原因としてのTTとIDAとを区別できることを示す。 式3で説明するクリアランス確率を計算するために使用される投影距離(Δ)の概略図である。血球がMCHC線上に投影され、クリアランス確率は、この線に沿った投影点から閾値(v)までの距離の関数である。 健常者の適合(濃いグレー)および経験的(薄いグレー)定常状態体積−ヘモグロビン含有量確率分布の比較を示す3次元グラフである。図6Aは、MCHC線を通る垂直面に投影された図を示す(図1参照)。図6Bは、起点に向かって見た90度回転図である。 図6A及び6Bは、健常者の適合(濃いグレー)および経験的(薄いグレー)定常状態体積−ヘモグロビン含有量確率分布の比較を示す3次元グラフである。図6Aは、MCHC線を通る垂直面に投影された図を示す(図1参照)。図6Bは、起点に向かって見た90度回転図である。 200を超える別個のシミュレーションによる最適化されたパラメータ値および残差平方和目的関数により判定された適合度の、一連の6つのヒストグラムである。目的関数のより小さい値は、より良い適合を意味する。図7は、すべてのモデルパラメータが、十分に定義された最適近傍を有することを示す。目的関数については式5を参照されたい。 図2に示す「健常者」および「IDA」患者についてのP0.85の分布を示す箱型図である。これらの結果に基づき、P0.85の閾値0.121が選択されて、図3Dに示す患者の独立集合の検査に使用された。本図は、潜伏IDAの閾値を特定するための方法の例である。 図8B〜8Dは、図2に示す定常状態の健常者(各図の上パネル)およびIDA患者(各図の下パネル)についてのP0.85(8B)、P0.75(8C)、およびP0.90(8D)の分布を示すヒストグラムである。 図8B〜8Dは、図2に示す定常状態の健常者(各図の上パネル)およびIDA患者(各図の下パネル)についてのP0.85(8B)、P0.75(8C)、およびP0.90(8D)の分布を示すヒストグラムである。 図8B〜8Dは、図2に示す定常状態の健常者(各図の上パネル)およびIDA患者(各図の下パネル)についてのP0.85(8B)、P0.75(8C)、およびP0.90(8D)の分布を示すヒストグラムである。 図9A〜9Cは、表1A、Bに挙げる単一RBC体積およびヘモグロビン含有量動力学の異なる関数形式の3つについての速度場を示す図である。黒の斜線はMCHCに等しい一定のヘモグロビン濃度を示す。 図10A及び10Bは、表1A、Bで説明する異なるクリアランス関数の2つについてのクリアランス確率を示す3次元グラフである。図10Aは関数形式AおよびCを示し、図10Bは関数形式Bを示す。MCHC線は黒で示される。 20人の健常者についての表1A、Bからの関数形式Bに基づくモデルパラメータと、3種類の軽度貧血を有する患者、すなわち、慢性疾患(ACD)に伴う貧血患者11人、サラセミア形質の患者33人、および鉄欠乏性貧血の患者27人についてのモデルパラメータの変化の様子とを示す、一連の6つの箱型図である。各箱の上縁および下縁は、75パーセンタイルおよび25パーセンタイルに位置する。中央値は箱内部の水平線により示される。垂直線は、箱からの距離が四分位距離の1.5倍以内であるデータ点まで延びる。より極値のデータ点はプラス(+)記号として示される。高速動力学はβ、低速動力学はα、不規則変動はD、およびクリアランス閾値はvにより特徴付けられる。 20人の健常者についての表2からの関数形式Cに基づくモデルパラメータと、3種類の軽度貧血を有する患者、すなわち、慢性疾患(ACD)に伴う貧血患者11人、サラセミア形質(TT)の患者33人、および鉄欠乏性貧血(IDA)の患者27人の変化の様子とを示す、一連の6つの箱型図である。各箱の上縁および下縁は、75パーセンタイルおよび25パーセンタイルに位置する。中央値は箱内部の水平線により示される。垂直線は、箱からの距離が四分位距離の1.5倍以内であるデータ点まで延びる。より極値のデータ点はプラス(+)記号として示される。高速動力学はβ、低速動力学はα、不規則変動はD、およびクリアランス閾値はvにより特徴付けられる。 図13A及び13Bは、本発明の一部の実施形態による、MCHC線への投影時に体積およびヘモグロビン含有量がこの線に沿った平均投影位置の閾値百分率よりも低いサンプル中の、赤血球分画を判定するための例示的な方法ステップを示すブロック図である。 図13A及び13Bは、本発明の一部の実施形態による、MCHC線への投影時に体積およびヘモグロビン含有量がこの線に沿った平均投影位置の閾値百分率よりも低いサンプル中の、赤血球分画を判定するための例示的な方法ステップを示すブロック図である。 本発明で使用する例示的なコンピューティング装置を示すブロック図である。 患者の血液サンプルについての体積およびヘモグロビン含有量測定からのモデルパラメータを推測するプログラムのための例示的なコンピュータコードを示す図である。 式4の数値近似を使用して定常状態RBC分布を計算する、呼び出された関数「ssRBC」のための例示的なコンピュータコードを示す図である。 被験者の血液サンプルについてP0.85を計算するプログラム「calculateCompensation」のための例示的なコンピュータコードを示す図である。この関数は、他の関数をも呼び出すが、パラメータ計算コードから独立して実行され得る。
赤血球は、ほとんど解明されていないプロセスにより循環から除去される。これらおよ
びその他のRBC成熟およびクリアランスプロセスの根底をなす特定の分子および細胞機
構は未知であるが、このようなプロセスは病気の状態において変化し、この変化によって
これらの病気を有する患者に種々の差を生じさせる可能性がある。例えば、体積またはヘ
モグロビン濃度のような特定の細胞特性が、循環における赤血球寿命中に変化する場合、
何らかの病気を有する患者の赤血球の母集団にわたるこの変化の平均速度が、病気のない
患者の赤血球の母集団にわたるこの変化の平均速度と異なる可能性がある。
健康な成人では、1日に約2.5×1011の新しいRBCが骨髄から末梢循環へ放出され、ほぼ同じ数が除去される。したがって、循環する母集団を有する血球は変化し続けるが、健常者(および軽度の病気を有する患者)では、母集団の特性が非常に安定している。臨床において、血液中の赤血球の体積分画(ヘマトクリット)、平均RBC体積(MCV)、RBC体積(RDW)の変動係数、および平均赤血球ヘモグロビン質量(MCH)等の母集団特性が、全血球計算(CBC)において定期的に測定される(1)。最近では、非常に若い(時間または日数が)循環RBC(網状赤血球)を特定し特徴付けることができるようになっている(2)。RBCは、骨髄から放出された後、数日で、体積およびヘモグロビンが急速に減少する(3)。この高速相に続いて、より緩やかに減少する、はるかに長い時間があり(4〜7)、この間に体積およびヘモグロビンが同時制御される(8)。図1Aを参照されたい。
網状赤血球および全循環RBC(図1A)の確率分布の比較は、体積とヘモグロビン含
有量との相関が、血球の成熟とともに、網状赤血球母集団の初期相関係数約0.40から
全母集団の約0.85まで増加することを示す。したがって、関与する分子機構の多くは
未知であるが、図1Aおよび1Bの等濃度線に示すように、ヘモグロビン濃度が母平均赤
血球ヘモグロビン濃度(MCHC)に向かう傾向があるように、平均RBCが成熟するこ
とが明らかである。この同時制御の結果として、ヘモグロビン濃度の変動が、体積および
ヘモグロビン含有量の変動よりも小さくなる(8)。
本明細書で開示されるパラメータは、TTおよびIDA患者のRBCが、正常状態にあ
るときよりも、絶対的および相対的に、はるかに小さい体積および低いヘモグロビン含有
量で末梢に残ることを示す。この残存は、このような貧血の非効率的な赤血球生成に反応
したクリアランスの代償性遅延を反映し得る。したがって、RBCクリアランスの要因と
なる挙動を変化させ得る機構が存在しなければならない。TTおよびIDAのRBCクリ
アランスを健常者またはACD患者のものと比較すると、要因を特定する新しい手段を提
供することができる。vの変動が、健常者の
の変動よりもはるかに小さく(14)、これは要因がMCHC線上の位置と非常に相関し
ていることを提示する。
貧血
貧血は、血流中に存在するヘモグロビンの量が正常よりも少ない状態である。貧血を、
赤血球の大きさ(MCVまたは平均赤血球体積として測定される)に基づき、主に3つの
群に細分することができる。正球性貧血は、赤血球が正常な大きさ(すなわち、80〜1
00fL)であるときに診断される。大球性貧血(通常、B12欠乏により生じる)は、
赤血球の大きさが正常よりも大きいとき、すなわち、MCV>100fLであるときに診
断される。小球性貧血は、赤血球が正常よりも小さいとき、すなわち、MCV<80fL
のときである。
貧血予備軍の診断/貧血発症の予測
本明細書で示すモデルおよび方法は、クリアランス遅延の兆候を探すことにより、明ら
かな臨床的貧血が発症する前に、潜伏貧血または代償性貧血の患者を特定するための可能
な方法を提供する。この可能性は、正常CBCを持っていたが、少なくとも30日および
90日未満後に別の正常CBCまたは臨床的なIDAが生じた患者の独立集合において検
査された。各患者のサンプルについて、全赤血球の(v、h)座標がMCHC線上に投影
され、確率密度が、平均(P0.85)の85%未満でこの線に沿って統合された。図3
A〜Cは、1患者のCBCが正常から潜伏IDAおよび最終的にIDAへ進む様子を示す
。中央パネル(図3B)に示すCBCは、臨床的には目立たないものであるが、P0.8
が異常であり、発症してはいないが51日の間に治療することになる貧血を予測してい
る。図3Dは、健康なままである患者の20の正常CBCと30〜90日後にIDAを発
症した患者の20の正常CBCについてのP0.85の値を示す。P0.85の値は、7
5%の感度および100%の特異性でIDAを予測した。実施例3「鉄欠乏性貧血の予測
」を参照されたい。現在の標準治療は、すべてのCBCが「正常」であるため、この母集
団で0%の感度を有する。このモデルに基づく予測は、異なる源および異なる時間点から
の情報の複数の測定および種類の統合に依拠することの多い統計的回帰手法と対照的に、
一時の単一CBC測定のみに依拠する(15)。
したがって、説明した方法を使用して、貧血予備軍を診断することができ、または発症
の少なくとも30〜90日前、例えば、被験者がまだ正常なCBCを有する(例えば、被
験者のHCT、MCV、RDW、MCH、MCHC、RBC、およびHGBレベルがすべ
て正常範囲内にある)ときに、被験者の貧血の発症を予測することができる。
正常CBC値は、分析を行う研究所によって変化し得るが、一般に、正常値は以下のと
おりである。
方法は、被験者のサンプル中の各赤血球について、赤血球体積およびHGB含有量(例
えば、HGB質量または濃度)を判定するステップと、サンプル中の全赤血球にわたる体
積およびHGB含有量値の分布を計算するステップと、例えば図5に示すような各赤血球
体積およびヘモグロビンの変換を計算するステップと、サンプル中の全赤血球についての
このような変換の分布を計算するステップと、例えば、図3A〜Dおよび8A〜Dに示す
ような、変換された体積およびHGB含有量の閾値レベルより変換が低い赤血球の百分率
を判定するステップとを含む。一部の実施形態では、方法が、例えば、被験者の赤血球を
含むサンプル中で、母平均赤血球ヘモグロビン濃度(MCHC)を判定するステップと、
サンプル中の全赤血球にわたる体積およびHGB含有量値の分布を計算するステップと、
例えば図5に示すような各赤血球体積およびヘモグロビンの変換を計算するステップと、
サンプル中の全赤血球についてのこのような変換の分布を計算するステップと、変換が平
均変換の閾値百分率よりも低いサンプル中の赤血球分画を判定して、サンプル分画を提供
するステップと、サンプル分画を基準分画と比較するステップとを含み、基準分画よりも
低いサンプル分画があれば、被験者が貧血、例えばIDAを発症する危険があることを示
す。閾値百分率および基準レベルの値を、当技術で公知の方法を使用して選択することが
でき、この値は感度および特異性を最大化する値を表すことができる。一部の実施形態で
は、閾値百分率が母平均MCHCの70%、75%、80%、85%、90%、または9
5%である。一部の実施形態では、基準分画が約0.10、0.11、0.12、0.1
3、0.14、または0.15である。当業者は、最適な閾値百分率および基準分画を、
公知の統計的な方法を使用して容易に特定することができるだろう。
鑑別診断 サラセミア対鉄欠乏性貧血
小球性貧血の鑑別診断は、鉄欠乏性貧血(IDA)、サラセミア形質(TT)慢性疾患
に伴う貧血(ACD)、および他の原因による貧血を含む。
ACDは、通常、正常値の下限よりも20%未満までのヘマトクリットの減少と、正常
または高いフェリチンと、低いまたは正常な総鉄結合能とを伴う。減少したRBCクリア
ランス閾値(v)が、貧血を補う適応生理学反応を表す場合には、ACDについてここ
で見られる正常vがおそらく予想され、貧血自体が適応生理学反応を表し得る(18)
。慢性疾患を有する被験者に貧血があれば、通常、ACDの診断を提案する。通常ACD
と関連する慢性疾患として、自己免疫疾患、例えば、クローン病、全身性エリテマトーデ
ス、関節リウマチ、および潰瘍性大腸炎;癌、例えば、リンパ腫およびホジキン病;慢性
腎臓病;肝臓病、例えば、肝硬変;慢性感染症、例えば、細菌性心内膜炎、骨髄炎(骨感
染症)、HIV/エイズ、B型肝炎、またはC型肝炎がある。例えば、GardnerおよびBen
z Jr.、「Anemia of chronic diseases.」 In: Hoffman他編、Hematology: Basic Princi
ples and Practice.、第5版、Philadelphia、Pa: Elsevier Churchill Livingstone、200
8、37章を参照されたい。
IDAは通常、ヘマトクリットの減少、低いMCV、および低いフェリチンを伴う。軽
度のIDAは通常、正常値の下限よりも20%未満までのヘマトクリットの減少を伴う。
IDAの患者は、正常ヘマトクリットを含む正常MCVの既存の証拠を示し得る。
TTは通常、高ヘモグロビンA2分画、または1もしくは複数のアルファグロビン遺伝
子突然変異の存在と関連する。TTは、通常、正常値の下限よりも20%未満までのヘマ
トクリットの減少と、低いMCV、および正常なフェリチンを伴うことが多い。TTは、
世界中で最も一般にスクリーニングされている疾患の1つであるが、既存の診断方法は非
常に高価であるか、受け入れがたいほど診断精度が低く、偽陽性率が30%にのぼる(1
9、20)。本明細書に記載のモデルおよび方法は、IDAとTTと、場合によっては他
の原因の貧血とを区別する、新しく、おそらくはより正確な方法を提供する。
図2は、Dが、小球性貧血の最も一般的な2つの原因であるIDAとTTとを区別す
ることを示す。10の訓練サンプルを分析することにより、Dの例示的な閾値値0.0
045が設定された。その後、軽度IDAまたはTTの診断を明確に行うことのできる5
0の独立した患者サンプルが分析され、Dが計算された。図4は、このD閾値が98
%の診断精度を有しており、IDAの22ケースのうち22ケース、およびTTの28ケ
ースのうち27ケースを正確に特定し、他の公開された手法よりも6〜41%優れている
ことを示す(20)。実施例4「小球性貧血の鑑別診断」を参照されたい。
したがって、本明細書に記載の方法は、小球性貧血の被験者において、鉄欠乏性貧血(IDA)とサラセミア形質(TT)と、場合によっては他の原因の貧血とを、被験者のサンプル中の平均ヘモグロビン含有量の変動の大きさの値(D を判定するステップと、Dを基準値と比較するステップとを含む方法により、鑑別診断するステップを含む。基準値よりも高いDがあれば、被験者がIDAであることを示し、基準値よりも低いDがあれば、被験者がTTであることを示す。
赤血球生成促進剤(例えば、エリスロポエチン(EPO))および/または鉄補給剤療
法の最適化
本明細書に記載の方法を使用して、EPOまたは鉄補給剤療法を最適化することができ
る。例えば、方法は、EPO治療を受けている被験者のサンプル中の正規化臨界体積、v
を判定するステップを含み得る。
一部の実施形態では、EPOまたは鉄補給剤の投与量を調節して、vまたはDを所
望レベル、例えば、閾値レベルよりも高いレベルもしくは低いレベル、または所望範囲内
のレベルで維持する。例えば、本明細書で言及するように、正常な健常者のvは母平均
体積
の約80%、または、一般的な90fLのMCVについて約72fLである。したがって
、一部の実施形態では、方法が、EPOまたは鉄補給剤の投与量を調節して、被験者の
の70%よりも高いv、被験者の
の75%よりも高いv、被験者の
の75〜85%よりも高いvを維持するステップを含む。一部の実施形態では、方法が
、治療を調節して70%よりも高いvまたは70〜75%のvを維持するステップを
含む。一部の実施形態では、vまたはDを時間とともに監視して、被験者に対するE
POまたは鉄補給剤の投与量を調節して、vまたはDを選択された閾値よりも高くも
しくは低く、または所与の範囲内で維持する。
GI疾患の危険性が高い被験者の特定
IDAまたは他の原因不明の貧血は、結腸癌(16)等の消化器疾患および小児期の栄
養不良(17)を含む深刻な疾患の最初の発現であることが多い。IDAに関連する他の
GI疾患として、結腸直腸(例えば、結腸)癌、消化管潰瘍(胃潰瘍、消化性潰瘍、盲腸
潰瘍)、憩室炎、腸管虚血、胃癌、胃炎、食道炎、GIポリープ、炎症性腸疾患(クロー
ン病、潰瘍性大腸炎)、およびセリアック病がある。この方法により貧血を初期に検出ま
たは予測することによって、このような疾患により早期に対応することができる。したが
って、例えば、GI疾患(例えば、大腸内視鏡検査を使用して)または栄養評価について
、スクリーニングまたは勧告追加スクリーニングのために、本明細書の方法を使用するこ
とができる。適切な患者の精密検査、例えば大腸内視鏡検査をより早期に開始することが
できるようになり、また適切な治療、例えば鉄補給剤の処方をより早期に開始することが
できるようになる。
血液ドーピングの検出
特に有酸素運動における運動選手のパフォーマンスは、存在するRBCの数によって影響され得る。RBCが多いほど、肺から運動中の筋肉への酸素の運搬および供給能力が大きくなる。赤血球の一時的な急増は、血液ドーピングと呼ばれる輸血、またはエリスロポエチン、類似体、および模倣薬等の赤血球生成促進剤の使用によって得られる。血液ドーピングは、検出が困難なことで知られている。例えば、JelkmannおよびLundby、Blood、2011年9月1日、118(9)、2395-404、Segura他、「Current strategic approaches for the detection of blood doping practices」、Forensic Sci Int. 2011年8月31日、電子ジャーナルを参照されたい。本方法を使用して、血液ドーピングまたは赤血球生成促進剤の使用の可能性を示す、RBC体積およびHGB含有量の異常分布の存在を検出することができる。例えば、一部の実施形態では、平均ヘモグロビン含有量の変動の大きさの値(D 、および平均赤血球体積(D)の変動の大きさの値の一方または両方の値を、血液ドーピングまたは赤血球生成促進剤の使用の疑いのある被験者のサンプルにおいて判定し、基準値よりも高いDおよび/またはDがある場合、血液ドーピングの存在もしくは使用または赤血球生成促進剤の使用を示す。
サンプル、システム、ソフトウェア、および試験方法
本明細書に記載の方法は、RBCを損なうことなく保存する公知の採取方法(例えば、
適量の真空(抜出し)と、RBCを実質的な溶血なしに採取可能な十分な大きさの針、例
えば、少なくとも25g以上の針とを使用する採血)を用いて得られた末梢血液サンプル
を使用して実施される。採取後24時間、12時間、または6時間以内に測定が行われる
ことが好ましい。RBC体積(例えば、低角度(2°〜3°)の散乱検出を使用)とヘモ
グロビン質量または濃度(例えば、高角度(5°〜15°)の散乱検出を使用)の両方を
測定可能な当技術で公知の方法または装置を使用して、網状赤血球およびCBCの測定を
行うことができる。例示的な方法は、米国特許出願公開第2011/0,178,716
号、米国特許出願公開第2011/0,164,803号、米国特許出願公開第2011
/0,149,061号、米国特許出願公開第2011/0,077,871号、米国特
許出願公開第2011/0,070,606号、および米国特許出願公開第2011/0
,070,210号に記載されている。
一部の実施形態では、血液分析器、例えば、手動、半自動、または自動の血液分析器を
使用して測定が行われる。血液分析器の例は、当技術で公知であり、例えば、参照により
本明細書に組み込まれている、米国特許第5,017,497号、米国特許第5,266
,269号、米国特許第5,378,633号、米国特許第5,631,165号、米国
特許第5,812,419号、米国特許第6,228,652号、米国特許第6,524
,858号、米国特許第6,320,656号、米国特許第7,324,194号、米国
特許第7,981,681号、米国特許出願公開第2008/0,153,170号、米
国特許出願公開第US2008/0,158,561号、米国特許出願公開第S2008
/0,268,494号、米国特許出願公開第US2011/0,178,716号、米
国特許出願公開第2011/0,077,871号、および米国特許出願公開第2011
/0,070,606号に記載されている。本方法で有用な血液分析器は、当技術で公知
のあらゆる検出方法、例えば、フローサイトメトリー、光学もしくは画像ベース分析また
はインピーダンスベースを使用することができる。本方法で有用な血液分析器は、通常、
CBCの全パラメータを測定可能なものである。すなわち、分析器は、少なくとも赤血球
(RBC)の赤血球体積(CV)と、血球ヘモグロビン濃度(CHC)または血球ヘモグ
ロビン質量(CH)とを判定することができるべきである。
血液分析器のいくつかのモデルが、例えば、Abbott Laboratories
(Abbott Park、IL、United States)(例えばCell−D
yn Sapphire)、およびSiemens(Deerfield、IL、Uni
ted States)(例えば、Advia120または2120自動血液分析器)か
ら市販されている。他の製造業者として、Beckman Coulter、Inc.(
Fullerton、CA、United States)、東亞医用電子(株)(日本
、神戸)、Constitution Medical(Boston、MA)、および
HORIBA ABX Inc(Irvine、CA、United States)が
ある。
本発明はまた、臨床的なサンプルを測定するための検出モジュールと、コンピューティング装置とを備えた血液分析器システムを提供する。このコンピューティング装置は、検出モジュールと通信し、(i)母平均赤血球ヘモグロビン濃度(MCHC)、ならびに、サンプル分画を判定するための、変換された体積およびヘモグロビン含有量がすべての変換された体積およびヘモグロビン含有量の平均の閾値百分率よりも低いサンプル中の赤血球分画、(ii)ヘモグロビン含有量減少率の母集団内赤血球間での変動の大きさの値(D 、(iii)赤血球体積減少率の母集団内赤血球間での変動の大きさの値(D 、および/または(iv)正規化臨界体積(v)もしくは他のクリアランス閾値の1または複数を、検出モジュールの出力に基づいて判定するためのプログラミングを備える。一部の実施形態では、コンピューティング装置が、血液分析器からの入力を使用する別個のコンピュータである。一部の実施形態では、コンピューティング装置が血液分析器に組み込まれるか、または血液分析器の一部となっている。
一般に、検出モジュールは、分析のために赤血球、例えば全血を含むサンプルを保持す
るように構成された分析チャンバを備える。一部の実施形態では、検出モジュールが、サ
ンプルを分析するように構成されたフローサイトメータを備える。一部の実施形態では、
検出モジュールが、サンプルを分析するように構成された光学または画像分析器を備える
本明細書では、(i)母平均赤血球ヘモグロビン濃度(MCHC)、ならびに、サンプル分画を判定するための、変換された体積およびヘモグロビン含有量がすべての変換された体積およびヘモグロビン含有量の平均の閾値百分率よりも低いサンプル中の赤血球分画、(ii)ヘモグロビン含有量減少率の母集団内赤血球間での変動の大きさの値(D 、(iii)赤血球体積減少率の母集団内赤血球間での変動の大きさの値(D 、および/または(iv)正規化臨界体積(v)もしくは他のクリアランス閾値の1または複数を計算するためのプログラミングを備えたコンピュータ読取可能な媒体も提供される。
一部の実施形態では、血液分析器からのデータ(例えば、赤血球(RBC)の赤血球体積(CV)、平均赤血球体積(MCV)、血球ヘモグロビン濃度(CHC)または血球ヘモグロビン含有量(CH)、平均赤血球ヘモグロビン濃度(MCHC)、および平均赤血球ヘモグロビン含有量(MCH)、ならびにこれらの母集団統計)がコンピューティング装置によって受信され、コンピューティング装置は、(i)母平均赤血球ヘモグロビン濃度(MCHC)、ならびに、サンプル分画を判定するための、変換された体積およびヘモグロビン含有量がすべての変換された体積およびヘモグロビン含有量の平均の閾値百分率よりも低いサンプル中の赤血球分画、(ii)ヘモグロビン含有量減少率の母集団内赤血球間での変動の大きさの値(D 、(iii)赤血球体積減少率の母集団内赤血球間での変動の大きさの値(D 、および/または(iv)正規化臨界体積(v)もしくは他のクリアランス閾値の1または複数を計算するために、本明細書に記載のアルゴリズムを含むプログラミングを実行する。

一部の実施形態では、コンピューティング装置が血液分析器からの入力を使用する別個
のコンピュータである。一部の実施形態では、コンピューティング装置が血液分析器に組
み込まれるか、または血液分析器の一部となっている。
プログラミングは、物理的記憶装置または伝送媒体において提供され得る。命令を受信
するコンピューティング装置(例えば、血液分析器システムの一部である別個の装置また
は情報処理モジュール)は、アルゴリズムを実行し、かつ/または本方法から入手したデ
ータを処理する。コンピュータ読取可能な記憶媒体の例として、コンピュータ内部または
外部にあってもなくても、フロッピーディスク、磁気テープ、CD−ROM、ハードディ
スクドライブ、ROMもしくは集積回路、光磁気ディスク、またはPCMCIAカード等
のコンピュータ読取可能なカード等がある。ファイルに含まれる情報を、コンピュータ読
取可能な媒体に記憶することができ、「記憶」とは、後日、ローカルまたは遠隔ネットワ
ーク上のコンピュータによりアクセス可能および検索可能となるように情報を記録するこ
とを意味する。一部の実施形態では、本明細書に記載の方法が、サンプルが実施される(
run)度に自動的に実行される。
図14は、本明細書に記載した技術とともに使用可能なコンピュータ装置700および
携帯型コンピュータ装置750の例を示す。コンピューティング装置700は、例えば、
ラップトップ、デスクトップ、ワークステーション、パーソナルデジタルアシスタント、
サーバ、ブレードサーバ、メインフレーム、およびその他の適切なコンピュータを含む種
々の形のデジタルコンピュータを表すものとする。コンピューティング装置750は、例
えば、パーソナルデジタルアシスタント、携帯電話、スマートフォン、およびその他の同
様のコンピューティング装置を含む種々の形の携帯型装置を表すものとする。ここに示す
構成要素、それらの接続および関係、ならびに機能は例にすぎず、本書類に記載し、かつ
/または特許請求した技術の実施を限定するものではない。
コンピューティング装置700は、ラップトップ、デスクトップ、ワークステーション
、パーソナルデジタルアシスタント、サーバ、ブレードサーバ、メインフレーム、および
、血液分析器システムまたは装置に組み込まれるコンピュータを含むその他の適切なコン
ピュータを含む種々の形のデジタルコンピュータを表すものとする。コンピューティング
装置750は、例えば、パーソナルデジタルアシスタント、携帯電話、スマートフォン、
およびその他の同様のコンピューティング装置を含む種々の形の携帯型装置を表すものと
する。ここに示す構成要素、それらの接続および関係、ならびに機能は例にすぎず、本書
類に記載し、かつ/または特許請求した技術の実施を限定するものではない。
コンピューティング装置700は、プロセッサ702、メモリ704、記憶装置706
、メモリ704および高速拡張ポート710に接続する高速ユーザインターフェース70
8、ならびに低速バス714および記憶装置706に接続する低速ユーザインターフェー
ス712を備える。構成要素702、704、706、708、710、712のそれぞ
れは、種々のバスを使用して相互接続され、共通のマザーボードに取り付けられるか、ま
たは必要に応じて他の方法で取り付けられ得る。プロセッサ702は、コンピューティン
グ装置700内で実行するための命令を処理することができ、この命令には、メモリ70
4または記憶装置706に記憶された、例えば、高速ユーザインターフェース708に結
合されたディスプレイ716を含む外部入出力装置にGUIのための図形情報を表示する
旨の命令が含まれる。他の実施では、複数のメモリおよび複数のタイプのメモリとともに
、複数のプロセッサおよび/または複数のバスを必要に応じて使用することができる。ま
た、複数のコンピューティング装置700を接続することができ、各装置が必要な動作の
一部を提供する(例えば、サーババンク、ブレードサーバ群、またはマルチプロセッサシ
ステムとして)。
メモリ704は、コンピューティング装置700内で情報を記憶する。一実施では、メ
モリ704が1または複数の揮発性メモリユニットである。別の実施では、メモリ704
が1または複数の不揮発性メモリユニットである。メモリ704は、例えば、磁気ディス
クまたは光ディスクを含む別の形のコンピュータ読取可能な媒体であってもよい。
記憶装置706は、コンピューティング装置700に大容量記憶を提供することができ
る。一実施では、記憶装置706がコンピュータ読取可能な媒体であるか、またはコンピ
ュータ読取可能な媒体を含むことができる。コンピュータ読取可能な媒体としては、例え
ば、フロッピーディスク装置、ハードディスク装置、光ディスク装置、テープ装置、フラ
ッシュメモリもしくは他の同様の固体メモリ装置、またはストーレジエリアネットワーク
内の装置または他の構成を含む多くの装置がある。コンピュータプログラム製品を、情報
担体で明確に具体化することができる。コンピュータプログラム製品は、実行時に、例え
ば前述した方法を含む1または複数の方法を実施する命令を含むことができる。情報担体
は、例えば、メモリ704、記憶装置706、プロセッサ702のメモリ等を含む、コン
ピュータまたは機械読取可能な媒体である。
高速コントローラ708は、コンピューティング装置700の帯域幅集約動作を管理し
、低速コントローラ712は低帯域幅集約動作を管理する。このような機能の割当ては一
例にすぎない。一実施では、高速コントローラ708がメモリ704、ディスプレイ71
6(例えば、図形プロセッサまたはアクセラレータを通して)、および種々の拡張カード
(図示せず)を受けることのできる高速拡張ポート710に結合される。この実施では、
低速コントローラ712が記憶装置706および低速拡張ポート714に結合される。種
々の通信ポート(例えば、USB、Bluetooth(登録商標)、イーサネット、無
線イーサネット)を含み得る低速拡張ポートを、例えば、キーボード、ポインティングデ
バイス、スキャナ、または例えばネットワークアダプタを通した、例えばスイッチもしく
はルータを含むネットワークデバイスを含む1または複数の入出力装置に結合することが
できる。
コンピューティング装置700は、図示するようにいくつかの異なる形で実施可能であ
る。例えば、コンピューティング装置700を、標準サーバ720として、またはこのよ
うなサーバ群において複数回実施することができる。また、コンピューティング装置70
0をラックサーバシステム724の一部として実施することもできる。加えて、または代
わりに、コンピューティング装置700を、例えば、ラップトップコンピュータ722を
含むパーソナルコンピュータ内で実施することができる。一部の例では、コンピューティ
ング装置700の構成要素を、例えば、装置750を含む携帯型装置(図示せず)の他の
構成要素と組み合わせることができる。このような装置はそれぞれ、コンピューティング
装置700、750の1または複数を含むことができ、システム全体を、互いに通信する
複数のコンピューティング装置700、750から構成することができる。
コンピューティング装置750は、他の構成要素のうち、プロセッサ752、メモリ7
64、ならびに例えば、ディスプレイ754、通信ユーザインターフェース766、およ
びトランシーバ768を含む入出力装置を備える。装置750は、例えば、マイクロドラ
イブまたは他の装置を含む記憶装置を備えて、追加の記憶容量を設けることができる。構
成要素750、752、764、754、766、768のそれぞれは、種々のバスを使
用して接続され、構成要素のうちの複数を共通のマザーボードに取り付けるか、または必
要に応じて他の方法で取り付けることができる。
プロセッサ752は、メモリ764に記憶された命令を含む、コンピューティング装置
750内の命令を実行することができる。プロセッサを、別個の複数のアナログおよびデ
ジタルプロセッサを備えるチップのチップセットとして実施することができる。プロセッ
サは、例えば、ユーザインターフェースの制御、装置750により実行されるアプリケー
ション、および装置750による無線通信を含む、例えば装置750の他の構成要素の調
整を行うことができる。
プロセッサ752は、制御ユーザインターフェース758およびディスプレイ754に
結合されたディスプレイユーザインターフェース756を通してユーザと通信することが
できる。ディスプレイ754は、例えば、TFT LCD(薄膜トランジスタ液晶ディス
プレイ)もしくはOLED(有機発光ダイオード)ディスプレイ、または他の適切なディ
スプレイ技術とすることができる。ディスプレイユーザインターフェース756は、ディ
スプレイ754を駆動して図形および他の情報をユーザに提示するための適切な回路を備
えることができる。制御ユーザインターフェース758は、ユーザからのコマンドを受け
、このコマンドをプロセッサ752へ送信するために変換する。加えて、外部ユーザイン
ターフェース762は、装置750と他の装置との近領域通信を可能にするように、プロ
セッサ742と通信することができる。外部ユーザインターフェース762は、例えば、
一部の実施における無線通信、または他の実施における無線通信を行うことができ、複数
のユーザインターフェースを使用することもできる。
メモリ764は、コンピューティング装置750内の情報を記憶する。メモリ764を
、1もしくは複数のコンピュータ読取可能な媒体、1もしくは複数の揮発性メモリユニッ
ト、1もしくは複数の不揮発性メモリユニットの1または複数として実施することができ
る。拡張メモリ774を設けて、例えば、SIMM(シングルインラインメモリモジュー
ル)カードユーザインターフェースを含み得る拡張ユーザインターフェース772を通し
て装置750に接続することもできる。このような拡張メモリ774は、装置750のた
めの追加の記憶空間を提供することができ、または装置750のためのアプリケーション
もしくは他の情報を記憶することもできる。すなわち、拡張メモリ774は、前述したプ
ロセスを実施または補う旨の命令を含むことができ、確実な情報を含むこともできる。し
たがって、例えば、拡張メモリ774を装置750のセキュリティモジュールとして設け
ることができ、装置750の確実な使用を可能にする命令を、拡張メモリ774にプログ
ラムすることができる。加えて、例えば、特定情報をSIMMカード上にハッキング不可
能に配置することを含む追加の情報とともに、SIMMカードを通して確実なアプリケー
ションを提供することができる。
メモリは、例えば、後述するようにフラッシュメモリおよび/またはNVRAMメモリ
を含むことができる。一実施では、コンピュータプログラム製品が、情報担体で明確に具
体化される。コンピュータプログラム製品は、実行時に、例えば前述した方法を含む1ま
たは複数の方法を実施する命令を含むことができる。情報担体は、例えば、トランシーバ
768または外部ユーザインターフェース762によって受信可能な、例えば、メモリ7
64、拡張メモリ774、および/またはプロセッサ752のメモリを含む、コンピュー
タまたは機械読取可能な媒体である。
装置750は、必要な場合にデジタル信号処理回路を含み得る通信ユーザインターフェ
ース766を通して無線通信することができる。通信ユーザインターフェース766は、
例えば、特にGSM音声通話、SMS、EMS、もしくはMMSメッセージ、CDMA、
TDMA、PDC、WCDMA、CDMA2000、またはGPRSを含む種々のモード
またはプロトコルによる通信を行うことができる。このような通信は、例えば、高周波ト
ランシーバ768を通して行うことができる。加えて、例えば、Bluetooth(登
録商標)、WiFi、または他のこのようなトランシーバ(図示せず)の使用を含む短距
離通信を行うことができる。加えて、GPS(全地球測位システム)レシーバモジュール
770は、装置750で実行されているアプリケーションにより必要に応じて使用可能な
追加のナビゲーション関連および位置関連無線データを装置750に提供することができ
る。
装置750は、ユーザから話される情報を受信して、使用可能なデジタル情報に変換す
ることのできる音声コーデック760を使用して音声通信することができる。音声コーデ
ック760は、同様に、例えば、装置750の送受話器内のスピーカを通して、ユーザに
対して可聴音を発生することができる。このような音は、音声通話からの音を含むことが
でき、録音された音(例えば、音声メッセージ、音楽ファイル等)を含むことができ、ま
た装置750で動作するアプリケーションにより発生する音を含むことができる。
コンピューティング装置750を、図示したように、複数の異なる形で実施することが
できる。例えば、コンピューティング装置750を携帯電話780として実施することが
できる。コンピューティング装置750を、スマートフォン782、パーソナルデジタル
アシスタント、またはその他の同様の携帯型装置の一部として実施することもできる。
本明細書に記載のシステムおよび技術の種々の実施を、デジタル電子回路、集積回路、
特別に設計されたASIC(特定用途向け集積回路)、コンピュータハードウェア、ファ
ームウェア、ソフトウェア、および/またはこれらの組合せにおいて実現することができ
る。これらの種々の実施は、少なくとも1つのプログラム可能プロセッサを含むプログラ
ム可能システムで実行可能かつ/または解釈可能な1または複数のコンピュータプログラ
ムにおける実施を含むことができる。このシステムは、専用または汎用とすることができ
、記憶システム、少なくとも1つの入力装置、および少なくとも1つの出力装置からデー
タおよび命令を受信し、かつこれらに対してデータおよび命令を送信するように結合され
る。
このようなコンピュータプログラム(プログラム、ソフトウェア、ソフトウェアアプリ
ケーション、またはコードとしても知られる)は、プログラム可能プロセッサのための機
械命令を含み、高レベル手続型および/またはオブジェクト指向プログラミング言語、な
らびに/またはアセンブリ/機械言語において実施可能である。本明細書で使用されるよ
うに、機械読取可能な媒体およびコンピュータ読取可能な媒体という用語は、機械命令を
受信する機械読取可能な媒体を含むプログラム可能なプロセッサに機械命令および/また
はデータを与えるために使用されるコンピュータプログラム製品、装置および/またはデ
バイス(例えば、磁気ディスク、光ディスク、メモリ、プログラム可能論理回路(PLD
))を指す。
ユーザとの対話を行うために、本明細書に記載のシステムおよび技術を、情報をユーザ
に表示するためのディスプレイ装置(例えば、CRT(陰極線管)またはLCD(液晶デ
ィスプレイ)モニタ)、ユーザがコンピュータに入力可能なキーボードおよびポインティ
ングデバイス(例えば、マウスまたはトラックボール)を有するコンピュータで実施する
ことができる。他の種類の装置を使用して、ユーザとの対話を行ってもよい。例えば、ユ
ーザに与えられるフィードバックは、感覚フィードバック(例えば、視覚フィードバック
、聴覚フィードバック、もしくは触覚フィードバック)の形とすることができ、ユーザか
らの入力を、音響、音声、または触覚入力を含む形で受信することができる。
本明細書に記載のシステムおよび技術を、演算システムで実施することができ、この演
算システムは、バックエンドコンポーネント(例えば、データサーバとして)、ミドルウ
ェアコンポーネント(例えば、アプリケーションサーバ)、フロントエンドコンポーネン
ト(例えば、本明細書に記載のシステムおよび技術の実施とユーザが対話する際に通す、
グラフィカルユーザインターフェースまたはウェブブラウザを有するクライアントコンピ
ュータ)、またはこのようなバックエンドコンポーネント、ミドルウェアコンポーネント
、もしくはフロントエンドコンポーネントの組合せを含む。システムのコンポーネントを
、デジタルデータ通信(例えば、通信ネットワーク)の形または媒体により相互接続する
ことができる。通信ネットワークの例として、ローカルエリアネットワーク(LAN)、
広域ネットワーク(WAN)、およびインターネットがある。
演算システムは、クライアントおよびサーバを含むことができる。クライアントおよび
サーバは、一般に互いにリモートであり、通常は通信ネットワークを通じて対話する。ク
ライアントおよびサーバの関係は、各コンピュータで実行され、互いにクライアント−サ
ーバ関係を有するコンピュータプログラムによって発生する。
一部の実施では、本明細書に記載のエンジンが、単一または複合エンジンと別個であっ
ても、組み合わされていても、組み込まれていてもよい。図示したエンジンは、本明細書
に記載のシステムを、図示したソフトウェアアーキテクチャに限定するものではない。
実施例
本発明をさらに以下の実施例で説明する。この実施例は、特許請求の範囲に記載された
本発明の範囲を限定するものではない。
モデル開発
生体内の個々のRBCの寿命中における体積およびヘモグロビン制御は非常に複雑で、
解明が困難である。RBCの大きな母集団の平均挙動を解明するほうが扱いやすい。この
母集団レベル挙動を見抜くために、RBC母集団の動力学を説明するRBC成熟およびク
リアランスのモデルが開発された。モデルは、時間(t)にわたる平均RBCの体積(v
)およびヘモグロビン(h)動力学を、決定論的還元(f)と、特定の関数形式が変動し
得る不規則変動(ζ)とに分解する。式1に一例を示し、式中、vおよびhは母平均
によってスケーリングされ、tは平均赤血球齢
によってスケーリングされる。先のレポート(3〜5、7)からのデータに基づき、2つ
のパラメータが、決定論的コンポーネント、すなわちRBCがMCHC線に近くなるまで
その効果が優位である高速変化(β)と、低速変化(α)とに挿入された。図1Bおよび
式1に示すように、平均ゼロおよび分散が拡散テンソル2Dにより与えられた、ガウスま
たは同様の分布を持つ不規則変数として、不規則変動をモデリングすることができる。
一般的な逆問題および特別なヒト病態生理学(9)のように、この問題は、fの異なる
関数形式が生体内動力学を再現するという意味で不良設定である。fの正確な関数形式は
重要ではない。このモデルの挙動は高速および低速の決定論的動力学および不規則変動の
定性的組合せに依拠した。表1A〜1Bは、モデルで使用可能な異なるタイプの関数形式
を示す。
図9A〜Cは、fの形式A〜Cに対応する速度場を示し、図10A、Bは、クリアランス
関数(d)形式AおよびBを示す。τは母集団内の赤血球の平均齢である。Δの定義を以
下に述べる。
パラメータが、決定論的進化(f)およびクリアランス関数(d)の異なる関数形式に
ついて推定された。定性的結果はこれらの異なる関数形式について一致しており、これら
の結果が生体内病態生理学の特徴を表し、データの過剰適合がないことを表すことを提示
する。関数形式の詳細は、表1A〜Bに示され、形式A〜Cについての推定は、図2、1
1、および12の箱型図に示される。
本モデルでは、典型的な個々の赤血球についての体積およびヘモグロビン含有量の不規
則変動および決定論的散逸または還元が、ポテンシャル(10)中のブラウン運動をモデ
リングするために通常使用されるランジュバン方程式により説明された。これにより、R
BCの全循環母集団の動力学を、フォッカープランク方程式(10、11)により近似可
能な時間依存性体積−ヘモグロビン確率分布(P(v,h,t))についてのマスタ方程
式によって説明することができる。式2は、この体積−ヘモグロビン分布についての確率
密度のドリフト(f)、拡散(D)、出生(b)、および死滅(d)を説明する。
出生死滅プロセスは、常に母集団に加えられ母集団から除去されるRBCの要因となる
。健常および軽度の病気の状態では、加えられる赤血球の総数は、除去される赤血球の総
数に等しい。
RBC除去の正確な要因および機構は完全には解明されていない(12)が、図1Aに示
すような経験的測定により、MCHC線に沿った閾値(v)があり、大部分のRBCが
MCHC線を超えて除去されたことが提示される。
経験的RBC分布の観察に基づき、RBCクリアランスの確率が、RBC体積およびヘ
モグロビン含有量の関数としてモデリングされる。体積−ヘモグロビン含有量平面におけ
る各RBCの位置がMCHC線上に投影され、クリアランス確率(d)が様々に定義され
、例えば、この線状のこの投影点から閾値vまでの距離のシグモイド(図10A)また
はステップ(図10B)関数として定義される。図1、5および表1Bを参照されたい。
式3は、この関係を定量化する。
CBC測定は人によって異なるが、これらの動力学プロセスが生体内定常状態Pに達
することを示す健常者(13)については大きく変化しない。ここで
所与のパラメータ集合について、病態生理学的範囲外の体積およびヘモグロビン含有量
でゼロになる確率の境界条件と、経験的に測定された網状赤血球分布に等しい初期条件と
を有する一次空間導関数
および二次空間導関数
についての有限差分近似を使用して、式2を数値的に解決することができる。体積−ヘモ
グロビン含有量平面は、一定のメッシュ幅で離散化され、各メッシュ赤血球中の確率密度
に等しい変数のベクトル(P)として表される。ここで報告されるシミュレーション結果
は、体積軸に沿った1.8fLおよびヘモグロビン含有量軸に沿った1.8pgのメッシ
ュ幅により実行された。このメッシュ幅は、経験的体積およびヘモグロビン含有量測定の
解析解と比較可能であった。より小さいメッシュ幅(1.2fLおよび1.2pg)を使
用してすべての結果が確認された。空間導関数の風上有限差分近似を使用して、対流寄与
(f)が数値的にモデリングされた。結果として得られる常微分方程式の線形システムが
、MATLABode15s積分器を使用して積分され、定常状態(P)に達するまで
反復された。
数値問題(P)についての定常状態分布は、解析的にも判定された。進化(J+L)
項およびクリアランス(d)項の数値近似は線形演算子であり、出生プロセスの積分スケ
ーリングは、平均齢の2倍の逆数に等しい定数である。
線形演算子を反転させて、式4に示すように、平均赤血球齢により示される定常状態分布
族を生じさせることができる。
解析および有限差分手法により得られた定常状態分布間の差はごくわずかであった。
パラメータ(α、β、D、およびv)の適切な選択により、本明細書に記載のモデル
は、健常者のRBC母集団の観察された分布を正確に再現する。
検証
モデルが健常者のRBC母集団の動力学を貧血患者のRBC母集団の動力学と区別する
ことができるか否かを検査するために、原因疾患が異なる3種類の最もよく見られる貧血
、すなわち炎症性疾患である慢性疾患に伴う貧血(ACD)、遺伝性疾患であるサラセミ
ア形質(TT)、および栄養性疾患である鉄欠乏性貧血(IDA)を有する患者について
、CBCおよび網状赤血球測定を行った(14)。RBC母集団特性が安定して見え、準
定常状態仮定が適正な各貧血の軽度のケースが、複数の見かけの健常対照群として選択さ
れた。
血液サンプルおよびCBC結果は、Partners Healthcare Ins
titutional Review Boardにより認可された研究プロトコルの下
で、三次医療成人病院の臨床検査室から得られた。網状赤血球およびCBC測定が、Si
emens Advia 2120自動血液分析器で、採取後6時間以内に行われた(2
1)。
IDAは、正常値の下限よりも20%未満までのヘマトクリットの軽度の減少、低いM
CV、低いフェリチン、および正常ヘマトクリットを含む正常MCVの既存の証拠として
定義される。急性疾患、急性出血、輸血後6か月以内、並行入院、慢性炎症性疾患、また
は異常ヘモグロビン症の患者は除外された。
TTは、高ヘモグロビンA2分画、またはアルファグロビン遺伝子突然変異の存在につ
いてのヘテロ接合体として定義されるとともに、正常値の下限よりも20%未満までのヘ
マトクリットの減少と、低いMCV、および正常なフェリチンとして定義された。急性疾
患、急性出血、輸血後6か月以内、並行入院、慢性炎症性疾患、またはさらなる異常ヘモ
グロビン症の患者は除外された。
ACDは、正常値の下限よりも20%未満までのヘマトクリットの減少と、正常または
高いフェリチンと、低いまたは正常な総鉄結合能として定義された。急性疾患、急性出血
、輸血後6か月以内、並行入院、慢性炎症性疾患、または異常ヘモグロビン症の患者は除
外された。
各患者のサンプルに対して、患者について観察された定常状態を再現する最適なパラメ
ータ集合(α、β、D、およびv)が特定された。シミュレーションした定常状態分布
と測定したCBC分布との最小二乗適合を使用して、最良適合を特定した。同一の患者に
ついて繰返し検査を行うことができるが、適合パラメータの変動はCBC測定の解析変動
によって説明されることがわかった。
勾配または非勾配最適化方法を使用して、パラメータ空間の最適近傍を各患者について
特定した。患者の経験的に測定された網状赤血球分布および初期のランダムに選択された
パラメータ集合を、始点として使用した。式4を使用して、これにより生じる定常状態R
BC分布を計算した。この計算された分布(P)を、その後、経験分布(PCBC)と
比較した。パラメータ推定の質が、式5に示すように、離散化分布に対する正規化残差平
方和に等しい目的関数を演算することにより定量化された。式中iおよびjは、離散化体
積−ヘモグロビン平面の赤血球の指標を表す。
一般に、目的関数は、2つの分布間の相違(または類似)の程度の測定値を提供する。
目的関数の値が閾値条件を満たすか否かに基づいて、パラメータ値を調節することができ
る。例えば、目的関数の値が閾値よりも高くなるか、または低くなるまでパラメータ値を
調節することができる。目的関数が相違の測定値(例えば平均二乗差、絶対差の和等)を
表す場合、目的関数は、パラメータを調節することによって減少することが求められる。
あるいは、目的関数が類似の測定値(例えば相関係数、相互情報等)を表す場合、目的関
数は、パラメータを調節することによって増加することが求められる。この例では、この
目的関数を減少させるために新しいパラメータ値が選択された。MATLABの勾配に基
づく(lsqnonlin関数)、ならびに非勾配に基づく(fminsearchおよ
びpatternsearch関数)最適化アルゴリズムを使用して、最適パラメータを
検索した。すべてのパラメータが非負に制限され、均一な初期パラメータ空間を定義して
、1000日よりも長いか、または5日よりも短い平均赤血球齢を除外した。その後、ラ
テン超方格サンプリングを使用して、この空間から初期パラメータを選択した。1つの最
適化制限が課され、αを平均赤血球齢が5日よりも長くなるのに十分な小ささに限定し、
平均赤血球齢が1000日よりも短くなるのに十分な大きさに限定する。図6は、モデル
が、シミュレーションされた定常状態確率分布と測定された定常状態確率分布とを比較す
ることにより、この健常者についてのPCBCを正確に再現することを示す。MCHC線
に沿って投影されると、この患者についての経験分布はモード近くでわずかに高い密度を
有し、モードのいずれかの側までの領域でわずかに低い密度を有していた。
図7は、1人の患者に対する200を超える最適化から得られたパラメータの最小の局
所解を示す。一部の結果は、軸の範囲よりも高い局所解を有した。すべてのシミュレーシ
ョン間の最良適合は、全パラメータについての値の小近傍を形成し、パラメータ推定プロ
セスが、この患者の血液サンプルについての明確に定義された最適近傍に達したことを示
す。
健常者および貧血患者についての適合パラメータが図2に示され、中央値適合パラメー
タが表2に挙げられる。
健常者から得た最良適合パラメータと貧血患者から得た最良適合パラメータとの間には
明白な差があり、異なる貧血疾患は異なる特性パラメータ集合を有した。例えば、健常者
およびACD患者は、高いβおよびβと低いαとを有し、すなわち、低速相中よりも
高速相中に体積およびヘモグロビンが比較的多く失われた。これに対して、TTおよびI
DA患者は、高速相中よりも低速相中に体積およびヘモグロビンが比較的多く失われた。
ACD患者は、健常者に対してわずかに高いDとわずかに低いDとを示したが、TT
はDの実質的な減少に伴うより大きなDの増加に関連した。IDA患者は健常者と同
様のDを有し、大部分の患者が健常者の10倍を超える、Dの劇的な増加を示した。
健常者およびACD患者の正規化臨界体積vは、
の約80%、または約72fLである。TTまたはIDA患者の大半は、通常、
および
が減少した。図2は、
および
の絶対的な減少に加えて、これらの患者のvがさらに減少下し、異なる患者によっては
るかに大きなばらつきがあったことを示す。
鉄欠乏性貧血の予測
本実施例は、平均よりも起点に近いMCHC線に沿って投影する赤血球の拡大母集団に
基づいて、代償性または潜伏IDAを、場合によっては、現在可能であるよりも少なくと
も90日早く予測することができるという仮説を検査した。投影動作を図5に示す。まず
、体積(v)およびヘモグロビン(h)を有する各赤血球のMCHC線に沿った投影位置
(u)が判定された。
MCHC線に沿った閾値は、平均投影位置の比率(φ)
として定義された。次に、起点とこの閾値との間に位置する投影された赤血球分画が以下
のように計算された。式中、fMCHCは、MCHC線上の位置の関数としての投影され
た赤血球の確率密度である。
図2で使用された定常状態CBCについての異なる閾値の識別効率を比較することによ
り、φ=85%の閾値が選択された。図8A〜Bに示す85%の閾値は、他の閾値(70
%、75%(図8Cに示す)、80%、90%(図8D)、および100%を含む)と比
較したときに最大の分離を生じさせた。P0.85の閾値0.121が、この訓練集合に
基づいて選択された。
すべて正常な、新しい独立した患者の40のCBCが特定された。これらの正常CBC
の20が、30〜90日後にさらなる正常CBCを有した患者からのものであり、これら
の正常CBCの20が、90日未満の後にIDAを発症した患者からのものである。2つ
のCBC間に急性出血または鉄補給のある患者は除外された。IDAの定義は、実施例1
で前述した。図3A〜Dは、P0.85の閾値0.121が、この独立した被験群の中で
75%の感度および100%の特異性を有するIDAを予測することができることを示す
小球性貧血の鑑別診断
小球性貧血の最も一般的な2つの原因を区別する際のDの診断精度を評価するために
、図4に示すように、パラメータを10の訓練ケース、すなわちIDAの5ケース、TT
の5ケースに適合させた。ケースの定義は、実施例2で前述した。TTの訓練集合は、1
.7×10−15〜2.3×10−15のDを有した。IDAの訓練集合は、0.00
9〜0.043のDを有した。0.0045に等しい閾値、IDA訓練集合の最小D
およびTTの訓練集合の最大Dの平均が選択された。新しい独立した50のケース、す
なわちIDAの22ケース、TTの28ケースがこれにより特定された。22のIDAケ
ースすべてが正確に分類され、28のTTケースのうち27ケースが正確に分類され、総
合診断精度は98%であった。この精度は、4つの通常挙げられる判別関数の精度よりも
優れている(20)。すなわち、Green&King式の精度は92%、Micro/
Hypoは84%、Mentzerは68%、およびEngland&Fraserは5
7%である。このような他の判別関数のさらなる詳細については参考文献(20)を参照
されたい。
(参考文献)

他の実施形態
本発明をその詳細な説明とともに説明したが、前述した説明は例示的なものであり、添
付の特許請求の範囲によって定義される本発明の範囲を限定するものではないことを理解
されたい。その他の態様、利点、変更が、以下の特許請求の範囲に含まれる。
700、750 コンピューティング装置
702、742、752 プロセッサ
704、764 メモリ
706 記憶装置
708 高速ユーザインターフェース
710 高速拡張ポート
712 低速ユーザインターフェース
714 低速バス
720 標準サーバ
722 ラップトップコンピュータ
724 ラックサーバシステム
754 ディスプレイ
756 ディスプレイユーザインターフェース
758 制御ユーザインターフェース
760 音声コーデック
762 外部ユーザインターフェース
766 通信ユーザインターフェース
768 トランシーバ
772 拡張ユーザインターフェース
774 拡張メモリ
780 携帯電話
782 スマートフォン

Claims (6)

  1. 被験者のエリスロポエチン(EPO)または他の赤血球生成促進剤(ESA)の投与量を最適化するための方法であって、
    EPOまたはESA治療を受けている前記被験者のサンプル中の正規化臨界体積(v)を判定するステップであって、前記正規化臨界体積は、母平均赤血球ヘモグロビン濃度(MCHC)線に沿った閾値を示し、当該MCHC線は、前記サンプル中の複数の赤血球の体積及びヘモグロビン含有量によって、少なくとも部分的に規定される、正規化臨界体積(v )を判定するステップと、
    前記v を基準と比較するステップとを含む方法。
  2. 前記被験者の前記母平均体積
    の75%よりも大きいv、または78〜82%のvを維持するように、前記EPOの投与量が調節される、請求項1に記載の方法。
  3. 前記サンプルは、前記被験者の全血を含む、請求項1又は2に記載の方法。
  4. 体積およびヘモグロビン含有量座標が、すべての体積およびヘモグロビン含有量座標の最小二乗線形適合を表す線上に投影されるように、個々の赤血球の体積およびヘモグロビン含有量が変換される、請求項1〜のいずれか一項に記載の方法。
  5. 前記v を判定するステップは、複数のパラメータを含む数値モデルに基づきv を判定するステップを含み、前記複数のパラメータは、前記v 、前記複数の赤血球の体積及びヘモグロビン含有量の低速動力学(α)の変化を示すパラメータ、前記複数の赤血球の体積及びヘモグロビン含有量の高速動力学(β)の変化を示すパラメータ、及び前記複数の赤血球の体積及びヘモグロビン含有量中の変動の大きさ(D)を示すパラメータ、を少なくとも含む、請求項1に記載の方法。
  6. 前記EPOまたはESAの投与量が前記比較に基づき調節される、請求項1に記載の方法。
JP2017144114A 2010-09-16 2017-07-26 診断のための赤血球動力学 Active JP6522059B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US38335710P 2010-09-16 2010-09-16
US61/383,357 2010-09-16

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2015189343A Division JP6185030B2 (ja) 2010-09-16 2015-09-28 診断のための赤血球動力学

Publications (2)

Publication Number Publication Date
JP2018009992A JP2018009992A (ja) 2018-01-18
JP6522059B2 true JP6522059B2 (ja) 2019-05-29

Family

ID=45832278

Family Applications (3)

Application Number Title Priority Date Filing Date
JP2013529382A Active JP5868982B2 (ja) 2010-09-16 2011-09-16 診断のための赤血球動力学
JP2015189343A Active JP6185030B2 (ja) 2010-09-16 2015-09-28 診断のための赤血球動力学
JP2017144114A Active JP6522059B2 (ja) 2010-09-16 2017-07-26 診断のための赤血球動力学

Family Applications Before (2)

Application Number Title Priority Date Filing Date
JP2013529382A Active JP5868982B2 (ja) 2010-09-16 2011-09-16 診断のための赤血球動力学
JP2015189343A Active JP6185030B2 (ja) 2010-09-16 2015-09-28 診断のための赤血球動力学

Country Status (6)

Country Link
US (2) US9938557B2 (ja)
EP (2) EP2616818B1 (ja)
JP (3) JP5868982B2 (ja)
CN (1) CN103430028A (ja)
IL (2) IL225275A (ja)
WO (1) WO2012037524A2 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103430028A (zh) 2010-09-16 2013-12-04 通用医疗公司 供诊断用的红细胞动力学
US20170108487A1 (en) * 2014-06-05 2017-04-20 The General Hospital Corporation Predicting morbidity associated with red blood cell volume variance
US20160161510A1 (en) * 2014-12-05 2016-06-09 C A Casyso Ag Blood Testing System Result Interpreter Interface and Methods
WO2017106461A1 (en) 2015-12-15 2017-06-22 The General Hospital Corporation Methods of estimating blood glucose and related systems
US11293852B2 (en) 2016-04-07 2022-04-05 The General Hospital Corporation White blood cell population dynamics
EP3586145A1 (en) 2017-02-21 2020-01-01 Abbott Diabetes Care Inc. Methods, devices, and systems for physiological parameter analysis
WO2021035038A1 (en) * 2019-08-20 2021-02-25 The General Hospital Corporation Single-cell modeling of clinical data to determine red blood cell regulation
JPWO2022102734A1 (ja) * 2020-11-13 2022-05-19
KR102517894B1 (ko) * 2021-02-23 2023-04-04 연세대학교 산학협력단 당뇨발에 의한 하지 절단 수술의 치료 반응성을 예측하는 방법

Family Cites Families (46)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU683749A1 (ru) * 1973-12-14 1979-09-05 Л. А. Тиунов, А. И. Клиорнн и Т. С. Колосова Способ количественного определени скорости распада гемоглобина в крови
NL8601000A (nl) 1986-04-21 1987-11-16 Jan Greve T H Twente Afdeling Het gebruik van gepolariseerd licht in stromingscytometrie.
JPH07119757B2 (ja) 1991-03-31 1995-12-20 日本光電工業株式会社 血球計数装置
WO1993016384A1 (en) 1992-02-07 1993-08-19 Abbott Laboratories Method for accurately enumerating and sensitively qualifying heterogeneous cell populations in cytolytic processing conditions
US5369014A (en) * 1992-08-17 1994-11-29 Beth Israel Hospital Assoc. In-vitro method for determining a surreptitious use of exogenous erythropoiesis stimulating agents by a normal living subject
US5956501A (en) 1997-01-10 1999-09-21 Health Hero Network, Inc. Disease simulation system and method
JP3703858B2 (ja) * 1993-09-30 2005-10-05 シスメックス株式会社 帰属度判別装置
US5631165A (en) 1994-08-01 1997-05-20 Abbott Laboratories Method for performing automated hematology and cytometry analysis
US5812419A (en) 1994-08-01 1998-09-22 Abbott Laboratories Fully automated analysis method with optical system for blood cell analyzer
JPH11326315A (ja) 1998-05-21 1999-11-26 Sysmex Corp β−サラセミアの鑑別方法
US6030838A (en) * 1998-09-24 2000-02-29 Emt & Associates, Inc. Hematological parameter
US6228652B1 (en) 1999-02-16 2001-05-08 Coulter International Corp. Method and apparatus for analyzing cells in a whole blood sample
CA2367780A1 (en) 1999-03-31 2000-10-05 Bayer Corporation Single channel, single dilution detection method
WO2001077140A2 (en) 1999-10-19 2001-10-18 Srienc Friederich Measurement of nutrient uptake in cells and methods based thereon
US6320656B1 (en) 2000-02-18 2001-11-20 Idexx Laboratories, Inc. High numerical aperture flow cytometer and method of using same
US6784981B1 (en) 2000-06-02 2004-08-31 Idexx Laboratories, Inc. Flow cytometry-based hematology system
EP1425589B1 (en) 2001-09-14 2007-12-05 F. Hoffmann-La Roche Ag Differential diagnosis of disorders of iron metabolism
US6916662B2 (en) * 2003-01-31 2005-07-12 Abbott Laboratories Performance improvement for hematology analysis
US8614204B2 (en) 2003-06-06 2013-12-24 Fibrogen, Inc. Enhanced erythropoiesis and iron metabolism
FR2883972B1 (fr) 2005-03-31 2007-11-16 C2 Diagnostics Sa Procede pour l'analyse d'un echantillon de sang et appareil et reactif pour sa mise en oeuvre
US7788045B2 (en) 2005-09-01 2010-08-31 Meditasks, Llc Systems and method for homeostatic blood states
US7586589B2 (en) * 2005-09-24 2009-09-08 Beckman Coulter, Inc. Methods of determination of responsiveness to erythropoietin treatment
US8481323B2 (en) 2005-10-28 2013-07-09 Hewlett-Packard Development Company, L.P. Systems and methods for measuring glycated hemoglobin
JP4925364B2 (ja) * 2006-01-20 2012-04-25 ベックマン コールター, インコーポレイテッド 鉄欠乏を検出する方法
US7674622B2 (en) 2006-12-22 2010-03-09 Abbott Laboratories, Inc. Method for determination of nucleated red blood cells and leukocytes in a whole blood sample in an automated hematology analyzer
US7804594B2 (en) 2006-12-29 2010-09-28 Abbott Laboratories, Inc. Method and apparatus for rapidly counting and identifying biological particles in a flow stream
JP4949898B2 (ja) 2007-03-09 2012-06-13 シスメックス株式会社 血球分析装置
NZ555163A (en) 2007-05-14 2010-05-28 Fonterra Co Operative Group Methods of immune or hematological enhancement, inhibiting tumour formation or growth, and treating or preventing cancer, cancer symptoms, or the symptoms of cancer treatments
JP5025371B2 (ja) 2007-07-31 2012-09-12 シスメックス株式会社 血液分析装置
CN101158691A (zh) * 2007-11-16 2008-04-09 东南大学 糖化血红蛋白含量的电化学检测方法及检测装置
EP2085486A1 (en) 2008-02-01 2009-08-05 Assistance Publique - Hopitaux de Paris Methods and kits for the rapid determination of patients at high risk of death during septic shock
CN102027368B (zh) 2008-03-21 2014-02-26 艾博特健康公司 利用红细胞内含有的血红蛋白的本征色素沉着来确定血样的红细胞指数的方法及设备
US9602777B2 (en) 2008-04-25 2017-03-21 Roche Diagnostics Hematology, Inc. Systems and methods for analyzing body fluids
US8911669B2 (en) 2009-08-24 2014-12-16 Abbott Laboratories Method for flagging a sample
JP5366747B2 (ja) 2009-09-30 2013-12-11 シスメックス株式会社 分析装置
WO2011057744A1 (en) 2009-11-13 2011-05-19 Roche Diagnostics Gmbh Differential diagnosis of iron deficiencies based on hepcidin and mean hemoglobin content per reticulocyte
EP2519820B1 (en) 2009-12-31 2013-11-06 Abbott Point Of Care, Inc. Method and apparatus for determining mean cell volume of red blood cells
US8906308B2 (en) 2010-01-15 2014-12-09 Abbott Laboratories Method for determining volume and hemoglobin content of individual red blood cells
CN103430028A (zh) 2010-09-16 2013-12-04 通用医疗公司 供诊断用的红细胞动力学
US9064301B2 (en) 2011-04-14 2015-06-23 Abbott Point Of Care, Inc. Method and apparatus for compressing imaging data of whole blood sample analyses
BR112015000018B1 (pt) 2012-07-05 2021-11-09 Beckman Coulter, Inc Sistema e método automatizado para determinar uma condição de células brancas do sangue em uma amostra biológica
US20140324460A1 (en) 2012-09-26 2014-10-30 Health Diagnostic Laboratory, Inc. Method for determining and managing total cardiodiabetes risk
US9351670B2 (en) 2012-12-31 2016-05-31 Abbott Diabetes Care Inc. Glycemic risk determination based on variability of glucose levels
US20170108487A1 (en) 2014-06-05 2017-04-20 The General Hospital Corporation Predicting morbidity associated with red blood cell volume variance
AU2016229162A1 (en) 2015-03-06 2017-09-21 iRepertoire, Inc. Method for measuring a change in an individual's immunorepertoire
CA3022798A1 (en) 2016-05-04 2017-11-09 LabThroughput LLC System and method for distinguishing blood components

Also Published As

Publication number Publication date
JP2013541712A (ja) 2013-11-14
JP2016035467A (ja) 2016-03-17
EP3246709A2 (en) 2017-11-22
US20130236566A1 (en) 2013-09-12
IL253807A0 (en) 2017-09-28
JP5868982B2 (ja) 2016-02-24
EP2616818A2 (en) 2013-07-24
JP6185030B2 (ja) 2017-08-23
EP2616818A4 (en) 2014-04-02
US9938557B2 (en) 2018-04-10
IL225275A (en) 2017-09-28
IL253807B (en) 2018-07-31
EP3246709A3 (en) 2018-02-14
US20180187235A1 (en) 2018-07-05
IL225275A0 (en) 2013-06-27
WO2012037524A2 (en) 2012-03-22
EP2616818B1 (en) 2017-03-15
JP2018009992A (ja) 2018-01-18
WO2012037524A3 (en) 2012-06-21
US11319571B2 (en) 2022-05-03
CN103430028A (zh) 2013-12-04

Similar Documents

Publication Publication Date Title
JP6522059B2 (ja) 診断のための赤血球動力学
KR102094326B1 (ko) 위장암의 위험을 평가하는 방법 및 시스템
Song et al. Neutrophil to lymphocyte ratio is associated with in-hospital mortality in older adults admitted to the emergency department
EP3152322B1 (en) Estimating risk of death using the clearance volume of red blood cells as biomarker
JP7197795B2 (ja) 機械学習プログラム、機械学習方法および機械学習装置
US11885733B2 (en) White blood cell population dynamics
CN111524594A (zh) 目标人群血液系统恶性肿瘤筛查系统
JP2023548253A (ja) 医学検査値分析を実行するためのコンピュータ実装された方法、及び、装置
Horne et al. Extreme erythrocyte macrocytic and microcytic percentages are highly predictive of morbidity and mortality
US20220293210A1 (en) Single-cell modeling of clinical data to determine red blood cell regulation
US10274412B2 (en) Flow cytometry data segmentation result evaluation systems and methods
WO2022210606A1 (ja) 認知症の将来の発症リスクの評価方法
Yördan et al. Hybrid AI-Based Chronic Kidney Disease Risk Prediction
CN113496769B (zh) 基于病理组织的三联综合肿瘤分析系统和应用
IL280496A (en) Artificial intelligence for predicting laboratory test results
US20240127951A1 (en) Methods and systems for predicting baseline creatinine values
Guan et al. A functional perspective on the conditional covariance comparison problem in dementia analysis
Faria Predicting the personalized duration of phlebotomy treatments in patients with Hereditary Hemochromatosis: a model based approach
WO2024056468A1 (en) <u style="single">METHODS AND SYSTEMS FOR PREDICTING INTENSIVE CARE UNIT MORTALITY
JP2024525499A (ja) 機械学習を使用した全血算に基づく異常検出
ATAR et al. IS RED CELL DISTRIBUTION WIDTH CORRELATED WITH MORTALITY IN GERIATRIC PATIENTS UNDERGOING HIP FRACTURE SURGERY: A PROSPECTIVE OBSERVATIONAL STUDY.
Shareef et al. An alternative model for assessing mortality risk in Stevens Johnson syndrome/toxic epidermal necrolysis using a random forests classifier: A pilot study

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180405

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180515

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20180809

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181015

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: 20190402

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190423

R150 Certificate of patent or registration of utility model

Ref document number: 6522059

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250