重組人BMP-2誘導(dǎo)C3H10T1/2間質(zhì)干細(xì)胞定向成骨分化早期基因表達(dá)譜分析

左長清 汪宗桂 鐘月春 盧旱云 戴忠 吳鐵 崔燎

1.廣東醫(yī)學(xué)院藥理學(xué)教研室，廣東東莞523808；2.廣東醫(yī)學(xué)院生物化學(xué)教研室，廣東東莞523808

重組人BMP-2誘導(dǎo)C3H10T1/2間質(zhì)干細(xì)胞定向成骨分化早期基因表達(dá)譜分析

左長清1汪宗桂2鐘月春1盧旱云1戴忠1吳鐵1崔燎1

1.廣東醫(yī)學(xué)院藥理學(xué)教研室，廣東東莞523808；2.廣東醫(yī)學(xué)院生物化學(xué)教研室，廣東東莞523808

目的研究間質(zhì)干細(xì)胞早期定向成骨分化基因表達(dá)譜，為研究基因?qū)υ缙诔晒嵌ㄏ蚍只{(diào)控機制提供實驗基礎(chǔ)。方法分別提取重組人骨形成蛋白2(rhBMP.2)誘導(dǎo)組和對照組C3H10T1/2細(xì)胞總RNA，進(jìn)行擴(kuò)增標(biāo)記后，與ArraySTAR小鼠基因芯片雜交，應(yīng)用生物信息學(xué)軟件GeneSpring和GATHER對基因芯片數(shù)據(jù)進(jìn)行分析。應(yīng)用STRING在線軟件對差異表達(dá)基因構(gòu)建蛋白互作網(wǎng)絡(luò)并進(jìn)行網(wǎng)絡(luò)分析。結(jié)果C3H10T1/2早期成骨分化中，主要富集發(fā)育、器官形成等分子功能本體以及細(xì)胞因子.細(xì)胞因子受體作用信號通路。成骨分化1 d和4 d均上調(diào)表達(dá)基因42個，下調(diào)表達(dá)基因45個。網(wǎng)絡(luò)分析研究表明：Egfr、Cxcl12等信號分子參與調(diào)控rhBMP.2誘導(dǎo)成骨分化。結(jié)論篩選的差異表達(dá)基因和信號分子對早期成骨分化調(diào)控具有重要作用，為進(jìn)一步全面解析早期成骨定向分化提供實驗基礎(chǔ)。

早期成骨分化；基因表達(dá)譜；間質(zhì)干細(xì)胞；生物信息學(xué)

成骨細(xì)胞主要由骨髓間質(zhì)干細(xì)胞(mesenchymal stem cells，MSCs)分化而來，在其發(fā)育分化的整個過程中，不同的基因嚴(yán)格按照特定的時間順序開啟和關(guān)閉，表現(xiàn)為與分化的階段相一致，支持并且控制著成骨細(xì)胞的分化與成熟。因此，骨髓間質(zhì)干細(xì)胞定向成骨分化調(diào)控研究在國內(nèi)外受到高度重視，促進(jìn)間質(zhì)干細(xì)胞成骨定向分化與成熟，將是治療骨質(zhì)疏松癥的新的有效方法和手段。但是，目前關(guān)于間質(zhì)干細(xì)胞早期成骨定向分化方面的調(diào)控機制尚未完全闡明。本實驗室采用重組人骨形成蛋白2(rhBMP.2)誘導(dǎo)小鼠間質(zhì)干細(xì)胞C3H10T1/2定向成骨分化細(xì)胞模型，采用小鼠全基因組芯片篩選鑒定差異表達(dá)基因譜并進(jìn)行生物信息學(xué)分析，探討rhBMP.2成骨分化調(diào)控的差異表達(dá)基因集，并進(jìn)一步通過網(wǎng)絡(luò)生物學(xué)分析，為探討間質(zhì)干細(xì)胞早期成骨定向分化提供實驗依據(jù)。

1 材料與方法

1.1 材料

小鼠間質(zhì)干細(xì)胞C3H10T1/2購自中國科學(xué)院上海細(xì)胞庫。rhBMP.2(R&D systems)，ArraySTAR mouse lncRNA microarray(8×60 K，上海康成生物，芯片包含31423條蛋白編碼mRNA和25376條lncRNA，本文僅分析編碼蛋白mRNA表達(dá)變化)。BCIP/NBT堿性磷酸酶染色試劑盒(Sigma)。

1.2 實驗方法

1.2.1 堿性磷酸酶染色鑒定早期成骨定向分化C3H10T1/ 2接種于6孔板中，細(xì)胞生長至約90%融合，rhBMP.2誘導(dǎo)組加入含200 ng/mL rhBMP.2培養(yǎng)基，對照組加入不含rhBMP.2的培養(yǎng)基，每2～3天更換一次培養(yǎng)基。當(dāng)成骨誘導(dǎo)分化7 d后，rhBMP.2誘導(dǎo)組、對照組的單層細(xì)胞用PBS清洗2次，然后用體積分?jǐn)?shù)為0.7乙醇室溫固定15 min，BCIP/NBT避光染色30 min，用蒸餾水漂洗細(xì)胞3次，倒置顯微鏡下觀察并拍照記錄。

1.2.2 C3H10T1/2早期成骨定向分化基因芯片檢測及其生物信息學(xué)分析C3H10T1/2成骨定向分化1 d和4 d后，rhBMP.2誘導(dǎo)組、對照組采用TRIzol提取總RNA，按照說明擴(kuò)增、標(biāo)記熒光、純化、與芯片雜交、洗脫。用Agilent Scanner G2505C掃描每張芯片的的信號強度，應(yīng)用Agilent Feature Extraction軟件(version 11.0.1.1)分析芯片圖像文件，用GeneSpring GX v11.5.1軟件(Agilent)進(jìn)行數(shù)據(jù)處理分析，根據(jù)基因表達(dá)水平篩選出差異2倍以上的基因并采用GATHER軟件［1］(http：//gather.genome.duke.edu/)進(jìn)行生物信息學(xué)分析富集基因本體(gene ontology，GO)和生物通路。

1.2.3 差異表達(dá)基因蛋白互作網(wǎng)絡(luò)分析STRING (http：//string.db.org/)軟件［2］從各種資源中整合了已知的和預(yù)測的蛋白質(zhì)相互作用(包括直接和間接作用)來生成網(wǎng)絡(luò)，其數(shù)據(jù)來自于基因組、高通量實驗、共表達(dá)、已有知識4種資源。本研究采用此數(shù)據(jù)庫資源對上述差異表達(dá)基因進(jìn)行蛋白互作網(wǎng)絡(luò)構(gòu)建。

2 結(jié)果

2.1 間質(zhì)干細(xì)胞C3H10T1/2成骨分化鑒定

間質(zhì)干細(xì)胞C3H10T1/2經(jīng)rhBMP.2成骨誘導(dǎo)分化7 d后，堿性磷酸酶染色顯示：rhBMP.2誘導(dǎo)組呈藍(lán)紫色(圖1)，這表明rhBMP.2能誘導(dǎo)C3H10T1/2間質(zhì)干細(xì)胞定向成骨分化。此細(xì)胞模型可用于后續(xù)實驗研究。

圖1 C3H10T1/2間質(zhì)干細(xì)胞成骨誘導(dǎo)7 d堿性磷酸酶染色（100×）

2.2 C3H10T1/2基因表達(dá)譜生物信息學(xué)分析

通過對C3H10T1/2成骨分化1 d和4 d芯片結(jié)果進(jìn)行分析，篩選出兩個時間點表達(dá)均升高2倍以上的基因42個；表達(dá)均降低2倍以上的基因45個(表1)。采用生物軟件進(jìn)行分子功能基因本體和通路分析發(fā)現(xiàn)，顯著性差異分子功能本體為發(fā)育、器官形成等，并與細(xì)胞因子.細(xì)胞因子受體作用通路相關(guān)。

表1 C3H10T1/2細(xì)胞成骨分化早期差異表達(dá)基因集

2.3 一致性差異表達(dá)基因編碼蛋白相互作用網(wǎng)絡(luò)分析

采用STRING網(wǎng)絡(luò)資源對差異表達(dá)基因編碼蛋白進(jìn)行網(wǎng)絡(luò)分析，構(gòu)建了部分基因的基因表達(dá)調(diào)控與蛋白互作網(wǎng)絡(luò)，通過對網(wǎng)絡(luò)進(jìn)行分析(圖2)，發(fā)現(xiàn)經(jīng)rhBMP.2誘導(dǎo)后差異表達(dá)的Egfr、Cxcl12等信號分子處于蛋白網(wǎng)絡(luò)重要交叉點，可能與信號通路交談有關(guān)，說明在rhBMP.2誘導(dǎo)C3H10T1/2間質(zhì)干細(xì)胞早期成骨分化過程中，Egfr、Cxcl12等信號分子參與調(diào)控BMP.2信號。

3 討論

間質(zhì)干細(xì)胞成骨分化是一個多基因參與的，多步驟的復(fù)雜的生物過程。機體一些重要的轉(zhuǎn)錄因子如Runx2、Dlx2等［3.4］以及眾多的信號通路在成骨分化起到非常重要調(diào)控作用［5］。但是，關(guān)于BMP.2誘導(dǎo)間質(zhì)干細(xì)胞早期成骨定向分化相關(guān)調(diào)控機制以及不同信號間的交談分子，目前尚未完全闡明。

基因芯片技術(shù)提供了一個用于同時研究一個細(xì)胞或組織全基因組表達(dá)的較好平臺。本研究對rhBMP.2誘導(dǎo)C3H10T1/2成骨分化1 d和4 d后的基因表達(dá)譜進(jìn)行分析，提取了rhBMP.2誘導(dǎo)后差異表達(dá)的基因共87個，其中上調(diào)基因42個，下調(diào)基因45個。通過文獻(xiàn)分析發(fā)現(xiàn)，多個基因已經(jīng)報道與成骨分化緊密相關(guān)，它們在成骨定向、分化、轉(zhuǎn)分化過程中發(fā)揮著重要作用，如Frizzled家族成員FZD9［6］、同源蛋白家族基因Dlx2［7］、轉(zhuǎn)錄因子ID4［8］、細(xì)胞外基質(zhì)相關(guān)基因Ecm1［9］等。其他差異基因，將成為筆者研究BMP.2成骨信號的重要切入點。同時，進(jìn)行分子功能相關(guān)基因本體分析發(fā)現(xiàn)，在早期成骨分化過程，富集得分最高的為發(fā)育、器官形成、細(xì)胞因子.細(xì)胞因子受體作用相關(guān)信號，這可能與BMP.2誘導(dǎo)成骨分化密切相關(guān)。

復(fù)雜網(wǎng)絡(luò)理論作為一門新興學(xué)科，為系統(tǒng)水平上研究生物網(wǎng)絡(luò)提供了新的理論依據(jù)和平臺。2000年Jeong等［10］在《Nature》上第一次發(fā)表利用復(fù)雜網(wǎng)絡(luò)理論研究代謝網(wǎng)絡(luò)的拓?fù)涮匦缘恼撐模源艘院螅脧?fù)雜網(wǎng)絡(luò)理論研究各種生物網(wǎng)絡(luò)迅速發(fā)展［11.12］，主要集中在探索活細(xì)胞中所有分子的行為和它們之間的相互作用，以及從局部到整體，生物體如何形成了極其復(fù)雜的動力學(xué)機制。蛋白.蛋白互作網(wǎng)絡(luò)和轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò)是生物信息調(diào)控的主要實現(xiàn)方式，在整個生命過程中有著非常重要的作用，是決定細(xì)胞命運的關(guān)鍵因素。通過對差異基因構(gòu)建蛋白互作和轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò)，本研究發(fā)現(xiàn)EGFR、CXCL12等信號分子在rhBMP.2誘導(dǎo)間質(zhì)干細(xì)胞早期成骨分化網(wǎng)絡(luò)中處于非常重要的網(wǎng)絡(luò)節(jié)點。研究表明，EGFR信號參與骨原始細(xì)胞群維持，成骨分化平衡調(diào)控［13］，并與前體成骨細(xì)胞活性和增殖緊密相關(guān)［14］。最近的研究也表明，CXCL12/ CXCR4信號參與骨形成和骨吸收的調(diào)控［15］，在BMP.9誘導(dǎo)的間質(zhì)干細(xì)胞成骨分化的早期和中期，CXCL12/ CXCR4信號軸也發(fā)揮重要作用［16］。本研究提示EGFR、CXCL12信號可能與BMP.2信號“交談”，參與調(diào)控BMP.2信號，因此進(jìn)一步探索這些分子在成骨分化平衡調(diào)控中的作用機制，具有重要的意義。同時，趨化因子(C.C模體)配體成員CCL2、CCL7等在網(wǎng)絡(luò)中也處于非常重要交叉節(jié)點，是否參與BMP.2信號調(diào)控成骨分化，筆者正在研究中。

綜上所述，本研究通過對C3H10T1/2成骨分化早期基因表達(dá)譜進(jìn)行分析，找到了一些受rhBMP.2調(diào)控表達(dá)的基因，這將是BMP.2成骨信號研究的重要切入點。同時，通過網(wǎng)絡(luò)生物學(xué)分析發(fā)現(xiàn)，EGFR、CXCL12等可能參與BMP.2信號調(diào)控，控制著間質(zhì)干細(xì)胞骨向分化。

圖2 STRING顯示成骨早期部分差異基因蛋白互作網(wǎng)絡(luò)圖

［1］Chang JT，Nevins JR.GATHER：a systems approach to interpreting genomic signatures［J］.Bioinformatics，2006，22 (23)：2926.2933.

［2］Franceschini A，Szklarczyk D，F(xiàn)rankild S，et al.STRING v9.1：protein.protein interaction networks，with increased coverage and integration［J］.Nucleic Acids Res，2013，41 (Database issue)：808.815.

［3］Karsenty G.Minireview：transcriptional control of osteoblast differentiation［J］.Endocrinology，2001，142(7)：2731.2733.

［4］Lian JB，Stein GS，Javed A，et al.Networks and hubs for the transcriptional control of osteoblastogenesis［J］.Rev Endocr Metab Disord，2006，7(1.2)：1.16.

［5］Zuo C，Huang Y，Bajis R，et al.Osteoblastogenesis regulation signals in bone remodeling［J］.Osteoporos Int，2012，23(6)：1653.1663.

［6］Albers J，Schulze J，Beil FT，et al.Control of bone formation by the serpentine receptor Frizzled.9［J］.J Cell Biol，2011，192(6)：1057.1072.

［7］Harris SE，Guo D，Harris MA，et al.Transcriptional regulation of BMP.2 activated genes in osteoblasts using gene expression microarray analysis：role of Dlx2 and Dlx5 transcription factors［J］.Front Biosci，2003，8：1249.1265.

［8］Tokuzawa Y，Yagi K，Yamashita Y，et al.Id4，a new candidate gene for senile osteoporosis，acts as a molecular switch promoting osteoblast differentiation［J］.PLoS Genet，2010，6(7)：1001019.

［9］Mongiat M，F(xiàn)u J，Oldershaw R，et al.Perlecan protein core interacts with extracellular matrix protein 1(ECM1)，a glycoprotein involved in bone formation and angiogenesis［J］. J Biol Chem，2003，278(19)：17491.17499.

［10］Jeong H，Tombor B，Albert R，et al.The large.scale organizationofmetabolicnetworks［J］.Nature，2000，407(6804)：651.654.

［11］Barabasi AL，Oltvai ZN.Network biology：understanding the cell's functional organization［J］.Nat Rev Genet，2004，5(2)：101.113.

［12］Rual JF，Venkatesan K，Hao T，et al.Towards a proteome. scale map of the human protein.protein interaction network［J］.Nature，2005，437(7062)：1173.1178.

［13］Zhu J，Shimizu E，Zhang X，et al.EGFR signaling suppresses osteoblast differentiation and inhibits expression of master osteoblastic transcription factors Runx2 and Osterix［J］.J Cell Biochem，2011，112(7)：1749.1760.

［14］Chandra A，Lan S，Zhu J，et al.Epidermal growth factor receptor(EGFR)signaling promotes proliferation and survival in osteoprogenitors by increasing early growth response 2(EGR2)expression［J］.J Biol Chem，2013，288 (28)：20488.20498.

［15］Shahnazari M，Chu V，Wronski TJ，et al.CXCL12/CXCR4 signaling in the osteoblast regulates the mesenchymal stem cell and osteoclast lineage populations［J］. FASEB J，2013，27(9)：3505.3513.

［16］Liu C，Weng Y，Yuan T，et al.CXCL12/CXCR4 signal axis plays an important role in mediating bone morphogenetic protein 9.induced osteogenic differentiation of mesenchymal stem cells［J］.Int J Med Sci，2013，10(9)：1181.1192.

Identification and analysis of gene expression profiles for early osteoblastic differentiation on C3H10T1/2 mesenchymal stem cells induced by rhBMP-2

ZUO Changqing1WANG Zonggui2ZHONG Yuechun1LU Hanyun1DAI Zhong1WU Tie1CUI Liao1l.Department of Pharmocology,Guangdong Medical College,Guangdong Province,Dongguan523808,China; 2.Department of Biochemistry,Guangdong Medical College,Guangdong Province,Dongguan523808,China

Objective To investigate the gene expression profiles associated with early osteoblastic differentiation on mesenchymal stem cells,and to provide an experimental basis for the study of gene regulation mechanism for early osteoblastic differentiation.Methods The total RNA was extracted from C3H10T1/2 mesenchymal stem cells between control and rhBMP.2 osteoblastic differentiation group,and then cRNA was labeled with fluorescence labels and hybridized with ArraySTAR mouse genechip.The gene expression profiles were analyzed with bioinformatics software GeneSpring and GATHER.The protein interaction network of differentially expressed genes was constructed and analyzed by online software STRING.Results The molecular function gene ontology（GO）including development and organogenesis and cytokine.cytokine receptor interaction signaling pathway were mainly enriched during C3H10T1/2 early osteogenic differentiation.A total of 42 differentially expressed genes were upregulated and 45 genes were downregulated from rhBMP.2 induced osteoblast 1 and 4 days.Protein network analysis showed that Egfr,Cxcl12 were very important signaling molecules which involved in rhBMP.2 induced osteoblast.Conclusion The screening differentially expressed genes and signaling molecules have an important role for regulation of early osteoblastic differentiation and provide further experimental basis for comprehensive analysis of early osteoblastic differentiation.

Early osteoblastic differentiation;Gene expression profiles;Mesenchymal stem cells;Bioinformatics

R966；R963；R336；R329

A

1673-7210（2014）02（a）-0025-04

2013.11.28本文編輯：程銘)

國家自然科學(xué)基金項目(編號81101357)；廣東省中醫(yī)藥局課題(編號20112145)；廣東省東莞市高等院校科研機構(gòu)科技項目(編號2011108102029)；廣東省湛江市科技攻關(guān)項目(編號2010C3104003)。

左長清(1977.)，男，博士；研究方向：干細(xì)胞成骨成脂分化調(diào)控，生物信息學(xué)。