CN108642160B - 检测胎儿地中海贫血致病基因的方法和试剂盒 - Google Patents
检测胎儿地中海贫血致病基因的方法和试剂盒 Download PDFInfo
- Publication number
- CN108642160B CN108642160B CN201810473083.4A CN201810473083A CN108642160B CN 108642160 B CN108642160 B CN 108642160B CN 201810473083 A CN201810473083 A CN 201810473083A CN 108642160 B CN108642160 B CN 108642160B
- Authority
- CN
- China
- Prior art keywords
- thalassemia
- dna
- gene
- detection
- sequencing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 208000002903 Thalassemia Diseases 0.000 title claims abstract description 44
- 108090000623 proteins and genes Proteins 0.000 title claims abstract description 27
- 230000001605 fetal effect Effects 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 title abstract description 18
- 230000001717 pathogenic effect Effects 0.000 title abstract description 17
- 238000002360 preparation method Methods 0.000 claims abstract description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 4
- 238000003745 diagnosis Methods 0.000 claims 1
- 210000003754 fetus Anatomy 0.000 abstract description 33
- 230000008774 maternal effect Effects 0.000 abstract description 29
- 238000012163 sequencing technique Methods 0.000 abstract description 29
- 230000008775 paternal effect Effects 0.000 abstract description 28
- 238000001514 detection method Methods 0.000 abstract description 26
- 238000005516 engineering process Methods 0.000 abstract description 21
- 239000000523 sample Substances 0.000 abstract description 20
- 230000002068 genetic effect Effects 0.000 abstract description 18
- 230000035772 mutation Effects 0.000 abstract description 15
- 238000012165 high-throughput sequencing Methods 0.000 abstract description 11
- 210000004369 blood Anatomy 0.000 abstract description 8
- 239000008280 blood Substances 0.000 abstract description 8
- 238000013461 design Methods 0.000 abstract description 8
- 238000012216 screening Methods 0.000 abstract description 8
- 206010064571 Gene mutation Diseases 0.000 abstract description 5
- 108020004414 DNA Proteins 0.000 description 81
- 102000054766 genetic haplotypes Human genes 0.000 description 49
- 238000010276 construction Methods 0.000 description 18
- 239000000047 product Substances 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 9
- 238000012252 genetic analysis Methods 0.000 description 9
- 239000000243 solution Substances 0.000 description 9
- 208000005980 beta thalassemia Diseases 0.000 description 8
- 210000005259 peripheral blood Anatomy 0.000 description 8
- 239000011886 peripheral blood Substances 0.000 description 8
- 201000006288 alpha thalassemia Diseases 0.000 description 7
- 238000003766 bioinformatics method Methods 0.000 description 6
- 239000003153 chemical reaction reagent Substances 0.000 description 6
- 238000000746 purification Methods 0.000 description 6
- 238000002864 sequence alignment Methods 0.000 description 6
- 238000001617 sequential probability ratio test Methods 0.000 description 6
- 101150071666 HBA gene Proteins 0.000 description 5
- 102100021519 Hemoglobin subunit beta Human genes 0.000 description 5
- 238000002669 amniocentesis Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000011144 upstream manufacturing Methods 0.000 description 5
- 238000007400 DNA extraction Methods 0.000 description 4
- 230000003321 amplification Effects 0.000 description 4
- 230000002759 chromosomal effect Effects 0.000 description 4
- 239000012634 fragment Substances 0.000 description 4
- 238000003199 nucleic acid amplification method Methods 0.000 description 4
- 101150013707 HBB gene Proteins 0.000 description 3
- 238000012408 PCR amplification Methods 0.000 description 3
- 206010043391 Thalassaemia beta Diseases 0.000 description 3
- 238000003793 prenatal diagnosis Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 102100027685 Hemoglobin subunit alpha Human genes 0.000 description 2
- 108091005904 Hemoglobin subunit beta Proteins 0.000 description 2
- 101001009007 Homo sapiens Hemoglobin subunit alpha Proteins 0.000 description 2
- 210000000349 chromosome Anatomy 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 238000007403 mPCR Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000002773 nucleotide Substances 0.000 description 2
- 125000003729 nucleotide group Chemical group 0.000 description 2
- 238000009609 prenatal screening Methods 0.000 description 2
- 206010000234 Abortion spontaneous Diseases 0.000 description 1
- 108700028369 Alleles Proteins 0.000 description 1
- 206010008805 Chromosomal abnormalities Diseases 0.000 description 1
- 208000031404 Chromosome Aberrations Diseases 0.000 description 1
- 208000032170 Congenital Abnormalities Diseases 0.000 description 1
- 206010016852 Foetal damage Diseases 0.000 description 1
- 108091092584 GDNA Proteins 0.000 description 1
- 206010018910 Haemolysis Diseases 0.000 description 1
- 208000026350 Inborn Genetic disease Diseases 0.000 description 1
- 206010056254 Intrauterine infection Diseases 0.000 description 1
- 108091028043 Nucleic acid sequence Proteins 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 210000004381 amniotic fluid Anatomy 0.000 description 1
- 208000007502 anemia Diseases 0.000 description 1
- 208000036878 aneuploidy Diseases 0.000 description 1
- 231100001075 aneuploidy Toxicity 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001574 biopsy Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007698 birth defect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000002939 deleterious effect Effects 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 238000007847 digital PCR Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 238000011049 filling Methods 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 102000054767 gene variant Human genes 0.000 description 1
- 208000016361 genetic disease Diseases 0.000 description 1
- 230000008588 hemolysis Effects 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 210000001161 mammalian embryo Anatomy 0.000 description 1
- 238000004949 mass spectrometry Methods 0.000 description 1
- 208000015994 miscarriage Diseases 0.000 description 1
- 230000036438 mutation frequency Effects 0.000 description 1
- 230000000869 mutational effect Effects 0.000 description 1
- 238000001821 nucleic acid purification Methods 0.000 description 1
- 108020004707 nucleic acids Proteins 0.000 description 1
- 102000039446 nucleic acids Human genes 0.000 description 1
- 150000007523 nucleic acids Chemical class 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 108090000765 processed proteins & peptides Proteins 0.000 description 1
- 239000012264 purified product Substances 0.000 description 1
- 210000000278 spinal cord Anatomy 0.000 description 1
- 208000000995 spontaneous abortion Diseases 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000002054 transplantation Methods 0.000 description 1
Images
Classifications
-
- 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/6869—Methods for sequencing
-
- 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
-
- 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/156—Polymorphic or mutational markers
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Organic Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Genetics & Genomics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Immunology (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Biotechnology (AREA)
- Biochemistry (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种检测胎儿地中海贫血致病基因的方法，包括如下步骤：(1)筛选SNP位点，SNP位点用于设计扩增基因组地贫基因的引物池和捕获孕妇血浆中游离DNA地贫基因的探针；(2)提取孕妇血浆中游离DNA和父‑母‑已出生同胞三人的全血基因组DNA并构建相应的DNA文库，进行模板制备和富集；(3)对步骤(2)游离DNA和全血基因组DNA文库测序；(4)构建地贫基因上SNP位点的单倍体基因型，结合游离DNA和全血基因组DNA文库测序信息，分析父源遗传情况和母源遗传情况，以确定胎儿SNP位点相应的基因型。本发明通过目标区域捕获结合高通量测序技术，实现对地贫无创产前检测；所需的样本量少；该方法在对父源突变检测的基础上，可以实现对胎儿母源基因突变的检测。
Description
技术领域
本发明涉及基因检测技术领域，尤其是一种检测胎儿地中海贫血致病基因的方法和试剂盒。
背景技术
地中海贫血(简称地贫)是全世界的高发遗传病之一，是由于人体α，β-珠蛋白基因突变或者缺失而导致α，β-珠蛋白肽链合成速率的不平衡，从而引起的溶血性贫血。地贫常见的两种类型是α-地贫和β-地贫，α-地贫相关基因为HBA2和HBA1，β-地贫相关基因为HBB。地贫集中在热带和亚热带地区，多发生于地中海沿岸国家，其次为中东、印度、巴基斯坦、东南亚、中国南方和北非，美国是移民国家，其发生率也较高。我国长江以南各省是地贫高发区，广东、广西、海南、台湾和香港等地受累人口超过2亿。目前地贫尚无根治的方法，只能使用输血治疗或脊髓移植，其价格昂贵且彻底治愈几率不高，因此产前筛查和诊断对预防该种出生缺陷疾病具有重大意义。
目前临床报道产前检测地贫基因突变的方法有几种：传统的羊膜穿刺、腹绒毛活检等，这些有创伤性的操作可能导致胎儿损伤、流产或宫内感染等。1997年，卢煜明发现孕妇外周血浆中存在来自于胎儿的游离DNA，这一发现为通过孕妇外周血无创产前诊断和检测胚胎基因型提供了新的可能。已有报道利用孕妇血浆游离DNA与质谱、数字PCR、COLD-PCR等技术对地贫已成功进行无创检测。但这些方法是基于单一位点的PCR技术，由于血浆DNA含量低及高度片段化的特点，很可能会带来PCR的假阴性。除此以外，这些方法只能够检测与母亲不同的特异性父源突变是否存在，而无法进一步确定母源突变的存在情况。
随着高通量测序技术的应用，孕妇血浆DNA无创产前检测胚胎染色体异常(如21、18、13染色体非整倍体变异)以及部分大的拷贝数变异已经成功应用于临床产前筛查和诊断。进一步地，通过高通量测序技术实现检测胎儿的家系中(指胎儿的父亲和母亲)基因型和单倍体型后，对孕妇血浆游离DNA的目标区域高通量测序，通过RHDO和SPRT等生物信息学方法实现对胎儿是否携带地贫相关基因变异的检测。高通量测序技术具有准确度高、通量大、样本需求量少、可以产生大量数据的特点。
发明内容
基于上述问题，本发明的目的在于克服上述现有技术的不足之处而提供一种高精确度的无创产前胎儿地中海贫血致病基因的检测方法，该方法不仅能够检测与母亲不同的特异性父源突变是否存在，还可以进一步确定母源突变的是否存在。
为实现上述目的，本发明采取的技术方案包括以下几个方面：
作为本发明的第一个方面，本发明提供一种检测胎儿地中海贫血致病基因的方法，包括如下步骤：(1)筛选SNP位点，所述SNP位点用于设计扩增基因组地贫基因的引物池和捕获孕妇血浆中游离DNA地贫基因的探针；(2)提取孕妇血浆中游离DNA和父-母-已出生同胞(指胎儿的姐姐或哥哥)三人的全血基因组DNA并构建相应的DNA文库，进行模板制备和富集；(3)对步骤(2)所述游离DNA和全血基因组DNA文库测序；(4)构建所述地贫基因上SNP位点的单倍体基因型，结合所述游离DNA和全血基因组DNA文库测序信息，分析父源遗传情况和母源遗传情况，以确定胎儿所述SNP位点相应的基因型。
应当说明的是，该方法基于多重PCR和目标区域捕获技术，能对α&β-地中海贫血进行无创产前检测，即针对α-地贫的设计区域(chr16:60022-2233661)：下游(10Kb)+HBA基因+上游(10Kb)，以及针对β-地贫血设计区域(chr11:3236690-7177304)：下游(10Kb)+HBB基因+上游(10Kb)，利用dbSNP数据库筛选出SNP位点195个α-地贫SNP位点和275个β-地贫SNP位点，根据筛选的SNP位点分别设计地贫引物与探针，通过多重PCR技术对父-母-已育小孩基因组DNA进行文库构建，在通过目标区域捕获技术对母亲游离DNA进行文库构建；将所有构建好的文库进行高通量测序，将得到的测序数据与参考序列人类基因组hg19作比对，通过扩增基因组DNA地贫基因的引物池位点信息进行序列比对，利用GATK(3.4版本)得到父-母-已育小孩三人核心家系基因组DNA的SNP位点，以测序深度大于30×(倍)，测序碱基质量值大于Q20为筛选条件，筛选多态位点；再根据孟德尔遗传定律，把筛选得到的SNP基因型信息构建染色体单倍体型，并结合家系样本表型和已知突变型信息将致病突变关联到单倍体型上；通过与捕获游离DNA地贫基因的探针位点信息进行序列比对，得到母亲外周血游离DNA各SNP位点等位碱基的频率，通过父源和母源生物信息学分析胎儿遗传的父母单倍体型，从而确定胎儿基因型。
优选地，所述步骤(2)和(3)中涉及DNA提取试剂，DNA纯化试剂，DNA文库构建试剂，DNA高通量测序试剂；DNA提取体系中提取游离DNA的血浆量不少于2mL，DNA纯化试剂为AgencourtAMPure XP-核酸纯化试剂盒。DNA文库构建试剂中基因组DNA模板不少于25ng，捕获血浆游离DNA文库模板不少于100ng；高通量测序体系选用Ion Torrent Proton测序平台进行测序。
优选地，所述父源遗传分析中，筛选父亲杂合母亲纯合的基因型位点和父亲和母亲均为纯合且基因型不同的位点，在上述两类位点中，再筛选出孕妇血浆在上述位点中为杂合的位点，进行父源遗传分析。
优选地，所述母源遗传分析中，在父源遗传情况分析基础上，筛选出父亲纯合母亲杂合的基因型位点，根据单倍型相对剂量(relative haplotype dosage，RHDO)分析结合序贯概率比试验(sequential probability ratio test，SPRT)进行母源遗传分析。
优选地，所述SNP位点的筛选条件为：在中国人千人基因组南方人和北方人数据库中同时MAF>＝20％，SNP位点上下100bp序列是特异区域，在基因组上没有同源性，并且45％<GC<70％，连续3个碱基不同。
更优选地，所述SNP位点信息见表15。
优选地，所述步骤(4)中，单倍体基因型的构建、父源遗传情况和母源遗传情况分析均使用NCBI数据库的hg19基因作为参考序列。
优选地，所述步骤(4)中，单倍体基因型的构建包括以下步骤：
1)构建父-母-已育小孩外周血基因组DNA的文库：
包括目的基因(即HBA、HBB)的多重PCR扩增；引物的部分消化；扩增产物的接头连接；连接产物的纯化；最后的文库片段PCR扩增及纯化；2)将上述构建的文库进行高通量测序，将人类基因组hg19作为参考序列，通过与扩增基因组DNA地贫基因的引物池位点信息进行序列比对，利用GATK(3.4版本)得到父-母-已育小孩三人核心家系基因组DNA的SNP位点，以测序深度大于30×，测序碱基质量值Q大于20为筛选条件，筛选多态位点，再根据孟德尔遗传定律，把筛选得到的SNP基因型信息构建染色体单倍体型，并结合家系样本表型和已知突变型信息将致病突变关联到单倍体型上。
其中，碱基质量值，简称Q值，是碱基读取错误率p的取整映射结果；目前，Q值与p值的对应关系有两种计算方式：标准Sanger变体计算方式以及早期Solexa管道(如IlluminaGenome Analyzer软件)计算方式。以标准Sanger变体方式为例，Q值被称为Phred qualityscore，它与p值之间的对应关系用以下公式计算得到：测序错误率与Q值的简明对应方式如下：Q13的碱基错误率为5％，Q20的碱基错误率为1％，Q30的碱基错误率为0.1％，Q40的碱基错误率为0.01％；地贫扩增引物池，包括了扩增HBA基因和其上下游2Mb染色体区域的1354对引物和扩增HBB基因和其上下游2Mb染色体区域的361对引物；应当说明的是，1354对引物和361对引物扩增出的目的片段都覆盖了相应的SNP位点(SNP位点信息参见表15)，为了获得足够的测序深度(30×)，所以设计了一共1715对引物，当为了获得更大的测序深度时，可以设计更多对引物，只要引物扩增出的目的片段覆盖表15的SNP位点即可；本领域普通技术人员在得知表15的SNP位点信息后，能够根据NCBI数据库中的HBA和HBB核苷酸序列以及引物设计网站或引物设计软件设计出所需的引物对(其中包括所述的1715对引物)。
优选地，所述步骤(4)中，父源遗传情况分析包括如下步骤：
1)构建母体外周血游离DNA的文库：
包括游离DNA的末端补平；接头连接；连接产物的扩增及纯化；纯化后产物的浓缩和探针捕获杂交；捕获后产物的富集；最后的富集产物PCR扩增及纯化；
2)将上述构建的文库进行高通量测序，将人类基因组hg19作为参考序列，通过捕获游离DNA地贫基因的探针位点信息进行序列比对，得到母亲外周血游离DNA各SNP位点等位碱基的频率，其中，挑选基因型为父亲杂合(MN，M和N代表碱基A,T,C,G)母亲纯合(MM)的位点和父亲和母亲均为纯合且基因型不同的位点(父亲MM，母亲NN)，在上述两类位点中，再筛选出孕妇血浆基因型为MN的位点，以此确定胎儿是否遗传了父源致病单倍体型。
优选地，所述步骤(4)中，母源遗传情况分析包括如下步骤：在父源遗传情况分析基础上，挑选出基因型为父亲纯合(MM)母亲杂合(MN)的位点，根据RHDO-SPRT确定胎儿是否遗传了母源的致病单倍体型，从而推断胎儿是否遗传了有害的突变型以及是否致病。
优选地，所述SPRT曲线根据以下公式计算：
其中，
所述Upper boundary与Lower bundary分别指的是上界和下界；在此检测方法中，把单倍体型1定义为上界，把单倍体型2定义为下界；根据SPRT的算法，计算筛选出的SNP位点累积的深度和累积的变异频率；如果筛选出的SNP位点都在上下界之间，那么，则不能分辨出胎儿遗传了母亲的哪一条单倍体型；如果累积的SNP位点中有一个或一个以上超过了上界，则可以判断胎儿遗传了母亲的单倍体型1；如果累积的SNP位点中有一个或一个以上超过了下界，则可以判断胎儿遗传了母亲的单倍体型2；
N是用于此分类的所有SNP位点的统计数；q0是胎儿遗传了母亲单倍体2的情况下，统计出reads来自于单倍体型1的比例；q1是胎儿遗传了母亲单倍体1的情况下，统计出reads来自于单倍体型1的比例。
作为本发明的第二个方面，本发明提供一种检测胎儿地中海贫血致病基因的试剂盒，所述试剂盒包括上述扩增基因组地贫基因的引物池和用于捕获孕妇血浆中游离DNA地贫基因的探针。
作为本发明的第三个方面，本发明还提供上述的试剂盒在制备诊断胎儿地中海贫血疾病的试剂盒中的应用。
综上所述，本发明的有益效果为：
1、本发明通过目标区域捕获结合高通量测序技术，实现对地贫无创产前检测；
2、本发明的方法所需的样本量少，仅需要5mL的外周血即可实现对样本的基因型准确检测；
3、本发明的方法具有无创性，通过孕妇外周血即可实现对胎儿基因型的准确检测；
4、本发明的方法在对父源突变检测的基础上，可以实现对胎儿母源基因突变的检测。
附图说明
图1为家系1通过生物信息学分析得到的RHDO-SPRT图；
图2为家系2通过生物信息学分析得到的RHDO-SPRT图；
图3为家系3通过生物信息学分析得到的RHDO-SPRT图；
图4为家系4通过生物信息学分析得到的RHDO-SPRT图；
图5为家系5通过生物信息学分析得到的RHDO-SPRT图。
具体实施方式
本发明提供一种无创产前胎儿α/β型地贫基因突变检测方法，具体包括地贫引物与探针设计，单倍型构建，父源遗传情况分析和母源遗传情况分析，其能够准确、高效地检测胎儿地贫基因突变，其结果与羊水穿刺检测结果一致，但安全性、无创性和高效性方面明显优于羊水穿刺检测。
为更好的说明本发明的目的、技术方案和优点，下面将结合附图和具体实施例对本发明作进一步说明。
样本说明：以下实施例所用的5个家系均收集自中山大学第三附属医院，分别编号为家系1-5，基因型如表1所示，其中包括了α-地贫的常见缺失型别--SEA/αα杂合和-α3.7/αα杂合及--SEA/-α3.7复合杂合和常见的突变型别αQSα/αα和αCSα/αα，还包括了β-地贫的常见突变型别β0/βN(IVS-2-654)、β0/βN((CD41-42(_CTTT))和β+/βN((-28(A>G))，还包括了野生型αα/αα&βN/βN。
表1 5个地贫家系基因型
实施例1单倍型(即单倍体基因型)的构建
(一)地贫引物与探针设计
针对α-地贫的设计区域(chr16:60022-2233661)：下游(10Kb)+HBA基因+上游(10Kb)，以及针对β-地贫血设计区域(chr11:3236690-7177304)：下游(10Kb)+HBB基因+上游(10Kp)，再利用dbSNP数据库筛选出SNP位点195个α-地贫SNP位点和275个β-地贫SNP位点。其中SNP位点筛选条件：SNP位点在千人基因组中国人南方人和北方人数据库中同时MAF>＝20％；SNP位点上下100bp序列是特异区域，在基因组上没有同源性；并且45％<GC<70％SNP位点没有连续3个相同碱基的情况。根据上述筛选的SNP位点，通过Ampliseq Designer网站(https://www.ampliseq.com/login/login.action)设计扩增引物对，通过安捷伦SureDesign网站(https://earray.chem.agilent.com/suredesign
/index.htm？sessiontimeout＝true)设计捕获探针并从该网站订购，设计的地贫引物和探针只需覆盖到筛选的SNP位点即可，引物序列和探针并非唯一。在表2中，列举了20对引物序列(10对HBA引物和10对HBB引物)和其对应的SNP位点作为示例。SNP位点详细信息见表15。
表2筛选的SNP位点和设计的对应的引物序列的示例
(二)检测步骤：
样本1-1,1-2,1-3；2-1,2-2，2-3；3-1,3-2,3-3；4-1,4-2,4-3；5-1,5-2,5-3用来单倍型的构建。具体实验步骤如下：
(1)DNA提取
按照DNeasy Blood&Tissue Kit(50)(厂家：Qiagen，货号：69504)试剂盒说明书对200μL的外周血进行基因组DNA提取。DNA溶液可置于(-20+5)℃保存。
(2)基因组文库构建
取25ng gDNA，按照Ion AmpliseqTM library Kit2.0-96LV(厂家：LifeTechnologies，货号:4480441)试剂盒说明书进行文库构建。文库溶液可置于(-20+5)℃保存。
(3)将文库稀释至50PM，将3个文库等量混合，从3个文库混合液中取2μL按照IonPGMTM HiQTM OT2Reagents 200(厂家：Life Technologies，货号：A26428)和Ion PGMTMHiQTMOT2Solutions 200(厂家：Life Technologies，货号：A26429)的试剂盒说明书进行模板制备和模板富集。
(4)按照Ion PGMTM HiQTM Sequencing 200Reagents(厂家：Life Technologies，货号：A26431)、Ion PGMTM HiQTM Sequencing 200Solutions(厂家：Life Technologies，货号：A26430)和Ion PGMTM Sequencing Nucleotides(厂家：Life Technologies，货号：A26432)的试剂盒说明书进行上机测序。
(三)单倍型的构建
测序得到平均5M reads的数据量，将人类基因组hg19作为参考序列，通过扩增基因组DNA地贫基因的引物池位点信息进行序列比对，分析目标区域的平均测序深度值，测序深度值(Sequencing Depth)≥500，进行下一步分析。利用GATK(3.4版本)将父-母-已育小孩同一染色体上均有信息的位点筛选出来，再以测序深度大于30×，测序碱基质量值大于Q20为筛选条件，筛选多态位点。碱基质量值，简称Q值，是碱基读取错误率p的取整映射结果。目前，Q值与p值的对应关系有两种计算方式：标准Sanger变体计算方式以及早期Solexa管道(如Illumina Genome Analyzer软件)计算方式。以标准Sanger变体方式为例，Q值被称为Phred quality score，它与p值之间的对应关系用以下公式计算得到：测序错误率与Q值的简明对应方式如下：Q13的碱基错误率为5％，Q20的碱基错误率为1％，Q30的碱基错误率为0.1％，Q40的碱基错误率为0.01％。再根据孟德尔遗传定律，把筛选得到的SNP基因型信息构建染色体单倍体型，并结合家系样本表型和已知突变型信息将致病突变关联到单倍体型上。
实施例2父源遗传情况分析
(一)检测步骤：
样本1-4，2-4，3-4，4-4和5-4用来父源和母源的遗传分析。具体实验步骤如下：
(1)DNA提取
按照QIAamp Circulating Nucleic Acid Kit(厂家：Qiagen，货号：55114)试剂盒说明书对2ml血浆进行游离DNA提取。DNA溶液可置于(-20+5)℃保存。
(2)文库构建
按照Ion Plus Fragment Library Kit(厂家：Life Technologies，货号：4471252)和sureselect Target Enrichment Kit(厂家：安捷伦，货号：51906514)试剂盒说明书对DNA进行文库构建。文库溶液可置于(-20+5)℃保存。
(3)将文库稀释至50PM，从1个文库稀释液中取2μL按照Ion PGMTM HiQTMOT2Reagents 200(厂家：Life Technologies，货号：A26428)和Ion PGMTM HiQTMOT2Solutions 200(厂家：Life Technologies，货号：A26429)的试剂盒说明书进行模板制备和模板富集。
(4)按照Ion PGMTM HiQTM Sequencing 200Reagents(厂家：Life Technologies，货号：A26431)、Ion PGMTM HiQTM Sequencing 200Solutions(厂家：Life Technologies，货号：A26430)和Ion PGMTMSequencing Nucleotides(厂家：Life Technologies，货号：A26432)的试剂盒说明书进行上机测序。
(二)数据分析
测序得到平均10M reads的数据量，将人类基因组hg19作为参考序列，通过捕获游离DNA地贫基因的探针位点信息进行序列比对，分析目标区域的平均测序深度值，测序深度值(Sequencing Depth)≥500，进行下一步分析。利用GATK(3.4版本)软件将筛选出孕妇血浆的变异位点。其中，挑选基因型为父亲杂合(MN，M和N代表碱基A,T,C,G)母亲纯合(MM)的位点和父亲和母亲均为纯合且基因型不同的位点(父亲MM，母亲NN)，在上述两类位点中，再筛选出孕妇血浆基因型为MN的位点，以此确定是否遗传了父源致病单倍体型。
实施例3母源遗传情况分析
在父源遗传情况分析基础上，进行母源遗传情况分析。挑选出基因型为父亲纯合(MM)母亲杂合的位点(MN)，根据RHDO-SPRT确定是否遗传了母源的致病单倍体型，其中，SPRT曲线根据以下公式计算：
其中，
Upper boundary与Lower bundary分别指的是上界和下界。发明人在此检测方法中把单倍体型1定义为上界，把单倍体型2定义为下界。用SPRT的算法，计算筛选出的位点对应深度的总和变异频率的总和。如果筛选出的位点都在上下界之间，那么，则不能分辨出胎儿遗传了母亲的哪一条单倍体型；如果累积的位点中有一个或一个以上超过了上界，则可以判断胎儿遗传了母亲的单倍体型1；如果累积的位点中有一个或一个以上超过了下界，则可以判断胎儿遗传了母亲的单倍体型2。
其中，N是用于此分类的所有位点的统计数；q0是胎儿遗传了母亲单倍体2的情况下，统计出reads来自于单倍体型1的比例；q1是胎儿遗传了母亲单倍体1的情况下，统计出reads来自于单倍体型1的比例。
实施例4 5个地贫家系检测结果
(1)地贫检测结果
通过实施例2的方法对家系1-5进行检测，家系1-5的单倍型构建结果如表3-7所示。根据孟德尔遗传定律，并结合家系样本表型和已知突变型信息将致病突变关联到单倍体型上(例如：父亲基因型为MN，母亲基因型为MN，已育小孩基因型为NN(M为致病位点，N为正常位点。M和N代表碱基A,T,C,G)，已育小孩分别遗传了父亲的单倍型1和母亲的单倍型1，则父亲致病位点M在单倍型2上，母亲致病位点M在单倍型2上)。家系1-5的父亲和母亲致病位点关联到单倍型上的结果如表8所示。
表3家系1的单倍型构建结果
表4家系2的单倍型构建结果
表5家系3的单倍型构建结果
表6家系4的单倍型构建结果
表7家系5的单倍型构建结果
表8家系1-5的父亲和母亲致病位点与单倍型的关联结果
通过实施例2和3的方法对家系1-5进行检测，表9-13分别为家系1-5的胎儿父源遗传情况列表。从表9-13可以看出家系1-5的胎儿遗传了父亲单倍体型1；图1-图5分别为家系1-5的RHDO-SPRT图。从图1-图5可以看出家系1-5的胎儿遗传了母亲的单倍体型1。综上所述，家系1-5胎儿的地贫检测结果如表14所示。
表9家系1的胎儿父源遗传情况
表10家系2的胎儿父源遗传情况
表11家系3的胎儿父源遗传情况
表12家系4的胎儿父源遗传情况
表13家系5的胎儿父源遗传情况
表14家系1-5胎儿无创地贫检测结果
结果显示，本发明的检测方法对α&β-地贫的临床样本的检测结果(表14)与原始羊水检测型别一致，表明本发明的方法能够准确、高效地实现对α&β-地贫无创产前检测检测，其结果与羊水穿刺检测分型一致，但安全性、无创性和高效性方面明显优于羊水穿刺检测。
表15 SNP位点信息(其中起始位置相对终止位置减1)
最后所应当说明的是，以上实施例仅用以说明本发明的技术方案而非对本发明保护范围的限制，尽管参照较佳实施例对本发明作了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的实质和范围。
SEQUENCE LISTING
<110> 广州市达瑞生物技术股份有限公司
<120> 检测胎儿地中海贫血致病基因的方法和试剂盒
<130> 2018
<160> 40
<170> PatentIn version 3.5
<210> 1
<211> 22
<212> DNA
<213> 人工序列
<400> 1
ccaaggtgcc aacctcagaa tt 22
<210> 2
<211> 22
<212> DNA
<213> 人工序列
<400> 2
gatgtccacg gagtcgaagt ac 22
<210> 3
<211> 19
<212> DNA
<213> 人工序列
<400> 3
cctgcccagt tatgggtgt 19
<210> 4
<211> 23
<212> DNA
<213> 人工序列
<400> 4
agacaccaag caccttgtga ata 23
<210> 5
<211> 21
<212> DNA
<213> 人工序列
<400> 5
caagacgatg agggtcctgt t 21
<210> 6
<211> 21
<212> DNA
<213> 人工序列
<400> 6
ctgtgcccac ataccagtgt a 21
<210> 7
<211> 22
<212> DNA
<213> 人工序列
<400> 7
ccaggctgca gattcctgaa at 22
<210> 8
<211> 22
<212> DNA
<213> 人工序列
<400> 8
agggctggaa agctcattct tc 22
<210> 9
<211> 23
<212> DNA
<213> 人工序列
<400> 9
ctcagctgat tcctcatcct tca 23
<210> 10
<211> 22
<212> DNA
<213> 人工序列
<400> 10
tgtgcctcct tgcgttcatt ta 22
<210> 11
<211> 22
<212> DNA
<213> 人工序列
<400> 11
cgtcacccta tttgcatgct ag 22
<210> 12
<211> 22
<212> DNA
<213> 人工序列
<400> 12
ggatcacatg gctgacaaaa gc 22
<210> 13
<211> 25
<212> DNA
<213> 人工序列
<400> 13
cttatcgttg atcctgtttg gtcct 25
<210> 14
<211> 27
<212> DNA
<213> 人工序列
<400> 14
cccaatagat gggttcaaaa gtaaacc 27
<210> 15
<211> 25
<212> DNA
<213> 人工序列
<400> 15
tctcagactg ctctactctt ctgtg 25
<210> 16
<211> 22
<212> DNA
<213> 人工序列
<400> 16
gtgatgaaga gccaaggtca gt 22
<210> 17
<211> 21
<212> DNA
<213> 人工序列
<400> 17
gcacagctgg tcaacagcta t 21
<210> 18
<211> 22
<212> DNA
<213> 人工序列
<400> 18
cctgtaagcc aagaagcctt ct 22
<210> 19
<211> 21
<212> DNA
<213> 人工序列
<400> 19
ctccactctc ctctccctca a 21
<210> 20
<211> 20
<212> DNA
<213> 人工序列
<400> 20
cctgagcaca aacccgtaag 20
<210> 21
<211> 23
<212> DNA
<213> 人工序列
<400> 21
ccaggtgatt catttggcct tga 23
<210> 22
<211> 25
<212> DNA
<213> 人工序列
<400> 22
ctacagccac atggaaggac aaggt 25
<210> 23
<211> 24
<212> DNA
<213> 人工序列
<400> 23
tgttcttagt cctgggagat gcat 24
<210> 24
<211> 24
<212> DNA
<213> 人工序列
<400> 24
cacacacaca gtctgtctca tctc 24
<210> 25
<211> 24
<212> DNA
<213> 人工序列
<400> 25
gctgtgtggc cttgagagaa taac 24
<210> 26
<211> 24
<212> DNA
<213> 人工序列
<400> 26
agagcaaaat gcggaaaaca atcc 24
<210> 27
<211> 24
<212> DNA
<213> 人工序列
<400> 27
gtgcctggcc aaagagaaat ttat 24
<210> 28
<211> 30
<212> DNA
<213> 人工序列
<400> 28
gatggcacct aaagtttttc ttaaaactcc 30
<210> 29
<211> 22
<212> DNA
<213> 人工序列
<400> 29
agcaaatgcc ctatcatgcc aa 22
<210> 30
<211> 29
<212> DNA
<213> 人工序列
<400> 30
gtgagggagt ataagaaata agcatgaca 29
<210> 31
<211> 29
<212> DNA
<213> 人工序列
<400> 31
gccagtaagt gttactattc tctcattcc 29
<210> 32
<211> 26
<212> DNA
<213> 人工序列
<400> 32
tagttatatt gcccacaatt gtggct 26
<210> 33
<211> 24
<212> DNA
<213> 人工序列
<400> 33
gacttttcct ccacttaagc aggt 24
<210> 34
<211> 30
<212> DNA
<213> 人工序列
<400> 34
tggtttctac acagtcatta aagaagcatt 30
<210> 35
<211> 22
<212> DNA
<213> 人工序列
<400> 35
gcttgaatgt gcagagaccc aa 22
<210> 36
<211> 28
<212> DNA
<213> 人工序列
<400> 36
ccttcaagcc acagtcatat tcttagat 28
<210> 37
<211> 25
<212> DNA
<213> 人工序列
<400> 37
agataaatga atgcctgcca actgt 25
<210> 38
<211> 25
<212> DNA
<213> 人工序列
<400> 38
cggtttgttt cgttttcagt aggga 25
<210> 39
<211> 24
<212> DNA
<213> 人工序列
<400> 39
ccattgaacc agaactcaat ccct 24
<210> 40
<211> 29
<212> DNA
<213> 人工序列
<400> 40
ccaccataga gaaagagaaa aaggaatca 29
Claims (2)
1.一种检测胎儿地中海贫血致病基因的试剂盒，其特征在于，所述试剂盒包括扩增基因组地贫基因的引物池，所述引物池由SEQ ID NO:1～SEQ ID NO:40的序列所示。
2.如权利要求1所述的试剂盒在制备诊断胎儿地中海贫血疾病的试剂盒中的应用。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810473083.4A CN108642160B (zh) | 2018-05-16 | 2018-05-16 | 检测胎儿地中海贫血致病基因的方法和试剂盒 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810473083.4A CN108642160B (zh) | 2018-05-16 | 2018-05-16 | 检测胎儿地中海贫血致病基因的方法和试剂盒 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108642160A CN108642160A (zh) | 2018-10-12 |
| CN108642160B true CN108642160B (zh) | 2022-03-11 |
Family
ID=63756300
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810473083.4A Active CN108642160B (zh) | 2018-05-16 | 2018-05-16 | 检测胎儿地中海贫血致病基因的方法和试剂盒 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108642160B (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112996926A (zh) * | 2018-12-07 | 2021-06-18 | 深圳华大生命科学研究院 | 一种靶基因文库的构建方法、检测装置及其应用 |
| CN111321210B (zh) * | 2018-12-13 | 2021-09-14 | 贵阳优乐复生医学检验所有限公司 | 一种无创产前检测胎儿是否患有遗传疾病的方法 |
| CN112080558B (zh) * | 2019-06-13 | 2024-03-12 | 杭州贝瑞和康基因诊断技术有限公司 | 同时检测hba1/2和hbb基因突变的试剂盒和方法 |
| CN110373458B (zh) * | 2019-06-27 | 2020-05-19 | 东莞博奥木华基因科技有限公司 | 一种地中海贫血检测的试剂盒及分析系统 |
| CN111192633A (zh) * | 2020-01-07 | 2020-05-22 | 深圳市早知道科技有限公司 | 一种预测地中海贫血病表型的方法及终端设备 |
| CN111471752A (zh) * | 2020-03-27 | 2020-07-31 | 广州辉园苑医药科技有限公司 | 一种用于检测hbb基因位点突变的引物及试剂盒 |
| CN111500574B (zh) * | 2020-05-07 | 2023-11-17 | 和卓生物科技（上海）有限公司 | 一种用于检测遗传性耳聋的探针组合及其应用 |
| CN111951890B (zh) * | 2020-08-13 | 2022-03-22 | 北京博昊云天科技有限公司 | 染色体和单基因病同步产前筛查的设备、试剂盒和分析系统 |
| CN112126677B (zh) * | 2020-11-25 | 2021-02-19 | 北京迈基诺基因科技股份有限公司 | 耳聋单倍型基因突变无创检测方法 |
| CN112575077A (zh) * | 2020-12-23 | 2021-03-30 | 东莞市妇幼保健院 | 一种胎儿显性遗传病新发突变的无创基因检测方法及应用 |
| CN112831555B (zh) * | 2021-02-01 | 2022-12-20 | 人和未来生物科技（长沙）有限公司 | 一种检测地中海贫血基因的参考品及其制备方法与应用 |
| CN112708674B (zh) * | 2021-03-29 | 2021-06-29 | 北京贝瑞和康生物技术有限公司 | 同时检测hba1/2和hbb基因位点多种突变的方法和试剂盒 |
| CN115019888A (zh) * | 2022-07-14 | 2022-09-06 | 苏州贝康医疗器械有限公司 | 一种基于外周血游离dna高通量测序的组织特异性基因标志物的筛选系统及其应用 |
| CN116959574A (zh) * | 2023-07-28 | 2023-10-27 | 苏州贝康医疗器械有限公司 | 一种针对于无创产前诊断领域的新型捕获探针及其设计方法和应用 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106755486A (zh) * | 2017-01-19 | 2017-05-31 | 人和未来生物科技（长沙）有限公司 | 无创产前胎儿β 型地贫基因突变检测文库构建方法、检测方法和试剂盒 |
| CN107099609A (zh) * | 2017-06-19 | 2017-08-29 | 东莞市妇幼保健院 | 一种微滴式数字pcr检测低比例痕量目标基因的遗传分析方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104232778B (zh) * | 2014-09-19 | 2016-08-17 | 天津华大基因科技有限公司 | 同时确定胎儿单体型及染色体非整倍性的方法及装置 |
| EP3294906B1 (en) * | 2015-05-11 | 2024-07-10 | Natera, Inc. | Methods for determining ploidy |
| CN107217095B (zh) * | 2017-06-15 | 2021-06-04 | 广东腾飞基因科技股份有限公司 | 人类亲权鉴定用多重pcr引物组及检测方法 |
-
2018
- 2018-05-16 CN CN201810473083.4A patent/CN108642160B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106755486A (zh) * | 2017-01-19 | 2017-05-31 | 人和未来生物科技（长沙）有限公司 | 无创产前胎儿β 型地贫基因突变检测文库构建方法、检测方法和试剂盒 |
| CN107099609A (zh) * | 2017-06-19 | 2017-08-29 | 东莞市妇幼保健院 | 一种微滴式数字pcr检测低比例痕量目标基因的遗传分析方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108642160A (zh) | 2018-10-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108642160B (zh) | 检测胎儿地中海贫血致病基因的方法和试剂盒 | |
| CN104232777B (zh) | 同时确定胎儿核酸含量和染色体非整倍性的方法及装置 | |
| CN103874767B (zh) | 对核酸样本中预定区域进行基因分型的方法和系统 | |
| KR101966262B1 (ko) | 대규모 병렬 게놈 서열분석을 이용한 태아 염색체 이수성의 진단 방법 | |
| CN104232778B (zh) | 同时确定胎儿单体型及染色体非整倍性的方法及装置 | |
| TWI727938B (zh) | 血漿粒線體ｄｎａ分析之應用 | |
| WO2024027569A1 (zh) | 一种不依赖于先证者的单体型构建方法 | |
| CN105648045B (zh) | 确定胎儿目标区域单体型的方法和装置 | |
| EP3564391A1 (en) | Method, device and kit for detecting fetal genetic mutation | |
| CN104531883A (zh) | Pkd1基因突变的检测试剂盒及检测方法 | |
| WO2015042980A1 (zh) | 确定染色体预定区域中snp信息的方法、系统和计算机可读介质 | |
| HK1231503A1 (zh) | 基因突变体及其应用 | |
| CN106906220B (zh) | 一种突变的col4a5基因及其应用 | |
| CN111508561B (zh) | 同源序列和同源序列中串联重复序列的检测方法、计算机可读介质和应用 | |
| CN111321210A (zh) | 一种无创产前检测胎儿是否患有遗传疾病的方法 | |
| CN104232649B (zh) | 基因突变体及其应用 | |
| CN108531580B (zh) | C5orf42基因突变体及其应用 | |
| CN111334513A (zh) | 一种非综合征型唇裂相关的低频/罕见突变及其检测方法 | |
| CN105838720A (zh) | Ptprq基因突变体及其应用 | |
| CN109097465B (zh) | Clip3基因的snp位点的应用 | |
| CN114292908A (zh) | 检测地中海贫血基因罕见突变的试剂、方法和试剂盒 | |
| CN106636351B (zh) | 一种与乳腺癌相关的snp标记及其应用 | |
| CN119193822B (zh) | 一种评价遗传多发性骨软骨瘤发病风险相关的snp标志物及其应用 | |
| CN112442503A (zh) | Kcnq1基因突变体及其应用 | |
| US20180179595A1 (en) | Fetal haplotype identification |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |