基于生物标志物的诊断青光眼的产品及其应用

摘要

本发明公开了基于生物标志物的诊断青光眼的产品及其应用，所述生物标志物选自DFNA5、MYH11和/或KHDRBS1。相比正常对照，DFNA5、KHDRBS1在青光眼患者中呈现显著上调，MYH11在青光眼患者中呈现显著下调，诊断效能结果显示，DFNA5、MYH11和/或KHDRBS1应用于青光眼的诊断具有较高的准确性、敏感性和特异性。

权利要求

1.诊断青光眼的产品在制备诊断青光眼的工具中的应用，其特征在于，所述产品包括检测样本中生物标志物DFNA5、MYH11和/或KHDRBS1水平的试剂。

2.根据权利要求1所述的应用，其特征在于，所述样本中DFNA5、MYH11和/或KHDRBS1水平通过测量样本中DFNA5、MYH11和/或KHDRBS1的蛋白水平或mRNA水平来确定。

3.根据权利要求2所述的应用，其特征在于，所述样本中DFNA5、MYH11和/或KHDRBS1的蛋白水平通过使用免疫染色、免疫荧光、蛋白质印迹或ELISA测量。

4.根据权利要求2所述的应用，其特征在于，所述样本中的DFNA5、MYH11和/或KHDRBS1的mRNA水平通过使用微阵列、RNA‑seq、原位杂交、RNA‑scope以及常规的半定量或定量RT‑PCR测量。

5.根据权利要求1‑4任一项所述的应用，其特征在于，所述产品还包括处理样本的试剂。

6.检测样本中生物标志物的试剂在制备诊断青光眼的产品中的应用，其特征在于，所述生物标志物包括DFNA5、MYH11和/或KHDRBS1。

7.根据权利要求1所述的应用，其特征在于，所述试剂包括测序技术、核酸杂交技术、核酸扩增技术、蛋白免疫技术检测生物标志物水平的试剂。

8.根据权利要求7所述的应用，其特征在于，所述试剂选自：特异性识别所述生物标志物的探针；特异性扩增所述生物标志物的引物；或特异性结合所述生物标志物的抗体。

9.根据权利要求6‑8任一项所述的应用，其特征在于，所述样本选自组织、血液。

说明书

基于生物标志物的诊断青光眼的产品及其应用

技术领域

[0001] 本发明涉及生物医学领域，涉及基于生物标志物的诊断青光眼的产品及其应用。

背景技术

[0002] 青光眼(glaucoma)是一组以视网膜神经节细胞(Retinal ganglion cells,RGCs)；死亡为重要特征，以特征性视野缺损和视神经萎缩为共同表现的疾病；作为目前全
球第二大致盲性眼病，青光眼导致了超过7000万人受累，严重威胁着人类的视觉健康
（Casson RJ, Chidlow G, Wood JP, et al. Definition of glaucoma: clinical and
experimental concepts[J]. Clip Exp Ophthalmol, 2012, 40(4):341‑349.）。

[0003] 青光眼在临床上分为先天性、原发性以及继发性三大类，其中原发性青光眼是指病因机制尚不明确的一类青光眼，根据前房角的状态为开放或关闭，原发性青光眼又可分
为原发性开角型青光眼(Primary open‑angle glaucoma, POAG)和原发性闭角型青光眼
(Primary angle‑closure glaucoma, PACG)(Quigley HA .Glaucoma[J]. Lancet, 2011,
377(9774):1367‑1377.)；POAG的临床特征表现为眼压虽升高，但前房角始终开放，是青光
眼最常见的临床类型，全球大约有70%的青光眼类型为POAG (Beidoe G, Mousa SA.
Current primary open‑angle glaucoma treatments and future directions[J]. Clin
Ophthalmol, 2012, 6(1):1699‑1707.。)

[0004] POAG作为青光眼的一种重要类型，其特征是获得性视神经萎缩、RGC及其轴突丢失且眼压升高时房角始终保持开放（中华医学会眼科学分会青光眼学组.我国原发性青光眼
诊断和治疗专家共识(2014年) [J].中华眼科杂志，2014,05(00).），是慢性进行性发展的
视神经病变。眼压是青光眼患者唯一可控的因素，有效的控制是青光眼的一种重要管理手
段（Mantravadi A V, Vadhar N. Glaucoma[J]. Prim Care, 2015,42(3): 437‑449.）。但
正常眼压型POAG的存在和POAG隐匿发病的特点会导致部分患者漏诊，确诊时通常已存在视
神经不可逆的损伤。为最大可能的挽救POAG患者的视功能，早期诊断具有重要的意义。

发明内容

[0005] 为了弥补现有技术的不足，本发明基于遗传因素在青光眼的发生发展中的作用，研究与青光眼发生发展相关的生物标志物，从而为青光眼的诊断和治疗提供新的手段。

[0006] 本发明第一方面提供了一种诊断青光眼的产品，所述产品包括检测样本中生物标志物DFNA5、MYH11和/或KHDRBS1水平的试剂。

[0007] 进一步，所述样本中DFNA5、MYH11和/或KHDRBS1水平通过测量样本中DFNA5、MYH11和/或KHDRBS1的蛋白水平或mRNA水平来确定。

[0008] 进一步，所述样本中DFNA5、MYH11和/或KHDRBS1的蛋白水平通过使用免疫染色、免疫荧光、蛋白质印迹或ELISA测量。

[0009] 进一步，所述样本中的DFNA5、MYH11和/或KHDRBS1的mRNA水平通过使用微阵列、RNA‑seq、原位杂交、RNA‑scope以及常规的半定量或定量RT‑PCR测量。

[0010] 进一步，所述产品还包括处理样本的试剂。

[0011] 本发明第二方面提供了本发明第一方面所述的产品在制备诊断青光眼的工具中的应用。

[0012] 本发明第三方面提供了检测样本中生物标志物的试剂在制备诊断青光眼的产品中的应用，所述生物标志物包括DFNA5、MYH11和/或KHDRBS1。

[0013] 进一步，所述试剂包括测序技术、核酸杂交技术、核酸扩增技术、蛋白免疫技术检测生物标志物水平的试剂。

[0014] 进一步，所述试剂选自：特异性识别所述生物标志物的探针；特异性扩增所述生物标志物的引物；或特异性结合所述生物标志物的抗体。

[0015] 进一步，所述样本包括但不限于组织或流体，如组织、血液、血浆、血清、淋巴液、尿液、浆膜腔液、脊髓液、滑液、房水、泪液、唾液、或其组分或经过处理的材料。

[0016] 优选地，所述样本选自组织、血液。

[0017] 本发明的有益效果：

[0018] 本发明通过检测DFNA5、MYH11和/或KHDRBS1的表达水平，可以实现青光眼的诊断，增加检测的敏感度，提高检测能力和效率，提升指导青光眼的临床治疗的能力。

附图说明

[0019] 图1显示DFNA5基因差异表达图；

[0020] 图2显示MYH11基因差异表达图；

[0021] 图3显示KHDRBS1基因差异表达图；

[0022] 图4显示DFNA5基因诊断青光眼的ROC曲线图；

[0023] 图5显示MYH11基因诊断青光眼的ROC曲线图；

[0024] 图6显示KHDRBS1基因诊断青光眼的ROC曲线图；

[0025] 图7显示DFNA5+MYH11基因诊断青光眼的ROC曲线图；

[0026] 图8显示DFNA5+KHDRBS1基因诊断青光眼的ROC曲线图；

[0027] 图9显示MYH11+KHDRBS1基因诊断青光眼的ROC曲线图；

[0028] 图10显示DFNA5+MYH11+KHDRBS1联合诊断青光眼的ROC曲线图。

具体实施方式

[0029] 以下将对本发明进一步详细说明，应理解，所述用语旨在描述目的，而非限制本发明。

[0030] 术语 “丰度” 、 “水平” 和 “含量” 在本文可互换使用，是指定量的含量（例如重量或摩尔）、半定量的含量、相对含量（例如，等级内的重量%或摩尔%）、浓度等。因此，这些术
语涵盖样品中疾病治疗生物标志物的绝对或相对含量或浓度。

[0031] 术语“和/或”是指并且包括一个或多个相关联的所列项目的任何和所有可能的组合，以及在备选方案（或）中解释时缺少组合。

[0032] 在本发明中，术语“生物标志物”意指化合物，优选是基因，与来自具有第二表型（例如没有疾病）的受试者或一组受试者的生物样品相比，它在来自具有第一表型（例如患
有疾病）的受试者或一组受试者的生物样品中差异地存在（即增加或减少）。术语“生物标志
物”通常是指一种基因的存在/浓度/含量或两种或更多种基因的存在/浓度/含量。

[0033] 术语 “生物标志物值” 或“生物标志物水平”是指针对受试者的至少一种相应的生物标志物测量或推导的值，并且其通常至少部分指示取自受试者的样品中的生物标志物
的丰度或浓度。因此，生物标志物值可以是测量的生物标志物值，其是针对受试者测量的生
物标志物值，或可替换地可以是推导的生物标志物值，其是推导自一个或多个测量的生物
标志物值的值，例如，通过将函数应用于一个或多个测量的生物标志物值推导的。生物标志
物值可以是任何合适的形式，这取决于其中确定值的方式。例如，可以使用高通量技术，诸
如测序平台、阵列和杂交平台、质谱、免疫测定、免疫荧光、流式细胞术或这些技术的任意组
合来确定生物标志物值。在一个优选的实例中，生物标志物值涉及表达产物或其他可测量
分子的丰度或活性水平，使用如定量RT‑PCR、测序等技术来定量。在这种情况中，如本领域
技术人员已知的，生物标志物值可以是扩增含量或循环次数的形式，其是样品内生物标志
物浓度的对数表示。在其他优选实例中，使用含有表达产物的细胞的免疫荧光来定量生物
标志物值。

[0034] 在本发明中，所述生物标志物包括DFNA5、MYH11和/或KHDRBS1。

[0035] 在本发明中，DFNA5（gasdermin E，gene ID：1687）包括DFNA5基因及其同源物，突变，和同等型。该术语涵盖全长，未加工的DFNA5，以及源自细胞中加工的任何形式的DFNA5。
该术语涵盖DFNA5的天然发生变体（例如剪接变体或等位变体）。

[0036] MYH11（myosin heavy chain 11，gene ID：4629）包括MYH11基因及其同源物，突变，和同等型。该术语涵盖全长，未加工的MYH11，以及源自细胞中加工的任何形式的MYH11。
该术语涵盖MYH11的天然发生变体（例如剪接变体或等位变体）。

[0037] KHDRBS1（KH RNA binding domain containing, signal transduction associated 1 ，gene ID：10657）包括KHDRBS1基因及其同源物，突变，和同等型。该术语涵
盖全长，未加工的KHDRBS1，以及源自细胞中加工的任何形式的KHDRBS1。该术语涵盖
KHDRBS1的天然发生变体（例如剪接变体或等位变体）。

[0038] gene ID可在https://www.ncbi.nlm.nih.gov/gene/获得。

[0039] 术语“引物”是指与DNA链配对时能够在合适的聚合剂存在下引发引物延伸产物合成的寡核苷酸。引物优选是单链的，以实现最大的扩增效率，但也可以是双链的。引物必须
足够长，以在聚合剂存在下引发延伸产物的合成。引物的长度取决于许多因素，包括应用、
待使用的温度、模板反应条件、其他试剂和引物来源。例如，根据靶序列的复杂性，引物可以
为至少约5、6、7、8、9、10、11、12、13、14、15、16、17、18、19、20，21、22、23、24、25、26、27、28、
29、30、35、40、50、75、100、150、200、300、400、500到引物3′ 端比模板序列短一个碱基的长
度，以允许核酸链的延伸，尽管引物的5′ 端可以在长度上延伸超过模板序列的3′ 端。在某
些实施方案中，引物可以是大的多核苷酸，如约35个核苷酸至几千个碱基或更多。可以选择
引物，以与模板上被设计成与其杂交的并且用作合成起始位点的序列 “基本上互补”。“基
本上互补”表示引物足够互补以与靶多核苷酸杂交。理想地，引物不含与设计成与其杂交的
模板的错配，但这不是必需的。例如，非互补核苷酸残基可以连接至引物的5′ 端，而引物序
列的其余部分与模板互补。或者，可以将非互补核苷酸残基或一段非互补核苷酸残基链散
布在引物中，只要该引物序列与模板序列具有足够的互补性以与其杂交，从而形成用于合
成引物延伸产品的模板。

[0040] 术语 “探针” 是指结合另一个分子的特定序列或子序列或其他部分的分子。除非另外指出，否则术语 “探针” 通常是指通过互补碱基配对结合另一个核酸（在本文也称为 
“靶多核苷酸” ）的核酸探针。探针可以结合缺少与探针的完整序列互补性的靶多核苷酸，
这取决于杂交条件的严格性。探针可以直接或间接标记并且包括它们范围内的引物。

[0041] 术语“抗体”在本文中以最广泛的含义使用，且包括单克隆抗体（ 例如，全长或完整的单克隆抗体）、多克隆抗体、多价抗体、多特异性抗体（ 例如，双特异性抗体，只要它们
显示希望的生物学活性），且也可以包括某些抗体片段（ 如本文更详细地描述）。抗体可以
是人、人源化和/或亲和力成熟的抗体。

[0042] 生物标志物核酸的检测

[0043] 在一些实施方案中，通过测定生物标志物核酸转录物水平来评估生物标志物。在说明性的基于核酸的测定中，根据标准方法（Sambrook等，1989，上文；和Ausubel等，1994，
上文）从生物样品中所含的细胞分离核酸。核酸通常是分级后的或全细胞RNA。在一些实施
方案中，通过模板依赖性核酸扩增技术来扩增核酸。多种模板依赖性方法可用于扩增给定
模板样品中存在的疾病治疗生物标志物序列。示例性核酸扩增技术是聚合酶链式反应（称
为PCR）。

[0044] 在某些有利的实施方案中，模板依赖性扩增涉及实时的转录物定量。例如，可以使用实时PCR技术来定量RNA或DNA。通过确定完成相同数目的循环并在其线性范围内的PCR反
应中靶DNA扩增产物的浓度，可以确定原始DNA混合物中特定靶序列的相对浓度。如果DNA混
合物是从不同组织或细胞分离的RNA合成的cDNA，则可以针对各个组织或细胞确定衍生该
靶序列的特定mRNA的相对丰度。实时PCR通常使用本领域可用的任何PCR仪器进行。通常，用
于实时PCR数据收集和分析的仪器包括热循环仪，用于荧光激发和发射收集的光学器件，和
任选的计算机和数据采集与分析软件。

[0045] 在某些实施方案中，使用印迹技术来定量靶核酸，这是本领域技术人员公知的。Southern印迹涉及使用DNA作为靶标，而Northern印迹涉及使用RNA作为靶标。各自提供不
同类型的信息，尽管在许多方面，cDNA印迹与印迹或RNA物质类似。简而言之，使用探针来靶
向已经固定于合适基质（常常是硝基纤维素滤器）上的DNA或RNA物质。不同的物质应在空间
上分开，以利于分析。这通常是通过核酸物质的凝胶电泳，然后 “印迹” 到滤器上来完成
的。随后，在促进变性和再杂交的条件下，将印迹的靶标与探针（通常是标记的）一起孵育。
去除未结合的探针，完成检测。在检测/定量之后，可以将在给定受试者中观察到的结果与
本文统计学显著的参照组或对照受试者群进行比较。这样，可以将检测到的疾病生物标志
物核酸的量与疾病的进展或严重程度相关联。

[0046] 芯片杂交利用连接固体基底的生物标志物特异性寡核苷酸，其可以由设计成微阵列的颗粒状固相组成，如尼龙滤器、载玻片或硅片。微阵列在本领域是已知的并且由表面组
成，序列对应于基因产物（如cDNA）的探针可以在所述表面上的已知位置特异性杂交或结
合，以检测生物标志物基因表达。

[0047] 杂交复合物的定量是本领域公知的并且可以通过几种方法中的任何一种来实现。这些方法通常是基于标记或标志物的检测，如本领域标准使用的任何放射性、荧光、生物或
酶标签或标记。可以将标记应用于寡核苷酸探针或源自生物样品的RNA。

[0048] 通常，mRNA定量可与校准曲线一起合适地进行，以实现准确的mRNA测定。此外，优选通过与对照样品比较来对源自生物样品的转录物进行定量，所述样品的特征在于所检查
的转录物的已知表达模式。

[0049] 生物标志物蛋白的检测

[0050] 在一些实施方案中，通过证明蛋白质的存在（分离的或一种或在细胞中） ，或通过生物标记物的一种或多种已知功能特性，在蛋白质表达水平上评估疾病治疗生物标记物。
例如，用于DFNA5特异性或MYH11特异性或KHDRBS1特异性蛋白检测中的抗DFNA5、MYH11、
KHDRBS1抗体是本领域已知的，可商购，并且也可以由本领域技术人员容易地生产。抗体和
抗原‑抗体复合物可以通过本领域几种熟知的测定法进行检测，包括免疫荧光测定法，免疫
组织化学，荧光激活细胞分选（FACS）分析，酶联免疫吸附测定法（ELISA） ，放射免疫测定法
（RIA） ，光发射免疫分析和蛋白质印迹分析。

[0051] 在特定的实施方案中，进行免疫荧光或免疫细胞化学来检测蛋白质。细胞，如患病组织细胞，可以通过本领域已知的方法来分离或富集。细胞的分离或富集是指其中特定细
胞（例如，患病组织细胞）的百分比相对于富集程序前样品中的百分比增加的过程。纯化是
富集的一个实例。在其他实施方案中，可以将针对细胞上的表面标志物的抗体连接到固体
TM
支持物上，以进行分离。用于分离的程序可包括使用抗体磁珠（例如Miltenyi 珠）的磁分
离、亲和色谱法、使用连接固体基质的抗体 “淘选” 或任何其他方便的技术，例如Laser 
Capture Microdissection。提供特别精确的分离的其他技术包括FACS。一旦将细胞沉积在
载玻片上，就可以将其固定，并用标记的抗体探测以检测疾病诊断生物标志物。

[0052] 疾病生物标志物特异性的抗体可以直接缀合荧光标志物，包括荧光素、FITC、罗丹明、Texas Red、Cy3、Cy5、Cy7和其他荧光标志物，并在配备有合适滤镜的荧光显微镜下观察
滤器。抗体也可以与酶缀合，酶在添加适当的底物后开始反应，从而在细胞上提供带有检测
到的生物标记蛋白的有色沉淀。然后可以通过标准光学显微镜观察载玻片。或者，可以将疾
病诊断生物标志物特异性的一抗进一步结合缀合至可检测部分的二抗。

[0053] 免疫组织化学原理上与免疫荧光或免疫细胞化学非常相似，然而，例如，与细胞悬浮液相反，用疾病治疗生物标志物特异性的抗体探测组织标本。固定并处理活检标本，并任
选将其包埋在石蜡中，如果需要，切片，以提供随后用乙酰肝素酶特异性抗体探测的细胞或
组织载玻片。或者，可将冷冻组织低温恒温器上切片，然后进行抗体探测，从而避免了固定
诱导的抗原掩盖。将抗体，如在免疫荧光或免疫细胞化学中的，偶联至荧光或酶联的可检测
部分，并通过针对免疫荧光描述的方法用于探测组织切片，然后根据所用的检测方法通过
荧光或共聚焦显微镜观察。反应产物形成后，如果利用了酶促可检测的部分，则可以通过标
准光学显微镜观察反应产物沉淀物的可视化。

[0054] 在其他实施方案中，使用如ELISA和RIA这样的测定，其遵循类似的原理用于检测特异性抗原。作为说明性实例，可以通过通常用125I放射性标记的DFNA5、MYH11或KHDRBS1
特异性抗体使用RIA来测量DFNA5、MYH11或KHDRBS1。在ELISA测定中，DFNA5、MYH11或
KHDRBS1特异性抗体与酶化学连接。DFNA5、MYH11或KHDRBS1特异性捕获抗体被固定在固体
支持物上。然后在其中阻断非特异性结合的条件下，将未标记的样本，例如来自生物样品的
蛋白质提取物，与固定的抗体一起温育，并通过洗涤除去未结合的抗体和/或蛋白质。结合
的DFNA5、MYH11或KHDRBS1通过第二个DFNA5、MYH11或KHDRBS1特异性标记抗体来检测。在
RIA中，通过测量放射性直接测量抗体结合，而在ELISA中，结合是通过将无色底物转化为有
色反应产物的反应来检测的，作为连接的酶的活性的函数。因此，可以通过分光光度法容易
地检测到变化。

[0055] 蛋白质生物标志物表达也可以通过发光免疫测定来检测。与ELISA和RIA极为相似，在发光免疫测定中，将待测试的生物样品/蛋白质提取物固定在固体支持物上，并用特
异性标记（标记的抗体）进行探测。该标记又是发光的，并在结合时发光，以作为特异性识别
的指示。发光标记包括在通过电磁辐射、电化学激发或化学激活而在激活时发光的物质，并
且可以包括荧光和磷光物质、闪烁体和化学发光物质。标记可以是催化反应系统的一部分，
如酶、酶片段、酶底物、酶抑制剂、辅酶或催化剂；色原系统的一部分，如荧光团、染料、化学
发光剂、发光剂或敏化剂；可分散颗粒（可以是非磁性或磁性的）、固体支持物、脂质体、配
体、受体、半抗原放射性同位素等。

[0056] Western印迹分析是评估生物样品中疾病生物标志物多肽含量的另一种方法。来自细胞（例如患病组织细胞）的生物样品的蛋白质提取物在变性电离环境中溶解，并将等分
试样施加到聚丙烯酰胺凝胶基质上。随着向阳极迁移，蛋白质会基于分子大小特性进行分
离。然后将抗原转移到硝酸纤维素膜、PVDF或尼龙膜上，然后进行膜阻断以最小化非特异性
结合。用与可检测部分直接偶联的抗体探查膜，或随后用含有可检测部分的二抗探查膜。通
常，将辣根过氧化物酶或碱性磷酸酶与抗体偶联，并使用发色或发光底物使活性可视化。

[0057] 在特定的实施方案中，使用允许同时检测和/或定量大量蛋白质的蛋白质捕获阵列。例如，滤膜上的低密度蛋白质阵列，现在有可能使用蛋白质阵列来分析体液中蛋白质的
表达，如健康或患病受试者的血清中以及药物治疗前后的受试者中的蛋白质表达。示例性
蛋白质捕获阵列包括包含空间寻址的抗原结合分子的阵列，通常称为抗体阵列，其可以促
进限定蛋白质组或亚蛋白质组的多种蛋白质进行广泛的平行分析。已经显示抗体阵列具有
所需的特异性和可接受背景的特性。

[0058] 诊断产品

[0059] 本发明提供了一种诊断青光眼的产品，所述产品包括检测样本中本发明所述的生物标志物的试剂；并且可包括使用所述产品评估受试者是否患有或易患青光眼的说明书。

[0060] 诊断产品还可以任选包括合适的用于检测标记的试剂、阳性和阴性对照、洗涤溶液、印迹膜、微滴定平板、稀释缓冲液等。例如，基于蛋白质的检测诊断产品可以包括（i）至
少一种疾病生物标志物多肽；和（ii）特异性地结合疾病治疗生物标志物多肽的抗体。或者，
基于核酸的测试剂盒可以包括（i）疾病治疗生物标志物多核苷酸；和（ii）与疾病治疗生物
标志物多核苷酸特异性地杂交的引物或探针。还可以包括适用于扩增核酸的酶，包括各种
TM
聚合酶（逆转录酶、Taq、Sequenase 、DNA连接酶等），脱氧核糖核苷酸和缓冲液，以提供扩增
所需的反应混合物。此类试剂盒通常还将以合适的方式包含用于每种单独的试剂和酶以及
每种引物或探针的不同容器。

[0061] 可使用本文所述的能够检测样本生物标志物的任何形式的样本测定。通常，所述测定将定量样本中生物标志物至一定的程度，例如它们的浓度或量是高于还是低于预定阈
值。此类试剂盒可采取测试条、浸杆、盒、药筒、基于芯片或基于珠粒的阵列、多孔板或一系
列容器等的形式。提供一种或多种试剂以检测所选样本生物标志物的存在和/或浓度和/或
量。可将受试者的样本直接分配到测定中，或从存储的或先前获得的样品中间接分配到测
定中。

[0062] 下面结合附图和实施例对本发明作进一步详细的说明。以下实施例仅用于说明本发明而不用于限制本发明的范围。实施例中未注明具体条件的实验方法，通常按照常规条
件，或按照制造厂商所建议的条件。

[0063] 实施例 与青光眼诊断相关的基因标志物

[0064] 1、数据和预处理

[0065] 从GEO数据库下载青光眼及其对照的数据集GSE27276的基因表达数据，并利用注释文件对其注释，多个探针对应同一个基因的取平均值作为其表达量，然后获得基因表达
矩阵文件。

[0066] 2、差异表达分析

[0067] 使用R软件中的“limma”包进行差异表达分析，差异基因的筛选标准为Pvalue<0.05。

[0068] 分析结果显示，DFNA5、KHDRBS1呈现显著上调，MYH11显著下调，其表达情况如图1‑3所示，其中，*：P<0.05；**：P<0.01；***：P<0.001。

[0069] 3、诊断效能分析

[0070] 使用R包“pROC”ROC曲线，分析差异表达基因作为检测变量的AUC值、敏感性和特异性，判断其诊断效能。在判断每个基因的诊断效能时，直接使用基因的表达量进行分析。首
先调用pROC包，然后读入目的基因构建的表达量矩阵，运行绘制ROC曲线的命令，其中命令
采用for循环，同时涉及添加AUC，thres（阈值）、smooth（拟合曲线）的命令。在判断基因联合
的诊断效能时，首先是使用glmnet对基因进行logistics回归，利用建立的Logistic回归模
型，使用predict（预测）函数观察某个预测变量在各个水平对结果概率的影响，计算预测概
率，绘制预测结果的ROC曲线。

[0071] 结果如表1和图4‑10所示，从表中可以看出DFNA5、MYH11、KHDRBS1及其组合应用于青光眼的诊断具有较高的准确性，尤其是三者的组合，具有较高的准确性、敏感性以及特异
性。

[0072] 表1 差异表达基因诊断效能分析

[0073]

[0074] 以上结合附图详细描述了本申请的优选实施方式，但是，本申请并不限于上述实施方式中的具体细节，在本申请的技术构思范围内，可以对本申请的技术方案进行多种简
单变型，这些简单变型均属于本申请的保护范围。

[0075] 另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合，为了避免不必要的重复，本申请对各种可
能的组合方式不再另行说明。

[0076] 此外，本申请的各种不同的实施方式之间也可以进行任意组合，只要其不违背本申请的思想，其同样应当视为本申请所公开的内容。