張 倩 王俊奎 羅 芳 和東遷 陳麗堯 盧童童 吳少飛 張瑞雪 陶金忠

(寧夏大學(xué)農(nóng)學(xué)院，銀川 750021)

剩余采食量(RFI)是衡量反芻動(dòng)物飼料轉(zhuǎn)化率的主要指標(biāo)，是指畜禽實(shí)際采食量和預(yù)期采食量的差值，可以用來(lái)有效評(píng)價(jià)動(dòng)物飼料轉(zhuǎn)化率、群體代謝水平以及選育標(biāo)準(zhǔn)[1]。關(guān)于RFI的相關(guān)研究涉及到諸多方面，如干物質(zhì)采食量(DMI)、飼料消化率、機(jī)體合成與分解代謝、運(yùn)動(dòng)量以及體溫調(diào)節(jié)等[2]。Patience[3]發(fā)現(xiàn)育成豬的每日干物質(zhì)采食量中34%的能量用于機(jī)體維持需要。從節(jié)約飼養(yǎng)成本的角度來(lái)看，降低動(dòng)物的維持需要，能夠提高飼料的轉(zhuǎn)化效率。Arthur等[4]和Richardson等[5]用超聲波掃描技術(shù)測(cè)定了安格斯牛皮下脂肪厚度，結(jié)果發(fā)現(xiàn)低RFI的個(gè)體體脂含量更低，體蛋白質(zhì)含量更髙，且背膘厚度與RFI顯著相關(guān)。RFI的高低反映了動(dòng)物的飼料利用率，即RFI越高的個(gè)體，其飼料轉(zhuǎn)化效率越低、代謝越快、采食次數(shù)越多、活動(dòng)量越大[6-7]。能量代謝的途徑是從葡萄糖開(kāi)始，由己糖激酶和葡萄糖-6-磷酸酶介導(dǎo)生成葡萄糖-6-磷酸[8]，再由磷酸己糖異構(gòu)酶的催化形成果糖-6-磷酸，之后生成丙酮酸進(jìn)入三羧酸循環(huán)(TCA)供能。飼料轉(zhuǎn)化效率和體內(nèi)能量代謝是否相關(guān)，目前還不清楚。

反芻動(dòng)物的RFI是根據(jù)個(gè)體體重和生長(zhǎng)[平均日增重(ADG)]的需求估算的。RFI在表型上與體重和ADG無(wú)關(guān)，但與采食量有關(guān)，因此對(duì)RFI的表型測(cè)定是一項(xiàng)亟待解決的問(wèn)題，因?yàn)樗枰獪y(cè)量每個(gè)個(gè)體2～3個(gè)月的采食量，需要投入巨大的時(shí)間、金錢、物力和人力，并且受到記錄方法和設(shè)備儀器等的限制[9]。

代謝組學(xué)是一種運(yùn)用高通量、高靈敏度的現(xiàn)代分析技術(shù)[10]，是以相對(duì)分子質(zhì)量在1 000以下的小分子代謝物為研究對(duì)象，包括內(nèi)源性和外源性分子，可以反映個(gè)體、器官組織及細(xì)胞代謝所產(chǎn)生的小分子代謝圖譜，能更為直接地反映當(dāng)前代謝的復(fù)雜網(wǎng)絡(luò)，將代謝途徑同生物機(jī)制相聯(lián)系，從而對(duì)疾病的診斷、發(fā)病機(jī)制的研究以及疾病的治療和預(yù)防提供新的技術(shù)和手段[11]。代謝組學(xué)研究揭示的不僅僅是遺傳學(xué)或基因組學(xué)提供的信息，更揭示了動(dòng)物遺傳學(xué)和生理表型之間的關(guān)系[12]。Clemmons等[13]研究表明，泛酸鹽、谷氨酰胺、肉堿和半胱氨酸等差異代謝物與黑安格斯的營(yíng)養(yǎng)代謝密切相關(guān)，可以以泛酸鹽、谷氨酰胺、肉堿和半胱氨酸等差異代謝物作為生物標(biāo)志物來(lái)解釋黑安格斯的RFI。運(yùn)用代謝組學(xué)研究反芻動(dòng)物RFI，不但可以提高反芻動(dòng)物的飼料轉(zhuǎn)化效率，還可以識(shí)別血漿中的潛在生物標(biāo)志物，但目前在羊上沒(méi)有看到相關(guān)的研究。因此，本試驗(yàn)擬采用液相色譜-串聯(lián)質(zhì)譜(LC-MS/MS)技術(shù)分析不同RFI下灘羊血漿中能量代謝相關(guān)代謝物的變化，篩選出對(duì)灘羊RFI有影響作用的差異代謝物，以探究在不同RFI下灘羊能量代謝情況，并為使用差異代謝物作為標(biāo)志性物質(zhì)預(yù)測(cè)RFI提供參考。

1 材料與方法

1.1 試驗(yàn)動(dòng)物

于寧夏某灘羊保種場(chǎng)選用345只[公羊163只、母羊182只(在飼養(yǎng)過(guò)程中會(huì)淘汰生病、體重不合格的羊只)]具有詳細(xì)系譜信息、出生日期相近、生長(zhǎng)狀況良好、體重[(30±5)kg]相近的6月齡灘羊，統(tǒng)一飼養(yǎng)管理，每天早晚各投喂1次顆粒料，控制采食時(shí)間范圍，自由飲水。試驗(yàn)分為過(guò)渡期15 d，預(yù)試期10 d，正試期50 d。準(zhǔn)確記錄每只羊每天的剩料量，并于正試期第1、10、20、30、40、50天晨飼前空腹稱重，定期對(duì)圈舍進(jìn)行清掃消毒。

1.2 試驗(yàn)分組

飼養(yǎng)階段結(jié)束后，利用Koch等[14]提出的預(yù)測(cè)RFI模型計(jì)算個(gè)體全期RFI值(第1～50天的RFI)，按照群體RFI的平均值(mean)和標(biāo)準(zhǔn)差(SD)將試驗(yàn)羊群分成極低RFI組(RFImean+0.5SD，16只)，采集血液，其中采樣灘羊的RFI情況如表1所示。然后按照性別和RFI高低將所采集的灘羊血液樣品分為3組(極低RFI母羊vs極高RFI母羊組、極低RFI公羊vs極高RFI公羊組和極低RFI羊vs極高RFI羊組)。

表1 采樣灘羊的RFI數(shù)據(jù)Table 1 RFI data of sampling Tan sheep

1.3 樣品采集

飼養(yǎng)試驗(yàn)結(jié)束后，對(duì)所選的極低和極高RFI灘羊進(jìn)行頸部采血，將10 mL血液置于抗凝管(含肝素鈉)中。對(duì)采集后血液樣品進(jìn)行3 000 r/min離心10 min，將上層血漿分裝，保存于-80 ℃冰箱。

1.4 儀器與試劑

表2 所需試劑Table 2 Required reagents

1.5 樣品預(yù)處理

因部分樣品發(fā)生溶血現(xiàn)象，故共檢測(cè)到26份血漿樣品，其中公羊、母羊各檢測(cè)13份。將血漿樣品從-80 ℃冰箱中取出，在冰上解凍后，渦旋10 s混勻，4 ℃、3 000 r/min離心5 min；取上清液50 μL加入到對(duì)應(yīng)的已編號(hào)的離心管中，加入150 μL純甲醇內(nèi)標(biāo)提取液；渦旋5 min后4 ℃、12 000 r/min離心10 min；取140 μL上清液裝入對(duì)應(yīng)進(jìn)樣瓶?jī)?nèi)襯管中，用于LC-MS/MS分析。

1.6 LC-MS/MS分析

電噴霧離子源(electrospray ionization，ESI)溫度450 ℃，質(zhì)譜電壓5 500 V(正)、-4 500 V(負(fù))，離子源氣體Ⅰ(GSⅠ)40 psi，氣體Ⅱ(GSⅡ)55 psi，氣簾氣(curtain gas，CUR)35 psi，碰撞誘導(dǎo)電離(collision-activated dissociation，CAD)參數(shù)設(shè)置為中。在三重四極桿(QTRAP)中，每個(gè)離子對(duì)是根據(jù)優(yōu)化的去簇電壓(declustering potential，DP)和碰撞能(collision energy，CE)進(jìn)行掃描檢測(cè)[15]。

1.7 數(shù)據(jù)處理與分析

將質(zhì)譜檢測(cè)的原始數(shù)據(jù)在R Studio軟件平臺(tái)，采用XCMS軟件對(duì)數(shù)據(jù)進(jìn)行特征峰提取，峰對(duì)齊，保留時(shí)間(retention time，RT)校正以及數(shù)據(jù)過(guò)濾，提取峰面積。將經(jīng)XCMS軟件提取的數(shù)據(jù)進(jìn)行完整性檢查，對(duì)組內(nèi)缺失值進(jìn)行刪除或補(bǔ)充，刪除極值及組內(nèi)缺失值超過(guò)50%的代謝物離子峰，數(shù)據(jù)以內(nèi)標(biāo)校正后進(jìn)行總峰面積歸一化處理，以確保各樣品間和代謝物之間的可平行比較。對(duì)數(shù)據(jù)進(jìn)行總峰面積歸一化后，使用SIMCA-P 14.1軟件中對(duì)數(shù)據(jù)進(jìn)行pareto-scaling處理后，進(jìn)行多元統(tǒng)計(jì)分析，包括主成分分析(principal component analysis,PCA)和正交偏最小二乘判別分析(orthogonal partial least squares discrimination analysis，OPLS-DA)；同時(shí)，在建模過(guò)程中對(duì)模型數(shù)據(jù)進(jìn)行置換檢驗(yàn)(permutation test)并計(jì)算變量投影重要度(variable important for the projection，VIP)。在此基礎(chǔ)上對(duì)數(shù)據(jù)進(jìn)行單變量統(tǒng)計(jì)分析，包括Student’st檢驗(yàn)和變異倍數(shù)(fold change，F(xiàn)C)分析，根據(jù)各代謝物P值和FC值篩選差異代謝物。

2 結(jié)果與分析

2.1 灘羊血漿樣品質(zhì)控(QC)分析

將QC樣品在正、負(fù)離子檢測(cè)模式下的質(zhì)譜總離子流(TIC)圖進(jìn)行譜圖疊加比較(圖1)，結(jié)果表明，代謝物檢測(cè)總離子流的曲線重疊性高，即保留時(shí)間和峰強(qiáng)度均一致，表明質(zhì)譜對(duì)同一樣品不同時(shí)間檢測(cè)時(shí)，信號(hào)穩(wěn)定性較好。整個(gè)試驗(yàn)過(guò)程中儀器誤差引起的變異較小，數(shù)據(jù)質(zhì)量可靠。

圖1 QC樣品在正離子(A)和負(fù)離子模式(B)下TIC重疊圖譜Fig.1 TIC overlapping patterns of QC samples in positive (A)and negative ion modes (B)

2.2 多元統(tǒng)計(jì)分析

PCA模型參數(shù)R2X、Q2分別表示模型對(duì)X變量的解釋率和模型預(yù)測(cè)能力。通過(guò)對(duì)不同組樣品數(shù)據(jù)進(jìn)行PCA(圖2)可知，極低RFI母羊vs極高RFI母羊組、極低RFI公羊vs極高RFI公羊組、極低RFI羊vs極高RFI羊組的模型解釋率(R2X)分別為0.235、0.351、0.145，聚類較好，有一定的分離趨勢(shì)，說(shuō)明代謝物檢測(cè)系統(tǒng)的穩(wěn)定性較好。

ELF與EHF分別表示極低RFI母羊與極高RFI母羊，ELM與EHM分別表示極低RFI公羊與極高RFI公羊，EL與EH分別表示極低RFI羊與極高RFI羊。下圖同。圖A、B、C分別為極低RFI母羊vs極高RFI母羊組、極低RFI公羊vs極高RFI公羊組和極低RFI羊vs極高RFI羊組PCA得分圖。ELF and EHF denoted very low RFI ewes and very high RFI ewes,respectively;ELM and EHM denoted very low RFI rams and very high RFI rams,respectively;EL and EH denoted very low RFI sheep and very high RFI sheep,respectively.The same as below.Figures A,B and C represented PCA score plots of very low RFI ewes vs very high RFI ewes group,very low RFI rams vs very high RFI rams group and very low RFI sheep vs very high RFI sheep group,respectively.圖2 3組灘羊PCA得分圖Fig.2 PCA score plots of Tan sheep in three groups

為了提高模型的有效性，直接區(qū)分組間變異進(jìn)行了OPLS-DA，結(jié)果見(jiàn)圖3。在OPLS-DA得分圖下，各組樣品區(qū)分良好。由置換檢驗(yàn)?zāi)Ｐ?圖4)可知，極低RFI母羊vs極高RFI母羊組Q2截距=-0.432(Q2截距<0)，極低RFI公羊vs極高RFI公羊組Q2截距=-0.346(Q2截距<0)，極低RFI羊vs極高RFI羊組Q2截距=-0.516(Q2截距<0)，表明OPLS-DA模型未過(guò)度擬合，說(shuō)明各樣品間存在差異。

圖A、B、C分別為極低RFI母羊vs極高RFI母羊組、極低RFI公羊vs極高RFI公羊組和極低RFI羊vs極高RFI羊組OPLS-DA得分圖。Figures A,B and C represented OPLS-DA score plots of very low RFI ewes vs very high RFI ewes group,very low RFI rams vs very high RFI rams group and very low RFI sheep vs very high RFI sheep group,respectively.圖3 3組灘羊OPLS-DA得分圖Fig.3 OPLS-DA score plots of Tan sheep in three groups

R2 intercept和Q2 intercept表示R2與Q2回歸直線與Y軸的截距，Q2 intercept<0表明OPLS-DA模型未過(guò)度擬合。橫坐標(biāo)代表隨機(jī)分組的Y與原始分組的Y的相關(guān)性，縱坐標(biāo)代表R2與Q2得分。圖A、B、C分別為極低RFI母羊vs極高RFI母羊組、極低RFI公羊vs極高RFI公羊組和極低RFI羊vs極高RFI羊組置換檢驗(yàn)分析圖。R2 intercept and Q2 intercept represented the Y-intercept of R2 and Q2 regression lines.Q2 intercept<0 indicated that the OPLS-DA model was not overfitting.The abaxial axis represented the correlation between the random group Y and the original group Y,and the vertical axis represented the scores of R2 and Q2.Figures A,B and C represented replacement test analysis diagrams of very low RFI ewes vs very high RFI ewes group,very low RFI rams vs very high RFI rams group and very low RFI sheep vs very high RFI sheep group,respectively.圖4 3組灘羊OPLS-DA置換檢驗(yàn)Fig.4 OPLS-DA replacement test of Tan sheep in three groups

2.3 差異代謝物篩選

本試驗(yàn)結(jié)合多元統(tǒng)計(jì)分析，利用OPLS-DA計(jì)算所得VIP，以VIP>1、P<0.05及FC>1.30或<0.77為標(biāo)準(zhǔn)篩選差異代謝物，并對(duì)其進(jìn)行鑒定，結(jié)果如表3所示。本試驗(yàn)在極低RFI公羊vs極高RFI公羊組中共篩選鑒定出2種差異代謝物，為D-木酮糖-5-磷酸和L-鳥(niǎo)氨酸鹽酸鹽，均下調(diào)；在極低RFI母羊vs極高RFI母羊組中共篩選鑒定出4種差異代謝物，其中3種上調(diào)，上調(diào)物質(zhì)分別為腺嘌呤、檸檬酸、L-天冬氨酸，1種下調(diào)，下調(diào)物質(zhì)為L(zhǎng)-瓜氨酸；在極低RFI羊vs極高RFI羊組中共篩選鑒定出4種差異代謝物，其中2種上調(diào)，上調(diào)物質(zhì)分別為腺嘌呤、檸檬酸，2種下調(diào)，下調(diào)物質(zhì)分別為L(zhǎng)-蘇氨酸、D-木酮糖-5-磷酸。

表3 已鑒定差異代謝物Table 3 Identified differential metabolites

2.4 差異代謝物聚類分析

將篩選出的差異代謝物利用Metabo Analyst 5.0在線分析平臺(tái)進(jìn)行層次聚類，結(jié)果如圖5所示。結(jié)果表明，多數(shù)樣品能很好地聚集在一起，但個(gè)別樣品仍然存在分離情況，表明組內(nèi)有差異，但具有相似變化的代謝物可以聚集在同一族內(nèi)。因此，進(jìn)一步采用受試者工作特征(ROC)曲線對(duì)篩選出的差異代謝物進(jìn)行進(jìn)一步分析，篩選出ROC曲線下面積(AUC)>0.90的差異代謝物。圖6顯示，極低RFI母羊vs極高RFI母羊組中檸檬酸AUC達(dá)到0.90以上，P<0.05，說(shuō)明檸檬酸是識(shí)別母羊組的主要指標(biāo)，但不能作為識(shí)別灘羊整體的主要指標(biāo)。利用液相色譜-四級(jí)桿-飛行時(shí)間質(zhì)譜(LC-Q/TOF-MS)峰強(qiáng)箱式圖對(duì)極低RFI母羊vs極高RFI母羊組的檸檬酸含量進(jìn)行進(jìn)一步分析，結(jié)果如圖7所示，檸檬酸能明顯分開(kāi)極高RFI母羊和極低RFI母羊。

圖A、B、C分別為極低RFI母羊vs極高RFI母羊組、極低RFI公羊vs極高RFI公羊組和極低RFI羊vs極高RFI羊組聚類熱圖。Figures A,B and C represented clustering heat maps of very low RFI ewes vs very high RFI ewes group,very low RFI rams vs very high RFI rams group,and very low RFI sheep vs very high RFI sheep group,respectively.圖5 3組灘羊差異代謝物聚類熱圖Fig.5 Clustering heat maps of differential metabolites of Tan sheep in three groups

圖A、B依次為極低RFI羊vs極高RFI羊組和極低RFI母羊vs極高RFI母羊組ROC曲線。AUC：ROC曲線下面積。The ROC curves of very low RFI sheep vs very high RFI sheep group and very low RFI ewes vs very high RFI ewes group were shown successively in figures A and B.AUC:area under curve.圖6 極低RFI羊vs極高RFI羊組和極低RFI母羊vs極高RFI母羊組ROC曲線Fig.6 ROC curves of very low RFI sheep vs very high RFI sheep group and very low RFI ewes vs very high RFI ewes group

圖7 極低RFI母羊vs極高RFI母羊組峰強(qiáng)度箱式圖Fig.7 Peak intensity boxplot of very low RFI ewe vs very high RFI ewe group

3 討 論

Elolimy等[16]研究結(jié)果表明性別的差異會(huì)影響瘤胃上皮基因溶質(zhì)載體家族9成員A1(SLC9A1)、缺氧誘導(dǎo)因子1亞單位A(HIF1A)和烏頭酸酶2(ACO2)的表達(dá)，并且公牛相較于母牛來(lái)說(shuō)mRNA的表達(dá)更高，而羥基丁酸脫氫酶1(BDH1)、溶質(zhì)載體家族9成員A2(SLC9A2)和丙酮酸脫氫酶E1亞單位A1(PDHA1)的表達(dá)更低；盧果等[17]通過(guò)代謝組學(xué)的方法研究36組人尿液樣品，研究表明由于性別的不同從而分離出不同的生物標(biāo)志物，女性尿液中檸檬酸的含量顯著高于男性；陳建軍[18]研究結(jié)果表明抑郁癥診斷由于性別差異標(biāo)志物也有所差異，甘醇酸(glycolate)、次黃嘌呤(hypoxanthine)可作為女性潛在的生物標(biāo)志物，酪氨酸(tyrosine)、檸檬酸鹽(citrate)和琥珀酸鹽(succinate)可作為男性潛在的生物標(biāo)志物，為篩選抑郁癥性別差異性診斷標(biāo)示物提供了進(jìn)一步研究數(shù)據(jù)；更有研究表明性別差異會(huì)影響藥物代謝[19]。本試驗(yàn)利用靶向代謝組學(xué)的方法對(duì)3組灘羊進(jìn)行血漿代謝物分析，在極低RFI母羊vs極高RFI母羊組中，共篩選出4種差異代謝物，分別為腺嘌呤、檸檬酸、L-天冬氨酸和L-瓜氨酸；而在極低RFI公羊vs極高RFI公羊組中，共篩選出2種差異代謝物，分別為L(zhǎng)-鳥(niǎo)氨酸鹽酸鹽和D-木酮糖-5-磷酸，結(jié)果顯示，公、母羊血漿代謝物種類不同是由于性別因素造成的。因此，性別的差異可能會(huì)導(dǎo)致血漿代謝物的差異。

在能量代謝過(guò)程中檸檬酸是三羧酸循環(huán)中的重要物質(zhì)。Kim等[20]研究發(fā)現(xiàn)，檸檬酸可顯著降低高脂飲食小鼠肝臟中脂肪的含量。Cheng等[21]研究發(fā)現(xiàn)，人類連續(xù)7 d服用500 mg/d劑量的檸檬酸，可顯著提高骨骼肌中糖原合成率。此外，檸檬酸還可用于脂肪和膽固醇合成，有利于物質(zhì)消化和營(yíng)養(yǎng)吸收，維持機(jī)體生長(zhǎng)發(fā)育[22]。本研究發(fā)現(xiàn)極低RFI母羊vs極高RFI母羊組和極低RFI羊vs極高RFI羊組篩選的差異代謝物中檸檬酸均上調(diào)，即檸檬酸表達(dá)量增加，進(jìn)入能量代謝途徑的相關(guān)物質(zhì)含量也相應(yīng)升高，能量轉(zhuǎn)化的效率也不斷增加，釋放更多的能量。低RFI灘羊檸檬酸表達(dá)量升高，因此低RFI灘羊產(chǎn)生的能量相較于高RFI灘羊更多。這說(shuō)明低RFI動(dòng)物在生長(zhǎng)潛力方面優(yōu)于高RFI動(dòng)物。在本研究中極低RFI母羊vs極高RFI母羊組檸檬酸的AUC大于0.9，說(shuō)明檸檬酸是識(shí)別母羊組的主要指標(biāo)，有望成為母灘羊RFI的生物標(biāo)志物。

腺嘌呤核苷酸(三羧酸循環(huán)過(guò)程中的輔助因子，不在代謝網(wǎng)絡(luò)圖里表現(xiàn))是機(jī)體重要的能量物質(zhì)[23]，L-天冬氨酸是生成草酰乙酸的前體物質(zhì)。腺嘌呤和L-天冬氨酸在極低RFI母羊vs極高RFI母羊組中均上調(diào)。腺嘌呤表達(dá)量增加，意味著整個(gè)能量代謝的進(jìn)程加快，釋放更多的能量；L-天冬氨酸表達(dá)量增加，能夠提高整個(gè)代謝途徑中能量的轉(zhuǎn)化效率，影響著整個(gè)能量代謝進(jìn)程。

瓜氨酸同組氨酸、脯氨酸、鳥(niǎo)氨酸和精氨酸一樣都是能量通路中轉(zhuǎn)化為谷氨酸的前體物質(zhì)，本試驗(yàn)極低RFI母羊vs極高RFI母羊組中L-瓜氨酸下調(diào)，雖然對(duì)谷氨酸的合成代謝有影響作用，但無(wú)法直接影響整個(gè)能量代謝過(guò)程。因腺嘌呤、L-天冬氨酸和瓜氨酸的AUC均小于檸檬酸，因此不能作為影響母羊組的生物標(biāo)志物。

D-木酮糖-5-磷酸是磷酸戊糖代謝途徑中的代謝物之一，通過(guò)磷酸戊糖途徑生成葡萄糖進(jìn)入三羧酸循環(huán)。本研究發(fā)現(xiàn)在極低RFI公羊vs極高RFI公羊組中D-木酮糖-5-磷酸下調(diào)。D-木酮糖-5-磷酸下調(diào)是由于核酮糖-5-磷酸在生成D-木酮糖-5-磷酸的同時(shí)生成了核糖-5-磷酸(圖8)，D-木酮糖-5-磷酸進(jìn)入磷酸戊糖途徑的含量減少，但核糖-5-磷酸是用于合成核苷酸的主要前體物，雖然D-木酮糖-5-磷酸表達(dá)量的下降并不能直接影響到能量轉(zhuǎn)化效率，但對(duì)磷酸戊糖途徑有一定的影響。L-鳥(niǎo)氨酸鹽酸鹽下調(diào)，與上述母羊L-瓜氨酸下調(diào)表述一致。但由于D-木酮糖-5-磷酸和L-鳥(niǎo)氨酸鹽酸的AUC均較低，故不能用作公羊組的生物標(biāo)志物。

圖8 差異代謝物代謝網(wǎng)絡(luò)圖Fig.8 Metabolic network map of differential metabolites

4 結(jié) 論

本研究發(fā)現(xiàn)不同RFI公與母灘羊的血漿代謝物種類不同，可能是由于性別因素造成的，因此性別也會(huì)影響灘羊能量代謝。在極低RFI母羊vs極高RFI母羊組和極低RFI羊vs極高RFI羊組中，極低RFI母羊和極低RFI羊的檸檬酸均上調(diào)，而檸檬酸作為三羧酸循環(huán)的關(guān)鍵代謝物，可成為影響母灘羊RFI的重要生物標(biāo)志物。