摘要：以生產(chǎn)上收集的一個(gè)平菇白色變異菌株（Po50）、黑色出發(fā)菌株（Po51）和一個(gè)白色生產(chǎn)菌株（Po9）為材料，采用拮抗、酯酶同工酶和全基因組重測(cè)序方法對(duì)3個(gè)菌株進(jìn)行鑒定。拮抗和酯酶同工酶結(jié)果表明，白色變異菌株和黑色出發(fā)菌株之間無(wú)明顯拮抗和酶譜差異，與生產(chǎn)平菇白色菌株差異較大；進(jìn)一步通過(guò)全基因組重測(cè)序技術(shù)鑒定結(jié)果表明，3個(gè)菌株鑒定出大量的SNP（單核苷酸多態(tài)性位點(diǎn)）、InDel（插入缺失位點(diǎn)）、SV（結(jié)構(gòu)變異位點(diǎn)）。平菇白色變異菌株P(guān)o50檢測(cè)到SNP、InDel突變總數(shù)為1 504 391個(gè)，生產(chǎn)菌株P(guān)o9檢測(cè)到SNP、InDel突變總數(shù)為1 371 818個(gè)，黑色出發(fā)菌株P(guān)o51檢測(cè)到SNP、InDel突變總數(shù)為1 501 877個(gè)，通過(guò)生物信息手段分析白色變異菌株黑色出發(fā)菌株以及對(duì)照菌株P(guān)o9基因組間的結(jié)構(gòu)差異，獲得遺傳變異圖譜，可以對(duì)變異菌株進(jìn)行快速精準(zhǔn)鑒定。

關(guān)鍵詞：平菇；變異菌株；全基因組重測(cè)序技術(shù)；酯酶同工酶

中圖分類號(hào)：S646.1⁺40.1 文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1002-1302（2024）18-0041-09

收稿日期：2023-09-20

基金項(xiàng)目：河北省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系食用菌創(chuàng)新團(tuán)隊(duì)建設(shè)專項(xiàng)資金（編號(hào)：HBCT2023090207）；河北省農(nóng)業(yè)科技成果轉(zhuǎn)化項(xiàng)目（編號(hào)：202360101010014）。

作者簡(jiǎn)介：夏會(huì)楠（1998—），女，河北承德人，碩士，從事食用和藥用真菌研究。E-mail：2199797807@qq.com。

通信作者：鄭素月，博士，教授，碩士生導(dǎo)師，主要從事食用菌遺傳育種研究。E-mail：zhengsuyue@sina.com。

平菇（Pleurotus ostreatus）栽培廣泛，含有人體所需的大量生物活性物質(zhì)、礦物質(zhì)和維生素，具有很高的營(yíng)養(yǎng)價(jià)值。由于其栽培技術(shù)簡(jiǎn)單、生長(zhǎng)周期短、收益快、適應(yīng)性強(qiáng)，為食用菌生產(chǎn)主栽品種之一。平菇產(chǎn)業(yè)中存在的主要問(wèn)題是育種工作落后，品種單一，生產(chǎn)缺乏優(yōu)良品種，因此平菇新品種種質(zhì)資源選育十分重要，經(jīng)常規(guī)鑒定后在分子層面進(jìn)行精準(zhǔn)分析，有利于更方便快捷地選育新品種。聶興華等以野生板栗為試驗(yàn)材料，通過(guò)重測(cè)序技術(shù)確定46份野生板栗的分類，發(fā)現(xiàn)野生板栗和栽培板栗屬于同種，為我國(guó)野生板栗提供了更準(zhǔn)確的邊界范圍^［1］。沈秀芬等通過(guò)對(duì)5個(gè)香菇進(jìn)行重測(cè)序，分析這些菌株中的插入/缺失位點(diǎn)開發(fā)InDel標(biāo)記，可將44個(gè)香菇菌株分為4個(gè)亞群^［2］。侯炳豪等以鐵觀音茶樹為試驗(yàn)材料，通過(guò)重測(cè)序技術(shù)對(duì)鐵觀音無(wú)性繁殖后代進(jìn)行遺傳差異研究，根據(jù)重測(cè)序數(shù)據(jù)得出在全基因組和編碼區(qū)范圍內(nèi)，各個(gè)樣本之間總變異位點(diǎn)數(shù)均差別不大，各樣本中涉及突變的基因數(shù)基本一致，推測(cè)鐵觀音品種的種性在無(wú)性繁殖過(guò)程中總體保持穩(wěn)定^［3］。

全基因組重測(cè)序技術(shù)對(duì)已知基因序列的物種進(jìn)行不同個(gè)體間的基因組重測(cè)序，并在此基礎(chǔ)上對(duì)個(gè)體進(jìn)行差異性分析^［4］，通過(guò)序列對(duì)比，在平菇基因組水平上開發(fā)大量分子標(biāo)記如鑒定食用菌中大量的SNP單核苷酸多態(tài)性位點(diǎn)、InDel插入缺失位點(diǎn)、SV結(jié)構(gòu)變異位點(diǎn)^［5-6］，通過(guò)生物信息手段，分析平菇白色變異菌株、栽培菌株P(guān)o9以及黑色出發(fā)菌株基因組間的結(jié)構(gòu)差異，獲得遺傳變異圖譜。全基因組重測(cè)序技術(shù)有助于快速發(fā)現(xiàn)重要性狀的遺傳變異以及對(duì)變異菌株精準(zhǔn)鑒定^［7-8］。

1 材料與方法

1.1 試驗(yàn)材料

試驗(yàn)于2023年2—5月在河北工程大學(xué)食用菌研究室內(nèi)進(jìn)行，在平菇系統(tǒng)選育過(guò)程中，發(fā)現(xiàn)變異白色平菇菌株1個(gè)，變異的白色子實(shí)體和對(duì)應(yīng)的出發(fā)菌株形態(tài)見圖1。對(duì)采集的白色變異菌株子實(shí)體和黑色出發(fā)菌株子實(shí)體進(jìn)行組織分離，獲得純菌種，白色變異菌株編號(hào)為Po50，出發(fā)黑色菌株編號(hào)為Po51。生產(chǎn)上栽培的白色平菇編號(hào)為Po9，以上材料均由河北工程大學(xué)食用菌實(shí)驗(yàn)室分離保藏。

1.2 試驗(yàn)方法

1.2.1 拮抗和酯酶同工酶測(cè)定 以平菇89、平菇99、雙抗黑平以及生長(zhǎng)上栽培的白平菇Po9作為對(duì)照菌株，對(duì)白色變異菌株和黑色出發(fā)菌株進(jìn)行拮抗和酯酶同工酶測(cè)定^［9］。

1.2.2 全基因組重測(cè)序技術(shù)

對(duì)平菇白色變異菌株P(guān)o50、黑色出發(fā)菌株P(guān)o51和生產(chǎn)對(duì)照白色菌株P(guān)o9進(jìn)一步進(jìn)行基因組重測(cè)序技術(shù)，文庫(kù)構(gòu)建和重測(cè)序數(shù)據(jù)均由北京諾禾致源科技股份有限公司提供。

1.2.3 SNP/InDel檢測(cè)及注釋統(tǒng)計(jì)

采用GATK 4.0.9.0軟件進(jìn)行MarkDuplicates、BaseRecalibrator處理后，再用GATK軟件中的HaplotypeCaller進(jìn)行SNP、InDel檢測(cè)，并對(duì)變異位點(diǎn)進(jìn)行質(zhì)控標(biāo)記。

2 結(jié)果與分析

2.1 白色變異菌株親緣關(guān)系鑒定

拮抗測(cè)定結(jié)果表明，平菇白色菌株P(guān)o50與出發(fā)黑色菌株P(guān)o51無(wú)拮抗作用，二者均與平菇99無(wú)拮抗作用（圖2）。酯酶同工酶圖譜顯示平菇白色變異菌株與原黑色菌株及標(biāo)準(zhǔn)99菌株酯酶同工酶圖譜相同，與2個(gè)標(biāo)準(zhǔn)菌株（平菇89、雙抗黑平、平菇黑色）和生產(chǎn)栽培白平菇Po9有較大差異（圖3），可以初步確定出發(fā)黑色菌株為平菇99的同物異名，而白色變異菌株在遺傳上與黑色出發(fā)菌株無(wú)顯著差異。

2.2 全基因組重測(cè)序技術(shù)

2.2.1 測(cè)序數(shù)據(jù)質(zhì)量評(píng)估及對(duì)比

利用Illumina PE150測(cè)序，由于原始測(cè)序數(shù)據(jù)可能包含低質(zhì)量序列、接頭序列等，為了保證信息分析結(jié)果的可靠性，需要過(guò)濾這些雜質(zhì)，從而得到高質(zhì)量clean data，共獲得2 020 Mb的平菇白色變異菌株、2 019 Mb的平菇黑色出發(fā)菌株、2 967 Mb的Po9基因組原始數(shù)據(jù)（raw data）。經(jīng)質(zhì)量評(píng)估和過(guò)濾雜質(zhì)，最終平菇白色變異菌株獲得1 763 Mb高質(zhì)量數(shù)據(jù) clean data、平菇黑色出發(fā)菌株獲得1 779 Mb高質(zhì)量數(shù)據(jù)clean data、Po9獲得2 519 Mb高質(zhì)量數(shù)據(jù)clean data，采用測(cè)序平臺(tái)對(duì)樣品進(jìn)行測(cè)序，數(shù)據(jù)統(tǒng)計(jì)見表1。

參考基因組對(duì)比結(jié)果將有效的測(cè)序數(shù)據(jù)通過(guò)BWA對(duì)比參考基因組得到最初對(duì)比結(jié)果，利用最初對(duì)比結(jié)果進(jìn)行mapping率、duplication等的統(tǒng)計(jì)。重測(cè)序數(shù)據(jù)與參考數(shù)據(jù)樣品對(duì)比平菇白色變異菌株P(guān)o50檢測(cè)到的reads數(shù)目為11 754 734個(gè)，黑色平菇Po51共檢測(cè)到reads數(shù)目為11 862 522個(gè)，Po9菌株共檢測(cè)到16 799 756個(gè) reads。整體對(duì)比率分別為80.48、81.11、75.83，覆蓋率分別為90.33%、90.32%、83.08%（表2）。

2.2.2 SNP及InDel突變檢測(cè)

隨著食用菌基因組的公布，通過(guò)全基因組重測(cè)序以及生物技術(shù)挖掘平菇基因組SNP標(biāo)記簡(jiǎn)便快捷，SNP主要是指在基因組水平上由單個(gè)核苷酸的變異所引起的DNA序列多態(tài)性。SNP所表現(xiàn)的多態(tài)性只涉及到單個(gè)堿基的變異，這種變異可由單個(gè)堿基的轉(zhuǎn)換（transition）或顛換（transversion）所引起。將所測(cè)得的數(shù)據(jù)采用GATK軟件進(jìn)行MarkDuplicates、BaseRecalibrator處理后，在利用GATK軟件中的HaplotypeCaller進(jìn)行SNP檢測(cè)，最后在對(duì)檢測(cè)的結(jié)果進(jìn)行質(zhì)量過(guò)濾。最終得到SNP、InDel統(tǒng)計(jì)結(jié)果（表3），分別檢測(cè)到突變總數(shù)為1 504 391、1 501 877、1 371 818個(gè)。

2.2.3 SNP注釋

采用使用ANNOVAR對(duì)檢測(cè)到的SNP進(jìn)行注釋，可對(duì)變異組合進(jìn)行多層次的組合篩選。SNP在各染色體上分布統(tǒng)計(jì)見圖4，Po50在染色體3503156.1上SNP數(shù)目最多，其次是染色體3503157.1、3503160.1、3503158.1，在染色體3503167.1、3503168.1、3503169.1、3503170.1、3503171.1、009905.1上SNP數(shù)目最少。Po9在染色體3503156.1、3503157.1上 SNP數(shù)目最多，在染色體3503167.1、3503168.1、3503169.1、3503170.1、3503171.1、009905.1上SNP數(shù)目最少。通過(guò)突變頻譜分析可以直觀看出點(diǎn)突變包含6種類型：T∶A→G∶C，T∶A→C∶G，T∶A→A∶T，C∶G→T∶A，C∶G→G∶C 和 C∶G→A∶T，在Po51、Po50、Po9中各種突變類型的比例存在某種突變類型的偏好性，其中T∶A→C∶G、C∶G→A∶T類型的突變較多，突變頻譜分布見圖5。

根據(jù)ANNOVAR的注釋結(jié)果，SNP在基因組各區(qū)域分布統(tǒng)計(jì)見圖6，SNP變異主要集中于發(fā)生在外顯子區(qū)域，比_{e71d63d306b1d50358201d0558687f6d4af4b42aab7079ad363ff0e68f703025}例在45.33%～45.66%，發(fā)生在基因下游的SNP變異比例在14.51%～14.61%，基因上游、基因區(qū)間、內(nèi)含子區(qū)域、ncRNA區(qū)域、基因剪切區(qū)域比例分別占16.86%、5.57%～5.61%、17.23%～17.52%、0、0.07%。其中，在外顯子區(qū)域發(fā)生同義突變的SNP比例為71.11%～71.24%，非同義突變?yōu)?8.52%～28.64%，使基因轉(zhuǎn)錄提前終止突變?yōu)?.17%，失去終止密碼子突變?yōu)?.03%（圖7）。

2.2.4 InDel注釋及統(tǒng)計(jì)

利用 ANNOVAR 對(duì)檢測(cè)出的 InDel 進(jìn)行注釋，InDel突變檢測(cè)與注釋方法與SNP檢測(cè)注釋方法一致。 InDel在基因組各區(qū)域分布統(tǒng)計(jì)見圖8，InDel變異主要集中于發(fā)生在基因上游外顯子區(qū)域，比例在31.90%～32.07%；發(fā)生在基因下游的 InDel變異比例在27.29%～27.49%，基因區(qū)間、內(nèi)含子區(qū)域、ncRNA區(qū)域、基因剪切區(qū)域比例分別占8.27%～8.28%、21.48%～21.49%、0.02%、0.39%?？傮w來(lái)說(shuō)，InDel突變?cè)?/p>

全基因組范圍內(nèi)的表現(xiàn)與SNP突變基本一致。在外顯子區(qū)域內(nèi)，移碼缺失突變所占比例最高，占33.98%～34.20%，移碼插入突變占外顯子區(qū)域突變總數(shù)的27.48%～28.24%，非移碼插入突變占17.56%～18.63%（圖9）。

在基因外顯子區(qū)域的SNP（InDel）位點(diǎn)中，部分位點(diǎn)導(dǎo)致基因編碼中氨基酸改變，而引起基因產(chǎn)物的突變，可能會(huì)影響平菇基因的生物學(xué)功能。

3 平菇差異基因的GO功能富集分析

使用BLAST2GO v2.5對(duì)所測(cè)Po50、Po9基因進(jìn)行GO功能注釋（圖10），Po50共注釋到5 760個(gè)基因，注釋到37個(gè)類別，主要包括生物過(guò)程（biological process，BP）、細(xì)胞組成（cellular component，CC）、分子功能（molecular function，MF）3個(gè)類別，其中CC數(shù)量最多，注釋到12個(gè)條目，總數(shù)為2 112個(gè)，最多的集中在細(xì)胞膜，占總數(shù)的24.5%，數(shù)量達(dá)到517個(gè)，其次是細(xì)胞和細(xì)胞部分（GO：0009696），分別占總數(shù)的13.4%和13.1%，數(shù)量為283個(gè)及276個(gè)；MF注釋8條，共有1 466個(gè)基因，最多的是催化活性部分（GO：0000155），占總數(shù)的51.9%，數(shù)量為761個(gè)，黏合物（GO：0046872）較少于活性催化部分，占總數(shù)的38.6%，數(shù)量為566個(gè)。其次是運(yùn)轉(zhuǎn)活性、結(jié)構(gòu)分子活性、活性轉(zhuǎn)錄因子，分別占總數(shù)的4.6%、2.0%、1.4%；氧化活性、分子轉(zhuǎn)換器活性最少，分別占總數(shù)的0.3%、0.1%。

Po9共注釋到5 386個(gè)基因，注釋到38個(gè)類別，主要包括生物過(guò)程、細(xì)胞組成、分子功能3個(gè)類別，其中CC數(shù)量最多，注釋到12個(gè)條目，總數(shù)為1 979個(gè)，占3個(gè)類別總數(shù)的43.8%。最多的集中在細(xì)胞膜（GO：0016021），占總數(shù)的25.0%，數(shù)量達(dá)到494個(gè)：細(xì)胞膜組成部分（GO：000021）少于細(xì)胞膜，占總數(shù)的24.0%，數(shù)量為474個(gè)，其次是細(xì)胞、細(xì)胞部分。MF注釋8條，共有1 392個(gè)基因，最多的是催化活性部分（GO：0000155），占總數(shù)的52.1%，數(shù)量為725個(gè)，黏合物（GO：0046872）較少于活性催化部分，占總數(shù)的38.1%，數(shù)量為530個(gè)。其次是轉(zhuǎn)錄活性，占總數(shù)的4.5%，數(shù)量為63個(gè)。在生物過(guò)程中Po9較Po50有1個(gè)基因注釋到細(xì)胞聚集類別中，生物過(guò)程及分子功能Po9涉及的基因數(shù)量，與Po50基本一致。

分別篩選出Po50、Po9的差異基因，將差異基因向GO數(shù)據(jù)庫(kù)映射，進(jìn)行分析，主要富集在生物過(guò)程、分子功能和細(xì)胞組成相關(guān)的3個(gè) GO 類別中，Po50、Po9在生物過(guò)程中參與細(xì)胞過(guò)程上調(diào)表達(dá)基因分別有372個(gè)和349個(gè)，代謝過(guò)程上調(diào)表達(dá)基因分別有357個(gè)和376個(gè)；在細(xì)胞組分過(guò)程中細(xì)胞膜上調(diào)表達(dá)基因分別有517個(gè)和494個(gè)，細(xì)胞膜組成過(guò)程中上調(diào)表達(dá)基因分別有498個(gè)和474個(gè)；在分子功能中催化活性上調(diào)表達(dá)基因分別有761個(gè)和725個(gè)，黏合物上調(diào)表達(dá)基因分別有566個(gè)和494個(gè)（圖11）。

3.1 不同差異基因indel網(wǎng)狀圖

Po50選擇10個(gè)顯著的代謝通路如圖12所示：肌動(dòng)蛋白細(xì)胞骨架調(diào)控（regulation of actin cytoskeleton）由PC9H_000055等102個(gè)基因控制；志賀菌病的致病過(guò)程（shigellosis）由PC9H_000055等90個(gè)基因控制；致病性大腸桿菌感染（pathogenic escherichia colii nfection）由PC9H_000055等92個(gè)差異基因控制；錯(cuò)配修復(fù)（mismatch repair）由PC9H__{bMkmRpjPmsYcWVy01BjqIt0/ln7g3cTLk74CKl8Qbm8=}000717等21個(gè)差異基因控制。附著連接（adherens junction）由PC9H_000055等67個(gè)差異基因控制；表皮細(xì)胞細(xì)菌入侵（bacterial invasion of epithelial cells）由PC9H_000055等71個(gè)差異基因控制；癌癥中的膽堿代謝（choline metabolism in cancer）由PC9H_000055等76個(gè)差異基因控制，其中還有緊密連接（tight junction）、沙門氏菌感染（salmonella infection）病毒性心肌炎（viral myocarditis）是個(gè)顯著通路關(guān)系網(wǎng)狀圖。

Po9選擇10個(gè)顯著的代謝通路：沙門氏菌感染由PC9H_000055等68個(gè)差異基因控制；肌動(dòng)蛋白細(xì)胞骨架調(diào)控由PC9H_000055等101個(gè)差異基因控制；志賀菌病的致病過(guò)程由PC9H_000055等84個(gè)基因控制；附著連接由PC9H_000055等65個(gè)差異基因控制；病毒性心肌炎由PC9H_000271等30個(gè)差異基因控制；癌癥中的膽堿代謝由PC9H_000055等71個(gè)差異基因控制；致病性大腸桿菌感染，由PC9H_000249等89個(gè)差異基因控制；表皮細(xì)胞細(xì)菌入侵由PC9H_000249等67個(gè)差異基因控制；趨化因子信號(hào)通路由PC9H_000249等67個(gè)差異基因控制。

3.2 不同糙皮側(cè)耳差異基因的KEGG代謝通路富集分析

將糙皮側(cè)耳差異基因通過(guò)KEGG數(shù)據(jù)庫(kù)進(jìn)行Pathway分析，選取KEGG數(shù)據(jù)庫(kù)里的真菌類，選取顯著的通路制作散點(diǎn)圖。對(duì)所得到的轉(zhuǎn)錄本進(jìn)行驗(yàn)證和注釋并篩選出的差異基因進(jìn)行富集通路注釋。將生物代謝通路分為6個(gè)類別：細(xì)胞過(guò)程（cellular processes）、生物體系統(tǒng)（organismal systems）、環(huán)境信息處理（environmental information processing）、遺傳信息處理（genetic information processing）、人類疾?。╤uman diseases）、新陳代謝（metabolism）。這些差異基因被分配到7個(gè)第1層級(jí)通路途徑第2層44個(gè)KEGG通路途徑，其中差異基因大部分富集在新陳代謝、疾病、生物體系統(tǒng)以及細(xì)胞過(guò)程上，其次是生物系統(tǒng)以及環(huán)境信息處理上，其中氨基酸和核苷酸糖代謝、緊密連接、白細(xì)胞骨架調(diào)節(jié)和附著連接在新陳代謝、細(xì)胞過(guò)程占主導(dǎo)地位（圖13）。

4 討論與結(jié)論

全基因組重測(cè)序已經(jīng)廣泛應(yīng)用到人類、動(dòng)物、植物的基因組和轉(zhuǎn)錄組分析中^［9-12］，然而全基因組重測(cè)序技術(shù)在食用菌中應(yīng)用較少，通常是應(yīng)用重測(cè)序技術(shù)進(jìn)行ISSR分子標(biāo)記的開發(fā)，在SNP以及InDel標(biāo)記研究較少^［13-15］。

本試驗(yàn)通過(guò)對(duì)白色變異菌株進(jìn)行全基因組重測(cè)序從基因組水平精準(zhǔn)快速地挖掘了大量的單核苷酸多態(tài)性位點(diǎn)SNP，插入缺失位點(diǎn)InDel分析不同個(gè)體間的結(jié)構(gòu)差異。Illumina PE150測(cè)序平菇白色變異菌株獲得 1 763 Mb高質(zhì)量數(shù)據(jù)clean data、平菇黑色出發(fā)菌株獲得 1 779 Mb高質(zhì)量數(shù)據(jù) clean data、Po9獲得2 519 Mb高質(zhì)量數(shù)據(jù)。在進(jìn)行全基因組測(cè)序時(shí)還會(huì)檢測(cè)到低質(zhì)量的數(shù)據(jù)，這些低質(zhì)量數(shù)據(jù)會(huì)影響基因組的分析，因此基因型對(duì)于基因組數(shù)據(jù)開發(fā)SNP以及InDel十分重要，所以重測(cè)序數(shù)據(jù)質(zhì)量評(píng)估至關(guān)重要。重測(cè)序數(shù)據(jù)與參考數(shù)據(jù)樣品對(duì)比平菇白色變異菌株P(guān)o50檢測(cè)到的reads數(shù)目為11 754 734個(gè)，黑色平菇Po51共檢測(cè)到reads數(shù)目為11 862 522個(gè)，Po9菌株共檢測(cè)到16 799 756個(gè) reads。分別檢測(cè)到Po51 SNP突變總數(shù)1 504 391、1 501 877、1 371 818個(gè)。Po50在染色體503 156.1上SNP數(shù)目最多。通過(guò)突變頻譜分析發(fā)現(xiàn)其中T∶A→C∶G、C∶G→A∶T類型的突變較多，SNP變異主要集中于發(fā)生在外顯子區(qū)域。InDel變異主要集中于發(fā)生在基因上游外顯子區(qū)域。Po50、Po9測(cè)基因進(jìn)行GO功能注釋，Po50共注釋到5 760個(gè)基因，注釋到38個(gè)類別，其中CC數(shù)量最多，注釋到12個(gè)條目，總數(shù)為2 433個(gè)。將糙皮側(cè)耳差異基因通過(guò)KEGG數(shù)據(jù)庫(kù)進(jìn)行Pathway分析，可將生物代謝通路分為6個(gè)類別：細(xì)胞過(guò)程、生物體系統(tǒng)、環(huán)境信息處理、遺傳信息處理、人類疾病、新陳代謝。這些差異基因被分配到7個(gè)第1層級(jí)通路途徑第2層44個(gè)KEGG通路途徑^［16-18］。本試驗(yàn)通過(guò)重測(cè)序技術(shù)對(duì)平菇白色變異菌株進(jìn)行鑒定，然而定位到導(dǎo)致顏色變異的基因還有待進(jìn)一步研究。

參考文獻(xiàn)：

［1］聶興華，張 煜，劉 松，等. 基于基因組重測(cè)序的野生板栗遺傳特征和分類地位研究［J］. 園藝學(xué)報(bào)，2023，50（8）：1622-1636.

［2］沈秀芬，章爐軍，張美彥，等. 利用InDel標(biāo)記分析中國(guó)香菇菌株的遺傳多樣性與群體結(jié)構(gòu)［J］. 菌物學(xué)報(bào)，2021，40（9）：2266-2281.

［3］侯炳豪，高 婷，魏月德，等. 基于高深度基因組重測(cè)序的‘鐵觀音’茶樹無(wú)性繁殖后代遺傳變異研究［J］. 園藝學(xué)報(bào)，2023，50（7）：1505-1517.

［4］Zhao W，F(xiàn)an L，Wu W J，et al. Re-sequencing and transcriptomic analysis reveal differential expression patterns and sequence variation in glucosyltransferase gene related to anthocyanin biosynthesis in walnut （Juglans regia L.）［J］. Scientia Horticulturae，2023，317：112077.

［5］Zhang Q P，Zhang Y P，Liu W S，et al. Re-sequencing and morphological data revealed the genetics of stone shell and kernel traits in apricot［J］. Frontiers in Plant Science，2023，14：1196754.

［6］Kim J Y，Hwang J E，Eo S H，et al. Development of InDel markers for interspecific hybridization between hill pigeons and feral pigeons based on whole-genome re-sequencing［J］. Scientific Reports，2022，12（1）：22618.

［7］Shan T F，Li Y Q，Pang S J. Identification of a genomic region linked with sex determination of Undaria pinnatifida （Alariaceae） through genomic resequencing and genetic linkage analyses of a segregating gametophyte family［J］. Journal of Phycology，2023，59（1）：193-203.

［8］Cheng Q，Sun L，Qiao H，et al. Loci underlying leaf agronomic traits identified by re-sequencing celery accessions based on an assembled genome［J］. iScience，2022，25（7）：104565.

［9］劉曉雪，王 強(qiáng)，張彬彬，等. 基于酯酶和ISSR技術(shù)的平菇單核菌株遺傳多樣性分析［J］. 江蘇農(nóng)業(yè)科學(xué)，2022，50（9）：27-32.

［10］ShresthaS，F(xiàn)u Y Q，Michael V N，et al.ftXyIQGdctfetEbRlIEYW7VWtLLRoIhM1OdUDfBZkTI= Whole genome re-sequencing and bulk segregant analysis reveals chromosomal location for Papaya ringspot virus W resistance in squash［J］. Frontiers in Plant Science，2022，13：848631.

［11］Zhao W，Zhang Y P，Zhang J P，et al. QTL mapping by whole genome re-sequencing and analysis of candidate genes for salt tolerance in linseed （Linum usitatissmum L.）［J］. Oil Crop Science，2022，7（2）：80-85.

［12］刁興旺，吳莉君，何 紅，等. 芒果炭疽病抗感品種全基因組重測(cè)序分析［J］. 江蘇農(nóng)業(yè)科學(xué)，2022，50（23）：55-61.

［13］Qing J，Meng Y D，He F，et al. Whole genome re-sequencing reveals the genetic diversity and evolutionary patterns of Eucommia ulmoides［J］. Molecular Genetics and Genomics：MGG，2022，297（2）：485-494.

［14］Song Z，Zhang Z，Dong J C，et al. Mapping immature fruit colour‐related genes via bulked segregant analysis combined with whole‐genome re‐sequencing in pepper （Capsicum annuum）［J］. Plant Breeding，2022，141（2）：277-285

［15］黃平仙，高永明，劉乃新，等. 基于全基因組重測(cè)序技術(shù)分析甜菜InDel標(biāo)記［J］. 中國(guó)糖料，2020，42（3）：1-6.

［16］Wassana K，Pumipat T，Orarat M，等. 基于SNPs全基因組測(cè)序技術(shù)對(duì)泰國(guó)辣椒地方品種遺傳多樣性分析和辣椒素含量關(guān)聯(lián)分析［J］. 辣椒雜志，2019，17（2）：37-46.

［17］宋海巖，孫淑霞，李 靖，等. 基于SSR標(biāo)記檢測(cè)與重測(cè)序技術(shù)的3個(gè)李品種鑒定與遺傳背景簡(jiǎn)析［J］. 中國(guó)南方果樹，2023，52（3）：94-101.

［18］李 淦. 基于全基因組重測(cè)序解析中國(guó)甜柿遺傳多樣性［D］. 武漢：華中農(nóng)業(yè)大學(xué)，2022.