辛靜 辛雅萱 董章宏 瞿紹宏 姚國(guó)瓊 李鳳 屈亞亞 唐軍榮

摘要：【目的】分析云南火焰蘭轉(zhuǎn)錄組中SSR分布及其序列分布特征，為大量開(kāi)發(fā)EST-SSR引物及開(kāi)展火焰蘭屬乃至蘭科植物的SSR遺傳多樣性分析、種質(zhì)資源評(píng)價(jià)、遺傳圖譜構(gòu)建及遺傳育種提供理論依據(jù)。【方法】以云南火焰蘭無(wú)菌苗幼嫩葉片為材料，利用高通量測(cè)序技術(shù)進(jìn)行轉(zhuǎn)錄組測(cè)序，使用Trinity對(duì)測(cè)序結(jié)果進(jìn)行de novo組裝，并以MISA對(duì)獲得的Unigene進(jìn)行SSR位點(diǎn)搜索，然后利用Excel 2010對(duì)云南火焰蘭轉(zhuǎn)錄組中SSR出現(xiàn)頻率、平均分布距離、基元類型及其重復(fù)類型組成等進(jìn)行統(tǒng)計(jì)分析?！窘Y(jié)果】從77888條去冗余的Unigenes序列中搜索到5051個(gè)SSR位點(diǎn)，其中有414個(gè)屬于復(fù)合型SSR位點(diǎn)，實(shí)際以4637個(gè)SSR位點(diǎn)進(jìn)行分析，SSR出現(xiàn)頻率為5.95%，平均分布距離13.53 kb。SSR位點(diǎn)中，二核苷酸為主要重復(fù)基元類型，數(shù)量最多，占SSR總數(shù)的52.21%，其次是三核苷酸重復(fù)基元，占39.42%，四核苷酸重復(fù)基元～六核苷酸重復(fù)基元數(shù)量較少。其中，二核苷酸重復(fù)基元中的主要重復(fù)類型為CT/GA，占SSR總數(shù)的20.64%，三核苷酸重復(fù)基元的主要重復(fù)類型為CTT/GAA，占SSR總數(shù)的4.70%。SSR基元各重復(fù)類型的重復(fù)次數(shù)以5～8次居多，占SSR總數(shù)的73.35%，其中，占比最多的是6次重復(fù)，占SSR總數(shù)的26.98%，其次是5次和7次，分別占SSR總數(shù)的23.08%和14.02%，且SSR數(shù)量隨各基元重復(fù)次數(shù)的增加而降低。云南火焰蘭轉(zhuǎn)錄組中SSR基序的平均長(zhǎng)度是18.28 bp，其中大小為12～20 bp的SSR基序占SSR總數(shù)的77.49%，21～30 bp基序占SSR總數(shù)的17.88%，大于30 bp的基序占SSR總數(shù)的4.63%，且SSR分布頻率隨基序長(zhǎng)度的增加而下降?！窘Y(jié)論】云南火焰蘭轉(zhuǎn)錄組SSR中低級(jí)基元類型較豐富，開(kāi)發(fā)出多態(tài)性高的SSR引物潛力較高，可用于云南火焰蘭乃至火焰蘭屬植物的遺傳多樣性及種質(zhì)資源評(píng)價(jià)、遺傳圖譜構(gòu)建及遺傳育種等研究。

關(guān)鍵詞： 云南火焰蘭;轉(zhuǎn)錄組;SSR;分布;序列特征;重復(fù)基元

中圖分類號(hào)： S682.31? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：2095-1191（2020）07-1634-08

Abstract：【Objective】The analysis of SSR distribution and sequence distribution characteristics in the transcriptome of Renanthera imschootiana Rolfe provided a theoretical basis for the development of? EST-SSR primers， SSR genetic diversity analysis， germplasm resource evaluation， genetic map construction and genetic breeding of Renanthera and even Orchidaceae. 【Method】The young leaves of aseptic R. imschootiana seedlings were used as material to perform transcriptome sequencing by high throughput sequencing. The de novo assembly of the sequencing results was conducted by Trinity， SSR loci search was performed using MISA. Excel 2010 was used to analyze the occurrence frequency， average distribution distance， primitive type and the composition of repetition of SSR in the transcriptome of R.imschootiana. 【Result】A total of 5051 SSR loci were searched from 77888 remove redundant Unigenes sequences， among them， 414 were compound SSR loci， then 4637 SSR loci were used as actual research object. SSR occurrence frequency was 5.95% and avera-ge distribution distance was 13.53 kb. In these SSR loci， dinucleotide repeat units were the main type， accounting for 52.21%， followed by the trinucleotide repeat units（39.42%）. The proportion of tetranucleotide repeat units to hexanucleotide repeat units was low. The main repeat motif types in dinucleotide was CT/GA ， accounting for 20.64% of the total SSR， the main repeat motif types in trinucleotide was CTT/GAA， accounting for 4.70% of the total SSR. The types of SSR repeat units were more of 5-8 repetition times， accounting for 73.35% of the total SSR. The proportion of 6-repetition was the largest with 26.98% proportion of the total SSR， followed by 5-repetition（23.08%） and 7-repetition（14.02%）. As the number of primitive repeats increased， SSR reduced. The average length of SSR sequence motif in the transcriptome of R.imschootiana was 18.28 bp， among them， the SSR sequence motif of 12-20 bp was accounting for 77.49% in the total SSR， the SSR motif of 21-30 bp was accounting for 17.88%， and SSR motif over 30 bp was accounting for 4.63% in the total SSR the distribution frequency of SSR reduced with increasing of motif length. 【Conclusion】The transcriptome SSR of R.imschootiana is rich in low-level primitive types and has a high potential to develop SSR primers with high polymorphism， which can be used to study the genetic diversity and germplasm resources evaluation， genetic map construction and genetic breeding of R.imschootiana and even the Renanthera plants.

Key words： Renanthera imschootiana Rolfe; transcriptome; SSR; distribution; sequence characteristics; repeat primitive

Foundation item： Joint Special Project of Agricultural Basic Research in Yunnan[2017FG001（-019）]

0 引言

【研究意義】云南火焰蘭（Renanthera imschoo-tiana Rolfe）為蘭科（Orchidaceae）火焰蘭屬（Renanthera）附生植物，因其花色紅而艷麗，猶如火焰而得名，主要生長(zhǎng)于云南元江縣（吳征鎰和洪德元，1999;中國(guó)科學(xué)院昆明植物研究所，2003）。野生云南火焰蘭的分布范圍狹窄、數(shù)量極少，且近年來(lái)遭到嚴(yán)重采挖，致使其野生種質(zhì)資源瀕臨滅絕，現(xiàn)已被列入《中國(guó)生物多樣性紅色名錄——高等植物卷》的極危（CR）物種和《瀕危野生動(dòng)植物種國(guó)際貿(mào)易公約》附錄I物種（李娟等，2018）。利用高通量測(cè)序數(shù)據(jù)，統(tǒng)計(jì)并分析云南火焰蘭SSR分布及其特征，以期開(kāi)發(fā)出大量SSR引物，對(duì)云南火焰蘭種質(zhì)資源評(píng)價(jià)、保護(hù)利用及遺傳育種具有重要意義。【前人研究進(jìn)展】目前，有關(guān)云南火焰蘭的研究報(bào)道主要涉及離體快速繁殖（林丹妮等，2006）、種子無(wú)菌萌發(fā)和試管育苗（林丹妮等，2008）、雜交育種（曹雯靜等，2014;陳和明等，2019）、種質(zhì)資源評(píng)價(jià)（姚婧，2015）及種群現(xiàn)狀分析（李娟等，2018）等，而有關(guān)其分子生物學(xué)方面的研究鮮見(jiàn)報(bào)道。SSR分子標(biāo)記具有多態(tài)性高、標(biāo)記數(shù)量多、重復(fù)性好、覆蓋面廣等優(yōu)點(diǎn)，已廣泛應(yīng)用于分子標(biāo)記輔助選擇、遺傳多樣性分析和遺傳圖譜構(gòu)建等領(lǐng)域（王大瑋等，2018;張金旺等，2018）。根據(jù)SSR的來(lái)源可將其分為基因組SSR和EST-SSR兩大類型（程小毛和黃曉霞，2011），前者是利用cDNA文庫(kù)構(gòu)建和基因組測(cè)序的方法獲得，該方法步驟繁瑣、工作量大且成本高;后者則通過(guò)搜索公共數(shù)據(jù)庫(kù)中的EST（表達(dá)序列標(biāo)簽）或轉(zhuǎn)錄組測(cè)序數(shù)據(jù)進(jìn)行開(kāi)發(fā)，需要大量EST數(shù)據(jù)支持（趙能等，2017）。葉鵬等（2019a，2019b）先后對(duì)藥用植物金鐵鎖（Psammosilene tunicoides）和云南金花茶（Camellia fascicularis）轉(zhuǎn)錄組SSR位點(diǎn)信息及序列特征進(jìn)行分析，結(jié)果發(fā)現(xiàn)二者的轉(zhuǎn)錄組中SSR出現(xiàn)頻率均較高，基元類型較豐富，重復(fù)頻次也較高，從而表現(xiàn)出高的多態(tài)性潛力。此外，研究人員還對(duì)云南松（Pinus yunnanensis）（蔡年輝等，2015）、杜仲（Eucommia ulmoides）（申響保等，2016）、四球茶（Camellia tetracocca）（黎瑞源等，2017）、木荷（Schima superba）（林艷等，2018）、牛樟（Cinnamomum kanehirae）（郭鶯等，2018）、芒果（Mangifera indica）（唐玉娟等，2018）、三七（Panax notoginseng）（揭應(yīng)碧等，2018）和牡丹（賀丹等，2019）等植物進(jìn)行轉(zhuǎn)錄組SSR位點(diǎn)信息分析，為深入開(kāi)展植物遺傳結(jié)構(gòu)和遺傳多樣性及分子標(biāo)記輔助育種提供了前期基礎(chǔ)數(shù)據(jù)?！颈狙芯壳腥朦c(diǎn)】云南火焰蘭在分子水平的研究基礎(chǔ)非常薄弱，至今未見(jiàn)有關(guān)云南火焰蘭轉(zhuǎn)錄組SSR分布及其序列特征的研究報(bào)道。【擬解決的問(wèn)題】對(duì)云南火焰蘭進(jìn)行轉(zhuǎn)錄組測(cè)序，并分析轉(zhuǎn)錄組中SSR數(shù)量、出現(xiàn)頻率及分布距離等，明確云南火焰蘭高多態(tài)性SSR序列的分布特征，為大量開(kāi)發(fā)EST-SSR引物及開(kāi)展火焰蘭屬乃至蘭科植物的SSR遺傳多樣性分析、種質(zhì)資源評(píng)價(jià)、遺傳圖譜構(gòu)建及遺傳育種等提供理論依據(jù)。

1 材料與方法

1. 1 試驗(yàn)材料

供試云南火焰蘭種植于云南元江縣，采集長(zhǎng)勢(shì)好且無(wú)病蟲(chóng)害云南火焰蘭母株上的果莢，采用種子無(wú)菌萌發(fā)成苗后，取其幼嫩葉片，迅速放入液氮中保存?zhèn)溆谩?/p>

1. 2 轉(zhuǎn)錄組測(cè)序

將云南火焰蘭幼嫩葉片采用干冰保存后送至廣州基迪奧生物科技有限公司，后續(xù)RNA提取、文庫(kù)構(gòu)建及測(cè)序工作均由該公司完成，測(cè)序平臺(tái)為Illumina HiSeqTM 4000。

1. 3 序列拼接及SSR位點(diǎn)搜索

使用Trinity對(duì)云南火焰蘭轉(zhuǎn)錄組測(cè)序結(jié)果進(jìn)行de novo組裝，并使用MISA（http：//pgrc.ipk-gatersleben.de/misa/）對(duì)轉(zhuǎn)錄組中所有的Unigenes進(jìn)行SSR位點(diǎn)搜索，參數(shù)設(shè)為二核苷酸重復(fù)次數(shù)至少6次，三核苷酸重復(fù)次數(shù)至少5次，四核苷酸～六核苷酸重復(fù)次數(shù)至少4次。

1. 4 統(tǒng)計(jì)分析

利用Excel 2010對(duì)云南火焰蘭轉(zhuǎn)錄組中SSR出現(xiàn)頻率、平均分布距離、基元類型及其重復(fù)類型組成等進(jìn)行統(tǒng)計(jì)分析（張震等，2018;楊帥等，2020），相關(guān)計(jì)算公式如下：

SSR出現(xiàn)頻率=SSR數(shù)量/Unigene總數(shù)

SSR位點(diǎn)平均分布距離=SSR總序列長(zhǎng)度/SSR數(shù)量

2 結(jié)果與分析

2. 1 云南火焰蘭轉(zhuǎn)錄組中SSR出現(xiàn)頻率和分布距離

通過(guò)對(duì)云南火焰蘭轉(zhuǎn)錄組進(jìn)行組裝，共獲得77888條去冗余的Unigenes，長(zhǎng)度為201～16505 bp，平均長(zhǎng)度805 bp，總長(zhǎng)度62740868 bp，GC含量42.73%。從中搜索出5051個(gè)SSR位點(diǎn)，其中，有414個(gè)屬于復(fù)合型SSR位點(diǎn)，實(shí)際以4637個(gè)SSR位點(diǎn)進(jìn)行分析。SSR出現(xiàn)頻率為5.95%，平均分布距離13.53 kb（表1），即云南火焰蘭轉(zhuǎn)錄組中平均13.53 kb就出現(xiàn)1個(gè)SSR位點(diǎn)，說(shuō)明SSR較豐富，為后期SSR引物開(kāi)發(fā)提供了大量的基礎(chǔ)數(shù)據(jù)。

2. 2 云南火焰蘭轉(zhuǎn)錄組中SSR基元類型

由表1可知，SSR位點(diǎn)中，二核苷酸為主要重復(fù)基元類型，數(shù)量最多，占SSR總數(shù)的52.21%，其次為三核苷酸重復(fù)基元，占SSR總數(shù)的39.42%;四核苷酸～六核苷酸重復(fù)基元數(shù)量相對(duì)較少;在SSR的平均分布距離方面，以五核苷酸重復(fù)基元最高，為1206.56 kb，以二核苷酸重復(fù)基元最低，僅為25.92 kb;在重復(fù)基元類型的出現(xiàn)頻率方面，以二核苷酸重復(fù)基元最高，為3.11%;三核苷酸重復(fù)基元次之，為2.35%，以五核苷酸重復(fù)基元最低，僅占0.07%?？梢?jiàn)，云南火焰蘭存在較多種類型的SSR基元，是其SSR引物開(kāi)發(fā)的重要保證。

2. 3 云南火焰蘭轉(zhuǎn)錄組中SSR堿基組成及比例

云南火焰蘭轉(zhuǎn)錄組中的SSR堿基組成及比例如表2所示。云南火焰蘭轉(zhuǎn)錄組SSR中，二核苷酸～六核苷酸重復(fù)基元的重復(fù)類型分別有6、30、66、43和105種，共有250種。二核苷酸中重復(fù)基元各基序所占比例排序?yàn)镃T/GA（957個(gè)，占SSR總數(shù)的20.64%）>TC/AG（851個(gè)，占SSR總數(shù)的18.35%）>AC/TG（241個(gè)，占SSR總數(shù)的5.20%）>CA/GT（188個(gè)，占SSR總數(shù)的4.05%）>AT/TA（172個(gè)，占SSR總數(shù)的3.71%）>CG/GC（12個(gè)，占SSR總數(shù)的0.26%），分別占雙堿基重復(fù)類型的比例為39.53%、35.15%、9.95%、7.77%、7.10%和0.50%（圖1-A）。三核苷酸重復(fù)基元最多的為CTT/GAA（218個(gè)），占SSR總數(shù)的4.70%，占三堿基重復(fù)類型的11.93%;其次是AAG/TTC（178個(gè)），占SSR總數(shù)的3.84%，占三堿基重復(fù)類型的9.74%;以ATA/TAT最少（6個(gè)），占SSR總數(shù)的0.13%，占三堿基重復(fù)類型的0.33%（圖1-B）。四核苷酸重復(fù)基元最多的為AAAT/TTTA（20個(gè)），占SSR總數(shù)的0.43%，占四堿基重復(fù)類型的10.10%（圖1-C）。五核苷酸重復(fù)基元最多的是ATGCA/TACGT（2個(gè)）和CCCTC/GGGAG（2個(gè)），均占SSR總數(shù)的0.04%，占五堿基重復(fù)類型的3.85%（圖1-D）。六核苷酸重復(fù)基元最多的是CCCTCT/GGGAGA（7個(gè)），占SSR總數(shù)的0.15%，占六堿基重復(fù)類型的5.07%（圖1-E）。綜上所述，云南火焰蘭轉(zhuǎn)錄組SSR堿基組成中，低級(jí)基元類型較豐富，可在該種或該屬植物中檢出較多的遺傳多態(tài)性。

2. 4 云南火焰蘭轉(zhuǎn)錄組中各基元重復(fù)次數(shù)

SSR基元重復(fù)類型的重復(fù)次數(shù)以5～8次居多，共計(jì)3401個(gè)，占SSR總數(shù)的73.35%，其中，占比最多的是6次重復(fù)，共計(jì)1251個(gè)，占SSR總數(shù)的26.98%;其次是5次和7次，分別為1070個(gè)和650個(gè)，占SSR總數(shù)的23.08%和14.02%。三核苷酸重復(fù)基元中重復(fù)類型最多的是5次，共993個(gè)，且SSR數(shù)量隨各基元重復(fù)次數(shù)的增加而降低（圖2）。

2. 5 云南火焰蘭轉(zhuǎn)錄組中SSR基序大小

云南火焰蘭轉(zhuǎn)錄組中大小為12～20 bp的SSR基序共計(jì)3593個(gè)，占SSR總數(shù)的77.49%，說(shuō)明從該轉(zhuǎn)錄組信息開(kāi)發(fā)出多態(tài)性高的SSR引物潛力較高。21～30 bp基序有829個(gè)，占SSR總數(shù)的17.88%;大于30 bp的基序有215個(gè)，占SSR總數(shù)的4.63%。SSR基序的平均長(zhǎng)度是18.28 bp，二核苷酸～六核苷酸的平均長(zhǎng)度分別是18.05、17.84、18.06、21.06和26.91 bp（表1），且SSR的分布頻率隨著基序長(zhǎng)度的增加而呈下降趨勢(shì)（圖3）。

3 討論

本研究對(duì)云南火焰蘭轉(zhuǎn)錄組測(cè)序獲得的77888條Unigenes進(jìn)行SSR位點(diǎn)搜索，共獲得二核苷酸重復(fù)基元～六核苷酸重復(fù)基元的SSR位點(diǎn)4637個(gè)，SSR出現(xiàn)頻率5.95%，平均分布距離13.53 kb。前人研究發(fā)現(xiàn)，從墨蘭（Cymbidium sinense）轉(zhuǎn)錄組測(cè)序獲得的51764條Unigenes中搜索到4419個(gè)SSR位點(diǎn)（陳起馨，2016）;從金釵石斛（Dendrobium nobile）轉(zhuǎn)錄組測(cè)序獲得的207283條Unigenes中搜索到32709個(gè)SSR位點(diǎn)（李清等，2017）;從寒蘭（Cymbidium kanran）轉(zhuǎn)錄組測(cè)序獲得的68699條Unigenes中搜索到9837個(gè)SSR位點(diǎn)（張亞楠等，2017）?？梢?jiàn)，云南火焰蘭的SSR位點(diǎn)數(shù)量與其親緣關(guān)系較近的物種均存在一定差異。就SSR出現(xiàn)頻率而言，本研究云南火焰蘭SSR出現(xiàn)頻率（5.95%）高于同科的鐵皮石斛（Dendrobium officinale）（5.11%）（肖冬長(zhǎng)等，2012），卻低于同科的金釵石斛（15.78%）（李清等，2017）、寒蘭（14.32%）（張亞楠等，2017）和天麻（Gastrodia elata）（8.16%）（周天華等，2017）的SSR出現(xiàn)頻率。但就平均分布距離來(lái)看，本研究云南火焰蘭SSR平均分布距離（13.53 kb）高于同科的墨蘭（11.45kb）（陳起馨，2016）、金釵石斛（3.75 kb）（李清等，2017）、寒蘭（6.25 kb）（張亞楠等，2017）和白及（Bletilla striata）（4.35 kb）（陳紅波，2019），與燈盞花（Erigeron breviscapus）（15.32 kb）（陳茵等，2014）和馬尾松（Pinus massoniana）（13.83 kb）（陳曉明等，2018）相當(dāng)?？梢?jiàn)，云南火焰蘭轉(zhuǎn)錄組中SSR出現(xiàn)頻率及其分布距離處于中等水平，且不同物種轉(zhuǎn)錄組測(cè)序信息存在明顯差異，一方面可能是其本身基因組差異的客觀表現(xiàn);另一方面可能是測(cè)序組織器官、轉(zhuǎn)錄組中Unigene數(shù)量和長(zhǎng)度、SSR搜索工具及設(shè)定篩選參數(shù)等差異引起（李俊仁等，2018）。

二核苷酸和三核苷酸是多種植物轉(zhuǎn)錄組中SSR的主要基元類型，但各基元類型占比有所差異（李翠婷等，2014;葉鵬等，2019a，2019b）。本研究中SSR以二核苷酸重復(fù)基元居多，占SSR總數(shù)的52.21%;其次是三核苷酸重復(fù)基元，占39.42%。在云南火焰蘭轉(zhuǎn)錄組中，二核苷酸重復(fù)基元以CT/GA為主，占SSR總數(shù)的20.64%;三核苷酸重復(fù)基元中，以CTT/GAA為主，占SSR總數(shù)的4.70%。這與蝴蝶蘭（Phalaenopsis spp.）（李冬梅等，2011）、石斛（Dendrobium）（邱道壽等，2013）和天麻（周天華等，2017）等植物的重復(fù)類型相同，但在杉木（Cunninghamia lanceolata）（文亞峰等，2015）、思茅松（P. kesiya）（趙能等，2017）和金鐵鎖（葉鵬等，2019a）等植物中以三核苷酸重復(fù)基元類型數(shù)量最多，其次是其他核苷酸重復(fù)基元類型，造成這種差異的原因可能與物種的不同有關(guān)。

SSR多態(tài)性是判斷其有效性的主要依據(jù)，而其主要影響因素是SSR長(zhǎng)度（李翠婷等，2014）。研究表明，當(dāng)SSR長(zhǎng)度大于20 bp時(shí)多態(tài)性較高，長(zhǎng)度為12～20 bp時(shí)呈多態(tài)性中等水平，長(zhǎng)度小于12 bp時(shí)多態(tài)性很低（Temnykh et al.，2001;楊華等，2011）。云南火焰蘭轉(zhuǎn)錄組中，有4422條SSR基序大小主要集中在12～30 bp，占SSR總數(shù)的95.37%，其中，基序大小12～20 bp的SSR有3593條，占SSR總數(shù)的77.49%，呈中等多態(tài)性;長(zhǎng)度大于20 bp的基序有1275條，占SSR總數(shù)的27.50%，呈較高多態(tài)性。Gao等（2003）研究發(fā)現(xiàn)，重復(fù)次數(shù)和SSR分子標(biāo)記多態(tài)性間呈正相關(guān)。Dreisigacker等（2004）研究表明，低級(jí)基元SSR多態(tài)性普遍高于高級(jí)基元SSR。本研究表明，云南火焰蘭轉(zhuǎn)錄組SSR中低級(jí)基元二核苷酸和三核苷酸長(zhǎng)度大于20 bp的基序有1022條，占SSR總數(shù)的22.04%，推測(cè)云南火焰蘭SSR具有較高的多態(tài)性。因此，該推論可作為潛在的SSR引物序列信息，極大地提高引物設(shè)計(jì)和開(kāi)發(fā)的目的性，為后續(xù)SSR引物的大量開(kāi)發(fā)提供基礎(chǔ)數(shù)據(jù)。

4 結(jié)論

云南火焰蘭轉(zhuǎn)錄組SSR中低級(jí)基元類型較豐富，開(kāi)發(fā)出多態(tài)性高的SSR引物潛力較高，可用于云南火焰蘭乃至火焰蘭屬植物的遺傳多樣性及種質(zhì)資源評(píng)價(jià)、遺傳圖譜構(gòu)建及遺傳育種等研究。

參考文獻(xiàn)：

蔡年輝，許玉蘭，徐楊，鄧麗麗，周麗，王大偉，何承忠，段安安. 2015. 云南松轉(zhuǎn)錄組SSR的分布及其序列特征[J]. 云南大學(xué)學(xué)報(bào)（自然科學(xué)版），37（5）：770-778. [Cai N H，Xu Y L，Xu Y，Deng L L，Zhou L，Wang D W，He C Z，Duan A A. 2015. The distribution and characteristics of SSR in trsnscriptome of Pinus yunnanensis[J]. Journal of Yunnan University（Natural Science），37（5）：770-778.]

曹雯靜，李虹俠，唐源江，吳坤林. 2014. 4個(gè)火焰蘭屬內(nèi)及屬間雜交新組合F1代核型分析[J]. 亞熱帶植物科學(xué)，43（1）：19-23. [Cao W J，Li H X，Tang Y J，Wu K L. 2014. Karyotype analysis of four intragenus and intergenus hybridized combinations F1 generation of Renanthera[J]. Subtropical Plant Science，43（1）：19-23.]

陳紅波. 2019. 基于白及轉(zhuǎn)錄組數(shù)據(jù)的EST-SSR分子標(biāo)記開(kāi)發(fā)、驗(yàn)證及應(yīng)用[D]. 遵義：遵義醫(yī)科大學(xué). [Chen H B. 2019. The development，validation and utilization of EST-SSR markers based on the transcriptome data in Bletilla striata[D]. Zunyi：Zunyi Medical University.]

陳和明，呂復(fù)兵，肖文芳，李佐，蔣明殿. 2019. 蝴蝶蘭與火焰蘭遠(yuǎn)緣雜交育種初探[J]. 中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，24（8）：60-71. [Chen H M，Lü F B，Xiao W F，Li Z，Jiang M D. 2019. Preliminary study on the intergeneric hybridization breeding between Phalaenopsis and Renanthera[J]. Journal of China Agricultural University，24（8）：60-71.]

陳起馨. 2016. 蘭花分子圖譜的構(gòu)建[D]. 廣州： 華南農(nóng)業(yè)大學(xué). [Chen Q X. 2016. A linkage map construction of Cymbidium using SSR marker[D]. Guangzhou：South China Agricultural University.]

陳曉明，李魁鵬，陳博雯，劉青華，周志春. 2018. 馬尾松轉(zhuǎn)錄組SSR序列特征分析及其分子標(biāo)記開(kāi)發(fā)[J]. 分子植物育種，16（19）：6407-6414. [Chen X M，Li K P，Chen B W，Liu Q H，Zhou Z C. 2018. Characterization analysis of SSR sequences in the transcriptome of Pinus massoniana and its molecular marker development[J]. Molecular Plant Breeding，16（19）：6407-6414.]

陳茵，李翠婷，姜倪皓，楊生超，陳軍文，楊建文，張廣輝. 2014. 燈盞花轉(zhuǎn)錄組中SSR位點(diǎn)信息分析及其多態(tài)性研究[J]. 中國(guó)中藥雜志，39（7）：1220-1224. [Chen Y，Li C T，Jiang N H，Yang S C，Chen J W，Yang J W，Zhang G H. 2014. SSR information in Erigeron breviscapus transcriptome and polymorphism analysis[J]. China Journal of Chinese Materia Medica，39（7）：1220-1224.]

程小毛，黃曉霞. 2011. SSR標(biāo)記開(kāi)發(fā)及其在植物中的應(yīng)用[J]. 中國(guó)農(nóng)學(xué)通報(bào)，27（5）：304-307. [Cheng X M，Huang X X. 2011. Development and application of SSR markers in plants[J]. Chinese Agricultural Science Bulletin，27（5）：304-307.]

郭鶯，孟紅巖，林文珍，林志楷，劉黎卿. 2018. 牛樟EST-SSR標(biāo)記的開(kāi)發(fā)及遺傳多態(tài)性分析[J]. 熱帶作物學(xué)報(bào)，39（8）：1561-1569. [Guo Y，Meng H Y，Lin W Z，Lin Z K，Liu L Q. 2018. Developing EST-SSR markers in Cinnamomum kanehirae and analyzing the genetic polymorphism[J]. Chinese Journal of Tropical Crops，39（8）：1561-1569.]

賀丹，吳芳芳，張佼蕊，謝棟博，李小康，劉藝平，栗燕，何松林. 2019. 牡丹轉(zhuǎn)錄組SSR信息分析及其分子標(biāo)記開(kāi)發(fā)[J]. 江蘇農(nóng)業(yè)學(xué)報(bào)，35（6）：1428-1433. [He D，Wu F F，Zhang J R，Xie D B，Li X K，Liu Y P，Li Y，He S L. 2019. Analysis of SSR information in transcriptome and development of molecular markers in Paeonia suffruticosa[J]. Jiangsu Journal of Agricultural Science，35（6）：1428-1433.]

揭應(yīng)碧，盧迎春，宋婉玲，韋坤華，張廣輝，楊生超. 2018. 三七基因組SSR位點(diǎn)分析和多態(tài)性引物開(kāi)發(fā)[J]. 熱帶作物學(xué)報(bào)，39（7）：1351-1358. [Jie Y B，Lu Y C，Song W L，Wei K H，Zhang G H，Yang S C. 2018. Analysis of SSR loci information in Panax notoginseng genomes and polymorphic primers development[J]. Chinese Journal of Tropical Crops，39（7）：1351-1358.]

李翠婷，張廣輝，馬春花，孟珍貴，陳軍文，楊生超. 2014. 野三七轉(zhuǎn)錄組中SSR位點(diǎn)信息分析及其多態(tài)性研究[J]. 中草藥，45（10）：1468-1472. [Li C T，Zhang G H，Ma C H，Meng Z G，Chen J W，Yang S C. 2014. Analysis on SSR loci information in transcriptome of Panax vienamensis var. fuscidiscus and its polymorphism[J]. Chinese Traditional and Herbal Drugs，45（10）：1468-1472.]

李冬梅，呂復(fù)兵，朱根發(fā)，洪彥彬，陳小平，梁炫強(qiáng)，操軍喜，孫映波. 2011. 蝴蝶蘭EST資源的SSR信息分析[J]. 廣東農(nóng)業(yè)科學(xué)，38（3）：117-120. [Li D M，Lü F B，Zhu G F，Hong Y B，Chen X P，Liang X Q，Cao J X，Sun Y B. 2011. Analysis of SSR information in EST resource of Phalaenopsis spp.[J]. Guangdong Agricultural Science，38（3）：117-120.]

李娟，石明，呂亞媚，楊錦超，杜凡. 2018. 極危植物云南火焰蘭的發(fā)現(xiàn)及種群現(xiàn)狀[J]. 林業(yè)科學(xué)研究，31（3）：9-14. [Li J，Shi M，Lü Y M，Yang J C，Du F. 2018. Discovery of an extremely endangered species Renanthera imschootiana and its population situation[J]. Forest Research，31（3）：9-14.]

李俊仁，陳秀珍，湯小婷，何瑞，詹若挺. 2018. 穿心蓮轉(zhuǎn)錄組SSR位點(diǎn)信息分析[J]. 中國(guó)中藥雜志，43（12）：2503-2508. [Li J R，Chen X Z，Tang X T，He R，Zhan R T. 2018. SSR loci information analysis in transcriptome of Andrographis paniculata[J]. China Journal of Chinese Materia Medica，43（12）：2503-2508.]

李清，李標(biāo)，郭順星. 2017. 金釵石斛轉(zhuǎn)錄組SSR位點(diǎn)信息分析[J]. 中國(guó)中藥雜志，42（1）：63-69. [Li Q，Li B，Guo S X. 2017. SSR information in transcriptome of Dendro-bium nobile[J]. China Journal of Chinese Materia Medica，42（1）：63-69.]

黎瑞源，邢輝，申鐵. 2017. 四球茶轉(zhuǎn)錄組SSR位點(diǎn)信息分析[J]. 安徽農(nóng)業(yè)大學(xué)學(xué)報(bào)，44（4）：558-562. [Li R Y，Xing H，Shen T. 2017. Analysis of SSR information in Camellia tetracocca Zhang transcriptome[J]. Journal of Anhui Agricultural University，44（4）：558-562.]

林丹妮，曾宋君，陳之林，段俊. 2006. 云南火焰蘭的離體快速繁殖[J]. 植物生理學(xué)通訊，42（1）：77. [Lin D N，Zeng S J，Chen Z L，Duan J. 2006. In vitro propagation of Renanthera imschootiana Rolfe[J]. Plant Physiology Communication，42（1）：77.]

林丹妮，陳之林，段俊，吳坤林，曾慶君. 2008. 云南火焰蘭種子無(wú)菌萌發(fā)和試管育苗[J]. 熱帶亞熱帶植物學(xué)報(bào)，16（1）：83-88. [Lin D N，Chen Z L，Duan J，Wu K L，Zeng Q J. 2008. Seed germination in vitro of Renanthera inschootiana Rolfe[J]. Journal of Tropical and Subtropical Botany，16（1）：83-88.]

林艷，何紫迪，毛積鵬，蔣開(kāi)彬，王金榜，黃少偉. 2018. 木荷基因組SSR位點(diǎn)開(kāi)發(fā)及初步分析[J]. 熱帶作物學(xué)報(bào)，39（9）：1766-1771. [Lin Y，He Z D，Mao J P，Jiang K B，Wang J B，Huang S W. 2018. Development and preliminary analysis of SSR locus in Schima superba genome[J]. Chinese Journal of Tropical Crops，39（9）：1766-1771.]

邱道壽，鄭希龍，蔡時(shí)可，鄭錦榮，羅煥明，張蕾，鄧瑞云，李武，劉曉津. 2013. 石斛SSR標(biāo)記的開(kāi)發(fā)及可轉(zhuǎn)移性分析[J]. 植物科學(xué)學(xué)報(bào)，31（5）：500-509. [Qiu D S，Zheng X L，Cai S K，Zheng J R，Luo H M，Zhang L，Deng R Y，Li W，Liu X J. 2013. Development and transfer analysis of SSR in Dendrobium[J]. Plant Science Journal，31（5）：500-509.]

申響保，閆文德，梁小翠，孫志強(qiáng). 2016. 杜仲轉(zhuǎn)錄組微衛(wèi)星特征分析[J]. 經(jīng)濟(jì)林研究，34（4）：47-51. [Shen X B，Yan W D，Liang X C，Sun Z Q. 2016. Characteristic analysis of microsatellites in Eucommia ulmoides transcriptome[J]. Nonwood Forest Research，34（4）：47-51.]

唐玉娟，黃國(guó)弟，羅世杏，周俊岸，莫永龍，李日旺，趙英，張宇，宋恩亮，寧琳. 2018. 芒果2個(gè)不同花分化時(shí)期轉(zhuǎn)錄組分析[J]. 南方農(nóng)業(yè)學(xué)報(bào)，49（7）：1257-1264. [Tang Y J，Huang G D，Luo S X，Zhou J A，Mo Y L，Li R W，Zhao Y，Zhang Y，Song E L，Ning L. 2018. Transcriptome of Mangifera indica L. in two different flower bud differentiation stages[J]. Journal of Southern Agriculture，49（7）：1257-1264.]

王大瑋，保云瑩，唐紅燕，段安安，蔡年輝，許玉蘭，周軍，李思廣. 2018. 思茅松SSR-PCR反應(yīng)體系優(yōu)化研究[J]. 西南林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），38（5）：34-37. [Wang D W，Bao Y Y，Tang H Y，Duan A A，Cai N H，Xu Y L，Zhou J，Li S G. 2018. Optimization of SSR-PCR reaction system in Pinus kesiya var. langbianensis[J]. Journal of Southwest Forestry University（Natural Science），38（5）：34-37.]

文亞峰，韓文軍，周宏，徐剛標(biāo). 2015. 杉木轉(zhuǎn)錄組SSR挖掘及EST-SSR標(biāo)記規(guī)?；_(kāi)發(fā)[J]. 林業(yè)科學(xué)，51（11）：40-49. [Wen Y F，Han W J，Zhou H，Xu G B. 2015. SSR mining and development of EST-SSR markers for Cunninghamia lanceolata based on transcriptome sequences[J]. Scientia Silvae Sinicae，51（11）：40-49.]

吳征鎰，洪德元. 1999. 中國(guó)植物志（第19卷）[M]. 北京：科學(xué)出版社. [Wu Z Y，Hong D Y. 1999. Flora of China（Vol. 19）[M]. Beijing：Science Press.]

肖冬長(zhǎng)，張智俊，管雨. 2012. 鐵皮石斛微衛(wèi)星SSR設(shè)計(jì)與應(yīng)用[J]. 生物技術(shù)通報(bào)，（7）：88-92. [Xiao D C，Zhang Z J，Guan Y. 2012. Microsatellite markers primer designing and screening from Dendrobium officinale[J]. Biotechno-logy Bulletin，（7）：88-92.]

楊華，陳琪，韋朝領(lǐng)，史成穎，方從兵，宛曉春. 2011. 茶樹(shù)轉(zhuǎn)錄組中SSR位點(diǎn)的信息分析[J]. 安徽農(nóng)業(yè)大學(xué)學(xué)報(bào)，38（6）：882-886. [Yang H，Chen Q，Wei C L，Shi C Y，F(xiàn)ang C B，Wan X C. 2011. Analysis on SSR information in Camellia sinensis transcriptome[J]. Journal of Anhui Agricultural University，38（6）：882-886.]

楊帥，楊喜翠，何方，李魯華，李振華，任明見(jiàn)，徐如宏. 2020. 貴紫麥1號(hào)轉(zhuǎn)錄組EST-SSR標(biāo)記發(fā)掘與應(yīng)用[J]. 河南農(nóng)業(yè)科學(xué)，49（5）：15-30. [Yang S，Yang X C， He F，Li L H， Li Z H， Ren M J，Xu R H. 2020. Exploration and app-lication of EST-SSR markers in Guizimai No.1 transcriptome[J]. Journal of Henan Agricultural Sciences，49（5）：15-30.]

姚婧. 2015. 火焰蘭種質(zhì)資源鑒定評(píng)價(jià)及組織培養(yǎng)技術(shù)的研究[D]. 海口：海南大學(xué). [Yao J. 2015. Identification and evaluation of germplasm resources and study on in vitro culture technique of Renanthera spp.[D]. Haikou：Hai-nan University.]

葉鵬，李顯煌，陳杰，瞿紹宏，董章宏，雷瀚，劉成，辛培堯，楊生超. 2019a. 金鐵鎖轉(zhuǎn)錄組SSR的分布及其序列特征[J]. 云南農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)），34（5）：808-814. [Ye P，Li X H，Chen J，Qu S H，Dong Z H，Lei H，Liu C，Xin P Y，Yang S C. 2019a. The distribution and characte-ristics of SSR in transcriptome of Psammosilene tunicoides[J]. Journal of Yunnan Agricultural University（Natural Science），34（5）：808-814.]

葉鵬，李顯煌，唐軍榮，李斌，張貴良，劉成，雷瀚，辛培堯. 2019b. 云南金花茶轉(zhuǎn)錄組SSR的分布及其序列特征[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào)，39（9）：86-91. [Ye P，Li X H，Tang J R，Li B，Zhang G L，Liu C，Lei H，Xin P Y. 2019b. Distribution and characteristics of SSR in transcriptome of Camellia fascicularis[J]. Journal of Central South University of Forestry & Technology，39（9）：86-91.]

張金旺，李代麗，尚峰男，阿茹娜，鐵龍，王君，姜鵬. 2018. 15個(gè)北方楊樹(shù)栽培品種倍性檢測(cè)與SSR遺傳分析[J]. 西南林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），38（5）：27-33. [Zhang J W，Li D L，Shang F N，A R N，Tie L，Wang J，Jiang P. 2018. Ploidy detection and genetic analysis of SSR mar-kers of 15 Populus cultivars in North China[J]. Journal of Southwest Forestry University（Natural Science），38（5）：27-33.]

張亞楠，楊柏云，熊冬金，羅火林. 2017. 寒蘭轉(zhuǎn)錄組SSR信息分析及其分子標(biāo)記開(kāi)發(fā)[J]. 南昌大學(xué)學(xué)報(bào)（理科版），41（3）：249-254. [Zhang Y N，Yang B Y，Xiong D J，Luo H L. 2017. Analysis on SSR information based on transcriptome and development of molecular markers in Cymbi-diun kanran[J]. Journal of Nanchang University（Natural Science），41（3）：249-254.]

張震，許彥明，陳永忠，李志鋼，王湘南，陳隆升，彭邵鋒，馬力，王瑞，李美群，唐煒，彭映赫. 2018. 油茶轉(zhuǎn)錄組測(cè)序與SSR特征分析[J]. 西南林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)），38（6）：63-68. [Zhang Z，Xu Y M，Chen Y Z，Li Z G，Wang X N，Chen L S，Peng S F，Ma L，Wang R，Li M Q，Tang W，Peng Y H. 2018. Transcriptome sequencing and analysis of SSR characteristics of Camellia oleifera[J]. Journal of Southwest Forestry University（Natural Science），38（6）：63-68.]

趙能，原曉龍，繆?？?，王毅，陳偉，李江，吳濤，王娟. 2017. 思茅松轉(zhuǎn)錄組SSR分析及標(biāo)記開(kāi)發(fā)[J]. 生物技術(shù)通報(bào)，33（5）：71-77. [Zhao N，Yuan X L，Miao F J，Wang Y，Chen W，Li J，Wu T，Wang J. 2017. Development of SSR molecular markers based on transcriptome data of Pinus kesiya var. langbianensis[J]. Biotechnology Bulletin，33（5）：71-77.]

中國(guó)科學(xué)院昆明植物研究所. 2003. 云南植物志（第14卷）[M]. 北京：科學(xué)出版社. [Kunming Institute of Botany，Chinese Academy of Sciences. 2003. Flora of Yunnan（Vol.14）[M]. Beijing： Science Press.]

周天華，丁家璽，田偉，王家. 2017. 天麻基因組微衛(wèi)星特征分析與分子標(biāo)記開(kāi)發(fā)[J]. 西北植物學(xué)報(bào)，37（9）：1728-1735. [Zhou T H，Ding J X，Tian W，Wang J. 2017. Genomic microsstellite characteristic analysis and molecular marker development for Gastrodia elata BI.[J]. Acta Botanica Boreali-Occidentalia Sinica，37（9）：1728-1735.]

Dreisigacker S，Zhang P，Warburton M L，Ginkel M V，Hoi-sington D，Bohn M，Melchinger A E. 2004. SSR and pedi-gree analyses of genetic diversity among CIMMYT wheat lines targeted to different megaenvironments[J]. Crop Science，44（2）：381-388.

Gao L F，Tang J F，Li H W，Jia J Z. 2003. Analysis of microsatellites in major crops assessed by computational and experimental approaches[J]. Molecular Breeding，12（3）：245-261.

Temnykh S，Declerck G，Lukashova A，Lipovich L，Cartinhour S，McCouch S. 2001. Computational and experimental analysis of microsatellites in rice（Oryza sativa L.）：frequency，length variation，transposon associations，and genetic marker potential[J]. Genome Research，11（8）：1441-1452.

（責(zé)任編輯 陳 燕）