蒙古黃芪全基因組SSR分子標記開發(fā)與居群遺傳多樣性分析
2025-04-07 00:00:00高慧霞梁云慧姚妙卓張鵬飛劉亞令
草地學報 2025年3期
摘要：為了解析蒙古黃芪［Astragalus membranaceus"（Fisch.）"Bge. var. mongholicus"（Bge.）"Hsiao］的遺傳多樣性及遺傳結(jié)構(gòu)，本研究基于蒙古黃芪全基因組數(shù)據(jù)進行了SSR分子標記的開發(fā)及引物的設(shè)計。結(jié)果顯示全基因組內(nèi)共檢索到1349個SSR位點，從單核苷酸重復至六核苷酸重復共有43種重復類型，其中二、三核苷酸重復最多，分別為6.45%和85.69%。通過對SSR分子標記有效性和多態(tài)性的篩選，得到15對引物多態(tài)信息含量變化范圍為0.076～0.462，具有良好的多態(tài)性。利用篩選出的引物對來自7個不同居群共120份蒙古黃芪樣本進行遺傳多樣性的分析，得到Shannon’s信息指數(shù)平均值為0.535，Nei’s基因多樣性指數(shù)平均值為0.422，渾源居群蒙古黃芪的遺傳多樣性最高，五寨居群蒙古黃芪的遺傳多樣性最低；居群間分化系數(shù)平均值為0.164，與AMOVA結(jié)果一致，蒙古黃芪的遺傳分化主要發(fā)生在居群內(nèi)；UPGMA聚類分析將不同居群的蒙古黃芪分別聚為一支，具有一定的地域性。通過蒙古黃芪遺傳多樣性分析結(jié)果驗證了本研究開發(fā)的SSR分子標記可用于核心種質(zhì)構(gòu)建、良種選育與不同物種分子鑒定等研究。
關(guān)鍵詞：蒙古黃芪；基因組；SSR；遺傳多樣性
中圖分類號：S567.239 """""""文獻標識碼：A """""""文章編號：1007-0435（2025）03-0696-11
Development of Genome-Wide SSR Molecular Markers and Analysis of Population Genetic Diversity in Astragalus mongolicus
GAO Hui-xia1，"LIANG Yun-hui1，"YAO Miao-zhuo1，"ZHANG Peng-fei2，"LIU Ya-ling1，3*
（1.College of Life Science，"Shanxi Agricultural University，"Taigu，"Shanxi Province 030801，"China；"2.College of Horticulture，"Shanxi Agricultural University，"Taigu，"Shanxi Province 030801，"China；"3.Shanxi Key Lab. for Modernization of TCVM，"College of Veterinary Medicine，"Shanxi Agricultural University，"Taigu，"Shanxi Province 030801，"China）
Abstract：To elucidate the genetic diversity and structure of Astragalus membranaceus"（Fisch.）"Bge. var. mongholicus"（Bge.）"Hsiao，"this study developed SSR （Simple Sequence Repeat）"molecular markers and designed primers leveraging the whole-genome sequence data. The outcomes revealed a total of 1349 SSR loci within the genome，"encompassing 43 distinct repeat types spanning from mononucleotide to hexanucleotide，"with dinucleotide （6.45%）"and trinucleotide （85.69%）"repeats being the most prevalent. After rigorous screening for efficacy and polymorphism，"15 primer pairs were selected，"exhibiting a polymorphism information content range of 0.076 to 0.462，"indicating robust polymorphism. These optimized primers were then utilized to analyze the genetic diversity of 120 A. mongholicus"samples from seven distinct populations. The average Shannon's diversity index was 0.535，"and the average Nei's gene diversity was 0.422. Notably，"the Hunyuan population demonstrated the highest genetic diversity，"whereas the Wuzhai population exhibited the lowest. The average fixation index among populations was 0.164，"concurring with AMOVA results，"suggesting that the primary genetic differentiation occurred within populations. UPGMA cluster analysis segregated the A. mongholicus"samples from different populations into distinct clusters，"reflecting geographical patterns. This study validated the utility of the developed SSR markers for applications such as core germplasm development，"elite cultivar breeding，"and molecular identification of related species.
Key words：A. mongholicus；Genome；SSR；Genetic diversity
黃芪為豆科黃芪屬植物蒙古黃芪［Astragalus membranaceus"（Fisch.）"Bge. var. mongholicus"（Bge.）"Hsiao］和膜莢黃芪［Astragalus membranaceus"（Fisch.）"Bge.］的干燥根［1］，具有補氣固表、利水消腫等功效，已有2000多年的藥用歷史［2］。蒙古黃芪為主流品種，除作藥用還被廣泛應用于食品保健方面［3］。隨著市場需求量的不斷增大，各地野生蒙古黃芪資源被無序采挖，加之生長環(huán)境不斷遭到人為破壞導致其資源利用量持續(xù)減少［4］，現(xiàn)如今僅有少部分野生資源分布在山西、陜西、河北等地區(qū)［5-6］。為了緩解野生資源利用量的壓力，人工栽培蒙古黃芪的面積及產(chǎn)量逐年增加，但市場調(diào)查發(fā)現(xiàn)，在栽培過程中存在種源混雜嚴重的現(xiàn)象［7］，加重了蒙古黃芪種間基因交流或遺傳漂變的發(fā)生，造成原始種質(zhì)“道地性”的喪失，影響到優(yōu)良種質(zhì)基因的篩選和育種工作的順利進行。因此，研究蒙古黃芪遺傳多樣性與遺傳結(jié)構(gòu)，為其資源保護與可持續(xù)利用提供科學依據(jù)。
簡單重復序列（Simple sequence repeats，"SSR）是一段由單核苷酸至六核苷酸基元串聯(lián)重復組成的DNA序列，長度通常為200 bp之內(nèi)，廣泛存在于真核生物基因組中。SSR序列具有共顯性、高多態(tài)性、易于操作等優(yōu)點，已成功應用于苦蕎［Fagopyrm tataricum （L.）"Gaertn.］［8］、防風［Saposhnikovia divaricatanbsp;（Turcz.）"Schischk.］［9］、紫蘇［Perilla frutescensnbsp;（L.）"Britt.］［10］等多種藥用植物遺傳多樣性分析、親緣關(guān)系鑒定以及指紋圖譜的構(gòu)建等方面［11-12］。在蒙古黃芪的SSR開發(fā)與應用中，賀潤麗等［13］利用蒙古黃芪轉(zhuǎn)錄組數(shù)據(jù)開發(fā)了19對SSR分子標記，可以準確反映20份蒙古黃芪和膜莢黃芪樣本間的差異；劉亞令等［14］利用同科植物甘草的10對SSR引物對380份蒙古黃芪和膜莢黃芪的遺傳多樣性和遺傳分化進行了分析，并發(fā)現(xiàn)蒙古黃芪的遺傳多樣性高于膜莢黃芪。相較于使用轉(zhuǎn)錄組數(shù)據(jù)開發(fā)SSR分子標記，全基因組序列覆蓋范圍更廣，不僅可以在基因表達區(qū)進行開發(fā)，還可以在非編碼區(qū)域開發(fā)得到更具特異性和穩(wěn)定性的標記，更好地反映物種的遺傳變異［15-16］。因此，本試驗基于蒙古黃芪全基因組數(shù)據(jù)進行SSR位點的開發(fā)與引物設(shè)計，并對來自7個不同居群共120份蒙古黃芪樣本進行遺傳多樣性的分析，為后續(xù)蒙古黃芪分布格局的解析、種質(zhì)資源的保護及分子輔助育種等研究提供有效的分子標記工具。
1 材料與方法
1.1 試驗材料
本試驗蒙古黃芪樣本采自山西、寧夏、內(nèi)蒙古、甘肅、河北與陜西6個省份，共7個居群120份樣本（表1）。將采集的新鮮嫩葉于干燥硅膠中保存，帶回實驗室進行蒙古黃芪基因組DNA的提取。
1.2 DNA的提取
將干燥的蒙古黃芪葉片放入高通量組織研磨儀中進行研磨，采用改良的CTAB法［17］提取基因組DNA，利用核酸蛋白儀檢測DNA樣本的濃度，測定A260/A280的值。利用1%瓊脂糖凝膠電泳檢測DNA樣本的質(zhì)量。
1.3 基因組SSR位點開發(fā)及引物設(shè)計
蒙古黃芪樣本質(zhì)檢合格后，首先利用磁珠富集、純化大片段DNA（大于15 k）。隨后，對片段化的DNA進行損傷及末端修復，純化后使用SQK-LSK109試劑盒中進行連接反應，完成文庫構(gòu)建。使用Qubit對構(gòu)建好的DNA文庫進行定量檢測。接著，將純化后的文庫使用PromethION P48測序儀進行測序。然后使用Oxford Nanopore GUPPY（version 4.0.2）軟件將高通量測序得到的原始圖像數(shù)據(jù)文件經(jīng)堿基識別后轉(zhuǎn)化為fast5原始數(shù)據(jù)，結(jié)果以FASTQ文件格式存儲。隨后，使用Oxford Nanopore GUPPY（version 3.0）軟件對原始數(shù)據(jù)進行質(zhì)量評估，過濾掉測序質(zhì)量值Qlt;7的低質(zhì)量reads，得到有效數(shù)據(jù)。最后，使用NextDenovo軟件進行組裝得到完整的蒙古黃芪基因組序列。
使用MISA（http：//pgrc.ipk-gatersleben.de/misa/misa.html）對蒙古黃芪基因組序列進行SSR位點的分析與檢索，參數(shù)設(shè)定為：單核苷酸重復10次，二核苷酸重復6次，三核苷酸至六核苷酸各重復5次，復合型SSR中2個簡單完整型SSR之間的距離小于100 bp。使用Primer 3.0（https：//Primer3.ut.ee）對SSR位點進行篩選并設(shè)計對應的引物，引物設(shè)計原則：GC含量在40%～60%之間，Tm值在50℃～60℃之間，上下游引物Tm差值不超過5℃，產(chǎn)物大小在500 bp之內(nèi)，引物長度在18～23 bp之間。隨機選取均勻分布于蒙古黃芪染色體并且產(chǎn)物大小在500 bp之內(nèi)的20對SSR引物送至生工生物工程（上海）股份有限公司合成。
1.4 SSR-PCR擴增
PCR擴增反應體系的總體積為10 μL，其中10×PCR Buffer（含Mg2+）1 μL，上下游引物各0.25 μL，2.5 mmol·L-1"dNTP 0.8 μL，5 U·μL-1"Taq酶0.25 μL，模板DNA 3 μL，ddH2O 4.45 μL。擴增的程序為：95℃預變性5 min；94℃變性40 s，55℃退火50 s，72℃延伸45 s，35個循環(huán)；72℃再延伸10 min。
擴增后的產(chǎn)物采用6%聚丙烯酰胺凝膠電泳進行檢測。
1.5 數(shù)據(jù)處理
采用人工方式進行讀帶，將電泳條帶轉(zhuǎn)化為0/1矩陣，使用PowerMarkerV3.25軟件計算SSR引物的多態(tài)性信息含量（PIC）、等位基因數(shù)（Na）、有效等位基因數(shù)（Ne）、Shannon’s信息指數(shù)（I）、期望雜合度（He）、觀察雜合度（Ho）。使用GenAlex 6.503軟件計算蒙古黃芪不同居群的等位基因數(shù)（Na）、有效等位基因數(shù)（Ne）、Shannon’s信息指數(shù)（I）、Nei’s基因多樣性（H）、期望雜合度（He）、觀察雜合度（Ho）、居群內(nèi)近交系數(shù)（Fis）、總居群近交系數(shù)（Fit）、居群間遺傳分化系數(shù)（Fst）、基因流（Nm）等遺傳多樣性參數(shù)；使用GenAlex 6.503軟件蒙古黃芪居群進行分子方差檢驗（AMOVA）；使用STRUCTURE軟件進行遺傳結(jié)構(gòu)的分析，將居群分組數(shù)設(shè)定為2～7，每個K值運算重復10次，將結(jié)果在StructureSelector（https：//lmme.ac.cn/StructureSelector/）中可視化進行蒙古黃芪最佳居群分組數(shù)的預測；使用MEGA 7.0［18］軟件中的UPGMA（非加權(quán)算術(shù)平均值聚類法）進行聚類分析，繪制居群聚類圖。
2 結(jié)果與分析
2.1 蒙古黃芪全基因組中SSR位點的分布
2.1.1 基因組中SSR在染色體上的分布特點 本試驗對蒙古黃芪全基因組數(shù)據(jù)進行SSR位點的檢索（表2），得到19 822條序列共存在1349個SSR位點。SSR位點發(fā)生頻率為5.83%，出現(xiàn)頻率為6.81%，平均分布距離為15 421 bp。其中有1011條序列僅包含1個SSR位點，144條序列包含1個以上SSR位點，以復合物形式存在的SSR位點有97個。蒙古黃芪全基因組8條染色體序列總長度約為1.40 Gb，其中Chr3包含SSR位點最多為318個，Chr6包含SSR位點最少為97個（圖1）。
2.1.2 SSR位點的重復類型與重復次數(shù) 蒙古黃芪基因組中共得到六種不同的SSR位點重復類型（圖2），主要集中在二、三核苷酸重復類型上。其中單核苷酸重復為77個，占全部SSR比例為5.71%；二核苷酸重復為87個，所占比例為6.45%；三核苷酸重復最多為1156個，所占比例為85.69%；四核苷酸重復最少為2個，所占比例為0.15%；五核苷酸重復為5個，所占比例為0.37%；六核苷酸重復為22個，所占比例為1.63%。
所有SSR類型重復次數(shù)中，5次重復的最多有634個，所占比例為47.00%；其次6次重復為333個，所占比例為24.68%；14次重復的最少有7個，所占比例為0.52%；≥15次重復為47個，所占比例為3.48%?？傮w上，當重復次數(shù)為5時，SSR位點數(shù)量為最多，隨著重復次數(shù)的增大，SSR位點數(shù)量呈現(xiàn)遞減趨勢。在單核苷酸重復中，10次重復的SSR數(shù)量最多（44個），所占比例為57.14%；二核苷酸重復中，6次重復的數(shù)量最多（41個），所占比例為47.13%；三核苷酸重復中，5次重復的數(shù)量最多（613個），所占比例為53.03%；四核苷酸重復中，僅有6次重復和11次重復，數(shù)量均為1；五、六核苷酸重復均為5次重復最多，數(shù)量和所占比例分別為3，18和60%，81.82%。
圖1 SSR位點在蒙古黃芪染色體上的分布情況
Fig.1 Distribution of SSR loci on chromosome of A. mongholicus
圖2 SSR位點重復次數(shù)
Fig.2 Repeat times of SSR loci
2.1.3 SSR位點基元重復類型及數(shù)量 對蒙古黃芪基因組SSR核苷酸重復基元的類型進行分析，結(jié)果如圖3所示。在1349個SSR位點共有43種重復基元類型，從單核苷酸重復基元至六核苷酸重復分別為77，87，1156，2，5和22。單核苷酸重復類型為A/T（68個，出現(xiàn)頻率為5.04%），C/G（9個，0.67%）（圖3A）；二核苷酸重復類型為AC/GT（10個，0.74%），AG/CT（58個，4.30%），AT/AT（19個，1.41%）（圖3B）；三核苷酸重復類型有10類（圖3C），其中AAC/GTT（235個，17.42%），AAG/CTT（262個，19.42%）重復類型為最多；四核苷酸重復類型為ACAT/ATGT（2個，0.22%）；五核苷酸重復類型有5類（圖3D），六核苷酸重復類型有22類（圖3E），重復數(shù)目均為1個，出現(xiàn)頻率均為0.07%。
2.2 蒙古黃芪基因組DNA提取
本試驗利用CATB法對7個不同居群共120份蒙古黃芪樣本的基因組DNA進行了提取，使用1%瓊脂糖凝膠電泳檢測（圖4），結(jié)果顯示樣品條帶清晰明亮，符合SSR-PCR擴增體系對DNA質(zhì)量的要求。
2.3 蒙古黃芪SSR引物多態(tài)性檢測與篩選
2.3.1 SSR引物多態(tài)性篩選 本實驗將合成的20對SSR引物對120份蒙古黃芪基因組DNA進行SSR-PCR擴增，通過6%聚丙烯酰胺凝膠電泳對引物進行多態(tài)性的篩選（圖5）。結(jié)果得到20對引物中有15對引物可擴增出目標條帶，且條帶清晰明亮，具有較好的重復性和豐富的多態(tài)性，可用于后續(xù)蒙古黃芪植物遺傳多樣性試驗中（表3）。剩余5對引物為不能擴增出條帶或條帶單一不具有多態(tài)性。
2.3.2 SSR引物多態(tài)性信息 本實驗利用PowerMarkerV3.25軟件計算SSR引物的Na，Ne，I，Ho，He，PIC（表4）。結(jié)果得到15對SSR引物序列包含多態(tài)性位點的Na變化范圍為1.714～2.000，平均值為1.933個；Ne變化范圍為1.099～1.923，平均值為1.655；I變化范圍為0.178～0.670，平均值為0.535；Ho變化范圍為0.059～0.789，平均值為0.371；He變化范圍為0.088～0.478，平均值為0.368；PIC變化范圍為0.176～0.498，平均值為0.405，有14對引物呈現(xiàn)中度多態(tài)性，1對引物呈現(xiàn)低度多態(tài)性。表明本試驗篩選得到的15對引物具有良好的多態(tài)性。
2.4 蒙古黃芪遺傳多樣性分析
2.4.1 蒙古黃芪遺傳多樣性分析 使用GenAlex 6.503對7個居群120份蒙古黃芪樣本進行遺傳多樣性的分析（表5）。不同居群的蒙古黃芪的Na范圍為1.733～2.000，均值為1.933；Ne范圍為1.564～1.739，均值為1.655；I范圍為0.485～0.573，均值為0.535；H范圍為0.384～0.456，均值為0.422；Ho范圍為0.318～0.467，均值為0.371；He范圍為0.339～0.418，均值為0.381。由此可得，蒙古黃芪遺傳多樣性水平最高的居群為渾源，最低的居群為五寨。其中，野生蒙古黃芪不同居群I的均值為0.559，栽培蒙古黃芪不同居群I的均值為0.526，表明野生蒙古黃芪遺傳多樣性高于栽培居群。
2.4.2 蒙古黃芪遺傳分化分析 使用GenAlex 6.503對7個居群120份蒙古黃芪樣本進行遺傳分化分析。結(jié)果如表6所示。Fis的范圍為-0.674～0.805，均值為0.008，表明蒙古黃芪居群雜合子不足，特別是位點HQ-19表現(xiàn)雜合子不足最明顯；Fit的范圍為-0.624～0.880，均值為0.148，同樣表明蒙古黃芪居群雜合子不足；Fst的范圍為0.015～0.391，均值為0.164；Nm的范圍為0.389～16.362，均值為3.880，表明蒙古黃芪居群間遺傳分化較小，基因流水平較高。
2.4.3 蒙古黃芪AMOVA分析 使用GenAlex 6.503軟件對7個居群的蒙古黃芪進行AMOVA分析，結(jié)果表明（表7），蒙古黃芪居群間的遺傳變異率為13%，居群內(nèi)遺傳變異率為87%。居群間遺傳變異率小于居群內(nèi)遺傳變異率，說明蒙古黃芪的遺傳變異主要發(fā)生在居群內(nèi)。
2.4.4 蒙古黃芪遺傳結(jié)構(gòu)分析 使用STRUCTURE軟件基于貝葉斯模型對7個居群蒙古黃芪進行遺傳結(jié)構(gòu)分析，當ΔK值達到最大時對應的K值為模擬的最佳居群數(shù)。將分組數(shù)設(shè)定為2～7，每個K值運算重復10次，得到K值與ΔK值的關(guān)系曲線圖（圖6），結(jié)果表明當K=4時，ΔK=2.820達到最大值，此時可以將7個蒙古黃芪居群分為4組（圖7），每個居群都有四種顏色的參與，說明蒙古黃芪基因具有復雜的遺傳背景，每個居群樣本都顯示有較為強烈的基因漸滲。
2.4.5 蒙古黃芪不同居群聚類分析 使用GenAlex 6.503計算蒙古黃芪不同居群間的遺傳距離，利用MEGA 7.0構(gòu)建UPGMA系統(tǒng)發(fā)育樹。結(jié)果表明（圖8），7個不同居群的蒙古黃芪樣本優(yōu)先按照地理進行聚類，不同居群的樣本之間出現(xiàn)少量混雜。與STRUCTURE結(jié)果一致，7個居群分為四組，第一組為隴西居群和子洲居群聚為一支；第二組為渾源居群與隆德居群聚為一支；第三組為主要為內(nèi)蒙古居群先與承德居群聚為一支，再與五寨居群聚為一支，其中有7個隴西樣本、5個子洲樣本、4個渾源樣本和1個隆德樣本混雜；第四組為內(nèi)蒙古居群中2個樣本聚為一支。
3 討論
本試驗基于蒙古黃芪全基因組測序數(shù)據(jù)進行SSR位點的檢索，共得到1349個SSR位點分布在1155條序列中，發(fā)生頻率為5.83%，出現(xiàn)頻率為6.81%，平均分布距離為15 421 bp。蒙古黃芪基因組SSR數(shù)量低于丹參（Salvia miltiorrhiza Bunge）［19］、罌粟（Papaver somniferum"L.）［20］等植物，可能是由于物種的基因組大小和結(jié)構(gòu)以及在SSR開發(fā)過程中設(shè)定的篩選標準不同。蒙古黃芪基因組大小為1.40 Gb，設(shè)定SSR篩選標準為1～6核苷酸重復依次為10，6，5，5，5，5，而丹參基因組大小為2.40 Gb，設(shè)定2～6核苷酸重復依次為6，4，3，3，2，罌粟基因組大小為2.22 Gb，設(shè)定1～6核苷酸重復依次為10，6，4，3，3，3。另外，物種獨特的遺傳背景導致基因組發(fā)生不同程度的突變也會對包含的SSR位點信息造成不同程度的改變。
本試驗通過對蒙古黃芪基因組中SSR位點信息進行分析，發(fā)現(xiàn)SSR位點重復類型較為豐富，1～6核苷酸均有分布，二、三核苷酸重復最多，分別占總SSR位點的6.75%，85.69%，表明蒙古黃芪可能起源較晚但進化水平較高，與綠豆［Vigna radiata"（L.）"Wilczek］［21］、黃芩（Scutellaria baicalensis Georgi）［22］、三七［Panax notoginseng （Burk.）"F. H. Chen］［23］等物種相似。二、三核苷酸重復頻率相較四、五、六核苷酸重復頻率高，可能是由于二、三核苷酸形成相對容易，結(jié)構(gòu)穩(wěn)定且存在突變偏好性，從而導致重復積累的數(shù)量較多。SSR核苷酸重復基元結(jié)構(gòu)共有43種，其中，優(yōu)勢基元主要為A/T，AG/CT，AAG/CTT，與蒙古黃芪轉(zhuǎn)錄組數(shù)據(jù)中SSR堿基偏好性相同［24］，也與楊梅（Myrica rubar"Sieb. et Zucc.）［25］、甘草（Glycyrrhiza uralensis"Fisch.）［26］、天麻（Gastrodia elata"Bl.）［27］等物種的結(jié)果一致，均以A/T為主要重復單元，堿基重復類型偏好性可能由以下因素造成：一是基因組DNA中CpG島的甲基化，通過脫氨基作用，甲基化的胞嘧啶C很可能轉(zhuǎn)變?yōu)樾叵汆奏［28］；二是G/C堿基的穩(wěn)定性，難以解鏈導致A/T堿基偏好性；三是自然選擇、遺傳背景、地理分布等因素，這些因素都有可能對物種基因組的結(jié)構(gòu)產(chǎn)生一定的影響，從而導致特定的SSR重復基元類型占據(jù)主導地位［29-30］。
利用篩選出的15對SSR引物對7個不同居群共120份蒙古黃芪樣本進行遺傳多樣性分析，得到蒙古黃芪的I，h高于巴戟天（Morinda officinalis F. C."How）［31］、廣金錢草［Desmodium styraci folium （Osb.）"Merr.］［32］、紫花苜蓿（Medicago sativa L.）［33］等藥用植物，表明蒙古黃芪的遺傳多樣性較為豐富。其中，野生蒙古黃芪居群的平均遺傳多樣性高于栽培居群，說明野生蒙古黃芪具有更豐富的遺傳變異類型和更強的適應環(huán)境變化能力［34］。另外，子洲地區(qū)的野生蒙古黃芪居群遺傳多樣性略低于隆德和內(nèi)蒙古地區(qū)的栽培居群的遺傳多樣性，可能是因為隆德和內(nèi)蒙古地區(qū)在栽培過程中種源復雜導致其遺傳多樣性略高［35］。蒙古黃芪的Fst為0.164，居群間分化較小。AMOVA分析同樣表明蒙古黃芪居群間遺傳分化主要來源于居群內(nèi)，Nm為3.880，與劉亞令［36］、王敖［37］等研究一致，即蒙古黃芪居群間遺傳分化水平較低，基因流水平較高，這可能與蒙古黃芪具有自交不親和的特性有關(guān)［38］。同時蒙古黃芪所處自然環(huán)境相對一致，在一定程度上也降低了居群間遺傳分化。對蒙古黃芪居群結(jié)構(gòu)分析發(fā)現(xiàn)，不同居群間的樣本出現(xiàn)少量混雜的現(xiàn)象，這可能是由于蒙古黃芪在栽培過程中種源的不確定及種苗引種栽培或運輸?shù)冗^程中促進基因在不同居群間流動，從而降低了地理因素對居群遺傳結(jié)構(gòu)的影響，導致不同居群的部分樣本相互混雜。
4 結(jié)論
本試驗基于蒙古黃芪基因組進行SSR分子標記位點的開發(fā)，共得到1349個SSR位點。選取20對SSR引物合成后進行篩選，得到15對多態(tài)性良好的SSR引物。利用篩選出的引物對來自7個不同居群共120份蒙古黃芪樣本進行遺傳多樣性及遺傳結(jié)構(gòu)的分析，結(jié)果得到蒙古黃芪具有較豐富的遺傳多樣性，其中渾源居群的遺傳多樣性最高，五寨居群的遺傳多樣性最低；AMOVA分析表明蒙古黃芪遺傳分化主要發(fā)生在居群內(nèi)；UPGMA聚類分析表明蒙古黃芪優(yōu)先以地理聚類。蒙古黃芪遺傳多樣性分析結(jié)果驗證了本研究開發(fā)的SSR分子標記可用于蒙古黃芪核心種質(zhì)構(gòu)建、良種選育與物種鑒定等研究。
參考文獻
［1］"國家藥典委員會. 中華人民共和國藥典：一部［M］. 北京：中國醫(yī)藥科技出版社：2020，315-316
［2］"汪娟. 黃芪的藥理作用研究進展［J］. 醫(yī)療裝備，2018，31（14）：202-203
［3］"王軍. 苦參藥用部位發(fā)育解剖暨黃芪道地藥材的年限鑒別及其次生代謝產(chǎn)物積累研究［D］. 合肥：安徽中醫(yī)藥大學，2013：17
［4］"李洪，王彧超，王瑞軍，等. 黃芪種質(zhì)資源研究與利用［J］. 安徽農(nóng)學通報，2020，26（20）：34，57
［5］"張?zhí)m濤，郭寶林，朱順昌，等. 黃芪種質(zhì)資源調(diào)查報告［J］. 中藥材，2006，29（8）：771-773
［6］"王劍. 黃芪種質(zhì)資源與優(yōu)良品種的選育研究［J］. 農(nóng)業(yè)災害研究，2022，12（7）：173-175
［7］"耿雅萍，馬艷珠，王芳，等. 基于psbA-trnH序列對藥用黃芪及其混偽品的鑒別［J］. 山西農(nóng)業(yè)科學，2020，48（5）：664-668
［8］"蔡齊宗，王佳蕊，陳慶富，等. 苦蕎全基因組SSR位點鑒定及分子標記開發(fā)［J］. 河南農(nóng)業(yè)大學學報，2022，56（3）：392-400
［9］"李爽，寇佩雯，許祎珂，等. 基于轉(zhuǎn)錄組的防風SSR分子標記開發(fā)及應用［J］. 中國中藥雜志，2024，49（3）：717-727
［10］"向依，潘金衛(wèi)，李慧琳，等. 基于SSR分子標記的紫蘇種質(zhì)資源遺傳多樣性及群體結(jié)構(gòu)分析［J］. 種子，2024，43（5）：40-47，2
［11］"馬驄毓，韓重陽，馬賽男，等. 基于SSR分子標記的13個白三葉（Trifolium repens"L.）品種指紋圖譜構(gòu)建［J］. 草地學報，2021，29（9）：1892-1899
［12］"陳志祥，羅小燕，李拴林，等. 基于SSR標記的木豆種質(zhì)資源遺傳多樣性與群體結(jié)構(gòu)分析［J］. 草地學報，2021，29（5）：904-911
［13］"賀潤麗，韓毅麗，王芳，等. 蒙古黃芪轉(zhuǎn)錄組SSR標記開發(fā)及遺傳多樣性研究［J］. 中國中藥雜志，2018，43（9）：1838-1843
［14］"劉亞令，耿雅萍，解瀟冬，等. 基于SSR分子標記的藥用黃芪遺傳多樣性與遺傳結(jié)構(gòu)分析［J］. 草地學報，2019，27（5）：1154-1162
［15］"GUPTA P K，BALYAN H S，SHARMA P C，et al. Microsatellites in plants："a new class of molecular markers［J］. Current Science，1996，70（1）：45-54
［16］"李建明，李洪杰，柴守誠，等. 普通小麥Genomic-SSR和EST-SSR分子標記遺傳差異［J］. 安徽農(nóng)業(yè)科學，2007，35（26）：8173-8175，8183
［17］"CAI C N，XIAO J H，CI X Q，et al. Genetic diversity of Horsfieldia tetratepala （Myristicaceae），"an endangered plant species with extremely small populations to China："implications for its conservation［J］. Plant Systematics and Evolution，2021，307（4）：50
［18］"KUMAR S，STECHER G，TAMURA K. MEGA7："Molecular evolutionary genetics analysis version 7.0 for bigger datasets［J］. Molecular Biology and Evolution，2016，33（7）：1870-1874
［19］"馮園園，郭林林，金花，等. 丹參基因組SSR的開發(fā)及其多態(tài)性分析［J］. 時珍國醫(yī)國藥，2019，30（4）：970-972
［20］"凌立貞，杜毛毛，張書東. 基于罌粟基因組的SSR位點鑒定及其特征分析［J］. 六盤水師范學院學報，2023，35（6）：99-105
［21］"葉衛(wèi)軍，陳圣男，楊勇，等. 綠豆SSR標記的開發(fā)及遺傳多樣性分析［J］. 作物學報，2019，45（8）：1176-1188
［22］"齊琳潔，龍平，蔣超，等. 黃芩基因組SSR分子標記的開發(fā)及遺傳多樣性分析［J］. 藥學學報，2015，50（4）：500-505
［23］"揭應碧，盧迎春，宋婉玲，等. 三七基因組SSR位點分析和多態(tài)性引物開發(fā)［J］. 熱帶作物學報，2018，39（7）：1351-1358
［24］"賀潤麗，韓毅麗，王芳，等. 蒙古黃芪轉(zhuǎn)錄組SSR標記開發(fā)及遺傳多樣性研究［J］. 中國中藥雜志，2018，43（09）：1838-1843
［25］"張淑文，梁森苗，鄭錫良，等. 楊梅基因組SSR引物的開發(fā)與應用［J］. 園藝學報，2019，46（1）：149-156
［26］"詹海仙，王穎莉，杜晨暉，等. 基于甘草全基因組序列的SSR分子標記開發(fā)［J］. 分子植物育種，2020，18（18）：6093-6100
［27］"周天華，丁家璽，田偉，等. 天麻基因組微衛(wèi)星特征分析與分子標記開發(fā)［J］. 西北植物學報，2017，37（9）：1728-1735
［28］"SCHORDERET D F，GARTLER S M. Analysis of CpG suppression in methylated and nonmethylated species［J］. Proceedings of the National Academy of Sciences，1992，89（3）：61-957
［29］"馬名川，劉龍龍，劉璋，等. 苦蕎全基因組SSR位點特征分析與分子標記開發(fā)［J］. 作物雜志，2021，（1）：38-46
［30］"柳娟，岳春雷，朱弘，等. 日本莢蒾轉(zhuǎn)錄組SSR位點分析及多態(tài)性SSR標記開發(fā)［J］. 分子植物育種，2022，20（21）：7183-7192
［31］"關(guān)婉，張丹純，丁曉霞，等. 廣東省巴戟天DNA條形碼及遺傳多樣性分析［J］. 世界科學技術(shù)-中醫(yī)藥現(xiàn)代化，2021，23（10）：3765-3771
［32］"何博，許勇，楊旻，等. 不同地區(qū)廣金錢草遺傳多樣性的ISSR分析［J］. 中國藥師，2016，19（5）：884-887
［33］"強海平，余國輝，劉海泉，等. 基于SSR標記的中美紫花苜蓿品種遺傳多樣性研究［J］. 中國農(nóng)業(yè)科學，2014，47（14）：2853-2862
［34］"王雪鳳，梁穎，劉德旺，等. 內(nèi)蒙古地區(qū)栽培及野生蒙古黃芪遺傳多樣性研究［J］. 中國草地學報，2018，40（1）：42-48
［35］"張茹. 藥用黃芪ISSR遺傳多樣性分析［D］. 太谷：山西農(nóng)業(yè)大學，2014：33-34
［36］"劉亞令，耿雅萍，解瀟冬，等. 基于SSR分子標記的藥用黃芪遺傳多樣性與遺傳結(jié)構(gòu)分析［J］. 草地學報，2019，27（05）：1154-1162
［37］"王敖. 蒙古黃耆和膜莢黃耆居群遺傳多樣性研究［D］. 北京：中央民族大學，2013：83
［38］"馮學金. 蒙古黃芪傳粉特性研究［J］. 中國現(xiàn)代中藥，2011，13（4）：24-26
（責任編輯""付宸）
引用格式：高慧霞， 梁云慧， 姚妙卓，"等.蒙古黃芪全基因組SSR分子標記開發(fā)與居群遺傳多樣性分析[J].草地學報，2025，33（3）：696-706
Citation：GAO Hui-xia， LIANG Yun-hui， YAO Miao-zhuo， et al.Development of Genome-Wide SSR Molecular Markers and Analysis of Population Genetic Diversity in"Astragalus mongolicus[J].Acta Agrestia Sinica，2025，33（3）：696-706
基金項目：國家自然科學基金（32070378）；山西農(nóng)業(yè)大學生物育種工程項目（YZGC136）；山西省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項資金（2021-11）資助
作者簡介：高慧霞（2001-），女，漢族，河南項城人，碩士研究生，主要從事藥用植物遺傳學研究，E-mail："ghxokng@163.com；*通信作者Author for correspondence，"E-mail："lylzpf@126.com
