基于RAPD的MCID法快速鑒定不結球白菜品種和遺傳多樣性分析

摘要：不結球白菜是長江中下游地區蔬菜周年供應的主栽品種，在種質資源保護和開發利用過程中，品種鑒定的難度隨著市場上品種的增多而增加。利用基于DNA分子標記的人工繪制品種鑒別圖（manual cultivar identification diagram，MCID）方法對筆者所在單位收集和保存的34份不結球白菜種質材料進行鑒定，同時進行NTSYS-pc遺傳多樣性分析。結果表明，利用篩選出的5對RAPD引物擴增出的多態性譜帶即可構建34份不結球白菜的品種鑒定圖，該圖可以快速地篩選出某品種鑒定所需要的引物及需要參考的特征性譜帶，具有較高的可操作性和實用性。聚類結果表明，34份不結球白菜種質材料與品種的聚類模式基本符合其農藝性狀表現和遺傳背景。研究結果可為不結球白菜品種鑒定和種質創新提供了科學依據，對不結球白菜種質資源的保護、研究及品種的早期鑒定具有重要意義。

關鍵詞：不結球白菜；RAPD；人工繪制植物品種鑒別圖（MCID）；品種鑒定

中圖分類號：S634.303 文獻標志碼：A

文章編號：1002-1302（2024）17-0162-07

收稿日期：2023-09-22

基金項目：蘇州市科技局農業科技創新項目（編號：SNG2020065）；江蘇省省級種質資源庫項目（編號：太湖地區特色作物JS-ZW-K18）；蘇州市級種質資源保護項目。

作者簡介：劉照坤（1988—），男，山東鄆城人，碩士，高級農藝師，主要從事不結球白菜種質資源保護、新品種選育及栽培技術研究。E-mail：saaslzk@qq.com。

不結球白菜（Brassica campestris ssp. chinensis），別稱小白菜或青菜，為十字花科蕓薹屬一年或二年生蔬菜作物，原產于我國。生產中選用適宜品種的不結球白菜，可做到四季供應，在蔬菜周年供應中占有十分重要的地位［1］。近年來，不結球白菜在日本、美國及歐洲一些國家被廣泛引種，已逐漸成為世界性蔬菜。

不結球白菜栽培歷史悠久，起源于我國長江中下游太湖地區一帶，在長期自然選擇和人工選擇的演化過程中，形成豐富的遺傳資源［2］。伴隨不結球白菜栽培面積的迅速擴大，品種日益增多，不結球白菜優良種質材料與市售品種在某些植物學特征特性、風味特征上表現出很高的相似性。傳統的形態學鑒定方法難以有效鑒別品種類型，同時品種命名方式多樣，往往存在同名異種或同種異名現象，這增加了種質資源保護和品種鑒定的難度。如市場上常把棵型直立、葉色深綠的不結球白菜種子作為蘇州青品種進行銷售，但田間表現和品質與典型蘇州青相差較大，擾亂了種業市場，影響不結球白菜種業健康發展。因此，本研究嘗試建立不結球白菜不同品種材料的快速可靠鑒定技術體系，對于遺傳資源的研究保護、苗期鑒定具有重要意義。

早期的品種鑒定方法如形態學鑒定和同工酶鑒定等易受外部環境、取材部位、人為因素等影響，而DNA分子標記技術建立在分子水平上，不受外界環境、作物發育階段、取樣部位等因素的影響，且檢出的多態性位點是無限的，鑒定結果具有高度的可靠性和可重復性，鑒別力強。在眾多分子標記技術中，基于聚合酶鏈式反應（polymerase chain reaction，PCR）技術的隨機擴增多態性DNA（random amplified polymorphic DNA，RAPD）技術以分析程序簡便快捷、周期性短、所需費用低等優勢，已被普遍應用于指紋圖譜構建、基因快速定位、遺傳多樣性分析和品種分類等研究［3］。林建麗等利用RAPD技術分析21份白菜類蔬菜種質資源的基因組DNA，通過系統聚類分析將其分為2類10組，結果表明，蕪菁和白菜的親緣關系較遠［4］。韓建明等對國內外不同類型的不結球白菜的遺傳多樣性進行RAPD分析，結果發現，普通白菜的遺傳多樣性最大，不同生態區域遺傳多樣性存在差異［5］。但目前關于利用RAPD分子標記技術鑒定不結球白菜種質資源與品種的研究較少。

人工繪制品種鑒別圖（manual cultivar identification diagram，MCID）法是近幾年廣泛應用的基于RAPD分子標記技術的人工繪制品種鑒別示意圖的方法，該方法可根據PCR擴增后的多態性條帶，實現品種間的區分鑒定。不同于利用計算機軟件繪制數字化指紋圖譜并用統計學軟件進行聚類分析的方法，MCID法更直觀顯現出區分品種的引物，且工作量小，實用性強。目前，該方法已成功在番茄、冬瓜、蘿卜、青花菜、石榴、枇杷、茶樹、蘋果、梨、桃、葡萄等植物上被應用［6-16］，重復性極好，然而，該技術在不結球白菜品種鑒定上仍未見報道。

本研究通過MCID法，對34個不結球白菜種質材料與品種進行鑒定，同時進行遺傳多樣性分析，旨在為不結球白菜的品種鑒定與種質資源利用提供科學依據。

1 材料與方法

1.1 材料

試驗材料為不結球白菜葉片，2022年11月采集于蘇州市農業科學院蔬菜研究所基地，選取每個品種的幼葉，置于冰盒內帶回實驗室用液氮進行速凍后，置于-80 ℃冰箱中保存待用。34個不結球白菜品種及來源見表1。

1.2 方法

1.2.1 DNA提取

利用天根生化科技（北京）有限公司植物基因組DNA提取試劑盒（DP305）提取34個不結球白菜葉片的DNA，利用NanoDropND-2 000 超微量蛋白核酸測定儀（Eppendorf公司，德國）測定DNA樣品濃度，用1.0%瓊脂糖凝膠電泳檢測DNA樣品質量。

1.2.2 引物篩選

試驗所用引物由生工生物工程（上海）股份有限公司合成，采用退火溫度梯度法對35個長度為11個堿基的隨機引物進行篩選，從中得到5個引物（表2）。引物篩選條件：進行連續3次梯度PCR， 選擇譜帶清晰、PCR產物重復出現、質量好的譜帶指紋。

1.2.3 PCR擴增與檢測

PCR反應體系（25 μL）：2×PCR Master Mix（上海萊楓生物科技有限公司）12.5 μL，10 μmol/L引物各1.0 μL，90 ng模板DNA，最后用雙蒸水補至25 μL。反應程序：94 ℃預變性5 min；94 ℃變性30 s，35～45 ℃復性70 s，72 ℃ 延伸2 min，42個循環；72 ℃延伸10 min，4 ℃保存。擴增產物在2%瓊脂糖凝膠中電泳30～50 min，紫外燈下觀察拍照。

1.2.4 DNA指紋圖譜統計分析以及品種鑒別的人工繪圖法

采用PCR擴增產物清晰穩定的引物對34個不結球白菜品種進行鑒定，鑒定方法參考Lin等的MCID法［14］。首先統計某一引物PCR擴增的34個品種的電泳圖，在相同的分子量片段處統計特征性譜帶的有無，將該處有擴增譜帶的品種歸類為一組，沒有譜帶的歸類為另一組，根據該引物其他特征性譜帶進行分組或利用更多的引物逐步進行鑒別；隨后按此方式根據其他引物擴增的多態性譜帶對各小組品種進行鑒別，直至所有品種被單獨鑒別開為止；最后，根據鑒別結果人工繪制鑒別圖，并在MCID中標記每一步所用的引物及多態性譜帶大小（bp）。

1.2.5 數據統計與分析

對選中的每個引物進行3次PCR擴增檢測，選取穩定的條帶進行數據統計與分析。擴增分離結果采用0、1統計方法：同一分子量處，有條帶記為1，無條帶記為0，在Excel中統計擴增出的特異性條帶數目，獲得數據庫。利用NTSYS-pc 2.10e軟件得到遺傳相似系數矩陣，進行UPGMA聚類分析，構建聚類圖，分析品種間的遺傳多樣性。

2 結果與分析

2.1 DNA提取

提取34個不結球白菜葉片的DNA，采用核酸蛋白測定儀測定DNA樣品濃度，確定達到試驗要求后，利用1%瓊脂糖凝膠進行電泳檢測，部分結果見圖1，可以看出，DNA條帶清晰，符合后續試驗要求。

2.2 不結球白菜資源品種鑒定

本試驗采用5條引物擴增出的特異條帶成功地將34個不結球白菜品種逐一區分。首先依據引物Y3在350、300、180 bp處擴增出的3條特異條帶的有無（圖2），將34個不結球白菜分為7組，特異性條帶的有無分別用“+”和“-”表示（圖3）。其中品種5通過180 bp（-）、300 bp（+）、350 bp（+）直接從34個品種中被鑒別出來；品種33通過180 bp（+）、300 bp（-）、350 bp（+）也與其他品種區別開來；第3組為180 bp（-）、300 bp（-）、350 bp（+），包括品種12、32；第4組為180 bp（+）、300 bp（+）、350 bp（-），包括編號為2、6、7、8、13、15、16、26、28、29的10個不結球白菜品種；第5組為180 bp（-）、300 bp（+）、350 bp（-），包括品種19、20、21；第6組為180 bp（-）、300 bp（-）、350 bp（-），包括編號為3、9、18、22、23、24、27、30、31、34的10個品種；第7組為180 bp（+）、300 bp（-）、350 bp（-），包括編號為1、4、10、11、14、17、25的7個品種。

繼續利用更多穩定、特異的引物對未區分開的品種進行鑒定，其中第3組的2個品種可以通過引物Y22在1 000 bp處擴增出的特異片段的有無區別開來。第4組的10個品種可以利用引物Y15在480、550、720 bp處擴增出的3條特異條帶的有無分為4個亞組，第1亞組為480 bp（-）、550 bp（+）、720 bp（+），包括編號為13、26和29的3個品種，繼續利用引物Y34在270 bp處擴增出的條帶的有無把品種26鑒定出來，余下2個品種13和29可以通過引物Y22在720 bp處擴增出的條帶的有無區分開；第2亞組為480 bp（+）、550 bp（+）、720 bp（+），品種2被單獨鑒別出來；第3亞組為480 bp（-）、550 bp（-）、720 bp（+），包括編號8、15、16和28的品種，繼續用Y34擴增出的270 bp條帶分為2組，其中品種8利用270 bp條帶被單獨區分開，

其余3個品種無270 bp條帶被分為一組，繼續用引物Y22擴增出的2個特異條帶可將3個品種區分開來，通過大小為580、1 000 bp的2個片段條帶的有無將15、16、28號品種分別區分開；第4亞組為 480 bp（-）、550 bp（-）、720 bp（-），包括品種6、7號，用引物Y34擴增出的270 bp條帶進行區分，270 bp（+）為6號，無此條帶為7號；同理，第5、6、7組分別用引物Y15、Y22、Y34、Y18等擴增出的特異條帶依次進行鑒別，直至將每個品種都區分開。最后，根據每一步獲得的鑒別結果人工繪制不結球白菜品種鑒定圖，將每一步用到的引物及所利用的多態性譜帶的大小同時標記到相應位置上，即MCID圖。

2.3 品種鑒定結果的驗證

為進一步驗證MCID法鑒定不結球白菜品種的可靠性及可利用性，隨機從不同組中選取編號為2、5、8、15、16的5個品種進行驗證，從MCID中可以看到，通過引物Y3、Y15、Y34和Y22擴增的特異條帶可以區別這5個品種。驗證結果如圖4所示，可以看出，通過Y3引物在180 bp和350 bp處擴增產生的特異條帶的有無可將5號品種與其他4個品種區別出來，通過Y15引物擴增出的480 bp和550 bp特異條帶的有無將2號品種區分出來，進一步用Y34引物擴增出的270 bp條帶將8號品種分離出來，最后利用Y22引物擴增出的580 bp條帶將品種15和16區分開。以上隨機驗證結果證明，優化的RAPD技術體系具有穩定性和可靠性，且MCID具有較高的可行性與實用性。

2.4 聚類分析

對34份不結球白菜種質資源和品種的遺傳多樣性進行分析。由聚類結果（圖5）可知，各品種間遺傳相似系數的變化范圍是0.71～1.00。在遺傳距離0.74處，34份不結球白菜種質資源和品種被分為3個類群，類群Ⅰ包括蘇州青、上海青、矮腳黃、耐熱類青菜等23個品種，數量最多；類群Ⅱ包括蘇州地方特色品種黃葉香青菜和中八葉等4個品種；類群Ⅲ包含黃心烏、五月慢、高梗白、紫色青菜等7個品種，分類結果大致符合品種自身的特征特性，表明各類群存在一定的親緣關系，各品種間存在著遺傳多樣性。

早薹蘇州青和四九菜心因其抽薹性早，遺傳關系較近。蘇州青類型中，A4蘇州青和KT-1蘇州青由收集的蘇州園區外跨塘地方品種提純復壯而來，蘇青14號是以蘇青2號為母本選育而來的品種，親緣關系較近，聚類在一起。AX02蘇州青棵型直立，葉柄稍扁，與上海青、黑矮青遺傳距離較近。大葉蘇州青是蘇州青類型中較耐熱的品種，葉色較其他蘇州青稍淺，與上海青、華王青梗菜等聚類在一起，在今后育種中可作為耐熱材料進行目標選育親本。B25蘇州青為蘇州青類型中抽薹稍晚的材料，M14青梗菜與四月慢主要農藝性狀表現相似，應為近似品種，親緣關系較近；但耐抽薹材料春綠3號、五月慢2個農藝性狀相似的近緣品種卻被聚合在類群Ⅲ，表明四月慢、五月慢可能來源不同，遺傳關系稍遠。類群Ⅱ包括蘇州地方特色品種香青菜中的3個類型，結果表明，黑葉香青菜雖與黃葉香青菜、繡花筋香青菜主要分布區域不同，但親緣關系仍相近；中八葉棵型踏地、葉柄扁平、葉面皺縮呈泡狀，聚類結果顯示與香青菜可能有近緣關系。類群Ⅲ中高梗白（腌菜）與紫色青菜田間表現均為葉片開展、葉柄細長、粗纖維較多；黃心烏和黑心烏來源于安徽，遺傳關系較近。

3 討論與結論

不結球白菜屬異花授粉植物，在自然環境下通過天然雜交形成符合各地消費習慣的地方品種群，長江中下游地區作為我國不結球白菜重要的起源地和消費地區，在長期的繁育過程中形成了大量的優異地方種質資源和品種，品種鑒定是種質資源研究、利用的前提。形態學標記易受生長環境、田間觀測期以及人為主觀因素影響，難以對近似品種進行有效鑒定，因此，準確實用的植物品種鑒定成為農業生產中亟待解決的重要技術難題。分子標記技術因不受環境和人為主觀因素影響， 克服了傳統

鑒定技術的缺點，而更加準確、可靠。目前，已有隨機擴增多態性DNA（random amplified polymorphic DNA，RAPD）［5］、簡單重復序列（simple sequence repeats，SSR）［17-18］、簡單重復序列見擴增（inter-simple sequence repeat，ISSR）［19］等多種DNA分子標記技術應用于不結球白菜品種間遺傳多樣性分析和品種鑒定，為不結球白菜的遺傳進化、品種分類提供了重要的分子證據，但是利用DNA分子標記獲得的指紋圖譜只適用于少數品種的區分鑒定，用NTSYS-pc軟件可計算遺傳距離和遺傳相似性系數，根據UPGMA方法可構建聚類樹狀圖，但無法轉化成鑒定品種的實用工具。本研究利用MCID方法將RAPD分子標記信息轉化為可視化、實用性強的圖表，可有效鑒定大量品種資源，且鑒定結果重復性好。

筆者利用基于RAPD分子標記的MCID法，用5對引物成功將34份已知的不結球白菜種質資源與品種進行鑒定，從MCID中可以明晰不同資源之間進行區分依據的引物和特異條帶，表現出較強的實用性。該方法類似圖書館檢索系統，便于查閱鑒定品種所需的相關信息，查找對應的那對引物擴增的幾個條帶的有無，進而進行針對性的品種鑒定。MCID法提高了引物的利用率，可以用少量的引物將大量的資源與品種區分開來。在已知品種的情況下，向前回溯查找相應引物，根據特征性條帶的大小和有無等信息，可清晰完整地展現所有品種如何被區分開。當有新的不結球白菜品種或種質資源加入時，可直接利用現有的5對引物對該品種進行PCR擴增分析，如發現與現有34個不結球白菜有1個或多個差異性譜帶，則可以把該品種補充到該MCID中去；若發現新的品種難以用5對引物與現有34個品種區分開，則需要添加新的引物對新添加品種進行區分鑒定，該方法工作量小、效率高且長期有效，且可以使不結球白菜MCID品種鑒定系統更加完善，鑒定范圍進一步擴大，成為不結球白菜品種鑒定的工具之一。

統計5對引物在34個不結球白菜種質資源與品種上擴增的特異性條帶，利用相關軟件對其進行遺傳多樣性分析，從聚類分析圖中發現，華王青梗菜與皇冠青梗菜的相似系數接近于1，難以區分；蘇州青類型的資源與品種聚類在類群Ⅰ的第1亞組，遺傳關系較近，且與上海青類型在同一類群；蘇州地方品種香青菜資源和中八葉被分在了類群Ⅱ，但在親緣關系上，黃葉香青菜與黑葉香青菜的親緣關系亦較近，繡花筋香青菜次之，中八葉較遠。研究結果表明，基于RAPD的MCID法可預測不結球白菜品種間的親緣關系或遺傳距離，推測結果符合品種自身特征特性。

本研究利用MCID法對不結球白菜品種進行區分鑒定，該方法具有快速、簡便、可重復的優點，可推測品種間親緣關系遠近和遺傳多樣性。研究結果表明，MCID不僅能為植物品種或材料鑒定提供幫助，還能為種質資源研究、品種權保護、實際生產提供幫助，為不結球白菜種質資源保護、鑒定和利用提供理論依據。

參考文獻：

［1］侯喜林，李 英，劉同坤.不結球白菜遺傳育種與分子生物學研究進展［J］. 南京農業大學學報，2022，45（5）：864-873.

［2］葉靜淵. 我國油菜的名實考訂及其栽培起源［J］. 自然科學史研究，1989，8（2）：158-165.

［3］朱旭東，上官凌飛，孫 欣，等. DNA標記在植物品種鑒定上的應用現狀［J］. 中國農學通報，2014（30）：234-240.

［4］林建麗，朱正歌，王建成，等. 白菜類蔬菜種質資源親緣關系的RAPD分析［J］. 華北農學報，2008，23（增刊2）：183-189.

［5］韓建明，侯喜林，徐海明，等. 不結球白菜種質資源遺傳多樣性RAPD分析［J］. 南京農業大學學報，2008，31（3）：31-36.

［6］Cao X，Wu Z，Zhou R，et al. A novel random amplified polymorphic DNA-based strategy for genetic diversity analysis and identification of tomatoes［J］. Genetics and Molecular Research，2015，14（1）：1650-1661.

［7］葉新如，劉建汀，李永平，等. 基于EST-SSR標記的MCID法鑒定冬瓜種質資源［J］. 核農學報，2021，35（4）：780-788.

［8］許園園，劉 哲，婁麗娜，等. 基于MCID法的蘿卜品種快速鑒定［J］. 江蘇農業學報，2016，32（6）：1384-1389.

［9］王葉群，袁金超，吳莎莎，等. 基于RAPD的MCID法快速鑒定青花菜品種和遺傳多樣性分析［J］. 分子植物育種，2019，17（23）：7778-7786.

［10］Zhang Y P，Tan H H，Cao S Y，et al. A novel strategy for identification of 47 pomegranate （Punica granatum） cultivars using RAPD markers［J］. Genetics and Molecular Research，2012，11（3）：3032-3041.

［11］孫 鈞，李曉穎，徐紅霞，等. 基于genic-SSR標記的MCID法鑒定浙江白沙枇杷地方種質資源［J］. 果樹學報，2018，35（5）：539-547.

［12］黃丹娟，馬建強，馬春雷，等. 基于SSR熒光標記的MCID快速鑒定13個湖北茶樹品種［J］. 南方農業學報，2018，49（6）：1045-1052.

［13］Wang W Y，Wang K，Liu F Z，et al. An efficient identification of 68 apple cultivars using a cultivar identification diagram （CID） strategy and RAPD markers［J］. Korean Journal of Horticultural Science and Technology，2012，30（5）：549-556.

［14］Lin J，Wang X C，Chang Y H，et al. Development of a novel and efficient strategy for practical identification of Pyrus spp.（Rosaceae） cultivars using RAPD fingerprints［J］. Genetics and Molecular Research，2011，10（2）：932-942.

［15］Yu M L，Wang W Y，Ma R J，et al. An improved strategy based on RAPD markers efficiently identified 95 peach cultivars［J］. Genetics and Molecular Research，2012，11（2）：1158-1168.

［16］張中海，張安世，高登濤，等. 巨峰系葡萄主栽品種的SRAP遺傳多樣性分析與MCID法快速鑒定［J］. 東北農業科學，2022，47（4）：38-42.

［17］李瑞峰，高 鵬，朱子成，等. 基于形態學標記及SSR標記的甜瓜主栽品種分類鑒定研究［J］. 中國蔬菜，2014（6）：20-27.

［18］劉冬媛，王新華，劉 楊，等. 不結球白菜品種（系）遺傳多樣性與親緣關系的SSR分析［J］. 分子植物育種，2014，12（3）：499-508.

［19］Wang A H，Zhou G L，Lin C F，et al. Genetic diversity study of Brassica campestris L.ssp.chinensis Makino based on ISSR markers［J］. Caryologia，2017，70（1）：48-54.