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

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

·中藥農業(yè)·

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

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

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

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

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

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

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

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

1 材料與方法

1.1 材料

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

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

1.2 調查方法

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

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

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

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

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

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

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

2 結果與分析

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

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

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

表1 丹參參試材料在各試點的有效成分含量

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

表2 丹參參試材料在各試點的單根重

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

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

2.2.1 丹酚酸B含量 從圖1A適應性分析可知，GGE雙標圖將4個試驗點分為兩組，河南方城和陜西商洛為第一組；北京海淀和北京密云為第二組。第一組中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適應性分析可知，GGE雙標圖將4個試驗點分為三組，北京海淀和河南方城為第一組；陜西商洛為第二組；北京密云為第三組。第一組無表現(xiàn)優(yōu)良的品系；第二組H24和HNCK表現(xiàn)最好；第三組H16表現(xiàn)最好。H15和H20在各試驗點表現(xiàn)都不好。從圖1D可知，H24含量最高，其次是HNCK、SXCK、H4和H16，H15和H20的含量低于總體水平。從穩(wěn)定性看，SXCK和H20穩(wěn)定性最好，其次為H4和H24，其他材料穩(wěn)定性差。

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

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

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

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

2.3 基于GGE雙標圖分析試驗點的區(qū)分力和代表性

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

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

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

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

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

3 討論

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

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

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

[1] 國家藥典委員會.中華人民共和國藥典：一部[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] 楊新杰，萬德光，林貴兵，等.不同種源丹參異地引種后質量分析[J].成都醫(yī)學學院學報，2011，6(4)：291-295.

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

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

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

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

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

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

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

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

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

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

[15] 羅俊，張華，鄧祖湖，等.應用GGE雙標圖分析甘蔗品種(系)的產量和品質性狀[J].作物學報，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] 嚴威凱.雙標圖分析在農作物品種多點試驗中的應用[J].作物學報，2010，36(11)：1805-1819.

[19] 王磊，程本義，鄂志國，等.基于GGE雙標圖的水稻區(qū)試品種豐產性、穩(wěn)產性和適應性評價[J].中國水稻科學，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] 許乃銀.農作物品種試驗數(shù)據(jù)管理與分析[M]，北京：中國農業(yè)科學技術出版社，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

國家自然科學基金(81303212)；高等學校博士學科點專項科研基金(20131106110033)；中組部“萬人計劃”科技領軍人才(2015)

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

10.13313/j.issn.1673-4890.2017.6.015

2016-11-17)

*[