膽囊收縮素A型受體基因序列變異和表達與略陽烏雞飼料轉化率的關聯性研究

王哲鵬，周雯馨，賀俊錫，虎巧燕，趙家悅

膽囊收縮素A型受體基因序列變異和表達與略陽烏雞飼料轉化率的關聯性研究

王哲鵬，周雯馨，賀俊錫，虎巧燕，趙家悅

西北農林科技大學動物科技學院，陜西楊凌 712100

【背景】食欲和飼料轉化率的提高構成了現代肉雞快速生長的生理基礎。膽囊收縮素A型受體（CCKAR）是介導膽囊收縮素（CCK）飽食信號和消化調節作用的關鍵受體。鑒于CCK信號在能量平衡調節中的重要作用，CCKAR的表達水平與雞的生長性能密切相關。然而，CCKAR水平調節雞生長的確切機制尚不清楚。【目的】通過研究闡明CCKAR基因表達水平和序列變異與雞飼料轉化率的關系，為雞生長性狀遺傳調控研究和飼料轉化率選育奠定基礎。【方法】測定略陽烏雞12—16周齡料肉比（公雞n=62，母雞n=67）和26—30周齡料蛋比（n=150）。用PCR-RFLP法檢測CCKAR基因內兩個同義突變（g.73206714T＞C和g.73209189C＞T）的基因型。以性別為固定效應，用雙因素方差分析檢驗g.73206714T＞C和g.73209189C＞T與料肉比的關聯性，用單因素方差分析檢驗2個SNP與料蛋比的關聯性。用JASPAR數據庫預測g.73206714T＞C對轉錄因子結合位點產生的影響。以g.73206714T＞C為標記，選取10只雜合子（CT）公雞，用PCR-RFLP法檢測CCKAR基因在g.73206714T＞C位點等位基因間的表達差異。選取30只帶有采食量和料肉比記錄的公雞，用qPCR檢測CCKAR基因在下丘腦、腦垂體、胰腺、膽囊和十二指腸中的表達水平，分析CCKAR表達水平與采食量和料肉比的相關性。【結果】性別效應對略陽烏雞料肉比有顯著影響（=29.44，<0.0001），公雞的料肉比（5.41±0.76）顯著低于母雞（6.62±1.00）。在兩個同義突變中，只有g.73206714T＞C與料肉比顯著關聯（=8.44，=0.0004），其中，TT型的料肉比（5.18±0.63）顯著低于CC（6.37±1.1）和CT型（6.17±0.95）。CC型和CT型的料肉比無顯著差異。g.73206714T＞C與性別的互作效應對料肉比的影響不顯著。G.73206714T＞C對28個轉錄因子結合位點有潛在影響。等位基因特異性表達試驗顯示，T等位基因的轉錄本豐度顯著高于C等位基因，表明CCKAR在等位基因間的表達差異很可能是構成g.73206714T＞C與料肉比關聯性的遺傳基礎。略陽烏雞26—30周齡的平均料蛋比為2.85±0.41。g.73206714T＞C和g.73209189C＞T的基因型效應對料蛋比的影響均不顯著。在5個組織中，CCKAR基因在胰腺中的表達水平最高，且與采食量呈負相關（=-0.45,=0.017）。但是，在其他組織中，CCKAR的表達水平與采食量的關聯性均不顯著。對于料肉比性狀，CCKAR基因在胰腺（=-0.41，=0.03）和腦垂體（=-0.57,=0.0018）中的表達水平越高，料肉比越低。在其他組織中，這種負相關性均不顯著。【結論】研究發現CCKAR基因表達水平上調有利于提高雞的飼料轉化率。G.73206714T＞C的T等位基因是一個與CCKAR基因高表達活性和高飼料轉化率顯著關聯的分子標記。

雞；生長；CCKAR；飼料轉化率；食欲

0 引言

【研究意義】鑒定影響雞飼料轉化率的候選基因對提高飼料的利用率，降低養殖成本，節約飼料用糧，緩解人畜爭糧的矛盾具有重要意義。【前人研究進展】膽囊收縮素（CCK）是一種主要由小腸Ⅰ型內分泌細胞分泌的肽類激素，具有誘導食欲調控中樞產生飽腹感，促進膽囊收縮、胰酶分泌，抑制嗉囊排空、胃酸分泌、調節胃腸蠕動等多種生理功能[1]。膽囊收縮素A型受體（CCKAR）是介導CCK信號的關鍵受體[2]。在人類中，CCKAR轉錄本異常剪切與肥胖和膽結石顯著關聯[3]。在大鼠中，由基因突變引起的CCKAR缺失將導致貪吃、肥胖和糖尿病[4]。在雞中，CCKAR基因表達下調能顯著提高生長速度[5]。CCKAR拮抗劑（對CCK與受體的結合產生競爭性抑制）處理能刺激食欲、提高采食量[6-7]，而CCKAR類似物處理會減少采食量[8]，表明CCKAR水平對哺乳動物能量平衡的調節極有可能通過食欲調控途徑來實現。但是，CCKAR拮抗劑、CCK類似物和CCK釋放促進劑處理并不能改變雞的采食量，而是對膽汁和胰淀粉酶的分泌、嗉囊的排空有明顯調節作用[9-10]。Covasa等發現CCKAR拮抗劑只有在高劑量（＞90 μg·kg-1BW）添加的情況下才能提高雞的采食量，但CCKAR拮抗劑預處理并不能有效阻斷內源CCK誘導的厭食效應[11]。這些研究對CCK信號在雞中的食欲調節作用提出質疑，但支持CCK信號與消化液分泌、嗉囊排空等消化相關的生理活動存在密切聯系。【本研究切入點】鑒于CCK信號對消化的調節作用及CCKAR的信號轉導功能，本研究推測CCKAR水平可能通過營養物質消化途徑影響雞的飼料轉化率。Xu等發現CCKAR在低剩余采食量組特異表達，支持CCKAR是影響雞飼料轉化率的一個候選基因[12]。【擬解決的關鍵問題】測定了略陽烏雞12—16周齡的料肉比和26—30周齡的料蛋比，檢測了CCKAR基因在下丘腦、腦垂體、胰腺、膽囊和十二指腸中的表達水平，研究了CCKAR基因表達水平及序列變異與飼料轉化率的關聯性，旨在闡明CCKAR基因遺傳差異與飼料轉化率的聯系，為從育種和CCKAR活性干預角度提高雞的飼料轉化率奠定理論基礎。

1 材料與方法

1.1 略陽烏雞的飼養和飼料轉化率測定

略陽烏雞是產自陜西省略陽縣的地方雞種。具有類似于絲羽烏雞的真皮超黑色素（dermal hyperpigmentation）特征，是我國一種體形較大的烏雞品種[13]。由于缺乏系統選育，略陽烏雞在產蛋和生長性能方面存在較大群體變異性，表明一些調控生產性能的功能突變很可能處于分離狀態[14]。這為功能突變的鑒定創造了條件。

在略陽烏雞蛋用系育種核心群三世代41個半同胞家系中隨機選取62只公雞和67只母雞，測定料肉比。在四世代56個半同胞家系中隨機選取150只母雞，測定料蛋比。料肉比和料蛋比測定所用樣本為同一批次孵化，在同一棟雞舍飼養。籠養育雛，全價料飼喂，自由采食、飲水。育雛期采用18—22 h光照，育成期采用自然光照，產蛋期采用16 h恒定光照。在60 d時將雞從育雛舍轉入產蛋雞舍，單籠飼養。三世代在2019年2月25日出雛，2020年1月8日淘汰，四世代在2019年12月9日出雛，2020年12月21日淘汰。三世代雞的性別用PCR法鑒定，用5′-GTAGTCTGTGTCTTATCATAGC-3′和5′-GGT CTCTGTACTTGGCATT-3′ 引物檢測個體是否攜帶W染色，用5′-GCCATGCATGTCTAAGTAC-3′和5′-G CACAGACAGTACCATCG -3′ 擴增18 s rRNA片段，確保PCR成功進行[15]。

以個體為單位，測定料肉比和料蛋比。2只雞之間隔1個空籠位，避免互相采食對方的飼料。每周測定初始料重、初始體重，周末料重、周末體重和產蛋量，按（期末料重-初始料重）/（期末體重-初始體重）計算料肉比，按（期末料重-初始料重）/試驗期內產蛋量計算料蛋比。料肉比測定周期為12—16周齡，測定時間為2019年5月22日至2019年6月19日；料蛋比測定周期為26—30周齡，測定時間為2020年6月8日至2020年7月5日。略陽烏雞養殖和料肉比、料蛋比測定均在陜西略陽龍佳農業發展有限公司完成。

1.2 CCKAR序列變異基因型檢測

在飼料轉化率測定結束后，翅靜脈采血0.5 mL，用ACD（血：ACD=3﹕1）抗凝。ACD的配方為0.48 g檸檬酸，1.32 g檸檬酸三鈉，1.47 g葡萄糖，用雙蒸水溶解后，定容至100 mL。取20 μL抗凝血，用Blood Genome DNA Extraction試劑盒（Takara）按試劑盒說明提取基因組DNA。用1 %瓊脂糖凝膠電泳檢測DNA的完整性，用NanoDrop? 2000（Thermo）分光光度計測定DNA濃度，將DNA濃度調至50 ng·μL-1備用。

g.73206714T＞C （dbSNP登錄號：rs315131298）和g.73209189C＞T（dbSNP登錄號：rs317682933）是分布位于CCKAR基因第3和第5外顯子內的兩個同義突變。研究發現g.73206714T＞C和g.73209189C＞T與雞的生長速度顯著關聯[5]。本研究用PCR-RFLP法檢測這兩個SNP的基因型。PCR體系由1 μL DNA（50 ng·μL-1），0.5 μL正鏈引物（10 μmol·L-1），0.5 μL反鏈引物（10 μmol·L-1），10 μL預染2 × Taq Plus MasterMix (康為) 和8 μLddH2O組成。g.73206714T＞C位點的PCR擴增引物為5′-AGCACGGAGACC AAGACCCA-3′和5′-GCAAATGGCACTGTACCGCT -3′，g.73209189C＞T位點的PCR擴增引物為5′-ATGC TGATGGTGATAGTGGT-3′和5′-GGTGAAAGTGGC TAGAAAAC-3′。PCR擴增條件為94℃預變性3 min，（94℃變性30 s, 60℃退火30 s, 72℃延伸20 s）×33個循環，72℃延伸5 min。PCR在伯樂T100TM熱循環儀上進行。g.73206714T＞C擴增產物用MnlI內切酶（NEB）消化，g.73209189C＞T擴增產物用MboII內切酶（NEB）消化。消化反應體系為1 μL 10× CutSmart Buffer, 1.5 μLPCR產物, 0.6 μL（3 U）內切酶和6.4 μL ddH2O。消化條件為37℃水浴消化1 h。用3.0%（wt/vol）瓊脂糖凝膠電泳分離消化產物，經0.1%（vol/vol）的溴乙錠染色，用Gel DocTMXR+凝膠成像系統（伯樂）拍照記錄結果。電泳條件為90 V，40 min。部分PCR-RFLP基因型檢測結果用Sanger測序進行驗證。

1.3 g.73206714T＞C位點對轉錄因子擾動作用預測

在JASPAR（http://jaspar.genereg.net）數據庫中，選取Vertebrata門類下1 011個轉錄因子，在Scan窗口提交FASTA格式序列，分別搜索T等位基因可結合的轉錄因子和C等位基因可結合轉錄因子，分析g.73206714T＞C對轉錄因子-DNA結合產生的潛在影響。搜索轉錄因子的Relative profile score 閾值設定為80 %。

1.4 CCKAR基因表達量檢測

在參與料肉比測定的62只公雞中隨機選取30只公雞，在16周齡宰殺，采集下丘腦、腦垂體、胰腺、膽囊和十二指腸，剪碎，置RNAstore（康為）4 ℃浸泡過夜，置-80℃冰箱保存備用。在采集膽囊時，需劃破膽囊，棄膽汁，用1×PBS緩沖液（Gibco）漂洗3次后，保存備用。在距幽門約5 cm處采集十二指腸。取約100 mg組織，用1 mL TRNzol Universal 總 RNA 提取試劑（天根），按TRNzol操作說明提取總RNA。用1 %瓊脂糖電泳檢測RNA的完整性，用NanoDrop? 2000（Thermo）分光光度計測定RNA的濃度。取1 μg RNA，用PrimeScriptTM1st Strand cDNA synthesis試劑盒（Takara）按試劑盒說明將RNA反轉錄為cDNA。用qPCR檢測CCKAR基因在各組織中的表達量。qPCR體系為2 μL cDNA，10 μL 2×TB Green Advantage qPCR Premix （Takara）, 0.4 μL正鏈引物（10 μmol·L-1），0.4 μL 反鏈引物（10 μmol·L-1）和7.2 μL ddH2O。qPCR在LightCycler? 96 Instrument實時熒光定量PCR儀（Roche）上進行，擴增條件為95℃ 30 s,（95℃5 s, 60℃20 s 72℃ 15 s）×40個循環。擴增反應結束后執行熔解曲線反應程序，檢測qPCR的特異性，溶解曲線反應程序為95℃ 10 s，65℃ 60 s和97℃ for 1 s。LBr被用作持家基因，校正由cDNA濃度和cDNA模板加樣量差異引起的實驗誤差。CCKAR基因擴增引物：5′-TACAGCAAGCTGGTC CCTTT-3′和5′-AATGAAATGAGGCCATACGC-3′， LBr擴增引物：5′-AGGGCTTCTTGTCACAGGTT-3′和5′-ATCACGAGCTTCTCGATGAACA-3′。目標和持家基因每個樣本均執行3個技術重復。用ΔCt（CtLBr- CtCCKAR）表示CCKAR基因的表達水平，即CCKAR基因表達水平越高，ΔCt值越大。

1.5 g.73206714T＞C位點等位基因特異性表達水平檢測

在30只采集組織樣的公雞中，有10只公雞g.73206714T＞C的基因型為CT。取10只公雞下丘腦、腦垂體、胰腺、膽囊和十二指腸cDNA，執行RT-PCR，獲得包含g.73206714T＞C位點的擴增片段。RT-PCR擴增體系為cDNA模板2 μL，正鏈引物0.5 μL（10 μmol·L-1），反鏈引物0.5 μL（10 μmol·L-1），2× Taq MasterMix（康為）10 μL，ddH2O 7μL。RT-PCR所用引物序列為5′-CTCATTCCCAACCTGCTGAA-3′和5′-GCAAATGGCACTGTACCGCT-3′。RT-PCR產物用MnlI（NEB）內切酶消化。RT-PCR擴增條件，酶切體系、條件，酶切產物電泳檢測方法見材料方法1.2。

1.6 統計分析

在統計分析前，以下四分位數（first quantile）-1.5 ×四分位距（interquantile range）為下限，以上四分位數（third quantile）+ 1.5 ×四分位距（interquantile range）為上限，刪除采食量、增重、產蛋量和CCKAR表達量數據中的異常值。g.73206714T＞C和g.73209189C＞T與料肉比的關聯性用雙因素方差分析進行檢驗，模型為FCRi= μ + sexi+ genotypei+ sexi× genotypei+ ei，FCRi為個體i的料肉比，μ為群體均值，genotypei為基因型效應，sexi×genotypei為性別與基因型的互作效應，ei為誤差項；料蛋比方差分析模型為FCRi= μ + genotypei+ ei，各項含義與料肉比模型相同。CCKAR基因在各組織間的表達差異用單因素方差分析進行檢驗，統計模型為ΔCt = μ + tissue + e。如果主效應統計檢驗的結果為顯著，則進一步執行各水平間最小二乘均數多重比較，多重比較校正方法為Tukey-Kramer。顯著水平為P<0.05。方差分析用SAS University Edition的GLM過程完成。CCKAR表達水平與采食量和料肉比的相關性用R語言Hmisc軟件包分析。

2 結果

2.1 略陽烏雞飼料轉化率描述統計

性別效應對略陽烏雞料肉比有顯著（F=29.44,<0.0001）影響，公雞料肉比5.41±0.76顯著低于母雞（6.62±1.00；圖1）。略陽烏雞26到30周齡平均料蛋比為2.85±0.41。略陽烏雞同期飼料轉化率不僅低于商業蛋雞和肉雞，而且存在較大群體變異性，公雞料肉比變異范圍在3.6—7.4之間（變異系數（C.V.）=14.0 %），母雞在4.8—9.1之間（C.V.= 15.1%），料蛋比變異范圍在2.0—4.4之間（C.V.=14.3%；圖2）。

圖1 略陽烏雞12—16周齡料肉比分布

2.2 CCKAR序列變異與略陽烏雞飼料轉化率的關聯性

本研究分析了CCKAR基因內兩個SNP與略陽烏雞料肉比和料蛋比的關聯性。g.73206714T＞C的基因型效應對略陽烏雞12—16周齡料肉比有顯著影響（= 8.44，=0.0004），TT型料肉比（5.18±0.63）顯著低于CC（6.37±1.1）和CT型（6.17±0.95），CC型和CT型的料肉比無顯著差異（圖2）。g.73206714T＞C與性別的互作效應對料肉比無顯著影響（=0.72，=0.49）。g.73206714T＞C位點的基因型效應對略陽烏雞26—30周齡的料蛋比無顯著影響（=2.16，=0.12；圖2）。在略陽烏雞樣本中，本研究在g.73209189C＞T位點主要檢測到CC和CT兩種基因型。g.73209189C＞T的基因型效應對料肉比（=0.67，=0.41）和料蛋比（=0.48，=0.62）均無顯著影響（圖2）。在料肉比性狀上，g.73209189C＞T與性別也不存在顯著的互作效應（=3.20，=0.07）。

箱線圖上方小寫字母標出多重比較結果，相同字母或未標注字母代表基因型間差異不顯著，不同字母代表基因型間存在差異顯著（P<0.05）

2.3 g.73206714T＞C位點等位基因特異性表達分析結果

關聯分析發現g.73206714T＞C顯著關聯到略陽烏雞料肉比。g.73206714T＞C是CCKAR基因第三外顯子上的一個同義突變。如果g.73206714T＞C直接參與飼料轉化效率調控，g.73206714T＞C可能通過影響CCKAR基因的表達發揮作用。經轉錄因子預測，g.73206714T＞C的T等位基因會影響RORC、Myb、Tcf15、Mycn、HOXA4、GSX2、hoxc4、YY2與DNA的結合，而C等位基因會影響NR5A1、Esrrg、VDR、Dmrt1、HLTF、Nr2e3、ZEB1、FOXK1、SOX10、Foxo1、DMRTA2、Sox5、LBX1、Foxj2、FOXO3、FOXC1、FOXO4、FOXO6、FOXI1、C/EBPα與DNA的結合。

為了進一步闡明g.73206714T＞C與CCKAR基因表達水平間的聯系，本研究通過等位基因特異性表達試驗檢測CCKAR基因在T和C等位基因間的表達差異。本研究以基因組DNA擴增產物為對照，控制由酶切效率產生的潛在誤差。在理論上，以基因組DNA為模板進行PCR擴增時，CT雜合子將產生數量相當的C等位基因和T等位基因擴增產物。PCR-RFLP結果顯示，CT產生138 bp（T等位基因擴增產物）和109 bp（C等位基因擴增產物）兩種亮度接近的電泳條帶（圖3），表明PCR-RFLP能夠準確檢測等位基因特異性擴增子豐度。在對10只雜合仔雞的RT-PCR（模板為cDNA）產物酶切后，除了個別樣本外，多數樣本T等位基因轉錄本豐度均明顯高于C等位基因，說明CCKAR基因在g.73206714T＞C位點存在等位基因表達不平衡性，其中，T等位基因與CCKAR基因高表達活性關聯。

以g.73206714T＞C為標記，選取10只雜合子（CT）個體，用MnlI內切酶消化各組織RT-PCR產物。消化產物用3%瓊脂糖凝膠電泳分離。T等位基因將產生138 bp條帶，C等位基因將產生109 bp和29 bp兩種條帶，29 bp條帶在3%瓊脂糖電泳中不可見。基因組DNA擴增產物被用作對照，驗證MnlI檢測等位基因特異性轉錄本豐度差異的準確性

2.4 CCKAR基因表達水平與料肉比的相關性

為了進一步分析CCKAR基因表達水平與料肉比的關系，本研究檢測了CCKAR基因在下丘腦、腦垂體、胰腺、膽囊和十二指腸中的表達量。CCKAR基因在胰腺中表達水平最高（圖4）。CCKAR基因表達水平除了在十二指腸和下丘腦間無顯著差異外，在其他組織間均存在顯著差異。CCKAR基因在胰腺和腦垂體中的表達水平與料肉比呈現出顯著負相關（圖5）。CCKAR基因在膽囊中的表達水平與料肉比存在潛在的負相關。在十二指腸和下丘腦中，CCKAR表達水平與料肉比間的負相關均未達到統計顯著水平。CCKAR基因在胰腺中的表達水平與采食量存在顯著負相關關系，但在其他組織中，二者的相關關系均不具有統計學意義。

箱線圖上方的字母標出了多重比較的結果，相同字母表示組織間差異不顯著，不同字母表示組織間差異顯著（P<0.05）

3 討論

鑒于CCK信號在食欲調控和能量平衡調節中的作用[16-17]，CCKAR的表達水平與動物的采食量、生長速度及一些代謝疾病的發生密切相關[3-5]。在雞中，CCKAR基因在快速生長型雞中的表達水平顯著低于慢速生長型雞[5]。研究人員推測CCKAR基因低表達促生長的機制很可能與飽食信號閾值改變有關[5]。除了影響食欲外，CCK信號對消化的調節作用表明，CCKAR的表達水平很可能與飼料利用率也有密切聯系。本研究發現CCKAR基因在胰腺、腦垂體和膽囊中的表達水平越高，雞的料肉比越低，為上述推測提供了直接證據。進而，在CCKAR基因第三外顯子上發現同義突變g.73206714T＞C不僅與CCKAR基因的表達活性有關，而且與料肉比顯著關聯。這些結果支持CCKAR基因是一個調節雞飼料轉化率的候選基因，而該基因內一些影響表達的功能變異很可能構成了雞飼料轉化率差異的分子基礎。

3.1 CCKAR基因表達水平影響飼料轉化率的潛在機制

基于相關分析的結果，本研究推測CCKAR水平可能通過以下3種途徑影響飼料轉化率。第一、消化液分泌途徑。CCK作為一種腸肽激素，能刺激胰酶和膽汁分泌[1]。CCKAR拮抗劑能有效阻斷CCK對胰酶分泌的刺激作用，說明CCKAR是介導CCK消化調節作用的特異性受體，CCKAR水平直接影響消化液的分泌[9]。本研究發現CCKAR基因在胰腺中的表達水平與料肉比呈中等強度的負相關，支持CCKAR水平可以通過影響消化液的分泌，影響飼料利用率。第二、脂代謝途徑。將飼料轉化為脂肪是極不經濟的一種做法。相反，提高脂肪的分解效率，抑制脂肪的沉積將有助于提高飼料利用率[18]。膽汁可以通過脂滴乳化促進脂肪的消化，也可以通過白色脂肪棕色化促進脂肪組織的分解[19]。外源CCK注射證明，膽汁的分泌量隨CCK劑量的增加而增加[9]。所以，CCKAR表達水平的升高有可能通過提高膽汁的分泌量，促進脂代謝，提高飼料轉化率。此外，生長激素可以促進脂代謝和同化作用[20]。外源注射CCK能顯著提高生長激素的分泌量[21]。本研究發現CCKAR基因在腦垂體中的表達水平與料肉比存在負相關性。據此推測CCKAR表達上調可能對生長激素的分泌有促進作用，進而提高飼料轉化率。第三、采食量途徑。研究表明，家禽的采食量越低，飼料利用率越高[12,22-23]。Zeng等[23]發現CCK在下丘腦中的表達水平越高，鴨的飼料轉化率越高。本研究發現CCKAR基因在胰腺中的表達水平不僅與采食量呈負相關，而且與料肉比也呈負相關。這些結果支持CCK系統也可能通過抑制雞的采食量，提高飼料轉化率。

ADFI_16wk代表16周齡日平均采食量，FCR_16wk代表16周齡料肉比。圖上方的數字是相關系數和相關系數顯著性檢驗的P值

3.2 CCKAR基因表達水平對雞食欲的調節作用

在腸-腦軸食欲調控系統中，CCK作為一種小腸分泌的飽食信號分子，它通過與迷走神經中的受體結合向食欲調控中樞（下丘腦弓狀核）傳遞飽食信號，也可以直接與食欲調控神經元表面的受體結合產生飽食信號，調節動物的食欲和采食量[24-25]。CCKAR作為CCK飽食信號轉導的關鍵分子，其表達水平與采食量密切相關[5]。本研究在胰腺中檢測到CCKAR表達水平與采食量呈負相關，但是在食欲調控關鍵組織（下丘腦和十二指腸）中均未檢測到這種負相關性。鑒于CCKAR基因在胰腺中表達水平最高，本研究推測胰腺很可能是發出CCK飽食信號的關鍵外周器官。

3.3 g.73206714T＞C對CCKAR基因表達的調控機制

本研究發現CCKAR基因的表達水平與雞飼料轉化率顯著相關，而鑒定表達調節突變是從分子層面揭示CCKAR對飼料轉化率遺傳調控機制的關鍵。g.73206714T＞C是一個與雞的體重顯著關聯的突變[5]。本研究進一步發現這個變異不僅與料肉比顯著關聯，而且與CCKAR基因的表達水平也密切相關。g.73206714T＞C對28個轉錄因子與DNA的結合產生潛在影響。這些轉錄因子的功能涉及脂代謝、性別決定、神經、軟骨、肌肉等組織的發育等多個生物學過程[26-31]。其中，FOXO1, C/EBPα和FOXK1是3個潛在的CCKAR表達調節蛋白。FOXO1通過調節食欲調控神經肽在下丘腦弓狀核中的表達來調控動物的食欲，而C/EBPα和FOXK1通過調節能量代謝相關蛋白FTO、PPARγ、IGF-1、leptin、insulin等的表達，參與能量平衡調節[29-33]。盡管目前并無直接證據表明這些轉錄因子參與CCKAR的表達調控，但是鑒于CCK系統在鳥類食欲和能量平衡調節中的作用[16]，CCKAR很可能也是受這些轉錄因子調控的靶基因。

3.4 g.73206714T＞C與飼料轉化率關聯性對雞育種的意義

通過檢測C等位基因在不同生長性能品種中的頻率分布，研究人員發現C等位基因頻率與雞生長性能存在顯著正相關性[5]。這表明C等位基因很可能是受到正向選擇的一個功能突變。然而，對飼料轉化率而言，本研究發現C等位基因并不是一個優勢等位基因。C等位基因有更低的表達活性，而CCKAR基因低表達將減弱雞飽食信號閾值，增大采食量[5]。這些結果提示，如果以C等位基因為標記，對雞生長性能進行輔助選擇，極有可能是以提高雞采食量的方式促進雞的生長[34]。這種育種策略并不利于提高飼料利用率，降低養殖成本。相反，以犧牲一定生長速度為代價改進品種的飼料轉化率，可能具有更大的實際意義。

4 結論

本研究發現提高CCKAR基因表達水平更有利于改進雞的飼料轉化率，而g.73206714T＞C的T等位基因是一個與CCKAR基因的高表達活性和高飼料轉化率關聯的分子標記。這些發現為雞的飼料轉化率遺傳調控機制研究提供線索，為從育種和CCKAR活性干預角度改進雞的飼料轉化率提供思路。

[1] REHFELD J F. Cholecystokinin-from local gut hormone to ubiquitous messenger. Frontiers in Endocrinology, 2017, 8: 47. doi:10.3389/fendo. 2017.00047.

[2] DUFRESNE M, SEVA C, FOURMY D. Cholecystokinin and gastrin receptors. Physiological Reviews, 2006, 86(3): 805-847. doi:10.1152/ physrev.00014.2005.

[3] MILLER L J, HOLICKY E L, ULRICH C D, WIEBEN E D. Abnormal processing of the human cholecystokinin receptor gene in association with gallstones and obesity. Journal of the American Chemical Society, 1995, 109(4): 1375-1380. doi:10.1016/0016-5085 (95)90601-0.

[4] Moran T H, Bi S. Hyperphagia and obesity in OLETF rats lacking CCK-1 receptors. Philosophical transactions Royal Society Biological sciences, 2006, 361(1471):1211-1218.

[5] DUNN I C, MEDDLE S L, WILSON P W, WARDLE C A, LAW A S, BISHOP V R, HINDAR C, ROBERTSON G W, BURT D W, ELLISON S J, MORRICE D M, HOCKING P M. Decreased expression of the satiety signal receptor CCKAR is responsible for increased growth and body weight during the domestication of chickens. American Journal of Physiology Endocrinology and Metabolism, 2013, 304(9): E909-E921. doi:10.1152/ajpendo.00580. 2012.

[6] BEGLINGER C, DEGEN L, MATZINGER D, D'AMATO M, DREWE J. Loxiglumide, a CCK-A receptor antagonist, stimulates calorie intake and hunger feelings in humans. American Journal of Physiology Regulatory, Integrative and Comparative Physiology, 2001, 280(4): R1149-R1154. doi:10.1152/ajpregu.2001.280.4.r1149.

[7] EBENEZER I S. Effects of intracerebroventricular administration of the CCK(1) receptor antagonist devazepide on food intake in rats. European Journal of Pharmacology, 2002, 441(1/2): 79-82. doi:10. 1016/s0014-2999(02)01485-1.

[8] CHRISTOFFERSEN B ?, SKYGGEBJERG R B, BUGGE A, KIRK R K, VESTERGAARD B, ULDAM H K, FELS J J, PYKE C, SENSFUSS U, SANFRIDSON A, CLAUSEN T R. Long-acting CCK analogue NN9056 lowers food intake and body weight in obese G?ttingen Minipigs. International Journal of Obesity (2005), 2020, 44(2): 447-456. doi:10.1038/s41366-019-0386-0.

[9] CHOI Y H, FURUSE M, SATOH S, OKUMURA J. Endogenous cholecystokinin is not a major regulator of food intake in the chicken. Journal of Comparative Physiology B, Biochemical, Systemic, and Environmental Physiology, 1994, 164(6): 425-429. doi:10.1007/ bf00714578.

[10] Furuse M, Chol Y H, Yang S I, Kita K, Okumura J. Enhanced release of cholecystokinin in chickens fed diets high in phenylalanine or tyrosine. Comparative Biochemistry and Physiology A Comparative Physiology, 1991, 99(3):449-451.

[11] COVASA M, FORBES J M. Effects of the CCK receptor antagonist MK-329 on food intake in broiler chickens. Frontiers in Genetics, 1994, 48(2): 479-486. doi:10.1016/0091-3057(94)90556-8.

[12] XU Z, JI C, ZHANG Y, ZHANG Z, NIE Q, XU J, ZHANG D, ZHANG X. Combination analysis of genome-wide association and transcriptome sequencing of residual feed intake in quality chickens. BMC Genomics, 2016, 17: 594. doi:10.1186/s12864-016-2861-5.

[13] Dang L P, Liu R F, Zhao W Y, Zhou W X, Min Y N, Wang Z P. Investigating structural impact of a valine to isoleucine substitution on anti-Müllerian hormone in silico and genetic association of the variant and AMH expression with egg production in chickens. Journal of Integrative Agriculture, 2020, 19:1635-1643.

[14] 黨李蘋, 周雯馨, 劉瑞芳, 白云, 王哲鵬. 略陽烏雞體重和產蛋數性狀遺傳參數估計. 中國農業科學, 2020(17): 3620-3628.

DANG L P, ZHOU W X, LIU R F, BAI Y, WANG Z P. Estimation of genetic parameters of body weight and egg number traits of Lueyang black-boned chicken. Scientia Agricultura Sinica, 2020(17): 3620-3628. (in Chinese)

[15] DHANASEKARAN S, RAJ G D, VIGNESH A R, SELVAN S T, PRAKASH B, PERUMAL P, ARIVUDAINAMBI S, BABU T G. Gender identification in Chicken () by PCR using whole blood and dried blood spot on filter paper as template: without prior DNA isolation. bioRxiv, 2016, doi:10.1101/046888.

[16] HONDA K, SANEYASU T, KAMISOYAMA H. Gut hormones and regulation of food intake in birds. The Journal of Poultry Science, 2017, 54(2): 103-110. doi:10.2141/jpsa.0160100.

[17] Wren A M, Bloom S R. Gut hormones and appetite control. Gastroenterology, 2007, 132(6):2116-2130.

[18] ZUIDHOF M J, SCHNEIDER B L, CARNEY V L, KORVER D R, ROBINSON F E. Growth, efficiency, and yield of commercial broilers from 1957, 1978, and 2005. Poultry Science, 2014, 93(12): 2970-2982. doi:10.3382/ps.2014-04291.

[19] VELAZQUEZ-VILLEGAS L A, PERINO A, LEMOS V, ZIETAK M, NOMURA M, POLS T W H, SCHOONJANS K. TGR5 signalling promotes mitochondrial fission and beige remodelling of white adipose tissue. Nature Communications, 2018, 9(1): 245. doi:10.1038/ s41467-017-02068-0.

[20] VIJAYAKUMAR A, YAKAR S, LEROITH D. The intricate role of growth hormone in metabolism. Frontiers in Endocrinology, 2011, 2: 32. doi:10.3389/fendo.2011.00032.

[21] CALOGERO A E, NICOLOSI A M, MONCADA M L, CONIGLIONE F, VICARI E, POLOSA P, D'AGATA R. Effects of cholecystokinin octapeptide on the hypothalamic-pituitary-adrenal axis function and on vasopressin, prolactin and growth hormone release in humans. Neuroendocrinology, 1993, 58(1): 71-76. doi:10. 1159/000126514.

[22] JIN S, YANG L, ZANG H, XU Y, CHEN X, CHEN X, LIU P, GENG Z. Influence of free-range days on growth performance, carcass traits, meat quality, lymphoid organ indices, and blood biochemistry of Wannan Yellow chickens. Poultry Science, 2019, 98(12): 6602-6610. doi:10.3382/ps/pez504.

[23] ZENG T, CHEN L, DU X, LAI S J, HUANG S P, LIU Y L, LU L Z. Association analysis between feed efficiency studies and expression of hypothalamic neuropeptide genes in laying ducks. Animal Genetics, 2016, 47(5): 606-609. doi:10.1111/age.12457.

[24] Owyang C, Heldsinger A. Vagal control of satiety and hormonal regulation of appetite. Journal of neurogastroenterology and motility, 2011, 17:338-348.

[25] PERRY B, WANG Y. Appetite regulation and weight control: the role of gut hormones. Nutrition & Diabetes, 2012, 2: e26. doi:10.1038/ nutd.2011.21.

[26] SMITH C A, ROESZLER K N, HUDSON Q J, SINCLAIR A H. Avian sex determination: what, when and where? Cytogenetic and Genome Research, 2007, 117(1/2/3/4): 165-173. doi:10.1159/ 000103177.

[27] SMITS P, LI P, MANDEL J, ZHANG Z, DENG J M, BEHRINGER R R, DE CROMBRUGGHE B, LEFEBVRE V. The transcription factors L-Sox5 and Sox6 are essential for cartilage formation. Developmental Cell, 2001, 1(2): 277-290. doi:10.1016/s1534-5807(01)00003-x.

[28] SHIKATANI E A, CHANDY M, BESLA R, LI C C, MOMEN A, EL-MOUNAYRI O, ROBBINS C S, HUSAIN M. C-myb regulates proliferation and differentiation of adventitial Sca1+ vascular smooth muscle cell progenitors by transactivation of myocardin. Arteriosclerosis, Thrombosis, and Vascular Biology, 2016, 36(7): 1367-1376. doi:10. 1161/atvbaha.115.307116.

[29] Peng S, Li W, Hou N, Huang N. A review of FoxO1-Regulated metabolic diseases and related drug discoveries. Cells, 2020, 9(1):184.

[30] SAKAGUCHI M, CAI W, WANG C H, CEDERQUIST C T, DAMASIO M, HOMAN E P, BATISTA T, RAMIREZ A K, GUPTA M K, STEGER M, WEWER ALBRECHTSEN N J, SINGH S K, ARAKI E, MANN M, ENERB?CK S, KAHN C R. FoxK1 and FoxK2 in insulin regulation of cellular and mitochondrial metabolism. Nature Communications, 2019, 10(1): 1582. doi:10.1038/s41467- 019-09418-0.

[31] Ren W, Guo J, Jiang F, Lu J, Ding Y, Li A, Liang X, Jia W. CCAAT/enhancer-binding protein α is a crucial regulator of human fat mass and obesity associated gene transcription and expression. BioMed Research International, 2014, 2014:406909.

[32] STAIGER J, LUEBEN M J, BERRIGAN D, MALIK R, PERKINS S N, HURSTING S D, JOHNSON P F. C/EBPbeta regulates body composition, energy balance-related hormones and tumor growth. Carcinogenesis, 2009, 30(5): 832-840. doi:10.1093/carcin/bgn273.

[33] MILLWARD C A, HEANEY J D, SINASAC D S, CHU E C, BEDERMAN I R, GILGE D A, PREVIS S F, CRONIGER C M. Mice with a deletion in the gene for CCAAT/enhancer-binding protein beta are protected against diet-induced obesity. Diabetes, 2007, 56(1): 161-167. doi:10.2337/db06-0310.

[34] Pauwels J, Coopman F, Cools A, Michiels J, Fremaut D, De Smet S, Janssens G P. Selection for growth performance in broiler chickens associates with less diet flexibility. PLoS One, 2015, 10(6):e0127819.

Association of Levels of Cholecystokinin A Receptor Expression and Sequence Variants with Feed Conversion Efficiency of Lueyang Black-Boned Chicken

WANG ZhePeng, ZHOU WenXin, HE JunXi, HU QiaoYan, ZHAO JiaYue

College of Animal Science and Technology, Northwest A&F University, Yangling 712100, Shaanxi

【Background】Increased appetite and feed conversion efficiency forms a physiological basis underlying high growth of modern commercial broiler chickens. Cholecystokinin A receptor (CCKAR) is a key receptor for mediating CCK signals and involving in regulation of appetite and digestion. In view of the important role of CCK signals in the control of energy balance, the levels of CCKAR expression are closely related to the growth of chickens. However, the mechanism by which the abundance of CCKAR affects the growth of chickens remains unclear. 【Objective】 The aim of this study was to elucidate the relationship of levels of CCKAR expression and sequence variants with feed conversion efficiency, which would promote the understanding of genetic basis of the growth traits and breeding for improvement of feed efficiency in the chicken. 【Method】The feed to meat conversion ratio of Lueyang black-boned chickens (rooster n=62, pullet n=67) from 12 to 16 weeks of age and the feed to egg conversion ratio from 26 to 30 weeks of age (n=150) were individually measured. Two synonymous mutations (g.73206714T＞C and g.73209189C＞T) in the CCKAR gene were genotyped by using PCR-RFLP. Association of two SNP with the feed to meat conversion ratio was tested by using two-way analysis of variance with sex treated as a fixed effect. Association of two SNP with the feed to egg conversion ratio was tested by using one-way analysis of variance. Disrupting effect of g.73206714T＞C on transcription factor binding sites was predicted in silico by using JASPAR database. Ten birds that were heterozygotes (CT) at the g.73206714T＞C locus were selected to analyze imbalance of CCKAR expression between T and C alleles. The difference of allele-specific expression was detected by using PCR-RFLP. The Levels of CCKAR expression were detected by qPCR in the hypothalamus, pituitary, pancreas, gallbladder and duodenum of Lueyang black-boned chickens (n=30). The correlation of levels of CCKAR expression with food intake and feed to meat conversion ratio was analyzed. 【Result】 Sex had a significant (=29.44，<0.0001) effect on the feed to meat conversion ratio, and the feed efficiency (5.41±0.76) of roosters was significantly superior to pullets (6.62± 1.00). Of two SNP, only g.73206714T＞C was significantly associated with the feed to meat conversion ratio. The feed to meat conversion ratio (5.18±0.63) of the TT genotype was significantly lower than ones of CC (6.37±1.1) and CT (6.17±0.95) genotypes. The feed to meat conversion ratio had no significant difference between CC and CT. There was no significant interaction effect between genotypes of g.73206714T＞C and sex. G.73206714T＞C potentially disrupted binding sites of 28 transcription factors. Allele-specific expression assay showed that the abundance of transcripts with the T allele was higher than one of transcripts with the C allele, supporting that differential expression of CCKAR gene between alleles might likely contribute to genetic association of g.73206714 with the feed to meat conversion ratio. Average feed to egg conversion ratio of Lueyang black-boned chickens was 2.85±0.41 at the 26-30 weeks of age. Both g.73206714T＞C and g.73209189C＞T had no significant effect on the feed to egg conversion ratio. The levels of CCKAR expression in the pancreas were the highest among five tissues, which showed a negative correlation (=-0.45,=0.017) with feed intake. But there was no significant correlation between levels of CCKAR expression and feed intake in other tissues. In the pancreas (=-0.41,=0.03) and pituitary (=-0.57,=0.0018), the higher levels of CCKAR expression was, the lower the feed to meat conversion ratio was. But in other tissues, there was no the negative correlation.【Conclusion】This study found that increased levels of CCKAR was associated with high feed efficiency. The T allele of g.73296714T＞C was associated with high expression activity of CCKAR gene, and with high feed efficiency.

chicken; growth; CCKAR; feed conversion efficiency; appetite

10.3864/j.issn.0578-1752.2022.22.017

2021-09-03；

2021-11-24

陜西省重點研發計劃（2021NY-028）、略陽烏雞蛋用型配套系選育（WJYJY-2021-3）

王哲鵬，Tel：15619295726；E-mail：wangzhepeng-001@163.com

（責任編輯 林鑒非）