一种利用扩增DNA片段长度多态性测定短链RNA的方法.pdf

本发明公开了一种利用扩增DNA片段长度多态性测定短链RNA的方法，属于分子生物技术领域，包括以下步骤：首先利用至少两种与待测短链RNA相比无天然同源序列的合成小RNA作为测定的内标，将这些合成小RNA以不同分子数混合，形成动态小RNA标准分子梯度；再以等量的所述动态小RNA标准与所述待测短链RNA混合，经由RNA反转录、cDNA加尾、PCR同步扩增及PCR产物DNA长度多态性片段的荧光定量分析，测得所述待测短链RNA扩增产生的DNA片段荧光强度在所述动态小RNA标准荧光强度梯度上的相对比例，本方法可以实现对短链RNA尤其是miRNA的绝对定量，且具有高特异性、高灵敏度、结果客观、重复性好，可在100-106拷贝数范围内准确定量，也适用于RNA粗提物。

1.  一种利用扩增DNA片段长度多态性定量测定短链RNA的方法，其特征在于，包括以下步骤：
首先利用至少两种与待测短链RNA相比无天然同源序列的合成小RNA作为测定的内标，将这些作为内标的所述合成小RNA以不同分子数混合，形成动态小RNA标准分子梯度；
再以等量的所述动态小RNA标准与所述待测短链RNA混合，经由RNA反转录、cDNA加尾、PCR同步扩增及PCR产物DNA长度多态性片段的荧光定量分析,测得所述待测短链RNA扩增产生的DNA片段荧光强度在所述动态小RNA标准荧光强度梯度上的相对比例，从而实现所述待测短链RNA的相对定量。

2.  如权利要求1所述的一种利用扩增DNA片段长度多态性定量测定miRNA的方法，其特征在于，测得所述待测短链RNA在所述动态小RNA标准梯度上的相对强度比例后，以待测短链RNA序列为模板合成短链RNA参照物，测定不同分子数量的待测RNA参照物在上述动态小RNA标准梯度上的相对比例，由此获得待测短链RNA参照物的分子数目-相对强度的校正曲线，再利用所述待测短链RNA相对强度比例通过所述校正曲线，计算出样品中待测短链RNA的绝对分子数量。

3.  如权利要求1或2所述的一种利用扩增DNA片段长度多态性定量测定短链RNA的方法，其特征在于，所述待测短链RNA为miRNA或siRNA。

4.  如权利要求1或2所述的一种利用扩增DNA片段长度多态性定量测定短链RNA的方法，其特征在于，所述RNA反转录过程中，所采用的引物为欧米茄引物。

5.  如权利要求1或2所述的一种利用扩增DNA片段长度多态性定量测定短链RNA的方法，其特征在于，所述RNA反转录过程中，所采用的引物为茎环引物。

6.  如权利要求5所述的一种利用扩增DNA片段长度多态性定量测定短链RNA的方法，其特征在于，所述茎环引物为经过长度编码的茎环引物。

7.  如权利要求6所述的一种利用扩增DNA片段长度多态性定量测定短链RNA的方法，其特征在于，所述长度编码的方法为：在所述茎环引物的PCR靶位点与探针序列之间加上不同碱基数目，并调整引物5’端的碱基序列使之维持茎环的二级结构不变。

8.  如权利要求2所述的一种利用扩增DNA片段长度多态性定量测定短链RNA的方法，其特征在于，所述校正曲线符合对数回归的方式，表述为：aX^b，其中a、b为常数项，经对不同分子数目的所述合成小RNA的实际测定值而确定。

9.  如权利要求1所述的一种利用扩增DNA片段长度多态性定量测定短链RNA的方法，其特征在于，所述与待测短链RNA相比无天然同源序列的合成小RNA的种数为三种。

一种利用扩增DNA片段长度多态性测定短链RNA的方法
技术领域
本发明涉及分子生物技术领域，特别涉及一种用扩增DNA片段长度多态性定量测定短链RNA的方法。 
背景技术
microRNA (即miRNA) 是一类调控基因表达的短单链小分子核糖核酸，由22个左右核苷酸（nt）组成。miRNA参与细胞的各种活动，影响并参与维持细胞的生理和分化状态，特定miRNA表达量的变化还能反映生物的生理、病理特征，预示其作为鉴定多种疾病的生物标记的可能性。例如肝细胞所特有的miR-122，在血液中的水平变化就可以揭示肝脏受损的情形。miR-122为肝脏细胞所独有，并在每个人类肝细胞中有13万个，具有目前用于肝损伤鉴定的生物标记物所欠缺的专一性及灵敏度。尽管miRNA在这些病变中所起的作用还有待明确，机制有待进一步的揭示，但这已成为目前miRNA研究的一个最重要的方向。其中miRNA定量检测方法的客观性、灵敏度及标准化是将miRNA研究成果推向临床疾病诊断应用的瓶颈和关键。 
已知人类拥有大约1千种以上的miRNA，与各种特定生理病理现象相关的miRNA表达量的动态范围很大，从个位数级别到上万，加上miRNA拥有超短的长度、相似的序列、二级构象、末端修饰及长度差异性等特性，给miRNA的纯化和准确定量带来很大的挑战。 
目前应用较多的对miRNA进行定量测定的主流技术包括定量反转录PCR（RT-qPCR）及同时测定多种miRNA的能力的miRNA深度测序法（miR-seq）和microArray杂交测定法等。其它也有很多基于RNA/DNA杂交技术的高通量技术平台，如NanoString nCounter、nanopore 技术，但受制于这些方法本身的不完善，缺乏应有的灵敏度和准确性，距离实际应用尚需时日。 
上述miRNA测定方法实际应用的主要问题是灵敏度，如其中最灵敏的NanoString nCounter法的测定低限为1 amol即约6x10⁵个分子以上的水平，而miR-seq和microArray的测定低限为0.1 fmol即6x10⁷个分子以上的水平，显然不满足大多数分子诊断的需求。 
其中miR-seq法，对miRNA的鉴定不依赖于miRNA序列，可以对miRNA新种的发现有不可替代的作用。而且其能统计出miRNA的相对读数，而有很多研究中利用它作了很多病理诊断应用方面的探索，发现了很多有益的苗头。但是，miR-seq的样品处理需要利用RNA Ligase、DNA Ligase制备适用的cDNA库，并经PCR扩增后才能上机检出，作到相对定量。冗长的样本制备过程及多种修饰酶的利用，破坏了待测目标数量与检测读值的线性相关性，影响miRNA测定的客观性。因此，如果没有更好的制备方法，miR-seq还不能作为miRNA定量测定的方法。 
反转录实时荧光PCR定量技术（RT-qPCR）是现代分子生物诊断方法中对RNA检出最灵敏、最客观可靠的方法，而且方法快捷，不需要特殊的设备。传统的RT-qPCR测定对象是长链的RNA，如mRNA，所用primers长度通常在20nt以上。但是，成熟miRNA的总长度才22nt左右,测定miRNA的RT-qPCR必须经过特殊的修饰，使之能适合qPCR。除茎环引物可作miRNA的RT以外，锁核苷酸引物(LNA)及多聚A聚合酶加尾-寡聚dT引发的RT也被应用于miRNA的定量测定。以利用TaqMan探针、茎环引物的反转录-qPCR方法为代表的实时荧光PCR定量法具备有测定单拷贝miRNA分子的能力，已经广泛的用于miRNA的机制研究、疾病诊断、治疗预后等的研究中。因其是所有方法中最灵敏、客观、容易应用，是验证从microArray、miR-seq方法获得的miRNA表达谱的优选方法。尽管不断的技术改进和对工艺的强化不断地催生新的miRNA qPCR检测方法，但是，至少有两大障碍阻止了这些技术在临床应用中的接受度：方法质量的评估及方法的标准化。这两个障碍严重影响到方法结果在不同批次、操作人员及实验室之间的重复性，即测定结果存在不一致不稳定的缺点。 
对液态的生物标本而言，miRNA拷贝数的绝对定量分析可以揭示miRNA在一定体积中的数量，是一种更容易理解和符合自然的单位，尤其适合液体样本。拷贝数的绝对定量分析有助于更准确地理解miRNA与疾病之间的关系，更能精细的揭示病变程度及预后。Bissel等用10倍稀释的合成miRNA等摩尔混合物作模板，利用microArray杂交测定法做出miRNA拷贝数的效正曲线，做到了对miRNA的绝对定量。并实际测定大鼠肝脏RNA中miR-122的含量为每10pg 约53000个拷贝。很遗憾的是该法的最低测定的限度为1 amol，也即6x10⁵拷贝，不能满足临床诊断的实际需求。实际上，因为不同microArray平台技术的原理均以DNA/RNA杂交来实现的，杂交信号没有扩增，受杂交的热力学平衡的制约，在满足测定的专一性的前提下，方法灵敏度很难突破amol的低限。另外一个可能带来差异的因素是miRNA的标记反应受miRNA的末端序列、二级构象及所用的标记方法所限，很难做到等比例地标记所有种类的miRNA，对不同miRNA的测定带来歧视性偏差。 
实用的miRNA定量测定方法必须能作到灵敏、准确和可标准化。临床样品具多样性，及成分的不确定性。一个简单的例子是来自不同个体的血清所含的成分对反转录反应及PCR反应可能没有作用、抑制或促进。现行的miR-seq、RT-qPCR法及microArray法用到RNA ligase、Reverse Transcriptase及Taq polymerase对RNA样品的杂质很敏感，加上这些技术设计原理是试图直接得到样品中待测目标的数量，而这种测定方式静态地假定纯化RNA为测定对象，没有设计抗衡RNA样品中不纯物影响的机制，使其所得的结果受样品来源影响而呈现随机性，增加了其标准化的难度。因而从严格意义上来说，绝对定量是不可能用现有测定方法来实现的。对于血液等生物样本来讲，从不同个体血清纯化出的RNA中所含的杂质对RT反应及PCR扩增效率具有抑制甚或促进的作用，加大了使用传统RT-qPCR方法进行定量测定的难度。另外样本的处理和纯化，一般采用盐酸胍苯酚等蛋白质变性剂抽提方法得到纯度很高的RNA，但由于miRNA的平均大小仅为22nt，需要高浓度乙醇长时间沉淀以避免丢失，使纯度和产率这对矛盾更加突出，而用亲和树脂柱的纯化方法一般能保证RNA的纯度但尚需证明纯化后RNA的完整性，即miRNA的丢失和选择性丢失的问题。从病人获得的病理标本通常是少量而不可重复的，因此，这些问题一直困扰着将miRNA定量检测技术推向临床应用的努力。 
发明内容
本发明的发明目的在于：针对上述存在的问题，提供一种改善短链RNA定量测定的灵敏度、准确度和标准化的新方法。 
本发明采用的技术方案是这样的：一种利用扩增DNA片段长度多态性定量测定miRNA的方法，包括以下步骤： 
首先利用至少两种与待测短链RNA相比无天然同源序列的合成小RNA作为测定的内标，将这些作为内标的所述合成小RNA以不同分子数混合，形成动态小RNA标准分子梯度；
再以等量的所述动态小RNA标准与所述待测短链RNA混合，经由RNA反转录、cDNA加尾、PCR同步扩增及PCR产物DNA长度多态性片段的荧光定量分析,测得所述待测短链RNA扩增产生的DNA片段荧光强度在所述动态小RNA标准荧光强度梯度上的相对比例，从而实现所述待测短链RNA的相对定量。
在本发明中，首先在样品中加入不同分子数目作为动态小RNA标准的合成小RNA与待测短链RNA混合，同一反应中进行同步检测，由同步测得的分子标准量为动态定量的标准尺度，获得待测短链RNA的相对峰值。为了不与样品中的待测短链RNA的序列相混，所以作为内标的合成小RNA要求与待测短链RNA无天然同源序列。杂质或影响反应的因素对样品中的待测短链RNA及动态小RNA标准分子均有相同程度的影响，这种设计思路明显地区别于现行的方法，使本发明的方法拥有足够的灵敏度和准确性的同时，兼具有对样品不纯成分的一定包容能力，有益于不同操作人员及实验室之间不同样品测定结果的比较，容易实施标准化。从而可以抗衡RNA尤其是短链RNA纯化、样本质量及仪器、操作、耗材等外在因素对RNA尤其是短链RNA测定结果的影响。 
作为优选的技术方案，测得所述待测短链RNA在所述动态小RNA标准梯度上的相对强度比例后，再以待测短链RNA序列为模板合成的短链RNA参照物，该短链RNA参照物的序列与待测短链RNA的序列相同，其不同分子数量，在上述动态小RNA标准为梯度上的相对比例而获得的分子数目-相对强度校正曲线，利用所述待测短链RNA相对强度比例通过所述校正曲线，计算出样品中待测短链RNA的绝对分子数量。 
对比用已知数量的合成目标RNA与相同小分子RNA标准的相对测定值作出的校正曲线而得到确定的待测RNA的绝对数目，从而可以获得待测RNA的绝对定量结果。 
作为优选的技术方案，所述待测短链RNA为miRNA或siRNA。 
本发明的方法可用于对各种生物样本或纯化RNA中miRNA、siRNA及RNA的种类（定性）及含量（定量）的测定，siRNA具有与miRNA一样的物理、化学特性，上述的方法同样适用于siRNA的定性定量测定。本发明也适用于鉴定RNA中的一段或几段较小的序列，从而达到各种不同RNA的定性、定量测定的目的，尤其适合于miRNA这类短链RNA。 
越来越多的研究已经揭示miRNA及小RNA表达水平变异对疾病的诊断及预后的联系，从病人获得的病理标本通常是少量而不可重复的，miRNA及小RNA的纯化方法还不能满足现行鉴定方法的需求，本发明可利用抵消miRNA样本中不纯物的影响而对待测miRNA进行精确定量，可以满足通过对miRNA的相对或绝对定量测定达到临床诊断、预后、刑侦等的需求。miRNA指纹谱（miRNA signature）是几个到几十个与某生命现象相关的miRNA在特定样本中的表达水平，其改变规律可揭示相关的生理或病理变化。本发明的实施，可以在一个或几个反应里，方便、快速、灵敏和准确地鉴定miRNA指纹谱，有助于miRNA指纹谱的研究和在临床诊断中的应用。 
本发明设计了通过DNA片段长度多态性荧光相对定量分析来完成RNA尤其是miRNA的分子数绝对定量检测的方法（miRNA-derived Fragment Length Polymorphism Assay，本发明中称为miRFLP定量分析法）。miRFLP定量分析法利用3种无天然同源序列的合成小RNA作为测定的内标，将这三种内标小RNA以不同分子数混合，形成动态小RNA标准分子梯度，再以等量的动态小RNA标准与待测短链RNA混合，经由RNA反转录、cDNA加尾、PCR同步扩增及PCR产物的荧光片段长度多态性分析等四个步骤测得待测短链RNA产生的DNA片段荧光强度在动态小RNA标准荧光梯度上的相对比例。用合成的待测短链RNA参照物分子数与其测得的荧光相对强度即可绘出短链RNA绝对数量的换算曲线，由此计算出待测短链RNA的绝对含量。 
作为优选的技术方案，所述RNA反转录过程中，所采用的引物为欧米茄引物。 
本发明所述的欧米茄引物，尤其指申请号为PCT/CN2013/070525的PCT专利申请中所公开的欧米茄引物。欧米茄引物作为本发明miRFLP定量分析反应的引物探针。欧米伽引物含有一个极稳定的二级茎环结构，可避免引物二聚体的形成，提高末端引发精确度，且将引物有效地分隔为探针区及编码区。通过加减编码区的碱基数目就可以方便地对不同引物进行长度编码，并可在同一反转录反应中将多种miRNA转换为长度不一的cDNA片段,经PCR同步扩增进行荧光片段长度多态性分析，解决了单一反应中检测多种miRNA的关键技术问题。 
miRFLP反应由miRNA反转录、cDNA加尾、PCR同步扩增及PCR产物的荧光片段长度多态性分析等四个步骤组成，见图1。反应的第一步是将miRNA和欧米伽引物混合杂交，miRNA与互补的探针形成配对，然后用反转录酶以未配对的miRNA 3’末端为模板合成cDNA。除去RNA后，新合成的cDNA与含有共同PCR靶位的3’寡居核苷酸引物杂交并用DNA聚合酶以cDNA的3’末端为起始补齐DNA配对后的单链缺口。经此组合后，正确装配的cDNA拥有相同的5’及3’末端序列，可经一对荧光PCR引物进行同步等比扩增，扩增产物通过毛细管凝胶电泳分离后测定不同长度PCR片段的荧光强度，完成miRFLP分析图谱。 
作为优选的技术方案，所述RNA反转录过程中，所采用的引物为茎环引物。 
作为进一步优选的技术方案，所述茎环引物为经过长度编码的茎环引物。 
作为更进一步优选的技术方案，所述长度编码的方法为：在所述茎环引物的PCR靶位点与探针序列之间加上不同碱基数目，并调整引物5’端的碱基序列使之维持茎环的二级结构不变。 
茎环结构的引物探针（Applied Biosystems, Inc, 如PCT/CN2013/070525所提及的）将多个小RNA同时定量地转换为cDNA，经与上述miRFLP反应路径相似的cDNA产物加尾、反应后经PCR同步扩增进行荧光片段长度多态性分析。实验结果证实，只要加上长度可变的编码区，茎环引物也可用于本发明的设计规划的实施，并在测定有限数目的miRNA目标时，可以与欧米伽引物互换或同时使用。 
作为优选的技术方案，所述校正曲线符合对数回归的方式，表述为：aX^b，其中a、b为常数项，经对不同分子数目的所述合成小RNA的实际测定值而确定。 
作为优选的技术方案，所述与待测短链RNA相比无天然同源序列的合成小RNA的种数为三种。由于荧光强度和荧光物质的数量不是线性的关系，需要至少3种来确定一元二次方程的回归，而采用三种合成小RNA，既可以达到定量测定的要求，同时检测过程也较经济。 
miRFLP定量分析法中的每一步中间产物的实现和信号的扩增都是以线性或接近线性定量的方式进行的，因此扩增后DNA片段的荧光强度与目标miRNA的分子数成线性定量的关系。但受DNA分析仪荧光探头的线性定量范围的客观限制，仪器探头的荧光强度与荧光物质的数量间的响应曲线适用于实测得到的待测目标的荧光强度与其分子数之间的换算。为确定仪器探头的响应曲线，我们选择了一份具有多个杂峰片段的PCR产物，经系列稀释后将样品用ABI 3730xl型DNA分析仪进行荧光定量分析，分析确定了稀释倍数与实测的荧光强度的相关性在25–25000 FU的测定范围内均符合一元二次方程的回归关系，R²平均值大于0.9999。本方法中所用的动态miRNA标准和样本目标荧光信号的扩增是以接近线性的定量方式同步进行的，一元二次回归方程适用于待测miRNA对动态miRNA标准的相对荧光强度的计算。现代的DNA测序仪具有极高的敏感度和动态探测范围（ABI 3730DNA分析仪：5-31000 FU）。不同厂家或不同型号的DNA荧光定量分析仪均可采用上述方法对适用的荧光响应曲线及范围进行标定，遴选出最适合亦即最精确的回归方式。 
由于各种外在因素对测定方法的影响，所测得的miRNA值还不是miRNA的绝对分子数目，而是待测miRNA在动态miRNA标准的荧光强度梯度上的比值，比值的大小反映出待测miRNA分子数目与动态miRNA标准分子数的比例。因待测miRNA与动态miRNA标准在同一反应中检测，同步扩增，排除了外在因素的干扰，动态miRNA标准的分子数目用量是恒定的，待测miRNA相对荧光强度比值就能客观地反映不同反应中待测miRNA的多寡。利用系列稀释的已知分子数的合成待测miRNA为参照物与动态miRNA标准的比较可以作出待测miRNA的相对比值与待测miRNA分子数目的校准曲线，将待测miRNA相对比值转换为其miRNA绝对拷贝数。用二倍或三倍稀释的合成miRNA参照物还可以用来准确地衡量方法的测定范围及在每个浓度点的大致误差范围，可以帮助提高所得的数据分析结果的可信度和标准化。用IBM SPSS Statistics 20分析软件对10个实际测得的三倍系列稀释的数据进行回归模型的分析，我们发现miRNA的相对荧光强度与分子数之间符合对数方程的回归规律(power regresssion)，拟合优度R²在适用的测定范围内大于0.99。校正miRNA参照物易于标准化，可用于校正不同操作员之间的操作误差，让不同实验室得出的数据具有可比性。 
由于不同miRNA的探针及3’适配寡聚核苷酸具有不同的Tm值，即使在相同的反应条件下，其热动力学的影响也导致不同miRNA拥有不同的校准曲线。由于没有能通用的miRNA的绝对定量校准曲线，每个待测miRNA均需有自己的绝对定量标准曲线。本发明适用于多个miRNA的同时鉴定，多种miRNA合成参照物也可以混合使用，大大减少所需反应的数目。例如完成9个不同miRNA的校正曲线，以三倍稀释制作一个覆盖50-5x10⁵个分子测定范围的标准曲线需要测定10个不同浓度的标准miRNA参照物，如果使用标准的RT-qPCR测定法，仅完成标准曲线的制作就需要至少90个反应才行，而用本发明的多目标同步测定法可将9种合成miRNA混合后作为模板，9种miRNA的绝对定量校准曲线仅需10个反应就可完成。 
综上所述，由于采用了上述技术方案，本发明的有益效果是：本发明的miRFLP定量分析法与现行方法相比，至少有如下优点： 
1.绝对定量：得益于与样本RNA在同一反应管中反应形成的动态小RNA标准分子，消除复杂样品以及RNA纯度对测定的干扰，可对待测miRNA进行绝对定量；
2.高特异性：可鉴别只有1个碱基差别或同一家族的miRNA；
3.高灵敏度：对细胞单拷贝miRNA测定的上样量范围为 0.1μg-1 ng 总RNA；
4.测定范围：在100-10⁶ 拷贝数范围内准确定量，超出此范围可作相对定量；
5.结果客观：设立了内控做为RT或PCR反应失败导致阴性结果的指标，对基因片段的甑别可以精确到±0.2个碱基；
6.重复性好：可控的合成miRNA标准参照保证操作人员及实验室之间的重复性和可比性，有益于方法的标准化。
7.适用于RNA粗提物：可直接从血清、血浆等病人样本中准确客观地测出miRNA的表达水平，重复好，结果可靠。 
附图说明
图1A和图1B是miRFLP定量分析法反应流程及测定原理图； 
图2是随机选取的具有多个杂峰的荧光PCR产物经25倍稀释后的DNA片段长度荧光分析图谱；
图2中，纵坐标FU为荧光强度，“□”标注的为随机选择的DNA片段；
图3是长度为62.00nt及91.22nt(B)的DNA片段在不同稀释度下的荧光测定值与稀释倍数关系符合一元二次方程回归曲线；
图3中，“A”为62.00nt的DNA片段，“B”为91.22nt的DNA片段，纵坐标FU为荧光强度单位；
图4A和图4B是动态miRNA标准和不同浓度miRNA的miRFLP分析图谱；
图4A和图4B中，“std 1”、“std 2”及“std 3”分别代表动态miRNA标准（std 1、std 2 和std 3）的荧光DNA片段， “1”、“2”、“3”、“4”分别代表反应1、反应2、反应3、反应4，“a”代表miR-92a的荧光DNA片段，“b”代表miR-92b的荧光DNA片段，纵坐标为荧光强度；
图5A和图5B是不同分子数miR-92a和miR-92b的miRFLP测定图谱及标准曲线；
图5A和图5B中，“1”为只含有动态miRNA标准的空白对照测定结果，“2-10”分别为加入了2.5x10⁵个至244个合成miR-92a分子和等量的miR-92b为待测目标，图5B为不同miRNA拷贝数与其相对荧光强度的点对应图；
图6是miRFLP方法测定miR-25家族成员矫正曲线的重复性及测定范围结果图；
图6中，横坐标为miRNA分子数，纵坐标为相对荧光强度；
图7是Let-7家族miRFLP定量分析谱对单个Let-7成员的测定结果图；
图7中，从上至下分别为let-7b、let-7c、let-7d及let-7g的测定结果，“std 1”、“std 2”及“std 3”代表的含义与图4A中的“std 1”、“std 2”及“std 3”一致；
图8是用茎环引物的miRFLP分析方法对miR-92b及miR-25的测定结果；
图9带有以碱基数目为编码的茎环引物；
图9中，A：为探测Standard 1 RNA设计的茎环引物，在5’端PCR靶位（箭头）与探针之间放置了4个碱基（方框）作为标记。茎环dG = -13.71。
B：为探测Standard 2 RNA设计的茎环引物，在5’端PCR靶位（箭头）与探针之间放置了7个碱基（方框）作为标记。茎环dG = -13.09。 
具体实施方式
下面对本发明作详细的说明。 
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。 
实施例1 miRFLP定量分析法反应流程及测定原理 
本实施例以欧米伽引物为例，具体阐述miRFLP定量分析法反应流程及测定原理。
miRFLP反应由miRNA反转录、cDNA加尾、PCR同步扩增及PCR产物的荧光片段长度多态性分析等四个步骤组成，如图1A和图1B所示。反应的第一步是将miRNA和欧米伽引物混合杂交，miRNA与互补的探针形成配对，然后用反转录酶以未配对的miRNA3’末端为模板合成cDNA。除去RNA后，新合成的cDNA与含有共同PCR靶位的3’寡居核苷酸引物杂交并用DNA聚合酶以cDNA的3’末端为起始补齐DNA配对后的单链缺口。经此组合后，正确装配的cDNA拥有相同的5’及3’末端序列，可经一对荧光PCR引物进行同步等比扩增，扩增产物通过毛细管凝胶电泳分离后测定不同长度PCR片段的荧光强度，完成miRFLP分析图谱，图1B为整个miRFLP定量分析的最后一步。 
  
实施例2 PCR荧光标记产物等倍稀释液荧光强度的测定试验
荧光定量分析仪测出的荧光强度与待测荧光物质的数量成正比关系，这种关系与仪器所配置的荧光探头有关，即不同的荧光探头拥有不同的荧光响应曲线。ABI的Prizma 310型DNA测序仪测得的荧光强度在5-7000 FU的范围里与待测的荧光物质数量成线性定量的关系，当荧光强度值大于7000FU时改变为抛物线。不同仪器的荧光探头响应特性不一样，可以影响到测出的荧光强度与待测荧光量的回归关系。本实施例对ABI 3730xl型DNA分析仪的荧光响应曲线用2倍系列稀释的PCR产物进行了标定。实验随机选取了一个具有多个不同片度且荧光强度变化较大的荧光PCR产物，用1xTE以25、50、100、200及400倍稀释后用ABI 3730xl型DNA分析仪分析。图2是该PCR产物经25倍稀释的样品荧光扫描图，从中选取11个峰值介于253-25000FU的DNA片段，统计出它们在不同稀释度的荧光读值（见表1），经IBM SPSS Statistics 20 统计软件对稀释倍数和荧光强度测定值的对应关系进行回归曲线的优化评估，显示一元二次方程回归曲线适用于描述所有片度的荧光强度与稀释倍数的关系，回归的拟合优度R平方值均大于0.999。因此，经ABI 3730xl型DNA分析仪测得的DNA片段的荧光强度与荧光物质数量的对应关系可以用一元二次方程回归曲线来精确计算，而不是简单的线性定量关系。重复实验证明同样的样品在不同批次的测定绝对值可以有所不同，但这种一元二次方程的回归模型并不受影响。图3显示了在低丰度片段（A：30FU-312FU）及高丰度片段（B：2052FU-24949FU）经系列稀释后的回归曲线图。不同厂家或不同型号的DNA荧光定量分析仪均可采用上述方法对适用的荧光响应曲线及范围进行标定，遴选出最适合亦即最精确的回归方式。
表1：用于图2的PCR产物经25、50、100、200和400倍稀释后用ABI 3730xl型DNA分析仪在不同稀释倍数下测得的荧光值及经IBM SPSS Statistics 20统计软件作一元二次方程回归统计后的拟合优度R平方： 
实施例3 miR-92a和miR-92b的miRFLP定量分析测定试验
miRFLP分析法将待测miRNA与动态miRNA标准混合均匀，经过miRNA反转录、cDNA加尾修饰及荧光PCR同步扩增处理，最后用DNA测序仪进行DNA片段长度和荧光定量的分析。具体做法，首先配制 4 μl杂交母液：2 μl 5xRT 缓冲液、1 μl 10 mM MgSO₄、1 μl 动态miRNA标准的混合物，该混合物包括分子数为3×10⁶的“标准1 RNA”（即图4A中的“std1”，下同）、分子数为×10⁵的“标准2 RNA”（即图4A中的“std2”，下同）、分子数为3×10⁴的“标准3 RNA”（即图4A中的“std3”，下同），动态miRNA标准的RNA序列见表2，经Sangers miRbase20版数据库所做的序列搜索证明无同源序列）。用4 μl杂交母液同2 μl待测RNA样本混合，加入2 μl 10 nM的欧米伽探针混合物（含：针对std1，2，3 及miR-92a和miR-92b的探针，探针序列见表2）。混合均匀后进行杂交，杂交反应条件为：55C 10分钟，然后以每分钟一度的速度从55C降到20C。加入2 μl 反转录酶混合液（1 μl MMTV 反转录酶，0.5 μl dNTP，0.5 μl 水，Takara），混匀后25C保温30分钟，37C10分钟，85C变性5分钟。取5 μl RT产物加入15 μl 3’ 适配引物缓冲液（ 10 μl JumpStart Taq ReadyMix，2 μl 100 nM 3’寡聚核苷酸适配引物, 1 μl 0.1 ug/μl RNase A，2 μl无菌水， Sigma)，反应条件：95C 2分钟，60c 10 分钟，重复5次55c 1分钟到30c 5分钟的循环，20c保温30分钟,以每分钟1C的速度从42c升到68c, 72c保温5分钟。取5 μl加入由15 μl JumpStart Taq ReadyMix（Sigma），0.5 uM PCR引物和10 μl水配成25 μl的PCR反应液。进行荧光PCR扩增，反应条件为：95℃ 2分钟，40次经95℃ 10秒，68℃ 3分钟，72℃ 30秒的循环。通用PCR引物为：[5Fam]GTGCTGAGTCACGAGGTATTCTA 及 CACCGACAGGAGACCTGTTCT（购自GenScript）。反应完成后，将PCR产物以1:20或1:50的稀释度稀释后用ABI 3730xl型DNA分析仪进行荧光片段长度多态性分析。miRFLP分析图谱结果见图4A和图4B所示；
图4A中，反应1为只含有动态miRNA分子标准的空白对照测定结果；反应2中加入了2.5x10⁵个miR-92a分子和等量的miR-92b为待测目标，图4B中，反应3中加入了3.1x10⁴个miR-92a和等量的miR-92b为待测目标，反应4中加入了78 pg H1299细胞总RNA。
在图4A和图4B的miRFLP分析图谱中，不同长度的DNA片段代表不同的miRNA分子（峰注上面的数值为片段长度，nt)，其荧光强度反映出该miRNA的相对数量（峰注下面的数值为荧光强度，FU），同一反应中不同动态miRNA标准测得的荧光强度形成与其分子数相应的荧光强度梯度。不同反应中按比例放入的各种动态miRNA标准的绝对数目是一致的，但因诸多原因，如操作、仪器、耗材、试剂、测定条件等的不同质性，均影响到荧光强度的绝对测定值，并显示出不小的差异（如反应1和反应3中的std 1的峰值)。而因待测目标的稀释倍数较小，反应2中用于RNA保存的保存液浓度要高于反应3，影响到RT或PCR反应的效率，引起动态miRNA标准荧光强度梯度的变动。因待测目标与动态miRNA标准反应前已经混匀，这些影响反应效率的因素对待测目标和动态miRNA标准的影响是相同的，待测目标产生的荧光强度在动态标准荧光强度梯度的相对比例可以用来客观地衡量待测目标与动态miRNA标准之间的分子数比。 
表2： miRFLP定量分析法测定miR-92a及miR-92b的欧米伽引物（其碱基序列分别为SEQ ID NO.1、SEQ ID NO. 2；采用的欧米欧米伽引物的碱基序列均为SEQ ID NO.3；采用的3’ 寡聚核苷酸适配引物（即表2中的“3’Adapter”）的碱基序列分别为SEQ ID NO.4、SEQ ID NO.5）、动态miRNA标准（Std 1，2，3）(其碱基序列分别为SEQ ID NO.6、SEQ ID NO.7、SEQ ID NO.8；分别采用的欧米欧米伽引物的碱基序列分别为SEQ ID NO.9、SEQ ID NO.10、SEQ ID NO.11，采用的3’ 寡聚核苷酸适配引物的碱基序列均为SEQ ID NO.12)及成分构成： 表中：dG, dC, dA and dT的平均分子量（308.95道尔顿）计为1nt， 5' 荧光基团分子量: 474.5 

实施例4 miRFLP定量分析校准曲线及细胞总RNA中miR-92a和miR-92b的定量测定试验
miRFLP定量分析反应对同一反应里的小RNA的反转录、修饰及扩增反应均以线性等比的方式进行，动态标准荧光强度梯度可以排除外在因素对测定的影响，但探针的Tm、PCR片段长度及构象的差异仍然能影响DNA片段的荧光强度。通过设计和反应条件的优化可以把后二者的影响降低到可以忽略的程度，但探针Tm对待测目标产生的荧光强度的影响不可忽视。因而需要用已知数目的待测目标作出分子数目与相对荧光强度之间的校准曲线才能将待测目标的相对荧光强度转换为待测目标的实际分子数。用系列稀释的合成miRNA作为待测目标，可获得相对荧光强度与目标分子数的转换曲线。
本实施例利用等量混合的合成miR-92a及miR-92b核酸小分子以2倍的稀释度稀释后为待测目标，遵循实施例3中的miR-92的miRFLP分析法的测定条件，测定了9个不同稀释度的合成miR-92a和miR-92b的相对荧光强度。实验组设置了对3.12ng、0.312ng、0.0312ng A549细胞RNA，0.222ng Hela细胞RNA及0.25ng H1299细胞RNA的测定反应。每个测定反应均设置三次重复实验，以判定测定的误差范围。图5A中 左边的数字显示每个反应中miR-92a和miR-92b的分子数目。右边为miRFLP测定图谱，显示miR-92a及miR-92b荧光片段的绝对荧光强度随待测分子数的降低而逐渐下降。用Microsoft Excel表格对动态miRNA标准的荧光强度和对应的分子数作回归统计，求得最佳的一元二次方程，并在此方程中代入待测目标荧光强度即得出待测目标相对于动态小RNA标准的的相对荧光强度。图5B为不同miRNA拷贝数与其相对荧光强度的点对应图，其趋势符合乘方回归的规律，miR-92b及miR-92a的分子数校准回归曲线分别为：0.0001×^7.4237 和 0.0013×^6.4765，拟合优度R平方值大于0.98。表3列出了三种细胞RNA的荧光强度测定结果以及利用各自反应中动态小RNA标准的荧光强度将测得的miR-92b及miR-92a的荧光强度转换成相对荧光强度（RFU）值。分别利用上述miR-92a、miR-92b的分子校准曲线的乘方回归方程，可以得出miRNA在细胞RNA中的含量及各自的误差范围。对3个10倍稀释的A549细胞RNA的测定结果显示在0.03ng-3ng 的RNA上样范围内，miR-92b及miR-92a的测定误差分别为10.13%和12.63%，显示出miRFLP分析法对miRNA的定量测定受样品上样量的影响很小，可靠性大。 
表3利用miRFLP方法测定细胞总RNA中miR-92a及miR-92b含量。细胞总RNA浓度由Qubic 2.0荧光定量仪测定。 
实施例5 miRFLP定量分析法的定量测定范围及误差验证试验
利用等量混合的合成miR-2、miR-92a及miR-92b核酸小分子以3倍的稀释度稀释后为待测目标，遵循实施例3中的miR-92的miRFLP分析测定条件，利用列于表5中的引物及动态小RNA分子标准测定9个不同稀释度的合成miR-25、miR-92a和miR-92b的相对荧光强度。每个反应中所含合成miRNA分子数分别为：250,000、83,333、27,778、9,259、3,086、1,029、343、114、38个等量混合的miR-25、miR-92a和miR-92b合成小RNA以及空白对照。每个测定反应均设置三次重复实验，以判定测定的误差范围。用Microsoft Excel表格对每个反应中动态小RNA标准的荧光强度和对应的分子数作回归统计，求得最佳的一元二次方程，并在此方程中代入待测目标荧光强度即得出待测目标相对于动态小RNA标准的的相对荧光强度。图6显示不同miRNA分子数与其测定的相对荧光强度的散点对应图及对数回归的拟合优度。miR-92b、miR-92a及miR-25的分子数校准回归曲线分别为：4.3965×^3.4854、9.9139×^3.0082及4.4131×^3.3164，拟合优度R平方值大于0.98。表4列出了三组重复实验所得不同分子数水平的相对荧光强度的误差范围（C.V） 显示出miR-25的miRFLP分析法的有效定量检测范围为38-250,000个，在38个分子的水平上的检测误差范围为104%，并可在此检测范围外对待测目标作出相对定量的测定。miRNA作为理想的生物标记物的一个重要原因是由于受生理病理的影响，miRNA数量的变动范围很大，具有很高的诊断灵敏性。miRFLP分析法显示出的miRNA定量测定范围、检测误差水平及客观性，均超出目前所有的鉴定方法，满足临床对miRNA分析的要求，具有很好的应用前景。
表4 miRFLP分析法测定miR-25家族各个成员在不同分子数水平的误差范围 

表5 miRFLP分析法测定miR-25家族成员的欧米伽引物（Omega Primer）、动态miRNA标准（Std 1，2，3）及3’ 寡聚核苷酸适配引物（3’Adapter）的碱基序列及成分构成表

实施例6 miRFLP定量分析法的方法专一性验证试验
遵循实施例3中的miRFLP分析测定条件，利用列于表6中的引物及动态miRNA分子标准测定的部分let-7 miRFLP分析谱分别测定let-7b、let-7c、let-7d及let-7g的相对荧光强度。每个反应中所加入的作为模板的合成miRNA分子数量均为12500个。
其中，let-7b的序列（SEQ ID NO.13）为：ugagguaguag guugugugguu  
      let-7c的序列（SEQ ID NO.14）为：ugagguaguag guuguaugguu
      let-7d的序列（SEQ ID NO.15）为：agagguaguag guugcauaguu
      let-7g的序列（SEQ ID NO.16）为：ugagguaguag uuuguacaguu
所得的miRFLP图谱见图7。从图7中可以看出，miRFLP测定法具有很高的测定专一性，可以区分单个碱基差异的miRNA分子。用12500个Let-7不同成员作模板，用Let-7家族miRFLP测定谱对单个Let-7成员分别进行独立的测试，结果证实了Let-7成员间的最大交叉反应在Let-7b和Let-7f之间，并小于5%，与ABI公司的茎环引物的专一性相似。但不同的是，本发明中交叉反应的测定是用含有能探测所有Let-7成员的欧米伽混合探针，也即在多目标测定的反应中实现的，而非在单一引物中进行，这种对专一性的评估方式更贴近实际应用。
表6： miRFLP定量分析法测定部分let-7家族成员（let-7b、let-7c、let-7d、let-7g）的欧米伽引物（碱基序列分别为SEQ ID NO.17、SEQ ID NO.18、SEQ ID NO.19、SEQ ID NO.20，采用的3’ 寡聚核苷酸适配引物（3’Adapter）的碱基序列均为SEQ ID NO.21）动态miRNA标准（Std 1，2，3）及成分构成 

实施例7利用茎环引物的miRFLP定量分析测定方法
利用表7列出的茎环引物及动态小RNA分子标准，按照实施例3中的miR-92的miRFLP分析测定条件，测合成miR-25和miR-92b的相对荧光强度。分析图谱正确地显示出代表目标miRNA的DNA片段，并且片段的荧光强度与目标miRNA的使用量成正比例，说明经过一定的优化，茎环引物也可用于miRNA的miRFLP分析法。图8中的miRFLP分析谱显示利用茎环引物对miR-92b及miR-25的定量测定结果。图8中，代表miR-92b的正确DNA片段为94.93nt。右图为miR-25的miRFLP分析谱，代表miR-25的正确DNA片段长度为80.74nt。茎环引物miRFLP分析法测出的信号荧光强度与不同miRNA模板的数量之间关系与欧米伽引物miRFLP分析法结果相似，茎环引物也可用来代替欧米伽引物以同样的方式测出miRNA的绝对定量值。
茎环引物的原始设计用途是用来启动miRNA的反转录反应，初衷是利用qPCR扩增循环数测定目标miRNA的浓度，用荧光标记的杂交探针对具有互补序列的PCR产物进行定量测定。因此并没有考虑利用PCR片段长度的多态性对目标RNA的鉴定的方式。茎环引物可以用于miRNA的反转录，也可应用于miRFLP分析法作为cDNA合成的引物。在传统茎环引物的PCR靶位点和探针之间引入不同碱基数目为编码，可以在同一反应中区分不同的miRNA目标，满足对多个miRNA的同时检测的需要。图9例举了利用碱基数编码的、针对miR-25和miR-92b的长度多态的茎环引物探针的设计方案。长度编码的茎环引物也可与不同长度的3’寡聚适配引物组合，增加miRNA的可测数目。图9中，引物A的PCR片段长度为67.17nt,引物B的PCR片段长度为70.08nt，可以在miRFLP分析谱上区分开来。 
表7 miRFLP分析法测定miR-25及miR-92b的茎环引物、动态miRNA标准（Std 1，2，3）及3’ 寡聚核苷酸适配引物（3’Adapter）的碱基序列及成分构成 
。