雙三倍體銀鯽超數(shù)微小染色體的基因組特征

石 倩 丁 苗 汪 洋 周 莉 王忠衛(wèi) 張曉娟李 志 李熙銀 * 桂建芳 *

(1. 大連海洋大學水產(chǎn)與生命學院, 大連 116023; 2. 中國科學院水生生物研究所淡水生態(tài)與生物技術國家重點實驗室,水產(chǎn)品種創(chuàng)制與高效養(yǎng)殖重點實驗室, 湖北洪山實驗室, 中國科學院種子創(chuàng)新研究院, 武漢 430072;3. 中國科學院大學, 北京 100049)

超數(shù)染色體(又叫超數(shù)B染色體、B染色體或者額外染色體), 是指除了常規(guī)染色體(又叫常規(guī)A染色體或A染色體: 包含常染色體和性染色體)之外的非必須染色體[1]。超數(shù)染色體最早于1907年在半翅目昆蟲中被記載[2], 且預計存在于約15%的真核生物中[3]。由于超數(shù)染色體不遵循孟德爾遺傳定律,所以其數(shù)量在同一群體不同個體之間也會存在差異[4]。超數(shù)染色體通常被認為起源于常規(guī)染色體,并且能夠積累細胞器來源的DNA序列[4]。超數(shù)染色體發(fā)現(xiàn)之初, 被認為不含有關鍵基因且沒有功能,然而, 隨著測序技術和基因功能研究手段的迅速發(fā)展, 越來越多的證據(jù)顯示超數(shù)染色體上的活性基因能夠發(fā)揮重要的生物學功能[5,6], 且為物種的演化提供了額外的基因組材料[7], 例如: 在小麥病原真菌(Zymoseptoria tritici)中, 超數(shù)染色體能增強其對小麥的致病性[8]; 在黑麥(Secale cereale)中, 超數(shù)染色體能夠增加常規(guī)染色體的重組率并增強其對逆境條件的適應性[9,10]; 在玉米(Zea maysL.)中, 超數(shù)染色體能夠增強其對壞死病的抵抗能力[11]; 在洞穴魚(Astyanax mexicanus)中, 超數(shù)染色體含有性別決定基因, 具有雄性決定作用[5]。

多倍體鯽復合種包括有性生殖四倍體鯽(Carassius auratus)和單性雌核生殖六倍體銀鯽(Carassius gibelio), 廣泛分布于亞歐大陸及鄰近島嶼的淡水水系[12,13]。全基因組解析揭示四倍體鯽是雙二倍體(AABB,n=100), 包含兩套二倍體基因組, 每套二倍體基因組源于不同祖先; 六倍體銀鯽是雙三倍體(AAABBB,n≈150), 包含兩套三倍體基因組, 每套三倍體基因組源自不同祖先[14—16]。約82—96萬年前, 雙三倍體銀鯽由祖先雙二倍體鯽經(jīng)同源三倍化而形成, 多倍化后的雙三倍體銀鯽通過雌核生殖克服三個同源染色體不能正常配對和均等分離的生殖障礙[14]。雌核生殖是指母本產(chǎn)生不減數(shù)卵子, 不減數(shù)的卵子在同域有性生殖物種精子的刺激下進行胚胎發(fā)育, 最終產(chǎn)生與母本遺傳背景一致的全雌子代[17,18]。雙三倍體銀鯽通過抑制減數(shù)第一次分裂, 從而產(chǎn)生不減數(shù)的卵子[14], 并利用同水域的雙二倍體鯽或其他物種的精子進行雌核生殖[17,18]。有意思的是, 雌核生殖銀鯽不同于其他單性脊椎動物, 在自然群體中含有少部分雄性個體[19]。

雙二倍體鯽含有100條常規(guī)染色體, 不含有超數(shù)染色體; 而雙三倍體銀鯽含有約150條染色體,除了常規(guī)染色體之外還有一些超數(shù)染色體, 由于銀鯽的超數(shù)染色體大小比常規(guī)染色體偏小, 我們將其稱為超數(shù)微小染色體。銀鯽雌性個體含有約9個超數(shù)微小染色體, 而雄性個體比雌性個體多出了額外的3—4個雄性特異超數(shù)微小染色體, 且這些雄性特異的超數(shù)微小染色體具有遺傳雄性決定作用[20,21]。隨后, 我們通過熒光顯微切割、體外擴增和測序, 對這些雄性特異超數(shù)微小染色體進行了解析, 并篩選到具有雄性特異或者雄性偏向表達的基因片段[20]。然而, 銀鯽雌性中超數(shù)微小染色體的基因組特征及其在單性雌核生殖方式演化中的作用仍不清楚。

在本研究中, 我們利用銀鯽超數(shù)微小染色體富集的重復序列和染色體熒光原位雜交, 鑒定了銀鯽雌性基因組中的潛在超數(shù)微小染色體序列。我們發(fā)現(xiàn)銀鯽超數(shù)微小染色體含有所有常染色體的同源序列、大量重復序列以及完整活性基因。此外,大部分銀鯽中擴張的減數(shù)分裂相關基因家族在超數(shù)微小染色體上存在擴增拷貝。這些結果不僅解析了雙三倍體銀鯽的超數(shù)微小染色體基因組特征,也為超數(shù)染色體在單性生殖演化中的作用提供了創(chuàng)新見解。

1 材料與方法

1.1 實驗材料

本實驗所使用的雙二倍體鯽(C. auratus)和雙三倍體銀鯽(C. gibelio)來自于國家水生生物種質資源庫(National Aquatic Biological Resource Center,NABRC), 中國科學院水生生物研究所。

1.2 中期染色體制備

雙三倍體銀鯽和雙二倍體鯽有絲分裂中期染色體的制備采用植物血球凝集素體內(nèi)誘導腎細胞制片技術, 具體參考朱華平等[22]描述的方法進行。

1.3 染色體熒光原位雜交(Fluorescence in situ hybridization, FISH)

FISH實驗參照丁苗等[20]描述的方法進行, 但方法有略微改動。FISH探針為超數(shù)微小染色體探針[20]和酵母人工染色體(Bacterial artificial chromosome,BAC)質粒探針。BAC質粒的接種和提取參考彭金霞等[23]描述的方法進行。其中, BAC質粒探針的制備按照DIG-Nick Translation Mix試劑盒(Roche公司)的說明書進行。染色體玻片于65℃放置3h后,移至70℃變性液(35 mL去離子甲酰胺, 10 mL滅菌水和5 mL 2×SSC)中變性3min。將玻片放入預冷的70%、90%和100%的乙醇中依次脫水5min, 自然晾干玻片。與此同時, 將探針混合液(50%去離子甲酰胺, 20×SSC, 10% SDS, 0.25 mg/mL 鮭精子DNA,50%硫酸葡聚糖, 100 ng BAC質粒探針和5 ng超數(shù)微小染色體探針), 于100℃的水中變性10min, 而后將探針混合液立即置于冰上; 將探針變性混合液孵育染色體玻片, 37℃避光雜交24h。

在雜交完畢后, 玻片放置于20%去離子甲酰胺(35 mL滅菌水, 10 mL去離子甲酰胺和5 mL 20×SSC)中洗10min。隨后于2×SSC(含1% Tween)中洗5min, 之后放于1×PBST中洗3次, 每次5min; 1×PBS中洗3次, 每次5min。接著于每張玻片中加200 μL羊抗Dig稀釋液(羊抗Dig﹕1×PBS =1﹕100), 37℃濕盒中避光孵育1h。然后將玻片于1×PBS中洗3次, 每次5min, 自然晾干后, 在每張玻片用150 μL(2 μg/mL)的DAPI染色, 蓋上蓋玻片, 37℃濕盒避光染色1h。最后將玻片放于1×PBS中洗3次, 每次5min, 加30 μL的抗淬滅劑封片, 借助共聚焦顯微鏡鏡檢。

1.4 DNA提取和PCR

剪取適量大小的尾鰭, 用干凈濾紙擦干鰭條表面黏液后, 將樣品放于1.5 mL的EP管中, 然后加入500 μL 核酸裂解液和10 μL(10 mg/mL)蛋白酶K, 放于55℃水浴鍋中1h; 待樣品裂解充分后, 14000 r/min離心5min, 將上清液轉入新的1.5 mL EP管中, 加入200 μL 蛋白沉淀劑, 上下顛倒混勻, 置于冰上10min;待樣品充分沉淀后, 14000 r/min離心5min, 吸取上清液于新的1.5 mL EP管中, 加入500 μL預冷的異丙醇, 上下顛倒混勻(可見白色絮狀物); 14000 r/min離心5min, 棄上清, 加入500 μL 80%乙醇混勻; 14000 r/min離心5min, 棄上清, 室溫放置待白色沉淀變?yōu)橥该魑? 加入100 μL滅菌水, 混勻以便充分溶解DNA。采用NanoDrop測濃度, 將樣品DNA濃度調(diào)為100 ng/μL,調(diào)整好濃度的樣品存于4℃待用。

PCR模板為雙三倍銀鯽和雙二倍體鯽基因組DNA, 反應體系為20 μL, 包括10 μL的Taqreaction Mix, 0.5 μL的上下游引物, 1 μL的模板和8 μL的滅菌水。擴增Cg-Ca-CL1、bmb-CgB、ccna2-CgB和nusap1-CgB,這些序列的引物在表1中, PCR反應程序為: 95℃ 3min; 95℃ 15s, 60℃ 20s, 72℃ 20s,35個循環(huán); 72℃ 5min, 12℃ 5min。待PCR程序結束后, PCR產(chǎn)物用1.5%的瓊脂糖凝膠檢測, 檢測時采用的Marker為TaKaRa DL2000。

表1 引物序列Tab. 1 Primer sequences

1.5 共線性分析

以已發(fā)表的銀鯽基因組及基因組BAC序列為參考數(shù)據(jù)庫,Cg-Ca-CL1為查詢序列, 借助 BLAST軟件[24], 使用默認參數(shù)比對, 獲取超數(shù)微小染色體基因組序列, 再將鑒定到的超數(shù)微小染色體序列分別和已發(fā)表的雙三倍體銀鯽常規(guī)A染色體序列(GenBank ID: PRJNA546443)及雙二倍體鯽基因組(GenBank ID: PRJNA546444)通過Minimap2[25]進行

基因組共線性分析, 結果借助Circos[26]展示。

1.6 重復序列注釋和基因通路分析

使用Repeatmask[27]軟件進行重復序列分析。對超數(shù)染色體序列上的注釋基因進行KEGG通路分析和GO通路分析(https://david.ncifcrf.gov/)。

1.7 熒光定量PCR(Real-time fluorescence quantitative PCR, qPCR)

qPCR采用SYBR Green染料法根據(jù)操作說明進行基因表達量的檢測。qPCR反應體系為20 μL, 包括10 μL的 SYBR Mix, 0.5 μL的上下游引物(表1),1 μL的模板和8 μL的滅菌水。反應程序為: 95℃預變性, 95℃ 15s, 60℃ 20s, 72℃ 20s, 40個循環(huán)。每個反應進行3次重復, 采用2-ΔΔCt計算獲得基因的相對表達量。

2 結果

2.1 雙三倍體銀鯽超數(shù)微小染色體序列鑒定

在前期研究中篩選到雙三倍體銀鯽中含量最高的重復序列Cg-Ca-CL1, 利用該重復序列為探針,通過染色體FISH能夠特異性標記銀鯽的超數(shù)微小染色體[20]。雙三倍體銀鯽由祖先雙二倍體鯽經(jīng)同源三倍化而形成[14], 利用Cg-Ca-CL1探針在雙二倍體鯽中進行染色體FISH分析, 發(fā)現(xiàn)雙二倍體鯽不含有超數(shù)微小染色體, 且常染色體和性染色體上也沒有明顯熒光信號(圖1A)。由于Cg-Ca-CL1在雙二倍體鯽中的含量低, 通過FISH在雙二倍體鯽中檢測不到Cg-Ca-CL1信號, 但是利用PCR擴增及其產(chǎn)物序列測定, 能夠在雙二倍體鯽中檢測到Cg-Ca-CL1(圖1B)。在雙二倍體鯽中,Cg-Ca-CL1擴增條帶單一, 說明該序列不是以串聯(lián)重復形式存在; 而在Cg-Ca-CL1在雙三倍體銀鯽中PCR擴增出不同長度條帶, 說明Cg-Ca-CL1發(fā)生了擴增且以串聯(lián)重復的形式存在(圖1B)。這些結果揭示同源三倍化之后,Cg-Ca-CL1在雙三倍體銀鯽的超數(shù)微小染色體上發(fā)生了擴增和富集。因此, 我們擬通過該重復序列并結合染色體熒光原位雜交來篩選雙三倍體銀鯽超數(shù)微小染色體序列。

圖1 重復序列Cg-Ca-CL1分析Fig. 1 Analysis of repetitive sequence Cg-Ca-CL1

對雙三倍體銀鯽雌性基因組(GenBank: PRJNA 546443)進行了重復序列Cg-Ca-CL1搜尋, 總共發(fā)現(xiàn)18個未組裝的scaffold和12個組裝好的染色體含有該重復序列, 包含Cg-Ca-CL1的拷貝數(shù)從5到3157不等(圖1C)。這18個scaffold由于未能組裝到染色體, 可能是潛在超數(shù)微小染色體序列。隨后, 在這尾測序個體構建的酵母人工染色體文庫中, 篩選到500個 BAC質粒含有Cg-Ca-CL1。拿這些BAC質粒序列與上述的18個scaffold進行序列比對, 最終鑒定出85個BAC質粒序列能比對上其中的15個scaffold(序列相似度大于95%)。在15個scaffold中, 分別選擇一個能匹配的BAC質粒為探針(表2), 通過FISH進行染色體定位分析。FISH結果顯示, 這15個BAC質粒的信號均位于超數(shù)微小染色體上, 而在常規(guī)染色體上沒有明顯信號, 因此, 可以推測這些BAC質粒所對應的scaffold序列是超數(shù)微小染色體序列(圖2)。

圖2 BAC質粒的染色體定位Fig. 2 Chromosomal localization of BAC plasmids

表2 超數(shù)微小染色體序列概況Tab. 2 Summary of sequences from supernumerary microchromosomes

2.2 超數(shù)微小染色體含有常規(guī)染色體同源序列

由于雙三倍體銀鯽由祖先雙二倍體鯽經(jīng)同源三倍化而形成[14], 同時選取了雙三倍體銀鯽及雙二倍體鯽為對照, 通過序列比對來揭示雙三倍體銀鯽超數(shù)微小染色體的序列特征。雙三倍體銀鯽中未組裝的序列不能確定是常規(guī)染色體還是超數(shù)微小染色體序列, 所以我們舍棄了這部分未組裝的序列,只用組裝到染色體的50條常規(guī)染色體序列來做超數(shù)微小染色體同源性分析。雙二倍體鯽中由于不含有超數(shù)微小染色體(圖1A), 其未組裝的序列都是來源于常規(guī)染色體, 所以雙二倍體鯽所有序列都用來做超數(shù)微小染色體同源性分析。我們發(fā)現(xiàn)銀鯽5.48%的超數(shù)微小染色體序列與銀鯽組裝的常規(guī)染色體具有同源性, 且同源序列在每個常規(guī)染色體上都有分布, 其中含有超數(shù)微小染色體同源序列最多的常規(guī)染色體依次是chr1B、chr4A和chr16A(圖3A)。同時, 與雙二倍體鯽基因組進行比對時發(fā)現(xiàn), 22.28%的超數(shù)微小染色體序列在雙二倍體鯽基因組中能找到同源序列, 其中含有超數(shù)微小染色體同源序列最多的是未組裝的序列及chr25B、chr21B和chr22B號染色體(圖3B)。由于雙三倍體銀鯽中部分屬于常規(guī)染色體但是又未能組裝到染色體的序列沒有用來做同源分析, 所以導致銀鯽常規(guī)染色體上同源序列比例(5.48%)遠低于鯽常規(guī)染色體(22.28%)。從上述的結果來看, 銀鯽超數(shù)微小染色體含有常規(guī)染色體同源序列, 暗示著銀鯽的超數(shù)微小染色體可能起源于常規(guī)染色體。

圖3 超數(shù)微小染色體同源序列特征Fig. 3 Homologous sequence characterization of supernumerary microchromosomes

2.3 超數(shù)微小染色體重復序列擴增

雙三倍體銀鯽超數(shù)微小染色體、雙三倍體銀鯽組裝的常規(guī)染色體及雙二倍體鯽常規(guī)染色體的GC含量分別為38.06%、37.61%和37.63%, 三者沒有顯著差異(圖4A)。隨后通過重復序列鑒定發(fā)現(xiàn)超數(shù)微小染色體的重復序列占所有序列的比重為78.29%, 比雙三倍體銀鯽組裝的常規(guī)染色體(47.93%)和雙二倍體鯽常規(guī)染色體(47.56%)的重復序列占比明顯偏高(圖4A)。相比雙三倍體銀鯽組裝的常規(guī)染色體(衛(wèi)星序列: 8.20%; TE: 69.45%)及雙二倍體鯽常規(guī)染色體(衛(wèi)星序列: 7.59%; TE: 72.53%), 超數(shù)微小染色體重復序列中衛(wèi)星序列(21.61%)的占比有所增加而TE的占比有所下降(50.67%; 圖4B)。此外, 超數(shù)微小染色體衛(wèi)星序列中Y染色體類型占比增加到52.05%(圖4C), TE中DNA類型占比減少到42.20%(圖4D)。 這些結果說明重復序列尤其是衛(wèi)星序列在雙三倍體銀鯽超數(shù)微小染色體上發(fā)生了大量擴增。

圖4 重復序列分析Fig. 4 Analysis of repetitive elements

2.4 超數(shù)微小染色體注釋基因通路分析

超數(shù)微小染色體序列上總共注釋到了318個完整基因, 其中包含141個已知基因和177個未知基因(圖5A)。本研究使用已知基因進行了KEGG(Kyoto Encyclopedia of Genes and Genomes)和GO(Gene Ontology)通路分析。在KEGG分析中, 213個基因能夠匹配到6個主要類別中的31個通路。包含最多基因數(shù)的通路是癌癥: 概述(16個基因); 膜運輸(14個基因); 癌癥: 特定類型以及全局和總覽圖(13個基因) (圖5B)。此外, 在GO分析中, 分別有401個、266個和235個基因能夠匹配到生物學過程、細胞學組分及分子生物學功能這3個主要類別中的47個通路。其中, 含有注釋基因組最多的通路依次是結合通路(83個基因), 細胞通路(72個基因)和細胞部分通路(53個基因; 圖5C)。

圖5 超數(shù)微小染色體注釋基因通路分析Fig. 5 Pathway analysis of annotated genes from supernumerary microchromosomes

2.5 超數(shù)微小染色體上減數(shù)分裂相關基因的鑒定及表達分析

前期的研究通過多個品系的比較基因組解析發(fā)現(xiàn)與雙二倍體鯽相比, 雙三倍體銀鯽中有13個基因家族發(fā)生了擴張, 其中有8個是減數(shù)分裂相關基因家族, 這些基因的擴張可能與銀鯽進行雌核生殖特殊的卵子發(fā)生方式有關[14]。本研究發(fā)現(xiàn)在這8個擴增的減數(shù)分裂相關基因家族中, 有6個在超數(shù)微小染色體序列中都能檢測到擴增拷貝。紡錘體組成相關基因incenp在HiC_sacffold_51和HiC_sacffold_53分別有1個拷貝; 細胞周期相關基因fbxo5在HiC_sacffold_57上有2個拷貝; 核膜組裝相關基因lap2在HiC_sacffold_56和HiC_sacffold_204上分別有1個拷貝;bmb在HiC_sacffold_64上有1個拷貝;ccna2在HiC_sacffold_53上有2個拷貝;nusap1在HiC_sacffold_166上有1個拷貝。其中,incenp、fbxo5和lap2這3個基因除了在雙三倍體銀鯽超數(shù)微小染色體序列上存在拷貝, 還分別在雙三倍體銀鯽的染色體chr3B、chr18B和chr9A上存在1個擴增拷貝; 而bmb、ccna2和nusap1這三個基因則只在雙三倍體銀鯽超數(shù)染色體序列上存在擴增拷貝。

本研究針對這6個基因在雙三倍體銀鯽中擴增的拷貝[incenp-Cg(B+chr3B)、fbxo5-Cg(B+chr18B)、lap2-Cg(B+chr9A)、bmb-CgB、ccna2-CgB和nusap1-CgB]設計了特異的PCR引物(表1), 分別在雙三倍體銀鯽出膜后30d、70d、160d、200d和1y(性成熟時期)的雌性性腺中進行了qPCR檢測。我們發(fā)現(xiàn)這些基因拷貝在銀鯽卵子發(fā)生過程中幾乎均有轉錄表達(圖6), 暗示著超數(shù)微小染色體上擴增的這些基因具有活性, 可能在銀鯽特殊的卵子發(fā)生過程中發(fā)揮作用。

圖6 雙三倍體銀鯽特異擴增減數(shù)分裂相關基因的轉錄活性檢測Fig. 6 Transcriptional activity detection of meiosis-related genes that are specifically expanded in amphitriploid C. gibelio

2.6 超數(shù)微小染色體減數(shù)分裂相關基因SCAR標記開發(fā)

此外,bmb-CgB、ccna2-CgB和nusap1-CgB這些超數(shù)微小染色體特異拷貝引物還可以用來開發(fā)超數(shù)微小染色體特異的SCAR(sequence characterized amplified regions)標記。我們在不含超數(shù)微小染色體的雙二倍體鯽和含有超數(shù)微小染色體的雙三倍體銀鯽中進行了PCR檢測, 發(fā)現(xiàn)這3對引物在雙三倍體銀鯽雌性和雄性個體中均能擴增出特異條帶, 而在雙二倍體鯽雌性和雄性個體中均擴增不出條帶(圖7)。因此, 這些SCAR標記還可以作為區(qū)分雙二倍體鯽和雙三倍體銀鯽的分子標記。

3 討論

性別在自然界中普遍存在, 廣義的有性生殖是指發(fā)生減數(shù)分裂以及隨后核子融合[28]。越來越多的證據(jù)表明有性生殖是真核生物的共源性狀[29,30]。然而在自然界中的確存在少數(shù)的生物類群進行單性生殖, 其中脊椎動物約有100種[31]。我們通常所說的單性生殖包括孤雌生殖、雌核生殖和雜種生殖[32]。孤雌生殖是雌性親本產(chǎn)生未減數(shù)的卵子, 這些卵子獨自發(fā)育成與母本遺傳背景一致的個體, 不需要雄性[33]。雌核生殖同樣也是雌性親本產(chǎn)生不減數(shù)的卵子, 但是卵子須要外源精子刺激胚胎發(fā)育,形成與母本遺傳背景一致的個體[17,18]。行雜種生殖的生物, 雌性產(chǎn)生減數(shù)的卵子, 但是卵子只含有母本的單倍體基因組, 父本的基因組在形成卵子的時候被排除[34]。單性生殖沒有減數(shù)分裂同源重組,導致有害突變積累及阻礙遺傳多樣性產(chǎn)生, 通常被認為是演化的死胡同。值得注意的是, 有些單性生殖生物類群存在的時間已經(jīng)遠超出其預測的滅絕時間, 且展示出很高的遺傳多樣性和很強的環(huán)境適應性[13,35,36]。因此, 單性動物對于回答生態(tài)學和演化生物學領域的基礎科學問題具有重要意義。

超數(shù)染色體由于其非必需的特質, 最初被認為是沒有或者很少有生物學功能的。隨著研究的拓展和深入, 逐漸發(fā)現(xiàn)超數(shù)染色體對于所在的生物類群能夠發(fā)揮重要的生物學功能和演化適應性[37]: 比如在真菌中, 超數(shù)染色體發(fā)現(xiàn)與致病性有關[38]; 在植物中, 超數(shù)染色體能夠增強植物的耐受性[39]; 在許多動物中, 超數(shù)染色體具有性別決定功能[40,41],近期在洞穴魚的超數(shù)染色體上還鑒定了性別決定基因[5]。我們前期的研究也發(fā)現(xiàn), 在銀鯽雄性中, 雄性特異的超數(shù)微小染色體上含有在性別決定關鍵時期雄性特異或者雄性偏向表達的基因[20], 且與銀鯽的雄性發(fā)生緊密相關[21]。在本文的研究中, 發(fā)現(xiàn)銀鯽中擴張的8個減數(shù)分裂相關基因, 有6個在超數(shù)微小染色體上存在擴增拷貝。其中incenp、fbxo5和lap2這3個基因在常規(guī)染色體和超數(shù)微小染色體上都存在擴增拷貝, 且熒光定量檢測引物沒有辦法區(qū)分常規(guī)染色體和超數(shù)微小染色體上的擴增拷貝,因此這3個基因在卵巢發(fā)育不同時期的轉錄表達是常規(guī)染色體和超數(shù)微小染色體擴增拷貝共同表達的結果。另外3個基因(bmb、ccna2和nusap1)則只在超數(shù)微小染色體上存在擴增拷貝, 且這3個基因在卵巢發(fā)育的不同階段都具有轉錄表達。這些結果暗示著超數(shù)微小染色體在雙三倍體銀鯽特殊的單性雌核生殖過程中可能發(fā)揮作用。

此外, 雙二倍體鯽和雙三倍體銀鯽共同組成了鯽復合種, 雙三倍體銀鯽是由祖先雙二倍體鯽經(jīng)同源三倍化而形成, 他們廣泛的分布于亞歐大陸及附近島嶼的淡水水系, 有著部分重疊的生態(tài)分布[12,13]。雙二倍體鯽和雙三倍體銀鯽是我國重要的大宗淡水養(yǎng)殖魚類, 在外觀和形態(tài)上沒有明顯差異, 通常我們是通過流式細胞儀對其血細胞的DNA含量進行測定來對二者進行區(qū)分[42]。通過流式細胞儀進行區(qū)分的流程比較繁瑣, 樣品處理不當還容易造成結果不準確。本研究中開發(fā)的超數(shù)微小染色體特異SCAR標記, 能夠很好地區(qū)分雙二倍體鯽和雙三倍體銀鯽, 為兩者的快速鑒定提供了一條快捷和準確的途徑。

總之, 本研究鑒定了雙三倍體銀鯽超數(shù)微小染色體序列, 揭示了銀鯽超數(shù)微小染色體含有常規(guī)染色體同源序列且積累了大量重復序列, 發(fā)現(xiàn)超數(shù)微小染色體上含有活性基因, 證明減數(shù)分裂相關基因通過超數(shù)微小染色體發(fā)生擴張, 開發(fā)了超數(shù)微小染色體特異的SCAR標記。本研究不僅解析了銀鯽超數(shù)微小染色體的基因組特征, 也為超數(shù)染色體在單性生物演化中的作用提供了創(chuàng)新見解。