不同生理階段邊雞卵巢轉(zhuǎn)錄譜的構(gòu)建及卵泡發(fā)育相關(guān)基因的分析

張俊珍，李彩娥，劉 博，李建慧，張 蒙，李亞莉，常強(qiáng)強(qiáng)

(山西農(nóng)業(yè)大學(xué)動物科學(xué)學(xué)院，太谷 030801)

家禽的就巢性影響著其產(chǎn)蛋性能，就巢性越強(qiáng)的品種產(chǎn)蛋性能越低。卵泡的發(fā)育與環(huán)境、營養(yǎng)和激素調(diào)控有關(guān)系，決定著家禽就巢行為的發(fā)生和維持，對產(chǎn)蛋性能具有重要影響。卵泡發(fā)育取決于細(xì)胞的增殖、分化和卵黃質(zhì)大量沉積，進(jìn)而促進(jìn)了卵泡質(zhì)量的快速增長[1]。在家禽生產(chǎn)過程中，一些與卵泡發(fā)育相關(guān)的基因參與調(diào)控家禽的就巢性。隨著家禽選育程度的不斷提高，多基因調(diào)控研究在家禽遺傳改良中得到廣泛的應(yīng)用。

邊雞是肉蛋兼用的中國地方品種，主要分布在內(nèi)蒙古自治區(qū)與山西省北部相毗連的長城內(nèi)外一帶，具有蛋重大、肉質(zhì)好、適應(yīng)性強(qiáng)、耐粗飼和抗寒等優(yōu)點。由于邊雞具有很強(qiáng)的就巢性，產(chǎn)蛋量很少，限制了該品種資源的開發(fā)和利用。因此，提高邊雞產(chǎn)蛋性能已受到相關(guān)家禽育種研究者的廣泛關(guān)注。家禽育種研究者通常以卵巢作為研究家禽就巢性的試驗材料，從分子水平上研究與就巢性相關(guān)基因的差異表達(dá)。江紅霞等[2]對克氏原螯蝦的卵巢和肝胰腺進(jìn)行轉(zhuǎn)錄組測序，共獲得53 006個單基因(unigene)，卵巢與肝、胰對比分析，發(fā)現(xiàn)卵巢中有20 382個差異性表達(dá)基因(differentially expressed genes, DEGs)，通過GO功能分類分析發(fā)現(xiàn)，部分DEGs被注釋到繁殖(reproduction)、繁殖過程(reproduction process)、免疫系統(tǒng)過程(immune system process)和生長(growth)GO條目。孟金柱等[3]對貴州白山羊的優(yōu)勢卵泡和從屬卵泡進(jìn)行轉(zhuǎn)錄組測序，篩選出233個上調(diào)表達(dá)基因，并對這些基因進(jìn)行GO、KEGG分析，最終篩選出4個與卵泡發(fā)育有關(guān)的基因。郝慶玲等[4]對山西省文水縣海福特母牛的不同生理階段發(fā)育卵泡進(jìn)行轉(zhuǎn)錄組分析，共獲得608個差異表達(dá)基因，其中有9個參與牛卵泡發(fā)育的信號通路調(diào)控，并篩選出11個直接參與信號通路調(diào)節(jié)及細(xì)胞增殖分化的基因。朱志明等[5]對山麻鴨開產(chǎn)期和產(chǎn)蛋高峰期卵巢組織轉(zhuǎn)錄組測序，發(fā)現(xiàn)存在1 929 個 差異表達(dá)基因，其中上調(diào)基因 989 個，下調(diào)基因 940 個，并將這些基因在Nr、GO、KEGG數(shù)據(jù)庫進(jìn)行分析比較，獲得了差異表達(dá)基因的功能、分類和代謝通路。Shen等[6]使用轉(zhuǎn)錄組測序來分析靖黃雞三種類型的卵泡細(xì)胞中circRNA和mRNA，利用生物信息學(xué)分析確定了14 502個circRNA，其中5 622個廣泛分布在卵泡細(xì)胞發(fā)育的所有階段，差異表達(dá)分析表明，在與生殖相關(guān)的途徑中富集了不同表達(dá)的circRNA和mRNA，包括TGF-β信號通路、卵母細(xì)胞減數(shù)分裂和血管平滑肌收縮，為卵泡發(fā)育提供了理論依據(jù)。董新龍等[7]利用轉(zhuǎn)錄組對京海黃雞卵巢組織進(jìn)行分析，發(fā)現(xiàn)了4 431個新轉(zhuǎn)錄本，這些新轉(zhuǎn)錄本比對到 GO、Nr、KEGG數(shù)據(jù)庫的數(shù)量分別為 917、1 795、422。與這些數(shù)據(jù)庫的成功比對，為進(jìn)一步證實新基因的存在提供了基礎(chǔ)；此外，經(jīng)過與KEGG數(shù)據(jù)庫比對，發(fā)現(xiàn)絲裂原活化蛋白激酶信號通路(MAPK signaling pathway)排在富集通路的首位。雖然轉(zhuǎn)錄組測序已經(jīng)在畜禽不同生理階段卵巢組織差異表達(dá)基因的研究中廣泛應(yīng)用，但是目前尚無對于邊雞卵巢組織轉(zhuǎn)錄組差異表達(dá)基因的研究。本研究為了挖掘邊雞卵巢卵泡發(fā)育過程中起調(diào)節(jié)作用的遺傳因子，揭示卵巢發(fā)育期、開產(chǎn)期、產(chǎn)蛋高峰期和就巢期卵巢組織的差異表達(dá)基因，篩選與卵泡發(fā)育相關(guān)的候選基因，利用轉(zhuǎn)錄組測序技術(shù)構(gòu)建不同日齡的邊雞卵巢的轉(zhuǎn)錄組文庫，對邊雞在不同生理狀態(tài)下卵巢組織做轉(zhuǎn)錄組差異分析、差異表達(dá)基因GO功能富集分析和KEGG分析，將差異表達(dá)基因進(jìn)行功能注釋，旨在研究與邊雞卵巢發(fā)育相關(guān)的基因，以揭示邊雞就巢性強(qiáng)的遺傳機(jī)制。

1 材料與方法

1.1 試驗動物和卵巢樣本采集

試驗用邊雞在標(biāo)準(zhǔn)化飼養(yǎng)條件下進(jìn)行飼養(yǎng)管理，每次取材前進(jìn)行連續(xù)7 d的產(chǎn)蛋觀察：產(chǎn)下第一枚蛋起，表示性成熟；產(chǎn)蛋率達(dá)到70%以上，代表產(chǎn)蛋高峰；320日齡后產(chǎn)蛋率下降為30%以下，觀察記錄一些健康的停產(chǎn)并表現(xiàn)出就巢行為的個體。隨機(jī)選取70日齡(卵巢發(fā)育期)、165日齡(開產(chǎn)期)、220日齡(產(chǎn)蛋高峰期)和330日齡(就巢期)健康、體重相近且符合取樣條件的邊雞個體，每個日齡5只，頸靜脈放血法處死，采集卵巢組織，迅速投入液氮中，轉(zhuǎn)運后置于-80 ℃超低溫冰箱中保存?zhèn)溆谩?/p>

1.2 cDNA文庫構(gòu)建和測序

采用TRIzol試劑(美國Invitrogen)提取不同日齡邊雞的卵巢組織總RNA，用Agilent 2100生物分析儀進(jìn)行總RNA品質(zhì)檢測。用帶有poly-T oligo附著的磁珠富集含有poly A的mRNA，使用逆轉(zhuǎn)錄酶和隨機(jī)引物將mRNA反轉(zhuǎn)錄合成第一條cDNA鏈，使用DNA polymerase I和RNase H合成第二條cDNA鏈。通過PCR擴(kuò)增方法純化和富集cDNA，構(gòu)建邊雞卵巢cDNA文庫。在BGISEQ-500平臺進(jìn)行測序(華大基因公司，深圳，中國)。

1.3 生物信息學(xué)分析

過濾低質(zhì)量片段，去除帶接頭污染和未知堿基含量比例過高的片段獲得clean reads。Clean reads與參考基因組比對預(yù)測新基因和轉(zhuǎn)錄本，檢測基因剪接事件，識別卵巢樣本間的差異表達(dá)基因(DEGs)，進(jìn)行聚類分析和功能注釋。

1.4 新轉(zhuǎn)錄本預(yù)測

利用HISAT軟件將clean reads定位到雞參考基因組?；蚪M定位后，用StringTie[8]重建轉(zhuǎn)錄本，基于基因組注釋信息，用Cuffcompare[9]識別新的轉(zhuǎn)錄本，并使用CPC[10]預(yù)測新轉(zhuǎn)錄本的編碼能力。

1.5 基因表達(dá)、DEGs檢測、GO功能富集和KEGG分析

利用Bowtie2[11]將新的編碼轉(zhuǎn)錄本與轉(zhuǎn)錄本合并，將clean reads映射成完整的參考基因組后，利用RSEM[12]計算基因表達(dá)水平。使用Cor計算卵巢樣本間的Pearson相關(guān)性，使用R軟件里的 hclust 函數(shù)進(jìn)行層次聚類分析，圖形的繪制采用了R軟件中的ggplot2包。

采用PossionDis[13]算法檢測不同日齡邊雞卵巢間的差異表達(dá)基因。根據(jù)GO注釋結(jié)果，對DEGs進(jìn)行功能分類和生物通路分類，同時使用R軟件中的phyper函數(shù)進(jìn)行GO的功能和KEGG富集，用超幾何分布函數(shù)計算出P值，計算每個P值的錯誤發(fā)現(xiàn)率(FDR)，當(dāng)FDR≤0.01時則被認(rèn)為顯著富集。

1.6 實時熒光定量PCR(qPCR)驗證

為了驗證轉(zhuǎn)錄組測序的結(jié)果，從已知參與卵泡發(fā)育的差異表達(dá)基因中隨機(jī)選擇11個基因，通過qPCR進(jìn)行驗證。將上述從測序同樣的個體樣本中提取的不同日齡邊雞卵巢的總RNA用MMLV逆轉(zhuǎn)錄酶(TaKaRa，大連，中國)將其逆轉(zhuǎn)錄為cDNA。通過SYBR方法進(jìn)行qPCR驗證(引物見表1)。擴(kuò)增條件為：預(yù)變性95 ℃ 30 s；(變性95 ℃ 5 s、退火55 ℃ 30 s、延伸72 ℃ 10 s)×40循環(huán)；熔解過程(95 ℃ 1 min、55 ℃ 30 s、55 ℃ 30 s)。內(nèi)參基因β-actin與目的基因的樣本要在同樣的條件下進(jìn)行qPCR，通過2-△△CT法計算目的基因的相對表達(dá)量，結(jié)果以“平均值±標(biāo)準(zhǔn)差(Means±SD)”表示，用Graph-Pad Prism 8.0軟件分析數(shù)據(jù)，采用one-way ANOVA進(jìn)行單因素方差分析。每個樣本做4次技術(shù)重復(fù)。

表1 引物序列

2 結(jié) 果

2.1 轉(zhuǎn)錄組測序原始數(shù)據(jù)及轉(zhuǎn)錄本分析

通過在BGISEQ-500平臺上對不同日齡邊雞卵巢樣本進(jìn)行測序，構(gòu)建基因表達(dá)譜，每個樣品約產(chǎn)生11.04 Gb的序列。共鑒定出18 891個基因，其中已知基因18 063個，預(yù)測的新基因828個。此外，共檢測出21 299個新的轉(zhuǎn)錄本，其中14 376個屬于已知基因且具有蛋白編碼潛力的新轉(zhuǎn)錄本，925個屬于新基因的具有蛋白編碼潛力的新轉(zhuǎn)錄本，其余5 998個為長鏈非編碼RNA。

2.2 不同生理階段邊雞卵巢基因表達(dá)及相關(guān)性分析

通過RSEM分析每組樣本的基因表達(dá)水平，計算每兩個樣本之間的Pearson相關(guān)性并對所有樣本進(jìn)行層次聚類分析。結(jié)果表明，在165與220 d邊雞卵巢基因表達(dá)聚類最近，二者間的Pearson值為0.983，說明有顯著相關(guān)性；而70和330 d與165和220 d邊雞卵巢基因表達(dá)均未聚類，且相互間相關(guān)性相對較低(圖1)。

樣品間彼此靠的越近代表樣品的表達(dá)譜越相似The closer the samples are to each other, the more similar the expression profiles of the samples are圖1 不同日齡邊雞卵巢表達(dá)基因的層次聚類和Pearson相關(guān)性Fig.1 Hierarchical clustering and Pearson correlation of genes expressed among the ovaries of Bian chicken at different ages

2.3 不同日齡邊雞卵巢基因表達(dá)譜中的差異表達(dá)基因

通過PossionDis計算差異表達(dá)基因水平，70、165、220和330 d 邊雞4個發(fā)育階段的卵巢轉(zhuǎn)錄組數(shù)據(jù)兩兩對比，形成6組差異表達(dá)的mRNA(圖2)。通過比較發(fā)現(xiàn)，70 d邊雞卵巢轉(zhuǎn)錄組與其他階段差異最大，獲得的差異表達(dá)基因也相對較多，而且下調(diào)基因的數(shù)量大于上調(diào)基因；而165與330 d 相比、220與 330 d相比，差異表達(dá)基因的上調(diào)數(shù)量大于下調(diào)基因的數(shù)量；165與220 d相比，差異表達(dá)基因的總量最少。

以差異倍數(shù)>1.5，F(xiàn)DR<0.01為標(biāo)準(zhǔn)，使用DESeq2R包對轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行篩選，找出了在不同生理階段與發(fā)育卵泡數(shù)變化相關(guān)的差異表達(dá)基因40個(表2)。結(jié)果表明，有18個基因參與卵泡發(fā)育的整個過程，有22個基因僅在220 d前表達(dá)。其中35個基因的表達(dá)量與發(fā)育卵泡數(shù)成正相關(guān)，F(xiàn)DX1L在70、220、165和330 d邊雞卵巢中表達(dá)量依次遞減，NDRG4在70、165、220和330 d邊雞卵巢中表達(dá)量隨日齡增加呈遞減趨勢，MC3R和PRLHR在70和330 d邊雞卵巢中表達(dá)量高于165和220 d的表達(dá)量，VIP僅在70、165和220 d邊雞卵巢中表達(dá)。

圖2 不同日齡邊雞卵巢中的差異表達(dá)基因Fig.2 Differentially expressed genes in the ovary of Bian chickens at different ages

表2 不同日齡邊雞卵巢中可能與卵泡發(fā)育相關(guān)的差異表達(dá)基因及其表達(dá)量(RPKM值)

(轉(zhuǎn)下頁 Carried forward)

(轉(zhuǎn)下頁 Carried forward)

2.4 不同日齡邊雞卵巢差異表達(dá)基因GO和KEGG富集分析

GO分為三大類，即生物過程、細(xì)胞組分、分子功能。將4個混合樣本相互組合形成6個對比組，對篩選獲得的差異表達(dá)基因進(jìn)行GO功能富集分析并繪制差異表達(dá)基因GO富集柱狀圖(圖3)。結(jié)果表明：70與165 d相比，共有1 748個差異基因，其中有615個 得到功能注釋；70與220 d相比，有1 704個差異基因，其中596個得到功能注釋；70 與330 d相比，有2 013個差異基因，其中有622個得到功能注釋； 165與220 d相比，有656個差異基因，其中169個 得到功能注釋；165與330 d相比，有1 541個差異基因，其中有480個得到功能注釋；220與 330 d相比，有1 593個差異基因，其中有474個得到功能注釋。差異基因主要集中在“生物過程”，其次是“細(xì)胞組分”，最后是“分子功能”。對比6個組合的差異基因功能富集GO term，其差異基因在“生物過程”分類中，細(xì)胞過程(cellular process)所占比例最多，其次是生物調(diào)控(biological regulation)和生物過程調(diào)控(regulation of biological process)，在“細(xì)胞組分”分類中，細(xì)胞(cell)占的比例最大，其次是細(xì)胞部分(cell part)和細(xì)胞膜等；在“分子功能”分類中，結(jié)合(binding)占的比例最多，其次是催化活性(catalytic activity)。此外，可以看到“生物過程”分類中與發(fā)育繁殖相關(guān)的有發(fā)育過程 (developmental process)、繁殖(reproduction)、繁殖過程(reproductive process)。在70與165 d相比，共涉及272個與發(fā)育繁殖相關(guān)的差異表達(dá)基因；70與220 d相比，共涉及269個與發(fā)育繁殖相關(guān)的差異表達(dá)基因；70與330 d相比，共涉及262個與發(fā)育繁殖相關(guān)的差異表達(dá)基因；165與 220 d相比，共涉及57個與發(fā)育繁殖相關(guān)的差異表達(dá)基因；165與330 d相比，共涉及194個與發(fā)育繁殖相關(guān)的差異表達(dá)基因；220 d 與330 d相比，共涉及188個與發(fā)育繁殖相關(guān)的差異表達(dá)基因。

通過對卵泡發(fā)育的不同生理階段的邊雞卵巢差異表達(dá)基因KEGG通路比較均可發(fā)現(xiàn)信號轉(zhuǎn)導(dǎo)、免疫系統(tǒng)以及信號分子與相互作用路徑上富集最多(表3)。

2.5 差異表達(dá)基因mRNA的實時定量PCR(qPCR)驗證

為進(jìn)一步驗證轉(zhuǎn)錄組測序的可靠性，本研究通過qPCR隨機(jī)驗證了11個差異基因的表達(dá)。結(jié)果表明，所選基因的mRNA表達(dá)趨勢與轉(zhuǎn)錄組測序數(shù)據(jù)一致(圖4)，可以證明轉(zhuǎn)錄組測序的結(jié)果可靠。

3 討 論

家禽的就巢行為與卵巢發(fā)育有直接聯(lián)系，就巢行為的發(fā)生會引起產(chǎn)蛋性能的下降。轉(zhuǎn)錄組測序以它檢測范圍廣、測序通量高、檢測成本低、分辨率好和靈敏度高等特點，已經(jīng)被廣泛應(yīng)用[14]，在金定鴨的卵巢[15]、京海黃雞卵巢[16]中都有應(yīng)用。本研究對邊雞卵巢組織的轉(zhuǎn)錄組測序，共鑒定出18 891個基因，其中已知基因18 063個，預(yù)測的新基因828個。此外，共檢測出21 299個新的轉(zhuǎn)錄本。本研究挖掘了邊雞卵巢組織轉(zhuǎn)錄組信息，為開展邊雞就巢性狀相關(guān)基因的研究及分子調(diào)控機(jī)制研究奠定基礎(chǔ)。

對所有樣本間的相關(guān)性進(jìn)行層次聚類分析，用pearson法計算每兩個樣本間所有基因表達(dá)量的相關(guān)系數(shù)，相關(guān)系數(shù)越接近1，表明這兩個樣本間的基因表達(dá)的相似度越高。本研究發(fā)現(xiàn)，在165和220 d之間的相關(guān)系數(shù)為0.983，基因表達(dá)相似度最高，說明卵巢發(fā)育成熟后，開始排卵和排卵高峰兩個樣本卵巢組織基因表達(dá)相關(guān)性最強(qiáng)，與層次聚類分析表達(dá)譜結(jié)果一致。從邊雞不同生理階段卵巢組織差異表達(dá)基因數(shù)量顯示，存在上調(diào)基因和下調(diào)基因，165與220 d相比，差異表達(dá)基因的總量最少(656個)，與樣本間的相關(guān)性進(jìn)行層次聚類分析結(jié)果一致，這些差異表達(dá)基因可能促進(jìn)卵泡發(fā)育。

本研究使用DESeq2R包對轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行篩選，以樣本間的相關(guān)性和不同生理階段卵泡數(shù)為篩選基礎(chǔ)，找出了40個可能與卵泡發(fā)育相關(guān)的差異表達(dá)基因，qPCR隨機(jī)驗證所選基因的mRNA表達(dá)趨勢與轉(zhuǎn)錄組測序數(shù)據(jù)一致，證明了轉(zhuǎn)錄組測序的結(jié)果可靠。其中35個基因的表達(dá)量與發(fā)育卵泡數(shù)成正相關(guān)，5個基因表達(dá)量與發(fā)育卵泡數(shù)不成正相關(guān)(鐵氧還蛋白-1L(FDX1L)、N-Myc下游調(diào)控基因-4(NDRG4)、黑素皮質(zhì)素受體(MC3R)、催乳素釋放激素受體(PRLHR)和血管活性腸肽(VIP))。FDX1L基因編碼一種183個氨基酸的線粒體蛋白Fdx2，它在鐵-硫簇生物發(fā)生過程中起重要作用[17]。它們通常通過半胱氨酸的巰基側(cè)鏈與假體蛋白結(jié)合，并作為電子受體或供體，在關(guān)鍵的細(xì)胞活動中發(fā)揮重要作用[18]。本試驗中發(fā)現(xiàn)，F(xiàn)DX1L在70、220、165和330 d邊雞卵巢中表達(dá)量依次遞減，表明該基因可能參與卵巢發(fā)育和調(diào)控排卵過程，作用機(jī)

圖4 qPCR驗證不同日齡邊雞卵巢中的差異表達(dá)基因Fig.4 Differentially expressed genes in the ovary of Bian chicken at different ages verified by the qPCR

理需進(jìn)一步研究探討。NDRG4在多種癌癥中被認(rèn)為與腫瘤細(xì)胞增殖能力相關(guān)[19-20]，NDRG4 表達(dá)下調(diào)后卵巢癌細(xì)胞增殖能力增高、凋亡減少，其機(jī)制可能與調(diào)節(jié) Bcl-2/Bax 表達(dá)及 P38 /MAPK 信號通路活性有關(guān)[21-22]，NDRG4通過Akt/CREB激活促進(jìn)肌原性分化[23]，支持激活細(xì)胞外信號調(diào)節(jié)激酶(ERK)促進(jìn)神經(jīng)元分化，并減少參與神經(jīng)元生長的Elk-1的激活[24]?，F(xiàn)已有報道稱Bcl-2與鵝就巢相關(guān)[25]，雖在本試驗中未發(fā)現(xiàn)Bcl-2在不同卵泡發(fā)育階段的邊雞卵巢中的表達(dá)量呈差異性，但發(fā)現(xiàn)NDRG4在邊雞卵巢中的表達(dá)量隨日齡增加呈遞減趨勢，NDRG4基因下調(diào)可能促進(jìn)卵巢機(jī)能衰老，與上述研究中NDRG4基因下調(diào)后卵巢癌細(xì)胞增殖能力增高、凋亡減少的信號通路可能相關(guān)，卵巢機(jī)能衰老與卵巢癌細(xì)胞增殖能力的相關(guān)性需要進(jìn)一步研究探討。黑素皮質(zhì)素受體(melanocortin receptors，MCRs)是G-蛋白偶聯(lián)受體(G protein-coupled receptors，GPCRs)超家族成員之一[26]。人、鼠和雞MCRs由5個成員組成，分別是MC1R、MC2R、MC3R、MC4R和MC5R，其中，MC3R基因不僅在外周組織中有一定表達(dá)，而且在大腦中呈現(xiàn)高表達(dá)，因此，又被稱為神經(jīng)性黑素皮質(zhì)素受體基因[27-28]。現(xiàn)有研究發(fā)現(xiàn)，MC3R主要參與調(diào)控動物的體重和能量代謝的過程[29]，也有研究發(fā)現(xiàn)，MC3R在人的心、乳腺、骨骼肌和腎以及睪丸和卵巢中表達(dá)[30]，在這些組織中表達(dá)很可能與免疫功能相關(guān)[31]。有關(guān)MC3R表達(dá)對禽類繁殖性能的研究尚未見報道，本試驗中發(fā)現(xiàn)，MC3R基因在邊雞卵巢70和330 d表達(dá)量高于165和220 d的表達(dá)量，推測MC3R可能起到抑制卵泡發(fā)育的作用，能夠促進(jìn)邊雞卵巢就巢性，這一發(fā)現(xiàn)為進(jìn)一步研究卵泡發(fā)育及排卵調(diào)控機(jī)制奠定基礎(chǔ)。李國輝等[32]發(fā)現(xiàn)，催乳素(PRL)和神經(jīng)肽Y(NPY)基因作為雞繁殖性能的候選基因，為多基因聚合育種技術(shù)在家禽遺傳改良方面的應(yīng)用。在本試驗中雖未發(fā)現(xiàn)PRL，但發(fā)現(xiàn)催乳素釋放激素受體(PRLHR)在卵泡發(fā)育期即165和220 d低表達(dá)，分別為1.14和1.16，而在未發(fā)育期70 d和就巢期330 d表達(dá)量相對較高，分別為3.13和2.72，PRLHR與NPY受體(NPYR)結(jié)合，從而影響PRL釋放，進(jìn)而調(diào)控邊雞的就巢性。Hosseinpour 等[33]發(fā)現(xiàn)，血管活性腸肽(VIP)及其受體(VIPR-1)的SNP與火雞產(chǎn)蛋性能相關(guān)，但未報道與就巢性相關(guān)。本試驗發(fā)現(xiàn)，在70、165和220 d的邊雞卵巢中VIP呈差異性表達(dá)，但在就巢期即330 d的邊雞卵巢中未發(fā)現(xiàn)VIP表達(dá)，該基因表達(dá)可能參與卵巢發(fā)育，與Hosseinpour 等[33]的結(jié)果一致。

本研究通過GO功能分析，發(fā)現(xiàn)邊雞不同日齡卵巢組織差異表達(dá)基因功能主要集中在“生物過程”。在“生物過程”中，與發(fā)育繁殖相關(guān)的通路有發(fā)育過程、繁殖和繁殖過程，富集在這些通路的差異表達(dá)基因可能參與卵泡的發(fā)育和排卵調(diào)控機(jī)制；而邊雞卵巢成熟開始排卵期間，即165與220 d相比，與發(fā)育繁殖相關(guān)的差異表達(dá)基因僅57個且數(shù)量最少，與本研究中卵巢基因差異表達(dá)相關(guān)性分析結(jié)果一致，這些少量的差異表達(dá)基因可能參與發(fā)育成熟卵泡數(shù)量的調(diào)控；220與330 d相比，共涉及188個與發(fā)育繁殖相關(guān)的差異表達(dá)基因，這些差異表達(dá)基因可能與邊雞就巢性的調(diào)控相關(guān)。對不同生理階段的邊雞卵巢差異表達(dá)基因KEGG通路比較，均可發(fā)現(xiàn)這些差異基因在信號轉(zhuǎn)導(dǎo)、免疫系統(tǒng)以及信號分子與相互作用的路徑上富集最多。通過揭示這些差異表達(dá)基因的功能，為多基因聚合育種改良邊雞就巢性提供了可能。

4 結(jié) 論

本研究通過構(gòu)建基因表達(dá)譜，獲取邊雞卵巢組織轉(zhuǎn)錄組信息，篩選出40個與卵巢發(fā)育相關(guān)的差異表達(dá)基因，通過GO功能富集和KEGG分析，獲得了差異表達(dá)基因的功能分類和通路，為開展邊雞卵泡發(fā)育相關(guān)基因的研究及分子調(diào)控機(jī)制研究奠定基礎(chǔ)。