表達(dá)譜芯片篩選獺兔毛色相關(guān)候選差異表達(dá)基因

施麗娟，秦立志，王文洲，閆曉榮，翁巧琴，吳信生*

(1.揚(yáng)州大學(xué)動(dòng)物科學(xué)與技術(shù)學(xué)院，揚(yáng)州 225009；2.江蘇省動(dòng)物繁育與分子設(shè)計(jì)重點(diǎn)實(shí)驗(yàn)室，揚(yáng)州 225009；3.浙江省余姚市欣農(nóng)兔業(yè)有限公司，余姚 315400)

表達(dá)譜芯片篩選獺兔毛色相關(guān)候選差異表達(dá)基因

施麗娟1，2，秦立志1，2，王文洲1，2，閆曉榮1，2，翁巧琴3，吳信生1，2*

(1.揚(yáng)州大學(xué)動(dòng)物科學(xué)與技術(shù)學(xué)院，揚(yáng)州 225009；2.江蘇省動(dòng)物繁育與分子設(shè)計(jì)重點(diǎn)實(shí)驗(yàn)室，揚(yáng)州 225009；3.浙江省余姚市欣農(nóng)兔業(yè)有限公司，余姚 315400)

本研究旨在了解青紫藍(lán)色獺兔和白色獺兔毛色形成的分子遺傳學(xué)機(jī)理。試驗(yàn)分為3組，每組包含全同胞青紫藍(lán)獺兔和白色獺兔各1只。利用安捷倫基因芯片技術(shù)分析不同分組之間差異基因表達(dá)情況，并通過實(shí)時(shí)熒光定量PCR方法對部分差異顯著基因進(jìn)行驗(yàn)證。差異表達(dá)分析結(jié)果顯示，3組青紫藍(lán)色獺兔對白色獺兔分別篩選出545、1 680和2 092個(gè)差異表達(dá)基因 (Differentially expressed genes，DEGs)，共有172個(gè)共同基因；GO(Gene Ontology)分類結(jié)果顯示，這些基因主要參與細(xì)胞分化、增殖、凋亡、發(fā)育、離子轉(zhuǎn)運(yùn)、免疫應(yīng)答等生物學(xué)過程。通過實(shí)時(shí)熒光定量PCR技術(shù)對Creb1、Wnt2、Gnaq和Gna14基因在組織樣中的表達(dá)量進(jìn)行檢測，結(jié)果與表達(dá)譜芯片一致。通過表達(dá)譜芯片結(jié)果結(jié)合生物信息學(xué)分析和文獻(xiàn)報(bào)道，篩選出了Creb1、TH、Wnt2、Gnaq、Dvl1、Mapk12等對毛色形成可能具有重要影響的候選基因。

獺兔；毛色；差異表達(dá)基因；基因芯片

獺兔的毛皮柔軟美觀、輕便保暖，一直以來都備受人們的喜愛，特別是其自然形成的多種多樣的毛色，更能滿足不同消費(fèi)者的審美需求。培育多色獺兔既可以增加獺兔的特色附加價(jià)值，也可以提高養(yǎng)兔業(yè)的經(jīng)濟(jì)效益。在獺兔的選育工作中，一般通過雜交方法，組合不同的基因型，培育出多種色型的獺兔。青紫藍(lán)獺兔是獺兔的一個(gè)品種，全身被毛基部為瓦藍(lán)色，中段為珍珠灰色，毛尖部為黑色。頸部毛色略淺于體側(cè)部，背部毛色較深，腹部毛色呈淺藍(lán)或白色。因?yàn)槠浔幻愃朴谇嘧纤{(lán)兔，故得名。毛色遺傳作為可以利用的遺傳標(biāo)記，在人、鼠、豬等多種哺乳動(dòng)物中都有研究[1-2]。哺乳動(dòng)物的毛色是遺傳基礎(chǔ)復(fù)雜的質(zhì)量性狀，而后天的環(huán)境也會(huì)對毛色有一定的調(diào)控作用[3]。D.C.Bennett等[4]研究表明，小鼠體內(nèi)控制毛色形成的多個(gè)通路中包含大約127個(gè)基因。哺乳動(dòng)物的皮毛之所以出現(xiàn)從白到黑的多種顏色，主要是由體內(nèi)黑色素細(xì)胞產(chǎn)生的真黑素和棕黑素之間的相對數(shù)量與分布決定的[5-6]。高文玉[7]通過對獺兔毛色遺傳資料的收集整理，發(fā)現(xiàn)影響獺兔毛色的基因至少涉及10個(gè)系統(tǒng)或座位，其可通過等位基因和不同位點(diǎn)基因間的相互作用決定獺兔毛色的具體表達(dá)。許多研究已證實(shí)，Mc1r、Mitf、Kit、Asip、Mlph、Pmel-17等基因參與了黑色素的產(chǎn)生，并與哺乳動(dòng)物毛色的形成過程相關(guān)[8-13]。

自1991年基因芯片概念被首次提出開始[14]，基因芯片技術(shù)因其自動(dòng)化、高通量、微型化檢測等優(yōu)點(diǎn)，迅速成為生命科學(xué)界的研究熱點(diǎn)[15]，隨著各種相關(guān)技術(shù)的發(fā)展，基因芯片已成為“后基因組時(shí)代”基因功能分析研究的最重要技術(shù)之一?；蛐酒夹g(shù)是一項(xiàng)由分子生物學(xué)、微電子學(xué)、物理學(xué)、化學(xué)和計(jì)算機(jī)學(xué)等多學(xué)科交叉融合而成的高新技術(shù)，高通量篩選的特點(diǎn)使其在基因表達(dá)分析、疾病診斷和治療、新藥發(fā)現(xiàn)等眾多領(lǐng)域得到廣泛應(yīng)用[16]。本研究利用家兔全基因表達(dá)譜芯片技術(shù)篩選青紫藍(lán)色獺兔和白色獺兔背部皮膚組織差異表達(dá)基因，探索對毛色形成可能具有重要影響的功能基因，初步揭示相關(guān)基因的差異表達(dá)模式，為進(jìn)一步篩選調(diào)控獺兔體內(nèi)黑色素生成、多種多樣毛色形成相關(guān)的主效基因以及毛色的進(jìn)一步選育改良工作提供參考。

1 材料與方法

1.1 樣本采集

選取14日齡青紫藍(lán)色和白色獺兔各3只(由浙江省余姚市欣農(nóng)兔業(yè)有限公司提供)，來自于3對獺兔全同胞個(gè)體，其中每對全同胞個(gè)體均含有青紫藍(lán)色獺兔和白色獺兔各1只。取背中部靠近脊椎處約 1.5 cm2皮膚樣，剪去表層毛發(fā)后將其剪成小塊放入去 RNA 酶指形管，于液氮中保存?zhèn)溆谩?/p>

1.2 基因芯片

本試驗(yàn)采用Agilent 兔全基因組4*44K 芯片(design ID：020908)，購自上海伯豪生物芯片有限公司。

1.3 RNA提取和純化

采用Trizol法提取樣品的總RNA，所得總RNA 經(jīng)Agilent Bioanalyzer 2100(Agilent technologies，US)電泳質(zhì)檢合格，并使用RNeasy mini kit(QIAGEN，Germany)和RNase-Free DNase Set(QIAGEN，Germany)純化總RNA。

1.4 熒光標(biāo)記

采用Agilent表達(dá)譜芯片配套試劑盒提取試驗(yàn)樣品總RNA，按照Low Input Quick Amp Labeling Kit、One-Color(Agilent technologies，US)的要求，對樣品總RNA中的mRNA進(jìn)行放大和熒光標(biāo)記，并用RNeasy mini kit(QIAGEN，Germany)對標(biāo)記后的cRNA進(jìn)行純化。

1.5 芯片雜交

按照Agilent表達(dá)譜芯片所提供的雜交標(biāo)準(zhǔn)流程和配套試劑盒，將1.65 μg的cRNA樣在65℃滾動(dòng)雜交爐中10 r·min-1滾動(dòng)雜交17 h，并在洗缸中洗片，洗片所用的試劑為Gene Expression Wash Buffer Kit(Agilent technologies，US)。

1.6 芯片的洗滌及掃描

從雜交爐中取出已經(jīng)完成雜交的芯片，用 GE Wash Buffer 洗滌芯片，并用 Agilent Microarray Scanner 進(jìn)行掃描。用 Feature Extraction software 10.7 讀取數(shù)據(jù)，最后采用 Gene Spring Software 11.0 進(jìn)行歸一化處理，所用的算法為Quantile。

1.7 數(shù)據(jù)分析

將讀取的數(shù)據(jù)進(jìn)行歸一化處理后，以 Ratio≥2 或Ratio≤0.5 篩選差異表達(dá)基因。對于篩選出的未知基因，根據(jù)其登錄號在 NCBI 上搜索其功能注釋。

1.8 熒光定量PCR驗(yàn)證

利用與芯片同一批組織樣提取RNA進(jìn)行熒光定量PCR驗(yàn)證，隨機(jī)選取3個(gè)芯片結(jié)果中的差異表達(dá)基因，以GAPDH作為內(nèi)參基因，采用 SYBR GreenⅠ法進(jìn)行熒光定量試驗(yàn)，以驗(yàn)證芯片結(jié)果的可靠性，預(yù)混體系UltraSYBR Mixture (With ROX) 購自北京康為世紀(jì)生物科技有限公司。引物利用Oligo7軟件根據(jù)芯片探針相對應(yīng)的核酸序列設(shè)計(jì)(表1)。采用TaKaRa 公司的 PrimeSYB?Perfect Real Time 試劑盒制備cDNA，標(biāo)準(zhǔn)曲線由 cDNA 梯度稀釋后建立。每個(gè)樣品設(shè)置3個(gè)重復(fù)，利用ABI7500熒光定量PCR儀(Applied Biosystems，US)分析PCR過程中的熒光信號值，并獲得各樣品和基因的Ct值，以2-ΔΔCt方法計(jì)算基因表達(dá)的差異倍數(shù)。

表1 實(shí)時(shí)熒光定量PCR的引物序列

Table 1 Primer sequences for quantitative real-time PCR

基因名稱Genename引物序列(5'→3')Primersequence片段大小/bpThefragmentlengthGAPDHF:CGGCACGAAGCTGGTGAAAGA R:GCGATCCAGGAAGGTGAGGAAGT293Wnt2F:ACTTCAGGAAAACGGGCGAT R:CTGCCTCTCGGTCCCTGATA115Creb1F:AGACTTCAGCACCTGCCATC R:TGTCCATCAGTGGTCTGTGC140GnaqF:AGGCTCATGCACAATTAGTTCG R:TTGTTGCGTAGGCAGGTAGG212Gna14F:GGAAGAGCACCTTTATCAAGCAA R:TCTGCTCACACACGTAAGGGAT165

2 結(jié) 果

2.1 差異表達(dá)基因

本試驗(yàn)利用兔的表達(dá)譜芯片初步得到不同毛色全同胞獺兔皮膚組織表達(dá)譜，以每組的白色獺兔基因表達(dá)值為對照，青紫藍(lán)色獺兔基因的表達(dá)值分別與其對比，按照差異倍數(shù)≥2或≤0.5的標(biāo)準(zhǔn)篩選差異基因，3組分別篩選出了545、1 680和2 092個(gè)差異表達(dá)基因，經(jīng)韋恩分析后發(fā)現(xiàn)3組共有差異表達(dá)基因172個(gè)(圖1)。在第1組中，111個(gè)基因在青紫藍(lán)色獺兔皮膚組織中表達(dá)上調(diào)，61個(gè)基因表達(dá)下調(diào)；在第2組中，107個(gè)基因在青紫藍(lán)色獺兔皮膚組織中表達(dá)上調(diào)，65個(gè)基因表達(dá)下調(diào)；在第3組中，92個(gè)基因在青紫藍(lán)色獺兔皮膚組織中表達(dá)上調(diào)，80個(gè)基因在青紫藍(lán)色獺兔皮膚組織中表達(dá)下調(diào)。172個(gè)基因中有89個(gè)在3組中表達(dá)趨勢一致，其中62個(gè)基因在青紫藍(lán)色獺兔皮膚組織中表達(dá)上調(diào)，27個(gè)基因在青紫藍(lán)色獺兔皮膚組織中表達(dá)下調(diào)。表2列出了其中的20個(gè)差異表達(dá)基因。

2.2 GO分類結(jié)果

篩選得到的差異基因經(jīng)過 GO(http://www.geneontology.org/)進(jìn)行功能分類，結(jié)果見表3。這些差異基因參與生物學(xué)過程的調(diào)控、細(xì)胞、代謝等生物學(xué)過程(圖2)。其中生物學(xué)過程的調(diào)控、生物學(xué)調(diào)節(jié)、細(xì)胞過程、代謝過程占10%以上，多細(xì)胞有機(jī)體過程、發(fā)育過程、定位、建立定位、一體化過程都占5%～10%，其他的生物學(xué)過程都占5%以下。

圖1 3組差異基因的Venn.分析Fig.1 Venn.analysis of the DEGs in 3 groups

2.3 相關(guān)差異表達(dá)基因信號通路富集

172個(gè)共同差異表達(dá)基因中，大部分基因處于已研究發(fā)現(xiàn)的KEGG、Biocarta信號通路，他們存在于多條KEGG、Biocarta信號通路中，例如Creb1基因處于絲裂素活化蛋白激酶信號通路(MAPKinase Signaling Pathway)等12個(gè)Biocarta信號通路，以及黑素生成(Melanogenesis)等5個(gè)KEGG信號通路。本試驗(yàn)找出了一些與色素形成相關(guān)的信號通路以及相關(guān)的差異基因(表4)。

表2 3組中表達(dá)趨勢一致的20個(gè)差異基因

Table 2 20 common differentially expressed genes in 3 groups

基因名稱GenesnameGenBank登錄號GenBankaccessionnumber差異倍數(shù)Foldchange第1組Group1第2組Group2第3組Group3Creb1NC_013675.13.382.352.68GnaqNC_013669.13.752.723.73THAF493546.13.072.812.66Gna14NC_013669.12.122.252.01Camk2aNC_013671.12.272.453.36Mapk12NM_013871.35.392.562.85Actn2NC_013684.12.052.032.02Dvl1NM_010091.32.362.122.56ST1A8NC_013674.13.163.273.18TtnNC_013675.15.405.493.34Reg4NC_013681.12.064.062.85Htr5bNC_013675.14.162.052.65Zfp11NM_172462.42.072.160.28Igsf1NC_013690.12.063.353.95Adrbk2NM_177078.32.073.802.23SctrNC_013675.10.320.290.12CgaNC_013680.10.310.460.08Wnt2NM_001171042.10.310.410.38Dnahc6NM_001370.10.110.460.42EgfrNC_013678.10.470.460.36

表3 分子生物學(xué)過程相關(guān)基因GO分類

Table 3 GO term to gene in biological process

基因功能分類Thefunctionalclassificationofgenes基因數(shù)NumberofgenesGO分類號GONo.生物學(xué)過程的調(diào)節(jié)Regulationofbiologicalprocess1320050789生物學(xué)過程的負(fù)調(diào)控Negativeregulationofbiologicalprocess350048519生物學(xué)過程的正調(diào)控Positiveregulationofbiologicalprocess380048518生物學(xué)調(diào)節(jié)Biologicalregulation1350065007一體化過程Multi-organismprocess700051704建立定位Establishmentoflocalization850051234定位Localization900051179刺激應(yīng)答Responsetostimulus420050896運(yùn)動(dòng)Locomotion320040011生長Growth40040007發(fā)育過程Developmentalprocess900032502繁殖Reproduction30000003免疫系統(tǒng)過程Immunesystemprocess480002376代謝過程Metabolicprocess1330008152細(xì)胞生理過程Cellularprocess1380009987解剖結(jié)構(gòu)形成Anatomicalstructureformation200010926凋亡Death190016265生殖過程Reproductiveprocess80022414生物黏附Biologicaladhesion120022610多細(xì)胞有機(jī)體過程Multicellularorganismalprocess1120032501

圖2 分子生物學(xué)過程相關(guān)基因GO分類Fig.2 GO term to gene in biological process

表4 色素形成相關(guān)的信號通路

Table 4 Signal pathway related to pigment formation

通路名稱Nameofpathway差異表達(dá)基因Differentiallyexpressedgenes黑素生成MelanogenesisWnt2、Gnaq、Dvl1、Creb1黑素細(xì)胞發(fā)展和色素通路MelanocytedevelopmentandpigmentationpathwayCreb1酪氨酸代謝TyrosinemetabolismTH絲裂素活化蛋白激酶信號通路MAPkinasesignalingpathwayMapk12、Creb1

2.4 熒光定量PCR驗(yàn)證差異基因

為驗(yàn)證芯片結(jié)果的可靠性，本研究從芯片雜交結(jié)果中選取4 個(gè)基因進(jìn)行熒光定量 PCR驗(yàn)證。芯片和實(shí)時(shí)定量 PCR 的結(jié)果比較見表5。結(jié)果顯示，Creb1、Wnt2、Gnaq和Gna14 各基因的表達(dá)譜芯片和熒光定量兩個(gè)試驗(yàn)所得的表達(dá)倍數(shù)雖然不同，但二者的表達(dá)趨勢基本相同(圖3)，因此認(rèn)為芯片的篩查結(jié)果可靠。

3 討 論

表達(dá)譜芯片是基因芯片的一個(gè)重要分支，其方法已經(jīng)相對成熟，并研制出許多物種的商品化表達(dá)譜芯片產(chǎn)品[17]。楊靜等[18]利用表達(dá)譜芯片技術(shù)篩選子宮頸上皮內(nèi)瘤變發(fā)展成子宮頸癌危險(xiǎn)性相關(guān)的差異表達(dá)基因發(fā)掘了ANKRD20A1、OLFM4、SNORD等13個(gè)在子宮頸癌發(fā)生的早期階段發(fā)揮重要作用的基因。李挺[19]通過表達(dá)譜芯片技術(shù)找到了新西蘭白兔與福建黃兔背肌和腿肌中共有的Tnnc1、Myo5b、Tpm3等7個(gè)影響兔肌肉生長發(fā)育的基因。本研究利用Agilent兔全基因組表達(dá)譜芯片篩選出了青紫藍(lán)色獺兔相對于白色獺兔的差異表達(dá)基因，并且選取部分基因進(jìn)行熒光定量PCR驗(yàn)證，結(jié)果顯示，芯片結(jié)果與熒光定量結(jié)果中基因表達(dá)量基本一致，說明利用表達(dá)譜芯片技術(shù)篩選毛色性狀相關(guān)候選基因是可靠的途徑。

本試驗(yàn)利用GO功能分類對篩選的差異表達(dá)基因進(jìn)行生物學(xué)過程分析，結(jié)果顯示，差異表達(dá)基因主要參與細(xì)胞黏附、增殖、分化及凋亡等細(xì)胞過程，以及發(fā)育、離子轉(zhuǎn)運(yùn)等生物學(xué)過程。本研究通過表達(dá)譜芯片結(jié)果結(jié)合生物信息學(xué)分析和文獻(xiàn)報(bào)道，篩選出TH、Creb1、Wnt2、Dvl1、Gnaq、Mapk12等與獺兔毛色形成過程相關(guān)性較大的候選基因，與前人的研究結(jié)果一致。本研究中酪氨酸羥化酶基因(Tyrosine hydroxylase，TH)在3組不同毛色獺兔表達(dá)譜芯片中的差異表達(dá)倍數(shù)均大于2，提示該基因可能與獺兔毛色形成相關(guān)，其編碼的酪氨酸羥化酶是一種以酪氨酸為底物，并催化酪氨酸形成多巴的酶，是多巴和兒茶酚胺生物合成過程中的關(guān)鍵步驟和第一個(gè)限速酶[20]。 韓宇等[21]認(rèn)為，TH基因控制家蠶體內(nèi)黑色素的形成。本研究中cAMP反應(yīng)元件結(jié)合蛋白基因(cAMP response element binding protein，Creb1)在3組不同毛色獺兔表達(dá)譜芯片中也具有差異表達(dá)。cAMP反應(yīng)元件結(jié)合蛋白是G蛋白耦合受體信號(GPCR)中最重要的一個(gè)核轉(zhuǎn)錄因子。卜今[22]的研究表明，丹皮酚能通過活化JNK/SAPK途徑，激活下游的酶性蛋白，抑制磷酸化Creb，進(jìn)而下調(diào)Mitf基因，減少Tyr基因的表達(dá)，抑制黑素合成。Creb1基因是該Creb家族基因中非常重要的一員，E.R.Price等[23]研究發(fā)現(xiàn)，在正常的黑色素細(xì)胞中，色素沉著信號從Mc1r基因傳遞到Mitf基因需要Creb1基因的參與。在篩選得到的差異基因中，Wnt2是Wnt家族的成員之一，Wnt/β-catenin通路是一條保守的信號通路，廣泛存在各種動(dòng)物細(xì)胞內(nèi)。Wnt家族是一個(gè)非常龐大的多基因家族，主要參與細(xì)胞的增殖與分化過程。Wnt可以與細(xì)胞表面的Frizzled受體結(jié)合，通過GPCR激活dishevelled(DVL)蛋白從而調(diào)節(jié)β-catenin[24]。Dvl在Wnt信號通路激活過程中，介導(dǎo)Wnt信號傳遞到不同的途徑，起著信號分流作用。G.M.Boland等[25]研究發(fā)現(xiàn)，過量表達(dá)Dvl1基因能夠使β-catenin在胞質(zhì)中穩(wěn)定積累。于秋菊等[26]研究表明，β-catenin在棕色羊駝皮膚組織中的基因表達(dá)量高于白色羊駝，可能與羊駝毛色形成具有相關(guān)性。K.Pham等[27]研究發(fā)現(xiàn)，Wnt2在惡性黑色素瘤中過表達(dá)。Gnaq基因編碼的Gnaq蛋白是異源三聚體G蛋白中a亞基下游中的亞基基因，其功能是將信息從細(xì)胞表面受體傳遞到細(xì)胞內(nèi)，在黑素細(xì)胞內(nèi)皮素B受體中發(fā)揮重要作用，存在于包括皮膚在內(nèi)的各種組織中。M.A.Davies等[28]研究指出，黑色素細(xì)胞中Gnaq蛋白的突變體和其他一些基因協(xié)作誘導(dǎo)小鼠黑色素瘤的形成。C.D.Van Raamsdonk等[29]研究表明，Gnaq超效等位基因的突變會(huì)引起真皮內(nèi)黑色素細(xì)胞增加從而導(dǎo)致皮膚色素?cái)U(kuò)散。Mapk12基因編碼的絲裂素活化蛋白激酶(MAPK)參與生物體內(nèi)的多種受體信號傳遞途徑，C.D.Van Raamsdonk等[30]研究表明，MAPK信號通路作用于黑色素小體的轉(zhuǎn)運(yùn)，可能在眼色素瘤和色素痣中發(fā)揮作用。

表5 熒光定量PCR與表達(dá)譜芯片的基因表達(dá)比較

Table 5 The comparison of gene expression level between quantification RT-PCR and chip

基因Genes組別Groups表達(dá)譜芯片差異倍數(shù)Foldchangechip實(shí)時(shí)熒光定量PCR差異倍數(shù)qRT-PCRfoldchangeCreb113.384.3522.352.2832.683.23Wnt210.310.8320.410.3330.380.42Gnaq13.752.0522.723.8533.732.44Gna1412.121.8822.252.5632.011.92

圖3 芯片與熒光定量PCR結(jié)果比對Fig.3 The comparisons on analysis result of microarray and qRT-PCR

4 結(jié) 論

本研究利用基因芯片技術(shù)篩選出同一窩出生的青紫藍(lán)色獺兔和白色獺兔中差異表達(dá)基因，共篩選出172個(gè)基因，參與了細(xì)胞的增殖、分化、凋亡，信號傳遞、離子轉(zhuǎn)運(yùn)、催化酶的活性、蛋白結(jié)合等多種細(xì)胞生物學(xué)過程。其中Creb1、Gnaq、Mapk12、TH、Wnt2、Dvl1基因的差異表達(dá)，表明這些基因可能是對獺兔毛色形成具有重要影響的候選基因。

[1] MARIE M，VERA S D，RICARDO M，et al.The developmental role of agouti in color pattern evolution [J].Science，2011，331(25)：1062-1065.

[2] STURM R A.Molecular genetics of human pigmentation diversity[J].HumMolGenet，2009，18(1)：9-17.

[3] TEASDELE R D，STURM R A，BOX N F.Human pigmentation genes：Identification，structure and consequences of polymorphic variation[J].Gene，2001，277：49-62.

[4] BENNETT D C，LAMOREUX M L.The color loci of mice-a genetic century[C].Pigment cell research / sponsored by the European Society for Pigment Cell Research and the International Pigment Cell Society，2003，16(4)：333-344.

[5] CONE R D，LU D，KOPPULA S，et al.The melanocortin receptors：agonists，antagonists，and the hormonal control of pigmentation[J].RecentProgHormRes，1996，51：287-317，318.

[6] HARTMEYER M，SCHOLZEN T，BECHER E，et al.Human dermal microvascular endothelial cells express the melanocortin receptor type 1 and produce increased levels of IL-8 upon stimulation with alpha-melanocyte-stimulating hormone[J].JImmunol，1997，15(4)：1930-1937.

[7] 高文玉．關(guān)于家兔毛色遺傳規(guī)律的探討分析[J]．中國農(nóng)學(xué)通報(bào)，2012，28(8)：19-23． GAO W Y.Study and analysis on the genetic regulation of rabbit hair color[J].ChineseAgriculturalScienceBulletin，2012，28(8)：19-23.(in Chinese)

[8] KIJAS J M H，WALES R，TORNSTEN A，et al.Melanocortin receptor 1 (MC1R) mutations and coat color in pigs[J].Genetics，1998，150：1177-1185.

[9] BAUER G L，PRAETORIUS C，BERGSLEINSDOTTIR K，et al.The role of MITF phosphorylation sites during coat color and eye development in mice analyzed by bacterial artificial chromosome transgene rescue[J].Genetics，2009，183(2)：581-594.

[10] MARKLUND S，MOLLER M，SANDBERG K，et al.Close association between sequence polymorphism in the KIT gene and the roan coat color in horses[J].MammGenome，1999，10(3)：283-288.

[11] MUNDY N I，KELLY J.Investigation of the role of the agouti signaling protein gene (ASIP) in coat color evolution in primates[J].MammGenome，2006，17(12)：1205-1213.

[12] PHILIPP U，HAMANN H，MECKLENBURG L，et al.Polymorphisms within the canine MLPH gene are associated with dilute coat color in dogs[J].BMCGenet，2005，6：34.

[13] BRUNBERG E，ANDERSSON L，COTHRAN G，et al.A missense mutation in PMEL17 is associated with the Silver coat color in the horse[J].BMCGenet，2006，7：46.

[14] FODOR S P A，READ J L，PIRRUN G M C.Light directed spatially addressable parallel chemical synthesis[J].Science，1991，251：767-773.

[15] CHEN C P，CHIANG M C，WANG T H，et al.Microvillus inclusion disease：prenatal ultrasound findings，molecular diagnosis and genetic counseling of congenital diarrhea[J].TaiwanJObstetGyne，2010，49(4)：487-494.

[16] 安利峰，勝 利．基因芯片及其應(yīng)用進(jìn)展[J]．西北民族學(xué)院學(xué)報(bào)，2002，23(2)：40-42． AN L F，SHENG L.Gene core and its application progress[J].JournalofNorthwestMinoritiesUniversity，2002，23(2)：40-42.(in Chinese)

[17] POLATI R，MENINI M，ROBOTTI E，et al.Proteomic changes involved in tenderization of bovine longissimus dorsi muscle during prolonged ageing[J].FoodChem，2012，135(3)：2052-2069.

[18] 楊 靜，馮家成，孟 曉，等．基因表達(dá)譜芯片在預(yù)測子宮頸上皮內(nèi)瘤變進(jìn)展風(fēng)險(xiǎn)中應(yīng)用[J]．臨床與實(shí)驗(yàn)病理學(xué)雜志，2013，29(1)：53-57． YANG J，F(xiàn)ENG J C，MENG X，et al.The role of gene expression profile chip in the prediction of cervical intraepithelial neoplasia to cervical cancer[J].ChineseJournalofClinicalandExperimentalPathology，2013，29(1)：53-57.(in Chinese)

[19] 李 挺．利用表達(dá)譜芯片技術(shù)篩選獺兔肌肉發(fā)育相關(guān)基因[D]．揚(yáng)州：揚(yáng)州大學(xué)，2013． LI T.Screening genes for muscle of rabbit by cDNA microarray[D].Yangzhou：Yangzhou University，2013.(in Chinese)

[20] LEVITT M，SPECTOR S，SJOERDSMA A，et al.Elucidation of the rate-limiting step in norepinephrine biosynthesis in the perfused guinea-pig heart[J].JPharmacolExpTher，1965，148：1-8.

[21] 韓 宇，劉 春，魏麗婉，等．家蠶酪氨酸羥化酶基因BmTh的表達(dá)及功能[J]．昆蟲學(xué)報(bào)，2010，53(9)：949-956. HAN Y，LIU C，WEI L W，et al.Expression profile and function of tyrosine hydroxylase gene BmTh in the silkworm[J].ActaEntomologicaSinica，2010，53(9)：949-956.(in Chinese)

[22] 卜 今．丹皮酚抑制黑素合成相關(guān)分子機(jī)制的研究及KC-MC共培養(yǎng)體系中細(xì)胞因子對黑素合成的反饋調(diào)節(jié)[D]．北京：北京協(xié)和醫(yī)院，2008． PU J.A study on the molecular mechanism of inhibition on melanogenesis induced by paeonol and the feedback regulation on melanin synthesis by cytokines in KC-MC co-culture system[D].Beijing：Peking Union Medical Hospital，2008.(in Chinese)

[23] PRICE E R，HORSTMANN M A，WELLS A G，et al.Fisher DE alpha-melanocyte-stimulating hormone signaling regulates expression of microphthalmia，a gene deficient in Waardenburg syndrome[J].JBiolChem， 1998，273：33042-33047.

[24] DIJKSTERHUIS J P，PETERSEN J，SCHULTE G.WNT/Frizzled signalling：receptor-ligand selectivity with focus on FZD-G protein signalling and its physiological relevance：IUPHAR review 3[J].BrJPharmacol，2014，171(5)：1195-1209.

[25] BOLAND G M，PERKINS G，HALL D J，et al.Wnt 3 apromotes proliferation and suppresses osteogenic differentiation of adult human mesenchymal stem cell[J].CellBiochem，2004，93(6)：1210-1230.

[26] 于秀菊，董常生，范闊海，等．β-catenin在不同毛色羊駝皮膚中的表達(dá)和定位[J]．畜牧獸醫(yī)學(xué)報(bào)，2010，41(3)：335-340． YU X J，DONG C S，F(xiàn)AN K H，et al.Study on the expression and localization of β-catenin in alpaca skin of different colors[J].ActaVeterinariaetZootechnicaSinica，2010，41(3)：335-340.(in Chinese)

[27] PHAM K，MILOVANOVIC T，BARR R J，et al.Wnt ligand expression in malignant melanoma：pilot study indicating correlation with histopathological features[J].MolPathol，2003，56(5)：280-285.

[28] DAVIES M A，GERSHENWALD J E.Targeted therapy for melanoma：a primer [J].SurgOncolClinNAm，2011，20(1)：165-180.

[29] VAN RAAMSDONK C D，F(xiàn)ITCH KR，F(xiàn)UCHS H，et al.Effects of G-protein mutations on skin color[J].NatGenet，2004，36：961-968.

[30] VAN RAAMSDONK C D，BEZROOKOVE V，GREEN G，et al.Frequent somatic mutations of GNAQ in uveal melanoma and blue naevi[J].Nature，2009，457(7229)：599-602.

(編輯 程金華)

Screening Candidate Differentially Expressed Genes of Rex Rabbits’ Fur Color with Gene Chips

SHI Li-juan1，2，QIN Li-zhi1，2，WANG Wen-zhou1，2，YAN Xiao-rong1，2，WENG Qiao-qing3，WU Xin-sheng1，2*

(1.CollegeofAnimalScienceandTechnology，YangzhouUniversity，Yangzhou225009，China；2.JiangsuKeyLaboratoryofAnimalGenetics&BreedingandMolecularDesign，Yangzhou225009，China；3.ZhejiangYuyaoXinnongRabbitIndustryCo.，Ltd.，Yuyao315400，China)

In order to explore the molecular genetic mechanism of fur color formation between chinchilla rex rabbits and white rex rabbits.The experiments were divided into 3 groups，every group contained one chinchilla rex rabbit and one white rex rabbit，which came from full-sib family.The Agilent’s rabbit gene expression microarray was used to determine the differentially expression genes between different fur color groups of rex rabbits，and the expression patterns of selected differentially expressed genes were further identified by quantitative real-time PCR.The results showed that 545，1 680 and 2 092 differentially expressed genes associated with fur color of rex rabbits were screened out from 3 groups of chinchilla rex rabbits and white rex rabbits，respectively.A total of 172 mutual genes were found among these 3 rex rabbits groups.GO (Gene Ontology) analysis results revealed that the differentially expressed genes were related to cell differentiation，proliferation，apoptosis，growth，ionic transport altered and immune response，etc.The quantitative real-time PCR showed gene expression trend ofCreb1，Wnt2，GnaqandGna14 genes was consistent with that of gene expression microarray data，indicating that the result of gene chip was reliable.Combined with bioinformatic analysis and related literatures，some genes were potentially related to fur color formation，such asCreb1，TH，Wnt2，Gnaq，Dvl1，Mapk12，and so on.

rex rabbit；fur color；DEGs；gene chip

10.11843/j.issn.0366-6964.2015.04.009

2014-09-05

現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)資金(CARS-44-1)；江蘇省高校優(yōu)勢學(xué)科建設(shè)工程資助項(xiàng)目(PAPD-2014-134)

施麗娟(1990-)，女，江蘇南通人，碩士生，主要從事動(dòng)物遺傳育種與繁殖研究，E-mail:msbdmy@126.com

*通信作者：吳信生，教授，博士，E-mail:xswu@yzu.edu.cn

Q81；S829.1

A

0366-6964(2015)04-0568-08