歐金梅，楊亞湉，錢程程，吳 瑞，王 瑞，黃璐琦

? 藥材與資源?

基于葉綠體基因的藥用梅群體遺傳學(xué)研究

歐金梅1, 2, 3，楊亞湉1，錢程程1，吳 瑞1，王 瑞4，黃璐琦2*

1. 安徽中醫(yī)藥大學(xué)，安徽 合肥 230012 2. 中國(guó)中醫(yī)科學(xué)院 中藥資源中心，北京 100700 3. 安徽省中醫(yī)藥科學(xué)院 中藥資源保護(hù)與開發(fā)研究所，安徽 合肥 230012 4. 蕪湖職業(yè)技術(shù)學(xué)院，安徽 蕪湖 241006

對(duì)藥用梅開展群體遺傳學(xué)研究，揭示其遺傳變異情況、單倍型地理分布格局，推測(cè)其起源和多樣化中心。采用葉綠體基因片段（A-H、16、S-G、1b）對(duì)我國(guó)23個(gè)地區(qū)的藥用梅野生群體進(jìn)行測(cè)序，開展群體遺傳多樣性、遺傳結(jié)構(gòu)和歷史動(dòng)態(tài)分析。藥用梅共檢測(cè)出20個(gè)單倍型，單倍型多樣性（d）為0.790，核苷酸多樣性（π）為0.001 71。AMOVA分析顯示，藥用梅在cpDNA水平上的遺傳變異主要發(fā)生在群體間，表明該物種具有明顯的譜系地理結(jié)構(gòu)，群體間基因流主要以花粉流為主。中性檢驗(yàn)和失配分布分析表明藥用梅群體歷史是基本穩(wěn)定的，沒有經(jīng)歷明顯的快速擴(kuò)張事件。單倍型中H2是分布最廣泛的共享單倍型，H1、H15位于單倍型網(wǎng)絡(luò)圖中心位置。藥用梅具有明顯的譜系地理結(jié)構(gòu)，群體遺傳性較高，四川、福建作為藥用梅主產(chǎn)區(qū)，也是其起源和多樣化中心。

梅；葉綠體基因；群體遺傳學(xué)；譜系地理學(xué)；A-H；16；S-G；1b

梅Sieb. et Zucc.，薔薇科李屬多年生木本，我國(guó)特有的藥食同源植物，其果實(shí)、花作中藥烏梅、梅花使用，均收錄于《中國(guó)藥典》2020年版中，烏梅可斂肺澀腸、生津安蛔，具有抗癌、抗氧化、抗肝纖維化等藥理活性[1-2]。梅在我國(guó)已有2000多年的藥用歷史，四川、云南、福建等地為藥用梅的主產(chǎn)區(qū)，其自然群體分布于長(zhǎng)江流域和整個(gè)江南地區(qū)[3]。課題組調(diào)查發(fā)現(xiàn)，藥用梅主產(chǎn)區(qū)栽培品種近百余種，產(chǎn)地盲目引種造成品種多而復(fù)雜，對(duì)藥材質(zhì)量和臨床療效產(chǎn)生很大影響[4-5]。藥用梅優(yōu)良種質(zhì)選育是目前亟待解決的問題。

葉綠體基因（cpDNA）具有單親遺傳、相對(duì)保守的特點(diǎn)。被子植物多為母系遺傳，依靠種子流進(jìn)行基因交流，其基因編碼區(qū)和非編碼區(qū)的變異速率差異較大，能真實(shí)反映出物種的遺傳分化、譜系結(jié)構(gòu)和擴(kuò)散情況[6-7]。近年來，基于cpDNA的分子譜系地理學(xué)早期主要集中研究珍稀瀕危植物的進(jìn)化歷史和保護(hù)，如十齒花、綿參、紅豆杉等[8-11]。近年來，分子譜系地理學(xué)在藥用植物遺傳多樣性、地理分布格局、瀕危機(jī)制等方面研究廣泛深入，如姜丹[12]應(yīng)用3個(gè)cpDNA序列4-A、AJ和AH對(duì)28個(gè)黃芩野生居群和1個(gè)甘肅黃芩野生居群進(jìn)行分子譜系研究，表明黃芩和甘肅黃芩之間存在顯著的遺傳分化，種間沒有基因滲透和雜交。張玉君[13]利用H-A、L-F、K序列對(duì)31個(gè)山藥居群進(jìn)行分析，得出山藥的遺傳多樣性主要在群體間，并推測(cè)出栽培薯蕷和參薯的馴化起源。

遺傳多樣性是生物多樣性的重要組成部分，是物種能夠長(zhǎng)期生存且保持進(jìn)化能力的物質(zhì)基礎(chǔ)。野生群體大都保留著較高的遺傳多樣性[14-16]，是栽培育種的種質(zhì)資源庫(kù)。本研究利用課題組前期從梅的葉綠體全基因組中篩選出的4對(duì)引物S-G、A-H、1b、16，用于藥用梅野生群體的遺傳多樣性和遺傳結(jié)構(gòu)分析，為藥用梅的種質(zhì)資源保護(hù)和品種選育提供分子依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

本研究共采集到23個(gè)地區(qū)的藥用梅野生群體樣品，分布于中國(guó)四川、云南、貴州、安徽、福建、西藏、江西地區(qū)，經(jīng)安徽中醫(yī)藥大學(xué)俞年軍教授鑒定為梅Sieb. et Zucc.，憑證標(biāo)本保存在安徽中醫(yī)藥大學(xué)中藥資源中心。每個(gè)群體采集10個(gè)以上個(gè)體（生境中不足10個(gè)個(gè)體全部采集），個(gè)體之間的間隔不小于10 m，采集新鮮健康葉片，用硅膠干燥后儲(chǔ)存在?70 ℃下超低溫冰箱備用。采集信息見表1。

表1 野生樣本采集信息

Table 1 Detailed information of samples in this study

序號(hào)群體編號(hào)數(shù)量采集地點(diǎn)經(jīng)度緯度海拔/m 1AHSW 4安徽寧國(guó)市仙霞鎮(zhèn)E119°18'27''N30°25'11''220 2AHHS 4安徽黃山市湯口鎮(zhèn)E118°08'47''N30°04'59''550 3FJSD11福建上杭縣臨城鎮(zhèn)E116°26'53''N24°57'40''230 4FJXY 5福建永泰葛嶺鎮(zhèn)E119°00'95''N25°83'55''150 5GZZJ12貴州赫章縣六曲河鎮(zhèn)E104°04'23''N27°16'15''1760 6GZNZ12貴州赫章縣農(nóng)莊村E104°51'25''N27°08'58''1720 7GZTJ10貴州赫章縣孫家院鄉(xiāng)E104°40'07''N27°08'57''1530 8GZML13貴州荔波縣茂蘭鎮(zhèn)E108'01'26''N25°18'09''570 9GZBK18貴州荔波縣茂蘭鎮(zhèn)E108'00'20''N25°16'33''570 10GZWZ10貴州荔波縣茂蘭鎮(zhèn)E108'01'24''N25°19'45''610 11JXBS 4江西景德鎮(zhèn)瑤里鎮(zhèn)E117°35'58''N29°33'35''150 12JXML10江西景德鎮(zhèn)瑤里鎮(zhèn)E117°38'38''N29°32'24''210 13SCRZ22四川寶興縣磽磧鄉(xiāng)E102°43'47''N30°41'60''2150 14SCCJ20四川達(dá)州市蔡家坡村E107°29°51''N30°59'60''510 15SCDZ10四川德昌縣德州鎮(zhèn)E102°10'26''N27°24'22''1375 16SCLD 6四川甘孜藏族自治區(qū)瀘定縣E102°12'58''N29°41'52''1680 17SCHL15四川會(huì)理縣E102°12'54''N26°34'25''1811 18SCXY10四川平武縣平通鎮(zhèn)E104°40'23''N32°05'13''1040 19SCMZ10四川馬邊縣梅子壩鄉(xiāng)E103°21'45''N 28°53'32''1274 20SCSH10四川馬邊縣三河口鄉(xiāng)E103°22'32''N 28°53'32''1377 21YNTC 5云南寧蒗縣大興鎮(zhèn)E100°53'11''N27°15'46''2370 22YNDY 5云南麗江得一公司種質(zhì)資源圃E100°15'09''N26°49'31''2390 23XZLZ10西藏林芝縣E94°25'04''N29°23'05''2940

1.2 基因組DNA提取

樣品葉片采用改良CTAB法提取基因組總DNA。經(jīng)1%瓊脂糖凝膠電泳檢測(cè)質(zhì)量和濃度后，調(diào)整濃度為50 ng/μL，保存于?20 ℃冰箱待用。

1.3 cpDNA引物篩選及群體擴(kuò)增

采用課題組前期篩選的A-H、16、S-G、1b序列（表2），每個(gè)群體選10個(gè)樣品（不足10個(gè)全選）進(jìn)行擴(kuò)增。用1%的瓊脂糖凝膠電泳檢測(cè)擴(kuò)增產(chǎn)物，送至上海生工生物工程有限公司（上海，中國(guó)）進(jìn)行雙向測(cè)序。

PCR反應(yīng)體系：模板DNA 1 μL，正反引物10 μmol/L各0.5 μL，MIX酶5 μL，補(bǔ)加ddH2O至總體系為10.00 μL。擴(kuò)增程序：95 ℃預(yù)變性4 min，95 ℃變性30 s，55 ℃退火30 s，72 ℃延伸1 min，共35次循環(huán)，72 ℃延伸6 min。擴(kuò)增產(chǎn)物置于?20 ℃冰箱保存。

表2 葉綠體片段引物信息

Table 2 Information of polymorphic cpDNA prime pairs

編號(hào)引物名稱引物序列（5’-3’）退火溫度/℃序列長(zhǎng)度/bp 1psbAGGAAGTTATGCATGAACGTAATG55380 trnHTGGTGGATTCACAATCCACTG 2rps16-1FTTTGCCTTGAATCATTGGGT55420 rps16-1DCGATGTGGTAGAAAGCAACG 3trnSCTTTAGTCCACTCAGCCATCT55672 trnGCTATACCCGCTACAATCC 4ycf1b FTCTCGACGAAAGTCCGGTTGTTGTGA55890 ycf1b RATACATGTCAAAGTGATGGAAAA

1.4 群體單倍型和遺傳多樣性分析

測(cè)序的原始數(shù)據(jù)，用GENEIOUS 10.2.6軟件進(jìn)行編輯、比對(duì)和校正，最后將4個(gè)cpDNA片段按照相同的樣品名拼接成一條序列。群體中具有多態(tài)性位點(diǎn)的序列數(shù)目（）可以簡(jiǎn)單反映群體的遺傳多樣性大小。利用DnaSP 5.10軟件[17]進(jìn)行單倍型個(gè)數(shù)（Hap）、核苷酸多態(tài)性（π）、單倍型多樣性（d）和多態(tài)性位點(diǎn)的統(tǒng)計(jì)。在考慮缺失的情況下將群體單倍型分布以及共享情況，利用POPART v.4.8.4軟件[18]進(jìn)行單倍型網(wǎng)絡(luò)圖構(gòu)建，在ArcGIS10.5軟件中繪制單倍型地理分布圖。利用RAxML8.2.4軟件，以最大似然法（maximum likelihood estimation，ML）構(gòu)建系統(tǒng)進(jìn)化樹。

1.5 群體遺傳結(jié)構(gòu)分析

基于cpDNA單倍型數(shù)據(jù)，利用DnaSP5.10軟件計(jì)算遺傳分化系數(shù)（st和st）?；?000次置換進(jìn)行檢驗(yàn)，通過比較st和st的大小，對(duì)藥用梅群體結(jié)構(gòu)進(jìn)行分析。利用ARLEQUIN v3.5.2.2軟件[19]計(jì)算群體間遺傳分化系數(shù)（st），采用AMOVA分析藥用梅群體內(nèi)和群體間遺傳變異程度，通過10 000次重復(fù)來檢測(cè)遺傳組成的顯著性水平。

1.6 群體歷史動(dòng)態(tài)分析

采用DnaSP 5.10軟件進(jìn)行單倍型的失配分布分析，繪制失配分布曲線，以推測(cè)物種水平、譜系地理分子水平間群體是否經(jīng)歷過動(dòng)態(tài)擴(kuò)張。同時(shí)進(jìn)行Tajinma’s D[20]和Fu's Fs[21]中性檢驗(yàn)，以推測(cè)葉綠體片段是否符合中性進(jìn)化模式。

2 結(jié)果與分析

2.1 cpDNA序列測(cè)序結(jié)果

采用4個(gè)葉綠體片段對(duì)23個(gè)群體193個(gè)野生個(gè)體進(jìn)行擴(kuò)增后，經(jīng)拼接、編輯和比對(duì)后共得到長(zhǎng)度為1 945 bp的核苷酸矩陣，共檢測(cè)到27個(gè)變異位點(diǎn)，其中有19個(gè)突變，8個(gè)插入或缺失。A-H片段為334 bp（1～334 bp），有1個(gè)單點(diǎn)突變，4個(gè)插入或缺失；S-G片段為355 bp（335～689 bp），有6個(gè)單點(diǎn)突變，2個(gè)插入或缺失；16片段為499 bp（690～1188 bp），有7個(gè)單點(diǎn)突變，1個(gè)插入或缺失；1b片段為757 bp（1189～1945 bp），有5個(gè)單點(diǎn)突變，1個(gè)插入或缺失，具體見表3。

2.2 單倍型及遺傳多樣性

共檢測(cè)出20個(gè)單倍型（H1～H20）。H2是分布最廣泛的共享單倍型，出現(xiàn)我國(guó)福建、四川和貴州的10個(gè)居群中；H1、H5和H14分布次之，其余均為特有單倍型。從居群分布來看，江西景德鎮(zhèn)（JXBS）、安徽黃山（AHHS）、福建永泰（FJXY）、四川瀘定（SCLD）、云南寧蒗（YNTC）和云南麗江（YNDY）各含一個(gè)特有單倍型，分別為H6～8、H11、H12、H18。貴州赫章（GZTJ）和西藏林芝（XZLZ）各有2個(gè)特有單倍型。四川達(dá)州（SCCJ）和四川馬邊（SCMZ）居群各有3個(gè)單倍型，福建上杭（FJSD）和四川德昌（SCDZ）各有4個(gè)單倍型，表明這4個(gè)居群的單倍型類型比較豐富。從單倍型網(wǎng)絡(luò)圖（圖1）上看出，H1、H15位于網(wǎng)絡(luò)圖的中心位置，推測(cè)H1、H15可能是藥用梅的原始單倍型。H1單倍型分布在福建、貴州和四川居群，H15單倍型分布在四川德昌居群，因此，福建、貴州和四川地區(qū)的居群可能是藥用梅的原始單倍型。

表3 藥用梅的葉綠體片段變異位點(diǎn)和單倍型分類

Table 3 Chloroplast DNA sequences polymorphisms detected in P.mume

核苷酸位點(diǎn)trnH-psbArps16trnS-trnGycf1b 9411410711718025528420328030634334634747965372414589637 H1TAAAACACTAATTTTCGTG H2TAAAACACTAATTTTCGTA H3TAAAACACTAATTTTTGTG H4TAAAACACTAAATTTCGTG H5TAAAACACTAATTTTCGAG H6TAAAACATTAATTTTCGAG H7TAAATCACTAATTTTCGAG H8TAAAACATTAATTTACTTG H9TAAAACACTAATTTACGTG H10TAAAACACTAATTTACTTG H11TAATATACTAAATCTTGTG H12TTAAACACTAAAACTCGTG H13GTAAATACTAAATCTCGTG H14GTAAATACAAAATCTCGTG H15TTAAATACTAAATCTCGTG H16TTAAATACTACATCTCGTG H17TAAAATACTACAACTCGTG H18TTCAATACTAAATCTCGTG H19TTATATACTAAAACTCGTG H20TAAAATGCTGATTCTCGTG

從單倍型網(wǎng)絡(luò)圖中可以反映出單倍型之間的親緣關(guān)系。以祖先單倍型H1、H15為中心，H1經(jīng)過一步變異演化為5個(gè)單倍型（H2～5和H9）；H15經(jīng)過一步演化成4個(gè)單倍型（H13、H16、H18和H19），這些單倍型經(jīng)過一步或多步變異又得到其他的單倍型，其中H20單倍型在分化進(jìn)程中經(jīng)歷了H1和H5的多次變異，H19、H12單倍型經(jīng)歷了H4的2步變異和H15的一步變異，表明這些單倍型在進(jìn)化過程中群體間可能經(jīng)歷了基因交流。以日本晚櫻Matsum.和稠李L.作為外類群，利用Mega軟件基于ML法構(gòu)建了20個(gè)單倍型的系統(tǒng)發(fā)育樹（圖2）。結(jié)果顯示ML樹分為2支，以H7、H5、H6、H3、H2、H1、H9、H8、H10和H4為一支，H16、H11、H13、H14、H18、H12、H19、H15、H17和H20為一支，與單倍型網(wǎng)絡(luò)圖結(jié)果基本一致。從單倍型分布區(qū)域來看，云貴高原的西藏、云南居群分為一支，華東地區(qū)的福建、江西、安徽和貴州（GZTJ除外）居群聚在一起。而四川居群則分散在2個(gè)分支中。四川是烏梅的主產(chǎn)區(qū)，分布的單倍型分布最廣，分別滲透到西南的云貴高原和華東地區(qū)，說明四川可能是藥用梅的起源中心。

圖1 基于cpDNA的藥用梅野生群體單倍型網(wǎng)絡(luò)圖

圖2 基于20個(gè)單倍型的cpDNA序列構(gòu)建的系統(tǒng)發(fā)育樹(ML)

2.3 群體遺傳結(jié)構(gòu)

DnaSP分析結(jié)果顯示，在群體水平上，藥用梅群體總的H為0.790，π為0.001 71，表明整個(gè)野生群體遺傳多樣性較高。在單個(gè)居群中，福建上杭，貴州赫章、荔波，四川達(dá)州、德昌和馬邊地區(qū)的遺傳多樣性較高，其中以福建上杭（FJSD）居群含有4個(gè)不同的單倍型而呈現(xiàn)出單倍型多樣性最高（0.800）。安徽、江西等15個(gè)居群僅有一個(gè)單倍型，其單倍型多樣性和核苷酸多態(tài)性均為0，結(jié)果見表4。

遺傳變異分析表明，群體間的遺傳分化系數(shù)st為0.673（0.12），st為0.891（0.03），st值大于st值，且群體間基因流（m）為0.03。通過AMOVA分析群體內(nèi)和群體間的變異情況，結(jié)果顯示群體間的遺傳分化系數(shù)（st）為0.891，群體間的遺傳變異占89.10%，群體內(nèi)的遺傳變異占10.90%。以上數(shù)據(jù)表明，藥用梅野生群體間存在譜系地理結(jié)構(gòu)，地理區(qū)域相近的群體分布的單倍型親緣關(guān)系也較近，群體的遺傳變異主要來自于群體間，在葉綠體基因水平上幾乎不存在基因流。

表4 基于4個(gè)葉綠體片段聯(lián)合的藥用梅野生群體遺傳多樣性

Table 4 Genetic diversity of P.mume based on four cpDNA fragments

居群樣本數(shù)單倍型Hdπ/(×10?3)居群樣本數(shù)單倍型Hdπ/(×10?3) AHSW 4H500SCRZ10H200 AHHS 4H700SCCJ10H2(5)、H3(2)、H4(3)0.7780.76 FJSD10H1(4)、H2(2)、H9(2)、H10(2)0.8000.65SCDZ10H13(6)、H14(2)、H15(2)0.6220.46 FJXY 5H800SCLD 6H1100 GZZJ10H200SCHL10H1400 GZNZ10H1(3)、H2(7)0.4670.24SCXY10H200 GZTJ10H19(7)、H20(3)0.4670.12SCMZ10H1(2)、H2(5)、H5(3)0.6890.54 GZML10H200SCSH10H200 GZBK10H200YNTC 5H1200 GZWZ10H1(4)、H2(6)0.5330.28YNDY 5H1800 JXBS 4H600XZLZ10H16(8)、H17(2)00 JXML10H500

2.4 群體歷史動(dòng)態(tài)

利用DnaSP5.10進(jìn)行中性檢驗(yàn)，Tajinma’s D的值為?0.142（＞0.10），未達(dá)到顯著水平，符合中性進(jìn)化，進(jìn)化方式為負(fù)向選擇。Fu and Li’D*和Fu and Li’F*的值分別為1.240（＞0.10）和0.836（＞0.10），均未達(dá)到顯著水平，符合中性進(jìn)化。結(jié)果表明，藥用梅野生群體在物種水平上未經(jīng)歷快速擴(kuò)張等歷史事件。

藥用梅野生群體的失配分布分析結(jié)果見圖3，可以看出，實(shí)際觀測(cè)的擴(kuò)張曲線呈現(xiàn)雙峰，而期望擴(kuò)張曲線則比較平滑，實(shí)際觀測(cè)曲線與期望觀測(cè)曲線不吻合，表明野生藥用梅群體歷史是穩(wěn)定的，沒有經(jīng)歷明顯的過快速擴(kuò)張事件，分析結(jié)果與中性檢驗(yàn)結(jié)果一致。

圖3 藥用梅群體的失配分布曲線

3 討論

3.1 遺傳多樣性

本研究利用葉綠體基因A-H、16、S-G、1b聯(lián)合序列對(duì)藥用梅野生群體進(jìn)行分析，結(jié)果顯示，藥用梅總的d為0.790，π為0.001 71，高于平均葉綠體DNA單倍型遺傳多樣性（0.670）[22]，表明整個(gè)野生群體遺傳多樣性較高，群體間的遺傳變異高于群體內(nèi)。藥用梅野生群體在我國(guó)主要分布在長(zhǎng)江流域[3]，在四川西南部、云南西北部和西藏東部交界的地區(qū)呈連續(xù)型分布，浙江、安徽、江西、福建等地有零散分布。20個(gè)單倍型中，H2是分布最廣泛的共享單倍型，H1、H15位于網(wǎng)絡(luò)圖中心位置，福建、四川地區(qū)野生群體的單倍型類型比較豐富，貴州、福建、四川地區(qū)是藥用梅的原始單倍型?；?個(gè)cpDNA片段構(gòu)建的單倍型網(wǎng)絡(luò)圖和系統(tǒng)進(jìn)化樹結(jié)果基本一致，將23個(gè)野生群體劃分為2支，在云貴高原周緣地區(qū)呈連續(xù)型分布的野生群體分在一組，貴州和華東地區(qū)的安徽、江西、福建分在一組，四川群體分散在2個(gè)分支中，單倍型分別滲透到云貴高原和華東地區(qū)。

根據(jù)溯祖理論，位于網(wǎng)狀進(jìn)化樹中心的單倍型和分支是最原始的，是其它單倍型和亞分支的祖先，其集中分布的地區(qū)是該物種的起源中心。結(jié)合單倍型構(gòu)建的進(jìn)化樹結(jié)果，推斷出四川、福建作為藥用梅主產(chǎn)區(qū)，也是其起源和多樣化中心。

3.2 遺傳結(jié)構(gòu)

葉綠體聯(lián)合片段分析結(jié)果顯示，群體間的遺傳分化系數(shù)st為0.673，st為0.891，st值大于st值，表明藥用梅的野生群體間存在譜系地理結(jié)構(gòu)，地理區(qū)域近的群體間的親緣關(guān)系也較近。

被子植物群體間的基因交流主要依靠花粉和種子基因流進(jìn)行傳播，核基因可以通過花粉和種子兩種方式進(jìn)行交流[23-24]，對(duì)于母系遺傳的葉綠體基因來說只能依靠種子流來傳播。群體間的基因流N為0.03，群體間的遺傳變異（89.10%）遠(yuǎn)大于群體內(nèi)（10.90%）。參考其它木本植物，中國(guó)櫻桃的野生群體間的分化低于群體內(nèi)[25]，槲樹、楊梅、杜梨等的野生群體間的分化高于群體內(nèi)[23,26-27]，推測(cè)藥用梅群體間分化遠(yuǎn)遠(yuǎn)大于群體內(nèi)，表明藥用梅群體間的基因流受阻，可能是因?yàn)樗幱妹返墓麑?shí)和種子較大，不適宜進(jìn)行遠(yuǎn)距離擴(kuò)散造成的。中性檢驗(yàn)和失配分析，常用于推測(cè)群體經(jīng)歷的歷史事件及群體擴(kuò)張時(shí)間[28-29]。在本研究中，中性檢驗(yàn)和失配分析結(jié)果一致，表明藥用梅野生群體歷史是基本穩(wěn)定的，沒有經(jīng)歷明顯的快速擴(kuò)張事件。

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

[1] 張華月, 李琦, 付曉伶. 烏梅化學(xué)成分及藥理作用研究進(jìn)展 [J]. 上海中醫(yī)藥雜志, 2017, 51(S1): 296-300.

[2] 楊琳, 陳金文, 王興海, 等. 烏梅本草考證 [J]. 中國(guó)現(xiàn)代中藥, 2019, 21(2): 247-251.

[3] 陳俊愉. 中國(guó)梅花 [M]. ?？? 海南出版社, 1996: 236.

[4] Ou J M, Wang R, Li X L,. Comparative analysis of free amino acids and nucleosides in different varieties of mume fructus based on simultaneous determination and multivariate statistical analyses [J]., 2020, 2020: 4767605.

[5] 歐金梅, 王瑞, 程慶兵, 等. ICP-MS法測(cè)定不同產(chǎn)地烏梅無機(jī)元素含量 [J]. 中草藥, 2020, 51(2): 482-489.

[6] 袁慶軍, 黃璐琦, 郭蘭萍, 等. 展望分子譜系地理學(xué)在道地藥材研究中的應(yīng)用 [J]. 中國(guó)中藥雜志, 2009, 34(16): 2007-2011.

[7] 李菁, 向婷婷, 許成, 等. 分子譜系地理學(xué)在藥用植物基因水平的研究進(jìn)展 [J]. 中草藥, 2015, 46(8): 1243-1246.

[8] Yuan Q J, Zhang Z Y, Peng H,. Chloroplast phylogeography of(Dipentodontaceae) in southwest China and northern Vietnam [J]., 2008, 17(4): 1054-1065.

[9] 拉瓊, 張勇群, 扎西次仁. 青藏高原特有植物分子譜系地理學(xué)研究現(xiàn)狀 [J]. 西藏大學(xué)學(xué)報(bào): 自然科學(xué)版, 2013, 28(1): 12-15.

[10] Zhang Y Z, Ma J, Yang B X,. The complete chloroplast genome sequence ofvar.(Taxaceae): Loss of an inverted repeat region and comparative analysis with related species [J]., 2014, 540(2): 201-209.

[11] 程蓓蓓. 中國(guó)紅豆杉屬分子譜系地理學(xué)與遺傳多樣性研究 [D]. 北京: 中國(guó)林業(yè)科學(xué)研究院, 2016.

[12] 姜丹. 黃岑道地性的遺傳和化學(xué)物質(zhì)基礎(chǔ)研究 [D]. 北京: 北京中醫(yī)藥大學(xué), 2018.

[13] 張玉君. 基于葉綠體分子譜系地理學(xué)的山藥的系統(tǒng)進(jìn)化和品種鑒定研究 [D]. 北京: 北京中醫(yī)藥大學(xué), 2018.

[14] Xu J J, Zang F Q, Wu Q C,. Analysis of the genetic diversity and molecular phylogeography of the endangered wild rose () in China based on chloroplast genes [J]., 2021, 28: e01653.

[15] Tóth E G, K?b?lkuti Z A, Pedryc A,. Evolutionary history and phylogeography of Scots pine (L.) in Europe based on molecular markers [J]., 2017, 28(4): 637-651.

[16] 石甜, 莫忠妹, 吳敏, 等. 藥食同源植物薤白的譜系地理學(xué)研究[J/OL].植物研究, [2022-01-03]. http://kns.cnki.net/kcms/detail/23.1480.s.20211130.1419.002.html.

[17] Librado P, Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data [J]., 2009, 25(11): 1451-1452.

[18] Leigh J W, Bryant D. Popart: Full-feature software for haplotype network construction [J]., 2015, 6(9): 1110-1116.

[19] Schneider S, Roessli D, Excoffier L. Arlequin: a software for population genetics data analysis [J]., 2000, 2: 2496-2497.

[20] Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism [J]., 1989, 123(3): 585-595.

[21] Fu X Y. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection [J]., 1997, 147(2): 915-925.

[22] Petit R J, Hampe A. Some evolutionary consequences of being a tree [J]., 2006, 37: 187-214.

[23] 劉路賢. 栽培楊梅譜系地理、馴化起源及品種(系)間關(guān)系的研究 [D]. 杭州: 浙江大學(xué), 2016.

[24] 徐興興. 矮牡丹的遺傳多樣性及栽培牡丹起源研究 [D]. 北京: 北京林業(yè)大學(xué), 2018.

[25] Chen T, Wang X R, Tang H R,. Genetic diversity and population structure of Chinese Cherry revealed by chloroplast DNAQ-16 intergenic spacers variation [J]., 2013, 60(6): 1859-1871.

[26] Yang J, Vázquez L, Chen X D,. Development of chloroplast and nuclear DNA markers for Chinese oaks (subgenus) and assessment of their utility as DNA barcodes [J]., 2017, 8: 816.

[27] Zong Y, Sun P, Liu J,. Chloroplast DNA-based genetic diversity and phylogeography ofbetulaefolia (Rosaceae) in Northern China [J]., 2014, 10(3): 739-749.

[28] Slatkin M, Hudson R R. Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations [J]., 1991, 129(2): 555-562.

[29] Roger AR, Harpending H. Population growth makes waves in the distribution of pairwise genetic differences [J]., 1992, 9(3): 552-569.

Analysis of population genetics of medicinal plum based on chloroplast gene

OU Jin-mei1, 2, 3, YANG Ya-tian1, QIAN Cheng-cheng1, WU Rui1, WANG Rui4, HUANG Lu-qi2

1. Anhui University of Chinese Medicine, Hefei 230012, China 2. National Resource Center for Chinese Meteria Medica, China Academy of Chinese Medical Sciences, Beijing 100700, China 3. Institute of Chinese Medicine Resources Protection and Development, Anhui Academy of Chinese Medicine, Hefei 230012, China 4. Wuhu Institute of Technical, Wuhu 241006, China

To study the population genetics of, reveal its genetic variation, haplotype geographical distribution pattern, and speculate its origin and diversification center.Chloroplast gene fragments (A-H,16,S-G,1b) were used to sequence the wild populations of medicinal plum from 23 regions in China, to carry out population genetic diversity, genetic structure and historical dynamic analysis.Twenty haplotypes were detected inpopulations, with a haplotype diversity (d) of 0.790 and a nucleotide diversity (π) of 0.001 71. AMOVA analysis showed that the genetic variation at the cpDNA level mainly occurred among populations, indicating that the species had an obvious phylogeographical structure inoccurred mainly between populations, which indicated that the species had a distinct distinct genealogical geographic structure and that interpopulation the gene flow was mainly dominated by pollen flow. Neutral test and mismatch distribution analysis showed that the population history was largely stable and did not experience significant rapid expansion event. Among the haplotypes, H2 was the most widely distributed shared haplotype, while H1 and H15 were located in the center of the haplotype network diagram.has obvious pedigree geographical structure with high population heritability, and the province of Sichuan and Fujian, as the main producing areas of, are also the centers of their origin and diversification.

Sieb. et Zucc.; chloroplast gene; population genetics; phylogeography;A-H;16;S-G;1b

R282.12

A

0253 - 2670(2022)17 - 5469 - 07

10.7501/j.issn.0253-2670.2022.17.023

2021-12-06

名貴中藥資源可持續(xù)利用能力建設(shè)項(xiàng)目（2060302-2201-23）；中央財(cái)政林業(yè)科技推廣示范資金項(xiàng)目（Z175070050002）

歐金梅（1980—），女，博士，副教授，研究方向?yàn)橹兴庂Y源與質(zhì)量評(píng)價(jià)。E-mail: ojm@ahtcm.edu.cn

黃璐琦，研究員，中國(guó)工程院院士。E-mail: huangluqi01@126.com

[責(zé)任編輯 時(shí)圣明]