一种用于检测高原适应性的SNP标志物及试剂盒的制作方法

一种用于检测高原适应性的snp标志物及试剂盒
技术领域
1.本发明涉及一种用于检测高原适应性的snp标志物及试剂盒，属于生物技 术领域。


背景技术：

2.近年来，随着高原地区经济的快速发展和旅游业的盛行，高原移居人口日 益增多。但是，高原的低氧环境给人们进入高原地区带来了挑战。非高原世居 人群，在进入高原低氧环境下，极易诱发急、慢性高原病。因此，预测个体高 原环境适应能力，将有望为高原地区的旅游、工作等提供科学指导。高原适应 性具有明显的遗传学基础。世居高原的藏族人群已适应了青藏高原上的低氧和 高寒等恶劣环境，并世代生存，将优势遗传基因代代相传，是世界高原人群中 居住海拔最高、适应历史最长的群体。人类基因组脱氧核糖核酸 (deoxyribonucleic acid，dna)序列中最常见的多态性包括单核苷酸多态性 (single nucleotide polymorphism，snp)，是机体表型差异和环境适应性的遗传 学基础。因此，有望基于snp信息，实现人体高原适应性的预测。
3.随着基于snp芯片技术、高通量测序技术和“中性理论检验”的分子进化 分析的逐步发展，人类高原适应性(high altitude adaptation，haa)相关基因 组学研究也在大步向前。多项报道称发现了epas1基因和egln1基因的多态性 是藏族人群能够适应高原地区低氧环境且免于罹患高原疾病的一个重要遗传学 基础。但rs1204168所在基因mcur1与高原适应性的关系及其作为高原适应性 的标志物尚未见报道，同时在已有的专利中尚无涉及rs1204168的应用。


技术实现要素：

4.为克服现有技术存在的缺陷，本发明的目的之一在于提供一种用于检测高 原适应性的snp标志物。
5.本发明的目的之二在于提供一种用于检测高原适应性的试剂盒。
6.为实现本发明的目的，提供以下技术方案。
7.一种用于检测高原适应性的snp标志物，所述snp标志物为rs1204168， rs号表示该位点在dbsnp数据库中的编号，基于人类基因组第19版本，rs1204168 位于人类第6号染色体13,836,672位置，所述rs1204168为cc或tc的基因型 是人体具有高原适应性的指示。
8.进一步地，sanger测序方法所使用的前向引物序列为核苷酸序列表中数字 标识符<210>所表示的序列标识符1所述的核苷酸序列，即核苷酸序列表中的 seq id no.1(简称seq id no.1)；后向引物序列为核苷酸序列表中的seq idno.2(简称seq id no.2)；利用所述前向引物和后向引物可以使用sanger测序 方法有效地检测本发明所述snp标志物rs1204168。
9.进一步地，sequenom基因分型方法所使用的前向引物序列为核苷酸序列表 中的seq id no.3(简称seq id no.3)，后向引物序列为核苷酸序列表中的seqid no.4(简称
seq id no.4)，延伸引物序列为核苷酸序列表中的seq id no.5 (简称seq id no.5)；利用所述前向引物、后向引物和延伸引物可以使用 sequenom基因分型方法有效地检测本发明所述snp标志物rs1204168。
10.一种用于检测高原适应性的试剂盒，所述试剂盒包括本发明所述的一种用 于检测高原适应性的snp标志物rs1204168，所述rs1204168为cc或tc的基 因型是人体具有高原适应性的指示。
11.进一步地，所述试剂盒还包括sanger测序方法所使用的序列为seq id no. 1的前向引物，以及序列为seq id no.2的后向引物；利用所述前向引物和后 向引物可以使用sanger测序方法有效地检测本发明所述snp标志物rs1204168。
12.进一步地，所述试剂盒还包括sequenom基因分型方法所使用的序列为seqid no.3的前向引物，序列为seq id no.4的后向引物和序列为seq id no.5 的延伸引物；利用所述前向引物、后向引物和延伸引物可以使用sequenom基因 分型方法有效地检测本发明所述snp标志物rs1204168。
13.进一步地，所述试剂盒的检测样本为被试人体的血液、皮肤或皮下组织。
14.有益效果
15.1.本发明提供了一种用于检测高原适应性的snp标志物，所述snp标志 物为rs1204168，可以将rs1204168为cc或tc的基因型作为人体具有高原适应 性的指示。
16.2.本发明提供了一种用于检测高原适应性的snp标志物，利用前向引物和 后向引物可以使用sanger测序方法有效地检测所述snp标志物rs1204168。
17.3.本发明提供了一种用于检测高原适应性的snp标志物，利用前向引物、 后向引物和延伸引物可以使用sequenom基因分型方法有效地检测所述snp标 志物rs1204168。
18.4.本发明提供了一种用于检测高原适应性的试剂盒，所述试剂盒包括所述 snp标志物rs1204168以及其使用sanger测序方法或sequenom基因分型方法检 测所使用的引物，可以通过所述试剂盒来检测rs1204168的基因型来有效地检测 人体是否具有高原适应性。
附图说明
19.图1为实施例1的藏族人群与汉族人群的f
st
结果曼哈顿图。
20.图2为实施例3中对实施例2的研究对象检测的rs1204168基因型频率图。
21.图3为实施例4的rs1204168基因型与血红蛋白浓度图。
具体实施方式
22.下面结合附图和实施例详细描述本发明，所述实施例的示例在附图中示出。 通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为 对本发明的限制。
23.若未特别指明，实施例中所采用的技术手段为本领域技术人员所熟知的常 规手段，可以参照《分子克隆实验指南》第三版或相关产品进行。所采用的试 剂盒产品也均为可商业获得的。未详细描述的各种过程和方法是本领域中公知 的常规方法，所用试剂的来源、商品名以及有必要列出其组成成分者，均在首 次出现时表明，其后所用相同试剂如无特殊说明，均与首次标明的内容相同。 实施例1全基因组进化分析研究
24.研究对象：48例藏族人和50例汉族人的基因组dna。48例藏族人的基因 组dna，分别来源于48个没有亲缘关系的藏族人的外周血，所述藏族人分别居 住在四川省甘孜藏族自治州海拔2,600米到4,200米的5个城市，青海慢性高原 病评分(qinghai chronic mountain sickness score)小于5分，健康无高原病。 50例汉族人的基因组dna，分别来源于50个没有亲缘关系的汉族人的外周血， 所述汉族人居住在广西省南宁市，平均海拔79米，具体信息如表1所示。
25.将所述研究对象通过全基因组测序检测rs1204168位点的基因型信息。
26.对所述研究对象进行进化分析研究，使用固定系数(fixation index，f
st
) 检验方法和基于加性模型的逻辑斯蒂(logistic)回归分析方法。其中，f
st
使用 vcftools软件计算，用于检验藏族人和汉族人的遗传距离。逻辑斯蒂回归以 rs1204168基因型为自变量，设定tt＝0，tc＝1，cc＝2，表型(民族)为因 变量，校正性别、年龄，使用plink软件1.07版本计算，用于检验rs1204168 在藏族人和汉族人之间的频率差异。
27.表1用于全基因组snp频率差异研究的人群信息总表
[0028][0029]
进化分析研究结果：
[0030]
全基因组平均f
st
可评估两个人群之间的遗传距离，本实施例发现藏族与汉 族之间的遗传距离，即全基因组平均f
st
是0.011。snp的f
st
值越大，提示这 些snp越可能受到正向选择(positive selection)。本实施例中取f
st
排名前1,500 的snp，作为潜在受到正向选择的snp阈值。这些snp的f
st
值(f
st
≥0.235) 均高于全基因组f
st
值平均值以上5个标准差(f
st
>0.146)。其中，除去排名第 一的snp，其f
st
值(f
st
＝0.565)高于平均值以上20个标准差(f
st
>0.551)， 新信号中排名第一的是6p23区域(index snp rs1204168，f
st
＝0.464)，该自然 选择信号区域包含14个f
st
非常显著的snps(f
st
>0.42)，并且，该区域及附 近基因尚未报道与高原适应性相关。由f
st
检验的结果显示，在大多数的snp 上，藏族人和汉族人的基因组遗传距离较小，如图1所示；但是在rs1204168上， 观测到rs1204168的f
st
值为0.464，提示
rs1204168是有强正向选择可能性的snp，即rs1204168的c等位基因与藏族人适应高原环境的正向选择相关的可能性大。
[0031]
逻辑斯蒂回归方法的结果显示，rs1204168是藏族人和汉族人的频率差异较大的snp，其统计学检验or＝0.08759，p＝1.11
×
10
‑7。藏族人携带的较高频率的变异是rs1204168的c等位基因。
[0032]
综合以上结果证明rs1204168的c等位基因与个体的高原适应性密切相关。
[0033]
实施例2利用sanger测序方法对研究对象的rs1204168进行基因分型
[0034]
研究对象：来自实施例1的48例藏族人基因组dna样本。利用3130xlgeneticanalyzer(abi，fostercity，ca)对所述dna样本的pcr扩增产物进行sanger测序，以便获得测序结果；基于测序结果，确定所述dna样本的rs1204168基因型。sanger测序方法中，用于针对rs1204168的基因分型的前向引物和后向引物如下：
[0035]
前向引物序列为：tatgagaccaccgagatttg(seqidno.1)，
[0036]
后向引物序列为：gccttcacattgtttccc(seqidno.2)。
[0037]
测序结果显示：所述dna样本分别呈现cc基因型和tt基因型，说明sanger测序可以有效地区分rs1204168的不同基因型。
[0038]
实施例3rs1204168在独立人群中的频率
[0039]
研究对象：672例藏族人和270例汉族人的基因组dna样本。藏族人均世居西藏拉萨市，平均海拔3,650米。汉族人来自全国多个省市，海拔均在2,500米以下。所有个体间均不存在亲缘关系，具体信息见表2。所有个体均已签署知情同意书，且本研究经辐射医学研究所伦理委员会批准执行。
[0040]
针对rs1204168设计引物，通过sequenom检测获得rs1204168突变情况。sequenom检测的具体步骤如下：
[0041]
针对所述研究对象的基因组dna样本中待分型的rs1204168，设计pcr引物以及延伸引物。这一过程使用massarraydesign3.0软件实施。每次实验每个样本需使用约15ng的基因组dna。使用pcr实验进行dna模板扩增，扩增产物和延伸引物被用于snp位点的特异性单碱基延伸反应。将反应产物经自动点样仪转移至spectrochip芯片(sequenom，加利福尼亚州，美国)。在sequenom核酸质谱平台(sequenom，加利福尼亚州，美国)进行基因型检测，下机数据，使用massarraytyper软件(sequenom，加利福尼亚州，美国)进行分析。sequenom检测中，用于针对rs1204168基因分型的引物如下：
[0042]
前向引物：acgttggatgtgtgaatggtggcccttcag(seqidno:3)；
[0043]
后向引物：acgttggatgaagtgtgtttctggcttctc(seqidno:4)；
[0044]
延伸引物：ggcttctctcttggaaag(seqidno:5)。
[0045]
表2用于实施例3的人群信息表
[0046][0047]
结果：通过计数可得rs1204168的c等位基因频率在藏族人群中为75％， 在汉族人群中为47％。即，rs1204168的c等位基因频率在藏族人群和汉族人群 中存在差异，如图2所示。
[0048]
实施例4rs1204168的基因型与血红蛋白浓度的相关性
[0049]
血红蛋白浓度指单位体积(l)血液内所含血红蛋白的量。血红蛋白是一种 含色素的结合蛋白质，是红细胞主要成分，能与氧结合，运输氧和二氧化碳。 血红蛋白浓度升高分为生理性和病理性。前者常见于住在高原地区的人，由于 高原地区相对来说氧浓度较低，所以长期住在那里的人，机体会代偿性生成较 多的血红蛋白。但血红蛋白浓度增减与红细胞浓度增减一致，当机体失代偿后 血红蛋白浓度超过阈值则会导致机体罹患红细胞增多症，即血红蛋白浓度的过 度升高与高原疾病的发生密切相关。
[0050]
研究对象：使用的是实施例3中的270个人的基因组dna样本。针对 rs1204168设计引物，通过sequenom检测获得rs1204168的基因型。血红蛋白浓 度，使用hemocue hb 201

分析仪(恩厄尔霍尔姆，瑞典)检测。rs1204168的 基因型与血红蛋白浓度的相关性，采用皮尔森相关性分析评估时，rs1204168基 因型中，设定tt＝0，tc＝1，cc＝2。当p<0.05时认为具有显著的统计学意 义。
[0051]
结果：为了探索rs1204168能否作为生物标记物，使用rs1204168的c等位 基因频率与血红蛋白浓度进行了相关性分析。发明人证明了rs1204168的c等位 基因频率与血红蛋白浓度呈显著负相关：r＝
‑
0.141，p＝0.0199，如图3所示， 图中n代表人数。说明c等位基因频率越高，血红蛋白浓度越低，进一步说明 如个体的rs1204168位点的基因型为cc或者tc，则认为其具有一定的高原适 应性。