CN108398516A

CN108398516A - 一种流苏石斛的质量检测方法

Info

Publication number: CN108398516A
Application number: CN201710064173.3A
Authority: CN
Inventors: 赵田; 张国强
Original assignee: Beijing Lanbiao Yicheng Technology Co Ltd
Current assignee: BEIJING LANBIAO YICHENG TECHNOLOGY Co.,Ltd.
Priority date: 2017-02-04
Filing date: 2017-02-04
Publication date: 2018-08-14
Anticipated expiration: 2037-02-04
Also published as: CN108398516B

Abstract

本发明公开流苏石斛的质量检测方法，包括：1)以ITS‑26SE和ITS‑17SE为引物，进行测序，以鉴定待测石斛药材的品种；2)对样本容量为n的样本进行色谱检测，获取以化学小分子成分夏佛托苷和/或柚皮素作为参照成分的检测数据；3)对样本分别进行指纹图谱检测，获取流苏石斛的全化学成分的指纹图谱峰面积值；4)以色谱数据中夏佛托苷和/或柚皮素的成分含量值作为响应变量，将指纹图谱中的其它成分的峰面积值作为自变量建立分析模型，通过Lasso方法筛选变量建立化学小分子成分的相关特征指纹图谱模型。通过一代测序及特征指纹图谱，精确地鉴别和控制流苏石斛药材的质量。

Description

一种流苏石斛的质量检测方法

技术领域

本发明属于中药成份检测分析领域，具体涉及一种流苏石斛的质量检测方法。

背景技术

流苏石斛(学名：Dendrobium fimbriatum Hook.)：茎粗壮，斜立或下垂，质地硬，圆柱形或有时基部上方稍呈纺锤形，长50cm～85cm，直径0.8cm～1.5cm，不分枝，具多节，干后淡黄色或淡黄褐色，节间长3.5cm～4.8cm，具多数纵槽；叶基部具紧抱于茎的革质鞘。总状花序长5cm～15cm，疏生6朵～12朵花；花序轴较细，多少弯曲；花序柄长2cm～4cm，基部被数枚套叠的鞘；鞘革质，筒状，位于基部的最短，长约0.3cm，顶端的最长，达1cm；花苞片膜质，卵状三角形，长0.3cm～0.5cm，先端锐尖；花梗和子房浅绿色，长2.5cm～3cm；花金黄色，质地薄，开展，稍具香气；中萼片长圆形，长1.3cm～1.8cm，宽0.6cm～0.8cm，先端钝，边缘全缘，具5条脉；侧萼片卵状披针形，与中萼片等长而稍狭窄，先端钝，基部歪斜，全缘，具5条脉；萼囊近圆形，长约0.3cm；花瓣长圆状椭圆形，长1.2cm～1.9cm，宽0.7cm～1.0cm，先端钝，边缘微啮蚀状，具5条脉；唇瓣比萼片和花瓣的颜色深，近圆形，长1.5cm～2cm，基部两侧具紫红色条纹并且收狭为长约0.3cm的爪，边缘具复流苏，唇瓣具一个新月形横生的深紫色斑块，上面密布短绒毛；蕊柱黄色，长约0.2cm，具长约0.4cm的蕊柱足；药帽黄色，圆锥形，光滑，前端边缘具细齿。流苏石斛的茎可入药，味甘，微寒。益胃生津，滋阴清热，可用于阴伤津亏，口干烦渴，食少干呕，病后虚热，目暗不明。外用治跌打损伤，骨折伤筋；内服治咽喉痒，咳嗽。

石斛一直被人们视为珍贵的中草药，具有十分重要的滋补功效。在临床上，石斛被用于治疗多种疾病，具有增强免疫力、抗氧化、降血糖和抑制癌症等药理功效。由于人为长期无节制采挖及不合理利用石斛，其野生资源日趋减少，市场上出现了一些以假乱真、以次充好的现象。此外，由于石斛品种较多，其品种间的杂交使得其近缘的种存在性状交叉现象，分类区别比较困难。因此，有必要建立石斛的特征指纹图谱对石斛的药材质量进行评价。

色谱指纹图谱是一种综合的、可量化的鉴别手段，作为一种全局分析的全景模式，反映的是样品的整体情况。但在色谱指纹图谱分析过程中，很多数据都是高维的，即数据包含很多属性或特征，比如有关流苏石斛色谱指纹图谱，就能更好地对流苏石斛进行描述，但在实际应用中对高维数据直接进行操作将会面临“维数灾难”的问题，“维数灾难”会导致建模过程所需要的样本数随着维数升高而呈指数级增长。面对高维数据，常规的最小二乘方法不再适用，为了提高模型的可解释性和预测的准确度，变量选择变得很重要。如何高效地从众多的变量中筛选出对因变量有重要作用的若干个变量，是在对指纹图谱进行分析时亟需解决的问题。

目前国家药典采用定量测定中药材中某一活性成分或有效成分即小分子成分的含量的高低来评价其质量。但研究证明，中药的疗效是来自其多种“活性成分”之间的协同作用，甚至是被普遍公认的有效的“活性成分”与“非活性成分”之间的协同作用或“生克作用”才能达到中药的疗效，而不是某一活性成分单独作用的结果。在中医理论指导下的中药，任何一种活性成分均不能全面反映中医用药所体现的整体疗效。

发明内容

本发明提供了一种流苏石斛的质量检测方法，通过一代测序鉴定石斛药材的品种，利用Lasso方法筛选变量建立流苏石斛中的化学小分子成分的相关特征指纹图谱模型，通过一代测序和相关特征指纹图谱模型准确评价流苏石斛的药材的质量。

本发明的目的是通过以下技术方案实现的：

一种流苏石斛的质量检测方法，包括：

1)以ITS-26SE：5’GAATTCCCCGGTTCGCTCGCCGTTAC 3’；

ITS-17SE：5’ACGAATTCATGGTCCGGTGAAGTGTTCG 3’为引物，进行PCR扩增测序，以鉴定待测石斛药材的品种为流苏石斛样本；

2)对样本容量为n的流苏石斛样本进行色谱检测，获取以化学小分子成分夏佛托苷和/或柚皮素作为参照成分的检测数据；

3)对样本分别进行指纹图谱检测，获取流苏石斛的全化学成分的指纹图谱峰面积值；

4)以色谱数据中夏佛托苷和/或柚皮素的成分含量值作为响应变量，将指纹图谱中的其它成分的峰面积值作为自变量建立分析模型，通过Lasso(The Least AbsoluteShrinkage and Selection Operator)方法筛选变量建立化学小分子成分的相关特征指纹图谱模型，其基本模型为：

y＝X^Tβ+ε

其中，y为响应变量，y＝(y₁,y₂,...,y_n)^T；X为矩阵，X＝(x₁,x₂,...,x_n)^T；E(ε)＝0；Var(ε)＝σ²I_n；ε为模型的随机误差项；σ是随机误差项的标准差；n为样本量；I_n是一个n×n的单位阵。

假定随机项服从古典假定，即：

(1)随机项具有零均值，E(ε_i|x_i)＝0；

(2)随机项具有同方差，Var(ε_i|x_i)＝σ²；

(3)随机项无序列相关性，Cov(ε_i,ε_j)＝0,i≠j；

(4)ε服从正态分布，ε_i～N(0,σ²)。

随机项的方差矩阵是一个对角线为σ²，其他地方为0的方阵，如下所示：

其中，I_n是一个n×n的单位阵，n为数据的样本量，

进一步地，所述Lasso方法是通过式Ⅰ计算实现的：

在式Ⅰ中，n为样本量；p^*为变量数；p为样本的维数；y＝(y₁,y₂,...,y_n)^T∈Rⁿ为响应变量；x＝(x₁,x₂,...,x_n)^T为n×p的设计矩阵，包含对响应变量有影响的所有候选自变量；λ为调整参数；为惩罚函数；β₀的含义为公式的截距项，也就是当所有自变量x为0时响应变量y的值；β_j的含义是自变量x_j的系数，即自变量x_j对响应变量y的影响程度。

进一步地，所述λ的选择方法为K折交叉验证法：

K-fold CV：

其中，K为5或10。

进一步地，所述分析模型为取CV值最小的子模型。

进一步地，所述λ的选择遵循GCV准则，所述GCV准则定义为：

其中，SSE_k是含有k个变量的CV子模型的残差平方和，df＝trace{P(λ)}；trace表示矩阵的迹。在线性代数中，一个n×n的矩阵A的主对角线(从左上方至右下方的对角线)上各个元素的总和被称为矩阵A的迹(或迹数)，一般记作tr(A)。也就是说，df等于矩阵P(λ)中主对角线上所有元素的和。

进一步地，所述分析模型为取GCV值最小的子模型。

进一步地，当色谱数据呈现超高维情形时，首先采用以下SIS(SureIndependenceScreening)方法筛选变量，再利用Lasso方法处理；

SIS：M_γ＝{1≤i≤p:|ω_i|是前|γn|个比较大的}

其中，M*＝{1≤i≤p:β_i≠0}表示真模型中非零系数的下标集；s＝|M_*|表示非零系数的个数；ω＝(ω₁,ω₂,...,ω_p)^T＝X^Ty；对于任意给定的γ∈(0,1)，ω的p个元素按绝对值从大到小排列并且定义；此时|γn|＜n，选取M_γ中下标对应的自变量，是超高维降到d(d≤n)维；其中，d＝n或者d＝[n/log n]。

本发明还提供了上述方法在流苏石斛质量控制中的应用。

与现有技术相比，本发明至少具有以下优点：

(a)本发明中先以ITS-26SE和ITS-17SE为引物测定石斛药材的特征序列，以确定该石斛药材为流苏石斛品种；其后以夏佛托苷和/或柚皮素的色谱数据作为自变量，将上述检测数据中的其它成分的色谱数据作为因变量建立小分子成分与指纹图谱线性回归模型，使石斛药材的质量评价更为精确；

(b)本发明采用Lasso方法对流苏石斛指纹图谱进行变量选择，有效的解决了“维数灾难”的问题；

(c)本发明对原有指纹图谱进行降维，建立流苏石斛化学小分子成分的相关特征指纹图谱，对单一化学小分子成分的含量解释针对性和适用性更强；

(d)本发明通过流苏石斛化学小分子成分的相关特征指纹图谱，实现化学小分子成分含量的关联性分析，能有效鉴别和控制流苏石斛药材的质量；

(e)当色谱数据呈现超高维情形时，首先采用SIS方法进行降维，再利用Lasso方法处理。

附图说明

图1为流苏石斛全成分指纹图谱色谱图；

图2为夏佛托苷化学小分子物质标准样品的色谱图；

图3为柚皮素化学小分子物质标准样品的色谱图；

图4为流苏石斛中夏佛托苷化学小分子的色谱图；

图5为流苏石斛中柚皮素化学小分子的色谱图。

注：1号峰夏佛托苷；2号峰柚皮素。

具体实施方式

下面结合附图和实施例对本发明作进一步详述，以下实施例只是描述性的，不是限定性的，不能以此限定本发明的保护范围。

实施例1 流苏石斛的一代测序

一代测序引物序列：

ITS-26SE：5’GAATTCCCCGGTTCGCTCGCCGTTAC 3’；

ITS-17SE：5’ACGAATTCATGGTCCGGTGAAGTGTTCG 3’。

扩增测序参数为：98℃变性2min后进行PCR循环，PCR循环参数为98℃20s；52℃30s；68℃1min，38个循环，68℃7min，扩增结束后设置4℃保温，并进行一代分子测序。

通过一代测序，鉴别待测石斛药材的品种为流苏石斛。

实施例2 流苏石斛的提取方法

取流苏石斛干燥样品，用粉碎机粉碎，过药典筛(孔径0.335mm)，精密称取石斛粉末1.000g(称量误差不能超过0.2％)，置于100ml锥形瓶中，分别加入50mL 75％甲醇(V水：V甲醇＝25:75)，室温下超声30min后取出，过滤，滤液旋蒸浓缩至干，用75％甲醇溶剂(V水：V甲醇＝25:75)溶解，最后转移至10ml容量瓶中定容，摇匀，以0.45μm微孔滤膜过滤，即得流苏石斛样品溶液。

实施例3 流苏石斛提取物的色谱检测方法

①对照品溶液制备

分别精密称取夏佛托苷4.10mg和柚皮素4.08mg，分别置于10ml容量瓶中，加75％(V/V)甲醇溶解稀释，摇匀，作为储备液。于4℃冰箱内冷藏备用。

再分别精密吸取一定量的对照品储备溶液，利用75％甲醇稀释，准确配制夏佛托苷和柚皮素混合对照品溶液。通过不同的稀释比例，稀释制备成分7个浓度点。注入高效液相色谱仪。

②流苏石斛的小分子成分含量测定样品提取处理方法：

取本品粉末(过三号筛)1.00g，精密称定，置100ml容量瓶中，精密加甲醇-水(75:25)50ml，超声处理(功率250W，频率40kHz)30分钟，放冷，过滤，滤液旋蒸浓缩至干，用5ml甲醇-水(75:25)溶解，上清液过0.45μm微孔滤膜，取续滤液，即得。

③流苏石斛小分子成分含量测定色谱条件：

色谱条件：

含量测定色谱条件：Grace Allitima C18色谱柱(250mm*4.6mm，5μm)；流动相采用二元梯度洗脱系统，A相：0.2％醋酸-水，B相：乙腈；梯度洗脱程序如表1；以波长290nm测定柚皮素，以波长334nm测定夏佛托苷；参比波长为500nm，柱温30℃；流速1.0mL/min，进样量20μL。

指纹图谱色谱条件：Grace Allitima C18色谱柱，优选为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；

图2为夏佛托苷(1号峰)标准样品的检测图谱图；

图3为柚皮素(2号峰)标准样品的检测图谱图；

图4为流苏石斛中夏佛托苷(1号峰)的检测图谱；

图5为流苏石斛中柚皮素(2号峰)的检测图谱。

实施例4 流苏石斛夏佛托苷的相关特征指纹图谱的建立

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.流苏石斛夏佛托苷的相关特征指纹图谱的建立方法

第一步：计算所有协变量x与y的相关系数；

第二步：将相关系数的绝对值从大到小进行排列，选出前2√n个协变量，记为x_1,x_2,…,x_p；

第三步：将y与x_1,x_2,…,x_p进行线性回归，采用Lasso方法，进行变量选择。

Lasso(Least Absolute Shrinkage and Selection Operator)函数的第一部分表示模型拟合的优良性，第二部分可以视为惩罚。该方法把小的系数往0压缩，一旦某个系数被压缩到0，对应的变量就被删除。就好像用“筛子”过滤，把影响小的变量一次就筛掉了。λ越小，模型中的变量越多λ越大，收缩量越大，选出的变量就越少。而Lasso方法是一种连续的、有序的过程，方差较小。当调节参数足够大时，惩罚项具有将其中某些系数的估计值强制设定为0的作用，因而Lasso方法可以进行变量选择，能够得到稀疏模型。

当自变量为p，样本量为n，当p>>n，首先采用SIS方法降维，再采用Lasso方法筛选变量。

SIS：M_γ＝{1≤i≤p:|ω_i|是前|γn|个比较大的}

其中，M_*＝{1≤i≤p:β_i≠0}表示真模型中非零系数的下标集；s＝|M_*|表示非零系数的个数；ω＝(ω₁,ω₂,...,ω_p)^T＝X^Ty；对于任意给定的γ∈(0,1)，ω的p个元素按绝对值从大到小排列并且定义；此时|γn|＜n，选取M_γ中下标对应的自变量，使超高维降到d(d≤n)维；其中，d＝n或者d＝[n/logn]。

Lasso方法筛选的线性模型为：

其中，y_i为第i个响应变量，y＝(y₁,y₂,...,y_n)；X_i是P_n×1阶的协变量，X＝(x₁,x₂,...,x_n)'；ε_i是均值为0，方差为σ²的i.i.d的随机误差项，E(ε)＝0，Var(ε)＝σ²I_n。

假定随机项服从古典假定，即：

(1)随机项具有零均值，E(ε_i|x_i)＝0；

(2)随机项具有同方差，Var(ε_i|x_i)＝σ²；

(3)随机项无序列相关性，Cov(ε_i,ε_j)＝0,i≠j；

(4)ε服从正态分布，ε_i～N(0,σ²)。

其中，I_n是一个n×n的单位阵，n为数据的样本量，

为了同时进行变量选择和对参数进行估计，Lasso方法通过惩罚最小二乘目标函数式Ⅰ的最小化来实现。

其中，y＝(y₁,y₂,...,y_n)^T∈Rⁿ为响应变量向量。以石斛数据为例，每种石斛有两个响应变量序列(夏佛托苷(μg/g)、柚皮素(μg/g))。响应变量受自变量的影响，一般情况下y为连续变量。

x＝(x₁,x₂,...,x_n)^T为n×p的设计矩阵，包含对响应变量有影响的所有候选自变量。

p为样本的维数，n为样本容量。在石斛数据中，维数p远大于样本容量n，因此最小二乘估计不再适用，需要采用变量选择的方法来进行模型估计。

β＝(β₁,β₂,...,β_p)^T是一个p维的参数。

为惩罚函数，λ为调整参数。在变量选择方法中，模型拟合的优良程度与对于入选变量个数惩罚的力度之间的平衡通过不同的准则来体现，而这里是通过直接选取调节参数来实现的，不同的λ值对应不同的惩罚力度。λ越大，压缩的程度越强，最后估计得到的非零参数越少，选择λ最常见的方法是K折交叉验证法：

K-fold CV：

一般，K可取为5或10。

GCV准则是CV准则中的K取n时的一种近似情形，定义为：

其中，SSE_k是含有k个变量的CV子模型的残差平方和，df＝trace{P(λ)}。对于最终最优模型的选择，可以取CV值或GCV值最小的子模型。

采用线性模型，由于初始自变量p＝431个，样本量n＝6个，p>>n，故变量筛选是较为重要的工作。先经SIS降维后，采用Lasso方法筛选变量。

Lasso方法筛选所得自变量有2个，R-square为0.9972。

Lasso方法筛选所得自变量及其相应系数如下表3，其中第一列为入选自变量编号，第二列为相应系数，第三列为系数方差，第四列为检验P值。

表3 Lasso方法所得入选自变量及其相应系数

x_model(选择的自变量X)	β(自变量X的系数)	Varbeta(β的方差)	p-value(P值)
				常数项	1.643352	0.44608	0.00023
342	0.040709	0.00515	2.69E-15
				72	0.003219	0.000778	3.5E-05

表3的结果显示，以夏佛托苷为响应变量的筛选结果：第1列为采用Lasso方法选出的变量，即选出了变量342和72，其所对应的p值(列4)均小于显著性水平0.05，具有显著性差异。上述自变量的含义：流苏石斛n批次样本(n不小于10)的指纹图谱根据保留时间对齐后的指纹图谱峰面积值。

列2为每个变量具体对应的β参数值。β值为正说明该变量对石斛夏佛托苷小分子存在正向影响；β值为负说明该变量对石斛夏佛托苷小分子存在负向影响。β值的绝对值大小显示的是该变量对石斛夏佛托苷小分子影响程度的大小。具体来说，在表3中，变量342、72对石斛夏佛托苷小分子的影响均为正，其中变量342的正向影响较大。

实施例5 流苏石斛柚皮素的相关特征指纹图谱的建立

1.对照品溶液制备

精密称定柚皮素4.08mg，置于10ml容量瓶，加75％甲醇溶解稀释，摇匀，作为储备液。于4℃冰箱内冷藏备用。再精密吸取一定量的对照品储备溶液，加75％甲醇稀释，准确配制柚皮素对照品溶液。通过不同的稀释比例，稀释制备成分7个浓度点，注入高效液相色谱仪。表4为得到的对照品柚皮素的线性关系。

表4 对照品柚皮素线性关系表

2.柚皮素相关特征指纹图谱的建立

采用线性模型，由于初始自变量p＝431个，样本量n＝8个，p>>n，故变量筛选是较为重要的工作。本发明采用Lasso方法筛选变量。

Lasso方法筛选所得自变量有3个，R-square为0.7414。

Lasso方法筛选所得自变量及其相应系数如下表5，其中第一列为入选自变量编号，第二列为相应系数，第三列为系数方差，第四列为检验P值。

表5 Lasso方法所得入选自变量及其相应系数

x_model(选择的自变量X)	β(自变量X的系数)	Varbeta(β的方差)	p-value(P值)
				常数项	7.347645	0.592511	2.59E-35
184	-0.00186	0.001227	0.129436
				89	0.003206	0.00139	0.021133
217	-0.01622	0.007275	0.025758

表5的结果显示，以柚皮素为自变量结果：采用Lasso方法选出的变量即为列1所示，即筛选出了变量184、89、217，该变量的p值(列4)均小于显著性水平0.05，具有显著性差异。上述自变量的含义：流苏石斛的n批次样本(n不小于10)的指纹图谱根据保留时间对齐后的指纹图谱峰面积值。

列2给出的是每个变量具体对应的β参数值。β值为正说明该变量对石斛柚皮素小分子存在正向影响；β值为负说明该变量对石斛柚皮素小分子存在负向影响。β值的绝对值大小显示的是该变量对石斛柚皮素小分子影响程度的大小。具体来说，在表5中，变量89对石斛柚皮素小分子的影响为正；变量184、217对石斛柚皮素小分子的影响为负，其中变量217对石斛柚皮素小分子的影响较大。

根据一代分子数据的结果显示该该技术可清晰区地把流苏石斛与其他石斛区分开来，具有可鉴定种的作用，一代测序分子序列可以作为该种鉴定指标。同时，流苏石斛的特征指纹图谱的两种化学小分子含量在种内各样品表现稳定，并能明确与其他石斛区分开来，也可以作为该种的鉴定指标之一。因此，一代数据和指纹图谱均可作为鉴定该种的指标，并且两者有以下关联关系：(1)当一代数据鉴定样品为流苏石斛时，该样品的指纹图谱具有特定的特征，即由一代数据即可鉴定出该种，并可知该种的指纹图谱特征；(2)当指纹图谱鉴定样品为流苏石斛时，也可以推知其一代数据。因此两者互为充分必要条件，并且由于该种的指纹图谱特征在各代表样品中表现稳定，可确定该种的药效成份含量。因此，一代测序数据及指纹图谱均可作为该种的鉴定指标，并可作为药材的评价指标。若在品种鉴定与质量评价中分别采用一代测序和指纹图谱进行测定时，可将两种方法的鉴定结果结合在一起，使得对石斛药材的鉴定结果更加的准确、可靠。

实施例7 流苏石斛的一代测序及小分子成分的相关特征指纹图谱在石斛药材质量评价中的应用

1.利用以下一代测序引物序列及扩增测序参数测定待测石斛药材的序列，将其鉴定为流苏石斛药材。

一代测序引物序列：

ITS-26SE：5’GAATTCCCCGGTTCGCTCGCCGTTAC 3’；

ITS-17SE：5’ACGAATTCATGGTCCGGTGAAGTGTTCG 3’；

扩增测序参数为：98℃变性2min后进行PCR循环，PCR循环参数为98℃20s；52℃30s；68℃1min，38个循环，68℃7min，扩增结束后设置4℃保温。

通过以上测序，确定待测石斛药材为流苏石斛品种。

2.供试品测定样品的制备

精密称取流苏石斛粉末，置100ml容量瓶中，每1g样品精密加体积比为75:25的甲醇-水50ml，以250W功率、40kHz频率超声处理30分钟，冷却后过滤，将滤液旋蒸浓缩至干，每1g流苏石斛粉末对应地用5ml体积比为75:25的甲醇-水溶解，上清液过0.45μm微孔滤膜，取续滤液，即得流苏石斛的小分子成分含量测定样品。

3.色谱检测

色谱条件：

色谱柱：Grace Allitima C18色谱柱(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。

样品制备方法：

在检测时以波长280nm测定全成分指纹色谱图，将得到的全成分指纹图谱与图1为对照的指纹图谱进行相似度比对；相似度大于0.85为质量合格。

中药指纹图谱是分析仪器检测得到的各种反映中药材、半成品和中成药(或植物药)所含复杂化学物质成分分布的量化特征关联药效活性控制为特点,从宏观上整体反映中药材、半成品和中成药(或植物药)中所含化学物质成分的种类、数量和含量特征，并能量化揭示潜在复杂的生物活性信息特征的图谱。

本发明首先通过一代测序确定石斛的种类，再通过测定石斛小分子成分含量数据与石斛指纹图谱全部峰面积进行关联性研究，通过关联数据建模，找出石斛内在的关联性，这样即可全面评价石斛质量，且能有效精确地鉴别和控制流苏石斛药材的质量，使得分析结果更加可靠，避免其他种类的石斛药材的干扰。

以上，仅为本发明较佳的具体实施方式，但发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。

Claims

1.一种流苏石斛的质量检测方法，其特征在于，包括：

1)以ITS-26SE：5’GAATTCCCCGGTTCGCTCGCCGTTAC 3’和

ITS-17SE：5’ACGAATTCATGGTCCGGTGAAGTGTTCG 3’为引物，进行PCR扩增测序，以鉴定待测石斛药材的品种为流苏石斛；

2)对样本容量为n的样本进行色谱检测，获取以化学小分子成分夏佛托苷和/或柚皮素作为参照成分的检测数据；

y＝X^Tβ+ε

2.根据权利要求1所述的方法，其特征在于，所述Lasso方法是通过式Ⅰ计算实现的：

3.根据权利要求2所述的方法，其特征在于，所述λ的选择方法为K折交叉验证法：

K-fold CV：

其中，K为5或10。

4.根据权利要求3所述的方法，其特征在于，所述分析模型为取CV值最小的子模型。

5.根据权利要求2所述的方法，其特征在于，所述λ的选择遵循GCV准则，所述GCV准则定义为：

其中，SSE_k是含有k个变量的CV子模型的残差平方和，df＝trace{P(λ)}。

6.根据权利要求5所述的方法，其特征在于，所述分析模型为取GCV值最小的子模型。

7.根据权利要求2所述的方法，当色谱数据呈现超高维情形时，首先采用以下SIS(SureIndependence Screening)方法筛选变量，再利用Lasso方法处理；

SIS：M_γ＝{1≤i≤p:|ω_i|是前|γn|个比较大的}

其中，M*＝{1≤i≤p:β_i≠0}表示真模型中非零系数的下标集；s＝|M_*|表示非零系数的个数；ω＝(ω₁,ω₂,...,ω_p)^T＝X^Ty；对于任意给定的γ∈(0,1)，ω的p个元素按绝对值从大到小排列并且定义；此时|γn|＜n，选取M_γ中下标对应的自变量，使超高维降到d(d≤n)维；其中，d＝n或者d＝[n/log n]。