英文名:Long Non-coding RNA Derived from lncRNA– mRNA Co-expression Networks Modulates the Locust Phase Change

期刊：Genomics Proteomics Bioinformatics

IF：7.051（1區）

通訊作者單位：中科院、河北大學

研究背景

長鏈非編碼RNA（lncRNA）調節基因表達、動物行為等各種生物學過程。盡管蛋白質編碼基因、microRNA和神經肽在亞洲飛蝗的表型可塑性調節中起著重要作用，但有關lncRNA在此過程中功能研究較少。本文應用高通量RNA-seq來比較蝗蟲型變時程中lncRNA和mRNA的表達模式。結果顯示lncRNA在型變的早期階段反應更快。功能注釋表明，早期改變的lncRNA在分散和群居階段采用了不同的途徑來應對種群密度的變化。篩選了分散和群居階段的網絡中兩個重疊的中樞lncRNA基因座進行功能驗證。本文進一步證明LNC1010057為潛在的蝗蟲型變因子。這項工作為深入了解蝗蟲型變的分子機制并擴大lncRNA在動物行為中的作用范圍提供了重要的數據。

材料方法

實驗材料：散居化處理為將群居型蝗蟲單獨飼養0h、4h、8h和16h后取出手機蝗蟲大腦進行測序；

群居化處理為將10只散居型蝗蟲與20只群居型蝗蟲飼養于一個小籠子（10厘米×10厘米×10厘米）中，于0h、4h、8h和16h后取出，收集蝗蟲大腦進行測序。在同一時間點采集樣品的三個生物學重復，提取RNA，建立cDNA文庫并進行RNA-seq。對lncRNA和mRNA進行差異表達分析及qRT-PCR驗證，對不同時間點的樣本進行STEM(ShortTime-seriesExpressionMiner)分析，并構建lncRNA–mRNA共表達網絡的構建。對目標中心節點lncRNA進行RNAi，采集其行為錄像數據進行行為數據分析。

實驗結果

1、蝗蟲lncRNA的外顯子比mRNA少但更長

研究發現大約78％的lncRNA包含2個外顯子，而mRNA中包含的外顯子數量為1至120（圖1C）。因此，lncRNA的外顯子明顯長于mRNA的外顯子（平均長度1142bpvs.263bp，圖1D，左）。同時，lncRNA的內含子明顯短于mRNA的內含子（平均長度：10086bpvs.12442bp，圖1D，右）。表達水平分析表明，lncRNA的總體表達水平明顯低于mRNA的表達水平（平均值為0.6vs.1.9；圖1E）。但是，表達特異性分析表明，lncRNA的表達受時間限制的程度要高于mRNA（平均值為0.672對0.436圖1F）。基于相對于mRNA的lncRNA基因組位置，蝗蟲lncRNA被分類為基因間、重疊區、有義內含子、反義內含子、有義外顯子和反義外顯子。蝗蟲77％以上的lncRNA是長基因間的ncRNA（lincRNA，圖1G）。這些結果表明，蝗蟲lncRNA與mRNA在結構和表達上有很大不同。蝗蟲lncRNAs更長，具有更少但更長的外顯子和更短的內含子。此外，lncRNA的表達模式顯示出比mRNA更高的時間特異性。

圖1lncRNA和mRNA之間的不同結構和表達

2、群居型蝗蟲中特異性表達的lncRNA比在散居型蝗蟲中表達的更多

在散居型和群居型蝗蟲大腦中共表達9722個lncRNA（73.9％），而散居型蝗蟲中特異性表達962個lncRNA（7.3％），群居型蝗蟲中特定表達2469個lncRNA（18.8％）（圖2A，頂部）。分別在散居型和群居型蝗蟲中特異性表達了479和1090個mRNA（3.1％和7.1％）（圖2A，底部）。這些結果表明，與mRNA相比，在兩個蝗蟲相型特異性表達的lncRNA的百分比更高，而在群居蝗蟲中比散居型蝗蟲中表達的lncRNAs更多。與散居型蝗蟲相比，群居型蝗蟲中335個lncRNA和779個mRNA的表達水平下調，而313個lncRNA和261個mRNA的表達上調[倍數變化（FC）>2和P<0.05；圖2B]。lincRNA在下調和上調的lncRNA中所占的比例高（分別為81.2％和87.8％），其次分別是反義外顯子和有義內含子lncRNA（圖2B）。在表達中按FC排名的前10個lncRNA基因顯示在圖2C中。這表明散居型和群居型蝗蟲大腦中表達的lncRNA的數量及其表達水平明顯不同。

圖2lncRNA在蝗蟲型變中顯示出不同的表達變化模式

3、lncRNA對種群密度變化的快速應答

為了進一步分析lncRNA和mRNA的表達模式，使用了短時間序列表達挖掘器（ShortTime-seriesExpressionMiner，STEM）對基于表達的轉錄本進行聚類。在CS期間，將214個lncRNA和242個mRNA分別分為8個和9個重要的表達譜（圖3D）。其中lncRNA表達譜15以及mRNA表達譜15、13和16直到群居化處理后8或16小時才改變。這些配置文件稱為后期更改的配置文件（模式d）。與后期更改的配置文件相反，早期更改的配置文件以4小時時間點開始的轉錄表達變化為特征。根據4小時后的表達變化，將早期變化的輪廓細分為模式a（早期變化）、模式b（早期中變化）和模式c（可持續變化）。在IG過程中，269個lncRNA和639個mRNA分別聚集成6個和7個顯著表達譜。所有lncRNA配置文件均為早期變化（圖3E）。在mRNA表達譜中，有6個是早期改變的，而譜12是晚期改變的。在CS和IG期間，lncRNA的聚類分布圖的數量與mRNA的分布無明顯差異。但是，通過計算譜圖中的轉錄本數量，發現早期改變的lncRNA的百分比比CS中的mRNA的百分比高61.1％（88.3％vs.54.8％，圖3F）。同樣，IG中早期改變的lncRNA的百分比高于IG中的mRNA（100％vs.81.3％，圖3F）。因此，lncRNA比mRNA對散居化和群居化處理的反應更快。在CS中，早期改變的lncRNA的表達變化速率快于285個mRNA的表達變化速率（0.55vs.0.24）和IG（0.77vs.0.39）（圖3G）。因此，早期改變的lncRNA的比例和表達變化率高于早期改變的mRNA。因此，蝗蟲lncRNA比mRNA對種群密度的變化更敏感。

圖3lncRNA對散居化和群居化處理的快速反應

4、lncRNA參與CS和IG早期變化的不同途徑

在CS中，自噬調節、剪接體和肌醇磷酸代謝途徑在上位和至少兩個模塊中均過代表（圖4C）。因此，CS中早期改變的lncRNA可能在調節自噬、RNA剪接和信號轉導途徑中發揮重要作用。lncRNA對自噬的調控參與了群居化的初期和中期，而對剪接體和磷酸肌醇代謝的調控則參與了整個過程。與神經遞質釋放有關的突觸小泡循環途徑的lncRNA僅在模式a（早期改變）模塊中被過代表（圖4C）。該結果表明，早期改變的lncRNA可能在CS期間調節神經系統中具有重要作用。此外，一些lncRNA與已知的型變相關基因密切相關。例如，早期改變的lncRNALNC531328.2、LNC1425451.2、LNC1088763.7和LNC492755.1與基因NPYR、NPF1a和Vat1相關（圖4A）。持續變化的lncRNALNC494161.1和LNC1065048.14與多巴胺途徑中的基因Ebony和Vat1相關。在CS網絡中，程度值前5％的lncRNA被視為hublncRNA（圖4A）。計算了lncRNA基因座的程度值，在CS中，基于度數排名前5％的10個lncRNA基因座被鑒定為hublncRNA。其中，四個基因表達上調，其他六個基因表達下調（圖4D）。在IG期間，包括谷氨酸能突觸、多巴胺能突觸和膽堿能突觸途徑在內的與突觸有關的途徑得以豐富。這些途徑中的大多數都富含早期改變的模塊。同時，多巴胺能突觸和膽堿能突觸途徑僅參與模式a（早期改變）的模塊（圖4E）。IG中也豐富了信號轉導途徑，例如鈣信號傳導、Ras信號傳導、胰島素信號傳導和MAPK信號傳導途徑。功能注釋的結果表明，與CS相比，IG中早期變化的lncRNA參與的突觸相關和信號處理途徑更多。

圖4早期改變的lncRNA參與CS和IG的不同途徑

5、LNC1010057可能調節蝗蟲型變

為為了驗證LNC1010057和LNC992414在蝗蟲型變中的功能，進行了一系列分子生物學實驗。首先，克隆了在LNC1010057和LNC992414基因座中鑒定出的長的轉錄本。LNC1010057由重復的A、B和C元素組成，這些元素依次分布并重復4次，但C元素重復三次（圖5A）。序列比對證明LNC992414在5’端與LNC1010057共享相似的重復序列，但是它們是從不同的基因組基因座轉錄而來的（圖5A）。盡管LNC992414的定量表達水平可以通過特異性引物檢測，但LNC1010057的表達水平為重復元件的總表達水平。如預期的那樣，LNC1010057和LNC992414的表達模式極為相似。在CS期間，兩種lncRNA的表達持續增加，并且在16h幾乎增加了四倍（圖5B）。在IG期間，4h后表達水平顯著下降，之后保持相對穩定（圖5B）。LNC1010057和LNC992414之間的序列結構和表達模式的相似性表明它們是同源的lncRNA。但是，實時qPCR分析表明，在散居型和群居型蝗蟲的大腦中，LNC1010057的表達水平比LNC992414的表達水平高約1000倍（圖5C）。其次，為了測試LNC1010057和LNC992414是否參與蝗蟲型變調節，在通過RNAi抑制蝗蟲大腦中的表達后進行了行為分析。顯著降低LNC1010057的表達水平（7470對3764610）和LNC992414（1.9vs.1.0）（圖5D）后行為分析表明，群居型蝗蟲其行為顯著改變為散居狀態（圖5E）。此外，多個與型有關的行為參數發生了變化。總移動距離（TDM）和總移動持續時間（TDMV）顯著降低（圖5F），但是，移動速度沒有區別。同時群居蝗蟲的特定喜好行為顯著下降（58.3對-13.5圖5F）。但LNC992414表達降低并未引起從群居狀態到散居狀態的轉變（圖5G和H）。與對照相比，與型相關的行為參數（包括運動和特定物種的優選行為）沒有差異（圖5I）。LNC992414的RNAi實驗不會引起行為變化，因此排除了LNC992414調節蝗蟲型變的可能性。這些結果表明是LNC1010057而不是LNC992414潛在地調節了蝗蟲的型變。

圖5LNC1010057潛在地調節了蝗蟲的型變

總結

lncRNA被證實是生物過程的關鍵調控因子，調節包括mRNA轉錄、穩定性、翻譯和翻譯后修飾等多種生物途徑。之前研究表明昆蟲lncRNA參與了例如殺蟲劑抗性、繁殖力和腺體凋亡等生物過程，但是尚未有證據證明lncRNA可以調節非模型昆蟲的行為。蝗蟲是世界范圍內的一種農業害蟲，表現出顯著的表型可塑性。為了鑒定與其型變相關的lncRNA，本文系統地分析了蝗蟲lncRNA的表達并注釋其功能，證實與mRNAs相比lncRNAs顯示對型變更敏感的響應，并證實了其中一個lncRNA可調節型變相關行為。本研究揭示了lncRNA在蝗蟲型變中的重要作用以及表型中蛋白編碼基因和lncRNA之間的相互作用。深入解析蝗蟲散居型和群居型轉變的分子機制，為可持續治理蝗蟲的新策略和新方法的開發提供了基礎。

研究背景

細胞自噬是一種非常保守的細胞內降解系統，近期的研究表明小鼠中細胞自噬現象表現出了生物鐘節律，nr1d1基因缺陷能導致小鼠骨骼肌中自噬活動的紊亂。然而生物鐘調節細胞自噬的機理還不清楚。斑馬魚的細胞自噬活動是否表現出生物鐘節律還沒有報道，本研究旨在以斑馬魚為實驗材料研究nr1d1基因在調節細胞自噬的生物鐘節律方面的作用。

實驗設計

實驗材料：2種材料：野生型WT和nr1d1突變型斑馬魚
2個時間點：CT2和CT14（分別為受精后122和134小時）;
無生物學重復
測序方法：百邁客Illumina HiSeq2000 PE125

結果展示

1.表型描述及關鍵基因在幼魚和肝臟中的表達情況

首先用TEM（透射電鏡）對一天內不同時間（ZT0，ZT6，ZT12，ZT18）的肝臟切片進行觀察，統計一天內自噬體數量的變化（紅色箭頭處，圖1，A），并進行統計（圖1，B），研究表明斑馬魚肝臟內自噬體的數量符合生物鐘節律。

TEM觀測及自噬體數量統計結果

選取了幾個以前研究過的與自噬體形成、移動、延伸、融合、降解有關的基因進行RT-PCR分析，如maplc3b，bnip3，becn1等。材料為受精后72-96小時的斑馬魚幼體和肝臟材料。由圖可見這些基因都表現出明顯的生物鐘節律。圖中白色底色表示白天，灰色底色表示黑夜。

圖2 關鍵基因的RT-PCR檢測

2.用TALEN技術構建nrid1突變斑馬魚，并進行表型研究和關鍵基因表達分析

用TEM觀察WT和nr1d1突變斑馬魚肝臟中的自噬體數量，研究表明自噬體數量在nr1d1中明顯高于WT。

圖3 不同樣品斑馬魚肝臟自噬體觀測及統計

用RT-PCR研究在WT和nr1d1突變斑馬魚中關鍵基因的表達差異，nr1d1突變體中基因表達模式與WT中明顯不同。

圖4 差異基因的RT-PCR驗證

3.用ChIP和熒光素酶研究Nr1d1直接調控的下游基因

體外熒光素酶實驗表明Nr1d1蛋白能直接結合ulk1a基因的promoter區域，體內ChIP實驗進一步證明Nr1d1蛋白能夠結合ulk1a基因的promoter區域。說明Nr1d1能直接調控ulk1a基因的表達。

圖5 熒光素和chip結果

4. 野生型WT和nr1d1突變型斑馬魚轉錄組分析

為了進一步說明nr1d1基因的功能，文章選取了野生型WT和nr1d1突變型斑馬魚進行轉錄組研究，并且取了2個時間點（CT2和CT14）。研究表明在CT2階段兩者有975個差異表達基因，在CT14階段有1360個差異表達基因。與WT相比，在nr1d1突變體中上調的基因數量要大于下調的基因數量。

同時，選取部分關鍵的差異表達基因進行RT-PCR驗證，結果與轉錄組分析結果一致。

圖6 差異基因韋恩圖

圖7 差異基因聚類圖

圖8 關鍵基因的RT-PCR驗證

研究結論

通過研究WT野生型和nr1d1突變型斑馬魚的基因表達情況，表明Nr1d1基因在生物鐘調節和細胞自噬過程中起著重要作用。生物鐘相關基因直接受Nr1d1調控。此外，通過轉錄組分析表明Nr1d1還調控脅迫、代謝、信號傳導、翻譯后修飾等生理過程。

圖9 Nr1d1在生物鐘和細胞自噬過程中的作用模式圖

亮點

1.首次利用斑馬魚研究生物鐘和細胞自噬的關系；
2.發現Nr1d1在生物鐘調節和細胞自噬過程中的關鍵作用；
3.轉錄組研究進一步證明Nr1d1的作用，并揭示Nr1d1其它生理功能。
4.雖然轉錄組的內容在這篇文章中比重不大，但卻揭示了除了生物鐘調節和細胞自噬，Nr1d1基因還有很多其它的功能，顯著提升了文章的水平，并為后續研究打好了基礎。

研究背景

小球藻為單細胞可食用綠藻，能夠合成蝦青素和油脂；蝦青素是一種結構獨特的酮式類胡蘿卜素，是自然界已知最強的抗氧化劑。含油脂高的藻類可作為生物質能源。

實驗設計

材料選擇：小球藻 取種子細胞培養到生長對數期的后期，分成3組，用于做不同處理；
處理方式：對照組T1：不做處理；
高光處理組T2：光照90μmol光子?2 s?1，6h；
葡萄糖誘導組T3：30g/L葡萄糖處理， 6h；
測序方法：Hiseq2500，4G/樣；

技術路線：

研究結果

1. 蝦青素和脂類定量分析

對照和葡萄糖處理組相比，葡萄糖誘導下的小球藻中，蝦青素含量96h后增加了30倍，而處理組中96h后只增加了8倍。脂類和TFA含量也有類似的情況。對照組中這三種成分增加較少可能是培養基中營養成分限制的原因，在氮饑餓的情況下，導致在相同器官中蝦青素和脂肪酸產量增加。

2. RNA-SEQ分析

三種不同的處理構建三個文庫，測序獲得clean data 10.63G，Trinty組裝獲得32931個unigenes。52.8%的unigenes在NR,Swiss-Prot,KEGG,GO,COG獲得注釋。

3. 差異表達基因分析

計算RPKM，做差異表達分析。T1與T2差異表達基因數為580，T1與T2差異表達基因數為1489。在光處理和葡萄糖處理兩種情況下，差異基因功能注釋表明這些基因都是與翻譯、核糖體、蛋白過程、脅迫和碳固定相關的。在光處理下，卟啉和葉綠素代謝、光刺激應答是被富集的。在葡萄糖處理情況想，葡萄糖、能量產生和轉化、脂肪合成等都是被富集的；而這些基因在光合作用、核糖體、蛋白運輸通路中是減少的。在碳源充足的情況下，生長和能量存儲基因是上調表達的，而光合作用是下調表達的。

4. RT-PCR驗證

作者選了10個與蝦青素和三酰甘油合成有關的unigenes進行了RT-PCR驗證，其中9個基因的表達趨勢與RNA-seq結果一致。只有PDS兩種結果不同，可能是RNA-seq結果不準確，因為前人的研究成果中，葡萄糖是上調PDS表達的，與本文中RT-PCR結果一致。

5. 蝦青素合成

在轉錄組測序結果中，發現MEP通路中的8個基因都被鑒定到，但是在MVA通路中只有2個基因被檢測到，其中HMG表達不顯著，ACAT表達量下降。這些結果表明小球藻中蝦青素的合成通路中IPP和DMAPP是MEP通路合成的。其中，發現一個可能編碼β-胡蘿卜素激酶（BKT2）與數據庫中的BKT1比對率僅為56%。功能分析發現BKT2可以轉化β-胡蘿卜素成為角黃素。下圖為蝦青素合成通路圖。在葡萄糖處理情況下，與蝦青素合成相關的6個基因上調差異表達。在小球藻中，葡萄糖誘導蝦青素的合成，伴隨著BKT1，CHYb,PDS在24h內上調表達。為研究轉錄組水平基因表達量與產物含量的關系，分析了5個限速基因的表達量，這些基因在葡萄糖處理72h或92h表達量達到最高，相應的，蝦青素的積累量也在96h達到穩定水平。

6. 三酰甘油的合成

藻類和植物中脂肪酸的生物合成主要在葉綠體中進行。乙酰-CoA作為合成的起始物質。基于轉錄組測序，小球藻質體中脂肪酸合成基因被鑒定到。在葡萄糖處理情況下，這些基因大都表達量上調。另外這次測序結果中還發現幾個基因對應的不同轉錄本。由于BC和SAD是脂肪酸合成的限速酶基因，做了進一步分析，發現BC和SAD表達量從6h到72h表達量持續上升，對應的三酰甘油含量從24h到96h持續增加。

三酰甘油（TAG）合成有兩條通路：乙酰-CoA依賴的通路和乙酰-CoA不依賴的通路；轉錄組數據中兩條通路中的基因都被鑒定到。

文章亮點

1：材料選擇+處理條件；
小球藻關注的成分具有重要的經濟價值，處理條件是可以促進經濟成分積累，研究更加有利用價值；

2：實驗數據+前人結果，深入驗證；
測序結果中重點關注關鍵成分（蝦青素和TAG）合成通路中的每一個基因，對其進行差異分析，并做不同時期的表達量分析；解釋清楚在處理情況下，物質合成與基因表達的關系。針對本文的研究結果與前人研究成果比較分析。

在瘤胃中的短鏈脂肪酸（SCFAs）對反芻動物的生長和健康起著關鍵的作用，微生物G蛋白偶聯受體（GPR）和微生物脫乙酰化酶（HDAC）可能是這些影響的主要調節通路。該研究主要是選用不同比例非纖維碳水化合物攝入（15-30%）的山羊模型，通過轉錄組測序和16srRNA測序聯合研究瘤胃上皮細胞和瘤胃微生物協同的響應機制。通過研究發現，瘤胃中微生物來源的SCFAs對上皮細胞的生長和代謝受到GPR和HDAC調節網絡的影響，通過對這些調控機制和飲食組成的相互關系的理解，可更深入的了解瘤胃的新陳代謝機制，更好的促進牧業的發展。

樣品取材：共6只雄性山羊，隨機分為兩組，第一組飼養為65%的干草+35%的飼料（MC組），另一組為90%的干草+10%的飼料（LC組）；喂養28d后進行取樣。

微生物多樣性樣品：飼養第28天，在清晨喂養0h、2h、5h以及8h進行取樣，瘤胃內容物用4層紗布取樣 ，收集瘤胃液15ml，-20℃保存用于提取。

轉錄調控測序取樣：取10cm2瘤胃組織，用冷卻的PBS進行清洗，去除肌肉層，取上皮組織切塊于TRIzol保存液保存。液氮速凍5min后，-80保存用于RNA提取。

短鏈脂肪酸（SCFAs）濃度測定：上述樣品加入5%HgCl2用于濃度測定。

測序分析：

瘤胃內容物微生物多樣性測序：細菌V3+V4區、Illumina Miseq 平臺測序

瘤胃上皮細胞轉錄組測序測序： PE125 測序、Illumina Hiseq2500 平臺測序

飼養不同處理組（MC和LC）、喂養取樣不同時間點間瘤胃內總短鏈脂肪酸（SCFA）、PH以及短鏈脂肪酸（SCFA）主要成分的變化情況進行分析，如下圖。

?圖1.1?總SCFA變化情況

圖1.2?PH變化情況

圖1.3?SCFA主成分變化情況

2、瘤胃微生物多樣性分析：

2.1細菌群落結構分析

在細菌門水平，一共有14個原核的門被鑒定，主要是厚壁菌門（35.7-30.1%），擬桿菌門（26.6-43.6%）、互養菌門（11-7%），且通過門水平維恩圖分析發現，在兩組（MC和LC）間門水平組成相同，比較分析兩組間微生物變化情況發現（如圖：FIG2），MC組相比于LC組，互養菌門增長了57%，無壁菌門增長了330%。而黏膠球形菌門降低了63%，纖維桿菌門降低了65%；

圖2.1?兩組間門水平Venn圖分析

圖2.2 兩組間門水平OTU分析

在細菌屬水平。共鑒定出75個屬水平細菌，在兩組間共有70個屬，在MC組有三個特有的屬，在LC水平有2個特有的屬，普氏菌屬（10.4-17.9%）在兩組間都為主要屬水平的菌。

圖3.1 兩組間屬水平Venn圖分析

圖3.2?兩組間屬水平OTU分析

2.2 微生物群落的多樣性和豐富性

在門水平，MC組的擬桿菌門和厚壁菌門顯著性高于Lc組（P＜0.05），互養菌門顯著性低于Lc組（如下圖4），通過最大似然法（ML）計算27個豐度大于1% 的OTU進行分析顯示，在MC組顯著增加的OTU屬于子單胞菌科，疣微菌科、韋榮球菌科、疣微菌門，相反的，58%（7/12）顯著降低的OTU屬于普雷沃氏菌科（如下圖5）；

?圖4

圖5

3.瘤胃表達譜分析：

3.1 瘤胃上皮細胞GPRs和HDACs的表達譜分析

RNA-Seq測序方法用于研究山羊瘤胃微生物G蛋白偶聯受體（GPRs）和微生物脫乙酰化酶（HDACs）表達情況以及調控網絡，過濾得到127M clean data，平均有83%數據比對到NCBI山羊參考基因組，得到73個GPR家族成員和11個HDAC家族成員比對到山羊基因組，轉錄組數據分析發現有20個GPR家族成員和7個HDAC家族成員在瘤胃上皮細胞中表達。通過比較LC組中基因的表達發現GPR1，87，89A，155顯著上調（P＜0.05），在MC組中GPR107，游離脂肪酸受體4（FFAR4, also known as GPR120），羥基烴酸受體2（HCAR2, also known as GPR109A）顯著上調（P＜0.05），通過組內基因表達分析，在兩組內都發現GPR家族中GPR87表達*高，HDAC家族中HDCA1的表達*高

3.2 GPRs和HDACs保守序列分析

為了研究GPR和HDAC在進化過程中差異表達的保守序列，通過用所有基因進行ML進化樹分析，所有的 GPRS和HDACS在脊椎動物門都是高度保守的。研究GPR家族內成員的保守序列發現，在GPR家族成員沒有超過30%的相似性，在HDAC家族成員之間也發現了同樣的結果。

圖6

3.3 差異表達的GPRS和HDACs基因相關KEGG通路分析

為了研究GPR和HDAC共表達網絡的生物學意義，改研究注釋了相關的信號通路和KEGG功能，發現顯著差異的鄰近基因被分類到瘤胃上皮細胞的生理調節過程，因此發現了兩類和上皮細胞生長調節有關系的功能網絡，其中一類是上皮細胞生長網絡（包括細胞凋亡，增殖和分化，見下圖7），另外一類是上皮細胞的代謝網絡（包括輔酶因子和維生素代謝，能量代謝、氨基酸代謝等，見下圖8）；

圖7：上皮細胞生長網絡圖

圖8：上皮細胞代謝網絡圖

在上皮細胞生長網絡中，基因GPR1，89和155的表達上調導致LAMTOR3蛋白，蛋白激酶和ELK1表達上調，這三種蛋白都存在于MAPK信號通路上，和細胞的增殖分化有關。此外，GPR1的 表達增加和酸性酰胺酶（ASAH1）的下調相關，ASAH1存在于調節細胞凋亡、生存和增殖的信號通路上。對于HDACs家族基因HDAC4,5,6和10基因的上調和五個調節細胞凋亡，生存和增殖基因的下調有關。

在上皮細胞代謝網絡中，GPR1、87和89A基因的上調表達與輔酶2、4以及4L（ND2，4，4L），ATP5H，NDUFA4有關，HDAC1的下調表達ND6的上調表達有關。以上所有相關的酶都和能量代謝的氧化磷酸化有關。

4.轉錄調控和16s rRNA數據聯合分析?

GPR和HDAC的表達、SCFAS主成分濃度比例以及屬水平細菌相對豐度關系分析

通過CCA相關性分析發現，顯著上調的GPRs和HDACs與8個顯著增加的微生物屬（56個 OTUS）是呈正相關，顯著下調的GPRs和HDACs和9個顯著減少的細菌屬呈正相關（63個OTU），然而HDAC1與丁酸鹽的比例呈負相關。

該研究得到的結果揭示，SCFA調節的GPR和HDAC共調控網絡存在于瘤胃的上皮細胞，該共調節網絡作用可作用于動物敏感√確地接受微生物區的信號調節，這些調控網絡調節上皮細胞各種生理學過程，尤其是細胞的生長和新陳代謝，在促進動物的生長和維持上皮細胞的完整性起到關鍵的作用。此外，這些調控網絡重要的調控機制對于共生的細菌來說，主要通過調節上皮細胞生理過程提高細菌在宿主中的共生條件。通過了解宿主動物和微生物區之間互相的調節機制，加上相應飲食中的外界干預調節因素，可以可持續性的更好的提高動物的健康和生長。