利用数字PCR鉴定B淋巴细胞的Ig基因分布的方法与流程

本发明属于生物
技术领域：
，涉及一种利用数字pcr鉴定b淋巴细胞的ig基因分布的方法。
背景技术：
：数字pcr(digitalpcr，dpcr)技术是一种新的核酸检测和定量方法，采用绝对定量，不依赖标准曲线或参考基因，比传统realtimepcr有更高的灵敏度。通过将微样品进行倍数稀释和分液，直至每个样品中所含有的待测分子数不会超过一个，再将所有样本在相同条件下进行pcr扩增，能够测量只含有一个目标分子的样品。可以用染料法或者taqman探针法进行pcr扩增，根据反应的荧光信号进行结果分析。由于其绝对定量的方法原理及高度灵敏检测特性，dpcr能够应用于不同的领域：拷贝数变异、突变检测、基因相对表达研究、二代测序结果验证、mirna表达分析、单细胞基因表达分析等。其优势主要表现在：1、对于稀有目标分子及突变的检测，能够精确灵敏进行。例如疾病无创(自血液、脱落细胞中)体外诊断筛查技术、体液(血液、尿液、粪便等)中肿瘤核酸标记物高分辨率检测。由于体液(血、尿液、粪便等)中或无创样本(自血液、脱落细胞中)，所含目标检测分子极少，而正常体细胞中背景信号极强，而dpcr技术能够通过微量化有效的减少背景的干扰，同时极其灵敏检测到低含量的目标分子。目前已经应用dpcr的技术有：肿瘤靶向药物基因检测、癌肿分子标记检测、无创产前等技术。2、病原微生物高灵敏度的定性及定量检测。在常规检测方法下，有极少部分病原微生物结果呈现为“弱阳性”。主要是由于这些常规技术在低拷贝检测时，由于方法限制无法检出。而通过dpcr技术则能够对该部分样本进行高灵敏度及特异性的检测，同时能够对病原微生物载量进行定量分析。3、实验室标准品绝对含量检测。分子生物学实验中，常需要对各种中间产物及其他核酸分子进行定量。而这些定量过程中使用的标准品含量通常是通过紫外发光或染料法测后倍比稀释而来，过程中的操作会极大影响标准品的准确性及后续实验。通过dpcr的高灵敏度绝对定量，可以对标准品含量进行准确检测或对样本直接进行绝对定量。成熟的b淋巴细胞抗体基因由vdj基因重排后组成，如图1所示，重排后的vdj基因后面是恒定区基因c。v基因根据功能分为fr1、cdr1、fr2、cdr2、fr3五个区域，v基因部分区域和d基因(轻链无d基因)、j基因部分区域、随机插入片段(n)组成cdr3区域。对基因组dna而言，j基因和c基因之间存在内含子，对表达的成熟信使rna而言，j基因和c基因的内含子被剪切，j基因和c基因如图1连接在一起。人免疫球蛋白(immunoglobulin，ig)的重链胚系基因定位于第14号染色体长臂，由v、d、j、c四种基因片段组成，包含40个vh基因片段、25个d基因片段、6个jh基因片段和9个ch基因片段。b细胞接触不同抗原刺激后，igc区基因的表达情况发生特异性改变，相继出现产生igm以外的ig如igg、iga或ige。ig基因将发生类别转换或者同种型转换产生多种不同类别的ig。ig的种类总共是五种，分别是iga、igm、igg、igd、ige，c基因决定ig的类别。目前，利用免疫组库技术可以扩增b淋巴细胞ig基因中的vc区域。并且对于b淋巴细胞ig重链基因来说，扩增vc区域比vj区域更有价值，c区可以反应五种ig的比例及变化情况。在免疫组库建库技术中通常必须设计多对v引物和c引物进行多重pcr或者设计多条c引物进行5’race建库，这样pcr过程中会产生偏好性，影响后续分析的准确性，特别是对五种ig比例的分析。无法确定是否产生偏好性，也无法得到具体是哪些引物产生偏好性，无法针对性地进行引物比例调整，无法验证ig基因数据分析结果是否可靠，缺乏一种评价pcr偏好性的技术手段。技术实现要素：本发明的目的是提供一种利用数字pcr鉴定b淋巴细胞的ig基因分布的方法。本发明所提供的利用数字pcr鉴定b淋巴细胞的5种ig的c区基因的拷贝数分布的方法，可包括如下步骤：(a1)根据5种ig的c区基因序列设计五套引物，每套引物中包含反转录引物、用于扩增c区基因的扩增引物和taqman探针。其中，所述反转录引物的序列与所述扩增引物中的下游引物序列既可以相同，也可以不同。(a2)取来自同一个受检者等量的5份rna作为待测样本，每份rna对应一套引物，采用所述反转录引物将所述rna反转录成cdna；采用所述扩增引物和所述taqman探针对所述cdna做数字pcr检测，从而获得5份cdna样本中相应ig的c区基因的原始拷贝数，进而计算所述受检者b淋巴细胞的5种ig的c区基因的拷贝数分布。在本发明的一个实施例中，所述五套引物具体如下：(b1)针对ige的c区基因：所述反转录引物的序列如序列表中序列1所示；所述扩增引物为由序列表中序列1和序列6所示的两条单链dna分子组成的引物对；所述taqman探针的序列如序列表中序列15所示；(b2)针对igm的c区基因：所述反转录引物的序列如序列表中序列2所示；所述扩增引物为由序列表中序列2和序列7所示的两条单链dna分子组成的引物对；所述taqman探针的序列如序列表中序列12所示；(b3)针对igd的c区基因：所述反转录引物的序列如序列表中序列3所示；所述扩增引物为由序列表中序列3和序列8所示的两条单链dna分子组成的引物对；所述taqman探针的序列如序列表中序列14所示；(b4)针对iga的c区基因：所述反转录引物的序列如序列表中序列4所示；所述扩增引物为由序列表中序列4和序列9所示的两条单链dna分子组成的引物对；所述taqman探针的序列如序列表中序列11所示；(b5)针对igg的c区基因：所述反转录引物的序列如序列表中序列5所示；所述扩增引物为由序列表中序列5和序列10所示的两条单链dna分子组成的引物对；所述taqman探针的序列如序列表中序列13所示。所述taqman探针5’端标记有荧光基团，如fam。所述taqman探针3’端标记有淬灭基团，如bqh1。在所述方法中，作为所述待测样本的rna具体可为从所述受检者外周血中提取的rna。所述方法在如下中的应用也属于本发明的保护范围：鉴定利用免疫组库测序法检测b淋巴细胞的5种ig的c区基因的拷贝数分布时所采用的多重pcr扩增引物是否存在扩增偏好性。本发明还提供了一种鉴定利用免疫组库测序法检测b淋巴细胞的5种ig的c区基因的拷贝数分布时所采用的多重pcr扩增引物是否存在扩增偏好性的方法。本发明所提供的鉴定利用免疫组库测序法检测b淋巴细胞的5种ig的c区基因的拷贝数分布时所采用的多重pcr扩增引物是否存在扩增偏好性的方法，可包括如下步骤：(c1)利用上述方法(即利用数字pcr鉴定b淋巴细胞的5种ig的c区基因的拷贝数分布的方法)检测受检者b淋巴细胞的5种ig的c区基因的拷贝数分布；(c2)利用免疫组库测序法检测所述受检者b淋巴细胞的5种ig的c区基因的拷贝数分布；(c3)根据(c1)和(c2)的检测结果按照如下确定利用免疫组库测序法检测b淋巴细胞的5种ig的c区基因的拷贝数分布时所采用的多重pcr扩增引物是否存在扩增偏好性：若(c1)所得所述受检者b淋巴细胞的5种ig的c区基因的拷贝数分布情况与(c2)所得所述受检者b淋巴细胞的5种ig的c区基因的拷贝数分布情况一致，则利用免疫组库测序法检测b淋巴细胞的5种ig的c区基因的拷贝数分布时所采用的多重pcr扩增引物不存在扩增偏好性；反之，则利用免疫组库测序法检测b淋巴细胞的5种ig的c区基因的拷贝数分布时所采用的多重pcr扩增引物存在扩增偏好性。前文所述两种方法在如下中的应用也属于本发明的保护范围：对利用免疫组库测序法检测b淋巴细胞的5种ig的c区基因的拷贝数分布时所采用的原本存在扩增偏好性的多重pcr扩增引物的用量和/或配比进行调整，以纠正原本存在的扩增偏好性，使偏好性达到最小，从而准确反应原始样本中b淋巴细胞的5种ig的c区基因表达的分布比例。本发明还请求保护如上所述的五套引物，以及这五套引物在利用数字pcr鉴定b淋巴细胞的5种ig的c区基因的拷贝数分布情况中的应用。本发明提出一种可以验证扩增c基因时，使用多对引物的pcr是否存在扩增偏好性的技术。该方法分别使用一条引物做rna反转，一对引物扩增cdna的ig基因，过程不存在多对引物扩增效率互相影响的偏好性。并且使用探针法dpcr精确定量cdna的拷贝数，数据结果可靠。可以准确的反应c基因中5种ig的c区基因的表达量，从而验证高通量免疫组库多重pcr中c区引物扩增的准确性，使测序的数据分析结果更加可信。提供了一种有效、准确度高的方法验证测序数据中5种ig的c区基因的表达量是否存在偏好性，为多重pcr中的引物是否存在偏好性提供有力的证据。ig比例的准确性有助于疾病免疫反应、疫苗评估等研究。附图说明图1为vdj基因重排。图2为dpcr鉴定b淋巴细胞ig基因分布的技术流程图。图3为微滴分析仪分析后5中ig基因的cdna拷贝数。各图依次为igha、ighm、ighg、ighd、ighe。图4为微滴式分析仪中显示的拷贝数形式。图5为免疫组库测序法得到的ig的c区基因的拷贝数分布图。图6为免疫组库测序法得到的ig的c区基因的拷贝数分布图。图7为数字pcr鉴定ig的c区基因的拷贝数分布图。具体实施方式下述实施例中所使用的实验方法如无特殊说明，均为常规方法。下述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。实施例1、利用数字pcr鉴定b淋巴细胞的ig基因分布一、技术流程如图2所示，首先根据5种ig的c区基因序列设计五套引物，每套引物中包含反转录引物、用于扩增c区基因的扩增引物和taqman探针；其中，所述反转录引物的序列与所述扩增引物中的下游引物序列相同。引物序列如表1至表3所示。表15种ig的c区基因序列的反转录引物ighercccacgcacccgagac(序列1)ighmrgtgaggtggctgcgtacttgc(序列2)ighdrccacgcatttgtactcgc(序列3)ighargggcagggcacagtcacatcc(序列4)ighgrggaaggtgtgcacgccgctggtc(序列5)表25种ig的c区基因序列的扩增引物ighefgtcttccccttgacccgct(序列6)ighmfctgtgagaattccccgtcgga(序列7)ighdfcatcatatcagggtgcagacac(序列8)ighafgtggtcrtcgcmtgcctggt(序列9)ighgfgtcttccccctggcrccct(序列10)ighercccacgcacccgagac(序列1)ighmrgtgaggtggctgcgtacttgc(序列2)ighdrccacgcatttgtactcgc(序列3)ighargggcagggcacagtcacatcc(序列4)ighgrggaaggtgtgcacgccgctggtc(序列5)表35种ig的c区基因序列的taqman探针ddhaccgctttcgctccaggtcacact(序列11)5'-fam；3'-bqh1ddhmaagccccgggtrctgctga(序列12)5'-fam；3'-bqh1ddhgaggactacttcccmgaaccggtg(序列13)5'-fam；3'-bqh1ddhdtaccacccaacgtccgtgact(序列14)5'-fam；3'-bqh1ddheccctccaatgccacctccgtgact(序列15)5'-fam；3'-bqh1然后抽取人的外周血1-2ml，提取rna。分别取5份等量的rna样本，实施例中以0.5μg为例。分别用表1中的c区特异性反转录引物反转为cdna。取相同体积的cdna做dpcr。利用bio-rad型号qx200tm的数字pcr仪，由微滴生成仪和微滴分析仪组成。采用油包水微滴式乳化技术，在一个反应管内乳化成上万个独立的pcr反应器，从而实现数字pcr技术。一个样本产生20000个微滴，每个微滴体积为1nl，仪器一次能检测96个样本，灵敏度万分之一。使用bio-rad微滴式数字pcr(dpcr)探针法试剂盒dropletpcrsupermix(货号1863024)，体系如下表4，放到微滴生成仪上自动化生成微滴。pcr运行条件为95℃10分钟，94℃30秒，56℃1分钟，98℃10分钟，4℃保存，第2和第3步循环40次，温度上升和下降梯度为2℃/s。把样本放到微滴分析仪上，用quantasoft进行分析，最后会显示每个孔内的cdna分子拷贝数。表4探针法dpcr体系试剂组分体积(μl)mix11引物(5μm)2.2探针(10μm)1.1cdna样本1-7.7水补到22μl总计22二、结果图3为微滴分析仪分析后5中ig基因的cdna拷贝数，图示依次为igha、ighm、ighg、ighd、ighe。由于人体中占大部分的免疫球蛋白是igha、ighm、ighg，ighd、ighe含量很少。于是，取相同体积的cdna中，igha、ighm、ighgcdna稀释了20倍，ighd、ighe没有稀释。样本浓度太高可能会导致ddpcr失败，可能微滴数太多，结果也会不准确。表5为5种ig基因的表达量，并按各自的拷贝数计算了5种ig基因表达量的比例分布，蓝色的点为有荧光信号的微滴数。微滴式分析仪中显示的拷贝数形式可见图4。表55种ig重链c区基因的cdna拷贝数目标区域20×cdna浓度(拷贝数/μl)cdna浓度(拷贝数/μl)百分比igha7871574070％ighm71.214246％ighg218436020％ighd无8254％ighe无14.90％另外，本发明的发明人用上述的dna样本进行免疫组库建库测序，图5和图6为免疫组库测序法得到的ig的c区基因的拷贝数分布图，图7为数字pcr鉴定ig的c区基因的拷贝数分布图。图5和图6是不同的两套引物，图6中使用的这套引物为文献lineagestructureofthehumanantibodyrepertoireinresponsetoinfluenzavaccinationningjiangetal.scitranslmed5,171ra19(2013)报道的，对ig重链的c基因分布无明显偏好的引物。图5与图7比较，图5使用的引物存在明显的扩增偏好，无法反映原始的ig基因的表达量。图6与图7比较，两者ig的c区基因的拷贝数的分布结果比较一致。<110>深圳华大基因研究院<120>利用数字pcr鉴定b淋巴细胞的ig基因分布的方法<130>gncln171539<160>15<170>patentinversion3.5<210>1<211>16<212>dna<213>人工序列<220><223><400>1cccacgcacccgagac16<210>2<211>21<212>dna<213>人工序列<220><223><400>2gtgaggtggctgcgtacttgc21<210>3<211>18<212>dna<213>人工序列<220><223><400>3ccacgcatttgtactcgc18<210>4<211>21<212>dna<213>人工序列<220><223><400>4gggcagggcacagtcacatcc21<210>5<211>23<212>dna<213>人工序列<220><223><400>5ggaaggtgtgcacgccgctggtc23<210>6<211>19<212>dna<213>人工序列<220><223><400>6gtcttccccttgacccgct19<210>7<211>21<212>dna<213>人工序列<220><223><400>7ctgtgagaattccccgtcgga21<210>8<211>22<212>dna<213>人工序列<220><223><400>8catcatatcagggtgcagacac22<210>9<211>20<212>dna<213>人工序列<220><221>misc_feature<222>(7)..(7)<223>r为g或a<220><221>misc_feature<222>(12)..(12)<223>m为a或c<400>9gtggtcrtcgcmtgcctggt20<210>10<211>19<212>dna<213>人工序列<220><221>misc_feature<222>(15)..(15)<223>r为g或a<400>10gtcttccccctggcrccct19<210>11<211>23<212>dna<213>人工序列<220><223><400>11ccgctttcgctccaggtcacact23<210>12<211>19<212>dna<213>人工序列<220><221>misc_feature<222>(12)..(12)<223>r为g或a<400>12aagccccgggtrctgctga19<210>13<211>23<212>dna<213>人工序列<220><221>misc_feature<222>(14)..(14)<223>m为a或c<400>13aggactacttcccmgaaccggtg23<210>14<211>21<212>dna<213>人工序列<220><223><400>14taccacccaacgtccgtgact21<210>15<211>24<212>dna<213>人工序列<220><223><400>15ccctccaatgccacctccgtgact24当前第1页12