一组右眼等效球镜数量性状相关的SNP标志物及其应用的制作方法

一组右眼等效球镜数量性状相关的snp标志物及其应用
技术领域
1.本发明属于生物医学、基因检测技术领域，具体一组右眼等效球镜数量性状相关的snp标志物及其应用。


背景技术：

2.眼科疾病是指发生在眼部区域的疾病，常见眼科疾病有近视、散光、白内障、青光眼、中心浆液性视网膜病变、干眼症、交感性眼炎、夜盲症、弱视、沙眼、糖尿病视网膜病变、结膜炎、老花眼、色盲、视网膜色素变性、视网膜中央动脉阻塞、视网膜脱落、远视、针眼、雪盲症、霰粒肿、飞蚊症等。在调节放松的状态下，平行光线经眼球屈光系统后聚焦在视网膜之前，称为近视。
3.屈光度是指屈光不正状态下所产生的度数大小，当屈光度发生时，一般多代表近视、散光、远视或是弱视等临床表现，近视、远视通常为球镜(s)，散光为柱镜(c)，目前采用等效球镜度数(spherical equivalent refraction，se)评价屈光度，等效球镜度数＝球镜度数 散光度数*1/2。屈光度≤-0.50d可作为近视诊断的循证共识阈值。更具体地，近视的定义是当调节放松时，眼睛的等效球镜屈光不正≤-0.5d的情况。高度近视的定义是当调节放松时，眼睛的等效球镜屈光不正≤-6.00d的情况。低度近视的定义是当调节放松时，眼睛的等效球镜屈光不正≤-0.5且》-6.00d的情况。
4.近视作为一种全球高发性疾病近年来受到国内外越来越多的关注，预计到2050年，近视将影响到全球近50％的人口，同时高度近视也将影响到全球近10％的人口。
5.同时，近视也是目前影响我国青少年人群视觉健康的主要眼病，其中高度近视是引起我国眼病患者眼盲和低视力的主要原因之一，将给患者及其家庭带来严重的经济负担。对于孩子视力问题，预防很重要，而早发现、早干预、早治疗，是防止近视度数加深的关键。少年儿童时期的眼病对视力发育危害极大，及早发现和诊断影响视力的非屈光不正性眼病非常重要，许多的眼病如果不能得到及时的发现和治疗，将造成眼睛的终生残疾。虽然通过药物、光学治疗和行为矫正可以部分减缓其加深速度，但我们距离能够逆转过去几十年的趋势还有很长的路要走。这使得近视及其相关并发症成为了研究的重中之重。高度近视的发生涉及多因素复杂过程，发病机制仍不清楚，当下的各种光学矫正手段以及手术治疗方法，如屈光手术和巩膜加固手术，不能从根本上阻止及延缓高度近视眼底病变的发展，可以说目前缺乏行之有效的高度近视的治疗措施。由此可见，寻找更有效的方法对高度近视患者和高危人群进行早期检测、风险预测和早期干预具有重要的临床意义。
6.单核苷酸多态性(single nucleotide polymorphisms，snp)是指在基因组水平上单个碱基的转换、颠换或由碱基的插入、缺失引起核苷酸变异导致的dna序列多态性，是人类可遗传的变异中最常见的一种，占所有已知多态性的80％以上，人群中变异频率高于1％，这也是区别于点突变的一个重要因素。在遗传学分析中，snps具有高频、稳定和易分析等特点。研究发现，snps与近视的发生发展具有较强的关联性。snps检测主要方法有时间飞行质谱(maldi-tofms)技术、荧光定量pcr技术、基因芯片技术、变形高效液相色谱等。
7.随着高通量snp检测技术方法的出现，作为数量最多且易于批量检测的多态标记，snp在连锁分析与基因定位，包括复杂疾病的基因定位、关联分析、个体和群体对环境致病因子与药物的研究中将发挥愈来愈重要的作用。随着技术的发展，snp检测费用也越来越经济，snp标志物将成为新一代分子标记。


技术实现要素：

8.本发明纳入大量6-18岁的受试人群，通过标准操作程序(sop)采集符合标准的口腔黏膜样本，提取dna，对单核苷酸多态性进行分型，筛选寻找到与右眼等效球镜高度相关的高特异性和敏感性的snp，进一步可研制出相应便于临床应用的试剂盒，为判断受试者右眼等效球镜提供数据支持，可在眼科问题发生早期进行干预，降低眼科疾病恶化的可能性，提高生活质量。
9.为实现本发明的技术目的，本发明提供以下技术方案：
10.第一方面，本发明提供了检测rs2547319的试剂在判断右眼等效球镜中的应用。
11.另一方面，本发明提供了一组右眼等效球镜相关的snp位点组合，所述snp位点组合包含以下任意一组：
12.(1)rs10889315、rs75144949、rs149213974、rs754615、rs45583335、rs2529438、rs2072236、rs2240403、rs78155900、rs114097201、rs2811736、rs2297722、rs2297723、rs117519669、rs7157977、rs10143899、rs3759779、rs143477571、rs2303221、rs34693726、rs3803752、rs142921938、rs3746295、rs2547319、rs2547318、rs4911402(26)；
13.(2)rs4648600、rs1061624、rs2275365、rs35909785、rs3103778、rs2457829、rs10889315、rs78587889、rs12568050、rs12564283、rs528123098、rs139842833、rs74953908、rs1077827、rs12492484、rs75144949、rs3732755、rs3135115、rs16837029、rs149213974、rs17001890、rs56382032、rs12651031、rs2736100、rs7724759、rs754615、rs5369、rs150157174、rs130068、rs45583335、rs78342561、rs2529438、rs1049402、rs2072236、rs2527878、rs2240403、rs78155900、rs114097201、rs2680903、rs59612871、rs3818584、rs16909677、rs10983755、rs2248077、rs2811736、rs2297722、rs2297723、rs9987876、rs7904554、rs907609、rs117519669、rs3750996、rs3802799、rs2276296、rs139692587、rs904628、rs11551723、rs144796593、rs78648016、rs7157977、rs10143899、rs3759779、rs7142098、rs143477571、rs3736911、rs11849022、rs77481241、rs3212112、rs9282879、rs55958706、rs3212056、rs34100926、rs114770841、rs3742368、rs118146919、rs2289702、rs2303221、rs34693726、rs3744566、rs9907420、rs375644567、rs185576254、rs3803752、rs17822735、rs17222523、rs34818467、rs12946397、rs3809724、rs181969864、rs7236574、rs142921938、rs3746295、rs151309111、rs533816037、rs2302948、rs182509214、rs2547319、rs2547318、rs2281209、rs4911402、rs138894461、rs5756223(102个)；
14.(3)rs16824585、rs4648600、rs17033266、rs1061624、rs202140107、rs2291804、rs2275365、rs139241751、rs35909785、rs6676052、rs3103778、rs116893711、rs116861599、rs2457829、rs10889315、rs17855078、rs78587889、rs140141416、rs4656197、rs45445497、rs143355122、rs12568050、rs12564283、rs10482、rs1929861、rs75401237、rs28517482、
rs1264208208、rs2631840、rs188480146、rs528123098、rs61732269、rs367931009、rs139842833、rs79908358、rs79358798、rs1434087、rs4894048、rs13421990、rs13419085、rs143196218、rs139713392、rs6734083、rs150887565、rs4663795、rs74953908、rs143001518、rs6437368、rs147838565、rs1077827、rs144623840、rs78233496、rs9876921、rs6622、rs2279909、rs547922447、rs300977、rs12492484、rs75144949、rs3732755、rs3135115、rs16837029、rs74787220、rs527903313、rs149213974、rs17001890、rs79769348、rs10025654、rs76996680、rs2306369、rs79482316、rs56382032、rs12651031、rs2736100、rs35130836、rs3749684、rs7724759、rs754615、rs45625534、rs3805625、rs137898880、rs147518598、rs202007699、rs2270900、rs75164992、rs17133512、rs76720430、rs5369、rs150157174、rs7756481、rs117298661、rs117644703、rs3094672、rs4713420、rs130068、rs67196728、rs185990098、rs9349593、rs143580489、rs3757302、rs3218602、rs9495370、rs45583335、rs78342561、rs117358940、rs9689983、rs2529438、rs1049402、rs2072236、rs2527878、rs2240403、rs781980620、rs56083327、rs2245368、rs2430307、rs3748126、rs36028884、rs78155900、rs114097201、rs13438494、rs117211799、rs115734214、rs148474510、rs200506987、rs8190955、rs8178175、rs2680903、rs117676215、rs115507207、rs369806130、rs78240711、rs59612871、rs117169590、rs3739523、rs145293869、rs139212809、rs3818584、rs16909677、rs141549766、rs147026086、rs188197626、rs13321、rs10983755、rs2301576、rs7038042、rs11849、rs2275161、rs2248077、rs2811736、rs2297722、rs2297723、rs9987876、rs7904554、rs148582275、rs3740168、rs149647444、rs7100474、rs7079648、rs202166096、rs60420182、rs36104328、rs7098255、rs3814163、rs2297785、rs3736924、rs3736923、rs141296329、rs117520659、rs57285449、rs112687793、rs73401681、rs79948463、rs10902170、rs10902171、rs867839、rs907608、rs907609、rs117519669、rs3750996、rs214086、rs3802799、rs76215382、rs147852038、rs2276296、rs143930161、rs183113103、rs139692587、rs10792769、rs7102675、rs116878689、rs35466364、rs142908193、rs7966459、rs7135745、rs34181394、rs1921038、rs201814780、rs11068531、rs904628、rs11551723、rs7983667、rs144796593、rs140526815_rs77878407、rs78646508、rs34523940、rs9521906、rs78648016、rs7157977、rs10143899、rs3759779、rs7142098、rs41309252、rs535020964、rs3742770、rs3742768、rs116881452、rs143477571、rs61992274、rs142973720、rs530110551、rs3742464、rs3736911、rs11849022、rs77481241、rs3212112、rs9282879、rs55958706、rs3212056、rs34100926、rs114770841、rs3742368、rs144857261_rs373032548、rs189734192、rs188393827、rs118146919、rs562829514、rs1805025、rs8192352、rs2289702、rs150921132、rs61869621、rs1802752、rs113355794、rs7198064、rs2291407、rs9941128、rs9941210、rs9941215、rs9931527、rs12708649、rs13306653、rs11643517、rs5723、rs7196736、rs2303221、rs2304483、rs149645127、rs11550470、rs17232091、rs150889000、rs16954357、rs371756524_rs557159990、rs192103504、rs115314145、rs34693726、rs2277651、rs2277652、rs71360269、rs71360270、rs3744566、rs3744568、rs9907420、rs375644567、rs185576254、rs2158460、rs3803752、rs17822735、rs17222523、rs34818467、rs12946397、rs3809724、rs35476409_rs397978887、
rs3744214、rs181969864、rs4789773、rs28561791、rs7236574、rs76597875、rs61749945、rs118033234、rs1052696、rs1052692、rs148018996、rs34488963、rs3745379、rs3745385、rs142921938、rs3746295、rs151309111、rs139271842、rs140566569、rs533816037、rs78237760、rs79621739、rs75990051、rs62107869、rs305974、rs1055099、rs41290128、rs76151738、rs16980091、rs2302948、rs182509214、rs2547319、rs2547318、rs16982743、rs2284385、rs2281209、rs4911402、rs80158178、rs7509972、rs150511237、rs2231391、rs2231390、rs762228、rs72564613、rs191574093、rs138894461、rs5999924、rs5750146、rs5756223、rs3747168、rs148345118、rs910796、rs910797、rs910798、rs738712、rs16445、rs3922872、rs56248708、rs9628315、rs8137790(342个)。
15.在本发明中，snp(单核苷酸多态性)是指dna中的单个碱基位置，受试者可以是纯合的或杂合的。本发明的snp位点以“rs
‑”
方式命名，本领域技术人员能够根据上文的rs-命名，从适合的数据库和相关的信息系统如单核苷酸多态性数据库(dbsnp)中确定其确切的位置、核苷酸序列。对于snp，每个位点上具有两种以上等位基因，例如，具有两种等位基因。
16.优选地，所述snp位点组合是根据来着受试者样本的进行检测，而后根据检测结果判断受试者当下右眼等效球镜，同时预测其发生眼科问题的风险。
17.本文所用“样本”是指包含来自或衍生自人类患者的核酸(特别是dna)的任何样品，例如来自患者的体液(血液、唾液、尿液等)、活体组织切片、组织和/或废物。因此，组织活体组织切片、粪便、痰液、唾液、血液、淋巴液等可以容易地筛选snp，基本上任何含有适当核酸的目标组织都可以如此。在一个实施方案中，样本是口腔上皮细胞。这些样本通常由患者在知情同意之后通过标准医学实验室方法获取。样品可以为从患者直接取得的形式，或者可以至少部分地加工(纯化)以除去至少一些非核酸材料；
18.优选地，所述样本是口腔拭子(口腔黏膜样本)；
19.另一方面，本发明还提供了一种检测右眼等效球镜的试剂盒，所述试剂盒中包括检测前述snp位点组合的试剂。
20.优选地，所述检测snp位点组合的试剂包括但不限于以下方法检测snp时所使用的试剂：taqman探针法、测序法、芯片法、飞行质谱仪(maldi-tofms)检测、限制性片段长度多态性法(pcr-rflp)、单链构象多态性法(pcr-sscp)、等位基因特异性pcr(as-pcr)、snapshot法、snplex法、变性高效液相色谱法(dhplc)、变性梯度凝胶电泳法(dgge)。本领域技术人员可选择任一种或几种方法来检测snp位点，只要可以实现snp位点的检测。
21.在本发明中“检测snp的制剂”是指特异地结合到包括在所述snp标志物组合物的snp标志物或snp，以便于识别，或者可以检测所述snp并扩增的制剂，“可以扩增snp标志物的制剂”是指反复复制包括在所述snp标志物组合物的snp标志物或snp，并增加其数量的制剂，例如，特异地扩增包括所述snp的多核苷酸的引物或者可以特异结合的探针，但不受此限定。
22.用于所述snp扩增的引物可以是适当缓冲剂中的妥当条件(例如，4种不同核苷三磷酸、dna、rna聚合酶或逆转录酶等聚合剂)以及在妥当温度下，可以作为柱状显示dna合成的起点，起作用的单链寡核苷酸，根据使用目的，所述引物的妥当长度会有不同，但通常是15至30个核苷酸。短的引物分子通常为了形成柱状和稳定的混合体，需要更低的温度。所述引物序列不需要与包括所述snp的多核苷酸形成完整的互补性，但其互补性需要达到可以
与包括所述snp的多核苷酸混合的程度。
23.在一种实施方式中，所述检测snp位点组合的试剂是扩增前述snp位点组合的引物组和/或探针组。
24.在一种实施方式中，所述试剂盒还可以包括相应pcr技术所需要的常用试剂，如dntps，mgcl2，双蒸水，荧光探针等，这些常用试剂都是本领域技术人员熟知，另外还有标准品和对照(如基因型标准品和空白对照等)。
25.另一方面，本发明提供了一种右眼等效球镜的检测系统，所述检测系统包括使用上述snp位点组合的检测结果判断受试者右眼等效球镜的计算装置。
26.优选地，所述诊断系统还可以包括snp位点的检测装置，例如测序仪等。
27.另一方面，本发明提供了使用以上snp位点组合构建检测右眼等效球镜模型的方法。
28.优选地，所述方法使用的是逻辑回归的方法。
29.优选地，本发明还包括10倍交叉验证的步骤。所述“10倍交叉验证”具体指利用sklearn.model_selection函数，选择0到100中的10个随机数字作为随机种子，分别将样本划分为10组训练集和验证集，进行10倍交叉验证，取10次结果的平均数作为模型结果。
30.另一方面，本发明还提供了一种检测右眼等效球镜的模型，所述模型根据前述snp位点组合的检测结果判断出受试者当下右眼等效球镜。
31.优选地，所述模型是经前述方法所构建得到。
32.优选地，所述模型可以是公式、诺莫图、或其他方便受试者操作的方式；根据此模型，可以直接计算得到受试者当下的右眼等效球镜。
33.另一方面，本发明提供了检测前述snp位点组合的试剂、上述试剂盒、系统、模型在制备检测右眼等效球镜的产品中的应用。
34.优选地，所述检测上述snp位点组合的试剂包括但不限于以下方法检测snp时所使用的试剂：taqman探针法、测序法、芯片法、飞行质谱仪(maldi-tofms)检测、限制性片段长度多态性法(pcr-rflp)、单链构象多态性法(pcr-sscp)、等位基因特异性pcr(as-pcr)、snapshot法、snplex法、变性高效液相色谱法(dhplc)、变性梯度凝胶电泳法(dgge)。本领域技术人员可选择任一种或几种方法来检测snp位点，只要可以实现snp位点的检测。
35.具体地，所述试剂包括但不限于引物、探针、芯片等。
36.另一方面，本发明提供了检测右眼等效球镜的方法，所述方法包括根据本发明所述snp位点组合的检测结果判断受试者右眼等效球镜的步骤，
37.具体地，所述方法可以包括以下步骤：
38.1)收集受试者样本，优选地，所述样本是口腔拭子(口腔黏膜样本)；
39.2)对样本进行snp检测，优选地，所述检测还可以包括提取dna的步骤；
40.3)根据2)的检测结果，判断受试者的右眼等效球镜。
41.所述判断可以是手动地、自动地、或它们组合地来执行或完成所选任务；可以根据检测结果手动计算结果，或输入前述系统自动地得到计算结果。
42.本发明所述“等效球镜(spherical equivalent refraction，se)”等于球镜度数 散光度数*1/2。屈光度≤-0.50d可作为近视诊断的循证共识阈值。更具体地，近视的定义是当调节放松时，眼睛的等效球镜屈光不正≤-0.5d的情况。高度近视的定义是当调节放松
时，眼睛的等效球镜屈光不正≤-6.00d的情况。低度近视的定义是当调节放松时，眼睛的等效球镜屈光不正≤-0.5且》-6.00d的情况。
43.本发明所述单核苷酸多态性(single nucleotide polymorphisms，snp)是指在基因组上单个核苷酸的变异，包括转换、颠换、缺失和插入，形成的遗传标记，其数量很多，多态性丰富。在人类基因组中大概每1000个碱基就有一个snp，人类基因组上的snp总量大概是3
×
106个。因此，snp成为第三代遗传标志，人体许多表型差异、对药物或疾病的易感性等等都可能与snp有关。
44.在本发明的整篇说明书中使用的术语“检测”意指判定受试者当下的右眼等效球镜。
45.具体地说，本发明解决问题的技术方案包括：
46.(1)建立统一标准的标本库和数据库：以标准操作程序(sop)采集符合标准的血液样本，系统收集完整的人口学资料和临床资料；
47.(2)基因型检测：纳入符合规定的样本，利用高密度snp芯片，在全基因范围内找出与右眼等效球镜相关的snp；
48.(3)对筛选出的阳性关联snp，进一步在另外的样本中进行检测，以验证其应用于临床诊断的可重复性和稳定性；
49.(4)右眼等效球镜检测试剂盒的研制：根据与右眼等效球镜有显著关联的snp开发检测试剂盒。
附图说明
50.图1本发明的所筛选到的snp组合在不同阈值下判断右眼等效球镜时的auc值变化曲线。
具体实施方式
51.下面结合实施例对本发明做进一步的说明，以下所述，仅是对本发明的较佳实施例而已，并非对本发明做其他形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更为同等变化的等效实施例。凡是未脱离本发明方案内容，依据本发明的技术实质对以下实施例所做的任何简单修改或等同变化，均落在本发明的保护范围内。
52.实施例1样本的收集、测序及数据分析
53.样本的纳入标准：
54.(1)最差眼的等效球镜小于等于-6.00d；
55.(2)年龄6岁-18岁；
56.(3)无其他眼部疾病、先天性遗传病、其他全身性疾病和身体异常；
57.(4)没有眼部手术史；
58.(5)中国人；
59.(6)自愿参加且签署知情同意书。
60.本发明是对2019年-2020年与温州医科大学附属眼视光医院与温州市政府合作开展的温州市全市中小学生近视筛查和干预项目的样本信息进行了整理，召集了符合上述纳入标准的10348(9852
‑‑
8961)样本进行全外显子组的测序，记录他们两眼的球镜(sph)、柱
dna technologies,san diego,usa)的标准流程，对预文库进行液相杂交捕获操作。
78.5、获得外显子文库：捕获产物洗脱回收后，进行pcr扩增和纯化即获得外显子文库。文库用qpcr方法进行定量，并使用agilent 2100对条带大小做检测。
79.6、illumina novaseq6000测序：使用illumina novaseq6000测序系统，对外显子文库进行150pe测序，原始图像使用casava v1.82软件进行碱基识别，生产原始测序数据。
80.7、比对人类参考基因组：使用burrows-wheeler aligner(bwa)工具将测序片段比对到人类参考基因组(ucsc hg19)，使用picard v1.57去除pcr重复序列。采用软件gatk进行变异检测，并对测序深度、覆盖深度及均一性等进行统计。
81.8、数据质量控制
82.对于样本：删除口腔拭子取样量不足，检出率小于90％，平均覆盖度小于10，平均基因型存在概率小于65％，基因型杂合率平均值偏离
±
4标准差(sd)，染色体异常、性别异常，具有亲缘关系和非东亚人群的样本。对于snp：将未通过vqsr质量控制的全外显子组检测结果，基因型检出率小于90％，哈代温伯格p值小于1e-06，变异allele数目(ac)为0的snp。
83.9、统计分析方法
84.通过snptest中的贝叶斯检验(bayesian tests)开展高度近视全外显子组数量性状研究，鉴定潜在的与高度近视临床状态显著相关的snp和基因。数量性状的bayesian tests是通过使用期望基因型(-method expected)的线性模型的共轭先验公式进行的。对于加性模型，公式为：
85.yi＝βgi ei,ei～n(0,σ2)。
86.其中yi是第i个个体的残差表型。残差表型是通过减去一个基线项和估计任何指定的协变量来计算的。gi为第i个个体的预期基因型的加性编码，野生纯合型＝“0”，杂合型＝“1”，突变纯合型＝“2”。
87.σ2是模型的误差方差。并且σ2～ig(a,β),β～n(m
β
,v
β
σ2)，默认的a＝3,β＝2,m
β
＝0,v
β
＝0.02。将年龄和性别作为协变量进行校正。
88.实验结果
89.1、snps选择：
90.利用snptest获取了8961个高度近视样本中共89095个常见变异与数量性状右眼等效球镜的全基因组关联关系。
91.通过4种阈值分别筛选snps进行建模：log10(bf)》3；log10(bf)》2；log10(bf)》1.5；log10(bf)》1，共获得342个snps。
92.2、原始snps基因型谱获取：
93.利用python提取vcf中对应所需snps(共342个snps)，选择加性模型(未突变设置为0，突变一个等位为1，突变两个等位为2)，获取原始的snp基因型谱。
94.3、snps基因型谱补缺失：
95.计算原始snps基因型谱中每个样本中snps的缺失数目，如果缺失》5％(n_缺失/n_allsamples)，去除对应的样本。对于缺失情况《5％的样本，利用人工神经网络进行补缺失。
96.人工神经网络共三层，包括输入层，隐藏层以及输出层：输入层节点数为len(label_data)，也就是完全没有缺失的snp；隐藏层设置20个节点；输出层设置3个节点；设
置学习率为0.003，训练迭代次数设置为20次。
97.4、分类模型：logistic regression
98.1)根据数量性状右眼等效球镜，取不同阈值，分别建立模型，并考察auc值，取auc最大的模型的阈值作为最终分类模型的阈值。
99.规则如下：
100.①
遍历右眼等效球镜的所有值，分别作为阈值将样本分为两组；
101.②
每组样本数量不得少于100
102.2)利用逻辑回归对样本进行分类，参数如下：
[0103][0104]
尽管它的名字是逻辑回归，但它是一种用于分类而不是回归的线性模型。
[0105]
逻辑回归默认情况下会应用正则化。这才机器学习中很常见，其一个优点就是提高了数值的稳定性，没有正则化相当于把c设置为非常高的值。
[0106]
作为一个优化问题，二进制类l2惩罚逻辑回归最小化以下成本函数：
[0107][0108]
相似的，l1正则化逻辑回归解决以下优化问题：
[0109][0110]
利用sklearn.model_selection函数，选择0到100中的10个随机数字作为随机种子，分别将样本划分为10组训练集和验证集，进行10倍交叉验证，取10次结果的平均数作为模型结果，以保证模型的客观性。
[0111]
各模型的cutoff值及其auc值如下：
[0112]
表1、各模型在不同cutoff值下的auc变化
[0113]
[0114][0115]
以logbf》3为标准筛选到的snp是rs2547319；
[0116]
以logbf》2为标准筛选到的snp是以下26个：
[0117]
rs10889315、rs75144949、rs149213974、rs754615、rs45583335、rs2529438、rs2072236、rs2240403、rs78155900、rs114097201、rs2811736、rs2297722、rs2297723、rs117519669、rs7157977、rs10143899、rs3759779、rs143477571、rs2303221、rs34693726、rs3803752、rs142921938、rs3746295、rs2547319、rs2547318、rs4911402；
[0118]
以logbf》1.5为标准筛选到的snp是以下102个：
[0119]
rs4648600、rs1061624、rs2275365、rs35909785、rs3103778、rs2457829、rs10889315、rs78587889、rs12568050、rs12564283、rs528123098、rs139842833、rs74953908、rs1077827、rs12492484、rs75144949、rs3732755、rs3135115、rs16837029、rs149213974、rs17001890、rs56382032、rs12651031、rs2736100、rs7724759、rs754615、rs5369、rs150157174、rs130068、rs45583335、rs78342561、rs2529438、rs1049402、rs2072236、rs2527878、rs2240403、rs78155900、rs114097201、rs2680903、rs59612871、rs3818584、rs16909677、rs10983755、rs2248077、rs2811736、rs2297722、rs2297723、rs9987876、rs7904554、rs907609、rs117519669、rs3750996、rs3802799、rs2276296、rs139692587、rs904628、rs11551723、rs144796593、rs78648016、rs7157977、rs10143899、rs3759779、rs7142098、rs143477571、rs3736911、rs11849022、rs77481241、rs3212112、rs9282879、rs55958706、rs3212056、rs34100926、rs114770841、rs3742368、rs118146919、rs2289702、rs2303221、rs34693726、rs3744566、rs9907420、rs375644567、rs185576254、rs3803752、rs17822735、rs17222523、rs34818467、rs12946397、rs3809724、rs181969864、rs7236574、rs142921938、rs3746295、rs151309111、rs533816037、rs2302948、rs182509214、rs2547319、rs2547318、rs2281209、rs4911402、rs138894461、rs5756223；
[0120]
以logbf》1为标准筛选到的snp是以下342个：
[0121]
rs16824585、rs4648600、rs17033266、rs1061624、rs202140107、rs2291804、rs2275365、rs139241751、rs35909785、rs6676052、rs3103778、rs116893711、rs116861599、rs2457829、rs10889315、rs17855078、rs78587889、rs140141416、rs4656197、rs45445497、rs143355122、rs12568050、rs12564283、rs10482、rs1929861、rs75401237、rs28517482、rs1264208208、rs2631840、rs188480146、rs528123098、rs61732269、rs367931009、rs139842833、rs79908358、rs79358798、rs1434087、rs4894048、rs13421990、rs13419085、rs143196218、rs139713392、rs6734083、rs150887565、rs4663795、rs74953908、rs143001518、rs6437368、rs147838565、rs1077827、rs144623840、rs78233496、rs9876921、rs6622、rs2279909、rs547922447、rs300977、rs12492484、rs75144949、rs3732755、rs3135115、rs16837029、rs74787220、rs527903313、rs149213974、rs17001890、rs79769348、rs10025654、rs76996680、rs2306369、rs79482316、rs56382032、rs12651031、rs2736100、rs35130836、rs3749684、rs7724759、rs754615、rs45625534、rs3805625、rs137898880、rs147518598、rs202007699、rs2270900、rs75164992、rs17133512、rs76720430、rs5369、rs150157174、rs7756481、rs117298661、rs117644703、rs3094672、rs4713420、rs130068、rs67196728、rs185990098、rs9349593、rs143580489、rs3757302、rs3218602、rs9495370、rs45583335、rs78342561、rs117358940、rs9689983、rs2529438、rs1049402、rs2072236、rs2527878、rs2240403、rs781980620、rs56083327、rs2245368、rs2430307、rs3748126、rs36028884、rs78155900、rs114097201、rs13438494、rs117211799、rs115734214、rs148474510、rs200506987、rs8190955、rs8178175、rs2680903、rs117676215、rs115507207、rs369806130、rs78240711、rs59612871、rs117169590、rs3739523、rs145293869、rs139212809、rs3818584、rs16909677、rs141549766、rs147026086、rs188197626、rs13321、rs10983755、rs2301576、rs7038042、rs11849、rs2275161、rs2248077、rs2811736、rs2297722、rs2297723、rs9987876、rs7904554、
rs148582275、rs3740168、rs149647444、rs7100474、rs7079648、rs202166096、rs60420182、rs36104328、rs7098255、rs3814163、rs2297785、rs3736924、rs3736923、rs141296329、rs117520659、rs57285449、rs112687793、rs73401681、rs79948463、rs10902170、rs10902171、rs867839、rs907608、rs907609、rs117519669、rs3750996、rs214086、rs3802799、rs76215382、rs147852038、rs2276296、rs143930161、rs183113103、rs139692587、rs10792769、rs7102675、rs116878689、rs35466364、rs142908193、rs7966459、rs7135745、rs34181394、rs1921038、rs201814780、rs11068531、rs904628、rs11551723、rs7983667、rs144796593、rs140526815_rs77878407、rs78646508、rs34523940、rs9521906、rs78648016、rs7157977、rs10143899、rs3759779、rs7142098、rs41309252、rs535020964、rs3742770、rs3742768、rs116881452、rs143477571、rs61992274、rs142973720、rs530110551、rs3742464、rs3736911、rs11849022、rs77481241、rs3212112、rs9282879、rs55958706、rs3212056、rs34100926、rs114770841、rs3742368、rs144857261_rs373032548、rs189734192、rs188393827、rs118146919、rs562829514、rs1805025、rs8192352、rs2289702、rs150921132、rs61869621、rs1802752、rs113355794、rs7198064、rs2291407、rs9941128、rs9941210、rs9941215、rs9931527、rs12708649、rs13306653、rs11643517、rs5723、rs7196736、rs2303221、rs2304483、rs149645127、rs11550470、rs17232091、rs150889000、rs16954357、rs371756524_rs557159990、rs192103504、rs115314145、rs34693726、rs2277651、rs2277652、rs71360269、rs71360270、rs3744566、rs3744568、rs9907420、rs375644567、rs185576254、rs2158460、rs3803752、rs17822735、rs17222523、rs34818467、rs12946397、rs3809724、rs35476409_rs397978887、rs3744214、rs181969864、rs4789773、rs28561791、rs7236574、rs76597875、rs61749945、rs118033234、rs1052696、rs1052692、rs148018996、rs34488963、rs3745379、rs3745385、rs142921938、rs3746295、rs151309111、rs139271842、rs140566569、rs533816037、rs78237760、rs79621739、rs75990051、rs62107869、rs305974、rs1055099、rs41290128、rs76151738、rs16980091、rs2302948、rs182509214、rs2547319、rs2547318、rs16982743、rs2284385、rs2281209、rs4911402、rs80158178、rs7509972、rs150511237、rs2231391、rs2231390、rs762228、rs72564613、rs191574093、rs138894461、rs5999924、rs5750146、rs5756223、rs3747168、rs148345118、rs910796、rs910797、rs910798、rs738712、rs16445、rs3922872、rs56248708、rs9628315、rs8137790。
[0122]
再将以上模型进行优化，筛选最佳阈值(cutoff)及比绘制roc曲线(如表2、图1所示)
[0123]
表2、各模型最佳阈值及其auc值
[0124]
模型阈值aucsnp数量模型1-11.1250.59以logbf》3为标准筛选到的1个模型2-11.1250.67以logbf》2为标准筛选到的26个模型3-11.250.77以logbf》1.5为标准筛选到的102个模型4-11.1250.80以logbf》1为标准筛选到的342个
[0125]
也即，使用以logbf》3为标准筛选到的1个snp位点为标志物时，可有效判断受试者的右眼等效球镜是否在-11.125之上；
[0126]
以logbf》2为标准筛选到的26个为标志物时，可有效判断受试者的右眼等效球镜是否在-11.125之上；
[0127]
以logbf》1.5为标准筛选到的102个为标志物时，就可有效判断受试者的右眼等效球镜是否在-11.25之上；
[0128]
以logbf》1为标准筛选到的342个为标志物时，就可有效判断受试者的右眼等效球镜是否在-11.125之上。