利用GBS技術研究240份寬皮柑橘的系統(tǒng)演化

王小柯，江東,2，孫珍珠

（1西南大學柑橘研究所，重慶 400712；2中國農(nóng)業(yè)科學院柑橘研究所，重慶 400712）

利用GBS技術研究240份寬皮柑橘的系統(tǒng)演化

王小柯1，江東1,2，孫珍珠1

（1西南大學柑橘研究所，重慶 400712；2中國農(nóng)業(yè)科學院柑橘研究所，重慶 400712）

【目的】GBS（genotyping-by-sequencing）是一種高效而經(jīng)濟的SNP（單核苷酸多態(tài)性）發(fā)掘和基因分型技術。采用GBS技術對240份寬皮柑橘進行基因分型，以闡明一些野生寬皮柑橘和地方品種的遺傳背景，為其起源和演化研究提供更可靠的證據(jù)?！痉椒ā窟x用國家柑橘種質重慶資源圃保存的具有廣泛遺傳多樣性和地理起源的240份寬皮柑橘作為材料，利用EcoR I限制性內切酶消化基因組DNA后構建GBS文庫；然后進行Illumina HiSeq PE150二代測序獲得短讀序列，通過BWA軟件將序列映射到克里曼丁參考基因組上，再利用SAMTOOLS軟件鑒定SNP位點。依據(jù)SNP的基因分型結果，采用鄰近法構建系統(tǒng)演化樹，并進行主成分分析?！窘Y果】利用GBS簡化基因組測序技術對240份寬皮柑橘進行測序，共獲得96.3 Gb的測序數(shù)據(jù)，平均每個樣本測序數(shù)據(jù)為401.26 Mb，經(jīng)過測序深度為4X、Miss0.2、次要等位基因頻率（MAF）＞0.01的篩選條件過濾，最后共獲得了114 200個高質量的SNP位點。主成分分析結果顯示240份寬皮柑橘被分為4大類，其中溫州蜜柑亞群、野生寬皮柑橘亞群可明顯區(qū)分于其他寬皮柑橘。利用系統(tǒng)演化樹可將240份寬皮柑橘劃分到11個類群中。系統(tǒng)演化樹和主成分分析都揭示了不同地理來源和特定形態(tài)的寬皮柑橘在遺傳水平上存在明顯的差異，比如來源于日本的溫州蜜柑、歐美的克里曼丁橘及其雜種后代，以及中國南方的野生寬皮柑橘由于地理分布不同而形成了較為獨特的類型，彼此間能夠相互區(qū)分開。進化樹結果表明中國南、北不同地域的寬皮柑橘可能存在不同的演化路徑，南嶺山脈及南方地區(qū)的野生寬皮柑橘、酸橘和目前南方地區(qū)栽培的砂糖橘存在較近的起源演化聯(lián)系，而北方寬皮柑橘的演化卻與寬皮柑橘中的古老地方品種存在緊密聯(lián)系。人工雜交育種、長期的人工選擇和馴化形成了不同類型的寬皮柑橘，同時也導致寬皮柑橘遺傳多樣性的增加。另外，一些寬皮柑橘資源中的可疑親本也通過GBS技術得以準確鑒定。本研究表明沃柑與金諾橘有較近的親緣關系?！窘Y論】GBS技術用于柑橘種質資源的基因分型高效可靠，建立的系統(tǒng)演化樹可以對240份寬皮柑橘進行準確劃分，與用植物形態(tài)學劃分的結果高度吻合。另外，GBS技術用于資源材料的準確鑒定和親緣關系的研究，可為柑橘植物新品種權的保護提供可靠的技術支撐。

GBS；寬皮柑橘；系統(tǒng)演化；品種鑒定

0 引言

【研究意義】中國是世界寬皮柑橘的起源中心，寬皮柑橘種類豐富多樣，至今仍有野生寬皮柑橘的分布[1]，這為研究寬皮柑橘的起源和演化提供了重要材料。現(xiàn)有的研究表明寬皮柑橘是柑橘屬下的基本種之一，但由于寬皮柑橘與柑橘屬下的其他種易發(fā)生種間雜交，導致其種類繁多，不僅增加了鑒定難度，其亞種或類群的劃分一直以來也存在較大的爭議[2]。通過基因分型對寬皮柑橘種質進行甄別和鑒定不僅是資源收集、保存、評價和利用的基本要求，也為寬皮柑橘品種類群的劃分、起源演化的研究、柑橘植物新品種權的保護和認定提供更加科學、充分的依據(jù)。在長期的栽培、馴化歷史過程中，來源于不同地域的寬皮柑橘可能保留了明顯的具有地域特征的遺傳差異。因此，對不同地理來源的寬皮柑橘進行遺傳多樣性研究，有助于了解這些材料的遺傳背景，從而為寬皮柑橘的育種提供可利用的基因資源。【前人研究進展】目前對寬皮柑橘種質的鑒定除利用植物形態(tài)學特征進行分辨外，也采用分子標記技術[3-4]和基因測序技術[5]。儲春榮等[3]利用SSR、Indel等分子標記證實了蘇州地區(qū)的黃皮橘、朱紅橘是兩類重要且較為古老的橘類資源，長江上游流域的紅橘是與朱紅橘不同的另一古老橘類品種；高恒錦等[4]利用類似的方法證實了起源于中國南嶺山脈的野生寬皮柑橘是一類獨特的原始野生寬皮柑橘類型，表明原始野生柑橘可能是寬皮柑橘中的初生種。WU等[5]通過對柑橘全基因組進行比較分析后，發(fā)現(xiàn)克里曼丁的基因組中有柚類基因的滲入，而來自中國的 1份野生寬皮柑橘與其他寬皮柑橘存在明顯的遺傳差異。近來隨著高通量測序技術的發(fā)展，基因分型的成本持續(xù)降低，GBS技術作為第二代深度測序基礎上發(fā)展起來的簡化基因組測序技術，通過采用酶切加標簽的方法，使多樣本高通量平行測序得以實現(xiàn)[6-7]。這不僅大大降低了基因測序的成本，也使對大樣本全基因組的基因分型成為可能，對深入了解種質資源的遺傳背景和系統(tǒng)演化具有重要意義[8-10]。同時GBS獲得的短讀序列可通過有參或無參基因組的形式進行拼接組裝，進而獲得高密度的SNP標記，利用這些SNP標記或者開發(fā)的bin標記可進行遺傳圖譜構建[11-12]、遺傳圖譜加密[13-15]、全基因組關聯(lián)分析（GWAS）[16-18]及基因組的輔助組裝等研究。目前GBS技術作為基因分型的重要手段，已經(jīng)在遺傳圖譜構建[11-12,19-20]、遺傳選擇[10]、基因多樣性研究[21-23]、種質鑒定[24-25]及品種識別[26-27]等領域得到廣泛的應用?！颈狙芯壳腥朦c】寬皮柑橘中的古老地方品種、野生資源在中國極其豐富，但國外對這些資源的研究很少涉及，因此，中國的野生和地方寬皮柑橘品種的遺傳多樣性值得深入研究。基于GBS技術的低成本、高通量和不依賴參考基因組等特點[28]，將其用于寬皮柑橘種質資源的基因分型是可行的?！緮M解決的關鍵問題】利用 GBS技術對240份寬皮柑橘的野生資源、地方品種以及國內外培育的雜種材料進行基因分型，通過系統(tǒng)演化和遺傳多樣性分析，以闡明這些寬皮柑橘的遺傳背景及起源和演化。

1 材料與方法

試驗于2016年在中國農(nóng)業(yè)科學院柑橘研究所進行。

1.1 試驗材料

本研究采用的240份寬皮柑橘材料均采自國家果樹種質重慶柑橘資源圃，這些寬皮柑橘種質來自于不同時期、不同地區(qū)的資源收集活動，包括起源于中國南方地區(qū)的野生寬皮柑橘15份，地方品種77份，選育品種140份，遺傳材料8份。在選育品種中包括有清見雜種材料11份，愛媛30號及愛媛28號雜種材料19份。這240份寬皮柑橘材料中來源于日本的溫州蜜柑有40份，來源于歐美的寬皮柑橘材料32份，尼泊爾材料5份，其余為國內材料，具體材料名單見電子版附表1。

1.2 試驗方法

1.2.1 DNA提取和GBS文庫構建 2016年春季采集所有寬皮柑橘資源的嫩葉，利用植物基因組DNA提取試劑盒（Magen HiPure Plant DNA Mini Kit，Guangzhou，China）提取DNA。提取的樣品DNA送諾禾致源進行DNA質檢、建庫和測序，利用Qubit? 2.0熒光測定計（Invitrogen，Carlsbad，USA）檢測核酸濃度，同時在1%的瓊脂糖凝膠上100 V電泳40 min檢測DNA純度，高質量的DNA用于GBS文庫構建和測序。每個樣本取1.5 μg DNA用于文庫構建，首先應用限制性內切酶EcoR I對基因組進行酶切，酶切后的片段兩端利用T4 DNA連接酶（NEB）加上適配體接頭，對每個樣品進行擴增，然后對樣品進行混合，電泳回收350—400 bp區(qū)間的DNA條帶，割取的片段進行純化，純化后產(chǎn)物用于測序，測序反應在Illumina HiSeq測序平臺上進行雙末端150 bp的測序。GBS分析的重復性、可靠性檢測使用電子版附表1中56號（1-1）、57號（專有橘）和58號（1-2）3份種質材料作為生物學重復。

1.2.2 單核苷酸多態(tài)性（SNP）鑒定 每條測序數(shù)據(jù)按照條碼劃分到對應的樣本中，對每條序列進行嚴格的質控篩選，去除原始測序序列中的接頭序列和低質量的短讀序列，包括單端測序中含有的N數(shù)量超過該條序列長度的10%，或者單端測序中低質量（Q≤5）堿基數(shù)超過該條序列長度50%的序列。剩下的高質量序列采用BWA（Burrows-Wheeler Aligner）程序將其匹配到克里曼丁參考基因組上（ftp://ftp.ncbi.nlm.nih. gov/genomes/all/GCF_000493195.1_Citrus_clementina_ v1.0）。SNP檢測采用Samtools軟件進行，對齊的匹配文件用Samtools軟件轉換為BAM文件后進行變異調用，Samtools收集BAM文件中的匯總信息，計算可能基因型的似然值，再利用Bcftools應用先驗值進行變異調用，采用貝葉斯模型檢測群體中的多態(tài)性位點獲得VCFs文件。

1.2.3 群體進化樹分析 采用SNP基因分型數(shù)據(jù)估測240份寬皮柑橘的系統(tǒng)演化，根據(jù)群體遺傳學特征的共同點或差異推斷出它們的親緣關系遠近。利用個體SNPs的檢測值計算群體間的遺傳距離。兩個個體i和j之間的p-距離通過如下公式計算：

式中，L表示高質量SNPs區(qū)域長度，dij表示兩個個體為不同基因型值的賦值表。運用 TreeBest（http:// treesoft.sourceforge.net/treebest.shtml）軟件計算距離矩陣，以此為基礎，通過鄰接法（neighbor-joining method）構建系統(tǒng)進化樹。引導值（bootstrap values）經(jīng)過達1 000次計算獲得。

1.2.4 主成分分析 通過GCTA（http://cnsgenomics. com/software/gcta/pca.html）軟件計算特征向量以及特征值，并利用R軟件繪制PCA分布圖。

2 結果

2.1 測序質量

GBS測序得到的原始圖像數(shù)據(jù)文件經(jīng)堿基識別分析轉化為原始測序序列，240個寬皮柑橘樣本的總測序數(shù)據(jù)量為96.31 Gb，去除低質量的序列后，產(chǎn)生的高質量序列數(shù)據(jù)量為96.3 Gb，平均每個樣本數(shù)據(jù)量為401.26 Mb。測序質量較高（Q20≥93.78%、Q30≥85.19%），GC分布正常，群體樣本與克里曼丁參考基因組的平均比對率為93.79%，表明本試驗的樣本與參考基因組克里曼丁有較高的相似性。序列平均覆蓋度為17.82%（至少有一個堿基的覆蓋），至少有4個堿基的平均覆蓋度為8.60%。利用Samtools軟件檢測后共獲得765 252個SNP位點，經(jīng)過條件為測序深度4×、Miss0.2、次要等位基因頻率（MAF）＞0.01的過濾后，最后共獲得114 200個高質量的SNP位點用于群體進化樹分析。

2.2 主成分分析

本研究選擇的 240份寬皮柑橘涵蓋了寬皮柑橘中的大部分類型，利用PCA分析方法對240份寬皮柑橘材料的親緣關系進行了分析，根據(jù)第一主成分和第二主成分可將240份寬皮柑橘分為4個亞群，分別為類群Ⅰ（溫州蜜柑亞群，包括溫州蜜柑、早熟蜜橘和諾瓦橘），類群Ⅱ（野生寬皮柑橘亞群，包括細皮狗屎柑、大坑野橘、聶都野橘、莽山橘、粗皮狗屎柑、道縣野橘），類群Ⅲ（包括大浦、早香和橘橙22-53），類群Ⅳ（其余寬皮柑橘及其雜種品種，包括酸橘、年橘、椪柑、紅橘、朱紅橘、地中海橘、早橘、橘橙或橘柚雜種、清見雜種、克里曼丁橘等）。具體結果見圖1。

圖1 240份寬皮柑橘的主成分分析圖Fig. 1 The principal component analysis of 240 mandarins

2.3 群體進化樹分析

群體進化樹的結果表明，240份寬皮柑橘種質聚為一個大類，支持寬皮柑橘是柑橘屬下一個基本種的說法。來源于不同地理區(qū)域的寬皮柑橘在遺傳上存在明顯的差異，比如來源于日本的溫州蜜柑、來源于歐美的克里曼丁橘及其雜種與中國的絕大部分寬皮柑橘品種間存在明顯的遺傳差異，表明地理來源的不同可能是引起寬皮柑橘遺傳差異的重要原因之一。根據(jù)進化樹對240份寬皮柑橘進一步細分后，可將這些材料劃分到11個寬皮柑橘亞群中（具體分類亞群見電子版附表1），這11個亞群分別是野生寬皮柑橘及其衍化種（類群Ⅰ，包括小果類型的酸橘、粗皮狗屎柑、大坑野橘、道縣野橘、莽山橘、聶都野橘、年橘、砂糖橘、八月橘、馬水橘、青皮蜜橘等），古老栽培品種（類群Ⅱ，包括土橘、建柑、椪柑、紅橘、地中海橘、槾橘等），朱紅橘（類群Ⅲ，包括朱紅橘、南橘、料紅、南豐蜜橘等），早橘（類群Ⅳ，包括本地早、城固冰糖橘、旺蒼皺皮柑、克里曼丁×本地早雜種等），歐美橘橙雜種（類群Ⅴ，包括王柑、默科特、恩科爾、麗紅、金諾、無核沃柑等），橘柚雜種（類群Ⅵ，包括明尼奧拉、清峰、漢源黃果柑、湘慈43號等），清見×椪柑雜種群（類群Ⅶ），愛媛系列雜種（類群Ⅷ，包括愛媛28號雜種、愛媛30號雜種），克里曼丁橘類（類群Ⅸ，主要為克里曼丁橘及其雜種后代，如費爾柴爾德、弗來蒙特橘、秋輝、帕森橘等以及其他類型雜種），溫州蜜柑（類群Ⅹ），清見雜種（類群Ⅺ，清見，清見與天草的雜交種，清見與明尼奧拉的雜交種）。進化樹的分類結果與利用植物形態(tài)學分類的結果是高度吻合的。

2.3.1 系統(tǒng)演化 從群體進化樹的結果來看，中國南、北不同地域的寬皮柑橘可能存在不同的演化路徑。來源于南嶺山脈及南方地區(qū)的野生寬皮柑橘、寬皮柑橘中小果類型的酸橘、華南地區(qū)廣泛種植的砂糖橘、八月橘、馬水橘、青皮蜜橘等存在較近的親緣關系，這些材料在地理分布上也主要位于華南地區(qū)，而分布于長江下游流域的小果寬皮柑橘中的乳橘類型，如南豐蜜橘、東華蜜橘卻與朱紅橘有較近的親緣關系，而小果類中的本地早卻與古老寬皮柑橘中的旺蒼皺皮柑有較近的親緣關系。結合這些品種的地理分布，推測南方地區(qū)栽植的砂糖橘等類型可能是從小果類的野生寬皮柑橘演化為酸橘后，再由其演化為現(xiàn)代栽培寬皮柑橘品種。同時還發(fā)現(xiàn)南方地區(qū)的紅橘和北方橘區(qū)的朱紅橘之間也存在明顯的差異，這表明即便都是紅皮橘類型，南北地理氣候差異也造成了兩個不同的類群。另外從進化樹上看，中國古老栽培品種旺蒼皺皮柑與城固冰糖橘遺傳距離較近，兩者有明顯的親緣關系，并且可能是本地早的早期親本來源之一，說明中國北方寬皮柑橘的演化有原始古老柑類品種的參與。來源于歐洲地中海地區(qū)的柳葉橘（阿凡娜橘）也形成了古老栽培品種中較為明顯的一個分支，與來源于東亞地區(qū)的土橘、椪柑、紅橘、朱紅橘等彼此分開，進一步證實了不同地理來源的品種間存在較明顯的遺傳差距。

圖2 240份寬皮柑橘群體進化樹Fig. 2 The phylogenetic tree of 240 mandarins

2.3.2 遺傳多樣性 從群體進化樹的結果可發(fā)現(xiàn)，寬皮柑橘與柑橘屬其他種的種間雜交，會造成遺傳多樣性的增加，進而形成獨特的亞群。比如來源于歐美的王柑、默科特等品種形成一個獨特的亞群，這些品種多數(shù)來源于橘橙品種間的天然雜種，由于具有橙類的遺傳基因，因此這個亞群的品種多數(shù)具有果皮硬、成熟晚等特點。而一些遺傳背景更為復雜的橘橙或橘柚雜種則形成了新的亞群，如清峰、湘慈43、漢源黃果柑等。通過人工雜交獲得的雜種后代也會導致遺傳多樣性的增加，而人工選擇的偏好性使具有相似特征的品種形成新的品種亞群，比如清見與不同椪柑品種的雜交后代形成了一個新的亞群（圖2中的類群Ⅶ），這些材料在遺傳距離上與椪柑相距較遠；而愛媛 30號、愛媛28號人工雜交后代也形成了新的分支，但這些雜種后代與母本愛媛30號、愛媛28號的遺傳距離相對較近。來源于歐美的克里曼丁橘及其雜種后代在進化樹上也形成了一個明顯的分支，表明克里曼丁橘應是一個遺傳背景較為復雜的寬皮柑橘類型，其參與了許多現(xiàn)代寬皮柑橘品種的形成，如費爾柴爾德、弗來蒙特橘、福瓊橘、凱旋柑、秋輝、苔絲、南香等品種中都具有該品種的基因滲入。

2.3.3 未知親本材料的鑒定 GBS技術不僅對研究寬皮柑橘的親緣關系具有重要價值，同時也為未知材料親本的鑒定提供了重要線索。在本研究中加入的56號（1-1）和58號（1-2）兩個生物學重復，其聚類結果顯示兩份重復材料均能聚在一起，表明通過 GBS得到的SNP數(shù)據(jù)和最后的分類結果均是穩(wěn)定可靠的。并且56號和58號兩份材料與愛媛28號和天草的遺傳距離都很近，這兩份材料本身來源于愛媛28號的雜交后代，而愛媛28號又是天草和南香的雜交后代，由此可見，利用GBS數(shù)據(jù)可以準確推導出這些材料的直接親本。另外，本研究也對沃柑的親本進行了探索，以往沃柑被認為是坦普爾和丹西紅橘的雜種后代[29]，但根據(jù)葉片、果實等形態(tài)的鑒定結果表明沃柑不像是從丹西紅橘衍生而來。本研究利用GBS技術發(fā)現(xiàn)沃柑與金諾的親緣關系十分密切，遺傳距離遠低于金諾與韋爾金的遺傳距離，而金諾和韋爾金屬于姊妹系，都是王柑和柳葉橘的雜交后代，因此，沃柑極有可能是王柑和柳葉橘的雜種后代，或者就是金諾的雜交后代。這與國外的報道是一致的，根據(jù)第12屆國際柑橘學會會議報道，利用SSR技術也證明了沃柑是金諾的雜交后代[30]。

同時 GBS技術對了解種質材料的遺傳背景具有重要意義。盛田溫州蜜柑具有果皮極光滑、不易浮皮等特點，以往被認為是宮川溫州蜜柑的早生芽變材料，但從本研究結果來看盛田溫州蜜柑與宮川溫州蜜柑的遺傳距離較遠，不像是宮川溫州蜜柑的芽變材料，極有可能是種子的實生繁殖后代。而大浦溫州蜜柑以往被認為是山崎早生溫州蜜柑芽變而來，但研究結果表明本試驗的大浦溫州蜜柑材料應為雜種來源，其芽變來源的可能性較小。

3 討論

中國是世界寬皮柑橘的起源演化中心，野生寬皮柑橘類型多樣，道縣野橘的發(fā)現(xiàn)表明南嶺山脈是南方寬皮柑橘的起源中心[4,31]，本研究結果進一步表明野生寬皮柑橘與南方地區(qū)的酸橘存在緊密聯(lián)系，可能從其中演化出了酸橘，進而再演化出了砂糖橘、馬水橘等現(xiàn)代栽培的寬皮柑橘品種。

寬皮柑橘由于栽培歷史久遠、地理分布廣泛、種間易發(fā)生雜交，因而遺傳類型極其多樣，對其劃分存在諸多爭議。美國的斯文格認為寬皮柑橘僅僅是柑橘屬下的一個基本種[32]，日本的田中長三郎將寬皮柑橘對應的蜜柑區(qū)進行了較為細致的劃分，包括5個亞區(qū)共 36個種[33]，而中國的曾勉教授根據(jù)植物學特性結合地理分布將寬皮柑橘分為柑和橘兩個亞區(qū)[34]，認為柑亞區(qū)為雜種起源?，F(xiàn)代分子生物學證據(jù)已經(jīng)表明，寬皮柑橘是柑橘屬下的一個基本種[35]，本研究支持該結論。從240份寬皮柑橘的群體進化樹分析來看，將寬皮柑橘劃分為一個種是合理的，但如何對諸多的寬皮柑橘資源進行有效管理，開展更細致的劃分也是必要的。本研究中的系統(tǒng)演化和主成分分析都表明，溫州蜜柑、野生寬皮柑橘、椪柑、克里曼丁橘彼此間存在明顯的遺傳差異，溫州蜜柑、克里曼丁橘是具有明顯雜種來源的寬皮柑橘類型，而椪柑是一個栽培寬皮柑橘類型，斯文格將其作為寬皮柑橘的代表種是不恰當?shù)?，原產(chǎn)中國南方地區(qū)的原始野生寬皮柑橘才真正具有寬皮柑橘代表種的地位，其可能衍生出了酸橘，進而衍生出了現(xiàn)代栽培寬皮柑橘砂糖橘等類型。

利用 GBS技術對寬皮柑橘進行更細致的分類研究，與用植物形態(tài)學分類的結果是吻合的，因而傳統(tǒng)的寬皮柑橘按照植物學特征進行分類劃分具有一定的合理性[36]，比如本研究中的酸橘、椪柑、紅橘、朱紅橘、地中海橘、克里曼丁橘等亞群都能準確地得以劃分。但由于寬皮柑橘與橙類、柚類的種間雜交易于發(fā)生，新的種間雜交還將不斷促進新類型的產(chǎn)生，因此，利用 GBS技術對柑橘資源進行基因分型有助于材料的準確劃分和科學管理。

GBS研究結果還表明，地理隔離、雜交育種、人工選擇是造成資源特異性和多樣性的內在動力，來源于不同地理區(qū)域的資源材料往往會帶有明顯的遺傳印記，這與育種習慣、人工選擇的偏好性有一定關系。人工的長期栽培和選擇使得一些品種群形成明顯類群，如椪柑、紅橘、朱紅橘、地中海橘、克里曼丁橘等栽培品種都形成了明顯的亞群。同時骨干親本在育種中的應用也促進了品種群的形成，比如歐美地區(qū)常利用王柑、默科特等為骨干親本開展雜交，其后代品種往往具有晚熟、皮較硬的特征；而東亞地區(qū)的日本、韓國等常采用清見和椪柑進行雜交，獲得一些肉質細嫩化渣、高糖的品種，進而形成了有明顯遺傳特征的亞群。這些都在本研究中得到充分的證實。

4 結論

本研究利用GBS技術對240份寬皮柑橘資源進行了簡化基因組的基因分型，建立的系統(tǒng)演化樹可以準確地對240份寬皮柑橘進行劃分，其結果與用植物形態(tài)學的劃分結果存在高度吻合，表明GBS技術可高效用于資源材料的準確鑒定和親緣關系的研究，可作為資源材料分子鑒定的重要手段?；贕BS技術具有高通量、低成本等特點，今后可用GBS技術建立涵蓋所有在圃柑橘種質資源的簡化基因組的基因分型數(shù)據(jù)，為今后品種的鑒別和品種權保護提供依據(jù)。

[1] 賀善文. 柑橘類種質資源中心問題的初步探討. 園藝學報, 1979,6(1): 19-25. HE S W. A preminary study of the native citrus in central China. Acta Horticulturae Sinica, 1979, 6(1): 19-25. (in Chinese)

[2] 周志欽. 柑橘分類研究進展一文獻綜述. 園藝學報, 1993, 20(3): 243-250.

ZHOU Z Q. Advances in citrus taxonomy: Literature review. Acta Horticulturae Sinica, 1993, 20(3): 243-250. (in Chinese)

[3] 儲春榮, 江東, 高恒錦, 陳紹彬, 周坤杰. 蘇州地區(qū)寬皮柑橘遺傳多樣性分析. 中國南方果樹, 2016, 45(3): 1-8. CHU C R, JIANG D, GAO H J, CHEN S B, ZHOU K J. Genetic diversity of Mandarin germplasm in Suzhou district. South China Fruits, 2016, 45(3): 1-8. (in Chinese)

[4] 高恒錦, 儲春榮, 王小柯, 陳紹彬, 晏承泉, 閆樹堂. 45份寬皮柑橘野生和地方資源遺傳多樣性分析. 中國南方果樹, 2016, 45(2): 1-9.

GAO H J, CHU C R, WANG X K,CHEN S B, YAN C Q, YAN S T. The genetic diversity of landrace mandarin germplasms. South China Fruits, 2016, 45(2): 1-9. (in Chinese)

[5] WU G A, PROCHNIK S, JENKINS J, SALSE J, HELLSTEN U, MURAT F, PERRIER X, RUIZ M, SCALABRIN S, TEROL J, et al. Sequencing of diverse mandarin, pummelo and orange genomes reveals complex history of admixture during citrus domestication. Nature Biotechnology, 2014, 32(7): 656-662.

[6] ELSHIRE R J, GLAUBITZ J C, SUN Q, POLAND J A, KAWAMOTO K, BUCKLER E S, MITCHELL S E. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE, 2011, 6(5): e19379.

[7] SONAH H, BASTIEN M, IQUIRA E, TARDIVEL A, LEGARE G, BOYLE B, NORMANDEAU E, LAROCHE J, LAROSE S, JEAN M, BELZILE F. An improved genotyping by sequencing (GBS) approach offering increased versatility and efficiency of SNP discovery and genotyping. PLoS ONE, 2013, 8(1): e54603.

[8] GLAUBITZ J C, CASSTEVENS T M, LU F, HARRIMAN, J., ELSHIRE R J, SUN Q. TASSEL-GBS: A high capacity genotyping by sequencing analysis pipeline. PLoS ONE, 2014, 9(2): e90346.

[9] POLAND J A, RIFE T W. Genotyping-by-sequencing for plant breeding and genetics. Plant Genome, 2012, 5(3): 92-102.

[10] POLAND J, ENDELMAN J, DAWSON J, RUTKOSKI J, WU S Y, MANES Y, DREISIGACKER S, CROSSA J, SANCHEZ-VILLEDA H, SORRELLS M, JANNINK J L. Genomic selection in wheat breeding using genotyping-by-sequencing. Plant Genome, 2012, 5(3): 103-113.

[11] POLAND J A, BROWN P J, SORRELLS M E, JANNINK J L. Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLoS ONE, 2012, 7(2): e32253.

[12] WARD J A, BHANGOO J, FERNANDEZ-FERNANDEZ F, MOORE P, SWANSON J D, VIOLA R, VELASCO R, BASSIL N, WEBER C A, SARGENT D J. Saturated linkage map construction in Rubus idaeus using genotyping by sequencing and genome-independent imputation. BMC Genomics, 2013, 14(1): 2.

[13] SPINDEL J, WRIGHT M, CHEN C, COBB J, GAGE J, HARRINGTON S, LORIEUX M, AHMADI N, MCCOUCH S. Bridging the genotyping gap: using genotyping by sequencing (GBS) to add high-density SNP markers and new value to traditional bi-parental mapping and breeding populations. Theoretical and Applied Genetics, 2013, 126(11): 2699-2716.

[14] LU F, ROMAY M C, GLAUBITZ J C, BRADBURY P J, ELSHIRE R J, WANG T, LI Y, LI Y, SEMAGN K, ZHANG X, HERNANDEZ A G, MIKEL M A, SOIFER I, BARAD O, BUCKLER E S. Highresolution genetic mapping of maize pan-genome sequence anchors. Nature Communications, 2015, 6: 6914.

[15] ZHOU Z, ZHANG C, ZHOU Y, HAO Z, WANG Z, ZENG X, DI H, LI M, ZHANG D, YONG H, ZHANG S, WENG J, LI X. Genetic dissection of maize plant architecture with an ultra-high density bin map based on recombinant inbred lines. BMC Genomics, 2016, 17: 178.

[16] ROMAY M C, MILLARD M J, GLAUBITZ J C, PEIFFER J A, SWARTS K L, CASSTEVENS T M, ELSHIRE R J, ACHARYA C B, MITCHELL S E, FLINT-GARCIA S A, MCMULLEN M D, HOLLAND J B, BUCKLER E S, GARDNER C A. Comprehensive genotyping of the USA national maize inbred seed bank. Genome Biology, 2013, 14(6): R55.

[17] MORRIS G P, RAMU P, DESHPANDE S P, HASH C T, SHAH T, UPADHYAYA H D, RIERA-LIZARAZU O, BROWN P J, ACHARYA C B, MITCHELL S E, HARRIMAN J, GLAUBITZ J C, BUCKLER E S, KRESOVICH S. Population genomic and genome-wide association studies of agroclimatic traits in sorghum. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(2): 453-458.

[18] LI X, LI X, FRIDMAN E, TESSO T T, YU J. Dissecting repulsion linkage in the dwarfing gene Dw3 region for sorghum plant height provides insights into heterosis. Proceedings of the National Academy of Sciences of the United States of America, 2015, 112(38): 11823-11828.

[19] GARDNER K M, BROWN P, COOKE T F, CANN S, COSTA F,BUSTAMANTE C, VELASCO R, TROGGIO M, MYLES S. Fast and cost-effective genetic mapping in apple using next-generation sequencing. G3-Genes Genomes Genetics, 2014, 4(9): 1681-1687.

[20] GUAJARDO V, SOLIS S, SAGREDO B, GAINZA F, MUNOZ C, GASIC K, HINRICHSEN P. Construction of high density sweet cherry (Prunus avium L.) linkage maps using microsatellite markers and SNPs detected by genotyping-by-sequencing (GBS). PLoS ONE, 2015, 10(5): e0127750.

[21] LIN M, CAI S, WANG S, LIU S, ZHANG G, BAI G. Genotypingby-sequencing (GBS) identified SNP tightly linked to QTL for pre-harvest sprouting resistance. Theoretical and Applied Genetics, 2015, 128(7): 1385-1395.

[22] LU F, LIPKA A E, GLAUBITZ J, ELSHIRE R, CHERNEY J H, CASLER M D, BUCKLER E S, COSTICH D E. Switchgrass genomic diversity, ploidy, and evolution: Novel insights from a network-based SNP discovery protocol. PLoS Genetics, 2013, 9(1): e1003215.

[23] BAJAJ D, DAS S, BADONI S, KUMAR V, SINGH M, BANSAL K C, TYAGI A K, PARIDA S K. Genome-wide high-throughput SNP discovery and genotyping for understanding natural (functional) allelic diversity and domestication patterns in wild chickpea. Scientific Reports, 2015, 5: 12468.

[24] WONG M M, GUJARIA-VERMA N, RAMSAY L, YUAN H Y, CARON C, DIAPARI M, VANDENBERG A, BETT K E. Classification and characterization of species within the genus lens using genotyping-by-sequencing (GBS). PLoS ONE, 2015, 10(3): e0122025.

[25] WU B, ZHONG G Y, YUE J Q, YANG R T, LI C, LI Y J, ZHONG Y, WANG X, JIANG B, ZENG J W, ZHANG L, YAN S T, BEI X J, ZHOU D G. Identification of pummelo cultivars by using a panel of 25 selected SNPs and 12 DNA segments. PLoS ONE, 2014, 9(4): e94506.

[26] LOMBARDI M, MATERNE M, COGAN N O, RODDA M, DAETWYLER H D, SLATER A T, FORSTER J W, KAUR S. Assessment of genetic variation within a global collection of lentil (Lens culinaris Medik.) cultivars and landraces using SNP markers. BMC Genetics, 2014, 15: 150.

[27] CABEZAS J A, IBANEZ J, LIJAVETZKY D, VELEZ D, BRAVO G, RODRIGUEZ V, CARRENO I, JERMAKOW A M, CARRENO J, RUIZ-GARCIA L, THOMAS M R, MARTINEZ-ZAPATER J M. A 48 SNP set for grapevine cultivar identification. BMC Plant Biology, 2011, 11: 153.

[28] 黎裕, 李英慧, 楊慶文, 張錦鵬, 張金梅, 邱麗娟. 基于基因組學的作物種質資源研究: 現(xiàn)狀與展望. 中國農(nóng)業(yè)科學, 2015, 48(17): 3333-3353.

LI Y, LI Y H, YANG Q W, ZHANG J P, ZHANG J M, QIU L J. Genomics-based crop germplasm research: Advances and perspectives. Scientia Agricultura Sinica, 2015, 48(17): 3333-3353. (in Chinese)

[29] 江東, 曹立. 晚熟高糖雜柑品種‘沃柑’在重慶的引種表現(xiàn). 中國南方果樹, 2011, 40(5): 33-34.

JIANG D, CAO L. The phenotype of ‘Or’ (late-maturing high-sugar hybrid) after cultivating in Chongqing. South China Fruits, 2011, 40(5): 33-34. (in Chinese)

[30] BARRY G H, JR F G G, CHEN C, ROOSE M L, FEDERICI C T, MCCOLLUM G T. Investigating the parentage of ‘orri’ and ‘fortune’mandarin hybrids. Acta Horticulturae, 2015, 36(1065): 449-456.

[31] 廖振坤, 張秋明, 劉衛(wèi)國, 丁偉平, 張玲. 南嶺山脈寬皮柑橘近緣野生種親緣關系鑒定. 湖南農(nóng)業(yè)大學學報(自然科學版), 2006, 32(4): 385-388.

LIAO Z K, ZHANG Q M, LIU W G, DING W P, ZHANG L. Identification of relative relationships of wild relatives of eucitrus originated from Nanling mountains by AFLP analysis. Journal of Hunan Agricultural University (Natural Sciences), 2006, 32(4): 385-388. (in Chinese)

[32] SWINGLE W T. The botany of Citrus and its wild relatives. Citrus Industry, 1967: 190-430.

[33] TANAKA T. Fundamental discussion of Citrus classification. Studia Citrologica, 1977, 14: 1-6.

[34] 曾勉. 對柑橘分類的認識體會和整理的意見. 中國果樹, 1962(2): 31-37.

ZENG M. The understanding and opinion of the citrus classification. China Fruits, 1962(2): 31-37. (in Chinese)

[35] 謝讓金, 周志欽, 鄧烈. 真正柑橘果樹類植物基于AFLP分子標記的分類與進化研究. 植物分類學報, 2008, 46(5): 682-691.

XIE R J, ZHOU Z Q, DENG L. Taxonomic and phylogenetic relationships among the genera of the True Citrus Fruit Trees Group (Aurantioideae, Rutaceae) based on AFLP markers. Journal of Systematics and Evolution (formerly Acta Phytotaxonomica Sinica), 2008, 46(5): 682-691. (in Chinese)

[36] 周開隆, 葉蔭民. 中國果樹志·柑橘卷. 北京: 中國林業(yè)出版社, 2010.

ZHOU K L, YE Y M. China Fruit’s Monograph: Citrus Volume. Beijing: China Forestry Publishing Press, 2010. (in Chinese)

（責任編輯 趙伶俐）

Study on Phylogeny of 240 Mandarin Accessions with Genotyping-by-Sequencing Technology

WANG XiaoKe1, JIANG Dong1,2, SUN ZhenZhu1
(1Citrus Research Institute, Southwest University, Chongqing 400712;2Citrus Research Institute of Chinese Academy of Agricultural Sciences, Chongqing 400712)

【Objective】Genotyping-by-sequencing (GBS) is an economic technique to discover SNP and genotype myriad of crop germplasms in an effective way. The aim of this study is to clarify the classification and evolution of 240 mandarin germplasms by using GBS. 【Method】 A total of 240 mandarin germplasms conserved in the National Citrus Germplasm Repository in Chongqing, with widely genetic diversity and geographic origin, were selected as trial materials. GBS library was constructed withgenomic DNAs after digested with EcoR I restriction endonuclease and sequenced on Illumina HiSeq PE150, then the sequences were mapped to the clementine (Citrus clementina hort. ex Tanaka) reference genome by using BWA, and SNPs were called with the SAMTOOLS pipeline. With the SNPs genotyping data, a phylogenetic tree was built by using Neighbor-joining method and a principal component analysis (PCA) was carried out. 【Result】By using GBS, a total of 96.3 Gb of sequences were generated from the 240 mandarin germplasms, and each sample produced 401.26 Mb in average. After screening with parameter of dp4, miss0.2 and minor alleles frequency (MAF)＞0.01, a total of high quality 114 200 SNP sites were retained. The PCA analysis showed that the 240 mandarin germplasms could be divided into 4 groups, in which satsuma sub-group and wild mandarin sub-group could be clearly separated from other mandarin accessions. With phylogenetic analysis, the 240 mandarin germplasms could be divided into 11 groups. Both the phylogenetic analysis and PCA suggested that the genetic variations were presented in mandarin germplasms with different geographical origins and morphological characteristics. For example, satsuma mandarin (Citrus unshiu Macf.) derived from Japan, clementine (Citrus clementina Hort.ex. Tanaka) and its offspring from Europe and America, as well as wild mandarins from China could be clearly distinguished based on phylogenetic tree, moreover the phylogenetic tree showed that the mandarin germplasms derived from the Southern and Northern of China have unique evolutionary route. The wild mandarins distributed in Nanling Mountain and southern China have a closer phylogenetic relationship with sour mandarin (Citrus sunki Hort.ex Tanaka) and Shatangju mandarin, which are widely cultivated in southern China nowadays, whereas the evolution of mandarins in the northern of China were related to some primitive and old cultivars. In addition, hybrid breeding, long-term artificial selection and domestication led to the subdivision formation and increased the genetic diversity of mandarin. Besides, results of this study showed that GBS has a potential advantage to identify some mandarin accessions with suspicious parents. For example, the phylogenetic tree clearly shows that “Or” tangor has a very close relationship with Kinnow mandarin. 【Conclusion】GBS technology provides an effective and high feasible approach to assistant the taxonomic classification of 240 mandarin accessions, the classification results are in accordance with the conclusion based on morphological method. Meanwhile GBS also is a powerful tool for germplasm identification, and can be applied in the new cultivars identification and intellectual property protection.

GBS; mandarin; phylogenetic evolution; germplasm identification

2016-11-16；接受日期：2017-02-10

國家“十二五”科技支撐計劃（2013BAD01B04）、重慶市科委重點項目（cstc2016shms-ztzx80004）

聯(lián)系方式：王小柯，Tel：18375638987；E-mail：wangxiaoke9191@163.com。通信作者江東，Tel：13983194771；E-mail：jiangdong@cric.cn