陳萌，魏建和，金鉞，隋春，楊成民

(中國醫(yī)學(xué)科學(xué)院 北京協(xié)和醫(yī)學(xué)院 藥用植物研究所 教育部重點(diǎn)實(shí)驗(yàn)室 瀕危藥材繁育國家工程實(shí)驗(yàn)室，北京 100193)

·中藥農(nóng)業(yè)·

應(yīng)用GGE雙標(biāo)圖分析丹參雜種一代品系的穩(wěn)定性和適應(yīng)性△

陳萌，魏建和，金鉞，隋春，楊成民*

(中國醫(yī)學(xué)科學(xué)院 北京協(xié)和醫(yī)學(xué)院 藥用植物研究所 教育部重點(diǎn)實(shí)驗(yàn)室 瀕危藥材繁育國家工程實(shí)驗(yàn)室，北京 100193)

目的：評(píng)估丹參雜種一代品系的穩(wěn)定性和適應(yīng)性。方法：采用GGE雙標(biāo)圖模型，分析5個(gè)丹參雜種一代品系和2個(gè)生產(chǎn)對(duì)照種質(zhì)開展一年多點(diǎn)試驗(yàn)數(shù)據(jù)，評(píng)價(jià)了參試材料單根重和有效成分含量的表現(xiàn)及穩(wěn)定性，同時(shí)考察了4個(gè)試驗(yàn)點(diǎn)的區(qū)分力和代表性。結(jié)果：在北京海淀，H4的單根重和丹酚酸B適應(yīng)性好，穩(wěn)定性高；在北京密云，H4和H16的單根重和有效成分的適應(yīng)性好，穩(wěn)定性高；在陜西商洛，H24的單根重和有效成分適應(yīng)性好，穩(wěn)定性高；在河南方城，H24的丹酚酸B含量高，適應(yīng)性好。通過分析環(huán)境的區(qū)分力和代表性可以確定丹參雜種一代的選育區(qū)及篩選區(qū)，北京海淀和河南方城分別作為選育高產(chǎn)、高丹酚酸B含量為目標(biāo)和高脂溶性成分為目標(biāo)的選育地點(diǎn)。品系需要在具有區(qū)分力的環(huán)境作進(jìn)一步淘汰，在北京密云淘汰產(chǎn)量低的品系，在河南方城淘汰脂溶性成分含量低的品系；在陜西商洛淘汰脂溶性成分含量低的品系。結(jié)論：應(yīng)用GGE雙標(biāo)圖評(píng)價(jià)分析方法可以較好地分析丹參雜種一代品系的穩(wěn)定性和適應(yīng)性，可為確定篩選出優(yōu)良品種、確定丹參選育試驗(yàn)地點(diǎn)和科學(xué)的建立丹參區(qū)域試驗(yàn)點(diǎn)提供科學(xué)參考。

丹參雜種一代品系；區(qū)域試驗(yàn)；適應(yīng)性；穩(wěn)定性；GGE雙標(biāo)圖

丹參為唇形科鼠尾草屬多年生草本植物丹參SalviamiltiorrhizaBge.干燥根及根莖[1]。丹參中含有水溶性酚酸類和脂溶性丹參酮類等藥效成分[2]，具有活血化瘀、鎮(zhèn)靜安神、涼血消癰等功效，尤其對(duì)心血管類疾病有顯著的療效[3]，是我國種植面積和使用量最大的常用藥材之一。但種質(zhì)混雜、無良種、各產(chǎn)地產(chǎn)量和有效成分含量差異大是丹參生產(chǎn)上面臨的主要問題。楊新杰等[4]提出了丹參質(zhì)量性狀受到遺傳和環(huán)境飾變的雙重作用，其中以遺傳主導(dǎo)型為主，環(huán)境飾變型為輔。因此，選育品質(zhì)優(yōu)良的新品種可以從種質(zhì)源頭上保障丹參生產(chǎn)的穩(wěn)產(chǎn)、高產(chǎn)、高含量，目前已選育出11個(gè)新品種[5]。新品種的穩(wěn)定性評(píng)價(jià)一直是農(nóng)作物審定和推廣的重要依據(jù)，我國的主要農(nóng)作物新品種是通過國家級(jí)或省級(jí)區(qū)域試驗(yàn)站進(jìn)行第三方評(píng)價(jià)[6]，但藥用植物新品種目前尚不具備這些條件。如何科學(xué)評(píng)價(jià)環(huán)境對(duì)藥用植物新品種(或品系)的影響及藥用植物新品種(或品系)產(chǎn)地適應(yīng)性和穩(wěn)定性，發(fā)現(xiàn)適合丹參產(chǎn)區(qū)推廣的優(yōu)良品種或品系，一直缺乏合理、規(guī)范、科學(xué)的分析方法。

基因型主效應(yīng)及其與環(huán)境互作(genotype main effect plus genotype-environment interaction，GGE)雙標(biāo)圖是研究基因型與環(huán)境互作以及不同環(huán)境下作物品種(或品系)產(chǎn)量、品質(zhì)穩(wěn)定性的新方法，該模型采用環(huán)境中心化后數(shù)據(jù)，結(jié)果中只含有與品種評(píng)價(jià)有關(guān)的G(品種效應(yīng))和GE(品種與環(huán)境互作效應(yīng))，是分析多點(diǎn)試驗(yàn)數(shù)據(jù)的理想工具。穩(wěn)定性是基于品種與環(huán)境互作效應(yīng)的大小做出的評(píng)價(jià)：互作效應(yīng)越大，說明該品種對(duì)環(huán)境越敏感，表現(xiàn)就越不穩(wěn)定，反之亦然[7]。迄今，GGE雙標(biāo)圖法已應(yīng)用于燕麥[8-9]、煙草[10]、棉花[11]、小麥[12]農(nóng)作物品種或品系的區(qū)域試驗(yàn)中，主要用于鑒別各地點(diǎn)的最佳品種及區(qū)域劃分、品種的平均表現(xiàn)及穩(wěn)定性、試驗(yàn)地點(diǎn)的區(qū)分能力和代表性等。并且GGE雙標(biāo)圖法也可應(yīng)用于一年多點(diǎn)試驗(yàn)數(shù)據(jù)分析，同樣達(dá)到較高的精確度[13-15]。目前GGE雙標(biāo)圖法尚未應(yīng)用于藥用植物的品種和品系分析中，尤其對(duì)于易受環(huán)境影響的次生代謝產(chǎn)物的適應(yīng)性和穩(wěn)定性評(píng)價(jià)應(yīng)更具意義。

課題組在丹參新品種選育過程中發(fā)現(xiàn)不同品系在不同產(chǎn)區(qū)的區(qū)域試驗(yàn)中的產(chǎn)量和有效成分含量存在顯著差異，受環(huán)境影響明顯。本研究以5個(gè)丹參F1代品系及2個(gè)生產(chǎn)對(duì)照為材料，應(yīng)用GGE-Biplot軟件，采用GGE雙標(biāo)圖模型分析了試驗(yàn)材料的表現(xiàn)，以揭示出一年多點(diǎn)試驗(yàn)不同雜種一代品系的產(chǎn)量和品質(zhì)性狀的表現(xiàn)及其穩(wěn)定性，同時(shí)評(píng)價(jià)試點(diǎn)的代表性和區(qū)分力，為篩選出優(yōu)良生態(tài)適合性的丹參雜種一代新品種，以及確定丹參選育試驗(yàn)地點(diǎn)和建立丹參區(qū)域試驗(yàn)點(diǎn)提供參考。

1 材料與方法

1.1 材料

參試材料H4、H24、H15、H16、H20，均為中國醫(yī)學(xué)科學(xué)院藥用植物研究所選育的雄性不育系DPS101與5個(gè)四代自交系雜交得到的F1代品系，其中H4和H24已于2014年分別被鑒定為新品種“中丹藥植1號(hào)”和“中丹藥植2號(hào)”[16]。對(duì)照為河南和陜西的生產(chǎn)種質(zhì)為(HNCK和SXCK)。分別在北京海淀、北京密云、河南方城和陜西商洛設(shè)置試驗(yàn)地點(diǎn)。

2013年9月在北京海淀中國醫(yī)學(xué)科學(xué)院藥用植物研究所實(shí)驗(yàn)基地統(tǒng)一育苗，2014年3月移栽到試驗(yàn)點(diǎn)；隨機(jī)區(qū)組設(shè)計(jì)，三次重復(fù)，每重復(fù)種植30株苗，行距60 cm，株距25 cm。2014年11月中旬采挖。

1.2 調(diào)查方法

每重復(fù)隨機(jī)取15株完整根系，自然陰干后稱重。所有樣品粉碎后過4號(hào)篩，采用文獻(xiàn)中的方法[17]測定丹酚酸B、丹參酮I、隱丹參酮和丹參酮ⅡA的含量。

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

1.3.1 GGE雙標(biāo)圖 研究基因型與環(huán)境互作及其不同環(huán)境下作物品種產(chǎn)量穩(wěn)定性時(shí)，某一性狀的試驗(yàn)數(shù)據(jù)為表現(xiàn)型(P)，可以分解為P=M+E+G+GE，其中M為多點(diǎn)試驗(yàn)環(huán)境總體平均值，E為環(huán)境主效應(yīng)，G為品種主效應(yīng)，GE為品種-環(huán)境互作效應(yīng)。通過數(shù)據(jù)中心化，可以從數(shù)據(jù)中移除M、E，形成一個(gè)環(huán)境中心化的兩向表。再將雙向表進(jìn)行特征值分解，得到各品種及環(huán)境的主成分表，以第一主成分(PC1)為橫坐標(biāo)，第二主成分(PC2)為縱坐標(biāo)，得到GGE雙標(biāo)圖[18]。

1.3.2分析方法 參試品系的適應(yīng)性分析：將同一方向上距離原點(diǎn)最遠(yuǎn)的品系連接起來形成一個(gè)多邊形，其他品系也包括在該多邊形內(nèi)，通過原點(diǎn)向各多邊形的各邊作垂線，稱為平均線，把多邊形分成了幾個(gè)扇形，位于頂點(diǎn)的品系是扇形區(qū)域內(nèi)所有試點(diǎn)表現(xiàn)最佳的品系。

參試品系的表現(xiàn)及其穩(wěn)定性分析：以圓圈代表平均環(huán)境，帶箭頭的直線表示平均環(huán)境軸，所指的方向?yàn)槠废翟谒协h(huán)境的近似平均表現(xiàn)；與平均環(huán)境軸垂直并通過原點(diǎn)的帶雙箭頭的直線代表品種與環(huán)境的互作的傾向性，越偏離平均環(huán)境軸表示品系越不穩(wěn)定。

試驗(yàn)點(diǎn)的代表性和區(qū)分力分析：圖中連接原點(diǎn)和環(huán)境變量的直線稱為向量，2 個(gè)環(huán)境向量間夾角的余弦近似于它們之間的遺傳相關(guān)系數(shù)，夾角小于90°表示正相關(guān)，大于90°表示負(fù)相關(guān)，接近90°表示無相關(guān)。正相關(guān)表明兩環(huán)境對(duì)品系排序相似，負(fù)相關(guān)表明兩環(huán)境對(duì)品系排序相左。環(huán)境向量的長度是試驗(yàn)點(diǎn)對(duì)品系區(qū)分能力的度量，圖中的圓圈代表平均環(huán)境，它的位置取決于各試驗(yàn)環(huán)境坐標(biāo)的平均值，帶箭頭的直線通過雙坐標(biāo)的原點(diǎn)和平均環(huán)境，試驗(yàn)點(diǎn)向量與平均環(huán)境向量的角度是其對(duì)目標(biāo)環(huán)境的代表性的度量，角度越小，代表性越強(qiáng)。

1.3.3 數(shù)據(jù)處理 數(shù)據(jù)采用Excel和SPSS 20 統(tǒng)計(jì)軟件對(duì)數(shù)據(jù)進(jìn)行平均值和方差分析。GGE雙標(biāo)圖構(gòu)建采用軟件GGE biplot(http：//www.ggebiplot.com)。

2 結(jié)果與分析

2.1 丹參參試材料的單根重和有效成分含量表現(xiàn)

2.1.1 丹酚酸B含量 從表1可知，丹酚酸B含量排名前3位的分別為H24、H4和HNCK，顯著高于其他樣品。H24和H4分別比2個(gè)對(duì)照的平均值分別提高5.91%和5.62%。4個(gè)試驗(yàn)點(diǎn)中5個(gè)品系進(jìn)行均值分析表明北京密云的丹酚酸B含量最高，其次為北京海淀。H15(在4個(gè)地區(qū))和H20(在北京海淀、密云)的丹酚酸B的含量小于3.0%，未達(dá)到《中華人民共和國藥典》標(biāo)準(zhǔn)。

2.1.2 丹參酮I、隱丹參酮和丹參酮ⅡA總含量 從表1可知，丹參酮I、隱丹參酮和丹參酮ⅡA總含量排名前3位的分別為H24、H16和H4，其中H24顯著高于其他樣品。H24、H16和H4比2個(gè)對(duì)照平均值分別提高30.74%、8.40%和-0.02%。4個(gè)試驗(yàn)點(diǎn)中5個(gè)品系進(jìn)行均值分析表明北京密云平均含量最高，其次為北京海淀。在4個(gè)地區(qū)種植的樣品均達(dá)到《中華人民共和國藥典》的標(biāo)準(zhǔn)，即丹參酮I、隱丹參酮和丹參酮ⅡA的總量不得低于0.25%[1]。

表1 丹參參試材料在各試點(diǎn)的有效成分含量

2.1.3 單根重 從表2可知，5個(gè)材料的平均單根重顯著高于2個(gè)對(duì)照。其中排名前3位的分別為H16、H4和H24，比2個(gè)對(duì)照的平均值均增產(chǎn)50.0%以上。4個(gè)試驗(yàn)點(diǎn)中，5個(gè)品系進(jìn)行均值分析表明北京海淀的平均單根重最高，其次為陜西商洛。

表2 丹參參試材料在各試點(diǎn)的單根重

2.2 基于GGE雙標(biāo)圖分析丹參參試材料的適應(yīng)性和穩(wěn)定性分析

橫坐標(biāo)代表第一主成分(PC1)，縱坐標(biāo)代表第二主成分(PC2)，兩者表示雙標(biāo)圖的擬合度，即兩向表中總變異的百分?jǐn)?shù)。PC1+PC2越大，擬合度越高，表示雙標(biāo)圖能更好的近似實(shí)際數(shù)據(jù)。當(dāng)擬合度大于80%時(shí)，得到的關(guān)系或規(guī)律接近真實(shí)[18]。以圖1A為例解讀擬合度信息：PC1解釋了81.7%數(shù)據(jù)變異平方和，PC2解釋了16.6%數(shù)據(jù)變異平方和。PC1和PC2一共解釋了98.3%的數(shù)據(jù)平方和(Sum=98.3%)。本實(shí)驗(yàn)中三個(gè)評(píng)價(jià)指標(biāo)的擬合度均在80%以上。

2.2.1 丹酚酸B含量 從圖1A適應(yīng)性分析可知，GGE雙標(biāo)圖將4個(gè)試驗(yàn)點(diǎn)分為兩組，河南方城和陜西商洛為第一組；北京海淀和北京密云為第二組。第一組中H24、SXCK表現(xiàn)最好；第二組中H4表現(xiàn)最好，H15和H20在第二組中表現(xiàn)都不好。從圖1B分析丹酚酸B含量，H24含量最高，其次是SXCK，H20和H15的含量低于總體水平。穩(wěn)定性看，HNCK、H24、H4、H16和H15的穩(wěn)定性好，H20、SXCK穩(wěn)定性差。

2.2.2 丹參酮I、隱丹參酮和丹參酮ⅡA總含量 從圖1C適應(yīng)性分析可知，GGE雙標(biāo)圖將4個(gè)試驗(yàn)點(diǎn)分為三組，北京海淀和河南方城為第一組；陜西商洛為第二組；北京密云為第三組。第一組無表現(xiàn)優(yōu)良的品系；第二組H24和HNCK表現(xiàn)最好；第三組H16表現(xiàn)最好。H15和H20在各試驗(yàn)點(diǎn)表現(xiàn)都不好。從圖1D可知，H24含量最高，其次是HNCK、SXCK、H4和H16，H15和H20的含量低于總體水平。從穩(wěn)定性看，SXCK和H20穩(wěn)定性最好，其次為H4和H24，其他材料穩(wěn)定性差。

注：A.丹參參試品系的適應(yīng)性試驗(yàn)中丹酚酸B含量；B.丹參參試品系的穩(wěn)定性試驗(yàn)中丹酚酸B含量；C.丹參參試品系的適應(yīng)性試驗(yàn)中丹參酮I、隱丹參酮和丹參酮ⅡA總含量；D.丹參參試品系的穩(wěn)定性試驗(yàn)中丹參酮I、隱丹參酮和丹參酮ⅡA總含量；E.丹參參試品系的適應(yīng)性試驗(yàn)中單根重；F.丹參參試品系的穩(wěn)定性試驗(yàn)中單根重。圖1 GGE雙標(biāo)圖分析丹參參試品系的生態(tài)適應(yīng)性和穩(wěn)定性

2.2.3 單根重 從圖1E適應(yīng)性分析可知，GGE雙標(biāo)圖將4個(gè)試驗(yàn)點(diǎn)分為兩組，北京密云、陜西商洛和北京海淀為第一組；河南方城為第二組。第一組H4和H16表現(xiàn)最好，H20和H24表現(xiàn)一般；第二組H15表現(xiàn)最好。HNCK和SXCK在各試驗(yàn)點(diǎn)表現(xiàn)都不好。從圖1F可知，H4單根重最高，其次是H16、H24、H20和H15，HNCK和SXCK的產(chǎn)量低于總體水平。從穩(wěn)定性看，H20和H24的穩(wěn)定性好，其他品系和對(duì)照的穩(wěn)定性差。

從表3可知，在北京海淀H4的單根重和丹酚酸B適應(yīng)性好，穩(wěn)定性高，適合在該地區(qū)推廣種植；在北京密云H4的單根重和有效成分的適應(yīng)性好，穩(wěn)定性高；在陜西商洛H24的單根重和有效成分適應(yīng)性好，穩(wěn)定性高；在河南方城H24適合該地區(qū)種植。

表3 基于GGE 雙標(biāo)圖的各地區(qū)產(chǎn)量和有效成分含量的表現(xiàn)

2.3 基于GGE雙標(biāo)圖分析試驗(yàn)點(diǎn)的區(qū)分力和代表性

區(qū)分力和代表性是對(duì)試驗(yàn)點(diǎn)評(píng)價(jià)的重要指標(biāo)，理想的試驗(yàn)地點(diǎn)應(yīng)當(dāng)具備兩個(gè)條件，一是對(duì)參試品種有較強(qiáng)的區(qū)分能力，二是對(duì)目標(biāo)生態(tài)區(qū)有較強(qiáng)的代表性。沒有區(qū)分能力的試驗(yàn)點(diǎn)是無用的，有區(qū)分能力但沒有代表性的試驗(yàn)點(diǎn)可用于淘汰不穩(wěn)定的品種，但不能用于選擇優(yōu)良品種。只有既有區(qū)分力又有代表性的試驗(yàn)點(diǎn)才能用來有效的選擇高產(chǎn)、高含量、穩(wěn)產(chǎn)的品種。

2.3.1 丹酚酸B含量 從圖2A可知，河南方城的區(qū)分力最強(qiáng)，其次為北京海淀。就代表性而言，4個(gè)地區(qū)間存在正相關(guān)，其中北京密云和北京海淀為緊密正相關(guān)。陜西商洛的代表性最強(qiáng)，其次為北京海淀，河南方城的最差。綜合而言，北京海淀的區(qū)分力和代表性都較好，可用于篩選優(yōu)良品系，河南方城的區(qū)分力強(qiáng)，可用于淘汰低丹酚酸B含量的品系。

2.3.2 丹參酮I、隱丹參酮和丹參酮ⅡA總含量 從圖2B可知，陜西商洛的區(qū)分力最強(qiáng)，河南方城和北京密云相當(dāng)，北京海淀最差。就代表性而言，除陜西商洛和北京密云間呈負(fù)相關(guān)外，其他兩地區(qū)間呈正相關(guān)。北京海淀的代表性最強(qiáng)，陜西商洛最差。綜合而言，河南方城的區(qū)分力和代表性都較好，可用于篩選優(yōu)良品系，陜西商洛的區(qū)分力強(qiáng)，可用于淘汰低脂溶性成分的品系。

2.3.3 單根重 從圖2C可知，北京密云的區(qū)分力最強(qiáng)，其他3個(gè)地區(qū)的區(qū)分力相當(dāng)。就代表性而言，4個(gè)地區(qū)間存在正相關(guān)，其中北京密云和陜西商洛為緊密正相關(guān)。北京海淀的代表性最強(qiáng)，河南方城最差。綜合而言，北京海淀的區(qū)分力和代表性都較好，可用于篩選優(yōu)良品系，北京密云的區(qū)分力強(qiáng)，可用于淘汰低產(chǎn)的品系。

注：A.丹酚酸B含量；B.丹參酮I、隱丹參酮和丹參酮ⅡA總含量；C.單根重。圖2 GGE雙標(biāo)圖分析試驗(yàn)點(diǎn)的代表性和區(qū)分力

3 討論

GGE雙標(biāo)圖法和傳統(tǒng)的算術(shù)平均值法的算法有所不同，相應(yīng)的結(jié)果也有區(qū)別，如SXCK的丹酚酸B含量顯著低于H4、H24和HNCK(表1)，而GGE分析結(jié)果則將SXCK排在第2位，僅次于H24(圖1B)。同樣的結(jié)果出現(xiàn)在水稻品系評(píng)價(jià)的試驗(yàn)中[19]。這是由于雙標(biāo)圖中的計(jì)算用到了參試品種在實(shí)驗(yàn)點(diǎn)的所有數(shù)據(jù)，遠(yuǎn)遠(yuǎn)多于重復(fù)值的平均值的數(shù)據(jù)，GGE雙標(biāo)圖中給出的估算值使得品種在每一地點(diǎn)的平均產(chǎn)量或含量偏向該品種總的平均值，顯得更加精確和可靠。

環(huán)境對(duì)農(nóng)作物的影響主要體現(xiàn)在產(chǎn)量上，產(chǎn)量的高低不影響農(nóng)作物的可用性。不同于農(nóng)作物，藥用植物是一種以有效成分為經(jīng)濟(jì)性狀的特殊的作物，且含量容易受到環(huán)境的影響[20]。有效成分含量是否達(dá)到《中華人民共和國藥典》的標(biāo)準(zhǔn)，直接決定藥材的可用性，因此對(duì)品系有效成分含量的評(píng)價(jià)要優(yōu)先于對(duì)產(chǎn)量的評(píng)價(jià)。嚴(yán)威凱等[18]特別指出(作物的)產(chǎn)量的穩(wěn)定性只有和高產(chǎn)相結(jié)合才有意義。這提示我們在確定各地區(qū)適宜栽培品種時(shí)，應(yīng)綜合考慮各性狀的適應(yīng)性和穩(wěn)定性，只關(guān)注適應(yīng)性或只關(guān)注穩(wěn)定性的做法是不可取的，如在河南方城和陜西商洛，SXCK的丹酚酸B的適應(yīng)性雖好(圖1A)，但穩(wěn)定性差(圖1B)，說明該種質(zhì)混雜度高；又如H20的脂溶性成分含量穩(wěn)定性好(圖2D)，但在4個(gè)試驗(yàn)點(diǎn)含量都不高(圖2C)，不適合在這4個(gè)試驗(yàn)點(diǎn)作為優(yōu)良品種種植。當(dāng)H4和H24分別種植在北京密云和陜西商洛時(shí)，單根重和有效成分含量表現(xiàn)好，穩(wěn)定性高，適合在當(dāng)?shù)赝茝V種植。同時(shí)對(duì)那些在某一特定環(huán)境下，某些性狀具有突出表現(xiàn)的品種，可以用于特殊用途的品種推廣，如H4種植在北京海淀時(shí)，單根重和丹酚酸B的表現(xiàn)好，穩(wěn)定性高，可用于水溶性成分的提取。

通過分析環(huán)境的區(qū)分力和代表性可以確定丹參選育區(qū)及區(qū)域篩選區(qū)，如選育高產(chǎn)、高丹酚酸B含量為目標(biāo)，可以將北京海淀作為選育地點(diǎn)；如選育高脂溶性成分為目標(biāo)，可以將河南方城作為選育地點(diǎn)。當(dāng)選育的雜種一代新品系需要在具有代表性的環(huán)境作進(jìn)一步淘汰不良品系時(shí)，在北京密云淘汰產(chǎn)量低的品系，在河南方城淘汰脂溶性成分含量低的品系；在陜西商洛淘汰脂溶性成分含量低的品系。值得注意的是目前的結(jié)論是基于一年多點(diǎn)的數(shù)據(jù)得出初步的提示性結(jié)果，由于品種生態(tài)劃分的結(jié)論對(duì)原種和品種推薦有長期的影響[21]，要為在不同地區(qū)建立丹參選育基地或區(qū)域試驗(yàn)基地提供可靠的科學(xué)依據(jù)，還需進(jìn)行多年試驗(yàn)數(shù)據(jù)方可確定。

[1] 國家藥典委員會(huì).中華人民共和國藥典：一部[S].北京：中國醫(yī)藥科技出版社，2015：76-77.

[2] Zhou L，Zuo Z，Chow MS.Danshen：an overview of its chemistry，pharmacology，pharmacokinetics and clinic use[J].Journal of clinical pharmacology，2005，45(12)：1345-1359.

[3] Zhou R，He LF，Li YJ，et al.Cardioprotective effect of water and ethanol extract ofSalviamiltiorrhizain an experimental model of myocardial infarction [J].Journal of ethnopharmacology，2012，139(2)：440-446.

[4] 楊新杰，萬德光，林貴兵，等.不同種源丹參異地引種后質(zhì)量分析[J].成都醫(yī)學(xué)學(xué)院學(xué)報(bào)，2011，6(4)：291-295.

[5] 楊成民，魏建和，隋春，等.我國中藥材新品種選育進(jìn)展與建議[J].中國現(xiàn)代中藥，2012，15(9)：727-737.

[6] 張群遠(yuǎn)，孔繁玲，廖琴，等.作物品種區(qū)域試驗(yàn)的評(píng)價(jià)體系及其評(píng)價(jià)方法[J].農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究，2000，16(2)：81-86.

[7] 嚴(yán)威凱，盛慶來，胡躍高，等.GGE疊圖法-分析品種×環(huán)境互作模式的理想方法[J].作物學(xué)報(bào)，2001，27(1)：21-28.

[8] 焦偉紅，劉景輝，齊冰潔，等.用GGE雙標(biāo)圖分析燕麥品種(系)農(nóng)藝和品質(zhì)性狀[J].中國農(nóng)業(yè)通報(bào)，2011，27(1)：24-29.

[9] 張志芬，付曉峰，劉俊青，等.用GGE雙標(biāo)圖分析燕麥區(qū)域試驗(yàn)品系產(chǎn)量穩(wěn)定性及試點(diǎn)代表性[J].作物學(xué)報(bào)，2010，36(8)：1377-1385.

[10] 趙潔宏，謝升東，王軼，等.GGE 雙標(biāo)圖在中間香型煙葉特色彰顯度分析中的應(yīng)用[J].中國煙草學(xué)報(bào)，2013，19(4)：28-34.

[11] 許乃銀，李健.利用GGE雙標(biāo)圖和綜合選擇指數(shù)劃分棉花品種生態(tài)區(qū)[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2014，22(9)：1113-1121.

[12] 楊進(jìn)文，朱俊剛，王曙光，等.用GGE雙標(biāo)圖及隸屬函數(shù)綜合分析山西小麥地方品種抗旱性[J].應(yīng)用生態(tài)學(xué)報(bào)，2013，24(4)：1031-1038.

[13] 王潔，廖琴，胡小軍，等.北方稻區(qū)國家水稻品種區(qū)域試驗(yàn)精確度分析[J].作物學(xué)報(bào)，2010，36(11)：1870-1876.

[14] 常磊，韓凡香，柴守璽，等.我國旱地小麥區(qū)域試驗(yàn)精確度及其環(huán)境綜合評(píng)價(jià)[J].應(yīng)用生態(tài)學(xué)報(bào)，2013，24(10)：2814-2820.

[15] 羅俊，張華，鄧祖湖，等.應(yīng)用GGE雙標(biāo)圖分析甘蔗品種(系)的產(chǎn)量和品質(zhì)性狀[J].作物學(xué)報(bào)，2013，39(1)：142-152.

[16] Chen M，Yang CM，Sui C，et al，“Zhongdanyaozhi No.1” and “Zhongdanyaozhi No.2”，are hybrid cultivars ofSalviamiltiorrhizawith high yield and active compounds content[J].PlosOne，2016，11(9)：0162691.

[17] Lu LL，Hou S，Zheng TT，et al，Analysis and comparison of the levels of bioactive components in roots ofSalviamiltiorrhizaBunge from different germplasms by high-performance liquid chromatography[J].Australian Journal of Crop Science.2013，7(1)：152-158.

[18] 嚴(yán)威凱.雙標(biāo)圖分析在農(nóng)作物品種多點(diǎn)試驗(yàn)中的應(yīng)用[J].作物學(xué)報(bào)，2010，36(11)：1805-1819.

[19] 王磊，程本義，鄂志國，等.基于GGE雙標(biāo)圖的水稻區(qū)試品種豐產(chǎn)性、穩(wěn)產(chǎn)性和適應(yīng)性評(píng)價(jià)[J].中國水稻科學(xué)，2015，29(4)：408-416.

[20] He CE，Wei JH，Jin Y，Chen SL，et al.Bioactive components of the root ofSalviamiltiorrhizae：Changes release to harvest time and germplasm line[J].Industrial Crops and Products，2010，32(3)：313-317.

[21] 許乃銀.農(nóng)作物品種試驗(yàn)數(shù)據(jù)管理與分析[M]，北京：中國農(nóng)業(yè)科學(xué)技術(shù)出版社，2016：133.

GGEBiplotsinStabilityandAdaptationEvaluationforHybridF1SalviamiltiorrhizaLines

CHENMeng，WEIJianhe，JINYue，SUIChun，YANGChenmin*

(InstituteofMedicinalPlantDevelopment，ChineseAcademyofMedicalSciencesandPekingUnionMedicalCollege，MinistryofEducationKeyLaboratory，NationalEngineeringLaboratoryofEndangeredMedicinalBreeding，Beijing100193，China)

Objective：To evaluate the stability and adaptation of single root weight and active compounds content of hybridSalviamiltiorrhizaBge F1 lines.Methods：GGE-biplot had been applied to analyze the stability of yield and quality traits，representativeness and discrimination of testing site，with 5 F1 hybrid lines and 2 production germplasm planted in four different testing sites.Results：H4 had better performance and stability at yield and salvianolic acid B content in Haidian of Beijing；H4 and H16 had better performance and stability at yield and active compounds content in Miyun of Beijing；H24 had better performance and stability at yield and active compounds content in Shangluo of Shanxi；H24 had better performance and stability at yield and lipophilic ingredients content in Fangcheng of Henan.The breeding region and selection region ofS.miltiorrhizaF1 lines could be modified through the analysis of testing site’s representativeness and discrimination.Through the analysis of testing site’s representativeness，the breeding region ofS.miltiorrhizacould be determined：Haidian of Beijing and Fangcheng of Henan had better representativeness at salvianolic acid B and lipophilic ingredients content，respectively.When planted in Miyun of Beijing cultivars with low yield could be eliminated，while Fangcheng of Henan and Shangluo of Shanxi could eliminate cultivar with low lipophilic ingredients content.Conclusion：GGE-biplot was a rational method for stability and adaptation evaluation ofS.miltiorrhizaF1 lines，and could assist us to select good culture and suitable test site.It provides scientific basis for the establishment of breeding site and selection site ofS.miltiorrhiza.

SalviamiltiorrhizaBge hybrid F1 lines；regional trial；stability；adaptation；GGE-biplot

國家自然科學(xué)基金(81303212)；高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(20131106110033)；中組部“萬人計(jì)劃”科技領(lǐng)軍人才(2015)

] 楊成民，副研究員，研究方向：中藥材栽培育種研究；E-mail：cmyang@implad.ac.cn

10.13313/j.issn.1673-4890.2017.6.015

2016-11-17)

*[