王祝年 邱燕連 晏小霞 陳永鋒 葉才華

摘 要 對(duì)海南島72份益智種質(zhì)的遺傳多樣性進(jìn)行系統(tǒng)評(píng)價(jià)。利用相關(guān)分析、主成分分析探究各表型性狀的相互關(guān)系，結(jié)合表型性狀及ISSR標(biāo)記研究資源的遺傳多樣性。結(jié)果表明：不同益智種質(zhì)的各種性狀均存在變異，所調(diào)查的14個(gè)性狀可歸類為營(yíng)養(yǎng)器官、產(chǎn)量、果實(shí)形態(tài)等3個(gè)部分，各部分平均變異系數(shù)表現(xiàn)為CV產(chǎn)量（0.464）>CV營(yíng)養(yǎng)器官（0.143）>CV果實(shí)形態(tài)（0.121）。性狀相關(guān)分析中，91對(duì)性狀組合中有34對(duì)呈極顯著相關(guān)。主成分分析前7個(gè)主成分積累貢獻(xiàn)率達(dá)到89.29%，前3個(gè)貢獻(xiàn)率最大的主成分積累貢獻(xiàn)率達(dá)到64.86%，單果體積因子第一，產(chǎn)量因子第二，單果形態(tài)因子第三。經(jīng)基于表型性狀的聚類分析，72份樣本可分為9組。利用20條ISSR引物對(duì)68份益智種質(zhì)材料基因組DNA進(jìn)行研究，共擴(kuò)增出466條條帶，其中446條為多態(tài)性帶，占總條帶的95.71%；當(dāng)相似系數(shù)為0.762時(shí)，68份樣本被分為6組，其中42.65%的樣品被分到第Ⅲ組。依據(jù)ISSR標(biāo)記的聚類與表型性狀的表現(xiàn)并不完全一致，這可能是由于DNA結(jié)構(gòu)差異不一定在表型性狀上得到相應(yīng)表現(xiàn)，本研究所調(diào)查的表型性狀只是益智表型性狀的一部分，而ISSR標(biāo)記所分析的位點(diǎn)則覆蓋整個(gè)基因組。

關(guān)鍵詞 海南島；益智；表型性狀；ISSR；聚類分析

中圖分類號(hào) Q949.71+8.33 文獻(xiàn)標(biāo)識(shí)碼 A

益智（Alpinia oxyphylla Miquel）屬姜科山姜屬植物，為中國(guó)四大南藥之一，主產(chǎn)于中國(guó)海南，廣東、云南、廣西、福建亦有少量栽培。益智以干燥成熟果實(shí)入藥，已有悠久的藥用歷史，其性溫味辛，具有暖腎固精縮尿、溫脾止瀉攝唾的功能[1]。此外，益智還是極有開發(fā)價(jià)值的香料植物，已被廣泛應(yīng)用于保健食品和藥膳等，屬藥、香、食同源植物，具有廣泛的研究和應(yīng)用價(jià)值，已開發(fā)出系列產(chǎn)品。

近年來，學(xué)者們陸續(xù)對(duì)益智進(jìn)行了相關(guān)研究，取得了很大的進(jìn)展，研究范圍包括化學(xué)成分[2-4]、藥理[5-7]、栽培技術(shù)[8]、組織培養(yǎng)[9]、遺傳多樣性[10-11]、食品開發(fā)[12-14]等多方面的內(nèi)容，但關(guān)于益智遺傳育種方面的研究尚未見報(bào)道。海南島益智主要為農(nóng)戶自發(fā)引種栽培，各地間種質(zhì)資源交流頻繁，遺傳變異復(fù)雜，遺傳背景不明，尚未有專門的機(jī)構(gòu)對(duì)其進(jìn)行深入研究，基于此，本課題以海南島益智種質(zhì)資源為材料，結(jié)合形態(tài)標(biāo)記和ISSR分析，開展海南島益智種質(zhì)資源的遺傳多樣性研究，旨在為益智的遺傳改良和生產(chǎn)應(yīng)用以及海南島益智種質(zhì)資源的保護(hù)和合理利用提供科學(xué)依據(jù)。

表型變異是基因型與環(huán)境共同作用的結(jié)果，可在一定程度上反映生物遺傳變異的程度，是理解物種適應(yīng)機(jī)制的重要方法[15]。ISSR（inter-simple sequence repeats）簡(jiǎn)單序列重復(fù)區(qū)間擴(kuò)增多態(tài)性是由Zietkiewicz等[16]于1994年提出的一種以PCR（polymerase chain reactiong，聚合酶鏈?zhǔn)椒磻?yīng)）擴(kuò)增為基礎(chǔ)的分析標(biāo)記技術(shù)，被廣泛應(yīng)用于居群遺傳學(xué)、品種鑒定、物種分類以及物種的遺傳多樣性研究中[17-19]。

1 材料與方法

1.1 材料

本課題于2014年5～6月對(duì)海南島5個(gè)地區(qū)8個(gè)居群72份益智種質(zhì)材料性狀進(jìn)行觀測(cè)。各個(gè)采集地的地理位置、經(jīng)緯度、海拔高度、生境特點(diǎn)和樣本數(shù)見表1，由于各種原因，部分樣品被他人提前收走，未能獲得數(shù)據(jù)，導(dǎo)致各居群樣本數(shù)不一致。居群中除SY1（三亞鳳凰嶺居群）為純野生居群外，其余均為栽培居群。

1.2 方法

1.2.1 表型性狀的數(shù)據(jù)處理 每叢為一份種質(zhì)材料（或稱樣品），各表型性狀皆以測(cè)量該性狀的最長(zhǎng)（寬）處為準(zhǔn)，株高、莖桿直徑、果穗重、果序長(zhǎng)、座果率等性狀每叢測(cè)10個(gè)重復(fù)并取其平均值；每份樣品取30粒正常發(fā)育的單果，用游標(biāo)卡尺測(cè)量果長(zhǎng)、果寬。采集地點(diǎn)的地理因子用手持GPS定位儀記錄。

1.2.2 ISSR分析 從上述72份種質(zhì)材料中，取68份（分子標(biāo)記材料缺HM2-3、LM1-1、QS1-3、QS2-6）材料；采用改良的CTAB法[20]提取益智嫩葉總DNA，并利用紫外分光光度計(jì)（島津UV-2700）檢測(cè)DNA的純度和濃度，用10 g/L瓊脂糖凝膠電泳檢測(cè)DNA質(zhì)量，最后將每一份樣本DNA稀釋至50 ng/μL，于-20 ℃保存?zhèn)溆谩SSR引物為哥倫比亞大學(xué)（UBC）公布的100條引物序列，由美國(guó)生命技術(shù)公司廣州合成部（life technologies）合成。從確定的反應(yīng)體系中篩選出20條多態(tài)性好及擴(kuò)增條帶清晰的引物用于后續(xù)試驗(yàn)。ISSR反應(yīng)體系為：10×PCR buffer（Mg2+plus）3.0 μL，模板DNA（50 ng/μL） 1.0 μL，引物（100 μmol/L）1.0 μL，dNTPs（2.5 mmol/L）1.5 μL，Taq酶（5 U/μL）0.2 μL，加ddH2O補(bǔ)足至20 μL，總反應(yīng)體積為20 μL。PCR反應(yīng)程序?yàn)椋?4 ℃預(yù)變性5 min；94 ℃變性30 s，55 ℃退火45 s，72 ℃延伸2 min，共36個(gè)循環(huán)；72 ℃延伸10 min，4 ℃保存。PCR擴(kuò)增反應(yīng)在德國(guó)Biometra公司生產(chǎn)的T1 PCR擴(kuò)增儀上進(jìn)行，PCR反應(yīng)結(jié)束后取7 μL PCR產(chǎn)物與1 μL 6×loading Buffer混勻，點(diǎn)入15 g/L瓊脂糖凝膠（每100 mL加入5 μL GoldViewTM DNA染料）中，電泳緩沖液為0.5×TBE，電壓90 V，電泳約40～50 min；電泳結(jié)束后在凝膠成像系統(tǒng)中觀察并拍照。

1.3 數(shù)據(jù)處理

用Excel軟件對(duì)所獲得的表型性狀數(shù)據(jù)進(jìn)行統(tǒng)計(jì)，以變異系數(shù)CV（標(biāo)準(zhǔn)差/平均值）作為各形態(tài)性狀變異的測(cè)度；運(yùn)用SPSS19.0軟件對(duì)形態(tài)數(shù)據(jù)經(jīng)標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化后進(jìn)行系統(tǒng)聚類。以擴(kuò)增條帶在相對(duì)遷移位置的有無，記為“1”或“0”，用NTSYS-PC2.1對(duì)擴(kuò)增條帶數(shù)據(jù)進(jìn)行分析，通過類平均法（UPGMA）對(duì)68份種質(zhì)進(jìn)行遺傳相似性聚類。

2 結(jié)果與分析

2.1 表型性狀多樣性

樣品形態(tài)性狀指標(biāo)的分析結(jié)果（表2）表明，72份益智材料的14個(gè)性狀均存在一定程度的變異，變異系數(shù)為0.089～0.821。14個(gè)性狀可歸類為營(yíng)養(yǎng)器官、產(chǎn)量、果實(shí)形態(tài)等3個(gè)部分。營(yíng)養(yǎng)器官部分包括株高、莖桿直徑2個(gè)性狀。株高的變異幅度呈中等大小，從138～385 cm不等；莖桿直徑變異幅度較株高小。產(chǎn)量部分主要包含莖桿總數(shù)、結(jié)果桿數(shù)、單株產(chǎn)量、果穗重、果序長(zhǎng)、座果率、單果重等7個(gè)性狀。單株產(chǎn)量是所有性狀中變異程度最大的，變異系數(shù)為0.821，同時(shí)受到其余6個(gè)性狀的影響。果實(shí)形態(tài)包括果長(zhǎng)、果寬、單果大?。üL(zhǎng)×果寬）、單果形態(tài)（果長(zhǎng)/果寬）及果頂端形態(tài)等5個(gè)性狀，在野外調(diào)查及觀測(cè)中，常見到果實(shí)有大、中、小及長(zhǎng)、圓之分。按大小可粗略將益智果實(shí)分為大果、中果、小果，按形態(tài)可分為圓果、橢圓果，果頂端形態(tài)有鈍和尖2種形態(tài)。圓果型果實(shí)的果長(zhǎng)/果寬接近于1，但絕不小于1，果粒之間排列緊密，互相擠壓；橢圓果型果實(shí)的果長(zhǎng)/果寬大于1，果粒之間排列較為疏松，果粒之間不互相擠壓；果長(zhǎng)/果寬越大，果實(shí)越長(zhǎng)，排列越疏松。在所調(diào)查到的樣品中，果長(zhǎng)/果寬范圍為1.041～1.965，比值小的極近于球形，比值大的果長(zhǎng)較果寬長(zhǎng)近2倍。由于果頂端形態(tài)為質(zhì)量性狀，因此未計(jì)算其變異系數(shù)，3個(gè)部分平均變異系數(shù)表現(xiàn)為CV產(chǎn)量（0.464）>CV營(yíng)養(yǎng)器官（0.143）>CV果實(shí)形態(tài)（0.121），果實(shí)形態(tài)的變異較小，說明其指標(biāo)有一定的遺傳穩(wěn)定性。

2.2 益智表型性狀的相關(guān)性

72份益智樣品14個(gè)性狀間表現(xiàn)出不同程度的相關(guān)性（表3）。91對(duì)性狀組合中有34對(duì)具有顯著相關(guān)關(guān)系。單株產(chǎn)量分別與莖桿總數(shù)、結(jié)果桿數(shù)、果穗重、座果率呈極顯著正相關(guān)，果穗重分別與果序長(zhǎng)、單果大?。üL(zhǎng)×果寬）呈極顯著正相關(guān)，結(jié)果桿數(shù)分別與莖桿總數(shù)、座果率呈極顯著正相關(guān)，果序長(zhǎng)分別與單果重、單果長(zhǎng)呈極顯著正相關(guān)，單果大小分別與單果重、果序長(zhǎng)、果長(zhǎng)、果寬呈極顯著正相關(guān)，單果形態(tài)（果長(zhǎng)/果寬）分別與果頂端形態(tài)、果長(zhǎng)呈極顯著正相關(guān)，與果寬呈極顯著負(fù)相關(guān)，果頂端形態(tài)與果寬呈極顯著負(fù)相關(guān)。

2.3 主成分分析

對(duì)14個(gè)表型性狀進(jìn)行主成分分析的結(jié)果（表4）表明，前7個(gè)主成分積累貢獻(xiàn)率達(dá)到89.29%，基本可以包含14個(gè)性狀所含有的信息。第一主成分貢獻(xiàn)率為24.49%，對(duì)它影響最大的性狀包括單果大?。?.486 687）、單果重（0.451 349）、果長(zhǎng)（0.382 861）、果寬（0.340 678），該成分與單果的體積有關(guān)，可以稱為單果體積因子；第二主成分貢獻(xiàn)率為21.60%，對(duì)它影響最大的性狀包括單株產(chǎn)量（0.513 016）、果穗重（0.424 571）、結(jié)果桿數(shù)（0.389 807）、莖桿總數(shù)（0.350 040）、果序長(zhǎng)（0.295 520）、座果率（0.222 168），該成分與產(chǎn)量有關(guān)，可以稱之為產(chǎn)量因子；第三主成分貢獻(xiàn)率為18.76%，對(duì)它影響最大的性狀包括單果形態(tài)（果長(zhǎng)/果寬）（0.575 507）、果頂端形態(tài)（0.453 547）、單果寬（-0.411 844）、單果長(zhǎng)（0.368 941），該成分與形態(tài)有關(guān)，可以稱之為單果形態(tài)因子；第四、五、六、七主成分貢獻(xiàn)率較小，分別為7.23%、6.72%、5.65%、4.84%，對(duì)第四、六主成分影響最大的性狀都是株高，對(duì)第五主成分影響最大的性狀為莖桿直徑，株高與莖桿直徑都為營(yíng)養(yǎng)器官，因此可以將第四、五、六3個(gè)主成分統(tǒng)稱為營(yíng)養(yǎng)器官因子；對(duì)第七主成分影響最大的為果序長(zhǎng)，屬于產(chǎn)量因子。

2.4 益智種質(zhì)資源表型性狀的聚類分析

用SPSS19.0軟件對(duì)72個(gè)供試材料進(jìn)行聚類分析，獲得基于表型性狀的樹狀圖（圖1）。在距離值為13時(shí)，72個(gè)樣品被分為9組，其中組Ⅰ與組Ⅱ相近，組Ⅴ、組Ⅵ、組Ⅶ相近，組Ⅷ、組Ⅸ相近，各個(gè)采集地的樣品分散在各個(gè)組里。經(jīng)分析各組信息后，將單果分為大果、中果、小果，小果單果重為0.60～1.20 g；中果單果重為1.20～1.65 g，大果單果重為1.65～2.10 g；產(chǎn)量分為3個(gè)級(jí)別，果穗重與單株產(chǎn)量接近總體平均值的設(shè)為中等，遠(yuǎn)低于總體平均值的設(shè)為產(chǎn)量較低，遠(yuǎn)高于總體平均值的設(shè)為產(chǎn)量較高。據(jù)實(shí)際測(cè)量及觀察的結(jié)果可設(shè)定，當(dāng)果長(zhǎng)/果寬在1.0～1.2時(shí)，果近于圓形；當(dāng)比值在1.2～1.4時(shí)，果略呈橢圓形；當(dāng)比值在1.4～1.6時(shí)，果呈橢圓形；當(dāng)比值在大于1.6，果呈長(zhǎng)橢圓形。各組特點(diǎn)見表5。

2.5 ISSR分析

取68份益智種質(zhì)材料的分子材料并用20條引物對(duì)其進(jìn)行PCR擴(kuò)增，部分引物擴(kuò)增結(jié)果見圖2，大部分?jǐn)U增片段大小在200～1 500 bp。20個(gè)ISSR引物共擴(kuò)增出466條譜帶，其中多態(tài)性條帶為446條，多態(tài)性比例（PPB）為95.71%（表6），表明益智種質(zhì)材料具有豐富的遺傳多樣性，ISSR標(biāo)記可很好地揭示供試益智種質(zhì)材料的遺傳差異和親緣關(guān)系。采用UPGMA方法對(duì)數(shù)據(jù)進(jìn)行聚類分析，得到益智種質(zhì)材料的ISSR聚類圖（圖3），由圖3可見，68份樣品相似系數(shù)在0.72～0.86，表明群體遺傳多樣性水平不高，來自各個(gè)群體的樣本基本上按照采集地點(diǎn)聚在一起。在相似系數(shù)0.762水平上，68份樣品被分為6組。第Ⅰ組有12份樣品，包括DZ1的10份樣品和HM1的2份樣品。第Ⅱ組有18份樣品，包含HM1的7份樣品，HM2的6份樣品，DZ1的2份樣品，LM1的2份樣品，LM2的1份樣品，本組各居群樣品比較分散。第Ⅲ組有29份樣品，包含LM1的1份樣品，LM2的8份樣品，QS1的10份樣品，QS2的10份樣品。第Ⅲ組在相似系數(shù)0.772水平上又可分為4個(gè)亞組，分別為A、B、C、D 4個(gè)亞組，A亞組包含LM1的1份樣品；B亞組有14份樣品，包含LM2的8份樣品，QS2的1份樣品，QS1的5份樣品；C亞組包含QS1的3份樣品；D亞組有11份樣品，包含QS1的2份樣品，QS2的9份樣品。各居群中以LM2及QS2的樣品較為集中，QS1的樣品較為分散。第Ⅳ組有2份樣品，包含HM2的1份樣品和LM1的1份樣品。第Ⅴ組為HM2的2份樣品。第Ⅵ類為SY1的5份樣品。從總體看，各樣品主要按來源聚在一起，其中以SY1、DZ1、LM2、QS2的樣品較為集中，LM1、QS1、HM1、HM2的樣品較為分散。SY1的樣品（野生）全部集中在一個(gè)類群中，并與其他樣品（野生撫育）最早在相似系數(shù)0.72水平上被分成2大類群，可見，SY1的野生居群與其余栽培居群具有較大的遺傳差異。

3 討論與結(jié)論

在野外調(diào)查及測(cè)量中發(fā)現(xiàn)，益智各表型性狀均存在各種變異，如株高、桿圍、葉形、果序長(zhǎng)等均呈連續(xù)型的變異，較為明顯且易于分辨及歸納的為果實(shí)性狀，常見到果實(shí)有大、中、小及長(zhǎng)、圓之分，果序長(zhǎng)短差異較大，但在同一叢中各性狀表現(xiàn)均較為一致。本研究利用SPSS19.0軟件對(duì)收集到的益智形態(tài)數(shù)據(jù)進(jìn)行聚類分析，72個(gè)樣品被分為9組，將表型性狀較為接近的材料很好地聚在一起。

主成分分析可以將多個(gè)相關(guān)指標(biāo)綜合成較少的相互獨(dú)立或不相關(guān)的新指標(biāo)，并保持絕大部分的信息量。本研究經(jīng)主成分分析綜合得到7個(gè)主成分，積累貢獻(xiàn)率達(dá)到89.29%，其中前3個(gè)貢獻(xiàn)率最大的主成分積累貢獻(xiàn)率達(dá)到64.85%，單果體積因子第一，產(chǎn)量因子第二，單果形態(tài)因子第三。

本研究中，利用ISSR分子標(biāo)記對(duì)68份益智種質(zhì)材料進(jìn)行分析，采用UPGMA方法對(duì)數(shù)據(jù)進(jìn)行聚類分析，結(jié)果顯示68份樣品相似系數(shù)在0.72～0.86，來自各個(gè)群體的樣本基本上按照采集地點(diǎn)聚在一起，可以看出，ISSR分子標(biāo)記聚類具有很高的靈敏度，能夠準(zhǔn)確地將來自不同采集地點(diǎn)的樣本聚在一起。

本研究中，依據(jù)ISSR標(biāo)記的聚類與表型性狀的聚類結(jié)果不完全一致，在表型分析中散布在各組的SY1種群，在ISSR標(biāo)記聚類中全部聚成一組，而在ISSR標(biāo)記中親緣關(guān)系較近的DZ1-6和LM1-4在表型分析中則相距較遠(yuǎn)。

本研究首次將益智種質(zhì)資源的ISSR分子標(biāo)記與表型性狀的聚類進(jìn)行對(duì)比分析，旨在探索分子標(biāo)記聚類方法與表型性狀的聚類是否一致，并通過2種方法的結(jié)合對(duì)海南島益智種質(zhì)資源進(jìn)行更為合理的鑒定與分類，然而2種聚類方法的結(jié)果并不完全一致。其原因可能在于：（1）所采用的72份益智種質(zhì)資源材料均為自然授粉的地方品種， 遺傳背景較復(fù)雜；（2）表型性狀是基因與環(huán)境間、基因與基因間互作的結(jié)果，分子標(biāo)記揭示的是物種在DNA序列上的差異，因此，DNA結(jié)構(gòu)差異不一定在表型性狀上得到相應(yīng)表現(xiàn)，即使能夠表現(xiàn)，其表現(xiàn)的形式也可能因?yàn)榛蜷g的互作或調(diào)控而呈現(xiàn)多種多樣，這也是許多作物中表型分類與分子標(biāo)記分類不完全一致的主要原因[21-23]；（3）也可能是由于本研究所調(diào)查的表型性狀只是益智所有表型性狀的一部分，而ISSR標(biāo)記所分析的位點(diǎn)則覆蓋整個(gè)益智基因組；（4）試驗(yàn)研究結(jié)果也受到研究者習(xí)慣的影響，不同的研究者得出的結(jié)果可能不完全一致。同時(shí)，本研究表明，當(dāng)供試材料形態(tài)相近時(shí)，表型聚類分辨能力有限，而ISSR分子標(biāo)記聚類則具有更高的靈敏度，能夠準(zhǔn)確地將來自不同采集地點(diǎn)的樣本聚在一起。因此，運(yùn)用分子標(biāo)記技術(shù)研究益智地方群體種質(zhì)的遺傳多樣性，可以不受環(huán)境、發(fā)育時(shí)期和器官等的限制，從基因組水平上揭示其遺傳變異的程度，且在早期就可進(jìn)行新品種的鑒定。

參考文獻(xiàn)

[1] 國(guó)家藥典委員會(huì). 中華人民共和國(guó)藥典-一部-（2010年版）[S]. 北京： 中國(guó)醫(yī)藥科技出版社， 2010： 273.

[2] Itokawa H， Aiyama R， Ikuta A. A pungent diarylheptanoid from Alpinia oxyphylla[J]. Phytochemistry， 1981， 20（4）： 769-771.

[3] Itokawa H， Aiyama R， Ikuta A. A pungent principle from Alpinia oxyphylla[J]. Phytochemistry， 1982， 21（1）： 241-243.

[4] 羅秀珍， 余競(jìng)光， 徐麗珍， 等. 中藥益智化學(xué)成分的研宄[J]. 藥學(xué)學(xué)報(bào)， 2000， 35（3）： 204-207.

[5] 易美華， 薛獻(xiàn)明， 肖 紅， 等. 益智提取物對(duì)油脂抗氧化作用研究[J]. 海南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2002， 20（1）：28-32.

[6] 易美華， 肖 紅， 尹學(xué)瓊， 等. 益智提取物對(duì)超氧陰離子自由基清除作用研究[J]. 中國(guó)食品學(xué)報(bào)， 2002， 2（4）：21-24.

[7] Lin R J， Yen C M， Chou T H， et al. Antioxidant， anti-adipocyte differentiation， antitumor activity and anthelmintic activities against Anisakis simplex and Hymenolepis nana of yakuchinone A from Alpinia oxyphylla[J]. BMC Complementary and Alternative Medicine， 2013， 13： 237.

[8] 楊福孫， 甘炳春， 李榕濤， 等. 野生撫育益智主要性狀與產(chǎn)量的回歸模型及相關(guān)分析[J]. 中國(guó)農(nóng)學(xué)通報(bào)， 2010， 26（2）： 272-276.

[9] 朱文麗， 劉小濤， 莫 饒， 等. 益智的組織培養(yǎng)與快速繁殖[J]. 植物生理學(xué)通訊， 2005， 41（3）： 335.

[10] Wang H Y， Liu X J， Wen M F， et al. Analysis of the genetic diversity of natural populations of Alpinia oxyphylla Miquel using inter-simple sequence repeat markers[J]. Crop Science， 2012， 52（4）： 1 767-1 775.

[11] 劉曉靜. 益智自然種群的演化和遺傳多樣性評(píng)價(jià)[D]. ?？冢?海南大學(xué)， 2008.

[12] 陳東良， 李平顏， 許 燕， 等. 天然益智鷓鴣茶及其研制加工[J]. 熱帶農(nóng)業(yè)工程， 2008， 32（1）： 8-10.

[13] 陽(yáng)辛鳳， 利美蓮. 益智與益智酒抗氧化活性的研究[J]. 華南熱帶農(nóng)業(yè)大學(xué)學(xué)報(bào)， 2001， 7（3）： 20-23.

[14] 洪 雁， 李遠(yuǎn)志， 盧昌阜， 等. 益智果脯無硫加工工藝研究[J]. 現(xiàn)代食品科技， 2012， 28（11）： 1 526-1 529.

[15] 姜建國(guó)， 嚴(yán) 媛， 宋冬林. 鹽藻Psy基因保守序列的克隆及同源性分析[J]. 華南理工大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2004， 32（10）： 76-79.

[16] Zietkiewicz E， Rafalskia， Labuda D. Genome fingerprinting by simple sequence repeat（SSR）-anchored polymerase chain reaction amplification[J]. Genomics， 1994， 20（2）： 176-183.

[17] Hadrys H， Balick M， Schierwater B. Applications of random amplified polymorphic DNA（RAPD）in molecular ecology[J]. Molecular Ecology， 1992， 1（1）： 55-63.

[18] Shimamoto Y， Hasegawa A， Abe J， et al. Glycine soja germplasm in Japan： isozyme and chloroplast DNA variation[J]. Soybean Genetics Newsletter， 1992， 19： 73-77.

[19] Yu H， Kiang Y T. Genetic variation in South Korean natural populations of wild soybean（Glycine soja）[J]. Euphytica， 1993， 68（3）： 213-221.

[20] 黃曉丹， 張?jiān)瀑F， 應(yīng)鐵進(jìn). 高質(zhì)量植物基因組DNA的提取[J]. 植物生理學(xué)通訊， 2006， 42（2）： 311-314.

[21] 龐朝友， 杜雄明， 馬峙英. 棉屬種間雜交基因漸滲系SSR標(biāo)記及其表型性狀的聚類分析[J]. 作物學(xué)報(bào)， 2006， 32（9）： 1 371-1 378.

[22] Schut J W， Qi X， Stam P. Association between relationship measures based on AFLP markers， pedigree data and morphological traits in barley[J]. Theoretical and Applied Genetics， 1997， 95（7）： 1 161-1 168.

[23] 陳學(xué)軍， 程志芳， 陳勁楓， 等. 辣椒種質(zhì)遺傳多樣性的RAPD和ISSR及其表型數(shù)據(jù)分析[J]. 西北植物學(xué)報(bào)， 2007， 27（4）： 662-670.