基于轉(zhuǎn)錄組序列的萬壽菊SSR標(biāo)記開發(fā)及應(yīng)用

摘要：萬壽菊（Tagetes erecta）是國內(nèi)外廣泛栽植的一年生草本花卉，具有重要的觀賞、藥用和經(jīng)濟(jì)價值?；谵D(zhuǎn)錄組測序結(jié)果挖掘萬壽菊簡單重復(fù)序列（simple sequence repeat，SSR）位點(diǎn)信息，并對萬壽菊屬的2個品種、24個品系進(jìn)行遺傳多樣性分析和聚類分析。結(jié)果表明，轉(zhuǎn)錄組測序共挖掘到10 182 個SSR位點(diǎn)，出現(xiàn)頻率為16.07%，平均6.15 kb出現(xiàn)1個SSR位點(diǎn)。SSR位點(diǎn)共出現(xiàn)259種基序重復(fù)類型，其中三核苷酸是主要重復(fù)類型（48.51%），其次是二核苷酸（27.05%），二核苷酸、三核苷酸的優(yōu)勢重復(fù)基元分別為AG/CT（37.22%）、ATC/ATG（28.63%）。篩選出10個多態(tài)性位點(diǎn)用于萬壽菊資源的遺傳多樣性分析，平均期望雜合度（He）、多態(tài)性信息含量（PIC）分別為0.361、0.364。26份萬壽菊材料的平均Shannon息指數(shù)（I）、期望雜合度（He）分別為0.642、0.475，遺傳距離介于0.027～0.913之間，表明供試材料的遺傳多樣性較豐富。非加權(quán)配對算數(shù)平均法（UPGMA）聚類分析結(jié)果表明，26份萬壽菊材料在相似系數(shù)0.69處可分為兩大類，其中花色相同的萬壽菊屬品種（系）聚為一類。綜上，萬壽菊轉(zhuǎn)錄組中SSR的出現(xiàn)頻率較高、類型豐富、多態(tài)性高，為萬壽菊種質(zhì)資源鑒別與利用、遺傳圖譜構(gòu)建和分子標(biāo)記輔助選擇等提供了備用分子標(biāo)記，研究結(jié)果為今后萬壽菊核心種質(zhì)的選擇、新品種的選育奠定了理論基礎(chǔ)。

關(guān)鍵詞：萬壽菊；簡單重復(fù)序列（SSR）；多態(tài)性引物；遺傳多樣性

中圖分類號：S682.1⁺10.32""文獻(xiàn)標(biāo)志碼：A

文章編號：1002-1302（2024）22-0023-08

萬壽菊（Tagetes erecta）為菊科萬壽菊屬一年生草本植物，因其具有花色艷麗、花大、花期長、適應(yīng)性強(qiáng)等優(yōu)點(diǎn)而被廣泛應(yīng)用于園林景觀布置中。萬壽菊花朵富含葉黃素，是食用色素及保健品的天然植物來源^［1-3］。此外，萬壽菊的根部富含噻吩類化合物，可用于生產(chǎn)生物殺蟲劑、替代化學(xué)藥劑，有顯著的生態(tài)效益^［4］。

分子標(biāo)記是進(jìn)行植物遺傳多樣性研究及品種鑒定的有效工具之一。簡單重復(fù)序列（simple sequence repeat，SSR）為共顯性標(biāo)記，具有數(shù)量多、特異性好、多態(tài)性信息豐富、易于檢測、可重復(fù)性好等優(yōu)點(diǎn)，被廣泛應(yīng)用于種質(zhì)資源遺傳多樣性研究、遺傳關(guān)系與遺傳結(jié)構(gòu)分析、遺傳圖譜構(gòu)建、分子標(biāo)記輔助育種等方面^［5］。在萬壽菊屬植物分子標(biāo)記開發(fā)與利用方面的研究取各了很多成果。例如，何燕紅利用（簡單序列重復(fù)區(qū)間擴(kuò)增多態(tài)性（ISSR）、相關(guān)序列擴(kuò)增多態(tài)性（SRAP）、擴(kuò)增片段長度多態(tài)性（AFLP）標(biāo)記結(jié)合混合分組分析法（BSA）技術(shù)，篩選得到與萬壽菊雄性不育基因緊密連鎖的1個SRAP標(biāo)記和4個AFLP標(biāo)記^［6］。牟丹利用ISSR、隨機(jī)擴(kuò)增多態(tài)性（RAPD）、SRAP標(biāo)記篩選出與萬壽菊雄性不育相關(guān)的分子標(biāo)記，并轉(zhuǎn)化成穩(wěn)定的特定序列擴(kuò)增（SCAR）標(biāo)記用于萬壽菊雄性不育、雄性可育植株的苗期鑒定^［7］。侯志強(qiáng)等利用開發(fā)的萬壽菊SRAP標(biāo)記對70份萬壽菊雜交后代優(yōu)選株系進(jìn)行遺傳多樣性和家系分析，為育種親本選配提供信息^［8］。Panwar等利用ISSR、通用引物（URP）標(biāo)記分析了22個萬壽菊基因型的遺傳多樣性^［9］。張華麗等進(jìn)行了菊花（Chrysanthemum morifolium）表達(dá)序列標(biāo)簽SSR（EST-SSR）標(biāo)記的開發(fā)及其對萬壽菊的可轉(zhuǎn)移性分析，結(jié)果表明，14對有效引物中有8對可以在萬壽菊中擴(kuò)增出清晰穩(wěn)定的條帶，可轉(zhuǎn)移率為57.14%^［10］。Majumder等利用ISSR分子標(biāo)記對誘變獲得的萬壽菊突變體進(jìn)行變異分析^［11］。汪歡笑基于表型性狀、SSR標(biāo)記分析了萬壽菊種質(zhì)的遺傳多樣性，構(gòu)建了當(dāng)前的萬壽菊初級核心種質(zhì)^［12］。

筆者所在課題組前期利用8對SSR引物對9個萬壽菊雄性不育兩用系材料進(jìn)行了遺傳多樣性分析^［13］。楊潔利用4對SSR引物構(gòu)建了6個萬壽菊屬品種（系）的指紋圖譜，為萬壽菊屬品系的特異性、真實(shí)性鑒定奠定了理論基礎(chǔ)^［14］。但是，目前關(guān)于萬壽菊SSR標(biāo)記開發(fā)與利用的報(bào)道還相對較少。因此，本研究以筆者所在課題組前期工作中對萬壽菊雄性不育兩用系2-2可育、不育花蕾的轉(zhuǎn)錄組數(shù)據(jù)為基礎(chǔ)，開發(fā)SSR標(biāo)記并對26份萬壽菊種質(zhì)資源進(jìn)行遺傳多樣性分析，以期為萬壽菊核心種質(zhì)的選擇、新品種的選育提供理論依據(jù)。

1"材料與方法

1.1"試驗(yàn)材料

本研究于2023年3—12月在青海大學(xué)高原花卉研究中心分子實(shí)驗(yàn)室進(jìn)行，以該中心提供的26份萬壽菊品種（系）為材料，其中高稈品系10份，中稈品種2份，矮稈品（種）系14份（表1）。其中，雪域3號、雪域1號為品種，D-1×C-19、D-1×C-21、YS、G10為雜交后代，C-19、C-21、C-13、B-4為常規(guī)品系，其余16份材料為雄性不育兩用系。各品種（系）于盛花期隨機(jī)選取3個單株（雄性不育兩用系取可育單株）的幼嫩葉片，用液氮冷凍后置于-80 ℃保存?zhèn)溆谩?/p>

1.2"DNA提取及質(zhì)量檢測

采用改良十六烷基三甲基溴化銨（CTAB）法提取供試材料的植物基因組DNA，用1%瓊脂糖凝膠電泳、核酸蛋白儀檢測DNA的質(zhì)量與濃度，將樣品于-20 ℃ "冰箱中保存?zhèn)溆?sup>［15］。

1.3"SSR位點(diǎn)的鑒別

筆者所在課題組前期委托廣州基迪奧生物科技有限公司對萬壽菊雄性不育兩用系2-2的可育株、不育株花蕾進(jìn)行二代轉(zhuǎn)錄組測序及數(shù)據(jù)組裝^［16］。使用MISA程序查找SSR位點(diǎn)，查找標(biāo)準(zhǔn)如下：二核苷酸、三核苷酸、四核苷酸、五核苷酸、六核苷酸的重復(fù)次數(shù)分別≥7、≥5、4、4、4次。

1.4"引物設(shè)計(jì)與篩選

用Primer Premier 5.0設(shè)計(jì)引物，具體標(biāo)準(zhǔn)如下：引物長度為18～25 bp，GC含量為40%～60%，退火溫度（Tm）為51～59 ℃，預(yù)計(jì)產(chǎn)物長度為200～300 bp。隨機(jī)選取60對SSR引物，以萬壽菊雄性不育兩用系2-2基因組DNA為模板進(jìn)行PCR擴(kuò)增。擴(kuò)增體系及反應(yīng)程序參照楊潔等的方法^［14］，總量為25 μL的體系中包含12.5 μL 2×Taq PCR Master Mix、1 μL模板DNA、各1 μL上下游引物、9.5 μL滅菌雙蒸水。反應(yīng)程序如下：94 ℃ 2 min；94 ℃ 30 s，51～60 ℃ 30 s，72 ℃ 1 min，38個循環(huán)；72 ℃ 10 min，4 ℃ 保存^［14］。首先用1%瓊脂糖凝膠電泳檢測擴(kuò)增產(chǎn)物，初步篩選出條帶清晰、明亮且擴(kuò)增產(chǎn)物符合預(yù)期大小的引物用于多態(tài)性引物篩選。

1.5"多態(tài)性引物的篩選及數(shù)據(jù)分析

利用初篩獲得的SSR引物對26份萬壽菊供試材料進(jìn)行PCR擴(kuò)增，經(jīng)瓊脂糖凝膠電泳檢測后，用6%變性聚丙烯酰胺凝膠電泳進(jìn)行引物的多態(tài)性分析，在220 V電壓條件下電泳90 min，用硝酸銀染色法顯色后觀察結(jié)果并拍照記錄。

計(jì)算SSR發(fā)生頻率（包含SSR位點(diǎn)的序列數(shù)/序列總數(shù)）、SSR出現(xiàn)頻率（識別的SSR位點(diǎn)數(shù)/序列總數(shù)）和平均分布距離（總堿基數(shù)/識別的SSR位點(diǎn)數(shù)）。用軟件Popgene 1.32計(jì)算有效等位基因數(shù)（Na）、香農(nóng)信息指數(shù)（I）、觀察雜合度（Ho）、期望雜合度（He）、多態(tài)信息含量（PIC）、遺傳距離。用軟件NTSYS-pc 2.2進(jìn)行UPGMA聚類分析。

2"結(jié)果與分析

2.1"萬壽菊轉(zhuǎn)錄組中SSR位點(diǎn)的分布與特點(diǎn)

萬壽菊雄性不育兩用系轉(zhuǎn)錄組的總長度為

在10 182個SSR位點(diǎn)中，三核苷酸為主要重復(fù)類型，共有4 939個，占總SSR位點(diǎn)的48.51%，其次分別是二核苷酸、四核苷酸、五核苷酸、六核苷酸，分別占總SSR位點(diǎn)的27.05%、11.57%、6.71%、6.17%。從重復(fù)次數(shù)看，以5次重復(fù)的SSR位點(diǎn)最多，共3 213個，占SSR位點(diǎn)總數(shù)的31.56%，其次分別是6、4、7次重復(fù)，分別占25.50%、17.87%、12.60%，20次以上重復(fù)的SSR位點(diǎn)只有52個，占0.51%（表2）。

2.2"萬壽菊SSR基序重復(fù)類型及頻率特征

10 182個SSR位點(diǎn)中包含259種基序重復(fù)類型，其中二核苷酸、三核苷酸、四核苷酸、五核苷酸、六核苷酸的重復(fù)類型分別有3、10、30、66、150 種。二核苷酸重復(fù)基序主要有AC/GT、AG/CT、AT/AT等3種（圖1），其中AG/CT重復(fù)類型最多，有1 025個，占二核苷酸重復(fù)總數(shù)的37.22%，占SSR位點(diǎn)總數(shù)的10.06%。三核苷酸中出現(xiàn)最多的重復(fù)類型是ATC/ATG（1 414個），占三核苷酸總數(shù)的28.63%，占SSR位點(diǎn)總數(shù)的13.89%。四核苷酸中出現(xiàn)最多的重復(fù)類型為ACAT/ATGT（205個），占四核苷酸重復(fù)總數(shù)的17.40%，占SSR位點(diǎn)總數(shù)的2.01%。五核苷酸重復(fù)類型中，每種重復(fù)基序類型的數(shù)量均低于100個，其中出現(xiàn)最多的重復(fù)類型是AATTC/AATTG（75個）， 占五核苷酸重復(fù)總數(shù)10.98%， 占SSR位點(diǎn)總數(shù)的0.74%。六核苷酸重復(fù)基元中出現(xiàn)最多的是AAGATG/ATCTTC（39個），占六核苷酸重復(fù)總數(shù)的6.21%，占SSR位點(diǎn)總數(shù)的0.38%。

2.3"萬壽菊SSR標(biāo)記篩選及多態(tài)性分析

引物初篩結(jié)果表明，60對引物中有41對可以擴(kuò)增出清晰明亮的條帶。用這41對引物對26份萬壽菊屬品種（系）基因組DNA進(jìn)行PCR擴(kuò)增，結(jié)果表明，編號為U104、U113、U203、U335、U461、U1205、U1203、U1757、U2041、U2261的10個SSR標(biāo)記能穩(wěn)定擴(kuò)增出目標(biāo)條帶且具有多態(tài)性（圖2），多態(tài)性位點(diǎn)比例為24.39%。

10個SSR多態(tài)性標(biāo)記在26份萬壽菊材料中的有效基因數(shù)平均1.614個，有效等位基因數(shù)最少的為1.122個（U1205），最多的為1.997個（U335）。Shannon’s多樣性指數(shù)（I）的分布范圍為0.221～0.692，平均為0.539；觀測雜合度（Ho）、期望雜合度（He）的變化范圍分別為0.039～0.654、0.109～0.499，平均分別為0.273、0.361；多態(tài)性信息含量（PIC）的變化范圍為0.114～0.503，平均為0.364，其中U335的多態(tài)性信息含量大于0.5（表3）。

2.4"基于SSR標(biāo)記的萬壽菊遺傳多樣性分析

26份萬壽菊屬品種（系）的有效等位基因數(shù)（Ne）范圍為1.600～2.000個，平均值為1.841個；Shannon^，s多樣性指數(shù)（I）的范圍為0.562～0.693，平均值為0.642；期望雜合度（He）在0.395～0.526之間，平均值為0.475（表4）。上述結(jié)果說明，26份萬壽菊供試材料的遺傳多樣性豐富。遺傳距離介于0.027～0.913之間（表5），其中D-1和5-2間的遺傳距離最大（0.913），C-16與B-4、C-13間的遺傳距離最?。?.027），說明D-1和5-2間的親緣關(guān)系最遠(yuǎn)，C-16與B-4、C-13這2個品系間的親緣關(guān)系最近。

UPGMA聚類分析結(jié)果表明，當(dāng)相似系數(shù)為0.69時，26份萬壽菊材料可分為兩大類（圖3）。結(jié)合表型性狀可知，這兩大類的劃分與花色具有相關(guān)性，其中第一大類（A）由15個品種（系）組成，花色以淺黃色為主，第二大類（B）由11個品種（系）組成，花色以橙黃色為主。當(dāng)相似系數(shù)為0.77時，第一大類可進(jìn)一步分為2個亞類，其中第1亞類（A1）由12個品系組成，大部分為矮稈，第二亞類（A2）包括3個品系，均為高稈。第二大類（B）進(jìn)一步也可分為2個亞類，第1亞類（B1）包括9個品種（系），第2亞類（B2）僅包含C-20、2-2等2個品系。

3"結(jié)論與討論

萬壽菊有重要的觀賞、藥用和經(jīng)濟(jì)價值，是國內(nèi)外廣泛栽植的一年生草本花卉。本研究從萬壽菊轉(zhuǎn)錄組中挖掘到10 182個SSR位點(diǎn)，其出現(xiàn)頻率為16.07%。在已有報(bào)道中，不同植物材料中SSR的出現(xiàn)頻率相差較大，張華麗等從萬壽菊雄性不育系S-121中挖掘出13 849個SSR位點(diǎn)，出現(xiàn)頻率為28.29%；盛文濤等在甜葉菊（Stevia rebaudiana）的5 646個unigenes中挖掘出88個SSR位點(diǎn)，出現(xiàn)頻率僅為1.60%。上述結(jié)果說明，SSR的出現(xiàn)頻率可能與物種本身基因組有關(guān)，也可能受到挖掘SSR位點(diǎn)的方式及篩選標(biāo)準(zhǔn)的影響^［17-20］。

本研究挖掘的SSR位點(diǎn)中包含259種基序重復(fù)類型，其中主導(dǎo)重復(fù)類型為二核苷酸、三核苷酸，分別占SSR位點(diǎn)總數(shù)的27.05%、48.51%。AG/CT、ATC/ATG在二核苷酸、三核苷酸重復(fù)類型總數(shù)中的占比最大，分別為37.22%、28.63%，這與張華麗等的研究結(jié)果^［17］一致。此外，金雞菊（Coreopsis lanceolata）、向日葵（Helianthus annuus）、黃纓菊（Xanthopappus subacaulis）、黃花蒿（Artemisia annua）中開發(fā)的SSR位點(diǎn)也表現(xiàn)出相同的主導(dǎo)重復(fù)類型^［21-24］。在乒乓菊（Dendranthema morifolium cv. pompon）、非洲菊（Gerbera jamesonii）的SSR位點(diǎn)中，單核苷酸為主要重復(fù)類型，并且A/T占比最大，與萬壽菊基序重復(fù)類型存在較大差異^［25-26］。由此可以看出，同科植物的SSR重復(fù)類型沒有表現(xiàn)出明顯相同的規(guī)律，但二核苷酸、三核苷酸是最常見的SSR重復(fù)類型^［27］。

本研究從隨機(jī)選取的60個SSR標(biāo)記中篩選出10個（24.39%）多態(tài)性標(biāo)記，PIC平均值為0.361。當(dāng)0.25lt;PIClt;0.5時，表示標(biāo)記具有中度多態(tài)性，本研究選擇的5個標(biāo)記（U203、U335、U1757、U2041、U2261）的PIC值高于平均值，表明其變異程度高、多態(tài)性好，可以在萬壽菊種質(zhì)資源評價和品種鑒定等研究中作為重點(diǎn)標(biāo)記加以應(yīng)用^［28］。利用10個多態(tài)性標(biāo)記對26份萬壽菊屬品種（系）進(jìn)行遺傳多樣性分析，發(fā)現(xiàn)供試材料具有豐富的遺傳多樣性， 其Shannon^，s信息指數(shù)（I）、 期望雜合度（He）的平均值分別為0.642、0.475，高于張華麗等報(bào)道的20份萬壽菊自交系的Shannon^，s信息指數(shù)（0.408）、期望雜合度（0.275）^［17］。UPGMA聚類結(jié)果顯示，當(dāng)相似系數(shù)為0.69時，26份供試材料可分為花色不同的兩大類，這與風(fēng)信子（Hyacinthus orientalis）的研究結(jié)果^［29］一致?？傮w來看，萬壽菊轉(zhuǎn)錄組中SSR出現(xiàn)頻率較高、類型豐富且多態(tài)性高，為萬壽菊種質(zhì)資源鑒別與利用、遺傳圖譜構(gòu)建和分子標(biāo)記輔助選擇等提供了候選分子標(biāo)記，遺傳多樣性分析為今后萬壽菊的新品種選育奠定了重要的理論基礎(chǔ)。

參考文獻(xiàn)：

［1］史金寶. 萬壽菊不同品種生長發(fā)育動態(tài)及農(nóng)藝性狀比較［D］. 烏魯木齊：新疆農(nóng)業(yè)大學(xué)，2022：3-9.

［2］曾晨峰，董"朕，劉"華，等. 萬壽菊葉黃素對蛋雞蛋品質(zhì)、血清生化指標(biāo)、 免疫功能和抗氧化能力的影響［J］. 飼料工業(yè)，2024，45（2）：49-55.

［3］李"峰，蔣高華，芮曉慶. 萬壽菊中葉黃素的提取工藝研究進(jìn)展［J］. 山東化工，2023，52（3）：98-99，103.

［4］王新國，徐漢虹，趙善歡. 殺蟲植物萬壽菊的研究進(jìn)展［J］. 西安聯(lián)合大學(xué)學(xué)報(bào)，2002，5（2）：5-10.

［5］李福榮. 萬壽菊雄性不育的遺傳與應(yīng)用研究［D］. 北京：北京林業(yè)大學(xué)，2005：6-10.

［6］何燕紅. 萬壽菊雄性不育性狀的遺傳分析及其育種應(yīng)用［D］. 武漢：華中農(nóng)業(yè)大學(xué)，2010：3-7.

［7］牟"丹. 萬壽菊雄性不育兩用系苗期分子標(biāo)記及應(yīng)用［D］. 西寧：青海大學(xué)，2016：9-12.

［8］侯志強(qiáng)，唐"楠，冶有林，等. 萬壽菊雜交后代優(yōu)選株系的SRAP-PCR遺傳多樣性分析［J］. 基因組學(xué)與應(yīng)用生物學(xué)，2017，36（10）：4295-4306.

［9］Panwar S，Singh K P，Banyal N，et al. DNA fingerprinting in African marigold （Tagetes erecta L.） genotypes using ISSR and URP markers［J］. Indian Journal of Horticulture，2018，75（1）：105.

［10］張華麗，王"濤，宋利娜，等. 菊花EST-SSR標(biāo)記開發(fā)及其對萬壽菊的可轉(zhuǎn)移性研究［J］. 核農(nóng)學(xué)報(bào)，2016，30（1）：42-49.

［11］Majumder J，Singh S K，Kumari M，et al. Variability and correlation studies on induced mutants of marigold （Tagetes erecta L.） for different traits and assessing them using molecular markers［J］. Plant Tissue Culture and Biotechnology，2018，28（2）：223-236.

［12］汪歡笑. 基于表型性狀和SSR標(biāo)記多樣評價構(gòu)建萬壽菊初級核心種質(zhì)［D］. 秦皇島：河北科技師范學(xué)院，2023：10-13.

［13］宋江琴，唐"楠，唐道城，等. 基于表型性狀和SSR標(biāo)記的9份萬壽菊種質(zhì)遺傳多樣性分析［J］. 種子，2021，40（10）：6-11，19.

［14］楊"潔. 萬壽菊‘雪域3號’品種特征及指紋圖譜構(gòu)建［D］. 西寧：青海大學(xué)，2021：10-15.

［15］李"巖，于海琳，馬曉彤，等. 對CTAB法提取菠菜基因組DNA實(shí)驗(yàn)的優(yōu)化設(shè)計(jì)［J］. 化工管理，2023（14）：15-17.

［16］Tang N，Liu W，Zhang W H，et al. Integrative analysis of transcriptomic and proteomic changes related to male sterility in Tagetes erecta［J］. Physiology and Molecular Biology of Plants，2020，26（10）：2061-2074.

［17］張華麗，叢日晨，王茂良，等. 基于萬壽菊轉(zhuǎn)錄組測序的SSR標(biāo)記開發(fā)［J］. 園藝學(xué)報(bào)，2018，45（1）：159-167.

［18］盛文濤，鄧建蘭，饒友生，等. 甜葉菊EST序列特征分析與分子標(biāo)記開發(fā)［J］. 分子植物育種，2019，17（15）：5045-5050.

［19］Biswas M K，Chai L J，Mayer C，et al. Exploiting BAC-end sequences for the mining，characterization and utility of new short sequences repeat （SSR） markers in citrus［J］. Molecular Biology Reports，2012，39（5）：5373-5386.

［20］Jarne P，Lagoda P J. Microsatellites，from molecules to populations and back［J］. Trends in Ecology amp; Evolution，1996，11（10）：424-429.

［21］孫浩男，李"蓉，李明陽，等. 栽培金雞菊EST-SSR分子標(biāo)記開發(fā)與親緣關(guān)系分析［J］. 草地學(xué)報(bào)，2022，30（6）：1430-1440.

［22］石"煜. 向日葵EST-SSR標(biāo)記的開發(fā)及應(yīng)用［D］. 呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2016：31-34.

［23］鄭長遠(yuǎn). 黃纓菊SSR分子標(biāo)記開發(fā)及遺傳多樣性研究［D］. 西寧：青海師范大學(xué)，2023：45-51.

［24］王"英，蔡漢強(qiáng)，賈偉章. 黃花蒿EST資源的SSR信息分析及標(biāo)記開發(fā)［J］. 中國中藥雜志，2012，37（5）：570-574.

［25］周"琳，蔡友銘，張永春，等. 基于RNA-Seq技術(shù)的乒乓菊轉(zhuǎn)錄組分析［J］. 分子植物育種，2020，18（18）：5917-5924.

［26］林發(fā)壯，李錦燁，安慧珍，等. 基于非洲菊轉(zhuǎn)錄組測序的SSR位點(diǎn)信息分析［J］. 福建農(nóng)業(yè)科技，2020（11）：1-6.

［27］余興華，楊艷梅，李玉祥，等. 基于八角轉(zhuǎn)錄組序列的SSR分子標(biāo)記引物開發(fā)［J］. 西南農(nóng)業(yè)學(xué)報(bào)，2022，35（6）：1261-1268.

［28］宋江琴. 萬壽菊矮稈兩用系及雜交后代遺傳多樣性分析［D］. 西寧：青海大學(xué)，2021：32-35.

［29］何元浩，胡鳳榮. 基于AFLP分子標(biāo)記的30個風(fēng)信子品種遺傳多樣性分析［J］. 分子植物育種，2023，21（9）：2955-2963.