連云港地區(qū)野生靈芝生物學(xué)特性及基因組特征分析

摘要：野生靈芝資源在連云港境內(nèi)有著廣泛的分布，本研究使用全基因組重測序技術(shù)對具有不同形態(tài)特征和生境特點的靈芝基因組進行測序及分析，以期解析連云港地區(qū)野生靈芝的生物學(xué)特性及基因組特征。結(jié)果表明，從不同區(qū)域采集的野生靈芝種質(zhì)資源菌絲具有鎖狀聯(lián)合，具備形成子實體的條件，生物學(xué)特性存在較大差異；通過主成分分析將12株野生靈芝菌株種質(zhì)資源聚為3類，其中第Ⅲ類具有更豐富的遺傳多樣性；比較單核苷酸多態(tài)性（SNP）、插入或者缺失（InDel）在12株靈芝菌株基因組中的分布情況發(fā)現(xiàn)，第1條染色體中的SNP最多，占比為46.07%，SNP中轉(zhuǎn)換顛換比（Ts/Tv）為2.03，說明變異以轉(zhuǎn)換為主，插入或缺失導(dǎo)致的基因差異在第1條染色體上最多，占總InDel數(shù)的44.31%。研究結(jié)果說明，連云港地區(qū)靈芝具有豐富的遺傳多樣性，可以用于進一步挖掘更多的功能基因和進行品種選育。

關(guān)鍵詞：靈芝；生物學(xué)特性；全基因組重測序；主成分分析

中圖分類號：S567.3+1文獻標(biāo)識碼：A文章編號：1000-4440（2024）02-0223-10

Analysis of biological characteristics and genomic characteristics of wild Ganoderma lucidum in Lianyungang

JI Wei1，2，SU Wen-ying1，LIU Xiao-mei1 ，REN Li-kai1，HU Shu-yun3，SUN Xiao-xiao3，CHEN Ke-long2

（1.Lianyungang Academy of Agricultural Sciences， Lianyungang 222000， China;2.Qinghai Normal University， Xining 810000， China;3.Lianyungang Gardening Vegetable Guidance Station， Lianyungang 222000， China）

Abstract：Wild Ganoderma lucidum resources are extensively found in Lianyungang. In this study， we utilized whole-genome resequencing technology to examine and analyze the genomes of Ganoderma lucidum with different morphological and habitat characteristics. The purpose of this study was to unveil the biological and genomic characteristics of wild Ganoderma lucidum in Lianyungang. The results showed that the mycelium of wild Ganoderma lucidum germplasm resources collected from different regions had clamp connection， and had the conditions to form fruiting bodies， and the biological characteristics were quite different. Through principal component analysis， 12 wild Ganoderma lucidum germplasm resources were grouped into three categories， and the third category had more abundant genetic diversity. Analyzing the distribution of single nucleotide polymorphism （SNP） and insertion or deletion （InDel） in the genome of 12 strains of Ganoderma lucidum， it was found that the number of SNP on chromosome 1 was the largest， accounting for 46.07%， and the conversion and inversion ratio （Ts/Tv） in the SNP was 2.03， indicating that the variation was mainly conversion. The genetic difference caused by insertion or deletion was the highest on chromosome 1， accounting for 44.31% of the total InDel number. These results showed that Ganoderma lucidum in Lianyungang had rich genetic diversity， which could be used to further explore more functional genes and breed varieties.

Key words：Ganoderma lucidum;biological characteristics;whole genome resequencing;principal component analysis

靈芝是一種大型藥食兩用真菌[1]，是馳名中外的珍稀中藥材之一[2]，其藥理成分豐富，且無毒副作用，有效成分包括靈芝多糖[3]、多肽[4]、三萜類[5]及16種氨基酸[6]等。靈芝能降低中樞神經(jīng)系統(tǒng)的興奮性，有一定鎮(zhèn)痛作用[7]，此外還有解毒[8]、降血糖[9]、抗輻射[10]、提高免疫力[11]、治療哮喘[12]和抗腫瘤[13]等作用。靈芝還具有極高的文化觀賞價值[14]，可將之培育成外形美觀的觀賞盆景[15]。

連云港位于江蘇北部，境內(nèi)山海齊觀，河湖、丘陵、灘涂、濕地、海島俱備[16]，具有豐富的生物多樣性，在菌物資源方面，連云港境內(nèi)已報道的大型真菌資源有50余種[17]，其中野生靈芝資源豐富、類型多樣，當(dāng)?shù)卮迕穸嗖杉吧`芝作為普通農(nóng)產(chǎn)品或中藥材在市場上銷售，但相關(guān)研究甚少。本研究團隊近幾年對連云港地區(qū)野生靈芝資源開展了一系列育種研究，在采集、鑒定、馴化栽培、栽培條件優(yōu)化及抗逆等方面做了一部分工作，探索并實踐了代料栽培、椴木栽培和近地保護方式的仿野生栽培，但沒有科學(xué)有效地從遺傳學(xué)角度選擇性地對這部分資源進行研究、開發(fā)和保護，對連云港地區(qū)不同區(qū)域采集的靈芝資源的生物學(xué)特性、基因組特征和遺傳多樣性尚不清楚。因此，本研究擬選取具有代表性的野生靈芝菌株進行分離、鑒定，對不同菌株間菌落生長速度、液體發(fā)酵生物量、液體發(fā)酵pH值和人工栽培吃料速度等生物學(xué)特性進行研究，首次采用全基因組重測序方法對不同靈芝菌株之間基因組特征及遺傳多樣性進行分析，以期明確其在分子水平上的遺傳多態(tài)性與親緣關(guān)系，為后期通過分子標(biāo)注技術(shù)與靈芝生物學(xué)特性、表型性狀評價相結(jié)合創(chuàng)造優(yōu)質(zhì)種質(zhì)資源及挖掘關(guān)鍵功能基因提供理論基礎(chǔ)。

1材料與方法

1.1試驗材料

本項目組于2022年6-9月在江蘇省連云港市范圍內(nèi)采集獲得野生靈芝標(biāo)本，根據(jù)形態(tài)特征、生境特點選取具有代表性的標(biāo)本12個，野生靈芝標(biāo)本信息如表1所示。將采集的野生靈芝進行組織分離得到菌株，保存于連云港市農(nóng)業(yè)科學(xué)院的食用菌菌種庫中。

試驗試劑：馬鈴薯葡萄糖水培養(yǎng)基，購自青島海博生物公司；馬鈴薯葡萄糖瓊脂（PDA）培養(yǎng)基，購自北京陸橋技術(shù)股份有限公司；蛋白胨，購自生工生物工程（上海）股份有限公司；酵母浸粉、瓊脂粉，購自北京奧博星生物技術(shù)有限責(zé)任公司。

試驗儀器：CX43顯微鏡，購自日本Olympus公司；超凈工作臺，購自德國Airtech公司；PCR儀、MD-550離心機，購自德國Eppendorf公司。

1.2試驗方法

1.2.1靈芝菌株基因組的提取及分子鑒定采用真菌基因組DNA快速提取試劑盒提取基因組DNA，采用真菌通用引物ITS1（5′-TCCGTAGGTGAACCTGCGG-3′）和ITS4（5′-TCCTCCGCTTATTGATATGC-3′）進行PCR擴增。擴增條件如下：94 ℃預(yù)變性5 min；94 ℃變性30 s，60 ℃ 退火30 s，72 ℃延伸1 min，35個循環(huán)；72 ℃延伸10 min。對PCR產(chǎn)物進行測序，用MEGA 7.0對所測內(nèi)轉(zhuǎn)錄間隔區(qū)（ITS）序列進行比對，用鄰接（Neighbor joining）法構(gòu)建系統(tǒng)發(fā)育樹[18]。

1.2.2靈芝菌株生物學(xué)特性比較分析

1.2.2.1菌絲體干質(zhì)量的測定量取50 ml靈芝發(fā)酵液倒入離心管中，10 000 r/min離心5 min，將離心后獲得的菌絲體轉(zhuǎn)移到已干燥且確保恒質(zhì)量的稱量瓶中，80 ℃烘干，稱量，計算菌絲體干質(zhì)量。每個試驗設(shè)2個重復(fù)。

1.2.2.2pH值的測定在超凈臺中，無菌操作取15 ml發(fā)酵液于50 ml燒杯中，用pH計進行測定，測定前pH計需要校準(zhǔn)，每個試驗設(shè)2個重復(fù)。

1.2.2.3菌落直徑及菌落生長速率的測定用打孔器定量，將各菌株接種于直徑90 mm的培養(yǎng)皿中，用馬鈴薯葡萄糖瓊脂（PDA）培養(yǎng)基在25 ℃避光培養(yǎng)5 d，測量菌落直徑。菌落生長速率=菌落直徑/培養(yǎng)天數(shù)。每個試驗設(shè)2個重復(fù)。

1.2.2.4菌絲吃料長度及吃料速率的測定將各菌株接種于直徑為12 cm的菌棒中，于25 ℃避光培養(yǎng)13 d，測量吃料長度。吃料速率=吃料長度/培養(yǎng)天數(shù)。每個試驗設(shè)3個重復(fù)。

1.2.3靈芝全基因組重測序分析全基因組重測序試驗流程包括DNA樣品檢測、文庫構(gòu)建、文庫質(zhì)控和上機測序等。

1.2.3.1樣品檢測為保證文庫構(gòu)建質(zhì)量，對基因組DNA進行檢測，待樣品合格后進行文庫構(gòu)建。檢測標(biāo)準(zhǔn)如下：瓊脂糖凝膠電泳檢測結(jié)果顯示基因組DNA主帶完整清晰，且無降解或RNA污染；Nanodrop檢測所得OD260/OD280值為1.8～2.2，無蛋白質(zhì)或肉眼可見的雜質(zhì)污染；Qubit 3.0檢測DNA樣品質(zhì)量濃度大于40 ng/μl，DNA總量大于2 μg。

1.2.3.2文庫的構(gòu)建樣品基因組DNA檢測合格后，嚴格按照NEBNextUltraTMⅡDNA Library Prep Kit for Illumina中提供的標(biāo)準(zhǔn)流程進行文庫構(gòu)建，主要試驗步驟如下：將基因組DNA用bioruptorUCD-200處理成200～500 bp的片段；對片段化的DNA進行末端修復(fù)并加A尾巴，然后連接測序接頭；根據(jù)預(yù)期的文庫選擇特定比例的AMPure XP Beads進行目標(biāo)片段選擇，純化去除接頭污染；通過PCR富集目標(biāo)DNA片段并用AMPure XP Beads進行純化，即完成測序文庫的構(gòu)建。

1.2.3.3文庫質(zhì)控文庫構(gòu)建完成后，對其進行質(zhì)量檢測，檢測結(jié)果達到要求后方可進行上機測序，檢測方法如下：用Qubit 3.0進行初步定量；用Agilent 2100對文庫的插入片段大?。↖nsert size）進行檢測，Insert size符合預(yù)期且無接頭污染才可進行下一步試驗；用德國 ANALYTIKJENA（耶拿）QTOWER實時熒光定量PCR儀對文庫的有效濃度進行準(zhǔn)確定量，即有效濃度＞2 nmol/L為合格文庫。

1.2.3.4上機測序按照目標(biāo)下機數(shù)據(jù)量對文庫進行Pooling（集中），用Illumina HiSeq平臺對DNA分子兩端分別測序150 bp。

1.2.4統(tǒng)計方法用SPSS、Excel 軟件對數(shù)據(jù)進行處理分析，用Duncan’s新復(fù)極差法對組間差異顯著性進行分析，Plt;0.05表示差異具有統(tǒng)計學(xué)意義，用Photoshop軟件處理圖片，用Origin 9.1軟件對所得數(shù)據(jù)作圖，用MEGA 7.0對內(nèi)轉(zhuǎn)錄間隔區(qū)序列構(gòu)建系統(tǒng)發(fā)育樹，用SIMCA 14.1軟件進行主成分分析。

2結(jié)果與分析

2.1野生靈芝菌株的分離及保存

本研究采集的野生靈芝詳見圖1，生境植被為麻櫟、板栗等闊葉樹種，幼嫩靈芝子實體呈近球形或近半圓形，已形成菌蓋的靈芝厚度不一，菌蓋上輪紋明顯，大部分幼嫩子實體為淡黃色，菌蓋邊緣為淡黃色，本研究分離到1株白肉靈芝和1株樹舌靈芝，詳見圖1B、圖1H。

經(jīng)分離獲得靈芝菌株12株（圖2A至圖2L），在馬鈴薯葡萄糖瓊脂（PDA）培養(yǎng)基中，菌落近圓形，菌落周邊菌絲表面平整，菌絲色澤為淡黃色或白色，菌絲均勻、濃密，無氣生菌絲。通過奧林帕斯顯微鏡觀察發(fā)現(xiàn)菌絲分支密度不高（圖2M），且菌絲粗壯（圖2O），有隔膜，可見鎖狀聯(lián)合（圖2N），無雜菌。孢子呈球形（圖2P），直徑為4.51～13.52 μm。將分離獲得的菌株用試管斜面保存于4 ℃冰箱。

2.2靈芝菌株基因組DNA的提取及分子鑒定

使用真菌基因組DNA快速提取試劑盒分別提取分離獲得的12株靈芝菌株基因組DNA，基因組DNA的電泳結(jié)果見圖3A。將提取的基因組DNA低溫保存，用于ITS分子鑒定和全基因組重測序；通過ITS測序，使用鄰接法構(gòu)建系統(tǒng)發(fā)育樹（圖3B）。從進化樹可以看出，12株野生靈芝被聚成3大類，其中第Ⅰ類有4株，占總株數(shù)的33.3%，第Ⅱ類有3株，占總株數(shù)的25.0%，第Ⅲ類有5株，占總株數(shù)的41.7%。由分支情況可知，每個分類的近交程度較低，其中第Ⅲ類表征遺傳多樣性的參數(shù)都高于其他2類（Ⅰ、Ⅱ），表明其遺傳性更豐富。

2.3靈芝菌株生物學(xué)特性的比較分析

對12株靈芝菌株進行生物學(xué)特性測試，由圖4可知，不同菌株之間的菌落直徑、菌落生長速率呈現(xiàn)較明顯的分化，其中菌株lynk002的菌落直徑、菌落生長速率最高，分別為8.39 cm、1.68 cm/d，其次是菌株lynk010，其菌落直徑、菌落生長速率均顯著高于菌株lynk005、lynk006、lynk007和lynk008（P＜0.05），菌株lynk007的菌落直徑、菌落生長速率最低。對不同菌株液體發(fā)酵生物量進行統(tǒng)計發(fā)現(xiàn)，菌株lynk010的菌絲體干質(zhì)量最高，為13.12 g/L，其次為菌株lynk002，菌株lynk003的菌絲體干質(zhì)量最低，菌株lynk010的菌絲體干質(zhì)量顯著高于菌株lynk001、lynk003、lynk004、lynk005、lynk006、lynk007、lynk008、lynk011和lynk012（P＜0.05）。對不同菌株液體發(fā)酵的pH值進行統(tǒng)計發(fā)現(xiàn)，菌株lynk005的pH值最高，為5.29，其次為菌株lynk007，pH值為5.27，菌株lynk002的pH值最低，為4.45，菌株lynk005的pH值顯著高于菌株lynk002、lynk008、lynk009、lynk010和lynk012（P＜0.05）。對不同菌株菌絲的吃料長度、吃料速率進行統(tǒng)計發(fā)現(xiàn)，菌株lynk010菌絲的吃料長度、吃料速率分別為7.27 cm、0.56 cm/d，其次是菌株lynk002，菌株lynk010菌絲的吃料長度、吃料速率顯著高于菌株lynk005、lynk006、lynk007和lynk008菌絲（P＜0.05）。

2.4測序數(shù)據(jù)統(tǒng)計分析

用Illumina HiSeq平臺對有效數(shù)據(jù)進行測序，結(jié)果顯示，高質(zhì)量的有效讀取數(shù)總計為263 448 358條，總堿基數(shù)為79 034 507 400 bp，平均每條測序有效讀取長度為300 bp，有效讀取數(shù)、總堿基數(shù)最多的是菌株lynk009，最少的是菌株lynk008，測序G+C含量為48.04%～58.45%，質(zhì)量數(shù)大于20的堿基所占比例（Q20）、質(zhì)量數(shù)大于30的堿基所占比例（Q30）的平均值分別為94.55%、88.95%，說明測序出錯率低，質(zhì)量符合要求，建庫成功，可用于后續(xù)SNP標(biāo)記挖掘，詳細數(shù)據(jù)統(tǒng)計結(jié)果見表2。

2.5比對率統(tǒng)計分析

以Liu等[19]報道的靈芝單核菌株DH-8的全基因組為參考基因組進行比對，由表3可以看出，在唯一位點比對率方面，菌株lynk002、lynk003和lynk008對應(yīng)的測序數(shù)據(jù)唯一位點比對率較低，其他菌株的唯一位點比對率為20.67%～33.40%，平均唯一位點比對率為22.63%；多位點比對率為0.81%～5.68%，平均比對率為3.68%。本研究采集的靈芝菌株與已報道的靈芝菌株DH-8基因組的比對率較低，說明本研究中的菌株與已報道的菌株之間的遺傳距離較遠。

2.6覆蓋度、覆蓋深度的統(tǒng)計分析

由圖5可以看出，在覆蓋率方面，菌株lynk002、lynk008較低，其他菌株的覆蓋率較接近，覆蓋率為34.23%～36.78%，平均覆蓋率為30.40%。在覆蓋深度方面，菌株lynk002、lynk003的覆蓋深度相對較低，分別為18.80、19.04，菌株lynk009的覆蓋深度最高，為84.49，平均覆蓋深度為50.68。

2.7變異類型分析

單核苷酸多態(tài)性（SNP）的類型及插入缺失標(biāo)記（InDel）長度統(tǒng)計結(jié)果詳見圖6。在本研究中，共獲得了705 742個變異位點，其中598 864個為SNP，106 878個為InDel。SNP數(shù)量統(tǒng)計結(jié)果顯示，401 001個SNP是轉(zhuǎn)換類型（A/G和C/T，為嘌呤之間或嘧啶之間的交換），197 863個SNP是顛換類型（A/C、A/T、C/G和G/T，為嘌呤和嘧啶之間的交換），轉(zhuǎn)換顛換比（Ts/Tv）為2.03。InDel數(shù)量統(tǒng)計結(jié)果顯示，長度小于或等于10 bp的InDel占比達到97.97%，長度大于10 bp的InDel占比為2.03%，長度為2 bp的InDel數(shù)量最多，為35 104個。

2.8變異在基因組上的分布

如表4所示，5條染色體中SNP總數(shù)量為598 864個，在染色體Chr01中，SNP最多，數(shù)量為275 900個，占比為46.07%，同時染色體Chr01對應(yīng)的SNP密度也最高；插入或缺失導(dǎo)致的基因差異在染色體Chr01上表現(xiàn)得較明顯，InDel數(shù)量為47 355個，占總InDel數(shù)量的44.31%，同時對應(yīng)的InDel密度也最高。

2.9主成分分析

為總體分析靈芝菌株之間的遺傳關(guān)系，對測序數(shù)據(jù)進行主成分分析。如圖7所示，貢獻度排前2位的主成分1、主成分2的特征值分別為58.4%、28.2%，合計貢獻率為86.6%。根據(jù)菌株分布情況，大致可以分為3簇（或3類），第Ⅰ簇位于第1象限，共3個菌株，為lynk002、lynk003和lynk008，其中l(wèi)ynk002和lynk008分布得更集中；第Ⅱ簇位于第3象限，共3個菌株，為lynk007、lynk010和lynk012，其中l(wèi)ynk007和lynk010分布得相對較近；第Ⅲ簇位于第2象限和第4象限交匯處，共6個菌株，為lynk001、lynk004、lynk005、lynk006、lynk009和lynk011。

3討論

有研究發(fā)現(xiàn)，除金針菇、黑木耳等少數(shù)單核菌絲可以形成子實體外，多數(shù)食用菌只有雙核菌絲才能形成子實體，雙核菌絲體與單核菌絲體相比，具有分支快、菌絲健壯和穩(wěn)定性好等優(yōu)勢，因此在食用菌生產(chǎn)中，菌株大多選擇雙核菌絲[20]。本研究分離的靈芝菌株具有鎖狀聯(lián)合，而鎖狀聯(lián)合是雙核菌絲的鑒定標(biāo)準(zhǔn)，因此本研究采集的菌株具備形成子實體的條件。菌落生長速率、液體發(fā)酵菌絲體干質(zhì)量和菌絲吃料速率等參數(shù)反映了菌絲的活力，而且呈正相關(guān)，菌株活力也關(guān)系著菌株優(yōu)劣、后續(xù)出菇好壞、栽培時間和子實體產(chǎn)量等[21]。韓鵬等[22]以羊肚菌菌絲生長速度、長勢作為菌株優(yōu)劣的評價標(biāo)準(zhǔn)，董玉蘭等[23-24]分別結(jié)合菌絲萌發(fā)活力、生物量對白靈菇、桑黃菌株的活力進行評價。本研究結(jié)果顯示，在連云港不同區(qū)域采集的野生靈芝種質(zhì)資源，菌株之間的生物學(xué)特性存在較大差異，菌落生長速率、液體發(fā)酵菌絲體干質(zhì)量和菌絲吃料速率最高值分別是最低值的1.49倍、1.77倍和1.40倍，差異程度達到顯著水平（P＜0.05），三者之間呈正相關(guān)趨勢，而與發(fā)酵液的pH值呈負相關(guān)趨勢。

Brown等[25]對32個花生品種進行了全基因組重測序，對品種間基于系譜和基于SNP的遺傳相似性進行對比，發(fā)現(xiàn)基于全基因組重測序數(shù)據(jù)獲得的結(jié)果更加準(zhǔn)確，具有高度相關(guān)性的品種緊密地聚在一起，可見重測序在探索品種的遺傳多樣性、鑒定未知來源材料方面具有不可比擬的優(yōu)越性。郭敏杰等[26]發(fā)現(xiàn)，基因型數(shù)據(jù)更能準(zhǔn)確地反映品種內(nèi)在遺傳基礎(chǔ)。本研究通過ITS聚類分析，將12株野生靈芝種質(zhì)資源聚成3大類，通過主成分分析，根據(jù)菌株在坐標(biāo)上的分布情況，同樣將這12個菌株分為3類，其中聚類分析的第2類與主成分分析的第Ⅱ簇完全吻合，而聚類分析的第Ⅰ、Ⅲ類與主成分分析的第Ⅰ、Ⅲ簇存在一定出入，聚類分析中第Ⅰ類中的lynk006、lynk009在主成分分析中聚為第Ⅲ簇，聚類分析中第Ⅲ類中的lynk002在主成分分析中聚為第Ⅰ簇，由于主成分分析數(shù)據(jù)基于全基因組重測序數(shù)據(jù)，該方法更可靠，上述2種分析方法的第Ⅲ類表征遺傳多樣性的參數(shù)都高于其他2類（Ⅰ、Ⅱ），說明第Ⅲ類具有更豐富的遺傳多樣性。

目前基于二代測序的全基因組重測序技術(shù)已作為主流方法應(yīng)用于很多物種SNP標(biāo)記的開發(fā)[27]。Shen等[28]揭示了四倍體棉花的SNP和InDel多樣性，Zhang等[29]確定了辣椒第1花節(jié)的候選基因，但靈芝這方面的研究還未見相關(guān)報道。本研究所測序列的Q20、Q30占比較高，說明測序出錯率很低，GC分布正常，數(shù)據(jù)量達到預(yù)期目標(biāo)，本研究采集的靈芝菌株與已報道的靈芝菌株DH-8的基因組比對率較低，說明本研究中的菌株與已報道的菌株之間遺傳距離較遠。12株靈芝中的SNP數(shù)量為598 864個，InDel數(shù)量為106 878個，在第1條染色體中SNP最多，占比46.07%，SNP中轉(zhuǎn)換顛換比（Ts/Tv）為 2.03，說明變異以轉(zhuǎn)換為主，插入或缺失導(dǎo)致的基因差異在第1條染色體上表現(xiàn)得較為明顯，占總InDel數(shù)的44.31%，表明可能有更多的功能基因在第1條染色體上，第1條染色體可以用于進一步挖掘更多功能基因。

4結(jié)論

本研究對連云港范圍內(nèi)野生靈芝生物學(xué)特性進行分析，使用全基因組重測序技術(shù)對靈芝基因組進行測序，并進行序列特征分析和主成分分析。結(jié)果表明，從不同區(qū)域采集的野生靈芝種質(zhì)資源在菌落生長速率、液體發(fā)酵液生物量和菌絲吃料速率等生物特性方面展現(xiàn)出較大差異；通過主成分分析將12株野生靈芝種質(zhì)資源聚為3類，其中第Ⅲ類具有更豐富的遺傳多樣性。在12株靈芝基因組中，第1條染色體的SNP最多，變異以轉(zhuǎn)換為主，插入或缺失導(dǎo)致的基因差異在第1條染色體上表現(xiàn)得明顯，說明可能有更多功能基因在第1條染色體上，第1條染色體可以用于進一步挖掘更多的功能基因。本研究結(jié)果可為后期分子標(biāo)注技術(shù)與靈芝生物學(xué)特性、表型性狀評價相結(jié)合，創(chuàng)造優(yōu)質(zhì)種質(zhì)資源和挖掘優(yōu)質(zhì)基因提供科學(xué)的理論基礎(chǔ)。

參考文獻：

[1]胡惠萍，劉遠超，莫偉鵬，等. 兩株西藏白肉靈芝菌株特性初探[J]. 食用菌學(xué)報，2017，24（1）：50-54.

[2]戴玉成，曹云，周麗偉，等. 中國靈芝學(xué)名之管見[J]. 菌物學(xué)報，2013，32（6）：947-952.

[3]SOHRETOGLU D， HUANG S. Ganoderma lingzhi polysaccharides as an anti-cancer agent[J]. Anti-Cancer Agents in Medicinal Chemistry，2018，18（5）：667-674.

[4]王朝川. 靈芝成分及功能的研究現(xiàn)狀[J]. 中國果菜，2018，38（8）：45-47，53.

[5]MENG J， WANG S Z， HE J Z， et al. Ganoderic acid A is the effective ingredient of Ganoderma triterpenes in retarding renal cyst development in polycystic kidney disease[J]. International Journal of Sports Medicine，2020，41（6）：782-790.

[6]李丹妮，朱長俊，鄭雨晴，等. 靈芝種植及初加工研究進展[J]. 南方農(nóng)業(yè)，2021，15（8）：233-235.

[7]崔黎明，劉利，潘麗艷，等. 靈芝孢子誤診為寄生蟲蟲卵3例分析[J]. 中國實驗診斷學(xué)，2016，20（7）：1191-1192.

[8]蒙田秀，楊力龍，龔志強，等. 壯族藥白鶴靈芝萘醌類化學(xué)成分及其藥理作用研究進展[J]. 中國實驗方劑學(xué)雜志，2020，26（10）：213-219.

[9]王夢晨，張雪漣，陳向東，等. 靈芝三萜與靈芝多糖抗腫瘤作用及其機制研究進展[J]. 中國實驗方劑學(xué)雜志，2022，28（5）：234-241.

[10]CHEN L H， ABUDUMIJITI A， LI M. Protective effect of Ganoderma （Lingzhi） on radiation and chemotherapy[J]. Advances in Experimental Medicine and Biology，2019，1182：119-142.

[11]REN L， ZHANG J， ZHANG T H. Immunomodulatory activities of polysaccharides from Ganoderma on immune effector cells[J]. Food Chemistry，2021，340：127933.

[12]莊靜. 試論中藥靈芝補肺湯在支氣管哮喘病例治療中的臨床應(yīng)用效果觀察[J]. 中國醫(yī)藥指南，2017，15（21）：183-184.

[13]沈滟惠，溫亞格，唐慶九，等. 靈芝抗腫瘤活性及其機制研究進展[J]. 海峽藥學(xué)，2020，32（4）：8-12.

[14]蘭慧瑩，毛可紅，陳俏彪，等. 工藝靈芝的制作技術(shù)芻議[J]. 食藥用菌，2019，27（6）：427-429.

[15]唐才明，張燕，陳嬌嬌，等. CO2緩釋劑在盆栽靈芝中應(yīng)用初探[J]. 食用菌，2019，41（2）：55-56， 65.

[16]豆芬芬. 連云港港口與城市耦合協(xié)調(diào)發(fā)展研究[D]. 徐州：中國礦業(yè)大學(xué)，2014.

[17]林祥天，邵世光，呂繼江，等. 連云港地區(qū)食用菌和毒菌的研究[J]. 江蘇預(yù)防醫(yī)學(xué)，1997（4）：13-14.

[18]邢煒，洪豆，徐碧林，等. 黃岡大別山一株野生靈芝的鑒定與馴化栽培[J]. 中國食用菌，2021，40（10）：28-33.

[19]LIU Y C， HUANG L H， HU H P， et al. Whole-genome assembly of Ganoderma leucocontextum （Ganodermataceae， Fungi） discovered from the Tibetan Plateau of China[J]. G3 Genes|Genomes|Genetics，2021，11（12）：jkab337.

[20]鮑大鵬. 食用菌雜交育種中的科學(xué)問題[J]. 食用菌學(xué)報，2020，27（4）：1-24.

[21]張先吉，花紀(jì)，顏俊清，等. 幾株主栽靈芝菌株液體菌種培養(yǎng)基優(yōu)化及其活力評價[J]. 食品科技，2020，45（4）：19-25.

[22]韓鵬，陳青君，賀國強，等. 2種野生羊肚菌分離、鑒定與菌絲培養(yǎng)條件[J]. 北京農(nóng)學(xué)院學(xué)報，2018，33（1）：37-42.

[23]董玉蘭，李書生，張麗萍，等. 白靈菇高活力液體菌種搖瓶發(fā)酵培養(yǎng)條件的優(yōu)化[J]. 北方園藝，2014（6）：143-146.

[24]曲德輝，楊焱，張赫男，等. 桑黃菌株活力評價及優(yōu)良菌株篩選[J]. 菌物學(xué)報，2016，35（10）：1187-1198.

[25]BROWN N， BRANCH W D， JOHNSON M， et al. Genetic diversity assessment of Georgia peanut cultivars developed during ninety years of breeding[J]. The Plant Genome，2021，14（3）：e20141.

[26]郭敏杰，鄧麗，李玉榮，等. 基于主成分和聚類分析的冀花高油酸花生品種綜合評價[J]. 中國油料作物學(xué)報，2022，44（6）：1210-1217.

[27]OGDEN R， GHARBI K， MUGUE N， et al. Sturgeon conservation genomics：SNP discovery and validation using RAD sequencing[J]. Molecular Ecology，2013，22（11）：3112-3123．

[28]SHEN C， JIN X， ZHU D， et al. Uncovering SNP and indel variations of tetraploid cottons by SLAF-seq[J]. BMC Genomics，2017，18：247.

[29]ZHANG X F， WANG G Y， CHEN B， et al. Candidate genes for first flower node identified in pepper using combined SLAF-seq and BSA[J]. PLoS One，2018，13（3）：e0194071.

（責(zé)任編輯：徐艷）