趙 瀅，艾 軍，楊義明，王振興，劉迎雪，何 偉，劉海雙

（中國農業(yè)科學院特產研究所，長春 130112）

0 引 言

山葡萄（Vitis amurensis Rupr.）屬多年生木質藤本，是葡萄屬中最抗寒的一個種[1-2]。休眠期枝條的抗寒性具有較高的遺傳力，是葡萄抗寒育種的重要親本材料[3-4]。但田間觀察發(fā)現(xiàn)不同山葡萄種質抗寒能力存在差異，而目前適于山葡萄休眠期枝條抗寒性的鑒定方法尚不夠完善，對不同山葡萄種質資源抗寒差異的了解仍不清晰[5]，影響了抗寒種質的進一步發(fā)掘與利用。

低溫半致死溫度（LT50）可用來反映植物的抗寒性，實驗室常用相對電導率法擬合計算低溫半致死溫度來鑒定葡萄屬的抗寒性[6-7]。然而，在實際測定中發(fā)現(xiàn)由于溫度、浸泡時間以及空氣中 CO2的溶解等因素都會對實際測定的電導值產生影響，造成測定結果重復性較差[8-10]。而采用隸屬函數(shù)法對作物抗寒性進行綜合評價的結果較為全面、準確[11-14]，但對于大量種質抗寒性的篩選則需要尋找更為簡便且準確的鑒定方法。

2，3，5-氯化三苯基四氮唑（TTC）染色法也可用于作物抗寒性鑒定[15-18]，具有直觀且操作簡便等優(yōu)點，但目前關于TTC染色法測定枝條組織活力的研究報道中，均是通過觀察枝條橫截面剪口著色深淺程度進行判定[19]。在實際操作中發(fā)現(xiàn)該方法存在一些不足，主要是染色生成紅色的三苯基甲臢（TTF）經常在失去組織活性的枝條橫截面剪口外側堆積，造成觀測誤差；另外，枝條橫切面著色深淺程度目測也難以準確判斷，需量化。

本研究將以40份山葡萄種質1 a生休眠期枝條為試材，對浸染枝段進行縱切，以枝段縱切面染色面積為判斷標準，可避開原有TTC染色法橫截面上出現(xiàn)TTF堆積問題；同時探索TTC染色圖像的可視化評估，以及量化獲得的染色指數(shù)能否配合 Logistic方程在山葡萄抗寒性鑒定上進行應用，并與隸屬函數(shù)法獲得的綜合抗寒性評價結果及田間鑒定結果進行比較，以此來驗證這種以枝段縱切面染色面積為判定標準的抗寒性鑒定新方法的可行性和準確性，以期建立起適于大量山葡萄種質休眠期抗寒性鑒定的簡便、可靠方法，從而判斷不同山葡萄種質抗寒性強弱，最終為山葡萄抗寒評價研究及抗寒育種材料選擇提供有效鑒定方法和參考依據。

1 材料與方法

1.1 試驗材料

供試材料采自國家果樹種質左家山葡萄圃，種質圃地處吉林省吉林市左家鎮(zhèn)，位于北緯 44°04′，東經 126°05′，海拔190 m。年均氣溫3.6 ℃，絕對低溫-39.8 ℃。圃地為南向緩坡地，土壤為暗棕森林土。圃中每份資源為扦插繁殖的自根苗，籬架生長，株行距1.0 m × 2.5 m，柱間距6 m。樹體肥水管理采用常規(guī)管理模式[20]。40份試材包括純種山葡萄種質（V. amurensis）38份，山歐雜種（（V.amurensis×V.vinifera）×V.amurensis）2 份。38 份山葡萄種質的休眠期枝條于2017年1月剪取，當年1月最低氣溫為-28.2 ℃；山歐雜種‘北冰紅’和‘左優(yōu)紅’為對照品種，其休眠枝條于深秋10月底未埋土防寒前剪取，然后放入冷庫中待處理。所用枝條均為長勢和粗度相對一致（約0.8 cm粗）且健康成熟的1年生枝條，架面南側1.0～1.8 m范圍內取樣。每份種質取50根枝條，每根60 cm長。

1.2 材料處理

材料處理參照 Zhang等的方法[21]，略有修改。將采集的枝條帶回實驗室，剪成20 cm長的小段，裝于自封袋中，-15 ℃冰箱中貯藏。使用高低溫交變試驗箱（型號：GDJW-225，北京雅士林試驗設備有限公司）對枝條進行低溫處理，處理溫度為-20、-25、-30、-35、-40、-45、-50 ℃，按4.0 ℃/h的速度降溫，至設定溫度后保持12 h，再按4.0 ℃/h的速度升至20 ℃后，取出靜置2 h。去除表皮，避開芽眼，以備用于抗寒性評價及指標測定等研究。

1.3 抗寒性的實驗室鑒定

1.3.1 TTC染色圖像可視化評估鑒定枝條抗寒性

不同低溫處理后，從各枝條上剪取0.5 cm小段，每份種質取10段，放入10 mL 0.5 % TTC染液中，于30 ℃恒溫培養(yǎng)箱中避光染色4 d后取出，然后將枝段縱切，蒸餾水沖洗 3次，濾紙吸干水分，放入連有 WinRHIZOTM圖像分析軟件的 LA2400掃描儀中對枝段縱切面進行掃描；利用圖像分析軟件計算縱切面染色面積并確定染色等級；同時根據染色等級計算染色指數(shù)，以染色指數(shù)反映各試材的受凍程度。

染色分級標準如下（圖1）：4級，縱切面90 %以上的面積呈紅色；3級，縱切面70%～90%的面積呈紅色；2級，縱切面40%～70%的面積呈紅色；1級，縱切面5%～40%的面積呈紅色；0級，縱切面5%以下的面積呈紅色。染色指數(shù)= Σ（染色級別×枝段數(shù)）/（4×枝段總數(shù)）。

圖1 染色分級標準Fig.1 Standard of staining identification

根據不同山葡萄種質各低溫處理下枝段染色指數(shù)的變化情況，配合 Logistic方程建立回歸方程，擬合溫度拐點來確定其低溫半致死溫度（LT50），具體參照解越等方法[22]。

1.3.2 隸屬函數(shù)法綜合評價枝條抗寒性

采用隸屬函數(shù)法對山葡萄種質抗寒性進行綜合評價。根據最適冷凍溫度下枝條 7項抗寒指標（包括相對電導率、丙二醛、脯氨酸、可溶性糖、可溶性蛋白、萌枝率、凍害指數(shù)）的平均隸屬度值將山葡萄種質抗寒性分為5個級別，分級標準參照張劍俠等的方法[5]，Ⅴ級，平均隸屬度值＞ 0.70～1.00；Ⅳ級＞ 0.60～0.70；Ⅲ級＞0.40～0.60；Ⅱ級＞ 0.30～0.40；Ⅰ級 0～0.30。

1） 生理指標測定

相對電導率測定參照 Verslues等[23]的方法。丙二醛含量測定采用硫代巴比妥酸法[24]；脯氨酸含量測定采用茚三酮比色法[25]；可溶性糖含量測定采用蒽酮比色法[26]；可溶性蛋白含量測定采用考馬斯亮藍G-250染色法[27]。

2） 水培萌枝率調查

每份試材選取帶芽眼的枝段共30根，10根為1次重復，共 3次重復。經-35 ℃最適冷凍溫度處理（處理方法同1.2）后插入水中，1個月后調查水培萌枝率，參照Nanami等[28]方法。培養(yǎng)期間每 5 d換一次水。萌枝率=萌芽枝數(shù)/總枝數(shù)×100 %。

3） 凍害指數(shù)調查

枝段經最適冷凍溫度處理后，觀察橫切面組織在低溫處理后的褐變情況來確定凍害級別，凍害級別及冷害指數(shù)計算參照呂中偉等的方法[29]。

1.4 抗寒性的田間鑒定

通過調查每份試材的田間萌芽率來對其田間抗寒性進行鑒定。田間萌芽率調查參照劉崇懷等[30]方法，本試驗于2017年5月初調查每份試材在資源圃中經自然越冬后的萌芽情況，統(tǒng)計萌芽率。萌芽率=單株芽眼萌芽數(shù)/單株總芽眼數(shù)×100%。

1.5 數(shù)據統(tǒng)計分析

采用SAS9.0軟件對數(shù)據進行方程及相關性分析，利用 Duncan多重比較進行檢驗；作圖采用 Sigmaplot10.0軟件。

2 結果與分析

2.1 TTC染色圖像可視化評估鑒定山葡萄種質抗寒性

山葡萄經不同低溫處理后染色指數(shù)與冷凍溫度呈典型的倒‘S’型曲線，說明可以利用非線性回歸分析結合Logistic方程擬合溫度拐點，確定低溫半致死溫度（LT50）；低溫半致死溫度越低，抗寒性越強[31]，‘雙優(yōu)’和‘北冰紅’低溫半致死溫度分別為-38.63和-29.75 ℃（表1），可見 2個品種抗寒性差異明顯；圖像可視化評估所求得的染色指數(shù)也能將‘雙優(yōu)’和‘北冰紅’的抗寒性有效分開（圖2a）；同時不同低溫處理下枝條的水培萌芽情況與2個品種抗寒能力相符（圖2b；表2）。說明TTC染色圖像可視化評估配合 Logistic方程鑒定山葡萄種質抗寒性具有一定的可行性，40份山葡萄種質抗寒性鑒定結果詳見表1。

表1 40份山葡萄種質抗寒性鑒定結果Table 1 Identification of cold resistance in 40 V. amurensis germplasms

表2 低溫處理下2個山葡萄品種枝條水培萌枝率Table 2 Shoots sprout rate of in two varieties of V. amurensis in hydroponic culture under cold temperatures %

通過對表1中同一溫度下40份山葡萄種質的染色指數(shù)與其低溫半致死溫度的相關性分析表明，-20、-25、-30、-35和-40 ℃處理下的枝段染色指數(shù)與其低溫半致死溫度之間均呈極顯著（P＜0.000 1）負相關，其中-20 ℃下，相關系數(shù)r=-0.632 55，P＜0.000 1；-25 ℃下，r=-0.791 16，P＜0.000 1；-30 ℃下，r=-0.832 20，P＜0.000 1；-35 ℃下，r=-0.861 38，P＜0.000 1；-40 ℃下，r=-0.726 03，P＜0.000 1?？梢妰烧咴?35 ℃下的相關性最強，因此本研究將-35 ℃作為最適冷凍溫度，進一步利用隸屬函數(shù)法對最適冷凍溫度下40山葡萄種質枝條的抗寒性進行綜合評價。

圖2 染色指數(shù)及水培萌枝率在冷凍溫度下的變化Fig.2 Change of staining index and hydroponic culture sprout rate under cold temperatures

2.2 隸屬函數(shù)法綜合評價山葡萄種質的抗寒性

枝條經-35 ℃最適冷凍溫度處理后，對40份山葡萄種質抗寒指標進行相關性分析結果表明，指標變量間存在一定的相關性（表 3），說明多個指標提供的抗寒信息有交叉重疊現(xiàn)象，進一步利用主成分分析對所測指標進行綜合評價。

表3 最適冷凍溫度下山葡萄種質抗寒性指標的相關性分析Table 3 Correlation analysis of cold resistance indexes in V.amurensis germplasms under optimum freezing temperature

以累計貢獻率大于85%為原則選擇主成分，其中第1主成分中特征向量絕對值較大的是相對電導率、丙二醛及可溶性蛋白含量，貢獻率達41.48%；第2主成分則是萌枝率和凍害指數(shù)，貢獻率達25.70%；第3主成分中可溶性糖含量特征向量絕對值較大，貢獻率為15.53%；第4主成分的游離脯氨酸含量特征向量絕對值較大，貢獻率為10.90%。4個獨立主成分的累計貢獻率為93.61%（表4），利用隸屬函數(shù)法對以上指標進行綜合分析，以評價山葡萄種質的綜合抗寒性。

表4 抗寒指標主成分分析Table 4 Principal components analysis for cold resistance indexes

依據平均隸屬度值將40份山葡萄種質抗寒性分為5個級別，從Ⅴ級到Ⅰ級抗寒性依次降低。其中Ⅴ級抗寒種質1份；Ⅳ級抗寒種質10份；Ⅲ級抗寒種質24份；Ⅱ級抗寒種質3份；Ⅰ級抗寒2份為山歐雜種（表5）。結果表明，40份山葡萄種質抗寒性存在明顯差異。整體來看，純種山葡萄抗寒性 ＞ 山歐雜種（對照）。

表5 隸屬函數(shù)法綜合評價山葡萄種質抗寒性Table 5 Comprehensive analysis of cold resistances in V. amurensis germplasms

2.3 田間萌芽率

本試驗連續(xù)3 a對40份山葡萄種質自然越冬后枝條的萌芽情況進行了調查，發(fā)現(xiàn)同一資源各年份間的田間萌芽率無明顯差異。但由于田間抗寒性鑒定數(shù)據的年份要與實驗室抗寒性鑒定數(shù)據的年份相對應，所以本試驗田間萌芽率為2017年的調查結果（表6）。圖3為2017年1—5月的氣象數(shù)據，最低氣溫出現(xiàn)在1月下旬，最低為-28.2 ℃。經自然越冬后，40份山葡萄種質因抗寒性存在差異，其田間萌芽率也表現(xiàn)出不同（表6）。

2.4 相關性分析

利用SAS相關性分析結果表明，TTC染色圖像可視化評估量化獲得的染色指數(shù)配合Logistic方程求得的40份山葡萄種質低溫半致死溫度與其抗寒平均隸屬度、田間萌芽率均呈極顯著負相關（P＜0.000 1），相關系數(shù)分別為－0.883 03（圖4a）和－0.876 70（圖4b），說明該方法與隸屬函數(shù)法對山葡萄種質抗寒性鑒定結果較為一致；且染色指數(shù)擬合計算出的低溫半致死溫度與實際田間萌芽率的線性關系為 y=－5.576 5x–127.36，為驗證模型準確性，另從源圃中隨機選取了15份山葡萄種質，利用該線性關系方程計算出了 15份種質田間萌芽率的理論值（表7），與田間萌芽率的實際值相關性為－0.985 84，說明TTC染色圖像可視化評估配合Logistic方程求得的山葡萄種質低溫半致死溫度能夠反映其在田間的實際抗寒性，且建模模型準確。

表6 山葡萄種質田間萌芽率調查結果Table 6 Investigation on germination rate of V. amurensis germplasms in field

圖3 2017年氣象數(shù)據Fig.3 Meteorological data in 2017

圖4 抗寒性鑒定結果的相關性分析Fig.4 Correlation analysis of results in cold resistances identification

表7 抗寒性數(shù)據模型驗證Table 7 Proof of cold resistance data model

3 結 論

1）TTC染色枝段圖像的可視化評估實現(xiàn)了染色面積的準確量化，量化得出的染色指數(shù)與冷凍溫度呈典型的倒‘S’型曲線，說明可以利用非線性回歸分析結合Logistic 方程擬合溫度拐點，確定低溫半致死溫度。

2）該方法量化得出的染色指數(shù)配合Logistic方程確定的40份山葡萄種質低溫半致死溫度與其抗寒性綜合評價的平均隸屬度值及田間萌芽率值均呈極顯著負相關（P＜0.000 1），相關系數(shù)分別為-0.883 03和-0.876 70。

3）染色指數(shù)擬合計算出的低溫半致死溫度與實際田間萌芽率的線性關系為y= - 5.576 5 x-127.36。利用該線性關系方程計算了建模數(shù)據以外隨機選取的山葡萄種質田間萌芽率的理論值，與田間萌芽率的實際值相關性為-0.985 84，說明TTC染色圖像可視化評估配合Logistic方程求得的山葡萄種質低溫半致死溫度能夠反映其在田間的實際抗寒性。

4）40份山葡萄種質低溫半致死溫度在-39.41～-29.75 ℃范圍內，種質抗寒性存在明顯差異；總體而言，該研究的純種山葡萄低溫半致死溫度的平均值為-37.25 ℃，山歐雜種山葡萄低溫半致死溫度的平均值為-29.81 ℃，可見純種山葡萄的抗寒性強于山歐雜種山葡萄。

[1] 賀普超， 晁無疾. 我國葡萄屬野生種質資源的抗寒性分析[J]. 園藝學報，1982，9(3)：17－21.He Puchao，Chao Wuji. Studies on the cold hardiness of eight species of vitis L.growing wild in China[J]. Acta Horticulturae Sinica, 1982, 9(3): 17－21. (in Chinese with English abstract)

[2] 沈育杰，史貴文，徐浩，等. 我國山葡萄種質資源研究與利用[J]. 特產研究，1992(3)：29－32.Shen Yujie, Shi Guiwen, Xu Hao, et al. Research and utilization of Vitis amurensis germplasm resources in China[J]. Special Wild Economic Animal and Plant Research,1992(3): 29－32. (in Chinese with English abstract)

[3] Welling A, Moritz T, Palvae T, et al. Independent activation of cold acclimation by low temperature and short photoperiod in hybrid aspen[J]. Plant Physiology, 2002, 129(4): 1633－1641.

[4] Junttila O, Kaurin L. Environmental control of cold acclimation in Salix pentandra[J]. Scandinavian Journal of Forest Research, 2008, 5(1): 195－204.

[5] 張劍俠，吳行昶，楊亞州. 葡萄種間雜種的抗寒性評價[J].北方園藝，2014(13)：1－6.Zhang Jianxia, Wu Xingchang, Yang Yazhou. Cold-resistance evaluation in grape interspecific hybrids[J]. Northern Horticulture, 2014(13): 1－6. (in Chinese with English abstract)

[6] Deans J D, Billington F J, Harvey F J. Assessment of frost damage to leafless stem tissues of Quercus petraea: A reappraisal of the method of relative conductivity[J]. Forestry,1995, 68: 25－34.

[7] 張倩，劉崇懷，郭大龍，等. 5個葡萄種群的低溫半致死溫度與其抗寒適應性的關系[J]. 西北農林科技大學學報(自然科學版)，2013，41(5)：149－154.Zhang Qian, Liu Chonghuai, Guo Dalong, et al. Relationship between LT50 and cold adaptability of five grape varieties[J].Journal of Northwest A & F University (Nat. Sci. Ed.), 2013,41(5): 149－154. (in Chinese with English abstract)

[8] Zhang M I N, Willison J H M. An improved conductivity method for the measurement of frost hardiness[J]. Canadian Journal of Botany, 1987, 65: 710－715.

[9] Stergios B G, Howell G S. Evaluation of viability tests for cold stressed plants[J]. Journal of the American Society for Horticultural Science, 1973, 98(4): 325－330.

[10] 安徽農業(yè)大學. 一種茶樹抗寒性早期鑒定方法：201110045014.1[P]. 2012-08-29.

[11] 張文娥，王飛，潘學軍. 應用隸屬函數(shù)法綜合評價葡萄種間抗寒性[J]. 果樹學報，2007，24(6)：849－853.Zhang Wene, Wang Fei, Pan Yiingjun. Comprehensive evaluation on cold hardiness of Vitis species by subordinate function[J]. Journal of Fruit Science, 2007, 24(6): 849－853.(in Chinese with English abstract)

[12] 劉杜玲，張博勇，孫紅梅，等. 早實核桃不同品種抗寒性綜合評價[J]. 園藝學報，2015，42(3)：545－553.Liu Duling, Zhang Boyong, Sun Hongmei, et al. Comprehensive evaluation on drought resistance of early fruiting walnut cultivars[J]. Acta Horticulturae Sinica, 2015, 42(3): 545－553. (in Chinese with English abstract)

[13] 令凡，焦健，李朝周，等. 不同油橄欖品種對低溫脅迫的生理響應及抗寒性綜合評價[J]. 西北植物學報，2015，35(3)：508－515.Ling Fan, Jiao Jian, Li Chaozhou, et al. Physiological response and comprehensive evaluation of cold resistance under cold stress for different varieties of Olea europaea[J].Acta. Bot. Boreal-Occident Sin, 2015, 35(3): 508－515. (in Chinese with English abstract)

[14] 劉海卿，孫萬倉，劉自剛，等. 北方寒旱區(qū)白菜型冬油菜抗寒性與抗旱性評價及其關系[J]. 中國農業(yè)科學，2015，48(18)：3743－3756.Liu Haiqing, Sun Wancang, Liu Zigang, et al. Evaluation of drought resistance and cold resistance and research of their relationship at seedling stage of winter rapeseed (Brassica campestris L.) in cold and arid regions in North China[J].Scientia Agricultura Sinica, 2015, 48(18): 3743－3756. (in Chinese with English abstract)

[15] 張海晶，張金梅，辛霞，等. 馬鈴薯莖尖超低溫保存流程TTC活力響應[J]. 植物遺傳資源學報，2015，16(3)：555－560.Zhang Haijing, Zhang Jinmei, Xin Xia, et al. Study on potato shoot tips vitality responses during cryopreservation by TTC staining method[J]. Journal of Plant Genetic Resources, 2015,16 (3) :555－560. (in Chinese with English abstract)

[16] 張啟翔. 梅花及其雜交種根系抗寒性研究初報[J]. 北京林業(yè)大學學報，1992，14(增刊4)：83－86.Zhang Qixiang. A preliminary study on cold resistance of root system of Mei and its hybrid cultivars[J]. Journal of Beijing Forestry University, 1992, 14 (Supp.4): 83－86. (in Chinese with English abstract)

[17] Thomas F M, Ahlers U. Effects of excess nitrogen on frost hardiness and freezing injury of above-ground tissue in you6ng oaks (Quercus petraea and robur)[J]. New phytologist,2010, 144(1): 73－83

[18] 董麗，黃亦工，賈麥娥，等. 北京園林主要常綠闊葉植物抗凍性及其測定方法[J]. 北京林業(yè)大學學報，2002，24(3)：70－73.Dong Li, Huang Yigong, Jia Maie, et al. Freezing tolerance and methods of evaluation of introduced evergreen broad-leaf plants in open gardens in Beijing of China[J]. Journal of Beijing Forestry University, 2002, 24 (3): 70－73. (in Chinese with English abstract)

[19] 田景花，王紅霞，高儀，等. 核桃屬植物休眠期的抗寒性鑒定[J]. 園藝學報，2013，40(6)：1051－1060.Tian Jinghua, Wang Hongxia, Gao Yi, et al. Assessment of freezing tolerance of Juglans Germplasms by using annual dormant branches[J]. Acta Horticulturae Sinica, 2013, 40 (6):1051－1060. (in Chinese with English abstract)

[20] 艾軍，沈育杰. 山葡萄規(guī)范化栽培與釀酒技術[M]. 北京：中國農業(yè)出版社，2017.

[21] Zhang J，Wu X，Niu R，et al. Cold-resistance evaluation in 25 wild grape species[J]. Vitis, 2012, 51(4):153－160.

[22] 解越，張敏，梁飛俠，等. 電導法配合Logistic方程鑒定西葫蘆和黃瓜果實抗寒性[J]. 農業(yè)工程學報，2015，31(7)：292－295.Xie Yue，Zhang Min，Liang Feixia, et al. Identification of cold resistance of squash and cucumber fruits by electrical conductivity and Logistic equation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(7): 292－295. (in Chinese with English abstract)

[23] Verslues P E, Agarwal M, Katiyar-Agarwal S, et al. Methods and concepts in quantifying resistance to drought, salt and freezing, abiotic stress that affect plant water status[J]. Plant Journal, 2006, 45(4): 523－539.

[24] Yang S H, Wang J, Li S H, et al. The effects of UV-B radiation on photosynthesis in relation to Photosystem II photochemistry, thermal dissipation and antioxidant defenses in winter wheat (Triticum aestivum L.) seedlings at different growth temperatures[J]. Functional Plant Biology, 2007,34(10): 907－917.

[25] 李合生. 植物生理生化實驗原理和技術[M]. 北京：高等教育出版社，2000.

[26] 高俊鳳. 植物生理學實驗指導[M]. 北京：高等教育出版社，2006.

[27] Bradford M M. A rapid and sensitive method for the quantification of microgram quantities of protein utilizing the principle of protein-dye binding[J]. Analytical Biochemistry,1976, 72, 248－254.

[28] Nanami S, Kawaguchi H, Tateno R, et al. Sprouting traits and population structure of co-occurring Castanopsis species in an evergreen broad- leaved forest in southern China[J].Ecological Research, 2004, 19(3): 341－348.

[29] 呂中偉，王鵬，許領軍，等. 河南省葡萄凍害調查及應對措施[J]. 山西果樹，2011(1)：37－38.Lu Zhongwei, Wang Peng, Xu Lingjun, et al. Investigation and countermeasures of grape freeze injury in Henan[J].Shanxi Fruits, 2011(1): 37－38. (in Chinese with English abstract)

[30] 劉崇懷，沈育杰，陳俊. 葡萄種質資源描述規(guī)范和數(shù)據標準[M]. 北京: 中國農業(yè)出版社，2006.

[31] 于立洋，李政，韓佩堯，等. 8個新疆野蘋果優(yōu)良無性系抗寒性比較[J]. 核農學報，2017，31(9)：1827－1835.Yu Liyang, Li Zheng, Han Peiyao, et al. Comparison of cold hardiness of 8 Malus sieversii Clones[J]. Journal of Nuclear Agricultural Sciences, 2017, 31(9):1827－1835. (in Chinese with English abstract)