基于高通量測(cè)序的思茅松微衛(wèi)星位點(diǎn)的特征分析

鄧麗麗，李德龍 ，蔡年輝 ，周 軍 ，孫 琪 ，唐紅燕 ，王大瑋 ，許玉蘭

（1. 西南林業(yè)大學(xué) 云南省高校林木遺傳改良與繁育重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650224；2. 普洱市林業(yè)科學(xué)研究所，云南 普洱 665099）

基于高通量測(cè)序的思茅松微衛(wèi)星位點(diǎn)的特征分析

鄧麗麗1，李德龍1，蔡年輝1，周 軍1，孫 琪1，唐紅燕2，王大瑋1，許玉蘭1

（1. 西南林業(yè)大學(xué) 云南省高校林木遺傳改良與繁育重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650224；2. 普洱市林業(yè)科學(xué)研究所，云南 普洱 665099）

以思茅松針葉為材料，采用Illumina Hiseq2000平臺(tái)測(cè)序，得到121 882條無冗余的序列，對(duì)這些序列進(jìn)行SSR位點(diǎn)搜索，共獲得3534個(gè)SSR位點(diǎn)，出現(xiàn)頻率為2.9%，分布平均距離25kb。所搜索的SSR以單核苷酸重復(fù)類型最多，其次為三核苷酸和二核苷酸，而四、五、六核苷酸重復(fù)類型較少（＜1%）。單核苷酸重復(fù)類型中以A/T基元較豐富（43.58%）；二核苷酸中出現(xiàn)頻率最高的是AT/TA基元（14.04%），AG/CT次之，AC/GT和CG/CG較低；三核苷酸重復(fù)類型中AAG/CTT最多，AGC/CTG和AGG/CCT次之，以CCG/CGG、ACG/CGT和ACT/AGT較少（＜1%）；四、五、六核苷酸類型中各重復(fù)基元相差不大，均較少。SSR數(shù)量隨對(duì)應(yīng)的重復(fù)類型重復(fù)次數(shù)的增加而降低，也隨重復(fù)區(qū)段堿基長(zhǎng)度的增加而降低。

思茅松；高通量測(cè)序；微衛(wèi)星

思茅松Pinus kesiyavar.langbianensis(A. Chev.)Gaussen是卡西亞松P.kesiyaKoyle ex Gordon的一個(gè)地理變種[1]，在我國(guó)分布于云南省亞熱帶、準(zhǔn)熱帶地區(qū)，是主要用材樹種[2]。自20世紀(jì)80年代以來對(duì)思茅松的遺傳改良做了許多工作，取得了階段性的成果和成熟的經(jīng)驗(yàn)[3]。目前，有關(guān)思茅松苗木培育[4]、森林資源管理[5]、遺傳變異[6]等方面開展多方位的研究，而從分子水平開展研究不多，集中在實(shí)驗(yàn)方法的建立[7-9]，利用分子標(biāo)記對(duì)種質(zhì)資源的評(píng)價(jià)研究更為鮮見，僅有少數(shù)利用同工酶及RAPD開展天然種群及種子園遺傳多樣性的研究[10-11]，迄今為止對(duì)思茅松的遺傳多樣性和遺傳結(jié)構(gòu)的了解仍十分欠缺，缺乏系統(tǒng)的遺傳學(xué)、群體遺傳學(xué)等基礎(chǔ)研究工作，僅局限于部分地區(qū)的有限群體或種子園，而且所選用的檢測(cè)手段以表型或顯性標(biāo)記為主，信息含量較少。因此，從DNA水平上深入研究思茅松群體遺傳結(jié)構(gòu)與分化就顯得更為重要，可為該種質(zhì)資源的保護(hù)與利用提供科學(xué)指導(dǎo)。

微衛(wèi)星是一種共顯性的分子標(biāo)記，因其多態(tài)性豐富、重復(fù)性高、覆蓋面廣、易于檢測(cè)、穩(wěn)定性好等優(yōu)點(diǎn)，廣泛應(yīng)用于遺傳多樣性分析、親緣關(guān)系分析、品種鑒定和輔助選擇育種等方面[12-19]。盡管微衛(wèi)星在許多在很多領(lǐng)域中得到廣泛應(yīng)用，但存在微衛(wèi)星的分離問題[20]。傳統(tǒng)的方法是先構(gòu)建文庫，再采用探針對(duì)所建立的文庫進(jìn)行篩選[21]，該方法耗時(shí)、耗力，且效率低，后來研究者相繼提出一些改進(jìn)的方法，如選擇性雜交、磁珠富集等，這些方法很快就取代傳統(tǒng)的耗時(shí)耗力的建庫掃描的方法[20]。近年來隨著高通量測(cè)序技術(shù)的發(fā)展，EST序列開發(fā)SSR分子標(biāo)記也越來越多，是一個(gè)比較有效的途徑[22-24]，與基因組SSR相比，EST-SSR具有通用性好、開發(fā)簡(jiǎn)便等特點(diǎn)，且來源于基因編碼序列，獲得基因表達(dá)信息，可開展功能基因的直接鑒定[25]，這方面的研究也越來越受到人們的關(guān)注。但是，目前關(guān)于思茅松基于EST序列的SSR研究未見報(bào)道。鑒于此，本研究基于Illumina Hiseq 2000平臺(tái)對(duì)思茅松轉(zhuǎn)錄組測(cè)序，搜索并分析SSR位點(diǎn)信息與分布特征，為思茅松微衛(wèi)星分子標(biāo)記的開發(fā)提供信息，以豐富思茅松分子標(biāo)記信息，為分子水平上深入研究思茅松遺傳變異奠定基礎(chǔ)，同時(shí)也可為近緣種特種SSR標(biāo)記開發(fā)及其遺傳分析提供便利。

1 材料與方法

1.1 研究材料

采集思茅松植株當(dāng)年生幼嫩針葉，干冰保存，核酸提取、測(cè)序（Illumina Hiseq2000）與組裝（Trinity軟件）由北京華諾時(shí)代科技有限公司完成。

1.2 SSR位點(diǎn)的搜索與統(tǒng)計(jì)

使 用 Micro Satellite（MISA,http://pgrc.ipkgatersleben. de/misa/）搜索轉(zhuǎn)錄組序列中潛在的1～6 bp的SSR，搜索標(biāo)準(zhǔn)為：?jiǎn)魏塑账嶂辽僦貜?fù)10次、二至四核苷酸至少重復(fù)5次、五核苷酸和六核苷酸至少重復(fù)4次。采用Excel軟件統(tǒng)計(jì)SSR的數(shù)量、出現(xiàn)頻率、分布距離與密度、重復(fù)類型、基元組成、重復(fù)區(qū)段長(zhǎng)度變化等，分析思茅松轉(zhuǎn)錄組SSR的分布特征。

2 結(jié)果與分析

2.1 思茅松轉(zhuǎn)錄組微衛(wèi)星序列的分布豐度與距離

組裝獲得的思茅松轉(zhuǎn)錄組共121 882條無冗余的轉(zhuǎn)錄本序列，按標(biāo)準(zhǔn)搜索1～6 bp的SSR，共搜索到3 534個(gè)SSR位點(diǎn)，其中復(fù)合型SSR的數(shù)量為96個(gè)，分布于3205條序列上，發(fā)生頻率（含SSR序列數(shù)/總序列數(shù)）為2.63%，其中2 924條序列含有單個(gè)SSR，281條序列含有2個(gè)或2個(gè)以上的SSR，總的SSR序列出現(xiàn)頻率（檢出的SSR個(gè)數(shù)/總序列數(shù)）為2.90%，平均距離為25.12 kb，即每25.12 kb出現(xiàn)1個(gè)SSR位點(diǎn)（表1）。

表1 思茅松轉(zhuǎn)錄組SSR位點(diǎn)的分布Table 1 The distribution of SSR repeat types in Piuns kesiya var. langbianensis transcriptome

在所搜索到的SSR中單核苷酸重復(fù)類型所占比例較高（44.74%），三核苷酸重復(fù)類型占總數(shù)的29.68%，二核苷酸重復(fù)類型占總數(shù)的23.80%，而四、五、六核苷酸重復(fù)單元所占比例均較低，分別為0.03%、0.01%和0.01%；若不考慮單核苷酸的情況下，二、三核苷酸占總數(shù)的96.77%。相應(yīng)地各重復(fù)單元類型的出現(xiàn)頻率、分布密度與SSR含量的變化規(guī)律相同，表現(xiàn)為：?jiǎn)魏塑账幔救塑账幔径塑账幔舅暮塑账幔疚搴塑账幔玖塑账?；不同重?fù)類型核苷酸的平均距離各不相同，其順序依次為：?jiǎn)魏塑账幔既塑账幔级塑账幔妓暮塑账幔嘉搴塑账幔剂塑账?，其中單核苷酸與六核苷酸重復(fù)類型的分布距離差異達(dá)113倍，即該轉(zhuǎn)錄組序列中每出現(xiàn)113個(gè)單核苷酸重復(fù)類型才出現(xiàn)一個(gè)六核苷酸重復(fù)類型的SSR。由表也可以看出，不同重復(fù)類型SSR位點(diǎn)的出現(xiàn)頻率或數(shù)量越高，其平均距離越小，分布密度越大。

2.2 思茅松轉(zhuǎn)錄組微衛(wèi)星序列的重復(fù)基元類型和比例

思茅松轉(zhuǎn)錄組微衛(wèi)星序列各重復(fù)基元統(tǒng)計(jì)如表2所示。

不同核苷酸重復(fù)類型的重復(fù)基元所占的比例差異較大，單核苷酸重復(fù)類型中以A/T基元居多，占43.58%，而C/G較少，僅占總數(shù)的1.16%；二核苷酸重復(fù)類型中各基元所占比例依次為：AT/AT（14.04%）＞AG/CT（8.09%）＞AC/GT（1.61%）＞CG/CG（0.06%）；三核苷酸重復(fù)類型中AAG/CTT較多（6.28%），其次是AGC/CTG（5.69%）和AGG/CCT（4.73%），而以ACT/AGT、ACG/CGT和CCG/CGG較低（＜1%）；四核苷酸重復(fù)類型中共搜索到20種重復(fù)基元，各重復(fù)基元的比例較低，最高的僅為0.14%；五、六核苷酸中出現(xiàn)的基元種類少、比例低。四、五、六核苷酸中各重復(fù)基元出現(xiàn)的頻率均比較低，因此表2中未列出具體的基元類型。

表2 思茅松轉(zhuǎn)錄組SSR各重復(fù)基元及其比例Table 2 The relative percentage of SSR repeat motif in Piuns kesiya var. langbianensis transcriptome

2.3 思茅松轉(zhuǎn)錄組微衛(wèi)星序列重復(fù)次數(shù)

重復(fù)基元的重復(fù)次數(shù)變異引起重復(fù)片段序列長(zhǎng)度的變異，思茅松轉(zhuǎn)錄組微衛(wèi)星序列各重復(fù)類型的SSR重復(fù)次數(shù)如圖1所示。

圖1 思茅松轉(zhuǎn)錄組各微衛(wèi)星重復(fù)類型重復(fù)次數(shù)分布頻率Fig.1 Frequency of SSR repeat types for the various number of repeats in Piuns kesiya var. langbianensis transcriptome

由圖1可知，各重復(fù)類型的重復(fù)次數(shù)差異較大，波動(dòng)于5～24次，總體表現(xiàn)為重復(fù)次數(shù)隨重復(fù)單元中核苷酸數(shù)量的增加而減少，其中單核苷酸重復(fù)10～23次，二核苷酸重復(fù)6～13次、其余核苷酸重復(fù)5～9次，不同核苷酸重復(fù)類型主導(dǎo)的重復(fù)次數(shù)也不相同。由圖也可以看出，SSR出現(xiàn)的頻率隨重復(fù)次數(shù)的增加而降低；與此同時(shí)，隨著SSR重復(fù)類型堿基數(shù)的增加，重復(fù)次數(shù)減少。總體來看（見圖2），SSR的重復(fù)次數(shù)以5～10次占多數(shù)，占77.31%，11～15次的占18.08%，而重復(fù)次數(shù)在15次以上的不足5%，表現(xiàn)為隨著重復(fù)次數(shù)的增加，SSR數(shù)量呈降低的趨勢(shì)。

2.4 思茅松轉(zhuǎn)錄組微衛(wèi)星序列長(zhǎng)度分布

思茅松轉(zhuǎn)錄組SSR重復(fù)片段的長(zhǎng)度波動(dòng)于10～76 bp，其中以10～24 bp居多（見圖3）。由圖3和表1可知，單核苷酸重復(fù)類型的長(zhǎng)度和二核苷酸重復(fù)類型的長(zhǎng)度均變化于10～24 bp，平均長(zhǎng)度分別為11.75 bp和14.53 bp；三核苷酸重復(fù)類型變化于15～27 bp，平均長(zhǎng)度為16.35 bp；四核苷酸重復(fù)類型長(zhǎng)度20～24 bp，另有一個(gè)重復(fù)較高的，重復(fù)19次（76 bp），平均長(zhǎng)度為22.24 bp；五核苷酸重復(fù)類型長(zhǎng)度25～40 bp，平均26.67 bp；六核苷酸30～66 bp，平均37.29 bp。綜合來看，各重復(fù)區(qū)段堿基片段長(zhǎng)度變化表現(xiàn)為隨著重復(fù)類型堿基數(shù)的增加，片段的平均長(zhǎng)度也隨之增加，即從單核苷酸至六核核苷酸，片段的平均長(zhǎng)度呈遞增的趨勢(shì)（表1）。從圖3也可以看出，除六核苷酸外，其余核苷酸的重復(fù)類型表現(xiàn)為SSR出現(xiàn)的頻率隨重復(fù)區(qū)段堿基片段長(zhǎng)度的增加而降低，微衛(wèi)星重復(fù)區(qū)段片段長(zhǎng)度與其對(duì)應(yīng)的SSR數(shù)量或比例成相反的變化趨勢(shì)，即重復(fù)區(qū)段堿基片段短的，其SSR數(shù)目較多。從全部核苷酸來看，90%以上的重復(fù)區(qū)段片段長(zhǎng)度為10～20 bp，而大于30 bp的SSR不足1%。

圖2 思茅松轉(zhuǎn)錄組微衛(wèi)星重復(fù)次數(shù)分布頻率Fig.2 Frequency of repeat number of SSR in Piuns kesiya var. langbianensis transcriptome

圖3 思茅松轉(zhuǎn)錄組各SSR重復(fù)類型片段長(zhǎng)度分布頻率Fig.3 Frequency of SSR repeat types for the various sequence length in Piuns kesiya var. langbianensis transcriptome

3 結(jié)論與討論

思茅松轉(zhuǎn)錄組中SSR的發(fā)生頻率為2.63%，平均分布距離為25.12 kb，39個(gè)SSR/Mb，與云南松Pinus yunnanensis轉(zhuǎn)錄組中SSR的分布及發(fā)生頻率較為接近[26]，而與楊樹Populus（14.83%）、桉樹Eucalyptus（14.99%）的發(fā)生頻率相比較低[27-28]，這種差異大小可能與SSR搜索標(biāo)準(zhǔn)、數(shù)據(jù)庫大小以及物種有關(guān)[29]，也可能與含有微衛(wèi)星的基因的表達(dá)豐度有關(guān)[27]。思茅松轉(zhuǎn)錄組序列所檢測(cè)的SSR中，單核苷酸重復(fù)類型所占比例接近一半，另一半主要是由三核苷酸重復(fù)類型和二核苷酸重復(fù)類型組成，這在紅松Pinus koraiensis轉(zhuǎn)錄組的研究中也提到，單核苷酸重復(fù)類型占總的46.90%[30]，楊樹轉(zhuǎn)錄組檢測(cè)到的核苷酸也表現(xiàn)出相同的規(guī)律，即單核苷酸＞三核苷酸＞二核苷酸[31]，但在云南松轉(zhuǎn)錄組中以三核苷酸較多[26]，這種占優(yōu)勢(shì)的重復(fù)類型在不同植物中有所差異，許玉蘭等[32]對(duì)多物種的統(tǒng)計(jì)可知以二、三核苷酸較為集中。若不考慮單核苷酸的情況下，在2～6 bp的SSR中，思茅松也以2～3 bp重復(fù)類型為主，這在地中海松Pinus halepensis、馬尾松Pinus massoniana和火炬松Pinus taeda等研究中也有類似的報(bào)道[22,33-35]。不同重復(fù)類型的SSR多態(tài)水平存在一定的差異[36]，在研究報(bào)道中，擴(kuò)增出多態(tài)性位點(diǎn)的引物重復(fù)單元以二、三核苷酸重復(fù)為主[30,37-38]，李淑嫻等[27]在桉樹研究中也發(fā)現(xiàn)微衛(wèi)星長(zhǎng)度的變異速率與其長(zhǎng)度之間存在一定的關(guān)系，表現(xiàn)為低級(jí)基元SSR理論多態(tài)性較高[39]。因此，思茅松轉(zhuǎn)錄組含有較多低級(jí)基元單、二、三核苷酸，多態(tài)性SSR分子標(biāo)記的開發(fā)具有很大的潛力。

不同重復(fù)類型中各重復(fù)基元所在比例存在差異，在單核苷酸重復(fù)類型中，以A/T占絕多數(shù)，二核苷酸重復(fù)類型中，AT/TA占多數(shù)，表現(xiàn)出AT優(yōu)勢(shì)。這種現(xiàn)象在云南松[26]、紅松[30]、火炬松和馬尾松[33]中均有研究報(bào)道。AT含量較多可能與其能量有關(guān)，因?yàn)锳-T間為兩氫鍵，C-G間為三氫鍵，打破AT所需的能量低于CG[40]。但是，在三核苷酸重復(fù)類型中以AAG/CTT、AGC/CTG和AGG/CCT較多，分別為6.28%、5.69%和4.73%，與其它的松樹報(bào)道的比較相似，如紅松中以AGC/CTG、AGG/CCT和AAG/CTT三種較多，分別占總數(shù)的7.63%、6.66%和6.55%[30]；云南松中也以AGC/CTG、AAG/CTT和AGG/CCT[26]。在前期對(duì)多樹種的統(tǒng)計(jì)中也可看出以AAG、AGC和AGG較多，而在如檉柳Tamarixspp.、橡膠樹Hevea brasiliensis、白樺Betula platyphylla等中還出現(xiàn)少量其它優(yōu)勢(shì)基元，如AAT、AGA、TCA等[32]，表明不同植物SSR基元組成存在一定的差異，但一些優(yōu)勢(shì)基元如AAG、AGC和AGG在多種植物中均存在。因此，推測(cè)這些重復(fù)基元可能在EST序列中的存在較為普遍，可能是優(yōu)勢(shì)的蛋白和DNA家族[36]。CG重復(fù)單元可能與某些特定的功能相關(guān)[41]，在思茅松轉(zhuǎn)錄組中發(fā)現(xiàn)一定數(shù)量的CG，如二核苷酸CG，三核苷酸ACG、AGC等，這些重復(fù)類型是否能加其它的生命活動(dòng)需要進(jìn)一步的研究。但是CG的含量很少，有研究認(rèn)為CG重復(fù)類型少是因?yàn)榛蚪MDNA中的CpG甲基化，極易發(fā)生突變，較少的GC是維持熱力學(xué)穩(wěn)定必須的因素[42]，這些將為思茅松SSR引物的開發(fā)提供啟示。

微衛(wèi)星位點(diǎn)多態(tài)性主要原因是基元重復(fù)數(shù)和堿基數(shù)不同而形成的序列長(zhǎng)度多態(tài)性[43]，微衛(wèi)星核心序列重復(fù)次數(shù)與該位點(diǎn)等位基因數(shù)存在較強(qiáng)的正相關(guān)，一般重復(fù)次數(shù)多，其變異性越大，該位點(diǎn)的等位基因數(shù)越多[41]，思茅松轉(zhuǎn)錄組SSR豐度隨各對(duì)應(yīng)重復(fù)類型堿基數(shù)的增加而減少，這在其它的研究中也有類似的報(bào)道，如唐古特紅景天Rhodiola algida中核苷酸類型表現(xiàn)為長(zhǎng)度變化與其對(duì)應(yīng)的重復(fù)類型的堿基長(zhǎng)度成反比，同種重復(fù)類型中，微衛(wèi)星的長(zhǎng)度與對(duì)應(yīng)的SSR數(shù)成反比[41]。云南松中也出現(xiàn)隨著重復(fù)次數(shù)的增加、重復(fù)片段的增長(zhǎng)，SSR出現(xiàn)的數(shù)量降低[26]。由此表明，這種重復(fù)類型特征不是思茅松所特有。一般來說，重復(fù)區(qū)段堿基片段長(zhǎng)度隨著重復(fù)次數(shù)的不同而發(fā)生變化，而SSR長(zhǎng)度變異是產(chǎn)生SSR多態(tài)性的主要原因[39]，不同長(zhǎng)度微衛(wèi)星的分布及變化的分析有利于多態(tài)性高的SSR標(biāo)記的獲得[41]。從思茅松轉(zhuǎn)錄組SSR分析來看，重復(fù)次數(shù)為5～24次，其中單核苷酸因容易發(fā)生錯(cuò)配排除外[43]，其余的核苷酸重復(fù)類型重復(fù)次數(shù)也有5～13次，甚至有高達(dá)19次的。從片段長(zhǎng)度來看，90%的長(zhǎng)度≤20 bp，按照Temnykh等[39]的提出SSR長(zhǎng)度≤20 bp時(shí)多態(tài)性高，當(dāng)＞20 bp時(shí)，仍以低級(jí)基元類型為主（72.94%），這部分SSR可能也具有較多的多態(tài)性潛能。由此推測(cè)，思茅松轉(zhuǎn)錄組挖掘的3 534個(gè)SSR位點(diǎn)大部分具有多態(tài)性潛能，可用于分子標(biāo)記的開發(fā)。當(dāng)然，本研究?jī)H針對(duì)SSR位點(diǎn)分布特征進(jìn)行分析，對(duì)于引物的開發(fā)以及功能標(biāo)記的研究還需入進(jìn)一步深入。

4 小 結(jié)

思茅松轉(zhuǎn)錄組SSR出現(xiàn)頻率較高，類型豐富，且分布密度也較大。在搜索到的3 534個(gè)SSR位點(diǎn)中，以單核苷酸占多數(shù)，其次為三、二核苷酸重復(fù)類型，而四、五、六核苷酸重復(fù)類型的SSR數(shù)量較少。單核苷酸和二核苷酸重復(fù)類型中，以A/T、AT/TA基元出現(xiàn)的比例較高，CG含量相對(duì)較低。SSR位點(diǎn)出現(xiàn)的頻率越高，其平均距離越小；各重復(fù)單元片段長(zhǎng)度變隨著重復(fù)單元堿基數(shù)的增加而增加；SSR出現(xiàn)的頻率隨對(duì)應(yīng)重復(fù)類型重復(fù)區(qū)段堿基片段長(zhǎng)度的增加而降低；同樣地，隨著重復(fù)次數(shù)的增加，SSR數(shù)量呈降低的趨勢(shì)。結(jié)合其它研究，評(píng)價(jià)思茅松轉(zhuǎn)錄組所挖掘的3 534個(gè)SSR位點(diǎn)具有較高多態(tài)性潛能，可用于分子標(biāo)記的開發(fā)。

[1]中國(guó)科學(xué)院昆明植物研究所. 云南植物志: 第四卷[M]. 北京:科學(xué)出版社, 1986.

[2]中國(guó)樹木志編委會(huì). 中國(guó)主要樹種造林技術(shù)[M]. 北京: 中國(guó)林業(yè)出版社, 1981.

[3]許玉蘭, 段安安. 思茅松遺傳改良研究現(xiàn)狀及育種策略的探討[J]. 西南林學(xué)院學(xué)報(bào), 2004, 24(4): 68-72.

[4]徐玉梅, 楊德軍, 史富強(qiáng), 等. 思茅松人工林不同空間配置模式下不同密度生長(zhǎng)試驗(yàn)與分析[J]. 南京林業(yè)大學(xué)學(xué)報(bào): 自然科學(xué)版, 2015, 39(1): 162-166.

[5]歐光龍, 胥 輝, 王俊峰, 等. 思茅松天然林林分生物量混合效應(yīng)模型構(gòu)建[J]. 北京林業(yè)大學(xué)學(xué)報(bào), 2015, 37(3): 101-111.

[6]李帥鋒, 蘇建榮, 劉萬德, 等. 思茅松天然群體種實(shí)表型變異[J]. 植物生態(tài)學(xué)報(bào), 2013, 37(11): 998-1009.

[7]王曙光, 楊文武, 普曉蘭. 思茅松大配子體DNA的提取及ISSR反應(yīng)體系的建立[J]. 四川林業(yè)科技, 2012, 33(2): 90-95.

[8]姚鵬強(qiáng), 李 培, 王曙光, 等. 思茅松AFLP分子標(biāo)記體系的建立[J]. 西南林業(yè)大學(xué)學(xué)報(bào), 2013, 32(6): 30-33.

[9]魏 博, 汪元超, 王大瑋, 等. 思茅松SRAP-PCR反應(yīng)體系的建立與優(yōu)化[J]. 江蘇農(nóng)業(yè)科學(xué), 2014, 42(3): 27-30.

[10]陳少瑜, 趙文書, 王 炯. 思茅松天然種群及其種子園的遺傳多樣性[J]. 福建林業(yè)科技, 2002, 29(3): 1-5.

[11]朱云鳳, 陳少瑜, 郝佳波, 等. 思茅松種子園遺傳多樣性的RAPD分析[J]. 西部林業(yè)科學(xué), 2015, 43(2): 141-146.

[12]陳贏男, 張新葉, 戴曉港. 利用微衛(wèi)星標(biāo)記鑒別油茶品種[J].經(jīng)濟(jì)林研究, 2014, 32(4): 140-143.

[13]Sanchez M, Ingrouille M J, Cowan R S,et al. Spatial structure and genetic diversity of natural populations of the Caribbean pine,Pinus caribaeavar.bahamensis(Pinaceae), in the Bahaman archipelago[J]. Botanical Journal of the Linnean Society, 2014,17(3): 359-383.

[14]管俊嬌, 張建華, 張 惠, 等. 辣椒雜交種純度鑒定的SSR核心引物篩選[J]. 生物技術(shù)通報(bào), 2014, 30(9): 97-101.

[15]魏琦琦, 林 青, 賈寶光, 等. 棗轉(zhuǎn)錄組序列的微衛(wèi)星特征分析[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2015, 35(6): 93-97.

[16]黃少勇, 張智俊, 羅淑萍, 等. 山核桃EST-SSR引物的篩選及通用性分析[J]. 經(jīng)濟(jì)林研究, 2013, 31(3): 10-15.

[17]盧代華, 王 劍, 伏榮桃, 等. 水稻對(duì)稻曲病不同抗性品種遺傳多樣性SSR分析[J]. 西南農(nóng)業(yè)學(xué)報(bào), 2015, 28(2): 457-464.

[18]閆麗君, 王紅霞, 黃芳芳, 等. 伯樂樹SSR-PCR反應(yīng)體系的優(yōu)化[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2014, 34(12): 83-86.

[19]Mason A S. SSR Genotyping[J]. Methods in Molecular Biology,2015, 1245: 77-89.

[20]Zane L, Bargelloni L, Patarnello T. Strategies for microsatellite isolation: a review[J]. Molecular Ecology, 2002, 11(1):1-16.

[21]Rassmann K, Schlotterer C, Tautz D. Isolation of simple sequence loci for use in polymerase chain reaction-based DNA fingerprinting[J]. Electrophoresis, 1991, 12(2-3): 113-118.

[22]楊維澤, 許宗亮, 楊紹兵, 等. 三種植物EST-SSR引物在滇重樓上的通用性分析[J]. 西南農(nóng)業(yè)學(xué)報(bào), 2014, 27(4): 1686-1690.

[23]Pinosio S, González-Martínez S C, Bagnoli F,et al. First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species,Pinus halepensisMill[J]. Molecular Ecology Resources, 2014, 14(4):846-856.

[24]Liu J J, Hammett C. Development of novel polymorphic microsatellite markers by technology of next generation sequencing in western white pine[J]. Conservation Genetic Resource, 2014, 6(3): 647-648.

[25]Ellis J R, Burke J M. EST-SSRs as a resource for population genetic analyses[J]. Heredity, 2007, 99(2): 125-132.

[26]蔡年輝, 許玉蘭, 徐 楊, 等. 云南松轉(zhuǎn)錄組SSR的分布及其序列特征[J]. 云南大學(xué)學(xué)報(bào): 自然科學(xué)版,2015,37(5):770-778.

[27]李淑嫻, 張新葉, 王英亞, 等. 桉樹EST序列中微衛(wèi)星含量及相關(guān)特征[J]. 植物學(xué)報(bào), 2010, 45(3): 363-371.

[28]張新葉, 宋叢文, 張亞東, 等. 楊樹EST-SSR標(biāo)記的開發(fā)[J].林業(yè)科學(xué), 2009, 45(9): 53-59.

[29]Nicot N, Chiquet V, Gandon B,et al.Study of simple sequence repeat (SSR) markers from wheat expressed sequence tags(ESTs)[J]. Theoretical & Applied Genetics, 2004, 109(4): 800-805.

[30]張 振, 張含國(guó), 莫 遲, 等. 紅松轉(zhuǎn)錄組SSR分析及ESTSSR標(biāo)記開發(fā)[J]. 林業(yè)科學(xué), 2015, 51(8): 114-120.

[31]雷淑蕓, 張發(fā)起, Khan G, 等. 利用高通量測(cè)序分析青藏高原地區(qū)青楊的SSR和SNP特征[J]. 林業(yè)科學(xué)研究, 2015, 28(1):37-43.

[32]許玉蘭, 蔡年輝, 康向陽, 等. EST-SSR標(biāo)記的開發(fā)及其在木本植物中的分布特點(diǎn)[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2012, 28(4): 1-7.

[33]Bai TD, Xu L, Xu M,et al. Characterization of masson pinePinus massonianaLamb. microsatellite DNA by 454 genome shotgun sequencing[J]. Tree Genetics & Genomes, 2014, 10(2):429-437.

[34]Pinosio S, González-Martínez SC, Bagnoli F,et al.First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species,Pinus halepensisMill[J]. Molecular Ecology Resources, 2014, 14(4):846-856.

[35]Bérubé Y, Zhuang J, Rungis D,et al. Characterization of ESTSSRs in loblolly pine and spruce[J]. Tree Genetics & Genomes,2007, 3(3): 251-259.

[36]張得芳, 李淑嫻, 夏 濤. 薔薇科6個(gè)屬植物EST-SSR特征分析[J]. 植物研究, 2014, 34(6): 810-815.

[37]Pan F, Shihui N, Huwei Y,et al. Development and Characterization of 25 EST-SSR markers inPinus sylvestrisvar.mongolica(Pinaceae)[J]. Applications in Plant Sciences, 2014,2(1): 1300057.

[38]Xiang X, Zhang Z, Wang Z,et al.Transcriptome sequencing and development of EST-SSR markers inPinus dabeshanensis, an endangered conifer endemic to China[J]. Molecular Breeding,2015, 35(8): 1-10.

[39]Temnykh S, Declerck G, Lukashova A,et al. Computational and experimental analysis of microsatellites in rice (Oryza sativaL.):Frequency, length variation, transposon associations, and genetic marker potential[J]. Genome Research, 2001, 11(8): 1441-1452.

[40]童治軍, 肖炳光. 3種煙草基因組SSR位點(diǎn)信息分析和標(biāo)記開發(fā)[J]. 西北植物學(xué)報(bào), 2014, 34(8): 1549-1558.

[41]雷淑蕓, 高慶波, 付鵬程, 等. 基于Solexa高通量測(cè)序的唐古特紅景天Rhodiola algida微衛(wèi)星信息分析[J]. 植物研究,2014, 34(6): 829-834.

[42]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): 957-961.

[43]王 森, 張 震, 姜倪皓, 等. 半夏轉(zhuǎn)錄組中的SSR位點(diǎn)信息分析[J]. 中藥材, 2014, 37(9): 1566-1569.

Characteristic analysis of microsatellite inPinus kesiyavar.langbianensisby using high-throughput sequencing

DENG Li-li1, LI De-long1, CAI Nian-hui1, ZHOU Jun1, SUN Qi1, TANG Hong-yan2, WANG Da-wei1, XU Yu-lan1
(1. Key Laboratory for Forest Genetic and Tree Improvement & Propagation in Universities of Yunnan Province, Southwest Forestry University, Kunming 650224, Yunnan, China；2. Puer City Institute of Forestry Sciences, Puer 665099, Yunnan, China)

Needle sample ofPinus kesiyavar.langbianensiswas sequenced with high-throughput sequencing technology (Illumina Hiseq 2000). In total 121 882 transcripts ofPinus kesiyavar.langbianensis.were screened using MISA software. 3534 SSRs were identi fied, and the frequency of SSR was 2.9% and mean distribution density was 25 kb. The mono-nucleotide repeats were the highest,followed by tri-nucleotide and di-nucleotide repeat types. The tetra-nucleotide, petra-nucleotide and hexa-nucleotide were all less than 1%.Among the mononucleotide repeats, the A/T repeats motifs were the highest frequency (43.58%). AT/AT was the most frequent repeat motifs in di-nucleotide repeats (14.04%), followed by AG/CT. AC/GT and CG/CG were minimum. AAG/CTT repeats motifs was the highest in tri-nucleotide followed by AGC/CTG and AGG/CCT, while CCG/CGG, ACG/CGT and ACT/AGT were lowest (＜1%). The repeat motifs were very few in tetra-nucleotide, petra-nucleotide and hexa-nucleotide and had no obvious differences. The SSR number of decreased with the increased number of repeats and length of repeats.

Pinus kesiyavar.langbianensis; high-throughput sequencing; microsatellite

S791.259

A

1673-923X(2016)10-0072-06

10.14067/j.cnki.1673-923x.2016.10.013

2016-01-19

西南林業(yè)大學(xué)科技創(chuàng)新基金（15092）；國(guó)家自然科學(xué)基金項(xiàng)目（31500536）；西南林業(yè)大學(xué)科研啟動(dòng)基金項(xiàng)目（01102-111436）

鄧麗麗，碩士研究生

許玉蘭，副教授，博士；E-mail：xvyulan@163.com

鄧麗麗，李德龍，蔡年輝，等. 基于高通量測(cè)序的思茅松微衛(wèi)星位點(diǎn)的特征分析[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2016, 36(10):72-77, 93.

[本文編校：吳 彬]