大規(guī)模群體解析豬日增重及達(dá)百千克體重日齡的潛在因果基因

摘 要： 旨在探究影響豬達(dá)100 kg體重日齡（age at 100 kg，AGE）和達(dá)100 kg體重平均日增重（average daily gain at 100 kg，ADG）的候選基因。本研究采集了來自大白、長白和杜洛克3個品種的共計4 593頭健康成年豬的耳組織作為試驗材料，其中公豬2 563頭，母豬2 030頭。通過記錄豬只的日齡和體重，計算校正AGE和ADG。通過酚氯仿法提取樣品DNA，利用“中芯一號”50K SNP芯片對樣本進(jìn)行基因分型，并對質(zhì)控后的SNPs進(jìn)行基因型填充，將SNPs位點數(shù)從4萬填充至800萬。隨后，利用GEMMA混合線性模型對AGE和ADG性狀進(jìn)行填充前后的全基因組關(guān)聯(lián)分析，使用BEDTools在顯著位點上、下游1 Mb范圍查找候選基因。同時，結(jié)合PigGTEx中的34個組織的表達(dá)量數(shù)量性狀位點數(shù)據(jù)，使用R軟件進(jìn)行共定位分析，挖掘與GWAS信號共享同一因果變異的基因，通過GWAS和共定位分析，確定AGE和ADG性狀的候選基因。GWAS分析結(jié)果顯示，AGE與ADG性狀的顯著SNP為1號染色體上的1_270827213。在該位點上、下游1 Mb范圍內(nèi)，共鑒定到32個基因。共定位分析結(jié)果顯示，對于AGE性狀，有10個基因的eQTL信號與GWAS信號共定位。對于ADG性狀，有11個基因的eQTL信號與GWAS信號共定位。最終，確定了CRAT、GPR107和USP20這3個基因作為AGE的候選基因，確定了CRAT、GPR107、USP20、FNBP1、PTGES和HMCN2這6個基因作為ADG的候選基因。本研究為豬品種改良提供了分子標(biāo)記，對豬生長相關(guān)性狀功能基因挖掘奠定基礎(chǔ)。

關(guān)鍵詞： 豬；AGE；ADG；全基因組關(guān)聯(lián)分析；候選基因

中圖分類號：

S828.2"""" 文獻(xiàn)標(biāo)志碼：A"" ""文章編號： 0366-6964（2025）03-1100-10

收稿日期：2024-07-26

基金項目：國家重點研發(fā)項目資金（2021YFD1200801）

作者簡介：黃雅妮（2000-），女，江西上饒人，博士生，主要從事動物遺傳育種與繁殖研究，E-mail：yani720@163.com

*通信作者：肖石軍，主要從事地方豬新品種培育和種豬育種新技術(shù)推廣應(yīng)用，E-mail：shjx_jxau@hotmail.com；張志燕，主要從事數(shù)量遺傳學(xué)研究，E-mail：bioducklily@hotmail.com

Large-scale Population Analysis of Potential Causal Genes for Daily Weight Gain and

Age at 100 kg in Pigs

HUANG" Yani1， TANG" Xi1， LI" Jingquan1， WEI" Jiacheng1， WU" Zhenfang2， LI" Xinyun3，

XIAO" Shijun^1*， ZHANG" Zhiyan^1*

（1.National Key Laboratory for Swine Genetic Improvement and Germplasm Innovation，

Jiangxi Agricultural University， Nanchang 330045，" China; 2.College of Animal Science，

South China Agricultural University， Guangzhou 510642，" China; 3.School of Animal

Science and Technology， Huazhong Agricultural University， Wuhan 430070，" China）

Abstract：" The aim of this study was to explore candidate genes influencing age at 100 kg （AGE） and average daily weight gain at 100 kg （ADG） in pigs. Ear tissue samples were collected from 4 593 healthy adult pigs， including 2 563 boars and 2 030 sows， representing three breeds： Large White， Landrace， and Duroc. The age and body weight of the pigs were recorded， and the corrected AGE and ADG were calculated accordingly. The DNA from the samples was extracted using the phenol-chloroform method. Genotyping was performed using the \"CC-1\" 50K SNP chip， and SNP markers were imputed from 40 000 to 8 million after quality control. Subsequently， genome-wide association study was performed using the GEMMA mixed linear model on AGE and ADG traits. Candidate genes were identified by searching within a 1 Mb window upstream and downstream of significant SNP loci using BEDTools. Additionally， the study integrated expression quantitative trait locus data from 34 tissues in the PigGTEx， and colocalization analysis was conducted using R software to identify genes that shared causal variants with the GWAS signals. Through GWAS and colocalization analysis， candidate genes for AGE and ADG traits were identified. The significant SNP associated with both AGE and ADG traits is 1_270827213 on chromosome 1. Within a 1 Mb region upstream and downstream of this SNP， 32 candidate genes were identified. Colocalization analysis showed that， there were 10 genes had eQTL signals that colocalized with the GWAS signal for AGE， while 11 genes had eQTL signals colocalizing with the GWAS signal for ADG. Ultimately， the genes CRAT， GPR107 and USP20 were identified as candidate genes for AGE， and CRAT， GPR107， USP20， FNBP1， PTGES and HMCN2 were identified as candidate genes for ADG. This study provides molecular markers for pig breed improvement and lays the foundation for the functional gene discovery related to growth traits in pigs.

Keywords： pigs; AGE; ADG; genome-wide association analysis; candidate genes

*Corresponding authors： XIAO Shijun， E-mail：shjx_jxau@hotmail.com; ZHANG Zhiyan， E-mail：bioducklily@hotmail.com

2023年，我國全年生豬出欄量已達(dá)到72 662萬頭^［1］，年消費豬肉超過了5 800萬噸^［1］，是世界上最大的生豬生產(chǎn)和消費國。生豬養(yǎng)殖對國民生活和消費具有重大影響。在生豬養(yǎng)殖中，豬的生長性能對生豬生產(chǎn)具有重要的經(jīng)濟(jì)效應(yīng)，其中，100 kg體重日齡（age at 100 kg，AGE）和100 kg體重平均日增重（average daily gain at 100 kg，ADG）是衡量豬生長速率的重要經(jīng)濟(jì)性狀，對養(yǎng)殖場的生產(chǎn)效率以及經(jīng)濟(jì)效益至關(guān)重要。西方商品豬的AGE平均值為154～167 d，ADG為610～820 g·d^-1［2-4］，這兩個生長性狀通常受到遺傳與非遺傳因素的影響，包括豬的品種、飼養(yǎng)方式以及營養(yǎng)水平等。由于ADG與AGE均具有中等遺傳力^［5-8］，表明可以通過遺傳方法進(jìn)行改善。

全基因組關(guān)聯(lián)研究（genome-wide association study，GWAS）是一種鑒別與性狀和疾病相關(guān)遺傳變異的有效方法。隨著測序技術(shù)的進(jìn)步，GWAS加快了對重要性狀和疾病的遺傳解析步伐^［9］。自2005年以來，GWAS廣泛應(yīng)用于數(shù)量性狀位點（quantitative trait locus，QTL）的鑒定和復(fù)雜性狀候選基因的定位。此外，隨著高密度單核苷酸多態(tài)性（single nucleotide polymorphism，SNP）芯片的發(fā)展和分析成本的降低，GWAS在家畜中得到了廣泛的應(yīng)用。然而GWAS發(fā)現(xiàn)的許多關(guān)聯(lián)位點在編碼區(qū)域之外，了解這些關(guān)聯(lián)的生物學(xué)基礎(chǔ)的方法之一就是探索細(xì)胞中間表型（特別是基因表達(dá)量）的GWAS信號與生物體復(fù)雜性狀的GWAS信號是否位于同一位點（即共定位），這將有助于這些性狀的生物學(xué)解釋^［10］。同時，當(dāng)前公共數(shù)據(jù)庫中已有大量關(guān)于基因表達(dá)量的GWAS匯總數(shù)據(jù)，即表達(dá)數(shù)量性狀位點（expression quantitative trait locus，eQTL）數(shù)據(jù)。通過共定位分析可以探索eQTL與GWAS信號是否由相同的因果變異驅(qū)動，從而加強基因與表型之間的關(guān)聯(lián)性。

本研究利用“中芯一號”50K芯片對大白、長白和杜洛克3個商業(yè)品種共計4 593頭豬進(jìn)行基因分型，針對AGE和ADG兩個重要的生長性狀開展GWAS研究，結(jié)合顯著位點附近的基因以及GWAS與eQTL的共定位分析來探究AGE與ADG性狀的候選基因。本研究鑒定到的候選位點和候選基因?qū)ωi基因組選育和生豬生產(chǎn)提供重要的理論基礎(chǔ)。

1 材料與方法

1.1 樣品及數(shù)據(jù)來源

本研究的所有個體均來自加大種豬有限公司，共4 593頭健康成年豬（公豬2 563頭、母豬2 030頭），包含3個品種，即大白豬（2 576頭）、長白豬（1 298頭）和杜洛克豬（719頭）。采集測定個體耳組織，使用酚氯仿法對每個個體的DNA進(jìn)行提取，將質(zhì)控后合格的DNA樣本利用“中芯一號”50K芯片進(jìn)行基因型分型，共獲得57 466個SNPs位點。

按照種豬生產(chǎn)性能測定規(guī)程（NYT 822—2019）要求對表型進(jìn)行測定并記錄。AGE校正公式如下：

AGE=S1+（WT-W1）×S1-AW1

式中，AGE為校正達(dá)100 kg日齡（d）；S1為出生到結(jié)測當(dāng)天的自然天數(shù)（d）；WT為目標(biāo)體重，本研究中該值為100（kg）；W1為結(jié)測當(dāng)天稱量的實際體重（kg）；A為校正參數(shù)。

AGE校正參數(shù)見表1。

ADG100校正方法如下：

按上式計算達(dá)30 kg體重日齡，按GB/T 8170對計算結(jié)果進(jìn)行修約：

AGE30=S+［30-W］×B

式中，AGE30為達(dá)到30 kg體重日齡（d）；S 為從出生到開測當(dāng)天的自然天數(shù)（d）；30為設(shè)定的開測體重（kg）；W為開測當(dāng)天稱量的實際重量（kg）；B為校正參數(shù)（杜洛克為1.536，長白為1.565，大白為1.550）。

按照上式計算測定日期日增重：

ADG=70×1 000（AGE-AGE30）

式中，ADG為測定期日增重（g）；70為目標(biāo)體重（100 kg）與開測體重（30 kg）之差（kg）；1 000為計量單位由千克換算為克。

1.2 數(shù)據(jù)處理

為了提高全基因組關(guān)聯(lián)分析的準(zhǔn)確性，本研究對表型數(shù)據(jù)和基因型數(shù)據(jù)均進(jìn)行了處理。針對AGE和ADG的表型數(shù)據(jù)，去除偏離數(shù)據(jù)3個標(biāo)準(zhǔn)差以上表型對應(yīng)的個體，最終獲得4 501頭豬的AGE數(shù)據(jù)和ADG數(shù)據(jù)，隨后，使用R（v4.4.0）軟件建立線性回歸模型對品種和性別進(jìn)行校正。使用 PLINK^［11］（v.1.90b7.2）對基因型數(shù)據(jù)進(jìn)行質(zhì)控，去除次等位基因頻率小于0.01的SNPs位點、基因型缺失率大于0.01的SNPs位點、沒有染色體位置信息的SNPs位點，以及性染色體上的SNPs位點。最終保留44 766個高質(zhì)量SNPs用于GWAS研究以及基因型填充。

1.3 基因型填充

對質(zhì)控后的基因型數(shù)據(jù)，使用Shapeit5^［12］（v5.1.1）構(gòu)建單倍型，以提高填充準(zhǔn)確性。隨后使用Beagle^［13］軟件以本實驗室Tong等^［14］構(gòu)建的千豬單倍型為參考面板對芯片位點進(jìn)行基因型填充。使用Plink再次對填充后的SNPs位點進(jìn)行質(zhì)控，去除實際等位基因與最可能等位基因相關(guān)性小于0.85的SNPs位點，去除次等位基因頻率小于0.01的SNPs、缺失率大于0.01的SNPs，剩余SNPs用于后續(xù)分析。使用R（v.x64 4.1.1）中的CMplot包（v.4.0）（https：//github.com/YinLiLin/CMplot）繪制SNPs密度圖。

1.4 全基因組關(guān)聯(lián)分析

使用PLINK對填充后的SNPs進(jìn)行主成分分析（principal component analysis，PCA），并在R中進(jìn)行可視化，觀察群體分層情況。隨后，使用GEMMA^［15］（v0.98.3）的混合線性模型對兩個性狀分別進(jìn)行GWAS分析，本研究確定顯著關(guān)聯(lián)的閾值為5×10^-8［16］。使用CMplot包生成了曼哈頓圖和Q-Q圖，使用PLINK估計Top SNP周圍區(qū)域中SNP對之間的連鎖不平衡（linkage disequilibrium，LD），并通過LDBlockshow（v1.40）^［17］進(jìn)行可視化。使用非參數(shù)檢驗法Kruskal-Wallis test對顯著位點的基因型和表型進(jìn)行單因素方差分析，以確定顯著位點不同基因型對應(yīng)的表型是否存在顯著差異。

1.5 候選基因鑒定

使用LocusZooms-like（v2.1）（https：//github.com/Geeketics/LocusZooms）對Top SNP上、下游200 kb范圍的基因進(jìn)行可視化。同時，使用BEDTools^［18］（v2.31.1）工具提取GWAS顯著位點上、下游1 Mb范圍的基因作為候選基因進(jìn)行進(jìn)一步分析，基因注釋文件使用Ensemble網(wǎng)站（https：//asia.ensembl.org/index.html）中的豬Sus scrofa 11.1基因組版本。然后在DAVID網(wǎng)站^［19］對候選基因進(jìn)行Gene Ontology （GO）分析，并在R中可視化。

1.6 共定位分析

本研究使用貝葉斯共定位方法來鑒定特定區(qū)域內(nèi)eQTL與GWAS信號之間的的共享因果變異，使用R包coloc^［20］（v5.2.3）中的coloc.abf函數(shù)對兩個性狀分別進(jìn)行GWAS與eQTL共定位分析，其中eQTL數(shù)據(jù)來自PigGTEx^［21］發(fā)布的34種組織的Cis-eQTL數(shù)據(jù)集，包括脂肪、肝臟、肌肉等。共定位區(qū)域根據(jù)顯著位點所在的1號染色體上的基因進(jìn)行確定，在eQTL數(shù)據(jù)集中篩選每個基因的Top SNP位點，提取Top SNP位點上下游1 Mb的區(qū)域與GWAS進(jìn)行共定位分析。該分析強調(diào)了5種后驗概率（PP），本研究認(rèn)為PP.H4gt;0.75，表明GWAS與eQTL之間存在共同的因果變異^［22］。

2 結(jié) 果

2.1 表型描述性統(tǒng)計分析

對AGE和ADG的表型數(shù)據(jù)進(jìn)行描述性統(tǒng)計（表2），AGE和ADG的表型平均值為160.07 d和615.70 g·d^-1，這一結(jié)果與前人研究報道相一致。

2.2 SNP在基因組分布及PCA分析

本研究對4 593頭豬進(jìn)行單倍型構(gòu)建以及基因型填充，填充后SNPs位點由44 766個增加到8 206 050個（圖1A），填充準(zhǔn)確性高達(dá)0.99，表明基因型填充策略準(zhǔn)確高效，可用于進(jìn)一步的關(guān)聯(lián)分析。由于本研究采用的群體來自3個不同的品種，將填充后的SNPs用于PCA分析，提取前兩個主成分進(jìn)行繪圖，PCA圖顯示3個品種能夠很好地分離（圖1B）。

2.3 AGE與ADG的全基因組關(guān)聯(lián)分析

對填充前和填充后的SNPs分別進(jìn)行AGE及ADG性狀的全基因組關(guān)聯(lián)分析。其中，兩個性狀填充前的GWAS結(jié)果均沒有顯著關(guān)聯(lián)位點（圖2A）。填充后GWAS結(jié)果顯示兩個性狀均有一個共同的顯著SNP位點1_270827213（圖2B），該位點在AGE性狀GWAS中的P值為1.23×10^-8，效應(yīng)值為-1.246，在ADG性狀GWAS中的P值為5.29×10^-9，效應(yīng)值為4.839（表3），表明該SNP有助于縮短達(dá)100kg體重日齡，增加100kg體重平均日增重。LDBlockshow圖顯示該顯著位點與附近的SNP位點存在強連鎖關(guān)系（圖2C）。接著對該位點上、下游200 kb范圍基因進(jìn)行可視化，LocusZooms-like圖顯示Top SNP位點落在ABL1基因上，該基因與其附近的ASS1基因以及EXCOS2基因連鎖程度較高（圖3A）。提取顯著位點1_270827213的基因型，與表型進(jìn)行單因素方差分析。結(jié)果顯示，在兩個性狀中，該位點不同基因型對應(yīng)的表型存在顯著差異（Plt;0.05，圖3B）。使用BEDToos工具獲得顯著位點上下游1 Mb范圍的基因列表，共檢測到32個基因，包括ABL1、AIF1L、ASS1、EXOSC2、FAM78A、FIBCD1、FUBP3、HMCN2、LAMC3、NCS1、NUP214、PLPP7、PRDM12、QRFP、NTMT1、ASB6、PRRX2、TOR1B、PTGES、FNBP1、PRRC2B、P0MT1、UCK1、RAPGEF1、MED27、USP20、C9orf78、PRRT1B、TOR1A、GPR107、NCS1和SNORD62。對這32個基因進(jìn)行GO富集分析（圖3C），結(jié)果表明這些基因富集到ERK1 and ERK2 cascade（ERK1和ERK2激酶級聯(lián)反應(yīng)）以及negative regulation of ERK1 and ERK2 cascade（ERK1和ERK2激酶級聯(lián)反應(yīng)負(fù)調(diào)控）條目，該級聯(lián)反應(yīng)參與多種過程的調(diào)節(jié)，包括細(xì)胞黏附、細(xì)胞周期進(jìn)程、細(xì)胞遷移、細(xì)胞存活、分化、代謝、增殖和轉(zhuǎn)錄。此外，還和kinesin binding（肌動蛋白結(jié)合）以及Substrate adhesion-dependent cell spreading（底物黏附依賴性細(xì)胞擴(kuò)展）條目有關(guān)，這些生物學(xué)過程對于遷移、增殖、分化以及維持組織結(jié)構(gòu)至關(guān)重要。

2.4 共定位分析確定潛在靶基因

為了確定風(fēng)險變異的潛在靶基因，本研究使用coloc R包進(jìn)行了兩個性狀填充后GWAS與eQTL共定位分析。AGE性狀的GWAS與34個組織eQTL分析結(jié)果顯示，有10個基因的PP.H4gt;0.75（表4），且共定位基因主要在血液、大腦和肌肉組織中。ADG性狀的GWAS與34個組織eQTL共定位分析結(jié)果顯示，有11個基因的PP.H4gt;0.75（表5），其中5個基因存在于肌肉組織中，而日增重性狀與肌肉組織的生長發(fā)育密切相關(guān)。6個基因（CRAT1、TOR1A、GPR107、EXOSC2、ENDOG和USP20）在兩個性狀的共定位分析中均顯示出顯著的共定位。

2.5 確定潛在候選基因

將顯著變異位點上、下游1 Mb范圍鑒定到的基因與共定位分析鑒定到的基因取交集，確定了4個AGE的候選基因（CRAT、EXOSC2、GPR107和USP20），以及8個ADG的候選基因（CRAT、EXOSC2、GPR107、USP20、FNBP1、PTGES、HMCN2和ASS1）以查找其生物學(xué)功能。

3 討 論

3.1 GWAS分析確定AGE以及ADG的顯著SNP位點及候選基因

本研究對4 593頭豬使用“中芯一號”50K芯片進(jìn)行基因分型，經(jīng)過表型校正處理、基因型質(zhì)控、基因型填充、PCA分析等一系列分析，使用混合線性模型對大規(guī)模群體豬的AGE和ADG性狀進(jìn)行全基因組關(guān)聯(lián)分析。GWAS結(jié)果顯示，兩個性狀均存在一個相同的顯著SNP位點1_270827213，該位點在前人的研究中并未發(fā)現(xiàn)^{［3，5，23-24］}。對該位點附近的基因進(jìn)行定位，檢測到32個候選基因，將這32個基因作為第一條證據(jù)鏈。接下來，利用PigGTEx發(fā)布的34個組織的eQTL數(shù)據(jù)集，進(jìn)行eQTL與GWAS共定位分析，將PP.H4gt;0.75的基因作為第二條證據(jù)鏈，結(jié)合兩條證據(jù)鏈確定AGE和ADG性狀的候選基因。

3.2 本研究鑒定的候選基因與前人研究結(jié)果的比較

目前，已有多項研究嘗試對ADG和AGE兩個性狀進(jìn)行遺傳解析，然而大多數(shù)研究使用有限數(shù)量的豬^［23］（樣本量lt;1 000）或者未填充的芯片數(shù)據(jù)^［24-25］進(jìn)行GWAS研究。竇騰飛等^［24］對1 805大白豬的50K芯片分型結(jié)果進(jìn)行AGE性狀的GWAS分析，確定了8個候選基因（DUSP6、MDF1、MYF5、MYF6、CSRP2、TAF8、FOXP4和ALX1），Ruan等^［5］對3 945頭加系杜洛克豬的50K芯片分型結(jié)果進(jìn)行了ADG和AGE兩個性狀的加權(quán)單步GWAS，并報道了6個ADG和AGE的候選基因（FAM135B、SLC27A6、ADRB2、ZFAT、NFIL3和ROR2）。Zhou等^［25］對美系杜洛克和加系杜洛克豬的50K芯片分型結(jié)果進(jìn)行單群體GWAS以及meta GWAS分析，并報道了兩個ADG的候選基因，STC1（美系杜洛克）和PHLPP1（加系杜洛克）。Jung等^［26］使用貝葉斯GWAS對韓國商業(yè)豬（杜洛克、長白和約克夏）進(jìn)行ADG和90 kg體重日齡（AGE90）的GWAS研究，鑒定了ADG重要候選基因有FAM177B（長白）、MGAT4C（約克夏）、BMP2（約克夏）、TRIB2（杜洛克）、RBMS3（杜洛克）等，以及AGE90的重要候選基因有LACTBL1（杜洛克）、HAO1（約克夏）、FAM177B（長白）等。Tang等^［3］對4個品種豬（杜洛克、約克夏、長白和皮特蘭）統(tǒng)一進(jìn)行ADG和AGE性狀的GWAS研究，其中9個候選基因被鑒定為與生長和代謝有關(guān)（TRIB3、GRP、TTR、CNR1、GLP1R、BRD2、HCRTR2、SEC11C和ssc-mir-122）。由此可見，不同研究方法和遺傳背景下鑒定到的候選基因差異較大。在本研究中并未發(fā)現(xiàn)與以上研究重合的候選基因，原因可能是分析方法和品種品系的不同，但本研究關(guān)聯(lián)到的顯著SNP位點和候選基因仍然可以為后續(xù)研究提供參考，且本研究中鑒定到的HMCN2在另一項豬的生長性狀GWAS研究中也"" 被確定為ADG性狀的候選基因^［27］。

3.3 重要的候選基因生物學(xué)功能

CRAT基因作為肉堿?；D(zhuǎn)移酶家族的一員，位于線粒體基質(zhì)中，參與能量穩(wěn)態(tài)和脂肪代謝^［28］。它通過將乙酰輔酶A轉(zhuǎn)化為乙酰肉堿酯，調(diào)節(jié)線粒體和細(xì)胞內(nèi)碳運輸^［29］，從而影響豬的生長發(fā)育。

HMCN2基因?qū)儆诩?xì)胞外基質(zhì)蛋白的纖維蛋白家族，在組織黏附和細(xì)胞遷移中發(fā)揮重要作用^［30］。在一項針對豬生長性狀的GWAS研究中，也將該基因確定為豬ADG性狀的候選基因^［27］，暗示其對生長性狀的影響。PTGES基因是前列腺素E2的誘導(dǎo)型合成酶^［31］，參與脂肪酸代謝和新陳代謝。該基因在脂肪前細(xì)胞中表達(dá)，并且在脂肪生成過程中始終檢測到其mRNA和蛋白質(zhì)水平，可能通過調(diào)節(jié)脂肪生成進(jìn)而對豬的生長產(chǎn)生影響。USP20基因是泛素特異性肽酶，穩(wěn)定膽固醇生物合成途徑中的限速酶HMGCR。其缺失或抑制可顯著減少飲食引起的體重增加，降低血清和肝中的脂質(zhì)水平，提高胰島素敏感性并增加能量消耗^［32］，表明該基因也可能通過調(diào)節(jié)脂質(zhì)水平進(jìn)而調(diào)整體重變化。FNBP1基因是一種肌動蛋白骨架相關(guān)蛋白，對許多基本的細(xì)胞過程至關(guān)重要，例如細(xì)胞形態(tài)，運動和胞質(zhì)分裂，在內(nèi)吞過程中與肌動蛋白細(xì)胞骨架重組協(xié)調(diào)膜變形^［33］。這一機制可能影響細(xì)胞的分裂進(jìn)而影響豬的生長相關(guān)性狀。GPR107是G蛋白偶聯(lián)受體，廣泛表達(dá)于多種組織。它與神經(jīng)抑素相互作用^［34］，在代謝功能、激素分泌等方面發(fā)揮重要作用^［35］，從而可能通過多方面的調(diào)節(jié)影響豬的生長發(fā)育。

綜上所述，以上6個基因在不同的代謝途徑和細(xì)胞過程中的作用，共同決定了豬的生長發(fā)育特性。然而，接下來還需要進(jìn)一步增加樣本數(shù)量提高發(fā)現(xiàn)的能力，并對候選基因開展生物學(xué)驗證，確定這些基因?qū)ιL性狀的實際影響。

4 結(jié) 論

本研究基于4 593頭豬的基因型數(shù)據(jù)，采用50K芯片進(jìn)行基因分型，對AGE和ADG兩個性狀進(jìn)行GWAS研究，發(fā)現(xiàn)了一個在AGE和ADG中均顯著的SNP位點1_270827213。在顯著位點附近1 Mb范圍檢測到32個候選基因，進(jìn)一步通過GWAS與eQTL共定位分析，確定了3個AGE的候選基因（CRAT、GPR107、USP20），6個ADG的候選基因（CRAT、GPR107、USP20、FNBP1、PTGES、HMCN2）。本研究中新檢測到的顯著SNP以及候選基因可應(yīng)用于豬的選擇育種，以加快豬的生長速度，進(jìn)而提高經(jīng)濟(jì)效益。

參考文獻(xiàn)（References）：

［1］ 張海峰，王祖力，陳澤芳，等.2023年全球生豬產(chǎn)業(yè)發(fā)展情況及2024年的趨勢［J］.豬業(yè)科學(xué)，2024，41（2）：32-36.

ZHANG H F，WANG Z L，CHEN Z F，et al.The global pig industry development in 2023 and trends for 2024［J］.Swine Industry Science，2024，41（2）：32-36.（in Chinese）

［2］ JIANG Y，TANG S，WANG C，et al.A genome-wide association study of growth and fatness traits in two pig populations with different genetic backgrounds［J］.J Anim Sci，2018，96（3）：806-816.

［3］ TANG Z S，XU J Y，YIN L L，et al.Genome-wide association study reveals candidate genes for growth relevant traits in pigs［J］.Front Genet，2019，10：302.

［4］ HERRERA-CCERES W，SNCHEZ J P.Selection for feed efficiency using the social effects animal model in growing Duroc pigs：evaluation by simulation［J］.Genet Sel Evol，2020，52（1）：53.

［5］ RUAN D L，ZHUANG Z W，DING R R，et al.Weighted single-step GWAS identified candidate genes associated with growth traits in a duroc pig population［J］.Genes （Basel），2021，12（1）：117.

［6］ 許 迪，顏 港，張 帥，等.大白豬生長性狀影響因素分析及遺傳參數(shù)估計［J］.中國畜牧獸醫(yī)，2024，51（1）：193-202.

XU D，YAN G，ZHANG S，et al.Influencing factors analysis and genetic parameters estimation of growth traits in yorkshire pigs［J］.China Animal Husbandry amp; Veterinary Medicine，2024，51（1）：193-202.（in Chinese）

［7］ 張笑科，廖偉莉，陳信佑，等.杜洛克豬生長性狀全基因組關(guān)聯(lián)分析及候選基因鑒定［J］.畜牧獸醫(yī)學(xué)報，2023， 54（5）：1868-1876.

ZHANG X K，LIAO W L，CHEN X Y，et al.Genome-wide association study for identifying candidate genes of growth traits in Duroc pigs［J］.Acta Veterinaria et Zootechnica Sinica，2023，54（5）：1868-1876.（in Chinese）

［8］ 賀婕妤，王斌虎，廖 柱，等.長白和大白豬主要生長性狀的遺傳參數(shù)估計［J］.畜牧獸醫(yī)學(xué)報，2021，52（8）：2115-2123.

HE J Y，WANG B H，LIAO Z，et al.Estimation of genetic parameters of main growth traits in landrace and large white pigs［J］.Acta Veterinaria et Zootechnica Sinica，2021，52（8）：2115-2123.（in Chinese）

［9］ LOOS R J F，YEO G S H.The genetics of obesity：from discovery to biology［J］.Nat Rev Genet，2022，23（2）：120-133.

［10］ NICA A C，MONTGOMERY S B，DIMAS A S，et al.Candidate causal regulatory effects by integration of expression QTLs with complex trait genetic associations［J］.PLoS Genet，2010，6（4）：e1000895.

［11］ PURCELL S，NEALE B，TODD-BROWN K，et al.PLINK：a tool set for whole-genome association and population-based linkage analyses［J］.Am J Hum Genet，2007，81（3）：559-575.

［12］ HOFMEISTER R J，RIBEIRO D M，RUBINACCI S，et al.Accurate rare variant phasing of whole-genome and whole-exome sequencing data in the UK Biobank［J］.Nat Genet，2023，55（7）：1243-1249.

［13］ BROWNING B L，ZHOU Y，BROWNING S R.A one-penny imputed genome from next-generation reference panels［J］.Am J Hum Genet，2018，103（3）：338-348.

［14］ TONG X K，CHEN D，HU J C，et al.Accurate haplotype construction and detection of selection signatures enabled by high quality pig genome sequences［J］.Nat Commun，2023，14（1）：5126.

［15］ ZHOU X，STEPHENS M.Genome-wide efficient mixed-model analysis for association studies［J］.Nat Genet，2012，44（7）： 821-824.

［16］ JANNOT A S，EHRET G，PERNEGER T.Plt;5×10^-8 has emerged as a standard of statistical significance for genome-wide association studies［J］.J Clin Epidemiol，2015，68（4）：460-465.

［17］ DONG S S，HE W M，JI J J，et al.LDBlockShow：a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files［J］.Brief Bioinform，2021，22（4）：bbaa227.

［18］ QUINLAN A R，HALL I M.BEDTools：a flexible suite of utilities for comparing genomic features［J］.Bioinformatics， 2010，26（6）：841-842.

［19］ SHERMAN B T，HAO M，QIU J，et al.DAVID：a web server for functional enrichment analysis and functional annotation of gene lists （2021 update）［J］.Nucleic Acids Res，2022，50（W1）：W216-W221.

［20］ GIAMBARTOLOMEI C，VUKCEVIC D，SCHADT E E，et al.Bayesian Test for colocalisation between Pairs of genetic association studies using summary statistics［J］.PLoS Genet，2014，10（5）：e1004383.

［21］ TENG J Y，GAO Y H，YIN H W，et al.A compendium of genetic regulatory effects across pig tissues［J］.Nat Genet， 2024，56（1）：112-123.

［22］ GUI J X，YANG X Y，TAN C，et al.A cross-tissue transcriptome-wide association study reveals novel susceptibility genes for migraine［J］.J Headache Pain，2024，25（1）：94.

［23］ 鄒若男，時坤鵬，張志勇，等.杜洛克豬生長性狀全基因組關(guān)聯(lián)分析［J/OL］.中國畜牧雜志， （2024-04-18）.https：//doi.org/10.19556/j.0258-7033.20230629-08.

ZOU R N，SHI K P，ZHANG Z Y，et al.Genome-wide association analysis of growth traits in Duroc pigs［J/OL］.Chinese Journal of Animal Science，（2024-04-18）.https：//doi.org/10.19556/j.0258-7033.20230629-08.（in Chinese）

［24］ 竇騰飛，吳姿儀，白利瑤，等.全基因組關(guān)聯(lián)分析鑒定大白豬生長性狀遺傳變異及候選基因［J］.中國畜牧雜志，2023，59（8）：264-272.

DOU T F，WU Z Y，BAI L Y，et al.Genome-wide association analysis identifies genetic variations and candidate genes for growth traits in Large White pigs［J］.Chinese Journal of Animal Science，2023，59（8）：264-272.（in Chinese）

［25］ ZHOU S P，DING R R，MENG F M，et al.A meta-analysis of genome-wide association studies for average daily gain and lean meat percentage in two Duroc pig populations［J］.BMC Genomics，2021，22（1）：12.

［26］ JUNG J H，LEE S M，OH S H.A genome-wide association study on growth traits of Korean commercial pig breeds using Bayesian methods［J］.Anim Biosci，2024，37（5）：807-816.

［27］ 王文斌.基于低深度重測序挖掘豬經(jīng)濟(jì)性狀相關(guān)突變位點［D］.武漢：華中農(nóng)業(yè)大學(xué)，2023.

WANG W B.Mining of mutation sites related to economic traits of pigs based on low depth resequencing［D］.Wuhan： Huazhong Agricultural University，2023.（in Chinese）

［28］ VIOLANTE S，IJLST L，RUITER J，et al.Substrate specificity of human carnitine acetyltransferase：implications for fatty acid and branched-chain amino acid metabolism［J］.Biochim Biophys Acta，2013，1832（6）：773-779.

［29］ MUOIO D M，NOLAND R C，KOVALIK J P，et al.Muscle-specific deletion of carnitine acetyltransferase compromises glucose tolerance and metabolic flexibility［J］.Cell Metab，2012，15（5）：764-777.

［30］ FEITOSA N M，ZHANG J L，CARNEY T J，et al.Hemicentin 2 and Fibulin 1 are required for epidermal-dermal junction formation and fin mesenchymal cell migration during zebrafish development［J］.Dev Biol，2012，369（2）：235-248.

［31］ SAMPEY A V，MONRAD S，CROFFORD L J.Microsomal prostaglandin E synthase-1：the inducible synthase for prostaglandin E2［J］.Arthritis Res Ther，2005，7（3）：114-117.

［32］ LU X Y，SHI X J，HU A，et al.Feeding induces cholesterol biosynthesis via the mTORC1-USP20-HMGCR axis［J］.Nature， 2020， 588（7838）：479-484.

［33］ TSUJITA K，SUETSUGU S，SASAKI N，et al.Coordination between the actin cytoskeleton and membrane deformation by a novel membrane tubulation domain of PCH proteins is involved in endocytosis［J］.J Cell Biol，2006，172（2）：269-279.

［34］ YOSTEN G L C，REDLINGER L J，SAMSON W K.Evidence for an interaction of neuronostatin with the orphan G protein-coupled receptor，GPR107［J］.Am J Physiol Regul Integr Comp Physiol，2012，303（9）：R941-R949.

［35］ SAMSON W K，ZHANG J V，AVSIAN-KRETCHMER O，et al.Neuronostatin encoded by the somatostatin gene regulates neuronal，cardiovascular，and metabolic functions［J］.J Biol Chem，2008，283（46）：31949-31959.

（編輯 郭云雁）