段修軍 孫國波 張蕾 陸艷鳳 王麗華 田培余 劉嘉晟 楊孟廣

摘要： 為挖掘黑羽番鴨肉質(zhì)風味相關候選基因，通過對黑羽番鴨胸肌、腿肌中揮發(fā)性風味物質(zhì)進行測定，并通過Illumina HiSeq2500高通量測序進行轉(zhuǎn)錄組對比分析，結合參考基因組對所獲得的序列進行序列比對、基因注釋和差異表達等分析，篩選出差異表達基因并進行GO富集分析。通過熒光定量 PCR（Real-time PCR，qRT-PCR）方法檢測4個差異候選基因表達水平，驗證測序結果的可靠性。結果表明，黑羽番鴨胸肌中揮發(fā)性風味物質(zhì)辛醛、2，3-辛二酮含量顯著高于腿肌，3-羥基-2-丁酮含量顯著低于腿肌（ P < 0.05）。通過參考基因組比對和差異表達分析，初步獲得614個差異表達基因，171個差異基因在黑羽番鴨腿肌組織上調(diào)表達，443個差異基因在黑羽番鴨腿肌組織下調(diào)表達。結合GO分析和KEGG富集分析，最終獲得了20個候選功能基因，它們分別參與機體氨基酸形成、糖代謝以及脂肪代謝等生物過程，其中10個差異候選基因在黑羽番鴨胸肌中參與肌內(nèi)脂肪代謝過程，這些基因可能通過形成肌內(nèi)揮發(fā)性風味物質(zhì)進而影響肌肉風味。qRT-PCR 驗證結果表明，篩選出的差異候選基因表達趨勢與轉(zhuǎn)錄組測序中表達的趨勢相似，說明測序結果可靠。在黑羽番鴨胸肌、腿肌中揮發(fā)性風味物質(zhì)對比基礎上，結合RNA-Seq技術對黑羽番鴨肉質(zhì)風味前體物質(zhì)特異性形成機制進行綜合分析，篩選出ACOT11、APOA1、FABP3等脂肪代謝相關的差異基因可能與黑羽番鴨胸、腿肌肉質(zhì)風味形成相關。

關鍵詞： 黑羽番鴨; 揮發(fā)性風味物質(zhì); 脂肪代謝; 轉(zhuǎn)錄組測序

中圖分類號： S834+.89?? 文獻標識碼： A?? 文章編號： 1000-4440（2022）03-0739-09

Identification of candidate genes related to meat flavor in black Muscovy duck based on RNA-Seq

DUAN Xiu-jun， SUN Guo-bo， ZHANG Lei， LU Yan-feng， WANG Li-hua， TIAN Pei-yu， LIU Jia-sheng， ?YANG Meng-guang

（Jiangsu Agri-animal Husbandry Vocational College， Taizhou 225300， China）

Abstract： To explore the candidate genes related to meat flavor in black Muscovy duck， volatile flavor compounds in breast and leg muscles of black Muscovy duck were determined， and transcriptome analysis was performed by Illumina HiSeq2500 high-throughput sequencing. Based on reference genome， the obtained sequences were analyzed by sequence alignment， gene annotation and differential expression analysis， and differentially expressed genes were screened out and GO enrichment analysis was performed. Moreover， to verify the sequencing results， qRT-PCR was used to detect the expression levels of four differential candidate genes. Results showed that the contents of octanal and 2，3-octanedione in breast muscle were significantly higher than those in leg muscle， and the content of 3-hydroxy-2-butanone in breast muscle was significantly lower than that in leg muscle （ P < 0.05）. Through genome comparison and differential expression analysis， 614 differentially expressed genes were preliminarily obtained， 171 differentially expressed genes were up-regulated and 443 differentially expressed genes were down-regulated in leg muscles of black Muscovy duck. By GO and KEGG enrichment analysis， 20 functional candidate genes were finally obtained. These candidate genes were involved in biological processes， such as amino acid formation， glucose metabolism and fat metabolism. Among them， ten differentially expressed candidate genes were involved in fat metabolism in breast muscle of black Muscovy duck， and these genes may affect muscle flavor by forming intramuscular volatile flavor compounds. Results of qRT-PCR verification showed that the expression trends of screened differentially expressed candidate genes were similar to those in transcriptome sequencing， indicating that the sequencing results were reliable. Based on the comparation of volatile flavor compounds in breast and leg muscles of black Muscovy duck， the specific formation mechanism of meat flavor precursors in black Muscovy duck was comprehensively analyzed， differentially expressed genes related to fat metabolism such as ACOT11， APOA1， FABP3 may be related to the formation of meat flavor in breast and leg muscles of black Muscovy duck.

Key words： black Muscovy duck; volatile flavor compounds; fat metabolism; transcriptome sequencing

番鴨（ Cairna moschata ）發(fā)源于南美洲和中美洲的熱帶，因其鴨嘴的基部和眼圈周圍具有黑色或紅色的肉瘤被稱為瘤頭鴨，又因其鴨肉瘦肉含量高，肉質(zhì)鮮香，具有野禽風味的特點被稱為麝香鴨 [1] 。番鴨按照羽色可以分為白羽、黑羽和花羽3種類群，其中黑羽番鴨因羽色和肉質(zhì)特異性深受消費者喜愛，但是其體型較小，生長速度較慢 [2-3] 。為進一步提高黑羽番鴨的生產(chǎn)性能，江蘇農(nóng)牧科技職業(yè)學院自2000年開始，利用法國番鴨和福建黑羽番鴨進行選育，構建了羽色、脛、爪全黑的優(yōu)質(zhì)黑羽番鴨類群，其具有生長速度快、體型適中、脂肪含量低、肉質(zhì)鮮美等優(yōu)點 [4-5] 。

孫國波等 [3] 對公、母黑羽番鴨不同階段胸、腿肌肉質(zhì)進行對比分析，結果表明黑羽番鴨膽固醇含量較低，氨基酸豐富且含量高，公母鴨胸、腿肌營養(yǎng)成分存在一定差異性。王錦鋒等 [6] 研究結果表明，黑羽番鴨胸肌的水分比腿肌高，腿肌粗蛋白質(zhì)含量比胸肌高。吉文林等 [7] 對黑羽番鴨屠宰性能和肉質(zhì)的研究結果表明，黑羽番鴨公鴨胸、腿肌的剪切力顯著高于母鴨，公鴨肌肉中（胸、腿肌）pH 值顯著高于母鴨（胸、腿肌）。林謙等 [8] 對不同羽色（黑羽、白羽以及黑白花羽）的公母番鴨進行屠宰性能測定對比分析，發(fā)現(xiàn)雄性黑羽番鴨的屠宰性能較好。錢建中等 [4] 研究結果表明，黑羽番鴨脂肪含量低于北京鴨、櫻桃谷鴨，且黑羽番鴨公鴨胸肌剪切力和腿肌失水率顯著高于母鴨。上述關于黑羽番鴨肉用性能的相關研究僅局限于個體肉質(zhì)分析水平，對于黑羽番鴨胸、腿肌肉質(zhì)性狀的形成機制仍不清楚。

鴨肉的風味主要由鴨肉中的揮發(fā)性風味物質(zhì)，即鴨肉中能產(chǎn)生揮發(fā)性肉香味的化合物成分，經(jīng)加熱后形成的特殊風味，是消費者對鴨肉食品產(chǎn)生的一系列生理反應和感官評定，是影響鴨肉食用口感的一個重要因素 [3，9] 。從基因組宏觀層面對黑羽番鴨肉質(zhì)性狀的調(diào)控機制進行深度挖掘，可以挖掘黑羽番鴨肉質(zhì)風味形成的關鍵基因。本研究以黑羽番鴨為研究對象，通過對黑羽番鴨胸肌、腿肌中揮發(fā)性風味物質(zhì)進行測定，結合RNA-Seq轉(zhuǎn)錄組測序，對黑羽番鴨肉質(zhì)風味前體物質(zhì)特異性沉積機制進行系統(tǒng)、全面的綜合分析，為培育優(yōu)質(zhì)肉用黑羽番鴨提供參考。

1 材料與方法

1.1 試驗材料

本試驗所用的黑羽番鴨由江蘇農(nóng)牧科技職業(yè)學院培育，選擇具有相同遺傳背景、同批孵化、同批出雛的黑羽番鴨，所用個體均在同一環(huán)境下飼養(yǎng)至300日齡。

1.2 樣品采集

在黑羽番鴨大群體中隨機抽樣6只（公鴨3只，母鴨3只），按照國家實驗動物處理行為準則進行屠宰，分別采取同一部位的胸肌、腿肌組織樣本。將同一樣本分為2組：一組迅速放入-80 ℃ 冰柜中速凍，以供測定分析揮發(fā)性風味物質(zhì)含量;另一組放入-80 ℃ 冰箱保存，以供進行轉(zhuǎn)錄組測序。

1.3 試驗方法

1.3.1 揮發(fā)性風味物質(zhì)測定 參考王衛(wèi)等 [10] 方法，利用氣-質(zhì)聯(lián)用儀（Trace DSQ GC-MS）進行黑羽番鴨胸肌、腿肌中揮發(fā)性風味物質(zhì)測定，并獲得各成分相對含量。

1.3.2 RNA提取及轉(zhuǎn)錄組測序 通過Trizol法提取各肌肉樣本組織總RNA，經(jīng)RNA質(zhì)檢后進行轉(zhuǎn)錄組測序 [11] 。

1.3.3 測序數(shù)據(jù)處理及分析 結合測序獲得的數(shù)據(jù)，對比參考基因組，獲得黑羽番鴨胸、腿肌肉組織樣本的特異轉(zhuǎn)錄組序列信息，完成基因組定位分析。在上述比對結果的基礎上，使用htseq-count（v0.6.0）軟件 [12] 統(tǒng)計每個轉(zhuǎn)錄本的表達水平，而后對所有基因進行對比分析，選取FDR（False discovery rate，錯誤發(fā)現(xiàn)率）小于0.05且差異倍數(shù)（FC）大于等于2的基因作為候選差異基因。最后通過GO（Gene ontology，基因本體）進行差異基因功能注釋分析，以及KEGG（Kyoto encyclopedia of genes and genomes， 京都基因與基因組百科全書）進行通路富集分析 [13] 。

1.3.4 qRT-PCR檢測 隨機抽取方法1.3.2中差異表達基因4個，以反轉(zhuǎn)錄得到的黑羽番鴨胸、腿肌肉組織（樣本號1～ 6）的cDNA為模板進行qRT-PCR（表1）檢測，檢驗上述基因在黑羽番鴨胸、腿肌肉組織的表達規(guī)律與測序分析結果是否一致。

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

qRT-PCR結果分析采用 2 -△△ Ct? 方法，SPASS 20.0軟件對試驗所有數(shù)據(jù)進行 t 檢驗，以 P<? 0.05視為差異顯著， P<? 0.01視為差異極顯著 [14] 。

2 結果與分析

2.1 黑羽番鴨胸肌、腿肌中揮發(fā)性風味物質(zhì)含量

對黑羽番鴨胸肌、腿肌中的揮發(fā)性風味物質(zhì)進行分離鑒定，通過NIST數(shù)據(jù)庫進行檢索，而后結合參考文獻[15]，篩選出可能與黑羽番鴨肉質(zhì)風味相關的7種揮發(fā)性風味物質(zhì)。由表2可知，黑羽番鴨胸肌中辛醛、2，3-辛二酮含量顯著高于腿肌，3-羥基-2-丁酮含量顯著低于腿肌，而正戊醇、甲基庚烯酮、壬醛以及正庚醇在胸、腿肌中含量無顯著差異。

2.2 黑羽番鴨胸肌、腿肌轉(zhuǎn)錄組測序分析

共構建300日齡公、母黑羽番鴨的胸肌、腿肌各3個組織樣品的轉(zhuǎn)錄組文庫，對測序初始文庫進行質(zhì)量控制篩選，去除帶接頭、無法確定堿基信息以及低質(zhì)量的測序數(shù)據(jù)（reads），獲得的篩選數(shù)據(jù)結果如表3所示。300日齡公、母黑羽番鴨的胸肌、腿肌12個樣品測序后獲得39 527 132～ 54 977 176 條初始數(shù)據(jù)，而后經(jīng)去除接頭以及低質(zhì)量數(shù)據(jù)篩選，獲得37 824 878～ 52 542 526 條有效數(shù)據(jù)。過濾后有效數(shù)據(jù)占原始數(shù)據(jù)的比例大于94%，表明本試驗構建的黑羽番鴨胸肌、腿肌12個樣品文庫質(zhì)量較好。另外，過濾后有效數(shù)據(jù)中準確率在99.9%以上的堿基總數(shù)占93%以上，進一步說明本試驗測序結果質(zhì)量較好。綜上所述，本試驗構建的公、母黑羽番鴨的胸肌、腿肌12個組織轉(zhuǎn)錄組文庫質(zhì)量較好，可用于后續(xù)差異基因功能分析，保證了研究結果的可靠性。

2.3 黑羽番鴨腿肌與胸肌差異基因篩選及功能分析

通過FPKM值 [9] 計算轉(zhuǎn)錄組測序所捕獲的基因參考表達值，以黑羽番鴨腿肌組織為對照，黑羽番鴨胸肌組織為試驗組，計算所捕獲基因的差異倍數(shù)（Fold change值）。為聚焦關鍵基因，篩選黑羽番鴨腿肌對比黑羽番鴨胸肌基因表達差異值|log2 ??Fold change |≥1（ P < 0.05）的差異基因作為黑羽番鴨胸肌、腿肌組織的候選差異基因，而后根據(jù)差異基因的 Fold change 對數(shù)值以及統(tǒng)計學顯著程度繪制候選差異基因的火山圖（圖1），以標示黑羽番鴨胸肌、腿肌組織差異基因的分布情況。最終篩選出614個差異基因，其中171個差異基因上調(diào)表達（黑羽番鴨腿肌對比黑羽番鴨胸?。?，443個差異基因下調(diào)表達（黑羽番鴨腿肌對比黑羽番鴨胸?。?。

2.4 黑羽番鴨腿肌與胸肌差異基因GO功能富集分析

為探究黑羽番鴨腿肌與黑羽番鴨胸肌差異基因的功能，結合GO數(shù)據(jù)庫對上述篩選獲得的614個顯著差異表達基因（|log2 ??Fold change |≥1（ P < 0.05）進行基因功能富集分類，差異基因分別富集至分子功能（Molecular function，MF）、細胞組分（Cellular component，CC）以及生物過程（Biological process，BP）3個方面。為進一步聚焦614個顯著差異基因在黑羽番鴨胸、腿肌肉組織風味物質(zhì)形成中的作用，篩選在上述3類分組中前10的功能亞類（圖2）。在MF功能中，黑羽番鴨胸、腿肌肉組織的差異表達基因在無機分子實體跨膜轉(zhuǎn)運蛋白活性（Inorganic molecular entity transmembrane transporter activity）與底物特異性通道活性（Substrate-specific channel activity）中所占比例最大，分別聚集到33和24個差異表達基因。在BP功能中，細胞通訊的調(diào)控（Regulation of cell communication）與信號調(diào)節(jié)（Regulation of signaling）富集到的差異基因數(shù)目最多，各有60個差異基因。在CC功能中，差異基因在質(zhì)膜組成成分（Integral component of plasma membrane）中所占的比例最高（46個差異基因），信號調(diào)節(jié)（Regulation of signaling）次之（36個差異基因）。

2.5 黑羽番鴨腿肌與胸肌差異基因KEGG通路注釋

為深入解析差異基因在黑羽番鴨胸肌、腿肌風味物質(zhì)形成中的作用，通過KEGG信號通路數(shù)據(jù)庫對上述篩選獲得的614個顯著差異表達基因進行通路注釋，共富集到109條差異表達的信號通路。以差異表達基因在此Pathway全部基因中所占的比例（GeneRatio）和 p 值進行篩選。GeneRatio越大，表示富集的程度越大; p 值越接近于0，表示富集越顯著。由此篩選出富集最顯著的10條差異信號通路（圖3），包括焦點黏連（Focal adhesion）信號通路、PI3K-Akt（PI3K-Akt signaling pathway）信號通路以及Hippo（Hippo signaling pathway）信號通路等。

2.6 黑羽番鴨風味物質(zhì)相關的差異表達基因篩選

結合上述測序結果、差異基因篩選及功能分析，獲得多個與風味物質(zhì)形成有關的基因存在顯著差異，表4列舉出黑羽番鴨胸肌、腿肌轉(zhuǎn)錄組中與風味物質(zhì)相關的差異表達前20個基因（TOP20）及其功能描述。根據(jù)GO分析以及KEGG功能注釋，其中基因LOC101793529、LOC101794721、NIM1K和SHMT1與氨基酸形成有關，基因 GPI 、 LDHB ?、PFKFB3、 PKM 參與糖代謝生物過程，剩余大部分差異基因與機體脂肪、脂肪酸代謝相關，例如基因ACOT11、ACSBG1、FABP3、ELOVL6和LOC101790670等。

為了進一步研究這些差異基因之間的互作關系，使用STRING（v.11.5）數(shù)據(jù)庫構建蛋白質(zhì)互作網(wǎng)絡分析圖。在黑羽番鴨胸肌、腿肌組織風味物質(zhì)相關差異基因TOP20中，共鑒定出15種已知編碼蛋白質(zhì)，其中13種蛋白質(zhì)被發(fā)現(xiàn)參與2組蛋白質(zhì)相互作用（圖4）。

2.7 黑羽番鴨腿肌與胸肌差異表達基因熒光定量驗證

為驗證黑羽番鴨胸、腿肌轉(zhuǎn)錄組差異基因測序結果，分別選擇表4中黑羽番鴨腿肌高表達基因LOC101790670、PFKFB3和胸肌高表達基因LOC101792626、FABP3以 GADPH 為內(nèi)參基因進行qRT-PCR 驗證。qRT-PCR結果（圖5）顯示：上述4個差異基因在黑羽番鴨胸肌和腿肌組織中的上、下調(diào)表達規(guī)律與本試驗通過黑羽番鴨胸肌和腿肌組織轉(zhuǎn)錄組測序獲得的上、下調(diào)表達規(guī)律一致，證實本試驗轉(zhuǎn)錄組測序獲得的結果較為準確，所篩選出的差異基因可信，可用于后續(xù)的功能驗證分析。

3 討 論

黑羽番鴨生長速度快、體型適中，其鴨肉中瘦肉含量高、肉質(zhì)鮮美，又因其具有獨特的野禽風味而深受廣大消費者喜愛 [1-2] 。本課題組前期研究結果表明，黑羽番鴨腿肌中膽固醇含量以及肌苷酸含量均低于胸肌 [3] 。而有關畜禽肉中風味物質(zhì)研究結果表明，膽固醇、肌苷酸和揮發(fā)性風味物質(zhì)等是形成其肉質(zhì)風味的主要來源 [15] ?；谇捌谘芯拷Y果，本研究進一步對黑羽番鴨胸肌、腿肌中的揮發(fā)性風味物質(zhì)進行對比分析，并結合RNA-Seq技術挖掘黑羽番鴨肉質(zhì)風味的調(diào)控基因，為進一步培育優(yōu)質(zhì)肉用黑羽番鴨以及探究肉質(zhì)風味調(diào)控機制提供參考。

禽肉中醛類、醇類及酮類含量不同可能是造成禽肉風味有差異的原因 [16] 。本試驗對黑羽番鴨胸肌、腿肌中揮發(fā)性風味物質(zhì)進行分離鑒定，篩選出可能與黑羽番鴨肉質(zhì)風味相關的7種揮發(fā)性風味物質(zhì)，其中辛醛、2，3-辛二酮含量在黑羽番鴨胸肌中顯著高于腿肌，3-羥基-2-丁酮含量顯著低于腿?。?P < 0.05），而正戊醇、甲基庚烯酮、壬醛以及正庚醇在胸肌和腿肌中含量無顯著差異。辛醛由油酸氧化而形成，具有堅果香、杏仁香味，對肉香味的貢獻較大，在黑羽番鴨胸肌中含量（1.13 μg/kg ）顯著高于腿?。?.04 μg/kg ）。2，3-辛二酮是脂質(zhì)氧化的另一重要產(chǎn)物，在含量低時有果香、奶香味，黑羽番鴨胸肌中含量（6.88 μg/kg ）顯著高于腿肌（1.25 μg/kg ），說明胸肌中酮類物質(zhì)含量較高，酮類物質(zhì)又可結合脂質(zhì)產(chǎn)生焦糖味 [17] ，從而有益于胸肌肉質(zhì)的風味提升。3-羥基-2-丁酮是某些風味物質(zhì)形成過程中的中間產(chǎn)物，具有甜香、奶香和油膩脂肪香味，黑羽番鴨胸肌中含量（1.06 μg/kg ）顯著低于腿?。?2.36 μg/kg ），但是由于其香味閾值高，需要在高含量下才能產(chǎn)生嗅覺反應 [15] ，因此對黑羽番鴨肉質(zhì)風味的整體貢獻不高。黑羽番鴨胸肌與腿肌揮發(fā)性風味物質(zhì)對比后，篩選出辛醛、2，3-辛二酮以及3-羥基-2-丁酮這3類含量存在顯著差異的物質(zhì)，這可能與黑羽番鴨胸肌與腿肌肉質(zhì)風味區(qū)別相關。

進一步對黑羽番鴨胸肌、腿肌進行轉(zhuǎn)錄組對比分析，篩選出614個差異基因，其中上調(diào)差異表達基因171個，下調(diào)差異表達基因443個。結合GO基因功能分析以及KEGG富集分析，篩選出黑羽番鴨胸肌、腿肌轉(zhuǎn)錄組中與風味物質(zhì)相關的TOP20顯著性差異基因。TOP20顯著性差異基因參與機體脂肪代謝相關的差異基因占絕大多數(shù)：與胸肌相比，黑羽番鴨腿肌中大多數(shù)基因ACOT11、ACSBG1、APOA1、APOD、ELOVL6、FABP3和 SCD 等均呈現(xiàn)下調(diào)表達，少數(shù)基因 CRAT 、LOC101790670和LOC101792317等呈上調(diào)表達。?；o酶A硫酯酶11（ACOT11）作為?；o酶A硫酯酶（ ACOT ）家族成員之一，是一類催化脂肪?；o酶A水解形成游離脂肪酸和輔酶A的酶，在脂類的合成與分解過程中發(fā)揮著巨大的作用，能減少能量消耗同時保存體內(nèi)熱量 [18] 。ACOT11基因能促進小鼠體內(nèi)游離脂肪酸濃度升高 [19] ，將該基因敲除后，小鼠脂肪酸氧化效率增加，脂質(zhì)沉積減少 [20] ，可以抵抗飲食導致的肥胖以及肝臟脂肪變性 [21] 。脂肪酸結合蛋白3（Fatty acid-binding protein 3，F(xiàn)ABP3）作為脂肪酸結合蛋白家族中的重要成員，可以通過長鏈脂肪酸的攝取、轉(zhuǎn)運等調(diào)節(jié)體內(nèi)脂肪和葡萄糖代謝過程，從而影響肌內(nèi)脂肪沉積 [22-23] 。相關研究結果表明，F(xiàn)ABP3的表達水平與30～ 90日齡的哈薩克綿羊肌肉組織中的脂肪含量呈顯著正相關 [24] ，與拜城油雞肌肉組織中的脂肪含量相關 [25] ，與6周齡北京鴨的腹脂率呈顯著負相關關系 [26] ;在前體脂肪細胞中過表達可以有效促進脂肪合成 [27-28] 。

上調(diào)表達基因表達的蛋白質(zhì)中，GPI（Glucose phosphate isomerase， 葡萄糖-6-磷酸異構酶）是一種多功能蛋白酶， 主要參與糖酵解和糖異生這2個生物過程，能夠催化葡萄糖-6-磷酸和果糖-6-磷酸之間的可逆反應 [29] ，目前研究多集中于醫(yī)學方面，GPI與非球型血紅細胞貧血癥、類風濕關節(jié)炎等疾病密切相關 [30-31] 。PFKFB3（6-phosphofructo-2-kinase/fructose-2，6-biphosphatase 3，6-磷酸果糖2-激酶/果糖-2，6-二磷酸酶3）參與生物體糖酵解反應，調(diào)節(jié)糖代謝通路和能量代謝過程 [32] ，研究證實，PFKFB3與腫瘤、血管鈣化、糖尿病等疾病的發(fā)生相關 [33-36] 。其余上調(diào)表達基因如LOC101790670、LOC101792317等，從GO功能注釋中可發(fā)現(xiàn)與脂肪代謝相關。但是，上述這些基因尚未有與肉質(zhì)風味相關的研究報道。

肉質(zhì)風味形成的主要原因之一是肌肉組織中含有大量的脂質(zhì)，其風味的差異與肌內(nèi)脂肪有著一定的關系。肌內(nèi)脂肪作為風味物質(zhì)的主要前體之一，在體內(nèi)合成與分解的過程中形成大量的揮發(fā)性化合物，如酮類、醇類、酸類等，是揮發(fā)性風味物質(zhì)的主要來源 [24] 。本試驗發(fā)現(xiàn)黑羽番鴨胸肌中揮發(fā)性風味物質(zhì)較腿肌中高，且通過轉(zhuǎn)錄組篩選獲得大量參與機體脂肪代謝相關的差異基因。通過蛋白質(zhì)網(wǎng)絡互作圖分析發(fā)現(xiàn)，這些差異基因的表達蛋白ACSBG1、APOA1、APOD、ELOVL6、FABP3、SCD和CRAT呈現(xiàn)出清晰的互作效應，這些基因可能共同作用于黑羽番鴨胸肌和腿肌的脂肪沉積過程，從而影響黑羽番鴨肉質(zhì)風味，其具體的調(diào)控機制有待結合黑羽番鴨胸肌和腿肌的脂肪含量及各類脂肪酸比例進行深入研究。

4 結 論

本研究在黑羽番鴨胸肌、腿肌中揮發(fā)性風味物質(zhì)對比基礎上，結合RNA-Seq技術對黑羽番鴨肉質(zhì)風味前體物質(zhì)特異性沉積機制進行了系統(tǒng)、全面的綜合分析，篩選出的ACOT11、APOA1、FABP3等脂肪代謝相關差異基因可能與黑羽番鴨肉質(zhì)風味形成相關。本研究結果為深入探討黑羽番鴨肉質(zhì)風味形成機理以及優(yōu)質(zhì)肉用黑羽番鴨的選育提供了參考依據(jù)。

參考文獻：

[1] 唐現(xiàn)文，張響英，陸艷鳳. ?MC4R 在黑羽番鴨肌肉組織中的表達研究[J].中國家禽，2020，42（7）：116-120.

[2] 孫國波，王 健，陸艷鳳，等. 番鴨羽色與產(chǎn)肉性狀的觀測研究[J].江蘇農(nóng)業(yè)科學，2019，47（2）：173-176.

[3] 孫國波，王日君，董 飚，等. 黑羽番鴨及繼代選育后代的胸腿肌營養(yǎng)成分測定[J].貴州農(nóng)業(yè)科學，2015，43（1）：104-107.

[4] 錢建中，段修軍，卞友慶，等. 不同性別黑羽番鴨屠宰性能、常規(guī)肉品質(zhì)分析[J].江蘇農(nóng)業(yè)科學，2014，42（2）：164-166.

[5] 董 飚，段修軍，孫國波，等. ?GH 基因?qū)谟鸱嗴w重和屠宰性能的遺傳效應分析[J].浙江農(nóng)業(yè)學報，2015，27（1）：22-27.

[6] 王錦鋒，包小成. 櫻桃谷鴨、白羽番鴨、黑羽番鴨和蘇牧鴨肉品質(zhì)的比較研究[J].畜牧與獸醫(yī)，2010，42（2）：40-42.

[7] 吉文林，段修軍，董 飚，等. 黑羽番鴨屠宰性能及肉品質(zhì)的研究[J].西南農(nóng)業(yè)學報，2013，26（2）：795-797.

[8] 林 謙，吳買生，蔣桂韜，等. 不同羽色和性別番鴨屠宰性能及肌肉成分比較研究[J].家畜生態(tài)學報，2014，35（1）：30-34.

[9] 章琳俐，李 麗，朱志明，等. 基于RNA-seq鑒定連城白鴨肉質(zhì)風味相關候選基因[J].農(nóng)業(yè)生物技術學報，2021，29（4）：711-722.

[10] 王 衛(wèi)，劉文龍，覃 瑤，等. 無抗放養(yǎng)模式有利于固始雞揮發(fā)性風味物質(zhì)的形成[J].湖南農(nóng)業(yè)科學，2017（6）：70-73.

[11] 張 蕾，章敬旗，章瑋月，等. 泰州鵝開產(chǎn)前后肝臟組織轉(zhuǎn)錄組對比分析[J].福建農(nóng)業(yè)學報，2020，35（9）：987-996.

[12] ROBINSON M D， MCCARTHY D J， SMYTH G K. edgeR： a Bioconductor package for differential expression analysis of digital gene expression data[J]. Bioinformatics， 2010，26（1）： 139-140.

[13] YU G C， WANG L G， HAN Y Y， et al. clusterProfiler： an R package for comparing biological themes among gene clusters[J]. Omics： a Journal of Integrative Biology，2012，16（5）：284-287.

[14] 唐現(xiàn)文，張 蕾，章敬旗，等.高郵鴨卵巢中雙黃蛋性狀相關lncRNA表達差異分析[J].甘肅農(nóng)業(yè)大學學報，2021，56（2）：40-49，60.

[15] 孫旺斌，付 琪，薛瑞林，等. 不同棗粉水平對陜北白絨山羊肉揮發(fā)性風味物質(zhì)的影響[J].動物營養(yǎng)學報，2021，33（10）：1-13.

[16] 崔小燕，茍鐘勇，蔣守群，等. 雞肉風味的形成機制與調(diào)控研究進展[J].動物營養(yǎng)學報，2019，31（2）：500-508.

[17] 李柯呈，徐寶才，姚 忠，等.干腌鹽量對南京鹽水鴨特征風味成分的影響[J].食品與發(fā)酵工業(yè)，2019，45（17）：98-104.

[18] 張 麗，強 俊，徐 跑，等. ?；o酶A硫酯酶11基因（ACOT11）及其家族的研究進展[J].江蘇農(nóng)業(yè)科學，2021，49（2）：12-18.

[19]? ZHANG Y Z， LI Y X， MICHELE W N， et al. Targeted deletion of thioesterase superfamily member 1 promotes energy expenditure and protects against obesity and insulin resistance[J].Proceedings of the National Academy of Sciences of the United States of America，2012，109（14）： 5417-5422.

[20]? TAKAYUKI O， KANAKO I， RYOTA M， et al. Chronic high-fat feeding impairs adaptive induction of mitochondrial fatty acid combustion-associated proteins in brown adipose tissue of mice[J]. Biochemistry and Biophysics Reports，2017，10：32-38.

[21] 胡曉維，靳聰飛，劉新峰，等. 牛ACOT11基因克隆及其編碼蛋白生物信息學分析[J].天津農(nóng)學院學報，2015，22（3）：1-6.

[22] 艾錦新，龍安炬，羅衛(wèi)星，等. 黔北麻羊FABP3基因多態(tài)性及其與生長性狀的相關性研究[J].畜牧與獸醫(yī)，2021，53（1）：13-18.

[23]? BI Y L， YUAN X Y， ZHU P F， et al. A novel long noncoding RNA， ENSGALG00000021686， regulates the intracellular transport of fatty acids by targeting the FABP3 gene in chicken[J]. Biochemical and Biophysical Research Communications，2020，528（4）：706-712.

[24] HUANG Z G， XIONG L， LIU Z S， et al. The develop-mental changes and effect on IMF content of H-FABP and PPAR-PPAR-γ mRNA expression on sheep muscle[J].Asia Cenetica Sinica，2006（6）：507-514.

[25] 沙爾山別克·阿不地力大，舒 婷，吳旭升，等. 拜城油雞 A-FABP、H-FABP和PPAR-PPAR-γ 基因mRNA表達的發(fā)育性變化及其與肌內(nèi)脂肪含量的相關性研究[J].中國家禽，2019，41（14）：5-11.

[26] 榮光輝，徐 琪，徐鐵山，等. 脂肪型和瘦肉型北京鴨 A-FABP 基因組織表達特征[J].中國畜牧雜志，2015，51（1）：5-9.

[27] YI BAO， WANG J G， WANG S， et al. Overexpression of Banna mini-pig inbred line fatty acid binding protein 3 promotes adipogenesis in 3T3-L1 preadipocytes.[J]. Cell Biology International， 2014， 38（8）：918-923.

[28] 葛 凱. 脂肪細胞大小及脂肪代謝調(diào)控巢湖鴨肌內(nèi)脂肪的分子機理[D].安徽農(nóng)業(yè)大學，2019.

[29] 韓 龍，杜翠紅.葡萄糖-6-磷酸異構酶研究進展[J].藥物生物技術，2012，19（6）：561-564.

[30] 白小娟，張麗卿.葡萄糖-6-磷酸異構酶在類風濕關節(jié)炎中的研究進展[J].山西醫(yī)藥雜志，2018，47（6）：668-670.

[31] 黃 鵬，唐 黎，汪愛平，等. ?GPI 基因突變致葡萄糖磷酸異構酶缺乏癥1例報告并文獻復習[J].臨床兒科雜志，2020，38（10）：785-788.

[32]? GRIESEL B A， MATSUZAKI S， BATUSHANSKY A， et al. PFKFB3-dependent glucose metabolism regulates 3T3-L1 adipocyte development[J]. The FASEB Journal， 2021， 35（7）： 1-19.

[33] 李 駿，張 倩，江東根，等. PFKFB3基因在前列腺癌中的表達及其對前列腺癌細胞糖酵解及生長的影響[J].中華腔鏡泌尿外科雜志（電子版），2020，14（6）：460-465.

[34] 饒利棟，鄧雪強，易 軒，等.PFKFB3通過調(diào)控PI3K/AKT信號通路促進骨肉瘤細胞的增殖及轉(zhuǎn)移[J].腫瘤，2021，41（2）：77-90.

[35] 牛 津，張 敏，王巧娟，等.PFKFB3抑制劑減輕高磷誘導的大鼠血管平滑肌細胞鈣化[J].心臟雜志，2021，33（2）：117-122.

[36] 袁 姿. 血清TLR4和PFKFB3水平與2型糖尿病及肥胖的相關性研究[D].蘇州：蘇州大學，2019.

（責任編輯：張震林）

收稿日期：2021-09-02

基金項目：泰州市科技支撐項目（TN201915）;江蘇省農(nóng)業(yè)重大新品種創(chuàng)制項目（PZCZ201737）

作者簡介：段修軍（1972-），男，湖南南縣人，碩士，教授，研究方向為動物遺傳育種與繁殖。（E-mail）sgb1981@126.com