核桃基因組學(xué)研究進(jìn)展

　　核桃具有極佳的經(jīng)濟(jì)價(jià)值與保健效果，富含人體必需的營養(yǎng)物質(zhì)和微量元素，在國內(nèi)國際市場均具有極大的消費(fèi)潛力。對(duì)核桃基因組進(jìn)行研究，獲得遺傳序列信息，有助于核桃重要農(nóng)藝性狀的定位與克隆，剖析復(fù)雜性狀的遺傳機(jī)制，對(duì)開發(fā)利用核桃優(yōu)質(zhì)種質(zhì)資源和培育、改良核桃品種具有重要意義。本文簡要統(tǒng)計(jì)了已進(jìn)行全基因組測序的核桃屬物種，并對(duì)核桃基因組學(xué)在測序技術(shù)上的更新和基因組學(xué)在核桃中的分子標(biāo)記開發(fā)、農(nóng)藝性狀定位、物種起源分化和多組學(xué)聯(lián)合等研究內(nèi)容進(jìn)行概述，對(duì)開發(fā)利用核桃優(yōu)質(zhì)種質(zhì)資源和培育、改良核桃品種具有重要意義。

 

 

　　核桃第1個(gè)參考基因組是美國加州大學(xué)測序組裝獲得，選擇物種為Chandler，使用了SOAP denovo 2和MaSuRCA兩種組裝技術(shù)，獲得667Mbp長的基因組序列，注釋到32498個(gè)基因。伴隨著核桃參考基因組的發(fā)布，對(duì)核桃的研究可以在全基因組范圍內(nèi)進(jìn)行，在該基因組數(shù)據(jù)的支持下對(duì)核桃的多酚合成途徑進(jìn)行分析，發(fā)現(xiàn)了第2個(gè)多酚氧化酶基因（JrPPO2），基因組序列的研究為核桃加速育種和促進(jìn)復(fù)雜性狀的遺傳剖析提供了重要的工具和方法。

 

　　Ning等對(duì)云南鐵核桃進(jìn)行了染色體水平的測序組裝，使用Illumina測序平臺(tái)和Nanopore測序平臺(tái)分別獲得了38.0Gb的短reads和66.31Gb的長reads，構(gòu)建Hi-C文庫獲得可組裝到染色體的數(shù)據(jù)，獲得的基因組大小為536.50Mbp，從基因組中預(yù)測了30387個(gè)蛋白編碼基因，其中99.8%已實(shí)現(xiàn)功能注釋。Huang等對(duì)山核桃進(jìn)行基因組測序，選用品種為美國的Pawnee和中國的ZAFU-1，測序深度為288 x和248 x，分別獲得187.22Gb和178.87Gb長的序列，基因組組裝長度為651.31Mbp和706.43Mbp，基因組測序后續(xù)分析可以為山核桃的氣候適應(yīng)性和抗逆性研究提供數(shù)據(jù)支持。

 

　　近年來，隨著結(jié)構(gòu)基因組學(xué)的發(fā)展，許多新的技術(shù)（如Hi-C、Chicago和光學(xué)圖譜等）被開發(fā)應(yīng)用于基因組的輔助組裝。Marrano等在核桃參考基因組Chandler v 1.0基礎(chǔ)上，使用納米孔長讀測序（Oxford Nanopore long-read sequencing）進(jìn)行深度序列覆蓋測序，結(jié)合染色體構(gòu)象捕獲技術(shù)（chromosome conformation capture）Hi-C和Illumina測序技術(shù)，組裝完成染色體級(jí)別的高度連續(xù)核桃基因組Chandler v 2.0，新基因組的組裝N50大小增加了84.4倍，與轉(zhuǎn)錄組結(jié)合挖掘到37554個(gè)基因，新的基因組極大地提高了基因預(yù)測的準(zhǔn)確性，平均基因長度高于之前基因注釋。Zhang等選擇雜合度低的西藏品種Zhongmucha-1進(jìn)行基因組測序，結(jié)合Hi-C、遺傳、物理圖譜獲得了核桃染色體水平的基因組序列，基因組大小為540Mbp，注釋了39432個(gè)基因，基于核桃的高質(zhì)量基因組圖譜，對(duì)黑核桃、野核桃和核桃楸等5個(gè)種的contig序列進(jìn)行排序和定向，得到了5個(gè)種的假染色體級(jí)別基因組序列（http: //xhhuanglab. cn/data/juglans. html）?；蚪M學(xué)數(shù)據(jù)獲取越來越快，而對(duì)基因組組裝序列的評(píng)估是復(fù)雜的，有研究發(fā)現(xiàn)約95%的RNA-seq讀數(shù)和97.25%的Illumina測序讀數(shù)可以排列到最終的裝配上，以此提出了一種基于基因含量進(jìn)化信息預(yù)期的方法，用于評(píng)估基因組組裝和注釋完整性——BUSCO。Zhang等使用BUSCO評(píng)估基因組組裝的完整性發(fā)現(xiàn)，大約94%的核心基因能夠被檢索到。

 

　　核桃的基因組高度雜合，組裝復(fù)雜，而種間雜種的基因組通常由親本中的單倍體基因組組成。Zhu等利用種間雜種的這一特性來避免雜合性，對(duì)親本為小黑核桃×核桃（J. microcarpa×J. regia）的雜交種進(jìn)行基因組測序，使用光學(xué)圖譜技術(shù)對(duì)雜交種構(gòu)建2個(gè)光學(xué)圖譜，再將光學(xué)圖譜上的排列序列重疊群填充到骨架（MS1-56_v0）中，填充后將雜交種的基因組組裝體分配到親本中，結(jié)合高密度遺傳連鎖圖譜填充親本的基因組間隙，生成了親本的基因組數(shù)據(jù)集JrSerr_v1.0和Jm31.01_v1.0。

 

　　對(duì)于核桃類的風(fēng)傳媒物種，對(duì)子代的研究使用來自母系遺傳的細(xì)胞器DNA會(huì)更加準(zhǔn)確，植物母系遺傳葉綠體基因組沒有重組和基因轉(zhuǎn)換的現(xiàn)象，廣泛用于系統(tǒng)發(fā)育研究。核桃的葉綠體基因組具有被子植物典型的四分體結(jié)構(gòu)，由大單拷貝區(qū)（LSC）和小單拷貝區(qū)（SSC）以及分隔的一對(duì)反向重復(fù)區(qū)（IRa與IRb）組成。Peng等使用高通量測序技術(shù)對(duì)核桃葉綠體基因組進(jìn)行全核苷酸序列測序發(fā)現(xiàn)，葉綠體基因組長度為160537bp，大拷貝區(qū)長度為90095bp，小拷貝區(qū)為18412bp，2對(duì)反向重復(fù)序列長26033bp，共編碼113個(gè)基因。Hu等對(duì)中國的5種核桃（核桃、鐵核桃、山核桃、麻核桃和核桃楸）的葉綠體基因組進(jìn)行測序，葉綠體基因組長度介于159714~160367bp之間，繪制了5種核桃的基因組變異圖譜，包括SNP、InDel和SV、SSR、重復(fù)序列等結(jié)構(gòu)變異的差異。根據(jù)5種核桃的葉綠體基因組數(shù)據(jù)，將這5種核桃資源劃分為2個(gè)已知區(qū)段核桃（胡桃）組和核桃楸組，為中國地區(qū)的核桃分類、系統(tǒng)發(fā)育提供了信息。Song等對(duì)美國核桃和中國核桃雜交的砧木進(jìn)行葉綠體測序，系統(tǒng)發(fā)育結(jié)果分析顯示雜交獲得的砧木與美國核桃的親緣關(guān)系較近。Yang等對(duì)胡桃科楓楊（Pterocarya stenoptera）的葉綠體基因組測序并分析其與核桃種之間的系統(tǒng)發(fā)育關(guān)系，與包括核桃在內(nèi)的14種近緣種蛋白序列分析，確定楓楊與核桃關(guān)系較近。

 

　　截至目前，已進(jìn)行全基因組測序的核桃屬物種共13個(gè)。隨著測序技術(shù)的快速發(fā)展，核桃種的測序深度逐漸加深，組裝數(shù)據(jù)也更加精準(zhǔn)，結(jié)合二代、三代測序平臺(tái)，組裝水平達(dá)到染色體級(jí)別（表1）。

 

 

　　基因組中含有大量的非編碼重復(fù)DNA，包括轉(zhuǎn)座子、反轉(zhuǎn)錄因子、串聯(lián)重復(fù)序列以及微衛(wèi)星位點(diǎn)等。微衛(wèi)星標(biāo)記（simple sequence repeat, SSR）共顯性好、多態(tài)性高，適用于高通量制圖、群體遺傳分析以及標(biāo)記輔助育種。Aziz等評(píng)估了12個(gè)核桃種中的SSR位點(diǎn)，通過鑒定核SSR（nuSSR）、葉綠體SSR（cpSSR）、線粒體SSR（mtSSR），比較它們?cè)诓煌蚪M中的頻率與分布，共驗(yàn)證了39000個(gè)SSR位點(diǎn)，顯示其他物種的末端序列在核桃種內(nèi)保持了90%~95%的一致性。Xu等利用泡核桃基因組數(shù)據(jù)開發(fā)了32個(gè)微衛(wèi)星位點(diǎn)，將這些位點(diǎn)用于3個(gè)自然群體的60個(gè)個(gè)體進(jìn)行檢測，檢測到這些位點(diǎn)具有高度的多態(tài)性，平均每個(gè)種群有4個(gè)等位基因，在核桃、山核桃、核桃楸等近緣種中有30個(gè)位點(diǎn)擴(kuò)增成功。

 

　　全基因組分子標(biāo)記的開發(fā)和基因分型是向分子育種過渡的第一步。單核苷酸多態(tài)性（SNPs）是個(gè)體間的單堿基差異，是基因組中最豐富的遺傳變異來源，具有數(shù)量大、分布廣、突變率低的特點(diǎn)。Liao等用第二代測序技術(shù)中的Roche 454測序技術(shù)對(duì)核桃進(jìn)行測序，總共檢測到49202個(gè)核苷酸變異，包括48165個(gè)SNPs和1037個(gè)InDels，這些SNPs和InDels將為核桃的遺傳研究提供有價(jià)值的標(biāo)記資源。Marrano等開發(fā)了一種高密度Axiom J. regia 700K SNP基因分型陣列，應(yīng)用這種基因分型工具，可以對(duì)核桃的關(guān)鍵性狀進(jìn)行遺傳剖析，使用這種SNP陣列從一組9.6Mbp的全基因組變異中獲得609k的高質(zhì)量SNP，數(shù)據(jù)表明獲得的SNP是先前重測序確定的，使用獲得的數(shù)據(jù)對(duì)超過1000株核桃進(jìn)行基因分型驗(yàn)證，55.7%的SNPs屬于多態(tài)性高的類型，親代與子代關(guān)系進(jìn)行鑒定也符合孟德爾遺傳定律。Arab等使用Axiom J. regia 700K SNP基因分型技術(shù)對(duì)伊朗地區(qū)的8個(gè)不同氣候區(qū)的95個(gè)核桃樣本進(jìn)行基因分型，其中53%的SNP屬于高分辨率的多態(tài)性。Ciarmiello等對(duì)18個(gè)地理來源不同的核桃品種的內(nèi)部轉(zhuǎn)錄間隔區(qū)（核糖體的ITS1和ITS2）中的一些SNPs進(jìn)行分析，對(duì)18個(gè)品種的ITS1-5.8s-ITS2的序列進(jìn)行比對(duì)發(fā)現(xiàn)序列分為2類，種間存在多態(tài)性，證明ITS1和ITS2區(qū)域進(jìn)行SNP標(biāo)記可用于核桃品種的鑒別。

 

 

　　隨著基因組測序技術(shù)的發(fā)展，與性狀基因相連的遺傳標(biāo)記逐漸被挖掘，具有豐富序列信息的標(biāo)記可以提高育種效率，使用這些標(biāo)記構(gòu)建遺傳圖譜，繪制高密度遺傳連鎖圖譜可以清晰目的基因與染色體之間的位置。Zhu等使用特異性長度擴(kuò)增片段測序（Specific Length Amplified Fragment sequencing，SLAF-seq）技術(shù)獲得153820個(gè)SLAF標(biāo)記，使用其中2577個(gè)多態(tài)性標(biāo)記，對(duì)F1群體構(gòu)建了含有16個(gè)連鎖群（linkage Group，LG）的高密度遺傳圖譜，總圖譜標(biāo)記覆蓋2457.82cM；根據(jù)區(qū)間映射（LOD＞3.0）檢測數(shù)量性狀，在LG14區(qū)間上的165.51~176.33cM檢測到1個(gè)炭疽病抗性QTL。

 

　　基于基因組對(duì)數(shù)量性狀進(jìn)行進(jìn)一步分析是數(shù)量性狀基因座定位（Quantitative Trait Loci，QTL）和全基因組關(guān)聯(lián)分析（Genome Wide Association Study，GWAS）手段對(duì)核桃基因組功能基因的注釋研究，QTL與GWAS結(jié)合使核桃從傳統(tǒng)育種走向基因組輔助育種。Marrano等使用QTL定位和GWAS關(guān)聯(lián)分析的方法探究影響核桃產(chǎn)量、果仁表皮顏色、物候期等5個(gè)性狀的因素，表型性狀數(shù)據(jù)結(jié)合SNP標(biāo)記構(gòu)建的遺傳圖譜，確認(rèn)相關(guān)性狀的基因座，分別在Chr1、Chr11染色體上發(fā)現(xiàn)了與核桃成熟期、產(chǎn)量與結(jié)果相關(guān)的基因組區(qū)域，在Chr6和Chr7上2個(gè)標(biāo)記關(guān)聯(lián)糖基轉(zhuǎn)化和氧化反應(yīng)。

 

　　Sideli等以核桃的果殼縫合強(qiáng)度為切入點(diǎn)，對(duì)核桃內(nèi)果皮縫合線進(jìn)行研究，使用700K SNP Array技術(shù)對(duì)464株核桃進(jìn)行SNP多態(tài)性分析，QTL作圖和GWAS結(jié)合分析，確定導(dǎo)致性狀變異的特定位點(diǎn)，控制核桃果殼縫合強(qiáng)度的主要QTL在LG05，該QTL的基因功能是切割和降解木葡聚糖聚合物，另在LG01和LG11上分別有1個(gè)次要QTL，也可能對(duì)核桃縫合線起調(diào)控作用。

 

　　對(duì)核桃的物候期和結(jié)果習(xí)性進(jìn)行研究有助于選育高產(chǎn)品種，Bernard等對(duì)全球各地的170份核桃材料進(jìn)行GWAS分析，并對(duì)與物候期相關(guān)的78個(gè)F1子代個(gè)體進(jìn)行基因組分型，結(jié)果發(fā)現(xiàn)在核桃第1條染色體上的SNP位點(diǎn)與核桃發(fā)育過程中的芽裂和雌花開花有很強(qiáng)的關(guān)聯(lián)性，該位點(diǎn)得到在該區(qū)域的QTL證實(shí)。果個(gè)大、出仁率高、易取仁是核桃商品出售的主要標(biāo)準(zhǔn)，Bernard等對(duì)核桃的果實(shí)特性包括性狀、大小、營養(yǎng)成分等品質(zhì)相關(guān)的25個(gè)性狀進(jìn)行全基因組關(guān)聯(lián)分析，結(jié)果發(fā)現(xiàn)60多個(gè)標(biāo)記與性狀關(guān)聯(lián)，包含2個(gè)與核桃果仁體積與重量關(guān)聯(lián)的SNP位點(diǎn)。

 

　　核桃仁澀皮顏色越淺越受到消費(fèi)者青睞，核桃仁的澀皮顏色在不同品種、不同發(fā)育階段表現(xiàn)不同。Sideli等對(duì)Chandler×Ldaho雙親后代的168株后代和本地的528株核桃進(jìn)行GWAS分析，發(fā)現(xiàn)與核桃仁澀皮顏色的表型在Chandler遺傳圖譜上有重疊區(qū)域，在Chr01等多個(gè)染色體上小效應(yīng)QTL位點(diǎn)，QTL與標(biāo)記錨定染色體定位于Chr01、Chr10、Chr14和Chr16上，包括轉(zhuǎn)錄因子MYB113在內(nèi)的12個(gè)基因被認(rèn)定為可能與核桃仁澀皮顏色色素沉著有關(guān)。

 

　　在核桃生產(chǎn)栽培中，與開花、落葉相關(guān)的性狀是植物生命周期的關(guān)鍵因素，也是決定核桃生產(chǎn)力的重要因素。Bükücü等對(duì)土耳其的188份核桃材料，通過全基因組關(guān)聯(lián)分析進(jìn)行基因型變異分析，結(jié)合3年與葉芽萌發(fā)、開花等13個(gè)開花、落葉相關(guān)的表型性狀數(shù)據(jù)，發(fā)現(xiàn)了與花序數(shù)量和果實(shí)結(jié)實(shí)等性狀相關(guān)的標(biāo)記，關(guān)聯(lián)分析結(jié)果顯示16個(gè)QTL位點(diǎn)至少與2個(gè)表型性狀相關(guān)，其中QTL05關(guān)聯(lián)數(shù)量最多，與7個(gè)性狀相關(guān)（開花期、萌芽期、雄花始花期、雌花始花期、雄花終花期、雌花終花期以及堅(jiān)果果形），這些鑒定出來的與開花、落葉相關(guān)的標(biāo)記位點(diǎn)與顯著性QTL為以后童期選擇提供了數(shù)據(jù)。

 

 

　　核桃的基因組測序長度在640M~990Mbp之間，測序結(jié)果的連續(xù)性和準(zhǔn)確性支持對(duì)基因組變異的注釋，基因組注釋得清晰對(duì)于了解核桃植物及其進(jìn)化分類都具有重要的意義。Stevens等對(duì)北加州黑核桃、黑核桃、小黑核桃和核桃進(jìn)行基因組深度重測序，重測序結(jié)果將核桃屬分為3類：①黑核桃組：J. nigra、J. microcarpa、J. hindsii；②核桃組：J. regia、J. sigillata；③核桃楸組：J. cathayensis。利用重測序數(shù)據(jù)對(duì)多酚氧化酶位點(diǎn)PPO的結(jié)構(gòu)和序列變化進(jìn)行深度分析，研究結(jié)果顯示PPO1位點(diǎn)在譜系上表現(xiàn)出加速分化和氨基酸置換過剩的現(xiàn)象，這些現(xiàn)象可能導(dǎo)致核桃和山核桃的馴化。Mu等使用全基因組限制位點(diǎn)相關(guān)的DNA測序方法（2b-RAD）對(duì)野生和栽培麻核桃種進(jìn)行研究，通過構(gòu)建物種進(jìn)化樹和測序結(jié)果得出麻核桃是核桃和核桃楸的雜交后代。

 

　　Zhang等對(duì)19種核桃種的80個(gè)個(gè)體，以及近緣種楓楊和化香樹（Platycarya strobilacea Sieb）進(jìn)行全基因組重測序，使用單拷貝核基因的系統(tǒng)發(fā)育網(wǎng)絡(luò)分析、全基因組位點(diǎn)模式概率等方法得出核桃及其地方品種泡核桃是美洲核桃和亞洲核桃的雜交后代，而灰核桃是亞洲核桃大量滲入美洲黑核桃基因組的結(jié)果，通過貝葉斯計(jì)算模型發(fā)現(xiàn)其在上新世晚期（3.45Mya）開始分化，而2種親本血統(tǒng)在歐洲已經(jīng)滅絕。Mu等基于RAD-seq技術(shù)和葉綠體基因組數(shù)據(jù)對(duì)核桃屬的系統(tǒng)發(fā)育進(jìn)行重新分析，分析結(jié)果在群體水平上支持胡桃亞科的拓?fù)浣Y(jié)構(gòu)，最終在化石記錄、系統(tǒng)發(fā)育和譜系分化時(shí)間等的證實(shí)下提出核桃木亞科起源于北美，通過白令海峽和北大西洋陸橋向歐亞大陸遷移。

 

　　通過對(duì)核桃基因組的研究預(yù)測過去的氣候生態(tài)信息，并更新核桃避難所的位置。Aradhya等對(duì)高加索地區(qū)、中亞、東亞、亞洲西南和東歐等317份核桃材料643種基因型進(jìn)行遺傳多樣性和群體結(jié)構(gòu)分析，分析核桃現(xiàn)今期、末次盛冰期（Last Glacial Maximum，LGM）、末次間冰期（Last Interglacia，LIG）的空間與時(shí)間分布，結(jié)果表明LGM期間核桃的分布范圍主要局限于南高加索地區(qū)的南緯地區(qū)、西藏西南部的中亞和南亞地區(qū)、印度東北部、錫金和不丹的喜馬拉雅地區(qū)以及中國東南部，在北緯30～45°的地理區(qū)域內(nèi)的避難所存活。

 

　　隨著末次盛冰期到來，迫使適應(yīng)溫暖的環(huán)北方植物群遷移，并被限制在亞洲東部、北美和歐亞大陸西南部的大型南方避難所。Bai等對(duì)亞洲南北2個(gè)地區(qū)的核桃進(jìn)行了多樣性分析，2個(gè)地區(qū)間存在一條干旱隔離帶，結(jié)果北部為核桃楸屬和日本核桃（Juglans ailantifolia Carr），南部為山核桃屬，推測主要原因?yàn)楸辈吭忌矫}的地域差異而分化，此后基本保持分離。Wang等利用12個(gè)SSR位點(diǎn)對(duì)中國東北和朝鮮半島的19個(gè)核桃群進(jìn)行了末次冰期后物種的擴(kuò)散模式探究，結(jié)果發(fā)現(xiàn)在末次冰期后物種擴(kuò)散的過程中，東北地區(qū)的核桃?guī)缀鯖]有遺傳多樣性缺失的現(xiàn)象，推測是由于在北方種群傳播過程中缺乏長距離的傳播機(jī)制，因而遺傳多樣性的維度沒有下降。Ebrahimi等對(duì)北美核桃、亞洲核桃和核桃的基因型和抗寒性進(jìn)行了評(píng)估，結(jié)合耐冷基因的表達(dá)量和系統(tǒng)發(fā)育分析得出，核桃與其他亞洲核桃核基因組聚為一組，且核桃的遺傳多樣性低于其他樣本，鑒定的耐寒冷基因最低，分子機(jī)理與生理數(shù)據(jù)高度相關(guān)，可以有效地表征核桃品種的耐寒性。

 

　　現(xiàn)在核桃群體種質(zhì)資源分布的多樣性除代表冰川后的擴(kuò)張、定居和種植外，還包括自然和人類選擇和馴化的復(fù)雜的相互作用。Bai等推算核桃物種演化時(shí)間，對(duì)11個(gè)溫帶核桃物種的基因組進(jìn)行兩兩測序的馬爾科夫合并方法，估算核桃有效種群大小（Ne）的變化軌跡，再使用G-PHOCS方法估算多個(gè)基因組間的分化時(shí)間，結(jié)果得出核桃屬在100萬年前開始分化，這一結(jié)果也被G-PHOCS對(duì)發(fā)散時(shí)間的估計(jì)所證實(shí)，但它們對(duì)冰期后的氣候反應(yīng)并不相似，Ne變化的時(shí)間和幅度不同，推測核桃種群的進(jìn)化歷程中不單單存在環(huán)境改變，與特定病原體的共同進(jìn)化交互作用等物種特異性因素可能發(fā)揮了關(guān)鍵作用。

 

　　Paola等使用14個(gè)微衛(wèi)星位點(diǎn)對(duì)91個(gè)歐洲核桃進(jìn)行遺傳多樣性評(píng)價(jià)，構(gòu)建歐洲核桃的種群歷史譜系和傳播途徑，結(jié)果表明在羅馬帝國時(shí)期就存在東歐、西歐的種間雜交，而在過去的5個(gè)世紀(jì)里，在歐洲東北部和西歐出現(xiàn)了人口規(guī)模的擴(kuò)張和隨后的下降，核桃在歐洲的實(shí)際分布是末次盛冰期后多個(gè)避難所的擴(kuò)張或收縮和過去5000年來人類對(duì)其開發(fā)的綜合作用的結(jié)果。

 

 

　　獲得高質(zhì)量的基因組，將使多倍體高重復(fù)的復(fù)雜物種的轉(zhuǎn)錄分析更加準(zhǔn)確。Hu等利用Illumina測序技術(shù)和De novo組裝獲得平均長度747bp的轉(zhuǎn)錄組數(shù)據(jù)，利用這些數(shù)據(jù)鑒定到63個(gè)新的SSR標(biāo)記，使用這些標(biāo)記對(duì)4個(gè)群體中的4~14個(gè)個(gè)體進(jìn)行標(biāo)記檢測，檢測到20個(gè)SSR標(biāo)記具有多態(tài)性且易于擴(kuò)增。Zhao等對(duì)核桃脂肪合成的3個(gè)發(fā)育時(shí)期（初期發(fā)育期、快速發(fā)育期、成熟發(fā)育期）進(jìn)行轉(zhuǎn)錄組測序分析，3個(gè)發(fā)育階段共獲得68.18Gb的數(shù)據(jù)，92%~94%的clean data可以比對(duì)到參考基因組，分析脂肪合成相關(guān)的關(guān)鍵基因的表達(dá)水平，確定ACCase、LACS和FAD7是與脂肪合成相關(guān)的關(guān)鍵基因。Quan等對(duì)核桃花芽分化不同時(shí)期進(jìn)行cDNA文庫測序及石蠟包埋組織切片方法觀察形態(tài)變化，雌花芽與葉芽間差異表達(dá)的DEGs有374個(gè)，構(gòu)建與花期有關(guān)的DEGs共表達(dá)網(wǎng)絡(luò)，鑒定出CRY2和NF-YA基因是調(diào)控花期的核心DEGs，通過cDNA測序?qū)颂页苫ㄟ^程的表達(dá)基因進(jìn)行篩選，為核桃的開花機(jī)理提供了理論依據(jù)。

 

　　將基因組數(shù)據(jù)和轉(zhuǎn)錄組數(shù)據(jù)補(bǔ)充現(xiàn)有的蛋白質(zhì)序列數(shù)據(jù)庫，利用這種蛋白基因組數(shù)據(jù)庫可以極大地改善MS/MS搜索結(jié)果。Cary等通過基因組測序和轉(zhuǎn)錄組、蛋白組的數(shù)據(jù)庫檢索光譜，并將提取的核桃進(jìn)行液相色譜-質(zhì)譜（LC-MS/MS）分析，最終在基因組數(shù)據(jù)庫中確定出10個(gè)獨(dú)特多肽段，轉(zhuǎn)錄組數(shù)據(jù)庫中的識(shí)別多肽數(shù)量增加了20%，為核桃等堅(jiān)果物種的蛋白數(shù)據(jù)庫奠定了基礎(chǔ)。Zaini等對(duì)早實(shí)、晚實(shí)核桃的種皮顏色進(jìn)行研究，揭示了在采收期中常見蛋白與特異蛋白的321個(gè)差異明顯的種皮膜層蛋白，蛋白質(zhì)組學(xué)分析檢測到4937個(gè)蛋白，涉及到抗氧化、次級(jí)代謝等途徑。

 

 

　　隨著新測序技術(shù)的不斷更新，基因組學(xué)相關(guān)研究實(shí)現(xiàn)快速發(fā)展，以高通量測序技術(shù)為基礎(chǔ)的各種組學(xué)研究廣泛應(yīng)用于核桃等植物的各個(gè)領(lǐng)域。測序技術(shù)的成熟加上測序成本的降低，越來越多的果樹物種的基因信息被破解。而果樹的全基因組測序工作仍需繼續(xù)推進(jìn)，核桃等果樹的功能基因組學(xué)研究遠(yuǎn)不及水稻、玉米等農(nóng)作物，主要原因仍為果樹是多年生植物，有性繁殖周期長，建立雜交群體至開花結(jié)果得到想要研究的性狀需要很長時(shí)間；另外對(duì)功能基因的定位與驗(yàn)證也受到完善的表達(dá)體系的限制，除蘋果等國內(nèi)主栽經(jīng)濟(jì)果樹外，大部分果樹仍需借助模式植物進(jìn)行驗(yàn)證。

 

　　隨著測序技術(shù)的不斷進(jìn)步，通過測序獲得了越來越多果樹的全基因組序列，同樣在轉(zhuǎn)錄組、蛋白組、代謝組等生物信息數(shù)據(jù)也在不斷積累，而建立一種包含基因組、轉(zhuǎn)錄組、蛋白組等相關(guān)數(shù)據(jù)的綜合性數(shù)據(jù)庫顯得尤為重要，目前在柑橘屬（Citrus Genome Database，https: //www. Citrusgenom-edb. org/）、梨屬（Pear Genome Project，http: //pear-genome. njau. edu. cn/）已擁有基因組數(shù)據(jù)庫，通過對(duì)綜合數(shù)據(jù)庫的生信數(shù)據(jù)進(jìn)行分析、儲(chǔ)存以及整合，可為果樹基礎(chǔ)研究和應(yīng)用研究提供可靠數(shù)據(jù)信息。對(duì)核桃的全基因組測序能夠推進(jìn)核桃的基因組學(xué)研究，豐富核酸序列和基因資源，為進(jìn)一步開發(fā)分子標(biāo)記、挖掘重要功能基因和解析生長發(fā)育機(jī)制等分子生物學(xué)研究提供重要基礎(chǔ)和依據(jù)，加快核桃分子輔助育種進(jìn)程。

 

　　聲  明：本文摘編自《中國果樹》2022年第2期“核桃基因組學(xué)研究進(jìn)展”（武鵬雨，劉婷婷，包建平，虎海防，馬凱，張銳）。