唐鑫鑫,鄭炬梅,駱 娜,營 凡,朱 丹,李 森,劉大偉,安炳星,文 杰, 趙桂蘋*,李和剛

(1.青島農(nóng)業(yè)大學動物科技學院,青島 266109;2.中國農(nóng)業(yè)科學院北京畜牧獸醫(yī)研究所,北京 100193; 3.佛山市高明區(qū)新廣農(nóng)牧有限公司,佛山 528000)

隨著肉雞養(yǎng)殖規(guī)模化和集約化的發(fā)展,商業(yè)化肉雞的選育目標逐漸變?yōu)樯L周期短、增重快、生產(chǎn)性能顯著提高,這給腿部骨骼尚未發(fā)育完全的肉雞帶來健康隱患的同時還導致機體代謝紊亂。這些問題使得我國肉雞腿病發(fā)生率顯著上升,整體可高達2.07%[1],嚴重影響了肉雞養(yǎng)殖的經(jīng)濟效益和動物福利[2-3]。在肉雞腿病中,跛足病約占60%,已被認為是造成禽肉行業(yè)經(jīng)濟損失的主要原因之一[4]。肉雞骨骼腿病發(fā)生原因錯綜復雜,如生長速度、機體吸收代謝、飼養(yǎng)管理方式、光照時間等均會導致腿病發(fā)生的風險。而大多數(shù)影響跛足和腿部畸形的非傳染性原因與骨代謝紊亂有關,伴隨著骨骼變小、骨密度和骨強度降低[5]。有研究發(fā)現(xiàn),雄性肉雞比雌性生長速度更快、腿病的發(fā)生率高出兩倍、行走能力也更差[6],也有研究表明肉雞腿病可遺傳[7-9],其遺傳力在0.10～0.40之間[10-12]。因此,通過遺傳選擇手段能有效降低肉雞腿病發(fā)生率,為肉雞育種工作打下基礎。

隨著基因芯片測序技術的日漸成熟,全基因組關聯(lián)分析(genome-wide association study, GWAS)技術作為一種識別遺傳區(qū)域和性狀間關聯(lián)的無假設方法,可用于鑒定與復雜疾病或性狀相關的單核苷酸多態(tài)性(single nucleotide polymorphism,SNP),挖掘與表型性狀變異相關的基因,已被廣泛使用于人類和畜禽的研究中[13-15]。在人類骨質(zhì)疏松癥的研究中,已經(jīng)通過GWAS發(fā)現(xiàn)了多個與骨骼生物學和代謝相關的基因,其中一些基因位點的突變或多態(tài)性變異與骨質(zhì)疏松癥的發(fā)生和進展密切相關[16-17];Guo等[18]通過對哈伯德肉雞樣本進行雙重基因分型測序(double-digest genotyping by sequencing, ddGBS),鑒定出5個顯著單核苷酸多態(tài)性(SNPs)位點和70個候選基因作為肉雞肢體內(nèi)外翻畸形(valgus-varus deformity, VVD)的潛在標記;Li等[19]將GWAS和選擇特征分析(固定指數(shù)值和核苷酸多樣性比率)相結(jié)合,發(fā)現(xiàn)與骨性狀相關的顯著SNPs主要集中在雞1、4和27號染色體上,確定了21個可能調(diào)控雞骨生長發(fā)育的候選基因。肉雞腿病發(fā)生一直以來是家禽生產(chǎn)中的重大難題,因其發(fā)生原因極其復雜,目前尚不清楚確切的發(fā)生原因和遺傳機制。因此,本研究采用我國自主研發(fā)的“京芯一號”55K SNP基因芯片,對我國自主培育的“廣明2號”白羽肉雞多個群體共2 331只白羽肉雞進行基因組測序,通過GWAS技術鑒定與肉雞腿病發(fā)生相關的遺傳位點和候選基因,闡明肉雞腿病發(fā)生的機理和遺傳基礎,為肉雞腿部疾病的發(fā)病機制提供更為科學合理的解釋。

1 材料與方法

1.1 試驗動物及樣品采集

本次試驗采用我國自主培育的“廣明2號”白羽肉雞A、B、C、D這4個資源群體,各有1 321、711、264、35只不同世代和批次共2 331只公雞個體,所有肉雞均在相同環(huán)境條件下以雙層單籠飼養(yǎng)條件管理,采用自由采食、飲水和一般免疫程序(中國肉雞飼養(yǎng)標準,NY 33-2004)。在肉雞自然生長狀態(tài)下,于42日齡時通過翅間靜脈穿刺采集2 mL血液于EDTA抗凝管中,-20 ℃保存用于后續(xù)DNA測序。

1.2 表型測定

有研究表明,肉雞脛骨彎曲更容易出現(xiàn)畸形和病變,增加跛足甚至腿病發(fā)生的風險[20]。于是本研究通過X光測定肉雞42日齡時腿部健康狀況及骨骼具體形態(tài),將脛骨形變角度≤10°個體作為健康組,定義為0;形變角度>10°個體作為患病組,定義為1,具體測定方法和儀器均參照Zheng等[21]的研究。

1.3 DNA提取與質(zhì)量檢測

本研究采用磁珠法提取肉雞血液基因組DNA。步驟如下:1)使用dsDNA Fragmentase對DNA樣品進行酶切,以實現(xiàn)片段化并修復酶切末端。2)通過2%瓊脂糖凝膠電泳檢測,確認在300～500 bp范圍內(nèi)存在明顯的亮帶,從而選取效果良好的樣品進行后續(xù)反應。3)使用連接酶將測序接頭與片段化DNA連接,并利用磁珠對連接產(chǎn)物進行純化。4)運用PCR對連接產(chǎn)物進行擴增,并使用磁珠對PCR產(chǎn)物進行片段篩選。篩選后的片段產(chǎn)物經(jīng)過Qubit熒光定量儀檢測濃度,2%瓊脂糖凝膠電泳確認片段大小,并通過Qsep400生物分析儀進一步核準文庫片段大小,僅對質(zhì)控合格的樣品進行后續(xù)文庫制備。

1.4 基因分型和質(zhì)量控制

采用中國農(nóng)業(yè)科學院北京畜牧獸醫(yī)研究所研發(fā)的“京芯一號”55K SNP芯片對2 331只肉雞質(zhì)控合格的血液樣品進行基因分型[22],由北京康普森生物技術有限公司進行測定。通過plink(Plink v1.90)軟件對所有個體分型后得到的基因組數(shù)據(jù)進行共有位點合并,再嚴格按照最小等位基因頻率≥5%、基因分型率≥90%的位點以及個體分型率≥90%的標準進行質(zhì)控,剔除不符合要求的個體和位點。在進行全基因組分析之前,位于性染色體(Z染色體和W染色體)中的SNPs被去除。

1.5 群體結(jié)構分析和固定效應檢驗

在建立模型之前,使用主成分分析(principal component analysis, PCA)進行群體分層檢驗,對可解釋方差最大的前2個主成分繪制散點圖。

將腿部健康表型(即健康或患病)作為因變量,世代、批次作為自變量,使用廣義線性回歸模型進行固定效應檢驗,使用二項分布和logit鏈接函數(shù)擬合邏輯回歸模型,并使用卡方檢驗評估世代和批次是否對腿病發(fā)生有顯著影響。本研究發(fā)現(xiàn),肉雞的批次和世代與腿病發(fā)生有顯著的相關性(P<0.05),因此將其作為協(xié)變量加入模型中。

1.6 全基因組關聯(lián)分析

參照此前對肉雞腿部健康表型的分類方式將樣本進行分類,采用plink軟件中的邏輯回歸模型進行GWAS分析。用“—dummy-coding”命令將固定效應轉(zhuǎn)換成0、1格式的二分類變量,將其和前兩個主成分作為協(xié)變量加入邏輯回歸模型中。使用邏輯回歸對腿部健康表型進行GWAS分析建立的模型如下:

Logit(P)=ln(P/(1-P))=β0+β1χ1,β2χ2……,βmχm

其中,P表示性狀發(fā)生的概率,取值范圍為0到1;χ1,χ2, ……,χm表示輸入的m個特征,是模型的自變量;β0,β1,β2,……,βm表示模型參數(shù),需要通過優(yōu)化算法學習得到合適的數(shù)值。

采用plink軟件中的參數(shù)“-indep-Pairwise 25 5 0.2”推測獨立檢驗的有效SNP數(shù)量,以確保SNP單獨具有統(tǒng)計顯著性且互相獨立,從而根據(jù)獨立的SNP數(shù)量調(diào)整全基因組顯著閾值和全基因組水平潛在顯著閾值[23]。計算得到有效獨立檢驗SNPs為7 468個,即GWAS分析的顯著性閾值為0.05/7 468=6.70×10-6,潛在顯著閾值為1/7 468=1.34×10-4。使用R(V4.2.3)軟件包“CMplot”(https://github.com/YinLiLin/CMplot)對GWAS結(jié)果進行曼哈頓圖和Q-Q圖繪制。

1.7 基因注釋及其表型解釋率

為了篩選可能與腿病發(fā)生相關的候選基因,選擇了顯著SNPs基因區(qū)段。利用Ensembl基因組數(shù)據(jù)庫(https://asia.ensembl.org/Sus_scrofa/Info/Index)獲取SNP的信息,對照Gullas(GRCg6a)參考基因組,并對GWAS所得SNPs進行注釋和基因本體(gene ontology,GO)和京都基因與基因組百科全書(kyoto encyclopedia of genes and genomes,KEGG)功能富集分析,以確定蛋白編碼基因的重要生物學功能。為衡量所注釋到的基因?qū)Ρ硇妥儺惖呢暙I程度,通過計算其表型解釋率(phenotypic variation explained,PVE):

beta為基因的效應大小或效應大小的估計值;MAF為最小等位基因頻率(minor allele frequency);se為基因效應的標準誤差(standard error);N為樣本數(shù)量。

識別出對特定表型特征具有較高遺傳影響的基因,以初步篩選影響肉雞腿病發(fā)生的重要候選基因。

2 結(jié) 果

2.1 腿部健康表型評價

X光測定可以精準評價肉雞腿部是處于健康狀況(圖1a)或病變狀態(tài)(圖1b),根據(jù)腿部骨骼形變程度和特征區(qū)分得到2 257只健康和74只腿病個體[24],總體腿病發(fā)生率為3.17%。各群體腿病發(fā)生具體情況統(tǒng)計如表1所示,A、B、C、D群體發(fā)病率分別為4.09%、0.70%、4.55%和8.57%,D群體由于個體數(shù)較少結(jié)果不予考慮,A群體和C群體腿病發(fā)生率均較高,而B群體腿病發(fā)生率較低。

a.健康個體(0)腿部骨骼形態(tài)X光圖像;b.腿病個體(1)腿部骨骼形態(tài)X光圖像 a. X-ray images of leg bone morphology in healthy individuals (0); b. X-ray images of leg bone morphology in leg disease individuals (1)圖1 X光測定評價肉雞腿部健康情況Fig.1 X-ray assessment to evaluate the health condition of broiler legs

表1 試驗群體肉雞腿病發(fā)生情況Table 1 Leg disease of broiler in experimental group

2.2 基因型數(shù)據(jù)分析

利用“京芯一號”55K SNP芯片對2 331只“廣明2號”白羽肉雞進行基因分型,基于最小等位基因頻率≥5%、基因分型率≥90%的位點和個體分型率≥90%的標準對基因型數(shù)據(jù)進行質(zhì)量控制,保留了2 330只肉雞個體和30 414個有效SNPs位點用于后續(xù)分析,各染色體有效位點密度圖見圖2。其中1號染色體上位點數(shù)量最多,有6 601個;16號染色體上位點數(shù)量最少,有134個。

圖2 基因分型質(zhì)控后染色體密度圖Fig.2 Chromosome density map after genotyping quality control

2.3 主成分分析

使用plink軟件對A、B、C、D群體進行主成分分析,結(jié)果(圖3)顯示這4個群體有明顯的結(jié)構分層,說明這4個白羽肉雞群體在人工選育過程中群體結(jié)構產(chǎn)生差異,且PC1和PC2分別可解釋方差的比例為28.4%和6.20%。因此,在后續(xù)全基因組關聯(lián)分析中將前2個主成分作為協(xié)變量加入模型中,以校正群體結(jié)構,降低群體分層對模型穩(wěn)定性的影響。

圖中的每個點對應一個家庭中的個體,紅色、綠色、藍色和紫色的圓點分別代表來自A、B、C、D群體的個體,代表4個群體 Each dot in the figure corresponds to an individual in a family, and the red, green, blue, and purple dots represent individuals from group A, B, C, and D, respectively, representing the 4 groups圖3 A、B、C、D群體結(jié)構主成分分析Fig.3 Principal component analysis of the structure of A, B, C and D groups

2.4 基于腿病表型的全基因組關聯(lián)分析

本研究基于肉雞腿部健康表型,將群體分為對照和患病群體進行GWAS分析,以獲得與腿病相關的位點和候選基因,結(jié)果用曼哈頓圖(圖4a)和QQ圖表示(圖4b)。QQ圖膨脹因子為1.018,表示分析結(jié)果的假陽性較低。曼哈頓圖顯示沒有與腿病發(fā)生相關的顯著SNP位點,在潛在顯著水平上共檢測到5個SNPs位點(P<9.90×10-5),分別位于6號染色體(rs7317511、rs7963346、rs11152911、rs11277299)和18號染色體(rs3239706)。在6號染色體上,最顯著的SNP是rs11277299(P=1.23×10-5)。

2.5 候選基因的SNP表型解釋率和基因頻率

所鑒定到的顯著SNP位點解釋的表型變異解釋(phenotypic variation explanation, PVE)范圍為0.65%～0.81%,這5個位點PVE值較其他位點高,說明對表型的貢獻較大。rs11277299位點的表型變異率最高,為0.81%,認為這幾個SNPs位點可能是重要的遺傳位點,在肉雞腿部健康發(fā)生發(fā)展調(diào)控中發(fā)揮重要作用。對上述5個顯著SNPs位點進行注釋,共得到TBCD、SIRT1和PBLD這3個候選基因(表2),通過計算各SNP位點在健康組(0)和患病組(1)間的比例,發(fā)現(xiàn)SIRT1基因中攜帶GG基因型的個體、PBLD基因中攜帶AA和GA基因型的個體、TBCD基因中攜帶CC基因型的個體患病比例較大,均有可能增加肉雞患腿病的風險(圖5)。

2.6 候選基因功能研究

對腿病表型注釋到的TBCD、SIRT1和PBLD這3個候選基因通過KOBAS網(wǎng)站(http://kobas.cbi.pku.edu.cn)進行GO和KEGG功能富集分析,結(jié)果顯示共有82個GO項達到了具有統(tǒng)計學意義的水平(P<0.05),這些GO條目主要涉及細胞分化、細胞生物合成、細胞因子結(jié)合和信號通路負調(diào)控等過程,轉(zhuǎn)化生長因子受體信號通路的負向調(diào)控GO條目(GO: 0030512)是GO富集分析結(jié)果中最顯著的GO條目,其中細胞質(zhì)過程(GO:0005737)在這3個基因中均顯著富集到(P<0.05)(圖6)。KEGG分析結(jié)果表明,僅有SIRT1基因富集到FoxO信號通路(gga04068)和細胞衰老信號通路(gga04218)并達到顯著水平(P<0.05),可能在腿部疾病和骨骼發(fā)育中起重要作用。

a.曼哈頓圖:水平黑色虛線表示全基因組潛在閾值線(P=1.34×10-4),黑色實線表示全基因組顯著閾值線(P=6.70×10-6);b.Q-Q圖 a. Manhattan plot: Horizontal black dashed line indicates the genome-wide potential threshold line (P=1.34×10-4), and the black solid line indicates the genome-wide significant threshold line (P=6.70×10-6); b; Q-Q plot圖4 全基因組水平上腿部健康表型性狀的曼哈頓圖和Q-Q圖Fig.4 Manhattan and Q-Q plots of leg health phenotypic traits at the genome-wide level

表2 白羽肉雞腿部健康性狀顯著相關的SNPs及其候選基因Table 2 SNPs and their candidate genes significantly associated with leg health traits in white feather broilers

a, b, c分別為顯著SNPs rs7317511、rs11277299和rs3239706的基因型頻率在健康組(0)和腿病組(1)間的比例,即不同等位基因?qū)ν炔拷】敌誀畹挠绊?a, b, c are the ratios of genotype frequencies of significant SNPs rs7317511, rs11277299, and rs3239706, respectively, between the healthy group (0) and the leg disease group (1), i.e., the effect of different alleles on leg health traits圖5 注釋到基因的顯著SNPs在健康和腿病組間的基因型比例Fig.5 Genotype proportions of significant SNPs annotated to genes between healthy and leg disease groups

圖6 候選基因GO富集分析Fig.6 GO enrichment analysis of candidate genes

3 討 論

隨著對現(xiàn)代商業(yè)化肉雞增重和生產(chǎn)性能的不斷選育,骨骼生長發(fā)育不足以支撐肉雞快速生長的體重導致多種類型腿病發(fā)生,腿病發(fā)生已引起養(yǎng)殖產(chǎn)業(yè)極大關注。本研究采用我國自主培育的“廣明2號”白羽肉雞,通過X光方法對已在人工選育過程中群體結(jié)構產(chǎn)生差異的A、B、C、D四個群體測定腿部骨骼健康狀況,對實踐應用有重要的指導意義。腿病個體在品系培育過程中經(jīng)歷了嚴格淘汰,可能導致已上市的配套系中相應性狀遺傳力降低,而本研究中特定的通過X光方法定義表型也可能對遺傳力估計造成影響。需要繼續(xù)探究腿病個體淘汰對整個品系性狀的影響且分析造成遺傳力較低的原因,進一步完善通過X光方法評價肉雞腿病發(fā)生的狀況,為改進肉雞品系培育和預防腿病發(fā)生提供借鑒。

本研究基于多個品系、批次和世代群體腿病發(fā)病率的統(tǒng)計,發(fā)現(xiàn)腿病發(fā)生率總體高達3.17%,A、C群體發(fā)病率分別高達4.09%、4.55%,這極度加大了家禽生產(chǎn)中的死淘率,造成了嚴重的經(jīng)濟損失,但群體內(nèi)世代和批次間的腿病發(fā)生率無顯著差異。而B群體腿病發(fā)生率較低,可進一步探究B群體腿部健康狀況遺傳機制,從遺傳選育角度入手改良腿部健康性狀。肉雞腿病發(fā)生受遺傳、營養(yǎng)、飼養(yǎng)管理、環(huán)境等多方面的影響,其中定位其主效遺傳因和位點是作為突破口進行遺傳改良以改善肉雞腿部疾病的關鍵[25]。

全基因組關聯(lián)分析(GWAS)作為揭示基因與疾病之間關聯(lián)性的有力研究方法,能為疾病的機制研究和治療提供新的方向和策略,已被廣泛應用于挖掘復雜疾病的易感基因[26-27]。在畜禽生產(chǎn)中,GWAS方法已被廣泛應用于經(jīng)濟性狀的研究中,如體重、產(chǎn)蛋性狀、免疫功能和骨骼等[28]。通過全基因組關聯(lián)研究和基因組選擇等手段能迅速加快腿部健康的遺傳改良進程,然而目前關于肉雞腿部健康表型相關遺傳變異位點的研究仍然十分缺乏,沒有大樣本用于鑒定與腿病相關的基因位點和候選基因。因此,本研究基于白羽肉雞腿部疾病發(fā)生的病例和對照樣本,將2 257例健康和74例腿病發(fā)生白羽肉雞通過“京芯一號”55K SNP芯片基因分型測序方法進行全基因組關聯(lián)分析,以期鑒定出與白羽肉雞腿部健康相關的候選基因和變異位點,為肉雞腿部健康遺傳改良提供參考。

本研究通過全基因組關聯(lián)分析定位到5個潛在SNPs位點與腿病發(fā)生顯著相關,其中4個位點位于6號染色體,1個位點位于18號染色體。PVE值可以作為評估基因位點對表型變異解釋能力的指標,本研究發(fā)現(xiàn)在肉雞腿部健康表型變異中,位點rs7317511、rs7963346、rs11152911、rs11277299和rs3239706的基因突變或多態(tài)性分別解釋了0.65%、0.66%、0.72%、0.81%和0.72%表型的變異,這些位點較其他位點PVE值大,一定程度上能作為評估肉雞腿部健康表型的變異位點的指標。其中,rs7317511注釋得到SIRT1基因,rs3239706注釋得到TBCD基因,rs11277299注釋得到PBLD基因,這3個已知候選基因相對來說對腿部健康表型具有較為重要的潛在意義。

骨骼作為雞的一個重要性狀,它的質(zhì)量與肉雞的產(chǎn)肉量和骨質(zhì)疏松癥有關。有研究采用高密度基因分型平臺檢測骨骼性狀的候選基因,發(fā)現(xiàn)位于不同染色體上的3個基因位點包含候選基因HTR2A、LPAR6、CAB39L、TRPC4、WNT9A、SPOP、NGFR、GIP和HOXB3與骨骼質(zhì)量和長度有關[29];Zhang等[30]對雞龍骨表型差異進行選擇消除分析,發(fā)現(xiàn)了10個具有強選擇信號且與骨骼相關的重要候選基因并富集到了MAPK和破骨細胞分化兩條重要的骨信號通路。重要的是,基于本研究注釋到的基因SIRT1、TBCD和PBLD具體功能及其如何參與機體骨骼代謝過程也已有多項研究。雞SIRT1基因cDNA全長約2 271 bp,位于6號染色體,可編碼756個氨基酸[31],它參與調(diào)節(jié)腫瘤、氧化應激、凋亡、衰老和炎癥等一系列細胞代謝過程[32-33]。有研究發(fā)現(xiàn),SIRT1敲除小鼠體內(nèi)骨量減少,成骨細胞與破骨細胞比例失衡,這表明該基因是通過直接作用于骨細胞來調(diào)控骨代謝過程[34],它的表達水平與骨礦物質(zhì)密度和骨脆性有關,其可能成為系統(tǒng)性骨代謝和骨相關疾病的潛在生物標志物[35],還能介導FoxO(gga04068)和細胞衰老(gga04218)等重要信號通路。FoxO是一類轉(zhuǎn)錄因子,參與調(diào)控多種細胞過程,而FoxO信號通路通過參與調(diào)控骨母細胞的增殖分化和吸收等過程影響骨細胞的功能和骨骼組織的形成[36-37]。細胞衰老信號通路主要包括端?？s短、DNA損傷響應、染色質(zhì)重塑、細胞周期調(diào)控和細胞死亡等多個分子機制。有研究發(fā)現(xiàn)了細胞衰老信號通路與骨骼生長發(fā)育過程的關聯(lián),即端?？s短和端粒酶的活性與骨質(zhì)疏松和骨折風險增加相關[38],但目前對于細胞衰老信號通路在骨骼生長發(fā)育中的具體功能還需要進一步研究證實。關于SIRT1基因?qū)趋来x和生長發(fā)育的調(diào)控機制和互作網(wǎng)絡在未來的研究中值得進一步關注,特別是GG基因型的肉雞。PBLD通過抑制NF-κB和上皮到間充質(zhì)轉(zhuǎn)化信號通路而發(fā)揮腫瘤抑制作用,控制炎癥的關鍵調(diào)節(jié),對骨關節(jié)炎等疾病都有著重要影響[39];TBCD主要參與微管相關的生物學過程,如細胞分裂、細胞軸向極化和細胞運輸?shù)萚40],尚未發(fā)現(xiàn)該基因與骨骼之間的聯(lián)系,但微管在骨骼生長和骨骼細胞的某些功能中也發(fā)揮著重要作用[41]。關于這3個基因如何具體影響骨骼生長發(fā)育的生理過程及其分子遺傳機理,以及攜帶不同基因型的肉雞腿部健康是否會有差異,還亟待進一步的理論和試驗研究進行探索。

本研究還發(fā)現(xiàn)了特定的潛在SNP位點rs7963346,其位置被定位到了lncRNA上,這個發(fā)現(xiàn)提供了進一步研究該SNP位點對腿病發(fā)生影響的線索。lncRNA是一類長鏈非編碼RNA,它的功能相對較為復雜和多樣,可以調(diào)控基因表達、細胞生物學過程等。可通過進一步試驗深入研究該位點所在的lncRNA的功能,還可通過構建SNP-lncRNA-disease網(wǎng)絡來探究該位點與腿病之間的關系,從而揭示其在腿部疾病發(fā)生發(fā)展機制中的作用,為肉雞腿病的預防、診斷和治療提供重要的理論依據(jù)。大多數(shù)疾病發(fā)生發(fā)展都是極其復雜的,受多種環(huán)境條件和基因控制。本研究在全基因組水平上發(fā)現(xiàn)多個顯著與肉雞腿病發(fā)生相關的位點,但這些位點未定位到相似的區(qū)域和基因。骨骼生長發(fā)育及病變過程受遺傳和環(huán)境因素共同影響調(diào)控,基因組的復雜性和多樣性決定了不同環(huán)境背景下基因的功能可能會呈現(xiàn)不同的效應。推測肉雞腿病作為復雜性狀,受到了多種效應的影響,進一步的功能注釋、環(huán)境相關研究以及更大規(guī)模的驗證研究將有助于深入理解這些位點在肉雞腿部健康中的潛在作用和復雜關系。SIRT1、PBLD和TBCD均是可能與腿病發(fā)生相關的候選基因,有望深入研究這些基因的潛在生物學意義,并為肉雞腿部健康研究提供新的見解。

4 結(jié) 論

本研究基于多個品系及世代的2 331只白羽肉雞個體,通過大量樣本數(shù)據(jù)基于55K SNP基因芯片利用GWAS方法研究白羽肉雞腿部健康的遺傳機制,挖掘與腿病發(fā)生相關的重要遺傳信息。研究發(fā)現(xiàn)6號染色體(rs7317511、rs7963346、rs11152911、rs11277299)和18號染色體(rs3239706)上的5個SNPs位點可能是影響肉雞腿病發(fā)生的變異位點,SIRT1、PBLD和TBCD可能是影響白羽肉雞腿病發(fā)生的重要候選基因。本試驗為研究肉雞腿病的發(fā)生機制和降低腿病發(fā)生率提供了重要遺傳參考。