侵染江蘇獼猴桃的北方根結線蟲（Meloidogyne hapla）形態學描述和分子特征分析

范亞磊　趙敏　鄧晟　劉瑞顯　姚煥釗　魏利輝　馮輝

摘要： 根結線蟲病長期以來對獼猴桃安全生產造成了嚴重威脅。調查發現江蘇省連云港市部分獼猴桃種植園根結線蟲病發生普遍，組織切片可見根部出現根結線蟲誘導的特異的巨細胞。形態學觀察和分子檢測結果顯示，盡管與其他群體存在一定差異，但仍確定本次分離自獼猴桃的根結線蟲（種群編號：CN43）為北方根結線蟲（Meloidogyne hapla）。基于內轉錄間隔區（ITS）和COⅡ-16S序列構建的遺傳進化樹顯示，CN43分離群體與已知的北方根結線蟲位于同一進化分支，而與其他種類的根結線蟲明顯區分開來。以ITS序列構建的中接網絡將34個北方根結線蟲分離群體劃分為10個單倍型，其中CN43屬于最大的單倍型群組;然而，單倍型分析不能確定各分離群體在地理分布上的關聯性。綜上，本研究首次報道了北方根結線蟲在中國獼猴桃上的侵染和發生，其分離群體的形態學描述和分子特征分析可為揭示獼猴桃根結線蟲的種群多樣性和種間變異提供參考。

關鍵詞： 獼猴桃;根結線蟲;北方根結線蟲;內轉錄間隔區（ITS）;COⅡ;進化分析

中圖分類號： S436.634 文獻標識碼： A 文章編號： 1000-4440（2021）01-0075-08

Morphological description and molecular characteristic analysis on Meloidogyne hapla which infect kiwifruit in Jiangsu province

FAN Ya-lei1，2， ZHAO Min1， DENG Sheng1， LIU Rui-xian3，YAO Huan-zhao4， WEI Li-hui1， FENG Hui1

（1.Institute of Plant Protection， Jiangsu Academy of Agricultural Sciences， Nanjing 210014，China;2.School of Horticulture and Plant Protection， Yangzhou University， Yangzhou 225009，China;3.Institute of Industrial Crops， Jiangsu Academy of Agricultural Sciences， Nanjing 210014，China;4.Plant Protection Station of Tongshan District of Xuzhou City， Xuzhou 221116，China）

Abstract： Root-knot diseases caused by Meloidogyne spp. have posed great threat to global kiwifruit industry for a long time. It was found through investigation that， the disease caused by root-knot nematode was common in part of the kiwifruit plantations in Lianyungang City， Jiangsu province. The specific giant cells were found in tissue slices of roots， which may be induced by root-knot nematodes. Results of morphological observation and molecular detection indicated that， the root-knot nematode （population code： CN43） separated from kiwifruits was confirmed to be Meloidogyne hapla， although there was difference compared with other polulations. Phylogenetic trees inferred from internal transcribed spacer （ITS） and partial COⅡ-16S sequences revealed that， isolate population of CN43 and M. hapla were within the same clade and could be separated from other root-knot nematode species obviously. Isolates of 34 M. hapla in NCBI database were divided into ten haplotypes using median-joining network constructed by ITS sequence， in which CN43 belonged to the largest haplotype group. However， haplotype analysis could hardly show the correlation of geographic distribution between the isolates. In conclusion， infection and occurrence of M. hapla in kiwifruit trees of China was reported for the first time. Analysis on the morphology and molecular characteristics of the isolated M. hapla populations can provide reference for revealing population diversity and intraspecific variability of root-knot nematodes in kiwifruit.

Key words： kiwifruit;root-knot nematode;Meloidogyne hapla;internal transcribed spacer（ITS）;COⅡ;phylogenetic analysis

獼猴桃（Actinidia deliciosa）是起源于中國的適合在溫熱帶地區栽培的重要果樹，目前已在20多個國家得到種植，截至2011年，中國以外地區獼猴桃收獲面積達到9.4×104 hm2[1]。然而，隨著獼猴桃產業和規模的不斷發展，獼猴桃根結線蟲病害問題日益凸顯。根結線蟲（Meloidogyne spp.）在全世界各獼猴桃種植區廣泛存在，其中南方根結線蟲（M. incognita）和北方根結線蟲（M. hapla）是主要病原線蟲種類。南方根結線蟲在巴西、智利、印度、美國、土耳其等國均有分布，北方根結線蟲則分布在巴西、智利、印度、意大利、西班牙、新西蘭、韓國等國家，在獼猴桃上存在和為害[2]。此外，其他根結線蟲種類如爪哇根結線蟲（M. javanica）、花生根結線蟲（M. arenaria）、埃塞俄比亞根結線蟲（M. ethiopica）[3]等也能侵染獼猴桃。

根結線蟲在中國主要獼猴桃產區長期為害，并造成嚴重的經濟損失。過去30年來，南方根結線蟲是中國獼猴桃上的優勢種群，在陜西[4]、河南[5-6]、湖北[7]、湖南[8]、貴州[9]、浙江[10]、福建[11-12]等獼猴桃產地普遍存在。部分獼猴桃產地還發現了爪哇根結線蟲[11]、花生根結線蟲[7]、M. actinidiae[13]和M. aberrans[2]等種類。然而，北方根結線蟲在中國獼猴桃上發生和為害的情況卻鮮有報道。

近年來，隨著中國農業產業結構調整和供給側結構性改革的推進，獼猴桃作為農民增收脫貧的重要產業在江蘇省得到快速發展。本研究于2019年7月對江蘇省多個獼猴桃種植基地根結線蟲病進行了調查，發現連云港部分果園的獼猴桃地上植株纖弱，葉片褪綠枯黃，坐果少，地下根系畸形，新根布滿大量珍珠狀根結，表現出典型的根結線蟲病的癥狀。通過病理解剖，從發病根部分離出根結線蟲。利用形態學觀察和分子檢測標記對其進行鑒定，確定侵染獼猴桃的線蟲種類為北方根結線蟲。同時，基于內轉錄間隔區（ITS）和COⅡ-16S序列分析了該群體的分子特征以及與其他分離群體間的進化關系。

1 材料與方法

1.1 線蟲分離

線蟲樣本采集自江蘇省連云港市贛榆區黑林鎮獼猴桃種植園。選擇長勢較弱、葉片發黃的獼猴桃植株，小心挖取根組織，并選取具有典型根結癥狀的根組織帶回實驗室進行后續分析。

1.2 形態學鑒定

用清水漂洗根組織以去除泥沙，在體式顯微鏡下，用挑針從根部根結處挑取根結線蟲雌蟲和雄蟲，二齡幼蟲直接由卵塊孵化獲得。線蟲殺死、固定、會陰花紋等標本制作參照文獻[13]。通過光學顯微鏡（DM2500，徠卡顯微系統有限公司產品）和體式顯微鏡（M205FA，徠卡顯微系統有限公司產品）觀察并拍攝雌蟲、雄蟲和幼蟲形態，測量de-Man值。

1.3 線蟲DNA提取和PCR擴增

DNA提取參考文獻[14]中的方法進行，略有修改。在解剖鏡下，挑取單條雌蟲至盛有20 μl裂解液（含1×PCR 緩沖液、1 μg蛋白酶K）的PCR管中，用鑷子將蟲體擠破使蟲液釋放，混合均勻后置于液氮中迅速冷凍20 min。隨后將樣品置于PCR儀上于65 ℃孵育30 min， 95 ℃ 5 min，最后于12 000 r/min離心1 min，吸取上清液用于PCR擴增。

核糖體ITS序列用引物TW81（5′-GTTTCCGTAGGTGAACCTGC-3′）和AB28（5′-ATATGCTTAA-GTTCAGCGGGT-3′）擴增，覆蓋3′-18S rRNA、ITS1、5.8S rRNA、ITS2和5′-28S rRNA序列[15]。25.0 μl反應體系含有2×Taq Master mix（北京康為世紀生物科技有限公司產品）12.5 μl，正反向引物各10 pmol，線蟲DNA粗提液5.0 μl。PCR反應條件為：94 ℃ 5 min;94 ℃ 30 s，55 ℃ 30 s，72 ℃ 60 s，共38個循環;72 ℃ 5 min。線粒體細胞色素c氧化酶Ⅱ亞基基因COⅡ序列應用C2F3（5′-GGTCAATGTTCAGAAATTTGTGG-3′）和MRH106-（5′-AATTTCTAAAGACTTTTCTTAGT-3′）擴增，覆蓋部分COⅡ編碼序列、tRNA-His序列以及部分16S rRNA序列[16]。反應體系同上，PCR反應條件為：94 ℃ 4 min;94 ℃ 1 min，50 ℃ 1 min，70 ℃ 2 min，共35個循環;72 ℃ 10 min。另外，利用南方根結線蟲SCAR（序列特異性擴增區域）標記Mi-F（5′-GTGAGGATTCAGCTCCCCAG-3′）和Mi-R（5′-ACGAGGAACATACTTCTCCGTCC-3′）[17]，北方根結線蟲SCAR標記Mh-F（5′-TGACGGCGGTGAGTGCGA-3′）和Mh-R（5′-TGACGGCGGTACCTCATAG-3′）[18]進一步驗證分離群體的種類。PCR產物經1.2%瓊脂糖凝膠電泳檢測后直接進行測序。

1.4 序列比對和進化分析

分別以江蘇分離群體ITS和COⅡ序列在美國國家生物技術信息中心（NCBI）數據庫進行BLAST（http：//www.ncbi.nlm.nih.gov/blast）比對，根據E值篩選相似度相近的序列以及合適外群序列構建貝葉斯進化樹。將新獲得的序列和下載序列通過MAFFT軟件（https：//www.ebi.ac.uk/Tools/msa/mafft）進行多序列比對，再通過Gblocks 0.91b（http：//www.phylogeny.fr/one_task.cgi？task_type=gblocks）進行策展分析并獲取保守序列。經DAMBE軟件運算發現，所有序列的核苷酸替換未達到飽和，適宜建樹。基于最小赤池信息準則（AIC）值，經Modeltest 3.7計算獲得ITS序列最優堿基替換模型為GTR+G，COⅡ最優堿基替換模型為K81uf+G。采用馬爾科夫鏈蒙特卡洛方法 （Markov chain Monte Carlo，MCMC）構建貝葉斯樹，每條馬爾科夫鏈運行8×105代。拋棄25%老化樣本進行聚合評估，剩余樣本用于后續分析構建一致性為50%的發育樹，并經FigTree 1.4進行優化。另對NCBI上北方根結線蟲各分離群體ITS序列進行重新分析，去除插入和缺失序列，最終確定34個群體樣本有11個核苷酸發生置換，基于該置換位點，利用Network v10構建單倍型中接網絡（Median-joining network）[19]。

2 結果與分析

2.1 受害獼猴桃根部癥狀和線蟲形態特征

受害獼猴桃植株長勢較弱，部分葉片褪綠枯黃，坐果量減少（圖1A）;地下根系不發達、畸形，新根布滿大量珍珠狀根結，表現出典型的根結線蟲病的癥狀（圖B）。利用酸性品紅進行組織染色可見根結處明顯出現近球形根結線蟲（圖1C）。通過病理石蠟切片，可見受害植株根內不同齡期的線蟲蟲態及根結線蟲特異誘導的巨細胞（圖1D、圖1E）。

在顯微鏡下對所分離的線蟲（標記為CN43分離群體）進行觀察，主要形態特征如下：雌蟲體型近球形，頸部短，口針纖細、基球小，錐部向背部稍彎曲，末端寬，口針基球圓形，與桿部有明顯界限;頭區大、無環紋（圖2A）;會陰花紋近圓形六邊形，背弓扁平，側線不明顯，線紋平滑形到波浪形，尾端有刻點（圖2B、圖2C）。雄蟲體型較短，頭冠圓、高，頭區無環紋，整個輪廓外突，與體區分界較明顯，口針椎體前端細尖，向后漸粗;口針基球球狀，與基桿界限明顯;尾端光滑鈍圓，側尾線口位于泄殖腔附近，交合刺略彎曲（圖2D、圖2E）。二齡幼蟲頭冠圓、窄，口針基球圓形，與基桿分界不明顯;尾細長，末端透明區薄，界限不明顯（圖2F、圖2G）。雌蟲、雄蟲和二齡幼蟲de-Man形態學參數與已報道的北方根結線蟲（M. hapla）特征值（表1）基本一致。

2.2 分子生物學特征和進化分析

由于在中國獼猴桃種植區南方根結線蟲普遍存在，因此首先用南方根結線蟲特異引物進行擴增，電泳未顯示出條帶，而以南方根結線蟲DNA為模板的陽性對照則獲得單一的995 bp左右的條帶。以ITS和COⅡ-16S通用引物PCR擴增分別獲得553 bp（NCBI登錄號為MT490918）和663 bp條帶（NCBI登錄號為MT536934）。利用北方根結線蟲特異引物PCR擴增獲得單一的610 bp條帶，與預期大小一致（圖3）。通過與NCBI數據庫登錄的序列進行比對分析，發現CN43分離群體ITS序列與其他北方根結線蟲分離群體相似度為85.19%～99.64%，其中與中國云南分離群體（NCBI登錄號：JX024147）相似度最高，為99.64%（E值=0），而與日本分離群體（NCBI登錄號：LC030358）相似度最低，為85.19%（E值=3e-155）。CN43分離群體COⅡ-16S序列與澳大利亞分離群體（NCBI登錄號：L76262）相似度最高，為96.69%（E值=0），而與埃塞俄比亞分離群體（NCBI登錄號：KP681262）相似度最低，為92.78%（E值=0）。

基于ITS和COⅡ-16S序列構建的貝葉斯進化樹顯示，盡管北方根結線蟲各分離群體位于不同較短的次級分支，但總體上CN43與北方根結線蟲劃歸為一個大的進化分支，與其他種類的線蟲位于不同分支（圖4、圖5）。基于NCBI數據庫中所有北方根結線蟲ITS序列進行重新分析，發現其中5個18S基因序列存在高度變異，推測這些遞交至NCBI上的序列可能不完整或不正確，因此將其排除，最終用于分析的序列共34個，包括中國分離群體20個，韓國和美國各3個，日本和澳大利亞各2個、新西蘭、塞爾維亞、西班牙、北愛爾蘭各1個。構建的中接網絡共獲得10個單倍型，每個單倍型代表1～21個分離群體。其中單倍型1出現頻率最高（21個群體），與其他9個單倍型相連，表明這些種群可能是其他單倍型的進化祖先。而中國[臺灣（NCBI登錄號：MK774795）、浙江（NCBI登錄號：KR535973）、云南（NCBI登錄號：MT249019）]和日本（NCBI登錄號：LC030361）、韓國（NCBI登錄號：MK188479）、澳大利亞（NCBI登錄號：AF516721）等群體存在單個獨立分支。由于無法排除測序誤差，將其分支用虛線表示;中接網絡圖顯示，各個單倍型和地理分布之間無明顯相關性（圖6）。

3 討論

根結線蟲在中國各個獼猴桃產區普遍發生。林尤劍等[11]報道福建省獼猴桃上的根結線蟲主要為爪哇根結線蟲和南方根結線蟲，其中爪哇根結線蟲為優勢種。姜鳳麗等[10]發現，浙江省獼猴桃根結線蟲主要為南方根結線蟲、爪哇根結線蟲和花生根結線蟲，其中南方根結線蟲為優勢種。陳文等[9]對貴州獼猴桃產區進行調查，基于雌蟲會陰花紋特征結合分子標記鑒定出所有分離群體均為南方根結線蟲。

近年來，獼猴桃種植作為江蘇部分地區農村振興和農民脫貧致富的重要產業得到迅速發展。然而，由于大量引種和管理措施不當，根結線蟲病逐漸成為制約當地獼猴桃安全生產的重要因素之一。筆者前期對江蘇省多個獼猴桃種植基地進行了調查，發現部分獼猴桃果園根結線蟲病發生較重，發病新根形成許多小珍珠般大小的根結，根結初與新根同色，后逐漸變為褐色，甚至壞死、腐爛;根際土壤蟲口密度最高可達每100 g土壤25～30條二齡幼蟲。通過SCAR分子標記進行檢測，基于特異的電泳條帶確定其主要為南方根結線蟲，表明南方根結線蟲是為害江蘇獼猴桃的主要線蟲種類。然而基于形態學觀察和分子生物學數據，證實其中1個來自連云港的分離群體（CN43），為北方根結線蟲——這是北方根結線蟲侵染中國獼猴桃的首次報道。

本研究分離的CN43群體與其他北方根結線蟲分離群體在形態特征上基本一致，例如雌蟲會陰花紋存在刻點以及二齡幼蟲尾部透明區界限不明顯等。盡管如此，仍與其他分離群體的de-Man參數值存在一定差異。例如，與美國夏威夷咖啡分離群體相比，CN43群體雌蟲、雄蟲SPI值、二齡幼蟲L值和W值較大，而雌蟲MEV-HE值和二齡幼蟲a值較小[20]。與國內記載的北方根結線蟲特征相比，CN43群體的雄蟲HH值、HW值和二齡幼蟲L值比遼寧沈陽紫花地丁分離群體大，而二齡幼蟲ST值、MEV-HE值和a值比遼寧沈陽紫花地丁分離群體小;雄蟲的HH值和a值比山東威海花生和濟南薄荷上的分離群體略小[21]。一般認為，北方根結線蟲存在A、B 2個細胞遺傳小種，其判定依據主要根據染色體數目和雄蟲頭部特征。盡管本次獲得的CN43群體雄蟲數目十分有限，然而這些雄蟲的口針基桿與基球連接處略微變粗，與北方根結線蟲A小種雄蟲描述較為一致[22]。宋志強[23]曾從江蘇連云港地區蔬菜上分離到北方根結線蟲，但缺乏形態學描述和分子特征分析，尚無法確定CN43群體與侵染當地蔬菜的北方根結線蟲為同一小種。

核糖體ITS和線粒體基因通常用于揭示根結線蟲中間變異特征[24-25]。盡管有研究結果表明ITS序列在根結線蟲屬內識別率較低[26]，然而本研究結果表明，ITS和COⅡ-16S序列能較好區分北方根結線蟲和其他種類根結線蟲。同時，序列比對和貝葉斯進化樹結果表明，CN43群體與NCBI登錄的北方根結線蟲均有高度相似性。由于NCBI數據庫中COⅡ數據十分有限，難以鑒別北方根結線蟲種間變異特征。因此，本研究采用ITS序列研究北方根結線蟲各個分離群體單倍型特征。基于中接網絡，確定CN43與其他20個分離群體歸為一類單倍型，然而由于潛在的測序錯誤，尚不能確定各分離群體在地理分布上的關聯性。新型分子標記的開發和應用有望消除這些誤差，以揭示其單倍型分布規律。

參考文獻：

[1] 張計育. 21世紀以來世界獼猴桃產業發展以及我國獼猴桃貿易與國際競爭力分析[J]. 中國農學通報， 2014， 30（23）：48-55.

[2] TAO Y， XU C， YUAN C， et al. Meloidogyne aberrans sp. nov.（Nematoda： Meloidogynidae）， a new root-knot nematode parasitizing kiwifruit in China[J]. PLoS One， 2017， 12（8）： e0182627.

[3] COFCEWICZ E， ALMEIDA M R， ABALLAY E， et al. Meloidogyne ethiopica， a major root-knot nematode parasitising Vitis vinifera and other crops in Chile[J]. Nematology， 2007， 9（5）： 633-639.

[4] 范學科，黨占平. 獼猴桃根結線蟲病的病原鑒定及其防治[J]. 陜西農業科學， 2007（6）： 71-72.

[5] 李淑君，喻 璋. 河南獼猴桃根結線蟲新種[J]. 河南農業大學學報， 1991， 3（6）： 71-72.

[6] XU J H， LIU P L， MENG Q P， et al. Characterisation of Meloidogyne species from China using isozyme phenotypes and amplified mitochondrial DNA restriction fragment length polymorphism[J]. European Journal of Plant Pathology， 2004， 110： 309-315.

[7] 楊清平，王立華，謝志斌，等. 湖北獼猴桃主要病害及其有機病害治理技術[J]. 湖北農業科學， 2014， 53（10）： 2307-2311.

[8] 方炎祖，羅桂菊，朱曉香，等. 湖南獼猴桃根結線蟲病害研究[J]. 湖南農業科學， 1991， 30（4）： 40-42.

[9] 陳 文，孫燕芳，吳石平，等. 貴州修文獼猴桃根結線蟲的發生種類與鑒定[J]. 西南農業學報， 2018， 31（1）： 84-88.

[10]姜鳳麗，邵桂英. 獼猴桃根結線蟲病的初步研究[J]. 浙江林學院學報， 1990（1）： 46-51.

[11]林尤劍，高日霞. 福建獼猴桃病害調查與鑒定[J]. 福建農業大學學報， 1995， 24（1）： 49-53.

[12]張紹升，林尤劍，高日霞. 福建獼猴桃根結線蟲病病原鑒定[J]. 福建農學院學報， 1993， 22（4）： 433-435.

[13]劉維志. 植物病原線蟲學[M]. 北京： 中國農業出版社， 2000： 384-400.

[14]馮 輝，魏利輝，陳懷谷，等. 細尖潛根線蟲（Hirschmanniella mucronata）江蘇分離群體形態學和分子特征描述[J]. 植物病理學報， 2016， 46（4）： 474-484.

[15]MAAFI Z T， SUBBOTIN S A， MOENS M. Molecular identification of cyst-forming nematodes （Heteroderidae） from Iran and a phylogeny based on ITS-rDNA sequences[J]. Nematology， 2003， 5（1）： 99-111.

[16]STANTON J， JUGALL A， MORITZ C. Nucleotide polymorphisms and an improved PCR-based mtDNA diagnostic for parthenogenetic root-knot nematodes （Meloidogyne spp.）[J]. Fundamental and Applied Nematology， 1997， 20（3）： 261-268.

[17]MENG Q， LONG H， XU J. PCR assays for rapid and sensitive identification of three major root-knot nematodes， Meloidogyne incognita， M. javanica and M. arenaria[J]. Acta Phytopathologcia Sinica， 2004，34：204-210.

[18]ZIJLSTRA C. Identification of Meloidogyne chitwoodi， M. fallax and M. hapla based on SCAR-PCR： a powerful way of enabling reliable identification of populations or individuals that share common traits[J]. European Journal of Plant Pathology， 2000， 106： 283-290.

[19]BANDELT H J， FORSTER P， RHL A. Median-joining networks for inferring intraspecic phylogenies[J]. Molecular Biology and Evolution， 1999， 16（1）： 37-48.

[20]HANDOO Z A， SKANTAR A M， CARTA L K， et al. Morphological and molecular evaluation of a Meloidogyne hapla population damaging coffee （Coffea arabica） in Maui， Hawaii[J]. Journal of Nematology， 2005， 37（2）： 136-145.

[21]趙洪海. 中國部分地區根結線蟲的種類鑒定和四種最常見種的種內形態變異研究 [D]. 沈陽：沈陽農業大學， 1999.

[22]劉維志. 植物線蟲志[M]. 北京： 中國農業出版社，2004： 372-380.

[23]宋志強. 江蘇省蔬菜根結線蟲的分布調查、種類鑒定及定量檢測技術研究 [D]. 南京：南京農業大學， 2013.

[24]POWERS T， MULLIN P G， HARRIS T S， et al. Incorporating molecular identication of Meloidogyne spp. into a large-scale regional nematode survey[J]. Journal of Nematology， 2005， 37： 226-235.

[25]POWERS T. Nematode molecular diagnostics： from bands to barcodes[J]. Annual Review of Phytopathology， 2004， 42： 367-383.

[26]林 宇，王金成，遲元凱，等. 基于GenBank分析28S（D2/D3）、18S和ITS序列作為根結線蟲條形碼標記的適用性[J]. 南京農業大學學報， 2013， 36（5）： 71-76.

（責任編輯：陳海霞）

收稿日期：2020-06-20

基金項目：國家自然科學基金項目（31871943）; 江蘇省農業科技自主創新基金項目[CX（18）2005]; 江蘇省農業科學院“小而特”學科建設項目[ZX（19）6005]

作者簡介：范亞磊（1994-），男，江蘇淮安人，碩士研究生，從事植物線蟲病害防治研究。

通訊作者：馮 輝，（Tel）025-84390769;（E-mail）fenghui@jaas.ac.cn