何 鑫，閆向前，孫華軍，張 琪

（河南省商丘市農(nóng)林科學院，河南 商丘 476000）

大豆新品種商豆1310的豐產(chǎn)性、穩(wěn)定性及適應(yīng)性分析*

何 鑫，閆向前，孫華軍，張 琪

（河南省商丘市農(nóng)林科學院，河南 商丘 476000）

為了更全面地了解大豆新品種商豆1310的生產(chǎn)特性，以2015年和2016年國家黃淮海大豆區(qū)域試驗結(jié)果為依據(jù)，運用變異系數(shù)法、高穩(wěn)系數(shù)法、回歸系數(shù)法及AMMI模型對商豆1310大豆品種的豐產(chǎn)性、穩(wěn)產(chǎn)性及適應(yīng)性進行分析。結(jié)果表明，商豆1310是一個產(chǎn)量高、豐產(chǎn)性和穩(wěn)產(chǎn)性均好，適應(yīng)范圍廣的品種，可以在生產(chǎn)上大面積推廣應(yīng)用。

大豆品種；商豆1310；豐產(chǎn)性；穩(wěn)產(chǎn)性；適應(yīng)性

高產(chǎn)大豆新品種商豆1310由河南省商丘市農(nóng)林科學院選育。該品種生育期101 d，有限結(jié)莢習性，卵圓葉，紫花，灰毛；株高78.8 cm，株型收斂，主莖16.0節(jié)，有效分枝數(shù)2.1個；單株有效莢數(shù)53.6個，單株粒數(shù)106.4粒，單株粒重19.8 g，百粒重19.2 g，籽粒橢圓形、黃色、微光澤、淺褐色臍；籽粒粗蛋白含量42.11%，粗脂肪含量20.49%。落葉性好，不裂莢，倒伏級別1級，穩(wěn)產(chǎn)性較好，豐產(chǎn)性好。

接種鑒定，抗花葉病毒病3號株系和7號株系，高感胞囊線蟲病1號生理小種。2015年參加黃淮海南區(qū)國家大豆品種區(qū)域試驗，平均產(chǎn)量3 574.05 kg/hm2，較對照中黃13增產(chǎn)10.07%，增產(chǎn)極顯著，居參試品種第1位。2016年續(xù)試，平均產(chǎn)量3 105.55 kg/hm2，較對照中黃13增產(chǎn)12.51%，增產(chǎn)極顯著，居參試品種第1位。2016年參加生產(chǎn)試驗，平均產(chǎn)量達到3 010.30 kg/hm2，7個試點均較對照中黃13增產(chǎn)，增產(chǎn)幅度為2.72%～24.64%，平均較對照中黃13增產(chǎn)11.22%，居生產(chǎn)試驗第1位。2017年通過國家品種審定委員會初審。

大豆品種的選育以豐產(chǎn)性作為主要鑒定指標，同時還要堅持在高產(chǎn)前提下的穩(wěn)產(chǎn)性及適應(yīng)性。作物新品種在生產(chǎn)上的豐產(chǎn)性、穩(wěn)產(chǎn)性和廣泛的適應(yīng)性是檢驗品種是否優(yōu)良的基本依據(jù)。而品種區(qū)域試驗對參試品種進行豐產(chǎn)性、穩(wěn)產(chǎn)性鑒定結(jié)果則是該品種審定的主要依據(jù)[1]。

為更好地了解商豆1310的生產(chǎn)性能，使其在生產(chǎn)上能因地制宜合理利用。運用變異系數(shù)法、回歸系數(shù)法、高穩(wěn)系數(shù)法及AMMI模型對2015—2016年國家黃淮海夏大豆品種南組區(qū)域試驗的數(shù)據(jù)進行分析，進一步研究商豆1310的豐產(chǎn)性、穩(wěn)產(chǎn)性和適應(yīng)性，以期為其審定和生產(chǎn)上的推廣應(yīng)用提供科學依據(jù)。

1 材料與方法

1.1 試驗材料

分析數(shù)據(jù)來源于2015—2016年國家黃淮海夏大豆品種區(qū)域試驗?zāi)掀珺組及2016年國家黃淮海夏大豆品種生產(chǎn)試驗。對照品種均為中黃13。2015年國家南片B組區(qū)域試驗參試品種有：商豆1310、周豆22號、菏豆29號、濟J12105、山寧17、濰豆8號、徐0212-3、徐9418-2、晉大78號、中黃13（CK）；2016年參試品種有：商豆1310、周豆22號、濰豆8號、菏豆29號、阜10-11、濉科8號、濉科20、菏豆28號、山寧18、徐0117-46、周豆25號、興豆3號、圣豆11號、中黃13（CK）。2年區(qū)域試驗地點相同，均設(shè)置12個試點，分別是安徽阜陽、安徽龍亢、安徽宿州、河南商丘、河南周口、河南駐馬店、江蘇灌云、江蘇淮安、江蘇徐州、山東菏澤、山東濟寧、山東臨沂。2016年生產(chǎn)試驗設(shè)置7個試點，分別是：安徽龍亢、安徽宿州、河南周口、河南駐馬店、江蘇徐州、山東菏澤、山東濟寧。

1.2 統(tǒng)計分析方法

品種的豐產(chǎn)性分析是以各參試品種產(chǎn)量間進行差異顯著性分析，品種間的多重比較運用LSD法，并以較對照增產(chǎn)百分比和增產(chǎn)試驗點比例來評價該品種（系）豐產(chǎn)性能。

變異系數(shù)法參照Francis等1978年提出的變異系數(shù)和平均產(chǎn)量結(jié)合的基因型分組法。變異系數(shù)計算公式為：

高穩(wěn)系數(shù)法由溫振民等提出[4]，計算公式為：

公式中HSCi為第i個參試品種的高穩(wěn)系數(shù)，i，SDi分別為第i個品種的平均產(chǎn)量和標準差，ck為對照品種的平均產(chǎn)量，HSCi值越大，表明該品種高產(chǎn)穩(wěn)產(chǎn)性越好。

回歸系數(shù)法根據(jù)作物的產(chǎn)量是由基因和環(huán)境共同作用結(jié)果的原理，以參試品種在各試點的平均產(chǎn)量作為因變量，以各試點所有參試品種的平均產(chǎn)量為自變量進行回歸分析，以回歸系數(shù)的大小來衡量品種的適應(yīng)性[5-6]。當回歸系數(shù)b=1時，表示該品種具有廣泛的適應(yīng)性；當b＞1時，說明該品種對環(huán)境變化的敏感性增大，其適應(yīng)性低于平均適應(yīng)性；當b＜1時，表示該品種對環(huán)境變化不敏感[7-8]，其適應(yīng)性高于平均適應(yīng)性。以b值為橫軸，產(chǎn)量為縱軸，以所有參試品種的回歸系數(shù)均值和產(chǎn)量均值為坐標原點將坐標圖劃成4個象限，得出散點圖。落在第Ⅰ象限的為高產(chǎn)適應(yīng)性差品種，落在第Ⅱ象限的則是高產(chǎn)適應(yīng)性好品種，落在第Ⅲ象限的是低產(chǎn)適應(yīng)性好品種，落在第Ⅳ象限的品種低產(chǎn)適應(yīng)性差品種。理想的品種應(yīng)該是回歸系數(shù)b接近于1，平均產(chǎn)量較高。

AMMI模型是將方差分析和主成分分析有機結(jié)合在一起的分析方法?；贏MMI模型，張澤等提出以一個品種（或基因型）在交互效應(yīng)主成分軸（IPCA）空間中與原點的歐氏距離（記為D）作為評價該品種穩(wěn)定性的指標[9]，Di值越小則品種越穩(wěn)定。Di值計算，采用吳為人改進的公式進行[10]：

式中，i表示第i個品種，c為統(tǒng)計測驗顯著的IPCA個數(shù)為第i個品種在第s個IPCA上得分的平方。ws為第s個IPCA的權(quán)重，以該IPCA所解釋的變異（平方和）占全部IPCA所解釋的變異的比例來表示。

所有數(shù)據(jù)處理采用Excel軟件和DPS數(shù)據(jù)處理系統(tǒng)進行。

2 結(jié)果與分析

2.1 商豆1310的豐產(chǎn)性分析

2015年商豆1310參加國家黃淮海夏大豆品種區(qū)域試驗，平均產(chǎn)量3 574.05 kg/hm2，較對照中黃13增產(chǎn)10.07%，較對照增產(chǎn)點次為91.70%，位居該組試驗第1位，增產(chǎn)極顯著，而且與產(chǎn)量位次為第2位的周豆22號產(chǎn)量間差異也達極顯著水平；2016年參加區(qū)域試驗，平均產(chǎn)量3 105.55 kg/hm2，較對照中黃13增產(chǎn)12.51%，較對照增產(chǎn)點次為100%，位居該組試驗第1位，增產(chǎn)極顯著水平（見表1）；2016年參加國家黃淮海夏大豆品種生產(chǎn)試驗，平均產(chǎn)量達到3 010.30 kg/hm2，較對照中黃13增產(chǎn)11.22%，居生產(chǎn)試驗第1位，各試驗點全部增產(chǎn)，增產(chǎn)幅度為2.72%～24.64%，其中產(chǎn)量最高試點達到3 599.85 kg/hm2（見表2）。經(jīng)過多年多點的試驗及示范試種，商豆1310產(chǎn)量表現(xiàn)穩(wěn)定，綜合性狀突出，具有很好的豐產(chǎn)性，是一個極具增產(chǎn)潛力的大豆新品種。

2.2 商豆1310不同年度、不同地點間豐產(chǎn)性、穩(wěn)產(chǎn)性及適應(yīng)性分析

對2015—2016年度國家黃淮海夏大豆品種區(qū)域試驗?zāi)掀珺組產(chǎn)量結(jié)果進行方差分析（見表3）。結(jié)果表明，2015年各參試品種（系）的品種主效、試點主效及品種×試點效應(yīng)均達到極顯著水平，商豆1310較對照增產(chǎn)顯著，居該組試驗第1位；2016年參試品種（系）的品種主效、試點主效及品種×試點效應(yīng)達極顯著水平，商豆1310與對照品種差異達極顯著。說明商豆1310在不同年度間、不同地點間及各試點內(nèi)小區(qū)間重復性較好，應(yīng)對不同氣候環(huán)境時的適應(yīng)性較強、豐產(chǎn)性好、產(chǎn)量表現(xiàn)穩(wěn)定。

2.3 商豆1310的高產(chǎn)穩(wěn)產(chǎn)性綜合分析

品種的穩(wěn)產(chǎn)性分析以變異系數(shù)來度量品種的穩(wěn)定性，變異系數(shù)越大，說明該品種穩(wěn)定性差。高穩(wěn)系數(shù)綜合了豐產(chǎn)性和穩(wěn)產(chǎn)性因素綜合評價品種，高穩(wěn)系數(shù)值越大，表明品種的高產(chǎn)穩(wěn)產(chǎn)性越好，高穩(wěn)系數(shù)值越小，則品種的高產(chǎn)穩(wěn)產(chǎn)性越差。對2015和2016年國家黃淮海夏大豆品種區(qū)域試驗?zāi)掀珺組各參試品種（系）產(chǎn)量的變異系數(shù)（CV）和高穩(wěn)系數(shù)（HSC）進行分析（見表1）。結(jié)果表明，2015年商豆1310的CV值和HSC值分別為5.97%和94.09%，低于對照中黃13的CV值（6.37%），較對照中黃13的HSC值（85.12%）大，分別位居第2、第1位；2016年商豆1310的CV值和HSC值分別為5.54%和96.61%，均位居第1位，且標準差位居第1位。運用Francis等提出的將變異系數(shù)和平均產(chǎn)量結(jié)合的基因型分組法，得到2年各參試品種（系）的散點圖（見圖1）。由圖1可知，商豆1310均落在第Ⅱ象限，為高產(chǎn)、穩(wěn)產(chǎn)品系。以上均說明商豆1310是一個豐產(chǎn)性和穩(wěn)產(chǎn)性均好的品種。

2.4 商豆1310的適應(yīng)性分析

2015年區(qū)域試驗，商豆1310的回歸系數(shù)（b）為0.998 0；2016年區(qū)域試驗，商豆1310的回歸系數(shù)為0.942 0。兩年的b值均小于1，且僅接近于1，說明商豆1310對環(huán)境變化不敏感，適應(yīng)性高于平均適應(yīng)性，且具有廣泛的適應(yīng)性。由參試品種（系）產(chǎn)量和回歸系數(shù)作出2年的散點圖（見圖2）可知，商豆1310均落在第Ⅱ象限，是高產(chǎn)適應(yīng)性好品種。加之該品種的平均產(chǎn)量明顯高于對照。綜上所述，商豆1310是一個高產(chǎn)而又適應(yīng)性廣的品種。

2.5 AMMI模型對商豆1310的穩(wěn)產(chǎn)性分析

AMMI模型分析表明（見表4），AMMI模型能較好的解釋各種變異。2015年基因、環(huán)境、基因與環(huán)境交互作用分別能解釋5.55%、81.46%、12.99%的平均產(chǎn)量變異，均達到極顯著水平。IPCA1～IPCA3均達極顯著水平，IPCA1、IPCA2、IPCA3分別解釋基因與環(huán)境交互作用的29.89%、25.82%、23.25%，三者累計解釋互作平方和的78.96%，占大部分；2016年基因、環(huán)境、基因與環(huán)境交互作用分別能解釋6.28%、67.84%、25.88%的平均產(chǎn)量變異，均達到極顯著水平。IPCA1～IPCA3均達極顯著水平，IPCA1、IPCA2、IPCA3分別解釋基因與環(huán)境交互作用的56.49%、14.66%、10.28%；三者累計解釋互作平方和的81.43%，占大部分。這說明AMMI模型能夠有效地分析基因與環(huán)境間互作。

根據(jù)公式（3）計算各參試品種（系）對應(yīng)的Di值（見表5）。Di值越小，離原點越近，說明品種穩(wěn)定性越好。2年結(jié)果表明，商豆1310的Di值，均位于第4位，對照中黃13的Di值均位于第3位，雖然商豆1310的Di值比對照小，但比較接近，且排名靠前，這能夠說明商豆1310的穩(wěn)產(chǎn)性較好。

表1 商豆1310 與其他參試品種（系）的產(chǎn)量、變異系數(shù)、高穩(wěn)系數(shù)和回歸系數(shù)

表2 商豆1310在國家黃淮海夏大豆品種生產(chǎn)試驗中的產(chǎn)量表現(xiàn)

表3 商豆1310不同年度、不同地點間方差分析結(jié)果及豐產(chǎn)性、穩(wěn)產(chǎn)性分析

圖1 參試品種（系）產(chǎn)量與變異系數(shù)散點圖

圖2 參試品種（系）產(chǎn)量與回歸系數(shù)散點圖

圖3 參試品種（系）產(chǎn)量與歐氏距離值散點圖

表4 AMMI模型對產(chǎn)量的分析

以各參試品種（系）產(chǎn)量為縱坐標，Di值為橫坐標，以所有參試品種的產(chǎn)量均值和Di均值為坐標原點作散點圖（見圖3）。由圖3可明顯看出，各參試品種（系）的產(chǎn)量類型。對照中黃13落在第Ⅲ象限為低產(chǎn)穩(wěn)產(chǎn)品種，商豆1310落在第Ⅱ象限為高產(chǎn)穩(wěn)產(chǎn)品種。

表5 參試品種（系）在顯著交互效應(yīng)軸上的得分、Di值及排序

3 討論與結(jié)論

運用與對照品種進行比較的方法分析品種的豐產(chǎn)性，簡單、直接、清晰明了，但不能反映品種的穩(wěn)產(chǎn)性[11]，例如單看2015年區(qū)域試驗表現(xiàn)，周豆22號產(chǎn)量居試驗第2位，但由于其標準差和變異系數(shù)較大，可以說明該品種是一個高產(chǎn)，但不穩(wěn)產(chǎn)的品種。

變異系數(shù)只是品種表型穩(wěn)定性的近似測度，主要反映了品種本身能對產(chǎn)量的影響，而對環(huán)境因素影響產(chǎn)量的程度體現(xiàn)不足，且很難反映出基因與環(huán)境間互作，變異系數(shù)可分析品種穩(wěn)產(chǎn)性，但體現(xiàn)不出品種的豐產(chǎn)性[12-13]。例如徐9418-2，標準差居第3位，變異系數(shù)居第4位，說明該品種穩(wěn)產(chǎn)性較好，但該品種較對照中黃13減產(chǎn)，不能體現(xiàn)品種的豐產(chǎn)性。但是變異系數(shù)法計算簡便，結(jié)合平均產(chǎn)量的散點圖可對參試品種（系）的產(chǎn)量潛力和穩(wěn)定性進行簡單的類型劃分，育種者可直觀的篩選出高產(chǎn)穩(wěn)產(chǎn)性的品種，此法適用于育種初期大量早代材料的多點無重復試驗[14]。

高溫系數(shù)法可以把目標品種的高產(chǎn)性和穩(wěn)產(chǎn)性結(jié)合起來綜合分析，通過參數(shù)HSCi的大小來衡量參試品種高產(chǎn)穩(wěn)產(chǎn)性的優(yōu)劣，計算方法簡便，結(jié)果量化便于排序，在產(chǎn)量試驗中具有較好的參考價值[12]，但是HSCi值高的品系存在標準差、變異系數(shù)大的現(xiàn)象，例如2015年區(qū)域試驗中，周豆22的標準差和變異系數(shù)分別位于第8，6位，但由于其產(chǎn)量高，居產(chǎn)量位第2位，通過計算其HSCi值位于第3位；2016年區(qū)域試驗中，濉科8號的標準差和變異系數(shù)雖較大，均位于第11位，因其產(chǎn)量較高位于第3位，使其HSCi值位于第5位。這說明了參試品系的穩(wěn)定性劣勢易被高產(chǎn)掩蓋。高穩(wěn)系數(shù)法還有不足之處是其忽視了品種產(chǎn)量高低受基因型和環(huán)境互作的影響，所以在實際運用中應(yīng)配合穩(wěn)定性參數(shù)進行分析效果更好[12,15]。

回歸系數(shù)法體現(xiàn)了環(huán)境對品種產(chǎn)量的影響，只能反映品種的適應(yīng)性。

AMMI模型綜合了方差分析和主成分分析的優(yōu)點，并可以性狀參數(shù)為指標對品種穩(wěn)定性進行評價，被廣泛應(yīng)用于各種動植物試驗分析中[16]，在區(qū)域試驗分析中具有較大應(yīng)用價值。AMMI模型不僅能有效揭示環(huán)境、基因及基因與環(huán)境互作效應(yīng)對品種穩(wěn)定性的影響，基于AMMI模型的歐氏距離值還可定量地分析品種的穩(wěn)產(chǎn)性表現(xiàn)，此法在分析時考慮了諸多對品種穩(wěn)定性具有影響的因素，分析過程嚴謹，因而在區(qū)域試驗品種穩(wěn)產(chǎn)性數(shù)據(jù)分析中更為可靠[12,17-18]。但是此方法在分析過程中對科技人員數(shù)據(jù)分析能力要求較高，基層農(nóng)技人員掌握起來有一定難度。采用AMMI模型分析參試品系的高產(chǎn)穩(wěn)產(chǎn)性時，借鑒Francis等的變異系數(shù)法將Di值與平均產(chǎn)量結(jié)合得到了各參試品系的散點圖，通過散點圖實現(xiàn)對參試品種的高產(chǎn)穩(wěn)產(chǎn)性綜合評價，如此一來不僅是試驗結(jié)果的可靠性得到保證，評價結(jié)果也較直觀更易于理解[12]。

綜上所述，在具體的分析過程中，不能局限于固定的方法和模型，應(yīng)結(jié)合育種目標和技術(shù)能力選擇合適的方法，并綜合參試品系在試驗中的具體表現(xiàn)，盡可能多的收集多年多點數(shù)據(jù)進行分析。

本試驗運用多種方法進行分析，結(jié)果表明：商豆1310具有較高的豐產(chǎn)性和穩(wěn)定性及較廣適應(yīng)性，是一個極具推廣潛力的大豆新品種，適合黃淮海南部地區(qū)大面積推廣種植。

[1] 蘇秋芹.花生新品種龍花163豐產(chǎn)性和穩(wěn)產(chǎn)性分析[J].中國農(nóng)學通報,2008,25(24):191-195.

[2] 穆培源,莊麗,張吉貞,等.作物品種穩(wěn)定性分析方法的研究進展[J].新疆農(nóng)業(yè)科學,2003,40(3):142-144．

[3] 王兵偉,時成俏,覃德斌,等.玉米品種豐產(chǎn)性和穩(wěn)定性分析方法比較[J].種子,2009,28(1):105-108.

[4] 溫振民,張永科.用高穩(wěn)系數(shù)法估算玉米雜交種高產(chǎn)穩(wěn)產(chǎn)性的探討[J].作物學報,1994,20(4):508-512.

[5] 劉正學,劉飛,李寶強,等.小麥新品種臨麥4號豐產(chǎn)性穩(wěn)產(chǎn)性及適應(yīng)性分析[J].中國農(nóng)學通報,2008,24(12):225-227.

[6] 李筠,王龍,晁林海,等.國審小麥連麥2號豐產(chǎn)性、穩(wěn)產(chǎn)性及適應(yīng)性分析[J].安徽農(nóng)業(yè)科學,2007,35(29):9191-9193.

[7]馮宜梅.玉米新品種武科8號的豐產(chǎn)穩(wěn)產(chǎn)性及適應(yīng)性分析[J].中國種業(yè),2015(7):30-32.

[8] 孫化軍,閆延梅,閆向前,等.大豆新品種商豆14號豐產(chǎn)性穩(wěn)產(chǎn)性及適應(yīng)性分析[J].農(nóng)業(yè)科技通訊,2014(9):104-107.

[9] 張澤,魯成,向仲懷.基于AMMI模型的家蠶品種穩(wěn)定性分析[J].中國農(nóng)業(yè)科學,1998,31(1):62-68.

[10] 吳為人.對基于AMMI模型的品種穩(wěn)定性分析方法的一點改進[J].遺傳,2000,22(1):31-32.

[11] 舒文濤,耿臻,楊青春,等.國審大豆新品種‘周豆23號’豐產(chǎn)性、穩(wěn)產(chǎn)性、適應(yīng)性分析[J].農(nóng)學學報,2016,6(12):8-10.

[12] 劉帆,袁玉順,楊俊青,等.大麥新品系高產(chǎn)穩(wěn)產(chǎn)性評價方法比較[J].云南農(nóng)業(yè)大學學報:自然科學,2016,31(4):597-603.

[13] 周青,范陽,徐淑霞,等.綜合分析大豆區(qū)試中品種的高產(chǎn)穩(wěn)產(chǎn)和適應(yīng)性[J].陜西農(nóng)業(yè)科學,2008(4):14-15.

[14] 穆培源,莊麗,張吉貞,等.作物品種穩(wěn)定性分析方法的研究進展[J].新疆農(nóng)業(yè)科學,2003,40(3):142-144.

[15] 劉占柱,姚丹,沈剛,等.幾種穩(wěn)定性分析法在大豆品種區(qū)試中的應(yīng)用與評析[J].吉林農(nóng)業(yè)科,2007,32(2):19-21,29.

[16] 劉博,衛(wèi)玲,樊云茜,等.基于AMMI模型的黃淮海夏大豆國家區(qū)試產(chǎn)量分析[J].中國農(nóng)學通報,2015,31(27):69-74.

[17] 陳勇生,鄧海華,劉福業(yè),等.對應(yīng)用幾種統(tǒng)計模型評價甘蔗品種穩(wěn)定性的初步比較[J].生物數(shù)學學報,2012,27(1):168-174.

[18] 何靄如,李觀康,陳勝勇,等.用AMMI模型分析甘薯品種產(chǎn)量性狀的穩(wěn)定性[J].安徽農(nóng)業(yè)大學學報,2014,41(3):430-434.

Analysis on High-yield Capacity,Stability and Adaptability of a New Soybean Variety Shangdou 1310

He Xin,Yan Xiangqian,Sun Huajun,Zhang Qi
(Henan Shangqiu Academy of Agriculture and Forestry Sciences,Shangqiu 476000,Henan,China)

The study aims to better understand the production characteristic of a new soybean variety Shangdou 1310.Based on the regional test results in Southern Huang-huai-hai region from 2015 and 2016,the yielding capacity,stability and adaptability of Shangdou 1310 were analyzed by variation coefficient method,high and stable coefficient(HSC),regression coefficient method,and AMMI model as well.The results showed that variety Shangdou 1310 had high yield,good stability and adaptability,and it can be widely planted for production in the Southern Huang-Huai region.

Soybean variety;Shangdou 1310;High yield;Stability;Adaptability

S565.1

A

1674-3547(2017)05-0037-09

2017-06-15

何鑫，女，碩士，主要從事大豆育種與栽培技術(shù)研究，E-mail:531863729@qq.com

現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項（CARS—004）