武鵬雨，劉婷婷，包建平，虎海防，馬 凱，張 銳

（1 塔里木大學(xué)園藝與林學(xué)學(xué)院，新疆阿拉爾843300）（2 新疆佳木果樹學(xué)國家長期科研基地）（3 新疆農(nóng)業(yè)科學(xué)院）

核桃屬（Juglans）隸屬于胡桃科（Juglandaceae），是世界四大堅果之一，是不可或缺的堅果和木本油料樹種。核桃屬有21 個二倍體種，均有2n＝32 條染色體[1]，根據(jù)系統(tǒng)發(fā)育關(guān)系和果實形態(tài)可分為四大類[2-4]，分別為核桃楸組（Cardiocaryon）的山核桃（Juglans cathayensisDode）、麻核桃（Juglans hopeiensisHu）、核桃楸（Juglans mandshuricaMaxim），核桃組（Juglans）的泡核桃（Juglans sigillataDode）和核桃（Juglans regiaL.），灰核桃組（Trachycaryon）的灰核桃（Juglans cinereaL.），黑核桃組（Rhysocaryon）的北加州黑核桃（Juglans hindsiiRehd）、黑核桃（Juglans nigraL.）以及小黑核桃（Juglans microcarpa）。核桃營養(yǎng)豐富，富含人體必需的ω-3 脂肪酸等營養(yǎng)物質(zhì)，并對延緩衰老，治療心臟病、糖尿病等有一定作用[5]。

我國核桃品種豐富，在品質(zhì)改良和品種選育過程中取得了顯著成績，對新品種的選育方法以雜交和優(yōu)選為主，也伴隨著培養(yǎng)周期長、種間雜交困難等問題，分子育種能有效縮短育種時間，快速篩選含有目標(biāo)性狀的植株，通過建立核桃屬植物的種質(zhì)資源庫及雜交后代群體等，開發(fā)應(yīng)用遺傳標(biāo)記，將分子標(biāo)記輔助育種技術(shù)、轉(zhuǎn)基因技術(shù)與傳統(tǒng)育種技術(shù)相結(jié)合，讓育種過程具有更好的目標(biāo)性[6]。

基因組學(xué)（Genomics）是對生物體內(nèi)核苷酸序列結(jié)構(gòu)、基因功能、分子遺傳進(jìn)化等機理的研究[7]。在模式植物擬南芥[8]中進(jìn)行了全基因組測序分析，獲得了120 Mbp 的基因組序列，編碼20 000 個基因，標(biāo)志著對植物學(xué)的研究進(jìn)入基因組學(xué)階段。隨著高通量測序技術(shù)的不斷完善，葡萄[9]、番木瓜[10]等國內(nèi)外果樹先后進(jìn)行基因組測序組裝，隨后越來越多的果樹完成基因組測序工作，果樹進(jìn)入基因組學(xué)篇章。核桃基因組學(xué)的研究通過對核桃的基因組進(jìn)行測序，定位核桃染色體中的序列信息，從而構(gòu)建核桃的種質(zhì)資源數(shù)據(jù)庫，獲取含有豐富遺傳信息的分子標(biāo)記，錨定相關(guān)功能基因和研究核桃物種的起源與演化。本文簡要統(tǒng)計了已進(jìn)行全基因組測序的核桃屬物種，并對核桃基因組學(xué)在測序技術(shù)上的更新和基因組學(xué)在核桃中的分子標(biāo)記開發(fā)、農(nóng)藝性狀定位、物種起源分化和多組學(xué)聯(lián)合等研究內(nèi)容進(jìn)行概述，對開發(fā)利用核桃優(yōu)質(zhì)種質(zhì)資源和培育、改良核桃品種具有重要意義。

1 核桃屬物種基因組測序研究

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

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

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

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

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

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

表1 核桃已測序基因組

2 基因組學(xué)在核桃上的應(yīng)用研究

2.1 分子標(biāo)記的開發(fā)

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

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

2.2 重要農(nóng)藝性狀的定位

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

基于基因組對數(shù)量性狀進(jìn)行進(jìn)一步分析是數(shù)量性狀基因座定位（Quantitative Trait Loci，QTL）和全基因組關(guān)聯(lián)分析（Genome Wide Association Study，GWAS）手段對核桃基因組功能基因的注釋研究，QTL 與GWAS 結(jié)合使核桃從傳統(tǒng)育種走向基因組輔助育種。Marrano 等[32]使用QTL 定位和GWAS關(guān)聯(lián)分析的方法探究影響核桃產(chǎn)量、果仁表皮顏色、物候期等5 個性狀的因素，表型性狀數(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 個標(biāo)記關(guān)聯(lián)糖基轉(zhuǎn)化和氧化反應(yīng)。

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

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

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

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

2.3 利用全基因組測序研究物種起源與演化

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

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

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

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

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

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

2.4 基因組與其他組學(xué)結(jié)合，挖掘基因表達(dá)信息

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

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

3 展 望

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

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