廣東黃心樹種質(zhì)表型與遺傳多樣性分析

丁釋豐 高筱鈺 吳林源 黃稚清 秦新生 馮志堅(jiān)

摘要：【目的】研究廣東地區(qū)黃心樹（Machilus gamblei）種質(zhì)資源表型與遺傳多樣性，為黃心樹種質(zhì)資源的保護(hù)、利用及新品種選育提供基礎(chǔ)數(shù)據(jù)?！痉椒ā窟x擇廣東地區(qū)具有代表性的黃心樹種群，從廣州王子山森林公園、從化石門森林公園、惠州南昆山風(fēng)景區(qū)等5個(gè)樣地取樣，測(cè)定黃心樹表型性狀指標(biāo)與RAPD分子標(biāo)記分析DNA片段多態(tài)性，探究黃心樹種質(zhì)資源表型與遺傳多樣性參數(shù)在群體間、群體內(nèi)的分布與特征?！窘Y(jié)果】測(cè)定的黃心樹15個(gè)表型性狀指標(biāo)中，有12個(gè)性狀指標(biāo)的變異系數(shù)（CV）低于其平均值，11個(gè)性狀指標(biāo)的遺傳多樣性指數(shù)（H'）高于其平均值，表明其表型性狀變異程度小，但變異類型豐富。高引物解析強(qiáng)度（Rp）引物OPAA-16、OPBH-10和OPB-17可有效評(píng)估黃心樹樣本間的遺傳多樣性。基于表型性狀對(duì)30份黃心樹樣本進(jìn)行聚類分析，以遺傳距離10為閾值，廣州王子山森林公園、從化石門森林公園與汕頭南澳島旅游區(qū)的樣本的表征親緣關(guān)系較接近。基于RAPD分子標(biāo)記對(duì)30份黃心樹樣本進(jìn)行聚類分析，在遺傳距離為15時(shí)，可將所有樣本分為5個(gè)類群，惠州南昆山風(fēng)景區(qū)、從化石門森林公園與廣州王子山森林公園的樣本遺傳距離較近，汕頭市南澳島旅游區(qū)與深圳馬巒山郊野公園的樣本也表現(xiàn)出相似的遺傳背景。對(duì)比廣東地區(qū)黃心樹表型性狀與RAPD分子標(biāo)記的聚類分析，發(fā)現(xiàn)二者類群中樣本情況一致性較差。【結(jié)論】廣東地區(qū)黃心樹種質(zhì)的表型與遺傳在種群間和種群內(nèi)部變異都較為豐富，其中惠州南昆山風(fēng)景區(qū)黃心樹種群的表型與遺傳多樣性最為豐富，樹種遺傳改良、資源開發(fā)利用等相關(guān)工作可于此處開展。

關(guān)鍵詞： 黃心樹;種質(zhì)資源;表型多樣性;遺傳多樣性;廣東

中圖分類號(hào)： S792.24? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：2095-1191（2021）05-1291-09

Abstract：【Objective】To study the phenotype and genetic diversity of Machilus gamblei King ex Hook. f. germplasm in Guangdong and provide the basic data for conservation， effective utilization of M. gamblei gemplasm and the breeding of new varieties. 【Method】Five typical populations of M. gamblei in Guangdongwere selected， such as Guangzhou Wangzishan Forest Park， Conghua Shimen National Forest Park， Huizhou Nankunshan Resort and so on. Phenotypic traits were determined and RAPD markerswere used in this study to analyze DNA fragment diversity， to explore the distribution and characteristics of phenotypicand genetic diversity between the populations and inside the populations. 【Result】Among the 15 phenotypic traits of M. gamblei， the coefficient variation（CV） of 12 phenotypic traits were lower than its average; and the genetic diversity index（H'） of 11 traits were higher than the average. The result showed that M. gamblei variation degree of phenotypic traits was low， but the variation types were abundant. The high resolving power（Rp） primers， OPAA-16， OPBH-10 and OPB-17， could effectively evaluate the genetic diversity among different samples of M. gambleiin Guangdong. Based on phenotypic traits， 30 samples of M. gamblei were clustered when genetic distance was 10， and the characterization relationship of samples which from Guangzhou Wangzishan Forest Park， Conghua Shimen Forest Park and Shantou Nanao Island Resorts was close. At the genetic distance of 15， based on the RAPD markers， 30 samples could be divided into 5 groups， samples of Huizhou Nankunshan Resorts， Shimen Forest Park and Wangzishan Forest Park had close genetic distance， samples of Nanao Island Resorts and Shenzhen Maluanshan Country Park showed similar genetic background. By comparing the cluster of phenotypic traits and RAPD markers， poor consistency of samples in the two groups was found. 【Conclusion】The phenotypic traits and genetic diversity of M. gamblei in Guangdong has rich variation between the populations and inside the populations. Samples from Huizhou Nankunshan Resort show the richest variation in both phenotypic traits and genetic diversity， which would be helpful for germplasm resources species improvement， effective utilization and other related researches.

Key words： Machilus gamblei King ex Hook. f.; germplasm resource; phenotypic diversity; genetic diversity; Guangdong

Foundation item： National Natural Science Foundation of China（31870699）; Forestry Science and Technology Innovation Project of Guangdong（2014KJCX005-02，2017KJCX024）

0 引言

【研究意義】種質(zhì)資源的遺傳多樣性是改良植物遺傳特性及培育新品種的重要基礎(chǔ)。黃心樹（Machilus gamblei King ex Hook. f.）系樟科（Lauraceae）潤(rùn)楠屬（Machilus）常綠喬木，分布于我國(guó)華南和西南地區(qū)，在印度、尼泊爾及越南北部也有分布。黃心樹樹姿優(yōu)美，花芽和新葉多彩艷麗，枝葉濃密青翠呈層狀伸展，材質(zhì)優(yōu)良且適應(yīng)性強(qiáng)，在園林、建筑、醫(yī)學(xué)和香料等領(lǐng)域具有重要的生態(tài)和經(jīng)濟(jì)價(jià)值，是開發(fā)前景較為廣闊的樹種（李陽(yáng)等，2014;Chetia et al.，2017）。近年來，人類的過度砍伐與開發(fā)，黃心樹的生境不斷遭受破壞，加之植物大小年花果情況差異大、自身結(jié)實(shí)能力較低及自然更新能力不高等自身原因，黃心樹等潤(rùn)楠屬植物生存境況令人堪憂（黃秋生等，2014）。廣東地區(qū)是黃心樹的主要產(chǎn)地之一，對(duì)該區(qū)域內(nèi)黃心樹種質(zhì)資源表型與遺傳多樣性進(jìn)行研究，有利于黃心樹種質(zhì)資源鑒定，為種質(zhì)資源保存提供遺傳背景和實(shí)驗(yàn)依據(jù)，對(duì)黃心樹種質(zhì)資源進(jìn)一步開發(fā)和地方特色品種選育工作有著重要的指導(dǎo)意義。【前人研究進(jìn)展】作為飼養(yǎng)琥珀蠶（Antheraea assama Ww.）的經(jīng)濟(jì)性作物（Chetia et al.，2017），國(guó)外對(duì)于黃心樹的相關(guān)研究開展較早，但大多集中于種植栽培方面。Das等（2010）、Chattopadhyay等（2014）、Boruah（2020）針對(duì)黃心樹炭疽病、灰枯病等病蟲害的發(fā)病原理進(jìn)行了相應(yīng)的研究。Devi等（2014）通過元素示蹤技術(shù)證實(shí)了黃心樹對(duì)于PM2.5及重金屬等有害物質(zhì)有吸收和富集作用，能夠有效凈化環(huán)境。Gogoi等（2009）對(duì)4個(gè)黃心樹同源四倍體及其二倍體親本的質(zhì)量性狀和數(shù)量性狀進(jìn)行研究，證實(shí)了四倍體性狀更具穩(wěn)定性，并提出有必要結(jié)合分子技術(shù)分析黃心樹的不同形態(tài)性狀。Rahman等（2012）應(yīng)用RAPD分子標(biāo)記對(duì)印度Assam Goalpara地區(qū)的黃心樹種質(zhì)資源遺傳多樣性進(jìn)行了研究，結(jié)果表明其遺傳多樣性較為豐富。Kumar（2016）針對(duì)黃心樹表型特征的研究，推測(cè)其表型特征的多樣性可能主要受氣候條件的影響?！颈狙芯壳腥朦c(diǎn)】目前，國(guó)內(nèi)對(duì)于黃心樹的研究較少且局限于栽培、用途及生態(tài)效益等方面（楊澤雄等，2016;馮力等，2020），在種質(zhì)的收集、識(shí)別以及評(píng)估的相關(guān)研究還處于空白狀態(tài)?！緮M解決的關(guān)鍵問題】以廣東地區(qū)黃心樹種群為試驗(yàn)材料，通過測(cè)定其表型與遺傳多樣性分析遺傳結(jié)構(gòu)的差異，探究黃心樹遺傳多樣性參數(shù)在群體間及群體內(nèi)的分布與特征，為樹種遺傳改良、種質(zhì)資源保護(hù)及資源深度開發(fā)利用提供參考依據(jù)。

1 材料與方法

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

在實(shí)地調(diào)查及2年物候觀測(cè)的基礎(chǔ)上，于2018年9月—2019年3月從廣東省廣州王子山森林公園、從化石門森林公園、惠州南昆山風(fēng)景區(qū)、深圳馬巒山郊野公園及汕頭南澳島旅游區(qū)黃心樹種群內(nèi)進(jìn)行取樣，選擇生長(zhǎng)發(fā)育正常、無嚴(yán)重缺陷和無明顯病蟲害的植株（表1）。為保證取樣的均勻性，樣株間水平距離50 m以上，海拔相差5 m以上，隨機(jī)摘取樣株樹冠中部當(dāng)年生枝5～10個(gè)和頂端輪生葉10～15片，用于表型性狀測(cè)定;于次年展葉期（2—3月），于同一樣株擇取健康幼葉若干，用密封袋貯藏，加入變色硅膠干燥，置于4 ℃低溫冷藏，用于種質(zhì)遺傳多樣性分析。

1. 2 試驗(yàn)方法

1. 2. 1 表型多樣性測(cè)定 選取黃心樹葉部表型性狀作為探討其親緣關(guān)系的性狀子集，共選擇15個(gè)表型性狀特征進(jìn)行分析（表2）。黃心樹表型性狀可分為質(zhì)量性狀和數(shù)量性狀，為準(zhǔn)確、客觀地用數(shù)碼差別來反映質(zhì)量性狀間的差異程度，同一質(zhì)量性狀有多種狀態(tài)者采用有序多態(tài)編碼，同一特征各狀態(tài)間差別基本一致者，采用順序序列間隔相等的編碼數(shù)值;數(shù)量形狀為可直接測(cè)量的量，經(jīng)掃描儀掃描后，應(yīng)用Digimizer application 4.6.0測(cè)定。

1. 2. 2 遺傳多樣性測(cè)定 （1）DNA提取。選用上海生工生物有限公司Ezup柱式植物基因組DNA抽提試劑盒進(jìn)行樣本DNA提取，試驗(yàn)操作參考說明書進(jìn)行。黃心樹葉片組織存在高濃度的酚類物質(zhì)，為去除植物組織中的多酚與多糖，按照操作說明書加入三氯甲烷之后，每管用等體積酚·三氯甲烷（pH 8.0）反復(fù)混勻，于12000 r/min離心5 min，取上層清液，反復(fù)抽提3次，再按照說明書進(jìn)行下一步操作。（2）RAPD引物篩選。隨機(jī)選取5個(gè)黃心樹樣本對(duì)Operon公司OPA、OPB、OPC、OPO、OPAA及OPBH組別的60個(gè)10堿基RAPD隨機(jī)引物進(jìn)行篩選（Rahman et al.，2012），引物由上海美吉生物醫(yī)藥科技有限公司合成。（3）PCR反應(yīng)體系與擴(kuò)增程序。PCR反應(yīng)過程中退火溫度和退火時(shí)間對(duì)于擴(kuò)增結(jié)果的影響較大（焦鋒和樓程富，2000），為使試驗(yàn)結(jié)果更準(zhǔn)確，在正式試驗(yàn)前設(shè)置不同的退火溫度和退火時(shí)間，進(jìn)行預(yù)實(shí)驗(yàn)，以確定最佳擴(kuò)增條件。通過相關(guān)預(yù)試驗(yàn)，最終決定PCR擴(kuò)增反應(yīng)在20 μL體系中進(jìn)行（10 μL 2×Taq Plus Master Mix，1.60 μL引物，1 μL樣本DNA，7.4 μL ddH2O）。PCR擴(kuò)增程序：94 ℃預(yù)變性2 min，94 ℃ 45 s，35 ℃ 45 s，72 ℃ 90 s，進(jìn)行35個(gè)循環(huán)，72 ℃延伸6 min。（4）電泳條件：采用2.0%瓊脂糖凝膠電泳分離擴(kuò)增產(chǎn)物，電壓70 V，電泳時(shí)間150 min。

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

1. 3. 1 表型多樣性數(shù)據(jù)分析 表型指標(biāo)數(shù)據(jù)用IBM SPSS Statistics 22.0進(jìn)行處理，將對(duì)原始數(shù)據(jù)矩陣進(jìn)行Z分?jǐn)?shù)標(biāo)準(zhǔn)化處理，并用標(biāo)準(zhǔn)化后的矩陣求變異系數(shù)（Coefficient of variation，CV），利用Shannon-Weaver遺傳多樣性指數(shù)（Genetic diversity index，H'）來衡量性狀遺傳多樣性大小（劉志齋等，2008;李華鋒等，2016）。

1. 3. 2 遺傳多樣性數(shù)據(jù)分析 綜合每條引物電泳的條帶信息（Bands information，iB），計(jì)算每條引物的條帶多態(tài)率（Percentage of polymorphic bands，PPB）與引物解析強(qiáng)度（Resolving power，Rp）。依據(jù)PCR擴(kuò)增產(chǎn)物電泳后的條帶，采取1/2賦值記帶，將在瓊脂糖凝膠上出現(xiàn)DNA片段的記為1，未出現(xiàn)DNA片段的記為2;所得數(shù)據(jù)運(yùn)行IBM SPSS Statistics 22.0進(jìn)行分析，采用Jaccard法進(jìn)行聚類分析，計(jì)算遺傳相似系數(shù)（Genetic similarity coefficient，GS）及遺傳距離。

2 結(jié)果與分析

2. 1 黃心樹種質(zhì)表型多樣性分析

2. 1. 1 變異系數(shù)與遺傳多樣性指數(shù)分析 廣東地區(qū)黃心樹種質(zhì)資源15個(gè)表型性狀的CV值為28.33%，其中幼枝毛被情況和冬芽芽鱗毛被情況的CV值最小，均為0%;葉尖偏向的CV值最大，其次是幼枝皮孔數(shù)量和葉面積指數(shù)，分別為118.42%、93.04%和29.59%。表型性狀的H'為1.1764，其中葉片長(zhǎng)度的H'最大，為1.6029;葉尖長(zhǎng)度、葉長(zhǎng)/葉尖、左葉基角和葉基角比（左/右）的H'也均超過1.5000;幼枝毛被情況和冬芽芽鱗毛被情況的H'最小，均為0（表3）。

2. 1. 2 基于表型性狀的聚類分析 將黃心樹15個(gè)表型性狀數(shù)據(jù)經(jīng)標(biāo)準(zhǔn)化處理后，對(duì)30份樣本進(jìn)行聚類分析。以遺傳距離15為閾值，樣本被分為3個(gè)類群，分別包含14、5和11份樣本。如圖1所示，以遺傳距離10為閾值，可將所有樣本分為5個(gè)類群，類群Ⅰ、Ⅱ、Ⅲ、Ⅳ和Ⅴ分別包含9、5、5、4和7份樣本。類群具體特點(diǎn)：類群Ⅰ中以從化石門森林公園和汕頭南澳島旅游區(qū)的樣本為主;深圳馬巒山郊野公園樣本多集中于類群Ⅲ中;類群Ⅴ中以廣州王子山森林公園和汕頭南澳島旅游區(qū)的樣本為主。此外，所有類群中均有惠州南昆山風(fēng)景區(qū)的樣本，而其他樣地的樣本分別只在3個(gè)類群中有所分布。

2. 2 引物擴(kuò)增產(chǎn)物分析

2. 2. 1 引物篩選結(jié)果 有研究表明，當(dāng)不同樣本間表型形態(tài)差異較大時(shí)，種類較少的引物就能夠評(píng)估不同個(gè)體之間的遺傳多樣性（Clegg et al.，2002）。選取WZS-1、SM-1、NKS-1、MLS-1和NAD-1等5個(gè)樣本對(duì)OPA、OPB、OPC、OPO和OPAA及OPBH組別的60個(gè)RAPD隨機(jī)引物進(jìn)行篩選，根據(jù)擴(kuò)增條帶的清晰度、穩(wěn)定性及多態(tài)性共篩選出10條RAPD引物（表4）。

2. 2. 2 引物條帶分析 廣東地區(qū)黃心樹種質(zhì)資源在表型性狀上表現(xiàn)出一定的差異性，通過RAPD分子標(biāo)記揭示其種質(zhì)遺傳物質(zhì)的多態(tài)性水平也較高。用10條RAPD引物對(duì)30份樣本進(jìn)行PCR擴(kuò)增，共計(jì)擴(kuò)增出51個(gè)條帶，大小為100～2500 bp，平均每條引物可擴(kuò)增出5.1個(gè)條帶，PPB為50.00%～85.71%（表5）。其中，多態(tài)性條帶共擴(kuò)增出34個(gè)，占總條帶數(shù)的66.67%，平均每條引物擴(kuò)增出3.4個(gè)多態(tài)性條帶，平均多態(tài)性比率達(dá)64.81%。Rp是衡量引物品種辨別能力和引物效率的重要指標(biāo)，本研究中使用的單一引物不能夠區(qū)分所有的樣本DNA的多樣性，但高Rp的有效引物（OPAA-16、OPBH-10和OPB-17）可用于評(píng)估不同樣本間的遺傳多樣性（圖2）。

2. 3 黃心樹種質(zhì)遺傳多樣性分析

2. 3. 1 遺傳相似系數(shù)分析 基于RAPD分子標(biāo)記的圖譜矩陣結(jié)果計(jì)算30份黃心樹樣本，結(jié)果顯示遺傳相似系數(shù)為0.36～0.96，平均值為0.75。遺傳相似系數(shù)矩陣熱圖（圖3）顯示了群體間NAD-6與MLS-3、NAD-5與NKS-3及WZS-5與NKS-1的遺傳相似系數(shù)最大，均為0.96，表明這些樣本的遺傳背景相似，親緣關(guān)系較近;NKS-5與MLS-4、MLS-6的GS值最小，分別為0.38和0.37，表明其親緣關(guān)系最遠(yuǎn)。此外，矩陣熱圖還顯示惠州南昆山風(fēng)景區(qū)的黃心樹樣本的遺傳相似系數(shù)平均值最低，僅0.69，表明在所有樣地中惠州南昆山風(fēng)景區(qū)的黃心樹種群遺傳多樣性程度最高;廣州王子山森林公園的黃心樹樣本的遺傳相似系數(shù)平均值最高，其次為汕頭南澳島旅游區(qū)，群體遺傳相似系數(shù)平均值分別為0.86和0.84，表明其群體內(nèi)遺傳多樣性程度較低。

2. 3. 2 基于RAPD分子標(biāo)記的聚類分析 基于RAPD分子標(biāo)記的PCR擴(kuò)增結(jié)果，對(duì)廣東地區(qū)30份黃心樹樣本進(jìn)行聚類分析（圖4）。在遺傳距離為20時(shí)，可將所有樣本分為3個(gè)類群，分別為包含12份樣本和16份樣本的2個(gè)主類群，以及包含2份樣本的次類群。在遺傳距離為15時(shí)，可將所有樣本分為5個(gè)類群，類群Ⅰ、Ⅱ、Ⅲ、Ⅳ和Ⅴ分別包含4、8、10、6和2份樣本;類群具體特點(diǎn)：類群Ⅰ以惠州南昆山風(fēng)景區(qū)、廣州王子山森林公園的樣本為主，類群Ⅱ主要包含從化石門森林公園、廣州王子山森林公園的樣本，類群Ⅲ的樣本以汕頭市南澳島旅游區(qū)、深圳馬巒山郊野公園為主。將上述結(jié)果與廣東地區(qū)黃心樹表型性狀聚類結(jié)果對(duì)比，發(fā)現(xiàn)二者類群中樣本情況不完全一致，推測(cè)其表型多樣性與遺傳多樣性發(fā)生機(jī)制不同。

3 討論

對(duì)廣東地區(qū)30份黃心樹種質(zhì)樣本的15個(gè)表型性狀的CV和H'進(jìn)行分析，結(jié)果表明30份樣本表型性狀CV與H'一致性較低。CV反映了性狀差異的離散程度，其大小反映的是性狀變異范圍大小（王海平等，2014）;H'是反映種質(zhì)資源間多樣性的另一個(gè)重要指標(biāo)，其大小反映了性狀多樣性的豐富程度（趙香娜等，2008），因此，CV和H'間不存在相關(guān)關(guān)系。植物種群的空間結(jié)構(gòu)特征意味著植物必須適應(yīng)不同的環(huán)境條件，并因生境的差異產(chǎn)生種群的生態(tài)分化，而植物營(yíng)固生長(zhǎng)的習(xí)性決定了植物在環(huán)境中的表型性狀特征具有很大的可塑性（胡啟鵬等，2008）。本研究發(fā)現(xiàn)，黃心樹15個(gè)表型性狀指標(biāo)中，有12個(gè)性狀指標(biāo)的CV低于其平均值，11個(gè)性狀指標(biāo)的H'高于其平均值;證明廣東地區(qū)黃心樹表型性狀變異程度小，但變異類型豐富，這一特性有利于黃心樹的資源配置。

本研究使用10條RAPD引物有效評(píng)估了廣東地區(qū)黃心樹種質(zhì)樣本的遺傳多樣性，其中高Rp的引物OPAA-16、OPBH-10和OPB-17可用于評(píng)估黃心樹不同樣本間的基因遺傳多樣性?；赗APD分子標(biāo)記所揭示廣東地區(qū)黃心樹種質(zhì)遺傳多樣性程度較高，但相較于RAPD分子標(biāo)記對(duì)浙江潤(rùn)楠、黑殼楠、舟山新木姜子及香樟等其他樟科樹種種質(zhì)遺傳多樣性的研究而言（宋愛云等，2003;Wang et al.，2005;張守鋒，2012;駱嘉言等，2017），廣東地區(qū)黃心樹種質(zhì)的遺傳變異程度并不算高。基于RAPD分子標(biāo)記的遺傳相似矩陣表明，惠州南昆山風(fēng)景區(qū)黃心樹種群表型及遺傳多樣性最為豐富，可為種質(zhì)資源的綜合利用、開發(fā)提供原材料，樹種遺傳改良、資源開發(fā)利用等相關(guān)工作可于此處開展。

聚類分析發(fā)現(xiàn)，廣東地區(qū)黃心樹樣本的表型與遺傳在群落間和群落內(nèi)均存在一定的多樣性。植物表型性狀特征的體現(xiàn)是生境條件與遺傳物質(zhì)共同作用的產(chǎn)物，研究發(fā)現(xiàn)生境條件類似的黃心樹種群在表型與遺傳多樣性上相似性程度較高。廣州王子山森林公園、從化石門森林公園與惠州南昆山風(fēng)景區(qū)的黃心樹樣本在表型與遺傳上的變異較其他樣地的樣本都更為接近，這3個(gè)樣地間的地理距離較為接近，并有著相似的氣候條件，結(jié)合潤(rùn)楠屬植物在我國(guó)的地理分布（胡曉敏等，2011;胡文強(qiáng)等，2017），推斷上述3個(gè)樣地的黃心樹種群親緣關(guān)系最為接近，之前或?qū)儆谕环N群，造成如今片狀分布的結(jié)果可能是由于人類活動(dòng)導(dǎo)致的限制因子或生態(tài)障礙所致。同時(shí)也發(fā)現(xiàn)，深圳馬巒山郊野公園與汕頭南澳島旅游區(qū)的黃心樹樣本的親緣關(guān)系也較為接近。現(xiàn)有的研究成果已經(jīng)肯定了不同生態(tài)小環(huán)境可導(dǎo)致不同群體遺傳結(jié)構(gòu)的顯著差異（Taylor and Aarssen，1990;Gehring and Linhart，1992）。若同種植物的生活史相同，且生活史格局具有同一性，那么小生境差異形成的自然選擇作用將進(jìn)一步引起種群間的歧化過程。汕頭南澳島旅游區(qū)和深圳馬巒山郊野公園均具有海、陸屬性動(dòng)態(tài)而復(fù)雜的自然環(huán)境，面臨著空氣、土壤中的鹽分及臺(tái)風(fēng)氣候等對(duì)植物造成的傷害（侯夢(mèng)瑩等，2019），這些遺傳上的變異雖然不一定表現(xiàn)在表型性狀特征上，但一定程度上可推斷黃心樹生態(tài)適應(yīng)的方向和途徑與其環(huán)境綜合作用中的主導(dǎo)因子密切相關(guān)。

此外，聚類分析發(fā)現(xiàn)同一樣地的黃心樹樣本并未形成單一類群，可能與其分布范圍廣、分布區(qū)內(nèi)環(huán)境因素多樣有關(guān)。地理隔離導(dǎo)致基因流動(dòng)的隔離是導(dǎo)致群體間遺傳分化的因子之一，但地理隔離在種群分化過程中的相對(duì)貢獻(xiàn)尚不明確（He et al.，2016;Zhang et al.，2016）。地理距離與群體間遺傳分化程度有所關(guān)聯(lián)，但造成特定的遺傳差異需結(jié)合各項(xiàng)生長(zhǎng)因子（葉俊偉等，2017），因此不能完全以地理距離劃分物種的遺傳變異方向異同（楊曉霞等，2016）。

4 結(jié)論

基于黃心樹表型性狀測(cè)定與RAPD分子標(biāo)記技術(shù)所揭示的廣東地區(qū)黃心樹種質(zhì)在群落間和群落內(nèi)的表型與遺傳多樣性程度較高。其中，惠州南昆山風(fēng)景區(qū)黃心樹種群表型與遺傳多樣性最為豐富，可為種質(zhì)資源的綜合利用、開發(fā)提供原材料，樹種遺傳改良、資源開發(fā)利用等相關(guān)工作可于此處開展。

參考文獻(xiàn)：

馮力，陳斯，夏尚文，王博. 2020. 雪災(zāi)干擾后林冠開闊度對(duì)黃心樹幼苗更新的影響[J]. 生態(tài)學(xué)雜志，39（3）：786-793. doi：10.13292/j.1000-4890.202003.024. [Feng L，Chen S，Xia S W，Wang B. 2020. Effects of canopy openness on seedling regeneration of Machilus bombycina after a snow damage[J]. Chinese Journal of Ecology，39（3）：786-793.]

侯夢(mèng)瑩，李芊芊，袁甜甜，王文卿，范志陽(yáng)，洪賜在. 2019. 南方濱海地區(qū)鹽霧沉降的時(shí)空分布——以福建古雷半島為例[J]. 生態(tài)學(xué)雜志，38（8）：2524-2530. doi：10.13292/j. 1000-4890.201908.020. [Hou M Y，Li Q Q，Yuan T T，Wang W Q，F(xiàn)an Z Y，Hong C Z. 2019. Spatiotemporal distribution of salt spray deposition in the coastal areas of South China： A case study of the Gulei Peninsula， Fujian[J]. Chinese Journal of Ecology，38（8）：2524-2530.]

胡啟鵬，郭志華，李春燕，馬履一. 2008. 植物表型可塑性對(duì)非生物環(huán)境因子的響應(yīng)研究進(jìn)展[J]. 林業(yè)科學(xué)，44（5）：135-142. doi：10.11707/j.1001-7488.20080525. [Hu Q P，Guo Z H，Li C Y，Ma L Y. 2008. Advance at phenotypic plasticity in plant responses to abiotic factors[J]. Scientia Silvae Sinicae，44（5）：135-142.]

胡文強(qiáng)，溫暖玲，肖玉，鐘任資，茹正忠. 2017. 短序潤(rùn)楠的地理分布及潛在分布區(qū)估計(jì)[J]. 廣東農(nóng)業(yè)科學(xué)，44（2）：82-85. doi：10.16768/j.issn.1004-874X.2017.02.013. [Hu W Q，Wen N L，Xiao Y，Zhong R Z，Ru Z Z. 2017. Geographic distribution and potential distribution estimation of Machilus breviflora[J]. Guangdong Agricultural Scien-ces，44（2）：82-85.]

胡曉敏，董安強(qiáng)，王發(fā)國(guó)，陳林，翟俊文，陳志明，邢福武. 2011. 廣東南嶺大東山浙江潤(rùn)楠群落物種多樣性與區(qū)系地理成分分析[J]. 植物科學(xué)學(xué)報(bào)，29（3）：265-271. doi：10.3724/SP.J.1142.2011.30265. [Hu X M，Dong A Q，Wang F G，Chen L，Zhai J W，Chen Z M，Xing F W. 2011. Species diversity and floristic components of Mac-hilus chekiangensis community in Dadongshan，Nanling National Nature Reserve，Guangdong[J]. Plant Science Journal，29（3）：265-271.]

黃秋生，劉光華，賈賢，廖景輝，孫京臣. 2014. 廣東省潤(rùn)楠屬植物資源現(xiàn)狀與保育對(duì)策[J]. 林業(yè)與環(huán)境科學(xué)，30（6）：72-76. doi：10.3969/j.issn.1006-4427.2014.06.014. [Huang Q S，Liu G H，Jia X，Liao J H，Sun J C. 2014. Current resource status and conservation strategy of Machilus plants in Guangdong[J]. Forestry and Environmental Science，30（6）：72-76.]

焦鋒，樓程富. 2000. RAPD技術(shù)應(yīng)用中的一些問題及對(duì)策[J]. 西北農(nóng)業(yè)學(xué)報(bào)，9（4）：98-102. doi：10.7606/j.issn.1004- 1389.2000.4.027. [Jiao F，Lou C F. 2000. Problems and its counter-measure in application of RAPD[J]. Acta Agriculturae Boreali-occidentalis Sinica，9（4）：98-102.]

李華鋒，滕杰，楊家干，楊亞杰，晏嫦妤，黃亞輝. 2016. 連南栽培型古茶樹資源葉片表型性狀遺傳多樣性及聚類分析[J]. 中國(guó)農(nóng)學(xué)通報(bào)，32（36）：109-114. doi：10.11924/j.issn. 1000-6850.casb16060018. [Li H F，Teng J，Yang J G，Yang Y J，Yan C Y，Huang Y H. 2016. Genetic diversity and cluster analysis of leaf phenotypic traits of Liannan cultivated ancient tea[J]. Chinese Agricultural Science Bulletin，32（36）：109-114.]

李陽(yáng)，江廣渝，王海洋. 2014. 4種樟科園林樹種揮發(fā)性物質(zhì)殺菌能力測(cè)定及有效成分分析[J]. 西南師范大學(xué)學(xué)報(bào)（自然科學(xué)版），39（6）：29-34. doi：10.13718/j.cnki.xsxb.2014. 06.008. [Li Y，Jiang G Y，Wang H Y. 2014. On analysis of bactericidal ability and effective component of four garden trees[J]. Journal of Southwest China Normal University（Natural Science Edition），39（6）：29-34.]

劉志齋，郭榮華，石云素，蔡一林，曹墨菊，宋燕春，王天宇，黎裕. 2008. 中國(guó)玉米地方品種核心種質(zhì)花期相關(guān)性狀的表型多樣性研究[J]. 中國(guó)農(nóng)業(yè)科學(xué)，41（6）：1591-1602. doi：10.3864/j.issn.0578-1752.2008.06.004. [Liu Z Z，Guo R H，Shi Y S，Cai Y L，Cao M J，Song Y C，Wang T Y，Li Y. 2008. Phenotypic diversity of flowering-related traits of maize landraces from the core collection preserved in China National Genebank[J]. Scientia Agricultura Sinica，41（6）：1591-1602.]

駱嘉言，王英，鐘文翰，劉海沖，蔣火俊. 2017. 樟科16個(gè)樹種木材的RAPD與ISSR分子鑒別[J]. 浙江農(nóng)林大學(xué)學(xué)報(bào)，34（5）：942-948. doi：10.11833/j.issn.2095-0756.2017.05. 023. [Luo J Y，Wang Y，Zhong W H，Liu H C，Jiang H J. 2017. Establishment of RAPD and ISSR markers for wood identification of 16 species of Lauraceae[J]. Journal of Zhejiang A & F University，34（5）：942-948.]

宋愛云，陳輝，董林水. 2003. RAPD分子標(biāo)記在鑒定香樟優(yōu)選株和普通株中的應(yīng)用[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào)，9（3）：263-265. doi：10.3321/j.issn：1006-687X.2003.03.011. [Song A Y，Chen H，Dong L S. 2003. Appliying RAPD technology to analysis of normal and excellent types of Cinnamomum camphora（L） Prest[J]. Chinese Journal of Applied and Environmental Biology，9（3）：263-265.]

王海平，李錫香，沈鏑，邱楊，宋江萍，張曉輝，Simon P W. 2014. 基于表型性狀的中國(guó)大蒜資源遺傳多樣性分析[J]. 植物遺傳資源學(xué)報(bào)，15（1）：24-31. doi：10.13430/j.cnki.jpgr.2014.01.004. [Wang H P，Li X X，Shen D，Qiu Y，Song J P，Zhang X H，Simon P W. 2014. Evaluation on genetic diversity of Garlic（Allium sativum L.） clones in China based on morphological characters[J]. Journal of Plant Genetic Resources，15（1）：24-31.]

楊曉霞，冷平生，鄭健，胡增輝，劉學(xué)婭，楊曉紅，竇德泉. 2016. 暴馬丁香不同種源種子和幼苗的表型性狀變異及其與地理—?dú)夂蛞蜃拥南嚓P(guān)性[J]. 植物資源與環(huán)境學(xué)報(bào)，25（3）：80-89. doi：10.3969/j.issn.1674-7895.2016.03. 10. [Yang X X，Leng P S，Zheng J，Hu Z H，Liu X Y，Yang X H，Dou D Q. 2016. Variation of phenotypic traits of seed and seedling of Syringa reticulata subsp. amurensis from different provenances and their correlations with geographic-climatic factors[J]. Journal of Plant Resources and Environment，25（3）：80-89.]

楊澤雄，燕騰，彭一航，鞏合德. 2016. 哀牢山黃心樹幼苗人工控制實(shí)驗(yàn)生長(zhǎng)特征[J]. 西部林業(yè)科學(xué)，45（1）：88-92. doi：10.16473/j.cnki.xblykx1972.2016.01.015. [Yang Z X，Yan T，Peng Y H，Gong H D. 2016. Seedling growth of Machilus gamblei in artificial controlled experiment in Ailao Mountain[J]. Journal of West China Forestry Scien-ce，45（1）：88-92.]

葉俊偉，張陽(yáng)，王曉娟. 2017. 中國(guó)—日本植物區(qū)系中的譜系地理間斷及其形成機(jī)制[J]. 植物生態(tài)學(xué)報(bào)，41（9）：1003-1019. doi：10.17521/cjpe.2016.0388. [Ye J W，Zhang Y，Wang X J. 2017. Phylogeographic breaks and the mechanisms of their formation in the Sino-Japanese floristic region[J]. Chinese Journal of Plant Ecology，41（9）：1003-1019.]

張守鋒. 2012. 上海地區(qū)不同表型香樟的遺傳多樣性分析[J]. 上海交通大學(xué)學(xué)報(bào)（農(nóng)業(yè)科學(xué)版），30（5）：1-4. doi：10.3969/J.ISSN.1671-9964.2012.05.001. [Zhang S F. 2012. Analysis of genetic diversity of different phenotypes Cinnamomum camphora in Shanghai[J]. Journal of Shanghai Jiaotong University（Agricultural Science），30（5）：1-4.]

趙香娜，李桂英，劉洋，岳美琪，張璞. 2008. 國(guó)內(nèi)外甜高粱種質(zhì)資源主要性狀遺傳多樣性及相關(guān)性分析[J]. 植物遺傳資源學(xué)報(bào)，（3）：302-307. doi：10.13430/j.cnki.jpgr. 2008.03.007. [Zhao X N，Li G Y，Liu Y，Yue M Q，Zhang P. 2008. Genetic diversity and correlation analysis of main agronomic characters in domestic and foreign sweet sorghum germplasm[J]. Journal of Plant Genetic Resources，（3）：302-307.]

Boruah J. 2020. Effect of 0.1% HgCl2 on surface sterilization of Som（Persea bombycina King.） explant during tissue culture—A major host plant of muga silkworm[J]. International Journal of Current Microbiology and Applied Sciences，9（7）：954-958. doi：10.20546/ijcmas.2020.907. 111.

Chattopadhyay S，Sangma C D，Tikader A，Rajan R K，Bindroo B B. 2014. Assessment of som（Persea bombycina） clones for resistance against Phyllosticta persae[J]. Journal of the Science of Food & Agriculture，39（3）：259-264. doi：10. 1590/S1982-56762014000300010.

Chetia H，Kabiraj D，Singh D，Mosahari，P V，Das S，Sharma P，Neog K，Sharma S，Jayaprakash P，Bora U. 2017. De novo transcriptome of the muga ilkworm，Antheraea assamensis（Helfer）[J]. Gene，611：54-65. doi：10.1016/j.gene.2017.02.021.

Clegg S M，Degnan S M，Kikkawa J，Moritz C，Estoup A，Owens I P F. 2002. Genetic consequences of sequential founder events by an island-colonizing bird[J]. Procee-dings of the National Academy of Sciences，99（12）：8127-8132. doi：10.1073/pnas.102583399.

Das R，Chutia M，Das K，Jha D K. 2010. Factors affecting sporulation of Pestalotiopsis disseminata causing grey blight disease of Persea bombycina Kost.，the primary food plant of muga silkworm[J]. Crop Protection，29（9）：963-968. doi：10.1016/j.cropro.2010.05.012.

Devi G，Bhattacharyya K G，Mahanta L B，Devi A. 2014. Trace metal composition of PM2.5，soil，and Machilus bombycina leaves and the effects on Antheraea assama silk worm rearing in the oil field area of Northeastern India[J]. Water Air & Soil Pollution，225（3）：1884. doi：10.1007/s11270-014-1884-2.

Gehring J L，Linhart Y B. 1992. Population structure and genetic differentiation in native and introduced populations of Deschampsia caespitosa（Poaceae.） in the colorrado alpine[J]. American Journal of Botany，79（12）：1337-1343. doi：10.2307/2445131.

Gogoi S N，Barua M，Chakravorty R. 2009. Studies on morphological characters of induced autotetraploid som[Persea bombycina（King ex. Hook f.）] Kost genotypes[J]. Indian Journal of Genetics & Plant Breeding，69（2）：145-151. doi：10.1007/s00251-009-0372-0.

He S L，Wang Y S，Li D Z，Yi T S. 2016. Environmental and historical determinants of patterns of genetic differentiation in wild soybean（Glycine soja Sieb. et Zucc）[J]. Scien-tific Reports，6：131-140. doi：10.1038/srep22795.

Kumar M. 2016. Studies on the foliar development of Som Machilus bombycena King; A primary host plant of Muga silkworm Antheraea assama Ww. in Bageshwar District of Uttarakhand[J]. International Journal of Environmental Sciences，5（4）：184-187.

Rahman A，Tanti B，Sarma G C，Kalita J. 2012. Genetic diversity of Persea gamblei from goalpara district of Assam，India[J]. Advances in Bioscience and Biotechnology，3（1）：20-24. doi：10.4236/abb.2012.31004.

Taylor D R，Aarssen L W. 1990. Complex competitive relationships among genotypes of three perennial grasses：Implications for species coexistence[J]. American Naturalist，136（3）：305-327.

Wang Z S，An S Q，Liu H，Leng X，Zheng J W，Liu Y H. 2005. Genetic structure of the endangered plant Neolitsea sericea（Lauraceae） from the Zhoushan archipelago using RAPD markers[J]. Annals of Botany，95（2）：305-313. doi：10.1093/aob/mci027.

Zhang Y H，Wang I J，Comes H P，Peng H，Qiu Y X. 2016. Contributions of historical and contemporary geographic and environmental factors to phylogeographic structure in a Tertiary relict species，Emmenopterys henryi（Rubiaceae）[J]. Scientific Reports，6：65-95. doi：10.1038/srep24041.

（責(zé)任編輯 鄧慧靈）