CN103602735B - 利用高通量测序精确测定线粒体dna高频和低频突变的方法 - Google Patents
利用高通量测序精确测定线粒体dna高频和低频突变的方法 Download PDFInfo
- Publication number
- CN103602735B CN103602735B CN201310544351.4A CN201310544351A CN103602735B CN 103602735 B CN103602735 B CN 103602735B CN 201310544351 A CN201310544351 A CN 201310544351A CN 103602735 B CN103602735 B CN 103602735B
- Authority
- CN
- China
- Prior art keywords
- mitochondrial dna
- dna
- frequency
- flux sequence
- library
- 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
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
- C12Q2531/00—Reactions of nucleic acids characterised by
- C12Q2531/10—Reactions of nucleic acids characterised by the purpose being amplify/increase the copy number of target nucleic acid
- C12Q2531/125—Rolling circle
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Organic Chemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Biotechnology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明属于基因组学技术领域，具体为一种利用高通量测序精确测定线粒体DNA高频和低频突变的方法。本发明通过设计一组线粒体DNA特异的引物，采用滚环复制从总DNA中富集线粒体DNA，并通过特定的限制性内切酶进一步富集线粒体DNA，然后构建测序文库并进行高通量测序，并开发特定的数据分析流程，从而达到测定线粒体基因组中高频率及低频率变异的目的。其最低可检测频率低至0.3%的突变位点。
Description
技术领域
本发明属于基因组学技术领域，具体涉及一种利用高通量测序（下一代测序、第二代测序、深度测序）技术测定真核生物线粒体DNA中高频率及低频率突变的方法。
背景技术
线粒体是细胞中提供能量的重要细胞器，在生物体能量转换和新陈代谢中处于核心地位。人体细胞中线粒体DNA(脱氧核糖核酸)的突变或缺失会引起氧化磷酸化和能量供应的异常，并进而引起近150种人类疾病的发生，比如心血管疾病、糖尿病、耳聋、肿瘤、神经退行性疾病及衰老等。
线粒体DNA是核外遗传物质，可以独立编码参与线粒体功能的一些重要蛋白。人类线粒体DNA为大小16571碱基的闭合双链分子，编码22种tRNA(转运核糖核酸)、2种rRNA(核糖体RNA)和氧化磷酸化复合物中的13个亚基。线粒体DNA的变异包括大片段的***或缺失，小片段的***或缺失，点突变。这些变异既可以出现在编码区，也可以出现在非编码区。体内活性氧(ROS,reactiveoxygenspecies)的异常增多是引起线粒体DNA碱基变异或点突变的主要原因。线粒体DNA上的某些区域由于受到活性氧引起的自由基的持续攻击，产生点突变的频率明显高于其它区域，则称之为突变热点。比如亮氨酸转运RNA基因的点突变就是线粒体DNA的一个病原学突变热点。研究表明，一些编码基因的点突变与某些特定的疾病有明显的相关性。比如，在色素性视网膜炎、一些神经疾病和Leigh综合征中，均检测到ATP合成酶中的亚基6上的点突变。并且，这些疾病的严重程度与该基因的点突变频率密切相关。疾病发生早期线粒体DNA的低频突变已经发生，但目前尚无有效的高通量检测方法。
线粒体基因组与核基因组在遗传上的重要不同是：在一个细胞中，核基因组只有两个拷贝，一个来自父亲，一个来自母亲；而线粒体基因组几乎完全来自母亲，并且有上千个拷贝（通常一个细胞有1000-2000个线粒体，一个线粒体有2-10个线粒体DNA拷贝）。这种拷贝数的不同也导致了在DNA突变检测上的策略不同。对于基因组上某一位置的点突变只可能出现三种情况：均无突变，有一个拷贝发生突变，或者两个拷贝均发生突变。而对于线粒体DNA来说，情况就要复杂的多，从一个线粒体DNA发生突变，到上千个线粒体DNA均发生突变，其中的突变频率可以是0.1%~100%中的任何一个比例。而要精确检测到这种突变频率的变化，特别是低频率(0.1%~10%)的点突变，用传统的桑格测序(Sangersequencing)是无法进行精确测定的。
随着高通量测序技术（第二代测序、下一代测序、深度测序）的飞速发展，单个碱基的测序费用迅速下降，测序的覆盖度也得以不断加深，使得之前所不能实现的目标有了新的可能。从理论上讲，要检测较低频率的突变，比如说0.4%，则至少需要250倍的覆盖度，加上每个点突变至少被检测到两次的标准以提高准确性，则至少需要500倍的覆盖度。深度测序的技术正好能够满足这样的覆盖度要求，因此使得检测线粒体DNA的低频突变成为可能。目前已有几个方法测定线粒体的突变，比如从细胞中分离线粒体，再提取线粒体DNA并构建文库进行深度测序。这种方法由于涉及到线粒体细胞器的分离，操作比较复杂繁琐，尤其是在样品多的时候工作量很大。另一种深度测序方法基于聚合酶链式反应(PCR)扩增，针对每个样品需要分10个以上管子进行PCR扩增，纯化后连接再打断，然后再构建高通量测序文库，并进行深度测序。这种方法虽然能从总DNA中扩增线粒体DNA片段，但每个样品需要分别做10个以上的不同引物对的扩增，还需要分别纯化，耗时耗力，并且不同扩增片段之间的效率也不尽相同，因此在线粒体基因组上的覆盖度均一性也较差。
为了克服以上方法的缺点，本发明提出了一种更简单、更高效的线粒体DNA深度测序方法，并取名为mitoRCA-seq。这种方法通过设计一组线粒体DNA特异的引物，采用滚环复制（RollingCircleAmplification,RCA）从总DNA中富集线粒体DNA，并通过特定的限制性内切酶进一步富集线粒体DNA，然后构建测序文库并进行高通量测序，并开发特定的数据分析流程，从而达到测定线粒体基因组中高频率及低频率点突变以及其它结构变异的目的。
发明内容
本发明的目的在于提供一种利用高通量测序技术测定真核生物线粒体DNA中高频率及低频率突变的方法。
本发明提供的测定真核生物线粒体基因组中高频率及低频率变异的方法，采用高通量测序技术。本方法的关键之处在于设计一组（38条，具体序列见序列表）线粒体DNA特异的引物，采用滚环复制从总DNA中富集线粒体DNA，并通过特定的限制性内切酶进一步富集线粒体DNA，然后构建测序文库并进行高通量测序，并开发特定的数据分析流程，从而达到测定线粒体基因组中高频率及低频率变异的目的（见图1）。由于每个样品在接头连接步骤中加入了条形码，因此可以将多个样品混合在一个通道中进行测序，随后在数据分析时可根据条形码将不同样品分开，使得单个样品的费用降低。
本发明方法为一种线粒体DNA高通量测序文库的构建新方法(简称：mitoRCA-seq)。起始材料为总DNA，终产物为扩增后的双链DNA测序文库。该文库经高通量测序（比如Illumina的GAIIx测序平台或HiSeq2000/HiSeq2500测序平台）后获得线粒体DNA的上百万个短序列基础数据，再经过生物信息学分析获得线粒体DNA高频率及低频率点突变及其它基因组变异等详细信息。采用本发明方法最低可检测频率低至0.3%的突变位点。
本发明方法的具体步骤如下：
（1）设计线粒体DNA特异性引物SEQIDNO.1-SEQIDNO.38；利用高保真的滚环复制从总DNA中直接扩增线粒体DNA；
（2）通过特定的限制性内切酶，进一步将扩增的线粒体DNA酶切成2条特异条带，以进一步富集线粒体DNA，这2条条带的长度加在一起为线粒体DNA的全长；
（3）回收此2条条带并片段化；连接接头后进行低循环PCR或无PCR文库构建；构建的文库可以单端测序，也可以双端测序；
（4）该文库经高通量测序获得高质量数据，再经过生物信息学分析获得真核生物线粒体DNA所有位点的点突变频率和其它类型的结构变异。
本发明方法的核心在于通过两步策略富集线粒体DNA，其一为采用多对线粒体DNA特异性引物，其二为通过特定的限制性内切酶将滚环复制扩增的产物切出2条特异条带。本发明在用这种方法对来自果蝇、小鼠、大鼠以及人类的总DNA进行了测试并获得成功。因此这种方法不管从原理上还是实际操作上均可适用于所有含有线粒体的真核生物，这些真核生物包括但不限于人类、猩猩、猕猴、大鼠、小鼠、果蝇、斑马鱼、线虫、拟南芥、水稻、棉花、油菜、酵母及疟原虫。此外，这一技术的核心在于线粒体DNA的富集步骤，而测序步骤既使用于Illumina测序平台（包括GAIIx、HiSeq2000、HiSeq2500、MiSeq），也适用于454测序平台和ABI的ABISOLiD4、ABI5500系列测序平台以及IonTorrent和IonProton测序平台。
本发明对这一方法的可重复性进行了验证（图2），结果表明三次不同技术重复均表现了很好的重复性。第一次和第二次的低频点突变相关系数达到了0.998(RCA1-1andRCA1-2)，而第一次和第三次相关系数也达到了0.996(RCA1-1andRCA2)，第二次和第三次相关系数为0.997(RCA1-2andRCA2)。我们随后利用Polg基因突变小鼠和野生型小鼠对这一方法(mitoRCA-seq)进行了进一步验证。Polg这一基因编码线粒体DNA聚合酶，Polg基因的突变导致了聚合酶复制修复功能的丧失，因而比野生型有更多的低频点突变。我们利用mitoRCA-seq测定了野生型和Polg突变小鼠的线粒体DNA所有位点的高频和低频突变，获得了很好的覆盖度（图3），并且发现Polg突变小鼠中确实含有更多的位点发生了低频点突变（图4）。这一结果也被另一个生物学重复实验所证实。以上的实验证实了mitoRCA-seq确实是一种有效并且可靠的方法用于检测线粒体DNA的突变，不管是高频率的，还是低频率的。
附图说明
图1是mitoRCA-seq方法示意图。橙色圆环表示线粒体DNA分子。
图2是mitoRCA-seq三次用于同一样品的点突变频率的相关性结果。用Polg突变小鼠肝脏样品对mitoRCA-seq的技术可重复性进行了验证，结果显示第一次和第二次的低频点突变相关系数(R)达到了0.998(RCA1-1andRCA1-2)，而第一次和第三次相关系数也达到了0.996(RCA1-1andRCA2)，第二次和第三次相关系数为0.997(RCA1-2andRCA2)。这一结果表明mitoRCA-seq具有很好的技术可重复性。
图3是野生型小鼠和Polg突变小鼠肝脏mitoRCA-seq所测定的测序覆盖度和点突变的位点和频率。从测序覆盖度上可以看出这一方法具有较好的均一性。而不同位点的点突变检测也可以清晰地看出。覆盖度图的纵坐标表示log2转化的测序深度，而点突变频率图中的纵坐标表示0%-100%的频率。野生型小鼠肝脏样品表示野生型小鼠生物学重复1的肝脏，Polg突变小鼠肝脏样品表示Polg突变小鼠生物学重复1的肝脏。小鼠线粒体基因组注释Ensembl在图的最下方。
图4是野生型小鼠和Polg突变小鼠大脑和肝脏突变位点数的统计图。野生型小鼠重复1和2指生物学重复。Polg突变小鼠重复1和2亦指生物学重复。左侧为大脑数据，右侧为肝脏数据。每条曲线的具体样品信息按照指示箭头所指。这一结果显示Polg突变小鼠中确实含有更多的位点发生了低频点突变。
具体实施方式
所有实施例中的有关DNA和RNA基本操作均参考《分子克隆：实验室操作指南》（金冬雁等译，科学出版社，北京（1993））和《精编分子生物学指南》（颜子颖等译，科学出版社，北京（1998））。分子操作中所用的酶和特殊试剂在之后括号中均注明了相应的公司。
实施例1：总DNA的提取
总DNA(包括人、猩猩、猕猴、大鼠、小鼠、果蝇、斑马鱼、线虫、拟南芥、水稻、棉花、油菜、酵母及疟原虫)的提取采用Qiagen公司的QIAampDNABloodMinikit试剂盒。不同的细胞或者组织按照Qiagen公司建议的流程来提取。通常来说，组织取样为20毫克，加入200μlAL缓冲液。匀浆采用钢珠，在2ml的圆底离心管中进行，采用的仪器为Qiagen公司的TissueLyser匀浆机，参数为每秒50下，共匀浆2分钟。柱子纯化后，洗脱到200μl纯水中。采用Nanodropspectrophotometer测定DNA浓度。100ng的总DNA用于构建mitoRCA-seq文库。
实施例2：采用滚环复制进行线粒体DNA扩增
100ng总DNA用于滚环扩增。将100ng总DNA加入50ul反应体系中，该体系中还包括：1xPhi29DNA聚合酶缓冲液(NEB),0.2μg/mlBSA,1mMdNTP(NEB)and25μM线粒体DNA特异性引物（小鼠引物详见序列表）。混匀这一体系并在95°C加热3分钟，然后冷却至室温，随后加入1μlPhi29DNA聚合酶(10unit/μl,NEB)。混匀后在37°C反应16小时。随后加热至65°C并保持10分钟以灭活该聚合酶。Qiagen公司的REPLI-gMitochondrialDNAkit试剂盒也可以用于滚环复制。
实施例3：限制性酶切
不同的物种需要选择不同的限制性内切酶。原则是选择1-2个限制性内切酶能够使得对应的线粒体DNA产生2条可以在琼脂糖凝胶电泳上分离开的条带。我们选择EcoRV(NEB)来产生9.5kb和6.8kb的2条条带。对于果蝇来说，NdeI(NEB)和EcoRV(NEB)可用于产生10kb和9.8kb的条带。对于人的样品的话，SacI(NEB)可用于产生6.9kb和9.6kb的2个条带。具体的限制性酶切体系按照每个酶的建议用量来实施。
实施例4：双链DNA的片段化、补平和3端加腺苷酸
我们从0.5%的琼脂糖凝胶电泳上切下实施例3中的2条条带，并采用ZYMO公司的Zymocleanlargefragmentgelpurificationkit试剂盒回收条带。片段化采用CovarisS2机器，根据最终的峰值大小300bp调参数。纯化并补平的DNA用Klenow(exo-)DNA聚合酶(Epicentre)和200μMdATP加上3端腺苷酸突出。该反应体系在37℃保持30分钟，然后用ZYMOResearch公司的clean&concentrator-5试剂盒纯化并洗脱至7ul纯水中。
实施例5：接头的连接反应
在实施例4中得到的7ulDNA中加入1ul10倍T4DNA连接酶缓冲液、3pmol的Illumina公司接头（或者SOLiD平台接头）和1ul高浓度T4DNA连接酶(NEB;2000units/μl)。混匀后在室温放置30分钟，随后用ZYMOResearch公司的clean&concentrator-5试剂盒纯化并洗脱至20ul纯水中。将20ul连接并纯化后的产物跑2%的琼脂糖凝胶电泳，切出250-450bp片段并用ZYMOResearch公司的DNARecovery试剂盒纯化并洗脱至20ul纯水中。
实施例6：低循环PCR扩增及无PCR的文库构建
连接好接头的双链DNA可在以下的50ulPCR体系中进行扩增：20ul的实施例5中获得的DNA模版、1倍浓度的Finnnzymes公司的HF缓冲液、1nmoldNTP、25pmol的Illumina正向PCR引物和25pmol的Illumina反向引物以及0.5μl的Phusion高保真DNA聚合酶(Finnzymes)。热循环的参数如下：98°Cfor30s;12cyclesof98°Cfor10s,67°Cfor30sand72°Cfor30s;72°Cfor10min;holdat10°C。取出20ulPCR产物跑2%琼脂糖凝胶电泳，切出300-500bp产物并用ZYMOResearch公司的DNARecovery试剂盒纯化并洗脱至20ul纯水中。用Qubit荧光仪(Invitrogen)测定文库的浓度。加入终浓度为0.1%的Tween-20用于测序前的样品保存。用Illumina公司的GAIIx或HiSeq2000/2500或MiSeq进行双端的测序，每端测定50个碱基。454测序文库、ABISOLiD4文库、ABI5500文库及ABIIonTorrent文库采用对应的测序仪测序。文库的构建除接头不同需用对应的PCR引物外，其余均完全一样。对于无PCR步骤的文库构建，可以按照NuGEN公司的EncoreRapidLibrarySystems试剂盒操作流程来进行。
实施例7：生物信息学分析
混合样品的原始数据通过构建文库时所用的条形码进行分开。对于来源不同物种的样品，我们用序列比对软件bwa(采用默认参数)将其映射到对应物种的最新线粒体DNA基因组版本上。分析所用到的工具包括samtools,bamtools以及bedtools，均为开放源软件。我们将所有具有唯一映射位置的短序列用于下游的点突变频率分析和其它类型突变的分析（比如小的***或缺失）。
SEQUENCELISTING
<110>复旦大学
<120>利用高通量测序精确测定线粒体DNA高频和低频突变的方法
<130>001
<160>38
<170>PatentInversion3.3
<210>1
<211>21
<212>DNA
<213>
<400>1
aacaaagcaaagcactgaaaa21
<210>2
<211>21
<212>DNA
<213>
<400>2
gcgtacttcattgctcaattc21
<210>3
<211>21
<212>DNA
<213>
<400>3
cttcagcaaaccctaaaaagg21
<210>4
<211>21
<212>DNA
<213>
<400>4
ttttatgttgagcttgaacgc21
<210>5
<211>21
<212>DNA
<213>
<400>5
attccaattctccaggcatac21
<210>6
<211>21
<212>DNA
<213>
<400>6
cctttcgtactgggagaaatc21
<210>7
<211>21
<212>DNA
<213>
<400>7
tttgatcaacggaccaagtta21
<210>8
<211>21
<212>DNA
<213>
<400>8
ataaagaataacgcgaatggg21
<210>9
<211>21
<212>DNA
<213>
<400>9
actattcggagctttacgagc21
<210>10
<211>21
<212>DNA
<213>
<400>10
ggccaggaggataattattga21
<210>11
<211>21
<212>DNA
<213>
<400>11
caactgaagcagcaacaaaat21
<210>12
<211>21
<212>DNA
<213>
<400>12
ggggtagggttattgtgctta21
<210>13
<211>21
<212>DNA
<213>
<400>13
taagcacaataaccctacccc21
<210>14
<211>21
<212>DNA
<213>
<400>14
ggtgggtagactgttcatcct21
<210>15
<211>21
<212>DNA
<213>
<400>15
aggctttggaaactgacttgt21
<210>16
<211>21
<212>DNA
<213>
<400>16
tgtcaagggatgagttggata21
<210>17
<211>21
<212>DNA
<213>
<400>17
tatccaactcatcccttgaca21
<210>18
<211>21
<212>DNA
<213>
<400>18
tagcagtcgtagttcaccagg21
<210>19
<211>21
<212>DNA
<213>
<400>19
catagggcaccaatgatactg21
<210>20
<211>21
<212>DNA
<213>
<400>20
cagctcatagtggaatggcta21
<210>21
<211>21
<212>DNA
<213>
<400>21
tcctattcccatcctcaaaac21
<210>22
<211>21
<212>DNA
<213>
<400>22
tgttggtacgaggctagaatg21
<210>23
<211>21
<212>DNA
<213>
<400>23
aaaccacataaatcaagccct21
<210>24
<211>21
<212>DNA
<213>
<400>24
tgtggatattaggtgagagcg21
<210>25
<211>21
<212>DNA
<213>
<400>25
tacaagctctgcacgtctacc21
<210>26
<211>21
<212>DNA
<213>
<400>26
atgaagcgtctaaggtgtgtg21
<210>27
<211>21
<212>DNA
<213>
<400>27
actgctaattgccctcatctt21
<210>28
<211>21
<212>DNA
<213>
<400>28
gtgtgagggttggaggttaat21
<210>29
<211>21
<212>DNA
<213>
<400>29
cttcaaatggtcttcccactt21
<210>30
<211>21
<212>DNA
<213>
<400>30
ttgatgtttgggtctgagtgt21
<210>31
<211>21
<212>DNA
<213>
<400>31
caacttttcattggctgagaa21
<210>32
<211>21
<212>DNA
<213>
<400>32
taattagtagggctcaggcgt21
<210>33
<211>21
<212>DNA
<213>
<400>33
actccaacatcatcaacctca21
<210>34
<211>21
<212>DNA
<213>
<400>34
tcaaggtggctttgtctactg21
<210>35
<211>21
<212>DNA
<213>
<400>35
cccatatattggaacaaccct21
<210>36
<211>21
<212>DNA
<213>
<400>36
agcttatatgcttggggaaaa21
<210>37
<211>21
<212>DNA
<213>
<400>37
aagaagaaggagctactcccc21
<210>38
<211>21
<212>DNA
<213>
<400>38
gagttttggttcacggaacat21
Claims (2)
1.一种利用高通量测序精确测定线粒体DNA高频和低频突变的方法，其特征在于具体步骤如下：
（1）设计小鼠线粒体DNA特异性引物SEQIDNO.1-SEQIDNO.38；利用高保真的滚环复制从小鼠的总DNA中直接扩增线粒体DNA；
（2）通过特定的限制性内切酶，进一步将扩增的线粒体DNA酶切成2条特异条带，以进一步富集线粒体DNA，这2条条带的长度加在一起为线粒体DNA的全长；
（3）回收此2条条带并片段化；连接接头后进行低循环PCR或无PCR文库构建；构建的文库为单端测序或双端测序文库；
（4）该文库经过高通量测序获得高质量数据，再经过生物信息学分析获得小鼠线粒体DNA所有位点的点突变频率和其它类型的结构变异；
其最低可检测频率低至0.3%的突变位点。
2.根据权利要求1所述的方法，其特征在于构建的文库进行高通量测序的平台为IlluminaGAIIx平台、IlluminaHiSeq2000、IlluminaHiSeq2500、IlluminaMiSeq平台。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310544351.4A CN103602735B (zh) | 2013-11-06 | 2013-11-06 | 利用高通量测序精确测定线粒体dna高频和低频突变的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310544351.4A CN103602735B (zh) | 2013-11-06 | 2013-11-06 | 利用高通量测序精确测定线粒体dna高频和低频突变的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103602735A CN103602735A (zh) | 2014-02-26 |
| CN103602735B true CN103602735B (zh) | 2016-03-30 |
Family
ID=50120998
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310544351.4A Active CN103602735B (zh) | 2013-11-06 | 2013-11-06 | 利用高通量测序精确测定线粒体dna高频和低频突变的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103602735B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105653899B (zh) * | 2014-09-30 | 2018-02-09 | 深圳华大基因研究院 | 同时确定多种样本的线粒体基因组序列信息的方法和*** |
| CN104805191B (zh) * | 2015-03-31 | 2018-02-13 | 江汉大学 | 一种测试纯系玉米新品种的特异性、一致性与稳定性的方法 |
| CN105653896B (zh) * | 2016-01-22 | 2019-02-12 | 北京圣谷同创科技发展有限公司 | 高通量测序突变检测结果验证方法 |
| CN105574365B (zh) * | 2016-01-22 | 2018-10-26 | 北京圣谷同创科技发展有限公司 | 高通量测序突变检测结果的统计学验证方法 |
| CN105907748B (zh) * | 2016-05-10 | 2017-10-13 | 广州嘉检医学检测有限公司 | 一种基于高通量测序的线粒体基因组文库及其构建方法 |
| CN106599616B (zh) * | 2017-01-03 | 2019-05-31 | 上海派森诺医学检验所有限公司 | 基于duplex-seq的超低频突变位点检测分析方法 |
| CN110310702B (zh) * | 2018-03-16 | 2021-03-23 | 深圳华大基因科技服务有限公司 | 一种基因组测序组装结果修复的方法、装置和存储介质 |
| CN112867800B (zh) * | 2018-12-17 | 2024-07-12 | Illumina公司 | 用于制备测序文库的方法和手段 |
| CN111621552B (zh) * | 2019-06-13 | 2022-05-06 | 中国科学院广州生物医药与健康研究院 | 检测mtDNA突变的方法和*** |
| WO2021142769A1 (zh) * | 2020-01-17 | 2021-07-22 | 深圳华大智造科技有限公司 | 同步进行dna正义链和反义链测序的方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103173441A (zh) * | 2013-02-05 | 2013-06-26 | 深圳华大基因研究院 | 线粒体全基因组dna扩增、引物、测序及突变检测 |
-
2013
- 2013-11-06 CN CN201310544351.4A patent/CN103602735B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103173441A (zh) * | 2013-02-05 | 2013-06-26 | 深圳华大基因研究院 | 线粒体全基因组dna扩增、引物、测序及突变检测 |
Non-Patent Citations (2)
| Title |
|---|
| Rolling Circle Amplification of Complete Nematode Mitochondrial Genomes;Sha Tang等;《Journal of Nematology》;20051231;第37卷;236-241 * |
| 基于新一代测序技术的mtDNA突变检测方法学的建立;李薇薇等;《现代生物医学进展》;20130510;第13卷(第13期);2433-2436 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103602735A (zh) | 2014-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103602735B (zh) | 利用高通量测序精确测定线粒体dna高频和低频突变的方法 | |
| US11996168B2 (en) | Systems and methods for determining relative abundances of biomolecules | |
| CN104894651B (zh) | 微量起始dna的高通量测序文库构建方法及其所构建的高通量测序文库 | |
| CN108300716B (zh) | 接头元件、其应用和基于不对称多重pcr进行靶向测序文库构建的方法 | |
| CN103060924B (zh) | 微量核酸样本的文库制备方法及其应用 | |
| EP3555305B1 (en) | Method for increasing throughput of single molecule sequencing by concatenating short dna fragments | |
| CN102409042B (zh) | 一种高通量基因组甲基化dna富集方法及其所使用标签和标签接头 | |
| CN105463585A (zh) | 基于单链dna分子构建测序文库的方法及其应用 | |
| CN105400776A (zh) | 寡核苷酸接头及其在构建核酸测序单链环状文库中的应用 | |
| Heidrich et al. | Investigating CRISPR RNA biogenesis and function using RNA-seq | |
| CN112359093B (zh) | 血液中游离miRNA文库制备和表达定量的方法及试剂盒 | |
| CN108517567B (zh) | 用于cfDNA建库的接头、引物组、试剂盒和建库方法 | |
| CN106350590A (zh) | 用于高通量测序的dna文库构建方法 | |
| Liu et al. | Research progress in molecular biology related quantitative methods of MicroRNA | |
| CN102586456A (zh) | 一种多重竞争性pcr检测拷贝数变异的方法 | |
| KR20210112350A (ko) | 다중 복제수 변이 검출 및 대립 유전자 비율 정량화를 위한 정량적 앰플리콘 서열분석 | |
| CN104093854A (zh) | 表征组合物中的rna的方法和试剂盒 | |
| CN105506747B (zh) | 富集原始转录本信息的rna文库的构建方法及其应用 | |
| CN105899682B (zh) | 用邻接标签序列复制dna用于基因组和表观基因组测序 | |
| CN104450872A (zh) | 一种高通量多样本多靶点单碱基分辨率的甲基化水平检测方法 | |
| CN114875118B (zh) | 确定细胞谱系的方法、试剂盒和装置 | |
| Gingrich et al. | Effect of RNA degradation on the data quality in quantitative PCR and microarray experiments | |
| WO2020259303A1 (zh) | 一种快速构建rna 3'端基因表达文库的方法 | |
| CN102943109A (zh) | 一种基于串联多内对照的检测拷贝数变异的方法 | |
| CN107904297B (zh) | 用于微生物多样性研究的引物组、接头组和测序方法 |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190505 Address after: 12th Floor, 926 Yishan Road, Xuhui District, Shanghai, 2003 Patentee after: Hei biotech (Shanghai) Co., Ltd. Address before: 200433 No. 220, Handan Road, Shanghai, Yangpu District Patentee before: Fudan University |
|
| TR01 | Transfer of patent right |