基于COS Marker分析柑橘屬及其近緣、遠緣屬植物的遺傳與進化

王炯，龔桂芝，彭祝春，李一兵，王艷杰，洪棋斌

（西南大學柑桔研究所/中國農業科學院柑桔研究所/國家柑桔工程技術研究中心，重慶 400712）

基于COS Marker分析柑橘屬及其近緣、遠緣屬植物的遺傳與進化

王炯，龔桂芝，彭祝春，李一兵，王艷杰，洪棋斌

（西南大學柑桔研究所/中國農業科學院柑桔研究所/國家柑桔工程技術研究中心，重慶 400712）

【目的】研究柑橘屬及其近緣、遠緣屬種質資源的遺傳多樣性及親緣關系，為柑橘類果樹的遺傳和進化關系提供新的視點，也為遺傳育種以及種質資源深入的收集、保護及評價利用提供新工具和參考。【方法】下載Phytozome網站公布的全部克里邁丁Conserved ortholog sequences，并在NCBI中與測序甜橙序列比對，獲得具有Gap的序列，從中設計60對COS引物，篩選并利用多態性標記對柑橘屬及其近緣、遠緣屬植物進行擴增檢測，通過Structure、PowerMarker和GenAlEx等軟件進行數據分析。【結果】從設計引物中篩選出譜帶清晰、穩定性好、具有多態性且分布在柑橘不同連鎖群的引物25對，對45份柑橘屬及其近緣、遠緣屬種質資源進行擴增檢測，共得到584條譜帶，平均每對引物約為23.36條；基因多樣性為0.26—0.88，平均0.49；PIC為0.26—0.87，平均0.48。群體遺傳結構分析顯示，在K=9時，群體結構圖清晰地揭示了柑橘屬及其近、遠緣屬植物間的遺傳組成關系，可區分為寬皮柑橘類群、柚類群、枸櫞類群、大翼橙類群、金柑類群、枳及澳沙檬澳圓檬類群、花椒類群、黃皮酒餅簕九里香類群以及印度野橘和莽山野柑雜合類群，類群關系與聚類分析結果總體比較一致。UPGMA聚類分析顯示，COS Marker能夠明確地區分柑橘屬及其近緣、遠緣屬種質資源，柑橘屬在較高分類層次與金柑屬、枳屬、澳沙檬屬、澳圓檬屬、花椒屬以及黃皮、酒餅簕、九里香等近、遠緣植物分離；在柑橘屬內，枸櫞和大翼橙類首先聚類，然后與寬皮柑橘、柚類等其他柑橘屬植物聚類；寬皮柑橘類群里，道縣野橘、岑溪酸橘等半野生寬皮柑橘與椪柑等栽培寬皮柑橘可以明確區分，而莽山野柑和印度野橘不在本類群。【結論】根據甜橙和克里邁丁基因測序信息設計的COS引物，在柑橘近緣、遠緣屬植物中均能獲得有效擴增，能夠有效區分柑橘屬及其近緣、遠緣屬種質資源；莽山野柑和印度野橘的遺傳組成表現雜合，聚類分析中與寬皮柑橘距離比較遠，可能并非寬皮橘類群最原始類型。

柑橘；COS Marker；親緣關系；種質資源；遺傳與進化

0 引言

【研究意義】蕓香科有約150屬900余種植物，進化層次差異大，地域分布范圍廣，既有栽培價值高、種類豐富、進化程度高的柑橘屬，也有種類相對單一的枳屬、澳沙檬屬等，以及進化層次較低的黃皮、九里香等。由于無融合生殖、屬種間易于雜交、無性變異、栽培歷史悠久等影響，學術界對柑橘的分類及系統演化關系并未形成被廣泛接受的結論。通過開發新標記，開展柑橘屬及其近、遠緣屬種質資源的遺傳多樣性研究，闡明其親緣及進化關系，將為柑橘類果樹的遺傳和進化關系提供新的視點，也為遺傳育種以及種質資源深入的收集、保護及評價利用提供新工具和有價值的參考。【前人研究進展】SWINGLE和REECE[1]根據雄蕊的數目和汁胞構造將多蕊橘屬、澳沙檬屬、澳指檬屬、枳屬、金柑屬以及柑橘屬等6屬植物歸為真正柑橘類果樹，并將柑橘屬劃分為16個種，形成了一個比較有影響力的柑橘分類系統。BARRETT和RHODES[2]采用更多的形態性狀進行數值分析，提出SWINGLE分類的柑橘亞屬中只有柚、枸櫞和寬皮柑橘3個真正生物學種，其余都是雜種起源或后生變異類型。分子標記對研究柑橘及近緣屬親緣關系提供了新的技術和視角。隨著分子生物學的快速發展，不同的分子標記技術不斷開發出來，并被廣泛應用到柑橘及近緣屬的分類和種質資源鑒定等研究中。RFLP、RAPD、AFLP[3-4]、SSR[5-6]、ISSR[7]、SRAP[8]等標記在柑橘的分類、種質資源鑒定和遺傳多樣性分析中得到應用。分子標記在研究柑橘及近緣屬親緣關系中，得到了與BARRETT和RHODES[2]性狀數值分類法相一致的結論，如柑橘屬由柚類、枸櫞和寬皮橘3個種構成，但也提出了許多不一致的論斷。例如，當使用不同的分子標記和聚類策略時，金柑屬與柑橘屬的聚類結果有很大的變化，有支持金柑屬獨立地位的[4,9]，也有認為金柑屬從屬于柑橘屬的[3,10]。導致不一致的原因，既與研究采用材料的代表性和全面性有關，也受到所采用的分子標記本身的影響。RFLP標記為共顯性標記，但操作復雜、費用高，難以對眾多代表性資源類型獲得大量廣泛覆蓋的標記；AFLP標記為顯性標記，有成簇聚集在基因組部分區域的傾向；RAPD理論上能夠覆蓋全部基因組，但其同樣為顯性標記，穩定性和產物的檢測和解讀存在不足；SSR標記比較高效穩定，但在遺傳背景差異大時，難以確定其差異的系譜來源的一致性，進而影響對親緣關系較遠的資源類型的相互關系的準確判斷。COS Marker 是基于保守直序同源序列Conserved ortholog sequence（COS）開發的一種比較新型的分子標記。COS概念首次由FULTON等[11]提出，通過對擬南芥基因組序列和番茄EST數據庫進行計算比對，共鑒別出大約130 000個長度大于 1 000 bp的 EST序列，將其命名為Conserved ortholog sequence（COS），據此設計COS引物并成功應用于茄科植物的圖譜比較。【本研究切入點】COS是高度保守且拷貝數低的基因序列，以PCR為基礎的COS Marker為親緣關系快速鑒定、遺傳圖譜構建和基因組比較等提供了有效的工具，已在不同植物中得到迅速應用，如松屬[12-13]、菊屬[14]、谷屬[15]、煙草屬[16]、薔薇屬[17]和蕓苔屬[18]等。而COS Marker在柑橘中的應用尚未見報道。【擬解決的關鍵問題】基于公布的柑橘基因組測序結果，開發一套分布于柑橘不同染色體的COS Marker，并將其運用到真正柑橘果樹類的 5個近緣屬植物和花椒、黃皮等5個遠緣屬植物的分類和進化關系研究。

1 材料與方法

試驗于 2015年在中國農業科學院柑桔研究所進行。

1.1 試驗材料

45份柑橘屬及其近、遠緣屬材料（表1）均來源于國家果樹種質重慶柑桔圃。

表1 用于COS Marker分析的45份材料Table 1 Forty-five Plant materials for COS Markers analysis

1.2 試驗方法

1.2.1 總 DNA的提取與檢測 采取試驗材料的新鮮幼嫩葉片帶回實驗室，清水洗凈表面，去除葉片主葉脈后稱取1.0 g，采用CTAB法提取DNA，參照文獻[5]并對其方法略做改動：研缽中加入少量PVPP與試驗材料一起研磨；CTAB溶液中加入體積分數為0.2%的β-巰基乙醇；體積比為24∶1的氯仿/異戊醇抽提2次，第1次抽提時加入5 mol·L-1的醋酸鉀，冰上放置5 min后10 000 r/min離心10 min；預冷的異丙醇沉淀DNA；1/10體積4 mol·L-1醋酸鈉和2倍體積預冷的無水乙醇沉淀DNA。200 μL超純水溶解DNA，經NanoDrop 2000微量紫外分光光度計檢測DNA濃度和 1%瓊脂糖凝膠電泳檢測譜帶完整度后，配制擴增檢測的工作模板。

1.2.2 COS引物的開發 從Phytozome v10.0網站主頁的Citrus clementina v1.0中收集下載全部已知測序克里邁丁保守直序同源序列。為了提高所設計開發引物的多態性，在 NCBI（National Center for Biotechnology Information）上利用nucleotide blast將全部下載的克里邁丁保守直序同源序列與已知測序甜橙序列進行比對，獲得在啟動子區或編碼區內克里邁丁和甜橙具有一個或多個 Gap的保守直序同源序列275條。按目的序列應滿足要求：Gap位置應處于目的序列中；均勻分布在不同染色體上，且平均每條染色體上有6—7條目的序列，運用NCBI網站的Pick Primers來設計所有滿足目的序列要求的引物，且所設計開發的引物應盡量保證其具有唯一性，從275條具Gap的序列中挑選并成功設計60對引物。

按照引物類型COS，染色體（Chromosome，C）編號以及在該染色體上的設計的引物序號進行引物命名。1.2.3 COS引物的篩選及產物的檢測 開發的60對COS引物均委托英濰捷基（上海）貿易有限公司進行合成。應用琯溪蜜柚、8號椪柑、梨橙2號、檸檬、金豆、飛龍枳、九里香、酒餅簕等8種代表性材料進行COS引物篩選，從中篩選出譜帶清晰、穩定性好、具有多態性且分布在不同連鎖群的25對COS引物。PCR反應體系為優化的 15 μL反應體系：1×PCR Buffer，1.5 mmol·L-1的Mg2+，0.2 mmol·L-1的dNTPs， 0.33 nmol·L-1的上、下游引物，1 U的Taq DNA聚合酶，DNA模板約75 ng。

PCR擴增反應在Biometra PCR儀上進行。擴增程序: 94℃預變性4 min；94℃變性45 s，64℃退火40 s，72℃延伸45 s，每次循環降低0.5℃，10個循環；94℃變性45 s，58℃退火40 s，72℃延伸 40 s，22個循環；72℃延伸5 min；最后PCR產物4℃保存。

取PCR擴增產物2 μL上樣，8%非變性聚丙烯酰胺凝膠、恒定電壓140 V，電泳約2 h（根據目的片段大小，電泳時間長短略有調整），銀染檢測并照相。1.2.4 數據處理 根據PCR擴增后電泳譜帶的有無及清晰程度，按照二倍體基因型格式讀取譜帶，之后根據相關分析軟件的要求作相應格式的轉換。

應用Structure軟件基于數學模型進行群體結構分析。設置MCMC（Markov Chain Monte Carlo）開始時內循環burn-in迭代為10 000次，再將burn-in迭代后的MCMC設為20 000次，亞群數量（K）設置為2—10，每個K值運行10次，運行完畢后在每個K值中取Ln P(D)最大和最小兩個數值，分別繪制10個最大值和10個最小值的折線圖，找出拐點，并參考系統聚類分析結果來判定合適的亞群數，選出合適的K值。當無明顯拐點時，參照并修改 EVANNO等[19]、強海平等[20]、劉志齋等[21]和GARCIA等[22]的△K算法選取合適的K值。

采用PowerMarker軟件進行聚類分析，計算等位基因（Allele）、基因多樣性（Gene Diversity）、多態信息含量（polymorphism information content，PIC）。采用Nei’ 1983方法計算遺傳距離、UPGMA法聚類，并對結果500次bootstrap重復抽樣，利用MEGA繪制遺傳聚類圖。

運用GenAlEx軟件進行主坐標軸分析。

圖1 引物COSC1-12的COS Marker擴增結果Fig. 1 The COS Marker amplified profile using the primer COSC1-12

2 結果

圖2 不同K值對應的Ln P(D)Fig. 2 Ln P(D) of different K values

2.1 COS遺傳多樣性及COS 引物的分析鑒定效果

用篩選的引物對45份試驗材料進行擴增，均能擴增出穩定的譜帶，范圍100—700 bp。25對引物獲得72個標記位點，擴增出584條譜帶，平均每對引物約為23.36條。擴增譜帶最多的引物為COSC4-3，46條，COSC3-8擴增譜帶最少，10條。25對COS引物的平均 PIC為 0.48，變化范圍 0.26（COSC3-9）—0.87（COSC3-11），其中有11個COS引物的PIC≥0.50（表2）。針對所研究的45份試驗材料，其平均基因多樣性為0.49，變化范圍處于0.26（COSC3-9）—0.88（COSC3-11），其中COSC3-11的PIC和基因多樣性值最高，分別為0.87和0.88，COSC3-9的PIC和基因多樣性值則最低，均為0.26。表2為25對COS引物的遺傳多樣性統計結果。

根據電泳譜帶結果分析，COSC1-12引物能夠將澳沙檬、澳圓檬區分開（圖1）；COSC1-12和COSC3-7引物結合能夠將所有大翼橙類有效分開；COSC2-8引物能夠有效區分黃皮、酒餅簕、九里香等遠緣屬；COSC2-8和COSC3-7引物結合能夠將5種金柑屬材料分開；COSC4-8引物能夠將枳屬兩種材料區分開；COSC4-9引物能夠將3種柚類區分開；COSC9-8引物能夠將花椒屬的兩個種類分離開。由于寬皮橘材料較多，一兩對引物不能有效區分，需結合多對引物才能達到滿意的分類效果。從聚類分析結果看，所篩選的引物能夠有效區分柑橘屬及近遠緣屬。

2.2 群體遺傳結構分析

將數據按照 Structure軟件數據處理要求進行處理，找出拐點對應的K值，發現在K=6、8、9時明顯出現拐點，參考UPMGA聚類分析結果，當K=9時，所有材料分為9個類群是較為合適的。

圖3 Structure群體遺傳結構分析圖Fig. 3 The population structure map of 45 materials

應用Structure軟件得到的群體遺傳結構圖（圖3）中，將45份材料分為9個類群，分別標記為A、B、C、D、E、F、G、H、I，整體符合目前柑橘屬及其近、遠緣屬的分類。晚白柚（1）、琯溪蜜柚（2）、嶺南沙田柚（3）構成柚類一族（A），其遺傳組成單一；黃皮（41）、酒餅簕（42）、九里香（43）構成遠緣類群（B），黃皮和酒餅簕遺傳組成單一，而九里香與花椒屬（44、45）共有少量基因，顯示出它們之間可能有一定的親緣關系；花椒（44）和野花椒（45）構成外類群（C），遺傳組成單一；野生香櫞（22）、合江佛手（23）、檸檬（24）、塔希提萊檬（25）作為枸櫞類組成一類（D），其中野生香櫞和合江佛手的遺傳組成無雜合現象，而檸檬和來檬則顯示有外緣基因滲入，與其明顯的雜種來源一致[23-24]；安江紅橘等寬皮柑橘和塔羅科血橙（20）、梨橙2號（21）等16個柑橘品種組成最多的一類（E），除甜橙、溫州蜜柑、印度酸橘、岑溪酸橘等明顯的雜種外，其余寬皮柑橘的遺傳組成未表現雜合現象；2個甜橙均顯示約3/4的寬皮橘基因以及1/4的其他基因來源，這與甜橙測序的分析結果比較一致[25]，但關于另1/4的其他基因來源則有不同；紅河大翼橙（26）、元江宜昌橙08-7（27）、

元江宜昌橙08-9（28）、金刀峽宜昌橙（29）、馬蜂柑（30）、馬來西亞大翼橙（31）構成另一個柑橘屬類群（F），除馬來西亞大翼橙外，其余大翼橙的遺傳組成無雜合現象；金豆（32）、羅紋（33）、金彈（34）、羅浮（35）、四季橘（36）作為金柑屬組成一類（G），除四季橘顯示出有外緣基因滲入外，其余材料遺傳組成無雜合現象；莽山野柑（17）和印度野橘（18）單獨成為一類（H），與其他寬皮橘有較大差異，遺傳來源較為獨特，但二者的遺傳組成表現出雜合現象，印度野橘有較少量的枸櫞血緣，莽山野柑有較多的寬皮橘血緣，顯示二者可能有比較獨特的遺傳來源，但并非該類群最原始類型；近緣屬中的旺蒼大葉枳（37）、飛龍枳（38）與澳沙檬（39）、澳圓檬（40）組成一個類群（I），除澳沙檬、澳圓檬顯示極少的外緣基因外，遺傳組成也未表現雜合現象。

表2 COS引物及其揭示的遺傳多樣性統計Table 2 Summary of COS Markers and genetic diversity revealed

2.3 聚類分析

運用PowerMarker軟件對45份試驗材料聚類分析，構建UPGMA聚類樹（圖4）。第一條分割線將材料分成A、B兩個類群：A類群為近遠緣屬植物；B類群為柑橘屬植物。第二條分割線將柑橘屬及其近遠緣屬進一步劃分為C、D、E、F、G、H六個類群：C類群為外緣屬植物，包括黃皮（41）、酒餅簕（42）、九里香（43）、花椒（44）、野花椒（45）；四季橘（36）與其他金柑屬植物聚成D類群，但它與其他金柑屬植物遺傳距離較遠，Structure群體遺傳結構分析顯示四季橘有金柑屬血緣并占主導地位，說明四季橘與金柑屬植物親緣關系較近；E類群為枳屬（37、38）和澳沙檬（39）、澳圓檬（40）構成的近緣屬類群；枸櫞類植物和大翼橙類植物構成F類群；G類群包括興津溫州蜜柑（12）、兩個甜橙（20、21）；H類群為柚類、寬皮柑橘類植物類群。第三條分割線將H類群進一步劃分為H1、H2兩個類群：莽山野柑（17）、印度野橘（18）、柚類（1、2、3）構成H1類群；其余寬皮橘構成H2類群，在H2類群中，又可細分為以安江紅橘（4）、萬州少核紅橘（5）、砂糖橘（9）、土橘（6）等為代表的栽培種和以旺蒼皺皮柑（13）、細皮狗屎柑（14）等為代表的野生、半野生寬皮橘兩類。聚類分析結果顯示，COS Marker能夠明確的區分柑橘屬及其近、遠緣屬種質資源，并對成員眾多、關系復雜的寬皮類群進行有效區分。

2.4 群體主坐標軸分析

采用GenAlEx軟件基于遺傳距離數據對45份材料進行的PCoA分析結果見圖5。可以看到在第一主坐標軸區，顯示綠色方形部分的寬皮柑橘和藍色三角形部分甜橙類塔羅科血橙（20）、梨橙2號（21）關系較近；在第二主坐標軸區，柚類和莽山野柑（17）、印度野橘（18）關系較近，與聚類分析結果一致；第三主坐標軸區的紫色菱形部分的大翼橙類和黃色圓形部分的枸櫞類關系較近，聚類分析也顯示兩者有較近的親緣關系；第四主坐標軸區的淺藍色部分的金柑屬聚在一起，四季橘（36）又遠離于其他金柑屬，這與 Structure群體遺傳結構分析結果、聚類分析結果一致；近遠緣屬的枳屬（37—38）、澳沙檬（39）、澳圓檬（40）、黃皮（41）、酒餅簕（42）、九里香、（43）花椒屬（44—45）等與柑橘屬距離較遠，與聚類分析結果一致。

3 討論

3.1 柑橘及其近、遠緣屬植物的分離分化時間

一個比較明確的時間點對認知植物的分類與進化關系非常有價值。地質事件和化石提供了比較準確的分離分化時間估計[26]。在柑橘屬及其近、遠緣屬植物起源和分布的“桶形”區域內[1]，發生了古大陸的分離和漂移以及印度板塊與亞歐板塊碰撞而帶來喜馬拉雅運動等重大地質事件，對柑橘屬及其近、遠緣屬植物的分離分化和分布產生了重要影響，也為研究其分離分化提供了比較客觀的參考時間點。

本研究所采用的材料，按SWINGLE和REECE[1]的分類，包括了親緣關系很近的5個真正柑橘果樹類植物，也包括了關系較近、與真正柑橘類位于相同亞族（subtribe）的酒餅簕，還有與真正柑橘類位于相同亞科但不同族的黃皮和九里香，以及與真正柑橘類不同亞科的花椒和野花椒。這些材料中，采樣分析的遠緣屬材料均來自中國，但花椒屬在世界范圍內廣泛分布，九里香屬、黃皮屬和酒餅簕局限在中國、東南亞、澳大利亞等地，分布區域與真正柑橘類非常一致；5個近緣屬中，中國是柑橘屬的主要起源地，是金柑屬和枳屬的唯一起源地，而澳大利亞和新幾內亞地區是澳指檬、澳沙檬和多蕊橘3個屬的唯一起源地。由于古大陸的分離和漂移而造成的廣闊海洋分隔，澳洲 3屬與中國的3屬以及九里香等遠緣植物間難以完成遺傳物質的交流，如果它們是同步由共同祖先分離分化并獨立演進形成，反映到分支或聚類關系上，6個柑橘近緣屬植物應該首先屬內聚類，然后澳洲3屬與中國起源3屬各自聚類，再在近緣屬聚類后與其他遠緣屬植物聚類。然而在本研究中各近、遠緣屬卻展示了不同的遺傳一致性和聚類層次。

圖4 基于COS Marker數據構建的UPGMA聚類樹（500次bootstrap檢驗）Fig. 4 An UPGMA tree of 45 different materials derived from bootstrap analysis (500 replications) of COS Marker data

采用COS Marker能夠將柑橘屬及其近、遠緣屬植物進行有效區分，但遺傳結構分析顯示枳屬和澳指檬屬與澳沙檬屬有非常一致的遺傳組成，黃皮、酒餅簕和九里香也顯示高度一致的遺傳組成；UPGMA聚類分析則顯示枳屬和澳指檬屬與澳沙檬屬在屬一級的分類中關系更近，它們先聚類后再與柑橘屬、金柑屬及其他屬構成較高級分類單元，花椒、九里香、酒餅簕和黃皮等遠緣植物彼此在分類上更接近。FEDERICI等[27]也得到類似結果，在其應用RFLP數據進行非雜種來源的材料聚類分析中，顯示枳屬與多蕊橘屬和澳指檬屬的M. warburgiana具有較近關系，它們先聚類后再與澳沙檬及澳指檬類型聚類；在應用RFLP數據進行所有材料聚類分析中，顯示枳屬與澳沙檬及多蕊橘屬和澳指檬屬的M. warburgiana聚類后，再與其他澳指檬類型聚類。

圖5 基于COS Marker數據的主坐標軸分析圖Fig. 5 PCoA plot constructed with COS Marker data

枳屬與3個澳洲近緣屬有更近似的遺傳組成和親緣關系，超越與其有長期共同地域分布的柑橘屬和金柑屬的關系，一個重要的分類學意義就是表明枳屬與澳洲3近緣屬分離分化的時間應該可能晚于它們與柑橘屬和金柑屬的分離分化。真正柑橘類的共同祖先可能是先分化形成柑橘屬原始類型、金柑屬原始類型以及枳屬與澳洲3近緣屬的共同原始類型，然后在長期的生殖隔離作用下，枳屬與澳洲3近緣屬的共同原始類型分別演化，形成了形態和適應生境差異甚大的枳、澳沙檬和多蕊橘等類型。真正柑橘類的最近共同祖先在分化形成不同分支類型時，理論上應該仍處于未完全割裂的陸地環境，花粉、種子或果實擴散能夠順利完成，對應的大地質背景應該是中國所在的歐亞古陸與澳大利亞所在的岡瓦納古陸尚未徹底分離。從古植物學方面看，在歐洲早始新世到晚上新世的第三紀地層中就有大量花椒屬種子化石被發現，在北美早中新世也發現有該屬種子化石，在中國山東臨朐中新世地層曾報道有花椒屬的葉片化石。在廣東和廣西第四紀地層中找到黃檗屬及酒餅簕屬的葉片化石[28]。SWINGLE和REECE[1]曾認為有充足的理由推測，澳沙檬是由澳指檬經過2千萬至3千萬年的漫長干旱適應演化形成，結合澳洲大陸分離的重大地質事件和被子植物出現的時期，認為真正柑橘類植物的產生應該在6 500萬年前的晚白堊紀。

枳屬與澳洲3近緣屬的相似遺傳組成和相近親緣關系，以及地理分布特點和區域內的地質事件為估計真正柑橘類植物的分離分化提供了參考，而印度野橘和莽山野柑則為柑橘屬內的分離分化提供了依據。印度野橘和莽山野柑是兩個地域來源明確、分布區域有限、遺傳特性獨特的野生原始柑橘種。前者分布在印度東北的阿薩姆邦區域，被認為是栽培柑橘的原始類型[29-30]；而后者則分布在中國中部的湖南莽山區域，因其性狀的原始，也被認為是寬皮柑橘的原始類型[31]。本研究中COS Marker能夠將二者進行有效區分，但遺傳結構分析顯示二者有較高的相似性，且相似成分非常獨特，該成分僅在甜橙、檸檬、來檬等明顯雜種來源的其他類型中低比例出現，UPGMA聚類分析也顯示印度野橘和莽山野柑親緣關系比較近，超過與其他寬皮柑橘或枸櫞類的相近程度。分布區域遙遠且其間橫亙了阻斷普通生物交流的喜馬拉雅山的兩個野生柑橘種，有近似的遺傳組成和親緣關系，表明其在起源上可能有獨特的關系。從地理上可知，印度東北的阿薩姆邦等區域實際上位于喜馬拉雅山南麓，本身是中國所在歐亞板塊的一部分，而非印度主體所在的印度板塊的組成部分，只是由于印度板塊向歐亞板塊的伏沖擠壓，強烈的“喜馬拉雅造山運動”，才被喜馬拉雅山分隔開來。印度野橘和莽山野柑可能是在其共同祖先的棲息地被“喜馬拉雅造山運動”分隔后，經過第四紀冰川，才在各自現今的區域獨自演化為迄今的類型。“喜馬拉雅造山運動”達到使普通生物交流困難的時段應該可以作為印度野橘和莽山野柑的分離分化參考時間點，可能在2千萬年至3千萬年前。

3.2 莽山野柑、印度野橘與柑橘的起源

處于相對隔絕的地域，未受或少受人工選擇或栽培柑橘的影響，野生柑橘對柑橘及近緣屬的起源和進化研究具有重要作用。莽山野柑和印度野橘是目前尚存但瀕危的野生柑橘，被認為是栽培柑橘的原始類型[29-31]。廖振坤等[31]根據汁囊結構、花粉壁的外壁雕紋特征比較及AFLP聚類結果顯示，莽山野柑是比道縣野橘更原始的類型，并且認為南嶺山脈存在莽山野柑→道縣野橘→栽培種的寬皮橘的演化路徑。JOSE等[32]研究發現寬皮橘多達38個左右的indels，而莽山野柑僅具有4個indels，并且莽山野柑與宜昌橙之間僅有10%的異質一致性，認為莽山野柑與寬皮橘、宜昌橙親緣關系較遠，但其聚類分析顯示莽山野柑與宜昌橙有較近的親緣關系，王福生等[33]也得到兩者具有較近的親緣關系的結果。而劉勇等[34]的SSR-AFLP聚類分析顯示莽山野柑與宜昌橙和大翼橙親緣關系較遠。MALIK[29]、LASKAR[35]等認為印度野橘是最原始的柑橘類型，是柑橘栽培種的祖先之一。劉通等[8]的研究涉及莽山野柑與印度野橘，通過SSR-SRAP聚類分析顯示，在遺傳相似系數為0.725處，它們和柚類聚為一個分支，表明三者有較近的親緣關系。而FEDERICI[27]認為印度野橘可能是卡西大翼橙或酸橙的雜種，但其RFLP和RAPD標記則顯示印度野橘優先與枸櫞聚類，NICOLOSI等[36]通過cpDNA聚類分析也顯示印度野橘和枸櫞優先聚為一類。

本研究中的遺傳結構分析顯示莽山野柑和印度野橘有較為獨特的相似遺傳組成，但與其他寬皮橘存在較大差異。莽山野柑有近半的寬皮柑橘遺傳組成，且其遺傳組成表現雜合現象。而印度野橘則含有較少量的枸櫞血緣，這與FEDERICI等[27]、NICOLOSI等[36]結果類似。UPGMA聚類分析也顯示印度野橘與莽山野柑優先聚類而遠離其他寬皮橘和枸櫞等類型，與劉通等[8]SSR-SRAP聚類結果一致。因此筆者認為莽山野柑和印度野橘是比較獨特原始的柑橘類型，但不應該是寬皮柑橘的直接祖先。

4 結論

根據甜橙和克里邁丁基因測序信息設計的 COS引物，在柑橘近緣、遠緣屬植物中均能獲得有效擴增，能夠有效區分柑橘屬及其近、遠緣屬種質資源。莽山野柑和印度野橘有獨特的遺傳組成，應該是比較獨特的柑橘類型，其遺傳組成表現雜合現象，聚類分析中與寬皮柑橘距離比較遠，應該不是寬皮柑橘的直接祖先。

[1] SWINGLE W T, REECE P C. The botany of citrus and its wild relatives//REUTHER W, WEBBER H J, BATCHELOR L D. The Citrus Industry. Vol. 1. 2nd ed. Berkley: University of California Press. 1967: 190-422.

[2] BARRETT H C, RHODES A M. A numerical taxonomic study of affinity relationships in cultivated citrus and its close relatives. Systematic Botany, 1976, 1: 105-136.

[3] 謝讓金, 周志欽, 鄧烈. 真正柑橘果樹類植物基于 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)

[4] PANG X M, HU C G, DEND X X. Phylogenetic relationships within citrus and its related genera as inferred from AFLP markers. Genetic Resources and Crop Evolution, 2007, 54: 429-436.

[5] 龔桂芝, 洪棋斌, 彭祝春, 江東, 向素瓊. 枳屬種質遺傳多樣性及其近緣屬植物親緣關系的 SSR和 cpSSR分析. 園藝學報, 2008, 38(12): 1742-1750. GONG G Z, HONG Q B, PENG Z C, JIANG D, XIANG S Q. Genetic diversity of poncirus and its phylogenetic relationships with relatives as revealed by nuclear and chloroplast SSR. Acta Horticulturae Sinica, 2008, 38(12): 1742-1750. (in Chinese)

[6] 雷天剛, 何永睿, 彭愛紅, 許蘭珍, 鄒修平, 姚利曉, 李金強, 陳善春. 貴州部分柑橘地方品種遺傳背景的 SSR 分析. 西南大學學報(自然科學版), 2015, 37(11): 31-35. LEI T G, HE Y R, PENG A H, XU L Z, ZOU X P, YAO L X, LI J Q, CHEN S C. Genetic background analysis of several Guizhou local citrus varieties based on SSR markers. Journal of Southwest University (Natural Science Edition), 2015, 37(11): 31-35. (in Chinese)

[7] 廖芳蕾, 陳民管, 桑丹, 陳文榮, 郭衛東. 佛手種質資源遺傳多樣性的ISSR分析. 園藝學報, 2013, 40(11): 2222-2228. LIAO F L, CHEN M G, SANG D, CHEN W R, GUO W D. Genetic diversity analysis of fingered citron germplasm by ISSR. Acta Horticulturae Sinica, 2013, 40(11): 2222-2228. (in Chinese)

[8] 劉通, 鄧崇嶺, 程玉芳, 李秋景, 陳傳武, 劉冰浩, 伊華林. 利用SSR和 SRAP技術分析廣西柑橘種質遺傳多樣性. 華中農業大學學報, 2016, 35(2): 23-29. LIU T, DENG C L, CHENG Y F, LI Q J, CHEN C W, LIU B H, YI H L. Analyzing genetic diversity of citrus germplasm in Guangxi province with SSR and SRAP markers. Journal of Huazhong Agricultural University, 2016, 35(2): 23-29. (in Chinese)

[9] PENJOR T, ANAI T, NAGANO Y, MATSUMOTO R, YAMAMOTO M. Phylogenetic relationships of citrus and its relatives based on rbcL gene sequences. Tree Genet Genomes, 2010, 6: 931-939.

[10] BARKLEY N A, ROOSE M L, KRUEGER R R, FEDRICI C T. Assessing genetic diversity and population structure in a citrus germplasm collection utilizing simple sequence repeat markers (SSRs). Theoretical and Applied Genetics, 2006, 112: 1519-1531.

[11] FULTON T M, HOEVENR V, EANNETTAN T, TANKSLEY S D. Identifiation, analysis, and utilization of conserved ortholog set markers for comparative genomics in higher plants. Plant Cell, 2002, 14: 1457-1467.

[12] KRUTOVSKY K V, TROGGIO M, BROWN G R, JERMSTAD K D, NEALE D B. Comparative mapping in the Pinaceae. Genetics, 2004, 168(1): 447-461.

[13] LIEWLAKSANEEYANAWIN C, ZHUNG J, TANG M, FARZANEH N, LUENG G, CULLIS C, FINDLAY S, RITLAND C E, BOHLMANN J, RITLAND K. Identification of COS markers in the Pinaceae. Tree Genetics & Genomes, 2008, 5(1): 247-255.

[14] CHAPMAN M A, CHANG J, WEISMAN D, KESSELIR V, BURKE J M. Universal markers for comparative mapping and phylogenetic analysis in the Asteraceae (Compositae). Theoretical and Applied Genetics, 2007, 115: 747-755.

[15] QURAISHI U M, ABROUK M, BOLOT S, PONT C, THROUDE M, GUILHOT N, CONFOLENT C, BORTOLINI F, PRAUD S, MURIGEUX A, CHARMET G, SALSE J. Genomics in cereals:from genome-wide conserved orthologous set (COS) sequences to candidate genes for trait dissection. Functional & Integrative Genomics, 2009, 9: 473-484.

[16] WU F N, NANCY T E, XU Y M, PLIESKE J, GANAL M, POZZI C, BAKAHER N, TANKSLEY S D. COSII genetic maps of two diploid nicotiana species provide a detailed picture of synteny with tomato and insights into chromosome evolution in tetraploid N. tabacum. Theoretical and Applied Genetics, 2010, 120: 809-827.

[17] BUSHAKRA J M, SARGENT D J. CABRERA A, CROWHURST R, GIRONA E L, VELASCO R, SYMONDS V V, KNAAP E, TROGGIO M, GARDINER S E. CHAGNE D. Rosaceae conserved orthologous set (RosCOS) markers as a tool to assess genome synteny between Malus and Fragaria. Tree Genetics & Genomes, 2011, 8: 643-658.

[18] JEONG Y M, CHUNG W H, CHUNG H, NAMSHIN K, PARK B S, LIM K B, YU H J, MUN J H. Comparative analysis of the radish genome based on a conserved ortholog set (COS) of Brassica. Theoretical and Applied Genetics, 2014, 127: 1975-1989.

[19] EVANNO G, REGNAUT S, GOUDET J. Detecting the number of clusters of individuals using the software STRUCTURE:a simulation study. Molecular Ecology, 2005, 14(8): 2611-2620.

[20] 強海平, 余國輝, 劉海泉, 高洪文, 劉貴波, 趙海明, 王贊. 基于SSR標記的中美紫花苜蓿品種遺傳多樣性研究. 中國農業科學, 2014, 47(14): 2853-2862. QIANG H P, YU G H, LIU H Q, GAO H W, LIU G B, ZHAO H M, WANG Z. Genetic diversity and population structure of Chinese and American Alfalfa (Medicago Sativa. L) germplasm assessed by SSR markers. Scientia Agricultura Sinica, 2014, 47(14): 2853-2862.(in Chinese)

[21] 劉志齋, 吳迅, 劉海利, 李永祥, 李清超, 王鳳格, 石云素, 宋燕春,宋偉彬, 趙久然, 賴錦盛, 黎裕,王天宇. 基于40個核心SSR標記揭示的 820份中國玉米重要自交系的遺傳多樣性與群體結構. 中國農業科學, 2012, 45(11): 2107-2138. LIU Z Z, WU X, LIU H L, LI Y X, LI Q C, WANG F G, SHI Y S, SONG Y C, SONG W B, ZHAO J R, LAI J S, LI Y, WANG T Y. Genetic diversity and population structure of important Chinese maize inbred lines revealed by 40 core simple sequence repeats (SSRs). Scientia Agricultura Sinica, 2012, 45(11): 2107-2138. (in Chinese)

[22] GARCIA A L, LURO F, OLLITRAULT P, NAVARRO L. Genetic diversity and population structure analysis of mandarin germplasm by nuclear, chloroplastic and mitochondrial markers. Tree Genetics & Genomes, 2015, 11: 123-137.

[23] 梁國魯. 部分柑桔屬及其近緣屬Giemsa C-帶帶型研究. 遺傳學報, 1988, 15(6): 409-415. LIANG G L. Study on the Giemsa C-banding Patterns of Some Citrus and Its Related Genera. Acta Genetica Sinica, 1988, 15(6): 409-415. (in Chinese)

[24] 熊光明. 應用 AFLP分子標記對柑橘屬(Cirtus)植物進行鑒別與系統分類研究[D]. 重慶: 西南農業大學, 2002. XIONG G M. Study on citrus Identification and Taxonomy by using AFLP markers [D]. Chongqing: Southwest Agricultural University, 2002. (in Chinese)

[25] XU Q, CHEN L L, RUAN X A, CHEN D J, ZHU A D, CHEN C L, BERTRAND D, JIAO W B, HAO B H, MATTHEW P L, CHEN J J, GAO S, XING F, LAN H, CHANG J W, GE X H, LEI Y, HU Q,MIAO Y, WANG L, XIAO S X, BISWAS M K, ZENG W F, GUO F, CAO H B, YANG X M, XU X W, CHENG Y J, XU J, LIU J H, OSCAR, LUO J H, TANG Z H, GUO W W, KUANG H H, ZHANG H Y, MIKEAL L R, NIRANJAN N, DENG X X, RUAN Y J. The draft genome of sweet orange (Citrus sinensis). Nature Genetics, 2013, 45: 59-66.

[26] DONOGHUE P C J, BENTON M J. Rocks and clocks: Calibrating the tree of life using fossils and molecules. Trends in Ecology & Evolution, 2007, 22: 424-431 .

[27] FEDERICI C T, FANG D Q, SCORA R W, ROOSE M L. Phylogenetic relationships within the genus Citrus （Rutaceae）and related genera as revealed by RFLP and RAPD analysis. Theoretical and Applied Genetics, 1998, 96: 812-822.

[28] 吳征鎰, 陳心啟. 中國植物志. 北京: 科學出版社, 1997, 43(2): 1. WU Z Y, CHEN X Q. Flora Republicae Popularis Sinicae. Beijing: Science Press, 1997, 43(2): 1. (in Chinese)

[29] MALIK S K, CHAUDHURY R, DHARIWAL O P, RAJWANT K K. Collection and characterization of Citrus indica Tanaka and C. macroptera Montr: wild endangered species of northeastern India. Genetic Resources and Crop Evolution, 2006, 53: 1485-1493.

[30] TANAKA T. Citrologia semi centennial commemoration papers on Citrus studies. Osaka: Citrologia Supporting Foundation, 1961: 114.

[31] 廖振坤, 張秋明, 劉衛國, 丁偉平, 張玲. 南嶺山脈寬皮柑橘近緣野生種親緣關系鑒定. 湖南農業大學學報(自然科學版), 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 Edition), 2006, 32(4): 385-388. (in Chinese)

[32] JOSE C C, ROBERTO A, VICTORIA I, JAVIER T, MANUEL T, JOAQUIN D. A phylogenetic analysis of 34 chloroplast genomes elucidates the relationships between wild and domestic species within the genus citrus. Molecular Biology and Evolution, 2015, 32(8): 2015-2035.

[33] 王福生, 江東. 應用cpSSR和EST-SSR標記進行柑橘特異種質資源遺傳背景研究. 園藝學報, 2010, 37(3): 465-474. WANG F S, JIANG D. Studies on genetic background of important germplasm resources among citrus based on cpSSR and ESTSSR marker. Acta Horticulturae Sinica, 2010, 37(3): 465-474. (in Chinese)

[34] 劉勇, 孫中海, 劉德春, 吳波, 陶建軍. 柚類種質資源 AFLP與SSR遺傳多樣性分析. 中國農業科學, 2005, 38(11): 2308-2315. LIU Y, SUN Z H, LIU D C, WU B, TAO J J. Assessment of the genetic diversity of pummelo germplasms using AFLP and SSR markers. Scientia Agricultura Sinica, 2005, 38(11): 2308-2315. (in Chinese)

[35] LASKAR M A, HYNNIEWTA M, RAO C S. In vitro propagation of citrus indica Tanaka—An endangered progenitor species. Indian Journal of Biotechnology, 2009, 8: 311-316.

[36] NICOLOSI E, DENG Z N, GENTILE A, MALFA S L, CONTINELLA G, TRIBULATO E. Citrus phylogeny and genetic origin of important species as investigated by molecular markers. Theoretical and Applied Genetics, 2000, 100: 1155-1166.

（責任編輯 趙伶俐）

Genetic and Phylogenetic Relationships Among Citrus and Its Close and Distant Relatives Based on COS Marker

WANG Jiong, GONG GuiZhi, PENG ZhuChun, LI YiBing, WANG YanJie, HONG QiBin
（Citrus Research Institute, Southwest University/Chinese Academy of Agricultural Sciences, National Citrus Engineering Research Center, Chongqing 400712）

【Objective】To reveal genetic diversity and phylogenetic relationships of Citrus and its close and distant relatives with newly developed conserved ortholog sequence (COS) Marker.【Method】According to the published C. sinensis and C. clementina genomic sequencing information, COS Markers were developed and screened. Polymorphic markers were used to amplify and detect Citrus and its close and distant relatives. Structure, Power Marker and GenAlEx were used to analyze the relevantdata.【Result】A total of 60 COS primers were designed, and 25 primer pairs distributed in different linkage groups were selected based on clear and good DNA polymorphic bands. In 45 samples from Citrus and its relatives, 25 COS markers produced 584 bands, with an average of about 23.36 bands per marker, and genetic diversity ranged from 0.26 to 0.88, with an average of 0.49; PIC value was between 0.26 and 0.87, with an average of 0.48. When K = 9, a clear genetic structure relationship of citrus and its relatives were revealed. They were divided into C. reticulata group, C. grandis group, C. medica group, Papeda group, Fortunella group, P. trifoliate and Eremocitrus and M. australis group, Z. bungeanum group, Clausena lansium and Atalantia buxifolia and Murraya paniculata group as well as C. indica and C. mangshanensis hybrid groups. This result was globally consistent with that in the UPGMA clustering analysis. UPGMA clustering analysis revealed that COS Marker could effectively distinguish Citrus from its close and distant germplasm resources. At higher taxonomy level, Citrus was first separated from P. trifoliate, Eremocitrus, M. australis, Z. bungeanum, Clausena lansium, Atalantia buxifolia and Murraya paniculata; In Citrus, C. medica and Papeda group were first clustered, and then they were clustered with C. reticulata and C. grandis and other plants of Citrus. In C. reticulata group, Daoxian yeju and Cenxi suanju and other half wild C. reticulata were clearly distinguished from ponkan and the rest cultivated C. reticulata, but C. mangshanensis and C. indica were not in the group.【Conclusion】COS primers based on Citrus genomic information could get effective amplification and distinguish Citrus from its close and distant relatives effectively. The genetic makeup of C. mangshanensis and C. indica were heterozygous and different from C. reticulata, so they may not be the most primitive types of C. reticulata.

Citrus; COS marker; phylogenetic relationship; germplasm; genetics and evolution

2016-05-06；接受日期：2016-07-15

國家科技支撐計劃（2013BAD02B02）、重慶市基礎與前沿研究計劃（cstc2014jcyjA80037）

聯系方式：王炯，Tel：18375633275；E-mail：18375633275@163.com。通信作者洪棋斌，Tel：023-68349025；E-mail：hongqb@sina.com