摘 要：【目的】在園林綠化苗木生產(chǎn)過程中，更早篩選出具有高觀賞價(jià)值、適合在北方地區(qū)推廣種植的優(yōu)良植株，是實(shí)現(xiàn)沙棗良種選育的關(guān)鍵?！痉椒ā繉?duì)106份沙棗Elaeagnus angustifolia 2年生幼苗為試驗(yàn)材料的表型性狀、生長情況進(jìn)行連續(xù)2年的觀察記錄，同時(shí)進(jìn)行相關(guān)性分析、主成分分析和聚類分析探究其表型性狀的多樣性和觀賞特性?！窘Y(jié)果】結(jié)果表明：1）106份沙棗植株具有豐富的表型變異性，地徑年生長量的變異系數(shù)最大，為57.56%，分枝角度的多樣性指數(shù)最高，為2.023 2；2）相關(guān)性分析表明，株高年生長量與地徑年生長量存在極顯著正相關(guān)（P<0.01），葉長與株高、地徑年生長量呈顯著正相關(guān)（P<0.05）；主成分分析表明，前4個(gè)主成分累計(jì)貢獻(xiàn)率達(dá)74.449%，其中葉片大小，苗木生長勢相關(guān)性狀指標(biāo)影響較大；3）供試的106份沙棗幼苗可通過聚類分析劃分為4大類群，其中第Ⅰ類群的種質(zhì)主要特征為葉片小，生長速度緩慢；第Ⅱ類群的種質(zhì)主要特征為葉片狹長，葉表顏色偏深灰綠，生長勢居中；第Ⅲ類群的種質(zhì)主要特征為葉片大且偏卵圓形、樹冠峭立、葉色鮮亮，生長勢居中；第Ⅳ類群的種質(zhì)主要特征為葉片大，樹冠張開，生長速度極快。4）采用層次分析法綜合分析發(fā)現(xiàn)，葉色、株形等8個(gè)性狀可作為沙棗植株觀賞性綜合評(píng)價(jià)的特征性指標(biāo)。BF1、BF19、BF10為排名前3位綜合性狀表現(xiàn)優(yōu)良的沙棗單株?！窘Y(jié)論】基于沙棗表型性狀分析開展觀賞性綜合評(píng)價(jià)，為沙棗實(shí)生苗選種、早期選育沙棗優(yōu)良植株及開發(fā)利用提供理論參考。

關(guān)鍵詞：沙棗；2年生植株；表型多樣性；觀賞性評(píng)價(jià)

中圖分類號(hào)：S687.9 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1673-923X（2024）08-0036-13

基金項(xiàng)目：國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2020YFD1000500）；北京林業(yè)大學(xué)建設(shè)世界一流學(xué)科和特色發(fā)展引導(dǎo)專項(xiàng)資金項(xiàng)目（2019XKJS0324）；北京園林綠化增彩延綠科技創(chuàng)新工程（2019-KJC-02-10）；北京市園林綠化局《沙棗培育技術(shù)規(guī)程》制定標(biāo)準(zhǔn)項(xiàng)目（20211245）。

Phenotypic diversity analysis and comprehensive evaluation of ornamental properties of Elaeagnus angustifolia biennial plant

WANG Zixu， ZENG Zhihan， WANG Xiujun， LI Zihang， LI Qingwei

（College of Horticulture， Beijing Forestry University， Beijing 100083， China）

Abstract：【Objective】In the production process of providing seedlings for garden greening， the key to achieving the selection and breeding of excellent Elaeagnus angustifolia plants is to screen out plants with high ornamental value that are suitable for promotion in northern regions at an early stage.【Method】Phenotypic traits and growth conditions of 106 E. angustifolia 2-year-old seedlings were recorded for a two-year continuous observation. Correlation analysis， principal component analysis， and cluster analysis were used to explore the diversity and ornamental characteristics of their phenotypic traits.【Result】The results show： 1） The 106 E. angustifolia plants exhibited rich phenotypic variation， with the largest variation coefficient in ground diameter annual growth （57.56%） and the highest diversity index in branch angle （2.023 2）； 2） Correlation analysis showed a highly significant positive correlation （P<0.01） between plant height annual growth and ground diameter annual growth， and significant positive correlations （P<0.05） between leaf length and plant height， ground diameter annual growth； Principal component analysis indicated that the first four principal components accounted for 74.449% of the total variance， with leaf size and seedling growth vigor traits having a larger impact； 3） The 106 E. angustifolia seedlings tested could be divided into four major groups through cluster analysis， with the main characteristics of the first group being small leaves， and slow growth； the second group’s characteristics being long and narrow leaves， dark grey-green leaf colour， and moderate growth vigor； the third group’s characteristics being large and ovoid leaves， erect canopy， bright leaf color， and moderate growth vigor； the fourth group’s characteristics being large leaves， open canopy， and fast growth； 4） Using the analytic hierarchy process for comprehensive analysis， it was found that leaf color， plant shape， and other eight traits could be used as characteristic indicators for the comprehensive evaluation of the ornamental qualities of E. angustifolia. BF1， BF19， BF10 were the top three E. angustifolia individuals with excellent comprehensive trait performance.【Conclusion】A theoretical reference for the selection of E. angustifolia seedlings was provided considering to the early breeding of quality plants， and their development and utilization， based on the analysis of phenotypic traits and comprehensive evaluation of ornamental qualities.

Keywords： Elaeagnus angustifolia； biennial plant； phenotypic diversity； ornamental evaluation

沙棗Elaeagnus angustifolia為胡頹子科胡頹子屬的直立落葉喬木或小喬木[1]，具有耐旱、耐鹽堿等特點(diǎn)，是西北荒漠地區(qū)造林綠化和鹽堿地改良的優(yōu)良樹種。目前全世界大約有80多個(gè)沙棗品種，中國是沙棗主產(chǎn)地之一[2]，近年來沙棗的引種工作在山西、山東、河南、河北、北京等多個(gè)地區(qū)進(jìn)行，取得了實(shí)質(zhì)性進(jìn)展[3]。沙棗憑借其獨(dú)特的株形，顏色（葉背銀白）和氣味，成為“北京園林綠化增彩延綠科技創(chuàng)新工程”引進(jìn)的重要彩色葉樹種，其應(yīng)用有效提升了首都園林景觀的綠化層次及人居環(huán)境質(zhì)量。早期由于缺乏對(duì)不同種源地沙棗表型變異機(jī)理的系統(tǒng)了解，加上植被退化、土壤鹽漬化加重、沙漠化等環(huán)境因素，導(dǎo)致甘肅、新疆、寧夏等沙棗種源地沙棗種質(zhì)資源受到嚴(yán)重威脅。因此，更早篩選出具有高觀賞價(jià)值、適合在北方地區(qū)推廣種植的優(yōu)良植株，是科學(xué)保護(hù)沙棗種質(zhì)資源多樣性的關(guān)鍵。

作為一種高生態(tài)、高觀賞和經(jīng)濟(jì)價(jià)值的多功能樹種，有學(xué)者對(duì)沙棗進(jìn)行了廣泛的研究。這些研究主要集中在非生物脅迫[4]、藥理價(jià)值[5]、食用價(jià)值[6]及活性成分提取[7]，而對(duì)沙棗表型性狀的研究較少。例如，有學(xué)者[8]采用性狀鑒別、顯微鑒別的方法對(duì)成熟沙棗葉進(jìn)行觀察，認(rèn)為受環(huán)境影響，不同產(chǎn)地沙棗葉片形態(tài)存在著一定差異；也有前人對(duì)甘肅[9]、寧夏[10]兩地沙棗成熟植株進(jìn)行形態(tài)學(xué)研究，認(rèn)為葉片相關(guān)性狀遺傳多樣性指數(shù)高于株形。但從幼苗形態(tài)學(xué)、生長等差異入手，研究其表型變異規(guī)律及其引種選育尚未見報(bào)道。

幼苗期對(duì)逆境脅迫十分敏感，是植物生長過程中的關(guān)鍵時(shí)期，幼苗對(duì)環(huán)境的適應(yīng)性[11]直接影響了植物種質(zhì)資源保存更新以及遺傳多樣性變化，因此發(fā)掘沙棗幼苗階段的優(yōu)良表型性狀具有重要意義。近年來，有關(guān)幼苗表型性狀的研究已有一些報(bào)道，張啟雷等[12]從26個(gè)一年生奇楠沉香品系中篩選出了適合在潮州推廣種植優(yōu)良品系。文珊娜等[13]對(duì)灰木蓮幼苗葉片形態(tài)性狀進(jìn)行觀測，發(fā)現(xiàn)灰木蓮幼苗葉形態(tài)與經(jīng)緯度密切相關(guān)。對(duì)比前人研究，目前對(duì)沙棗種質(zhì)的觀賞性分析中大多未引入定性定量的評(píng)價(jià)體系[14]。層次分析法則有效的將個(gè)人主觀判斷轉(zhuǎn)化為層級(jí)客觀定量，目前已建立了桃、梅等植物的觀賞性綜合評(píng)價(jià)體系。因此，本研究以2019年9月收集的北京106份沙棗2年生幼苗為材料，首先通過分析不同種質(zhì)資源之間在生長、適應(yīng)性等早期生長差異及相關(guān)性，其次進(jìn)行主成分分析、聚類分析，最后結(jié)合層次分析法建立沙棗的觀賞性綜合評(píng)價(jià)體系，并篩選出具有高觀賞價(jià)值、適合園林應(yīng)用的沙棗優(yōu)良植株，為更早完成沙棗的良種選育提供理論依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

從甘肅地區(qū)收集同一批1年生沙棗實(shí)生苗500株，苗木規(guī)格一致，株高均為50 cm。栽植于北京市黃垡苗圃沙棗實(shí)生苗繁育基地和北京市白浮沙棗實(shí)生苗繁育基地，經(jīng)過一年的觀察記錄，最終選育出106株長勢優(yōu)良的沙棗2年生實(shí)生苗。其中北京市黃垡苗圃沙棗實(shí)生苗繁育基地（39°34′N，116°20′E）56份，編號(hào)為DX1-DX56；北京市白浮沙棗實(shí)生苗繁育基地（40°10′N，116°15′E）50份，編號(hào)為BF1-BF50（表1）。

1.2 評(píng)價(jià)指標(biāo)

1.2.1 表型性狀觀測

在沙棗植株的東南西北各個(gè)方向隨機(jī)選取10枚葉片進(jìn)行表型性狀觀測，分為定性描述性狀）和定量描述性狀。定性性狀包括葉表顏色、葉背顏色和樹冠圓整度。由于沙棗沒有種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)，因此，采用分級(jí)賦值法[15]進(jìn)行記錄（表2）。葉表和葉背顏色用英國皇家園林協(xié)會(huì)植物比色卡（RHS）測定，并記錄顏色色值編號(hào)，色值以每株樣株的10枚葉片中呈現(xiàn)最多的葉表顏色、葉背顏色值為代表。

定量性狀包括葉長（葉片基部至葉尖的距離）、葉寬（葉片最寬處的寬度）、葉形指數(shù)（葉長與葉寬的比值）、株高年生長量（在整個(gè)年齡期間每年植株頂部至地表的距離的增長量）、地徑年生長量（在整個(gè)年齡期間每年苗木距地面30 cm處直徑的增長量）和分枝角度（主枝與主干間的張開角度）。每株樣株各指標(biāo)均重復(fù)測量10次，結(jié)果取平均值。用游標(biāo)卡尺（精度0.01 cm）測量葉長和葉寬；用卷尺（精度0.1 cm）測量株高、地徑；用量角器（精度1）測量分枝角度。

1.3 層次分析法

本研究綜合考慮了沙棗幼苗觀賞性影響因素和北京地區(qū)沙棗應(yīng)用需求的實(shí)際情況，建立沙棗的綜合評(píng)價(jià)模型（表3），具體方法可參照Tang[16]。一般用CI來表示判斷矩陣的一致性，CI越小則判斷矩陣越接近完全一致性，并采用CR為一致性指標(biāo)，其值為CI/RI，若CR<0.1，則認(rèn)為該判斷矩陣合理。

1.4 數(shù)據(jù)分析

定性性狀的分布頻度以某一級(jí)別性狀出現(xiàn)的樣株數(shù)與總樣株數(shù)的百分比表示；定量性狀利用Excel 2019軟件計(jì)算各觀賞性狀指標(biāo)的最大值、最小值、平均值、標(biāo)準(zhǔn)差、變異系數(shù)（CV）和Shannon-weaver信息指數(shù)（H′）[17]。在主成分分析和聚類分析前對(duì)利用平均數(shù)（X）和標(biāo)準(zhǔn)差（S）對(duì)定量性狀進(jìn)行標(biāo)準(zhǔn)化分級(jí)處理（表4）。利用SPSS 22.0軟件，進(jìn)行主成分分析及Q型聚類分析。

2 結(jié)果與分析

2.1 沙棗種質(zhì)資源表型多樣性定性、定量分析

對(duì)供試材料定性、定量性狀的多樣性指數(shù)進(jìn)行計(jì)算（表5）。結(jié)果顯示，定性性狀的多樣性指數(shù)范圍在1.020 5～1.223 5之間，其中，葉背顏色多樣性指數(shù)最高，為1.223 5，樹冠圓整度多樣性指數(shù)最低，為1.020 5。葉背顏色多樣性指數(shù)最高是因?yàn)槿~背顏色變化大小分為5個(gè)級(jí)別，其中灰黃綠色（191B）顏色比例最大，占52.83%，灰橄欖綠色（147B）比例最小，為0.94%。然而，樹冠圓整度只有3個(gè)級(jí)別，因此該性狀多樣性指數(shù)相對(duì)較低，其中樹冠圓整的頻率最低，占16.98%，而定量性狀中分枝角度多樣性指數(shù)最大，為2.023 2，推測是分枝角度過大導(dǎo)致株形雜亂。定量性狀的多樣性指數(shù)范圍在1.917 0～2.023 2之間，其中分枝角度和地徑年生長量兩個(gè)性狀多樣性指數(shù)均超過2.000，說明在沙棗定量性狀變異中，株形多樣性程度高于葉形（表6）。沙棗的6個(gè)定量性狀均具有較高的表型多樣性，并且沙棗定量性狀多樣性指數(shù)均明顯高于定性性狀。

選取106份供試沙棗材料的6個(gè)定量性狀的變異系數(shù)進(jìn)行分析，發(fā)現(xiàn)株高、地徑年生長量的變異系數(shù)明顯大于其余性狀，變異系數(shù)分別為43.73%、56.09%；葉形指數(shù)的變異系數(shù)最小，為13.62%，變異系數(shù)范圍為13.62%～56.09%，均值為29.14%（表6）。綜合分析各性狀的變異系數(shù)、變異范圍、平均值和標(biāo)準(zhǔn)差，發(fā)現(xiàn)與葉片相關(guān)的表型性狀變異系數(shù)均較小，而與株形相關(guān)的表型性狀變異系數(shù)均較大。

2.2 沙棗表型的相關(guān)性

為了探究沙棗8個(gè)表型性狀間的關(guān)系，對(duì)所觀察的8個(gè)性狀進(jìn)行Pearson相關(guān)性分析，并繪制相關(guān)系數(shù)。結(jié)果表明，各性狀之間存在一定的相關(guān)性（表7），其中葉片長度為葉片寬度、葉形指數(shù)等其余四個(gè)性狀之間都存在顯著相關(guān)性，葉長不僅與葉寬和葉形指數(shù)呈極顯著正相關(guān)（P<0.01），且與株高、地徑年生長量呈顯著正相關(guān)（P<0.05）；葉片寬度的相關(guān)性僅次于葉片長度，其與葉形指數(shù)、分枝角度與株形均呈極顯著負(fù)相關(guān)（P<0.01），與株高年生長量均呈顯著正相關(guān)（P<0.05）。株高年生長量與地徑年生長量呈極顯著正相關(guān)（P<0.01）。因此，可根據(jù)葉片的長寬判斷株高年生長量，株高年生長量判斷地徑年生長量或通過分枝角度判斷株形趨勢。

2.3 沙棗表型的主成分分析

為了反映沙棗表型基本特征，對(duì)沙棗各性狀進(jìn)行主成分分析，結(jié)果顯示，前4個(gè)主成分的累計(jì)貢獻(xiàn)率達(dá)74.448%，并且其特征值均大于1，可以客觀反映沙棗表型性狀的基本特征（表8）。第一主成分的貢獻(xiàn)率最高，為24.201%，其中葉片寬度、葉片長度、株高年生長量的特征向量絕對(duì)值較大（大于0.53），主要反映了沙棗的葉片大小特征和苗木生長勢；第二主成分的貢獻(xiàn)率次之，為22.197%，其中葉長、地徑年生長量、分枝角度的特征向量絕對(duì)值較大（大于0.59），主要反映了沙棗的葉片大小、苗木生長勢、枝形相關(guān)特征；第三主成分的貢獻(xiàn)率偏低，為15.524%，其中株形、葉形指數(shù)的值均大于0.60，主要反映了沙棗的葉形特征和樹冠特征；第四主成分的貢獻(xiàn)率最低，為12.526%，其中葉色的值較大，為0.644，主要反映了沙棗葉片顏色的特征。

2.4 沙棗表型的聚類分析

為了明確各沙棗優(yōu)良單株的表型多樣性，采用系統(tǒng)聚類平方歐式距離法對(duì)106份沙棗種質(zhì)資源的8個(gè)表型性狀進(jìn)行聚類分析。結(jié)果顯示，在歐式距離10處，可將106株樣株分為4大類群，第Ⅰ類群包含19個(gè)樣株，第Ⅱ類群包含36個(gè)樣株，第Ⅲ類群包含16個(gè)樣株，第Ⅳ類群包含35個(gè)樣株（圖1）。比較每一類群各性狀的均值結(jié)果后，發(fā)現(xiàn)3個(gè)類群中不存在所有性狀均高于或低于其他類群的分組（表9～10）。第Ⅰ類群葉片長寬，株高、地徑年生長量的分級(jí)均值最低，而分枝角度較大，表現(xiàn)為葉片小，生長速度緩慢，可做盆景、切枝等；第Ⅱ類群葉片長度和葉表顏色的分級(jí)均值最高，表現(xiàn)為葉片狹長，葉表顏色偏深灰綠，生長勢居中，可用于觀葉；第Ⅲ類群樹冠圓整度，葉背顏色，葉片長寬的分級(jí)均值最高，而分枝角度均值最低，表現(xiàn)為葉片大且偏卵圓形、樹冠峭立、株形美觀、葉色鮮亮，生長勢居中，可用于園林綠化觀賞。第Ⅳ類群分枝角度和株高、地徑年生長量均最高，表現(xiàn)為葉片大，樹冠張開，生長速度極快，經(jīng)過適當(dāng)整形修剪后可做行道樹、園景樹（表10）。

2.5 供試沙棗樣株觀賞性評(píng)價(jià)分析

層次判斷矩陣的計(jì)算結(jié)果均顯示CR<0.10，與一致性檢驗(yàn)的要求相符（表11）。根據(jù)沙棗的自身特點(diǎn)，參考前人研究及專家的意見，建立3個(gè)約束層，葉片性狀權(quán)重系數(shù)最高（0.525 7），生長速度最低（0.141 6）。從指標(biāo)層權(quán)重結(jié)果可以看出，葉色、株形和分枝角度為權(quán)重排名前三的因子，權(quán)重分別為0.332 8、0.222 5和0.111 3。（表12）。觀測并統(tǒng)計(jì)供試植株的相關(guān)性狀并參考賦值標(biāo)準(zhǔn)進(jìn)行加權(quán)總分計(jì)算（表13），劃分為3個(gè)等級(jí)，對(duì)沙棗種質(zhì)資源進(jìn)行評(píng)價(jià)，等級(jí)越高，該沙棗單株的綜合表現(xiàn)越優(yōu)良（表14）。其中，得分≥3.3的單株樣本劃分為等級(jí)Ⅰ，此等級(jí)共有20株沙棗單株；該等級(jí)中排名第1的BF1的綜合得分為4.493 0，其株形、葉色、葉形和生長勢這4個(gè)性狀綜合觀賞性評(píng)價(jià)最優(yōu)。此等級(jí)植株的觀賞期更長，可觀賞性更強(qiáng)，即可用于景觀營造，或保留植株進(jìn)行針對(duì)性地培育及推廣。得分在2.8～3.3之間的單株樣本為等級(jí)Ⅱ，此等級(jí)共43株沙棗單株；這些品種部分性狀指標(biāo)表現(xiàn)最佳，如葉色、株形和生長勢等，但其他指標(biāo)表現(xiàn)一般，可以根據(jù)植株的突出優(yōu)勢進(jìn)行培育及應(yīng)用。得分≤2.8的43株單株樣本為等級(jí)Ⅲ，葉形、株形均表現(xiàn)一般，景觀應(yīng)用價(jià)值較低，因此可以發(fā)掘如盆景、切枝等其他方面的應(yīng)用價(jià)值。

3 討論與結(jié)論

通過對(duì)幼苗表型性狀多樣性進(jìn)行分析評(píng)價(jià)，可預(yù)測沙棗植株有價(jià)值的遺傳信息，為沙棗良種選育和遺傳改良奠定基礎(chǔ)[18]。群體內(nèi)的表型多樣性變異反映了種內(nèi)多樣性，一般情況下，變異系數(shù)大于10%，表明植物的表型變異豐富，群體對(duì)環(huán)境適應(yīng)的適應(yīng)范圍廣[19-20]。本研究中，沙棗的6個(gè)數(shù)量形狀的變異系數(shù)范圍為13.62%～57.56%，均值為29.57%，說明沙棗枝葉表型性狀變異豐富，能更好地適應(yīng)生長環(huán)境，利于自身生長、種群繁衍，這可能是沙棗除了用于觀葉外被廣泛應(yīng)用于造林綠化和鹽堿地改良的原因之一。但曾郅涵等對(duì)本引種群體進(jìn)行實(shí)生地植株的表型多樣性調(diào)查[9]，表型性狀變異系數(shù)范圍為10.29%～104.77%，遠(yuǎn)高于本研究結(jié)果，這可能是沙棗引種到適生區(qū)以外，不同環(huán)境對(duì)表型形狀產(chǎn)生顯著影響[21]，環(huán)境影響機(jī)制還有待探究。此外，沙棗多樣性指數(shù)分析表明，株型性狀的多樣性指數(shù)較葉型性狀高，除少數(shù)性狀如葉色、株形的變異系數(shù)和多樣性指數(shù)之間變化趨勢不一致外，多數(shù)性狀變異系數(shù)與多樣性指數(shù)間呈正相關(guān)。這與前人對(duì)長柄扁桃Amygdalus pedunculata的研究結(jié)論一致[22]，但與譚萬慶等[23]、馬茹暄等[24]認(rèn)為兩者之間具有負(fù)相關(guān)的關(guān)系的結(jié)論相反，具體原因有待探究。

本研究對(duì)8個(gè)性狀進(jìn)行相關(guān)性研究，結(jié)果表明存在9對(duì)顯著的相關(guān)關(guān)系，多數(shù)為葉片與葉片、株形與株形間，葉片與株形之間無顯著相關(guān)，與曾郅涵研究張掖沙棗相關(guān)性中發(fā)現(xiàn)葉寬與分枝數(shù)呈顯著正相關(guān)的結(jié)論不完全一致[9]。在這9對(duì)性狀中，株高年生長量與地徑年生長量呈極顯著正相關(guān)（P<0.01），與孫鵬等[25]對(duì)柿屬植物相關(guān)性研究的結(jié)論一致。沙棗的分枝角度與株形呈極顯著負(fù)相關(guān)，說明分枝角度過大會(huì)對(duì)樹冠圓整度產(chǎn)生不利影響。葉片相關(guān)性狀中，葉長、葉寬、葉形指數(shù)、株高、地徑年生長量之間的相關(guān)性均較高，說明在沙棗樹種資源的評(píng)價(jià)與培育的過程中，可通過調(diào)整葉形指數(shù)以選育適合觀葉的優(yōu)良單株，也可通過葉長、葉寬預(yù)測沙棗長勢，這也與陳文星等[26]對(duì)新樟屬植物葉片相關(guān)性的研究結(jié)果一致。因此，在沙棗良種早期篩選時(shí)，可通過一種性狀間接選擇另一種性狀，提高育種效率。而葉片性狀間顯著的相關(guān)性可能是由于沙棗自身多個(gè)遺傳基因具有多樣性，由這些基因共同調(diào)控定量性狀所致[26]。

表型性狀主成分分析結(jié)果表明，前4個(gè)主成分累計(jì)貢獻(xiàn)率達(dá)74.449%，4個(gè)主成分中葉片性狀、株型性狀，苗木生長勢均能反映沙棗表型性狀變異的大部分信息，其中葉片性狀對(duì)分組的貢獻(xiàn)率大于其他性狀，這與楊艷等[18]對(duì)南方黑楊種質(zhì)資源居群表型多樣性研究結(jié)果相似。聚類分析中部分依地理距離而聚類，第Ⅰ、Ⅳ類群根據(jù)試驗(yàn)區(qū)地理位置歸類，但也有極少數(shù)聚集了另一試驗(yàn)區(qū)植株，可能是由于種源間相互引種，導(dǎo)致相同試驗(yàn)區(qū)沙棗種質(zhì)資源變異較大。而Ⅱ、Ⅲ類群材料整體比較分散，根據(jù)不同性狀的相似性歸類，說明沙棗這兩個(gè)類群間表型性狀的變異是不連續(xù)的，具有隨機(jī)變異的特點(diǎn)，推測是由于經(jīng)過一年的移栽，在相似的生態(tài)環(huán)境長期作用下，沙棗表型性狀具有較強(qiáng)可塑性，兩試驗(yàn)區(qū)種源趨于相似，因而聚在一起[27]。表型性狀易受到環(huán)境因素影響，通過表型多樣性反映遺傳多樣性具有一定的局限性[19]，因此，在后續(xù)研究中需保持沙棗的種質(zhì)資源純化，以確保沙棗單株鑒別質(zhì)量。同時(shí)，進(jìn)一步擴(kuò)大沙棗種質(zhì)資源的調(diào)查范圍，引進(jìn)更多地區(qū)的優(yōu)良沙棗種質(zhì)資源，并最終建立沙棗核心種質(zhì)資源庫。結(jié)合形態(tài)標(biāo)記與生理生態(tài)特征、微觀結(jié)構(gòu)和分子標(biāo)記等手段進(jìn)行分析，以便更好地了解、開發(fā)利用和保護(hù)沙棗種質(zhì)資源。

本研究通過評(píng)價(jià)初步篩選出生長勢優(yōu)良、可觀賞性強(qiáng)、應(yīng)用前景良好的3種優(yōu)良沙棗種質(zhì)資源：BF1、BF19、BF10，建議在生產(chǎn)上進(jìn)行有針對(duì)性地培育和推廣應(yīng)用。這一評(píng)價(jià)結(jié)果為沙棗在我國北方地區(qū)開展引種造林綠化及深層次的開發(fā)利用奠定基礎(chǔ)。但是對(duì)于具有多功能的木本植物而言，僅憑借觀賞形狀來評(píng)價(jià)其在林業(yè)生產(chǎn)應(yīng)用的功能是不準(zhǔn)確的[28]。本研究僅能代表沙棗早期生長情況，后續(xù)研究中還需要對(duì)供試材料進(jìn)行長期追蹤記錄，以完善現(xiàn)有的北京地區(qū)沙棗觀賞性綜合評(píng)價(jià)體系。

在林木漫長的生長周期中高效記錄和評(píng)價(jià)優(yōu)良觀賞性狀是良種選育的關(guān)鍵。106份沙棗植株具有豐富的表型變異，本研究基于沙棗葉片性狀、生長速度和株型性狀，建立北京地區(qū)沙棗觀賞性綜合評(píng)價(jià)體系，最終篩選出具有高觀賞價(jià)值的BF1、BF19、BF10植株，適合在北京地區(qū)推廣種植。

參考文獻(xiàn)：

[1] SUN M， LIN Q. A revision of Elaeagnus L. （Elaeagnaceae） in Mainland of China[J]. Journal of Systematics and Evolution， 2010，48（5）：356-390.

[2] HAMIDPOUR R， HAMIDPOUR S， HAMIDPOUR M， et al. Russian olive （Elaeagnus angustifolia L.）： from a variety of traditional medicinal applications to its novel roles as active antioxidant， anti-inflammatory， anti-mutagenic and analgesic agent[J]. Journal of Traditional and Complementary Medicine， 2017，7（1）：24-29.

[3] 王天，林祥娜，劉軍，等.沙棗的主要活性物質(zhì)及其在食品中的應(yīng)用[J].中國果菜，2023，43（8）：42-48. WANG T， LIN X N， LIU J， et al. Main functional substances of Elaeagnus angustifolia and its application in food[J]. China Fruit Vegetable，2023，43（8）：42-48.

[4] CHANG W， ZHANG Y， PING Y， et al. Label-free quantitative proteomics of arbuscular mycorrhizal Elaeagnus angustifolia seedlings provides insights into salt-stress tolerance mechanisms[J]. Frontiers in Plant Science，2023，13：1-19.

[5] HAN J X， CHEN X Y， LIU W， et al. Triterpenoid saponin and lignan glycosides from the traditional medicine Elaeagnus angustifolia flowers and their cytotoxic activities[J]. Molecules， 2020，25（3）：1-9.

[6] WANG Y X， LI H， LI X Z， et al. Widely targeted metabolomics analysis of enriched secondary metabolites and determination of their corresponding antioxidant activities in Elaeagnus angustifolia var. orientalis （L.） Kuntze fruit juice enhanced by Bifidobacterium animalis subsp. Lactis HN-3 fermentation[J]. Food Chemistry，2022，374：1-14.

[7] AMIRI TEHRANIZADEH Z， BARATIAN A， HOSSEINZADEH H. Russian olive （Elaeagnus angustifolia） as a herbal healer[J]. Bioimpacts，2016，6（3）：155-167.

[8] 夏米斯巴奴·艾則孜，巴依爾太，孫蕓.新疆沙棗葉的鑒別與化學(xué)成分初步研究[J].新疆醫(yī)科大學(xué)學(xué)報(bào)，2019，42（4）：529-532，540. SHAMISBANU AIZEZI， BAIRTAI， SUN Y. Preliminary study on identification and chemical composition of Elaeagnus angustifolia L. leaves in Xinjiang[J]. Journal of Xinjiang Medical University，2019，42（4）：529-532，540.

[9] 曾郅涵，李慶衛(wèi)，遆羽靜.甘肅省張掖市野生沙棗的表型多樣性及觀賞性分析[J].植物資源與環(huán)境學(xué)報(bào)，2021，30（4）：12-21. ZENG Z H， LI Q W， TI Y J. Analysis on phenotypic diversity and ornamental of wild Elaeagnus angustifolia from Zhangye city of Gansu province[J]. Journal of Plant Resources and Environment， 2021，30（4）：12-21.

[10] 范義昌，柴珊珊，張曼曼，等.寧夏沙棗種質(zhì)資源表型多樣性分析[J].北方園藝，2018（23）：37-43. FAN Y C， CHAI S S， ZHANG M M， et al. Phenotypic genetic diversity of Elaeagnus angustifolia resources from Ningxia[J]. Northern Horticulture，2018（23）：37-43.

[11] SHI R R， ZHU Z， SHI N R， et al. Phenotypic diversity analysis in Elaeagnus angustifolia populations in Gansu province， China[J]. Forests，2023，14（6）：1-20.

[12] 張啟雷，劉小金，李小飛，等.奇楠沉香種質(zhì)資源評(píng)價(jià)及優(yōu)良品系早期生長分析[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2023， 43（11）：36-43. ZHANG Q L， LIU X J， LI X F， et al. Evaluation of Aquilaria sinensis ‘Qi Nan’ germplasm resources and early growth analysis of excellent clones[J]. Journal of Central South University of Forestry Technology，2023，43（11）：36-43.

[13] 文珊娜，仲崇祿，姜清彬，等.灰木蓮種源幼苗葉片性狀表型多樣性分析[J].植物研究，2017，37（2）：288-297. WEN S N， ZHONG C L， JIANG Q B， et al. Phenotypic diversity analysis of seedling leaf traits of Manglietia conifera dandy[J]. Bulletin of Botanical Research，2017，37（2）：288-297.

[14] 喬雨軒，申瀟瀟，焦雪輝，等.豫中地區(qū)46份桃種質(zhì)的觀賞性綜合評(píng)價(jià)[J].果樹學(xué)報(bào)，2024，41（2）：216-228. QIAO Y X，SHEN X X，JIAO X H，et al. Comprehensive evaluation of the ornamental value of 46 peach germplasm in the central Henan[J]. Journal of Fruit Science，2024，41（2）：216-228.

[15] WANG X L， CAO Z L， GAO C J， et al. Strategy for the construction of a core collection for Pinus yunnanensis Franch. to optimize timber based on combined phenotype and molecular marker data[J]. Genetic Resources and Crop Evolution， 2021，68（8）：3219-3240.

[16] TANG X H. Ecological garden landscape with nanotechnology design and application effect evaluation considering geotechnical characteristics[J]. Applied Nanoscience，2023，13（3）：2027-2035.

[17] NIE M， GAO L X， YAN J H， et al. Population variation of invasive Spartina alterniflora can differentiate bacterial diversity in its rhizosphere[J]. Plant Ecology，2010，209（2）：219-226.

[18] 楊艷，湯玉喜，唐潔，等.南方型黑楊種質(zhì)資源表型及生長性狀遺傳多樣性分析[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2019，39（7）： 31-36. YANG Y， TANG Y X， TANG J， et al. Genetic diversity of phenotypic and growth characters of southern type of Populus deltoides[J]. Journal of Central South University of Forestry Technology， 2019，39（7）：31-36.

[19] 吉德娟，張得芳，于倩.柴達(dá)木盆地唐古特白刺的表型多樣性[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2021，41（10）：57-66. JI D J， ZHANG D F， YU Q. Phenotypic diversity of Nitraria tangutorum Bobr. in Qaidam basin[J]. Journal of Central South University of Forestry Technology，2021，41（10）：57-66.

[20] 張強(qiáng)英，張啟翔，程堂仁.榆葉梅野生群體表型多樣性的研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2012，32（5）：155-160. ZHANG Q Y， ZHANG Q X， CHENG T R. Study on phenotypic diversity of Prunus triloba wild populations[J]. Journal of Central South University of Forestry Technology，2012，32（5）：155-160.

[21] 蔡金峰，郁萬文，汪貴斌，等.不同產(chǎn)地苦楝果實(shí)和果核表型多樣性分析[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2021，41（11）：1-9. CAI J F， YU W W， WANG G B， et al. Analysis on phenotypic variations of fruit and core traits of Melia azedarach from different seed sources[J]. Journal of Central South University of Forestry Technology，2021，41（11）：1-9.

[22] 柳江群，尹明宇，左絲雨，等.長柄扁桃天然種群表型變異[J].植物生態(tài)學(xué)報(bào)，2017，41（10）：1091-1102. LIU J Q， YIN M Y， ZUO S Y， et al. Phenotypic variations in natural populations of Amygdalus pedunculata[J]. Chinese Journal of Plant Ecology，2017，41（10）：1091-1102.

[23] 譚萬慶，彭麗平，劉政安，等.鳳丹栽培群體表型多樣性分析及綜合評(píng)價(jià)[J].經(jīng)濟(jì)林研究，2022，40（3）：180-191. TAN W Q， PENG L P， LIU Z A， et al. Phenotypic diversity analysis and comprehensive evaluation in Paeonia ostii cultivated population[J]. Non-wood Forest Research，2022，40（3）：180-191.

[24] 馬茹暄，于笑，紀(jì)若璇，等.蒙古蕕種質(zhì)資源表型性狀多樣性分析[J].干旱區(qū)資源與環(huán)境，2023，37（9）：164-171. MA R X， YU X， JI R X， et al. Diversity of phenotypic traits in Caryopteris mongholica germplasm resources[J]. Journal of Arid Land Resources and Environment，2023，37（9）：164-171.

[25] 孫鵬，傅建敏，周道順，等.柿屬植物種子表型多樣性及其萌發(fā)與幼苗生長研究[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2013，33（11）： 19-25. SUN P， FU J M， ZHOU D S， et al. Seed phenotypic diversity and germination growth of Diospyros plants[J]. Journal of Central South University of Forestry Technology，2013，33（11）：19-25.

[26] 陳文星，楊琳懿，唐軍榮，等.10個(gè)新樟屬植物葉片表型性狀多樣性分析[J].中南林業(yè)科技大學(xué)學(xué)報(bào)，2023，43（12）：83-93. CHEN W X， YANG L Y， TANG J R， et al. Analysis of the diversity on leaf phenotypic traits of 10 Neocinnamomum genus[J]. Journal of Central South University of Forestry Technology，2023，43（12）： 83-93.

[27] 楊曉霞，冷平生，鄭健，等.暴馬丁香不同種源種子和幼苗的表型性狀變異及其與地理-氣候因子的相關(guān)性[J].植物資源與環(huán)境學(xué)報(bào)，2016，25（3）：80-89. YANG X X， LENG P S， ZHENG J， et al. Variation of phenotypic traits of seed and seedling of Syringa reticulata subsp. amurensis from different provenances and their correlations with geographicclimatic factors[J]. Journal of Plant Resources and Environment， 2016，25（3）：80-89.

[28] 李冬林，梁珍海，金雅琴，等.8種梾木觀賞性及經(jīng)濟(jì)學(xué)性狀的綜合評(píng)價(jià)[J].經(jīng)濟(jì)林研究，2018，36（2）：127-132. LI D L， LIANG Z H， JIN Y Q， et al. A comprehensive evaluation on ornamental and economic characters for eight dogwoods[J]. Non-wood Forest Research，2018，36（2）：127-132.

[本文編校：羅 列]