粗葉榕轉(zhuǎn)錄組SSR位點分析及分子標記開發(fā)

摘""要：五指毛桃（Fici"Hirtae"Radix）是一種中藥材，其來源于粗葉榕（Ficus"hirta"Vahl.）的干燥根部。開發(fā)利用藥用植物粗葉榕SSR分子標記，為其種質(zhì)遺傳多樣性分析、遺傳圖譜構(gòu)建、輔助育種等提供理論依據(jù)及數(shù)據(jù)支撐。以粗葉榕全緣、三裂、五裂、七裂4種典型葉片為試材，經(jīng)轉(zhuǎn)錄組測序后進行位點特征分析及分子標記開發(fā)。結(jié)果表明：基于34"375條轉(zhuǎn)錄組unigene序列，篩選到16"851個SSR位點，分布于12"094條unigene，分布頻率為35.18%，平均每2638"bp有1個SSR位點，分布密度為3.68"kb。共觀察到5種不同核苷酸基序類型，其中最常見的是重復6次的基序，占SSR總數(shù)的20.92%，且長度基本小于29"bp。雙堿基重復序列占主導地位，其次是三堿基和四堿基重復類型。其中，（AG/CT）n是最常見的基元，占所有SSR的26.05%；其次是（AT/AT）n基元，占18.69%。隨機選取100條EST-SSR引物，以粗葉榕4種不同葉型的基因組DNA擴增目標產(chǎn)物，其中89條引物在粗葉榕4個種質(zhì)間表現(xiàn)出多態(tài)性。選取9對條帶清晰的引物，對不同地區(qū)的31個粗葉榕種質(zhì)進行DNA擴增和指紋圖譜構(gòu)建，這些引物組合多態(tài)性條帶比例均值達72.04%，并對31個種質(zhì)進行了有效鑒別，成功建立了DNA指紋圖譜。研究結(jié)果證明粗葉榕轉(zhuǎn)錄組SSR位點信息豐富，利用性高，可為后續(xù)研究提供數(shù)據(jù)支撐。

關(guān)鍵詞：粗葉榕；五指毛桃；轉(zhuǎn)錄組；EST-SSR；DNA指紋圖譜中圖分類號：S567.19""""""文獻標志碼：A

Comprehensive"Analysis"of"SSR"Loci"Within"the"Ficus"hirta"Vahl."Transcriptome"and"the"Developmental"Application"of"Molecular"Markers

YAN"Qing1，"QIN"Yanan1，3，"ZHANG"Hualin1，"LI"Wenxiu1，"YE"Yana1，"HUANG"Anying2，"LUO"Ping1*

1."Zhanjiang"Experimental"Station，"Chinese"Academy"of"Tropical"Agricultural"Sciences"/"Zhanjiang"Rubber"Understory"Economic"Engineering"Technology"Research"Center，"Zhanjiang，"Guangdong"524013，"China;"2."Research"Institute"of"Fast-growing"Trees，"Chinese"Academy"of"Forestry，"Zhanjiang，"Guangdong"524300，"China;"3."College"of"Tropical"Crops，"Yunnan"Agricultural"University，"Pu’er，"Yunnan"650201，"China

Abstract："Fici"Hirtae"Radix，"a"Chinese"herbal"medicine，"is"prepared"from"the"dry"roots"of"Ficus"hirta"Vahl."This"study"aimed"to"provide"theoretical"and"empirical"support"for"the"genetic"diversity"analysis，"genetic"map"construction，"and"assisted"breeding"of"F."hirta"Vahl."by"developing"SSR"molecular"markers."Four"typical"leaves"of"F."hirta"Vahl."were"selected"as"the"samples，"and"site"feature"analysis"and"molecular"marker"development"were"performed"after"transcriptome"sequencing."A"total"of"16"851"SSR"loci"were"identified"from"34"375"transcriptome"unigene"sequences，"distributed"across"12"094"unigene"sequences，"resulting"in"a"distribution"frequency"of"35.18%."The"average"frequency"was"one"SSR"locus"every"2638"bp，"with"a"density"of"3.68"kb."Five"distinct"nucleotide"motifs"were"observed，nbsp;with"the"motif"repeating"six"times"being"the"most"prevalent，"constituting"20.92%"of"all"SSRs，"and"typically"less"than"29"bp"in"length."Di-nucleotide"repeats"were"predominant，"followed"by"tri-and"tetra-nucleotide"repeats."（AG/CT）n"was"the"most"abundant"repeat"motif，"representing"26.05%"of"all"SSRs，"while"（AT/AT）n"was"the"second"most"common，"accounting"for"18.69%."100"EST-SSR"primers"were"randomly"selected"to"amplify"target"product"using"genomic"DNA"from"four"different"leaf"types"of"F."hirta"Vahl."89"primers"exhibited"polymorphism"across"the"germplasms."Nine"primer"pairs，"distinguished"by"the"high"definition"and"clarity，"were"selected"to"perform"DNA"amplification"and"construct"DNA"fingerprints"for"31"specimens"sourced"from"different"regions."The"average"percentage"of"polymorphic"bands"observed"across"these"primers"was"72.04%."The"substantial"degree"of"polymorphism"facilitated"the"unambiguous"identification"of"each"of"the"31"germplasm"samples，"and"robust"DNA"fingerprint"profiles"were"successfully"established"for"each"specimen."This"study"demonstrated"that"SSR"loci"within"the"transcriptome"of"F."hirta"Vahl."were"rich"in"information"and"highly"utilizable，"offering"valuable"data"support"for"future"research"endeavors.

Keywords："Ficus"hirta"Vahl.;"Fici"Hirtae"Radix;"transcriptome;"EST-SSR;"DNA"fingerprinting

DOI："10.3969/j.issn.1000-2561.2025.05.007

微衛(wèi)星標記，又稱簡單重復序列（simple"sequence"repeat，"SSR），是一種在植物遺傳多樣性研究、基因定位、目標性狀定位以及標記輔助選擇育種中被廣泛使用的標記系統(tǒng)。它具有突變率高、分布廣、密度大的優(yōu)點，已經(jīng)廣泛應用于種質(zhì)資源鑒別、遺傳多樣性分析以及指紋圖譜的構(gòu)建等領域[1]。SSR主要是由串聯(lián)重復單元構(gòu)成，其核心單位含1～6個核苷酸，因重復次數(shù)和程度的差異，呈現(xiàn)出高度的多態(tài)性。SSR標記具有共顯性特征，穩(wěn)定性強且等位基因多樣性豐富，根據(jù)序列獲得方式不同，可分為基因組SSR（genomic"SSR，"g-SSR）和轉(zhuǎn)錄組SSR（expressed"sequence"tags"SSR，"EST-SSR）[2]。與g-SSR相比，EST-SSR并不需要對整個基因組進行測序，而是位于基因編碼區(qū)內(nèi)，具有方便快捷、效率高、成本低廉及與表型性狀相關(guān)性好等優(yōu)點。傳統(tǒng)的SSR開發(fā)方法難度大、成本高、勞動密集，高通量測序技術(shù)的進步和測序成本的降低減少了轉(zhuǎn)錄組測序的成本，推動了通過轉(zhuǎn)錄組測序開發(fā)SSR標記的快速發(fā)展與廣泛應用。轉(zhuǎn)錄組測序后生成的序列數(shù)據(jù)中包含大量的有效基因，為評估物種多樣性提供了新方法。目前，通過轉(zhuǎn)錄組數(shù)據(jù)開發(fā)EST-SSR標記已在黃連、冬蟲夏草、白木香、川貝母、穿心蓮等諸多植物中得到廣泛應用[3]，極大促進了植物的育種進程。

粗葉榕（Ficus"hirta"Vahl.），屬于?？崎艑俚臒o花果亞屬，是一種木本植物，主要分布于亞洲的亞熱帶地區(qū)，中國的熱帶及亞熱帶地區(qū)均有分布[4]。中藥材五指毛桃（Fici"Hirtae"Radix）藥食同源，來源于粗葉榕的干燥根，又稱五指牛奶、南芪、土北芪、五爪龍等，具有益氣、健脾、潤肺、祛風濕等功效，在海南、廣東等地被廣泛應用，藥用和食用價值較高[5]。粗葉榕的高遺傳多樣性與低群體分化決定了其更高的適應性與進化潛力[6]，在長期適應不同生境及遺傳漂移的過程中形成了復雜的變型。從已收集的種質(zhì)資源來看，粗葉榕不同種質(zhì)間表型差異大，僅通過表型鑒定無法準確及時對其進行區(qū)分，同時不同種質(zhì)間遺傳背景不清晰，遺傳距離不明了，極大限制了植物新品種選育，影響其育種進程。因此，表型評價結(jié)合分子標記技術(shù)無疑是一個極佳的解決手段，通過轉(zhuǎn)錄組測序開發(fā)粗葉榕SSR標記，對不同種質(zhì)資源進行精準鑒定，對后續(xù)粗葉榕資源評價及基因組選擇育種具有重要意義。

當前對粗葉榕的報道多集中于傳統(tǒng)分類學[7]、化學成分[8-9]、繁育[10]與藥理作用[11]等方面，少有品種選育與精準鑒別方面的研究，缺少豐富的分子標記資源和準確的DNA指紋圖譜信息。國內(nèi)外對大豆[12]、玉米[13]、水稻[14]等作物的研究結(jié)果表明，SSR標記可有效檢測位點多態(tài)性，對育種工作具有較好的推動作用。而關(guān)于粗葉榕轉(zhuǎn)錄組SSR分子標記開發(fā)利用方面的研究報道較少，僅有LU等[6]基于同屬植物無花果（Ficus"carica"Linn.）的9對SSR引物對18個粗葉榕居群開展的遺傳多樣性分析。截至目前，已有學者開展了全緣、三裂、五裂、七裂4種葉型五指毛桃果實的氨基酸組成比較及營養(yǎng)價值評價[15]，尚無通過不同葉裂類型粗葉榕轉(zhuǎn)錄組測序開發(fā)SSR標記的研究。本研究以不同葉裂粗葉榕葉片轉(zhuǎn)錄組測序為基礎，利用生物信息學方法分析其SSR位點序列分布特征、堿基重復類型等，開發(fā)了4種典型葉型粗葉榕中均具有的多態(tài)性的葉片SSR標記，并在31種不同來源種質(zhì)中進行驗證。為粗葉榕種質(zhì)遺傳多樣性分析、遺傳圖譜構(gòu)建、輔助育種等提供理論及數(shù)據(jù)支撐。

1""材料與方法

1.1""材料

供試材料為河源全緣、河源三裂、河源五裂、河源七裂4個粗葉榕典型種質(zhì)（圖1）及31份不同來源粗葉榕種質(zhì)（表1），經(jīng)中國熱帶農(nóng)業(yè)科學院湛江實驗站羅萍研究員鑒定為粗葉榕。樣品采自中國熱帶農(nóng)業(yè)科學院湛江實驗站粗葉榕種質(zhì)圃，該種質(zhì)圃位于廣東省東南部，年均氣溫為23.8"℃，年均降水量為1396.3～1759.4"mm。選取生長良好，已分化且分化明顯的不同種質(zhì)新鮮葉片各3片，液氮快速冷凍后置于?80"℃冰箱保存，備用。

1.2""方法

1.2.1""轉(zhuǎn)錄組測序及組裝""采集典型種質(zhì)河源全緣、河源三裂、河源五裂、河源七裂4個種質(zhì)各3份葉片，提取RNA后分別反轉(zhuǎn)錄構(gòu)建cDNA文庫，然后利用Illumina"HiSeq"TM"2000平臺測序。使用Trinity軟件對原始數(shù)據(jù)進行過濾、重組、去除冗余和同源轉(zhuǎn)錄本聚類，從而得到unigene序列。

1.2.2""DNA提取""取不同葉型及不同來源的12個粗葉榕種質(zhì)葉片分別研磨，采用改進的CTAB方法提取DNA，并通過1%瓊脂糖凝膠電泳檢查其完整性。使用時檢測濃度與純度，最后將其稀釋到一致的濃度后再使用。

1.2.3""轉(zhuǎn)錄組SSR位點鑒別及EST-SSR引物設計""利用MISA"1.0軟件對粗葉榕的轉(zhuǎn)錄組數(shù)據(jù)查找SSR位點時，重復序列的要求為：二核苷酸≥6；三核苷酸≥5；四核苷酸、五核苷酸、六核苷酸均≥4。此外，如果2個SSR之間的距離小于100"bp，則認為是1個SSR。使用Primer"3.0軟件設計引物，設計預測產(chǎn)物長度為100～280"bp，序列長度為8～25"bp，GC含量為40%～60%，退火溫度為（60±5）℃，上下游引物的退火溫度之間的差異不超過5"℃。

1.2.4""EST-SSR引物篩選及遺傳多樣性分析""（1）EST-SSR引物篩選。從引物組中隨機選取100對長度超過20"bp的SSR序列，以全緣、三裂、五裂、七裂4種不同種質(zhì)粗葉榕葉片DNA為模板進行擴增。PCR擴增體系為10"μL：上下游引物各0.4"μL，1"μL"DNA，5"μL"Mix，3.2"μL"ddH2O。PCR擴增程序為：98"℃預變形5"min，98"℃變性10"s，60"℃退火15"s，72"℃延伸45"s，35個循環(huán)，72"℃延伸2"min。使用12%聚丙烯酰胺凝膠電泳對PCR擴增產(chǎn)物進行檢測，設定電壓為200"V，電流為120"A，電泳2"h。電泳結(jié)束后，進行銀染洗脫，記錄結(jié)果。

（2）EST-SSR引物多態(tài)性分析。選取來自廣東、廣西、福建和云南4個?。▍^(qū)）的31份粗葉榕種質(zhì)資源（表1）進行EST-SSR引物驗證。以上述種質(zhì)DNA為模板，對1.2.4-（1）篩選得到的9對引物進行PCR擴增，電泳后讀帶。將讀帶結(jié)果進行統(tǒng)計，采用Popgen"32軟件分析不同引物條帶數(shù)量、觀測等位基因數(shù)（observed"number"of"alleles，"Na）、有效等位基因數(shù)（effective"number"of"alleles，"Ne）、香農(nóng)指數(shù)（Shannon’s"information"index，"I）、Nei’s基因多樣性指數(shù)（Nei’s"gene"diversity"index，"H）及不同種質(zhì)間的遺傳距離；采用Genalex軟件分析期望雜合度（expected"heterozygosity，"He）與無偏期望雜合度（unbiased"expected"heterozygosity，"uHe）。將數(shù)據(jù)輸入到NTSYS軟件中，利用非加權(quán)組平均法在NTSYS軟件的相似性模塊中建立聚類樹，具體方法：從NTSYS軟件的相似性模塊中獲取相似性矩陣或距離矩陣，在Clustering模塊中得到聚類樹，然后將種質(zhì)相關(guān)資料輸入二維碼制作工具，從而為各粗葉榕資源創(chuàng)建二維碼形式的DNA分子身份標識[16-17]。

2""結(jié)果與分析

2.1""轉(zhuǎn)錄組的組裝

通過對河源全緣、河源三裂、河源五裂、河源七裂4種典型種質(zhì)進行轉(zhuǎn)錄組測序，共得到12組序列。對轉(zhuǎn)錄組組裝質(zhì)量統(tǒng)計分析，共得到34"375條unigene，其總長度為44"450"502"bp，平均長度為1293"bp；最長的unigene長度為

14"243"bp，最短的unigene長度為201"bp；GC含量為41.50%；N50長度為2186"bp，是平均長度的1.69倍；N50數(shù)量為6808條，是總組裝數(shù)量的19.81%，說明粗葉榕轉(zhuǎn)錄組的組裝質(zhì)量好。

2.2""粗葉榕葉片轉(zhuǎn)錄組中SSR位點分析

利用MISA軟件進行分析，在粗葉榕的34"375條轉(zhuǎn)錄組unigene序列數(shù)據(jù)中，12"094條unigene上含有16"851個SSR位點，含SSR位點的unigene數(shù)占unigene總數(shù)的35.18%，而SSR位點的數(shù)量則占unigene總數(shù)的49.02%，平均每2638"bp有1個SSR位點。此外，還發(fā)現(xiàn)有4756個unigene包含至少1個以上的SSR位點。

2.3""粗葉榕葉片轉(zhuǎn)錄組SSR重復類型分析

從表2可以看出，在粗葉榕的轉(zhuǎn)錄組SSR位點中，雙堿基和三堿基重復類型的頻率較高，分別占24.83%和10.71%。四堿基、五堿基和六堿基重復類型的頻率相對較低，分別僅占3.14%、1.48%和1.63%。平均每19.88"bp堿基長度出現(xiàn)1個雙堿基重復序列，平均每19.56"bp堿基長度出現(xiàn)1個雙堿基至四堿基，而每24.16"bp堿基長度時才可能會出現(xiàn)五、六堿基，可見雙堿基、六堿基與其他類型堿基間的平均分布長度差異較大，特別是雙堿基在長度分布上占有很高的優(yōu)勢地位。雙堿基和三堿基的總長度達到了20多萬的數(shù)量級，剩余類型堿基長度總和也有236.48"kb，但四堿基、五堿基、六堿基類型總長度只有45.38"kb左右。在總的unigene序列上平均每2.64"kb出現(xiàn)1個重復基序，各種類型的基序出現(xiàn)距離由雙堿基、三堿基、四堿基、五堿基到六堿基依次增加。

2.4""粗葉榕葉片轉(zhuǎn)錄組中SSR位點類型

在粗葉榕轉(zhuǎn)錄組中，單一出現(xiàn)頻率高于2%的主要有7種重復單元，其中包含雙堿基重復3種，即AC/GT重復、AG/CT重復、AT/AT重復，尤以AG重復數(shù)量最多，為4389個，比例達26.05%。三堿基重復基元重復次數(shù)較多的共8種，AAG/CTT重復最多，出現(xiàn)頻率為7.62%，ATC/ATG重復較少，頻率為1.41%；四堿基重復中，數(shù)量最多的是AAAT/ATTT，出現(xiàn)頻率為3.41%；五堿基重復數(shù)量最多的為AAAAT/ATTTT單元。其他類型出現(xiàn)頻率低于0.90%的重復單元共占總SSR的21.97%（圖2）。

2.5""粗葉榕葉片轉(zhuǎn)錄組SSR重復次數(shù)分析

利用MISA軟件對轉(zhuǎn)錄組中的所有unigene進行搜索，并分析SSR基元的重復次數(shù)。SSR重復單元的重復次數(shù)分布在4～13次之間，其中4～9次重復的有14"235個，占78.46%；10～13次重復的SSR位點有3908個，占21.54%。6次重復基元出現(xiàn)的頻率最高，占總SSR位點數(shù)的20.92%；其次為5、7次重復基元出現(xiàn)的頻率，分別占總SSR位點數(shù)的14.98%、13.49%。從SSR的重復次數(shù)來看，雙堿基主要集中在6～11次，占雙堿基分布總數(shù)的86.01%；三堿基主要集中在5～7次，占三堿基分布總數(shù)的81.30%；四堿基主要集中在4～6次，占四堿基分布總數(shù)的96.21%；五堿基主要集中在4～5次，占五堿基分布總數(shù)的95.51%；六堿基主要集中在4～5次，占六堿基總數(shù)的94.84%（圖3）。

2.6""粗葉榕葉片轉(zhuǎn)錄組SSR可用性評價

篩選的SSR基序長度分布在12～407"bp范圍內(nèi)，在該范圍內(nèi)，基序長度為12～19"bp的中度多態(tài)性位點共有7631個，占比為45.29%；基序長度≥20"bp的高度多態(tài)性位點有9220個，占54.71%（圖4）。說明可以從該葉片轉(zhuǎn)錄組中開發(fā)出具有高度多態(tài)性的SSR引物。

2.7""EST-SSR引物篩選與遺傳多樣性分析

2.7.1""EST-SSR引物篩選""利用16"851個SSR位點，最終設計得到12"336對EST-SSR引物。以全緣、三裂、五裂、七裂4種不同葉型粗葉榕DNA為模板進行驗證，隨機選擇100對EST-SSR引物進行擴增。結(jié)果顯示，有89對引物可擴增出目標產(chǎn)物，其中55對引物在4種不同葉型粗葉榕中表現(xiàn)出多態(tài)性，篩選出9對條帶清晰的引物（表3）。

2.7.2""EST-SSR引物遺傳多樣性分析""以31份不同粗葉榕種質(zhì)DNA為模板，對篩選出的9對EST-SSR引物進行擴增及遺傳多樣性分析，并構(gòu)建系統(tǒng)進化樹及DNA分子身份證。利用9對EST-SSR引物對31份不同粗葉榕種質(zhì)進行擴增，計算獲得遺傳多樣性指標和多態(tài)性信息含量（表4）。9對引物的觀測多態(tài)性條帶數(shù)在3～5條之間，平均多態(tài)性條帶數(shù)量為4.44條；觀測等位基因數(shù)（Na）在0.903～1.839之間，平均值為1.577；有效等位基因數(shù)（Ne）為1.273～1.656，平均值為1.476；香農(nóng)指數(shù)（I）介于0.235～0.543之間，平均值為0.405；期望雜合度（He）介于0.158～0.372之間，平均值為0.274；無偏期望雜合度（uHe）介于0.181～0.414之間，平均值為0.309；Nei’s基因多樣性指數(shù)（H）為0.158～0.372，平均值為0.274。多態(tài)性條帶比例分布范圍在41.94%～"90.32%之間，平均值為72.04%。結(jié)果表明，篩選的9對引物在31份種質(zhì)中的多態(tài)性較高。

31份粗葉榕種質(zhì)的遺傳相似系數(shù)在0.35～"0.90之間，種質(zhì)間均存在遺傳差異。CY1與FLT7的遺傳相似系數(shù)最?。?.35），遺傳差異最大；FLT1與FLT5的遺傳相似系數(shù)最大（0.90），遺傳差異最小?；贜TSYS軟件進行聚類分析構(gòu)建系統(tǒng)進化樹（圖5），結(jié)果表明，31個種質(zhì)被有效區(qū)分，且結(jié)果穩(wěn)定可重復，表明基于EST-SSR分子標記體系能對粗葉榕種質(zhì)資源進行遺傳多樣性分析，通過聚類分析能有效區(qū)分與鑒定不同種質(zhì)的遺傳多樣性。同時可構(gòu)建不同粗葉榕種質(zhì)指紋圖譜，為粗葉榕資源創(chuàng)建二維碼形式的DNA分子身份標識（圖5為部分種質(zhì)指紋圖譜）。

3""討論

本研究對4種不同葉型的粗葉榕葉片進行轉(zhuǎn)錄組測序，開發(fā)了SSR分子標記。研究顯示，4種葉型粗葉榕葉片轉(zhuǎn)錄組中的SSR位點總數(shù)為16"851個，出現(xiàn)頻率為49.02%，分布于12"094條unigene上，平均每2.64"kb有1個SSR位點。出現(xiàn)頻率高于同?？浦参锷洌?7.09%）[18]，并高于防風（15.51%）[19]、白三葉（13.25%）[20]、八角（18.85%）[21]等藥用植物，平均分布距離高于辣椒花（7.9"kb）[22]、苦蕎（3.04"kb）[23]、海帶（7.07"kb）[24]等植物，低于薏苡（0.20"kb）、三七（0.13"kb）等[25]。這種差異可能和物種之間SSR位點的表達差異有關(guān)，也可能是由于物種遺傳特性、SSR搜索標準的差異，數(shù)據(jù)量的大小以及選擇分析軟件的不同導致的[26-29]?？傮w來看粗葉榕轉(zhuǎn)錄組中SSR位點的出現(xiàn)頻率較高，平均分布距離較短，說明粗葉榕轉(zhuǎn)錄組中SSR位點的數(shù)量較多且種類豐富。

不同物種的SSR位點重復基元類型存在差異。主要體現(xiàn)在不同植物SSR位點數(shù)量、重復類型及占比有所不同。本研究中，粗葉榕轉(zhuǎn)錄組中的SSR基元類型豐富多樣，包括雙堿基到六堿基的重復，其中雙堿基和三堿基的重復占據(jù)主導地位，這2類基元類型構(gòu)成的SSR位點共計72.49%，而且以雙堿基重復位點數(shù)量最多，6次重復基元出現(xiàn)的頻率最高。蔣爽等[30]通過對枇杷轉(zhuǎn)錄組SSR重復類型進行分析，發(fā)現(xiàn)該轉(zhuǎn)錄組中的二核苷酸重復數(shù)量占總數(shù)的71.1%，以AT/TA為主，三核苷酸重復類型占15.0%，以AAG/CTT為主，雙堿基和三堿基重復為主要類型，共占總SSR位點的86.1%，這與本研究結(jié)果一致。本研究中雙堿基重復類型較多，占總SSR位點的50.65%，郭俊等[31]通過對油梨轉(zhuǎn)錄組SSR重復類型進行分析，也發(fā)現(xiàn)該轉(zhuǎn)錄組的重復以二核苷酸和三核苷酸為主，占SSR總數(shù)的87.15%。陳友吾等[32]對多花黃精轉(zhuǎn)錄組SSR位點分析發(fā)現(xiàn)，SSR位點中的主導類型是二核苷酸重復，占總SSR的53.14%，其次是三核苷酸重復，占比達33.31%，均與本研究結(jié)果相似。

評估分子標記的適用性需要考慮其多態(tài)性程度，而具有豐富多態(tài)性和高潛力的低級SSR重復基序是開發(fā)這些標記的關(guān)鍵條件[20，"28]。已有研究表明，當基序長度lt;12"bp時，SSR的多態(tài)性較低；而當基序長度≥12"bp時，隨著基序長度的增加，SSR的分布頻率則降低；基序長度為12～19"bp時，SSR的多態(tài)性顯著下降；基序長度≥20"bp時，多態(tài)性較高[33]。在本研究中，粗葉榕轉(zhuǎn)錄組中基序長度超過20"bp以上的SSR位點有9220個，占SSR總數(shù)的54.71%，這部分多態(tài)性較高的SSR具有良好的開發(fā)應用價值，大多可作為后續(xù)研究的基礎。隨機選擇100對引物進行驗證，結(jié)果顯示，89對引物可擴增出目標產(chǎn)物，有效性達89%；55對引物具有多態(tài)性，多態(tài)性引物比例達55%，說明該轉(zhuǎn)錄組多態(tài)性潛能高。55對多態(tài)性引物中挑選9對條帶清晰的引物進行引物遺傳多樣性分析，成功將來源不同的31個粗葉榕種質(zhì)進行合理區(qū)分并成功構(gòu)建指紋圖譜，且9對引物在31份種質(zhì)中多態(tài)性條帶比例高，均值達72.04%。

綜上所述，本研究對粗葉榕轉(zhuǎn)錄組測序后進行位點特征分析及分子標記開發(fā)，證實粗葉榕轉(zhuǎn)錄組SSR位點信息豐富，可利用性高。本研究結(jié)果可用于后續(xù)種質(zhì)資源評價及遺傳多樣性分析等研究，是對現(xiàn)有粗葉榕轉(zhuǎn)錄組數(shù)據(jù)的進一步發(fā)掘和利用，所獲得的分子標記對豐富榕屬植物的分子標記類型及開展特定性狀的輔助選擇具有重要意義。

參考文獻

[1]"錢潤，"周駿輝，"楊健，"黃璐琦，"袁媛."中藥材分子標記輔助育種技術(shù)研究進展[J]."中國中藥雜志，"2020，"45（20）："4812-4818.QIAN"R，"ZHOU"J"H，"YANG"J，"HUANG"L"Q，"YUAN"Y."Study"progress"on"molecular"marker-assisted"breeding"of"Chinese"medicinal"materials[J]."China"Journal"of"Chinese"Materia"Medica，"2020，"45（20）："4812-4818."（in"Chinese）

[2]"王娟娟，"趙明，"韓雨威，"吳蒙蒙，"劉瑞振，"沈超，"祁哲晨."微衛(wèi)星DNA標記開發(fā)技術(shù)進展及其在經(jīng)濟植物研究中的應用[J]."生命科學研究，"2016，"20（3）："260-266.WANG"J"J，"ZHAO"M，"HAN"Y"W，"WU"M"M，"LIU"R"Z，"SHEN"C，"QI"Z"C."Advances"in"development"of"microsatellite"DNA"markers"and"their"applications"in"economic"plant"research[J]."Life"Science"Research，"2016，"20（3）："260-266."（in"Chinese）

[3]"崔亞杰，"史荷鳳，"聶廣樓，"廖海迪，"李學剛，"龐蕾."黃連轉(zhuǎn)錄組SSR位點信息分析及分子標記開發(fā)[J]."中草藥，"2024，"55（19）："6713-6720.CUI"Y"J，"SHI"H"F，"NIE"G"L，"LIAO"H"D，"LI"X"G，"PANG"L."Site"information"analysis"and"molecular"marker"based"development"of"SSR"in"Coptis"chinensis"transcriptome[J]."Chinese"Traditional"and"Herbal"Drugs，"2024，"55（19）："6713-6720."（in"Chinese）

• 陸婧."中國榕屬綿毛榕組粗葉榕復合群的系統(tǒng)分類學研究[D]."上海："華東師范大學，"2015.LU"J."Phylogeny"and"taxonomy"of"the"‘hairy-fig’"complex"of"Ficus"sect."Eriosycea"（Moraceae）"in"China[D]."Shanghai："East"China"Normal"University，"2015."（in"Chinese）
• 馬雅靜，"劉煥，"史晶晶，"尚明英，"李耀利，"劉廣學，"徐風，"蔡少青."五指毛桃的質(zhì)量標準研究[J]."中草藥，"2017，"48（4）："782-791.MA"Y"J，"LIU"H，"SHI"J"J，"SHANG"M"Y，"LI"Y"L，"LIU"G"X，"XU"F，"CAI"S"Q."Study"on"quality"standard"of"Fici"Hirtae"Radix[J]."Chinese"Traditional"and"Herbal"Drugs，"2017，"48（4）："782-791."（in"Chinese）
• LU"Y，"CHEN"J"L，"CHEN"B，"LIU"Q"Q，"ZHANG"H"L，"YANG"L"Y，"CHAO"Z，"TIAN"E"W."High"genetic"diversity"and"low"population"differentiation"of"a"medical"plant"Ficus"hirta"Vahl.，"uncovered"by"microsatellite"loci："implications"for"conservation"and"breeding[J]."BMC"Plant"Biology，"2022，"22（1）："334.
• 楊麗薇，"杜麗敏，"林旭俊，"李善志，"李曉斌."粗葉榕種質(zhì)資源調(diào)查與收集淺析[J]."熱帶林業(yè)，nbsp;2021，"49（4）："54-56.YANG"L"W，"DU"L"M，"LIN"X"J，"LI"S"Z，"LI"X"B."Investigation"and"collection"of"Ficus"hirta"Vahl."resources[J]."Tropical"Forestry，"2021，"49（4）："54-56."（in"Chinese）
• 杜麗敏，"林旭俊，"楊麗薇，"黃國寧，"李曉斌，"郭霞."不同表型粗葉榕根部營養(yǎng)成分特征分析[J]."分子植物育種，"2024，"22（18）："6138-6145.DU"L"M，"LIN"X"J，"YANG"L"W，"HUANG"G"N，"LI"X"B，"GUO"X."Analysis"of"nutrient"components"in"roots"of"different"phenotypes"of"Ficus"hirta"Vahl.[J]."Molecular"Plant"Breeding，"2024，"22（18）："6138-6145."（in"Chinese）
• 鄭蓉蓉，"軋霽，"王文婧，"楊寒冰，"張慶文，"張曉琦，"葉文才."五指毛桃的化學成分研究[J]."中國中藥雜志，"2013，"38（21）："3696-3701.ZHENG"R"R，"YA"J，"WANG"W"J，"YANG"H"B，"ZHANG"Q"W，"ZHANG"X"Q，"YE"W"C."Chemical"studies"on"roots"of"Ficus"hirta[J]."China"Journal"of"Chinese"Materia"Medica，"2013，"38（21）："3696-3701."（in"Chinese）
• 易明濤，"李明蘭."粗葉榕組織培養(yǎng)常見問題及對策[J]."熱帶農(nóng)業(yè)工程，"2024，"48（2）："117-119.YI"M"T，"LI"M"L."Analysis"on"the"common"problems"and"countermeasures"of"Ficus"Hirta"Vahl."tissue"culture[J]."Tropical"Agricultural"Encineering，"2024，"48（2）："117-119."（in"Chinese）

[11]"黃國婷，"張青槐，"揭潔，"周刊."基于壯醫(yī)理論探討壯藥五指毛桃的藥理作用及臨床應用[J]."中國民族民間醫(yī)藥，"2024，"33（17）："61-66.HUANG"G"T，"ZHANG"Q"H，"JIE"J，"ZHOU"K."Based"on"the"theory"of"Zhuang"medicine，"the"pharmacological"effects"and"clinical"application"of"Zhuang"medicine"Ficus"hirta"Vahl."were"discussed[J]."Chinese"Journal"of"Ethnomedicine"and"Ethnopharmacy，"2024，"33（17）："61-66."（in"Chinese）

[12]"GAO"W"R，"MA"R"H，"LI"X，"LIU"J"Q，"JIANG"A"H，"TAN"P"T，"XIONG"G"X，"DU"C"Z，"ZHANG"J"J，"ZHANG"X"C，"FANG"X"M，"YI"Z"L，"ZHANG"J."Construction"of"genetic"map"and"QTL"mapping"for"seed"size"and"quality"traits"in"soybean"（Glycine"max"L.）[J]."International"Journal"of"Molecular"Sciences，"2024，"25（5）："2857

[13]"PATEL"R，"MEMON"J，"KUMAR"S，"PATEL"D"A，"SAKURE"A"A，"PATEL"M"B，"DAS"A，"KARJAGI"C"G，"PATEL"S，"PATEL"U，"ROYCHOWDHURY"R."Genetic"diversity"and"population"structure"of"maize"（Zea"mays"L.）"inbred"lines"in"association"with"phenotypic"and"grain"qualitative"traits"using"SSR"genotyping[J]."Plants"（Basel），"2024，"13（6）："823.

[14]"劉樹芳，"張先聞，"董麗英，"劉沛，"田維逵，"周伍民，"楊勤忠."云南不同稻區(qū)稻瘟病菌群體遺傳結(jié)構(gòu)的SSR分析[J]."中國農(nóng)學通報，"2023，"39（12）："123-130.LIU"S"F，"ZHANG"X"W，"DONG"L"Y，"LIU"P，"TIAN"W"K，"ZHOU"W"M，"YANG"Q"Z."SSR"analysis"of"population"genetic"structure"of"Magnaporthe"oryzae"from"different"rice-growing"regions"in"Yunnan[J]."Chinese"Agricultural"Science"Bulletin，"2023，"39（12）："123-130."（in"Chinese）

[15]"陳榮珠，"王小平，"高偉城，"羅曉鋒，"周建金，"林美珍."4種裂葉類型五指毛桃果實中氨基酸組成比較及營養(yǎng)價值評價[J]."中國食品添加劑，"2024，"35（1）："238-246.CHEN"R"Z，"WANG"X"P，"GAO"W"C，"LUO"X"F，"ZHOU"J"J，"LIN"M"Z."Comparative"analysis"of"amino"acids"composition"in"fruits"of"four"types"Ficus"Hirta"with"lobed"leaf"and"evaluation"of"their"nutritional"value[J]."China"Food"Additives，"2024，"35（1）："238-246."（in"Chinese）

[16]"胡超瓊."基于SCoT、IRAP和SSR分子標記的柿屬植物親緣關(guān)系分析[D]."楊陵："西北農(nóng)林科技大學，"2018.HU"C"Q."Genetic"relationship"analysis"of"persimmon"by"SCoT，"IRAP"and"SSR"molecular"markers[D]."Yangling："Northwest"A"amp;"F"University，"2018."（in"Chinese）

[17]"趙秋菊，"崔學強，"鄧杰玲，"黃昌艷，"歐美景，"李紅燕，"張自斌."金釵石斛居群遺傳多樣性分析及DNA指紋圖譜構(gòu)建[J]."中草藥，"2023，"54（19）："6434-6442.ZHAO"Q"J，"CUI"X"Q，"DENG"J"L，"HUANG"C"Y，"OU"M"J，"LI"H"Y，"ZHANG"Z"B."Genetic"diversity"analysis"and"DNA"fingerprint"construction"of"Dendrobium"nobile"populations[J]."Chinese"Traditional"and"Herbal"Drugs，"2023，"54（19）："6434-"6442."（in"Chinese）

[18]"姚劉慧，"唐翠明，"戴凡煒，"王振江，"羅國慶."基于桑樹轉(zhuǎn)錄組測序的SSR標記開發(fā)和引物篩選[J]."蠶業(yè)科學，"2015，"41（1）："18-27.YAO"L"H，"TANG"C"M，"DAI"F"W，"WANG"Z"J，"LUO"G"Q."Development"and"primer"screening"of"SSRnbsp;markers"based"on"mulberry"transcriptome"sequencing[J]."Science"of"Sericulture，"2015，"41（1）："18-27."（in"Chinese）

[19]"李爽，"寇佩雯，"許祎珂，"劉長樂，"孫曉春，"李鉑，"黃文靜."基于轉(zhuǎn)錄組的防風SSR分子標記開發(fā)及應用[J]."中國中藥雜志，"2024，"49（3）："717-727.LI"S，"KOU"P"W，"XU"Y"K，"LIU"C"L，"SUN"X"C，"LI"B，"HUANG"W"J."Development"and"application"of"SSR"markers"of"Saposhnikovia"divaricata"based"on"transcriptome[J]."China"Journal"of"Chinese"Materia"Medica，"2024，"49（3）："717-727."（in"Chinese）

[20]"張婷婷，"張鶴山，"宋康杰，"趙澤宇，"許本波，"劉洋."白三葉轉(zhuǎn)錄組SSR位點特征分析及引物開發(fā)[J]."草業(yè)科學，"2023，"40（9）："2266-2275.ZHANG"T"T，"ZHANG"H"S，"SONG"K"J，"ZHAO"Z"Y，"XU"B"B，"LIU"Y."Analysis"of"SSR"site"characteristics"of"Trifolium"repens"transcriptome"and"primer"development[J]."Pratacultural"Science，"2023，"40（9）："2266-2275."（in"Chinese）

[21]"余興華，"楊艷梅，"李玉祥，"陶永宏，"孟凡來，"滕娟."基于八角轉(zhuǎn)錄組序列的SSR分子標記引物開發(fā)[J]."西南農(nóng)業(yè)學報，"2022，"35（6）："1261-1268.YU"X"H，"YANG"Y"M，"LI"Y"X，"TAO"Y"H，"MENG"F"L，"TENG"J."Development"of"SSR"molecular"markers"based"on"transcriptome"sequencing"of"Illicium"verum[J]."Southwest"China"Journal"of"Agricultural"Sciences，"2022，"35（6）："1261-"1268."（in"Chinese）

[22]"聶智星，"傅鴻妃，"李高青，"郭賽賽，"鄭積榮."辣椒花藥轉(zhuǎn)錄組SSR位點信息分析及標記開發(fā)[J].nbsp;分子植物育種，"2023，"21（23）："7803-7810.NIE"Z"X，"FU"H"F，"LI"G"Q，"GUO"S"S，"ZHENG"J"R."SSR"site"information"analysis"and"marker"development"of"pepper"（Cap"sicum"annuum"L.）"anther"transcriptome[J]."Molecular"Plant"Breeding，"2023，"21（23）："7803-7810."（in"Chinese）

[23]"蔡齊宗，"王佳蕊，"陳慶富，"李洪有."苦蕎全基因組SSR位點鑒定及分子標記開發(fā)[J]."河南農(nóng)業(yè)大學學報，"2022，"56（3）："392-400.CAI"Q"Z，"WANG"J"R，"CHEN"Q"F，"LI"H"Y."Identification"of"SSR"loci"and"molecular"marker"development"in"the"whole"genome"of"tartary"buckwheat[J]."Journal"of"Henan"Agricultural"University，"2022，"56（3）："392-400."（in"Chinese）

[24]"林錦彬，"伍永久，"韓藝娟，"陳松彪，"唐唯其."海帶全基因組SSR分子標記開發(fā)[J]."分子植物育種，"2023，"21（12）："3964-3972.LIN"J"B，"WU"Y"J，"HAN"Y"J，"CHEN"S"B，"TANG"W"Q."Genome-wide"development"of"SSR"markers"in"Saccharina"japonica[J]."Molecular"Plant"Breeding，"2023，"21（12）："3964-"3972."（in"Chinese）

[25]"揭應碧，"盧迎春，"宋婉玲，"韋坤華，"張廣輝，"楊生超."三七基因組SSR位點分析和多態(tài)性引物開發(fā)[J]."熱帶作物學報，"2018，"39（7）："1351-1358.JIE"Y"B，"LU"Y"C，"SONG"W"L，"WEI"K"H，"ZHANG"G"H，"YANG"S"C."Analysis"of"SSR"loci"information"in"Panax"noto g inseng"genomes"and"polymorphic"primers"development[J]."Chinese"Journal"of"Tropical"Crops，"2018，"39（7）："1351-1358."（in"Chinese）

[26]"陳巍."三葉木通微衛(wèi)星序列分析及其在核心種質(zhì)構(gòu)建中的應用[D]."雅安："四川農(nóng)業(yè)大學，"2023.CHEN"W."Microsatellite"sequence"analysis"of"Akebia"trifoliata"and"application"in"core"collection"construction[D]."Ya’an："Sichuan"Agriculture"University，"2023."（in"Chinese）

[27]"王彬彬，"高妍夏，"郭欣慰，"紀宏亮，"王敬，"王暉."基于高通量測序的北柴胡根轉(zhuǎn)錄組SSR位點信息分析[J]."西北農(nóng)業(yè)學報，"2024，"33（4）："707-717.WANG"B"B，"GAO"Y"X，"GUO"X"W，"JI"H"L，"WANG"J，"WANG"H."SSR"loci"analysis"using"high-throughput"transcrip tome"sequencing"in"roots"of"Bupleurum"chinense"DC.[J]."Acta"Agriculturae"Boreali-occidentalis"Sinica，"2024，"33（4）："707-717."（in"Chinese）

[28]"周勝芳."石斛SSR分子標記的開發(fā)及應用研究[D]."南充："西華師范大學，"2023.ZHOU"S"F."Development"and"application"of"SSR"molecular"markers"for"the"Dendrobium[D]."Nanchong："China"West"Normal"University，"2023."（in"Chinese）

[29]"劉強，"趙麗，"郭虹霞，"張麗光，"鄧妍，"王陸軍，"王創(chuàng)云."玉米野生近緣種大芻草幼苗葉片的轉(zhuǎn)錄組SSR分析[J]."玉米科學，"2020，"28（2）："69-75.LIU"Q，"ZHAO"L，"GUO"H"X，"ZHANG"L"G，"DENG"Y，"WANG"L"J，"WANG"C"Y."Transcriptome"SSR"analysis"of"seedling"leaves"in"maize"wild"relative"species"teosinte[J]."Journal"of"Maize"Sciences，"2020，"28（2）："69-75."（in"Chinese）

[30]"蔣爽，"張學英，"安海山，"徐芳杰，"章加應."枇杷全基因組SSR標記開發(fā)及其多態(tài)性研究[J]."園藝學報，"2021，"48（5）："1013-1022.JIANG"S，"ZHANG"X"Y，"AN"H"S，"XU"F"J，"ZHANG"J"Y."Development"and"analysis"of"polymorphism"of"SSR"markers"in"the"whole"genome"of"loquat[J]."Acta"Horticulturae"Sinica，"2021，"48（5）："1013-1022."（in"Chinese）

[31]"郭俊，"朱婕，"謝尚潛，"張葉，"葉蓓蕾，"鄭麗燕，"凌鵬."油梨轉(zhuǎn)錄組SSR分子標記開發(fā)與種質(zhì)資源親緣關(guān)系分析[J]."園藝學報，"2020，"47（8）："1552-1564.GUO"J，"ZHU"J，"XIE"S"Q，"ZHANG"Y，"YE"B"L，"ZHENG"L"Y，"LING"P."Development"of"SSR"molecular"markers"based"on"transcriptome"and"analysis"of"genetic"relationship"of"germplasm"resources"in"avocado[J]."Acta"Horticulturae"Sinica，"2020，"47（8）："1552-1564."（in"Chinese）

[32]"陳友吾，"廖榮俊，"葉碧歡，"宋其巖，"楊陽，"李楠，"李海波."多花黃精轉(zhuǎn)錄組SSR位點分析及分子標記開發(fā)[J]."中草藥，"2020，"51（1）："182-189.CHEN"Y"W，"LIAO"R"J，"YE"B"H，"SONG"Q"Y，"YANG"Y，"LI"N，"LI"H"B."Analysis"on"SSR"loci"in"transcriptome"and"development"of"molecular"markers"in"Polygonatum"cyrtonema[J]."Chinese"Traditional"and"Herbal"Drugs，"2020，"51（1）："182-189."（in"Chinese）

[33]"TEMNYKH"S，"DECLERCK"G，"LUKASHOVA"A，"LIPOVI CH"L，"CARTINHOUR"S，"MCCOUCH"S."Computational"and"experimental"analysis"of"microsatellites"in"rice"（Oryza"sativa"L.）："frequency，"length"variation，"transposon"associations，"and"genetic"marker"potential[J]."Genome"Research，"2001，"11（8）："1441-1452.