CN101076806A - 败血症的早期检测 - Google Patents
败血症的早期检测 Download PDFInfo
- Publication number
- CN101076806A CN101076806A CNA2005800424703A CN200580042470A CN101076806A CN 101076806 A CN101076806 A CN 101076806A CN A2005800424703 A CNA2005800424703 A CN A2005800424703A CN 200580042470 A CN200580042470 A CN 200580042470A CN 101076806 A CN101076806 A CN 101076806A
- Authority
- CN
- China
- Prior art keywords
- septicemia
- analysis
- group
- data
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/68—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids
- G01N33/6863—Cytokines, i.e. immune system proteins modifying a biological response such as cell growth proliferation or differentiation, e.g. TNF, CNF, GM-CSF, lymphotoxin, MIF or their receptors
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/68—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids
- G01N33/6872—Intracellular protein regulatory factors and their receptors, e.g. including ion channels
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
- G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
- G16B40/10—Signal processing, e.g. from mass spectrometry [MS] or from PCR
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/158—Expression markers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2333/00—Assays involving biological materials from specific organisms or of a specific nature
- G01N2333/435—Assays involving biological materials from specific organisms or of a specific nature from animals; from humans
- G01N2333/705—Assays involving receptors, cell surface antigens or cell surface determinants
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2800/00—Detection or diagnosis of diseases
- G01N2800/26—Infectious diseases, e.g. generalised sepsis
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Immunology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- General Health & Medical Sciences (AREA)
- Biotechnology (AREA)
- Genetics & Genomics (AREA)
- Analytical Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Urology & Nephrology (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Hematology (AREA)
- Biomedical Technology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Cell Biology (AREA)
- Pathology (AREA)
- Microbiology (AREA)
- Theoretical Computer Science (AREA)
- Organic Chemistry (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biochemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- General Physics & Mathematics (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Artificial Intelligence (AREA)
- Databases & Information Systems (AREA)
- Public Health (AREA)
- Software Systems (AREA)
- Epidemiology (AREA)
- Bioethics (AREA)
- Evolutionary Computation (AREA)
- Signal Processing (AREA)
Abstract
本发明涉及一种检测感染早期征象尤其是鉴别最有可能发生败血症的个体的系统和方法。细胞因子和/或细胞活化标记物的特定组合的表达水平的测量,可选地结合采用预测算法,可以得到高的预测准确度。该方法适用于民用和军事环境。
Description
背景
尽管严重感染和败血症的诊断、治疗以及支持疗法有很大的改进,但是它们仍然是很重要的死因,经常导致在急性发病中生存下来的人们产生长期的体弱或致残。尽管在其它健康成人中突然、极其严重的感染相对较为稀少,但是在免疫缺陷的个体、重症监护的重症患者、烧伤病人和年幼的儿童中,其风险增加。在部分病例中,表面上能够治疗的感染导致败血症的发生;表现为进行性循环衰竭的对感染不规律、不适当的应答导致肾衰竭和呼吸衰竭、凝血异常、严重并且无应答的低血压,并使得大约30%的病例死亡。在北美人群中败血症的发病率每年为人群的大约0.3%(大约750000例),在老年人中死亡率增加至40%,在最严重的病例即脓毒性休克中,死亡率增加至50%(Angus等,2001,Crit Care Med 29:1303-1310)。
在感染了感染性微生物后,机体产生经典的炎症反应,首先为激活先天的非特异性免疫应答,然后为特异性的、获得性免疫应答。在细菌感染的情况中,在发展成为特异性的、抗原引发的免疫应答的很早之前,通过检测到细菌毒素,在发展到特异性的、抗原启动的免疫应答的很长一段时间以前,菌血症导致快速的发热(30-90分钟以内),并释放炎症细胞因子,例如白介素1(IL-1)和肿瘤坏死因子-α(TNF-α)。
在感染导致的革兰氏阴性菌血症中,例如伤寒、鼠疫、兔热病和布鲁氏菌病,或者在革兰氏阴性消化道微生物例如大肠杆菌、克氏杆菌、变形杆菌或假单胞菌导致的腹膜炎中,大部分是对来自细菌细胞壁的脂多糖(LPS)和其它组分的响应。循环LPS,尤其是其成分脂质A,激发了广泛的系统反应。其有可能与肝细胞的Kupffer细胞接触,首先导致IL-1的释放,并导致发热。循环系统中单核细胞和巨噬细胞的激活导致细胞因子的释放,例如IL-6、IL-12、IL-15、IL-18、TNF-α、巨噬细胞移动抑制因子(MIF)、以及细胞因子类分子,例如高迁移群B1(HMGB1),其依次激活嗜中性粒细胞、淋巴细胞和血管内皮上调细胞粘附分子,并诱导前列腺素、氮氧化物合成酶、和急性期蛋白。血小板活化因子(PAF)、前列腺素、白细胞三烯和凝栓质降解物的释放激活了血管内皮,调节血管紧张度,并激活外源性凝血级联。这些响应的不当调节导致并发败血症和脓毒性休克,即导致低血压的外周血管舒张,以及导致血栓形成和血管内凝血的异常凝血和血纤维蛋白溶解(Cohen,2002,Nature 420:885-891)。
LPS主要通过结合至血清LPS-结合蛋白(LBP)以及单核细胞和巨噬细胞表达的CD14而作用于细胞。一旦结合LPS和LBP的复合物,CD14与co受体、Toll类受体(TLR-4)和另一组分MD-2一起作用,形成信号复合物,并启动巨噬细胞的激活和细胞因子的释放(Palsson-McDermott&O’Neill,2004,Immunology 113:153-162)。Toll类受体家族是一组细胞表面受体,涉及一类作为先天免疫系统激活剂的细菌和真菌配体,包括革兰氏阳性细胞壁结构、鞭毛蛋白和表征细菌DNA的CpG重复单元。
在感染革兰氏阳性病原的情况下,脓毒性休克伴有外毒素的生成。例如,中毒性休克综合症,通常影响其它健康个体的脓毒性休克的一种特殊急性形式,是由于感染了金黄色葡萄球菌的特定株,其产生了被称为中毒性休克综合症毒素-1(TSST-1)的外毒素。一个类似的综合症由特定的化脓性链球菌A菌株的入侵感染引起的,其通常伴有链球菌化脓性肠毒素A(SPE-A)。某些革兰氏阳性外毒素(包括TSST-1)被认为是主要由于其超抗原性能而起作用。通过将抗原呈递细胞上的MHC II型分子与特定类型的T细胞受体交联,超抗原能够非特异性地刺激T淋巴细胞。通常,T细胞受体(TCR)-主要组织相容性复合物(MHC)的交互作用是高度特异的,只有携带特异性识别MHC呈递的短抗原衍生肽的TCR的T细胞才能够结合并被激活,以确保抗原特异性T细胞响应。超抗原规避了这种导致大规模的、不适当的T细胞激活的机制。然而,SPE-A并不是有效的超抗原,其应当涉及某些其它机制。
应当注意的是,临床败血症还可能来自某些病毒(例如,委内端拉马脑脊髓炎病毒,VEEV)和真菌的感染,在此类情况中,可能涉及其它的机制。
尽可能早地检测潜在严重感染的能力,尤其是在易感个体中预测败血症的发生显然是非常有利的。多年来,在尝试建立明确的确定临床疾病单元的标准方面投入了大量的努力,例如休克、败血症、脓毒性休克、中毒性休克和系统炎症反应综合症(SIRS)。类似地,基于所测量的一系列临床、化学、生化、免疫学和细胞学参数,以及大量的评分系统、所提出的不同的预测成功性和准确性,在设计可靠的预测模型方面进行了许多尝试。
根据美国胸内科医生学会(ACCP)和重症监护医学学会(SCCM)的1991年会议达成的共识,当患者具有以下一个以上症状时被认为具有SIRS:体温高于38℃或低于36℃,心率高于90/min,涉及呼吸速率高于20/min或PaCO2低于32mmHg的换气过度,白细胞计数高于12000个/ul或少于4000个/ul(Bone等,1992,Crit Care Med 20:864-874)。
“败血症”被定义为感染引起的SIRS。被接受的是,可以在没有感染的情况下例如烧伤、胰腺炎和其它疾病状态中发生SIRS。“感染”被定义为由致病或潜在致病的微生物侵入通常无菌的组织、液体或体腔而引发的病理过程。
“严重败血症”被定义为并发器官衰竭的败血症,其由Marshall等(1995,Crit Care Med 23:1638-1652)或Sequential Organ FailureAssessment(器官顺序衰竭评价,SOFA)得分(Ferreira等,2002,JAMA286:1754-1758)来确定。
“脓毒性休克”是指(成人中)败血症加上急性循环衰竭状态,其特征为其它原因不能解释的持续性动脉低血压。
为了评价重症监护患者中败血症的严重性,并进行合理的治疗规划,对败血症研究出了或者从更加通用的模型中改编得到了大量的临床严重程度模型。首先被广泛接受的系统为Acute Physiology andChronic Health Evaluation得分(急性生理和慢性健康评价得分,APACHE,以及其修定版APACHEII和III)(Knaus等,1985,CritCare Med 13:818-829;Knaus等,1991,Chest 100:1619-1636),和被广泛采用的通用预测模型Mortality Prediciton Model(死亡率预测模型,MPM)(Lemeshow等,1993,JAMA 270:2957-2963),以及Simplified Acute Physiology(SAPS)得分(简化的急性生理得分)(Le Gall等,1984,Crit Care Med 12:975-977)。对于更加严重的状况,包括败血症,研究出了更加专业的模型,例如Multiple OrganDysfunction Score(多器官障碍得分,MODS)(Marshall等,1995,Crit Care Med 23:1638-1652),Sequential Organ Failure Assessment(SOFA)得分(器官顺序衰竭评价得分)(Ferreira等,2002,JAMA286:1754-1758)和Logistical Organ Dysfunction Score(LODS)(Logistical器官障碍得分,Le Gall等,1996,JAMA276:802-810)。最近,提出了一个专门的模型PIRO(Levy等,2003,Intensive Care Med29:530-538)。所有这些模型采用了广泛的普通和专业临床手段的组合,试图得到反映患者状况严重程度和可能结果的有用评分。
除了上述的标准预测模型之外,对败血症和大量特异的血清标记物之间的关联进行了深入的研究,其目的是研究出特异性的诊断和预测检测,其中有:
C反应蛋白(CRP)是来源于肝的一种血清急性期蛋白,其是公知的非特异性炎症标记物。最近(Toh等,2003,Intensive Care Med29:55-61),被称为脂蛋白复合C-反应蛋白(LCCRP)的依赖于钙的CRP和极低密度脂蛋白(VLDL)的复合物显示在败血症中涉及影响凝血机制。尤其是,称为活化部分凝血活酶时间的常见化验在败血症病例中显现出特定的模式,其已经被提出作为快速诊断检测的基础。
TNF-α和IL-1是长久以来都知晓的在败血症中升高的原型急性炎症细胞因子(Damas等,1989,Critical Care Med 17:975-978),已经报道其是成人呼吸窘迫综合症的器官衰竭中的有用预测物,上述综合症是败血症的一个严重并发症(Meduni等,1995,Chest 107:1062-1073)。
已经提出将活性补体产物C3(C3a)和IL-6作为对宿主对微生物入侵响应的有用指示物,优于发热和白细胞计数(Groeneveld等,2001,Clin Diagn Lab Immunol 8:1189-1195)。在同一研究中,发现分泌型磷脂酶A2并不是很可靠的标记物。
前降钙素是降钙素的前肽前体,知道的是,其血清浓度响应于LPS而升高,并与IL-6和TNF-α水平相关。对其作为败血症预测物已经进行了评估(AI-Nawas等,1996,Eur J Med Res 1:331-333)。采用0.1ng/ml为阈值,其正确地鉴定出39%的败血症患者。然而,其他报道提示,其没有采用系列CRP测定可靠(Neely等,2004,J Burn CareRehab 25:76-80),但是其优于IL-6或IL-8(Harbarth等,Am J RespCrit Care Med 164:396-402)。
白细胞激活标记物的嗜中性粒细胞的表面表达的改变(例如CD11b,CD31,CD35,L-选择蛋白,CD16)已经作为SIRS的标记物,并被发现与IL-6和接下来的器官衰竭的发生相关(Rosenbloom等,1995,JAMA 274:58-65)。类似地,还检查了血小板表面抗原的表达例如CD63,CD62P,CD36和CD31,但是没有建立可靠的预测模型。
最后,单核细胞HLA-DR表达的下调显示为败血症预后较差的预测物,其可能指示单核细胞的去活化,损害了TNF-α的生成。用IFN-γ进行治疗显示在所述病例中是有利的(Docke等,1997,Nature Med 3:678-681)。
然而,尽管这些标记物中的多个与败血症相关,某些指示了状况的严重程度,但是仍然没有一个标记物或标记物的组合显示作为可靠的诊断检测,更不用说作为败血症发展的预测物。2001年国际败血症定义会议(2001 International Sepsis Definition Conference)得出结论:“采用生物标记物以诊断败血症是不成熟的”(Levy等,2003,Intensive Care Med 29:530-538)。
从随着时间变化的复杂的变量组中提取可靠的诊断模式和可靠的预测指示物要求非常复杂的分析方法,所述变量包括传统的临床观察资料,临床化学,生化、免疫和细胞学数据。对于医学诊断应用采用专家系统和人工智能,包括神经网络已经发展了一定时间(Place等,1995,Clinical Biochemistry 28:373-389;Lisboa,2002,Neural Networks15:11-39)。通过采用多元logistic回归和采用APACHE得分以及上述的1991 ACCP/SCCM SIRS标准的神经网络模型(Bone等,1992,CritCare Med 20:864-874),已经研究出了专门的系统以试图预测败血症患者的生存(Flanagan等,1996,Clinical Performance & Quality HealthCare4:96-103)。所述研究提示,尽管两种方法都可以给出良好的预测结果,但是神经网络系统对预选择的阈值不是很敏感(Rosenberg对大量研究的结果进行了综述,2002,Curr Opin Crit Care 8:321-330)。Brause等(2004,Journal fur Anasthesie und Intesivbehandlung11:40-43)提供了用于预测败血症的神经网络模型的例子。该模型(MEDAN)分析了一系列的标准临床手段,并将其结果与采用APACHE II、SOFA、SAPS II和MODS模型的结果比较。该研究得出以下结论,在可利用的标记物中,信息量最丰富的是收缩压和舒张压以及血小板计数。
神经网络为能够确定复杂数据系统的模式的非线性函数。其通过采用大量的数学函数来实现,从而使所述网络能够在噪声数据集中确定结构。其原因是,来自一个系统的数据基于数据内变量之间的关系可以产生模式。如果在已知的“训练”期,神经网络看见足够多的所述数据点的例子,则其能够“学习”该结构,然而在以后的数据点或测试数据中确定这些模式。这样,神经网络通过对已经看到的数据内存在的模式进行建模,能够预测或对未来例子进行分类。然后通过所述网络正确对测试数据进行预测或分类的能力,对该网络的性能进行评估。所述神经网络的平行处理能力依赖于其处理元素的结构,所述结构根据生物神经元的模式安排成相互作用。通过连接权重调节一个或多个输入,以改变所述处理元素内的刺激水平。所述处理元素的输出与其激活水平有关,该输出可以是非线性或不连续的。因此,神经网络的训练包括根据元素之间的转移函数、互联结构的细节以及该系统遵循的学习规则调整互连权重(Place等,1995,Clinical Biochemistry28:373-389)。所述系统已经被应用到大量的临床情况中,包括创伤患者的健康结局模型(Marble & Healy(1999)Art Intell Med 15:299-307)。
Warner等(1996,Ann Clin Lab Sci 26:471-479)描述了一预测败血症结局的多参数模型,其采用测量“脓毒性休克因子”(患者在收治时是否具有脓毒性休克症状,其看上去比较简单)、IL-6、可溶性IL-6受体(通过酶联免疫吸附试验测量)以及APACHE II得分作为多层、前馈神经网络模型的四输入算法的组分。然而,该系统对于仍然没有临床表现的个体来说是非预测性的,而且有疑义的是,通过细胞因子例如IL-6的时间血清水平的升高,如果不是结果的话,其诊断在临床上是非常明显的。
Dybowski等(1996,Lancet 347:1146-1150)采用分类和回归树(CART,Classification and Regression Tree)从157个可能的败血症预测规则中选择输入,然而采用运行遗传算法的神经网络以选择预测标记物的最佳组合。其包括许多常规的临床值和代表指示物,而不是血清或细胞表面标记物。然而,存在的问题是已经明确诊断败血症而且已经非常严重的患者的预后。
该遗传算法方法的进一步改进涉及采用人工免疫系统(ArtificialImmune Systems),所述系统的一个版本是人工免疫识别系统(AIRS,Artificial Immune Recognition System)(Timmis等,An overview ofartificial immune system.In:Paton,Bolouri,Holcombe,Parish和Tateson(eds.)“Commputation in cells and tissues:Perspectives and Toolsfor Thought”,Natural Computation Series,51-86页,Springer,2004;Timmis(L.N.De Castro和J.Timmis,Artificial Immune Systems:A newcomputational intelligence approach.Springer-Verlag,2002),其是受到应用到人工智能的生物免疫系统的克隆选择和亲和成熟过程的激发而改编的系统。
从免疫学上说,AIRS受到免疫系统的克隆选择理论的激发(F.Burnett.The clonal selection theory of acquired immunity.Cambrideg University Press,1959)。克隆选择理论试图解释如何通过匹配、克隆、突变和选择的过程制造出能够确定感染物的抗体。AIRS利用该过程,通过匹配、克隆和突变的过程,涉及一组能够用作未见数据项目的分类器的存储检测器。与其他免疫激发的方法不同,例如负选择,AIRS专门设计用于分类,尤其是一次就完成的监督学习。
美国专利申请2002/0052557描述了一种基于患者心率变化预测大量灾难性疾病发生的方法。再次,在对所述数据进行建模和分析的可能方法中包括神经网络。
国际专利申请WO00/52472描述了一用于小儿的快速检测方法,其基于血清或嗜中性粒细胞表面的CD11b或“CD11b复合物”(Mac-1,CR3)水平。该方法仅仅采用一种标记物,该标记物是一种有争论的、响应于炎症的嗜中性粒细胞活化的公知标记物。
分析所述复杂数据集的另一种方法是采用复杂的统计分析方法,例如logistic回归,所述复杂数据集的数据通常是定性和离散的,而不是定量和连续的。尽管采用定性的二元因变量的Logistic回归就选择显著变量而言并不能够充分地辨别,但是可以采用多变量方法。多元logisitc回归模型和神经网络的输出均是连续的变量,但是神经网络模型计算的似然性通常落在一个端点或另一个端点,很少有值位于中间范围。在临床中,这通常是有帮助的,能够给出更加明确的结果(Flanagan等,1996,Clinical Performance & Quality Health Care4:96-103)。
检测感染和/或败血症的最早征兆的能力就尽可能早地进行治疗而言具有显然的优势。状况严重程度和如果不治疗的可能结局的指示告知了关于治疗选择方案的结果。在脆弱的医院人群例如重症监护、烧伤或者免疫缺陷的人群,以及在发生严重感染和之后发生败血症的风险升高的其它人群中,这是非常重要的。对暴露的个体的最早期征兆进行检测的快速和可靠的手段将是有利的,在战场和平民区中采用或怀疑使用生物武器就是一个例子。
本发明的描述
感染后,免疫系统的细胞通过激活识别和对病原产生应答。其导致不同信使蛋白的生成(例如细胞因子和化学因子(chemokine)),以及细胞表面活化标记物和粘附因子的表达(附图1)。这些的生成有利于细胞之间的通信,使得针对某一特定因素产生协同的免疫响应。由于该炎性免疫响应对于感染的响应相对恒定,并且发生在疾病过程的非常早的阶段,因此监测所述标记物表达的改变可以预测败血症发展的早期阶段。本发明的目的是提供一种在临床症状和疾病发生之前早期检测严重感染的方法,优选在干预治疗的窗口期(附图1)。
第一方面,本发明描述了在临床征象明显显现之前的几天检测感染、SIRS或败血症的早期迹象的系统和方法。其还提供了能够预测所述状况的临床病程的时间的方法。该系统包括分析基于生物样品的一组或多组检验的结果,优选血液样本。可选地,也可以包括其它常规临床手段进行分析。
第一组检验包括通过逆转录聚合酶链反应(RT-PCR)放大细胞因子mRNA来确定血白细胞的一组化学因子和细胞因子的表达水平。该组包括CD178(FAS-L)、MCP-1(单核细胞趋化蛋白-1)、TNF-α、IL-1β、IL-6、IL-8、IL-10、INF-α和INF-γ。CD178编码并表达为II型膜蛋白,但是可以被认为是一种细胞因子,原因是其被金属蛋白酶裂解而释放可溶性同质三聚体、可溶性Fasl或sFasl。
第二组检验包括确定一组表面标记物的细胞表面表达,包括HLA-DR、CD54、CD11b、CD31、CD69、CD80、CD62L(L-选择蛋白)、CD25和CD97。其可以通过标准方法实现,优选通过特异性的、经标记的抗体或其片段标记细胞来实现,以及通过流式细胞术分析。
该分析中可以包括的第三组检验包括常规的临床数据,包括体温、心率、总白细胞计数和白细胞分类计数(单核细胞、淋巴细胞、粒细胞、嗜中性粒细胞)、血小板计数、血清肌酐、尿素、乳酸、碱过剩、pO2、HCO3-和C-反应蛋白。
对一组或多组检验结果进行分析,优选通过神经网络程序进行分析,该程序能够对被采集样品的患者发展成败血症进行是/否的预测。
在阳性预测结果的情况下,优选进行进一步的分析,预测明显临床征象和症状发生的时间。或者,采用多元logistic回归进行两种分析。
可以逐个或同时进行检验组的分析。优选的是,将临床数据作为PCR数据的补充数据进入神经网络。同时,可以采用所述神经网络处理流式细胞数据。通过神经网络进行处理仅仅需要一组数据,但是输入一组、两组或全部的三组数据有好处,其原因是,这些额外的例子有助于“训练”神经网络,改善程序输出的可信度。
另一方面,采用神经网络处理前记录临床数据,或所述数据的数据库可以用于训练神经网络,提高其预测能力。
所述方法可以作为对重症监护病人的常规监测的一部分,而采集常规血液样品用作其它目的。还可以监测其它有感染和/或败血症倾向的医院患者。所述倾向条件包括先天的或获得性免疫缺陷(包括HIV/AIDS)或免疫抑制(例如普通外科手术患者、移植接受者或接受激素治疗的患者)、糖尿病、淋巴瘤、白血病或其它恶性肿瘤、贯穿或污染伤、烧伤或腹膜炎。另一方面,本发明的方法用于筛选传染性疾病爆发期间的个体,或者筛选已经或者被怀疑已经暴露于感染病原的个体,不论其是偶然的或是生物恐怖主义或者在武装冲突期间采用生物武器而故意导致的。
因此,本发明提供了一种筛选生物样品以检测感染、SIRS或败血症的早期阶段的方法,包括以下步骤:通过RT-PCR检测一组提供信息的生物标记物的表达和/或通过流式细胞术检测一组提供信息的细胞表面生物标记物的表达和/或监测一组标准临床测量结果;通过神经网络或多元统计分析来分析检测结果;根据明显感染发生的可能性和疾病发展的时机对所述样品进行分类。优选地,所述生物样品为血样。进一步优选地,对RT-PCR和/或流式细胞术和/或标准临床测量的结果进行第一分析,以产生关于SIRS或发展成败血症的可能性的预测。在一备选实施例中,其表达为概率。在一备选实施例中,其表达为二元是/否结果。
可选地,当第一分析提示可能发生SIRS或败血症(定义为超过预定的任意概率阈值,或“是”预测结果),对所述结果进行第二次分析,以确定发展到明显临床征象的可能时机,或者给出关于所述临床疾病的可能严重性的指示。
优选地,该组提供信息的通过RT-PCR检测表达的生物标记物由至少选自以下列表的4个、优选至少6个、更优选至少7个、最优选所有标记物组成,所述列表由CD178(FsaL)、MCP-1(或者称为CCL-2)、TNFα、IL-1β、IL-6、IL-8、IL-10、INF-α和INF-γ组成。
或者,通过流式细胞术检测其表达的提供信息的细胞表面生物标记物组由至少1个、优选至少2个、最优选至少三个选自以下列表的标记物组成,所述列表由CD31、HLA-DR、CD54、CD11b、CD69、CD80、CD62L、CD25和CD97组成。
当对标准临床测量结果进行分析时,其由至少一个、优选至少三个、更优选至少五个选自以下列表的测量结果组成,所述列表由体温、心率、总白细胞计数和白细胞分类计数(单核细胞、淋巴细胞、粒细胞、嗜中性粒细胞)、血小板计数、血清肌酐、尿素、乳酸、碱过剩、pO2、HCO3-和C-反应蛋白组成。
在一高度优选的实施方式中,通过神经网络进行所述分析。最优选的是,其为多层感知神经网络。
优选地,所述神经网络能够在大于70%的病例中正确预测SIRS或败血症(在对照组中不能通过预防性治疗防止其发展的试验中确定)、优选在至少80%病例中,更优选在至少85%病例中,最优选在至少95%病例中。优选地是,SIRS或败血症能够在明显临床征象发生的至少前一天预测到,更优选至少前两天,更加优选至少前三天,最优选在诊断出SIRS或败血症前的三天以上。
在另一优选实施方式中,通过多元统计分析进行所述分析,优选包括主成分分析和/或判别函数分析。更优选地是,所述多元统计分析包括判别函数分析。
另一方面,本发明提供了一种筛选生物样品的系统以检测感染、SIRS或败血症的早期阶段,其包括:从所述样品获得的细胞中提取和纯化RNA的装置、热循环控制装置或其它通过逆转录聚合酶链反应(RT-PCR)放大所选择的RNA序列的装置、检测并定量所述RT-PCR的结果的装置、能够分析所述结果的经训练的基于计算机的神经网络和由此能够向操作者交流神经网络分析结论的显示装置。
注意:在该方面,可以采用判别函数分析对RT-PCR的结果进行分析,但是神经网络为优选的实施方式。
或者,系统包括:对从生物样品获得的细胞上的特异性细胞表面标记物进行标记并对所述标记物的表达进行量化的装置、检测或对所述标记和量化的结果进行定量的装置、能够分析所述结果的经训练的基于计算机的神经网络和由此能够向操作者交流神经网络分析结论的显示装置。优选地,所述标记是通过标记抗体或抗体片段进行的,所述量化是通过荧光激活细胞分类(FACS)或其它形式的流式细胞术进行的。
优选地,所述系统包括同时进行RT-PCR和FACS分析的装置,优选地,还可以包括分析其它常规临床检验结果,并将所有这些结果结合至该分析中。
在另一方面,本发明提供了根据上述方法的任意实施方式的分析,以制备用于诊断SIRS、败血症或感染的诊断装置、或者采用上述系统以制备诊断感染的诊断装置。
还提供的是根据上述方法的早期诊断SIRS、感染或败血症的一种方法。
本发明的详细描述
参考以下附图和例子对本发明进行详细描述。
附图1:感染后,免疫系统的细胞通过激活识别并对病原做出应答。其导致了不同信使蛋白(例如细胞因子和化学因子)的产生,以及活化标记物和粘附分子在细胞表面的表达。这些的生成有利于细胞之间的通信,使得针对某一特定因素产生协同的免疫响应。由于该炎性免疫响应对于感染的响应相对恒定,并且发生在疾病过程的非常早的阶段,因此监测所述标记物表达的改变可以预测败血症发展的早期阶段。理想的是,在临床症状和疾病发生之前,在干预治疗的窗口期进行上述监测。
附图2:通过流式细胞术测量的CD31在粒细胞上表达的图形。来自诊断为败血症(n=6)的患者的前三天血样,以及没有发生成为败血症的ICU患者(n=24)的血样。每个符号表示对一个患者进行的一次测量。
附图3:根据表4、模型2设计分析临床数据的神经网络。WCC=白细胞计数,CRP=C-反应蛋白。
附图4:用金黄色葡萄球菌体外血液感染之后获得的细胞因子模式的改变。来自三个自愿者的血液数据详细见例8。
附图5:体外金黄色葡萄球菌败血症模型的神经网络分析的结果。
实施例1:通过对细胞因子表达、细胞表面标记物和临床测量结果的神经网络分析对败血症进行预测。
研究设计和患者
对Queen Alexandra医院ICU部门研究败血症的发生,产生了91名患者的队列(Dstl/CR08631)。在他们在ICU期间,每天从这些患者中采集血样,一共有24名患者被诊断为发展成败血症。对被诊断为败血症的当天(0天)采集的样品、诊断为败血症之前向后推六天的样品(-6天)采用RT-PCR和流式细胞术分别分析活化标记物和细胞因子mRNA的表达。此外,记录标准医院数据和临床观察结果。以同样方式对来自对照组患者的样品进行处理以提供供传统统计分析用的数据。
根据通常采用的实验室技术进行RT-PCR。简言之,对于血液样品,取全血,在RNA稳定剂存在下裂解细胞。采用珠亲和结合分离出RNA,其通过离心(或适当的话,采用磁法)进行分离,用去氧核糖核酸酶DNase消化除去污染的DNA,对RNA进行RT-PCT。
荧光激活细胞分类(FACS)流式细胞术是本领域非常公知的,可以采用任何标准方法。
数据分析
生物系统的复杂性以及本研究所采用的标记物之间的复杂关系使得标准的线性数据分析方法得出不确定的结果。因此,不清楚所述数据是否存在某种模式,由此寻求一种更加有效的能够非线性建模的方法以应对所述数据集的复杂性。
为了分析,对发生败血症之前1至4天的患者的数据进行核对,并与没有发展为败血症的ICU患者构成的年龄/性别匹配的对照组进行比较。单个样品提供测量多达56个不同参数的数据,将变量的选择性组合输入多层感知神经网络(Proforma,Hanon Solution,Glasgow,Scotland).
采用后向传播算法用随机选择的70%经平衡的败血症和对照数据对每个网络进行训练,然后用剩余的30%数据进行检验。在该网络内对数据进行建模时进行五次尝试,每个模型对数据归纳的能力不同。最成功的模型为最能够对之前未能预见的患者进行正确分类的模型,所述患者来自败血症或非败血症对照组。
结果
表2显示了对测试患者进行正确分类或“评分”29/35(或者82.9%)的成功模型的例子。
表2:采用细胞因子mRNA变量(第1至4天)的分类显示
得分 | 样品正确预测百分数% | |
患者总数 | 29/35 | 82.9% |
对照患者 | 16/20 | 80.0% |
败血症患者 | 13/15 | 86.7% |
为了提高该模型的可信度,对其实施五次,每次采用随机选择的不同数据,以训练和检验网络。一旦完成,对每个模型的得分进行平均,以得到关于所述网络将患者分类到正确的败血症或非败血症对照组的能力的全面指示。
一系列5个数据集的平均预测准确性为大约80%,如表3所示。
表3:仅采用RT-PCR数据预测败血症的神经网络(5个方案的分类性能分析)
每中/发生 | % | 命中/预测 | % | 概率% | 提高率% | 比率 | |
模型1 | |||||||
状况 | 23/27 | 85.2 | N/A | N/A | 50.0 | 35.2 | 1.7∶1 |
败血症 | 9/13 | 69.2 | 9/9 | 100.0 | 48.1 | 51.9 | 2.1∶1 |
对照 | 14/14 | 100.0 | 14/18 | 77.8 | 51.9 | 25.9 | 1.5∶1 |
模型2 | |||||||
状况 | 29/35 | 82.9 | N/A | N/A | 50.0 | 32.9 | 1.7∶1 |
败血症 | 13/15 | 86.7 | 13/17 | 76.5 | 42.9 | 33.6 | 1.8∶1 |
对照 | 16/20 | 80.0 | 16/18 | 88.9 | 57.1 | 31.7 | 1.6∶1 |
模型3 | |||||||
状况 | 32/40 | 80.0 | N/A | N/A | 50.0 | 30.0 | 1.6∶1 |
败血症 | 20/22 | 90.9 | 20/26 | 76.9 | 55.0 | 21.9 | 1.4∶1 |
对照 | 12/18 | 66.7 | 12/14 | 85.7 | 45.0 | 40.7 | 1.9∶1 |
模型4 | |||||||
状况 | 21/26 | 80.8 | N/A | N/A | 50.0 | 30.8 | 1.6∶1 |
败血症 | 9/12 | 75.0 | 9/11 | 81.8 | 46.2 | 35.7 | 1.8∶1 |
对照 | 12/14 | 85.7 | 12/15 | 80.0 | 53.8 | 26.2 | 1.5∶1 |
模型5 | |||||||
状况 | 23/29 | 79.3 | N/A | N/A | 50.0 | 29.3 | 1.6∶1 |
败血症 | 12/15 | 80.0 | 12/15 | 80.0 | 51.7 | 28.3 | 1.5∶1 |
对照 | 11/14 | 78.6 | 11/14 | 78.6 | 48.3 | 30.3 | 1.6∶1 |
128/157=81.5%
79.3+85.2+80+82.9+80.8=408.2/5=81.64%
表4列出了采用不同变量组合而构造的一系列网络的平均预测准确值。
最成功的模型是采用细胞因子mRNA表达和流式细胞术数据的CD31表达百分率构建的(平均准确率为81.0%,表3,模型1),临床数据得分也较高(80.4%,表3,模型2)。
表4:神经网络分析的结果
模型 | 标记物 | 预测准确性(%) |
1 | FasL,MCP-1,TNF-α,IL-1β,IL-6,IL-8,IL-10 | 81.6 |
2 | 肌酐,单核细胞,CRP,淋巴细胞,体温,嗜中性粒细胞,白细胞计数 | 80.4 |
3 | FasL,MCP-1,IL-8,白细胞计数,体温,肌酐 | 79.0 |
4 | FasL,MCP-1,TNF-α,IL-1β,IL-6,IL-8&IL-10,%CD31&肌酐 | 78.7 |
5 | FasL,TNF-α,IL-1β,IL-6,IL-8&IL-10 | 78.1 |
6 | FasL.MCP-1,TNF-α,IL-1β,IL-6,IL-8,IL-10,肌酐,单核细胞,CRP,淋巴细胞,体温,嗜中性粒细胞,白细胞计数 | 76.0 |
7 | MCP-1,TNF-α,IL-1β,IL-6,IL-8,IL-10 | 76.0 |
8 | 血小板,HCO3,PO2,尿素,肌酐,心率 | 70.7 |
为了进一步检测我们的预测模型,我们采用临床症状发生之前的第1至4天获得的100%细胞因子数据对网络进行训练。然而,我们选择包括“第0天”败血症患者的检验数据和来自第-5、-6天以及第0、5、6天的数据,还从单独的自愿者研究中选择14名对照患者,他们中的7个在采样的9天内出现上呼吸道感染(URTI)的症状(Dstl/CR08631)。所述结果显示在下表5中。
表5:在其它组的预测中细胞因子mRNA模型(第-4至-1天)的性能
检测组 | 得分 | 样品正确预测百分率 |
第0天败血症 | 8/9 | 89% |
败血症第-5天 | 7/9 | 78% |
败血症第-6天 | 5/6 | 83% |
该表显示了,我们的由临床败血症发生前4天所表达的模式而构建的模型正确地确认89%的第0天败血症患者、78%的第-5天败血症患者和83%的第-6天败血症患者,或进行“评分”。总体而言,采用神经网络的分析产生了大量的针对败血症的预测模型。仅仅采用细胞因子数据构建的模型证明能够持续成功地区分将发展成为败血症的个体和不会发展成为败血症的个体。
实施例2:采用神经网络模型使得非败血症志愿者中没有假阳性结果
表6显示了通过细胞因子RT-PCR对志愿者组的检测结果,该组中没有人出现SIRS或败血症表现。
表6
姓名 | 命中/发生 | % | 命中/预测 | % | 概率 | 提高 | 比率 |
总计 | 13/13 | 100.0 | N/A | N/A | 50.0% | 50.0% | 2.0∶1 |
对照 | 13/13 | 100.0 | 13/13 | 100.0 | 100.0% | 0.0% | 1.0∶1 |
败血症 | 0/0 | N/A | 0/0 | 0.0 | 0.0% | 0.0% | N/A |
实施例3:采用临床数据,神经网络的败血症预测准确度超过90%。
采用下表7中描述的参数,采用临床数据而检验的神经网络模型见表4的模型2,并进一步在附图3中描述:
表7:分析临床数据的神经网络参数
从输入至隐藏的权重 | ||||||
输入 | 单元1 | 单元2 | 单元3 | 单元4 | 单元5 | |
1 | -3.04389 | 0.783085 | -8.14579 | -5.31918 | -11.699 | |
2 | -2.40492 | -3.28341 | 6.81119 | -2.22889 | 10.3357 | |
3 | 20.0835039 | -2.25136 | 15.2575 | -3.38677 | 11.4842 | |
4 | -16.2918 | 3.45666 | -5.86258 | 8.93773 | -1.57967 | |
5 | 12.7828 | 45.2618 | 6.36791 | 7.45053 | 17.1156 | |
6 | 34.3201 | 7.31471 | -18.392 | 16.806 | 10.2057 | |
7 | -3.31337 | -1.41443 | -10.2639 | -3.68324 | 1.4483 | |
从偏性至隐藏的权重 | ||||||
1 | -0.407035 | -9.38879 | 4.06257 | -12.6877 | -10.7582 | |
从隐藏至输出的权重 | ||||||
1 | -8.79483 | 8.82026 | ||||
2 | 8.79794 | -8.38064 | ||||
3 | 3.60783 | -3.62177 | ||||
4 | 3.00073 | -3.83778 | ||||
5 | 4.85609 | -4.85559 | ||||
从偏性至输出的权重 | ||||||
1 | -2.64015 | 2.65373 |
表8
名称 | 命中/发生 | % | 命中/预测 | % | 概率 | 提高 | 比率 |
状况 | 14/15 | 93.3 | N/A | N/A | 50.0% | 43.3% | 1.9∶1 |
败血症 | 8/8 | 100.0 | 8/9 | 88.9 | 53.3% | 35.6% | 17∶1 |
对照 | 6/7 | 85.7 | 6/6 | 100.0 | 46.7% | 53.3% | 2.1∶1 |
实施例4:采用人工免疫识别系统
表述
初始AIRS系统(A.Watins.An Artifical Immune RecognitionSystem.密西西比州立大学:MSc Thesis.,2001)采用简单的实值形状空间。近来,其它人将所述表述延伸至Hamming形状空间(J.Hanamaker和L.Boggess.The effect of distance metrics on AIRS.InProc.Of Congress on Evolutionary Computation(CEC).IEEE,2004)和自然语言(D.Goodman,L.Bogess和A.Waktins.“An investigation intothe source of power for AIRS,an artifcial immune classificationsystem”.In Proc.Int Joint Conference on Neural Networks,第1678-1683页.IEEE,2003)。AIRS保持了一组人工识别球(ArtificialRecognition Balls,ARBs),其含有被学习的数据矢量,刺激水平和大量资源。在训练期间,通过评价ARB中数据矢量对训练项目的亲和性来计算所述刺激水平,匹配越强,刺激越大。采用该刺激水平以规定ARB如何产生多个克隆,并如何影响ARB的生存。
亲和性测量
其依赖于所采用的表述。已经提出大量的亲和性测量以用于AIRS,包括Hamming距离、Euclidean距离等。在该研究中,同时采用Euclidean和Hamming距离量度,其中Euclidean的结果最佳。
免疫算法
实质上,AIRS涉及两个群体,一个存储库和一个ARB库C。其有独立的训练和检验相,检验相类似于K-最近领域分类器。在该训练相中,将训练数据相呈给M。然后随机启动该数据集,实验证据将建议AIRS对初始启动点不敏感。所述训练项目与集M的所有存储单元匹配,确定一个单元为较高匹配Mcmatch。然后克隆并突变该Mcmatch。克隆与刺激呈正比(刺激越高,克隆率越高),突变与刺激呈反比(刺激越高,突变率越低)。将这些克隆插入ARB库,C。然后将该训练项目呈给ARB库的成员,其中采用循环程序从而对新的侯选存储单元进行克隆和突变。通过群控制处理,生存由ARB可以要求的资源数量确定,产生新的侯选存储单元。该机制基于J.Timmis和M.Neal.A resource limited artificial immune system.Knowledge BasedSystemsm 14(3/4):121-130,2001中提出的资源分配算法。用该训练项目,将该新的候选单元与Mcmatch相比。如果所述侯选单元和MVCmatch之间的亲和性更高,则该存储单元被所述侯选单元取代。
对每个训练项目都实施该过程,存储集将含有大量能够用于分类的单元。以k-最近领域方式进行未见数据项的分类。
实验设置
尝试采用与将神经网络用于该数据集具有可比性的实验程序。对于所有研究,采用标记物:asL,MCP-1,TNF-α,IL-β,IL-6,IL-8和IL-10。然而,不可能完全精确地复制该数据集,其原因是在神经网络的研究过程中就数据的前处理而言其信息是不完整的。
实验1
在第一组实验中,采用从败血症发生前1至4天的患者所收集的数据,以及作为训练数据从对照组患者中收集的数据。尤其是,采用显示败血症征象患者的第1、2、3和4天的组合数据,随机收集对照组患者,以训练AIRS。一共采用了59个训练数据项目。为了检验AIRS,采用对对照组患者的随机收集,以及来自所有天的组合数据(排除已经在训练过程中用过的数据)。一共采用了34个检验数据点。AIRS的设置参见表9:
表9:AIRS的参数设置
参数 | 设置 |
Epochs克隆率 | 10 |
突变率 | 0.8 |
群的初始大小 | 5 |
亲和阈值 | 0.2 |
刺激阈值 | 0.8 |
源数量 | 200 |
实验2
对于我们的第二组实验,采用第0、5和6天的数据和对照患者对显示败血症征象的患者进行分类。采用与实验1想同的数据训练AIRS系统,同时采用相同的参数。
结果
所述结果与通过神经网络分析得到的结果不能直接相比,其原因是,很难从原始报告中严格确定在一段天数期间数据是如何首先组合,然后又如何分为训练和检验组的。因此,考虑所获得的结果时应当记住上述问题。
实验1
对AIRS算法进行独立运行,然后计算平均和标准偏差。发现,AIRS能够达到平均73(2.96)%的分类准确性。其与神经网络分析大约低10%(采用相同的标记物)。然而,直接比较时应当注意。
实验2
再次,对AIRS算法进行10次独立运行,然后计算平均值和标准偏差。此次,分析了败血症发作之前的天数(0、5和6)和对照组。再次,用第1至4天和对照组的数据对AIRS进行训练。所述结果参见表10。
表10:在其他组中预测(括号内为标准差)
检验组 | 准确度 |
败血症的第0天 | 83(7.6)% |
败血症的第5天 | 96(6.1)% |
败血症的第6天 | 91(8.2)% |
对照 | 70(12.8)% |
从表10可以看出,AIRS确定了高百分率的败血症患者(在第5和6天优于神经网络,但是应当再次注意上述可比性警告)。对照组修整的不是很好,低于预期结果,并且显著低于神经网络法。其可能是由于用于训练的败血症患者的数据多于非败血症患者的数据从而AIRS偏向于败血症患者。
结论
AIRS似乎能够事先确定潜在的败血症病例,在一定水平上与神经网络法具有可比性。
实施例5:采用CD31表达预测败血症研究设计和患者
参见例1。
流式细胞术
将血液收集至肝素钠容器中(HMJ&S,Chessington,Surrey),并在室温下运至实验室。将100ul等分血样与免疫荧光染料混合,其体积采用生产商(Beckman Coulter Limited,High Wycombe,Buckinghamshire,和Becton Dickinson UK Limited,Cowley,Oxford)所推荐的。通过CD3和CD4共染色鉴别T辅助细胞,通过同时染色CD3和CD8鉴别T细胞毒细胞。对这些细胞集落染色HLA-DR,CD25,CD54和CD69。通过CD19染色鉴别B细胞,并通过CD80,CD86,CD25,CD54检查。用CD56染色区分并用CD11b,CD25,CD54和CD69检查自然杀伤细胞。通过CD14染色选择单核细胞集落,这些细胞用CD11b,CD54,CD80,CD86和HLA-DR染色作为探针。采用门控作用鉴别粒细胞集落,其对CD11b,CD69,CD31,CD54和CD62L染色。将这些染料在室温下孵育20分钟。在每管中加入500ul Optilyse C(Beckman Coulter Limited),并立即进行涡旋混合。在室温下孵育样品10分钟,以裂解红细胞,然后加入500ul lsoton(Beckman Coulter Limited),以固定染料。立即对管进行涡旋混合,在室温下孵育10分中。然后在Beckman CoulterEpics XL System 2 Flow Cytometer中对细胞进行计数。
统计
采用11.0版本的SPSS软件包的二元Logistic回归模型对数据进行分析。该分析将对照组的免疫调节剂表达的均值与败血症患者的7个时间点的均值进行比较:诊断之前6天、诊断前的第5、4、3、2、1天,以及0天,诊断败血症的当天。当数据点缺失时,用该组的平均值代替,以保证可接受的n值。如果代替的数据点不涉及被该模型高亮显示为可能预测物的标记物时,仅仅报告模型结果。
结果
数据分析发现弱的证据,预测效果=0.114。CD31的表达减少指示,诊断前3天是败血症的可能预测物,p=0.037(n=6)。诊断前6天获得的结果是不确定的,其原因在于该时期的样本量很小(n=4)。诊断前5、4、2或1天或者诊断当天没有发现有统计学显著意义的预测物。
结论
对从发展成为败血症之前的患者和没有发生该疾病的患者获得的流式细胞数据进行了比较。采用诊断败血症之前的患者的结果构建组,对照组由没有发展成为败血症的年龄匹配的患者所测量的结果组成。由于数据的分布,显示中位数和第90以及第10百分数的条图很难解释,因此进行了统计分析。
当对原始数据进行做图时(参见附图2),可以看出,6个败血症患者中4个(66.6%)的CD31的表达低于对照组。可以看出,对照组(无败血症)的数据点在11.8%至100%之间分布,而诊断测量结果前三天内的6个数据点中的4个低于9%。因此,CD31可以在出现临床征象和症状之前的三天预测败血症的发生。其提示了,CD31可能是一种有用的预测标记物,尤其是结合其他提供信息的败血症标记物。
实施例6:多变量统计分析以预测败血症
介绍
对出现败血症症状前1-6天的患者收集的数据进行多变量数据分析程序。测量包括流式细胞术、PCR和经典临床观察资料。对数据矩阵采用主成分分析(PCA),对3类观察资料的每个进行单独考虑,并组合成完整的数据集。采用判别函数分析(DFA)以确定就变量的均值而言组之间是否不同,然后采用该变量预测组的从属关系(例如新的病例)。对来自PCA的结果和三类观察资料的每个以及组合成的完整数据集进行上述分析。
数据描述、处理和多变量方法
在PCA之前,通过产生概率密度函数对数据进行总结。由于PCA和DFA之前需要正态分布,因此,采用Johnson转换算法对非正态数据进行转换。
PCA为减少维数的方法,其试图将多变量数据矩阵分解为一些潜在变量,其由变量的线性组合构成,其解释了原始矩阵的方差的主要部分。这样,可以确定数据集内参数之间的关联(正或负)。
实质上,DFA与方差分析法(ANOVA)类似。可以将DFA问题转化成一元分析ANOVA问题。尤其是,人们可以询问就某一特定变量的均值而言两个组之间是否显著不同。然而,应当清楚的是,如果一个变量的均值在不同组间显著不同,则可以做出结论,该变量在不同组间不同。
在单个变量的情况下,变量在不同组间是否不同的最终显著检验称为F-test。F实质上为数据组间方差与总(平均)组内方差的比率。如果组间方差显著较大,则均值之间必定有显著差异。
当考虑多变量时,可以确定几个变量中的哪个对组间差异做出贡献。其导致总方差和协方差的矩阵;以及总组内方差和协方差的矩阵。通过多变量F-test比较这些矩阵,以确定组间是否存在显著差异(就所有变量而言)。该过程与多元方差分析或MANOVA相同。如同在MANOVA中,进行多元检验,如果有统计学显著性,检查变量中的哪些在组间具有显著不同的均值。这样,即使多个变量的计算更加复杂,但是主要推理仍然适用,即寻求组间不同的变量,如所观察到的均值差一样显而易见。
采用完整数据矩阵(包括替代的均值)对临床、流式细胞术和RT-PCR数据进行DFA,将含有替代均值的一个或多个参数的数据点排除。采用模拟模型对由模型1得到的模型进行PCA得分分析。后者的目的是确定转化的数据矩阵(PCA)是否可以对观察资料进行分类。
模型1观察数据的主成分分析(PCA)
i)基于临床数据的PCA模型
在给定模型中衍生和使用的PC数量通常确定为特征值>1.6的成分。对于临床数据PC满足该条件,并且解释了数据集中74.3%的方差。由于每个PC与其它PC正交(不相关),一个临床参数与特定PC之间的关联确定了该PC,并解释了该参数如何影响数据集的方差。
表12对每个参数的负荷进行了总结,从临床数据得到6个PC。负荷值>0.5指示一特定参数对给定PC的强贡献。源自该数据集的PC可以解释如下:
PC1 其由WCC、单核细胞、嗜中性粒细胞和血小板之间的强关联控制。肌酐和乳酸之间存在强关联。这些组具有负关联(PC1比例尺的相反端)。BXS和HCO3 -高度相关,对PC1和PC2的贡献相同。后面这些参数在PC1中与肌酐和乳酸相反。表2和附图2中总结了这些关联。
PC2 显示了由WCC、单核细胞、嗜中性粒细胞和血小板构成的组与由BXS和HCO3 -构成的组呈负关联。
PC3 该PC的特征在于体温、HR和CRP之间的强关联,如附图3所示
PC4 尽管对于该PC,许多参数显示显著性,但是仅仅有CRP明确与该PC相关,如附图4所示
PC5 在该组分中,pO2与尿素和MAP相反。
PC6 该PC专门解释了由淋巴细胞引入该数据集的方差。
在解释PC负荷时,相关联的临床参数提示,这些样本/物理参数的水平在属于PC1,2等的患者中升高或降低。其将在判别分析部分中进行。
ii)基于流式细胞数据的PCA模型
为了清楚和简便起见,对参数进行简写,如表13所示。9个PC解释了该数据集的80.5%方差,如表14的特征值矩阵所示,表15中总结了相关负荷。所测量的参数之间的关联如附图16所示,在衍生的PC中以下之间具有强关联:
fc1和fc3
fc5和fc6
fc7和fc8
fc9与fc10和fc11,
fc12和fc13
fc11和fc14
fc17与fc20和fc23
fc21与fc17,20和22
fc23与fc20和fc22
fc28和fc29
所述PC结构可以解释如下:
PC1 | CD3 CD4中的CD3 CD4 CD25,CD3 CD4中的CD3 CD4HLA-DR,CD3 CD8中的CD3 CD8 CD25,CD19中的CD19CD80,CD19中的CD19 CD86,CD14中的CD14 CD80,CD14中的CD14 CD86,CD19中的CD19 CD54,CD19中的CD19CD25和CD56中的CD56 CD54相关 |
PC2 | CD19中的CD19 CD86与CD14中的CD14 HLA-DR,CD14中的CD14 HLA-DR CD11B,CD14中的CD14 HLA-DRCD11B CD54相反 |
PC3 | CD3 CD4中的CD3 CD4 CD54,CD3 CD4中的CD3 CD4CD69,CD3 CD8中的CD3 CD8 CD54,CD3 CD8中的CD3CD8 CD69 |
PC4 | CD56中的CD56 CD59,CD69(%),CD11B CD69(%) |
PC5 | CD3 CD8 CD62L(%)中的CD3 CD8 CD54, |
PC6 | CD31(%) |
PC7 | 没有EV>0.5的显著组分 |
PC8 | CD54(%) |
PC9 | CD14中的CD14 CD11B |
通过仅仅考虑一对或一组中的一个参数,则可以除去9个参数,由此增加了该数据矩阵的y组分。然而,确定的是,仅仅对CD31(%)、CD54(%)、CD62L(%)、CD11B%)、CD69(%)和CD11B CD69(%)进行统计分析(fc24-fc29)。
所选择的流式细胞变量的特征值矩阵显示在表17中,表18总结了所构建的PCA模型的负荷。采用3个PC解释了该数据集的76.6%的方差。该PC模型显示:
PC1 | CD69(%)与CD11B CD69(%)相关 |
PC2 | CD31(%),CD62L(%)和CD11B(%)相关 |
PC3 | 由与CD54(%)相关的方差组成 |
iii)基于RT-PCR数据的PCA模型
表19指示,仅仅3个PC解释了RT-PCR数据的72.9%方差。该模型的负荷参见表20。变量与每个PC之间的关联见附图23和24,并显示了以下:
PC1 | Fas-l、MCP-1、TNF-α、IL-6和IL-8相关 |
PC2 | IL-1和IL10相关 |
PC3 | IL-1和IL10相反 |
IL-1和IL-10在PC2中相关而在PC3中相反是非常有意思的。其显示,在某些病人中,这些变量可能高度相关或者相反,可能提供了一种有效的区分患者的方式。
iv)基于组合的临床、流式细胞和RT-PCR数据的PCA模型
表21总结了包括在最终模型中的参数和关联矩阵的相关特征值。表21指示,9个PC的特征值大于1,其解释了数据方差的68.7%。该模型的负荷见表22。分析允许对PCA模型进行如下解释:
PC1 | 显示WCC、嗜中性粒细胞、APTR、HCO3-、BXS、血小板、CD69(%)和CD 11B CD69(%)呈正相关。该PC还使上述与相互关联的乳酸和肌酐相反 |
PC2 | CD69(%)、CD 11B CD69(%)、WCC、嗜中性粒细胞、单核细胞和INR相关联。这些参数与TNF-α相反 |
PC3 | PCR参数Fas-L、MCP-1、TNF-α、IL-6和IL8强相关 |
PC4 | CRP和IL-10相反 |
PC5 | 流式细胞参数CD31(%)、CD54(%)、CD62L(%)和CD69(%)相关 |
PC6 | CD62L(%)与HR相关 |
PC7 | 与体温相关 |
PC8 | 与IL-1相关 |
PC9 | 与PO2相关 |
模型2:基于观察资料和PCA得分数据的判别函数分析(DFA)
在下表23中解释了在DFA模型中所有模型定义中共有的术语,并显示了数值。
模型
分析的目的是建立如何能最佳地预测病例属于哪个组的“模型”。在以下讨论中,采用术语“在模型中”以指被包括在组员预测中的变量,如果它们没有被包括,则“不在模型中”。
●前向逐步分析
在逐步判别函数分析中,逐步建立判别模型。尤其是,在每步中,对所有变量进行检查并评价,以确定哪个对区分不同组的贡献最大。然后将该变量包括在该模型中。
●后向逐步分析
可以后向逐步分析;在此情况下,程序首先在模型中包括所有变量,然后在每一步中,除去对预测组员贡献最小的变量。这样,由于成功的判别函数分析,将在模型中仅仅保留“重要”的变量,也就是说,对区分不同组的贡献最大的变量。
●进入F,出去F
所述逐步法由各自进入的F和出去F值“指导”。一个变量的F值指示其在组间判别时具有的统计学显著性,即其是一个变量对预测组员的独特贡献的量度。总的来说,该程序继续选择将被保留在模型中的变量,只要这些变量的各自F值大于指定进入的F;并且,如果变量的显著性低于指定排出的F,则从模型中将这些变量排除(除去)。
●容许度
一个变量的容许值计算为该变量与模型中其它所有变量的1-R2。这样,容许度为各个变量的冗余度的量度。例如,容许值为0.10意味着该变量与该模型中的其它变量的冗余为90%。
●Wilks’λ
该参数给出了该模型的判别能力的量度,可以评价范围为0(完全判别)至1(不能判别)的值
●部分λ
这是与各个变量对该模型判别能力的独特贡献(正交测量)相关的Wilks’λ。
就这一点,逐步判别分析的结果的常见错误解释是取表面值的统计显著性值。当该程序确定在下步分析中哪个变量被包括或排除时,其实际计算每个被考虑的变量的贡献的显著性。因此,本质上,由于逐步法“挑挑选选”将被包括在模型中的变量,以产生最大区别,则所述逐步法利用了几率。这样,当采用逐步法时,应知晓显著性水平不反映实际的 错误率,即必须保持的错误拒绝HO(组间无区别的无效假设)的概率。
●典型相关分析(CCA)
这是评价变量之间关系的额外程序。尤其是,该操作允许对两组变量之间的关系进行说明。以下对表征该分析的参数进行详细说明。
●根的显著性(χ2检验)
术语根用于描述单个判别函数(DF)。衍生的DF的统计学显著性通过连续DF的χ2检验来检测。得到所有典型根的逐步减小的检验报告,该报告含有所有DF的显著性,然后在除去第一根之后,得到报告剩余根的显著性的第二条线,等等。由此得到解释的DF数量。
●判别函数系数
产生两个输出,一个作为原始系数,一个作为标化系数。此处,原始意味着所述系数可以与所观察的数据结合以计算(原始)判别函数得分。标化系数是指通常用于解释的系数,其原因是,它们与标准变量有关,因此是指可比的尺度。
●特征值
得到每个DF的特征值和每个函数所解释的方差的累积比例。该值在PCA和DFA中定义的方式相同。该值越大,DF解释的方差的量越大。
●因子结构系数
这些系数代表了变量和DF之间的关联,并通常采用以解释判别函数的“含义”。以类似于PCA的方式,对因子的解释应当基于因子结构系数。
●典型变量的均值
当获知变量如何参与不同组间的判别时,下一逻辑步骤是确定每个DF的判别性质。回答该问题的第一步是检查典型均值。对于给定DF,其典型均值越大,该DF的判别能力越大。
i)基于临床数据的PCA得分的DFA模型
a)含有替代均值
表19总结了该分析的结果。0.4的Wilk’λ值指示相对无效的分类模型。所述三个衍生的DF共解释了该数据集方差的89.9%,DF主要由PC1,3和4组成。因子结构系数指示:
●DF1由PC1和较少程度上PC4所解释的方差组成
●DF2专门由PC3解释的方差组成
●DF3由PC4解释的方差组成
标化系数确定了这些关联。典型变量的均值指示:
●DF1使第1,2和3天与第5和6天负相关
●DF2确定了对照组的观察结果
●DF3确定第2天的观察结果
该模型的分类及其与PC相关的判别性质的总结见表25。该模型的分类矩阵见表26。表26提示,可以对对照组和第6天的数据进行良好的分类,分别有80%和83%的观察结果被正确分类。然而,该模型的总分类能力较差,所有观察结果中仅有48%被正确分类。
b)排除替代均值
表24总结了该分析的结果。0.45的Wilk’λ值指示相对无效的分类模型。三个衍生的DF共解释了该数据集方差的95%,DF主要由PC1,3和5组成。因子结构系数指示:
●DF1由PC1和PC5之间的负关联所解释的方差组成
●DF2专门由PC3解释的方差组成
●DF3由PC3解释的方差组成,但是其低于DF2
标化系数确定了这些关联。典型变量的均值指示:
●DF1使第1,2和3天与第5和6天负相关
●DF2使第1和6天与对照组负相关
●DF3确定第6天的观察结果
该模型的分类及其与PC相关的判别性质的总结见表25。该模型的分类矩阵见表26。表26提示,可以对对照组和第6天的数据进行良好的分类,分别有83%和67%的观察结果被正确分类。然而,该模型的总分类能力较差,所有观察结果中仅有44%被正确分类,低于用均值代替的变量。
采用均值替代值和不采用均值替代值的类似预测效率证实了接下来包括这些值的进行DFA的方法。
ii)基于临床数据的转换值的DFA模型
为了提高观察资料的分类,从原始数据矩阵转化的变量值经过DFA。其主题是由于PCA为维素减少的方法,可能某些数据量丢失了,对PCA得分进行DFA导致预测力低的模型。
表27总结了该分析的结果。0.22的Wilk’λ值指示PCA得分分类模型的改进。五个衍生的DF共解释了该数据集的99%方差,DF由BXS、CRP、乳酸、尿素、体温、肌酐、嗜中性粒细胞、pO2和HCO3 -以及其它对观察资料的分类没有影响的临床变量组成。因子结构系数指示:
●DF1对与乳酸呈负相关的BXS和HCO3 -之间的关联进行分类
●DF2对显示BXS、CRP和HCO3 -之间呈高度关联的观察资料进行分类
●DF3对体温和肌酐具有负相关的样本进行分类
●DF4对体温和pO2具有负相关的样本进行分类
●DF5对尿素和嗜中性粒细胞具有负相关的样本进行分类
该模型的分类及其与临床变量相关的判别性质的总结见表31。该模型的分类矩阵见表32。表32提示了可以对对照组和第6天的数据进行良好的分类,分别有80%和83%的观察结果被正确分类。与PCA得分模型相比,第1,2和5天显著改善,但是该模型的总体分类能力较差,所有观察结果中有55%被正确分类。
iii)基于流式细胞的PCA得分的DFA模型
表33总结了该分析的结果。0.39的Wilk’λ值指示相对无效的分类模型。两个衍生的DF共解释了该数据集的71%方差,DF主要由PC1,3和5组成。因子结构系数指示:
●DF1由PC1和PC8所解释的方差组成
●DF2专门由PC5解释的方差组成
标化系数确定了这些关联。典型变量的均值指示:
●DF1使第1天与第5和6天负相关
●DF2确定第3天的观察结果
该模型的分类及其与PC相关的判别性质的总结见表34。该模型的分类矩阵见表35。表35提示了可以对对照组和第6天的数据进行良好的分类,两个组中均有66%的观察结果被正确分类。但是该模型的总体分类能力较差,所有观察结果中仅有44%被正确分类。
iv)基于流式细胞数据的DFA模型
该模型的分类及其与变量相关的判别性质的总结见表36。0.034的Wilk’λ值指示优良的分类模型。三个衍生的DF共解释了该数据集的74%方差,DF主要由fc7-8、fcl1、fc16、fc25、fc28、fc29组成。因子结构系数指示:
●DF1将fc7,16,28相关,而这些与fc8和25相反
●DF2将fc7,8,16和25相关,而这些与fc12相反
●DF3与fc11相关
表37总结了该信息。
这些关联被标化系数证实。典型变量的均值指示:
●DF1使第1天与第5和6天相反
●DF2使对照组与第3,4和5天相反,将对照组与第6天相关
●DF3使第2和3天与第4天相反
表38提示可以对所有组进行良好分类。该模型的总分类能力令人印象深刻,全部观察资料中的76.6%被正确分类。
对于PCA得分和转换的数据,RT-PCR的DFA模型较差,Wilk’λ值>0.8,它们被放弃,不再考虑。
v)基于组合的临床、流式细胞和RT-PCR数据的DFA模型
该模型的分类及其与变量相关的判别性质的总结见表39。0.0087的Wilk’λ值指示优良的分类模型。四个衍生的DF共解释了该数据集的89.2%方差,DF因子结构指示:BXS、fc25、fc22、fc11、体温、CRP、fc18、fc6、IL-6、INR、APTR、fc16、尿素、乳酸、Fas-L、fc13、fc1、fc3、MCP-1、fc28、IL-10、fc27、fc26、嗜中性粒细胞、fc14、WCC、fc29、血小板、pO2被包括在该模型中。所有其它参数不能满足逐步法标准,因此从模型中删除。
典型变量的均值指示:
●DF1使第1天与第5和6天相反
●DF2使对照组与第6天相关,将这些与第3,4和5天相反
●DF3使对照组与第5天相关,将这些与第2,3和6天相反
●DF4使第4和5天相反
表40显示了能够得到对于所有组的良好分类,最小准确分配率为76%。该模型的总分类能力令人印象深刻,全部观察资料中的86.9%被正确分类。
当采用该模型,对数据应用每个DF时,可以清楚地对患者组进行聚类,然后与其它组空间分离。
结论
PCA对所有类型的患者的测量变量之间具有高度关联。从分子生物学角度而言,这些关联中的多个是预料中的。某些PCA模型极大地减少了数据集的维素,但是所得到的得分不能空间地分离患者组。
对从PCA获得的得分的DFA显示令人失望的结果。当采用PCA得分时,该模型的判别能力为44-56%。这些模型的低判别能力可能是PCA过程中数据集维素的减少和大量细节丢失的结果。在DFA中采用经转换的变量得到改进很多的模型。临床和流式细胞模型的判别能力分别为55%和76%。当对完整数据集(临床、流式细胞和RT-PCR变量)进行DFA时,观察到86.9%的预测效率。因此,推荐测量并采用包括在后一模型中的变量(表36)以对怀疑对败血症易感的新患者进行分类。
该模型令人印象最深刻的特征是其对症状发作前6天的患者正确分类的能力。因此可以监测关键的判别变量,并确定阈值水平,此时必须给予医学治疗。采用表34显示的参数,可以从患者获取数据,并采用转换算法,将数据输入DFA模型。然后,其能够将患者分类至适当的组中,有效率接近90%。当用于临床实验室中,其可能具有很大的价值。
表11关联矩阵的特征值,以及对临床观察资料的相关统计
PC | 特征值 | 总方差的的百分率% | 累积特征值 | 累积方差% |
1 | 3.73 | 24.87 | 3.73 | 24.87 |
2 | 2.50 | 16.70 | 6.23 | 41.58 |
3 | 1.53 | 10.23 | 7.77 | 51.81 |
4 | 1.22 | 8.14 | 8.99 | 59.95 |
5 | 1.10 | 7.39 | 10.10 | 67.35 |
6 | 1.04 | 6.96 | 11.14 | 74.31 |
7 | 0.87 | 5.80 | 12.01 | 80.11 |
8 | 0.77 | 5.19 | 12.79 | 85.31 |
9 | 0.66 | 4.42 | 13.46 | 89.73 |
10 | 0.54 | 3.62 | 14.00 | 93.36 |
11 | 0.44 | 2.99 | 14.45 | 96.35 |
12 | 0.29 | 1.95 | 14.74 | 98.30 |
13 | 0.210 | 1.40 | 14.95 | 99.70 |
14 | 0.02 | 0.16 | 14.98 | 99.87 |
15 | 0.01 | 0.12 | 15.00 | 100.00 |
表12临床测量对每个PC的负荷(与显示95%CL的特征值>0.5相关)
临床测量 | PC1 | PC2 | PC3 | PC4 | PC5 | PC6 |
温度 | 0.14 | 0.08 | 0.60 | -0.36 | 0.14 | 0.40 |
心率 | 0.15 | -0.06 | 0.58 | -0.48 | 0.11 | -0.37 |
MAP | 0.22 | 0.21 | -0.29 | -0.48 | -0.54 | -0.14 |
WCC | 0.63 | -0.74 | 0.00 | 0.05 | 0.02 | -0.01 |
嗜中性粒细胞 | 0.58 | -0.75 | 0.01 | 0.17 | 0.00 | 0.05 |
淋巴细胞 | 0.35 | 0.00 | -0.14 | -0.16 | -0.19 | -0.62 |
单核细胞 | 0.62 | -0.66 | 0.08 | -0.05 | 0.01 | -0.01 |
血小板 | 0.53 | -0.09 | -0.37 | 0.05 | 0.22 | 0.04 |
CRP | -0.15 | 0.11 | 0.62 | 0.50 | 0.08 | -0.36 |
PO2 | -0.26 | 0.07 | -0.36 | 0.12 | 0.53 | -0.39 |
HCO3- | 0.76 | 0.47 | 0.16 | 0.23 | -0.13 | -0.07 |
BXS | 0.78 | 0.49 | 0.15 | 0.20 | -0.11 | -0.05 |
乳酸 | -0.55 | -0.46 | 0.16 | -0.12 | -0.08 | -0.22 |
尿素 | -0.30 | -0.17 | 0.09 | 0.44 | -0.58 | 0.04 |
肌酐 | -0.69 | -0.44 | 0.05 | 0.02 | -0.19 | -0.06 |
表13流式细胞数据分析中采用的缩写
缩写 | 流式细胞参数 |
fc1 | CD3 CD4中的CD3 CD4 CD25 |
fc2 | CD3 CD4中的CD3 CD4 HLA-DR |
fc3 | CD3 CD8中的CD3 CD8 CD25 |
fc4 | CD3 CD8中的CD3 CD8 HLA-DR |
fc5 | CD3 CD4中的CD3 CD4 CD54 |
fc6 | CD3 CD4中的CD3 CD4 CD69 |
fc7 | CD3 CD8中的CD3 CD8 CD54 |
fc8 | CD3 CD8中的CD3 CD8 CD69 |
fc9 | CD19中的CD19 CD80 |
fc10 | CD19中的CD19 CD86 |
fc11 | CD14中的CD14 CD80 |
fc12 | CD14中的CD14 CD86 |
fc13 | CD19中的CD19 CD54 |
fc14 | CD19中的CD19 CD25 |
fc15 | CD56中的CD56 CD54 |
fc16 | CD56中的CD56 CD69 |
fc17 | CD14中的CD14 CD54 |
fc18 | CD14中的CD14 HLA-DR |
fc19 | CD14中的CD14 CD11B |
fc20 | CD14中的CD14 CD54 HLA-DR |
fc21 | CD14中的CD14 HLA-DR CD11B |
fc22 | CD14中的CD14 HLA-DR CD11B CD54 |
fc23 | CD14中的CD14 CD11B CD54 |
fc24 | CD31(%) |
fc25 | CD54(%) |
fc26 | CD62L(%) |
fc27 | CD11B(%) |
fc28 | CD69(%) |
fc29 | CD11B CD69(%) |
表14关联矩阵的特征值,以及对流式细胞数据的相关统计
PC | 特征值 | 总方差的百分数% | 累积特征值 | 累积方差% |
1 | 7.16 | 24.68 | 7.16 | 24.68 |
2 | 3.81 | 13.12 | 10.96 | 37.80 |
3 | 2.63 | 9.07 | 13.59 | 46.88 |
4 | 2.43 | 8.39 | 16.03 | 55.27 |
5 | 1.90 | 6.54 | 17.93 | 61.81 |
6 | 1.62 | 5.59 | 19.55 | 67.41 |
7 | 1.46 | 5.05 | 21.01 | 72.45 |
8 | 1.21 | 4.19 | 22.23 | 76.64 |
9 | 1.10 | 3.78 | 23.32 | 80.42 |
10 | 0.87 | 3.00 | 24.19 | 83.42 |
11 | 0.85 | 2.93 | 25.04 | 86.36 |
12 | 0.71 | 2.44 | 25.75 | 88.80 |
13 | 0.63 | 2.17 | 26.38 | 90.97 |
14 | 0.53 | 1.82 | 26.91 | 92.78 |
15 | 0.39 | 1.36 | 27.30 | 94.14 |
16 | 0.35 | 1.21 | 27.65 | 95.35 |
17 | 0.34 | 1.16 | 27.99 | 96.51 |
18 | 0.31 | 1.07 | 28.30 | 97.58 |
19 | 0.20 | 0.88 | 28.50 | 98.26 |
20 | 0.12 | 0.42 | 28.62 | 98.68 |
21 | 0.10 | 0.34 | 28.72 | 99.02 |
22 | 0.08 | 0.27 | 28.79 | 99.29 |
23 | 0.07 | 0.23 | 28.86 | 99.52 |
24 | 0.05 | 0.18 | 28.91 | 99.70 |
25 | 0.04 | 0.13 | 28.95 | 99.83 |
26 | 0.02 | 0.08 | 28.98 | 99.91 |
27 | 0.01 | 0.05 | 28.99 | 99.96 |
28 | 0.01 | 0.03 | 29.00 | 99.99 |
29 | 0.00 | 0.01 | 29.00 | 100.00 |
表15所有流式细胞数据对每个PC的负荷(与显示95%CL的特征值>0.5相关)
PC1 | PC2 | PC3 | PC4 | PC5 | PC6 | PC7 | PC8 | PC9 | |
fc1 | 0.59 | -0.36 | 0.19 | 0.06 | -0.13 | -0.28 | 0.28 | -0.11 | -0.33 |
fc2 | 0.61 | -0.07 | -0.20 | -0.21 | 0.11 | 0.41 | 0.34 | -0.01 | -0.17 |
fc3 | 0.62 | -0.31 | -0.07 | 0.14 | 0.14 | -0.11 | 0.21 | 0.05 | -0.31 |
fc4 | 0.24 | -0.03 | -0.26 | -0.19 | 0.11 | 0.44 | 0.46 | -0.24 | -0.03 |
fc5 | -0.06 | -0.01 | 0.89 | -0.24 | -0.10 | 0.07 | 0.20 | 0.07 | -0.10 |
fc6 | -0.11 | -0.11 | 0.62 | -0.02 | -0.45 | 0.06 | 0.42 | 0.15 | -0.07 |
fc7 | 0.08 | 0.18 | 0.64 | -0.39 | 0.50 | 0.05 | -0.25 | -0.09 | -0.06 |
fc8 | 0.14 | 0.29 | 0.63 | -0.41 | 0.43 | 0.03 | -0.24 | -0.12 | -0.07 |
fc9 | 0.70 | -0.52 | 0.17 | 0.05 | -0.15 | -0.01 | -0.05 | -0.03 | 0.18 |
fc10 | 0.58 | -0.57 | 0.05 | 0.14 | -0.01 | -0.06 | -0.25 | -0.35 | 0.11 |
fc11 | 0.79 | -0.35 | -0.06 | -0.09 | 0.02 | 0.30 | 0.04 | 0.03 | 0.14 |
fc12 | 0.53 | -0.36 | 0.06 | 0.15 | -0.22 | -0.09 | -0.33 | -0.44 | 0.02 |
表15续
表16采用流式细胞变量的PCA模型的关联矩阵
表16续
fc16 | fc17 | fc18 | fc19 | fc20 | fc21 | fc22 | fc23 | fc24 | fc25 | fc26 | fc27 | fc28 | fc29 | |
fc1 | 0.06 | 0.41 | 0.04 | 0.06 | 0.16 | 0.06 | 0.16 | 0.40 | 0.36 | 0.07 | -0.24 | 0.02 | 0.03 | 0.06 |
fc2 | 0.47 | 0.31 | 0.27 | 0.13 | 0.32 | 0.25 | 0.29 | 0.29 | 0.26 | 0.23 | 0.15 | 0.20 | 0.00 | 0.03 |
fc3 | 0.12 | 0.38 | 0.04 | 0.04 | 0.18 | 0.05 | 0.19 | 0.40 | 0.30 | 0.22 | -0.05 | 0.18 | -0.07 | -0.04 |
fc4 | 0.27 | 0.00 | 0.06 | 0.08 | 0.07 | 0.12 | 0.13 | 0.05 | 0.12 | 0.08 | 0.09 | 0.10 | -0.01 | -0.01 |
fc5 | -0.01 | -0.06 | -0.02 | 0.04 | -0.10 | 0.00 | -0.07 | -0.03 | -0.01 | -0.09 | -0.20 | -0.25 | -0.05 | -0.05 |
fc6 | -0.07 | -0.09 | -0.02 | 0.01 | -0.12 | -0.02 | -0.11 | -0.08 | -0.08 | -0.09 | -0.29 | -0.32 | -0.04 | -0.04 |
fc7 | 0.11 | 0.03 | 0.03 | 0.07 | 0.03 | 0.05 | 0.07 | 0.06 | 0.12 | 0.00 | 0.10 | 0.05 | -0.02 | -0.04 |
fc8 | 0.10 | 0.09 | 0.16 | 0.09 | 0.17 | 0.20 | 0.22 | 0.14 | 0.17 | 0.03 | 0.12 | 0.03 | -0.02 | -0.03 |
fc9 | 0.15 | 0.34 | 0.06 | 0.10 | 0.12 | 0.11 | 0.17 | 0.41 | 0.15 | 0.20 | -0.17 | -0.04 | 0.07 | 0.09 |
fc10 | 0.11 | 0.26 | -0.06 | -0.06 | 0.04 | -0.05 | 0.04 | 0.28 | 0.09 | 0.12 | -0.10 | 0.09 | 0.01 | 0.04 |
fc11 | 0.52 | 0.37 | 0.21 | 0.15 | 0.25 | 0.22 | 0.27 | 0.44 | 0.17 | 0.28 | 0.03 | 0.11 | 0.01 | 0.04 |
fc12 | 0.07 | 0.28 | 0.13 | -0.07 | 0.17 | 0.14 | 0.19 | 0.31 | 0.09 | 0.04 | -0.23 | 0.11 | 0.01 | 0.03 |
fc13 | 0.20 | 0.32 | 0.08 | 0.19 | 0.14 | 0.10 | 0.12 | 0.33 | 0.24 | 0.14 | 0.09 | 0.07 | -0.02 | 0.02 |
fc14 | 0.23 | 0.38 | 0.07 | 0.11 | 0.16 | 0.10 | 0.20 | 0.44 | 0.18 | 0.37 | 0.01 | 0.07 | 0.02 | 0.04 |
fc15 | 0.07 | 0.44 | 0.12 | -0.01 | 0.24 | 0.09 | 0.20 | 0.44 | 0.07 | 0.26 | 0.01 | -0.16 | -0.07 | -0.04 |
fc16 | 1.00 | 0.06 | 0.30 | 0.05 | 0.18 | 0.19 | 0.08 | -0.01 | -0.02 | 0.06 | 0.11 | 0.16 | 0.39 | 0.41 |
fc17 | 0.06 | 1.00 | 0.43 | 0.04 | 0.77 | 0.33 | 0.60 | 0.84 | 0.34 | 0.24 | 0.15 | 0.24 | -0.22 | -0.18 |
fc18 | 0.30 | 0.43 | 1.00 | 0.16 | 0.82 | 0.90 | 0.74 | 0.40 | 0.13 | -0.01 | 0.07 | 0.11 | -0.10 | -0.09 |
fc19 | 0.05 | 0.04 | 0.16 | 1.00 | 0.11 | 0.42 | 0.39 | 0.42 | 0.41 | -0.02 | 0.12 | 0.15 | 0.05 | 0.08 |
fc20 | 0.18 | 0.77 | 0.82 | 0.11 | 1.00 | 0.72 | 0.87 | 0.67 | 0.28 | 0.03 | 0.22 | 0.20 | -0.21 | -0.18 |
fc21 | 0.19 | 0.33 | 0.90 | 0.42 | 0.72 | 1.00 | 0.86 | 0.50 | 0.25 | 0.01 | 0.06 | 0.09 | -0.06 | -0.05 |
fc22 | 0.06 | 0.60 | 0.74 | 0.39 | 0.87 | 0.86 | 1.00 | 0.77 | 0.38 | 0.05 | 0.19 | 0.18 | -0.16 | -0.13 |
fc23 | -0.01 | 0.84 | 0.40 | 0.42 | 0.67 | 0.50 | 0.77 | 1.00 | 0.46 | 0.20 | 0.14 | 0.24 | -0.18 | -0.14 |
fc24 | -0.02 | 0.34 | 0.13 | 0.41 | 0.28 | 0.25 | 0.38 | 0.46 | 1.00 | 0.15 | 0.22 | 0.37 | 0.08 | 0.13 |
fc25 | 0.06 | 0.24 | -0.01 | -0.02 | 0.03 | 0.01 | 0.05 | 0.20 | 0.15 | 1.00 | 0.05 | 0.16 | 0.03 | 0.03 |
fc26 | 0.11 | 0.15 | 0.07 | 0.12 | 0.22 | 0.06 | 0.19 | 0.14 | 0.22 | 0.05 | 1.00 | 0.20 | -0.17 | -0.14 |
fc27 | 0.16 | 0.24 | 0.11 | 0.15 | 0.20 | 0.09 | 0.18 | 0.24 | 0.37 | 0.16 | 0.20 | 1.00 | 0.05 | 0.09 |
fc28 | 0.39 | -0.22 | -0.10 | 0.05 | -0.21 | -0.06 | -0.16 | -0.18 | 0.08 | 0.03 | -0.17 | 0.05 | 1.00 | 0.98 |
fc29 | 0.41 | -0.18 | -0.09 | 0.08 | -0.18 | -0.05 | -0.13 | -0.14 | 0.13 | 0.03 | -0.14 | 0.09 | 0.98 | 1.00 |
表17关联矩阵的特征值,以及对所选择的流式细胞变量的相关统计(fc24-fc29)
PC | 特征值 | 总方差的百分数% | 累积特征值 | 累积方差% |
1 | 2.06 | 34.26 | 2.06 | 34.26 |
2 | 1.60 | 26.64 | 3.65 | 60.90 |
3 | 0.94 | 15.71 | 4.60 | 76.60 |
4 | 0.76 | 12.70 | 5.36 | 89.31 |
5 | 0.63 | 10.45 | 5.99 | 99.76 |
6 | 0.01 | 0.24 | 6.00 | 100.00 |
表18所选择的流式细胞数据对每个PC的负荷(与显示95%CL的特征值>0.5相关)
PC1 | PC2 | PC3 | |
fc24 | -0.23 | -0.72 | -0.13 |
fc25 | -0.10 | -0.40 | 0.87 |
fc26 | 0.20 | -0.60 | -0.39 |
fc27 | -0.19 | -0.72 | -0.06 |
fc28 | -0.98 | 0.14 | -0.05 |
fc29 | -0.98 | 0.09 | -0.08 |
表19关联矩阵的特征值,以及对RT-PCR数据的相关统计
PC | 特征值 | 总方差的百分数% | 累积特征值 | 累积方差% |
1 | 2.9 | 41.0 | 2.9 | 41.0 |
2 | 1.2 | 17.6 | 4.1 | 58.6 |
3 | 1.0 | 14.4 | 5.1 | 73.0 |
4 | 0.6 | 8.3 | 5.7 | 81.3 |
5 | 0.6 | 7.9 | 6.2 | 89.2 |
6 | 0.4 | 6.2 | 6.7 | 95.4 |
7 | 0.3 | 4.6 | 7.0 | 100.0 |
表20 RT-PCR数据对每个PC的负荷(与显示95%CL的特征值>0.5相关)
PC1 | PC2 | PC3 | |
Fas-L | 0.68 | -0.01 | 0.340 |
MCP-1 | 0.71 | -0.38 | -0.27 |
TNF-alpha | 0.77 | 0.21 | -0.25 |
IL-1 | 0.37 | 0.64 | 0.56 |
IL-6 | 0.72 | -0.22 | -0.20 |
II-8 | 0.70 | -0.28 | 0.389 |
II-10 | 0.34 | 0.698 | -0.48* |
●95%CL的临界显著性
表21关联矩阵的特征值,以及对组合的临床、RT-PCR和流式细胞变量的相关统计
特征值 | 总方差的百分数% | 累积特征值 | 累积% | |
1 | 4.550 | 15.167 | 4.550 | 15.167 |
2 | 3.861 | 12.870 | 8.411 | 28.037 |
3 | 3.295 | 10.982 | 11.706 | 39.019 |
4 | 2.028 | 6.760 | 13.734 | 45.779 |
5 | 1.816 | 6.052 | 15.549 | 51.831 |
6 | 1.520 | 5.065 | 17.069 | 56.896 |
7 | 1.234 | 4.113 | 18.303 | 61.010 |
8 | 1.153 | 3.844 | 19.456 | 64.854 |
9 | 1.139 | 3.797 | 20.595 | 68.650 |
10 | 0.958 | 3.192 | 21.553 | 71.842 |
11 | 0.931 | 3.104 | 22.484 | 74.946 |
12 | 0.906 | 3.021 | 23.390 | 77.968 |
13 | 0.825 | 2.748 | 24.215 | 80.716 |
14 | 0.740 | 2.466 | 24.955 | 83.182 |
15 | 0.685 | 2.284 | 25.640 | 85.466 |
16 | 0.611 | 2.036 | 26.251 | 87.502 |
17 | 0.559 | 1.863 | 26.810 | 89.365 |
18 | 0.542 | 1.808 | 27.352 | 91.173 |
19 | 0.480 | 1.601 | 27.832 | 92.774 |
表21续
20 | 0.441 | 1.469 | 28.273 | 94.244 |
21 | 0.371 | 1.235 | 28.644 | 95.479 |
22 | 0.333 | 1.108 | 28.976 | 96.588 |
23 | 0.299 | 0.998 | 29.276 | 97.586 |
24 | 0.250 | 0.833 | 29.526 | 98.419 |
25 | 0.175 | 0.585 | 29.701 | 99.003 |
26 | 0.139 | 0.462 | 29.840 | 99.465 |
27 | 0.114 | 0.381 | 29.954 | 99.846 |
28 | 0.020 | 0.068 | 29.974 | 99.914 |
29 | 0.018 | 0.059 | 29.992 | 99.973 |
30 | 0.008 | 0.027 | 30.000 | 100.000 |
表22组合的临床、RT-PCR和流式细胞数据对每个PC的负荷(与显示95%CL的特征值>0.5相关)
参数 | PC1 | PC2 | PC3 | PC4 | PC5 | PC6 | PC7 | PC8 | PC9 |
CD31(%) | -0.02 | -0.03 | 0.42 | 0.18 | 0.52 | -0.07 | -0.40 | -0.17 | 0.03 |
CD54(%) | -0.16 | 0.04 | 0.20 | -0.09 | 0.48 | 0.48 | 0.06 | 0.18 | -0.19 |
CD62L(%) | -0.41 | -0.09 | 0.20 | 0.28 | 0.16 | -0.21 | 0.01 | -0.01 | 0.39 |
CD11B(%) | -0.13 | -0.16 | 0.22 | 0.19 | 0.54 | -0.13 | -0.09 | -0.19 | 0.05 |
CD69(%) | 0.49 | -0.55 | 0.23 | -0.27 | 0.13 | -0.16 | 0.00 | 0.40 | -0.25 |
CD11B CD69(%) | 0.48 | -0.53 | 0.29 | -0.27 | 0.14 | -0.16 | -0.03 | 0.40 | -0.22 |
体温 | 0.20 | 0.06 | 0.30 | -0.35 | -0.02 | 0.37 | 0.47 | 0.17 | 0.30 |
HR | 0.15 | -0.12 | 0.13 | -0.35 | 0.20 | 0.65 | -0.03 | -0.32 | 0.00 |
MAP | 0.19 | 0.19 | -0.15 | 0.10 | -0.44 | 0.31 | -0.30 | 0.04 | -0.35 |
WCC | 0.57 | -0.64 | 0.14 | 0.38 | 0.01 | 0.06 | 0.02 | -0.08 | 0.08 |
嗜中性粒细胞 | 0.53 | -0.64 | 0.21 | 0.38 | -0.01 | -0.05 | 0.00 | -0.06 | 0.12 |
淋巴细胞 | 0.31 | -0.10 | -0.32 | 0.01 | -0.02 | 0.09 | -0.06 | -0.30 | -0.28 |
单核细胞 | 0.56 | -0.61 | 0.04 | 0.25 | -0.04 | 0.19 | 0.06 | -0.04 | 0.09 |
血小板 | 0.49 | -0.02 | -0.31 | 0.29 | 0.01 | -0.15 | 0.26 | 0.04 | -0.01 |
INR | -0.21 | -0.71 | 0.09 | 0.10 | -0.25 | 0.21 | -0.24 | -0.12 | -0.04 |
APTR | -0.50 | -0.49 | -0.04 | 0.07 | -0.33 | 0.11 | -0.22 | -0.08 | -0.08 |
CRP | -0.07 | -0.09 | 0.26 | -0.63 | 0.14 | -0.32 | -0.23 | -0.18 | 0.02 |
PO2 | -0.23 | 0.13 | -0.10 | 0.18 | 0.45 | -0.10 | 0.14 | -0.12 | -0.49 |
HCO3- | 0.79 | 0.32 | -0.07 | -0.25 | -0.03 | -0.12 | -0.21 | -0.06 | 0.01 |
BXS | 0.80 | 0.33 | -0.07 | -0.24 | -0.04 | -0.11 | -0.22 | -0.08 | 0.03 |
乳酸 | -0.50 | -0.41 | -0.02 | -0.15 | 0.03 | 0.06 | 0.22 | 0.08 | -0.08 |
尿素 | -0.25 | -0.15 | 0.21 | -0.12 | -0.24 | -0.32 | 0.16 | -0.08 | -0.19 |
肌酐 | -0.68 | -0.45 | 0.16 | -0.06 | -0.15 | -0.04 | -0.03 | 0.10 | -0.18 |
Fas-L | -0.05 | 0.12 | 0.77 | -0.08 | 0.01 | 0.11 | -0.16 | -0.05 | -0.05 |
MCP-1 | 0.12 | 0.19 | 0.61 | 0.07 | -0.36 | -0.09 | 0.20 | -0.32 | -0.13 |
TNF-alpha | 0.14 | 0.49 | 0.57 | 0.32 | 0.01 | 0.24 | 0.13 | 0.05 | -0.12 |
IL-1 | -0.22 | 0.28 | 0.20 | 0.20 | -0.12 | 0.11 | -0.45 | 0.49 | 0.16 |
IL-6 | 0.07 | 0.19 | 0.66 | 0.15 | -0.15 | -0.20 | 0.20 | -0.05 | -0.25 |
IL-8 | 0.02 | 0.16 | 0.73 | -0.05 | -0.27 | -0.04 | -0.03 | -0.01 | 0.16 |
IL-10 | 0.14 | 0.44 | 0.02 | 0.51 | 0.09 | 0.03 | -0.05 | 0.21 | -0.19 |
表23所有DFA模型中使用的模型定义
参数 | 值 |
引入变量 | 前向逐步 |
容许度 | 0.010 |
模型中包括的项 | 衍生 |
进入F | 1.00 |
除去F | 0.00 |
步骤数量 | 显著的 |
表24基于含有替代均值的临床数据的PCA得分的DFA模型总结
参数 | 值/项 | |||
Wilks’λ | <0.40 | |||
模型中包括的变量 | PC1,PC3,PC4 | |||
χ2 | 3个显著的DF | |||
累积比例 | DF1 | DF2 | DF3 | |
56.7% | 79.3% | 89.9% | ||
特征值 | 0.60 | 0.24 | 0.11 | |
因子结构 | PC1 | 0.63 | 0.02 | 0.15 |
PC3 | 0.09 | -0.87 | 0.39 | |
PC4 | 0.44 | 0.03 | -0.51 | |
标化等数 | PC1 | 0.78 | 0.02 | 0.13 |
PC3 | 0.12 | -0.89 | 0.35 | |
PC4 | 0.61 | 0.04 | -0.48 | |
典型变量的均值 | 对照 | -0.15 | 0.63 | 0.04 |
第1天 | 0.51 | -0.39 | 0.02 | |
第2天 | 0.57 | -0.17 | -0.64 | |
第3天 | 0.72 | 0.02 | 0.35 | |
第4天 | 0.21 | 0.02 | 0.43 | |
第5天 | -1.80 | -0.01 | -0.28 | |
第6天 | -1.33 | -1.14 | 0.37 |
表25基于含有替代均值的临床数据的PCA得分的判别DFA模型中变量组合的总结
DF | PC | 每个PC的组分 |
1 | 1 | WCC,单核细胞,嗜中性粒细胞和血小板,BXS和HCO3 -与肌酐和乳酸负相关 |
4 | CRP | |
2 | 3 | HR和CRP |
3 | 4 | CRP |
表26基于含有替代均值的临床数据的PCA得分的DFA模型的分类矩阵
组 | 正确百分数 | 对照 | 1天 | 2天 | 3天 | 4天 | 5天 | 6天 |
p=.28 | p=23 | p=.14 | p=.08 | p=.11 | p=.09 | p=.05 | ||
对照 | 80 | 24 | 4 | 0 | 1 | 0 | 0 | 1 |
1天 | 36 | 12 | 9 | 3 | 0 | 0 | 0 | 1 |
2天 | 40 | 5 | 2 | 6 | 0 | 1 | 0 | 1 |
3天 | 11 | 1 | 5 | 1 | 1 | 1 | 0 | 0 |
4天 | 8 | 7 | 3 | 0 | 0 | 1 | 1 | 0 |
5天 | 50 | 2 | 2 | 0 | 0 | 0 | 5 | 1 |
6天 | 83 | 1 | 0 | 0 | 0 | 0 | 0 | 5 |
总计 | 48 | 52 | 25 | 10 | 2 | 3 | 6 | 9 |
表27基于不含替代均值的临床数据的PCA得分的DFA模型的总结
参数 | 值/项 | |||
Wliks’λ | <0.45 | |||
模型中包含的变量 | PC1,PC3,PC5 | |||
χ2 | 3个显著的DF | |||
累积比例 | DF1 | DF2 | DF3 | |
52% | 89% | 95% | ||
特征值 | 0.49 | 0.35 | 0.05 | |
因子结构 | PC1 | 0.65 | 0.09 | 0.08 |
PC3 | -0.13 | 0.80 | 0.45 | |
PC5 | -0.61 | -0.32 | -0.07 | |
标化系数 | PC1 | 0.79 | 0.10 | 0.06 |
PC3 | -0.15 | 0.87 | 0.38 | |
PC5 | -0.73 | -0.34 | -0.06 | |
典型变量的均值 | 对照 | -0.18 | 0.56 | -0.10 |
第1天 | 0.54 | -0.59 | 0.07 | |
第2天 | 0.45 | -0.14 | 0.36 | |
第3天 | 0.81 | 0.40 | -0.19 | |
第4天 | 0.22 | -0.05 | -0.27 | |
第5天 | -1.61 | -0.00 | 0.33 | |
第6天 | -1.40 | -1.90 | -0.59 |
表28基于不含替代均值的临床数据的PCA得分的判别DFA模型中变量组合的总结
DF | PC | 每个PC的组分 |
1 | 1 | WCC,单核细胞,嗜中性粒细胞和血小板,BXS和HCO3 -与肌酐和乳酸负相关 |
4 | PO2与尿素和MAP相反 | |
2 | 3 | HR和CRP |
3 | 4 | CRP |
表29基于不含替代均值的临床数据的PCA得分的DFA模型的分类矩阵
组 | 正确百分率 | 对照 | 1天 | 2天 | 3天 | 4天 | 5天 | 6天 |
p=.34 | p=.21 | p=.13 | p=.06 | p=.10 | p=.09 | p=.03 | ||
对照 | 83 | 25 | 2 | 1 | 1 | 0 | 1 | 0 |
1天 | 47 | 9 | 9 | 1 | 0 | 0 | 0 | 0 |
2天 | 0 | 6 | 6 | 0 | 0 | 0 | 0 | 0 |
3天 | 0 | 3 | 3 | 0 | 0 | 0 | 0 | 0 |
4天 | 0 | 5 | 3 | 0 | 0 | 0 | 1 | 0 |
5天 | 25 | 4 | 2 | 0 | 0 | 0 | 2 | 0 |
6天 | 67 | 0 | 0 | 0 | 0 | 0 | 1 | 2 |
总计 | 44 | 52 | 25 | 2 | 1 | 0 | 5 | 2 |
表30基于含有替代均值的经转换的临床数据的DFA模型的总结
参数 | 值/项 | |||||
Wilks’λ | 0.25 | |||||
模型中包括的变量 | BXS,CRP,lactate,urea,temperature,creatinine,neuts,PO2,HCO3- | |||||
χ2 | 5个显著DF | |||||
累积比例 | DF1 | DF2 | DF3 | DF4 | DF5 | |
46% | 74% | 85% | 94% | 99% | ||
特征值 | 0.78 | 0.47 | 0.18 | 0.15 | 0.09 | |
因子结构 | BXS | -0.62 | 0.58 | -0.09 | 0.21 | 0.22 |
CRP | 0.12 | 0.63 | 0.16 | -0.03 | -0.31 | |
体温 | 0.17 | 0.03 | -0.56 | 0.62 | 0.11 | |
乳酸 | 0.65 | 0.03 | 0.05 | 0.05 | 0.01 | |
尿素 | -0.19 | -0.18 | -0.19 | 0.17 | -0.57 | |
肌酐 | 0.42 | -0.26 | 0.51 | 0.20 | -0.19 | |
HCO3- | -0.54 | 0.67 | -0.13 | 0.19 | 0.16 | |
PO2 | 0.16 | -0.05 | -0.04 | -0.67 | -0.31 | |
嗜中性粒细胞 | -0.14 | 0.05 | -0.24 | -0.18 | 0.70 | |
标化系数 | BXS | -1.45 | -0.78 | 0.84 | 0.09 | 1.00 |
CRP | 0.13 | 0.65 | 0.07 | -0.20 | -0.20 | |
体温 | 0.41 | -0.02 | -0.64 | 0.53 | 0.02 | |
乳酸 | 0.56 | 0.34 | -0.22 | -0.01 | 0.21 | |
尿素 | -0.52 | -0.22 | -0.51 | 0.09 | -0.60 | |
肌酐 | 0.13 | -0.11 | 0.93 | 0.57 | 0.10 | |
HCO3- | 1.02 | 1.45 | -0.53 | 0.32 | -0.88 | |
PO2 | 0.19 | 0.07 | -0.28 | -0.61 | -0.21 | |
嗜中性粒细胞 | -0.10 | 0.01 | -0.30 | -0.32 | 0.67 | |
典型变量的均值 | 对照 | -0.54 | -0.90 | -0.10 | -0.14 | -0.07 |
第1天 | -0.22 | 0.53 | -0.35 | -0.03 | 0.40 | |
第2天 | -0.19 | 0.84 | 0.20 | -0.58 | -0.36 | |
第3天 | -0.42 | 0.64 | 0.41 | 0.62 | -0.28 | |
第4天 | -0.34 | -0.09 | 0.43 | 0.61 | 0.06 | |
第5天 | 2.05 | -0.37 | 0.61 | -0.23 | 0.27 | |
第6天 | 1.95 | 0.03 | -1.05 | 0.47 | -0.54 |
表31基于含有替代均值的临床数据的DFA模型中变量组合的总
结
DF | 确定DF的临床变量 |
1 | 与乳酸呈负相关的BXS和HCO3 - |
2 | BXS,CRP和HCO3 -之间相关 |
3 | 体温和肌酐呈负相关 |
4 | 体温和PO2呈负相关 |
5 | 尿素和嗜中性粒细胞呈负相关 |
表32基于含有替代均值的临床数据的DFA模型的分类矩阵
组 | 正确百分数 | 对照 | 1天 | 2天 | 3天 | 4天 | 5天 | 6天 |
p=.28 | p=.23 | p=.14 | p=.08 | p=.11 | p=.09 | p=.05 | ||
对照 | 80 | 24 | 5 | 1 | 0 | 0 | 0 | 0 |
1天 | 56 | 4 | 14 | 5 | 0 | 1 | 0 | 1 |
2天 | 60 | 2 | 2 | 9 | 1 | 0 | 0 | 1 |
3天 | 11 | 1 | 4 | 2 | 1 | 1 | 0 | 0 |
4天 | 8 | 7 | 4 | 0 | 0 | 1 | 0 | 0 |
5天 | 50 | 1 | 2 | 1 | 0 | 0 | 5 | 1 |
6天 | 83 | 0 | 0 | 0 | 0 | 0 | 1 | 5 |
总计 | 55 | 39 | 31 | 18 | 2 | 3 | 6 | 8 |
表33基于含有替代均值的流式细胞数据的PCA得分的DFA模型的总结
参数 | 值/项 | ||
Wilks’λ | <0.39 | ||
模型中包括变量 | PC1,PC5,PC8 | ||
χ2 | 2个显著的DF | ||
累积比例 | DF1 | DF2 | |
46% | 71% | ||
特征值 | 0.49 | 0.26 | |
因子结构 | PC1 | 0.30 | -0.04 |
PC5 | -0.19- | -0.59 | |
PC8 | 0.419 | -0.14 | |
标化系数 | PC1 | 0.40 | -0.04 |
PC5 | -0.25 | -0.67 | |
PC8 | 0.56 | -0.16 | |
典型变量的均值 | 对照 | -0.21 | -0.35 |
第1天 | -0.71 | -0.12 | |
第2天 | -0.24 | 0.16 | |
第3天 | 0.14 | 1.44 | |
第4天 | 0.48 | 0.32 | |
第5天 | 0.88 | -0.28 | |
第6天 | 1.99 | -0.43 |
表34基于含有替代均值的流式细胞数据的PCA得分的DFA模型中变量组合的总结
DF | PC | 每个DF的组分 |
1 | 1 | fc1-3,9-15,17,20-23 |
8 | fc19 | |
2 | 5 | fc7+26 |
表35基于含有替代均值的流式细胞数据的PCA得分的DFA模型的分类矩阵
组 | 正确百分数 | 对照 | 1天 | 2天 | 3天 | 4天 | 5天 | 6天 |
p=.28 | p=.23 | p=.14 | p=.08 | p=.11 | p=.09 | p=.05 | ||
对照 | 66.6 | 20 | 7 | 0 | 1 | 0 | 1 | 1 |
1天 | 36.0 | 12 | 9 | 2 | 0 | 0 | 2 | 0 |
2天 | 13.3 | 4 | 8 | 2 | 1 | 0 | 0 | 0 |
3天 | 44.4 | 2 | 2 | 0 | 4 | 0 | 1 | 0 |
4天 | 33.3 | 3 | 2 | 1 | 0 | 4 | 1 | 1 |
5天 | 40.0 | 4 | 1 | 0 | 0 | 0 | 4 | 1 |
6天 | 66.6 | 1 | 1 | 0 | 0 | 0 | 0 | 4 |
总计 | 43.9 | 46 | 30 | 5 | 6 | 4 | 9 | 7 |
表36基于含有替代均值的流式细胞数据的DFA模型的总结
参数 | 值/项 | |||
Wilks’λ | <0.034 | |||
模型中包括的变量 | fc7-8,fc11,fc16,fc25, fc28,fc29 | |||
χ2 | 3个显著的DF | |||
累积比例 | DF1 | DF2 | DF3 | |
46% | 62% | 74% | ||
特征值 | 2.17 | 0.94 | 0.63 | |
因子结构 | fc7 | 0.06 | -0.16 | -0.10 |
fc8 | -0.06 | -0.19 | -0.07 | |
fc11 | 0.02 | 0.13 | -0.24 | |
fc16 | 0.05 | -0.14 | 0.14 | |
fc25 | -0.34 | -0.23 | -0.08 | |
fc28 | 0.04 | -0.02 | 0.02 | |
fc29 | 0.06 | -0.03 | 0.01 | |
标化系数 | fc7 | 1.84 | 0.20 | -1.03 |
fc8 | -1.66 | -1.00 | 0.32 | |
fc11 | 0.87 | 1.89 | -1.16 | |
fc16 | -0.73 | -1.33 | 0.80 | |
fc25 | -0.86 | -0.29 | 0.02 | |
fc28 | -2.34 | 0.65 | 0.29 | |
fc29 | 2.64 | 0.12 | -0.66 | |
典型变量的均值 | 对照 | 0.50 | -0.92 | -0.40 |
第1天 | 1.07 | -0.32 | 0.12 | |
第2天 | 0.22 | 0.07 | 1.10 | |
第3天 | 0.20 | 1.75 | 1.26 | |
第4天 | -0.39 | 1.49 | -1.34 | |
第5天 | -1.05 | 0.73 | -0.55 | |
第6天 | -5.27 | -1.06 | 0.48 |
表37基于含有替代均值的流式细胞数据的DFA模型中变量组合的总结
DF | 限定DF的临床变量 |
1 | fc7,16,28和29相关,这些与fc8和25相反 |
2 | fc7,8,16和25相关,这些与fc12相反 |
3 | 与fc11相关 |
表38基于含有替代均值的临床数据的DFA模型的分类矩阵
组 | 正确百分数 | 对照 | 1天 | 2天 | 3天 | 4天 | 5天 | 6天 |
p=.28 | p=.23 | p=.14 | p=.08 | p=.11 | p=.09 | p=.05 | ||
对照 | 80.0 | 24.0 | 2.0 | 0.0 | 1.0 | 2.0 | 1.0 | 0.0 |
1天 | 68.0 | 6.0 | 17.0 | 1.0 | 0.0 | 1.0 | 0.0 | 0.0 |
2天 | 73.3 | 3.0 | 1.0 | 11.0 | 0.0 | 0.0 | 0.0 | 0.0 |
3天 | 66.7 | 2.0 | 0.0 | 0.0 | 6.0 | 1.0 | 0.0 | 0.0 |
4天 | 91.7 | 1.0 | 0.0 | 0.0 | 0.0 | 11.0 | 0.0 | 0.0 |
5天 | 80.0 | 0.0 | 1.0 | 1.0 | 0.0 | 0.0 | 8.0 | 0.0 |
6天 | 83.3 | 0.0 | 0.0 | 0.0 | 0.0 | 1.0 | 0.0 | 5.0 |
总计 | 76.6 | 36.0 | 21.0 | 13.0 | 7.0 | 16.0 | 9.0 | 5.0 |
表39基于组合的临床、RT-PCR和流式细胞变量的DFA模型的总结
参数 | 值/项 | ||||
Wilks’λ | <0.0087 | ||||
模型中包括的变量 | fc7-8,fc11,fc16,fc25,fc28,fc29 | ||||
χ2 | 3个显著的DF | ||||
累积比例 | DF1 | DF2 | DF3 | DF4 | |
44% | 63% | 79% | 89% | ||
特征值 | 3.79 | 1.60 | 1.34 | 0.86 | |
因子结构 | BXS | -0.30 | 0.08 | -0.26 | 0.17 |
fc25 | 0.22 | -0.23 | -0.08 | 0.13 | |
fc22 | 0.22 | -0.03 | 0.08 | 0.28 | |
fc11 | -0.02 | 0.12 | 0.05 | 0.23 | |
体温 | 0.09 | -0.11 | -0.13 | 0.20 | |
CRP | 0.01 | 0.16 | -0.33 | -0.16 | |
fc18 | 0.15 | 0.01 | -0.12 | 0.06 | |
fc6 | -0.03 | 0.01 | -0.02 | 0.I2 | |
IL-6 | -0.11 | -0.06 | 0.11 | -0.07 | |
INR | 0.00 | 0.06 | 0.03 | 0.12 | |
APTR | 0.15 | 0.12 | 0.30 | -0.03 | |
fc16 | -0.07 | -0.09 | -0.11 | -0.07 | |
尿素 | -0.02 | -0.10 | 0.01 | 0.17 | |
乳酸 | 0.25 | 0.04 | -0.05 | -0.22 | |
Fas-L | 0.00 | 0.15 | -0.03 | 0.04 | |
fc13 | -0.04 | -0.01 | -0.07 | 0.11 | |
fc24 | 0.00 | -0.15 | 0.11 | 0.15 | |
fc1 | -0.01 | 0.12 | -0.12 | 0.22 | |
fc3 | -0.02 | 0.04 | -0.02 | 0.15 | |
MCP-1 | -0.06 | 0.04 | 0.02 | -0.02 | |
fc28 | -0.05 | -0.03 | -0.08 | -0.01 | |
Il-10 | 0.01 | -0.10 | 0.09 | -0.07 |
表39续.基于组合的临床、RT-PCR和流式细胞变量的DFA模型的总结
fc27 | 0.01 | -0.12 | 0.08 | -0.02 | |
fc26 | 0.04 | -0.10 | 0.19 | -0.06 | |
嗜中性粒细胞 | -0.11 | -0.07 | 0.03 | -0.04 | |
fc14 | 0.00 | 0.09 | 0.08 | 0.19 | |
WCC | -0.10 | -0.08 | 0.03 | -0.01 | |
fc29 | -0.07 | -0.04 | -0.08 | 0.03 | |
血小板 | -0.11 | 0.04 | 0.12 | -0.03 | |
PO2 | 0.05 | -0.07 | 0.03 | -0.27 | |
参数 | 值/项 | ||||
累积比例 | DF1 | DF2 | DF3 | DF4 | |
44% | 63% | 79% | 89% | ||
特征值 | 3.79 | 1.60 | 1.34 | 0.86 | |
标化系数 | BXS | -0.80 | 0.33 | -0.21 | -0.09 |
fc25 | 0.81 | -0.52 | -0.29 | 0.22 | |
fc22 | 0.46 | -0.58 | 0.05 | 0.68 | |
fc11 | 0.04 | 2.00 | 0.53 | 0.96 | |
体温 | 0.51 | -0.41 | -0.07 | 0.23 | |
CRP | 0.07 | 0.14 | -0.58 | -0.13 | |
fc18 | 0.64 | 0.72 | -0.53 | -0.11 | |
fc6 | 0.02 | -0.29 | 0.17 | 0.08 | |
IL-6 | -0.69 | -0.38 | 0.28 | -0.12 | |
INR | -0.90 | 0.03 | -0.99 | 0.94 | |
APTR | 0.82 | 0.62 | 1.28 | -0.44 | |
fc16 | -0.63 | -1.23 | -0.54 | -0.69 | |
尿素 | -0.13 | -0.38 | 0.06 | 0.58 | |
乳酸 | 0.14 | 0.12 | -0.06 | -0.68 | |
Fas-L | 0.20 | 0.92 | -0.04 | 0.05 | |
fc13 | -0.27 | -0.82 | -0.33 | -0.07 | |
fc24 | -0.46 | -0.72 | 0.11 | 0.15 | |
fc1 | 1.04 | 0.84 | 0.09 | 0.12 | |
fc3 | -0.67 | -0.46 | -0.41 | 0.06 | |
MCP-1 | 0.33 | 0.01 | -0.21 | -0.05 | |
fc28 | 3.20 | 0.89 | -2.48 | -0.21 | |
ll-10 | 0.22 | -0.03 | 0.26 | -0.16 | |
fc27 | 0.12 | 0.14 | 0.23 | -0.13 | |
fc26 | 0.15 | 0.26 | 0.27 | -0.01 | |
嗜中性粒细胞 | 0.83 | 0.74 | 0.27 | -1.34 | |
fc14 | -0.43 | -0.73 | 0.32 | -0.88 | |
WCC | -0.70 | -1.06 | 0.05 | 0.77 | |
fc29 | -2.73 | -0.34 | 2.71 | 0.56 | |
血小板 | -0.47 | -0.06 | 0.33 | 0.27 | |
PO2 | -0.17 | -0.12 | 0.20 | -0.45 | |
典型变量的均值 | 对照 | -0.39 | -1.04 | 1.24 | 0.23 |
第1天 | -1.80 | -0.51 | -0.50 | -0.19 | |
第2天 | -0.65 | 0.70 | -1.16 | -1.17 | |
第3天 | -0.40 | 1.12 | -1.19 | 0.45 | |
第4天 | 0.38 | 2.02 | 0.20 | 1.77 | |
第5天 | 3.33 | 1.36 | 1.41 | -1.45 | |
第6天 | 5.37 | -2.44 | -2.21 | 0.76 |
表40基于组合的临床、RT-PCR和流式细胞变量的DFA模型的
分类矩阵
组 | 正确百分数 | 对照 | 1天 | 2天 | 3天 | 4天 | 5天 | 6天 |
p=.28 | p=.23 | p=.14 | p=.08 | p=.11 | p=.09 | p=.05 | ||
对照 | 90.0 | 27.0 | 0.0 | 1.0 | 1.0 | 1.0 | 0.0 | 0.0 |
1天 | 76.0 | 3.0 | 19.0 | 3.0 | 0.0 | 0.0 | 0.0 | 0.0 |
2天 | 93.3 | 0.0 | 1.0 | 14.0 | 0.0 | 0.0 | 0.0 | 0.0 |
3天 | 77.8 | 1.0 | 0.0 | 1.0 | 7.0 | 0.0 | 0.0 | 0.0 |
4天 | 91.7 | 0.0 | 0.0 | 1.0 | 0.0 | 11.0 | 0.0 | 0.0 |
5天 | 90.0 | 0.0 | 0.0 | 1.0 | 0.0 | 0.0 | 9.0 | 0.0 |
6天 | 100.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 6.0 |
总计 | 86.9 | 31.0 | 20.0 | 21.0 | 8.0 | 12.0 | 9.0 | 6.0 |
例7:二元logistic回归分析以预测败血症
采用二元logistic回归模型以单独分析来自发展为败血症和呈现阳性微生物结果的ICU患者的RT-PCR、流式细胞结果和临床数据。该模型采用由年龄匹配的、没有诊断为败血症的作为对照组的ICU患者获得的结果。尽管该模型确定了大量的可能预测物,但是某些预测物似乎具有有限的用途,其原因在于,对有症状之前的败血症患者所获得的值在非败血症患者所获得的值的范围内。表36中列出了确实产生某些不同于非败血症数据的败血症前数据点的潜在预测标记物。然而,一旦组合,这些预测标记物仅仅从24个前败血症患者中确定8个。
表41通过二元logistic回归分析确定的潜在预测标记物的总结
时间点 | 总的预测效果 | 预测物 | 概率得分 |
诊断前5天 | p=0.10 | 碱过剩 | p=0.047 |
诊断前4天 | p=0.042 | IL-10 | p=0.058 |
诊断前3天 | p=0.095 | 血液碳酸氢盐 | p=0.021 |
p=0.114 | CD 31 | p=0.037 | |
p=0.055 | IL-10 | p=0.052 | |
诊断前2天 | p=0.10 | 血液碳酸氢盐 | p=0.017 |
C反应蛋白 | p=0.038 | ||
P=0.085 | TNF-α | p=0.022 |
发现可以从24名稍后发展成为具有明确感染的SIRS的患者中预测出8个败血症的标记物组合并没有构成诊断检测。尽管粒细胞上的CD31的预测能力看上出有希望(66%),但是该标记物仅仅在表现出临床症状前3天有效。仅仅基于CD31的检验不能构成诊断检测,其原因是,如果有效,需要更大的诊断窗口。其可以通过发现更多的标记物来实现。然而,该研究可能已经发现,某些标记物可以在以后构成诊断检测的一部分,但是应当谨慎。在20世纪八十年代中期,HLA-DR被认为可以预测感染和败血症的发生(Spilltler,A.&Roth,E.2003,Intensive Care Med,第29卷,第1211-1213页)。然而,近来的研究显示,该标记物的手术后水平不能预测SIRS、败血症或感染并发症的发生(Oczenski,W等2003,Intensive Care Med,第29卷,第1253-1257页,以及Perry,S等2003,Intensive Care Med,第29卷,第1245-1252页)。就HLA-DR而出现的所述相矛盾的报告显示为什么应谨慎使用本研究的结果。这些发现可能是因为区域因素,例如抗生素政策、诊断标准、临床实践、外科方法、治疗方案、环境因素以及患者的嗜好因素(Angus,D.&Wax,R.2001,Critical Care Medicine,第29卷,第7期(增刊),第109-116页)。需要进行涉及来自几个不同医院、优选分布在不同卫生机构的更多病人的更大研究,以进一步评价确定用于预测败血症的标记物的用处。
例8:败血症作为对生物武器反应的模型
背景和方法
由于所要求的发明的一个应用是早期检测由暴露于生物武器所导致的人为感染,因此,调查了败血症作为所述感染的模型的应用。预测假定的生物武器病原例如Burkholderia pseudomallei和弗朗西斯氏菌属tularensis产生严重的败血症(见表42),其出于显而易见的原因而难以建模。
然而,体外全血感染可以作为模型,并测量活化标记物表达和细胞因子应答。为了将该体外感染模型与体内的情况相比,选择金黄色葡萄球菌感染作为感染物模型,其可以直接与体内医源性感染数据相比。
表42
感染 | “不良”结果的形式 |
炭疽 | 败血症,脓毒性休克 |
布鲁氏菌病 | 败血症,脓毒性休克或慢性形式 |
鼻疽和类鼻疽 | 败血症,脓毒性休克或持续形式 |
鼠疫 | 败血症,脓毒性休克 |
兔热病 | 败血症,脓毒性休克 |
流行性斑疹伤寒 | 败血症 |
Q热 | 败血症 |
埃博拉和马尔堡出血热 | 败血症,脓毒性(中毒性感染,血量减少)休克 |
日本脑炎 | 败血症,脓毒性休克 |
天花 | 败血症,脓毒性休克 |
黄热病 | 败血症 |
来自25名健康志愿者的血液在体外用金黄色葡萄球菌感染,在感染后24小时和48小时如前所述测量以下活化标记物和细胞因子水平。
FACS
树突状细胞:CD54,CD97,CCR6,CCR7
NK细胞:CD25,CD44,CD62L,CD69,CD97
单核细胞:CD44,CD54,CD62L,CD69,CD97,CD107a
嗜中性粒细胞:CD44,CD62L,CD69,CD107a
实时RT PCR
IL-1β,IL-6,IL-8,IL-10,MCP-1,TNF-α,sFasL
将这些组的输入参数的每一个(即树突状细胞标记物、NK细胞标记物、单核细胞标记物、嗜中性粒细胞标记物、24小时的RT PCR数据、48小时的RT PCR数据)用于训练其自身的神经网络模型。采用随机选择的感染或非感染血样以用于训练(70%)或接下来的检验(30%)。神经网络分析的检验期基于其正确预测输入参数检验组来自感染或非感染样本的次数百分数而给出预测准确度。每组输入参数的检验重复5次。每次进行检验时,采用新的随机选择的70%感染和非感染样品构建新的神经网络。根据由5个随机输入数据组构建的5个神经网络计算得到的5个预测准确度,计算其均值,从而对每组输入参数,得到平均预测准确度。该方法类似于在败血症患者研究中所采用的。
结果
由RT PCR结果得到最一致的结果。附图4显示了从三个受试者得到的数据,其显示所述轮廓变化的某些不一样的模式。然而,当进行上述的神经网络分析时,所述算法得到良好的感染样品对未感染对照的识别水平(附图5)。
Claims (19)
1.一种筛选生物样品以检测感染、SIRS或败血症的早期阶段的方法,包括以下步骤:
a)通过RT-PCR检测一组提供信息的标记物的表达,和/或通过流式细胞术检测一组提供信息的细胞表面标记物的表达,和/或监测一组标准临床测量结果
b)通过神经网络或多变量统计分析来分析检测的结果
c)根据明显感染的发生几率和时间对所述样品进行分类。
2.如权利要求1所述的方法,其中结果的分析产生对临床SIRS或败血症发生概率的预测。
3.如权利要求1所述的方法,其中结果的分析产生对临床SIRS或败血症发生的二元是/否预测。
4.如权利要求2或3所述的方法,其中如果预测了临床SIRS或败血症的发生,则对结果进行第二次分析以确定所述临床疾病的可能时间和/或严重度。
5.如权利要求1至4任意所述的方法,其中通过RT-PCR检测一组提供信息的生物标记物的表达,所述生物标记物组由至少6个选自以下组的标记物组成,所述组由FasL、MCP-1、TNFα、IL-1β、IL-6、IL-8、IL-10、INF-α和INF-γ组成。
6.如权利要求5所述的方法,其中检测至少7个生物标记物的表达。
7.如权利要求1至4任意所述的方法,其中通过流式细胞术检测一组提供信息的生物标记物的表达,所述生物标记物组由至少2个选自以下组的标记物组成,所述组由CD31、HLA-DR、CD54、CD11b、CD62L、CD25、CD69、CD80和CD97组成。
8.如权利要求1至4任意所述的方法,其中被分析的标准临床测量结果选自由以下构成的组,温度、心率、白细胞总数和白细胞分类计数(单核细胞、淋巴细胞、粒细胞、嗜中性粒细胞)、血小板计数、血清肌酐、尿素、乳酸、碱过剩、pO2、HCO3 -和C-反应蛋白。
9.如权利要求1至8任意所述的方法,其中通过多变量统计分析进行的分析包括主成分分析和/或判别函数分析。
10.如权利要求9所述的方法,其中多变量统计分析包括判别函数分析。
11.如权利要求1至8任意所述的方法,其中通过神经网络进行分析。
12.如权利要求10所述的方法,其中神经网络为多层感知神经网络。
13.如权利要求9至11任意所述的方法,其中该分析能够正确预测大于80%病例的临床SIRS或败血症。
14.一种筛选生物样品以检测感染、SIRS或败血症的早期阶段的系统,其包括:从所述样品获得的细胞中提取和纯化RNA的装置、热循环控制装置或其它通过逆转录聚合酶链反应(RT-PCR)放大所选择的RNA序列的装置、检测并定量所述RT-PCR的结果的装置,能够分析所述结果的经训练的基于计算机的神经网络和由此能够向操作者交流神经网络分析结论的显示装置。
15.一种筛选生物样品以检测SIRS或败血症的早期阶段的系统,其包括:标记特定细胞表面标记物和对所述标记物的表达进行定量的装置、检测和量化所述标记和定量结果的装置、能够分析所述结果的经训练的基于计算机的神经网络和通过其可以将神经网络分析的结论传递给操作者的显示装置,优选地,所述标记通过标记的抗体或抗体碎片记进行,量化是通过荧光激活细胞分类(FACS)或其它形式的流式细胞术进行。
16.根据权利要求13或14所述的系统,包括同时进行RT-PCR和FACS分析的装置。
17.根据权利要求13至15任意所述的系统,进一步包括从标准临床测量结果获取和分析数据的装置。
18.根据权利要求1至12任意所述方法或采用权利要求13至16任意所述系统进行的分析,以制备用于诊断感染、SIRS或败血症的诊断装置。
19.权利要求1至12任意所述的方法或权利要求13至16任意所述系统的应用,以制备用于诊断感染、SIRS或败血症的诊断装置。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB0426982.5 | 2004-12-09 | ||
GBGB0426982.5A GB0426982D0 (en) | 2004-12-09 | 2004-12-09 | Early detection of sepsis |
PCT/GB2005/004755 WO2006061644A1 (en) | 2004-12-09 | 2005-12-09 | Early detection of sepsis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101076806A true CN101076806A (zh) | 2007-11-21 |
CN101076806B CN101076806B (zh) | 2012-01-11 |
Family
ID=34073429
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2005800424703A Expired - Fee Related CN101076806B (zh) | 2004-12-09 | 2005-12-09 | 败血症的早期检测 |
Country Status (8)
Country | Link |
---|---|
US (1) | US20080114576A1 (zh) |
EP (2) | EP1839218A1 (zh) |
JP (1) | JP2008522604A (zh) |
CN (1) | CN101076806B (zh) |
AU (1) | AU2005313114B2 (zh) |
CA (1) | CA2587270A1 (zh) |
GB (1) | GB0426982D0 (zh) |
WO (1) | WO2006061644A1 (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103649329A (zh) * | 2010-11-26 | 2014-03-19 | ImmuneXpress有限公司 | 诊断和/或筛选剂及其用途 |
CN108431832A (zh) * | 2015-12-10 | 2018-08-21 | 渊慧科技有限公司 | 利用外部存储器扩增神经网络 |
CN108604466A (zh) * | 2016-01-28 | 2018-09-28 | 贝克曼考尔特公司 | 感染检测和区分系统及方法 |
CN111374639A (zh) * | 2018-12-28 | 2020-07-07 | 陈柏龄 | 败血症的菌种预测系统与方法 |
CN111989032A (zh) * | 2018-06-15 | 2020-11-24 | 拜克门寇尔特公司 | 存在母细胞标志时确定败血症的方法 |
CN112447289A (zh) * | 2019-08-30 | 2021-03-05 | 希尔-罗姆服务公司 | 败血症监测系统 |
WO2021102776A1 (zh) * | 2019-11-28 | 2021-06-03 | 北京肿瘤医院(北京大学肿瘤医院) | 一种用于早期体腔内感染性并发症诊断的标记物及方法 |
CN114167062A (zh) * | 2020-09-10 | 2022-03-11 | 北京肿瘤医院(北京大学肿瘤医院) | 一种用于早期腹腔内感染性并发症诊断的标记物及方法 |
TWI803893B (zh) * | 2021-06-28 | 2023-06-01 | 中國醫藥大學 | 敗血症之人工智慧輔助醫療診斷方法及其系統 |
Families Citing this family (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9053222B2 (en) | 2002-05-17 | 2015-06-09 | Lawrence A. Lynn | Patient safety processor |
US8068990B2 (en) * | 2003-03-25 | 2011-11-29 | Hologic, Inc. | Diagnosis of intra-uterine infection by proteomic analysis of cervical-vaginal fluids |
US7191068B2 (en) * | 2003-03-25 | 2007-03-13 | Proteogenix, Inc. | Proteomic analysis of biological fluids |
AU2006236588A1 (en) | 2005-04-15 | 2006-10-26 | Becton, Dickinson And Company | Diagnosis of sepsis |
EP1951899B1 (en) * | 2005-11-25 | 2010-06-30 | The Provost, Fellows And Scholars Of The College Of The Holy And Undivided Trinity Of Queen Elizabeth Near Dublin | A method for detecting or monitoring sepsis by analysing cytokine mrna expression levels |
US8295392B2 (en) | 2005-12-14 | 2012-10-23 | Nec Corporation | Digital communication system, indoor unit, and outdoor unit |
GB0610078D0 (en) * | 2006-05-20 | 2006-06-28 | Secr Defence | Sepsis detection microarray |
EP1950310A1 (en) * | 2007-01-23 | 2008-07-30 | Charite-Universitätsmedizin Berlin | Method for risk prediction of a postoperative sepsis in a human |
GB0722582D0 (en) * | 2007-11-16 | 2007-12-27 | Secr Defence | Early detection of sepsis |
US8669113B2 (en) | 2008-04-03 | 2014-03-11 | Becton, Dickinson And Company | Advanced detection of sepsis |
WO2009137682A1 (en) | 2008-05-07 | 2009-11-12 | Lynn Lawrence A | Medical failure pattern search engine |
US11728041B2 (en) * | 2008-05-07 | 2023-08-15 | Lawrence A. Lynn | Real-time time series matrix pathophysiologic pattern processor and quality assessment method |
EP2166358A1 (en) * | 2008-09-17 | 2010-03-24 | Fundacio Institut de Recerca de l'Hospital Universitari Vall d'Hebron | Differential diagnostic biomarkers of stroke mimicking conditions and methods of use thereof |
US10460408B2 (en) * | 2009-01-31 | 2019-10-29 | Mayo Foundation For Medical Education And Research | Presentation of critical patient data |
JP5667353B2 (ja) * | 2009-09-25 | 2015-02-12 | シスメックス株式会社 | 血球計数装置、診断支援装置、診断支援方法及びコンピュータプログラム |
CA2788636C (en) * | 2010-03-02 | 2020-07-07 | F. Hoffmann-La Roche Ag | Il-6 detection based early diagnosis and prediction of systemic inflammatory response syndrome and sepsis in asymptomatic patients |
US9719842B2 (en) * | 2010-07-06 | 2017-08-01 | Norman A. Paradis | Method for the discovery, validation and clinical application of multiplex biomarker algorithms based on optical, physical and/or electromagnetic patterns |
US20140120174A1 (en) * | 2012-10-30 | 2014-05-01 | Lascco, Sa | Methods of prognosis and diagnosis of sepsis |
US10354429B2 (en) | 2012-11-14 | 2019-07-16 | Lawrence A. Lynn | Patient storm tracker and visualization processor |
US9953453B2 (en) | 2012-11-14 | 2018-04-24 | Lawrence A. Lynn | System for converting biologic particle density data into dynamic images |
WO2014117165A1 (en) * | 2013-01-28 | 2014-07-31 | Vanderbilt University | Method for differentiating sepsis and systemic inflammatory response syndrome (sirs) |
WO2014134526A1 (en) | 2013-02-28 | 2014-09-04 | Lynn Lawrence A | System for analysis and imaging using perturbation feature quanta |
KR101595784B1 (ko) * | 2013-05-31 | 2016-02-26 | 가톨릭대학교 산학협력단 | 패혈증 고위험군 예측 방법 및 시스템 |
SG11201510282PA (en) * | 2013-06-28 | 2016-01-28 | Acumen Res Lab Pte Ltd | Sepsis biomarkers and uses thereof |
US10260111B1 (en) * | 2014-01-20 | 2019-04-16 | Brett Eric Etchebarne | Method of detecting sepsis-related microorganisms and detecting antibiotic-resistant sepsis-related microorganisms in a fluid sample |
GB201402293D0 (en) * | 2014-02-11 | 2014-03-26 | Secr Defence | Biomarker signatures for the prediction of onset of sepsis |
US10761093B2 (en) | 2014-10-09 | 2020-09-01 | Texas Tech University System | Microdevice for cell separation utilizing activation phenotype |
WO2017165693A1 (en) * | 2016-03-23 | 2017-09-28 | Peach Intellihealth, Inc. | Use of clinical parameters for the prediction of sirs |
KR102287968B1 (ko) * | 2019-03-29 | 2021-08-09 | 연세대학교 산학협력단 | 균혈증 발병 위험도 예측 방법 및 이를 이용한 균혈증 발병 위험도 예측용 디바이스 |
DE102022114248A1 (de) | 2022-06-07 | 2023-12-07 | TCC GmbH | Verfahren sowie Vorhersagesystem zur Ermittlung der Eintrittswahrscheinlichkeit einer Sepsis eines Patienten |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6190872B1 (en) * | 1994-05-06 | 2001-02-20 | Gus J. Slotman | Method for identifying and monitoring patients at risk for systemic inflammatory conditions and apparatus for use in this method |
US6804551B2 (en) * | 1998-03-17 | 2004-10-12 | University Of Virginia Patent Foundation | Method and apparatus for the early diagnosis of subacute, potentially catastrophic illness |
WO2000010455A1 (en) * | 1998-08-24 | 2000-03-02 | Emory University | Method and apparatus for predicting the onset of seizures based on features derived from signals indicative of brain activity |
AU3862200A (en) | 1999-03-01 | 2000-09-21 | Board Of Trustees Of The Leland Stanford Junior University | Rapid assay for infection in small children |
DE10155600B4 (de) * | 2001-11-09 | 2009-08-27 | Oligene Gmbh | Nukleinsäure-Array |
US20040038201A1 (en) * | 2002-01-22 | 2004-02-26 | Whitehead Institute For Biomedical Research | Diagnostic and therapeutic applications for biomarkers of infection |
US7774143B2 (en) * | 2002-04-25 | 2010-08-10 | The United States Of America As Represented By The Secretary, Department Of Health And Human Services | Methods for analyzing high dimensional data for classifying, diagnosing, prognosticating, and/or predicting diseases and other biological states |
BR0316231A (pt) * | 2002-11-12 | 2005-10-04 | Becton Dickinson Co | Métodos para determinar o estado de sepsia para prognosticar o começo de sepsia e para diagnosticar a sìndrome de resposta inflamatória sistêmica em um indivìduo e para isolar um biomarcador, perfil biomarcador r kit |
CA2505785A1 (en) * | 2002-11-12 | 2004-05-27 | Becton, Dickinson And Company | Diagnosis of sepsis or sirs using biomarker profiles |
WO2004044555A2 (en) * | 2002-11-12 | 2004-05-27 | Becton, Dickinson And Company | Diagnosis of sepsis or sirs using biomarker profiles |
EP1673465A4 (en) * | 2003-09-29 | 2008-04-30 | Biosite Inc | PROCESS AND COMPOSITIONS FOR SEPSIS DIAGNOSIS |
TW200519384A (en) * | 2003-12-12 | 2005-06-16 | We Gene Technologies Inc | Sepsis detection chip and fabricating method thereof and method of detecting sepsis |
US20050196817A1 (en) * | 2004-01-20 | 2005-09-08 | Molecular Staging Inc. | Biomarkers for sepsis |
GB0610078D0 (en) * | 2006-05-20 | 2006-06-28 | Secr Defence | Sepsis detection microarray |
-
2004
- 2004-12-09 GB GBGB0426982.5A patent/GB0426982D0/en not_active Ceased
-
2005
- 2005-12-09 US US11/792,724 patent/US20080114576A1/en not_active Abandoned
- 2005-12-09 JP JP2007544987A patent/JP2008522604A/ja active Pending
- 2005-12-09 EP EP05825137A patent/EP1839218A1/en not_active Withdrawn
- 2005-12-09 CA CA002587270A patent/CA2587270A1/en not_active Abandoned
- 2005-12-09 AU AU2005313114A patent/AU2005313114B2/en not_active Ceased
- 2005-12-09 WO PCT/GB2005/004755 patent/WO2006061644A1/en active Application Filing
- 2005-12-09 CN CN2005800424703A patent/CN101076806B/zh not_active Expired - Fee Related
- 2005-12-09 EP EP11008209A patent/EP2458517A1/en not_active Withdrawn
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103649329A (zh) * | 2010-11-26 | 2014-03-19 | ImmuneXpress有限公司 | 诊断和/或筛选剂及其用途 |
CN108431832A (zh) * | 2015-12-10 | 2018-08-21 | 渊慧科技有限公司 | 利用外部存储器扩增神经网络 |
CN108604466A (zh) * | 2016-01-28 | 2018-09-28 | 贝克曼考尔特公司 | 感染检测和区分系统及方法 |
CN108604466B (zh) * | 2016-01-28 | 2022-10-11 | 贝克曼考尔特公司 | 感染检测和区分系统及方法 |
CN111989032A (zh) * | 2018-06-15 | 2020-11-24 | 拜克门寇尔特公司 | 存在母细胞标志时确定败血症的方法 |
CN111374639A (zh) * | 2018-12-28 | 2020-07-07 | 陈柏龄 | 败血症的菌种预测系统与方法 |
CN112447289A (zh) * | 2019-08-30 | 2021-03-05 | 希尔-罗姆服务公司 | 败血症监测系统 |
WO2021102776A1 (zh) * | 2019-11-28 | 2021-06-03 | 北京肿瘤医院(北京大学肿瘤医院) | 一种用于早期体腔内感染性并发症诊断的标记物及方法 |
CN114167062A (zh) * | 2020-09-10 | 2022-03-11 | 北京肿瘤医院(北京大学肿瘤医院) | 一种用于早期腹腔内感染性并发症诊断的标记物及方法 |
CN114167062B (zh) * | 2020-09-10 | 2024-02-02 | 北京肿瘤医院(北京大学肿瘤医院) | 一种用于早期腹腔内感染性并发症诊断的标记物及方法 |
TWI803893B (zh) * | 2021-06-28 | 2023-06-01 | 中國醫藥大學 | 敗血症之人工智慧輔助醫療診斷方法及其系統 |
Also Published As
Publication number | Publication date |
---|---|
WO2006061644A1 (en) | 2006-06-15 |
AU2005313114A1 (en) | 2006-06-15 |
EP1839218A1 (en) | 2007-10-03 |
CN101076806B (zh) | 2012-01-11 |
CA2587270A1 (en) | 2006-06-15 |
EP2458517A1 (en) | 2012-05-30 |
AU2005313114B2 (en) | 2012-01-19 |
JP2008522604A (ja) | 2008-07-03 |
GB0426982D0 (en) | 2005-01-12 |
US20080114576A1 (en) | 2008-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101076806A (zh) | 败血症的早期检测 | |
Zhang et al. | Defining inflammatory cell states in rheumatoid arthritis joint synovial tissues by integrating single-cell transcriptomics and mass cytometry | |
Katzenelenbogen et al. | Coupled scRNA-seq and intracellular protein activity reveal an immunosuppressive role of TREM2 in cancer | |
US6235504B1 (en) | Methods for identifying genomic equivalent markers and their use in quantitating cells and polynucleotide sequences therein | |
US7711492B2 (en) | Methods for diagnosing lymphoma types | |
CN101541976A (zh) | 利用血液白细胞中的基因表达标记物的急性感染的鉴别诊断 | |
US20150315652A1 (en) | Method for Determining Coronary Artery Disease Risk | |
US20060073496A1 (en) | Methods of identifying patients at risk of developing encephalitis following immunotherapy for Alzheimer's disease | |
CN1313891A (zh) | 使用基因表达特征监控疾病状态和治疗的方法 | |
CN1647067A (zh) | 数据解析装置及方法 | |
CN1335893A (zh) | 细胞表达特征的统计合并 | |
CN101378764A (zh) | 通过血液白细胞微阵列分析对系统性红斑狼疮的诊断、预后和疾病发展的监测 | |
CN1811447A (zh) | 核酸薄膜层析快速检测方法及其试纸条及其用途 | |
CN1268033A (zh) | 用来预测未来健康的系统 | |
CN1298486A (zh) | 评价化学和生物学分析的方法 | |
Yuan et al. | Single-cell microfluidics to study the effects of genome deletion on bacterial growth behavior | |
Kiss et al. | Interleukin-3 coordinates glial-peripheral immune crosstalk to incite multiple sclerosis | |
CN103097549A (zh) | 肾移植耐受的确定方法 | |
Regueira‐Iglesias et al. | Critical review of 16S rRNA gene sequencing workflow in microbiome studies: From primer selection to advanced data analysis | |
Foulkes et al. | Characterizing the relationship between HIV‐1 genotype and phenotype: prediction‐based classification | |
CN1809644A (zh) | 第17号染色体数量性状基因座精细定位及其在标记辅助选择中的应用 | |
WO2009063249A2 (en) | Early detection of sepsis | |
Chevignon et al. | De novo transcriptome assembly and analysis of the flat oyster pathogenic protozoa Bonamia ostreae | |
CN107960107A (zh) | 测量嵌合状态的方法 | |
Gill et al. | The TNF (− 308A) polymorphism is associated with microchimerism in transfused trauma patients |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120111 Termination date: 20131209 |