絨山羊胎兒期毛囊發(fā)生發(fā)育關(guān)鍵circRNA的鑒定及功能預(yù)測(cè)分析

尚方正 馬 榮 戎友俊 王 敏 潘劍鋒 梁麗麗 張燕軍* 李金泉

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué) 動(dòng)物科學(xué)學(xué)院，呼和浩特 010018；2.農(nóng)業(yè)部肉羊遺傳育種重點(diǎn)實(shí)驗(yàn)室 呼和浩特 010018；3.內(nèi)蒙古自治區(qū)動(dòng)物遺傳育種與繁殖重點(diǎn)實(shí)驗(yàn)室，呼和浩特 010018；4.內(nèi)蒙古自治區(qū)山羊遺傳育種工程技術(shù)研究中心，呼和浩特 010018)

circRNA是一類比較特殊的RNA，它沒(méi)有游離的5′帽子結(jié)構(gòu)和3′poly(A)結(jié)構(gòu)，并且對(duì)核酸酶不敏感[1]。circRNA的作用機(jī)制主要包括：調(diào)控親本基因的表達(dá)[2-3],與 RNA 結(jié)合蛋白相互作用[4],翻譯蛋白質(zhì)[5],作為競(jìng)爭(zhēng)性內(nèi)源RNA調(diào)控基因的表達(dá)[6-8]。絨山羊毛囊發(fā)育受多種信號(hào)分子的調(diào)控，特別是近年來(lái)較多的circRNA在絨山羊毛囊發(fā)生發(fā)育有了部分研究進(jìn)展。研究發(fā)現(xiàn)circRNA-1926 可與miR-148a/b-3p 競(jìng)爭(zhēng)性結(jié)合，促進(jìn)CDK19基因的表達(dá)和絨山羊次級(jí)毛囊干細(xì)胞的分化[9-10]；通過(guò)對(duì)成年的遼寧絨山羊和內(nèi)蒙古絨山羊生長(zhǎng)期皮膚樣本轉(zhuǎn)錄組測(cè)序數(shù)據(jù)分析，在遼寧絨山羊中鑒定出17個(gè)上調(diào)的circRNAs和15個(gè)下調(diào)的circRNAs，進(jìn)一步研究發(fā)現(xiàn)ciRNA128、circRNA6854、circRNA4154和circRNA3620在不同毛被遼寧絨山羊中存在顯著差異表達(dá)，揭示circRNA可能對(duì)羊絨細(xì)度有調(diào)控作用[11]。其他方面研究發(fā)現(xiàn)在湖羊皮膚毛囊中共鑒定出5 527個(gè)circRNA，其中在卷毛和直毛比對(duì)組中有114個(gè)差異表達(dá)的circRNA。circRNA在絨山羊胎兒期毛囊發(fā)生發(fā)育中的研究還相對(duì)匱乏。

絨山羊毛囊可分為初級(jí)毛囊和次級(jí)毛囊。初級(jí)毛囊與次級(jí)毛囊形態(tài)發(fā)生的啟動(dòng)處于胎兒發(fā)育的不同時(shí)期，張燕軍等[12-13]研究發(fā)現(xiàn)初級(jí)毛囊的啟動(dòng)早于次級(jí)毛囊。在胎兒期45～55 d，皮膚形成完整的表皮結(jié)構(gòu)，毛囊還未發(fā)生；在55～65 d，胎兒各部位初級(jí)毛囊開始發(fā)生；在65 d時(shí)，體側(cè)部可觀察到明顯的初級(jí)毛囊毛芽；在65～75 d，在胎兒各部位可觀察到次級(jí)毛囊原始體，次級(jí)毛囊開始發(fā)生；在75 d時(shí)，體側(cè)部可觀察到明顯的次級(jí)毛囊毛芽。

研究發(fā)現(xiàn)，定位于細(xì)胞質(zhì)中的外顯子circRNA主要和mRNA競(jìng)爭(zhēng)miRNA的靶標(biāo)結(jié)合位點(diǎn)，從而調(diào)控mRNA的表達(dá)[14]；而內(nèi)含子circRNA主要通過(guò)參與其親本基因的表達(dá)，如circ-ankrd52是來(lái)自于ANKRD52的第二內(nèi)含子，ci-ankrd52的敲低會(huì)導(dǎo)致ANKRD52 mRNA的表達(dá)顯著降低，表明ci-ankrd52 參與了其親本基因的的順式表達(dá)[15]。因此，本研究在課題組前期得到絨山羊胎兒期(45、55、65和75 d)4個(gè)時(shí)期皮膚組織circRNA表達(dá)譜的基礎(chǔ)上[16]，篩選鑒定絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)的外顯子circRNA和內(nèi)含子circRNA，并初步探究其功能，為后續(xù)解析circRNA在絨山羊胎兒期毛囊發(fā)生發(fā)育中的分子機(jī)制提供理論基礎(chǔ)。

1 材料與方法

1.1 樣品采集

于內(nèi)蒙古金萊牧業(yè)科技有限責(zé)任公司選取12只飼喂條件相同的內(nèi)蒙古絨山羊阿爾巴斯型三歲母羊分為四組，每組3只，進(jìn)行同期發(fā)情處理，記錄配種時(shí)間。采取胎兒期45、55、65和75 d各3只羊皮膚樣品。剖腹產(chǎn)后用焦碳酸二乙酯水清洗胎兒。使用無(wú)菌、無(wú)酶的一次性手術(shù)刀和鑷子快速采集體側(cè)部皮膚，放入冷凍儲(chǔ)存管中，編號(hào)，在液氮中快速冷凍，并儲(chǔ)存在-80 ℃冰箱中，用于qRT-PCR 試驗(yàn)。所有的胎兒皮膚標(biāo)本均按照國(guó)際動(dòng)物研究指導(dǎo)原則采集，并獲得內(nèi)蒙古農(nóng)業(yè)大學(xué)科學(xué)研究和學(xué)術(shù)道德委員會(huì)和內(nèi)蒙古農(nóng)業(yè)大學(xué)生物醫(yī)學(xué)研究倫理的批準(zhǔn)(批準(zhǔn)文號(hào)[2020]056)。

1.2 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)circRNA的篩選

基于全轉(zhuǎn)錄組測(cè)序得到的circRNA在絨山羊胎兒期(45、55、65和75 d)4個(gè)時(shí)期表達(dá)譜RNA-Seq數(shù)據(jù)被提交到SRA數(shù)據(jù)庫(kù)，注冊(cè)號(hào)(SRR13306949、SRR13306948、SRR13306947、SRR13306946、SRR13306945、SRR13306944、SRR13306943、SRR13306942、SRR13306941、SRR13306940、SRR13306939、SRR13306938)及各比對(duì)組中差異表達(dá)的circRNA(log2(Fold Change)絕對(duì)值≥1且P<0.05)。因此結(jié)合毛囊在不同時(shí)期的發(fā)生發(fā)育研究結(jié)果，篩選絨山羊毛囊發(fā)生發(fā)育關(guān)鍵circRNA。

1.3 總RNA的提取及反轉(zhuǎn)錄DNA制備

總RNA的提取按照Trizol Reagent試劑操作說(shuō)明進(jìn)行。將提取的總RNA用NanoDrop 2000紫外分光光度計(jì)測(cè)定總RNA的質(zhì)量濃度。利用1%的瓊脂糖凝膠電泳檢測(cè)RNA的完整性。將總RNA等量分成兩份，一份不做處理(供內(nèi)參定量使用)，另一份經(jīng)RNase R(RNA外切酶)37 ℃ 25 min，70 ℃ 10 min處理去除線性RNA。按照Takara的 PrimeScriptRT reagent Kit with gDNA Eraser反轉(zhuǎn)錄試劑盒置于PCR儀中37 ℃ 15 min；85 ℃ 5 s；4 ℃∞獲得 cDNA模板。反轉(zhuǎn)錄后的cDNA放置于-20 ℃冰箱保存?zhèn)溆谩?/p>

1.4 引物設(shè)計(jì)

以circRNA的全長(zhǎng)序列為模板設(shè)計(jì)絨山羊毛囊發(fā)育相關(guān)circRNA的引物，β-actin作為內(nèi)參基因，相關(guān)引物信息詳見表1，引物由廣州銳博生物有限公司合成。

表1 circRNA qRT-PCR引物序列信息Table 1 circRNA qRT-PCR primer sequence information

1.5 circRNA 的實(shí)時(shí)熒光定量PCR驗(yàn)證

qRT-PCR反應(yīng)體系為20 μL，其中TB Green premix EXTapⅡ10 μL，cDNA 2 μL，上下游引物各0.8 μL，ROX 0.4 μL，ddH2O 6 μL。使用LightCycler?96 PCR擴(kuò)增儀進(jìn)行擴(kuò)增，擴(kuò)增程序?yàn)椋侯A(yù)變性：95 ℃ 30 s；PCP反應(yīng)：95 ℃ 10 s，62 ℃ 30 s，72 ℃ 10 s，40個(gè)循環(huán)；熔解：95 ℃ 10 s，60 ℃ 60 s，97 ℃ 1 s；冷卻：37 ℃ 30 s。利用2-ΔΔCt進(jìn)行計(jì)算。

1.6 外顯子circRNA調(diào)控網(wǎng)絡(luò)的構(gòu)建

運(yùn)用TargetScan和miRanda(分析TargetScan和miRanda值的大小，TargetScan閾值為50，該值越大說(shuō)明兩者存在互作的可能性越大；miRanda閾值為-10，該值越小說(shuō)明兩者存在互作的可能性越大)，初步得到circRNA-miRNA靶向關(guān)系，進(jìn)一步對(duì)靶向miRNA的靶基因進(jìn)行預(yù)測(cè)，將富集到毛囊發(fā)育相關(guān)重要通路WNT、TGF-β和Notch等信號(hào)通路中的靶基因構(gòu)成的調(diào)控網(wǎng)絡(luò)作為毛囊發(fā)生發(fā)育相關(guān)circRNA-miRNA-mRNA重要的調(diào)控網(wǎng)絡(luò)。并用Cytoscape將其可視化。

1.7 內(nèi)含子circRNA宿主基因的功能分析

綜合分析22個(gè)差異表達(dá)內(nèi)含子circRNA宿主基因在KEGG的富集情況，探究其富集的通路是否有與絨山羊胎兒期毛囊發(fā)育相關(guān)的通路。

1.8 生物信息學(xué)分析及相關(guān)數(shù)據(jù)分析

運(yùn)用Hisat2和Tophat-fusion軟件進(jìn)行參考基因組比對(duì)分析，運(yùn)用CIRCExplorer和CIRI鑒定識(shí)別circRNA，得到circRNA其序列全長(zhǎng)及來(lái)源基因，進(jìn)一步在NCBI查找相應(yīng)circRNA的結(jié)構(gòu)組成(Ensembl)。最后使用Cytoscape將circRNA-miRNA-mRNA調(diào)控網(wǎng)絡(luò)可視化。

運(yùn)用SPSS 18.0進(jìn)對(duì)相關(guān)數(shù)據(jù)進(jìn)行Kruskal-Wallis秩和檢驗(yàn)及Nemenyi檢驗(yàn)，其中P<0.05表示差異顯著，P<0.01表示差異極顯著，P>0.05表示差異不顯著。

2 結(jié)果與分析

2.1 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)circRNA的篩選

在絨山羊胎兒期45、55、65和75 d皮膚中分別鑒定出2 910、3 710、3 228和3 059個(gè)circRNA。其中在絨山羊胎兒D75組與D45組中上調(diào)59個(gè)，下調(diào)33個(gè)circRNA；胎兒D65組與D45組中上調(diào)96個(gè)，下調(diào)63個(gè)circRNA；胎兒D55組與D45組中，上調(diào)76個(gè)，下調(diào)42個(gè)circRNA。將上述3個(gè)比對(duì)組中均差異表達(dá)的circRNA作為毛囊發(fā)生發(fā)育相關(guān)的circRNA(圖1)。共篩選到4個(gè)絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)的circRNA，包括circRNA2049、circRNA2225、circRNA3411和circRNA5681。

圖1 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)circRNA的篩選Fig.1 Screening of circRNAs related to hair follicle development in cashmere goats during fetal period

2.2 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)circRNA的鑒定

circRNA2049長(zhǎng)度為536 bp，由山羊ATRNL1基因第2～5外顯子的全部序列組成，為外顯子circRNA；circRNA2225長(zhǎng)度為791 bp，由山羊SERINC5基因第3～8外顯子的全部序列組成，為外顯子circRNA；circRNA5681序列長(zhǎng)度為265 bp，由山羊基因102 184 877第3～5外顯子的全部序列組成，為外顯子circRNA；circRNA3411長(zhǎng)度為333 bp，由山羊102 173 234基因第175～178外顯子的全部序列組成，為外顯子circRNA(表2)。

表2 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)circRNA的結(jié)構(gòu)組成Table 2 The structural composition of circRNA related to hair follicle development in the fetal period of cashmere goats

表2(續(xù))

2.3 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)circRNA差異表達(dá)分析

本研究運(yùn)用qRT-PCR檢測(cè)了circRNA2049、circRNA2225、circRNA3411和circRNA5681在絨山羊胎兒期(45、55、65和75 d)皮膚組織中的表達(dá)量。Kruskal-Wallis秩和檢驗(yàn)表明上述circRNA在4個(gè)時(shí)期中均差異表達(dá)(P<0.05)，Nemenyi組間多重比較發(fā)現(xiàn)circRNA2049、circRNA2225和circRNA5681在D65vsD55差異顯著(P<0.05)(圖2)。表明circRNA2049、circRNA2225和circRNA5681可能調(diào)控初級(jí)毛囊的發(fā)生；circRNA3411在D75vsD45差異顯著(P<0.05)，推測(cè)circRNA3411可能調(diào)控次級(jí)毛囊的發(fā)生。

*P<0.05表示差異顯著。* P<0.05 indicates significant difference.圖2 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)circRNA qRT-PCR驗(yàn)證分析Fig.2 Validation analysis of circRNA qRT-PCR related to hair follicle development in the fetal period of cashmere goats

2.4 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)circRNA靶向miRNA的預(yù)測(cè)

本研究運(yùn)用TargetScan和miRanda預(yù)測(cè)了篩選到的4個(gè)外顯子circRNA靶向的miRNA(TargetScan score percentile>50且miRanda energy<-10)。其中circRNA2049有19個(gè)靶向的miRNA，circRNA2225有8個(gè)靶向miRNA，circRNA3411有4個(gè)靶向miRNA，circRNA5681有1個(gè)靶向miRNA(表3)。

表3 TargetScan和miRanda 預(yù)測(cè)circRNA靶向的miRNATable 3 Prediction of miRNAs targeted circRNA by Targetscan and miRanda

表3(續(xù))

2.5 絨山羊胎兒期毛囊發(fā)生發(fā)育外顯子circRNA調(diào)控網(wǎng)路構(gòu)建

本研究基于上述circRNA靶向miRNA的預(yù)測(cè)結(jié)果，同樣利用TargetScan和miRanda預(yù)測(cè)靶向miRNA的靶基因。并結(jié)合mRNA的富集情況，篩選富集在Notch signaling pathway 和TGF-β signaling pathway與毛囊發(fā)生發(fā)育相關(guān)的重要信號(hào)通路中的mRNA構(gòu)建了調(diào)控絨山羊毛囊發(fā)生發(fā)育的調(diào)控網(wǎng)絡(luò)。其中DLL4位于NOTCH信號(hào)通路的起始段，同時(shí)研究表明DLL4與羊毛囊發(fā)育相關(guān)。因此circRNA2049-chi-miR-27b-3p-DLL4，circRNA2225-chi-miR-145-5p-DLL4和circRNA9181-chi-miR-145-5p-DLL4成為后續(xù)絨山羊胎兒期毛囊發(fā)生發(fā)育的重要調(diào)控網(wǎng)絡(luò)(圖3)。

綠色三角形為circRNA;粉色四邊形為miRNA;藍(lán)色圓形為mRNA。The green triangle was circRNA; The pink quadrilateral was miRNA; Blue circle was mRNA.圖3 絨山羊毛囊發(fā)生發(fā)育重要circRNA-miRNA-mRNA調(diào)控網(wǎng)絡(luò)Fig.3 The important circRNA-miRNA-mRNA regulatory network of hair follicle development in cashmere goats

2.6 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)內(nèi)含子circRNA宿主基因功能分析

本研究為探究?jī)?nèi)含子circRNA宿主基因的功能，分析了6個(gè)對(duì)照組中22個(gè)差異表達(dá)的ciRNA(表4)。其中ciRNA450的宿主基因100 860 901富集在Wnt信號(hào)通路中；ciRNA587的宿主基因102 189 934富集在Mapk信號(hào)通路中(表5)，而Wnt 信號(hào)通路和Mapk信號(hào)通路是毛囊發(fā)生發(fā)育相關(guān)的重要通路，推測(cè)ciRNA可能通過(guò)調(diào)控其宿主基因的表達(dá)來(lái)影響毛囊的發(fā)生發(fā)育。

表4 絨山羊胎兒期毛囊發(fā)生發(fā)育差異表達(dá)的內(nèi)含子circRNA統(tǒng)計(jì)Table 4 Statistics of intron circRNA differentially expressed during hair follicle development in the fetal period of cashmere goats

表5 絨山羊胎兒期毛囊發(fā)生發(fā)育內(nèi)含子circRNA宿主基因富集的KEGG通路Table 5 KEGG pathway enriched in circRNA host genes in the development of fetal hair follicles in cashmere goats

3 討 論

3.1 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)circRNA結(jié)構(gòu)

毛囊的發(fā)生發(fā)育受到多種信號(hào)通路分子的調(diào)控，除了眾多編碼RNA外，非編碼也是調(diào)控毛囊發(fā)生發(fā)育的重要因子，包括miRNA、lncRNA和circRNA，但circRNA在毛囊的發(fā)生發(fā)育中調(diào)控作用的研究還相對(duì)匱乏。因此，本研究根據(jù)課題組前期研究結(jié)果，篩選出4個(gè)絨山羊毛囊發(fā)生發(fā)育相關(guān)的circRNA，進(jìn)一步運(yùn)用生物信息學(xué)分析出這4個(gè)長(zhǎng)度分別為536、791、265和333 bp，這與Li等[17]在綿羊circRNA測(cè)序中鑒定到circRNA長(zhǎng)度在200～400 bp的結(jié)果基本一致。同時(shí)，篩選的4個(gè)circRNA均由不同個(gè)數(shù)的外顯子組成，均為外顯子circRNA。與Li等[18]在circRNA的鑒定中發(fā)現(xiàn)外顯子circRNA占80%以上的研究結(jié)果相符。qRT-PCR表明篩選的4個(gè)circRNA均為差異表達(dá)的circRNA，組間差異分析發(fā)現(xiàn)部分結(jié)果與測(cè)序結(jié)果不一致，但整體表達(dá)趨勢(shì)一致，推測(cè)可能與實(shí)驗(yàn)誤差有關(guān)。

3.2 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)內(nèi)含子circRNA功能預(yù)測(cè)

全部由內(nèi)含子環(huán)化而成的circRNA稱內(nèi)含子circRNA，這些ciRNAs通常定位于細(xì)胞核中，而它們的親本基因主要定位在細(xì)胞質(zhì)中。研究表明，環(huán)狀RNA可以通過(guò)與線性轉(zhuǎn)錄本競(jìng)爭(zhēng)性剪接，從而影響宿主基因的表達(dá)[4]。Zhang等[19]研究發(fā)現(xiàn)內(nèi)含子環(huán)狀RNA能順式調(diào)控宿主基因的表達(dá)。本研究對(duì)前期測(cè)序結(jié)果中22個(gè)差異表達(dá)的內(nèi)含子環(huán)狀RNA的宿主基因富集情況進(jìn)行整合，發(fā)現(xiàn)部分宿主基因富集在Wnt和Mapk信號(hào)通路中，其中100 860 901基因富集于Wnt信號(hào)通路中。Wu等[20]篩選了Wnt信號(hào)通路中Wnt10A，驗(yàn)證發(fā)現(xiàn)對(duì)絨山羊胎兒期毛囊發(fā)生發(fā)育有調(diào)控作用，馬森等[21]研究發(fā)現(xiàn)，Wnt/β-catenin 通路與絨山羊絨毛周期性再生相關(guān)，并且能夠促進(jìn)絨毛的快速生長(zhǎng)，另外發(fā)現(xiàn)Wnt蛋白是胎兒毛囊發(fā)育啟動(dòng)的最早信號(hào)之一[22]；102 189 934基因是本研究中富集到MAPK通路中的基因。研究顯示，Tβ4能夠影響β-catenin的表達(dá)，通過(guò)誘導(dǎo)VEGF及MMP-2的表達(dá)，使得MAPK中P38的激活，進(jìn)而影響了毛囊分布和毛干數(shù)量[23]；另外研究發(fā)現(xiàn)，調(diào)控毛囊周期的一些蛋白受MAPK信號(hào)激活，腫瘤壞死因子(TNF)作為MAPK信號(hào)中的重要成員之一，對(duì)毛囊的細(xì)胞凋亡有重要作用[24]。因此，推測(cè)內(nèi)含子circRNA可能通過(guò)調(diào)控宿主基因的表達(dá)來(lái)影響毛囊的發(fā)生發(fā)育。

3.3 絨山羊胎兒期毛囊發(fā)生發(fā)育相關(guān)外顯子circRNA功能預(yù)測(cè)

全部由外顯子環(huán)化而成的circRNA為外顯子circRNA。外顯子circRNA主要定位于細(xì)胞質(zhì)中，近年來(lái)，有大量的研究發(fā)現(xiàn)定位于細(xì)胞質(zhì)中的circRNA可以和mRNA競(jìng)爭(zhēng)miRNA的靶標(biāo)結(jié)合位點(diǎn)，從而調(diào)控mRNA的表達(dá)[25]。研究表明Notch 信號(hào)和毛囊的形態(tài)發(fā)生有直接的關(guān)系[26]。本研究篩選到了絨山羊胎兒期毛囊發(fā)生發(fā)育關(guān)鍵基因Notch2和DLL4，而先前研究中也發(fā)現(xiàn)Notch2對(duì)調(diào)節(jié)羊絨細(xì)度具有潛在的重要意義[27]，DLL4與羊毛囊發(fā)育相關(guān)[28]。TGF-β信號(hào)通路同樣是毛囊發(fā)育相關(guān)的重要通路，本研究中篩選到TGF-β信號(hào)通路中的FBLN5對(duì)毛發(fā)生長(zhǎng)有影響[29]，本課題組Han等[30]研究表明chi-miR-199a-5p通過(guò)調(diào)控TGF-β信號(hào)通路中的TGF-β2來(lái)調(diào)控次級(jí)毛囊的發(fā)生發(fā)育通路，另外研究發(fā)現(xiàn)TGF-β1對(duì)小鼠毛囊退行期有調(diào)控作用，敲除TGF-β1小鼠生長(zhǎng)期向退行期的轉(zhuǎn)變明顯延遲[31]，TGF-β2可以抑制毛干的延長(zhǎng)并誘導(dǎo)毛囊退行期的形態(tài)學(xué)改變[32]，TGF-β2可能通過(guò)調(diào)控內(nèi)源性的細(xì)胞凋亡引起毛囊角質(zhì)形成細(xì)胞的凋亡，進(jìn)而抑制毛發(fā)的生長(zhǎng)[33]。本研究預(yù)測(cè)關(guān)鍵circRNA靶向miRNA以及靶向miRNA的靶基因，同時(shí)部分靶基因富集在Notch和TGF-β信號(hào)通路中。當(dāng)然，circRNA對(duì)Notch和TGF-β信號(hào)通路調(diào)控的研究也相對(duì)較多，但主要集中在人類疾病的研究上。研究發(fā)現(xiàn)，下調(diào)的環(huán)狀RNA hsa_circ0067301通過(guò)miR-141/Notch信號(hào)通路調(diào)節(jié)子宮內(nèi)膜異位癥上皮-間質(zhì)轉(zhuǎn)化[34]；CircNFIX可通過(guò)Notch信號(hào)通路介導(dǎo)miR-34a-5p，調(diào)控NOTCH1和Notch信號(hào)通路，促進(jìn)膠質(zhì)瘤的發(fā)生[35]；CircAPLP2通過(guò)靶向miR-101-3p激活Notch信號(hào)通路調(diào)控結(jié)直腸癌的增殖和轉(zhuǎn)移[36]。circRNA010567通過(guò)靶向TGF-β1抑制miR-141促進(jìn)心肌纖維化[37]；circRNAcESRP1通過(guò)分泌miR-93-5p 抑制TGF-β通路在小細(xì)胞肺癌化療敏感性中起著關(guān)鍵作用，提示cESRP1可能是一個(gè)有價(jià)值的預(yù)后生物標(biāo)志物和潛在的小細(xì)胞肺癌治療靶點(diǎn)[38]；hsa_circ0001085可能通過(guò)hsa-miR-196b-5p間接調(diào)節(jié)TGF-β信號(hào)通路從而對(duì)前列腺癌細(xì)胞起調(diào)節(jié)作用[39]。因此，認(rèn)為外顯子circRNA主要通過(guò)circRNA-miRNA-mRNA途徑調(diào)控毛囊的發(fā)生發(fā)育。

4 結(jié) 論

本研究篩選了絨山羊毛囊發(fā)生發(fā)育關(guān)鍵circRNA并發(fā)現(xiàn)在山羊胎兒期毛囊發(fā)育各階段存在差異；外顯子circRNA主要通過(guò)競(jìng)爭(zhēng)性結(jié)合miRNA以circRNA-miRNA-mRNA調(diào)控網(wǎng)絡(luò)來(lái)調(diào)控毛囊的發(fā)生發(fā)育；22個(gè)內(nèi)含子circRNA的宿主基因可能通過(guò)Wnt和Mapk等信號(hào)通路調(diào)控毛囊的發(fā)生發(fā)育。本研究為解析circRNA在絨山羊毛囊發(fā)生發(fā)育中的分子機(jī)制奠定了基礎(chǔ)，也為絨山羊分子輔助育種提供了理論依據(jù)。