本发明属于中药成份检测分析领域，具体涉及一种鼓槌石斛的指纹图谱的分析建模方法。
背景技术：
：石斛是一种常用的滋补中药，主要为石斛属植物。石斛属植物是兰科植物中最大的一个属，包括铁皮石斛、金钗石斛、密花石斛、鼓槌石斛等多个品种。在全球范围内，约有1100多种石斛，其中在我国发现的有76种。石斛的药用历史悠久。早在古代，石斛在《神农本草经》中就被列为滋补上品，长久以来，随着时代的发展，石斛一直被人们视为珍贵的中草药，具有十分重要的滋补功效。在临床上，石斛被用于治疗多种疾病，具有增强免疫力、抗氧化、降血糖和抑制癌症等药理功效。由于人为长期无节制采挖及不合理利用石斛，其野生资源日趋减少，市场上出现了一些以假乱真、以次充好的现象。此外，由于石斛品种较多，其品种间的杂交使得其近缘的种存在性状交叉现象，分类区别比较困难。因此，有必要建立石斛的特征指纹图谱对石斛鲜品、干品的药材质量进行评价。鼓槌石斛(DendrobiumchrysotoxumLindley)：茎直立，肉质，纺锤形，长6cm～30cm，中部直径1.5cm～5cm，具2节～5节间，具多数圆钝的条棱，干后金黄色，近顶端具2枚～5枚叶；叶不下延为抱茎的鞘。总状花序从近茎顶端发出。斜出或稍下垂，长达20cm；花序轴粗壮，疏生多数花；花序柄基部具4枚～5枚鞘；花苞片小，膜质，卵状披针形，长0.2cm～0.3cm，先端急尖；花梗和子房黄色，长达5cm；花质地厚，金黄色，稍带香气；中萼片长圆形，长1.2cm～2cm，中部宽0.5cm～0.9cm，先端稍钝，具7条脉；侧萼片与中萼片近等大；萼囊近球形，宽约0.4cm；花瓣倒卵形，等长于中萼片，宽约为萼片的2倍，先端近圆形，具约10条脉；唇瓣颜色比萼片和花瓣深，近肾状圆形，长约2cm，宽2.3cm，先端浅2裂，基部两侧多少具红色条纹，边缘波状，上面密被短绒毛；唇瓣通常呈倒V型隆起，有时具U形的栗色斑块；蕊柱长约0.5cm；药帽淡黄色，尖塔状。色谱指纹图谱是一种综合的、可量化的鉴别手段，作为一种全局分析的全景模式，反映的是样品的整体情况，而对单一化学小分子成分的含量解释不具有针对性。迄今为止，有关药用石斛HPLC指纹图谱的专利有1项(CN200810218494.5)，成金乐等对石斛中的酚类和黄酮类成分进行了HPLC指纹图谱分析。此方法适合于酚类和黄酮类这两类成分的质量控制。另一方面，在色谱指纹图谱分析过程中，很多数据都是高维的，即数据包含很多属性或特征，比如有关鼓槌石斛色谱指纹图谱，就能更好地对鼓槌石斛进行描述，但在实际应用中对高维数据直接进行操作将会面临“维数灾难”的问题，“维数灾难”会导致建模过程所需要的样本数随着维数升高而呈指数级增长。面对高维数据，常规的最小二乘方法不再适用，为了提高模型的可解释性和预测的准确度，变量选择变得很重要。如何高效地从众多的变量中筛选出对自变量有重要作用的若干个变量，是在对指纹图谱进行分析时亟需解决的问题。另外，在现有的石斛检测标准中都是检测其中的多糖等大分子药效成分，但是，检测多糖等大分子成分需要的技术、仪器要求高，国内现有的仪器设备通常无法准确检测大分子成分，无法准确评价石斛的质量。技术实现要素：本发明建立了鼓槌石斛化学小分子成分的相关特征指纹图谱，通过该图谱实现对鼓槌石斛药材相关化学小分子成分含量的关联性分析。采用化学小分子成分的相关特征指纹图谱获取反映药用鼓槌石斛相关化学小分子成分含量的信息。本发明的目的是通过以下技术方案实现的：一种鼓槌石斛的化学小分子成分的相关特征指纹图谱的建立方法，所述相关特征指纹图谱的建立方法包括：1)对样本容量为n的样本分别进行色谱检测，获取以化学小分子成分夏佛托苷和/或柚皮素作为参照成分的检测数据，所述n取值为不小于10的整数；2)对样本分别进行指纹图谱检测，获取以鼓槌石斛的全化学成分的指纹图谱峰面积值；3)以色谱数据中夏佛托苷和/或柚皮素的成分含量值作为响应变量，将指纹图谱中其它成分的峰面积值作为自变量建立分析模型，进行主成分提取和回归计算，建立小分子成分与指纹图谱线性回归模型，并得到回归方程；4)通过变量相关系数指标(CorrelationCoefficient)表征自变量对回归建模的贡献程度，对应的相关系数越大的自变量的重要程度越大，获得和小分子成分密切相关的指纹图谱特征峰群，建立化学小分子成分的相关特征指纹图谱。在根据本发明的一个实施方案中，步骤1)所述的色谱检测还包括鼓槌石斛样品处理方法：所述样品处理方法包括：精密称取鼓槌石斛粉末，置100ml容量瓶中，每1g样品精密加体积比为75∶25的甲醇-水50ml，以250W功率、40kHz频率超声处理30分钟，冷却后过滤，将滤液旋蒸浓缩至干，每1g鼓槌石斛粉末对应地用5ml体积比为75∶25的甲醇-水溶解，上清液过0.45μm微孔滤膜，取续滤液，即得鼓槌石斛的小分子成分含量测定样品。在根据本发明的一个实施方案中，步骤1)所述的色谱检测的检测条件包括：含量测定色谱条件：GraceAllitimaC18色谱柱；流动相采用二元梯度洗脱系统，A相：0.2％醋酸-水，B相：乙腈；以波长290nm测定柚皮素，以波长334nm测定夏佛托苷，参比波长500nm，柱温30℃；流速1.0mL/min，进样量20μL；洗脱梯度为：指纹图谱色谱条件：GraceAllitimaC18色谱柱，优选为250mm×4.6mm，5μm规格的色谱柱；流动相：A相：0.4％乙酸+20mmol/L醋酸铵水溶液，B相：乙腈；梯度洗脱：0～12min：2％～15％B相，12～35min：15％～24％B相，35～45min：24％～36％B相，45～60min：36％～75％B相，60～80min：75％～95％B相；流速1.0mL/min；柱温30℃；进样量20μL；检测波长280nm。在根据本发明的一个实施方案中，步骤2)中所述的建立分析模型包括以下步骤：a)当色谱数据呈现超高维情形时，采用以下SIS(SureIndependenceScreening)方法通过式l将超高维降到d(d≤n)维；其中，d＝n或者Mγ＝{1≤i≤p：|ωi|是前|γn|个比较大的}式I其中，M*＝{1≤i≤p：βi≠0}表示真模型中非零系数的下标集；s＝|M*|表示非零系数的个数；ω＝(ω1，ω2，...，ωp)T＝XTy；当p＞n且对某个ξ有logp＝O(nξ)，M*与Mγ之间有以下联系：对于任意给定的γ∈(0，1)，ω的p个元素按绝对值从大到小排列并且定义；此时|γn|＜n，选取Mγ中下标对应的自变量，达到降维的目的；b)在SIS方法降维之后，采用SCAD(SmoothClippedAbsoluteDeviation)方法筛选变量。在根据本发明的一个实施方案中，所述SCAD方法是通过式II计算实现的：其中，惩罚函数pλ(β)的导函数p′λ(β)满足；a＞2，β＞0；优选地，d＝3.7。y＝(y1，y2，...，yn)T∈Rn为响应变量向量。以石斛数据为例，每种石斛有两个响应变量序列(夏佛托苷(μg/g)、柚皮素(μg/g))。响应变量受自变量的影响，一般情况下y为连续变量。x＝(x1，x2，...，xn)T为n×p的设计矩阵，包含对响应变量有影响的所有候选自变量。p为样本的维数，n为样本容量。在石斛数据中，维数p远大于样本容量n，因此最小二乘估计不再适用，需要采用变量选择的方法来进行模型估计。β＝(β1，β2，...，βp)T是一个p维的参数。为惩罚函数，λ为调整参数。在变量选择方法中，模型拟合的优良程度与对于入选变量个数惩罚的力度之间的平衡通过不同的准则来体现，而这里是通过直接选取调节参数来实现的，不同的λ值对应不同的惩罚力度。λ越大，压缩的程度越强，最后估计得到的非零参数越少，选择λ最常见的方法是K折交叉验证法：K-foldCV：一般，K可取为5或10。GCV准则是CV准则中的K取n时的一种近似情形，定义为：其中，SSEk是含有k个变量的CV子模型的残差平方和，df＝trace{P(λ)}。对于最终最优模型的选择，可以取CV值或GCV值最小的子模型。惩罚函数pλ(β)的导函数p′λ(β)满足其中参数a＞2，β＞0，可以通过二维广义交叉验证(GeneralizedCross-Validation)的方法来估计，但是计算比较复杂。通过贝叶斯的观点结合模拟分析，Fan提出a取3.7。SCAD方法满足无偏性、稀疏性和连续性，且在选取了合适的正则化参数后，SCAD方法得到的估计量满足Oracle性质。本发明还提供了上述方法在鼓槌石斛质量控制中的应用。在上述鼓槌石斛质量控制中所述鼓槌石斛可以为鲜品或干品。与现有技术相比，本发明至少具有以下优点：(a)本发明采用SCAD方法，对鼓槌石斛指纹图谱进行变量选择，有效的解决了“维数灾难”的问题；(b)本发明对原有指纹图谱进行降维，建立鼓槌石斛化学小分子成分的相关特征指纹图谱，对单一化学小分子成分的含量解释针对性和适用性更强；(c)本发明通过鼓槌石斛化学小分子成分的相关特征指纹图谱，实现化学小分子成分含量的关联性分析，能有效鉴别、控制和评价鼓槌石斛药材的质量。附图说明图1为鼓槌石斛指纹图谱色谱图；图2a为鼓槌石斛中夏佛托苷(1号峰)的检测图谱；图2b为鼓槌石斛中柚皮素(2号峰)的检测图谱。具体实施方式下面结合附图和实施例对本发明作进一步详述，以下实施例只是描述性的，不是限定性的，不能以此限定本发明的保护范围。实施例1鼓槌石斛的提取方法取鼓槌石斛干燥样品，用粉碎机粉碎，过药典筛(孔径0.335mm)，精密称取石斛粉末1.000g(称量误差不能超过0.2％)，置于100ml锥形瓶中，分别加入50mL75％甲醇(V水∶V甲醇＝25∶75)，室温下超声30min后取出，过滤，滤液旋蒸浓缩至干，用75％甲醇溶剂(V水∶V甲醇＝25∶75)溶解，最后转移至10ml容量瓶中定容，摇匀，以0.45μm微孔滤膜过滤，即得鼓槌石斛样品溶液。实施例2鼓槌石斛提取物的色谱检测方法1.鼓槌石斛的小分子成分含量测定样品提取处理方法：取本品粉末(过三号筛)1.00g，精密称定，置100ml容量瓶中，精密加甲醇-水(75∶25)50ml，超声处理(功率250W，频率40kHz)30分钟，放冷，过滤，滤液旋蒸浓缩至干，用5ml甲醇-水(75∶25)溶解，上清液过0.45μm微孔滤膜，取续滤液，即得。2.对照品溶液制备分别精密称定夏佛托苷4.10mg、柚皮素4.08mg，分别置于10ml容量瓶，加75％甲醇溶解稀释，摇匀，作为储备液。于4℃冰箱内冷藏备用。再分别精密吸取一定量的对照品储备溶液，加75％甲醇稀释同，分别准确配制夏佛托苷、柚皮素对照品溶液。通过不同的稀释比例，稀释制备成分7个浓度点。注入高效液相色谱仪。3.鼓槌石斛小分子成分含量测定色谱条件：色谱条件在根据本发明的一个实施方案中，步骤1)所述的色谱检测的检测条件包括：含量测定色谱条件：GraceAllitimaC18色谱柱；流动相采用二元梯度洗脱系统，A相：0.2％醋酸-水，B相：乙腈；以波长290nm测定柚皮素，以波长334nm测定夏佛托苷，参比波长500nm，柱温30℃；流速1.0mL/min，进样量20μL；洗脱梯度如表1所示。指纹图谱色谱条件：GraceAllitimaC18色谱柱，优选为250mm×4.6mm，5μm规格的色谱柱；流动相：A相：0.4％乙酸+20mmol/L醋酸铵水溶液，B相：乙腈；梯度洗脱：0～12min：2％～15％B相，12～35min：15％～24％B相，35～45min：24％～36％B相，45～60min：36％～75％B相，60～80min：75％～95％B相；流速1.0mL/min；柱温30℃；进样量20μL；检测波长280nm。表1洗脱梯度表图1所示为鼓槌石斛全成分指纹图谱色谱图，检测波长280nm；图2a为鼓槌石斛中夏佛托苷(1号峰)的检测图谱；图2b为鼓槌石斛中柚皮素(2号峰)的检测图谱。实施例3鼓槌石斛中与夏佛托苷的相关特征指纹图谱的建立1.石斛样品溶液制备取石斛干燥样品，用粉碎机粉碎，过药典筛(孔径0.335mm)，精密称取石斛粉末1.000g(称量误差不能超过0.2％)，置于100ml锥形瓶中，分别加入50mL75％甲醇(V水∶V甲醇＝25∶75)，室温下超声30min后取出，过滤，滤液旋蒸浓缩至干，用75％甲醇溶剂(V水∶V甲醇＝25∶75)溶解，最后转移至10ml容量瓶中定容，摇匀，以0.45μm微孔滤膜过滤，即得。表2为得到的对照品夏佛托苷的线性关系。表2对照品夏佛托苷线性关系表2.鼓槌石斛夏佛托苷的相关特征指纹图谱的建立方法采用线性模型，由于初始自变量p＝488个，样本量n＝13个，p＞＞n，故变量筛选是较为重要的工作。本研究采用SCAD方法筛选变量。SCAD方法筛选所得自变量有5个，R-square为0.7495。SCAD方法筛选所得自变量及其相应系数如表3所示，其中第一列为入选自变量编号，第二列为相应系数，第三列为系数方差，第四列为检验p值。表3：SCAD方法所得入选自变量及其相应系数以夏佛托苷为响应变量的筛选结果：SCAD方法的R方能达到0.74以上，入选自变量p值均小于0.05。上述自变量的含义：鼓槌石斛n批次样本(n不小于10)的指纹图谱根据保留时间对齐后的指纹图谱峰面积值，如X9即是样本容量中的n批次样本(n不小于10)在相同保留时间下的第9项对齐的指纹图谱峰面积值。综合来看，可采用5个自变量如表4所示(按相关系数程度从大到小排)。表45个筛选后的自变量X334X9X338X327X82实施例4鼓槌石斛柚皮素的相关特征指纹图谱的建立1.对照品溶液制备精密称定柚皮素4.08mg，置于10ml容量瓶，加75％甲醇溶解稀释，摇匀，作为储备液。于4℃冰箱内冷藏备用。再精密吸取一定量的对照品储备溶液，加75％甲醇稀释，准确配制柚皮素对照品溶液。通过不同的稀释比例，稀释制备成分7个浓度点，注入高效液相色谱仪。表5为得到的对照品的线性关系。表5对照品柚皮素线性关系表2.鼓槌石斛柚皮素相关特征指纹图谱的建立采用线性模型，由于初始自变量p＝488个，样本量n＝13个，p＞＞n，故变量筛选是较为重要的工作。本研究采用SCAD方法筛选变量。SCAD方法筛选所得自变量有3个，R-square为0.8063。SCAD方法筛选所得自变量及其相应系数如表6所示，其中第一列为入选自变量编号，第二列为相应系数，第三列为系数方差，第四列为检验P值。表6：SCAD方法所得入选自变量及其相应系数以柚皮素为响应变量结果：SCAD方法的R方能达到0.80以上，入选自变量p值均小于0.05。综合来看，可采用3个自变量如表7所示(按相关系数程度从大到小排)。表7：3个筛选后的自变量X210X105X277上述自变量的含义：齿瓣石斛的n批次样本(n不小于10)的指纹图谱根据保留时间对齐后的指纹图谱峰面积值，如X210即是n批次样本在相同保留时间下的第210项对齐的指纹图谱峰面积值。实施例5鼓槌石斛化学小分子成分的相关特征指纹图谱在鼓槌石斛质量评价中的应用1.供试品测定样品的制备精密称取鼓槌石斛粉末，置100ml容量瓶中，每1g样品精密加体积比为75∶25的甲醇-水50ml，以250W功率、40kHz频率超声处理30分钟，冷却后过滤，将滤液旋蒸浓缩至干，每1g鼓槌石斛粉末对应地用5ml体积比为75∶25的甲醇-水溶解，上清液过0.45μm微孔滤膜，取续滤液，即得鼓槌石斛的小分子成分含量测定样品；2.色谱检测色谱条件：色谱柱：GraceAllitimaC18色谱柱(250mm×4.6mm，5μm)；流动相：A相：0.4％乙酸+20mmol/L醋酸铵水溶液，B相：乙腈；梯度洗脱：0～12min：2％～15％B相，12～35min：15％～24％B相，35～45min：24％～36％B相，45～60min：36％～75％B相，60～80min：75％～95％B相；流速，1.0mL/min；柱温30℃；进样量20μL；检测波长280nm。样品制备方法：取本品粉末(过三号筛)1.00g，精密称定，置100ml容量瓶中，精密加甲醇-水(75∶25)50ml，超声处理(功率250W，频率40kHz)30分钟，放冷，过滤，滤液旋蒸浓缩至干，用5ml甲醇-水(75∶25)溶解，上清液过0.45μm微孔滤膜，取续滤液，即得。在检测时以波长280nm测定全成分指纹图谱，将得到的全成分指纹图谱与图1为对照的指纹色谱图进行相似度比对；相似度大于0.85为质量合格。以上，仅为本发明较佳的具体实施方式，但发明的保护范围并不局限于此，任何熟悉本
技术领域：
的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。当前第1页1 2 3