中國地方品種黃雞線粒體DNA D-loop遺傳多樣性研究

黃勛和，翁茁先，李威娜，王慶，何丹林，羅威，張細權，杜炳旺

中國地方品種黃雞線粒體DNA D-loop遺傳多樣性研究

1嘉應學院/廣東省山區特色農業資源保護與精準利用重點實驗室/廣東省五華三黃雞科技創新中心，廣東梅州 514015；2華南農業大學動物科學學院，廣州 510642

【目的】地方品種黃雞是中國重要的家禽遺傳資源，系統評估其遺傳多樣性水平，為提高利用率和制定科學的保護策略提供依據。【方法】對新測的694份樣本和下載的589份樣本，合計28個中國南方地方品種黃雞共1 283份的線粒體DNA片段（D-loop，519 bp）遺傳多樣性進行整合分析，構建單倍型中介網絡圖，通過單倍型地理分布格局、主坐標分析、分子變異分析和中性檢測，確認其內部的群體遺傳結構和群體歷史。通過分析亞洲和太平洋地方雞線粒體DNA地理分布格局，推測中國南方地方品種黃雞稀有單倍型類群的來源?！窘Y果】從1 283份樣本檢測到101個變異位點，其中92個為多態位點。定義了169種單倍型，歸屬于6個單倍型類群A-E和G，其中A-C和E為優勢單倍型類群，D和G為稀有單倍型類群，占總體樣本比例分別為15.43%、49.26%、18.55%、16.37%、0.31%和0.08%。河南省單倍型類群最豐富（6個）；廣東省單倍型數量最多（61），浙江?。?9）和海南省（12）較少。單倍型類群A和B在28個黃雞品種中均有分布；D只分布于淮南麻黃雞、固始雞、寧都黃雞和霞煙雞中；G只分布于固始雞中。從單倍型類群D的分布頻率和單倍型數量來看，中國南方地方品種黃雞單倍型類群D可能來源于東南亞地區。從單倍型類群G的地理分布和中介網絡圖來看，河南和南亞的單倍型類群G可能來源于中國西南地區。中國南方地方品種黃雞總體單倍型多樣性和核苷酸多樣性分別為0.903±0.005和0.01269±n.d.。河南省和湖南省的單倍型多樣性和核苷酸多樣性最高，分別為0.916±0.011、0.01358±0.00039和0.913±0.012、0.01345±0.00042；海南省的單倍型多樣性最低（0.736±0.076），江西省的核苷酸多樣性最低（0.00981±0.00072）。在不同地方品種間，淮南麻黃雞、固始雞和鄖陽大雞的單倍型多樣性最高，江漢雞、淮南麻黃雞和黃郎雞的核苷酸多樣性最高，洪山雞、廣西三黃雞和蕭山雞的單倍型多樣性和核苷酸多樣性最低。文昌雞與河田雞的遺傳關系最近，與歷史記載相吻合。地方品種黃雞群體遺傳分化不明顯，分子變異主要來源于品種內。中性檢測顯示黃雞在品種水平上近期未經歷明顯種群擴張，但單倍型類群A、B和E經歷明顯的擴張歷史。【結論】結果表明地方品種黃雞總體上保留著較高的遺傳多樣性水平，處于較好的保種狀態，但洪山雞、蕭山雞和廣西三黃雞應加強保護。地方品種黃雞存在雜交現象，整體上經歷了群體擴張歷史。東南亞和西南地方雞對中國南方地方品種黃雞有一定的遺傳貢獻。

黃雞；線粒體DNA D-loop；遺傳變異；遺傳分化；群體歷史

0 引言

【研究意義】中國獨特的地理環境和多樣的風俗文化形成了眾多形態各異的地方雞品種，僅在《中國畜禽遺傳資源志：家禽志》記載的就有107個，滿足了不同地區、不同消費習慣及不同消費層次的需求[1]。地方品種黃雞是指產于當地、歷史悠久、品種純正、具有“三黃”（毛黃、腳黃、皮黃）或黃麻羽的美麗外觀、獨特的肉質風味、廣受消費者歡迎的地方品種，早在公元540年前北魏時期的農學專著《齊民要術》上就有記載[1-2]。地方品種黃雞是祖先留給我們的寶貴財富，是我國重要的畜禽遺傳資源和極具產業優勢的特色肉雞品種，也是我國畜牧業可持續發展的種源基礎。在當今世界各國把種質資源挖掘保護作為國家戰略高度之際，誰占有資源，誰就占領了制高點[1]?？茖W認識和正確評估遺傳多樣性水平是今后制定家禽遺傳資源保護規劃和開展利用工作的重要基礎和依據。【前人研究進展】線粒體DNA（mitochondrial DNA，mtDNA）因母系遺傳、核苷酸突變速率快、無重組等特點，被廣泛用于動物遺傳多樣性和群體母系歷史研究[3-4]。Wu等[5]分析了中國14個地方雞品種的mtDNA控制區（D-loop，片段長度360 bp）遺傳多樣性，結果顯示這些品種處于較高的遺傳變異水平，不存在明顯的遺傳結構。基于黃郎雞mtDNA D-loop（片段長度524 bp）的研究顯示其處于較高的遺傳變異水平，起源于中國北方和西南地區，同時受鄰省地方雞的影響[6]。賈曉旭等[7]評估了華東地區11個地方雞品種（含白耳黃雞、崇仁麻雞、仙居雞、瑯琊雞等黃雞）的mtDNA D-loop全序列（片段長度1 232 bp）核苷酸多樣性，探討了品種的母系起源，指出部分品種可能受到歐美高產品系的遺傳滲入?；趍tDNA D-loop全序列（片段長度1 232 bp）研究表明江蘇省鹿苑雞、溧陽雞和太湖雞等3個黃雞的遺傳多樣性較低，主要單倍型類群為A和B[8]。黃勛和等[9]分析了廣東省及鄰省共12個黃雞品種的mtDNA D-loop（片段長度520 bp）核苷酸多樣性，結果顯示廣東省黃雞的遺傳多樣性水平較高，品種的形成受到鄰省和北方家雞，以及東南亞紅原雞的影響?！颈狙芯壳腥朦c】然而，這些基于mtDNA D-loop的地方雞遺傳多樣性研究所用的序列長短不同，不同研究結果存在的差異使得不同研究比較難以進行平行比較分析，同時研究材料（品種）數量較少及覆蓋面較低，不利于摸清中國地方品種黃雞的遺傳多樣性水平及其群體歷史?！緮M解決的關鍵問題】本研究整合中國南方28個地方品種黃雞樣本及mtDNA D-loop數據，以統一的分析序列長度大尺度對中國南方地方品種黃雞的遺傳多樣性進行整體評估，探討品種的群體歷史，為中國家禽遺傳資源保護和創新利用提供科學依據。

1 材料與方法

1.1 試驗材料

新采集26個中國南方地方品種黃雞694份樣本，同時從GenBank（https://www.ncbi.nlm.nih.gov/ nuccore）數據庫下載17個中國南方地方品種黃雞589條mtDNA D-loop序列數據，合計28個地方品種共1 283份樣本（表1）。五華三黃雞之外的27個品種均入選了《中國畜禽遺傳資源志：家禽志》[1]，五華三黃雞入選了《中國禽類遺傳資源》[10]。采集地為保種場。根據系譜記錄采集無血緣關系的個體，肱靜脈取血100 μL，置于70%乙醇于-80℃保存。采用標準的酚氯仿方法提取基因組DNA，-20℃保存備用。

1.2 PCR擴增與序列測定

mtDNA D-loop的PCR擴增引物為：L16750，5￠-AGGACTACGGCTTGAAAAGC-3￠[11]；H522，5￠-ATG TGCCTGACCGAG-GAACC-AG-3￠[12]。反應體系為30 μL，含3 μL10×PCR Buffer（含Mg2+），2.4 μLdNTP mixture（2.5 mmol·L-1），正、反向引物（20 μmol·L-1）各0.3 μL，1 μL模板DNA，0.3 μLDNA聚合酶（寶生物，大連）（5 U/ μL）。PCR擴增條件：94℃預變性4min；35個循環（94℃變性30 s，63℃復性1 min，72℃延伸50 s）；最后72℃延伸10 min。PCR擴增產物用1.5%瓊脂糖凝膠電泳檢測到目標條帶后，送廣州艾基生物技術有限公司雙向測序。

1.3 序列分析

使用軟件ChromasPro 2.1.8（Technelysium，澳大利亞）檢查測序質量；以紅原雞mtDNA基因組全序列（NC_007235）作為參考序列，用MEGA 7.0[13]比對序列后，基于序列長度的分布情況和方便不同研究結果的平行比較分析，本文統一截取序列長度為519 bp（對應NC_007235核苷酸序列位置nps：1—519）作進一步分析。提取線粒體的SNP位點，利用Structure 2.3.4[14]分析品種間的群體遺傳結構。使用DnaSP 6.12.03[15]定義單倍型，提取變異位點信息（variable sites），計算單倍型多樣性（haplotype diversity）、核苷酸多樣性（nucleotide diversity）和錯配分布分析（mismatch distribution analysis）。利用Arlequin 3.5[16]以省份和品種為劃分依據，分別計算群體間分化指數（genetic differences among population，ST）、組間分化指數（genetic differences among groups defined a priori，CT）、組內群體間分化指數（genetic differences among population within groups，SC），進行群體內和群體間的分子方差分析（Analysis of Molecular Variance，AMOVA）以及中性檢測（Tajima’s和Fu's檢驗）。ST數值用于主坐標分析（principal coordinate analysis，PCoA），該計算由GenAlEx 6.503[17]執行。應用軟件MitoToolPy（http://www.mitotool.org/）獲得序列變異信息，然后劃入特定的單倍型類群[18-19]。中介網絡圖（median-joining network）由NETWORK 10.2.0.0[20]繪制。

表1 樣品信息

樣品/數據來源：1：本研究；2：Miao et al，2013；3：Gao et al，2017；4：黃勛和等，2016；5：黃勛和等，2018；6：Liu et al，2006

Sample or data references: 1: This study; 2: Miao et al, 2013; 3: Gao et al, 2017; 4: Huang et al, 2016; 5: Huang et al, 2018; 6: Liu et al, 2006

2 結果

2.1 遺傳多樣性

獲得了1 283條519 bp的mtDNA D-loop序列，檢測到101個變異位點，占總分析位點的19.46%，其中92個為多態位點，定義了169種單倍型。28個黃雞群體的總體單倍型多樣性為0.903±0.005，總體核苷酸多樣性為0.01269±n.d.。河南省和湖南省黃雞的單倍型多樣性和核苷酸多樣性均排在10個省自治區中的前二位，而海南省的單倍型多樣性最低，江西省的核苷酸多樣性最低（表2）。在地方品種中，淮南麻黃雞、固始雞和鄖陽大雞的單倍型多樣性位列28個品種的前三位，江漢雞、淮南麻黃雞和黃郎雞的核苷酸多樣性位列前三位，而洪山雞、廣西三黃雞和蕭山雞的單倍型多樣性和核苷酸多樣性均排在后三位。值得注意的是，蕭山雞的核苷酸多樣性遠低于其他黃雞品種。

表2 28個中國南方地方品種黃雞的線粒體DNA D-loop遺傳多樣性

*：在＜0.05水平顯著差異；**：在＜0.01水平顯著差異*: Significant at＜0.05;**: Significant at＜0.01

2.2 單倍型地理分布格局

169種單倍型劃分為6個單倍型類群A-E和G。單倍型類群B為優勢單倍型類群，占總樣本數的49.26%，其次為C（18.55%）、E（16.37%）和A（15.43%），D和G分別只檢測到4個（0.31%）和1個（0.08%）（表3）。廣東省黃雞的單倍型數量最多（61），河南省次之（41），浙江?。?3）和海南?。?2）較少。河南省黃雞的單倍型類群組成最為豐富，覆蓋6個單倍型類群，安徽省、江西省和廣西壯族自治區有5個單倍型類群，而其他省份只有4種單倍型類群（表3）。中介網絡圖顯示單倍型類群A、B、C和E為優勢單倍型類群，在10個省份的黃雞品種均有分布；并呈星狀分布模式，即每個單倍型類群由一個主導單倍型和一些衍生單倍型組成（圖1）。單倍型類群A和B在28個黃雞品種中均有分布。除寧都黃雞和霞煙雞外，單倍型類群C在其他26個雞品種均有分布；淮北麻雞、洪山雞、寧都黃雞、蕭山雞缺失單倍型類群E。單倍型類群D和G為稀有單倍型類群，D分布于淮南麻黃雞、固始雞、寧都黃雞和霞煙雞中，G只分布于固始雞中。

連接點數字表示核苷酸轉換的位置，圓的大小對應單倍型頻率，不同省份用不同顏色標注

表3 28個中國南方地方品種黃雞線粒體DNA D-loop的單倍型及其在品種中的分布

1)：單倍型數量及該品種或地區的獨享型單倍型（括號數值）Numbers in this bracket indicate the private haplotypes to a given breed or region

為進一步探討黃雞稀有單倍型類群D和G的來源，我們分析了已發表的8 076份亞洲（含本研究數據）和太平洋地區的地方雞線粒體D-loop數據，定義了1 345個單倍型類群D個體和288個單倍型類群G個體。中國和東亞的單倍型類群D的分布頻率最低，分別為1.55%和1.94%；太平洋和東南亞地區的最高，分別為74.38%和54.00%；其中東南亞的單倍型數最為豐富（附表1）。在4 953份中國地方雞樣品中，檢測到76個單倍型類群D，主要分布于斗雞、黑雞、烏雞和藏雞中，其中藏雞分布頻率最高，達17.73%（50/282）（附表2）。定義的亞洲單倍型類群G個體有288個，主要分布在云南（164），其次是四川（46）、南亞（11）和東南亞（11），而廣西（4）、貴州（2）和河南（1）只有少量分布。值得注意的是，藏雞的單倍型類群G分布頻率較高，達17.38%（49/282）。中介網絡圖顯示河南和南亞的單倍型類群G可能來源于西南地區（附圖1）。

2.3 遺傳分化

主坐標分析顯示品種間的遺傳距離主要與地理位置相關，但浙江、廣西壯族自治區和湖南的黃雞品種遺傳距離較為分散，如蕭山雞與仙居雞、霞煙雞與其他廣西黃雞、黃郎雞與桃源雞和東安雞，而河田雞與文昌雞的遺傳距離最為接近（圖2）。分子變異分析顯示遺傳變異主要來源于品種內，其次是品種間，均顯著差異；省份之間的變異很小，差異不顯著（表4）。在群體遺傳結構分析中，= 2時的Δ值遠大于= 3或其他值（附圖2）。Δ值與群體分化信號的強度呈正相關，推測= 2為最優聚類。然而，黃雞品種內部分個體雖含有混合組分，但品種間群體分化不明顯，未形成明顯的群體遺傳結構（圖3）。

2.4 群體歷史

中性檢測顯示只有河南、湖北和廣東黃雞的Fu’s數值達到顯著水平，但Tajima’s檢驗數值均不顯著；在品種水平上，只有蕭山雞的Tajima’s檢驗數值達到顯著水平，而其他品種數值均不顯著（表2）?；谑》荩▓D4）和品種（附圖3）的黃雞核苷酸錯配分布圖均為多峰，不符合種群擴張的單峰模式。優勢單倍型類群A、B和D中性檢測數值達到顯著水平（表2），并且核苷酸錯配分布圖為單峰，但單倍型類群C不顯著（圖4）。

A：基于省區；B：基于品種 A: Based on provinces; B: Based on breeds

表4 28個黃雞品種的群體遺傳結構分子變異分析

CT、SC和ST分別反映組間、群體間和群體內的遺傳變異

CT,SCandSTrepresent genetic differences among groups defined a priori, among populations within group, and within populations, respectively

每條豎線表示每個個體被分配到給定聚類中的一個，每個品種用黑色線條分開。GS：固始雞；ZY：正陽三黃雞；HB：淮北麻雞；HN：淮南麻黃雞；HS：洪山雞；JH：江漢雞；YY：鄖陽大雞；TY：桃源雞；HL：黃郎雞；DA：東安雞；BE：白耳黃雞；CR：崇仁麻雞；ND：寧都黃雞；XJ：仙居雞；XS：蕭山雞；HT：河田雞；HY：惠陽胡須雞；QY：清遠麻雞；WH：五華三黃雞；XH：杏花雞；ZS：中山沙欄雞；HX：懷鄉雞；GXP：廣西麻雞；XY：霞煙雞；YA：瑤雞；GXY：廣西三黃雞；WC：文昌雞

橫坐標為成雙配對差異，縱坐標為頻率；實線表示觀察值，虛線表示預測值

3 討論

本研究應用mtDNA D-loop分子標記分析了中國南方28個地方品種黃雞的遺傳多樣性和單倍型地理分布格局，探討了品種間的遺傳分化和群體歷史。總體而言，這些黃雞品種保持著較高的遺傳多樣性水平，但洪山雞、蕭山雞和廣西三黃雞核苷酸多樣性明顯較低。洪山雞規模較小，開發利用不充分；而廣西三黃雞分布較廣，開發應用程度較高；蕭山雞近期群體數量急劇減少，在Tajima’s檢驗結果也得以體現，這些因素可能與這3個品種遺傳多樣性較低有關[1]。因此，這3個品種尤其是蕭山雞，保種工作應引起足夠的重視，應采取提升保種群體數量、家系選育等措施，以逐漸提高其遺傳多樣性水平。

不同的樣本來源、樣本數量或分子標記長度差異會對研究結果造成一定的影響。如廣西地方雞，雖然在LIAO等[21]的研究中核苷酸多樣性均高于本研究，但當取相同D-loop長度（nps 189—509）分析時，除了廣西三黃雞之外，其他3個品種均高于LIAO等[21]的研究。在江西地方雞遺傳多樣性研究中，GAO等[22]采用D-loop全序列分析，其遺傳多樣性均低于本研究，這可能與D-loop的突變位點分布頻率的差異有關，因為D-loop的突變位點大部分集中核苷酸序列167—446的位置；而當取相同長度分析時，仍然低于本研究的結果，如白耳黃雞（0.00739）、崇仁麻雞（0.00770），這可能與GAO等[22]的樣本來源有關。類似的差異還有賈曉旭等[7]的研究等。這些差異不利于不同研究的平行比較分析。因此，本研究統一了D-loop分子標記長度，從整體上對中國南方主要的地方品種黃雞進行了遺傳多樣性水平評估，分析結果對制定針對性的保種策略，提高保種利用率具有積極的參考作用。

單倍型類群地理格局分析有助于追溯動物的群體歷史[3,23-24]。本研究來源于28個地方品種黃雞的1 283份樣本共檢測到6個單倍型類，其中A、B、C和E是黃雞常見單倍型類群，在黃雞品種中較均勻分布，而D（0.31%）和G（0.07%）則是稀有單倍型類群，這與有關研究結果相一致[7,9,18,21-23]。單倍型類群D在中國的分布頻率較低（1.53%），主要分布于斗雞和藏雞群體[4,18,23]，但在東南亞、太平洋等海島地區，以及馬達加斯加等部分地區分布廣泛且單倍型豐富，甚至可達89.13%[24-26]，并且中性檢測顯示東南亞地區家雞群體曾經歷了顯著的歷史擴張。雖然太平洋等海島地區的家雞單倍型類群D分布頻率高于其他地區，但其主要來源于東南亞[18,26-27]。中國沿海地區如廣東、福建等地自古以來與東南亞交流頻繁?；诖耍茢帱S雞單倍型類群D可能是通過商貿等途徑從東南亞地區攜帶而來。有趣的是，在河南固始雞中檢測到1例單倍型類群G的個體，另外兩個研究中也有報道[18,28]。該單倍型集中分布在四川、云南地方雞和藏雞群體中，并且單倍型類群G在西南地區尤其是云南普遍分布，而在中國其他地區的地方雞品種中均未檢測到。因此，固始雞的單倍型類群G可能是從西南地方雞種中獲得，但具體途徑和時間還有待查證。

分子變異分析顯示變異主要來自于品種內，品種間的遺傳變異也達到顯著水平，說明品種的遺傳特性明顯。但黃雞的遺傳結構不清晰，地區之間分化也不顯著。一方面，這是持續多年保種努力的結果，說明黃雞種群處于穩定較好的保種狀態。國家于2005年12月通過了《中華人民共和國畜牧法》，并于2015年4月進行了修正，將49個家禽品種列為國家級保護品種，97個列為省級保護品種。廣東省早在20世紀50年代就開展了地方雞的種質資源調查工作，以及后續的一系列復查、再核查、反復論證和補充調查等大量工作，挖掘出一批優質雞品種并建立了相應的保種場，制定了科學的保護策略，取得了較好的效果[9,29-30]。另一方面，由于地方雞沒有嚴格的特定分布/飼養范圍，不同品種間可能存在混養/雜交現象，該假設得到了來自28個品種的核苷酸錯配分布（典型雙峰分布）和遺傳結構分析的支持。并且不同品種建立保種場的時間以及重視程度也存在差異，因此地方品種黃雞的遺傳分化不明顯。商品雞在新品種（配套系）的應用對黃雞品種的遺傳特性也具有一定的影響。比如，嶺南黃雞有隱性白洛克的血統，而嶺南黃雞在華南地區的大量飼養對地方雞有一定的滲透作用[31]；雞60K SNP芯片檢測研究發現中國地方雞存在廣泛的商品雞基因滲透現象[32]。此外，由于大部分地方品種黃雞的分布范圍和飼養規模有限，以及悠久的品種形成歷史和復雜的基因交流，因此未檢測到以品種/省份為單元的明顯擴張現象。然而，單倍型類群A、B和E錯配分布分析和中性檢測表明地方品種黃雞整體上經歷了顯著的群體擴張歷史，形成了目前豐富的地方品種資源分布格局，但地方品種黃雞群體歷史的具體細節還需進一步研究。

由于不同分子標記的分辨率不同，對結果解析度也存在差異。D-loop突變速率較快，基因覆蓋度低，遺傳分化解析度不如微衛星分子標記或全基因組。本研究的線粒體DNA群體遺傳分化不明顯，但相同的樣本在微衛星或全基因組研究中卻可以明確地劃分群體遺傳結構。HUANG等[31]運用微衛星分子標記將12個華南家雞分為3個明顯的群體遺傳結構，主要與地理位置有關。HUANG等[33]全基因組重測序分析了10個傳統黃雞的遺傳結構，主成分分析（Principal component analysis）、混合分析（ADMIXTURE）和鄰接系統發生樹分析（Neighbor-joining phylogenetic tree）均顯示這些雞品種可明顯分為北部、中部和南部，與微衛星標記分析的趨勢一致[31,34]。在主坐標分析中，文昌雞與河田雞、五華三黃雞等客家地區的品種關系較近，與全基因組分析結果一致[33]，也與歷史記錄相符[1]。據記載，文昌雞早在清代由廣東、福建地區的客家移民帶入雜交選育而成。因此，遺傳分析結果與歷史記錄相結合，是對了解品種形成歷史的有效方法之一。

4 結論

中國南方地方品種黃雞線粒體DNA D-loop數據整合分析顯示，大部分黃雞群體保留著較高的遺傳多樣性和群體間遺傳變異水平，處于穩定的保種狀態，但洪山雞、蕭山雞和廣西三黃雞應加強保種工作。中國南方地方品種黃雞遺傳組分混雜，缺乏明顯的遺傳結構，反映了其復雜的雜交起源歷史。東南亞和中國西南地方雞對中國南方地方品種黃雞有一定的遺傳貢獻。

[1] 國家畜禽遺傳資源委員會. 中國畜禽遺傳資源志-家禽志. 北京: 中國農業出版社, 2011.

China National Commission of Animal Genetic Resources. Animal Genetic Resources in China: Poultry. Beijing: Chinese Agriculture Press, 2011. (in Chinese)

[2] 孟祥兵, 徐廷生, 杜炳旺, 滕小華. 雞藝——中國古代養雞智慧附書法藝術[M]. 北京: 中國農業出版社, 2017.

Meng X B, Xu T S, Du B W, Teng X H. Chicken art - ancient Chinese wisdom on chicken-raising with an appendix of the Chinese art of calligraphy. Beijing: Chinese Agriculture Press, 2017. (in Chinese)

[3] Lan D, Hu YD, Zhu Q, Liu Y P. Mitochondrial DNA study in domestic chicken. Mitochondrial DNA Part A: DNA Mapping, Sequencing, and Analysis, 2017, 28: 25-29.

[4] Zhang L, Zhang P, Li Q Q, Gaur U, Liu Y P, Zhu Q, Zhao X L, Wang Y, Yin H D, Hu Y D, Liu A P, Li D Y. Genetic evidence from mitochondrial DNA corroborates the origin of Tibetan chickens. PLoS ONE, 2017, 12: e0172945.

[5] WU Y P, HUO J H, XIE J F, LIU L X, WEI Q P, XIE M G, KANG Z F, JI H Y, MA Y H. Phylogeography and origin of Chinese domestic chicken. Mitochondrial DNA, 2014, 25(2): 126-130. doi:10.3109/ 19401736.2013.786704.

[6] 黃勛和, 翁茁先, 李威娜, 陳潔波, 鐘福生, 唐壽桂, 趙宇旗. 基于線粒體DNAD-loop序列的黃郎雞遺傳多樣性與品種起源研究. 湖南農業大學學報(自然科學版), 2016(1): 75-80.

HUANG X H, WENG Z X, LI W N, CHEN J B, ZHONG F S, TANG S G, ZHAO Y Q. Genetic diversity and breed origin of Huanglang chicken inferred from mitochondrial DNA D-loop sequence. Journal of Hunan Agricultural University (Natural Sciences), 2016(1): 75-80. (in Chinese)

[7] 賈曉旭, 唐修君, 樊艷鳳, 陸俊賢, 黃勝海, 葛慶聯, 高玉時, 韓威. 華東地區地方雞品種mtDNA控制區遺傳多樣性. 生物多樣性, 2017, 25(5): 540-548. doi:10.17520/biods.2017012.

JIA X X, TANG X J, FAN Y F, LU J X, HUANG S H, GE Q L, GAO Y S, HAN W. Genetic diversity of local chicken breeds in East China based on mitochondrial DNA D-loop region. Biodiversity Science, 2017, 25(5): 540-548. doi:10.17520/biods.2017012. (in Chinese)

[8] JIA X X, LU J X, TANG X J, FAN Y F, HUANG S H, GE Q L, GAO Y S. Genetic diversity of Jiangsu native chicken breeds assessed with the mitochondrial DNA D-loop region. British Poultry Science, 2018, 59(1): 34-39. doi:10.1080/00071668.2017.1395391.

[9] 黃勛和, 余哲琪, 翁茁先, 何丹林, 易振華, 李威娜, 陳潔波, 張細權, 杜炳旺, 鐘福生. 廣東省地方雞線粒體遺傳多樣性與母系起源. 生物多樣性, 2018, 26(3): 238-247. doi:10.17520/biods.2017259.

HUANG X H, YU Z Q, WENG Z X, HE D L, YI Z H, LI W N, CHEN J B, ZHANG X Q, DU B W, ZHONG F S. Mitochondrial genetic diversity and maternal origin of Guangdong indigenous chickens. Biodiversity Science, 2018, 26(3): 238-247. doi:10.17520/ biods.2017259. (in Chinese)

[10] 陳國宏. 中國禽類遺傳資源. 上海: 上?？茖W技術出版社, 2004.

CHEN G H. Poul Try Genetic Resources in China. Shanghai: Shanghai Scientific & Technical Publishers, 2004. (in Chinese)

[11] FumihitoA, Miyake T, Sumi S, Takada M, Ohno S, Kondo N. One subspecies of the red jungle fowl () suffices as the matriarchic ancestor of all domestic breeds. Proceedings of the National Academy of Sciences, USA, 1994, 91: 12505-12509.

[12] 傅衍, 牛冬, 羅靜, 阮暉, 何國慶, 張亞平. 中國家雞的起源探討. 遺傳學報, 2001, 28(5): 411-417.

FU Y, NIU D, LUO J, RUAN H, HE G Q, ZHANG Y P. Studies of the origin of Chinese domestic fowls. Acta Genetica Sinica, 2001, 28(5): 411-417. (in Chinese)

[13] Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Molecular Biology and Evolution, 2016, 33: 1870-1874.

[14] Pritchard J K, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics, 2000, 155(2): 945-959.

[15] RozasJ, Ferrer-Mata A, Sánchez-DelBarrio J C, Guirao-Rico S, Librado P, Ramos-Onsins S E, Sánchez- Gracia A. DnaSP 6: DNA sequence polymorphism analysis of large datasets. Molecular Biology Evolution, 2017, 34: 3299-3302.

[16] Excoffier L, Lischer H E. Arlequin suite ver. 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources, 2010, 10: 564-567.

[17] Peakall R, Smouse P E. GENALEX6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes, 2006, 6: 288-295.

[18] Miao Y W, Peng M S, Wu G S, Ouyang Y N, Yang Z Y, Yu N, Liang J P, Pianchou G, Beja-Pereira A, Mitra B, Palani-chamy M G, Baig M, Chaudhuri T K, Shen Y Y, Kong Q P, Murphy R W, Yao Y G, Zhang Y P. Chicken domestic-ation: An updated perspective based on mitochondrial gen-omes. Heredity, 2013, 110: 277-282.

[19] Peng M S, Fan L, Shi N N, Ning T, Yao Y G, Murphy R W, Wang W Z, Zhang Y P. DomeTree: a canonical toolkit for mitochondrial DNA analyses in domesticated animals. Molecular Ecology Resources, 2015, 15: 1238-1242.

[20] Bandelt H J, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Molecular Biology and Evolution, 1999, 16: 37-48.

[21] Liao Y Y, Mo G D, Sun J L, Wei F Y, Liao D J. Genetic diversity of Guangxi chicken breeds assessed with microsa-tellites and the mitochondrial DNA D-loop region. Molecu-lar Biology Reports, 2016, 43: 415-425.

[22] Gao Y S, Jia X X, Tang X J, Fan Y F, Lu J X, Huang S H, Tang M J. The genetic diversity of chicken breeds from Jiangxi, assessed with BCDO2 and the complete mitochondrial DNA D-loop region. PLoS ONE, 2017, 12: e0173192.

[23] Huang XH, Wu YJ, Miao YW, Peng MS, Chen X, He DL, Suwannapoom C, Du BW, Li XY, Weng ZX, Jin SH, Song JJ, Wang MS, Chen JB, Li WN, Otecko NO, Geng ZY, Qu XY, Wu YP, Yang XR, Jin JQ, Han JL, Zhong FS, Zhang XQ, Zhang YP. Was chicken domesticated in north-ern China? New evidence from mitochondrial genomes. Science Bulletin, 2018, 63: 743-746.

[24] Herrera M, Thomson V A, Wadley J J, Piper P J, Sulandari S, Dharmayanthi A B, Kraitsek S, Gongora J, Austin J J. East African origins for Madagascan chickens as indicated by mitochondrial DNA. Royal Society Open Science, 2017(4): 160787.

[25] Dancause K N, Vilar M G, Steffy R, Lum J K. Characterizing genetic diversity of contemporary pacific chickens using mitochondrial DNA analyses. PloS ONE, 2011, 6(2): e16843.

[26] Langford S M, Kraitsek S, Baskerville B, Ho S Y W, Gongora J. Australian and Pacific contributions to the genetic diversity of Norfolk Island feral chickens. BMC Genetics, 2013, 14: 91.

[27] Thomson V A, Lebrasseur O, Austin J J, Hunt T L, Burney D A, Denham T, Rawlence N J, Wood J R, Gongora J, Flink G L, Linderholm A, Dobney K, Larson G, Cooper A. Using ancient DNA to study the origins and dispersal of ancestral Polynesian chickens across the Pacific. Proceedings of the National Academy of Sciences, USA, 111(13): 4826-4831.

[28] Liu Y P, Wu G S, Yao Y G, Miao Y W, Luikart G, Baig M, Beja-Pereira A, Ding Z L, Palanichamy M G, Zhang Y P. Multiple maternal origins of chickens: Out of the Asian jungles. Molecular Phylogenetics and Evolution, 2006, 38: 12-19.

[29] 瞿浩, 舒鼎銘, 楊純芬. 廣東優質雞種質資源的現狀及發展建議. 廣東農業科學, 2004, 31(6): 4-7.

Qu H, Shu D M, Yang C F. The actuality and develop-ment proposal in Guangdong quality chicken germplasm resources. Guangdong Agricultural Sciences, 2004, 31(6): 4-7. (in Chinese)

[30] 張細權, 梁少東, 張德祥. 廣東優質雞新品種選育回顧. 中國家禽, 2006, 28(1): 9-11.

Zhang X Q, Liang S D, Zhang D X. Review on new breed selection of quality chicken in Guangdong Province. China Poultry, 2006, 28(1): 9-11. (in Chinese)

[31] Huang X H, Zhang J F, He D L, Zhang X Q, Zhong F S, Li W N, Zheng Q M, Chen J B, Du B W. Genetic diversity and population structure of indigenous chicken breeds in South China. Frontiers of Agricultural Science and Engineering, 2016, 3: 97-101.

[32] Zhang C Y, Lin D, Wang Y Z, PENG D Z, LI H F, FEI J, CHEN K W, YANG N, HU X X, ZHAO Y Q, LI N. Widespread introgression in Chinese indigenous chicken breeds from commercial broiler. Evolutionary applications, 2019, 12(3): 610-621.

[33] Huang X H, Otecko N O, Peng M S, Weng Z X, Li W N, Chen J B, Zhong M, Zhong F S, Jin S H, Geng Z Y, Luo W, He D L, Ma C, Han J L, Ommeh S C, Zhang Y P, Zhang X Q, Du B W. Genome-wide genetic structure and selection signatures for color in 10 traditional Chinese yellow-feathered chicken breeds. BMC Genomics, 2020, 21: 316.

[34] Qu L J, Li X Y, Xu G F, Chen K W, Yang H J, Zhang L C, Wu G Q, Hou Z C, Xu G Y, Yang N. Evaluation of genetic diversity in Chinese indigenous chicken breeds using mic-rosatellite markers. Science in China Series C: Life Science, 2006, 49: 332-341.

Genetic Diversity of Indigenous Yellow-Feathered Chickens in Southern China Inferred from Mitochondrial DNA D-Loop Region

1Jiaying University/Guangdong Provincial Key Laboratory of Conservation and Precision Utilization of Characteristic Agricultural Resources in Mountainous Areas/Guangdong Science and Technology Innovation Centre of Wuhua Yellow Chicken, Meizhou 514015, Gugangdong;2College of Animal Sciences, South China Agricultural University, Guangzhou 510642

【Objective】Indigenous yellow-feathered chicken is an important poultry genetic resource in China. The genetic diversity of these chickens was evaluated to facilitate scientific protection policies and to improve their use.【Method】 A total of 1283 mitochondrial DNA fragments (mtDNA D-loop, 519 bp) were investigated including 694and 589 previously published units from 28 indigenous yellow-feathered chicken breeds from southern China. These data were subsequently used to calculate the genetic diversity, to construct median-joining networks of haplotypes, and to investigate geographical distribution patterns of haplotypes by using principal coordinate analysis, analysis of molecular variance, and neutral test. The origin of rare haplogroups was inferred from the geographical distribution pattern of mitochondrial DNA haplogroups of Asian and Pacific indigenous chickens.【Result】One hundred and one mutation sites were detected in 1 283 samples, of which 92 were polymorphic. One hundred and sixty-nine haplotypes, belonging to haplogroups A-E and G, were defined. The predominant haplogroups were A (15.43%), B (49.26%), C (18.55%), and E (16.37%), while D (0.31%) and G (0.08%) were the rare haplogroups within the total samples. Chickens from Henan Province covered all six haplogroups; Chickens from Guangdong Province had the largest number of haplotypes (64), whereas those from Zhejiang (19) and Hainan (12) Provinces had the lowest. Haplogroups A and B occurred in all 28 breeds. Haplogroup D existed only in Huainan yellow, Gushi, Ningdu yellow and Xiayan chickens, whereas haplogroup G only occurred in Gushi chickens. The distribution frequency of haplogroup D and the number of haplotypes indicated that haplogroup D of indigenous yellow-feathered chickens in southern China could originated from Southeast Asia. The evidence from geographical distribution and median-joining networks of haplogroup G indicated that haplogroup G of chickens from Henan Province and South Asia could originated from Southwest China. The indigenous yellow-feathered chickens in southern China had a total haplotype diversity of 0.903±0.005 and nucleotide diversity of 0.01269±n.d. Indigenous yellow-feathered chickens from Henan and Hunan Provinces had the highest level of haplotype diversity of 0.916±0.011 and 0.01358±0.00039, respectively, and the nucleotide diversity of 0.913±0.012 and 0.01345±0.00042, respectively. The lowest level of haplotype diversity was found in Hainan Province (0.736±0.076), and the lowest level of nucleotide diversity was found in Jiangxi Province (0.00981±0.00072). Huainan yellow, Gushi and Yunyang da chickens retained the highest level of haplotype diversity; Jianghan, Huainan yellow and Huanglang chickens had the highest level of nucleotide diversity, whereas the haplotype diversity and nucleotide diversity in Hongshan, Guangxi yellow and Xiaoshan chickens were the lowest. Wenchang chickens were genetically close to Hentian chickens, which was consistent with their breeding history. The population genetic structure of yellow-feathered chickens was less clear, and the molecular variance component of within-population was significantly higher than others. The neutral test indicated that the yellow-feathered chickens in southern Chinadid not sustain obvious population expansion at the breed level, with the exception of haplogroups A, B and E.【Conclusion】The results suggested that yellow-feathered chickens had a good conservation status with high level of genetic diversity, although Hongshan, Xiaoshan and Guangxi yellow chickens required further protection. Hybridization between chickens was common in yellow-feathered chickens, which had experienced population expansion. Indigenous chickens from Southeast Asia and Southwest China had substantial genetic contributions to indigenous yellow-feathered chickens in southern China.

yellow-feathered chicken; mitochondrial DNA D-loop; genetic variation; genetic divergence; demographic history

10.3864/j.issn.0578-1752.2022.22.016

2021-09-13；

2021-12-16

廣東省自然科學基金（2014A030307018）、廣東省公益研究與能力建設項目（2016A030303068）、廣東省高等教育“沖補強”提升計劃重點建設學科（農業資源與環境）建設項目（粵教科函[2018]181號）、梅州市應用型科技專項資金項目（2019B0201003、2019B0201006）

黃勛和，E-mail：hxh826@126.com。通信作者杜炳旺，E-mail：dudu903@163.com

（責任編輯 林鑒非）