加入放歸個(gè)體后的栗子坪保護(hù)區(qū)大熊貓小種群遺傳結(jié)構(gòu)研究

陳子君， 黃蜂，2， 何可, 馮菲菲,古曉東,楊旭煜,戴強(qiáng), 齊敦武,楊志松

(1.西華師范大學(xué) 西南野生動(dòng)植物資源保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室， 四川　南充　637002；2.四川栗子坪國(guó)家級(jí)自然保護(hù)區(qū)管理處， 四川　石棉　625400；3.四川省野生動(dòng)物資源調(diào)查保護(hù)管理站， 四川　成都　610081；4.中國(guó)科學(xué)院成都生物研究所， 四川　成都　610041；5.成都大熊貓繁育研究基地, 四川　成都　610086)

?

加入放歸個(gè)體后的栗子坪保護(hù)區(qū)大熊貓小種群遺傳結(jié)構(gòu)研究

陳子君1， 黃蜂1，2， 何可1, 馮菲菲1,古曉東3,楊旭煜3,戴強(qiáng)4, 齊敦武5,楊志松1

(1.西華師范大學(xué) 西南野生動(dòng)植物資源保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室， 四川南充637002；2.四川栗子坪國(guó)家級(jí)自然保護(hù)區(qū)管理處， 四川石棉625400；3.四川省野生動(dòng)物資源調(diào)查保護(hù)管理站， 四川成都610081；4.中國(guó)科學(xué)院成都生物研究所， 四川成都610041；5.成都大熊貓繁育研究基地, 四川成都610086)

小相嶺山系是現(xiàn)有大熊貓種群數(shù)量最少的山系之一。國(guó)家林業(yè)局?jǐn)M通過(guò)在該山系的栗子坪國(guó)家級(jí)自然保護(hù)區(qū)放歸大熊貓“瀘欣”、“淘淘”、“張想”來(lái)復(fù)壯大熊貓小種群。為評(píng)估加入3只放歸個(gè)體后的小種群遺傳結(jié)構(gòu)，2012～2013年在研究區(qū)域采集大熊貓糞便樣品92份，鑒定出14只個(gè)體(含放歸個(gè)體)。放歸個(gè)體加入后，測(cè)得放歸個(gè)體與本地大熊貓小種群親緣關(guān)系較遠(yuǎn)。種群年齡組成為青年組8只、成年組4只、未知年齡個(gè)體2只。平均觀測(cè)雜合度(Ho)為0.557，平均期望雜合度(He)為0.581，遺傳多樣性處于中等偏上水平。5個(gè)位點(diǎn)上共計(jì)增加6個(gè)等位基因，平均等位基因數(shù)(MNA)達(dá)到3.667。遺傳基因簇由1簇增加為4簇，出現(xiàn)新的基因頻率。以上結(jié)果顯示放歸個(gè)體加入后該區(qū)大熊貓小種群遺傳多樣性有所提高，放歸成效初顯。

大熊貓；小種群；遺傳多樣性；放歸；栗子坪國(guó)家級(jí)自然保護(hù)區(qū)

放歸是促進(jìn)小種群復(fù)壯，恢復(fù)破碎化棲息地的有效方法[1-2]。國(guó)內(nèi)外已有不少成功案例，國(guó)內(nèi)的如普氏野馬(Equus ferus przewalskii)、麋鹿(Elaphurus davidianus )、揚(yáng)子鱷(Alligator sinensis)等。國(guó)外的如加拿大盤(pán)羊(Ovis canadensis)、美洲野牛(Bison bison)、美洲黑熊(Ursus americanus)等。然而，也有少數(shù)放歸活動(dòng)因近交[3]、遺傳多樣性丟失[4]或者人類(lèi)干擾、天敵捕食[5]而成為失敗的嘗試。因此，及時(shí)分析放歸成效，總結(jié)經(jīng)驗(yàn)、教訓(xùn)，對(duì)放歸工作的進(jìn)一步開(kāi)展有著極為重要的意義。

小相嶺山系是現(xiàn)有大熊貓種群數(shù)量最少的山系之一，種群生存力分析表明，如果與其他種群沒(méi)有遺傳交流或者無(wú)外來(lái)個(gè)體補(bǔ)充，該小種群極有可能在80年內(nèi)會(huì)趨于滅絕[6]。為提高該區(qū)大熊貓小種群生存能力，促進(jìn)大熊貓小種群的復(fù)壯，(截止到本研究開(kāi)展的2013年)國(guó)家林業(yè)局在小相嶺山系的栗子坪國(guó)家級(jí)保護(hù)區(qū)先后放歸了3只大熊貓(Ailuropoda melanoleuca)：“瀘欣”已放歸超過(guò)6年，“淘淘”已放歸超過(guò)3年，“張想”也已放歸超過(guò)2年。然而，加入放歸個(gè)體后，栗子坪自然保護(hù)區(qū)大熊貓小種群遺傳結(jié)構(gòu)是否有所改善，目前仍不清楚。本研究利用野外調(diào)查收集的大熊貓糞便，通過(guò)糞便咬節(jié)測(cè)量、糞便DNA性別鑒定以及微衛(wèi)星分析，確定了放歸地大熊貓種群的年齡結(jié)構(gòu)、性比和種群遺傳結(jié)構(gòu)，并從種群遺傳上對(duì)放歸成效進(jìn)行了評(píng)估。這一評(píng)估結(jié)果有助于確定將來(lái)大熊貓放歸和小種群復(fù)壯項(xiàng)目的方向和重點(diǎn)，對(duì)其他瀕危珍稀野生動(dòng)物的小種群復(fù)壯工作也有很好的參考價(jià)值。

1　研究對(duì)象及區(qū)域

圖1　研究區(qū)域區(qū)位示意圖

栗子坪自然保護(hù)區(qū)位于四川石棉縣境內(nèi)，地理位置為E 102°10′33″～102°29′07″，N 28°51′02″～29°08′42″；南北長(zhǎng)23 km，東西寬17.8 km，總面積47 940 hm2。區(qū)內(nèi)最高海拔為4 551 m，最低海拔1 330 m，最大相對(duì)高差達(dá)3 221 m，平均相對(duì)高差大于2 000 m。保護(hù)區(qū)是小相嶺山系大熊貓種群的核心分布區(qū)和關(guān)鍵棲息地，也是小相嶺山系大熊貓種群遺傳基因交流的關(guān)鍵走廊帶和重要集中分布區(qū)。

本項(xiàng)目研究區(qū)域位于栗子坪自然保護(hù)區(qū)內(nèi)公益海保護(hù)站周邊，面積約為23 251.45 hm2。“瀘欣”、“淘淘”、“張想”先后放歸于該區(qū)(詳見(jiàn)圖1和表1)，遺傳分析所用糞便也采集自這一區(qū)域。這一區(qū)域的大熊貓棲息地與周邊大熊貓棲息地存在較為嚴(yán)重的隔離[7-8]，第四次全國(guó)大熊貓調(diào)查也認(rèn)為該區(qū)域大熊貓種群是一個(gè)相對(duì)獨(dú)立的局域種群[9]。

表1　放歸大熊貓個(gè)體信息

2　研究方法

2．1野外采樣

野外樣品于2012年4月至2013年12月采集自保護(hù)區(qū)內(nèi)公益海保護(hù)站周邊(見(jiàn)圖1)。 研究區(qū)域內(nèi)調(diào)查樣線的抽樣強(qiáng)度為1 條/200 hm2,采用非損傷性遺傳取樣方法在野外收集糞便樣品[10]。 在采樣同時(shí)，對(duì)所有的樣品采集點(diǎn)均進(jìn)行GPS定位并記錄。

2．2大熊貓年齡組劃分

測(cè)量采集到的大熊貓糞便的咬節(jié)長(zhǎng)度，采用“離差平方和法”對(duì)咬節(jié)長(zhǎng)度進(jìn)行聚類(lèi)，參照胡錦矗的年齡組劃分標(biāo)準(zhǔn)[11]確定個(gè)體年齡段：幼年組(小于1.5 歲)，咬節(jié)值小于30 mm，糞便直徑30 mm以下；青年組(大于等于1.5歲，并且小于5.5 歲),咬節(jié)值在30～35 mm(包含30 mm) 之間, 糞便直徑30～55 mm(包含30 mm)；成年組(大于等于5.5歲，并且小于20 歲)，咬節(jié)長(zhǎng)度在35～40 mm(包含35 mm)之間，糞便直徑55～65 mm (包含55 mm)；老年組(大于等于20 歲)，咬節(jié)長(zhǎng)度在40 mm(包含40 mm)以上，糞便直徑在65 mm(包含65 mm)以上。

2．3糞便微衛(wèi)星DNA提取與分析

使用QIAGEN公司的QIAamp DNA Stool kit試劑盒進(jìn)行大熊貓糞便DNA的提取，性別鑒定PCR、微衛(wèi)星PCR以及結(jié)果判定參考胡義波[12]、詹祥江[10]的方法，我們挑選了12對(duì)微衛(wèi)星引物(Ame -μ10, Ame -μ11, Ame -μ22, Ame -μ13, Ame -μ15, Ame -μ24, Ame -μ26, Ame -μ27以及AY161179, AY161195, AY161213, AY161217)來(lái)對(duì)大熊貓糞便DNA進(jìn)行微衛(wèi)星擴(kuò)增[12]。

因大熊貓糞便中DNA含量較低，為避免因引物選取、PCR擴(kuò)增、數(shù)據(jù)讀取等引起的誤差，我們利用MICRO-CHECKER軟件[13]檢測(cè)基因分型數(shù)據(jù)中是否存在無(wú)效等位基因(Null alleles)、大片段等位基因丟失(Large allele dropout)及鬼影峰誤判(Stutter scoring error)等分型錯(cuò)誤。利用GIMLET 1．3．3軟件[14]計(jì)算了12個(gè)座位的聯(lián)合P(ID)值，以確定選擇的微衛(wèi)星引物數(shù)量是否能將影子效應(yīng)(Shadow effect)[15]減少到一個(gè)不顯著的水平。

2．4個(gè)體鑒定與種群數(shù)量估計(jì)

利用MStools插件[16]，尋找基因分型數(shù)據(jù)中相匹配的基因型。判斷為同一個(gè)個(gè)體的標(biāo)準(zhǔn)是：所有座位上的基因型都相同，或者只有一個(gè)座位上的一個(gè)等位基因存在差別[17]。由于我們采樣分析是針對(duì)放歸成效的研究，糞便采集也集中在放歸大熊貓主要活動(dòng)的區(qū)域(依據(jù)放歸大熊貓佩戴的GPS項(xiàng)圈確定)，可以認(rèn)為出現(xiàn)在此區(qū)域外或是很少出現(xiàn)在該區(qū)域的個(gè)體被捕捉概率很低，因此我們選擇CAPWIRE軟件里TIRM(Two Innate Rates Model)模型[10]進(jìn)行種群數(shù)量估計(jì)。

2．5種群遺傳多樣性及親緣關(guān)系分析

利用ARLEQUIN 3．5[18]和GENEPOP 3．4[19]計(jì)算期望雜合度(He, expected heterozygosity)、表觀雜合度(Ho, observed heterozygosity)、等位基因數(shù)(A)，測(cè)算總體和各個(gè)位點(diǎn)是否符合哈溫平衡(Hardy-Weinberg)。利用FSTAT 2．9．3．2[20]來(lái)檢驗(yàn)位點(diǎn)是否存在連鎖不平衡(Linkage disequilibrium)。利用STRUCTURE方法[21]檢測(cè)研究區(qū)域所有個(gè)體中存在的遺傳基因簇。利用Kingroup v2．08[22]來(lái)測(cè)算個(gè)體間的親緣關(guān)系遠(yuǎn)近，其中個(gè)體間的親緣關(guān)系值(pairwise relatedness, r)參考Wang[23]的研究。

3　結(jié)果

3．1種群數(shù)量、性比及年齡結(jié)構(gòu)

本研究從92 份樣品(含糞便樣品及放歸個(gè)體血樣)中得到52個(gè)可靠基因型。TIRM模型分析表明，研究區(qū)域大熊貓種群數(shù)量(含3 只放歸個(gè)體)為25 只(95%置信范圍：14～28)。種群數(shù)量最低估計(jì)值(lowCI)為14 只(與實(shí)驗(yàn)室分析鑒定結(jié)果一致)，因此研究區(qū)域中大熊貓種群最保守的數(shù)量應(yīng)該是14 只。14只個(gè)體分布在公益海保護(hù)站區(qū)域，海拔集中在2 291～3 510 m(見(jiàn)圖2)，分別是8 只雌性、6只雄性，雌雄比例為4∶3。其中“瀘欣”、“淘淘”和“張想”為放歸個(gè)體，年齡有準(zhǔn)確記錄，無(wú)需通過(guò)糞便咬節(jié)確定年齡，另有兩只個(gè)體未能獲得有效咬節(jié)數(shù)據(jù)，無(wú)法判斷年齡，其余個(gè)體的年齡段均根據(jù)咬節(jié)法測(cè)得。結(jié)果表明，本地種群的年齡組成為：青年組6 只，成年組3 只，未知年齡個(gè)體2只(見(jiàn)表2)。從年齡結(jié)構(gòu)上來(lái)看，青年組個(gè)體數(shù)明顯多于成年組個(gè)體數(shù)，年齡組成為增長(zhǎng)型。 加入放歸個(gè)體后的大熊貓種群年齡組成為：青年組8 只、成年組4只，未知年齡個(gè)體2只(見(jiàn)表2)，與放歸前比較，年齡組結(jié)構(gòu)的增長(zhǎng)型特征更為明顯。

表2　本研究大熊貓種群年齡組組成

3．2本地小種群遺傳結(jié)構(gòu)

MICRO-CHECKER檢驗(yàn)表明，本研究基因分型數(shù)據(jù)結(jié)果未受無(wú)效等位基因、大片段等位基因丟失及鬼影帶誤判等分型錯(cuò)誤的顯著干擾。GIMLET分析顯示，12 個(gè)微衛(wèi)星座位的聯(lián)合區(qū)分率很高，兩只個(gè)體具有相同基因型的無(wú)偏概率P(ID)為3．286×10-9，即使為雙胞胎，判斷錯(cuò)誤的概率P(sib) 也僅為0．04778%，小于0．01[15]。哈溫平衡檢驗(yàn)結(jié)果顯示研究區(qū)域本地大熊貓小種群未偏離哈溫平衡(P>0．05)。連鎖不平衡分析結(jié)果顯示本研究采用的12 個(gè)位點(diǎn)，對(duì)當(dāng)?shù)卮笮茇埿》N群而言不存在連鎖遺傳。12個(gè)微衛(wèi)星位點(diǎn)分析結(jié)果顯示研究區(qū)域本地大熊貓小種群的平均觀測(cè)雜合度(Ho)為0．635，平均期望雜合度(He)為0．608。平均等位基因數(shù)(MNA)為3．455(見(jiàn)表3)。STRUCTURE 分析結(jié)果顯示K=1時(shí)(圖3A)，LnP(D)具有最大的平均值，圖3B表明研究區(qū)域本地大熊貓小種群整體上是一個(gè)遺傳基因簇。

表3　研究區(qū)域本地大熊貓小種群(代號(hào)A)與加入放歸個(gè)體后的種群(代號(hào)B)遺傳多樣性信息

圖3　LnP(D)值及遺傳基因簇圖圖A:本地小種群LnP(D)值；圖C:放歸個(gè)體加入后的種群LnP(D)值；圖B圖D分別為K=1、K=4的遺傳基因簇圖

3．3加入放歸個(gè)體后小種群遺傳結(jié)構(gòu)

3．3．1放歸個(gè)體和本地大熊貓親緣關(guān)系

親緣關(guān)系分析結(jié)果顯示(見(jiàn)表4)，“瀘欣”、“淘淘”、“張想”和整個(gè)本地小種群的親緣關(guān)系都較遠(yuǎn)，與其它13 只個(gè)體的親緣系數(shù)(r)平均值平均為-0．07、-0．21、-0．23。

雖然“瀘欣”和lzp6、lzp7、lzp10、lzp11、“張想”的親緣系數(shù)(r)>0，但是這五只個(gè)體均為雌性，不存在因親緣關(guān)系近而帶來(lái)的近交風(fēng)險(xiǎn)。“淘淘”和lzp2、lzp7、“張想”的親緣系數(shù)(r)分別為0．05、0．10、0．10，但是r<0．125，“淘淘”和這三只個(gè)體的親緣關(guān)系水平低于雙重表親(Double first cousins)，不存在明顯的因親緣關(guān)系較近而帶來(lái)的近交風(fēng)險(xiǎn)。“張想”和lzp6、“淘淘”、“瀘欣”的親緣系數(shù)(r)分別為0．03、0．10、0．11，r<0．125，“張想”和這三只個(gè)體的親緣關(guān)系水平也低于雙重表親(Double first cousins)，不存在明顯的因親緣關(guān)系近而帶來(lái)的近交風(fēng)險(xiǎn)。因lzp7是雌性大熊貓，雖然它和“張想”親緣系數(shù)(r)為0．30，屬于親子/全同胞(parent-offspring/full sibs)的水平，但也不存在因親緣關(guān)系近的近交風(fēng)險(xiǎn)。

表4　研究區(qū)域大熊貓個(gè)體親緣關(guān)系

3．3．2加入放歸個(gè)體后的遺傳結(jié)構(gòu)

根據(jù)哈溫平衡檢驗(yàn)分析結(jié)果，加入放歸個(gè)體后的大熊貓小種群仍未偏離哈溫平衡(p>0．05)。連鎖不平衡分析結(jié)果顯示對(duì)于加入放歸個(gè)體后的大熊貓小種群，也不存在連鎖遺傳。12個(gè)微衛(wèi)星位點(diǎn)分析結(jié)果顯示：加入放歸個(gè)體后種群平均觀測(cè)雜合度(Ho)為0．557，平均期望雜合度(He)為0．581(見(jiàn)表3)。相比Zhang等[24]的研究，加入放歸個(gè)體后的遺傳多樣性仍處于中等偏上的水平。平均等位基因數(shù)(MNA)增加到3．667，表現(xiàn)在Ame-μ11(Locus1)、Ame-μ22(Locus3)、Ame-μ24(Locus7)、Ame-μ13(Locus10)和Ame-μ27(Locus12)5個(gè)位點(diǎn)上，共計(jì)增加6 個(gè)等位基因。 STRUCTURE分析結(jié)果表明，K=4時(shí)(見(jiàn)圖3C)，LnP(D)具有最大的平均值，遺傳基因簇由一簇變?yōu)榱怂拇?見(jiàn)圖3D)，這表明放歸個(gè)體加入后，出現(xiàn)了新的基因頻率。

4　討 論

三只放歸個(gè)體加入當(dāng)?shù)胤N群后，研究區(qū)域大熊貓小種群的等位基因數(shù)有所增加，同時(shí)也出現(xiàn)了新的基因頻率，這表明該區(qū)大熊貓小種群的遺傳多樣性確實(shí)呈增長(zhǎng)趨勢(shì)。遺傳多樣性的高低在一定程度上反映了大熊貓種群的進(jìn)化潛力[25]，因此放歸個(gè)體的加入有助于研究區(qū)域大熊貓小種群的種群遺傳多樣性的提高，將對(duì)該種群的復(fù)壯帶來(lái)積極影響。

本研究分析結(jié)果顯示放歸個(gè)體與研究區(qū)域本地大熊貓小種群親緣關(guān)系較遠(yuǎn)。楊波等[26]的研究表明：親緣關(guān)系越遠(yuǎn)，釋放個(gè)體與本地種群組成的混合種群內(nèi)近交率會(huì)越低，從而減少近交衰退的風(fēng)險(xiǎn)，更有利于改善野生種群遺傳水平。因此這三只放歸個(gè)體的加入，能有效提高研究區(qū)域本地大熊貓小種群的遺傳多樣性水平。

迄今，“瀘欣”、“淘淘”、“張想”已放歸分別超過(guò)6年、3年、2年，它們?cè)诖笮茇埛艢w地——栗子坪國(guó)家級(jí)自然保護(hù)區(qū)生存狀況良好。年齡組結(jié)構(gòu)分析顯示，加入放歸個(gè)體后，研究區(qū)域大熊貓小種群年齡結(jié)構(gòu)得到進(jìn)一步改善。此外，該區(qū)域的紅外相機(jī)監(jiān)測(cè)結(jié)果顯示，“瀘欣”已成功帶崽，表明放歸個(gè)體已經(jīng)初步融入當(dāng)?shù)胤N群，放歸工作取得初步成效。若其余兩只放歸個(gè)體也能成功繁殖，將進(jìn)一步促進(jìn)該地區(qū)大熊貓小種群的復(fù)壯。

目前，小相嶺山系大熊貓種群數(shù)量仍然很少，成功放歸三只大熊貓遠(yuǎn)不足以復(fù)壯當(dāng)?shù)匦》N群。針對(duì)這一問(wèn)題，該山系還需進(jìn)一步持續(xù)開(kāi)展大熊貓放歸工作。而在放歸個(gè)體的選擇上，需要比較放歸個(gè)體和當(dāng)?shù)匦》N群的遺傳背景，優(yōu)先選擇能有效提高當(dāng)?shù)匦》N群遺傳多樣性的個(gè)體。

[1]ZHU L F, ZHAN X J, MENG T, et al． Landscape features influence gene flow as measured by cost-distance and genetic analyses: a case study for giant pandas in the Daxiangling and Xiaoxiangling mountains [J]． BMC genetics, 2010(1):72．

[2]ZHU L F, ZHANG S N, GU X D, et al． Significant genetic boundaries and spatial dynamics of giant pandas occupying fragmented habitat across southwest China [J]． Molecular Ecology, 2011, 20(6):1122-1132．

[3]BIEBACH I, KELLER L F． Inbreeding in reintroduced populations: the effects of early reintroduction history and contemporary processes [J]． Conservation Genetics, 2010(2):527-538．

[4]ALCAIDE M．, NEGRO J J．, SERRANO D, et al． Captive breeding and reintroduction of the lesser kestrel Falco naumanni: a genetic analysis using microsatellites [J]． Conservation Genetics, 2010(11): 331-338．

[5]GREY-ROSS R, DOWNS C T, KIRKMAN K．． Reintroduction Failure of Captive-Bred Oribi (Ourebia ourebi) [J]． South African Journal of Wildlife Research, 2009, 39(1): 34-38．

[6]朱磊, 吳攀文, 張洪, 等． 小相嶺山系大熊貓種群生存力分析 [J]． 西華師范大學(xué)學(xué)報(bào)( 自然科學(xué)版), 2008, 29(2): 112-116．

[7]ZHU L F, ZHAN X J, WU H, et al． Conservation implications of drastic reductions in the smallest and most isolated populations of giant pandas [J]． Conservation Biology, 2010, 24(5): 1299-1306．

[8]青菁, 胥池, 楊彪, 等． 小相嶺山系大熊貓廊道規(guī)劃[J]． 生態(tài)學(xué)報(bào), 2016, 36(4): 1-9．

[9]四川省林業(yè)廳． 四川的大熊貓[M]． 成都：四川科學(xué)技術(shù)出版社，2015.

[10]詹祥江．利用非損傷性遺傳取樣研究大熊貓的種群數(shù)量和擴(kuò)散模式 [D]． 中國(guó)科學(xué)院動(dòng)物研究院, 2006．

[11]胡錦矗． 從野外大熊貓的糞便估計(jì)年齡及其種群年齡結(jié)構(gòu)的研究 [J]． 獸類(lèi)學(xué)報(bào), 1987, 7 (2): 81-84．

[12]胡義波． 涼山山系大熊貓的種群歷史、景觀與空間遺傳格局研究 [D]． 北京：中國(guó)科學(xué)院研究生院, 2008．

[13]VAN OOSTERHOUT C, HUTCHINSON W F, WILLS D P M, et al． Micro-checker: software for identifying and correcting genotyping errors in microsatellite data [J]． Molecular Ecology Notes, 2004(3): 535-538．

[14]VALIERE N． GIMLET: a computer program for analysing genetic individual identification data [J]． Molecular Ecology Notes, 2002, 2: 377-379．

[15]MILLS L S, CITTA J J, LAIR K P, et al． Estimating animal abundance using noninvasive DNA sampling: Promise and pitfalls [J]． Ecological Applications, 2000, 10(1):283-294．

[16]PARK S D E． Trypanotolerance in west African cattle and the population genetic effects of selection [D]． Dublin：University of Dublin, 2001．

[17]BELLEMAIN E, SWENSON J E, TALLMON D, et al． Estimating population size of elusive animals with DNA from hunter-collected feces: Four methods for brown bears [J]． Conservation Biology, 2005, 19(1): 150-161．

[18]EXCOFFIER L, LAVAL G, SCHNEIDER S． Arlequin (version 3．0): an integrated software package for population genetics data analysis [J]． Evolutionary Bioinformatics Online, 2005, 1: 47-50．

[19]RAYMOND M, ROUSSET F． Genepop (Version1．2): Population Genetics Software for Exact Tests and Ecumenicism [J]． Journal of Heredity, 1995,86(3):248-249．

[20]GOUDET J． FSTAT, a program to estimate and test gene diversities and fixationindices (version 2．9．3)． 2001． Updated from Goudet J． FSTAT (Version 1．2): A computer program to calculate F-statistics [J]． Journal of Heredity, 1995, 86: 485-486．

[21]PRITCHARD J K, STEPHENS M, DONNELLY P． Inference of population structure using multilocus genotype data [J]． Genetics, 2000, 155(2): 945-959．

[22]KONOVALOV D A, MANNING C, HENSHAW M T． Kingrouf: a program for pedigree relationship reconstruction and kin group assignments using genetic markers [J]． Molecular Ecology Notes, 2004(4): 779-782．

[23]WANG J． COANCESTRY: a program for simulating, estimating and analysing relatedness and inbreeding coefficients [J]． Molecular Ecology Resources, 2011, 11(1): 141-145．

[24]ZHANG B W, LI M, ZHANG Z J, et al． Genetic Viability and Population History of the Giant Panda, Putting an end to the "Evolutionary dead end"? [J]． Molecular Biology and Evolution, 2007, 24(18): 1801-1810．

[25]張澤鈞, 張陜寧, 魏輔文, 等． 移地與圈養(yǎng)大熊貓野外放歸的探討 [J]． 獸類(lèi)學(xué)報(bào), 2006, 26 (3): 292-299．

[26]楊波, 楊承忠, 涂飛云, 等． 大熊貓野化放歸中的遺傳學(xué)分析 [J]． 四川動(dòng)物, 2013, 32(1):149-155．

[Abstract]The population of giant panda in Xiaoxiangling is one of the smallest . In order to rejuvenate this small panda population, the State Forestry Administration has tried to release three giant pandas(Ailuropoda melanoleuca)： “Luxin”, “Taotao” and “Zhangxiang” in the Liziping National Nature Reserve. To assess the genetic structure of this small population after reintroduction, 92 stools of giant pandas were collected in research area from 2012 to 2013, and 14 individuals(including three reintroduced individuals)were identified. After reintroduced the three individuals, measurement proved that the kinship between the reintroduced individuals and local small population was far. The population age composition was: eight individuals in youth group，four individuals in adult group and two individuals in unknown group.The mean observed heterozygosity (Ho) was 0.557. The mean expected heterozygosity (He) was 0.581. Genetic diversity was in the medium level.Five sites increased six alleles.And the mean allele number (MNA) achieved 3.667. What’s more ,the genetic cluster increased to four meaning that new gene frequency had appeared. All these results suggest that the genetic diversity of the small giantpanda population in this region is improved. The reintroduction is beginning to bear results.

[Key words]giant panda(Ailuropoda melanoleuca)； small population； genetic diversity； reintroduction； Liziping National Nature Reserve

[責(zé)任編輯劉景平]

The Genetic Structure of Small Giant Panda Population in the Liziping Nature Reserve after Reintroduction

CHEN Zi-jun1, HUANG Feng1,2, HE Ke1, FENG Fei-fei1, GU Xiao-dong3, YANG Xu-yu3,DAI Qiang4, QI Dun-wu5, YANG Zhi-song1

(1.Key Laboratory of Southwest China Wildlife Resource Conservation, China West Normal University,Nanchong, Sichuan 637002, China;2.Sichuan Liziping Nature Reserve Management Office, Shimian, Sichuan 625400, China;3.Sichuan Station of Wildlife Survey and Management, Chengdu, Sichuan 610081, China;4.Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu, Sichuan 610041, China;5.Chengdu Research Base of Giant Panda Breeding, Chengdu, Sichuan 610086, China)

Q953

A

1672-9021(2016)02-0001-07

陳子君(1989-)，女，四川綿陽(yáng)人，西華師范大學(xué)生命科學(xué)學(xué)院碩士研究生，主要研究方向：大熊貓保護(hù)生物學(xué)；楊志松(1975-)，男，四川南充人，博士，西華師范大學(xué)生命科學(xué)學(xué)院副教授，碩士生導(dǎo)師，主要研究方向：大熊貓保護(hù)生物學(xué)、分子系統(tǒng)地理學(xué)、 鳥(niǎo)類(lèi)生態(tài)學(xué)。

香港海洋公園保育基金資助項(xiàng)目(GP07_1213)。

2016-03-01