CN117887812A - 高通量测序质控用文库及其制备方法和应用 - Google Patents
高通量测序质控用文库及其制备方法和应用 Download PDFInfo
- Publication number
- CN117887812A CN117887812A CN202410289937.9A CN202410289937A CN117887812A CN 117887812 A CN117887812 A CN 117887812A CN 202410289937 A CN202410289937 A CN 202410289937A CN 117887812 A CN117887812 A CN 117887812A
- Authority
- CN
- China
- Prior art keywords
- sequence
- library
- quality control
- seq
- insert
- 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.)
- Granted
Links
- 238000003908 quality control method Methods 0.000 title claims abstract description 88
- 238000012165 high-throughput sequencing Methods 0.000 title claims abstract description 47
- 238000002360 preparation method Methods 0.000 title abstract description 19
- 239000002773 nucleotide Substances 0.000 claims description 40
- 125000003729 nucleotide group Chemical group 0.000 claims description 40
- 238000000034 method Methods 0.000 claims description 16
- 238000012408 PCR amplification Methods 0.000 claims description 12
- 239000012634 fragment Substances 0.000 claims description 12
- 238000011144 upstream manufacturing Methods 0.000 claims description 12
- 238000003780 insertion Methods 0.000 claims description 9
- 230000037431 insertion Effects 0.000 claims description 9
- 238000012163 sequencing technique Methods 0.000 abstract description 86
- 108090000623 proteins and genes Proteins 0.000 abstract description 3
- 239000013612 plasmid Substances 0.000 description 26
- 108020004414 DNA Proteins 0.000 description 17
- 238000007847 digital PCR Methods 0.000 description 11
- 238000006243 chemical reaction Methods 0.000 description 8
- 238000002156 mixing Methods 0.000 description 8
- 239000011324 bead Substances 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 241000282414 Homo sapiens Species 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000007689 inspection Methods 0.000 description 4
- 238000000746 purification Methods 0.000 description 4
- 238000010257 thawing Methods 0.000 description 4
- 241000196324 Embryophyta Species 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 244000005700 microbiome Species 0.000 description 3
- 238000007481 next generation sequencing Methods 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 238000013112 stability test Methods 0.000 description 3
- GUAHPAJOXVYFON-ZETCQYMHSA-N (8S)-8-amino-7-oxononanoic acid zwitterion Chemical compound C[C@H](N)C(=O)CCCCCC(O)=O GUAHPAJOXVYFON-ZETCQYMHSA-N 0.000 description 2
- 238000011529 RT qPCR Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000003153 chemical reaction reagent Substances 0.000 description 2
- 238000010367 cloning Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- LSNNMFCWUKXFEE-UHFFFAOYSA-M Bisulfite Chemical compound OS([O-])=O LSNNMFCWUKXFEE-UHFFFAOYSA-M 0.000 description 1
- 102000053602 DNA Human genes 0.000 description 1
- 102000016928 DNA-directed DNA polymerase Human genes 0.000 description 1
- 108010014303 DNA-directed DNA polymerase Proteins 0.000 description 1
- 241000588724 Escherichia coli Species 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 108091028043 Nucleic acid sequence Proteins 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000005119 centrifugation Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000009509 drug development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002547 new drug Substances 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
- XEBWQGVWTUSTLN-UHFFFAOYSA-M phenylmercury acetate Chemical compound CC(=O)O[Hg]C1=CC=CC=C1 XEBWQGVWTUSTLN-UHFFFAOYSA-M 0.000 description 1
- 239000013600 plasmid vector Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000007480 sanger sequencing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000002798 spectrophotometry method Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明涉及基因测序技术领域,尤其涉及一种高通量测序质控用文库及其制备方法和应用。本发明提供的高通量测序质控用文库在每个碱基位置的碱基占比均是均衡的。将本发明提供的高通量测序质控用文库与待测样本文库混合后共同上机测序,可消除因待检样本文库碱基序列的不平衡性所导致的测序质量降低,提高对不平衡文库测序的准确性。
Description
技术领域
本发明涉及基因测序技术领域,尤其涉及一种高通量测序质控用文库及其制备方法和应用。
背景技术
高通量测序技术又称下一代测序技术(Next-generation sequencing,NGS),一次能对几十万到几百万条核酸分子进行序列测定,是目前发展最为迅速的技术之一,其已深入生命科学的各个领域,在肿瘤检测、新药研发等方面前景广阔。高通量测序技术通量大、灵敏度高等优点已使其成为目前临床分子检测中最重要的应用平台之一。
高通量测序仪在测序时会采用前几十个测序Cycle信号来确定簇的位置和密度等信息,这些信息影响了测序信号采集、测序碱基质量。碱基不平衡文库(即每个测序Cycle中A、G、C、T四种碱基的含量远远偏离25%)在测序时会导致测序仪荧光信号采集失真,大大降低碱基识别的准确性,进而影响测序结果的准确性。尤其在甲基化DNA测序中,常通过亚硫酸氢盐使DNA中未甲基化的C转化为U,进而在PCR过程中转化为T,这种转化使待测序的DNA中C碱基占比大大降低,很容易导致测序中碱基的不平衡。
为解决这一问题,目前主要采用混合质控文库与待检样本文库共同上机测序的方式。但是现有质控文库主要通过非人源DNA(如PhiX噬菌体DNA)来制备,通过片段筛选方式获得一定长度范围的文库,最终混入待检样本文库中一同测序。这种方法的主要问题是:制备过程中涉及片段筛选、获得的质控文库非固定序列容易导致批间存在差异。
有鉴于此,特提出本发明。
发明内容
为解决上述技术问题,本发明提供了一种高通量测序质控用文库及其制备方法和应用,可有效降低因待测文库碱基不平衡导致的碱基测序质量下降问题。
具体的,本发明的技术方案如下:
第一方面,本发明提供了一种高通量测序质控用文库,由包括序列插入片段的模板DNA经PCR扩增得到;所述序列插入片段包括序列插入片段a、序列插入片段b、序列插入片段c和序列插入片段d;所述序列插入片段a的核苷酸序列如SEQ ID NO.1所示;所述序列插入片段b的核苷酸序列如SEQ ID NO.2所示;所述序列插入片段c的核苷酸序列如SEQ IDNO.3所示;所述序列插入片段d的核苷酸序列如SEQ ID NO.4所示。
序列插入片段a:
CGAGGATTCAACAAGGTGATCGAAGTCGACAGATCCAGAGAGACGGGCTTCAAAGCTGCGTGACGACGCTTGCGAGTCCGTATCAATATCCTCACAATAAGCACACGTGACCGTTGGTTGAACAGCACAGGACGAGCTGACCAGAAGCACTATTATATCGCTTAACGGCTCTTGAGCCAGTGTGCGTTACCTTGCAGCAATCGAGGCCGTCCGTTAATTCCTCTTGCATTCATATCGCGTATTGTTGTCTCTGTACGCGCTTACTTGGATCAGGATGACATAGCTTCTTACAGGAGCGTC
序列插入片段b:
ACTTAGACTGCAGCTAGTCGATGGTCGACTGATCGTGACTCCTGATATCATGCCTAACGCGATGTGTAGCATGCTAGTACGTGACTGCATGCATGCTAGCTAGTCACGTCAGTCATCAGTGCATCGTACTCGATGCTGTCGTCAGTAGCTACAGGACAGATCGCTGACGATGCATAGTCACTGATCGATGGCATGCATGCATAGCATGTCAGTCGTGACTGAGACTACGATGCTAGCATGCAGTCGACAGTCAGCTAGTAGACTGACTGCACTATCTCACCTAGCGTAGCTCACGTACAA
序列插入片段c:
GACATCGAGTGTTGCTACTAGCCTAGTCTATCGAAGCTGCTACTCACAGGACTTCGGATACTCACTGTAGCATACCAGTACGCTGGCTGATAGCTGCGTGATCGTGTCAGTACGCATGCACTGCATGTGCATGCTAGACTAGGCTCTAGAGTCACGTGATGACTGACTAGACGTCTAAGTACACGTACGTTAGCATCGCTGACTATAACAGTAGACTCGATGAGACTGACGCAGCTATCAGCACGACGTCAGTCGATACGACGACGTATGTTCTCACTGTAGCTACAGCTGGTATCGTCG
序列插入片段d:
TTGCCTCGACTGCTACCAGCTATCCAATGGCTCGTATCTAGTGATCTGACGTGGATCTATACGTACACTAGCATGTCAGTACAGTCAGCGAGTGATGCCTCGTAGTAACTGTAAGCACACTGTGTACGTATCTACTACAGTATTCGCTTGCGGCTCGCTCAGAGCTTATCGAACGCTGTCGACTAACGCAAGCATGTATGCGTCTCGTAGTACACGCGAGACTCGAGTCGATGCGATGACTGCAACTAGAGACATGCTATCGTGACACCAGGACGGAGTGGCTAGAGCAGATCTCATAGT
上述4种序列插入片段的长度均为300bp,其各自的4种碱基AGCT含量均相同(四种碱基各占25%),且这4种特定序列插入片段对应的相同位置上DNA碱基(A/G/C/T)均不相同。4种特定序列插入片段均无法比对到人、常见动植物和微生物,可有效防止数据污染。
优选地,所述高通量测序质控用文库包括4组目标特定序列;所述目标特定序列包含双端接头、双端Index和核心区序列。所述核心区序列由所述序列插入片段扩增得到;4组目标特定序列的核心区序列长度相同,且在序列插入片段a、序列插入片段b、序列插入片段c和序列插入片段d上对应的位置相同;所述核心区序列的长度为100-300bp。
本发明以上述4种序列插入片段a/b/c/d或包含上述4种序列插入片段a/b/c/d的质粒或长链DNA分子为模板,经PCR扩增,可获得含有4种目标特定序列的文库。目标特定序列文库由4种不同序列的文库组成,文库两端包含双端接头、双端index及特定序列(特定序列来自上述序列插入片段a/b/c/d)。4种目标特定序列文库的核心区的长度均相同(100~300bp,优选150~250bp,优选250bp)。4种目标特定序列文库的序列均是单一且固定的。4种目标特定序列各自的4种碱基AGCT含量均相同(四种碱基各占25%),且这4种目标特定序列对应的相同位置上DNA碱基(A/G/C/T)均不相同,因此可使用在同一个测序Cycle中,这4种文库测出的序列各不相同。4种目标特定序列文库需与待检样本文库的长度接近,避免差异过大,以使目标特定序列文库与待检样本文库在测序仪上进行无差异反应(例如通过桥式扩增生成簇的过程)。另外,本发明提供的高通量测序质控用文库含有双端标签(index)序列,测序数据下机后可直接进行数据拆分,用于质控数据的详细评估(包括但不限于质控数据量在总数量中的占比、测序质控用文库测序结果与已知测序结果的对比等)。
举例说明本发明涉及的4种目标特定序列文库组成结构,包含但不限于下述结构序列:
illumina测序平台测序质控用文库结构序列:
5'-AATGATACGGCGACCACCGAGATCTACAC
-Index2(CGCATGTC)
-ACACTCTTTCCCTACACGACGCTCTTCCGATCTT
-序列插入片段a/b/c/d
-AAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
-Index1(TCATTCGA)
-ATCTCGTATGCCGTCTTCTGCTTG-3'。
华大测序平台测序质控用文库结构序列:
5'-CTCTCAGTACGTCAGCAGTT
-Index2(CGCATGTC)
-CAACTCCTTGGCTCACAGAACGACATGGCTACGATCCGACTTT
-序列插入片段a/b/c/d
-AAAGTCGGAGGCCAAGCGGTCTTAGGAAGACAA
-Index1(TCATTCGA)
-CTGATAAGGTCGCCATGC-3'。
本发明优选等摩尔浓度混合4种目标特定序列文库(根据文库的质量浓度和分子量计算摩尔浓度并按照等摩尔浓度混合),同时采用数字PCR进行质检。质检合格后,加入一定浓度的poly(dN)(N=A/C/T/G碱基),制成测序质控用文库。Poly(dN)优选Poly(dA),Poly(dA)的浓度为0.1~100ng/μL,优选1~10ng/μL,优选1ng/μL;dA碱基数为5~100,优选10~50,优选18~30。
将测序质控用文库与待检样本文库按一定比例混合后通过高通量测序仪测序。测序质控用文库在总文库(测序质控用文库与待检样本文库之和)中的占比为1~50%,优选5~20%,优选5~10%,优选5%。
第二方面,本发明提供了所述高通量测序质控用文库的制备方法,由包括序列插入片段的模板DNA经PCR扩增得到;所述序列插入片段包括序列插入片段a、序列插入片段b、序列插入片段c和序列插入片段d;所述序列插入片段a的核苷酸序列如SEQ ID NO.1所示;所述序列插入片段b的核苷酸序列如SEQ ID NO.2所示;所述序列插入片段c的核苷酸序列如SEQ ID NO.3所示;所述序列插入片段d的核苷酸序列如SEQ ID NO.4所示。
本发明提供的方法可快速、稳定制备出含有目标特定序列且固定长度的可拆分测序质控用文库。
优选地,用于PCR扩增的引物对的上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ ID NO.6所示。此时,4种文库所含目标特定序列长度为300bp,仅需一轮PCR和一轮纯化即可完成制备。
作为另一种可选的实施方式,用于PCR扩增的引物对包括如下(1)、(2)、(3)、(4)和(5)所示引物对:
(1)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ IDNO.7所示;
(2)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ IDNO.8所示;
(3)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ IDNO.9所示;
(4)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ IDNO.10所示;
(5)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ IDNO.11所示。
当用于PCR扩增的引物对选自(1)、(2)、(3)、(4)和(5)时,获得的4种文库所含目标特定序列长度为160bp,共需两轮PCR和纯化完成制备。其中,(1)、(2)、(3)和(4)所示引物对用于第一轮PCR扩增,(5)所示引物对用于第二轮PCR扩增。
基于相同的原理,其他未列出的能特异性针对本发明所述序列插入片段a/b/c/d,在相同碱基位置进行PCR扩增,且能获得相同技术效果(即4种目标特定序列各自的4种碱基AGCT含量均相同,且这4种目标特定序列对应的相同位置上DNA碱基均不相同)的引物组合也在本发明的保护范围之内。
第三方面,本发明提供了所述高通量测序质控用文库或者所述制备方法制备得到的高通量测序质控用文库在高通量测序中的应用。
优选地,本发明所述应用将所述高通量测序质控用文库与待检样本文库混合,然后通过高通量测序仪测序;其中,所述高通量测序质控用文库中的目标特定序列与待测样本序列长度相同或相差±100bp以内。
第四方面,本发明还提供了一种高通量测序文库,包括所述高通量测序质控用文库,还包括待测样本文库;所述高通量测序质控用文库在所述高通量测序文库中的占比为1-50%,优选为5-20%,更优选5-10%,更优选为5%。
本发明优选根据数字PCR得到的拷贝浓度进行分析,若所述高通量测序质控用文库中的4种目标特定序列文库两两之间的拷贝浓度的相对偏差不超过±10%即可判定为合格。如果不合格,则根据数字PCR结果进一步调整文库浓度,直至合格。
有益效果:
本发明提供了一种高通量测序质控用文库及其制备方法和应用。本发明提供的高通量测序质控用文库在每个碱基位置的碱基占比均是均衡的。将本发明提供的高通量测序质控用文库与待测样本文库混合后共同上机测序,可消除因待检样本文库碱基序列的不平衡性所导致的测序质量降低,提高对不平衡文库测序的准确性。
具体实施方式
本发明提供了一种高通量测序质控用文库及其制备方法和应用。所述高通量测序质控用文库由包括序列插入片段的模板DNA经PCR扩增得到。
在本发明提供的一种具体实施方式中,所述高通量测序质控用文库的制备及应用方法包括以下步骤:
(1)合成一组4种含有特定序列插入片段a/b/c/d的质粒:质粒上的插入片段包括测序文库的双端接头、双端index及特定序列插入片段。其中4种特定序列插入片段的长度均为300bp,4种特定序列插入片段所含的4种碱基(A/G/C/T)的比例均相同(四种碱基各占25%),且这4种特定序列插入片段对应的相同碱基位置上的碱基(A/G/C/T)均不相同。4种特定序列插入片段均无法比对到人、常见动植物和微生物,可有效防止数据污染。
(2)通过上述质粒获得一组4种含有目标特定序列的文库:这组文库包括双端接头、双端index及核心区序列(核心区序列来自上述质粒含有的特定序列插入片段)。4种文库的核心区长度均相同(100~300bp,优选150~250bp,优选250bp),且与待检样本文库的长度接近。
(3)等摩尔浓度混合上述4种文库:根据质量浓度和分子量计算摩尔浓度并按照等摩尔浓度混合,同时采用数字PCR进行质检。质检合格后,加入一定浓度的poly(dN)(N=A/C/T/G碱基),制成测序质控用文库。Poly(dN)优选Poly(dA),Poly(dA)的浓度为0.1~100ng/μL,优选1~10ng/μL,优选1ng/μL;dA碱基数为5~100,优选10~50,优选18~30。
(4)将测序质控用文库与待检样本文库按一定比例混合后通过高通量测序仪测序。测序质控用文库在总文库(测序质控用文库与待检样本文库之和)中的占比为1~50%,优选5~20%,优选5~10%,优选5%。
本发明提供的高通量测序质控用文库通过利用AGCT四种平衡碱基以平衡整体待测样本文库的碱基平衡性,可有效解决由于待测文库碱基不平衡可能导致的测序准确性降低问题,提高测序质量。相较于现有技术,本发明仅采用4种含有目标特定序列且固定长度的测序质控用文库。最重要的是这4种目标特定序列中各自的4种碱基AGCT含量均相同(四种碱基各占25%),且这4种目标特定序列对应的相同碱基位置上DNA碱基(A/G/C/T)均不相同,因此可使在同一个测序Cycle中,这组质控文库测出的AGCT碱基含量均相同(四种碱基各占25%),从而避免由于待测样本文库碱基不平衡而导致的测序质量下降,大大改善整体测序的质量。本发明在碱基不平衡的文库样本测序时,可以加入较大比例测序质控用文库,以此来抵消待测文库样本的碱基不平衡性。在正常待测文库中也可以少量的加入测序质控用文库,以此来监控测序质量。同时,本发明提供的高通量测序质控用文库还具有稳定性高,便于使用和保存,制备方法简单快捷、稳定可靠,且无需文库片段筛选、批间差异小等优点。
由于不同测序平台对应的测序文库的接头和index可能不同,因此以下以illumina测序平台为例说明实施方式。
该制备方法主要步骤如下:
S1、采用基因合成法,合成4个含有特定序列的插入片段(即序列插入片段a/b/c/d)的大肠杆菌克隆质粒(或长链DNA分子)。
以质粒为例说明实施方式:质粒载体为克隆性质粒,包括但不限于pUC57质粒。质粒上的插入片段包括测序文库的双端接头、双端index及特定序列插入片段a/b/c/d。质粒插入片段为:5'-AATGATACGGCGACCACCGAGATCTACAC-Index2(CGCATGTC)-ACACTCTTTCCCTACACGACGCTCTTCCGATCTT-序列插入片段a/b/c/d-AAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC-Index1(TCATTCGA)-ATCTCGTATGCCGTCTTCTGCTTG-3'。4种特定序列插入片段a/b/c/d均无法比对到人、常见动植物和微生物,可有效防止数据污染。通过Sanger测序法确认质粒插入片段序列合成正确。获得的4种特定序列文库质粒对应编号为1/2/3/4。
S2、采用分子生物学常规方法进行质粒提取纯化(例如使用商品化试剂盒)、测定质粒的质量浓度(例如使用分光光度计、荧光计等)。
S3、通过PCR获得目标特定序列文库:
I、当4种文库所含目标特定序列长度为300bp时,仅需一轮PCR和一轮纯化即可完成制备。
(1)PCR
将制备目标特定序列文库的上述4种质粒以及下表中试剂取出,室温解冻,完全融化后轻柔振荡混匀,短暂离心后置于冰盒上备用。
在0.2mL的离心管中,按照下表分别配制4管PCR反应体系。其中质粒模板用量为:1pg~50ng,优选1pg~1ng,优选10pg。使用高保真DNA聚合酶,以确保DNA序列在扩增过程中的保真性。
表 1 PCR反应体系
4种质粒模板均可使用相同的PCR引物:
PCR正向引物(P5):AATGATACGGCGACCACCGAG
PCR反向引物(P7):CAAGCAGAAGACGGCATACGAGA
表 2 PCR反应条件
注:模板用量越少,循环数应对应增加。
(2)根据产品说明书,使用磁珠(VAHTS DNA Clean Beads)纯化PCR产物。
II、当4种文库所含目标特定序列长度为100~300bp时,需两轮PCR和纯化完成制备。下文以160bp的目标特定序列为例进行说明。
(1)第1轮PCR
反应体系同表1,4种质粒模板均可使用相同的正向引物和各自对应的反向引物:
PCR正向引物(P5):AATGATACGGCGACCACCGAG
PCR反向引物(对应质粒1):GTTAAGCGATATAATAGTGCGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTT
PCR反向引物(对应质粒2):CAGCGATCTGTCCTGTAGCTGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTT
PCR反向引物(对应质粒3):TCAGTCATCACGTGACTCTAGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTT
PCR反向引物(对应质粒4):AGCTCTGAGCGAGCCGCAAGGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTT
反向引物中的下划线序列对应于质粒上的特定序列,且均位于4个质粒对应的相同位置,以确保PCR扩增出的产物长度一致,且在4个PCR产物中的特定序列的每个碱基位置都各不相同。
PCR程序同表2。
(2)根据产品说明书,使用磁珠(VAHTS DNA Clean Beads)纯化第1轮PCR产物。
(3)第2轮PCR
反应体系同表1,PCR模板为4种纯化后的第1轮PCR产物1μL。
4种质粒模板均可使用相同的正向引物和反向引物:
PCR正向引物(P5):AATGATACGGCGACCACCGAG
PCR反向引物:CAAGCAGAAGACGGCATACGAGAT-Index1(TCGAATGA)-GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTT
PCR程序同表2。循环数可减少至5~10轮。
(4)根据产品说明书,使用磁珠(VAHTS DNA Clean Beads)纯化第2轮PCR产物。
S4、制备测序质控用文库:
(1)通过分光光度计或荧光计对4种目标特定序列文库的质量浓度进行测定,并根据下述公式计算出对应的摩尔浓度:
摩尔浓度(nM)=质量浓度(ng/μL)÷分子量(g/mol)10^-6
根据上述公式计算出的4种目标特定序列文库的摩尔浓度,统一稀释到相同的摩尔浓度,摩尔浓度为:5~200nM,优选5~100nM,优选10nM。
(2)通过数字PCR对稀释后的文库进行质检。以Thermo QuantStudio 3D数字PCR平台为例:
表 3 数字PCR反应体系
数字PCR反应体系中PCR正向引物(P5):AATGATACGGCGACCACCGAG;PCR反向引物(P7):CAAGCAGAAGACGGCATACGAGA。
表 4 数字PCR反应条件
根据数字PCR得到的拷贝浓度进行分析,若4种目标特定序列文库两两之间的拷贝浓度的相对偏差不超过±10%即判定为合格。如果不合格,则根据数字PCR结果进一步调整文库浓度,直至合格。
(3)向已满足上述等摩尔浓度混合标准的4种目标特定序列文库中,加入一定浓度的poly(dN),制成测序质控用文库。Poly(dN)优选Poly(dA),Poly(dA)的浓度为0.1~100ng/μL,优选1~10ng/μL,优选1ng/μL;dA碱基数为8~100,优选10~50,优选18~30。
S5、将测序质控用文库与待检样本文库按一定比例混合后通过高通量测序仪测序。测序质控用文库在总文库(测序质控用文库与待检样本文库之和)中的占比为1~50%,优选5~20%,优选5~10%,优选5%。
S6、下机后数据可通过文库标签将测序质控用文库数据拆分出来,单独进行分析,不影响待检样本文库的结果分析。同时,也可通过序列比对,将测序质控用文库的已知序列从测序结果中剔除,避免影响待检样本测序结果分析。
为了便于理解本发明,下面将对本发明实施例中的技术方案进行清晰完整的描述,但是不能把它们理解为对本发明保护范围的限定。实施例中未注明具体条件的,按照常规条件或制造商建议的条件进行。所用试剂和仪器未注明生产厂家的,均为可以通过市售购买获得的常规产品。
实施例1
本实施例按照具体实施方式所述进行illumina测序平台质控用文库制备,选择3天制备,每天制备3个批次,每个批次制备1管质控用文库。摩尔浓度稀释至10nM。
将制备好的测序质控用文库使用qPCR方法(KAPA Library Quantification Kit)进行定量。结果显示:所有制备出的质控文库的实测浓度与理论浓度的相对偏差均不超过10%,9批次质控文库的浓度的相对标准偏差(CV)仅为2.41%,表明按本发明方法制备出的测序质控用文库制备的批间稳定性好,具体结果见下表。
表 5 测序质控用文库多批次制备定量浓度
将上述制备好的9个批次测序质控用文库分别与9组待测文库样本,均按5%的比例混合后在illumina Nextseq CN500测序仪上完成测序,针对第1组待测文库样本设置对照组:即不加入测序质控用文库、仅对该组待测文库样本进行测序。结果显示:对照组测序数据中,簇密度(Cluster Density)出现异常,仅为146 K/mm2(正常范围在170~230 K/mm2),碱基质量值≥Q30占比也较低;对应的该组待测文库样本以及其他组待测样本文库在加入测序质控用文库的所在上机批次均符合要求。同时数据下机后使用其特有Index标签进行拆分后,测序质控用文库实际产出占比与测序质控用文库掺入比例相差极少(相对偏差在5%以内),得到的整体下机数据的簇密度(Cluster Density),簇通过率(Cluster PassingFilter),碱基质量值≥Q30占比均符合要求。结果表明:按本发明方法制备出的测序质控用文库制备的批间差异很小,且用其与待检样本文库混合上机测序得到的数据质量高。具体结果见下表。
表 6 测序质控用文库使用效果
实施例2
本实施例按照具体实施方式所述进行制备测序质控用文库,摩尔浓度为10nM,其中加入Poly(dA)的浓度为1ng/μL。设立对照组,对照组中质控文库摩尔浓度为10nM,但不加入poly(dA)。两组测序质控用文库同时做实时稳定性测试(保存于-30~-10℃,在第0/4/8/12/16/20/24个月检测文库浓度)、加速破坏稳定性测试(在反复冻融10次、室温(20~30℃)放置2天后检测文库浓度)。共制备3个批次测序质控用文库及对应的对照组进行测试。使用qPCR方法(KAPA Library Quantification Kit)进行定量,定量结果显示,按照具体实施方式所述制备出的测序质控用文库,无论是实时稳定性,还是加速破坏稳定性,各次检测结果间无明显差异;对照组(不加入poly(dA))在对应各次测试中的稳定性均表现较差。结果表明:按本发明方法制备出的测序质控用文库制备的实时稳定性、加速破坏稳定性(反复冻融、室温放置)均良好。具体结果见下表。
表 7 测序质控用文库实时稳定性测试
表 8 测序质控用文库加速破坏稳定性测试
以上所述实施例仅表达了本发明的几种实施方式,便于具体和详细地理解本发明的技术方案,但并不能因此而理解为对本发明保护范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。
Claims (10)
1.高通量测序质控用文库,其特征在于,由包括序列插入片段的模板DNA经PCR扩增得到;所述序列插入片段包括序列插入片段a、序列插入片段b、序列插入片段c和序列插入片段d;所述序列插入片段a的核苷酸序列如SEQ ID NO.1所示;所述序列插入片段b的核苷酸序列如SEQ ID NO.2所示;所述序列插入片段c的核苷酸序列如SEQ ID NO.3所示;所述序列插入片段d的核苷酸序列如SEQ ID NO.4所示。
2.根据权利要求1所述高通量测序质控用文库,其特征在于,所述高通量测序质控用文库包括4组目标特定序列;所述目标特定序列包含双端接头、双端Index和核心区序列;所述核心区序列由所述序列插入片段扩增得到;4组目标特定序列的核心区序列长度相同,且在序列插入片段a、序列插入片段b、序列插入片段c和序列插入片段d上对应的位置相同;所述核心区序列的长度为100-300bp。
3.根据权利要求2所述高通量测序质控用文库,其特征在于,还包括poly(dN),N=A/C/T/G碱基。
4.根据权利要求3所述高通量测序质控用文库,其特征在于,所述poly(dN)选用poly(dA),所述poly(dA)的浓度为0.1-100mg/μL,dA碱基数为5-100。
5.权利要求1-4任意一项所述高通量测序质控用文库的制备方法,其特征在于,由包括序列插入片段的模板DNA经PCR扩增得到;所述序列插入片段包括序列插入片段a、序列插入片段b、序列插入片段c和序列插入片段d;所述序列插入片段a的核苷酸序列如SEQ ID NO.1所示;所述序列插入片段b的核苷酸序列如SEQ ID NO.2所示;所述序列插入片段c的核苷酸序列如SEQ ID NO.3所示;所述序列插入片段d的核苷酸序列如SEQ ID NO.4所示。
6.根据权利要求5所述制备方法,其特征在于,用于PCR扩增的引物对的上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ ID NO.6所示。
7.根据权利要求5所述制备方法,其特征在于,用于PCR扩增的引物对包括如下(1)、(2)、(3)、(4)和(5)所示引物对:
(1)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ ID NO.7所示;
(2)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ ID NO.8所示;
(3)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ ID NO.9所示;
(4)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ ID NO.10所示;
(5)上游引物核苷酸序列如SEQ ID NO.5所示,下游引物核苷酸序列如SEQ ID NO.11所示。
8.权利要求1-4任意一项所述高通量测序质控用文库或者权利要求5-7任意一项所述制备方法制备得到的高通量测序质控用文库在高通量测序中的应用。
9.根据权利要求8所述应用,其特征在于,将所述高通量测序质控用文库与待检样本文库混合,然后通过高通量测序仪测序;其中,所述高通量测序质控用文库中的目标特定序列与待测样本序列长度相同或相差±100bp以内。
10.高通量测序文库,其特征在于,包括权利要求1-4任意一项所述高通量测序质控用文库,还包括待测样本文库;所述高通量测序质控用文库在所述高通量测序文库中的占比为1-50%;所述高通量测序质控用文库中的4种目标特定序列两两之间的拷贝浓度相对偏差≤±10%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410289937.9A CN117887812B (zh) | 2024-03-14 | 高通量测序质控用文库及其制备方法和应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410289937.9A CN117887812B (zh) | 2024-03-14 | 高通量测序质控用文库及其制备方法和应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117887812A true CN117887812A (zh) | 2024-04-16 |
CN117887812B CN117887812B (zh) | 2024-07-09 |
Family
ID=
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108203847A (zh) * | 2016-12-16 | 2018-06-26 | 深圳华大智造科技有限公司 | 用于二代测序质量评估的文库、试剂及应用 |
CN110628890A (zh) * | 2019-11-07 | 2019-12-31 | 中国人民解放军军事科学院军事医学研究院 | 测序质控标准品及其应用与产品 |
WO2020191521A1 (zh) * | 2019-03-22 | 2020-10-01 | 深圳华大智造科技有限公司 | 核酸序列、rna目标区域测序文库的构建方法及应用 |
CN113621609A (zh) * | 2021-09-15 | 2021-11-09 | 深圳泛因医学有限公司 | 文库构建引物组及其在高通量检测中的应用 |
CN114703266A (zh) * | 2022-04-06 | 2022-07-05 | 杭州杰毅医学检验实验室有限公司 | 一种内参、包含有该内参的试剂盒及其应用 |
CN116179664A (zh) * | 2023-02-10 | 2023-05-30 | 青岛锐翌精准医学检验有限公司 | 基于内参确定微生物的高通量检测方法和系统及试剂盒 |
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108203847A (zh) * | 2016-12-16 | 2018-06-26 | 深圳华大智造科技有限公司 | 用于二代测序质量评估的文库、试剂及应用 |
WO2020191521A1 (zh) * | 2019-03-22 | 2020-10-01 | 深圳华大智造科技有限公司 | 核酸序列、rna目标区域测序文库的构建方法及应用 |
CN110628890A (zh) * | 2019-11-07 | 2019-12-31 | 中国人民解放军军事科学院军事医学研究院 | 测序质控标准品及其应用与产品 |
CN113621609A (zh) * | 2021-09-15 | 2021-11-09 | 深圳泛因医学有限公司 | 文库构建引物组及其在高通量检测中的应用 |
CN114703266A (zh) * | 2022-04-06 | 2022-07-05 | 杭州杰毅医学检验实验室有限公司 | 一种内参、包含有该内参的试剂盒及其应用 |
CN116179664A (zh) * | 2023-02-10 | 2023-05-30 | 青岛锐翌精准医学检验有限公司 | 基于内参确定微生物的高通量检测方法和系统及试剂盒 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110129415B (zh) | 一种ngs建库分子接头及其制备方法和用途 | |
CN106048009B (zh) | 一种用于超低频基因突变检测的标签接头及其应用 | |
CN111808854B (zh) | 带有分子条码的平衡接头及快速构建转录组文库的方法 | |
WO2019062289A1 (zh) | 一种探针及其适用于高通量测序的对目标区域进行富集的方法 | |
CN114107459B (zh) | 一种基于寡核苷酸链杂交标记的高通量单细胞测序方法 | |
CN110628890A (zh) | 测序质控标准品及其应用与产品 | |
WO2022099794A1 (zh) | 一种质检寡核苷酸交叉污染的方法 | |
CN112795654A (zh) | 用于生物体融合基因检测与融合丰度定量的方法及试剂盒 | |
CN116024308A (zh) | 肺癌相关基因高通量扩增子文库的制备方法、多重pcr引物对及应用 | |
CN110724731A (zh) | 一种在多重pcr体系内加入内参定量核酸拷贝数的方法 | |
CN117887812B (zh) | 高通量测序质控用文库及其制备方法和应用 | |
CN110872609A (zh) | 对小rna分子精准建库和测序的方法以及应用 | |
WO2021047175A1 (zh) | 基因拷贝数定量分析 | |
CN117887812A (zh) | 高通量测序质控用文库及其制备方法和应用 | |
CN115807056B (zh) | 一种bcr或tcr重排序列模板池及其应用 | |
CN103789414A (zh) | 17个x染色体短串联重复序列的复合扩增试剂盒 | |
CN114277114B (zh) | 一种扩增子测序添加唯一性标识符的方法及应用 | |
CN116064818A (zh) | 检测igh基因重排及超突变的引物组、方法和系统 | |
CN117867086B (zh) | 高通量测序文库定量用标准品及其制备方法和应用 | |
CN109266723A (zh) | 稀有突变检测方法、其试剂盒及应用 | |
CN117867090B (zh) | Dna片段化水平定量用标准品及其制备方法和应用 | |
CN114807302B (zh) | 扩增子文库构建方法及用于地中海贫血突变型与缺失型基因检测的试剂盒 | |
CN114277096B (zh) | 鉴别地中海贫血αααanti4.2杂合型和HKαα杂合型的方法和试剂盒 | |
CN116626189B (zh) | 一种基于LC-MS对On-DNA化学反应中DNA损伤的分析方法 | |
US20220298505A1 (en) | Method for constructing pacbio sequencing library |
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 |