基于ITS2 序列及二級結構對紫蘇與其變種、紫蘇子與其混偽品的鑒別

鄭夢迪， 賀紫涵， 張 春， 汪興軍

（西安醫學院 藥學院，陜西 西安 710021）

紫蘇Perilla frutescens（L.） Britt.為唇形科紫蘇屬植物，是一種常見的藥食兩用植物資源[1]，據2015版《中國藥典》記載，紫蘇干燥成熟的果實入藥稱為紫蘇子，具有降氣化痰、止咳平喘、潤腸通便的作用；其干燥葉（或帶嫩枝）為紫蘇葉入藥，具有解表散寒、行氣和胃的功效[2]。 由于紫蘇屬下變種較多，歷代本草文獻中對于紫蘇及其變種用途及分類的記載混亂，其植物命名和分類爭議存在已久[3-4]。 魏晉時期《名醫別錄》最早記載紫蘇入藥，分為“蘇”和“荏”，“蘇”即現代所使用的紫蘇，“荏”即指白蘇[3]。南宋王介編繪的《履巉巖本草》中記載了野蘇，經考證為紫蘇的近似種。 《本草綱目》中除對紫蘇和白蘇的記載， 另外還收錄了一種花紫蘇，“其葉細齒，紐如剪成之狀， 色香莖子并無異者， 人稱回回蘇”[3]。《中國植物志》對紫蘇屬劃分為紫蘇P. frutescens（原變種）、白蘇P. frutescens var. frutescens、回回蘇P. frutescens var. crispa（Thunb.） Hand.-Mazz.、耳齒紫蘇P. frutescens var. auriculato-dentata C. Y. Wu.et Hsuan 和野生紫蘇P. frutescens var. purpurascens（Hayata） H. W. Li[1]。然而由于紫蘇屬植物在我國廣泛人工栽培，屬下不同變種植物的表型特征非常相似，并且表型往往受到植物自身生長周期和外界環境因素的影響， 從而造成傳統鑒定方法易產生混淆，進而危及用藥安全。 另一方面，紫蘇子在加工食品、功能性飲品、保健品及化妝品等領域的應用逐漸增多， 但產量低不能滿足增長的市場需求量；種苗錯誤種植、混合種植、混合收購、人為摻入等現象屢見不鮮，再加之藥用紫蘇子經過炮制后外觀特征丟失，常有急性子、沙苑子、白芥子、油菜子和菟絲子等種子混入， 從而造成傳統方法對紫蘇子與其混偽品的鑒別困難[5-6]。

DNA 條形碼（DNA barcoding）技術是利用生物基因組上一段具有足夠變異， 且易于PCR 擴增的片段，可用來區分不同物種的分子生物學技術，是對傳統物種鑒定方法的有效補充[7]。 其中ITS2（Internal transcribed spacer 2）條形碼在植物分類鑒別中的穩定性和適用性已得到廣泛驗證[8-11]。 ITS2片段位于5.8S 和28S 核rDNA 基因之間。 rDNA 轉錄形成rRNA，RNA 的一級結構是指A、U、C、G 4 種核苷酸的排列組合；RNA 的二級結構是指其自身回折形成莖環結構，不配對堿基構成環區，莖干區由配對堿基構成[12-13]；而在莖干區，有時可發生補償性堿基替換（Compensatory base changes，簡稱CBC），即為了維持配對結構，一側堿基發生突變會引起另一側堿基的響應配對性突變，如C-G 配對中，胞嘧啶C 突變為腺嘌呤A， 原配對中的G 相應突變為U，即由C-G 變化為A-U，從而維持二級結構的穩定性。 研究表明，發生CBC 的兩個個體之間存在生殖隔離，CBC 與雜交不親和有很強的相關性，即ITS2 二級結構中的CBC 是一種對真核生物進行分類的有效方法[14-15]。

近年來對紫蘇與其變種的分子鑒別僅停留在核酸序列層面[16]，并且尚無對紫蘇子與其混偽品的DNA 條形碼鑒別研究。 因此作者擬采用ITS2 序列結合二級結構，對紫蘇與其變種、紫蘇子與其混偽品進行分析，為中藥材紫蘇葉與紫蘇子在分子水平上的鑒別提供科學依據，保障臨床用藥的合理性和有效性。

1 材料與方法

1.1 材料

本研究收集紫蘇植物材料3 份，分別來自浙江衢州、陜西西安、桂林全州；白蘇植物1 份，來源于江西贛州； 于GenBank 獲得ITS 序列， 包括紫蘇3條、白蘇10 條、野生紫蘇2 條、回回蘇3 條、耳齒紫蘇3 條、3 個日本紫蘇品種 （文中記為日本紫蘇I、II、III）各1 條、外類群枸骨5 條；共33 份樣本信息見表1。 另外，紫蘇子與其混偽品樣本信息見表2，包括紫蘇子基源植物樣本2 份，分別來源于浙江衢州和陜西西安； 菟絲子和油菜子基源植物樣本各1份， 分別來自廣東茂名和陜西西安；GenBank 下載ITS 序列包括紫蘇子3 條、 急性子4 條、 白芥子5條、沙苑子3 條、菟絲子3 條、油菜子2 條；共24 份樣本。 所有植物材料經西安醫學院藥學院生藥教研室汪興軍老師鑒定，并保存于-80 ℃冰箱。

表1 紫蘇與其變種的材料種類和來源的詳細資料Table 1 Detailed information of materials types and sources of P. frutescens and its varieties

表2 紫蘇子與其混偽品的材料種類和來源的詳細資料Table 2 Samples of Perillae fructus and its adulterants

1.2 方法

1.2.1 DNA 提取 用體積分數75%乙醇擦拭新鮮植物樣本葉或莖的表面， 取約150 mg 組織于1.5 mL 滅菌離心管中，迅速放入液氮中冷凍。 在液氮環境下，用滅菌小木棒研磨植物組織至粉末狀。 再按照DP305 植物基因組DNA 提取試劑盒說明書操作提取DNA。

1.2.2 PCR 擴增、 測序和校正 選擇ITS2 通用引物進行序列擴增[17]，ITS2F（5′-ATGCGATACTTGGT GTGAAT-3′）和ITS2R（5′-GACGCTTCTCCAGACT ACAAT-3′）。RCR 反應條件：94 ℃5 min，94 ℃30 s、56 ℃30 s、72 ℃45 s，循環40 次，72 ℃2 min。PCR反應體系：DNA 模板1 μL， 引物ITS2F 和ITS2R 各1 μL，DNA 聚合酶2×Taqmix 12.5 μL，ddH2O 補充至25 μL。 PCR 產物由北京擎科澤西生物科技有限公司完成Sanger 測序，每個樣本重復3 次測序。 通過Contig Express 軟件去除每條序列的兩端引物區，檢查測序峰圖并進行人工校正，最終得到全部序列。

1.2.3 GenBank 序列下載及ITS2 注釋 通過GenBank（https：//www.ncbi.nlm.nih.gov/）分別下載紫蘇與其變種、紫蘇子與其各混偽品的ITS 序列，并通過BLAST 相似性搜索排除可疑序列。將所得序列連同測序序列，基于隱馬爾可夫模型HMM[18]除去位于兩端的5.8S 和28S 序列， 得到所有樣本的ITS2 序列。 使用Clustal omega 在線軟件對所有收集到的序列進行比對和矯正。

1.2.4 序列分析、 遺傳距離計算及NJ 樹構建MEGA7.0 軟件比對序列矩陣，統計GC 含量、保守位點和變異位點等信息； 基于Kimura 2-parameter（K2P）模型計算各物種種內和種間遺傳距離，采用鄰接法（Neighbor joining，NJ）構建系統發育樹，設置參數Bootstrap 自展支持率1000 次重復檢驗各分支。

1.2.5 ITS2 二級結構預測及PNJ 樹構建 利用ITS2 Database （http：//its2.bioapps.biozentrum.uni -wuerzburg.de/）中Predict 功能預測各物種ITS2 的二級結構。 下載各樣本的核酸序列和二級結構聯合矩陣，通過4Sale 軟件比對聯合矩陣，進行序列和二級結構的比對。比對完成后對ITS2 一級序列和二級結構系統發育樹進行構建，4Sale[19]導出的比對結果輸入ProfDistS[20]軟件，基于距離法利用剖面鄰接Profile neighbor-joining （PNJ） 構建系統發育樹，并自展分析1000 次來評價分支系統發育關系的可靠性。

2 結果與分析

2.1 紫蘇與其變種ITS2 變異位點分析

紫蘇與其變種的ITS2 序列長度相近， 集中在231 bp 和233 bp。 3 個日本紫蘇品種的ITS2 序列GC 含量較高在81.20%～82.00%， 而紫蘇與其他變種的序列GC 含量為64.20%～65.40%（見表3）。 紫蘇I（序列號KC011247）在123 bp 處有T→G 一個變異位點，即紫蘇種內存在2 種單倍型；回回蘇和耳齒紫蘇分別在222 bp 和22 bp 處各有一個變異位點。野生紫蘇在22、103 bp 處發生變異，共兩個變異位點。白蘇無變異位點。3 個日本紫蘇與其他變種間序列變異較大。 變異位點信息如表4 所示。

表3 紫蘇與其變種的ITS2 序列差異分析Table 3 ITS2 sequence analysis of P. frutescens and its varieties

表4 紫蘇與其變種的變異位點分析Table 4 Variation sites of P. frutescens and its varieties

2.2 紫蘇與其變種遺傳距離分析

基于K2P 模型計算的遺傳距離顯示（見表5），紫蘇的平均種內遺傳距離為0.002， 其余各變種的種內遺傳距離均為0.000。 白蘇與紫蘇的種間遺傳距離為0.001， 紫蘇的種內遺傳距離大于紫蘇與白蘇種間遺傳距離， 因此不能將白蘇和紫蘇區分開。白蘇、 紫蘇與其他變種種間的遺傳距離為0.005～0.033，大于各物種的種內遺傳距離0.000～0.002；即白蘇、紫蘇分別與其他變種及其他各物種（不含構骨）的種間最小遺傳距離大于各物種種內最大遺傳距離。 因此，從遺傳距離水平，ITS2 序列可將紫蘇、白蘇與其他變種區分開。

表5 紫蘇與其變種的遺傳距離Table 5 Genetic distance table of P. frutescens and its varieties

2.3 紫蘇與其變種二級結構及CBC 模型分析

紫蘇與其變種的ITS2 二級結構均由一環 （主環）四臂（螺旋區）組成（見圖1），并且中心環和莖臂的夾角在不同物種間較為相近。 其中紫蘇包含3 個內環和4 個發夾環，與白蘇、耳齒紫蘇和回回蘇相比在臂IV 上區少了2 個內環，與野生紫蘇、日本紫蘇I 相比在臂III 和臂IV 區上各少2 個內環， 與日本紫蘇II 相比，臂IV 區少了2 個內環。 白蘇包含有4 個發夾環和5 個內環，白蘇比野生紫蘇在臂III 區少1 個內環， 白蘇和耳齒紫蘇的臂IV 區的內環的位置也不同。 在臂IV 區處回回蘇的第1 個內環比耳齒紫蘇更靠近中心環。 回回蘇相比野生紫蘇在臂III 上少1 個內環， 并且臂IV 區的第1 個內環更靠近中心環。 即紫蘇與白蘇及其他變種之間，ITS2 二級結構在內環的數量和位置均存在差異。CBC（補償性堿基替換）模型的統計結果顯示，回回蘇分別與白蘇、耳齒紫蘇、野生紫蘇在臂Ⅳ區具有1 個CBC（203/222 T-A→G-C）和1 個hCBC（204/221 T-A→T-G）（見圖2），基于CBC 模型顯示回回蘇分別與白蘇、耳齒紫蘇、野生紫蘇符合雜交不親和現象。 綜上所述，ITS2 二級結構對于紫蘇與其變種的鑒別具有指導意義。

圖1 紫蘇與其變種的ITS2 二級結構圖Fig. 1 ITS2 secondary structure prediction diagram of P. frutescens and its varieties

圖2 紫蘇與其變種的ITS2 二級結構CBC 模型Fig. 2 CBC of ITS2 secondary structure of P. frutescens and its varieties

2.4 紫蘇與其變種NJ 樹分析

由ITS2 序列構建的NJ 樹可以看出，日本紫蘇與我國產紫蘇與其變種明顯區分開（見圖3）。 從分支情況和支持率可見，紫蘇和白蘇樣本不能有效的區分開，共聚為一個大分支，二者親緣關系最近；回回蘇、耳齒紫蘇和野生紫蘇形成姐妹分支，其中野生紫蘇與耳齒紫蘇親緣關系較近。 冬青科植物枸骨Ilex cornuta 作為外類群。

圖3 紫蘇與其變種的NJ 樹Fig. 3 NJ tree of P. frutescens and its varieties

2.5 紫蘇子與其混偽品的序列和遺傳距離分析

紫蘇子與其混偽品的ITS2 序列差異較大（見表6）， 序列長度在191～233 bp，GC 含量在50.50%～64.90%，物種間也存在大量變異位點。 如表7 所示，基于K2P 模型計算紫蘇子與其混偽品的遺傳距離，急性子、白芥子、油菜子種內遺傳距離均為0.000，紫蘇子種內平均遺傳距離為0.003， 種內最大遺傳距離的是混偽品沙苑子0.013。 種間最小遺傳距離是白芥子和沙苑子0.627； 紫蘇子與其混偽品白芥子的遺傳距離最小，而紫蘇子與其混偽品菟絲子的遺傳距離最大。所有物種的種間最小遺傳距離（0.627）大于種內最大遺傳距離（0.013），并且紫蘇子與其混偽品的種間最小遺傳距離（0.634）遠大于紫蘇子種內遺傳距離， 表明基于K2P 遺傳距離，ITS2 序列能夠區分紫蘇子與其混偽品（急性子、沙苑子、白芥子、菟絲子和油菜子）。

表6 紫蘇子與其混偽品的ITS2 序列分析Table 6 ITS2 sequence analysis of Perillae fructus and its adulterants

表7 紫蘇子與其混偽品的遺傳距離Table 7 Genetic distance table of Perillae fructus and its adulterants

2.6 紫蘇子與其混偽品的二級結構預測及分析

紫蘇子與其混偽品的二級結構符合被子植物一環（主環）四臂（螺旋區）的特征，其中臂III 最長，臂IV 變異最大，其次是臂I，臂II 最保守，螺旋上又有大小和數量不等的莖環（見圖4）。 紫蘇子與其混偽品在環的大小和位置、旋臂的角度上均有較大區別，其中急性子、沙苑子、白芥子和油菜子在臂I 和臂IV 之間有一段富含嘌呤的單鏈結構。 因此ITS2二級結構對紫蘇子與其混偽品的物種鑒別具有一定的參考價值。

圖4 紫蘇子與其混偽品的ITS2 二級結構圖Fig. 4 Prediction diagram of ITS2 secondary structure of Perillae fructus and its adulterants

2.7 紫蘇子與其混偽品的NJ 樹、PNJ 樹分析

由紫蘇子與其混偽品的NJ 樹中可見， 紫蘇子單獨聚為一支，支持率高達99%，表現出良好的單系性（見圖5）。常見混偽品沙苑子、白芥子、急性子、菟絲子和油菜子也分別聚為一支。ITS2 一級序列與二級結構聯合矩陣構建PNJ 樹顯示，各物種均能分支區分開與NJ 樹結果吻合，并且PNJ 樹結果比NJ樹顯示出更多的系統發育關系信息位點。

圖5 紫蘇子與其混偽品的NJ 樹和PNJ 樹Fig. 5 NJ tree and PNJ tree of Perillae fructus and its adulterants

3 結語

我國是紫蘇的原產地，紫蘇具有很高的食用價值和藥用價值，是藥食同源的代表植物。 紫蘇作為食材是美味的天然調味品，早在兩漢時期，已有將活鯉魚切成生魚片搭配紫蘇食用，既可去腥又可保持魚肉的鮮美。 并且紫蘇在我國用藥歷史悠久，紫蘇葉、 紫蘇梗和紫蘇子被收錄為3 種不同藥材，分別具有治療風寒感冒、調節腸胃、促進氣血循環、減輕疼痛、解咳嗽和氣喘、緩解便秘的作用。 另外，詩經中“荏染柔木，君子樹之”，已經用紫蘇指代人們美好的情感，后更用“荏苒冬春謝，寒暑忽流易”，將紫蘇與流逝的時光聯系起來。 可見，對中國人來說，紫蘇不僅是一種食物、一種藥材，更是一種中國文化精神象征。

紫蘇屬下變種多，且變種易受到環境因素的影響而改變，品種的鑒定依賴于具備植物學和生藥學等專業知識且經驗豐富的學者。 近年來，我國紫蘇市場銷勢活躍，產銷兩旺，批量成交增加，出口訂單上升。 而前幾年我國紫蘇產量少，市場供應缺口逐年加大，短期難以緩解，紫蘇的緊缺已使價格大幅上漲，導致市面上出現紫蘇子與其混偽品、紫蘇與其變種存在不合理的銷售問題。

作者采用ITS2 序列以及二級結構分別對紫蘇與其變種、紫蘇子與其混偽品進行鑒別，結果顯示，僅僅依靠ITS2 序列不能將紫蘇與白蘇兩個種區分開， 但二者的ITS2 二級結構存在內環數量的差異；紫蘇、白蘇分別與其他變種均可有效區分，符合種間最小遺傳距離大于種內最大遺傳距離。 在ITS 序列、K2P 遺傳距離、ITS2 二級結構、 基于序列的NJ樹和基于一級結構和二級結構聯合矩陣的PNJ 樹均顯示紫蘇子與其混偽品可以明顯區分。 因此，建議ITS2 可以作為紫蘇子與其混偽品的DNA 條形碼，而對紫蘇與白蘇的鑒別或建議合并的假設還有待進一步研究。