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

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

Info

Publication number
JP5868982B2
JP5868982B2 JP2013529382A JP2013529382A JP5868982B2 JP 5868982 B2 JP5868982 B2 JP 5868982B2 JP 2013529382 A JP2013529382 A JP 2013529382A JP 2013529382 A JP2013529382 A JP 2013529382A JP 5868982 B2 JP5868982 B2 JP 5868982B2
Authority
JP
Japan
Prior art keywords
red blood
blood cell
volume
hemoglobin
fraction
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
JP2013529382A
Other languages
English (en)
Other versions
JP2013541712A (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 JP2013541712A publication Critical patent/JP2013541712A/ja
Application granted granted Critical
Publication of JP5868982B2 publication Critical patent/JP5868982B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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)
  • Urology & Nephrology (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medicinal Chemistry (AREA)
  • Biotechnology (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Microbiology (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • General Physics & Mathematics (AREA)
  • Food Science & Technology (AREA)
  • Cell Biology (AREA)
  • Organic Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Pharmacology & Pharmacy (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • General Chemical & Material Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Wood Science & Technology (AREA)
  • Zoology (AREA)
  • Diabetes (AREA)
  • General Engineering & Computer Science (AREA)
  • Biophysics (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号および第HL091331号の下で、政府の支援により行われた。政府は本発明における確実な権利を有する。
本発明は、貧血または貧血予備軍の患者を特定し、サラセミア形質貧血を鉄欠乏性貧血および他の原因の貧血と区別するための方法に関する。この方法を使用して、治療、例えば、鉄補給剤および/またはエリスロポエチン(EPO)等の赤血球生成促進剤による療法を最適化することもできる。また、例えば前記方法に使用するための装置およびシステムについても記載する。
循環するヒト赤血球(RBC)の数、大きさ、ヘモグロビン濃度、その他の特徴を制御するシステムについては、ほとんど解明されていない。RBCは骨髄から放出された後、未知の機構(Lew VL他、(1995)、Blood 86、334-341、Waugh RE他、(1992)、Blood 79、1351-1358)により、体積および総ヘモグロビン含有量が減少し、約120日後に未知の要因に反応して除去される。
統計物理学による理論および病院の臨床検査室からのデータ(d’Onofrio G他、(1995)Blood 85、818-823)の一部に基づき、本発明は、RBC成熟およびクリアランスのマスタ方程式モデルを使用する。モデルは、貧血患者を正確に特定し、サラセミア形質貧血を鉄欠乏性貧血と区別する。驚いたことに、このモデルは、貧血が臨床的に検出可能になる数週間前に、多くの貧血予備軍患者を特定する。
一態様では、被験者が鉄欠乏性貧血(IDA)を発症する危険を判定する方法は、被験者の赤血球を含むサンプル中で、母平均赤血球ヘモグロビン濃度(MCHC)を判定するステップと、赤血球体積およびヘモグロビン含有量をMCHC線上に投影することにより、各赤血球体積およびヘモグロビン含有量を指標に変換するステップとを含む。また方法は、変換された体積およびヘモグロビン含有量指標がMCHC線上の平均投影位置の閾値百分率よりも低いサンプル中の赤血球分画を判定して、サンプル分画を提供するステップと、サンプル分画を基準分画と比較するステップとを含む。基準分画よりも高いサンプル分画があれば、被験者がIDAを発症する危険があることを示す。
別の態様では、消化器(GI)疾患をスクリーニングするために被験者をスクリーニングまたは選択する方法は、被験者の赤血球を含むサンプル中で、母平均赤血球ヘモグロビン濃度(MCHC)を判定するステップと、赤血球体積およびヘモグロビン含有量をMCHC線上に投影することにより、各赤血球体積およびヘモグロビン含有量を指標に変換するステップとを含む。また方法は、変換された体積およびヘモグロビン含有量指標がMCHC線上の平均投影位置の閾値百分率よりも低いサンプル中の赤血球分画を判定するステップと、サンプル分画を基準分画と比較するステップとを含む。方法は、GI評価のために被験者を選択するか、またはサンプル分画が基準分画よりも低い場合にさらにスクリーニングするステップをさらに含む。
別の態様では、本開示は、小球性貧血を有する被験者、例えば、慢性疾患のない被験者の鉄欠乏性貧血(IDA)とサラセミア形質(TT)とを鑑別診断する方法を特徴とする。方法は、被験者の赤血球のサンプル中で、赤血球間のヘモグロビン含有量減少率(D)の変動の大きさの値を判定するステップと、Dを基準値と比較するステップとを含む。基準値よりも高いDがあれば、被験者がIDAであるかIDAになる可能性が高いことを示し、基準値よりも低いDがあれば、被験者がTTであるかTTになる可能性が高いことを示す。
別の態様では、健常者の血液ドーピングまたは赤血球生成促進剤の存在または使用を検出するための方法は、ヘモグロビン含有量減少率(D)の変動の大きさ、および被験者のサンプル中の赤血球体積減少率(D)の変動の大きさの一方または両方の値を判定するステップを含む。基準値よりも高いDおよび/またはDがあれば、血液ドーピングの存在もしくは使用または赤血球生成促進剤の使用を示す。
別の態様では、被験者のエリスロポエチン(EPO)または他の赤血球生成促進剤(ESA)の投与量を最適化するための方法は、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集合のパラメータに基づいて赤血球の定常状態分布を推定するステップとを含む。また、方法は、赤血球の推定定常状態分布と経験分布との相違の程度を表す目的関数を計算するステップを含む。
前述した方法、システム、コンピュータ読取可能な記憶装置の実施は、以下の特徴の組合せを含むことができる。
被験者の母平均体積
Figure 0005868982
の75%よりも大きいv、または78〜82%のvを維持するように、EPOの投与量を調節することができる。被験者の母平均体積
Figure 0005868982
の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の変動が、健常者の
Figure 0005868982
の変動よりもはるかに小さく(14)、これは要因がMCHC線上の位置と非常に相関していることを提示する。
貧血
貧血は、血流中に存在するヘモグロビンの量が正常よりも少ない状態である。貧血を、赤血球の大きさ(MCVまたは平均赤血球体積として測定される)に基づき、主に3つの群に細分することができる。正球性貧血は、赤血球が正常な大きさ(すなわち、80〜100fL)であるときに診断される。大球性貧血(通常、B12欠乏により生じる)は、赤血球の大きさが正常よりも大きいとき、すなわち、MCV>100fLであるときに診断される。小球性貧血は、赤血球が正常よりも小さいとき、すなわち、MCV<80fLのときである。
貧血予備軍の診断/貧血発症の予測
本明細書で示すモデルおよび方法は、クリアランス遅延の兆候を探すことにより、明らかな臨床的貧血が発症する前に、潜伏貧血または代償性貧血の患者を特定するための可能な方法を提供する。この可能性は、正常CBCを持っていたが、少なくとも30日および90日未満後に別の正常CBCまたは臨床的なIDAが生じた患者の独立集合において検査された。各患者のサンプルについて、全赤血球の(v、h)座標がMCHC線上に投影され、確率密度が、平均(P0.85)の85%未満でこの線に沿って統合された。図3A〜Cは、1患者のCBCが正常から潜伏IDAおよび最終的にIDAへ進む様子を示す。中央パネル(図3B)に示すCBCは、臨床的には目立たないものであるが、P0.85が異常であり、発症してはいないが51日の間に治療することになる貧血を予測している。図3Dは、健康なままである患者の20の正常CBCと30〜90日後にIDAを発症した患者の20の正常CBCについてのP0.85の値を示す。P0.85の値は、75%の感度および100%の特異性でIDAを予測した。実施例3「鉄欠乏性貧血の予測」を参照されたい。現在の標準治療は、すべてのCBCが「正常」であるため、この母集団で0%の感度を有する。このモデルに基づく予測は、異なる源および異なる時間点からの情報の複数の測定および種類の統合に依拠することの多い統計的回帰手法と対照的に、一時の単一CBC測定のみに依拠する(15)。
したがって、説明した方法を使用して、貧血予備軍を診断することができ、または発症の少なくとも30〜90日前、例えば、被験者がまだ正常なCBCを有する(例えば、被験者のHCT、MCV、RDW、MCH、MCHC、RBC、およびHGBレベルがすべて正常範囲内にある)ときに、被験者の貧血の発症を予測することができる。
正常CBC値は、分析を行う研究所によって変化し得るが、一般に、正常値は以下のとおりである。
Figure 0005868982
方法は、被験者のサンプル中の各赤血球について、赤血球体積およびHGB含有量(例えば、HGB質量または濃度)を判定するステップと、サンプル中の全赤血球にわたる体積およびHGB含有量値の分布を計算するステップと、例えば図5に示すような各赤血球体積およびヘモグロビンの変換を計算するステップと、サンプル中の全赤血球についてのこのような変換の分布を計算するステップと、例えば、図3A〜Dおよび8A〜Dに示すような、変換された体積およびHGB含有量の閾値レベルより変換が低い赤血球の百分率を判定するステップとを含む。一部の実施形態では、方法が、例えば、被験者の赤血球を含むサンプル中で、母平均赤血球ヘモグロビン濃度(MCHC)を判定するステップと、サンプル中の全赤血球にわたる体積およびHGB含有量値の分布を計算するステップと、例えば図5に示すような各赤血球体積およびヘモグロビンの変換を計算するステップと、サンプル中の全赤血球についてのこのような変換の分布を計算するステップと、変換が平均変換の閾値百分率よりも低いサンプル中の赤血球分画を判定して、サンプル分画を提供するステップと、サンプル分画を基準分画と比較するステップとを含み、基準分画よりも低いサンプル分画があれば、被験者が貧血、例えばIDAを発症する危険があることを示す。閾値百分率および基準レベルの値を、当技術で公知の方法を使用して選択することができ、この値は感度および特異性を最大化する値を表すことができる。一部の実施形態では、閾値百分率が母平均MCHCの70%、75%、80%、85%、90%、または95%である。一部の実施形態では、基準分画が約0.10、0.11、0.12、0.13、0.14、または0.15である。当業者は、最適な閾値百分率および基準分画を、公知の統計的な方法を使用して容易に特定することができるだろう。
鑑別診断 サラセミア対鉄欠乏性貧血
小球性貧血の鑑別診断は、鉄欠乏性貧血(IDA)、サラセミア形質(TT)慢性疾患に伴う貧血(ACD)、および他の原因による貧血を含む。
ACDは、通常、正常値の下限よりも20%未満までのヘマトクリットの減少と、正常または高いフェリチンと、低いまたは正常な総鉄結合能とを伴う。減少したRBCクリアランス閾値(v)が、貧血を補う適応生理学反応を表す場合には、ACDについてここで見られる正常vがおそらく予想され、貧血自体が適応生理学反応を表し得る(18)。慢性疾患を有する被験者に貧血があれば、通常、ACDの診断を提案する。通常ACDと関連する慢性疾患として、自己免疫疾患、例えば、クローン病、全身性エリテマトーデス、関節リウマチ、および潰瘍性大腸炎;癌、例えば、リンパ腫およびホジキン病;慢性腎臓病;肝臓病、例えば、肝硬変;慢性感染症、例えば、細菌性心内膜炎、骨髄炎(骨感染症)、HIV/エイズ、B型肝炎、またはC型肝炎がある。例えば、GardnerおよびBenz Jr.、「Anemia of chronic diseases.」 In: Hoffman他編、Hematology: Basic Principles and Practice.、第5版、Philadelphia、Pa: Elsevier Churchill Livingstone、2008、37章を参照されたい。
IDAは通常、ヘマトクリットの減少、低いMCV、および低いフェリチンを伴う。軽度のIDAは通常、正常値の下限よりも20%未満までのヘマトクリットの減少を伴う。IDAの患者は、正常ヘマトクリットを含む正常MCVの既存の証拠を示し得る。
TTは通常、高ヘモグロビンA2分画、または1もしくは複数のアルファグロビン遺伝子突然変異の存在と関連する。TTは、通常、正常値の下限よりも20%未満までのヘマトクリットの減少と、低いMCV、および正常なフェリチンを伴うことが多い。TTは、世界中で最も一般にスクリーニングされている疾患の1つであるが、既存の診断方法は非常に高価であるか、受け入れがたいほど診断精度が低く、偽陽性率が30%にのぼる(19、20)。本明細書に記載のモデルおよび方法は、IDAとTTと、場合によっては他の原因の貧血とを区別する、新しく、おそらくはより正確な方法を提供する。
図2は、Dが、小球性貧血の最も一般的な2つの原因であるIDAとTTとを区別することを示す。10の訓練サンプルを分析することにより、Dの例示的な閾値値0.0045が設定された。その後、軽度IDAまたはTTの診断を明確に行うことのできる50の独立した患者サンプルが分析され、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は母平均体積
Figure 0005868982
の約80%、または、一般的な90fLのMCVについて約72fLである。したがって、一部の実施形態では、方法が、EPOまたは鉄補給剤の投与量を調節して、被験者の
Figure 0005868982
の70%よりも高いv、被験者の
Figure 0005868982
の75%よりも高いv、被験者の
Figure 0005868982
の75〜85%よりも高いvを維持するステップを含む。一部の実施形態では、方法が、治療を調節して70%よりも高いvまたは70〜75%のvを維持するステップを含む。一部の実施形態では、vまたはDを時間とともに監視して、被験者に対するEPOまたは鉄補給剤の投与量を調節して、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−Dyn Sapphire)、およびSiemens(Deerfield、IL、United 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に接続する高速ユーザインターフェース708、ならびに低速バス714および記憶装置706に接続する低速ユーザインターフェース712を備える。構成要素702、704、706、708、710、712のそれぞれは、種々のバスを使用して相互接続され、共通のマザーボードに取り付けられるか、または必要に応じて他の方法で取り付けられ得る。プロセッサ702は、コンピューティング装置700内で実行するための命令を処理することができ、この命令には、メモリ704または記憶装置706に記憶された、例えば、高速ユーザインターフェース708に結合されたディスプレイ716を含む外部入出力装置にGUIのための図形情報を表示する旨の命令が含まれる。他の実施では、複数のメモリおよび複数のタイプのメモリとともに、複数のプロセッサおよび/または複数のバスを必要に応じて使用することができる。また、複数のコンピューティング装置700を接続することができ、各装置が必要な動作の一部を提供する(例えば、サーババンク、ブレードサーバ群、またはマルチプロセッサシステムとして)。
メモリ704は、コンピューティング装置700内で情報を記憶する。一実施では、メモリ704が1または複数の揮発性メモリユニットである。別の実施では、メモリ704が1または複数の不揮発性メモリユニットである。メモリ704は、例えば、磁気ディスクまたは光ディスクを含む別の形のコンピュータ読取可能な媒体であってもよい。
記憶装置706は、コンピューティング装置700に大容量記憶を提供することができる。一実施では、記憶装置706がコンピュータ読取可能な媒体であるか、またはコンピュータ読取可能な媒体を含むことができる。コンピュータ読取可能な媒体としては、例えば、フロッピーディスク装置、ハードディスク装置、光ディスク装置、テープ装置、フラッシュメモリもしくは他の同様の固体メモリ装置、またはストーレジエリアネットワーク内の装置または他の構成を含む多くの装置がある。コンピュータプログラム製品を、情報担体で明確に具体化することができる。コンピュータプログラム製品は、実行時に、例えば前述した方法を含む1または複数の方法を実施する命令を含むことができる。情報担体は、例えば、メモリ704、記憶装置706、プロセッサ702のメモリ等を含む、コンピュータまたは機械読取可能な媒体である。
高速コントローラ708は、コンピューティング装置700の帯域幅集約動作を管理し、低速コントローラ712は低帯域幅集約動作を管理する。このような機能の割当ては一例にすぎない。一実施では、高速コントローラ708がメモリ704、ディスプレイ716(例えば、図形プロセッサまたはアクセラレータを通して)、および種々の拡張カード(図示せず)を受けることのできる高速拡張ポート710に結合される。この実施では、低速コントローラ712が記憶装置706および低速拡張ポート714に結合される。種々の通信ポート(例えば、USB、Bluetooth(登録商標)、イーサネット、無線イーサネット)を含み得る低速拡張ポートを、例えば、キーボード、ポインティングデバイス、スキャナ、または例えばネットワークアダプタを通した、例えばスイッチもしくはルータを含むネットワークデバイスを含む1または複数の入出力装置に結合することができる。
コンピューティング装置700は、図示するようにいくつかの異なる形で実施可能である。例えば、コンピューティング装置700を、標準サーバ720として、またはこのようなサーバ群において複数回実施することができる。また、コンピューティング装置700をラックサーバシステム724の一部として実施することもできる。加えて、または代わりに、コンピューティング装置700を、例えば、ラップトップコンピュータ722を含むパーソナルコンピュータ内で実施することができる。一部の例では、コンピューティング装置700の構成要素を、例えば、装置750を含む携帯型装置(図示せず)の他の構成要素と組み合わせることができる。このような装置はそれぞれ、コンピューティング装置700、750の1または複数を含むことができ、システム全体を、互いに通信する複数のコンピューティング装置700、750から構成することができる。
コンピューティング装置750は、他の構成要素のうち、プロセッサ752、メモリ764、ならびに例えば、ディスプレイ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によって受信可能な、例えば、メモリ764、拡張メモリ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は母平均
Figure 0005868982
によってスケーリングされ、tは平均赤血球齢
Figure 0005868982
によってスケーリングされる。先のレポート(3〜5、7)からのデータに基づき、2つのパラメータが、決定論的コンポーネント、すなわちRBCがMCHC線に近くなるまでその効果が優位である高速変化(β)と、低速変化(α)とに挿入された。図1Bおよび式1に示すように、平均ゼロおよび分散が拡散テンソル2Dにより与えられた、ガウスまたは同様の分布を持つ不規則変数として、不規則変動をモデリングすることができる。
Figure 0005868982
一般的な逆問題および特別なヒト病態生理学(9)のように、この問題は、fの異なる関数形式が生体内動力学を再現するという意味で不良設定である。fの正確な関数形式は重要ではない。このモデルの挙動は高速および低速の決定論的動力学および不規則変動の定性的組合せに依拠した。表1A〜1Bは、モデルで使用可能な異なるタイプの関数形式を示す。
Figure 0005868982
Figure 0005868982
図9A〜Cは、fの形式A〜Cに対応する速度場を示し、図10A、Bは、クリアランス関数(d)形式AおよびBを示す。τは母集団内の赤血球の平均齢である。Δの定義を以下に述べる。
パラメータが、決定論的進化(f)およびクリアランス関数(d)の異なる関数形式について推定された。定性的結果はこれらの異なる関数形式について一致しており、これらの結果が生体内病態生理学の特徴を表し、データの過剰適合がないことを表すことを提示する。関数形式の詳細は、表1A〜Bに示され、形式A〜Cについての推定は、図2、11、および12の箱型図に示される。
本モデルでは、典型的な個々の赤血球についての体積およびヘモグロビン含有量の不規則変動および決定論的散逸または還元が、ポテンシャル(10)中のブラウン運動をモデリングするために通常使用されるランジュバン方程式により説明された。これにより、RBCの全循環母集団の動力学を、フォッカープランク方程式(10、11)により近似可能な時間依存性体積−ヘモグロビン確率分布(P(v,h,t))についてのマスタ方程式によって説明することができる。式2は、この体積−ヘモグロビン分布についての確率密度のドリフト(f)、拡散(D)、出生(b)、および死滅(d)を説明する。
Figure 0005868982
出生死滅プロセスは、常に母集団に加えられ母集団から除去されるRBCの要因となる。健常および軽度の病気の状態では、加えられる赤血球の総数は、除去される赤血球の総数に等しい。
Figure 0005868982
RBC除去の正確な要因および機構は完全には解明されていない(12)が、図1Aに示すような経験的測定により、MCHC線に沿った閾値(v)があり、大部分のRBCがMCHC線を超えて除去されたことが提示される。
経験的RBC分布の観察に基づき、RBCクリアランスの確率が、RBC体積およびヘモグロビン含有量の関数としてモデリングされる。体積−ヘモグロビン含有量平面における各RBCの位置がMCHC線上に投影され、クリアランス確率(d)が様々に定義され、例えば、この線状のこの投影点から閾値vまでの距離のシグモイド(図10A)またはステップ(図10B)関数として定義される。図1、5および表1Bを参照されたい。式3は、この関係を定量化する。
Figure 0005868982
CBC測定は人によって異なるが、これらの動力学プロセスが生体内定常状態Pに達することを示す健常者(13)については大きく変化しない。ここで
Figure 0005868982
所与のパラメータ集合について、病態生理学的範囲外の体積およびヘモグロビン含有量でゼロになる確率の境界条件と、経験的に測定された網状赤血球分布に等しい初期条件とを有する一次空間導関数
Figure 0005868982
および二次空間導関数
Figure 0005868982
についての有限差分近似を使用して、式2を数値的に解決することができる。体積−ヘモグロビン含有量平面は、一定のメッシュ幅で離散化され、各メッシュ赤血球中の確率密度に等しい変数のベクトル(P)として表される。ここで報告されるシミュレーション結果は、体積軸に沿った1.8fLおよびヘモグロビン含有量軸に沿った1.8pgのメッシュ幅により実行された。このメッシュ幅は、経験的体積およびヘモグロビン含有量測定の解析解と比較可能であった。より小さいメッシュ幅(1.2fLおよび1.2pg)を使用してすべての結果が確認された。空間導関数の風上有限差分近似を使用して、対流寄与(f)が数値的にモデリングされた。結果として得られる常微分方程式の線形システムが、MATLABode15s積分器を使用して積分され、定常状態(P)に達するまで反復された。
数値問題(P)についての定常状態分布は、解析的にも判定された。進化(J+L)項およびクリアランス(d)項の数値近似は線形演算子であり、出生プロセスの積分スケーリングは、平均齢の2倍の逆数に等しい定数である。
Figure 0005868982
線形演算子を反転させて、式4に示すように、平均赤血球齢により示される定常状態分布族を生じさせることができる。
Figure 0005868982
解析および有限差分手法により得られた定常状態分布間の差はごくわずかであった。
パラメータ(α、β、D、およびv)の適切な選択により、本明細書に記載のモデルは、健常者のRBC母集団の観察された分布を正確に再現する。
検証
モデルが健常者のRBC母集団の動力学を貧血患者のRBC母集団の動力学と区別することができるか否かを検査するために、原因疾患が異なる3種類の最もよく見られる貧血、すなわち炎症性疾患である慢性疾患に伴う貧血(ACD)、遺伝性疾患であるサラセミア形質(TT)、および栄養性疾患である鉄欠乏性貧血(IDA)を有する患者について、CBCおよび網状赤血球測定を行った(14)。RBC母集団特性が安定して見え、準定常状態仮定が適正な各貧血の軽度のケースが、複数の見かけの健常対照群として選択された。
血液サンプルおよびCBC結果は、Partners Healthcare Institutional Review Boardにより認可された研究プロトコルの下で、三次医療成人病院の臨床検査室から得られた。網状赤血球およびCBC測定が、Siemens Advia 2120自動血液分析器で、採取後6時間以内に行われた(21)。
IDAは、正常値の下限よりも20%未満までのヘマトクリットの軽度の減少、低いMCV、低いフェリチン、および正常ヘマトクリットを含む正常MCVの既存の証拠として定義される。急性疾患、急性出血、輸血後6か月以内、並行入院、慢性炎症性疾患、または異常ヘモグロビン症の患者は除外された。
TTは、高ヘモグロビンA2分画、またはアルファグロビン遺伝子突然変異の存在についてのヘテロ接合体として定義されるとともに、正常値の下限よりも20%未満までのヘマトクリットの減少と、低いMCV、および正常なフェリチンとして定義された。急性疾患、急性出血、輸血後6か月以内、並行入院、慢性炎症性疾患、またはさらなる異常ヘモグロビン症の患者は除外された。
ACDは、正常値の下限よりも20%未満までのヘマトクリットの減少と、正常または高いフェリチンと、低いまたは正常な総鉄結合能として定義された。急性疾患、急性出血、輸血後6か月以内、並行入院、慢性炎症性疾患、または異常ヘモグロビン症の患者は除外された。
各患者のサンプルに対して、患者について観察された定常状態を再現する最適なパラメータ集合(α、β、D、およびv)が特定された。シミュレーションした定常状態分布と測定したCBC分布との最小二乗適合を使用して、最良適合を特定した。同一の患者について繰返し検査を行うことができるが、適合パラメータの変動はCBC測定の解析変動によって説明されることがわかった。
勾配または非勾配最適化方法を使用して、パラメータ空間の最適近傍を各患者について特定した。患者の経験的に測定された網状赤血球分布および初期のランダムに選択されたパラメータ集合を、始点として使用した。式4を使用して、これにより生じる定常状態RBC分布を計算した。この計算された分布(P)を、その後、経験分布(PCBC)と比較した。パラメータ推定の質が、式5に示すように、離散化分布に対する正規化残差平方和に等しい目的関数を演算することにより定量化された。式中iおよびjは、離散化体積−ヘモグロビン平面の赤血球の指標を表す。
Figure 0005868982
一般に、目的関数は、2つの分布間の相違(または類似)の程度の測定値を提供する。目的関数の値が閾値条件を満たすか否かに基づいて、パラメータ値を調節することができる。例えば、目的関数の値が閾値よりも高くなるか、または低くなるまでパラメータ値を調節することができる。目的関数が相違の測定値(例えば平均二乗差、絶対差の和等)を表す場合、目的関数は、パラメータを調節することによって減少することが求められる。あるいは、目的関数が類似の測定値(例えば相関係数、相互情報等)を表す場合、目的関数は、パラメータを調節することによって増加することが求められる。この例では、この目的関数を減少させるために新しいパラメータ値が選択された。MATLABの勾配に基づく(lsqnonlin関数)、ならびに非勾配に基づく(fminsearchおよびpatternsearch関数)最適化アルゴリズムを使用して、最適パラメータを検索した。すべてのパラメータが非負に制限され、均一な初期パラメータ空間を定義して、1000日よりも長いか、または5日よりも短い平均赤血球齢を除外した。その後、ラテン超方格サンプリングを使用して、この空間から初期パラメータを選択した。1つの最適化制限が課され、αを平均赤血球齢が5日よりも長くなるのに十分な小ささに限定し、平均赤血球齢が1000日よりも短くなるのに十分な大きさに限定する。図6は、モデルが、シミュレーションされた定常状態確率分布と測定された定常状態確率分布とを比較することにより、この健常者についてのPCBCを正確に再現することを示す。MCHC線に沿って投影されると、この患者についての経験分布はモード近くでわずかに高い密度を有し、モードのいずれかの側までの領域でわずかに低い密度を有していた。
図7は、1人の患者に対する200を超える最適化から得られたパラメータの最小の局所解を示す。一部の結果は、軸の範囲よりも高い局所解を有した。すべてのシミュレーション間の最良適合は、全パラメータについての値の小近傍を形成し、パラメータ推定プロセスが、この患者の血液サンプルについての明確に定義された最適近傍に達したことを示す。
健常者および貧血患者についての適合パラメータが図2に示され、中央値適合パラメータが表2に挙げられる。
Figure 0005868982
健常者から得た最良適合パラメータと貧血患者から得た最良適合パラメータとの間には明白な差があり、異なる貧血疾患は異なる特性パラメータ集合を有した。例えば、健常者およびACD患者は、高いβおよびβと低いαとを有し、すなわち、低速相中よりも高速相中に体積およびヘモグロビンが比較的多く失われた。これに対して、TTおよびIDA患者は、高速相中よりも低速相中に体積およびヘモグロビンが比較的多く失われた。ACD患者は、健常者に対してわずかに高いDとわずかに低いDとを示したが、TTはDの実質的な減少に伴うより大きなDの増加に関連した。IDA患者は健常者と同様のDを有し、大部分の患者が健常者の10倍を超える、Dの劇的な増加を示した。健常者およびACD患者の正規化臨界体積vは、
Figure 0005868982
の約80%、または約72fLである。TTまたはIDA患者の大半は、通常、
Figure 0005868982
および
Figure 0005868982
が減少した。図2は、
Figure 0005868982
および
Figure 0005868982
の絶対的な減少に加えて、これらの患者のvがさらに減少下し、異なる患者によってはるかに大きなばらつきがあったことを示す。
鉄欠乏性貧血の予測
本実施例は、平均よりも起点に近いMCHC線に沿って投影する赤血球の拡大母集団に基づいて、代償性または潜伏IDAを、場合によっては、現在可能であるよりも少なくとも90日早く予測することができるという仮説を検査した。投影動作を図5に示す。まず、体積(v)およびヘモグロビン(h)を有する各赤血球のMCHC線に沿った投影位置(u)が判定された。
Figure 0005868982
MCHC線に沿った閾値は、平均投影位置の比率(φ)
Figure 0005868982
として定義された。次に、起点とこの閾値との間に位置する投影された赤血球分画が以下のように計算された。式中、fMCHCは、MCHC線上の位置の関数としての投影された赤血球の確率密度である。
Figure 0005868982
図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.009〜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は57%である。このような他の判別関数のさらなる詳細については参考文献(20)を参照されたい。
(参考文献)
Figure 0005868982

Figure 0005868982
他の実施形態
本発明をその詳細な説明とともに説明したが、前述した説明は例示的なものであり、添付の特許請求の範囲によって定義される本発明の範囲を限定するものではないことを理解されたい。その他の態様、利点、変更が、以下の特許請求の範囲に含まれる。
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 (11)

  1. 被験者が鉄欠乏性貧血(IDA)を発症する危険を判定する方法であって、
    前記被験者の赤血球を含むサンプル中で、母平均赤血球ヘモグロビン濃度(MCHC)を判定するステップと、
    赤血球体積およびヘモグロビン含有量をMCHC線上に投影することにより、各赤血球体積およびヘモグロビン含有量を指標に変換するステップと、
    変換された体積およびヘモグロビン含有量指標が前記MCHC線上の平均投影位置の閾値百分率よりも低い前記サンプル中の赤血球分画を判定して、サンプル分画を提供するステップと、
    前記サンプル分画を基準分画と比較するステップとを含み、
    前記基準分画よりも高いサンプル分画があれば、前記被験者がIDAを発症する危険があることを示す方法。
  2. 消化器(GI)疾患をスクリーニングするために被験者をスクリーニングまたは選択する方法であって、
    前記被験者の赤血球を含むサンプル中で、母平均赤血球ヘモグロビン濃度(MCHC)を判定するステップと、
    赤血球体積およびヘモグロビン含有量をMCHC線上に投影することにより、各赤血球体積およびヘモグロビン含有量を指標に変換するステップと、
    変換された体積およびヘモグロビン含有量指標が前記MCHC線上の平均投影位置の閾値百分率よりも低い前記サンプル中の赤血球分画を判定するステップと、
    前記サンプル分画を基準分画と比較するステップと、
    前記サンプル分画が前記基準分画よりも低い場合、GI評価のために前記被験者を選択するか、またはさらにスクリーニングするステップとを含む方法。
  3. 前記閾値百分率が、前記MCHC線上の前記母平均投影の70%、75%、80%、85%、90%、または95%である、請求項1または2に記載の方法。
  4. 前記基準分画が、約0.10、0.11、0.12、0.13、0.14、または0.15である、請求項1または2に記載の方法。
  5. 前記サンプルは、前記被験者の全血を含む、請求項1〜4のいずれか一項に記載の方法。
  6. 命令を記憶するためのメモリと、
    動作を行う旨の前記記憶された命令を実行可能な1または複数のプロセッサまたは処理装置と
    を有するコンピューティング装置を備え、前記動作が、
    サンプル中の各赤血球の体積およびヘモグロビン濃度または含有量を表すデータを受信するステップと、平均赤血球ヘモグロビン濃度(MCHC)を表す線上に赤血球体積およびヘモグロビン含有量を投影することにより、各赤血球体積およびヘモグロビン含有量を指標に変換するステップと、
    変換された体積およびヘモグロビン含有量指標が平均投影の閾値百分率よりも低い前記サンプル中の赤血球分画を計算するステップと、
    前記分画を表す出力を提供するステップとを含む、システム。
  7. 前記サンプル中の平均赤血球ヘモグロビン濃度(MCHC)を表すデータを受信する旨の命令をさらに含む、請求項に記載のシステム。
  8. 赤血球(RBC)の赤血球体積(CV)、平均赤血球体積(MCV)、血球ヘモグロビン濃度(CHC)、平均赤血球ヘモグロビン濃度(MCHC)、および平均赤血球ヘモグロビン含有量(MCH)の1または複数、例えば、すべてと、これらの母集団統計とを検出するように構成された検出モジュールをさらに備える、請求項に記載のシステム。
  9. 1または複数のプロセッサまたは処理装置により実行可能な動作を行う旨の命令を格納するように構成された、1または複数の機械読取可能な記憶装置であって、前記動作が、
    各赤血球の体積およびヘモグロビン濃度を表すデータを受信するステップと、平均赤血球ヘモグロビン濃度(MCHC)を表す線上に赤血球体積およびヘモグロビン含有量を投影することにより、各赤血球体積およびヘモグロビン含有量を指標に変換するステップと、
    変換された体積およびヘモグロビン含有量指標が平均投影の閾値百分率よりも低い前記サンプル中の赤血球分画を計算するステップと、
    前記分画を表す出力を提供するステップとを含む記憶装置。
  10. 前記サンプル中の平均赤血球ヘモグロビン濃度(MCHC)を表すデータを受信する旨の命令をさらに含む、請求項に記載の機械読取可能な記憶装置。
  11. (i)ヘモグロビン含有量減少率(D)の変動の大きさ、
    (ii)赤血球体積減少率(D)の変動の大きさ、
    (iii)正規化臨界体積またはクリアランス閾値(v)、
    (iv)低速相体積およびヘモグロビン含有量の平均減少率(α)、
    (v)高速相体積の平均減少率(β)、ならびに
    (vi)高速相ヘモグロビン含有量の平均減少率(β
    のうちの1または複数を計算する旨の命令をさらに含む、請求項6〜8のいずれか一項に記載のシステム、または請求項9または10に記載の機械読取可能な記憶装置。
JP2013529382A 2010-09-16 2011-09-16 診断のための赤血球動力学 Active JP5868982B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US38335710P 2010-09-16 2010-09-16
US61/383,357 2010-09-16
PCT/US2011/052038 WO2012037524A2 (en) 2010-09-16 2011-09-16 Red blood cell dynamics for diagnosis

Related Child 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
JP2013541712A JP2013541712A (ja) 2013-11-14
JP5868982B2 true JP5868982B2 (ja) 2016-02-24

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 After (2)

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

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
WO2012037524A2 (en) 2010-09-16 2012-03-22 The General Hospital Corporation Red blood cell dynamics for diagnosis
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
US10955423B2 (en) 2015-12-15 2021-03-23 The General Hospital Corporation Methods of estimating blood glucose and related systems
WO2017177192A1 (en) 2016-04-07 2017-10-12 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
WO2022102734A1 (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 日本光電工業株式会社 血球計数装置
EP0626067A4 (en) 1992-02-07 1997-06-25 Abbott Lab METHOD FOR ACCURATELY COUNTING AND SENSITIVELY CLASSIFYING HETEROGENEOUS CELL POPULATIONS UNDER 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
AU775260B2 (en) 1999-03-31 2004-07-29 Bayer Corporation Single channel, single dilution detection method
AU5901501A (en) 1999-10-19 2001-10-23 Nicholas R. Abu-Absi 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
JP2005503559A (ja) 2001-09-14 2005-02-03 エフ.ホフマン−ラ ロシュ アーゲー 鉄代謝異常の判別診断
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
US7361510B2 (en) * 2006-01-20 2008-04-22 Beckman Coulter, Inc. Methods of detection of iron deficiency
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
WO2009117664A2 (en) 2008-03-21 2009-09-24 Abbott Point Of Care, Inc. Method and apparatus for determining red blood cell indices of a blood sample utilizing the intrinsic pigmentation of hemoglobin contained within the red blood cells
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
WO2012037524A2 (en) 2010-09-16 2012-03-22 The General Hospital Corporation Red blood cell dynamics for diagnosis
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
EP2870458B1 (en) 2012-07-05 2020-09-09 Beckman Coulter, Inc. Method and apparatus for determining white blood cell counts
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
US10529442B2 (en) 2015-03-06 2020-01-07 iRepertoire, Inc. Method for measuring a change in an individual's immunorepertoire
WO2017192739A1 (en) 2016-05-04 2017-11-09 LabThroughput LLC System and method for distinguishing blood components

Also Published As

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

Similar Documents

Publication Publication Date Title
JP6185030B2 (ja) 診断のための赤血球動力学
Burke Genetic tests: clinical validity and clinical utility
KR102094326B1 (ko) 위장암의 위험을 평가하는 방법 및 시스템
Pérez-Persona et al. New criteria to identify risk of progression in monoclonal gammopathy of uncertain significance and smoldering multiple myeloma based on multiparameter flow cytometry analysis of bone marrow plasma cells
JP2019522774A (ja) 試料に対する光学測定の実施
Pedersen et al. cyCombine allows for robust integration of single-cell cytometry datasets within and across technologies
EP3152322B1 (en) Estimating risk of death using the clearance volume of red blood cells as biomarker
Schwarz et al. Diagnostic value of the neutrophil myeloperoxidase index in horses with systemic inflammation
Spielmann et al. Accuracy and precision of hemoglobin point‐of‐care testing during major pediatric surgery
Anoke et al. Approaches to treatment effect heterogeneity in the presence of confounding
Sehgal et al. Reference range evaluation of complete blood count parameters with emphasis on newer research parameters on the complete blood count analyzer Sysmex XE-2100
Ghosh Semiparametric inferences for association with semi‐competing risks data
WO2017177192A1 (en) White blood cell population dynamics
Goag et al. A simplified mortality score using delta neutrophil index and the thrombotic microangiopathy score for prognostication in critically ill patients
Sekhar et al. Optimizing hemoglobin thresholds for detection of iron deficiency among reproductive-age women in the United States
US20220293210A1 (en) Single-cell modeling of clinical data to determine red blood cell regulation
Lewerin et al. Holotranscobalamin is not influenced by decreased renal function in elderly men: the MrOS Sweden study
Yu et al. Leukemia can be Effectively Early Predicted in Routine Physical Examination with the Assistance of Machine Learning Models
Cook Comment: measures to summarize and compare the predictive capacity of markers
Chauhan et al. Z score analysis: A novel approach to interpretation of an erythrogram
US20240127951A1 (en) Methods and systems for predicting baseline creatinine values
CN113496769B (zh) 基于病理组织的三联综合肿瘤分析系统和应用
Shen et al. The use of statistics in heart rhythm research: a review
Zekavat et al. A computational model of 1, 5‐AG dynamics during pregnancy
Thuluvath et al. A model to predict inhospital mortality in patients with cirrhosis, ascites and hyponatremia

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140916

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150415

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150526

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20150826

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150928

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160106

R150 Certificate of patent or registration of utility model

Ref document number: 5868982

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

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