張會永,李國輝,薛 倩,周成浩,殷建玫,蘇一軍,夏樹立,韓 威
(1.江蘇省家禽科學研究所,揚州 225125;2.江蘇省家禽科學研究所科技創(chuàng)新有限公司,揚州 225125;3.天津市農業(yè)科學院畜牧獸醫(yī)研究所,天津 300381)
大圍山微型雞(Daweishan Mini chicken,WX)俗稱“香雞”、“金雞”、“嬌雞”,原產地為云南省屏邊縣,屬于肉蛋兼用和競技觀賞型地方品種,2011年收錄于《中國畜禽遺傳資源志·家禽志》[1]。在長期的自然選擇和人工選擇過程中,大圍山微型雞形成了獨具一格的外貌特征和種質特性,具有體型小、胸腿肌發(fā)達、靈活好斗、耐粗飼、抗病性強等特點,是中國珍稀的家禽遺傳資源。
針對大圍山微型雞的研究主要集中在表型、微衛(wèi)星技術及線粒體DNA(mtDNA)的研究。李正田等[2]對310日齡不同性別大圍山微型雞的屠宰性能進行了測定分析,為品種的選育、開發(fā)利用提供了基礎資料。佟薈全等[3]對大圍山微型雞生長曲線分析與擬合比較分析中發(fā)現,大圍山微型雞生長曲線拐點為6.94周齡。賈俊靜等[4]利用微衛(wèi)星標記法對大圍山微型雞遺傳多樣性進行分析,結果表明其平均雜合度和群體平均多態(tài)信息含量較低。賈曉旭等[5]基于mtDNA D-loop區(qū)全序列變異對大圍山微型雞遺傳多樣性及其起源進化關系進行研究,結果表明其遺傳多樣性較低,大圍山微型雞在遺傳來源組成上具有4個血統。歐陽依娜等[6]對大圍山微型雞、云龍矮腳雞和拉伯高腳雞mtDNA D-loop遺傳多樣性分析的研究發(fā)現,大圍山微型雞由7個母系血統組成。
隨著分子技術的不斷發(fā)展,簡化基因組測序(RAD-Seq)技術以高通量、成本低、操作簡便、靈敏度高等特點[7],一次運行可以獲得0.45~200 Gb的數據,克服了微衛(wèi)星標記及mtDNA遺傳信息覆蓋整個雞基因組范圍極其微小的缺點(FAO推薦用于雞遺傳多樣性分析的微衛(wèi)星標記約30個,線粒體DNA的D-loop區(qū)域全長也僅約1 200 bp[8]),本試驗利用RAD-Seq技術對大圍山微型雞及其他20個品種雞進行測序分析,從分子水平上探討大圍山微型雞的遺傳進化,為品種保護、評價和開發(fā)利用提供科學依據。
大圍山微型雞保種群體與其他18個地方雞品種,包括北京油雞(BY)、東鄉(xiāng)綠殼蛋雞(DX)、鹿苑雞(LY)、河南斗雞(DJ)、瓢雞(PJ)、蕭山雞(XS)、文昌雞(WC)、邊雞(BJ)、狼山雞(LS)、白耳黃雞(BE)、茶花雞(CH)、仙居雞(XJ)、安義瓦灰雞(WH)、固始雞(GS)、大骨雞(DG)、惠陽胡須雞(HX)、金湖烏鳳雞(JH)、藏雞(ZZ)及2個國外引進品種隱性白羽肉雞(RW)、安卡雞(AK),每個品種按照家系選擇30個樣本(公雞10只、母雞20只),個體間無親緣關系。在同一飼養(yǎng)條件下,300日齡時翅下靜脈采血1.5 mL/只,檸檬酸鈉(ACD)抗凝,-20 ℃保存?zhèn)溆谩?/p>
1.2.1 DNA提取與質量檢測 采用常規(guī)酚-氯仿法提取21個品種雞血液樣本基因組DNA,-20 ℃保存?zhèn)溆?。使用紫外分光光度計對提取的DNA進行濃度和純度的測定。
1.2.2 構建pair-end文庫 檢測合格的DNA樣品,采用ddRAD建庫方式構建長度范圍在300~500 bp的pair-end文庫:取基因組DNA 500 ng,加入0.6 UEcoRⅠ(NEB公司)、T4DNA連接酶(NEB公司)、ATP(NEB公司)和EcoRⅠ接頭(含區(qū)分樣品的Index序列)在37 ℃反應3 h,65 ℃退火1 h。然后加限制性內切酶NlaⅢ(NEB公司)和NlaⅢ接頭在37 ℃反應3 h。反應結束后在PCR儀中65 ℃放置30 min使內切酶失活。使用1.2%瓊脂糖凝膠電泳選擇400~600 bp回收酶切產物。使用Qubit3.0(Life Technology公司)對回收產物進行DNA定量,等量混合樣品。 使用IlluminaTruSeq試劑盒對混合產物進行DNA文庫構建。
1.2.3 RAD-Seq簡化基因組測序與SNP質控 在pair-end文庫構建完成后,在IlluminaHiSeq X平臺進行RAD-Seq簡化基因組測序。SNP的檢測使用GATK4.1.2.0和samtoolsl.9軟件工具包,質控條件:堿基測序正確識別率(QPhred (20),Q20)>95%;SNP深度>60%,SNP檢出率>70%,最小等位基因頻率(minorallele frequency,MAF)>0.05。
應用PopGen32軟件進行平均雜合度(heterozygosity,Ho)、核苷酸多樣度(Pi)、近交系數(Fis)和群體分化指數(Fst)計算。應用Admixture1.3.0軟件進行群體聚類分析,采用PLINK1.90軟件進行選擇信號檢測,基于Fst檢驗方法進行選擇信號分析(以100 kb為窗口、10 kb為步長),Z/W性染色體位點不進行分析。 利用DAVID在線軟件(https:∥david.ncifcrf.gov/)對受選擇基因進行GO功能注釋和KEGG信號通路分析。
通過數據質控,在大圍山微型雞群體中鑒定出331 892個SNPs,主要分布于1、2、3、4、5號和Z染色體上(圖1),SNP位點數量與染色體長度總體上呈正相關。
圖1 SNP在染色體上的分布Fig.1 Distribution of SNP on chromosomes
21個品種的平均雜合度、核苷酸多樣度和近交系數見表1,大圍山微型雞的平均雜合度為0.219;核苷酸多樣度為0.245,與其他品種相比,呈中度多態(tài);近交系數為0.107。
表1 各品種雞的平均雜合度、核苷酸多樣度和近交系數Table 1 Ho,Pi and Fis of chicken breeds
大圍山微型雞與其他18個地方雞品種的Fst見表2,大圍山微型雞與瓢雞的Fst最低(0.0929),親緣關系最近;與河南斗雞的Fst最高(0.2179),親緣關系最遠。對21個雞種進行遺傳結構聚類分析發(fā)現,大圍山微型雞與瓢雞、茶花雞和藏雞聚合在一起(圖2),親緣關系較近。
表2 大圍山微型雞與其他品種雞間的群體分化指數Table 2 Fst between WX and other chicken breeds
圖2 大圍山微型雞及其他品種雞的遺傳進化樹Fig.2 Phylogenetic tree of WX and other chicken breeds
本研究基于Fst檢驗方法,大圍山微型雞分別與其他18個地方雞品種比較,以100 kb為窗口、10 kb為step進行滑動的Fst選擇信號分析,篩選Fst值前5%的位點作為選擇位點,注釋相應的基因,以≥10個品種受到選擇的基因作為選取標準,共篩選出200個受到選擇的基因。GO功能分析結果顯示,這些受選擇基因主要富集在對刺激的應答、信號傳導、運動等生物過程(圖3)。KEGG通路富集分析結果表明,這些受選擇基因主要富集在新陳代謝、肌動蛋白細胞骨架的調節(jié)、RIG-Ⅰ樣受體、Notch、MAPK、Apelin等信號通路(圖4)。 根據GO功能和KEGG通路富集分析結果,對與大圍山微型雞種質特性相關基因進行初步篩選發(fā)現,NTRK2、SEMA3A、NTF3、NRG1基因與刺激反應信號傳導等生物過程相關;NTRK2、GNAQ、GNA14、CRHBP基因與攻擊行為相關;PDE4D、SEMA3A、AGGF1基因與心血管系統相關;ALDH1A1、PDE8B、ECH1、PTPRD基因與代謝相關;SEMA3A、NPR2基因與骨骼發(fā)育相關,TULP3、FOXP2基因與顱骨形成相關;PRPF3、GLI3、DIAPH3基因與視力、聽覺感知相關;TRIM14、MALT1基因與免疫相關。
圖3 受選擇基因的GO功能分析Fig.3 GO function analysis of selected genes
圖4 受選擇基因的KEGG分析Fig.4 KEGG analysis of selected genes
開展地方雞保種群遺傳多樣性分析對保種工作至關重要,保種工作的核心工作之一就是確保種群世代傳遞過程中遺傳多樣性的穩(wěn)態(tài)。賈靜雯等[4]用微衛(wèi)星方法測得的大圍山微型雞平均雜合度為0.1737,賈曉旭等[5]基于mtDNA D-loop區(qū)分析大圍山微型雞的核苷酸多樣度為0.00443,上述研究說明大圍山微型雞的遺傳多樣性相對匱乏。本研究結果顯示,大圍山微型雞的平均雜合度處于中等水平,核苷酸多樣度處于中度多態(tài)。與前人研究不一致的原因可能是本研究對大圍山微型雞群體的采樣更有代表性,囊括了其主要類型(白羽、白麻、黑麻、黃麻、黑花等,單冠、豆冠),采樣個體表型上有一定差異。
大圍山微型雞的近交系數與其他18個地方雞品種相比,處于中上水平,這提示在今后的保種工作中,應嚴格控制近交水平的過快增加,保持遺傳多樣性的相對穩(wěn)定,必要時可從原產地引進一定數量的群體(或其他類群)進行血緣更新。群體遺傳結構分析中與大圍山微型雞親緣關系較近的3個地方雞品種,從地理位置上看,大圍山微型雞與瓢雞、茶花雞同為云南省地方雞品種,藏雞在云南省西北部迪慶藏族自治州也有分布,具有明顯的地域特征。
本研究基于Fst檢驗方法對大圍山微型雞種質特性相關基因進行篩選,共篩選出200個受選擇的基因,進行了GO功能和KEGG通路分析。結果表明,與其他18個地方雞品種及2個引入肉雞品種相比,受選擇候選基因顯著富集于不同的生物學通路。
自然選擇是進化的原動力,大圍山微型雞公雞長期處于半野生半家養(yǎng)狀態(tài),為了保護種群應對野外的危險,因而具有較強的戰(zhàn)斗力,與其他種群公雞相比,它體型較小,成年公雞的平均體重僅為1 kg左右,但斗性強。本研究發(fā)現,大圍山微型雞有著發(fā)達的神經系統,與刺激反應信號傳導等生物學過程相關的基因NTRK2、SEMA3A、NTF3、NRG1等受到了選擇。研究表明,NTRK2編碼酪氨酸激酶受體B (TrkB),不僅可促進神經細胞生存,增加突觸可塑性及神經發(fā)生[9-10];SEMA3A是神經軸突發(fā)育的一個關鍵信號蛋白,通過結合神經纖毛蛋白1(NRP1)、神經纖毛蛋白2(NRP2)和叢狀蛋白A(LXNA)復合物受體參與軸突排斥﹑樹突分支、突觸形成和神經元遷移過程[11];NTF3是早期神經發(fā)生期間神經元發(fā)育的關鍵介質,在維持突觸功能和突觸形成中起重要作用,在神經營養(yǎng)蛋白選擇性地結合特定的酪氨酸受體激酶存儲信號在各種信號通路中的作用,對神經元的存活和軸突投射至關重要[12-14];NRG1影響神經膠質細胞的發(fā)育與遷移、突觸可塑性、髓鞘形成[15-17]。這些基因使得大圍山微型雞反應迅速、行動敏捷。
5-羥色胺(5-HT)、多巴胺(DA)已被證實是調控攻擊行為的主要神經遞質[18],本研究發(fā)現,與此相關的基因NTRK2、GNAQ、GNA14、CRHBP受到了選擇。研究表明,NTRK2對5-羥色胺和多巴胺能神經元的生長發(fā)育、分化再生具有維持和促進作用[19],GNAQ、GNA14通過磷脂酶C激活多巴胺受體信號通路[20-21]。李佳玉等[22]研究發(fā)現,CRHBP基因rs10062367位點多態(tài)性可能與攻擊行為有關。本研究中還發(fā)現與心血管系統相關的基因PDE4D、SEMA3A、AGGF1等受到了選擇。PDE4D能夠調解心肌細胞收縮和心率[23],SEMA3A能夠降低室性心速過快[24],AGGF1在誘導血管生成、血管發(fā)育以及血管損傷后新內膜形成發(fā)揮作用[25]。發(fā)達的心血管系統和充足的能量供應為大圍山微型雞打斗的爆發(fā)力和持久力、野外逃避危險的能力提供了保障。多個代謝通路(糖酵解和糖異生、甘油磷脂代謝、甘油酯代謝等)相關的基因ALDH1A1、PDE8B、ECH1、PTPRD受到了選擇,這可能為大微山微型雞打斗爆發(fā)力和持久力提供充足的能量。
大圍山微型雞骨細,而較強的斗性行為要求必須有堅硬的骨骼,本研究發(fā)現與骨骼發(fā)育相關的基因SEMA3A、NPR2等、與顱骨形成相關的基因TULP3、FOXP2,以及與視力和聽覺感知相關的基因PRPF3、GLI3、DIAPH3等也受到了選擇。研究表明,SEMA3A可促進成骨細胞和抑制破骨細胞的生成,并可調節(jié)骨重塑[26-27];NPR2能夠編碼利納肽B型受體(NPR-B),NPR-B與C型利納肽結合促進骨骼基質的合成和刺激軟骨的增殖分化[28]。研究發(fā)現,FOXP2基因能夠調控顱骨形成和長骨骨重塑[29];TULP3是正常胚胎發(fā)育所必需的,它的破壞導致胚胎顱面區(qū)域形態(tài)缺陷[30],這些基因可能使得大圍山微型雞形成了獨特的頭部外貌特征。PRPF3在視網膜的發(fā)育過程中起重要作用[31]。GLI3可通過調節(jié) Wnt/β-catenin信號傳導來控制晶狀體和視網膜的發(fā)育[32]。這些基因可能使大圍山微型雞能夠面對野外生存危險保持足夠的機警。在免疫調節(jié)方面,TRIM14和MALT1等基因受到了選擇。TRIM14能夠正向調控抗病毒天然免疫[33],MALT1蛋白可以調控免疫細胞的生存和增殖分化,在神經系統損傷等多種疾病發(fā)揮作用[34]。這些基因可能與大圍山微型雞的抗病能力有關。本研究初步篩選出大圍山微型雞的特性及其相關基因,這些基因所在的通路及其作用還需要進一步的研究。
大圍山微型雞遺傳多樣性呈中度多態(tài),與瓢雞親緣關系最近(Fst為0.0929),并篩選出了200個受到選擇的基因,這些基因在代謝、信號傳導、顱骨發(fā)育等多個方面發(fā)揮作用。結果為大圍山微型雞種質特性的深入挖掘、保護與開發(fā)利用提供了科學依據。