不同產區艾葉ISSR遺傳多樣性分析

余孟娟，趙博宇，董誠明，李 詢，邢 冰，李 曼

不同產區艾葉ISSR遺傳多樣性分析

余孟娟，趙博宇，董誠明*，李 詢，邢 冰，李 曼

河南中醫藥大學藥學院，鄭州 河南 450046

對不同產區艾葉進行遺傳多樣性分析，為艾葉規范化種植及新品種選育工作奠定基礎。采用ISSR分子標記技術對36個不同產區的艾葉樣品進行擴增，然后利用POPGENEv1.32軟件進行遺傳學分析，得出它們的遺傳特征系數，用Ntsyspc2.1對艾葉ISSR分子標記結果進行UPGMA聚類分析。篩選出16條引物，共擴增出條帶402條，其中多態性條帶398條，多態性百分比為99%；在對栽培種和野生種之間遺傳距離分析后發現，兩者之間的遺傳一致度為0.989 1，遺傳距離為0.011 0；對按緯度劃分的4個區域艾資源進行遺傳結構與遺傳變異分析，4個區域間基因多樣性指數（diversity index，t）為0.192 7、遺傳多樣性指數（genetic diversity index，s）為0.164 2、遺傳變異指數（genetic variation index，st）為0.147 9、基因流（gene flow，m）為2.881 8，NTSYSpc軟件對4個區域進行系統聚類分析，可將4個區位劃分為2大類，且4區位相鄰2個之間均具有一定的遺傳相似性。艾葉具有極其豐富的遺傳多樣性，栽培種和野生種之間親緣關系接近、遺傳背景差異較小。艾葉其遺傳特性隨著緯度的變化呈現出一種遞進式遺傳擴散，其中河南南部居群多樣性最為豐富，其次為河南中西部和河南中部，北部居群遺傳多樣性最低。

艾葉；ISSR分子標記；遺傳多樣性；栽培；野生；不同緯度

艾葉為菊科植物艾LevletVant的干燥葉，性溫，味苦、辛，有溫經止血、散寒止痛之功效。在臨床上可分為外用和內服2類，其中外用以火灸著名，而內服用途甚廣。我國艾分布范圍較廣，從東北到華南均有生長[1]。國內外對艾葉成分含量、藥理藥效及化學成分的提取均有比較深入的研究[2-3]，但利用分子生物學手段對不同產地艾葉種質資源遺傳多樣性的探索卻較少[4-5]。遺傳多樣性的本質是遺傳物質變異，探索遺傳多樣性是研究植物進化和親緣關系的基礎[6]，故將現代分子生物技術應用于艾規范化種植及輔助育種進而發掘優良資源或改良資源顯得尤為重要[7]。王惠君等[5]在2015年首次將分子標記方法應用到艾葉的種質資源多樣性研究領域，但其卻未能收集到北艾的種質資源，分析其遺傳多樣性，在《國藥提要》《中國藥學大詞典》以及國外的《泰西本草名疏》中均以北艾作為艾葉的正品基原，本實驗所用樣品包括北艾，填補了北艾種質資源多樣性研究的空白。

簡單重復序列間擴增（inter-simple sequence repeat，ISSR）是基于微衛星系列發展來的分子標記技術，具有快速高效、引物特異性強、多態性好等優點，目前被廣泛應用于品種鑒定、物種分類、居群遺傳學研究等領域，是研究植物遺傳多樣性常用的分子生物學方法[8-10]。因此本研究擬采用ISSR分子標記技術對我國傳統產區的36份艾葉從分子水平分析其遺傳多樣性，旨在從DNA水平上探明艾葉資源遺傳多樣性現狀；野生艾葉和栽培艾葉遺傳多樣性差距以及緯度與艾葉遺傳多樣性的關系，為今后艾的種質資源利用和品種選育等奠定基礎[11-12]]，因而具有非常重要的現實意義。

1 材料和儀器

1.1 材料

查閱《中國植物志》《中藥志》《河南植物志》《河南志》等相關資料和文獻，了解艾葉生境、分布、栽培等信息，以北緯38°20′40″到北緯33°03′6″劃分采樣區域，南起湖北蘄春，北至河北安國，基本覆蓋了全國艾葉主要產地。經河南中醫藥大學董誠明教授鑒定為艾Levl.et Vant.，取樣部位均為其幼嫩葉片，總計36份樣品（表1），野生種和栽培種各有18種。所有材料在采樣前均經過了去雜處理，以確保種質純度。采集后封裝于有干燥劑的自封袋內，放置入冰盒，并及時送至實驗室用錫箔紙包裹浸入液氮，迅速冷凍后放入于?80 ℃冰箱保存備用。

表1 樣品生長環境信息

1.2 試劑與儀器

植物組織基因組DNA抽提試劑盒（生工生物公司，編號B518261）、SanPrep柱式DNA膠回收試劑盒（生工生物公司，編號B518131）、Taq Plus DNA 聚合酶（BBI，編號B600090）、滅菌水。

SW-CJ-1D型潔凈工作臺（江蘇蘇潔凈化設備廠）、PCR-96型PCR反應擴增儀（BBI）、G508009型高速微量離心機（100～14 800 r/min，生工生物公司）、DYY-6C型電泳儀（北京六一公司）、FR980型凝膠成像系統（上海復日科技有限公司）。

2 方法

2.1 DNA提取與檢測

采用試劑盒（生工生物，批號B518261）法提取樣品葉片基因組DNA。

2.2 ISSR引物篩選及PCR反應體系確定

以提取的艾基因組DNA為模板，從哥倫比亞大學（UBC）公布的100條引物[13]中篩選出16條引物，均能擴增出清晰完整、多態性和穩定性均較好的條帶。通過一系列實驗優化，最終確定本實驗的ISSR-PCR反應體系為（50 μL）：10×PCR Buffer 5 μL，引物F 2 μL，引物R 2 μL，Dntp 2 μL,Taq Plus DNA Polymerase 0.5 μL，最后用滅菌水加至50 μL。擴増反應程序：95 ℃預變性5 min；95 ℃變性30 s，52 ℃退火30 s，72 ℃延伸50 s，進行35個循環；72 ℃延伸8 min；4 ℃保存。

2.3 ISSR-PCR擴增及檢測

將需要使用的反應體系取出，對所有樣品基因組DNA擴增，將擴增后的產物用2%的瓊脂糖凝膠在150 V 100 mA的電泳條件下電泳20 min，置于全自動凝膠成像分析系統上成像，調節至最佳光圈及亮度后拍照保存。

2.4 數據處理與分析

將實驗中所做出來的瓊脂糖凝膠膠板拍照，對膠板中的電泳條數進行判讀。統計原則以瓊脂糖凝膠電泳上的每1個條帶當成1個分子標記，其所代表的意義為1個基因位點。對于不同產地的艾葉DNA樣品在瓊脂糖凝膠電泳統計時按照有條帶為1，無條帶為0；相同相對分子質量中，以較為清晰的記為1不清晰的記為0，模糊條帶記為0的原則，在EXCEL中記錄成0、1矩陣形式。用Ntsyspc2.1對艾葉ISSR分子標記統計好的0、1矩陣數據做UPGMA聚類分析。用POPGENEv1.32軟件對艾葉ISSR分子標記統計好的0、1矩陣數據進行遺傳學分析，得出它們的遺傳特征系數。

3 結果與分析

3.1 引物擴增結果分析

以16條引物對36個不同產地的艾葉資源進行ISSR-PCR擴增（表2），共擴增出402條清晰的條帶，其中多態性條帶有398條，多態性百分比為99%。16條引物擴增出的有效條帶數范圍為17～31條，均值為24.875條。引物UBC840擴增出的條帶最多，為33條，引物UBC822擴增條帶最少，為17條（圖1～3）。ISSR擴增結果表明艾種質資源具有較為豐富的遺傳多樣性。

表2 ISSR引物擴增結果

圖1 引物892瓊脂糖凝膠擴增結果

圖2 引物822瓊脂糖凝膠擴增結果

圖3 引物834瓊脂糖凝膠擴增結果

3.2 不同產區艾野生種質與栽培種質遺傳多樣性分析

對ISSR引物電泳結果進行統計后，通過POPGENEv1.32軟件對艾種質資源栽培種與野生種的遺傳背景進行分析，從2個資源的遺傳多樣性指數來看，遺傳多樣性指數栽培艾要略高于野生品，但遺傳多樣性差距并不顯著。等位基因數（number of alleles，a）、有效等位基因數（effective number of allele，e）、Nei’s 氏基因多樣性指數（Nei’s gene diversity index，h）、Shannon’s 信息指數（Shannon's information index，）等指標結果見表3。對栽培種和野生種之間遺傳距離分析后發現，兩者之間的遺傳一致度接近于1.000 0，為0.989 1，遺傳距離較小，為0.011 0，表明兩者親緣關系接近，說明栽培種來自于野生種的馴化，且馴化時間并不很長，因此遺傳背景差異較小[14]。對人工栽培艾與野生艾種質資源進行遺傳結構與遺傳變異分析可以發現，有2.77%的遺傳分化存在于2個亞群間，在采樣過程中能夠觀察到同一采集點的每株艾的表型均略有不同，生物的表觀多樣性主要由基因決定，基因差異明顯則表觀形態多樣性會更加豐富，由于艾的人工馴化時間短，栽培品種并非嚴格意義上遺傳背景高度純合的品種，而是多為一些混合品種，故種內遺傳分化程度高則使得在同一樣點單株表觀形態特性呈現多樣性。栽培種與野生種之間基因流（gene flow，m）為17.520 8，由于艾的繁殖更多是通過無性繁殖途徑實現，此處的基因流更多是反映野生種與栽培種兩者之間遺傳多樣性趨同的一種關系，即兩者遺傳背景差距并不大[15-17]。

表3 人工栽培艾與野生艾種質資源遺傳多樣性

3.3 不同產區（不同緯度）艾遺傳多樣性分析

按照地理緯度，將本次采集的緯度差異相對較小的產地劃分為4個大的區域（表4），從這些樣本中共擴增出396條多態性較為豐富的條帶，多態性位點占比98.51%，該擴增結果表明不同緯度間艾具有較好的遺傳多樣性。通過POPGENEv1.32軟件對按表4方法劃分的不同區域艾種質資源進行分析，結果如表5所示，從遺傳多樣性指數來看，4個區域的艾葉之間具有顯著的遺傳多樣性。對4個區域資源進行遺傳結構與遺傳變異分析后發現，4個區域間艾種群基因多樣性指數（gene diversity index，t）為0.192 7、遺傳多樣性指數（genetic diversity index，s）為0.164 2、遺傳變異指數（genetic variation index，st）為0.147 9、m為2.881 8。由此表明在地理區劃上4個居群間的艾品種具有一定的基因交流。在地理區劃上4個居群間的艾葉具有一定的基因交流。對4個區域間遺傳距離分析后發現（表6），它們之間的遺傳一致度接近于1，遺傳距離較小，接近于0，表明這4個居群間親緣關系極近。再次證明了4個居群具有一定的聯系。利用PopGen對4個區域進行聚類，由聚類圖（圖4）可知，中西部ZXB與南部NB親緣關系更為接近，被聚成第1類。中部ZB次之，與第1類單獨聚成第2類。北部BB遺傳分化最強，單獨與第2類聚為一類。按地理劃分的聚類結果顯示出不同區域間由南向北依次具有進階相似性，即不同產區艾在緯度從南向北相互兩兩具有聯系，這表明各地區艾的親緣關系受自身生長的緯度和相鄰地域間的距離影響。

表4 不同產區不同緯度區劃

表5 不同產區區劃艾種質資源遺傳變異結果

表6 不同緯度區域遺傳距離相似度

圖4 不同緯度區域聚類結果

4 討論

以篩選出的16條引物對36個不同產地的艾葉資源進行擴增，多態性百分比為99%。有效條帶數范圍為17～31條，均值為24.875條。擴增結果表明艾種質資源具有較為豐富的遺傳多樣性。在采集樣品的過程中也能夠看出田間艾草臨株之間的差異性，表明即便在同一地點生長的不同植株間也存在一定的差異[18-19]。植物遺傳多樣性的成因較為復雜，植物在不同環境下與不同異栽培物之間存在協同共存，通過蟲媒等繁殖途徑，與其他異栽培物進行接觸，導致其出現了較高的遺傳多樣性。

從栽培種與野生種2個類型來看，其兩種生存狀態下的艾均具有豐富的遺傳多樣性。兩者之間的遺傳一致度接近于1，遺傳距離也很小，趨近于0，這表明人類對于艾的馴化歷史相對較短，人工栽培種和野生種之間尚未分化完全。由于艾的繁殖更多是通過無性繁殖途徑實現，因此栽培種與野生種之間存在的m為17.520 8時更多是反映野生種與栽培種兩者之間遺傳多樣性趨同的一種關系。即兩者遺傳背景差距并不大，親緣關系接近，有利于植物間育種。本研究在將緯度差異相對較小的采集地劃分為的4個大的區域時，發現其遺傳多樣性整體呈現出一種遞進式的遺傳擴散。由聚類結果表明，相鄰區域間艾葉兩兩聯系，說明艾葉遺傳距離與緯度有關，從側面也能反映出不同種質艾葉間基因交流程度較大。其中，南部居群多樣性最為豐富，中西部和中部分別居于第2和第3位置，北部居群遺傳多樣性在4個居群中最低。根據艾葉4個居群的遺傳分析可以推斷出，艾的分布是由南部向中部、北部依次擴散的。性狀表型是基因型和環境間互做的結果，河南中西部屬于伏牛山脈分布區域，其氣候多樣性強，生態類型復雜多變，由此為不同產區艾的生存和變異提供了物質基礎，而北部居群所在區域受北部大陸季風影響較強，氣候類型較為單一，這也是其遺傳類型較少的主要原因[20-21]。

綜上可得，艾葉的遺傳多樣性與地理距離具有一定的相關性，種質資源遺傳多樣性是育種的重要內容之一，種質資源遺傳多樣性越豐富越有利于新品種選育，種質親緣關系越近越有利于植物間育種。本研究通過分子標記技術為艾葉規范化種植及新品種選育工作奠定了一定的基礎。

利益沖突 所有作者均聲明不存在利益沖突

[1] 宋川. 艾葉生藥及化學成分的研究 [D]. 昆明: 云南中醫學院, 2013.

[2]田璐. 艾葉化學成分分析及其抗炎功效研究 [D]. 廣州: 暨南大學, 2017.

[3]于曉. 艾葉化學成分及一測多評分析方法的研究 [D]. 濟南: 山東中醫藥大學, 2016.

[4]焦倩. 艾葉和野艾葉指紋圖譜及多成分含量測定研究 [D]. 石家莊: 河北醫科大學, 2018.

[5]王惠君, 盧誠, 黎明, 等. 傳統中藥艾種質資源的ISSR遺傳多樣性分析 [J]. 貴州農業科學, 2015, 43(8): 1-5.

[6]李冰敏. 分子標記在園林植物育種中的應用 [J]. 現代園藝, 2019(15): 148-149.

[7]Nasim N, Sandeep I S, Sahoo A,. Population genetic structure and diversity analysis in economically important(Forssk.) Kuntze accessions employing ISSR and SSR markers [J]., 2020, 143: 111894.

[8]李璐, 董誠明, 朱畇昊, 等. 基于ISSR的連翹遺傳多樣性研究 [J]. 中草藥, 2019, 50(15): 3673-3680.

[9]Shahabzadeh Z, Mohammadi R, Darvishzadeh R,. Genetic structure and diversity analysis of tall fescue populations by EST-SSR and ISSR markers [J]., 2020, 47(1): 655-669.

[10]姜武, 吳志剛, 陶正明, 等. 基于ISSR和SRAP標記的梔子種質遺傳多樣性研究 [J]. 中草藥, 2019, 50(2): 510-516.

[11]Sudha G S, Ramesh P, Sekhar A C,. Genetic diversity analysis of selected Onion (L.) germplasm using specific RAPD and ISSR polymorphism markers [J]., 2019, 17: 110-118.

[12]俞文灝, 吳保鋒, 劉勇波. 生境破碎化對動植物遺傳多樣性的影響研究進展 [J]. 應用與環境生物學報, 2019, 25(3): 743-749.

[13]周延清. DNA分子標記技術在植物研究中的應用 [M]. 北京: 化學工業出版社, 2005.

[14]梁文匯, 劉凱, 黃開順, 等. 肉桂家系遺傳背景的ISSR分析 [J]. 廣西林業科學, 2016, 45(1): 35-39.

[15]趙琬玥, 冷秋思, 屈燕, 等. 全緣葉綠絨蒿() 的ISSR遺傳多樣性分析 [J]. 分子植物育種, 2019, 17(20): 6891-6899.

[16]邢婀娜, 徐詩濤, 任明迅. 海南特有毛花馬鈴苣苔小尺度局域種群的年齡結構與基因流 [J]. 熱帶生物學報, 2018, 9(1): 37-46.

[17]刁毅, 周麗娟. 地木耳遺傳多樣性的ISSR分析 [J]. 北方園藝, 2018(14): 143-147.

[18]宋鳳. 基于SCoT分子標記的甘草屬植物自然居群的遺傳多樣性與基因漸滲研究 [D]. 石河子: 石河子大學, 2017.

[19]趙博宇, 董誠明. 河南艾不同種質資源形態特征比較 [J]. 中醫學報, 2019, 34(5): 1024-1029.

[20]彭國清, 唐紹清. 喀斯特耐陰灌木淡黃金花茶的精細空間遺傳結構與基因流 [J]. 生態學報, 2017, 37(21): 7313-7323.

[21]吳蕊. 紫牡丹引種馴化及遠緣雜交后代早期鑒定研究 [D]. 北京: 中國農業科學院, 2011.

ISSR analysis on genetic diversity offrom different habitats

YU Meng-juan, ZHAO Bo-yu, DONG Cheng-ming, LI Xun, XING Bing, LI Man

College of Pharmacy, Henan University of Traditional Chinese Medicine, Zhengzhou 450046, China

The genetic diversity analysis ofin different producing areas will lay the foundation for the standardized planting and selection of new varieties ofISSR molecular marker technology was used to amplify artemisia sylvestris samples from 36 different production areas, and then genetic analysis was performed using POPGENE v1.32 software to obtain their genetic characteristic coefficients. Ntsyspc2.1 was used to perform ISSR molecular marker results UPGMA cluster analysis.16 primers were selected, and a total of 402 bands were amplified, including 398 polymorphic bands, and the percentage of polymorphism was 99%.. After analyzing the genetic distance between the wild species and the wild species, it was found that the genetic identity between the two is 0.989 1 and the genetic distance is 0.011 0. The diversity index (t) is 0.192 7, the genetic diversity index (s) is 0.164 2, the genetic variation index (st) is 0.147 9, and the gene flow (m) is 2.881 8. NTSYSpc software performs a systematic cluster analysis on 4 regions. The 4 loci are divided into 2 categories, and there are certain genetic similarities between the two adjacent 4 loci.Artemisia oleifera is extremely rich in genetic diversity, the genetic relationship between cultivated and wild species is close, and the genetic background is small. The genetic characteristics ofexhibit a progressive genetic spread with changes in latitude. Among them, southern Henan has the most abundant population diversity, followed by central and western Henan and central Henan, and northern population has the lowest genetic diversity.

Artemisiae Argyi FoliumISSR molecular marker; genetic diversity; cultivated germplasm; wild germplasm; different latitudes

R286.12

A

0253 - 2670(2023)13 - 4306 - 06

10.7501/j.issn.0253-2670.2023.13.023

2023-02-06

國家自然科學基金項目（81603232）

余孟娟（1997—），女，碩士研究生，從事中藥資源保護與開發研究。E-mail: 1600213809@qq.com

董誠明（1963—），E-mail: dcm663@sina.com

[責任編輯 時圣明]