摘 要：為篩選適宜河南地區(qū)種植的農(nóng)藝性狀優(yōu)良的朝天椒新品種，并為辣椒品種改良提供理論依據(jù)，對38個朝天椒品種的果實性狀進(jìn)行差異性分析、相關(guān)性分析、主成分分析和聚類分析。結(jié)果表明，品種間果實性狀差異較大，變異系數(shù)范圍為9.24%～53.03%，其中二氫辣椒素含量變異程度最大；果實性狀間存在相關(guān)性，其中果實長度和果實橫徑與單果質(zhì)量、果實長度與心室數(shù)、果皮厚度與烘干所需時間、辣椒素含量與二氫辣椒素含量均呈極顯著正相關(guān)。主成分分析結(jié)果表明，共篩選出4個主成分，方差累積貢獻(xiàn)率為70.287%，能反映辣椒農(nóng)藝性狀的大部分信息。通過表型聚類和液相芯片基因型數(shù)據(jù)聚類分析，明確了供試材料間的親緣關(guān)系。綜上，品種DKP-8418（P31）、好農(nóng)515紅綠（P28）和DO1844歐蘭德（P24）綜合表現(xiàn)優(yōu)良，適宜在河南地區(qū)推廣種植。

關(guān)鍵詞：朝天椒；農(nóng)藝性狀；主成分分析；聚類分析

中圖分類號：S641.3 文獻(xiàn)標(biāo)志碼：A 文章編號：1673-2871（2024）08-067-09

Fruit traits evaluation and cluster analysis of different Chaotian pepper varieties

Lü Mengyao1， WANG Panqiao1， LI Jisuo2， LIU Keke1， LI Ying1， LI Chen1， JIA Zhiqi1

（1. College of Horticulture， Henan Agricultural University， Zhengzhou 450046， Henan， China; 2. Henan Qidu Agricultural Technology Co.， Ltd， Zhengzhou 450000， Henan， China）

Abstract： In order to screen new varieties of Chaotian pepper with good agronomic characters suitable for planting in Henan， and provide theoretical basis for the improvement of pepper varieties， the difference analysis， correlation analysis， principal component analysis and cluster analysis of fruit characters of 38 Chaotian pepper varieties were carried out. The results showed that the variation coefficient of fruit traits was 9.24%-53.03%， and the variation degree of dihydrocapsaicin character was the largest. There was significant positive correlation between fruit length， fruit transverse diameter and single fruit mass， fruit length and ventricular number， peel thickness and time of drying， capsaicin content and dihydrocapsaicin content. The results of principal component analysis showed that four principal components were selected， and the cumulative contribution rate of variance was 70.287%， which could reflect most of the information of capsicum agronomic traits. Through phenotypic clustering and liquid chip genotypic data clustering analysis， the relationship between the tested materials was clarified. In summary， the varieties DKP-8418（P31）， Haonong 515 Red Green（P28）and DO1844 Olander（P24）showed excellent performance and were suitable for popularization and planting in Henan.

Key words： Pepper; Agronomic traits; Principal component analysis; Clustering analysis

朝天椒（Capsicum annuum L. var. conoides）是世界上重要的蔬菜作物之一，屬于茄科辣椒屬一年或多年生植物，是辣椒的一個重要亞種。我國朝天椒栽培歷史悠久，產(chǎn)量和種植面積均居世界第一位，總產(chǎn)量占世界總產(chǎn)量的70%以上[1]。由于朝天椒口感獨特、營養(yǎng)豐富，所以常被用于制作各種調(diào)味品，是制作辣椒面、辣椒醬等調(diào)味品的重要原料。朝天椒富含辣椒素，是辣椒中一種重要的生理活性物質(zhì)，具有抗氧化、抗炎、抗菌、抗病毒、調(diào)節(jié)免疫等生理功能，在醫(yī)藥和食品工業(yè)中有著廣泛的應(yīng)用[2]。

近年來，河南省朝天椒生產(chǎn)規(guī)?；匾呀?jīng)形成，在柘城、臨潁、內(nèi)黃等主產(chǎn)區(qū)，朝天椒已成為當(dāng)?shù)氐闹鲗?dǎo)產(chǎn)業(yè)。河南地區(qū)目前種植的朝天椒品種單一、辣味低，不能滿足市場需求，市場急需綜合性狀優(yōu)良的朝天椒品種[3]。朝天椒果實的脫水速度影響朝天椒的收獲采摘和機械化采收，是衡量朝天椒果實特性的重要指標(biāo)，但對該性狀的評價研究較少。另外，目前朝天椒育種以傳統(tǒng)常規(guī)方法為主，分子育種水平較低。基于SNP的靶向測序基因型檢測（genotyping by target sequencing，GBTS）技術(shù)，也稱為液相芯片技術(shù)，因為其高通量和低成本，在各個農(nóng)作物育種中的應(yīng)用越來越廣泛[4]，辣椒液相芯片的開發(fā)使用，有助于深入了解辣椒品種間的親緣關(guān)系，提高育種效率，對提高我國辣椒生物育種水平具有重要意義。

王丹丹等[5]對辣椒主要農(nóng)藝性狀相關(guān)性的研究表明，單果質(zhì)量與商品果橫徑呈極顯著正相關(guān)。張國儒等[6]將272份辣椒資源進(jìn)行聚類分析，發(fā)現(xiàn)在歐式距離7.5處將全部種質(zhì)資源分為6類。楊志剛等[7]采用聚類分析法對92份辣椒資源的色價值及性狀指標(biāo)進(jìn)行分析，篩選出了用于干制紅辣椒的優(yōu)良種質(zhì)材料17份。盡管對辣椒農(nóng)藝性狀已有相關(guān)方面的研究，但隨著辣椒種質(zhì)資源的不斷創(chuàng)新和辣椒遺傳育種的需要以及河南辣椒新種質(zhì)的不斷引進(jìn)，仍需對辣椒種質(zhì)資源進(jìn)行評價研究。筆者對來源于22家公司及科研單位的38個朝天椒品種進(jìn)行綜合性狀評價，對辣椒素含量、單果質(zhì)量、果實長度、果實橫徑、果肉厚度等主要果實性狀進(jìn)行相關(guān)性分析、主成分分析以及液相芯片基因型聚類分析，以期篩選出適合河南地區(qū)露地栽培的優(yōu)良朝天椒品種。

1 材料與方法

1.1 材料

試驗材料來自22家育種單位提供的38個朝天椒新品種，均為雜交1代，品種來源信息見表1。

1.2 試驗設(shè)計

試驗材料于2022年種植于在河南省鄭州市惠濟區(qū)中原地展試驗基地，露地栽培，5月1日用72孔穴盤播種育苗，6月5日定植，每份材料20株，株行距30 cm×50 cm，3次重復(fù)，完全隨機分布。于11月12日上午進(jìn)行果實采集。

1.3 果實性狀及脫水速度評價

1.3.1 朝天椒果實性狀評價 從每個品種測試果中取10個鮮果，對測試品種的單果質(zhì)量、心室數(shù)、種子數(shù)量進(jìn)行觀察、測量、統(tǒng)計，用游標(biāo)卡尺測量果實的長度，使用神眇電子光學(xué)數(shù)字顯微鏡測量果肉厚度、果實橫徑[8]。參考麥玉連[9]的方法測定辣椒素和二氫辣椒素含量。

1.3.2 果實脫水速度評價 每個品種隨機取10個果實，先稱取鮮椒的總質(zhì)量，再放入烘箱60 ℃ 烘烤24 h，轉(zhuǎn)入45 ℃后每24小時稱質(zhì)量1次，直到3次稱質(zhì)量穩(wěn)定，作為果實干質(zhì)量，并記錄所需時間d（天），即烘干所需時間。

1.4 數(shù)據(jù)處理

采用Excel對數(shù)據(jù)進(jìn)行統(tǒng)計與整理；采用 SPSS 進(jìn)行相關(guān)性分析、主成分分析和聚類分析；利用與中國農(nóng)業(yè)大學(xué)園藝學(xué)院、河南七度農(nóng)業(yè)科技有限公司和石家莊博瑞迪生物技術(shù)有限公司合作開發(fā)的辣椒40 K液相芯片進(jìn)行基因型分析和聚類分析。

2 結(jié)果與分析

2.1 各品種果實性狀的變異系數(shù)分析

變異系數(shù)反映變量之間的相對離散程度，對38個試驗品種的其中10個指標(biāo)進(jìn)行了變異系數(shù)分析，從表2可以看出，在10個指標(biāo)中，心室數(shù)的變異系數(shù)最小，二氫辣椒素含量的變異系數(shù)最大。變異系數(shù)從小到大依次為：心室數(shù)（9.24%）＜果實橫徑（12.90%）＜干鮮比（16.18%）＜果實長度（18.58%）＜烘干所需時間（21.20%）＜單果質(zhì)量（25.26%）＜果肉厚度（26.32%）＜種子數(shù)量（27.00%）＜辣椒素含量（48.08%）＜二氫辣椒素含量（53.03%），38個朝天椒品種間性狀差異較大，具有較為豐富的遺傳多樣性，適宜篩選優(yōu)異朝天椒品種。

2.2 各品種果實性狀的相關(guān)性分析

對38個試驗品種的其中10個指標(biāo)進(jìn)行相關(guān)性分析，結(jié)果由表3所示，果實長度、果實橫徑與單果質(zhì)量呈極顯著正相關(guān)，果肉厚度、烘干所需時間與單果質(zhì)量呈顯著正相關(guān)，果肉厚度與烘干所需時間呈極顯著正相關(guān)，果實長度與心室數(shù)呈極顯著正相關(guān)；辣椒素含量與二氫辣椒素含量呈極顯著正相關(guān)。

2.3 主成分分析及綜合評價

為充分反映各因素中起主導(dǎo)作用的綜合指標(biāo)，進(jìn)一步明確農(nóng)藝性狀對38個品種的影響，對9個農(nóng)藝性狀進(jìn)行主成分分析，結(jié)果見表4，共篩選出特征值大于1的4個主成分，方差貢獻(xiàn)率分別為25.473%、20.484%、12.680%、11.649%，方差累積貢獻(xiàn)率達(dá)到了70.287%，能綜合反映辣椒農(nóng)藝性狀的大部分信息。其中，第一主成分的特征值是2.293，以單果質(zhì)量的特征向量值最大，其次為果實橫徑、果實長度和心室數(shù)，是果實性狀因子；第二主成分的主要影響因子為辣椒素含量和二氫辣椒素含量，是辣度因子；在第三主成分中，果肉厚度的特征向量值最大，是果肉厚度因子；第四主成分種子數(shù)量的特征向量值最大，是單果種子數(shù)量因子。

根據(jù)因子得分系數(shù)矩陣及其所對應(yīng)的主成分，可以計算出4個主成分因子得分。Fl=0.525X1+0.426X2+0.497X3+0.191X4+0.416X5+0.226X6+0.080X7-0.070X8-0.169X9；F2=0.089X1+0.205X2-

0.018X3+0.057X4+0.105X5-0.168X6+0.191X7+0.663X8+0.656X9；F3=0.184X1-0.364X2+0.013X3+0.737X4-0.221X5+0.219X6-0.374X7+0.222X8+0.002X9；

F4=0.068X1-0.320X2+0.152X3-0.220X4-0.269X5+0.685X6+0.524X7+0.086X8+0.091X9。

其中，X1～X9分別表示單果質(zhì)量、果實長度、果實橫徑、果肉厚度、心室數(shù)、種子數(shù)量、干鮮比、二氫辣椒素含量和辣椒素含量9個性狀指標(biāo)。

以各主成分特征值比率為權(quán)重，構(gòu)建38個品種性狀綜合評價得分函數(shù)。F=0.362F1+0.292F2+0.180F3+0.165F4。其中，F(xiàn)為38個辣椒品種9個性狀的綜合評價得分。不同辣椒品種性狀的綜合得分及排名如表5所示，F(xiàn)值越大，表明該品種的綜合品質(zhì)性狀越好。其中，編號P31（DKP-8418）綜合得分122.04，排名第1；編號P28（好農(nóng)515紅綠）綜合得分121.04，排名第2；編號P24（DO1844歐蘭德）綜合得分110.49，排名第 3。

2.4 基于果實性狀的聚類分析

在主成分分析的基礎(chǔ)上，用SPSS軟件采用平方歐氏距離法對9個農(nóng)藝性狀進(jìn)行聚類分析，在歐氏距離為20處，將38個品種分為3類，分別包含22、14和2個品種（圖1）。由圖1和表6可以看出，第I類群共有22個品種，占試驗品種總數(shù)的57.9%，此類群二氫辣椒素含量最高，辣椒素含量次之，除干鮮比處，其他農(nóng)藝性狀值均最低，主要有P1、P3、P4、P5、P6等；第II類群共有14個品種，占試驗品種總數(shù)的36.8%，此類群心室數(shù)和果肉厚度均為最高，主要有P2、P7、P9、P15等；第Ⅲ類群共有2個品種，占試驗品種總數(shù)的5.3%，此類群單果質(zhì)量、果實長度、果實橫徑、心室數(shù)、種子數(shù)量、干鮮比和辣椒素含量均最大或最高，主要有P8和P28。

2.5 基于液相芯片的基因型聚類分析

筆者利用液相芯片對38個朝天椒全基因組范圍內(nèi)的SNP進(jìn)行檢測，結(jié)果得到最小等位基因頻率（MAF）為0.013 16，群體觀察雜合度（Ho）最大值為1.000 0，最小值為0.026 3；期望雜合度（He）最大值為0.500 00，最小值0.025 97。選取的38個朝天椒品種的表型數(shù)據(jù)大部分品種間存在一定的遺傳差異。為了探索該芯片所包含的基因型是否能夠應(yīng)用于親緣關(guān)系的評價，對38個朝天椒品種進(jìn)行了聚類分析，結(jié)果如圖2所示，在歐式距離為 2 時，將38個品種分為3類，第I類群，共有7個品種，占試驗品種總數(shù)的18.4%，均來自不同單位；第II類群共有18份種質(zhì)，占試驗品種總數(shù)的47.4%，其中有 4 個品種來源于河南農(nóng)業(yè)大學(xué)園藝學(xué)院，有 3 個品種來源于昆明市華興種業(yè)有限公司，京研益農(nóng)（北京）種業(yè)科技有限公司和湖南湘研種業(yè)有限公司分別有2個品種，其他均來自不同公司；第Ⅲ類群共有13個品種，占試驗品種總數(shù)的34.2%，其中有 3 個品種來源于河南紅綠辣椒種業(yè)有限公司，其他均來自不同公司。部分公司存在交叉現(xiàn)象，同時被歸到不同的類群中，例如京研益農(nóng)（北京）種業(yè)科技有限公司的4個品種被歸到I、II和Ⅲ類群中；河南農(nóng)業(yè)大學(xué)園藝學(xué)院和北京大一韓日國際種業(yè)有限公司的品種被歸到I和II類群中；歐蘭德種業(yè)有限公司和內(nèi)蒙古天紅椒業(yè)股有限公司的品種被歸到I和Ⅲ類群中；沈陽市辣伙伴農(nóng)業(yè)有限公司的品種被歸到II和Ⅲ類群中。

3 討論與結(jié)論

農(nóng)藝性狀的觀測調(diào)查是辣椒新品種培育和篩選的重要手段，通過變異度與多樣性指數(shù)分析，能夠了解不同品種性狀的差異程度[10]。韓婭楠等[11]研究表明，24個農(nóng)藝性狀中株幅、果梗長度、果肉厚度的多樣性指數(shù)較大，性狀遺傳較豐富；8個品質(zhì)性狀中，辣椒素和二氫辣椒素含量的變異系數(shù)較大，分別為106.98%和112.44%。在本研究的10個農(nóng)藝指標(biāo)中，變異系數(shù)最大的是二氫辣椒素含量，其次是辣椒素含量，與韓婭楠等[11]、李寧等[12]、詹永發(fā)等[13]的研究結(jié)果一致。說明這38個品種農(nóng)藝性狀的變異潛力較大，有助于篩選出優(yōu)異的朝天椒品種。

相關(guān)性分析能反映不同性狀間的關(guān)系，為制定育種目標(biāo)和新品種選育提供參考。袁娟偉等[14]對275份辣椒種質(zhì)資源表型性狀進(jìn)行相關(guān)性分析，結(jié)果表明，單果質(zhì)量與果實橫徑、果肉厚度、心室數(shù)呈顯著正相關(guān)。筆者在本研究中的結(jié)果表明，單果質(zhì)量與果實長度、果實橫徑和果肉厚度呈極顯著正相關(guān)，與袁娟偉等[14]的研究結(jié)果一致。同時本研究中果肉厚度與烘干所需時間呈極顯著正相關(guān)，果實長度與心室數(shù)呈極顯著正相關(guān)。很多作物的農(nóng)藝性狀間存在相關(guān)性，在進(jìn)行育種工作時可以利用性狀間的相關(guān)性對一些性狀進(jìn)行選擇。因此，在今后辣椒的高產(chǎn)育種中，可以將果實形狀、果肉厚度、單果質(zhì)量等作為主要的參考性狀。

主成分分析采取降維的方法，利用幾個綜合因子代表原來眾多的變量，使這些綜合因子盡可能地反映原來變量的信息量[15]。王春勇等[16]采用主成分分析、聚類分析等方法對17個鮮加工型辣椒品種農(nóng)藝性狀進(jìn)行分析，結(jié)果表明，前4個主成分特征值均大于1，累積貢獻(xiàn)率達(dá)80.748%，主成分綜合得分由高到低是嬌鳳、龍椒17號、園椒35等。王丹丹等[5]對14個辣椒品種的11個農(nóng)藝性狀進(jìn)行主成分分析，結(jié)果表明，前3個主成分的累積貢獻(xiàn)率達(dá)84.865%，在歐氏距離為10時，將其分為5大類。筆者對38個朝天椒品種的9個主要農(nóng)藝性狀進(jìn)行主成分分析，結(jié)果表明，特征值大于1的主成分有4個，累積貢獻(xiàn)率達(dá)70.287%，第一主成分以單果質(zhì)量的特征向量值最大；第二主成分主要影響因子是辣椒素含量和二氫辣椒素含量，第三主成分中，果肉厚度的特征向量值最大，第四主成分中種子數(shù)量的向量值最大；主成分綜合得分前三名由高到低依次是P31、P28、P24。

聚類分析常用于類群間的相互關(guān)系研究，特別是作物品種資源的分類與遺傳多樣性分析。任朝輝等[17]、鞏雪峰等[18]通過聚類分析將多份種質(zhì)材料進(jìn)行分類，并根據(jù)辣椒素含量篩選出優(yōu)質(zhì)類群。步洪鳳等[19]對26個辣椒品種的13個主要農(nóng)藝性狀進(jìn)行相關(guān)性和聚類分析，表明果形和產(chǎn)量是影響聚類結(jié)果的重要因素。在本研究中，基于果實性狀調(diào)查結(jié)果，將38個朝天椒品種分為3類。第Ⅰ類群二氫辣椒素含量最高；第Ⅱ類群心室數(shù)和果肉厚度最高；第Ⅲ類群單果質(zhì)量、果實長度、果實橫徑、心室數(shù)、種子數(shù)量、干鮮比和辣椒素含量最高。筆者在本研究中利用與中國農(nóng)業(yè)大學(xué)共同開發(fā)的辣椒40 K液相芯片，對38個朝天椒品種進(jìn)行基因型分析，用NTSYS軟件對朝天椒基因型數(shù)據(jù)進(jìn)行聚類分析，從聚類結(jié)果可以看出，各個品種在分子水平上有較好的多態(tài)性，說明各個公司的育種材料來源廣泛，種質(zhì)資源的遺傳多樣性較豐富，但個別單位也存在品種遺傳多樣性有限、遺傳基礎(chǔ)狹窄的問題。因此在辣椒育種中應(yīng)該注重選用外來資源來拓展我國栽培辣椒的遺傳基礎(chǔ)。另外，辣椒45 K液相芯片還可以用于朝天椒遺傳圖譜構(gòu)建、種質(zhì)資源鑒定、功能基因定位、品種定向選育、全基因組關(guān)聯(lián)分析等諸多方面，能夠大幅度提高朝天椒育種效率和分子育種水平[20]。

綜上所述，不同朝天椒品種間果實性狀差異較大，變異系數(shù)范圍為9.24%～53.03%，其中辣椒素性狀變異程度最大；果實性狀間存在相關(guān)性，其中果實長度、果實橫徑與單果質(zhì)量，果實長度與心室數(shù)，果皮厚度與烘干所需時間，辣椒素含量與二氫辣椒素含量均呈極顯著正相關(guān)。主成分分析結(jié)果表明，品種DKP-8418（P31）、好農(nóng)515紅綠（P28）和DO1844歐蘭德（P24）綜合表現(xiàn)優(yōu)良，適宜在河南地區(qū)推廣種植。

參考文獻(xiàn)

[1] 周悅，劉一彤，白麗，等．中國辣椒出口貿(mào)易結(jié)構(gòu)與國際競爭力分析[J]．中國蔬菜，2023（3）：1-7．

[2] 鄒學(xué)校，馬艷青，戴雄澤，等．辣椒在中國的傳播與產(chǎn)業(yè)發(fā)展[J]．園藝學(xué)報，2020，47（9）：1715-1726．

[3] 馬德俊，吳雪梅．露地朝天椒新品種引進(jìn)試驗簡報[J]．蔬菜，2022（6）：72-73．

[4] 徐云碧，楊泉女，鄭洪建，等．靶向測序基因型檢測（GBTS）技術(shù)及其應(yīng)用[J]．中國農(nóng)業(yè)科學(xué)，2020，53（15）：2983-3004．

[5] 王丹丹，師建華，李燕，等．基于主成分與聚類分析的辣椒主要農(nóng)藝性狀評價[J]．中國瓜菜，2021，34（2）：47-53．

[6] 張國儒，唐亞萍，楊濤，等．272份辣椒資源的農(nóng)藝性狀分析[J]．新疆農(nóng)業(yè)科學(xué)，2021，58（12）：2300-2311．

[7] 楊志剛，胡栓紅，王勇，等．辣椒種質(zhì)色價及主要果實性狀指標(biāo)的聚類分析[J]．西北農(nóng)林科技大學(xué)學(xué)報（自然科學(xué)版），2015，43（5）：149-155．

[8] 李錫香，張寶璽．辣椒種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)[M]．北京：中國農(nóng)業(yè)出版社，2006．

[9] 麥玉連，李惠華，孫勝枚．高效液相色譜法測定辣椒及其制品中辣椒素類物質(zhì)[C]//廣東省食品學(xué)會．健康食品研發(fā)與產(chǎn)業(yè)技術(shù)創(chuàng)新高峰論壇暨2022年廣東省食品學(xué)會年會論文集，2023：223-227．

[10] 桂敏，龍洪進(jìn)，鐘秋月，等．引進(jìn)辣椒種質(zhì)資源果實性狀的多元統(tǒng)計分析[J]．西南農(nóng)業(yè)學(xué)報，2014，27（2）：772-776．

[11] 韓婭楠，行文妍，史艷艷，等．加工型辣椒種質(zhì)資源遺傳多樣性分析[J]．山東農(nóng)業(yè)科學(xué)，2024，56（1）：34-42．

[12] 李寧，王飛，姚明華，等．國內(nèi)外辣椒種質(zhì)資源表型性狀多樣性及相關(guān)性分析[J]．辣椒雜志，2015，13（1）：8-13．

[13] 詹永發(fā)，楊紅，涂祥敏，等．辣椒品種資源的遺傳多樣性和聚類分析[J]．貴州農(nóng)業(yè)科學(xué)，2010，38（11）：12-15．

[14] 袁娟偉，賈利，唐菁，等．275份辣椒種質(zhì)資源表型性狀的遺傳多樣性分析[J]．安徽農(nóng)業(yè)大學(xué)學(xué)報，2024，51（2）：244-255．

[15] 常曉軻，董曉宇，韓婭楠，等．基于主成分分析的不同朝天椒品種品質(zhì)綜合評價[J]．中國瓜菜，2023，36（3）：42-47．

[16] 王春勇，王耐紅，閆穎，等．17個鮮加工型辣椒品種農(nóng)藝性狀分析[J]．種子，2023，42（7）：57-63．

[17] 任朝輝，田旭芳，廖衛(wèi)琴，等．不同辣椒種質(zhì)資源的品質(zhì)性狀評價[J]．西南農(nóng)業(yè)學(xué)報，2020，33（9）：1884-1891．

[18] 鞏雪峰，陳鑫，李紅，等．加工辣椒無損品質(zhì)檢測方法及典型相關(guān)性分析[J]．江蘇農(nóng)業(yè)科學(xué)，2022，50（14）：170-176．

[19] 步洪鳳，張忠武，鄧正春，等．26個辣椒品種主要農(nóng)藝性狀分析與測評[J]．熱帶作物學(xué)報，2019，40（7）：1290-1296．

[20] LI Z X，JIA Z Q，LI J S，et al．Development of a 45 K pepper GBTS liquid-phase gene chip and its application in genome-wide association studies[J]．Frontiers in Plant Science，2024，15：1405190.